Abstract

In a future ocean, dissolved organic carbon (DOC) release by seaweed has been considered a pathway for organic carbon that is not incorporated into growth under carbon dioxide (CO2) enrichment/ocean acidification (OA). To understand the influence of OA on seaweed DOC release, a 21-day experiment compared the physiological responses of three seaweed species, two which operate CO2 concentrating mechanisms (CCMs), Ecklonia radiata (C. Agardh) J. Agardh and Lenormandia marginata (Hooker F. and Harvey) and one that only uses CO2 (non-CCM), Plocamium cirrhosum (Turner) M.J. Wynne. These two groups (CCM and non-CCM) are predicted to respond differently to OA dependent on their affinities for Ci (defined as CO2 + bicarbonate, HCO3). Future ocean CO2 treatment did not drive changes to seaweed physiology—growth, Ci uptake, DOC production, photosynthesis, respiration, pigments, % tissue carbon, nitrogen, and C:N ratios—for any species, regardless of Ci uptake method. Our results further showed that Ci uptake method did not influence DOC release rates under OA. Our results show no benefit of elevated CO2 concentrations on the physiologies of the three species under OA and suggest that in a future ocean, photosynthetic CO2 fixation rates of these seaweeds will not increase with Ci concentration.

Abbreviations

    Abbreviations
     
  • CCM

    carbon dioxide concentrating mechanism

  •  
  • CO2

    carbon dioxide

  •  
  • Ci

    dissolved inorganic carbon

  •  
  • DOC

    dissolved organic carbon

  •  
  • HCO3-

    bicarbonate

  •  
  • OA

    ocean acidification

Introduction

Over the past 250 years, CO2 concentrations in the atmosphere, driven by anthropogenic emissions, have increased from ~280 to 400 ppmv, of which ∼30% has been absorbed by the surface ocean (Doney et al., 2009; Mcleod et al., 2011; Hurd et al., 2018; IPCC, 2022). The ongoing absorption of atmospheric CO2 by the surface ocean has caused seawater pH to decrease by 0.1 unit since the industrial revolution and altered the relative abundance of the different species of Ci in seawater (termed OA). By the end of this century, it is projected that the pH of the ocean will drop by a further 0.4 units and simultaneously HCO3, the most abundant species of Ci in seawater (concentrations of ∼2000 μM), will increase by 20% and CO2 by 260% (Koch et al., 2013; Stockdale et al., 2016; Young and Gobler, 2016). The changes in concentrations of dissolved HCO3 and CO2, can independently and synergistically influence the physiology of seaweeds (Hurd et al., 2019).

Seaweeds (marine macroalgae) are dominant marine autotrophs in coastal oceans around the world from three distinct evolutionary lineages—the Rhodophyta, Chlorophyta, and Ochrophyta. They provide essential services within the marine environment such as habitat creation, food provision, and nutrient cycling (e.g. CO2, oxygen, and nitrogen) (Steneck et al., 2002; Graham, 2004; Christie et al., 2009; Egan et al., 2013; Roleda and Hurd, 2019). All seaweeds take up CO2 from seawater by passive diffusion, however, the majority of species also operate CO2 concentrating mechanisms (CCMs) to facilitate the use of HCO3, which is abundant in seawater compared to CO2 (Andria et al., 1999; Raven et al., 2008; Fernández et al., 2014). HCO3 is converted to CO2, which is fixed in the chloroplast by RUBISCO and incorporated into amino acids, polysaccharides, and cell walls (Giordano et al., 2005; Raven et al., 2008; Fernández et al., 2014; Hurd et al., 2014). Around 35% of studied seaweeds (mostly phylum Rhodophyta) take up only CO2 and are termed CO2-only users or non-CCMs (Raven et al., 2002; Cornwall et al., 2015; Kubler and Dudgeon, 2015). CCMs have higher energetic requirements than diffusive CO2 uptake and for some species, under CO2 enrichment or low light, CCM activity can be reduced in favour of diffusive CO2, termed “down-regulation” of HCO3 use (Cornwall et al., 2015, 2017; Britton et al., 2019).

Seaweeds are predicted to respond differently to ocean acidification (OA) based on their Ci uptake method and affinity for Ci including: low and high affinity non-CCMs, low and high affinity CCMs with either down regulation or no down regulation of the CCM, and HCO3 only use (Hepburn et al., 2011; Cornwall and Hurd, 2019; Hurd et al., 2019; van der Loos et al., 2019). Various physiological responses of seaweeds to OA have been recorded, ranging from beneficial to negative and neutral (Cornwall et al., 2017; Britton et al., 2019; van der Loos et al., 2019). It is predicted that non-CCM seaweeds with a high affinity for CO2 will not be affected by future ocean CO2 concentrations while non-CCM seaweeds with a low affinity for CO2 may benefit due to the increase in seawater CO2 available for photosynthesis (Hurd et al., 2019; van der Loos et al., 2019). Seaweed CCMs are pH sensitive and observed responses are likely due to differences in CCM effectiveness under pH changes (Ober and Thornber, 2017; van der Loos et al., 2019). When enriched with CO2, some seaweed species show elevated photosynthetic rates but no increase in growth, suggesting an additional pathway for the photosynthetically fixed CO2, which is not used for growth (Olischlager and Wiencke, 2013; Gordillo et al., 2016; Kang and Kim, 2016; Bender-Champ et al., 2017; Reidenbach et al., 2017; Kang and Chung, 2018).

A carbon overflow mechanism (called the carbon overflow hypothesis in older literature), where dissolved organic carbon (DOC) is released into the surrounding seawater to maintain intracellular carbon and nitrogen stoichiometry is one such pathway for fixed CO2 not used in seaweed growth (Gordillo et al., 2001; Iñiguez et al., 2016b; Weigel and Pfister, 2021). Discrepancies between photosynthesis and growth rates of seaweeds under future ocean CO2 concentrations are considered a response of species-specific re-organization of the cellular carbon budget, whereby seaweeds either store the excess carbon or release it into the surrounding seawater as DOC (Iñiguez et al., 2016a). Elevated DOC release rates under OA have been observed for seaweeds on the Great Barrier Reef (Diaz-Pulido and Barron, 2020), in Spain (Gordillo et al., 2001), in the Artic, and Antarctica (Iñiguez et al., 2016a, 2017). This overflow of DOC may be enhanced if seaweeds are limited for growth by other essential nutrients, particularly nitrogen (Gordillo et al., 2001; Riebesell, 2004; Thornton, 2014; Iñiguez et al., 2016a; Diaz-Pulido and Barron, 2020).

