Accepted Manuscript Impacts of sea cucumber farming on biogeochemical characteristics in the Yellow River estuary, Northern China Jing Fu, Hisashi Yokoyama, Baoshan Cui, Jin Zhou, Jiaguo Yan, Xu Ma, Shozo Shibata PII:
S1474-7065(16)30261-3
DOI:
10.1016/j.pce.2016.12.006
Reference:
JPCE 2543
To appear in:
Physics and Chemistry of the Earth
Received Date: 20 October 2016 Revised Date:
8 December 2016
Accepted Date: 16 December 2016
Please cite this article as: Fu, J., Yokoyama, H., Cui, B., Zhou, J., Yan, J., Ma, X., Shibata, S., Impacts of sea cucumber farming on biogeochemical characteristics in the Yellow River estuary, Northern China, Physics and Chemistry of the Earth (2017), doi: 10.1016/j.pce.2016.12.006. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. 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.
ACCEPTED MANUSCRIPT
Impacts of sea cucumber farming on biogeochemical characteristics in the Yellow River
2‐
estuary, Northern China
3‐
Jing Fu1 , Hisashi Yokoyama2, Baoshan Cui3, Jin Zhou4, Jiaguo Yan3, Xu Ma3, Shozo
4‐
Shibata1
RI PT
1‐
5‐
Japan 2
Field Science Education and Research Center, Kyoto University, Kyoto 6068502, Japan
3
State Key Laboratory of Water Environment Simulation, School of Environment, Beijing
8‐ 9‐ 10‐
Normal University, Beijing 100875, P.R. China 4
12‐
East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences,
Shanghai 200090, P.R. China
TE D
11‐
M AN U
7‐
Graduate School of Global Environmental Studies, Kyoto University, Kyoto 6068501,
SC
1
6‐
13‐ 14‐
EP
15‐
Present address. Graduate School of Global Environmental Studies, Kyoto University,
17‐
Yoshida-Honmachi, Sakyo-ku, Kyoto, 606-8501, Japan. Tel.: +81-80-9472-5668.
18‐
E-mail address:
AC C
16‐
Corresponding author.
19‐
20‐
Email:
[email protected] (Fu J.)
21‐
1‐ ‐
ACCEPTED MANUSCRIPT
22‐
Abstract To investigate the potential environmental effects of pond farming for Apostichopus
24‐
japonicas in Yellow River estuary, we examined discrepancies of distance-based typical
25‐
pollution indicators (TOC, TN, NO3-, NH4+, NO2- and PO43-) and biochemical tracers (δ13C
26‐
and δ15N) in water column and sediment, as well as dietary characteristics of dominant
27‐
macrobenthos between farming and non-farming areas. The results revealed that studied
28‐
variables in water column showed no uniform spatial differences. Meanwhile, those in
29‐
sediments displayed similar decrease tendencies from farming pond to the adjacent tidal
30‐
flat, which was considered to represent the environmental effects of farming. Biochemical
31‐
tracers (δ13C and δ15N) in both water column and sediment confirmed the origin of organic
32‐
matters from the aquaculture waste. The detectable dispersion distance of aquaculture waste
33‐
was restricted to an area within 50 m distance as determined by most variables in sediment
34‐
(TOC, TN, NO3- and NH4+), particularly by C/N ratio and δ13C with which origins of the
35‐
wastes were traced. Bayesian mixing models indicated that in the farming area BMA had a
36‐
larger contribution, while POM(marine) showed a smaller contribution to the diets of Helice
37‐
tridens and Macrophthalmus abbreviates compared to those in the non-farming area. The
38‐
overall results showed that pond farming for Apostichopus japonicus in the Yellow River
39‐
estuary altered the local environment to a certain extent. For methodological consideration,
40‐
sediment biogeochemical characteristics as a historical recorder much more effectively
41‐
reflected aquaculture waste accumulation, and stable isotope approaches are efficient in
42‐
tracing the origin and extent of various allogenous sources.
AC C
EP
TE D
M AN U
SC
RI PT
23‐
43‐ 44‐
Keywords: sea cucumber farming; Apostichopus japonicus; environmental impact;
45‐
aquaculture waste; stable isotope; Yellow River estuary
46‐
2‐ ‐
ACCEPTED MANUSCRIPT
47‐
1. Introduction Coastal aquaculture has rapidly grown in recent years with an increasing demand for
49‐
seafood production (Naylor et al., 2000; Tanner and Fernandes, 2010). The expansion of
50‐
coastal aquaculture, however, has raised widespread concerns on environmental costs and
51‐
sustainable aquaculture techniques (Pusceddu et al., 2007; Tanner and Fernandes, 2010;
52‐
Troell et al., 2003).
SC
RI PT
48‐
Aquaculture farming release large amounts of dissolved and particulate nutrients
53‐
originating from excretion process of reared organisms and uneaten feed and fecal materials
55‐
(Holby and Hall, 1991; Lupatsch and Kissil, 1998). The redundant particulate organic
56‐
carbon and nitrogen loads result in organic enrichment and sediment biogeochemical
57‐
characteristics alternation (Brooks and Mahnken, 2003; Holby and Hall, 1991; Sarà et al.,
58‐
2004). The excess dissolved inorganic nutrients that easily dispersed in the water and
59‐
deposited in the sediment are the most concerned contaminant contributed to coastal
60‐
eutrophication (Paerl, 2006; Smith et al., 1999; Wu et al., 2014).
TE D
M AN U
54‐
In recent years, sea cucumber farming has increased tremendously encouraged by the
EP
61‐
strong and growing market demand (Chen, 2003; FAO, 2008; Han et al., 2016) for its huge
63‐
profit and high nutritional and pharmaceutical values (Bordbar et al., 2011). Consequently,
64‐
much attention had been recently focused on its culturing manner for better production and
65‐
development (Dong et al., 2010; Liu et al., 2010; Paltzat et al., 2008; Yokoyama, 2013).
66‐
Although a variety of studied related to sea cucumber industry, such as hatchery production,
67‐
culturing manner, sociological issues of stocking, were conducted in previous studies (Han
68‐
et al., 2016; Ren et al., 2012), there existed a relative paucity of details concerning the
AC C
62‐
3‐ ‐
ACCEPTED MANUSCRIPT
environmental effects of sea cucumber farming activities (Li et al., 2014; Ren et al., 2010;
70‐
Zheng et al., 2009). One of the most important reason probably lies in that most sea
71‐
cucumbers are deposit feeders. They always are considered as the “environmental cleaners
72‐
or scavengers (Yang et al., 1999) and therefore no exogenous nutrients were introduced.
73‐
Actually, the exogenous nutrients input mainly depend on culturing manners. In most areas
74‐
of China, commercial feeds which mainly composed of fish meal and soybean meal
75‐
(accounting for 15% and 30% of total weight, respectively) were commonly used in sea
76‐
cucumber ponds. The excessive feeds always make a large amount of dissolved and
77‐
particulate waste which result in nutrients enrichment of water and sediment in the pond
78‐
and adjacent environment.‐
M AN U
SC
RI PT
69‐
Sea cucumber farming in China increased tremendously along coasts, and the yield
79‐
reached 205,791 tons (dry weight) in 2015 (China Fisheries Yearbook, 2016). For the
81‐
present study area Yellow River estuary, sea cucumber production exponentially increased
82‐
in the past ten years. The farming ponds area increased from 21.3 ha in 2003 to 21734 ha in
83‐
2014 with an annual production of 11241 tons (dry weight) annually. Among the sea
84‐
cucumber species cultured (nearly ten species) around the world at present, Apostichopus
85‐
japonicus was the largest species in production, which was mainly artificially cultured in
86‐
China (Han et al., 2016; Purcell et al., 2012). The annual hatchery production of juvenile of
87‐
this species was approximately 6 billion in China, and the sum exceeded the total
88‐
production of all other reared sea cucumber species (Purcell et al., 2012).‐The rapid growth
89‐
of Apostichopus japonicus farming promoted a number of studies focusing on genetics (Li
90‐
et al., 2009), energetic (Liu et al., 2009), neurology (SEO and LEE, 2011) and feeding
AC C
EP
TE D
80‐
4‐ ‐
ACCEPTED MANUSCRIPT
ecology and physiology (Sun et al., 2013) etc. However, few studies reported the impacts
92‐
of Apostichopus japonicus farming on pond sediment (Ren et al., 2012; Zheng et al., 2009)
93‐
or coastal waters (Kang and Xu, 2016), no comprehensive investigation of water, sediment
94‐
and biological characteristics was conducted to date .
RI PT
91‐
Environmental variables selected in this research were all universally recognized by
95‐
previous studies. Some of them were specially emphasized as highly informative
97‐
descriptors of heterotrophic loadings from aquaculture, such as dissolved inorganic
98‐
nutrients (DIN) like nitrite (NO3-), nitrate (NO2-), ammonium (NH4+) and phosphate (PO43-)
99‐
in water (Sarà, 2007), available phosphorus (AP) and acid volatile sulphide (AVS) in
M AN U
SC
96‐
sediment (Yokoyama, 2003a). Other studied variables like total organic carbon (TOC) and
101‐
total nitrogen (TN) has been widely used as important indices to represent organic matter
102‐
contents in water and sediment (Gao et al., 2012). Similarly, C/N ratio, δ13C and δ15N were
103‐
commonly utilized in recent years as effective tools to trace and verify the origin, extent
104‐
and fate of nutrient sources in the environment and ecosystem (Andrews et al., 1998b;
105‐
Sakamaki et al., 2010). The combination of these biogeochemical indicators will
106‐
sufficiently describe and interpret the local environment variations (Kollongei and Lorentz,
107‐
2014) induced by aquaculture activities.
AC C
EP
TE D
100‐
The main purposes of this study were to (1) examine the impacts of sea cucumber pond
108‐ 109‐
discharge on local tidal flat by evaluating spatial differences (between farming and
110‐
non-farming area or among stations with a gradient distance to cultured pond) of
111‐
environmental variables in water and sediment; (2) detect differences in feeding strategies
112‐
of dominant macrobenthos between farming and non-farming areas. We hypothesized that
5‐ ‐
ACCEPTED MANUSCRIPT
significant spatial differences exist for contents of selected environmental variables by
114‐
aquaculture activities with varying degrees. Besides, stable isotopes serve as an effective
115‐
tool for revealing the source and fate of aquaculture induced organic matters in the
116‐
environment.
RI PT
113‐
117‐
2. Materials and methods
119‐
2.1 Study area
SC
118‐
This study was conducted in the north part (37°50′ – 38°10′ N, 118°20′ – 119°00′ E) of
121‐
the Yellow River estuary (Fig. 1b) which is located in the northeast of Shandong Province,
122‐
the People’s Republic of China (Fig. 1a). The study site is characterized by a
123‐
warm-temperate continental monsoonal climate, with mean annual temperature and annual
124‐
mean precipitation of 11.9 oC and 640 mm, respectively (Gao et al., 2014).
M AN U
120‐
TE D
A total area of 66.67 km2 was constructed for sea cucumber ponds along the north coast
125‐
of Yellow River estuary. Each pond was constructed with 3 to 6 ha in size and 1.5 to 3m in
127‐
depth. Apostichopus japonicus are cultured in ponds with a density of 10 individuals m-2.
128‐
Commercial feed (fish meal 15%, soybean meal 20%, peanut flour 10% and seaweed meal
129‐
40%) is added to ponds 1-2 times daily. Water exchange is managed uniformly twice a
130‐
month when spring tide occurs. Pond influent water is taken from one tidal creek using
131‐
pumping systems. Pond effluents are released into unified drainage channels then
132‐
discharged into adjacent wetland directly through a sluice gate.
AC C
EP
126‐
133‐
The tidal flat (C-tidal flat) adjacent to the sea cucumber pond area is separated from sea
134‐
cucumber ponds by a 7 m width road. It is dominated by salt marsh plant Suaeda salsa and
6‐ ‐
ACCEPTED MANUSCRIPT
two macrobenthic species, Helice tridens (omnivore) and Macrophthalmus abbreviates
136‐
(omnivore) which are widely distributed in this locality. In order to eliminate potential
137‐
effects of farming activities, we investigated the environments of a tidal flat (N-tidal flat) in
138‐
Natural reserve area (non-farming area) which is 15 km departing from the C-tidal flat and
139‐
was strictly considered with similar natural characteristics (such as sediment type, tidal
140‐
characteristics, elevation).
141‐
2.2 Field sampling and laboratory procedures
M AN U
142‐
SC
RI PT
135‐
Field sampling was conducted in the sea cucumber pond area, C-tidal flat and N-tidal
143‐
flat from May 1 to May 21, 2016. In the C-tidal flat, there were no other allochthonous
145‐
nutrients loading except aquaculture effluent. In the sea cucumber pond area, two stations
146‐
(Fig. 1c) were positioned in the sea cucumber pond and effluent channel, respectively. Two
147‐
gradient-based sampling were set up in C-tidal flat and N-tidal flat. Six stations (Fig. 1c) in
148‐
the C-tidal flat creek were set at 0 m, 50 m, 100 m, 200 m, 500 m and 800 m distances from
149‐
the sluice gate using an electronic distance meter (Nikon COOLSHOT 40i Laser
150‐
Rangefinders) and a GPS receiver (Garmin GPS map64s). Four stations (Fig. 1d) in the
151‐
N-tidal flat creek were set at 0 m, 50 m, 100 m and 200 m distances from the reclaimed
152‐
land. Other two stations (CM and NM) were set at the coastal area of C-tidal flat and
153‐
N-tidal flat (Fig. 1c, d) for water analysis (organic analysis only).
AC C
EP
TE D
144‐
154‐
At each station (except for CM and NM stations), both water and sediment samples
155‐
were collected for analyses of TOC, TN, C/N ratio, δ13C and δ15N, and dissolved inorganic
156‐
nutrients (NO3-, NO2-, NH4+ and PO43-/ AP). The chemical items analyzed in this study
7‐ ‐
ACCEPTED MANUSCRIPT
were shown with their abbreviations in Table 1. Dominant consumers and primary
158‐
producers were sampled at 0 m, 50 m and 100 m stations in each tidal flat. All of these
159‐
samplings were conducted in the ebb tides.
RI PT
157‐
At each station, three replicate water samples were collected by a bucket from the
160‐
surface layer of the water for the analyses of particulate organic matter (POM) and
162‐
dissolved inorganic nutrients. Besides, marine POM (POM(marine)) in CM and NM stations
163‐
were collected at the sea surface of each adjacent offshore. POM samples were collected by
164‐
filtering 0.4 to 3 L surface seawater through a 125 µm mesh (to remove zooplankton)
165‐
followed by glass-fiber filters (Whatman GF/F, φ47 mm). To remove carbonates, POM
166‐
samples were treated with 1.2 N HCl overnight and then oven dried at 60 oC for stable
167‐
isotope analysis. The remained filtered water samples were treated immediately for nutrient
168‐
analysis using SAN++ continuous flow analyzer (SKALAR, Netherlands).
TE D
M AN U
SC
161‐
To detect the biogeochemical characteristics of deposited sediment, three replicate
170‐
surface layer (0-1 cm) sediments were collected by using a 48 mm (inner diameter) acrylic
171‐
tube. An aliquot of the samples was used for the analysis of acid volatile sulphide (AVS)
172‐
concentration to present sediment condition affected by aquaculture (Yokoyama, 2010).
173‐
The remainder was divided into two groups (group A and group B). Group A was used for
174‐
elemental (TOC & TN) and stable isotope analysis (δ13C & δ15N), group B was used for the
175‐
analyses of dissolved inorganic nutrients in sediment. Sediment samples of group A were
176‐
weighed, dried at 60 oC for 72 h,then weighed again to determine the water content, and
177‐
ground to a fine powder. Prior to analyses, samples for carbon isotope analysis were soaked
178‐
in 1.2 N HCl overnight to remove carbonates, filtered on a Nuclepore polycarbonate
AC C
EP
169‐
8‐ ‐
ACCEPTED MANUSCRIPT
track-etch membrane filter (pore size = 0.2 µm) and dried again. Aliquots of samples not
180‐
soaked in 1.2 N HCl were used for nitrogen isotopic analysis. Group B samples used for
181‐
dissolved inorganic nutrients analysis were air-dried at room temperature and passed
182‐
through a 2 mm mesh sieve. Samples for NO3- and NH4+ analysis (5.0 g) were extracted
183‐
with 2 mol l-1KCL (25 ml), and samples for AP (1.2 g) were extracted with 0.5 mol l-1
184‐
NaHCO3 (25 ml) for determining the exchangeable nitrogen and phosphorous. After shock
185‐
extraction for 30 min, the samples were filtered and analyzed by SAN++ continuous flow
186‐
analyzer to determine the concentrations of NO3-, NH4+ and AP.
M AN U
SC
RI PT
179‐
Benthic microalgae (BMA) samples (n = 3) were extracted from the surface sediment in
187‐
each tidal flat (Couch, 1989; Riera and Richard, 1996). The sediment was brought back to
189‐
the laboratory within 30 min and then spread on a flat tray (1 cm depth) and covered by a
190‐
nylon screen (63 µm mesh). A 4 to 5 mm layer of combusted sand (acid-washed and burned,
191‐
125-500 µm) was spread over the nylon screen. The sand was kept wet using filtered
192‐
seawater and then cultured under fluorescent light for 12 to 15 h. The top 2 mm of sand was
193‐
then removed and filtered through a 63 µm fiber sieve by washing with distilled water.
194‐
Water containing the microalgae was filtered onto GF/F (φ25 mm) filters, and oven dried at
195‐
60 oC for stable isotope analysis.
AC C
EP
TE D
188‐
196‐
The dominant salt marsh plant Suaeda salsa was collected by hand. Leaves from 10
197‐
plants were taken as a single sample. Plant debris was collected by hand nets. The dominant
198‐
macrobenthic animals (two crabs) were caught by a long-handled scoop around the
199‐
sampling station within 5 m distance. Each sample was then preserved on ice and returned
200‐
to the laboratory for sorting and sampling. Muscle tissue of crab was taken for animal
9‐ ‐
ACCEPTED MANUSCRIPT
isotope samples. All animal and plant samples were first treated with 1.2 N HCl to remove
202‐
carbonates and then oven dried at 60 oC and ground to fine powder for isotope analysis.
203‐
The elemental and isotope analysis of samples were analyzed using a Delta XP plus
204‐
isotope ratio mass spectrometer linked to a Conflo IIII interface (Thermo Fisher Scientific,
205‐
USA) with a Flash EA1112 elemental analyzer
SC
Stable isotope ratios were expressed in the typical estuary notation (δ13C or δ15N),
206‐ 207‐
RI PT
201‐
defined by the following equation:
δ13C or δ15N = (Rsample/Rstandard – 1) ×103
209‐
and R = 13C/12C or 15N/14N
M AN U
208‐
The standard reference materials were Pee Dee Belemnite (PDB) and atmospheric N2
210‐ 211‐
for carbon and nitrogen isotopes, respectively.
213‐
2.3 Data analysis
TE D
212‐
Principle components analysis (PCA) (Manly, 1986) was conducted in R V3.3.1 to
215‐
explore the relationships among stations that driven by sediment and water variables using
216‐
the function prcomp and princomp of R package stats. To evaluate spatial differences in
217‐
N-tidal flat, C-tidal flat, effluent channel and sea cucumber pond, a multivariate technique
218‐
was applied using the PRIMER 6.0 ecological software package developed by the
219‐
Plymouth Marine Laboratory. Ordination was computed via hierarchical cluster analysis
220‐
(CLUSTER) based on the Euclidean distance of studied environmental variables and
221‐
isotopic compositions. Seriation test and SIMPER (Similarity percentages) analysis were
222‐
performed to analyze the community structure. One-way ANOSIM was performed to test
AC C
EP
214‐
10‐ ‐
ACCEPTED MANUSCRIPT
223‐
the difference in the studied areas. Multiple comparisons of each environmental variable between stations were conducted
225‐
to detect spatial differences between farming and non-farming areas or among stations with
226‐
a gradient distance to the cultured pond. As no gradient change was observed in the N-tidal
227‐
flat along the distance (0 m, 50 m, 100 m, 200 m), the data from different sampling stations
228‐
were calculated together to represent a reference location expressed as N-tidal flat.
SC
RI PT
224‐
Bonferroni test (Rice, 1989) was used when data is normally distributed (Shapiro and
229‐
Wilk, 1965) and homoscedastic (Bartlett, 1937). While Kruskal-Wallis test (Kruskal and
231‐
Wallis, 1952) was performed when the data is not normally distributed or homoscedastic,
232‐
followed by Gao post hoc test (Gao et al., 2008; Konietschke et al., 2015) to compare
233‐
differences among stations for the selected environmental variable. Statistical significance
234‐
level was considered at 5%. Statistical analyses were performed using the function
235‐
Wilk-Shapiro test, Bartlett’s test, Bonferroni test and Kruskal test of the R package stats
236‐
and the function Gao of the R package nparcomp (Ihaka and Gentleman, 1996; Shirley,
237‐
1977) and oneway_test function of the R package coin (Ihaka and Gentleman, 1996) in R
238‐
V3.3.1.
EP
TE D
M AN U
230‐
Carbon and nitrogen stable isotopes of primary producers (food sources) and consumers
AC C
239‐ 240‐
were initially investigated using biplot diagrams. Because no differences in stable isotope
241‐
composition of consumers and primary producers were found along the distance gradient in
242‐
N-tidal flat and C-tidal flat, each of them was calculated together in each tidal flat. A t-test
243‐
was performed to test significant differences in stable isotope values of consumers and
244‐
primary producers. The Bayesian stable isotope mixing model, SIAR v 4.0 (Stable Isotope
11‐ ‐
ACCEPTED MANUSCRIPT
Analysis in R) (Parnell et al., 2010) was used to determine proportional contributions of
246‐
potential food sources to dominant consumers. Contributions of dietary sources were
247‐
represented as mean and 95% credibility interval (Ci). The fractionations δ13C and δ15N for
248‐
crustacean (muscle, acid treated) calculated from raw data presented in Yokoyama et al.
249‐
(2005) were δ13C = 1.93 ± 0.40‰, δ15N = 3.97 ± 0.49‰, respectively.
RI PT
245‐
3. Results
252‐
3.1 Environmental variables
M AN U
251‐
SC
250‐
Each sampling location (sea cucumber pond, effluent channel, C-tidal flat and N-tidal
253‐
flat) were clearly separated from each other by PCA analysis (Fig. 2a) with seventeen water
255‐
and sediment variables (Table 1). The first three principal components (PC1, PC2 and PC3)
256‐
explained 73% of the total variance. The main contributions to PC1 included sediment
257‐
variables TOC(S), TN(S), CN(S), δ13C(S), NO3-(S), NH4+ (S), PO43- (S) (with large positive
258‐
coefficients), δ15N(S) and δ13C(POM) (with negative coefficient). Vector plots of sediment
259‐
variables and δ13C(POM) showed that PC1 represented axes of increasing organic
260‐
accumulation and nutrients concentration in sediment and decreasing δ15N(S) and δ13C(POM)
261‐
from tidal flat to effluent channel and then followed by sea cucumber pond. As for PC2
262‐
(Fig. 2b), POM variables (except for δ13C(POM)) and PO43-(water) in water were closely related
263‐
to PC2 as indicated by large positive coefficients (TOC(POM), TN(POM), CN(POM), δ15N(POM))
264‐
and negative coefficient (PO43-(water)), respectively. PC2 represented axes of increasing
265‐
organic accumulation in POM and decreasing PO43-(water) concentration from sea cucumber
266‐
pond area to N-tidal flat. The main contribution to PC3 (Fig. 2c) was from DIN(water) that
AC C
EP
TE D
254‐
12‐ ‐
ACCEPTED MANUSCRIPT
showed highly positive coefficients (NO3-(water), NO2-(water), NH4+(water)). As the vector plots
268‐
of DIN(water) that PC3 represented, however, no obvious gradient variation was found across
269‐
stations.
RI PT
267‐
The cluster analysis (Fig.3) showed that studied localities could be distinguished into
271‐
three groups (sea cucumber pond, the combination of C-tidal flat and effluent channel,
272‐
N-tidal flat) with the threshold Euclidean distance more than 7.0 (Fig. 3). The results of the
273‐
pairwise test in ANOSIM showed that significant differences were present between N tidal
274‐
flat and sea cucumber pond (R = 1, P < 0.01), N tidal flat and the combination of C tidal
275‐
flat and effluent channel (R = 0.529, P < 0.01), sea cucumber pond and the combination of
276‐
C tidal flat and effluent channel (R = 0.991, P < 0.01). With the subsequent SIMPER
277‐
analysis, the primary contributors to three groups were TOC(s) (89.60% contribution),
278‐
TOC(water) (63.16% contribution) and TOC(s) (97.82% contribution).
TE D
M AN U
SC
270‐
279‐
3.2 Water and sediment biogeochemical characteristics
281‐
3.2.1 Water biogeochemical characteristics
EP
280‐
TOC(POM) and TN(POM) concentrations (Figs. 4a, b) detected in farming area (mean ± SD
282‐
= 0.07 ± 0.02 and 0.05 ± 0.2 mg l-1, respectively) were significantly lower than those in
284‐
non-farming area (1.4 ± 0.6 and 0.2 ± 0.1 mg l-1, respectively). C/N ratios of POM (Fig. 4c)
285‐
ranged from 4.7 to 9.1 and showed significant differences between sea cucumber pond and
286‐
N-tidal flat (F = 9.350, P < 0.001). δ13C(POM) in sea cucumber pond (-20.6 ± 1.0‰),
287‐
effluent channel (-19.9 ± 1.2‰) and C-marine stations (-21.1 ± 0.1‰) showed significantly
288‐
lower values than those in C-tidal flat (-16.7 ± 1.6‰), N-tidal flat (-16.9 ± 1.2‰) and
AC C
283‐
13‐ ‐
ACCEPTED MANUSCRIPT
N-marine (-18.4 ± 0.4‰) stations (Fig. 4d). δ15N(POM) values in the sea cucumber pond
290‐
(-2.0 ± 0.2‰) was significantly lower than those in other stations (F = 38.06, P < 0.001).
291‐
No remarkable differences were found within C-tidal flat stations (0.9 ± 0.6‰).
RI PT
289‐
NO3- and NH4+ showed no significant differences (F = 1.90, P > 0.05; F = 0.62, P >
293‐
0.05) among stations (Figs. 4f, h), whereas significant differences were present in NO2- and
294‐
PO43- (X² = 21.85, DF = 8, P < 0.01; X² = 28, DF = 8, P < 0.001) (Fig. 4g, i). The NO2-
295‐
concentration in sea cucumber pond (53 ± 2 µg l-1) was clearly higher than that in N-tidal
296‐
flat (40 ± 9 µg l-1), but no distinct differences were found among stations in the farming
297‐
area or between C-tidal flat and N-tidal flat. The PO43- concentrations in sea cucumber
298‐
pond (32 ± 1 µg l-1), effluent channel (33 ± 1 µg l-1) and C-tidal flat (29 ± 3 µg l-1) were
299‐
clearly higher than that in N-tidal flat (20 ± 2 µg l-1), but no distinct differences were found
300‐
within stations in the farming area.
TE D
301‐ 302‐
M AN U
SC
292‐
3.2.2 Sediment biogeochemical characteristics AVS concentration detected in sea cucumber pond was 0.1 ± 0.5 mg g-1 (dry sediment)
304‐
higher than those in the effluent channel, C-tidal flat and N-tidal flat (all less than 0.01 mg
305‐
g-1). For TOC(S) and TN(S), significant differences (X² = 25.61, DF = 8, P < 0.01; X² =
306‐
31.246, DF = 8, P < 0.001) were observed among stations (Figs. 5a, b). TOC(S) and TN(S) in
307‐
the sea cucumber pond (5.9 ± 1 mg g-1 and 0.8 ± 0.3 mg g-1, respectively) and effluent
308‐
channel (2.95 ± 0.43 mg g-1 and 0.51 ± 0.08 mg g-1, respectively) were significantly higher
309‐
than those in the C-tidal flat (1.2 ± 0.3 mg g-1 and 0.3 ± 0.1 mg g-1, respectively) and
310‐
N-tidal flat (1.8 ± 0.2 mg g-1 and 0.3 ± 0 mg g-1, respectively). Both TOC(S) and TN(S)
AC C
EP
303‐
14‐ ‐
ACCEPTED MANUSCRIPT
311‐
contents showed decreasing trends from sea cucumber pond to 0 m station, and no gradient
312‐
variations along the dispersal distance in C-tidal flat. C/N ratios (Fig. 5c) in the pond, effluent channel, 0 m station and N-tidal flat had
314‐
significant differences from other stations in C-tidal flat (50 m, 100 m, 200 m and 800 m).
315‐
A decreasing trend was found from pond to 0 m. δ13C(S) and δ15N(S) values of sea cucumber
316‐
pond were remarkably different from the effluent channel, C-tidal flat and N-tidal flat
317‐
stations (Fig. 5d, e). A decreasing trend was found from pond to 0 m station. δ13C(S) of 0 m
318‐
station showed a large range of standard deviation and no significant differences with the
319‐
effluent channel and other stations in C-tidal flat and N-tidal.
M AN U
SC
RI PT
313‐
No significant differences (Fig. 5f) among stations were found for NO3-(S) (X² = 11.982,
320‐
DF = 8, P > 0.05). Both NO3-(S) and NH4+(S) were found a decreasing trend from pond to 0
322‐
m station. The concentration of NH4+(S) ranged from 1.8 to 18.1 mg kg-1 (Fig. 5g) and
323‐
varied across stations (X² = 27.671, DF = 8, P < 0.001). The highest NH4+(S) concentration
324‐
found in sea cucumber pond (15.8 ± 2.0 mg kg-1) was significantly higher than those from
325‐
effluent channel (8.6 ± 2.9 mg kg-1), C-tidal flat (5.7 ± 2.6 mg kg-1) and N-tidal flat stations
326‐
(4.3 ± 0.5 mg kg-1). NH4+(S) content in effluent channel showed no significant differences
327‐
with C-tidal flat stations (except for 500m), but was obviously higher than that in N-tidal
328‐
flat. A gradient decrease was observed from pond to 0 m station. AP contents (Fig. 5h) in
329‐
pond (8.2 ± 0.9 mg kg-1) and effluent channel (10.4 ± 3.9 mg kg-1) were significantly higher
330‐
than those in both C-tidal flat (3.6 ± 0.6 mg kg-1) and N-tidal flat (3.4 ± 2.1 mg kg-1)
331‐
stations (F = 7.00, P < 0.001).
AC C
EP
TE D
321‐
332‐
15‐ ‐
ACCEPTED MANUSCRIPT
333‐
3.3 Feeding strategies of dominant macrobhenthos δ13C value of H. tridens in the C-tidal flat (-13 ± 0.6‰) was significantly higher (t-test,
334‐
P < 0.01) than that in the N-tidal flat (-15.5 ± 3.3‰) (Table. 2). δ15N value of H. tridens in
336‐
C-tidal flat (4.3 ± 1.1‰) was significant than that in N-tidal flat (5.8 ± 0.4‰). As for the
337‐
crab M. abbreviates, δ13C values were similar (t-test, P > 0.05) between C-tidal flat (-12.1 ±
338‐
0.6‰) and N-tidal flat (-12.3 ± 0.4‰), whereas δ15N values in C-tidal flat (2.1 ± 0.7‰)
339‐
was significantly lower (t-test, P < 0.01) compared to that in N-tidal flat (3.8 ± 0.9‰).
SC
RI PT
335‐
M AN U
In according with the presumptive trophic shift of crustaceans (1.9 ± 0.4‰ in δ13C and
340‐
4.0 ± 0.5‰ in δ15N) presented by Yokoyama et al. (2005), the expected δ13C and δ15N
342‐
values of the H. tridens diet would be -17.4 ± 3.7‰ and 1.8 ± 0.9‰, respectively in N-tidal
343‐
flat, and -14.9 ± 1.0‰ and 0.3 ± 1.6‰, respectively in C-tidal flat. The expected δ13C and
344‐
δ15N values of M. abbreviates diet would be -17.4 ± 3.7‰ and 1.8 ± 0.9‰, respectively in
345‐
N-tidal flat, and -14.9 ± 1.0‰ and 0.3 ± 1.6‰, respectively in C-tidal flat.
TE D
341‐
Among the primary producers, salt marsh plant Suaeda salsa had the lowest δ13C
347‐
(around -29.7 ± 1.3 to -30.4 ± 0.5‰ ) and highest δ15N (around 6.1 ± 0.9 to 6.3 ± 1.5‰)
348‐
values. There were large differences in isotope compositions between Suaeda salsa and
349‐
hypothetical diet of the two crabs. The δ13C values of hypothetical diet for two crabs were
350‐
between POM (POM(marine), POM(tidal flat) and POM(pond)) and BMA (Fig. 6).
AC C
EP
346‐
SIAR analysis result (Fig. 7a, b) revealed that in N-tidal flat, the contributions of
351‐ 352‐
POM(marine) to the diet of H. tridens (25%, Ci 0 – 49%) was similar to those of BMA (27%,
353‐
Ci 8 – 45%) and plant debris (25%, Ci 8 – 42%) (Fig. 7a), while POM(tidal flat) (19%, Ci 0 –
354‐
37%) and salt marsh plant Suaeda salsa (4%, Ci 0 – 9%) made minor contributions. In
16‐ ‐
ACCEPTED MANUSCRIPT
C-tidal flat, POM(marine) made a smaller contribution to the diet of H. tridens (3%, Ci 0 –
356‐
7%), while BMA made a greater contribution (45%, Ci, 36 – 54%) (Fig. 7b). Besides, the
357‐
contributions of POM(tidal flat) (15%, Ci 0 – 32%), plant debris (18%, Ci 0 – 36%) and
358‐
Suaeda salsa (1%, Ci 0 – 2%) were relatively low. In addition, POM(pond) contributed 19%
359‐
(Ci 2 – 34%) of the diet of H. tridens in C-tidal flat.
RI PT
355‐
As for crab M. abbreviates (Fig. 7c, d), BMA overwhelmingly dominated the diet in
SC
360‐
N-tidal flat (51%, Ci 39 – 63%) and C-tidal falt (54%, Ci 37 – 69%). But same with H.
362‐
tridens, POM(marine) showed the much lower contribution to the diet of M. abbreviates in
363‐
C-tidal flat (6%, Ci 0 – 16%) than that in N-tidal flat(22%, Ci 0 – 44%). Likewise, POM(tidal
364‐
flat)
365‐
contribution to the diet of M. abbreviates in both N-tidal flat and C-tidal falt. But, POM(pond)
366‐
contributed a relatively lower amount (9%, Ci 0 – 22%) to the diet of M. abbreviates than it
367‐
contributed to H. tridens in C-tidal flat.
M AN U
361‐
368‐
TE D
and plant debris made small contributions, and Suaeda salsa made a minimal
4. .Discussion
370‐
4.1 Differences of distance-based environmental variables and their implications
AC C
EP
369‐
Overall, this study intended to evaluate the environmental impacts of sea cucumber
371‐ 372‐
farming based on two kinds of analyses. The first one was the detection of discrepancies
373‐
among distance-based pollution indicators in the water column and sediment. The second
374‐
one was the utilization of biochemical tracers (δ13C and δ15N) to trace the origin of organic
375‐
matters. Both methods had been widely confirmed as reliable tools to evaluate
376‐
environmental effects of aquaculture (Gowen and Bradbury, 1987; Sarà, 2007; Sarà et al.,
17‐ ‐
ACCEPTED MANUSCRIPT
2004; Wai et al., 2011). A total of eight typical environmental variables (TOC, TN, C/N,
378‐
AVS, NO3-, NH4+, NO2-, PO43- and AP) were studied in this paper. For the studied
379‐
parameters in the water column, they showed no uniform variation tendencies. For example,
380‐
TOC(POM), TN(POM), NH4+(water) and NO3-(water) contents were similar between farming and
381‐
non-farming areas, or among stations within the farming area (such as pond and its adjacent
382‐
tidal flat). NO2-(water) and PO43-(water) contents in the farming area were higher than those in
383‐
N-tidal flat, while they did not differ among farming areas (Fig. 4). In our study,
384‐
distribution patterns of studied parameters may have many potential reasons. One possible
385‐
reason lies in that eutrophication, related to the rapid development of mariculture
386‐
(Bouwman et al., 2013) or other anthropogenic activities, are widely present in the Yellow
387‐
River estuary and its adjacent waters (Chen et al., 2013; Kang and Xu, 2016; Zhang et al.,
388‐
2008). High background values hinder the statistical differences.
TE D
M AN U
SC
RI PT
377‐
In pervious studies, NO3-(water) was recognized as the typical dominant part of DIN in
390‐
the studied coastal waters (Gong, 2012; Zhang et al., 2008). But our study found the same
391‐
result as Kang (2016), that NH4+(water) represented the relative major proportion of DIN than
392‐
NO3-(water) in the mariculture zone. This result corresponded well to the meta-analysis of
393‐
Sarà (2007) that NH4+(water) appears to be the most affected compound by aquaculture
394‐
loadings than NO2-(water), NO3-(water) and followed by PO43-(water). However, discrepancies
395‐
among distribution patterns of studied parameters prevented us from drawing much more
396‐
clear conclusions.
AC C
EP
389‐
Despite the distinct distribution patterns of nutrients in the water column, NO3-(S),
397‐ 398‐
NH4+(S), AP(S), TOC(S), TN(S) and AVS contents in sediment showed a uniform
18‐ ‐
ACCEPTED MANUSCRIPT
distance-based differences. They had significantly higher values in farming area than those
400‐
in non-farming areas. Within the farming area, a decrease transect from pond to 0 m station
401‐
was detected. These findings clearly demonstrated the presence of aquaculture waste from
402‐
the sea cucumber farming, and additionally implied a restricted dispersion distance
403‐
(approximately 50 m from farming pond to 0 m station in C-tidal flat) of aquaculture waste
404‐
that accumulated in sediment. It was also worth noting that although AVS content in sea
405‐
cucumber pond in this study was below the threshold (0.2 mg/g) in some rules ensuring a
406‐
sustainable aquaculture (Uede, 2008; Yokoyama, 2003b; Yokoyama et al., 2007), the
407‐
content was much higher than those reported from other natural sediments. Similar cases
408‐
were also found in TOC(S) and TN(S) concentration in sea cucumber pond. This phenomenon
409‐
were also considered as results of sea cucumber aquaculture.
411‐
4.2 Detection of origins of organic matters and their implications
TE D
410‐
M AN U
SC
RI PT
399‐
In the water column, the result of C/N ratios(water) and isotopic composition analysis
412‐
revealed that sea cucumber pond had different organic origins from adjacent tidal flat and
414‐
non-farming area. For example, C/N ratios(POM) detected in pond, 0 m, 50 m, 100 m and
415‐
500 m stations were close to the Redfield ratio (C/N=106:16) for marine plankton
416‐
(Andrews et al., 1998a; Burkhardt et al., 1999; Redfield, 1963), which might indicate that
417‐
the main contributor of POM in these stations was marine phytoplankton . The relatively
418‐
higher values of C/N ratios(POM) in the effluent channel, 200 m, 800 m and N-tidal flat
419‐
suggested that POM resulted from a mixture of marine phytoplankton and other organic
420‐
matter sources. Low C/N ratios(POM) detected in N-marine and C-marine areas might be
AC C
EP
413‐
19‐ ‐
ACCEPTED MANUSCRIPT
attributed to the contribution of marine zooplankton with low C/N ratio (Båmstedt, 1986;
422‐
Walve and Larsson, 1999). The lowest δ13C(POM) and δ15N(POM) values found in sea
423‐
cucumber pond maybe due to the allogenous sources. Terrestrial plant (soybean and peanut)
424‐
and seaweed flour (seaweed cultured in offshore) introduced by aquaculture (commercial
425‐
feeds) had low δ13C and δ15N values, and the presence of allogenous sources resulted in low
426‐
isotopic values. In this case, the distribution pattern of δ13C(POM) and δ15N(POM) values
427‐
indicated the environmental impact of aquaculture.
SC
RI PT
421‐
M AN U
In the sediments, high δ13C(S) and low δ15N(S) values were detected in sea cucumber
428‐
pond. Similar to the reason of low δ15N(POM) compositions in the water column, low δ15N(S)
430‐
values should be due to the allogenous organic matter input. δ13C(S) in pond sediment was
431‐
approximately equal to δ13C(POM) of POM , which suggested that they had homogenous
432‐
resources. In the adjacent and natural tidal flat, deposition of organic matter from salt marsh
433‐
plant (δ13C = 29.74 ± 1.31, δ15N = 6.1 ± 0.91) contributed to the sediment with much lower
434‐
δ13C and high δ15N values. Therefore, high δ13C and low δ15N values in the pond were a
435‐
reflection of the disturbance of aquaculture waste.
EP
TE D
429‐
436‐
4.3 Effects of pond farming waste on food sources for dominant macrobenthos
AC C
437‐
SIAR analysis results revealed that BMA was the most important food resources for
438‐ 439‐
H. tridens and M. abbreviates in both farming and non-farming area. However, for both
440‐
species, contributions of BMA to the diet in the C-tidal flat were bigger than those in
441‐
N-tidal flat. This founding might imply that nutrient enrichment of surface sediment may
442‐
enhance the primary production of BMA which furthermore enabled a higher contribution
20‐ ‐
ACCEPTED MANUSCRIPT
to the consumer's diet. Meanwhile, the contribution of POM(marine) in C-tidal flats was lower
444‐
than that in N-tidal flats and the POM(pond) showed opposite tendencies. This differences
445‐
indicated that sea cucumber farming activities might modify the dietary characteristics of
446‐
benthic animals to a certain extent. The similar modifications of food composition of
447‐
dominant benthic animals by aquaculture were confirmed in previous studies of other
448‐
aquaculture type, such as fish farming (Mazzola and Sarà, 2001; Yokoyama and Ishihi,
449‐
2007) and shrimp pond farming (Kon et al., 2009).
M AN U
450‐ 451‐
SC
RI PT
443‐
4.4 Implications for operation and management of sea cucumber farming The overall conclusion of the present study was that sea cucumber farming slightly
452‐
altered the local environment, however, this modification was restricted to a limited area
454‐
(within 50 m). Up to present, there existed a relative paucity of details concerning the
455‐
environmental effects of sea cucumber farming (Li et al., 2014; Ren et al., 2010; Zheng et
456‐
al., 2009). Therefore, it’s hard to make a comparison of our results with the previous
457‐
studies. In this study, the short dispersal distance might be attributed to the ideal operation
458‐
mode and management of sea cucumber pond farming in the studied region. Firstly, the
459‐
main feed for Apostichopus japonicus is composed of much seaweed flour (40%) and less
460‐
protein feeds (fishmeal = 15%) in comparison with the shrimp feed (around 45%) (Dierberg
461‐
and Kiattisimkul, 1996; Paez-Osuna, 2001; Trott et al., 2004), that reduced the
462‐
accumulation of organic matter from the waste feed. Besides, the intense of sea cucumber
463‐
farming activities was moderate, with low densities (10 individuals m-2) that far below the
464‐
typical densities (20 individuals m-2) practiced by Chinese farmers (Dong et al., 2010; Kang
AC C
EP
TE D
453‐
21‐ ‐
ACCEPTED MANUSCRIPT
and Xu, 2016). In addition, waste waters from sea cucumber ponds, connected by internal
466‐
pipes, were gathered in large channels and drained intermittently by valves, which may
467‐
result in a higher deposition before the discharge. The modification in sediment was
468‐
restricted to 50 m distance, which indicates that farming pond and effluent channel were
469‐
heavily disturbed. Consequently, high natural feeds content, low stocking densities, and line
470‐
effluent channels reduced the deleterious effects of sea cucumber farming on the adjacent
471‐
environment. Besides, deposited sludge should be more effiently removed physical efforts
472‐
or by filters like salt marsh plant, bivalves, microalgae.
M AN U
SC
RI PT
465‐
For the management of aquaculture activities, one of the challenge is to develop
474‐
practical, cost-effective monitoring techniques that reflect the impact and processes
475‐
(Burford et al., 2003). The present study provides meaningful information about
476‐
environmental impacts of Apostichopus japonicus farming. The unique characteristics of
477‐
Apostichopus japonicus farming discharges indicated that the predominant effects of sea
478‐
cucumber farming discharge were in the sediment processes, rather than water column
479‐
processes. In this case, environmental variables should be unequally treated.
EP
TE D
473‐
In conclusion, the present study considered that pond farming for Apostichopus
481‐
japonicus in the Yellow River estuary altered local environment to a certain extent, as
482‐
determined by the discrepancies among distance-based environmental parameters in
483‐
sediment (TOC(S), TN(S) and AVS, NO3-(S), NH4+(S), AP(S)) as well as isotopic signatures
484‐
(δ13C and δ15N). Such kind of modification was restricted to 50 m distance (covers farming
485‐
pond and effluent channel). Apostichopus japonicus farming also modified dietary
486‐
characteristics of some dominant macroinvertebrates. For methodological consideration,
AC C
480‐
22‐ ‐
ACCEPTED MANUSCRIPT
sediment biogeochemical characteristics as a historical recorder much more effectively
488‐
reflected both organic and dissolved inorganic nutrients induced by aquaculture.
RI PT
487‐
489‐ Acknowledgements
491‐
This study was supported financially by National Key Basic Research Programme of China
492‐
(2013CB430406). We acknowledge experimental support by Center for Ecological Research of Kyoto
493‐
University. We are grateful to the great support from the lab of Aquatic Environmental Biology, Kyoto
494‐
University. We thank Dr. Edouard Lavergne for the idea of statistical analyses in environment variables.
495‐
And thank Prof. Yoh Yamashita, Dr. Bunei Nishimura and Dr. Zou Yuxuan for their assistance in
496‐
different aspects of this study.
498‐ 499‐
M AN U
497‐
SC
490‐
References
501‐ 502‐ 503‐ 504‐ 505‐ 506‐ 507‐ 508‐ 509‐ 510‐ 511‐ 512‐ 513‐ 514‐ 515‐ 516‐ 517‐ 518‐ 519‐ 520‐ 521‐ 522‐ 523‐ 524‐ 525‐ 526‐ 527‐ 528‐
Andrews, J., Greenaway, A., Dennis, P., 1998a. Combined carbon isotope and C/N ratios as indicators of source and fate of organic matter in a poorly flushed, tropical estuary: Hunts Bay, Kingston Harbour, Jamaica. Estuarine, Coastal and Shelf Science 46, 743-756. Andrews, J.E., Greenaway, A.M., Dennis, P.F., 1998b. Combined Carbon Isotope and C/N Ratios as Indicators of Source and Fate of Organic Matter in a Poorly Flushed, Tropical Estuary: Hunts Bay, Kingston Harbour, Jamaica. Estuarine, Coastal and Shelf Science 46, 743-756. Båmstedt, U., 1986. Chemical composition and energy content. The biological chemistry of marine copepods. Clarendon Press, Oxford, 1-58. Bartlett, M.S., 1937. Properties of sufficiency and statistical tests. Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 268-282. Bordbar, S., Anwar, F., Saari, N., 2011. High-value components and bioactives from sea cucumbers for functional foods—a review. Marine drugs 9, 1761-1805. Bouwman, L., Beusen, A., Glibert, P.M., Overbeek, C., Pawlowski, M., Herrera, J., Mulsow, S., Yu, R., Zhou, M., 2013. Mariculture: significant and expanding cause of coastal nutrient enrichment. Environmental Research Letters 8, 044026. Brooks, K.M., Mahnken, C.V., 2003. Interactions of Atlantic salmon in the Pacific northwest environment: II. Organic wastes. Fisheries Research 62, 255-293. Burford, M.A., Costanzo, S.D., Dennison, W.C., Jackson, C.J., Jones, A.B., McKinnon, A.D., Preston, N.P., Trott, L.A., 2003. A synthesis of dominant ecological processes in intensive shrimp ponds and adjacent coastal environments in NE Australia. Marine Pollution Bulletin 46, 1456-1469. Burkhardt, S., Zondervan, I., Riebesell, U., 1999. Effect of CO2 concentration on C: N: P ratio in marine phytoplankton: A species comparison. Limnology and Oceanography 44, 683-690. Chen, J., 2003. Overview of sea cucumber farming and sea ranching practices in China. SPC beche-de-mer Information Bulletin 18, 18-23. Chen, P., Liu, S., Zhang, G., Li, L., Gu, P., Cao, X., 2013. Monthly variation of nutrient concentrations and fluxes in the lower Huanghe River: under the influence of artifical floods. Acta Oceanol Sin 35, 59-71. Couch, C.A., 1989. Carbon and nitrogen stable isotopes of meiobenthos and their food resources. Estuarine, Coastal and Shelf Science 28, 433-441.
AC C
EP
TE D
500‐
23‐ ‐
ACCEPTED MANUSCRIPT
EP
TE D
M AN U
SC
RI PT
Dierberg, F.E., Kiattisimkul, W., 1996. Issues, impacts, and implications of shrimp aquaculture in Thailand. Environmental Management 20, 649-666. Dong, S., Liang, M., Gao, Q., Wang, F., Dong, Y., Tian, X., 2010. Intra specific effects of sea cucumber (Apostichopus japonicus) with reference to stocking density and body size. Aquaculture Research 41, 1170-1178. FAO, 2008. Sea cucumbers, a global review of fisheries and trade. Food and Agriculture Organization, Rome, Italy. Gao, H., Bai, J., He, X., Zhao, Q., Lu, Q., Wang, J., 2014. High temperature and salinity enhance soil nitrogen mineralization in a tidal freshwater marsh. PloS one 9, e95011. Gao, X., Alvo, M., Chen, J., Li, G., 2008. Nonparametric multiple comparison procedures for unbalanced one-way factorial designs. Journal of Statistical Planning and Inference 138, 2574-2591. Gao, X., Yang, Y., Wang, C., 2012. Geochemistry of organic carbon and nitrogen in surface sediments of coastal Bohai Bay inferred from their ratios and stable isotopic signatures. Mar Pollut Bull 64, 1148-1155. Gong, Y., 2012. Influencing factor of nutrient transport in the lower reach of Yellow River. Doctoral Dissertation. Qingdao: Ocean University of China. 1–110. Gowen, R.J., Bradbury, N.B., 1987. THE ECOLOGICAL IMPACT OF SALMONID FARMING IN COASTAL WATERS - A REVIEW. Oceanography and Marine Biology 25, 563-575. Han, Q., Keesing, J.K., Liu, D., 2016. A Review of Sea Cucumber Aquaculture, Ranching, and Stock Enhancement in China. Reviews in Fisheries Science & Aquaculture 24, 326-341. Holby, O., Hall, P., 1991. Chemical fluxes and mass balances in a marine fish cage farm. II. Phosphorus. Mar. Ecol. Progr. Ser. 70: 263-272. Ihaka, R., Gentleman, R., 1996. R: a language for data analysis and graphics. Journal of computational and graphical statistics 5, 299-314. Kang, P., Xu, S., 2016. The impact of mariculture on nutrient dynamics and identification of the nitrate sources in coastal waters. Environmental science and pollution research international 23, 1300-1311. Kollongei, K.J., Lorentz, S.A., 2014. Connectivity influences on nutrient and sediment migration in the Wartburg catchment, KwaZulu-Natal Province, South Africa. Physics and Chemistry of the Earth, Parts A/B/C 67, 12-22. Kon, K., Kawakubo, N., Aoki, J.I., Tongnunui, P., Hayashizaki, K.I., Kurokura, H., 2009. Effect of shrimp farming organic waste on food availability for deposit feeder crabs in a mangrove estuary, based on stable isotope analysis. Fisheries Science 75, 715-722. Konietschke, F., Placzek, M., Schaarschmidt, F., Hothorn, L.A., 2015. nparcomp: An R Software Package for Nonparametric Multiple Comparisons and Simultaneous Confidence Intervals. 2015 64, 17. Kruskal, W.H., Wallis, W.A., 1952. Use of ranks in one-criterion variance analysis. Journal of the American statistical Association 47, 583-621. Li, J., Dong, S., Gao, Q., Wang, F., Tian, X., Zhang, S., 2014. Total organic carbon budget of integrated aquaculture system of sea cucumber Apostichopus japonicus, jellyfish Rhopilema esculenta and shrimp Fenneropenaeus chinensis. Aquaculture Research 45, 1825-1831. Li, Q., Chen, L., Kong, L., 2009. A genetic linkage map of the sea cucumber, Apostichopus japonicus (Selenka), based on AFLP and microsatellite markers. Animal genetics 40, 678-685. Liu, G., Yang, H., Liu, S., 2010. Effects of rearing temperature and density on growth, survival and development of sea cucumber larvae, Apostichopus japonicus (Selenka). Chinese Journal of Oceanology and Limnology 28, 842-848. Liu, Y., Dong, S., Tian, X., Wang, F., Gao, Q., 2009. Effects of dietary sea mud and yellow soil on growth and energy budget of the sea cucumber Apostichopus japonicus (Selenka). Aquaculture 286, 266-270. Lupatsch, I., Kissil, G.W., 1998. Predicting aquaculture waste from gilthead seabream (Sparus aurata) culture using a nutritional approach. Aquatic Living Resources 11, 265-268. Manly, B.F., 1986. Randomization and regression methods for testing for associations with geographical, environmental and biological distances between populations. Researches on Population Ecology 28, 201-218. Mazzola, A., Sarà, G., 2001. The effect of fish farming organic waste on food availability for bivalve molluscs (Gaeta Gulf, Central Tyrrhenian, MED): stable carbon isotopic analysis. Aquaculture 192, 361-379. Naylor, R.L., Goldburg, R.J., Primavera, J.H., Kautsky, N., Beveridge, M.C.M., Clay, J., Folke, C.,
AC C
529‐ 530‐ 531‐ 532‐ 533‐ 534‐ 535‐ 536‐ 537‐ 538‐ 539‐ 540‐ 541‐ 542‐ 543‐ 544‐ 545‐ 546‐ 547‐ 548‐ 549‐ 550‐ 551‐ 552‐ 553‐ 554‐ 555‐ 556‐ 557‐ 558‐ 559‐ 560‐ 561‐ 562‐ 563‐ 564‐ 565‐ 566‐ 567‐ 568‐ 569‐ 570‐ 571‐ 572‐ 573‐ 574‐ 575‐ 576‐ 577‐ 578‐ 579‐ 580‐
24‐ ‐
ACCEPTED MANUSCRIPT
EP
TE D
M AN U
SC
RI PT
Lubchenco, J., Mooney, H., Troell, M., 2000. Effect of aquaculture on world fish supplies. Nature 405, 1017-1024. Paerl, H.W., 2006. Assessing and managing nutrient-enhanced eutrophication in estuarine and coastal waters: Interactive effects of human and climatic perturbations. Ecological Engineering 26, 40-54. Paez-Osuna, F., 2001. The environmental impact of shrimp aquaculture: Causes, effects, and mitigating alternatives. Environmental Management 28, 131-140. Paltzat, D.L., Pearce, C.M., Barnes, P.A., McKinley, R.S., 2008. Growth and production of California sea cucumbers (Parastichopus californicus Stimpson) co-cultured with suspended Pacific oysters (Crassostrea gigas Thunberg). Aquaculture 275, 124-137. Parnell, A.C., Inger, R., Bearhop, S., Jackson, A.L., 2010. Source partitioning using stable isotopes: coping with too much variation. PloS one 5, e9672. Purcell, S.W., Hair, C.A., Mills, D.J., 2012. Sea cucumber culture, farming and sea ranching in the tropics: Progress, problems and opportunities. Aquaculture 368, 68-81. Pusceddu, A., Fraschetti, S., Mirto, S., Holmer, M., Danovaro, R., 2007. Effects of intensive mariculture on sediment biochemistry. Ecological Applications 17, 1366-1378. Redfield, A.C., 1963. The influence of organisms on the composition of sea-water. The sea, 26-77. Ren, Y., Dong, S., Qin, C., Wang, F., Tian, X., Gao, Q., 2012. Ecological effects of co-culturing sea cucumber Apostichopus japonicus (Selenka) with scallop Chlamys farreri in earthen ponds. Chinese journal of Oceanology and Limnology 30, 71-79. Ren, Y., Dong, S., Wang, F., Gao, Q., Tian, X., Liu, F., 2010. Sedimentation and sediment characteristics in sea cucumber Apostichopus japonicus (Selenka) culture ponds. Aquaculture research 42, 14-21. Rice, W.R., 1989. Analyzing tables of statistical tests. Evolution 43, 223-225. Riera, P., Richard, P., 1996. Isotopic determination of food sources of Crassostrea gigas along a trophic gradient in the estuarine bay of Marennes-Oleron. Estuarine Coastal and Shelf Science 42, 347-360. Sakamaki, T., Shum, J.Y.T., Richardson, J.S., 2010. Watershed Effects on Chemical Properties of Sediment and Primary Consumption in Estuarine Tidal Flats: Importance of Watershed Size and Food Selectivity by Macrobenthos. Ecosystems 13, 328-337. Sarà, G., 2007. A meta-analysis on the ecological effects of aquaculture on the water column: Dissolved nutrients. Marine Environmental Research 63, 390-408. Sarà, G., Scilipoti, D., Mazzola, A., Modica, A., 2004. Effects of fish farming waste to sedimentary and particulate organic matter in a southern Mediterranean area (Gulf of Castellammare, Sicily): a multiple stable isotope study (δ13C and δ15N). Aquaculture 234, 199-213. SEO, J.Y., LEE, S.M., 2011. Optimum dietary protein and lipid levels for growth of juvenile sea cucumber Apostichopus japonicus. Aquaculture Nutrition 17, e56-e61. Shapiro, S.S., Wilk, M.B., 1965. An analysis of variance test for normality (complete samples). Biometrika 52, 591-611. Shirley, E., 1977. A non-parametric equivalent of Williams' test for contrasting increasing dose levels of a treatment. Biometrics, 386-389. Smith, V.H., Tilman, G.D., Nekola, J.C., 1999. Eutrophication: impacts of excess nutrient inputs on freshwater, marine, and terrestrial ecosystems. Environmental pollution 100, 179-196. Sun, Z., Gao, Q., Dong, S., Shin, P.K.S., Wang, F., 2013. Seasonal changes in food uptake by the sea cucumber Apostichopus japonicus in a farm pond: Evidence from C and N stable isotopes. Journal of Ocean University of China 12, 160-168. Tanner, J.E., Fernandes, M., 2010. Environmental effects of yellowtail kingfish aquaculture in South Australia. Aquaculture Environment Interactions 1, 155-165. Troell, M., Halling, C., Neori, A., Chopin, T., Buschmann, A., Kautsky, N., Yarish, C., 2003. Integrated mariculture: asking the right questions. Aquaculture 226, 69-90. Trott, L.A., McKinnon, A.D., Alongi, D.M., Davidson, A., Burford, M.A., 2004. Carbon and nitrogen processes in a mangrove creek receiving shrimp farm effluent. Estuarine Coastal and Shelf Science 59, 197-207. Uede, T., 2008. Validity of acid volatile sulfide as environmental index and experiment for fixing its standard value in aquaculture farms along the coast of Wakayama Prefecture, Japan. Nippon Suisan Gakkaishi 74,
AC C
581‐ 582‐ 583‐ 584‐ 585‐ 586‐ 587‐ 588‐ 589‐ 590‐ 591‐ 592‐ 593‐ 594‐ 595‐ 596‐ 597‐ 598‐ 599‐ 600‐ 601‐ 602‐ 603‐ 604‐ 605‐ 606‐ 607‐ 608‐ 609‐ 610‐ 611‐ 612‐ 613‐ 614‐ 615‐ 616‐ 617‐ 618‐ 619‐ 620‐ 621‐ 622‐ 623‐ 624‐ 625‐ 626‐ 627‐ 628‐ 629‐ 630‐ 631‐ 632‐
25‐ ‐
ACCEPTED MANUSCRIPT
M AN U
SC
RI PT
402-411. Wai, T.C., Leung, K.M.Y., Wu, R.S.S., Shin, P.K.S., Cheung, S.G., Li, X.Y., Lee, J.H.W., 2011. Stable isotopes as a useful tool for revealing the environmental fate and trophic effect of open-sea-cage fish farm wastes on marine benthic organisms with different feeding guilds. Marine Pollution Bulletin 63, 77-85. Walve, J., Larsson, U., 1999. Carbon, nitrogen and phosphorus stoichiometry of crustacean zooplankton in the Baltic Sea: implications for nutrient recycling. Journal of Plankton Research 21, 2309-2321. Wu, H., Peng, R., Yang, Y., He, L., Wang, W., Zheng, T., Lin, G., 2014. Mariculture pond influence on mangrove areas in south China: Significantly larger nitrogen and phosphorus loadings from sediment wash-out than from tidal water exchange. Aquaculture 426-427, 204-212. Yang, H.-s., Zhou, Y., Wang, J., Zhang, T., Wang, P., He, Y.-c., Zhang, F.-s., 1999. A modelling estimation of carrying capacities for Chlamys farreri, Laminaria japonica and Apostichopus japonicus in Sishilivan Bay, Yarntai, China. Journal of fishery sciences of China/Zhongguo shui chan ke xue yan jiu yuan bian ji chu ban 7, 27-31. Yearbook, C.F., 2016. China Bureau of Fisheries Beijing: China Agriculture Press China. Yokoyama, H., 2003a. Environmental quality criteria for fish farms in Japan. Aquaculture 226, 45-56. Yokoyama, H., 2003b. Environmental quality criteria for fish farms in Japan. Aquaculture 226, 45-56. Yokoyama, H., 2010. Monitoring, assessment and management of fish farm environments in Japan. Reviews in Aquaculture 2, 154-165. Yokoyama, H., 2013. Growth and food source of the sea cucumber Apostichopus japonicus cultured below fish cages — Potential for integrated multi-trophic aquaculture. Aquaculture 372-375, 28-38. Yokoyama, H., Ishihi, Y., 2007. Variation in food sources of the macrobenthos along a land–sea transect: a stable isotope study. Marine Ecology Progress Series 346, 127-141. Yokoyama, H., Nishimura, A., Inoue, M., 2007. Macrobenthos as biological indicators to assess the influence of aquaculture on Japanese coastal environments, Ecological and Genetic Implications of Aquaculture Activities. Springer, pp. 407-423. Yokoyama, H., Tamaki, A., Harada, K., Shimoda, K., Koyama, K., Ishihi, Y., 2005. Variability of diet-tissue isotopic fractionation in estuarine macrobenthos. Marine Ecology Progress Series 296, 115-128. Zhang, J., Liu, S., Zhang, Q., Liu, Y., 2008. Nutrient distribution and eutrophication assessment for the adjacent waters of the Yellow River Estuary. Marine Science Bulletin 27, 66-71. Zheng, Z., Dong, S., Tian, X., Wang, F., Gao, Q., Bai, P., 2009. Sediment water Fluxes of Nutrients and Dissolved Organic Carbon in Extensive Sea Cucumber Culture Ponds. CLEAN–Soil, Air, Water 37, 218-224.
TE D
633‐ 634‐ 635‐ 636‐ 637‐ 638‐ 639‐ 640‐ 641‐ 642‐ 643‐ 644‐ 645‐ 646‐ 647‐ 648‐ 649‐ 650‐ 651‐ 652‐ 653‐ 654‐ 655‐ 656‐ 657‐ 658‐ 659‐ 660‐ 661‐ 662‐ 663‐ 664‐
AC C
EP
665‐
26‐ ‐
ACCEPTED MANUSCRIPT
666‐
Figure captions
667‐ Fig. 1. Map of the study area showing sampling stations in sea cucumber farming area (sea cucumber
669‐
pond, effluent channel, stations with a gradient distance from the pond) and Natural reserve area (0m, 50
670‐
m, 100 m and 200 m from land).
671‐ 672‐
Fig. 2. Principal component analysis (PCA) based on the seventeen sediment and water parameters.
673‐
SC
Fig. 3. Dendrogram showing cluster relationship between samples. Fig. 4. Comparisons of water parameters. Fig. 5. Comparisons of sediment parameters.
M AN U
674‐ 675‐ 676‐ 677‐ 678‐ 679‐ 680‐ 681‐ 682‐ 683‐ 684‐ 685‐ 686‐
RI PT
668‐
Fig. 6. Stable isotope biplots of δ13C and δ15N for dominant consumers and their potential food sources in N-tidal flat (a) and C-tidal flat (b). Fig. 7. Proportionate contribution of potential food sources to the diet of Helice tridents and Macrophthalmus abbreviates by SIAR boxplots.
TE D
687‐ 688‐ 689‐ 690‐
EP
691‐ 692‐ 693‐
AC C
694‐ 695‐ 696‐ 697‐ 698‐ 699‐ 700‐ 701‐ 702‐
27‐ ‐
ACCEPTED MANUSCRIPT
Table 1. Chemical items analyzed in the present study and their abbreviations. Water
Particulate organic matter in the pond area
Sediment
POM(pond)
in the tidal flat area POM(tidal flat) in the coastal area
POM(marine) TOC(POM)
TOC(S)
Total nitrogen
TN(POM)
TN(S)
Carbon to nitrogen ratio
C/N ratio(POM) C/N ratio(S)
Stable carbon isotope ratio
δ13C(POM)
δ13C(S)
Stable nitrogen isotope ratio δ15N(POM)
δ15N(S)
Nitrate
NO3-(water)
Nitrite
NO2-(water)
Ammonium
NH4+(water)
Phosphate Acid volatile sulfides
PO43-(water)
M AN U
Total organic carbon
NO3-(S)
NH4+(S)
AP(S) AVS
AC C
EP
TE D
706‐
28‐ ‐
RI PT
Analyzing items
SC
703‐ 704‐ 705‐
ACCEPTED MANUSCRIPT
Table 2 Carbon and nitrogen isotopic compositions of consumers and their potential food sources in the tidal flat of Natural Reserve Area and sea cucumber pond area N-tidal flat C-tidal flat δ13C (‰) δ15N (‰) n δ13C (‰) δ15N (‰) n Consumers / producers a a b Helice tridens -15.5 ± 3.3 5.8 ± 0.4 12 -13.0 ± 0.6 4.3 ± 1.1b 21 Macrophthalmus abbreviates -12.3 ± 0.4a 3.8 ± 0.9a 10 -12.1 ± 0.6a 2.1 ± 0.7b 7 a a b POM (marine) -18.4 ± 0.3 1.5 ± 0.3 3 -21.1 ± 0.1 2.4 ± 0.3b 3 a a a b POM (tidal flat) -16.9 ± 1.2 2.7 ± 0.4 12 -16.7 ± 1.6 0.9 ± 0.6 18 POM (pond) -20.6 ± 0.9 -2.0 ± 0.2 3 BMA -8.1 ± 0.4a 2.2 ± 0.4a 3 -10.7 ± 1.5a 3.3 ± 0.4b 3 a a a a Salt marsh plant -30.4 ± 0.5 6.3 ± 1.5 9 -29.7 ± 1.3 6.1 ± 0.9 9 Plant debris -25.8 ± 0.9a 0.7 ± 0.4a 3 -17.7 ± 1.6b -0.1 ± 0.4a 3 The mean ± SD values are given for food items with the number of samples (n) ≥ 3. BMA: benthic microalgae. Different letters (a, b) on stable isotope values indicate significant difference (P < 0.05) by t-test.
M AN U
TE D EP AC C
711‐ 712‐ 713‐ 714‐
SC
RI PT
707‐ 708‐ 709‐ 710‐
29‐ ‐
EP
Fig. 1. Map of the study area showing sampling stations (sample numbers from 1 to 36) in four areas. Sea cucumber pond area: sea cucumber pond (pond), effluent channel (EC); C-tidal flat area: 0 m to 800 m; N-tidal flat area: 0 m to 200 m; farming coastal area (CM); non-farming coastal area (NM).
AC C
715‐ 716‐ 717‐ 718‐ 719‐ 720‐ 721‐ 722‐
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
30‐ ‐
Fig. 2. Principal component analysis (PCA) based on the seventeen sediment and water parameters. (a) 3D scores plot (PC1, PC2 and PC 3) of the environmental variables to the samples. (b) PC1 and PC2. (c) PC1 and PC3.
AC C
723‐ 724‐ 725‐ 726‐ 727‐ 728‐ 729‐ 730‐ 731‐ 732‐ 733‐ 734‐ 735‐ 736‐ 737‐
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
31‐ ‐
ACCEPTED MANUSCRIPT
EP
TE D
Fig. 3. Dengrogram showing cluster relationship between samples. Different sampling locations were shown in different colors of solid circle (green = N-tidal flat, blue = C-tidal flat, yellow = effluent channel, red = sea cucumber pond). See Fig. 1 for the sample numbers.
AC C
739‐ 740‐ 741‐ 742‐ 743‐ 744‐ 745‐ 746‐ 747‐ 748‐ 749‐ 750‐ 751‐ 752‐ 753‐ 754‐ 755‐ 756‐ 757‐ 758‐ 759‐ 760‐ 761‐ 762‐ 763‐ 764‐ 765‐ 766‐ 767‐ 768‐ 769‐ 770‐
M AN U
SC
RI PT
738‐
32‐ ‐
AC C
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
Fig. 4. Comparisons of water parameters. Different letters above bars indicate significant differences (p<0.05) by Bonferroni and Gao test. Boxes represent 25%-75% percentiles; range bars represent the 5% and 95% percentiles, solid dots in the boxes represent mean values, hollow dots represent values outside of 95% confidences interval.
771‐
33‐ ‐
772‐ 773‐ 774‐ 775‐ 776‐ 777‐ 778‐ 779‐ 780‐ 781‐
AC C
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
Fig. 5. Comparisons of sediment parameters. Different letters above bars indicate significant differences (p< 0.05) by Bonferroni test and Gao test. Boxes represent 25%-75% percentiles; range bars represent the 5% and 95% percentiles, solid dots in the boxes represent mean values, hollow dots represent values outside of 95% confidences interval.
34‐ ‐
ACCEPTED MANUSCRIPT
784‐ 785‐ 786‐ 787‐ 788‐ 789‐ 790‐ 791‐ 792‐ 793‐ 794‐
AC C
EP
TE D
M AN U
SC
RI PT
782‐ 783‐
Fig. 6. Stable isotope biplots of δ13C and δ15N for dominant consumers and their potential food sources in N-tidal flat (a) and C-tidal flat (b). Grey solid dots: expected isotope composition of the diet of each crab. Error bars represent SD.
35‐ ‐
ACCEPTED MANUSCRIPT
797‐ 798‐ 799‐ 800‐ 801‐ 802‐ 803‐ 804‐
AC C
EP
TE D
M AN U
SC
RI PT
795‐ 796‐
Fig. 7. Proportionate contribution of potential food sources to the diet of Helice tridents and Macrophthalmus abbreviates by SIAR boxplots. Plant: salt marsh plant. Debris: plant debris. Dark gray, gray and light gray boxes show 50, 75, 95 credibility interval respectively.
805‐
36‐ ‐
ACCEPTED MANUSCRIPT
Highlights
Sediment biogeochemical characteristics are effective indices for estimation of sea cucumber pond aquaculture impact Dispersion distance of aquaculture waste is restricted within a distance
Stable isotope approaches are efficient in tracing aquaculture waste in environment
AC C
EP
TE D
M AN U
SC
RI PT