Philippa Kaur

Biskaborn, B. K. et al. Permafrost is warming at a global scale. Nat. Commun. 10, 264 (2019).Article 
ADS 

Google Scholar 
Murton, J. B. What and where are periglacial landscapes? Permaf. Periglac. Process. 32, 186–212 (2021).Article 

Google Scholar 
Woodcroft, B. J. et al. Genome-centric view of carbon processing in thawing permafrost. Nature 560, 49–54 (2018).Article 
ADS 
CAS 

Google Scholar 
Reyes, F. & Lougheed, V. L. Rapid nutrient release from permafrost thaw in Arctic aquatic ecosystems. Arct. Antarct. Alp. Res. 47, 35–48 (2015).Article 

Google Scholar 
Fouché, J., Christiansen, C. T., Lafrenière, M. J., Grogan, P. & Lamoureux, S. F. Canadian permafrost stores large pools of ammonium and optically distinct dissolved organic matter. Nat. Commun. 11, 4500 (2020).Article 
ADS 

Google Scholar 
Alley, N. F., Hore, S. B. & Frakes, L. A. Glaciations at high-latitude Southern Australia during the Early Cretaceous. Aust. J. Earth Sci. 67, 1045–1095 (2020).Article 
ADS 

Google Scholar 
Hore, S. B., Hill, S. M. & Alley, N. F. Early Cretaceous glacial environment and paleosurface evolution within the Mount Painter Inlier, northern Flinders Ranges, South Australia. Aust. J. Earth Sci. 67, 1117–1160 (2020).Article 
ADS 
CAS 

Google Scholar 
Rodríguez-López, J. P. et al. Glacial dropstones in the western Tethys during the late Aptian–early Albian cold snap: Palaeoclimate and palaeogeographic implications for the mid-Cretaceous. Palaeogeogr. Palaeoclimatol. Palaeoecol. 452, 11–27 (2016).Article 

Google Scholar 
Schneider, S. et al. Macrofauna and biostratigraphy of the Rollrock Section, northern Ellesmere Island, Canadian Arctic Islands e a comprehensive high latitude archive of the Jurassic–Cretaceous transition. Cret. Res. 114, 104508 (2020).Article 

Google Scholar 
Jeans, C. V. & Platten, I. M. The erratic rocks of the Upper Cretaceous Chalk of England: how did they get there, ice transport or other means? Acta Geol. Pol. 71, 287–304 (2021).
Google Scholar 
Wu, C. & Rodríguez-López, J. P. Cryospheric processes in Quaternary and Cretaceous hyper-arid oases. Sedimentology 68, 755–770 (2021).Article 

Google Scholar 
Grasby, S. E., McCune, G. E., Beauchamp, B. & Galloway, J. M. Lower Cretaceous cold snaps led to widespread glendonite occurrences in the Sverdrup Basin, Canadian High Arctic. GSA Bull. 129, 771–787 (2017).Article 
CAS 

Google Scholar 
Galloway, J. M. et al. Finding the VOICE: organic carbon isotope chemostratigraphy of the Late Jurassic–Early Cretaceous of Arctic Canada. Geol. Mag. 1–15 https://doi.org/10.1017/S0016756819001316 (2019).Rogov, M. et al. Database of global glendonite and ikaite records throughout the Phanerozoic. Earth Syst. Sci. Data 13, 343–356 (2021).Article 
ADS 

Google Scholar 
Price, G. D. The evidence and implications of polar ice during the Mesozoic. Earth–Sci. Rev. 48, 183–210 (1999).Article 
ADS 

Google Scholar 
Savidge, R. A. Evidence of early glaciation of southeastern Beringia. Can. J. Earth Sci. 57, 199–226 (2020).Article 
ADS 

Google Scholar 
Wang, Y. et al. Relict sand wedges suggest a high altitude and cold temperature during the Early Cretaceous in the Ordos Basin, North China. Int. Geol. Rev. https://doi.org/10.1080/00206814.2022.2081938 (2022).Nelson, D. A., Cottle, J. M., Bindeman, I. N. & Camacho, A. Ultra-depleted hydrogen isotopes in hydrated glass record Late Cretaceous glaciation in Antarctica. Nat. Commun. 13, 5209 (2022).Article 
ADS 
CAS 

Google Scholar 
Yang, W.-B. et al. Isotopic evidence for continental ice sheet in mid-latitude region in the supergreenhouse Early Cretaceous. Sci. Rep. 3, 2732 (2013).Article 

Google Scholar 
Gao, T. et al. Accelerating permafrost collapse on the eastern Tibetan Plateau. Environ. Res. Lett. 16, 054023 (2021).Article 
ADS 

Google Scholar 
Huang, Y. B. The origin and evolution of the desert in southern Ordos in early Cretaceous: Constraint from Magnetostratigraphy of Zhidan Group and magnetic susceptibility of its sediment. Doctoral Dissertation. Lanzhou University (2010).Ma, J. Sedimentary Basin Analysis of the Cretaceous Ancient Desert in the Ordos Basin. Master’s thesis, China University of Geosciences (2020).Wu, C. H., Rodríguez-López, J. P. & Santosh, M. Plateau archives of lithosphere dynamics, cryosphere and paleoclimate: the formation of Cretaceous desert basins in east Asia. Geosci. Front. 13, 101454 (2022).Article 
CAS 

Google Scholar 
Zhu, R. X., Chen, L., Wu, F. Y. & Liu, J. L. Timing, scale and mechanism of the destruction of the North China Craton. Sci. China Earth Sci. 54, 789–797 (2011).Article 
ADS 
CAS 

Google Scholar 
Rodríguez-López, J. P., Clemmensen, L. B., Lancaster, N., Mountney, N. P. & Veiga, G. D. Archean to Recent aeolian sand systems and their preserved successions: current understanding and way forward. Sedimentology 61, 1487–1534 (2014).Article 

Google Scholar 
Murton, J. B. in Encyclopedia of Quaternary Science Vol. 3 (eds Elias, S. A. & Mock, C. J.) 436–451 (Elsevier, Amsterdam, 2013).Rodríguez-López, J. P., Van Vliet-Lanöe, B., López-Martínez, J. & Martín-García, R. Scouring by rafted ice and cryogenic pattern ground preserved in a Palaeoproterozoic equatorial proglacial lagoon succession, eastern India, Nuna supercontinent. Mar. Pet. Geol. 123, 104766 (2021).Article 

Google Scholar 
Murton, J. B., Worsley, P. & Gozdzik, J. Sand veins and wedges in cold aeolian environments. Quat. Sci. Rev. 19, 899–922 (2000).Article 
ADS 

Google Scholar 
Kovács, J., Fábián, S. A., Schweitzer, F. & Varga, G. A relict sand-wedge polygon site in north-central Hungary. Permafr. Periglac. Process. 18, 379–384 (2007).Article 

