Multi-substance indicators based on caged Gammarus

May 2, 2019 - ... Hydrobiidae densities with increasing levels of metal contamination are identified relative to Baetidae, Chironomidae, and Hydropsyc...
0 downloads 0 Views 501KB Size
Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC

Environmental Measurements Methods

Multi-substance indicators based on caged Gammarus bioaccumulation reveal the influence of chemical contamination on stream macroinvertebrate abundances across France Benjamin Alric, Olivier Geffard, Andre Chandesris, Martial Ferreol, Adeline François, Olivier Perceval, Jeremy Piffady, Bertrand Villeneuve, and Arnaud Chaumot Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b01271 • Publication Date (Web): 02 May 2019 Downloaded from http://pubs.acs.org on May 3, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

Environmental Science & Technology

1

Multi-substance indicators based on caged Gammarus

2

bioaccumulation reveal the influence of chemical contamination on

3

stream macroinvertebrate abundances across France

4

Benjamin Alric†, O. Geffard†, A. Chandesris‡, M. Ferréol‡, A. François†, O. Perceval§, J.

5

Piffady‡, B. Villeneuve‡, A. Chaumot*†

6 7 8

†Irstea,

UR RiverLy, Laboratoire d’écotoxicologie, centre Lyon-Villeurbanne, 5 rue de la

9

Doua CS 20244, F-69625, Villeurbanne, France

10

‡Irstea,

UR RiverLy, Laboratoire d’hydrobiologie quantitative, centre Lyon-Villeurbanne, 5

11

rue de la Doua CS 20244, F-69625, Villeurbanne, France

12

§Agence

13

Vincennes, France

française pour la biodiversité, site de Vincennes, 5 Square Felix Nadar, 94300

14 15

*Corresponding author: Arnaud Chaumot, Irstea, UR RiverLy, Laboratoire

16

d’écotoxicologie, centre de Lyon-Villeurbanne, 5 rue de la Doua CS 20244, F-69625,

17

Villeurbanne, France. [email protected]

18 19 20 21 22 23

Running tittle: Biota for assessing pollutants in the ecosystem

24

ACS Paragon Plus Environment

1

Environmental Science & Technology

25

Abstract

26

Most of anthropogenic stressors affecting freshwater systems are qualitatively known.

27

But the quantitative assessment of contaminant exposure and effects to aquatic

28

communities is still difficult, limiting the understanding of consequences on aquatic

29

ecosystem functioning and the implementation of effective management plans. Here,

30

multi-substance indicators based on caged gammarid bioaccumulated contamination

31

data are proposed (for metals and persistent organic pollutants, POPs) to map the

32

bioavailable contamination level of freshwater ecosystems at a large spatial scale. We

33

assessed the ability of these indicators to highlight the relationships between chemical

34

exposure gradients and alteration in the abundance of macroinvertebrate populations,

35

on a data set of 218 watercourses distributed throughout France. We identified spatial

36

regional heterogeneities in the levels of bioavailable contamination of metals (18

37

compounds) and POPs (43 compounds). Besides, a degradation of Gammaridae,

38

Ephemeridae, and Hydrobiidae densities with increasing levels of metal contamination

39

are identified relative to Baetidae, Chironomidae, and Hydropsychidae. We show here

40

that active biomonitoring allows the establishment of multi-substance indicators of

41

bioavailable contamination, which reliably quantify chemical exposure gradients in

42

freshwater ecosystems. Our ability to identify species-specific responses to chemical

43

exposure gradients demonstrates the promising possibility to further decipher the

44

effects of chemical contamination on macroinvertebrate assemblages through this type

45

of indicators.

Page 2 of 32

46

ACS Paragon Plus Environment

2

Page 3 of 32

Environmental Science & Technology

47

Introduction

48

Human-made chemical contaminants released into aquatic environments have attracted

49

widespread attention during recent years owing to the increasing anthropogenic

50

loading1 as well as their persistence, biological accumulation and potential toxic

51

effects,2,3 which can compromise the future provisioning of vital ecosystem services.4,5

52

Given that water is regarded as the most essential of natural resources,6,7 a forward-

53

looking system of chemical status assessment for surface water bodies has been

54

implemented in Europe and embodied in the Water Framework Directive (WFD8). The

55

measurement of contaminant concentrations in biota is increasingly recommended for

56

assessing the environmental exposure of biological communities to chemicals.9,10

57

Bioaccumulated concentration data collected in biomonitoring surveys provide time-

58

integrated records of the fraction of contaminants potentially toxic for organisms, i.e.,

59

bioavailable and bioaccumulative fraction.11,12 More recently, active biomonitoring

60

(caged organism transplantation) has been proposed as an alternative approach to the

61

field collection of wild organisms, which may improve spatial and time comparisons of

62

contamination levels. Active biomonitoring obviously allows one to avoid the potential

63

lack of native species in studied sites, but also to minimize some confounding factors

64

known to influence the accumulation of contaminants in organisms (e.g., variability in

65

the exposure time, age and size of sampled organisms).13,14 Among freshwater sentinel

66

species, recent interest has been focused on the genus Gammarus because of its

67

widespread distribution throughout the Northern Hemisphere,15,16 and its important

68

ecological role as a trophic link between the base of food webs, as shredder, and higher-

69

order consumers, as food for fish, amphibians and birds.17,18 Also, an efficient

70

accumulation of metal and organic contaminants has been reported in gammarids at

ACS Paragon Plus Environment

3

Environmental Science & Technology

71

environmental exposure levels,19,20 making it a suitable candidate for identifying the

72

sources and origins of contamination in freshwater ecosystems.10,21–23

73

Page 4 of 32

Evaluating the effects of contamination in aquatic ecosystems requires tools that

74

quantitatively link the exposure to chemical contaminants to the impairment of

75

biological communities. Relationships between the contamination level of individual

76

metal measured in native organisms and the abundance of macroinvertebrates known

77

to be sensitive to metals have been described in the literature. This was reported at the

78

scale of single watersheds affected by occasional and local chemical contaminations

79

(e.g., mining discharge),24–26 but also at larger spatial scales (Flanders, Belgium) with a

80

lower ecological quality when levels of bioaccumulated individual pollutants in native

81

eels were elevated.9 Recently, active biomonitoring with caged gammarids was

82

successfully carried out at the French national scale to identify the possible link between

83

exposure to three metal contaminants and native gammarid densities.27 All these

84

findings highlight the effect of bioavailable chemical contaminations on aquatic

85

populations and communities. Nevertheless, in terms of management, the establishment

86

of such relationships obtained separately for each contaminant can be inefficient to

87

quantify the link between contamination gradient and ecological degradation because

88

the exposure to multiple contaminants seems to be the rule rather than the

89

exception.28,29

90

Thus, this study aimed to propose multi-substance indicators, based on the

91

aggregation of concentrations accumulated by caged gammarid biomonitors, to

92

demonstrate their relevance for mapping the contamination level of aquatic systems and

93

to describe responses of macroinvertebrates to chemical pressure. This work benefited

94

from caged gammarid contamination data at the French national scale acquired in the

95

WFD regulatory context by regional water agencies, considering 18 metals and 43

ACS Paragon Plus Environment

4

Page 5 of 32

Environmental Science & Technology

96

persistent organic pollutants (POPs), and records of benthic macroinvertebrate

