Reviewer #1:

'Comment on egusphere-2025-2864', Anonymous Referee #1, 6 Aug 2025

Fecal pellets act as important particulate organic matter components, and play vital role in marine carbon cycling. This manuscript presented the downward fluxes of zooplankton fecal pellets in the southern South China Sea, based on a long time sediment trap mooring system. The results are very helpful to understand the importance of fecal pellets in organic carbon export and cycling. However, there are some defects through the manuscript. Below are my major advice:

Reply: We sincerely appreciate your effort reviewing our manuscript and your positive feedback on our research topic. We have read through all your comments carefully and have made related modifications to improve the quality of our manuscript. We highly appreciate your time and consideration.

Comment #1:

The potential impacts of monsoon, typhoon and eddies to fecal pellets export fluxes are discussed. However, through the discussion section, their importance is only qualitatively discussed, and no direct quantitative analysis is done. E.g., in section 4.1, the authors focused on comparing the MLD during monsoon and non-monsoon periods, and then the relationships nitrate and Chl a with MLD. I advise the authors quantify the effects of the above events in enhancing the FPN and FPC using statistical analysis/test.

Reply: We sincerely thank the reviewer for this valuable and insightful comment. We acknowledge that the initial version of our manuscript primarily provides a qualitative discussion of the effects of the winter monsoon, typhoons and mesoscale eddies on fecal pellet fluxes, without a direct quantitative assessment.

In response, we have now incorporated several additional quantitative approaches in Section 4.4 (Lines 497-519) [see Supplement for revised manuscript annotated] to better evaluate the influence of these dynamic events on fecal pellet numerical flux (FPN) and fecal pellet carbon flux (FPC).

1. Correlation among hydrological parameters

As mentioned in the comment, the relationship between MLD, nitrate, and Chl a concentration is intriguing yet not quantified in our manuscript. Thus, we conducted Pearson correlation analyses among key hydrological parameters, including mixed layer depth (MLD), sea level anomaly (SLA), sea surface temperature (SST), zooplankton biomass (ZMC), and calculated all correlation matrix. The results are presented in Figure 9 in Line 513.

2. Relationships between fecal pellet fluxes and hydrological parameters

In addition, we applied Mantel tests to further access the potential relationships between fecal pellet fluxes (FPN, FPC) and hydrological parameters. We found significant correlations between fecal pellets and SLA, ZMC, and other components of the BCP, such as opal, carbonate, and terrigenous fluxes. These results reinforce our previous hypothesis that upper-ocean dynamics strongly modulate particle export, and that various BCP components tend to increase simultaneously under such influences. The results, also shown in Figure 9, are discussed in Lines 498-504.

3. General linear model (GLM):

To quantify the relative contributions of upper-ocean dynamic events, we employed a general linear model (GLM). Our results indicate that The East Asian monsoon (EAM) exerts a demonstrated positive influence, explaining approximately 42 % of FPN and over 32 % of FPC variability. The model outputs are summarized in Table 2 (Line 516), with interpretations provided in Lines 505-512. However, we note that the limited sample size (n = 13) may constrain model robustness and introduce potential overfitting, which we discuss in Lines 510-512, emphasizing the need for future studies with larger datasets.

Comment #2:

The language and logic of the manuscript must be improved before submitting and published to such high quality journal.

Reply: We sincerely thank the reviewer for this valuable suggestion. To meet the standards of such a high-quality journal, we have made comprehensive improvements to both the language and logic structure of the manuscript.

1. Language:

We carefully revised the entire manuscript to enhance readability and academic tone. The language was refined to be more professional and concise, with redundant expressions removed. We also adjusted sentence structures to improve clarity and coherence. Since revisions were made throughout the text, we provide one sample for illustration:

Previous version (Section 2.1):

‘The South China Sea (SCS), the largest semi-enclosed marginal sea in the western Pacific, spans 3.5 × 10⁶ km² with an average depth of 1140 m (Wang and Li, 2009). Our study focuses on the southern SCS near the Sunda Shelf (Fig. 1a), where the East Asian monsoon system (EAM) dominates surface wind fields and ocean circulation (Shaw and Chao, 1994).’ 

Revised version (Lines 81-85):