Google Scholar 
Fábián, S. Á. et al. Distribution of relict permafrost features in the Pannonian Basin, Hungary. Boreas 43, 722–732 (2014).Article 

Google Scholar 
Williams, G. E. Proterozoic (pre-Ediacaran) glaciation and the high obliquity, low-latitude ice, strong seasonality (HOLIST) hypothesis: principles and tests. Earth–Sci. Rev. 87, 61–93 (2008).Article 
ADS 

Google Scholar 
Williams, G. E., Schmidt, P. W. & Young, G. M. Strongly seasonal Proterozoic glacial climate in low palaeolatitudes: radically different climate system on the pre-Ediacaran Earth. Geosci. Front. 7, 555–571 (2016).Article 

Google Scholar 
Van Vliet-Lanoë, B. Deformations in the active layer related with ice/soil wedge growth and decay in present day Arctic. Paleoclimate implications. Ann. Soc. Géol. Nord. 13, 81–95 (2005).
Google Scholar 
Remillard, A. M. et al. Chronology and palaeoenvironmental implications of the ice-wedge pseudomorphs and composite wedge casts on the Magdalen Islands (eastern Canada). Boreas 44, 658–675 (2015).Article 

Google Scholar 
Murton, J. B. Thermokarst sediments and sedimentary structures, Tuktoyaktuk Coastlands, western Arctic Canada. Glob. Planet. Change 28, 175–192 (2001).Article 
ADS 

Google Scholar 
Harris, C., Murton, J. B. & Davies, M. C. R. An analysis of mechanisms of ice-wedge casting based on geotechnical centrifuge modelling. Geomorphology 71, 328–343 (2005).Article 
ADS 

Google Scholar 
Houmark-Nielsen, M. et al. Early and Middle Valdaian glaciations, ice-dammed lakes and periglacial interstadials in northwest Russia: new evidence from the Pyoza River area. Glob. Planet. Change 31, 215–237 (2001).Article 
ADS 

Google Scholar 
Murton, J. B. & Kolstrup, E. Ice-wedge casts as indicators of palaeotemperatures: precise proxy or wishful thinking? Prog. Phys. Geogr. 27, 155–170 (2003).Article 

Google Scholar 
Harry, D. G. & Gozdzik, J. S. Ice wedges: growth, thaw transformation, and palaeoenvironmental significance. J. Quat. Sci. 3, 39–55 (1988).Article 

Google Scholar 
Wolfe, S. A., Morse, P. D., Neudorf, C. M., Kokelj, S. V., Lian, O. B. & O’Neill, H. B. Contemporary sand wedge development in seasonally frozen ground and paleoenvironmental implications. Geomorphology 308, 215–229 (2018).Article 
ADS 

Google Scholar 
Murton, J. B. & Bateman, M. D. Syngenetic sand veins and anti-syngenetic sand wedges, Tuktoyaktuk Coastlands, western Arctic Canada. Permafr. Periglac. Process. 18, 33–47 (2007).Article 

Google Scholar 
Obu, J., Westermann, S., Kääb, A., & Bartsch, A. Ground Temperature Map, 2000–2016, Northern Hemisphere Permafrost (Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, PANGAEA, 2018)Obu, J. et al. Northern Hemisphere permafrost map based on TTOP modelling for 2000–2016 at 1 km2 scale. Earth–Sci. Rev. 193, 299–316 (2019).Article 
ADS 

Google Scholar 
Hock, R. et al. in IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (eds Pörtner, H.-O. et al.) 131–202 (Cambridge University Press, Cambridge, UK and New York, NY, USA, 2019).Mackay, J. R. The origin of hummocks, western arctic coast, Canada. Can. J. Earth Sci. 17, 996–1006 (1980).Article 
ADS 

Google Scholar 
Kokelj, S. V., Burn, C. R. & Tarnocai, C. The structure and dynamics of earth hummocks in the subarctic forest near Inuvik, Northwest Territories, Canada. Arct. Antarct. Alp. Res. 39, 99–109 (2007).Article 

Google Scholar 
Rodríguez-López, J. P., Meléndez, N., de Boer, P. L., Soria, A. R. & Liesa, C. L. Spatial variability of multicontrolled aeolian supersurfaces in central-erg and marine erg-margin systems. Aeolian Res. 11, 141–154 (2013).Article 
ADS 

Google Scholar 
Lunt, D. J. et al. Palaeogeographic controls on climate and proxy interpretation. Clim. Past 12, 1181–1198 (2016).Article 

Google Scholar 
Cheng, G., Bai, Y. & Sun, Y. Paleomagnetic study on the tectonic evolution of the Ordos Block, North China. Seismol. Geol. 10, 81–87 (1988).
Google Scholar 
Zheng, Z. et al. The apparent polar wander path for the North China Block since the Jurassic. Geophys. J. Int. 104, 29–40 (1991).Article 
ADS 

Google Scholar 
Malinverno, A., Hildebrandt, J., Tominaga, M. & Channell, J. E. T. M-sequence geomagnetic polarity time scale (MHTC12) that steadies global spreading rates and incorporates astrochronology constraints. J. Geophys. Res. 117, B06104 (2012).ADS 

Google Scholar 
Zachos, J. C., Shackleton, N. J., Revenaugh, J. S., Pälike, H. & Flower, B. P. Climate response to orbital forcing across the Oligocene–Miocene boundary. Science 292, 274–278 (2001).Article 
ADS 
CAS 

Google Scholar 
Li, M. et al. Astronomical tuning of the end-Permian extinction and the Early Triassic Epoch of South China and Germany. Earth Planet. Sci. Lett. 441, 10–25 (2016).Article 
ADS 
CAS 

Google Scholar 
Westall, F. The nature of fossil bacteria: a guide to the search for extraterrestial live. J. Geophys. Res. 104, 437–16,451 (1999).
Google Scholar 
Yang, H., Chen, Z.-Q. & Papineau, D. Cyanobacterial spheroids and other biosignatures from microdigitate stromatolites of Mesoproterozoic Wumishan Formation in Jixian, North China. Precambrian Res. 368, 106496 (2022).Article 
ADS 
CAS 

Google Scholar 
Kremer, B., Kazmierczak, J., Łukomska-Kowalczyk, M. & Kempe, S. Calcification and silicification: fossilization potential of cyanobacteria from stromatolites of Niuafo’ou’s caldera lakes (Tonga) and implications for the early fossil record. Astrobiology 12, 535–548 (2012).Article 
ADS 
CAS 

Google Scholar 
Astafieva M. M. et al. Fossil Bacteria and Other Microorganisms in Terrestrial Rocks and Astromaterials (Paleontological Institute Russian Academy of Science, Moscow, 2011).Rozanov, A. Y. & Zavarzin, G. A. Bacterial paleontology. Vestn. Akad. Med. Nauk 67, 241–245 (1997).
Google Scholar 
Perez-Mon, C., Stierli, B., Plötze, M. & Frey, B. Fast and persistent responses of alpine permafrost microbial communities to in situ warming. Sci. Total Environ. 807, 150–720 (2022).Article 

