Philippa Kaur

Using a targeted metabolomics approach, we showed that late copepodite stages of the keystone Arctic copepod Calanus glacialis experience important changes in several central energetic pathways following exposure to decreasing pH. These findings shed light on the physiological changes underpinning the effects of OA on fitness related traits such as ingestion rate and metabolic rate previously observed in this species17,18,20.Cellular energy metabolismCellular energy production was altered consistently in both stage CIV and CV, with concentrations of higher energy adenosine phosphates (ATP and ADP) increasing, and concentrations of the lower energy, less-phosphorylated AMP decreasing, with decreasing seawater pH. Moreover, Phospho-L-arginine, which in crustaceans functions as phosphagen in the replenishment of ATP from ADP during transient energy demands32, increased significantly in stage CV. These changes strongly suggest that exposure to low pH affects energy production and expenditure in both developmental stages, although with nuanced differences.NAD+ increased significantly in stage CIV. NAD+ is an essential redox carrier receiving electrons from oxidative processes in the glycolysis, the TCA cycle, and fatty acid oxidation to form NADH. A high NAD+/NADH ratio facilitates higher rates of these reactions and thus potentially higher rates of ATP production (unfortunately, the LC-HRMS could not detect NADH). But most importantly, the produced NADH serves as electron donors to ATP synthesis in the oxidative phosphorylation. For every ATP produced in the oxidative phosphorylation one NADH is oxidised back to NAD+. High rates of ATP production in the oxidative phosphorylation would therefore amass NAD+, as observed in stage CIV. Conversely, ATP production in the glycolysis and TCA cycle consumes NAD+ (9 NAD+ per 4 ATP) and glycolytic ATP production would decrease the NAD+ concentration with decreasing pH.Heterotrophic organisms generally face a trade-off between rate and yield of ATP production. The efficient low rate/high yield production in the TCA cycle/oxidative phosphorylation may prevail under certain circumstances, whereas under other circumstances, the less efficient high rate/low yield production in the glycolysis may predominate33. Because glycolysis and oxidative phosphorylation compete for ADP, the one dominate over the other in terms of rates depending on the substrate being metabolised. In stage CIV copepodites, the TCA cycle pathway was enriched in the MetPA, and metabolites associated with glycolysis and the TCA cycle showed significant changes in their concentrations at decreasing seawater pH. Glucose, the entry point to glycolysis, increased significantly with decreasing pH. High levels of blood glucose (hyperglycemia) have been observed as a general stress response in decapod crustaceans34. Copepods have no circulatory system (although they have a dorsal heart) but might nevertheless react similarly on the cellular level. Along with the significant increase in glucose, lactate decreased significantly with pH in stage CIV. Lactate is an inevitable end product of glycolysis, because lactate dehydrogenase has the highest Vmax of any enzymes in the glycolytic pathway and the Keq for pyruvate to lactate is far in the direction of lactate35. Accordingly, although the glycolysis was not enriched in the MetPA, conceivably because none of its intermediate metabolites were included in the analyses (the protocol did not allow for it), we hypothesise that stage CIV copepodites experience a general down-regulation of glycolysis under decreasing pH. Alternatively, the amassing of glucose and depletion of lactate could also indicate increased gluconeogenesis. Gluconeogenesis occurs during starvation to replenish glycogen stores and ingestion rates did decrease in stage CIV20. But again, we did not target any intermediates in our analyses, and thus cannot firmly conclude on this.Phosphofructokinase-1 is a key regulatory enzyme of glycolysis36. This enzyme is allosterically inhibited by ATP and activated by AMP, and interestingly this regulation is augmented by low pH37,38. Thus, phosphofructokinase-1 could be key to the down-regulation of glycolysis we hypothesise. The fact that we found increasing oxygen consumption with decreasing pH in stage CIV copepodites from the same experiment adds further momentum to this line of thought20. It seems that stage CIV copepodites might experience the so-called Pasteur effect—a decrease in glycolysis at increased levels of oxygen uptake—when exposed to decreasing pH39. Although ATP and AMP were significantly affected also in stage CV, glucose, pyruvate and lactate did not change with decreasing pH, which perhaps indicate absence of the down-regulation of glycolysis we hypothesise for stage CIV. There is, nevertheless, one indication that down-regulation may in fact occur also in this developmental stage. Alpha-glycerophosphate decreased significantly with decreasing pH in stage CV. This molecule is an intermediate in the transfer of electrons from NADH produced by glycolysis in the cytosol to the oxidative phosphorylation in the mitochondria, and decreased concentrations could result from down-regulation of the glycolysis also in stage CV copepodites.The TCA cycle was enriched for stage CIV and most of the measured TCA cycle metabolites (alpha-ketoglutarate, succinate, fumarate, and malate) showed increasing concentrations at decreasing pH. Trigg et al.40 observed a similar increase in concentrations of TCA cycle-related metabolites in the Dungeness crab, Cancer magister (Dana, 1852), at decreased pH and concluded that TCA cycle activity is upregulated under OA. Since NAD+ is the product of the transport of electrons from the TCA cycle to the oxidative phosphorylation in the mitochondria,  the increase in NAD+ concentration we observed in stage CIV could reflect an increase in the flow of electrons from the TCA cycle to the oxidative phosphorylation, and by extension an increase in the energy production by the TCA cycle and the oxidative phosphorylation. There is negative feedback from the TCA cycle to glycolysis through inhibition of phosphofructokinase-1 by citrate, a metabolite of the TCA cycle38. Unfortunately, we did not target citrate in our targeted approach to specifically test this hypothesis, but the amassing of NAD+ do provide additional support to the idea that glycolysis is down-regulated at decreasing pH. Again, there is a less clear picture of how cellular energy metabolism is affected by decreasing pH in stage CV when compared to stage CIV. There was no clear pattern of regulation of TCA metabolites, and the TCA cycle was not enriched in the MetPA. Nevertheless, alpha-ketoglutarate concentrations did increase with decreasing pH in CVs.The glyoxylate/dicarboxylate cycle was also enriched in the pathway analysis, but this is probably also a result of the increases in concentrations of alpha-ketoglutarate, succinate, fumarate, and malate, and we are unable to distinguish it from the TCA cycle based on the set of metabolites analysed.Conclusively, lowered glycolysis due to inhibition of phosphofructokinase-1 and upregulation of the TCA cycle and oxidative phosphorylation at low pH in stage CIV appear plausible causes for the changes in ATP, ADP and AMP concentrations we observed. Alongside these effects, down-regulation of transcription of genes involved in the glycolysis were also present in nauplii of C. glacialis exposed to 35–38 days of low pH conditions16. On the other hand, studies on the acclimatisation and adaptation to OA in another calanoid copepod species, Pseudocalanus acuspes (Giesbrecht, 1881), showed no increase in expression of mitochondrial genes at pHT 7.54, which would have been expected if the TCA cycle or oxidative phosphorylation is upregulated41. Interestingly, De Wit et al.41 also showed natural selection in a large fraction of mitochondrial genes under OA conditions. Even evolutionarily conserved sequences, such as cytochrome oxidase subunit I, were under selection and it was hypothesised that the mitochondrial function of oxidative phosphorylation is a target for natural selection in copepods at low pH41.Besides its role in the transfer of energy from the mitochondria to the cell, ATP is also used to fuel cell homeostasis and active cellular acid–base regulation by activation of ATP-dependent enzymes involved in osmo-ionic- and acid–base regulation. In crustaceans, acid–base status is linked to ion regulation, and is maintained primarily through ion transport mechanisms moving acid and/or base equivalents between the extracellular fluid and the ambient water42. One prominent process in this respect is regulation by Na+/K+-ATPase42,43. While this regulation takes place in the gills of decapod crustaceans43, it is located in the maxillary glands and other specialised organs on the swimming legs of copepods44. Any extensive ATPase mediated pH regulation could have manifested itself by decreasing ATP concentrations, but this is contrary to what we report here. Interestingly, while the pCO2-sensitive isopod Cymodoce truncata (Leach, 1814) is able to maintain its cellular ATP concentration at the expense of the concentration of carbonate anhydrase (an enzyme involved in the cellular transformation of water and CO2 to bicarbonate ions and H+ prior to the ATPase mediated transport of H+ across the cell membrane), the pCO2-tolerant isopod Dynamene bifida (Torelli, 1930) upregulates ATP with no functional compromise to CA concentrations45. Finally, C. glacialis nauplii have shown upregulation of Na+/H+-antiporters independent of ATPase as a response to OA16, which one could hypothesise also may be the case in the copepodites. Arctic populations of the amphipod Gammarus setosus also do not experience increased ATPase activity during OA conditions46. It seems that C. glacialis faces OA without any ATP dependent acid/base regulation activity.Glycolysis is the first step of catabolism of carbohydrates for the production of energy. When down-regulating glycolysis the copepods may be increasingly dedicated to catabolism of amino acids e.g. through oxidative deamination of glutamate and/or catabolism of fatty acids through beta-oxidation to produce the energy they require21. Both lead to the production of molecules entering the TCA cycle and ultimately the oxidative phosphorylation for energy production in the mitochondria.Amino acid metabolismOf the free amino acids which were significantly affected by decreasing pH, the majority decreased in concentration, for both stage CIV and CV copepodites. This could be an indication of changes in protein synthesis at decreasing pH. Supporting this idea, biosynthesis of aminoacyl-tRNA was indicated as significantly enriched in the MetPA in both stage CIV and CV. Aminoacyl-tRNA partakes in the elongation of the protein amino acid chain during protein synthesis and the enrichment was most likely due to the changes in concentration of the many amino acids tested. One probable cause of protein synthesis is the increased demands of enzymes needed to handle stress at low pH, including for example enzymes involved in acid–base- and osmo-regulation or regulation of energy production. Increased protein synthesis caused by OA conditions has been observed in larvae of the purple sea urchin Strongylocentrotus purpuratus (O.F. Müller, 1776), where in vivo rates of protein synthesis and ion transport increased ∼50%47. Costs of protein synthesis are high and have shown to constitute a major part of copepod metabolic demand48 and we did observe significant increases in metabolic rate in copepodite stage CIV from the same experiment20 giving further credit to the idea that protein synthesis was upregulated.An alternate but not mutually exclusive explanation is that the copepods experience increased amino acid catabolism under OA. Glutamate increased in stage CIV accompanied by a significant increase in alpha-ketoglutarate in both stage CIV and CV. Alpha-ketoglutarate is part of the metabolic pathway of glutamine, glutamate and arginine in which glutamate acts as an intermediate in catabolism of these amino acids when it is deaminated to alpha-ketoglutarate to enter the TCA cycle49. Glutamate metabolism (in conjunction with alanine and aspartate metabolism) was significantly enriched in the MetPA in both stage CIV and CV, and these changes could be taken as an indication of a shift towards amino acid catabolism with decreasing pH. The key enzyme catalysing the oxidative deamination of glutamate is glutamate dehydrogenase (GDH), which functions in both directions: deamination of glutamate to form alpha-ketoglutarate or formation of glutamate from alpha-ketoglutarate. Studies on the ribbed mussel, Modiolus dernissus (Dillwyn, 1817), have shown that the balance of this action is strongly pushed towards deamination when pH decreases from 8.0 to 7.550. GDH is activated by ADP, and one could argue that the increase in ADP we observed would work against this shift, but ADP activates GDH mainly in the glutamate forming direction51. The other measured amino acids enter the TCA cycle at different positions we unfortunately could not target in our analyses. Glutamate also partakes in the arginine biosynthesis pathway in which it is transformed to ornithine to enter the urea cycle. Arginine biosynthesis was enriched in the MetPA and it is therefore possible that decreasing pH also changes amino acid catabolism to increase urea excretion. Decreasing pH has a similar depressing effect on amino acid concentration in the gills of the shore crab Carcinus maenas (Linnaeus, 1758) which also has been interpreted as a sign of increased protein catabolism52. Hammer and colleagues52 argued that this increase in catabolism served to buffer H+ by supplying nitrogen to NH4 formation in the cells. All in all, we hypothesise that increased amino acid catabolism, possibly driven by changes in GDH activity, and the down-regulation of glycolysis by inhibition of phosphofructokinase-1 may be major drivers of a shift from carbohydrate metabolism towards catabolism of amino acids.D-glutamine/D-glutamate metabolism was highly enriched in the MetPA in both developmental stages. Several studies show enriched D-glutamine/D-glutamate metabolism in crustaceans [e.g. 53], but they offer no explanation of its function or the reason why it is enriched. While D-glutamate act in neurotransmission, this action is evolutionarily restricted to ctenophores, and biochemical measurements of D-amino acid concentrations have shown absence of D-glutamate in crustaceans54,55.We observed no changes in concentrations of 8-oxy-2-deoxyguanosine, a product of DNA oxidation. Furthermore, regulation of cellular response to oxidative stress is down-regulated in C. glacialis nauplii16, and OA may not induce oxidative stress in C. glacialis.Fatty acid metabolismBesides their importance in energy storage as wax esters, fatty acids are involved in many central processes in cells, most prominently through their function as cell membrane building blocks. Many fatty acids are obtained from the diet but some longer chain fatty acids, such as 20:1n-9 are synthesised de novo in copepods56. Stage CV copepodites experienced increases in most of the targeted free fatty acids (18 of 21) with decreasing pH. Only one of those 18 increased significantly, but since the direction of change were the same in all, we argue that the pattern of change does merit consideration. Conspicuous exceptions were eicosapentaenoic acid (EPA) 20:5n-3 and docosahexaenoic acid (DHA) 22:6n-3, which both decreased significantly. The only other study (to our knowledge) of metabolomic effects of environmental changes in copepods showed the exact same response to starvation in a mix of C. finmarchicus and C. helgolandicus stage CV copepodites, with most fatty acids increasing while EPA and DHA decreased in concentration57. EPA and DHA are key marine polyunsaturated fatty acids (PUFAs) exclusively produced by marine algae. They contribute a major fraction of the fatty acids of cell membrane phospholipids58, and zooplankton reproductive production is highly dependent on especially EPA59. EPA and DHA are key for cell membrane fluidity, which for calanoid copepods is especially important during diapause in the deep during copepodite stage CV60. They have also been linked to diapause buoyancy control, and are selectively metabolized in diapausing copepodites61. The importance of EPA and DHA for cell membrane integrity may be central for the changes we observed. Glycerol-3-phosphate, the precursor for the glycerol backbone of cell membrane phospholipids also decreased significantly and it seems decreasing pH could affect cell membrane turnover.Changing fatty acid concentration could be due to either a change in lipid intake from feeding or increased fatty acid catabolism. While ingestion rates decreased in stage CIV, they were unchanged in stage CV with decreasing pH20. Also, Thalassiosira weissflogii (Grunow) G.Fryxell & Hasle, 1977, the diatom we fed to the copepods, is rich in 16:0, 16:1n-7 and EPA59. The concentrations of 16:0 and 16:1n-7 increased, whereas EPA concentration decreased. If fatty acid concentrations reflected feeding, we would have seen increased concentrations of all three. We therefore believe that the general increases in concentrations of free fatty acids were caused by increasing catabolism of the wax esters stored in stage CV. It may be that due to the metabolic reconfiguration to enter hibernation, stage CV copepodites are already committed to the catabolism of fatty acids through beta-oxidation, and stored wax esters are being hydrolysed to increase the availability of free fatty acids for energy production. Mayor and colleagues57 arrived at the same conclusion. We hypothesise that stress due to low pH increases the organism’s energetic demands, but carbohydrates are not used to accommodate these demands due to the down-regulation of the glycolysis, rather demands are met by hydrolysing and metabolising wax esters in stage CV. The further ramifications of future OA could therefore be a less efficient build-up of wax esters so important for hibernation in this species.Finally, besides their importance for cell membrane fluidity, EPA and DHA are important precursors for eicosanoid endocrine hormones. These hormones are important regulators of, among other processes, ion flux62. As mentioned above, acid base regulation is coupled to osmoregulation in crustaceans42, and the decrease in concentrations of these two specific fatty acids, when all other fatty acid concentrations increased might represent an indication for changing endocrine hormone production to counter adverse whole-organism effects of OA.Changes in metabolite concentrations cannot be directly translated into changes in the rate of the processes they are involved in. However, they do pin-point processes which are affected by the imposed environmental changes. Also, in our analyses we targeted a limited range of molecules. In that respect OA could inflict changes in other important metabolic pathways we did not investigate. The absence of specific biochemical pathways in our analyses and discussion should therefore not be taken as indication that these are not implicated in this species responses to OA.From our previously published study on copepodites from the same incubations, we know that high pCO2/low pH conditions have detrimental effects on the balance between energy input (ingestion) and energy expenditure (metabolism) in stage CIV copepodites but not in stage CV copepodites20. The effects we report here help in this sense to shed light on the metabolic origin of the rather severe effects on energy balance we observed in stage CIV copepodites and the difference in response between stage CIV and CV20. Copepods develop through six nauplii and five copepodite stages before maturation, and while previous studies show negligible effects in stage CV and adults17,18,20, any effects in any developmental stage along the way will affect the fitness of the individual and the recruitment to the population as a whole. In addition, the enhanced fatty acid metabolism observed in stage CV needs further investigation, to determine the magnitude of the fitness implications of the energy diverted away from energy storage for hibernation. More

Odum, E. P. The strategy of ecosystem development. Science 164, 262–270 (1969).Article 
ADS 
CAS 

Google Scholar 
Gorham, E., Vitousek, P. M. & Reiners, W. A. The regulation of element budgets over the course of terrestrial ecosystem succession. Annu. Rev. Ecol. Syst. 10, 53–84 (1979).Article 
CAS 

Google Scholar 
Corman, J. R. et al. Foundations and frontiers of ecosystem science: Legacy of a classic paper (Odum 1969). Ecosystems 22, 1160–1172. https://doi.org/10.1007/s10021-018-0316-3 (2019).Article 

Google Scholar 
Santín, C. et al. Towards a global assessment of pyrogenic carbon from vegetation fires. Glob. Change Biol. 22, 76–91. https://doi.org/10.1111/gcb.12985 (2016).Article 
ADS 

Google Scholar 
Kominoski, J. S., Gaiser, E. E. & Baer, S. G. Advancing theories of ecosystem development through long-term ecological research. Bioscience 68, 554–562. https://doi.org/10.1093/biosci/biy070 (2018).Article 

Google Scholar 
Balch, J. K., Bradley, B. A., D’Antonio, C. M. & Gómez-Dans, J. Introduced annual grass increases regional fire activity across the arid western USA (1980–2009). Glob. Change Biol. 19, 173–183. https://doi.org/10.1111/gcb.12046 (2013).Article 
ADS 

Google Scholar 
Abatzoglou, J. T. & Kolden, C. A. Climate change in Western US deserts: Potential for increased wildfire and invasive annual grasses. Rangeland Ecol. Manag. 64(5), 471–478 (2011).Article 

Google Scholar 
Shi, H. et al. Historical cover trends in a sagebrush steppe ecosystem from 1985 to 2013: Links with climate, disturbance, and management. Ecosystems 21, 913–929. https://doi.org/10.1007/s10021-017-0191-3 (2018).Article 

Google Scholar 
Seyfried, M. S. & Wilcox, B. P. Scale and the nature of spatial variability: Field examples having implications for hydrologic modeling. Water Resour. Res. 31, 173–184. https://doi.org/10.1029/94WR02025 (1995).Article 
ADS 

Google Scholar 
Gasch, C. K., Huzurbazar, S. V. & Stahl, P. D. Description of vegetation and soil properties in sagebrush steppe following pipeline burial, reclamation, and recovery time. Geoderma 265, 19–26. https://doi.org/10.1016/j.geoderma.2015.11.013 (2016).Article 
ADS 

