Ecology

Ritchie, H. & Roser, M. Urbanization. Our World in Data (2018). https://ourworldindata.org/urbanization.Chen, H., Weersink, A., Beaulieu, M., Lee, Y. N. & Nagelschmitz, K. A historical review of changes in farm size in canada. Tech. Rep., University of Guelph, Institute for the Advanced Study of Food and and Agricultural Policy (2019).Gudelj, I. & White, K. Spatial heterogeneity, social structure and disease dynamics of animal populations. Theor. Popul. Biol. 66, 139–149 (2004).CAS 
MATH 
Article 

Google Scholar 
Augustin, N., Mugglestone, M. A. & Buckland, S. T. An autologistic model for the spatial distribution of wildlife. J. Appl. Ecol. 339–347 (1996).Karlsson, E. K., Kwiatkowski, D. P. & Sabeti, P. C. Natural selection and infectious disease in human populations. Nat. Rev. Genetics 15, 379–393 (2014).CAS 
Article 

Google Scholar 
Fornaciari, A. Environmental microbial forensics and archaeology of past pandemics. Microbiol. Spect. 5, 5–1 (2017).Article 

Google Scholar 
Thèves, C., Crubézy, E. & Biagini, P. History of smallpox and its spread in human populations. Microbiol. Spect. 4, 4–4 (2016).Article 

Google Scholar 
Coltart, C. E., Lindsey, B., Ghinai, I., Johnson, A. M. & Heymann, D. L. The ebola outbreak, 2013–2016: old lessons for new epidemics. Philos. Trans. R. Soc. B Biol. Sci. 372, 20160297 (2017).Article 

Google Scholar 
Colizza, V., Barrat, A., Barthelemy, M., Valleron, A.-J. & Vespignani, A. Modeling the worldwide spread of pandemic influenza: Baseline case and containment interventions. PLOS Med. 4, e13 (2007).Article 

Google Scholar 
Lüthy, I. A., Ritacco, V. & Kantor, I. N. One hundred years after the “Spanish” flu. Medicina 78, 113–118 (2018).
Google Scholar 
Zhang, Y., Zhang, A. & Wang, J. Exploring the roles of high-speed train, air and coach services in the spread of COVID-19 in China. Transport Policy 94, 34–42 (2020).Article 

Google Scholar 
Coelho, M. T. P. et al. Global expansion of COVID-19 pandemic is driven by population size and airport connections. PeerJ 8, e9708 (2020).Article 

Google Scholar 
Tompkins, D. M., Carver, S., Jones, M. E., Krkošek, M. & Skerratt, L. F. Emerging infectious diseases of wildlife: A critical perspective. Trends Parasitol. 31, 149–159 (2015).Article 

Google Scholar 
Soulsbury, C. D. & White, P. C. Human-wildlife interactions in urban areas: A review of conflicts, benefits and opportunities. Wildl. Res. 42, 541–553 (2015).Article 

Google Scholar 
VanderWaal, K. L. et al. Network analysis of cattle movements in Uruguay: Quantifying heterogeneity for risk-based disease surveillance and control. Prevent. Vet. Med. 123, 12–22 (2016).Article 

Google Scholar 
Rossi, G. et al. The potential role of direct and indirect contacts on infection spread in dairy farm networks. PLOS Comput. Biol. 13, e1005301 (2017).Article 

Google Scholar 
Stoddard, S. T. et al. The role of human movement in the transmission of vector-borne pathogens. PLOS Neg. Trop. Dis. 3, e481 (2009).Article 

Google Scholar 
Cosner, C. Models for the effects of host movement in vector-borne disease systems. Math. Biosci. 270, 192–197 (2015).MathSciNet 
MATH 
Article 

Google Scholar 
Scherer, P.C. Infection on the move: individual host movement drives disease persistence in spatially structured landscapes. Dr. rer. nat. thesis, Universität Potsdam (2019).Riley, S. Large-scale spatial-transmission models of infectious disease. Science 316, 1298–1301 (2007).ADS 
CAS 
Article 

Google Scholar 
Dougherty, E. R., Seidel, D. P., Carlson, C. J., Spiegel, O. & Getz, W. M. Going through the motions: Incorporating movement analyses into disease research. Ecol. Lett. 21, 588–604 (2018).Article 

Google Scholar 
Daversa, D., Fenton, A., Dell, A., Garner, T. & Manica, A. Infections on the move: How transient phases of host movement influence disease spread. Proc. R. Soc. B Biol. Sci. 284, 20171807 (2017).Article 

Google Scholar 
MacArthur, R. H. & Wilson, E. O. The Theory of Island Biogeography (Princeton University Press, 2001).Book 

Google Scholar 
Kobayashi, K. & Okumura, M. The growth of city systems with high-speed railway systems. Ann. Region. Sci. 31, 39–56 (1997).Article 

Google Scholar 
VanderWaal, K., Perez, A., Torremorrell, M., Morrison, R. M. & Craft, M. Role of animal movement and indirect contact among farms in transmission of porcine epidemic diarrhea virus. Epidemics 24, 67–75 (2018).Article 

Google Scholar 
Hanski, I. Metapopulation dynamics. Nature 396, 41–49 (1998).ADS 
CAS 
Article 

Google Scholar 
Colizza, V. & Vespignani, A. Epidemic modeling in metapopulation systems with heterogeneous coupling pattern: Theory and simulations. J. Theor. Biol. 251, 450–467 (2008).ADS 
MathSciNet 
MATH 
Article 

Google Scholar 
Wang, L. & Li, X. Spatial epidemiology of networked metapopulation: An overview. Chin. Sci. Bull. 59, 3511–3522 (2014).Article 

Google Scholar 
Ruxton, G. D. Low levels of immigration between chaotic populations can reduce system extinctions by inducing asynchronous regular cycles. Proc. R. Soc. Lond. Seri. B Biol. Sci. 256, 189–193 (1994).ADS 
Article 

Google Scholar 
Earn, D. J. D., Rohani, P. & Grenfell, B. T. Persistence chaos and synchrony in ecology and epidemiology. Proc. R. Soc. Lond. Seri. B Biol. Sci. 265, 7–10 (1998).CAS 
Article 

Google Scholar 
Rosenzweig, M. L. Paradox of enrichment: Destabilization of exploitation ecosystems in ecological time. Science 171, 385–387 (1971).ADS 
CAS 
Article 

Google Scholar 
Hilker, F. M. & Schmitz, K. Disease-induced stabilization of predator-prey oscillations. J. Theor. Biol. 255, 299–306 (2008).ADS 
MathSciNet 
MATH 
Article 

Google Scholar 
Brown, J. H. & Kodric-Brown, A. Turnover rates in insular biogeography: Effect of immigration on extinction. Ecology 58, 445–449 (1977).Article 

Google Scholar 
Philipson, T. Economic epidemiology and infectious diseases. Handb. Health Econ. 1, 1761–1799 (2000).Article 

Google Scholar 
Murdoch, W. W., Briggs, C. J. & Nisbet, R. M. Consumer-Resource Dynamics, Monographs in Population Biology Vol. 36 (Princeton University Press, 2003).
Google Scholar 
Murdoch, W. W. & Oaten, A. Predation and population stability. In Advances in Ecological Research, vol. 9, 1–131 (Elsevier, 1975).Bolker, B. & Grenfell, B. T. Space, persistence and dynamics of measles epidemics. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 348, 309–320 (1995).ADS 
CAS 
Article 

Google Scholar 
Keeling, M. J. & Rohani, P. Estimating spatial coupling in epidemiological systems: a mechanistic approach. Ecol. Lett. 5, 20–29 (2002).Article 

Google Scholar 
Arino, J. Spatio-temporal spread of infectious pathogens of humans. Infect. Dis. Model. 2, 218–228 (2017).
Google Scholar 
Keeling, M. J. & Rohani, P. Modeling Infectious Diseases in Humans and Animals (Princeton University Press, 2011).MATH 
Book 

Google Scholar 
Wilson, E. B. & Worcester, J. The spread of an epidemic. Proc. Nat. Acad. Sci. 31, 327 (1945).ADS 
CAS 
Article 

Google Scholar 
Rushton, S. & Mautner, A. The deterministic model of a simple epidemic for more than one community. Biometrika 42, 126–132 (1955).MathSciNet 
MATH 
Article 

Google Scholar 
Lourenço, J. & Recker, M. Natural, persistent oscillations in a spatial multi-strain disease system with application to dengue. PLOS Comput. Biol. 9, e1003308 (2013).ADS 
Article 

Google Scholar 
Wikramaratna, P. S., Pybus, O. G. & Gupta, S. Contact between bird species of different lifespans can promote the emergence of highly pathogenic avian influenza strains. Proc. Natl. Acad. Sci. 111, 10767–10772 (2014).ADS 
CAS 
Article 

Google Scholar 
Xiao, Y., Zhou, Y. & Tang, S. Modelling disease spread in dispersal networks at two levels. Math. Med. Biol. J. IMA 28, 227–244 (2011).MathSciNet 
MATH 
Article 

Google Scholar 
Arino, J., Ducrot, A. & Zongo, P. A metapopulation model for malaria with transmission-blocking partial immunity in hosts. J. Math. Biol. 64, 423–448 (2012).MathSciNet 
MATH 
Article 

Google Scholar 
De Roos, A. M., Mccauley, E. & Wilson, W. G. Mobility versus density-limited predator-prey dynamics on different spatial scales. Proc. R. Soc. Lond. Ser. B Biol. Sci. 246, 117–122 (1991).ADS 
Article 

Google Scholar 
Dey, S., Goswami, B. & Joshi, A. Effects of symmetric and asymmetric dispersal on the dynamics of heterogeneous metapopulations: Two-patch systems revisited. J. Theor. Biol. 345, 52–60 (2014).ADS 
MathSciNet 
MATH 
Article 

Google Scholar 
Anderson, R. M., Jackson, H. C., May, R. M. & Smith, A. M. Population dynamics of fox rabies in Europe. Nature 289, 765–771 (1981).ADS 
CAS 
Article 