The goal of this study was to examine the DOC overflow hypothesis of seaweeds with differing Ci uptake strategies. We tested three seaweed species with known Ci uptake strategies for their physiological responses—growth, Ci uptake, DOC production, photosynthesis, respiration, pigments, % tissue carbon, nitrogen, and C:N ratios—to future ocean CO2 concentrations under a 21-day OA experiment. This is the first study to report DOC release values for a non-CCM seaweed, and only the second to study the effect of OA on the DOC release of temperate seaweeds (Gordillo et al., 2001). The study species were, Ecklonia radiata (Ochrophyta, CCM), Lenormandia marginata (Rhodophyta, CCM), and Plocamium cirrhosum (Rhodophyta, non-CCM)—note the Ci uptake method of each species has been previously confirmed using a pH drift experiment and isotope δ13C analysis by Cornwall et al. (2015). We hypothesized that if the non-CCM seaweed, P. cirrhosum, was limited by CO2 under current ocean concentrations we would observe an increase in CO2 uptake, photosynthesis, and growth in the OA treatment. For the CCM seaweeds, E. radiata and L. marginata, it was hypothesized that there would be no increase in photosynthetic rates or growth but instead an increase in the Ci uptake, in particular CO2 taken up from the seawater under OA—as found by Britton et al. (2016) for E. radiata. We predicted an increase in the DOC release rate for all three seaweed species under OA due to the hypothesized change in Ci uptake and re-organization of the cellular carbon budget. The effect of OA on pigment content, % dry tissue carbon, nitrogen, and C:N ratios for all three species was unknown but considered linked with growth, DOC release, and Ci uptake.

Methods

Seaweed collection

Seaweeds (n = 16 individuals each for E. radiata, P. cirrhosum, and L. marginata) were collected using SCUBA at a depth of 8–10 m from Coal Point, Bruny Island, Tasmania, Australia (43.33 S, 147.32 E) on 12 July 2021. Specimens of each species were pressed using a seaweed press, confirmed using field guide books (Scott, 2013), and lodged at the Tasmanian State Herbarium. Juveniles of E. radiata were used because preliminary studies revealed no difference in the DOC release rates (μmol C gDW−1 h−1) by whole mature blades (blade and stipe) and whole juvenile E. radiata (see Supplementary Figure S1). Mature individuals of P. cirrhosum and L. marginata were collected which were an appropriate size (∼ 50 cm2) for the experimental chambers (650 ml). Seaweeds were placed in an insulated container with seawater for transport to the laboratory 2 h away.

Pre-experimental treatment of seaweeds and equipment cleaning procedure

Once in the laboratory, seaweeds were left in the insulated container overnight with freshly filtered seawater and aeration under 80 μmol photons m−2 s−1 on a 12:12 light/dark cycle in a cool room set to 11°C (similar to conditions at the collection site). Prior to the experiment, seaweeds were cleaned of epiphytes using KimtechTM wipes and tweezers. For L. marginata, calcifying invertebrates could not be removed due to potential tissue damage and covered ∼1% of the seaweed blade on each individual. All glassware, tubes, and syringes used in the experiment were cleaned using Decon-90 (2% v/v) and HCL (10% v/v). Glass DOC vials (Shimadzu TOC certified) and glass filter paper (Whatman GF/F, 0.7 μm) were combusted in a muffler overnight (8 h at 400°C) to remove any residual organic carbon (Yoro et al., 1999). Bottles for total alkalinity samples (HDPE, Nalgene, 60 ml, Thermo ScientificTM) were rinsed five times each in distilled water and left to dry before use.

Experimental culture conditions

On day 1 of the 21-day experiment, 16 individual seaweeds from each species were randomly assigned to one of two pH treatments (“current” and “future,” n = 8 for each species |$\times$| treatment combination ∴ 48 seaweeds total). Mean target pHT was 8.0 for the current treatment and 7.6 for the future treatment. The pHT of both treatments fluctuated on a diurnal cycle to reflect the natural daily pH fluctuation within a seaweed bed, increasing 0.2 units every 4 h during the day and decreasing 0.2 units every 4 h at night (see Supplementary Figure S2; Britton et al., 2019). This gave the treatments a range of 7.85–8.15 pH units (CO2 at 400 ppm) and 7.50–7.80 pH units (CO2 at 850 ppm) for current and future treatments, respectively. The experimental set-up was situated in a cool room set to 11°C (ambient seawater at the time of seaweed collection) and was measured in the chambers throughout the experiment as 13°C ± 0.5 SD. The experiment was conducted in an automated culture system described in Britton et al. (2019). Perspex culture chambers (650 ml, total number of chambers = 48) were placed on stirrer plates underneath full spectrum LED lights (Zeus-70 LED ZEAL) on a diurnal cycle set to a daily maximum of 80 μmol photons m−2 s−1 (12:12, light/dark). This irradiance was typical of the site depth (Britton et al., 2020). Target pHT was achieved using mass flow controllers (FMA5418A and FMA5402A, Omega Engineering, USA) that controlled the ratio of air and CO2 in a gas mix introduced to seawater through mini module membrane contactors (3MTM Liqui-CelTM MM-15.5 Series) from a single header tank using peristaltic pumps. The pH of the incoming seawater was maintained within 0.03 units of the pHT using an automated spectrophotometric feedback system (McGraw et al., 2010). Chambers were filled with seawater corresponding to their designated pH treatment (current or future), which was replaced every 4 h with fresh seawater. The pH of the seawater in each culture chamber was monitored every 2 h using a separate spectrophotometric system to assess how the seaweeds altered the pH of the seawater between seawater changes (for the carbonate chemistry parameters of each treatment see Supplementary Table S1). All pHT measurements were measured between 12.55 and 16.28°C.

Growth rates

To determine growth rates of seaweeds, all individuals were blotted dry using KimWipesTM and weighed for wet weight at the start (day 1), and end (day 21) of the experiment. Growth rates according to wet weight were calculated using equations from Kain (1982).

For E. radiata, a small hole (0.5 cm2) was punched 2.5 cm above the stipe (just above the meristematic region) using a single, steel hole punch for determination of lateral blade elongation (note that blade length and elongation were not measured for P. cirrhosum and L. marginata due to their apical meristematic growth from numerous thalli resulting in this measurement to be unsuitable) (Dean and Hurd, 2007). Photos of seaweeds were taken on day 1 on a light board for initial surface area measurements. At the completion of the 21-day experiment, seaweeds were re-photographed for final determination of surface area, blade elongation, and blade erosion. Images were analysed for surface area using ImageJ (1.52i). Blade length and blade elongation (using a hole punch method) were calculated using the following equation from Dean and Hurd (2007):
where for blade length (BTotL, cm d−1), L1 is the distance from blade/stipe junction to blade tip on day 1, L2 is the distance from blade/stipe junction to blade tip on day 21, and T is the duration of growth period (21 d). For blade elongation (BElgn, cm d−1), L1 is the distance from blade/stipe junction to hole on day 1, L2 is the distance from blade/stipe junction to hole on day 21, and T is the duration of growth period (21 d). Erosion was calculated using the following equation from Dean and Hurd (2007):
where BElgn is the blade elongation as calculated using the hole punch method above (cm d−1) and BTotL is the total blade length (cm d−1) as calculated above.