97

abundance, in 218 French river sites. Thanks to multi-substance indicators developed in

98

this study, we sought in a first step to establish a global snapshot of metal and POP

99

bioavailable contamination of watercourses at the national scale, and in a second step, to

100

evaluate the response of resident macroinvertebrates to bioavailable chemical

101

contamination. According to our previous study,27 we thus assumed that a multi-

102

substance indicator could allow us to describe quantitative relationships between metal

103

or POP environmental contamination gradients and native gammarid abundances.

104

Finally, taking advantage of such an integrated indicator of metal contamination, we

105

inspected the covariation between the estimated metal bioavailable contamination

106

gradient and the distribution of abundances of five additional widely distributed

107

invertebrate families (Ephemeridae, Hydrobiidae, Baetidae, Chironomidae, and

108

Hydropsychidae), known to present contrasting toxicological sensitivities toward metals

109

in ecotoxicological surveys.

110 111

Materials and methods

112

Chemical contamination and ecological bioindication datasets

113

This study was based on two datasets gathering information on (a) the bioavailable

114

chemical contamination in 218 wadeable river sites distributed throughout France, and

115

(b) the abundances of benthic macroinvertebrates in these sites. Chemical and biological

116

sampling dates (one sampling/year) cover the period 2009–2016 and the sampled sites

117

belong to two major surveying networks established by French regional water agencies

118

implementing the WFD requirements: the surveillance monitoring network (SMN) and

119

the operational monitoring network. The former aims to assess long-term changes in the

120

overall surface water status within each catchment and sub-catchment of the river basin

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 32

121

district, while the latter aims to follow the status of French waterbodies identified as

122

being at risk of failing to meet good quality objectives. Hence, available ecological data

123

were obtained by a normalized protocol, applied once a year to sample benthic

124

macroinvertebrates for this purpose.30 This protocol is based on twelve sample units

125

distributed along different predefined habitats in reaches by regarding both their

126

relative coverage but also their hosting capacity (i.e., the capacity of a given substrate to

127

host a rich and diverse invertebrate assemblage). These sample units, performed with a

128

normalized Surber net (sampling area 0.05 m2, mesh size 500 µm), are distributed

129

evenly into three groups according to their habitat representativeness. The first, group

130

A, focused on ‘marginal habitats’ (i.e., with an individual share of less than 5% coverage)

131

with a high hosting capacity, while the other eight sample units, split in two groups (B

132

and C), are taken from ‘major habitats’ (i.e., with an individual share of at least 5%

133

coverage). The last two groups are distinguished by their hosting capacity and the

134

relative coverage of major mesohabitats. From these samples, benthic

135

macroinvertebrate assemblages are characterized by identifying and counting

136

individuals at the genus level.

137

Bioavailable chemical contamination data (contamination levels in caged

138

Gammarus fossarum) were acquired either through French regional water agencies in

139

the WFD context, or through additional research programs carried out by our own

140

laboratory (ANR CESA program GAMMA021 02; ANR Blanc program Multistress0004-

141

03). The methodology used for caged organism transplantation is presented in previous

142

studies.27,31 Briefly, it consists in transplanting during 7 days 20 individuals G. fossarum,

143

which have previously been acclimated 3 weeks in the laboratory under controlled

144

housing parameters. To limit the influence of body size and sex on contaminants

145

accumulation, size-calibrated male individuals were exposed. To avoid the influence of

ACS Paragon Plus Environment

6

Page 7 of 32

Environmental Science & Technology

146

starvation on survival rate and contamination uptake, alder leaves were supplied into

147

cages. Exposed gammarids came from a source population, commonly employed by our

148

laboratory, and characterized by very low concentration levels for both metals and

149

POPs. Hence, a dataset of 313 sampling events at 218 sites was obtained for metals and a

150

dataset of 196 sampling events at 165 sites for POPs. There were 136 common

151

samplings between metal and POP data sets.

152 153

Physical state and land use data

154

The surface water bodies studied in the present work are categorized and characterized

155

according to System B of Annex II of the WFD by considering a climatic, hydrogeological

156

and landform hydro-ecological classification.32 By comparing the distribution of our

157

study sites in the different classes of physical state and land use types with the one

158

determined on the SMN (𝑛 = 1501 sites) we can evaluate how far our set of sites

159

encompass the diversity of rivers over the whole metropolitan French territory. The

160

physical state of rivers was based on the Strahler number and the geological nature of

161

the watershed soils. Digital land cover information, obtained from the CORINE Land

162

Cover 2006 database (CLC2006),33 was used in order to identify land use types in the

163

watersheds. On the basis of river catchment delimitations estimated beforehand with a

164

digital field model, land use cover proportions were calculated at the scale of entire

165

watersheds using ArcInfo analytical tools (ArcGis 10.5 software). We focused on two

166

groups of land use categories likely to induce chemical contamination: urbanization of

167

watersheds (i.e., urban areas, industrial and commercial areas, roads and highways,

168

mines, dumps and construction sites, non-agricultural artificial green spaces), and

169

intensive agriculture (i.e., arable lands, permanent crops, orchards, vineyards, annual

170

crops associated with permanent crops, cropping and field systems).

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 32

171 172

Multi-substance indicators of bioavailable contamination

173

In order to represent the overall contamination pressure on aquatic ecosystems,

174

indicators of integrated bioavailable contamination (IBC) are proposed, one for metals

175

(IBCmetals) and one for POPs (IBCPOPs). The IBC calculated for each site studied takes into

176

account contamination levels observed in caged organisms in regard to a threshold

177

value determined for each compound (BBAC: bioavailable background assessment

178

concentrations), as defined by Besse et al. (2013)31 and Ciliberti et al. (2017)27. First, the

179

BBACs for all elements considered in the present study (Table 1–2) were determined by

180

implementing the same statistical approach described by Besse et al. (2013)31 with the

181

present national database. From this, the indicators IBCmetals and IBCPOPs computed a

182

given site only consider the compounds whose concentrations recorded in G. fossarum

183

significantly exceed their respective BBAC. Second, the concentrations were normalized

184

according to an average contamination level to avoid scale effects among compounds.

185

Hence, the IBC indicator is given as: 𝑛

∑𝑘 = 1 186

𝐼𝐵𝐶 =

𝑐𝑘 > 𝐵𝐵𝐴𝐶𝑘 𝜑𝑘 𝑛

187

where 𝑐𝑘 is the concentration of compound k measured in gammarids (only when the

188

measured concentrations exceed their respective BBAC), 𝜑𝑘 is the average

189

contamination level defined as the mean between BBAC value of compound k and the

190

90th percentile of concentrations of the compound measured over the whole dataset,

191

and 𝑛 is the number of analyzed compounds. A square-root transformation of compound

192

contributions was applied to reduce the asymmetry of the distribution of concentration

193

values between sites.

ACS Paragon Plus Environment

8

Page 9 of 32

Environmental Science & Technology

194 195

Link between chemical bioavailable contamination and macroinvertebrate

196

abundances

197

To inspect whether IBC indicators could predict effects on population densities at a

198

national scale, macroinvertebrate abundances recorded at corresponding sites should

199

be previously weighted by an ecoregional reference abundance level (i.e., carrying

200

capacity in undisturbed watercourses) in order to make the comparison possible

201

