Philippa Kaur

1.Bonan, G. B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).CAS 
PubMed 
Article 
ADS 

Google Scholar 
2.Pugh, T. A. M. et al. Role of forest regrowth in global carbon sink dynamics. Proc. Natl. Acad. Sci. U. S. A. 116, 4382–4387 (2019).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
3.Allen, C. D. et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 259, 660–684 (2010).Article 

Google Scholar 
4.Williams, C. A., Gu, H., MacLean, R., Masek, J. G. & Collatz, G. J. Disturbance and the carbon balance of US forests: A quantitative review of impacts from harvests, fires, insects, and droughts. Glob. Planet. Change 143, 66–80 (2016).Article 
ADS 

Google Scholar 
5.Kurz, W. A. et al. Mountain pine beetle and forest carbon feedback to climate change. Nature 452, 987–990 (2008).CAS 
PubMed 
Article 
ADS 

Google Scholar 
6.Lovett, G. M. et al. Nonnative forest insects and pathogens in the United States: Impacts and policy options. Ecol. Appl. 26, 1437–1455 (2016).PubMed 
PubMed Central 
Article 

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

Google Scholar 
8.Nave, L. E. et al. Reforestation can sequester two petagrams of carbon in US topsoils in a century. Proc. Natl. Acad. Sci. U. S. A. 115, 2776–2781 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
9.Millar, C. I., Stephenson, N. L. & Stephens, S. L. Climate change and forests of the future: Managing in the face of uncertainty. Ecol. Appl. 17, 2145–2151 (2007).PubMed 
Article 

Google Scholar 
10.McCarthy, J. K., Dwyer, J. M. & Mokany, K. A regional-scale assessment of using metabolic scaling theory to predict ecosystem properties. Proc. Biol. Sci. 286, 20192221 (2019).PubMed 
PubMed Central 

Google Scholar 
11.Woodall, C. W., Miles, P. D. & Vissage, J. S. Determining maximum stand density index in mixed species stands for strategic-scale stocking assessments. For. Ecol. Manag. 216, 367–377 (2005).Article 

Google Scholar 
12.Reineke, L. H. Perfecting a stand-density index for even-aged forests. J. Agric. Res. 46, 627–638 (1933).
Google Scholar 
13.Long, J. N. A practical approach to density management. For. Chron. 61, 23–27 (1985).Article 

Google Scholar 
14.Domke, G. et al. Forests. In Second State of the Carbon Cycle Report (SOCCR2): A Sustained Assessment Report (eds Cavallaro, N., Shrestha, G., Birdsey, R., Mayes, M. A., Najjar, R. G., Reed, S. C., Romero-Lankao, P. & Zhu, Z.) 365–398 (US Global Change Research Program, 2018).15.Yoda, K., Kira, T., Ogawa, H. & Hozumi, K. Self-thinning in overcrowded pure stands under cultivated and natural conditions. J. Biol. Osaka City Univ. 14, 106–129 (1963).
Google Scholar 
16.Drew, T. J. & Flewelling, J. W. Stand density management: An alternative approach and its application to Douglas-fir plantations. For. Sci. 25, 518–532 (1979).
Google Scholar 
17.Bechtold, W. A. & Patterson, P. L. The Enhanced Forest Inventory and Analysis Program: National Sampling Design and Estimation Procedures. SRS GTR-80. USDA Forest Service, Southern Research Station, Asheville, North Carolina, USA. (2005). https://doi.org/10.2737/SRS-GTR-80.18.McGill, B. J., Enquist, B. J., Weiher, E. & Westoby, M. Rebuilding community ecology from functional traits. Trends Ecol. Evol. 21, 178–185 (2006).PubMed 
Article 

Google Scholar 
19.Andrews, C., Weiskittel, A., D’Amato, A. W. & Simons-Legaard, E. Variation in the maximum stand density index and its linkage to climate in mixed species forests of the North American Acadian Region. For. Ecol. Manag. 417, 90–102 (2018).Article 

Google Scholar 
20.Nagel, L. M. et al. Adaptive silviculture for climate change: A national experiment in manager–scientist partnerships to apply an adaptation framework. J. For. 115, 167–178 (2017).
Google Scholar 
21.Pretzsch, H. & Biber, P. A re-evaluation of the Reineke’s rule and stand density index. For. Sci. 51, 304–320 (2005).
Google Scholar 
22.Condés, S. et al. Climate influences on the maximum size-density relationship in Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) stands. For. Ecol. Manag. 385, 295–307 (2017).Article 

Google Scholar 
23.Ducey, M. J., Woodall, C. W. & Bravo-Oviedo, A. Climate and species functional traits influence maximum live tree stocking in the Lake States, USA. For. Ecol. Manag. 386, 51–61 (2017).Article 

Google Scholar 
24.Zhao, D., Bullock, B. P., Montes, C. R. & Wang, M. Rethinking maximum stand basal area and maximum SDI from the aspect of stand dynamics. For. Ecol. Manag. 475, 118462 (2020).Article 

Google Scholar 
25.Weiskittel, A. R. & Kuehne, C. Evaluating and modeling variation in site-level maximum carrying capacity of mixed-species forest stands in the Acadian Region of northeastern North America. For. Chron. 95, 171–182 (2019).Article 

Google Scholar 
26.Pretzsch, H. & del Río, M. Density regulation of mixed and mono-specific forest stands as a continuum: A new concept based on species-specific coefficients for density equivalence and density modification. For. Int. J. For. Res. 93, 1–15 (2020).
Google Scholar 
27.Senf, C., Buras, A., Zang, C. S., Rammig, A. & Seidl, R. Excess forest mortality is consistently linked to drought across Europe. Nat. Commun. 11, 6200 (2020).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
28.Woodall, C. W., Perry, C. H. & Miles, P. D. The relative density of forests in the United States. For. Ecol. Manag. 226, 368–372 (2006).Article 

Google Scholar 
29.Venturas, M. D., Todd, H. N., Trugman, A. T. & Anderegg, W. R. L. Understanding and predicting forest mortality in the western United States using long-term forest inventory data and modeled hydraulic damage. New Phytol. 230, 1896–1910 (2020).PubMed 
Article 

Google Scholar 
30.Higuera, P. E. & Abatzoglou, J. T. Record-setting climate enabled the extraordinary 2020 fire season in the western United States. Glob. Change Biol. 27, 1–2 (2021).Article 
ADS 

Google Scholar 
31.Peters, M. P. & Iverson, L. R. Projected drought for the conterminous United States in the 21st century. In Effects of Drought on Forests and Rangelands in the United States (eds Vose, J. M., Peterson, D. L., Luce, C. H. & Patel-Weynand, T.) vol. Gen. Tech. Rep. WO-98 19–39 (USDA Forest Service, 2019).32.Coulston, J. W., Woodall, C. W., Domke, G. M. & Walters, B. F. Refined forest land use classification with implications for United States national carbon accounting. Land Use Policy 59, 536–542 (2016).Article 

Google Scholar 
33.Wear, D. N. & Coulston, J. W. From sink to source: Regional variation in U.S. forest carbon futures. Sci. Rep. 5, 16518 (2015).PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
34.Senf, C., Sebald, J. & Seidl, R. Increasing canopy mortality affects the future demographic structure of Europe’s forests. One Earth 4, 749–755 (2021).Article 

Google Scholar 
35.Morin, X., Fahse, L., Scherer-Lorenzen, M. & Bugmann, H. Tree species richness promotes productivity in temperate forests through strong complementarity between species. Ecol. Lett. 14, 1211–1219 (2011).PubMed 
Article 

Google Scholar 
36.Griscom, B. W. et al. Natural climate solutions. Proc. Natl. Acad. Sci. U. S. A. 114, 11645–11650 (2017).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
37.Gunn, J. S., Ducey, M. J. & Belair, E. Evaluating degradation in a North American temperate forest. For. Ecol. Manag. 432, 415–426 (2019).Article 

Google Scholar 
38.Domke, G. M., Oswalt, S. N., Walters, B. F. & Morin, R. S. Tree planting has the potential to increase carbon sequestration capacity of forests in the United States. Proc. Natl. Acad. Sci. U. S. A. https://doi.org/10.1073/pnas.2010840117 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
39.King, D. I. & Schlossberg, S. Synthesis of the conservation value of the early-successional stage in forests of eastern North America. For. Ecol. Manag. 324, 186–195 (2014).Article 

Google Scholar 
40.Stephens, S. L. et al. Forest restoration and fuels reduction: Convergent or divergent?. Bioscience 71, 85–101 (2020).
Google Scholar 
41.Berner, L. T., Law, B. E., Meddens, A. J. H. & Hicke, J. A. Tree mortality from fires, bark beetles, and timber harvest during a hot and dry decade in the western United States (2003–2012). Environ. Res. Lett. 12, 065005 (2017).Article 
ADS 

Google Scholar 
42.Stanke, H., Finley, A. O., Domke, G. M., Weed, A. S. & MacFarlane, D. W. Over half of western United States’ most abundant tree species in decline. Nat. Commun. 12, 451 (2021).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
43.Weiskittel, A. R., Gould, P. J. & Temesgen, H. Sources of variation in the self-thinning boundary line for three species with varying levels of shade tolerance. For. Sci. 55, 84–93 (2009).
Google Scholar 
44.Ducey, M. J. & Knapp, R. A. A stand density index for complex mixed species forests in the northeastern United States. For. Ecol. Manag. 260, 1613–1622 (2010).Article 

Google Scholar 
45.Kurz, W. A., Stinson, G., Rampley, G. J., Dymond, C. C. & Neilson, E. T. Risk of natural disturbances makes future contribution of Canada’s forests to the global carbon cycle highly uncertain. Proc. Natl. Acad. Sci. U. S. A. 105, 1551–1555 (2008).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
46.Seidl, R., Schelhaas, M.-J. & Lexer, M. J. Unraveling the drivers of intensifying forest disturbance regimes in Europe. Glob. Change Biol. 17, 2842–2852 (2011).Article 
ADS 

Google Scholar 
47.Nelson, M. D. et al. Defining the United States land base: A technical document supporting the USDA Forest Service 2020 RPA assessment. In Gen. Tech. Rep. NRS-191, Vol. 191, 1–70 (2020).48.Patterson, P. L. & Reams, G. A. Combining panels for forest inventory and analysis estimation. Gen. Tech. Rep. SRS-80. Asheville, NC: US Department of Agriculture, Forest Service, 79–84 (2005).49.Bailey, R. G. Delineation of ecosystem regions. Environ. Manag. 7, 365–373 (1983).Article 
ADS 

Google Scholar 
50.Salas-Eljatib, C. & Weiskittel, A. R. Evaluation of modeling strategies for assessing self-thinning behavior and carrying capacity. Ecol. Evol. 8, 10768–10779 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
51.Geraci, M. Linear quantile mixed models: The lqmm package for Laplace quantile regression. J. Stat. Softw. 57(13), 1–29. http://www.jstatsoft.org/v57/i13/ (2013).
Google Scholar 
52.R Development Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).
Google Scholar 
53.Wang, T., Hamann, A., Spittlehouse, D. & Carroll, C. Locally downscaled and spatially customizable climate data for historical and future periods for North America. PLoS ONE 11, e0156720 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
54.Omernik, J. M. & Griffith, G. E. Ecoregions of the conterminous United States: Evolution of a hierarchical spatial framework. Environ. Manag. 54, 1249–1266 (2014).Article 
ADS 

Google Scholar 
55.De’ath, G. Boosted trees for ecological modeling and prediction. Ecology 88, 243–251 (2007).PubMed 
Article 

Google Scholar 
56.Long, J. N. & Daniel, T. W. Assessment of growing stock in uneven-age stands. West. J. Appl. For. 11, 59–61 (1990).Article 

Google Scholar 
57.Yang, L. et al. A new generation of the United States National Land Cover Database: Requirements, research priorities, design, and implementation strategies. ISPRS J. Photogramm. Remote Sens. 146, 108–123 (2018).Article 
ADS 

Google Scholar  More

1.Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377 (2017).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
2.Duarte, C. M. et al. Rebuilding marine life. Nature 580, 39–51 (2020).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
3.Darling, E. S. et al. Relationships between structural complexity, coral traits, and reef fish assemblages. Coral Reefs 36, 561–575 (2017).ADS 
Article 

Google Scholar 
4.McWilliam, M., Chase, T. J. & Hoogenboom, M. O. Neighbor diversity regulates the productivity of coral assemblages. Curr. Biol. 28, 3634–3639 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
5.Graham, N. A. J., Jennings, S., MacNeil, M. A., Mouillot, D. & Wilson, S. K. Predicting climate-driven regime shifts versus rebound potential in coral reefs. Nature 518, 94–97 (2015).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
6.Hughes, T. P. et al. Global warming transforms coral reef assemblages. Nature 556, 492–496 (2018).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
7.Cornwall, C. E. et al. Global declines in coral reef calcium carbonate production under ocean acidification and warming. Proc. Natl. Acad. Sci. 118, e2015265118 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
8.Gardner, T. A. Long-term region-wide declines in caribbean corals. Science 301, 958–960 (2003).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
9.De’ath, G., Fabricius, K. E., Sweatman, H. & Puotinen, M. The 27—year decline of coral cover on the Great Barrier Reef and its causes. Proc. Natl. Acad. Sci. USA https://doi.org/10.1073/pnas.1208909109 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
10.Madin, J. S. et al. Cumulative effects of cyclones and bleaching on coral cover and species richness at Lizard Island. Mar. Ecol. Prog. Ser. 604, 263–268 (2018).ADS 
Article 

Google Scholar 
11.Dietzel, A., Bode, M., Connolly, S. R. & Hughes, T. P. Long-term shifts in the colony size structure of coral populations along the Great Barrier Reef: Demographic change in Australia’s corals. Proc. R. Soc. B Biol. Sci. 287, 20201432 (2020).Article 

Google Scholar 
12.Claar, D. C. et al. Dynamic symbioses reveal pathways to coral survival through prolonged heatwaves. Nat. Commun. 11, 1–10 (2020).ADS 
Article 
CAS 

Google Scholar 
13.Claar, D. C. & Baum, J. K. Timing matters: Survey timing during extended heat stress can influence perceptions of coral susceptibility to bleaching. Coral Reefs 38, 559–565 (2019).ADS 
Article 

Google Scholar 
14.Edmunds, P. J. Vital rates of small reef corals are associated with variation in climate. Limnol. Oceanogr. 66, 901–913 (2021).ADS 
Article 

Google Scholar 
15.Hall, T. E. et al. Stony coral populations are more sensitive to changes in vital rates in disturbed environments. Ecol. Appl. 31, 1–11 (2021).Article 

Google Scholar 
16.Madin, J. S., Baird, A. H., Dornelas, M. & Connolly, S. R. Mechanical vulnerability explains size-dependent mortality of reef corals. Ecol. Lett. 17, 1008–1015 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
17.Edmunds, P. J. & Riegl, B. Urgent need for coral demography in a world where corals are disappearing. Mar. Ecol. Prog. Ser. 635, 233–242 (2020).ADS 
Article 

