Flow and nutrient transport in intermittent rivers: a modelling case-study on the Vène River using SWAT 2005

Abstract

Intermittent rivers have a specific hydrological behaviour which also influences water quality dynamics. The objective of this work was to model the flow and water quality dynamics of a coastal Mediterranean intermittent river using the Soil and Water Assessment Tool (SWAT 2005). Flow, sediment, nitrogen and phosphorus transport were simulated on the Vène experimental catchment, France. The model was sequentially calibrated at sub-catchment scale and validated both at sub-catchment and catchment scales. A procedure for building the data records for the point sources is presented. The results indicate that, while the model produces good results for flow simulation, its performance for sediment transport is less satisfactory. This in turn impacts on the nutrient transport module. The reasons behind these shortcomings are analysed, taking into account the length of the data records, their distribution and the equations used in the SWAT model. The need for a thorough multi-objective model validation is illustrated.

Citation Chahinian, N., Tournoud, M.-G., Perrin, J.-L. & Picot, B. (2011) Flow and nutrient transport in intermittent rivers: a modelling case-study on the Vène River using SWAT 2005. Hydrol. Sci. J. 56(2), 268–287.     

Numerical modelling of morphodynamics—Vilaine Estuary

The main objective of this paper is to develop a method to simulate long-term morphodynamics of estuaries dominated by fine sediments, which are subject to both tidal flow and meteorologically induced variations in freshwater run-off and wave conditions. The method is tested on the Vilaine Estuary located in South Brittany, France. The estuary is subject to a meso–macrotidal regime. The semi-diurnal tidal range varies from around 2.5 to 5 m at neap and spring, respectively. The freshwater input is controlled by a dam located approximately 8 km from the mouth of the estuary. Sediments are characterised as mostly fines, but more sandy areas are also found. The morphology of the estuary is highly influenced by the dam. It is very dynamic and changes in a complicated manner with the run-off from the dam, the tide and the wave forcing at the mouth of the estuary. Extensive hydrodynamic and sediment field data have been collected in the past and provide a solid scientific basis for studying the estuary. Based on a conceptual understanding of the morphodynamics, a numerical morphological model with coupled hydrodynamic, surface wave and sediment transport models is formulated. The numerical models are calibrated to reproduce sediment concentrations, tidal flat altimetry and overall sediment fluxes. Scaling factors are applied to a reference year to form quasi-realistic hydrodynamic forcing and river run-off, which allow for the simulations to be extended to other years. The simulation results are compared with observed bathymetric changes in the estuary during the period 1998–2005. The models and scaling factors are applied to predict the morphological development over a time scale of up to 10 years. The influence of the initial conditions and the sequence of external hydrodynamic forcing, with respect to the morphodynamic response of the estuary, are discussed.     

The Multiscale TROPIcal CatchmentS critical zone observatory M-TROPICS dataset I: The Nyong River Basin,Cameroon

Simultaneously acquiring time series of climate, hydrology and hydrochemical data over decades on river systems is pivotal to understand the complex interactions involving rock, soil water, air and biota in the Critical Zone, to build integrated modelling and to propose predictive scenarios. Among the Critical Zone Observatories (CZOs) implemented in the past 25 years, only a few are located in the humid Tropics despite the importance of these regions in terms of population density, fast-changing land use, biodiversity hotspots, biomass stock on continents, size of river systems, etc. Since 1994, weathering and erosion processes and fluxes have been investigated at both local (experimental watershed) and regional scales in the Nyong River Basin (Cameroon) which belongs to the Critical Zone Observatories network named Multiscale TROPIcal CatchmentS (M-TROPICS). The data shared by M-TROPICS in Cameroon are: (1) rainfall; (2) air temperature, air relative humidity, wind speed and direction, and global radiation; (3) stream and river water level; (4) pH, electrical conductivity, water temperature and suspended particulate matter (SPM) concentration; (5) major ion, alkalinity and dissolved organic carbon (DOC) concentrations. The dataset already contributed to describe the water partitioning in these tropical humid watersheds, to better understand the factors controlling chemical weathering and physical erosion in tropical ecosystems, particularly the role of organic matter. The dataset also contributed to calculate elemental weathering fluxes and saprolite production rate and to propose denudation rates on tropical cratonic landscapes. Hydrological modelling allowed quantification of the geographical water sources contributing to streamflow. DOC data were used to determine greenhouse-gas emissions and carbon budgets from African inland waters. However, long-term solute concentrations at the outlet of a small tributary of the Nyong River exhibit non-stationary behaviour over the last 26 years. The processes governing those fluctuations are not yet fully understood and might be related to changes in the hydrological regime, land-cover and land-use. The latter highlights the need for longer time-series and continued support for CZOs particularly in the humid tropics.     

