Ecology

Tusting, L. S. et al. Mapping changes in housing in sub-Saharan Africa from 2000 to 2015. Nature 568, 391–394 (2019).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lozano, R. et al. Measuring progress from 1990 to 2017 and projecting attainment to 2030 of the health-related Sustainable Development Goals for 195 countries and territories: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 392, 2091–2138 (2018).
Google Scholar 
Kassebaum, N. J. et al. Global, regional, and national levels and causes of maternal mortality during 1990–2013: A systematic analysis for the Global Burden of Disease Study 2013. Lancet 384, 980–1004 (2014).PubMed 
PubMed Central 

Google Scholar 
Dhiman, S. Are malaria elimination efforts on right track? An analysis of gains achieved and challenges ahead. Infect. Dis. Poverty 8, 1–19 (2019).
Google Scholar 
WHO. World Malaria Report 2018 (WHO, 2018).
Google Scholar 
Janko, M. M. et al. The links between agriculture, Anopheles mosquitoes, and malaria risk in children younger than 5 years in the Democratic Republic of the Congo: A population-based, cross-sectional, spatial study. Lancet Planet. Health 2, e74–e82 (2018).PubMed 
PubMed Central 

Google Scholar 
Gerland, P. et al. World population stabilization unlikely this century. Science 346, 234–237 (2014).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Jayne, T. S., Chamberlin, J. & Headey, D. D. Land pressures, the evolution of farming systems, and development strategies in Africa: A synthesis. Food Policy 48, 1–17 (2014).
Google Scholar 
Ordway, E. M., Asner, G. P. & Lambin, E. F. Deforestation risk due to commodity crop expansion in sub-Saharan Africa. Environ. Res. Lett. 12, 044015 (2017).ADS 

Google Scholar 
Chaves, L. S. M. et al. Global consumption and international trade in deforestation-associated commodities could influence malaria risk. Nat. Commun. 11, 1–10 (2020).
Google Scholar 
Adenle, A. A., Azadi, H. & Manning, L. The era of sustainable agricultural development in Africa: Understanding the benefits and constraints. Food Rev. Int. 34, 411–433 (2018).
Google Scholar 
Ijumba, J. N. & Lindsay, S. W. Impact of irrigation on malaria in Africa: Paddies paradox. Med. Vet. Entomol. 15, 1–11 (2001).CAS 
PubMed 

Google Scholar 
Warra, A. A. & Prasad, M. N. V. African perspective of chemical usage in agriculture and horticulture—their impact on human health and environment. In Agrochemicals, Detection Treatment and Remediation 401–436 (Elsevier, 2020).
Google Scholar 
Fornace, K. M., Diaz, A. V., Lines, J. & Drakeley, C. J. Achieving global malaria eradication in changing landscapes. Malar. J. 20, 1–14 (2021).
Google Scholar 
Burkett-Cadena, N. D. & Vittor, A. Y. Deforestation and vector-borne disease: Forest conversion favors important mosquito vectors of human pathogens. Basic Appl. Ecol. 26, 101–110 (2018).PubMed 

Google Scholar 
Lindblade, K. A., Walker, E. D., Onapa, A. W., Katungu, J. & Wilson, M. L. Land use change alters malaria transmission parameters by modifying temperature in a highland area of Uganda. Trop. Med. Int. Heal. 5, 263–274 (2000).CAS 

Google Scholar 
Yasuoka, J. & Levins, R. Impact of deforestation and agricultural development on anopheline ecology and malaria epidemiology. Am. J. Trop. Med. Hyg. 76, 450–460 (2007).PubMed 

Google Scholar 
Guerra, C. A., Snow, R. W. & Hay, S. I. A global assessment of closed forests, deforestation and malaria risk. Ann. Trop. Med. Parasitol. 100, 189–204 (2006).CAS 
PubMed 

Google Scholar 
Laporta, G. Z., de Prado, P. I. K. L., Kraenkel, R. A., Coutinho, R. M. & Sallum, M. A. M. Biodiversity can help prevent malaria outbreaks in tropical forests. PLoS Negl. Trop. Dis. 7, e2139 (2013).PubMed 
PubMed Central 

Google Scholar 
Baeza, A., Santos-Vega, M., Dobson, A. P. & Pascual, M. The rise and fall of malaria under land-use change in frontier regions. Nat. Ecol. Evol. 1, 0108 (2017).
Google Scholar 
Patz, J. A., Graczyk, T. K., Geller, N. & Vittor, A. Y. Effects of environmental change on emerging parasitic diseases. Int. J. Parasitol. https://doi.org/10.1016/S0020-7519(00)00141-7 (2000).Article 
PubMed 

Google Scholar 
Sogoba, N. et al. Spatial analysis of malaria transmission parameters in the rice cultivation area of Office du Niger, Mali. Am. J. Trop. Med. Hyg. 76, 1009–1015 (2007).PubMed 

Google Scholar 
Mwangangi, J. M. et al. Anopheles larval abundance and diversity in three rice agro-village complexes Mwea irrigation scheme, central Kenya. Malar. J. 9, 1–10 (2010).
Google Scholar 
Diuk-Wasser, M. A. et al. Patterns of irrigated rice growth and malaria vector breeding in Mali using multi-temporal ERS-2 synthetic aperture radar. Int. J. Remote Sens. 27, 535–548 (2006).CAS 
PubMed 
PubMed Central 

Google Scholar 
Briët, O. J. T., Dossou-Yovo, J., Akodo, E., Van De Giesen, N. & Teuscher, T. M. The relationship between Anopheles gambiae density and rice cultivation in the savannah zone and forest zone of Côte d’Ivoire. Trop. Med. Int. Heal. 8, 439–448 (2003).
Google Scholar 
Klinkenberg, E., McCall, P. J., Wilson, M. D., Amerasinghe, F. P. & Donnelly, M. J. Impact of urban agriculture on malaria vectors in Accra, Ghana. Malar. J. 7, 1–9 (2008).
Google Scholar 
Keiser, J. et al. Effect of irrigation and large dams on the burden of malaria on a global and regional scale. Am. J. Trop. Med. Hyg. 72, 392–406 (2005).PubMed 

Google Scholar 
Kyei-Baafour, E. et al. Impact of an irrigation dam on the transmission and diversity of Plasmodium falciparum in a seasonal malaria transmission area of Northern Ghana. J. Trop. Med. 2020, 1–8 (2020).
Google Scholar 
Kibret, S. Time to revisit how dams are affecting malaria transmission. Lancet Planet. Heal. 2, e378–e379 (2018).
Google Scholar 
Kibret, S., Lautze, J., McCartney, M., Nhamo, L. & Yan, G. Malaria around large dams in Africa: Effect of environmental and transmission endemicity factors. Malar. J. 18, 1–12 (2019).
Google Scholar 
Kibret, S., Wilson, G. G., Ryder, D., Tekie, H. & Petros, B. Malaria impact of large dams at different eco-epidemiological settings in Ethiopia. Trop. Med. Health 45, 1–14 (2017).
Google Scholar 
Keiser, J., Singer, B. H. & Utzinger, J. Reducing the burden of malaria in different eco-epidemiological settings with environmental management: A systematic review. Lancet Infect. Dis. 5, 695–708 (2005).PubMed 

Google Scholar 
Ijumba, J. N., Shenton, F. C., Clarke, S. E., Mosha, F. W. & Lindsay, S. W. Irrigated crop production is associated with less malaria than traditional agricultural practices in Tanzania. Trans. R. Soc. Trop. Med. Hyg. 96, 476–480 (2002).CAS 
PubMed 

Google Scholar 
Ijumba, J. N., Mosha, F. W. & Lindsay, S. W. Malaria transmission risk variations derived from different agricultural practices in an irrigated area of northern Tanzania. Med. Vet. Entomol. 16, 28–38 (2002).CAS 
PubMed 

Google Scholar 
Klinkenberg, E. et al. Malaria and irrigated crops, Accra, Ghana. Emerg. Infect. Dis. 11, 1290–1293 (2005).PubMed 
PubMed Central 

Google Scholar 
USAID. The DHS Program—DHS Methodology. https://www.dhsprogram.com/What-We-Do/Survey-Types/DHS-Methodology.cfm (1984).Siraj, A. S. et al. Altitudinal changes in malaria incidence in highlands of Ethiopia and Colombia. Science 343, 1154–1158 (2014).ADS 
CAS 
PubMed 

Google Scholar 
Pascual, M., Ahumada, J. A., Chaves, L. F., Rodó, X. & Bouma, M. Malaria resurgence in the East African highlands: Temperature trends revisited. Proc. Natl. Acad. Sci. U. S. A. 103, 5829–5834 (2006).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lindblade, K. A. et al. Sustainability of reductions in malaria transmission and infant mortality in Western Kenya with use of insecticide-treated bednets 4 to 6 years of follow-up. J. Am. Med. Assoc. 291, 2571–2580 (2004).CAS 

Google Scholar 
Mordecai, E. A. et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol. Lett. 16, 22–30 (2013).PubMed 

Google Scholar 
Weiss, D. J. et al. Re-examining environmental correlates of Plasmodium falciparum Malaria endemicity: A data-intensive variable selection approach. Malar. J. 14, 1–18 (2015).
Google Scholar 
Bauhoff, S. & Busch, J. Does deforestation increase malaria prevalence? Evidence from satellite data and health surveys. World Dev. 127, 104734 (2020).
Google Scholar 
Austin, K. F., Bellinger, M. O. & Rana, P. Anthropogenic forest loss and malaria prevalence: A comparative examination of the causes and disease consequences of deforestation in developing nations. AIMS Environ. Sci. 4, 217–231 (2017).
Google Scholar 
Kabaria, C. W., Gilbert, M., Noor, A. M., Snow, R. W. & Linard, C. The impact of urbanization and population density on childhood Plasmodium falciparum parasite prevalence rates in Africa. Malar. J. 16, 1–10 (2017).
Google Scholar 
Herrera, D. et al. Upstream watershed condition predicts rural children’s health across 35 developing countries. Nat. Commun. 8, 1–8 (2017).CAS 

Google Scholar 
Van Ittersum, M. K. et al. Can sub-Saharan Africa feed itself?. Proc. Natl. Acad. Sci. U. S. A. 113, 14964–14969 (2016).PubMed 
PubMed Central 

Google Scholar 
Ickowitz, A. Shifting cultivation and deforestation in tropical Africa: Critical reflections. Dev. Change 37, 599–626 (2006).
Google Scholar 
Kar, N. P., Kumar, A., Singh, O. P., Carlton, J. M. & Nanda, N. A review of malaria transmission dynamics in forest ecosystems. Parasit. Vectors 7, 1–12 (2014).
Google Scholar 
Phalan, B. et al. Crop expansion and conservation priorities in tropical countries. PLoS ONE 8, e51759 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Afrane, Y. A. et al. Does irrigated urban agriculture influence the transmission of malaria in the city of Kumasi, Ghana?. Acta Trop. 89, 125–134 (2004).PubMed 

Google Scholar 
De Silva, P. M. & Marshall, J. M. Factors contributing to urban malaria transmission in sub-saharan Africa: A systematic review. J. Trop. Med. 2012, 1–10 (2012).
Google Scholar 
Kibret, S., Wilson, G. G., Tekie, H. & Petros, B. Increased malaria transmission around irrigation schemes in Ethiopia and the potential of canal water management for malaria vector control. Malar. J. 13, 1–12 (2014).
Google Scholar 
Dongus, S. et al. Urban agriculture and Anopheles habitats in Dar es Salaam, Tanzania. Geospat. Health 3, 189–210 (2009).PubMed 

Google Scholar 
Bhatt, S. et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature 526, 207–211 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Georganos, S. et al. Modelling the wealth index of demographic and health surveys within cities using very high-resolution remotely sensed information. Remote Sens. 11, 2543 (2019).ADS 

Google Scholar 
Pascual, M. & Baeza, A. What happens when forests fall?. Elife 10, e67863 (2021).PubMed 
PubMed Central 

Google Scholar 
MacDonald, A. J. & Mordecai, E. A. Amazon deforestation drives malaria transmission, and malaria burden reduces forest clearing. Proc. Natl. Acad. Sci. U. S. A. 116, 22212–22218 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Valle, D. & Clark, J. Conservation efforts may increase malaria burden in the Brazilian Amazon. PLoS ONE 8, e57519 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Luis, A. D., Kuenzi, A. J. & Mills, J. N. Species diversity concurrently dilutes and amplifies transmission in a zoonotic host–pathogen system through competing mechanisms. Proc. Natl. Acad. Sci. U. S. A. 115, 7979–7984 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Miller, E. & Huppert, A. The effects of host diversity on vector-borne disease: The conditions under which diversity will amplify or dilute the disease risk. PLoS ONE https://doi.org/10.1371/journal.pone.0080279 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Yamana, T. K. & Eltahir, E. A. B. Incorporating the Effects of Humidity in a Mechanistic Model of Anopheles Gambiae Mosquito Population Dynamics in the Sahel Region of Africa. http://www.parasitesandvectors.com/content/6/1/235. https://doi.org/10.1186/1756-3305-6-235 (2013).Wielgosz, B., Kato, E. & Ringler, C. Agro-ecology, household economics and malaria in Uganda: Empirical correlations between agricultural and health outcomes. Malar. J. 13, 1–11 (2014).
Google Scholar 
Asale, A., Duchateau, L., Devleesschauwer, B., Huisman, G. & Yewhalaw, D. Zooprophylaxis as a control strategy for malaria caused by the vector Anopheles arabiensis (Diptera: Culicidae): A systematic review. Infect. Dis. Poverty 6, 1–14 (2017).
Google Scholar 
Halliday, F., Rohr, J. & Laine, A.-L. Biodiversity loss underlies the dilution effect of biodiversity. Ecol. Lett. https://doi.org/10.1101/2020.04.20.050377 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Pienkowski, T., Dickens, B. L., Sun, H. & Carrasco, L. R. Empirical evidence of the public health benefits of tropical forest conservation in Cambodia: A generalised linear mixed-effects model analysis. Lancet Planet. Health 1, e180–e187 (2017).PubMed 

