Discharge Driven Nitrogen Dynamics in a Mesoscale River Basin As

Jul 22, 2016 - The purpose of simulating the discharge for all investigated subcatchments of the Bode River basin is to provide information on nitrate...
0 downloads 11 Views 2MB Size
Subscriber access provided by STEPHEN F AUSTIN STATE UNIV

Article

Discharge driven nitrogen dynamics in a mesoscale river basin as constrained by stable isotope patterns Christin Mueller, Matthias Zink, Luis Samaniego, Ronald Krieg, Ralf Merz, Michael Rode, and Kay Knoeller Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b01057 • Publication Date (Web): 22 Jul 2016 Downloaded from http://pubs.acs.org on July 25, 2016

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

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

Page 1 of 29

Environmental Science & Technology

1

Discharge driven nitrogen dynamics in a mesoscale

2

river basin as constrained by stable isotope patterns

3

Christin Mueller1*, Matthias Zink2, Luis Samaniego2, Ronald Krieg1, Ralf Merz1, Michael Rode3,

4

Kay Knöller1

5

[1] Helmholtz Center for Environmental Research – Department Catchment Hydrology - Stable

6

Isotope Group, Theodor-Lieser-Straße 4, D-06120 Halle (Saale), Germany

7

[2] Helmholtz Center for Environmental Research – Department Computational Hydrosystems -

8

Permoserstraße 15, 04318 Leipzig, Germany

9

[3] Helmholtz Center for Environmental Research – Department Aquatic Ecosystem Analysis -

10

Brückstraße 3a, 39114 Magdeburg, Germany

11 12 13 14 15 16 17

*Corresponding Author. Christin Mueller Address: Helmholtz Centre for Environmental Research, UFZ Theodor-Lieser-Str. 4 06120 Halle (Saale) – GERMANY Phone: +49 345 558 5490 E-Mail: [email protected]

18 19

KEYWORDS: stable isotope pattern, regional process dynamics, mesoscale hydrological model (mhM), TERENO

20

ACS Paragon Plus Environment

Environmental Science & Technology

21

ABSTRACT: Nitrate loads and corresponding dual-isotope signatures were used to evaluate large

22

scale N dynamics and trends in a river catchment with a strong anthropogenic gradient (forest

23

conservation areas in mountain regions, and intensive agriculturally used lowlands). The Bode

24

River catchment with an area of 3200 km2 in the Harz Mountains and central German lowlands

25

was investigated by a two years monitoring program including 133 water sampling points each

26

representing a sub-catchment. Based on discharge data either observed or simulated by the

27

mesoscale Hydrological Model (mhM) a load based interpretation of hydrochemical and isotope

28

data was conducted. Nitrate isotopic signatures in the entire catchment are influenced by (I) the

29

contribution of different nitrogen sources, (II) by variable environmental conditions during the

30

formation of nitrate, and (III) by a minor impact of denitrification. For major tributaries, a

31

relationship between discharge and nitrate isotopic signatures is observed. This may in part be

32

due to the fact, that during periods of higher hydrologic activity a higher wash out of isotopically

33

lighter nitrate formed by bacterial nitrification processes of reduced or organic soil nitrogen

34

occurs. Beyond that, in-stream denitrification seems to be more intense during periods of low

35

flow.

36

INDRODUCTION

37

Nitrogen cycling in river systems impacted by human activity. Nitrogen is an essential

38

nutrient for the synthesis of nucleic acids and proteins, for the growth of biomass and thus for

39

agriculture production. On the other hand, it is a major multiplier on eutrophication of rivers,

40

lakes and coastal ecosystems. The emission of nitrate has increased in the last decade triggered

41

by three main causes: (1) the application of organic or inorganic fertilizer; (2) the widespread

42

cultivation of rice, legumes, and other crops that are capable of converting atmospheric N2 to

43

organic nitrogen and (3) the combustion of fossil fuels which transforms gaseous N2 and fossil N

ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

Environmental Science & Technology

44

to reactive nitrogen oxides1-3. To control the impact of nitrate and its sustainable mitigation in

45

aquatic systems, it is crucial to understand and quantify sources as well as biochemical processes

46

which (permanently) remove nitrate4-8. Isotope signatures of nitrate can help identify and quantify

47

potential nitrate sources and characterize biogeochemical N-transformation processes9-14.

48

Detailed information about relevant biochemical processes of natural and anthropogenic nitrate

49

sources in the nitrogen cycle can be found in the support information part.

50

Microbiological isotope fractionation in N-compounds. The two major nitrate consuming

51

bacterial processes in aquatic systems, NO3- denitrification and assimilation, result in significant

52

isotopic fractionation rates (Table S1, supporting information). Nitrate which contains light

53

isotopes (14N and

16

54

isotopes (15N and

18

55

processes, which is measured in the streams, is generally more enriched in heavy isotopes.

57

Impact of environmental conditions on the N-cycle in a river system. Nitrate process

58

dynamics in river catchments are influenced by climate related variables and their changes,

59

especially increasing stream temperature and discharge variability due to changed precipitation

60

patterns15. Similarly, denitrification rates increase with higher water temperatures2,7. Boyacioglu

61

et al.15 found that especially in wet periods and with the canalization of streams, denitrification

62

rates decrease. However, discharge variability has a higher impact on denitrification than

63

seasonal temperature changes15. In contrast, Wexler et al.16 observed microbial denitrification in

64

temperate forest catchments in groundwater only, whereby surface water showed no evidence for

65

the occurrence of bacterial denitrification. They also found that during small rain events

66

precipitation mobilizes nitrate that makes up about 34 % of NO3- in the stream water. This nitrate

67

is obviously derived from the microbial nitrification of soil-nitrogen. Generally, N-inputs by

O) is generally more rapidly consumed by bacteria than nitrate with heavy O). Thereby, the residual nitrate pool during microbiological fractionation

ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 29

68

atmospheric deposition exceed the export from catchments. This implies an accumulation in the

69

soil zone or a loss as gaseous N2 during denitrification17, 18.

70

OBJECTIVES

71

Past research related to the quantitative assessment of the nitrogen balance in catchments is

72

mainly focused on small catchments6,16 or does not include the advantage of stable isotope

73

observations19. In this study, we clarify potentially occurring nitrogen transformation processes

74