Measurement of total dissolved phosphorus in small volumes of iron rich interstitial water

It is shown that sorption of orthophosphate to iron compounds, formed during persulfate digestion, can cause a significant underestimation of total dissolved phosphorus in interstitial waters rich in iron and poor in phosphorus. Labelling the samples with carrier free³²PO₄ before digestion allows to correct for these losses.     

Cadmium deposition and mobility in the sediments of an acidic oligotrophic lake

A dated core from the profoundal zone in a pristine oligotrophic acidic lake was analyzed for Cd as well as for Al, Ca, Fe, Mg, Mn, Pb, Ti and total carbon and nitrogen. Overlying water and porewater samples were also obtained on six occasions at the same site, and yielded vertical profiles of pH and dissolved Cd, Ca, Fe, Mg, Mn, sulfide, SO₄⁻², organic and inorganic carbon concentrations. These extensive porewater and sediment geochemical data were used, together with information on infaunal benthos, to decipher the sedimentary record of Cd contamination. Depth variation of sediment Ca concentrations indicate that the lake suffered from progressive acidification starting about 1950. The present-day accumulation rate of Cd (J_acc^Cd = 5.4 ± 0.4 × 10⁻¹¹ mol cm⁻² yr⁻¹) in the sediments is the sum of the flux of Cd deposited with settling particles (J_S^Cd = 3.3 ± 0.2 × 10⁻¹¹ mol cm⁻² yr⁻¹) and the fluxes of dissolved Cd across the sediment-water interface due to molecular diffusion (J_D^Cd = 1.8 ± 0.3 × 10⁻¹¹ mol cm⁻² yr⁻¹), bioturbation (J_B^Cd = 1.1 ± 0.2 × 10⁻¹⁴ mol cm⁻² yr⁻¹) and bioirrigation (J_I^Cd = 0.27 ± 0.05 × 10⁻¹¹ mol cm⁻² yr⁻¹). Biological mixing of the sediments was negligible. The shape of the vertical profile of total Cd concentration with depth in the sediment appears to be determined more by its input history than by post-depositional mobilization and redistribution in the sediment column.     

An experimental approach to investigate mercury species transformations under redox oscillations in coastal sediments

This work describes a laboratory experiment designed to unravel mercury species reactivity in superficial coastal sediments oscillating between oxic and anoxic conditions. The experimental set-up has been applied to a sediment slurry from the Arcachon Bay (France) to follow the evolution of both naturally occurring (i.e. endogenous) and isotopically enriched added mercury species (i.e. exogenous, 1??IHg and 2?1MMHg) at environmental levels. The transformation and partition between the different phases (aqueous, solid and gaseous) of the endogenous and exogenous mercury species (inorganic Hg (IHg), monomethyl Hg (MMHg), elemental Hg (Hg?) and dimethyl Hg (DMHg)) have been investigated by isotopic speciation methods throughout the experiment. The results demonstrate that the experimental approach is able to promote sediment redox oscillations and to simultaneously follow the biogeochemical fate of naturally occurring or added mercury species. Experimentally driven redox transition events were found to significantly enhance the aqueous Hg species concentrations, while the MMHg burden is not greatly affected. Indeed, during the anoxic-oxic transition, while aqueous endogenous IHg and MMHg exhibited a two-fold increase, aqueous exogenous IHg and MMHg increased 7 and 4 times, respectively. Transient increases of the net IHg methylation were recorded during the redox transitions with the largest increase of the MMHg contents (factor 1.8) observed during the oxic-anoxic transition. High resolution in situ redox experiments have not been performed up to now, therefore the developed experimental set-up provides novel insights in both the influence of redox conditions on Hg methylation/demethylation and adsorption/desorption processes and kinetics in superficial sediments.     

