the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 28 Jul 2026
Contrasting physical and biological drivers of oxygen change in the northern California Current System
Abstract. Eastern Boundary Upwelling Systems like the California Current System are among the ocean’s most productive regions, shaped by distinct surface and subsurface currents and water masses. They are also vulnerable to deoxygenation, with several systems experiencing extreme low-oxygen events and associated mass mortality in recent decades. Determining potential drivers is complicated by seasonal dynamics, multiple water sources, and biological consumption. In the northern California Current System, oxygenated subarctic water mixes with oxygen depleted water from the south. We analyse three 20-year records, including one synthesized from local data, to provide valuable insights into subsurface oxygen variability. These long coastal time series are rare. We show that on the slope in spring, subsurface oxygen variations are physically driven by the relative contributions of these water masses, and their salinity and temperature ratios (a property called `spiciness') can be used to estimate oxygen levels. Seasonally, decreasing spiciness provides evidence of springtime oxygen renewal on the shelf, before respiration consumes oxygen over summer. Since 2003, mid-depth isopycnals show increased southern water contributions in spring and summer, leading to lower oxygen in spring, but not summer, suggesting summer-time oxygen consumption may have decreased with time. On the densest isopycnal near the shelf floor (~ 200 m), where oxygen is already lowest, no changes in spiciness or oxygen have been observed. Future research should examine the use of temperature and salinity ratios to estimate oxygen levels in the wider California Current System, and how changing source water conditions and biological activity may affect the risk of low oxygen in coastal waters.
- Preprint
(6522 KB) - Metadata XML
-
Supplement
(2246 KB) - BibTeX
- EndNote
Status: open (until 22 Oct 2026)
- RC1: 'Comment on egusphere-2026-4197', Anonymous Referee #1, 07 Sep 2026 reply
-
RC2: 'Comment on egusphere-2026-4197', Anonymous Referee #2, 15 Oct 2026
reply
General Summary
This manuscript presents an analysis of the subsurface oxygen variability over a 20-year period in the northern California Current system. The authors demonstrate how physical advection (tracked by spiciness) and local biological consumption drive local oxygen concentrations in different seasons and in different isopycnal surfaces. Results suggest that seasonal oxygen variations on the shelf and slope environments are primarily driven by the relative contributions of subarctic and southern/subtropical source waters, but oxygen signatures get further modified by local subsurface respiration. While the overall premise of the study is an interesting and novel contribution, I believe the analysis would benefit from additional quantitative analyses to support some of the claims. Specifically, I think accounting for oxygen trends in source waters and incorporating Apparent Oxygen Utilization (AOU) can help clarify if some of the signals are indeed due to locally enhanced or decreased subsurface respiration. I believe this manuscript would be suitable for publication after additional complementary analyses and adjustments to some of the figures to better depict uncertainties.
Main Comments
Seasonality of oxygen in the northern CCS
- Lines 121-135: The authors point to the lack of oxygen increase on the shelf on the 26.0 and 26.3 isopycnals as a respiratory signal in the early EUW season. This is accompanied by a decrease in spiciness, which would lead to an expected increase in oxygen. The hypothesis here is that this lack of oxygen increase is driven by enhanced local respiration, perhaps due to an increased delivery of nutrients onto the shelf during upwelling/isopycnal shoaling stimulating a corresponding increase in primary production and oxygen demand during remineralization.
To explore this hypothesis more fully, the authors should consider incorporating Apparent Oxygen Utilization (AOU) to track if oxygen consumption is indeed increasing seasonally during the transition from LDW to EUW. AOU can be easily calculated from the variables already presented. If discrete nutrient bottle data (e.g., nitrate) are available for these time series sites, then this could further corroborate your hypothesis. Any AOU anomaly should be stoichiometrically mirrored by a nitrate anomaly, which would provide evidence that the lack of oxygen increase is biologically driven rather than an artifact from mixing with a third water mass.
Trends in spiciness and oxygen, 2003-2023
I think the study would benefit from more robust quantitative analysis to support the statements about observed trends (or lack thereof) in seasonal subsurface oxygen variability on mid-depth isopycnals being caused by changes in respiration over time.
- Lines 154-159: The authors point to the oxygen decline in upper isopycnals during EUW is related to a change in relative source water contributions, indicated by an increase in spice on the shelf and slope. Because the shelf oxygen is not fully compensated by the mixing expected from the spice relationship, an increase in local respiration during EUW is again hypothesized as an important driver of oxygen loss. This deviation from expectation is calculated using the static spice-oxygen relationship presented in Figure 2. However, as shown in the source water timeseries (Figure S6), the subarctic end-member (P26) shows significant decline in baseline oxygen in the upper two isopycnals (Table S3). Because the source water is arriving at the study site pre-conditioned with less oxygen over the 20-year period, the empirical mixing line will underestimate dissolved oxygen decline. In this case, source water deoxygenation must be factored into the “expected” oxygen calculation before the deviation can be appropriately attributed to local biological drawdown. To ensure this is a robust local signature and not an artifact of offshore deoxygenation, the authors could define the oxygen-spice relationship in discrete time bins (e.g., 5-year windows) and compute time-varying residuals as deviations from the expected relationship. Further, similarly computing trends in AOU and oxygen saturation for the upper isopycnals could corroborate local changes and explicitly separate the physical solubility drivers from biological drivers. For example, if there is a robust trend in AOU and not in oxygen saturation, then the authors can conclude there is a biological influence.
- Lines 160-172: The analysis in the LUW season, the authors hypothesize a decrease in subsurface respiration over time based on an insignificant trend in oxygen with the corresponding significant trend in spice. However this is slightly contradicted by Line 170 and Figure S10, where the trends in seasonal oxygen consumption between spring and summer are not statistically significant. I would suggest removing this hypothesis, or alternatively exploring whether this is actually a phenological shift on biological oxygen demand. If primary production is shifting earlier in the year (into the EUW window), the earlier consumption of organic matter would leave less available for respiration later in the season during LUW. This shift in seasonal timing could help explain both the EUW increase in oxygen decline and LUW lack of oxygen decline on observations without requiring a long-term decrease in subsurface respiration.
Specific Comments
- Line 25: most productive compared to the rest of the CCS?
- Line 31 (and throughout): an equal sign is missing after σt
- Figure 1: Could you make the study area box with the black outline a different color? It was a little hard to see at first.
- Line 68-69: Can you clarify what you mean in this statement “When data exist in every year, there are 20 data points in each seasonal, regional timeseries.”? Figure 1c caption indicates that data from early downwelling are too sparse to be included, and Figure S7 does not support 20 data points per (dynamical) season. Can you reconcile this?
- Line 70: Please also provide a physical definition of spiciness in addition “spicy” and “minty”. It could be useful to introduce in reference to track the mixing of different water masses along a constant density surface.
- Figure 2: I suggest you use another color for CalCOFI instead of red to accommodate red/green colorblindness.
- Figure 3: Can you clarify what you mean by “slightly offset for legibility”?
- Figures 4, S4, S5, S6, S9, S10: The confidence intervals look incorrectly plotted. The trend lines should fully be within the confidence intervals.
- Figure S1: Panels (d) and (e) show the same plot, just (e) with shading. Panel (d) can be removed. Similar to other figures, the confidence interval shading in panel (e) looks incorrect.
Citation: https://doi.org/10.5194/egusphere-2026-4197-RC2 - Lines 121-135: The authors point to the lack of oxygen increase on the shelf on the 26.0 and 26.3 isopycnals as a respiratory signal in the early EUW season. This is accompanied by a decrease in spiciness, which would lead to an expected increase in oxygen. The hypothesis here is that this lack of oxygen increase is driven by enhanced local respiration, perhaps due to an increased delivery of nutrients onto the shelf during upwelling/isopycnal shoaling stimulating a corresponding increase in primary production and oxygen demand during remineralization.
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 212 | 115 | 37 | 364 | 77 | 45 | 35 |
- HTML: 212
- PDF: 115
- XML: 37
- Total: 364
- Supplement: 77
- BibTeX: 45
- EndNote: 35
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This manuscript investigates the drivers of oxygen changes in the coastal northern Pacific using a novel approach based on the spiciness-oxygen relationship. The method and results are clearly laid-out. This contribution will be suitable for publication once further validation of some of the results is provided.
Main comments
1. In a recent publication by the same authors, the temporal evolution in spiciness is discussed in much the same way as in the present manuscript. Most of the figures in Maier et al (2025) are reproduced in the supplementary material of the present manuscript or as panels of the main figures. Since these results have already been published, I recommend that for these figures the authors refer to the published paper instead of reproducing them. Part of the results and of the discussion of trends could probably also be removed from the present manuscript.
2. The analysis is based on annually averaged observations. Two pieces of information are missing to support the robustness of the results: the spread of the values underlying each annual average (standard deviation), and the number of observations contributing to it. If the spread is much larger than the obtained trends, the trends may not be significant. Without this information, it is also not possible to assess whether there is a strong spatial or temporal variability within the considered regions and time frames, which could affect the results if there are biases in the observations (towards specific locations or dates in certain years). Therefore, I suggest that the authors:
(i) compute the standard deviations associated with the annual averages, and possible the associated uncertainty by also taking into account the number of observations;
(ii) display them as error bars on figure 2, 3, 4;
(iii) modify the statistical tests so that they account for these standard deviations. For instance, one approach is to weight each data point by the inverse of its variance.
These validations appear essential to confirm the validity of the results.
3. While most of the analysis are qualitative, the proposed method seems to offer possibilities for more a quantitative treatment.
First, the mixing ratio of the two parent water masses could be quantified from the distance along the spiciness axis of the mixing line (in Fig. 1).
Second, I believe that the vertical offset in oxygen between individual points and the mixing line would provide a quantitative estimate of departure from parental water properties, and hence of respiration.
Third, as discussed at L186-187 referring to other studies, a trend in parent water masses properties would affect the results. These trends could be quantified based on the dataset used here. I suggest adding a figure showing these trends possibly in the supplementary material, and then discussing the effect of the trends on the results. Maybe such trends could provide hints into the unexplained trends discussed at L186-187.
Then, I would encourage the authors suggest to use these quantitative analysis to provide more quantitative results in the abstract and conclusion, instead of general conclusions such as the one at L201.
4. The proposed method assumes that mixing is isopycnal. I suggest that the authors discuss the potential role of diapycnal mixing, as mixing could change the properties as waters move between the slope and the shelf. For instance, the following studies seems to indicate a significant role for diffusion and vertical mixing in shelf-sea interactions in that region: https://doi.org/10.1029/2010JC006720, https://doi.org/10.1080/07055900.1989.9649345.
5. I suggest deepening the interpretation of the results (L93-100, L147, L163) proposing a physical/climatic hypothesis to explain the observed change in water mass composition, for instance in relation to known changes in circulation, ventilation, or other mechanisms at play in the region. Hypothesis are provided for the change in respiration, but not in composition. Relatedly, multi-annual variability appears to be visible on Fig. 4 in addition to the long-term trends. It could be acknowledged, and suggests a role for natural or climatic variability.
6. I would like to clarify one part of the discussion. In Table S3, the spiciness trend is weaker in LUW than in EUW, as also mentioned at L166-167. However, Fig. 4 shows the opposite. Please clarify this inconsistency. From the figure, the difference in the spiciness trend between the two appears significant, and would reduce the O2 trend, bringing it close to the level of detectability, as observed.
Specific comments
Throughout, these appears to be a formatting problem with the density, with a missing equal sign or space. See for example L31.
L7-10 : The text from “We show” to “over summer” would benefit from clarification. It is not clear what the authors mean by “Seasonally” in the second sentence.
L25: Please clarify the northern CCS are the most productive in terms of what and compared to what.
L26: Please replace “extreme low” with a concentration.
L26: “a mix of” may be more accurate than “subject to”.
Fig. 1 inset: a colorbar is missing for spiciness
Fig. 1 caption: A colon may be needed so that the text reads “Inset: mean ...”
L61 : The reader wonders if the depths provided in parenthesis are for a specific season. Actually, I suggest adding a figure showing density and oxygen profiles in the four locations. This would help the reader interpret the water mass dynamics.
Fig. 2: I am not convinced that panels b-c, e-f, and h-i are necessary.
L132: I suggest comparing these rates with existing estimates, and commenting on whether they are realistic.
L134: Is the higher productivity related to terrestrial nutrient input ? If yes, this could be mentioned.
L136: This first sentence is strongly redundant with sentence “However...” at L137.
Fig. 3 caption: Please clarify what is meant by “dots slightly offset for legibility”. Is it difficult to distinguish offsets from genuine differences.
L160-165: These sentences are difficult to follow; I suggest rewriting them for clarity.
L161: The comma after the closing parenthesis appears unnecessary.
Fig. 4: For the shading showing the bootstrap confidence intervals, a bootstrap interval should contain the trend line. Here, the trend line lies at one edge of the shaded range or even outside it. Please check the calculation or the plotting. A further indication of a possible problem is that the shading appears to be symmetrical about the zero-slope line rather than about the fitted trend. This would normally imply that the trends are not statistically significant at the corresponding confidence level.
L192: The same result would, I believe, also be obtained if the contributions of the two water masses changed by equal amounts in opposite directions.
Fig. S2: Panels d and e are redundant, I suggest removing panel d.
Fig. S10: I am always uncomfortable when trends that are not statistically significant are displayed on figures. I would suggest removing the trends from this figure, or maybe showing them with a very faint dotted line.