Google Scholar 
Huber, D. P. et al. Vegetation and precipitation shifts interact to alter organic and inorganic carbon storage in desert soils. Ecosphere 10, e02655. https://doi.org/10.1002/ecs2.2655 (2019).Article 

Google Scholar 
Knight, D. H., Jones, G. P., Reiners, W. A. & Romme, W. H. Mountains and Plains: The Ecology of Wyoming Landscapes (Yale University Press, 2014).
Google Scholar 
Patton, N. R., Lohse, K. A., Seyfried, M. S., Godsey, S. E. & Parsons, S. Topographic controls on soil organic carbon on soil mantled landscapes. Sci. Rep. 9, 6390. https://doi.org/10.1038/s41598-019-42556-5 (2019).Article 
ADS 
CAS 

Google Scholar 
Schwabedissen, S. G., Lohse, K. A., Reed, S. C., Aho, K. A. & Magnuson, T. S. Nitrogenase activity by biological soil crusts in cold sagebrush steppe ecosystems. Biogeochemistry 134, 57–76. https://doi.org/10.1007/s10533-017-0342-9 (2017).Article 
CAS 

Google Scholar 
You, Y. et al. Biological soil crust bacterial communities vary along climatic and shrub cover gradients within a sagebrush steppe ecosystem. Front. Microbiol. 12, 2365. https://doi.org/10.3389/fmicb.2021.569791 (2021).Article 

Google Scholar 
Burke, I. C., Reiners, W. A. & Olson, R. K. Topographic control of vegetation in a mountain big sagebrush steppe. Vegetation 84, 77–86 (1989).Article 

Google Scholar 
Poulos, M. J., Pierce, J. L., Flores, A. N. & Benner, S. G. Hillslope asymmetry maps reveal widespread, multi-scale organization. Geophys. Res. Lett. 39, 6. https://doi.org/10.1029/2012GL051283 (2012).Article 

Google Scholar 
Smith, T. & Bookhagen, B. Climatic and biotic controls on topographic asymmetry at the global scale. J. Geophys. Res.: Earth Surf. 126, e2020JF005692. https://doi.org/10.1029/2020JF005692Received22 (2021).Article 
ADS 

Google Scholar 
Seyfried, M., Link, T., Marks, D. & Murdock, M. Soil temperature variability in complex terrain measured using fiber-optic distributed temperature sensing. Vadose Zone J. 15, 6. https://doi.org/10.2136/vzj2015.09.0128 (2016).Article 

Google Scholar 
Chambers, J. C. et al. Resilience and resistance of sagebrush ecosystems: Implications for state and transition models and management treatments. Rangel. Ecol. Manage. 67, 440–454. https://doi.org/10.2111/REM-D-13-00074.1 (2014).Article 

Google Scholar 
Chambers, J. C. et al. Operationalizing resilience and resistance concepts to address invasive grass-fire cycles. Front. Ecol. Evol. 7, 2369. https://doi.org/10.3389/fevo.2019.00185 (2019).Article 

Google Scholar 
Boehm, A. R. et al. Slope and aspect effects on seedbed microclimate and germination timing of fall-planted seeds. Rangel. Ecol. Manage. 75, 58–67. https://doi.org/10.1016/j.rama.2020.12.003 (2021).Article 

Google Scholar 
Sankey, J. B., Germino, M. J., Sankey, T. T. & Hoover, A. N. Fire effects on the spatial patterning of soil properties in sagebrush steppe, USA: A meta-analysis. Int. J. Wildl. Fire 21, 545–556. https://doi.org/10.1071/WF11092 (2012).Article 

Google Scholar 
Fellows, A., Flerchinger, G., Seyfried, M. S. & Lohse, K. A. Rapid recovery of mesic mountain big sagebrush gross production and respiration following prescribed fire. Ecosystems 21, 1283–1294. https://doi.org/10.1007/s10021-017-0218-9 (2018).Article 

Google Scholar 
Vega, S. P. et al. Interaction of wind and cold-season hydrologic processes on erosion from complex topography following wildfire in sagebrush steppe. Earth Surf. Process. Landforms https://doi.org/10.1002/esp.4778 (2019).Article 

Google Scholar 
Xie, J., Li, Y., Zhai, C., Li, C. & Lan, Z. CO2 absorption by alkaline soils and its implication to the global carbon cycle. Environ. Geol. 56, 953–961 (2009).Article 
ADS 
CAS 

Google Scholar 
Stanbery, C., Pierce, J. L., Benner, S. G. & Lohse, K. On the rocks: Quantifying storage of inorganic soil carbon on gravels and determining pedon-scale variability. CATENA 157, 436–442. https://doi.org/10.1016/j.catena.2017.06.011 (2017).Article 
CAS 

Google Scholar 
Stanbery, C. et al. Controls on the presence and concentration of soil inorganic carbon in a semi-arid watershed. CATENA https://doi.org/10.2139/ssrn.4267018 (2023).Article 

Google Scholar 
Cerling, T. E. & Quade, J. Stable carbon and oxygen isotopes in soil carbonates. Geophys. Monogr. 78, 217–231 (1993).ADS 

Google Scholar 
Tappa, D. J., Kohn, M. J., McNamara, J. P., Benner, S. G. & Flores, A. N. Isotopic composition of precipitation in a topographically steep, seasonally snow-dominated watershed and implications of variations from the global meteoric water line. Hydrol. Process. 30, 4582–4592. https://doi.org/10.1002/hyp.10940 (2016).Article 
ADS 
CAS 

Google Scholar 
Salomons, W., Goudie, A. & Mook, W. G. Isotopic composition of calcrete deposits from Europe, Africa and India. Earth Surf. Process. 3, 43–57. https://doi.org/10.1002/esp.3290030105 (1978).Article 
CAS 

Google Scholar 
Salomons, W. & Mook, W. G. In Handbook of Environmental Isotope Geochemistry (eds P. Fritz & J. Fontes) Ch. 6, 241–269 (Elsevier, 1986).Bodí, M. B. et al. Wildland fire ash: Production, composition and eco-hydro-geomorphic effects. Earth Sci. Rev. 130, 103–127. https://doi.org/10.1016/j.earscirev.2013.12.007 (2014).Article 
ADS 
CAS 

Google Scholar 
Kéraval, B. et al. Soil carbon dioxide emissions controlled by an extracellular oxidative metabolism identifiable by its isotope signature. Biogeosciences 13, 6353–6362. https://doi.org/10.5194/bg-13-6353-2016 (2016).Article 
ADS 
CAS 

Google Scholar 
Goforth, B. R., Graham, R. C., Hubbert, K. R., Zanner, C. W. & Minnich, R. A. Spatial distribution and properties of ash and thermally altered soils after high-severity forest fire, southern California. Int. J. Wildland Fire 14, 343–354 (2005).Article 

Google Scholar 
Glossner, K. L. et al. Long-term suspended sediment and particulate organic carbon yields from the Reynolds Creek Experimental Watershed and Critical Zone Observatory. Hydrol. Process. 36, e14484. https://doi.org/10.1002/hyp.14484 (2022).Article 
CAS 

Google Scholar 
Seyfried, M. S. et al. Reynolds creek experimental watershed and critical zone observatory. Vadoze Zone J. 17, 180129. https://doi.org/10.2136/vzj2018.07.0129 (2018).Article 
CAS 

Google Scholar 
McIntyre, D. H. Cenozoic geology of the Reynolds Creek Experimental Watershed, Owyhee County, Idaho (Idaho Bureau of Mines and Geology, 1972).Earth Resources Observation and Science (EROS) Center, U. Image of the week: Burned Area Analysis for the Soda Fire, Idaho, https://eros.usgs.gov/media-gallery/image-of-the-week/burned-area-analysis-the-soda-fire-idaho (2015).Jenny, H. Factors of Soil Formation (McGraw-Hill, 1941).Book 

Google Scholar 
Kormos, P. R. et al. 31 years of hourly spatially distributed air temperature, humidity, and precipitation amount and phase from Reynolds Critical Zone Observatory. Earth Syst. Sci. Data 10, 1197–1205. https://doi.org/10.5194/essd-10-1197-2018 (2018).Article 
ADS 

Google Scholar 
Thomas, G. W. In Methods in Soil Analysis. Part 3. Chemical Methods (ed Sparks, D. L. ) (Soil Science Society of America and American Society of Agronomy, 1996).Brodie, C. R. et al. Evidence for bias in C and N concentrations and δ13C composition of terrestrial and aquatic organic materials due to pre-analysis acid preparation methods. Chem. Geol. 282, 67–83. https://doi.org/10.1016/j.chemgeo.2011.01.007 (2011).Article 
ADS 
CAS 

Google Scholar 
Patton, N. P., Lohse, K. A., Seyfried, M. S., Will, R. & Benner, S. G. Lithology and coarse fraction adjusted bulk density estimates for determining total organic carbon stocks in dryland soils. Geoderma 337, 844–852. https://doi.org/10.1016/j.geoderma.2018.10.036 (2019).Article 
ADS 
CAS 

Google Scholar 
McGuire, L. A., Rasmussen, C., Youberg, A. M., Sanderman, J. & Fenerty, B. Controls on the Spatial distribution of near-surface pyrogenic carbon on hillslopes 1 year following wildfire. J. Geophys. Res.: Earth Surf. 126, e2020JF005996. https://doi.org/10.1029/2020JF005996 (2021).Article 
ADS 

Google Scholar 
Jiménez-González, M. A. et al. Spatial distribution of pyrogenic carbon in Iberian topsoils estimated by chemometric analysis of infrared spectra. Sci. Total Env. 790, 148170. https://doi.org/10.1016/j.scitotenv.2021.148170 (2021).Article 
CAS 

Google Scholar 
Baldock, J. A. et al. Quantifying the allocation of soil organic carbon to biologically significant fractions. Soil Res. 51, 561–576. https://doi.org/10.1071/SR12374 (2013).Article 
CAS 

Google Scholar 
Sanderman, J. et al. Soil organic carbon fractions in the Great Plains of the United States: An application of mid-infrared spectroscopy. Biogeochemistry 156, 97–114. https://doi.org/10.1007/s10533-021-00755-1 (2021).Article 
CAS 

Google Scholar 
Sherrod, L. A., Dunn, G., Peterson, G. A. & Kolberg, R. L. Inorganic carbon analysis by modified pressure-calcimeter method. Soil Sci. Soc. Am. J. 66, 299–305 (2002).Article 
ADS 
CAS 

Google Scholar 
Mikutta, R., Kleber, M., Kaiser, K. & Jahn, R. Review. Soil Sci. Soc. Am. J. 69, 120–135. https://doi.org/10.2136/sssaj2005.0120 (2005).Article 
ADS 
CAS 

Google Scholar 
Risk, D., Nickerson, N., Creelman, C., McArthur, G. & Owens, J. Forced Diffusion soil flux: A new technique for continuous monitoring of soil gas efflux. Agric. For. Meteorol. 151, 1622–1631. https://doi.org/10.1016/j.agrformet.2011.06.020 (2011).Article 
ADS 

Google Scholar 
Fierer, N. & Schimel, J. P. Effects of drying–rewetting frequency on soil carbon and nitrogen transformations. Soil Biol. Biochem. 34, 777–787. https://doi.org/10.1016/S0038-0717(02)00007-X (2002).Article 
CAS 

Google Scholar 
Dane, J. H., Topp, G. C. & Campbell, G. S. In Methods of Soil Analysis: Physical Methods. Vol. 4 (ed SSSA) 721–738 (2002). More

Takoutsing, B. et al. Assessment of soil health indicators for sustainable production of maize in smallholder farming systems in the highlands of Cameroon. Geoderma 276, 64–73 (2016).Article 
ADS 
CAS 

Google Scholar 
Wenzel, W. W. et al. Soil and land use factors control organic caron status and accumulation in agricultural soils of Lower Austria. Geoderma 409, 115595. https://doi.org/10.1016/j.geoderma.2021.115595 (2022).Article 
ADS 
CAS 

Google Scholar 
Rinot, O., Levy, G. J., Steinberger, Y., Svoray, T. & Eshel, G. Soil health assessment: A critical review of current methodologies and a proposed new approach. Sci. Total. Environ. 648, 1484–1491 (2019).Article 
ADS 
CAS 

Google Scholar 
Veum, K. S., Sudduth, K. A., Kremer, R. J. & Kitchen, R. (2017) Sensor data fusion for soil health assessment. Geoderma 305, 53–61 (2017).Article 
ADS 

Google Scholar 
Nunes, M. R., Van Es, H. M., Schindelbeck, R., Ristow, A. J. & Ryan, M. No-till and cropping system diversification improve soil health and crop yield. Geoderma 328, 30–43 (2018).Article 
ADS 
CAS 

Google Scholar 
Chhipaa, V., Stein, A., Shankar, H., George, K. J. & Alidoost, F. Assessing and transferring soil health information in a hilly terrain. Geoderma 343, 130–138 (2019).Article 
ADS 

Google Scholar 
Oliver, D. P., Bramley, R. G. V., Riches, D., Porter, I. & Edwards, J. A review of soil physical and chemical properties as indicators of soil quality in Australian viticulture. Aust. J. Grape Wine Res. 19, 129–139 (2013).Article 
CAS 

Google Scholar 
Riches, D. et al. Review: soil biological properties as indicators of soil quality in Australian viticulture. Aust. J. Grape Wine Res. 19, 311–323 (2013).CAS 

Google Scholar 
Ritz, K., Black, H. I. J., Campbell, C. D., Harris, J. A. & Wood, C. Selecting biological indicators for monitoring soils: a framework for balancing scientific opinion to assist policy development. Ecol. Ind. 9, 1212–1221 (2009).Article 
CAS 

Google Scholar 
Zhuo, Z., Kirchner, I., Pfahl, S. & Cubasch, U. Climate impact of volcanic eruptions: the sensitivity to eruption season and latitude in MPI-ESM ensemble experiments. Atmos. Chem. Phys. 21, 13425–13442 (2021).Article 
ADS 
CAS 

Google Scholar 
Griffiths, B. S., Bonkowski, M., Roy, J. & Ritz, K. Functional stability, substrate utilisation and biological indicators of soils following environmental impacts. Appl. Soil Ecol. 16(1), 49–61 (2001).Article 

Google Scholar 
Avidano, L., Gamalero, E., Cossa, G. P. & Carraro, E. Characterization of soil health in an Italian polluted site by using microorganisms as bioindicators. Appl. Soil Ecol. 30(1), 21–33 (2005).Article 

Google Scholar 
Pattison, A. B. et al. Development of key soil health indicators for the Australian banana industry. Appl. Soil Ecol. 40(1), 155–164 (2008).Article 

Google Scholar 
Damsma, K. M., Rose, M. T. & Cavagnaro, T. R. Landscape scale survey of indicators of soil health in grazing systems. Soil Res. 53(2), 154–167 (2015).Article 

Google Scholar 
Fine, A. K., van Es, H. M. & Schindelbeck, R. R. Statistics, scoring functions, and regional analysis of a comprehensive soil health database. Soil Sci. Soc. Am. J. 81(3), 589–601 (2016).Article 

Google Scholar 
Roper, W. R., Osmond, D. L., Heitman, J. L., Wagger, M. G. & Reberg-Horton, S. C. Soil Health indicators do not differentiate among agronomic management systems in North Carolina soils. Soil Sci. Soc. Am. J. 81(4), 828–843 (2016).Article 

Google Scholar 
Li, Z. et al. Rapid diagnosis of agricultural soil health: A novel soil health index based on natural soil productivity and human management. J. Environ. Manage. 277, 111402. https://doi.org/10.1016/j.jenvman.2020.111402 (2021).Article 

Google Scholar 
Oren, A. & Steinberger, Y. Catabolic profiles of soil fungal communities along a geographic climatic gradient in Israel. Soil Biol. Biochem. 40, 2578–2587 (2008).Article 
CAS 

Google Scholar 
Yu, J., Glazer, N. & Steinberger, Y. Carbon utilization, microbial biomass, and respiration in biological soil crusts in the Negev Desert. Biol. Fert. Soils 50, 285–293 (2014).Article 
CAS 

Google Scholar 
Van der Putten, W. H. et al. Plant–soil feedbacks: the past, the present and future challenges. J. Ecol. 101, 265–276 (2013).Article 

Google Scholar 
Sherman, C. & Steinberger, Y. Microbial functional diversity associated with plant litter decomposition along a climatic gradient. Microb. Ecol. 64, 399–415 (2012).Article 
CAS 

Google Scholar 
Dwivedi, V. & Soni, P. A review on the role of soil microbial biomass in eco-restoration of degraded ecosystem with special reference to mining areas. J. Appl. Nat. Sci. 3(1), 151–158 (2011).Article 

Google Scholar 
Barreiro, A., Martín, A., Carballas, T. & Díaz-Raviña, M. Long-term response of soil microbial communities to fire and fire-fighting chemicals. Biol. Fertil. Soils 52, 963–975 (2016).Article 
CAS 

Google Scholar 
Soil Science Division Staff. Soil survey manual. In USDA Handbook 18 (ed. Ditzler, C., Scheffe, K. & Monger, H.C.). (Washington, G. P. O., 2017).Campbell, C. D., Chapman, S. J., Cameron, C. M., Davidson, M. S. & Potts, J. M. A rapid microtiter plate method to measure carbon dioxide evolved from carbon substrate amendments so as to determine the physiological profiles of soil microbial communities by using whole soil. Appl. Environ. Microbiol. 69, 3593–3599 (2003).Article 
ADS 
CAS 

Google Scholar 
Anderson, J. P. E. & Domsch, K. H. Physiological method for quantitative measurement of microbial biomass in soils. Soil Biol. Biochem. 10, 215–221 (1978).Article 
CAS 

Google Scholar 
Creamer, R. E., Stone, D., Berry, P. & Kuiper, I. Measuring respiration profiles of soil microbial communities across Europe using MicroResp™ method. Appl. Soil Ecol. 97, 36–43 (2016).Article 

Google Scholar 
Oren, A. & Steinberger, Y. Coping with artifacts induced by CaCO3–CO2–H2O equilibria in substrate utilization profiling of calcareous soils. Soil Biol. Biochem. 40, 2569–2577 (2008).Article 
CAS 

Google Scholar 
Zak, J. C., Willig, M. R., Howard, D. L. & Wildman, G. Functional diversity of microbial communities: A quantitative approach. Soil Biol. Biochem. 26(9), 1101–1108 (1994).Article 

Google Scholar 
Hotelling, H. The most predictable criterion. J. Educ. Psychol. 26, 139–142 (1935).Article 

Google Scholar 
Morrison, D. F. Multivariate Statistical Methods 2nd edn. (McGraw-Hill, 1976).MATH 

Google Scholar 
Rencher, A. C. Methods of Multivariate Analysis (Wiley, Uk, 1995).MATH 

Google Scholar 
IBM Corp. Released 2020. IBM SPSS Statistics for Windows, Version 27.0. (Armonk, NY: IBM Corp., 2020)R Core Team. A language and environment for statistical computing (R Foundation for Statistical Computing, 2021).
Google Scholar 
Bartoń K. MuMIn: Multi-Model Inference. R package version 1.46.0, https://CRAN.R-project.org/package=MuMIn, 2022.Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).Book 
MATH 

Google Scholar 
Johnson, J. B. & Omland, K. S. Model selection in ecology and evolution. Trends Ecol. Evol. 19, 101–108 (2004).Article 

