Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review

Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review

Journal Pre-proof Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review Andrea Castaño-Sánchez, Grant C. Hose, Ana S...

4MB Sizes 1 Downloads 53 Views

Journal Pre-proof Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review Andrea Castaño-Sánchez, Grant C. Hose, Ana Sofia P.S. Reboleira PII:

S0045-6535(19)32662-1

DOI:

https://doi.org/10.1016/j.chemosphere.2019.125422

Reference:

CHEM 125422

To appear in:

ECSN

Received Date: 4 September 2019 Revised Date:

18 November 2019

Accepted Date: 19 November 2019

Please cite this article as: Castaño-Sánchez, A., Hose, G.C., Reboleira, A.S.P.S., Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review, Chemosphere (2019), doi: https://doi.org/10.1016/j.chemosphere.2019.125422. This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2019 Published by Elsevier Ltd.

CRediT author statement Andrea Castaño-Sánchez: Conceptualization, Methodology, Visualization, Data Curation, Writing – Original Draft, Investigation, Formal Analysis. Grant Hose: Visualization, Writing – Review & Editing. Ana Sofia Reboleira: Conceptualization, Methodology, Data Curation, Resources, Investigation, Formal Analysis, Funding Acquisition, Project Administration, Supervision.

1

Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review

2 3

Andrea Castaño-Sáncheza; Grant C. Hoseb; Ana Sofia P.S. Reboleiraa *

4

a

5

2100 Copenhagen, Denmark

6

b

Department of Biological Sciences, Macquarie University, NSW 2109 Sydney, Australia

7

*

Corresponding author: [email protected]

Natural History Museum of Denmark, University of Copenhagen, Universitetsparken 15,

8 9

Abstract

10

How anthropogenic stressors affect biodiversity is a central question in a changing world.

11

Subterranean ecosystems and their biodiversity are particularly vulnerable to change, yet,

12

these species are frequently neglected in analyses of global biodiversity and assessments of

13

ecological status and risk. Are these hidden species affected by anthropogenic stressors? Do

14

they survive outside of the current thermal limits of their ecosystems? These and other

15

important questions can be addressed with ecotoxicological testing, relating contaminants and

16

temperature resistance of species with measured environmental concentrations and climatic

17

data. Ecotoxicological knowledge specific to subterranean ecosystems is crucial for

18

establishing thresholds for their protection, but such data are both scarce and scattered. Here,

19

we review the existing ecotoxicological studies of these impacts to subterranean-adapted

20

species. An effort that includes 167 measured endpoints and presents a database containing

21

experimentally derived species’ tolerance data for 28 contaminants and temperature, for 46

22

terrestrial and groundwater species, including fungi and animals. The lack of standard data

23

among the studies is currently the major impediment to evaluate how stressors affect

24

subterranean-adapted species and how differently they respond from their relatives at surface. 1

25

Improving understanding of ecotoxicological effects on subterranean-adapted species will

26

require extensive analysis of physiological responses to a wide range of untested stressors,

27

standardization of testing protocols and evaluation of exposures under realistic scenarios.

28 29

Keywords

30

Subsurface ecosystems; human impacts; aquifer; ecotoxicological assays; stygofauna;

31

troglofauna.

32 33

1. Introduction

34

While most of the world's biodiversity is found at or very near the surface, numerous species

35

also live deep in the underground (Danielopol et al., 2003; Humphreys, 2006). The

36

subterranean domain harbors a peculiar ecosystem, containing narrowly-distributed species

37

with unique morphological and physiological traits, providing an aphotic habitat both in air

38

and water filled spaces below the surface (Humphreys, 2006; Dole-Olivier et al., 2009; Culver

39

and Pipan, 2019). As ecosystems, they range from shallow subterranean habitats, to deep

40

caves or meso and macrofissures that may be inundated by fresh, brackish or marine waters

41

(Culver and Pipan, 2019). As depth below ground increases, the influence of surface

42

processes decreases, and environmental conditions become increasingly constant (stable

43

temperature, high humidity and permanent darkness) (Culver and Pipan, 2019). These stable

44

conditions support diverse communities of highly adapted organisms (mainly invertebrates)

45

that have colonized this environment at different times (Gibert and Deharveng, 2002). These

46

subterranean ecosystems constitute one of the most poorly known and unprotected natural

47

resources of our planet (Pipan et al., 2010).

2

48

Access to subterranean environments is often limited to artificial (wells and boreholes) and

49

natural points (caves and springs) (Sket, 2018), which is reflected in the relatively limited

50

sampling of subterranean environments globally (Schneider and Culver, 2004; Culver et al.,

51

2006). Studies of more accessible ecotone habitats such as the aquatic hyporheic zone and the

52

terrestrial mesovoid shallow substrate (MSS) have provided valuable insights although the

53

extent of knowledge of deeper subterranean habitats is patchy. Indeed, the difficulty of

54

accessing subterranean ecosystems combined with the strong geographical biases on species

55

distribution have been recently described as the “Racovitzan impediment” (Ficetola et al.,

56

2019).

57

The infiltration of contaminants into the subsurface can be rapid and difficult to manage,

58

making subterranean ecosystems particularly vulnerable to pollution (Marmonier et al., 2013,

59

2018). Despite worldwide recognition that subterranean ecosystems are both the most

60

important sources of freshwater for human consumption, and critically endangered

61

ecosystems, initiatives like the EU Water Framework Directive (EU, 2000), the Groundwater

62

Directive (EU-GWD, 2006), the Safe Drinking Water Act of the US EPA (HDR, 2017), the

63

guidelines for drinking-water quality from the World Health Organization (WHO, 2011), or

64

the Australian Drinking Water Guidelines (NRMMC, 2011), focus only on water quality and

65

the need to achieve a good physicochemical status of groundwater, neglecting protection to its

66

endemic biodiversity (Reboleira et al., 2013). Subterranean biota play a key role in functions

67

that provide ecosystem services, including the maintenance of water quality and the support

68

of groundwater dependent ecosystems, such as springs and rivers (Griebler et al., 2014).

69

Therefore, the evaluation of the condition of subterranean ecosystems should consider both

70

abiotic and biotic components. Furthermore, the fact that terrestrial subterranean ecosystems

71

are intimately linked with the groundwater cycle is also neglected and scientific information 3

72

concerning the effect of pollution in these ecosystems is needed for their protection.

73

Several authors suggested that the evolutionary adaptations of subterranean biota make them

74

particularly sensitive to anthropogenic impacts (Griebler et al., 2010; Griebler and Avramov,

75

2015; Jiménez Valverde et al., 2017). The primary, global human impact in subterranean

76

ecosystems is the extraction of water leading to the depletion of aquifers. Around 50% of the

77

world’s population relies on aquifers for drinking water (World Water Assessment

78

Programme, 2003), and this is likely to increase under the predicted global warming and

79

climate change scenarios (Danielopol et al., 2003; Kløve et al., 2014). In parallel with the

80

over extraction of water, subterranean ecosystems are threatened by anthropogenic pollutants

81

(Danielopol et al., 2003; Marmonier et al., 2013).

82

Ecotoxicological dose-response bioassays are an important line of evidence for environmental

83

risk assessments (ERA) because they provide experimental evidence of cause and effect. A

84

large number of bioassays have been developed to assess the toxicity of contaminants to

85

plants, animals and microbes in terrestrial, marine and freshwater environments (Breitholtz et

86

al., 2006). However, the development of bioassays for subterranean organisms has been

87

limited. The first reported dose-response toxicity tests using subterranean species was

88

conducted in 1976 (Bosnak and Morgan, 1981), and since then, a small number of studies

89

have been done using aquatic subterranean organisms. While there is an extensive body of

90

literature on bioassays for superficial soil biota (Sochová et al., 2006; Van Gestel, 2012;

91

Duarte et al., 2017), there have been no bioassays for assessing the toxicity of pollutants using

92

terrestrial subterranean-adapted organisms.

93

Here we provide a critical review of the ecotoxicological effects of anthropogenic stressors on

94

subterranean organisms, and examine the state of the art of ecotoxicological bioassays with

95

subterranean organisms, including those associated with climate change. We synthesized 4

96

available data on the ecotoxicological effects known for subterranean species, focusing on the

97

type of stressors, studied organisms and experimental conditions, and provide future

98

perspectives for ecotoxicology on subterranean ecosystems. This constitutes a fundamental

99

starting point for the development of standard methods for estimating the effect of stressors

100

and climate change in these hidden ecosystems, definitely contributing to their protection.

101 102 103

2. The subterranean ecosystem 2.1 Environmental constraints

104

Conditions in subterranean ecosystems are typified by the absence of light, high air humidity,

105

oligotrophic (low carbon and nutrient) conditions (Schmidt and Hahn, 2012), and stable

106

temperatures (relative to surface environments); overall, physical and chemical conditions are

107

relatively stable compared to surface systems (Gibert and Deharveng, 2002). This suite of

108

conditions are common across all types of subterranean ecosystems (Culver and Pipan, 2019).

109

Subterranean ecosystems occur in highly complex and heterogeneous geological settings

110

(Mösslacher and Notenboom, 1999) including unconsolidated (e.g. alluvial sediments), and

111

consolidated rocks (e.g. karst). The connectivity within and between subterranean habitats,

112

and between them and the surface is determined by the connectivity among voids of the

113

matrix. Habitats can thus range from those that are highly and homogeneously porous and

114

well connected to the surface and adjoining habitats, to those that are highly complex and

115

fragmented with isolated pockets, and limited vertical and lateral connectivity. This

116

architecture of spatial partitioning and high isolation presents a pattern similar to island

117

biogeography, and has resulted in strong spatial heterogeneity and patchiness in species

118

distributions (Dole-Olivier et al., 2009; Malard et al., 2009).

5

119

The size of voids within the geological matrix, and the environmental conditions have been

120

primary selection forces for the species that have successfully colonized subterranean

121

ecosystems (e.g. Korbel and Hose, 2015). The size of the voids in the subterranean matrix

122

limit the size of biota that can inhabit it (Korbel et al., 2019). As a consequence, vertebrates

123

are relatively uncommon and biota is dominated by microbes (such as bacteria, fungi and

124

protozoa) and invertebrates (macro and meiofauna) (Gibert and Deharveng, 2002; Griebler

125

and Lueders, 2009; Pipan and Culver, 2017).

126 127

2.2 Simplified trophic chain

128

The subterranean food web is often strongly truncated, basally and apically, due to the

129

absence of light and consequently of photosynthetic producers, and potentially low abundance

130

of strict predators (Fig. 1) (Gibert and Deharveng, 2002). With the general absence of

131

vertebrates, invertebrates are often the highest trophic level. In the absence of photosynthetic

132

organisms, subterranean ecosystems are largely reliant on carbon infiltrating from the surface,

133

with few exceptions where there is significant primary production from chemolithotrophic

134

bacteria (Sarbu et al., 1996; Kumaresan et al., 2014; Brankovits et al., 2017). As a

135

consequence, carbon is often limited and microbial diversity and activity are typically low

136

(Griebler and Lueders, 2009). Microbial diversity found in the subsurface is represented by

137

Archaea, Bacteria, Protozoa and Fungi, which remain largely unknown (Griebler and Lueders,

138

2009). The scarcity of carbon and nutrients constrains population densities (biomass) and has

139

forced species to adapt and cope with a low energy environment by decreasing metabolic

140

rates and enhancing food detection through chemoreception and well-developed sensory

141

appendages (Di Lorenzo et al., 2015a). Also, functional redundancy, i.e., several species that

142

fulfill the same role (food web) within the ecosystem might occur (Gibert and Deharveng, 6

143

2002; Korbel and Hose, 2011; Saccó et al., 2019a), however specific feeding traits have been

144

observed in the same trophic level, even within the same genus (Arnedo et al., 2008; Hutchins

145

et al., 2014; Francois et al., 2015).

146 147

2.3 Ecological classification of biota

148

The environmental constraints of the subterranean domain provide strong selective pressures

149

and have led to similar morphological, physiological and life-history traits across independent

150

evolutionary lineages of organisms (Notenboom et al., 1994; Hancock et al., 2005;

151

Humphreys, 2006). These specialized traits have been named troglomorphisms or

152

“troglomorphic syndrome” (Christiansen, 2012). The ecological classification of subterranean

153

species is based on their degree of dependence on the subterranean ecosystems (White et al.,

154

2019). Species that complete their whole life-cycle in the subterranean environment and

155

typically exhibit typical morpho-physiological traits to the subterranean environment (e.g.,

156

lack of pigment and ocular regression, body and appendages elongation, hypertrophy of

157

sensory organs, low metabolic rate and prolonged life span with low reproductive output) are

158

referred to as troglobiont (terrestrial) or stygobiont (aquatic). The suffix -phile

159

(troglophile/stygophile) is used for those species with affinities to subterranean ecosystems,

160

but can also complete their life-cycle at the surface. Finally, the suffix -xene

161

(trogloxene/stygoxene) is used to describe species that are present sporadically in

162

subterranean ecosystems but are unable to establish a subterranean population; they are

163

considered occasional inhabitants (Sket, 2008).

164 165

3. Anthropogenic stressors on the subterranean environment

7

166

Anthropogenic pressures have surpassed natural forces as drivers of change in Earth’s

167

environment (Steffen et al., 2007). Subsoil and groundwater contamination, groundwater

168

pumping and mining activities threaten subterranean ecosystems and the provision of

169

ecosystem services (Danielopol et al., 2003; Larned, 2012) (Fig. 2). Moreover, the changes to

170

the atmosphere and surface environments associated with climate change (warming effect)

171