‘Mooring TJ-S (6.72ºN, 110.76ºE, 1630 m water depth) is located in the southern SCS near the Sunda Shelf (Fig. 1a), where the East Asian monsoon system (EAM) exerts dominant control on local wind fields and ocean circulation (Shaw and Chao, 1994).’

This revision removed redundant background information and strengthens focus on the study site.

2. Logic:

We also reorganized sections of the manuscript to improve logical flow and strengthen discussions. We now improved the logic of BCP introduction (Lines 19-29), refined the description of methodological approaches (Lines 33-41), and enhanced the logical flow in describing hydrological data sources (Lines 144-175). We have now strengthened the logic in the introductory part of the discussion section (Lines 296-323), and we also reorganized paragraphs and removed redundant information in the further discussion to maintain a clear focus on major arguments (Lines 325-338, 396-402).

Comment #3:

The manuscript only simply presented the data sources, such as wind field, primary production…. But, how the trap system pretreated before using? In other words, how the bottles cleaned and did any chemicals (e.g. HgCl₂ or NaCl solution) added to each bottle? Since the storing durations of the particles in the bottles are different and the degradation percent in each bottle should be considered. In addition, POC is used in Figure 4, but how was it analyzed is not mentioned in the method section. 

Reply: We thank the reviewer for this valuable comment. In response, we have expanded Section 2.2 (Lines 111-142) to include detailed descriptions of the pretreatment and handling of the sediment trap system as well as procedures for POC measurement. Specifically, we have added information on the use of NaCl-buffered HgCl₂ solution prior to deployment as a preservative to minimize biological degradation during sample storage (Lines 114-117). We also added detailed steps involved in sample processing after retrieval, including sample storage conditions, rinsing and splitting of subsamples, and drying procedures (Lines 122-127). The analytical method for total mass flux (TMF) and POC determination are also described in Lines 128-130.

Comment #4:

There are some mistakes/typos in certain figures: There is no panel d in Figure 4. The y-axis names are repeated in Figure S7, the right one should be wrong. Base on the contents, Figure 3a and Figure S7d are repeated. Figure 3a and 3b, 3d and 3e are also repeated. Figure 5 and Figure 7 are partially repeated and the logic/order in Figure 5 are bad (they should be presented in the order of physical parameters, then biochemical parmeters, and raw data and then reanalized data). Additionally, the panels should be cited for the first time in the main text sequentially.

Reply: We thank the reviewer for this careful examination and for pointing out the errors in our figures. In response, we have thoroughly checked and revised all the figures in both the main text and the supplementary materials: Figure 3 has been reorganized to avoid data repetition: Fig. 3a now presents the total FPN flux, while Fig. 3b now shows the individual contributions of different fecal pellet types (Line 237). All panels in Figure 4 have been carefully reviewed and corrected (Line 256). The orders of hydrological parameters in Figure 5 have been completely reorganized to improve logic flow. As suggested, physical parameters and raw data are now presented first, followed by reanalyzed data and biogeochemical parameters (Line 289). In Figure S7, panel d has been removed to eliminate redundancy (Line 110, supplementary materials). The right y-axis is now correctly labeled as representing the mean biovolume of the single pellet (the line plot), while the left y-axis represents the biovolume distribution (bar plot). Additionally, we have ensured that all panels are now cited sequentially in the main text according to their first appearance.

Comment #5:

Please add note for the bold numbers in Table 1. I think they the average values. Please also present the standard deviations.

Reply: We thank the reviewer for this helpful suggestion. We have now clarified in the caption of Table 1 (Lines 234-235) that the bold numbers represent average values. In addition, we have now included the corresponding standard deviations to provide a more comprehensive statistical description.

Comment #6:

Lines 25-26: “sequestering over 1300 Pg C” Is this the integrating value? If yes, what’s the time scale?

Reply: We thank the reviewer for pointing out this issue. The 1300 Pg C value is assimilated by an ensemble numerical model (Nowicki et al., 2022), and it is an integrating value with an average sequestration time of  127₁₂₂¹³³ years. Detailed calculation methodology is described in section 2.6.2 of their research by using the equation below:

dc/dt = Ac + J_DIC - Q_surf^C