Google Scholar 
Kiryushin, V. I. The management of soil fertility and productivity of agrocenoses in adaptive-landscape farming systems. Eurasian Soil Sci. 52, 1137–1145 (2019).Article 
ADS 
CAS 

Google Scholar 
Hermans, S. M. et al. Using soil bacterial communities to predict physic-chemical variables and soil quality. Microbiome 8, 79 (2020).Article 
CAS 

Google Scholar 
Fierer, N. & Jackson, R. B. The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 103, 626–631 (2006).Article 
ADS 
CAS 

Google Scholar 
Frey, S. D., Drijber, R., Smith, H. & Melillo, J. M. Microbial biomass, functional capacity, and community structure after 12 years of soil warming. Soil Biol. Biochem. 40, 2904–2907 (2008).Article 
CAS 

Google Scholar 
Powlson, D. S., Brookes, P. C. & Christensen, B. T. Measurement of soil microbial biomass provides an early indication of changes in total soil organic matter due to straw incorporation. Soil Biol. Biochem. 19, 159–164 (1987).Article 
CAS 

Google Scholar 
Brookes, P. C. The use of microbial parameters in monitoring soil pollution by heavy metals. Biol. Fertil. Soils. 19, 269–279 (1995).Article 
CAS 

Google Scholar 
Hermans, S. M. et al. Bacteria as emerging indicators of soil condition. Appl. Environ. Microbiol. 83, e02826-e2916. https://doi.org/10.1128/AEM.02826-16 (2016).Article 

Google Scholar 
Liddicoat, C. et al. Can bacterial indicators of a grassy woodland restoration inform ecosystem assessment and microbiota-mediated human health?. Environ. Int. 129, 105–117 (2019).Article 

Google Scholar 
Jeanne, T., Parent, S. -É. & Hogue, R. Using a soil bacterial species balance index to estimate potato crop productivity. PLoS ONE 14, e0214089. https://doi.org/10.1371/journal.pone.0214089 (2019).Article 
CAS 

Google Scholar 
Taylor, B. R. & Parkinson, D. Respiration and mass loss rates of aspen and pine leaf litter decomposing in laboratory microcosms. Can. J. Bot. 66, 1948–1959 (1988).Article 

Google Scholar 
Wardle, D. A. & Parkinson, D. Response of the soil microbial biomass to glucose, and selective inhibitors, across a soil moisture gradient. Soil Biol. Biochem. 22, 825–834 (1990).Article 
CAS 

Google Scholar 
Baldrian, P., Merhautová, V., Petránková, M., Cajthaml, T. & Snajdr, J. Distribution of microbial biomass and activity of extracellular enzymes in a hardwood forest soil reflect soil moisture content. Appl. Soil Ecol. 46, 177–182 (2010).Article 

Google Scholar 
Holland, T. C. et al. The response of soil biota to water availability in vineyards. Pedobiol. Int. J. Soil Biol. 56, 9–14 (2013).
Google Scholar 
Yu, J. & Steinberger, Y. Vertical distribution of microbial-community functionality under the canopies of Zygophyllum dumosum and Hammada scoparia in the Negev Desert. Microb. Ecol. 62, 218–227 (2011).Article 

Google Scholar 
Wardle, D. A. & Parkinson, D. Interaction between microclimatic variables and the soil microbial biomass. Biol. Fertil. Soils. 9, 273–280 (1990).Article 

Google Scholar  More

Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Change 3, 919–925 (2013).Article 
ADS 

Google Scholar 
Wiens, J. J. Climate-related local extinctions are already widespread among plant and animal species. PLoS Biol. 14, e2001104 (2016).Article 

Google Scholar 
Bellard, C., Bertelsmeier, C., Leadley, P., Thuiller, W. & Courchamp, F. Impacts of climate change on the future of biodiversity. Ecol. Lett. 15, 365–377 (2012).Article 

Google Scholar 
Assis, J., Serrão, E. A., Duarte, C. M., Fragkopoulou, E. & Krause-Jensen, D. Major expansion of marine forests in a warmer Arctic. Front. Mar. Sci. 9, 850368 (2022).Article 

Google Scholar 
Assis, J. et al. Major shifts at the range edge of marine forests: The combined effects of climate changes and limited dispersal. Sci. Rep. 7(44348), 1–10 (2017).CAS 

Google Scholar 
O’Leary, J. K. et al. The resilience of marine ecosystems to climatic disturbances. BioScience. https://doi.org/10.1093/biosci/biw161 (2017).Article 

Google Scholar 
Steneck, R. S. et al. Kelp forest ecosystems: Biodiversity, stability, resilience and future. Environ. Conserv. 29, 436–459 (2002).Article 

Google Scholar 
Filbee-Dexter, K. & Scheibling, R. E. Detrital kelp subsidy supports high reproductive condition of deep-living sea urchins in a sedimentary basin. Aquat. Biol. 23, 71–86 (2014).Article 

Google Scholar 
Filbee-Dexter, K. Ocean forests hold unique solutions to our current environmental crisis. One Earth https://doi.org/10.1016/j.oneear.2020.05.004 (2020).Article 

Google Scholar 
Krumhansl, K. A. & Scheibling, R. E. Production and fate of kelp detritus. Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps09940 (2012).Article 

Google Scholar 
Edwards, M. S. & Hernández-Carmona, G. Delayed recovery of giant kelp near its southern range limit in the North Pacific following El Niño. Mar. Biol. 147, 273–279 (2005).Article 

Google Scholar 
Cavanaugh, K. C., Reed, D. C., Bell, T. W., Castorani, M. C. N. & Beas-Luna, R. Spatial variability in the resistance and resilience of giant kelp in southern and Baja California to a multiyear heatwave. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00413 (2019).Article 

Google Scholar 
Butler, C. L., Lucieer, V. L., Wotherspoon, S. J. & Johnson, C. R. Multi-decadal decline in cover of giant kelp Macrocystis pyrifera at the southern limit of its Australian range. Mar. Ecol. Prog. Ser. 653, 1–18 (2020).Article 
ADS 

Google Scholar 
Martínez, B. et al. Distribution models predict large contractions of habitat-forming seaweeds in response to ocean warming. Divers. Distrib. 24, 1350–1366 (2018).Article 

Google Scholar 
Bell, T. W., Allen, J. G., Cavanaugh, K. C. & Siegel, D. A. Three decades of variability in California’s giant kelp forests from the Landsat satellites. Remote Sens. Environ. 238, 110811 (2020).Article 
ADS 

Google Scholar 
Mann, M. E. & Emanuel, K. A. Atlantic Hurricane trends linked to climate change. Eos 87, 233–241 (2006).Article 
ADS 

Google Scholar 
Jensen, J. R., Estes, J. E. & Tinney, L. Remote sensing techniques for kelp surveys. Photogramm. Eng Remote Sens. 46, 743–755 (1980).
Google Scholar 
Cavanaugh, K. C. et al. A review of the opportunities and challenges for using remote sensing for management of surface-canopy forming kelps. Front. Mar. Sci. https://doi.org/10.3389/fmars.2021.753531 (2021).Article 

Google Scholar 
Cavanaugh, K. C., Siegel, D. A., Reed, D. C. & Dennison, P. E. Environmental controls of giant-kelp biomass in the Santa Barbara Channel, California. Mar. Ecol. Prog. Ser. 429, 1–17 (2011).Article 
ADS 

Google Scholar 
Kadhim, M. A. & Abed, M. H. Convolutional neural network for satellite image classification. Stud. Comput. Intell. 830, 165–178 (2020).Article 

Google Scholar 
Segal-Rozenhaimer, M., Li, A., Das, K. & Chirayath, V. Cloud detection algorithm for multi-modal satellite imagery using convolutional neural-networks (CNN). Remote Sens. Environ. 237, 111446 (2020).Article 
ADS 

Google Scholar 
Canonico, G. et al. Global observational needs and resources for marine biodiversity. Front. Mar. Sci. 6, 367 (2019).Article 

Google Scholar 
LeCun, Y., Bengio, Y. & Hinton, G. Deep learning. Nature 521, 436–444 (2015).Article 
ADS 
CAS 

Google Scholar 
Yu, L. & Gong, P. Google Earth as a virtual globe tool for Earth science applications at the global scale: Progress and perspectives. Int. J. Remote Sens. 33, 3966–3986 (2012).Article 

Google Scholar 
Guirado, E., Tabik, S., Rivas, M. L., Alcaraz-Segura, D. & Herrera, F. Whale counting in satellite and aerial images with deep learning. Sci. Rep. 9, 14259 (2019).Article 
ADS 

Google Scholar 
Borowicz, A. et al. Aerial-trained deep learning networks for surveying cetaceans from satellite imagery. PLoS ONE 14, 1–15 (2019).Article 

Google Scholar 
Lorencin, I., Anđelić, N., Mrzljak, V. & Car, Z. Marine objects recognition using convolutional neural networks. Nase More 66, 112–119 (2019).Article 

Google Scholar 
Ridge, J. T., Gray, P. C., Windle, A. E. & Johnston, D. W. Deep learning for coastal resource conservation: Automating detection of shellfish reefs. Remote Sens. Ecol. Conserv. 6, 431–440 (2020).Article 

Google Scholar 
Wang, Y. et al. Machine learning-based ship detection and tracking using satellite images for maritime surveillance. J. Ambient Intell. Smart Environ. 13, 361–371 (2021).Article 

Google Scholar 
Han, Q., Yin, Q., Zheng, X. & Chen, Z. Remote sensing image building detection method based on Mask R-CNN. Complex Intell. Syst. https://doi.org/10.1007/s40747-021-00322-z (2021).Article 

Google Scholar 
Girshick, R. Fast R-CNN. In 2015 IEEE International Conference on Computer Vision (ICCV) 1440–1448. https://doi.org/10.1109/ICCV.2015.169 (2015).Ren, S., He, K., Girshick, R. & Sun, J. Faster R-CNN: Towards real-time object detection with region proposal networks. IEEE Trans Pattern Anal. Mach. Intell. 39, 28 (2017).Article 

Google Scholar 
Shelhamer, E., Long, J. & Darrell, T. Fully convolutional networks for semantic segmentation. IEEE Trans. Pattern Anal. Mach. Intell. 39, 3431–3440 (2017).Article 

Google Scholar 
He, K., Gkioxari, G., Dollár, P. & Girshick, R. Mask R-CNN. In Proceedings of the IEEE international Conference on Computer Vision (2017).Arafeh-Dalmau, N. et al. Extreme Marine Heatwaves alter kelp forest community near its equatorward distribution limit. Front. Mar. Sci. 6, 1–18 (2019).Article 
ADS 

Google Scholar 
Nie, X., Duan, M., Ding, H., Hu, B. & Wong, E. K. Attention Mask R-CNN for ship detection and segmentation from remote sensing images. IEEE Access 8, 9325–9334 (2020).Article 

Google Scholar 
Abdulla, W. Mask R-CNN for object detection and instance segmentation on Keras and TensorFlow. GitHub Repository (2017).Fragkopoulou, E. et al. Global biodiversity patterns of marine forests of brown macroalgae. Glob. Ecol. Biogeogr. https://doi.org/10.1111/geb.13450 (2022).Article 

Google Scholar 
Markham, B. L., Storey, J. C., Williams, D. L. & Irons, J. R. Landsat sensor performance: History and current status. IEEE Trans. Geosci. Remote Sens. https://doi.org/10.1109/TGRS.2004.840720 (2004).Article 

Google Scholar 
Gorelick, N. et al. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 202, 18–27 (2017).Article 
ADS 

Google Scholar 
Aghamohamadnia, M. & Abedini, A. A morphology-stitching method to improve Landsat SLC-off images with stripes. Geodesy Geodyn. 5, 27–33 (2014).Article 

Google Scholar 
Houskeeper, H. F. et al. Automated satellite remote sensing of giant kelp at the Falkland Islands (Islas Malvinas). PLoS ONE 17, e0257933 (2022).Article 
CAS 

Google Scholar 
Mantha, K. B. et al. From Fat Droplets to Floating Forests: Cross-Domain Transfer Learning Using a PatchGAN-Based Segmentation Model (2022).Finger, D. J. I., McPherson, M. L., Houskeeper, H. F. & Kudela, R. M. Mapping bull kelp canopy in northern California using Landsat to enable long-term monitoring. Remote Sens. Environ. 254, 112243 (2021).Article 
ADS 

Google Scholar 
Siegel, D. A., Wang, M., Maritorena, S. & Robinson, W. Atmospheric correction of satellite ocean color imagery: The black pixel assumption. Appl. Opt. 39, 3582–3591 (2000).Article 
ADS 
CAS 

Google Scholar 
Loisel, H., Nicolas, J. M., Sciandra, A., Stramski, D. & Poteau, A. Spectral dependency of optical backscattering by marine particles from satellite remote sensing of the global ocean. J. Geophys. Res. Oceans https://doi.org/10.1029/2005JC003367 (2006).Article 

Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).Article 
ADS 
CAS 

Google Scholar 
Dutta, A. & Zisserman, A. The VIA annotation software for images, audio and video. In MM 2019: Proceedings of the 27th ACM International Conference on Multimedia. https://doi.org/10.1145/3343031.3350535 (2019).Pfister, C. A., Berry, H. D. & Mumford, T. The dynamics of Kelp Forests in the Northeast Pacific Ocean and the relationship with environmental drivers. J. Ecol. 106, 1520–1533 (2018).Article 

Google Scholar 
Cavanaugh, K. C., Cavanaugh, K. C., Bell, T. W. & Hockridge, E. G. An automated method for mapping giant kelp canopy dynamics from UAV. Front. Environ. Sci. 8, 587354 (2021).Article 

Google Scholar 
Castorani, M. C. N. et al. Connectivity structures local population dynamics: A long-term empirical test in a large metapopulation system. Ecology 96, 3141–3152 (2015).Article 

Google Scholar 
Irmak, E. Implementation of convolutional neural network approach for COVID-19 disease detection. Physiol. Genom. 52, 590–601 (2020).Article 
CAS 

Google Scholar 
Assis, J., Araújo, M. B. & Serrão, E. A. Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic. Glob. Change Biol. 24, 1365–2486 (2017).
Google Scholar 
Cao, C. et al. An improved faster R-CNN for small object detection. IEEE Access 7, 106838–106846 (2019).Article 

Google Scholar 
Konar, J., Khandelwal, P. & Tripathi, R. Comparison of various learning rate scheduling techniques on convolutional neural network. In 2020 IEEE International Students’ Conference on Electrical, Electronics and Computer Science, SCEECS 2020. https://doi.org/10.1109/SCEECS48394.2020.94 (2020).LeCun, Y., Bottou, L., Bengio, Y. & Haffner, P. Gradient-based learning applied to document recognition. Proc. IEEE 86, 2278–2324 (1998).Article 

Google Scholar 
Johnson, J. W. Automatic nucleus segmentation with mask-RCNN. Adv. Intell. Syst. Comput. 944, 399–407 (2020).
Google Scholar 
Lin, T. Y. et al. Microsoft COCO: Common objects in context. In Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), Vol. 8693 LNCS (2014).McKnight, P. E. & Najab, J. Mann-Whitney U Test. Corsini Encycl. Psychol. https://doi.org/10.1002/9780470479216.corpsy0524 (2010).Article 

Google Scholar 
R Development Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2021).
Google Scholar 
Haklay, M. & Weber, P. OpenStreet map: User-generated street maps. IEEE Pervasive Comput. 7, 12–18 (2008).Article 

Google Scholar 
Wäldchen, J. & Mäder, P. Machine learning for image based species identification. Methods Ecol. Evol. https://doi.org/10.1111/2041-210X.13075 (2018).Article 
MATH 

Google Scholar 
Weinstein, B. G. A computer vision for animal ecology. J. Anim. Ecol. https://doi.org/10.1111/1365-2656.12780 (2018).Article 

Google Scholar 
Chilson, C. et al. Automated detection of bird roosts using NEXRAD radar data and Convolutional Neural Networks. Remote Sens. Ecol. Conserv. 5, 20–32 (2019).Article 

Google Scholar 
O’Gara, S. & McGuinness, K. Comparing data augmentation strategies for deep image classification. Ir. Mach. Vis. Image Process. Conf. https://doi.org/10.21427/148b-ar75 (2019).Article 

Google Scholar 
Li, W., Chen, C., Zhang, M., Li, H. & Du, Q. Data augmentation for hyperspectral image classification with deep CNN. IEEE Geosci. Remote Sens. Lett. 16, 593–597 (2019).Article 
ADS 

Google Scholar 
Bharati, P. & Pramanik, A. Deep learning techniques—R-CNN to Mask R-CNN: A survey. In Computational Intelligence in Pattern Recognition (eds Das, A. K. et al.) 657–668 (Springer, 2020).Chapter 

Google Scholar 
Li, A. S., Chirayath, V., Segal-Rozenhaimer, M., Torres-Perez, J. L. & van den Bergh, J. NASA NeMO-Net’s convolutional neural network: Mapping marine habitats with spectrally heterogeneous remote sensing imagery. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 13, 5115–5133 (2020).Article 
ADS 

Google Scholar 
Hamilton, S. L., Bell, T. W., Watson, J. R., Grorud-Colvert, K. A. & Menge, B. A. Remote sensing: generation of long-term kelp bed data sets for evaluation of impacts of climatic variation. Ecology 101, e03031 (2020).Article 

Google Scholar 
Bell, T. W., Cavanaugh, K. C. & Siegel, D. A. Remote monitoring of giant kelp biomass and physiological condition: An evaluation of the potential for the Hyperspectral Infrared Imager (HyspIRI) mission. Remote Sens. Environ. 167, 218–228 (2015).Article 
ADS 

Google Scholar 
Schroeder, S. B., Dupont, C., Boyer, L., Juanes, F. & Costa, M. Passive remote sensing technology for mapping bull kelp (Nereocystis luetkeana): A review of techniques and regional case study. Glob. Ecol. Conserv. https://doi.org/10.1016/j.gecco.2019.e00683 (2019).Article 

Google Scholar 
Kristollari, V. & Karathanassi, V. Convolutional neural networks for detecting challenging cases in cloud masking using Sentinel-2 imagery. Remote Sens. Geoinf. Environ. https://doi.org/10.1117/12.2571111 (2020).Article 

Google Scholar 
Wilson, M. J. & Oreopoulos, L. Enhancing a simple MODIS cloud mask algorithm for the landsat data continuity mission. IEEE Trans. Geosci. Remote Sens. 51, 723–731 (2013).Article 
ADS 

Google Scholar 
Zhuge, X. Y., Zou, X. & Wang, Y. A fast cloud detection algorithm applicable to monitoring and nowcasting of daytime cloud systems. IEEE Trans. Geosci. Remote Sens. 55, 6111–6119 (2017).Article 
ADS 

Google Scholar 
Lin, T. Y. et al. Feature pyramid networks for object detection. In Proceedings: 30th IEEE Conference on Computer Vision and Pattern Recognition, CVPR 2017 (2017).Jacox, M. G. et al. Impacts of the 2015–2016 El Niño on the California Current System: Early assessment and comparison to past events. Geophys. Res. Lett. https://doi.org/10.1002/2016GL069716 (2016).Article 

Google Scholar 
Chavez, F. P. et al. Biological and chemical consequences of the 1997–1998 El Niño in central California waters. Prog. Oceanogr. https://doi.org/10.1016/S0079-6611(02)00050-2 (2002).Article 

Google Scholar 
Tegner, M. J. & El Dayton, P. K. Niño effects on Southern California kelp forest communities. Adv. Ecol. Res. 17, 243–279 (1987).Article 

