Philippa Kaur

Authors and AffiliationsJacob Blaustein Center for Scientific Cooperation, The Jacob Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Midreshet Ben-Gurion, IsraelGopal MuraliMitrani Department of Desert Ecology, The Swiss Institute for Dryland Environments and Energy Research, The Jacob Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Midreshet Ben-Gurion, IsraelGopal Murali & Uri RollDepartment F.-A. Forel for Aquatic and Environmental Sciences, Faculty of Science, University of Geneva, Geneva, SwitzerlandTakuya IwamuraDepartment of Forest Ecosystems and Society, College of Forestry, Oregon State University, Corvallis, OR, USATakuya IwamuraSchool of Zoology, Tel Aviv University, Tel Aviv, IsraelShai MeiriThe Steinhardt Museum of Natural History, Tel Aviv University, Tel Aviv, IsraelShai MeiriAuthorsGopal MuraliTakuya IwamuraShai MeiriUri RollCorresponding authorCorrespondence to
Gopal Murali. More

In 2008, I had just begun volunteering at Equilibrio Azul — a non-profit marine-research and -conservation organization based in Quito — when colleagues discovered a hawksbill sea turtle (Eretmochelys imbricata) nesting at La Playita beach in Ecuador. The eastern Pacific population of hawksbill sea turtles is one of the most endangered in the world and was considered functionally extinct in the region before this turtle and others were observed.That discovery was a tipping point for hawksbill research in Ecuador and throughout the Pacific Ocean. Since 2008, we’ve found about 20 nests each year at La Playita, and one season, we documented 50.We have tagged 11 adult females with satellite transmitters. Previously, most of our understanding of these turtles had been based on observations in the Caribbean, where the reptiles are strictly coral-reef dwellers. But Ecuador’s reefs are mostly rocky, with patches of coral, and we were surprised to see females migrate south to mangroves, mainly for food.
Women in science
In this image, we have just attached a transmitter to a baby turtle — a first for hawksbill turtles this young and in the eastern Pacific region. We did not know much about hawksbills at this young age. It is tricky working with baby turtles, because they grow very fast, and the transmitters, which give us location data, can easily fall off. We’ve used cement to glue the devices to the shells of six newborns so far. The longest the transmitters have lasted is three months and the shortest period was only six weeks — but the devices provided our first insights into the ‘lost years’ of sea-turtle biology.Our findings have overturned assumptions that neonates were just carried along by currents. Instead, we found that one-day-old turtles can swim against the current. They aim for a specific direction — north by northwest — as they learn to dive and swim. We tracked one-year-old hawksbills to Costa Rican waters, a journey of roughly 2,000 kilometres, before we lost their signal.Cristina Miranda is a scientific coordinator at Equilibrio Azul in Quito, Ecuador. Interview by Virginia Gewin. More

Pachauri, R. K. et al. Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC, 2014).Fu, F. X., Warner, M. E., Zhang, Y., Feng, Y. & Hutchins, D. A. Effects of Increased temperature and CO2 on photosynthesis, growth, and elemental ratios in marine Synechococcus and Prochlorococcus (cyanobacteria) 1. J. Phycol. 43, 485–496 (2007).
Google Scholar 
Schippers, P., Lürling, M. & Scheffer, M. Increase of atmospheric CO2 promotes phytoplankton productivity. Ecol. Lett. 7, 446–451 (2004).
Google Scholar 
Raven, J. A., Gobler, C. J. & Hansen, P. J. Dynamic CO2 and pH levels in coastal, estuarine, and inland waters: Theoretical and observed effects on harmful algal blooms. Harmful Algae 91, 101594 (2020).CAS 

Google Scholar 
Gobler, C. J. et al. Ocean warming since 1982 has expanded the niche of toxic algal blooms in the North Atlantic and North Pacific oceans. Proc. Natl Acad. Sci. 114, 4975–4980 (2017).CAS 

Google Scholar 
Rost, B., Richter, K. U., Riebesell, U. & Hansen, P. J. Inorganic carbon acquisition in red tide dinoflagellates. Plant, Cell Environ. I 29, 810–822 (2006).CAS 

Google Scholar 
Honjo, T. Harmful Red Tides of Heterosigma akashiwo. NOAA Technical Report NMFS. 111, 27–32 (1992).Rensel, J. J. & Haigh, N. Fraser river sockeye salmon marine survival decline and harmful blooms of Heterosigma akashiwo. Harmful Algae 10, 98–115 (2010).
Google Scholar 
Herndon, J. & Cochlan, W. P. Nitrogen utilization by the raphidophyte Heterosigma akashiwo: growth and uptake kinetics in laboratory cultures. Harmful Algae 6, 260–270 (2007).
Google Scholar 
Haley, S. T., Alexander, H., Juhl, A. R. & Dyhrman, S. T. Transcriptional response of the harmful raphidophyte Heterosigma akashiwo to nitrate and phosphate stress. Harmful Algae 68, 258–270 (2017).CAS 

Google Scholar 
Wang, Z.-h, Liang, Y. & Kang, W. Utilization of dissolved organic phosphorus by different groups of phytoplankton taxa. Harmful Algae 12, 113–118 (2011).CAS 

Google Scholar 
Ji, N. et al. Metatranscriptome analysis reveals environmental and diel regulation of a Heterosigma akashiwo (raphidophyceae) bloom. Environ. Microbiol. 20, 1078–1094 (2018).CAS 

Google Scholar 
Zhang, H. et al. Functional differences in the blooming phytoplankton Heterosigma akashiwo and Prorocentrum donghaiense revealed by comparative metaproteomics. Appl. Environ. Microbiol. 85, e01425–01419 (2019).CAS 

Google Scholar 
Redfield, A. C. The biological control of chemical factors in the environment. Am. Scientist 46, 230A–221 (1958).
Google Scholar 
Liefer, J. D. et al. The macromolecular basis of phytoplankton C: N: P under nitrogen starvation. Front. Microbiol. 10, 763 (2019).
Google Scholar 
Matsumoto, K., Tanioka, T. & Rickaby, R. Linkages between dynamic phytoplankton C: N: P and the ocean carbon cycle under climate change. Oceanography 33, 44–52 (2020).
Google Scholar 
Thrane, J. E., Hessen, D. O. & Andersen, T. Plasticity in algal stoichiometry: Experimental evidence of a temperature‐induced shift in optimal supply N: P ratio. Limnol. Oceanogr. 62, 1346–1354 (2017).CAS 

Google Scholar 
Toseland, A. et al. The impact of temperature on marine phytoplankton resource allocation and metabolism. Nat. Clim. Change 3, 979–984 (2013).CAS 

Google Scholar 
Mittler, R., Finka, A. & Goloubinoff, P. How do plants feel the heat? Trends Biochem. Sci. 37, 118–125 (2012).CAS 

Google Scholar 
Dingman, J. E. & Lawrence, J. E. Heat-stress-induced programmed cell death in Heterosigma akashiwo (Raphidophyceae). Harmful Algae 16, 108–116 (2012).
Google Scholar 
Whitten, S. T., García-Moreno E, B. & Hilser, V. J. Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins. Proc. Natl Acad. Sci. 102, 4282–4287 (2005).CAS 

Google Scholar 
Casey, J. R., Grinstein, S. & Orlowski, J. Sensors and regulators of intracellular pH. Nat. Rev. Mol. Cell Biol. 11, 50–61 (2010).CAS 

Google Scholar 
Kim, H., Spivack, A. J. & Menden-Deuer, S. pH alters the swimming behaviors of the raphidophyte Heterosigma akashiwo: implications for bloom formation in an acidified ocean. Harmful Algae 26, 1–11 (2013).CAS 

Google Scholar 
Hennon, G. M., Williamson, O. M., Limón, M. D. H., Haley, S. T. & Dyhrman, S. T. Non-linear physiology and gene expression responses of harmful alga Heterosigma akashiwo to rising CO2. Protist 170, 38–51 (2019).CAS 

Google Scholar 
Xu, H., Jaynes, J. & Ding, X. Combining two-level and three-level orthogonal arrays for factor screening and response surface exploration. Statistica Sin. 24, 269–289 (2014).
Google Scholar 
Boyd, P. W. & Ellwood, M. J. The biogeochemical cycle of iron in the ocean. Nat. Geosci. 3, 675–682 (2010).CAS 

Google Scholar 
Sterner, R. W. & Elser, J. J. in Ecological Stoichiometry (Princeton university press, 2002).Liu, H. C., Liao, H. T. & Charng, Y. Y. The role of class A1 heat shock factors (HSFA1s) in response to heat and other stresses in Arabidopsis. Plant Cell Environ. 34, 738–751 (2011).CAS 

Google Scholar 
Geider, R. J. & La Roche, J. J. Redfield revisited: variability of C [ratio] N [ratio] P in marine microalgae and its biochemical basis. Eur. J. Phycol. 37, 1–17 (2002).
Google Scholar 
Loladze, I. & Elser, J. J. The origins of the Redfield nitrogen‐to‐phosphorus ratio are in a homoeostatic protein‐to‐rRNA ratio. Ecol. Lett. 14, 244–250 (2011).
Google Scholar 
Hennige, S. J., Coyne, K. J., MacIntyre, H., Liefer, J. & Warner, M. E. The photobiology of Heterosigma akashiwo. Photoacclimation, diurnal periodicity, and its ability to rapidly exploit exposure to high light. J. Phycol. 49, 349–360 (2013).CAS 

Google Scholar 
Collier, J. L. & Grossman, A. A small polypeptide triggers complete degradation of light‐harvesting phycobiliproteins in nutrient‐deprived cyanobacteria. EMBO J. 13, 1039–1047 (1994).CAS 

Google Scholar 
Gordillo, F. J., Jimenez, C., Figueroa, F. L. & Niell, F. X. Influence of elevated CO2 and nitrogen supply on the carbon assimilation performance and cell composition of the unicellular alga Dunaliella viridis. Physiologia Plant. 119, 513–518 (2003).CAS 

