Ecology

FAO (Food and Agricultural Organization). Food and Agricultural Organization of the United Nations. http://www.faostat.org (2020).Nouatin, G., Kougbadi, S. & Afouda, L. Analyse des contraintes de la production rizicole et les stratégies développées par les femmes de la commune de Gogounou. Ann. Sci. Agro. 12, 45–59 (2009).
Google Scholar 
Totin, E. et al. Barriers and opportunities for innovation in rice production in the inland valleys of Benin. NJAS-Wagen J. Life Sci. 60–63, 57–66 (2012).Article 

Google Scholar 
Tanaka, A., Saito, K., Azoma, K. & Kobayashi, K. Factors affecting variation in farm yields of irrigated lowland rice in southern-central Benin. Eur. J. Agron. 44, 46–53 (2013).Article 

Google Scholar 
Nonvide, G. M. A., Sarpong, D. B., Kwadzo, G.T.-M., Anim-Somuah, H. & Amoussouga, G. F. Farmers’ perceptions of irrigation and constraints on rice production in Benin: A stakeholder-consultation approach. Int. J. Water Resour. Dev. 34, 1001–1021 (2018).Article 

Google Scholar 
Seye, B., Arouna, A., Sall, S. N. & Ndiaye, A. A. Determinants de l’adoption des semences certifiees de varietes ameliorees du riz au Benin. J. Rech. Sci. Univ. Lomé 18, 93–106 (2016).
Google Scholar 
Chiambo, P. J., Coelho, J. P., Soares, F. B. & Salumbo, A. Characterization of rice production system in Camacupa and Catabola municipalities of the province of Bié in Angola. DRJAFS. 7, 250–263 (2019).
Google Scholar 
Kleinhenz, V., Chea, S. & Hun, N. Survey of rice cropping systems in Kampong Chhnang Province, Cambodia. Rice Sci. 20, 154–164 (2013).Article 

Google Scholar 
Loko, Y. L. E. et al. On-Farm Management of Rice Diversity, Varietal Preference Criteria, and Farmers’ Perceptions of the African (Oryza glaberrima Steud.) Versus Asian Rice (Oryza sativa L.) in the Republic of Benin (West Africa): Implications for Breeding and Conservation. Econ. Bot. 75, 1–29 (2021)Bello, O. L., Baiyegunhi, L. J. S. & Danso-Abbeam, G. Productivity impact of improved rice varieties’ adoption: Case of smallholder rice farmers in Nigeria. Econ. Innov. New Tech. https://doi.org/10.1080/10438599.2020.1776488 (2020).Article 

Google Scholar 
Gnacadja, C., Azokpota, P., Moreira, J. & Sie, M. Perceptions des producteurs et consommateurs sur le riz africain (Oryza glaberrima). Int. J. Biol. Chem. Sci. 11, 2778–2792 (2017).Article 

Google Scholar 
Yokouchi, T. & Saito, K. Factors affecting farmers’ adoption of NERICA upland rice varieties: The case of a seed producing village in central Benin. Food Sec. 8, 197–209 (2016).Article 

Google Scholar 
Chandio, A. A. & Yuansheng, J. Determinants of adoption of improved rice varieties in northern Sindh, Pakistan. Rice Sci. 25, 103–110 (2018).Article 

Google Scholar 
Dagnelie P. Statistiques théoriques et appliquées [Theoretical and applied statistics]. Paris, France: De Boeck & Larcier SA. (1998).Adebo, H. O. et al. Ethnobotanical Knowledge of Jute (Corchorus olitorius L.) in Benin. Eur. J. Med. Plants. 26, 1–11 (2018).INSAE (Institut National de la Statistique et de l’Analyse Économique). Principaux indicateurs sociodémographiques et économiques (RGPH-4, 2013). https://insae.bj/statistiques/enquetes-et-recensements (2016).MAEP (Ministère de l’Agriculture, de l’Elevage et de la Pêche). Plan Stratégique de Relance du Secteur Agricole. Repport, MAEP, Bénin. http://extwprlegs1.fao.org/docs/pdf/ben149176.pdf (2008).Dansi, A. et al. Traditional leafy vegetables and their use in the Benin Republic. Genet Resour Crop Evol. 55, 1239–1256 (2008).Article 

Google Scholar 
Keuls, M. The use of the “studentized range” in connection with an analysis of variance. Euphytica 1, 112–122 (1952).Article 

Google Scholar 
Xie, Y. knitr: A General-Purpose Package for Dynamic Report Generation in R. R package version 129 (2020).R Core Team. R: A Language and Environment for Statistical Computing Vienna: R Foundation for Statistical Computing. https://www.gbif.org/fr/tool/81287/r-a-language-and-environment-for-statistical-computing (2018).Fox, J. & Weisberg, S. An R Companion to Applied Regression, Third edition Sage, Thousand Oaks CA. https://socialsciences.mcmaster.ca/jfox/Books/Companion/ (2019).Oksanen, J. et al. Vegan: Community Ecology Package R package version 25-6. https://cran.r-project.org/package=vegan (2019).Robinson, D., Hayes, A. & Couch, S. Broom: Convert Statistical Objects into Tidy Tibbles R package version 070. https://cran.r-project.org/package=broom (2020).Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).Book 

Google Scholar 
Kinkingninhoun-Mêdagbé, F. M., Diagne, A., Simtowe, F., Agboh-Noameshie, A. R. & Adégbola, P. Y. Gender discrimination and its impact on income, productivity, and technical efficiency: evidence from Benin. Agric. Human Values. 27, 57–69 (2010).Article 

Google Scholar 
Adétonah S. et al. Analysis of gender and governance of value chain-based systems on rice and vegetable crops in southern Benin and Mali. Open .J Soc. Sci. 3, 134–141 (2015).Le Vido, A. A. riz africain (Oryza glaberrima Steudel) dans l’agrosystème des Fon du plateau d’Abomey (Bénin) au XIXè siècle: essai d’approche historique. Rev. Iv. Hist. 20, 59–76 (2012).
Google Scholar 
Sakurai, T. Intensification of rainfed lowland rice production in West Africa: Present status and potential green revolution. Dev. Econ. 44, 232–251 (2006).Article 

Google Scholar 
Feng, S. Land rental, off-farm employment and technical efficiency of farm households in Jiangxi Province, China. NJAS 55, 363–378 (2008).
Google Scholar 
Affholder, F., Poeydebat, C., Corbeels, M., Scopel, E. & Tittonell, P. The yield gap of major food crops in family agriculture in the tropics: Assessment and analysis through field surveys and modelling. Field Crops Res. 143, 106–118 (2013).Article 

Google Scholar 
Anang, B. T. & Awuni, J. A. Effect of training on small-scale rice production in northern Ghana. APSTRACT. 12, 13–20 (2018).Article 

Google Scholar 
Nonvide, G. M. A. A re-examination of the impact of irrigation on rice production in Benin: An application of the endogenous switching model. Kasetsart J. Soc. Sci. 40, 657–662 (2019).
Google Scholar 
Osawe, O. W., Akinyosoye, V. O., Omonona, B. T., Okoruwa, V. O. & Salman, K. K. Productivity differentials in rice production systems: Evidence from rice farmers in five agroecological zones in Nigeria. J. Nutraceut. Food Sci. 2, 1–18 (2017).
Google Scholar 
Yabi, I. & Afouda, F. Extreme rainfall years in Benin (West Africa). Quat. Int. 262, 39–43 (2012).Article 

Google Scholar 
Goulart, R. Z., Reichert, J. M. & Rodrigues, M. F. Cropping poorly-drained lowland soils: Alternatives to rice monoculture, their challenges and management strategies. Agric. Syst. 177, 102715. https://doi.org/10.1016/jagsy2019102715 (2020).Dobermann, A. & Fairhurst, T. Rice: Nutrient Disorders & Nutrient Management Handbook Series, Potash & Phosphate Institute (PPI), Potash & Phosphate Institute of Canada (PPIC) and International Rice Research Institute, Philippine, 191 (2000).Roder, W., Maniphone, S. & Keoboulapha, B. Pigeon pea for fallow improvement in slash-and-burn systems in the hills of Laos?. Agrofor. Syst. 39, 45–57 (1997).Article 

Google Scholar 
Van Campenhout, B. The role of information in agricultural technology adoption: Experimental evidence from rice farmers in Uganda. Econ. Dev. Cult. Change 69, 1239–1272 (2021).Article 

Google Scholar 
Castillo, J., Kirk, G. J. D., Rivero, M. M. J., Dobermann, A. & Haefele, M. The nitrogen economy of rice-livestock systems in Uruguay. Glob. Food Sect. 30, 100566; https://doi.org/10.1016/j.gfs.2021.100566 (2021).Paman, U., Inaba, S. & Uchida, S. The mechanization of small-scale rice farming: Labor requirements and costs. Eng. Agric. Environ. Food 7, 122–126 (2014).Article 

Google Scholar 
Mesfin, A. H. & Zemedu, L. Choices of varieties and demand for improved rice seed in Fogera district of Ethiopia. Rice Sci. 25, 350–356 (2018).Article 

Google Scholar 
Chandio, A. A., Jiang, Y., Gessesse, A. T. & Dunya, R. The nexus of agricultural credit, farm size and technical efficiency in Sindh, Pakistan: A stochastic production frontier approach. J. Saudi Soc. Agric. Sci. 18(3), 348–354 (2019).
Google Scholar 
Naseem, A., Mhlanga, S., Diagne, A., Adegbola, P. Y. & Midingoyi, G. S. Economic analysis of consumer choices based on rice attributes in the food markets of West Africa—The case of Benin. Food Sect. 5, 575–589 (2013).Article 

Google Scholar 
Demont, M., Fiamohé, R. & Kinkpé, T. Comparative advantage in demand and the development of rice value chains in West Africa. World Dev. 96, 578–590 (2017).Article 

Google Scholar 
Zannou, A., Kpenavoun, C. S,. Saliou, I. O. & Biaou, G. Technical efficiency of irrigated rice seed farmers in Koussin-Lélé, Benin Republic. J. Dev. Agric. Econ. 10, 28–37 (2018).Ouédraogo, M. & Dakouo, D. Evaluation de l’adoption des variétés de riz NERICA dans l’Ouest du Burkina Faso. Afr. J. Agric. Resour. Econ. 12, 1–16 (2017).Article 

Google Scholar 
Bruce, A. K. K., Donkoh, S. A. & Ayamga, M. Improved rice variety adoption and its effects on farmers’ output in Ghana. J. Dev. Agric. Econ. 6, 242–248 (2014).Article 

Google Scholar 
Beke, T. E. Institutional constraints and adoption of improved rice varieties: Econometric evidence from Ivory Coast. RAEStud. 92, 117–141 (2011).Article 

Google Scholar 
Hagos, A. & Zemedu, L. Determinants of improved rice varieties adoption in Fogera district of Ethiopia. Sci. Technol. Arts Res. J. 4, 221–228 (2015).Reardon, T., Stamoulis, K. & Pingali, P. Rural nonfarm employment in developing countries in an era of globalization. Agric. Econ. 37, 173–183 (2007).Article 