will result in changes to the subsurface quality (Dijck et al., 2006) and its annual thermal

172

regime, but the lack of knowledge makes impossible to determine the magnitude and direction

173

to climate related trends for subterranean ecosystems (Kundzewicz and Döll, 2009).

174

The main cause of subterranean pollution is by infiltration from surface due to industry, urban

175

discharge (urban water disposals), transport, agriculture and farming (Marmonier et al., 2013).

176

Contaminants frequently detected in subterranean ecosystems include metals, pesticides,

177

fertilizers, emerging contaminants, and volatile organic compounds (VOCs) (Fig. 3, Table

178

S.1). Measured environmental concentrations (MECs) at which these contaminants are

179

measured in groundwater can be found in databases such as EEA (2019) and the National

180

Water Quality Monitoring (2019), among others. However, similar information for sediment

181

contamination in the subterranean terrestrial and aquatic compartments are still missing. The

182

impact of contaminants depends on the source, physicochemical characteristics, mobility and

183

behavior in soil or the aqueous environment and associated hazards (Stuart et al., 2012). In the

184

following sections we examine the current extent of common contaminants in subterranean

185

ecosystems.

186 187

3.1 Metals

188

Metal (and metalloid) contamination of subterranean and surface environments has received

189

more attention than other contaminants because of the long history of metals use (Momodu 8

190

and Anyakora, 2010; Kumar et al., 2012) (Fig. 3A). The toxicity, biological activity,

191

bioavailability, solubility, mobility, reactivity and elimination of a metal can vary depending

192

on its speciation form (Santos et al., 2002), which itself depends on the physicochemical

193

characteristics of the environment (such as pH, redox potential, temperature, complexation

194

with other dissolved constituents, sorption and ion-exchange capacity of the geological

195

materials, and organic matter content) (Evanko and Dzombak, 1997; Hashim et al., 2011).

196

Accordingly, contamination by metals in karst areas needs special consideration, since the

197

alkalinity (> 7 pH) and the carbonates may remove metals from the dissolved phase (Vesper,

198

2019). The formation of metal-carbonate precipitates and secondary minerals reduces the

199

metal solubility and transport. The resulting chemical species will have different toxicity and

200

bioavailability (White et al., 2019).

201

Anthropogenic sources of metals in subterranean environments include point sources such as

202

mining activities (in which contamination may continue long after mining activities have

203

ceased (Hutton and Symon, 1986; Nriagu, 1989; Duruibe et al., 2007)), biosolid and sludge

204

applications to land, and waste water discharges, as well as diffuse sources such as urban and

205

agricultural runoff (Evanko and Dzombak, 1997). Diffuse metal contamination is a significant

206

problem in some agricultural areas (Vodela et al., 1997; Wongsasuluk et al., 2014) since

207

metals are common constituents of pesticides and fertilizers (Santos et al., 2002).

208 209

3.2 Pesticides

210

Pesticides detected in subterranean ecosystems (Fig. 3B) include synthetic chemicals and

211

natural substances that can be classified according to their specific biological activity on

212

target species (Van der Hoff and Van Zoonen, 1999; Ahmed, 2001). As a result, this group

213

includes chemicals with a wide variety of uses, properties and potential impacts on the 9

214

environment (Stuart et al., 2012). Over 500 compounds are currently registered worldwide as

215

pesticides, or metabolites of pesticides (Van der Hoff and Van Zoonen, 1999; Ahmed, 2001;

216

Loos et al., 2010).

217

Pesticides are considered a diffuse source of pollution for subterranean environments since

218

they are typically applied over broad geographical areas (Stuart et al., 2012). Most pesticides

219

infiltrate by percolating through the soil profile, or in the case of caves, mobilization of

220

contaminated soil from the surface. The mobility of a pesticide in the environment depends

221

heavily on its chemical properties, particularly its solubility in water and its affinity for

222

binding to soil and organic matter (Arias-Estévez et al., 2008). Degradation of pesticides is

223

governed by both abiotic and biotic processes, with the degradation products also biologically

224

active and often equally or more toxic than the parent compound (Arias-Estévez et al., 2008;

225

Stuart et al., 2012).

226

Studies have shown that shallow groundwater is more vulnerable to contamination from

227

pesticides than deep groundwater, because of its proximity to the surface (Haarstad and

228

Ludvigsen, 2007; Close and Skinner, 2012; Toccalino et al., 2014). Moreover, Stuart et al.

229

(2012) emphasize the importance of the rapid transport routes (like drains, sinkholes and

230

fractures), which are especially common in karst areas, that by-pass the main natural pesticide

231

attenuation areas in the vadose zone (Jiménez-Sánchez et al., 2008).

232

3.3 Fertilizers

233

The largest source of nutrients in groundwater is nitrogen and phosphorous – based fertilizers

234

used in agriculture (Fig. 3C). These can be synthetic chemical fertilizers or organic-waste

235

fertilizers (e.g. compost and blood or bone meal) which contain essential and nonessential

236

elements as metals (Khan et al., 2018). Both nitrogen and phosphorous are essential macro-

237

nutrients for plant growth and are thus critical for primary production in the biosphere; 10

238

however, when present at concentrations exceeding natural background levels they can cause

239

eutrophication problems (Gruber and Galloway, 2008).

240

Most of the world’s phosphate fertilizers come from phosphate rocks which also frequently

241

contain hazardous elements such as Cd, Pb, Hg, U, Cr and As. The application of phosphate

242

fertilizers can promote the contamination of agricultural soils by introducing these hazardous

243

metals and enhancing the mobility of metals already in the soils (Campos, 2002; Dissanayake

244

and Chandrajith, 2009; Khan et al., 2018).

245

The production of ammonia by the reaction of nitrogen and hydrogen, known as the Haber-

246

Bosh process, has enabled humans to approximately double the input of available nitrogen to

247

the Earth’s land surface, and greatly increase food production through the use of nitrogenous

248

fertilizers. However, this increase in nitrogen in a readily bioavailable form has also created

249

environmental issues (Galloway et al., 2008; Gruber and Galloway, 2008). Nitrogen is a

250

particular threat to groundwaters because of its high solubility in its various forms.

251

Background concentrations of nitrate in groundwater are typically below 2 mg/L,

252

consequently concentrations over 3 mg/L are considered indicative of anthropogenic inputs

253

(Burkart and Stoner, 2008).

254

Factors such as climate, soil properties, groundwater fluxes and management practices

255

influence the fertilizer load entering subterranean ecosystems (Wang et al., 2010). However,

256

the relationships between cropping systems, fertilizer demand, soil nitrate and phosphate

257

accumulation and groundwater fertilizers concentrations are still not fully understood (Ju et

258

al., 2006). Other nutrient sources in groundwater, include leakage from sewage systems and

259

septic tanks, industrial spillages, landfills, river or channel infiltration, and fertilizers used in

260

gardens (Wakida and Lerner, 2005). As for pesticides, shallow water-table aquifers and karst

11

261

systems are more vulnerable to nutrient contamination (Bijay-Singh et al., 1995; Boyer and

262

Pasquarell, 1995; Jiménez-Sánchez et al., 2008).

263 264

3.4 Volatile organic compounds (VOCs)

265

Volatile organic compounds (VOCs) are among the most frequently detected contaminants in

266

groundwaters (Squillace et al., 2002; EEA, 2007), usually at concentrations in the order of

267

µg/L (Fig. 3D). They are carbon-based compounds that evaporate readily at normal

268

temperatures (Fan et al., 2009) and include some oxygenated compounds, such as ethers or

269

aromatic hydrocarbons from mineral oil (e.g. benzene, toluene), and chlorinated aliphatic

270

hydrocarbons (Rivett et al., 2011).

271

Due to their volatility, VOCs enter the atmosphere in large quantities from anthropogenic

272

(e.g. biomass burning, oil and gas industry emissions) and biogenic sources (e.g. natural

273

wetland emission and domestic ruminants) (Atkinson, 2000; Fan et al., 2009). The

274

atmospheric deposition is the main non-point/diffusive source of VOCs in subterranean

275

ecosystems whereas point sources are primarily industrial or waste disposal sites where the

276

use and production or inadequate storage or disposal contribute VOCs to the soil or

277

groundwater (Pankow et al., 1997; Fan et al., 2009; Stephenson et al., 2013). Once in the soil,

278

VOCs can undergo diffusive transport through the unsaturated zone as a gas-phase or

279

infiltrate into the groundwater. In groundwater, VOCs may act as a non-aqueous phase liquid

280

(NAPL) due to their limited solubility in water, allowing them to remain in a separate phase.

281

Aromatic hydrocarbons often form light non-aqueous-phase liquids (LNAPL) and chlorinated

282

aliphatic hydrocarbons form denser non-aqueous-phase liquids (DNAPLs) (White et al., 2019;

283

Rivett et al., 2011). Once in groundwater, the DNAPLs, dissolve very slowly creating plumes

284

that persist over decades (Loop and White, 2001; Ellis and Rivett, 2007). 12

285

Field studies have shown changes in subterranean communities in response to VOC

286

contamination (e.g. Stephenson et al., 2013), but few laboratory-based studies have been

287

conducted. Due to their physical and (bio) chemical features such as volatility, limited

288

solubility in water, biodegradability, as well as sorption and bioaccumulation potential, VOCs

289

“present difficulties in the execution and interpretation of standard aquatic toxicity tests”

290

(Avramov et al., 2013).

291 292

3.5 Emerging contaminants

293

Emerging contaminants are substances that may have been present for a long time in the

294

environment, but whose presence is now being detected continuously at subtoxic

295

concentrations in groundwater (Fig. 3E) and can cause unexpected consequences in humans

296

and ecosystems (Daughton, 2005). Emerging contaminants include pharmaceuticals, personal

297

care

298

perfluorooctanesulfonic acid compounds (PFOS) (see Barceló, 2003).

299

The behavior of emerging contaminants in the subsoil, their solubility or ability to travel

300

through the aqueous environment, and persistence in the subsurface environment is generally

301

poorly known (Stuart et al., 2012). They occur mostly as mixtures, meaning that organisms

302

are usually exposed to complex multi-component mixtures over their whole life cycle

303

(Barceló, 2003).

products,

UV

filters/sunscreen,

food

additives,

flame/fire

retardants

and

304 305

3.6 Temperature

306

Subterranean-adapted organisms are typically exposed to small thermal fluctuations in their

307

environment. They are considered mostly ectothermic with a thermal tolerance is within the

308

range of conditions experienced in their natural habitat (Addo-Bediako et al., 2000), however 13

309

some stygobionts exhibit eurythermal characteristics (Issartel et al., 2005). Increasing

310

temperatures associated with climate change may affect the quality of subterranean

311

environments in two ways (Kløve et al., 2014). Firstly, expected changes in recharge rates

312

may affect the mobilization of compounds in the unsaturated zone, and, secondly, changes to

313

thermal regimes have important implications for the rates of temperature-dependent reactions

314

that directly affect many biogeochemical processes (e.g. redox reactions linked to nitrogen

315

and carbon cycles, and the speciation and mobilization of contaminants). These effects all

316

pose a risk to subterranean biodiversity.

317 318

4. Effects of anthropogenic stressors on subterranean organisms

319

Ecotoxicological testing is a useful tool for assessing the tolerance of organisms to

320

contaminants that threaten subterranean environments (Avramov et al., 2013). Bioassays are

321

based on dose/exposure-response assessment, which quantify the relationship between dose or

322

degree of exposure to a substance and the toxic effect or response (EC, 2003). They span a

323

continuum of biological complexity and ecological realism, from small, rapid and highly

324

reproducible tests on populations under laboratory conditions to long-term, large scale,

325

complex field-based community studies (Breitholtz et al., 2006). Bioassays vary in duration,

326

typically considered as either: a) acute assays – short-term exposure test, used to observed the

327

effects of toxicants under a single event exposure, where lethality is the most common

328

endpoint; or b) chronic assays – longer-term tests, with (usually) lower and continuous

329

contaminant exposure, mostly used to estimate sub-lethal endpoints related to the organism’s

330

life span (population effect) (Van Leeuwen et al., 1996).

331

The toxicity values reported in ecotoxicological bioassays, together with the environmental

332

contaminant exposure assessment are key lines of evidence in ERA (EC, 2003). Importantly, 14

333

the toxicity data used in ERA should reflect the species likely to be affected (EC, 2003;

334

Breitholtz et al., 2006; Kulkarni et al., 2013). The current lack of such data is a critical

335

knowledge gap for determining ecological thresholds values for subterranean ecosystems (Di

336

Lorenzo et al., 2015b).

337

Ecotoxicity data reported here for subterranean fauna were obtained from the available

338

literature using: i) Google Scholar, Web of Science and Scopus, with the following strings

339

applied to the topic field: “stygo*”, “troglo*”, “subterranean”, “groundwater”, “thermal

340

tolerance” “sensitivity” and “ecotox*”; and ii) EXPLORE of the US EPA ECOTOX database

341

(https://cfpub.epa.gov/ecotox/explore.cfm). A detailed literature survey was done to review all

342

ecotoxicological data available for subterranean organisms, targeting specifically: 1) organism

343

information (species and ecological classification); 2) organism origin and manipulation

344

before the test start (specimens origin and stage, acclimation to the laboratory condition and

345

test condition); 3) contaminant information (name, class); 3) type of test performed (test

346

media, temperature, aeration, light, feeding and duration) and; 4) evaluation of the estimated

347

effect results (end point and toxicity measured). Data are summarized in Table S.2 for

348

contaminants and Table S.3 for temperature (supplementary material).

349 350