between geographically distant watercourses. Following the procedure previously

202

applied by Ciliberti et al. (2017)27, macroinvertebrate abundances recorded at the sites

203

defined in the French reference monitoring network (𝑛 = 392 with 3,127 sampling

204

events performed between 2005 and 2013) were grouped according to the physico-

205

chemical typology used in France. This typology classes streams based on alkalinity,

206

altitude, and nutrient conditions.34 Six classes make up this typology (hard water

207

streams: H1, H2, H3; and soft water streams: S1, S2 and S3). Reference levels of

208

abundance for a given taxon in H1, H2, H3, S1, S2, S3 streams were determined like the

209

95th percentile of the abundances recorded in the reference sites belonging to the six

210

classes respectively. The macroinvertebrate abundances recorded in our study sites

211

were then weighted by these reference levels which stand for a proxy of the optimal

212

carrying capacity of the respective physico-chemical class of each site. Data for six

213

macroinvertebrate families were collected and analyzed; Baetidae, Chironomidae,

214

Ephemeridae, Gammaridae, Hydrobiidae, and Hydropsychidae. Abundance values were

215

expressed at the family scale as a total of the genus level abundances in sampling units

216

of groups B and C so as to ensure fine representative sampling of the site. Following our

217

previous study,27 Gammaridae abundances were confronted either to IBCmetals or to

218

IBCPOPs. In a last step, taking advantage of IBCmetals as an integrated indicator of metal

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 32

219

contamination, the relationships between Baetidae, Chironomidae, Ephemeridae,

220

Hydrobiidae, and Hydropsychidae abundances and IBCmetals were assessed. These widely

221

distributed macroinvertebrate families were selected because they are regarded to

222

present contrasting toxicological sensitivities to metals in the ecotoxicological

223

literature.35–37

224

Quantile regression analyses were used to estimate the relationships between

225

IBC and changes in the density of macroinvertebrates. This type of regression analysis

226

developed to estimate rates of changes in all parts of the distribution of a response

227

variable is well suited to examine relationships between an environmental driver and a

228

descriptor of biological communities in cases where other influencing factors are

229

unmeasured and unaccounted for.38 Since other factors than chemical contamination

230

may also affect the macroinvertebrate abundances (e.g., habitat impairment, life history

231

disruption, loss or excess of food availability, etc.), a quantile regression model was

232

constructed to estimate the abundance of macroinvertebrate as a function of IBCmetals or

233

IBCPOPs by using the quantreg package39 of the R software. Quantile regression

234

coefficients were estimated for percentiles (𝜏) from 80% to 95% by increments of 5%.

235

Linear and curvilinear (exponential) models were built to each percentile and selection

236

procedure based on Akaike Information Criteria (AIC) was used to determine which

237

model (linear or curvilinear) best fits the data.

238 239

Results

240

Study site dataset

241

Out of 218 sites sampled and listed in the database, 211 sites were sufficiently

242

contextualized to be classified in the 6-class national general physicochemical classes

243

(H1: 𝑛 = 24, H2: 𝑛 = 70, H3: 𝑛 = 47, S1: 𝑛 = 8, S2: 𝑛 = 23, S3: 𝑛 = 39). These 211 sites

ACS Paragon Plus Environment

10

Page 11 of 32

Environmental Science & Technology

244

then were grouped in relation to watercourse size, geological nature of the watershed

245

soils, and land use types (Figure S1 of the Supporting Information, SI). The comparison

246

with sites of the entire SMN (𝑛 = 1,501 sites) shows that our set of sites was diverse and

247

represented with reliability the hydrographic network at the national scale. The

248

proportions of sites allocated to the different classes in the variables of agricultural land

249

cover and watercourse size conform with those calculated for the sites of the French

250

SMN (agricultural land cover: 𝜒2 = 1.11, 𝑝 = 0.57; watercourse size: 𝜒2 = 0.87, 𝑝 = 0.93).

251

In our database, the proportions of sites in the urbanization intensity classes are slightly

252

different than in the national SMN (𝜒2 = 7.72, 𝑝 = 0.02). Sites with a medium and high

253

urban land cover are, respectively, underrepresented and overrepresented (44% and

254

42% vs. 63% and 26% in the national SMN, respectively). The highest number of shifts

255

in the distribution of the sites among the classes is observed for the variable geology (𝜒2

256

= 18.52, 𝑝 = 9.50 × 10 ―5), with an underrepresentation of mixed profiles in our set of

257

sites (2% vs. 18% in the national SMN, respectively) while the proportion of sites with a

258

sedimentary profile is higher (80% vs. 56% in the national SMN, respectively).

259 260

National scale snapshot of bioavailable contamination in French watercourses

261

supplied by IBCmetals and IBCPOPs

262

The list of metal and POP compounds considered, respectively, in IBCmetals and IBCPOPs, is

263

presented in Tables 1 and 2. Considering metals, 17 compounds are quantified in more

264

than 90% of the sites, while one compound [mercury (Hg)] is quantified in only 45% of

265

sites (Table 1). For POPs, there is a more variable quantification among compounds,

266

with 15, 14, and 14 compounds quantified, respectively, on average in 90%, 62%, and

267

33% of sites (Table 2). The accumulation capacity of G. fossarum is assessed by the ratio

268

of maximal to minimal recorded concentrations. All compounds exhibited a ratio higher

ACS Paragon Plus Environment

11

Environmental Science & Technology

269

than 2, with values higher than 100 for some metal contaminants [chromium (Cr) and

270

nickel (Ni)]], some DDT isomers (4,4’-DDE), some PAHs (benzo(a)pyrene,

271

benzo(e)pyrene, chrysene, perylene), and some PCB congeners (n° 101, 138, 153 and

272

180).

273

Page 12 of 32

Overall, the geographic distribution of our study sites covers the whole territory,

274

but with an underrepresentation in the central area of France (Figure 1). Mapping

275

shows that weak and high metal contamination levels exhibited a homogeneous spatial

276

repartition at the territory scale (Figure 1a). Nonetheless, a striking case is the western

277

part of the Loire-Bretagne (LB) basin, where all sites exhibited contamination values

278

higher than the median of the national IBCmetals values, and 75% of them are higher than

279

the 75th percentile of the national IBCmetals values. Similar observations were made for

280

the hydro-ecoregion of Cévennes (situated on the south-western border of the Rhône-

281

Méditerranée-Corse (RMC) basin; 100% of IBCmetals values > median, 57% of IBCmetals

282

values > 75th percentile) and Vosges (situated on the eastern border of the Rhin-Meuse

283

(RM) basin; 100% of IBCmetals values > median, 83% of IBCmetals values > 75th percentile).

284

To a lesser extent, high values of IBCmetals are measured in the eastern area of the hydro-

285

ecoregion of Massif Central (situated on the eastern border of the LB basin; 62.5% of

286

IBCmetals values > median and 40% of them > 75th percentile). For POPs, 85% of sites

287

with a contamination value higher than the 75th percentile of the national IBCPOPs values

288

are located in the eastern part of France, RMC and RM basins, and north of France, in the

289

Artois-Picardie (AP) basin (Figure 1b). A spatial heterogeneity between metal and POP

290

contamination levels is shown by the lack of a correlation between the spatial

291

repartition of IBCmetals and IBCPOPs (𝑐𝑜𝑟 = ―0.036, 𝑝 = 0.655). For instance, sites