Google Scholar 
Khush, G. S. Modern varieties—Their real contribution to food supply and equity. GeoJourna 35, 275–284 (1995).Bannor, R. K., Kumar, G. A. K., Oppong-Kyeremeh, H. & Wongnaa, C. A. Adoption and impact of modern rice varieties on poverty in Eastern India. Rice Sci. 27, 56–66 (2020).Article 

Google Scholar 
Anik, A. R. & Salam, M. D. A. Determinants of adoption of improved onion variety in Banglades. J. Agric. Environ. Int. Dev. 109, 71–88 (2015).
Google Scholar 
Nascente, A. S. & Kromocardi, R. Genotype selection and addition of fertilizer increases grain yield in upland rice in Suriname. Acta Amazon. 47, 185–194 (2017).Article 

Google Scholar 
Hossain, M. G., Sabiruzzaman, M., Islam, S., Ohtsuki, F. & Lestrel, P. E. Effect of craniofacial measures on the cephalic index of Japanese adult female students. Anthropol. Sci. 118, 117–121 (2010).Article 

Google Scholar 
Nonvide, G. M. A. Identification of factors affecting adoption of improved rice varieties among smallholder farmers in the municipality of Malanville, Benin. J. Agric. Sci. Technol. 22, 305–316 (2020).
Google Scholar  More

Jones, H. P. et al. Invasive mammal eradication on islands results in substantial conservation gains. Proc. Natl Acad. Sci. USA 113, 4033–4038 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
DIISE. The database of island invasive species eradications, developed by island conservation, coastal conservation action laboratory UCSC, IUCN SSC Invasive Species Specialist Group. http://diise.islandconservation.org (2020).Doherty, T. S., Glen, A. S., Nimmo, D. G., Ritchie, E. G. & Dickman, C. R. Invasive predators and global biodiversity loss. Proc. Natl Acad. Sci. USA 113, 11261–11265 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Simberloff, D. et al. Yes we can! Exciting progress and prospects for controlling invasives on islands and beyond. West. North Am. Nat. 78, 942 (2018).
Google Scholar 
Kappes, P. J., Bond, A. L., Russell, J. C. & Wanless, R. M. Diagnosing and responding to causes of failure to eradicate invasive rodents. Biol. Invasions 21, 2247–2254 (2019).
Google Scholar 
Browett, S. S., O’Meara, D. B. & McDevitt, A. D. Genetic tools in the management of invasive mammals: recent trends and future perspectives. Mammal. Rev. 50, 200–210 (2020).
Google Scholar 
Burgess, B. T., Irvine, R. L., Howald, G. R. & Russello, M. A. The promise of genetics and genomics for improving invasive mammal management on islands. Front. Ecol. Evol. 9, 704809 (2021).
Google Scholar 
Abdelkrim, J., Pascal, M., Calmet, C. & Samadi, S. Importance of assessing population genetic structure before eradication of invasive species: examples from insular Norway rat populations. Conserv. Biol. 19, 1509–1518 (2005).
Google Scholar 
Adams, A. L., van Heezik, Y., Dickinson, K. J. M. & Robertson, B. C. Identifying eradication units in an invasive mammalian pest species. Biol. Invasions 16, 1481–1496 (2014).
Google Scholar 
Sjodin, B. M. F., Irvine, R. L., Ford, A. T., Howald, G. R. & Russello, M. A. Rattus population genomics across the Haida Gwaii archipelago provides a framework for guiding invasive species management. Evol. Appl. 13, 889–904 (2019).
Google Scholar 
Russello, M. A., Smith-Vidaurre, G. & Wright, T. F. In Naturalized Parrots of the World: Distribution, Ecology, and Impacts of the World’s Most Colorful Colonizers (ed. Pruett-Jones, S.) Ch. 4 (Princeton Univ. Press, 2021).Russell, J. C. et al. Survivors or reinvaders? Using genetic assignment to identify invasive pests following eradication. Biol. Invasions 12, 1747–1757 (2010).
Google Scholar 
Amos, W., Nichols, H. J., Churchyard, T. & Brooke, MdeL. Rat eradication comes within a whisker! A case study of a failed project from the South Pacific. R. Soc. Open Sci. 3, 160110 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sjodin, B. M. F., Puckett, E. E., Irvine, R. L., Munshi-South, J. & Russello, M. A. Global origins of invasive brown rats (Rattus norvegicus) in the Haida Gwaii archipelago. Biol. Invasions 23, 611–623 (2021).
Google Scholar 
Gaston, A. J., Golumbia, T., Martin, J.-L. & Sharpe, S. Lessons from the islands: introduced species and what they tell us about how ecosystems work. In Proc. from the Research Group on Introduced Species 2002 Symposium 103–116 (Canadian Wildlife Society, 2008).Parks Canada. Gwaii Haanas National Park Reserve, National Marine Conservation Area Reserve, and Haida Heritage Site. https://www.pc.gc.ca/en/pn-np/bc/gwaiihaanas/nature/faune-wildlife (2019).Stockton, S. A., Allombert, S., Gaston, A. J. & Martin, J.-L. A natural experiment on the effects of high deer densities on the native flora of coastal temperate rain forests. Biol. Conserv. 126, 118–128 (2005).
Google Scholar 
Allombert, S., Stockton, S. & Martin, J.-L. A natural experiment on the impact of overabundant deer on forest invertebrates. Conserv. Biol. 19, 1917–1929 (2005).
Google Scholar 
Allombert, S., Gaston, A. J. & Martin, J.-L. A natural experiment on the impact of overabundant deer on songbird populations. Biol. Conserv. 126, 1–13 (2005).
Google Scholar 
Chollet, S., Maillard, M., Schörghuber, J., Grayston, S. J. & Martin, J. Deer slow down litter decomposition by reducing litter quality in a temperate forest. Ecology 102, e03235 (2021).Deagle, G. Traditional west coast native medicine. Can. Fam. Physician 34, 4 (1988).
Google Scholar 
Bellis, K. X. T., Peet, R. T., Irvine, R. L., Howald, G. & Alsop, G. J. In Island Invasives: Scaling up to Meet the Challenge Ch. 3 (eds Veitch, C. R., Clout, M. N., Russell, J. C. & West, C. J.) 494–496 (IUCN, 2019).Stroh, N., Baltzinger, C. & Martin, J.-L. Deer prevent western redcedar (Thuya plicata) regeneration in old-growth forests of Haida Gwaii: Is there a potential for recovery? Ecol. Manag 255, 3973–3979 (2008).
Google Scholar 
MacArthur, R. H. & Wilson, E. O. An equilibrium theory of insular zoogeography. Evolution 17, 373–387 (1963).
Google Scholar 
Catchen, J., Hohenlohe, P. A., Bassham, S., Amores, A. & Cresko, W. A. Stacks: an analysis tool set for population genomics. Mol. Ecol. 22, 3124–3140 (2013).PubMed 
PubMed Central 

Google Scholar 
Foll, M. & Gaggiotti, O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics 180, 977–993 (2008).PubMed 
PubMed Central 

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

Google Scholar 
Excoffier, L., Smouse, P. E. & Quattro, J. M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131, 479–491 (1992).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wilson, G. A. & Rannala, B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics 163, 1177–1191 (2003).PubMed 
PubMed Central 

Google Scholar 
Mussmann, S. M., Douglas, M. R., Chafin, T. K. & Douglas, M. E. BA3‐SNPs: contemporary migration reconfigured in BayesAss for next‐generation sequence data. Methods Ecol. Evol. 10, 1808–1813 (2019).
Google Scholar 
Loiselle, B. A., Sork, V. L., Nason, J. & Graham, C. Spatial genetic structure of a tropical understory shrub Psychotria officinalis (Rubiaceae). Am. J. Bot. 82, 1420 (1995).
Google Scholar 
Weir, B. S. & Cockerham, C. C. Estimating F-statistics for the analysis of population structure. Evolution 38, 1358–1370 (1984).CAS 
PubMed 
PubMed Central 

Google Scholar 
Goudet, J. Hierfstat, a package for r to compute and test hierarchical F-statistics. Mol. Ecol. Notes 5, 184–186 (2005).
Google Scholar 
Bellard, C., Cassey, P. & Blackburn, T. M. Alien species as a driver of recent extinctions. Biol. Lett. 12, 20150623 (2016).PubMed 
PubMed Central 

Google Scholar 
Robinette, W. L. Mule deer home range and dispersal in Utah. J. Wildl. Manag 30, 335 (1966).
Google Scholar 
Bunnell, F. L. & Harestad, A. S. Dispersal and dispersion of black-tailed deer: models and observations. J. Mammal. 64, 201–209 (1983).
Google Scholar 
Anderson, A. & Wallmo, O. Odocoileus hemionus. Am. Soc. Mammal. 219, 1–9 (1984).
Google Scholar 
Quigley, D. T. G. & Moffatt, S. Sika-like deer Cervus nippon Temminck, 1838 observed swimming out to sea at Greystones, Co. Wicklow: Increasing deer population pressure? Bull. Ir. Biogeogr. Soc. 38, 251–262 (2014).Colson, K. E., Brinkman, T. J., Person, D. K. & Hundertmark, K. J. Fine-scale social and spatial genetic structure in Sitka black-tailed deer. Conserv. Genet. 14, 439–449 (2013).
Google Scholar 
Powell, J. H. et al. Microsatellites indicate minimal barriers to mule deer Odocoileus hemionus dispersal across Montana, USA. Wildl. Biol. 19, 102–110 (2013).
Google Scholar 
Dolman, P. M. & Wäber, K. Ecosystem and competition impacts of introduced deer. Wildl. Res. 35, 202 (2008).
Google Scholar 
Edge, K.-A., Crouchley, D., McMurtrie, P., Willans, M. J. & Byrom, A. In Island Invasives: Eradication and Management (eds Veitch, C. R., Clout, M. N. & Towns, D. R.) 166–171 (IUCN, 2011).Hess, S. C., Muise, J. & Schipper, J. Anatomy of an eradication effort: removing Hawaii’s illegally introduced axis deer. Wildl. Prof. 9, 26–29 (2015).Masters, P., Markopoulos, N., Florance, B. & Southgate, R. The eradication of fallow deer (Dama dama) and feral goats (Capra hircus) from Kangaroo Island, South Australia. Australas. J. Environ. Manag. 25, 86–98 (2018).
Google Scholar 
Macdonald, N., Nugent, G., Edge, K.-A. & Parkes, J. P. In Island Invasives: Scaling up to Meet the Challenge Ch. 2 (eds Veitch, C. R., Clout, M. N., Russell, J. C. & West, C. J.) 256–260 (IUCN, 2019).Keitt, B. et al. In Island Invasives: Eradication and Management (eds Veitch, C. R., Clout, M. N. & Towns, D. R.) 74–77 (IUCN, 2011).Gerber, L. R. Conservation triage or injurious neglect in endangered species recovery. Proc. Natl Acad. Sci. USA 113, 3563–3566 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Baird, N. A. et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 3, e3376 (2008).PubMed 
PubMed Central 

Google Scholar 
Lemay, M. A. & Russello, M. A. Genetic evidence for ecological divergence in kokanee salmon. Mol. Ecol. 24, 798–811 (2015).CAS 
PubMed 

Google Scholar 
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
Rochette, N. C. & Catchen, J. M. Deriving genotypes from RAD-seq short-read data using Stacks. Nat. Protoc. 12, 2640–2659 (2017).CAS 
PubMed 

Google Scholar 
Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Beaumont, M. A. & Balding, D. J. Identifying adaptive genetic divergence among populations from genome scans. Mol. Ecol. 13, 969–980 (2004).CAS 
PubMed 

Google Scholar 
Lischer, H. E. L. & Excoffier, L. PGDSpider: an automated data conversion tool for connecting population genetics and genomics programs. Bioinformatics 28, 298–299 (2012).CAS 
PubMed 

Google Scholar 
Meirmans, P. G. & Tienderen, Van, P. H. Genotype and Genodive: two programs for the analysis of genetic diversity of asexual organisms. Mol. Ecol. Notes 4, 792–794 (2004).
Google Scholar 
Zheng, X. et al. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 28, 3326–3328 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol. Ecol. 14, 2611–2620 (2005).CAS 
PubMed 

Google Scholar 
Rosenberg, N. A. distruct: a program for the graphical display of population structure: PROGRAM NOTE. Mol. Ecol. Notes 4, 137–138 (2004).
Google Scholar 
Vekemans, X. & Hardy, O. J. New insights from fine-scale spatial genetic structure analyses in plant populations. Mol. Ecol. 13, 921–935 (2004).CAS 
PubMed 

Google Scholar 
Burgess, B. T., Irvine, R. L. & Russello, M. A. Population genomics of Sitka black-tailed deer supports invasive species management and ecological restoration on islands. Dryad https://doi.org/10.5061/dryad.q2bvq83mq (2022). More

Slade, C. & Law, B. The other half of the coastal State Forest estate in New South Wales; The value of informal forest reserves for conservation. Aust. Zool. 39, 359–370. https://doi.org/10.7882/AZ.2016.011 (2017).Article 

Google Scholar 
Munks, S. A., Chuter, A. E. & Koch, A. J. ‘Off-reserve’ management in practice: Contributing to conservation of biodiversity over 30 years of Tasmania’s forest practices system. For. Ecol. Manag. 465, 117941. https://doi.org/10.1016/j.foreco.2020.117941 (2020).Article 

Google Scholar 
Lande, R. Demographic models of the northern spotted owl (Strix occidentalis caurina). Oecologia 75, 601–607 (1988).ADS 
CAS 
Article 

Google Scholar 
Franklin, C. M. A., Macdonald, S. E. & Nielsen, S. E. Can retention harvests help conserve wildlife? Evidence for vertebrates in the boreal forest. Ecosphere 10(3), e02632 (2019).Article 

Google Scholar 
McAlpine, C. A. et al. Conserving koalas: A review of the contrasting regional trends, outlooks and policy challenges. Biol. Conserv. 192, 226–236. https://doi.org/10.1016/j.biocon.2015.09.020 (2015).Article 

Google Scholar 
Kavanagh, R. P. & Stanton, M. A. Koalas use young Eucalyptus plantations in an agricultural landscape on the Liverpool Plains, New South Wales. Ecol. Manag. Restor. 13, 297–305. https://doi.org/10.1111/emr.12005 (2012).Article 

Google Scholar 
Matthews, A., Lunney, D., Gresser, S. & Maitz, W. Movement patterns of koalas in remnant forest after fire. Aust. Mammal. 38, 91–104. https://doi.org/10.1071/AM14010 (2016).Article 

Google Scholar 
McAlpine, C. A. et al. The importance of forest area and configuration relative to local habitat factors for conserving forest mammals: A case study of koalas in Queensland, Australia. Biol. Conserv. 132, 153–165. https://doi.org/10.1016/j.biocon.2006.03.021 (2006).Article 

Google Scholar 
Beyer, H. L. et al. Management of multiple threats achieves meaningful koala conservation outcomes. J. Appl. Ecol. 55, 1966–1975. https://doi.org/10.1111/1365-2664.13127 (2018).Article 

Google Scholar 
Kavanagh, R. P., Stanton, M. A. & Brassil, T. E. Koalas continue to occupy their previous home-ranges after selective logging in Callitris–Eucalyptus forest. Wildl. Res. 34, 94–107. https://doi.org/10.1071/WR06126 (2007).Article 

Google Scholar 
Kavanagh, R. P., Debus, S., Tweedie, T. & Webster, R. Distribution of nocturnal forest birds and mammals in north-eastern New South Wales: Relationships with environmental variables and management history. Wildl. Res. 22, 359–377. https://doi.org/10.1071/WR9950359 (1995).Article 

Google Scholar 
Roberts, P. Associations Between Koala Faecal Pellets and Trees at Dorrigo, M.Sc. Thesis (University of New England, 1998).
Google Scholar 
Smith, A. P. Koala conservation and habitat requirements in a timber production forest in north-east New South Wales. In Conservation of Australia’s Forest Fauna (ed. Lunney, D.) 591–611 (Royal Zoological Society of New South Wales, 2004).Chapter 

Google Scholar 
Radford Miller, S. Aspects of the ecology of the koala, Phascolarctos cinereus, in a tall coastal production forest in north eastern New South Wales. PhD thesis (Southern Cross University, 2012).Law, B. S. et al. Passive acoustics and sound recognition provide new insights on status and resilience of an iconic endangered marsupial (koala Phascolarctos cinereus) to timber harvesting. PLoS One 13(10), e0205075. https://doi.org/10.1371/journal.pone.0205075 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Ellis, W. et al. Koala habitat use and population density: Using field data to test the assumptions of ecological models. Aust. Mammal. 35, 160–165. https://doi.org/10.1071/AM12023 (2013).Article 

Google Scholar 
Ashman, K. R., Rendall, A. R., Symonds, M. R. E. & Whisson, D. Understanding the role of plantations in the abundance of an arboreal folivore. Landsc. Urban Plan. 193, 103684. https://doi.org/10.1016/j.landurbplan.2019.103684 (2020).Article 

Google Scholar 
Cristescu, R. H., Rhodes, J., Frere, C. & Banks, P. B. Is restoring flora the same as restoring fauna? Lessons learned from koalas and mining rehabilitation. J. Appl. Ecol. 50(2), 423–431. https://doi.org/10.1111/1365-2664.12046 (2013).Article 

Google Scholar 
Chandler, R. B. & Royle, J. A. Spatially explicit models for inference about density in unmarked or partially marked populations. Ann. Appl. Stat. 7(2), 936–954. https://doi.org/10.1214/12-AOAS610 (2013).MathSciNet 
Article 
MATH 

Google Scholar 
Law, B., Gonsalves, L., Burgar, J., Brassil, T., Kerr, I., Wilmott, L., Madden, K., Smith, M., Mella, V., Crowther, M., Krockenberger, M., Rus, A., Pietsch, R., Truskinger, A., Eichinski, P. & Roe, P. Validation of spatial count models to estimate koala Phascolarctos cinereus density from acoustic arrays. Wildl. Res. (in press).MacKenzie, D. I. et al. Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence (Elsevier, 2006).MATH 

Google Scholar 
Smith, M. Behaviour of the Koala, Phascolarctos cinereus (Goldfuss), in Captivity III. Vocalisations. Wildl. Res. 7, 13–34. https://doi.org/10.1071/WR9800013 (1980).Article 

Google Scholar 
Ellis, W. et al. Koala bellows and their association with the spatial dynamics of free-ranging koalas. Behav. Ecol. 22, 372–377. https://doi.org/10.1093/beheco/arq216 (2011).Article 

Google Scholar 
Ellis, W. et al. The role of bioacoustic signals in koala sexual selection: Insights from seasonal patterns of associations revealed with gps-proximity units. PLoS One 10(7), e0130657. https://doi.org/10.1371/journal.pone.0130657 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Martin, R. W. Overbrowsing and decline of a population of the koala, Phascolarctos cinereus, in Victoria II. Population condition. Aust. Wildl. Res. 12, 367–375 (1985).ADS 
Article 

Google Scholar 
Penn, A. M. et al. Demographic forecasting in koala conservation. Conserv. Biol. 14(3), 629–638. https://doi.org/10.1046/j.1523-1739.2000.99385.x (2000).Article 

Google Scholar 
Watchorn, D. J. & Whisson, D. A. Quantifying the interactions between koalas in a high-density population during the breeding period. Aust. Mammal. 42(1), 28–37. https://doi.org/10.1071/AM18027 (2019).Article 

Google Scholar 
Crowther, M. S. et al. Comparison of three methods of estimating the population size of an arboreal mammal in a fragmented rural landscape. Wildl. Res. 48, 105–114. https://doi.org/10.1071/WR19148 (2020).Article 

Google Scholar 
Witt, R. R. et al. Real-time drone derived thermal imagery outperforms traditional survey methods for an arboreal forest mammal. PLoS One 15(11), e0242204. https://doi.org/10.1371/journal.pone.0242204 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Law, B.S, Gonsalves, L., Burgar, J., Brassil, T., Kerr I. & O’Loughlin C. Fire severity and its local extent are key to assessing impacts of Australian mega-fires on koala (Phascolarctos cinereus) density. Glob. Ecol. Biogeogr. 00, 1–13. https://doi.org/10.1111/geb.13458 (2022).Hynes, E. F., Whisson, D. A. & Di Stefano, J. Response of an arboreal species to plantation harvest. For. Ecol. Manag. 490, 119092. https://doi.org/10.1016/j.foreco.2021.119092 (2021).Article 

Google Scholar 
Law, B., Gonsalves, L., Burgar, J., Brassil, T., Kerr, I., O’Loughlin, C., Eichinski, P. & Roe, P. Regulated timber harvesting does not reduce koala density in north-east forests of New South Wales. Unpubl. Report to NSW (Natural Resources Commission, 2021).Phillips, S. Aversive behaviour by koalas (Phascolarctos cinereus) during the course of a music festival in northern New South Wales, Australia. Aust. Mammal. 38(2), 158–163. https://doi.org/10.1071/AM15006 (2016).Article 

Google Scholar 
Fedrowitz, K. et al. Can retention forestry help conserve biodiversity? A meta-analysis. J. Appl. Ecol. 51, 1669–1679. https://doi.org/10.1111/1365-2664.12289 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Mori, A. S. & Kitagawa, R. Retention forestry as a major paradigm for safeguarding forest biodiversity in productive landscapes: A global meta-analysis. Biol. Conserv. 175, 65–73. https://doi.org/10.1016/j.biocon.2014.04.016 (2014).Article 

Google Scholar 
Law, B. et al. Development and field validation of a regional, management-scale habitat model: A koala Phascolarctos cinereus case study. Ecol. Evol. 7, 7475–7489. https://doi.org/10.1002/ece3.3300 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Phillips, S., Wallis, K. & Lane, A. Quantifying the impacts of bushfire on populations of wild koalas (Phascolarctos cinereus): Insights from the 2019/20 fire season. Ecol. Manag. Restor. 22, 80–88. https://doi.org/10.1111/emr.12458 (2021).Article 