4.1 Contaminants tested on subterranean organisms

351

Subterranean organisms have been tested mostly in acute assays, under conditions of short,

352

often high-concentration exposures and starvation conditions with mortality as end points

353

(Table S.2). Among all studies there was considerable variability in the protocols used. Test

354

organisms were field-collected (rather than laboratory cultures), mostly from Europe, but also

355

from North Africa, western USA and southeastern Australia (Fig. 4). Pre-test acclimation

356

periods ranged from a few hours to eight months and only in some cases the organisms were 15

357

under starvation. The life stage, size or sex of the organism tested was not always recorded

358

and the number of tested organisms per concentration varied between 6 and 30. The medium

359

used in some cases was aerated or renewed and varied from environmental to artificial water,

360

(except for fungi tested in agar) and the temperature ranged from 10.5 to 22 °C. Some studies

361

added food and most were conducted in darkness. The end point measured was mortality

362

(with the exception of growth inhibition in a fungus) expressed most commonly as the lethal

363

(LC50) or effective (EC50) concentrations for half of the population, or as lethal concentration

364

values which affect a specified proportion (e.g. 10, 25, 100%) of the population. The exposure

365

periods (24 h to 28 d) were always shorter than the organism’s life cycle (Table S.2).

366

The majority of toxicity studies with subterranean organisms tested the effects of metals,

367

followed in number by pesticides, fertilizers, chlorine and volatile organic compounds. In all

368

experiments, the contaminants were tested individually (Table S.2).

369

To date most species tested have been groundwater (stygobiont and stygophile) crustaceans:

370

the copepods Bryocamptus echinatus, Diacyclops belgicus, Parastenocaris germanica, two

371

undescribed species of cyclopoids “Cyclopioda_sp1 and Cyclopoida_sp2” and one

372

harpaticoid “Harpaticoida_sp1”; the isopods Caecidotea bicrenata, C. stygia, Proasellus

373

assaforensis, P. cavaticus, P. lusitanicus, and Typhlocirolana haouzensis; the amphipods

374

Metacrangonyx spinicaudatus, Niphargus aquilex and N. rhenorhodanesis; and the decapod

375

Orconectes australis. Non-crustacean tested taxa include the stygobiont Trichodrilus tenuis

376

(annelid worm) and the stygophile fungi Penicillium sp. and Rhodotorula minuta (Table S.2).

377

The sensitivity of any terrestrial subterranean species to toxicants remains unknown.

378 379

4.1.1 Metals

16

380

Metals are the most tested compounds, representing 69% of available data. Acute toxicity has

381

only been tested in subterranean organisms for seven metals and metalloids (As, Cd, Cr, Cu,

382

Ni, Pb and Zn) (Fig. 5). The various studies reported either of two end points: lethal

383

concentration measured at different percentage (e.g., IC25, LC10/LC50) or toxicity data for

384

more than one period (e.g., test duration ranged from 24 h to 28 d) (Table S.2.). The

385

sensitivity responses of the different tested taxa to each metal are variable and no taxon was

386

consistently more or less sensitive, nor was any metal (Fig. 5, Table S.2).

387

The results reflect large variability in the methodological test conditions (see Table S.2),

388

especially regarding the end points and length of the experiment, which confound

389

comparisons of responses between taxa. Bearing this in mind, some patterns in sensitivity are

390

still evident within and between taxonomic groups. For example, the response of stygobiont

391

copepods to Cr and Zn was similar among the three species tested, while stygobiont isopods

392

showed a wide range of sensitivity to Cr, Cu and Zn among the five species tested. There was

393

more than one order of magnitude variation in sensitivity to Cu and Cr at the genus level in

394

Proasellus spp (48-h EC50 = 1.12 mg Cr/L and 6.21 mg Cu/L for P. lusitanicus, and 48-h

395

EC50 = 17.9 mg Cr/L and > 52 mg Cu/L for P. assaforensis), and within the genus

396

Caecidotea, C. bicrenata (96-h LC50 = 2.20 mg Cd/L) was considerably more tolerant to Cd

397

than C. stygia (96-h LC50 = 0.29 mg Cd/L). However, C. bicrenata and C. stygia showed

398

similar sensitivity to Cu (96-h LC50 = 2.2 and 2.3 mg Cu/L, respectively).

399

A wide range of sensitivity was observed across groundwater species, especially for As (from

400

24-h EC50 = 4.4 mg As/L for Rhodotorula minuta to 24-h IC25 = 330 mg As/L for Penicillium

401

sp.). Sensitivities to Cd (from 96-h LC50 = 0.29 mg Cd/L for Caecidotea stygia, to 48-h LC50

402

= 150 mg Cd/L for Typhlocirolana haouzensis), and Zn (from 48-h LC50 = 0.45 mg Zn/L for

403

Metacrangonyx spinicaudatus to 96-h LC50 = 118-180 mg Zn/L for Niphargus aquilex) were 17

404

similarly variable. In contrast, the range of response to Cu and Cr were relatively narrow

405

(from 24-h IC25 = 0.02 mg Cu/L and 0.13 mg Cr/L for Penicillium sp, to 48-h EC50 = 17.9 mg

406

Cr/L and > 52 mg Cu/L for P. assaforensis). The effects of Pb and Ni have been even scarcely

407

tested (Fig. 5, Table S.2).

408

In the absence of data from chronic exposure, some authors expanded the exposure time using

409

the acute test conditions until 28 d, to account for the potential slower response in animals

410

with low metabolic rates (e.g., Canivet et al., 2001; Hose et al., 2016). The major weakness of

411

this approach lies in the fact that extending length of the acute test, will likely affect the

412

results, because it is performed under starvation conditions.

413 414

4.1.2 Pesticides

415

Assays with pesticides represent 16% of all assays using subterranean organisms. Ten

416

different pesticides have been tested (α-endosulfan, 3,4-dichlorophenol, aldicarb, ArianeTM,

417

chlorpyrifos, desethylatrazine, Imazamox, penthaclorophenol, S-metolachlor and thiram) in

418

three species of copepods and one amphipod. The range of the sensitivity observed among the

419

pesticides is broad (e.g. 96-h LC50 = 0.003 mg thiram/L for P. germanica, up to 96-h LC50 =

420

199 mg Imazamox/L for Diacyclops belgicus). Aldicarb was the only pesticide tested in more

421

than one species and the sensitivity was similar for both (B. echinatus 96-h LC50 = 2.71 mg/L,

422

and P. germanica 96-h LC50 = 2.99 mg/L) (Table S.2).

423 424

4.1.3 Fertilizers

425

Ammonium chloride (NH4Cl), ionized ammonia (NH4⁺), and urea (CH4N2O) are the

426

fertilizers that have been tested. Ammonium chloride was more toxic to the amphipod M.

427

spinicaudatus (48-h LC50 = 1.4 mg NH4/L) than to the isopod T. haouzensis (48-h LC50 = 96 18

428

mg NH4/L). Similar LC50 values were reported for the stygophile copepods B. echinatus (96-h

429

LC50 = 14.61 mg/L) and D. belgicus (96-h LC50 = 16.33 mg/L) exposed to ionized ammonia.

430

Di Lorenzo et al. (2014) suggest that the LC50 value for urea (96-h LC50 = 3.14 g/L) was the

431

least toxic of the fertilizers that have been tested (Table S.2).

432 433

4.1.4 Total residuals chlorine

434

Total residual chlorine was tested with two different species. The isopod C. bicrenata (96-h

435

LC50 = 0.11 mg/L) and the decapod Orconectes australis australis (24-h LC50 = 2.7-3.39

436

mg/L) (Table S.2).

437 438

4.1.5 VOCs

439

The toxicity of toluene was studied using the amphipod Niphargus inopinatus. The 96-h LC50

440

values (49-63.9 mg/L) were in the same range of magnitude than the LC50 measured after 13-

441

18 d (39-6-47.8 mg/L) (Table S.2).

442 443

4.1.6 Mixtures

444

Three studies have reported the effect of a toxicant mixture on subterranean organisms.

445

Canivet and Gibert (2002) tested the amphipod N. rhenorhodanensis in a complex mixture of

446

metals using maturated lead secondary smelting slags, where the dilutions were obtained

447

mixing fractions of slags from different percolation stages with river water. Di Marzio et al.

448

(2018) exposed copepods to mixtures of the pesticide Imamox and NH4⁺, and the toxicity of

449

the mixture was synergistic (96-h EC50 = 2.71 mg/L), i.e., it was more toxic than the exposure

450

to the two compounds individually (Imamox 96-h LC50 = 199.23 mg/L; NH4⁺ 96-h LC50

451

16.33 mg/L). No interaction was found between S-metolachlor and desethylatrazine on N. 19

=

452

rhenorhodanensis, with high survival recorded in the highest concentration tested (Maazouzi

453

et al., 2016) (Table S.2).

454 455

4.2 Temperature tested on subterranean organisms

456

Thermal tolerance was studied in terrestrial and aquatic subterranean organisms, however

457

there was considerable variability in the methods used (Table S.3). The end points measured

458

were mainly mortality expressed in upper and lower thermal tolerance, while few included

459

sublethal studies on oxygen consumption as a function of temperature and heat shock proteins

460

synthetization (Table S.3). Overall, the experimental approaches used either fixed

461

temperatures or ramping (increasing or decreasing) temperatures, with between 1 h to three

462

months exposure, or until specimens’ death. Results were expressed in days, temperatures,

463

standard respiration rates or survival percentage for a tested temperature (Table S.3).

464

Four groundwater species have been tested for thermal tolerance. The copepod D. belgicus

465

did not express a change in metabolic activity at temperatures up to 1.1 ºC above its actual

466

maximum environmental temperature (Di Lorenzo and Galassi, 2017). The amphipod N.

467

rhenorhodanensis survived within the range of -2 to 28 ºC. They also accumulated

468

cryoprotective molecules and maintained locomotory activity and aerobic metabolism even at

469

-2 ºC (Issartel et al., 2005). The isopod Proasellus valdensi had wide thermal tolerance, with

470

90% survival between 2 to 16 ºC (Mermillod-Blondin et al., 2013). Only the isopods P. “n.

471

sp. 1” and P. “n. sp. 2” were stenothermic, i.e., surviving within a very narrow thermal niche

472

breadth (Mermillod-Blondin et al., 2013).

473

Thermal tolerance has been tested in 17 terrestrial species (Table S.3). Two beetles from the

474

genus Neobathyscia (N. mancinii and N. pasai) and four beetle imagoes: Trapezodirus

475

arcticollis, Macharoscelis infernus, and two species of Troglocharinus (T. fonti and T. ferreri) 20

476

were tested. The beetles tolerated temperatures from 1 to 20 ºC, but none of the four species

477

survived more than 24 h at 25 ºC. Beetles were exposed to gradually increasing temperatures

478

for less than 24 h, and demonstrated a broad range temperature tolerance. Troglocharinus

479

fonti and T. ferreri had a lower thermal limit (LTL) of -2.5 °C and an upper thermal limit

480

(UTL) of  50.7 °C, while the LTL50 was -11.72°C for N. mancinii and -16.96 °C for N. pasai,

481

and their UTL50 was 28 °C. Sublethal assays reported the synthesis of heat shock proteins

482

(HSPs) by N. mancinii and N. pasai at a temperature close to their UTL50. However, this

483

seems to be unspecific for subterranean species, as the synthesis of these HTPs not always

484

occur in stenothermic organisms, e.g., artic sea organisms (Bernabò et al., 2011). The viability

485

of one lineage of subterranean beetles to cope with climate change was also studied using

486

bioclimatic models. Around 60% of the species were predicted to have suitable conditions

487

(within their physiological tolerance range) under different predicted scenarios to 2080, but

488

most of the species were predicted to be exposed to temperatures and rates of change that they

489

have never experienced through their evolutionary history (Sánchez-Fernández et al., 2016).

490

Thirteen Collembola belonging to three categories of adaptation to subterranean ecosystems

491

(nine troglophiles and four troglobionts) from the western Carpathians, were exposed to

492

gradually increasing and decreasing temperatures (ramping rate of 0.15 ºC/min) before and

493

after being exposed for 1 h to the test temperature (7 to 14 temperatures in a step of 0.6 ºC). In

494

all species, the UTL50 (expressed as heat 50% lethal dose) values were higher than the annual

495

maximum temperature in the respective cave, in which the wider thermal tolerant species was

496

Hypogastura crassaegranulata (from a cold lethal dose, LDc50 = -6.9 °C to a heat lethal dose,

497

LDh50 =  36.6 °C), the subtroglophile Tetrodontophora bielanensis was the most cold-sensitive

498

species (LDc50 = -4.4 °C) and troglobiont Pseudosinella paclti the most heat-sensitive species

21

499

(LDh50 =  31.3 °C), results which outline that cold resistance in Collembola is negatively

500

correlated with species body length. There were no significant differences in upper thermal

501

values between different ecological groups. However, it appears that species belonging to

502

ecological groups less associated with the stable cave environments (trogloxenes,

503

troglophiles) had a wider range of temperature tolerance than the more cave-adapted species

504

(troglophiles or troglobionts), and all troglobionts from deeper caves, were markedly heat-

505

sensitive (Raschmanova et al., 2018).

506 507

5. Discussion

508

Data on the responses of subterranean species to anthropogenic stressors come mostly from

509

measuring the mortality of field-collected organisms. However, experimental conditions

510

varied widely between studies which complicates comparisons of results. The selection of test

511

organisms highlights a strong geographic bias, and the rather ad hoc nature of current research

512

in this field, where the primary criterion for collection site and test taxa was the availability of

513