Evidence for millennial-scale climatic events in the sedimentary infilling of a macrotidal estuarine system,the Seine estuary (NW France)

High-resolution sedimentological and rock magnetic analyses from sediment cores collected in the Seine estuary record changes in coastal sedimentary dynamics linked to climatic variations during the late Holocene. Using AMS ¹⁴C and paleomagnetic data we present a first attempt in developing a reliable age model on macrotidal estuarine archives, with a decadal resolution. Correlations between sedimentary successions from the outer Seine estuary document the main sedimentary infilling phases of the system during the last 3000 years. Between 3000 and 1150 cal. BP sedimentary patterns reveal that sequence deposition and preservation are predominantly controlled by marine and tidal hydrodynamics while severe storm events are recorded at ca. 2700 and 1250 cal. BP in the outermost estuary. Conversely, the Medieval Warm Period (MWP; 900–1200 AD) is characterized by a drastic waning of the influence of marine hydrodynamics on sediment preservation. Pronounced episodes of Seine river floods indicate a much stronger impact of continental inputs on sedimentary patterns during this period. The onset of the Little Ice Age marks a diminishing influence of continental inputs while tidal and open marine hydrodynamics again exerted a primary control on the sedimentary evolution of the system during 1200–2003 AD. Coastal sedimentary dynamics as preserved within sedimentary successions appear to have been largely influenced by changes in storminess during the last 3000 years. We have matched the preservation of MWP Seine river floods, as revealed by sedimentological and rock magnetic proxy data, to a prolonged interval of weakened storminess in Normandy permitting the preservation of estuarine flood deposits within a context of reduced coastal erosion in northern Europe. The preservation of sedimentary successions in the Seine estuary is therefore maximal when climate conditions resembled those of the preferred low phase of the NAO on multidecadal timescales such as during 800–1200 AD (MWP). In contrast, increased removal and transport of estuarine sediments occur when winter storm activity greatly intensified over northwestern France. We report four prominent centennial-scale periods of stronger storminess, occurring with a pacing of ～1500 years, which are likely to be related to the last four Bond's Holocene cold events. Our results documenting a close link between coastal sedimentary dynamics, millennial-scale variations in Holocene climate and North Atlantic atmospheric circulation are fairly consistent with other records from Scandinavia, central Greenland and southern Europe.     

In situ adsorption of mercury,methylmercury and other elements by iron oxyhydroxides and organic matter in lake sediments

Samples of authigenic material, sediment overlying water and oxic surface sediment (0–0.5-cm depth) from a perennially oxygenated lacustrine basin were analysed to investigate which solid phases are important for binding a suite of trace elements (Ag, As, Ca, Cd, Cu, Hg, In, methylmercury (MeHg), Mg, Mo, Pb, Sb and Zn). The authigenic material, which was collected with inert Teflon sheets deployed for several years across the sediment–water interface, contained mainly poorly crystallized Fe oxyhydroxides and natural organic matter, presumably humic substances derived from the watershed. Manganese oxyhydroxides were not present in the collected authigenic material due to the slightly acidic condition (pH = 5.6) of the lake that prevents the formation and recycling of these compounds. Conditional equilibrium constants for the adsorption of cationic (K_Fe–M) and anionic (K_Fe–A) trace elements onto the authigenic Fe oxyhydroxides were estimated from their concentrations in the authigenic material and in bottom water samples. These field-derived values of K_Fe–M and K_Fe–A were compared with those predicted by the surface complexation model, using laboratory-derived intrinsic adsorption constants and the water composition at the study site. Equilibrium constants (K_POM–M) were also calculated for the adsorption of the cationic trace elements onto the humic substances contained in the diagenetic material. The field-derived values of K_POM–M were compared to those predicted by the speciation code WHAM 6 for the complexation of the trace elements by dissolved humic substances in the lake. Combining the results of the present study with those on the distributions of trace elements in the porewater and solid-phase sediments reported in previous studies at the same site, it was determined whether the trace elements bind preferentially to Fe oxyhydroxides or natural organic matter in oxic sediments. The main inferences are that the anionic trace elements As, Mo and Sb, as well as the cationic metal Pb are preferentially bound to the authigenic Fe oxyhydroxides whereas the other trace elements, and especially Hg and MeHg, are preferentially bound to the humic substances.     