where A is the transport matrix and c is DIC from respired organic matter. J_DIC is the source of DIC from organic matter respiration, and Q_surf is a matrix operator that removes carbon from the system once it reaches the surface layer. J_DIC is obtained by solving the model equations for each carbon sequestration pathway in the same way as the export pathways are partitioned, using a diagnostic linear model to avoid issues associated with nonlinearities in the original model equations (Nowicki et al., 2022).  

We have realized that using this number may be not proper to merely demonstrate the importance of BCP in this paragraph. We have now reorganized the whole paragraph into: “The ocean plays a pivotal role in regulating global carbon sink, absorbing 2.9 ± 0.4 Gt C annually through coupled physical and biological mechanisms, mitigating increasing anthropogenic carbon dioxide (CO₂) emissions (Friedlingstein et al., 2025). Central to this uptake lies the biological carbon pump (BCP), which converts massive dissolved CO₂ in the surface ocean into particulate organic carbon (POC) via phytoplankton photosynthesis (Falkowski, 2012; Boyd and Trull, 2007; Nowicki et al., 2022).” (Line 20-24)

Comment #7:

Lines 37-38: The authors mentioned “These studies”, but did not cite any references, please add the relevant references. In addition, why the contribution of fecal pellets to POC can > 100%?

Reply: We thank the reviewer for the careful examination and for pointing out the missing citations. Accordingly, we have now added relevant references after this sentence (Lines 33-41), including Shatova et al., 2012; Siegel et al., 2014; Turner et al., 2015; Stamieszkin et al., 2015; Estapa et al., 2017; Li et al., 2022; Countryman et al., 2022; Terrats et al., 2023; Darnis et al., 2024.

Regarding the FPC/POC ratio occasionally exceeding 100 %, this phenomenon arises primarily from methodological rather than ecological causes. Theoretically, FPC is a component of total POC, thus the actual FPC/POC ratio should not surpass 100 %. However, FPC and POC are quantified using different approaches. POC is directly measured from bulk samples, while FPC flux is estimated from pellet volumes converted to carbon contents using carbon-volume conversion factors. Uncertainties in pellet size measurement, image recognition, and the choice of conversion factors can lead to over estimations of FPC. In some cases, especially when fecal pellets dominate the local BCP process, these uncertainties may produce calculated FPC/POC ratios close to or slightly above 100 %, though the true ratio cannot exceed this limit.

Comment #8:

Line 110: Based on the location of the mooring station and depth, terrigenous particles may be not the only component of the fine particles. I mean biological originated organic/inorganic particles, including detritus of phytoplankton, partially degraded organic particles, etc, may contribute more to this fraction. 

Reply: We thank the reviewer for raising this insightful comment. We fully agree that, given the location and depth of the TJ-S mooring station, the fine-particle fraction in our sediment trap samples likely includes not only terrigenous materials, but also biologically derived components such as phytoplankton detritus and partially-degraded organic and inorganic materials. We apologize for not clearly describing all the relevant pretreatment procedures in the original version.

In the revised manuscript, we have expanded Section 2.2 (Lines 122-130) to include more details on sample pretreatment, processing, and the procedures used to separate FPs from other components, thereby clarifying the composition and handling of the fine-particle fraction.

At station TJ-S, our retrieved trap samples were divided into terrigenous particles, opal, carbonate, and organic carbon. The total mass flux (TMF) was calculated from sample dry weight (mg) normalized to the trap collection area and sampling duration, while POC and component-specific mass fluxes (including opal flux, carbonate flux, and terrigenous flux) were derived by the measured TMF percentage (%). Before zooplankton numeration, these fluxes were removed to ensure accurate quantification of fecal pellet contributions.

Comment #9:

Lines 155-156: Why the authors say spherical pellets are “slightly smaller”? What’s the criteria? Based on the average biovolume, it’s much smaller than cylindrical ones, but larger than the ellipsoidal ones (not similar). Please do statistical significance and present p values. 

Reply: We thank the reviewer for the careful examination and for pointing out this problem. We agree that, based on the average biovolume data, our previous statement was not precise.  In response, we have now conducted ANOVA tests to compare the biovolumes of fecal pellets with different shapes. The results indicate that cylindrical pellets were significantly larger than both elliptical (p < 0.001) and spherical (p < 0.001) pellets, while spherical pellets were also larger than elliptical pellets (p < 0.05). We have now corrected the expressions and added these results in our revised manuscript (Lines 192-194).