Google Scholar 
18.Hughes, T. P. et al. Ecological memory modifies the cumulative impact of recurrent climate extremes. Nat. Clim. Chang. 9, 40–43 (2019).ADS 
Article 

Google Scholar 
19.Pratchett, M. et al. Spatial, temporal and taxonomic variation in coral growth—Implications for the structure and function of coral reef ecosystems. Oceanogr. Mar. Biol. Ann. Rev. 53, 215–295 (2015).
Google Scholar 
20.Cantin, N. E. & Lough, J. M. Surviving coral bleaching events: Porites growth anomalies on the great barrier reef. PLoS ONE 9, e88720 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
21.Linares, C., Pratchett, M. S. & Coker, D. J. Recolonisation of Acropora hyacinthus following climate-induced coral bleaching on the Great Barrier Reef. Mar. Ecol. Prog. Ser. 438, 97–104 (2011).ADS 
Article 

Google Scholar 
22.Victor, S., Golbuu, Y., Yukihira, H. & Van Woesik, R. Acropora size-frequency distributions reflect spatially variable conditions on coral reefs of Palau. Bull. Mar. Sci. 85, 149–157 (2009).
Google Scholar 
23.Wilson, S. K., Robinson, J. P. W., Chong-Seng, K., Robinson, J. & Graham, N. A. J. Boom and bust of keystone structure on coral reefs. Coral Reefs 38, 625–635 (2019).ADS 
Article 

Google Scholar 
24.Pratchett, M. S., McWilliam, M. J. & Riegl, B. Contrasting shifts in coral assemblages with increasing disturbances. Coral Reefs 39, 783–793 (2020).Article 

Google Scholar 
25.Loya, Y. et al. Coral bleaching: The winners and the losers. Ecol. Lett. 4, 122–131 (2001).Article 

Google Scholar 
26.Van Woesik, R., Sakai, K., Ganase, A. & Loya, Y. Revisiting the winners and the losers a decade after coral bleaching. Mar. Ecol. Prog. Ser. 434, 67–76 (2011).ADS 
Article 

Google Scholar 
27.McWilliam, M., Pratchett, M. S., Hoogenboom, M. O. & Hughes, T. P. Deficits in functional trait diversity following recovery on coral reefs. Proc. R. Soc. B Biol. Sci. 287, 20192628 (2020).Article 

Google Scholar 
28.Marshall, P. A. & Baird, A. H. Bleaching of corals on the Great Barrier Reef: Differential susceptibilities among taxa. Coral Reefs 19, 155–163 (2000).Article 

Google Scholar 
29.Graham, N. A. J., Cinner, J. E., Norström, A. V. & Nyström, M. Coral reefs as novel ecosystems: Embracing new futures. Curr. Opin. Environ. Sustain. 7, 9–14 (2014).Article 

Google Scholar 
30.Sully, S., Burkepile, D. E., Donovan, M. K., Hodgson, G. & van Woesik, R. A global analysis of coral bleaching over the past two decades. Nat. Commun. 10, 1–5 (2019).CAS 
Article 

Google Scholar 
31.Gilmour, J. P., Smith, L. D., Heyward, A. J., Baird, A. H. & Pratchett, M. S. Recovery of an isolated coral reef system following severe disturbance. Science 340, 69–71 (2013).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
32.Hughes, T. P. et al. Global warming impairs stock–recruitment dynamics of corals. Nature 568, 387–390 (2019).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
33.Vercelloni, J. et al. Forecasting intensifying disturbance effects on coral reefs. Glob. Chang. Biol. 26, 2785–2797 (2020).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
34.Team, R. C. R: A Language and Environment for Statistical Computing. (2020).35.Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378 (2017).Article 

Google Scholar 
36.Evans, R. D. et al. Early recovery dynamics of turbid coral reefs after recurring bleaching events. J. Environ. Manag. 268, 110666 (2020).Article 

Google Scholar 
37.Carlot, J. et al. Juvenile corals underpin coral reef carbonate production after disturbance. Glob. Chang. Biol. 27, 2623–2632 (2021).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
38.Bellwood, D. R. et al. Coral reef conservation in the Anthropocene: Confronting spatial mismatches and prioritizing functions. Biol. Conserv. 236, 604–615 (2019).Article 

Google Scholar 
39.Baird, A., Emslie, M. & Lewis, A. Extended periods of coral recruitment on the Great Barrier Reef. In Proc. 12th Int. Coral Reef Symp. (2012).40.Foster, N. L., Baums, I. B. & Mumby, P. J. Sexual vs. asexual reproduction in an ecosystem engineer: The massive coral Montastraea annularis. J. Anim. Ecol. 76, 384–391 (2007).PubMed 
Article 
PubMed Central 

Google Scholar 
41.Edmunds, P. J. Patterns in the distribution of juvenile corals and coral reef community structure in St. John, US Virgin Islands. Mar. Ecol. Prog. Ser. 202, 113–124 (2000).ADS 
Article 

Google Scholar 
42.Hughes, T. P., Linares, C., Dakos, V., van de Leemput, I. A. & van Nes, E. H. Living dangerously on borrowed time during slow, unrecognized regime shifts. Trends Ecol. Evol. 28, 149–155 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
43.Hughes, T. P. et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science 359, 80–83 (2018).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
44.Wismer, S., Tebbett, S. B., Streit, R. P. & Bellwood, D. R. Spatial mismatch in fish and coral loss following 2016 mass coral bleaching. Sci. Total Environ. 650, 1487–1498 (2019).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
45.Wismer, S., Tebbett, S. B., Streit, R. P. & Bellwood, D. R. Young fishes persist despite coral loss on the Great Barrier Reef. Commun. Biol. 2, 1–7 (2019).Article 

Google Scholar 
46.Abràmoff, M. D., Hospitals, I., Magalhães, P. J. & Abràmoff, M. Image processing with ImageJ. Biophotonics Int. 11, 36–42 (2004).
Google Scholar  More

1.Clapham, P. J. Humpback whale: Megaptera novaeangliae. In Encyclopedia of Marine Mammals 489–492 (Elsevier, 2018).Chapter 

Google Scholar 
2.Corkeron, P. J. & Connor, R. C. Why do baleen whales migrate?. Mar. Mamm. Sci. 15, 1228–1245 (1999).Article 

Google Scholar 
3.Geijer, C. K. A., Notarbartolo di Sciara, G. & Panigada, S. Mysticete migration revisited: Are Mediterranean fin whales an anomaly?. Mamm. Rev. 46, 284–296 (2016).Article 

Google Scholar 
4.Baker, C. S. & Herman, L. M. Aggressive behavior between humpback whales (Megaptera novaeangliae) wintering in Hawaiian waters. Can. J. Zool. 62, 1922–1937 (1984).Article 

Google Scholar 
5.Herman, L. M. The multiple functions of male song within the humpback whale (Megaptera novaeangliae) mating system: Review, evaluation, and synthesis. Biol. Rev. 92, 1795–1818 (2017).PubMed 
Article 

Google Scholar 
6.Palsbøll, P. J., Clapham, P. J., Mattila, D. K. & Vasquez, O. Composition and dynamics of humpback whale competitive groups in the West Indies. Behaviour 122, 182–194 (1992).Article 

Google Scholar 
7.Payne, R. S. & Mcvay, S. Songs of Humpback Whales. Science 173, 585–597 (1971).CAS 
PubMed 
Article 
ADS 

Google Scholar 
8.Kroodsma, D. E. & Byers, B. E. The function (s) of bird song. Am. Zool. 31, 318–328 (1991).Article 

Google Scholar 
9.Garland, E. C. et al. Dynamic horizontal cultural transmission of humpback whale song at the Ocean Basin Scale. Curr. Biol. 21, 687–691 (2011).CAS 
PubMed 
Article 

Google Scholar 
10.Noad, M. J. & Cato, D. H. Swimming speeds of singing and non-singing humpback whales during migration. Mar. Mamm. Sci. 23, 481–495 (2007).Article 

Google Scholar 
11.Smith, J. N., Goldizen, A. W., Dunlop, R. A. & Noad, M. J. Songs of male humpback whales, Megaptera novaeangliae, are involved in intersexual interactions. Anim. Behav. 76, 467–477 (2008).Article 

Google Scholar 
12.Ross-Marsh, E., Elwen, S., Prinsloo, A., James, B. & Gridley, T. Singing in South Africa: Monitoring the occurrence of humpback whale (Megaptera novaeangliae) song near the Western Cape. Bioacoustics 30, 163–179 (2020).Article 

Google Scholar 
13.Stimpert, A. K., Peavey, L. E., Friedlaender, A. S. & Nowacek, D. P. Humpback whale song and foraging behavior on an Antarctic feeding ground. PLoS ONE 7, e51214 (2012).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
14.Vu, E. T. et al. Humpback whale song occurs extensively on feeding grounds in the western North Atlantic Ocean. Aquat. Biol. 14, 175–183 (2012).Article 

Google Scholar 
15.McSweeney, D., Chu, K., Dolphin, W. & Guinee, L. North Pacific humpback whale songs: A comparison of southeast Alaskan feeding ground songs with Hawaiian wintering ground songs. Mar. Mamm. Sci. 5, 139–148 (1989).Article 

Google Scholar 
16.Kowarski, K., Evers, C., Moors-Murphy, H., Martin, B. & Denes, S. L. Singing through winter nights: Seasonal and diel occurrence of humpback whale (Megaptera novaeangliae) calls in and around the Gully MPA, offshore eastern Canada. Mar. Mamm. Sci. 34, 169–189 (2018).Article 

Google Scholar 
17.Clark, C. W. & Clapham, P. J. Acoustic monitoring on a humpback whale (Megaptera novaeangliae) feeding ground shows continual singing into late spring. Proc. R. Soc. B 271, 1051–1057 (2004).PubMed 
PubMed Central 
Article 

Google Scholar 
18.International Whaling Commission. Report on the workshop on the comprehensive assessment of Southern Hemisphere humpback whales. J. Cetac. Res. Manag. 3, 1–50 (2011).
Google Scholar 
19.International Whaling Commission. Annex H: Report of the Sub-Committee on Other Southern Hemisphere Whale Stocks. (2016).20.Garland, E. C. et al. Humpback whale song on the Southern Ocean feeding grounds: Implications for cultural transmission. PLoS ONE 8, e79422 (2013).PubMed 
PubMed Central 
Article 
ADS 
CAS 

Google Scholar 
21.Gabriele, C. & Frankel, A. The occurrence and significance of humpback whale songs in Glacier Bay, Southeastern Alaska. Arctic Res. USA 16, 42–47 (2002).
Google Scholar 
22.Payne, R. & Guinee, L. Humpback whales (Megaptera novaeangliae) songs as an indicator of stocks. In Communication and Behavior of Whales (ed. Payne, R.) 333–358 (Westview Press, 1983).
Google Scholar 
23.Payne, K. & Payne, R. Large scale changes over 19 years in songs of humpback whales in Bermuda. Z. Tierpsychol. 68, 89–114 (1985).Article 

Google Scholar 
24.Winn, H. et al. Song of the humpback whale—population comparisons. Behav. Ecol. Sociobiol. 8, 41–46 (1981).Article 

Google Scholar 
25.Winn, H. & Winn, L. The song of the humpback whale Megaptera novaeangliae in the West Indies. Mar. Biol. 47, 97–114 (1978).Article 

Google Scholar 
26.Cholewiak, D. M., Sousa-Lima, R. S. & Cerchio, S. Humpback whale song hierarchical structure: Historical context and discussion of current classification issues. Mar. Mamm. Sci. 29, E312–E332 (2013).Article 

Google Scholar 
27.Kowarski, K., Moors-Murphy, H., Maxner, E. & Cerchio, S. Western North Atlantic humpback whale fall and spring acoustic repertoire: Insight into onset and cessation of singing behavior. J. Acoust. Soc. Am. 145, 2305–2316 (2019).PubMed 
Article 
ADS 

Google Scholar 
28.Magnúsdóttir, E. E. & Lim, R. Subarctic singers: Humpback whale (Megaptera novaeangliae) song structure and progression from an Icelandic feeding ground during winter. PLoS ONE 14, e0210057 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
29.Magnúsdóttir, E. E. et al. Humpback whale (Megaptera novaeangliae) song unit and phrase repertoire progression on a subarctic feeding ground. J. Acoust. Soc. Am. 138, 3362–3374 (2015).PubMed 
Article 
ADS 

Google Scholar 
30.Mattila, D. K., Guinee, L. N. & Mayo, C. A. Humpback whale songs on a North Atlantic feeding ground. J. Mammal. 68, 880–883 (1987).Article 

Google Scholar 
31.Teschke, K., Pehlke, H., Deininger, M., Jerosch, K. & Brey, T. Scientific Background Document in Support of the Development of a CCAMLR MPA in the Weddell Sea (Antarctica)–Version 2016. (2016).32.Gridley, T., Silva, M., Wilkinson, C., Seakamela, S. & Elwen, S. H. Song recorded near a super-group of humpback whales on a mid-latitude feeding ground off South Africa. J. Acoust. Soc. Am. 143, 298–304 (2018).Article 
ADS 

Google Scholar 
33.Spreen, G., Kaleschke, L. & Heygster, G. Sea ice remote sensing using AMSR-E 89-GHz channels. J. Geophys. Res.-Oceans 113, C02S03 (2008).Article 
ADS 

Google Scholar 
34.Tynan, C. T. & Thiele, D. Report on Antarctic ice edge definition by the ad hoc working group on ice data collection in the Antarctic. Paper: SC/55/19, submitted to the Scientific Committee of the International Whaling Commission (2003).35.Dice, L. R. Measures of the amount of ecologic association between species. Ecology 26, 297–302 (1945).Article 

Google Scholar 
36.Kohonen, T. Median strings. Pattern Recogn. Lett. 3, 309–313 (1985).Article 
ADS 

Google Scholar 
37.Schall, E. et al. Multi-year presence of humpback whales in the Atlantic sector of the Southern Ocean but not during El Niño. Commun. Biol. 4, 1–7 (2021).Article 

Google Scholar 
38.Ritschard, M. & Brumm, H. Zebra finch song reflects current food availability. Evol. Ecol. 26, 801–812 (2012).Article 

Google Scholar 
39.Darling, J. D., Acebes, J. M. V., Frey, O., Urbán, R. J. & Yamaguchi, M. Convergence and divergence of songs suggests ongoing, but annually variable, mixing of humpback whale populations throughout the North Pacific. Sci. Rep. 9, 1–14 (2019).ADS 