Google Scholar 
Bartsch, I. et al. Changes in kelp forest biomass and depth distribution in Kongsfjorden, Svalbard, between 1996–1998 and 2012–2014 reflect Arctic warming. Polar Biol. 39, 2021–2036 (2016).Article 

Google Scholar 
Simonson, E. J., Scheibling, R. E. & Metaxas, A. Kelp in hot water: I. Warming seawater temperature induces weakening and loss of kelp tissue. Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps11438 (2015).Article 

Google Scholar 
Oliver, E. C. J. et al. Projected marine heatwaves in the 21st century and the potential for ecological impact. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00734 (2019).Article 

Google Scholar  More

Monte Carlo analysisSeaweed production costs and net costs of climate benefits were estimated on the basis of outputs of the biophysical and technoeconomic models described below. The associated uncertainties and sensitivities were quantified by repeatedly sampling from uniform distributions of plausible values for each cost and economic parameter (n = 5,000 for each nutrient scenario from the biophysical model, for a total of n = 10,000 simulations; see Supplementary Figs. 14 and 15)47,48,49,50,51,52. Parameter importance across Monte Carlo simulations (Fig. 3 and Supplementary Fig. 9) was determined using decision trees in LightGBM, a gradient-boosting machine learning framework.Biophysical modelG-MACMODS is a nutrient-constrained, biophysical macroalgal growth model with inputs of temperature, nitrogen, light, flow, wave conditions and amount of seeded biomass30,53, that we used to estimate annual seaweed yield per area (either in tons of carbon or tons of dry weight biomass per km2 per year)33,34. In the model, seaweed takes up nitrogen from seawater, and that nitrogen is held in a stored pool before being converted to structural biomass via growth54. Seaweed biomass is then lost via mortality, which includes breakage from variable ocean wave intensity. The conversion from stored nitrogen to biomass is based on the minimum internal nitrogen requirements of macroalgae, and the conversion from biomass to units of carbon is based on an average carbon content of macroalgal dry weight (~30%)55. The model accounts for farming intensity (sub-grid-scale crowding) and employs a conditional harvest scheme, where harvest is optimized on the basis of growth rate and standing biomass33.The G-MACMODS model is parameterized for four types of macroalgae: temperate brown, temperate red, tropical brown and tropical red. These types employed biophysical parameters from genera that represent over 99.5% of present-day farmed macroalgae (Eucheuma, Gracilaria, Kappahycus, Sargassum, Porphyra, Saccharina, Laminaria, Macrocystis)39. Environmental inputs were derived from satellite-based and climatological model output mapped to 1/12-degree global resolution, which resolves continental shelf regions. Nutrient distributions were derived from a 1/10-degree resolution biogeochemical simulation led by the National Center for Atmospheric Research (NCAR) and run in the Community Earth System Model (CESM) framework35.Two nutrient scenarios were simulated with G-MACMODS and evaluated using the technoeconomic model analyses described below: the ‘ambient nutrient’ scenario where seaweed growth was computed using surface nutrient concentrations without depletion or competition, and ‘limited nutrient’ simulations where seaweed growth was limited by an estimation of the nutrient supply to surface waters (computed as the flux of deep-water nitrate through a 100 m depth horizon). For each Monte Carlo simulation in the economic analysis, the technoeconomic model randomly selects either the 5th, 25th, 50th, 75th or 95th percentile G-MACMODS seaweed yield map from a normal distribution to use as the yield map for that simulation. Figures and numbers reported in the main text are based on the ambient-nutrient scenario; results based on the limited-nutrient scenario are shown in Supplementary Figures.Technoeconomic modelAn interactive web tool of the technoeconomic model is available at https://carbonplan.org/research/seaweed-farming.We estimated the net cost of seaweed-related climate benefits by first estimating all costs and emissions related to seaweed farming, up to and including the point of harvest at the farm location, then estimating costs and emissions related to the transportation and processing of harvested seaweed, and finally estimating the market value of seaweed products and either carbon sequestered or GHG emissions avoided.Production costs and emissionsSpatially explicit costs of seaweed production ($ tDW−1) and production-related emissions (tCO2 tDW−1) were calculated on the basis of ranges of capital costs ($ km−2), operating costs (including labour, $ km−2), harvest costs ($ km−2) and transport emissions per distance travelled (tCO2 km−1) in the literature (Table 1, Supplementary Tables 1 and 2); annual seaweed biomass (tDW km−2, for the preferred seaweed type in each grid cell), line spacing and number of harvests (species-dependent) from the biophysical model; as well as datasets of distances to the nearest port (km), ocean depth (m) and significant wave height (m).Capital costs were calculated as:$$c_{cap} = c_{capbase} + left( {c_{capbase} times left( {k_d + k_w} right)} right) + c_{sl}$$
(1)
where ccap is the total annualized capital costs per km2, ccapbase is the annualized capital cost per km2 (for example, cost of buoys, anchors, boats, structural rope) before applying depth and wave impacts, kd and kw are the impacts of depth and waviness on capital cost, respectively, each expressed as a multiplier between 0 and 1 modelled using our Monte Carlo method and applied only to grid cells with depth >500 m and/or significant wave height >3 m, respectively, and csl is the total annual cost of seeded line calculated as:$$c_{sl} = c_{slbase} times p_{sline}$$
(2)
where cslbase is the cost per metre of seeded line, and psline is the total length of line per km2, based on the optimal seaweed type grown in each grid cell.Operating and maintenance costs were calculated as:$$c_{op} = c_{ins} + c_{lic} + c_{lab} + c_{opbase}$$
(3)
where cop is the total annualized operating and maintenance costs per km2, cins is the annual insurance cost per km2, clic is the annual cost of a seaweed aquaculture license per km2, clab is the annual cost of labour excluding harvest labour, and copbase is all other operating and maintenance costs.Harvest costs were calculated as:$$c_{harv} = c_{harvbase} times n_{harv}$$
(4)
where charv is the total annual costs associated with harvesting seaweed per km2, charvbase is the cost per harvest per km2 (including harvest labour but excluding harvest transport), and nharv is the total number of harvests per year.Costs associated with transporting equipment to the farming location were calculated as:$$c_{eqtrans} = c_{transbase} times m_{eq} times d_{port}$$
(5)
where ceqtrans is total annualized cost of transporting equipment, ctransbase is the cost to transport 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons and dport is the ocean distance to the nearest port in km.The total production cost of growing and harvesting seaweed was therefore calculated as:$$c_{prod} = frac{{left( {c_{cap}} right) + left( {c_{op}} right) + left( {c_{harv}} right) + (c_{eqtrans})}}{{s_{dw}}}$$
(6)
where cprod is total annual cost of seaweed production (growth + harvesting), ccap is as calculated in equation (1), cop is as calculated in equation (3), charv is as calculated in equation (4), ceqtrans is as calculated in equation (5) and sdw is the DW of seaweed harvested annually per km2.Emissions associated with transporting equipment to the farming location were calculated as:$$e_{eqtrans} = e_{transbase} times m_{eq} times d_{port}$$
(7)
where eeqtrans is the total annualized CO2 emissions in tons from transporting equipment, etransbase is the CO2 emissions from transporting 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons and dport is the ocean distance to the nearest port in km.Emissions from maintenance trips to/from the seaweed farm were calculated as:$$e_{mnt} = left( {left( {2 times d_{port}} right) times e_{mntbase} times left( {frac{{n_{mnt}}}{{a_{mnt}}}} right)} right) + (e_{mntbase} times d_{mnt})$$
(8)
where emnt is total annual CO2 emissions from farm maintenance, dport is the ocean distance to the nearest port in km, nmnt is the number of maintenance trips per km2 per year, amnt is the area tended to per trip, dmnt is the distance travelled around each km2 for maintenance and emntbase is the CO2 emissions from travelling 1 km on a typical fishing maintenance vessel (for example, a 14 m Marinnor vessel with 2 × 310 hp engines) at an average speed of 9 knots (16.67 km h−1), resulting in maintenance vessel fuel consumption of 0.88 l km−1 (refs. 28,56).Total emissions from growing and harvesting seaweed were therefore calculated as:$$e_{prod} = frac{{(e_{eqtrans}) + (e_{mnt})}}{{s_{dw}}}$$
(9)
where eprod is total annual emissions from seaweed production (growth + harvesting), eeqtrans is as calculated in equation (7), emnt is as calculated in equation (8) and sdw is the DW of seaweed harvested annually per km2.Market value and climate benefits of seaweedFurther transportation and processing costs, economic value and net emissions of either sinking seaweed in the deep ocean for carbon sequestration or converting seaweed into usable products (biofuel, animal feed, pulses, vegetables, fruits, oil crops and cereals) were calculated on the basis of ranges of transport costs ($ tDW−1 km−1), transport emissions (tCO2-eq t−1 km−1), conversion cost ($ tDW−1), conversion emissions (tCO2-eq tDW−1), market value of product ($ tDW−1) and the emissions avoided by product (tCO2-eq tDW−1) in the literature (Table 1). Market value was treated as globally homogeneous and does not vary by region. Emissions avoided by products were determined by comparing estimated emissions related to seaweed production to emissions from non-seaweed products that could potentially be replaced by seaweed (including non-CO2 greenhouse gas emissions from land use)24. Other parameters used are distance to nearest port (km), water depth (m), spatially explicit sequestration fraction (%)57 and distance to optimal sinking location (km; cost-optimized for maximum emissions benefit considering transport emissions combined with spatially explicit sequestration fraction; see ‘Distance to sinking point calculation’ below). Each Monte Carlo simulation calculated the cost of both CDR via sinking seaweed and GHG emissions mitigation via seaweed products.For seaweed CDR, after the seaweed is harvested, it can either be sunk in the same location that it was grown, or be transported to a more economically favourable sinking location where more of the seaweed carbon would remain sequestered for 100 yr (see ‘Distance to optimal sinking point’ below). Immediately post-harvest, the seaweed still contains a large amount of water, requiring a conversion from dry mass to wet mass for subsequent calculations33:$$s_{ww} = frac{{s_{dw}}}{{0.1}}$$
(10)
where sww is the annual wet weight of seaweed harvested per km2 and sdw is the annual DW of seaweed harvested per km2.The cost to transport harvested seaweed to the optimal sinking location was calculated as:$$c_{swtsink} = c_{transbase} times d_{sink} times s_{ww}$$
(11)
where cswtsink is the total annual cost to transport harvested seaweed to the optimal sinking location, ctransbase is the cost to transport 1 ton of material 1 km on a barge, dsink is the distance in km to the economically optimized sinking location and sww is the annually harvested seaweed wet weight in t km−2 as in equation (10).The costs associated with transporting replacement equipment (for example, lines, buoys,anchors) to the farming location and hauling back used equipment at the end of its assumed lifetime (1 yr for seeded line, 5–20 yr for capital equipment by equipment type) in the sinking CDR pathway were calculated as:$$c_{eqtsink} = left( {c_{transbase} times left( {2 times d_{sink}} right) times m_{eq}} right) + (c_{transbase} times d_{port} times m_{eq})$$
(12)
where ceqtsink is the total annualized cost to transport both used and replacement equipment, ctransbase is the cost to transport 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons, dsink is the distance in km to the economically optimized sinking location and dport is the ocean distance to the nearest port in km. We assumed that the harvesting barge travels from the farming location directly to the optimal sinking location with harvested seaweed and replaced (used) equipment in tow (including used seeded line and annualized mass of used capital equipment), sinks the harvested seaweed, returns to the farm location and then returns to the nearest port (see Supplementary Fig. 16). These calculations assumed the shortest sea-route distance (see Distance to optimal sinking point).The total value of seaweed that is sunk for CDR was therefore calculated as:$$v_{sink} = frac{{left( {v_{cprice} – left( {c_{swtsink} + c_{eqtsink}} right)} right)}}{{s_{dw}}}$$
(13)
where vsink is the total value (cost, if negative) of seaweed farmed for CDR in $ tDW−1, vcprice is a theoretical carbon price, cswtsink is as calculated in equation (11), ceqtsink is as calculated in equation (12) and sdw is the annually harvested seaweed DW in t km−2. We did not assume any carbon price in our Monte Carlo simulations (vcprice is equal to zero), making vsink negative and thus representing a net cost.To calculate net carbon impacts, our model included uncertainty in the efficiency of using the growth and subsequent deep-sea deposition of seaweed as a CDR method. The uncertainty is expected to include the effects of reduced phytoplankton growth from nutrient competition, the relationship between air–sea gas exchange and overturning circulation (hereafter collectively referred to as the ‘atmospheric removal fraction’) and the fraction of deposited seaweed carbon that remains sequestered for at least 100 yr. The total amount of atmospheric CO2 removed by sinking seaweed was calculated as:$$e_{seqsink} = k_{atm} times k_{fseq} times frac{{tC}}{{tDW}} times frac{{tCO_2}}{{tC}}$$
(14)
where eseqsink is net atmospheric CO2 sequestered annually in t km−2, katm is the atmospheric removal fraction and kfseq is the spatially explicit fraction of sunk seaweed carbon that remains sequestered for at least 100 yr (see ref. 57).The emissions from transporting harvested seaweed to the optimal sinking location were calculated as:$$e_{swtsink} = e_{transbase} times d_{sink} times s_{ww}$$
(15)
where eswtsink is the total annual CO2 emissions from transporting harvested seaweed to the optimal sinking location in tCO2 km−2, etransbase is the CO2 emissions (tons) from transporting 1 ton of material 1 km on a barge (tCO2 per t-km), dsink is the distance in km to the economically optimized sinking location and sww is the annually harvested seaweed wet weight in t km−2 as in equation (10). Since the unit for etransbase is tCO2 per t-km, the emissions from transporting seaweed to the optimal sinking location are equal to (e_{mathrm{transbase}} times d_{mathrm{sink}} times s_{mathrm{ww}}), and the emissions from transporting seaweed from the optimal sinking location back to the farm are equal to 0 (since the seaweed has already been deposited, the seaweed mass to transport is now 0). Note that this does not yet include transport emissions from transport of equipment post-seaweed-deposition (see equation 16 below and Supplementary Fig. 16).The emissions associated with transporting replacement equipment (for example, lines, buoys, anchors) to the farming location and hauling back used equipment at the end of its assumed lifetime (1 yr for seeded line, 5–20 yr for capital equipment by equipment type)28,41 in the sinking CDR pathway were calculated as:$$e_{eqtsink} = left( {e_{transbase} times left( {2 times d_{sink}} right) times m_{eq}} right) + (e_{transbase} times d_{port} times m_{eq})$$
(16)
where eeqtsink is the total annualized CO2 emissions in tons from transporting both used and replacement equipment, etransbase is the CO2 emissions from transporting 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons, dsink is the distance in km to the economically optimized sinking location and dport is the ocean distance to the nearest port in km. We assumed that the harvesting barge travels from the farming location directly to the optimal sinking location with harvested seaweed and replaced (used) equipment in tow (including used seeded line and annualized mass of used capital equipment), sinks the harvested seaweed, returns to the farm location and then returns to the nearest port. These calculations assumed the shortest sea-route distance (see Distance to optimal sinking point).Net CO2 emissions removed from the atmosphere by sinking seaweed were thus calculated as:$$e_{remsink} = frac{{left( {e_{seqsink} – left( {e_{swtsink} + e_{eqtsink}} right)} right)}}{{s_{dw}}}$$
(17)
where eremsink is the net atmospheric CO2 removed per ton of seaweed DW, eseqsink is as calculated in equation (14), eswtsink is as calculated in equation (15), eeqtsink is as calculated in equation (16) and sdw is the annually harvested seaweed DW in t km−2.Net cost of climate benefitsSinkingTo calculate the total net cost and emissions from the production, harvesting and transport of seaweed for CDR, we combined the cost and emissions from the sinking-pathway cost and value modules. The total net cost of seaweed CDR per DW ton of seaweed was calculated as:$$c_{sinknet} = c_{prod} – v_{sink}$$
(18)
where csinknet is the total net cost of seaweed for CDR per DW ton harvested, cprod is the net production cost per DW ton as calculated in equation (6) and vsink is the net value (or cost, if negative) per ton seaweed DW as calculated in equation (13).The total net CO2 emissions removed per DW ton of seaweed were calculated as:$$e_{sinknet} = e_{remsink} – e_{prod}$$
(19)
where esinknet is the total net atmospheric CO2 removed per DW ton of seaweed harvested annually (tCO2 tDW−1 yr−1), eremsink is the net atmospheric CO2 removed via seaweed sinking annually as calculated in equation (17) and eprod is the net CO2 emitted from production and harvesting of seaweed annually as calculated in equation (9). For each Monte Carlo simulation, locations where esinknet is negative (that is, net emissions rather than net removal) were not included in subsequent calculations since they would not be contributing to CDR in that location under the given scenario. Note that these net emissions cases only occur in areas far from port in specific high-emissions scenarios. Even in such cases, most areas still contribute to CO2 removal (negative emissions), hence costs from locations with net removal were included.Total net cost was then divided by total net emissions to get a final value for cost per ton of atmospheric CO2 removed:$$c_{pertonsink} = frac{{c_{sinknet}}}{{e_{sinknet}}}$$
(20)
where cpertonsink is the total net cost per ton of atmospheric CO2 removed via seaweed sinking ($ per tCO2 removed), csinknet is total net cost per ton seaweed DW harvested as calculated in equation (18) ($ tDW−1) and esinknet is the total net atmospheric CO2 removed per ton seaweed DW harvested as calculated in equation (19) (tCO2 tDW−1).GHG emissions mitigationInstead of sinking seaweed for CDR, seaweed can be used to make products (including but not limited to food, animal feed and biofuels). Replacing convention products with seaweed-based products can result in ‘avoided emissions’ if the emissions from growing, harvesting, transporting and converting seaweed into products are less than the total greenhouse gas emissions (including non-CO2 GHGs) embodied in conventional products that seaweed-based products replace.When seaweed is used to make products, we assumed it is transported back to the nearest port immediately after being harvested. The annualized cost to transport the harvested seaweed and replacement equipment (for example, lines, buoys, anchors) was calculated as:$$c_{transprod} = frac{{left( {c_{transbase} times d_{port} times left( {s_{ww} + m_{eq}} right)} right)}}{{s_{dw}}}$$
(21)
where ctransprod is the annualized cost per ton seaweed DW to transport seaweed and equipment back to port from the farm location, ctransbase is the cost to transport 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons, dport is the ocean distance to the nearest port in km, sww is the annual wet weight of seaweed harvested per km2 as calculated in equation (10) and sdw is the annual DW of seaweed harvested per km2.The total value of seaweed that is used for seaweed-based products was calculated as:$$v_{product} = v_{mkt} – left( {c_{transprod} + c_{conv}} right)$$
(22)
where vproduct is the total value (cost, if negative) of seaweed used for products ($ tDW−1), vmkt is how much each ton of seaweed would sell for, given the current market price of conventional products that seaweed-based products replace ($ tDW−1), ctransprod is as calculated in equation (21) and cconv is the cost to convert each ton of seaweed to a usable product ($ tDW−1).The annualized CO2 emissions from transporting harvested seaweed and equipment back to port were calculated as:$$e_{transprod} = frac{{left( {e_{transbase} times d_{port} times left( {s_{ww} + m_{eq}} right)} right)}}{{s_{dw}}}$$
(23)
where etransprod is the annualized CO2 emissions per ton seaweed DW to transport seaweed and equipment back to port from the farm location, etransbase is the CO2 emissions from transporting 1 ton of material 1 km on a barge, meq is the annualized equipment mass in tons, dport is the ocean distance to the nearest port in km, sww is the annual wet weight of seaweed harvested per km2 as calculated in equation (10) and sdw is the annual DW of seaweed harvested per km2.Total emissions avoided by each ton of harvested seaweed DW were calculated as:$$e_{avprod} = e_{subprod} – left( {e_{transprod} + e_{conv}} right)$$
(24)
where eavprod is total CO2-eq emissions avoided per ton of seaweed DW per year (including non-CO2 GHGs using a GWP time period of 100 yr), esubprod is the annual CO2-eq emissions avoided per ton seaweed DW by replacing a conventional product with a seaweed-based product, etransprod is as calculated in equation (23) and econv is the annual CO2 emissions per ton seaweed DW from converting seaweed into usable products. esubprod was calculated by converting seaweed DW to caloric content58 for food/feed and comparing emissions intensity per kcal to agricultural products24, or by converting seaweed DW into equivalent biofuel content with a yield of 0.25 tons biofuel per ton DW59 and dividing the CO2 emissions per ton fossil fuel by the seaweed biofuel yield.To calculate the total net cost and emissions from the production, harvesting, transport and conversion of seaweed for products, we combined the cost and emissions from the product-pathway cost and value modules. The total net cost of seaweed for products per ton DW was calculated as:$$c_{prodnet} = c_{prod} – v_{product}$$
(25)
where cprodnet is the total net cost per ton DW of seaweed harvested for use in products, cprod is the net production cost per ton DW as calculated in equation (6) and vproduct is the net value (or cost, if negative) per ton DW as calculated in equation (22).The total net CO2-eq emissions avoided per ton DW of seaweed used in products were calculated as:$$e_{prodnet} = e_{avprod} – e_{prod}$$
(26)
where eprodnet is the total net CO2-eq emissions avoided per ton DW of seaweed harvested annually (tCO2 tDW−1 yr−1), eavprod is the net CO2-eq emissions avoided by seaweed products annually as calculated in equation (24) and eprod is the net CO2 emitted from production and harvesting of seaweed annually as calculated in equation (9). For each Monte Carlo simulation, locations where eprodnet is negative (that is, net emissions rather than net emissions avoided) were not included in subsequent calculations since they would not be avoiding any emissions in that scenario.Total net cost was then divided by total net emissions avoided to get a final value for cost per ton of CO2-eq emissions avoided:$$c_{pertonprod} = frac{{c_{prodnet}}}{{e_{prodnet}}}$$
(27)
where cpertonprod is the total net cost per ton of CO2-eq emissions avoided by seaweed products ($ per tCO2-eq avoided), cprodnet is total net cost per ton seaweed DW harvested for products as calculated in equation (25) ($ tDW−1) and eprodnet is total net CO2-eq emissions avoided per ton seaweed DW harvested for products as calculated in equation (26) (tCO2 tDW−1).Parameter ranges for Monte Carlo simulationsFor technoeconomic parameters with two or more literature values (see Supplementary Table 1), we assumed that the maximum literature value reflected the 95th percentile and the minimum literature value represented the 5th percentile of potential costs or emissions. For parameters with only one literature value, we added ±50% to the literature value to represent greater uncertainty within the modelled parameter range. Values at each end of parameter ranges were then rounded before Monte Carlo simulations as follows: capital costs, operating costs and harvest costs to the nearest $10,000 km−2, labour costs and insurance costs to the nearest $1,000 km−2, line costs to the nearest $0.05 m−1, transport costs to the nearest $0.05 t−1 km−1, transport emissions to the nearest 0.000005 tCO2 t−1 km−1, maintenance transport emissions to the nearest 0.0005 tCO2 km−1, product-avoided emissions to the nearest 0.1 tCO2-eq tDW−1, conversion cost down to the nearest $10 tDW−1 on the low end of the range and up to the nearest $10 tDW−1 on the high end of the range, and conversion emissions to the nearest 0.01 tCO2 tDW−1.We extended the minimum range values of capital costs to $10,000 km−2 and transport emissions to 0 to reflect potential future innovations, such as autonomous floating farm setups that would lower capital costs and net-zero emissions boats that would result in 0 transport emissions. To calculate the minimum value of $10,000 km−2 for a potential autonomous floating farm, we assumed that the bulk of capital costs for such a system would be from structural lines and flotation devices, and we therefore used the annualized structural line (system rope) and buoy costs from ref. 41 rounded down to the nearest $5,000 km−2. The full ranges used for our Monte Carlo simulations and associated literature values are shown in Supplementary Table 1.Distance to optimal sinking pointDistance to the optimal sinking point was calculated using a weighted distance transform (path-finding algorithm, modified from code in ref. 60) that finds the shortest ocean distance from each seaweed growth pixel to the location at which the net CO2 removed is maximized (including impacts of both increased sequestration fraction and transport emissions for different potential sinking locations) and the net cost is minimized. This is not necessarily the location in which the seaweed was grown, since the fraction of sunk carbon that remains sequestered for 100 yr is spatially heterogeneous (see ref. 57). For each ocean grid cell, we determined the cost-optimal sinking point by iteratively calculating equations (11–20) and assigning dsink as the distance calculated by weighted distance transform to each potential sequestration fraction (0.01–1.00) in increments of 0.01. Except for transport emissions, the economic parameter values used for these calculations were the averages of unrounded literature value ranges; we assumed that the maximum literature value reflected the 95th percentile and the minimum literature value represented the 5th percentile of potential costs or emissions, or for parameters with only one literature value, we added ±50% to the literature value to represent greater uncertainty within the modelled parameter range. For transport and maintenance transport emissions, we extended the minimum values of the literature ranges to zero to reflect potential net-zero emissions transport options and used the mean values of the resulting ranges. The dsink that resulted in minimum net cost per ton CO2 for each ocean grid cell was saved as the final dsink map, and the associated sequestration fraction value that the seaweed is transported to via dsink was assigned to the original cell where the seaweed was farmed and harvested (Supplementary Fig. 19). If the cost-optimal location to sink using this method was the same cell where the seaweed was harvested, then dsink was 0 km and the sequestration fraction was not modified from its original value (Supplementary Fig. 18).Comparison of gigaton-scale sequestration area to previous estimatesPrevious related work estimating the ocean area suitable for macroalgae cultivation13 and/or the area that might be required to reach gigaton-scale carbon removal via macroalgae cultivation13,19,36 has yielded a wide range of results, primarily due to differences in modelling methods. For example, Gao et al. (2022)36 estimate that 1.15 million km2 would be required to sequester 1 GtCO2 annually when considering carbon lost from seaweed biomass/sequestered as particulate organic carbon (POC) and refractory dissolved organic carbon (rDOC), and assume that the harvested seaweed is sold as food such that the carbon in the harvested seaweed is not sequestered. The area (0.31 million km2) required to sequester 1 GtCO2 in our study assumes that all harvested seaweed is sunk to the deep ocean to sequester carbon.Additionally, Wu et al.19 estimates that roughly 12 GtCO2 could be sequestered annually via macroalgae cultivation in approximately 20% of the world ocean area (that is, 1.67% ocean area per GtCO2), which is a much larger area per GtCO2 than our estimate of 0.085% ocean area. This notable difference arises for several reasons (including differences in yields, which in Wu et al. are around 500 tDW yr−1 in the highest-yield areas, whereas yields in our cheapest sequestration areas from G-MACMODS average 3,400 tDW km−2 yr−1) that arise from differences in model methodology. First, Wu et al. model temperate brown seaweeds, while our study considers different seaweed types, many of which have higher growth rates, and uses the most productive seaweed type for each ocean grid cell. The G-MACMODS seaweed growth model we use also has a highly optimized harvest schedule, includes luxury nutrient uptake (a key feature of macroalgal nutrient physiology) and does not directly model competition with phytoplankton during seaweed growth. Finally, tropical red seaweeds (the seaweed type in our cheapest sequestration areas) grow year-round, while others, such as the temperate brown seaweeds modelled by Wu et al., only grow seasonally. These differences all contribute to higher productivity in our model, leading to a smaller area required for gigaton-scale CO2 sequestration compared with Wu et al.Conversely, the ocean areas we model for seaweed-based CO2 sequestration or GHG emissions avoided are much larger than the 48 million km2 that Froehlich et al.13 estimate to be suitable for macroalgae farming globally. Although our maps show productivity and costs everywhere, the purpose of our modelling was to evaluate where different types of seaweed grow best and how production costs and product values vary over space, to highlight the lowest-cost areas (which are often the highest-producing areas) under various technoeconomic assumptions.Comparison of seaweed production costs to previous estimatesAlthough there are not many estimates of seaweed production costs in the scientific literature, our estimates for the lowest-cost 1% area of the ocean ($190–$2,790 tDW−1) are broadly consistent with previously published results: seaweed production costs reported in the literature range from $120 to $1,710 tDW−1 (refs. 40,41,61,62), but are highly dependent on assumed seaweed yields. For example, Camus et al.41 calculate a cost of $870 tDW−1 assuming a minimum yield of 12.4 kgDW m−1 of cultivation line (equivalent to 8.3 kgDW m−2 using 1.5 m spacing between lines). Using the economic values from Camus et al. but with our estimates of average yield for the cheapest 1% production cost areas (2.6 kgDW m−2) gives a much higher average cost of $2,730 tDW−1. Contrarily, van den Burg et al.40 calculate a cost of $1,710 tDW−1 using a yield of 20 tDW ha−1 (that is, 2.0 kg m−2). Instead assuming the average yield to be that from our lowest-cost areas (that is, 2.6 kgDW m−2 or 26 tDW ha−1) would decrease the cost estimated by van den Burg et al. (2016) to $1,290 tDW−1. Most recently, Capron et al.62 calculate an optimistic scenario cost of $120 tDW−1 on the basis of an estimated yield of 120 tDW ha−1 (12 kg m−2; over 4.5 times higher than the average yield in our lowest-cost areas). Again, instead assuming the average yield to be that in our lowest-cost areas would raise Capron et al.’s production cost to $540 tDW−1 (between the $190–$880 tDW−1 minimum to median production costs in the cheapest 1% areas from our model; Fig. 1a,b).Data sourcesSeaweed biomass harvestedWe used spatially explicit data for seaweed harvested globally under both ambient and limited-nutrient scenarios from the G-MACMODS seaweed growth model33.Fraction of deposited carbon sequestered for 100 yrWe used data from ref. 57 interpolated to our 1/12-degree grid resolution.Distance to the nearest portWe used the Distance from Port V1 dataset from Global Fishing Watch (https://globalfishingwatch.org/data-download/datasets/public-distance-from-port-v1) interpolated to our 1/12-degree grid resolution.Significant wave heightWe used data for annually averaged significant wave height from the European Center for Medium-range Weather Forecasts (ECMWF) interpolated to our 1/12-degree grid resolution.Ocean depthWe used data from the General Bathymetric Chart of the Oceans (GEBCO).Shipping lanesWe used data of Automatic Identification System (AIS) signal count per ocean grid cell, interpolated to our 1/12-degree grid resolution. We defined a major shipping lane grid cell as any cell with >2.25 × 108 AIS signals, a threshold that encompasses most major trans-Pacific and trans-Atlantic shipping lanes as well as major shipping lanes in the Indian Ocean, the North Sea, and coastal routes worldwide.Marine protected areas (MPAs)We used data from the World Database on Protected Areas (WDPA) and defined an MPA as any ocean MPA >20 km2.Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

Azam, F. et al. The ecological role of water-column microbes in the sea. Marine Ecol. Prog. Ser. 10, 257–263 (1983).Fuhrman, J. A. & Caron D. A. in Manual of Environmental Microbiology (eds Yates, M. V. et al.) 4.2.2–4.2.2.-34 (ASM Press, 2016).Gasol, J. M. & Kirchman, D. L. Microbial Ecology of the Oceans (John Wiley & Sons, 2018).Fuhrman, J. A. et al. A latitudinal diversity gradient in planktonic marine bacteria. Proc. Natl Acad. Sci. 105, 7774–7778 (2008).Article 
ADS 
CAS 

Google Scholar 
Gilbert, J. A. et al. The seasonal structure of microbial communities in the Western English Channel. Environ. Microbiol. 11, 3132–3139 (2009).Article 
CAS 

Google Scholar 
Gilbert, J. A. et al. Defining seasonal marine microbial community dynamics. ISME J. 6, 298–308 (2012).Article 
CAS 

Google Scholar 
Hatosy, S. M. et al. Beta diversity of marine bacteria depends on temporal scale. Ecology 94, 1898–1904 (2013).Article 

Google Scholar 
Ward, C. S. et al. Annual community patterns are driven by seasonal switching between closely related marine bacteria. ISME J. 11, 1412–1422 (2017).Article 

Google Scholar 
Fuhrman, J. A. et al. Annually reoccurring bacterial communities are predictable from ocean conditions. Proc. Natl Acad. Sci. 103, 13104–13109 (2006).Article 
ADS 
CAS 

Google Scholar 
Gonzalez, J. M., Sherr, E. B. & Sherr, B. F. Size-selective grazing on bacteria by natural assemblages of estuarine flagellates and ciliates. Appl. Environ. Microbiol. 56, 583–589 (1990).Article 
ADS 
CAS 

Google Scholar 
Guixa-Boixereu, N., Vaque, D., Gasol, J. M. & Pedros-Alio, C. Distribution of viruses and their potential effect on bacterioplankton in an oligotrophic marine system. Aquat. Microb. Ecol. 19, 205–213 (1999).Article 

Google Scholar 
Šimek, K. et al. Shifts in bacterial community composition associated with different microzooplankton size fractions in a eutrophic reservoir. Limnol. Oceanogr. 44, 1634–1644 (1999).Article 
ADS 

Google Scholar 
Hewson, I., Vargo, G. & Fuhrman, J. Bacterial diversity in shallow oligotrophic marine benthos and overlying waters: effects of virus infection, containment, and nutrient enrichment. Microb. Ecol. 46, 322–336 (2003).Article 
CAS 

Google Scholar 
Schwalbach, M. S., Hewson, I. & Fuhrman, J. A. Viral effects on bacterial community composition in marine plankton microcosms. Aquat. Microb. Ecol. 34, 117–127 (2004).Article 

Google Scholar 
Winter, C., Smit, A., Herndl, G. J. & Weinbauer, M. G. Linking bacterial richness with viral abundance and prokaryotic activity. Limnol. Oceanogr. 50, 968–977 (2005).Article 
ADS 

Google Scholar 
Chow, C.-E. T., Kim, D. Y., Sachdeva, R., Caron, D. A. & Fuhrman, J. A. Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists. ISME J. 8, 816–829 (2014).Article 
CAS 

Google Scholar 
Suzuki, S. et al. Comparison of community structures between particle-associated and free-living prokaryotes in tropical and subtropical Pacific Ocean surface waters. J. Oceanogr. 73, 383–395 (2017).Article 
CAS 

Google Scholar 
Milici, M. et al. Diversity and community composition of particle‐associated and free‐living bacteria in mesopelagic and bathypelagic Southern Ocean water masses: evidence of dispersal limitation in the Bransfield Strait. Limnol. Oceanogr. 62, 1080–1095 (2017).Article 
ADS 

Google Scholar 
D’ambrosio, L., Ziervogel, K., MacGregor, B., Teske, A. & Arnosti, C. Composition and enzymatic function of particle-associated and free-living bacteria: a coastal/offshore comparison. ISME J. 8, 2167–2179 (2014).Article 

Google Scholar 
Rieck, A., Herlemann, D. P., Jürgens, K. & Grossart, H.-P. Particle-associated differ from free-living bacteria in surface waters of the Baltic Sea. Front. Microbiol. 6, 1297 (2015).Article 

Google Scholar 
Yung, C.-M., Ward, C. S., Davis, K. M., Johnson, Z. I. & Hunt, D. E. Insensitivity of diverse and temporally variable particle-associated microbial communities to bulk seawater environmental parameters. Appl. Environ. Microbiol. 82, 3431–3437 (2016).Article 
ADS 
CAS 

Google Scholar 
Buchan, A., LeCleir, G. R., Gulvik, C. A. & González, J. M. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698 (2014).Article 
CAS 

Google Scholar 
Duret, M. T., Lampitt, R. S. & Lam, P. Prokaryotic niche partitioning between suspended and sinking marine particles. Environ. Microbiol. Rep. 11, 386–400 (2019).Article 
CAS 

Google Scholar 
Crespo, B. G., Pommier, T., Fernández‐Gómez, B. & Pedrós‐Alió, C. Taxonomic composition of the particle‐attached and free‐living bacterial assemblages in the Northwest Mediterranean Sea analyzed by pyrosequencing of the 16S rRNA. Microbiologyopen 2, 541–552 (2013).Article 
CAS 

Google Scholar 
Mestre, M., Borrull, E., Sala, M. & Gasol, J. M. Patterns of bacterial diversity in the marine planktonic particulate matter continuum. ISME J. 11, 999–1010 (2017).Yeh, Y. C. et al. Comprehensive single‐PCR 16S and 18S rRNA community analysis validated with mock communities, and estimation of sequencing bias against 18S. Environ. Microbiol. 23, 3240–3250 (2021).Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).Article 
CAS 