Comment #10:

Lines 176-178: I am confused here. The contributions of ellipsoidal, spherical and cylindrical to FPN are 48%, 38% and 15%, respectively, so how can we got > 90% for ellipsoidal + spherical?

Reply: We thank the reviewer for the careful examination and for pointing out the inconsistency. We apologize for the inaccurate expression in the original text. The values of 48 %, 38 % and 15 % represent annual averages, whereas during specific periods, the combined contribution of ellipsoidal and spherical pellets indeed exceeded 90 % (e.g., 93 % in UP01 during September, and 92 % in UP08 during January). To clarify this point and avoid confusion, we have revised this sentence into:  ‘During specific seasons, ellipsoidal and spherical pellets together accounted for over 90 % of total FPN.’ (Lines 227-228)

Comment #11:

Lines 244-249: What’s the relationship of these sentences and pellets export? I do not think they are necessary in discuss how phytoplankton blooms. I also suggest to re-organize this section to make it simple and clear.

Reply: We thank the reviewer for this constructive suggestion. We agree that, the discussions of the plankton bloom hypothesis was not directly relevant to fecal pellet export, and the logic of this section required improvement. In response, we have now removed redundant information unrelated to fecal pellet export (e.g., general theories of phytoplankton blooms), and reorganized this section to improve clarity, coherence and focus on the topic. (Lines 296-323)

Comment #12:

Line 252: “microbial remineralization and degradation”, please delete “and degradation”, because remineralization is generally equal to degradation.

Reply: We thank the reviewer for the careful examination. We agree with this suggestion, and we have now uniformly use the term ‘remineralization’ instead of ‘degradation’ in our revised manuscript. The revised sentence reads:  ‘Processes during sedimentation, such as microbial remineralization, can strongly regulate the transform efficiency but remain difficult to quantify due to the lack of downward sediment trap samples.’ (Lines 311-312)

Reviewer #2:

'Comment on egusphere-2025-2864', Anonymous Referee #2, 14 Sep 2025

The study by Wu et al., uses time-series data collected from a sediment trap deployed in the southern South China Sea to provide a year-round assessment of fecal pellet–mediated carbon fluxes and to evaluate how these fluxes respond to upper-ocean dynamics. The findings highlight that, in addition to the seasonal monsoon-driven pattern, other physical processes—such as typhoons, eddies, and episodes of high wind speed—also contribute significantly to the variability of deep fecal pellet carbon fluxes.

Overall, the manuscript is well written and logically organized. The figures are presented effectively, and the results strongly support the main points. I do not have major concerns. My suggestions are mainly minor, relating to figure design, methodological description, and expansion of the discussion to improve clarity and depth. Please see my detailed comments below:

Reply: We sincerely appreciate your effort reviewing our manuscript and your positive feedback on the manuscript. We have read through all your comments carefully and have made related modifications to improve the quality of our manuscript. We highly appreciate your time and consideration.

Specific Comments:

Comment #1:

Line 40: Please add relevant references after each study approach cited, to demonstrate how these methods have improved our understanding of fecal pellet–mediated carbon fluxes.

Reply: We thank the reviewer for this helpful suggestion. Accordingly, we have added the relevant references to support each study approach mentioned to strengthen this statement (Lines 33-41) (see Supplement for the revised manuscript annotated). We now have included references for study approaches including in situ observations (Shatova et al., 2012; Turner et al., 2015; Li et al., 2022; Wang et al., 2023; Cao et al., 2024; Darnis et al., 2024), satellite observations (Siegel et al., 2014), Bio-Argo float profiling (Estapa et al., 2017; Terrats et al., 2023), and numerical modeling (Stamieszkin et al., 2015; Countryman et al., 2022).

Comment #2:

Line 75: The abbreviation SCS has already been defined in a previous section; repetition is unnecessary.

Reply: We thank the reviewer for this helpful comment. We agree that this abbreviation is already discussed and is not necessarily presented. To avoid unnecessary repetition, we have now revised Section 2.1 by removing the redundant abbreviation of SCS and replacing it with a clearer description of the mooring station TJ-S (Line 81-85).

The revised text now reads: ‘Mooring TJ-S (6.72ºN, 110.76ºE, 1630 m water depth) is located in the southern SCS near the Sunda Shelf (Fig. 1a), where the East Asian monsoon system (EAM) exerts dominant control on local wind fields and ocean circulation (Shaw and Chao, 1994).’