Immediately after growth measurements were taken, seaweed tissue was preserved for subsequent assays. For pigment content, tissue (0.1 g) was frozen at −80°C and, for % carbon, nitrogen, and C:N ratios, tissue (0.1 g) was oven dried at 60°C for three nights.

Ci uptake

Seawater samples for determination of total alkalinity (AT) were taken from each incubation chamber (48 chambers total) on day 19 of the experiment. Upon refreshment of seawater in each incubation chambers, all chambers were sealed, and an initial sample of seawater (60 ml) was collected in a pre-rinsed alkalinity bottle. Final seawater samples were taken after 4 h, prior to the next seawater refreshment. AT was measured using the open-cell titration method using a Metrohm 862 Compact Titrosampler coupled with an Aquatrode Plus with PT1000 240 temperature sensor following methods outlined in Ferderer et al. (2022) and Dickson et al. (2007). The change in AT measured on day 19 and pH over 4 h on day 17 during the experiment was used to calculate Ci uptake by seaweeds using the CRAN package “seacarb” in the software program R version 4.1.3 (R Core Team, 2022). An issue with the spectrophotometric system prevented use of pH measurements taken at the same time as AT. Ci uptake was standardized with seaweed wet weight (gWW−1) and time (h).

DOC seawater sample analysis

To estimate the rate of DOC release from each replicate seaweed (n = 8 for each species |$\times$| treatment combination ∴ 48 seaweeds total) over a 24 h period, seawater samples were taken at 8 a.m. (day 20) immediately after each chamber flush and final samples were taken at 8 a.m. on day 21—chambers were not flushed over this incubation time (note that pH would have fluctuated in each chamber over this period with metabolic activity. The measured pH fluctuation was ± 0.01–0.08 over 2 h). DOC samples were analysed at the Queensland Chemistry Centre using an automated total organic carbon analyser (Analytica Jena Multi N/C 3100) via combustion at 720°C over a platinum catalyst in accordance with method 5310 D (Standard Methods Committee of the American Public Health Association 2017). DOC release rates were determined as the difference between initial and final samples and were standardized to incubation time (h), culture volume (L), and seaweed biomass (grams of dry weight, gDW). The rate of DOC production was expressed as μmol C gDW−1 h−1.

Photosynthetic and respiration rates

On the final day of the experiment (day 21), the change in oxygen concentration was measured in each culture chamber (48 total); initial and final dissolved oxygen concentrations (mmol O2 g−1 h−1) were measured under experimental irradiance over 3 h for net photosynthesis. The lights were turned off and algae incubated for a further 3 h for dark respiration. Dissolved oxygen concentration was measured using a portable oxygen meter (HachVR, HQ40D Portable Multi Meter) and net photosynthesis/respiration were standardized to fresh weight (gWW) and expressed as μmol O2 g−1 h−1.

Seawater nutrient analysis

To assess nutrient concentrations, seawater samples were taken from both the header tank and a random subset of chambers on the last day of the experiment (day 21). Seawater samples (12 ml) were taken using pre-cleaned syringes, filtered (Whatman GF/F 0.7 μM) and stored in 12 ml polyethylene nutrient tubes at −20°C for later analysis. A QuickChem® 8000 Automated Ion Analyser (LaChat Instruments) was used to determine concentrations of nitrate (NO3, measured as NO3 + nitrite [NO2]) and ammonium (NH4+) following methods outlined in Diamond (2008) and Liao (2008).

Percent dry tissue carbon, nitrogen, C:N ratios, and δ13C isotopes

Oven dry seaweed tissue (n = 8 for each species |$\times$| treatment combination ∴ 48 seaweeds total) was analysed for % tissue nitrogen, % tissue carbon, and C:N ratios (calculated as mol:mol). To confirm and study changes in the Ci uptake method of each seaweed species, stable carbon isotopes (δ13C) were determined using a NA1500 elemental analyser coupled to a Thermo Scientific Delta V Plus via a Conflo IV. Tissue was combusted at 1020°C and reduced at 650°C. Carbon stable isotope ratios were normalized using NBS-22, NBS-19, and LSVEC.

Pigment content

For the red seaweeds, L. marginata and P. cirrhosum, phycobilin and total chlorophyll contents were determined using preserved algal material (0.1 ± 0.01 g) following methods outlined in Pritchard et al. (2013). Under dim light, seaweed tissue was ground in a mortar and pestle using liquid nitrogen, transferred to polyethylene centrifuge tubes (15 ml) with 5 ml of ice-cold 0.1 M phosphate buffer (pH 6.8), and left in a fridge at 4°C for overnight extraction of phycobiliproteins. The following day, the solution was centrifuged and the supernatant decantated off and absorbance measured at 564, 592, 455, 618, and 645 nm using a scanning spectrophotometer (Halo RB-10, Dynamica). Phycobilin content were calculated using equations from Beer and Eshel (1985). Chlorophyll content was extracted from the remaining pellets with 5 ml of 100% ethanol. Solution was vortexed and left overnight at 4°C. Solution was again centrifuged and the supernatant absorbance was measured at 629, 649, and 665 nm. Chlorophyll content were calculated using equations from Ritchie (2006).

For the brown seaweed E. radiata, fucoxanthin and chlorophyll contents were extracted from preserved tissue (0.1 ± 0.01 g) using methods outlined in Seely et al. (1972). Frozen tissue was placed in centrifuge tubes (15 ml) with 4 ml of dimethyl-sulfoxide (DMSO) and left to extract for 10 min. Decanted DMSO extract was measured for absorbance using a scanning spectrophotometer (Halo RB-10, Dynamica) at wavelengths of 665, 631, 582, and 480 nm. After removal of DMSO, 6 ml of ice-cold acetone (90% v/v) was added to the tubes with seaweed tissue and left to extract for 30 min. Tubes were agitated sporadically throughout this extraction time and visual observation ensured no remaining pigments in the seaweed tissue prior to decanting off acetone solution. Acetone solution was then measured for absorbance at wavelengths 664, 631, 581, and 470 nm. Fucoxanthin and chlorophyll content were calculated using equations described in Seely et al. (1972).

Statistical analysis

Statistical analyses were performed in the software program R version 4.1.3 (R Core Team, 2022) and graphs were created using SigmaPlot (Systat Software Inc). One-way ANOVAs were used to determine the effect of CO2 treatment on E. radiata blade elongation and erosion, as well as pigment content for each species individually (due to different pigments). To compare the parameters for growth, Ci uptake, DOC production, photosynthesis, respiration, % tissue carbon, nitrogen, and C:N ratios, two-way ANOVAs were conducted with the fixed factors Treatment (2 levels, current, and future), Species (3 levels, E. radiata,L. marginata, andP. cirrhosum) and the Treatment|$\times $|Species interaction. ANOVA assumptions of normality of residuals and homoscedasticity were assessed using diagnostic plots of model residuals and data were transformed where necessary using the package, bestNormalize (Peterson, 2021); see Supplementary Table S2 for data transformations. Where Species was statistically significant and/or an interaction was observed between Treatment|$\times $|Species (p < 0.05), a Tukey’s HSD multiple comparison test was run to distinguish which means were statistically significant from others and are discussed in the results section.

Results

Growth rates

There was no difference between CO2 treatments for growth rates (ANOVA: F(1,39) = 0.07, p = 0.793; Figure 1a) of all three species studied and no interaction between Treatment|$\times $|Species. Species influenced growth with E. radiata growth rates being higher than L. marginata and P. cirrhosum (ANOVA: F(2,39) = 94.967, p < 0.001) but were highly variable under both current (0.022–0.054 gWW−1 d−1) and future CO2 concentrations (0.018–0.053 gWW−1 d−1). L. marginata showed negligible and negative growth rates for both CO2 concentrations, and P. cirrhosum exhibited no to little growth (-0.001–0.019 gWW−1 d−1) under both CO2 concentrations (Figure 1a). Similarly, no difference was observed between current and future treatments for E. radiata blade elongation and distal erosion (elongation ANOVA: F(1,13) = 0.461, p = 0.509 and erosion ANOVA: F(1,13) = 1.145, p = 0.304; Figure 1b).

(a) Growth rates (gWW−1 d−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, dark grey bars) and future (850 ppm, pH 7.50–7.80, white bars) CO2 concentrations, and (b) blade elongation and erosion (cm d−1) of E. radiata, measured using the hold punch method, under current and future fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way and one-way ANOVAs and post hoc test results see Supplementary Tables S2 and S3.
Figure 1.

(a) Growth rates (gWW−1 d−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, dark grey bars) and future (850 ppm, pH 7.50–7.80, white bars) CO2 concentrations, and (b) blade elongation and erosion (cm d−1) of E. radiata, measured using the hold punch method, under current and future fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way and one-way ANOVAs and post hoc test results see Supplementary Tables S2 and S3.

Ci uptake and carbon stable isotopes

CO2 treatment had no effect on Ci uptake (CO2 + HCO3 combined) (ANOVA: F(1, 33) = 0.216, p = 0.645; Figure 2a). Note that air bubbles were inevitably created in each sealed chamber due to the initial sampling of the seawater and could have slightly reduced the Ci in the future treatment due to outgassing however this was considered minimal given the short incubation period and small headspace to volume ratio. Species did influence Ci uptake, with P. cirrhosum taking up more Ci than L. marginata (ANOVA: F(2, 33) = 3.698, p = 0.036; Figure 2a).

(a) Ci uptake (CO2 + HCO3−) and (b) δ13C isotope values for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white bars) and future (850 ppm, pH 7.50–7.80, grey bars) fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.
Figure 2.

(a) Ci uptake (CO2 + HCO3) and (b) δ13C isotope values for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white bars) and future (850 ppm, pH 7.50–7.80, grey bars) fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.

To determine the Ci uptake method, carbon stable isotopes were analysed with CO2 treatment having no effect on carbon stable isotopes (δ13C) (ANOVA: F(1, 39) = 0.460, p = 0.502; Figure 2b). There was however an effect of species with P. cirrhosum having more negative ‰ values than E. radiata and L. marginata (ANOVA: F(2, 39) = 461.579, p < 0.001; Figure 2b).

DOC production

There was no difference between future and current ocean CO2 concentrations for DOC release rates (ANOVA: F(1, 37) = 3.664, p > 0.0633; Figure 3). DOC release was highly variable for both treatments for all three species studied with no influence of species on DOC production (ANOVA: F(2, 37) = 0.894, p = 0.418; Figure 3) or the interaction between Treatment|$\times $|Species.

Net DOC production (μmol C gDW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white circles) and future (850 ppm, pH 7.50–7.80, grey circles) fluctuating CO2 concentrations. Circles are replicate seaweeds (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.
Figure 3.

Net DOC production (μmol C gDW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white circles) and future (850 ppm, pH 7.50–7.80, grey circles) fluctuating CO2 concentrations. Circles are replicate seaweeds (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.

Photosynthesis and respiration

CO2 treatment had no effect on photosynthetic rates (ANOVA: F(1, 39) = 0.011, p = 0.916; Figure 4a) but species did (ANOVA: F(2, 39) = 7.503, p = 0.001; Figure 4a) with L. marginata exhibiting the highest rates of oxygen production (μmol O2 gWW−1 h−1) of the three species studied. Respiration rates were not influenced by CO2 treatment (ANOVA: F(1, 39) = 0.249, p = 0.668; Figure 4b), species (ANOVA: F(2, 39) = 0.710, p = 0.585; Figure 4b), or the interaction between Treatment|$\times $|Species.

(a) Net photosynthetic rates (μmol O2 gWW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, dark grey bars) and future (850 ppm, pH 7.50–7.80, white bars) CO2 concentrations, and (b) respiration rates (μmol O2 gWW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current and future fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.
Figure 4.

(a) Net photosynthetic rates (μmol O2 gWW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, dark grey bars) and future (850 ppm, pH 7.50–7.80, white bars) CO2 concentrations, and (b) respiration rates (μmol O2 gWW−1 h−1) for E. radiata, L. marginata, and P. cirrhosum under current and future fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). For two-way ANOVA and post hoc test results see Supplementary Tables S2 and S3.

Percent tissue carbon and nitrogen and pigment content

CO2 treatment had no effect on % dry tissue carbon (ANOVA: F(1, 39) = 0.299, p = 0.588; Figure 5a), % nitrogen content (ANOVA: F(1, 39) = 25.2, p = 0.0937; Figure 5b), or C:N (ANOVA: F(1, 39) = 0.921, p = 0.343; Figure 5c). For the brown seaweed E. radiata, there were higher pigment contents of both chlorophyll a and fucoxanthin under current compared to future CO2 treatment (chlorophyll ANOVA: F(1, 11) = 6.212, p = 0.029 and fucoxanthin ANOVA: F(1, 11) = 4.433, p = 0.059; Figure 5d). CO2 treatment had no influence on L. marginata chlorophyll and phycobilin content (chlorophyll ANOVA: F(1, 11) = 0.571, p = 0.466 and phycobilin ANOVA: F(1, 11) = 0.547, p = 0.475; Figure 5d) or P. cirrhosum chlorophyll and phycobilin content (chlorophyll ANOVA: F(1, 11) = 0.172, p = 0.686 and phycobilin ANOVA: F(1, 11) = 1.013, p = 0.336; Figure 5d).

(a) Percent tissue carbon, (b) % tissue nitrogen, (c) C:N ratio, and (d) pigment content for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white bars) and future (850 ppm, pH 7.50–7.80, grey bars) fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). Significant differences between current and future CO2 indicated by asterisks (*). For one way and two-way ANOVAs and post hoc test results see Supplementary Tables S2 and S3.
Figure 5.

(a) Percent tissue carbon, (b) % tissue nitrogen, (c) C:N ratio, and (d) pigment content for E. radiata, L. marginata, and P. cirrhosum under current (400 ppm, pH 7.85–8.15, white bars) and future (850 ppm, pH 7.50–7.80, grey bars) fluctuating CO2 concentrations. Bars are mean ± SE (n = 8). Significant differences between current and future CO2 indicated by asterisks (*). For one way and two-way ANOVAs and post hoc test results see Supplementary Tables S2 and S3.

Discussion

The physiologies of E. radiata, L. marginata, and P. cirrhosum, did not respond to future compared to current ocean CO2 treatment, regardless of Ci uptake method (CCM vs. non-CCM). Although negligible growth rates were observed for both red seaweeds, photosynthetic rates were positive suggesting that photosynthate (organic carbon created via photosynthesis) was used for other cellular processes—or released as DOC—rather than invested in new tissue growth (Hurd et al., 2014). The isotope values indicated that L. marginata and E. radiata were both using CCMs to take up HCO3 during the experiment (δ13C between −10 and −30‰) and P. cirrhosum only used CO213C more negative than −30‰). The only statistical significance for treatment was for E. radiata chlorophyll and fucoxanthin contents which were reduced under OA, likely due to re-arrangement of the light harvesting system at higher CO2 concentrations by which excess pigment molecules are metabolized, termed the “pigment economy” hypothesis (Gordillo et al., 1999, 2015).

Various physiological responses of seaweeds to OA have been recorded (see Table 2 in van der Loos et al. [2019]), and our results show that for E. radiata, L. marginata, and P. cirrhosum there was no benefit of OA on seaweed growth and photosynthesis, i.e., a neutral response (Cornwall et al., 2017; Britton et al., 2019). The neutral response to OA observed in this study may have been a physiological mechanism to counteract the increase of H+ caused by elevated pCO2, which could have driven homeostatic functions such as the synthesis of H+ transport-related proteins, rather than invest metabolic energy into growth (Roleda et al., 2012; Xu et al., 2017). Combined, these findings suggest E. radiata, L. marginata, and P. cirrhosum were not taking up, or metabolizing, more Ci in the future compared to the current ocean treatment.

Relationships between DOC release and OA

These results are just the second published measurements of the effect of OA on DOC release rates by temperate (cool water) seaweeds (the first being, Gordillo et al., 2001). All three species in our study released DOC at similar rates suggesting that Ci uptake strategy did not influence the rate of DOC release—a novel finding, which has not been previously reported in the literature. The rate of DOC release for E. radiata, L. marginata, and P. cirrhosum was not influenced by CO2 treatment, in contrast to the results of similar studies on the Great Barrier Reef, Australia (Diaz-Pulido and Barron, 2020), Spain (Gordillo et al., 2001), and in the Arctic and Antarctica (Iñiguez et al., 2016a, 2017), which all found elevated DOC release under future ocean CO2 concentrations. These previous studies suggest that DOC was released as an overflow under CO2 enrichment to prevent a build-up of photosynthetic products, which were not used for growth, thereby maintaining intracellular C:N balance when environmental conditions were unfavourable (such as nitrogen limitation, discussed below) (Gordillo et al., 2001; Iñiguez et al., 2016a, 2017; Diaz-Pulido and Barron, 2020). In our study, no increase in DOC release under OA was, therefore, likely a result of the non-response of growth, photosynthesis, and Ci uptake for all three seaweeds, causing the intracellular C:N balance to similarly remain unaffected by the OA treatment.

Interactions between OA and nitrogen

In a previous study, the combination of nitrogen enrichment and OA increased growth rates of Ulva rigida in Spain, suggesting that the physiological benefits of OA were dependent on the enhancement effect of CO2 on nitrogen assimilation (Gordillo et al., 2001). In relation to this, DOC released as a % of the primary production was highest under nitrogen limitation and CO2 enrichment suggesting release as a mechanism of overflow during unfavourable environmental conditions (Gordillo et al., 2001). Ambient NO3 concentrations in Tasmania (maintained throughout our experiment with chamber flushing) were 1000-fold lower than used by Gordillo et al., (2001) in Spain (0.00043 mM compared to 0.25–5 mM). The non-response to future ocean CO2 concentrations, therefore, may be explained by the closely linked metabolism of carbon and nitrogen (Turpin et al., 1990; Raven et al., 2004; Roleda and Hurd, 2019). Both red seaweeds were nitrogen sufficient C:N ∼10 mol:mol, however E. radiata, showed nitrogen limitation with C:N ratios ∼20 mol:mol, suggesting that DOC is not released under nitrogen limitation for this species. In contrast, in a parallel study we found that E. radiata, L. marginata, and P. cirrhosum released DOC when intracellular C:N values were lowest annually at ≤10 mol:mol for L. marginata and P. cirrhosum and ∼15 mol:mol for E. radiata, indicating nitrogen sufficiency (Paine et al., in Review). The overflow of DOC from these three seaweed species is, therefore, likely caused by nitrogen sufficiency preventing the “shutting off” of photosynthesis, causing DOC to be released as an overflow to prevent a detrimental build-up of photosynthetic products around RUBISCO (Thornton, 2014; Paine et al., 2021).

Influence of surface epibionts, and the microbiome on responses to OA

There was high variability in DOC release rates with both positive and negative production rates. Although not tested here, negative release rates were likely a result of DOC consumption by the associated seaweed macro- and micro-biome (Haas et al., 2011; Nelson et al., 2013). Lenormandia marginata is a seaweed, which is highly colonized by animal epibionts, usually a combination of spiroid worms (covering ∼1% of the seaweed blade) and bryozoans, which likely influenced the species response to current and future CO2 concentrations. Although seaweeds were cleaned using tissue, the calcifying invertebrates could not be removed because of the risk of tissue damage. Although not tested, the high variability in L. marginata growth, photosynthesis, and respiration results (in particular for the future compared to current CO2 treatment) may be attributed to respiration, or lysis, of the associated surface epibionts under OA. Additionally, there was a high standard deviation for DOC release by L. marginata (current SD ± 14.41 and future SD ± 24.82) with numerous replicates showing evidence of DOC uptake (as observed by negative DOC release rates) potentially occurring as a result of consumption by the macro- and micro-biome (Haas et al., 2011; Nelson et al., 2013). This relationship between seaweed and its associated biome is an area of particular research interest at the moment due to uncertainty in how this relationship is beneficial to the seaweed growth in the marine environment (Egan et al., 2013).