Google Scholar 
Gupta, S., Ferguson, N. & Anderson, R. Chaos persistence, and evolution of strain structure in antigenically diverse infectious agents. Science 280, 912–915 (1998).ADS 
CAS 
Article 

Google Scholar 
Holland, M. D. & Hastings, A. Strong effect of dispersal network structure on ecological dynamics. Nature 456, 792–794 (2008).ADS 
CAS 
Article 

Google Scholar 
McCann, K., Hastings, A. & Huxel, G. R. Weak trophic interactions and the balance of nature. Nature 395, 794–798 (1998).ADS 
CAS 
Article 

Google Scholar 
Singh, A. & Gakkhar, S. Synchronization of chaos in a food web in ecological systems. World Acad. Sci. Eng. Technol. 70, 94–98 (2010).
Google Scholar 
Gotelli, N. J. Metapopulation models: The rescue effect, the propagule rain, and the core-satellite hypothesis. American Naturalist 138, 768–776 (1991).Article 

Google Scholar 
Heino, M., Kaitala, V., Ranta, E. & Lindström, J. Synchronous dynamics and rates of extinction in spatially structured populations. Proc. R. Soc. Lond. Ser. B Biol. Sci. 264, 481–486 (1997).ADS 
Article 

Google Scholar 
Molofsky, J. & Ferdy, J.-B. Extinction dynamics in experimental metapopulations. Proc. Natl. Acad. Sci. 102, 3726–3731 (2005).ADS 
CAS 
Article 

Google Scholar 
Saxena, G., Prasad, A. & Ramaswamy, R. Amplitude death: The emergence of stationarity in coupled nonlinear systems. Phys. Rep. 521, 205–228 (2012).ADS 
Article 

Google Scholar 
Majhi, S. & Ghosh, D. Amplitude death and resurgence of oscillation in networks of mobile oscillators. Europhys. Lett. 118, 40002 (2017).ADS 
Article 

Google Scholar 
Shen, C., Chen, H. & Hou, Z. Mobility and density induced amplitude death in metapopulation networks of coupled oscillators. Chaos 24, 043125 (2014).MATH 
Article 

Google Scholar 
Karnatak, R., Ramaswamy, R. & Feudel, U. Conjugate coupling in ecosystems: Cross-predation stabilizes food webs. Chaos Solitons Fractals 68, 48–57 (2014).ADS 
MathSciNet 
MATH 
Article 

Google Scholar 
Bolker, B. M. & Grenfell, B. T. Chaos and biological complexity in measles dynamics. Proc. R. Soc. Lond. Ser. B Biol. Sci. 251, 75–81 (1993).ADS 
CAS 
Article 

Google Scholar 
Olsen, L. F., Truty, G. L. & Schaffer, W. M. Oscillations and chaos in epidemics: A nonlinear dynamic study of six childhood diseases in Copenhagen, Denmark. Theor. Popul. Biol. 33, 344–370 (1988).MathSciNet 
CAS 
MATH 
Article 

Google Scholar 
Lundberg, P., Ranta, E., Ripa, J. & Kaitala, V. Population variability in space and time. Trends Ecol. Evolut. 15, 460–464 (2000).CAS 
Article 

Google Scholar 
Dekker, A. Realistic social networks for simulation using network rewiring. In International Congress on Modelling and Simulation, 677–683 (2007).Milgram, S. The small world problem. Psychol. Today 2, 60–67 (1967).
Google Scholar 
Sallaberry, A., Zaidi, F. & Melançon, G. Model for generating artificial social networks having community structures with small-world and scale-free properties. Soc. Netw. Anal. Min. 3, 597–609 (2013).Article 

Google Scholar 
Olesen, J. M., Bascompte, J., Dupont, Y. L. & Jordano, P. The modularity of pollination networks. Proc. Natl. Acad. Sci. 104, 19891–19896 (2007).ADS 
CAS 
MATH 
Article 

Google Scholar 
Stouffer, D. B. & Bascompte, J. Compartmentalization increases food-web persistence. Proc. Natl. Acad. Sci. 108, 3648–3652 (2011).ADS 
CAS 
Article 

Google Scholar 
Girvan, M. & Newman, M. E. Community structure in social and biological networks. Proc. Natl. Acad. Sci. 99, 7821–7826 (2002).ADS 
MathSciNet 
CAS 
MATH 
Article 

Google Scholar 
Krause, A. E., Frank, K. A., Mason, D. M., Ulanowicz, R. E. & Taylor, W. W. Compartments revealed in food-web structure. Nature 426, 282–285 (2003).ADS 
CAS 
Article 

Google Scholar 
Rezende, E. L., Albert, E. M., Fortuna, M. A. & Bascompte, J. Compartments in a marine food web associated with phylogeny, body mass, and habitat structure. Ecol. Lett. 12, 779–788 (2009).Article 

Google Scholar 
Pastor-Satorras, R. & Vespignani, A. Epidemics and immunization in scale-free networks. In Handbook of Graphs and Networks, 111–130 (Wiley Online Library, 2002).Lloyd-Smith, J. O., Schreiber, S. J., Kopp, P. E. & Getz, W. M. Superspreading and the effect of individual variation on disease emergence. Nature 438, 355–359 (2005).ADS 
CAS 
Article 

Google Scholar 
Shirley, M. D. & Rushton, S. P. The impacts of network topology on disease spread. Ecol. Complex. 2, 287–299 (2005).Article 

Google Scholar 
Keeling, M. J. & Eames, K. T. Networks and epidemic models. J. R. Soc. Interface 2, 295–307 (2005).Article 

Google Scholar 
Godfrey, S. S., Bull, C. M., James, R. & Murray, K. Network structure and parasite transmission in a group living lizard the gidgee skink, Egernia stokesii. Behav. Ecol. Sociobiol. 63, 1045–1056 (2009).Article 

Google Scholar 
VanderWaal, K. L., Atwill, E. R., Hooper, S., Buckle, K. & McCowan, B. Network structure and prevalence of Cryptosporidium in Belding’s ground squirrels. Behav. Ecol. Sociobiol. 67, 1951–1959 (2013).Article 

Google Scholar 
Proulx, S. R., Promislow, D. E. & Phillips, P. C. Network thinking in ecology and evolution. Trends Ecol. Evolut. 20, 345–353 (2005).Article 

Google Scholar 
Craft, M. E. & Caillaud, D. Network models: An underutilized tool in wildlife epidemiology?. Interdiscip. Perspect. Infect. Dis. 2011, (2011).Bajardi, P. et al. Human mobility networks, travel restrictions, and the global spread of 2009 h1n1 pandemic. PloS one 6, e16591 (2011).ADS 
CAS 
Article 

Google Scholar 
Gog, J. R. et al. Seven challenges in modeling pathogen dynamics within-host and across scales. Epidemics 10, 45–48 (2015).Article 

Google Scholar 
Cen, X., Feng, Z. & Zhao, Y. Emerging disease dynamics in a model coupling within-host and between-host systems. J. Theor. Biol. 361, 141–151 (2014).ADS 
MathSciNet 
MATH 
Article 

Google Scholar 
Meakin, S. R. & Keeling, M. J. Correlations between stochastic epidemics in two interacting populations. Epidemics 26, 58–67 (2019).Article 

Google Scholar 
Machado, G. et al. Identifying outbreaks of porcine epidemic diarrhea virus through animal movements and spatial neighborhoods. Sci. Rep. 9, 1–12 (2019).
Google Scholar 
Tonkin, J. D. et al. The role of dispersal in river network metacommunities: Patterns, processes, and pathways. Freshwater Biol. 63, 141–163 (2018).Article 

Google Scholar 
Pedersen, T. L. tidygraph: a tidy API for graph manipulation (2019). R package version 1.1.2.Rackauckas, C. & Nie, Q. Differentialequations.jl–a performant and feature-rich ecosystem for solving differential equations in Julia. J. Open Res. Softw. 5 (2017).Rackauckas, C. & Nie, Q. Confederated modular differential equation APIS for accelerated algorithm development and benchmarking. Adv. Eng. Softw. 132, 1–6 (2019).Article 

Google Scholar 
Bezanson, J., Edelman, A., Karpinski, S. & Shah, V. B. Julia: A fresh approach to numerical computing. SIAM Rev. 59, 65–98 (2017).MathSciNet 
MATH 
Article 

Google Scholar 
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).MATH 
Book 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2020). More

Zhu, X. China’s karst tiankeng and its value for science and tourism. Sci. Technol. Rev. 19, 60–63 (2001).
Google Scholar 
Zhu, X. et al. A brift study on karst tiankeng. Carsol. Sin. 22, 51–65 (2003).
Google Scholar 
Zhu, X. & Waltham, T. Tiankeng: Definition and description. Cave Karst Sci. 32, 75–79 (2005).
Google Scholar 
Zhu, X. & Chen, W. Tiankengs in the karst of China. Cave Karst Sci. 32, 55–56 (2005).
Google Scholar 
Shui, W., Chen, Y., Wang, Y., Su, Z. & Zhang, S. Origination, study progress and prospect of karst tiankeng research in China. Acta Geogr. Sin. 70, 431–446 (2015).
Google Scholar 
Alexander, K. Cave un-roofing as a large-scale geomorphic process. Carsol. Sin. 4, 1–11 (2006).
Google Scholar 
Palmer, A. & Palmer, M. Hydraulic processes in the origin of tiankengs. Speleogenesis Evol. Karst Aquifers 4, 8 (2006).
Google Scholar 
Waltham, T. Collapse processes at the tiankengs of Xingwen. Cave Karst Sci. 32, 107–110 (2005).
Google Scholar 
White, W. & White, E. Size scales for closed depression landforms: The place of tiankengs. Cave Karst Sci. 32, 111–118 (2005).
Google Scholar 
Chen, W., Zhu, X., Zhu, D. & Ma, Z. Karst geological relics and development of Xiaozhai Tiankeng and Tianjinxia Fissure Gorge, Fengjie, Chongqing. J. Mountain Sci. 22, 22–29 (2004).CAS 