Google Scholar 
Satoh, E., Watanabe, M. M. & Fujii, T. Photoperiodic regulation of cell division and chloroplast replication in Heterosigma akashiwo. Plant Cell Physiol. 28, 1093–1099 (1987).
Google Scholar 
Ashworth, J. et al. Genome-wide diel growth state transitions in the diatom Thalassiosira pseudonana. Proc. Natl Acad. Sci. 110, 7518–7523 (2013).CAS 

Google Scholar 
Thangaraj, S. & Sun, J. J. E. M. Transcriptomic reprogramming of the oceanic diatom Skeletonema dohrnii under warming ocean and acidification. Environ. Microbiol. 23, 980–995 (2021).CAS 

Google Scholar 
Nakajima, K., Tanaka, A. & Matsuda, Y. SLC4 family transporters in a marine diatom directly pump bicarbonate from seawater. Proc. Natl Acad. Sci. 110, 1767–1772 (2013).CAS 

Google Scholar 
Kranz, S. A. et al. Low temperature reduces the energetic requirement for the CO2 concentrating mechanism in diatoms. N. Phytologist 205, 192–201 (2015).CAS 

Google Scholar 
Ralston, A. & Shaw, K. Gene expression regulates cell differentiation. Nat. Educ. 1, 127–131 (2008).
Google Scholar 
Lobo, I. Environmental influences on gene expression. Nat. Educ. 1, 39 (2008).
Google Scholar 
Suzuki, N. et al. Respiratory burst oxidases: the engines of ROS signaling. Curr. Opin. Plant Biol. 14, 691–699 (2011).CAS 

Google Scholar 
Saidi, Y., Finka, A. & Goloubinoff, P. Heat perception and signalling in plants: a tortuous path to thermotolerance. N. Phytologist 190, 556–565 (2011).CAS 

Google Scholar 
Saidi, Y. et al. The heat shock response in moss plants is regulated by specific calcium-permeable channels in the plasma membrane. Plant Cell 21, 2829–2843 (2009).CAS 

Google Scholar 
Zhang, W. et al. Molecular and genetic evidence for the key role of AtCaM3 in heat-shock signal transduction in Arabidopsis. Plant Physiol. 149, 1773–1784 (2009).CAS 

Google Scholar 
Li, S. et al. Functional characterization of Arabidopsis thaliana WRKY39 in heat stress. Mol. Cells 29, 475–483 (2010).CAS 

Google Scholar 
Sangwan, V., Örvar, B. L., Beyerly, J., Hirt, H. & Dhindsa, R. S. Opposite changes in membrane fluidity mimic cold and heat stress activation of distinct plant MAP kinase pathways. Plant J. 31, 629–638 (2002).CAS 

Google Scholar 
Reddy, A. S., Ali, G. S., Celesnik, H. & Day, I. S. Coping with stresses: roles of calcium-and calcium/calmodulin-regulated gene expression. Plant Cell 23, 2010–2032 (2011).CAS 

Google Scholar 
Meiri, D. & Breiman, A. J. Arabidopsis ROF1 (FKBP62) modulates thermotolerance by interacting with HSP90. 1 and affecting the accumulation of HsfA2‐regulated sHSPs. Plant J. 59, 387–399 (2009).CAS 

Google Scholar 
Mishkind, M., Vermeer, J. E., Darwish, E. & Munnik, T. J. Heat stress activates phospholipase D and triggers PIP2 accumulation at the plasma membrane and nucleus. Plant J. 60, 10–21 (2009).CAS 

Google Scholar 
Zheng, S. Z. et al. Phosphoinositide‐specific phospholipase C9 is involved in the thermotolerance of Arabidopsis. Plant J. 69, 689–700 (2012).CAS 

Google Scholar 
Pincus, D. et al. BiP binding to the ER-stress sensor Ire1 tunes the homeostatic behavior of the unfolded protein response. PLoS Biol. 8, e1000415 (2010).
Google Scholar 
Sugio, A., Dreos, R., Aparicio, F. & Maule, A. J. The cytosolic protein response as a subcomponent of the wider heat shock response in Arabidopsis. Plant Cell 21, 642–654 (2009).CAS 

Google Scholar 
Vasseur, F., Pantin, F. & Vile, D. J. P. Cell & Environment. Changes in light intensity reveal a major role for carbon balance in Arabidopsis responses to high temperature. Plant Cell Environ. 34, 1563–1576 (2011).CAS 

Google Scholar 
Paroutis, P., Touret, N. & Grinstein, S. The pH of the secretory pathway: measurement, determinants, and regulation. J. Physiol. 19, 207–215 (2004).CAS 

Google Scholar 
Forgac, M. Vacuolar ATPases: rotary proton pumps in physiology and pathophysiology. Nat. Rev. Mol. Cell Biol. 8, 917–929 (2007).CAS 

Google Scholar 
Cipriano, D. J. et al. Structure and regulation of the vacuolar ATPases. Biochim. et. Biophys. Acta -Bioenerg. 1777, 599–604 (2008).CAS 

Google Scholar 
Abad, M. F. C., Di Benedetto, G., Magalhães, P. J., Filippin, L. & Pozzan, T. Mitochondrial pH monitored by a new engineered green fluorescent protein mutant. J. Biol. Chem. 279, 11521–11529 (2004).CAS 

Google Scholar 
McCORMACK, J. G., Halestrap, A. P. & Denton, R. M. Role of calcium ions in regulation of mammalian intramitochondrial metabolism. Physiol. Rev. 70, 391–425 (1990).CAS 

Google Scholar 
Garlid, K. D., Sun, X., Paucek, P. & Woldegiorgis, G. in Methods in enzymology Vol. 260 331–348 (Elsevier, 1995).Yamada, E. W. & Huzel, N. J. J. B. Calcium-binding ATPase inhibitor protein of bovine heart mitochondria. Role in ATP synthesis and effect of calcium. Biochemistry 28, 9714–9718 (1989).CAS 

Google Scholar 
Moreno-Sánchez, R. Inhibition of oxidative phosphorylation by a Ca2+-induced diminution of the adenine nucleotide translocator. Biochim. et. Biophys. Acta -Bioenerg. 724, 278–285 (1983).
Google Scholar 
Matsuyama, S., Llopis, J., Deveraux, Q. L., Tsien, R. Y. & Reed, J. Changes in intramitochondrial and cytosolic pH: early events that modulate caspase activation during apoptosis. Nat. Cell Biol. 2, 318–325 (2000).CAS 

Google Scholar 
Sunda, W. G., Price, N. M. & Morel, F. M. Trace metal ion buffers and their use in culture studies. Algal Cultur. Tech. 4, 35–63 (2005).
Google Scholar 
Sun, J. et al. Effects of changing pCO2 and phosphate availability on domoic acid production and physiology of the marine harmful bloom diatom Pseudo‐nitzschia multiseries. Limnol. Oceanogr. 56, 829–840 (2011).CAS 

Google Scholar 
Pierrot, D., Lewis, E. & Wallace, D. J. MS Excel Program Developed for CO2 System Calculations ORNL/CDIAC‐105 (US Dept. of Energy, Oak Ridge, TN, 2006).Wilbur, K. M. & Anderson, N. G. Electrometric and colorimetric determination of carbonic anhydrase. J. Biol. Chem. 176, 147–154 (1948).CAS 

Google Scholar 
Solórzano, L. & Sharp, J. H. Determination of total dissolved phosphorus and particulate phosphorus in natural waters 1. Limnol. Oceanogr. 25, 754–758 (1980).
Google Scholar 
Myklestad, S. M., Skånøy, E. & Hestmann, S. J. Sensitive and rapid method for analysis of dissolved mono-and polysaccharides in seawater. Mar. Chem. 56, 279–286 (1997).CAS 

Google Scholar 
Pakulski, J. D. & Benner, R. J. An improved method for the hydrolysis and MBTH analysis of dissolved and particulate carbohydrates in seawater. Mar. Chem. 40, 143–160 (1992).CAS 

Google Scholar 
Folch, J. & Lees, M. & Sloane Stanley, G. H. A simple method for the isolation and purification of total lipids from animal tissues. J. Biol. Chem. 226, 497–509 (1957).CAS 

Google Scholar 
Pande, S., Khan, R. P. & Venkitasubramanian, T. Microdetermination of lipids and serum total fatty acids. Anal. Biochem. 6, 415–423 (1963).CAS 

Google Scholar 
Lowry, O., Rosebrough, N., Farr, A. L. & Randall, R. J. Protein measurement with the Folin phenol reagent. J. Biol. Chem. 193, 265–275 (1951).CAS 

Google Scholar 
Berdalet, E., Roldán, C., Olivar, M. P. & Lysnes, K. Quantifying RNA and DNA in planktonic organisms with SYBR Green II and nucleases. Part A. Optimisation of the assay. Sci. Mar. 69, 1–16 (2005).CAS 

Google Scholar 
Chomoczynski, P. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloro-form extraction. Anal. Biochem. 162, 156–159 (1987).
Google Scholar 
Sañudo-Wilhelmy, S. A. et al. The impact of surface-adsorbed phosphorus on phytoplankton Redfield stoichiometry. Phycol. Res. 432, 897–901 (2004).
Google Scholar 
Dyhrman, S. T. Nutrients and their acquisition: phosphorus physiology in microalgae. Physiol. Microalgae 155–183 (2016).Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).CAS 

Google Scholar 
Waterhouse, R. M. et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543–548 (2018).CAS 

Google Scholar 
Pruitt, K. D., Tatusova, T. & Maglott, D. R. J. N. A. R. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 35, D61–D65 (2007).CAS 

Google Scholar 
Kanehisa, M. & Goto, S. J. N. A. R. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000).CAS 

Google Scholar 
Mitchell, A. L. et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 47, D351–D360 (2019).CAS 

Google Scholar 
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21 (2014).
Google Scholar 
Wang, L., Feng, Z., Wang, X., Wang, X. & Zhang, X. J. B. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26, 136–138 (2010).
Google Scholar  More

Saatchi, S. S. et al. Benchmark map of forest carbon stocks in tropical regions across three continents. Proc. Natl Acad. Sci. 108, 9899–9904 (2011).Article 
CAS 

Google Scholar 
Baccini, A. et al. Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nat. Clim. Change 2, 182–185 (2012).Article 
CAS 