specimens, the accessibility of sampling points and proximity to the laboratory. The effect of

514

contaminants has not been tested with terrestrial fauna, while the effect of temperature was

515

observed with both terrestrial and aquatic subterranean organisms. The great majority of taxa

516

tested were crustaceans, which corresponds with them being often the most frequently

517

encountered and abundant taxa in groundwater (Sket, 2018). There were more data on the

518

ecotoxicological effects of metals to subterranean organisms than for other toxicants, with test

519

endpoints mostly related to organism sensitivity, however some data on bioaccumulation were

520

also reported for the troglophile millipede Apfelbeckia insculpta (Vranković et al., 2017), the

521

stygobiont amphipods N. rhenorhodanesis (Plénet, 1999; Canivet et al., 2001; Canivet and

522

Gibert, 2002) and N. montelianus (Krupa and Guidolin, 2003), the cave fish Clarias 22

523

gariepinus (Du Preez and Wepener, 2016), and the cave salamander Proteus anguinus

524

(Pezdirc et al., 2011). According to EU (2003), bioaccumulation may be used for secondary

525

poisoning assessment to understand the impacts in the trophic chain, which may be

526

particularly pertinent for subterranean ecosystems, where the trophic chains are simplified.

527

Currently, there is no indication that subterranean organisms are more or less sensitive than

528

surface organisms to pollutants. This is a major knowledge gap that limits environmental risk

529

assessments and the setting of specific thresholds for the protection of subterranean

530

ecosystems. To address this central question, focus should be placed in studying the

531

comparative responses of species of the same genus that have surface closely related species,

532

(e.g., across European asellids), in a similar way that was done for surface vs. subterranean

533

populations of the same species (Jemec et al., 2017).

534

Following the initial research focus on metal toxicity, more recent ecotoxicological

535

approaches with subterranean species have followed general trends of research in surface

536

ecosystems, which have focused increasingly on mixtures, pesticides, VOCs, and temperature

537

increases, where sublethal end points have also been used (Guillén et al., 2012; Gavrilescu et

538

al., 2015). However, the progress of such research in subterranean ecosystems is limited by:

539

a) the lack of knowledge of the concentrations of emerging contaminants and fate of these

540

contaminants in subterranean environments (Stuart et al., 2012), especially in the sediment of

541

both aquifer and terrestrial compartment; b) the limited knowledge of subterranean biota and

542

ecosystems structure; and c) the lack of standardized protocols for subterranean ecotoxicology

543

assays (Gibert and Deharveng, 2002; Danielopol et al., 2004; Hose, 2005; Griebler et al.,

544

2010; Reboleira et al. 2013; Di Lorenzo et al., 2014, 2019). As a consequence, there are

545

relatively few toxicity data available compared to the large amount available for surface

546

ecosystems. 23

547

The three main constraints on the use of subterranean organisms for estimating the effects of

548

contaminants in the subsurface are:

549

i.

No standard criteria for choosing species and selection of model/target species. The

550

well-known “short-range” endemic distribution among all taxonomic groups of

551

subsurface-adapted organisms, leads to an absence of widely distributed species that

552

can be used as model/target species for ecotoxicology studies. Furthermore, the

553

subterranean ecosystem is difficult to access and the populations are usually small

554

(Danielopol et al., 2003; Hose, 2005; Gibert et al., 2009; Daam et al., 2010;

555

Reboleira et al., 2013).

556

ii.

The obstacles to developing stable laboratory cultures with subterranean organisms.

557

The limited knowledge of life history and ecological requirements of subterranean

558

fauna and microbes, and typically low reproductive rates, remain an impediment to

559

their successful culturing in the laboratory. Other traits of subterranean organisms,

560

such as slow metabolism, long life cycles, and low thermal tolerance also complicate

561

the rearing of organisms under laboratory conditions (Thulin and Hahn, 2008;

562

Avramov et al., 2013).

563

iii.

Taxonomic

impediments.

Due

to

morphological

convergence,

identifying

564

subterranean species is challenging and time consuming, requiring proper taxonomic

565

training (Ficetola et al., 2019). The absence of a commercial stock of subterranean

566

species requires the use of field collected animals for assays, in which the species

567

identification based on morphology is not always easy to be achieved (Curini-Galletti

568

et al., 2012; Fonseca et al., 2014; Daam et al., 2010; Hose et al., 2016).

569

Despite these challenges, the management of anthropogenic pressures on subterranean

570

ecosystems requires comprehensive multidisciplinary policy framework that recognizes and 24

571

protects both human and ecological values including ecosystems services (Artigas et al.,

572

2012; Saccò et al., 2019b). Ecotoxicological testing of subterranean organisms remains

573

critical for the protection of subterranean ecosystems, and the standardization of testing

574

protocols, including terrestrial species, is urgently needed in order to progress this research

575

field (Di Lorenzo et al., 2019). Ecotoxicological data are a critical line of evidence in ERA,

576

yet testing approaches lag behind other lines of evidence, such as field biomonitoring data, for

577

which indicators and protocols are established (Danielopol et al., 2004; Danielopol and

578

Griebler, 2008; Griebler et al., 2010; Korbel and Hose, 2011, 2017) or the use of sentinel

579

organisms in situ for real life exposure assessment, which is a valuable first step toward

580

including ecotoxicological data in biological indices (Marmonier et al., 2013, 2018). Even the

581

legislation is ahead of the science; the need for robust scientific data specifically for

582

subterranean ecosystems is clearly articulated in some legislations and environmental

583

agencies. For example, the EU-GWD states that “Research should be conducted in order to

584

provide better criteria for ensuring groundwater ecosystem quality and protection. Where

585

necessary, the findings obtained should be taken into account when implementing or revising

586

this Directive” (EU-GWD, 2006).

587

There is a growing appreciation globally of the need to protect subterranean ecosystems

588

(Mammola et al., 2019) through conservation planning and risk assessments that are specific

589

to subterranean ecosystems (Boulton, 2005, 2009; Linke et al., 2019; Di Lorenzo et al., 2019).

590

Subterranean ecosystems and their fauna need urgent protection from anthropogenic activities

591

because:

592

i.

Subterranean fauna have a high risk of extinction due to fewer possibilities for

593

recolonization and recovery (Culver et al., 2000; Pipan et al., 2010); they are

594

apparently more vulnerable to contaminants due to their narrow geographic ranges, 25

595

high patterns of endemism, and presumably low population size, low densities of

596

species, and low reproductive output (Di Lorenzo et al., 2014).

597

ii.

The low metabolic rates of subterranean organisms relative to phylogenetically

598

related surface species will likely influence toxicant uptake, metabolism and

599

depuration rates, which consequently may lead to different expression of toxicity

600

effects relative to surface species (Hose, 2005, 2007; Avramov et al., 2013; Di

601

Lorenzo et al., 2014, 2015b; Di Marzio et al., 2018). As a consequence, some authors

602

propose increasing the duration of acute toxicity tests with subterranean organisms

603

since a delayed manifestation of toxic effect is expected (Avramov et al., 2013; Hose

604

et al., 2016);

605

iii.

The biotic composition and structure, the abiotic conditions (and their interactions) in

606

subterranean ecosystems differ markedly from those in surface ecosystems.

607

Consequently, it is expected that subterranean biota will also have different

608

sensitivities towards stressors at the community level. The truncated trophic chains

609

and functional redundancy that are common in subterranean communities mean that

610

ecotoxicological testing using realistic approaches, specifically for subterranean

611

ecosystems are needed (Mösslacher and Notenboom, 1999; Daam et al., 2010).

612

There is presently no scientific evidence that subterranean organisms are overall more or less

613

sensitive than surface organisms to pollutants. This paucity of ecotoxicological data for

614

subterranean fauna has led to the use of sensitivity data for surface species as surrogate in the

615

development of environmental quality and protection thresholds. This approach is undesirable

616

because subterranean ecosystems have substantially different environmental conditions, food

617

chains and species traits to surface ecosystems. The development of standardized protocols

618

for testing and, with that, more data specific to subterranean taxa are needed to implement the 26

619

ERA process (Mösslacher and Notenboom, 1999; Hose, 2005; Humphreys, 2007; Daam et al.,

620

2010).

621

The priorities for supporting ERA in subterranean ecosystems should be to: 1) characterize

622

and quantify reliable data on contaminants’ concentrations found in the subterranean

623

ecosystems, including the emerging ones, which are currently out of regular monitoring

624

programs (Artigas et al., 2012); 2) prioritize substances for testing on subterranean species,

625

using existing information and sensitivity data on closely related surface species as surrogates

626

(Artigas et al., 2012, Di Lorenzo et al., 2019); and 3) select model species among the

627

subterranean ones, which should be a representative organism for each trophic level from both

628

compartment terrestrial and aquatic (EC, 2003). The criteria of species selection should be

629

based on knowledge of its morphology, biology. ecological function and routes of exposure

630

(Breithltz et al., 2006). Reporting information about the organism’s feeding behavior is also

631

desirable to enable comparisons between organisms from different trophic levels.

632

It is commonplace in ecotoxicological studies using subterranean species whose results are

633

compared to those for surface species exposed to the same toxicants. Such comparisons

634

provide useful context for the sensitivity of subterranean taxa particularly for comparisons

635

with model test species and those for which standardized protocols are available (e.g.

636

Notenboom et al., 1991; Avramov et al., 2013; Reboleira et al., 2013; Di Lorenzo et al., 2014;

637

Hose et al., 2016) or comparisons with field collected species from the same area of study or

638

other species more closely related (e.g. Bosnak and Morgan, 1981; Boutin et al., 1995; Di

639

Marzio et al., 2009; Di Lorenzo et al., 2014; Maazouzi et al., 2016; Raschmanova et al.,

640

2019). These comparisons have been inconsistent in demonstrating that either group is more

641

or less sensitive than the other. Bosnak and Morgan (1981), Boutin et al. (1995), Avramov et

642

al. (2013), Reboleira et al. (2013) and Maazouzi et al. (2016) observed that surface freshwater 27

643

species were more sensitive than groundwater species to the tested compounds; in contrast,

644

Barr (1976), Bosnak and Morgan (1981) and Di Lorenzo et al. (2014) observed the opposite

645

tendency. Hose et al. (2016) found little congruence between 96 h response data for Daphnia

646

magna and the response of the groundwater species. Interestingly, even subterranean

647

organisms from the same genus differed considerably in their sensitivity to the tested

648

compounds (Reboleira et al., 2013), and it was the species with the highest degree of

649

adaptation to subterranean life (P. lusitanicus) that was more sensitive to pollutants than the

650

species less adapted to subterranean ecosystems (P. assaforensis). For the temperature

651

tolerance, Raschmanova et al. (2019) observed that less adapted species to the subterranean

652

environment showed a wider range of temperature tolerance, but they did not observe

653

significant differences in UTL50 values between different ecological groups while the cold

654

resistance was significantly related to the body length within the same taxon (Collembola).

655

However, comparisons of surface and subterranean species remains confounded by

656

phylogenetic effects and paucity of data such that there are currently insufficient data for a

657

meaningful broad-scale comparison (Hose, 2005).

658 659

6. Future perspectives

660

The environmental and ecological characteristics of subterranean ecosystems demand

661

protocols to assess the impact of contaminants that are specific to the unique traits of the

662

biota. In parallel, a greater understanding of contaminant concentrations and fate in

663

subterranean ecosystems, and deeper knowledge of subterranean ecology are needed for a

664

holistic assessment of risk to subterranean ecosystems. Indeed, these represent the main

665

knowledge needs to progress a framework for environmental assessment of subterranean

666

ecosystems. The key steps to fulfill these knowledge needs are outlined below. 28

i)

667

Greater knowledge of subterranean ecology:

668

Key steps include:

669



Developing standardized methods for sample collection and monitoring to inform risk

670

assessments (Danielopol et al., 2003; Tomlinson et al., 2007; Humphreys, 2007, 2009;

671

Boulton, 2009; Larned, 2012);

672



Understanding the behavior and life cycles of subterranean fauna in order to develop laboratory cultures and enable chronic assays with sublethal end points (Larned, 2012);

673 674



Understanding the structure of subterranean trophic chains to underpin ecotoxicological

675

assays at the community level (e.g. understanding the role of subterranean fauna in the

676

cycling of carbon and nutrients) (Gibert and Deharveng, 2002; Korbel and Hose, 2011);

677



Increasing the knowledge of biodiversity, the relationship between biodiversity and

678

ecosystems functioning (Gibert et al., 2009), and the role of microorganisms and

679

invertebrates in the provision of ecosystem services (Larned, 2012; Hose and Stumpp,

680

2019).

681

ii)

Greater knowledge of contaminant fate and distribution in subterranean ecosystems

682 683

It is important to determine the current extent of contaminants in the subsurface, and to

684

predict the likely presence of contaminants into the future, which requires an understanding of

685

the fate and movement of contaminants in the subsurface (Stuart et al., 2012). Key steps

686

include:

687



Identifying actual and possible new subsurface pollutants, including the development

688

of new analytical techniques and technologies for in situ and real time analysis in

689

groundwater and sediments;

29

690



Increasing monitoring and assessment to determine the environmental occurrence of contaminants;

691 692



Characterizing the sources and identifying pathways of contaminant entry to the

693

subterranean environment, including connectivity with surface ecosystems (e.g.,

694

Bugnot et al., 2019);

695



Defining and quantifying the fate (persistence) and transport (dispersion) of contaminants in subterranean environments.

696 697

iii)

Identify potential ecological effects of contaminant exposure

698

In order to protect subterranean ecosystems and their fauna, it will be necessary to develop