Google Scholar 
Needham, D. M. et al. Dynamics and interactions of highly resolved marine plankton via automated high-frequency sampling. ISME J. 12, 2417 (2018).Article 
CAS 

Google Scholar 
McNichol, J., Berube, P. M., Biller, S. J. & Fuhrman, J. A. Evaluating and improving small subunit rRNA PCR primer coverage for bacteria, archaea, and eukaryotes using metagenomes from global ocean surveys. Msystems 6, e00565–00521 (2021).Article 
CAS 

Google Scholar 
Chow, C. E. T. & Fuhrman, J. A. Seasonality and monthly dynamics of marine myovirus communities. Environ. Microbiol. 14, 2171–2183 (2012).Article 

Google Scholar 
Filée, J., Tétart, F., Suttle, C. A. & Krisch, H. Marine T4-type bacteriophages, a ubiquitous component of the dark matter of the biosphere. Proc. Natl Acad. Sci. 102, 12471–12476 (2005).Article 
ADS 

Google Scholar 
Pagarete, A. et al. Strong seasonality and interannual recurrence in marine myovirus communities. Appl. Environ. Microbiol. 79, 6253–6259 (2013).Article 
ADS 
CAS 

Google Scholar 
Comeau, A. M. & Krisch, H. M. The capsid of the T4 phage superfamily: the evolution, diversity, and structure of some of the most prevalent proteins in the biosphere. Mol. Biol. Evolution 25, 1321–1332 (2008).Article 
CAS 

Google Scholar 
Needham, D. M. et al. Short-term observations of marine bacterial and viral communities: patterns, connections and resilience. ISME J. 7, 1274–1285 (2013).Article 
CAS 

Google Scholar 
Needham, D. M., Sachdeva, R. & Fuhrman, J. A. Ecological dynamics and co-occurrence among marine phytoplankton, bacteria and myoviruses shows microdiversity matters. ISME J. 11, 1614–1629 (2017).Article 

Google Scholar 
Ahlgren, N. A., Perelman, J. N., Yeh, Y. C. & Fuhrman, J. A. Multi‐year dynamics of fine‐scale marine cyanobacterial populations are more strongly explained by phage interactions than abiotic, bottom‐up factors. Environ. Microbiol. 21, 2948–2963 (2019).Article 
CAS 

Google Scholar 
Ren, J., Ahlgren, N. A., Lu, Y. Y., Fuhrman, J. A. & Sun, F. VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data. Microbiome 5, 1–20 (2017).Article 