Google Scholar 
Xu, L. et al. Changes in global terrestrial live biomass over the 21st century. Sci. Adv. 7, eabe9829 (2021).Article 
CAS 

Google Scholar 
Betts, R. A. et al. The role of ecosystem-atmosphere interactions in simulated Amazonian precipitation decrease and forest dieback under global climate warming. Theor. Appl. Climatol. 78, 157–175 (2004).Article 

Google Scholar 
Cox, P. M. et al. Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability. Nature 494, 341–344 (2013).Article 
CAS 

Google Scholar 
Rammig, A. et al. Estimating the risk of Amazonian forest dieback. N. Phytol. 187, 694–706 (2010).Article 
CAS 

Google Scholar 
Huntingford, C. et al. Towards quantifying uncertainty in predictions of Amazon ‘dieback’. Philos. Trans. R. Soc. B Biol. Sci. 363, 1857–1864 (2008).Article 

Google Scholar 
Galbraith, D. et al. Multiple mechanisms of Amazonian forest biomass losses in three dynamic global vegetation models under climate change. N. Phytol. 187, 647–665 (2010).Article 

Google Scholar 
Kumar, D., Pfeiffer, M., Gaillard, C., Langan, L. & Scheiter, S. Climate change and elevated CO2 favor forest over savanna under different future scenarios in South Asia. Biogeosciences 18, 2957–2979 (2021).Article 
CAS 

Google Scholar 
Huntingford, C. et al. Simulated resilience of tropical rainforests to CO2-induced climate change. Nat. Geosci. 6, 268–273 (2013).Article 
CAS 

Google Scholar 
Brienen, R. J. W. et al. Forest carbon sink neutralized by pervasive growth-lifespan trade-offs. Nat. Commun. 11, 4241 (2020).Article 
CAS 

Google Scholar 
Koch, A., Hubau, W. & Lewis, S. L. Earth system models are not capturing present-day tropical forest carbon dynamics. Earths Future 9, e2020EF001874 (2021).Article 
CAS 

Google Scholar 
Negrón-Juárez, R. I., Koven, C. D., Riley, W. J., Knox, R. G. & Chambers, J. Q. Observed allocations of productivity and biomass, and turnover times in tropical forests are not accurately represented in CMIP5 Earth system models. Environ. Res. Lett. 10, 064017 (2015).Article 

Google Scholar 
Fleischer, K. et al. Amazon forest response to CO2 fertilization dependent on plant phosphorus acquisition. Nat. Geosci. 12, 736–741 (2019).Article 
CAS 

Google Scholar 
Terrer, C. et al. Nitrogen and phosphorus constrain the CO2 fertilization of global plant biomass. Nat. Clim. Change 9, 684–689 (2019).Article 
CAS 

Google Scholar 
Malhi, Y. et al. Exploring the likelihood and mechanism of a climate-change-induced dieback of the Amazon rainforest. Proc. Natl Acad. Sci. 106, 20610–20615 (2009).Article 
CAS 

Google Scholar 
Zelazowski, P., Malhi, Y., Huntingford, C., Sitch, S. & Fisher, J. B. Changes in the potential distribution of humid tropical forests on a warmer planet. Philos. Trans. R. Soc. Math. Phys. Eng. Sci. 369, 137–160 (2011).
Google Scholar 
Huang, L. et al. Drought dominates the interannual variability in global terrestrial net primary production by controlling semi-arid ecosystems. Sci. Rep. 6, 24639 (2016).Article 
CAS 

Google Scholar 
Castanho, A. D. A. et al. Potential shifts in the aboveground biomass and physiognomy of a seasonally dry tropical forest in a changing climate. Environ. Res. Lett. 15, 034053 (2020).Article 

Google Scholar 
Santoro, M. & Cartus, O. ESA Biomass Climate Change Initiative (Biomass_cci): global datasets of forest above-ground biomass for the years 2010, 2017 and 2018, v3. NERC EDS Centre for Environmental Data Analysis https://doi.org/10.5285/5f331c418e9f4935b8eb1b836f8a91b8 (2021).Baccini, A. et al. Tropical forests are a net carbon source based on aboveground measurements of gain and loss. Science 358, 230–234 (2017).Article 
CAS 

Google Scholar 
Harris, N. L. et al. Global maps of twenty-first century forest carbon fluxes. Nat. Clim. Change 11, 234–240 (2021).Article 

Google Scholar 
Gatti, L. V. et al. Amazonia as a carbon source linked to deforestation and climate change. Nature 595, 388–393 (2021).Article 
CAS 

Google Scholar 
Qin, Y. et al. Carbon loss from forest degradation exceeds that from deforestation in the Brazilian Amazon. Nat. Clim. Change 11, 442–448 (2021).Article 

Google Scholar 
Brienen, R. J. W. et al. Long-term decline of the Amazon carbon sink. Nature 519, 344–348 (2015).Article 
CAS 

Google Scholar 
Hubau, W. et al. Asynchronous carbon sink saturation in African and Amazonian tropical forests. Nature 579, 80–87 (2020).Article 
CAS 

Google Scholar 
Phillips, O. L. et al. Drought sensitivity of the amazon rainforest. Science 323, 1344–1347 (2009).Article 
CAS 

Google Scholar 
Ross, C. W. et al. Woody-biomass projections and drivers of change in sub-Saharan Africa. Nat. Clim. Change 11, 449–455 (2021).Article 

Google Scholar 
Larjavaara, M., Lu, X., Chen, X. & Vastaranta, M. Impact of rising temperatures on the biomass of humid old-growth forests of the world. Carbon Balance Manag. 16, 31 (2021).Article 

Google Scholar 
Romps, D. M., Seeley, J. T., Vollaro, D. & Molinari, J. Projected increase in lightning strikes in the United States due to global warming. Science 346, 851–854 (2014).Article 
CAS 

Google Scholar 
Gora, E. M., Bitzer, P. M., Burchfield, J. C., Gutierrez, C. & Yanoviak, S. P. The contributions of lightning to biomass turnover, gap formation and plant mortality in a tropical forest. Ecology 102, e03541 (2021).Article 

Google Scholar 
Magnabosco Marra, D. et al. Windthrows control biomass patterns and functional composition of Amazon forests. Glob. Change Biol. 24, 5867–5881 (2018).Article 

Google Scholar 
Negrón-Juárez, R. I. et al. Windthrow variability in central amazonia. Atmosphere 8, 28 (2017).Article 

Google Scholar 
Silva Junior, C. H. L. et al. Persistent collapse of biomass in Amazonian forest edges following deforestation leads to unaccounted carbon losses. Sci. Adv. 6, 40 (2020).Article 

Google Scholar 
Yin, Y. et al. Fire decline in dry tropical ecosystems enhances decadal land carbon sink. Nat. Commun. 11, 1900 (2020).Article 
CAS 

Google Scholar 
Koch, A. & Kaplan, J. O. Tropical forest restoration under future climate change. Nat. Clim. Change 12, 279–283 (2022).Article 

Google Scholar 
Wang, S. et al. Recent global decline of CO2 fertilization effects on vegetation photosynthesis. Science 370, 1295–1300 (2020).Article 
CAS 

Google Scholar 
Case, M. F. & Staver, A. C. Fire prevents woody encroachment only at higher-than-historical frequencies in a South African savanna. J. Appl. Ecol. 54, 955–962 (2017).Article 
CAS 

Google Scholar 
Mau, A. C., Reed, S. C., Wood, T. E. & Cavaleri, M. A. Temperate and tropical forest canopies are already functioning beyond their thermal thresholds for photosynthesis. Forests 9, 47 (2018).Article 

Google Scholar 
Kolby Smith, W. et al. Large divergence of satellite and Earth system model estimates of global terrestrial CO2 fertilization. Nat. Clim. Change 6, 306–310 (2016).Article 
CAS 

Google Scholar 
Martens, C. et al. Large uncertainties in future biome changes in Africa call for flexible climate adaptation strategies. Glob. Change Biol. 27, 340–358 (2021).Article 
CAS 

Google Scholar 
Doughty, C. E. & Goulden, M. L. Are tropical forests near a high temperature threshold?. J. Geophys. Res. Biogeosciences 113, G00B07 (2008).Article 

Google Scholar 
Doughty, C. E. & Goulden, M. L. Seasonal patterns of tropical forest leaf area index and CO2 exchange. J. Geophys. Res. Biogeosciences 113, G00B06 (2008).Article 

Google Scholar 
Langenbrunner, B., Pritchard, M. S., Kooperman, G. J. & Randerson, J. T. Why does amazon precipitation decrease when tropical forests respond to increasing CO2? Earths Future 7, 450–468 (2019).Article 

Google Scholar 
Harris, I., Osborn, T. J., Jones, P. & Lister, D. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Sci. Data 7, 109 (2020).Article 

Google Scholar 
Maurer, E. P., Brekke, L., Pruitt, T. & Duffy, P. B. Fine-resolution climate projections enhance regional climate change impact studies. EOS Trans. Am. Geophys. Union 88, 504–504 (2007).Article 

Google Scholar 
Reclamation. Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections: Release of Hydrology Projections, Comparison with preceding Information, and Summary of User Needs. https://gdo-dcp.ucllnl.org/downscaled_cmip_projections/techmemo/BCSD5HydrologyMemo.pdf (2014).Silva de Miranda, P. L. et al. Using tree species inventories to map biomes and assess their climatic overlaps in lowland tropical South America. Glob. Ecol. Biogeogr. 27, 899–912 (2018).Article 

Google Scholar 
Beguería, S., Vicente-Serrano, S. M., Reig, F. & Latorre, B. Standardized precipitation evapotranspiration index (SPEI) revisited: parameter fitting, evapotranspiration models, tools, datasets and drought monitoring. Int. J. Climatol. 34, 3001–3023 (2014).Article 

Google Scholar 
Wright, M. N. & Ziegler, A. ranger: a fast implementation of random forests for high dimensional data in C++ and R. J. Stat. Softw. 77, 1–17 (2017).Article 

Google Scholar 
Middleton, N., Thomas, D. & UNEP. World Atlas of Desertification (Arnold, 1997).Staver, A. C., Archibald, S. & Levin, S. A. The global extent and determinants of Savanna and forest as alternative biome states. Science 334, 230–232 (2011).Article 
CAS 