Google Scholar 
Yue, Y., Wang, K., Zhang, W., Chen, H. & Wang, M. Relationships between soil and environment in Peak-Cluster Depression areas of karst region based on canonical correspondence analysis. Environ. Sci. 29, 1400–1405 (2008).
Google Scholar 
Huang, B., Cai, W., Xue, Y. & Zhu, X. Research on tourism resource characteristics of tiankeng group in Dashiwei, Guangxi. Geogr. Geo-Inf. Sci. 20, 109–112 (2004).CAS 

Google Scholar 
Zhu, X. Discovery of erosional tiankeng in Houping, Wulong and its value of science and tourism. Carsol. Sin. 2, 93–98 (2006).
Google Scholar 
Yuan, D. The development of modern karstology in China. Geol. Rev. 52, 733–736 (2006).CAS 

Google Scholar 
Gunn, J. Turloughs and tiankengs: Distinctive doline forms. Cave Karst Sci. 32, 83–84 (2005).
Google Scholar 
Klimchou, A. Cave un-roofing as a large-scale geomorphic process. Cave Karst Sci. 32, 93–98 (2005).
Google Scholar 
Zhu, X., Chen, W. & Erin, L. Wulong karst systems and as an indicator of local tectonic uplift. Carsol. Sin. 26, 119–125 (2007).
Google Scholar 
Shui, W. & Wang, X. Geological expedition and analysis on formation and evolvement of erosive Karst Tiankeng: A case study of Xingwen World Geopark. Adv. Mater. Res. 250–253, 2002–2006 (2011).Article 

Google Scholar 
Su, S., Huang, K. & Ma, B. Diversity study on pteridophyte flora in the area of Dashiwei Tiankeng group of Leye County. Hubei Agric. Sci. 51, 948–950 (2012).
Google Scholar 
Huang, K. & Su, S. Resource investigation and application research of pteridophyte flora resource in the area of Dashiwei Tiankeng Group. Anhui Agric. Sci. Bull. 21, 74–80 (2015).CAS 

Google Scholar 
Su, Y., Xue, Y., Fan, B., Mo, F. & Feng, H. Plant community structure and species diversity in Liuxing tiankeng of Guangxi. Acta Botan. Boreali-Occiden. Sin. 36, 2300–2306 (2016).
Google Scholar 
Li, W., Xiang, Y., Du, Y. & Wu, X. Underground forest communities in Zhanyi, Yunnan Province. For. Sci. Technol. 20–25 (2001).Jian, X. et al. Interspecific relationships of grassland plant community’s dominant species in moderate-degraded tiankeng of Yunnan, China. Chin. J. Appl. Ecol. 29, 1–14. https://doi.org/10.13287/j.1001 (2018).Article 

Google Scholar 
Chen, W., Zhu, D. & Zhu, X. Characteristics and evaluation of karst landscape in tiankeng-difeng scenery site, Fengjie, Chongqing. Geogr. Geo-Inf. Sci. 20, 80–83 (2004).CAS 

Google Scholar 
Wang, J. & Guo, C. Comparison between the positive and negative topographic ecosystem in karst mountainous areas and its bearing capability. Guizhou Agric. Sci. 35, 85–87 (2007).MathSciNet 

Google Scholar 
Tony, W. Tiankengs of the world, outside China. Speleogenesis Evol. Karst Aquifers 4, 1–12 (2006).
Google Scholar 
Su, Y., Tang, Q., Mo, F. & Xue, Y. Karst tiankengs as refugia for indigenous tree flora amidst a degraded landscape in southwestern China. Sci. Rep. 7, 1–10 (2017).ADS 
Article 
CAS 

Google Scholar 
Bátori, Z. et al. A comparison of the vegetation of forested and non-forested solution dolines in Hungary: A preliminary study. Biologia 69, 1339–1348 (2014).Article 
CAS 

Google Scholar 
Bátori, Z. et al. The conservation value of karst dolines for vascular plants in woodland habitats of Hungary: Refugia and climate change. Int. J. Speleol. 43, 15–26 (2014).Article 

Google Scholar 
Bátori, Z. et al. Importance of karst sinkholes in preserving relict, mountain, and wet-woodland plant species under sub-Mediterranean climate: A case study from southern Hungary. J. Cave Karst Stud. Natl. Speleol. Soc. Bull. 74, 127–134 (2012).Article 

Google Scholar 
Bátori, Z. et al. Large-and small-scale environmental factors drive distributions of cool-adapted plants in karstic microrefugia. Ann. Bot. 119, 301–309 (2017).PubMed 
Article 

Google Scholar 
Vilisics, F. et al. Small scale gradient effects on isopods (Crustacea: Oniscidea) in karstic sinkholes. Biologia 66, 499–505 (2011).Article 

Google Scholar 
Dolinar, B. & Vreš, B. Pregled flore Mišje doline in zgornjega porečja Rašice (Dolenjska, Slovenija). Hladnikia 30, 3–37 (2012).
Google Scholar 
Raschmanová, N., Miklisová, D., Ľubomír, K. & Šustr, V. Community composition and cold tolerance of soil Collembola in a collapse karst doline with strong microclimate inversion. Biologia 70, 802–811 (2016).Article 

Google Scholar 
Macarthur, R. & Wilson, E. The Theory of Island Biogeography (Princeton University Press, 1967).
Google Scholar 
Kreft, H. & Jetz, W. Global patterns and determinants of vascular plant diversity. Proc. Natl. Acad. Sci. U. S. A. 104, 5925 (2007).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhou, F. & Yu, S. Study on the plant diversity of island-like habitats in Karst Mountain Areas. Guizhou For. Sci. Technol. 40, 18–22 (2012).CAS 

Google Scholar 
Hu, F., Lou, Q. & Sun, Y. Community composition and species diversity of different island habitat on Karst Mountainous in Central Guizhou. Guizhou Sci. 29, 23–28 (2011).
Google Scholar 
Culver, D. Karst environment. Z. Geomorphol. Suppl. 60, 103–117 (2016).Article 

Google Scholar 
Daily, G. Restoring value to the world’s degraded lands. Science 269, 350–354 (1995).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bartgis, R. The endangered sedge Scirpus ancistrochaetus and the flora of sinkhole ponds in Maryland and West Virginia. Castanea 57, 46–51 (1992).
Google Scholar 
Yu, X., Li, Y. & Ma, Z. A preliminary study on flora diversity of karst microhabitat in Shilin Park, Yunnan, China. J. Mountain Sci. 25, 438–447 (2007).
Google Scholar 
Dang, G., Feng, H., Tang, Q., Mo, F. & Xue, Y. New recorded species in Guangxi, China. J. Guangxi Normal Univ. (Nat. Sci. Edit.) 34, 147–150 (2016).
Google Scholar 
Eigenbrod, F., Gonzalez, P., Dash, J. & Steyl, I. Vulnerability of ecosystems to climate change moderated by habitat intactness. Glob. Change Biol. 21, 275–286 (2015).ADS 
Article 

Google Scholar 
Cornell, H. & Lawton, J. Species interactions, local and regional processes, and limits to the richness of ecological communities: A theoretical perspective. J. Anim. Ecol. 61, 1–12 (1992).Article 

Google Scholar 
Helmus, M., Mahler, D. & Losos, J. Island biogeography of the Anthropocene. Nature 513, 543–546 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Yuan, T., Zhang, H., Ou, Z. & Tan, Y. Effects of topography on the diversity and distribution pattern of ground plants in karst montane forests in Southwest Guangxi, China. Chin. J. Appl. Ecol. 25, 2803–2810 (2014).
Google Scholar 
Wen, L. et al. The succession characteristics and its driving mechanisms of plant community in karst region, Southwest China. Acta Ecol. Sin. 35, 5822–5833 (2015).Article 

Google Scholar 
Zhang, Z., Hu, G. & Ni, J. Erratum to: Effects of topographical and edaphic factors on the distribution of plant communities in two subtropical Karst forests, Southwestern China. J. Mt. Sci. 10, 337–338 (2013).Article 

Google Scholar 
Du, H. et al. Plant community characteristics and its coupling relationships with soil in depressions between karst hills, North Guangxi, China. Chin. J. Plant Ecol. 37, 197–208 (2013).Article 

Google Scholar 
Liu, S., Zhang, B., Yang, Q., Hu, C. & Su, C. Species composition and diversity of plant communities in Xiaoyanwan Garden of Xingwen Karst National Geopark, Sichuan Province. Subtrop. Plant Sci. 38, 37–40 (2009).
Google Scholar 
Tan, C. The preliminary discussion about Haifeng’s wetland ecosystem. For. Sci. Technol. 1–8 (2002).Ma, K. & Liu, Y. Measurement of biodiversity: The measurement of α diversity. Chin. Biodivers. 2, 231–239 (1995).
Google Scholar 
Pielou, E. The measurement of diversity in different types of biological collections. J. Theor. Biol. 13, 131–144 (1966).ADS 
Article 

Google Scholar 
Simpson, E. Measurement of diversity. Nature 163, 688 (1949).ADS 
MATH 
Article 

Google Scholar  More