Google Scholar 
Roux, S., Enault, F., Hurwitz, B. L. & Sullivan, M. B. VirSorter: mining viral signal from microbial genomic data. PeerJ 3, e985 (2015).Article 

Google Scholar 
Ignacio-Espinoza, J. C., Ahlgren, N. A. & Fuhrman, J. A. Long-term stability and Red Queen-like strain dynamics in marine viruses. Nat. Microbiol. 5, 265–271 (2020).Article 
CAS 

Google Scholar 
Brum, J. R. et al. Patterns and ecological drivers of ocean viral communities. Science 348, (2015).Brown, M. V. et al. Global biogeography of SAR11 marine bacteria. Mol. Syst. Biol. 8, 595 (2012).Article 
ADS 

Google Scholar 
Johnson, Z. I. et al. Niche partitioning among Prochlorococcus ecotypes along ocean-scale environmental gradients. Science 311, 1737–1740 (2006).Article 
ADS 
CAS 

Google Scholar 
Zwirglmaier, K. et al. Global phylogeography of marine Synechococcus and Prochlorococcus reveals a distinct partitioning of lineages among oceanic biomes. Environ. Microbiol. 10, 147–161 (2008).
Google Scholar 
Martiny, A. C., Tai, A. P., Veneziano, D., Primeau, F. & Chisholm, S. W. Taxonomic resolution, ecotypes and the biogeography of Prochlorococcus. Environ. Microbiol. 11, 823–832 (2009).Article 

Google Scholar 
Bond, N. A., Cronin, M. F., Freeland, H. & Mantua, N. Causes and impacts of the 2014 warm anomaly in the NE Pacific. Geophys. Res. Lett. 42, 3414–3420 (2015).Article 
ADS 

Google Scholar 
Di Lorenzo, E. & Mantua, N. Multi-year persistence of the 2014/15 North Pacific marine heatwave. Nat. Clim. Change 6, 1042–1047 (2016).Article 
ADS 

Google Scholar 
Traving, S. J. et al. Prokaryotic responses to a warm temperature anomaly in northeast subarctic Pacific waters. Commun. Biol. 4, 1–12 (2021).Article 

Google Scholar 
Peña, M. A., Nemcek, N. & Robert, M. Phytoplankton responses to the 2014–2016 warming anomaly in the northeast subarctic Pacific Ocean. Limnol. Oceanogr. 64, 515–525 (2019).Article 
ADS 

Google Scholar 
Yang, B., Emerson, S. R. & Peña, M. A. The effect of the 2013–2016 high temperature anomaly in the subarctic Northeast Pacific (the “Blob”) on net community production. Biogeosciences 15, 6747–6759 (2018).Article 
ADS 
CAS 

Google Scholar 
Cavole, L. M. et al. Biological impacts of the 2013–2015 warm-water anomaly in the Northeast Pacific: winners, losers, and the future. Oceanography 29, 273–285 (2016).Article 

Google Scholar 
Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).Article 
CAS 

Google Scholar 
Needham, D. M. & Fuhrman, J. A. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat. Microbiol. 1, 16005 (2016).Article 
CAS 

Google Scholar 
Grossart, H. P., Levold, F., Allgaier, M., Simon, M. & Brinkhoff, T. Marine diatom species harbour distinct bacterial communities. Environ. Microbiol. 7, 860–873 (2005).Article 
CAS 

Google Scholar 
Teeling, H. et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611 (2012).Article 
ADS 
CAS 

Google Scholar 
Chafee, M. et al. Recurrent patterns of microdiversity in a temperate coastal marine environment. ISME J. 12, 237–252 (2018).Article 

Google Scholar 
Teeling, H. et al. Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. elife 5, e11888 (2016).Article 

Google Scholar 
Unfried, F. et al. Adaptive mechanisms that provide competitive advantages to marine bacteroidetes during microalgal blooms. ISME J. 12, 2894–2906 (2018).Article 
CAS 

Google Scholar 
Francis, T. B. et al. Changing expression patterns of TonB-dependent transporters suggest shifts in polysaccharide consumption over the course of a spring phytoplankton bloom. ISME J. 15, 2336–2350 (2021).Thingstad, T. F. & Lignell, R. Theoretical models for the control of bacterial growth rate, abundance, diversity and carbon demand. Aquat. Microb. Ecol. 13, 19–27 (1997).Article 

Google Scholar 
Thingstad, T. F., Våge, S., Storesund, J. E., Sandaa, R.-A. & Giske, J. A theoretical analysis of how strain-specific viruses can control microbial species diversity. Proc. Natl Acad. Sci. 111, 7813–7818 (2014).Article 
ADS 
CAS 

Google Scholar 
Thingstad, T. F., Pree, B., Giske, J. & Våge, S. What difference does it make if viruses are strain-, rather than species-specific? Front. Microbiol. 6, 320 (2015).Article 

Google Scholar 
Prokopowich, C. D., Gregory, T. R. & Crease, T. J. The correlation between rDNA copy number and genome size in eukaryotes. Genome 46, 48–50 (2003).Article 
CAS 

Google Scholar 
Zhu, F., Massana, R., Not, F., Marie, D. & Vaulot, D. Mapping of picoeucaryotes in marine ecosystems with quantitative PCR of the 18S rRNA gene. FEMS Microbiol. Ecol. 52, 79–92 (2005).Article 
CAS 

Google Scholar 
Sintes, E. & Del Giorgio, P. A. Feedbacks between protistan single-cell activity and bacterial physiological structure reinforce the predator/prey link in microbial foodwebs. Front. Microbiol. 5, 453 (2014).Article 

Google Scholar 
Del Giorgio, P. A. et al. Bacterioplankton community structure: protists control net production and the proportion of active bacteria in a coastal marine community. Limnol. Oceanogr. 41, 1169–1179 (1996).Article 
ADS 

Google Scholar 
Andersson, A., Larsson, U. & Hagström, Å. Size-selective grazing by a microflagellate on pelagic bacteria. Marine Ecol. Prog. Ser. 33, 51–57 (1986).Pernthaler, J. Predation on prokaryotes in the water column and its ecological implications. Nat. Rev. Microbiol. 3, 537–546 (2005).Article 
CAS 

Google Scholar 
Baltar, F. et al. Marine bacterial community structure resilience to changes in protist predation under phytoplankton bloom conditions. ISME J. 10, 568–581 (2016).Article 

Google Scholar 
Suzuki, M. T. Effect of protistan bacterivory on coastal bacterioplankton diversity. Aquat. Microb. Ecol. 20, 261–272 (1999).Article 

Google Scholar 
Yokokawa, T. & Nagata, T. Growth and grazing mortality rates of phylogenetic groups of bacterioplankton in coastal marine environments. Appl. Environ. Microbiol. 71, 6799–6807 (2005).Article 
ADS 
CAS 

Google Scholar 
Eren, A. M. et al. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol. Evolution 4, 1111–1119 (2013).Article 

Google Scholar 
Coleman, M. L. & Chisholm, S. W. Code and context: Prochlorococcus as a model for cross-scale biology. Trends Microbiol. 15, 398–407 (2007).Article 
CAS 

Google Scholar 
Scanlan, D. J. et al. Ecological genomics of marine picocyanobacteria. Microbiol. Mol. Biol. Rev. 73, 249–299 (2009).Article 
CAS 

Google Scholar 
Moore, L. R., Rocap, G. & Chisholm, S. W. Physiology and molecular phylogeny of coexisting Prochlorococcus ecotypes. Nature 393, 464–467 (1998).Article 
ADS 
CAS 

Google Scholar 
Rusch, D. B., Martiny, A. C., Dupont, C. L., Halpern, A. L. & Venter, J. C. Characterization of Prochlorococcus clades from iron-depleted oceanic regions. Proc. Natl Acad. Sci. 107, 16184–16189 (2010).Article 
ADS 
CAS 

Google Scholar 
Larkin, A. A. et al. Persistent El Niño driven shifts in marine cyanobacteria populations. PloS ONE 15, e0238405 (2020).Article 
CAS 

Google Scholar 
Arandia‐Gorostidi, N. et al. Warming the phycosphere: differential effect of temperature on the use of diatom‐derived carbon by two copiotrophic bacterial taxa. Environ. Microbiol. 22, 1381–1396 (2020).Article 

Google Scholar 
Arandia‐Gorostidi, N., Huete‐Stauffer, T. M., Alonso‐Sáez L, G. & Morán, X. A. Testing the metabolic theory of ecology with marine bacteria: different temperature sensitivity of major phylogenetic groups during the spring phytoplankton bloom. Environ. Microbiol. 19, 4493–4505 (2017).Article 

Google Scholar 
Fagan, A. J., Moreno, A. R. & Martiny, A. C. Role of ENSO conditions on particulate organic matter concentrations and elemental ratios in the Southern California Bight. Front. Mar. Sci. 6, 386 (2019).Article 

Google Scholar 
Chang, C. W. et al. Reconstructing large interaction networks from empirical time series data. Ecol. Lett. 24, 2763–2774 (2021).Article 

Google Scholar 
Lie, A. A., Kim, D. Y., Schnetzer, A. & Caron, D. A. Small-scale temporal and spatial variations in protistan community composition at the San Pedro Ocean Time-series station off the coast of southern California. Aquat. Microb. Ecol. 70, 93–110 (2013).Article 

Google Scholar 
Yeh, Y.-C., Needham, D. M., Sieradzki, E. T. & Fuhrman, J. A. Taxon disappearance from microbiome analysis reinforces the value of mock communities as a standard in every sequencing run. MSystems 3, e00023–00018 (2018).Article 

Google Scholar 
Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010).Article 
CAS 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).Article 
CAS 

Google Scholar 
Guillou, L. et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D579–D604 (2013).
Google Scholar 
Amir, A. et al. Deblur rapidly resolves single-nucleotide community sequence patterns. MSystems 2, (2017).Decelle, J. et al. Phyto REF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy. Mol. Ecol. Resour. 15, 1435–1445 (2015).Article 
CAS 

Google Scholar 
Amin, S. et al. Interaction and signalling between a cosmopolitan phytoplankton and associated bacteria. Nature 522, 98–101 (2015).Article 
ADS 
CAS 

Google Scholar 
Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280 (2001).Article 
ADS 

Google Scholar 
Hill, M. O. & Gauch, H. G. J. Detrended correspondence analysis: an improved ordination technique. Vegetatio 42, 47–58 (1980).Ter Braak, C. J. Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis. Ecology 67, 1167–1179 (1986).Article 

Google Scholar 
Peres-Neto, P. R., Legendre, P., Dray, S. & Borcard, D. Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology 87, 2614–2625 (2006).Article 

Google Scholar  More

Rice, D. W. Marine mammals of the world: systematics and distribution. In The Society for Marine Mammalogy (ed. Rice, D. W.) 231 (Allen Press, 1998).
Google Scholar 
Best, P. B. External characters of southern minke whales and the existence of a diminutive form. Sci. Rep. Whales Res. Inst. 36, 1–33 (1985).
Google Scholar 
Acevedo, J. et al. Occurrence of dwarf minke whales (Balaenoptera acutorostrata subsp.) around the Antarctic Peninsula. Polar Biol. 34, 313–318 (2011).Article 

Google Scholar 
Risch, D., Norris, T., Curnock, M. & Friedlaender, A. Common and Antarctic minke whales: Conservation status and future research directions. Front. Mar. Sci. 6, 247. https://doi.org/10.3389/fmars.2019.00247 (2019).Article 

Google Scholar 
International Whaling Commission (IWC). Report of the scientific committee. J. Cetacean Res. Manag. 14, 102 (2013).
Google Scholar 
Matsuoka, K. et al. Overview of minke whale sightings surveys conducted on IWC/IDCR and SOWER Antarctic cruises from 1978/79 to 2000/01. J. Cetacean Res. Manag. 5, 173–201 (2003).
Google Scholar 
Glover, K. A. et al. Migration of Antarctic minke whales to the Arctic. PLoS One 5, e15197. https://doi.org/10.1371/journal.pone.0015197 (2010).Article 
ADS 
CAS 

Google Scholar 
Williams, R., Brierley, A., Friedlaender, A. & Scheidat, M. Densitiy of Antarctic minke whales in Weddell Sea from helicopter survey data. Ecology 63, IA14 (2011).
Google Scholar 
Williams, R. et al. Counting whales in a challenging, changing environment. Sci. Rep. 4, 4170. https://doi.org/10.1038/srep04170 (2014).Article 
CAS 

Google Scholar 
Shabangu, F. W., Findlay, K. & Stafford, K. M. Seasonal acoustic occurrence, diel vocalizing patterns and bioduck call-type composition of Antarctic minke whales off the west coast of South Africa and the Maud Rise Antarctica. Mar. Mamm. Sci. 36, 658–675 (2019).Article 

Google Scholar 
Kasamatsu, F., Nishiwaki, S. & Ishikawa, H. Breeding areas and southbound migrations of southern minke whales Balaenoptera acutorostrata. Mar. Ecol. Prog. Ser. 119, 1–10 (1995).Article 
ADS 

Google Scholar 
Tamura, T. & Konishi, K. Food habit and prey consumption of Antarctic minke whale Balaenoptera bonaerensis in the JARPA research area. J. Northwest Atl. Fish. Sci. 42, 13–25 (2009).Article 

Google Scholar 
Perrin, W. F., Mallette, S. D. & Brownell, R. L. Minke whales. In Encyclopedia of Marine Mammals (eds Perrin, W. F. et al.) 608–613 (Academic Press, 2018).Chapter 

Google Scholar 
Taylor, R. J. F. An unusual record of three species of whale being restricted to pools in Antarctic sea-ice. Proc. R. Soc. Lond. 129, 325–331 (1957).
Google Scholar 
Ensor, P. H. Minke whales in the pack ice zone, East Antarctica, during the period of maximum annual ice extent. Rep. Int. Whal. Commn 39, 219–225 (1989).
Google Scholar 
Scheidat, M. et al. Cetacean surveys in the Southern Ocean using icebreaker-supported helicopters. Polar Biol. 34, 1513–1522 (2011).Article 

Google Scholar 
Meirelles, A. C. O. & Furtado-Neto, M. A. A. Stranding of an Antarctic minke whale, Balaenoptera bonaerensis Burmeister, 1867, on the northern coast of South America. Lat. Am. J. Aquat. Mamm. 3, 81–82 (2004).Article 

Google Scholar 
Juri, E., Valdivia, M., Simoes-Lopes, P. C. & Le Bas, A. A note on minke whales (Cetacea: Balaenopteridae) in Uruguay: Strandings review. JCRM 21, 135–140 (2020).Article 

Google Scholar 
Williamson, G. R. Minke whales off Brazil. Sci. Rep. Whales Res. Inst. 27, 37–59 (1975).
Google Scholar 
Pastene, L. A. & Goto, M. Genetic characterization and population genetic structure of the Antarctic minke whale Balaenoptera bonaerensis in the Indo-Pacific region of the Southern Ocean. Fish Sci. 82, 873–886 (2016).Article 
CAS 

Google Scholar 
Balbuena, J. A., Aznar, F. J., Fernández, M. & Raga, J. A. Parasites as indicators of social structure and stock identity of marine mammals. Dev. Mar. Biol. 4, 133–139 (1995).
Google Scholar 
Kuramochi, T., Araki, J., Uchida, Moriyama, N., Takeda, Y., Hayashi, N., Wakao, H., Machida, M. & Nagasawa, K. Summary of parasite and epizoit investigations during JARPN surveys 1994–1999, with reference to stock structure analysis for the western North Pacific minke whales. IWC Scientific Committee Workshop to Review the Japanese Whaling Programme under Special Permit for North Pacific Minke Whales (JARPN) SC/F2K/J19 (2000).Kaliszewska, Z. A. et al. Population histories of right whales (Cetacea: Eubalaena) inferred from mitochondrial sequence diversities and divergences of their whale lice (Amphipoda: Cyamus). Mol. Ecol. 14, 3439–3456 (2005).Article 
CAS 

Google Scholar 
Ólafsdóttir, D. & Shinn, A. P. Epibiotic macrofauna on common minke whales, Balaenoptera acutorostrata Lacépède, 1804 Icelandic waters. Parasit. Vectors 6, 1–10 (2013).Article 

Google Scholar 
Matthews, C. J., Ghazal, M., Lefort, K. J. & Inuarak, E. Epizoic barnacles on Arctic killer whales indicate residency in warm waters. Mar. Mamm. Sci. 36, 1010–1014 (2020).Article 

Google Scholar 
Flach, L., Van Bressem, M. F., Pitombo, F. & Aznar, F. J. Emergence of the epibiotic barnacle Xenobalanus globicipitis in Guiana dolphins after a morbillivirus outbreak in Sepetiba Bay Brazil. Estuar. Coast. Shelf Sci. 263, 107632. https://doi.org/10.1016/j.ecss.2021.107632 (2021).Article 

Google Scholar 
Ten, S., Raga, J. A. & Aznar, F. J. Epibiotic fauna on cetaceans worldwide: A systematic review of records and indicator potential. Front. Mar. Sci. 9, 846558. https://doi.org/10.3389/fmars.2022.846558 (2022).Article 

Google Scholar 
Liouville, J. Cétacés de l’Antarctique. Paris: Deuxième Expédition Antarctique Française (1908–1910) (1913).Ohsumi, S., Masaki, Y. & Kawamura, A. Stock of the Antarctic minke whale. Sci. Rep. Whales Res. Inst. 22, 75–125 (1970).
Google Scholar 
Ohsumi, S. Find of marlin spear from the Antarctic minke whales. Sci. Rep. Whales Res. Inst. 25, 237–239 (1973).
Google Scholar 
Ivashin, M. V. External Parasites on Lesser Rorquals in the Antarctic 125–127 (Naukova Dumka, 1975).
Google Scholar 
Berzin, A. A. & Vlasova, L. P. Fauna of the Cetacea Cyamidae (Amphipoda) of the world ocean. Investig. Cet. 13, 149–164 (1982).
Google Scholar 
Best, P. B. Seasonal abundance, feeding, reproduction, age and growth in minke whales off Durban (with incidental observations from the Antarctic). Rep. Int. Whal. Commn 32, 759–786 (1982).
Google Scholar 
Avdeev, V. V. Parasitic amphipods of the family Cyamidae and the problem of Cetacea origin. Biol. Morja 4, 27–33 (1989).
Google Scholar 
Bushuev, S. G. A study of the population structure of the southern minke whale (Balaenoptera acutorostrata Lacepede) based on morphological and ecological variability. Rep. Int. Whal. Commn 40, 317–324 (1990).
Google Scholar 
Sedlak-Weinstein, E. Preliminary report of parasitic infestation of the minke whale Balaenoptera acutorostrata taken during the 1988/89 Antarctic expedition. Unpublished paper (1990).Dailey, M. D. & Vogelbein, W. Parasite fauna of 3 species of Antarctic whales with reference to their use as potential stock indicators. Fish. Bull. 89, 355–365 (1991).
Google Scholar 
Nemoto, T., Best, P. B., Ishimaru, K. & Takano, H. Diatom films on whales in South African waters. Sci. Rep. Whales Res. Inst. 32, 97–103 (1980).
Google Scholar 
Donovan, G. A review of IWC stock boundaries. Rep. Int. Whal. Commn 13, 39–68 (1991).
Google Scholar 
Lester, R. J. G. & MacKenzie, K. The use and abuse of parasites as stock markers for fish. Fish. Res. 97, 1–2 (2009).Article 