Google Scholar 
40.Schall, E. et al. Large-scale spatial variabilities in the humpback whale acoustic presence in the Atlantic sector of the Southern Ocean. R. Soc. Open Sci. 7, 201347 (2020).PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
41.Van Opzeeland, I., Van Parijs, S., Kindermann, L., Burkhardt, E. & Boebel, O. Calling in the cold: Pervasive acoustic presence of humpback whales (Megaptera novaeangliae) in Antarctic coastal waters. PLoS ONE 8, 1–7 (2013).
Google Scholar 
42.Craig, A. S., Herman, L. M., Gabriele, C. M. & Pack, A. A. Migratory timing of humpback whales (Megaptera novaeangliae) in the central north Pacific varies with age, sex and reproductive status. Behaviour 140, 981–1001 (2003).Article 

Google Scholar 
43.Dawbin, W. Temporal segregation of humpback whales during migration in southern hemisphere waters. Mem. Qld. Mus. 42, 105–138 (1997).
Google Scholar 
44.Magnúsdóttir, E., Rasmussen, M., Lammers, M. & Svavarsson, J. Humpback whale songs during winter in subarctic waters. Polar Biol. 37, 427–433 (2014).Article 

Google Scholar 
45.Bombosch, A. et al. Predictive habitat modelling of humpback (Megaptera novaeangliae) and Antarctic minke (Balaenoptera bonaerensis) whales in the Southern Ocean as a planning tool for seismic surveys. Deep Sea Res. Part 1 91, 101–114 (2014).Article 

Google Scholar 
46.Thiele, D. et al. Seasonal variability in whale encounters in the Western Antarctic Peninsula. Deep Sea Research (Part II, Topical Studies in Oceanography) 51, 2311–2325 (2004).Article 
ADS 

Google Scholar 
47.Brown, M. R., Corkeron, P. J., Hale, P. T., Schultz, K. W. & Bryden, M. M. Evidence for a sex-segregated migration in the humpback whale (Megaptera novaeangliae). Proc. R. Soc. Lond. B 259, 229–234 (1995).CAS 
Article 
ADS 

Google Scholar 
48.McDonald, M. A., Mesnick, S. L. & Hildebrand, J. A. Biogeographic characterisation of blue whale song worldwide: Using song to identify populations. J. Cetac. Res. Manage. 8, 55–65 (2006).
Google Scholar 
49.Thomisch, K. et al. Spatio-temporal patterns in acoustic presence and distribution of Antarctic blue whales Balaenoptera musculus intermedia in the Weddell Sea. Endanger. Species Res. 30, 239–253 (2016).Article 

Google Scholar 
50.Oleson, E. M., Širović, A., Bayless, A. R. & Hildebr, J. A. Synchronous seasonal change in fin whale song in the North Pacific. PLoS ONE 9, e115678 (2014).PubMed 
PubMed Central 
Article 
ADS 
CAS 

Google Scholar 
51.Simon, M., Stafford, K. M., Beedholm, K., Lee, C. M. & Madsen, P. T. Singing behavior of fin whales in the Davis Strait with implications for mating, migration and foraging. J. Acoust. Soc. Am. 128, 3200–3210 (2010).PubMed 
Article 
ADS 

Google Scholar 
52.Stafford, K. M. et al. Spitsbergen’s endangered bowhead whales sing through the polar night. Endanger. Species Res. 18, 95–103 (2012).Article 

Google Scholar 
53.Risch, D. et al. Minke whale acoustic behavior and multi-year seasonal and diel vocalization patterns in Massachusetts Bay, USA. Mar. Ecol. Prog. Ser. 489, 279–295 (2013).Article 
ADS 

Google Scholar 
54.Brenowitz, E. A., Margoliash, D. & Nordeen, K. W. An introduction to birdsong and the avian song system. J. Neurobiol. 33, 495–500 (1997).CAS 
PubMed 
Article 

Google Scholar 
55.Tobias, J., Gamarra-Toledo, V., García-Olaechea, D., Pulgarin, P. & Seddon, N. Year-round resource defence and the evolution of male and female song in suboscine birds: Social armaments are mutual ornaments. J. Evol. Biol. 24, 2118–2138 (2011).CAS 
PubMed 
Article 

Google Scholar 
56.Vu, E. T., Clark, C., Catelani, K., Kellar, N. M. & Calambokidis, J. Seasonal blubber testosterone concentrations of male humpback whales (Megaptera novaeangliae). Mar. Mam. Sci. 31, 1258–1264 (2015).Article 

Google Scholar 
57.Yamada, K. & Soma, M. Diet and birdsong: Short-term nutritional enrichment improves songs of adult Bengalese finch males. J. Avian Biol. 47, 865–870 (2016).Article 

Google Scholar 
58.Casagrande, S., Pinxten, R., Zaid, E. & Eens, M. Positive effect of dietary lutein and cholesterol on the undirected song activity of an opportunistic breeder. PeerJ 4, e2512 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
59.Weinrich, M. Humpback whale competitive groups observed on a high-latitude feeding ground. Mar. Mamm. Sci. 11, 251–254 (1995).Article 

Google Scholar 
60.Chittleborough, R. The breeding cycle of the female humpback whale, Megaptera nodosa (Bonnaterre). Mar. Freshw. Res. 9, 1–18 (1958).Article 

Google Scholar 
61.Chittleborough, R. Studies on the ovaries of the humback whale, Megaptera nodosa (bonnaterre), on the western Australian coast. Mar. Freshw. Res. 5, 35–63 (1954).Article 

Google Scholar 
62.Cerchio, S., Jacobsen, J. K. & Norris, T. F. Temporal and geographical variation in songs of humpback whales, Megaptera novaeangliae: Synchronous change in Hawaiian and Mexican breeding assemblages. Anim. Behav. 62, 313–329 (2001).Article 

Google Scholar 
63.Garland, E. C. et al. Quantifying humpback whale song sequences to understand the dynamics of song exchange at the ocean basin scale. J. Acoust. Soc. Am. 133, 560–569 (2013).PubMed 
Article 
ADS 

Google Scholar 
64.Allen, J. A., Garland, E. C., Dunlop, R. A. & Noad, M. J. Cultural revolutions reduce complexity in the songs of humpback whales. Proc. R. Soc. B 285, 20182088 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
65.Stevick, P. T. et al. Migrations of individually identified humpback whales between the Antarctic Peninsula and South America. J. Cetac. Res. Manag. 6, 109–113 (2004).
Google Scholar 
66.Engel, M. et al. Mitochondrial DNA diversity of the Southwestern Atlantic humpback whale (Megaptera novaeangliae) breeding area off Brazil, and the potential connections to Antarctic feeding areas. Conserv. Genet. 9, 1253–1262 (2008).CAS 
Article 

Google Scholar 
67.Amaral, A. R. et al. Population genetic structure among feeding aggregations of humpback whales in the Southern Ocean. Mar. Biol. 163, 1–13 (2016).Article 

Google Scholar 
68.Rekdahl, M. L. et al. Culturally transmitted song exchange between humpback whales (Megaptera novaeangliae) in the southeast Atlantic and southwest Indian Ocean basins. R. Soc. Open Sci. 5, 172305 (2018).PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
69.Darling, J. D. & Sousa-Lima, R. S. Songs indicate interaction between humpback whale (Megaptera novaeangliae) populations in the western and eastern South Atlantic Ocean. Mar. Mamm. Sci. 21, 557–566 (2005).Article 

Google Scholar 
70.Razafindrakoto, Y., Cerchio, S., Collins, T., Rosenbaum, H. & Ngouessono, S. Similarity of humpback whale song from Madagascar and Gabon indicates significant contact between South Atlantic and southwest Indian Ocean populations. PLoS ONE 8, e79422 (2009).
Google Scholar 
71.Zerbini, A. et al. Migration and summer destinations of humpback whales (Megaptera novaeangliae) in the western South Atlantic Ocean. J. Cetac. Res. Manag. 3, 113–118 (2011).
Google Scholar 
72.Rosenbaum, H. C., Maxwell, S. M., Kershaw, F. & Mate, B. Long-range movement of humpback whales and their overlap with anthropogenic activity in the South Atlantic Ocean. Conserv. Biol. 28, 604–615 (2014).PubMed 
Article 

Google Scholar 
73.Filun, D. et al. Frozen verses: Antarctic minke whales (Balaenoptera bonaerensis) call predominantly during austral winter. R. Soc. Open Sci. 7, 192112 (2020).PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
74.Rettig, S. et al. In International Confence and Exhibition on Underwater Acoustics. (eds Papadakis, J. & Bjorno, L.) 1669–1674.75.Baumgartner, M. F. & Mussoline, S. E. A generalized baleen whale call detection and classification system. J. Acoust. Soc. Am. 129, 2889–2902 (2011).PubMed 
Article 
ADS 

Google Scholar 
76.Klinck, H. et al. Long-range underwater vocalizations of the crabeater seal (Lobodon carcinophaga). J. Acoust. Soc. Am. 128, 474–479 (2010).PubMed 
Article 
ADS 

Google Scholar 
77.Risch, D. et al. Mysterious bio-duck sound attributed to the Antarctic minke whale (Balaenoptera bonaerensis). Biol. Lett. 10, 20140175 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
78.Schall, E. & Van Opzeeland, I. Calls produced by Ecotype C killer whales (Orcinus orca) off the Eckstrom iceshelf, Antarctica. Aquat. Mamm. 43, 117–126 (2017).Article 

Google Scholar 
79.Van Opzeeland, I. et al. Acoustic ecology of Antarctic pinnipeds. Mar. Ecol. Prog. Ser. 414, 267–291 (2010).Article 
ADS 

Google Scholar 
80.Dunlop, R. A., Cato, D. H. & Noad, M. J. Non-song acoustic communication in migrating humpback whales (Megaptera novaeangliae). Mar. Mamm. Sci. 24, 613–629 (2008).Article 

Google Scholar 
81.Stimpert, A. K., Au, W. W., Parks, S. E., Hurst, T. & Wiley, D. N. Common humpback whale (Megaptera novaeangliae) sound types for passive acoustic monitoring. J. Acoust. Soc. Am. 129, 476–482 (2011).PubMed 
Article 
ADS 

Google Scholar 
82.Cavalieri, D., Parkinson, C., Gloersen, P. & Zwally, H. Sea Ice Concentrations from Nimbus-7 SMMR and DMSP SSM/I-SSMIS Passive Microwave Data, Version 1 (NASA National Snow and Ice Data Center Distributed Active Archive Center, 1996).
Google Scholar 
83.Greene, C. A., Gwyther, D. E. & Blankenship, D. D. Antarctic mapping tools for MATLAB. Comput. Geosci. 104, 151–157 (2017).Article 
ADS 

Google Scholar 
84.Greene, C. A. Daily Antarctic Sea Ice Concentration (2020).85.Schall, E., Roca, I. & Van Opzeeland, I. Acoustic metrics to assess humpback whale song unit structure from the Atlantic sector of the Southern ocean. J. Acoust. Soc. Am. 149, 4649–4658 (2021).PubMed 
Article 
ADS 

Google Scholar 
86.Dalla Rosa, L., Secchi, E., Maia, Y. G., Zerbini, A. & Heide-Jørgensen, M. Movements of satellite-monitored humpback whales on their feeding ground along the Antarctic Peninsula. Polar Biol. 31, 771–781 (2008).Article 

Google Scholar 
87.Zann, R. & Cash, E. Developmental stress impairs song complexity but not learning accuracy in non-domesticated zebra finches (Taeniopygia guttata). Behav. Ecol. Sociobiol. 62, 391–400 (2008).Article 

Google Scholar 
88.Woodgate, J. L., Mariette, M. M., Bennett, A. T., Griffith, S. C. & Buchanan, K. L. Male song structure predicts reproductive success in a wild zebra finch population. Anim. Behav. 83, 773–781 (2012).Article 

Google Scholar 
89.Boogert, N. J., Giraldeau, L.-A. & Lefebvre, L. Song complexity correlates with learning ability in zebra finch males. Anim. Behav. 76, 1735–1741 (2008).Article 

Google Scholar 
90.Templeton, C. N., Laland, K. N. & Boogert, N. J. Does song complexity correlate with problem-solving performance in flocks of zebra finches?. Anim. Behav. 92, 63–71 (2014).Article 

Google Scholar 
91.R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2018). https://www.R-project.org/.92.Suzuki, R., Terada, Y. & Shimodaira, H. pvclust: Hierarchical Clustering with P-values via Multiscale Bootstrap Resampling. R Package Version 2.2–0 (2019).93.Garland, E. C. et al. Improved versions of the Levenshtein distance method for comparing sequence information in animals’ vocalisations: Tests using humpback whale song. Behaviour 149, 1413–1441 (2012).Article 

Google Scholar 
94.Van der Loo, M. P. The stringdist package for approximate string matching. R J. 6, 111–122 (2014).Article 

Google Scholar 
95.Pawlowicz, R. M_Map: A Mapping Package for MATLAB v. Version 1.4m. www.eoas.ubc.ca/~rich/map.html (2020). More

1.Takemura, A., Rahman, M. S. & Park, Y. J. External and internal controls of lunar-related reproductive rhythms in fishes. J. Fish Biol. 76, 7–26 (2010).CAS 
PubMed 
Article 

Google Scholar 
2.Brockington, S. & Clarke, A. The relative influence of temperature and food on the metabolism of a marine invertebrate. J. Exp. Mar. Bio. Ecol. 258, 87–99 (2001).CAS 
Article 

Google Scholar 
3.Kelly, M. S. Environmental parameters controlling gametogenesis in the echinoid Psammechinus miliaris. J. Exp. Mar. Bio. Ecol. 266, 67–80 (2001).Article 

Google Scholar 
4.Muthiga, N. A. The reproductive biology of a new species of sea cucumber, Holothuria (Mertensiothuria) arenacava in a Kenyan marine protected area: The possible role of light and temperature on gametogenesis and spawning. Mar. Biol. 149, 585–593 (2006).Article 

Google Scholar 
5.Emilio, L. et al. Is the Orton’s rule still valid? Tropical sponge fecundity, rather than periodicity, is modulated by temperature and other proximal cues. Hydrobiologia 815, 187–205 (2018).Article 

Google Scholar 
6.St.Gelais, A. T., Chaves-Fonnegra, A., Moulding, A. L., Kosmynin, V. N. & Gilliam, D. S. Siderastrea siderea spawning and oocyte resorption at high latitude. Invertebr. Reprod. Dev. 60, 212–222 (2016).Article 

Google Scholar 
7.Zhadan, P. M., Vaschenko, M. A. & Ryazanov, S. D. Assessing the effect of environmental factors on the spawning activity of the sea urchin Strongylocentrotus intermedius through video recording observations. Mar. Ecol. Prog. Ser. 588, 101–119 (2018).CAS 
Article 
ADS 

Google Scholar 
8.Grange, L. J., Tyler, P. A., Peck, L. S. & Cornelius, N. Long-term interannual cycles of the gametogenic ecology of the Antarctic brittle star Ophionotus victoriae. Mar. Ecol. Prog. Ser. 278, 141–155 (2004).Article 
ADS 