Google Scholar 
ESRI Data & Maps. World Continents Version 10.3. (2015).Uribe, M. R. et al. Net loss of biomass predicted for tropical biomes in a changing climate. Dryad https://doi.org/10.7280/D1D124 (2023). More

Hibbing, M. E., Fuqua, C., Parsek, M. R. & Peterson, S. B. Bacterial competition: surviving and thriving in the microbial jungle. Nat. Rev. Microbiol 8, 15–25 (2010).CAS 

Google Scholar 
Cragg, G. M. & Newman, D. J. Natural products: a continuing source of novel drug leads. Biochim. Biophys. Acta 1830, 3670–3695 (2013).CAS 

Google Scholar 
Blin, K., Kim, H. U., Medema, M. H. & Weber, T. Recent development of antiSMASH and other computational approaches to mine secondary metabolite biosynthetic gene clusters. Brief. Bioinform. 20, 1103–1113 (2019).CAS 

Google Scholar 
Paoli, L. et al. Biosynthetic potential of the global ocean microbiome. Nature 607, 111–118 (2022).Gavriilidou, A. et al. Compendium of specialized metabolite biosynthetic diversity encoded in bacterial genomes. Nat. Microbiol. 7, 726–735 (2022).CAS 

Google Scholar 
Scherlach, K. & Hertweck, C. Triggering cryptic natural product biosynthesis in microorganisms. Org. Biomol. Chem. 7, 1753–1760 (2009).CAS 

Google Scholar 
Gilly, W. F., Beman, J. M., Litvin, S. Y. & Robison, B. H. Oceanographic and biological effects of shoaling of the oxygen minimum zone. Ann. Rev. Mar. Sci. 5, 393–420 (2013).
Google Scholar 
Schmidtko, S., Stramma, L. & Visbeck, M. Decline in global oceanic oxygen content during the past five decades. Nature 542, 335–339 (2017).ADS 
CAS 

Google Scholar 
Naqvi, S. W. A. et al. Marine hypoxia/anoxia as a source of CH 4 and N 2 O. Biogeosciences 7, 2159–2190 (2010).ADS 
CAS 

Google Scholar 
Scranton, M. I., Sayles, F. L., Bacon, M. P. & Brewer, P. G. Temporal changes in the hydrography and chemistry of the Cariaco Trench. Deep-Sea Res. Part A. Oceanogr. Res. Pap. 34, 945–963 (1987).ADS 
CAS 

Google Scholar 
Taylor, G. T. et al. Chemoautotrophy in the redox transition zone of the Cariaco Basin: a significant midwater source of organic carbon production. Limnol. Oceanogr. 46, 148–163 (2001).ADS 
CAS 

Google Scholar 
Scranton, M. I., Astor, Y., Bohrer, R., Ho, T.-Y. & Muller-Karger, F. Controls on temporal variability of the geochemistry of the deep Cariaco Basin. Deep-Sea Res. I Oceanogr. Res. Pap. 48, 1605–1625 (2001).ADS 
CAS 

Google Scholar 
Scranton, M. I. et al. Interannual and subdecadal variability in the nutrient geochemistry of the Cariaco Basin. Oceanography 27, 148–159 (2014).
Google Scholar 
Dalsgaard, T., Thamdrup, B., Farías, L. & Revsbech, N. P. Anammox and denitrification in the oxygen minimum zone of the eastern South Pacific. Limnol. Oceanogr. 57, 1331–1346 (2012).ADS 
CAS 

Google Scholar 
Canfield, D. E. et al. A cryptic sulfur cycle in oxygen-minimum–zone waters off the Chilean coast. Science 330, 1375–1378 (2010).ADS 
CAS 

Google Scholar 
Schlosser, C. et al. H 2 S events in the Peruvian oxygen minimum zone facilitate enhanced dissolved Fe concentrations. Sci. Rep. 8, 1–8 (2018).
Google Scholar 
Rapp, I. et al. Controls on redox-sensitive trace metals in the Mauritanian oxygen minimum zone. Biogeosciences 16, 4157–4182 (2019).ADS 
CAS 

Google Scholar 
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21 (2014).
Google Scholar 
Blin, K. et al. antiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res. 49, W29–W35 (2021).Edgcomb, V. P. et al. Comparison of Niskin vs. in situ approaches for analysis of gene expression in deep Mediterranean Sea water samples. Deep-Sea Res. II: Top. Stud. Oceanogr. 129, 213–222 (2016).ADS 
CAS 

Google Scholar 
Cabello-Yeves, P. J. et al. The microbiome of the Black Sea water column analyzed by shotgun and genome centric metagenomics. Environ. Microbiome 16, 1–15 (2021).
Google Scholar 
Suter, E. A., Pachiadaki, M., Taylor, G. T., Astor, Y. & Edgcomb, V. P. Free‐living chemoautotrophic and particle‐attached heterotrophic prokaryotes dominate microbial assemblages along a pelagic redox gradient. Environ. Microbiol. 20, 693–712 (2018).CAS 

Google Scholar 
Mestre, M. et al. Spatial variability of marine bacterial and archaeal communities along the particulate matter continuum. Mol. Ecol. 26, 6827–6840 (2017).CAS 

Google Scholar 
Li, J. et al. Characterization of particle-associated and free-living bacterial and archaeal communities along the water columns of the South China Sea. Biogeosciences 18, 113–133 (2021).ADS 
CAS 

Google Scholar 
Mestre, M., Borrull, E., Sala, M. M. & Gasol, J. M. Patterns of bacterial diversity in the marine planktonic particulate matter continuum. ISME J. 11, 999–1010 (2017).
Google Scholar 
Pelve, E. A., Fontanez, K. M. & DeLong, E. F. Bacterial succession on sinking particles in the ocean’s interior. Front. Microbiol. 8, 2269 (2017).
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).CAS 

Google Scholar 
Sinninghe Damsté, J. S., Rijpstra, W. I. C., Geenevasen, J. A. J., Strous, M. & Jetten, M. S. M. Structural identification of ladderane and other membrane lipids of planctomycetes capable of anaerobic ammonium oxidation (anammox). FEBS J. 272, 4270–4283 (2005).
Google Scholar 
Fuchsman, C. A., Staley, J. T., Oakley, B. B., Kirkpatrick, J. B. & Murray, J. W. Free-living and aggregate-associated Planctomycetes in the Black Sea. FEMS Microbiol. Ecol. 80, 402–416 (2012).CAS 

Google Scholar 
Scherlach, K. & Hertweck, C. Mining and unearthing hidden biosynthetic potential. Nat. Commun. 12, 1–12 (2021).
Google Scholar 
Letzel, A.-C., Pidot, S. J. & Hertweck, C. A genomic approach to the cryptic secondary metabolome of the anaerobic world. Nat. Prod. Rep. 30, 392–428 (2013).CAS 

Google Scholar 
Navarro-Muñoz, J. C. et al. A computational framework to explore large-scale biosynthetic diversity. Nat. Chem. Biol. 16, 60–68 (2020).
Google Scholar 
Mungan, M. D. et al. ARTS 2.0: feature updates and expansion of the antibiotic resistant target seeker for comparative genome mining. Nucleic Acids Res. 48, W546–W552 (2020).CAS 

Google Scholar 
Alanjary, M. et al. The antibiotic resistant target seeker (ARTS), an exploration engine for antibiotic cluster prioritization and novel drug target discovery. Nucleic Acids Res. 45, W42–W48 (2017).CAS 

Google Scholar 
Waters, A. L., Hill, R. T., Place, A. R. & Hamann, M. T. The expanding role of marine microbes in pharmaceutical development. Curr. Opin. Biotechnol. 21, 780–786 (2010).CAS 

Google Scholar 
Long, R. A. & Azam, F. Antagonistic interactions among marine pelagic bacteria. Appl. Environ. Microbiol. 67, 4975–4983 (2001).ADS 
CAS 

Google Scholar 
Graça, A. P., Calisto, R. & Lage, O. M. Planctomycetes as novel source of bioactive molecules. Front. Microbiol. 7, 1241 (2016).
Google Scholar 
Murphy, C. L. et al. Genomes of novel Myxococcota reveal severely curtailed machineries for predation and cellular differentiation. Appl. Environ. Microbiol. 87, e01706–e01721 (2021).CAS 

Google Scholar 
Pachiadaki, M. G. et al. Charting the complexity of the marine microbiome through single-cell genomics. Cell 179, 1623–1635 (2019).CAS 

Google Scholar 
Charlesworth, J. C. & Burns, B. P. Untapped resources: biotechnological potential of peptides and secondary metabolites in archaea. Archaea 2015, 282035 (2015).Wang, S. & Lu, Z. in Biocommunication of Archaea (ed. Witzany, G.) 67–101 (Springer, 2017).Castelle, C. J. et al. Genomic expansion of domain archaea highlights roles for organisms from new phyla in anaerobic carbon cycling. Curr. Biol. 25, 690–701 (2015).CAS 

Google Scholar 
McInnes, L., Healy, J. & Melville, J. Umap: uniform manifold approximation and projection for dimension reduction. Journal of Open Source Software 3, 861 (2018).Rattray, J. E. et al. A comparative genomics study of genetic products potentially encoding ladderane lipid biosynthesis. Biol. Direct 4, 1–16 (2009).
Google Scholar 
Orakov, A. et al. GUNC: detection of chimerism and contamination in prokaryotic genomes. Genome Biol. 22, 1–19 (2021).
Google Scholar 
Choudoir, M. J., Pepe-Ranney, C. & Buckley, D. H. Diversification of secondary metabolite biosynthetic gene clusters coincides with lineage divergence in Streptomyces. Antibiotics 7, 12 (2018).
Google Scholar 
Li, Y. & Rebuffat, S. The manifold roles of microbial ribosomal peptide–based natural products in physiology and ecology. J. Biol. Chem. 295, 34–54 (2020).CAS 

Google Scholar 
Ma, L. & Payne, S. M. AhpC is required for optimal production of enterobactin by Escherichia coli. J. Bacteriol. 194, 6748–6757 (2012).CAS 

