the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Building your own mountain: the effects, limits, and drawbacks of cold-water coral ecosystem engineering
Anna-Selma van der Kaaden
Sandra R. Maier
Siluo Chen
Laurence H. De Clippele
Evert de Froe
Theo Gerkema
Johan van de Koppel
Furu Mienis
Christian Mohn
Max Rietkerk
Karline Soetaert
Dick van Oevelen
Download
- Final revised paper (published on 22 Feb 2024)
- Preprint (discussion started on 30 May 2023)
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2023-949', Anonymous Referee #1, 21 Jun 2023
The approach of this paper fits in line with several other publications/manuscripts of this group in recent years. It applies statistics to morphological and environmental parameters trying to explain the distribution of cold-water corals on coral mounds or coral mounds in general, in this case referring to the Logachev Mounds at the flank of the Rockall Bank. Actually, most of the data presented here have been published before by these authors, the new aspect here is a (another) different statistical approach to evaluate the data. And I have to admit right away that I am not a specialist in statistics and that it is very difficult for me (if possible at all) to follow the statistical approach applied here (which presumably is due to my ignorance on this topic). That means, I cannot assess the credibility of the statistical work presented here. Nevertheless, probably I am a common “consumer” for such a study with a reasonable background in cold-water coral research that enables me to assess the interpretations and the conclusions presented here.
For many years we already see the self-enhancing effect of coral mound formation. Thus, I applaud the authors for picking this up and trying to go into the details how this process works, which so far was simply seen as a direct link between increasing mound size results in stronger hydrodynamics that result in a better food supply to the corals – with these processes affecting in particular the mound summits. This presence of this general pattern is also described in the conclusions of this ms (although Fig. 2 shows that it does not fully apply to all the mounds investigated here). What is new, is the finding presented here that reef growth is controlled by downward water velocities in 250-300 m water depth in February. Actually, that these downward (and upward) fluxes occur around the mounds has already been presented before (e.g., Sotaert et al. 2016). However, that these control reef growth seems to be a clear overinterpretation of the data. Although this parameter is the strongest contributor to RDA axis 1 (see Fig. 4), it still is only one among several others with a similar weight. Furthermore, RDA axis 1 only explains 27% of the variance, with a total of ~55% being not explained at all by the parameters considered. So, there are probably other parameters involved and the available data do not justify to give the overall control to the downward velocities in February. The mechanistic explanation behind this “control” is seen in “deep winter mixing and deep-water cascading triggering the downward flux of POM (no reference for the POM link!) which is strongest near the winter-spring transition (here: data for February), which fits to the early onset of the phytoplankton bloom (March/April) at Rockall Bank”. Mohn and White (2007), referenced for the latter statement, however, show that the spring bloom over Rockall Bank begins in May. Thus, does this mean that the downward fluxes in February predominantly transport “old” POM to the corals or no POM, thereby exerting “strong control” on the coral reef cover. How would this help the corals during the food-limited winter months as stated in line 328? And what happens later in spring, under spring bloom conditions, that most likely provide high amounts of fresh POM. How do the corals react on this? Here it becomes obvious that other parameters (some are mentioned on page 16) need to be included in this study to get a better grip on the controls of benthic cover distribution (and to explain more of the remaining 55% of the variance). Seeing these points, a hydrodynamically stimulated food supply in winter (whereas the model only shows vertical velocities and gives no hint on food supply) as a major factor controlling reef cover (see line 414) is not evidenced by these data.
Consequently, in its present stage this manuscript does not provide any new, well documented information on the controls on benthic cover in the Logachev Mounds. Thus, this ms definitely needs to consider further parameters, add more process understanding, provide more information about the meaning of its MEMs, and discuss its new ideas (rather than just backing them up by references that do not really fit).
Further points to consider are listed below:
Coral distribution: In line 193 it is stated that the live coral coverage typically increases towards the mound summit. The use of the term “typically” already hints to a somewhat weak signal. And indeed, in Fig. 2 in transect 1, highest live coral coverage is at the lower slope and in transects 2 and 5 it is along the upper slopes. In transects 3, 4, and 6 live coral coverage is continuously very low with maybe some more corals at the summit in transect 3. Thus, a clear signal reflecting the statement above is only seen in transect 7, thus, 1 out of 7. One might include here transects 2 and 5 with highest coral coverage along the upper slopes, but these do not show any further increase towards the summit. So, basically, the statement repeated above is somewhat misleading – as already seen by the authors letting them put in “typically”. In the conclusions, it is again referred to that only mound summits are topped with dense live corals (line 395).
Interestingly, the RDA reveals that the most important single variable to explain high vs. low reef cover between 600 and 1100 m water depth is the “February downward velocity between 250 m and 300 m water depth”. In addition, for the 2nd RDA axis, differentiating areas with high vs. low coral rubble cover, the spatial variable within-transect MEM3 dominates. The 4th RDA axis, finally differentiating between live and dead coral cover, however, only explains <1% of the variability and is driven by within-transect MEMs 6 and 7 and by “upward velocities in February”. Thus, in addition to the up- and downward velocities, the MEMs are important variables controlling the benthic cover. However, besides one very vague sentence (line 299-301), no hint is given in which way the MEMs exert any control. If these are so important, a clearer definition of the MEMs and an in-depth discussion how (!) these affect the benthic cover would be needed. Fig. A1 is of no help in this context: “colours and bubble size represent the values belonging to the individual annotated video frames” – which values? which of the 1200 video frames? how do colours and bubble sizes related to any numbers? I am sorry, but I cannot read anything out of this figure. This figure is also referred to in line 350ff, where in combination with Fig. 4 evidence is presented that high coral reef cover is driven by MEM 3 and 4 and by strong bottom currents. However, Fig. 4 points to downwelling in February, MEM2 and depth as most prominent parameters. It also does not become clear to what extent MEM3 and 4 are linked to specific flanks of the mounds … This part is quite confusing.
This continues in lines 354 ff. While the first sentence is not to follow as the described settings are not obvious in Fig. 4, it does not become clear at all, why these suggest strong upwelling around the mound flanks … Most of the contribution of the parameters described in the first sentence to the RDA axes are really minor. Thus, to me also this reads as a clear overinterpretation of the data. And the following sentence (line 356) is pure speculation (especially as the last two references only refer to the benefit high primary productivity can have for the corals).
Mound engineering effect: a key finding here is that the mound engineering effect results in a deceleration on mound summits (line 256), which is picked up again in the discussion (line 332). So, why do the corals prefer the mound summits? Available current meter data usually show highest current speeds at or close to the summit. So, what is the meaning of your assessment that the mound engineering leads to reduced current speeds at the summits? And in Fig. 7a it is shown that bottom current speed significantly increases with mound height. But this does not refer to the summits? From line 333 this continues with a description of a hydrodynamic regime on coral mounds also seen in other studies. However, neither in Genin et al. nor in Soetaert et al., this pattern becomes obvious. See, for instance, points 1) and 2), for which neither of the two other paper present any evidence, and the same applies partly to the other points. These are to see to some extent in the modelling study by van der Kaaden et al. (2021), however, field observations do not show such a pattern (e.g., Genin et al. 1986, Dorschel et al., 2007). So, the pattern described here might be introduced as a new conceptual idea, but then it needs more discussion about the involved processes, and not simply the hint to some other papers (which actually do not show this pattern).
One comment on the methodology: what does not become clear from the text, is how the hydrodynamic model data with grid cells of 250 m in the horizontal are combined with the video frame data. According to Fig. 2, the individual transects are between ~500 and ~1500 m long; only transect 3 is much longer with ~3100 m, but has a less dense video frame density. Thus, most transects correspond to 2 to 6 grid cells in the hydrodynamic model, with seeing the small-scale typical summit and slope settings probably being partly very difficult to discern within 250 m grid cells. On the contrary, the video frame density usually is much, much higher, with a total of ~1200 frames. Thus, ~1200 frames, which in parts show quite some variability along a transect, are linked to ~30 grid cells, that seemingly cannot always discern between summits, slopes, etc. – how this is done becomes – at least for me – not clear from chapter 2.3. Chapter 2.2.2 explains that for each frame the nearest grid cell has been chosen for correlation, so, there might be a few more than the 30 grid cells mentioned above. But nevertheless, the question remains …
Obviously, the authors like their own phrases very much. For instance, the hint to better conditions for the corals “higher up in the water column” occurs at least 4 times, e.g., with the statement that the corals would benefit from being higher in the water column as there, POM quality and quantity is higher (only giving a reference for the POM distribution, but not for the link to the corals). This is followed up at the beginning of the discussion where coral performance again is linked to water depth – described as a result of this study. However, these statements are in contrast to one of the main conclusions referring to the “Massenerhebung effect” stating that not the position in the water column, but the position on the mound is a key parameter controlling coral growth. In addition, also the title of the ms is more a nice wording than a good description of the content of the paper. While the word “mountain” in the first part is only picked up in the paper with reference to terrestrial mountains, the second part of the title goes far beyond its actual content.
Table 1 caption: which colours?
Line 177: with (!) coral mounds
Citation: https://doi.org/10.5194/egusphere-2023-949-RC1 -
AC2: 'Reply on RC1', Anna van der Kaaden, 16 Aug 2023
Dear reviewer 1,
Thank you for your review of our manuscript. From your comments, we infer that, unfortunately, our main conclusions did not come across clearly enough. We never intended to state that downward velocities in winter resulting from broad-scale non-engineered processes are the only factor controlling reef growth. We added to the revised abstract (line 24): “We find that downward velocities in winter, related to non-engineered environmental factors, e.g., deep winter mixing and dense water cascading, cause substantial differences in reef cover at the broadest spatial scale (20-30 km). Such hydrodynamic processes that stimulate the food supply towards the corals in winter seem more important for the reefs than cold-water coral mound engineering or than similar hydrodynamic processes in summer. While the ecosystem engineering effect of cold-water corals is frequently discussed, our results highlight also the importance of non-engineered environmental processes.
Therefore, we conclude that (line 443): “We showed that downward velocities in winter, likely driven by non-engineered processes as deep winter mixing and dense water cascading, control reef growth to such an extent that reef cover differs substantially at the scale of the Logachev cold-water coral mound province.”
And that (line 457): “When studying cold-water coral ecosystems and their response to global change the effects of broad-scale, non-engineered environmental processes should not be underestimated, and it is especially important to gain a better understanding of how benthic-pelagic coupling will change in winter and to quantify how much alternative food supply mechanisms and coral mound engineering could compensate a lower food availability.”
That hydrodynamics are important for cold-water corals has been shown, but we extend this knowledge by showing whether it is the engineered hydrodynamics or non-engineered hydrodynamics that exert the major influence on cold-water corals. Moreover, in contrast with earlier studies (e.g., Soetaert et al., 2016), we showed the statistical correlation between cold-water coral cover and hydrodynamic variables. Our analysis with live and dead coral framework further integrates long-term processes, whereas Soetaert et al. (2016) studied processes on the relatively short term (i.e., tidal processes and spring-neap tide).
We added a paragraph (4.3) on the mechanisms by which the hydrodynamic zones might affect cold-water coral growth.
Indeed, we explain ‘only’ 45.1% of the variance in benthic cover, which actually means that our RDA performed well in comparison to similar deep-sea studies. The results of an RDA are only reliable when the explanatory variables are not too much correlated to one another. Therefore, we used the forward selection procedure to identify the most important explanatory variables contributing to explaining the variation in benthic cover. Adding more variables would increase the percentage of variance explained but would make the analysis unreliable. Other variables that might influence benthic cover are anthropogenic influences, small-scale differences in environmental variables, interspecific interactions, or random effects, but these could unfortunately not be included in this study.
To the discussion (line 294) we added: “By explaining 45.1% of the variation in benthic cover, our set of variables performed well in comparison to similar deep-sea studies using RDA: e.g., Vad et al. (2020) explained 15.1% of the variation in benthic cover and Kazanidis et al. (2021) explained 26.7% of macrofauna variation. Henry et al. (2013) explained 65% of the variation in reef fauna using a combination of MEMs, hydrography, and bathymetric variables. So, the inclusion of both spatial factors (MEMs) and hydrodynamic variables appears to benefit the analysis. Especially within-transect MEMs described a relatively large part of the variation in our benthic cover data that was not explained by the environmental variables. Henry et al. (2013) similarly found that community composition was structured mainly by broad-scale environmental forcings, but that fine-scale MEMs described a significant fraction of the variation in community composition. This within-transect variation in benthic cover can be related to spatial patterns of fauna (Henry et al., 2013), such as self-organized regular spatial patterns in coral reefs (van der Kaaden et al., in press) and to the typical reef zonation on cold-water coral mounds (section 4.2). “
Mohn and White (2007) describe Chl-a levels over the Rockall Bank that are enhanced throughout the year and remain high during winter as compared to the surrounding water. In their figure 3, Chl-a levels over the Rockall Bank are higher in March than in February. Duineveld et al. (2007) further report that a fresh phytoplankton signal reaches the cold-water coral mounds at the end of February, signalling the start of the phytoplankton bloom. We added the reference to Duineveld et al. (2007) and changed the sentence to (line 327) to: “Both processes cause an accelerated downward flow of organic matter, that peaks near the winter-spring transition, at the onset of the phytoplankton bloom that is particularly early (March/April) at the Rockall Bank (Mohn and White, 2007; Duineveld et al., 2007).”
Below, I will reply in detail to your comments, displaying your comments in italics followed by our reply.
- Further points to consider are listed below:
Coral distribution: In line 193 it is stated that the live coral coverage typically increases towards the mound summit. The use of the term “typically” already hints to a somewhat weak signal. And indeed, in Fig. 2 in transect 1, highest live coral coverage is at the lower slope and in transects 2 and 5 it is along the upper slopes. In transects 3, 4, and 6 live coral coverage is continuously very low with maybe some more corals at the summit in transect 3. Thus, a clear signal reflecting the statement above is only seen in transect 7, thus, 1 out of 7. One might include here transects 2 and 5 with highest coral coverage along the upper slopes, but these do not show any further increase towards the summit. So, basically, the statement repeated above is somewhat misleading – as already seen by the authors letting them put in “typically”. In the conclusions, it is again referred to that only mound summits are topped with dense live corals (line 395).
We phrased the sentence as we did because we thought that clear trends were not obvious from the raw data. We now rephrased the sentence to (line 190): “There is variability in benthic cover within- and between transects (Fig. 2), but clear trends are not obvious. For instance, transect 3, located on the lower southern flank of the northernmost mound, and transect 6, located at the western flank of the highest mound (Fig. 1), have a relatively low live and dead coral cover overall. Transect 7, across the highest mound in the area, shows a relatively high dead coral cover and low sediment and rubble cover overall.”
We also rephrased the sentence in the conclusion (line 450): “…, leading to a typical reef zonation whereby mound upper flanks and summits have the densest reef cover.”
- Interestingly, the RDA reveals that the most important single variable to explain high vs. low reef cover between 600 and 1100 m water depth is the “February downward velocity between 250 m and 300 m water depth”. In addition, for the 2nd RDA axis, differentiating areas with high vs. low coral rubble cover, the spatial variable within-transect MEM3 dominates. The 4th RDA axis, finally differentiating between live and dead coral cover, however, only explains <1% of the variability and is driven by within-transect MEMs 6 and 7 and by “upward velocities in February”. Thus, in addition to the up- and downward velocities, the MEMs are important variables controlling the benthic cover. However, besides one very vague sentence (line 299-301), no hint is given in which way the MEMs exert any control. If these are so important, a clearer definition of the MEMs and an in-depth discussion how (!) these affect the benthic cover would be needed. Fig. A1 is of no help in this context: “colours and bubble size represent the values belonging to the individual annotated video frames” – which values? which of the 1200 video frames? how do colours and bubble sizes related to any numbers? I am sorry, but I cannot read anything out of this figure. This figure is also referred to in line 350ff, where in combination with Fig. 4 evidence is presented that high coral reef cover is driven by MEM 3 and 4 and by strong bottom currents. However, Fig. 4 points to downwelling in February, MEM2 and depth as most prominent parameters. It also does not become clear to what extent MEM3 and 4 are linked to specific flanks of the mounds … This part is quite confusing.
MEMs merely describe a certain spatial scale of spatial separation between the observation points, so they do not exert any control on reef cover. We clarified this by changing section 2.2.3 (line 136): “In the MEM-approach, the study area is divided into groups (i.e., a map with 1 group, 2 groups, 3 groups, etc.) and each annotated video frame is assigned a value that represents its proximity to the centre of the group in which it is located (Borcard et al., 2011). Moran’s Eigenvector Map 1 thus represents variations on the broadest spatial scale of the entire study area (~25 km) and the highest MEMs (i.e., when the study area is divided in the highest number of groups) describe variation at the finest spatial scale of the difference between individual annotated video frames (~4m). We visualized the MEM-procedure in Fig. A1 (Appendix A). MEMs merely describe a certain spatial scale and thus do not exert any control on CWC growth, their specific values are therefore meaningless with respect to CWC growth.”
And (line 146): “At the broadest spatial scale (signalling differences between transects) datapoints were classified according to their east-to-west position (MEM 1) or according to their north-to-south position (MEM 2). At finer scales (signalling differences within transects) datapoints were classified by MEMs according to their position along a mound flank (MEM 3, 4, and 5), or according to regular intervals along transects (MEM 6 to MEM 15).”
To clarify the discussion section 4.2, we changed the sentence in line 360: “…we found that high cold-water coral reef and rubble cover correlated to the spatial scale from mound foot-to-summit (MEM 3 and 4)…”. We further removed the reference to the supplementary figure in the discussion.
We will clarify the supplementary figure illustrating the MEM values: showing the values of MEM 1 to 7, along three example transects for each annotated video frame.- This continues in lines 354 ff. While the first sentence is not to follow as the described settings are not obvious in Fig. 4, it does not become clear at all, why these suggest strong upwelling around the mound flanks … Most of the contribution of the parameters described in the first sentence to the RDA axes are really minor. Thus, to me also this reads as a clear overinterpretation of the data. And the following sentence (line 356) is pure speculation (especially as the last two references only refer to the benefit high primary productivity can have for the corals).
Indeed, the sentence was ambiguous. We changed it to (line 363): “High dead coral cover correlated to strong upward velocities in February (section 3.2) and high coral rubble cover to strong upward velocities in August (Fig. 4b), whereas high live coral and sediment cover correlated to weak upward velocities in February (section 3.2).”
The fourth RDA axis contributes significantly, albeit explaining a smaller fraction of the variance than the other axes. We moved the sentence of the upwelling of nutrients to the new discussion section 4.3 and added references on the nutrient upwelling (line 391): “Nutrient-rich deep water is pushed up towards the lower mound flanks (Cyr et al., 2016) and the community associated to the dead coral framework releases substantial amounts of nutrients (Maier et al., 2021). These nutrients are likely transported by the predominantly upward water motions towards the ocean surface (Soetaert et al., 2016; Findlay et al., 2014; de Froe et al., 2022) where they stimulate primary productivity that, in turn, benefits the coral reefs (Eisele et al., 2011; Wienberg et al., 2022, 2020).”
- Mound engineering effect: a key finding here is that the mound engineering effect results in a deceleration on mound summits (line 256), which is picked up again in the discussion (line 332). So, why do the corals prefer the mound summits? Available current meter data usually show highest current speeds at or close to the summit. So, what is the meaning of your assessment that the mound engineering leads to reduced current speeds at the summits? And in Fig. 7a it is shown that bottom current speed significantly increases with mound height. But this does not refer to the summits? From line 333 this continues with a description of a hydrodynamic regime on coral mounds also seen in other studies. However, neither in Genin et al. nor in Soetaert et al., this pattern becomes obvious. See, for instance, points 1) and 2), for which neither of the two other paper present any evidence, and the same applies partly to the other points. These are to see to some extent in the modelling study by van der Kaaden et al. (2021), however, field observations do not show such a pattern (e.g., Genin et al. 1986, Dorschel et al., 2007). So, the pattern described here might be introduced as a new conceptual idea, but then it needs more discussion about the involved processes, and not simply the hint to some other papers (which actually do not show this pattern).
Horizontal bottom currents are only decelerated when mound summits are flat. We write (line 257): “…a deceleration on some mound summits”. In the discussion we refer to it again (line 351): “…, 3) flow deceleration at the summit if the summit is flat,…”.
Figure 7a does not show that horizontal bottom currents increase with mound height, but that the absolute effect on bottom current speed increases with mound height. To clarify we added to the caption of figure 7: “Note that the y-axes indicate the effect of mound presence on the hydrodynamic variables, not the actual values.”
The description of the hydrodynamic regime is based on Genin et al. (1986) where it is eloquently explained: “The distribution of corals on wide peaks can be explained by prediction of physical theory. When water is upwelled above wide topography, vortex lines are compressed and anticyclonic motion is induced, because of conservation of potential vorticity. The resulting flow-field is a combination of the velocity induced by this anticyclonic component and the free-stream velocity. Thus, looking downstream, flow acceleration occurs on the left side of the topographic feature (in the Northern Hemisphere) and deceleration above the centre and right side of the peak. […] The overall mean velocity pattern would therefore comprise alternating periods of acceleration on the peak’s edges and recurrent deceleration trend over the centre. This postulated velocity pattern is consistent with the lower densities of corals observed at the centre of the wide peaks of Jasper Seamount. Furthermore, the combination of tidal fluctuations and a general southward flow in the vicinity of the seamount could result in a stronger average acceleration on the eastern side of the seamount.”
The upwelling along the flanks and downwelling at the summit can be seen in Soetaert et al. (2016) figure 3B where the downwelling is also explicitly linked to the organic matter flux (their figure 3F). We substantiated the references to the hydrodynamics and rephrased (line 345): “These findings are in line with hydrodynamic theory of waterflow passed an object (e.g., Juva et al., 2020; Dewey et al., 2005; Baines, 1995; Genin et al., 1986) and with other studies reporting the hydrodynamics around coral mounds (e.g., Genin et al., 1986; de Froe et al., 2022; Findlay et al., 2013; Davies et al., 2009; Juva et al., 2020).”
- One comment on the methodology: what does not become clear from the text, is how the hydrodynamic model data with grid cells of 250 m in the horizontal are combined with the video frame data. According to Fig. 2, the individual transects are between ~500 and ~1500 m long; only transect 3 is much longer with ~3100 m, but has a less dense video frame density. Thus, most transects correspond to 2 to 6 grid cells in the hydrodynamic model, with seeing the small-scale typical summit and slope settings probably being partly very difficult to discern within 250 m grid cells. On the contrary, the video frame density usually is much, much higher, with a total of ~1200 frames. Thus, ~1200 frames, which in parts show quite some variability along a transect, are linked to ~30 grid cells, that seemingly cannot always discern between summits, slopes, etc. – how this is done becomes – at least for me – not clear from chapter 2.3. Chapter 2.2.2 explains that for each frame the nearest grid cell has been chosen for correlation, so, there might be a few more than the 30 grid cells mentioned above. But nevertheless, the question remains …
Indeed, the hydrodynamic simulations have a rather coarse resolution which may cause us to overestimate the importance of the broad-scale non-engineered hydrodynamics in the RDA. We did not stress this enough in our manuscript, so to section 4.1 we added (line 332): “This finding suggests that broad-scale environmental processes exert a stronger control on cold-water coral cover than engineered environmental variables. However, the rather coarse resolution of the hydrodynamic simulations might cause us to overestimate the importance of broad-scale environmental variables. Variance partitioning showed that 8.8% of the variation in benthic cover was explained solely by within-transect MEMs, compared to only 2.4% by between-transect MEMs. This signals that the variation at the within-transect scale might be related to environmental variables at a finer resolution. Such fine-scale variation in hydrodynamic variables is more likely to be caused by the interaction between water flow and fine-scale topography, i.e., ecosystem engineering. Since MEMs 3 to 5, that describe the mound foot-to-summit spatial scale, play a role in explaining the variation in benthic cover, cold-water coral mound ecosystem engineering might cause self-organized variation in benthic cover. “
- Obviously, the authors like their own phrases very much. For instance, the hint to better conditions for the corals “higher up in the water column” occurs at least 4 times, e.g., with the statement that the corals would benefit from being higher in the water column as there, POM quality and quantity is higher (only giving a reference for the POM distribution, but not for the link to the corals). This is followed up at the beginning of the discussion where coral performance again is linked to water depth – described as a result of this study. However, these statements are in contrast to one of the main conclusions referring to the “Massenerhebung effect” stating that not the position in the water column, but the position on the mound is a key parameter controlling coral growth. In addition, also the title of the ms is more a nice wording than a good description of the content of the paper. While the word “mountain” in the first part is only picked up in the paper with reference to terrestrial mountains, the second part of the title goes far beyond its actual content.
Depth explained 11.3% of the variation in benthic cover and was an important variable in the RDA axis 1 and 2, so that depth is an important variable explaining benthic cover is a finding of our study. As cold-water coral reef cover can be more or less dense, this result is not in contrast with the massenerhebung effect. To the discussion (4.2) we added (line 384): “Since depth is also an important variable explaining benthic cover, it is likely that the position of cold-water coral reefs on coral mounds is determined by the feedbacks between the mound and the hydrodynamics, but that shallower mounds have a denser reef cover than deeper mounds.”
- Table 1 caption: which colours?
We removed the reference to colours from the caption.
- Line 177: with (!) coral mounds
Yes, thank you.
Citation: https://doi.org/10.5194/egusphere-2023-949-AC2 - Further points to consider are listed below:
-
AC2: 'Reply on RC1', Anna van der Kaaden, 16 Aug 2023
-
RC2: 'Comment on egusphere-2023-949', Anonymous Referee #2, 17 Jul 2023
Review ms from Van der Kaaden et al, Building your own mountain: the effects, limits, and drawbacks of cold-water coral ecosystem engineering
The work of Van der Kaaden et al, in Biogeosciences egusphere-2023-949 provides a study of the Logachev cold-water coral mound province to investigate their ecosystem engineering ability at a local scale. To do so, they use 1) hydrodynamics modelling to test the effect of presence/absence of the coral mounds in currents and energy dissipation, and 2) statistical analysis (redundancy analysis) based on video transect images and environmental variables (both hydrodynamics and spatial variables from Moran’s Eigenvector Maps).
Major comments:
Self-structuration of the CWC reefs and their environmental drivers is an important question to assess especially in a global change context where sustainable management and future remediation strategies are to be defined, and the ms from Van der Kaaden et al, has the advantage to try to tackle this issue. The ms is well written and properly referenced. However, I feel like the ms in its current form needs further work before it can be published in Biogeosciences.
One of the difficulty I have with the ms is that some of the conclusions seem speculative, and the authors need to better assess their conclusions, their reliability to be more assertive.
For starter, their RDA is based on variables that fail to explain 55% of the variance, leaving a significant part of the variance in the CWC benthic cover unexplained: among the principal contributors are depth, and hydrodynamic parameters (which had already been demonstrated see Wagner et al, 2011, De Clippele et al, 2018, Mohn et al, 2014, Cyr 2016, Van der Kaaden et al, 2012 – all referenced in the ms), and the first four MEMS.
I find it hard to identify what the MEMs can realistically represent. Especially since mainly the larger ones are significant contributors to the RDA 2 main axis (MEM 1 to 4) which therefore corresponds to what? Several km scale? If MEM1 is 25km entire study area, .MEM4 is down to which spatial scale? Figure Appendix A1 lacks a color and a bubble scale to help interpret more easily these spatial variables. Especially since the conclusions of the authors are made on the mound scales – mound being in the order of hundreds of metres to a few kilometres. Maybe provide scalograms to indicate the portion of variance (R2) explained by each spatial scale (See Dray et al, 2012 https://doi.org/10.1890/11-1183.1). Regarding this spatial scale consideration, two remarks:
- Hydrodynamic simulations: the grid size being 250m, how confident are the authors to mimic and interpret the effect of mound removal when mound sizes are “in the order of hundreds meters to a few kilometers (Wheeler et al, 2006 https://link.springer.com/article/10.1007/s00531-006-0130-6)” mainly “up to 100 m high, 1–2 km long “(Van Weering et al., 2003)? Is 2 to 4 (max 8) grid points enough to draw hard conclusions? Did the authors perform any sensitivity analysis (maybe testing the minimum mound size to be have an effect on the hydrodynamics variables they consider)? Just to give the reader the sense on how the hydrodynamics effect of the CWC presence act.
- If the four first MEMS are the main spatial non-hydrodynamic environmental drivers of the CWC benthic cover, and therefore represent spatial variations “between transect” and “within transect MEM3 and 4” of several hundred meters, how can the authors draw conclusions at the mound scale? Although I am rather seduced by the masserhebung explanation and the analogy to mountain specie zonation, I fail to reconcile the scales at stake there.
Detailed comments:
Line 179 add positive in “i.e.a feedback”. As feedbacks can be both ways…
Figure 3: Images #81 and #52 are blurry. Maybe the authors could select other more focused images to make their point alhtough I do understand the need to select “end-members” pictures in terms of reef cover and rubble abundance.
Figure 8b) Side view – it is rather difficult to picture a convection loop with descending currents and ascending currents concomitant at that scale of a mound… Do the authors have any measurements (or model output) of such co-existing opposite vertical currents?
Line 320-330. Only transects 2, 3, 6 and 7 reach 600m, transects 1, 4 and 5 are deeper: the deep water mixing extending down to 600m won’t reach these mounds. Please suggest additional mechanisms for food supply to those.
Line 354 – 355. The authors state “High dead coral cover correlated to strong upward velocities in February and high coral rubble cover to strong upward velocities in August, whereas high live coral and sediment cover correlated to weak upward velocities in February (Fig. 4).”. Although I do see correlation for high coral rubble cover to strong upward velocities in August, I failed to see the two others. Both high dead coral and live coral cover seem correlated to strong downward February (Wdown February being on the same side of the axis that Dead and Live coral). Correlation between sediment and weak upward velocities in Feb also seems ambiguous.
Line 360-363. It is worth mentioning that hydrodynamics (upwelling and downwelling) and tidal cycles impact also the biogeochemical properties of the seawater (e.g. inorganic carbon and pH) above the reef and therefore the ecosystem functioning of the reef (Findlay et al, 2013 https://doi.org/10.1111/gcb.12256 Jiang et al, 2020 https://doi.org/10.1016/j.ocemod.2019.101555 , Findlay et al, 2014 https://doi.org/10.1038/srep03671). Maybe the authors should elaborate on that and integrate it to the discussion.
Citation: https://doi.org/10.5194/egusphere-2023-949-RC2 -
AC1: 'Reply on RC2', Anna van der Kaaden, 16 Aug 2023
Dear reviewer 2,
Thank you for your helpful comments. The results of RDA are only reliable when the explanatory variables are not too much correlated to one another. Therefore, we used the forward selection procedure to identify the most important explanatory variables contributing to explaining the variation in benthic cover. Adding more variables would increase the percentage of variance explained but would make the analysis unreliable. Compared to similar deep-sea studies using an RDA on benthic cover, our analysis actually performed better or comparable. Other factors that might influence benthic cover are anthropogenic influences, small-scale hydrodynamics, interspecific interactions, or random effects, but these could unfortunately not be included in this study.
That hydrodynamics are important for cold-water corals has been shown before, but we extend this knowledge by showing whether it is the ‘engineered’ hydrodynamics or ‘non-engineered’ hydrodynamics that exert the major influence on cold-water corals. In contrast to earlier studies (e.g., Wagner et al. 2011, De Clippele et al. 2017, Cyr et al. 2016, and van der Kaaden et al. 2021), we identified significant statistical correlations between cold-water coral cover and hydrodynamic variables.
To the discussion (line 293) we now added: “By explaining 45.1% of the variation in benthic cover, our set of variables performed well in comparison to similar deep-sea studies using RDA: e.g., Vad et al. (2020) explained 15.1% of the variation in benthic cover and Kazanidis et al. (2021) explained 26.7% of macrofauna variation. Henry et al. (2013) explained 65% of the variation in reef fauna using a combination of MEMs, hydrography, and bathymetric variables. So, the inclusion of both spatial factors (MEMs) and hydrodynamic variables appears to benefit the analysis. Especially within-transect MEMs described a relatively large part of the variation in our benthic cover data that was not explained by the environmental variables. Henry et al. (2013) similarly found that community composition was structured mainly by broad-scale environmental forcings, but that fine-scale MEMs described a significant fraction of the variation in community composition. This within-transect variation in benthic cover can be related to spatial patterns of fauna (Henry et al., 2013), such as self-organized regular spatial patterns in coral reefs (van der Kaaden et al., in press) and to the typical reef zonation on cold-water coral mounds (section 4.2).”
Below, I will reply in detail to your comments, displaying your comments in italics followed by our reply.
- I find it hard to identify what the MEMs can realistically represent. Especially since mainly the larger ones are significant contributors to the RDA 2 main axis (MEM 1 to 4) which therefore corresponds to what? Several km scale? If MEM1 is 25km entire study area, .MEM4 is down to which spatial scale? Figure Appendix A1 lacks a color and a bubble scale to help interpret more easily these spatial variables. Especially since the conclusions of the authors are made on the mound scales – mound being in the order of hundreds of metres to a few kilometres. Maybe provide scalograms to indicate the portion of variance (R2) explained by each spatial scale (See Dray et al, 2012 https://doi.org/10.1890/11-1183.1). Regarding this spatial scale consideration, two remarks:
Hydrodynamic simulations: the grid size being 250m, how confident are the authors to mimic and interpret the effect of mound removal when mound sizes are “in the order of hundreds meters to a few kilometers (Wheeler et al, 2006 https://link.springer.com/article/10.1007/s00531-006-0130-6)” mainly “up to 100 m high, 1–2 km long “(Van Weering et al., 2003)? Is 2 to 4 (max 8) grid points enough to draw hard conclusions? Did the authors perform any sensitivity analysis (maybe testing the minimum mound size to be have an effect on the hydrodynamics variables they consider)? Just to give the reader the sense on how the hydrodynamics effect of the CWC presence act.
In Van der Kaaden et al. (2021), we investigated the effect of cold-water coral mounds of different sizes on the hydrodynamics and showed that mounds influence the hydrodynamics already from about 25m height. By correlating the height of all Logachev mounds to their effect on the hydrodynamics we here show that all hydrodynamic variables are significantly affected by cold-water coral mound presence, except for downward velocities in February. For this correlation, the scale of the hydrodynamic simulations is of lesser importance as we calculated the mean values for each mound.
In the RDA, the scale of the hydrodynamic simulations is more important and might lead to an overestimation of the importance of broad-scale environmental processes over finer-scale processes. Indeed, we did not stress this enough in our manuscript (also mentioned by reviewer 1), so to section 4.1 we added (line 332): “This finding suggests that broad-scale environmental processes exert a stronger control on cold-water coral cover than engineered environmental variables. However, the rather coarse resolution of the hydrodynamic simulations might cause us to overestimate the importance of broad-scale environmental variables. Variance partitioning showed that 8.8% of the variation in benthic cover was explained solely by within-transect MEMs, compared to only 2.4% by between-transect MEMs. This signals that the variation at the within-transect scale might be related to environmental variables at a finer resolution. Such fine-scale variation in hydrodynamic variables is more likely to be caused by the interaction between water flow and fine-scale topography, i.e., ecosystem engineering. Since MEMs 3 to 5, that describe the mound foot-to-summit spatial scale, play a role in explaining the variation in benthic cover, cold-water coral mound ecosystem engineering might cause self-organized variation in benthic cover.”
- If the four first MEMS are the main spatial non-hydrodynamic environmental drivers of the CWC benthic cover, and therefore represent spatial variations “between transect” and “within transect MEM3 and 4” of several hundred meters, how can the authors draw conclusions at the mound scale? Although I am rather seduced by the masserhebung explanation and the analogy to mountain specie zonation, I fail to reconcile the scales at stake there.
MEM 3, MEM 4, and MEM 5 describe the variation at the scale from mound foot to mound summit. MEM 3 is the most important explanatory variable explaining the variation in rubble cover versus non-rubble cover (RDA axis 2) and MEM 4 plays a role in explaining the variation in coral reef versus non-reef cover (RDA axis 1). Between-transect MEMs 1 and 2 describe more of the variation in the RDA, but variation partitioning shows that a larger share of the variance in benthic cover is explained solely by within-transect MEMs (i.e., MEM 3 and 4 and to a lesser extent 5, 6, and 7). This indicates that the mound foot-to-summit variation in benthic cover is, for quite a significant part (i.e., almost as much as depth), not explained by environmental variables included in this study. Spatial variation that cannot be related to environmental variables can be self-induced, as with a massenerhebung effect. It is not unreasonable to think that a massenerhebung effect might happen at cold-water coral mounds, which is why we propose that a massenerhebung effect exists for cold-water coral mounds.
To clarify the scale described by the different MEMs, we changed section 2.2.3 (line 146): “At the broadest spatial scale (signalling differences between transects) datapoints were classified according to their east-to-west position (MEM 1) or according to their north-to-south position (MEM 2). At finer scales (signalling differences within transects) datapoints were classified by MEMs according to their position along a mound flank (MEM 3, 4, and 5), or according to regular intervals along transects (MEM 6 to MEM 15).”
And in line 360: “…we found that high cold-water coral reef and rubble cover correlated to the spatial scale from mound foot-to-summit…” instead of “…the spatial scale of a mound…”.
We will clarify the supplementary figure illustrating the MEM values: showing the values of MEM 1 to 7, along three example transects for each annotated video frame.
We further worded our conclusions more carefully (line 382): “We here hypothesize that this massenerhebung effect, …, applies to cold-water coral mounds too.” Instead of: “We here present evidence…”
- Detailed comments:
- Line 179 add positive in “i.e.a feedback”. As feedbacks can be both ways…
Yes, that is true, but we intend to refer to a feedback (writing that the effect of the mound on the hydrodynamics increases with increasing mound height). At some sides of the mound bottom currents are decelerated, so we refer to both a positive and negative feedback.
- Figure 3: Images #81 and #52 are blurry. Maybe the authors could select other more focused images to make their point alhtough I do understand the need to select “end-members” pictures in terms of reef cover and rubble abundance.
We used higher-resolution images from a different transect for figure 3.
- Figure 8b) Side view – it is rather difficult to picture a convection loop with descending currents and ascending currents concomitant at that scale of a mound… Do the authors have any measurements (or model output) of such co-existing opposite vertical currents?
Indeed, such ‘upwelling’ and ‘downwelling’ can be seen in measurements and model output (e.g., Soetaert et al. 2016, Davies et al. 2009, Cyr et al. 2016), and is known from theory of waterflow around an obstacle (e.g., Juva et al. 2020). Of course, the up- and downward motion do not necessarily exist at the same time, but can be the result of tidal processes. We clarified in the caption of figure 8: “The up- and downward water motions do not necessarily happen simultaneously but can result from tidal motions.”
- Line 320-330. Only transects 2, 3, 6 and 7 reach 600m, transects 1, 4 and 5 are deeper: the deep water mixing extending down to 600m won’t reach these mounds. Please suggest additional mechanisms for food supply to those.
Deep-water mixing brings organic matter from the ocean surface some hundred meters down into the water column. This would increase the organic matter quality and quantity at all depths not only within the mixed layer. So, the mechanisms that we propose are relevant also for the deeper mounds.
- Line 354 – 355. The authors state “High dead coral cover correlated to strong upward velocities in February and high coral rubble cover to strong upward velocities in August, whereas high live coral and sediment cover correlated to weak upward velocities in February (Fig. 4).”. Although I do see correlation for high coral rubble cover to strong upward velocities in August, I failed to see the two others. Both high dead coral and live coral cover seem correlated to strong downward February (Wdown February being on the same side of the axis that Dead and Live coral). Correlation between sediment and weak upward velocities in Feb also seems ambiguous.
Indeed, this was unclear. The reference to figure 4 concerned only the coral rubble distinction from RDA axis 2. We changed the sentence (line 363): “High dead coral cover correlated to strong upward velocities in February (section 3.2) and high coral rubble cover to strong upward velocities in August (Fig. 4b), whereas high live coral and sediment cover correlated to weak upward velocities in February (section 3.2).”
- Line 360-363. It is worth mentioning that hydrodynamics (upwelling and downwelling) and tidal cycles impact also the biogeochemical properties of the seawater (e.g. inorganic carbon and pH) above the reef and therefore the ecosystem functioning of the reef (Findlay et al, 2013 https://doi.org/10.1111/gcb.12256 Jiang et al, 2020 https://doi.org/10.1016/j.ocemod.2019.101555 , Findlay et al, 2014 https://doi.org/10.1038/srep03671). Maybe the authors should elaborate on that and integrate it to the discussion.
Yes, we like this suggestion very much, so we added a section (4.3) on “Ecosystem engineering effects on cold-water corals and climate change”, where we discuss the mechanisms by which the hydrodynamic zones might affect cold-water coral growth, including temperature, oxygen, and pH.
Citation: https://doi.org/10.5194/egusphere-2023-949-AC1 - I find it hard to identify what the MEMs can realistically represent. Especially since mainly the larger ones are significant contributors to the RDA 2 main axis (MEM 1 to 4) which therefore corresponds to what? Several km scale? If MEM1 is 25km entire study area, .MEM4 is down to which spatial scale? Figure Appendix A1 lacks a color and a bubble scale to help interpret more easily these spatial variables. Especially since the conclusions of the authors are made on the mound scales – mound being in the order of hundreds of metres to a few kilometres. Maybe provide scalograms to indicate the portion of variance (R2) explained by each spatial scale (See Dray et al, 2012 https://doi.org/10.1890/11-1183.1). Regarding this spatial scale consideration, two remarks: