Keywords

1 Introduction

The implementation of the Paris agreement to keep global warming below 1.5–2.0 °C is accompanied by a variety of measures such as the compilation of the NIR to monitor the emission of CO2 and other greenhouse gases in order to review the realization of climate pledges at national levels (e.g. UNEP 2019). The National Greenhouse Gas Inventory Report (NIR) splits CO2 emissions into four sectors: (I) energy, (II) industrial process and product use (IPPU), (III) waste, as well as (IV) agriculture, forestry, and other land use (AFOLU). In order to achieve a comparability of greenhouse gas emissions among different countries, the IPCC provides guidelines to quantify greenhouse gas emissions within these sectors (IPCC 2006). According to the South African, Namibian, and German NIR’s the first three sectors are net sources of CO2, while AFOLU acts as a CO2 sink to the atmosphere in all three countries (Table 25.1, German Environment Agency 2020; Department of Environmental Affairs 2019; Government of Namibia 2021).

Table 25.1 CO2 Emissions from South Africa, Namibia, and Germany according to the NIRs of these countries

AFOLU quantifies CO2 emissions caused by land use and land cover changes and is based on the quantification of net changes in carbon stocks. It considers six land-use categories. Among them are wetlands including also coastal ecosystems such as mangrove forests, tidal marshes, and seagrass meadows (IPCC 2014). Carbon stocks of these coastal wetland reservoirs are named ‘blue carbon’ due to the colour of the ocean (e.g. Nellemann et al. 2009). Despite being linked by name to the marine environment, blue carbon ignores the CO2 uptake by phytoplankton and its storage within the ocean’s biological carbon pump.

In terms of carbon fixation, phytoplankton in the ocean are as productive as terrestrial plants with a global rate of about 50 Pg C year−1 (Field et al. 1998), but the storage of carbon differs in systems dominated by terrestrial plants and marine phytoplankton. Terrestrial plants are the world’s largest reservoir of living biomass with a reservoir size of 450 Pg C, whereas phytoplankton represent, with 1–2 Pg C, a relatively small carbon stock (Bar-On et al. 2018; Falkowski et al. 1998). Instead of building-up a huge biomass, carbon fixed by phytoplankton is exported into the deep ocean where the vast majority of the exported biomass is respired and stored as dissolved inorganic carbon (DIC). This results in a strong gradient between low DIC concentrations in surface waters where carbon is fixed via photosynthesis into biomass, and high DIC concentrations at greater depth where the exported organic matter is remineralized. The transfer of DIC into the deep sea favours the CO2 uptake by the ocean by decreasing the CO2 concentration and therewith the partial pressure of CO2 (pCO2) in surface waters. A hypothetical collapse of the biological carbon pump is assumed to increase the atmospheric CO2 concentrations by 200–300 ppm, representing a doubling of the pre-industrial atmospheric CO2 concentration (Heinze et al. 2015).

The CO2 uptake efficiency of the biological carbon pump is linked to the global thermohaline circulation of the ocean (Heinze et al. 1991; Tschumi et al. 2011) and believed to have played a significant role in controlling atmospheric CO2 concentrations during glacial/interglacial transitions (Broecker and Barker 2007; Bauska et al. 2016; Schmitt et al. 2012). Even though it is widely assumed that the biological carbon pump also responds to the current climate change (e.g. Laufkötter et al. 2017; Duce et al. 2008; DeVries and Deutsch 2014; Riebesell et al. 2007) and fisheries (Bianchi et al. 2021), neither the magnitude nor the direction of change is predictable (Passow and Carlson 2012; Laufkötter et al. 2017; Laufkötter and Gruber 2018). Considering their potential impact on the CO2 concentration in the atmosphere and the Paris goals, it appears as a necessity to strongly reduce these uncertainties. This means that variations of the CO2 uptake by the biological carbon pump should be quantified and respective drivers should be identified, which, in turn, provides background to the discussion on responsibilities (national versus the international community). Here we use data obtained within the BMBF (Federal German Ministry for Education and Research) funded project TRAFFIC (Trophic tRAnsfer eFFICiency) to develop concepts that help to describe the status of the biological carbon pump and to quantify changes of the CO2 uptake within the BUS.

2 Study Area: The Benguela Upwelling System

The BUS stretches from the Angola Benguela Frontal Zone (ABFZ) at ~15°S to Cape Agulhas (~35°S, Fig. 25.1). The south-easterly trade winds emanating from the interplay between the South Atlantic Anticyclone and the continental low-pressure trough, cause the emergence of distinct upwelling cells along the coast (e.g. Kämpf and Chapman 2016; Veitch et al. 2009; Shannon and Nelson 1996). The Lüderitz Upwelling Cell at 26°40′S is the strongest of these cells and divides the BUS into a northern (NBUS) and southern (SBUS) subsystem (e.g. Hutchings et al. 2009; Duncombe Rae 2005; Shannon 1985).

Fig. 25.1
A map of the South Atlantic Ocean currents. The arrows indicate the direction of the currents. The major currents indicated here are the Benguela current flows northward, the South Atlantic current flows eastward, and the Angola current which flows southward.

Southeast Atlantic Ocean with its main currents as adopted and simplified from Verheye et al. (2016) and Hardman-Mountford et al. (2003) as well as annual mean sea surface temperatures (Smith et al. 2008). ABFZ (Angola Benguela Frontal Zone) and ‘Subtropical Convergence’ (dotted line) mark major oceanographic fronts. SECC South Equatorial Counter Current, PUC Poleward Under Current

Sub-Antarctic Mode Water mainly feeds upwelling (Marinov et al. 2006; Sarmiento et al. 2004), but its properties vary due to regional differences at sites of its formation and mixing with other water masses. The South Atlantic Central Water (SACW) is comprised of Sub-Antarctic Mode Water that originates as mixture of mainly Antarctic Intermediate Water and Subtropical Mode Water (Souza et al. 2018; Karstensen and Quadfasel 2002; McCartney 1977). It is subducted beneath subtropical surface waters north of the Sub-Antarctic Front at about 36°S–54°S and is fed via the South Atlantic Current into the Benguela Current (e.g. Gordon 1981; Donners et al. 2005). In the southern part of the SBUS, the Benguela Current follows the shelf break towards the equator, while the Benguela Jet, also known as the Cape Jet, flushes the SBUS shelf and carries Indian Ocean Central Water into the SBUS (Durgadoo et al. 2017; Tim et al. 2015). The mixture between South Atlantic Central Water and Indian Ocean water creates a new water mass, which is referred to as the Eastern South Atlantic Water (ESACW, e.g. Liu and Tanhua 2021). The ESACW as well as the SACW which largely encompasses the SBUS along with the Benguela Current, feed the complex equatorial current system. It comprises a number of east- and westward flowing currents and undercurrents (e.g. Pitcher et al. 2021) and finally feed the Angola Current from where modified SACW enters the NBUS as a poleward undercurrent at the ABFZ (Fig. 25.1). During its voyage through the south and equatorial Atlantic Ocean the SACW gets enriched in nutrients and depleted in oxygen (Mohrholz et al. 2008).

South of Lüderitz, the interplay between the Benguela Current, the Benguela Jet, and coastal upwelling creates a complex and temporally varying frontal system starting with the Oceanic Front (e.g. Veitch et al. 2010; Hardman-Mountford et al. 2003). This front is considered as the outer boundary of the SBUS (Barange and Pillar 1992; Shannon and Nelson 1996) whereas the Shelf Break Front and Upwelling Front develop in response to the variable nature of upwelling across the wide shelf regions (Shannon 1985).

During the last two decades, satellite data indicate no significant trends in productivity (Demarcq 2009; Lamont et al. 2019; Verheye et al. 2016), although there is a tendency towards an intensification of upwelling in the SBUS (Lamont et al. 2018, 2019; Tim et al. 2015; Wang et al. 2015; Sydeman et al. 2014). Associated with this trend are shifts in the ecosystem structure such as an overall decrease in zooplankton abundance and a shift from larger to smaller zooplankton species (Lamont et al. 2019; Verheye et al. 2016; Jarre et al. 2015; Hutchings et al. 2012).

3 Background

3.1 Nutrient Recycling and Productivity

The recycling efficiency of nutrients strongly influences the productivity whereas in general one can distinguish between two different nutrient recycling machineries in pelagic ecosystems. One recycling machinery operates in the seasonal thermocline, while the other one is located in the sunlit surface ocean (euphotic zone) where it directly affects primary production. Eppley and Peterson (1979) divided primary production (CO2 fixation via photosynthesis) into new and regenerated production. Upwelled nutrients impel new production, while the recycling of nutrients which have been introduced from the dark deep ocean into the euphotic zone fuels the regenerated production. The ratio between new production and primary production defines the f-ratio. An increasing f-ratio indicates a higher contribution of new production to the primary production and hence, a less efficient recycling of nutrients in the euphotic zone. Vice versa, a decreasing f-ratio reflects an enhanced regenerated production and thus a more intense recycling of nutrients. Nevertheless, nutrient recycling does not prevent the loss of nutrients via the export of organic matter from the euphotic zone. Over an annual cycle this so-called export production is assumed to equal new production. Even though export production represents a loss of nutrients from the euphotic zone, the exported organically bound nutrients are not necessarily lost to the pelagic ecosystem as export production drives a secondary recycling machinery in the seasonal thermocline.

The seasonal thermocline defines the subsurface layer from which water is introduced into the euphotic zone on the seasonal times scale which means that the more exported organic matter is remineralized within this subsurface layer, the more formerly exported nutrients return into the euphotic zone and can be recycled (Rixen et al. 2019a). In a coastal upwelling system which is characterized by onshore flowing subsurface waters (upwelling source waters), upwelling at the coast and offshore flowing surface waters, the secondary recycling machinery is part of the nutrient trapping system provoking new production by increasing nutrient concentrations in the upwelling source waters (Dittmar and Birkicht 2001; Tyrrell and Lucas 2002; e.g. Barange and Pillar 1992; Flynn et al. 2020; Rixen et al. 2021). Hereby, the utilization of upwelled nutrients and the associated development of plankton blooms in the offshore flowing upwelled waters initiate the nutrient trapping by exporting the formerly upwelled and now organically bound nutrients into the subsurface waters. This reduces their advection into the open ocean along with offshore flowing surface waters and keeps them within the upwelling system. Although these nutrients increase new production they come from the remineralization of organic matter as those nutrients which drive the regenerated production in the euphotic zone. Hence, these subsurface nutrients are also called regenerated nutrients. Since in addition to nutrients, CO2 is also released during the remineralization of organic matter, primary production driven by the utilization of regenerated nutrients consumes the associated regenerated CO2 and creates no need for additional CO2 as long as the Redfield carbon to nutrient uptake and remineralization ratio is constant. Hence, even though recycling of nutrients in the euphotic zone and the seasonal thermocline contributes to an elevated productivity of upwelling systems through their impact on regenerated and new production, respectively, it hardly affects the CO2 storage of the biological carbon pump. This differs with the utilization of preformed nutrients.

3.2 Preformed Nutrients and the CO2 Uptake Efficiency of Biological Carbon Pump

In contrast to regenerated nutrients which are released during the remineralization of exported organic matter, physical processes carry preformed nutrients into the deep ocean. This also occurs during the Sub-Antarctic Mode Water formation in winter when the lack of light prevents the utilization of regenerated nutrients in surface waters. During the subduction of mode water, the biologically unused preformed nutrients are transported into the deep ocean (Knox and McElroy 1984; Sarmiento and Toggweiler 1984; Siegenthaler and Wenk 1984; Duteil et al. 2012; Ito and Follows 2005). The regenerated CO2 formerly associated with the regenerated nutrients, is, in turn, released from the biological carbon pump. Vice versa, the biological carbon pump takes up CO2 by the utilization of preformed nutrients and their retransformation into regenerated nutrients at sites such as the BUS where upwelling introduces mode waters into the euphotic zone. To which extent the CO2 release and uptake affects the CO2 flux across the air–sea interface depends on its influence on the solubility pump. It could absorb CO2 released from the biological carbon pump or favour its emission into the atmosphere.

3.3 Solubility Pump and Upwelling from a Carbon Cycling Perspective

The solubility pump is an abiotically driven carbon pump which responds to the partial pressure difference of CO2 (ΔpCO2) between the surface ocean and the atmosphere. The net flux of CO2 follows the pressure gradient, which means that CO2 invades the ocean if the pCO2 in the atmosphere exceeds the oceanic pCO2 and vice versa the ocean emits CO2 into the atmosphere if the pCO2 in the ocean is higher than it is in the atmosphere. The CO2 concentration and its temperature- and salinity-controlled solubility determines the pCO2 in the ocean (Weiss 1974). Hence, CO2 concentrations link the biological carbon pump and the solubility pump, which gains strength in cold waters due to an increased solubility of CO2 at low temperatures. If, for example, ocean currents carry warm tropical waters with a pCO2 exceeding those in the atmosphere into polar regions, the pCO2 could fall below those in the atmosphere. Consequently, the water takes up CO2, simply due to decreasing sea surface temperatures. Vice versa, if water masses which are formed in the Southern Ocean such as the SACW, upwell in the BUS, they warm up and release CO2 that was previously taken up in the Southern Ocean. Hence, the opposing impacts of the biological carbon pump and the solubility pump on the CO2 concentrations during the Sub-Antarctic Mode Water formation and upwelling of mode waters in the BUS are parts of one system.

4 Data and Methods

In order to study the biological carbon pump within the framework of TRAFFIC, two cruises were carried out with the German research vessels Meteor (M153: February, 19th–March 31st 2019) and Sonne (SO285: August 20th–November 2nd 2021, Fig. 25.2). Mooring operations were additionally carried out during the RV Sonne cruise SO283, the RS Algoa cruise ALG 269, and other pre-TRAFFIC cruises with the German RVs Meteor and M.S. Merian. In addition to sediment trap samples, we analysed water samples to characterize source water masses and the chemical composition of plankton. This allows us to quantify the supply and utilization of preformed nutrients as well as Redfield ratios, which are crucial to translate nutrient consumption into carbon uptake. Furthermore, satellite data were evaluated and linked to sediment trap data in order to estimate the reliability of export production rates. Satellite-derived sea surface temperature (SST) and primary production rates (PP) were downloaded from the MODIS-Aqua website in August 2021 (https://oceandata.sci.gsfc.nasa.gov/MODIS-Aqua/Mapped/Monthly/4km/sst/) and the ocean primary production website in August 2020 (http://www.science.oregonstate.edu/ocean.productivity/).

Fig. 25.2
A bathymetric map of the coast of California exhibits the major features of the ocean floor, continental slope, and the abyssal plain. The major locations indicated here are some major cities and ports along the coast of California.

Working area and sampling sites during the two TRAFFIC cruises M153 and SO285. EEZ = 200-mile Exclusive Economic Zones, CT indicates Cape Town. Please note that the timing of the cruises indicates the time within the working area

4.1 Particulate Matter

Particulate matter was collected by moored and drifting sediment traps, as well as through the filtration of water samples and by net catches.

A sediment trap mooring was deployed off Walvis Bay (14.04°E/23.02°S) and equipped with a Hydrobios MST-12 sediment trap at water-depths between 65 m and 75 m (~74 m above the sea floor). Since the moored TRAFFIC traps have yet not fully been recovered, we focus on the pre-TRAFFIC deployment periods. This includes seven deployment periods, which are in part not time-coherent. The first deployment started in December 2009 and the last deployment ended in November 2017. Furthermore, four short-term drifting sediment trap systems were deployed in the NBUS and SBUS for about 41–83 h during the cruise M153 (Fig. 25.2). These systems were equipped with four to five single sediment traps and a Hydrobios MST-12 trap at the bottom of the drifter at water-depths of down to 500 m. In total, 24 drifter samples were collected and analysed in Germany.

According to common sediment trap processing procedures, sediment trap samples were split into fractions of >1 mm and <1 mm (Haake et al. 1993; Honjo et al. 2008). The >1 mm fraction contained larger swimmers and the <1 mm fraction was assumed to represent the gravitationally driven export of particles (= particle flux). In addition, the 24 drifter samples were additionally macroscopically and microscopically examined and the picked zooplankton was divided into three genera, namely copepods, amphipods, and euphausiids. Furthermore, we obtained zooplankton and fish samples from net catches at the deployment stations during the cruise M153 (Fig. 25.2). These samples were fully homogenized and analysed.

Water samples (5 L and 30 L) obtained from the CTD rosette water sampler were filtered (WHATMAN GF/F; ~0.7 μm; 47 mm diameter) and rinsed with deionized water. Subsequently, filters were dried, shipped to Germany, and analysed in accordance to the sediment trap samples. A detailed description of the sediment trap sample procedure including the biogeochemical analysis is given by Haake et al. (1993) and Rixen et al. (1996). The analysis included the determination of total carbon and nitrogen and organic carbon. The phosphorus content was additionally determined in zooplankton samples according to a method modified from Aspila et al. (1976) and Grasshoff et al. (1999).

4.2 Water Sampling

The profiling SeabirdElectronics (SBE) 911plus CTD system was equipped with a Digiquartz pressure sensor, a SBE3 temperature and SBE4 conductivity sensor within a double sensor setup, and a SBE43 dissolved oxygen sensor. Additionally, water samples were taken from the CTD rosette water sampler to analyse DIC concentrations as well as the total alkalinity (TA) and nutrient concentrations. Therefore, we used a cavity ringdown spectrometer (Picarro G2201-I, 1510CFIDS2047_v1.0) attached to a Liaison A0301 (DIC), a VINDTA 2C analyser (MARIANDA, Kiel, Germany, TA), and a continuous flow injection system (Skalar SAN plus System). Detection limits for nitrate/nitrite (NOx) and phosphate were 0.08 μM and 0.07 μM, respectively (Flohr et al. 2014). Both the Picarro and the VINDTA 2C analyser were calibrated using Certified Reference Material (CRM, batch #177) provided by A. Dickson (Scripps Institution of Oceanography, La Jolla, CA, USA). The accuracy was ±12 μmol kg−1 and ± 4.3 μmol kg−1 for DIC and TA, respectively.

5 Results and Discussion

5.1 CO2 Concentrations

Within the BUS, CO2 concentrations were high in a narrow belt along the coast (Emeis et al. 2018) where ESACW and SACW enriched in nutrients and CO2 reach the surface. Plankton blooms developing in response to the nutrient input consume nutrients and CO2. Thereby they decrease the concentration of both components in the offshore flowing surface water until nitrate is consumed and phosphate concentrations drop to values of ~0.2 μmol kg−1 (Flohr et al. 2014). This so-called excess phosphate is assumed to be the result of a nutrient uptake which follows the Redfield ratio and a relative enrichment of phosphate over nitrate due to anaerobic processes within the Oxygen Minimum Zone (OMZ) on the shelf and surface sediments (Flohr et al. 2014; Mashifane 2021; Nagel et al. 2013; Goldhammer et al. 2010). However, after the consumption of nitrate within offshore flowing surface waters, a balance between nitrate consumption and release during the respiration of organic matter prevents nitrate accumulation in the euphotic zone. Since such a regenerated nutrient cycle hardly affects CO2 concentrations, the variability of pCO2 is strongly reduced at approximately 340 km from the coast. This can be seen in the BUS (Siddiqui et al. 2023) as well as in the California Upwelling System (Chavez and Messié 2009) and marks in addition to nitrate concentrations below the detection limit, the diminishing influence of upwelling on the pelagic ecosystem. At these outskirts of the upwelling-driven ecosystem, CO2 concentrations reflect the balance between the CO2 uptake via the utilization of preformed nutrients (biological carbon pump) and the warming of the upwelling water (solubility pump). Vice versa, changing CO2 concentrations at such sites indicate varying intensities of these two pumps within the upwelling system. Hence, moored CO2 observing buoys, including oxygen, nitrate, and pH sensors could help to detect long-term changes in the balance between the biological carbon pump and solubility pump. Since the region influenced by upwelling varies and is fragmented by mesoscale eddies and filaments (e.g. Rubio et al. 2009), it is suggested to deploy such CO2 observing buoys along transects perpendicular to the coast encompassing the region in which nitrate depletion and low variability of pCO2 indicate the diminishing effect of upwelling on pelagic ecosystems.

5.2 Export Production

Export/new production can be used in conjunction with satellite-derived primary production rates to characterize the recycling efficiency in the euphotic zone as indicated by the f-ratio. Combined with the supply of preformed nutrients along with upwelling source water masses, it is also a parameter that allows us to estimate the CO2 uptake by the biological carbon pump as we will discuss in the following sections. However, so far there are no reliable methods to measure export production directly. Due to uncertainties involved in methods commonly applied to determine export production, estimates of global mean export production rates reveal a wide range with values between 1.8 and 27.5 Pg C year−1 (Lutz et al. 2007; del Giorgio and Duarte 2002; Honjo et al. 2008). A widely accepted global mean export production is centred around 10 Pg C year−1 (Emerson 2014: and references therein) resulting in a mean area normalized export production rate of 28 g C m−2 year−1 (area of the ocean 360 1012 m2) and a global mean f-ratio of 0.2.

5.2.1 Export Production: A Top-Down Approach

In contrast to the bottom-up approach in which nutrient concentrations are multiplied by upwelling velocities to quantify new/export production (e.g. Chavez and Messié 2009; Messié et al. 2009; Waldron et al. 2009), other methods address the issue of export production from a top-down approach by looking at primary production. Thereby, export production is described as a simple function of primary production. Eppley and Peterson (1979) introduced one of the first functions which assumes that export/new production (POCExport) contributes 50% (f-ratio 0.5) to the total primary production (PP) in high productive systems (see Eq. 25.1). Laws et al. (2000) and Henson et al. (2011) additionally considered impacts of sea water temperatures (SST) on the link between export and primary production as shown in Eqs. (25.2) and (25.3).

$$ {\textrm{POC}}_{\textrm{Export}}{=}_{0.5\cdot \textrm{PP},\kern0.75em \textrm{if}\ \textrm{PP}>200}^{0.0025\cdot {\textrm{PP}}^2,\kern0.75em \textrm{if}\ \textrm{PP}<200} \vspace*{-6pt} $$
(25.1)
$$ {\textrm{POC}}_{\textrm{Export}}=\left(-0.02\cdot \textrm{SST}+0.63\ \right)\cdot \textrm{PP} \vspace*{-6pt} $$
(25.2)
$$ {\textrm{POC}}_{\textrm{Export}}=0.23\cdot {\exp}^{\left(-0.08\cdot \textrm{SST}\right)}\cdot \textrm{PP} $$
(25.3)

In the BUS, SSTs vary in general between <20 °C at the ABFZ and about 10 °C at sites where upwelling source waters reach the surface. Within this temperature range, f-ratios (= POCExport/PP) as derived from Eqs. (25.2) and (25.3) are far below 0.5, with means of 0.3 ± 0.06 and 0.07 ± 0.02, respectively (Fig. 25.3). This implies a much stronger nutrient recycling in the euphotic zone and a lower export production at the same primary production. We used these equations in addition to satellite-derived primary production rates and SST to calculate export production rates, which were compared to organic carbon fluxes measured by a sediment trap deployed in the NBUS (Figs. 25.2 and 25.4).

Fig. 25.3
A graph for the f-ratio versus sea surface temperature reveals 4 lines for the E q 25.1 from (10, 0.5) to (20, 0.5), E q 25.2 from (10, 0.43) to (20, 0.23), E q 25.3 from (10, 0.1) to (20, 0.04), and global mean from (10, 0.2) to (20, 0.2). Values are approximate.

f-ratio versus sea surface temperatures according to Eqs. (25.1), (25.2), and (25.3) as well as in comparison to the global mean f-ratio (dotted line) and those determined by our SBUS drifter study (grey shaded area)

Fig. 25.4
A multi-line graph exhibits data for milligram centimeter negatively squared day superscript negative 1 versus years from 2010 to 2018. The line for E q 25.1 plots the highest peak in 2012 at 3500.

Export production rates derived from Eqs. (25.1)–(25.3) (top-down approaches) as well as from bottom-up approaches (this study and Messié et al. 2009) which are based on upwelling velocity, nutrient concentrations in source waters, and a fixed carbon to nutrient ratio. Open circles indicate organic carbon fluxes measured by moored sediment trap off Walvis Bay in the NBUS at water-depth between 60 m and 74 m (Trap, see Fig. 25.2) while black circles indicate the corrected sediment trap data. The sediment trap data have been corrected for remineralization below the surface mixed layer by assuming that 78% of the exported organic matter is decomposed between the base of the surface mixed layer and the trap (Trap 78%)

The sediment trap was deployed in a water-depth of 60–75 m close to the Namibian coast within the narrow belt along the coast where CO2 concentrations and primary production are high. According to satellite data, the mean primary production rate over the period of the sediment trap experiment from 2009 to 2018 was about 3000 mg m−2 day−1. This is high, but within the upper range of primary production rates measured in situ during expeditions (Barlow et al. 2009; Wasmund et al. 2005). Export production rates derived from the three equations vary (Eq. 25.1 = 1794 ± 426 mg m−2 day−1; Eq. 25.2 = 1101 ± 295 mg m−2 day−1) due to the different f-ratios, whereas results obtained from Eq. (25.3) (229 ± 64 mg m−2 day−1) are similar to organic carbon fluxes measured by sediments traps (136 ± 90 mg m−2 day−1). However, since the sediment trap was deployed at a water-depth of 60–75 m within the thermocline where exported organic matter is respired as indicated by low oxygen concentrations, trap results do not reflect export production rates but only the export of organic matter at water-depth 60–75 m. Hence, the difference between the export production and the measured organic carbon flux indicates the organic matter remineralization between the base of the euphotic zone and the trap-depth. Data obtained in the SBUS allowed us to study organic matter decomposition at such low water-depths in more detail.

During the cruise M153 (February, 19th–March 31st 2019) a drifter was deployed close to the South African coast, with sediment traps attached at four different water-depths (20, 30, 50, and 75 m) at one of our 24 hour-stations. At this station in situ primary production was also determined by using measured photosynthesis-light-curves (PE curves) in combination with incubation experiments at in situ temperatures and different light intensities. Furthermore, CTD casts were conducted four times a day in order to capture a diurnal cycle and during the cruise SO285 (August 20th–November 2nd 2021) this site was revisited (Fig. 25.5). The drifter traps showed that the total flux of organic carbon was highest just below the surface mixed layer at a water-depth of 20 m (1588 mg m−2 day−1) and decreased step-wise with increasing water-depth. At a water-depth of 75 m, the ‘20 m-flux’ had already decreased by 78% to 358 mg m−2 day−1. Measured in situ primary production reached values of about 2430 mg m−2 day−1 and were similar to monthly mean satellite-derived primary production rates of 3428 mg m−2 day−1 (March, 2019). Primary production of 2430–3428 mg m−2 day−1 and an export production of 1588 mg m−2 day−1 result in a f-ratio of 0.65 and 0.37, respectively. f-ratios used in Eq. (25.1) and obtained by Eq. (25.2) partly fall in this range whereas f-ratios derived from Eq. (25.3) seem to be too low (Fig. 25.3).

Fig. 25.5
2 multi-line graphs. A, water depth versus C o r g per P P for temperature plots fluctuating and increasing curves. B, water depth versus organic carbon fluxes for the oxygen plots fluctuating curves for the time at 12, 16, 2, 8, and 5 hours.

Temperature and oxygen profiles obtained at the near shore drifter site in the SBUS (see Fig. 25.2) at different times of the day during the cruise M153 (March 2019, station 7) and during the cruise SO285 (September/October 2021, station 43) at 06:00 in the morning (S: 6:00). The red circles indicate the ratio between exported organic matter and primary production (Corg/PP, a) and organic carbon fluxes determined by the drifter study in four different water-depth. The dark and light grey shades indicate the depth of the surface mixed layer during the cruises M153 and SO285

Assuming that 78% of the export production is decomposed between the base of the surface mixed layer and a water-depth of about 75 m implies that the measured mean organic flux of 136 ± 90 mg m−2 day−1 at our NBUS trap site represents only 78% of the export production which, in turn, suggests an export production of 618 ± 410 mg m−2 day−1. If one applies this 78%-correction to all measured sampling intervals, it shows that there are several periods at which the corrected organic carbon export rates are similar to those derived from Eq. (25.1) and partly even Eq. (25.2) (Fig. 25.4). This by far does not solve the dilemma of determining export production rates nor the accuracy of sediment trap measurements, but merely increases the reliability of Eq. (25.1) and (25.2) as well as our drifter experiment.

The revisit of the SBUS drifter site during the cruise SO285 showed that the seasonally varying mixed layer depth is another factor that needs to be considered when comparing shallow water sediment trap results and satellite data. Even though a similar situation was observed during the cruises M153 and SO285 at this site (Fig. 25.5), the comparison of data obtained during these cruises also shows some pronounced differences. For instance, the surface mixed layer was about 14 m deeper, while bottom water oxygen concentrations were higher in September 2021 (SO285) than in March 2019 (M153). The intensity of the bottom water OMZ is known to reveal a pronounced seasonal cycle with decreasing oxygen concentrations during the main upwelling season in austral summer and increasing oxygen concentrations during the austral winter (Pitcher et al. 2014). Hence, the bottom water OMZ was more intense in March at the end of the upwelling season than in September at the beginning of the upwelling season. Vertical mixing is, in turn, assumed to supply oxygen to the intermediate water layer so that this layer maintains nearly constant oxygen concentrations over the years. Nevertheless, spatial variations of the Upwelling Front cause temporally limited events during which the oxygen concentration drop by up to ~150 μM (Rixen et al. 2021).

However, the 14 m winter deepening of the surface mixed layer could have strongly influenced results obtained by sediment traps which are deployed at a fixed water-depth. Assuming a trap-depth of 75 m, a winter deepening of the surface mixed layer by 14 m could have reduced the distance between the trap and the base of surface mixed layer by 30% from 60 m in March to 42 m in September. Considering the rapid decomposition of organic matter as observed in the drifter experiment, such a shortening of the distance between the trap and the base of surface mixed layer should have had a significant impact on the organic carbon fluxes. Hence, measurements of the mixed layer depth should be integrated into shallow water sediment trap studies to correct fluxes for changes in the distance between base of surface mixed layer and the sediment trap. Respective information could, e.g., be obtained by combining sediment trap studies with the deployment of CO2 observing buoys as discussed before.

5.2.2 Export Production: A Bottom-Up Approach

In accordance with Eppley and Peterson (1979) and their Eq. (25.1), Messié et al. (2009) determined with their bottom-up approach a f-ratio of 0.5 for the BUS and f-ratios of 0.4–0.7 for the other eastern boundary upwelling systems. This agrees with results obtained from a comprehensive regional model study suggesting a mean f-ratio of 0.58 (Emeis et al. 2018; Schmidt and Eggert 2016). This regional model was also used to quantify the amount of water (0.9 × 106 m3 s−1) which upwells in the NBUS off Namibia between 16°S and 28°S (Müller et al. 2014). Following the bottom-up approach we primarily used the M153 and SO285 data to determine the mean nutrient concentrations in the upwelling source waters (Table 25.2), which, subsequently multiplied with upwelling velocities derived from the model, results in the new/export production rate for the BUS.

Table 25.2 Mean properties of the upwelling source water between 18–28°S during the cruises M153 and SO285. Please note that nutrient concentrations have been measured in discrete samples during the cruise M153 and nitrate concentrations by the Seascan Suna Deep Nitrate Sensor during the cruise SO285. Other nutrient data from cruise SO285 are not yet available since this cruise was still on-going while this paper was written, n = number of samples and Std = standard deviation

SACW and ESACW, the principal source waters in the BUS, can be identified in θ-S-diagrams as they fall on distinct straight lines within a density range of 27.3 to 26.4 (Fig. 25.6). The comparison between cruises M153 and SO285 shows pronounced differences in the surface water properties. Summer warming decreased the density of surface water and increased the density gradient between surface and upwelling source waters in March, whereas winter cooling had an opposing effect (Fig. 25.6). It decreased the density gradient between surface waters and upwelling source water masses. In contrast to the surface water, the physical properties of upwelling source waters with this density horizon showed hardly any variations except in concentrations of dissolved oxygen (Table 25.2). In the NBUS, they were lower in March than in September due to an enhanced inflow of oxygen-poor SACW between March and May (Mohrholz et al. 2008).

Fig. 25.6
4 graphs with a dual vertical axis for the potential temperature and oxygen versus a horizontal axis for salinity plots data for N BUS M 153, S BUS M 153, N BUS S O 285, and S BUS S O 285 in parts A to D.

θ-S-diagrams including concentrations of dissolved oxygen derived from CTD casts in NBUS and SBUS during the cruises M153 (February/March 2019) and SO285 (September/October 2021). The red and blue straight lines indicate the θ-S characteristics of the ESACW and SACW according to the equation given by Flohr et al. (2014)

In comparison to oxygen the seasonal variability of nitrate was low (Table 25.2). A mean nitrate concentration of 23.5 μmol kg−1 (= (22.59 + 24.39)/2) and a phosphate concentration of 1.75 μmol kg−1 multiplied by the mean upwelling velocity of 0.9 × 106 m3 s−1 (~0.93 × 106 kg s−1) suggest nutrient inputs into the surface layer of 52.5 × 109 mol phosphate and 718.5 × 109 mole nitrate. Multiplied by a carbon to nitrogen and carbon to phosphorus ratio of 7.6 and 106, which will be discussed in the following sections, these nutrient inputs suggest a new/export production rate of 65–67 Tg C year−1. Considering an area of 0.16 × 1012 m2 these values imply export production rates of 1113–1147 mg m−2 day−1. which is within the range of export production rates derived from Eq. 25.2 (1101 ± 295 mg m−2 day−1) and below those (1416 mg m−2 day−1) obtained by Messié et al. (2009). Multiplied by an area of 0.16 × 1012 m2 the export production of 1416 mg m−2 day−1 amounts to 83 Tg C year−1 which we will consider as upper estimate for the export production in the BUS between 18°S and 28°S in the following discussion.

5.3 Redfield Ratio

In 1934, Redfield introduced the concept of a constant carbon to nutrient ratio in the year 1934 and provided further support in a couple of follow-up publications (e.g. Redfield 1934, 1958; Redfield et al. 1963). He discovered that changes in the concentration of DIC, nitrate, phosphate, and oxygen in the water column reflect the stoichiometry of plankton. Today a C/N/P/-O2 ratio of 106/16/1/−138 is considered as the classical Redfield ratio. During the last 60 years, the number of observations increased enormously showing that the Redfield ratios derived from correlations of dissolved components reveal a remarkedly low variability (Anderson and Sarmiento 1994; Takahashi et al. 1985) in comparison to those derived from the elementary composition of plankton (Boyd and Trull 2007; Planavsky 2014; Martiny et al. 2013). Hence, we are again confronted with a discrepancy between bottom-up (dissolved components) and top-down (plankton stoichiometry) approaches.

5.3.1 Redfield Ratio: A Top-Down Approach

A main factor controlling the Redfield ratio is the ratio between cyanobacteria and eukaryotic plankton as cyanobacteria often show enhanced and highly variable Redfield ratios (Karl et al. 1997; Bertilsson et al. 2003; Sanudo-Wilhelmy et al. 2004). Hence, C/N/P ratios in oligotrophic subtropical gyres dominated by cyanobacteria are high, with C/N/P values of 195:28:1 (Martiny et al. 2013; Teng et al. 2014). However, in addition to variations in plankton community structure, environmental changes such as temperature and pH variations as well as nutrient stress could influence the carbon to nutrient ratios in individual plankton clades (Boyd and Trull 2007; Planavsky 2014; Martiny et al. 2013; Geider and La Roche 2002; Riebesell et al. 2007). Hence, organic matter in nutrient-rich environments and at high latitudes often reveal lower Redfield ratios (78/13/1, C/N = 6) than in nutrient-depleted regions (138/18/1, C/N = 7.7) at lower latitudes (Martiny et al. 2013; Teng et al. 2014).

This compilation also shows that in comparison to C/P ratios, the variability of C/N ratios is quite low (6–7.7) and hardly deviates from the classical Redfield ratio of 6.6. We obtained a similar result from plankton samples which were collected within the chlorophyll maximum during the cruise M153. The chlorophyll maximum is the zone within surface waters where the highest chlorophyll concentration occurs. C/N ratios determined from these samples vary between 5.6 and 8.6, with a mean of about 6.7 in the NBUS and SBUS. Zooplankton, fish, squid, and jelly fish samples collected by net catches during the cruise M153 reveal, in turn, lower C/N ratios of about 4.7 ± 0.7 (n = 37), suggesting an enrichment of nitrogen within higher trophic levels of the food web (Table 25.3).

Table 25.3 Mean C/N and C/P ratios with corresponding standard deviations (std) and n number of samples measured. The samples were obtained from net catches at locations marked as ‘M153 Zooplankton’ in Fig. 25.2

Our zooplankton samples include euphausiids, amphipods, and decapods, but excluded copepods. Copepods are abundant in the BUS as well as in our drifter samples and are known to reveal highly variable C/N ratios of 3.6–10.4 (Bode et al. 2015; Schukat et al. 2014). In order to study the role of zooplankton in more detail, we also deployed drifting sediment traps during the cruise M153 as mentioned before. The results show that the zooplankton carbon contributes on average 65 ± 21% to the total organic carbon collected by the traps with a mean C/N ratio of 7.8 ± 1.1. Considering the total organic carbon flux (zooplankton and the <1 mm fractions) results in a mean C/N ratio of 8.1 ± 1.8. This is similar to the mean C/N ratios of 8.7 determined by other sediment trap experiments in the region including traps deployed at water-depths down to 2500 m (Vorrath et al. 2018) and those C/N ratios measured in the moored NBUS traps (8.4 ± 1.3). Compared to phytoplankton (5.6–8.6, mean: 6.7), these slightly enhanced C/N ratios could be caused by a preferential decomposition of organically bound nitrogen in the water column or inputs of organic matter from land and resuspended sediments. On the other hand, we discovered in line with an earlier study along the European continental margin (Antia 2005) a preferential leaching of phosphorus containing organic components, which raised C/P ratios to values of >400. Hence, C/P ratios measured in trap samples have been ignored in this discussion. Along the European continental margin where traps had been deployed at water-depths between 600 and 3200 m, leaching of organically bound nitrogen raised the C/N ratio on average from 8.1 to 11.3 (Antia 2005). Even though a mean C/N ratio of 8.1 ± 1.8 as measured by our traps at a water-depth of 60–75 m implies a comparably low leaching of organically bound nitrogen, it could suffice to explain the slightly enhanced C/N ratio in the trapped exported organic matter. To which extent C/N ratios derived from concentrations of DIC, phosphate, and nitrate reflect those of particulate matter will be discussed in the next section.

Fig. 25.7
4 graphs with scattered dots around an increasing diagonal line plot data for the N BUS in parts A and B and S BUS in parts C and D for D I C versus P O 4.

DIC versus phosphate concentrations as well as phosphate versus nitrate concentrations measured in samples obtained in NBUS and SBUS during the cruise M153

5.3.2 Redfield Ratio: A Bottom-Up Approach

In the NBUS and SBUS, nitrate and phosphate concentrations are correlated, while the slope derived from the regression equation suggests mean N/P ratios of 13.7–14 (Fig. 25.7). The slightly lower N/P ratio in the NBUS (Fig. 25.7b) is most likely a consequence of nitrate reduction (Tyrrell and Lucas 2002). This occurs in the NBUS OMZ (e.g. Nagel et al. 2013), which is more intense than the SBUS OMZ. However, C/N and C/P ratios are more difficult to determine since the solubility pump accounts for the majority of the carbon dissolved in ocean waters. The y-axis intercept of the regression line derived from the correlation between phosphate and DIC concentrations indicates the DIC background concentrations of 2144 μmol kg−1 in the NBUS and 2085 μmol kg−1 in SBUS. These differences could be caused by the dissolution of carbonates, as well as varying pCO2 disequilibria between surface waters and the atmosphere and the uptake of anthropogenic CO2 during the formation of ESACW and SACW. The slopes, in turn indicate C/P ratios of 81.1 and 86.7 in the NBUS and SBUS, respectively.

Fig. 25.8
A graph of D I C corrected versus P O 4 plots 3 diagonal lines for C over P =93 and redfield 106 from (0.0, 0) to (2.8, 250) and (2.6, 250), respectively. Several scattered dots are around the diagonal lines.

Corrected DIC versus phosphate concentrations from the NBUS and SBUS. The red line marks the regression equation (DICcorrected = 92.9 × PO4 + 2.129, r = 0.782, n = 617) and the blue line shows show the Redfield ratio

Well-accepted and often used methods (Sabine et al. 1999; Gruber et al. 2019) associated with the quasi-conservative tracer ∆C* can be used to quantify the influence of the pCO2 disequilibria between ocean and atmosphere, the dissolution of carbonate and the uptake of anthropogenic CO2 on the DIC concentration (Gruber et al. 1996). But since these methods already rely on the assumption of a constant Redfield ratio as derived from Anderson and Sarmiento (1994) of 117/16/1/170 (C/N/P/−O2), they are not applicable here where we aim to determine Redfield ratios in the BUS. Nevertheless, deviation from the minimum TA was considered as a consequence of carbonate dissolution and the corresponding release of DIC (DICdissolution = (TAmeasured – TAminimum)/2) as well as the DIC background concentrations were subtracted from the measured DIC concentrations (DICcorrected = DICmeasured – DICdissolution – DICbackground). These corrections reduced the differences between NBUS and SBUS, and the slope of the resulting regression equation implied a C/P ratio of 93 (Fig. 25.8). C/P ratios varying between 81.1 and 93.0 combined with N/P ratios of 13.7–14.0 suggest a mean Redfield ratio of 87 ± 6/14 ± 0.2/1. The resulting C/N of about 6.3 falls slightly below the classical Redfield C/N ratio of 6.6 and C/N ratios determined by us in the chlorophyll maximum with a mean of about 6.7. Considering the chlorophyll maximum C/N ratios instead of the N/P ratio derived from nitrate and phosphate concentrations results in a mean Redfield ratio of 87 ± 6/13 ± 3/1. The differences between these two approaches are low and the overall relatively low Redfield ratios are similar to those derived from plankton in nutrient-rich environments and high latitudes (78/13/1). However, in February 2011, we determined with a bottom-up approach a Redfield ratio of about 100 ± 1/16 ± 2/1 off Namibia and also located a region on the Namibian shelf where anaerobic processes in association with nutrient fluxes across the sediment water interface reduced the Redfield ratio to 66/10/1 (Flohr et al. 2014). Hence, Redfield ratios seem to vary interannually and could be influenced by sediment water interactions within pronounced bottom water OMZ. Nevertheless, regarding the CO2 uptake by the biological carbon pump, our bottom-up and top-down approach based on the M153 data support Redfield’s findings and imply that the consumption of upwelled nitrate and phosphate adsorbs as much CO2 as is released during the respiration of exported organic matter. This, in turn, poses the question as to how the carbon export could influence the CO2 uptake by the biological carbon pump.

5.4 CO2 Sequestration

In order to estimate carbon sequestration by the biological carbon pump, we applied the bottom-up approach to calculate export production and used the Redfield ratio to convert phosphate supply into surface waters into carbon export production rates. In principle, this approach is similar to those used by Waldron et al. (2009) for the SBUS but in contrast to Waldron et al. (2009), we used only preformed and not the total nutrient concentrations of the source waters. As discussed before, preformed nutrients have been detached from CO2 so that their utilization acts as a CO2 sink. Since the preformed phosphate concentrations contribute about 30% to the total phosphate concentration (Table 25.2), the CO2 uptake by the biological carbon pump amounts to about one third of the total export production. In the BUS between 18°S and 28°S which approximately covers the Namibian part of the BUS, we estimated a total export production of 65–83 Tg C year−1, whereas Waldron et al. (2009) suggested a total export production of about 7–39 Tg C year−1 for the South African part of the BUS. Assuming that utilization of preformed nutrients accounts for 30% of the total export production suggests a CO2 uptake of 20–25 Tg C year−1 (mean = 22.5 ± 2.5 Tg C year−1) in the Namibian part of the BUS and 2–12 Tg C year−1 (mean = 7 ± 5 Tg C year−1) in the South African part of the BUS which is similar to the results of our previous study (Siddiqui et al. 2023).

In order to estimate possible impacts of changing Redfield ratio, we repeated the calculations but instead of the classical C/P ratio of 106, we used a C/P ratio of 87 ± 6 as determined by using the M153 data. The resulting effect can be calculated by the rule of three (export driven by preformed nutrients/106 × 87 ± 6) which lowered the mean CO2 uptake rates to 18.5 ± 3.3 Tg C year−1 and 6.0 ± 5.0 Tg C year−1 for the Namibian and the South African part of the BUS, respectively. This represents a decrease of the calculated CO2 uptake by the biological carbon pump of 17% (4 Tg C year−1) and 13% (1 Tg C year−1) off Namibia and South Africa.

The CO2 uptake rates by the biological carbon pump of 18.5 ± 3.3 Tg C year−1 and 6.0 ± 5.0 Tg C year−1 are in the same order of magnitude as the CO2 uptake by AFOLU of 30.6 Tg C year−1 (Namibia) and 7.5 Tg C year−1 (South Africa, see Table 25.1), but these carbon fluxes cannot be compared to each other since the latter represents human impacts on the CO2 uptake by terrestrial ecosystems. Even though marine ecosystems respond to global change (Lamont et al. 2019; Verheye et al. 2016; Jarre et al. 2015; Hutchings et al. 2012), fisheries (Bianchi et al. 2021) and in particular bottom trawling which will be discussed in the following section, effects on the CO2 uptake by the biological carbon pump via influences on the Redfield ratio have not yet been quantified. Nevertheless, the magnitude of the CO2 uptake rates by the biological carbon pump is on levels at which small changes due to, e.g., decreasing C/P ratios of 17% or 13%, respectively, could influence the CO2 uptake by AFOLU significantly.

5.5 Benthos

So far, we discussed the CO2 uptake by the biological carbon pump in the water column, but it also influences the burial of carbon in sediments. This represents a removal of carbon from the climate active carbon cycle and humans could favour and reduce it. An enhanced supply of clay and other minerals via dust inputs and river discharges due to soil as well as coastal erosion favours the organic carbon sedimentation by reducing the remineralization of organic matter in the water column and surface sediments mainly by two processes. Firstly, minerals incorporated into particles increase their sinking speed and the resulting accelerated transport of organic matter from the euphotic zone onto the sediments lowers the residence time of organic matter in the water columns and therewith its remineralization (e.g. Rixen et al. 2019b). Secondly, the adsorption of organic matter to mineral surfaces reduces its remineralization in the water column and surface sediments by protecting organic matter against bacterial attacks (e.g. Armstrong et al. 2002; Hedges 1977). Vice versa, bottom water trawling remobilizes the sedimentary organic matter and favours its return into the climate active carbon cycle. Sala et al. (2021) introduced an approach to quantify impacts of bottom trawling on the remobilization of carbon from marine sediments on a global scale. Assuming a continuous bottom trawling, Sala et al. (2021) suggested, for instance, a global long-term sedimentary organic carbon loss of about 158 Tg C year−1, and also broke it down to the BUS where it amounts to about 5 Tg C year−1. In comparison to AFOLU and our estimates of the CO2 uptake by the biological Compared to AFOLU and our estimates of the CO2 uptake by the biological carbon pump of about 24 Tg C/year-1 in the entire BUS (NBUS and SBUS), this is a substantial amount. It shows that humans can significantly affect the CO2 uptake of the biological carbon pump which still needs to be quantified.

6 Conclusion

In contrast to the often assumed and used constant global mean Redfield ratio, our results indicate varying Redfield ratios whereas the one determined by us during the TRAFFIC cruise M153 (87 ± 6/13 ± 3/1) is similar to those derived by other studies in nutrient-rich environments and high latitudes (78/13/1). CO2 uptake by the biological carbon pump based on the supply of preformed nutrients and the Redfield ratio amounts to 18.5 ± 3.3 Tg C year−1 and 6.0 ± 5.0 Tg C year−1 for the Namibian and South African parts of the BUS, respectively. These uptake rates are on a level at which small changes could influence the CO2 uptake by AFOLU significantly. Considering the release of sedimentary carbon by bottom trawling and that ecosystems in the BUS respond to global change and fisheries, it is quite likely that humans affect the CO2 uptake by the biological carbon pump already but uncertainties are still too large to quantify such impacts. Hence, we suggest to improve methods to estimate the supply of preformed nutrients, carbon export rates, and Redfield ratios and expand monitoring strategies by linking established observing methods such as remote sensing and sediment trap studies, and focus on CO2 observation at the outer rim of the upwelling system. This appears to be a key site to detect long-term changes in the balance between the biological carbon pump and solubility pump. Despite large uncertainties regarding the CO2 uptake by the biological carbon pump and its influence on the solubility pump, we propose to include parameters such as preformed nutrient supply, carbon export rates, Redfield ratios, and CO2 concentrations measured at specific key sites into the NIR.