Conclusions

In conclusion, growth, photosynthetic rates, and Ci uptake were not increased in the future compared to current CO2 concentration treatment for either, E. radiata, L. marginata, or P. cirrhosum, regardless of Ci uptake strategy. This result may suggest all species had high affinities for Ci but further research is required to understand the putative responses of Ci uptake mechanisms to future ocean change. OA/CO2 enrichment had no effect on seaweed DOC release likely because overflow of photosynthate was not required as both photosynthesis and growth were similarly unaffected by OA. Our results show no benefit of elevated CO2 concentrations for any of the study seaweeds in a future ocean.

ACKNOWLEDGEMENTS

We thank Pam Quayle, Eloise Bennett, and Olivia Wynn for their help in the laboratory and field. We acknowledge and pay respects to the Palawa people of Lunawanna-allonah and Nipaluna, whose water and land upon which we work.

Funding

This study was funded by the Australian Research Council by Laureate Fellowship number FL160100131 (to PWB), discovery project number DP160103071 funding to GD-P and CLH, and DP200101467 to CLH. DB was funded by the Tasmanian Department of Natural Resources and Environment (AIRF: 2020–047).

Data availability

Available upon request.

Authors’ contributions

ERP: experimental design (lead), conducted experiments (lead), data analysis (lead), and drafting and editing of manuscript (lead); DB: experimental design and laboratory support (supporting), data analysis (supporting), and editing of manuscript (supporting); MS: editing of manuscript (supporting); EAB: laboratory analysis support and editing of manuscript; GD-P: original concept (supporting) and editing of manuscript (supporting); PWB: editing of manuscript (supporting); CLH: original concept, experimental design (supporting), and editing manuscript (supporting)

References

Andria
J. R.
,
Vergara
J. J.
,
Perez-Llorens
J. L.
1999
.
Biochemical responses and photosynthetic performance of Gracilaria sp. (Rhodophyta) from Cadiz, Spain, cultured under different inorganic carbon and nitrogen levels
.
European Journal of Phycology
,
34
:
497
504
.

Beer
S.
,
Eshel
A.
1985
.
Determining phycoerythrin and phycocyanin concentrations in aqueous crude extracts of red algae
.
Marine and Freshwater Research
,
36
:
785
792
.

Bender-Champ
D.
,
Diaz-Pulido
G.
,
Dove
S.
2017
.
Effects of elevated nutrients and CO2 emission scenarios on three coral reef macroalgae
.
Harmful Algae
,
65
:
40
51
.

Britton
D.
,
Cornwall
C. E.
,
Revill
A. T.
,
Hurd
C. L.
,
Johnson
C. R.
2016
.
Ocean acidification reverses the positive effects of seawater pH fluctuations on growth and photosynthesis of the habitat-forming kelp, Ecklonia radiata
.
Scientific Reports
,
6
:
26036
.

Britton
D.
,
Mundy
C. N.
,
McGraw
C. M.
,
Revill
A. T.
,
Hurd
C. L.
,
Norkko
J.
2019
.
Responses of seaweeds that use CO2 as their sole inorganic carbon source to ocean acidification: differential effects of fluctuating pH but little benefit of CO2 enrichment
.
ICES Journal of Marine Science
,
76
:
1860
1870
.

Britton
D.
,
Schmid
M.
,
Noisette
F.
,
Havenhand
J. N.
,
Paine
E. R.
,
McGraw
C. M.
,
Revill
A. T.
et al.
2020
.
Adjustments in fatty acid composition is a mechanism that can explain resilience to marine heatwaves and future ocean conditions in the habitat-forming seaweed Phyllospora comosa (Labillardière) C.Agardh
.
Global Change Biology
,
26
:
3512
3524
.

Christie
H.
,
Norderhaug
K. M.
,
Fredriksen
S.
2009
.
Macrophytes as habitat for fauna
.
Marine Ecology Progress Series
,
396
:
221
233
.

Cornwall
C. E.
,
Hurd
C. L.
2019
.
Variability in the benefits of ocean acidification to photosynthetic rates of macroalgae without CO2-concentrating mechanisms
.
Marine and Freshwater Research
,
71
:
275
280
.

Cornwall
C. E.
,
Revill
A. T.
,
Hall-Spencer
J. M.
,
Milazzo
M.
,
Raven
J. A.
,
Hurd
C. L.
2017
.
Inorganic carbon physiology underpins macroalgal responses to elevated CO2
.
Scientific Reports
,
7
:
46297
.

Cornwall
C. E.
,
Revill
A. T.
,
Hurd
C. L.
2015
.
High prevalence of diffusive uptake of CO2 by macroalgae in a temperate subtidal ecosystem
.
Photosynthesis Research
,
124
:
181
190
.

Dean
P. R.
,
Hurd
C. L.
2007
.
Seasonal growth, erosion rates, and nitrogen and photosynthetic ecophysiology of Undaria pinnatifida (Heterokontophyta) in southern New Zealand
.
Journal of Phycology
,
43
:
1138
1148
.

Diamond
D.
2008
.
Determination of nitrate/nitrite in brackish or seawater by flow injection analysis
.
ICES Document 31-107-04-1-A
.
19
pp.

Diaz-Pulido
G.
,
Barron
C.
2020
.
CO2 enrichment stimulates dissolved organic carbon release in coral reef macroalgae
.
Journal of Phycology
,
56
:
1039
1052
.

Dickson
A.
,
Chris
S.
,
Christian
J. R.
2007
.
Guide to best practices for ocean CO2 measurements
.
North Pacific Marine Science Organization
,
Sidney
.
3
pp.

Doney
S. C.
,
Fabry
V. J.
,
Feely
R. A.
,
Kleypas
J. A.
2009
.
Ocean acidification: the other CO2 problem
.
Annual Review of Marine Science
,
1
:
169
192
.

Egan
S.
,
Harder
T.
,
Burke
C.
,
Steinberg
P.
,
Kjelleberg
S.
,
Thomas
T
.
2013
.
The seaweed holobiont: understanding seaweed–bacteria interactions
.
FEMS Microbiology Reviews
,
37
:
462
476
.

Ferderer
A.
,
Chase
Z.
,
Kennedy
F.
,
Schulz
K.
,
Bach
L.
2022
;
Assessing the influence of ocean alkalinity enhancement on a coastal phytoplankton community
.
Biogeosciences
,
19
:
5375
5399
.