Study areaThe study was conducted on two rivers in north-eastern Amazonia sensu lato, including the Guiana Shield and the Amazon River drainage (Fig. 2). The climate of the entire study area is homogeneous and the region is covered by dense, uniform lowland primary rainforest51. The altitude is in the range of 0–860 m a.s.l. The regional climate is equatorial, and the annual rainfall ranges from 3600 mm in the northeast to 2000 mm in the southwest. The Maroni River is 612 km long from its source to its estuary, and its watershed covers a surface of >68,000 km2 in Suriname and French Guiana. The Oyapock River (length, 404 km; area, 26,800 km2) is located in the state of Amapa in Brazil and in French Guiana.The foregoing river basins host nearly 400 freshwater fish species and more than 180 mammal species52,53. Most of the mammal species have a large distribution range, covering the entire study area53. The fish species have a less homogeneous distribution and a distinct upstream-downstream composition gradient54,55. Here, only large rivers were considered and most fish species were widespread over the whole study area. As habitat availability increases with river size, species richness is expected to increase upstream to dowsntream31,32. The Oyapock and Maroni river basins are among the last remaining wilderness areas on Earth17. Nevertheless, ecological disturbances are increasing there because of a growing human population and the development of small-scale gold mining activity. These disturbances have caused limited but diffuse deforestation23,56. The deforested areas currently comprise 0.67% of all Maroni and Oyapock catchments.SamplingEnvironmental DNA (eDNA) was collected from water samples at 74 locations (hereafter, sites) along the main channel and the large tributaries of the Maroni and Oyapock rivers (Fig. 2). Thirty-seven sites were sampled at each river basin. The minimum and maximum distances between adjacent sites were 1.07 and 50.20 km, respectively. The mean and median distances between adjacent sites were 10.18 and 9.14 km, respectively, and the standard deviation (SD) was 7.79 km. The sites were located from sea level to 157 m a.s.l. At all sites, the river was wider than 20 m and deeper than 1 m (Strahler orders 4–8; Supplementary Fig. 5). The physicochemical properties of the water slightly varied among sites. The temperature, pH, and conductivity were in the ranges of 28.4–33.2 °C, 6.5–7.6, and 16.9–54.6 µS/cm, respectively, at all sites except two estuarine locations where the conductivity was relatively high because of seawater incursion (Supplementary Data 2).The eDNA samples were collected during the dry seasons (October–November) of 2017 and 2018 for Maroni and Oyapock, respectively. At both rivers, the sites were sequentially sampled from downstream to upstream at a rate of 1–4 sites per day depending on the distance and travel time between sites. Following the protocol of ref. 45, we collected the eDNA by filtering two replicates of 34 L of water per site. A peristaltic pump (Vampire Sampler; Buerkle GmbH, Bad Bellingen, Germany) and single-use tubing were used to pump the water into a single-use filtration capsule (VigiDNA, pore size 0.45 μm; filtration surface 500 cm2, SPYGEN, Bourget-du-Lac, France). The tubing input was placed a few centimetres below the water surface in zones with high water flow as recommended by Cilleros et al.43. Sampling was performed in turbulent areas with rapid hydromorphologic units to ensure optimal eDNA homogeneity throughout the water column. To avoid eDNA cross-contamination among sites, the operator remained on emerging rocks downstream from the filtration area. At the end of filtration, the capsule was voided, filled with 80 mL CL1 preservation buffer (SPYGEN), and stored in the dark up to one month before the DNA extraction. No permits were required for the eDNA sampling and the access to all sites was legally permitted. The study complies with access and benefit permits ABSCH-IRCC-FR-246820-1 and ABSCH-IRCC-FR-245902-1, authorizing collection, transport and analysis of all environmental DNA samples used in this study.Laboratory procedures and bioinformatic analysesFor the DNA extraction, each filtration capsule was agitated on an S50 shaker (Ingenieurbüro CAT M. Zipperer GmbH, Ballrechten-Dottingen, Germany) at 800 rpm for 15 min, decanted into a 50 mL tube, and centrifuged at 15,000 × g and 6 °C for 15 min. The supernatant was removed with a sterile pipette, leaving 15 mL of liquid at the bottom of the tube. Subsequently, 33 mL of ethanol and 1.5 mL of 3 M sodium acetate were added to each 50 mL tube, and the mixtures were stored at −20 °C for at least one night. The tubes were then centrifuged at 15,000 × g and 6 °C for 15 min, and the supernatants were discarded. Then, 720 µL of ATL buffer from a DNeasy Blood & Tissue Extraction Kit (Qiagen, Hilden, Germany) was added. The tubes were vortexed, and the supernatants were transferred to 2 mL tubes containing 20 µL proteinase K. The tubes were then incubated at 56 °C for 2 h. DNA extraction was performed using a NucleoSpin Soil kit (Macherey-Nagel GmbH, Düren, Germany) starting from step six of the manufacturer’s instructions. Elution was performed by adding 100 µL of SE buffer twice. After the DNA extraction, the samples were tested for inhibition by qPCR following the protocol in ref. 57. Briefly, quantitative PCR was performed in duplicate for each sample. If at least one of the replicates showed a different Ct (Cycle threshold) than expected (at least 2 Cts), the sample was considered inhibited and diluted 5-fold before the amplification.For the fish, the “teleo” primers58 (forward: 3ʹ-ACACCGCCCGTCACTCT-5ʹ; reverse: 3ʹ-CTTCCGGTACACTTACCATG-5ʹ) were used as they efficiently discriminated local fish species43,45. For the mammals, the 12S-V5 vertebrate marker59 (forward: 3ʹ-TAGAACAGGCTCCTCTAG-5ʹ; reverse: 3ʹ-TTAGATACCCCACTATGC-5ʹ) was used as it also effectively distinguishes local mammal species44,60. The DNA amplifications were performed in a final volume of 25 μL containing 1 U AmpliTaq Gold DNA Polymerase (Applied Biosystems, Foster City, CA, USA), 0.2 μM of each primer, 10 mM Tris-HCl, 50 mM KCl, 2.5 mM MgCl2, 0.2 mM of each dNTP, and 3 μL DNA template. Human blocking primer was added to the mixture for the “teleo”58 (5′-ACCCTCCTCAAGTATACTTCAAAGGAC-C3-3′) and the “12S-V5” primers61 (5′-CTATGCTTAGCCCTAAACCTCAACAGTTAAATCAACAAAACTGCT-C3-3′) at final concentrations of 4 μM and 0.2 μg/μL bovine serum albumin (BSA; Roche Diagnostics, Basel, Switzerland). Twelve PCR replicates were performed per field sample. The forward and reverse primer tags were identical within each PCR replicate. The PCR mixture was denatured at 95 °C for 10 min, followed by 50 cycles of 30 s at 95 °C, 30 s at 55 °C for the “teleo” primers and 50 °C for the 12S-V5 primers, 1 min at 72 °C, and a final elongation step at 72 °C for 7 min. This step was conducted in a dedicated room for DNA amplification that is kept under negative air pressure and is physically separated from the DNA extraction rooms maintained under positive air pressure. The purified PCR products were pooled in equal volumes to achieve an expected sequencing depth of 500,000 reads per sample before DNA library preparation.For the fish analyses, 10 libraries were prepared using a PCR-free library protocol (https://www.fasteris.com/metafast) at Fasteris, Geneva, Switzerland. Four libraries were sequenced on an Illumina HiSeq 2500 (2 × 125 bp) (Illumina, San Diego, CA, USA) with a HiSeq SBS Kit v4 (Illumina), three were sequenced on a MiSeq (2 × 125 bp) (Illumina) with a MiSeq Flow Cell Kit Version3 (Illumina), and three libraries were sequenced on a NextSeq (2 × 150 bp + 8) (Illumina) with a NextSeq Mid kit (Illumina). The libraries run on the NextSeq were equally distributed in four lanes. Sequencing was performed according to the manufacturer’s instructions at Fasteris. For the mammal analyses, eight libraries were prepared using a PCR-free library protocol (https://www.fasteris.com/metafast) at Fasteris. Two libraries were sequenced on an Illumina HiSeq 2500 (2 × 125 bp) (Illumina) using a HiSeq Rapid Flow Cell v2 and a HiSeq Rapid SBS Kit v2 (Illumina), three libraries were prepared on a MiSeq (2 × 125 bp) (Illumina) with a MiSeq Flow Cell Kit Version3 (Illumina), and three libraries were prepared using a NextSeq (2 × 150 bp + 8) (Illumina) and a NextSeq Mid kit (Illumina). The libraries run on the NextSeq were equally distributed in four lanes. As different sequencing platforms were used (MiSeq and NextSeq for the Maroni and HiSeq 2500 and MiSeq for the Oyapock; Supplementary Fig. 6 and Supplementary Data 3), the possible influences of the platforms on the sequencing results were verified. To this end, we compared the differences in species numbers between the sample replicates assigned to the same platform (accounting for replicate effect only) against those of the sample replicates assigned to different platforms (accounting for replicate and platform effects). As there were more sites with their two replicates sequenced with the same platform than sites with their replicates sequenced with different platforms (see Supplementary Fig. 6), sites with replicates on the same platform were randomly selected for the comparisons. We repeated this procedure 50 times. The number of species between replicates sequenced on the same platform and those sequenced on different platforms did not differ for >98.5% of all fish and mammal samples (Supplementary Fig. 7 and Supplementary Note 2). Similar to these results, a previous study on 16 S rRNA amplicon has shown that the samples were not influenced by the Illumina sequencing platform used62.To monitor for contaminants, 13 negative extraction controls were performed for each of the primers (“teleo” and “12S-V5”); one control was amplified twice. All of them were amplified and sequenced by the same methods as the samples and in parallel to them. Therefore, for the negative extraction controls, 168 amplifications were prepared with the “teleo” primers (13 negative controls; one amplified and sequenced twice) and 156 amplifications with the “12S-V5” primers (13 negative controls). Fourteen negative PCR controls (ultrapure water; 12 replicates) were amplified and sequenced in parallel to the samples. Eight were amplified with the “teleo” primers and six were amplified with the “12S-V05” primers. Thus, for the PCR negative controls, there were 96 amplifications with the “teleo” primers and 72 amplifications with the Vert01 primers. Sequencing information for the controls is shown in Supplementary Data 3c.An updated version of the reference database from ref. 43 was used. There were 265 Guianese species for the fish analyses (ref. 47). The GenBank nucleotide database was consulted, but it contained little information on the Guianese fish species. Most of the sequences were derived from ref. 43. For the mammal analyses, the vertebrate database was built using ecoPCR software63 from the releases 134 and 138 of the European Nucleotide Archive (ENA), for the Maroni and Oyapock river samples, respectively. The two releases were compared, and it was established that the new mammal species added to each version did not originate from French Guiana. Hence, the results were not influenced by the EMBL release number. The relevant metabarcoding fragment was extracted from this database with ecoPCR63 and OBITools64. Therefore, the reference database comprised the local database of French Guianese mammals60, which references 576 specimens from 164 species as well as all available vertebrate species in EMBL.The sequence reads were analyzed with the OBITools package according to the protocol described by Valentini et al.58. Briefly, the forward and reverse reads were assembled with the illuminapairedend programme. The ngsfilter programme was then used to assign the sequences to each sample. A separate dataset was created for each sample by splitting the original dataset into several files with obisplit. Sequences shorter than 20 bp or occurring less than 10 times per sample were discarded. The obiclean program was used to identify amplicon sequence variants (ASVs) that have likely arisen due to PCR or sequencing errors. It uses the information of sequence counts and sequence similarities to classify whether a sequence is a variant (“internal”) of a more abundant (“head”) ASV64. After this step, we matched the ASV with the reference database to obtain the taxonomic assignation for each ASV. Sequences labelled by the obiclean programme as ‘internal’’ and probably corresponding to PCR errors were discarded. The ecotag programme was then used for taxonomic assignment of molecular operational taxonomic units (MOTUs). The taxonomic assignments from ecotag were corrected to avoid overconfidence in assignments. Species-level assignments were validated only for ≥98% sequence identity with the reference database. Sequences below this threshold were discarded.Measuring disturbance intensity using GIS dataIn riverine systems, the disturbances may accumulate because of hydrologic connectivity, which is the downstream transfer of matter and pollutants4. Hence, the upstream sub-basin drainage network was considered to determine the size of the upstream sub-basin affecting local biodiversity (Fig. 1). The sub-basins were delineated by applying a flow accumulation algorithm to the SRTM global 30 m digital elevation model65. Deforestation was measured over 14 upstream spatial extents with radii of 0.5, 1.5, 3, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, and 90 km for each sampling site. Then, these 14 upstream spatial extents were intersected with the sub-basin drainage network. In addition, mammals and fish can also be affected by disturbances other than those mediated by hydrologic connectivity. Thus, deforestation was also measured upstream and downstream from the eDNA sampling sites using the same foregoing 14 radii.At each sampling site, deforestation intensity was quantified for each of the 14 spatial extents. We summed upstream (only accounting for disturbances mediated by river hydrologic connectivity) or upstream and downstream (not only considering disturbances mediated by hydrologic connectivity) deforested surfaces from Landsat satellite image datasets. Forest loss surfaces were obtained from the Global Forest Change dataset66. The Global Forest Change dataset identifies areas deforested between 2001 and 2017 on a 30 m spatial scale. To incorporate deforested areas prior to 2000, tree canopy cover data for that year were also used. Except for river courses, all pixels with More

We express our gratitude to Lukáš Falteisek, Richard Dorrell, Jan Petrášek, Stanislav Volsobě, Kateřina Schwarzerová and Jana Krtková for constructive discussions. English has been kindly corrected by William Bourland. Furthermore, we thank to Dovilė Barcytė, William Bourland, Antonio Calado, Dora Čertnerová, Yana Eglit, Ivan Fiala, Martina Hálová, Miroslav Hyliš, Dagmar Jirsová, Petr Kaštánek, Viktorie Kolátková, Alena Kubátová, Alexander Kudryavtsev, Frederik Leliaert, Julius Lukeš, Jan Mach, Joost Mansour, Jan Mourek, Yvonne Němcová, Fabrice Not, Vladimír Scholtz, Alastair Simpson, Pavel Škaloud, Jan Šťastný, Róbert Šuťák, Daria Tashyreva, Dana Savická, Jan Šobotník, Zdeněk Verner, Jan Votýpka for kindly providing cultures and taxonomic identifications. More

Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237–240 (1998).Article 

Google Scholar 
Guidi, L. et al. Plankton networks driving carbon export in the oligotrophic ocean. Nature 532, 465–470 (2016).Article 

Google Scholar 
Ptacnik, R. et al. Diversity predicts stability and resource use efficiency in natural phytoplankton communities. Proc. Natl.Acad. Sci. USA 105, 5134–5138 (2008).Article 

Google Scholar 
Corcoran, A. A. & Boeing, W. J. Biodiversity increases the productivity and stability of phytoplankton communities. PLoS ONE 7, e49397 (2012).Article 

Google Scholar 
Arteaga, L., Pahlow, M. & Oschlies, A. Global patterns of phytoplankton nutrient and light colimitation inferred from an optimality-based model. Glob. Biogeochem. Cycles 28, 648–661 (2014).Article 

Google Scholar 
Lewis, M., Hebert, D., Harrison, W. G., Platt, T. & Oakey, N. S. Vertical nitrate fluxes in the oligotrophic ocean. Science 234, 870–873 (1986).Article 

Google Scholar 
McGillicuddy, D. J. J. et al. Eddy/wind interactions stimulate extraordinary mid-ocean plankton blooms. Science 316, 1021–1026 (2007).Article 

Google Scholar 
Duce, R. A. et al. Impacts of atmospheric anthropogenic nitrogen on the open ocean. Science 320, 893–897 (2008).Article 

Google Scholar 
Tang, W. et al. Revisiting the distribution of oceanic N2 fixation and estimating diazotrophic contribution to marine production. Nat. Commun. 10, 831 (2019).Article 

Google Scholar 
Letscher, R. T., Primeau, F. & Moore, J. K. Nutrient budgets in the subtropical ocean gyres dominated by lateral transport. Nat. Geosci. 9, 815–819 (2016).Article 

Google Scholar 
Righetti, D., Vogt, M., Gruber, N., Psomas, A. & Zimmermann, N. E. Global pattern of phytoplankton diversity driven by temperature and environmental variability. Sci. Adv. 5, eaau6253 (2019).Article 

Google Scholar 
Ibarbalz, F. M. et al. Global trends in marine plankton diversity across kingdoms of life. Cell 179, 1084–1097 (2019).Article 

Google Scholar 
Lévy, M., Franks, P. J. S. & Smith, K. S. The role of submesoscale currents in structuring marine ecosystems. Nat. Commun. 9, 4758 (2018).Article 

Google Scholar 
Dutkiewicz, S. et al. Dimensions of marine phytoplankton diversity. Biogeosciences 17, 609–634 (2020).Article 

Google Scholar 
Gove, J. M. et al. Near-island biological hotspots in barren ocean basins. Nat. Commun. 7, 10581 (2016).Article 

Google Scholar 
Doty, M. S. & Oguri, M. The island mass effect. ICES J. Mar. Sci. 22, 33–37 (1956).Article 

Google Scholar 
Bell, J. D. et al. Planning the use of fish for food security in the Pacific. Mar. Policy 33, 64–76 (2009).Article 

Google Scholar 
Bakker, D. C., Nielsdóttir, M. C., Morris, P. J., Venables, H. J. & Watson, A. J. The island mass effect and biological carbon uptake for the subantarctic Crozet Archipelago. Deep Sea Res. Pt II 54, 2174–2190 (2007).Article 

Google Scholar 
Heywood, K. J., Stevens, D. P. & Bigg, G. R. Eddy formation behind the tropical island of Aldabra. Deep Sea Res. Pt I 43, 555–578 (1996).Article 

Google Scholar 
Palacios, D. M. Factors influencing the island-mass effect of the Galapagos archipelago. Geophys. Res. Lett. 29, 2134 (2002).Article 

Google Scholar 
Gilmartin, M. & Revelante, N. The ‘island mass’ effect on the phytoplankton and primary production of the Hawaiian Islands. J. Exp. Mar. Biol. Ecol. 16, 181–204 (1974).Article 

Google Scholar 
Signorini, S. C., McClain, C. R. & Dandonneau, Y. Mixing and phytoplankton bloom in the wake of the Marquesas Islands. Geophys. Res. Lett. 26, 3121–3124 (1999).Article 

Google Scholar 
Messié, M., Radenac, M.-H., Lefèvre, J. & Marchesiello, P. Chlorophyll bloom in the western Pacific at the end of the 1997-98 El Niño: the role of the Kiribati Islands. Geophys. Res. Lett. 33, L14601 (2006).Article 

Google Scholar 
Messié, M. & Radenac, M.-H. Seasonal variability of the surface chlorophyll in the western tropical Pacific from SeaWiFS data. Deep Sea Res. Pt I 53, 1581–1600 (2006).Article 

Google Scholar 
Le Borgne, R., Dandonneau, Y. & Lemasson, L. The problem of the island mass effect on chlorophyll and zooplankton standing crops around Mare (Loyalty Islands) and New Caledonia. Bull. Mar. Sci. 37, 450–459 (1985).
Google Scholar 
Messié, M. et al. The delayed island mass effect: how islands can remotely trigger blooms in the oligotrophic ocean. Geophys. Res. Lett. 47, e2019GL085282 (2020).Article 

Google Scholar 
Dandonneau, Y. & Charpy, L. An empirical approach to the island mass effect in the south tropical Pacific based on sea surface chlorophyll concentrations. Deep Sea Res. Pt A 32, 707–721 (1985).Article 

Google Scholar 
Shiozaki, T., Kodama, T. & Furuya, K. Large-scale impact of the island mass effect through nitrogen fixation in the western South Pacific Ocean. Geophys. Res. Lett. 41, 2907–2913 (2014).Article 

Google Scholar 
Caputi, L. et al. Community-level responses to iron availability in open ocean plankton ecosystems. Glob. Biogeochem. Cycles 33, 391–419 (2019).Article 

Google Scholar 
Martinez, E., Rodier, M., Pagano, M. & Sauzède, R. Plankton spatial variability within the Marquesas archipelago, South Pacific. J. Mar. Syst. 212, 103432 (2020).Article 

Google Scholar 
Behrenfeld, M. J. & Falkowski, P. G. Photosynthetic rates derived from satellite-based chlorophyll concentration. Limnol. Oceanogr. 42, 1–20 (1997).Article 

Google Scholar 
Laws, E. A., Ducklow, H. & McCarthy, J. J. Temperature effects on export production in the open ocean. Glob. Biogeochem. Cycles 14, 1231–1246 (2000).Article 

Google Scholar 
Messié, M. & Chavez, F. P. A global analysis of ENSO synchrony: the oceans’ biological response to physical forcing. J. Geophys. Res. 117, C09001 (2012).
Google Scholar 
Luo, Y.-W., Lima, I. D., Karl, D. M., Deutsch, C. A. & Doney, S. C. Data-based assessment of environmental controls on global marine nitrogen fixation. Biogeosciences 11, 691–708 (2014).Article 

Google Scholar 
Messié, M. & Chavez, F. P. Seasonal regulation of primary production in eastern boundary upwelling systems. Prog. Oceanogr. 134, 1–18 (2015).Article 

Google Scholar 
Mouw, C. B. et al. A consumer’s guide to satellite remote sensing of multiple phytoplankton groups in the global ocean. Front. Mar. Sci. 4, 41 (2017).Article 

Google Scholar 
Alvain, S., Moulin, C., Dandonneau, Y. & Bréon, F. M. Remote sensing of phytoplankton groups in case 1 waters from global SeaWiFS imagery. Deep Sea Res. Pt I 52, 1989–2004 (2005).Article 

Google Scholar 
Rêve-Lamarche, A.-H. et al. Ocean color radiance anomalies in the North Sea. Front. Mar. Sci. https://doi.org/10.3389/fmars.2017.00408 (2017).Alvain, S., Loisel, H. & Dessailly, D. Theoretical analysis of ocean color radiances anomalies and implications for phytoplankton groups detection in case 1 waters. Opt. Express 20, 1070–1083 (2012).Article 

Google Scholar 
Mackey, D. J., Blanchot, J., Higgins, H. W. & Neveux, J. Phytoplankton abundances and community structure in the equatorial Pacific. Deep Sea Res. Pt II 49, 2561–2582 (2002).Article 

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

Google Scholar 
Martiny, A. C., Kathuria, S. & Berube, P. M. Widespread metabolic potential for nitrite and nitrate assimilation among Prochlorococcus ecotypes. Proc. Natl. Acad. Sci. USA 106, 10787–10792 (2009).Article 

Google Scholar 
Vallina, S. M. et al. Global relationship between phytoplankton diversity and productivity in the ocean. Nat. Commun. 5, 4299 (2014).Article 

Google Scholar 
Dai, S. et al. The seamount effect on phytoplankton in the tropical western Pacific. Mar. Environ. Res. 162, 105094 (2020).Article 

Google Scholar 
Leitner, A. B., Neuheimer, A. B. & Drazen, J. C. Evidence for long-term seamount-induced chlorophyll enhancements. Sci. Rep. 10, 12729 (2020).Article 

Google Scholar 
Bowen, B. W., Rocha, L. A., Toonen, R. J. & Karl, S. A. The origins of tropical marine biodiversity. Trends Ecol. Evol. 28, 359–366 (2013).Article 

Google Scholar 
Worm, B., Lotze, H. K. & Myers, R. A. Predator diversity hotspots in the blue ocean. Proc. Natl. Acad. Sci. USA 100, 9884–9888 (2003).Article 

Google Scholar 
Block, B. A. et al. Tracking apex marine predator movements in a dynamic ocean. Nature 475, 86–90 (2011).Article 

Google Scholar 
Harrison, A.-L. et al. The political biogeography of migratory marine predators. Nat. Ecol. Evol. 2, 1571–1578 (2018).Article 

Google Scholar 
Pompa, S., Ehrlich, P. R. & Ceballos, G. Global distribution and conservation of marine mammals. Proc. Natl. Acad. Sci. USA 108, 13600–13605 (2011).Article 

Google Scholar 
Wessel, P. & Smith, W. H. F. A global, self-consistent, hierarchical, high-resolution shoreline database. J. Geophys. Res. 101, 8741––8743 (1996).Article 

Google Scholar 
Nunn, P. D., Kumar, L., Eliot, I. & McLean, R. F. Classifying Pacific islands. Geosci. Lett 3, 7 (2016).Article 

Google Scholar 
Hasegawa, D., Lewis, M. R. & Gangopadhyay, A. How islands cause phytoplankton to bloom in their wakes. Geophys. Res. Lett. 36, L20605 (2009).Article 

Google Scholar 
Platt, T. & Sathyendranath, S. Oceanic primary production: estimation by remote sensing at local and regional scales. Science 241, 1613–1620 (1988).Article 

Google Scholar 
Hasegawa, D., Yamazaki, H., Ishimaru, T., Nagashima, H. & Koike, Y. Apparent phytoplankton bloom due to island mass effect. J. Mar. Syst. 69, 238–246 (2008).Article 

Google Scholar 
Silsbe, G. M., Behrenfeld, M. J., Halsey, K. H., Milligan, A. J. & Westberry, T. K. The CAFE model: a net production model for global ocean phytoplankton. Glob. Biogeochem. Cycles 30, 1756–1777 (2016).Article 

Google Scholar 
Ben Mustapha, Z., Alvain, S., Jamet, C., Loisel, H. & Dessailly, D. Automatic classification of water-leaving radiance anomalies from global SeaWiFS imagery: application to the detection of phytoplankton groups in open ocean waters. Remote Sens. Environ. 146, 97–112 (2014).Article 

Google Scholar 
Alvain, S., Moulin, C., Dandonneau, Y. & Loisel, H. Seasonal distribution and succession of dominant phytoplankton groups in the global ocean: a satellite view. Glob. Biogeochem. Cycles 22, GB3001 (2008).Article 

Google Scholar 
Bray, J. R. & Curtis, J. T. An ordination of the upland forest communities of Southern Wisconsin. Ecol. Monogr. 27, 325–349 (1957).Article 

Google Scholar 
Pielou, E. The measurement of diversity in different types of biological collections. J. Theor. Biol. 13, 131–144 (1966).Article 

Google Scholar 
Shannon, C. E. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).Article 

