Subscriber access provided by PEPPERDINE UNIV
Article
CO2 outgassing from an urbanized river system fueled by wastewater treatment plant effluents Tae Kyung Yoon, Hyojin Jin, Most Shirina Begum, Namgoo Kang, and Ji-Hyung Park Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b02344 • Publication Date (Web): 22 Aug 2017 Downloaded from http://pubs.acs.org on August 23, 2017
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 32
Environmental Science & Technology
CO2 outgassing from an urbanized river system fueled by wastewater treatment plant effluents Tae Kyung Yoon†,‡, Hyojin Jin†, Most Shirina Begum†, Namgoo Kang§,∥, Ji-Hyung Park†,*
†
Department of Environmental Science and Engineering, Ewha Womans University, Seoul 03760,
Republic of Korea §
Center for Gas Analysis, Korea Research Institute of Standards and Science, Daejeon 34113,
Republic of Korea
∥Science
of Measurement, University of Science and Technology, KRISS Campus, Daejeon
34113, Republic Korea
1
ACS Paragon Plus Environment
Environmental Science & Technology
1
ABSTRACT
2
Continuous underway measurements were combined with a basin-scale survey to examine
3
human impacts on CO2 outgassing in a highly urbanized river system in Korea. While the partial
4
pressure of CO2 (pCO2) was measured at 15 sites using syringe equilibration, three cruises
5
employing an equilibrator were done along a 30-km transect in the Seoul metropolitan area. The
6
basin-scale survey revealed longitudinal increases in surface water pCO2 and dissolved organic
7
carbon (DOC) in the downstream reach. Downstream increases in pCO2, DOC, fluorescence
8
index, and inorganic N and P reflected disproportionately large contributions from wastewater
9
treatment plant (WWTP) effluents carried by major urban tributaries. Cruise transects exhibited
10
strong localized peaks of pCO2 up to 13,000 µatm and 13CO2 enrichment along the confluences
11
of tributaries at average flow, whereas CO2 pulses were dampened by increased flow during the
12
monsoon period. Fluctuations in pCO2 along the eutrophic reach downstream of the confluences
13
reflected environmental controls on the balance between photosynthesis, biodegradation, and
14
outgassing. The results underscore WWTP effluents as an anthropogenic source of nutrients,
15
DOC, and CO2 and their influences on algal blooms and associated C dynamics in eutrophic
16
urbanized river systems, warranting further research on urbanization-induced perturbations to
17
riverine metabolic processes and carbon fluxes.
18
2
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32
Environmental Science & Technology
19
INTRODUCTION
20
As the terrestrial biosphere “breathes”,1 streams, rivers, and lakes release carbon dioxide (CO2)
21
to the atmosphere during transit of organic and inorganic carbon (C) to the oceans, constituting
22
an important link in the global carbon cycle.2-8 Recent global estimates of CO2 outgassing from
23
inland waters range from < 1 to 3.8 Pg C y-1, which can result in a significant loss of the annual
24
terrestrial sink (~2.5 Pg C y-1) of CO2 emitted from anthropogenic sources.4-7, 9 Large
25
uncertainties in current estimates result from a combination of factors including sparse spatial
26
coverage of direct measurements of the partial pressure of CO2 (pCO2) in inland waters and
27
difficulties in estimating gas transfer velocities.6, 7, 9 Although direct measurements of pCO2 have
28
identified previously unrecognized inland water C sources such as wetlands and estuarine
29
marshes,10-12 pCO2 data are still sparse for anthropogenically modified river systems worldwide,
30
except for well-studied rivers in North America13, 14 and Europe.2, 15 Particularly sparse are in-
31
situ measurements in rapidly urbanizing watersheds across Asia, so obtaining pCO2 data in Asian
32
rivers has recently been suggested as a top priority in reducing large uncertainty in estimating the
33
global riverine CO2 outgassing.16
34
Anthropogenic perturbations such as dams, river channeling, and water pollution can alter
35
inland water C fluxes substantially.7 Although anthropogenic changes in water quality and
36
biological communities have been studied extensively,17, 18 few studies have investigated altered
37
C fluxes in urbanized river systems from the perspective of human impacts on aquatic ecosystem
38
processes.19 Some recent studies conducted in urbanized watersheds have identified river
39
impoundment and pollution as key controls on the quantity and quality of dissolved organic
40
matter (DOM).20-22 However, only a small number of studies have linked the altered flux and
3
ACS Paragon Plus Environment
Environmental Science & Technology
41
composition of DOM to direct measurements of pCO2 in urban rivers.14, 23, 24 Large spatial
42
heterogeneity in pCO2 and DOM in urbanized river systems has been attributed to infrastructures
43
such as wastewater treatment plants (WWTPs), dams, weirs, channels, and urban stormwater
44
facilities that may represent discontinuities in the river continuum, altering natural riverine
45
processes and ecosystem resilience.23, 25, 26 Discontinuities in the river continuum can either
46
decrease or increase the rate of riverine CO2 outgassing via enhanced photosynthesis27 or organic
47
matter biodegradation.21 Few studies have addressed these discontinuities and associated large
48
spatial variations in the riverine C fluxes in urbanized watersheds using approaches with
49
adequate spatial resolution.
50
The Han River basin is a highly modified river system draining the middle of the Korean
51
Peninsula including the densely populated Seoul metropolitan area with a population > 25
52
million along the lower reach (Figure 1). Longitudinal variations in dominant land use, along
53
with a cascade of dams along the middle reach, make the Han River basin an ideal venue to
54
study anthropogenic perturbations to hydrologic and biogeochemical flows. A basin-scale survey
55
of pCO2, DOM, and nutrients from the forested headwater stream to the estuary of the Han River
56
was combined with continuous underway measurements of pCO2 employing an automated
57
equilibrator system8 along the downstream reach in the Seoul metropolitan area to explore how
58
river impoundment and pollution alter riverine CO2 outgassing and organic matter transport in
59
the highly urbanized watershed. We hypothesized that urban tributaries delivering high loads of
60
CO2, DOM, and nutrients derived from WWTP effluents to the downstream river could not only
61
enhance CO2 outgassing in the receiving water but also alter the rates of the primary production
62
and organic matter biodegradation along the downstream reach. We expected that enhanced
63
primary production, particularly during algal blooms, could significantly lower the level of pCO2
4
ACS Paragon Plus Environment
Page 4 of 32
Page 5 of 32
Environmental Science & Technology
64
via planktonic uptake of CO2, resulting in large longitudinal variations in pCO2 along the
65
eutrophic downstream reach.
66 67
Figure 1 Study sites in in the Han River basin: (a) a basin-scale transect (~300 km); (b) an
68
underway transect (30 km) and an urban tributary transect (JN; ~35 km) in the Seoul
69
metropolitan area. Three major urban streams (AY: Anyang River, JN: Jungnang River, TC: Tan
70
River), four wastewater treatment plants (WWTPs), and two submerged weirs are marked on the
71
Seoul metropolitan area map (b).
72
73
74
MATERIALS AND METHODS Site description. The 515-km long Han River consists of the North Han and South Han
75
branches and the lower Han River. The river drains an area of 26,142 km2 along the middle of
76
the Korean Peninsula and flows westward into the Yellow Sea (Figure 1; more detailed
77
descriptions of the study site and methods provided in SI).28 The dominant land cover in the Han
78
River basin varies longitudinally from the highly forested upstream reach to increasingly
79
urbanized areas along the mid- to downstream reaches. The Han River flowing through the Seoul 5
ACS Paragon Plus Environment
Environmental Science & Technology
80
metropolitan area is channelized with two submerged weirs at the east (RKM 39) and west end
81
(RKM 70) of Seoul (Figure 1b). Three major urban tributaries feed into the lower reach along the
82
city of Seoul: the Tan River (TC), Jungnang River (JN), and the Anyang River (AY) enter the
83
mainstem at RKM 68, 64, and 49, respectively (Figure 1b). Three WWTPs discharge at the rate
84
of 1.5 million m3 d-1 to JN, with the bulk (1.3 million m3 d-1) discharged from the WWTP located
85
near the mouth. Four WWTPs located within Seoul release effluents from tertiary treatments
86
including modified Ludzack Ettinger (MLE) and anaerobic-anoxic/oxic process (A2O) at the rate
87
of 4.3 million m3 d-1.29
88
Basin-scale survey. A basin-scale survey along the North Han branch and the lower reach
89
was conducted at 15 locations on 4–6 May 2015 at average river flow (183 m3 s-1 at RKM 57;
90
Figure S1). The sites were selected at intervals from a forested headwater stream (RKM 299)
91
through the middle and lower reaches impounded by dams or weirs to the estuarine reach (RKM
92
23) (Figure 1a, Table S1). A 35-km transect of the urban tributary JN, together with WWTP
93
effluents near the confluence, was surveyed at eight sampling points from the headwater to the
94
confluence on 13 November 2015 and 12–13 May 2016. In addition to air temperature and
95
barometric pressure measured with a portable sensor (Watchdog 1650 Micro Station, Spectrum
96
Technologies, USA), we measured in situ water temperature, pH, electrical conductivity, and DO
97
using a portable pH meter (Orion 5-Star Portable, Thermo Scientific, USA). The surface water
98
pCO2 at 20 cm depth was measured by a manual headspace equilibration method.8, 30, 31 The
99
equilibrated headspace air sample, as well as an ambient air sample collected for equilibration,
100
was injected into a GC (7890A, Agilent, USA) fitted with a Supelco Hayesep Q 12 ft 1/8”
101
column for the measurement of CO2 concentration. Grab water samples collected 20 cm below
102
the surface were filtered (GF/F, Whatman) at the laboratory and then analyzed for dissolved
6
ACS Paragon Plus Environment
Page 6 of 32
Page 7 of 32
Environmental Science & Technology
103
organic C (DOC; TOC-VCPH, Shimadzu, Japan), fluorescence excitation-emission matrices
104
(EEMs; F7000, Hitachi, Japan), UV absorbance (8453, Agilent, USA), major ions (883 Basic IC
105
plus, Metrohm, Switzerland), and chlorophyll a (Chl-a).32 Fluorescence and UV absorbance data
106
were used to calculate fluorescence index (FI), as an optical index that can distinguish the origin
107
of DOM between terrestrial/allochthonous sources (1.3–1.5) and microbial/autochthonous
108
sources (1.7–2.0).33-35
109
Continuous underway measurements. Continuous underway measurements of pCO2 were
110
performed on three cruises along the lower reach (RKM 70 to 39; Fig. 1), twice at average flow
111
(212 m3 s-1 on both 11 May 2015 and 10 June 2016,) and once at high flow (757 m3 s-1 on 28
112
July 2015). The boat speed was maintained at ∼10 km h-1 based on the prior tests and other
113
studies.8, 36 To compare spatial differences in pCO2 during a selected day period between three
114
cruise transects, the measurements were conducted within four hours in the afternoon when the
115
mean pCO2 is close to the daily mean value based on our year-long continuous pCO2
116
measurements at RKM 53 (unpublished data). pCO2 was continuously measured using a spray-
117
type equilibrator8, 10, 37, 38 connected to an infrared gas analyzer (IRGA; LI820, Li-Cor, USA).
118
The equilibrated headspace air in the equilibrator chamber was circulated through a channel
119
including an air filter, a desiccant (Drierite) column, and the IRGA. The data were logged every
120
1 s on a laptop. Basic water quality parameters were measured in the water from an additional
121
water pump channel onboard, using a portable multiparameter meter (Orion 5-Star Portable,
122
Thermo Scientific, USA or 6820 V2, YSI Inc., USA). Air temperature and barometric pressure
123
were recorded in a micro-logger (Watchdog 1650 Micro Station, Spectrum Technologies, USA).
124
The cruises were tracked by a GPS tracking unit (Montana 650, Garmin Ltd., USA). The gas
125
samples of the manual headspace and the water samples were collected onboard and at the urban
7
ACS Paragon Plus Environment
Environmental Science & Technology
126
streams near the confluences, and later analyzed in the laboratory with the same methods used
127
for the basin-scale field survey. In the June 2016 transect, stable C isotope ratios of CO2 (δ13CO2)
128
in the gas samples were analyzed at the UC Davis Stable Isotope Facility using a GasBench-
129
IRMS (ThermoScientific, Bremen, Germany). The actual water saturated pCO2 in the
130
equilibrator was calculated from the pCO2 measured in dry air, barometric pressure, water
131
temperature, and equilibrium water vapor (refer to the supporting information; Eq. S1).39
132
Aqueous pCO2 in the equilibrator was then corrected to in situ temperature using a temperature
133
effect coefficient. To compensate the time lags of the measured pCO2 resulting from the
134
residence time of the water circuit and the response time of the sensors, the logged values were
135
corrected according to Eq. S2.36
136
CO2 outgassing estimation and downstream mass balance of CO2 and DOC. The CO2
137
flux at the water-air interface (fCO2) was estimated for the mainstem section along the city of
138
Seoul based on the differences in measured pCO2 between water and air and an estimated gas
139
transfer velocity (K). Our observations and available data were insufficient to apply empirical K
140
models,40 so a K value was adopted from a study conducted in the Hudson River (5 cm h-1)41
141
based on the similarities in river hydrology and climatic conditions. The downstream evolution
142
of CO2 and DOC was simulated using a simple mass balance model adapted from the modeling
143
approach used for estimating CO2 fluxes from the Amazon river and wetlands.12 The model
144
consists of inputs of CO2 and DOC from the upper reach (RKM 69) and from the urban
145
tributaries (TC and JN); and CO2 outgassing from the section from RKM 69 to 48 (Eq. S3, S4).
146
The proportion (termed “residual”) of the total C measured at RKM 62 and 48 that was not
147
accounted for by the upstream and/or tributary inputs of CO2 and DOC and CO2 outgassing from
148
the given section was assumed to represent either CO2 production (i.e., gain) or consumption (i.e.,
8
ACS Paragon Plus Environment
Page 8 of 32
Page 9 of 32
Environmental Science & Technology
149
loss) resulting from organic matter biodegradation or photosynthesis, respectively. All data
150
processing, statistical analyses, and modeling were conducted using R.42
151
RESULTS AND DISCUSSION
152
Spatial pattern of pCO2 and DOC across the Han River basin. The basin-wide survey
153
revealed distinct downstream increases in pCO2, DOC, and FI, particularly from the point where
154
the river flows into the city boundary of Seoul (Figures 2a–c). The surface water pCO2 ranged
155
from 100 to 730 µatm across the upper and middle reaches and began to exceed 1,500 µatm in
156
the upper part of the lower reach (Figure 2a). Relatively low pCO2 values in the vicinity of the
157
dams (288 µatm; N = 7) may indicate an enhanced algal uptake of CO2 in the impounded reaches.
158
Compared to the upstream values, DOC and FI increased 1.5–2-fold in the highly urbanized
159
lower reach (Figure 2b, c). The three major urban tributaries exhibited distinctively higher values
160
of pCO2, DOC, and FI than all the mainstem sites. This finding suggests that the urban tributaries
161
and WWTP effluents, in combination with in-stream metabolic processes in the downstream
162
reach, can significantly alter the fluxes of CO2 and DOC in the urbanized river system.
163
The transect study along the urban tributary JN revealed longitudinal increases in CO2, DOC,
164
and inorganic nutrients resulting in large exports into the confluence (Figures 2e-h). The
165
effluents from a WWTP located 4 km upstream of the JN confluence were highly enriched in
166
CO2, DOC, and inorganic nutrients (Figures 2e-h, Table S2). pCO2, DOC, and FI measured at the
167
most downstream location of the tributary were similar to those measured at the WWTP effluents,
168
indicating that the WWTP effluents account for the bulk of the tributary export of carbon and
169
other nutrients to the Han River. Along the tributary upstream of the WWTP, spanning from the
170
upper, rural area (> 28 km upstream from the confluence) to the middle, urbanized area (4–28
9
ACS Paragon Plus Environment
Environmental Science & Technology
171
km), pCO2 ranged from 1,000–3,000 µatm. Upon receiving the WWTP effluent, the pCO2 in the
172
receiving water increased 4-fold, approaching 8,000 µatm.
Page 10 of 32
173 174
Figure 2 Spatial variations in pCO2, dissolved organic carbon (DOC), and fluorescence index (FI)
175
along a basin-scale transect (a–d) and an urban tributary (Jungnang River; JN on Figure 1)
176
transect (e–h). RKM refers to the river km from the mouth. The measurements in the wastewater
177
treatment plant (WWTP) effluents are indicated by red box plots (e-g).
178
179
The spatial pattern of DOC along the JN transect was similar to the pCO2 pattern (Figure 2f).
180
The FI values indicated a gradual shift from the dominance of terrestrial DOM (FI < 1.5) in the
181
upper reach (Figure 2g) to the greater contribution of autochthonous DOM (FI > 1.7) in the
182
WWTP effluent and its downstream reach.21, 34 High loads of CO2 and DOC in the WWTP
10
ACS Paragon Plus Environment
Page 11 of 32
Environmental Science & Technology
183
effluents likely originate from the wastewater treatment processes including microbial oxidation
184
and biomass respiration that transform organic compounds to CO2 and DOM.43, 44 The WWTP
185
effluents strongly influenced other water quality variables, especially inorganic nutrients (Table
186
S2).
187
Continuous underway measurements of pCO2 along the lower reach. Three cruise
188
surveys along the lower reach of the Han River exhibited unique spatial patterns of pCO2, DOC,
189
and inorganic nutrients that reflect the material inputs from tributaries and WWTPs (Figures 3, 4;
190
Table S3). Mean pCO2 values of 1,066 µatm, 2,518 µatm, and 2,566 µatm, for May 2015, July
191
2015, and June 2016 transects, respectively, were all above the atmospheric saturation level (409
192
µatm), pointing to the role of the lower reach as a net source of CO2. The levels of
193
supersaturation are similar to or higher than the reported values for some temperate rivers,14, 45
194
but lower than those for subtropical and tropical rivers.23, 46, 47 Two transects at average flow
195
(May 2015 and June 2016) showed strong localized pulses of pCO2 along the tributary JN
196
confluence (Figures 3a,c) and noticeable downstream increases in DOC, FI, and inorganic N and
197
P concentrations (Figure 4; Table S3). pCO2 increased from < 100 µatm at a location upstream of
198
the TC confluence (RKM 69) to ~3,000 µatm in the downstream reach of the JN confluence
199
(RKM 62). Strong localized pulses of pCO2 up to ~13,000 µatm were observed in the JN stream
200
plume mixing with the mainstem along RKM 62–65, reflecting large inputs of CO2 delivered by
201
the tributary (Figures 3, S2). The observed localized pCO2 peaks along the tributary confluences
202
have rarely been observed by the previous continuous underway measurements conducted in
203
urbanized rivers and estuaries14, 45 or systems with lower proportions of urban land use.46-48 In a
204
rare case, Zhai et al.49 found exceptionally high pCO2 values approaching 5,000 µatm along an
205
upper reach of the Pearl River estuary with a polluted tributary.
11
ACS Paragon Plus Environment
Environmental Science & Technology
206 207
Figure 3 Underway measurements of pCO2 along the 30-km lower reach of the Han River
208
dissecting the Seoul metropolitan area. The gray lines indicate the measurements along the
209
confluences receiving the urban tributary discharge (a, c). The May 2015 measurements along
210
the mainstem were modified from a prior publication.8
211
12
ACS Paragon Plus Environment
Page 12 of 32
Page 13 of 32
Environmental Science & Technology
212 213
Figure 4 Urban tributary impacts on downstream water quality: pCO2 (a), dissolved organic
214
carbon (DOC; b), fluorescence index (FI; c), chlorophyll-a (Chl-a; d), inorganic nutrients (e–g),
215
and dissolved oxygen (DO; h). The upstream (76–68 km; N = 6) and downstream (62–48 km; N
216
= 12) reaches were compared with two urban streams (TC; N = 2–3) and (JN; N = 5–8). Asterisks
217
(*) and daggers (†) indicate significant (* P < 0.05, ** P < 0.01, and *** P < 0.001) differences
218
in space (before and after the confluence; α) and time (sampling season; β), respectively (Yijk = µ
219
+ αi + βj + (αβ)ij + εijk).
220
221
DOC concentrations increased by 45% and 8% across the transition from the up- to the
222
downstream location of the TC and JN confluences, respectively (Figure 4b; Table S3). FI,
223
inorganic N and P, EC, and Chl-a also showed general longitudinal increases, consistent with the
224
basin-scale survey. The longitudinal increases along the reach downstream of the two tributaries
225
TC and JN were greater at average flow than at high flow, indicating a homogenizing effect of
226
fast and high river flow during the monsoon period. In the May 2015 transect, high
13
ACS Paragon Plus Environment
Environmental Science & Technology
227
concentrations of PO43- and NH4+ peaked at RKM 53 near the JN confluence and rapidly
228
decreased along the downstream reach, likely due to a combination of dilution and rapid algal
229
uptake of labile nutrients during the algal bloom.
230
Observed longitudinal patterns can be explained by the mixing of the mainstem and tributary
231
waters, as indicated by downstream changes associated with distinctively higher pCO2, DOM,
232
and inorganic nutrients and lower DO and Chl-a in the tributaries (Figure 4). For example, pCO2
233
at TC (6240 ± 794 µatm) and JN (9699 ± 1442 µatm) was 80- and 124-fold higher, respectively,
234
than the value (78 µatm) measured at an upstream location (RKM 69). DOC concentrations were
235
2‒3 times higher at the two urban tributaries than the upstream value. The May 2015
236
measurements had particularly high values of pCO2 along the JN confluence and gradually
237
decreased moving downstream from ~3,000 µatm at RKM 62 to ~260 µatm at RKM 53.
238
The measurements of δ13CO2 during the June 2016 cruise exhibited a gradual longitudinal
239
enrichment of 13C along the reach downstream of two tributaries TC and JN (Figure 3c). δ13C
240
increased from −20.9 ‰ at RKM 76 to −16.7 ‰ at RKM 49. The observed values of δ13C fall in
241
the lower range of δ13CO2 and δ13DIC measured for various inland waters (−25 – 0‰).23, 46, 50, 51
242
Interpreting the observed variations in δ13C is challenging because δ13CO2 reflects the interplay
243
between the complex processes involving gas exchange, organic matter decomposition, and
244
dissolution of carbonates of either natural or anthropogenic origin.23, 50 The observed range is
245
similar to those for CO2 originating from the decomposition of organic matter (−25 – −15‰),
246
implying the dominant role of biogenic processes in the river or anthropogenic sources such as
247
WWTPs. 13C enrichment along the reach downstream of the tributaries also suggests that there
248
was a contribution from urban tributaries enriched in 13C (TC: –18.3 ‰, JN: –18.2 ‰, AY: –14.7
14
ACS Paragon Plus Environment
Page 14 of 32
Page 15 of 32
Environmental Science & Technology
249
‰) relative to the upstream value (–20.9 ‰). As suggested by Zeng et al.23, 13C enrichment in
250
DIC can result from a combination of processes including the preferential algal uptake of CO2
251
depleted in 13C, organic matter degradation, and the enhanced air-water gas exchange in the
252
impounded, slow-moving water facilitating the invasion of atmospheric CO2 enriched in 13C
253
(δ13C: –8 ‰) as observed at a location upstream of the TC confluence (RKM 69 and 78), in
254
which pCO2 was undersaturated.
255
In the July 2015 transect during the summer monsoon period, spatial variations in pCO2,
256
DOM, and inorganic nutrients were less pronounced than the measurements during the other
257
transects at lower flows (Figure 3b; Table S3). The higher river flow might have weakened the
258
effects of urban tributary inputs and the spatial variability along the mainstem. Most
259
measurements except pCO2 were lower in magnitude in the monsoon season than at lower flows
260
(Table S3). The lower FI values observed during the monsoon transect (1.41) than those obtained
261
during the low-flow transects (1.55 in 2015 and 1.52 in 2016) can be explained by the greater
262
contribution of DOM from terrestrial sources in the upper reaches. The decrease in FI resulting
263
from an increasing contribution of terrestrial DOM has been observed during the monsoon
264
season across major rivers in Korea.52
265
Downstream transformation of DOC and CO2. The mass balance model consisting of the
266
measurements of pCO2 and DOC in the mainstem and the two urban tributaries (TC and JN) and
267
the estimated rate of CO2 outgassing were used to track changes in major C components along
268
the sections from RKM 69 (upstream of TC confluence) through RKM 62 (downstream of JN
269
confluence) to RKM 48 (upstream of AY confluence) (Figure 5). In May 2015, two tributaries
270
TC and JN accounted for 108% and 47% of pCO2 and DOC concentrations measured at RKM 62,
271
respectively, with the WWTP effluent near the JN confluence alone contributing to 15
ACS Paragon Plus Environment
Environmental Science & Technology
272
approximately one-half of CO2 and one-fifth of DOC carried by the two tributaries. Although we
273
cannot rule out potential overestimation resulting from cumulative errors in estimating the
274
upstream and tributary contributions and CO2 evasion downstream of the confluence, the
275
extremely high tributary pCO2 relative to the very low upstream value indicates the dominant
276
role of the tributary carrying WWTP effluents in increasing pCO2 along the downstream reach at
277
relatively low flow (May 2015 and June 2016) compared to the monsoon period (July 2015). In
278
the monsoon season (July 2015), however, the contribution of the same WWTP effluents was
279
much smaller than at lower flows because of the dilution effect of high river flow. Outgassing
280
from the 14-km downstream transect was estimated to release 38% and 12% of CO2 measured at
281
RKM 62 at average and high flow, respectively. When the rate of outgassing is estimated for the
282
lower reach up to the weir preventing the tidal water intrusion (RKM 39), 55% and 23% of the
283
CO2 inflow from the WWTP near the JN confluence may be released within the 23-km
284
downstream transect (RKM 62 to 39) at average and high flow, respectively.
285
When the fraction associated with CO2 outgassing (fCO2) from the RKM 62-48 section was
286
subtracted from the total C concentration measured at RKM 62, the difference (DOCRKM62 +
287
CO2RKM62 – fCO2RKM48–62) was either higher (May 2015) or lower (July 2015 and June 2016) than
288
the combined concentration of DOC and CO2 measured at RKM 48 (DOCRKM48 + pCO2RKM48).
289
The residual found in May 2015 indicates the consumption of CO2 via planktonic uptake in the
290
14-km transect. In contrast, the increased C concentration at RKM 48 observed in July 2015 and
291
June 2016 can be explained by the production of CO2 in the same section. At high flow (July
292
2015) the concentration of Chl-a varied little across the transect, but increased at average flow
293
between RKM 62 and 48 by 31.4 and 11.8 µg L-1 in May 2015 and June 2016, respectively
294
(Table S3).
16
ACS Paragon Plus Environment
Page 16 of 32
Page 17 of 32
Environmental Science & Technology
295 296
Figure 5 Changes in CO2, dissolved organic carbon (DOC), chlorophyll-a (Chl-a) from a
297
mainstem location (RKM 69) upstream of the first tributary (TC) through a downstream location
298
(RKM 62) of the second tributary (JN) to a location (RKM 48) upstream of the third tributary
299
(AY). Estimates of the tributary input and CO2 outgassing, and residuals are indicated between
300
the measurements at three locations. Note that fCO2 from RKM 69 to RKM 62 was too low (