Fernández
P. A.
,
Hurd
C. L.
,
Roleda
M. Y.
,
Wernberg
T.
2014
.
Bicarbonate uptake via an anion exchange protein is the main mechanism of inorganic carbon acquisition by the giant kelp Macrocystis pyrifera (Laminariales, Phaeophyceae) under variable pH
.
Journal of Phycology
,
50
:
998
1008
.

Giordano
M.
,
Beardall
J.
,
Raven
J. A
.
2005
.
CO2 concentrating mechanisms in algae: mechanisms, environmental modulation, and evolution
.
Annual Review of Plant Biology
,
56
:
99
131
.

Gordillo
F. J. L.
,
Aguilera
J.
,
Wiencke
C.
,
Jiménez
C.
2015
.
Ocean acidification modulates the response of two Arctic kelps to ultraviolet radiation
.
Journal of Plant Physiology
,
173
:
41
50
.

Gordillo
F. J. L.
,
Carmona
R.
,
Viñegla
B.
,
Wiencke
C.
,
Jiménez
C.
2016
.
Effects of simultaneous increase in temperature and ocean acidification on biochemical composition and photosynthetic performance of common macroalgae from Kongsfjorden (Svalbard)
.
Polar Biology
,
39
:
1993
2007
.

Gordillo
F. J. L.
,
Jiménez
C.
,
Figueroa
F. L.
,
Niell
F. X.
1998
.
Effects of increased atmospheric CO2 and N supply on photosynthesis, growth and cell composition of the cyanobacterium Spirulina platensis (Arthrospira)
.
Journal of Applied Phycology
,
10
:
461
.

Gordillo
F. J. L.
,
Niell
F. X.
,
Figueroa
F. L.
2001
.
Non-photosynthetic enhancement of growth by high CO2 level in the nitrophilic seaweed Ulva rigida C. Agardh (Chlorophyta)
.
Planta
,
213
:
64
70
.

Graham
M. H.
2004
.
Effects of local deforestation on the diversity and structure of southern California giant kelp forest food webs
.
Ecosystems
,
7
:
341
357
.

Haas
A. F.
,
Nelson
C. E.
,
Kelly
L. W.
,
Carlson
C. A.
,
Rohwer
F.
,
Leichter
J. J.
,
Wyatt
A.
et al.
2011
.
Effects of coral reef benthic primary producers on dissolved organic carbon and microbial activity
.
PLoS ONE
,
6
:
e27973
.

Hepburn
C. D.
,
Pritchard
D. W.
,
Cornwall
C. E.
,
McLeod
R. J.
,
Beardall
J.
,
Raven
J. A.
,
Hurd
C. L.
2011
.
Diversity of carbon use strategies in a kelp forest community: implications for a high CO2 ocean
.
Global Change Biology
,
17
:
2488
2497
.

Hurd
C. L.
,
Beardall
J.
,
Comeau
S.
,
Cornwall
C. E.
,
Havenhand
J. N.
,
Munday
P. L.
,
Parker
L. M.
et al.
2019
.
Ocean acidification as a multiple driver: how interactions between changing seawater carbonate parameters affect marine life
.
Marine and Freshwater Research
,
71
:
263
274
.

Hurd
C. L.
,
Harrison
P. J.
,
Bischof
K.
,
Lobban
C. S.
2014
.
Seaweed Ecology and Physiology
,
Cambridge University Press
,
Cambridge
.
551
pp.

Hurd
C. L.
,
Lenton
A.
,
Tilbrook
B.
,
Boyd
P. W.
2018
.
Current understanding and challenges for oceans in a higher-CO2 world
.
Nature Climate Change
,
8
:
686
694
.

Iñiguez
C.
,
Carmona
R.
,
Lorenzo
M. R.
,
Niell
F. X.
,
Wiencke
C.
,
Gordillo
F. J. L.
2016a
.
Increased CO2 modifies the carbon balance and the photosynthetic yield of two common Arctic brown seaweeds: D esmarestia aculeata and Alaria esculenta
.
Polar Biology
,
39
:
1979
1991
.

Iñiguez
C.
,
Carmona
R.
,
Lorenzo
M. R.
,
Niell
F. X.
,
Wiencke
C.
,
Gordillo
F. J. L.
2016b
.
Increased temperature, rather than elevated CO2, modulates the carbon assimilation of the Arctic kelps Saccharina latissima and Laminaria solidungula
.
Marine Biology
,
163
:
248
.

Iñiguez
C.
,
Heinrich
S.
,
Harms
L.
,
Gordillo
F. J. L.
2017
.
Increased temperature and CO2 alleviate photoinhibition in Desmarestia anceps: from transcriptomics to carbon utilization
.
Journal of Experimental Botany
,
68
:
3971
3984
.

IPCC
.
2022
.
Climate Change 2022: mitigation of climate change
.
Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change
.

Kain
J. M.
1982
.
Morphology and growth of the giant kelp Macrocystis pyrifera in New Zealand and California
.
Marine Biology
,
67
:
143
157
.

Kang
E. J.
,
Kim
K. Y.
2016
.
Effects of future climate conditions on photosynthesis and biochemical component of Ulva pertusa (Chlorophyta)
.
Algae
,
31
:
49
59
.

Kang
J. W.
,
Chung
I. K
.
2018
.
The interactive effects of elevated CO2 and ammonium enrichment on the physiological performances of Saccharina japonica (Laminariales, Phaeophyta)
.
Ocean Science Journal
,
53
:
487
497
.

Koch
M.
,
Bowes
G.
,
Ross
C.
,
Zhang
X. H.
2013
.
Climate change and ocean acidification effects on seagrasses and marine macroalgae
.
Global Change Biology
,
19
:
103
132
.

Kubler
J. E.
,
Dudgeon
S. R.
2015
.
Predicting effects of ocean acidification and warming on algae lacking carbon concentrating mechanisms
.
PLoS ONE
,
10
:
e0132806
.

Liao
N.
2008
.
Determination of ammonia in brackish or seawater by flow injection analysis
.
ICES Document 31-107-06-1-B
.
20
pp.

McGraw
C. M.
,
Cornwall
C. E.
,
Reid
M. R.
,
Currie
K. I.
,
Hepburn
C. D.
,
Boyd
P.
,
Hurd
C. L.
et al.
2010
.
An automated pH-controlled culture system for laboratory-based ocean acidification experiments
.
Limnology and Oceanography: Methods
,
8
:
686
694
.

Mcleod
E.
,
Chmura
G. L.
,
Bouillon
S.
,
Salm
R.
,
Björk
M.
,
Duarte
C. M.
,
Lovelock
C. E.
et al.
2011
.
A blueprint for blue carbon: toward an improved understanding of the role of vegetated coastal habitats in sequestering CO2
.
Frontiers in Ecology and the Environment
,
9
:
552
560
.