Google Scholar 
Colwell, R. K., Mao, C. X. & Chang, J. Interpolating, extrapolating, and comparing incidence-based species accumulation curves. Ecology 85, 2717–2727 (2004).Article 

Google Scholar 
De Monte, S., Soccodato, A., Alvain, S. & d’Ovidio, F. Can we detect oceanic biodiversity hotspots from space? ISME J. 7, 2054–2056 (2013).Article 

Google Scholar 
Soccodato, A. et al. Estimating planktonic diversity through spatial dominance patterns in a model ocean. Mar. Geonom. 29, 9–17 (2016).Article 

Google Scholar 
Messié, M., Petrenko, A., Doglioli, A., Martinez, E. & Alvain, S. Data from: Basin-scale biogeochemical and ecological impacts of islands in the tropical Pacific Ocean (v1.0.0). Zenodo https://doi.org/10.5281/zenodo.6416130 (2022).Messié, M. Code for: Basin-scale biogeochemical and ecological impacts of islands in the tropical Pacific Ocean (v1.0.0). Zenodo https://doi.org/10.5281/zenodo.6494328 (2022). More

Curtis, P. G., Slay, C. M., Harris, N. L., Tyukavina, A. & Hansen, M. C. Classifying drivers of global forest loss. Science 361, 1108–1111 (2018).CAS 
Article 