Google Scholar 
Kramer, A. et al. California spotted owl habitat selection in a fire-managed landscape suggests conservation benefit of restoring historical fire regimes. For. Ecol. Manag. 479, 118576 (2021).Article 

Google Scholar 
Jones, G. M. et al. Megafire causes persistent loss of an old-forest species. Anim. Conserv. 24, 925–936. https://doi.org/10.1111/acv.12697 (2021).Article 

Google Scholar 
Hagens, S. V., Rendall, A. R. & Whisson, D. A. Passive acoustic surveys for predicting species’ distributions: Optimising detection probability. PLoS One 13(7), e0199396. https://doi.org/10.1371/journal.pone.0199396 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Law, B. et al. Using passive acoustic recording and automated call identification to survey koalas in the southern forests of New South Wales. Aust. Zool. 40, 477–486 (2019).Article 

Google Scholar 
Towsey, M., Planitz, B., Nantes, A., Wimmer, J. & Roe, P. A toolbox for animal call recognition. Bioacoustics 21, 107–125. https://doi.org/10.1080/09524622.2011.648753 (2012).Article 

Google Scholar 
Royle, J. A. & Dorazio, R. M. Parameter-expanded data augmentation for Bayesian analysis of capture–recapture models. J. Ornithol. 152(2), 521–537 (2012).Article 

Google Scholar 
Royle, J. A., Chandler, R. B., Sollmann, R. & Gardner, B. Spatial Capture–Recapture 1st edn. (Elsevier, 2014). https://doi.org/10.1016/B978-0-12-405939-9.00020-7.Book 

Google Scholar 
Clark, J. D. Comparing clustered sampling designs for spatially explicit estimation of population density. Popul. Ecol. 61(1), 93–101. https://doi.org/10.1002/1438-390X.1011 (2019).Article 

Google Scholar 
Sun, C. C., Fuller, A. K. & Royle, J. A. Trap configuration and spacing influences parameter estimates in spatial capture-recapture models. PLoS One 9(2), e88025. https://doi.org/10.1371/journal.pone.0088025 (2014).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Plummer, M. JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling. In Proceedings of the 3rd International Workshop on Distributed Statistical Computing, Vol. 124(125.10), pp. 1–10 (2003).Plummer, M. rjags: Bayesian graphical models using MCMC. R package version 4(6) (2016).Burgar, J. M., Stewart, F. E., Volpe, J. P., Fisher, J. T. & Burton, A. C. Estimating density for species conservation: Comparing camera trap spatial count models to genetic spatial capture-recapture models. Glob. Ecol. Conserv. 15, e00411. https://doi.org/10.1016/j.gecco.2018.e00411 (2018).Article 

Google Scholar 
Stewart-Oaten, A., Murdoch, W. W. & Parker, K. R. Environmental impact assessment: “Pseudoreplication” in time?. Ecology 67(4), 929–940. https://doi.org/10.2307/1939815 (1986).Article 

Google Scholar 
Stewart-Oaten, A. & Bence, J. R. Temporal and spatial variation in environmental impact assessment. Ecol. Monogr. 71(2), 305–339. https://doi.org/10.1890/0012-9615(2001)071[0305:TASVIE]2.0.CO;2 (2001).Article 

Google Scholar  More

IPCC. Climate change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2013).
Google Scholar 
Pachauri, R. K. et al. Climate change 2014: Synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change (IPCC, 2014).
Google Scholar 
Reay, D. S. et al. Global agriculture and nitrous oxide emissions. Nat. Clim. Change 2, 410–416 (2012).ADS 
CAS 

Google Scholar 
Jassal, R. S., Black, T. A., Roy, R. & Ethier, G. Effect of nitrogen fertilization on soil CH4 and N2O fluxes, and soil and bole respiration. Geoderma 162, 182–186 (2011).ADS 
CAS 

Google Scholar 
Hu, H. W., Chen, D. & He, J. Z. Microbial regulation of terrestrial nitrous oxide formation: Understanding the biological pathways for prediction of emission rates. FEMS Microbiol. Rev. 39, 729–749 (2015).CAS 
PubMed 

Google Scholar 
Bateman, E. J. & Baggs, E. M. Contributions of nitrification and denitrification to N2O emissions from soils at different water-filled pore space. Biol. Fertil. Soils 41, 379–388 (2005).CAS 

Google Scholar 
Yang, Y. D., Hu, Y. G., Wang, Z. M. & Zeng, Z. H. Variations of the nirS-, nirK-, and nosZ-denitrifying bacterial communities in a northern Chinese soil as affected by different long-term irrigation regimes. Environ. Sci. Pollut. Res. 25, 14057–14067 (2018).CAS 

Google Scholar 
Pan, Y., Ye, L., Ni, B. & Yuan, Z. Effect of pH on N2O reduction and accumulation during denitrification by methanol utilizing denitrifiers. Water Res. 46, 4832–4840 (2012).CAS 
PubMed 

Google Scholar 
Hallin, S., Philippot, L., Loffler, F. E., Sanford, R. A. & Jones, C. M. Genomics and ecology of novel N2O-reducing microorganisms. Trends Microbiol. 1485, 43–55 (2017).
Google Scholar 
Yang, L., Zhang, X. & Ju, X. Linkage between N2O emission and functional gene abundance in an intensively managed calcareous fluvo-aquic soil. Sci. Rep. 7, 43283 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Cui, P. Y. et al. Long-term organic and inorganic fertilization alters temperature sensitivity of potential N2O emissions and associated microbes. Soil Biol. Biochem. 93, 131–141 (2016).CAS 

Google Scholar 
Gerber, J. S. et al. Spatially explicit estimates of N2O emissions from croplands suggest climate mitigation opportunities from improved fertilizer management. Glob. Change Biol. 22, 3383–3394 (2016).ADS 

Google Scholar 
Shcherbak, I., Millar, N. & Robertson, G. P. Global metaanalysis of the nonlinear response of soil nitrous oxide (N2O) emissions to fertilizer nitrogen. Proc. Natl. Acad. Sci. 111, 9199–9204 (2014).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wang, J., Chadwick, D. R., Cheng, Y. & Yan, X. Global analysis of agricultural soil denitrification in response to fertilizer nitrogen. Sci. Total Environ. 616, 908–917 (2018).ADS 
PubMed 

Google Scholar 
Albanito, F. et al. Direct nitrous oxide emissions from tropical and sub-tropical agricultural systems—A review and modelling of emission factors. Sci. Rep. 7, 44235 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Wolsing, M. & Priemé, A. Observation of high seasonal variation in community structure of denitrifying bacteria in arable soil receiving artificial fertilizer and cattle manure by determining T-RFLP of nir gene fragments. FEMS Microbiol. Ecol. 48, 261–271 (2004).CAS 
PubMed 

Google Scholar 
Akiyama, H., McTaggart, I. P., Ball, B. C. & Scott, A. N2O, NO, and NH3 emissions from soil after the application of organic fertilizers, urea and water. Water Air Soil Pollut. 156, 113–129 (2004).ADS 
CAS 

Google Scholar 
Wang, Y. Y. et al. Responses of N2O reductase gene (nosZ)-denitrifer communities to long-term fertilization follow a depth pattern in calcareous purplish paddy soil. J. Integr. Agric. 16, 2597–2611 (2017).CAS 

Google Scholar 
Fernandez-Luqueno, F. et al. Emission of CO2 and N2O from soil cultivated with common bean (Phaseolus vulgaris L.) fertilized with different N sources. Sci. Total Environ. 407, 4289–4296 (2009).ADS 
CAS 
PubMed 

Google Scholar 
Yin, C. et al. Different denitrification potential of aquic brown soil in Northeast China under inorganic and organic fertilization accompanied by distinct changes of nirS-and nirK-denitrifying bacterial community. Eur. J. Soil Biol. 65, 47–56 (2014).CAS 

Google Scholar 
Harter, J. et al. Linking N2O emissions from biochar-amended soil to the structure and function of the N-cycling microbial community. ISME J. 8, 660–674 (2013).PubMed 
PubMed Central 

Google Scholar 
Hai, B. et al. Quantification of key genes steering the microbial nitrogen cycle in the rhizosphere of sorghum cultivars in tropical agroecosystems. Appl. Environ. Microbiol. 75, 4993–5000 (2009).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Huang, R. et al. Nitrous oxide emission and the related denitrifier community: A short-term response to organic manure substituting chemical fertilizer. Ecotoxicol. Environ. Saf. 192, 110291 (2020).CAS 
PubMed 

Google Scholar 
Xu, X. et al. NosZ clade II rather than clade I determine in situ N2O emissions with different fertilizer types under simulated climate change and its legacy. Soil Biol. Biochem. 150, 107974 (2020).CAS 

Google Scholar 
Henderson, S. L. et al. Changes in denitrifier abundance, denitrification gene mRNA levels, nitrous oxide emissions, and denitrification in anoxic soil microcosms amended with glucose and plant residues. Appl. Environ. Microbiol. 76, 2155–2164 (2010).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Palmer, K., Biasi, C. & Horn, M. A. Contrasting denitrififier communities relate to contrasting N2O emission patterns from acidic peat soils in arctic tundra. ISME J. 6, 1058–1077 (2012).CAS 
PubMed 

Google Scholar 
Dandie, C. E. et al. Abundance, diversity and functional gene expression of denitrifier communities in adjacent riparian and agricultural zones. FEMS Microbiol. Ecol. 77, 69–82 (2011).CAS 
PubMed 

Google Scholar 
Avrahami, S., Conrad, R. & Braker, G. Effect of soil ammonium concentration on N2O release and on the community structure of ammonia oxidizers and denitrifiers. Appl. Environ. Microbiol. 68, 5685–5692 (2002).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Yang, Y. J. et al. Compost supplementation with nitrogen loss and greenhouse gas emissions during pig manure composting. Bioresour. Technol. 297, 122435 (2019).PubMed 

Google Scholar 
Yang, Y. J. et al. Exploring the microbial mechanisms of organic matter transformation during pig manure composting amended with bean dregs and biochar. Bioresour. Technol. 313, 123647 (2020).CAS 
PubMed 

Google Scholar 
Yang, J. H., Wang, C. L. & Dai, H. L. Agricultural Soil Analysis and Environmental Monitoring (China Land Press, 2008) (in Chinese).
Google Scholar 
Wang, Q. R., Li, Y. C. & Klassen, W. Changes of soil microbial biomass carbon and nitrogen with cover crops and irrigation in a tomato field. J. Plant Nutr. 30, 623–639 (2007).CAS 

Google Scholar 
Moore, J. M., Klose, S. & Tabatabai, M. A. Soil microbial biomass carbon and nitrogen as affected by cropping systems. Biol. Fertil. Soils 31, 200–210 (2000).CAS 

Google Scholar 
Jones, D. L. & Willett, V. B. Experimental evaluation of methods to quantify dissolved organic nitrogen (DON) and dissolved organic carbon (DOC) in soil. Soil Biol. Biochem. 38, 991–999 (2006).CAS 

Google Scholar 
Ghani, A., Dexter, M. & Perrott, K. W. Hot-water extractable carbon in soils: A sensitive measurement for determining impacts of fertilization, grazing and cultivation. Soil Biol. Biochem. 35, 1231–1243 (2003).CAS 