699

new tools and/or adapt those that exist for other ecological domains to assess the chemical

700

and ecological state of subterranean ecosystems, considering also geological heterogeneity

701

and consequent regional endemicity patterns (Reboleira et al. 2013). Only when supported by

702

appropriate tools and data will robust environmental standards and legislation be established

703

to protect subterranean ecosystems. For example, the Water Framework Directive (EU, 2000)

704

and its Groundwater Directive (EU-GWD, 2006) are currently focused on setting the

705

threshold values for all pollutants that put groundwater bodies at risk of achieving a good

706

chemical state. However, to ensure the sustainability of subterranean ecosystems and

707

resources it is necessary to use an integrated approach that includes both the chemical and the

708

ecological status (Stuart et al., 2012). Key steps to achieve this include:

709



Setting appropriate ecotoxicological standards (protocols for culturing and protocols

710

for ecotoxicological assays) (Mösslacher and Notenboom, 1999; Gibert and

711

Deharveng, 2002; Danielopol et al., 2004; Hose, 2005; Humphreys, 2007; Di

712

Lorenzo et al., 2014, 2019), accounting for the specificities of these organisms and

30

713

following recommended guidelines as those for stygobiont crustaceans (Di Lorenzo

714

et al., 2019);

715



Establishing model or target subterranean species in order to generate a database of

716

chemical sensitivity to inform environmental risk assessments (Di Lorenzo et al.,

717

2019);

718



Working with realistic doses or mixtures and combining multilevel endpoints

719

(population, community and ecosystems) to assess the vulnerability of the

720

ecosystems from local to global scales;

721



Implementing an ERA for subterranean ecosystems that is based on relevant

722

scientific evidence and standard procedures for surface ecosystems (EMA, 2006,

723

2018), but modified in terms of requiring ecotoxicological data for at least three

724

different trophic levels to accommodate the simplified trophic chains in subterranean

725

ecosystems.

726

Finally, it is necessary to promote throughout society the importance and urgency of studying

727

and conserving subterranean resources to preserve their ecosystem services. It is also

728

necessary to create a science-policy dialogue to address management questions, ensure the

729

support and uptake of new research, and ultimately ensure the protection and sustainability of

730

subterranean ecosystems (Boulton, 2009).

731 732

Acknowledgements

733

This work was supported by a research grant (15471) from the VILLUM FONDEN.

734 735

Declarations of interest: none.

31

736 737

References

738

Addo-Bediako, A., Chown, S.L., Gaston, K.J., 2000. Thermal tolerance, climatic variability

739

and

740

https://doi.org/10.1098/rspb.2000.1065

741 742

latitude.

Proc.

R.

Soc.

Lond.

B.

Biol.

Sci.

267,

739–745.

Ahmed, F.E., 2001. Analyses of pesticides and their metabolites in foods and drinks. TrACTrend. Anal. Chem. 20, 649–661. https://doi.org/10.1016/S0165-9936(01)00121-2

743

Arias-Estévez, M., López-Periago, E., Martínez-Carballo, E., Simal-Gándara, J., Mejuto, J.C.,

744

García-Río, L., 2008. The mobility and degradation of pesticides in soils and the

745

pollution of groundwater resources. Agric. Ecosyst. Environ. 123, 247–260.

746

https://doi.org/10.1016/j.agee.2007.07.011

747

Arnedo, M.A., Oromí, P., Múrria, C., Macías-Hernández, N., Ribera, C. 2008. The dark side

748

of an island radiation: systematics and evolution of troglobitic spiders of the genus

749

Dysdera Latreille (Araneae: Dysderidae) in the Canary Islands. Invertebrate

750

Systematics, 21(6), 623–660.

751

Artigas, J., Arts, G., Babut, M., Caracciolo, A.B., Charles, S., Chaumot, A., Combourieu, B.,

752

Dahllöf, I., Despréaux, D., Ferrari, B., Friberg, N., Garric, J., Geffard, O., Gourlay-

753

Francé, C., Hein, M., Hjorth, M., Krauss, M., De, H.L., Lahr, J., Lehtonen, K.K.,

754

Lettieri, T., Liess, M., Lofts, S., Mayer, P., Morin, S., Paschke, A., Svendsen, C.,

755

Usseglio-Polatera, P., Van, N. den B., Vindimian, E., Williams, R., 2012. Towards a

756

renewed research agenda in ecotoxicology. Environ. Pollut. 160, 201–206.

757

https://doi.org/10.1016/j.envpol.2011.08.011

758 759

Atkinson, R., 2000. Atmospheric chemistry of VOCs and NOx. Atmos. Environ. 34, 2063– 2101. https://doi.org/10.1016/S1352-2310(99)00460-4 32

760

Avramov, M., Schmidt, S.I., Griebler, C., 2013. A new bioassay for the ecotoxicological

761

testing of VOCs on groundwater invertebrates and the effects of toluene on Niphargus

762

inopinatus.

763

https://doi.org/10.1016/j.aquatox.2012.12.023

Aquat.

Toxicol.

130–131,

1–8.

764

Bernabò, P., Latella, L., Jousson, O., Lencioni, V., 2011. Cold stenothermal cave-dwelling

765

beetles do have an HSP70 heat shock response. J. Therm. Biol. 36(3), 206–208.

766

https://doi.org/10.1016/j.jtherbio.2011.03.002

767 768 769 770

Barceló, D., 2003. Emerging pollutants in water analysis. TrAC-Trend. Anal. Chem. 22, 14– 16. https://doi.org/10.1016/S0165-9936(03)01106-3 Barr, T.C., 1976. Ecological Effects of Water Pollutants in Mammoth Cave: Final technical report to the national park service. University of Kentucky.

771

Bijay-Singh, Yadvinder-Singh, Sekhon, G.S., 1995. Fertilizer-N use efficiency and nitrate

772

pollution of groundwater in developing countries. J. Contam. Hydrol. 20, 167–184.

773

https://doi.org/10.1016/0169-7722(95)00067-4

774

Bosnak, A.D., Morgan, E.L., 1981. Acute toxicity of cadmium, zinc, and total residual

775

chlorine to epigean and hypogean isopods (Asellidae). NSS Bulletin 43, 12–18.

776

Boulton, A.J., 2009. Recent progress in the conservation of groundwaters and their dependent

777 778 779

ecosystems. Aquat. Conserv. 19, 731–735. https://doi.org/10.1002/aqc.1073 Boulton, A.J., 2005. Chances and challenges in the conservation of groundwaters and their dependent ecosystems. Aquat. Conserv. 15, 319–323. https://doi.org/10.1002/aqc.712

780

Boutin, C., Boulanouar, M., Yacoubi-Khebiza, M., 1995. A simple biological test for the

781

lethal toxicity of water and sediments of wells. Compared in vitro toxicity of some

782

heavy metals and ammonium to three genera of aquatic crustacea living in the wells.

783

Hydroécologie appliquée, 7 1–2, 91–109. 33

784

Boyer, D.G., Pasquarell, G.C., 1995. Nitrate Concentrations in Karst Springs in an

785

Extensively Grazed Area1. J. Am. Water. Resour. Assoc. 31, 729–736.

786

https://doi.org/10.1111/j.1752-1688.1995.tb03397.x

787

Brankovits, D., Pohlman, J. W., Niemann, H., Leigh, M. B., Leewis, M. C., Becker, K. W.,

788

Iliffe, T. M., Alvarez, F., Lehmann, M. F., Phillips, B., 2017. Methane-and dissolved

789

organic carbon-fueled microbial loop supports a tropical subterranean estuary

790

ecosystem. Nat. commun, 8(1), 1835. https://doi.org/10.1038/s41467-017-01776-x

791

Breitholtz, M., Rudén, C., Ove Hansson, S., Bengtsson, B.E., 2006. Ten challenges for

792

improved ecotoxicological testing in environmental risk assessment. Ecotox. Environ.

793

Safe. 63, 324–335. https://doi.org/10.1016/j.ecoenv.2005.12.009

794

Bugnot, A.B., Hose, G.C., Walsh, C., Floerl, O., French, K., Dafforn, K.A., Hanford, J.,

795

Lowe, L., Hahs, A., 2019. Urban impacts across realms: Making the case for inter-

796

realm

797

https://doi.org/10.1016/j.scitotenv.2018.08.134

monitoring

and

management.

Sci.

Total

Environ.

648,

711–719.

798

Burkart, M.R., Stoner, J.D., 2008. Chapter 7 - Nitrogen in Groundwater Associated with

799

Agricultural Systems, in: Hatfield, J.L., Follett, R.F. (Eds.), Nitrogen in the

800

Environment,

801

https://doi.org/10.1016/B978-0-12-374347-3.00007-X

802 803

second

ed.

Academic

Press,

San

Diego,

pp.

177–202.

Campos, V., 2002. Arsenic in groundwater affected by phosphate fertilizers at Sao Paulo, Brazil. Env. Geol. 44 (1), 83–87. https://doi.org/10.1007/s00254-002-0540-0

804

Canivet, V., Chambon, P., Gibert, J., 2001. Toxicity and Bioaccumulation of Arsenic and

805

Chromium in Epigean and Hypogean Freshwater Macroinvertebrates. Arch. Environ.

806

Contam. Toxicol. 40, 345–354. https://doi.org/10.1007/s002440010182

34

807

Canivet, V., Gibert, J., 2002. Sensitivity of epigean and hypogean freshwater

808

macroinvertebrates

809

Chemosphere 46, 999–1009. https://doi.org/10.1016/S0045-6535(01)00169-2

810 811 812

to

complex

mixtures.

I:

Laboratory

experiments.

Christiansen, K., 2012. Morphological adaptations, in: Culver, D.C., White, W.B. (Eds.), Encyclopedia of Caves. Elsevier, Amsterdam, pp 517–527. Close, M.E., Skinner, A., 2012. Sixth national survey of pesticides in groundwater in New

813

Zealand.

814

https://doi.org/10.1080/00288330.2012.707131

815

Part

New

Zeal.

J.

Mar.

Fresh

46,

443–457.

Culver, D. C, Deharveng, L., Bedos, A., Lewis, J.J., Madden, M., Reddell, J.R., Sket, B.,

816

Trontelj, P., White, D., 2006. The mid‐latitude biodiversity ridge in terrestrial cave

817

fauna. Ecography. 29(1), 120–128. https://doi.org/10.1111/j.2005.0906-7590.04435.x

818

Culver, D.C., Master, L.L., Christman, M.C., Hobbs, H.H., 2000. Obligate Cave Fauna of the

819

48

820

https://doi.org/10.1046/j.1523-1739.2000.99026.x

821 822 823

Contiguous

United

States.

Conserv.

Biol.

14,

386–401.

Culver, D. C., Pipan, T., 2019. The biology of caves and other subterranean habitats, second ed. Oxford University Press, Oxford. Culver, D.C., Trontelj, P., Zagmajster, M., Pipan, T., 2012. Paving the Way for Standardized

824

and Comparable Subterranean Biodiversity Studies. Subt. Biol. 10, 43–50.

825

https://doi.org/10.3897/subtbiol.10.4759

826

Curini-Galletti, M., Artois, T., Delogu, V., Smet, W.H.D., Fontaneto, D., Jondelius, U., Leasi,

827

F., Martínez, A., Meyer-Wachsmuth, I., Nilsson, K.S., Tongiorgi, P., Worsaae, K.,

828

Todaro, M.A., 2012. Patterns of Diversity in Soft-Bodied Meiofauna: Dispersal

829

Ability

830

https://doi.org/10.1371/journal.pone.0033801

and

Body

Size

Matter.

35

PloS

One

7,

e33801.

831

Daam, M.A., Silva, E., Leitão, S., Trindade, M.J., Cerejeira, M.J., 2010. Does the actual

832

standard of 0.1µg/L overestimate or underestimate the risk of plant protection products

833

to

834

https://doi.org/10.1016/j.ecoenv.2009.12.029

groundwater

ecosystems?

Ecotox.

Environ.

Safe.

73,

750–756.

835

Danielopol, D.L., Gibert, J., Griebler, C., Gunatilaka, A., Hahn, H.J., Messana, G.,

836

Notenboom, J., Sket, B., 2004. Incorporating ecological perspectives in European

837

groundwater

838

https://doi.org/10.1017/S0376892904001444

management

policy.

Environ.

Conserv.

31,

185–189.

839

Danielopol, D.L., Griebler, C., 2008. Changing Paradigms in Groundwater Ecology – from

840

the ‘Living Fossils’’ Tradition to the “New Groundwater Ecology.”’ Internat. Rev.

841

Hydrobiol. 93(4‐5), 565–577. https://doi.org/10.1002/iroh.200711045

842

Danielopol, D.L., Griebler, C., Gunatilaka, A., Notenboom, J., 2003. Present state and future

843

prospects

844

https://doi.org/10.1017/S0376892903000109

845 846

for

groundwater

ecosystems.

Environ.

Conserv.

30,

104–130.

Daughton, C.G., 2005. “Emerging” Chemicals as Pollutants in the Environment: a 21st Century Perspective. Renewable resources journal 6–23.

847

Deharveng, L., Stoch, F., Gibert, J., Bedos, A., Galassi, D., Zagmajster, M., Brancelj, A.,

848

Camacho, A., Fiers, F., Martin, P., Giani, N., Magniez, G., Marmonier, P., 2009.

849

Groundwater

850

https://doi.org/10.1111/j.1365-2427.2008.01972.x

biodiversity

in

Europe.

Freshwater

Biol.

54,

709–726.

851