Google Scholar 
Gibbs, H. K. et al. Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proc. Natl Acad. Sci. USA 107, 16732–16737 (2010).CAS 
Article 

Google Scholar 
Payn, T. et al. Changes in planted forests and future global implications. Ecol. Manag. 352, 57–67 (2015).Article 

Google Scholar 
Pendrill, F., Persson, U. M., Godar, J. & Kastner, T. Deforestation displaced: trade in forest-risk commodities and the prospects for a global forest transition. Environ. Res. Lett. 14, 055003 (2019).Article 

Google Scholar 
Hurni, K. & Fox, J. The expansion of tree-based boom crops in mainland Southeast Asia: 2001 to 2014. J. Land Use Sci. 13, 198–219 (2018).Article 

Google Scholar 
Vijay, V. et al. The impacts of oil palm on recent deforestation and biodiversity loss. PLoS ONE 11, e0159668 (2016).Heilmayr, R., Echeverría, C. & Lambin, E. F. Impacts of Chilean forest subsidies on forest cover, carbon and biodiversity. Nat. Sustain. 3, 701–709 (2020).Article 

Google Scholar 
le Maire, G., Dupuy, S., Nouvellon, Y., Loos, R. A. & Hakamada, R. Mapping short-rotation plantations at regional scale using MODIS time series: case of eucalypt plantations in Brazil. Remote Sens. Environ. 152, 136–149 (2014).Article 

Google Scholar 
Wang, M. M. H., Carrasco, L. R. & Edwards, D. P. Reconciling rubber expansion with biodiversity conservation. Curr. Biol. 30, 3825–3832 (2020).CAS 
Article 

Google Scholar 
Lewis, S. L., Wheeler, C. E., Mitchard, E. T. A. & Koch, A. Restoring natural forests is the best way to remove atmospheric carbon. Nature 568, 25–28 (2019).CAS 
Article 

Google Scholar 
Dave, R. et al. Second Bonn Challenge Progress Report: Application of the Barometer in 2018 (IUCN, 2019).Sloan, S., Meyfroidt, P., Rudel, T. K., Bongers, F. & Chazdon, R. The forest transformation: planted tree cover and regional dynamics of tree gains and losses. Glob. Environ. Change 59, 101988 (2019).Article 

Google Scholar 
Petersen, R. et al. Mapping Tree Plantations with Multispectral Imagery: Preliminary Results for Seven Tropical Countries (WRI, 2016).Erb, K.-H. et al. Land management: data availability and process understanding for global change studies. Glob. Change Biol. 23, 512–533 (2017).Article 

Google Scholar 
Souza, C. M. et al. Reconstructing three decades of land use and land cover changes in Brazilian biomes with Landsat Archive and Earth Engine. Remote Sens. 12, 2735 (2020).Article 