Google Scholar 
Huang, R. et al. Variation in N2O emission and N2O related microbial functional genes in straw- and biochar-amended and non-amended soils. Appl. Soil. Ecol. 137, 57–68 (2019).
Google Scholar 
Yang, Y. J. et al. Soil organic carbon transformation and dynamics of microorganisms under different organic amendments. Sci. Total Environ. 750, 141719 (2021).ADS 
CAS 
PubMed 

Google Scholar 
López-Fernández, S. et al. Effects of fertiliser type and the presence or absence of plants on nitrous oxide emissions from irrigated soils. Nutr. Cycl. Agroecosyst. 78, 279–289 (2007).
Google Scholar 
Wallenstein, M. D., Myrold, D. D., Firestone, M. & Voytek, M. Environmental controls on denitrifying communities and denitrification rates: insights from molecular methods. Ecol. Appl. 16, 2143–2152 (2006).PubMed 

Google Scholar 
Ciarlo, E., Conti, M., Bartoloni, N. & Rubio, G. Soil N2O emissions and N2O/(N2O+N2) ratio as affected by different fertilization practices and soil moisture. Biol. Fertil. Soils 44, 991–995 (2008).CAS 

Google Scholar 
Dandie, C. E. et al. Changes in bacterial denitrifier community abundance over time in an agricultural field and their relationship with denitrification activity. Appl. Environ. Microbiol. 74, 5997–6005 (2008).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Francis, C. A., O’Mullan, G. D., Cornwell, J. C. & Ward, B. B. Transitions in nirS-type denitrifier diversity, community composition, and biogeochemical activity along the Chesapeake Bay estuary. Front. Microbiol. 4, 237 (2013).PubMed 
PubMed Central 

Google Scholar 
Lin, W. et al. Evaluation of N2O sources after fertilizers application in vegetable soil by dual isotopocule plots approach. Environ. Res. 188, 109818 (2020).CAS 
PubMed 

Google Scholar 
Chen, M. M. et al. Nitrosospira cluster 3 lineage of AOB and nirK of Rhizobiales respectively dominated N2O emissions from nitrification and denitrification in organic and chemical N fertilizer treated soils. Ecol. Indic. 127, 107722 (2021).CAS 

Google Scholar 
Malghani, S., Kim, J., Lee, S. H., Yoo, G. Y. & Kang, H. Application of two contrasting rice-residue-based biochars triggered gaseous loss of nitrogen under denitrification-favoring conditions: a short-term study based on acetylene inhibition technique. Appl. Soil Ecol. 127, 112–119 (2018).
Google Scholar 
Sun, R., Guo, X., Wang, D. & Chu, H. Effects of long-term application of chemical and organic fertilizers on the abundance of microbial communities involved in the nitrogen cycle. Appl. Soil. Ecol. 95, 171–178 (2015).
Google Scholar 
Philippot, L., Andert, J., Jones, C. M., Bru, D. & Hallin, S. Importance of denitrifiers lacking the genes encoding the nitrous oxide reductase for N2O emissions from soil. Glob. Change Biol. 17, 1497–1504 (2011).ADS 

Google Scholar 
Chen, Z. et al. Differentiated response of denitrifying communities to fertilization regime in paddy soil. Microb. Ecol. 63, 446–459 (2012).PubMed 

Google Scholar 
Yoshida, M., Ishii, S., Otsuka, S. & Senoo, K. nirK-harboring denitrifiers are more responsive to denitrification-inducing conditions in rice paddy soil than nirS-harboring bacteria. Microbes Environ. 25, 45–48 (2010).PubMed 

Google Scholar 
Yin, C. et al. Denitrification potential under different fertilization regimes is closely coupled with changes in the denitrifying community in a black soil. Appl. Microbiol. Biotechnol. 99, 5719–5729 (2015).CAS 
PubMed 

Google Scholar 
Barrett, M. et al. Carbon amendment and soil depth affect the distribution and abundance of denitrifiers in agricultural soils. Environ. Sci. Pollut. Res. 23, 7899–7910 (2016).CAS 

Google Scholar 
Yoshida, M., Ishii, S., Otsuka, S. & Senoo, K. Temporal shifts in diversity and quantity of nirS and nirK in a rice paddy field soil. Soil Biol. Biochem. 41, 2044–2051 (2009).CAS 

Google Scholar 
Kandeler, E., Deiglmayr, K., Tscherko, D., Bru, D. & Philippot, L. Abundance of narG, nirS, nirK, and nosZ genes of denitrifying bacteria during primary successions of a glacier foreland. Appl. Environ. Microbiol. 72, 5957–5962 (2006).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar  More