Google Scholar 
Rivkina, E. et al. Earth’s perennially frozen environments as a model of cryogenic planet ecosystems. Permafr. Periglac. Process. 29, 246–256 (2018).Article 

Google Scholar 
Vishnivetskaya, T. A. et al. Insights into community of photosynthetic microorganisms from permafrost. FEMS Microbiol. Ecol. 96, fiaa229 (2020).Article 
CAS 

Google Scholar 
Hultman, J. et al. Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 521, 208–212 (2015).Article 
ADS 
CAS 

Google Scholar 
Choe, Y. H. et al. Comparing rock-inhabiting microbial communities in different rock types from a high arctic polar desert. FEMS Microbiol. Ecol. 94, fiy070 (2018).ADS 

Google Scholar 
Wu, X. et al. Comparative metagenomics of the active layer and permafrost from low-carbon soil in the Canadian High Arctic. Environ. Sci. Technol. 55, 12683–12693 (2021).Article 
ADS 
CAS 

Google Scholar 
Vickers, M. L. et al. The duration and magnitude of Cretaceous cold events: evidence from the northern high latitudes. Geol. Soc. Am. Bull. 131, 1979–1994 (2019).Article 
ADS 
CAS 

Google Scholar 
Lehmann, J. in Ammonoid Palaeobiology: From Macroevolution to Palaeogeography (eds Klug, C. De Baets, K., Kruta I. & Mapes, R. H.) 403–429 (Springer, Amsterdam, 2015).Keller, M. A. & Macquaker, J. H. S. in Studies by the U.S. Geological Survey in Alaska: US Geological Survey Professional Paper 1814-B Vol. 15 (ed Dumoulin, J. A.) 1–35 (US Geological Survey, US Department of The Interior, Reston, 2015).Cavalheiro, L. et al. Impact of global cooling on Early Cretaceous high pCO2 world during the Weissert Event. Nat. Commun. 12, 5411 (2021).Article 
ADS 
CAS 

Google Scholar 
McArthur, J. M. et al. Palaeotemperatures, polar ice-volume, and isotope stratigraphy (Mg/Ca, d18O, d13C, 87Sr/86Sr): the Early Cretaceous (Berriasian, Valanginian, Hauterivian). Palaeogeogr. Palaeoclimatol. Palaeoecol. 248, 391–430 (2007).Article 

Google Scholar 
Lini, A., Weissert, H. & Erba, E. The Valanginian carbon isotope event: a first episode of greenhouse climate conditions during the Cretaceous. Terra Nova 4, 374–384 (1992).Article 
ADS 

Google Scholar 
Li, X. et al. Carbon isotope signatures of pedogenic carbonates from SE China: rapid atmospheric pCO2 changes during middle–late Early Cretaceous time. Geol. Mag. 151, 830–849 (2014).Article 
ADS 
CAS 

Google Scholar 
O’Brien, Ch. L. et al. Cretaceous sea-surface temperature evolution: constraints from TEX86 and planktonic foraminiferal oxygen isotopes. Earth–Sci. Rev. 172, 224–247 (2017).Article 
ADS 

Google Scholar 
Price, G. D. et al. A high-resolution Belemnite geochemical analysis of early Cretaceous (Valanginian–Hauterivian) environmental and climatic perturbations. Geochem. Geophys. Geosyst. 19, 3832–3843 (2018).Article 
CAS 

Google Scholar 
Turetsky, M. R. et al. Carbon release through abrupt permafrost thaw. Nat. Geosci. 13, 138–143 (2020).Article 
ADS 
CAS 

Google Scholar 
Van der Kolk, D. A., Whalen, M. T., Wartes, M. A., Newberry, R. J. & McCarthy, P. in Arctic to the Cordillera: Unlocking the Potential. American Association of Petroleum Geologists Pacific Section Meeting, May 8–11, Anchorage, AK, USA, Search and Discovery Article 90125 (American Association of Petroleum Geologists, 2011).Walter Anthony, K. M. et al. 21st-century modeled permafrost carbon emissions accelerated by abrupt thaw beneath lakes. Nat. Commun. 9, 3262 (2018).Article 
ADS 

Google Scholar 
Cheng, F. et al. Alpine permafrost could account for a quarter of thawed carbon based on Plio-Pleistocene palaeoclimate analogue. Nat. Commun. 13, 1329 (2022).Article 
ADS 
CAS 

Google Scholar 
Brouillette, M. How microbes in permafrost could trigger a massive carbon bomb. Nature 591, 360–362 (2021).Article 
ADS 
CAS 

Google Scholar 
Murton, J. B. in Climate Change, Observed Impacts on Planet Earth, 3rd edn (ed Letcher, T.) 281–326 (Elsevier, Amsterdam, 2021).Schnyder, J., Ruffell, A., Deconinck, J. F. & Baudin, F. Conjunctive use of spectral gamma-ray logs and clay mineralogy in defining late Jurassic–early Cretaceous palaeoclimate change (Dorset, UK). Palaeogeogr. Palaeoclimatol. Palaeoecol. 229, 303–320 (2006).Article 

Google Scholar 
Li, M. et al. Astrochronology of the Anisian stage (Middle Triassic) at the guandao reference section, south china. Earth Planet. Sci. Lett. 482, 591–606 (2018).Article 
ADS 
CAS 

Google Scholar 
Li, M. et al. Palaeoclimate proxies for cyclostratigraphy: comparative analysis using a Lower Triassic marine section in South China. Earth–Sci. Rev. 189, 125–146 (2019).Article 
ADS 
CAS 

Google Scholar 
Li, M., Hinnov, L. & Kump, L. Acycle: time–series analysis software for palaeoclimate research and education. Comput. Geosci. 127, 12–22 (2019).Article 
ADS 
CAS 

Google Scholar 
Laskar, J. et al. A long–term numerical solution for the insolation quantities of the Earth. Astron. Astrophys. 428, 261–285 (2004).Article 
ADS 

Google Scholar  More

Lavelle, P., Blanchart, E., Martin, A., Spain, A. V. & Martin, S. Impact of soil fauna on the properties of soils in the humid tropics. In Myths and Science of Soils of the Tropics (eds Lal, R. & Sanchez, P.) 157–185 (Soil Science Society of America, 1992).
Google Scholar 
Eisenhauer, N. The action of an animal ecosystem engineer: Identification of the main mechanisms of earthworm impacts on soil microarthropods. Pedobiologia 53, 343–352 (2010).Article 