292

presenting higher values (> 90th percentile) of IBCmetals (respectively IBCPOPs) had null or

293

very low POP (or metals) contamination scores (Figure 2).

ACS Paragon Plus Environment

12

Page 13 of 32

Environmental Science & Technology

294 295

Co-variation between IBCmetals, IBCPOPs and gammarid population abundances

296

The relationships between the bioavailable chemical contaminations with the river-

297

typology-corrected abundance of gammarids are depicted in Figure 2. Curvilinear model

298

was the best fitted model of the relationships between the abundance of free-ranging

299

gammarids and the IBCmetals and ICBPOPs values. One significant reduction in the

300

abundance of free-ranging gammarids is apparent with the increase in values of metal

301

contaminations for all tested percentiles from 80% to 95% (Tables S1). The quantile

302

regression for 𝜏 = 90 is plotted on Figure 2a. A scatter plot revealed a similar trend

303

toward a decrease in the relative abundance of free-ranging gammarids as a function of

304

the increase in values of POP contaminations (Figure 2b), but the relationship in this

305

latter case with fewer study sites (n=196 sampling events) was not significant whatever

306

the considered percentiles (Tables S2).

307 308

Population abundances of six macroinvertebrate families toward the national

309

IBCmetals gradient

310

Figure 3 shows the distribution of river-typology-corrected abundances of six

311

macroinvertebrate families with a widespread distribution in French watercourses

312

(including gammarids) along the increasing gradient of IBCmetals values. For the two

313

families Ephemeridae and Hydrobiidae a significant negative relationship in quantile

314

regression analysis is detected between IBCmetals and corrected abundance as seen for

315

Gammaridae (Figure 3a, 3c, 3e). The best fitted model was curvilinear for Gammaridae

316

and Hydrobiidae, and linear for Ephemeridae (Table S1). For Hydrobiidae, regression

317

slopes were significant only for the highest percentiles (𝜏 = 90–95, Table S1). By

318

contrast, no significant relationship is detectable for Baetidae, Chironomidae, and

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 32

319

Hydropsychidae whatever the considered percentile (Figure 3b, 3d, 3f and Table S1). In

320

addition, the proportions of sites with zero abundance values are significantly different

321

among the six macroinvertebrate families in sites with the highest contamination scores,

322

for instance, the sites with IBCmetals > 0.39, i.e., the 90th percentile of the national

323

distribution (light gray/white areas in Figure 3, 𝜒2 = 35.08, 𝑝 = 1.452 × 10 ―6).

324 325

Discussion

326

Gammarid caging for monitoring the bioavailable contamination in freshwater

327

ecosystems

328

This study confirms and expands previous results that showed the relevance of caged G.

329

fossarum for monitoring chemical contamination in watercourses.27,31,40 Data from the

330

present study show that, during environmental exposure, a great diversity of

331

contaminants were accumulated in tissues of G. fossarum at concentrations exceeding

332

their respective BBAC (Tables 1 and 2). On the whole of contaminants quantified, 35

333

were identified like WFD priority substances. A total of 18 metals, whose eight had not

334

been detected in our previous study,31 were quantified. Also, 43 POPs were measured,

335

whose 24 compounds were different to Besse et al. (2013)31. For metal contaminants, in

336

common with our previous studies of the caged G. fossarum biomonitor, we report here

337

similar or higher frequencies of quantification (except for Hg27,31). For these same

338

metals, the ratios between maximal and minimal concentrations are higher in our study

339

than reported earlier,31 indicating that the enrichment of our database allows us to

340

cover a wider range of metal concentrations at the national scale. Regarding POPs, 19

341

contaminants are common to the regional study of Besse et al. (2013)31, and a difference

342

in the frequency of quantification is observed for only two contaminants (2,4’-DDT and

343

4,4’-DDD) for which we found lower frequencies, 14% and 46% instead of 100%,

ACS Paragon Plus Environment

14

Page 15 of 32

Environmental Science & Technology

344

respectively. This result suggests that the weaker quantification of these substances in

345

our national study, compared with the regional study, is explained by limited

346

contamination of the study sites by these contaminants rather than by low accumulation

347

in G. fossarum. These findings show the relevance of using the biota as a storage matrix

348

for metals and POPs to study the contamination levels, but also as an integration matrix

349

over time, giving a best representation of the contamination of aquatic ecosystems.

350 351

Spatial heterogeneity of bioavailable contamination revealed by IBCs at a national

352

scale

353

Numerous studies confirmed that contaminants co-occur in mixtures with many

354

chemical compounds in aquatic environments28,29,41 warranting the development of

355

multi-substance indicators. But the integration of many compounds into a single

356

indicator raises questions about how to manage the specific variability of each

357

compound, taking into account the background concentration levels of each and the

358

scale effect induced by ranges of contamination levels variable between compounds. By

359

considering in the IBC only the concentrations above the defined BBACs and by

360

weighting these concentrations according to the range of concentration values observed

361

at the national scale, we provided an approach allowing to graduate and compare

362

reliably the spatial and temporal levels of bioavailable contamination in streams.

363

The regional hotspots of high metal contamination observed in the LB basin, RM

364

basin, and RMC basin can be compared to the physicochemical characteristics of

365

watercourses. The water geochemistry influences the bioavailability of metal

366

contaminants (e.g., decrease in bioavailable of cadmium (Cd), nickel (Ni) and lead (Pb)

367

under raised water calcium concentration), which can reduce their bioaccumulation in

368

organisms.42–44 Also, biomonitoring integrates the influence of the water physico-

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 32

369

chemistry on the bioavailability of contaminants.45,46 Given that sites of these regions

370

belong to soft water classes of the physico-chemical typology, the low water hardness,

371

which makes metals relatively more bioavailable, explains the high metal bioavailable

372

contamination recorded thereinto. These findings reinforce the hypothesis that IBC,

373

because it is based on bioaccumulation data, exposes the vulnerability of communities

374

faced with exposure to metals in relation to national physicochemical typology.

375

Moreover, in the same hydro-ecoregion with soft water, profiles of contaminants vary

376

between study sites. For instance, in the western part of the LB basin, apart from the

377

cobalt (Co) and nickel (Ni) that were present (albeit at various concentrations) in all

378

study sites, cadmium (Cd) prevailed in five sites, while lead (Pb) and/or chromium (Cr)

379

prevailed in three other sites, and selenium (Se) in another site, suggesting different

380

sources of metal contamination. Similar reasoning can be followed to account for the

381

spatial heterogeneity in organic contamination levels. It can be argued that the

382

differentiated urbanization rates between regions explain the contrast in the

383

distribution of POP contamination among eastern and northern/western zones of the

384

territory, with contamination levels and urban land use higher in the latter (Figure S2),

385

while a more extensive urbanization releases a higher quantity of persistent organic

386

substances in aquatic ecosystems.47 Consequently, by converting a multitude of chemical

387

measurements into a single value, IBC indicators allow us to qualify sites relative to each

388

other and to identify those most subject to bioavailable contamination and/or the most

389

vulnerable, which can make them relevant prioritization tools for managers.

390 391

Relevance of multi-substance IBC indicators to decipher the effect of chemical

392

contaminations on macroinvertebrate populations

ACS Paragon Plus Environment

16