Di Lorenzo, T., Cifoni, M., Lombardo, P., Fiasca, B., Galassi, D.M.P., 2015b. Ammonium

852

threshold values for groundwater quality in the EU may not protect groundwater

853

fauna: evidence from an alluvial aquifer in Italy. Hydrobiologia 743, 139–150.

854

https://doi.org/10.1007/s10750-014-2018-y 36

855

Di Lorenzo, T., Di Marzio, W.D., Fiasca, B., Galassi, D.M.P., Korbel, K., Iepure, S., Pereira,

856

J., Reboleira, A.S.P.S., Schmidt, S., Hose, G.C., 2019. Recommendations for

857

ecotoxicity testing with stygobiotic species in the framework of groundwater

858

environmental

859

https://doi.org/10.1016/j.scitotenv.2019.05.030

risk

assessment.

Sci.

Total

Environ.

681,

292–304.

860

Di Lorenzo, T., Di Marzio, W.D., Sáenz, M.E., Baratti, M., Dedonno, A.A., Iannucci, A.,

861

Cannicci, S., Messana, G., Galassi, D.M.P., 2014. Sensitivity of hypogean and epigean

862

freshwater copepods to agricultural pollutants. Environ. Sci. Pollut. Res. 21, 4643–

863

4655. https://doi.org/10.1007/s11356-013-2390-6

864

Di Lorenzo, T., Di Marzio, W.D., Spigoli, D., Baratti, M., Messana, G., Cannicci, S., Galassi,

865

D.M.P., 2015a. Metabolic rates of a hypogean and an epigean species of copepod in an

866

alluvial aquifer. Freshwater Biol. 60, 426–435. https://doi.org/10.1111/fwb.12509

867

Di Lorenzo, T., Galassi, D., 2017. Effect of Temperature Rising on the Stygobitic Crustacean

868

Species Diacyclops belgicus: Does Global Warming Affect Groundwater Populations?

869

Water 9, 1–12. https://doi.org/10.3390/w9120951

870

Di Marzio, W.D., Castaldo, D., Pantani, C., Di Cioccio, A., Di Lorenzo, T., Sáenz, M.E.,

871

Galassi, D.M.P., 2009. Relative Sensitivity of Hyporheic Copepods to Chemicals. B.

872

Environ. Contam. Tox. 82, 488–491. https://doi.org/10.1007/s00128-008-9603-6

873

Dijck, S.J.E.V., Laouina, A., Carvalho, A.V., Loos, S., Schipper, A.M., Kwast, H.V. der,

874

Nafaa, R., Antari, M., Rocha, A., Borrego, C., Ritsema, C.J., 2006. Desertification in

875

northern Morocco due to effects of climate change on groundwater recharge, in:

876

Desertification in the Mediterranean Region. A Security Issue, NATO Security

877

Through Science Series. Springer, Dordrecht, pp. 549–577. https://doi.org/10.1007/1-

878

4020-3760-0_26 37

879

Dissanayake, C.B., Chandrajith, R., 2009. Phosphate mineral fertilizers, trace metals and

880

human health. J. Nantn. Sci. Foundation Sri Lanka. 37 (3), 153–165.

881

https://doi.org/10.4038/jnsfsr.v37i3.1219

882

Dole‐Olivier, M.-J., Castellarini, F., Coineau, N., Galassi, D.M.P., Martin, P., Mori, N.,

883

Valdecasas, A., Gibert, J., 2009. Towards an optimal sampling strategy to assess

884

groundwater biodiversity: comparison across six European regions. Freshwater Biol.

885

54, 777–796. https://doi.org/10.1111/j.1365-2427.2008.02133.x

886 887 888

Duarte, A., Cachada, A., Rocha-Santos, T., 2017. Soil Pollution. Elsevier/Academic Press. 312 pp. Duruibe, J.O., Ogwuegby, M.O.C., Egwurugwu, J.N., 2007. Heavy metal pollution and

889

human

biotoxic

effects.

Int.

890

https://doi.org/10.1016/j.proenv.2011.09.146

J.

Phys.

Sci.

2,

112–118.

891

Du Preez, G., Wepener, V., 2016. Influence of mining pollution on metal bioaccumulation

892

and biomarker responses incave dwelling fish, Clarias gariepinus. Bull. Environ.

893

Contam. Toxicol. 97(1), 1823. https://doi.org/10.1007/s00128-016-1799-2

894

EC, 2003. Technical guidance document in support of commission directive 93/67/EEC on

895

risk assessment for new notified substances and commission regulation (EC) No

896

1488/94 on risk assessment for exiting substances. European Commission,

897

Luxembourg.

898 899

EEA, 2019. European Environmental Agency's Waterbase-quality. http://www.eea.europa.eu (accessed 18 June 2019).

900

EEA, 2007. Progress in management of contaminated sites (CSI 015/LSI 003). European

901

environment agency May 2007 assessment. http://www.eea.europa.eu/data-and-

902

maps/indicators 38

903

Ellis, P.A., Rivett, M.O., 2007. Assessing the impact of VOC-contaminated groundwater on

904

surface water at the city scale. J. Contam. Hydrol., Issues in urban hydrology: The

905

emerging

906

https://doi.org/10.1016/j.jconhyd.2006.08.015

field

of

urban

contaminant

hydrology

91,

107–127.

907

EMA (European Medicines Agency), 2018. Guideline on assessing the environmental and

908

human health risks of veterinary medicinal products in groundwater. EMA/

909

CVMP/ERA/103555/2015.

910

EMA (European Medicines Agency), 2006. Guidelines on the environmental risk assessment

911

of medicinal products for human use. Doc. Ref. 627 risks of veterinary medicinal

912

products in groundwater. EMA/CVMP/ERA/103555/2015

913

EU, 2000. Directive 2000/60/EC of the european parliament and of the council of 23 October

914

2000 establishing a framework for Community action in the field of water policy. Off.

915

J. Eur. Union L327, 1–72.EU-GWD, 2006. Directive 2006/118/EC of the European

916

Parliament and of the Council on the protection of groundwater against pollution and

917

deterioration. Off. J. Eur. Union L372, 19–31.

918

Evanko, C.R., Dzombak, D.A., 1997. Remediation of metals-contaminated soils and

919

groundwater (Technology Evaluation Report TE-97-01). Ground-Water Remediation

920

Technologies Analysis Center, Carnegie Mellon University, Deparment of Civil and

921

Environmental Engineering (Pittsburgh, PA).

922

Fan, C., Wang, G.S., Chen, Y.C., Ko, C.H., 2009. Risk assessment of exposure to volatile

923

organic compounds in groundwater in Taiwan. Sci. Total Environ. 407, 2165–2174.

924

https://doi.org/10.1016/j.scitotenv.2008.12.015

39

925

Ficetola, G.F., Canedoli, C., Stoch, F., 2019. The Racovitzan impediment and the hidden

926

biodiversity

of

unexplored

environments.

927

https://doi.org/10.1111/cobi.13179

Conserv.

Biol.

33,

214–216.

928

Fonseca, V.G., Carvalho, G.R., Nichols, B., Quince, C., Johnson, H.F., Neill, S.P.,

929

Lambshead, J.D., Thomas, W.K., Power, D.M., Creer, S., 2014. Metagenetic analysis

930

of patterns of distribution and diversity of marine meiobenthic eukaryotes. Global

931

Ecol. Biogeogr. 23, 1293–1302. https://doi.org/10.1111/geb.12223

932

Francois, C. M., MermillodBlondin, F., Malard, F., Fourel, F., Lécuyer, C., Douady, C. J.,

933

Simon, L., 2016. Trophic ecology of groundwater species reveals specialization in a

934

lowproductivity

935

https://doi.org/10.1111/1365-2435.12484

environment.

Funct.

Ecol.

30(2),

262–273.

936

Galloway, J.N., Townsend, A.R., Erisman, J.W., Bekunda, M., Cai, Z., Freney, J.R.,

937

Martinelly, A., Seitzinger, S.P., Sutton, M.A., 2008. Transformation of the Nitrogen

938

Cycle: Recent Trends, Questions, and Potential Solutions. Science 320, 889–892.

939

https://doi.org/10.1126/science.1136674

940

Gavrilescu, M., Demnerová, K., Aamand, J., Agathos, S., Fava, F., 2015. Emerging pollutants

941

in the environment: present and future challenges in biomonitoring, ecological risks

942

and

943

https://doi.org/10.1016/j.nbt.2014.01.001

bioremediation.

New

Biotechnol.

32,

147–156.

944

Gibert, J., Culver, D.C., DoleOlivier, M.-J., Malard, F., Christman, M.C., Deharveng, L.,

945

2009. Assessing and conserving groundwater biodiversity: synthesis and perspectives.

946

Freshwater Biol. 54, 930–941. https://doi.org/10.1111/j.1365-2427.2009.02201.x

40

947

Gibert, J., Deharveng, L., 2002. Subterranean Ecosystems: A Truncated Functional

948

Biodiversity.

BioScience

52,

473–481.

949

3568(2002)052[0473:SEATFB]2.0.CO;2

https://doi.org/10.1641/0006-

950

Guillén, D., Ginebreda, A., Farré, M., Darbra, R.M., Petrovic, M., Gros, M., Barceló, D.,

951

2012. Prioritization of chemicals in the aquatic environment based on risk assessment:

952

Analytical, modeling and regulatory perspective. Sci. Total Environ. 440, 236–252.

953

https://doi.org/10.1016/j.scitotenv.2012.06.064

954 955 956 957

Griebler, C., Avramov, M., 2015. Groundwater ecosystem services: a review. Freshwater Sci. 34, 355–367. https://doi.org/10.1086/679903 Griebler, C., Lueders, T., 2009. Microbial biodiversity in groundwater ecosystems. Freshwater Biol. 54, 649–677. https://doi.org/10.1111/j.1365-2427.2008.02013.x

958

Griebler, C., Malard, F., Lefébure, T., 2014. Current developments in groundwater ecology—

959

from biodiversity to ecosystem function and services. Curr. Opin. Biotech. 27, 159–

960

167. https://doi.org/10.1016/j.copbio.2014.01.018

961

Griebler, C., Stein, H., Kellermann, C., Berkhoff, S., Brielmann, H., Schmidt, S., Selesi, D.,

962

Steube, C., Fuchs, A., Hahn, H.J., 2010. Ecological assessment of groundwater

963

ecosystems

964

https://doi.org/10.1016/j.ecoleng.2010.01.010

965 966



Vision

or

illusion?

Ecol.

Eng.

36,

1174–1190.

Gruber, N., Galloway, J.N., 2008. An Earth-system perspective of the global nitrogen cycle. Nature 451, 293–296. https://doi.org/10.1038/nature06592

967

Guillén, D., Ginebreda, A., Farré, M., Darbra, R.M., Petrovic, M., Gros, M., Barceló, D.,

968

2012. Prioritization of chemicals in the aquatic environment based on risk assessment:

969

Analytical, modeling and regulatory perspective. Sci. Total Environ. 440, 236–252.

970

https://doi.org/10.1016/j.scitotenv.2012.06.064 41

971

Haarstad, K., Ludvigsen, G.H., 2007. Ten Years of Pesticide Monitoring in Norwegian

972

Ground

Water.

Groundwater

Monit.

973

https://doi.org/10.1111/j.1745-6592.2007.00153.x

Remediat.

27,

75–89.

974

Hancock, P.J., Boulton, A.J., Humphreys, W.F., 2005. Aquifers and hyporheic zones:

975

Towards an ecological understanding of groundwater. Hydrogeol. J. 13, 98–111.

976

https://doi.org/10.1007/s10040-004-0421-6

977

Hashim, M.A., Mukhopadhyay, S., Sahu, J.N., Sengupta, B., 2011. Remediation technologies

978

for heavy metal contaminated groundwater. J. Environ. Manage. 92, 2355–2388.

979

https://doi.org/10.1016/j.jenvman.2011.06.009

980 981

HDR, 2017. Safe Drinking Water Act update, 15th edition. hdrinc.com/sdwa (accessed 18 June 2019).

982

Hose, G.C., 2005. Assessing the Need for Groundwater Quality Guidelines for Pesticides

983

Using the Species Sensitivity Distribution Approach. Hum. Ecol. Risk. Assess. 11,

984

951–966. https://doi.org/10.1080/10807030500257788

985

Hose, G.C., 2007. Response to Comments on Assessing the Need for Groundwater Quality

986

Guidelines for Pesticides Using the Species Sensitivity Distribution Approach. Hum.

987

Ecol. Risk. Assess. 13, 241–246. https://doi.org/10.1080/10807030601105217

988

Hose, G.C., Stumpp, C., 2019. Architects of the underworld: bioturbation by groundwater

989

invertebrates

influences

aquifer

hydraulic

990

https://doi.org/10.1007/s00027-018-0613-0

properties.

Aquat.

Sci.

81.

991

Hose, G.C., Symington, K., Lott, M., Lategan, M., 2016. The toxicity of arsenic(III),

992

chromium(VI) and zinc to groundwater copepods. Environ. Sci. Pollut. R. 23, 18704–

993

18713. https://doi.org/10.1007/s11356-016-7046-x

42

994 995

Humphreys, W.F., 2006. Aquifers: the ultimate groundwater-dependent ecosystems. Aust. J. Bot. 54, 115. https://doi.org/10.1071/BT04151

996

Humphreys, W.F., 2007. Comment on Assessing the Need for Groundwater Quality

997

Guidelines for Pesticides Using the Species Sensitivity Distribution Approach by Hose

998

(2005).

999

https://doi.org/10.1080/10807030601107551

1000 1001

Hum.

Ecol.

