Ecology

1.Liu, J. et al. Research progress in Sapindus L. germplasm resources. World For. Res. 30, 15–21 (2017).
Google Scholar 
2.Sun, C., Jia, L., Xi, B., Wang, L. & Weng, X. Natural variation in fatty acid composition of Sapindus spp. seed oils. Ind. Crops Prod. 102, 97–104 (2017).CAS 
Article 

Google Scholar 
3.Liu, J. et al. Variation in fruit and seed properties and comprehensive assessment of germplasm resources of the genus Sapindus. Sci. Silva Sin. 55, 44–54 (2019).
Google Scholar 
4.Xu, Y., Jia, L., Chen, Z. & Gao, Y. Advances on triterpenoid Saponin of Sapindus mukorossi. Chem. Bull. 081, 1078–1088 (2018).CAS 

Google Scholar 
5.Basu, A., Basu, S., Bandyopadhyay, S. & Chowdhury, R. Optimization of evaporative extraction of natural emulsifier cum surfactant from Sapindus mukorossi: Characterization and cost analysis. Ind. Crops Prod. 77, 920–931 (2015).CAS 
Article 

Google Scholar 
6.Mukhopadhyay, S., Hashim, M. A., Sahu, J. N., Yusoff, I. & Gupta, B. S. Comparison of a plant based natural surfactant with SDS for washing of As(V) from Fe rich soil. Journal of Environmental Sciences 25 (2013).7.Mukhopadhyay, S. et al. Ammonium-based deep eutectic solvents as novel soil washing agent for lead removal. Chem. Eng. J. 294, 316–322 (2016).CAS 
Article 

Google Scholar 
8.Mukherjee, S. et al. Optimization of pulp fibre removal by flotation using colloidal gas aphrons generated from a natural surfactant. J. Taiwan Inst. Chem. Eng. 53, 15–21 (2015).CAS 
Article 

Google Scholar 
9.Shinobu-Mesquita, C. et al. Cellular structural changes in candida albicans caused by the hydroalcoholic extract from Sapindus saponaria L. Molecules 20, 9405–9418 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
10.Smułek, W. et al. Sapindus saponins’ impact on hydrocarbon biodegradation by bacteria strains after short- and long-term contact with pollutant. Colloids Surf. B 142, 207–213 (2016).Article 
CAS 

Google Scholar 
11.Rodríguez-Hernández, D. et al. Highly potent anti-leishmanial derivatives of hederagenin, a triperpenoid from Sapindus saponaria L. Eur. J. Med. Chem. 124, 153–159 (2016).PubMed 
Article 
CAS 

Google Scholar 
12.Rodriguez-Hernández, D., Demuner, A. J., Barbosa, L. C. A., Csuk, R. & Heller, L. Hederagenin as a triterpene template for the development of new antitumor compounds. Eur. J. Med. Chem. 105, 57–62 (2015).PubMed 
Article 
CAS 

Google Scholar 
13.Singh, P. T. D. & Singh, M. M. Anti-Trichomonas activity of Sapindus saponins, a candidate for development as microbicidal contraceptive. J. Antimicrob. Chemother. 62, 526–534 (2008).PubMed 
Article 
CAS 

Google Scholar 
14.Muntaha, S. T. & Khan, M. N. Natural surfactant extracted from Sapindus mukurossi as an eco-friendly alternate to synthetic surfactant: A dye surfactant interaction study. J. Clean. Prod. 93, 145–150 (2015).CAS 
Article 

Google Scholar 
15.Sun, C., Jia, L., Ye, H. O., Gao, Y. & Weng, X. Geographic variation evaluating and correlation with fatty acid composition of economic characters of Sapindus spp. fruits. J. Beijing For. Univ. 12, 73–83 (2016).
Google Scholar 
16.Barry, C. & Cox, P. D. Biogeography: An Ecological and Evolutionary Approach (Blackwell, 1980).
Google Scholar 
17.Stocker, T. F. et al. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of IPCC the Intergovernmental Panel on Climate Change, 95–123. http://www.ipcc.ch/publications_and_data/publications_ipcc_fourth_assessment_report_wg1_report_the_physical_science_basis.htm (2014).18.Hatfield, J. L. & Prueger, J. H. Temperature extremes: Effect on plant growth and development. Weather Clim. Extremes 10, 4–10 (2015).Article 

Google Scholar 
19.Ray, R., Gururaja, K. V. & Ramchandra, T. V. Predictive distribution modeling for rare Himalayan medicinal plant Berberis aristata DC. J. Environ. Biol. 32, 725–730 (2011).PubMed 

Google Scholar 
20.Lawler, J. J. Climate change adaptation strategies for resource management and conservation planning. Ann. N. Y. Acad. Sci. 1162, 79–98 (2009).ADS 
PubMed 
Article 

Google Scholar 
21.Richard et al. Will plant movements keep up with climate change? Trends Ecol. Evol. (2013).22.Araújo, M. B. & Peterson, A. T. Uses and misuses of bioclimatic envelope modeling. Ecology 93, 1527–1539 (2012).PubMed 
Article 

Google Scholar 
23.Elith, J. & Leathwick, J. R. Species distribution models: Ecological explanation and prediction across space and time. Annu. Rev. Ecol. Evol. Syst. 40, 677–697 (2009).Article 

Google Scholar 
24.Zhu, G., Liu, G., Bu, W. & Gao, Y. Ecological niche modeling and its applications in biodiversity conservation. Biodiv. Sci. 1, 94–102 (2013).
Google Scholar 
25.Stockwell, D. & Peters, D. P. The GARP modelling system: Problems and solutions to automated spatial prediction. Int. J. Geogr. Inf. Sci. 13, 143–158 (1999).Article 

Google Scholar 
26.Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259 (2006).Article 

Google Scholar 
27.Beaumont, L. J. et al. Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species’ current and future distribution. Ecol. Model. 186, 251–270 (2005).Article 

Google Scholar 
28.Liaw, A. & Wiener, M. Classification and regression by random forest. R News 23, 1–10 (2002).
Google Scholar 
29.Elith, J., Leathwick, J. R. & Hastie, T. A working guide to boosted regression trees. J. Anim. Ecol. 77, 802–814 (2008).CAS 
PubMed 
Article 

Google Scholar 
30.Elith, J. et al. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17, 43–57 (2011).Article 

Google Scholar 
31.Li, G., Du, S. & Wen, Z. Mapping the climatic suitable habitat of oriental arborvitae (Platycladus orientalis) for introduction and cultivation at a global scale. Sci. Rep. 6, 30009 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
32.Merow, C., Smith, M. J. & Silander, J. A. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography 36, 1058–1069 (2013).Article 

Google Scholar 
33.Merow, C. et al. What do we gain from simplicity versus complexity in species distribution models?. Ecography 37, 1267–1281 (2014).Article 

Google Scholar 
34.Wang, W. et al. Assessment of potential habitat for firmiana danxiaensis, a plant species with extremely small populations in danxiashan national nature reserve based on maxent model. Scientia Silvae Sinicae (2019).35.Guo, Y., Guo, J., Shen, X., Wang, G. & Wang, T. Predicting the bioclimatic habitat suitability of Ginkgo biloba L. in China with field-test validations. Forests https://doi.org/10.3390/f10080705 (2019).Article 

Google Scholar 
36.Huang, Z. et al. Geographic distribution and impacts of climate change on the suitable habitats of Zingiber species in China. Ind. Crops Prod. 138, 111429 (2019).Article 

Google Scholar 
37.Rong, Z. et al. Modeling the effect of climate change on the potential distribution of Qinghai Spruce (Picea crassifolia Kom.) in Qilian Mountains. Forests https://doi.org/10.3390/f10010062 (2019).Article 

Google Scholar 
38.Mohammadi, S., Ebrahimi, E., Moghadam, M. S. & Bosso, L. Modelling current and future potential distributions of two desert jerboas under climate change in Iran. Ecol. Inform. 52, 7–13 (2019).Article 

Google Scholar 
39.Ramos, R. S., Kumar, L., Shabani, F. & Picanco, M. C. Risk of spread of tomato yellow leaf curl virus (TYLCV) in tomato crops under various climate change scenarios. Agric. Syst. 173, 524–535. https://doi.org/10.1016/j.agsy.2019.03.020 (2019).Article 

Google Scholar 
40.Sultana, S., Baumgartner, J. B., Dominiak, B. C., Royer, J. E. & Beaumont, L. J. Impacts of climate change on high priority fruit fly species in Australia. PLoS ONE 15, e0213820 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
41.Raffini, F. et al. From nucleotides to satellite imagery: Approaches to identify and manage the invasive pathogen Xylella fastidiosa and its insect vectors in Europe. Sustainability 12, 4508 (2020).CAS 
Article 

Google Scholar 
42.Melo-Merino, S. M., Reyes-Bonilla, H. & Lira-Noriega, A. Ecological niche models and species distribution models in marine environments: A literature review and spatial analysis of evidence. Ecol. Model. https://doi.org/10.1016/j.ecolmodel.2019.108837 (2020).Article 

Google Scholar 
43.Sterne, T. K., Retchless, D., Allee, R. & Highfield, W. Predictive modelling of mesophotic habitats in the north-western Gulf of Mexico. Aquat. Conserv. Mar. Freshw. Ecosyst. https://doi.org/10.1002/aqc.3281 (2020).Article 

Google Scholar 
44.Convertino, M., Annis, A. & Nardi, F. Information-theoretic portfolio decision model for optimal flood management. Environ. Model. Softw. 119, 258–274 (2019).Article 

Google Scholar 
45.Ardestani, E. G. & Mokhtari, A. Modeling the lumpy skin disease risk probability in central Zagros Mountains of Iran. Prev. Vet. Med. 176, 104887–104887. https://doi.org/10.1016/j.prevetmed.2020.104887 (2020).Article 
PubMed 

Google Scholar 
46.Hanafi-Bojd, A. A., Vatandoost, H. & Yaghoobi-Ershadi, M. R. Climate change and the risk of malaria transmission in Iran. J. Med. Entomol. 57, 50–64. https://doi.org/10.1093/jme/tjz131 (2020).Article 
PubMed 

Google Scholar 
47.Zhang, L., Jing, Z., Li, Z., Liu, Y. & Fang, S. Predictive modeling of suitable habitats for Cinnamomum Camphora (L.) presl using maxent model under climate change in China. Int. J. Environ. Res. Public Health https://doi.org/10.3390/ijerph16173185 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
48.Peng, L.-P. et al. Modelling environmentally suitable areas for the potential introduction and cultivation of the emerging oil crop Paeonia ostii in China. Sci. Rep. https://doi.org/10.1038/s41598-019-39449-y (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
49.Li, J., Fan, G. & He, Y. Predicting the current and future distribution of three Coptis herbs in China under climate change conditions, using the MaxEnt model and chemical analysis. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2019.134141 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
50.Zhang, K., Zhang, Y. & Tao, J. Predicting the potential distribution of Paeonia veitchii (Paeoniaceae) in China by incorporating climate change into a maxent model. Forests https://doi.org/10.3390/f10020190 (2019).Article 

Google Scholar 
51.Sun, C. et al. Genetic structure and biogeographic divergence among Sapindus species: An inter-simple sequence repeat-based study of germplasms in China. Ind. Crops Prod. 118, 1–10 (2018).CAS 
Article 

Google Scholar 
52.Mahar, K. S., Palni, L. M. S., Ranade, S. A., Pande, V. & Rana, T. S. Molecular analyses of genetic variation and phylogenetic relationship in Indian soap nut (Sapindus L.) and closely related taxa of the family Sapindaceae. Meta Gene, S2214540017300336 (2017).53.Li, J., Chang, H., Liu, T. & Zhang, C. The potential geographical distribution of Haloxylon across Central Asia under climate change in the 21st century. Agric. For. Meteorol. 275, 243–254. https://doi.org/10.1016/j.agrformet.2019.05.027 (2019).ADS 
Article 

Google Scholar 
54.Prevéy, J. S., Parker, L. E., Harrington, C. A., Lamb, C. T. & Proctor, M. F. Climate change shifts in habitat suitability and phenology of huckleberry (Vaccinium membranaceum). Agric. For. Meteorol. 280, 107803 (2020).ADS 
Article 

Google Scholar 
55.Adeyemi, T. O., Ogundipe, O. T. & Olowokudejo, J. D. Species distribution modelling of family Sapindaceae in West Africa. Int. J. Bot. 8, 45–49 (2012).Article 

Google Scholar 
56.Sun, C. et al. Association of fruit and seed traits of sapindus mukorossi germplasm with environmental factors in Southern China. Forests 8, 491 (2017).Article 

Google Scholar 
57.Sun, C. et al. Genetic diversity and association analyses of fruit traits with microsatellite ISSRs in Sapindus. J. For. Res. 30, 197–207 (2019).Article 
CAS 

Google Scholar 
58.Gao, Y. et al. Canopy characteristics and light distribution in Sapindus mukorossi Gaertn. are influenced by crown architecture manipulation in the hilly terrain of Southeast China. Sci. Hortic. 240, 11–22 (2018).Article 

Google Scholar 
59.Zhang, Y.-Q. et al. Spatio-temporal effects of canopy microclimate on fruit yield and quality of Sapindus mukorossi Gaertn. Sci. Hortic. 251, 136–149 (2019).CAS 
Article 

Google Scholar 
60.Liu, J. et al. Potential distribution and ecological characteristics of genus Sapindus in China based on MaxEnt model. Sci. Silvae Sin. 57, 1–12. https://doi.org/10.11707/j.1001-7488.20210501 (2021).Article 

Google Scholar 
61.Cabral, A. L., Sales, J. D. F., Barbosa, K. F., Rodrigues, A. A. & Filho, S. C. V. Dormancy breakage and germination in Sapindus saponaria L. seeds as a function of temperature and germination substrate. Semina 40, 3345–3358 (2019).CAS 

Google Scholar 
62.Jia, L. & Sun, C. Research progress of biodiesel tree Sapindus mukorossi. J. China Agric. Univ. 017, 191–196 (2012).
Google Scholar 
63.Wei, X., Dai, T., Liu, S. & Jia, L. Effects of formula fertilization on leaf nutrient dynamics and yield of Sapindus mukorossi Gaertn. J. Nanjing For. Univ. 42, 21–28 (2018).
Google Scholar 
64.Pal, A. K., Vaishnav, V., Meena, B., Pandey, N. & Rana, T. S. Adaptive fitness of Sapindus emarginatus Vahl populations towards future climatic regimes and the limiting factors of its distribution. Sci. Rep. 10, 1–11 (2020).Article 
CAS 

Google Scholar 
65.Jayasinghe, S. L. & Kumar, L. Modeling the climate suitability of tea [Camellia sinensis(L.) O. Kuntze] in Sri Lanka in response to current and future climate change scenarios. Agric. For. Meteorol. 272–273, 102–117 (2019).ADS 
Article 