Page 17 of 32

Environmental Science & Technology

393

Previous studies showed a sensitivity of macroinvertebrate populations to bioavailable

394

chemical contaminants, but they considered one contaminant at a time.48–50 Dealing with

395

water concentrations of multiple substances, some studies have used alternative metrics

396

to evaluate the ecological effect of multiple contaminants, for instance, by using the axis

397

scores from an ordination method to express a contamination gradient51 or by

398

developing models such as the cumulative criterion unit (CCU, which corrects for the

399

hardness of the water and therefore does correct the metal availability because of water

400

quality52). However, the use of data of contaminant concentrations measured in the

401

water column may limit the ability to establish a link to ecological effects. Ignoring

402

bioavailability could partly explain the weak relationships between CCU and the

403

macroinvertebrate richness.53 As such, we evaluated the relevance of the indicators,

404

which integrate over one week the bioavailable chemical contamination measured in

405

transplanted organisms (G. fossarum), as proxies of freshwater ecosystem chemical

406

contamination to which biological communities are exposed.

407

Analyzing the on-site abundances of Gammaridae against the IBCmetals or IBCPOPs

408

values shows a negative effect of bioavailable contamination on populations (Figure 2).

409

This reinforces the findings of our previous work carried out separately on three

410

individual metal contaminants27 and confirms that the ecosystem chemical

411

contamination such as that proposed through IBC indicators is reflected by a response of

412

the population abundance of this sensitive species. IBC indicators express the influence

413

of bioavailable contamination as a gradient, in contrast to Ciliberti et al. (2017)27, who

414

depicted a binary response of Gammaridae abundance in relation to a threshold of

415

contamination by each single metal element. It offers the opportunity to establish a

416

quantitative link between the chemical contamination exposure and the impairment of

417

populations. The exposure to chemical contamination can take various forms; most of

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 32

418

the contaminated study sites displayed a complex pattern of bioavailable contamination

419

with concentrations above BBACs for many substances, simultaneously for metals and

420

organics. There was no correlation among the spatial repartition of metal and organic

421

contamination with the sites presenting the higher values (>90th percentile) of IBCmetals

422

(IBCPOPs) and having no or very low organic (metal) contamination (Figure 2); this

423

shows that the patterns of Gammaridae abundance degradation with increasing

424

bioavailable contaminations of metals or POPs are independent at this national scale.

425

Regarding the five other macroinvertebrate families, contrasting patterns were

426

reported in their response to the metal contamination of ecosystems (expressed from

427

IBCmetals, Figure 3). Two families (Ephemeridae, Hydrobiidae), in addition to

428

Gammaridae, were found to be sensitive to metal contamination, as expressed by a

429

significant negative relationship in quantile regression analysis, while the other three

430

families seem tolerant to metal contamination (i.e., Baetidae, Chironomidae,

431

Hydropsychidae). Interestingly, the family-specific sensitivity patterns observed for the

432

six families are consistent with results of earlier experimental studies in which

433

Ephemeridae and Hydrobiidae were described with Gammaridae as being more

434

sensitive to metal toxicity relatively to Baetidae, Chironomidae, and Hydropsychidae.35–

435

37

436

metals is yet recognized.54–56 Therefore the relative tolerance observed here in several

437

populations should certainly be explained by ecological factors such as demographic

438

recovery or adaptation abilities, also already documented in experimental studies.57–59 A

439

similar pattern was hypothesized to explain the resilience of gammarid populations

440

exposed to insecticides in agricultural landscapes.60,61

441 442

This statement should be nuanced for Baetidae for which physiological sensitivity to

These results validate the relevance of the caging gammarid approach to decipher the effect of ambient chemical contaminations on macroinvertebrate

ACS Paragon Plus Environment

18

Page 19 of 32

Environmental Science & Technology

443

populations, by considering the complex link between the source of contamination and

444

population exposure at a large spatial scale. Further work is currently in progress to

445

validate the hypothesis of the effect of chemical contaminations, expressed through IBC

446

indicators, on macroinvertebrate community assemblages, expressed through

447

community metrics. Another prominent aspect is the potential of these indicators as

448

explicative variables in pressure-impact models62 to better understand the link between

449

pressure and ecological status, which is an essential step for developing effective river

450

management plans and shaping future environmental policy.63

451

In conclusion, providing tools allowing for a reliable assessment of the chemical

452

risk of aquatic ecosystems remains a challenging issue. The successful application of the

453

proposed multi-substance indicators over an entire country reinforces that resorting to

454

biota for large-scale monitoring of bioavailable contaminants in continental waters

455

constitutes a highly valuable approach. Because the bioavailable contamination is

456

expressed relative to reference thresholds defined at a large spatial scale, such

457

indicators allow us to reliably quantify levels of bioavailable contamination in aquatic

458

ecosystems. Our findings demonstrate that they are relevant in tracing the chemical

459

pressure to which some components of biological communities are subjected, for

460

instance, by highlighting species-specific responses to chemical exposure gradients. It

461

paves the way for further studies to examine the link between aquatic ecosystem

462

exposure to chemicals and potential changes in the structure of macroinvertebrate

463

communities and the functions they sustain.

464 465

Acknowledgments

466

We are grateful for financial support from AFB (the French National Agency Grant

467

agreement for action program 2016-2018 for Water and Aquatic Ecosystems). This

ACS Paragon Plus Environment

19

Environmental Science & Technology

Page 20 of 32

468

work benefited also from the support of two programs: ANR CESA program GAMMA021

469

02 and Blanc program Multistress0004-03.

470 471

Authors’ contributions

472

ArC, OG and OP conceived the idea and designed methodology; ArC, AF and OG collected

473

the chemical data; AnC collected the land use data; MF, JP and BV gathered data on

474

macroinvertebrates; BA analyzed the data; BA led the writing of the manuscript. BA, ArC

475

and OG contributed critically to the drafts and all authors gave final approval for

476

publication.

477 478

Supporting Information

479

Figures showing the distribution of the whole 211 study sites according to their physical

480

state and land use data, and into eastern and northern/western regions of the French

481

territory. Tables showing the fit statistics for linear and curvilinear quantile regression

482

models describing the river-typology-corrected abundance organisms as a function of

483

the chemical bioavailable contamination (expressed from IBCmetals and IBCPOPs) for six

484

families of macroinvertebrates.

485

If paper is accepted for publication, supplementary data will be found online on the

486

website of the journal.

ACS Paragon Plus Environment

20

Page 21 of 32

Environmental Science & Technology

References (1)

(2) (3) (4) (5) (6)

(7) (8) (9)

(10)

(11) (12) (13) (14)

