Abstract: Groundwater fluxes in tidal marshes largely control key ecosystem functions and services, such as vegetation growth, soil carbon sequestration, and nutrient cycling. In tidal marshes restored on formerly embanked agricultural land, groundwater fluxes are often limited as compared to nearby natural marshes, as a result of historical agricultural soil compaction. To improve the functioning of restored tidal marshes, knowledge is needed on how much certain design options can optimize soil-groundwater interactions in future restoration projects. Based on measured data on soil properties and tidally induced groundwater dynamics, we calibrated and evaluated a 2D vertical model of a creek-marsh cross-section, accounting for both saturated and unsaturated groundwater flow and solute transport in a variably saturated groundwater flow model. We found that model simulations of common restoration practices such as soil amendments (increasing the depth of porous soil on top of the compact layer) and creek excavation (increasing the creek density) increase the soil aeration depth and time, the drainage depth and the solute flux, and decrease the residence time of solutes in the porewater. Our simulations indicate that increasing the depth to the compact layer from 20 cm to 40 cm, or increasing the creek density from 1 creek to 2 creeks along a 50 m marsh transect (while maintaining the total creek cross-sectional area), in both cases more than doubles the volume of water processed by the marsh soil. We discuss that this may stimulate nutrient cycling. As such, our study demonstrates that groundwater modelling can support the design of marsh restoration measures aiming to optimize groundwater fluxes and related ecosystem services.