Google Scholar 
9.Balogh, R., Wolfe, K. & Byrne, M. Gonad development and spawning of the vulnerable commercial sea cucumber, Stichopus herrmanni, in the southern Great Barrier Reef. J. Mar. Biol. Assoc. United Kingdom 99, 487–495 (2019).Article 

Google Scholar 
10.Stenseth, N. C. et al. Studying climate effects on ecology through the use of climate indices: The North Atlantic Oscillation, El Niño Southern Oscillation and beyond. Proc. R. Soc. B Biol. Sci. 270, 2087–2096 (2003).Article 

Google Scholar 
11.Wood, S. et al. El Nino and coral larval dispersal across the eastern Pacific marine barrier. Nat. Commun. 7, 1 (2016).
Google Scholar 
12.Turner, J. The El Niño-Southern Oscillation and Antarctica. Int. J. Climatol. 24, 1–31 (2004).Article 

Google Scholar 
13.La, H. S. et al. Zooplankton and micronekton respond to climate fluctuations in the Amundsen Sea polynya, Antarctica.. Sci. Rep. 9, 1–7 (2019).CAS 
Article 
ADS 

Google Scholar 
14.Xuebin, Z. & Mcphaden, M. J. Eastern equatorial Pacific forcing of ENSO sea surface temperature anomalies. J. Clim. 21, 6070–6079 (2008).Article 
ADS 

Google Scholar 
15.Oliver, E. C. J. et al. Longer and more frequent marine heatwaves over the past century. Nat. Commun. 9, 1–12 (2018).CAS 
Article 

Google Scholar 
16.Ryan, J. P. et al. Causality of an extreme harmful algal bloom in Monterey Bay, California, during the 2014–2016 northeast Pacific warm anomaly. Geophys. Res. Lett. 44, 5571–5579 (2017).Article 
ADS 

Google Scholar 
17.Conde, A. & Prado, M. Changes in phytoplankton vertical distribution during an El Niño event. Ecol. Indic. 90, 201–205 (2018).Article 

Google Scholar 
18.Santidrián Tomillo, P. et al. The impacts of extreme El Niño events on sea turtle nesting populations. Clim. Change https://doi.org/10.1007/s10584-020-02658-w (2020).Article 

Google Scholar 
19.Wilson, S. K. et al. Climatic forcing and larval dispersal capabilities shape the replenishment of fishes and their habitat-forming biota on a tropical coral reef. Ecol. Evol. 8, 1918–1928 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
20.Welhouse, L., Lazzara, M., Keller, L., Tripoli, G. & Hitchman, M. Composite analysis of the effects of ENSO events on Antarctica. J. Clim. 29, 1797–1808 (2016).Article 
ADS 

Google Scholar 
21.Testa, J. W. et al. Temporal variability in Antarctic marine ecosystems: periodic fluctuations in the phocid seals. Can. J. Fish. Aquat. Sci. 48, 631–639 (1991).Article 

Google Scholar 
22.Román-González, A. et al. Analysis of ontogenetic growth trends in two marine Antarctic bivalves Yoldia eightsi and Laternula elliptica: Implications for sclerochronology. Palaeogeogr. Palaeoclimatol. Palaeoecol. 465, 300–306 (2017).Article 

Google Scholar 
23.Brown, M. et al. Long-term effect of photoperiod, temperature and feeding regimes on the respiration rates of Antarctic Krill (Euphausia superba). Open J. Mar. Sci. 3, 40–51 (2013).Article 

Google Scholar 
24.Ainley, D. G. et al. Decadal trends in abundance, size and condition of Antarctic toothfish in McMurdo Sound, Antarctica, 1972–2011. Fish Fish. 14, 343–363 (2013).Article 

Google Scholar 
25.Doney, S. C. et al. Climate Change Impacts on Marine Ecosystems. Ann. Rev. Mar. Sci. 4, 11–37 (2012).PubMed 
Article 

Google Scholar 
26.Peck, L. S. Antarctic Marine Biodiversity: Adaptations, Environments and Responses to Change. Oceanogr. Mar. Biol. An Annu. Rev. 56, 105–236 (2018).Article 

Google Scholar 
27.Peck, L. S. A Cold Limit to Adaptation in the Sea. Trends Ecol. Evol. 31, 13–26 (2016).PubMed 
Article 

Google Scholar 
28.Brockington, S., Peck, L. S. & Tyler, P. A. Gametogenesis and gonad mass cycles in the common circumpolar Antarctic echinoid Sterechinus neumayeri. Mar. Ecol. Prog. Ser. 330, 139–147 (2007).Article 
ADS 

Google Scholar 
29.Grange, L. J., Tyler, P. A. & Peck, L. S. Multi-year observations on the gametogenic ecology of the Antarctic seastar Odontaster validus. Mar. Biol. 153, 15–23 (2007).Article 

Google Scholar 
30.Brockington, S. The seasonal ecology and physiology of Sterechinus neumayeri (Echinodermata; Echinoidea) at Adelaide Island, Antarctica. PhD thesis The Open University. (2001).31.Bosch, I., Beauchamp, K. A., Steele, M. E. & Pearse, J. S. Development, metamorphosis, and seasonal abundance of embryos and larvae of the Antarctic sea urchin Sterechinus Neumayeri. Biol. Bull. 173, 126–135 (1987).PubMed 
Article 

Google Scholar 
32.Stanwell-Smith, D. & Peck, L. S. Temperature and embryonic development in relation to spawning and field occurrence of larvae of three Antarctic echinoderms. Biol. Bull. 194, 44–52 (1998).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
33.Fogt, R. L., Bromwich, D. H. & Hines, K. M. Understanding the SAM influence on the South Pacific ENSO teleconnection. Clim. Dyn. 36, 1555–1576 (2011).Article 

Google Scholar 
34.Kwok, R. & Comiso, J. C. Spatial patterns of variability in Antarctic surface temperature: Connections to the Southern Hemisphere Annular Mode and the Southern Oscillation. Geophys. Res. Lett. 29, 2–5 (2002).
Google Scholar 
35.Santamaría-Del-ángel, E. et al. Interannual climate variability in the west antarctic peninsula under austral summer conditions. Remote Sens. 13, 1 (2021).Article 

Google Scholar 
36.Montgomery, D. & Peck, E. Introduction to linear regression analysis. (Wiley, 1992).37.Halberg, F., Shankaraiah, K. & Giese, A. The chronobiology of marine invertebrates: methods of analysis. in Reproduction of marine invertebrates, Vol IX. General aspects: seeking unity in diversity 331–384 (The Boxwood Press, 1987).38.Loeb, V. J., Hofmann, E. E., Klinck, J. M., Holm-Hansen, O. & White, W. B. ENSO and variability of the antarctic peninsula pelagic marine ecosystem. Antarct. Sci. 21, 135–148 (2009).Article 
ADS 

Google Scholar 
39.White, W. B., Chen, S. C., Allan, R. J. & Stone, R. C. Positive feedbacks between the Antarctic Circumpolar Wave and the global El Niño-Southern Oscillation wave. J. Geophys. Res. C Ocean. 107, 29–31 (2002).
Google Scholar 
40.Saba, G. K. et al. Winter and spring controls on the summer food web of the coastal West Antarctic Peninsula. Nat. Commun. 5, 1–8 (2014).CAS 

Google Scholar 
41.Cavanagh, R. D. et al. A synergistic approach for evaluating climate model output for ecological applications. Front. Mar. Sci. 4, 1 (2017).Article 

Google Scholar 
42.Vergani, D. F., Labraga, J. C., Stanganelli, Z. B. & Dunn, M. The effects of El Niño-La Niña on reproductive parameters of elephant seals feeding in the Bellingshausen Sea. J. Biogeogr. 35, 248–256 (2008).Article 

Google Scholar 
43.Clark, G. F. et al. Light-driven tipping points in polar ecosystems. Glob. Chang. Biol. 19, 3749–3761 (2013).PubMed 
Article 
ADS 

Google Scholar 
44.Schneider, D. P., Okumura, Y. & Deser, C. Observed Antarctic interannual climate variability and tropical linkages. J. Clim. 25, 4048–4066 (2012).Article 
ADS 

Google Scholar 
45.Yuan, X. ENSO-related impacts on Antarctic sea ice: A synthesis of phenomenon and mechanisms. Antarct. Sci. 16, 415–425 (2004).Article 
ADS 

Google Scholar 
46.Loeb, V. J. & Santora, J. A. Population dynamics of Salpa thompsoni near the Antarctic Peninsula: Growth rates and interannual variations in reproductive activity (1993–2009). Prog. Oceanogr. 96, 93–107 (2012).Article 
ADS 

Google Scholar 
47.Moran, A. L., McAlister, J. S. & Whitehill, E. A. G. Eggs as energy: Revisiting the scaling of egg size and energetic content among echinoderms. Biol. Bull. 224, 184–191 (2013).CAS 
PubMed 
Article 

Google Scholar 
48.Gómez-Robles, E. & Saucedo, P. E. Evaluation of quality indices of the gonad and somatic tissues involved in reproduction of the pearl oyster Pinctada mazatlanica with histochemistry and digital image analysis. J. Shellfish Res. 28, 329–335 (2009).Article 

Google Scholar 
49.Gómez-Valdez, M., Ocampo, L., Carvalho-Saucedo, L. & Gutiérrez-González, J. Reproductive activity and seasonal variability in the biochemical composition of a pen shell, Atrina maura.. Mar. Ecol. Prog. Ser. 663, 99–113 (2021).Article 
ADS 
CAS 

Google Scholar 
50.Steinberg, D. K. et al. Long-term (1993–2013) changes in macrozooplankton off the Western Antarctic Peninsula. Deep. Res. Part I Oceanogr. Res. Pap. 101, 54–70 (2015).Article 
ADS 

Google Scholar 
51.Rozema, P. D. et al. Interannual variability in phytoplankton biomass and species composition in northern Marguerite Bay (West Antarctic Peninsula) is governed by both winter sea ice cover and summer stratification. Limnol. Oceanogr. 62, 235–252 (2017).Article 
ADS 

Google Scholar 
52.Starr, M., Himmelman, J. H. & Therriault, J. Direct coupling of marine invertebrate spawning with phytoplankton blooms. Science 247, 1071–1074 (1990).CAS 
PubMed 
Article 
ADS 

Google Scholar 
53.Harrington, L. H., Walker, C. W. & Lesser, M. P. Stereological analysis of nutritive phagocytes and gametogenic cells during the annual reproductive cycle of the green sea urchin, Strongylocentrotus droebachiensis.. Invertebr. Biol. 126, 202–209 (2007).Article 

Google Scholar 
54.Magniez, P. Reproductive cycle of the brooding echinoid Abatus cordatus (Echinodermata) in Kerguelen (Antarctic Ocean): changes in the organ indices, biochemical composition and caloric content of the gonads. Mar. Biol. 74, 55–64 (1983).CAS 
Article 

Google Scholar 
55.Pérez, A. F., Morriconi, E., Boy, C. & Calvo, J. Seasonal changes in energy allocation to somatic and reproductive body components of the common cold temperature sea urchin Loxechinus albus in a Sub-Antarctic environment. Polar Biol. 31, 443–449 (2008).Article 

Google Scholar 
56.Hernandez, E., Vázquez, O. A., Torruco, A. & Rahman, M. S. Reproductive cycle and gonadal development of the Atlantic sea urchin Arbacia punctulata in the Gulf of Mexico: changes in nutritive phagocytes in relation to gametogenesis. Mar. Biol. Res. 16, 177–194 (2020).Article 

Google Scholar 
57.Bronstein, O., Kroh, A. & Loya, Y. Reproduction of the long-spined sea urchin Diadema setosum in the Gulf of Aqaba – Implications for the use of gonad-indexes. Sci. Rep. 6, 1–11 (2016).Article 
CAS 

Google Scholar 
58.Alturkistani, H. A., Tashkandi, F. M. & Mohammedsaleh, Z. M. Histological Stains: A Literature Review and Case Study. Glob. J. Health Sci. 8, 72–79 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
59.Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).CAS 
PubMed 
Article 

Google Scholar 
60.Rueden, C. T. et al. Image J2: ImageJ for the next generation of scientific image data. BMC Bioinformatics 18, 1–26 (2017).Article 
ADS 

Google Scholar 
61.Lau, S. C. Y., Grange, L. J., Peck, L. S. & Reed, A. J. The reproductive ecology of the Antarctic bivalve Aequiyoldia eightsii (Protobranchia: Sareptidae) follows neither Antarctic nor taxonomic patterns. Polar Biol. 41, 1693–1706 (2018).Article 

Google Scholar 
62.Reed, A. J., Morris, J. P., Linse, K. & Thatje, S. Reproductive morphology of the deep-sea protobranch bivalves Yoldiella ecaudata, Yoldiella sabrina, and Yoldiella valettei (Yoldiidae) from the Southern Ocean. Polar Biol. 37, 1383–1392 (2014).Article 

Google Scholar 
63.Cleveland, W. S. Robust locally weighted regression and smoothing scatterplots. J. Am. Stat. Assoc. 74, 829–836 (1979).MathSciNet 
MATH 
Article 

Google Scholar 
64.Venables, H. J., Clarke, A. & Meredith, M. P. Wintertime controls on summer stratification and productivity at the western Antarctic Peninsula. Limnol. Oceanogr. 58, 1035–1047 (2013).Article 
ADS 

Google Scholar 
65.Clarke, A., Meredith, M. P., Wallace, M. I., Brandon, M. A. & Thomas, D. N. Seasonal and interannual variability in temperature, chlorophyll and macronutrients in northern Marguerite Bay, Antarctica.. Deep Res. Part II Top. Stud. Oceanogr. 55, 198–206 (2008).
Google Scholar 
66.Zuur, A., Ieno, E. N. & Smith, G. M. Analyzing Ecological Data. in Analyzing Ecological Data (ed. M. Gail, K. Krickeberg, J. Samet, A. Tsiatis, W. W.) 23–47 (Springer-Verlag New York, 2007).67.Burnham, K. P. & Anderson, D. R. Model selection and multimodel inference. A practical information-theoretical approach. Model Selection and Multimodel Inference (Springer, 2002). https://doi.org/10.1007/978-0-387-22456-5_768.Fisher, R., Wilson, S. K., Sin, T. M., Lee, A. C. & Langlois, T. J. A simple function for full-subsets multiple regression in ecology with R. Ecol. Evol. 8, 6104–6113 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
69.Wood, S. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. R. Stat. Soc. 73, 3–36 (2011).MathSciNet 
MATH 
Article 

Google Scholar 
70.De Leij, R., Peck, L. S. & Grange, L. J. R code and csv. files. https://doi.org/10.5061/dryad.6q573n5z1 (2021).71.Grange, L. J., Peck, L. S. & Tyler, P. A. Reproductive ecology of the circumpolar Antarctic nemertean Parborlasia corrugatus: No evidence for inter-annual variation. J. Exp. Mar. Bio. Ecol. 404, 98–107 (2011).Article 

