2.1 Impacts on CO2 fluxes
Saltmarshes are among the most productive ecosystems in the world with the capacity to store large amounts of C per unit area, referred as “Blue Carbon” 2. Our high-resolution eddy covariance (EC) measurements of CO2 fluxes revealed that the S. alterniflora marsh was a net sink of atmospheric CO2. Daily CO2 fluxes in the marsh exhibited distinct seasonal patterns, with peak CO2 fluxes, e.g., gross primary productivity (GPP), ecosystem respiration (Re), and net ecosystem CO2 exchange (NEE), occurring in summer when vegetation has grown vigorously (Fig. 1a). From late-September, CO2 fluxes gradually declined since both GPP and ER were inhibited by cool temperature. Because of the long green period of S. alterniflora, the saltmarsh switched from net CO2 source to net CO2 sink in early April and did not switch back to net CO2 source until December. Over an annual cycle, 65.2 Mg CO2 ha‒1 (1,778 g C m‒2) was accumulated as GPP and 47.4 Mg CO2 ha‒1 (1,293 g C m‒2) was released as Re in this saltmarsh, resulting in a net uptake of 17.8 Mg CO2 ha‒1 (485 g C m‒2) annually (Fig. 2a).
In contrast to S. alterniflora marsh, NEE in the mariculture ponds fluctuated between ‒34.1 and 215.3 mg CO2 m‒2 h‒1 and no apparent seasonal variations were observed, except for the pulse of CO2 emission in the initial drainage period (Fig. 1b). On an annual basis, mariculture ponds were net atmospheric CO2 sources ranging from 0.91 to 1.76 Mg CO2 ha‒1, and the highest CO2 emission was observed in the 1-year old mairculture pond (MP1) (Fig. 2a). The NEE in our ponds compared reasonably with observations from other mariculture systems which also receives feed inputs. Zhang et al.25 reported an annual CO2 emission of 4.69 Mg CO2 ha‒1 in a nearby crab-shrimp-clam polyculture system. Recent seasonal measurements from subtropical estuarine shrimp ponds also demonstrated culturing-period CO2 emissions ranging from 0.29 to 0.55 Mg CO2 ha‒1 (ref. 26). Aquaculture waters have been documented autotrophic and act as a sink of atmospheric CO2, because primary production of phytoplankton is generally higher than the water-column respiration26,27. However, phytoplankton is easily decomposed compared to plant biomass, resulting in minimal retention of primary production. Yang et al.26 reported a net water-column CO2 fixation of 2.02 Mg CO2 ha‒1 in a subtropical shrimp pond, which was ~ 1 order of magnitude lower than the CO2 uptake rate in S. alterniflora marsh. Therefore, the shift of CO2 sink to source following saltmarsh conversion can be primarily due to the loss of CO2 sequestrating capacity of primary producer. Specifically, prevalence of CO2 efflux in the ponds suggested that water-column CO2 uptake was largely counterbalanced and surpassed by other components of respiration, e.g., soil respiration26. It has been estimated that only < 20% of the feed C can be converted to biomass in aquaculture systems, the residuals are buried and then transformed to CO2 and CH4 by microbes17,28. However, we found that NEE was significantly higher in MP1 than in the 4-year (MP2) and 5-year (MP3) ponds (Fig. 2a). As C inputs from feeds and fertilizers were similar among ponds, the difference in NEE indicated that the decomposition of residual fish feed and feces can only be partly responsible for the CO2 emission. Decomposition of SOC, on the other hand, might account more for the CO2 efflux. This can be supported by the significant lower SOC contents in the older ponds relative to MP1 (Table 1). Losses of SOC have been widely reported in wetlands subjecting to reclamation and conversion to other land uses29,30. Arifanti et al.31 found that the SOC mineralization rate over 16 years following mangrove conversion to shrimp ponds was up to 24.5 Mg C ha‒1 yr‒1, accounting for 75.0% of total C losses. When coastal wetlands are converted into aquaculture ponds, the top layers of soil are excavated and exposed to oxygen which enhances the susceptibility of SOC to mineralization. Further, operation activities such as aeration, drainage and harvest could also destabilize the SOC, resulting in large amounts of CO2 emission during the initial culturing years32. Thereafter, subsequent losses were much slower probably due to the increased proportion of recalcitrant organic C fractions33,34, as exhibited in the relatively lower NEE in MP2 and MP3.
Table 1
Soil properties in the Spartina alterniflora saltmarsh and mariculture ponds.
Systems
|
Culture time
(year)
|
Area
(ha)
|
SOC
(g C kg–1)
|
TN
(g N kg–1)
|
DOC
(mg C kg–1)
|
NH4+
(mg N kg–1)
|
NO3–
(mg N kg–1)
|
WD
(cm)
|
DO
(mg L–1)
|
Salinity
(ppt)
|
Saltmarsh
|
–
|
–
|
15.8 ± 0.2a
|
1.04 ± 0.04b
|
135 ± 9a
|
2.31 ± 0.18b
|
0.14 ± 0.03b
|
0.12 ± 0.47b
|
–
|
28.2 ± 0.1b
|
MP1
|
1
|
3.87
|
13.5 ± 0.1b
|
0.82 ± 0.03a
|
47.0 ± 3.8b
|
24.4 ± 1.3a
|
1.41 ± 0.21a
|
88.2 ± 4.0a
|
6.83 ± 0.25a
|
30.8 ± 0.5a
|
MP2
|
4
|
4.83
|
11.9 ± 0.2c
|
0.78 ± 0.01a
|
36.3 ± 3.0c
|
26.8 ± 1.9a
|
1.35 ± 0.23a
|
88.9 ± 5.6a
|
8.12 ± 0.33b
|
29.0 ± 0.4b
|
MP3
|
5
|
4.47
|
11.1 ± 0.2d
|
0.77 ± 0.01a
|
35.0 ± 3.2c
|
28.2 ± 2.0a
|
1.23 ± 0.21a
|
82.4 ± 5.3a
|
8.89 ± 0.31a
|
28.6 ± 0.4b
|
F-value
|
–
|
–
|
104.14
|
30.17
|
55.09
|
115.75
|
18.40
|
598.30
|
12.07
|
19.34
|
P-value
|
–
|
–
|
< 0.001
|
< 0.001
|
< 0.001
|
< 0.001
|
< 0.001
|
< 0.001
|
< 0.001
|
< 0.001
|
MP, mariculture pond; SOC, soil organic carbon; TN, total nitrogen; DOC, dissolved organic carbon; NH4+, ammonium; NO3–, nitrate; WD, water depth; DO, dissolved oxygen. Different letters denote significant differences (P < 0.05) between groups based on ANOVA followed by Tukey’s HSD test. Values are means ± 1 standard errors during the entire year but not for SOC and TN. |
It is noteworthy that the activities of animals, bivalve molluscs in particular, can also accelerate SOC decomposition in mariculture systems. Declines of quantity and quality of SOC have been documented in the no-feeding shellfish farms in a subtropical estuary (Ao River estuary, China)35. Bioturbation and bioirrigation associated with the burrowing, feeding, and other activities of molluscs can mix, flush and suspend the sediments, resulting in oxygenation of sediments and elevated accessibility of organic matters to microorganisms, which favors the sediment respiration and decomposition of SOC (ref. 36).
2.2 Impacts on CH4 fluxes
Significant differences were also observed in the magnitude and seasonal pattern of CH4 fluxes between S. alterniflora marsh and mariculture ponds (Fig. 1c). In the saltmarsh, daily CH4 fluxes gradually increased from the start of the growth season in April, peaked on August, and then decreased to near zero during the wintertime. In the ponds, CH4 fluxes were low at the beginning of the culturing period, then steadily increased and reached the first peak in July-August. From mid-August, CH4 fluxes declined and were maintained low until draining. Drainage triggered the second peaks of CH4 emission and the highest CH4 flux was measured in MP1 (Fig. 1d). Annual CH4 emission in the S. alterniflora marsh was 144 kg CH4 ha–1, 83.8% of which occurred in growth season. Converting saltmarsh to mariculture ponds dramatically decreased CH4 emission by 80.3‒96.0%, with annual CH4 emission in the ponds ranging from 5.77 to 28.3 kg CH4 ha–1. Similar to NEE, CH4 emission exhibited a decreasing trend with culturing time, which was 3.91 times greater in MP1 than in MP3 (Fig. 2b). CH4 emission during culturing period in the ponds was 3.93–19.2 kg CH4 ha–1, contributing 62.2–68.1% annual budget. While CH4 emission in drainage period mainly occurred in the initial 10 days, contributing more than half of the emission (Figs. 1c and 2b).
Poffenbarger et al.37 proposed a salinity threshold of > 18 ppt for coastal wetlands that CH4 emission could be neglected, which has recently been adopted by IPCC (ref. 20). However, our results suggested that vegetated coastal saltmarsh could still emit appreciable CH4 albeit the high salinity (28.4 ppt). In a previous study, we found that 67.3% of CH4 produced in S. alterniflora marsh originated from ‘noncompetitive’ trimethylamine38, as methanogens are often outcompeted for acetate or H2 by sulfate reduction in marine environments. These methylamines can be produced during the decomposition of the cytoplasmic osmolytes (e.g., glycine betaine, choline, and proline) of marine plants39, supporting substantial CH4 production in S. alterniflora marsh. We had expected that converting S. alterniflora marsh to mariculture ponds would increase CH4 emission, as anthropogenic labile organic C inputs and continuous inundation may favor CH4 production in the ponds. In contrast to our expectations, this conversion significantly reduced CH4 emission (Fig. 2b), which could be related to the decreased total C inputs and substrate availability. In S. alterniflora marsh, annual gross C assimilation was 17.8 Mg C ha‒1, which was 18 times higher than the C inputs from feeds and fertilizers in mariculture ponds (0.73 Mg C ha‒1, Supplementary Table 1). Here, C inputs from primary production of phytoplankton and periphyton in the ponds were not quantified. However, previous study reported an annual net autotrophic C fixation of 182‒218 g C m‒1 in an analogous crab-shrimp-clam polyculture system adjacent to our sites40. Therefore, we consider that when primary productivity in the ponds was counted, C inputs for mariculture ponds were still much lower than that in S. alterniflora saltmarsh, which eventually decreased the supply of methanogenic substrates. We have speculated that, compared with plant residues, organic compounds in the feeds (e.g., starch and protein) can be more easily decomposed and become as substrates, especially “noncompetitive” methylamines17,23. Indeed, CH4 emission factor of applied/inputted C (EFC) in test ponds (1.43%) was 2.35 times that in S. alterniflora marsh (0.61%), indicating that the inputted C could be more efficiently converted to CH4 in mariculture ponds. However, the considerably lower CH4 emissions in the mariculture ponds suggested that the benefit of higher C decomposability was completely offset by the loss of primary production.
More important, conversion of S. alterniflora saltmarsh to mariculture ponds significantly reduced SOC by 14.6‒29.7%, and especially DOC by 65.2‒74.1% (Table 1). Although the importance of recent photosynthates as source of CH4 has been widely reported, SOC still accounted for a large fraction of CH4 production in wetlands41,42. For example, 40‒96% of emitted CH4 from paddy fields originated from SOC (ref. 42,43). Consistently, Kankaala and Bergström44 found that the SOC-derived CH4 contributed more than 90% to the total CH4 emission from a shallow boreal lake. Particularly, CH4 emission, SOC and DOC concentration in the ponds decreased with culturing age, and MP1 had the higher CH4 emission and organic C contents relative to MP2 and MP3. Hence, SOC decomposition partly, at least, supported CH4 production, since the application rates of feed and fertilizer were identical. Alongside with decreased substrate availability, changes in hydrological condition may also contribute to the reduction of CH4 emission. Compared with the intermittently flooded saltmarsh, the ponds were permanently inundated during the culturing period (Supplementary Fig. 2). Generally, permanent inundation can create ideal anaerobic condition favoring CH4 production45. However, the standing water can also act as a diffusion barrier which may suppress CH4 emission by increasing CH4 oxidation during transport46. Here, the maximal CH4 fluxes in the ponds were observed in the early drainage period (Fig. 2d). Given the low temperature of this stage (early winter) and aerobic/semi-anaerobic condition, we consider that the emitted CH4 was not currently produced but was a release of trapped CH4 produced in culturing period. Inevitably, prolonged transport process increased the chance of CH4 being oxidized. Nykänen et al.47 found that up to 97.5% of CH4 produced in a shallow small lake was oxidized before reaching the atmosphere. In contrast, CH4 transport could proceed quickly in saltmarsh via aerenchyma tissue of S. alterniflora, which may also help bypass oxidation in soil zones48.
2.3 Impacts on N2O fluxes
In S. alterniflora saltmarsh, N2O fluxes were at a low level and fluctuated within a narrow range between weak uptakes (‒24.1 µg N2O m‒2 h‒1) and emissions (27.7 µg N2O m‒2 h‒1), while uptake frequently occurred in summer and autumn (Fig. 1e). In mariculture ponds, basal manure application and the first supplemental fertilizing did not stimulate N2O fluxes. After the second supplemental fertilizing, a peak of N2O flux was observed, albeit the relatively lower rates. Similar to CH4 fluxes, a burst of N2O was observed immediately upon drainage (Fig. 1f). Annually, S. alterniflora marsh was a weak sink of atmospheric N2O (‒0.04 kg N2O ha–1) and 68.3% of the N2O uptakes occurred in growth season (Fig. 2c). However, converting saltmarsh to mariculture ponds produced a shift to N2O source. Annual N2O emission in the ponds ranged from 1.06 to 1.11 kg N2O ha–1. In spite of the shorter time and lower temperature, drainage period N2O emission was higher than that in culturing period. Particularly, the initial 10 days of drainage period contributed 22.1‒31.4% annual N2O budget (Fig. 2c).
Our annual N2O budget in the S. alterniflora marsh was comparable to the reported atmospheric N2O consumptions of 0.01‒0.14 kg N2O ha‒1 yr‒1 in coastal saltmarshes which receive limited anthropogenic N loading49,50. These data supported the view of productive and pristine saltmarshes as negligible sources or even small sinks of atmospheric N2O (ref. 49). However, converting S. alterniflora marsh to mariculture ponds fundamentally shifted the N2O sink to a net source of N2O. At least two explanations can be proposed. First, soil mineral N for nitrifiers and denitrifiers was considerably increased following the conversion (Table 1). Generally, N is the key limiting nutrient in pristine coastal wetlands. Here annual mean NO3‒ in the S. alterniflora marsh was only 0.14 mg N kg‒1, a value that was lower than the threshold (5 mg N kg‒1) for denitrification51 and was close to the threshold (~ 1 mg N kg‒1) for denitrifiers using N2O as the unique electron acceptor52. Hence, the lack of NO3‒ might be responsible for negative N2O fluxes in S. alterniflora marsh. Conversely, mineral N concentrations in the ponds were a magnitude higher (Table 1). The mariculture ponds received higher N inputs from feeds and fertilizers while only a small part of the added N (~ 25%) was converted to fish biomass, leading to large amount of residual active N accumulated and facilitating N2O production in the ponds16,17. Second, Chen et al.53 found that N2O production via denitrification relies on NO3‒/DOC ratio in soils. The relatively lower DOC concentration and higher NO3‒/DOC ratio in ponds could increase the proportion of N2O in denitrification products since high DOC content increases electron pressure, facilitating N2O reduction and resulted in a lower ratio of N2O/(N2O + N2) or even net N2O consumption, especially in anaerobic environments54. Senbayram et al.55 pointed out that when soil NO3‒ concentrations fell below 20 mg kg‒1, organic substrate addition would sharply reduce N2O emission in anaerobic soils. In aquaculture systems, the lack of organic electron donors for N2O reduction was also recognized as a major reason for N2O emissions56, whereas carbohydrate (starch) addition was noted to dramatically reduce N2O emission by 83.4% in an intensive tilapia culturing system57. Therefore, it is likely that the elevated NO3‒ and lowered DOC availability collectively increased N2O emission via denitrification following the conversion of saltmarsh to mariculture ponds.
As well as with the change in magnitude, N2O fluxes in mariculture ponds showed distinct temporal patterns from that in S. alterniflora marsh. Substantial N2O fluxes were triggered by pond draining. This finding agreed well with previous studies that observed drainage episodes as the critical period for N2O emission from aquaculture ponds21,22,58. As discussed above, abundant soil NO3‒ usually generates a higher ratio of N2O/(N2 + N2O) because NO3‒ is preferred over N2O as an electron acceptor59. The ratio also increases if some oxygen (O2) is present, because the inhibition of N2O reductase by O2 is stronger than the other reductases involved in denitrification60. On the other hand, the draining pond soil was still saturated initially and might cause optimum moisture condition (~ 80% water-filled pore space in soil) for a large proportion of N2O produced during denitrification58,59. Furthermore, with no standing water, N2O could diffuse more easily and escape from consumption compared to the waterlogging condition. Hence, we consider these factors have collectively contributed to the N2O emission pulse in the early stage of the drainage period.
2.4 Full GHG debt of mariculture ponds
Although aquaculture is often inferred as important sources of GHG to the atmosphere61,62, only a few studies could comprehensively account for the GHG emissions associating with the conversion of natural coastal wetlands. Here, we incorporated the direct fluxes of CO2, CH4, and N2O to evaluate their net global warming potential (nGWP) of each system. Annually, mariculture ponds exerted an overall positive net radiative forcing impact of 2.16 Mg CO2eq ha–1 over a 100-year time horizon, with CO2, CH4, and N2O accounting for 68.9%, 18.0%, and 13.2% of the nGWP, respectively (Fig. 2d). Among the ponds, MP1 held up the highest nGWP (3.22 Mg CO2eq ha–1 yr–1), primarily due to its higher NEE and CH4 emissions relative to the older ponds. In contrast to mariculture ponds, S. alterniflora marsh showed a negative nGWP of ‒13.8 Mg CO2eq ha–1 yr–1. This substantial negative value was mainly attributed to the high CO2 uptake by plants, despite 22.6% of its climatic cooling impact was offset by the simultaneous CH4 emission. Specifically, the negative nGWP of S. alterniflora marsh indicated that its conversion to mariculture ponds caused the loss of foregone GHG mitigation capacity by the saltmarsh63, yielding a full GHG debt of 15.9 Mg CO2eq ha–1 yr–1. This value was 7.36-fold higher than their direct GHG emissions from the mariculture ponds.
2.5 Comparison with other animal protein production sectors
We compared the direct GHG emissions from our mariculture ponds with those from other aquaculture systems in China, e.g., freshwater crab ponds (< 0.5 ppt) and brackish shrimp ponds (4.2 ppt) (ref. 17,26,64). We found that, in spite of the differences in C and N inputs and operation scenarios, nGWP and yield-scaled GWP (GHGI) in aquaculture systems decreased significantly along the gradient of salinity (Fig. 3a,b). Our estimates of nGWP and GHGI in the polyhaline mariculture ponds (28.4 ppt) were just 8.51‒13.7% and 1.38‒14.9% of those in freshwater and brackish ponds. This difference was mainly attributed to the minor CH4 emission in our ponds, given N2O and CO2 emissions were much comparable among three systems. Salinity can strongly impact CH4 production by regulating the availability of sulfate. In contrast to polyhaline water ponds where CH4 production was inhibited by presence of sulfate, lower salinity and sulfate availability in freshwater and brackish aquaculture ponds allowed intense CH4 production from H2/CO2 and acetate, whereas these “competitive” substrates can be easily decomposed from residues of feeds and submerged macrophytes17,23. This result substantiates the role of CH4 emission in GHG budgets of aquaculture systems and suggests that mariculture at high salinity could have substantial lower direct GHG-cost relative to freshwater and brackish aquaculture due to their virtual low CH4 emissions.
We also compared the full GHG debt in aquaculture systems with the estimates from other animal protein production sectors on a kg CO2eq kg‒1 protein basis (GHGI-protein) (Fig. 3c). Similar to nGWP and GHGI, the direct GHGI-protein in our mariculture ponds (6.48 kg CO2eq kg‒1 protein) was significantly lower than that in freshwater crab (185 kg CO2eq kg‒1 protein) and brackish shrimp (14.7 kg CO2eq kg‒1 protein) ponds, but was higher than that in a no-feeding oyster farm (0.13 kg CO2eq kg‒1 protein) (ref. 65). Notably, we found the GHGI-protein in aquaculture systems were much lower than those from terrestrial livestock production (39.5‒466 kg CO2eq kg‒1 protein), except for freshwater crab ponds. However, when the loss of foregone GHG mitigation capacity of pristine wetlands was counted, GHGI-protein in our mariculture ponds would increase to 46.7 CO2eq kg‒1 protein, which is higher than poultry but still lower than beef, small ruminants, or pork production. These results imply that mariculture can provide a low-GHG animal protein source relative to freshwater aquaculture and livestock production. Developing approaches to attenuate net emissions associated with the foregone GHG mitigation loss will be crucial for minimizing the climate impact of mariculture. Future aquaculture programs could preferentially engineer aquaculture systems by conversion of lowly productive sites (e.g., saline bare land, abandoned ponds, etc.) and carrying out it in high-salinity waters30,31, as reduced foregone C sequestration losses and CH4 emissions could be simultaneously achieved under such conditions.
Overall, our results show that the conversion of saltmarsh to mariculture ponds in coastal region of China creates a GHG debt of 15.9 Mg CO2eq ha–1 yr–1. Within them, the loss of foregone GHG mitigation capacity of pristine wetlands makes the largest contributions (86.4%) to overall GHG debt, whereas the other 15.6% arises from the direct emission of CO2, CH4, and N2O in the mariculture ponds. This finding contributes to growing recognition of the climatic roles of mariculture and highlights the importance of taking the foregone C sequestration losses associated with coast wetlands conversion into the full GHG accounting. Nonetheless, the GHG-cost of protein in mariculture was much lower than inland freshwater aquaculture and the major terrestrial livestock productions (such as beef, small ruminants, and pork), indicating mariculture can provide a low-GHG animal protein source. For initiatives aiming at reduction of GHG emissions, future aquaculture should primarily avoid conversion of the highly productive coastal wetlands, and further minimize CH4 emissions by using high-salinity waters.