Google Scholar 
Eisenhauer, N. & Eisenhauer, E. The “intestines of the soil”: The taxonomic and functional diversity of earthworms—A review for young ecologists. Preprint at https://doi.org/10.32942/osf.io/tfm5y (2020).Gates, G. E. Burmese earthworms, an introduction to the systematics and biology of megadrile oligochaetes with special reference to South-east Asia. Trans. Amer. Phil. Soc. 62, 1–326. https://doi.org/10.2307/1006214 (1972).Article 

Google Scholar 
Blakemore, R. J. Origin and means of dispersal of cosmopolitan Pontodrilus litoralis (Oligocaheta: Megascolecidae). Eur. J. Soil Biol. 443, S3–S8. https://doi.org/10.1016/j.ejsobi.2007.08.041 (2007).Article 

Google Scholar 
Seesamut, T., Sutcharit, C., Jirapatrasilp, P., Chanabun, R. & Panha, S. Morphological and molecular evidence reveal a new species of the earthworm genus Pontodrilus Perrier, 1874 (Clitellata, Megascolecidae) from Thailand and Peninsular Malaysia. Zootaxa 4496, 218–237. https://doi.org/10.11646/zootaxa.4496.1.18 (2018).Article 

Google Scholar 
Seesamut, T., Jirapatrasilp, P., Chanabun, R., Oba, Y. & Panha, S. Size variation and geographical distribution of the luminous earthworm Pontodrilus litoralis (Grube, 1855) (Clitellata, Megascolecidae) in Southeast Asia and Japan. Zookeys 862, 23–43. https://doi.org/10.3897/zookeys.862.35727 (2019).Article 

Google Scholar 
Seesamut, T., Jirapatrasilp, P., Sutcharit, C., Tongkerd, P. & Panha, S. Mitochondrial genetic population structure and variation of the littoral earthworm Pontodrilus longissimus Seesamut and Panha, 2018 along the coast of Thailand. Eur. J. Soil Biol. 93, 103091. https://doi.org/10.1016/j.ejsobi.2019.103091 (2019).Article 

Google Scholar 
Attrill, M. J. A testable linear model for diversity trends in estuaries. J. Anim. Ecol. 71, 262–269. https://doi.org/10.1046/j.1365-2656.2002.00593.x (2002).Article 

Google Scholar 
McLusky, D. S. & Elliott, M. The Estuarine Ecosystem: Ecology, Threats and Management 3rd edn. (Oxford University Press, 2004).Book 

Google Scholar 
Telesh, I. V. & Khlebovich, V. V. Principal processes within the estuarine salinity gradient: A review. Mar. Pollut. Bull. 61, 149–155. https://doi.org/10.1016/j.marpolbul.2010.02.008 (2010).Article 
CAS 

Google Scholar 
Owojori, O. J. & Reinecke, A. J. Effects of natural (flooding and drought) and anthropogenic (copper and salinity) stressors on the earthworm Aporrectodea caliginosa under field conditions. Appl. Soil Ecol. 44, 156–163. https://doi.org/10.1016/j.apsoil.2009.11.006 (2010).Article 

Google Scholar 
Guzyte, G., Sujetoviene, G. & Zaltauskaite, J. Effects of salinity on earthworm (Eisenia fetida). Environ. Eng. 8, 111 (2011).
Google Scholar 
Ganapati, P. N. & Subba Rao, B. V. S. S. R. Salinity tolerance of a littoral oligochaete, Pontodrilus bermudensis Beddard. Proc. Ind. Nat. Sci. Acad. 38, 350–354 (1972).
Google Scholar 
Subba Rao, B. V. S. S. R. Volume regulation in a euryhaline oligochaete, Pontodrilus bermudensis Beddard. Proc. Indian Acad. Sci. 87, 339–347 (1978).Article 

Google Scholar 
Owojori, O. J., Reinecke, A. J. & Rozanov, A. B. Effects of salinity on partitioning, uptake and toxicity of zinc in the earthworm Eisenia fetida. Soil Biol. Biochem. 40, 2385–2393. https://doi.org/10.1016/j.soilbio.2008.05.019 (2008).Article 
CAS 

Google Scholar 
Seesamut, T. et al. Occurrence of bioluminescent and nonbioluminescent species in the littoral earthworm genus Pontodrilus. Sci. Rep. 11, 8407 (2021).Article 
CAS 

Google Scholar 
Sivinski, J. & Forrest, T. Luminous defense in an earthworm. Fla. Entomol. 66, 517 (1983).Article 

Google Scholar 
Verdes, A. & Gruber, D. F. Glowing worms: Biological, chemical, and functional diversity of bioluminescent annelids. Integr. Comp. Biol. 57, 18–32. https://doi.org/10.1093/icb/icx017 (2017).Article 
CAS 

Google Scholar 
Shimomura, O. & Yampolsky, I. Bioluminescence: Chemical Principles and Methods 3rd edn. (World Scientific, 2019).Book 

Google Scholar 
Easton, E. G. Earthworms (Oligochaeta) from islands of the south-western Pacific, and a note on two species from Papua New Guinea. N. Z. J. Zool. 11, 111–128. https://doi.org/10.1080/03014223.1984.10423750 (1984).Article 

Google Scholar 
Shen, H.-P., Tsai, S.-C. & Tsai, C.-F. Occurrence of the earthworms Pontodrilus litoralis (Grube, 1855), Metaphire houlleti (Perrier, 1872), and Eiseniella tetraedra (Savigny, 1826) from Taiwan. Taiwania 50, 11–21 (2005).
Google Scholar 
Satheeshkumar, P., Khan, A. B. & Senthilkumar, D. Annelida, Oligochaeta, Megascolecidae, Pontodrilus litoralis (Grupe, 1985): First record from Pondicherry mangroves, southeast coast of India. Int. J. Zool. Res. 7, 406–409. https://doi.org/10.3923/ijzr.2011.406.409 (2011).Article 

Google Scholar 
Nguyen, T. T., Nguyen, D. A., Tran, T. T. B. & Blakemore, R. J. A comprehensive checklist of earthworm species and subspecies from Vietnam (Annelida: Clitellata: Oligochaeta: Almidae, Eudrilidae, Glossoscolecidae, Lumbricidae, Megascolecidae, Moniligastridae, Ocnerodrilidae, Octochaetidae). Zootaxa 4140, 1–92. https://doi.org/10.11646/zootaxa.4140.1.1 (2016).Article 

Google Scholar 
Chen, S.-Y., Hsu, C.-H. & Soong, K. How to cross the sea: Testing the dispersal mechanisms of the cosmopolitan earthworm Pontodrilus litoralis. R. Soc. Open Sci. 8, 202297. https://doi.org/10.1098/rsos.202297 (2021).Article 
ADS 

Google Scholar 
Smyth, K. & Elliott, M. Effects of changing salinity on the ecology of the marine environment. In Stressors in the Marine Environment (eds Solan, M. & Whiteley, N. M.) 161–175 (Oxford University Press, 2016).Chapter 