Google Scholar  More

1.Lowe, S., Browne, M., Boudjelas, S. & De Poorter, M. 100 of the world’s worst invasive alien species: A selection from the global invasive species database. In Encyclopedia of Biological Invasions 12 (The Invasive Species Specialist Group (ISSG), Species Survival Commission (SSC), World Conservation Union (IUCN), 2000). https://doi.org/10.1525/9780520948433-159.2.North American Invasive Species Network. The ten most important invasive species or invasive species assemblages in North America in 2015. https://www.bugwoodcloud.org/mura/naisn/assets/File/NAISNPRJan2015.pdf (2015).3.Keuling, O. et al. Eurasian wild boar Sus scrofa (Linnaeus, 1758). in Ecology, Conservation and Management of Wild Pigs and Peccaries (eds. Melleti, M. & Meijaard, E.) 202–233 (Cambridge University Press, 2017).4.Strickland, B. K., Smith, M. D. & Smith, A. L. Wild pig damage to resources. In Invasive Wild Pigs in North America: Ecology, Impacts, and Management (eds VerCauteren, K. C. et al.) 143–174 (RC Press, London, 2020).
Google Scholar 
5.Pimental, D. Environmental and economic costs of vertebrate species invasions into the United States. In Managing Vertebrate Invasive Species: Proceedings of an International Symposium (eds. Witmer, G. W., Pitt, W. C. & Fagerstone, K. A.) 2–8 (USDA National Wildlife Research Center, Fort Collins, CO, USA, 2007).6.Ditchkoff, S. S. & Bodenchuk, M. J. Management of wild pigs. In Invasive Wild Pigs in North America: Ecology, Impacts, and Management (eds VerCauteren, K. C. et al.) 175–198 (CRC Press, London, 2020).
Google Scholar 
7.Maher, C. R. & Lott, D. F. Definitions of territoriality used in the study of variation in vertebrate spacing systems. Anim. Behav. 49, 1581–1597 (1995).Article 

Google Scholar 
8.Bastille-Rousseau, G. et al. Multi-level movement response of invasive wild pigs (Sus scrofa) to removal. Pest Manag. Sci. 77, 85–95 (2021).CAS 
Article 

Google Scholar 
9.Boitani, L., Mattei, L., Nonis, D. & Corsi, F. Spatial and activity patterns of wild boars in Tuscany, Italy. J. Mammal. 75, 600–612 (1994).Article 

Google Scholar 
10.Ilse, L. M. & Hellgren, E. C. Resource partitioning in sympatric populations of collared peccaries and feral hogs in southern Texas. J. Mammal. 76, 784–799 (1995).Article 

Google Scholar 
11.Gabor, T. M., Hellgren, E. C., Bussche, R. A. V. D. & Silvy, N. J. Demography, sociospatial behaviour and genetics of feral pigs (Sus scrofa) in a semi-arid environment. J. Zool. 247, 311–322 (1999).Article 

Google Scholar 
12.Sparklin, B. D., Mitchell, M. S., Hanson, L. B., Jolley, D. B. & Ditchkoff, S. S. Territoriality of feral pigs in a highly persecuted population on Fort Benning, Georgia. J. Wildl. Manag. 73, 497–502 (2009).Article 

Google Scholar 
13.Beasley, J. C., Ditchkoff, S. S., Mayer, J. J., Smith, M. D. & VerCauteren, K. C. Research priorities for managing invasive wild pigs in North America. J. Wildl. Manag. 82, 674–681 (2018).Article 

Google Scholar 
14.Gray, S. M., Roloff, G. J., Montgomery, R. A., Beasley, J. C. & Pepin, K. M. Wild pig spatial ecology and behavior. In Invasive Wild Pigs in North America: Ecology, Impacts, and Management (eds VerCauteren, K. C. et al.) 33–56 (CRC Press, London, 2020).
Google Scholar 
15.Emlen, J. T. Defended area? A critique of the territory concept and of conventional thinking. Ibis 99, 352 (1957).
Google Scholar 
16.Kamath, A. & Wesner, A. B. Animal territoriality, property and access: A collaborative exchange between animal behaviour and the social sciences. Anim. Behav. 164, 233–239 (2020).Article 

Google Scholar 
17.ESRI. ArcGIS Pro. Environmental Systems Research Institute (2021).18.Mayer, J. J. Wild hog. In Ecology and Management of a Forested Landscape: Fifty Years on the Savannah River Site (eds Kilgo, J. C. & Blake, J. I.) 374–379 (Island Press, Washington, 2005).
Google Scholar 
19.Mayer, J. J., Edwards, T. B., Garabedian, J. E. & Kilgo, J. C. Sanitary waste landfill effects on an invasive wild pig population. J. Wildl. Manag. 85, 868–879 (2021).Article 

Google Scholar 
20.Royle, J. A., Chandler, R. B., Sollmann, R. & Gardner, B. Spatial Capture-Recapture (Academic Press, Cambridge, 2014).
Google Scholar 
21.Kranstauber, B., Kays, R., LaPoint, S. D., Wikelski, M. & Safi, K. A dynamic Brownian bridge movement model to estimate utilization distributions for heterogeneous animal movement. J. Anim. Ecol. 81, 738–746 (2012).Article 

Google Scholar 
22.Byrne, M. E., Guthrie, J. D., Hardin, J., Collier, B. A. & Chamberlain, M. J. Evaluating wild Turkey movement ecology: An example using first-passage time analysis. Wildl. Soc. Bull. 38, 407–413 (2014).Article 

Google Scholar 
23.Clontz, L. M., Pepin, K. M., VerCauteren, K. C. & Beasley, J. C. Behavioral state resource selection in invasive wild pigs in the Southeastern United States. Sci. Rep. 11, 6924 (2021).CAS 
Article 
ADS 

Google Scholar 
24.White, G. C. & Garrott, R. A. Analysis of Wildlife Radio-Tracking Data (Academic Press, Cambridge, 1990).
Google Scholar 
25.Potts, J. R., Harris, Stephen & Giuggioli, L. Quantifying behavioral changes in territorial animals caused by sudden population declines. Am. Nat. 182, E73–E82 (2013).Article 

Google Scholar 
26.Fieberg, J. & Kochanny, C. O. Quantifying home-range overlap: The importance of the utilization distribution. J. Wildl. Manag. 69, 1346–1359 (2005).Article 

Google Scholar 
27.R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing (2021).28.Schielzeth, H. & Forstmeier, W. Conclusions beyond support: Overconfident estimates in mixed models. Behav. Ecol. 20, 416–420 (2009).Article 

Google Scholar 
29.Kay, S. L. et al. Quantifying drivers of wild pig movement across multiple spatial and temporal scales. Mov. Ecol. 5, 14 (2017).Article 

Google Scholar 
30.Hurvich, C. M. & Tsai, C.-L. Regression and time series model selection in small samples. Biometrika 76, 297–307 (1989).MathSciNet 
Article 

Google Scholar 
31.Long, J. A., Nelson, T. A., Webb, S. L. & Gee, K. L. A critical examination of indices of dynamic interaction for wildlife telemetry studies. J. Anim. Ecol. 83, 1216–1233 (2014).Article 

Google Scholar 
32.Benhamou, S., Valeix, M., Chamaillé-Jammes, S., Macdonald, D. W. & Loveridge, A. J. Movement-based analysis of interactions in African lions. Anim. Behav. 90, 171–180 (2014).Article 

Google Scholar 
33.Brotherton, P. N. M., Pemberton, J. M., Komers, P. E. & Malarky, G. Genetic and behavioural evidence of monogamy in a mammal, Kirk’s dik–dik (Madoqua kirkii). Proc. R. Soc. Lond. B Biol. Sci. 264, 675–681 (1997).CAS 
Article 
ADS 

Google Scholar 
34.Burt, W. H. Territoriality and home range concepts as applied to mammals. J. Mammal. 24, 346–352 (1943).Article 

Google Scholar 
35.Cooper, N. W., Sherry, T. W. & Marra, P. P. Modeling three-dimensional space use and overlap in birds. Auk 131, 681–693 (2014).Article 

Google Scholar 
36.Millspaugh, J. J., Gitzen, R. A., Kernohan, B. J., Larson, M. A. & Clay, C. L. Comparability of three analytical techniques to assess joint space use. Wildl. Soc. Bull. 32, 148–157 (2004).Article 

Google Scholar 
37.Pepin, K. M. et al. Contact heterogeneities in feral swine: Implications for disease management and future research. Ecosphere 7, e01230 (2016).Article 

Google Scholar 
38.Yang, A. et al. Effects of social structure and management on risk of disease establishment in wild pigs. J. Anim. Ecol. 90, 820–833 (2021).Article 

Google Scholar 
39.Carpenter, F. L. Food abundance and territoriality: To defend or not to defend?. Am. Zool. 27, 387–399 (1987).Article 

Google Scholar 
40.Both, C. & Visser, M. E. Density dependence, territoriality, and divisibility of resources: From optimality models to population processes. Am. Nat. 161, 326–336 (2003).Article 

Google Scholar 
41.Doncaster, C. P. & Macdonald, D. W. Optimum group size for defending heterogenous distributions of resources: A model applied to red foxes, Vulpes vulpes, Oxford city. J. Theor. Biol. 159, 189–198 (1992).Article 
ADS 

Google Scholar 
42.Krause, J. & Ruxton, G. D. Living in Groups (University Press, Oxford, 2002).
Google Scholar 
43.Garabedian, J. E., Moorman, C. E., Peterson, M. N. & Kilgo, J. C. Effects of group size and group density on trade-offs in resource selection by a group-territorial central-place foraging woodpecker. Ibis 162, 477–491 (2020).Article 

Google Scholar  More

