Biomass Mapping for an Improved Understanding of the Contribution of Cold-Water Coral Carbonate Mounds to C and N Cycling

This study used a novel approach combining biological, environmental, and ecosystem function data of the Logachev cold-water coral carbonate mound province to predictively map coral framework (bio)mass. A more accurate representation and quantification of cold-water coral reef ecosystem functions such as Carbon and Nitrogen stock and turnover were given by accounting for the spatial heterogeneity. Our results indicate that 45% is covered by dead and only 3% by live coral framework. The remaining 51%, is covered by fine sediments. It is estimated that 75,034–93,534 tons (T) of live coral framework is present in the area, of which ∼10% (7,747–9,316 T) consists of Cinorg and ∼1% (411–1,061 T) of Corg. A much larger amount of 3,485,828–4,357,435 T (60:1 dead:live ratio) dead coral framework contained ∼11% (418,299–522,892 T) Cinorg and <1% (0–16 T) Corg. The nutrient turnover by dead coral framework is the largest, contributing 45–51% (2,596–3,626 T) C year–1 and 30–62% (290–1,989 T) N year–1 to the total turnover in the area. Live coral framework turns over 1,656–2,828 T C year–1 and 53–286 T N year–1. Sediments contribute between 1,216–1,512 T C year–1 and 629–919 T N year–1 to the area’s benthic organic matter mineralization. However, this amount is likely higher as sediments baffled by coral framework might play a much more critical role in reefs CN cycling than previously assumed. Our calculations showed that the area overturns 1–3.4 times the C compared to a soft-sediment area at a similar depth. With only 5–9% of the primary productivity reaching the corals via natural deposition, this study indicated that the supply of food largely depends on local hydrodynamical food supply mechanisms and the reefs ability to retain and recycle nutrients. Climate-induced changes in primary production, local hydrodynamical food supply and the dissolution of particle-baffling coral framework could have severe implications for the survival and functioning of cold-water coral reefs.

This study used a novel approach combining biological, environmental, and ecosystem function data of the Logachev cold-water coral carbonate mound province to predictively map coral framework (bio)mass. A more accurate representation and quantification of cold-water coral reef ecosystem functions such as Carbon and Nitrogen stock and turnover were given by accounting for the spatial heterogeneity. Our results indicate that 45% is covered by dead and only 3% by live coral framework. The remaining 51%, is covered by fine sediments. It is estimated that 75,034-93,534 tons (T) of live coral framework is present in the area, of which ∼10% (7,747-9,316 T) consists of C inorg and ∼1% (411-1,061 T) of C org . A much larger amount of 3,485,828-4,357,435 T (60:1 dead:live ratio) dead coral framework contained ∼11% (418,299-522,892 T) C inorg and <1% (0-16 T) C org . The nutrient turnover by dead coral framework is the largest, contributing 45-51% (2,596-3,626 T) C year −1 and 30-62% (290-1,989 T) N year −1 to the total turnover in the area. Live coral framework turns over 1,656-2,828 T C year −1 and 53-286 T N year −1 . Sediments contribute between 1,216-1,512 T C year −1 and 629-919 T N year −1 to the area's benthic organic matter mineralization. However, this amount is likely higher as sediments baffled by coral framework might play a much more critical role in reefs CN cycling than previously assumed. Our calculations showed that the area overturns 1-3.4 times the C compared to a soft-sediment area at a similar depth. With only 5-9% of the primary productivity reaching the corals via natural deposition, this study indicated that the supply of food largely depends on local hydrodynamical food supply mechanisms and the reefs ability to retain and recycle nutrients. Climate-induced changes in primary production, local hydrodynamical food supply and the dissolution of particle-baffling coral framework could have severe implications for the survival and functioning of cold-water coral reefs.

INTRODUCTION
Cold-water coral (CWC) carbonate mounds are important marine ecosystems . They are topographic seafloor structures that can be several hundreds of meters in height and have accumulated through successive periods of reef development, sedimentation and (bio)erosion over glacialinterglacial periods (Kenyon et al., 2003;Van Weering et al., 2003;Mienis et al., 2007;Roberts et al., 2009). They are hotspots of biomass and biodiversity and provide essential ecosystem functions through nutrient [Carbon (C) and Nitrogen (N)] cycling in a resource-limited deep sea (Henry and Roberts, 2007;van Oevelen et al., 2009;Armstrong et al., 2012). However, significant gaps remain in our understanding of the spatial distribution of their overall biomass and capacity to remineralise organic matter (OM) (De Clippele et al., 2021).
Cold-water coral reefs depend on OM produced at the ocean's surface to support their growth (Duineveld et al., 2004Kiriakoulakis et al., 2005). This OM can be transported to the reef from surface waters through deposition, tidal downwelling, nepheloid layers and deep-water advection Davies et al., 2009;Findlay et al., 2013;Mohn et al., 2014;Soetaert et al., 2016). In addition, when reefs (tens of meters high) accumulate over time to form large CWC carbonate mounds (hundreds of meters high), they can induce a "topographicallyenhanced carbon pump" (Soetaert et al., 2016). The mounds large size interrupts the currents, which creates downwelling events bringing OM from surface waters to the mound's summits and upper flanks (Guinan et al., 2009;Mohn et al., 2014;Rengstorf et al., 2014;Soetaert et al., 2016). Baffling of currents caused by the coral framework can also locally increase the POM concentration at the reefs (Soetaert et al., 2016).
The availability of this food is a major determinant controlling CWCs occurrence and the zonation of macrohabitats on the mounds (De Clippele et al., 2019;Maier et al., 2021). The mound bases are covered by sediments (bio-and siliciclastic sands), pebbles, cobbles and boulders (de Haas et al., 2009). Dense Lophelia pertusa patches characterize the summits of the carbonate mounds, while the flanks of the mounds are covered with patches of coral rubble, dead coral branches and living corals (Kenyon et al., 2003;Van Weering et al., 2003;de Haas et al., 2009;De Clippele et al., 2019;Maier et al., 2021). Dead coral framework is particularly biodiverse as it provides complex micro-and macrohabitats for diverse communities (Jonsson et al., 2004;Henry and Roberts, 2007). It is this living fauna (including e.g., anthozoans, hydroids, ophiuroids, and sponges) that contributes the most to a reef 's capacity to mineralize OM (de Froe et al., 2019;Maier et al., 2019Maier et al., , 2020De Clippele et al., 2021).
Knowing how much live and dead coral framework biomass is present on a CWC reef and their contribution toward OM mineralization is critical information to understand how well the reef is functioning. It also provides a baseline that can help us understand the extent of the potential effects of ocean acidification, warming and decreases in ocean O 2 levels on these vulnerable ecosystems (Hennige et al., 2014(Hennige et al., , 2015(Hennige et al., , 2020Roberts and Cairns, 2014;Sweetman et al., 2017). To estimate biomass and OM mineralization on CWC carbonate mounds, we apply the novel approach by De Clippele et al. (2021). This approach uses surface area measurements of the coral L. pertusa, extracted from high-definition (HD) video frames and combines this with biomass and respiration data. We hypothesize that this method allows to map live and dead coral framework at the CWC Logachev Mound province (LMP) and quantify the ecosystem function of this area.

METHODOLOGY Location
The LMP consists of a cluster of CWC carbonate mounds located on the south-eastern slope of Rockall Bank in the North-East Atlantic (Kenyon et al., 2003; Figure 1). The CWC carbonate mounds are between 5 and 360 m tall, up to a few kilometers long and located between 500 and 1,000 m depth (Kenyon et al., 2003;de Haas et al., 2009). The dominant current direction in the LMP is in a southwest direction, following from a clockwise circumventing flow around Rockall bank , while the local diurnal barotropic tide causes cross slope transport in a northwest-southeast direction White, 2007).

Biological Data
Eight HD video transects were recorded during the Changing Oceans 2012 expedition, RRS James Cook cruise 073 (Roberts, 2013), using the Remotely Operated Vehicle (ROV) Holland-1 (more details in De Clippele et al., 2019; Table 1). Using the software Photoshop CC 2018, video frames were extracted every 500th frame. The video frames were used to measure the surface area of live and dead coral framework (see Section "Biomass Estimation"). The remaining area (total area minus [dead + live] coral framework) was referred to as sediment. However, hard substrates such as pebbles, cobbles, boulders and lithified substrate can also be present (De Clippele et al., 2019;Maier et al., 2021). The ROV was equipped with two parallel pointers, marking a fixed distance of 10 cm on the video frames, which was used to scale the images.
In addition, data on the dry weight, C inorg and C org stock of the live and dead coral framework, collected with a NIOZ boxcorer (diameter: 50 cm; height 50 cm; surface area ∼0.2 m 2 ) was used. Six cores were collected during the 2017 R/V Pelagia research cruise and used to derive ex situ benthic O 2 and N flux measurements of the CWC community (see Table 1 in de Froe et al., 2019). Photographs were taken of the core surface after sampling and used to calculate the surface area (m 2 ): dry weight (kg) ratio. The photographs were scaled using the dimensions of the boxcorer.

Environmental Data
Particulate organic matter (POM) concentrations were obtained from a POM model with a resolution of 250 m × 250 m (Soetaert et al., 2016). This model provides values that represent the concentration of reactive freshly-produced organic matter available in the water column. These are below the actual measured values of POM concentration, which additionally  include refractory organic matter (Soetaert et al., 2016). In addition, terrain variables were extracted from bathymetry data provided by the Irish National Seabed Survey program (INSS) at a 20 m × 20 m resolution. The following topographic terrain variables were derived from the bathymetry data using the ArcGIS 10.1, ESRI Software and the Benthic Terrain Modeller (Wright et al., 2005): depth, slope, aspect (eastness and northness), rugosity (calculated at two spatial scales, using a square kernel window of, respectively, 3 pixels × 3 pixels and 9 pixels × 9 pixels) and bathymetric positioning index (BPI; calculated at two spatial scales using an annulus kernel window with inner and outer radius of, respectively, 3 × 6 and 6 × 9 cells). More information on these variables is provided in De Clippele et al., 2019.

Coral Presence Habitat
Our model area was defined by the habitat suitability model of CWC presence/absence produced by Rengstorf et al. (2014) and covers 253 km 2 . This habitat suitability model was chosen as particulate organic matter (POM) is used as an environmental variable to explain the spatial variability in coral biomass. The POM was calculated by Soetaert et al. (2016) who used the above mentioned habitat suitability model to study benthic respiration and the amount of food supplied to the LMP . Because the POM model assumed that OM deposition/uptake was increased by a constant factor in the presence of corals, we cannot compare coral-presence habitat with coral-absence habitat (Soetaert et al., 2016).

Biomass Estimation
Biomass is here defined as the live tissue of a specimen.
In this study we therefore refer to "(bio)mass" to indicate the differentiation between measuring mass and biomass for, respectively, live and dead coral framework. The approach by De Clippele et al. (2021), was adapted due to a difference in coral morphologies, i.e., the presence of coral thickets at the LMP rather than the globular colonies at the Mingulay Reef (Figure 2).
Here, to convert surface area to (bio)mass, (bio)mass data from boxcores collected at the LMP were used (de Froe et al., 2019). The steps are described in detail below and in Figure 3.

Image Analyses
Step 1 The video still frames from the HD videos (see Section "Biological Data") ( Figure 3A), were imported in Adobe Photoshop. Bad quality images or images that overlapped were excluded from analyses. In Photoshop, the laser-scale dots, live and dead coral framework were labeled each with a unique color aided by Photoshop's "quick selection tool" (Figure 3B) (van der Kaaden and De Clippele, 2021). The benthic surface area covered by dead and live coral framework was calculated in R. An R-function (van der Kaaden and De Clippele, 2021) was used to semiautomatically calculate the percentage cover and image size from the labeled images, from which the surface area in m 2 could be calculated. This faster method is an alternative to the method proposed in De Clippele et al. (2021) where the opensource software ImageJ2 (Greene et al., 1999;Rueden et al., 2017) was used.

Predictive Mapping
Step 2 The surface area data points were imported in ArcGIS and combined in 20 m sub-samples (x-axis). The length of 20 m was chosen as this length gave the most accurate representation of the coral framework variability in relation to the multibeam grid cell size. Then, the ArcGIS Extract Values to Points tool was used to extract the environmental variables (i.e., depth, BPI, slope, rugosity, eastness, northness, and POM) ( Figure 3C) associated with each sub-sample data point.
Step 3 The response (i.e., surface area) and explanatory (i.e., environmental) variables were then used to model a predictive map using the Random Forest approach ( Figure 3D) with the randomForest package in R (Breiman, 2001). This supervised classification methodology is referred to as a regression tree with a number of simple decision trees. Each tree is based on a bootstrapped sample of the response and explanatory data set. This group of simple trees vote for the most popular class, which is capable of predicting a response when presented with a set of explanatory variables (Cutler et al., 2007;Rogan et al., 2008). Random Forest modeling is commonly used to produce predictive maps (Baccini et al., 2008;Wei et al., 2010;Zhang, 2015;Conti et al., 2019;De Clippele et al., 2021) and provides similar results to approaches using logarithmic regression and Deep Neural Networks (Conti et al., 2019). Here, the training dataset contained one-third of the total data points. This Random Forest model can then be applied to a new set of the same response variables to create a predictive map of the unseen data of the whole area. Correlated environmental variables (<0.5) were removed prior to analyses, using the cor test in R. The importance of the environmental variables in predicting surface area was assessed by calculating the Mean Decrease in Accuracy for each variable to indicate their contribution to the model performance.
To evaluate the uncertainty of the model outputs, we first used a bootstrap technique to produce estimates of model uncertainty (Rowden et al., 2017). This is done by repeating the Random Forest model one hundred times, with the same model but with a replacement random sample of the training data each time. This results in 100 estimates from which the coefficient of variance (CV) was calculated to examine the models output stability. This provides a range of how the Random Forest model output varies and is measured as the standard deviation/mean × 100 (Wei et al., 2010).

Biomass Calculation
Step 4 From the predictive Lophelia reef maps (see step 3), the total amount of live and dead coral framework surface area for the Frontiers in Marine Science | www.frontiersin.org   (Table 3). From the known surface area of the boxcorer (see Section "Biological Data"), the dry weight of live and dead coral framework per boxcore sample was calculated (kg dry weight boxcore −1 ). The benthic surface area of live and dead coral framework from boxcore photographs (Supplementary Materials) was then used to calculate a conversion factor of live/dead coral framework surface area to live/dead coral framework dry mass for each boxcore sample. The mean conversion factor of 8.01 ± 5.52 kg m −2 and 33.78 ± 6.27 kg m −2 for live and dead coral framework, respectively, was then used to convert benthic surface area measurements from the HD video extracted frames to skeletal dry weight (kg dry weight) (Figure 3F), i.e., dry weight (kg) = video surface area (m 2 ) × conversion factor (kg m −2 ). From this, the biomass (live tissue mass) was calculated using the linear relationship between tissue dry weight and tissue and skeletal weight with the following equation from Hennige et al. (2014): TDW = 0.0415 (TWW + SDW) + 0.0849. With TDW = tissue dry weight (g) = SDW * 5%; TWW = tissue wet weight (g); SDW = skeletal dry weight (g) (De Clippele et al., 2021;Hennige et al., 2014). Sediment were calculated as the remainder (total area minus [dead + live] coral framework). The above conversion calculation was also used in the ArcGIS raster calculator tool to convert the surface area predictive map to a (bio)mass. The standard deviation of the conversion factors was used to calculate the absolute minimum, mean and maximum biomass (see the "Results" section and Supplementary Materials).

Carbon and Nitrogen Turnover and Stock
Step 5 The total biomass and sediment surface area data was used to calculate the yearly C and N turnover for the area, using O 2 consumption data reported in literature ( Table 2 and Figure 3G). Carbon and N turnover are here defined as the conversion of ingested food into biomass and loss by respiration as CO 2 and DIN. The C turnover is calculated from the total O 2 consumption assuming a respiratory quotient (C:O 2 ratio) of 1:1 (Glud, 2008) ( Figure 3H). This does not account for temporal changes in O 2 consumption that the coral might experience during its lifetime, as this data is currently not available. The live and dead coral framework were multiplied with their respective DIN release rates. More details on this calculation can be found in Step 3 in De Clippele et al. (2021). In Table 3 by de Froe et al. (2019), the percentage of C org and N org of live and dead coral framework per boxcore sample are reported. The C org and N org stock of live and dead coral framework per boxcore sample was obtained by multiplying the percentage C org and N org content of live and dead coral framework with their total dry weight (C or N kg −1 dry weight). These values were used to calculate a mean conversion factor for the mass calculated in step 3 to C org and N org ( Table 3). The C inorg stock was calculated by multiplying the CaCO 3 mass by 0.12 (Windholz et al., 1983;De Clippele et al., 2021). To account for uncertainties the absolute minimum, mean and maximum CN turnover and C stock were calculated (see the "Results" section and Supplementary Materials).

Predictive Maps
Our model predicts live coral framework covering 8 km 2 (3%) and dead covering 115 km 2 (45%) of the CWC habitat area. The remaining 130 km 2 (51%) is therefore considered to consist of sediment.
The environmental variables used in the mean live coral framework Random Forest model explained 65.54% of the variation in the data. The environmental variables that contributed most to explaining the spatial variability in the amount of live coral framework were BPI (inner cell radius 6 × outer cell radius 9), POM concentration and rugosity (9 × 9 cells) (Figure 4). The live coral framework biomass map (Figure 5) illustrated that the highest live coral biomass is located on the summits of the mounds. The study area contained a total live coral framework skeletal mass of 64,054 T C inorg (range: 62,280-77,635 T) and biomass of 13,117 T C org (range: 12,754-15,899 T). Our model results showed highest uncertainty at deeper depths and at the most eastern mounds (Figure 6).
The environmental variables used in the mean dead coral framework Random Forest model explained 54.21% of the variation in the data. The environmental variables that contributed most to explaining the spatial variability in the dead coral framework were BPI, depth and POM (Figure 4). The dead coral framework predictive map (Figure 7) showed that the highest mass is located on the northeast flanks and on the summits of the mounds, and that it decreases with depth. The area has a total mean dead framework skeletal mass of 2,875,706 T C inorg (range: 3,485,828-4,357,435 T) and variability was also here higher at depth and the most eastern mounds (Figure 8).

DISCUSSION
This study applied a new methodology to map live and dead coral framework biomass at the Logachev Mound Province. These biomass maps were used to estimate region-scale inorganic CN turnover, as well as the organic and inorganic CN standing stocks. Even though the reefs at the LMP occur in relatively deep and under food-limited waters compared to shallower inshore reefs (De Clippele et al., 2021) they contribute significantly to the global CN turnover and CN stock. This is as CWC mounds in the LMP, cover a large area and form big mounds due to their ability to persist over glacial-interglacial time scales. This work advances our growing knowledge of their significance to remineralise OM, a criteria used to define Ecologically or Biologically Significant Areas (EBSAs) (Titschack et al., 2015;Johnson et al., 2018).

Distribution of Dead and Live Coral Framework
Spatial differences in environmental conditions drive the small and large scale patterns in biomass observed at the LMP. This study showed that bathymetric positioning index (BPI) is the most important environmental predictor of both live and dead coral framework. This is as coral carbonate mounds form through periods of successive reef development . When reef growth dominates over erosion, many small reefs will cover the surface of the mound before they merge and continue the development cycle . These smaller reefs coincide with positive BPI values across the LMP (De Clippele et al., 2017). Our predictive biomass map also indicates that live coral framework is predominantly located on the summits of the mounds, which can largely be explained by the higher amount of available POM. In these relative food-limited waters, the supply of POM, i.e., their food, from surface water to the reefs is important to support the high metabolic C demand of the live coral framework (Davies et al., 2009;Roberts et al., 2009). Soetaert et al. (2016) suggest that the high elevation of the coral carbonate mounds induces downwelling and hence POM supply from the ocean surface, a concept described as topographically-enhanced carbon pump. The baffling created by the coral framework can locally increase the POM concentration measured on the reefs (Soetaert et al., 2016). These higher POM concentrations, in turn, can increase the biomass of the live coral framework. This positive feedback loop could explain why the prediction of live coral is strongly driven by POM concentration. The explanatory power of POM could increase even more if a higher resolution POM model would be used. Reduced POM supply linked with bio-and hydrodynamic erosion are the likely causes for the observed reduction in dead coral cover at greater depths . This is also shown by Maier et al. (2021), where more degraded coral framework is found at deeper depths and less degraded framework at more shallow depths. Live corals form complex three-dimensional frameworks, and increase the local terrain roughness (Jenness, 2002). Rugosity, a measure of  terrain roughness, therefore, represents the third most important variable explaining the prediction of live coral at the LMP. Our model predicted more dead than live coral framework in the Logachev Mound Province, which is supported by previous studies on CWC reefs (De Clippele et al., 2019de Froe et al., 2019;Maier et al., 2021). De Clippele et al. (2019) and Maier et al. (2021) calculated the percentage cover from ROV videos and found that dead coral framework covered 35-93% and live coral framework covered 3-25% of the transects at the LMP. This study found that dead coral framework surface area covered 45%, compared to live coral framework covering 3% of the area. Here, we report 60 times more dead framework mass than live coral mass in the area, which is twice as much as reported by earlier studies (27 times; de Froe et al., 2019). de Froe et al. (2019) based this difference in dry mass on collected boxcores, while here, the whole area is accounted for by means of predictive modeling. Other studies, such as Conti et al. (2019) have calculated the percentage of dead and live coral framework using video mosaic segmentation and classification approaches. They found that the Piddington Mound has 33-43% dead framework (rubble + dead coral framework), 2-3% live coral framework and 48-58% sediments (incl. dropstones) (Conti et al., 2019). While this is similar to what was found in our area, the Piddington mound has a spatial extent of approximately 40 m × 60 m and is one of the smaller mounds found in the Belgica Mound Province. The percentages of the substrates found will vary depending on the extent and the spatial heterogeneity of CWC reef area analyzed. The latter underlines the important contribution of representative sampling techniques and predictive models to more accurately represent CWC framework surface area coverage and biomass. Dead coral framework is important as it facilitates the high biodiversity typical of CWC reefs (Henry and Roberts, 2007) and contributes substantially more to reef fauna biomass and benthic fluxes (de Froe et al., 2019;De Clippele et al., 2021;Maier et al., 2021). Live corals protect themselves against colonization, for example by production of mucus (Freiwald, 2002;Wild et al., 2008;Buhl-Mortensen et al., 2010). In contrast, dead, unprotected coral framework is more easily colonized and provides the majority of micro-and macrohabitats in a CWC reef (Mortensen and Fosså, 2006;Buhl-Mortensen et al., 2010).

Oxygen Consumption and Nitrogen Release
Cold-water coral reefs are hotspots of O 2 consumption and N release, i.e., OM mineralization (van Oevelen et al., 2009;Cathalot et al., 2015;de Froe et al., 2019;De Clippele et al., 2021;Maier et al., 2021). The average C turnover at the LMP, which we derived from O 2 consumption measurements, was 1-3.4 times the global average for a soft-sediment area at the same depth (Glud, 2008). Dead coral framework contributed 49%, live coral framework 29%, and sediments 22% to the total C turnover of the area. At the same time, the reefs at the LMP released 1.9 times more DIN compared to adjacent soft-sediment grounds (de Froe et al., 2019). Here the majority of the DIN was released in the form of NO − 3 by both dead coral framework (64% of the total N turnover) and sediments (27% of the total N turnover), while NH + 4 was released by live coral framework (9% of the total N turnover). It should be noted that the partitioning of DIN release in NH 4 + and NO 3 − originates from the model assumption, i.e., that live cold-water corals typically release ammonium as metabolic end product (Khripounoff et al., 2014), while dead framework and sediment release mostly nitrate, due to the activity of nitrifying microorganisms (de Froe et al., 2019;Maier et al., 2021). Tidal induced upwelling of this nutrientrich reef water could promote new phytoplankton primary production in the surface waters, which in turn would increase OM export to the reef (Davies et al., 2009;Eisele et al., 2011;Hebbeln et al., 2014;Soetaert et al., 2016). Such a loop has been suggested for cold-water coral ecosystems at the shallower Porcupine Bank (White et al., 2005) and could explain how coldwater coral reefs are sustained in the relative resource-poor deep sea. If such a loop is also present at the deeper LMP remains to be determined.
However, it is important to note that the C turnover reported in this study (5.21-6.44 mmol C m −2 d −1 ) is 3-12 times lower than previously reported respiration measurements (van Oevelen et al., 2009;de Froe et al., 2019) (11-75 mmol C m −2 d −1 ). There are three reasons why we might observe this difference. Firstly, in contrast to previous studies, our study accounted for the spatial variability in the biomass of the coral framework across the whole region and revealed that on average 51% is covered by sediments, 45% by dead coral framework, and 3% by live coral framework. Respiration measurements collected by Aquatic Eddy-Correlation (AEC) or boxcore samples are not able to grasp the spatial heterogeneity of such a large area (de Froe et al., 2019;De Clippele et al., 2021). Our study therefore provides a more accurate representation of the total C and N turnover in the area, as it accounts for the spatial complexity.
Secondly, our calculations might be an underestimation as the physical structure of the coral framework baffles sediment (de Haas et al., 2009;de Froe et al., 2019). This type of sediment has a higher OM concentration and C turnover rate (5.32 ± 0.59 mmol O 2 m −2 d −1 , de Froe et al., 2019) compared to non-reef sediments adjacent to the reefs (de Froe et al., 2019;de Haas et al., 2009). Given that the area that is covered by live and dead coral framework, it will also contain baffled sediment which could significantly increase the total OM mineralization capacity of the area (de Haas et al., 2009). To illustrate the potential contribution of the baffled sediments to the LMP carbon turnover, we assume the area covered by live and dead coral framework is also covered by baffled sediments and add this to our calculation. In this scenario, the result suggests that communities associated with dead coral branches would contribute a mean of 33% (3,056 T C year −1 ), baffled sediments 33% (3,025 T C year −1 ), sediments 15% (1,389 T C year −1 ), and live coral framework 19% (1,793 year −1 ) to the total benthic C turnover at the LMP. If we recalculate the C turnover per square meter for this new scenario, a total mean of 8.34 mmol C m −2 d −1 is found instead of a mean of 5.64 mmol C m −2 d −1 . This indicates that cold-water coral carbonate mound sediments might play a much more important role in the C turnover of CWC reefs and carbonate mounds than previously thought.
Thirdly, our predictive maps indicate that in situ measurements by de Froe et al. (2019) were deployed in areas with relatively high coral framework cover (Figure 9). This might provide an overestimation that does not reflect the spatial variability of the C turnover in the area. The Oreo Mound, where we found a particularly high live coral framework cover (Figure 9A), showed the highest AEC O 2 flux of 45.3 mmol m −2 d −1 (de Froe et al., 2019). In contrast, the Haas Mound, which contained a lower live coral framework coverage near the AEC deployment site ( Figure 9A) showed an O 2 flux of 11.5-22.4 mmol m −2 d −1 (de Froe et al., 2019). This indicates that live coral patches are hotspots of metabolic activity within the reef and highlights the importance of understanding the spatial distribution of live and dead coral framework and sediments when planning research equipment deployments.
Similar to the Mingulay Reef, dead coral framework at the LMP contributes to the majority (49%) of the C turnover (De Clippele et al., 2021). This is expected as cold-water coral carbonate mounds consist of predominantly dead coral framework as discussed above. At the Mingulay Reef, the fauna associated with dead coral framework contributes ∼6 times more to C turnover compared to live coral framework. This is a larger difference compared to the LMP where the fauna associated with dead coral framework contributes only ∼1.7 times more. At Mingulay Reef, a higher biomass of fauna grows on the dead coral framework, hence turning over more C than the live corals (Kazanidis et al., 2016;De Clippele et al., 2021). This biomass difference may be caused by the reefs' shallower depth and the higher surface primary production above Mingulay Reef (0.048 g C org m −3 d −1 ), compared to the LMP (0.0067 g C org m −3 d −1 ) (Tyberghein et al., 2012;Assis et al., 2018;De Clippele et al., 2021) ( Table 5). This results in an annual PP over the LMP of 6,194 T C year −1 and 294 T C year −1 at the Mingulay Reef (Tyberghein et al., 2012;Assis et al., 2018;De Clippele et al., 2021).

How Much Organic Matter Is Required to Sustain the Deep Reefs?
From the annual C and N turnover of the coral presence habitat in the LMP area, we estimate a minimum annual C requirement of 5,763-7,124 T C year −1 and 973-3,194 T N year −1 . Using the parametrisation by Suess (1980), the amount of particulate organic carbon reaching the seafloor from the sea surface (annual primary production: 6,194 T C year −1 ) via deposition was estimated to be 511-322 T C year −1 for the shallowest (500 m) and deepest point (800 m) of LMP The reported values of carbon (C) turnover are the mean values. The total C turnover at the Mingulay Reef is based on data from live Lophelia pertusa, the sponge Spongosorites coralliophaga and Aquatic Eddy-Correlation (AEC) data of the dead framework. The AEC measurements capture the oxygen consumption of the dead framework and baffled sediment community. The results on C turnover for Logachev include data of live and dead L. pertusa, sediments and baffled sediments (see Section "Oxygen Consumption and Nitrogen Release"). reef habitat area, respectively. This indicates that almost the entire primary production, i.e., 91-95% (5,252-6,802 T C year −1 ), would have to be supplied through tidal downwelling, nepheloid layers, lateral deep-water advection and/or by the topographically-enhanced carbon pump (Duineveld et al., 2004;White et al., 2005;Mienis et al., 2007;Soetaert et al., 2016).
Our study suggests that the C requirement of the reef could be higher than the yearly PP over the area of 6,194 T C year −1 (Tyberghein et al., 2012;Assis et al., 2018) depending on seasonal PP and/or biomass variability. This could mean that the PP right above the area might not be sufficient to sustain the reef and highlights the importance of the supply of food trough advection from the wider area, bottom currents together with material retention and recycling of waste material on the reef, in particular during winter food limitation (Maier et al., 2020(Maier et al., , 2021. For example, studies have indicated that the reef could benefit from nitrification (re-utilization) of faunal-produced NH + 4 (Maier et al., 2020(Maier et al., , 2021 and utilization of dissolved OM, which is produced by the corals as mucus (Wild et al., 2004). The dependence of the reef 's function on these alternative supply mechanisms appears greater at the LMP compared to the Mingulay Reef (De Clippele et al., 2021) (Table 5), and is likely due to their location at greater depths with comparatively lower food flux. The supply of food needed to sustain the reef could be severely impacted by climate-induced changes in primary production, local hydrodynamical food supply, which could have severe implications for the survival and functioning of CWC reefs.

Conclusion
Biomass maps can guide sampling and monitoring expeditions and our current approach can be applied to other habitats, to provide large-scale maps of biomass, hotspots of metabolic activity and nutrient mineralization, in particular in the understudied, but large deep-sea realm. The predictive power of this approach can be improved by adding more coral surface area data, especially where the coefficient of variation of the map is higher. Additional local measurements on nutrient cycling, high resolution multibeam data (De Clippele et al., 2019) or more environmental variables (e.g., hydrodynamics) and the use of photo mosaics (Bodenmann and Thornton, 2017;Conti et al., 2019;Price et al., 2019) could further improve our understanding how complex habitats contribute to nutrient cycling. Biomass maps can also advice on the most optimal locations to collect AEC respiration data, to ensure a representative amount of habitat complexity is captured in the measurements (Rovelli et al., 2015;De Clippele et al., 2021). Alternatively, AEC deployments could be used to ground truth the biomass maps. Understanding how much dead and live coral framework is present in this area is especially important in deeper reefs such as the LMP, where ocean acidification threatens to dissolve dead coral framework (Hennige et al., 2015). If the dead coral framework dissolves, not only will the habitat of CWC reef organisms disappear (Kazanidis et al., 2016;Maier et al., 2021), but C and N demineralization, and sediment baffling will diminish. This may ultimately reduce primary production in surface waters, affecting the CO 2 being extracted from the atmosphere. Consequently, these effects will negatively impact the existence of the CWC mounds in the LMP and the overall functioning of the area.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.