Bindler, R.; Rydberg, J.; Renberg, I. Establishing natural sediment reference conditions for metals and the legacy of long-range and local pollution on lakes in Europe. J. Paleolimnol. 2011, 45, 519–531. http://doi.org/10.1007/s10933-0109425-5 Fleeger, J. W.; Carman, K. R.; Nisbet, R. M. Indirect effects of contaminants in aquatic ecosystems. Sci. Total Environ. 2003, 317, 207–233. https://doi.org/10.1016/S0048-9697 Beketov, M. A.; Kefford, B. J.; Schäfer, R. B.; Liess, M. Pesticides reduce regional biodiversity of stream invertebrates. Proc. Natl. Acad. Sci. USA 2013, 110, 11039‒11043. https://doi.org/10.1073/pnas.1305618110 Kremen, C. Managing ecosystem services: what do we need to know about their ecology? Ecol. Lett. 2005, 8, 468–479. https://doi.org/10.1111/j.14610248.2005.00751.x Grizzetti, B.; Lanzanova, D.; Liquete, C.; Reynaud, A.; Cardoso, A.C. Assessing water ecosystem services for water resource management. Environ. Sci. Policy 2016, 61, 194–203. https://doi.org/10.1016/j.ensci.2016.04.008 Vörösmarty, C. J.; Lévêque, C.; Revenga, C. Fresh water. In Hassan, R.; Scholes, R.; Ash, N. (Eds.), Ecosystems and Human Well-Being: Current States and Trends, Millennium Ecosystem Assessment Report. (pp 165–207), USA: Washington, DC, Island Press, 2005 Hanna, D. E. L.; Tomscha, S. A.; Dallaire, C. O.; Bennett, E. M. A review of riverine ecosystem service quantification: research gaps and recommendations. J. Appl. Ecol. 2018, 55, 1299–1311. https://doi.org/10.1111/1365-2664.13045 European Council. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for community action in the field of water policy, 72pp, 2000 Van Ael, E.; Belpaire, C.; Breine, J.; Geeraerts, C.; Van Thuyne, G.; Eulaers, R.; Blust, R. Bervoets, L. Are persistent organic pollutants and metals in eel muscle predictive for the ecological water quality? Environ. Pollut. 2014, 186, 165–171. https://doi.org/10.1016/j.envpol.2013.12.006 Shahid, N.; Becker, J. M.; Krauss, M.; Brack, W.; Liess, M. Pesticide body burden of the crustacean Gammarus pulex as a measure of toxic pressure in agricultural streams. Environ. Sci. Technol. 2018, 52, 7823–7832. https://doi.org/10.1021/acs.est.8b01751 Hare, L. Aquatic insects and trace metals: bioavailability, bioaccumulation, and toxicity. Crit. Rev. Toxicol. 1992, 22, 327–369. https://doi.org/10.3109/10408449209146312 Oikari, A. Caging techniques for field exposures of fish to chemical contaminants. Aquat. Toxicol. 2006, 78, 370–381. https://doi.org/10.1016/j.aquatox.2006.03.010 De Kock, W. C.; Kramer, K. J. M. Active biomonitoring (ABM) by translocation of bivalve molluscs. In Kramer K. J. M. (Ed.), Biomonitoring of Coastal Waters and Estuaries. (pp 51–84), CRC Press, Boca–Raton, 1994 Smolders, R.; Bervoets, L.; Wepener, V.; Blust, R. A conceptual framework for using mussels as biomonitors in whole effluent toxicity. Hum. Ecol. Risk Assess. 2003, 9, 741–760. https://doi.org/10.1080/713609965

ACS Paragon Plus Environment

21

Environmental Science & Technology

(15)

(16) (17) (18) (19) (20)

(21)

(22)

(23)

(24)

(25)

(26)

(27)

Page 22 of 32

Westram, A. M.; Jokela, J.; Baumgartner, C.; Keller, I. Spatial distribution of cryptic species diversity in European freshwater Amphipods (Gammarus fossarum) as revealed by pyrosequencing. PLoS One 2011, 6, e23879. https://doi.org/10.1371/journal.pone.0023879 Weiss, M.; Leese, F. Widely distributed and regionally isolated! Drivers of genetic structure in Gammarus fossarum in a human-impacted landscape. BMC Evol. Biol. 2016, 16, 153. https://doi.org/10.1186/s12862-016-0723-z MacNeil, C.; Dick, J. T. A., Elwood, R.W. The dynamics of predation on Gammarus spp. (Crustacean: Amphidopa). Biol. Rev. 1999, 74, 375–395. https://doi.org/10.1111/j.1469-185X.1999.tb00035.x Kelly, D.; Dick, J. T. A.; Montgomery, W. I. The functional role of Gammarus (Crustacea, Amphipoda): shredders, predators, or both? Hydrobiologia, 2002, 485, 199–203. https://doi.org/10.1023/A:102137045349 Bizzotto, E. C.; Villa, S.; Vighi, M. POP bioaccumulation in macroinvertebrates of alpine freshwater systems. Environ. Pollut. 2009, 157, 3192–3198. https://doi.org/10.1016/j.envpol.2009.06.001 Tlili, K.; Labadie, P.; Bourges, C.; Desportes, A.; Chevreuil, M. Bioaccumulation of polybrominated diphenyl ethers by the freshwater benthic amphipod Gammarux pulex. Arch. Environ. Contam. Toxicol. 2012, 63, 69–76. https://doi.org/10.1007/s00244-012-9751-x Besse, J. P.; Geffard, O.; Coquery, M. Relevance and applicability of active biomonitoring in continental waters under the Water Framework Directive. Trends in Anal. Chem. 2012, 36, 113–127. https://doi.org/10.1016/j.trac.2012.04.004 Inostroza, P. A.; Vera-Escalona, I.; Wild, R.; Norf, H.; Brauns, M. Tandem action of natural and chemical stressors in stream ecosystems: insights from a population genetic perspective. Environ. Sci. Technol. 2018, 52, 7962–7971. https://doi.org.10.1021/acs.est.8b01259 Munz, N. A.; Fu, Q.; Stamm, C.; Hollender, J. Internal concentrations in gammarids reveal increased risk of organic micropollutants in wastewater-impacted streams. Environ. Sci. Technol. 2018, 52, 10347–10358. https://doi.org.10.1021/acs.est.8b03632 Fialkowski, W.; Fialkowska, E.; Smith, B. D.; Rainbow, P. S. Biomonitoring survey of trace metal pollution in streams of a catchment draining a zinc and lead mining area of Upper Silesia, Poland using the amphipod Gammarus fossarum. Internat. Rev. Hydrobiol. 2003, 88, 187–200. https://doi.org/10.1002/iroh.200390014 Luoma, S. N.; Cain, D. J.; Rainbow, P. S. Calibrating biomonitors to ecological disturbance: a new technique for explaining metal effects in natural waters. Integr. Environ. Assess. Manag. 2010, 6, 199–209. https://doi.org/10.1897/IEAM_2009.067.1 Rainbow, P. S.; Hildrew, A. G.; Smith, B. D.; Geatches, T., Luoma, S. N. Caddisflies as biomonitors identifying thresholds of toxic metal bioavailability that affect the stream benthos. Environ. Pollut. 2012, 166, 196-207. https://doi.org/10.1016/j.envpol.2012.03.017 Ciliberti, A.; Chaumot, A.; Recoura-Massaquant, R.; Chandesris, A.; François A.; Coquery, M.; Ferréol, M.; Geffard, O. Caged Gammarus as biomonitors identifying thresholds of toxic metal bioavailability that affect gammarid densities at the French national scale. Water Res. 2017, 118, 131–140. https://doi.org/10.1016/j.waters.2017.04.031

ACS Paragon Plus Environment

22

Page 23 of 32

Environmental Science & Technology

(28) (29)

(30)

(31)

(32)

(33) (34) (35)

(36) (37) (38) (39) (40)

(41)