on a meso- to large scale river catchment considering hydrochemical and stable isotope

75

observations. The purpose of this study is to examine spatio and temporal variabilities of nitrate

76

stream loading related to land use and human impact on the N-cycle. This includes the

77

identification of different nitrogen sources and microbiological turnover processes that consume

78

or produce nitrate (Table S1 in supporting information) using

79

combined with water chemical parameters and hydrological modeling. Furthermore, we intended

80

to recognize and evaluate predominate catchment transformation zones on a large scale and the

81

impact of different hydrological conditions (base flow and high flow) on the catchment scale

82

nitrogen dynamics. Also, we want to investigate to what catchment size (stream order)

83

information on N-sources and N-turnover processes can be reliably differentiated by isotope

84

methods. With our study, we want to contribute a tool set that helps predict the nitrate export

85

from large river systems as a prerequisite for suggesting future safety precautions that guarantee a

86

good ecological status of the water bodies from headwater to coast.

87

STUDY AREA

88

The Bode catchment is a hydrologically and hydrochemically intensively studied region in

89

central Germany (Figure 1). It is part of the Terrestrial Environmental Observatories (TERENO)

90

network which is hosted within the Helmholtz Association of German Research Centres20. The

15

ACS Paragon Plus Environment

N and

18

O nitrate isotope data

Page 5 of 29

Environmental Science & Technology

91

Bode catchment is characterized by a spatial gradient of diverse nitrate sources. The southern

92

natural protected headwaters in the Harz Mountains represent the low nitrate concentrated region.

93

With increasing river length the impact of agricultural land use increases. Major tributaries can be

94

differentiated into four main groups after Mueller et al.21:

95



Firstly, the mountainous river headwaters with a discharge of about 28 % of total flow

96

discharging into the Rappbode drinking water reservoir including the Kalte Bode, Warme

97

Bode, Rappbode and Luppbode stream.

98



Secondly, transfer streams (Hassel, Selke and Holtemme with 27 % of the total

99

discharge) flow between protected areas dominated by forest and areas of agricultural

100

land use. Within the group of transfer streams, 21 sub-catchments close to the headwater

101

catchments are solely covered by forest. Further downstream, 9 sub-catchments are

102

characterized by an increasing impact of agricultural land use and a mixed forest and

103

agricultural zone (forest and agricultural land use covering each 50 % ± maximal 10 %).

104

The tributaries are dominated mainly by agricultural land use (>60 % of each sub-

105

catchment is covered by arable land) which is the case for 13 sub-catchments.

106



with 21 % of total discharge).

107 108 109