Google Scholar 
Shah, H. A., Huxley, P., Elmes, J. & Murray, K. A. Agricultural land-uses consistently exacerbate infectious disease risks in Southeast Asia. Nat. Commun. 10, 4299 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Parham, P. E. & Michael, E. Modeling the effects of weather and climate change on malaria transmission. Environ. Health Perspect. 118, 620–626 (2010).PubMed 

Google Scholar 
Krefis, A. C. et al. Modeling the relationship between precipitation and malaria incidence in children from a holoendemic area in Ghana. Am. J. Trop. Med. Hyg. 84, 285–291 (2011).PubMed 
PubMed Central 

Google Scholar 
Drakeley, C. J. et al. Altitude-Dependent and -Independent Variations in Plasmodium falciparum Prevalence in Northeastern Tanzania. J. Infect. Dis. 191, 1589–1598 (2005).PubMed 

Google Scholar 
Masuda, K. Length of maternal schooling and children’s risk of malaria infection: Evidence from a natural experiment in Uganda. BMJ Glob. Health 5, 4–11 (2020).
Google Scholar 
Ma, C. et al. Is maternal education a social vaccine for childhood malaria infection? A cross-sectional study from war-torn Democratic Republic of Congo. Pathog. Glob. Health 111, 98–106 (2017).PubMed 
PubMed Central 

Google Scholar 
Njau, J. D., Stephenson, R., Menon, M. P., Kachur, S. P. & McFarland, D. A. Investigating the important correlates of maternal education and childhood malaria infections. Am. J. Trop. Med. Hyg. 91, 509–519 (2014).PubMed 
PubMed Central 

Google Scholar 
Degarege, A., Fennie, K., Degarege, D., Chennupati, S. & Madhivanan, P. Improving socioeconomic status may reduce the burden of malaria in sub Saharan Africa: A systematic review and meta-analysis. PLoS ONE 14, 1–26 (2019).
Google Scholar 
Sonko, S. T. et al. Does socio-economic status explain the differentials in malaria parasite prevalence? Evidence from the Gambia. Malar. J. 13, 1–12 (2014).
Google Scholar 
Tusting, L. S. et al. Housing improvements and malaria risk in Sub-Saharan Africa: A multi-country analysis of survey data. PLoS Med. 14, 1–15 (2017).
Google Scholar 
Yang, D. et al. Drinking water and sanitation conditions are associated with the risk of malaria among children under five years old in sub-Saharan Africa: A logistic regression model analysis of national survey data. J. Adv. Res. 21, 1–13 (2020).PubMed 

Google Scholar 
Hay, S. I., Guerra, C. A., Tatem, A. J., Atkinson, P. M. & Snow, R. W. Urbanization, malaria transmission and disease burden in Africa. Nat. Rev. Microbiol. 3, 81–90 (2011).
Google Scholar 
Murray, C. J. L. et al. Global malaria mortality between 1980 and 2010: A systematic analysis. Lancet 379, 413–431 (2012).PubMed 

Google Scholar 
Nankabirwa, J. et al. Malaria in school-age children in Africa: An increasingly important challenge. Trop. Med. Int. Heal. 19, 1294–1309 (2014).
Google Scholar 
Okiro, E. A. et al. Age patterns of severe paediatric malaria and their relationship to Plasmodium falciparum transmission intensity. Malar. J. 8, 1–11 (2009).
Google Scholar 
Fullman, N., Burstein, R., Lim, S. S., Medlin, C. & Gakidou, E. Nets, spray or both? the effectiveness of insecticide-treated nets and indoor residual spraying in reducing malaria morbidity and child mortality in sub-Saharan Africa. Malar. J. 12, 1 (2013).
Google Scholar 
Agusto, F. B. et al. The impact of bed-net use on malaria prevalence. J. Theor. Biol. 320, 58–65 (2013).ADS 
PubMed 
MATH 

Google Scholar 
Hughes, R. A., Heron, J., Sterne, J. A. C. & Tilling, K. Accounting for missing data in statistical analyses: Multiple imputation is not always the answer. Int. J. Epidemiol. 48, 1294–1304 (2019).PubMed 
PubMed Central 

Google Scholar 
Beck-Johnson, L. M. et al. The effect of temperature on Anopheles mosquito population dynamics and the potential for malaria transmission. PLoS ONE 8, e79276 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Hien, D. F. D. S. et al. Plant-mediated effects on mosquito capacity to transmit human malaria. PLoS Pathog. 12, 1–17 (2016).
Google Scholar 
Donnelly, B., Berrang-Ford, L., Ross, N. A. & Michel, P. A systematic, realist review of zooprophylaxis for malaria control. Malar. J. 14, 1–16 (2015).
Google Scholar 
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).ADS 
CAS 
PubMed 

Google Scholar 
Townes, L. R., Mwandama, D., Mathanga, D. P. & Wilson, M. L. Elevated dry-season malaria prevalence associated with fine-scale spatial patterns of environmental risk: A case-control study of children in rural Malawi. Malar. J. 12, 1 (2013).
Google Scholar 
Brock, P. M. et al. Predictive analysis across spatial scales links zoonotic malaria to deforestation. Proc. R. Soc. B 286, 20182351 (2019).PubMed 
PubMed Central 

Google Scholar 
Zhang, Z., Manjourides, J., Cohen, T., Hu, Y. & Jiang, Q. Spatial measurement errors in the field of spatial epidemiology. Int. J. Health Geogr. 15, 1–12 (2016).
Google Scholar 
Rockström, J. et al. Managing water in rainfed agriculture: The need for a paradigm shift. Agric. Water Manag. 97, 543–550 (2010).
Google Scholar 
Rockström, J., Barron, J. & Fox, P. Water productivity in rain-fed agriculture: Challenges and opportunities for smallholder farmers in drought-prone tropical agroecosystems. Water Product. Agric. limits Oppor. Improv. 85199, 145–162. https://doi.org/10.1079/9780851996691.0145 (2009).Article 

Google Scholar 
Steinmann, P., Keiser, J., Bos, R., Tanner, M. & Utzinger, J. Schistosomiasis and water resources development: Systematic review, meta-analysis, and estimates of people at risk. Lancet. Infect. Dis. 6, 411–425 (2006).PubMed 

Google Scholar 
Halstead, N. T. et al. Agrochemicals increase risk of human schistosomiasis by supporting higher densities of intermediate hosts. Nat. Commun. 9, 837 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
Rohr, J. R. et al. Emerging human infectious diseases and the links to global food production. Nat. Sustain. 2, 445–456 (2019).PubMed 
PubMed Central 

Google Scholar 
Sokolow, S. H. et al. Reduced transmission of human schistosomiasis after restoration of a native river prawn that preys on the snail intermediate host. Proc. Natl. Acad. Sci. U. S. A. 112, 9650–9655 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rasolofoson, R. A., Hanauer, M. M., Pappinen, A., Fisher, B. & Ricketts, T. H. Impacts of forests on children’s diet in rural areas across 27 developing countries. Sci. Adv. 4, 1–10 (2018).
Google Scholar 
Doxsey-Whitfield, E. et al. Taking advantage of the improved availability of census data: A first look at the gridded population of the world, version 4. Pap. Appl. Geogr. 1, 226–234 (2015).
Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 850, 850–854 (2013).ADS 

Google Scholar 
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).
Google Scholar 
Hollister, M. J. Package ‘elevatr’ Title Access Elevation Data from Various APIs. (2018).Bontemps, S. et al. Consistent global land cover maps for climate modelling communities: Current achievements of the ESA’s land cover CCI. Proc. ESA Living Planet Symp. 13, 9–13 (2013).
Google Scholar 
Mahende, C. et al. Performance of rapid diagnostic test, blood-film microscopy and PCR for the diagnosis of malaria infection among febrile children from Korogwe District, Tanzania. Malar. J. 15, 1–7 (2016).
Google Scholar 
Stauffer, W. M. et al. Diagnostic performance of rapid diagnostic tests versus blood smears for malaria in US clinical practice. Clin. Infect. Dis. 49, 908–913 (2009).PubMed 

Google Scholar 
Yankson, R., Anto, E. A. & Chipeta, M. G. Geostatistical analysis and mapping of malaria risk in children under 5 using point-referenced prevalence data in Ghana. Malar. J. 18, 1–12 (2019).
Google Scholar 
Gatton, M. L. et al. Impact of Plasmodium falciparum gene deletions on malaria rapid diagnostic test performance. Malar. J. 19, 1–11 (2020).
Google Scholar 
Austin, K. F. Export agriculture is feeding malaria: A cross-national examination of the environmental and social causes of malaria prevalence. Popul. Environ. 35, 133–158 (2013).
Google Scholar 
Tyukavina, A. et al. Aboveground carbon loss in natural and managed tropical forests from 2000 to 2012. Environ. Res. Lett. 10, 074002 (2015).ADS 

Google Scholar 
Ayele, D. G., Zewotir, T. T. & Mwambi, H. G. Prevalence and risk factors of malaria in Ethiopia. Malar. J. 11, 1 (2012).
Google Scholar 
Acheson, E. S. & Kerr, J. T. Nets versus spraying: A spatial modelling approach reveals indoor residual spraying targets Anopheles mosquito habitats better than mosquito nets in Tanzania. PLoS ONE 13, 1–19 (2018).
Google Scholar 
Siraj, A. S. et al. Temperature and population density determine reservoir regions of seasonal persistence in highland malaria. Proc. R. Soc. B 282, 20151383 (2015).PubMed 
PubMed Central 

Google Scholar 
Ishengoma, D. S. et al. Trends of Plasmodium falciparum prevalence in two communities of Muheza district North-eastern Tanzania: Correlation between parasite prevalence, malaria interventions and rainfall in the context of re-emergence of malaria after two decades of progressive. Malar. J. 17, 1–10 (2018).
Google Scholar 
Weiss, D. J. et al. Air temperature suitability for Plasmodium falciparum malaria transmission in Africa 2000–2012: A high-resolution spatiotemporal prediction. Malar. J. 13, 1–11 (2014).
Google Scholar 
Watts, A. G. et al. Elevation as a proxy for mosquito-borne zika virus transmission in the Americas. PLoS ONE 12, 1–16 (2017).
Google Scholar 
Shah, H. A., Dritsaki, M., Pink, J. & Petrou, S. Psychometric properties of Patient Reported Outcome Measures (PROMs) in patients diagnosed with Acute Respiratory Distress Syndrome (ARDS). Health Qual. Life Outcomes 14, 15 (2016).PubMed 
PubMed Central 

Google Scholar 
Eneanya, O. A. et al. Environmental suitability for lymphatic filariasis in Nigeria. Parasites Vectors 11, 1–13 (2018).
Google Scholar 
Craney, T. A. & Surles, J. G. Model-dependent variance inflation factor cutoff values. Qual. Eng. 14, 391–403 (2002).
Google Scholar 
Anderson, D. & Burnham, K. Model Selection and Multimodel Inference (Springer, 2002).MATH 

Google Scholar 
Giam, X. & Olden, J. D. Quantifying variable importance in a multimodel inference framework. Methods Ecol. Evol. 7, 388–397 (2016).
Google Scholar 
Guo, G. & Zhao, H. Multilevel modeling for binary data. Annu. Rev. Sociol. 26, 441–462 (2000).
Google Scholar 
Li, B., Lingsma, H. F., Steyerberg, E. W. & Lesaffre, E. Logistic random effects regression models: A comparison of statistical packages for binary and ordinal outcomes. BMC Med. Res. Methodol. 11, 1–11 (2011).
Google Scholar 
Symonds, M. R. E. & Moussalli, A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav. Ecol. Sociobiol. 65, 13–21 (2011).
Google Scholar 
Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. Behav. Ecol. Sociobiol. 65, 23–35 (2011).
Google Scholar 
Shmueli, G. To explain or to predict?. Stat. Sci. 25, 289–310 (2010).MathSciNet 
MATH 

Google Scholar 
Hamlet, A. et al. The seasonal influence of climate and environment on yellow fever transmission across Africa. PLoS Negl. Trop. Dis. 12, 1–17 (2018).
Google Scholar 
Foley, J. A. et al. Global consequences of land use. Science 309, 570–574 (2005).ADS 
CAS 
PubMed 

Google Scholar 
Norton, E. C., Dowd, B. E. & Maciejewski, M. L. Marginal effects: Quantifying the effect of changes in risk factors in logistic regression models. JAMA 320, 84–85 (2018).PubMed 

Google Scholar 
RStudio Team. R Studio: Integrated Development for R (2015).Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Google Scholar  More