Risk.

Assess.

13,

236–240.

Humphreys, W.F., 2009. Hydrogeology and groundwater ecology: Does each inform the other? Hydrogeol. J. 17, 5–21. https://doi.org/10.1007/s10040-008-0349-3

1002

Hutchins, B.T., Schwartz, B.F., Nowlin, W.H., 2014. Morphological and trophic

1003

specialization in a subterranean amphipod assemblage. Freshwater biol. 59(12), 2447–

1004

2461. https://doi.org/10.1111/fwb.12440

1005

Hutton, M., Symon, C., 1986. The quantities of cadmium, lead, mercury and arsenic entering

1006

the U.K. environment from human activities. Sci. Total Environ. 57, 129–150.

1007

https://doi.org/10.1016/0048-9697(86)90018-5

1008

Issartel, J., Hervant, F., Voituron, Y., Renault, D., Vernon, P., 2005. Behavioural, ventilatory

1009

and respiratory responses of epigean and hypogean crustaceans to different

1010

temperatures.

1011

https://doi.org/10.1016/j.cbpb.2005.02.013

Comp.

Biochem.

Phys.

A

141,

1–7.

1012

Jemec, A., Škufca, D., Prevorčnik, S., Fišer, Ž., Zidar, P., 2017. Comparative study of

1013

acetylcholinesterase and glutathione S-transferase activities of closely related cave and

1014

surface Asellus aquaticus (Isopoda: Crustacea). PloS one. 12(5), e0176746.

1015

hpttps://doi.org//10.1271/journal.pone.0176746

1016

Jiménez-Sánchez, M., Stoll, H., Vadillo, I., López-Chicano, M., Domínguez-Cuesta, M.,

1017

Martín-Rosales, W., Meléndez-Asensio, M., 2008. Groundwater contamination in 43

1018

caves:

four

case

of

studies

in

1019

https://doi.org/10.5038/1827-806X.37.1.5

Spain.

Int.

J.

Speleol.

37,

53–66.

1020

JiménezValverde, A., Sendra, A., Garay, P., Reboleira, A.S.P.S., 2017. Energy and

1021

speleogenesis: Key determinants of terrestrial species richness in caves. Ecol. Evol. 7,

1022

10207–10215. https://doi.org/10.1002/ece3.3558

1023

Ju, X.T., Kou, C.L., Zhang, F.S., Christie, P., 2006. Nitrogen balance and groundwater nitrate

1024

contamination: Comparison among three intensive cropping systems on the North

1025

China

1026

https://doi.org/10.1016/j.envpol.2005.11.005

Plain.

Environ.

Pollut.

143,

117–125.

1027

Khan, M.N., Mobin, M., Abbas, Z.K., Alamri, S.A., 2018. Fertilizers and their contaminants

1028

in soils, surface and groundwater. Encyclopedia of the Anthropocene, In: Dominick A.

1029

DellaSala, and Michael I. Goldstein (Eds.) The Encyclopedia of the Anthropocene. 5,

1030

225–240. Oxford: Elsevier. https://doi.org/10.1016/B978-0-12-809665-9.09888-8

1031

Kløve, B., Ala-Aho, P., Bertrand, G., Gurdak, J.J., Kupfersberger, H., Kværner, J., Muotka,

1032

T., Mykrä, H., Preda, E., Rossi, P., Uvo, C.B., Velasco, E., Pulido-Velazquez, M.,

1033

2014. Climate change impacts on groundwater and dependent ecosystems. J.

1034

Hydrology 518, 250–266. https://doi.org/10.1016/j.jhydrol.2013.06.037

1035

Korbel, K., Hose, G.C., 2011. A tiered framework for assessing groundwater ecosystem

1036

health. Hydrobiologia 661, 329–349. https://doi.org/10.1007/s10750-010-0541-z

1037

Korbel, K., Hose G.C., 2015. Habitat, water quality, seasonality or site? Identifying

1038

environmental correlates of the distribution of groundwater biota. Freshwater Sci. 34,

1039

329–343. https://doi.org/10.1086/680038

44

1040

Korbel, K.L., Hose, G.C., 2017. The weighted groundwater health index: Improving the

1041

monitoring and management of groundwater resources. Ecol. Indic. 75, 164–181.

1042

https://doi.org/10.1016/j.ecolind.2016.11.039

1043

Korbel, K., Stephenson, S., Hose, G.C., 2019. Sediment size influences habitat selection and

1044

use

by

groundwater

macrofauna

and

1045

https://doi.org/10.1007/s00027-019-0636-1

meiofauna.

Aquat.

Sci.

81,

39.

1046

Krupa, O.C., Guidolin, L., 2003. Responses of Niphargus montellianus and Gammarus

1047

balcanicus (Crustacea, Amphipoda) from karst waters to heavy metal exposure. J.

1048

Phys. IV 107, 323–326. https://doi.org/10.1051/jp4:20030307

1049

Kulkarni, D., Gergs, A., Hommen, U., Ratte, H.T., Preuss, T.G., 2013. A plea for the use of

1050

copepods in freshwater ecotoxicology. Environ. Sci. Pollut. Res. 20, 75–85.

1051

https://doi.org/10.1007/s11356-012-1117-4

1052

Kumar, P.J.S., Delson, P.D., Babu, P.T., 2012. Appraisal of Heavy Metals in Groundwater in

1053

Chennai City Using a HPI Model. Bull. Environ. Contam. Toxicol. 89, 793–798.

1054

https://doi.org/10.1007/s00128-012-0794-5

1055

Kumaresan, D., Wischer, D., Stephenson, J., Hillebrand-Voiculescu, A., Murrell, J.C., 2014.

1056

Microbiology of Movile Cave a chemolithoautotrophic ecosystem. Geomicrobiol. J.

1057

31(3): 186–193. https://doi.org/10.1080/01490451.2013.839764

1058 1059

Kundzewicz, Z.W., Döll, P., 2009. Will groundwater ease freshwater stress under climate change? Hydrolog. Sci. J. 59, 665–675. https://doi.org/10.1623/hysj.54.4.665

1060

Larned, S.T., 2012. Phreatic groundwater ecosystems: research frontiers for freshwater

1061

ecology. Freshwater Biol. 885–906. https://doi.org/10.1111/j.1365-2427.2012.02769.x

45

1062

Lategan, M.J., Hose, G.C., 2014. Development of a groundwater fungal strain as a tool for

1063

toxicity

assessment.

Environ.

1064

https://doi.org/10.1002/etc.2748

Toxicol.

Chem.

33,

2826–2834.

1065

Lategan, M.J., Klare, W., Kidd, S., Hose, G.C., Nevalainen, H., 2016. The unicellular fungal

1066

tool RhoTox for risk assessments in groundwater systems. Ecotox. Environ. Safe. 132,

1067

18–25. https://doi.org/10.1016/j.ecoenv.2016.05.019

1068 1069

Latella, L., Bernabò, P., Lencioni, V., 2008. Distribution pattern and thermal tolerance in two cave dwelling Leptodirinae (Coleoptera, Cholevidae). Subt. Biol. 6, 81–86.

1070

Linke, S., Turak, E., Asmyhr, M., Hose, G.C., 2019. 3D conservation planning: Including

1071

aquifer protection in freshwater plans alters priorities without much additional effort.

1072

Aquat. Conserv. https://doi.org/10.1002/aqc.3129

1073

Loop, C.M., White, W.B., 2001. A Conceptual Model for DNAPL Transport in Karst Ground

1074

Water

Basins.

Groundwater

1075

6584.2001.tb00357.x

39,

119–127.

https://doi.org/10.1111/j.1745-

1076

Loos, R., Locoro, G., Comero, S., Contini, S., Schwesig, D., Werres, F., Balsaa, P., Gans, O.,

1077

Weiss, S., Blaha, L., Bolchi, M., Gawlik, B.M., 2010. Pan-European survey on the

1078

occurrence of selected polar organic persistent pollutants in ground water. Water Res.

1079

44(14), 4115–4126. https://doi.org/10.1016/j.watres.2010.05.032

1080

Maazouzi, C., Coureau, C., Piscart, C., Saplairoles, M., Baran, N., Marmonier, P., 2016.

1081

Individual and joint toxicity of the herbicide S-metolachlor and a metabolite,

1082

deethylatrazine on aquatic crustaceans: difference between ecological groups.

1083

Chemosphere 165, 118–125. https://doi.org/10.1016/j.chemosphere.2016.09.030

46

1084

Malard, F., Boutin, C., Camacho, A.I., Ferreira, D., Michel, G., Sket, B., Stoch, F., 2009.

1085

Diversity patterns of stygobiotic crustaceans across multiple spatial scales in Europe.

1086

Freshwater Biol. 54, 756–776. https://doi.org/10.1111/j.1365-2427.2009.02180.x

1087

Mammola, S., Cardoso, P., Culver ,D.C., Deharveng, L., Ferreira, R.L., Fišer, C., Galassi,

1088

D.P.M., Griebler, C., Halse, S., Humphreys, W.F., Isaia, M., Malard, F., Martinez, A.,

1089

Moldovan, O.T., Niemiller, M.L., Pavlek, M., Reboleira, A.S.P.S., Souza-Silva, M.,

1090

Teeling, E.C., Wynne, J.J., Zagmajster, M., 2019. Scientists’ warning on the

1091

conservation

1092

https://doi.org/10.1093/biosci/biz064

of

subterranean

ecosystems.

Bioscience

69(8),

641–650.

1093

Mammola, S., Giachino, P.M., Piano, E., Jones, A., Barberis, M., Badino, G., Isaia, M., 2016.

1094

Ecology and sampling techniques of an understudied subterranean habitat: the Milieu

1095

Souterrain Superficiel (MSS). Sci. Nat. 103, 11–12. https://doi.org/10.1007/s00114-

1096

016-1413-9

1097

Marmonier, P., Maazouzi, C., Baran, N., Blanchet, S., Ritter, A., Saplairoles, M., Dole-

1098

Olivier, M.J., Galassi, D.M.P., Eme, D., Dolédec, S., Piscart, C., 2018. Ecology-based

1099

evaluation of groundwater ecosystems under intensive agriculture: A combination of

1100

community analysis and sentinel exposure. Sci. Total Environ. 613–614, 1353–1366.

1101

https://doi.org/10.1016/j.scitotenv.2017.09.191

1102

Marmonier, P., Maazouzi, C., Foulquier, A., Navel, S., François, C., Hervant, F., Mermillod-

1103

Blondin, F., Vieney, A., Barraud, S., Togola, A., Piscart, C., 2013. The use of

1104

crustaceans as sentinel organisms to evaluate groundwater ecological quality. Ecol.

1105

Eng. 57, 118–132. https://doi.org/10.1016/j.ecoleng.2013.04.009

47

1106

Mathews, R., Bosnak, A.D., Tennant, D.S., Morgan, E.L., 1977. Mortality curves of blind

1107

cave crayfish (Orconectes australis australis) exposed to chlorinated stream water.

1108

Hydrobiologia 53, 107–111. https://doi.org/10.1007/BF00029288

1109

Meinel, V.W., Krause, R., 1988. Zur korrelation zwischen Zink und verschiedenen pH-

1110

Werten in ihrer toxischen Wirkung auf einige Grundwasser-Organismen. Zeitschrift

1111

für angewandte Zoologie 75, 159–182.Meinel, V.W., Krause, R., Musko, J., 1989.

1112

Experimente zur pH-Wert Abhängigen Toxizität von Kadmium bei einigen

1113

Grundwasserorganismen. Zeitschrift für angewandte Zoologie 76, 101–125.

1114

Mermillod-Blondin, F., Lefour, C., Lalouette, L., Renault, D., Malard, F., Simon, L., Douady,

1115

C., 2013. Thermal tolerance breadths among groundwater crustaceans living in a

1116

thermally

1117

https://doi.org/10.1242/jeb.081232

1118 1119 1120 1121 1122 1123

constant

environment.

J.

Exp.

Biol.

21(9),

1683–1694.

Momodu, M.A., Anyakora, C.A., 2010. Heavy Metal Contamination of Ground Water: The Surulere Case Study. Res. J. Environ. Earth. Sci. 2, 39–43. Mösslacher, F., Notenboom, J., 1999. Groundwater Biomonitoring. Environmental Science Forum; Trans. Tech. Publications, Switzeland 96, 119–140. National

Water

Quality

Monitoring

Council,

2019.

Water-quality

Portal.

http://www.waterqualitydata.us (accessed 18 June 2019).

1124

Notenboom, J., Cruys, K., Hoekstra, J., Van Beelen, P., 1991. Effect of Ambient Oxygen

1125

Concentration upon the Acute Toxicity of Chlorophenols and Heavy Metals to the

1126

Groundwater Copepod Parastenocaris germanica (Crustacea). Ecotox. Environ. Safe.

1127

131–143. https://doi.org/10.1016/0147-6513(92)90041-Z

1128

Notenboom, J., Boessenkool, J.J., 1992. Acute toxicity testing with the groundwater copepod

1129

Parastenocaris germanica (Crustacea). In proceeding first international conference on 48

1130

ground water ecology, Ameriacan Water Resources Association, Tampa, Florida 301–

1131

309.

1132

Notenboom, J., Plénet, S., Turquin, M.J., 1994. Groundwater contamination and its impact on

1133

groundwater animals and ecosystems, in: Gibert, J., Danielopol, D.L., Stanford, J.A.

1134

(Ed.), Groundwater Eecology. ed. Academic Press, Elsevier.

1135 1136

Nriagu, J.O., 1989. A global assessment of natural sources of atmospheric trace metals. Nature 338, 47–49. https://doi.org/10.1038/338047a0