Comment #3:

Lines 165–176: Cite the corresponding figures or tables when describing results.

Reply: We thank the reviewer for this helpful suggestion. We apologize for not citing corresponding figures or tables in our discussion sections. In response, we have revised the Results sections and added citations to corresponding figures and tables to improve clarity (Line 188-209).

Comment #4:

Line 130: Provide more details (including model data source, configuration, and validation) for the models used to derive the biogeochemical parameters (e.g., Global Ocean Eddy-Resolving Reanalysis, CMEMS Global Ocean Low and Mid Trophic Levels). Many biogeochemical models are available; please justify why these particular models were selected and comment on their performance in the SCS.

Reply: We thank the reviewer for this valuable suggestion. In response, we have expanded Section 2.3 to include additional details on the model and the data sources, including their configuration, data assimilation schemes, and their reported performance in the South China Sea (Lines 144-175). In addition, we have added a brief justification for selecting these particular products and cited recent studies that validate or employ them in this region (Lines 146-148).

Comment #5:

Figure 1: Consider embedding a small map showing the relative location of the SCS within the broader Pacific Ocean. This would help readers unfamiliar with the region.

Reply: Thank you for this helpful advice. We have now added a small global map showing the relative location of SCS in Figure 1 (Fig. 1b, Line 104).

Comment #6:

Line 80: The description that the mixed layer depth (MLD) is “relatively deep” in the SCS may not be appropriate. As shown in your results, there is a significant seasonal cycle, and the maximum MLD remains shallower than 60 m.

Reply: We thank the reviewer for this valuable advice. We agree with your concern that our previous description is not appropriate. To provide a more accurate summary of MLD variability in the southern SCS, we added relevant references (Qu et al., 2007; Thompson and Tkalich, 2014; Liang et al., 2019) and have rewrote this sentence accordingly. The new version reads: 

‘In this region, the mixed layer depth (MLD) is primarily controlled by air-sea heat fluxes and exhibits remarkable seasonal signals under monsoon forcing, deepening during June-August and December-February but remaining shallower than 60 m annually.’ (Line 90-93)

Comment #7:

Line 90: It is difficult to believe that diatoms contribute ~70% of the phytoplankton biomass in the open-ocean upper layer. Do you mean in the sunlit surface water, or are you referring to aggregates/sinking particles in the mesopelagic zone? Please clarify.

Reply: We thank the reviewer for pointing out this issue. We agree with your concern. The previously cited values refer to sunlit surface water contributions, and we apologize for the inaccurate expression in the original version. In response, we have added clarifications and relevant references (Lines 96-99) and revised the sentence as follows: ‘Diatoms and dinoflagellates dominate the surface phytoplankton community in the region, with contributions from picophytoplankton and cyanobacteria (Zhu et al., 2003; Ke et al., 2012, 2016; Wang et al., 2022).’ 

Comment #8:

Line 120: Please add methodological details on how POC flux was quantified from the sediment trap, since the FPC:POC flux ratio is discussed extensively later. In addition, the observation that the FPC:POC ratio remains relatively constant despite increases in total POC flux is intriguing. Does this reflect a true ecological signal—i.e., that different components contributing to POC flux increase proportionally—or could it be an artifact of the methodological approach? Some clarification on this point would strengthen the interpretation.

Reply: We thank the reviewer for this valuable comment. In response, we have expanded Section 2.2 to include additional methodological details on how POC flux was quantified from the sediment trap (Lines 126-130). Since POC flux and FPC flux are calculated independently, we infer that the relative constant FPC/POC ratio reflects a true ecological signal rather than methodological artifact, suggesting that different components contributing to POC flux increased proportionally during periods of increased pellet export.

Comment #9:

Line 125: The a in Chl a should be italicized. Please modify throughout the manuscript.

Comment #10:

Line 140: The p in p-value should be italicized. Please modify throughout the manuscript.

Reply: Thank you for this advice. We have corrected the term Chl a and p-value in the manuscript as suggested, and have ensured its consistency throughout the text. Specific corrections are listed as follows:

Chl a: Line 159, Line 345, Line 353, Line 445, Line 471, Line 488, Line 494.