Google Scholar 
Veiga, M. P. T., Gutierre, S. M. M., Castellano, G. C. & Freire, C. A. Tolerance of high and low salinity in the intertidal gastropod Stramonita brasiliensis (Muricidae): Behaviour and maintenance of tissue water content. J. Molluscan Stud. 82, 154–160. https://doi.org/10.1093/mollus/eyv044 (2016).Article 

Google Scholar 
Carley, W. W., Caracciolo, E. A. & Mason, R. T. Cell and coelomic fluid volume regulation in the earthworm Lumbricus terrestris. Comp. Biochem. Physiol. 74, 569–575 (1983).Article 

Google Scholar 
Sharif, F. et al. Salinity tolerance of earthworms and effects of salinity and vermi amendments on growth of Sorghum bicolor. Arch. Agron. Soil Sci. 62, 1169–1181. https://doi.org/10.1080/03650340.2015.1132838 (2016).Article 
CAS 

Google Scholar 
Wu, Z. et al. Effects of salinity on earthworms and the product during vermicomposting of kitchen wastes. Int. J. Environ. Res. Public Health 16, 4737. https://doi.org/10.3390/ijerph16234737 (2019).Article 
CAS 

Google Scholar 
Oglesby, L. C. Volume regulation in aquatic invertebrates. J. Exp. Zool. 215, 289–301 (1981).Article 
CAS 

Google Scholar 
Generlich, O. & Giere, O. Osmoregulation in two aquatic oligochaetes from habitats with different salinity and comparison to other annelids. Hydrobiologia 334, 251–261. https://doi.org/10.1007/BF00017375 (1996).Article 

Google Scholar 
Carregosa, V. et al. Tolerance of Venerupis philippinarum to salinity: Osmotic and metabolic aspects. Comp. Biochem. Physiol. A 171, 36–43. https://doi.org/10.1016/j.cbpa.2014.02.009 (2014).Article 
CAS 

Google Scholar 
Freitas, R. et al. The effects of salinity changes on the polychaete Diopatra neapolitana: Impacts on regenerative capacity and biochemical markers. Aquat. Toxicol. 163, 167–176. https://doi.org/10.1016/j.aquatox.2015.04.006 (2015).Article 
CAS 

Google Scholar 
Rivera-Ingraham, G. A. & Lignot, J. H. Osmoregulation, bioenergetics and oxidative stress in coastal marine invertebrates: Raising the questions for future research. J. Exp. Biol. 220, 1749–1760. https://doi.org/10.1242/jeb.135624 (2017).Article 

Google Scholar 
Munnoli, P. M. & Bhosle, S. Effect of soil cow dung proportion of vermicomposting. J. Sci. Ind. Res. 68, 57–60 (2009).
Google Scholar  More

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Google Scholar  More