Loos, R.; Gawlik, B. M.; Locoro, G.; Rimaviciute, E.; Contini, S.; Bidoglio, G. EU-wide survey of polar organic persistent pollutants in European river waters. Environ. Pollut. 2009, 157, 561–568. https://doi.org/10.1016/j;envpol.2008.06.020 Moschet, C.; Wittmer, I.; Simovic, J.; Junghans, M.; Piazzoli, A.; Singer, H.; Stamm, C.; Leu, C.; Hollender, J. How a complete pesticide screening changes the assessment of surface water quality. Environ. Sci. Technol. 2014, 48, 5423–5432. https://doi.org/10.1021/es500371t Mondy, C. P.; Villeneuve, B.; Archaimbault, V.; Usseglio-Polatera, P. A new macroinvertebrate-based multimetric index (I2M2) to evaluate ecological quality of French wadeable streams fulfilling the WFD demands: A taxonomical and trait approach. Ecol. Indic. 2012, 18, 452–467. https://doi.org/10.1016/j.ecolind.2011.12.013 Besse, J. P.; Coquery, M.; Lopez, C.; Chaumot, A.; Budzinski, H.; Labadie, P.; Geffard, O. Caged Gammarus fossarum (Crustacea) as a robust tool for the characterization of bioavailable contamination levels in continental waters: towards the determination of threshold values. Water Res. 2013, 47, 650–660. https://doi.org/10.1016/j.waters.2012.10.024 Wasson, J. G.; Chandersis, A.; Pella, H.; Sauquet, E.; Mengin, N. Scientific support for the implementation of the European Water Framework Directive. Typology of rivers in metropolitan France. Rapport, Ministère de l’Aménagement du Territoire et de l’Environnement, Cemagref Lyon BEA/LHQ, Lyon, 62pp, 2006 Büttner, G.; Kosztra, B. CLC2006 Technical Guidelines. European Environment Agency (Technical Report), 70 pp, 2007 Bougon, N.; Ferréol, M. Physico-Chimie soutenant la biologie. Typologie chimique nationale des cours d’eau. Irstea, 39 pp, 2011 Wogram, J.; Liess, M. Rank ordering of macroinvertebrate species sensitivity to toxic compounds by comparison with that of Daphnia magna. Bull. Environ. Contam. Toxicol. 2001, 67, 360–367. https://doi.org/10.1007/s00128-001-01338 Von der Ohe, P. C.; Liess, M. Relative sensitivity distribution of aquatic invertebrates to organic and metal compounds. Environ. Toxicol. Chem. 2004, 23, 150–156. https://doig.org/10.1897/02-577 Malaj, E.; Grote, M.; Schäfer, R. B.; Brack, W.; Von der Ohe, P.C. Physiological sensitivity of freshwater macroinvertebrates to heavy metals. Environ. Toxicol. Chem. 2012, 31, 1754–1764. https://doi.org:10.1002/etc.1868. Cade, B. S.; Terrell, J. W.; Schroeder, R. L. Estimating effects of limiting factors with regression quantiles. Ecology, 1999, 80, 311–323. https://doi.org/10.1890/0012-9658(1999)080[0311:EEOLFW]2.0.CO;2 Koenker, R., 2018. quantreg: Quantreg Regression. R package version 5.36. https://CRAN.R-project.org/package=quantreg Urien, N.; Lebrun, J. D.; Fechner, L. C.; Uher, E.; François, A.; Quéau, H.; Coquery, M.; Chaumot, A.; Geffard, O. Environmental relevance of laboratory-derived kinetic models to predict trace metal bioaccumulation in gammarids: field experimentation at a large spatial scale (France). Water Res. 2016, 95, 330–339. https://doi.org/10.1016/j.watres.2016.03.023 Busch, W.; Schmidt, S.; Kühne, R.; Schulze, T.; Krauss, M.; Altenburger, R. Micropollutants in European rivers: a mode of action survey to support the development of effect-based tools for water monitoring. Environ. Toxicol. Chem. 2016, 35, 1887–1899. https://doi.org/10.1002/etc.3460

ACS Paragon Plus Environment

23

Environmental Science & Technology

(42) (43) (44)

(45)

(46) (47)

(48) (49)

(50)

(51)

(52)

(53) (54)

Page 24 of 32

Qiu, J. W.; Xei, Z. C.; Wang, W. X. Effects of calcium on the uptake and elimination of cadmium and zinc in Asiatic clams. Arch. Environ. Contam. Toxicol. 2005, 48, 278–287. https://doi.org/10.1007/s00244-004-0001-8 Tan, Q. G.; Wang, W. X. The influences of ambient and body calcium on cadmium and zinc accumulation in Daphnia magna. Environ. Toxicol. Chem. 2008, 27, 1605– 1613. https://doi.org/10.1897/07-592.1 Pellet, B.; Geffard, O.; Lacour, C.; Kermoal, T.; Gourlay-Francé, C.; TusseauVuillemin, M. A model predicting waterborne cadmium bioaccumulation in Gammarux pulex: the effects of dissolved organic ligands, calcium, and temperature. Environ. Toxicol. Chem. 2009, 28, 2434–2442. https://doi.org/10.1897/09-042.1 McGeer, J. C.; Brix, K. V.; Skeaff, J. M.; DeForest, D. K.; Brigham, S. I.; Adams, W. J.; Green, A. Inverse relationship between bioconcentration factor and exposure concentration for metals: implications for hazard assessment of metals in the aquatic environment. Environ. Toxicol. Chem. 2003, 22, 1017–1037. https://doi.org/10.1022/etc.5620220509 Rainbow, P. S. Trace metal bioaccumulation: models, metabolic availability and toxicity. Environ. Intern. 2007, 33, 576–582. https://doi.org/10.1016/j.envint.2006.05.007 Dsikowitzky, L.; Schwarzbauer, J. Organic Contaminants from Industrial Wastewaters: Identification, Toxicity and Fate in the Environment. In Lichtfouse, E.; Schwarzbauer, J.; Robert, D. (Eds.), Pollutant Diseases, Remediation and Recycling. (pp 45–101), Switzerland: Springer International Publishing, 2013 Fialkowski, W.; Rainbow, P. S. The discriminatory power of two biomonitors of trace metal bioavailabilities in freshwater streams. Water Res. 2006, 40, 1805– 1810. https://doi.org/10.1016/j.watres.2006.02.035 Schmidt, T. S.; Clements, W. H.; Zuellig, R. E.; Mitchell, K. A.; Church, S. E.; Wanty, R. B., San Juan, C.A.; Adams, M.; Lamothe, P. J. Critical tissue residue approach linking accumulated metals in aquatic insects to population and community-level effects. Environ. Sci. Technol. 2011, 45, 7004–7010. https://doi.org/10.1021/es200215s Bervoets, L.; De Jonge, M.; Blust, R. Identification of threshold body burdens of metals for the protection of the aquatic ecological status using two benthic invertebrates. Environ. Pollut. 2016, 210, 76–84. https://doi.org/10.1016/j.envpol.2015.12.005 Kail, J.; Arle, J.; Jähnig, S. C. Limiting factors and threshold for macroinvertebrate assemblages in European rivers: Empirical evidence from three datasets on water quality, catchment urbanization, and river restoration. Ecol. Indic. 2012, 18, 63– 72. https://doi.org/10.1016/j.ecolind.2011.09.038 Clements, W. H.; Carlisle, D. M.; Lazorchak, J. M.; Johnson, P. C. Heavy metals structure benthic communities in Colorado mountain streams. Ecol. Appl. 2000, 10, 626–638. https://doi.org/10.1890/10510761(2000)010[0626:HMSBCI]2.0.CO;2 Hirst, H.; Jüttner, I.; Ormerod, S. J. Comparing the responses of diatoms and macroinvertebrates to metals in upland streams of Wales and Cornwall. Fresh. Biol. 2002, 47, 1752–1765. https://doi.org/10.1046/j.1365-2427.2002.00904.x Leland, H. V.; Fend, S. V.; Dudley, T. L.; Carter, J. L. Effects on species composition of benthic insects in a Sierra Nevada, California, stream. Fresh. Biol. 1989, 21, 163–179. https://doi.org/10.1111/j.1365-2427-tb01356.x