Google Scholar 
Davis, C. et al. The role of glutathione S-transferase GliG in gliotoxin biosynthesis in Aspergillus fumigatus. Chem. Biol. 18, 542–552 (2011).CAS 

Google Scholar 
Kautsar, S. A. et al. MIBiG 2.0: a repository for biosynthetic gene clusters of known function. Nucleic Acids Res. 48, D454–D458 (2020).
Google Scholar 
Wang, Y. et al. Phenazine-1-carboxylic acid promotes bacterial biofilm development via ferrous iron acquisition. J. Bacteriol. 193, 3606–3617 (2011).CAS 

Google Scholar 
Laursen, J. B. & Nielsen, J. Phenazine natural products: biosynthesis, synthetic analogues, and biological activity. Chem. Rev. 104, 1663–1686 (2004).CAS 

Google Scholar 
McParland, E. et al. Cycling of suspended particulate phosphorus in the redoxcline of the Cariaco Basin. Mar. Chem. 176, 64–74 (2015).CAS 

Google Scholar 
McRose, D. L. & Newman, D. K. Redox-active antibiotics enhance phosphorus bioavailability. Science 371, 1033–1037 (2021).ADS 
CAS 

Google Scholar 
Cundliffe, E. How antibiotic-producing organisms avoid suicide. Annu. Rev. Microbiol. 43, 207–233 (1989).ADS 
CAS 

Google Scholar 
Webber, M. A. & Piddock, L. J. V. The importance of efflux pumps in bacterial antibiotic resistance. J. Antimicrob. Chemother. 51, 9–11 (2003).CAS 

Google Scholar 
Vetting, M. W. et al. Pentapeptide repeat proteins. Biochemistry 45, 1–10 (2006).CAS 

Google Scholar 
Kauppinen, S., Siggaard-Andersen, M. & von Wettstein-Knowles, P. β-ketoacyl-ACP synthase I of Escherichia coli: nucleotide sequence of thefabB gene and identification of the cerulenin binding residue. Carlsberg Res. Commun. 53, 357–370 (1988).CAS 

Google Scholar 
Kloosterman, A. M., Shelton, K. E., van Wezel, G. P., Medema, M. H. & Mitchell, D. A. RRE-Finder: a genome-mining tool for class-independent RiPP discovery. mSystems 5, e00267–20 (2020).CAS 

Google Scholar 
Barry, S. M. & Challis, G. L. Mechanism and catalytic diversity of Rieske non-heme iron-dependent oxygenases. ACS Catal. 3, 2362–2370 (2013).CAS 

Google Scholar 
Benjdia, A., Balty, C. & Berteau, O. Radical SAM enzymes in the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs). Front. Chem. 5, 87 (2017).
Google Scholar 
Pandey, R. P., Parajuli, P. & Sohng, J. K. Metabolic engineering of glycosylated polyketide biosynthesis. Emerg. Top. Life Sci. 2, 389–403 (2018).CAS 

Google Scholar 
Argueta, E. A., Amoh, A. N., Kafle, P. & Schneider, T. L. Unusual non-enzymatic flavin catalysis enhances understanding of flavoenzymes. FEBS Lett. 589, 880–884 (2015).CAS 

Google Scholar 
Jarrett, J. T. Surprise! A hidden B12 cofactor catalyzes a radical methylation. J. Biol. Chem. 294, 11726–11727 (2019).CAS 

Google Scholar 
Byers, D. M. & Gong, H. Acyl carrier protein: structure–function relationships in a conserved multifunctional protein family. Biochem. Cell Biol. 85, 649–662 (2007).CAS 

Google Scholar 
D’Andrea, L. D. & Regan, L. TPR proteins: the versatile helix. Trends Biochem. Sci. 28, 655–662 (2003).
Google Scholar 
Ganesh, S. et al. Size-fraction partitioning of community gene transcription and nitrogen metabolism in a marine oxygen minimum zone. ISME J. 9, 2682–2696 (2015).CAS 

Google Scholar 
Ganesh, S., Parris, D. J., DeLong, E. F. & Stewart, F. J. Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. ISME J. 8, 187–211 (2014).CAS 

Google Scholar 
Fuchsman, C. A., Kirkpatrick, J. B., Brazelton, W. J., Murray, J. W. & Staley, J. T. Metabolic strategies of free-living and aggregate-associated bacterial communities inferred from biologic and chemical profiles in the Black Sea suboxic zone. FEMS Microbiol. Ecol. 78, 586–603 (2011).CAS 

Google Scholar 
Alldredge, A. L. & Cohen, Y. Can microscale chemical patches persist in the sea? Microelectrode study of marine snow, fecal pellets. Science 235, 689–691 (1987).ADS 
CAS 

Google Scholar 
Scranton, M. I. et al. Temporal variability in the nutrient chemistry of the Cariaco Basin. in Past and Present Water Column Anoxia. Nato Science Series: IV: Earth and Environmental Sciences, Vol. 64. (ed. Neretin, L.) 139–160 (Springer Dordrecht, 2006).Firn, R. D. & Jones, C. G. The evolution of secondary metabolism–a unifying model. Mol. Microbiol. 37, 989–994 (2000).CAS 

Google Scholar 
Junkins, E. N., McWhirter, J. B., McCall, L.-I. & Stevenson, B. S. Environmental structure impacts microbial composition and secondary metabolism. ISME Commun. 2, 1–10 (2022).
Google Scholar 
Penn, K. et al. Genomic islands link secondary metabolism to functional adaptation in marine Actinobacteria. ISME J. 3, 1193–1203 (2009).CAS 

Google Scholar 
Thaker, M. N. et al. Identifying producers of antibacterial compounds by screening for antibiotic resistance. Nat. Biotechnol. 31, 922–927 (2013).CAS 

Google Scholar 
Taylor, C. D. & Doherty, K. W. Submersible Incubation Device (SID), autonomous instrumentation for the in situ measurement of primary production and other microbial rate processes. Deep-Sea Res. Part A. Oceanogr. Res. Pap. 37, 343–358 (1990).ADS 
CAS 

Google Scholar 
Pachiadaki, M. G., Rédou, V., Beaudoin, D. J., Burgaud, G. & Edgcomb, V. P. Fungal and prokaryotic activities in the marine subsurface biosphere at Peru Margin and Canterbury Basin inferred from RNA-based analyses and microscopy. Front. Microbiol. 7, 846 (2016).
Google Scholar 
Frias-Lopez, J. et al. Microbial community gene expression in ocean surface waters. Proc. Natl Acad. Sci. USA 105, 3805–3810 (2008).ADS 
CAS 

Google Scholar 
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).CAS 

Google Scholar 
Stewart, F. J., Ulloa, O. & DeLong, E. F. Microbial metatranscriptomics in a permanent marine oxygen minimum zone. Environ. Microbiol. 14, 23–40 (2012).CAS 

Google Scholar 
Conroy, J. L. et al. Unprecedented recent warming of surface temperatures in the eastern tropical Pacific Ocean. Nat. Geosci. 2, 46–50 (2009).ADS 
CAS 

Google Scholar 
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).MathSciNet 
CAS 

Google Scholar 
Kang, D. D., Froula, J., Egan, R. & Wang, Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ 3, e1165 (2015).
Google Scholar 
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015).CAS 

Google Scholar 
Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 36, 1925–1927 (2020).Giovannoni, S. J., Britschgi, T. B., Moyer, C. L. & Field, K. G. Genetic diversity in Sargasso Sea bacterioplankton. Nature 345, 60–63 (1990).ADS 
CAS 

Google Scholar 
Eren, A. M. et al. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ 3, e1319 (2015).
Google Scholar 
Delmont, T. O. et al. Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes. Nat. Microbiol. 3, 804–813 (2018).CAS 

Google Scholar 
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).CAS 

Google Scholar 
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Google Scholar 
Team, R. C. R: a language and environment for statistical computing (R Foundation for Statistical Computing, 2013).Marçais, G. et al. MUMmer4: a fast and versatile genome alignment system. PLoS Comput. Biol. 14, e1005944 (2018).
Google Scholar 
Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv, 1303.3997v2 (2013).Ben Woodcroft. CoverM. https://github.com/wwood/CoverM (2022).Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 11, 1–11 (2010).
Google Scholar 
Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240 (2014).CAS 

Google Scholar 
Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).CAS 

Google Scholar 
Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100 (2018).CAS 

Google Scholar 
Bray, N. L., Pimentel, H., Melsted, P. & Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34, 525–527 (2016).CAS 

Google Scholar 
Konopka, T. umap. Uniform manifold approximation and projection (2018).Wickham, H., Chang, W. & Wickham, M. H. Package ‘ggplot2’. Create elegant data visualisations using the grammar of graphics. Version 2, 1–189 (2016).
Google Scholar 
Hahsler, M., Piekenbrock, M. & Doran, D. dbscan: fast density-based clustering with R. J. Stat. Softw. 91, 1–30 (2019).
Google Scholar 
Geller-McGrath, D. et al. Diverse secondary metabolites are expressed in particle-associated and free-living microorganisms of the permanently anoxic Cariaco Basin. https://github.com/d-mcgrath/cariaco_basin (2023). More

Oil palm plantations can supplant once biodiverse tropical forests. As planted areas expand, it is vital to plan landscapes to better balance biodiversity and oil palm production. Strategic ‘set-asides’ offer a key approach.In recent decades, oil palm has expanded spectacularly in some of the most biodiverse areas of the tropics, especially in Indonesia and Malaysia. This expansion has caused extensive deforestation (including loss of more than 2.1 million ha of primary forests in Borneo2, as well as other forests and agroforests), and management of plantations often relies heavily on clearing, herbicides and pesticides. This has generated many direct and indirect impacts on wildlife, ecosystems, climate and human communities3. Further expansion is ongoing, and global demand continues to rise4. More

Phalan, B. et al. Crop expansion and conservation priorities in tropical countries. PLoS ONE 8, e51759 (2013).Article 
CAS 

Google Scholar 
Poore, J. & Nemecek, T. Reducing food’s environmental impacts through producers and consumers. Science 360, 987–992 (2018).Article 
CAS 