Automatic feeder productivityTable 1 and Figs. 4, 5 and 6 show the automatic feeder productivity as affected by the different feed pellets sizes (1, 2 and 3 mm), air flow rates (10, 15 and 20 m3 min−1) and rotational speeds of screw (180, 360, 540, 720 and 900 rpm). The results indicate that the automatic feeder productivity increases with increasing feed pellets size, air flow rate and rotational speed of screw. It indicates that when the feed pellets size increased from 1 to 3 mm, the automatic feeder productivity significantly increased from 11.16 to 13.87 (by 19.54%) kg min−1. It also indicates that when the air flow rate increased from 10 to 20 m3 min−1, the automatic feeder productivity significantly increased from 11.02 to 14.03 (by 21.45%) kg min−1, while the automatic feeder productivity significantly increased from 3.33 to 21.46 (by 84.48%) kg min−1 when the rotational speed of screw increased from 180 to 900 rpm.Table 1 Automatic feeder productivity at different feed pellets sizes, air flow rates and rotational speeds of screw.Full size tableFigure 4Automatic feeder productivity at different feed pellet sizes and rotational speeds of screw.Full size imageFigure 5Automatic feeder productivity at different feed pellet sizes and air flow rates.Full size imageFigure 6Automatic feeder productivity at different rotational speeds of screw and flow rates.Full size imageIt could be noticed that increasing the feed pellets size from 1 to 3 mm, tends to increase the automatic feeder productivity from 3.04 to 3.79, 6.23 to 8.92, 11.86 to 14.10, 15.27 to 18.94 and 19.42 to 23.62 kg min−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively. The results also indicate that the automatic feeder productivity increased from 3.04 to 19.42, 3.16 to 21.36 and 3.79 to 23.62 kg min−1 at 1, 2 and 3 mm feed pellets sizes, respectively when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 4.From statistical analysis, there were no significant different between feed pellets sizes 1 and 2 on the automatic feeder productivity, meanwhile, there were significant differences between feed pellets size 3 and sizes 1 and 2 on the productivity. Regarding the effect of air flow rate, there were significant differences between air flow rates on the automatic feeder productivity, the same trend was happened with the effect of rotational speed of screw on productivity. The analysis showed also that the interaction between both ABC was non-significant. On the other hand, the interaction between the effect of both AB, AC and BC on the data was significant as shown in Table 1.Regarding the effect of feed pellet size and air flow rate on the automatic feeder productivity, the results indicate that the automatic feeder productivity increases with increasing the feed pellets size and flow rate. It increased from 9.53 to 12.37, 11.23 to 13.82 and 12.73 to 15.43 kg min−1 for 10, 15 and 20 m3 min−1 air flow rate, respectively, when the feed pellets size increased from 1 to 3 mm. The results also indicate that the automatic feeder productivity increased from 9.53 to 12.73, 11.16 to 13.92 and 12.37 to 15.43 kg min−1 at 1, 2 and 3 mm feed pellets size, respectively, when the air flow rate increased from 10 to 20 m3 min−1 as shown in Fig. 5.The results also indicate that the automatic feeder productivity increased from 2.26 to 4.54, 6.39 to 8.90, 11.76 to 14.56, 15.25 to 18.68 and 19.44 to 23.45 kg min−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively, when the air flow rate increased from 10 to 20 m3 min−1. The results also indicate that the automatic feeder productivity increased from 2.26 to 19.44, 3.19 to 21.50 and 4.54 to 23.45 kg min−1 at 10, 15 and 20 m3 min−1 air flow rate, respectively, when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 6.Multiple regression analysis was carried out to obtain a relationship between the automatic feeder productivity as dependent variable and different of feed pellets size, air flow rate and rotational speed of screw as independent variables. The best fit for this relationship is presented in the following equation:-$$ Pr_{actual} = – 8.457 + 1.354PS + 0.301FR + 0.025RS{text{ R}}^{{2}} = 0.98{ ,} $$
(13)
where PS is the feed pellets size, mm; FR is the air flow rate, m3 min−1; RS is the rotational speed of screw, rpm.This equation could be applied in the range of 1 to 3 mm feed pellets size, 10 to 20 m3 min−1 air flow rate and from 180 to 900 rpm of rotational speed of screw.Automatic feeder efficiencyTable 2, Figs. 7, 8 and 9 show the automatic feeder efficiency as affected by the different feed pellets sizes (1, 2 and 3 mm), air flow rates (10, 15 and 20 m3 min−1) and rotational speeds of screw (180, 360, 540, 720 and 900 rpm). The results indicate that, when the feed pellets size increased from 1 to 3 mm, the automatic feeder efficiency significantly increased from 65.30 to 82.14 (by 20.50%) %. It also indicates that when the air flow rate increased from 10 to 20 m3 min−1, the automatic feeder efficiency significantly increased from 62.58 to 85.07 (by 26.44%) %, while the automatic feeder efficiency significantly increased from 61.58 to 78.69 (by 21.74%) % when the rotational speed of screw increased from 180 to 900 rpm.Table 2 Automatic feeder efficiency at different feed pellets sizes, air flow rates and rotational speeds of screw.Full size tableFigure 7Automatic feeder efficiency at different feed pellet sizes and rotational speeds of screw.Full size imageFigure 8Automatic feeder efficiency at different feed pellet sizes and air flow rates.Full size imageFigure 9Automatic feeder efficiency at different rotational speeds of screw and air flow rates.Full size imageIt could be noticed that increasing the feed pellets size from 1 to 3 mm, tends to increase the automatic feeder efficiency from 55.79 to 69.41, 57.10 to 81.78, 72.48 to 86.13, 69.96 to 86.81 and 71.19 to 86.58% at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively. The results also indicate that the automatic feeder efficiency increased from 55.79 to 71.19, 57.98 to 78.29 and 69.41 to 86.58% at 1, 2 and 3 mm feed pellets sizes, respectively when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 7.The statistical analysis showed that the differences between the obtained data of automatic feeder efficiency due to the effect of feed pellets size (A) and air flow rate (B) were significant. Regarding the effect of rotational speed of screw, there were significant differences between rotational speeds of screw 1, 2 and 3, meanwhile, there were no significant differences between rotational speeds of screw 3, 4 and 5. The analysis showed also that the interaction between both ABC was non-significant. On the other hand, the interaction between the effect of both AB, AC and BC on the data was significant as shown in Table 2.Regarding the effect of feed pellet size and air flow rate on the automatic feeder productivity, the results indicate that the automatic feeder efficiency increases with increasing the feed pellets size and flow rate. It increased from 53.91 to 70.69, 65.23 to 81.19 and 76.78 to 94.54% for 10, 15 and 20 m3 min−1 air flow rate, respectively, when the feed pellets size increased from 1 to 3 mm. The results also indicate that the automatic feeder efficiency increased from 53.91 to 76.78, 63.14 to 83.89 and 70.69 to 94.54% at 1, 2 and 3 mm feed pellets size, respectively, when the air flow rate increased from 10 to 20 m3 min−1 as shown in Fig. 8.The results also indicate that the automatic feeder efficiency increased from 41.37 to 83.28, 58.53 to 81.54, 71.85 to 84.96, 69.88 to 85.59 and 71.27 to 85.98% at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively, when the air flow rate increased from 10 to 20 m3 min−1. The results also indicate that the automatic feeder efficiency increased from 41.37 to 71.27, 58.53 to 80.82 and 83.28 to 85.98% at 10, 15 and 20 m3 min−1 air flow rate, respectively, when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 9.Increasing the parameters seams to increase the productivity but regarding the efficiency, results show that the efficiency increases with increasing this parameter at (540 rpm) started to be constant and 720–900 rpm decreased in all treatments under study (Figs. 7, 9). It is concluded that efficiency with the parameters increased, became constant and decreased.Multiple regression analysis was carried out to obtain a relationship between the automatic feeder efficiency as dependent variable and different of feed pellets size, air flow rate and rotational speed of screw as independent variables. The best fit for this relationship is presented in the following equation:-$$ eta = 9.566 + 8.417PS + 2.249FR + 0.025RS{text{ R}}^{{2}} = 0.89{ ,} $$
(14)
where this equation could be applied in the range of 1 to 3 mm feed pellets size, 10 to 20 m3 min−1 air flow rate and from 180 to 900 rpm of rotational speed of screw.Specific energy consumptionTable 3, Figs. 10, 11 and 12 show the specific energy consumption of automatic feeder as affected by the different feed pellets sizes (1, 2 and 3 mm), air flow rates (10, 15 and 20 m3 min−1) and rotational speeds of screw (180, 360, 540, 720 and 900 rpm). The results indicate that the specific energy consumption of automatic feeder decreases with increasing feed pellets size, air flow rate and rotational speed of screw. It indicates that when the feed pellets size increased from 1 to 3 mm, the specific energy consumption of automatic feeder significantly decreased from 8.93 to 6.74 (by 24.52%) W h kg−1. It also indicates that when the air flow rate increased from 10 to 20 m3 min−1, the specific energy consumption of automatic feeder significantly decreased from 10.83 to 5.42 (by 49.95%) W h kg−1, while the specific energy consumption significantly decreased from 9.08 to 6.55 (by 27.86%) W h kg−1 when the rotational speed of screw increased from 180 to 900 rpm.Table 3 Specific energy consumption at different feed pellets sizes, air flow rates and rotational speeds of screw.Full size tableFigure 10Specific energy consumption at different feed pellet sizes and rotational speeds of screw.Full size imageFigure 11Specific energy consumption at different feed pellet sizes and air flow rates.Full size imageFigure 12Specific energy consumption at different rotational speeds of screw and air flow rates.Full size imageIt could be noticed that increasing the feed pellets size from 1 to 3 mm, tends to decrease the specific energy consumption from 9.87 to 7.94, 9.18 to 7.63, 9.14 to 7.30, 8.65 to 6.63 and 7.79 to 4.20 W h kg−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively. The results also indicate that the specific energy consumption decreased from 9.87 to 7.79, 9.42 to 7.65 and 7.94 to 4.20 W h kg−1 at 1, 2 and 3 mm feed pellets sizes, respectively when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 10.From statistical analysis, there were no significant differences between feed pellets sizes 1 and 2 on the specific energy consumption, meanwhile, there were significant differences between feed pellets size 3 and 1 and 2 on the specific energy consumption. Regarding the effect of air flow rate, there were significant differences between air flow rates and specific energy consumption. Regarding the effect of rotational speed of screw, there were significant differences between rotational speeds of screw 1, 2, 4 and 5 on the specific energy consumption, meanwhile, there were no significant differences between rotational speeds of screw 2 and 3 on the specific energy consumption. The analysis showed also that the interaction between both ABC was non-significant. On the other hand, the interaction between the effect of both AB, AC and BC on the data was significant as shown in Table 3.Regarding the effect of feed pellet size and air flow rate on the specific energy consumption, the results indicate that the specific energy consumption decreases with increasing the feed pellets size and flow rate. It decreased from 12.05 to 9.07, 8.81 to 6.56 and 5.92 to 4.59 W h kg−1 for 10, 15 and 20 m3 min−1 air flow rate, respectively, when the feed pellets size increased from 1 to 3 mm. The results also indicate that the specific energy consumption decreased 12.05 to 5.92, 11.37 to 5.75 and 9.07 to 4.59 W h kg−1 at 1, 2 and 3 mm feed pellets size, respectively, when the air flow rate increased from 10 to 20 m3 min−1 as shown in Fig. 11.The results also indicate that the specific energy consumption decreased from 12.31 to 6.18, 11.43 to 5.63, 11.21 to 5.63, 10.38 to 5.21 and 8.81 to 4.46 W h kg−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively, when the air flow rate increased from 10 to 20 m3 min−1. The results also indicate that the specific energy consumption decreased from 12.31 to 8.81, 8.75 to 6.37 and 6.18 to 4.46 W h kg−1 at 10, 15 and 20 m3 min−1 air flow rate, respectively, when the rotational speed of screw increased from 180 to 900 rpm as shown Fig. 12.
Multiple regression analysis was carried out to obtain a relationship between the specific energy consumption of automatic feeder as dependent variable and different of feed pellets size, air flow rate and rotational speed of screw as independent variables. The best fit for this relationship is presented in the following equation:-$$ SEC = 20.045 – 1.095PS – 0.541FR – 0.003RS{text{ R}}^{{2}} = 0.92 , {.} $$
(15)
This equation could be applied in the range of 1 to 3 mm feed pellets size, 10 to 20 m3 min−1 air flow rate and from 180 to 900 rpm of rotational speed of screw.Total costs of automatic feederTable 4, Figs. 13, 14 and 15 show the total cost of automatic feeder as affected by the different feed pellets sizes (1, 2 and 3 mm), air flow rates (10, 15 and 20 m3 min−1) and rotational speeds of screw (180, 360, 540, 720 and 900 rpm). The results indicate that the total cost of automatic feeder decreases with increasing feed pellets size, flow rate and rotational speed of screw. It indicates that when the feed pellets size increased from 1 to 3 mm, the total cost of automatic feeder significantly decreased from 0.15 to 0.11 (by 26.27%) EGP kg−1. It also indicates that when the air flow rate increased from 10 to 20 m3 min−1, the total cost of automatic feeder significantly decreased from 0.16 to 0.09 (by 43.75%) EGP kg−1, while the total cost of automatic feeder significantly decreased from 0.16 to 0.10 (by 37.50%) EGP kg−1 when the rotational speed of screw increased from 180 to 900 rpm.Table 4 Total cost of automatic feeder at different feed pellets sizes, air flow rate and rotational speeds of screw.Full size tableFigure 13Total cost of automatic feeder at different feed pellet sizes and rotational speeds of screw.Full size imageFigure 14Total cost of automatic feeder at different feed pellet sizes and air flow rates.Full size imageFigure 15Total cost of automatic feeder at different rotational speeds of screw and air flow rate.Full size imageIt could be noticed that increasing the feed pellets size from 1 to 3 mm, tends to decrease the total cost of automatic feeder from 0.18 to 0.14, 0.16 to 0.12, 0.15 to 0.11, 0.13 to 0.09 and 0.12 to 0.08 EGP kg−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively. The results also indicate that the total cost of automatic feeder decreased from 0.18 to 0.12, 0.17 to 0.10 and 0.14 to 0.08 EGP kg−1 at 1, 2 and 3 mm feed pellets sizes, respectively when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 13.From statistical analysis, there were no significant differences between feed pellets sizes 1 and 2 on the total cost of automatic feeder, meanwhile, there were significant differences between feed pellets size 3 and 1 and 2 on the total cost of automatic feeder. Regarding the effect of air flow rate, there were significant differences between air flow rates and specific energy consumption. Regarding the effect of rotational speed of screw, there were no significant differences between rotational speeds of screw 1 and 2, also 3 and 4 on the total cost of automatic feeder, meanwhile, there were significant differences between rotational speeds of screw 2 and 3 on the total cost of automatic feeder.Regarding the effect of feed pellet size and flow rate on the total cost of automatic feeder, the results indicate that the total cost of automatic feeder decreases with increasing the feed pellets size and air flow rate. It decreased from 0.18 to 0.13, 0.16 to 0.11 and 0.10 to 0.08 EGP kg−1 for 10, 15 and 20 m3 min−1 air flow rate, respectively, when the feed pellets size increased from 1 to 3 mm. The results also indicate that the total cost of automatic feeder decreased from 0.18 to 0.10, 0.16 to 0.10 and 0.13 to 0.08 EGP kg−1 at 1, 2 and 3 mm feed pellets size, respectively, when the air flow rate increased from 10 to 20 m3 min−1 as shown in Fig. 14.The results also indicate that the total cost of automatic feeder decreased from 0.22 to 0.11, 0.18 to 0.10, 0.16 to 0.10, 0.13 to 0.09 and 0.12 to 0.07 EGP kg−1 at 180, 360, 540, 720 and 900 rpm rotational speed of screw, respectively, when the air flow rate increased from 10 to 20 m3 min−1. The results also indicate that the total cost of automatic feeder decreased from 0.22 to 0.12, 0.16 to 0.11 and 0.11 to 0.07 EGP kg−1 for 10, 15 and 20 m3 min−1 air flow rate, respectively, when the rotational speed of screw increased from 180 to 900 rpm as shown in Fig. 15.
Multiple regression analysis was carried out to obtain a relationship between the total costs of automatic feeder as dependent variable and different of feed pellets size, air flow rate and rotational speed of screw as independent variables. The best fit for this relationship is presented in the following equation:$$ TC = 0.315 – 0.020PS – 0.006FR – 8.8 times 10^{ – 5} RS{text{ R}}^{{2}} = 0.87{,} $$
(16)
where: TC is the total cost of automatic feeder, EGP kg−1.This equation could be applied in the range of 1 to 3 mm feed pellets size, 10 to 20 m3 min−1 air flow rate and from 180 to 900 rpm of rotational speed of screw. More

