Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.     

Graphical tools based on Turc-Budyko plots to detect changes in catchment behaviour

Abstract

This paper aims at revisiting the use of non-dimensional representations for catchment change and model robustness analysis. As well as being helpful for catchment classification according to hydroclimatic conditions, the Turc-Budyko representation enables visualization of temporal changes in these conditions and/or modification of the hydrological behaviour of the catchment. This brings a new perspective to hydrological analysis, different from the classical time series plots, which helps when interpreting catchment functioning with respect to hydroclimatic constraints. These tools do not require statistical analyses of observed series and are therefore very simple to implement. Four case studies are considered here to illustrate the use of the Q/P = f(P/E_p) graph. When combined with the inter-annual Turc-Mezentsev formula, this visual framework enables anticipation of potential difficulties for models based on the same hypotheses, solely using the analysis of observed data.     

Determinants of sediment properties and organic matter in beach and dune environments based on boosted regression trees

Grain size properties and the variation of organic matter in coastal beach and dune environments are assumed to be controlled by the intensity of aeolian processes, time and the sediment source. However, assumptions are based on relatively limited empirical studies. In this study, we examined which environmental variables are the main predictors of multiple topsoil properties. To achieve this, we analysed an extensive dataset systematically collected across all beach zones and a large geographical area at the Finnish Baltic Sea coast characterized by post‐glacial land uplift. We included a comprehensive set of predictors in the analysis and applied boosted regression trees, a modern modelling technique particularly suited for analysis without prior assumptions of the data model. The results suggest that mean grain size and sorting are mainly determined by northing and fetch. Northing, disturbance and fetch predicted the variation of soil organic matter while litter cover was strongly related to disturbance. Based on the analyses, we were able to identify the main drivers of multiple topsoil properties on land uplift beaches. Parent material is suggested to determine sediment textural properties, which largely masks the effects of transient processes. Mean grain size and sorting are highly interdependent: grains become finer and sorting improves with increasing shore exposure. The intensity of momentary geomorphic processes controls the accumulation of litter whereas the slower accumulation of organic matter in the soil is influenced also by the static exposure setting. Skewness and kurtosis of the grain size distribution are mainly influenced by unmeasured processes, potentially relating to the geomorphological origin of the sediment. Copyright © 2014 John Wiley & Sons, Ltd.     

Hydro‐meteorological drivers and sources of suspended sediment flux in the pro‐glacial zone of the retreating Castle Creek Glacier,Cariboo Mountains,British Columbia,Canada

Glaciers are major agents of erosion that increase sediment load to the downstream fluvial system. The Castle Creek Glacier, British Columbia, Canada, has retreated ~1.0 km in the past 70 years. Suspended sediment concentration (SSC) and streamflow (Q) were monitored independently at five sites within its pro‐glacial zone over a 60 day period from July to September 2011, representing part of the ablation season. Meteorological data were collected from two automatic weather stations proximal to the glacier. The time‐series were divided into hydrologic days and the shape and magnitude of the SSC response to hydro‐meteorological conditions (‘cold and wet’, ‘hot and dry’, ‘warm and damp’, and ‘storm’) were categorized using principal component analysis (PCA) and cluster analysis (CA). Suspended sediment load (SSL) was computed and summarized for the categories. The distribution of monitoring sites and results of the multivariate statistical analyses describe the temporal and spatial variability of suspended sediment flux and the relative importance of glacial and para‐glacial sediment sources in the pro‐glacial zone. During the 2011 study period, ~ 60% of the total SSL was derived from the glacial stream and sediment deposits proximal to the terminus of the glacier; during ‘storm’ events, that contribution dropped to ~40% as the contribution from diffuse and point sources of sediment throughout the pro‐glacial zone and within the meltwater channels increased. While ‘storm’ events accounted for just 3% of the study period, SSL was ~600% higher than the average over the monitoring period, and ~20% of the total SSL was generated in that time. Determining how hydro‐meteorological conditions and sediment sources control sediment fluxes will assist attempts to predict how pro‐glacial zones respond to future climate changes. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of heavy metals and arsenic speciation discharged by the industrial activity on the Tinto-Odiel estuary, SW Spain

This study reports the annual amount of heavy metals discharged by industrial activity into the estuary of the Ría of Huelva (SW Spain). The findings showed that the discharged metals found in highest amounts were Fe (11 t y⁻¹), Zn (3.4 t y⁻¹) and Mo (0.88 t y⁻¹). There were other metals with high pollutant charge, such as Ti (232 kg y⁻¹), As (228 kg y⁻¹), Ni (195 kg y⁻¹), Pb (100 kg y⁻¹), Cr (39 kg y⁻¹) and Cd (33 kg y⁻¹). These results were compared with pollutants transported via the Tinto and Odiel rivers from abandoned mining activities in the Iberian Pyrite Belt (IPB), and it was deduced that the amounts spilled exclusively by industries were less than 1% in relation to the total discharge. Hence, the treatment of residues from the IPB should be the priority goal to improve water quality in the estuary.     

Spatial and temporal distribution of heavy metal concentrations in mussels (Mytilus edulis) from the Baie des Chaleurs, New Brunswick, Canada

Previous studies on heavy metal contamination of the Baie des Chaleurs focus only on industrial centers and overlooked the ecosystem as a whole. To fill this gap, the objective of this study is to establish a baseline of the spatio-temporal distribution of heavy metals in mussels from the Baie des Chaleurs based on the ecosystem approach. Our results show, for the first time, a cadmium contamination in mussels across the south coast of the Baie des Chaleurs and not only in industrial centers. Our results also confirm previous studies showing heavy metal contamination of the Belledune area. This study demonstrates that the use of the ecosystemic approach is essential to obtain a comprehensive picture of environmental contamination in marine ecosystems.     

Do benthic biofilters contribute to sustainability and restoration of the benthic environment impacted by offshore cage finfish aquaculture?

Benthic biofilters were deployed under a cage fish farm and in two reference locations to assess the influence of the farm on the biofilters and the surroundings, as well as to verify the usefulness of this technology as a mitigation tool. The biofilters underneath the farm recruited a fouling community practically identical to that of the control biofilters, which included a variety of trophic strategies. The former showed a higher ¹⁵N enrichment, indicating that fouling beneath the farm was benefiting from the farm waste. The waste retention efficiency was low (0.02 g N m⁻² month⁻¹) beneath the farm. Benthic biofilters aggregated demersal wild fish around and within them. Pelagic wild fish also frequently used the biofilters beneath the farm, forming compact shoals around them. The increased complexity of the habitat below the fish farm enhanced biodiversity, but this improvement did not lead to the recovery of the sediments around the biofilters.