Water temperature dynamics in a headwater forest stream: Contrasting climatic, anthropic and geological conditions create thermal mosaic of aquatic habitats.

The thermal regime of streams is a relevant driver of their ecological functioning. As this regime is presently submitted to numerous alterations (among others, impoundments, and climate change), it seems important to study both their effects and potential recovery from the latter. Thus, we investigated the surface and hyporheic water temperature along a small headwater stream with contrasting environmental contexts: forest landscape, open grassland landscape without riparian vegetation, several artificial run-of-the-river impoundments and one discharge point of a by-pass impoundment. The main objectives were to study the influence of these contrasting contexts on surface and subsurface water temperature at a local scale. Contrasting contexts were supposed to create effects on both surface and hyporheic thermal regimes at a local scale. Differences of thermal regimes between surface and hyporheos were expected, as well as between geological contexts. Sensors located at multiple stations allowed monitoring of stream and hyporheos temperature along the stream, while comparison with adjacent reference stream allowed for surface water thermal regime benchmark. Impoundments and landscapes significantly influenced stream thermal regime at a local scale (impoundments created up to +3.7°C temperature increase in average). Their effect on hyporheos thermal regime was less marked than the ones generated by solar radiation or geological features. Hyporheos thermal regime varies from stream one by temperature dynamics delay (up to 18h) and decrease (up to -7°C between surface and hyporheos temperature in average). These coupled effects create a mosaic of thermal habitats, which could be used for river biodiversity preservation and restoration.

Oxygen data assimilation for estimating micro-organism communities' parameters in river systems.

The coupling of high frequency data of water quality with physically based models of river systems is of great interest for the management of urban socio-ecosystems. One approach to exploit high frequency data is data assimilation which has received an increasing attention in the field of hydrology, but not for water quality modeling so far. We present here a first implementation of a particle filtering algorithm into a community-centered hydro-biogeochemical model to assimilate high frequency dissolved oxygen data and to estimate metabolism parameters in the Seine River system. The procedure is designed based on the results of a former sensitivity analysis of the model (Wang et al., 2018) that allows for the identification of the twelve most sensible parameters all over the year. Those parameters are both physical and related to micro-organisms (reaeration coefficient, photosynthetic parameters, growth rates, respiration rates and optimal temperature). The performances of the approach are assessed on a synthetic case study that mimics 66 km of the Seine River. Virtual dissolved oxygen data are generated using time varying parameters. This paper aims at retrieving the predefined parameters by assimilating those data. The simulated dissolved oxygen concentrations match the reference concentrations. The identification of the parameters depends on the hydrological and trophic contexts and more surprisingly on the thermal state of the river. The physical, bacterial and phytoplanktonic parameters can be retrieved properly, leading to the differentiation of two successive algal blooms by comparing the estimated posterior distribution of the optimal temperature for phytoplankton growth. Finally, photosynthetic parameters' distributions following circadian cycles during algal blooms are discussed.

Time-dependent global sensitivity analysis of the C-RIVE biogeochemical model in contrasted hydrological and trophic contexts.

Dissolved oxygen within water column is a key variable to characterize the water quality. Water quality modeling has been extensively developed for decades. However, complex biogeochemical cycles are described using a high number of parameters. Hence, parameters' uncertainty constitutes a major problem in the application of these models. Sensitivity analysis allows the identification of the most influential parameters in a model and a better understanding of the governing processes. This paper presents a time-dependent sensitivity analysis for dissolved oxygen using Morris and Sobol methods combined with a functional principal components analysis for dimension reduction. The aim of this study is to identify the most important parameters of C-RIVE model in different trophic contexts and to understand the biogeochemical functioning of river systems. The results indicate that the maintenance respiration of phytoplankton and the photosynthetic parameters (i.e. photosynthetic capacity, the maximal photosynthesis rate and light extinction coefficients) are the most influential parameters during algal blooms. When the river system becomes heterotrophic, the bacterial activities (moderate and high temperature) and the reaeration coefficients (low temperature) affect the most the dissolved oxygen concentration in the water column. An anthropogenic effect (ship navigation) on variation of dissolved oxygen concentration has been identified and the role of this anthropogenic effect evolves with hydrological and trophic conditions.

Estimation of the water quality of a large urbanized river as defined by the European WFD: what is the optimal sampling frequency?

Assessment of the quality of freshwater bodies is essential to determine the impact of human activities on water resources. The water quality status is estimated by comparing indicators with standard thresholds. Indicators are usually statistical criteria that are calculated on discrete measurements of water quality variables. If the time step of the measured time series is not sufficient to fully capture the variable's variability, the deduced indicator may not reflect the system's functioning. The goal of the present work is to assess, through a hydro-biogeochemical modeling approach, the optimal sampling frequency for an accurate estimation of 6 water quality indicators defined by the European Water Framework Directive (WFD) in a large human-impacted river, which receives large urban effluents (the Seine River across the Paris urban area). The optimal frequency depends on the sampling location and on the monitored variable. For fast varying compounds that originate from urban effluents, such as PO[Formula: see text], NH[Formula: see text] and NO[Formula: see text], a sampling time step of one week or less is necessary. To be able to reflect the highly transient character of bloom events, chl a concentrations also require a short monitoring time step. On the contrary, for variables that exert high seasonal variability, as NO[Formula: see text] and O 2, monthly sampling can be sufficient for an accurate estimation of WFD indicators in locations far enough from major effluents. Integrative water quality variables, such as O 2, can be highly sensitive to hydrological conditions. It would therefore be relevant to assess the quality of water bodies at a seasonal scale rather than at annual or pluri-annual scales. This study points out the possibility to develop smarter monitoring systems by coupling both time adaptative automated monitoring networks and modeling tools used as spatio-temporal interpolators.

Estimating River Conductance from Prior Information to Improve Surface-Subsurface Model Calibration.

Most groundwater models simulate stream-aquifer interactions with a head-dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill-posed and individual model parameters are likely to be poorly constrained. Ill-posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface-subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water-groundwater use.

Modelling the fate of nitrite in an urbanized river using experimentally obtained nitrifier growth parameters.

Maintaining low nitrite concentrations in aquatic systems is a major issue for stakeholders due to nitrite's high toxicity for living species. This study reports on a cost-effective and realistic approach to study nitrite dynamics and improve its modelling in human-impacted river systems. The implementation of different nitrifying biomasses to model riverine communities and waste water treatment plant (WWTP)-related communities enabled us to assess the impact of a major WWTP effluent on in-river nitrification dynamics. The optimal kinetic parameters and biomasses of the different nitrifying communities were determined and validated by coupling laboratory experiments and modelling. This approach was carried out in the Seine River, as an example of a large human-impacted river with high nitrite concentrations. The simulation of nitrite fate was performed at a high spatial and temporal resolution (Δt = 10 min, dx¯ = 500 m) including water and sediment layers along a 220 km stretch of the Seine River for a 6-year period (2007-2012). The model outputs were in good agreement with the peak of nitrite downstream the WWTP as well as its slow decrease towards the estuary. Nitrite persistence between the WWTP and the estuary was mostly explained by similar production and consumption rates of nitrite in both water and sediment layers. The sediment layer constituted a significant source of nitrite, especially during high river discharges (0.1-0.4 mgN h(-1) m(-2)). This points out how essential it is to represent the benthic layer in river water quality models, since it can constitute a source of nitrite to the water-column. As a consequence of anthropogenic emissions and in-river processes, nitrite fluxes to the estuary were significant and varied from 4.1 to 5.5 TN d(-1) in low and high water discharge conditions, respectively, over the 2007-2012 period. This study provides a methodology that can be applied to any anthropized river to realistically parametrize autochthonous and WWTP-related nitrifier communities and simulate nitrite dynamics. Based on simulation analysis, it is shown that high spatio-temporal resolution hydro-ecological models are efficient to 1) estimate water quality criteria and 2) forecast the effect of future management strategies. Process-based simulations constitute essential tools to complete our understanding of nutrient cycling, and to decrease monitoring costs in the context of water quality and eutrophication management in river ecosystems.

Pluri-annual sediment budget in a navigated river system: the Seine River (France).

This study aims at quantifying pluri-annual Total Suspended Matter (TSM) budgets, and notably the share of river navigation in total re-suspension at a long-term scale, in the Seine River along a 225 km stretch including the Paris area. Erosion is calculated based on the transport capacity concept with an additional term for the energy dissipated by river navigation. Erosion processes are fitted for the 2007-2011 period based on i) a hydrological typology of sedimentary processes and ii) a simultaneous calibration and retrospective validation procedure. The correlation between observed and simulated TSM concentrations is higher than 0.91 at all monitoring stations. A variographic analysis points out the possible sources of discrepancies between the variabilities of observed and simulated TSM concentrations at three time scales: sub-weekly, monthly and seasonally. Most of the error on the variability of simulated concentrations concerns sub-weekly variations and may be caused by boundary condition estimates rather than modeling of in-river processes. Once fitted, the model permits to quantify that only a small fraction of the TSM flux sediments onto the river bed (<0.3‰). The river navigation contributes significantly to TSM re-suspension in average (about 20%) and during low flow periods (over 50%). Given the significant impact that sedimentary processes can have on the water quality of rivers, these results highlight the importance of taking into account river navigation as a source of re-suspension, especially during low flow periods when biogeochemical processes are the most intense.

Modelling the fate of nonylphenolic compounds in the Seine River--part 2: assessing the impact of global change on daily concentrations.

This study aims at modelling the daily concentrations of nonylphenolic compounds such as 4-nonylphenol (4-NP), nonylphenol monoethoxylate (NP1EO) and nonylphenoxy acetic acid (NP1EC) within the Seine River downstream of Paris City for over a year, firstly in the present state (year 2010) and for years 2050 and 2100 in order to assess the consequences of global change on the fate of nonylphenolic compounds in the Seine river. Concentrations were first simulated for the year 2010 and compared to monthly measured values downstream of Paris. To achieve this goal, the hydrodynamic and biogeochemical model, ProSe, was updated to simulate the fate of 4-NP, NP1EO and NP1EC. The Seine upstream and Oise River (tributaries of the Seine River) concentrations are estimated according to concentrations-flow relationships. For Seine Aval wastewater treatment plant (SA-WWTP), the concentrations are considered constant and the median values of 11 campaigns are used. The biodegradation kinetics of 4-NP, NP1EO and NP1EC in the Seine River were deduced from the results of the companion paper. The Nash-Sutcliffe coefficient indicates a good efficiency to simulate the concentrations of 4-NP, NP1EC and NP1EO over an entire year. Eight scenarios were built to forecast the impacts of global warming (flow decrease), population growth (SA-WWTP flow increase) and optimisation of wastewater treatment (improvement of the quality of effluents) on annual concentrations of 4-NP, NP1EO and NP1EC at Meulan by 2050 and 2100. As a result, global warming and population growth may increase the concentrations of 4-NP, NP1EC and NP1EO, especially during low-flow conditions, while the optimisation of wastewater treatment is an efficient solution to balance the global change by reducing WWTP outflows.

Modelling the fate of nonylphenolic compounds in the Seine River--part 1: determination of in-situ attenuation rate constants.

Assessing the fate of endocrine disrupting compounds (EDCs) in the environment is currently a key issue for determining their impacts on aquatic ecosystems. The 4-nonylphenol (4-NP) is a well known EDC and results from the biodegradation of surfactant nonylphenol ethoxylates (NPnEOs). Fate mechanisms of NPnEO are well documented but their rate constants have been mainly determined through laboratory experiments. This study aims at evaluating the in-situ fate of 4-NP, nonylphenol monoethoxylate (NP1EO) and nonylphenolic acetic acid (NP1EC). Two sampling campaigns were carried out on the Seine River in July and September 2011, along a 28km-transect downstream Paris City. The field measurements are used for the calibration of a sub-model of NPnEO fate, included into a hydro-ecological model of the Seine River (ProSe). The timing of the sampling is based on the Seine River velocity in order to follow a volume of water. Based on our results, in-situ attenuation rate constants of 4-NP, NP1EO and NP1EC for both campaigns are evaluated. These rate constants vary greatly. Although the attenuation rate constants in July are especially high (higher than 1d(-1)), those obtained in September are lower and consistent with the literature. This is probably due to the biogeochemical conditions in the Seine River. Indeed, the July sampling campaign took place at the end of an algal bloom leading to an unusual bacterial biomass while the September campaign was carried out during common biogeochemical status. Finally, the uncertainties on measurements and on the calibration parameters are estimated through a sensitivity analysis. This study provides relevant information regarding the fate of biodegradable pollutants in an aquatic environment by coupling field measurements and a biogeochemical model. Such data may be very helpful in the future to better understand the fate of nonylphenolic compounds or any other pollutants at the basin scale.

Modeling nitrate fluxes at the catchment scale using the integrated tool CAWAQS.

Nitrates fluxes in the Grand Morin basin (1200 km(2)), that is subjected to intense agricultural pressure, are considered using in-stream observations (around 250 sampling days over 5 years) and physically based simulations using the CAWAQS model (CAtchment WAter Quality Simulator). In-stream nitrate concentration averaged 6 mg N L(-1), increasing by approximately 0.2 mg N L(-1) yr(-1) around this value (period 1991-1996). Our results show that, over the period of 1991-1996, the differences between in-stream observed nitrate concentrations and simulated nitrate concentrations result from nitrate losses at the basin scale. These losses are due to denitrification by transfer through wetlands, alluvial plains, the hyporheic zone, and by benthic processes in rivers. A mean annual mass balance at the basin scale indicates that 40% of the infiltration flux (3360 kg N km(-2) yr(-1)) is removed from the system via the river network, 40% is stored in aquifers and 20% is lost by denitrification (period 1991-1996).

Primary production in headwater streams of the Seine basin: the Grand Morin river case study.

Periphytic biomass has an important influence on the water quality of many shallow streams. The purpose of this paper is to synthesize the knowledge obtained on periphyton during the PIREN Seine research program. Periphyton was sampled using chl a measurements by acetone extraction and oxygen measurements with microelectrodes. The experiments reveal the presence of an important fixed biomass ranging between 123 and 850 mgchl a m(-2) and the mean gross production (photosynthesis) is shown to range between 180 and 315 mgC m(-2) h(-1). An independent approach was performed using the ProSe model, which simulates transport and biogeochemical processes in 22 km of the Grand Morin stream. A strong agreement between in situ measurements and the model results was obtained. The gross production obtained using ProSe is 220 mgC m(-2) h(-1) for the periphyton, which matches the experimental data. Although the net photosynthetic activity of the phytoplankton (0.84 gC gC(-1) d(-1)) is higher than the periphytic one (0.33 gC gC(-1) d(-1)), the absolute periphytic activity is greater since the mean biomass (3.4 gC m(-)(2)) is 10 times higher than the phytoplanktonic one (0.3 gC m(-2)), due to the short residence time of the water body (1.5d).

Modelling the impacts of Combined Sewer Overflows on the river Seine water quality.

To achieve the objectives of the European Water Framework Directive (EWFD), the Seine basin Water Authority has constructed a number of prospective scenarios forecasting the impact of planned investments in water quality. Paris and its suburbs were given special attention because of their impact on the river Seine. Paris sewer system and overflow control is of major concern in future management plans. The composition and fate of the urban effluents have been characterized through numerous in situ samplings, laboratory experiments and modelling studies. The PROSE model was especially designed to simulate the impact on the river of both permanent dry-weather effluents and of highly transient Combined Sewer Overflow (CSO). It was also used to represent the impact of Paris at large spatial and temporal scales. In addition to immediate effects on oxygen levels, heavy particulate organic matter loads that settle downstream of the outlets contribute to permanent oxygen consumption. Until the late 90s, the 50 km long reach of the Seine inside Paris was permanently affected by high oxygen consumption accounting for 112% of the flux upstream of the city. 20% of this demand resulted from CSO. However, the oxygenation of the system is strong due to high phytoplankton activity. As expected, the model results predict a reduction of both permanent dry-weather effluents and CSOs in the future that will greatly improve the oxygen levels (concentrations higher than 7.3 mgO(2) L(-1), 90% of the time instead of 4.0 mgO(2) L(-1) in the late 90s). The main conclusion is that, given the spatial and temporal extent of the impact of many CSOs, water quality models should take into account the CSOs in order to be reliable.

Assessment of nitrate pollution in the Grand Morin aquifers (France): combined use of geostatistics and physically based modeling.

The objective of this work is to combine several approaches to better understand nitrate fate in the Grand Morin aquifers (2700 km(2)), part of the Seine basin. cawaqs results from the coupling of the hydrogeological model newsam with the hydrodynamic and biogeochemical model of river ProSe. cawaqs is coupled with the agronomic model Stics in order to simulate nitrate migration in basins. First, kriging provides a satisfactory representation of aquifer nitrate contamination from local observations, to set initial conditions for the physically based model. Then associated confidence intervals, derived from data using geostatistics, are used to validate cawaqs results. Results and evaluation obtained from the combination of these approaches are given (period 1977-1988). Then cawaqs is used to simulate nitrate fate for a 20-year period (1977-1996). The mean nitrate concentrations increase in aquifers is 0.09 mgN L(-1)yr(-1), resulting from an average infiltration flux of 3500 kgN.km(-2)yr(-1).