Distribution and temporal variation of cadmium in the St. Lawrence river basin

Samples of raw water were collected at regular intevals at two transects in the St. Lawremce River and four of its tributaries from March to November 1991 and from April to June 1992. Water samples were analyzed for both the dissolved and the particulate phase for cadmium (Cd), organic carbon, iron and manganese. Mean dissolved Cd concentration was 10±5 ng/L and no spatial variability was observed. Higher concentrations were found during high flow periods, suggesting an uptake of cadmium by phytoplankton during summer. In addition, dissolved cadmium did not appear to be associated with either DOC, dissolved Fe or dissolved Mn. The mean particulate Cd concentration was 1.3±1.1 μg/g, with almost all stations presenting the same concentration except the Yamaska River, which had a concentration of 0.5±0.2 μg/g. Particulate Cd showed a negative correlation with suspended particulate matter and a positive correlation with particulate organic carbon and particulare manganese. Fifty-nine percent of the cadmium was found to be in the particulate phase. Partition coefficients for cadmium (Kd), organic carbon (Kc), iron (KdFe) and manganese (KdMn) were calculated for each sample. Log Kd varied from 3.9 to 5.9, with an average of 5.0±0.4. Log Kd decreased with increasing particulate, matter as did Log Kc and Log KdMn. No significant correlation was found between Log Kd and Log Kc, suggesting that the distribution of cadmium between the dissolved and the particulate phase is not influenced by the distribution of organic carbon. In contrast, positive correlations were observed between Log Kd, Log KdFe and Log KdMn. Cadmium distribution appears to be influenced by Fe and Mn distribution.     

Sequestration mechanisms and anthropogenic inputs of rhenium in sediments from Eastern Canada lakes

The concentrations of Re, as well as those of several other geochemical variables, were measured in dated sediment cores and in porewater samples from four lacustrine basins in Eastern Canada: one, perennially oxic, located 40 km from Québec City and three, seasonally anoxic, located within 25 km of non-ferrous metal smelters. The drainage basins of these lakes are uninhabited and have not been affected by human activity or wildfires. All of the depth profiles of dissolved Re indicate: higher Re concentrations in the water overlying the sediment than in the porewater; diffusion of Re across the sediment-water interface; a progressive decrease in porewater Re concentrations to reach minimum values of ∼0.5 pM within a 10-cm sediment depth interval. Modeling of these Re porewater profiles with a one-dimensional transport-reaction equation indicates that Re is removed from porewater within this depth interval. Based on thermodynamic predictions of Re speciation and of saturation states and on comparison of these predictions with sulfide porewater profiles, we infer that Re is removed from porewater by precipitation of rheniite (ReS_2(s)). The rate constant for the formation of ReS_2(s) in sediments is estimated from the modeling exercise to be 0.51 ± 0.64 × 10⁻²¹ mol cm⁻³ s⁻¹. Accumulation of sedimentary Re shows a strong authigenic component, as in anoxic marine sediments. Sharp increases in solid-phase Re during the last century are attributed to atmospheric deposition of anthropogenic Re deriving from coal burning and nearby smelter emissions.