FLUXNET dataThe data used in this study belong to the FLUXNET LaThuile9 and FLUXNET2015 Tier 1 and Tier 2 datasets10, which make up the global network of CO2, water vapour and energy flux measurements. We merged the two FLUXNET releases and retained the FLUXNET2015 (the most recent and with a robust quality check) version of the data when the site was present in both datasets. Croplands were removed to avoid the inclusion of sites that are heavily managed in the analysis (for example, fertilization and irrigation).The sites used cover a wide variety of climate zones (from tropical to Mediterranean to Arctic) and vegetation types (wetlands, shrublands, grasslands, savanna, evergreen and deciduous forests). It should be noted though that tropical forests are underrepresented in the FLUXNET database (Extended Data Figs. 1, 3).Sites were excluded in cases in which: (i) data on precipitation or radiation were not available or completely gap-filled; (ii) the calculation of functional properties failed because of low availability of measured data (see ‘Calculation of ecosystem functions from FLUXNET’); and (iii) fluxes showed clear discontinuities in time series indicating a change of instrumentation set-up (for example, changes in the height of the ultrasonic anemometer or gas analyser).The final number of sites selected was 203 (1,484 site years). The geographical distribution is shown in Extended Data Fig. 1, the distribution in the climate space is shown in Extended Data Fig. 2 and the fraction of sites for each climate classes is reported in Extended Data Fig. 3.For each site, we downloaded the following variables at half-hourly temporal resolution: (i) gross primary productivity (GPP, μmol CO2 m–2 s–1) derived from the night-time flux partitioning26 (GPP_NT_VUT_50 in FLUXNET 2015 and GPP_f in LaThuile), (ii) net ecosystem exchange (NEE, μmol CO2 m–2 s–1) measurements filtered using annual friction velocity (u*, m s−1) threshold (NEE_VUT_50 in FLUXNET 2015; NEE in LaThuile); (iii) latent heat (LE, W m−2) fluxes, which were converted to evapotranspiration (ET, mm); (iv) sensible heat (H, W m−2) fluxes; (v) air temperature (Tair, °C); (vi) vapour pressure deficit (VPD, hPa); (vii) global shortwave incoming radiation (SWin, W m−2); viii) net radiation (Rn, W m−2); (ix) ground heat flux (G, W m−2); (x) friction velocity u* (m s−1); and (xi) wind speed (u, m s−1). For the energy fluxes (H, LE) we selected the fluxes not corrected for the energy balance closure to guarantee consistency between the two FLUXNET datasets (in the LaThuile dataset energy fluxes were not corrected).The cumulative soil water index (CSWI, mm) was computed as a measure of water availability according to a previous report27. Half-hourly values of transpiration estimates (T, mm) were calculated with the transpiration estimation algorithm (TEA)28. The TEA has been shown to perform well against both model simulations and independent sap flow data28.For 101 sites, ecosystem scale foliar N content (N%, gN 100 g−1) was computed as the community weighted average of foliar N% of the major species at the site sampled at the peak of the growing season or gathered from the literature29,30,31,32. Foliar N% for additional sites was derived from the FLUXNET Biological Ancillary Data Management (BADM) product and/or provided by site principal investigators (Supplementary Table 1, Extended Data Fig. 1). It should be noted that this compilation of N% data might suffer from uncertainties resulting from the scaling from leaves to the eddy covariance footprint, the sampling strategy (including the position along the vertical canopy profile), the species selection and the timing of sampling. About 30% of the data comes from a coordinated effort that minimized these uncertainties29,30, and for the others we collected N% data that were representative for the eddy covariance footprint31,32.Maximum leaf area index (LAImax, m2 m−2) and maximum canopy height (Hc, m) were also collected for 153 and 199 sites, respectively, from the literature32,33, the BADM product, and/or site principal investigators.Earth observation retrievals of above-ground biomass (AGB, tons of dry matter per hectare (t DM ha−1)) were extracted from the GlobBiomass dataset34 at its original resolution (grid cell 100 × 100 m) for each site location. All the grid cells in a 300 × 300 m and 500 × 500 m window around each location were selected to estimate the median and 95th percentiles of AGB for each site. The median of AGB was selected to avoid the contribution of potential outliers to the expected value of AGB. The analysis further explored the contribution of higher percentiles in the local variation of AGB as previous studies have highlighted the contribution of older and larger trees in uneven stand age plots to ecosystem functioning35. According to the evaluation against AGB measured at 71 FLUXNET sites (Extended Data Fig. 10), we decided to use the product with median AGB values extracted from the 500 × 500 m window.A total of 94 sites have all the data on vegetation structure (N%, LAImax, Hc, and AGB).The list of sites is reported in Supplementary Table 1 along with the plant functional type (PFT), Köppen-Geiger classification, coordinates, and when available N%, LAImax, Hc and AGB.In this study we did not make use of satellite information, with the exception of the AGB data product. Future studies will benefit from new missions such as the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS), the fluorescence explorer (FLEX), hyperspectral, and radar and laser detection and ranging (LiDAR) missions (for example, Global Ecosystem Dynamics Investigation (GEDI)), to characterize a multivariate space of structural and functional properties.Calculation of ecosystem functions from FLUXNETStarting from half-hourly data, we calculated at each site a single value for each of the ecosystem functions listed below. For the calculations of functional properties we used, unless otherwise indicated, good-quality data: quality flag 0 (measured data) and 1 (good-quality gap-filled data) in the FLUXNET dataset.Gross primary productivity at light saturation (GPPsat)GPP at light saturation using photosynthetically active radiation as driving radiation and 2,000 μmol m−2 s−1 as saturating light. GPPsat represents the ecosystem-scale maximum photosynthetic CO2 uptake15,30,36. The GPPsat was estimated from half-hourly data by fitting the hyperbolic light response curves with a moving window of 5 days and assigned at the centre of the moving window30,37. For each site the 90th percentile from the GPPsat estimates was then extracted.Maximum net ecosystem productivity (NEPmax)This was computed as the 90th percentile of the half-hourly net ecosystem production (NEP = −NEE) in the growing season (that is, when daily GPP is higher than 30% of the GPP amplitude). This metric represents the maximum net CO2 uptake of the ecosystem.Basal ecosystem respiration (Rb and Rbmax)Basal ecosystem respiration at reference temperature of 15 °C was derived from night-time NEE measurements26. Daily basal ecosystem respiration (Rbd) was derived by fitting an Arrhenius type equation over a five-day moving window and by keeping the sensitivity to temperature parameter (E0) fixed as in the night-time partitioning algorithms26,38. Rbd varies across seasons because it is affected by short-term variations in productivity33,39, phenology40 and water stress41. For each site, the mean of the Rbd (Rb) and the 95th percentile (Rbmax) were computed. The calculations were conducted with the REddyProc R package v.1.2.2 (ref. 38).Apparent carbon-use efficiency (aCUE)The aCUE as defined in this study is the efficiency of an ecosystem to sequester the carbon assimilated with photosynthesis39. aCUE is an indication of the proportion of respired carbon with respect to assimilated carbon within one season. A previous report6 showed that little of the variability in aCUE can be explained by climate or conventional site characteristics, and suggested an underlying control by plant, faunal and microbial traits, in addition to site disturbance history. Daily aCUE (aCUEd) is defined as aCUEd = 1 − (Rbd/GPPd), where GPPd is daily mean GPP and Rbd is derived as described above. For each site, aCUE was computed as the median of aCUEd.Metrics of water-use efficiency (WUE)Various metrics of WUE are described below: stomatal slope or slope coefficient (G1), underlying water-use efficiency (uWUE), and water-use efficiency based on transpiration (WUEt). The three metrics were used because they are complementary, as shown in previous studies11,42.Stomatal slope or slope coefficient (G1)This is the marginal carbon cost of water to the plant carbon uptake. G1 is the key parameter of the optimal stomatal model derived previously43. G1 is inversely related to leaf-level WUE. At leaf level, G1 is calculated using nonlinear regression and can be interpreted as the slope between stomatal conductance and net CO2 assimilation, normalized for VPD and CO2 concentration43. A previous report42 showed the potential of the use of G1 at ecosystem scale, where stomatal conductance is replaced by surface conductance (Gs), and net assimilation by GPP. The methodology is implemented in the bigleaf R package44. The metric was computed in the following situations: (i) incoming shortwave radiation (SWin) greater than 200 W m−2; (ii) no precipitation event for the last 24 h45, when precipitation data are available; and (iii) during the growing season: daily GPP > 30% of its seasonal amplitude44.Underlying water-use efficiency (uWUE)The underlying WUE was computed following a previous method46. uWUE is a metric of water-use efficiency that is negatively correlated to G1 at canopy scale44:$${rm{uWUE}}=frac{{rm{GPP}}sqrt{{rm{VPD}}}}{{rm{ET}}}.$$uWUE was calculated using the same filtering that was applied for the calculation of G1. The median of the half-hourly retained uWUE values was computed for each site and used as a functional property.Water-use efficiency based on transpiration (WUEt)The WUE based on transpiration (T) was computed to reduce the confounding effect resulting from soil evaporation11,28:$${{rm{WUE}}}_{{rm{t}}}=frac{{rm{GPP}}}{T},$$where T is the mean annual transpiration calculated with the transpiration estimation algorithm (TEA) developed by in a previous study28 and GPP is the mean annual GPP.Maximum surface conductance (G
smax)Surface conductance (Gs) was computed by inverting the Penman–Monteith equation after calculating the aerodynamic conductance (Ga).Among the different formulations of Ga (m s–1) in the literature, we chose to use here the calculation of the canopy (quasi-laminar) boundary layer conductance to heat transfer, which ranges from empirical to physically based (for example, ref. 47). Other studies48,49 suggested an empirical relationship between Ga, the horizontal wind speed (u) and the friction velocity, u*:$${G}_{{rm{a}}}=frac{1}{(frac{u}{{u}^{* 2}}+6.2u{* }^{-0.67})}$$Gs (m s−1) is computed by inverting the Penman–Monteith equation:$${G}_{{rm{s}}}=frac{{{rm{LEG}}}_{{rm{a}}}gamma }{Delta ({R}_{{rm{n}}}-G-S)+rho {C}_{{rm{p}}}{G}_{{rm{a}}}{rm{VPD}}-{rm{LE}}(Delta +gamma )}$$where Δ is the slope of the saturation vapour pressure curve (kPa K−1), ρ is the air density (kg m−3), Cp is the specific heat of the air (J K−1 kg−1), γ is the psychrometric constant (kPa K−1), VPD (kPa), Rn (W m−2), G (W m−2) and S is the sum of all energy storage fluxes (W m−2) and set to 0 as not available in the dataset. When not available, G also was set to 0.Gs represents the combined conductance of the vegetation and the soil to water vapour transfer. To retain the values with a clear physiological interpretation, we filtered the data as we did for the calculation of G1.For each site, the 90th percentile of the half-hourly Gs was calculated and retained as the maximum surface conductance of each site (Gsmax). Gs was computed using the bigleaf R package44.Maximum evapotranspiration in the growing season (ETmax)This metric represents the maximum evapotranspiration computed as the 95th percentile of ET in the growing season and using the data retained after the same filtering applied for the G1 calculation.Evaporative fraction (EF)EF is the ratio between LE and the available energy, here calculated as the sum of H + LE (ref. 50). For the calculation of EF, we used the same filtering strategy as for G1. We first calculated mean daytime EF. We then computed  the EF per site as the growing season average of daytime EF. We also computed the amplitude of the EF in the growing season by calculating the interquartile distance of the distribution of mean daytime EF (EFampl).Principal component analysisA PCA was conducted on the multivariate space of the ecosystem functions. Each variable (ecosystem functional property, EFP) was standardized using z-transformation (that is, by subtracting its mean value and then dividing by its standard deviation). From the PCA results we extracted the explained variance of each component and the loadings of the EFPs, indicating the contribution of each variable to the component. We performed the PCA using the function PCA() implemented in the R package FactoMineR51.We justify using PCA over nonlinear methods because it is an exploratory technique that is highly suited to the analysis of the data volume used in this study, whereas other nonlinear methods applied to such data would be over-parameterized. For the same reason, PCA was used in previous work concerning the global spectrum of leaf and plant traits, and fluxes1,3,52.To test the significance of dimensionality of the PCA, we used a previously described methodology53. We used the R package ade4 (ref. 54) and evaluated the number of significant components of the PCA to be retained to minimize both redundancy and loss of information (Supplementary Information 2). We tested the significance of the PCA loadings using a combination of the bootstrapped eigenvector method55 and a threshold selected using the number of dimensions56 (Supplementary Information 2).Predictive variable importanceA random forests (RF) analysis was used to identify the vegetation structure and climate variables that contribute the most to the variability of the significant principal components, which were identified with the PCA analysis (see ‘Principal component analysis’). In the main text we refer to the results of this analysis as ‘predictive variable importance’ to distinguish this to the ‘causal variable importance’ described below.The analysis was conducted using the following predictor variables: as structural variables, N% (gN 100 g−1), LAImax (m2 m−2), AGB (t DM ha−1) and Hc (m); as climatic variables, mean annual precipitation (P, mm), mean VPD during the growing season (VPD, hPa), mean shortwave radiation (SWin, W m−2), mean air temperature (Tair, °C); and the cumulative soil water index (CSWI, −), as indicator of site water availability.We used partial dependencies of variables to assess the relationship between individual predictors and the response variable (that is, PC1, PC2 and PC3).The results from the partial dependency analysis can be used to determine the effects of individual variables on the response, without the influence of the other variables. The partial dependence function was calculated using the pdp R package57.The partial dependencies were calculated restricted to the values that lie within the convex hull of their training values to reduce the risk of interpreting the partial dependence plot outside the range of the data (extrapolation).Invariant causal regression models and causal variable importanceWe have quantified the dependence of the principal components on the different structural and climatic variables using nonlinear regression. Such dependencies can only be interpreted causally if the regression models are in fact causal regression models (see Supplementary Information 3 for a formal definition), which may not be the case if there are hidden confounders. To see whether the regression models allow for a causal interpretation, we use invariant causal prediction58. This method investigates whether the regression models are stable with respect to different patterns of heterogeneity in the data, encoded by different environments (that is, subsets of the original dataset). The rationale is that a causal model, describing the full causal mechanism for the response variable, should be invariant with respect to changes in the environment if the latter does not directly influence the response variable13,59. Other non-causal models may be invariant, too, but a non-invariant model cannot be considered causal.How to choose the environments is a modelling choice that must satisfy the following criteria. First, it should be possible to assign each data point to exactly one environment. Second, the environments should induce heterogeneity in the data, so that, for example, the predictor variables have different distributions across environments. Third, the environments must not directly affect the response variable, only via predictors, although the distribution of the response may still change between environments. The third criterion can be verified by expert knowledge and is assumed to hold for our analysis. In addition, if it is violated, then, usually, no set is invariant58, which can be detected from data.In our analysis, we assigned each data point (that is, each site) to one of two environments (two subsets of the original dataset): the first includes forest sites in North America, Europe or Asia; and the second includes non-forest and forest ecosystems from South America, Africa or Oceania, and non-forest ecosystems from North America, Europe or Asia (see Supplementary Information 3.1.3.1 for details). Our choice satisfies the method’s assumption that the distribution of the predictors is different between the two environments (that is, they induce heterogeneity in the data; see Supplementary Fig. 3.1). Environments that are too small or too homogeneous do not provide any evidence against the full set of covariates being a candidate for the set of causal predictors. Other choices of environments than the one presented here yield consistent results (Supplementary Information 3.2.1, Supplementary Fig. 3.4).For each subset of predictors, we test whether the corresponding regression model is invariant (yielding the same model fit in each environment). Although many models were rejected and considered non-invariant, the full model (with all the nine predictors and used in the predictive variable importance analysis) was accepted as invariant, establishing the full set of covariates as a reasonable candidate for the set of direct causal predictors. We used both RF (randomForest package in R60) and generalized additive models, GAM61 (mgcv package62 in R) to fit the models. Both methods lead to comparable results but with a better average performance of the RF: GAM led to slightly better results than RF for PC1, whereas for PC2 and PC3 RF showed a much better model performance (Supplementary Table 3.1, Supplementary Information 3.2.2). Therefore, in the main text we showed only the results from the RF (except for PC1).If, indeed, the considered regression models are causal, this allows us to make several statements. First, we can test for the existence of causal effects by testing for statistical significance of the respective predictors in the fitted models. Second, we can use the response curves of the fitted model to define a variable importance measure with a causal interpretation. In the main text we refer to this variable importance as ‘causal variable importance’. For details, see Supplementary Information 3.1.2. More formally, we considered the expected value of the predicted variables (the principal components) under joint interventions on all covariates (AGB, Hc, LAImax, N%, Tair, VPD, SWin, CSWI and P) at once, and then, to define the importance, we quantified how this expected value depends on the different covariates. We applied the same analysis to groups of vegetation structural and climate covariates (see ‘Groupwise variable importance’ in Supplementary Information 3.1.2.3, 3.2.3).The details of the methodology and the results are described in Supplementary Information 3, in which we also provide further details on the choice of environment variable and on the statistical tests that we use to test for invariance. An overview of the invariance-based methodology is shown in Supplementary Fig. 3.1.Land surface model runsWe run two widely used land surface models: Orchidee-CN (OCN) and Jena Scheme for Biosphere Atmosphere Coupling in Hamburg (JSBACH):OCNThe dynamic global vegetation model OCN is a model of the coupled terrestrial carbon and nitrogen cycles63,64, derived from the ORCHIDEE land surface model. It operates at a half-hourly timescale and simulates diurnal net carbon, heat and water exchanges, as well as nitrogen trace gas emissions, which jointly affect the daily changes in leaf area index, foliar nitrogen, and vegetation structure and growth. The main purpose of the model is to analyse the longer-term (interannual to decadal) implication of nutrient cycling for the modelling of land–climate interactions64,65. The model can run offline, driven by observed meteorological parameters, or coupled to the global circulation model.JSBACHJSBACH v.3 is the land surface model of the MPI Earth System Model66,67. The model operates at a half-hourly time step and simulates the diurnal net exchange of momentum, heat, water and carbon with the atmosphere. Daily changes in leaf area index and leaf photosynthetic capacity are derived from a prognostic scheme assuming a PFT-specific set maximum leaf area index and a set of climate responses modulating the seasonal course of leaf area index. Carbon pools are prognostic allowing for simulating the seasonal course of net land–atmosphere carbon exchanges.We selected OCN and JSBACH because they are widely used land surface models with different structures. JSBACH is a parsimonious representation of the terrestrial energy, water and carbon exchanges used to study the coupling of land and atmosphere processes in an Earth system model67. OCN has also been derived from the land surface model ORCHIDEE68, but it includes a more comprehensive representation of plant physiology, including a detailed representation of the tight coupling of the C and N cycling63. Both models contribute to the annual global carbon budget of the Global Carbon Project69 and have shown good performance compared to a number of global benchmarks. OCN was further used in several model syntheses focused on the interaction between changing N deposition and CO2 fertilization70,71,72. Both OCN and JSBACH can operate at a half-hourly timescale and simulate net and gross carbon exchanges, water and energy fluxes, and therefore are ideal for the extraction of ecosystem functional properties, as done with the eddy covariance data.The models were driven by half-hourly meteorological variables (shortwave and longwave downward flux, air temperature and humidity, precipitation, wind speed and atmospheric CO2 concentrations) observed at the eddy covariance sites. OCN was furthermore driven by N deposition fields73. Vegetation type, soil texture and plant available water were prescribed on the basis of site observations, but no additional site-specific parameterization was used. Both models were brought into equilibrium with respect to their ecosystem water storage and biogeochemical pools by repeatedly looping over the available site years. We added random noise (mean equal to 0 and standard deviation of 5% of the flux value) to the fluxes simulated by the models to mimic the random noise of the eddy covariance flux observations. An additional test conducted without noise addition showed only a marginal effect on the calculations of the functional properties and the ecosystem functional space.We used runs of the JSBACH and OCN model for 48 FLUXNET sites (Supplementary Table 1). The simulated fluxes were evaluated against the observation to assess the performance of the models at the selected sites. From the model outputs and from each site we derived the ecosystem functions using the same methodology described above. Then the PCA analysis was performed on the three datasets (FLUXNET, OCN and JSBACH) and restricted to the 48 sites used to run the models. We ran the models only on the subset of sites for which the information for the parameterization and high-quality forcing was available. However, the different ecosystem functions emerge from the model structure and climatological conditions. Therefore, even with a smaller set of site we can evaluate whether models reproduce the key dimensions of terrestrial ecosystem function by comparing the PCA results from FLUXNET and the model runs.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this paper. More