Google Scholar 
Miettinen, J. et al. Extent of industrial plantations on Southeast Asian peatlands in 2010 with analysis of historical expansion and future projections. GCB Bioenergy 4, 908–918 (2012).Article 

Google Scholar 
Vancutsem, C. et al. Long-term (1990–2019) monitoring of forest cover changes in the humid tropics. Sci. Adv. 7, eabe1603 (2021).Puyravaud, J.-P., Davidar, P. & Laurance, W. F. Cryptic destruction of India’s native forests. Conserv. Lett. 3, 390–394 (2010).Article 

Google Scholar 
Fagan, M. E. et al. Mapping pine plantations in the southeastern U.S. using structural, spectral, and temporal remote sensing data. Remote Sens. Environ. 216, 415–426 (2018).Article 

Google Scholar 
Tropek, R. et al. Comment on “High-resolution global maps of 21st-century forest cover change”. Science 344, 981 (2014).CAS 
Article 

Google Scholar 
Global Forest Resources Assessment 2020 (FAO, 2020).FAOSTAT Agricultural Statistics Database (FAO, 2019); http://faostat.fao.org/site/291/default.aspxCook-Patton, S. C. et al. Mapping carbon accumulation potential from global natural forest regrowth. Nature 585, 545–550 (2020).CAS 
Article 

Google Scholar 
Hurni, K., Schneider, A., Heinimann, A., Nong, D. H. & Fox, J. Mapping the expansion of boom crops in mainland Southeast Asia using dense time stacks of Landsat data. Remote Sens. 9, 320 (2017).Article 

Google Scholar 
Miettinen, J., Shi, C. & Liew, S. C. 2015 Land cover map of Southeast Asia at 250 m spatial resolution. Remote Sens. Lett. 7, 701–710 (2016).Article 

Google Scholar 
Torbick, N., Ledoux, L., Salas, W. & M. Zhao, M. Regional mapping of plantation extent using multisensor imagery. Remote Sens. 8, 236 (2016).Azizan, F. A., Kiloes, A. M., Astuti, I. S. & Abdul Aziz, A. Application of optical remote sensing in rubber plantations: a systematic review. Remote Sens. 13, 429 (2021).Article 

Google Scholar 
Bégué, A. et al. Remote sensing and cropping practices: a review. Remote Sens. 10, 99 (2018).Article 

Google Scholar 
Bey, A. & Meyfroidt, P. Improved land monitoring to assess large-scale tree plantation expansion and trajectories in Northern Mozambique. Environ. Res. Commun. 3, 115009 (2021).Jucker, T. et al. Topography shapes the structure, composition and function of tropical forest landscapes. Ecol. Lett. 21, 989–1000 (2018).Article 

Google Scholar 
Féret, J.-B. & Asner, G. P. Spectroscopic classification of tropical forest species using radiative transfer modeling. Remote Sens. Environ. 115, 2415–2422 (2011).Article 

Google Scholar 
Poortinga, A. et al. Mapping plantations in Myanmar by fusing Landsat-8, Sentinel-2 and Sentinel-1 data along with systematic error quantification. Remote Sens. 11, 831 (2019).Article 

Google Scholar 
Gutiérrez-Vélez, V. H. et al. High-yield oil palm expansion spares land at the expense of forests in the Peruvian Amazon. Environ. Res. Lett. 6, 044029 (2011).Article 

Google Scholar 
Descals, A. et al. High-resolution global map of smallholder and industrial closed-canopy oil palm plantations. Earth Syst. Sci. Data. 13, 1211–1231 (2021).Article 

Google Scholar 
Ordway, E. M., Naylor, R. L., Nkongho, R. N. & Lambin, E. F. Oil palm expansion and deforestation in Southwest Cameroon associated with proliferation of informal mills. Nat. Commun. 10, 114 (2019).CAS 
Article 

Google Scholar 
Heilmayr, R., Echeverría, C., Fuentes, R. & Lambin, E. F. A plantation-dominated forest transition in Chile. Appl. Geogr. 75, 71–82 (2016).Article 

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

Google Scholar 
Bond, W. J., Stevens, N., Midgley, G. F. & Lehmann, C. E. R. The trouble with trees: afforestation plans for Africa. Trends Ecol. Evol. 34, 963–965 (2019).Article 

Google Scholar 
Veldman, J. W. et al. Where tree planting and forest expansion are bad for biodiversity and ecosystem services. Bioscience 65, 1011–1018 (2015).Article 

Google Scholar 
Venter, O. et al. Global terrestrial Human Footprint maps for 1993 and 2009. Sci. Data 3, 160067 (2016).Article 

Google Scholar 
Fagan, M. E. A lesson unlearned? Underestimating tree cover in drylands biases global restoration maps. Glob. Change Biol. 26, 4679–4690 (2020).Bastin, J. F. et al. The extent of forest in dryland biomes. Science 356, 635–638 (2017).CAS 
Article 

Google Scholar 
Fagan, M. E., Reid, J. L., Holland, M. B., Drew, J. G. & Zahawi, R. A. How feasible are global forest restoration commitments? Conserv. Lett. 13, e12700 (2020).Article 

Google Scholar 
Malkamäki, A. et al. A systematic review of the socio-economic impacts of large-scale tree plantations, worldwide. Glob. Environ. Change 53, 90–103 (2018).Article 

Google Scholar 
Schwartz, N. B., Aide, T. M., Graesser, J., Grau, H. R. & Uriarte, M. Reversals of reforestation across Latin America limit climate mitigation potential of tropical forests. Front. For. Glob. Change 3, 85 (2020).Article 

Google Scholar 
Noojipady, P. et al. Managing fire risk during drought: the influence of certification and El Niño on fire-driven forest conversion for oil palm in Southeast Asia. Earth Syst. Dynam. 8, 749–771 (2017).Bullock, E. L., Woodcock, C. E., Souza, C. Jr. & Olofsson, P. Satellite-based estimates reveal widespread forest degradation in the Amazon. Glob. Change Biol. 26, 2956–2969 (2020).Article 

Google Scholar 
Sloan, S. & Sayer, J. A. Forest Ecology and Management Forest Resources Assessment of 2015 shows positive global trends but forest loss and degradation persist in poor tropical countries. Ecol. Manag. 352, 134–145 (2015).Article 

Google Scholar 
Heinrich, V. H. A. et al. Large carbon sink potential of secondary forests in the Brazilian Amazon to mitigate climate change. Nat. Commun. 12, 1785 (2021).CAS 
Article 

Google Scholar 
Potapov, P. et al. Mapping global forest canopy height through integration of GEDI and Landsat data. Remote Sens. Environ. 253, 112165 (2021).Article 

Google Scholar 
Bernal, B., Murray, L. T. & Pearson, T. R. H. Global carbon dioxide removal rates from forest landscape restoration activities. Carbon Balance Manag. 13, 22 (2018).CAS 
Article 

Google Scholar 
Li, W., Goodchild, M. F. & Church, R. An efficient measure of compactness for two-dimensional shapes and its application in regionalization problems. Int. J. Geogr. Inf. Sci. 27, 1227–1250 (2013).Article 

Google Scholar 
Asner, G. P. Cloud cover in Landsat observations of the Brazilian Amazon. Int. J. Remote Sens. 22, 3855–3862 (2001).Article 

Google Scholar 
Wilson, A. M. & Jetz, W. Remotely sensed high-resolution global cloud dynamics for predicting ecosystem and biodiversity distributions. PLoS Biol. 14, e1002415 (2016).Article 
CAS 

Google Scholar 
Gutiérrez-Vélez, V. H. & DeFries, R. Annual multi-resolution detection of land cover conversion to oil palm in the Peruvian Amazon. Remote Sens. Environ. 129, 154–167 (2013).Article 

Google Scholar 
Reiche, J. et al. Combining satellite data for better tropical forest monitoring. Nat. Clim. Change 6, 120–122 (2016).Article 

Google Scholar 
Erinjery, J. J., Singh, M. & Kent, R. Mapping and assessment of vegetation types in the tropical rainforests of the Western Ghats using multispectral Sentinel-2 and SAR Sentinel-1 satellite imagery. Remote Sens. Environ. 216, 345–354 (2018).Article 

Google Scholar 
Shimada, M. et al. New global forest/non-forest maps from ALOS PALSAR data (2007–2010). Remote Sens. Environ. 155, 13–31 (2014).Article 

Google Scholar 
Torres, R. et al. GMES Sentinel-1 mission. Remote Sens. Environ. 120, 9–24 (2012).Article 

Google Scholar 
Potapov, P. et al. The last frontiers of wilderness: tracking loss of intact forest landscapes from 2000 to 2013. Sci. Adv. 3, e1600821 (2017).Article 

Google Scholar 
World Database on Protected Areas User Manual 1.4 (UNEP-WCMC, 2016).AutoML: Automatic Machine Learning (H2O.ai, 2020); https://h2o-release.s3.amazonaws.com/h2o/rel-yau/5/docs-website/h2o-docs/automl.htmlHealey, S. P. et al. Mapping forest change using stacked generalization: an ensemble approach. Remote Sens. Environ. 204, 717–728 (2018).Article 

Google Scholar 
Lagomasino, D. et al. Measuring mangrove carbon loss and gain in deltas. Environ. Res. Lett. 14, 25002 (2019).Article 

Google Scholar 
Bunting, P. et al. The global mangrove watch—a new 2010 global baseline of mangrove extent. Remote Sens. 10, 1669 (2018).Article 

Google Scholar 
Pickens, A. H. et al. Mapping and sampling to characterize global inland water dynamics from 1999 to 2018 with full Landsat time-series. Remote Sens. Environ. 243, 111792 (2020).Article 

Google Scholar 
Olofsson, P. et al. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 148, 42–57 (2014).Article 

Google Scholar 
Stehman, S. V. Estimating area and map accuracy for stratified random sampling when the strata are different from the map classes. Int. J. Remote Sens. 35, 4923–4939 (2014).Article 

Google Scholar 
Olofsson, P. et al. Mitigating the effects of omission errors on area and area change estimates. Remote Sens. Environ. 236, 111492 (2020).Article 