MAG binning and general featuresFrom the four hot springs, we assembled four associated metagenomes and then binned a total of 42 MAGs. We recovered 12 MAGs from RB10 hot spring, 13 from RB13, 14 from RB32 and 3 from RB108. Out of these 42 MAGs, 7 were of high-quality, 25 of nearly-high quality, 9 of medium quality and 1 of low quality (Table 1) based on metagenomic standards26. The GC% was quite variable, ranging from 25.76 to 70.35% among all MAGs and between 32.15 and 69.21% only among the high- and near high-quality MAGs. With the exception of RB108 from which we only recovered bacterial MAGs, we retrieved both bacterial and archaeal MAGs in the other hot springs. Two thirds of the MAGs (26/42) were assigned to the domain Bacteria and the rest to the domain Archaea (16/42) (Table 2).Table 1 General characteristics of the 42 MAGs obtained from RB10, RB13, RB32 and RB108 samples.Full size tableTable 2 Classification of the MAGs based on the taxonomic classification of GTDB-Tk software (v2.1.0) and the Genome Taxonomy Database (07-RS207 release).Full size tableTaxonomic and phylogenomic analyses of MAGsThe taxonomic affiliation of the MAGs was investigated in detail through the workflow classify of GTDB-Tk (v 2.1.0; GTDB reference tree 07-RS207) (Table 2) and through de novo phylogenomic analyses (Fig. S1a–i). We also tried to classify MAGs on the basis of overall genome relatedness indices (OGRI), which is detailed in supplementary material (Text S1, Table S2, Fig. S2).De novo phylogenomic analyses globally confirmed the positioning of MAGs provided by GTDB-Tk, with high branching support. For Bacteria, GTDB-Tk analyses allowed us to place the MAGs in the following clades: six in the phylum Aquificota from the four different springs, comprising four MAGs belonging to the genus Hydrogenivirga (family Aquificaceae) (RB10-MAG07, RB13-MAG10, RB32-MAG07, RB108-MAG02), and two belonging to the family ‘Hydrogenobaculaceae’ (RB10-MAG12, RB32-MAG11) (Table 2, Fig. S1a). Their closest cultured relatives originated either from hot springs or from deep-sea hydrothermal vents27. Three MAGs from three geothermal springs belonged to the phylum Armatimonadota (RB10-MAG03, RB13-MAG04, RB32-MAG03) and had no close cultured relatives. Seven MAGs have been classified into the phylum Chloroflexota: three MAGs belonging to the genus Thermoflexus from three different springs (RB10-MAG04, RB13-MAG05, RB32-MAG02), one affiliating with the genus Thermomicrobium (RB32-MAG08), one falling into the family Ktedonobacteraceae (RB108-MAG03), one belonging to the class Dehalococcoidia (RB32-MAG04) and another one whose phylogenetic position is more difficult to assert because it is a MAG of medium quality (RB32-MAG14). Six MAGs from four various hot springs belonged to the phylum Deinococcota, and to the genera Thermus (RB10-MAG08, RB10-MAG11, RB13-MAG09, RB32-MAG10, RB108-MAG01) and Meiothermus (RB13-MAG13). One MAG belonged to the family ‘Sulfurifustaceae’ (RB13-MAG01), in the phylum Proteobacteria (Gamma-class). The MAG referenced as RB32-MAG13 was classified into the phylum ‘Patescibacteria’, in the class ‘Paceibacteria’, and was distantly related to MAGs originating from groundwater and from hot springs. Finally, two MAGs from two different springs belonged to the phylum WOR-3, in the Candidatus genus ‘Caldipriscus’ (RB32-MAG12, RB10-MAG09).For Archaea, almost all the MAGs reconstructed in this study, e.g. 15 of the 16 archaeal MAGs, belonged to the phylum Thermoproteota. Among them, four belonged to the genus Ignisphaera (RB10-MAG05, RB13-MAG08, RB13-MAG11, RB32-MAG05), three to the genus Infirmifilum (RB10-MAG06, RB13-MAG03, RB32-MAG09), two to the genus Zestosphaera (RB10-MAG02, RB13-MAG06), three to the family Acidilobaceae (RB10-MAG01, RB13-MAG02, RB32-MAG01) and two to the order Geoarchaeales (RB10-MAG10, RB32-MAG06). Additionally, one belonged to the family Thermocladiaceae (RB13-MAG07). Lastly, the MAG belonging to another phylum (RB13-MAG12) was affiliated with the ‘Aenigmatarchaeota’, class ‘Aenigmatarchaeia’, and was distantly related to MAGs from hot springs and from deep-sea hydrothermal vent sediments28,29.Out of these 42 MAGs, at least 19 MAGs corresponded to different taxa at the taxonomic rank of species or higher according to GTDB (Table 2). Eighteen of them belonged to lineages with several cultivated representatives including the species Thermus thermophilus. 13 new genomic species within the GTDB genera Hydrogenivirga, HRBIN17, Thermoflexus, SpSt-223, CADDYT01, Zestosphaera, Ignisphaera, Infirmifilum, Thermus, Thermus_A, Meiothermus_B, JAHLMO01 and Caldipriscus, and 6 putative new genomic genera belonging to the GTDB families Hydrogenobaculaceae, Acidilobaceae, WAQG01, Thermocladiaceae, Sulfurifustaceae and HR35 could be identified (Table 2). In addition, 9 MAGs belonged to lineages that are predominantly or exclusively known through environmental DNA sequences. Thus, these 42 MAGs comprised a broad phylogenetic range of Bacteria and Archaea at different levels of taxonomic organization, of which a large majority were not reported before.The approaches implemented here were not intended to describe the microbial diversity present in these sources in an exhaustive way or to compare them in a fine way, and cannot allow it because of a 2-year storage at 4 °C. This long storage has probably led to changes in the microbial communities and to the selective loss or enrichment of some taxa. As a result, no analysis of abundance or absence of taxa can be conducted from these metagenomes and the results are discussed taking this bias into account. However, they do provide an overview of the microbial diversity effectively present. If we compare the phylogenetic diversity of the MAGs found in the four hot springs, we can observe that 3 shared phyla (Deinococcota, Aquificota and Chloroflexota: phyla names according to GTDB), 2 shared families (Thermaceae and Aquificaceae), and one shared genus (Hydrogenivirga) were found among the four sources (Fig. 2). In addition, hot springs RB10, RB13 and RB32, that are geographically close ( More

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Waples, R. S. & Gaggiotti, O. INVITED REVIEW: What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol. Ecol. 15, 1419–1439 (2006).Article 
CAS 

Google Scholar 
Mendez, M., Rosenbaum, H. C., Subramaniam, A., Yackulic, C. & Bordino, P. Isolation by environmental distance in mobile marine species: Molecular ecology of franciscana dolphins at their southern range. Mol. Ecol. 19, 2212–2228 (2010).Article 
CAS 

Google Scholar 
De Meeûs, T. et al. Population genetics and molecular epidemiology or how to “débusquer la bête”. Infect. Genet. Evol. 7, 308–332 (2007).Article 

Google Scholar 
Durigan, M. et al. Population genetic analysis of Giardia duodenalis: Genetic diversity and haplotype sharing between clinical and environmental sources. MicrobiologyOpen 6, e00424 (2017).Article 

Google Scholar 
Amaral, A. R. et al. Seascape genetics of a globally distributed, highly mobile marine mammal: The short-beaked common dolphin (genus Delphinus). PLoS ONE 7, e31482 (2012).Article 
ADS 
CAS 

Google Scholar 
Mendez, M. et al. Molecular ecology meets remote sensing: Environmental drivers to population structure of humpback dolphins in the Western Indian Ocean. Heredity 107, 349–361 (2011).Article 
CAS 

Google Scholar 
de los Angeles Bayas-Rea, R., Félix, F. & Montufar, R. Genetic divergence and fine scale population structure of the common bottlenose dolphin (Tursiops truncatus, Montagu) found in the Gulf of Guayaquil. Ecuador. PeerJ 6, e4589 (2018).Article 

Google Scholar 
Natoli, A., Peddemors, V. M. & Rus Hoelzel, A. Population structure and speciation in the genus Tursiops based on microsatellite and mitochondrial DNA analyses. J. Evol. Biol. 17, 363–375 (2004).Article 
CAS 

Google Scholar 
Oliveira, L. R., Loizaga De Castro, R., Cárdenas-Alayza, S. & Bonatto, S. L. Conservation genetics of South American aquatic mammals: An overview of gene diversity, population structure, phylogeography, non-invasive methods and forensics. Mammal Rev. 42, 275–303 (2012).Article 

Google Scholar 
Vollmer, N. L. & Rosel, P. E. Fine-scale population structure of common bottlenose dolphins (Tursiops truncatus) in offshore and coastal waters of the US Gulf of Mexico. Mar. Biol. 164, 1–15 (2017).Article 

Google Scholar 
MacLeod, C. D. Global climate change, range changes and potential implications for the conservation of marine cetaceans: A review and synthesis. Endanger. Species Res. 7, 125–136 (2009).Article 

Google Scholar 
Hartl, D. L., Clark, A. G. & Clark, A. G. Principles of Population Genetics, Vol. 116 (Sinauer associates Sunderland, 1997).Thomas, C. D. et al. Extinction risk from climate change. Nature 427, 145 (2004).Article 
ADS 
CAS 

Google Scholar 
Reeves, R. R., Smith, B. D., Crespo, E. A. & Notarbartolo di Sciara, G. Dolphins, whales and porpoises: 2002–2010 conservation action plan for the world’s cetaceans, Vol. 58 (IUCN, 2003).Crespo, E. A. & Hall, M. A. In Marine Mammals, 463–490 (Springer, 2002).Crespo, E. A. et al. Direct and indirect effects of highseas fisheries on the marine mammal populations in the northern and central Patagonian coast. J. Northwest Atl. Fish. Sci. 22, 189–207 (1997).Article 

Google Scholar 
Harlin-Cognato, A. D., Markowitz, T., Würsig, B. & Honeycutt, R. L. Multi-locus phylogeography of the dusky dolphin (Lagenorhynchus obscurus): Passive dispersal via the west-wind drift or response to prey species and climate change?. BMC Evol. Biol. 7, 1–17 (2007).Article 

Google Scholar 
Hoelzel, A. Evolution of population genetic structure in marine mammal species. In Population genetics for animal conservation, 294–318 (Cambridge University Press, Cambridge, 2009).Fraser, C. I., Nikula, R., Ruzzante, D. E. & Waters, J. M. Poleward bound: Biological impacts of Southern Hemisphere glaciation. Trends Ecol. Evol. 27, 462–471 (2012).Article 

Google Scholar 
Louis, M. et al. Influence of past climate change on phylogeography and demographic history of narwhals, Monodon monoceros. Proc. R. Soc. B 287, 20192964 (2020).Article 
CAS 

Google Scholar 
Skovrind, M. et al. Circumpolar phylogeography and demographic history of beluga whales reflect past climatic fluctuations. Mol. Ecol. 30, 2543–2559 (2021).Article 

Google Scholar 
Foote, A. D. et al. Ancient DNA reveals that bowhead whale lineages survived Late Pleistocene climate change and habitat shifts. Nat. Commun. 4, 1–7 (2013).Article 

Google Scholar 
Crespo, E. A. et al. Status, population trend and genetic structure of South American fur seals, Arctocephalus australis, in southwestern Atlantic waters. Mar. Mamm. Sci. 31, 866–890 (2015).Article 

Google Scholar 
Feijoo, M., Lessa, E. P., De Castro, R. L. & Crespo, E. A. Mitochondrial and microsatellite assessment of population structure of South American sea lion (Otaria flavescens) in the Southwestern Atlantic Ocean. Mar. Biol. 158, 1857–1867 (2011).Article 

Google Scholar 
Túnez, J. I., Cappozzo, H. L., Nardelli, M. & Cassini, M. H. Population genetic structure and historical population dynamics of the South American sea lion, Otaria flavescens, in north-central Patagonia. Genetica 138, 831–841 (2010).Article 

Google Scholar 
Oliveira, L., Ott, P. H., Grazziotin, F. G., White, B. & Bonatto, S. In Paper (SC/S11/RW26) presented to the Southern Right Whale Assessment Workshop (Commission International Whaling).Loizaga de Castro, R., Dans, S. L. & Crespo, E. A. Spatial genetic structure of dusky dolphin, Lagenorhynchus obscurus, along the argentine coast: Preserve what scale?. Aquat. Conserv. Mar. Freshw. Ecosyst. 26, 173–183 (2016).Article 

Google Scholar 
Pimper, L. E., Goodall, R. N. P. & Remis, M. I. First mitochondrial DNA analysis of the spectacled porpoise (Phocoena dioptrica) from Tierra del Fuego, Argentina. Mamm. Biol. 77, 459–462 (2012).Article 

Google Scholar 
Pichler, F. B. et al. Origin and radiation of Southern Hemisphere coastal dolphins (genus Cephalorhynchus). Mol. Ecol. 10, 2215–2223 (2001).Article 
CAS 

Google Scholar 
Dawson, S. M. In Encyclopedia of Marine Mammals, 166–172 (Elsevier, 2018).Robineau, D., Goodall, R. N. P., Pichler, F. & Baker, C. S. Description of a new subspecies of Commerson’s dolphin, Cephalorhynchus commersonii (Lacépède, 1804), inhabiting the coastal waters of the Kerguelen Islands. Mammalia 71, 172–180 (2007).Article 

Google Scholar 
Crespo, E. A. et al. Cephalorhynchus commersonii, Commerson’s Dolphin. IUCN; The IUCN Red List of Threatened Species; 10-2017; 1-14 (2017).Goodall, R. Commerson’s dolphin Cephalorhynchus commersonii (Lacépède 1804). Handb. Mar. Mamm. 5, 241–267 (1994).
Google Scholar 
Coscarella, M. A. Ecologıa, comportamiento y evaluación del impacto de embarcaciones sobre manadas de tonina overa Cephalorhynchus commersonii en Bahıa Engano, Chubut (Universidad de Buenos Aires, Buenos Aires, 2005).Dellabianca, N. A. et al. Spatial models of abundance and habitat preferences of commerson’s and peale’s dolphin in southern patagonian waters. PLoS ONE 11, e0163441 (2016).Article 

Google Scholar 
Goodall, R. et al. Studies of Commerson’s dolphins, Cephalorhynchus commersonii, off Tierra del Fuego, 1976–1984. Report of the International Whaling Commission (Special Issue 9), 143–160 (1988).White, R. The Distribution of Seabirds and Marine Mammals in Falkland Islands Waters (Joint Nature Conservation Committee, 2002).Loizaga de Castro, R., Dans, S. L., Coscarella, M. A. & Crespo, E. A. Living in an estuary: Commerson’s dolphin (Cephalorhynchus commersonii (Lacépède, 1804)), habitat use and behavioural pattern at the Santa Cruz River, Patagonia, Argentina. Latin Am. J. Aquat. Res. 41, 985–991 (2013).Article 

Google Scholar 
Pedraza, S. Ecología poblacional de la tonina overa, Cephalorhynchus commersonii, (Lacépède, 1804) en el litoral patagónico. Unpublished PhD thesis, Universidad de Buenos Aires, Buenos Aires, Argentina (2008).Garaffo, G. V. et al. Modeling habitat use for dusky dolphin and Commerson’s dolphin in Patagonia. Mar. Ecol. Prog. Ser. 421, 217–227 (2011).Article 
ADS 

Google Scholar 
Cipriano, F., Hevia, M. & Iñíguez, M. Genetic divergence over small geographic scales and conservation implications for Commerson’s dolphins (Cephalorhynchus commersonii) in southern Argentina. Mar. Mamm. Sci. 27, 701–718 (2011).Article 
CAS 

Google Scholar 
Pimper, L. E., Baker, C. S., Goodall, R. N. P., Olavarría, C. & Remis, M. I. Mitochondrial DNA variation and population structure of Commerson’s dolphins (Cephalorhynchus commersonii) in their southernmost distribution. Conserv. Genet. 11, 2157–2168 (2010).Article 

Google Scholar 
O’Brien, S. J. A role for molecular genetics in biological conservation. Proc. Natl. Acad. Sci. 91, 5748–5755 (1994).Article 
ADS 
CAS 

Google Scholar 
Loizaga de Castro, R., Hoelzel, A. & Crespo, E. Behavioural responses of Argentine coastal dusky dolphins (Lagenorhynchus obscurus) to a biopsy pole system. Anim. Welf. 22, 13–23 (2013).Article 
CAS 

Google Scholar 
Elphinstone, M. S., Hinten, G. N., Anderson, M. J. & Nock, C. J. An inexpensive and high-throughput procedure to extract and purify total genomic DNA for population studies. Mol. Ecol. Notes 3, 317–320 (2003).Article 
CAS 

Google Scholar 
Bérubé, M. & Palsbøll, P. Identification of sex in cetaceans by multiplexing with three ZFX and ZFY specific primers. Mol. Ecol. 5, 283–287 (1996).Article 

Google Scholar 
Hoelzel, A., Hancock, J. & Dover, G. Evolution of the cetacean mitochondrial D-loop region. Mol. Biol. Evol. 8, 475–493 (1991).CAS 

Google Scholar 
Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).Article 
CAS 