1.Chesson, P. Mechanisms of maintenance of species diversity. Annu. Rev. Ecol. Syst. 31, 343–366 (2000).Article 

Google Scholar 
2.Wright, J. S. Plant diversity in tropical forests: A review of mechanisms of species coexistence. Oecologia 130, 1–14 (2002).ADS 
PubMed 
Article 

Google Scholar 
3.Janzen, D. H. Herbivores and the number of tree species in tropical forests. Am. Nat. 104, 501–528 (1970).Article 

Google Scholar 
4.Connell, J. On the role of natural enemies in preventing competitive exclusion in some marine animals and rain forest trees. Dyn. Popul. 298, 312 (1971).
Google Scholar 
5.Terborgh, J. W. Toward a trophic theory of species diversity. PNAS 112, 11415–11422 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
6.Johnson, D. J., Beaulieu, W. T., Bever, J. D. & Clay, K. Conspecific negative density dependence and forest diversity. Science 336, 904–907 (2012).ADS 
CAS 
PubMed 
Article 

Google Scholar 
7.LaManna, J. A. et al. Plant diversity increases with the strength of negative density dependence at the global scale. Science 356, 1389–1392 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
8.Chisholm, R. A. & Muller-Landau, H. C. A theoretical model linking interspecific variation in density dependence to species abundances. Theor. Ecol. 4, 241–253 (2011).Article 

Google Scholar 
9.Mangan, S. A. et al. Negative plant–soil feedback predicts tree-species relative abundance in a tropical forest. Nature 466, 752–755 (2010).ADS 
CAS 
PubMed 
Article 

Google Scholar 
10.Chisholm, R. A. & Fung, T. Comment on “Plant diversity increases with the strength of negative density dependence at the global scale”. Science 360, eaar4685 (2018).PubMed 
Article 
CAS 

Google Scholar 
11.Hülsmann, L. & Hartig, F. Comment on “Plant diversity increases with the strength of negative density dependence at the global scale”. Science 360, eaar2435 (2018).PubMed 
Article 
CAS 

Google Scholar 
12.Detto, M., Visser, M. D., Wright, S. J. & Pacala, S. W. Bias in the detection of negative density dependence in plant communities. Ecol. Lett. 22, 1923–1939 (2019).PubMed 
Article 

Google Scholar 
13.LaManna, J. A. et al. Response to Comment on “Plant diversity increases with the strength of negative density dependence at the global scale”. Science 360, eaar3824 (2018).PubMed 
Article 
CAS 

Google Scholar 
14.LaManna, J. A. et al. Response to Comment on “Plant diversity increases with the strength of negative density dependence at the global scale”. Science 360, eaar5245 (2018).PubMed 
Article 
CAS 

Google Scholar 
15.LaManna, J. A., Mangan, S. A. & Myers, J. A. Conspecific negative density dependence and why its study should not be abandoned. Ecosphere 12, e03322 (2021).Article 

Google Scholar 
16.Gaston, K. J. Global patterns in biodiversity. Nature 405, 220–227 (2000).CAS 
PubMed 
Article 

Google Scholar 
17.Mittelbach, G. G. et al. Evolution and the latitudinal diversity gradient: Speciation, extinction and biogeography. Ecol. Lett. 10, 315–331 (2007).PubMed 
Article 

Google Scholar 
18.Janzen, D. H. Why mountain passes are higher in the tropics. Am. Nat. 101, 233–249 (1967).Article 

Google Scholar 
19.Ricklefs, R. E. & He, F. Region effects influence local tree species diversity. PNAS 113, 674–679 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
20.Comita, L. S. et al. Testing predictions of the Janzen-Connell hypothesis: A meta-analysis of experimental evidence for distance- and density-dependent seed and seedling survival. J. Ecol. 102, 845–856 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
21.Currie, D. J. Energy and large-scale patterns of animal- and plant-species richness. Am. Nat. 137, 27–49 (1991).Article 

Google Scholar 
22.Grosso, S. D. et al. Global potential net primary production predicted from vegetation class, precipitation, and temperature. Ecology 89, 2117–2126 (2008).PubMed 
Article 

Google Scholar 
23.Chase, J. M. Stochastic Community Assembly Causes Higher Biodiversity in More Productive Environments. Science 27, (2010).24.O’Brien, E. M. Climatic gradients in woody plant species richness: Towards an explanation based on an analysis of Southern Africa’s woody flora. J. Biogeography 20, 181–198 (1993).Article 

Google Scholar 
25.McCain, C. M. & Grytnes, J.-A. Elevational Gradients in Species Richness. In eLS (American Cancer Society, 2010).26.Barry, R. G. Mountain Weather and Climate (Cambridge University Press, 2008).Book 

Google Scholar 
27.LaManna, J. A., Walton, M. L., Turner, B. L. & Myers, J. A. Negative density dependence is stronger in resource-rich environments and diversifies communities when stronger for common but not rare species. Ecol. Lett. 19, 657–667 (2016).PubMed 
Article 

Google Scholar 
28.Zhu, K., Woodall, C. W., Monteiro, J. V. D. & Clark, J. S. Prevalence and strength of density-dependent tree recruitment. Ecology 96, 2319–2327 (2015).PubMed 
Article 

Google Scholar 
29.Yao, J. et al. Abiotic niche partitioning and negative density dependence across multiple life stages in a temperate forest in northeastern China. J. Ecol. 108, 1299–1310 (2020).Article 

Google Scholar 
30.Leigh, E. G. et al. Why do some tropical forests have so many species of trees?. Biotropica 36, 447–473 (2004).
Google Scholar 
31.Terborgh, J. Enemies maintain hyperdiverse tropical forests. Am. Nat. 179, 303–314 (2012).PubMed 
Article 

Google Scholar 
32.Altman, J. et al. Linking spatiotemporal disturbance history with tree regeneration and diversity in an old-growth forest in northern Japan. PPEES 21, 1–13 (2016).
Google Scholar 
33.Kubota, Y., Hirao, T., Fujii, S., Shiono, T. & Kusumoto, B. Beta diversity of woody plants in the Japanese archipelago: The roles of geohistorical and ecological processes. J. Biogeogr. 41, 1267–1276 (2014).Article 

Google Scholar 
34.Mori, A. S. Local and biogeographic determinants and stochasticity of tree population demography. J. Ecol. 107, 1276–1287 (2019).Article 

Google Scholar 
35.Oohata, S. Distribution of tree species along the temperature gradient in the Japan archipelago (ii).: Life form and species distribution. Jap. J. Ecol. 40, 71–84 (1990).ADS 

Google Scholar 
36.Kira, T. A Climatological Interpretation of Japanese Vegetation Zones 21–30 (Springer, 1977).
Google Scholar 
37.Mori, A. S. Environmental controls on the causes and functional consequences of tree species diversity. J. Ecol. 106, 113–125 (2018).Article 

Google Scholar 
38.Suzuki, S. N., Ishihara, M. I. & Hidaka, A. Regional-scale directional changes in abundance of tree species along a temperature gradient in Japan. Glob. Chan. Biol. 21, 3436–3444 (2015).ADS 
Article 

Google Scholar 
39.Hara, M. Analysis of seedling banks of a climax beech forest: Ecological importance of seedling sprouts. Vegetatio 71, 67–74 (1987).
Google Scholar 
40.Homma, K. Effects of snow pressure on growth form and life history of tree species in Japanese beech forest. J. Veg. Sci. 8, 781–788 (1997).Article 

Google Scholar 
41.Gansert, D. Treelines of the Japanese Alps—altitudinal distribution and species composition under contrasting winter climates. Flora 199, 143–156 (2004).Article 

Google Scholar 
42.Hukusima, T. et al. New phytosociological classification of beech forests in Japan. Jpn. J. Ecol. 45, 79–98 (1995).
Google Scholar 
43.Matsui, T. et al. Probability distributions, vulnerability and sensitivity in Fagus crenata forests following predicted climate changes in Japan. J. Veg. Sci. 15, 605–614 (2004).Article 

Google Scholar 
44.Johnson, D. J., Condit, R., Hubbell, S. P. & Comita, L. S. Abiotic niche partitioning and negative density dependence drive tree seedling survival in a tropical forest. Proc. R. Soc. B 284, 20172210 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
45.Ishihara, M. I. et al. Forest stand structure, composition, and dynamics in 34 sites over Japan. Ecol. Res. 26, 1007–1008 (2011).Article 

Google Scholar 
46.Law, R. et al. Ecological information from spatial patterns of plants: Insights from point process theory. J. Ecol. 97, 616–628 (2009).Article 

Google Scholar 
47.Wright, S. J. et al. Reproductive size thresholds in tropical trees: Variation among individuals, species and forests. J. Trop. Ecol. 21, 307–315 (2005).Article 

Google Scholar 
48.Zhu, Y., Comita, L. S., Hubbell, S. P. & Ma, K. Conspecific and phylogenetic density-dependent survival differs across life stages in a tropical forest. J. Ecol. 103, 957–966 (2015).Article 

Google Scholar 
49.Ripley, B. D. Spatial point pattern analysis in ecology. In Develoments in Numerical Ecology (eds Legendre, P. & Legendre, L.) 407–429 (Springer, 1987).Chapter 

Google Scholar 
50.Wiegand, T. & Moloney, K. A. Handbook of Spatial Point-Pattern Analysis in Ecology (CRC Press, 2013).Book 

Google Scholar 
51.Loosmore, N. B. & Ford, E. D. Statistical inference using the G or K point pattern spatial statistics. Ecology 87, 1925–1931 (2006).PubMed 
Article 

Google Scholar 
52.R Core Team. R: A Language and Environment for Statistical Computing (2020).53.Baddeley, A. & Turner, R. spatstat: An R Package for Analyzing Spatial Point Patterns. J. Stat. Soft. 12, 1–42 (2005).Article 

Google Scholar 
54.Wills, C., Condit, R., Foster, R. B. & Hubbell, S. P. Strong density- and diversity-related effects help to maintain tree species diversity in a neotropical forest. PNAS 94, 1252–1257 (1997).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
55.Givnish, T. J. On the causes of gradients in tropical tree diversity. J. Ecol. 87, 193–210 (1999).Article 

Google Scholar 
56.Fibich, P., Vítová, A. & Lepš, J. Interaction between habitat limitation and dispersal limitation is modulated by species life history and external conditions: A stochastic matrix model approach. Comm. Ecol. 19, 9–20 (2018).Article 

Google Scholar 
57.Miyawaki, A. A vegetation ecological view of the Japanese archipelago. Bull. Inst. Environ. Sci. Technol. Yokohama Natl. Univ. 11, 85–101 (1984).
Google Scholar 
58.Mori, A. S. et al. Community assembly processes shape an altitudinal gradient of forest biodiversity. Glo. Ecol. Biogeogr. 22, 878–888 (2013).Article 

Google Scholar 
59.Grime, J. P. Plant Strategies, Vegetation Processes, and Ecosystem Properties (Wiley, 2001).
Google Scholar 
60.Brown, C., Law, R., Illian, J. B. & Burslem, D. F. R. P. Linking ecological processes with spatial and non-spatial patterns in plant communities. J. Ecol. 99, 1402–1414 (2011).Article 

Google Scholar 
61.Bastias, C. C. et al. Species richness influences the spatial distribution of trees in European forests. Oikos 129, 380–390 (2020).Article 

Google Scholar 
62.Hülsmann, L., Chisholm, R. A. & Hartig, F. Is variation in conspecific negative density dependence driving tree diversity patterns at large scales?. Trends Ecol. Evol. 36, 151–163 (2021).PubMed 
Article 

Google Scholar 
63.Damgaard, C. & Weiner, J. It’s about time: A critique of macroecological inferences concerning plant competition. Trends Ecol. Evol. 32, 86–87 (2017).PubMed 
Article 

Google Scholar 
64.Murata, I. et al. Effects of sika deer (Cervus nippon) and dwarf bamboo (Sasamorpha borealis) on seedling emergence and survival in cool-temperate mixed forests in the Kyushu Mountains. J. For. Res. 14, 296–301 (2009).Article 

Google Scholar 
65.Ackerly, D. D. et al. The geography of climate change: Implications for conservation biogeography. Divers. Distrib. 16, 476–487 (2010).Article 

Google Scholar  More