General characteristics of study soilsTable 2 presents the descriptive statistics regarding the soil characteristics. Significant changes were observed in the distribution of sand (110–850 g kg−1), silt (50–530 g kg−1), clay (100–610 g kg−1), and soil textural class (7 texture classes) showing the diversity of natural and human processes involved in the formation and development of these soils28. Almost all soil samples were alkaline (with reaction at a range of 7.4–8.1) and calcareous (with CCE at a range of 5.5–35%). The EC of some soils was  > 4 dS/m (about 7% of the soil samples), indicating the partial salinity of the study soils. The organic carbon and total N contents of the soils were, on average, 2% (0.8–3.1%) and 0.28% (0.05–0.51%), respectively, placing them within the range of the moderate class. Likewise, the mean CEC of the soil, which is an effective indicator of soil fertility and quality, was in the moderate class of 12–25 cmol kg−129. The CEC was found to be highly correlated with clay (r = 0.76 P  Pb (58 mg kg−1)  > Ni (55.4 mg kg−1)  > Cu (38.8 mg kg−1)  > Cd (0.88 mg kg−1). In most soil samples, these ranges are comparable with data reported for other urban soils around the world—e.g. Ref.30 in Poland, Ref.31 in China, and Ref.32 in Greece. The values of Cd, Cu, and Zn were below their acceptable ranges as per the international standards4 in all soil samples. Nonetheless, the Pb and Ni contents were higher than their acceptable ranges in 13.1% and 17.4% of the samples, respectively. Furthermore, the concentrations of the five elements were higher than their background values in all urban soil samples. This difference was considerable for Cd, Pb, and Ni. The heavy metals had CV in the order of Cd (53%)  > Pb (51%)  > Ni (46%)  > Zn (21%)  > Cu (18%). This CV variation implies great variations in Cd, Pb, and Ni, which is linked to anthropogenic activities33. The background values of the metals, estimated by the median absolute deviation method10,14, were 52.3, 18.7, 0.45, 29.1, and 30.8 mg kg−1 for Zn, Cu, Cd, Pb, and Ni, respectively.We compared the concentrations of the heavy metals between urban and non-urban soils and found significant increases in the concentration of the metals in most soil types (Fig. 2). The urban soils had 17–36%, 14–21%, 41–70%, 43–69%, and 13–24% higher Zn, Cu, Cd, Pb, and Ni contents than the non-urban soils. The effluent and waste entry from multiple food processing and storage units, dying plants, metal plating facilities, and plastic production in close proximity of the study area is believed to be the reason for the high concentration of these trace elements. Research in various parts of the world, e.g., Ref.34 in India, Ref.35 in Brazil, and Ref.36 in China, has documented that the facilities have introduced significant quantities of heavy metals to soils. However, traffic and agrochemicals also play a key role in the accumulation of heavy metals in this region10.Figure 2The comparison of the mean values of Zn (a), Cu (b), Cd (c), Pb (d), and Ni (e) between urban and non-urban soils in different soil types. Different letters indicate significant differences in metal content within each soil type at P  Ni  > Cu. These findings are comparable to the results reported by37 and12. The highest EF for all five elements was observed in the Fluvisols soil type, reflecting that this soil type had been exposed to element pollution induced by urban activities to a greater extent than the other soil types. In a study on the pollution potential of four soil types in Central Greece, Ref.38 reported different ranges of element pollution across different soil types.Figure 3The comparison of the mean enrichment factor of Zn (a), Cu (b), Cd (c), Pb (d), and Ni (e) between urban and non-urban soils in different soil types. Different letters indicate significant differences in enrichment factor within each soil type at P  Pb (1.89)  > Ni (1.86)  > Cu (1.73)  > Zn (1.51). Mean PI for non-urban soils followed the order Cd (1.5)  > Zn (1.4)  > Cu (1.33)  > Pb (1.31)  > Ni (1.29). Nearly 7% and 16% of the urban soils showed moderate pollution (MP, PI = 2–3) and high pollution classes (HP, PI  > 3) of PI for Cd and 39% and 4% showed the MP and HP class of PI for Pb, respectively. However, the PI class was low pollution (PI = 1–2) for all soil samples and soil types in the non-urban soils. The results on the pollution index indicate a widespread intensification of soil pollution in urban soils across all studied heavy metals.Table 3 The level and terminology of PI and Ei of the analyzed heavy metals in urban and non-urban soils.Full size tableEcological risk, Ei was similarly found to be significantly higher in the urban soils than in the non-urban soils, even though the concentration of all elements except Cd fell within the low-risk class (Ei ≤ 40) in both urban and non-urban soils (Table 3). The mean Ei for Cd was 58.7 (moderate-risk class) and 39.2 (low-risk class) in the urban and non-urban soils, respectively. This means that urban activities have enhanced the ecological risk class of Cd by one grade. Overall, Cd had the highest EF, PI, and Ei among all heavy metals and in all soil samples, indicating a greater risk potential by Cd than Zn, Cu, Pb, and Ni across the water-soil–plant-human domain. Elevated Cd pollution by anthropogenic activities has been widely reported in the literature10,12,39. Cadmium as a Group 1 carcinogen element40 can accumulate in plant tissue without exhibiting visual symptoms. Therefore, Cd generally transfers from soil to the food chain covertly. Cadmium pollution can also influence soil quality and reduce crop yields and grain quality3.Similar to EF, PI, and Ei, the mean ER was significantly elevated in all urban soil types than the non-urban soils (Fig. 4). Among different soil types, the ER magnitude was in the order of Fluvisols (66.6%)  > Regosols (66.1%)  > Cambisols (59.8%)  > Calcisols (47%). These results indicate that Fluvisols carry a higher ecological risk potential for heavy metal accumulations than other soil types. In the study region, Fluvisols due to higher fertility and productivity are subject to more intense and extensive agronomic operations than other soil types13. Heavy application of agrochemicals (e.g., pesticides, herbicides, insecticides, and chemical fertilizers), accelerate the heavy metal input to the Fluvisols. Widespread application of nitrogen fertilizers and subsequent reduction in average soil pH markedly increases the solubility of certain heavy metals (e.g., Zn, Cu, Cd) which can be another factor increasing the ecological risk of heavy metal contamination in Fluvisols41. In addition, these Fluvisols are located on the margin of open urban wastewater channels, which are sometimes used for irrigation. A combination of mentioned processes can be implicated for higher ER of Fluvisols than that of other soil types as for BF, PI, and Ei.Figure 4The comparison of the mean ecological risk of selected heavy metals between urban and non-urban soils in different soil types. Different letters indicate significant differences in ecological risk within each soil type at P  Cu  > Ni  > Cd  > Pb in the roots, partially differing from that of the grain—Zn  > Cu  > Pb  > Ni  > Cd. Heavy metals concentrations observed in the corn roots and grains are almost comparable with those reported by42 in China and43 in Peru.Table 4 Summary statistical attributes of the concentration of heavy metals in corn root (R) and grain (G) along with their BCF and TF.Full size tableThe accumulation of heavy metals in the edible parts of corn is of higher importance. In the present study, the concentrations of these metals were lower than the acceptable level in the corn grains based on international references44. So, the consumption of corns grown in the regions should not threaten human and animal health in the short term, but caution should be exercised in their long-term consumption because some of these elements, especially Cd and Pb, which have long decomposition half-lives, gradually accumulate in body organs, especially in kidneys and livers45. Besides, the ratio of Zn, Cu, Cd, Pb, and Ni of the corn grain to their acceptable standard concentration, known as the pollution index of crop heavy metals, Ref.12 was lower than 0.7 for most corn samples, indicating the unpolluted risk class.The mean concentrations of Cd, Pb, and Ni were 5, 3.1, and 9.2 times as great in the corn roots as in their grains. This observation exhibits a notable phytoremediatory function of corn roots through restriction of radial translocation of heavy metals to the xylems and eventually into the grains. A similar trend of heavy metal accumulation in different plant organs has been reported in previous observations46,47. Based on Kabata-Pendias4 and Adriano22, plant cells can use the defensive tools of the roots to cope with heavy metals, especially Cd and Pb—highly toxic metals to plant cytosols. Accordingly, plant cells can fix these elements in the root system by such approaches as precipitating on cell walls, storing in vacuoles, and/or chelating by phytochelatins, thereby alleviating their toxic effects and inhibiting their translocation to plant shoots. For Zn, Cu, and Cd metals, a significant correlation was observed between their concentration in corn roots and grains. But, a less significant correlation (P  Cu (0.17)  > Zn (0.12)  > Ni (0.02)  > Pb (0.01). This implies that Cd, and to a smaller extent Cu is taken up by corn roots from the soil more readily, but Pb and Ni are less absorbable. These results are consistent with the reports of48 and46. The greater value of BCF-Cd may be related to a combination of the specific factors e,g., Cd concentration and chemistry, as well as soil characteristics (e.g., soil texture, pH, and calcium carbonate content)4. As was already discussed, the examined soils were characterized by high alkaline (pH = 7.4–8.1) and calcareous properties (CCE = 5.5–35%) with a high concentration of Soluble salts (EC = 0.7–6.6 dS m−1). These characteristics can result in the formation of complex Cd ions, especially CdOH+, CdCl20, CdCl+, CdSO40, and CdHCO3+4,22. These ions are plant-available, resulting in a further increase in Cd BCF. Regarding Ni and Pb, the alkaline and calcareous properties of the soils may have motivated insoluble compounds such as NiHCO3+ and NiCO30 (for Ni) and Pb(OH)2, PbCO3, PbSO4, and PbO (for Pb)4,22. These compounds cannot be uptake by plant roots, which may have resulted in a significant decrease in the BCF of these metals versus the other analyzed elements.Like BCF, the heavy metals had TF of  Pb (0.21)  > Cd (0.2)  > Ni (0.15). This implies that Zn and Cu are translocated from roots to grains readily, about four times as great as the other metals, while Ni, Cd, and Pb are translocated in smaller concentrations.The comparison of BCF and TF of Cd showed that less than 30% of Cd, on average, accumulated in the corn roots were translocated to the grains. This states that Cd is immobilized by various mechanisms before it can find its way into the grains. Some of the important mechanisms include (i) the antagonistic effects of Cd with other equivalent elements, especially Zn, Fe, and Ca, in the vascular system of corn, which reduces its mobility in the corn root-stem-grain system22, (ii) Cd sequestration in active exchange sites on the cell wall in the corn root-stem pathway10, and (iii) the binding of Cd with some specific compounds, e.g., phytochelatins of root vacuoles, which immobilizes it before its translocation to grains4,22. Lin and Aarts52 remarked that Cd mostly tends to be trapped in root vacuoles, which reduces its translocation to the upper parts of the plants. In general, it was found that corn plants have a high potential to absorb and accumulate Cd in their roots and Zn in their grains, which is consistent with previous studies41. For the majority of heavy metals, the values of BCF and TF in different soil types were in the order of Fluvisols  > Regosols  > Cambisols  > Calcisols, indicating that the great variety of soil types for the uptake and translocation of heavy metals in the soil-root-grain of the corn (Fig. 5).Figure 5Effect of soil type on the mean bioconcentration factor (a) and translocation factor (b) of selected heavy metals in urban soils. Different letters indicate significant differences in bioconcentration and translocation factors among soil types for each metal at P  Zn  > Cu  > Pb  > Ni for children, differing from that for adults (Cu  > Cd  > Pb  > Zn  > Ni). The values of HQ was  1 in over 87% of the samples, implying the low non-carcinogenic risk of this metal for corn-consuming children in the study region53. Rapidly developing children’s nervous system are highly sensitive to environmental factors, including heavy metals, so even a relatively low concentration of Cd in children’s blood may irreversibly affect their mental growth and functioning54.The highest HI was observed in children (min = 1.16, max = 2.31, mean = 1.63) followed by women and men which was similar to the found pattern of HQ (Table 7). These data show a moderate non-carcinogenic health risk (1 ≤ HI  More

Characteristics of heavy metal concentrationsOn the basis of the soil sample collection and chemical analysis, the concentration data for heavy metals in the urban soils of Wuxi were obtained. Through the statistical analysis of the soil heavy metal concentration data (Table 1), on the whole, the concentration of each heavy metal is as follows: Mn  > Zn  > Cr  > Ni  > Pb  > Cu  > Co  > Be  > Cd. Among these, the concentration range of Cr was 64.5–99 mg kg-1, and the average concentration was 72.9 mg kg−1. The concentration range of Ni was 31.4–67.5 mg kg−1, and the average concentration was 38.2 mg kg−1. The concentration range of Cu was 19.8–37.2 mg kg−1, and the average concentration was 25.5 mg kg−1. The concentration range of Zn was 72.4–1146 mg kg−1, and the average concentration was 90.2 mg kg−1. The concentration range of Cd was 0.34–1.06 mg kg−1, and the average concentration was 0.51 mg kg−1. The concentration range of Pb was 25.6–66.4 mg kg−1, and the average concentration was 37.6 mg kg−1. The variation coefficients of urban soil heavy metal concentration in Wuxi is between 0.09 and 0.33, which is less than 1. The spatial fluctuation of urban soil heavy metal concentration in Wuxi is small, indicating that the sources may be the same or similar.Table 1 Statistics of the heavy metal concentrations and Pb isotope ratios in the urban soils of Wuxi city (unit of heavy metal: mg kg−1; CV: coefficient of variation).Full size tableBy analysing the spatial distributions of the urban soil heavy metal concentrations in Wuxi, several obvious spatial distribution characteristics are found (Fig. 2). First, the heavy metals have high values in the central area of Wuxi, due to where has a high population density and various industries. The central aggregation of Pb is more obvious. Due to the dense roads in the city centre, vehicle traffic, bus stop signs and gas stations are mostly concentrated here, which will lead to Pb contents in this area that are significantly higher than those in other areas. In addition to the heavy metal concentrations, such as those for Cu, Zn and Cr in the downtown area, there are also areas with high values in western Wuxi and low values in eastern Wuxi. This phenomenon may be related to the land use types in Wuxi. In the western area of Wuxi, most land use types are urban and construction land, and the soils in this area are greatly disturbed by human activities. In the eastern region of Wuxi, woodland and grassland account for a large proportion of the land use types, which are less disturbed by human activities.Figure 2Spatial distribution characteristics of heavy metals in the urban soils of Wuxi city (unit: mg kg−1) [the figure was generated by Yan Li using the ArcGIS 10.2 (http:// https://developers.arcgis.com/)].Full size imageSource analysis of heavy metalsExploring for heavy metal pollution from emission sources is an important prerequisite for the study of urban soil pollution. By analysing the sources of heavy metals in soil environments, we can accurately determine which industries are major sources28,29,30 and whether there is homologous pollution. This is not only a theoretical basis for the study of lake sediment pollution and to clarify the risks brought by different pollution sources to the urban soil environment but also provides important guides for local government control of specific polluting industries and pollutant emissions. Based on this, the correlations and significance of heavy metals in the urban soils of Wuxi were analysed (Table 2). Generally, a heavy metal pollution source will emit multiple heavy metals at the same time. If the pollution source has a large emission, the concentration of these heavy metals in the environment will show a high level; on the contrary, if the emission of this pollution source is small, the concentration of these heavy metals in the environment will show a low level10. The correlations between the heavy metals Zn, Cr, Ni, Pb, Cu and Cd are between 0.655–0.907 and show strong correlations and significance at a level of 0.01. The strong significant correlations between different heavy metals indicate that these heavy metals have similar emission sources and transmission routes, which also means that they have consistent sources.Table 2 Correlations of Heavy Metals in the Urban Soils of Wuxi City.Full size tableTo further determine which industries are the sources of the heavy metals found in the urban soil of Wuxi, we analysed the Pb isotope data. The variation range of 208Pb/206Pb in soil is 2.09–2.12, and the average value is 2.10. The variation range of 206Pb/207Pb in soil is 1.17–1.18, and the average value is 1.177 (Table 1). After consulting relevant literature and materials, the main pollution sources of heavy metals in cities in eastern China include coal combustion, oil combustion, factory emissions, municipal wastes and so on3. Therefore, we collected the corresponding Pb isotope data in the emissions of heavy metal pollution sources. By collecting and comparatively analysing the Pb isotope data of known pollution sources (Fig. 3), it was determined that the Pb isotopes of the urban soil heavy metals in the soils of Wuxi have distinct characteristics. First, the Pb isotope distributions in the soils of Wuxi are relatively concentrated, and the ranges of variation are relatively small, which indicate that these heavy metals may have the same source or similar sets of sources. Second, the Pb isotopes in the urban soils of Wuxi city have few similarities with those of the uncontaminated soils and granites in eastern China; in contrast, the Pb isotopes in the urban soils of Wuxi are distributed in areas that are associated with coal combustion, automobile exhaust and urban waste (supplementary materials). The urban soil heavy metals in Wuxi generally have similar pollution sources and are greatly affected by human activities such as coal combustion and automobile exhaust emissions. Wuxi has a developed industrial economy and large numbers of factories. In the production and processing activities, the combustion of energy and fuel and the incomplete utilization of raw materials will lead to the enrichment of pollutants in the surrounding environment. By comparing other studies30,31, the Pb isotope analysis results in this study well indicate the source of soil heavy metals in Wuxi and make up for the Pb isotope data in this area. In the process of urban development, we should develop and apply clean energy, reduce the utilization of petroleum fossil fuels, and control the enrichment of heavy metals and other pollutants in the soil from the source.Figure 3Comparison of the Pb isotope compositions in the urban soils of Wuxi city with known sources.Full size imageEcological risk analysisBy calculating the potential ecological risk index for the heavy metals in the urban soils of Wuxi, the risks of heavy metals in the Wuxi soils were evaluated (Table 3). According to previous studies21, an Ei value lower than 40 indicates that a heavy metal is in a low-risk state at this location, and Ei values greater than or equal to 40 indicate that a heavy metal represents a high-risk state at this location. The average value of the potential ecological risk index of soil heavy metal Cd in Wuxi is 80.3, which represents a high-risk state. The average distributions of the potential ecological risk indexes of the heavy metals Cr, Cu, Zn, Pb and Ni are 1.8, 4.3, 1.1, 5.5 and 4.8, respectively, which all indicate a low-risk state. The risk statuses of different heavy metals may show certain correlations in space, which may be mutually complementary or antagonistic. Examining the spatial interactions of different heavy metal compound pollutants in urban soils plays an important role in the prevention and control of urban heavy metal pollution. Based on this, we used the Lisa analysis method to explore the spatial correlations of the different heavy metal risks in the urban soils of Wuxi (Fig. 4). The Moran scatter diagram can be divided into four quadrants that correspond to four different spatial patterns. High means that the variable value is higher than the average value, and Low means that the variable value is lower than the average value. In the upper right quadrant (High–High), a high-value area is surrounded by high-value neighbours; in the upper left quadrant (Low–High), a low-value area is surrounded by high-value neighbours; in the lower left quadrant (LL), a low-value area is surrounded by low-value neighbours; and in the lower right quadrant (High–Low), a high-value area is surrounded by low-value neighbours. High-High and Low-Low indicate that the differences between the region and its surrounding areas are small; that is, the regions with higher or lower values are concentrated, while the Low–High and High–Low quadrants indicate that the variable values between a region and its surrounding areas are different to a certain extent.Table 3 Ecological risk and health risk analysis of heavy metals in the urban soils of Wuxi (Cr-E represents the ecological risk of metal element Cr; Ni-E represents the ecological risk of metal element Ni; Cu-E represents the ecological risk of metal element Cu; Zn-E represents the ecological risk of metal element Zn; Cd-E represents the ecological risk of metal element Cd; Pb-E represents the ecological risk of metal element Pb; ADDderm-C is the average exposure to skin contact pathways for child; ADDderm-A is the average exposure to skin contact pathways for adult; ADDing-C is the average daily exposure to intake pathway for child; ADDing-A is the average daily exposure to intake pathway for adult; HI-C is the total health risk caused by accumulation of heavy metals in multiple ways in the same environmental medium for child; HI-A is the total health risk caused by accumulation of heavy metals in multiple ways in the same environmental medium for adult).Full size tableFigure 4LISA analysis of the ecological risks from different heavy metals [the figure was generated by Yan Li using the ArcGIS 10.2 (http:// https://developers.arcgis.com/)].Full size imageIn this study, two main results were obtained from spatial correlation Lisa analysis between different heavy metals. One is a High-High area, which is mainly distributed in the central and western regions of Wuxi city, which is consistent with the spatial distribution of the urban soil heavy metal concentrations in Wuxi city and is strongly disturbed by human activities. The other is the insignificant area, in which there are also large numbers of factories and enterprises and in which the forestland and grassland are distributed at intervals, which leads to an insignificant spatial correlation of soil heavy metal contents. Based on the above analysis, the high-risk areas for heavy metals in the urban soils of Wuxi are mainly concentrated in the central and western regions, and the relevant management activities need to be given great attention. In the eastern region, sporadic high-risk areas are also present, which should also receive due attention. Moran’s I is a method to measure the interdependence and degree of objects or phenomena by constructing statistics on certain characteristics or attributes for a certain spatial unit in the study area and the surrounding spatial units. It can be used to describe the spatial characteristics of spatial units such as aggregation or outliers in the distribution of certain attributes and is a very important technology in spatial data analysis33,34. However, few studies have applied it to the spatial relationship analysis of different heavy metals in urban soil.Health risk analysisBy using the health risk assessment model that is recommended by the U.S. EPA, this study calculated the health risks of soil heavy metals to adults and children through skin contact and ingestion. For both adults and children, the risk of soil heavy metals through ingestion was much higher than that caused by skin exposure (Table 3). For children, the total health risk that was caused by soil heavy metals is 0.078, which is four times that of adults. This may be related to children’s habits. Most children like to play with sand and climb around on the ground. These behaviours greatly increase the frequency of children contacting the soil, which thus increases the health risk caused by heavy metals in the soil. To further explore the spatial characteristics of the health risks of heavy metals in the soils of Wuxi, this study provides spatial predictions of the health risk values of soil heavy metals (Fig. 5). The total health risk values of soil heavy metals for children and adults have similar spatial distribution characteristics. High health risk values appear in the central area of Wuxi and decrease in a ring-shaped pattern. This is similar to the development degree of the city. The downtown area of Wuxi is densely populated, the pedestrian flow is very large, and the health risk of soil heavy metals in this area is very high, which poses a very serious potential threat. The health risk values for the western region of Wuxi are high, and there is also a potential threat. When compared with western Wuxi, eastern Wuxi has a lower risk.Figure 5Health risk analysis of heavy metals in the urban soils of Wuxi [the figure was generated by Yan Li using the ArcGIS 10.2 (http:// https://developers.arcgis.com/)].Full size image More

IPCC. in Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al) (Cambridge University Press. In Press, 2021).Otto, F. E. L. et al. Toward an inventory of the impacts of human-induced climate change. Bull. Am. Meteorol. Soc. 101, E1972–E1979 (2020).
Google Scholar 
Stanners, D. et al. in Sustainability Indicators. A Scientific Assessment (eds Moldan, B., Hak, T. & Dahl, A. L.) 127–144 (Island Press, 2007).Cohen-Shacham, E. et al. Core principles for successfully implementing and upscaling Nature-based Solutions. Environ. Sci. Policy 98, 20–29 (2019).
Google Scholar 
Seddon, N. et al. Getting the message right on nature-based solutions to climate change. Glob. Chang. Biol. https://doi.org/10.1111/gcb.15513 (2021).Keesstra, S. et al. The superior effect of nature based solutions in land management for enhancing ecosystem services. Sci. Total Environ. 610-611, 997–1009 (2018).CAS 

Google Scholar 
Seddon, N. et al. Understanding the value and limits of nature-based solutions to climate change and other global challenges. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190120 (2020).
Google Scholar 
Gómez Martín, E., Máñez Costa, M. & Schwerdtner Máñez, K. An operationalized classification of Nature Based Solutions for water-related hazards: from theory to practice. Ecol. Econ. 167 https://doi.org/10.1016/j.ecolecon.2019.106460 (2020).Doswald, N. et al. Effectiveness of ecosystem-based approaches for adaptation: review of the evidence-base. Clim. Dev. 6, 185–201 (2014).
Google Scholar 
Chausson, A. et al. Mapping the effectiveness of nature-based solutions for climate change adaptation. Glob. Chang. Biol. https://doi.org/10.1111/gcb.15310 (2020).Rebelo, A. J., Holden, P. B., Esler, K. & New, M. G. Benefits of water-related ecological infrastructure investments to support sustainable land-use: a review of evidence from critically water-stressed catchments in South Africa. R. Soc. Open Sci. 8, 201402 (2021).
Google Scholar 
Berrang-Ford, L. et al. A systematic global stocktake of evidence on human adaptation to climate change. Nat. Clim. Change 11, 989–1000 (2021).
Google Scholar 
Griscom, B. W. et al. Natural climate solutions. Proc. Natl Acad. Sci. USA 114, 11645–11650 (2017).CAS 

Google Scholar 
Bastin, J.-F. et al. The global tree restoration potential. Science 365, 76–79 (2019).CAS 

Google Scholar 
Koch, A., Brierley, C. & Lewis, S. L. Effects of Earth system feedbacks on the potential mitigation of large-scale tropical forest restoration. Biogeosciences 18, 2627–2647 (2021).CAS 

Google Scholar 
Girardin, C. A. J. et al. Nature-based solutions can help cool the planet – if we act now. Nature 593, 191–194 (2021).CAS 

Google Scholar 
Sudmeier-Rieux, K. et al. Scientific evidence for ecosystem-based disaster risk reduction. Nat. Sustain. 4, 803–810 (2021).
Google Scholar 
Otto, F. E. L. Attribution of weather and climate events. Annu. Rev. Environ. Resources 42, 627–646 (2017).
Google Scholar 
Philip, S. et al. A protocol for probabilistic extreme event attribution analyses. Adv. Stat. Climatol. Meteorol. Oceanogr. 6, 177–203 (2020).
Google Scholar 
Herring, S. C., Christidis, N., Hoell, A., Hoerling, M. P. & Stott, P. A. Explaining extreme events of 2019 from a climate perspective. Bull. Amer. Meteorol. Soc. 102, S1–S112 (2021).Otto, F. E. L. et al. Challenges to understanding extreme weather changes in lower income countries. Bull. Am. Meteorol. Soc. https://doi.org/10.1175/bams-d-19-0317.1 (2020).Pall, P. et al. Anthropogenic greenhouse gas contribution to flood risk in England and Wales in autumn 2000. Nature 470, 382–385 (2011).CAS 

Google Scholar 
Kay, A. L., Crooks, S. M., Pall, P. & Stone, D. A. Attribution of Autumn/Winter 2000 flood risk in England to anthropogenic climate change: a catchment-based study. J. Hydrol. 406, 97–112 (2011).
Google Scholar 
Schaller, N. et al. Human influence on climate in the 2014 southern England winter floods and their impacts. Nat. Clim. Change 6, 627–634 (2016).
Google Scholar 
Wolski, P., Stone, D., Tadross, M., Wehner, M. & Hewitson, B. Attribution of floods in the Okavango basin, Southern Africa. J. Hydrol. 511, 350–358 (2014).
Google Scholar 
Ross, A. C. et al. Anthropogenic influences on extreme annual streamflow into Chesapeake Bay from the Susquehanna River. Bull. Am. Meteorol. Soc. 102, S25–S32 (2021).Mitchell, D. et al. Attributing human mortality during extreme heat waves to anthropogenic climate change. Environ. Res. Lett. 11, 074006 (2016).
Google Scholar 
Botai, C., Botai, J., de Wit, J., Ncongwane, K. & Adeola, A. Drought Characteristics over the Western Cape Province, South Africa. Water https://doi.org/10.3390/w9110876 (2017).Wolski, P. How severe is Cape Town’s “Day Zero” drought? Significance 15, 24–27 (2018).
Google Scholar 
Stafford, L., Shemie, D., Kroeger, T., Baker, T. & Apse, C. The Greater Cape Town Water Fund. Assessing the return on investment for Ecological Infrastructure restoration. Business case. (The Nature Conservancy, 2018).Otto, F. E. L. et al. Anthropogenic influence on the drivers of the Western Cape drought 2015–2017. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/aae9f9 (2018).Pascale, S., Kapnick, S. B., Delworth, T. L. & Cooke, W. F. Increasing risk of another Cape Town “Day Zero” drought in the 21st century. Proc. Natl. Acad. Sci. USA https://doi.org/10.1073/pnas.2009144117 (2020).Van Wilgen, B. W., Measey, J., Richardson, D. M., Wilson, J. R. & Zengeya, T. A. Biological Invasions in South Africa (Springer Nature, 2020).Le Maitre, D. et al. Impacts of plant invasions on terrestrial water flows in South Africa in Biological Invasions in South Africa (eds van Wilgen, B. W., Measey. J., Richardson, D. M., Wilson, J. R. & Zengeya, T. A.) 431–457 (Springer, 2020).Brown, A. E., Zhang, L., McMahon, T. A., Western, A. W. & Vertessy, R. A. A review of paired catchment studies for determining changes in water yield resulting from alterations in vegetation. J. Hydrol. 310, 28–61 (2005).
Google Scholar 
Dennedy-Frank, P. J. & Gorelick, S. M. Insights from watershed simulations around the world: watershed service-based restoration does not significantly enhance streamflow. Glob. Environ. Change https://doi.org/10.1016/j.gloenvcha.2019.101938 (2019).Calder, I. D. & Dye, P. Hydrological impacts of invasive alien plants. Land Use Water Resour. Res. 7, 1–12 (2001).
Google Scholar 
Trabucco, A., Zomer, R. J., Bossio, D. A., van Straaten, O. & Verchot, L. V. Climate change mitigation through afforestation/reforestation: a global analysis of hydrologic impacts with four case studies. Agric. Ecosyst. Environ. 126, 81–97 (2008).
Google Scholar 
Farley, K. A., Jobbagy, E. G. & Jackson, R. B. Effects of afforestation on water yield: a global synthesis with implications for policy. Glob. Change Biol. 11, 1565–1576 (2005).
Google Scholar 
Jackson, R. B. Trading water for carbon with biological carbon sequestration. Science 310, 1944–1947 (2005).Filoso, S., Bezerra, M. O., Weiss, K. C. B. & Palmer, M. A. Impacts of forest restoration on water yield: a systematic review. PLoS ONE 12, e0183210 (2017).
Google Scholar 
Sitzia, T., Campagnaro, T., Kowarik, I. & Trentanovi, G. Using forest management to control invasive alien species: helping implement the new European regulation on invasive alien species. Biol. Invasions 18, 1–7 (2015).
Google Scholar 
Richardson, D. M. & Rejmánek, M. Trees and shrubs as invasive alien species – a global review. Divers. Distrib. 17, 788–809 (2011).
Google Scholar 
Everard, M. et al. Can control of invasive vegetation improve water and rural livelihood security in Nepal? Ecosyst. Serv. 32, 125–133 (2018).
Google Scholar 
Everard, M. Can management of ‘thirsty’ alien trees improve water security in semi-arid India? Sci. Total Environ. 704, 135451 (2020).CAS 

Google Scholar 
Archer, S. R. et al. Woody plant encroachment: causes and consequences in Rangeland Systems Springer Series on Environmental Management (ed. Briske, D. D.) Chapter 2, 25–84 (2017).Wood, M. Bootstrapped confidence intervals as an approach to statistical inference. Organ. Res. Methods 8, 454–470 (2016).
Google Scholar 
Tan, S. H. The correct interpretation of confidence intervals. Proc. Singapore Healthc. 19 (2010).Coetsee, C., Gray, E. F., Wakeling, J., Wigley, B. J. & Bond, W. J. Low gains in ecosystem carbon with woody plant encroachment in a South African savanna. J. Trop. Ecol. 29, 49–60 (2012).
Google Scholar 
Stevens, N., Erasmus, B. F., Archibald, S. & Bond, W. J. Woody encroachment over 70 years in South African savannahs: overgrazing, global change or extinction aftershock? Philos. Trans. R. Soc. Lond. B Biol. Sci. https://doi.org/10.1098/rstb.2015.0437 (2016).Venter, Z. S., Cramer, M. D. & Hawkins, H. J. Drivers of woody plant encroachment over Africa. Nat. Commun. 9, 2272 (2018).CAS 

Google Scholar 
Forsyth, G. G., Le Maitre, D. C., Smith, J. & Lotter, D. Upper Berg River Catchment (G10A) Management Unit Control Plan. (Natural Resources Management (NRM) Department of Environmental Affairs, 2016).Dirmeyer, P. A., Balsamo, G., Blyth, E. M., Morrison, R. & Cooper, H. M. Land‐atmosphere interactions exacerbated the drought and heatwave over northern Europe during summer 2018. AGU Adv. 2, e2020AV000283 (2021).
Google Scholar 
Rejmánek, M., Richardson, D. M. & Pysek, P. Trees and shrubs as invasive alien species – 2013 update of the global database. Divers. Distrib. 19, 1093–1094 (2013).
Google Scholar 
Terrer, C. et al. A trade-off between plant and soil carbon storage under elevated CO2. Nature 591, 599–603 (2021).CAS 

Google Scholar 
Ziervogel, G. et al. Climate change impacts and adaptation in South Africa. Wiley Interdiscip. Rev. Clim. Change 5, 605–620 (2014).
Google Scholar 
Thomas, A. et al. Global evidence of constraints and limits to human adaptation. Reg. Environ. Change https://doi.org/10.1007/s10113-021-01808-9 (2021).Dow, K., Berkhout, F. & Preston, B. L. Limits to adaptation to climate change: a risk approach. Curr. Opin. Environ. Sustain. 5, 384–391 (2013).
Google Scholar 
Manning, J. & Goldblatt, P. Plants of the greater Cape Floristic Region 1: the Core Cape Flora., (South African National Biodiversity Institute, 2012).Nel, J. L. et al. Strategic water source areas for urban water security: Making the connection between protecting ecosystems and benefiting from their services. Ecosyst. Serv. 28, 251–259 (2017).
Google Scholar 
Wolski, P. What Cape Town learned from its drought. Bull. At. Sci. https://thebulletin.org/2018/04/what-cape-town-learned-from-its-drought/ (2018).D. W. S. Cape Town River Systems State of Dams on 2021-08-16. Department of Water and Sanitation. Republic of South Africa. https://www.dws.gov.za/Hydrology/Weekly/RiverSystems.aspx?river=CT (2021).Rebelo, A. J. et al. The hydrological benefits of restoration: a modelling study of alien tree clearing in four mountain catchments in South Africa. Preprint at J. Hydrol. https://doi.org/10.21203/rs.3.rs-1316834/v1.DWAF. The Assessment of Water Availability in the Berg Catchment (WMA 19) by Means of Water Resource Related Models: Report 9 (Groundwater Model): Volume 9 (Breede River Alluvium Aquifer Model). (Department of Water Affairs and Forestry, 2008).DWAF. The Assessment of Water Availability in the Berg Catchment (WMA 19) by Means of Water Resource Related Models: Report 9 (Groundwater Model): Volume 3 (Regional Conceptual Model). (Department of Water Affairs and Forestry, 2008).Blake, D., Mlisa, A. & Hartnady, C. Large scale quantification of aquifer storage and volumes from the Peninsula and Skurweberg Formations in the southwestern Cape. Water SA 36, 177–184 (2010).
Google Scholar 
Holden, P. B., Rebelo, A. J. & New, M. G. Mapping invasive alien trees in water towers: a combined approach using satellite data fusion, drone technology and expert engagement. Remote Sens. Appl.: Soc. Environ. https://doi.org/10.1016/j.rsase.2020.100448 (2021).Midgley, J. & Scott, D. The use of stable isotopes of water in hydrological studies in the Jonkershoek Valley. Water SA 20, 151–154 (1994).
Google Scholar 
Van Genuchten, M. T. A closed‐form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892–898 (1980).
Google Scholar 
Moriasi, D. N., Gitau, M. W., Pai, N. & Daggupati, P. Hydrologic and water quality models: performance measures and evaluation criteria. Trans. ASABE 58, 1763–1785 (2015).
Google Scholar 
Stone, D. A. et al. A basis set for exploration of sensitivity to prescribed ocean conditions for estimating human contributions to extreme weather in CAM5.1-1degree. Weather Clim. Extremes 19, 10–19 (2018).
Google Scholar 
Risser, M. D., Stone, D. A., Paciorek, C. J., Wehner, M. F. & Angélil, O. Quantifying the effect of interannual ocean variability on the attribution of extreme climate events to human influence. Clim. Dyn. 49, 3051–3073 (2017).
Google Scholar 
Jones, G. S., Stott, P. A. & Christidis, N. Attribution of observed historical near-surface temperature variations to anthropogenic and natural causes using CMIP5 simulations. J. Geophys. Res. Atmos. 118, 4001–4024 (2013).
Google Scholar 
Sun, L. et al. Drivers of 2016 record Arctic warmth assessed using climate simulations subjected to factual and counterfactual forcing. Weather Clim. Extremes 19, 1–9 (2018).CAS 

Google Scholar 
Guillod, B. P. et al. weather@home 2: validation of an improved global–regional climate modelling system. Geosci. Model Dev. 10, 1849–1872 (2017).
Google Scholar 
Massey, N. et al. weather@home—development and validation of a very large ensemble modelling system for probabilistic event attribution. Q. J. R. Meteorol. Soc. 141, 1528–1545 (2014).
Google Scholar 
Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).
Google Scholar 
Flato, G. et al. in Climate change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change 741–866 (Cambridge University Press, 2014).Hargreaves, G. H. & Samani, Z. A. Reference crop evapotranspiration from temperature. Appl. Eng. Agriculture 1, 96–99 (1985).
Google Scholar 
Cayan, D. R., Maurer, E. P., Dettinger, M. D., Tyree, M. & Hayhoe, K. Climate change scenarios for the California region. Clim. Change 87, 21–42 (2008).
Google Scholar 
Cannon, A. J., Sobie, S. R. & Murdock, T. Q. Bias correction of GCM precipitation by quantile mapping: how well do methods preserve changes in quantiles and extremes? J. Clim. 28, 6938–6959 (2015).
Google Scholar 
R Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing. https://www.R-project.org/, 2020).Paciorek, C. J., Stone, D. A. & Wehner, M. F. Quantifying statistical uncertainty in the attribution of human influence on severe weather. Weather Clim. Extremes 20, 69–80 (2018).
Google Scholar 
Tadono, T. et al. Generation of the 30 M-Mesh Global Digital Surface Model by Alos Prism. ISPRS – International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences XLI-B4, 157–162, https://developers.google.com/earth-engine/datasets/catalog/JAXA_ALOS_AW3D30_V3_2#description (2016).
Google Scholar 
Takaku, J., Tadono, T., Tsutsui, K. & Ichikawa, M. Validation of “Aw3d” Global Dsm Generated from Alos Prism. ISPRS Ann. Photogramm. III-4, 25–31 (2016).
Google Scholar 
Viviroli, D. Increasing dependence of lowland population on mountain water resources. Nat. Sustain. 3, 917–928 (2020).
Google Scholar 
Meybeck, M. A New typology for mountains and other relief classes: an application to global continental water resources and population distribution. Mt. Res. Dev. 21, 34–45 (2001).DWS. Surface water home. Department of Water and Sanitation. Republic of South Africa. https://www.dws.gov.za/Hydrology/Unverified/UnverifiedDataFlowInfo.aspx (2021). More

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain
the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in
Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles
and JavaScript. More