Google Scholar 
Ruzzante, D. E. et al. Validation of close-kin mark–recapture (CKMR) methods for estimating population abundance. Methods Ecol. Evol. 10, 1445–1453 (2019).Article 

Google Scholar 
Faircloth, B. C., Branstetter, M. G., White, N. D. & Brady, S. G. Target enrichment of ultraconserved elements from arthropods provides a genomic perspective on relationships among H ymenoptera. Mol. Ecol. Resour. 15, 489–501 (2015).Article 
CAS 

Google Scholar 
Faircloth, B. C. MSATCOMMANDER: Detection of microsatellite repeat arrays and automated, locus-specific primer design. Mol. Ecol. Resour. 8, 92–94 (2008).Article 
CAS 

Google Scholar 
Zhan, L. et al. MEGASAT: Automated inference of microsatellite genotypes from sequence data. Mol. Ecol. Resour. 17, 247–256 (2017).Article 
CAS 

Google Scholar 
Nei, M. Molecular Evolutionary Genetics (Columbia University Press, 1987).Librado, P. & Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).Article 
CAS 

Google Scholar 
Schneider, S., Roessli, D. & Excoffier, L. Arlequin: A software for population genetics data analysis, version 2.000. Genetics Biometry Laboratory, Department of Anthropology, University of Geneva, Switzerland (2000).Excoffier, L., Smouse, P. E. & Quattro, J. M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131, 479–491 (1992).Article 
CAS 