Since the bone discontinuities noted in the three vertebrae analyzed show no clear sign of bone overgrowth, it is pivotal to rule out the possibility that we are dealing with preservation damages before proposing an accurate diagnosis for the lesions. The close-up view examination of the abnormalities shows that their edges have clear signs of smoothing and rounding (Fig. 1), which represent important evidence of osteoblastic activity18,19. Additionally, the similar color of the cortical damage and normal bone can be used as secondary evidence to rule out post-mortem processes as a possible origin of the alterations, since recent destructive processes are lighter than the rest of the bone19. Therefore, as taphonomic processes can be ruled out, the pointed evidence strongly suggests that the discontinuities observed are of pathological origin. More specifically, these breaks found in all three vertebrae are indicative of bone fracture.Based on fracture analysis criteria applied here20, which consider the location and morphological pattern of the fractures, we classified the fractures noted in all vertebrae as traumas belonging to Type A (vertebral body compression), Group A2 (split fractures), and subgroup A2.1 (sagittal split fracture). This diagnosis implies that the traumatic episode was likely caused by a compressive force on the vertebral column, which split the vertebral bodies in the sagittal plane. This type of injury is considered stable—i.e., the fracture does not have a tendency to displace after reduction—and neurological deficit is uncommon20,22,23. Although stable traumas cause only moderate pain, without generating significant movement limitations20, the Eremotherium individual here analyzed died with unhealed bones, as there is no evidence of callus formation.The absence of other skeletal signs that point to the presence of another type of disease concomitantly to the fractures allows us to reject the possibility that they have been generated as a result of a pre-existing disease (e.g., infection, neoplasm). We also consider that the vertebral injuries were not caused by repetitive force (stress fractures) because this type of injury is commonly characterized as a nondisplaced line or crack in the bone, called hairline fracture3. Those refer to situations where the broken bone fragments are not visibly out of alignment and exhibit very little relative displacement21. Although the Eremotherium vertebrae fractures’ can be described as nondisplaced, they also have a noticeable gap between their edges that is mostly narrow with wider parts in the middle, something found in split fractures20 but that is not characteristic of hairline fractures. Lastly, the subgroup C1.2.1 (rotational sagittal split fracture) might be a source of confusion due to similar morphological pattern with subgroup A2.1 (sagittal split fracture). However, in subgroup C1.2.1 there are compressive and rotational forces acting simultaneously, producing total separation into two parts20, which clearly did not occur in the vertebrae analyzed here.In humans, compression fractures are most commonly caused by osteoporosis, although infection, neoplasm and trauma can also be etiological factors23,24,25. However, as aforementioned, the absence of other pathological skeletal marks is an important characteristic to take note as it serves to disregard the possibility of the fractures’ genesis to be secondary to another pathology. As such, in this case, osteoporosis, infection and neoplasm are unlikely etiologies. On the other hand, a compression fracture in a healthy individual is commonly generated after a severe traumatic event such as a fall from great height23,26. This scenario seems to better explain the origin of the vertebral fractures in the case of the Eremotherium ground sloth herein studied.The three fractured vertebrae were recovered in the Toca das Onças site (Fig. 2), a small cave considered as one of the richest paleontological sites of the Brazilian Quaternary15. Two complete skeletons of Eremotherium laurillardi and fragments belonging to at least thirteen other individuals, together with several other bones assigned to different smaller species are known to this cave14. It comprises of a single dry chamber that can only be entered through vertical entrances approximately 4.5 m high (Figs. 2b–d and 3). Two different hypotheses concerning the depositional process of Toca da Onças were previously proposed: (1) the animals climbed down into the cave in search of water14; or (2) due to the vertical character of the cave entrance, it could have functioned as a natural trap where animals accidentally fell into the cave15.Figure 2Location map of the Toca das Onças site and images of the cave. (a) Detail of the location, (b) cave entrance area view, (c) view from inside the cave, (d) Cave entrance detail. Scale bars 10 m in (b) and 5 m in (c). This figure was generated by Adobe Photoshop CS6 software (https://www.adobe.com/br/products/photoshop.html).Full size imageFigure 3Schematic representation of the Toca das Onças site. (a) Ground plan of the cave illustrating its morphology and dimension, (b) Cross-section illustrating the abyss-shaped entrance.Full size imageThe first hypothesis would indicate that the animal fell into the cave during an attempt to climb down. However, there is no report in the literature indicating that Eremotherium laurillardi could have been a climbing animal. In addition, the vertical morphology of the cave entrance would be a limiting factor for climbing behavior (see Fig. 3).Therefore, based on the type of fracture (compression sagittal split fracture) observed in the three vertebrae of Eremotherium as well as the inferred origin mechanism (fall from a great height), the presence of the individual here analyzed in the fossil accumulation of Toca das Onças is more likely explained by the second hypothesis. This idea is not particularly new as ‘entrapment due to fall’ has been described as a fossil accumulation mode to several other caves worldwide (e.g.,27,28). However, the use of bones fractures as an indicator of fossil accumulation mode is an interesting novelty. Of course, a detailed taphonomic investigation in the Toca das Onças still needs to be conducted in order to accurately interpret the formation of this important Quaternary fossil accumulation from Brazil.In sum, we suggest that the animal accidentally fell into the cave, fractured at least three sequential vertebrae (12th, 13th thoracic vertebrae and 1st lumbar vertebra) after the impact on the ground, survived for a while, but succumbed trapped inside the cave without food and water (Fig. 4). Other animals found in the cave, but without signs of bone fracture, may have fallen and not fractured their bones or not survived after the fall, especially the smaller ones. Finally, the proposal of falls to explain the unusual record of giant ground sloth fossils preserving much of its skeleton in caves, as reported for Toca das Onças site, contrasts with the better-documented pattern of skeletal accumulation via hydraulic action.Figure 4Artistic reconstruction of the suggested fall of the individual Eremotherium laurillardi into the cave. Artwork by Júlia d’Oliveira.Full size image More

Our compiled dataset consisted of 1674 records of marine mammal records after removing duplicate reports. It included 660 reports of sightings, 59 reports of induced mortalities or hunting records, 240 reports of incidental mortalities, 632 unique stranding records (live / dead), and 83 records which could not be categorised because of incomplete information.SightingsA total of 660 opportunistic sightings (number of individuals, ni = 3299) were recorded throughout the Indian coastline between 1748 and 2017 (Fig. 1a, 2a, 3a). Sighting data on the east coast (species = 18, ni = 1105) was mostly restricted to Odisha and Tamil Nadu (representing 97% of total east coast sightings). On the west coast (ni = 1297), Maharashtra (ni = 549), Gujarat (ni = 248) and Karnataka (ni = 307) contributed to highest sighting records (representing 85% of total west coast sightings). Sightings from the islands also contributed to 24.85% of the dataset (Andaman & Nicobar Islands = 24.37%, Lakshadweep = 0.48%). Highest incidence of sightings was for DFP (ni = 1894) followed by dugongs (ni = 959), BW (ni = 58) and SBW (ni = 17).Figure 1Marine mammal records obtained from data compiled between years 1748 – 2017 along the east coast, west coast and the islands of India for the groups i.e., baleen whales (BW), dolphins and finless porpoise (DFP), sperm and beaked whales (SBW) and dugongs, given as color-coded stacked bars where (a) sighting records—records where live animals were sighted (b) induced mortalities—records where animals were reported hunted or killed or were driven ashore, (c) incidental mortalities—records where animals were found dead after entanglement in fishing nets or being struck by vessels and (d) stranding records—records where dead or live animals were found washed ashore, or floating near shore or stranded alive and were attempted for rescue.Full size imageFigure 2Marine mammal records obtained every year from the data compiled between years 1748–2017 along Indian coastline given as cumulative numbers for each group i.e., baleen whales (BW), dolphins and finless porpoise (DFP), sperm and beaked whales (SBW) and dugongs, as color-coded stacked bars, where (a) sighting records—records where live animals were sighted (b) induced mortalities—records where animals were reported hunted or killed or were driven ashore, (c) incidental mortalities—records where animals were found dead after entanglement in fishing nets or being struck by vessels and (d) stranding records—records where dead or live animals were found washed ashore, or floating near shore or stranded alive and were attempted for rescue.Full size imageFigure 3Bubble plots showing distribution of marine mammal records obtained from data compiled between years 1748–2017 along the Indian coastline for each group i.e., baleen whales (BW), dolphins and finless porpoise (DFP), sperm and beaked whales (SBW) and dugongs, as color-coded stacked bars, where (a) sighting—records where live animals were sighted (b) induced mortalities—records where animals were reported hunted or killed or were driven ashore, (c) incidental mortalities—records where animals were found dead after entanglement in fishing nets or being struck by vessels and (d) strandings—records where dead or live animals were found washed ashore, or floating near shore or stranded alive and were attempted for rescue. Size of the bubble indicates number of individuals. These maps were created using ArcGIS 10.5 (https://desktop.arcgis.com/en/arcmap/10.3/map/working-with-layers/about-symbolizing-layers-to-represent-quantity.htm).Full size imageInduced mortalitiesA total of 59 incidences (ni = 102) were recorded of marine mammals being hunted/ captured between the years 1748–2017 (Fig. 1b, 2b, 3b). The total number of animals hunted/ captured deliberately is similar along east coast (ni = 33), west coast (ni = 29) and islands (ni = 36). Out of all marine mammal species, 90% of the animals hunted at the east coast were dugong D. dugon (ni = 30, all from Tamil Nadu). On the west coast, records of hunting incidences of finless porpoise Neophocaena phocaenoides were highest (79% of total records on west coast, Goa ni = 17, Kerala ni = 4, Karnataka and Maharashtra ni = 1). In the islands (i.e., Andaman and Nicobar Islands), 94% of the hunting records were of dugongs (ni = 34).Incidental mortalitiesA total of 240 net entanglements (ni = 1356) were reported along the Indian coast between the years 1748 and 2017 (Fig. 1c, 2c, 3c). Similar counts of individuals entangled along east (ni = 670) and west coast (ni = 654) were obtained with low reporting from the islands (ni = 26). Fourteen species were reported entangled from both east and west coast with only 4 species recorded from the islands. D. dugon was found to be most frequently entangled along the east coast (63 incidences, ni = 594, contributing to 56% of the total numbers on east coast), followed by Tursiops sp. (11 incidences, ni = 14, 9% of the east coast dataset). On the west coast, Tursiops sp. was the most frequently entangled (18 incidences, ni = 117, contributing to 18% of the west coast dataset), followed by N. phocaenoides (17 incidences, ni = 34, contributing to 17% of the dataset). The total number of DFP being entangled from west coast (ni = 623) were higher than east coast (ni = 68). More dugong individuals were entangled along east coast (i.e., from Tamil Nadu; ni = 594) as compared to the west coast (i.e., Gujarat; ni = 3) and Islands (i.e., Andaman and Nicobar; ni = 19). D. dugon was the most frequently entangled species in the islands (19 incidences, ni = 19, contributing to 79% of the total numbers in islands dataset) followed by false killer whale Pseudorca crassidens (3 incidences, ni = 5, contributing to 12% of the islands dataset). Very few BW or SBW (11 incidences, ni = 11) were recorded accidently entangled throughout the Indian coastline.StrandingsMarine mammals stranding reports consisted of 91.93% dead (ni = 581) and 8.07% live strandings (ni = 51) (Figs. 1d, 2d, 3d). Considering mass strandings as strandings with ni  > 2 (excluding mother and calf;33,34), 8.5% of all reports were mass strandings (21 strandings, ni = 1054). Most of the records did not have information about the sex of the stranded animal (83%), the age class (88%) or the state of decomposition of the carcass (53%). Highest strandings were reported of dugongs (strandings = 190, ni = 228), followed by BW (strandings = 178, ni =  = 190), DFP (strandings = 157, ni =  = 552) and SBW (strandings = 47, individuals = 48). There were 54 incidences (ni = 54, 9% of total stranding data) where the animal was not identified reliably to include in either of the groups.Species composition and frequencies of strandings were different on east coast, west coast and in the islands (Fig. 1, Table 1). Twenty-two species were reported as stranded on the east coast with D. dugon as the most frequently stranded species (83 incidences, ni = 107, ~ 29% of all records), followed by Indo-Pacific humpback dolphin Sousa chinensis, (31 incidences, ni = 108, ~ 10% of all records). On the west coast, out of 20 species reported as stranded, Balaenoptera musculus was most frequent (28 incidences, ni = 29, ~ 12% of all records) followed by N. phocaenoides (23 incidences, ni = 39, ~ 10% of all records). In the islands, 13 species were reported as stranded, D. dugon (93 incidences, ni = 102, contributing to 77% of the total animals found on the islands) followed by strandings of sperm whale Physeter macrocephalus (8 incidences, ni = 8, contributing to 6% of the data; Table 1).

a. Baleen whales

Table 1 Number of stranding events reported for marine mammals between 1748–2017 in India from the east coast, the west coast and Lakshadweep and Andaman & Nicobar archipelagos.Full size tableA total of 178 BW strandings (ni = 190) were reported. Most species were unidentified (east coast ni= 27, west coast ni = 58, islands ni = 4; i.e., 47% of the data). Identified strandings comprised of 6 species (see Table 1), some of which were later found to be misidentification (no confirmed evidence for common Minke Whale Balaenoptera acutorostrata, Sei Whale Balaenoptera borealis and Fin Whale Balaenoptera physalus from Indian waters; MMRCNI, 2018). Higher number of strandings occurred on the west coast (ni = 126), as compared to east coast (ni = 60). The east and west coast reported all six species of BW, whereas only three species stranded on the islands. B. borealis (misidentified) was the most stranded species across the east coast (12 incidences, ni = 12, contributing to 11% of the data) whereas blue whale Balaenoptera musculus was the most frequent across the west coast (28 incidences, ni = 29, contributing to 11% of the data). Baleen whale strandings were rare in the islands (4 incidences, ni = 4).Forty-seven SBW strandings (ni = 48) were reported along the Indian coast. More SBW stranded on the east coast (ni = 23) as compared to the west coast (ni = 13) and the islands (ni = 12). P. macrocephalus was most frequently reported (70% of all SBW records, east coast ni = 20, west coast ni = 6, islands ni = 8).There were 157 strandings (ni =552) of DFP belonging to 14 species. Twenty-one of these events were mass strandings (ni  > 2). The largest mass stranding event (ni = 147) occurred of short-finned pilot whale Globicephala macrorhynchus along the west coast (Tamil Nadu). Higher number of DFP strandings were recorded from east coast (ni = 418) as compared to west coast (ni = 83) and the islands (ni = 51; Table 1). East coast received a higher diversity of stranded DFP (number of species = 11) as compared to west coasts (number of species = 9) and the islands (number of species = 3). S. chinensis was the most frequently stranded species along the east coast (31 incidences, ni = 108, contributing to 33% of the data) whereas N. phocaenoides was the most frequent along the west coast (23 incidences, ni = 39, contributing to 37% of the data; Table 1).

d. Dugongs

The current distribution of dugongs in India is in the shallow coastal waters of Gujarat, Tamil Nadu and Andaman & Nicobar Islands37,38. There are 190 stranding events recorded between the years 1893 and 2017. The highest number of stranded dugongs were recorded from Tamil Nadu (ni = 107) closely followed by Andaman and Nicobar Islands (ni = 102) and few records from Gujarat (ni = 19).Temporal stranding patternsOur analysis of temporal trends for the last 42 years (1975–2017) showed that the mean number of strandings along the Indian coast was 11.25 ± SE 1.39 / year. The number of stranding reports show an increasing trend for two decades after 1975, dropping between 1995 and 2004. We observed a distinct rise in strandings post 2005 (18.23 ± SE 2.98 / year) with the highest reports from 2015–17 (27.66 ± SE 8.51/year) (Fig. 4).

a. Baleen whales

Figure 4A beanplot of decadal trends in marine mammal stranding in India from data compiled between years 1975–2017. Data prior to 1975 was discontinuous over the years to be considered for decadal trends. The data for last decade considered here includes only two years (2015–17) where increased reporting is evident. The bold horizontal lines indicate the mean number of strandings in each decade whereas the smaller horizontal lines indicate stranding numbers recorded for each year within the decade.Full size imageOn the west coast, mean stranding rate throughout the years (1975–2017) was 0.0010 ± SE 0.0014 strandings/km, and a steady rise was observed in rate of reported strandings after 2010. A seasonal trend was observed as well, with a peak in the month of September (sr = 0.0061 ± SE 0.0016 strandings/km), i.e., towards the end of monsoon season, and lowest strandings were recorded in the month of June (sr = 0.0016 ± SE 0.006 strandings/ km) (Fig. 5).Figure 5Temporal patterns (annual and monthly stranding rates / 100 km of coastline) in strandings of marine mammal records obtained from data compiled between years 1975–2017 along east and west coast of India for each group where (a) annual stranding rate and (b) monthly stranding rate for baleen whales (BW); (c) annual stranding rate and (d) monthly stranding rate for dolphins and finless porpoise (DFP); (e) annual stranding rate and (f) monthly stranding rate for sperm and beaked whales (SBW) and (g) annual stranding rate and (h) monthly stranding rate for dugongs.Full size imageThe mean stranding rate of BW on the east coast through 1975–2017 was 0.0013 ± SE 0.0017 strandings/km, but no specific trends were observed according to years or seasons. Stranding rates of BW did not differ between east and west coast (Mann–Whitney U test, U = 390, U standardized = -0.025, p value  > 0.05).The stranding rates of SBW differed significantly along both the coasts (Mann Whitney U test, U = 192, U standardized = 0.0, p value  More

Present fire-climate classificationTo identify the different regions of the planet with suitable climatic conditions for fire activity, we compare the global distribution of climate indicators based on temperature and precipitation, with satellite-derived GFED4 burned area data21 (Fig. 1). Starting from four general climates (Tr-tropical, Ar-arid, Te-temperate and Bo-boreal) based on the Köppen–Geiger climate classification main categories22, we create four fire-prone classes using climate thresholds to define the patterns observed in Fig. 1. Each category is characterised by the element that boosts fire activity during the FS: low precipitation, high temperatures or a combination of both. The classification is made by contrasting the probability distribution of the climatic variables at data points associated with high fire activity vs. points with low fire activity within the main Köppen-Geiger categories (see Threshold Selection in Methods section for a detailed explanation).Fig. 1: Burned area observations and climate drivers.a 1996–2016 maximum annual burned area (BAmax) and monthly burned area time series for selected regions. b Average monthly precipitation percentage from the annual total for the fire season (PPFS). c Average monthly temperature anomaly from the annual mean for the fire season (TAFS).Full size imageThe environmental conditions associated with fire occurrence emerge more clearly in this comparison, yielding the different threshold sets in Table 1 that determine the fire-prone months at any location (the selection method is detailed in the Methods section). We define those years with at least 1-month meeting the thresholds, as fire-prone years (FPY). Depending on the number of FPY at each location, the categories of Table 1 are sub-divided into recurrent (r), occasional (o) and infrequent (i) (see Methods). The average number of fire-prone months during the FPY is defined as the potential FS length (PFSL), i.e., the season with climatic characteristics prone to fire activity.Table 1 Fire classification defining criteria.Full size tableFigure 2a depicts the global map of the burned areas classified according to the selected thresholds (Table 1). Savanna fires are responsible for the largest proportion of burned area on the global scale21. The FS in these areas is longer than in other climates (see Supplementary Fig. 1) and, despite savanna fires being also dependent on ignition patterns and human policies and practices, the FS is tied to a pronounced seasonal cycle of precipitation23,24,25, with fire occurring mainly during the dry part of the cycle. Because of this, the Tropical – dry season fire class (Tr-ds) coincides with the distribution of the tropical savanna climate. In Fig. 2, boreal fires are represented as hot season fires (Bo-hs) due to the large positive temperature anomaly existing in those locations during the FS (Fig. 1c). In fact, temperature variations explain much of the variability in boreal burned area26,27. Temperate fires are classified as dry and hot season (Te-dhs) because they affect regions where the dry season coincides with the warm season (Fig. 1b, c). Here, high temperatures and precipitation seasonality determine fire activity and inter-annual burned area variability, e.g., in Western North America28,29,30,31 and Southern Europe32,33. Fire activity in arid regions occurs during warm months, but the relation with precipitation is more complex. The FS is associated with a hot season in cooler (MAT  27.5 °C), the FS starts right at the beginning of the dry season (e.g., the Sahel, Supplementary Fig. 12) while where MATs are more moderate, between 18.5 and 27.5 °C, it takes longer to develop (e.g., Central Australia and the Kalahari desert, Supplementary Figs. 12 and 13). Due to the dependency between fires and the existence of fuel in arid climates, we named this class Arid fuel limited (Ar-fl). A more in-depth discussion about the definition of this fire-climate class can be found in the section entitled Threshold selection for each climate of the Supplementary Information.Fig. 2: Fire-prone region classification.a With observed burned area data as a reference: not classified (NC, white) and misclassified (C, black) areas with BAmax = 0 ha, unclassified (NC, grey) and classified (Tr-ds, Ar-fl, Te-dhs and Bo-hs) areas with BAmax  > 0 ha. Each class is subdivided into three subcategories depending on the recurrence of the fire-prone conditions: recurrent (r), occasional (o) and infrequent (i). b Present (1996–2016) fire-prone climatic regions. c Future (2070–2099) fire-prone climatic regions with shaded grey representing a  0 ha) or fireless (BA = 0 ha). This reveals a two-way relation between fires and climate: fires take place under specific climatic conditions, and most places with these climatic conditions are indeed fire-prone, which supports our earlier hypothesis. Fire activity is controlled by weather, resources to burn and ignitions, as represented through the fire regime triangle12,20. On broad temporal scales and large spatial scales, temperature and precipitation have an important impact on fire because these climate variables influence vegetation type and the abundance, composition, moisture content, and structure of fuels34. Although ignitions may be driving fires to a greater extent than temperature or precipitation at specific locations or events35, they do not seem to limit fire activity at coarse spatial and temporal resolutions, implying that where fuels are sufficient and atmospheric conditions are conducive to combustion, the potential for ignition exists, either by lightning or human causes13,20. For all these reasons, we can identify specific climates that are prone to fires.The areas classified as fire-prone in Fig. 2b comprise 99.26% of the observed global mean annual burned area in Supplementary Fig. 2. This percentage is above 85% for all four general climates (Supplementary Fig. 20). The percentage of land area with non-zero burned area data classified as fire-prone is 91.22%. Considering for each location only the obtained FPY, the percentage of the observed burned area classified is 90.36%. Furthermore, the PFS obtained in the fire-climate classification (Fig. 3b) also correlates well with the timing of observed fire incidence, as globally 87.91% of the observed mean burned area occurs during the identified months of PFS at classified fire-prone locations.Fig. 3: Potential fire season.a Future minus present potential fire season length (PFSL) difference in months (ΔPFSL). b Present potential fire season. c Future potential fire season.Full size imageUnclassified regions (in grey in Fig. 2a) correspond for the most part to those with the least burned area or those where agricultural practices modify the climatic seasonality of fires. In addition, as the classification is conceived from a climatic point of view, locations with fire activity associated with specific meteorological conditions that are not appreciable at the monthly temporal resolution, are probably not well identified. For example, a week of extremely high temperatures could be almost unnoticeable in the monthly mean temperature, but not in fire activity. Similarly, months with the same total precipitation may have different fire activity if the precipitation falls concentrated in a few days or is distributed throughout the month. Furthermore, the short temporal sampling period of the burned area data could also be influencing our results. Locations with long fire cycles may not be well represented in the data.Future fire-climate classificationA future fire-climate classification map is derived by applying the thresholds obtained in the present fire-climate classification to future climatology variables from multiple coupled model intercomparison project phase 5 (CMIP5) global circulation model (GCM) outputs, considering the RCP8.5 scenario (the worst-case climate change scenario of the CMIP5). Two contrasting approaches can be taken for analysing future fire activity, one that considers quick vegetation adaptation to the new climatic conditions, and another that does not. These two approaches clearly diverge in the boreal regions, where the biome (mainly taiga) is strongly conditioned by the low temperatures and where future temperature changes at the end of the 21st century will have a greater amplitude. It is expected that the boreal forest of these areas will not be immediately replaced by a temperate mixed forest where the average annual temperature exceeds the range of values typical of the taiga biome. Terrestrial vegetation compositional and structural change could occur during the 21st century where vegetation disturbance is accelerated or amplified by human activity, but equilibrium states may not be reached until the 22nd century or beyond36.Based on the assumption that during the future period (2070–2099) the vegetation will not be fully adapted to the new climatic conditions, and since the present Köppen–Geiger climate classification (on which we base our Tr, Ar, Te and Bo categories) closely corresponds to the different existent biomes22, we analyse only the projected changes in the specific fire-climate classification variables, maintaining the general division of Tropical, Arid, Temperate and Boreal regions as is in present climate conditions. The future fire-climate classification is shown in Fig. 2c.We note that we determine future fire activity from relationships of the latter with the present climate; however, these relationships might not be stationary. Our approach does not contemplate possible future changes in precipitation frequency if they are not noticeable in monthly precipitation amounts. Areas with the rising incidence of extreme precipitation events due to global warming37 may experience an increase in monthly precipitation but a decrease in rainy days, which may lead us to consider the conditions there less favourable for fire activity than they actually will be.Future changes in global fire activityModelled future fire-prone regions experience significant variations with respect to the present (Fig. 2b, c). Due to global warming, the Bo-hs fire class pertaining to boreal forests would spread over a larger area, reaching most of Northern Scandinavia and undergoing a southward and northward expansion in Canada, Alaska and Russia. This category may experience a percentual expansion of 47.0% according to our results. This expansion is more accentuated for the combination of the highest recurrence subcategories Bo-hs-r and Bo-hs-o, reaching a value of 111.5%.The conjunction of Te-dhs-r and Te-dhs-o fire classes of midlatitudes also undergoes a considerable expansion of 24.5% in the area (Fig. 2b, c). The most remarkable changes are expected in Southern China and Southern Europe. A large part of Europe transitions from an infrequent fire category to a more frequent fire category with Csa and Csb Mediterranean climates38.The Tr-ds fire classes with frequent fire-prone conditions in the Tropics presents fewer spatial changes (Fig. 2b, c), with a spatial contraction of 6.3%. The most important differences are found in South America. Some of the climate model results considered here indicate also that some parts of the Eastern Amazon rainforest will move from a non-fire class to Tr-ds fire class, as other studies have suggested39.The Arid fire-prone classes Ar-fl-r and Ar-fl-o would increase its area by 5.0%. Projected changes in the extent of this class are very sensitive to changes in annual precipitation, conducive to vegetation and fuel reduction or increment, thus there is significant uncertainty in the proximity of desert regions (Fig. 2c).Clearer conclusions can be drawn from the FPY and PFSL calculation (Figs. 3 and 4). The number of months meeting the set of conditions in Table 1 yields the estimated PFSL (Fig. 3b), and the number of years with at least 1-month meeting the thresholds, the FPY. In the boreal regions, we obtain a general lengthening of the PFS. The PFS of these areas is conditioned by temperature, so the amplified warming of Artic zones40 is expected to make the FS longer. Notwithstanding, in certain parts of Eastern Asia, the intense warming is counterbalanced by an increase of the precipitation in certain warm months (see Supplementary Figs. 21 and 22), leading to a slight shortening of our estimated PFS. There is evidence, however, that temperature increases may lead to drier fuels in the future despite the precipitation increase, thus augmenting fire risk, as some investigations have shown for Canada41. Our results agree in general with several other studies that have previously pointed towards an increase of the FSL in boreal areas1,17,42, even when some suggest a more pronounced lengthening in more northerly latitudes1,17. In terms of the frequency of years with fire-prone conditions, the conclusions are even clearer. A general increase of the FPY is observed, especially for northerly latitudes, where the differences reach values of more than +4 years per decade (Fig. 4a). This possible increase in fire activity in boreal areas may result in significant peatland combustion and a release of the large quantities of soil carbon that they store into the atmosphere43. These greenhouse gas emissions may create a positive feedback loop, leading to a further increase in temperature, which in turn will enhance boreal wildfire incidence and more peatland burning.Fig. 4: Fire-prone years.a Future minus a present number of years with at least one month classified as fire-prone per decade (ΔFPY). b Present fire-prone years per decade. c Future fire-prone years per decade.Full size imageThe Te-dhs fire class, corresponding to temperate climates, would also experience a general lengthening of the PFS (Fig. 3). A future precipitation decline may be especially significant in Southern Europe (Supplementary Fig. 21), associated with an increased anticyclonic circulation yielding more stable conditions44, while the temperature rise would be quite homogeneous among all Te-dhs fire-climate class areas. The FS drought intensification around the Mediterranean, together with the general warming (Supplementary Fig. 21), would lead to a lengthening of the PFS of around 2 months (Fig. 3a), but summer months could also experience this precipitation decline (Supplementary Fig. 22), meaning that the FS would be more severe. The Western US, which has already experienced over the last decades the lengthening of the FS45 and the increase of large fires46 and extreme wildfire weather47,48 due to climate change, may also experience an FS lengthening by the end of the 21st century. Some authors18,48,49,50 have studied projected fire future changes from other points of view (occurrence of very large fires, wildfire potential, etc.), finding also a general increase of fire severity by the end of the century in some of these Te-dhs fire regions. The interannual recurrence of fire-prone conditions will significantly increase in countries like France, Italy or Eastern China (Fig. 4a).The PFSL of the Tropical Tr-ds fire-climate class presents slight differences between present and future values (Fig. 3). Some areas of the Northern African savanna may experience a shortening of the PFS, while Southern Africa shows a lengthening. A dipole pattern of wetting in tropical Eastern Africa and drying in Southern Africa51 could be the reason for these future changes. There is a contrasting influence of ENSO in present African fire patterns52, which suggests that the future pattern of precipitation variations in Central Africa may be associated with ENSO future changes under climate change conditions53. Although the quantification of ENSO changes in a warmer climate is still an issue that continues to be investigated, an expansion and strengthening of ENSO teleconnections is confirmed by some authors53,54,55. The general increase in precipitation along all seasons in western equatorial Africa would lead to a significant decrease in the recurrence of interannual fire-prone conditions (Fig. 4a).Our results show that fire-prone areas in Temperate and especially Boreal climates are projected to undergo the most significant expansion and lengthening of the potential FS at the end of the XXI century driven by rising temperatures. In the Tropics, little change is expected in these respects. Notwithstanding, global warming is likely to make fire risk more severe mostly everywhere, and in particular in some regions such as Mediterranean Europe and the Eastern Amazon, where an important decrease in precipitation is also predicted during the PFS. More favourable fire conditions will potentially increment fire activity and burned areas in many places. In others, especially in the Tropics, increasing suppression efforts and a cease to agricultural and pastoral practices like vegetation clearing by fire, replaced by more intensive farming, could counteract the impact of a warmer climate. A reduction of these human-caused fires in the Tropics could bring global burned area down2, despite rising trends elsewhere, given the vast contribution of Tropical fires to the burned areas at the global scale (Fig. 1). More

Yao C-L, Somero GN. The impact of ocean warming on marine organisms. Chin Sci Bull. 2014;59:468–79.
Google Scholar 
Frölicher TL, Fischer EM, Gruber N. Marine heatwaves under global warming. Nature. 2018;560:360–4.PubMed 

Google Scholar 
Bindoff NL, Cheung WWL, Kairo JG, Arístegui J, Guinder VA, Hallberg R, et al. Changing ocean, marine ecosystems, and dependent communities. Switzerland: Intergovernmental Panel on Climate Change (IPCC); 2019.Breitburg D, Levin LA, Oschlies A, Grégoire M, Chavez FP, Conley DJ, et al. Declining oxygen in the global ocean and coastal waters. Science. 2018;359:eaam7240.PubMed 

Google Scholar 
Collins M, Knutti R, Arblaster J, Dufresne J-L, Fichefet T, Friedlingstein P, et al. Long-term climate change: projections, commitments and irreversibility. 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. United Kingdom and New York, NY, USA: Cambridge; 2013.Mackenzie BR, Schiedek D. Daily ocean monitoring since the 1860s shows record warming of northern European seas. Glob Change Biol. 2007;13:1335–47.
Google Scholar 
Gruner DS, Bracken MES, Berger SA, Eriksson BK, Gamfeldt L, Matthiessen B, et al. Effects of experimental warming on biodiversity depend on ecosystem type and local species composition. Oikos. 2017;126:8–17.
Google Scholar 
Forsman A, Berggren H, Åström M, Larsson P. To what extent can existing research help project climate change impacts on biodiversity in aquatic environments? A review of methodological approaches. Multidiscipl Digital Publishing Inst. 2016;4:75.
Google Scholar 
HELCOM. Eutrophication in the Baltic Sea—An integrated thematic assessment of the effects of nutrient enrichment and eutrophication in the Baltic Sea region. Baltic Sea Environ Proc. 2009. Report No.: 115B.Carstensen J, Andersen JH, Gustafsson BG, Conley DJ. Deoxygenation of the Baltic Sea during the last century. Proc Natl Acad Sci USA. 2014;111:5628–33.CAS 
PubMed 
PubMed Central 

Google Scholar 
Broman E, Sjostedt J, Pinhassi J, Dopson M. Shifts in coastal sediment oxygenation cause pronounced changes in microbial community composition and associated metabolism. Microbiome. 2017;5:96.PubMed 
PubMed Central 

Google Scholar 
Schmidtko S, Stramma L, Visbeck M. Decline in global oceanic oxygen content during the past five decades. Nature. 2017;542:335–9.CAS 
PubMed 

Google Scholar 
Brewer PG, Peltzer ET. Depth perception: the need to report ocean biogeochemical rates as functions of temperature, not depth. Philos Trans R Soc Mathemat Phys Eng. 2017;375:20160319.
Google Scholar 
Laruelle GG, Cai W-J, Hu X, Gruber N, Mackenzie FT, Regnier P. Continental shelves as a variable but increasing global sink for atmospheric carbon dioxide. Nat Commun. 2018;9:454.PubMed 
PubMed Central 

Google Scholar 
Gilbert D, Rabalais NN, Díaz RJ, Zhang J. Evidence for greater oxygen decline rates in the coastal ocean than in the open ocean. Biogeosciences. 2010;7:2283–96.CAS 

Google Scholar 
Kauppi L, Norkko J, Ikonen J, Norkko A. Seasonal variability in ecosystem functions: quantifying the contribution of invasive species to nutrient cycling in coastal ecosystems. Marine Ecol Progr Series. 2017;572:193–207.CAS 

Google Scholar 
Lu X, Zhou F, Chen F, Lao Q, Zhu Q, Meng Y, et al. Spatial and seasonal variations of sedimentary organic matter in a subtropical bay: implication for human interventions. Int J Environ Res Public Health. 2020;17:1362.CAS 
PubMed Central 

Google Scholar 
Turner JT. Zooplankton fecal pellets, marine snow, phytodetritus and the ocean’s biological pump. Progr Oceanograph. 2015;130:205–48.
Google Scholar 
Gupta A, Gupta R, Singh RL. Microbes and environment. In: Singh R (eds) Principles and Applications of Environmental Biotechnology for a Sustainable Future. Applied Environmental Science and Engineering for a Sustainable Future. Springer, Singapore; 2017:43–84.American Society for Microbiology. Microbes and Climate Change: Report on an American Academy of Microbiology and American Geophysical Union Colloquium held in Washington, DC, in March 2016. Washington (DC): American Society for Microbiology; 2017.Sarmento H, Montoya JM, Vazquez-Dominguez E, Vaque D, Gasol JM. Warming effects on marine microbial food web processes: how far can we go when it comes to predictions? Philos Trans R Soc B Biol Sci. 2010;365:2137–49.
Google Scholar 
IPCC. Climate Change 2021: The physical science basis. Contribution of working group I to the sixth assessment report of the Intergovernmental Panel on Climate Change. Cambridge University Press (In Press); 2021.Moberg A, Humborg C. Second assessment of climate change for the Baltic Sea basin. Second assessment of climate change for the Baltic Sea basin. Berlin Heidelberg: Springer; 2008.
Google Scholar 
Humborg C, Geibel MC, Sun X, McCrackin M, Mörth C-M, Stranne C, et al. High emissions of carbon dioxide and methane from the coastal Baltic Sea at the end of a summer heat wave. Front Marine Sci. 2019;6:493.
Google Scholar 
Smith TP, Thomas TJH, García-Carreras B, Sal S, Yvon-Durocher G, Bell T, et al. Community-level respiration of prokaryotic microbes may rise with global warming. Nat Commun. 2019;10:5124.PubMed 
PubMed Central 

Google Scholar 
Broman E, Li L, Fridlund J, Svensson F, Legrand C, Dopson M. Spring and late summer phytoplankton biomass impact on the coastal sediment microbial community structure. Microbial Ecol. 2018;77:288–303.
Google Scholar 
Gao Y, Cornwell JC, Stoecker DK, Owens MS. Influence of cyanobacteria blooms on sediment biogeochemistry and nutrient fluxes. Limnol Oceanograph. 2014;59:959–71.CAS 

Google Scholar 
Sawicka JE, Brüchert V. Annual variability and regulation of methane and sulfate fluxes in Baltic Sea estuarine sediments. Biogeosciences. 2017;14:325–39.CAS 

Google Scholar 
Berner RA. A new geochemical classification of sedimentary environments. J Sediment Res. 1981;51:359–65.CAS 

Google Scholar 
Nealson KH. Sediment bacteria: who’s there, what are they doing, and what’s new? Ann Rev Earth Planet Sci. 1997;25:403–34.CAS 

Google Scholar 
EPA. Quality criteria for water. Washington D.C., USA: Office of Water Regulations and Standards; 1986.
Google Scholar 
Tamme R, Hiiesalu I, Laanisto L, Szava-Kovats R, Pärtel M. Environmental heterogeneity, species diversity and co-existence at different spatial scales. J Veget Sci. 2010;21:796–801.
Google Scholar 
Klier J, Dellwig O, Leipe T, Jürgens K, Herlemann DPR. Benthic bacterial community composition in the oligohaline-marine transition of surface sediments in the Baltic Sea based on rRNA analysis. Front Microbiol. 2018;9:236.PubMed 
PubMed Central 

Google Scholar 
Broman E, Sachpazidou V, Pinhassi J, Dopson M. Oxygenation of hypoxic coastal Baltic Sea sediments impacts on chemistry, microbial community composition, and metabolism. Front Microbiol. 2017;8:2453.PubMed 
PubMed Central 

Google Scholar 
Orlygsson J, Kristjansson JK. The family Hydrogenophilaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The Prokaryotes: Alphaproteobacteria and Betaproteobacteria. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. p. 859–68.Liu Z, Frigaard NU, Vogl K, Iino T, Ohkuma M, Overmann J, et al. Complete genome of Ignavibacterium album, a metabolically versatile, flagellated, facultative anaerobe from the phylum Chlorobi. Front Microbiol. 2012;3:185.CAS 
PubMed 
PubMed Central 

Google Scholar 
Watanabe M, Kojima H, Fukui M. Desulfoplanes formicivorans gen. nov., sp. nov., a novel sulfate-reducing bacterium isolated from a blackish meromictic lake, and emended description of the family Desulfomicrobiaceae. Int J Syst Evol Microbiol. 2015;65:1902–7.CAS 
PubMed 

Google Scholar 
Galushko A, Desulfocapsaceae JK. Bergey’s Manual of Systematics of Archaea and Bacteria. Hoboken, New Jersey: Wiley; 2015. p. 1–6.
Google Scholar 
Dyksma S, Bischof K, Fuchs BM, Hoffmann K, Meier D, Meyerdierks A, et al. Ubiquitous Gammaproteobacteria dominate dark carbon fixation in coastal sediments. ISME J. 2016;10:1939–53.CAS 
PubMed 
PubMed Central 

Google Scholar 
Ye Q, Wu Y, Zhu Z, Wang X, Li Z, Zhang J. Bacterial diversity in the surface sediments of the hypoxic zone near the Changjiang Estuary and in the east China Sea. Microbiologyopen. 2016;5:323–39.CAS 
PubMed 
PubMed Central 

Google Scholar 
Fahrbach M, Kuever J, Remesch M, Huber BE, Kampfer P, Dott W, et al. Steroidobacter denitrificans gen. nov., sp. nov., a steroidal hormone-degrading gammaproteobacterium. Int J Syst Evol Microbiol. 2008;58:2215–23.CAS 
PubMed 

Google Scholar 
Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, et al. Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraeota (phyl. nov.). Front Microbiol. 2017;8:682.PubMed 
PubMed Central 

Google Scholar 
Reyes C, Schneider D, Thürmer A, Kulkarni A, Lipka M, Sztejrenszus SY, et al. Potentially active iron, sulfur, and sulfate reducing bacteria in Skagerrak and Bothnian bay sediments. Geomicrobiol J. 2017;34:840–50.CAS 

Google Scholar 
Lovley DR, Roden EE, Phillips EJP, Woodward JC. Enzymatic iron and uranium reduction by sulfate-reducing bacteria. Marine Geol. 1993;113:41–53.CAS 

Google Scholar 
Funkey CP, Conley DJ, Reuss NS, Humborg C, Jilbert T, Slomp CP. Hypoxia sustains cyanobacteria blooms in the Baltic sea. Environ Sci Technol. 2014;48:2598–602.CAS 
PubMed 
PubMed Central 

Google Scholar 
Boden R, Hutt LP, Rae AW. Reclassification of Thiobacillus aquaesulis (Wood & Kelly, 1995) as Annwoodia aquaesulis gen. nov., comb. nov., transfer of Thiobacillus (Beijerinck, 1904) from the Hydrogenophilales to the Nitrosomonadales, proposal of Hydrogenophilalia class. nov. within the ‘Proteobacteria’, and four new families within the orders Nitrosomonadales and Rhodocyclales. Int J Syst Evol Microbiol. 2017;67:1191–205.CAS 
PubMed 

Google Scholar 
Howarth R, Unz RF, Seviour EM, Seviour RJ, Blackall LL, Pickup RW, et al. Phylogenetic relationships of filamentous sulfur bacteria (Thiothrix spp. and Eikelboom type 021N bacteria) isolated from waste water treatment plants and description of Thiothrix eikelboomii sp. nov., Thiothrix unzii sp. nov., Thiothrix fructosivorans sp. nov. and Thiothrix defluvii sp. nov. Int J Syst Evol Microbiol. 1999;49:1817–27.CAS 

Google Scholar 
Imhoff JF. The family Chromatiaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The Prokaryotes: Gammaproteobacteria. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. p. 151–78.Bižić M, Klintzsch T, Ionescu D, Hindiyeh MY, Günthel M, Muro-Pastor AM, et al. Aquatic and terrestrial cyanobacteria produce methane. Sci Adv. 2020;6:eaax5343.PubMed 
PubMed Central 

Google Scholar 
Rana K, Rana N, Singh B. Chapter 10 – Applications of sulfur oxidizing bacteria. In: Salwan R, Sharma V, editors. Physiological and Biotechnological Aspects of Extremophiles. London, UK: Academic Press; 2020. p. 131–6.
Google Scholar 
Zhuang W-Q, Yi S, Bill M, Brisson VL, Feng X, Men Y, et al. Incomplete Wood-Ljungdahl pathway facilitates one-carbon metabolism in organohalide-respiring Dehalococcoides mccartyi. Proc Natl Acad Sci USA. 2014;111:6419–24.CAS 
PubMed 
PubMed Central 

Google Scholar 
Roncarati D, Scarlato V. Regulation of heat-shock genes in bacteria: from signal sensing to gene expression output. FEMS Microbiol Rev. 2017;41:549–74.CAS 
PubMed 

Google Scholar 
Nagar SD, Aggarwal B, Joon S, Bhatnagar R, Bhatnagar S. A network biology approach to decipher stress response in bacteria using Escherichia coli as a model. OMICS. 2016;20:310–24.CAS 
PubMed 
PubMed Central 

Google Scholar 
Jonas K, Liu J, Chien P, Laub MT. Proteotoxic stress induces a cell-cycle arrest by stimulating lon to degrade the replication initiator DnaA. Cell. 2013;154:623–36.CAS 
PubMed 
PubMed Central 

Google Scholar 
Miss P. Oskarshamn power plant and Clab—Annual report over the radioecological environmental control under 2020. Reg.Nr.2021-02902. Made public 2021-03-21 (In Swedish). Oskarshamn, Sweden; 2021.Lindh MV, Figueroa D, Sjostedt J, Baltar F, Lundin D, Andersson A, et al. Transplant experiments uncover Baltic Sea basin-specific responses in bacterioplankton community composition and metabolic activities. Front Microbiol. 2015;6:223.PubMed 
PubMed Central 

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

Bhatt, S. et al. The global distribution and burden of dengue. Nature 496, 504–509 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gratz, N. G. Critical review of the vector status of Aedes albopictus. Med. Vet. Entomol. 18, 215–227 (2004).CAS 
PubMed 

Google Scholar 
Weaver, S. C., Charlier, C., Nasilakis, N. & Lecuit, M. Zika, chikungunya, and other emerging vector-borne viral diseases. Annu. Rev. Med. 69, 1–14 (2018).
Google Scholar 
Wilder-Smith, A. et al. Epidemic arboviral diseases: priorities for research and public health. Lancet Infect. Dis. 17, e101–e106 (2017).PubMed 

Google Scholar 
Felicetti, T., Manfroni, G., Cecchetti, V. & Cannalire, R. Broad-spectrum flavivirus inhibitors: a medicinal chemistry point of view. Chem. Med. Chem. 15, 2391–2419 (2020).CAS 
PubMed 

Google Scholar 
da Silveira, L. T. C., Bernardo, T. & Santos, M. Systemic review of dengue vaccine efficacy. BMC Inf. Dis. 19, 750 (2019).
Google Scholar 
Katzelnich, L. C. et al. Zika virus infection enhances future risk of severe dengue disease. Science 369, 1123–1128 (2020).ADS 

Google Scholar 
Rezza, G. & Weaver, S. C. Chikungunya as a paradigm for emerging viral diseases: evaluating disease impact and hurdles to vaccine development. PLoS Negl. Trop. Dis. 13, e0006919 (2019).PubMed 
PubMed Central 

Google Scholar 
Reiter, P. & Gubler, D. J. Surveillance and control of urban dengue vectors. In Dengue and dengue hemorrhagic fever (eds Gubler, D. J. & Kuno, G.) 425–462 (CAB International, 1997).
Google Scholar 
Moyes, C. L. et al. Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Negl. Trop. Dis. 11, e0005625 (2017).PubMed 
PubMed Central 

Google Scholar 
Bowman, L. R., Donegan, S. & McCall, P. J. Is dengue vector control deficient in effectiveness or evidence?: systematic review and meta-analysis. PLoS Negl. Trop. Dis. 10, e0004551 (2016).PubMed 
PubMed Central 

Google Scholar 
Erlanger, T. E., Keiser, J. & Utzinger, J. Effect of dengue vector control interventions on entomological parameters in developing countries: a systematic review and meta-analysis. Med. Vet. Entomol. 22, 203–221 (2008).CAS 
PubMed 

Google Scholar 
Banerjee, S., Aditya, G. & Saha, G. K. Household disposables as breeding habitats of dengue vectors: linking wastes and public health. Waste Manag. 33, 233–239 (2013).PubMed 

Google Scholar 
Chadee, D. D., Doon, R. & Severson, D. W. Surveillance of dengue fever cases using a novel Aedes aegypti population sampling method in Trinidad, West Indies: the cardinal points approach. Acta Trop. 104, 1–7 (2007).PubMed 

Google Scholar 
Barrera, R., Acevedo, V. & Amador, M. Role of abandoned and vacant houses on Aedes aegypti productivity. J. Med. Entomol. 104, 145–150 (2020).
Google Scholar 
Chadee, D. D. & Rahaman, A. Use of water drums by humans and Aedes aegypti in Trinidad. J. Vector Ecol. 25, 28–35 (2000).CAS 
PubMed 

Google Scholar 
Padmanabha, H., Soto, E., Mosquera, M., Lord, C. C. & Lounibos, L. P. Ecological links between water storage behaviors and Aedes aegypti production: implications for dengue vector control in variable climates. EcoHealth 7, 78–90 (2010).CAS 
PubMed 

Google Scholar 
Colton, Y. M., Chadee, D. D. & Severson, D. W. Natural skip oviposition of the mosquito Aedes aegypti indicated by codominant genetic markers. Med. Vet. Entomol. 17, 195–204 (2003).CAS 
PubMed 

Google Scholar 
Davis, T. J., Kaufman, P. E., Hogsette, J. A. & Kline, D. I. The effects of larval habitat quality on Aedes albopictus skip oviposition. J. Am. Mosq. Control Assoc. 31, 321–328 (2015).PubMed 

Google Scholar 
David, M. R., Lourenco-de-Oliveira, R. & de Freitas, R. M. Container productivity, daily survival rates and dispersal of Aedes aegypti mosquitoes in a high income dengue epidemic neighbourhood of Rio de Janeiro: presumed influence of differential urban structure on mosquito biology. Mem. Inst. Oswaldo Cruz 104, 927–932 (2009).PubMed 

Google Scholar 
Focks, D. A. & Chadee, D. D. Pupal survey: an epidemiologically significant surveillance method for Aedes aegypti: an example using data from Trinidad. Am. J. Trop. Med. Hyg. 56, 159–167 (1997).CAS 
PubMed 

Google Scholar 
Morrison, A. C. et al. Temporal and geographic patterns of Aedes aegypti (Diptera: Culicidae) production in Iquitos, Peru. J. Med. Entomol. 41, 1123–1142 (2004).PubMed 

Google Scholar 
Chadee, D. D. Oviposition strategies adopted by gravid Aedes aegypti (L.) (Diptera: Culicidae) as detected by ovitraps in Trinidad, West Indies (2002–2006). Acta Trop. 111, 279–283 (2009).CAS 
PubMed 

Google Scholar 
Chadee, D. D. Seasonal incidence and horizontal distribution patterns of oviposition by Aedes aegypti in an urban environment in Trinidad, West Indies. J. Am. Mosq. Control Asso. 8, 281–284 (1992).CAS 

Google Scholar 
Fay, R. W. & Eliason, D. A. A preferred oviposition site as a surveillance method for Aedes aegypti. Mosq. News 26, 531–535 (1966).
Google Scholar 
Johnson, B. J., Ritchie, S. A. & Fonseca, D. M. The state of the art of lethal oviposition trap-based mass interventions for arboviral control. Insects 8, 5 (2017).PubMed Central 

Google Scholar 
Eiras, A. E., Buhagiar, T. S. & Ritchie, S. A. Development of the gravid Aedes trap for the capture of adult female container-exploiting mosquitoes (Diptera: Culicidae). J. Med. Entomol. 51, 200–209 (2014).PubMed 

Google Scholar 
Mackay, A. J., Amador, M. & Barrera, R. An improvied autocidal gravid ovitrap for the control and surveillance of Aedes aegypti. Parasit. Vectors 6, 225 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Olson, K. E. & Blair, C. D. Arbovirus-mosquito interactions: RNAi pathway. Curr. Opin. Virol. 15, 119–126 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hapairai, L. K. et al. Lure-and-kill yeast interfering RNA larvicides targeting neural genes in the human disease vector mosquito Aedes aegypti. Sci. Rep. 7, 13223 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Mysore, K. et al. Yeast interfering RNA larvicides targeting neural genes induce high rates of Anopheles larval mortality. Malaria J. 16, 461 (2017).
Google Scholar 
Mysore, K. et al. Characterization of a broad-based mosquito yeast interfering RNA larvicide with a conserved target site in mosquito semaphorin-1a genes. Parasit. Vectors 12, 256 (2019).PubMed 
PubMed Central 

Google Scholar 
Mysore, K. et al. Characterization of a yeast interfering RNA larvicide with a target site conserved in the synaptotagmin gene of multiple disease vector mosquitoes. PLoS Negl. Trop. Dis 13, e0007422 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hapairai, L. K. et al. Evaluation of large volume yeast interfering RNA lure-and-kill ovitraps for attraction and control of Aedes mosquitoes. Med. Vet. Entomol. 35, 361–370 (2021).CAS 
PubMed 

Google Scholar 
Zeileis, A. & Hothorn, T. Diagnostic checking in regression relationships. R News 2, 7–10 (2002).
Google Scholar 
Braks, M. A. H., Honorio, N. A., Lourenco-de-Oliveira, R., Juliano, S. A. & Lounibos, L. P. Convergent habitat segregation of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in southeastern Brazil and Florida. J. Med. Entomol. 40, 785–794 (2003).PubMed 

Google Scholar 
Kumari, R., Kumar, K. & Chauhan, L. S. First dengue virus detection in Aedes albopictus from Delhi, India: Its breeding ecology and role in dengue transmission. Trop. Med. Int. Health 16, 949–954 (2012).
Google Scholar 
Apostol, B. L., Black, W. C. IV., Reiter, P. & Miller, B. R. Use of randomly amplified polymorphic DNA amplified by polymerase chain reaction markers to estimate the number of Aedes aegypti families at oviposition sites in San Juan, Puerto Rico. Am. J. Trop. Med. Hyg. 51, 89–97 (1994).CAS 
PubMed 

Google Scholar 
Corbet, P. S. & Chadee, D. D. An improved method for detecting substrate preferences shown by mosquitoes that exhibit ‘skip oviposition’. Physiol. Entomol. 18, 114–118 (1993).
Google Scholar 
Reinbold-Wasson, D. D. & Reiskind, M. H. Comparative skip-oviposition behavior among container breeding Aedes spp. mosquitoes (Diptera: Culicidae). J. Med. Entomol. https://doi.org/10.1093/jme/tjab084 (2021).Article 
PubMed 

Google Scholar 
Barrera, R. Spatial stability of adult Aedes aegypti populations. Am. J. Trop. Med. Hyg. 85, 1087–1092 (2011).PubMed 
PubMed Central 

Google Scholar 
Barrera, R., Amador, M., Ruiz-Valcarcel, J. & Acevedo, V. Factors modulating captures of gravid Aedes aegypti females. J. Am. Mosq. Control Assoc. 36, 66–73 (2020).PubMed 
PubMed Central 

Google Scholar 
Moura, M. B. C. M. et al. Spatio-temporal dynamics of Aedes aegypti and Aedes albopictus oviposition in an urban area of northeastern Brazil. Trop. Med. Int. Health 25, 1510–1521 (2020).PubMed 

Google Scholar 
Crawford, J. E. et al. Efficient production of male Wolbachia-infected Aedes aegypti mosquitoes enables large-scale suppression of wild populations. Nat. Biotechnol. 38, 482–492 (2020).CAS 

Google Scholar 
Lau, K. W. et al. Vertical distribution of Aedes mosquitoes in multiple story buildings in Selangor and Kuala Lumpur, Malaysia. Trop. Biomed. 30, 36–45 (2013).CAS 
PubMed 

Google Scholar 
Perich, M. J. et al. Field evaluation of a lethal ovitrap against dengue vectors in Brazil. Med. Vet. Entomol. 17, 205–210 (2003).CAS 
PubMed 

Google Scholar 
Serpa, L. L. N. et al. Study of the the distribution and abundance of the eggs of Aedes aegypti and Aedes albopictus according to the habitat and meteorlogical variables, municipality of Sao Sebastiao, Sao Paulo state, Brazil. Parasit. Vectors 6, 321 (2014).
Google Scholar 
Sithiprasasna, R. et al. Field evaluation of a lethal ovitrap for the control of Aedes aegypti (Diptera: Culicidae) in Thailand. J. Med. Entomol. 40, 455–462 (2003).PubMed 

Google Scholar 
Barrera, R., Amador, M., Munoz, J. & Acevedo, V. Integrated vector control of Aedes aegypti mosquitoes around target houses. Parasit. Vectors 11, 88 (2018).PubMed 
PubMed Central 

Google Scholar 
Naranjo, D. P. et al. Vector control programs in Saint Johns County, Florida and Guayas, Ecuador: Successes and barriers to integrated vector management. BMC Public Health 14, 674 (2014).PubMed 
PubMed Central 

Google Scholar 
Regis, L. N. et al. Sustained reduction of the dengue vector population resulting from an integrated control strategy applied in two Brazilian cities. PLoS ONE 8, e67682 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Stewart, A. T. M. et al. Community acceptance of yeast interfering RNA larvicide technology for control of Aedes mosquitoes in Trinidad. PLoS ONE 15, e0237675 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Winter, N. et al. Assessment of Trinidad community stakeholder perspectives on the use of yeast interfering RNA-baited ovitraps for biorational control of Aedes mosquitoes. PLoS ONE 16, e0252997 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Chadee, D. D. & Corbet, P. S. Seasonal incidence and diel patterns of oviposition in the field of the mosquito, Aedes aegypti (L.) (Diptera: Culicidae) in Trinidad, West Indies: a preliminary study. Ann. Trop. Med. Parasitol. 81, 151–161 (1987).CAS 
PubMed 

Google Scholar 
Edman, J. D. et al. Aedes aegypti (Diptera: Culicdae) movement influenced by availability of oviposition sites. J. Med. Entomol. 35, 578–583 (1998).CAS 
PubMed 

Google Scholar 
Reiter, P., Amador, M. A., Anderson, R. A. & Clark, G. G. Short report: dispersal of Aedes aegypti in an urban area after blood feeding as demonstrated by rubidium-marked eggs. Am. J. Trop. Med. Hyg. 52, 177–179 (1995).CAS 
PubMed 

Google Scholar 
Mysore, K. et al. Preparation and use of a yeast shRNA delivery system for gene silencing in mosquito larvae. Methods Mol. Biol. 1858, 213–231 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Chadee, D. D., Fat, F. H. & Persad, R. C. First record of Aedes albopictus from Trinidad, West Indies. J. Am. Mosq. Control Assoc. 19, 438–439 (2003).PubMed 

Google Scholar 
Clemons, A., Mori, A., Haugan, M., Severson, D. W. & Duman-Scheel, M. Culturing and egg collection of Aedes aegypti. Cold Spring Harb. Protoc. https://doi.org/10.1101/pdb.prot5507 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Stan Developmental Team. Stan Modeling Language Users Guide and Reference Manual, v.2.22.1 https://mc-stan.org (2020). More

Holmgren, M. et al. Extreme climatic events shape arid and semiarid ecosystems. Front. Ecol. Environ. 4, 87–95 (2006).
Google Scholar 
Ummenhofer, C. C. & Meehl, G. A. Extreme weather and climate events with ecological relevance: a review. Philos. Trans. R. Soc. B-Biol. Sci. 372, 20160135. https://doi.org/10.1098/rstb.2016.0135 (2017).Chesson, P. et al. Resource pulses, species interactions, and diversity maintenance in arid and semi-arid environments. Oecologia 141, 236–253 (2004).ADS 
PubMed 

Google Scholar 
McCluney, K. E. et al. Shifting species interactions in terrestrial dryland ecosystems under altered water availability and climate change. Biol. Rev. 87, 563–582 (2012).PubMed 

Google Scholar 
Reyer, C. P. O. et al. A plant’s perspective of extremes: Terrestrial plant responses to changing climatic variability. Glob. Change Biol. 19, 75–89 (2013).ADS 

Google Scholar 
Scheffer, M., Carpenter, S., Foley, J. A., Folke, C. & Walker, B. Catastrophic shifts in ecosystems. Nature 413, 591–596 (2001).ADS 
CAS 
PubMed 

Google Scholar 
Schwinning, S. & Sala, O. E. Hierarchy of responses to resource pulses in and and semi-arid ecosystems. Oecologia 141, 211–220 (2004).ADS 
PubMed 

Google Scholar 
Borer, E. T., Seabloom, E. W. & Tilman, D. Plant diversity controls arthropod biomass and temporal stability. Ecol. Lett. 15, 1457–1464 (2012).PubMed 

Google Scholar 
Kwok, A. B. C., Wardle, G. M., Greenville, A. C. & Dickman, C. R. Long-term patterns of invertebrate abundance and relationships to environmental factors in arid Australia. Austral Ecol. 41, 480–491 (2016).
Google Scholar 
Prugh, L. R. et al. Ecological winners and losers of extreme drought in California. Nat. Climate Change 8, 819–824 (2018).ADS 

Google Scholar 
Deguines, N., Brashares, J. S. & Prugh, L. R. Precipitation alters interactions in a grassland ecological community. J. Anim. Ecol. 86, 262–272 (2017).PubMed 

Google Scholar 
Ripple, W. J. et al. What is a trophic cascade?. Trends Ecol. Evol. 31, 842–849 (2016).PubMed 

Google Scholar 
Greenville, A. C., Wardle, G. M. & Dickman, C. R. Extreme climatic events drive mammal irruptions: regression analysis of 100-year trends in desert rainfall and temperature. Ecol. Evol. 2, 2645–2658 (2012).PubMed 
PubMed Central 

Google Scholar 
Molyneux, J., Pavey, C. R., James, A. I. & Carthew, S. M. Persistence of ground-dwelling invertebrates in desert grasslands during a period of low rainfall—Part 2. J. Arid. Environ. 157, 39–47 (2018).ADS 

Google Scholar 
Seymour, C. L., Simmons, R. E., Joseph, G. S. & Slingsby, J. A. On bird functional diversity: Species richness and functional differentiation show contrasting responses to rainfall and vegetation structure in an arid landscape. Ecosystems 18, 971–984 (2015).
Google Scholar 
Prather, C. M. et al. Invertebrates, ecosystem services and climate change. Biol. Rev. 88, 327–348 (2013).PubMed 

Google Scholar 
Del Toro, I., Ribbons, R. R. & Pelini, S. L. The little things that run the world revisited: a review of ant-mediated ecosystem services and disservices (Hymenoptera: Formicidae). Myrmecol. News 17, 133–146 (2012).
Google Scholar 
Gerlach, J., Samways, M. & Pryke, J. Terrestrial invertebrates as bioindicators: an overview of available taxonomic groups. J. Insect Conserv. 17, 831–850 (2013).
Google Scholar 
Doblas-Miranda, E., Sanchez-Pinero, F. & Gonzalez-Megias, A. Different microhabitats affect soil macroinvertebrate assemblages in a Mediterranean arid ecosystem. Appl. Soil Ecol. 41, 329–335 (2009).
Google Scholar 
Hadley, N. F. & Szarek, S. R. Productivity of desert ecosystems. Bioscience 31, 747–753 (1981).
Google Scholar 
Barnett, K. L. & Facey, S. L. Grasslands, invertebrates, and precipitation: A review of the effects of climate change. Front. Plant Sci. 7, 1196 (2016).PubMed 
PubMed Central 

Google Scholar 
Zhu, H. et al. Effects of altered precipitation on insect community composition and structure in a meadow steppe. Ecol. Entomol. 39, 453–461 (2014).
Google Scholar 
Palmer, C. M. Chronological changes in terrestrial insect assemblages in the arid zone of Australia. Environ. Entomol. 39, 1775–1787 (2010).PubMed 

Google Scholar 
Liu, R. T., Zhu, F. & Steinberger, Y. Ground-active arthropod responses to rainfall-induced dune microhabitats in a desertified steppe ecosystem, China. J. Arid Land 8, 632–646 (2016).
Google Scholar 
Mendelsohn, J., Jarvis, A., Roberts, C. & Robertson, T. Atlas of Namibia: A portrait of the land and its people. 3rd edn, (Sunbird Publishers, 2009).Theron, L. Temporal and spatial composition of arboreal insects along the Omaruru river, Namibia Magister scientiae thesis, University of the Free State Bloemfontein, (2010).Wagner, T. C., Richter, J., Joubert, D. F. & Fischer, C. A dominance shift in arid savanna: An herbaceous legume outcompetes local C4 grasses. Ecol. Evol. 8, 6779–6787 (2018).PubMed 
PubMed Central 

Google Scholar 
Wagner, T. C., Hane, S., Joubert, D. F. & Fischer, C. Herbaceous legume encroachment reduces grass productivity and density in arid rangelands. PLoS ONE 11, e0166743; https://doi.org/10.1371/journal.pone.0166743 (2016).Picker, M., Griffiths, C. & Weaving, A. Field Guide to Insects of Southern Africa. (Struik Nature, 2004).Scholtz, C. H. & Holm, E. Insects of Southern Africa. 2nd edn, (Protea Book House, 2008).Blaum, N., Seymour, C., Rossmanith, E., Schwager, M. & Jeltsch, F. Changes in arthropod diversity along a land use driven gradient of shrub cover in savanna rangelands: identification of suitable indicators. Biodivers. Conserv. 18, 1187–1199 (2009).
Google Scholar 
Franca, L. F., Figueiredo-Paixao, V. H., Duarte-Silva, T. A. & dos Santos, K. B. The effects of rainfall and arthropod abundance on breeding season of insectivorous birds, in a semi-arid neotropical environment. Zoologia-Curitiba. https://doi.org/10.3897/zoologia.37.e37716 (2020).Wagner, T. C., Uiseb, K. & Fischer, C. Rolling pits of Hartmann’s mountain zebra (Zebra equus hartmannae) increase vegetation diversity and landscape heterogeneity in the Pre-Namib. Ecol. Evol. 11, 13036–13051 (2021).PubMed 
PubMed Central 

Google Scholar 
Dormann, C. F. et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36, 27–46 (2013).
Google Scholar 
R Core Team. R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria, 2019).Oksanen, J., et al. vegan: Community Ecology Package. R package version 2.5-7. (2020).Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280 (2001).ADS 
PubMed 
PubMed Central 

Google Scholar 
Anderson, M. J. & Walsh, D. C. I. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing?. Ecol. Monogr. 83, 557–574 (2013).
Google Scholar 
Anderson, M. J. in Wiley StatsRef: Statistics Reference Online (eds N. Balakrishnan et al.) (2017).Stopher, K. V., Bento, A. I., Clutton-Brock, T. H., Pemberton, J. M. & Kruuk, L. E. B. Multiple pathways mediate the effects of climate change on maternal reproductive traits in a red deer population. Ecology 95, 3124–3138 (2014).
Google Scholar 
Bolker, B. M. et al. Generalized linear mixed models: A practical guide for ecology and evolution. Trends Ecol. Evol. 24, 127–135 (2009).PubMed 

Google Scholar 
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Google Scholar 
Pinheiro, J. C. & Bates, D. M. Mixed-Effects Models in S and S-PLUS. (Springer Verlag, 2000).Zhang, D. rsq: R-Squared and related measures. R package version 2.2. (2021).Barnes, A. D. et al. Direct and cascading impacts of tropical land-use change on multi-trophic biodiversity. Nat. Ecol. Evol. 1, 1511–1519 (2017).PubMed 

Google Scholar 
Henschel, J. R. Long-term population dynamics of Namib desert Tenebrionid beetles reveal complex relationships to pulse-reserve conditions. Insects 12, 804. https://doi.org/10.3390/insects12090804 (2021).Cloudsley-Thompson, J. L. The adaptational diversity of desert biota. Environ. Conserv. 20, 227–231 (1993).
Google Scholar 
Sømme, L. in Invertebrates in Hot and Cold Arid Environments 135–157 (Springer, 1995).Suttle, K. B., Thomsen, M. A. & Power, M. E. Species interactions reverse grassland responses to changing climate. Science 315, 640–642 (2007).ADS 
CAS 
PubMed 

Google Scholar 
Henschel, J., Klintenberg, P., Roberts, C. & Seely, M. Long-term ecological research from an arid, variable, drought-prone environment. Sécheresse 18, 342–347 (2007).
Google Scholar 
Cloudsley-Thompson, J. L. Adaptations of arthropoda to arid environments. Annu. Rev. Entomol. 20, 261–283 (1975).CAS 
PubMed 

Google Scholar 
Schuldt, A. et al. Belowground top-down and aboveground bottom-up effects structure multitrophic community relationships in a biodiverse forest. Sci. Rep. 7 (2017).Vidal, M. C. & Murphy, S. M. Bottom-up vs. top-down effects on terrestrial insect herbivores: a meta-analysis. Ecol. Lett. 21, 138–150 (2018).Báez, S., Collins, S. L., Lightfoot, D. & Koontz, T. L. Bottom-up regulation of plant community structure in an aridland ecosystem. Ecology 87, 2746–2754 (2006).PubMed 

Google Scholar 
Gibb, H. et al. Testing top-down and bottom-up effects on arid zone beetle assemblages following mammal reintroduction. Austral Ecol. 43, 288–300 (2018).
Google Scholar 
Coll, M. & Guershon, M. Omnivory in terrestrial arthropods: Mixing plant and prey diets. Annu. Rev. Entomol. 47, 267–297 (2002).CAS 
PubMed 

Google Scholar 
Karolyi, F., Hansal, T., Krenn, H. W. & Colville, J. F. Comparative morphology of the mouthparts of the megadiverse South African monkey beetles (Scarabaeidae: Hopliini): feeding adaptations and guild structure. PeerJ 4, e1597; https://doi.org/10.7717/peerj.1597 (2016).Greenslade, P. Survival of Collembola in arid environments: Observations in South Australia and the Sudan. J. Arid. Environ. 4, 219–228 (1981).ADS 

Google Scholar 
Fattorini, S. Effects of fire on tenebrionid communities of a Pinus pinea plantation: A case study in a Mediterranean site. Biodivers. Conserv. 19, 1237–1250 (2009).
Google Scholar 
Sanders, N. J., Moss, J. & Wagner, D. Patterns of ant species richness along elevational gradients in an arid ecosystem. Glob. Ecol. Biogeogr. 12, 93–102 (2003).
Google Scholar  More

Design parametersA new type of Debris-flow grille dam is proposed to be built with a height of 8 m. Column section 500 mm × 700 mm, spacing 5000 mm. The cross section of the beam is 400 mm × 300 mm, and the spacing is 4000 mm. The section steel adopts I-steel 45a, the spacing is 250 mm. The counterfort wall is 300 mm thick and 6500 mm high. Pile foundation adopts manual digging pile, pile by 1000 mm, 5000 mm deep. The concrete is C30; Stressed bar is HRB335; Stirrups is HRB300; Stay Cable is 3 (emptyset) s15.2. The design size of the anchor piers is shown in Fig. 12. In the Figure where (T = 2 times 10^{5} N); (L_{l} = 8500;{text{mm}}); (E_{l} = 1.95 times 10^{5} ;{text{N/mm}}^{2}); (A_{l} = 420;{text{mm}}); (D_{e} = 1000;{text{mm}}); (L_{m} = 1200;{text{mm}}); (E_{e} = 3.0 times 10^{4} ;{text{N/mm}}^{2}); (H = 1000;{text{mm}}); (mu = 0.2); (E = 20;{text{N/mm}}^{2}). The parameter of gully bed soil is shown in Table 1.Figure 12The parameters of anchor piers.Full size imageTable 1 The parameters of gully bed soil.Full size tableAnalysis of results(1) The effect of the elastic modulus and Poisson’s ratio of the surrounding soil on the displacement deformation of the anchor-pulling system.The elastic modulus (E) and Poisson’s ratio (mu) are important parameters for calculating the displacement deformation of soil. They have something to do with both the properties of materials and the stress level. To analyze the effect of the physical parameter variation of the surrounding soil on the displacement deformation of the anchor-pulling system, we can study changing the elastic modulus and Poisson’s ratio. The variation range of the elastic modulus is 15–45 N/mm2, and the variation range of Poisson’s ratio is 0.15–0.25.Figure 13 shows the variation curve in which the displacement deformation increases with the elastic modulus of the soil around the anchor pier. We can see that as the elastic modulus of the soil around the anchor pier increases, the displacement deformation decreases gradually. When the elastic modulus is in the range of 15–35 N/mm2, the curve is steep, and the decrease in deformation is apparent. After 35 N/mm2, the curve becomes smooth, and the decrease in deformation tends to be stable.Figure 13The effect of the elastic modulus E(15–45 N/mm2) of the surrounding soil on the displacement of the anchor-pulling system.Full size imageIn Fig. 14, the displacement deformation increases linearly with Poisson’s ratio of the soil around the anchor pier. However, the total impact is not large. From calculation, the variation of elastic modulus of the soil around the anchor pier has nothing to do with elastic deformation of the stayed cable ((S_{1} )), but mainly influences relative shear displacement between anchor piers and the surrounding soil ((S_{2} )) and the compression performance of the soil on the front of anchor piers ( (S_{3} )). where ((S_{2} )) accounted for 89% and (left( {S_{3} } right)) accounted for 11%. When the Poisson ratio increases, the displacement deformation also increases. Poisson’s ratio has the greatest influence on the relative shear displacement ((S_{2} )) of the anchor pier and soil, accounting for approximately 96.4%. The design parameters should be selected correctly during design. The influence of parameters on the deformation of anchor system is analyzed by using control variable method. The influence of a single variable on the results can be intuitively obtained. However, the elastic modulus E and Poisson ‘ s ratio (mu) of rock and soil are not independent. Therefore, Matlab is used to analyze the influence of the two aspects on the deformation of the tensile anchor system, and the results are shown in Fig. 15. It can be seen from Fig. 15 that the influence of elastic modulus E on the deformation of tensile anchor system is much greater than that of Poisson’s ratio (mu). And the variation of the curve is basically the same, so the interaction between the two is weak.Figure 14The effect of Poisson’s ratio (mu)(0.15–0.26) of the surrounding soil on the displacement of the anchor-pulling system.Full size imageFigure 15Influence of elastic modulus E (15–45 N/mm2) and Poisson’s ratio (mu left( {0.15 – 0.26} right)) on deformation of anchor system.Full size image(2) The effect of the design parameters of anchor piers on the displacement deformation of the anchor-pulling system.The design parameters of anchor piers include the equivalent width (D_{e}), length (L_{m}) and height (H). Different design parameters have varying effects on the displacement deformation of the anchor-pulling system. Keep other parameters unchanged and let ( D_{e} ) vary in 0.5–1.5 m, (L_{m}) vary in 0.6–2.0 m, and (H) vary in 0.5–1.5 m. Analyzing their effect on the displacement deformation of the anchor-pulling system, the results are shown in Figs. 16 and 17.Figure 16The effect of equivalent width (D_{e})(500–1500 mm) on the displacement of the anchor-pulling system.Full size imageFigure 17The effect of equivalent length (L_{m})(600–2000 mm) on the displacement of the anchor-pulling system.Full size imageAs illustrated in Figs. 16 and 17, the effects of the design parameters of the anchor piers on the displacement deformation of the anchor-pulling system are almost the same. As the size increases, the displacement deformation gradually decreases, and the front section decreases quickly, while the rear section becomes gradually smooth. Here, the equivalent width (D_{e}) and length (L_{m}) mainly affect the compression performance of the soil on the front of anchor piers (left( {S_{3} } right)). The anchor piers can be seen as rigid bodies where horizontal displacement takes place. Increasing the size means increasing the contact area between the anchor pier and soil body. With this increase, the compression performance of the soil on the front of the anchor piers decreases. However, the effect of the height (H) on the displacement deformation of the anchor-pulling system is the contribution to the relative shear displacement between the anchor piers and the surrounding soil ((S_{2} )). When (H) grows, ((S_{2} )) grows accordingly. However, theoretically, the larger the effect of the size, the better it is. Because of the constraint of topographic conditions, construction conditions and economic benefits in practical engineering, it is necessary to choose the best size. the anchor pier provides enough anchor force and saves all kinds of resources. The best design dimensions suggested are (D_{e}) = 1.2 m–1.8 m, (L_{m}) = 1.5 m–2.5 m, and (H) = 1.0 m–1.6 m.It can be seen from Fig. 18 that the width (D_{e}) and the height (L_{m}) of anchor pier influence each other greatly. When (D_{e}) is 600 mm, with the increase of (L_{m}), the deformation of tension anchor system will first decrease and then increase. When (D_{e}) is greater than 800 mm, with the increase of (L_{m}), the deformation of tension anchor system will continue to decrease. And with the increase of (L_{m}), the decreasing trend is more obvious. When (L_{m}) is 500 mm, with the increase of the height of the anchor pier (D_{e}), the deformation of the anchor system will increase first. When (L_{m}) is greater than 800 mm, with the increase of (D_{e}), the deformation of the anchor system will continue to decrease. But the decreasing trend is not much different.Figure 18Influence of Anchor Pier Width (D_{e} left( {500 – 1500;{text{mm}}} right)) and Anchor Pier Height (L_{m} left( {600 – 2000;{text{mm}}} right)) on Deformation of Anchorage System.Full size imageThe numerical validationThe establishment of the finite element modelWhen the finite element model of the anchor-pulling system and surrounding soil is created, the constitutive model of the surrounding soil uses the Mohr–Coulomb elastoplastic model. The anchor pier and surrounding soil use eight nodes as oparametric elements, such as solid45, of which the basic grid unit is cubic units. When the grid is divided, the grid between the anchor pier and the surrounding soil contact is dense. The LINK10 unit is used to simulate cables, which have a bilinear stiffness matrix. It can simulate not only tensile bar units but also compressed bar units. For example, when the pull-up option is used alone, if the unit is under pressure, its stiffness disappears, so it can be used to simulate the relaxation of cables or chains. This feature is very significant for the static problem of wire rope, which uses a unit to simulate the entire cable. It can also be used for dynamic analysis with inertial or damping effects when the needed relaxation unit should pay attention to its performance rather than its movement. The soil is homogeneous. The soil physical parameters and structure design parameters are consistent with the theoretical calculation parameters mentioned above. The tensile force of the cable is exerted on the nodes as a force. The top surface of the model is free, and the normal displacements of the remaining faces are constrained such that the displacements are zero. The contact of the anchor pier and surrounding soils is a rigid-flexible surface-to-surface contact element to reflect the interaction. The surface of the anchor pier is regarded as the “target” surface, and the surface of the soil body is regarded as the “contact” surface. The coefficient of friction and normal penalty stiffness are 0.35 and 0.15, respectively. The scope of interaction between the anchor pier and the surrounding soil in the model is taken as 15 m × 11 m × 12 m, referring to past experience in engineering and the research data of the effect scope that the related anchors have had on the soil. The values of the model geometric parameters and physical and mechanical parameters are the same as in “Design parameters” section. The finite element model is shown in Fig. 19.Figure 19Finite element model of the anchor-pulling system and surrounding soil.Full size imageResearch on finite element model gridIn order to verify the convergence of numerical simulation, the soil was divided into three different mesh sizes. Condition 1 is fine finite element meshing. The stress nephogram of condition 1 is shown in Fig. 20. Condition 2 is medium finite element mesh. The stress nephogram of condition 1 is shown in Fig. 21. Condition 3 is coarse finite element mesh. The stress nephogram of condition 1 is shown in Fig. 22. See Table 2 for specific grid division.Figure 20Condition 1 stress cloud diagram.Full size imageFigure 21Condition 1 stress cloud diagram.Full size imageFigure 22Condition 1 stress cloud diagram.Full size imageTable 2 Mesh size of three working conditions.Full size tableIt can be seen from the stress nephogram of the three working conditions that the thicker the grid is, the greater the displacement of the anchor system is. The maximum displacement difference between condition 2 and condition 3 is 2.6%; the maximum displacement of condition 1 is 17% different from that of condition 2. The finer the mesh, the more accurate the numerical simulation results. But with the increase in computing time. It can be seen from Table 2 that the maximum iteration of condition 1 is 10 times, and the result will converge. The maximum iterations of condition 2 and 3 only need 7 times, and the results can converge.The calculation resultsFigure 23 and Fig. 24 are the displacement nephograms of the soil around the anchor piers for 100 kN and 400 kN, respectively. The soil displacement increases with increasing load, the affected area will increase and become uniform, and the area under load will also increase. The soil within the range of 1–3 m around the anchor pier is greatly affected, accounting for 80% of the total force. The soil around the anchor pier should be reinforced, and the anchoring force should be enhanced in the design.Figure 23Displacement fringe of soil around the anchor piers for 100 kN.Full size imageFigure 24Displacement fringe of soil around the anchor piers for 400 kN.Full size imageIn order to further study the influence of anchorage pier size on the displacement and deformation of anchorage system, finite element models with different sizes are established by finite element method. The stress nephogram is shown in Figs. 25, 26 and 27.Figure 25Top 800 mm, bottom 800 mm anchor pier stress nephogram.Full size imageFigure 26Top 1000 mm, bottom 1000 mm anchor pier stress nephogram.Full size imageFigure 27Top 800 mm, bottom 1000 mm anchor pier stress nephogram.Full size imageFrom Figs. 25, 26 and 27, it can be seen that when the anchor pier is rectangular, the deformation of the tensile anchor system decreases with the increase of the size of the anchor pier, but the degree is small. When the anchor pier is trapezoidal, the material is small, but the deformation is more ideal than the rectangular. It can be seen that reasonable selection of anchor pier size is crucial, not blindly increase the size of anchor pier.Figure 28 shows that the displacement of the soil around the anchor pier increases with increasing load, and the added value is obvious at approximately 2–3 mm. Figure 29 shows that the increase in load has a great effect on the soil in front of the anchor pier. As the load increases, the compressive deformation of the soil gradually increases. As the distance from the anchor pier increases, the displacement of the soil decreases, and the scope of influence gradually decreases. The displacement of the soil tends to be stable beyond 4–5 m from the anchor pier.Figure 28The displacement of soil around anchor pier.Full size imageFigure 29The horizontal displacement of soil along cable axis.Full size imageComparison of theoretical calculation and numerical simulation results at the time of load variationTo verify the correctness of the theoretical calculation, we compare the theoretical calculation with numerical simulation results of displacement deformation of anchor-pulling system under different pulling force of stayed cable. The results are shown in Fig. 30, see Table 3 for data.Figure 30Comparison of theoretical calculation and numerical simulation results.Full size imageTable 3 Comparison between theoretical calculation and numerical simulation.Full size tableAs seen from Fig. 30, the theoretical and numerical simulation results are consistent, showing a linear growth trend. The slope difference of the two straight lines is approximately 5%, which meets the accuracy requirements of geotechnical engineering. As the restraint effect of the surrounding soil on the anchor pier is not fully considered, the theoretical calculation result is too large. The deformation of anchor (left( {S_{1} } right)) in displacement deformation is the same, and the relative shear displacement (left( {S_{2} } right)) of the anchor pier and the soil and the compressive deformation ((S_{3} )) of the soil at the front end of the anchor pier are 1.25 times and 1.08 times the numerical simulation results, respectively. The change in (left( {S_{2} } right)) in the calculation results is large and should be taken into account in the design. More