1137

NRMMC, 2011. Australian Drinking Water Guidelines Paper 6 National Water Quality

1138

Management Strategy. National Health and Medical Research Council, National

1139

Resource Management Ministerial Council, Commonwealth of Australia, Canberra.

1140

https://www.nhmrc.gov.au/sites/default/files/documents/reports/aust-drinking-water-

1141

guidelines.pdf (accessed 18 June 2019).

1142

Pankow, J.F., Thomson, N.R., Johnson, R.L., Baehr, A.L., Zogorski, J.S., 1997. The Urban

1143

Atmosphere as a Non-Point Source for the Transport of MTBE and Other Volatile

1144

Organic Compounds (VOCs) to Shallow Groundwater. Environ. Sci. Technol. 31,

1145

2821–2828. https://doi.org/10.1021/es970040b

1146

Pezdirc, M., Heath, E., Mali, L. B., and Bulog, B., 2011. PCB accumulation and tissue

1147

distribution in cave salamander (Proteus anguinus anguinus, Amphibia, Urodela) in

1148

the polluted karstic hinterland of the Krupa River, Slovenia. Chemosphere 84(7), 987–

1149

993. https://doi.org/10.1016/j.chemosphere.2011.05.026

1150

Pipan, T., Culver, D.C., 2017. The unity and diversity of the subterranean realm with respect

1151

to

invertebrate

body

size.

1152

https://doi.org/10.4311/2016LSC0119

J.

49

Cave

Kars

Stud.

79,

1–9.

1153

Pipan, T., Holt, N., Culver, D.C., 2010. How to protect a diverse, poorly known, inaccessible

1154

fauna: identification and protection of source and sink habitats in the epikarst. Aquat.

1155

Conserv. 20, 748–755. https://doi.org/10.1002/aqc.1148

1156

Plénet, S., 1999. Metal accumulation by an epigean and a hypogean freshwater amphipod:

1157

Consideration for water quality assessment. Water Environ. Res. 71, 1298–1309.

1158

https://doi.org/10.2175/106143096X122401

1159

Raschmanova et al., 2018. Testing the climatic variability hypothesis in edaphic and

1160

subterranean

Collembola

(Hexapoda).

1161

https://doi.org/10.1016/j.jtherbio.2018.11.004

J.

Therm.

Biol.

78,

391–400.

1162

Reboleira, A.S.P.S., Abrantes, N., Oromí, P., Gonçalves, F., 2013. Acute toxicity of copper

1163

sulfate and potassium dichromate on stygobiont Proasellus: general aspects of

1164

groundwater ecotoxicology and future perspectives. Water Air Soil Pollut. 224, 1550.

1165

https://doi.org/10.1007/s11270-013-1550-0

1166

Rivett, M.O., Wealthall, G.P., Dearden, R.A., McAlary, T.A., 2011. Review of unsaturated-

1167

zone transport and attenuation of volatile organic compound (VOC) plumes leached

1168

from

1169

https://doi.org/10.1016/j.jconhyd.2010.12.013

shallow

source

zones.

J.

Contam.

Hydrol.

123,

130–156.

1170

Rizzo, V., Sánchez-Fernández, D., Fresneda, J., Cieslak, A., Ribera, I., 2015. Lack of

1171

evolutionary adjustment to ambient temperature in highly specialized cave beetles.

1172

Evol. Biol. 15, 10. https://doi.org/10.1186/s12862-015-0288-2

1173

Saccò, M., Blyth, A.J., Humphreys, W.F., Kuhl, A., Mazumder, D., Smith, C., Grice, M.,

1174

2019a. Elucidating stygofaunal trophic web interactions via isotopic ecology. PLoS

1175

ONE 14(10), e0223982. https://doi.org/10.1371/journal.pone.022398

50

1176

Saccò, M., Blyth, A., Bateman, P. W., Hua, Q., Mazumder, D., White, N., Humphreys, W.F.,

1177

Laini, A., Griebler, C., Grice, K., 2019b. New light in the dark - a proposed

1178

multidisciplinary framework for studying functional ecology of groundwater fauna.

1179

Sci. Total Environ. 662, 963–977. https://doi.org/10.1016/j.scitotenv.2019.01.296

1180

Sánchez-Fernández, D., Rizzo, V., Cieslak, A., Faille, A., Fresneda, J., Ribera, I., 2016.

1181

Thermal niche estimators and the capability of poor dispersal species to cope with

1182

climate change. Sci. Rep. 6. https://doi.org/10.1038/srep23381

1183

Santos, A., Alonso, E., Callejón, M., Jiménez, J.C., 2002. Distribution of Zn, Cd, Pb and Cu

1184

Metals in Groundwater of the Guadiamar River Basin. Water Air Soil Pollut. 134,

1185

273–283. https://doi.org/10.1023/A:1014149107094

1186 1187

Sarbu, S., Kane, T.C., Kinkle, B.K., 1996. A chemoautotrophically based groundwater ecosystem. Science. 272, 1953–55. https://doi.org/10.1126/science.272.5270.1953

1188

Schmidt, S.I., Hahn, H.J., 2012. What is groundwater and what does this mean to fauna? – An

1189

opinion. Limnologica. 42(1), 1–6. https://doi.org/10.1016/j.limno.2011.08.002

1190

Schneider, K., Culver, D.C., 2004. Estimating subterranean species richness using intensive

1191

sampling and rarefaction curves in a high density cave region in West Virginia. J.

1192

Cave Kars Stud. 66, 39–45.

1193 1194 1195 1196

Sket, B., 2008. Can we agree on an ecological classification of subterranean animals? J. Nat. Hist. 42(21-22), 1549–1563. https://doi.org/10.1080/00222930801995762 Sket, B., 2018. Collecting and processing crustaceans of subterranean habitats. J. Crustacean. Biol. 38, 380–384. https://doi.org/10.1093/jcbiol/rux125

1197

Squillace, P.J., Scott, J.C., Moran, M.J., Nolan, B.T., Kolpin, D.W., 2002. VOCs, Pesticides,

1198

Nitrate, and Their Mixtures in Groundwater Used for Drinking Water in the United

1199

States. Environ. Sci. Technol. 36, 1923–1930. https://doi.org/10.1021/es015591n 51

1200 1201 1202

Sochová, I., Hofman, J., Holoubek, I., 2006. Using nematodes in soil ecotoxicology. Environ. Int. 32(3), 374–383. https://doi.org/10.1016/j.envint.2005.08.031 Steffen, W., Crutzen, P.J., McNeill, J.R., 2007. The Anthropocene: Are Humans Now

1203

Overwhelming

the

Great

Forces

of

Nature.

AMBIO

1204

https://doi.org/10.1579/0044-7447(2007)36[614:TAAHNO]2.0.CO;2

36,

614–621.

1205

Stephenson, S., Chariton, A.A., Holley, M.P., O’Sullivan, M., Gillings, M.R., Hose, G.C.,

1206

2013. Changes in Prokaryote and Eukaryote Assemblages Along a Gradient of

1207

Hydrocarbon Contamination in Groundwater. Geomicrobiol J. 30, 623–634.

1208

https://doi.org/10.1080/01490451.2012.746408

1209

Stuart, M., Lapworth, D., Crane, E., Hart, A., 2012. Review of risk from potential emerging

1210

contaminants

1211

https://doi.org/10.1016/j.scitotenv.2011.11.072

1212 1213

in

UK

groundwater.

Sci.

Total

Environ.

416,

1–21.

Thulin, B., Hahn, H.J., 2008. Ecology and living conditions of groundwater fauna. Swedish Nuclear Fuel and Waste Management Co Technical Report (TR-08-06).

1214

Toccalino, P.L., Gilliom, R.J., Lindsey, B.D., Rupert, M.G., 2014. Pesticides in Groundwater

1215

of the United States: DecadalScale Changes, 1993–2011. Groundwater 52, 112–125.

1216

Tomlinson, M., Boulton, A.J., Hancock, P.J., Cook, P.G., 2007. Deliberate omission or

1217

unfortunate oversight: Should stygofaunal surveys be included in routine groundwater

1218

monitoring programs? Hydrogeol. .J. 15, 1317–1320. https://doi.org/10.1007/s10040-

1219

007-0211-z

1220

Van der Hoff, G.R., van Zoonen, P., 1999. Trace analysis of pesticides by gas

1221

chromatography. J. Chromatogr. A 843, 301–322. https://doi.org/10.1016/S0021-

1222

9673(99)00511-7

52

1223 1224

Van Gestel, C. A., 2012. Soil ecotoxicology: state of the art and future directions. ZooKeys, 176, 275–296. https://doir.org/10.3897/zookeys.176.2275

1225

Van Leeuwen, C.J., Bro-Rasmussen, F., Feijtel, T.C.J., Arndt, R., Bussian, B.M., Calamari,

1226

D., Glynn, P., Grandy, N.J., Hansen, B., Van Hemmen, J.J., Hurst, P., King, N., Koch,

1227

R., Müller, M., Solbé, J.F., Speijers, G.A.B., Vermeire, T., 1996. Risk assessment and

1228

management of new and existing chemicals. Environ. Toxicol. Phar. 2, 243–299.

1229

https://doi.org/10.1016/S1382-6689(96)00072-5

1230

Vranković, J., Borković Mitić, S., Ilić, B., Radulović, M., Milošević, S., Makarov, S., Mitić,

1231

B., 2017. Bioaccumulation of metallic trace elements and antioxidant enzyme

1232

activities in Apfelbeckia insculpta (L. Koch, 1867) (Diplopoda: Callipodida) from the

1233

cave Hadzi-Prodanova

1234

https://doi.org/10.5038/1827-806X.46.1.1981

Pecina

(Serbia).

Int. J.

Speleol.

46

(1),

99–108.

1235

Vesper, D.J., 2019. Contamination of cave waters by heavy metals, in Culver, D.C., White,

1236

W.B., Pipan, T. (Eds.), Encyclopedia of caves. Academic Press, pp. 320–325.

1237

https://doi.org/10.1016/B978-0-12-814124-3.00035-2

1238

Vodela, J.K., Renden, J.A., Lenz, S.D., McElhenney, W.H., Kemppainen, B.W., 1997.

1239

Drinking

1240

trichloroethylene). 1. Interaction of contaminants with nutritional status on general

1241

performance and immune function in broiler chickens. Poult. Sci. 76, 1474–1492.

1242

https://doi.org/10.1093/ps/76.11.1474

1243 1244

water

contaminants

(arsenic,

cadmium,

lead,

benzene,

and

Wakida, F.T., Lerner, D.N., 2005. Non-agricultural sources of groundwater nitrate: a review and case study. Water Res. 39, 3–16. https://doi.org/10.1016/j.watres.2004.07.026

1245

Wang, Q., Li, F., Zhao, L., Zhang, E., Shi, S., Zhao, W., Song, W., Vance, M.M., 2010.

1246

Effects of irrigation and nitrogen application rates on nitrate nitrogen distribution and 53

1247

fertilizer nitrogen loss, wheat yield and nitrogen uptake on a recently reclaimed sandy

1248

farmland. Plant soil. 337(1–2), 325–339. https://doi.org/10.1007/s11104-010-0530-z

1249 1250

White, W.B., Culver, D.C., Pipan, T., 2019. Encyclopedia of Caves. Elsevier/Academic Press, Amsterdam; Boston.

1251

WHO, 2011. Guidelines for drinking-water quality. World Health Organization, Geneva.

1252

Wongsasuluk, P., Chotpantarat, S., Siriwong, W., Robson, M., 2014. Heavy metal

1253

contamination and human health risk assessment in drinking water from shallow

1254

groundwater wells in an agricultural area in Ubon Ratchathani province, Thailand.

1255

Environ. Geochem. Health. 36, 169–182. https://doi.org/10.1007/s10653-013-9537-8

1256 1257

World Water Assessment Programme, 2003. “Water for people, water for life.” The United Nations World Water Development Report. UNESCO-WWP, Paris.

1258

54

1

Figure captions

2

Figure 1. Conceptual model of surface vs. subterranean ecosystems, over time and depth,

3

pointing out the simplified food web, generally lower organic matter content, higher relative

4

humidity and small temperature variation in the subterranean ecosystems.

5

Figure 2. Anthropogenic stressors in the subterranean ecosystems, with contamination sources

6

and threatened ecosystem services.

7

Figure 3. Maximum Concentration Detected (MCD) of contaminants in European

8

groundwater, data from the European Environmental Agency's Waterbase-quality (EEA,

9

2019). A) Metals; B) Pesticides; C) Fertilizers; D) Volatile Organic Compounds; and E)

10

Emerging contaminants. Asterisks points out the compounds that have been tested in

11

subterranean organisms.

12

Figure 4. World distribution of the subterranean species tested for stressors. Circles represent

13

species tested for contaminants sensitivity and squares represent species tested for upper-limit

14

survival temperature.

15

Figure 5. Graphical summary of the ecotoxicological endpoints for metal exposure on

16

groundwater species. End points shown: LC50 (Lethal Concentration, 50%) and IC25

17

(Inhibition Concentration, 25%). * no toxic responses for the maximum tested concentrations,

18

i.e., LC50 value is higher.

1

Highlights: •

Impacts of anthropogenic stressors on subterranean ecosystems are poorly known



Interdisciplinary research is needed for subterranean ecosystems’ sustainability



Subterranean biota is neglected in legal protection frameworks



Contaminants and temperature variation have pernicious effects on subterranean biota



Standardization of reliable testing protocols is urgently needed

Declaration of interests ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. ☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: