Brian Four

Evelyne Arce

Michaël Danger

Juliette Gaillard

Marielle Thomas

Damien Banas

Catchment land use-dependent effects of barrage fishponds on the functioning of headwater streams

Keywords: Fishpond, headwater stream, land use litter breakdown, macroinvertebrate, hyphomycete, stream ecosystem functioning

come   L'archive ouverte pluridisciplinaire

Introduction __________________________________________________________________

Headwater streams, which may represent up to 80% of the total stream length [START_REF] Macdonald | Influence of headwater streams on downstream reaches in forested areas[END_REF], are essential for ensuring high water quality and good status of downstream ecosystems [START_REF] Alexander | The role of headwater streams in downstream water quality[END_REF][START_REF] Gomi | Understanding processes and downstream linkages of headwater systems[END_REF]. The functioning of streams, especially headwater streams, is closely related to their associated catchments [START_REF] Fisher | Energy flow in bear brook, New Hampshire: an integrative approach to stream ecosystem metabolism[END_REF][START_REF] Wallace | Multiple trophic levels of a forest stream linked to terrestrial litter inputs[END_REF]. Human activities in the catchments can significantly alter headwater streams (e.g. [START_REF] Bernot | Inter-regional comparison of land-use effects on stream metabolism[END_REF]. Intensive agriculture, in particular, has been identified as a major stressor for freshwater ecosystems [START_REF] Vörösmarty | Global threats to human water security and river biodi-versity[END_REF]). As shown by several studies, conventional crop production markedly affects surface water with, for example, increases in nutrient concentrations [START_REF] Hughes | Updating previous estimates of the load and source apportionment of nitrogen to waters in the UK[END_REF][START_REF] White | The sources of phosphorus in the waters of Great Britain[END_REF], sediment discharges [START_REF] Collins | Assessing the likelihood of catchments across England and Wales meeting Bgood ecological status^ due to sediment contributions from agricultural sources[END_REF] and pesticide concentrations [START_REF] Kreuger | Pesticides in stream water within an agricultural catchment in southern Swenden, 1990-1996[END_REF]. Streams draining agricultural land are also often characterized by loss of hydromorphological features of natural streams (e.g. meanders) and loss of habitat heterogeneity due to dredging and channelization [START_REF] Negishi | Effects of channelisation on stream habitat in relation to a spate and flow refugia for macroinver-tebrates in northern Japan[END_REF][START_REF] Pedersen | Effects of channelisation, riparian structure and catchment area on physical habitats in small lowland streams[END_REF]). These physical and chemical alterations may be a cause of biodiversity loss in streams [START_REF] Liess | The footprint of pesticide stress in communities-species traits reveal community effects of toxicants[END_REF][START_REF] Rasmussen | Pesticide effects on the structure and function of stream ecosystems[END_REF]) and affect fundamental ecological processes such as leaf litter decomposition and gross primary production [START_REF] Peters | Review on the effects of toxicants on freshwater ecosystem functions[END_REF][START_REF] Rasmussen | Impacts of pesticides and natural stressors on leaf litter decomposition in agricultural streams[END_REF][START_REF] Schäfer | Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe[END_REF][START_REF] Robinson | Nutrient addition accelerates leaf breakdown in an alpine Springbrook[END_REF].

Beyond these quite well-known effects of pollution, far less attention has been paid to other aspects of river integrity. Yet, due to their small sizes, numerous headwater streams suffer strong hydromorphological impact effects [START_REF] Elosegi | Effects of hydromorphological impacts on river ecosystem functioning: a review and suggestions for assessing ecological impacts[END_REF]). Yet, hydromorphology, defined as the complex interaction between water flow and channel form, is an essential status component for streams [START_REF] Poole | Stream hydrogeomorphology as a physical science basis for advances in stream ecology[END_REF]. Hydromorphology can have considerable effects on water quality, community structure and stream ecosystem functioning [START_REF] Elosegi | Effects of hydromorphological integrity on biodiversity and functioning of river ecosystems[END_REF]). Among the various catchment uses that could affect river continuum and significantly alter water flow, barrage fishponds (i.e. drainable ponds used as extensive fish production systems), made by building a small dam on headwater streams, represent a common hydromorphological alteration. Fishponds account for a large proportion of all surface water bodies in France (1200 km2), the Czech Republic (410 km2) and Germany (420 km2) (Le [START_REF] Quéré | La pisciculture d'étangs française[END_REF][START_REF] Pokorný | The restoration offish ponds in agricultural landscapes[END_REF]. In France, it is estimated that there are over 251,000 ponds [START_REF] Bartout | L'inventaire des plans d'eau français :outil d'une meilleure gestion des eaux de surface[END_REF], most of which are located on first-order streams. Although the ecosystem services provided by fishpond systems are increasingly recognized [START_REF] Bekefi | Multifunctional pond fish farmers in Hungary[END_REF][START_REF] Blayac | Perceptions of the services provided by pond fish farming in Lorraine (France)[END_REF][START_REF] Mathé | Local knowledge of pond fish-farming ecosystem services: management implications of stakeholders' perceptions in three different contexts (Brazil, France and Indonesia)[END_REF], the presence of dams can strongly affect the biodiversity and functioning of streams [START_REF] Bunn | Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity[END_REF][START_REF] Elosegi | Effects of hydromorphological impacts on river ecosystem functioning: a review and suggestions for assessing ecological impacts[END_REF][START_REF] González | Small dams decrease leaf litter breakdown rates in Mediterranean mountain streams[END_REF]. Due to increases in water retention time, sediment retention and particulate organic matter production, to changes in water chemistry or to the partial or total reduction of the aquatic organism migration [START_REF] Pringle | Exploring how disturbance is transmitted upstream: going against the the flow[END_REF][START_REF] Kunz | Impact of a large tropical reservoir on riverine transport of sediment, carbon, and nutrients to downstream wetlands[END_REF]Gonzales et al. 2013;[START_REF] Colas | Synergistic impacts of sediment contamination and dam presence on river functioning[END_REF], the effects of dams have most often been seen as negative for stream ecosystems. Yet, in the context of agricultural landscapes, some studies have shown that the presence of these fish production systems in streams can influence positively the flux of water, decreasing the suspended matter, and pesticide and nutrient contents of the downstream water [START_REF] Banas | Flux de matière en étangs piscicoles extensifs :Rétention, Sédimentation, Exportation[END_REF][START_REF] Banas | Discharge of sediments, nitrogen (N) and phosphorus (P) during the emptying of extensive fishponds: effect of rain-fall and management practices[END_REF][START_REF] Gaillard | Rôle des étangs de barrage à vocation piscicole dans la dynamique des micropolluants en têtes de bassins versants[END_REF]Gaillard et al. 2016a) except during the draining period when significant release of nutrients can be observed [START_REF] Banas | Assessment of sediment concentration and nutrients loads in effluents drained from extensivelymanaged fishponds in France[END_REF]). More recently, Gaillard et al. (2016b) measured a significant reduction in pesticide peak concentrations downstream from barrage fishponds (from 49 to 99%, depending on the molecule). From this, it could be said that barrage fishponds improve water quality for streams in agricultural catchment areas. As a Headwater streams, which may represent up to 80% of the total stream length (MacDonald and Coe 2007), are essential for ensuring high water quality and good status of downstream ecosystems [START_REF] Alexander | The role of headwater streams in downstream water quality[END_REF][START_REF] Gomi | Understanding processes and downstream linkages of headwater systems[END_REF]. The functioning of streams, especially headwater streams, is closely related to their associated catchments [START_REF] Fisher | Energy flow in bear brook, New Hampshire: an integrative approach to stream ecosystem metabolism[END_REF][START_REF] Wallace | Multiple trophic levels of a forest stream linked to terrestrial litter inputs[END_REF]. Human activities in the catchments can significantly alter headwater streams (e.g. [START_REF] Bernot | Inter-regional comparison of land-use effects on stream metabolism[END_REF]. Intensive agriculture, in particular, has been identified as a major stressor for freshwater ecosystems [START_REF] Vörösmarty | Global threats to human water security and river biodi-versity[END_REF]). As shown by several studies, conventional crop production markedly affects surface water with, for example, increases in nutrient concentrations [START_REF] Hughes | Updating previous estimates of the load and source apportionment of nitrogen to waters in the UK[END_REF][START_REF] White | The sources of phosphorus in the waters of Great Britain[END_REF], sediment discharges [START_REF] Collins | Assessing the likelihood of catchments across England and Wales meeting Bgood ecological status^ due to sediment contributions from agricultural sources[END_REF] and pesticide concentrations [START_REF] Kreuger | Pesticides in stream water within an agricultural catchment in southern Swenden, 1990-1996[END_REF]. Streams draining agricultural land are also often characterized by loss of hydromorphological features of natural streams (e.g. meanders) and loss of habitat heterogeneity due to dredging and channelization [START_REF] Negishi | Effects of channelisation on stream habitat in relation to a spate and flow refugia for macroinver-tebrates in northern Japan[END_REF][START_REF] Pedersen | Effects of channelisation, riparian structure and catchment area on physical habitats in small lowland streams[END_REF]). These physical and chemical alterations may be a cause of biodiversity loss in streams [START_REF] Liess | The footprint of pesticide stress in communities-species traits reveal community effects of toxicants[END_REF][START_REF] Rasmussen | Pesticide effects on the structure and function of stream ecosystems[END_REF]) and affect fundamental ecological processes such as leaf litter decomposition and gross primary production [START_REF] Peters | Review on the effects of toxicants on freshwater ecosystem functions[END_REF][START_REF] Rasmussen | Impacts of pesticides and natural stressors on leaf litter decomposition in agricultural streams[END_REF][START_REF] Schäfer | Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe[END_REF][START_REF] Robinson | Nutrient addition accelerates leaf breakdown in an alpine Springbrook[END_REF].

Beyond these quite well-known effects of pollution, far less attention has been paid to other aspects of river integrity. Yet, due to their small sizes, numerous headwater streams suffer strong hydromorphological impact effects [START_REF] Elosegi | Effects of hydromorphological impacts on river ecosystem functioning: a review and suggestions for assessing ecological impacts[END_REF]). Yet, hydromorphology, defined as the complex interaction between water flow and channel form, is an essential status component for streams [START_REF] Poole | Stream hydrogeomorphology as a physical science basis for advances in stream ecology[END_REF]. Hydromorphology can have considerable effects on water quality, community structure and stream ecosystem functioning [START_REF] Elosegi | Effects of hydromorphological integrity on biodiversity and functioning of river ecosystems[END_REF]). Among the various catchment uses that could affect river continuum and significantly alter water flow, barrage fishponds (i.e. drainable ponds used as extensive fish production systems), made by building a small dam on headwater streams, represent a common hydromorphological alteration. Fishponds account for a large proportion of all surface water bodies in France (1200 km2), the Czech Republic (410 km2) and Germany (420 km2) (Le [START_REF] Quéré | La pisciculture d'étangs française[END_REF][START_REF] Pokorný | The restoration offish ponds in agricultural landscapes[END_REF]. In France, it is estimated that there are over 251,000 ponds [START_REF] Bartout | L'inventaire des plans d'eau français :outil d'une meilleure gestion des eaux de surface[END_REF], most of which are located on first-order streams. Although the ecosystem services provided by fishpond systems are increasingly recognized [START_REF] Bekefi | Multifunctional pond fish farmers in Hungary[END_REF][START_REF] Blayac | Perceptions of the services provided by pond fish farming in Lorraine (France)[END_REF][START_REF] Mathé | Local knowledge of pond fish-farming ecosystem services: management implications of stakeholders' perceptions in three different contexts (Brazil, France and Indonesia)[END_REF], the presence of dams can strongly affect the biodiversity and functioning of streams [START_REF] Bunn | Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity[END_REF][START_REF] Elosegi | Effects of hydromorphological impacts on river ecosystem functioning: a review and suggestions for assessing ecological impacts[END_REF][START_REF] González | Small dams decrease leaf litter breakdown rates in Mediterranean mountain streams[END_REF]. Due to increases in water retention time, sediment retention and particulate organic matter production, to changes in water chemistry or to the partial or total reduction of the aquatic organism migration [START_REF] Pringle | Exploring how disturbance is transmitted upstream: going against the the flow[END_REF][START_REF] Kunz | Impact of a large tropical reservoir on riverine transport of sediment, carbon, and nutrients to downstream wetlands[END_REF]Gonzales et al. 2013;[START_REF] Colas | Synergistic impacts of sediment contamination and dam presence on river functioning[END_REF], the effects of dams have most often been seen as negative for stream ecosystems. Yet, in the context of agricultural landscapes, some studies have shown that the presence of these fish production systems in streams can influence positively the flux of water, decreasing the suspended matter, and pesticide and nutrient contents of the downstream water [START_REF] Banas | Flux de matière en étangs piscicoles extensifs :Rétention, Sédimentation, Exportation[END_REF][START_REF] Banas | Discharge of sediments, nitrogen (N) and phosphorus (P) during the emptying of extensive fishponds: effect of rain-fall and management practices[END_REF][START_REF] Gaillard | Rôle des étangs de barrage à vocation piscicole dans la dynamique des micropolluants en têtes de bassins versants[END_REF]Gaillard et al. 2016a) except during the draining period when significant release of nutrients can be observed [START_REF] Banas | Assessment of sediment concentration and nutrients loads in effluents drained from extensivelymanaged fishponds in France[END_REF]). More recently, Gaillard et al. (2016b) measured a significant reduction in pesticide peak concentrations downstream from barrage fishponds (from 49 to 99%, depending on the molecule). From this, it could be said that barrage fishponds improve water quality for streams in agricultural catchment areas. As a consequence, it could be expected that the influence of fishponds on the functioning of headwater streams might be highly dependent on the catchment land use.

To investigate headwater stream status, water or various biological indices mainly based on the sensitivity of some aquatic taxa (fish, diatoms, macroinvertebrates and macrophytes) to multiple stressors [START_REF] Birk | Three hundred ways to assess Europe's surface waters: an almost complete overview of biological methods to implement the Water Framework Directive[END_REF]. However, these indices may be considered unsatisfying due to major inconsistencies with the Water Framework Directive (WFD, European Union 2000; [START_REF] Roche | Les enjeux de recherche liés à la directive-cadre européenne sur l'eau[END_REF][START_REF] Mondy | A new macroinvertebrate-based multimetric index (I2M2) to evaluate ecological quality of French wadeable streams fulfilling the WFD demands: a taxonomical and trait approach[END_REF]. Recently, ecologists have emphasized the need to use both the structural and functional components of ecosystem communities to assess the ecological status of streams and the effects of various environmental stressors [START_REF] Clarke | Macroinvertebrate diversity in headwater streams: a review[END_REF][START_REF] Graça | The role of invertebrates on leaf litter decomposition in streamsareview[END_REF][START_REF] Maltby | Spray drift of pesticides and stream macroin-vertebrates: experimental evidence of impacts and effectiveness of mitigation measures[END_REF][START_REF] Tachet | Invertébrés d'eau douce[END_REF]. In order to understand in greater depth the influence that these stressors have on stream ecosystem functioning, leaf litter decomposition has been proposed as an integrative indicator of headwater stream functioning and developed as such by several authors [START_REF] Gessner | A case for using litter breakdown to assess functional stream integrity[END_REF][START_REF] Graça | The role of invertebrates on leaf litter decomposition in streamsareview[END_REF][START_REF] Graça | Methods to study litter decomposition: a practical guide[END_REF][START_REF] Graça | A conceptual model of litter breakdown in low order streams: litter breakdown in low order streams[END_REF][START_REF] Tank | A review of allochthonous organic matter dynamics and metabolism in streams[END_REF]. This parameter is a key ecosystem level process that integrates the activity of both microbial decomposers (mostly aquatic hyphomycetes) and aquatic macroinvertebrates. It has been successfully applied to the ecological monitoring of streams (e.g. [START_REF] Woodward | Continental-scale effects of nutrient pollution on stream ecosystem functioning[END_REF].

In order to investigate the influence of barrage fishponds on headwater streams along with the catchment land use, we selected five fishponds located on five headwater streams (two in agricultural landscapes and three in forested landscapes) and monitored (i) leaf litter processing (evaluated by the litter bag technique), (ii) fungal biomass and (iii) structural and functional metrics of macroinvertebrate communities associated with litter bags. The barrage fishponds were distributed along a gradient of land uses in the various catchments from extensive forest management to intensive cereal production. Since land use is well known to highly influence water quality of streams, we hypothesized that barrage fishponds could have differential effects on the functioning of headwater streams depending on the catchment land use.

2.

Material and Methods __________________________________________________________

Study area

he study was conducted in the Lorraine region (north-eastern France) located at the extreme east end of the Paris sedimentary drainage basin near the Vosges mountains (Fig. 1). The climate is semi-continental with an average annual air temperature of 10.7 °C, an average minimum temperature of about 2.2 °C in January and an average maximum temperature of 19.7 °C in July (30-year average, Château-Salins, Meteo France 2011). The average annual precipitation is about 800 mm, and it is well distributed throughout the year. During the investigation period (21 January to 17 March 2014), the average air temperature was 5.1 °C, which falls within the temperature range of relatively warm years. The cumulated daily precipitation during the whole investigation period was 92.2 mm, which falls within the precipitation range of relatively dry years. This region, with a total surface area of 23,547 km2, is characterized by a large cover of mixed deciduous forests (36%, where Quercus, Carpinus and Fagus are the dominant genera) and agricultural lands (27% arable lands and 20% pastures). Agricultural lands are often managed intensively using pesticides and chemical fertilizers [START_REF] Joulin | Les produits phytosanitaires utilisés par l'agriculture lorraine[END_REF]. Moreover, the Lorraine region, with a fish production of 1100 t (CRAL 2005), ranks third in freshwater fish production in France thanks to a high density of fishponds covering an area of about 7000 ha (Le [START_REF] Quéré | La pisciculture d'étangs française[END_REF]. Most of these fishponds were created along a headwater stream continuum in the Middle Ages and are commonly exploited extensively for fish production. Numerous fishponds are often exposed to pesticide pression as a result of agricultural land use [START_REF] Lazartigues | Pesticide pressure and fish farming in barrage pond in Northeastern France part I: site characterization and water quality[END_REF](Lazartigues et al. , 2013a;;Lazartigues et al. 2013b).

Study sites

To investigate the influence of fishponds on headwater stream functioning, we selected five barrage fishponds located on first-order headwater streams (according to the classification of [START_REF] Strahler | Quantitative analysis of watershed geomorphology[END_REF]) in a limited geographical area (60 km2) with homogenous geology (Triassic sedimentary deposits). We studied potential differences in leaf litter processing, fungal biomass and structural and functional metrics of macroinver-tebrate communities between sites located upstream (Up) and downstream (Down) from the fishponds. At each sampling point (i.e., Up and Down of the five fishponds), a 30-m-long reach was selected at approximately 50 m from the ponds to avoid the drawdown zone upstream as well as the direct effect of water fall from dams downstream. Selection of the fishponds was based on the following criteria: (i) major water inflows brought by one main tributary to the pond and (ü) similar land uses between upstream and downstream sampling points in each catchment area. In order to carry this field study, we selected sites representative of the field context where ponds are located either on agricultural catchments (but always with few percent of other land cover) or in forest (but always with few percent of agricultural area). Classification of our sites into two land use categories was done with an accurate protocol, on the basis of three different criteria. To define a catchment in the category 'Forest', it is necessary to meet the three following criteria:

1. River basin must be dominated by native deciduous forest.

2. Stream running through a significant distance of forest 3. Main land cover surrounding the sampling points must be also dominated by native tree species. In this case, organic matter (OM) available as food source for trophic webs is dominated by tree litters (i.e. leaves and branches, as observed for the three forested sites).

If these three criteria are met, the catchment was ranged in the category Forest. If not, the catchment was considered in the category 'Agriculture'. By applying these criteria, we categorized there three fishponds into 'Forest sites' (noted F1, F2 and F3) and two into 'Agriculture sites' (noted A1 and A2). In the forest catchments, the dominant tree species were Carpinus betulus, Fagus sylvatica and Quercus sp. Among arable lands, the cultivated surface areas and the proportion of crop variety varied annually as follows. The farmers mainly grew rapeseed (Brassica napus), wheat (Triticum aestivum) and barley (Hordeum vulgare) in a 3-year cultural rotation as commonly practised in the region [START_REF] Xiao | Modeling the spatial distribution of crop sequences at a large regional scale using landcover survey data: a case from France[END_REF]) (Table 1).