Abstract: Human activities related to urbanization and industrialization have established a vast territory of urban soil worldwide. On traditional soil maps, urban and industrial areas usually appear as blind spots as they were beyond the interest of national soil survey campaigns. Furthermore, these soil maps are likely already outdated with respect to urban soil due to rapid urban expansion in recent decades. This research aims to evaluate the use of land cover data to estimate the urban soil information gap considering the highly urbanized region of Flanders, Belgium, as a case study. The current extent and spatial distribution of anthropogenic urban soil (1) was estimated through reclassification of recently acquired (2012) exhaustive land cover data, discriminating three qualitative likelihood levels (high-intermediate-low) of anthropogenic influence by urbanization, and (2) compared with its occurrence as represented by the 'Technosols/Not Surveyed area' in the legacy soil map of Belgium, as this map unit best matches with the likelihood for anthropogenic urban soil at the time of the National Soil Survey conducted between end 1940s and mid 1970s. The proposed reclassification of the land cover map resulted in 16.3% and 16.7% of Flanders' total area that corresponds with a high and intermediate likelihood for anthropogenic urban soil, which highlights the underestimation of the anthropogenic urban soil extent as represented by the 'Technosol/Not Surveyed' unit in the legacy soil map (only 13.7%). Moreover, a more realistic spatial pattern of anthropogenic urban soil occurrence was obtained, providing an improved basis for urban soil spatial analysis studies. The produced anthropogenic urban soil likelihood map therefore presents a useful supporting tool for coordinating future soil surveys in urban environments.

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: To properly estimate and manage pesticide occurrence in urban rivers, it is essential, but often highly challenging, to identify the key pesticide transport pathways in association to the main sources. This study examined the concentration-discharge hysteresis behaviour (hysteresis analysis) for three pesticides and the parent-metabolite concentration dynamics for two metabolites at sites with different levels of urban influence in a mixed land use catchment (25 km(2)) within the Swiss Greifensee area, aiming to identify the dominant pesticide transport pathways. Combining an adapted hysteresis classification framework with prior knowledge of the field conditions and pesticide usage, we demonstrated the possibility of using hysteresis analysis to qualitatively infer the dominant pesticide transport pathway in mixed land-use catchments. The analysis showed that hysteresis types, and therefore the dominant transport pathway, vary among pesticides, sites and rainfall events. Hysteresis loops mostly correspond to dominant transport by flow components with intermediate response time, although pesticide sources indicate that fast transport pathways are responsible in most cases (e.g. urban runoff and combined sewer overflows). The discrepancy suggests the fast transport pathways can be slowed down due to catchment storages, such as topographic depressions in agricultural areas, a wastewater treatment plant (WWTP) and other artificial storage units (e.g. retention basins) in urban areas. Moreover, the WWTP was identified as an important factor modifying the parent-metabolite concentration dynamics during rainfall events. To properly predict and manage pesticide occurrence in catchments of mixed land uses, the hydrological delaying effect and chemical processes within the artificial structures need to be accounted for, in addition to the catchment hydrology and the diversity of pesticide sources. This study demonstrates that in catchments with diverse pesticide sources and complex transport mechanisms, the adapted hysteresis analysis can help to improve our understanding on pesticide transport behaviours and provide a basis for effective management strategies.(C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Large river catchments with mixed land use capture pesticides from many sources, and degradable pesticides are converted during downstream transport. Unravelling the contribution of pesticide source and the effect of degradation processes is a challenge in such areas. However, insight and understanding of the sources is important for targeted managen-lent, especially when water is abstracted from the river for drinking water production. The river Meuse is such a case. A long-term monitoring data set was applied in a modelling approach for assessing the contribution of waste water treatment plants (WWTPs) and tributaries (sub-basins) to surface water contamination, and to evaluate the effect of decay on the downstream concentrations of glyphosate and AMPA at the point of drinking water abstraction. The results show that WWTPs are important contributors for glyphosate and AMPA in large river catchments with mixed land uses. In the studied area, the river Meuse in the Netherlands, the relative contribution of WWTP effluents is above 29% for glyphosate and around 12% for AMPA. Local industries are found to be potentially big contributors of AMPA. Glyphosate entering the river system is gradually converted to AMPA and other degradation-productions, which results in downstream loads that are considerably lower than the sum of all influxes. In summer when the travel time is longer due to lower discharge, the first order decay of glyphosate in the river Meuse is estimated to result in about 50% reduction of the downstream glyphosate concentrations over a river stretch of 250 km. The contribution of glyphosate decay to the observed AMPA concentrations ranges between 2% and 10%. Contributions are sensitive to seasonal variations in discharge that influence the concentrations through dilution and degradation. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Monitoring and modelling tools may improve irrigation strategies in precision agriculture. We used non-invasive soil moisture monitoring, a crop growth and a soil hydrological model to predict soil water content fluctuations and crop yield in a heterogeneous sandy grassland soil under supplementary irrigation. The sensitivity of the soil hydrological model to hydraulic parameters, water stress, crop yield and lower boundary conditions was assessed after integrating models. Free drainage and incremental constant head conditions were implemented in a lower boundary sensitivity analysis. A time-dependent sensitivity analysis of the hydraulic parameters showed that changes in soil water content are mainly affected by the soil saturated hydraulic conductivity K-s and the Mualem-van Genuchten retention curve shape parameters n and alpha. Results further showed that different parameter optimization strategies (two-, three-, four- or six-parameter optimizations) did not affect the calculated water stress and water content as significantly as does the bottom boundary. In this case, a two-parameter scenario, where K-s was optimized for each layer under the condition of a constant groundwater depth at 135-140 cm, performed best. A larger yield reduction, and a larger number and longer duration of stress conditions occurred in the free drainage condition as compared to constant boundary conditions. Numerical results showed that optimal irrigation scheduling using the aforementioned water stress calculations can save up to 12-22 % irrigation water as compared to the current irrigation regime. This resulted in a yield increase of 4.5-6.5 %, simulated by the crop growth model.

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.

Abstract: Undoubtedly, the most important advance in the environmental regulatory monitoring of elements of the last decade is the widespread introduction of ICP-mass spectrometry (ICP-MS) due to standards developed by the European Committee for Standardization. The versatility of ICP-MS units as a tool for the determination of major, minor and trace elements (Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, K, Mg, Mn, Mo, Na, Ni, P, Pb, Sb, Se, Sn, Ti, V and Zn) in surface water, groundwater, river sediment, topsoil, subsoil, fine particulates and atmospheric deposition is illustrated in this paper. Ranges of background concentrations for major, minor and trace elements obtained from a regional case study (Flanders, Belgium) are summarized for all of these environmental compartments and discussed in the context of a harmonized implementation of European regulatory monitoring requirements. The results were derived from monitoring programs in support of EU environmental quality directives and were based on a selection of (non-polluted) background locations. Because of the availability of ICP-MS instruments nowadays, it can be argued that the main hindrance for meeting the European environmental monitoring requirements is no longer the technical feasibility of analysis at these concentration levels, but rather (i) potential contamination during sampling and analysis, (ii) too limited implementation of quality control programs, validating the routinely applied methods (including sampling and low level verification) and (iii) lack of harmonization in reporting of the chemical environmental status between the individual member states.

Abstract: The latest life-cycle assessment methods account for land use, due to the production, use and disposal of products and services, in terms of ecosystem damage. The process of brownfield remediation converts otherwise idle urban space into productive space. The value to ecosystems in this context is of course limited since the brownfield site remains urban. When evaluating brownfield remediation technologies, the availability of space on-site is dependent on the duration of time required by the remediation technology to reach the remediation target. Remediation technology alternatives tend to vary largely in terms of duration. Comparative life-cycle assessments of remediation technologies, to date, present the large variations between alternatives in terms of remediation duration but do not translate this into an impact or parameter. The restored subsurface zone is often defined as a functional unit, when in fact the surface area is the resource restored by the remediation service. The economic benefits of making land resources available are particularly important considerations in the context of brownfield remediation. The research proposes an innovative impact assessment approach that allows land to be considered as a finite resource. The method is applied in a comparative life-cycle assessment of two potential remediation scenarios for an idle brownfield in the Brussels region of Belgium. The results show that there is a trade-off between greenhouse gas emissions and land availability and that both are largely dependent on the efficiency of the contaminant extraction mechanism. The results also raise the question as to whether the economic valuation of land, like precious metals and fossil fuels, provides an accurate reflection of the true value of the resource. Considering land as a resource at the midpoint level is also relevant in other urban contexts where competition exists between different land-uses, where urban sprawl is detrimental to undeveloped areas and where urban intensification is a policy objective. (C) 2015 Elsevier Ltd. All rights reserved.

Abstract: Urban runoff can be a significant source of pesticides in urban streams. However, quantification of this source has been difficult because pesticide use by urban residents (e.g., on pavements or in gardens) is often unknown, particularly at the scale of a residential catchment. Proper quantification and characterization of pesticide loss via urban runoff require sound information on the use and occurrence of pesticides at hydrologically-relevant spatial scales, involving various hydrological conditions. We conducted a monitoring study in a residential area (9.5 ha, Flanders, Belgium) to investigate the use and loss of a widely-used herbicide (glyphosate) and its major degradation product (aminomethylphosphonic acid, AMPA). The study covered 13 rainfall events over 67 days. Overall, less than 0.5% of glyphosate applied was recovered from the storm drain outflow in the catchment. Maximum detected concentrations were 6.1 mu g/L and 5.8 mu g/L for glyphosate and AMPA, respectively, both of which are below the predicted no-effect concentration for surface water proposed by the Flemish environmental agency (10 mu g/L), but are above the EU drinking water standard (0.1 mu g/L). The measured concentrations and percentage loss rates can be attributed partially to the strong sorption capacity of glyphosate and low runoff potential in the study area. However, glyphosate loss varied considerably among rainfall events and event load of glyphosate mass was mainly controlled by rainfall amount, according to further statistical analyses. To obtain urban pesticide management insights, robust tools are required to investigate the loss and occurrence of pesticides influenced by various factors, particularly the hydrological and spatial factors. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: The use of temperature-time series measured in streambed sediments as input to coupled water flow and heat transport models has become standard when quantifying vertical groundwater-surface water exchange fluxes. We develop a novel methodology, called LPML, to estimate the parameters for 1-D water flow and heat transport by combining a local polynomial (LP) signal processing technique with a maximum likelihood (ML) estimator. The LP method is used to estimate the frequency response functions (FRFs) and their uncertainties between the streambed top and several locations within the streambed from measured temperature-time series data. Additionally, we obtain the analytical expression of the FRFs assuming a pure sinusoidal input. The estimated and analytical FRFs are used in an ML estimator to deduce vertical groundwater-surface water exchange flux and its uncertainty as well as information regarding model quality. The LPML method is tested and verified with the heat transport models STRIVE and VFLUX. We demonstrate that the LPML method can correctly reproduce a priori known fluxes and thermal conductivities and also show that the LPML method can estimate averaged and time-variable fluxes from periodic and nonperiodic temperature records. The LPML method allows for a fast computation of exchange fluxes as well as model and parameter uncertainties from many temperature sensors. Moreover, it can utilize a broad frequency spectrum beyond the diel signal commonly used for flux calculations.

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: A spatially distributed model for leaching of Cd from the unsaturated zone was developed for the Belgian-Dutch transnational Kempen region. The model uses as input land-use maps, atmospheric deposition data, and soil data and is part of a larger regional model that simulates transport of Cd in soil, groundwater, and surface water. A new method for deriving deposition from multiple sites was validated using soil data in different wind directions. Leaching was calculated for the period 1890 to 2010 using a reconstruction of metal loads in the region. The model was able to reproduce spatial patterns of concentrations in soil and groundwater and predicted the concentration in shallow groundwater adequately well for the purpose of evaluating management options. For 42% of the data points, measurements and calculations were within the same concentration class. The model was used for forecasting under a reference scenario, an autonomous development scenario including climate change, and a scenario with implementation of remediation measures. The impact of autonomous development (under the most extreme scenario of climatic change) amounted to an increase of 10% in cumulative Cd flux after 100 yr as compared with the reference scenario. The impact of remediation measures was mainly local and is less pronounced (i.e., only 3% change in cumulative flux at the regional scale). The integrated model served as a tool to assist in developing management strategies and prioritization of remediation of the wide-spread heavy metal contamination in the region.

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: 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.