Google Scholar 
66.He, X., Ning, X., Guo, Y. & Wei, H. Geographical distribution of Xanthoceras sorbifolia Bunge in China and predicting suitable area under the climate change scenario. Res. Agric. Modern. 40, 138–146 (2019).
Google Scholar 
67.GBIF.org (13 October 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.4d9kye.68.Veloz, S. D. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence-only niche models. J. Biogeogr. 36, 2290–2299 (2010).Article 

Google Scholar 
69.Boria, R. A., Olson, L. E., Goodman, S. M. & Anderson, R. P. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecol. Model. 275, 73–77 (2014).Article 

Google Scholar 
70.Brown, J. L. SDMtoolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 5, 694–700. https://doi.org/10.1111/2041-210x.12200 (2014).Article 

Google Scholar 
71.Beckmann, M. et al. glUV: A global UV-B radiation data set for macroecological studies. Methods Ecol. Evol. 5, 372 (2014).Article 

Google Scholar 
72.Mao, J. F. & Wang, X. R. Distinct niche divergence characterizes the homoploid hybrid speciation of Pinus densata on the Tibetan Plateau. Am. Nat. 177(4), 424–439 (2011).ADS 
PubMed 
Article 

Google Scholar 
73.Li, G., Xu, G., Guo, K. & Du, S. Mapping the global potential geographical distribution of black locust (Robinia Pseudoacacia L.) using herbarium data and a maximum entropy model. Forests 5, 2773–2792 (2014).Article 

Google Scholar 
74.Allouche, O., Tsoar, A. & Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43, 1223–1232 (2006).Article 

Google Scholar 
75.Elith, J. et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129–151 (2006).Article 

Google Scholar 
76.Phillips, S. J. Transferability, sample selection bias and background data in presence-only modelling: A response to Peterson et al. (2007). Ecography 31, 272–278 (2008).Article 

Google Scholar 
77.Araujo, M. B., Pearson, R. G., Thuiller, W. & Erhard, M. Validation of species–climate impact models under climate change. Glob. Change Biol. 11, 1504–1513 (2005).ADS 
Article 

Google Scholar 
78.Ikhumhen, H. O., Li, T., Lu, S. & Matomela, N. Assessment of a novel data driven habitat suitability ranking approach for Larus relictus specie using remote sensing and GIS. Ecol. Model. 432, 109221 (2020).Article 

Google Scholar 
79.Bosso, L., De Conno, C. & Russo, D. Modelling the risk posed by the zebra mussel Dreissena polymorpha: Italy as a case study. Environ. Manage. 60, 304–313 (2017).ADS 
PubMed 
Article 

Google Scholar  More

Botanical content of the Libri Picturati Brazilian collectionOur identifications of all plant illustrations are listed with their vernacular names, page numbers, and associated information on growth form, geographical origin, conservation and domestication status in Supplementary Dataset S1. From the entire collection of Brazilian plant illustrations in the Libri Picturati, we identified 198 taxa that are organized in the Theatrum, LP and MC as indicated in Supplementary Table S1. Between folios 729 and 731 of the Theatrum, an illustration of a tea plant (Camellia sinensis (L.) Kuntze) is glued, which was sent by Cleyer from Batavia (currently Jakarta, Indonesia), the headquarters of the Dutch East Indian Company. As it was inserted later in the Theatrum and not depicted in Brazil, we did not include it in our analysis. A few plants remained unidentified due to a lack of morphological characters, the limited quality of the drawing and/or the lack of references to written sources by Marcgrave or Piso7,8.Among the LP botanical watercolors, we identified 34 vascular plant species (38 images) with the Passifloraceae as the most represented family (five species, six images), followed by the Fabaceae (five species, five images). Among the MC plant drawings, we identified 26 vascular plant species (34 images) and the most represented families were the Cucurbitaceae (three species, seven images) and the Myrtaceae (three species, three images). Among the illustrated content of the Theatrum, we identified 162 vascular plant species (175 images) and one basidiomycete fungus (Copelandia cyanescens (Sacc.) Singer, Bolbitiaceae). Fungi were commonly placed within the plant kingdom until the mid-twentieth century. The most represented families among the illustrated content are the Fabaceae (22 species, 22 images), followed by the Solanaceae (10 species, 11 images), Lamiaceae (six species, six images) and Myrtaceae (six species, eight images). The Fabaceae is the most diverse plant family in the world41, while the Myrtaceae is one of the most rich-species woody plant family in the Atlantic Forest in Brazil42.Mentzel’s unfinished task: the intended botanical content of the Theatrum
The Theatrum also includes 206 empty folios, interleaved between 160 folios with plant illustrations (see example in Fig. 1). On most folios, vernacular names and references to the pages of the HNB and IURNM are written on the top center, often relating to one taxon, but sometimes referring to two taxa (Fig. 1). This occurs specially at the end of the collection, as if the maker had ended up with little space and somehow had to squeeze them in. Among these unillustrated folios, the vernacular plant names and references to Marcgrave and Piso’s sources allowed us to identify 196 vascular plant species (218 records) including five ferns from the families Drypteriaceae (one species), Polypodiaceae (one species) and Pteridaceae (three species); one alga (Sargassum tenuissimum (Endlicher & Diesling) Grunow, Phaeophyceae) and a marine sponge (Clathria cf. nicoleae Vieira de Barros, Santos & Pinheiro, Microcioniadae) (Supplementary Dataset S1). Considering that the study of spongiology (Porifera) did not develop until the mid-nineteenth century, these animal colonies must have been considered an aquatic plant because of the tree-like shape and the fact of living attached to the seabed. The most represented family that would correspond to the empty folios was the Fabaceae (29 unillustrated species, 33 records), followed by the Arecaceae (nine species, ten records), Solanaceae (nine species, nine records), and Asteraceae (seven species, seven records). Estimates of the intended botanical content (i.e., empty folios with references together with the illustrated folios) are shown in Table 1.Figure 1Similar vernacular names for related taxa and distinct taxa associated to the same vernacular name in the Theatrum Rerum Naturalium.Full size imageTable 1 Estimations of the botanical content of the Theatrum Rerum Naturalium, including empty and illustrated folios.Full size tableOn p. 139 of the Theatrum, the vernacular name Ambaibuna is written on an empty page without reference to Marcgrave’s or Piso’s books. The page with Ambaibuna is located between Ambaiba (p. 137), which corresponds to the illustration of Cecropia pachystachya Trécul, and a blank page with only the vernacular name Ambaitinga (p. 141), which corresponds to C. hololeuca Miq.7: 92,24 (Fig. 1). The Brazilian Cecropia species are known in Tupi-related languages as Ambauba, Ambauva or Umbaúba (https://dataplamt.org.br/), which are phonetically and morphologically similar to Ambaibuna. For those reasons, we initially assumed that Ambaibuna referred to a Cecropia species, but the same name Ambaibuna is later repeated together with the name Iito (p. 227) next to an illustration that represents a completely different tree species: Guarea guidonia (L.) Sleumer (Fig. 1). Furthermore, the name Ambaibuna is also written above the illustration of a grapevine, Vitis vinifera L. (p. 257), also unrelated to Cecropia (Fig. 1).Whether Ambaibuna was a generic name to designate several non-related species or represents a mistake by the author who wrote the names on the illustrations remains unknown. On the other hand, neither Marcgrave nor Piso mentioned Ambaibuna in their descriptions of the Brazilian flora. Aside from Marcgrave and Piso’s books7,8, it is yet to be determined which source(s) Mentzel relied on when arranging the botanical content of the Theatrum. It is nonetheless clear that he must have been confused by the similarity of some of the Tupi-related plant names. Unfortunately, Marcgrave was no longer present to help him match the illustrations, names and descriptions, because he died about 16 years before Mentzel started organizing the Brazilian plant illustrations.Origin of the exotic species in the Libri Picturati
The Libri Picturati collection depicts in its majority native Brazilian plants. Most of the species represented in the Theatrum are native from Brazil, but the proportion of native species is much lower in the MC and lowest in the LP, in which almost half of the illustrations represent introduced species (Fig. 2).Figure 2Proportion of native and introduced species in the Brazilian collection of the Libri Picturati: Theatrum Rerum Naturalium (Theatrum), Libri Principis (LP) and Miscellanea Cleyeri (MC).Full size imageThere are 35 species of exotic origin in the complete Brazilian collection of the Libri Picturati (Supplementary Table S2). These introduced species now occur in (sub-) tropical areas around the world. Most of the exotic plants originally came from other parts of the Americas, especially Mexico, the Caribbean and the Andes region (14 species); followed by those that originated in the African continent (10 species) and tropical Asia (nine species) (Supplementary Table S2). Most of the exotic American plants that were introduced to Brazil were domesticated and traded by indigenous groups long before the European colonization, such as papaya (Carica papaya L.), cotton (Gossypium barbadense L.), sweet potato (Ipomoea batatas (L.) Lam.), beans (Phaseolus vulgaris L.), guava (Psidium guajava L.) and maize (Zea mays L.)37. Most of the species of Asiatic origin were already naturalized or cultivated in Africa and introduced to Brazil by means of the Trans-Atlantic slave trade before the Dutch arrived, such as yams (Dioscorea alata L.), plantains (Musa × paradisiaca L.) and weeds like Abrus precatorius L. and Plumbago zeylanica L.43,44. Others were introduced from Europe by merchants and settlers, such as the Portuguese Jesuits, who incorporated them as remedies into their boticas (Jesuit pharmacies in the colonies). For example, the various Citrus and pomegranate fruits were not only planted as fruits but also used to expel roundworms and to combat cold fevers, respectively45: 88. Before their arrival to Brazil, the Portuguese and Dutch must have been familiar with some African plants, such as Aloe vera (L.) Burm.f., Ricinus communis L. and Tamarindus indica L. These useful plants were already known in Europe through Arabic and Greek medical texts, which knowledge was boosted by their translations into Latin during the High Renaissance45,46. Punica granatum L. was introduced into the Iberian Peninsula via ancient merchant routes in the Mediterranean47 and brough to Brazil by the Portuguese45. Grapes (Vitis vinifera L.) were already cultivated by the Portuguese in Pernambuco around 154248. Along the Atlantic coast, lemons, pomegranates and grape vines adapted to the new environmental conditions and thrived in the vicinities of Johan Maurits’ residence, as evidenced by the illustrations in the Theatrum and textual accounts6,7,40.The presence of these globally commodified plants is common today in Brazil as in many regions worldwide. Other species seem to have lost their popularity over time. The so-called Ethiopian, Guinean or Negro pepper, Xylopia aethiopica (Dunal) A.Rich., was present around the 1640s in northeast Brazil, as evidenced in the Libri Picturati by a painting with a fruiting branch with leaves named Piperis aethiopici spés (Fig. 3a). The first iconography of this aromatic tree in Europe is found in Matthioli’s commentaries on Dioscorides under the name of Piper aethiopicum49: 575 and its fruits were previously cited by the Persian polymath Avicenna (980–1037)30. In Europe, this African pepper was commonly used until southeast Asian spices gained popularity in the sixteenth century50. In the plantation societies of tropical America, X. aethiopica constituted a food crop for enslaved Africans in the early colonial period43: 135. Today, its fruits are used in aphrodisiac tonics51 and special dishes prepared for African deities (Orishas) in Cuba43: 90, but it is unclear whether the species grows in Brazil. Its current distribution range encompasses West, Central and Southern Africa (https://gbif.org/occurrence/map?taxon_key=3157151). The dry fruits are used in tropical Africa as a condiment, in rituals and as medicine to treat cough, bronchitis, rheumatism, malaria, amenorrhea and uterine fibroids52,53,54. There is an herbarium record in Brazil made by photographer and anthropologist Pierre Verger. The label on the specimen mentions ‘Brazil’ and ‘Plantas de Candomblé’ and it indicates that the voucher was deposited at the Herbarium Alexandre Real Costa (ALCB, according to Index Herbariorum: http://sweetgum.nybg.org/science/ih/, accessed 23 August 2021) in Bahia (Verger s/n, ALCB012478, available at ALCB, via Species Link: https://specieslink.net/search/, accessed 23 August 2021) Verger presumably collected this specimen in Bahia in 1967 while he was researching on ritual and medicinal plants used in Candomblé (http://inct.florabrasil.net/alcb-resgate/, accessed 2 June 2021)55. However, it seems to be a mixed collection, as the leaves are oppositely arranged and with long petioles, which is uncommon to Annonaceae30. In Brazil, the fruits of the Brazilian relative Xylopia aromatica (Lam.) Mart. have probably served as a good substitute for X. aethiopica, as they have a similar peppery taste and stomachic properties56: 3, and are more easily gathered from the cerrado savannahs or the Amazon rainforest. Voeks57 documented X. aethiopica seed powder as used in Candomblé rituals by Yoruba practitioners in Bahia. Nevertheless, there is no clear information whether X. aethiopica is cultivated in the Neotropics or imported; thus, the origin of the fruits, seeds or its powder in Brazil remains uncertain.Figure 3 (a) The African spice-producing tree Xylopia aethiopica depicted in the Theatrum Rerum Naturalium (p. 321); (b) The first record of the sunflower (Helianthus annuus) in Brazil (Theatrum: 555).Full size imageThe first reference to the sunflower (Helianthus annuus L.) in Brazil dates to the twentieth century, when it was introduced by European immigrants due to its economic value as an oil-producing crop58. Sunflowers are of North American and/or Mexican origin 59,60, and were introduced to Europe in the sixteenth century by the Spanish, as part of the Columbian exchange61. Merchants observed how native Americans benefited from this plant and exported the sunflower to Europe, where it was primarily valued as ornamental and later as a food crop, propelled by genetic improvement by the Russians in the 1800s59. Before the sunflower became a popular and well-stablished crop in the twentieth century, this plant was already encountered in northeast Brazil, as evidenced by the illustration in the Theatrum (Fig. 3b). Portuguese sailors may have played a role in its introduction to Brazilian territories or it could have been intentionally brought by merchants or Jesuits, although the latter paid more attention to medicinal plants45,62. We may also consider the Dutch as active agents in its introduction to their colonies in the northeast. A relevant female agent in the dissemination of the sunflower in the Netherlands was Christine Bertolf (1525–1590), who was acquainted with the Spanish court and keen of the rare plants that thrived in the Royal Botanical Garden in Madrid63. She spread textual and visual information about the sunflower, and possibly also its seeds, among her network of naturalists and collectors, including the Flemish botanists Dodoens and Clusius63. After Dodoens64: 295 depicted the first European sunflower in his herbal in 1568, images and descriptions of this species began to circulate in manuscripts of other naturalists and physicians in Europe (e.g., Matthioli65: preface, Fragoso66: title page, Monardes67: 109 and Clusius68: 14–15). Thus, by the seventeenth century, Dutch scholars and collectors of exotic naturalia were familiar with sunflowers, which possibly promoted its cultivation at Johan Maurits’ gardens with ornamental purposes.Interestingly, the sunflower is referenced as Camará-guaçú, an indigenous term from the macrolinguistic Tupi family. Camará, Kamará or Cambará is a generic name given to several unrelated species, such as Lantana camara L. (Verbenaceae) and Ageratum conyzoides L. (Asteraceae) (http://www.dataplamt.org.br/, accessed 2 June 2021), both found in the Theatrum (p. 341 and 343 respectively). According to Tibiriçá69, in Tupi caa means plant and mbaraá means illness, and according to Cherini70 Cambará means “leaf of rough bark”. Hence, Camara also refers to medicinal plants with rough leaves in general. Guaçú means big and miri small71, which matches with the larger inflorescence of H. annuus in contrast to the African weed Sida rhombifolia L. (Malvaceae), documented as Camara-miri in the HNB and “used by black people as a broom to sweep the houses of their masters”7: 110. According to the Tupi-based nomenclature associated to H. annuus in the Theatrum, Tupi indigenous groups were already familiar with the sunflower in Brazil around the 1640 s.Life forms and domestication status of the Libri Picturati plantsMost of the species in the Theatrum are tropical trees, followed by shrubs, herbs, and lianas (Fig. 4). Several are rainforest trees, such as Andira fraxinifolia Benth., Garcinia brasiliensis Mart. and Syagrus coronata (Mart.) Becc. The same trend was observed for the illustrations in the MC, with trees as the most often represented life forms, followed by shrubs, lianas and herbs. Typically, the LP contains much less trees, but more small herbs, shrubs and vines that were probably found in and around Mauritsstad (i.e., the former capital of Dutch Brazil, currently a part of the Brazilian city of Recife), such as Commelina erecta L. and Turnera subulata Sm., which commonly grow in disturbed landscapes.Figure 4Proportion of life forms of the species depicted in the Brazilian collection of the Libri Picturati: Theatrum Rerum Naturalium (Theatrum), Libri Principis (LP) and Miscellanea Cleyeri (MC).Full size imageAlthough the majority of the species depicted in the Theatrum and the MC are wild forest trees, some species are found both in the wild and cultivated, such as Psidium guineense Sw., which was part of the pre-Columbian anthropogenic forests or ‘indigenous landscape’ in Brazil37,38,72. Some trees were planted in or around Recife. Hancornia speciosa Gomes, known by its Tupi-based name Mangabiba or Mangaiba [Mangabeira]7: 121, was cultivated in Mauritsstad6: 242,40. The fruit of H. speciosa (Mangaba) was harvested in great amounts as it was a highly appreciated food7: 122. Seeds were collected to plant the tree, and Marcgrave gave details about the specific locations of varieties in different northeastern locations (Salvador, Sergipe and Olinda). H. speciosa was already selected and managed by indigenous groups before colonization37, yet wild populations of this tree are still found in the Brazilian rainforest and savannah (http://floradobrasil.jbrj.gov.br/reflora/floradobrasil/FB15558, accessed 4 June 2021).Domesticated plants are represented in higher proportions within the LP and the MC (Fig. 5), accounting mostly for introduced fruit species (Supplementary Dataset S1), such as Citrus spp., Musa x paradisiaca and Cocos nucifera L., which were cultivated in Maurits’ gardens in Recife40. The influence of the European colonization of Brazil is also visible by the presence of weeds from Asia and Africa among the illustrations in the Theatrum and the LP, such as Abrus precatorius L., Argemone mexicana L., Boerhavia coccinea Mill. and Plumbago zeylanica L. Some of these plants were introduced from Africa via the slave ships, while others may have dispersed naturally44. Guilandina bonduc L., an African scrambling shrub depicted as Inimboi in the Theatrum, was described by Piso7: 95 as “growing in abundance in sandy and dry forests of the coasts”. We categorized G. bonduc as a wild plant: its round seeds could have been brought by oceanic currents from West African shores and germinated in the coastal vegetation of Pernambuco and other South American areas73. However, G. bonduc may also have reached Brazil during the Trans-Atlantic slave trade, as the hard, grey seeds are used in the African game Oware and also used in bead ornaments74.Figure 5Domestication status of the species in the Brazilian collection of the Libri Picturati: Theatrum Rerum Naturalium (Theatrum), Libri Principis (LP) and Miscellanea Cleyeri (MC).Full size imagePlant parts represented in the Libri Picturati
The way plants are depicted in the Libri Picturati provides us with information about the level of botanical skills of the artists, and how closely they worked together with the naturalists in the Dutch colony. Some plants are represented by only loose parts or depicted sterile, while others show us different organs and reproductive stages, which greatly facilitated their taxonomic identification (Table 2).Table 2 Plant parts represented in the botanical illustrations of the Libri Picturati: Theatrum Rerum Naturalium (Theatrum), Libri Principis (LP) and Miscellanea Cleyeri (MC).Full size tableMost illustrations depict fertile plant species with flowers and fruits, often cut in half to show the seeds, which reveals a high level of botanical knowledge. Fertile plants are more common in the Theatrum, in a few occasions also showing their tubers, such as Spondias tuberosa Arruda, known as Umbi [Iva Umbu], of which the prominent tuber in the bottom front captures the attention of the observer (Fig. 6a). Likely associated to a scientific purpose, drawing some plant parts out of proportion corresponds to a pictorial style also observed in other iconographies. This is also the case in the Icones Plantarum Malabaricarum, which depicts plants from Ceylon (modern Sri Lanka) in the eighteenth century and often accentuates useful fruits, flowers or roots75.Figure 6 (a) Spondias tuberosa with the tuber painted in front of the branch with leaves, tiny white flowers and a detail of the immature (green) and mature (yellow) fruit in the back (Theatrum Rerum Naturalium: 261); (b) Infertile individual of Hippeastrum psittacinum (Theatrum: 389); c Ficus gomelleira leaf, probably picked from the ground (Theatrum: 157); (d) Flowering vine of Centrosema brasilianum (Libri Principis: 1); (e) Amphilophium crucigerum dry open fruit without seeds (Theatrum: 387).Full size imagePiso7: 78 indicated that roots [tubers] of S. tuberosa deserved special attention, because of the way they developed underground and their use as a refreshment [water reservoir] for feverish patients and exhausted travelers, as he experimented himself. He and Marcgrave7: 108 also described how its fruits were valued as food. This example not only provides textual and visual evidence of the field trips to the interior by these naturalists and their first-hand experiments, but also adds insights into the connectedness between artistic and scientific practices in seventeenth century Dutch Brazil. Currently S. tuberosa, known as Umbu or Umbuzeiro (https://dataplamt.org.br/), is an important economic and subsistence food resource for rural communities in semiarid regions of northeast Brazil76,77. Its specialized root system (xylopodia) bears tubers that store liquids, sugars and other nutrients and allow the survival of the tree during the dry seasons of the caatinga and central Brazilian savanna, where this species is endemic78. The water or sweet juice of these xylopodia is still used as an emergency thirst quencher in extreme arid areas of the Brazilian sertão79; also see https://www.youtube.com/watch?v=NyGNlrljAww, accessed 25 August 2012].In the Theatrum, a small proportion of plants are illustrated in their sterile stage, such as Hippeastrum psittacinum (Ker Gawl.) Herb. (Fig. 6b) or Ficus gomelleira Kunth & C.D.Bouché (Fig. 6c). Marcgrave7: 32 did not see the impressive flower of H. psittacinum as it is lacking in his observations25: 59. The Theatrum painting must have been made in the wet season in the interior of Pernambuco, when Marcgrave and the painter(s) encountered the lily with only leaves, before these fall off and make place for the mesmerizing flower25. Ficus gomelleira, depicted by a single oblong leaf with its characteristic pinnate venation (p. 157), is a large tree, up to 40 m tall (https://portal.cybertaxonomy.org/flora-guianas/node/3041, accessed 4 June 2021). It can be challenging to collect a branch, so the painter(s) or local assistants possibly picked a leaf from the ground (Fig. 6c).The LP contains mainly flowering plants (e.g., Ruellia cf. elegans Poir.), tendrillate vines (e.g., Centrosema brasilianum (L.) Benth. (Fig. 6d)) and cultivated crops, such as peanuts (Arachis hypogaea L.), pumpkins (Cucurbita pepo L.) or guava (Psidium guajava L.) (Supplementary Dataset S1). Compared to the MC and the LP, a smaller proportion of the illustrations display only flowers or fruits in the Theatrum. Yet, these deserve special attention as the reasons for only painting the reproductive organs in the three collections may differ. While in the MC and LP flowers or fruits represent species that are domesticated or more likely to be found in urban areas, such as Capsicum baccatum L. or Hancornia speciosa, the Theatrum contains more loose parts of native plants found in the rainforest. Amphilophium crucigerum (L.) L.G.Lohmann is a liana referenced by the Tupi-related name Iaruparicuraba (Theatrum: 387) and today known in Brazil as pente-de-macaco (https://dataplamt.org.br/) due to its large dehiscent fruit (c.17 cm long) that opens in two valves covered with soft spines, hence its name “monkey’s comb” (Fig. 6e). Its winged seeds are not present in the drawing, possibly one empty valve was gathered from the ground, and the seeds were already dispersed by the wind.In the MC, we also find some drawings of infertile structures, but these mostly belonged to species that were depicted on several folios. When assembling those folios, we observed the whole plant represented in its fertile stage: the watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai) is depicted with its leaves and fruit on folio 63 (verso) and its leaves on folio 64 (recto). In the case of Furcraea foetida (L.) Haw., whoever bounded the drawings in the MC collection did not realize that folios 63 (recto), 64 (verso) and double folio 68 formed together one entire plant (See Supplementary Fig. S1).On other occasions, the painters focused on painting the plant parts that were valuable to humans. Several rainforest trees were highly valued for its edible fruits or seeds, such as Hymenaea courbaril L.7: 101 or Lecythis pisonis L., of which the “seeds (also called chestnuts) were eaten raw or roasted”7: 128 and “were considered aphrodisiacs”7: 65. The fruit of Macoubea guianensis Aubl. was “appreciated for its sweetness by the indigenous peoples to eat during their travels, while Europeans used it to treat chest affections”8: 242. The fruit of Swartzia pickelii Killip ex Ducke was “not eaten unless it was cooked, from which the inhabitants made a wholesome delicacy for the stomach called Manipoy”8: 165. The same applies to the tomato-like fruits of the African eggplant Solanum aethiopicum L., which were “eaten cooked, after seasoning with oil and pepper; it has lemon taste”7: 24. While these plants are represented in the Theatrum only by their fruits (Supplementary Dataset S1), the tree branches or the whole plant are depicted in the written sources. The illustrations in the books were most likely made by Marcgrave, who aimed to describe and depict as many plant parts as possible, although compromising in aesthetic aspect. The painters, on the other hand, focused on the edible parts without sacrificing their aesthetics. In any case, the illustrations from the Theatrum and the woodcuts and descriptions in the HNB and IURNM often complemented each other and thus facilitated our identifications.Current conservation status of the Brazilian species in the Libri Picturati
In the past centuries, the Atlantic Forest and savannah regions of northeast Brazil have been severely affected by habitat loss and degradation due to the expansion of urbanization, intensive agriculture, farming and logging80,81. Several plant species that were abundant enough to be noted by European artists around 1640 are not common anymore today. According to the IUCN Red List, eight species in the Libri Picturati, seven in the Theatrum and one in the LP are currently experiencing population decline or are at risk of facing extinction (Supplementary Table S3). Several endemic plants from the northeast Atlantic rainforest and caatinga biomes appear in the illustrations. Four species in the Libri Picturati are currently CITES-listed and restricted to trade: the cacti Brasiliopuntia brasiliensis (Willd.) A.Berger, Cereus fernambucensis Lem., Epiphyllum phyllanthus (L.) Haw. and Melocactus violaceus subsp. margaritaceus N.P.Taylor The latter is an endemic cactus of the coastal sand dunes’ ecoregion in the Atlantic rainforest known as restinga, which is severely threatened by agricultural expansion and urbanization82.Some endemic species are classified as Least Concern by the IUCN or the CNC Flora (12 species), while others (13 species) have not been evaluated yet (Supplementary Dataset S1). The MC does not contain threatened species, but includes two endemic trees: Attalea compta Mart. and Eugenia cf. brasiliensis Lam., which are only found in the biodiversity hotspots of the Atlantic rainforest and the cerrado, both greatly affected by habitat loss23. The mangrove vegetation along the Brazilian coast has been severely affected by urbanization, pollution by industrial and domestic waste and climate change83,84, threatening the populations of the mangrove trees Avicennia schaueriana Stapf & Leechm. ex Moldenke and Laguncularia racemosa (L.) C.F.Gaertn. The occurrence of anthropogenic impacts and the lack of available data call for the implementation of more in-depth and continuous studies on the conservation status of these vulnerable populations.Linking the plant illustrations to the published works and Marcgrave’s herbariumA total of 357 different plant species is described in the HNB and IURNM (Supplementary Dataset S2). Because the Theatrum constitutes a larger number of illustrations, we found more taxa from the books and the herbarium represented in this source (102 out of 163, 63%). However, the largest overlap was found between the MC and the HNB / IURNM: 21 out of 26 taxa (81%) were also described in the books. A smaller overlap exists between the LP and the HNB / IURNM (18 out of 34, 53%). We counted 143 taxa at species level in Marcgrave’s herbarium (Supplementary Dataset S3) and we observed some of these preserved species in all three pictorial works, with the largest percentage of taxa in common with the MC (seven out of 26, 27%), probably because of its smaller number of images. Strikingly, a third of the species illustrated in the whole Brazilian collection of the Libri Picturati could not be ascribed to the species described by Marcgrave or Piso (61 out of 180, 34%). More

Acinetobacter isolates from the TGR regionIn the TGR region, we isolated 21 Acinetobacter strains (one A. johnsonii, one A. haemolyticus and 19 Acinetobacter sp. strains) (Table S1). Based on phylogenetic analysis after 16S rRNA gene sequencing, these 21 isolates belong to the genus Acinetobacter, exhibiting a similarity of 95.38–99.93% with known Acinetobacter strains in GenBank (Table S1). In phylogenetic tree (N-J) constructed with both isolated and known Acinetobacter strains, these 21 isolates branched deeply with three Acinetobacter clusters consisting of important clinical Acinetobacter species, such as A. johnsonii H10 (FJ009371), A. junii NH88-14 (FJ447529), A. baumannii ATCC19606T (HE651907), A. lwoffii DSM2403T (X81665) and A. haemolyticus TTH04-1 (KF704077) (Fig. 1). Five reference Acinetobacter strains were selected and used21. Currently, the genus Acinetobacter comprises 68 species with validly-published names (https://apps.szu.cz/anemec/Classification.pdf, May 25, 2021). Among the named species, A. baumannii is the most studied species associated with clinical infections followed by the non-A. baumannii species A. haemolyticus, A. junii, A. johnsonii, and A. lwofii21.Figure 1A phylogenetic tree of 16S rRNA gene sequences showing position of isolates among species of genus Acinetobacter. Both isolates from the TGR region (the bold fonts) and reference strains used to infect Caenorhabditis elegans (the red fonts) are shown.Full size imageEffect of different Acinetobacter strains isolated from the TGR region and reference strains on lifespan of nematodesL4-larvae were exposed to different Acinetobacter strains for 24-h. Totally 21 Acinetobacter strains isolated from the TGR region and 5 reference strains of Acinetobacter species were used for the lifespan analysis. Based on the comparison of lifespan curves, exposure to Acinetobacter strains of AC2, AC3, AC4, AC5, AC6, AC7, AC8, AC9, AC10, AC11, AC12, AC13, AC14, AC16, AC17, AC19, AC20, A. johnsonii H10, and A. haemolyticus TTH0-4 could not alter lifespan curve (Fig. 2). Similarly, Acinetobacter strains of AC2, AC3, AC4, AC5, AC6, AC7, AC8, AC9, AC10, AC11, AC12, AC13, AC14, AC16, AC17, AC19, AC20, A. johnsonii, and A. haemolyticus also could not influence mean lifespan (Fig. 2). Different from these, the lifespan curves of nematodes exposed to Acinetobacter strains of AC1, AC15, AC18, AC21, A. baumannii ATCC 19606T, A. junii NH88-14, and A. lwoffii DSM 2403T were significantly (P  More

1.Beyene, Y. et al. The characteristics and chronology of the earliest Acheulean at Konso, Ethiopia. Proc. Natl. Acad. Sci. USA 110, 1584–1591 (2013).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
2.Haslam, M. et al. Late Acheulean hominins at the marine isotope stage 6/5e transition in north-central India. Quat. Res. 75, 670–682 (2011).CAS 
Article 

Google Scholar 
3.Carotenuto, F. et al. Venturing out safely: The biogeography of Homo erectus dispersal out of Africa. J. Hum. Evol. 95, 1–12 (2016).MathSciNet 
CAS 
PubMed 
Article 

Google Scholar 
4.Benito, B. M. et al. The ecological niche and distribution of Neanderthals during the Last Interglacial. J. Biogeogr. 44, 51–61 (2017).Article 

Google Scholar 
5.Bae, C. J., Douka, K. & Petraglia, M. D. On the origin of modern humans: Asian perspectives. Science (80-) 358, 65 (2017).Article 
CAS 

Google Scholar 
6.Rogers, A. R., Harris, N. S. & Achenbach, A. A. Neanderthal-Denisovan ancestors interbred with a distantly related hominin. Sci. Adv. 6, 1–8 (2020).
Google Scholar 
7.Ruff, C. B., Trinkaus, E. & Holliday, T. W. Body mass and encephalization in Pleistocene Homo. Nature 387, 173–176 (1997).CAS 
PubMed 
Article 
ADS 

Google Scholar 
8.Rightmire, G. P. Brain size and encephalization in early to Mid-Pleistocene Homo. Am. J. Phys. Anthropol. 124, 109–123 (2004).PubMed 
Article 

Google Scholar 
9.Rightmire, G. P. Later Middle Pleistocene Homo. in Handbook of Palaeoanthropology (eds. Henke, W. & Tattersall, I.) 2221–2242 (Springer, 2015).10.Lahr, M. M. The complex landscape of recent human evolution. Science (80-) 372, 995 (2021).Article 
CAS 

Google Scholar 
11.Athreya, S. & Hopkins, A. Conceptual issues in hominin taxonomy: Homo heidelbergensis and an ethnobiological reframing of species. Am. J. Phys. Anthropol. https://doi.org/10.1002/ajpa.24330 (2021).Article 
PubMed 

Google Scholar 
12.Grun, R. et al. U-series and ESR analyses of bones and teeth relating to the human burials from Skhul. J. Hum. Evol. 49, 316–334 (2005).PubMed 
Article 

Google Scholar 
13.Harvati, K. et al. Apidima Cave fossils provide earliest evidence of Homo sapiens in Eurasia. Nature 571, 500–504 (2019).CAS 
PubMed 
Article 

Google Scholar 
14.Hershkovitz, I. et al. The earliest modern humans outside Africa. Science (80-) 359, 456–459 (2018).CAS 
Article 
ADS 

Google Scholar 
15.Kolobova, K. A. et al. Archaeological evidence for two separate dispersals of Neanderthals into southern Siberia. Proc. Natl. Acad. Sci. U. S. A. 117, 2879–2885 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
16.Chen, F. et al. A late Middle Pleistocene Denisovan mandible from the Tibetan Plateau. Nature 569, 409–412 (2019).CAS 
PubMed 
Article 
ADS 

Google Scholar 
17.Douka, K. et al. Age estimates for hominin fossils and the onset of the Upper Palaeolithic at Denisova Cave. Nature 565, 640–644 (2019).CAS 
PubMed 
Article 
ADS 

Google Scholar 
18.Rizal, Y. et al. Last appearance of Homo erectus at Ngandong, Java, 117,000–108,000 years ago. Nature 577, 381–385 (2020).CAS 
PubMed 
Article 

Google Scholar 
19.Brown, P. et al. A new small-bodied hominin from the Late Pleistocene of Flores, Indonesia. Nature 431, 1055–1061 (2004).CAS 
PubMed 
Article 
ADS 

Google Scholar 
20.Détroit, F. et al. A new species of Homo from the Late Pleistocene of the Philippines. Nature 568, 181–186 (2019).PubMed 
Article 
ADS 
CAS 

Google Scholar 
21.Wedage, O. et al. Specialized rainforest hunting by Homo sapiens ~45,000 years ago. Nat. Commun. 10, 739 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
22.Athreya, S. Modern human emergence in South Asia: A review of the fossil and genetic evidence. in Emergence and Diversity of Modern Human Behaviour in Paleolithic Asia (eds. Kaifu, Y., Izuho, M., Goebel, T., Sato, H. & Ono, A.). 61–79. (Texas A&M University Press, 2015).23.Patnaik, R. et al. New geochronological, paleoclimatological, and archaeological data from the Narmada Valley hominin locality, central India. J. Hum. Evol. 56, 114–133 (2009).PubMed 
Article 

Google Scholar 
24.Pappu, S. et al. Early Pleistocene presence of Acheulian hominins in South India. Science 331, 1596–1599 (2011).CAS 
PubMed 
Article 
ADS 

Google Scholar 
25.Paddayya, K. et al. Recent findings on the Achuelian of the Hunsgi and Baichbal valleys, Karnataka, with special reference to the Isampur excavation and its dating. Curr. Sci. 83, 641–647 (2002).
Google Scholar 
26.Blinkhorn, J. & Petraglia, M. D. Environments and cultural change in the Indian subcontinent: Implications for the dispersal of Homo sapiens in the Late Pleistocene. Curr. Anthropol. 58, S463–S479 (2017).Article 

Google Scholar 
27.Benito-calvo, A., Barfod, D. N., Mchenry, L. J. & De, I. The geology and chronology of the Acheulean deposits in the Mieso area (East-Central Ethiopia). J. Hum. Evolut. 2014, 26–38 (2014).Article 

Google Scholar 
28.De la Torre, I., Mora, R., Arroyo, A. & Benito-calvo, A. Acheulean technological behaviour in the Middle Pleistocene landscape of Mieso (East-Central Ethiopia ). J. Hum. Evol. 76, 1–25 (2014).PubMed 
Article 

Google Scholar 
29.Shipton, C. et al. Acheulean technology and landscape use at Dawadmi, central Arabia. PLoS ONE 13, e200497 (2018).
Google Scholar 
30.Scerri, E. M. L. et al. The expansion of later Acheulean hominins into the Arabian Peninsula. Sci. Rep. 8, 1–9 (2018).CAS 
Article 

Google Scholar 
31.Brooks, A. S. et al. Long-distance stone transport and pigment use in the earliest Middle Stone Age. Science (80-) 360, 90–94 (2018).CAS 
Article 
ADS 

Google Scholar 
32.Valladas, H. et al. Dating the Lower to Middle Paleolithic transition in the Levant: A view from Misliya Cave, Mount Carmel, Israel. J. Hum. Evol. 65, 585–593 (2013).PubMed 
Article 

Google Scholar 
33.Carmignani, L., Moncel, M. H. E., Fernandes, P. & Wilson, L. Technological variability during the Early Middle Palaeolithic in Western Europe: Reduction systems and predetermined products at the Bau de l’Aubesier and Payre (South-East France). PLoS ONE 12, e178550 (2017).Article 
CAS 

Google Scholar 
34.Blinkhorn, J. et al. The first directly dated evidence for Palaeolithic occupation on the Indian coast at Sandhav, Kachchh. Quat. Sci. Rev. 224, 105975 (2019).Article 

Google Scholar 
35.Blinkhorn, J., Achyuthan, H., Petraglia, M. & Ditchfield, P. Middle Palaeolithic occupation in the Thar Desert during the Upper Pleistocene: the signature of a modern human exit out of Africa?. Quat. Sci. Rev. 77, 233–238 (2013).Article 
ADS 

Google Scholar 
36.Petraglia, M. et al. Middle Paleolithic assemblages from the Indian subcontinent before and after the Toba super-eruption. Science 317, 114–116 (2007).CAS 
PubMed 
Article 
ADS 

Google Scholar 
37.Petraglia, M. D., Ditchfield, P., Jones, S., Korisettar, R. & Pal, J. N. The Toba volcanic super-eruption, environmental change, and hominin occupation history in India over the last 140,000 years. Quat. Int. 258, 119–134 (2012).Article 

Google Scholar 
38.Clarkson, C. et al. Human occupation of northern India spans the Toba super-eruption ~74,000 years ago. Nat. Commun. 11, 1–10 (2020).Article 
CAS 

Google Scholar 
39.Akhilesh, K. et al. Early Middle Palaeolithic culture in India around 385–172 ka reframes out of Africa models. Nature 554, 97–101 (2018).CAS 
PubMed 
Article 
ADS 

Google Scholar 
40.Dennell, R. The Palaeolithic Settlement of Asia. (Cambridge University Press, 2009).41.Chauhan, P. R. Human evolution in the center of the old world: An updated review of the South Asian Paleolithic. in Pleistocene Archaeology—Migration, Technology, and Adaptation (eds. Ono, R. & Pawlik, A.). 94625. (IntechOpen, 2020).42.Roberts, P., Blinkhorn, J. & Petraglia, M. D. A transect of environmental variability across South Asia and its influence on Late Pleistocene human innovation and occupation. J. Quat. Sci. 33, 285–299 (2018).Article 

Google Scholar 
43.Goudie, A., Allchin, B. & Hegde, K. The former extensions of the Great Indian Sand desert. Geogr. J. 139, 243–257 (1973).Article 

Google Scholar 
44.Blinkhorn, J. The gateway to the oriental zone: Environmental change and palaeolithic behaviour in the Thar Desert. Quat. Int. 596, 79–92 (2021).Article 

Google Scholar 
45.Gaillard, C., Misra, V. N., Rajaguru, S. N., Raju, D. R. & Raghavan, H. Acheulian occupation at Singi Talav in the Thar Desert: A preliminary report on 1981 excavation. Bull. Deccan Coll. Res. Inst. 44, 141–152 (1985).
Google Scholar 
46.Raghavan, H., Gaillard, C. & Rajaguru, S. N. Genesis of Calcretes from the Calc-pan site of Singi Talav near a micromorphological approach. Geoarchaeology 6, 151–168 (1991).Article 

Google Scholar 
47.Misra, V. N. & Rajaguru, S. N. Palaeoenvironments and prehistory of the Thar Desert, Rajasthan, India. in South Asian Archaeology 1985. 296–320. (Scandinavian Institute of Asian Studies Occasional Papers, 1989).48.Misra, V. Geoarchaeology of the Thar desert, northwest India. Mem. Soc. India 210–230 (1995).49.Gaillard, C. Contribution à la Connaissance du Paléolithique Inférieur-Moyen en Inde. (Université de Provence, 1993).50.Gaillard, C., Mishra, S., Singh, M., Deo, S. & Abbas, R. Lower and Early Middle Pleistocene Acheulian in the Indian sub-continent. Quat. Int. 223–224, 234–241 (2010).Article 

Google Scholar 
51.Gaillard, C., Mishra, S., Singh, M., Deo, S. & Abbas, R. Reply to: “Comment on ‘lower and early Middle Pleistocene Acheulian in the Indian sub-continent’” by P. Chauhan. Quat. Int. 223–224, 260–264 (2010).Article 

Google Scholar 
52.Kailath, A. J. et al. Electron spin resonance characterization of calcretes from Thar desert for dating applications. Radiat. Meas. 32, 371–383 (2000).CAS 
Article 

Google Scholar 
53.Chauhan, P. R. Comment on ‘Lower and Early Middle Pleistocene Acheulian in the Indian sub-continent’ by Gaillard et al. (2009) (Quaternary International). Quat. Int. 223–224, 248–259 (2010).54.Alexandre, A., Meunier, J. D., Lézine, A. M., Vincens, A. & Schwartz, D. Phytoliths: Indicators of grassland dynamics during the late Holocene in intertropical Africa. Palaeogeogr. Palaeoclimatol. Palaeoecol. 136, 213–229 (1997).Article 

Google Scholar 
55.Diester-Haass, L., Schrader, H. J. & Thiede, J. Sedimentological and paleoclimatological investigations of two pelagic ooze cores off Cape Barbas, North-West Africa. Meteorol. Forshungergebnisse 16, 19–66 (1973).
Google Scholar 
56.Twiss, P. C. Predicted world distribution of C3 and C4 grass phytoliths. in Phytolith Systematics: Emerging Issues (eds. Rapp, G. R. & Mullholland, S. C.). 113–128. (Springer, 1992).57.Breecker, D. O., Sharp, Z. D. & McFadden, L. D. Seasonal bias in the formation and stable isotopic composition of pedogenic carbonate in modern soils from central New Mexico, USA. Bull. Geol. Soc. Am. 121, 630–640 (2009).CAS 
Article 

Google Scholar 
58.Szabo, B. J., McKinney, C., Dalbey, T. S. & Paddayya, K. On the age of the Acheulian culture of the Hunsgi-Baichbal Valleys, Peninsular India. Bull. Deccan Coll. Postgrad. Res. Inst. 50, 317–321 (1990).
Google Scholar 
59.Atkinson, T. J., Brown, P. J., Powar, N. J. & Kale, V. S. The Acheulian horizon at Chirki-Nevasa and its chronological implications. Bull. Deccan Coll. Postgrad. Res. Inst. 36, 3–14 (1990).
Google Scholar 
60.Deo, S. G., Joglekar, J. J. & Rajaguru, S. N. Geomorphic context of two acheulian sites in semi-arid peninsular India: Inferring palaeoenvironment and chronology. Quat. Int. 480, 166–177 (2018).Article 

Google Scholar 
61.Blinkhorn, J. A new synthesis of evidence for the Upper Pleistocene occupation of 16R Dune and its southern Asian context. Quat. Int. 300, 282–291 (2013).Article 

Google Scholar 
62.Baskaran, M., Maratheb, A. R., Rajagurub, S. N. & Somayajulu, B. L. K. Geochronology of Palaeolithic cultures in the Hiran Valley, Saurashtra, India. J. Archaeol. Sci. 13, 505–514 (1986).Article 

Google Scholar 
63.Jain, M., Tandon, S. K. & Bhatt, S. C. Late Quaternary stratigraphic development in the lower Luni, Mahi and Sabarmati river basins, western India. J. Earth Syst. Sci. 113, 453–471 (2004).Article 
ADS 

Google Scholar 
64.Juyal, N., Chamyal, L. S., Bhandari, S., Bhushan, R. & Singhvi, A. K. Continental record of the southwest monsoon during the last 130 ka : evidence from the southern margin of the Thar Desert, India. Quat. Sci. Rev. 25, 2632–2650 (2006).Article 
ADS 

Google Scholar 
65.Ajithprasad, P. Early Middle Palaeolithic: A transition phase between the Upper Acheulian and Middle Palaeolithic cultures in the Orsang Valley, Gujarat. Man Environ. 30, 1–11 (2005).
Google Scholar 
66.Bednarik, R. G., Kumar, G., Watchman, A. & Roberts, R. G. Preliminary results of the EIP Project. Rock Art Res. 22, 147–197 (2005).
Google Scholar 
67.Blinkhorn, J., Achyuthan, H., Ditchfield, P. & Petraglia, M. Palaeoenvironmental dynamics and Palaeolithic occupation at Katoati, Thar Desert. India. Quat. Res. (United States) 87, 298 (2017).CAS 

Google Scholar 
68.Gaillard, C. & Rajaguru, S. N. Revisiting the Acheulian site of Singi Talav at Didwana (Rajasthan) 35 years. in Rethinking the Past: A Tribute to Professor V. N. Misra (ed. Deo, S. G.). 25–39. https://doi.org/10.7765/9780719098451.00013 (Indian Society for Prehistoric and Quaternary Studies, 2017).69.d’Errico, F., Gaillard, C. & Misra, V. N. Collection of non-utilitarian objects by Homo erectus in India. in Hominidae: Proceedings of the 2nd International Congress of Human Paleontology (ed. Giacobini, G.). 237–239. (Jaca Book, 1989).70.Key, A. J. M., Jarić, I. & Roberts, D. L. Modelling the end of the Acheulean at global and continental levels suggests widespread persistence into the Middle Palaeolithic. Hum. Soc. Sci. Commun. 8, 1–12 (2021).Article 

Google Scholar 
71.Wasson, R. J., Smith, G. I. & Agrawala, D. P. Late Quaternary sediments, minerals and inferred geochemical history of Didwana Lake, Thar Dessert, India. Palaeogeogr. Palaeoclimatol. Palaeoecol. 46, 345–372 (1984).CAS 
Article 

Google Scholar 
72.Jakher, G. R., Bhargava, S. C. & Sinha, R. K. Comparative limnology of Sambhar and Didwana lakes (Rajasthan, NW India). Saline Lakes 67, 245–256 (1990).Article 

Google Scholar 
73.Sinha, R. et al. Late Quaternary palaeoclimatic reconstruction from the lacustrine sediments of the Sambhar playa core, Thar Desert margin, India. Palaeogeogr. Palaeoclimatol. Palaeoecol. 233, 252–270 (2006).Article 

Google Scholar 
74.Sinha, R. & Raymahashay, B. C. Evaporite mineralogy and geochemical evolution of the Sambhar Salt Lake, Rajasthan, India. Sediment. Geol. 166, 59–71 (2004).CAS 
Article 
ADS 

Google Scholar 
75.Deotare, B. C. et al. Palaeoenvironmental history of Bap-Malar and Kanod playas of western Rajasthan, Thar desert. J. Earth Syst. Sci. 113, 403–425 (2004).Article 
ADS 

Google Scholar 
76.Roy, P. D., Sinha, R., Smykatz-Kloss, W., Singhvi, A. K. & Nagar, Y. C. Playas of the Thar Desert: Mineralogical and geochemical archives of Late Holocene climates. Asian J. Earth Sci. 1, 43–61 (2008).Article 
CAS 

Google Scholar 
77.Achyuthan, H., Kar, A. & Eastoe, C. Late Quaternary-Holocene lake-level changes in the eastern margin of the Thar Desert, India. J. Paleolimnol. 38, 493–507 (2007).Article 
ADS 

Google Scholar 
78.Dhir, R. P. et al. Multiple episodes of aggradation and calcrete formation in Late Quaternary aeolian sands, Central Thar Desert, Rajasthan, India. J. Asian Earth Sci. 37, 10–16 (2010).Article 
ADS 

Google Scholar 
79.Juyal, N., Chamyal, L., Bhandari, S., Bhushan, R. & Singhvi, A. Continental record of the southwest monsoon during the last 130 ka: Evidence from the southern margin of the Thar Desert, India. Quat. Sci. Rev. 25, 2632–2650 (2006).Article 
ADS 

Google Scholar 
80.Giosan, L. et al. Fluvial landscapes of the Harappan civilization. Proc. Natl. Acad. Sci. 109, E1688–E1694 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
81.Achyuthan, H., Quade, J., Roe, L. & Placzek, C. Stable isotopic composition of pedogenic carbonates from the eastern margin of the Thar Desert, Rajasthan, India. Quat. Int. 162, 50–60 (2007).Article 

Google Scholar 
82.Sharma, K., Bhatt, N., Shukla, A. D., Cheong, D. K. & Singhvi, A. K. Optical dating of late Quaternary carbonate sequences of Saurashtra, western India. Quat. Res. (United States) 87, 133–150 (2017).CAS 

Google Scholar 
83.Blinkhorn, J., Achyuthan, H., Jaiswal, M. & Singh, A. K. The first dated evidence for Middle-Late Pleistocene fluvial activity in the central Thar Desert. Quat. Sci. Rev. 250, 106656 (2020).Article 

Google Scholar 
84.Singhvi, A. K. et al. A ~200 ka record of climatic change and dune activity in the Thar Desert, India. Quat. Sci. Rev. 29, 3095–3105 (2010).Article 
ADS 

Google Scholar 
85.Blinkhorn, J., Achyuthan, H., Ditchfield, P. & Petraglia, M. Palaeoenvironmental dynamics and Palaeolithic occupation at Katoati, Thar Desert, India. Quat. Res. 87, 298–313 (2017).CAS 
Article 

Google Scholar 
86.Shipton, C. Hierarchical organization in the Acheulean to Middle Palaeolithic transition at Bhimbetka, India. Camb. Archaeol. J. 26, 601–618 (2016).Article 

Google Scholar 
87.Shipton, C. et al. Generativity, hierarchical action and recursion in the technology of the Acheulean to Middle Palaeolithic transition: A perspective from Patpara, the Son Valley, India. J. Hum. Evol. 65, 93–108 (2013).CAS 
PubMed 
Article 

Google Scholar 
88.Potts, R. et al. Environmental dynamics during the onset of the Middle Stone Age in eastern Africa. Science (80-) 360, 86–90 (2018).CAS 
Article 
ADS 

Google Scholar 
89.Tryon, C. A. & Mcbrearty, S. Tephrostratigraphy of the bedded tuff member (Kapthurin Formation, Kenya) and the nature of archaeological change in the later middle Pleistocene. Quatern. Res. 65, 492–507 (2006).Article 
ADS 

Google Scholar 
90.Tryon, C. A., McBrearty, S. & Texier, P.-J. Levallois lithic technology from the Kapthurin Formation, Kenya: Acheulian origin and Middle Stone Age diversity. African Archaeol. Rev. 22, 199–229 (2006).Article 

Google Scholar 
91.Zaidner, Y. & Weinstein-Evron, M. The emergence of the Levallois technology in the Levant: A view from the Early Middle Paleolithic site of Misliya Cave. Israel. J. Hum. Evol. 2020, 102785 (2020).Article 

Google Scholar 
92.Jacobs, G. S. et al. Multiple deeply divergent Denisovan ancestries in papuans. Cell 177, 1010-1021.e32 (2019).CAS 
PubMed 
Article 

Google Scholar 
93.Browning, S. R., Browning, B. L., Zhou, Y., Tucci, S. & Akey, J. M. Analysis of human sequence data reveals two pulses of archaic Denisovan admixture. Cell 173, 53-61.e9 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
94.Sankararaman, S. et al. The combined landscape of Denisovan and Neanderthal ancestry in present-day humans report the combined landscape of Denisovan and Neanderthal ancestry in present-day humans. Curr. Biol. 26, 1–7 (2016).Article 
CAS 

Google Scholar 
95.Holt, B. G. et al. An update of Wallace’s zoogeographic regions of the world. Science 339, 74–78 (2013).CAS 
PubMed 
Article 
ADS 

Google Scholar 
96.Blott, S. J. & Pye, K. Technical communication Gradistat : A grain size distribution and statistics package for the analysis of unconsolidated sediments. Earth Surf. Process. Landforms 26, 1237–1248 (2001).Article 
ADS 

Google Scholar 
97.Keeling, P. S. Some experiments on the low-temperature removal of carbonaceous material from clays. Clay Miner. Bull. 62, 155–158 (1962).Article 
ADS 

Google Scholar 
98.Ball, D. F. Loss-on-ignition as an estimate of organic matter and organic carbon in non-calcareous soils. J. Soil Sci. 15, 84–92 (1964).CAS 
Article 

Google Scholar 
99.Dean, W. E. Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: comparison with other methods. J. Sediment. Res. 44, 242–248 (1974).CAS 

Google Scholar 
100.Bengtsson, L. & Enell, M. Chemical analysis. in Handbook of Holocene Paleoecology and Paleohydrology (eds. Bengtsson, L., Enell, M. & Berglund, B. E.) 423–454 (Wiley, 1986).101.Juggins, S. rioja: Anlaysis of Quaternary science data. R Doc. 1–58 (2015).102.Team, R. D. C. R: A Language and Environment for Statistical Computing. http://www.r-project.org (2020).103.Parker, A. An index of weathering for silicate rocks. Geol. Mag. 107, 501–504 (1970).CAS 
Article 
ADS 

Google Scholar 
104.Nesbitt, H. W. & Young, G. M. Formation and diagenesis of weathering profiles. J. Geol. 2, 129–147 (1989).Article 
ADS 

Google Scholar 
105.Piperno, D. R. Phytolith Analysis: An Archaeological and Geological Perspective. (Academic Press, 1988).106.Twiss, P. A Curmudgeon’s view of grass phytolithology. in Phytoliths: Applications in Earth Sciences and Human History (ed. Meunier, J. D.) 7–25 (A.A. Balkema Publishers, 2001).107.Eksambekar, S. Contribution of the Study of Phytoliths to Bioarchaeology. (Deccan College PGRI, 2002).108.Durcan, J., King, E. G. & Duller, G. A. T. DRAC : Dose rate and age calculator for trapped charge dating DRAC : Dose rate and age calculator for trapped charge dating. Quat. Geochronol. 28, 54–61 (2015).Article 

Google Scholar 
109.Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).Article 

Google Scholar 
110.Kathayat, G. et al. Indian monsoon variability on millennial-orbital timescales. Sci. Rep. 6, 4–10 (2016).Article 
CAS 

Google Scholar 
111.Caley, T. et al. New Arabian Sea records help decipher orbital timing of Indo-Asian monsoon. Earth Planet. Sci. Lett. 308, 433–444 (2011).CAS 
Article 
ADS 

Google Scholar 
112.Clemens, S. C. & Prell, W. L. A 350,000 year summer-monsoon multi-proxy stack from the Owen Ridge, Northern Arabian Sea. Mar. Geol. 201, 35–51 (2003).Article 
ADS 

Google Scholar 
113.Bolton, C. T. et al. A 500,000 year record of Indian summer monsoon dynamics recorded by eastern equatorial Indian Ocean upper water-column structure. Quat. Sci. Rev. 77, 167–180 (2013).Article 
ADS 

Google Scholar  More

1.Schoock, M. Tractatus de turffis ceu cespitibus bituminosis: quo multa, ab aliis hactenus aut neglecta, aut minus diligenter examinata, accuratius aliquanto excutiuntur. [A treatise on peats as bituminous sods] (Cöllen, Groningen, 1658).2.Tanneberger, F. et al. Towards net zero CO2 in 2050: An emission reduction pathway for organic soils in Germany. Mires Peat 27, 1–17, https://doi.org/10.19189/MaP.2020.SNPG.StA.1951 (2021).Article 

Google Scholar 
3.Klimkowska, A. et al. Are we restoring functional fens? – The outcomes of restoration projects in fens re-analysed with plant functional traits. PLOS One 14, e0215645, https://doi.org/10.1371/journal.pone.0215645 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
4.Emsens, W.-J. et al. Recovery of fen peatland microbiomes and predicted functional profiles after rewetting. ISME J. 14, 1701–1712, https://doi.org/10.1038/s41396-020-0639-x (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
5.Lamers, L. P. M. et al. Ecological restoration of rich fens in Europe and North America: from trial and error to an evidence-based approach. Biol. Rev. 90, 182–203, https://doi.org/10.1111/brv.12102 (2015).PubMed 
Article 

Google Scholar 
6.Findorff, J. C. Beiträge und Fragmente zu einem Moorkatechismus/von Jürgen Christian Findorff und Anmerkungen von F. Brüne, K. Lilienthal und F. Overbeck. [Contributions and fragments to a peatland catechism from Jürgen Christian Findorff with remarks from F. Brüne, K. Lilienthal and F. Overbeck] (Stalling, Oldenburg, 1764, published 1937).7.Kirpotin, S. N. et al. Great Vasyugan Mire: How the world’s largest peatland helps addressing the world’s largest problems. Ambio 71, 618, https://doi.org/10.1007/s13280-021-01520-2 (2021).Article 

Google Scholar 
8.Scharlemann, J. P. W., Tanner, E. V. J., Hiederer, R. & Kapos, V. Global soil carbon: understanding and managing the largest terrestrial carbon pool. Carbon Manag. 5, 81–91, https://doi.org/10.4155/cmt.13.77 (2014).CAS 
Article 

Google Scholar 
9.Xu, J., Morris, P. J., Liu, J. & Holden, J. Hotspots of peatland-derived potable water use identified by global analysis. Nat. Sustain. 1, 246–253, https://doi.org/10.1038/s41893-018-0064-6 (2018).CAS 
Article 

Google Scholar 
10.Joosten, H., Tanneberger, F. & Moen, A. (eds). Mires and peatlands of Europe. Status, distribution and conservation (Schweizerbart Science Publishers, Stuttgart, 2017).11.Bonn, A., Allott, T., Evans, M., Joosten, H. & Stoneman, R. (eds). Peatland restoration and ecosystem services (Cambridge University Press, Cambridge, 2016).12.Leifeld, J. & Menichetti, L. The underappreciated potential of peatlands in global climate change mitigation strategies. Nat. Commun. 9, 1071, https://doi.org/10.1038/s41467-018-03406-6 (2018).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
13.Joosten, H. & Clarke, D. Wise use of mires and peatlands. Background and principles including a framework for decision-making (Internat. Mire Conservation Group, Totnes, 2002).14.Günther, A. et al. Prompt rewetting of drained peatlands reduces climate warming despite methane emissions. Nat. Commun. 11, 1644, https://doi.org/10.1038/s41467-020-15499-z (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
15.Erkens, G., van der Meulen, M. J. & Middelkoop, H. Double trouble: subsidence and CO2 respiration due to 1,000 years of Dutch coastal peatlands cultivation. Hydrogeol. J. 24, 551–568, https://doi.org/10.1007/s10040-016-1380-4 (2016).ADS 
CAS 
Article 

Google Scholar 
16.Tiemeyer, B. & Kahle, P. Nitrogen and dissolved organic carbon (DOC) losses from an artificially drained grassland on organic soils. Biogeosciences 11, 4123–4137, https://doi.org/10.5194/bg-11-4123-2014 (2014).ADS 
CAS 
Article 

Google Scholar 
17.Minayeva, T., Bragg, O. & Sirin, A. Peatland biodiversity and its restoration. In Peatland restoration and ecosystem services, (eds Bonn, A., Allott, T., Evans, M., Joosten, H. & Stoneman, R.) pp. 44–62 (Cambridge University Press, Cambridge, 2016).18.Janssen, J. A. M. et al. European red list of habitats (Publications Office of the European Union, 2016).19.Mrotzek, A., Michaelis, D., Günther, A., Wrage-Mönnig, N. & Couwenberg, J. Mass balances of a drained and a rewetted peatland: on former losses and recent gains. Soil Syst. 4, 16, https://doi.org/10.3390/soilsystems4010016 (2020).CAS 
Article 

Google Scholar 
20.Zerbe, S. et al. Ecosystem service restoration after 10 years of rewetting peatlands in NE Germany. Environ. Manag. 51, 1194–1209, https://doi.org/10.1007/s00267-013-0048-2 (2013).ADS 
Article 

Google Scholar 
21.Leifeld, J., Wüst-Galley, C. & Page, S. Intact and managed peatland soils as a source and sink of GHGs from 1850 to 2100. Nat. Clim. Change 9, 945–947, https://doi.org/10.1038/s41558-019-0615-5 (2019).ADS 
CAS 
Article 

Google Scholar 
22.IPCC. Global warming of 1.5 °C (IPCC, Geneva, Switzerland, 2018).23.Liu, H. & Lennartz, B. Hydraulic properties of peat soils along a bulk density gradient-A meta study. Hydrol. Process 33, 101–114, https://doi.org/10.1002/hyp.13314 (2019).ADS 
Article 

Google Scholar 
24.Franz, D., Koebsch, F., Larmanou, E., Augustin, J. & Sachs, T. High net CO2 and CH4 release at a eutrophic shallow lake on a formerly drained fen. Biogeosciences 13, 3051–3070, https://doi.org/10.5194/bg-13-3051-2016 (2016).ADS 
CAS 
Article 

Google Scholar 
25.Hemes, K. S., Chamberlain, S. D., Eichelmann, E., Knox, S. H. & Baldocchi, D. D. A biogeochemical compromise: the high methane cost of sequestering carbon in restored wetlands. Geophys. Res. Lett. 45, 6081–6091, https://doi.org/10.1029/2018GL077747 (2018).ADS 
CAS 
Article 

Google Scholar 
26.Zak, D. & Gelbrecht, J. The mobilisation of phosphorus, organic carbon and ammonium in the initial stage of fen rewetting (a case study from NE Germany). Biogeochemistry 85, 141–151, https://doi.org/10.1007/s10533-007-9122-2 (2007).CAS 
Article 

Google Scholar 
27.Emsens, W.-J., Aggenbach, C. J. S., Smolders, A. J. P., Zak, D. & van Diggelen, R. Restoration of endangered fen communities: the ambiguity of iron-phosphorus binding and phosphorus limitation. J. Appl. Ecol. 54, 1755–1764, https://doi.org/10.1111/1365-2664.12915 (2017).CAS 
Article 

Google Scholar 
28.Chytrý, M. et al. EUNIS Habitat Classification: Expert system, characteristic species combinations and distribution maps of European habitats. Appl. Veg. Sci. 23, 648–675, https://doi.org/10.1111/avsc.12519 (2020).Article 

Google Scholar 
29.Zak, D. et al. Unraveling the importance of polyphenols for microbial carbon mineralization in rewetted riparian peatlands. Front. Environ. Sci. 7; https://doi.org/10.3389/fenvs.2019.00147 (2019).30.Henneberg, A., Sorrell, B. K. & Brix, H. Internal methane transport through Juncus effusus: experimental manipulation of morphological barriers to test above- and below-ground diffusion limitation. N. Phytol. 196, 799–806, https://doi.org/10.1111/j.1469-8137.2012.04303.x (2012).CAS 
Article 

Google Scholar 
31.Agethen, S., Sander, M., Waldemer, C. & Knorr, K.-H. Plant rhizosphere oxidation reduces methane production and emission in rewetted peatlands. Soil Biol. Biochem 125, 125–135, https://doi.org/10.1016/j.soilbio.2018.07.006 (2018).CAS 
Article 

Google Scholar 
32.Vitt, D. H. Peatlands: ecosystems dominated by bryophytes. In Bryophyte biology, (eds Shaw, A. J. & Goffinet, B.), pp. 312–343 (Cambridge University Press, Cambridge, 2002).33.Wilson, D. et al. Greenhouse gas emission factors associated with rewetting of organic soils. Mires Peat, 1–28; https://doi.org/10.19189/MaP.2016.OMB.222 (2016).34.Kotowski, W., Thörig, W., van Diggelen, R. & Wassen, M. J. Competition as a factor structuring species zonation in riparian fens – a transplantation experiment. Appl. Veg. Sci. 9, 231–240, https://doi.org/10.1111/j.1654-109X.2006.tb00672.x (2006).Article 

Google Scholar 
35.Frantz, D. FORCE—Landsat+ Sentinel-2 Analysis Ready Data and Beyond. Remote Sens. 11, 1124, https://doi.org/10.3390/rs11091124 (2019).ADS 
Article 

Google Scholar 
36.Körner, C., Stöcklin, J., Reuther-Thiébaud, L. & Pelaez-Riedl, S. Small differences in arrival time influence composition and productivity of plant communities. N. Phytol. 177, 698–705, https://doi.org/10.1111/j.1469-8137.2007.02287.x (2008).Article 

Google Scholar 
37.Fritz, C., Campbell, D. I. & Schipper, L. A. Oscillating peat surface levels in a restiad peatland, New Zealand-magnitude and spatiotemporal variability. Hydrological Process. 22, 3264–3274, https://doi.org/10.1002/hyp.6912 (2008).ADS 
Article 

Google Scholar 
38.Loisel, J. & Yu, Z. Surface vegetation patterning controls carbon accumulation in peatlands. Geophys. Res. Lett. 40, 5508–5513, https://doi.org/10.1002/grl.50744 (2013).ADS 
CAS 
Article 

Google Scholar 
39.Koebsch, F. et al. The impact of occasional drought periods on vegetation spread and greenhouse gas exchange in rewetted fens. Proc. R. Soc. Lond. B Biol. Sci. 375, 20190685, https://doi.org/10.1098/rstb.2019.0685 (2020).CAS 
Article 

Google Scholar 
40.Fenner, N. & Freeman, C. Drought-induced carbon loss in peatlands. Nat. Geosci. 4, 895–900, https://doi.org/10.1038/ngeo1323 (2011).ADS 
CAS 
Article 

Google Scholar 
41.Hobbs, R. J., Higgs, E. S. & Hall, C. M. Defining novel ecosystems. In Novel Ecosystems. Intervening in the New Ecological World Order (eds Hobbs, R. J., Higgs, E. S. & Hall, C.) 1st ed., Vol. 24, pp. 58–60 (Wiley-Blackwell, s.l., 2013).42.Bonnett, S. A. F., Ross, S., Linstead, C. & Maltby, E. A review of techniques for monitoring the success of peatland restoration. Natural England Commissioned Reports, Number 086 and Natural England Technical Information Note TIN097. http://publications.naturalengland.org.uk/publication/46013 (2011).43.Schwieger, S. et al. Wetter is better: Rewetting of minerotrophic peatlands increases plant production and moves them towards carbon sinks in a dry year. Ecosystems 45, 279, https://doi.org/10.1007/s10021-020-00570-z (2020).CAS 
Article 

Google Scholar 
44.Beyer, F. et al. Drought years in peatland rewetting: rapid vegetation succession can maintain the net CO2 sink function. Biogeosciences 18, 917–935, https://doi.org/10.5194/bg-18-917-2021 (2021).ADS 
CAS 
Article 

Google Scholar 
45.Xu, H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 27, 3025–3033, https://doi.org/10.1080/01431160600589179 (2006).ADS 
Article 

Google Scholar 
46.Chytrý, M. et al. European Vegetation Archive (EVA): an integrated database of European vegetation plots. Appl. Veg. Sci. 19, 173–180, https://doi.org/10.1111/avsc.12191 (2016).Article 

Google Scholar 
47.Euro+Med 2006+. Euro+Med PlantBase – the information resource for Euro-Mediterranean plant diversity. https://www.europlusmed.org. Accessed 25-Feb-2021.48.Hodgetts, N. G. et al. An annotated checklist of bryophytes of Europe, Macaronesia and Cyprus. J. Bryol. 42, 1–116, https://doi.org/10.1080/03736687.2019.1694329 (2020).Article 

Google Scholar 
49.Jansen, F. & Dengler, J. Plant names in vegetation databases – a neglected source of bias. J. Veg. Sci. 21, 1179–1186, https://doi.org/10.1111/j.1654-1103.2010.01209.x (2010).Article 

Google Scholar 
50.Dufrene, M. & Legendre, P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol. Monogr. 67, 345–366 (1997).
Google Scholar 
51.Jurasinski, G. & Kreyling, J. Upward shift of alpine plants increases floristic similarity of mountain summits. J. Veg. Sci. 18, 711–718 (2007).Article 

Google Scholar 
52.R Core Team. R: a language and environment for statistical computing. R version 4.0.2 (R Foundation for Statistical Computing. URL http://www.R-project.org, Vienna, Austria, 2020).53.Tanneberger, F. et al. The peatland map of Europe. Mires Peat, 1–17; https://doi.org/10.19189/MaP.2016.OMB.264 (2017). More

1.Krebs, C. J. Of lemmings and snowshoe hares: the ecology of northern Canada. Proc. R. Soc. B: Biol. Sci. 278, 481–489 (2011).Article 

Google Scholar 
2.Oli, M. K. Population cycles in voles and lemmings: state of the science and future directions. Mammal Rev. 49, 226–239 (2019).Article 

Google Scholar 
3.Tallian, A. et al. Competition between apex predators? Brown bears decrease wolf kill rate on two continents. Proc. R. Soc. B: Biol. Sci. 284, 20162368 (2017).Article 

Google Scholar 
4.Henden, J.-A., Ims, R. A., Yoccoz, N. G., Hellström, P. & Angerbjörn, A. Strength of asymmetric competition between predators in food webs ruled by fluctuating prey: the case of foxes in tundra. Oikos 119, 27–34 (2010).Article 

Google Scholar 
5.Sih, A., Englund, G. & Wooster, D. Emergent impacts of multiple predators on prey. Trends Ecol. Evol. 13, 350–355 (1998).CAS 
PubMed 
Article 

Google Scholar 
6.Griffen, B. D. Detecting emergent effects of multiple predator species. Oecologia 148, 702–709 (2006).ADS 
PubMed 
Article 

Google Scholar 
7.Jones, H. P. et al. Severity of the effects of invasive rats on seabirds: a global review. Conserv. Biol. 22, 16–26 (2008).PubMed 
Article 

Google Scholar 
8.Dowding, J. E. & Murphy, E. C. The impact of predation by introduced mammals on endemic shorebirds in New Zealand: a conservation perspective. Biol. Conserv. 99, 47–64 (2001).Article 

Google Scholar 
9.Pascal, M., Lorvelec, O., Bretagnolle, V. & Culioli, J.-M. Improving the breeding success of a colonial seabird: a cost-benefit comparison of the eradication and control of its rat predator. Endangered Species Res. 4, 267–276 (2008).Article 

Google Scholar 
10.Ruffino, L. et al. Invasive rats and seabirds after 2,000 years of an unwanted coexistence on Mediterranean islands. Biol. Invasions 11, 1631–1651 (2009).Article 

Google Scholar 
11.Harper, G. A. Detecting predation of a burrow-nesting seabird by two introduced predators, using stable isotopes, dietary analysis and experimental removals. Wildl. Res. 34, 443–453 (2007).Article 

Google Scholar 
12.Gaze, P. The response of a colony of sooty shearwater (Puffinus griseus) and flesh-footed shearwater (P. carneipes) to the cessation of harvesting and the eradication of Norway rats (Rattus norvegicus). N. Z. J. Zool. 27, 375–379 (2000).Article 

Google Scholar 
13.Dowding, J. E. & Murphy, E. C. Ecology of ship rats (Rattus rattus) in a kauri (Agathis australis) forest in Northland, New Zealand. N. Z. J. Ecol., 19–27 (1994).14.Major, H., Bond, A., Jones, I. & Eggleston, C. Stability of a seabird population in the presence of an introduced predator. Avian Conserv. Ecol. 8 (2013).15.Bond, A. L. & Eggleston, C. J. Application of a non-invasive indexing method for introduced Norway rats (Rattus norvegicus) in the Aleutian Islands, Alaska. Biodiversity Conserv. 24, 2551–2563 (2015).Article 

Google Scholar 
16.Major, H. L., Jones, I. L., Byrd, G. V. & Williams, J. C. Assessing the effects of introduced Norway rats (Rattus norvegicus) on survival and productivity of Least Auklets (Aethia pusilla). Auk 123, 681–694 (2006).Article 

Google Scholar 
17.Eaton, M. et al. Birds of conservation concern 4: the population status of birds in the UK, Channel Islands and Isle of Man. British Birds 108, 708–746 (2015).
Google Scholar 
18.JNCC. Manx shearwater (Puffinus puffinus). https://jncc.gov.uk/our-work/manx-shearwater-puffinus-puffinus/ (2021).19.Bell, E. et al. The ground-based eradication of Norway rats (Rattus norvegicus) from the Isle of Canna, Inner Hebrides, Scotland. In Island Invasives: Eradication and Management (eds Veitch C. R. et al.) 269–274 (IUCN, Gland, Switzerland, 2011).20.Luxmoore, R., Swann, R. & Bell, E. Canna seabird recovery project: 10 years on. In Island invasives: scaling up to meet the challenge (eds Veitch C. R. et al.) 576–579 (IUCN, Gland, Switzerland, 2019).21.Mavor, R. A., Parsons, M., Heubeck, M. & Schmitt, S. Seabird numbers and breeding success in Britain and Ireland, 2004. (Peterborough, Joint Nature Conservation Committee (UK Nature Conservation, No. 29), 2005).22.Ratcliffe, N., Mitchell, I., Varnham, K., Verboven, N. & Higson, P. How to prioritize rat management for the benefit of petrels: a case study of the UK, Channel Islands and Isle of Man. Ibis 151, 699–708 (2009).Article 

Google Scholar 
23.Duron, Q., Shiels, A. B. & Vidal, E. Control of invasive rats on islands and priorities for future action. Conserv. Biol. 31, 761–771 (2017).PubMed 
Article 

Google Scholar 
24.Watt, L. & Sargent, I. Management Plan for Rum National Nature Reserve 2016–2026 (Scottish Natural Heritage, 2018).25.Lloyd, C., Tasker, M. L. & Partridge, K. The status of seabirds in Britain and Ireland. (A&C Black, 2010).26.Harris, M. Breeding biology of the Manx Shearwater Puffinus puffinus. Ibis 108, 17–33 (1966).Article 

Google Scholar 
27.Thompson, K. R. The ecology of the Manx Shearwater Puffinus puffinus on Rhum, west Scotland, Ph.D. thesis, University of Glasgow (1987).28.Perrins, C., Harris, M. & Britton, C. Survival of Manx shearwaters Puffinus puffinus. Ibis 115, 535–548 (1973).Article 

Google Scholar 
29.Clutton-Brock, T. H. & Ball, M. Rhum: the natural history of an island. (Edinburgh University Press, 1987).30.Virtanen, R., Edwards, G. R. & Crawley, M. J. Red deer management and vegetation on the Isle of Rum. J. Appl. Ecol. 39, 572–583 (2002).31.Corbet, G. Origin of the British insular races of small mammals and of the ‘Lusitanian’ fauna. Nature 191, 1037–1040 (1961).ADS 
Article 

Google Scholar 
32.Berry, R. History in the evolution of Apodemus sylvaticus (Mammalia) at one edge of its range. J. Zool. 159, 311–328 (1969).Article 

Google Scholar 
33.Meehan, A. P. Rats and mice: Their biology and control. (Rentokil Ltd, 1984).34.Carlisle, S. B. Rum rats: ecology and behaviour, Ph.D. thesis, Anglia Ruskin University (2019).35.Berry, R., Evans, I. & Sennitt, B. The relationships and ecology of Apodemus sylvaticus from the Small Isles of the Inner Hebrides, Scotland. J. Zool. 152, 333–346 (1967).Article 

Google Scholar 
36.SNH. Rum SSSI Citation, https://apps.snh.gov.uk/sitelink-api/v1/sites/1396/documents/1 (2010).37.Ambagis, J. A comparison of census and monitoring techniques for Leach’s storm petrel. Waterbirds 27, 211–215 (2004).Article 

Google Scholar 
38.Bicknell, T. W., Reid, J. B. & Votier, S. C. Probable predation of Leach’s Storm-petrel Oceanodroma leucorhoa eggs by St Kilda Field Mice Apodemus sylvaticus hirtensis. Bird Study 56, 419–422 (2009).Article 

Google Scholar 
39.Rodríguez, B., De León, L., Martín, A., Alonso, J. & Nogales, M. Status and distribution of breeding seabirds in the northern islets of Lanzarote, Canary Islands. Atlantic Seabirds 5, 41–56 (2003).
Google Scholar 
40.Seto, N. W. H. & Conant, S. The effects of rat (Rattus rattus) predation on the reproductive success of the Bonin Petrel (Pterodroma hypoleuca) on Midway Atoll. Colon. Waterbirds 19, 171–185 (1996).Article 

Google Scholar 
41.Walsh, P. et al. Seabird monitoring handbook for Britain and Ireland: a compilation of methods for survey and monitoring of breeding seabirds. (JNCC/RSPB/ITE/Seabird Group, 1995).42.Quy, R. J., Cowan, D. P. & Swinney, T. Tracking as an activity index to measure gross changes in Norway rat populations. Wildl. Soc. Bull. 21, 122–127 (1993).
Google Scholar 
43.Martin, A. & Richardson, M. Rodent eradication scaled up: clearing rats and mice from South Georgia. Oryx 53, 27–35 (2019).Article 

Google Scholar 
44.European Chemicals Agency. Technical Notes for Guidance on Product Evaluation Appendices to Chapter 7 Product Type 14 Efficacy Evaluation of Rodenticidal Biocidal Products. https://echa.europa.eu/documents/10162/16960215/bpd_guid_revised_appendix_chapter_7_pt14_2009_en.pdf (2009).45.European Chemicals Agency. Transitional Guidance on Efficacy Assessment for Product Type 14 Rodenticides, Helsinki, Finland. https://echa.europa.eu/documents/10162/23492134/tg_efficacy_pt14_superseded_en.pdf/f82f63de-6ca6-abed-53fe-8b7f11b53493 (2016).46.CRRU. CRRU UK Code of Best Practice. Campaign for Responsible Rodenticide Use, Leeds, UK. https://bpca.org.uk/write/MediaUploads/Documents/Codes%20of%20Best%20Practice/COBP-CRRU-Rodent-Control-and-Safe-Use-of-Rodenticides-2015.PDF (2015).47.Pelz, H.-J. et al. The genetic basis of resistance to anticoagulants in rodents. Genetics 170, 1839–1847 (2005).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
48.Rost, S. et al. Novel mutations in the VKORC1 gene of wild rats and mice – a response to 50 years of selection pressure by warfarin? BMC Genet. 10, 1–9 (2009).Article 
CAS 

Google Scholar 
49.Berny, P., Esther, A., Jacob, J. & Prescott, C. Development of resistance to anticoagulant rodenticides in rodents. In Anticoagulant rodenticides and wildlife. Emerging Topics in Ecotoxicology (Principles, Approaches and Perspectives), Vol. 5 (eds van den Brink, N. et al.) 259–286 (Springer, 2018).50.Fisher, P. Review of house mouse (Mus musculus) susceptibility to anticoagulant poisons. (Department of Conservation, Wellington, 2005).51.Lambert, M. S., Quy, R. J., Smith, R. H. & Cowan, D. P. The effect of habitat management on home-range size and survival of rural Norway rat populations. J. Appl. Ecol. 45, 1753–1761. https://doi.org/10.1111/j.1365-2664.2008.01543.x (2008).Article 

Google Scholar 
52.Taylor, K. D. & Quy, R. J. Long distance movements of a Common rat (Rattus norvegicus ) revealed by radio-tracking. Mammalia 42, 63–71 (1978).Article 

Google Scholar 
53.Schall, R. Estimation in generalized linear models with random effects. Biometrika 78, 719–727 (1991).MATH 
Article 

Google Scholar 
54.Engeman, R. M. Indexing principles and a widely applicable paradigm for indexing animal populations. Wildl. Res. 32, 203–210 (2005).Article 

Google Scholar 
55.Saad, S. M., Sanderson, R., Robertson, P. & Lambert, M. Effects of supplementary feed for game birds on activity of brown rats Rattus norvegicus on arable farms. Mammal Res. 66, 163–171 (2021).56.O’Hara, R. B. & Kotze, D. J. Do not log-transform count data. Methods Ecol. Evol. 1, 118–122 (2010).Article 

Google Scholar 
57.Watts, C. H. The foods eaten by wood mice (Apodemus sylvaticus) and bank voles (Clethrionomys glareolus) in Wytham Woods, Berkshire. J. Anim. Ecol. 37, 25–41 (1968).58.Zubaid, A. & Gorman, M. The diet of wood mice Apodemus sylvaticus living in a sand dune habitat in north-east Scotland. J. Zool. 225, 227–232 (1991).Article 

Google Scholar 
59.Bicknell, A. W. et al. Stable isotopes reveal the importance of seabirds and marine foods in the diet of St Kilda field mice. Sci. Rep. 10, 1–12 (2020).Article 
CAS 

Google Scholar 
60.Haarsma, A.-J. & Kaal, R. Predation of wood mice (Apodemus sylvaticus) on hibernating bats. Popul. Ecol. 58, 567–576 (2016).Article 

Google Scholar 
61.Taylor, K. D., Quy, R. J. & Gurnell, J. Comparison of three methods for estimating the numbers of common rats (Rattus norvegicus). Mammalia 45, 403–413 (1981).Article 

Google Scholar 
62.Lambert, M. et al. Validating activity indices from camera traps for commensal rodents and other wildlife in and around farm buildings. Pest Manag. Sci. 74, 70–77 (2018).CAS 
PubMed 
Article 

Google Scholar 
63.Pankhurst, S., Conlan, H., Rogers, M., Smith, A. Survey of Wood Mice (Apodemus sylvaticus) on the Isle of Rum, Inner Hebrides. (http://insight.cumbria.ac.uk/id/eprint/2191/1/Pankhurst_SurveyOfWoodmice.pdf, 2010).64.Harris, D. B. & Macdonald, D. W. Interference competition between introduced black rats and endemic Galápagos rice rats. Ecology 88, 2330–2344 (2007).PubMed 
Article 

Google Scholar 
65.Rayner, M. J., Hauber, M. E., Imber, M. J., Stamp, R. K. & Clout, M. N. Spatial heterogeneity of mesopredator release within an oceanic island system. Proc. Natl. Acad. Sci. 104, 20862–20865 (2007).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
66.MacKay, J. W., Russell, J. C. & Murphy, E. C. Eradicating house mice from islands: successes, failures and the way forward. In Managing Vertebrate Invasive Species: Proceedings of an International Symposium (eds Witmer, G. W. et al.) 294–304 (USDA/APHIS/WS, National Wildlife Research Center, Fort Collins, CO., 2007).67.Capizzi, D. A review of mammal eradications on Mediterranean Islands. Mammal Rev. 50, 124–135 (2020).Article 

Google Scholar 
68.Ballari, S. A., Kuebbing, S. E. & Nuñez, M. A. Potential problems of removing one invasive species at a time: a meta-analysis of the interactions between invasive vertebrates and unexpected effects of removal programs. PeerJ 4, e2029 (2016).69.Anon. No UK products authorised for control of wood mice, https://www.pestcontrolnews.com/rodenticide-alert-prompts-reminder-of-permitted-target-species/ (2020).70.Furness, R. Predation on ground-nesting seabirds by island populations of red deer Cervus elaphus and sheep Ovis. J. Zool. 216, 565–573 (1988).Article 

Google Scholar 
71.Duff, A. Population ecology of the Manx shearwater Puffinus puffinus on the island of Rum, Scotland, MSc thesis, Edinburgh Napier University (2011).72.Towns, D.R., West, C.J. & Broome, K.G. Purposes, outcomes and challenges of eradicating invasive mammals from New Zealand islands: an historical perspective. Wildl. Res. 40, 94–107 (2012). Article 

Google Scholar  More

Temporal variation and dynamic analysis of eucalyptus forestsData from the 1st-9th NFIs suggested that the total area of eucalyptus plantations had started to increase since 1973 (Table 1). In 1973–1976, Eucalyptus plantations only existed in Guangxi and Guangdong (including Hainan) in China with a total area of 23.0 × 104 hectares, taking up 0.38% of total forest area in China. The stock volume was 372.0 × 104 m3, about 0.04% of that in China. In 2014–2018, the eucalyptus plantation area increased to 546.74 × 104 hectares, about 24 times of that in 1973–1976. The growing stock has increased to 21,562.90 × 104 m3 in 2014–2018 which increased about 58 times from 1973–1976.Table 1 Eucalyptus plantation area and stand volume in different time periods by province.Full size tableThe stock volume per unit area of eucalyptus plantations did not increase significantly from 1973–2008, ranging from 14–30 m3/hectares, but it increased rapidly from 2009 to 2018, reaching 39.43 m3/hectares. This increase occurred because China started to focus and value the development of eucalyptus plantations. As a. result, plantations expanded rapidly, and the need for eucalyptus with greater trunk radius increased. The extended harvest cycle of eucalyptus plantations, not the rise of eucalyptus productivity, caused the increase in stock volume per unit area29,30,31.Based on CFLDM, the distribution of eucalyptus plantations in 2003 and 2016 are mapped (Fig. 1a,b). It suggests that the distribution of eucalyptus plantations extended from Leizhou Peninsula, Guangdong and Hainan Province to the north (Guangxi, Hunan, and Guizhou provinces), east (Fujian and Jiangxi provinces), and west (Yunan and Sichuan provinces). This is consistent with data from the NFIs. The widespread expansion of eucalyptus also leads to several regions with clustered plantations.Figure 1Distribution of eucalyptus in the south of China [(a) 2003; (b) 2016]. This figure was created by spatially overlaying spatial sample plots data from National Forest Inventory (NFI) and patch vectors data from China Forest-Land Database Map (CFLDM), (a) shows that the point data are from 6th NFIs eucalyptus sample plots and the polygon vector data are from the 2003 CFLDM; (b) shows that the point data are from 9th NFIs eucalyptus sample plots and the polygon vector data are from the 2016 CFLDM. The extents of eucalyptus plantations is mainly concerned with 11 provinces (e.g. Zhejiang, Fujian, Guangdong, Guangxi) in Southern China.Full size imageChanges in spatial distributionBased on the database of sample plots (including climate and elevation data) and sampled eucalyptus, we analyze the distribution of eucalyptus plantations, and how it is affected by elevation and climate conditions. It is found that most eucalyptus plantations are within the region of 110°23′–120°5′E and 18°21′–30°39′N. The annual mean temperature within this region ranges from 11 to 25 °C with an average of 19.5 °C, and the annual precipitation ranges from 600 to 2000 mm with an average of 1455 mm. Elevation in this region is 0–2500 m with an average of 338 m.To find out the most suitable conditions for eucalyptus growth and its plantation management, we classify this region based on their elevation and climate conditions. The classification is done separately and independently for each factor (i.e., elevation, temperature, and precipitation). In terms of elevation, the region is assigned to seven grades from below 300 m to 2100 m with an interval of 300 m in between (i.e., below 300 m, 300–600 m, 600–900 m, 900–1200 m, 1200–1500 m, 1500–1800 m, 1800–2100 m). Land with elevation above 2100 mm has limited eucalyptus plantations, and thus is not taken into consideration. Similar criterion is applied to the classification based on annual mean temperature and annual precipitation. The grades are from 11 to 25 °C within an interval of 2 °C for temperature and 600–2000 mm with an interval of 200 mm for precipitation.We examine how the eucalyptus forest area changes with these factors. This is done by plotting the eucalyptus plantation area within a certain group against the corresponding grade number (Fig. 2). Eucalyptus is mostly distributed below 300 m, reaching an area of 301.1 × 104 hectares and counting for 67.58% of the total eucalyptus plantation area in China. Eucalyptus occurs rarely in areas with elevation above 900 m.Figure 2Area of eucalyptus plantations in China based on grades defined in text.Full size imageEucalyptus is sensitive to temperature, and its distribution is limited within areas with annual mean temperature below 19 °C. Eucalyptus is mostly distributed within areas with annual mean temperature of 19–21 °C. Areas with annual mean temperature within this range have a total of 291.49 × 104 hectares eucalyptus plantations, approximately 65.43% of all eucalyptus plantations in China. This result is slightly different from previous studies3,5, which suggests that eucalyptus prefers areas with mean annual temperature above 20 °C.Eucalyptus has a high tolerance to annual precipitation. Eucalyptus plantations can be found in areas with annual precipitation ranging from 600–2000 mm. It should be noted that this is related to irrigation conditions in production and management. However, in areas with annual precipitation below 600 mm, Eucalyptus plantations are extremely rare. Areas with annual precipitation of 1400–1600 mm (and without considering other factors) have the largest portion of eucalyptus plantations, whose total area reaches 146.49 × 104 hectares. This accounts for about 32.94% of total eucalyptus plantations in China.Productivity analysis of eucalyptusVariability in mean productivityEucalyptus annual productivity for each province based on the 5th to 9th NFIs (no digitized data for the 1th to 4th NFIs) is calculated, which includes 3564 sample plots (Table 2). Among which, 769 sample plots had been harvested at the time of the survey, and to remove the influence of these 769 sample plots, their data were removed during the productivity of the eucalyptus age-productivity relationship graph (see Fig. 3), which shows that the period of maximum productivity for eucalyptus lasts for approximately 2–3 years. Its productivity declines rapidly after 10 years of growth. Therefore, the harvest cycle of eucalyptus is normally 4–5 years. After coppicing and growing for another 4–5 years, Eucalyptus will be harvested again, which will be followed by its replanting.Table 2 Basic sample plot statistics (quantity, mean and maximum annual productivity) of eucalyptus plantations by province from the 5th to 9th NFIs.Full size tableFigure 3Relationship between age and productivity of eucalyptus sample plots.Full size imageVariability in eucalyptus productivityEucalyptus productivity for each province based on the 5th to 9th NFIs is calculated and shown in Table 2. It can be seen that from 1994 to 2018, mean and maximum productivities of eucalyptus plantations have increased. This is especially show for Guangxi and Fujian Provinces during 2009–2018. The averaged productivity of eucalyptus plantations in China increased from 4.14 to 8.57 m3 hm−2 a−1 from 1994–1998 to 2014–2018, which can be explained by the improved management of eucalyptus plantations (e.g., high soil fertility for newly cultivated lands and improved ability for irrigation) and their expansion.Data from the 5th to 9th NFIs suggest that a lot of sampled plots were no longer used for growing eucalyptus before the next inventory (Table 3). There were 226, 273, 687, and 109 eucalyptus plots in the 5th, 6th, 7th, 8th inventories, respectively, and in the corresponding next NFI (6th, 7th, 8th, 9th), only 150, 179, 544, and 848 of these sample plots were left unabandoned. This suggests that 33.63%, 34.43%, 20.82%, and 22.63% of the plots were abandoned before the next inventory. New plots have been included in each inventory, but large portions of these plots were abandoned as well. There are 123, 508, and 552 new eucclyptus plots in the 6th, 7th, 8th inventories, and 76, 413, and 433 of them were left unabandoned in the next inventories. The land abandonment rates for them are 38.21%, 18.70%, 21.56%, respectively.Table 3 Quantity of newly-cultivated, retained, and abandoned sample plots during different NFIs.Full size tableWe examine how the productivity changes with time for eucalyptus plantations that have been operated for more than 20 years. From the 5th to 9th NFIs, we find 55 and 38 such (operating for more than 25 and 20 years, respectively; Table 4). It is found that the productivity of eucalyptus is relatively low in the first 5 years of its growing. The productivity increases in the 5th–10th years, and reaches its peak after 10–15 years of growing eucalyptus. For land that have been continuously growing eucalytptus for 15–25 years or more, the productivity decreases significantly. This is due to the decrease in soil fertility.Table 4 Relationship between Continuous planting time and Mean productivity of reserved and newly-cultivated eucalyptus sample plots during different NFIs (unit: m3 hm−2 a−1).Full size tableMost (~ 90%) of the sample plots have mean annual productivity below 10 m3 hm−2 a−1 (Fig. 4). The productivity reaches its peak after 10–15 years of growing eucalyptus (mean annual average: 7.0175 m3 hm−2 a−1), and starts to decrease afterwards. Statistical model (Table 5) is established between productivity and age of eucalyptus plots. The results suggest that eucalyptus productivity follows a consistent pattern: it increases with time until a peak and then decrease (Fig. 5), and this applies to the old, newly-included plots and their average. The statistical model also agrees to the observed data, suggesting that the productivity peak is reached 10–15 years after the planting of eucalyptus, and the productivity reaches its minimum or even zero after 50 years of growing eucalyptus continuously.Figure 4Distribution of mean annual productivity for sample plot from different NFIs.Full size imageTable 5 Productivity prediction model of multi-stage reserved and increase eucalyptus sample plots.Full size tableFigure 5Statistical model showing how mean annual productivity of eucalyptus sample plots changes with time.Full size imageSoil fertility variation of eucalyptus plantationsHow eucalyptus affects soil fertility is not well-studied. Here, based on 948 sample points from Tang32, which includes monitoring of soil fertility of eucalyptus plantations from 1993 to 2018, we report and study the temporal soil fertility variation for eucalyptus plantations. After 25 years of growing eucalyptus, acidification of the corresponding lands persists. The pH value changed to 4.63 in 2018, a 4.14% decrease compared to that in1993. The organic content within the soil reached its minimum of 17.98 g/kg in 2018, a decrease of 23.19% compared to 1993. Total nitrogen content of the investigated samples changed from 2.11 to 1.98 g/kg, and total phosphorus content decreased from 1.12 to 0.75 g/kg. The temporal variation of potassium does not change in a consistent pattern with time. Alkaline hydrolysis of nitrogen and available potassium content in 2018 are significantly lower than those in 1993. From more to less, the rank of soil fertility indicator affiliation polygon area is 1993  > 1998  > 2003  > 2013  > 2018  > 2008. The rank of soil fertility index is 1993  > 1998  > 2018  > 2003  > 2013  > 2008. It decreased first, and then increased. The minimum soil fertility (0.475) was reached in 2008 (22.51), which is smaller than that in 1993. The soil fertility decreases at the greatest rate after 15 years of growing eucalyptus. This argument from Tang32 is consistent with this work (Table 6). Soil fertility generally decreases with the age of eucalyptus plantations.Table 6 Evolutionary characteristics of soil chemical indicators in eucalyptus plantation forests.Full size tableIn addition, Parfitt et al.33 studied the variation of soil fertility of pine plantations in New Zealand for a period of 20 years, and found that long-term successive rotations lead to an increase of the soil C/N ratio. Carbon is lost at a speed much greater than nitrogen. Successive rotations of eucalyptus lead to environmental issues such as decrease in soil fertility and ecological diversity and soil erosion. These would limit the sustainable management of eucalyptus plantations34,35,36,37,38,39,40.Abandonment of sample plotsWe find that many sample plots were not used for growing eucalyptus anymore after each inventory. The abandonment rate is high, ranging from 18.7 to 38.21%. The 226 eucalyptus sample plots in the 5th inventory decreased to 103 (the others are abandoned) during the 7th inventory, and the land abandonment rate was 31.33%. In the 8th and 9th inventories, the abandonment rates are 30.10% and 23.61%, respectively. The cumulative land abandonment rates are 33.63%, 54.43%, 68.15%, and 75.66% after 5, 10, 15, and 20 years of growing eucalyptus, respectively (Table 7).Table 7 Quantity (rates) of retained and abandoned sample plots after certain periods of plantation management.Full size tableThere are a total of 1843 eucalyptus plots from the 5th to 9th NFIs. In the last NFI, there are 1282 sampled plots still growing eucalyptus, and the rest 561 plots are abandoned. The averaged land abandonment rates of these plots every 5 years are 23.92%, 24.26% (43.52% cumulatively), 32.10% (68.48% cumulatively), and 23.61% (75.66% cumulatively) over 5, 10, 15, and 20 years, respectively.These data suggest that the abandonment rate of eucalyptus plantations reaches its peak (about one third) after 15 years of operation. For other time intervals (i.e., 5, 10, and 20 years), the rate remains at around 25%. This is related to the management of eucalyptus plantations in the south of China: the first eucalyptus harvest cycle is about 6 years. The second generation of eucalyptus reproduces by division propagation (sprout naturally) with 4 years of harvest cycle, and the third generation follows the same pattern. These amount to 15–16 years-long period for plantation management. Eucalyptus requires stubble-cleaning after twice of division propagations (sprout reproduction), and needs to be re-planted. This is consistent with the timing of abandonment rate peak as stated above. It is highly likely that the eucalyptus plantations are abandoned due to the low soil fertility, and plantation managers or land owners decide to stop growing eucalyptus as a result.A simple statistical model (second-order polynomial) is established between eucalyptus plantations abandonment rate and time (Fig. 6), which suggests that all plantations will stop growing eucalyptus after 50 years, and the corresponding lands will be used for other purposes. The expansion of eucalyptus plantations relies on sustained cultivation of new lands (land reclamation). The total area of eucalyptus plantations reached 5,647,400 hectares in the 9th NFI, but only 4.29% of them (that) have been continuously growing eucalyptus since the 5th NFL (i.e., 24.34% of the plots from the 5th NFI are kept).Figure 6Statistical model showing how abandonment and replanting rates of eucalyptus plantations change with time.Full size imageThere are two main reasons to explain the loss of eucalyptus plantations. The land might be taken over for non-agricultural use (e.g., infrastructure and building construction), or they could be used for growing other crops. The latter is help for soil fertility restoration and soil microorganism readjustment. As most eucalyptus plantations in China are cultivated on lands with poor growing conditions, most of them were abandoned voluntarily by the land owner or plantation manager as stated earlier.After harvest, eucalyptus plantations could be reused for the continuation of eucalyptus growing or used for other purposes (e.g., growing other crops). The plots that were temporarily not used for growing eucalyptus could be used for re-growing it under certain conditions after a certain time period. We investigate the replanting rate of the 561 abandoned eucalyptus plantations, and study whether the abandoned plantations are used for growing other outcrops, and, if so, the corresponding tree species (Tables 8, 9, 10, 11). The 6th NFI data suggests that there are 76 plots abandoned after the 5th NFI. Their replanting rates are 2.63%, 7.89% (10.53% cumulatively), 0.00% (10.53% cumulatively), and 5.26% (15.79% cumulatively) within every 5 years, and after 5, 10, and 15 years of abandonment. For all the 561 abandoned plots, replanting rates are 9.09%, 5.53% (cumulatively 14.62% within 10 years), 0.53% (cumulatively 15.15% within 15 years), 5.26% (cumulatively 15.86% within 20 years) within every 5 years, and after 5, 10, and 15 years of abandonment. These suggest that about five sixth of the abandoned plots had not replanted eucalyptus for at least 20 years since abandonment.Table 8 Land use of eucalyptus plantation sample plots during different NFIs.Full size tableTable 9 Temporal change of tree species planted in sample plots.Full size tableTable 10 Replanting rate of eucalyptus plantations.Full size tableTable 11 Productivity of plantations that have replanted eucalyptus.Full size tableA simple statistical model is established between replanting rate and time (second-order polynomial) based on the current data (Fig. 5). It suggests that the replanting rates after 30 and 50 years are around 20% and 30%, respectively. These suggest that if the plantation management does not improve significantly, it would be difficult to maintain the current supply of eucalyptus and areal distribution of its plantations in the long term. It is necessary to rely on both land rotation and cultivation of new lands to maintain the current supply of eucalyptus.The NFI data suggests that very few eucalyputs plots are turned to non-plantation purposes. The exception is from the 6th NFI in which 34.21% of plots have been used for other purposes after harvest. This rate is below 20% for all other inventory data. A lot of abandoned eucalyptus plantations are still used as plantations, and they are for growing eucalyptus, and the rate of regrowing eucalyptus tends to remain low for a long period of time (below one sixth after 20 years based on current data). This is because eucalyptus grows fast with high productivity, and it has high demand for soil fertility and water. Land rotation is necessary after a few harvest cycles to restore the soil fertility, which would take relatively long period of time before the land becomes suitable to regrow euccalyptus. Among the 561 abandoned eucalyptus plots, broad-leaf and economic tree species are the most commonly planted species after stop growing eucalyptus (e.g., rubber tree and Lychee; 18.54% and 18.36%; Table 9). The greater variability of land use for the abandoned plots suggests greater management intensity. Afforestation with eucalyptus is dominated by short rotation period (harvest cycle). Frequently modifying tree species planted within plantations helps maintain a high productivity of the land.
Carbon storage and fixation of eucalyptus plantationsVariability in carbon storageBased on the 9th NFI data and Eqs. (2–8), we calculate the BEF of eucalyptus in each province (Table 12). The results suggest that the BEF ranges from 0.982–1.652 with a weighted average of 1.236 (weight determined by stock volume).Table 12 BEF of eucalyptus by province.Full size tableCalculation from Eqs. (9) and (10) suggests that the total carbon storage (excluding harvest volume) of eucalyptus in China is 2.40 TgC (1973–1976, 1Tg = 1012 g), 4.14 TgC (1977–1981), 2.73 TgC (1984–1988), 5.42 TgC (1989–1993), 9.73 TgC (1994–1998), 12.58 TgC (1999–2003), 28.90 TgC (2004–2008), 98.61 TgC (2009–2013), and 133.00 TgC (2014–2018) in different time periods in the past 45 years (Table 13).Table 13 Eucalyptus carbon density and storage by province.Full size tableThe carbon storage of eucalyptus increased rapidly in the past 45 years especially since the end of last century. This is due to the rapid expansion of eucalyptus plantations in China, and its carbon storage in 2014–2018 is 55.42 times of that in 1973–1976. The carbon density per square hectometer also increases from 5.22 MgC (1 Mg = 106 g) in 1973–1976 to 12.16 MgC in 2014–2018, about 2.33 times of the former.Carbon fixationThe mean annual productivity of eucalyptus is 8.57 m3 hm−2 a−1 in 2014–2018 based on the 9th NFI. This is a lot greater compared to other species widely planted in the same areas (Pinus massoniana Lamb.: 2.91 m3 hm−2 a−1; Cunninghamia lanceolata Lamb.: 3.93 m3 hm−2 a−1). Using the stock volume biomass method with BEF being 1.2336 (from previous calculation) and carbon storage coefficient of 0.5, the mean annual carbon fixed by eucalyptus is 5.29 t hm−2 a−1, which are about 2.95 and 2.18 times that of Pinus massoniana Lamb. (1.79 t hm−2 a−1) and Cunninghamia lanceolata Lamb. (2.42 t hm−2 a−1), respectively. This shows that eucalyptus is characterized by high biomass productivity and high carbon fixation capability. It thus plays an important role in maintaining the carbon balance in China. More

1.Laurance, W. F., Sayer, J. & Cassman, K. G. Agricultural expansion and its impacts on tropical nature. Trends Ecol. Evol. 29, 107–116 (2014).Article 

Google Scholar 
2.Bélanger, J. & Pilling, D. (eds) The State of the World’s Biodiversity for Food and Agriculture (FAO Commission on Genetic Resources for Food and Agriculture, 2019).3.Rockström, J. et al. Sustainable intensification of agriculture for human prosperity and global sustainability. Ambio 46, 4–17 (2017).Article 

Google Scholar 
4.Pretty, J. et al. Global assessment of agricultural system redesign for sustainable intensification. Nat. Sustain. 1, 441–446 (2018).Article 

Google Scholar 
5.Renard, D. & Tilman, D. National food production stabilized by crop diversity. Nature 571, 257–260 (2019).CAS 
Article 

Google Scholar 
6.Tilman, D. Benefits of intensive agricultural intercropping. Nat. Plants 6, 604–605 (2020).Article 

Google Scholar 
7.Gomez, A. A. & Gomez, K. A. Multiple Cropping in the Humid Tropics of Asia (IDRC, 1983).8.Li, L. et al. Diversity enhances agricultural productivity via rhizosphere phosphorus facilitation on phosphorus-deficient soils. Proc. Natl Acad. Sci. USA 104, 11192–11196 (2007).CAS 
Article 

Google Scholar 
9.Li, C. J. et al. Syndromes of production in intercropping impact yield gains. Nat. Plants 6, 653–660 (2020).Article 

Google Scholar 
10.Xu, Z. et al. Intercropping maize and soybean increases efficiency of land and fertilizer nitrogen use; a meta-analysis. Field Crops Res. 246, 107661 (2020).Article 

Google Scholar 
11.Zhu, Y. Y. et al. Genetic diversity and disease control in rice. Nature 406, 718–722 (2000).CAS 
Article 

Google Scholar 
12.Li, W. X. et al. Effects of intercropping and nitrogen application on nitrate present in the profile of an Orthic Anthrosol in Northwest China. Agric. Ecosyst. Environ. 105, 483–491 (2005).CAS 
Article 

Google Scholar 
13.Manevski, K., Borgesen, C. D., Andersen, M. N. & Kristensen, I. S. Reduced nitrogen leaching by intercropping maize with red fescue on sandy soils in North Europe: a combined field and modeling study. Plant Soil 388, 67–85 (2015).CAS 
Article 

Google Scholar 
14.Tilman, D., Reich, P. B. & Knops, J. M. H. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441, 629–632 (2006).CAS 
Article 

Google Scholar 
15.Roscher, C. et al. Identifying population- and community-level mechanisms of diversity–stability relationships in experimental grasslands. J. Ecol. 99, 1460–1469 (2011).Article 

Google Scholar 
16.Zhou, B. R. et al. Plant functional groups asynchrony keep the community biomass stability along with the climate change—a 20-year experimental observation of alpine meadow in eastern Qinghai-Tibet Plateau. Agric. Ecosyst. Environ. 282, 49–57 (2019).Article 

Google Scholar 
17.Schnabel, F. et al. Drivers of productivity and its temporal stability in a tropical tree diversity experiment. Glob. Change Biol. 25, 4257–4272 (2019).Article 

Google Scholar 
18.Pohl, M., Alig, D., Körner, C. & Rixen, C. Higher plant diversity enhances soil stability in disturbed alpine ecosystems. Plant Soil 324, 91–102 (2009).CAS 
Article 

Google Scholar 
19.Pérès, G. et al. Mechanisms linking plant community properties to soil aggregate stability in an experimental grassland plant diversity gradient. Plant Soil 373, 285–299 (2013).Article 

Google Scholar 
20.Gould, I. J., Quinton, J. N., Weigelt, A., De Deyn, G. B. & Bardgett, R. D. Plant diversity and root traits benefit physical properties key to soil function in grasslands. Ecol. Lett. 19, 1140–1149 (2016).Article 

Google Scholar 
21.Lange, M. et al. Plant diversity increases soil microbial activity and soil carbon storage. Nat. Commun. 6, 6707 (2015).CAS 
Article 

Google Scholar 
22.Dybzinski, R., Fargione, J. E., Zak, D. R., Fornara, D. & Tilman, D. Soil fertility increases with plant species diversity in a long-term biodiversity experiment. Oecologia 158, 85–93 (2008).Article 

Google Scholar 
23.Tisdall, J. M. & Oades, J. M. Organic matter and water-stable aggregates in soils. J. Soil Sci. 33, 141–163 (1982).CAS 
Article 

Google Scholar 
24.Amézketa, E. Soil aggregate stability: a review. J. Sustain. Agric. 14, 83–151 (1999).Article 

Google Scholar 
25.Six, J., Bossuyt, H., Degryze, S. & Denef, K. A history of research on the link between (micro)aggregates, soil biota, and soil organic matter dynamics. Soil Tillage Res. 79, 7–31 (2004).Article 

Google Scholar 
26.Tiemann, L. K., Grandy, A. S., Atkinson, E. E., Marin-Spiotta, E. & McDaniel, M. D. Crop rotational diversity enhances belowground communities and functions in an agroecosystem. Ecol. Lett. 18, 761–771 (2015).CAS 
Article 

Google Scholar 
27.Christensen, B. T. & Johnston, A. E. in Developments in Soil Science (eds Gregorich, E. G. & Carter, M. R.) 399–430 (Elsevier, 1997).28.Fornara, D. A. & Tilman, D. Soil carbon sequestration in prairie grasslands increased by chronic nitrogen addition. Ecology 93, 2030–2036 (2012).Article 

Google Scholar 
29.Cong, W. F. et al. Plant species richness promotes soil carbon and nitrogen stocks in grasslands without legumes. J. Ecol. 102, 1163–1170 (2014).CAS 
Article 

Google Scholar 
30.Zhang, W. F. et al. Closing yield gaps in China by empowering smallholder farmers. Nature 537, 671–674 (2016).CAS 
Article 

Google Scholar 
31.Wan, N. F. et al. Global synthesis of effects of plant species diversity on trophic groups and interactions. Nat. Plants 6, 503–510 (2020).Article 

Google Scholar 
32.Isbell, F. et al. Benefits of increasing plant diversity in sustainable agroecosystems. J. Ecol. 105, 871–879 (2017).Article 

Google Scholar 
33.Status of the World’s Soil Resources Main Report (FAO and ITPS, 2015).34.Jensen, E. S. & Williamson, S. in Replacing Chemicals with Biology: Phasing Out Highly Hazardous Pesticides with Agroecology (eds Watts, M. & Williamson, S.) 162–167 (Pestice Action Network International, 2015).35.Bender, S. F., Wagg, C. & van der Heijden, M. G. A. An underground revolution: biodiversity and soil ecological engineering for agricultural sustainability. Trends Ecol. Evol. 31, 440–452 (2016).Article 

Google Scholar 
36.Zhang, R. Z. et al. Response of the arbuscular mycorrhizal fungi diversity and community in maize and soybean rhizosphere soil and roots to intercropping systems with different nitrogen application rates. Sci. Total Environ. 740, 139810 (2020).CAS 
Article 

Google Scholar 
37.Bao, S. D. Analysis on Soil and Agricultural Chemistry (China Agricultural Press, 2005).38.Kemper, W. D. & Rosenau, R. C. Aggregate Stability and Size Distribution: Methods of Soil Analysis, Part 1 Physical and Mineralogical Methods (American Society of Agronomy, 1986).39.Willey, R. W. Evaluation and presentation of intercropping advantages. Exp. Agric. 21, 119–133 (1985).Article 

Google Scholar 
40.van Ruijven, J. & Berendse, F. Contrasting effects of diversity on the temporal stability of plant populations. Oikos 116, 1323–1330 (2007).Article 

Google Scholar 
41.The Price Division, National Development and Reform Commission of China The National Agricultural Products Cost–Benefit Compilation of Information (China Statistics Press, 2010–2019).42.R Core Team R: A Language and Environment for Statistical Computing v.3.6.2 (R Foundation for Statistical Computing, 2019). More