Review of: ‘Unexpected quasi-independence of colored dissolved organic matter

absorption from chlorophyll-! concentration in the Southern Ocean’ by Li et al.

Reviewer: Emmanuel Boss, UMaine.

This paper investigates the relationship between CDOM (the absorbing part of DOM) and chlorophyll using BGC-Argo data. A simple bio-optical model to decompose the contributions to the diffuse attenuation of the measured irradiance is used with additional inputs from the chlorophyll fluorometer and backscattering sensors on the float. For each wavelength a CDOM estimate is computed and compared to chlorophyll to investigate how they relate both in the SO as well as further north.

The subject of the paper is very interesting and pertinent as there is currently a debate regarding the skill of remote sensing retrievals in the SO. Understanding why they may not work, such as if assumption in them are wrong, would help advance the science we can do from remote sensing in the SO by removing biases.

The paper is relatively clearly written though the notation/formatting could be significantly improved to avoid confusion and make it easier to read (for example CDOM, the subject of the paper, is written as a_y, ay and a_y, see attached pdf). Such sloppiness in the notation of the subject of the paper does not reflect well on the author (two of which I know well and are excellent and careful scientists).

I have some significant comments that I feel, if addressed, will significantly improve this manuscript. I am also returning an annotated PDF.

1. I feel like there is no real closure in this paper. This could be done, for example by showing the degree of agreement between the different wavelengths when computing an a_y spectra. Since each radiometer is an independent sensor, it provides you an opportunity to test the quality of your retrieval for each profile. You only provide us the one estimate of the spectral slope of a_y for each region. It would be very informative to see it for each profile. It could be used by you as an independent quality check as CDOM slopes are well constraint in the literature. You could, for example, provide a histogram of a_y(380)/a_y(412), a_y(380)/a_y(490), and a_y(412)/a_y(490) (or their reciprocal), and compare to what is expected. You are making so many assumption (e.g. a_p to chl relationship, chl to Fchl relationship etc’) that it is hard to evaluate if your uncertainties are indeed as small as you think they are. This closure exercise will help determine it.
2. Alternatively, you could use kd(490) and Fchl to constrain better chlorophyll (using a Xing like method) and then see if a_y(380) and a_y(412) are better related.
3. The adjustment done to the chlorophyll and bbp is not well justified. It is hand waved. In addition, at low chl, Fchl could be significantly affected by CDOM (see another Xing paper) biasing your adjustment.
4. Table 1 – not clear if you assume no uncertainties at 412 and 490 when they are not provided.
5. You mix Discussion (or interpretation) into your Result while you provide two separate sections. I suggest, in particular because of the high uncertainties of your assumptions, that you clearly separate the two. Have the Result section be an objective description of what you found and use only the Discussion for speculation about the cause for your observation. Also, you may want to add a Summary section where you summarize the most important implication of what you found and suggest future work.
6. It is not clear what Fig. 2 represent and how it is linked to Eq. 4. I suspect it is the relative contribution of different terms to (a+bb), but this is not what is written.
7. The discussion relating Fig. 8 to DOC in the upper ocean is wrong. There is actually more DOC in low latitude compared to high latitudes (see Hansel, 2009, Oceanography). And, generally, CDOM and DOC do not correlate in the open ocean.
8. The discussion regarding case 1 and case 2 could be more nuanced. Bricaud and Morel have a paper from 1981 showing decorrelation between chl and a_y. When I asked Morel he told me that locally one may see a lack of correlation but when one uses a large dynamic range from many contrasting oceanic environment, a general relationship may emerge. When looking at your data I see a huge spread around your best fit lines, with the lines being a poor representation of most of the data (e.g. if I had to predict a_y from [chl], even not in the SO, the likely error could be very large, to the point of being useless for many applications. Maybe if you used a heat map for your regression plot (with color being the concentration of datapoint of a given value one could see better how the relationship derived relate to the statistical distribution of the points and appreciate that the fits are much better than I claim here).

Dear authors: I am often wrong. If you feel that my comments are off the mark, please contact me and I would be happy to change them if convinced.

Authors present analysis of differences in bio-optical properties between the low-latitude ocean and the Southern Ocean (SO) using many 1000's of BioARGO float profiles which include multispectral radiometers. Their focus is on relationships between colored dissolved organic matter (CDOM or a_y) and chlorophyll retrievals and show a rough independence between these two properties. From what I can tell, the analysis of deriving a_y spectra seems valid as supported by the sensitivity analyses presented. There is a lack of statistical rigor in the analyses presented and the discussion presented is neither insightful nor useful. I think that the paper could be a fine contribution to the literature once a couple of issues are resolved. These issues are delineated in the text that follows.

I found section 3.1 very confusing. First, I think the description of figure 2 is incorrect. The text and caption both state that what is plotted is the relative contributions of optical properties to equation 4, including a_y.  But, the right hand side of equation 4 is a_y. So, this does not make any sense. I think what is plotted in figure 2 are the relative contributions that the various optical properties make to K_d.

Figure 3 is fine, showing invariance of several of the approaches to the retrieved a_y value. However, I do not understand what Table 2 is trying to communicate. Line 230 states that "K_d is the main contributor to differences in a_y" and cites Table 2, but I do not see how the "dispersion of the mean a_y" from the sensitivity analysis shown in Table 2 demonstrates that.

Section 3.2 shows differences in the K_d vs. Chl relationships from the low-latitude ocean and the Southern Ocean. The analysis presented is completely qualitative and needs to be more quantitative. I suggest that the authors assess whether there are statistically significant differences among the various relationships presented. The colors selected in figure 5 are very difficult to see (especially the yellow and white text).

Section 3.3 (Figure 6) gets to differences in the relationships between a_y spectra and Chl concentrations. As in the last section, there is a lack of statistical rigor and the colors selected for the figure needs changing. The authors need to show that the differences commented upon are actually statistically significant. Most importantly, a test of independence between a_y and Chl would be that the derived power-law slope is not statistically different than zero. This result would prove the author's hypothesis in a compelling way.

As mentioned above, the discussion is very weak and does not get anywhere useful. The existing sections need to be further elaborated upon and expanded new points to be introduced. One suggestion would be to address previous work on why the SO bio-optical properties differ and satellite Chl values are underestimated compared with lower latitudes. Some researchers thought it was unusually low pigment normalized phytoplankton absorption spectra (Mitchell and others) or glacial flour changing particle backscattering (Dierssen) or differences in photochemical processes. This would help round out section 4.2. Also, you should read through Yamamoto et al. JGR (2024) which does a global analysis of CDOM sources and sinks that would help think about the processes that are creating and destroying CDOM and why the SO might be different.

Last, this is not the first time that Southern Ocean CDOM patterns have been addressed (see Ortega-Retuerta et al. 2010, there are likely others).