Thirdly, the northern agricultural lowlands are dominated by arable land (Großer Graben



Fourthly, other minor tributaries which are not considered in the sampling network or diffuse inflow (including groundwater) account for 22 % of the total discharge.

110

Further characterization of major streams, climatic conditions as well as landscape information

111

can be found in Mueller et al.21.

112

The spring water from the Harz Mountains is used for drinking water and electrical power supply

113

which is collected in the Rappbode reservoir. The LHW (Landesbetrieb für Hochwasserschutz ACS Paragon Plus Environment

Environmental Science & Technology

114

und Wasserwirtschaft) Saxony-Anhalt has 30 observation gauges which continuously register

115

water-levels along the Bode River and major tributaries (Figure 1).

116

METHODS

117

Monitoring approach. In 2012, we started our two-year monitoring campaign in the Bode River

118

catchment. It includes 133 sampling sites (Figure 1) which represent nearly undisturbed

119

mountainous regions, forests, settlement areas as well as anthropogenic impacted areas by waste

120

water and agricultural land use. Further information about regional, climatic and hydrological

121

aspects are given in Zacharias et al.20 and Mueller et al.21. Spatial high resolution sampling

122

campaigns were performed three times per year. Low resolution sampling campaigns including

123

25 representative test sites were sampled monthly (starting in 2013), between the high resolution

124

campaigns. All main tributaries flowing into the Bode River were included into the monitoring

125

network. Precipitation was collected monthly at the 25 representative test sites. During high and

126

low resolution sampling campaigns, water samples were collected at the outlets of the major sub-

127

catchments. Water samples were filtered in the field (0.45 µm cellulose acetate filter) and stored

128

in high density polyethylene (HDPE)-bottles in the refrigerator. Additionally, all water samples

129

were measured for pH-value, temperature and electrical conductivity in the field. The density of

130

sampling spots is much higher in the southern mountainous area due to the high number of small

131

creeks compared to the northern lowland area.

132

Laboratory analysis: nitrate and isotopic signatures. Laboratory analyses on all water samples

133

were conducted for stable isotope signatures of nitrate (δ15N and δ18O) and corresponding NO3

134

contents within one month after sampling. Concentration of nitrate in water was analyzed by ion

135

chromatography using a Dionex ICS-2000 combined with AS50. Nitrate isotopic signatures were

136

measured with the denitrifier method by using bacteria strains of Pseudomonas chlororaphis

ACS Paragon Plus Environment

Page 6 of 29

Page 7 of 29

Environmental Science & Technology

137

(ATCC #13985)23,24. Isotopic ratios are expressed in delta notation (equation 1). δ15NNO3 is

138

reported relative to atmospheric nitrogen (AIR), whereby δ18ONO3 is expressed relative to Vienna

139

Standard Mean Ocean Water (VSMOW), calculated as:

140

𝑅𝑠𝑎𝑚𝑝𝑙𝑒

𝛿𝑠𝑎𝑚𝑝𝑙𝑒 [‰] = (𝑅

𝑠𝑡𝑎𝑛𝑑𝑎𝑟𝑑

− 1) × 1000

[1]

141

Ratio R describes the relation between heavier and lighter isotopes. For the analytical

142

measurement a mass DELTA V Plus spectrometer was used in combination with a GasBench II

143

from Thermo Scientific. The standard deviation of the described analytical measurement

144

including sampling and preparation uncertainties is ±1.6 ‰ for δ18O and ±0.4 ‰ for δ15N.

145

Isotope measurements are conducted in duplicates and are repeated when the standard deviation

146

exceeds the quality control threshold stated above. For calibration of nitrogen and oxygen isotope

147

values, the reference nitrates IAEA-N3 (δ15N: +4.7 ‰ AIR; δ18O: +25.6 ‰ VSMOW), USGS-32

148

(δ15N: +180 ‰ AIR; δ18O: +25.7 ‰ VSMOW), USGS-34 (δ15N: -1.8 ‰ AIR; δ18O: -27.9 ‰

149

VSMOW), and USGS-35 (δ15N: +2.7 ‰ AIR; δ18O: +57.5 ‰ VSMOW) were used.

150

Regionalization of runoff: Hydrological modeling and nitrogen loads. The purpose of

151

simulating the discharge for all investigated sub-catchments of the Bode River basin is to provide

152

information on nitrate export from the respective sub-catchments and to understand the process

153

dynamics across scales. Large stream basins are the sum of its associated sub-catchments and can

154

be regionalized with different rainfall-runoff models, as explained and compared in e.g. Oudin et

155

al.25. The mesoscale Hydrologic Model mHM used in this study is based on accepted

156

hydrological conceptualizations and employs a multi-scale parameter regionalization (available at

157

www.ufz.de/mHM). It is predicated on numerical approximations of commanding hydrological

158

processes. Detailed information can be found in Samaniego et al.26. The software mHM is a

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 29

159

distributed hydrological model which treats grid cells as hydrological units. It is a conceptual

160

hydrological model similar to HBV27 or VIC28, which considers interception, snow accumulation

161

and melting, infiltration, soil water retention, recharge, groundwater storage and runoff

162

generation as main hydrological processes. The generated discharge consists of direct runoff,

163

slow and fast interflow and baseflow. It is routed through the study domain using the Muskingum

164

method after Tewolde et al.29.

165

A key feature of mHM is the estimation of effective parameters called Multiscale Parameter

166

Regionalization technique (MPR26). Effective model parameters are estimated using a two-step

167

approach. In a first step the parameters (e.g. porosity) are estimated on the resolution of the

168

morphological input data (i.e. 100 x 100 m). With the help of transfer functions (e.g. pedotransfer

169

function) and a couple of transfer parameters (e.g. coefficients of the pedotransfer function),

170

which are associated with terrain characteristics (e.g. sand and clay content), spatially distributed

171

parameter fields are estimated. In a second step, this parameter fields are upscaled to the

172

resolution of the hydrological simulation (i.e. 1 x 1 km2). Thus, mHM explicitly accounts for sub

173

grid variabilities and depends only on a couple of transfer or global parameters, which are space

174

and time invariant. They are estimated during an automated calibration. MPR enables mHM to

175

perform satisfying across different locations and scales26.

176

The objective of the calibration is to minimize the Nash Sutcliffe Efficiency (NSE)30 of discharge

177

and the logarithm of discharge using the Dynamical Dimensioned Search Algorithm31. The NSE

178

(equation 2) is calculated as

179

𝐸 = 1.0 −

∑𝑁 𝑖=1(𝑄𝑜𝑏𝑠 (𝑖)− 𝑄𝑠𝑖𝑚 (𝑖))² ∑𝑁 𝑄𝑜𝑏𝑠 (𝑖)− 𝑄̅𝑜𝑏𝑠 )² 𝑖=1

ACS Paragon Plus Environment

[2]

Page 9 of 29

Environmental Science & Technology

180

̅ obs the mean where Qobs denotes the observed discharge, Qsim the simulated discharge and Q

181

observed discharge for all observations N. The calibration is done on the spatial resolution of

182

2 km in the time period 2000-2009, whereas discharge simulations for analyzing nitrate loads are

183

conducted at 1 km spatial resolution. The calibration of mHM at a coarser resolution is reasoned

184

in the savings in computational time. Kumar et al.32 showed that this is a reasonable choice.

185

Since the Bode catchment is highly influenced by human interaction, the model has been fed by

186

measured discharge time series downstream the area of the biggest reservoir, the Rappbode dam.

187

Thus, the simulated upstream discharge time series has been neglected at this point. The feeding

188

or inflow gauge is located at gauge B1 shown in Figure 1. For calibration the gauges B3 and S4

189

are used (Figure 1). B3 is the gauge which is representing the headwaters of the Bode River and

190

advantageously, is still not influenced by the hydraulic connection to the neighboring catchment.

191

S4 is the outlet of the Selke River, which is the most natural behaving sub-catchment in the Bode

192

basin. Furthermore, it is used for detailed nitrate dynamics analysis within this study. The

193

estimated discharge at 1 x 1 km2 is evaluated at 20 out of 30 gauges, which were not or only

194

insignificantly influenced by human activity. All simulations are conducted with a spin up period

195

of 10 years to ensure that the initial conditions are reasonable.

196

RESULTS AND DISCUSSION

197

Discharge simulation. The discharge simulations conducted with the hydrological model mHM

198

could not be evaluated at all gauges available in the Bode catchment, because gauges downstream

199

of the Rappbode dam are influenced by upstream water regulations or uptake. The results of the

200

majority of the remaining gauges indicate a sufficient estimation of discharge with an average

201

NSE of 0.65 for all gauges in the validation period 1980-1999 (Table 1) according to Moriasi et

202

al.33. The gauge GG1 reveals low performances because a hydraulic connection from the Großer

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 29

203

Graben to the neighboring catchment does exist. This connection leads to the exchange of water

204

between both catchments. At the river Holtemme and its tributaries, many water regulations are

205

located in the upstream areas. In consequence, the downstream gauges are impacted by these

206

regulations. This is the reason for lower NSEs at the gauges HO1, HO3 and HO10. One of the

207

largest summer floods during decades took place in June 2013 in the Mulde, Saale and Elbe

208

catchment, of which the Bode is a tributary34. This extreme event is captured by the model

209

simulation, but the peak discharge is usually underestimated. One of the reasons is the magnitude

210

of the event, which has rarely been observed before. This induces errors in the observed

211

discharge time series, due to considerable uncertainties in the rating curve. Another reason could

212

be erroneous precipitation data, caused by the high wind speeds during the storm event occurring

213

in connection with the flood event and insufficient spatial density of rainfall gauges for that

214

event. The underestimation of the flood event in May and June of 2013 can be observed in

215

Figure 2 for the Selke tributary and the entire Bode catchment. The underestimation of discharge

216

at the Bode outlet (Figure 2) during this period is additionally influenced by supplementary water

217

entering the Bode catchment from the neighboring catchment via the hydraulic connection at the

218

tributary Großer Graben. Finally, the model simulations are sufficient for relating loads observed

219

at the sampling points with simulated discharge.

ACS Paragon Plus Environment

Page 11 of 29

Environmental Science & Technology

220

Table 1. Observed and simulated discharge at gauge stations (o - outlet, MQ - mean discharge,

221

NSE - Nash-Sutcliffe efficiency)

Landscape type

Intensive agriculture

Intensive agriculture and forest

Headwater

ID

222 223

*2

NSE

Time period

MQ [m³ s-1]

2000-2009

1980-1999

2012-2013

5.48

0.90

0.91

0.73

B2

Bode

1961-2013

B3

Bode

1893-2013

7.80

0.89

0.85

0.63

B4

Bode

1930-2013

12.97

0.81

0.79

0.68

B5 GG1

Bode (o) Großer Graben (o)

1988-2013

12.64

0.71

0.75

0.67

1986-2013

2.41

0.44*1

0.51*1

0.57*1

HO1

Holtemme

1961-2013

0.28

0.38

0.40

0.38

HO2

Holtemme

1970-2013

0.34

0.67

0.60

0.58

HO3

Holtemme

1983-2013

0.13

0.59

0.40

0.45

Holtemme

1971-2013

1.34

0.65

0.72

0.59

HO8

0.47

*2

0.49

*2

0.56*2

HO10

Holtemme (o)

1982-2013

1.55

S1

Selke

2006-2010

0.33

0.72

-

-

S2

Selke

1948-2013

1.13

0.64

0.68

0.56

S3

Selke

1920-2013

1.46

0.70

0.69

0.45

S4

Selke (o)

1980-2013

1.73

0.79

0.74

0.54

KB1

Kalte Bode (o)

1950-2013

0.72

0.57

0.52

0.42

WB1

Warme Bode

1960-2013

2.06

0.81

0.69

0.75

WB2

Warme Bode (o)

1950-2013

2.42

0.64

0.67

0.75

HA1

Hassel (o)

1968-2013

0.54

0.80

0.52

0.69

RB1

Rappbode (o)

1950-2013

0.75

0.73

0.71

0.76

Luppbode (o)

1996-2013

0.35

0.60

0.80

0.47

LB1 *1

River

hydraulic connection to the neighboring catchment impact of regulated tributary (Zillierbach dam)

224

Distribution of discharge and nitrate load on catchment scale. The average annual discharge

225

intensity on the Bode River outlet (B5) is about 8.3 m3s-1 in 2012 and 14.9 m3s-1 in 2013 with a

226

mean nitrate load of about 6.5 t per year. The distribution of mean discharge of major tributaries

227

of the Bode catchment over the years 2012 and 2013 is shown in Figure 3A. Discharge at the

228

Bode outlet is composed of four components of comparable size: (I) headwater tributaries, (II)

229

tributaries in mixed land use catchment with intensive agriculture and forest (transfer streams),

230

(III) tributaries in catchments with agricultural land use only and (IV) minor tributaries directly to

231

the Bode River or other unspecified inflows such as groundwater.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 29

232

The corresponding calculated nitrate loads are distributed differently as shown in Figure 3B. The

233

lowland agricultural tributaries (Großer Graben) and the transfer streams (Holtemme and Selke

234

River) feature the major nitrate impacts with 73 % of the total nitrate load. Tributaries upstream

235

the Rappbode dam in the high mountains transport low amounts of nitrate (6 % of the total nitrate

236

load). Most of the samples were taken in the agricultural regions with 56 samples representing

237

the major land cover type covering 75 % of the entire Bode catchment area. The forest region and

238

headwaters are characterized by samples from 37 sub-catchments and 24 sub-catchments,

239

respectively. The remaining 16 sub-catchments show a mixed land cover type of intensive

240

agriculture and forest (each 50 % ± maximal 10 %).

242

243

For further temporal and spatial process dynamics analysis on catchment scale, the discharge and

244

corresponding nitrate loads are weighted to the respective sub-catchment size. The area-weighted

245

or so-called specific discharge (Qspecific) and nitrate loads ((NO3-)specific, are calculated using

246

equation (3) and (4)), respectively.

247

248

𝑄𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐 = 𝐴

𝑄𝑠𝑖𝑚

𝑠𝑢𝑏 𝑏𝑎𝑠𝑖𝑛

(𝑁𝑂3− )𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐 =

[mm a-1]

𝑁𝑂3− 𝑐𝑜𝑛𝑐. ∗ 𝑄𝑠𝑖𝑚 𝐴𝑠𝑢𝑏 𝑏𝑎𝑠𝑖𝑛

[3]

[g m-2 a-1]

[4]

249

Temporal and spatial dependency of specific nitrate loads and isotope signatures on specific

250

discharge. To detect possible nitrate sources and sinks, the comparison of specific discharge and

251

nitrate loads between major land cover types along the river system can provide useful

252

information21. In general, the specific nitrate load increases with increasing specific discharge.

253

For agricultural catchments specific nitrate load is much higher than for non-agricultural

ACS Paragon Plus Environment

Page 13 of 29

Environmental Science & Technology

254

catchments for the whole range of specific discharges (Figure 4A). The specific nitrate loads in

255

agricultural areas are significantly higher compared to non-farming regions with a maximum

256

specific load of 70.3 g m-2 a-1 and a corresponding specific runoff of 2.2 mm a-1. The specific

257

nitrate loads for headwater, forest and mixed regions do not differ significantly between each

258

other. The strongest correlation between specific discharge and specific nitrate load is obvious for

259

agricultural catchments. For all other land use types, correlations are much weaker. This implies

260

that the nature of the nitrate source is different between agricultural and all other land use types.

261

The temporal variations of nitrate concentrations in (mainly lowland) tributaries are higher than

262

in the main stream itself. The highest coefficients of variation cv related to temporal nitrate load

263

variability of up to 40 % are obvious for agricultural tributaries in the northern sub-basins

264

(Figure 4B). A significant role plays the temporal differentiated release of fertilizer in the

265

agricultural areas as well as fluctuating contribution of urban waste water.

266

The relationship between mean isotopic signatures of nitrate (δ15N/δ18O) for each sub-catchment

267

and the proportion of agricultural land use is shown in Figure 5A and B. The nitrogen isotopic

268

signature generally increases with increasing agricultural portions (Figure 5A). The outlet of the

269

main river indicates a δ15NNO3 value of +12 ‰ ± 1 ‰. Headwaters and catchments with forest

270

land use show a variation range of δ15NNO3 from 0 ‰ to about 8 ‰ with low mean nitrate

271

concentrations of 4.4 mg L-1. These values are slightly higher than in other studies which report

272

δ15NNO3 in forested watersheds between -1 ‰ and +5 ‰35. For forested catchments with more

273

than 20 % of agricultural land use, a positive shift of the N isotope signatures up to +11 ‰

274

indicating the impact of N sources related to agricultural practice including waste water.

275

Generally, nitrate sources from agricultural land use are characterized by relatively elevated

276

nitrogen isotopic signatures up to +25 ‰13 which are obvious for agricultural dominated

277

tributaries (Figure 5C and E). Synthetic fertilizers show only insignificant impact on the isotopic ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 29

278

δ15N signature of the stream network, because they could only lower nitrogen isotopic values

279

towards typical signatures between -3 ‰ and +3 ‰ and corresponding oxygen isotopic signature

280

between +23 ‰ and +24 ‰13,37. Increasing inflow of urban or agricultural waste water into the

281

stream network (Figure 5E) is assumed to result in a higher pollution especially in the lower

282

Bode region36. On the contrary, the amount of waste water for the major Selke stream just

283

accounts 2.9 % of the entire nitrate loads36.

284

The oxygen isotopic signature of nitrate at the sub-catchment outlets varies between -0.3 and

285

+9.8 ‰. The main river is a mixture of these signatures and varies between +2.9 and +6.7 ‰.

286

Large δ18ONO3-variations in small catchments can be explained either by reactive isotope

287

fractionation processes during nitrate formation and their dependency on temperature and

288

hydrological conditions or by a minor yet variable direct input of atmospheric nitrate (δ18O-NO3

289

+75 ‰ for the Bode region after Mueller et al.21) into the stream system. Generally, the oxygen

290

isotopic signatures of nitrate do not show a significant correlation with major land use types

291

(Figure 5B). However, it can be noted that, especially in catchments with agricultural land use,

292

elevated N isotope signatures are accompanied by slightly elevated O isotope signatures (dual

293

nitrate plot in Figure 5E). This may indicate a certain impact of bacterial denitrification as a

294

nitrate consuming process occurring in groundwater feeding the springs or in the stream itself38.

295

Correlations between isotope signatures and nitrate concentration data do not provide additional

296

evidence for the occurrence of denitrification. This might be due to the fact that the variation of

297

the initial nitrate concentrations exceeds the concentration impact of denitrification.

298

Nitrate dynamics related to different hydrological conditions. Slightly elevated oxygen

299

isotope signatures suggest that only a small share of stream water nitrate in headwater catchments

300

directly derives from atmospheric deposition. Groundwater feeding the springs during base flow ACS Paragon Plus Environment

Page 15 of 29

Environmental Science & Technology

301

conditions indicate recycling (immobilization with subsequent mineralization) of atmospheric

302

nitrate within the soil zone. During high flow events (Figure 6A), oxygen isotopic signatures of

303

nitrate follow a similar pattern as during base flow conditions (Figure 6B). Elevated nitrate

304

concentrations during high flow, however, indicate a more intense mobilization of nitrified soil

305

nitrogen that originates either from the atmosphere or from the decay of organic material.

306

Catchment scale assessment of predominant nitrogen transformation zones. To assess nitrate

307

dynamics in a nearly natural environment we selected the Selke River, a tributary of the Bode, as

308

a special study site. Within this catchment we try to identify the influence of changing land use

309

along the river course on nitrate isotopic signatures (Figure 1). According to Jiang et al.37 and

310

Musolff et al.19, the mesoscale Selke River shows increasing inorganic nitrogen (IN as

311

ammonium, nitrate or nitrite) concentrations with rising runoff. Furthermore, the Selke region is

312

affected by high temporal and spatial variations of IN leaching, influenced by hydrological

313

transport processes37. The mean concentration of inorganic nitrogen along the Selke River

314

increases downstream at the discharge stations (S2 to S3 to S4) from 1.44 to 1.75 to 3.91 mg L-1

315

37

. Corresponding nitrate isotopic signatures show variations from +4 to +13 ‰ for δ15NNO3

316

(Figure 7A) and 0 to +5 ‰ for δ18ONO3 (Figure 7B), respectively. Generally, the nitrate isotopic

317

signature increases downstream the Selke River from gauges S2 to S3 to S4. This is due to the

318

increasing impact of agricultural land use. Coefficients of determination (R2) for the correlation

319

between measured discharge and N-isotope signatures of nitrate range between 0.30 and 0.47 for

320

Selke River stations. Lower R2 values for Bode River stations (0.16 to 0.20) indicate an even

321

weaker correlation. Even though weak, the observed correlation between discharge and N-isotope

322

signatures may in part be related to an advanced high flow mobilization of N from the

323

nitrification of reduced or organic soil nitrogen (as indicated by a positive correlation between

ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 29

324

discharge and N-concentrations, see: Jiang et al.37 and Musolff et al.19) that relatively less

325

actively emits nitrate during periods of low flow.

326

Even though the coefficients of determination (R2) for the correlation between measured

327

discharge and O-isotope signatures of nitrate are very low (0.12 to 0.49 for Selke River stations),

328

a relationship seems obvious in a way that higher discharge is accompanied by decreased δ18ONO3

329

values. Such trend might in part be due to the enhanced mobilization of isotopically lighter nitrate

330

formed by bacterial nitrification processes of reduced or organic soil nitrogen during periods of

331

higher hydrologic activity and the resulting higher wash out. During nitrification, the oxygen

332

isotope signature of nitrate is controlled by the incorporation of O-H2O and O-O2 (+23.5 ‰) into

333

the nitrate molecule. Normally, two O atoms are derived from ambient water while one O atom

334

comes from dissolved oxygen. However, due to variable kinetic and equilibrium fractionation

335

during O incorporation an additional uncertainty of up to 5 ‰ is introduced to the 2:1 ratio39.

336

Such different discharge scenarios and the related differences of the residence time in the reactive

337

zone, where nitrification is taking place, may have significant impact on the isotope fractionation

338

during O-incorporation. Consequently, oxygen isotope signatures of nitrate depend in part on the

339

discharge intensity (Figure 7B). The isotopic signature of NO3 in precipitation does not affect

340

significantly the isotopic signature of the Selke tributary because it would increase δ18ONO3

341

values21.

342

Beyond a different mobilization scheme during periods of higher flow, the observed correlation

343

pattern between isotope signatures, nitrate concentrations and discharge may in part be due to the

344

in-stream turnover of nitrate either by denitrification or by uptake. During low flow conditions, a

345

more intensive denitrification following oxygen depletion of stream water nitrate with biofilms

346

on the stream bed substrate would increase the impact of the nitrate turnover and result in

ACS Paragon Plus Environment

Page 17 of 29

Environmental Science & Technology

347

elevated isotope signatures of both nitrogen and oxygen in nitrate (Fig. 7 A, B) and lower nitrate

348

concentrations37. The slope of the regression line between oxygen and nitrogen isotope signatures

349

measured at three gauging stations of the Selke River is close to 0.5, a value that is considered to

350

be characteristic for denitrification in terrestrial aquatic systems13. For the two Bode River

351

gauging stations this correlation becomes very weak as mixing processes of nitrate from all

352

tributaries mask the potential isotopic effect of in-stream nitrate turnover.

353

As a consequence of the general dependency of the variability of isotope signatures on discharge

354

intensity, a reliable resolution of source and process dynamics is only possible for smaller sub-

355

catchments. While an evaluation of the isotope data set with respect to source delineation and

356

certain target processes such as denitrification is still possible for the Selke River catchment, an

357

increasing impact of mixing processes on the isotope signal in the main river does not allow for a

358

differentiated analysis of the N dynamics anymore. After all, our study reveals that the nitrate

359

degradation potential in groundwater or stream water is insufficient to control the nitrate export

360

from the Bode River watershed. Therefore, nitrate loads in the river system can only be regulated

361

by preventive measures especially by an optimized land use management.

362

AUTHOR CONTRIBUTIONS

363

The manuscript was drafted by Christin Mueller, who performed the sampling, partly laboratory

364

analysis and interpretation of the data under the direct supervision of Kay Knöller and Ralf Merz.

365

Michael Rode gave useful information for interpreting the isotope results and a model approach.

366

Ronald Krieg was continuously involved in the sample collection and characterization of related

367

parameters. Luis Samaniego and Matthias Zink developed and performed the mesoscale

368

Hydrological Model.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 29

369

ACKNOWLEDGMENT

370

For the assistance in the field and laboratory we would like to thank Julia Jäger, Benedikt Ahner

371

and Stefanie Kolbe. Furthermore, we give our sincere thanks to Martina Neuber, Silke Köhler,

372

Petra Blümel and Wolfgang Städter for running numerous chemical and isotope analyses. We

373

kindly acknowledge our data providers, the German Weather Service (DWD), the Joint Research

374

Center of the European Commission, the European Environmental Agency, the Federal Institute

375

for Geosciences and Natural Resources (BGR), the European Water Archive and the Global

376

Runoff data Centre. The project is supported by TERENO (Terrestrial Environmental

377

Observatories of the Helmholtz Association). The River Network is supported by ATKIS R DLM

378

1000 © Bundesamt für Kartographie und Geodäsie, 2003. Topographical information is based on

379

Geobasisinformation of the Vermessungs- und Katasterverwaltung by permission of Landesamt

380

für Landesvermessung und Datenverarbeitung Sachsen-Anhalt Gen.-Nr. I, Verm D/P/086/95.

381

SUPPORTING INFORMATION

382

Information about biochemical processes of natural and anthropogenic nitrate sources in the N-

383

cycle can be found in Supporting information (Text S1). Table S1 shows selected nitrogen

384

enrichment factors (ε) with relation to associated nitrate-oxygen ratios.

385

REFERENCES

386 387 388 389 390 391 392 393

1. Galloway, J. N.; Aber, J. D.; Erisman, J. W.; Seitzinger, S. P.; Howarth, R. W.; Cowling, E. B.; Cosby, B. J., The Nitrogen Cascade. Bioscience 2003, 53, (4), 341-356. 2. Zweimüller, I.; Zessner, M.; Hein, T., Effects of climate change on nitrate loads in a large river: the Austrian Danube as example. Hydrological Processes 2008, 22, (7), 1022-1036, DOI 10.1002/hyp.7000. 3. Pattinson, S. N.; Garcıá -Ruiz, R.; Whitton, B. A., Spatial and seasonal variation in denitrification in the Swale–Ouse system, a river continuum. Science of The Total Environment 1998, 210–211, (0), 289-305, DOI 10.1016/S0048-9697(98)00019-9.

ACS Paragon Plus Environment

Page 19 of 29

394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440

Environmental Science & Technology

4. Whitehead, P. G.; Wilby, R. L.; Battarbee, R. W.; Kernan, M.; Wade, A. J., A review of the potential impacts of climate change on surface water quality. Hydrological Sciences Journal 2009, 54, (1), 101-123, DOI 10.1623/hysj.54.1.101. 5. Wagenschein, D.; Rode, M., Modelling the impact of river morphology on nitrogen retention—A case study of the Weisse Elster River (Germany). Ecological Modelling 2008, 211, (1–2), 224-232, DOI 10.1016/j.ecolmodel.2007.09.009. 6. Mulholland, P. J.; Helton, A. M.; Poole, G. C.; Hall, R. O.; Hamilton, S. K.; Peterson, B. J.; Tank, J. L.; Ashkenas, L. R.; Cooper, L. W.; Dahm, C. N.; Dodds, W. K.; Findlay, S. E. G.; Gregory, S. V.; Grimm, N. B.; Johnson, S. L.; McDowell, W. H.; Meyer, J. L.; Valett, H. M.; Webster, J. R.; Arango, C. P.; Beaulieu, J. J.; Bernot, M. J.; Burgin, A. J.; Crenshaw, C. L.; Johnson, L. T.; Niederlehner, B. R.; O/'Brien, J. M.; Potter, J. D.; Sheibley, R. W.; Sobota, D. J.; Thomas, S. M., Stream denitrification across biomes and its response to anthropogenic nitrate loading. Nature 2008, 452, (7184), 202-205, DOI 10.1038/nature06686. 7. Birgand, F.; Skaggs, R. W.; Chescheir, G. M.; Gilliam, J. W., Nitrogen Removal in Streams of Agricultural Catchments—A Literature Review. Critical Reviews in Environmental Science and Technology 2007, 37, (5), 381-487, DOI 10.1080/10643380600966426. 8. von Schiller, D.; Martí, E.; Riera, J.; Ribot, M.; Argerich, A.; Fonollà, P.; Sabater, F., Inter-annual, Annual, and Seasonal Variation of P and N Retention in a Perennial and an Intermittent Stream. Ecosystems 2008, 11, (5), 670-687, DOI 10.1007/s10021-008-9150-3. 9. Hauck, R. D., Nitrogen Tracers in Nitrogen Cycle Studies—Past Use and Future Needs1. J. Environ. Qual. 1973, 2, (3), 317-327, DOI 10.2134/jeq1973.00472425000200030001x. 10. Heaton, T. H. E., Isotopic studies of nitrogen pollution in the hydrosphere and atmosphere: A review. Chemical Geology: Isotope Geoscience section 1986, 59, (0), 87-102, DOI 10.1016/0168-9622(86)90059-X.. 11. Durka, W.; Schulze, E.-D.; Gebauer, G.; Voerkeliust, S., Effects of forest decline on uptake and leaching of deposited nitrate determined from 15N and 18O measurements. Nature 1994, 372, (6508), 765-767, DOI 10.1038/372765a0. 12. Wankel, S. D.; Kendall, C.; Pennington, J. T.; Chavez, F. P.; Paytan, A., Nitrification in the euphotic zone as evidenced by nitrate dual isotopic composition: Observations from Monterey Bay, California. Global Biogeochemical Cycles 2007, 21, (2), DOI 10.1029/2006GB002723. 13. Kendall, C.; McDonnell, J. J. E., Tracing Nitrogen Sources and Cycling in Catchments. Isotope Tracers in Catchment Hydrology, ; Kendall, C., McDonnell, J. J., Eds. 1998, Elsevier Science B.V., Amsterdam, , pp 519ff. 14. Kendall, C.; Elliott, E. M.; Wankel, S. D., Tracing Anthropogenic Inputs of Nitrogen to Ecosystems. In Stable Isotopes in Ecology and Environmental Science, Blackwell Publishing Ltd: 2008; pp 375-449, DOI 10.1002/9780470691854.ch12.. 15. Boyacioglu, H.; Vetter, T.; Krysanova, V.; Rode, M., Modeling the impacts of climate change on nitrogen retention in a 4th order stream. Climatic Change 2012, 113, 18, DOI 10.1007/s10584-011-0369-1. 16. Wexler, S. K.; Goodale, C. L.; McGuire, K. J.; Bailey, S. W.; Groffman, P. M., Isotopic signals of summer denitrification in a northern hardwood forested catchment. Proceedings of the National Academy of Sciences 2014, 111, (46), 16413-16418, DOI 10.1073/pnas.1404321111. 17. Aber, J. D.; Goodale, C. L.; Ollinger, S. V.; Smith, M.-L.; Magill, A. H.; Martin, M. E.; Hallett, R. A.; Stoddard, J. L., Is Nitrogen Deposition Altering the Nitrogen Status of Northeastern Forests? Bioscience 2003, 53, (4), DOI 10.1641/00063568(2003)053[0375:indatn]2.0.co;2.

ACS Paragon Plus Environment

Environmental Science & Technology

441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488

Page 20 of 29

18. MacDonald, J. A.; Dise, N. B.; Matzner, E.; Armbruster, M.; Gundersen, P.; Forsius, M., Nitrogen input together with ecosystem nitrogen enrichment predict nitrate leaching from European forests. Global Change Biology 2002, 8, (10), 1028-1033, DOI 10.1046/j.13652486.2002.00532.x. 19. Musolff, A.; Schmidt, C.; Selle, B.; Fleckenstein, J. H., Catchment controls on solute export. Advances in Water Resources 2015, 86, 13, DOI 10.1016/j.advwatres.2015.09.026. 20. Zacharias, S.; Bogena, H.; Samaniego, L.; Mauder, M.; Fuß, R.; Pütz, T.; Frenzel, M.; Schwank, M.; Baessler, C.; Butterbach-Bahl, K.; Bens, O.; Borg, E.; Brauer, A.; Dietrich, P.; Hajnsek, I.; Helle, G.; Kiese, R.; Kunstmann, H.; Klotz, S.; Munch, J. C.; Papen, H.; Priesack, E.; Schmid, H. P.; Steinbrecher, R.; Rosenbaum, U.; Teutsch, G.; Vereecken, H., A Network of Terrestrial Environmental Observatories in Germany. Vadose Zone Journal 2011, 10 (3), (Soil Science Society of America, Madison, WI.), 19, DOI 10.2136/vzj2010.0139. 21. Mueller, C.; Krieg, R.; Merz, R.; Knoller, K., Regional nitrogen dynamics in the TERENO Bode River catchment, Germany, as constrained by stable isotope patterns. Isotopes Environ Health Stud 2015, 52 (1-2), 16-74, DOI 10.1080/10256016.2015.1019489. 22. EEA, Land Cover. European Environment Agency. 2012, Accessed November 2013. 23. Casciotti, K. L.; Sigman, D. M.; Hastings, M. G.; Bohlke, J. K.; Hilkert, A., Measurement of the oxygen isotopic composition of nitrate in seawater and freshwater using the denitrifier method. Anal Chem 2002, 74, (19), 4905-12, DOI 10.1021/ac020113w. 24. Sigman, D. M.; Casciotti, K. L.; Andreani, M.; Barford, C.; Galanter, M.; Böhlke, J. K., A Bacterial Method for the Nitrogen Isotopic Analysis of Nitrate in Seawater and Freshwater. Analytical Chemistry 2001, 73, (17), 4145-4153, DOI 10.1021/ac010088e. 25. Oudin, L.; Andréassian, V.; Perrin, C.; Michel, C.; Le Moine, N., Spatial proximity, physical similarity, regression and ungaged catchments: A comparison of regionalization approaches based on 913 French catchments. Water Resources Research 2008, 44, (3), W03413, DOI 10.1029/2007WR00624. 26. Samaniego, L.; Kumar, R.; Attinger, S., Multiscale parameter regionalization of a gridbased hydrologic model at the mesoscale. Water Resources Research 2010, 46, (5), W05523, DOI 10.1029/2008WR007327. 27. Bergström, S., Development and application of a conceptual runoff model for Scandinavian catchments. Bulletin Series A 1976, 52, (Norrköping, Sweden: SMHI). 28. Liang, X.; Lettenmaier, D. P.; Wood, E. F.; Burges, S. J., A simple hydrologically based model of land surface water and energy fluxes for general circulation models. Journal of Geophysical Research: Atmospheres 1994, 99, (D7), 14415-14428, DOI 10.1029/94JD00483. 29. Tewolde, M. H.; Smithers, J. C., Flood routing in ungauged catchments using Muskingum methods. Water SA 2006, 32, (3), 379-388. DOI 10.4314/wsa.v32i3.5263.. 30. Nash, J. E.; Sutcliffe, J. V., River flow forecasting through conceptual models part I — A discussion of principles. Journal of Hydrology 1970, 10, (3), 282-290, DOI 10.1016/00221694(70)90255-6. 31. Tolson, B. A.; Shoemaker, C. A., Dynamically dimensioned search algorithm for computationally efficient watershed model calibration. Water Resources Research 2007, 43, (1), W01413, DOI 10.1029/2005WR004723. 32. Kumar, R.; Livneh, B.; Samaniego, L., Toward computationally efficient large-scale hydrologic predictions with a multiscale regionalization scheme. Water Resources Research 2013, 49, (9), 5700-5714, DOI 10.1002/wrcr.20431. 33. Moriasi, D. N.; Arnold, J. G.; Liew, M. W. V.; Bingner, R. L.; Harmel, R. D.; Veith, T. L., Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Transactions of the ASABE 2007, Vol. 50(3), 15. DOI 10.13031/2013.23153. ACS Paragon Plus Environment

Page 21 of 29

489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506

Environmental Science & Technology

34. Conradt, T.; Roers, M.; Schroter, K.; Elmer, F.; Hoffmann, P.; Koch, H.; Hattermann, F. F.; Wechsung, F., Comparison of the extreme floods of 2002 and 2013 in the German part of the Elbe River basin and their runoff simulation by SWIM-live. Hydrol. Wasserbewirtsch. 2013, 57, (5), 241-245, DOI 10.5675/HyWa-2013,5-4. 35. Nestler, A.; Berglund, M.; Accoe, F.; Duta, S.; Xue, D.; Boeckx, P.; Taylor, P., Isotopes for improved management of nitrate pollution in aqueous resources: review of surface water field studies. Environmental Science and Pollution Research 2011, 18, (4), 519-533, DOI 10.1007/s11356-010-0422-z. 36. Rode, M.; Halbedel née Angelstein, S.; Anis, M. R.; Borchardt, D.; Weitere, M., Continuous In-Stream Assimilatory Nitrate Uptake from High-Frequency Sensor Measurements. Environmental Science & Technology 2016, 50, (11), 5685-5694. 37. Jiang, S.; Jomaa, S.; Rode, M., Modelling inorganic nitrogen leaching in nested mesoscale catchments in central Germany. Ecohydrology 2014, 7, (5), 1345-1362, DOI 10.1002/eco.1462. 38. Mariotti, A.; Mariotti, F.; Champigny, M.-L.; Amarger, N.; Moyse, A., Nitrogen Isotope Fractionation Associated with Nitrate Reductase Activity and Uptake of NO(3)− by Pearl Millet. Plant Physiology 1982, 69, (4), 880-884. 39. Buchwald, C.; Casciotti, K. L., Isotopic ratios of nitrite as tracers of the sources and age of oceanic nitrite. Nature Geosci 2013, 6, (4), 308-313, DOI 10.1038/ngeo1745.

ACS Paragon Plus Environment

Environmental Science & Technology

TOC TOC 142x80mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29

Environmental Science & Technology

Figure 1. Location of the Bode catchment with major sub-catchments and spatial distribution of stream discharge gauges. The Selke tributary is heightened with runoff stations and land use information22. Figure 1 570x413mm (72 x 72 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 2. Simulated versus observed discharge at the Selke tributary outlet (S4) and the Bode River outlet (B5). Seasonal (black arrow) and monthly (grey arrow) sampling campaigns are indicated on the temporal scale. Figure 2 371x292mm (72 x 72 DPI)

ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29

Environmental Science & Technology

Figure 3 A. Annual mean discharge distribution of major Bode River tributaries in 2012 and 2013; B. Distribution of calculated nitrate loads based on measured discharge and nitrate concentration Figure 3 709x384mm (72 x 72 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 4 A. Trend between specific discharge and specific nitrate load the stream differentiated between major land use types. B Relationship between calculated portion of agricultural land use for each subcatchment and the coefficient of variation cv in the nitrate concentration. Additionally, the spatial distribution of major land cover types and corresponding cv is shown for the entire Bode River catchment. Figure 4 340x511mm (72 x 72 DPI)

ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

Environmental Science & Technology

Figure 5. Relationship between mean δ15N (A) and δ18O (B) for nitrate and portion of agricultural land use; Relationship between specific nitrate loads and nitrate isotopic signatures δ15N (C) and δ18O (D) as well as mean dual nitrate isotope plot (E) with typical NO3- source signatures Figure 5 493x508mm (72 x 72 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 6. Contrasting hydrological conditions (high flow (A) and base flow (B)) with corresponding nitrate concentration and δ18O-NO3- of stream water and precipitation water related to major land-cover type and sub-catchment size. Figure 6 610x649mm (72 x 72 DPI)

ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29

Environmental Science & Technology

Figure 7. Relationship between measured nitrate isotopic signature (δ15N (A) and δ18O (B)) (S2-S3-S4, monthly measurements in 2013) and continuously measured discharge at stations along the tributary Selke River. For comparison, nitrate isotopic signatures and discharge measurements are shown for stations of the Bode River right downstream of the inflow of the Selke River into the Bode River (B4) and on the outlet of the Bode River (B5). Figure 7 415x533mm (72 x 72 DPI)

ACS Paragon Plus Environment