Fishpond areas were 4.4, 31.6, 4.9, 1.5 and 7.2 ha for A1, A2, F1, F2 and F3, respectively (Table 1). The volumes of these ponds were estimated at 44000, 568,000, 49,000,16,000 and 108,000 m3 [START_REF] Banas | Flux de matière en étangs piscicoles extensifs :Rétention, Sédimentation, Exportation[END_REF][START_REF] Gaillard | Rôle des étangs de barrage à vocation piscicole dans la dynamique des micropolluants en têtes de bassins versants[END_REF]Lazartigues 2013). The ratios between the pond and catchment surface areas, which can serve as a proxy for hydraulic retention time ØT), were 1:20, 1:11, 1:20, 1:42 and 1:21 for A1, A2, F1, F2 and F3, respectively (Table 1). Highest HRT was expected for A2 and lowest HRT was expected for F2. In the Lorraine region (n =105 sites), the median pond to surface ratio is 1:31 which means the studied sites were typical of ponds in the region (unpublished data).

All the fishponds under investigation are extensively managed for polyculture, including species with different diets and behaviour (Cyprinus carpio, Rutilus rutilus, Tinca tinca, Perca fluviatilis, Esox lucius and/or Sander lucioperca). Management operations include a 2-year cycle production with three steps. The first corresponds to filling of the ponds due to water inputs from small streams, surface water runoff and precipitations. The following step corresponds to a pseudo-balance phase during which fish are stocked and grown in ponds. It must be emphasized that additional feeding and direct use of fertilizers or pharmaceutical compounds are not added in the five studied fishponds, as is usually the case. Finally, the third step is drainage in late autumn or early spring. All of the water is discharged downstream to allow fish harvest every 1 or 2 years.

The riparian vegetation, basin surface area and land use characteristics of the ten sampling points are given in Table 1.

Water physical and chemical parameters

The experiment was carried out over 8 weeks (21 January 2014 corresponding to day 0 to 17 March 2014 corresponding to day 55). Water quality was monitored at all sampling points (upstream and downstream from the five fishponds). Water temperature was measured throughout the entire experimental period with temperature loggers ata 15min time step (HOBO Pendant data logger UA-001-64). Turbidity, pH, dissolved oxygen, conductivity and oxidationreduction potential were measured (n = 6-8) throughout the decomposition experiment (PONSEL ODEON X line tools). Water samples were collected at each litter sampling date (day 0, 14, 34 and 55). Samples were taken to the laboratory in a cool box for nitrate, nitrite, ammonium and soluble reactive phosphorus (SRP) analyses. Within 24 h, water samples were filtered (Whatman GF/F) and nutrient measurements were conducted on the filtered water. Nitrate concentrations were determined using the hydrazine reduction method (NF ISO 15923-1, 2014), nitrite concentrations using the sulphanilamide method (NF ISO 15923-1, 2014), ammonium concentrations using the indophenol blue spectrophotometric method (NF T 90-015-2, 2000) and SRP concentrations using the molybdate method (NF ISO 159231, 2014). All the water physical and chemical results are reported in Table 2.

Leaf decomposition

In autumn 2013, leaves were collected just after abscission from the same stand of maple trees (Acer pseudoplatanus) in the Vosges Mountains using a net hung between the trees 1 m above the ground. Maple leaves were chosen because they are supposed to have a medium to fast litter breakdown rate [START_REF] Lecerf | Perturbations anthropiques et fonctionnement écologique des cours d'eau de tête de bassin :étude du processus de décomposition des litières[END_REF][START_REF] Petersen | Leaf processing in a woodland stream[END_REF], an important property for our study because study sites can be temporarily dry out. Thus, the whole experiment must be conducted between the end of autumn and before summer drying. Like in these small streams a large part of the shredder communities implicated in the litter decay have annual cycles with larvae aquatic stage during the winter and adult terrestrial stage during the summer, we have chosen to conduct the study at the end of their aquatic stage and just before emergencies (i.e. from January to March 2014).

Leaf litter was air-dried in a room and stored in the dark under dry conditions before being used. Leaf petioles were removed, and then leaves were weighed in batches of 3 ± 0.02 g, moistened with deionized water, placed in 40 fine and 120 coarse mesh bags (0.5 and 10 mm mesh size, respectively, following the method described by [START_REF] Gessner | A case for using litter breakdown to assess functional stream integrity[END_REF]) and deployed at the ten sampling points on day 0. Coarse mesh bags allowed for shredder colonization whereas fine mesh bags excluded them and therefore only reflecting microbial activity (largely by microfungi) and leaching. Coarse mesh bags were used to determine the total leaf breakdown rate at each site, and fine mesh bags were used to determine fungal biomass.

Four coarse mesh bags were collected at each sampling point after 14, 34 and 55 days of exposure, and the fine mesh bags were collected after 34 days of exposure. All leaf bags were removed from the streams using 0.5 mm (mesh size) sieves to avoid invertebrate loss, stored individually in ziplock bags and returned to the laboratory in a cool box for processing. In the laboratory, leaves were carefully washed on a 0.5 mm sieve to remove sediment and exogenous organic matter and to collect invertebrates. Then, samples from coarse mesh bags were oven-dried at 105 °C for 48 h and weighed to the nearest 0.01 g and ground. In order to minimize the bias from sediment contamination, results were expressed in ash free dry mass (AFDM). Leaf AFDM was determined on 500 mg sub-samples of ground leaves ashed in a muffle furnace (550 °C for 4 h). As described by [START_REF] Abelho | From Litterfall to breakdown in streams: a review[END_REF], leaching of soluble compounds can account for more than 40% of initial dry mass loss. In order to minimize the effect of this initial leaching, four unexposed samples were placed in drinking water for 2 days, weighed and ashed in a muffle furnace to determine initial leaf AFDM. The leaf mass remaining in the bags exposed in streams was expressed as a ratio between samples and initial leached litter expressed in AFDM.

Fungal biomasssults

The role of fungal biomass was demonstrated to highly influence the microbial conditioning of leaves, an important driver of leaf decay into the streams [START_REF] Gessner | Importance of stream microfungi in controlling breakdown rates of leaf litter[END_REF]. Fungal biomass was measured after 34 days of exposure (2nd sampling date) corresponding to the maximum values of fungal biomass expected for this type of leaves in low-order streams [START_REF] Gessner | Importance of stream microfungi in controlling breakdown rates of leaf litter[END_REF]. A set of five 12-mmdiameter discs were cut from five random leaves of each fine mesh bag cleaned sample, avoiding the central veins, and frozen at 18 °C until processing for ergosterol content as a measure of fungal biomass [START_REF] Gessner | Ergosterol-to-biomass conversion factors for aquatic hyphomycetes[END_REF].

Frozen discs were freeze-dried and weighed before estimation of ergosterol content. Ergosterol extraction and quantification were performed following a method based on solid-phase extraction (Waters, Oasis HLB, 60 mg, 3 cm3, Milford, MA, USA) and high-performance liquid chromatography [START_REF] Gessner | Ergosterol-to-biomass conversion factors for aquatic hyphomycetes[END_REF]. Fungal biomass in leaves was expressed as mycelium mass per gram of dry leaf litter.

Benthic invertebrates

Invertebrates from the 34-day exposure coarse mesh bags retained on 0.5-mm sieves (in the field and in the laboratory after leaf wash) were preserved in 70% ethanol until being identified, counted and measured [START_REF] Tachet | Invertébrés d'eau douce[END_REF]). Identification and count were performed under a stereomicroscope (Nikon SMZ-800N) to genus or species when possible or to family, sub-family or tribe for some Diptera. Body length was measured to the nearest millimetres (from the first 50 individuals of each taxon per leaf bag). Biomass of taxa was evaluated by length-mass relationship [START_REF] Benke | Length-mass relationships for freshwater macroinvertebrates in North America with particular reference to the southeastern United States[END_REF][START_REF] Salins Méthot | Macroinvertebrate size-mass relationships: how specific should they be? Freshwater[END_REF]. A number of biotic metrics were calculated for each leaf bag (abundance, richness, Shannon's and Simpson's diversity (indices), Pielou's evenness (index) of benthic macroinvertebrates, richness and densities of Ephemeroptera, Plecoptera and Trichoptera (EPT) and Crustaceans). Organic pollution (OP) tolerance of taxa was also evaluated using saprobic value trait modalities from [START_REF] Tachet | Invertébrés d'eau douce[END_REF] and was expressed as the relative abundance of intolerant or tolerant taxa (i.e. xeno-saprobic and oligo-saprobic taxa or alpha meso-saprobic and polysaprobic taxa, respectively) in a leaf bag. Benthic macroinvertebrates were also assigned into functional feeding groups (FFGs: shredders, collector-gatherers, collector-filterers, scrapers, predators, parasite-piercers) according to their affinity score (i.e. 0-5) fixed by feeding habits as described by [START_REF] Tachet | Invertébrés d'eau douce[END_REF]. Since invertebrate biomass is known as one of the most important drivers of leaf litter breakdown in headwater streams, benthic macroinvertebrate communities assigned to each FFG (e.g. shredder) were expressed as relative biomass of FFG and as total biomass of FFG per gram of leaf litter AFDM remaining, in a leaf bag. It should be noted that ben-thic macroinvertebrate communities involved in the leaf decomposition process were only a part of the whole communities leaving in the river bed, and so, interpretations could just underline differences of this part of the whole communities.

Data analysis

Leaf litter breakdown rates were determined by fitting massloss data into negative exponential decay models (%AFDMt = ekt, where %AFDMt is the percentage of leaf litter AFDM remaining at timet (day) corrected by the initial AFDM (after leaching) and k (g day1) is the breakdown rate constant) as described by Benfield (2006). The best-fit model was identified by comparing the Akaike information criteria (AIC) [START_REF] Akaike | Information theory and an extension of the maximum likelihood principle[END_REF]) of five leaf litter breakdown models: (i) one parameter: common model with all the data; (ü) two parameters: location of the sites (upstream vs. downstream); (üi) two parameters: dominant land use in the catchment (agricultural vs. forested); (iv) four parameters: agriculture upstream, agriculture downstream, forest upstream and forest downstream or (v) ten parameters: the ten sampling points. According to this method, the best model is defined by the lowest AIC. Then, differences between fitted models were compared with nested model testing. Furthermore, in order to have an overview of the data distribution, the best regression model, the common model and coarse mesh bag leaf litter AFDM were plotted and analysed graphically. Comparison of the k values (corresponding to the various breakdown rates) among identified treatments of the best model was done using their 95% confidence intervals.

Fungal biomass and invertebrate metrics were compared by Kruskal-Wallis test or two-way ANOVA followed by Turkey's HSD tests [START_REF] Zar | Biostatistical Analysis[END_REF] to identify structural differences of communities among site locations (upstream vs. downstream from fishponds) and dominant land use (DLU) in the catchments (forested vs. agricultural). The Bray-Curtis dissimilarity matrix was computed [START_REF] Bray | An ordination of the upland forest communities of southern Wisconsin[END_REF] on the log-transformed abundance taxa data and then tested by a permutational multivariate analysis of variance (PERMANOVA) (McArdle and Anderson 2001) using location, DLU and their interaction to evaluate differences in the composition of benthic macroinver-tebrate communities. Non-metric multidimensional scaling (NMDS) ordination was then used to ordinate Bray-Curtis dissimilarities between samples, and results were graphically represented. Ellipses of each significantly different group identified with the PERMANOVA were overlaid on the graph representing the 95% confidence interval standard error from the centroid of each group.

The relative FFG biomass (expressed as percentage of total macroinvertebrate biomass) was used to investigate the macroinvertebrate distribution among functional groups. Besides, the mean FFG biomass per gram of leaf litter (g1 AFDM) was determined in an attempt to reflect the macroinvertebrate feeding activity potential from the communities.

Overall differences among groups (Location, DLU, Loc*DLU) were analysed with PERMANOVA (Bray-Curtis dissimilarity matrix on the log-transformed FFG data).

All data analyses were performed with R software (R development Core Team 2008) using 'vegan', 'ggplot2', 'nlstools' and 'lattice' packages. The significance level for all statistical analyses was set at 0.05.

3.

Results ______________________________________________________________________

Water physical and chemical parameters

The physical and chemical parameters (Table 2) showed that water was of very good quality according to the French reference document SEQ-Eau (MEDD and French water agency 2003).

Mean daily temperature during the experimental period was <8 °C. Temporal trends in the temperature raw data were identified by applying a locally weighted scatterpoint smoothing (lowess) (Figure in SI). Temperature slowly increased during the study period. There was no difference in temperature among sites, but differences were observed between upstream and downstream locations. Significant differences among physico-chemical parameters between upstream and downstream locations of each pond were detected by means of non-parametric paired Wilcoxon signed-rank test. In the forested sites, temperature appeared to be significantly lower upstream compared to downstream. In the agricultural sites, temperature was higher upstream compared to downstream at the beginning of the experiment. This tendency was inversed in the middle of the experiment. There was also a difference in temperature variation between upstream and downstream locations. In the forested sites, temperature variations was higher upstream (coefficient of variation = 41.5%) than downstream (CV = 30.5%). In the agricultural streams, temperature variation was lower upstream (CV = 26.7%) compared to downstream (CV = 41.8%).

Dissolved oxygen concentrations were generally above 10 mg/L. No temporal trend was observed during the study period. There was no difference in dissolved oxygen among sites, but differences were observed between upstream and downstream locations for A2 only. In this agricultural pond, lower dissolved oxygen concentrations were measured upstream compared to downstream. pH varied between 6.1 and 8.3 among sites. Slightly lower pH values were measured at the forest sites compared to agricultural sites. Significantly lower pH values were measured at the upstream location compared to the downstream location for F2 and A2. Turbidity measurements were in the range of 4 and 244 NTU. Turbidity was significantly higher upstream compared to downstream for A1, A2 and F2.

Conductivities were in the range of 83 and 1202 μS/cm, a subject to large temporal and spatial variation. Similarly, nitrate concentrations were in the range of 0.5 and 35.1 mg/L. Higher conductivities were associated with agricultural sites, presumably as a result of fertilizer application (beginning of March) or desorption from soil. In the agricultural sites, nitrate concentrations were significantly higher upstream compared to downstream.

Leaf litter breakdown

Comparison of the five AIC (corresponding to the five leaf litter breakdown models) highlighted two best-fit models (Table in SI). These models considered either four parameters (model (iv): upstream forested sites: Up F, downstream forested sites: Down F, upstream agricultural sites: Up A and downstream agricultural sites: Down A) or ten parameters (model (v): five streams × two locations) (AIC = 192.78 and 192.77,respectively). The two models did not show any significant statistical differences (ANOVA, p = 0.08, Table in SI). We therefore chose to retain the more parsimonious one, i.e. the model requiring the smallest number of parameters (four parameters, namely Up F, Down F, Up A and Down A). The values of remaining leaf litter were well distributed around the common model values (model (i)) for both upstream sites, while they were always above the mean model values for downstream agricultural sites and below the mean model values at the last date for downstream forested sites (Fig. 2). The four leaf litter breakdown rates (g day1) obtained were 0.027 (Up F), 0.021 (Up A), 0.032 (Down F) and 0.011 (Down A) (Fig. 3). The breakdown rate 95% interval confidences analysed revealed significant differences between downstream forested sites and agricultural sites (greater decomposition at forested sites). Significant differences between downstream agricultural sites and all other sites were measured (with two to three times lower breakdown rates for downstream agricultural sites).

Decomposition rates were also expressed in terms of grams (degree day)1 in place of per day so as to correct differences among sites. The four leaf litter breakdown rates (grams (degree day)1) obtained were 0.006 (Up F), 0.004 (Up A), 0.007 (Down F) and 0.002 (Down A). The breakdown rate 95% interval confidences analysed revealed significant differences between forested and agricultural sites (greater decomposition at forested sites) as well as between downstream and upstream agricultural sites.

Fungal biomass

Fungal biomass results are shown in Fig. 4. After 34 days of exposure in the streams, fungal biomass on maple leaves varied between 15.4 (Down A) and 79.3 (Down F) milligrams of mycelium gram of leaf dry mass among the studied groups. Fungal biomass was markedly affected by location (upstream vs. downstream) and dominant land use in the catchments (agricultural vs. forested) (two-way ANOVA, p = 0.03, p < 0.001 and interaction p < 0.001, respectively). Fungal biomass was similar at upstream forested and agricultural sites. Much lower values of fungal biomass were observed only at downstream agricultural sites. No significant difference was measured in forested sites between upstream (Up F) and downstream (Down F) locations (Fig. 4).

Macroinvertebrate communities

A total of 6792 benthic macroinvertebrates were sampled in coarse mesh bags after 34 days of experiment, resulting in 44 taxa distributed in 32 families. The structure of the communities was not significantly different among the four studied groups (Up F, Up A, Down F, Down A) for individual abundances, taxon richness, Shannon and Simpson diversities and Pielou evenness (p = 0.13, p = 0.86, p = 0.46, p = 0.47, p = 0.39, respectively) (Table 3). Significant differences in relative EPT abundances were found between both dominant land uses in the catchments (forest vs. agriculture) with around 30-fold more EPT in the forest communities than in the agricultural communities (p < 0.001), but no significant difference was observed between upstream and downstream locations. Significant differences in EPT richness were found between three groups (Up F, Down F, Up/Down A) with greater values at the upstream forested sites (p < 0.001) (Table 3). By contrast, significant differences in crustacean abundances were found between the studied groups with 15-fold lower relative abundances at the upstream forested sites than at the upstream agricultural sites (p < 0.001) but no differences were observed for the downstream sites (intermediate values). Furthermore, significant differences in the percentages of intolerant taxa (xeno-and oligo-saprobic taxa) and tolerant taxa (alpha meso-and poly-saprobic taxa) with regard to organic matter pollution only discriminated the upstream forested sites characterized by the highest relative abundance of intolerant taxa and the lowest relative abundance of tolerant taxa (p < 0.001, p = 0.003).

Results of the PERMANOVA on pairwise Bray-Curtis dissimilarity matrices of benthic macroinvertebrate abundances (named taxa composition) showed significant effects of land use (p = 0.001), location (p = 0.001) and their interaction (p = 0.002) (Table 4). The ordination of the studied groups' communities by NMDS showed four different community compositions from the overlap of confidence interval ellipses (Fig. 5).

The main feeding habits (FFG) of the benthic macroinvertebrate communities were not significantly different between the upstream and downstream sites (from barrage fishponds) as underlined by the results of PERMANOVA analyses on the composition and total biomass of functional feeding groups (Table 4). Moreover, these results highlighted functional differences according to the dominant land uses in the catchments. Our findings showed that shredders dominated the functional composition of all the communities with at least 50% of the biomass represented by shredder taxa (Fig. 6a). Furthermore, gatherers, filterers and parasites always represented less than 8% of the biomass of the different communities whatever the studied group (Fig. 6a). The percentage of shredder biomass decreased from forested sites (66-65%) to agricultural sites (58-50%) and scraper biomass displayed the inverse trend with values ranging from 11 to 13% in forested sites and from 22 to 18% in agricultural sites. The relative biomass of predators varied between 8 and 14% among the groups (Up F, Up A, Down F, Down A). Concerning the total biomass of communities spread in the FFGs, the results highlighted higher biomass for forested sites (Fig. 6b). This greater biomass in forested sites was well divided among the FFG with a (quasi) systematic higher biomass obtained in forested sites. The greatest difference was observed in the total biomass of shredders. The latter was almost 10-fold higher for forested sites than it was for agricultural sites. Similarly, the biomass of shredders showed an average value 2-fold higher for downstream than for upstream forested sites (Fig. 6b). Nevertheless, at downstream site F3, the shredder biomass was almost eight times lower than at both other downstream sites (F1 and F2) and subsequently associated with EPT density decrease (data not shown).

4.

Discussion ___________________________________________________________________

In order to evaluate the effects of barrage fishponds on headwater streams, we will first discuss the need to take into account the land uses within the studied catchments. On these grounds, we will then discuss the results we obtained for invertebrate communities, fungal biomass and leaf litter decomposition for a better understanding of fishpond effects on headwater stream functioning.

Effect of land use on headwater stream functioning

Looking only at the upstream sampling points of the five studied sites, we did not record any differences in benthic macroinvertebrate abundance, taxon richness, diversity or evenness between agricultural and forested sites. Similarly, shredders were dominant for both land uses suggesting that shredding activity is an important ecosystem process in running waters, especially in headwater streams where allochtho-nous leaf litter is abundant. However, we recorded significant differences in benthic macroinvertebrate community compositions depending on the catchment land use (forest vs. agriculture). We noted only 45% of common species, as well as 30-fold more EPT, 15-fold less crustaceans and almost 14% more taxa intolerant to organic pollution at upstream sites located in forested catchments compared to agricultural catchments. The quasi-systematic disappearance of the most sensitive taxa in agricultural catchments (Capnia sp., Nemoura sp., HabØphlebia sp., Oligostomis reticulata) indicated that land use considerably influences the stream community. Several studies have already shown the effects of agriculture on terrestrial [START_REF] Mazzia | The impact of management strategies in apple orchards on the structural and functional diversity of epigeal spiders[END_REF][START_REF] Frampton | The effects on terrestrial invertebrates of reducing pesticide inputs in arable crop edges: a meta-analysis[END_REF] and aquatic [START_REF] Castela | Evaluation of stream ecological integrity using litter decomposition and benthic invertebrates[END_REF][START_REF] Schäfer | Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe[END_REF]) macroinver-tebrate communities. These communities are often affected by pesticide inputs especially for the most sensitive taxa (e.g. [START_REF] Ippolito | Sensitivity assessment of freshwater macroinvertebrates to pesticides using biological traits[END_REF][START_REF] Lecerf | Perturbations anthropiques et fonctionnement écologique des cours d'eau de tête de bassin :étude du processus de décomposition des litières[END_REF][START_REF] Schäfer | Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe[END_REF]). During our study, Gaillard et al. (2016a, b) followed pesticides in three out of our five studied sites (i.e. F2, A1 and A2). They measured very low concentrations of pesticides upstream from F2 (a forested site, with maximum concentrations of tritosulfuron, 0.1 μg.L1). Forested sites F1 and F3 were not investigated in this study for pesticides. However, since F2 was by far the most exposed forested site to agricultural inputs (among the ponds under forest cover, F2 had the catchment with the highest arable land percentage and high nitrate concentrations confirmed agricultural inputs; Tables 1 and2), we could expect that pesticide levels in F1 and F3 would be lower than in F2. On the other hand, agricultural sites were characterized by high pesticide concentrations (with maximum MCPA concentrations of 26.5 μg L1 measured upstream from A1 and 8.26 μg L1 measured upstream from A2). As such, we cannot exclude that pesticides associated with conventional crop production could have an effect on the macroinvertebrate communities.

Beyond EPT, the absence of a relationship between macroinvertebrate taxon richness, diversity, evenness and land use was in agreement with other studies whose authors did not establish any correlation between pesticides (associated with agricultural land use) and those metrics [START_REF] Brock | On the choice of structural parameters and endpoints to indicate responses of freshwater ecosystems to pesticide stress[END_REF][START_REF] Maltby | Spray drift of pesticides and stream macroin-vertebrates: experimental evidence of impacts and effectiveness of mitigation measures[END_REF]. In addition, [START_REF] Morrissey | Neonicotinoid contamination of global surface waters and associated risk to aquatic invertebrates: a review[END_REF] showed that crustacean taxa are often less sensitive than EPT taxa to acute or chronic insecticide exposition, which could corroborate our results concerning the high crustacean abundance in streams from agricultural catchments. Thus, the aquatic macroinvertebrates of agricultural sites are probably less sensitive taxa to diverse agricultural perturbations, as also suggested by the results of relative abundance of intolerant organisms.

In the literature, community modifications related to anthropogenic environmental perturbations, which are themselves linked to conventional agricultural practices, were also often associated with changes in their functional feeding habits [START_REF] Flores | Effects of a fungicide (imazalil) and an insecticide (diazinon) on stream fungi and invertebrates associated with litter breakdown[END_REF][START_REF] Rasmussen | Impacts of pesticides and natural stressors on leaf litter decomposition in agricultural streams[END_REF][START_REF] Schäfer | Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe[END_REF]). [START_REF] Rasmussen | Impacts of pesticides and natural stressors on leaf litter decomposition in agricultural streams[END_REF] reported a decrease in leaf litter breakdown with increasing pesticide inputs. However, in our study, the results were not as straightforward. Indeed, the leaf litter breakdown rate did not vary significantly between forested and agricultural upstream sites when it was expressed in grams per day and was only lower when corrected in grams (degree day)1. Moreover, the higher temperatures observed in agricultural sites could have been due to riparian management (i.e. streams with little riparian shading in agricultural sites). Besides, shredder biomass showed significant differences between forested and agricultural sites suggesting potential feeding habit differences. We observed lower shredder biomass in the agricultural catchments, but most of it came from Gammaridae (Gammarus pulex and G. roeseli) compared to the forested sites where EPT dominated. Several studies showed that high leaf litter breakdown rates can be correlated with shredder communities dominated by Gammaridae [START_REF] Dangles | Species richness-decomposition relationships depend on species dominance: biodiversity and dominance in ecosystems[END_REF][START_REF] Piscart | Effects of intense agricultural practices on heterotrophic processes in streams[END_REF]. As already outlined, Gammaridae appeared to be better shredders than some EPT taxa. Such information is noteworthy given the fact that leaf litter fungal conditioning was similar for the various upstream sites whatever the land use. Despite the known effects of agriculture on microbial activity measured in several studies [START_REF] Bundschuh | Fungal composition on leaves explains pollutant-mediated indirect effects on amphipod feeding[END_REF][START_REF] Robinson | Nutrient addition accelerates leaf breakdown in an alpine Springbrook[END_REF][START_REF] Rasmussen | Impacts of pesticides and natural stressors on leaf litter decomposition in agricultural streams[END_REF][START_REF] Schäfer | Effects of pesticides monitored with three sampling methods in 24 sites on macroinvertebrates and microorganisms[END_REF], our results did not indicate any direct effects of agriculture on fungal biomass. However, nutrient losses in freshwater associated with soil fertilization by farmers were investigated in several laboratory and field studies showing that moderate increase in nutrients favoured microbial biomass and activity [START_REF] Colas | Synergistic impacts of sediment contamination and dam presence on river functioning[END_REF][START_REF] Ferreira | A metaanalysis of the effects of nutrient enrichment on litter decomposition in streams: Nutrient enrichment and litter decomposition[END_REF][START_REF] Robinson | Nutrient addition accelerates leaf breakdown in an alpine Springbrook[END_REF][START_REF] Suberkropp | Effect of dissolved nutrients on two aquatic hypho-mycetes growing on leaf litter[END_REF]). This is consistent with the nitrate concentrations measured at our sampling points which were correlated with land uses, with 2-to 16-fold higher values at upstream agricultural sites (Fig. 1 and Table 2).

Effect of barrage fishponds on the ecosystem functioning of headwater streams

Numerous studies have investigated the effects that daminduced water flow changes have on the fish and macroinvertebrate communities of downstream ecosystems [START_REF] Bunn | Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity[END_REF][START_REF] Bredenhand | Impact of a dam on benthic macro-invertebrates in asmall river in a biodiversity hotspot: Cape Floristic Region, South Africa[END_REF][START_REF] Casas | The effect of a headwater dam on the use of leaf bags by invertebrate communities[END_REF][START_REF] Martínez | Stream regulation by small dams affects benthic macroinver-tebrate communities: from structural changes to functional implications[END_REF]). Nevertheless, only few of them conjointly considered community modifications and leaf litter breakdown, a fundamental ecosystem function of headwater streams [START_REF] Casas | The effect of a headwater dam on the use of leaf bags by invertebrate communities[END_REF][START_REF] Martínez | Stream regulation by small dams affects benthic macroinver-tebrate communities: from structural changes to functional implications[END_REF]). In any case, to our knowledge, no research had been conducted yet on leaf litter breakdown and their associated communities when extensive fish production in pond dams was coupled with land uses in the river basins. Land uses can affect habitat quality, which is an important parameter for communities and their functional activities in streams. On these grounds, we hypothesized that barrage fishponds could have differential effects on the functioning of headwater streams depending on the catchment land use.

For example, as in forested catchments, ponds induce large aperture in the canopy and consequently change the access to solar energy, in this study, we observed an increase in downstream temperatures only in the forested catchments. This observation is consistent with another study whose authors showed an increase in water temperatures (1 °C, on average) between upstream and downstream sites (Tou c hart and [START_REF] Touchart | La gestion du risque thermique en étang[END_REF]. Concerning water quality (monitored during the step of growing fish, see BStudy sites^ section), our results indicated a decrease in nutrients and turbidity (used as an indicator of suspended matter concentration) from upstream to downstream sites (Table 2). The differences were even more significant when concentrations were at their highest at upstream sites (e.g. in agricultural catchments and F2). This is in agreement with observations made by several authors who reported suspended matter and nutrient retention in ponds [START_REF] Banas | Flux de matière en étangs piscicoles extensifs :Rétention, Sédimentation, Exportation[END_REF][START_REF] Passy | Restoration of ponds in rural landscapes: modelling the effect on nitrate contamination of surface water (the Seine River Basin, France)[END_REF]. For example, [START_REF] Passy | Restoration of ponds in rural landscapes: modelling the effect on nitrate contamination of surface water (the Seine River Basin, France)[END_REF] measured 27 to 56% nitrogen retention in ponds depending on the studied year. Gaillard et al. (2016a, b) observed the same trend for pesticide concentrations suggesting that fishponds can improve downstream water quality, an important parameter for macroinvertebrate and fungal communities.

The results of our investigation into the structure and composition of the benthic macroinvertebrate communities showed that upstream/downstream differences were more pronounced for fishponds in forested catchments than for those in agricultural catchments. For example, EPT richness and abundance exhibited greater decreases at forested sites (between Up F and Down F) than at agricultural sites (differences of 15.3 and 14.5%, respectively, for forested sites vs. 0.4 and 7.5% for agricultural sites; Table 3). The results on the organic tolerance status of our communities showed the same trend. The percentages of tolerant and intolerant individuals observed at downstream forested sites were close to those recorded at upstream and downstream agricultural sites. These results were in accordance with a decrease in or disappearance of the most sensitive taxa (Capnia sp., Habrophlebia sp. and O. reticulata) between upstream and downstream forested sites, whereas they were already absent in agricultural sites. Indeed, we observed more common species at the downstream sites (i.e. at both agricultural and forested downstream sites) than at the upstream sites (57.5 and 31%, respectively). A hypothesis is that alteration to the flow regimes of waters downstream from the fishponds could favour the homogenization of the benthic habitat, which in turn could have an effect on the whole benthic macroinvertebrate community structure and composition [START_REF] Tachet | Invertébrés d'eau douce[END_REF]). However, no substratum homogenization was visually observed at our downstream sites as compared with upstream sites (personal observation). This phenomenon might be related to substratum composition because fishponds in the Lorraine region are already located on homogenous and impermeable clay substratum [START_REF] Cral ; | Benthic algae stimulate leaf litter decomposition in detritus-based headwater streams: a case of aquatic priming effect?[END_REF]. By contrast, the higher macroinverte-brate community similarity and the higher densities of organic intolerant taxa observed for upstream and downstream agricultural sites are consistent with the results of other authors who reported non-cumulative effects of multiple stressors on co-tolerant species [START_REF] Christensen | Multiple anthropogenic stressors cause ecological surprises in boreal lakes[END_REF][START_REF] Kneitel | Trade-offs in community ecology: linking spatial scales and species coexistence[END_REF][START_REF] Vinebrooke | Impacts of multiple stressors on biodiversity and ecosystem functioning: the role of species co-tolerance[END_REF].

As an example, [START_REF] Kneitel | Trade-offs in community ecology: linking spatial scales and species coexistence[END_REF] suggested that ecological trade-offs can have synergetic interactions on the taxa since the exposure to one stressor can select for species or individuals tolerant to that stressor (e.g. in this study, pesticides or nutrients in agricultural catchments) but potentially to an additional stressor (e.g. in this study, fishpond effects due to increasing temperatures or autochthonous organic matter in the downstream rivers). Furthermore, as evidenced by Gaillard et al. (2016a, b), fishponds can have positive effects on pesticide concentrations in streams. Indeed, fishponds cause disruptions of stream continuum (i.e. of fine particulate organic matter (FPOM)/coarse particulate organic matter (CPOM) ratio, temperature and ecological continuity), which could have antagonist effects on the stream water quality that could diminish their impact on mac-roinvertebrate communities especially when the stream ecosystems are already disrupted by toxicants and nutrient runoff. [START_REF] Ward | The serial discontinuity concept of lotic ecosystems[END_REF] showed that headwater dams can increase the FPOM/CPOM ratio due to preferential sedimentation of the largest and heaviest particulates in lentic systems. Indeed, while determining the FFGs of these communities, we expected a shift in the main feeding habits of the downstream communities with more collector feeders and fewer shredders. Surprisingly, according to the PERMANOVA results, we did not measure any significant differences in the functional feeding composition and biomass of the communities between upstream and downstream sites. Shredders dominated the benthic macroinvertebrate communities regardless of the group studied (Up A, Up F, Down A or Down F). This finding suggests that coarse particulate organic matter remains the main food source from all our headwater stream ecosystems (despite land use modification and fishpond addition) and that community modifications have no direct effects on FFGs.

From our results on the influence of fishponds on the main feeding habits of macroinvertebrate communities estimated by calculating the percentage of each FFG in biomass (Fig. 6b), little effect was expected on the leaf litter breakdown rates downstream from the fishponds. However, we measured that leaf litter breakdown remained stable at downstream forested sites (in g day1 or in grams (degree day)1), but dropped at downstream agricultural sites (in g day1 or in grams (degree day)1) despite stability of the shredder biomass. Applying fungal biomass as an indicator of leaf litter microbial conditioning [START_REF] Gessner | Growth and production of aquatic hy-phomycetes in decomposing leaf litter[END_REF], our results have evidenced significant differences for downstream agricultural sites with three to four times less fungal biomass than in all other studied groups. Numerous studies reported close correlation between leaf conditioning by hyphomycetes and an increase in palatability or nutritional value of leaf litter (Gessner andChauvet 1997, 2002;[START_REF] Gulis | Interactions between stream fungi and bacteria associated with decomposing leaf litter at different levels of nutrient availability[END_REF]. Furthermore, authors have shown that shredders preferentially choose and feed on conditioned leaf litter with high palatability and high nutritional value [START_REF] Chung | Contribution of fungal biomass to the growth of the shredder, Pycnopsyche gentilis (Trichoptera: Limnephilidae)[END_REF][START_REF] Graça | The role of invertebrates on leaf litter decomposition in streamsareview[END_REF][START_REF] Graça | A conceptual model of litter breakdown in low order streams: litter breakdown in low order streams[END_REF][START_REF] Nelson | Gammarus-microbial interactions: a review[END_REF]. Thus, even if we have no significant difference in the biomass (relative and total) of shredders between upstream and downstream sites for both catchment types, the difference in shredding activity could be explained by lower palatability of leaf litter at downstream agricultural sites.

The decrease in fungal biomass observed at downstream agricultural sites suggests that fishponds affect relatively more the fungal community growth in agricultural sites than in forested sites. As highlighted in the first part of this discussion, we hypothesized that fungal community growth was sustained by high nutrient content in upstream agricultural sites. However, at downstream agricultural sites, values were 6-to 8-fold lower (for nitrates) than at upstream sites. These results could provide some explanation for such low fungal biomass at downstream agricultural sites. Moreover, based on the passive and short-distance dispersal of hyphomycetes [START_REF] Bärlocher | Reproduction and dispersal in aquatic hyphomycetes[END_REF][START_REF] Thomas | Changes in concentration of aquatic hyphomycete spores in Lees Creek, ACT, Australia[END_REF] and the low availability of growing substrate (leaf litter) especially in fishponds from agricultural sites (lack of riparian trees when compared to forested sites; Fig. 1), we therefore suggest that fishponds in agricultural landscapes could cause larger gaps for hyphomycete ecological continuity than in forested landscapes.

5.

Conclusion __________________________________________________________________

The barrage fishponds under investigation showed that they affect ecosystem functioning in headwater streams, but the effects depend on the dominant land use in the catchment. In forested landscapes, we observed that fishponds have significant effects on the structure and composition of benthic mac-roinvertebrates and moderate effects on ecosystem functioning. In agricultural landscapes, we found that the structure and composition of benthic macroinvertebrates are less affected but, on the other hand, that the ecosystem functioning is strongly affected by fishponds. Most studies dealing with the influence of water flow and quality alteration associated with dams or barrage fishponds are based essentially on physical and chemical water parameters or on community structure (mostly macroinvertebrates and fish).

Our study indicates that those parameters are not sufficient to estimate freshwater status especially when it concerns the effects of complex stressors such as dams or fishpond dams in this study. 329) 402 ( 425) 333 ( 207) 679 ( 444) 512 ( 320) 997 ( 389) 291 ( 221) 310 ( 351) 3060 ( 874) 3054 ( 698)

Temperature (°C) b 5.1 (2.

Figure 1 .Figure 2 .Figure 3 .Figure 4 .Figure 6 .

 12346 Figure 1. Location of the study area (a France; b Lorraine region) andland use of the river basins withstudy site locations

Table 1 . Fishpond, catchment and riparian characteristics of the ten sampling points under study in the Lorraine region

 1 Heavy Alnus, Heavy Alnus Heavy Alnus Heavy Alnus Heavy Corylus

	Parameter	F1		F2		F3		A1		A2
		Up	Down	Up	Down	Up	Down	Up	Down	Up	Down
	Fishpond
	characteristics Pond area (ha)	-	4.9	-	1.6	-	7.2	-	4.4	-	31.6
	Pond volumes (m 3 )	-	49,000	-	16,000	-	108,000	-	44,000	-	568,000
	Hydraulic retention	-	1:20	-	1:42	-	1:21	-	1:20	-	1:11
	time Catchment
	characteristics Catchment area (ha) 74.1	100.6	46.1	64.4	85.8	152.8	34.4	86.2	209.3	345.3
	Forest (%)	58.6	64.6	61.0	62.8	98.2	94.3	-	-	0.6	1.6
	Arable land (%)	4.2	3.1	39.0	25.4	1.8	1.0	49.9	39.3	44.7	38.6
	Permanent pastures	37.2	27.4	-	9.3	-	-	34.4	49.3	51.9	43.6
	(%) Oters (e.g. roads,	-	4.9	-	2.5	-	5.7	15.7	11.4	2.8	16.2
	buildings + pond)
	(%)
	Sampling point
	characteristics
	Land use near to the Deciduous Forest Deciduous Deciduous Deciduous Deciduous Deciduous forest Agricultural Agricultural land Agricultural land Agricultural land
	sampling point		forest	forest	forest	forest		land
	Tree cover	Heavy						Isolated	Partial	Weak	Partial
	Dominant genders in	Quercus, Fraxinus,	Fraxinus,	Carpinus,	Fraxinus,	Carpinus,	Quercus,	trees Carex, Typha,	Phragmites	Alnus, Urtica Urtica,	Phragmites,
	riparian vegetation	Fagus	Carpinus	Quercus	Carpinu	Quercus	Fagus	Phragmites		Prunus,	Prunus, Urtica
					s					Carex

Table 2 . Physical and chemical characteristics (mean (SD)) of the five study sites or the ten sampling points in the Lorraine region

 2

	Parameter	F1		F2		F3		A1		A2
		Up	Down	Up	Down	Up	Down	Up	Down	Up	Down
	Stream order e	Intermittent	Intermittent	Intermittent	Intermittent	Intermittent	1	1	1	1	1
	Mean outflow rate (m 3 j 1 ) a	416 (340)		506 (386)		754 (415)		301 (292)		3057 (787)
		430 (

  n = 6 for F1 and F3, n = 3 for each sampling point; 15 min recording for F2, A1, A2 b 15 min recording during all the studying periods c n = 8-12 for each study site and 4-6 for each sampling point d n = 8 for each study sites and 4 at each sampling point e As defined by[START_REF] Strahler | Quantitative analysis of watershed geomorphology[END_REF]

		1)		5.2 (1.8)		5.4 (1.9)		5.1 (2.0)		5.5 (1.7)
		4.5 (2.4)	5.5 (1.5)	4.8 (1.7)	6.0 (1.9)	4.8 (1.7)	5.6 (1.8)	5.2 (1.1)	4.9 (2.2)	5.4 (1.7)	5.6 (2.2)
	Conductivity (μS cm 1 ) c	212.0 (49.6)		545.2 (133.8)		166.8 (69.4)		824.9 (161.9)		938.9 (84.5)
		245.5 (64.2)	178.5 (1.3)	609.9 (169.8)	480.6 (23.8)	104.7 (26.1)	228.9 (16.7)	899 (177.1)	750.8 (114.5)	1003.9 (74.7)	873.9 (4)
	Turbidity (NTU) c	59.1 (27.1)		61.1 (51.9)		32.6 (15.8)		47.2 (35.4)		16.2 (9.9)
		68.9 (36.8)	49.2 (10.3)	88.1 (63.1)	34.1 (9.3)	46.9 (5.2)	18.2 (3.1)	63.1 (42.4)	31.3 (18.8)	22.9 (9.7)	9.5 (4)
	O2 (mg/L) c	11.3 (1.1)		10.7 (1.7)		10.3 (1.0)		11 (1.1)		11.8 (1.1)
		11.3 (1.4)	11.4 (1.0)	9.8 (1.9)	11.7 (0.9)	10.1 (0.7)	10.5 (1.2)	10.6 (1.2)	11.4 (0.9)	11.1 (1.2)	12.5 (0.8)
	pH c	7.6 (0.2)		7.6 (0.4)		6.4 (1.3)		7.8 (0.5)		8.1 (0.4)
		7.6 (0.3)	7.7 (0.1)	7.5 (0.3)	7.8 (0.4)	6.7 (0.4)	6.1 (1.9)	7.7 (0.6)	7.8 (0.5)	7.8 (0.3)	8.3 (0.3)
	Nitrate (mg L 1 ) d	1.8 (0.5)		5.9 (4.4)		0.8 (0.4)		19.6 (15.2)		12.8 (10.1)
		1.93 (0.84)	1.73 (0.26)	10.57 (3.35)	2.43 (1.05)	0.58 (0.08)	0.83 (0.30)	33.94 (0.98)	4.79 (0.81)	21.73 (1.08)	3.02 (0.69)
	Nitrite (mg L 1 ) d	0.01 (0.00)		0.01 (0.00)		0.01 (0.01)		0.04 (0.04)		0.02 (0.01)
		0.01 (0.00)	0.01 (0.00)	0.01 (0.00)	0.01 (0.00)	0.01 (0.00)	0.01 (0.00)	0.05 (0.02)	0.02 (0.00)	0.02 (0.01)	0.01 (0.00)
	Ammonia (mg L 1 ) d	0.09 (0.06)		0.02 (0.01)		0.04 (0.03)		0.33 (0.37)		0.05 (0.02)
		0.11 (0.07)	0.11 (0.07)	0.03 (0.00)	0.02 (0.01)	0.07 (0.02)	0.02 (0.01)	0.71 (0.26)	0.02 (0.02)	0.05 (0.02)	0.06 (0.01)
	SRP (mg L 1 ) d	0.04 (0.00)		0.03 (0.02)		0.04 (0.02)		0.04 (0.01)		0.03 (0.01)
		0.04 (0.00)	0.03 (0.00)	0.05 (0.01)	0.03 (0.01)	0.05 (0.01)	0.04 (0.02)	0.06 (0.00)	0.03 (0.01)	0.04 (0.01)	0.03 (0.01)
	SRP soluble reactive phosphorus

a

Table 3 . Effect of fishponds on streams for selected benthic invertebrate metrics and considering the dominant land use in the catchment

 3 Results indicate the mean value of each group, and the various letters indicate significantly different groups. % Intolerant taxa =community score of xeno-+ oligo-saprobic taxa according to Tachet et al. 2010 saprobic value classification. % Tolerant taxa =community score of α-meso-+ poly-saprobic taxa according to Tachet et al. 2010 saprobic value classification EPT Ephemeroptera + Plecoptera + Trichoptera *0.01 > p > 0.001; **p < 0.001; Kruskal test or two-way ANOVA followed by Tukey's post hoc test

		Forest		Agriculture
		Up	Down	Up	Down
	Individuals' abundance	260.5 311.7	186.1 210.1
	Taxa richness	8.8	9.7	8.7	8.6
	%EPT( abundance)**	72.9 a 57.6 a	2.2 b	1.8 b
	%EPT (richness)**	36.6 a 22.1 a,b	16.6 b	9.1 b
	% Crustaceans	4.1 a 18.9 a,b	62.2 c	26.4 b,c
	(abundance)** Shannon's diversity. H′	1.65	1.58	1.66	1.35
	Simpson's diversity (1-D)	0.56	0.48	0.54	0.45
	Pielou's equitability. J′	0.54	0.47	0.54	0.44
	% Of intolerant taxa*	49.0a 37.8b	35.4b 34.8b
	% Of tolerant taxa*	18.9a 28.2b	29.8b 31.5b

Table 4 . Results of the PERMANOVA analysis on taxa composition, as well as composition of the functional feeding groups (FFG) in relative biomass and in total biomass per dominant land use (forest vs. agriculture), location (upstream vs. downstream from the fishponds) and their interactions

 4

	Community parameter	Source of variation	df Sum of squares Mean squares Pseudo F p value Significance
	Taxa composition	Location	1	0.7636	0.76357	7.3368	0.001	***
		Land use near to the sampling point 1	1.7625	1.76247	16.9349	0.001	***
		Land use: location	1	0.5088	0.50878	4.8887	0.002	**
		Residuals	36 3.7467	0.10407	0.55248
		Total	39 6.7815	1.00000
	FFG in relative biomass	Location	1	0.03222	0.032222	1.2875	0.27	NS
		Land use	1	0.13376	0.133757	5.3444	0.004	**
		Land use: location	1	0.0272	0.027196	1.0866	0.324	NS
		Residuals	36 0.90099	0.025027	0.82345
		Total	39 1.09416	1.00000
	FFG in total biomass	Location	1	0.00871	0.00871	0.3654	0.753	NS
		Land use	1	0.66539	0.66539	27.9285	0.001	***
		Land use: location	1	0.02827	0.02827	1.1866	0.282	NS
		Residuals	36 0.85769	0.02382	0.54978
		Total	39 1.56005	1.00000
	Results are based on 999 permutations
	NS no significant
	PERMANOVA significance test: *** p < 0.001, **p < 0.05

Acknowledgements ____________________________________________________________

The authors are grateful to the 'Agence de l'Eau Rhin-Meuse' and the 'Zone Atelier Moselle' for their financial support to this project. We sincerely thank P. Hartmeyer, A.