Google Scholar 
Ten, S. et al. Epibiotic barnacles of sea turtles as indicators of habitat use and fishery interactions: an analysis of juvenile loggerhead sea turtles, Caretta caretta, in the western Mediterranean. Ecol. Indic. 107, 105672. https://doi.org/10.1016/j.ecolind.2019.105672 (2019).Article 

Google Scholar 
Calman, W. T. A whale-barnacle of the genus Xenobalanus from Antarctic Seas. Ann. Mag. Nat. Hist. 6, 165–166 (1920).Article 

Google Scholar 
Kato, H., Hiroyama, H., Fujise, Y. & Ono, K. Preliminary report of the 1987/88 Japanese feasibility study of the special permit proposal for Southern Hemisphere Minke Whales. Rep. int. Whal. Commn 39, 235–248 (1989).
Google Scholar 
International Whaling Commission (IWC). Report of the Intersessional Workshop to review data and results from special permit research on minke whales in the Antarctic, Tokyo, 7–8 December 2006. J. Cetacean Res. Manag. 10, 411–445 (2008).
Google Scholar 
Bush, A. O., Lafferty, K. D., Lotz, J. M. & Shostak, A. W. Parasitology meets ecology on its own terms: Margolis et al. revisited. J. Parasitol. 83, 575–583 (1997).Article 
CAS 

Google Scholar 
Kim, H., Chan, B., Kang, C., Kim, H. & Kim, W. How do whale barnacles live on their hosts? Functional morphology and mating-group sizes of Coronula diadema (Linnaeus, 1767) and Conchoderma auritum (Linnaeus, 1767) (Cirripedia: Thoracicalcarea). J. Crustac. Biol. 40, 808–824 (2020).Article 

Google Scholar 
Reiczigel, J. Confidence intervals for the binomial parameter: Some new considerations. Stat. Med. 22, 611–621 (2003).Article 

Google Scholar 
Kato, H. Migration strategy of southern minke whales to maintain high reproductive rate. Dev. Mar. Biol. 4, 465–480 (1995).
Google Scholar 
Zuur, A. F., Ieno, E. N., Walker, N. J., Saveliev, A. A. & Smith, G. M. Mixed effects models and extensions in ecology with R. In Statistics for Biology and Health (ed. Gail, M.) (Springer, 2009).MATH 

Google Scholar 
Fransen, C. H. J. M. & Smeenk, C. Whale-lice (Amphipoda: Cyamidae) recorded from The Netherlands. Zool. Meded. 65, 393–405 (1991).
Google Scholar 
Barton, N. A., Farewell, T. S. & Hallett, S. H. Using generalized additive models to investigate the environmental effects on pipe failure in clean water networks. NPJ Clean Water 3, 31. https://doi.org/10.1038/s41545-020-0077-3 (2020).Article 

Google Scholar 
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682. https://doi.org/10.1038/nmeth.2019 (2012).Article 
CAS 

Google Scholar 
Kane, E. A., Olson, P. A., Gerrodette, T. & Fiedler, P. Prevalence of the commensal barnacle Xenobalanus globicipitis on cetacean species in the eastern tropical Pacific Ocean, and a review of global occurrence. Fish. Bull. 106, 395–404 (2008).
Google Scholar 
Aznar, F. J., Balbuena, J. A. & Raga, J. A. Are epizoites biological indicators of a western Mediterranean striped dolphin die-off?. Dis. Aquat. Organ. 18, 159–163 (1994).Article 

Google Scholar 
Carrillo, J. M., Overstreet, R. M., Raga, J. A. & Aznar, F. J. Living on the edge: Settlement patterns by the symbiotic barnacle Xenobalanus globicipitis on small cetaceans. PLoS One 10, e0127367. https://doi.org/10.1371/journal.pone.0127367 (2015).Article 
CAS 

Google Scholar 
Moreno-Colom, P., Ten, S., Raga, J. A. & Aznar, F. J. Spatial distribution and aggregation of Xenobalanus globicipitis on the flukes of striped dolphins, Stenella coeruleoalba: An indicator of host hydrodynamics?. Mar. Mamm. Sci. 36, 897–914 (2020).Article 

Google Scholar 
Aznar, F. J. et al. Changes in epizoic crustacean infestations during cetacean die-offs: The mass mortality of Mediterranean striped dolphins Stenella coeruleoalba revisited. Dis. Aquat. Org. 67, 239–247 (2005).Article 
CAS 

Google Scholar 
Wood, S. N. & Augustin, N. H. GAMs with integrated model selection using penalized regression splines and applications to environmental modelling. Ecol. Modell. 157, 157–177 (2002).Article 

Google Scholar 
Wood, S. N. Generalized Additive Models: An Introduction with R (Chapman and Hall/CRC, 2017).Book 
MATH 

Google Scholar 
Bloch, D. et al. Short-term movements of long-finned pilot whales Globicephala melas around the Faroe Islands. Wildl. Biol. 9, 47–58 (2003).Article 

Google Scholar 
Beasley, I. et al. Stomach contents of long-finned pilot whales, Globicephala melas mass-stranded in Tasmania. PLoS One 14, e0206747. https://doi.org/10.1371/journal.pone.0206747 (2019).Article 
CAS 

Google Scholar 
Ohno, M. & Fujino, K. Biological investigation on the whales caught by the Japanese Antarctic whaling fleets, season 1950/51. Sci. Rep. Whales Res. Inst. 7, 125–188 (1952).
Google Scholar 
Clarke, R. The stalked barnacle Conchoderma, ectoparasitic on whales. Norsk Hvalfangst-Tidende 55, 153–168 (1966).
Google Scholar 
Christensen, I. First record of gooseneck barnacles (Conchoderma auritum) on a minke whale (Balaenoptera acutorostrata). ICES C. M. 1985/N:9 (1985).Bertulli, C. G., Cecchetti, A., Van Bressem, M. F. & Van Waerebeek, K. Skin disorders in common minke whales and white-beaked dolphins off Iceland, a photographic assessment. J. Mar. Anim. Ecol. 5, 29–40 (2012).
Google Scholar 
Knowlton, N. Sibling species in the sea. Annu. Rev. Ecol. Evol. Syst. 24, 189–216 (1993).Article 

Google Scholar 
Trontelj, P. & Fišer, C. Perspectives: Cryptic species diversity should not be trivialised. Syst. Biodivers. 7, 1–3 (2009).Article 

Google Scholar 
Norris, R. & Hull, P. The temporal dimension of marine speciation. Evol. Ecol. 26, 393–415 (2011).Article 

Google Scholar 
Rawson, P., Macnamee, R., Frick, M. & Williams, K. Phylogeography of the coronulid barnacle, Chelonibia testudinaria, from loggerhead sea turtles Caretta caretta. Mol. Ecol. 12, 2697–2706 (2003).Article 
CAS 

Google Scholar 
Cabezas, M. P., Cabezas, P., Machordom, A. & Guerra-García, J. M. Hidden diversity and cryptic speciation refute cosmopolitan distribution in Caprella penantis (Crustacea: Amphipoda: Caprellidae). J. Zool. Syst. Evol. 51, 85–99 (2013).Article 

Google Scholar 
Boyd, L. L., Zardus, J. D., Knauer, C. M. & Wood, L. D. Evidence for host selectivity and specialization by epizoic Chelonibia barnacles between hawksbill and green sea turtles. Front. Ecol. Evol. 9, 807237. https://doi.org/10.3389/fevo.2021.807237 (2021).Article 

Google Scholar 
Schell, D., Rowntree, V. & Pfeiffer, C. Stable-isotope and electron-microscopic evidence that cyamids (Crustacea: Amphipoda) feed on whale skin. Can. J. Zool. 78, 721–727 (2000).Article 

Google Scholar 
Iwasa-Arai, T. & Serejo, C. S. Phylogenetic analysis of the family Cyamidae (Crustacea: Amphipoda): A review based on morphological characters. Zool. J. Linn. Soc. 184, 66–94 (2018).Article 

Google Scholar 
Fraija-Fernández, N. et al. Living in a harsh habitat: Epidemiology of the whale louse, Syncyamus aequus (Cyamidae), infecting striped dolphins in the Western Mediterranean. J. Zool. 303, 199–206 (2017).Article 

Google Scholar 
Angot, M. Rapport scientifique sur les expeditions baleinieres autour de Madagascar (saisons 1949 et 1950). Mem. Inst. Sci. Madag. Ser. A 6, 439–486 (1951).
Google Scholar 
Newman, W. A. & Abbott, D. P. Cirripedia: The barnacles. In Intertidal Invertebrates of California (eds Morris, R. H. et al.) 504–535 (Stanford University Press, 1980).
Google Scholar 
Nogata, Y. & Matsumura, K. Larval development and settlement of a whale barnacle. Biol Lett. 2, 92–93 (2006).Article 

Google Scholar 
Hiro, F. The fauna of Akkeshi Bay. II. Cirripedia. J. Fac. Sci. Hokkaido Univ. 4, 213–229 (1935).
Google Scholar 
Rice, D. W. Progress report on biological studies of the larger Cetacea in the waters off California. Norsk Hvalfangst-Tid 52, 181–187 (1963).
Google Scholar 
Klinkhart, E. G. The beluga whale in Alaska. State Alsk. Dep. Fish 7, 11 (1966).
Google Scholar 
Nilsson-Cantell, C. A. Cirripedia Thoracica and Acrothoracica. MIOS 5, 1–133 (1978).
Google Scholar 
Scarff, J. E. Occurrence of the barnacles Coronula diadema, C. reginae and Cetopirus complanatus (Cirripedia) on right whales. Sci. Rep. Whales Res. Inst. 37, 129–153 (1986).
Google Scholar 
Kakuwa, Z., Kawakami, T. & Iguchi, K. Biological investigation on the whales caught by the Japanese Antarctic whaling fleets in the 1951–52 season. Sci. Rep. Whales Res. Inst. 8, 147–213 (1953).
Google Scholar 
Nishiwaki, M. Humpback whales in Ryukyuan waters. Sci. Rep. Whales Res. Inst. 14, 49–87 (1959).
Google Scholar 
Best, P. B. The presence of coronuline barnacles on a southern right whale Eubalaena australis. S. Afr. J. Mar. Sci. 11, 585–587 (1991).Article 

Google Scholar 
Mackintosh, N. A. & Wheeler, J. F. G. Southern blue and fin whales. Disc. Rep. 1, 257–540 (1929).
Google Scholar 
Nilsson-Cantell, C. A. Thoracic cirripedes collected in 1925–1927. Disc. Rep. 2, 223–260 (1930).
Google Scholar 
Nishiwaki, M. & Hayashi, K. Biological survey of fin and blue whales taken in the Antarctic season 1947–48 by the Japanese fleet. Sci. Rep. Whales Res. Inst. 3, 132–190 (1950).
Google Scholar 
Mizue, K. & Murata, T. Biological investigation on the whales caught by the Japanese Antarctic whaling fleets season 1949–50. Sci. Rep. Whales Res. Inst. 6, 73–131 (1951).
Google Scholar 
Nishiwaki, M. & Oye, T. Biological investigation on blue whales (Balaenoptera musculus) and Fin Whales (Balaenoptera physalus) caught by the Japanese Antarctic Whaling Fleets. Sci. Rep. Whales Res. Inst. 5, 91–167 (1951).
Google Scholar 
Tomilin, A. G. Cetacea. In Mammals of the U.S.S.R. and Adjacent Countries Vol. 9 (ed. Tomilin, A. G.) 717 (Akademii Nauk SSSR, 1957).
Google Scholar 
Cockrill, W. R. Pathology of the cetacea. A veterinary study on whales. Br. Vet. J. 116, 1–28 (1960).
Google Scholar 
Kawamura, A. Some consideration on the stock unit of sei whales by the aspect of ectoparasitic organisms on the body. Bull. Jpn. Soc. Fish. Oceanogr. 14, 38–43 (1969).
Google Scholar 
Fraija-Fernández, N., Hernández-Hortelano, A., Ahuir-Baraja, A. E., Raga, J. A. & Aznar, F. J. Taxonomic status and epidemiology of the mesoparasitic copepod Pennella balaenoptera in cetaceans from the western Mediterranean. Dis. Aquat. Org. 128, 249–258 (2018).Article 

Google Scholar 
Foster, B. A. & Willan, R. C. Foreign barnacles transported to New Zealand on an oil platform. N. Z. J. Mar. Freshw. Res. 13, 143–149 (1979).Article 

Google Scholar 
González, J. et al. Cirripedia of the Canary islands: Distribution and ecological notes. J. Mar. Biol. Assoc. U.K. 92, 129–141 (2012).Article 

Google Scholar 
Zettler, M. L. An example for transatlantic hitchhiking by macrozoobenthic organisms with a research vessel. Helgol. Mar. Res. 75, 4. https://doi.org/10.1186/s10152-021-00549-w (2021).Article 

Google Scholar 
Matthews, L. H. The humpback whale Megaptera novaeangliae. Disc. Rep. 17, 7–92 (1937).
Google Scholar 
Scheffer, V. B. Organisms collected from whales in the Aleutian Islands. Murrelet 20, 67–69 (1939).Article 

Google Scholar 
Symons, H. W. & Weston, R. D. Studies on the humpback whale (Megaptera nodosa) in the Bellinghausen Sea. Norsk Hvalfangsttid 47, 53–81 (1958).
Google Scholar 
Van Waerebeek, K., Reyes, J. C. & Alfaro, J. Helminth parasites and phoronts of dusky dolphins Lagenorhynchus obscurus (Gray, 1828) from Peru. Aquat. Mamm. 19, 159–169 (1993).
Google Scholar 
Fertl, D. Barnacles. In Encyclopedia of Marine Mammals (eds Perrin, W. F. et al.) 75–78 (Academic Press, 2002).
Google Scholar 
Cornwall, I. E. The barnacles of british Columbia. Br. Col. Prov. Mus. Dept. 7, 5–69 (1955).
Google Scholar 
Abaunza, P., Arroyo, N. L. & Preciado, I. A contribution to the knowledge on the morphometry and the anatomical characters of Pennella balaenopterae (Copepoda, Ciphonostomatoida, Pennellidae), with special reference to the buccal complex. Crustaceana 74, 193–210 (2001).Article 

Google Scholar 
Marcer, F. et al. Parasitological and pathological findings in fin whales Balaenoptera physalus stranded along Italian coastlines. Dis. Aquat. Org. 133, 25–37 (2019).Article 
CAS 

Google Scholar 
Turner, W. On Pennella balænopteræ: A crustacean, parasitic on a finner whale, Balaenoptera musculus. Earth. Environ. Sci. Trans. R. Soc. Edinb. 41, 409–434 (1905).Article 

Google Scholar 
Walker, W. A. & Hanson, M. B. Biological observations on Stejneger’s beaked whale, Mesoplodon stejnegeri, from strandings on Adak Alaska. Mar. Mamm. Sci. 15, 1314–1329 (1999).Article 

Google Scholar 
Delaney, M. A., Ford, J. K. B., Tang, K. & Gaydos, J. K. Mesoparasitic copepod (Pennella balaenopterae) infestation of a stranded offshore orca (Orcinus orca) in Southeast Alaska: Review of significance as a health indicator in cetaceans. In IAAAM 21–26 (2016).Suyama, S., Kakehi, S., Yanagimoto, T. & Chow, S. Infection of the pacific saury Cololabis saira (Brevoort, 1856) (Teleostei: Beloniformes: Scomberesocidae) by Pennella sp. (Copepoda: Siphonostomatoida: Pennellidae) south of the Subarctic Front. J. Crust. Biol. 40, 384–389 (2020).Article 

Google Scholar 
Rowntree, V. J. Feeding, distribution and reproductive behavior of cyamids (Crustacea: Amphipoda) living on humpback and right whales. Can. J. Zool. 74, 103–109 (1996).Article 

Google Scholar 
Leung, Y. M. Life cycle of Cyamus scammoni (Amphipoda: Cyamidae), ectoparasite of gray whale, with a remark on the associated species. Sci. Rep. Whales Res. Inst. 28, 153–160 (1976).
Google Scholar 
MacIntyre, R. J. Rapid growth in stalked barnacles. Nature 212, 637–638 (1966).Article 
ADS 

Google Scholar 
Rasmussen, T. Notes on the biology of the shipfouling gooseneck barnacle Conchoderma auritum Linnaeus, 1776 (Cirripedia; Lepadomorpha). Biol. Mar. 2, 37–44 (1980).
Google Scholar 
Dalley, R. & Crisp, D. J. Conchoderma: A fouling hazard to ships underway. Mar. Biol. Lett. 2, 141–152 (1981).
Google Scholar 
Dalley, R. The larval stages of the oceanic, pedunculate barnacle Conchoderma auritum (L) (Cirripedia, Thoracica). Crustaceana 46, 39–54 (1984).Article 

Google Scholar 
Foskolos, I., Provata, M. T. & Frantzis, A. First record of Conchoderma auritum (Cirripedia: Lepadidae) on Ziphius cavirostris (Cetacea: Ziphiidae) in Greece. Ann. Ser. Hist. 27, 29–34 (2017).
Google Scholar 
Lee, J. F., Friedlaender, A. S., Oliver, M. J. & DeLiberty, T. L. Behavior of satellite-tracked Antarctic minke whales (Balaenoptera bonaerensis) in relation to environmental factors around the western Antarctic Peninsula. Anim. Biotelem. 5, 23. https://doi.org/10.1186/s40317-017-0138-7 (2017).Article 

Google Scholar 
Darwin, C. A Monograph on the Subclass Cirripedia Vol. 1 (The Ray Society, 1851).
Google Scholar 
Tsikhon-Lukanina, V. A., Soldatova, I. N., Kuznetsova, I. A. & Il’in, I. I. Macrofouling community in the Strait of Tunisia (Sicily). Oceanology 16, 519–522 (1977).
Google Scholar 
Nilsson-Cantell, C. A. Cirripedien von der Stewart Insel und von Südgeorgien. Senckenbergiana 12, 210–213 (1930).
Google Scholar 
Slijper, E. J. Whales (Hutchinson, 1962).
Google Scholar 
Kaufman, G. D. & Forestell, P. H. Hawaii’s humpback whales, a complete whalewatching guide (Pacific Whale Foundation Press, 1986).
Google Scholar 
Dawbin, W. H. Baleen whales. In Whales, Dolphins and Porpoises (eds Harrison, R. & Bryden, M.) 44–65 (Facts on File, 1988).
Google Scholar 
Félix, F., Bearson, B. & Falconí, J. Epizoic barnacles removed from the skin of a humpback whale after a period of intense surface activity. Mar. Mamm. Sci. 22, 979–984 (2006).Article 

Google Scholar 
Towers, J. R. et al. Seasonal movements and ecological markers as evidence for migration of common minke whales photo-identified in the eastern North Pacific. J. Cetacean Res. Manag. 13, 221–229 (2013).
Google Scholar 
Iwasa-Arai, T. et al. The host-specific whale louse (Cyamus boopis) as a potential tool for interpreting humpback whale (Megaptera novaeangliae) migratory routes. J. Exp. Mar. Biol. Ecol. 505, 45–51 (2018).Article 

Google Scholar 
Lehnert, K. et al. Whale lice (Isocyamus deltobranchium & Isocyamus delphinii; Cyamidae) prevalence in odontocetes off the German and Dutch coasts – Morphological and molecular characterization and health implications. Int. J. Parasitol. 15, 22–30 (2021).
Google Scholar 
Dreyer, N. et al. How whale and dolphin barnacles attach to their hosts and the paradox of remarkably versatile attachment structures in cypris larvae. Org. Divers. Evol. 20, 233–249 (2020).Article 

Google Scholar 
Visser, I. N., Cooper, T. E. & Grimm, H. Duration of pseudo-stalked barnacles (Xenobalanus globicipitis) on a New Zealand Pelagic ecotype orca (Orcinus orca), with comments on cookie cutter shark bite marks (Isistius sp.); can they be used as biological tags?. Biol. Divers. 11, 1067–1086 (2020).
Google Scholar 
Van Waerebeek, K. & Reyes, J. C. A note on incidental fishery mortality of southern minke whales off western South America. Rep. Int. Whal. Commn 15, 521–523 (1994).
Google Scholar 
Félix, F. & Haase, B. A note on the northernmost record of the Antarctic minke whale (Balaenoptera bonaerensis) in the Eastern Pacific. J. Cetacean Res. Manag. 13, 191–194 (2013).
Google Scholar 
Esposito, C., Bichet, O. & Petit, M. First sightings of Antarctic minke whale (Balaenoptera bonaerensis) mother–calf pairs in French Polynesia. Aquat. Mamm. 47, 175–180 (2021).Article 

Google Scholar 
Karaa, S., Insacco, G., Bradai, M. N. & Scaravelli, D. Records of Xenobalanus globicipitis on Balaenoptera physalus and Stenella coeruleoalba in Tunisian and Sicilian waters. Nat. Rerum 1, 55–59 (2011).
Google Scholar 
Oliveira, J. B., Morales, J. A., González-Barrientos, R. C., Hernández-Gamboa, J. & Hernández-Mora, G. Parasites of cetaceans stranded on the Pacific Coast of Costa Rica. Vet. Parasitol. 182, 319–328. https://doi.org/10.1016/j.vetpar.2011.05.014 (2011).Article 
CAS 

Google Scholar 
Dı́az-Gamboa, R. E. Varamiento de orcas pigmeas (Feresa attenuata Gray 1874) en Yucatán: Reporte de caso. Bioagrociencias 8, 36–43 (2015).
Google Scholar 
IJsseldijk, L. L. et al. Beached bachelors: An extensive study on the largest recorded sperm whale Physeter macrocephalus mortality event in the north sea. PloS One 13, e0201221. https://doi.org/10.1371/journal.pone.0201221 (2018).Article 
CAS 

Google Scholar 
Guerrero-Ruiz, M. & Urbán, J. R. First report of remoras on two killer whales (Orcinus orca) in the Gulf of California Mexico. Aquat. Mamm. 26, 148–150 (2000).
Google Scholar 
Kautek, G., Van Bressem, M. F. & Ritter, F. External body conditions in cetaceans from La Gomera, Canary Islands Spain. J. Marine Anim. Ecol. 11, 4–17 (2008).
Google Scholar 
Bearzi, M. & Patonai, K. Occurrence of the barnacle (Xenobalanus globicipitis) on coastal and offshore common bottlenose dolphins (Tursiops truncatus) in Santa Monica Bay and adjacent areas California. Bull. S. Calif. Acad. Sci. 109, 37–44. https://doi.org/10.3160/0038-3872-109.2.37 (2010).Article 

Google Scholar 
Foote, A. D. et al. Genetic differentiation among North Atlantic killer whale populations. Mol. Ecol. 20, 629–641. https://doi.org/10.1111/j.1365-294X.2010.04957.x (2011).Article 

Google Scholar 
Toth, J. L., Hohn, A. A., Able, K. W. & Gorgone, A. M. Defining bottlenose dolphin (Tursiops truncatus) stocks based on environmental, physical and behavioral characteristics. Mar. Mamm. Sci. 28, 461–478. https://doi.org/10.1111/j.1748-7692.2011.00497.x (2012).Article 

Google Scholar 
Urian, K. W., Kaufmann, R., Waples, D. M. & Read, A. J. The prevalence of ectoparasitic barnacles discriminates stocks of Atlantic common bottlenose dolphins (Tursiops truncatus) at risk of entanglement in coastal gill net fisheries. Mar. Mamm. Sci. 35, 290–299. https://doi.org/10.1111/mms.12522 (2019).Article 

Google Scholar 
Siciliano, S. et al. Epizoic barnacle (Xenobalanus globicipitis) infestations in several cetacean species in South-Eastern Brazil. Mar. Biol. Res. 16, 1–13. https://doi.org/10.1080/17451000.2020.1783450 (2020).Article 

Google Scholar 
Whitehead, T. O., Rollinson, D. P. & Reisinger, R. R. Pseudostalked barnacles Xenobalanus globicipitis attached to killer whales Orcinus orca in South African waters. Mar. Biodivers. Rec. 45, 873–876. https://doi.org/10.1007/s12526-014-0296-2 (2014).Article 

Google Scholar 
Methion, S. & Dı́az López, B. First record of atypical pigmentation pattern in fin whale Balaenoptera physalus in the Atlantic ocean. Dis. Aquat. Org. 135, 121–125. https://doi.org/10.3354/dao03385 (2019).Article 

Google Scholar 
Herr, H., Burkhardt-Holm, P., Heyer, K., Siebert, U. & Selling, J. Injuries, malformations and epidermal conditions in cetaceans of the strait of Gibraltar. Aquat. Mamm. 46, 215–235. https://doi.org/10.1578/AM.46.2.2020.215 (2020).Article 

Google Scholar 
Herr, H. et al. Return of large fin whale feeding aggregations to historical whaling grounds in the southern ocean. Sci. Rep. 12, 9458. https://doi.org/10.1038/s41598-022-13798-7 (2022).Article 
ADS 
CAS 

Google Scholar 
Gruvel, J. A. Cirrhipèdes Provenant Des Campagnes Scientifiques De S.A.S. Le Prince De Monaco, (1885– 1913). In Résultas Des Campagnes Scientifiques Accomplies Sur Son Yacht Par Albert Ler (Monaco: Prince Souverain de Monaco) 1-88 (1920).Annandale, N. The rate of growth in Conchoderma and Lepas. Rec. Indian Mus. 3, 295 (1909).
Google Scholar 
Il’in, I. I., Kuznetsova, L. A. & Starostin, I. V. Oceanic fouling in the equatorial Atlantic. Oceanology 18, 597–599 (1978).
Google Scholar 
Eckert, K. L. & Eckert, S. A. Growth rate and reproductive condition of the barnacle Conchoderma virgatum on gravid leatherback sea turtles in Caribbean waters. J. Crust. Biol. 7, 682–690. https://doi.org/10.2307/1548651 (1987).Article 

Google Scholar 
Arroyo, N. L., Abaunza, P. & Preciado, I. The first naupliar stage of Pennella balaenopterae Koren and Danielssen 1877 (Copepoda: Siphonostomatoida, Pennellidae). Sarsia 87, 333–337. https://doi.org/10.1080/0036482021000155785 (2002).Article 

Google Scholar  More

Evans, N. P. et al. Quantification of drought during the collapse of the classic Maya civilization. Science 361, 498–501 (2018).Article 
ADS 
CAS 

Google Scholar 
Gill, R. B. The Great Maya Droughts: Water, Life, and Death (University of New Mexico Press, 2001).
Google Scholar 
Coe, M. D. The Maya (Thames and Hudson, 1993).
Google Scholar 
Douglas, P. M. J. et al. Drought, agricultural adaptation, and sociopolitical collapse in the Maya Lowlands. Proc. Natl. Acad. Sci. USA 112, 5607–5612 (2015).Article 
ADS 
CAS 

Google Scholar 
Haug, G. H. et al. Climate and the collapse of Maya civilization. Science 299, 1731–1735 (2003).Article 
ADS 
CAS 

Google Scholar 
Ford, A. & Nigh, R. Origins of the Maya forest garden: Maya resource management. J. Ethnobiol. 29, 213–236 (2009).Article 

Google Scholar 
Anderson, E. N. et al. Las Plantas de los Mayas: Etnobotánica en Quintana Roo, México (CONABIO-ECOSUR, 2005).
Google Scholar 
Fedick, S. L. Maya cornucopia: Indigenous food plants of the Maya lowlands. in The Real Business of Ancient Maya Economies (eds. Masson, M. A., Freidel, D. A. & Demarest, A. A.). 224–237 (University Press Florida, 2020).Ford, A. & Clarke, K. C. Linking the past and present of the ancient Maya: Lowland land use, population distribution, and density in the Late Classic period. in The Oxford Handbook of Historical Ecology and Applied Archaeology (eds. Isendahl, C. & Stump, D.) (Oxford Handbook of Historical Ecology and Applied Archaeology, 2015).Ford, A. & Nigh, R. The Maya Forest Garden: Eight Millennia of Sustainable Cultivation of the Tropical Woodlands (Routledge, 2016).Gómez-Pompa, A. On maya silviculture. Mexican Stud. (Estudios Mexicanos) 3, 1–17 (1987).Article 

Google Scholar 
Beach, T., Luzzadder-Beach, S., Krause, S. & Walling, S. ‘Mayacene’ floodplain and wetland formation in the Rio Bravo watershed of northwestern Belize. Holocene 25(10), 1612–1622 (2015).Pohl, M. D. et al. Early agriculture in the Maya lowlands. Lat. Am. Antiq. 7, 355–372 (1996).Article 

Google Scholar 
Fedick, S. L. The Managed Mosaic: Ancient Maya Agriculture and Resource Use (University of Utah Press, 1996).
Google Scholar 
Mueller, A. D. et al. Recovery of the forest ecosystem in the tropical lowlands of northern Guatemala after disintegration of Classic Maya polities. Geology 38, 523–526 (2010).Article 
ADS 

Google Scholar 
Hodell, D. A., Curtis, J. H. & Brenner, M. Possible role of climate in the collapse of Classic Maya civilization. Nature 375, 391–394 (1995).Article 
ADS 
CAS 

Google Scholar 
Islebe, G. A., Hooghiemstra, H., Brenner, M., Curtis, J. H. & Hodell, D. A. A Holocene vegetation history from lowland Guatemala. Holocene 6, 265–271 (1996).Article 
ADS 

Google Scholar 
Medina-Elizalde, M., Polanco-Martínez, J. M., Lases-Hernández, F., Bradley, R. & Burns, S. Testing the ‘tropical storm’ hypothesis of Yucatan Peninsula climate variability during the Maya Terminal Classic Period. Quat. Res. 86, 111–119 (2016).Aragón-Moreno, A. A., Islebe, G. A., Torrescano-Valle, N. & Arellano-Verdejo, J. Middle and late Holocene mangrove dynamics of the Yucatan Peninsula, Mexico. J. South Am. Earth Sci. 85, 307–311 (2018).Article 
ADS 

Google Scholar 
Aragón-Moreno, A. A., Islebe, G. A., Roy, P. D., Torrescano-Valle, N. & Mueller, A. D. Climate forcings on vegetation of the southeastern Yucatán Peninsula (Mexico) during the middle to late Holocene. Palaeogeogr. Palaeoclimatol. Palaeoecol. 495, 214–226 (2018).Article 

Google Scholar 
Kennett, D. J. et al. Development and disintegration of Maya political systems in response to climate change. Science 338, 788–791 (2012).Article 
ADS 
CAS 

Google Scholar 
Conde, C. et al. El Niño y la agricultura. in Los impactos de El Niño en México (ed. Magaña, V.). 103–135 (Dirección General de Protección Civil, Secretaría de Gobernación, México, 1999).Magaña, V. O., Vázquez, J. L., Pérez, J. L. & Pérez, J. B. Impact of El Niño on precipitation in Mexico. Geofísica Int. 42, 313–330 (2003).
Google Scholar 
Wahl, D., Byrne, R. & Anderson, L. An 8700 year paleoclimate reconstruction from the southern Maya lowlands. Quat. Sci. Rev. 103, 19–25 (2014).Article 
ADS 

Google Scholar 
Nooren, K. et al. Climate impact on the development of Pre-Classic Maya civilisation. Clim. Past 14, 1253–1273 (2018).Article 

Google Scholar 
Palomo-Kumul, J., Valdez-Hernández, M., Islebe, G. A., Cach-Pérez, M. J. & El Andrade, J. L. Niño-Southern oscillation affects the water relations of tree species in the Yucatan Peninsula. Mexico. Sci. Rep. 11, 10451 (2021).Article 
ADS 
CAS 

Google Scholar 
Rosenswig, R. M., VanDerwarker, A. M., Culleton, B. J. & Kennett, D. J. Is it agriculture yet? Intensified maize-use at 1000 cal BC in the Soconusco and Mesoamerica. J. Anthropol. Archaeol. 40, 89–108 (2015).Article 

Google Scholar 
Mueller, A. D. et al. Climate drying and associated forest decline in the lowlands of northern Guatemala during the late Holocene. Quat. Res. 71, 133–141 (2009).Article 

Google Scholar 
Aragón-Moreno, A. A., Islebe, G. A. & Torrescano-Valle, N. A ~3800-yr, high-resolution record of vegetation and climate change on the north coast of the Yucatan Peninsula. Rev. Palaeobot. Palynol. 178, 35–42 (2012).Article 

Google Scholar 
Carrillo-Bastos, A., Islebe, G. A. & Torrescano-Valle, N. 3800 Years of quantitative precipitation reconstruction from the Northwest Yucatan Peninsula. PLoS ONE 8, e84333 (2013).Article 
ADS 

Google Scholar 
Berglund, B. E. Human impact and climate changes—Synchronous events and a causal link?. Quat. Int. 105, 7–12 (2003).Article 

Google Scholar 
Vela-Peláez, A. A., Torrescano-Valle, N., Islebe, G. A., Mas, J. F. & Weissenberger, H. Holocene precipitation changes in the Maya forest, Yucatán peninsula. Mexico. Palaeogeogr. Palaeoclimatol. Palaeoecol. 505, 42–52 (2018).Article 
ADS 

Google Scholar 
Torrescano-Valle, N. & Islebe, G. A. Holocene paleoecology, climate history and human influence in the southwestern Yucatán Peninsula. Rev. Palaeobot. Palynol. 217, 1–8 (2015).Article 

Google Scholar 
Anselmetti, F. S., Hodell, D. A., Ariztegui, D., Brenner, M. & Rosenmeier, M. F. Quantification of soil erosion rates related to ancient Maya deforestation. Geology 35, 915–918 (2007).Article 
ADS 

Google Scholar 
Beach, T. et al. A review of human and natural changes in Maya Lowland wetlands over the Holocene. Quat. Sci. Rev. 28, 1710–1724 (2009).Article 
ADS 

Google Scholar 
Kerr, M. T. Holocene Precipitation Variability, Prehistoric Agriculture, and Natural and Human-Set Fires in Costa Rica (University of Tennessee, 2019).
Google Scholar 
Ebert, C. E., Peniche May, N., Culleton, B. J., Awe, J. J. & Kennett, D. J. Regional response to drought during the formation and decline of Preclassic Maya societies. Quat. Sci. Rev. 173, 211–235 (2017).Article 
ADS 

Google Scholar 
De la Barreda, B., Metcalfe, S. E. & Boyd, D. S. Precipitation regionalization, anomalies and drought occurrence in the Yucatan Peninsula, Mexico. Int. J. Climatol. 40, 4541–4555 (2020).Article 

Google Scholar 
Islebe, G. A. et al. Holocene Paleoecology and Paleoclimatology of south and south-eastern Mexico: A palynological approach. in Mexico´s Environmental Holocene and Anthropocene History (eds. Torrescano-Valle, N., Islebe, G. A. & Roy, P.) (Springer, 2019).Tuxill, J., Reyes, L. A., Moreno, L. L., Uicab, V. C. & Jarvis, D. I. All maize is not equal: Maize variety choices and Mayan foodways in rural Yucatan, Mexico. in Pre-Columbian Foodways: Interdisciplinary Approaches to Food, Culture, and Markets in Ancient Mesoamerica (eds. Staller, J. & Carrasco, M.) 467–486 (Springer, 2010).Torrescano-Valle, N., Ramírez-Barajas, P. J., Islebe, G. A., Vela-Pelaez, A. A. & Folan, W. J. Human influence versus natural climate variability. in The Holocene and Anthropocene Environmental History of Mexico: A Paleoecological Approach on Mesoamerica (eds. Torrescano-Valle, N., Islebe, G. A. & Roy, P. D.). 171–194 (Springer, 2019).Faegri, K. & Iversen, J. Textbook of Pollen Analysis (Wiley, 1989).
Google Scholar 
Ford, A. The Maya forest: A domesticated landscape. in The Maya World (eds. Hutson, S. R. & Ardren, T.). 519–539 (Routledge, 2020).Fedick, S. L. & Santiago, L. S. Large variation in availability of Maya food plant sources during ancient droughts. Proc. Natl. Acad. Sci. USA 119, 2115657118 (2022).Article 

Google Scholar 
Puleston, D. E. The role of ramón in Maya subsistence. in Maya Subsistence. 353–366 (Elsevier, 1982).Atran, S. et al. Itza Maya tropical agro-forestry [and comments and replies]. Curr. Anthropol. 34, 633–700 (1993).Article 

Google Scholar 
Dussol, L., Elliott, M., Michelet, D. & Nondédéo, P. Ancient Maya sylviculture of breadnut (Brosimum alicastrum Sw.) and sapodilla (Manilkara zapota (L.) P. Royen) at Naachtun (Guatemala): A reconstruction based on charcoal analysis. Quat. Int. 457, 29–42 (2017).Ebel, R., de Jesús Méndez Aguilar, M. & Putnam, H. R. Milpa: One sister got climate-sick. The impact of climate change on traditional Maya farming systems. Int. J. Sociol. Agric. Food (Online) 24, 175–199 (2018).
Google Scholar 
Hernández-González, O. & Vergara-Yoisura, S. Studies on the productivity of Brosimum alicastrum a tropical tree used for animal feed in the Yucatan Peninsula. Bothalia 22, 7 (2014).
Google Scholar 
Martínez-Ruiz, N. del R. & Larqué-Saavedra, A. Semilla de Ramón. in Alimentos Vegetales Autóctonos Iberoamericanos Subutilizados (eds. Sonia, S.-A. & Álvarez-Parrilla, E.). 177–192 (Fabro Editores, 2018).Hatfield, J. L. & Dold, C. Water-use efficiency: Advances and challenges in a changing climate. Front. Plant Sci. 10, 103 (2019).Article 

Google Scholar 
Basso, B. & Ritchie, J. T. Evapotranspiration in high-yielding maize and under increased vapor pressure deficit in the US Midwest. Agric. Environ. Lett. 3, 170039 (2018).Article 

Google Scholar 
Gregory, P. J., Simmonds, L. P. & Pilbeam, C. J. Soil type, climatic regime, and the response of water use efficiency to crop management. Agron. J. 92, 814–820 (2000).Article 

Google Scholar 
Moy, C. M., Seltzer, G. O., Rodbell, D. T. & Anderson, D. M. Variability of El Niño/Southern Oscillation activity at millennial timescales during the Holocene epoch. Nature 420, 162–165 (2002).Article 
ADS 
CAS 

Google Scholar 
Revelle, W. psych: Procedures for Psychological, Psychometric, and Personality Research. R package at https://CRAN.R-project.org/package=psych (2022).Wickham, H. & Bryan, J. readxl: Read Excel Files. R package at https://readxl.tidyverse.org/ (2022).Wei, T. et al. Package ‘corrplot’. Statistician 56, e24 (2017).
Google Scholar 
QGIS Development Team. QGIS Geographic Information System. QGIS Association at https://www.qgis.org (2022)Instituto Nacional de Estadistica Geographia e Informatica (INEGI). 1:1000000 Merida, Carta de Precipitacion. Merida, Yucatán, Mexico (1981). More