Google Scholar 
Database of Global Administrative Areas (GADM) v.3.6 (GADM, 2018); https://gadm.org/download_country_v3.htmlHijmans, R. J., Williams, E., Vennes, C. M. & Hijmans, M. R. J. Package ‘geosphere’ version 1.5-10. Spherical trigonometry (2017).Olson, D. M. et al. Terrestrial ecoregions of the world: a new map of life on Earth: a new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. Bioscience 51, 933–938 (2001).Article 

Google Scholar 
Mittermeier, R. A., Turner, W. R., Larsen, F. W., Brooks, T. M. & Gascon, C. in Biodiversity Hotspots: Distribution and Protection of Conservation Priority Areas (eds Zachos, F. E. & Habel, J. C.) 3–22 (Springer, 2011).Potapov, P. et al. The last frontiers of wilderness: tracking loss of intact forest landscapes from 2000 to 2013. Sci. Adv. 3, e1600821 (2017). More

Yuuka Fukui
Yuuka Fukui received Ph.D. degree from Keio University in 2012 under the supervision of Prof. Keiji Fujimoto. She was a JSPS research fellow (DC2) from 2010 to 2012. She joined the laboratory of Prof. Keiji Fujimoto at Keio university as a research associate in 2012 and was promoted to an assistant professor in 2017. Her research interests focus on the design and synthesis of polymeric materials (particles, membranes, porous structures) and organic-inorganic hybrid materials inspired from biological systems.About the award article: The authors reported a new technique to prepare nanoparticles from biomass-derived polymers, which will be utilized as an eco-friendly alternative to synthetic particulate plastics. Nanosized agarose gel particles were produced via sol-to-gel transition of agarose inside water nanodroplets prepared by W/O miniemulsion method. Subsequently, the water evaporation was carried out to generate xerogel nanoparticles (AgarX). The morphologies and crystal structure of AgarX were controlled by changing the pressure and temperature during the water evaporation. The resultant AgarX possessed high crystallinity and exhibited a water dispersibility and a water resistance.

Mikihiro Hayashi
Mikihiro Hayashi received his Ph.D. degree from Nagoya University (Prof. Yushu Matsushita group) in 2015. During his doctor course, he had been selected as a JSPS research fellow (DC2) and experienced researches in ESPCI Paris-Tech (Prof. Ludwik Leibler) and in Shanghai Jiao Tong University (Prof. Xinyuan Zhu). He then re-joined Ludwik Leibler’s group as a postdoc, and experienced another postdoc in Prof. Masatoshi Tokita in Tokyo institute of technology. In 2017, he became an assistant professor in Prof. Akinori Takasu group (Nagoya institute of technology), and currently manages his own laboratory as a PI. His research interest is the design of functional cross-linked materials.About the award article: the authors reported a preparation vitrimer-like elastomers with dynamic bond-exchangeable cross-links. A poly(ethyl acrylate)-based copolymer bearing random pyridine groups was synthesized, which was cross-linked by quaternization reaction with dibromo cross-linkers. In this system, the bond exchange was operated via trans-N-alkylation of the quaternized pyridine groups, showing useful sustainable functions, such as reprocessability, recyclability, and dissolution ability in some selective solvents.

Ryohei Ishige
Ryohei Ishige received his Ph.D. from Tokyo Institute of Technology in 2011 under the supervision of Prof. Junji Watanabe. He joined Prof. Atsushi Takahara’s laboratory at Kyushu University (2011–2013) and Prof. Yoshinobu Tsujii’s laboratory at Kyoto University (2013–2014). From 2014, he joined Prof. Shinji Ando’s laboratory at Tokyo Institute of Technology as an assistant professor and was promoted to an associate professor in 2021. His research interests are liquid-crystalline aromatic polymers and those structure-property relationships.About the award article: the authors developed a novel analytical technique integrating spectroscopies (infrared pMAIRS, and spectroscopic ellipsometry) and scattering methods (GI-WAXS), applied to the process where thin film polyimide, PI, is generated from linear poly(amic ester), PAE, precursors whose backbone consists of para-linkage. They revealed that PAE-based thin PI films form heterogeneous structure composed of non-oriented amorphous region and oriented ordered region which includes anisotropic nanopores causing structural birefringence. This method enables comprehensive evaluation of the evolution in complex hierarchical structures following chemical reactions for every noncrystalline thin film polymers.

Ryohei Kakuchi
Ryohei Kakuchi received his Ph.D. degree from the Hokkaido University in 2009 with a JSPS (Japan Society for Promotion of Science) research fellowship. After the Ph.D., he has made postdoctoral works in Germany from 2009 to 2014 and joined Kanazawa University as a research assistant professor in 2014. Based on the Leading Initiative for Excellent Young Researchers program, he was then appointed as an assistant professor (PI) at Gunma University in 2017. His research interests are the novel polymer synthesis based on unique organic transformation reactions including multicomponent reactions.About the award article: The authors proposed a new synthetic strategy to utilize wood-biomass sourced compounds in a green fashion. To achieve sustainable material chemistry, the intrinsic reactivity of lignin-derived poly(methacrylated vanillin) (PMV) was spotlighted because many multicomponent reactions employ aldehydes as a reactant. First, the Passerini three-component reaction (Passerini-3CR) of the PMV was revealed to proceed with >90% aldehyde conversions. Taking advantage of this high reactivity of the PMV, its immobilized cellulose fabric, a wood-biomass sourced organic hybrid, was revealed to accept the surface Passerini-3CR with amino acid derivatives, thereby demonstrating a fully bio-based material fabrication. More

Site descriptionThe respective sites were classified into five climatic regions in California, containing San Cruz and San Rose in region 1, Saint Helena and San Jose in region 2, Livermore and Cloverdale in region 3, Davis, Lodi and Fontana in region 4, Fresno and Bakerfield in region 5 (Fig. 1). There were differences in annual mean temperature among five climatic regions, ranging from 14.3°C to 18.6°C. In each region, the GHDs, quality of musts and wines, and wine tasting notes were recorded for 148 cultivars from 1935 to 1941. Meanwhile, in region 2, namely in Napa, the GHDs and must sugar content (in °Brix) were recorded for four representative cultivars (Cabernet Sauvignon, Chardonnay, Merlot and Sauvignon Blanc) during 1991–2018.Fig. 1The locations of five climatic regions for wine grape classed by Winkler index in California. The insert plot represents the distinct Winkler index (WI) during 1935–1941 in five climatic regions.Full size imageClimate dataThe climate data was collected from five stations for over one hundred year-period (1911–2018), including daily average, maximum and minimum temperature (Table 1). Climate data was retrieved from the National Oceanic and Atmospheric Administration (NOAA)’s National Centers for Environmental Information (NCEI). The database from which the data was retrieved was the “Global Historical Climatology Network – Daily (GHCN-Daily), Version 3” (https://www1.ncdc.noaa.gov/pub/data/ghcn/daily/by_station/)25,26. Table 1 showed the search codes and names of five stations in the website. The climate data of region 1 and region 5 were for the periods of 1911–2011 and 1938–2018, respectively.Table 1 Description of weather stations and time-span in five climatic regions.Full size tableBioclimatic indicesHere, we presented seven temperature-related indices to explore the changing climate in five climatic regions during the last 100 years. We compared the changes of these indices between the past (1935–1941) and current climate conditions (1991–2018). Thereafter, four indices were chosen to describe annual changes, including average, maximum, minimum temperature and diurnal temperature range (DTR). Furthermore, other indices were used to analyse growing season temperature (GST), Winkler index (WI) and Huglin index (HI) for the grape-growing season5,27,28. The equations used to calculate the bioclimatic indices of grape-growing season are:$$GST=frac{{sum }_{Apr1}^{Oct31}frac{{T}_{max}+{T}_{min}}{2}}{n}$$
(1)
$$WI={sum }_{Apr1}^{Oct31}left(frac{{T}_{max}+{T}_{min}}{2}-10right)$$
(2)
$$HI={sum }_{Apr1}^{Sep30}left(frac{{T}_{max}+{T}_{ave}}{2}-10right)times K$$
(3)
where Tmax, Tmin and Tave represent daily maximum, minimum and average temperatures, respectively. K is a length of day coefficient ranging from 1.02 to 1.06 between 40 and 50 of latitude in the northern hemisphere.Sample collection, harvest dates, quality of musts and wines measurementSample collection, harvest dates, quality of musts and wines measurement were detailed in the report of Amerine and Winkler24. Briefly, grape berries (22–220 kg) were picked in the morning from representative vines of variety collections or commercial vineyards by Amerine and Winkler24, as well as numerous vineyard owners. The harvest dates were recorded after picking. All grapes picked were crushed within 24 hours except for a few samples in 1935. The clear juice was taken after the coarse sediment had settled, in order to measure total soluble solids (°Brix), total acid (grams per 100 cc), and pH of must. The must was placed in an open oak fermenting tank. After fermentation, it was completed in a closed oak container. Then, the alcohol (percent by volume), extract (grams per 100 cc), tannin (grams per 100 cc), and fixed acid (grams per 100 cc) of wine were measured. The must °Brix was measured with a Brix hydrometer floating in a cylinder, must total acid was determined by titration with sodium hydroxide to a phenolphthalein end point, and must pH was measured with a quinhydrone electrode or a Beckman pH meter. In addition, wine alcohol was measured by the hydrometer and reported as percentage by volume, the extract and tannin of wine were measured by means of a special 0° to 8° Balling hydrometer and the Association of Official Agricultural Chemists method24. Note that the fixed acid of wine are equal to total acid minus volatile acid, where the total acid was measured by titration with phenolphthalein as an indicator while the volatile acid was determined also by titration with pretreated wines by method II of the Association of Official Agricultural Chemists24.Wine tasting notesThe purpose of wine tasting was to evaluate the cultivars based on the merits and defects of wine. The descriptive terms used for recording the results of the organoleptic examination contained appearance, color, odors, volatile acidity, total acidity, dryness, body, taste, smoothness and astringency, and general quality. More