Google Scholar 
Springmann, M. et al. Options for keeping the food system within environmental limits. Nature 562, 519–525 (2018).Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl Acad. Sci. USA 108, 20260–20264 (2011).Article 
CAS 

Google Scholar 
Searchinger, T. D., Wirsenius, S., Beringer, T. & Dumas, P. Assessing the efficiency of changes in land use for mitigating climate change. Nature 564, 249–253 (2018).Edwards, D. P. et al. Conservation of tropical forests in the Anthropocene. Curr. Biol. 29, R1008–R1020 (2019).Article 
CAS 

Google Scholar 
Newbold, T. et al. Global patterns of terrestrial assemblage turnover within and among land uses. Ecography 39, 1151–1163 (2016).Article 

Google Scholar 
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).Article 
CAS 

Google Scholar 
Gibbs, H. K. et al. Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proc. Natl Acad. Sci. USA 107, 16732–16737 (2010).Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).Article 
CAS 

Google Scholar 
Newbold, T. et al. A global model of the response of tropical and sub-tropical forest biodiversity to anthropogenic pressures. Proc. R. Soc. B 281, 20141371 (2014).Article 

Google Scholar 
Clough, Y. et al. Combining high biodiversity with high yields in tropical agroforests. Proc. Natl Acad. Sci. USA 108, 8311–8316 (2011).Article 
CAS 

Google Scholar 
Giam, X. Global biodiversity loss from tropical deforestation. Proc. Natl Acad. Sci. USA 114, 5775–5777 (2017).Article 
CAS 

Google Scholar 
van der Werf, G. R. et al. CO2 emissions from forest loss. Nat. Geosci. 2, 737–738 (2009).Article 

Google Scholar 
Harvey, C. A. et al. Climate‐smart landscapes: opportunities and challenges for integrating adaptation and mitigation in tropical agriculture. Conserv. Lett. 7, 77–90 (2014).Article 

Google Scholar 
Harris, N. L. et al. Baseline map of carbon emissions from deforestation in tropical regions. Science 336, 1573–1576 (2012).Article 
CAS 

Google Scholar 
Song, X.-P. et al. Global land change from 1982 to 2016. Nature 560, 639–643 (2018).Quezada, J. C., Etter, A., Ghazoul, J., Buttler, A. & Guillaume, T. Carbon neutral expansion of oil palm plantations in the Neotropics. Sci. Adv. 5, eaaw4418 (2019).Article 
CAS 

Google Scholar 
Oil Palm and Biodiversity: a Situation Analysis by the IUCN Oil Palm Task Force (International Union for Conservation of Nature, 2018). https://doi.org/10.2305/IUCN.CH.2018.11.enMeijaard, E. & Sheil, D. The moral minefield of ethical oil palm and sustainable development. Front. For. Glob. Change 2, 22 (2019).The Future of Food and Agriculture – Alternative Pathways to 2050 (FAO, 2018).Henders, S., Persson, U. M. & Kastner, T. Trading forests: land-use change and carbon emissions embodied in production and exports of forest-risk commodities. Environ. Res. Lett. 10, 125012 (2015).Article 

Google Scholar 
Donofrio, S., Rothrock, P. & Leonard, J. Supply Change: Tracking Corporate Commitments to Deforestation-Free Supply Chains (Forest Trends, 2017).Terrenoire, E., Hauglustaine, D. A., Gasser, T. & Penanhoat, O. The contribution of carbon dioxide emissions from the aviation sector to future climate change. Environ. Res. Lett. 14, 084019 (2019).Article 
CAS 

Google Scholar 
Parsons, S., Raikova, S. & Chuck, C. J. The viability and desirability of replacing palm oil. Nat. Sustain. 3, 412–418 (2020).Article 

Google Scholar 
Taheripour, F., Hertel, T. W. & Ramankutty, N. Market-mediated responses confound policies to limit deforestation from oil palm expansion in Malaysia and Indonesia. Proc. Natl Acad. Sci. USA 116, 19193–19199 (2019).Article 
CAS 

Google Scholar 
Laurance, W. F. et al. Improving the performance of the roundtable on sustainable palm oil for nature conservation. Conserv. Biol. 24, 377–381 (2010).Article 

Google Scholar 
Meijaard, E., Abrams, J. F., Juffe-Bignoli, D., Voigt, M. & Sheil, D. Coconut oil, conservation and the conscientious consumer. Curr. Biol. 30, R757–R758 (2020).Article 
CAS 

Google Scholar 
Driving Change With Sustainable Palm Oil (Roundtable on Sustainable Palm Oil, accessed August 2022). https://rspo.org/aboutGarrett, R. D., Carlson, K. M., Rueda, X. & Noojipady, P. Assessing the potential additionality of certification by the round table on responsible soybeans and the roundtable on sustainable palm oil. Environ. Res. Lett. 11, 045003 (2016).Article 

Google Scholar 
Mittermeier, R. A., Myers, N., Mittermeier, C. G. & Robles, G. Hotspots: Earth’s Biologically Richest and Most Endangered Terrestrial Ecoregions (Conservation International, 1999).Gaveau, D. L. et al. Rapid conversions and avoided deforestation: examining four decades of industrial plantation expansion in Borneo. Sci. Rep. 6, 32017 (2016).Article 
CAS 

Google Scholar 
Luke, S. H. et al. Riparian buffers in tropical agriculture: scientific support, effectiveness and directions for policy. J. Appl. Ecol. 56, 85–92 (2019).Article 

Google Scholar 
Mitchell, S. L. et al. Riparian reserves help protect forest bird communities in oil palm dominated landscapes. J. Appl. Ecol. 55, 2744–2755 (2018).Article 

Google Scholar 
Scriven, S. A. et al. Testing the benefits of conservation set-asides for improved habitat connectivity in tropical agricultural landscapes. J. Appl. Ecol. 56, 2274–2285 (2019).Article 

Google Scholar 
Deere, N. J. et al. Riparian buffers can help mitigate biodiversity declines in oil palm agriculture. Front. Ecol. Environ. 20, 459–466 (2021).Woodham, C. R. et al. Effects of replanting and retention of mature oil palm riparian buffers on ecosystem functioning in oil palm plantations. Front. Glob. Change 2, 29 (2019).Article 

Google Scholar 
Carlson, K. M. et al. Influence of watershed‐climate interactions on stream temperature, sediment yield, and metabolism along a land use intensity gradient in Indonesian Borneo. J. Geophys. Res. Biogeosci. 119, 1110–1128 (2014).Article 

Google Scholar 
Carlson, K. M. et al. Effect of oil palm sustainability certification on deforestation and fire in Indonesia. Proc. Natl Acad. Sci. USA 115, 121–126 (2018).Article 
CAS 

Google Scholar 
Fleiss, S. et al. Conservation set-asides improve carbon storage and support associated plant diversity in certified sustainable oil palm plantations. Biol. Conserv. 248, 108631 (2020).Article 

Google Scholar 
Wunder, S., Angelsen, A. & Belcher, B. Forests, livelihoods, and conservation: broadening the empirical base. World Dev. 64, S1–S11 (2014).Struebig, M. J. et al. Quantifying the biodiversity value of repeatedly logged rainforests: gradient and comparative approaches from Borneo. Adv. Ecol. Res. 48, 183–224 (2013).Article 

Google Scholar 
Shevade, V. S. & Loboda, T. V. Oil palm plantations in Peninsular Malaysia: determinants and constraints on expansion. PLoS ONE 14, e0210628 (2019).Article 
CAS 

Google Scholar 
Pirker, J., Mosnier, A., Kraxner, F., Havlík, P. & Obersteiner, M. What are the limits to oil palm expansion? Glob. Environ. Change 40, 73–81 (2016).Article 

Google Scholar 
Launching the RSPO Jurisdictional Approach (JA) Piloting Framework (Roundtable on Sustainable Palm Oil, accessed August 2022).Abram, N. K. et al. Synergies for improving oil palm production and forest conservation in floodplain landscapes. PLoS ONE 9, e95388 (2014).Article 

Google Scholar 
Othman, N. et al. Shift of paradigm needed towards improving human–elephant coexistence in monoculture landscapes in Sabah. Int. Zoo Yearb. 53, 161–173 (2019).Article 

Google Scholar 
Horton, A. J. et al. Can riparian forest buffers increase yields from oil palm plantations? Earths Future 6, 1082–1096 (2018).Article 

Google Scholar 
Ewers, R. M. et al. A large-scale forest fragmentation experiment: the Stability of Altered Forest Ecosystems Project. Phil. Trans. R. Soc. Lond. B 366, 3292–3302 (2011).Article 

Google Scholar 
Pfeifer, M. et al. Creation of forest edges has a global impact on forest vertebrates. Nature 551, 187–191 (2017).Ewers, R. M., Thorpe, S. & Didham, R. K. Synergistic interactions between edge and area effects in a heavily fragmented landscape. Ecology 88, 96–106 (2007).Article 

Google Scholar 
Deere, N. J. et al. High carbon stock forests provide co-benefits for tropical biodiversity. J. Appl. Ecol. 55, 997–1008 (2018).Article 
CAS 

Google Scholar 
Hemprich-Bennett, D. R. et al. Altered structure of bat–prey interaction networks in logged tropical forests revealed by metabarcoding. Mol. Ecol. 30, 5844–5857 (2021).Article 

Google Scholar 
Williamson, J. et al. Riparian buffers act as microclimatic refugia in oil palm landscapes. J. Appl. Ecol. 58, 431–442 (2021).Article 

Google Scholar 
Slade, E. M., Mann, D. J. & Lewis, O. T. Biodiversity and ecosystem function of tropical forest dung beetles under contrasting logging regimes. Biol. Conserv. 144, 166–174 (2011).Article 

Google Scholar 
Gray, R. E. J. et al. Movement of forest-dependent dung beetles through riparian buffers in Bornean oil palm plantations. J. Appl. Ecol. 59, 238–250 (2022).Woodman, S. M. et al. esdm: a tool for creating and exploring ensembles of predictions from species distribution and abundance models. Methods Ecol. Evol. 10, 1923–1933 (2019).Article 

Google Scholar 
Liu, C., Berry, P. M., Dawson, T. P. & Pearson, R. G. Selecting thresholds of occurrence in the prediction of species distributions. Ecography 28, 385–393 (2005).Article 