ACS Paragon Plus Environment

24

Page 25 of 32

Environmental Science & Technology

(55)

(56) (57) (58) (59) (60) (61) (62)

(63)

Cadmus, P.; Clements, W. H.; Williamson, J. L.; Ranville, J. F.; Meyer, J. S.; GutiérrezGinés, M. J. The use of field and mesocosm experiments to quantify effects of physical and chemical stressors in mining-contaminated streams. Environ. Sci. Technol. 2016, 50, 7825–7833. https://doi.org/10.1021/acs.est.6b01911 Mebane, C. A.; Schmidt, T. S.; Balistrieri, L. S. Larval aquatic insect responses to cadmium and zinc in experimental streams. Environ. Toxicol. Chem. 2017, 36, 749–762. https://doi.org/10.1002/etc.3599 Kiffney, P. M.; Clements, W. H. Effects of heavy metals on a macroinvertebrate assemblage from a Rocky Mountain stream in experimental microcosms. J. N. Am. Benthol. Soc. 1994, 13, 511–523. https://doi.org/10.2307/1467847 Mackay, R. J. Colonization by lotic macroinvertebrates: a review of processes and patterns. Can. J. Fish. Aquat. Sci. 1992, 49, 617–628. https://doi.org/10.1139/f92071 Clements, W. H. Benthic invertebrate community responses to heavy metals in the upper Arkansas River Basin, Colorado. J. N. Am. Benthol. Soc. 1994, 13, 30–44. https://doi.org/10.2307/1467263 Liess, M.; Schulz, R. Linking insecticide contamination and population response in an agriculture stream. Environ. Toxicol. Chem. 1999, 18, 1948–1955. https://doi.org/10.1002/etc.5620180913 Liess, M.; Von der Ohe, P. C. Analyzing effects of pesticides on invertebrate communities in streams. Environ. Toxicol. Chem. 2005, 24, 954–965. https://doi.org/10.1897/03-652.1 Villeneuve, B.; Piffady, J.; Valette, L.; Souchon, Y.; Usseglio-Polatera. Direct and indirect effects of multiple stressors on stream invertebrates across watershed, reach and site scales: a structural equation modelling better informing on hydromorphological impacts. Sci. Total Environ. 2018, 612, 660–671. https://doi.org/10.1016/j.scitotenv.2017.08.197 Hering, D.; Carvalho, L.; Argillier, C.; Beklioglu, M.; Borja, A.; Cardoso, A.C.; Duel, H.; Ferreira, T.; Globevnik, L.; Hanganu, J.; Hellsten, S.; Jeppesen, E.; Kodeš, V.; Solheim, A.L.; Nõges, T.; Ormerod, S.; Panagopoulos, Y.; Schmutz, S.; Venohr, M.; Birk, S. Mapping aquatic ecosystems and water resources under multiple stress – An introduction to the MARS project. Sci. Total Environ. 2015, 503–504, 10–21. http://doi.org/10.1016/j.scitotenv.2014.06.106

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 32

Figures Figure 1. Distribution throughout France of sites where active biomonitoring (caged G. fossarum) is applied to evaluate the level of (a) bioavailable metal contamination (𝑛 = 218 sites) and (b) bioavailable POP contamination (𝑛 = 165 sites). Dot color corresponds to the chemical contamination expressed either from IBCmetals values or IBCPOPs values. The color gradient (metals: from green to blue, POPs: from green to red) indicates growing level of chemical contamination. When IBCmetals equals zero (green points), all measured values for each contaminant are above their respective BBACs. Geographic delimitations on the map represent the boundaries of six French regional environmental agencies; Adour-Garonne (AG), Artois-Picardie (AP), LoireBretagne (LB), Rhin-Meuse (RM), Rhône-Méditerranée-Corse (RMC).

ACS Paragon Plus Environment

26

Page 27 of 32

Environmental Science & Technology

Figure 2. River-typology-corrected abundance of free-ranging Gammaridae (expressed as a percentage of the 95th percentile of the abundance observed on reference sites with the same physico-chemical typology) as a function of (a) the metal bioavailable contamination (IBCmetals) and (b) the POP bioavailable contamination (IBCPOPs). When IBCmetals or IBCPOPs equal zero (green points in both panels), all measured values for each contaminant are above their respective BBACs. In the top panel, dot colors correspond to the chemical contamination expressed from IBCPOPS. Black dots correspond to sampling sites for which data of POPs contamination are not available. In the bottom panel, dot colors correspond to the chemical contamination expressed from IBCmetals. Gray lines indicate regressions fitted at 90th percentile distributions. Significant relationships are shown with solid lines, and non-significant relationships are shown with dashed lines. According to the AIC value, curvilinear model was the best to fit the relationships between relative abundance of Gammaridae and IBCmetals or IBCPOPs values (Tables S1, S2). 𝑝 values were computed by bootstrapping with 200 replications.

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 32

Figure 3. River-typology-corrected abundance of free-ranging organism estimates as a function of the metal bioavailable contamination expressed from the computation of IBCmetals for six families of macroinvertebrates: (a) Gammaridae, (b) Chironomidae, (c) Hydrobiidae, (d) Hydropsychidae, (e) Ephemeridae, and (f) Baetidae. Points in the white part of the scatter plot exhibit a value of metal contamination higher than the 90th percentile of the distribution of IBCmetals values. Gray lines indicate regressions fitted at 90th percentile distributions. Significant relationships are shown with solid lines, and non-significant relationships are shown with dashed lines. According to the AIC value, curvilinear model was the best to fit the relationships between relative abundance of Gammaridae (a), Hydrobiidae (c), and Baetidae (f) and IBCmetals vlaues, while linear model was the best for Ephemeridae (e), Chironomidae (b), and Hydropsychidae (d) (Table S1). 𝑝 values were computed by bootstrapping with 200 replications.

ACS Paragon Plus Environment

28

Page 29 of 32

Environmental Science & Technology

Tables Table 1. Concentrations of metals, limits of quantification (LQ) and frequency of quantification (% data ≥ LQ) in caged G. fossarum after 7 days of exposure for the 218 sites. Concentrations and LQ are expressed in µg g-1 dw (dry weight). Compounds Concentration (µg g-1 dw) Frequency (%) Mean Med Min 90th Max LQ Agb 0.18 0.14 0.04 0.301 1.02 0.02 90 Al 600.24 463.40 36.10 1300 2311.00 2.00 100 B 3.99 3.54