Google Scholar 
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 9, 772 (2012).Article 
CAS 

Google Scholar 
Excoffier, L. & Lischer, H. E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour 10, 564–567 (2010).Article 

Google Scholar 
Bandelt, H.-J., Forster, P. & Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).Article 
CAS 

Google Scholar 
Fu, Y.-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147, 915–925 (1997).Article 
CAS 

Google Scholar 
Rogers, A. R. & Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569 (1992).CAS 

Google Scholar 
Peakall, R. & Smouse, P. E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6, 288–295 (2006).Article 

Google Scholar 
Mantel, N. The detection of disease clustering and a generalized regression approach. Can. Res. 27, 209–220 (1967).CAS 

Google Scholar 
Drummond, A. J. & Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7, 214 (2007).Article 

Google Scholar 
Harlin, A. D., Markowitz, T., Baker, C. S., Würsig, B. & Honeycutt, R. L. Genetic structure, diversity, and historical demography of New Zealand’s dusky dolphin (Lagenorhynchus obscurus). J. Mammal. 84, 702–717 (2003).Article 

Google Scholar 
Rambaut, A., Suchard, M., Xie, D. & Drummond, A. Tracer v1. 6. http://beast.bio.ed.ac.uk/Tracer (2014).Van Oosterhout, C., Hutchinson, W. F., Wills, D. P. & Shipley, P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4, 535–538 (2004).Article 

Google Scholar 
Rice, W. R. Analyzing tables of statistical tests. Evolution 43, 223–225 (1989).
Google Scholar 
Goudet, J. FSTAT, a program to estimate and test gene diversities and fixation indices, version 2.9. 3. http://www2.unil.ch/popgen/softwares/fstat.htm (2001).Waples, R. S. & Do, C. LDNE: A program for estimating effective population size from data on linkage disequilibrium. Mol. Ecol. Resour. 8, 753–756 (2008).Article 

Google Scholar 
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).Article 
CAS 

Google Scholar 
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).Article 
CAS 

Google Scholar 
Earl, D. A. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361 (2012).Article 

Google Scholar 
Queller, D. C. & Goodnight, K. F. Estimating relatedness using genetic markers. Evolution 43, 258–275 (1989).
Google Scholar 
Wilson, G. A. & Rannala, B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics 163, 1177–1191 (2003).Article 

Google Scholar 
Milinkovitch, M. C., Leduc, R., Tiedemann, R. & Dizon, A. In Marine Mammals: Biology and Conservation (ed Evans, P. G. H. & Raga, J. A.) 325–359 (Springer, 2002).Pichler, F. Population structure and genetic variation in Hector’s dolphin (Cephalorhynchus hectori), ResearchSpace@ Auckland (2001).Pichler, F. & Baker, C. Loss of genetic diversity in the endemic Hector’s dolphin due to fisheries-related mortality. Proc. R. Soc. Lond. Ser. B Biol. Sci. 267, 97–102 (2000).Article 
CAS 

Google Scholar 
Greenwood, P. J. Mating systems, philopatry and dispersal in birds and mammals. Anim. Behav. 28, 1140–1162 (1980).Article 

Google Scholar 
Chilvers, B. L. & Wilkinson, I. S. Philopatry and site fidelity of New Zealand sea lions (Phocarctos hookeri). Wildl. Res. 35, 463–470 (2008).Article 

Google Scholar 
Engelhaupt, D. et al. Female philopatry in coastal basins and male dispersion across the North Atlantic in a highly mobile marine species, the sperm whale (Physeter macrocephalus). Mol. Ecol. 18, 4193–4205 (2009).Article 
CAS 

Google Scholar 
Möller, L. M. & Beheregaray, L. B. Genetic evidence for sex-biased dispersal in resident bottlenose dolphins (Tursiops aduncus). Mol. Ecol. 13, 1607–1612 (2004).Article 

Google Scholar 
Jansen van Vuuren, B., Best, P., Roux, J. P. & Robinson, T. Phylogeographic population structure in the Heaviside’s dolphin (Cephalorhynchus heavisidii): Conservation implications. Anim. Conserv. 5, 303–307 (2002).Article 

Google Scholar 
Pérez-Alvarez, M. J. et al. Microsatellite markers reveal strong genetic structure in the endemic Chilean dolphin. PLoS ONE 10, e0123956 (2015).Article 

Google Scholar 
Hamner, R. M., Pichler, F. B., Heimeier, D., Constantine, R. & Baker, C. S. Genetic differentiation and limited gene flow among fragmented populations of New Zealand endemic Hector’s and Maui’s dolphins. Conserv. Genet. 13, 987–1002 (2012).Article 

Google Scholar 
Pichler, F., Dawson, S., Slooten, E. & Baker, C. Geographic isolation of Hector’s dolphin populations described by mitochondrial DNA sequences. Conserv. Biol. 12, 676–682 (1998).Article 

Google Scholar 
Kraft, S. et al. From settlers to subspecies: Genetic differentiation in commerson’s Dolphins between South America and the Kerguelen Islands. Front. Mar. Sci. 8, 782512 (2021).Article 

Google Scholar 
Grant, W. & Bowen, B. W. Shallow population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. J. Hered. 89, 415–426 (1998).Article 

Google Scholar 
Ponce, J. F., Rabassa, J., Coronato, A. & Borromei, A. M. Palaeogeographical evolution of the Atlantic coast of Pampa and Patagonia from the last glacial maximum to the Middle Holocene. Biol. J. Lin. Soc. 103, 363–379 (2011).Article 

Google Scholar 
Wright, S. Isolation by distance. Genetics 28, 114 (1943).Article 
CAS 

Google Scholar 
Meirmans, P. G. Nonconvergence in B ayesian estimation of migration rates. Mol. Ecol. Resour. 14, 726–733 (2014).Article 

Google Scholar  More

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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