Google Scholar 
Piccini, I. et al. Greenhouse gas emissions from dung pats vary with dung beetle species and with assemblage composition. PloS ONE 12, e0178077 (2017).Article 

Google Scholar 
Raine, E. H. & Slade, E. M. Dung beetle–mammal associations: methods, research trends and future directions. Proc. R. Soc. B 286, 20182002 (2019).Article 

Google Scholar 
Nichols, E., Gardner, T., Peres, C., Spector, S. & Network, S. R. Co‐declining mammals and dung beetles: an impending ecological cascade. Oikos 118, 481–487 (2009).Article 

Google Scholar 
Asner, G. P. et al. Mapped aboveground carbon stocks to advance forest conservation and recovery in Malaysian Borneo. Biol. Conserv. 217, 289–310 (2018).Article 

Google Scholar 
Jucker, T. et al. Estimating aboveground carbon density and its uncertainty in Borneo’s structurally complex tropical forests using airborne laser scanning. Biogeosciences 15, 3811–3830 (2018).Article 

Google Scholar 
Philipson, C. D. et al. Active restoration accelerates the carbon recovery of human-modified tropical forests. Science 369, 838–841 (2020).Article 
CAS 

Google Scholar 
Nunes, M. H. et al. Recovery of logged forest fragments in a human-modified tropical landscape during the 2015–16 El Niño. Nat. Commun. 12, 1526 (2021).Article 
CAS 

Google Scholar 
Woittiez, L. S., van Wijk, M. T., Slingerland, M., van Noordwijk, M. & Giller, K. E. Yield gaps in oil palm: a quantitative review of contributing factors. Eur. J. Agron. 83, 57–77 (2017).Article 

Google Scholar  More