p values: Line 180, Line 194, Line 195, Line 342, Line 343, Line 347, Line 348, Line 500, Line 501, Line 507.

Comment #11:

Figure 3: Use different color schemes to denote the two monsoon seasons to improve visualization.

Reply: Thank you for this advice. We have changed the color scheme of summer and winter monsoon in Figure 3 (Line 237), Figure 4a, 4b (Line 256), Figure 5 (Line 288), Figure 7k (Line 443), and Figure 8j (Line 492) to improve visualization.

Comment #12:

Figure 5a: Please clarify which direction the positive values represent.

Reply: We thank the reviewer for this helpful comment. We apologize for not adding clarifications of wind directions in the caption of Figure 5. Fig. 5a is the eastward component of 10 m wind, where positive values represent eastern winds and negative values for western winds. Similarly, Fig. 5b is the northward component, where positive values represent northern winds, and negative values for southern winds. To improve clarification, we have now revised the description of Figure 5 as follows: ‘(a) eastward component for 10 m wind (u), where positive values represent eastern winds, and negative values for western wind; (b) northward component for 10 m wind (v), where positive values represent northern winds, and negative values for southern winds.’ (Lines 289-291)

Comment #13:

Figure 5h: Report depth-integrated NPP in mg C m⁻² d⁻¹ to align with the units of POC flux. Also, the time unit is missing.

Precipitation data may be redundant and contribute little to your analysis; consider removing it.

Reply: Thank you for this helpful advice. We have changed the units of NPP in Figure 5j to align with POC flux, and added corresponding time units. We agree that the precipitation data is not used in our further discussion and analysis, and this data has been removed (Figure 5, Line 288).

Comment #14:

Figures 6g–6h: Place surface nitrate on the x-axis, since the aim is to examine how zooplankton respond to upper-ocean dynamics.

Reply: Thank you for this helpful advice. We agree that in our discussion, zooplankton biomass should be the dependent variable (y axis), since we aim to examine how they correspond to upper-ocean dynamics. Thus, the nutrients should be in x axis, as they represent independent variables. In response, we have revised Figure 6 and changed the x-axis of Figure 6a-6d to examine the respond of biochemical parameters to upper ocean dynamics. (Line 350)

Comment #15:

Line 295: Please clarify how you estimated the ~22-day time lag at your study site.

Reply: We thank the reviewer for this valuable comment. The estimated time lag between primary production and fecal pellet export at 500 m was inferred mainly from previous studies in the SCS (e.g., Bao et al., 2023), where biological responses were reported to occur within several days to weeks. In our study, each sediment-trap sample integrated over a 22-day period; therefore, we assumed that events occurring during one sampling interval (e.g., UP06) could influence both that sample and the subsequent one (e.g., UP07). This assumption guided our interpretation that WSP 2 and WSP 3 may correspond to PK 2 and PK 3. However, due to the lack of concurrent surface observations (e.g., phytoplankton or zooplankton analyses), the exact time lag at our station cannot be determined quantitatively.

Comment #16:

Line 335: Replace CMD with the more standard term deep chlorophyll maximum (DCM).

Reply: Thank you for this helpful advice. We have corrected the term deep chlorophyll maximum (DCM) in the manuscript as suggested, and have ensured its consistency throughout the text and figures.

Specific corrections are listed as follows:

DCM: Line 174, Line 293, Line 359, Line 362, Line 434, Line 446, Line 494, Line 525.

Comment #17:

Line 410: What is the typical fractional contribution of zooplankton fecal pellets to total POC flux in oligotrophic oceans? Additional discussion and comparison with previous studies would strengthen this section.

Reply: Thank you for this valuable comment. We agree with that including comparative discussion will strengthen our interpretation. Accordingly, we have reorganized this sentence as follows:

‘In the SCS, zooplankton fecal pellets make the most contribution to POC export in the southern region, with FPC/POC ratio ranging from 10.0 % to 42.6 %, reaching an average of 21.6 %. This value is larger than most oligotrophic regions including the central north Pacific (Wilson et al., 2008), the Sargasso Sea (Shatova et al., 2012), and the Mediterranean (Carroll et al., 1998), demonstrating the critical role of zooplankton fecal pellets in the unique carbon export process in the southern SCS.’ (Lines 530-535)