Nelson
C. E.
,
Goldberg
S. J.
,
Kelly
L. W.
,
Haas
A. F.
,
Smith
J. E.
,
Rohwer
F.
,
Carlson
C. A
.
2013
.
Coral and macroalgal exudates vary in neutral sugar composition and differentially enrich reef bacterioplankton lineages
.
The ISME Journal
,
7
:
962
979
.

Ober
G. T.
,
Thornber
C. S.
2017
.
Divergent responses in growth and nutritional quality of coastal macroalgae to the combination of increased pCO2 and nutrients
.
Marine Environmental Research
,
131
:
69
79
.

Olischlager
M.
,
Wiencke
C.
2013
.
Ocean acidification alleviates low-temperature effects on growth and photosynthesis of the red alga Neosiphonia harveyi (Rhodophyta)
.
Journal of Experimental Botany
,
64
:
5587
5597
.

Paine
E. R.
,
Schmid
M.
,
Boyd
P. W.
,
Diaz-Pulido
G.
,
Hurd
C. L.
2021
.
Rate and fate of dissolved organic carbon release by seaweeds: a missing link in the coastal ocean carbon cycle
.
Journal of Phycology
,
57
:
1375
1391
.

Peterson
R. A.
2021
.
Finding optimal normalizing transformations via bestNormalize
.
The R Journal
,
13
:
310
329
.

Pritchard
D. W.
,
Hurd
C. L.
,
Beardall
J.
,
Hepburn
C. D
.
2013
.
Survival in low light: photosynthesis and growth of a red alga in relation to measured in situ irradiance
.
Journal of Phycology
,
49
:
867
879
.

R Core Team
.
2022
.
R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria
. .

Raven
J. A.
,
Giordano
M.
,
Beardall
J.
2008
.
Insights into the evolution of CCMs from comparisons with other resource acquisition and assimilation processes
.
Physiologia Plantarum
,
133
:
4
14
.

Raven
J. A.
,
Handley
L. L.
,
Andrews
M.
2004
.
Global aspects of C/N interactions determining plant-environment interactions
.
Journal of Experimental Botany
,
55
:
11
25
.

Raven
J. A.
,
Johnston
A. M.
,
Kübler
J. E.
,
Korb
R.
,
McInroy
S. G.
,
Handley
L. L.
,
Scrimgeour
C. M.
et al.
2002
.
Mechanistic interpretation of carbon isotope discrimination by marine macroalgae and seagrasses
.
Functional Plant Biology
,
29
:
355
378
.

Reidenbach
L. B.
,
Fernandez
P. A.
,
Leal
P. P.
,
Noisette
F.
,
McGraw
C. M.
,
Revill
A. T.
,
Hurd
C. L.
et al.
2017
.
Growth, ammonium metabolism, and photosynthetic properties of Ulva australis (Chlorophyta) under decreasing pH and ammonium enrichment
.
PLoS ONE
,
12
:
e0188389
.

Riebesell
U.
2004
.
Effects of CO2 enrichment on marine phytoplankton
.
Journal of Oceanography
,
60
:
719
729
.

Ritchie
R. J.
2006
.
Consistent sets of spectrophotometric chlorophyll equations for acetone, methanol and ethanol solvents
.
Photosynthesis Research
,
89
:
27
41
.

Roleda
M. Y.
,
Hurd
C. L.
2019
.
Seaweed nutrient physiology: application of concepts to aquaculture and bioremediation
.
Phycologia
,
58
:
552
562
.

Roleda
M. Y.
,
Morris
J. N.
,
McGraw
C. M.
,
Hurd
C. L.
2012
.
Ocean acidification and seaweed reproduction: increased CO2 ameliorates the negative effect of lowered pH on meiospore germination in the giant kelp Macrocystis pyrifera (Laminariales, Phaeophyceae)
.
Global Change Biology
,
18
:
854
864
.

Scott
F. J.
2013
.
Rare marine macroalgae of Southern Australia
. In
Botany
.
University of Tasmania
,
Tasmania
.

Seely
G. R.
,
Duncan
M. J.
,
Vidaver
W. E.
1972
.
Preparative and analytical extraction of pigments from brown algae with dimethyl sulfoxide
.
Marine Biology
,
12
:
184
188
.

Steneck
R. S.
,
Graham
M. H.
,
Bourque
B. J.
,
Corbett
D.
,
Erlandson
J. M.
,
Estes
J. A.
,
Tegner
M. J.
2002
.
Kelp forest ecosystems: biodiversity, stability, resilience and future
.
Environmental Conservation
,
29
:
436
459
.

Stockdale
A.
,
Tipping
E.
,
Lofts
S.
,
Mortimer
R. J. G.
2016
.
Effect of ocean acidification on organic and inorganic speciation of trace metals
.
Environmental Science and Technology
,
50
:
1906
1913
.

Thornton
D. C. O.
2014
.
Dissolved organic matter (DOM) release by phytoplankton in the contemporary and future ocean
.
European Journal of Phycology
,
49
:
20
46
.

Turpin
D. H.
,
Weger
H. G.
,
Smith
R. G.
,
Plaxton
W. C.
,
Lin
M
.
1990
.
Interactions between respiration and nitrogen assimilation during photosynthesis
. In
Inorganic Nitrogen in Plants and Microorganisms: Uptake and Metabolism
, pp.
124
130
.. Ed. by
Ullrich
W. R.
,
Aparicio
P. J.
,
Fuggi
A.
,
Rigano
C.
Springer
,
Berlin, Heidelberg
.

van der Loos
L. M.
,
Schmid
M.
,
Leal
P. P.
,
McGraw
C. M.
,
Britton
D.
,
Revill
A. T.
,
Virtue
P.
et al.
2019
.
Responses of macroalgae to CO2 enrichment cannot be inferred solely from their inorganic carbon uptake strategy
.
Ecology and Evolution
,
9
:
125
140
.

Weigel
B. L.
,
Pfister
C. A.
2021
.
The dynamics and stoichiometry of dissolved organic carbon release by kelp
.
Ecology
,
102
:
e03221
.

Xu
Z.
,
Gao
G.
,
Xu
J.
,
Wu
H.
2017
.
Physiological response of a golden tide alga (Sargassum muticum) to the interaction of ocean acidification and phosphorus enrichment
.
Biogeosciences
,
14
:
671
681
.

Yoro
S.
,
Panagiotopoulos
C.
,
Sempéré
R.
1999
.
Dissolved organic carbon contamination induced by filters and storage bottles
.
Water Research
,
33
:
1956
1959
.

Young
C. S.
,
Gobler
C. J.
2016
.
Ocean acidification accelerates the growth of two bloom-forming macroalgae
.
PLoS ONE
,
11
:
e0155152
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
Handling Editor: Ivan Rodil
Ivan Rodil
Handling Editor
Search for other works by this author on:

Supplementary data