Plant growth model without environmental forcingThe model without environmental forcing closely follows the original description of the Thornley transport resistance (TTR) model29. A summary of the model parameters is provided in Supplementary Table 2. The shoot and root mass pools (MS and MR, in kg structural dry matter) change as a function of growth and loss (equations (1) and (2)). The litter (kL) and maintenance respiration (r) loss rates (in kg kg−1 d−1) are treated as constants. In the original model description29 r = 0. The parameter KM (units kg) describes how loss varies with mass (MS or MR). Growth (Gs and Gr, in kg d−1) varies as a function of the carbon and nitrogen concentrations (equations (3) and (4)). CS, CR, NS and NR are the amounts (kg) of carbon and nitrogen in the roots and shoots. These assumptions yield the following equations for shoot and root dry matter,$${mathrm{MS}}[t+1]={mathrm{MS}}[t]+{G}_{{mathrm{S}}}[t]-frac{({k}_{{mathrm{L}}}+r){mathrm{MS}}[t]}{1+frac{{K}_{{{M}}}}{{mathrm{MS}}[t]}},$$
(1)
$${mathrm{MR}}[t+1]={mathrm{MR}}[t]+{G}_{{mathrm{R}}}[t]-frac{({k}_{{mathrm{L}}}+r){mathrm{MR}}[t]}{1+frac{{K}_{{{M}}}}{{mathrm{MR}}[t]}},$$
(2)
where GS and GR are$${G}_{{mathrm{S}}}=gfrac{{mathrm{CS}}times {mathrm{NS}}}{{mathrm{MS}}},$$
(3)
$${G}_{{mathrm{R}}}=gfrac{{mathrm{CR}}times {mathrm{NR}}}{{mathrm{MR}}},$$
(4)
and g is the growth coefficient (in kg kg−1 d−1).Carbon uptake UC is determined by the net photosynthetic rate (a, in kg kg−1 d−1) and the shoot mass (equation (5)). Similarly, nitrogen uptake (UN) is determined by the nitrogen uptake rate (b, in kg kg−1 d−1) and the root mass. The parameter KA (units kg) forces both photosynthesis and nitrogen uptake to be asymptotic with mass. The second terms in the denominators of equations (5) and (6) model product inhibitions of carbon and nitrogen uptake, respectively; that is, the parameters JC and JN (in kg kg−1) mimic the inhibition of source activity when substrate concentrations are high,$${U}_{{mathrm{C}}}=frac{a{mathrm{MS}}}{left(1+frac{{mathrm{MS}}}{{K}_{{mathrm{A}}}}right)left(1+frac{{mathrm{CS}}}{{mathrm{MS}}times {J}_{{mathrm{C}}}}right)},$$
(5)
$${U}_{{mathrm{N}}}=frac{b{mathrm{MR}}}{left(1+frac{{mathrm{MR}}}{{K}_{{mathrm{A}}}}right)left(1+frac{{mathrm{NR}}}{{mathrm{MR}}times {J}_{{mathrm{N}}}}right)}.$$
(6)
The substrate transport fluxes of C and N (τC and τN, in kg d−1) between roots and shoots are determined by the concentration gradients between root and shoot and by the resistances. In the original model description29, these resistances are defined flexibly, but we simplify and assume that they scale linearly with plant mass,$${tau }_{{mathrm{C}}}=frac{{mathrm{MS}}times {mathrm{MR}}}{{mathrm{MS}}+{mathrm{MR}}}left(frac{{mathrm{CS}}}{{mathrm{MS}}}-frac{{mathrm{CR}}}{{mathrm{MR}}}right)$$
(7)
$${tau }_{{mathrm{N}}}=frac{{mathrm{MS}}times {mathrm{MR}}}{{mathrm{MS}}+{mathrm{MR}}}left(frac{{mathrm{NR}}}{{mathrm{MR}}}-frac{{mathrm{NS}}}{{mathrm{MS}}}right)$$
(8)
The changes in mass of carbon and nitrogen in the roots and shoots are then$${mathrm{CS}}[t+1]={mathrm{CS}}[t]+{U}_{{mathrm{C}}}[t]-{f}_{{mathrm{C}}}{G}_{{mathrm{s}}}[t]-{tau }_{{mathrm{C}}}[t]$$
(9)
$${mathrm{CR}}[t+1]={mathrm{CR}}[t]+{tau }_{{mathrm{C}}}[t]-{f}_{{mathrm{C}}}{G}_{{mathrm{r}}}[t]$$
(10)
$${mathrm{NS}}[t+1]={mathrm{NS}}[t]+{tau }_{{mathrm{N}}}[t]-{f}_{{mathrm{N}}}{G}_{{mathrm{s}}}[t]$$
(11)
$${mathrm{NR}}[t+1]={mathrm{NR}}[t]+{U}_{{mathrm{N}}}[t]-{f}_{{mathrm{N}}}{G}_{{mathrm{r}}}[t]-{tau }_{{mathrm{N}}}[t]$$
(12)
where fC and fN (in kg kg−1) are the fractions of structural carbon and nitrogen in dry matter.Adding environmental forcing to the plant growth modelIn this section, we describe how the net photosynthetic rate (a), the nitrogen uptake rate (b), the growth rate (g) and the respiration rate (r) are influenced by environmental-forcing factors. These environmental-forcing effects are described in equations (13)–(17) and summarized graphically in Extended Data Fig. 1. All other model parameters are treated as constants. Previous work that implemented the TTR model as a species distribution model30 is used as a starting point for adding environmental forcing. As in this previous work30, we assume that parameters a, b and g are co-limited by environmental factors in a manner analogous to Liebig’s law of the minimum, which is a crude but pragmatic abstraction. The implementation here differs in some details.Unlike previous work30, we use the Farquhar model of photosynthesis47,48 to represent how solar radiation, atmospheric CO2 concentration and air temperature co-limit photosynthesis35. We assume that the Farquhar model parameters are universal and that all vegetation in our study uses the C3 photosynthetic pathway. The Farquhar model photosynthetic rates are rescaled to [0,amax] to yield afqr. The effects of soil moisture (Msoil) on photosynthesis are represented as an increasing step function ({{{{S}}}}(M_{mathrm{soil}},{beta }_{1},{beta }_{2})=max left{min left(frac{M_{mathrm{soil}}-{beta }_{1}}{{beta }_{2}-{beta }_{1}},1right),0right}). This allows us to redefine a as,$$a={a}_{{mathrm{fqr}}} {{{{S}}}}(M_{mathrm{soil}},{beta }_{1},{beta }_{2})$$
(13)
The processes influencing nitrogen availability are complex, and global data products on plant available nitrogen are uncertain. We therefore assume that nitrogen uptake will vary with soil temperature and soil moisture. That is, the nitrogen uptake rate b is assumed to have a maximum rate (bmax) that is co-limited by soil temperature Tsoil and soil moisture Msoil,$$b={b}_{{mathrm{max}}} {{{{S}}}}({T}_{soil},{beta }_{3},{beta }_{4}) {{{{Z}}}}(M_{mathrm{soil}},{beta }_{5},{beta }_{6},{beta }_{7},{beta }_{8}).$$
(14)
In equation (14), we have assumed that the nitrogen uptake rate is a simple increasing and saturating function of temperature. We have also assumed that the nitrogen uptake rate is a trapezoidal function of soil moisture with low uptake rates in dry soils, higher uptake rates at intermediate moisture levels and lower rates once soils are so moist as to be waterlogged. The trapezoidal function is ({{{{Z}}}}(M_{mathrm{soil}},{beta }_{5},{beta }_{6},{beta }_{7},{beta }_{8})=max left{min left(frac{M_{mathrm{soil}}-{{{{{beta }}}}}_{5}}{{{{{{beta }}}}}_{6}-{{{{{beta }}}}}_{5}},1,frac{{{{{{beta }}}}}_{8}-M_{mathrm{soil}}}{{beta }_{8}-{beta }_{7}}right),0right}).The previous sections describe how the assimilation of carbon and nitrogen by a plant are influenced by environmental factors. The TTR model describes how these assimilate concentrations influence growth (equations (3) and (4)). In our implementation, we additionally allow the growth rate to be co-limited by temperature (soil temperature, Tsoil) and soil moisture (Msoil),$$g={g}_{{mathrm{max}}} {{{{Z}}}}({T}_{{mathrm{soil}}},{beta }_{9},{beta }_{10},{beta }_{11},{beta }_{12}) {{{{S}}}}(M_{mathrm{soil}},{beta }_{13},{beta }_{14}).$$
(15)
We use Tsoil since we assume that growth is more closely linked to soil temperature, which varies slower than air temperature. The respiration rate (r, equations (1) and (2)) increases as a function of air temperature (Tair) to a maximum rmax,$$r={r}_{{mathrm{max}}}{{{{S}}}}({T}_{{mathrm{air}}},{beta }_{15},{beta }_{16}).$$
(16)
The parameter r is best interpreted as a maintenance respiration. Growth respiration is not explicitly considered; it is implicitly incorporated in the growth rate parameter (g, equation (15)), and any temperature dependence in growth respiration is therefore assumed to be accommodated by equation (15).Fire can reduce the structural shoot mass MS as follows,$${mathrm{MS}}[t+1]={mathrm{MS}}[t](1-{{{{S}}}}(F,{beta }_{17},{beta }_{18})).$$
(17)
where F is an indicator of fire severity at a point in time (for example, burnt area) and the function S(F, β17, β18) allows MS to decrease when the fire severity indicator F is high. If F = 0, this process plays no role. This fire impact equation was used in preliminary analyses, but the data on fire activity did not provide sufficient information to estimate β17 and β18; we therefore excluded this process from the final analyses.We further estimate two additional β parameters (βa and βb) so that each site can have unique maximum carbon and nitrogen uptake rates. Specifically, we redefine a as ({a}^{{prime} }={beta }_{a} a) and b as ({b}^{{prime} }={beta }_{b} b).Data sources and preparationTo describe vegetation activity, we use the GIMMS 3g v.1 NDVI data26,27 and the MODIS EVI28 data. The GIMMS data product has been derived from the AVHRR satellite programme and controls for atmospheric contamination, calibration loss, orbital drift and volcanic eruptions26,27. The data provide 24 NDVI raster grids for each year, starting in July 1981 and ending in December 2015. The spatial resolution is 1/12° (~9 × 9 km). The EVI data used are from the MODIS programme’s Terra satellite; it is a 1 km data product provided at a 16-day interval. We use data from the start of the record (February 2000) to December 2019. The MODIS data product (MOD13A2) uses a temporal compositing algorithm to produce an estimate every 16 days that filters out atmospheric contamination. The EVI is designed to reduce the effects of atmospheric, bare-ground and surface water on the vegetation index28.For environmental forcing, we use the ERA5-Land data31,32 (European Centre for Medium-Range Weather Forecasts Reanalysis v. 5; hereafter, ERA5). The ERA5 products are global reanalysis products based on hourly estimates of atmospheric variables and extend from present back to 1979. The data products are supplied at a variety of spatial and temporal resolutions. We used the monthly averages from 1981 to 2019 at a 0.1° spatial resolution (~11 km). The ERA5 data provide air temperature (2 m surface air temperature), soil temperature (0–7 cm soil depth), surface solar radiation and volumetric soil water (0–7 cm soil depth). Fire data were taken from the European Space Agency Fire Disturbance Climate Change Initiative’s AVHRR Long-Term Data Record Grid v.1.0 product49. This product provides gridded (0.25° resolution) data of monthly global (from 1982 to 2017) burned area derived from the AVHRR satellite programme. As mentioned, the fire data did not enrich our analysis, and the analyses we present here therefore exclude further consideration of the fire data.All data were resampled to the GIMMS grid. The mean pixel EVI was then calculated for each GIMMS cell for each time point in the MODIS EVI data. We used linear interpolation on the NDVI, EVI and ERA5 environmental-forcing data to estimate each variable on a weekly time step. This served to set the time step of the TTR difference equations to one week and to synchronize the different time series.Site selectionThe GIMMS grid cells define the spatial resolution of our sample points. GIMMS grid cells are large (1/12°, ~9 km), meaning that most grid cells contain multiple land-cover types. We focused on wilderness landscapes, and our aim was to find multiple grid cells for the major ecosystems of the world. We used the following classification of ecosystem types to guide the stratification of our grid-cell selection: tropical evergreen forest (RF), boreal forest (BF), temperate evergreen and temperate deciduous forest (TF), savannah (SA), shrubland (SH), grassland (GR), tundra (TU) and Mediterranean-type ecosystems (MT).We used the following criteria to select grid cells. (1) Selected grid cells should contain relatively homogeneous vegetation. Small-scale heterogeneity was allowed (for example, catenas, drainage lines, peatlands) as long as many of these elements are repeated in the scene (for example, rolling hills were accepted, but elevation gradients were rejected). (2) Sites should have no signs of transformative human activity (for example, tree harvesting, crop cultivation, paved surfaces). We used the Time Tool in Google Earth Pro, which provides annual satellite imagery of the Earth from 1984 onwards, to ensure that no such activity occurred during the observation period (note that the GIMMS record starts in July 1981; however, it is likely that evidence of transformative activity between July 1981 and 1984 would be visible in 1984). Grid cells with extensive livestock holding on non-improved pasture were included. In some cases, a small agricultural field or pasture was present, and such grid cells were used as long as the field or pasture was small and remained constant in size. (3) Grid cells should not include large water bodies, but small drainage lines or ponds were accepted as long as they did not violate the first criterion. (4) Grid cells should be independent (neighbouring grid cells were not selected) and cover the major ecosystems of the world. Using these criteria, we were able to include 100 sites in the study (Extended Data Figs. 2 and 3 and Supplementary Table 4).State-space modelWe used a Bayesian state-space approach. Conceptually, the analysis takes the form,$$M[t]=f(M[t-1],{boldsymbol{beta}},{boldsymbol{theta}}_{t-1},{epsilon }_{t-1})$$
(18)
$${mathrm{VI}}[t]=m M[t]+eta .$$
(19)
Here M[t] is the plant growth model’s prediction of biomass (M = MS + MR) at time t, and ϵt−1 is the process error associated with the state variables. In the model, each underlying state variable (MS, MR, CS, CR, NS and NR) has an associated process error term. The function f(M[t − 1], β, θt−1, ϵt−1) summarizes that the development of M is influenced by the state variables MS, MR, CS, CR, NS and NR, the environmental-forcing data θt−1 and the β parameters. The observation equation (equation (19)) uses the parameter m to link the VI (vegetation index, either NDVI or EVI) observations to the modelled state M. The parameter η is the observation error. Equation (19) assumes that there is a linear relationship between modelled biomass (M) and VI, which is a simplification of reality50,51,52. The observation error η is structured by our simplification of the data products quality scores (coded Q = 0, 1, 2, with 0 being good and 2 being poor; Supplementary Table 3) to allow the error to increase with each level of the quality score. Specifically, we define η = e0 + e1 × Q.The model was formulated using the R package LaplacesDemon53. All β parameters are given vague uniform priors. The parameter m is given a vague normal prior (truncated to be >0). The process error terms are modelled using normal distributions, and the variances of the error terms are given vague half-Cauchy priors. The ex parameters are given vague normal priors. The model also requires the parameterization of M[0], the initial vegetation biomass; M[0] is given a vague uniform prior. We used the twalk Markov chain Monte Carlo (MCMC) algorithm as implemented in LaplacesDemon53 and its default control parameters to estimate the posterior distributions of the model parameters. We further fitted the model using DEoptim54,55, which is a robust genetic algorithm that is known to perform stably on high-dimensional and multi-modal problems56, to verify that the MCMC algorithm had not missed important regions of the parameter space. The models estimated with MCMC had significantly lower log root-mean-square error than models estimated with DEoptim (paired t-test NDVI analysis: t = –2.9806, degrees of freedom (d.f.) = 99, P = 0.00362; EVI analysis: t = –4.6229, d.f. = 99, P = 1.144 × 10–5), suggesting that the MCMC algorithm performed well compared with the genetic algorithm.Anomaly extraction and trend estimationWe use the ‘seasonal and trend decomposition using Loess’ (STL57) as implemented in the R58 base function stl. STL extracts the seasonal component s of a time series using Loess smoothing. What remains after seasonal extraction (the anomaly or remainder, r) is the sum of any long-term trend and stochastic variation. We estimate the trend in two ways. First, we estimate the trend by fitting a quadratic polynomial (r = a + bx + cx2) to the remainder (r is the remainder, x is time and a, b and c are regression coefficients). The use of polynomials allows the data to specify whether a trend exists, whether the trend is linear, cup or hat shaped and whether the overall trend is increasing or decreasing. As a second method, we estimate the trend by fitting a bent-cable regression to the remainder. Bent-cable regression is a type of piecewise linear regression for estimating the point of transition between two linear phases in a time series59,60. The model takes the form r = b0 + b1x + b2 q(x, τ, γ)60. Here r is the remainder, x is time, b0 is the initial intercept, b1 is the slope in phase 1, the slope in phase 2 is b2 − b1 and q is a function that defines the change point: (q(x,tau ,gamma )=frac{{(x-tau +gamma )}^{2}}{4gamma }I(tau -gamma < tau +gamma )+(x-tau )I(x > tau +gamma )); τ represents the location of the change point and γ the span of the bent cable that joins the two linear phases; I(A) is an indicator function that returns 1 if A is true and 0 if A is false. The bent-cable model allows the data to specify whether a trend exists and whether there has been a switch in the trend, thereby allowing the identification of whether the trend is linear, cup or hat shaped and whether the overall trend is increasing or decreasing. Both the polynomial and bent-cable models were estimated using LaplacesDemon’s53 Adaptive Metropolis MCMC algorithm and vague priors, although for the bent-cable model we constrained τ to be in the middle 70% of the time series and γ to be at most two years.The STL extraction of the seasonal components in the air temperature, soil temperature, soil moisture and solar radiation data (there is no stochasticity or seasonal trend in the CO2 data we used) allows us to simulate detrended time series d of these forcing variables as (d=bar{y}+s+{{{{N}}}}(mu ,sigma )) where N(μ, σ) is a normally distributed random variable with mean and standard deviation estimated from the remainder r (we verified that r was well described by the normal distribution), (bar{y}) is the mean of the data over the time series and s is the seasonal component extracted by STL. For CO2, the detrended time series is simply the average CO2 over the time series. More