Philippa Kaur

Phung, T. X., Nascimento, J. C. S., Novarro, A. J. & Wiens, J. J. Correlated and decoupled evolution of adult and larval body size in frogs. Proc. R. Soc. Lond. B 287, 20201474–10 (2020).
Google Scholar 
Hall, B. K. & Wake, M. H. The Origin and Evolution of Larval Forms (Gulf Professional Publishing, 1999).Hime, P. M. et al. Phylogenomics reveals ancient gene tree discordance in the amphibian tree of life. Syst. Biol. 70, 49–66 (2021).CAS 
PubMed 

Google Scholar 
Duellman, W. E. & Trueb, L. Biology of Amphibians (Johns Hopkins University Press, 1994).Nunes-de-Almeida, C. H., Batista Haddad, C. F. & Toledo, L. F. A revised classification of the amphibian reproductive modes. Salamandra 57, 413–427 (2021).
Google Scholar 
Vences, M. & Köhler, J. Global diversity of amphibians (Amphibia) in freshwater. Hydrobiologia 595, 569–580 (2007).
Google Scholar 
Blackburn, D. G. Evolution of vertebrate viviparity and specializations for fetal nutrition: a quantitative and qualitative analysis. J. Morphol. 276, 961–990 (2015).PubMed 

Google Scholar 
AmphibiaWeb. Electronic Database (University of California, Berkeley, CA, USA, 2019). https://amphibiaweb.org.Frost, D. R. Amphibian Species of the World: an Online Reference. Version 6.0 (date of access: 01.08.2019). Electronic Database (American Museum of Natural History, New York, USA, 2019). https://amphibiansoftheworld.amnh.org/.Bonett, R. M., Ledbetter, N. M., Hess, A. J., Herrboldt, M. A. & Denoël, M. Repeated ecological and life cycle transitions make salamanders an ideal model for evolution and development. Developmental Dynamics 251, 957–972 (2022).Salthe, S. N. Reproductive modes and the number and sizes of ova in the urodeles. Am. Midl. Naturalist 81, 467490 (1969).
Google Scholar 
Salthe, S. N. & Duellman, W. E. in Evolutionary Biology of the Anurans (ed. Vial, J. L.) 229–249 (University of Missouri Press Columbia, 1973).Haddad, C. & Prado, C. P. A. Reproductive modes in frogs and their unexpected diversity in the Atlantic Forest of Brazil. BioScience 55, 207–217 (2005).
Google Scholar 
Lutz, B. Ontogenetic evolution in frogs. Evolution 2, 29–39 (1948).CAS 
PubMed 

Google Scholar 
Crump, M. L. Anuran reproductive modes: evolving perspectives. J. Herpetol. 49, 1–16 (2015).
Google Scholar 
Schoch, R. Evolution of life cycles in early amphibians. Annu. Rev. Earth Planet. Sci. 37, 135–162 (2009).ADS 
CAS 

Google Scholar 
Meegaskumbura, M. et al. Patterns of reproductive-mode evolution in Old World tree frogs (Anura, Rhacophoridae). Zool. Scr. 44, 509–522 (2015).
Google Scholar 
Portik, D. M. & Blackburn, D. C. The evolution of reproductive diversity in Afrobatrachia: A phylogenetic comparative analysis of an extensive radiation of African frogs. Evolution 70, 2017–2032 (2016).PubMed 
PubMed Central 

Google Scholar 
San Mauro, D. et al. Life-history evolution and mitogenomic phylogeny of caecilian amphibians. Mol. Phylogenet. Evol. 73, 177–189 (2014).PubMed 

Google Scholar 
Pereira, E. B., Collevatti, R. G., de Carvalho Kokubum, M. N., de Oliveira Miranda, N. E. & Maciel, N. M. Ancestral reconstruction of reproductive traits shows no tendency toward terrestriality in leptodactyline frogs. BMC Evol. Biol. 15, 91 (2015).Gomez-Mestre, I., Pyron, R. A. & Wiens, J. J. Phylogenetic analyses reveal unexpected patterns in the evolution of reproductive modes in frogs. Evolution 66, 3687–3700 (2012).PubMed 

Google Scholar 
Wake, D. B. & Hanken, J. Direct development in the lungless salamanders: what are the consequences for developmental biology, evolution and phylogenesis? Int. J. Dev. Biol. 40, 859–869 (1996).CAS 
PubMed 

Google Scholar 
Dubois, A. Developmental pathway, speciation and supraspecific taxonomy in amphibians: 1. Why are there so many frog species in Sri Lanka? Alytes 22, 19–37 (2004).
Google Scholar 
Hedges, S. B., Duellman, W. E. & Heinicke, M. P. New World direct-developing frogs (Anura: Terrarana): molecular phylogeny, classification, biogeography, and conservation. Zootaxa 1737, 1–182 (2008).
Google Scholar 
Dugo-Cota, Á., Vilà, C., Rodríguez, A. & Gonzalez-Voyer, A. Ecomorphological convergence in Eleutherodactylus frogs: a case of replicate radiations in the Caribbean. Ecol. Lett. 22, 884–893 (2019).PubMed 

Google Scholar 
Simpson, G. G. The Major Features of Evolution (Columbia University Press, 1953).Vági, B., Végvári, Z., Liker, A., Freckleton, R. P. & Szekely, T. Parental care and the evolution of terrestriality in frogs. Proc. R. Soc. Lond. B 286, 20182737–10 (2019).
Google Scholar 
Furness, A. I. & Capellini, I. The evolution of parental care diversity in amphibians. Nat. Commun. 10, 1–12 (2019).CAS 

Google Scholar 
Furness, A. I., Venditti, C. & Capellini, I. Terrestrial reproduction and parental care drive rapid evolution in the trade-off between offspring size and number across amphibians. PLoS Biol. 20, e3001495 (2022).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wollenberg, K. C., Vieites, D. R., Glaw, F. & Vences, M. Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs. BMC Evol. Biol. 11, 217 (2011).PubMed 
PubMed Central 

Google Scholar 
Cayuela, H. et al. Determinants and consequences of dispersal in vertebrates with complex life cycles: a review of pond-breeding amphibians. Q. Rev. Biol. 95, 1–36 (2020).
Google Scholar 
Phillimore, A. B., Freckleton, R. P., Orme, C. D. L. & Owens, I. P. F. Ecology predicts large‐scale patterns of phylogenetic diversification in birds. Am. Nat. 168, 220–229 (2006).PubMed 

Google Scholar 
Chen, J.-M. et al. An integrative phylogenomic approach illuminates the evolutionary history of Old World tree frogs (Anura: Rhacophoridae). Mol. Phylogenet. Evol. 145, 106724 (2020).PubMed 

Google Scholar 
Zimkus, B. M., Lawson, L., Loader, S. P. & Hanken, J. Terrestrialization, miniaturization and rates of diversification in African Puddle Frogs (Anura: Phrynobatrachidae). PLoS ONE 7, e35118 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Moen, D. S. Improving inference and avoiding over-interpretation of hidden-state diversification models: specialized plant breeding has no effect on diversification in frogs. Evolution 76, 373–384 (2022).PubMed 

Google Scholar 
Eastman, J. M. & Storfer, A. Correlations of life-history and distributional-range variation with salamander diversification rates: evidence for species selection. Syst. Biol. 60, 503–518 (2011).PubMed 

Google Scholar 
Bonett, R. M., Steffen, M. A., Lambert, S. M., Wiens, J. J. & Chippindale, P. T. Evolution of paedomorphosis in plethodontid salamanders: ecological correlates and re-evolution of metamorphosis. Evolution 68, 466–482 (2014).PubMed 

Google Scholar 
Akaike, H. A new look at the statistical-model identification. IEEE Trans. Autom. Control 19, 716–723 (1974).ADS 
MathSciNet 
MATH 

Google Scholar 
Wagenmakers, E.-J. & Farrell, S. AIC model selection using Akaike weights. Psychon. Bull. Rev. 11, 192–196 (2004).PubMed 

Google Scholar 
Beaulieu, J. M., Oliver, J. C., O’Meara, B. & Boyko, J. R package ‘corHMM’: hidden markov models of character evolution. (2020).Pagel, M. & Meade, A. BayesTraits, version 4. University of Reading, Berkshire, UK. http://www.evolution.rdg.ac.uk (2022).Pagel, M. & Meade, A. Bayesian analysis of correlated evolution of discrete characters by reversible-jump Markov chain Monte Carlo. Am. Nat. 167, 808–825 (2006).PubMed 

Google Scholar 
Maddison, W., Midford, P. & Otto, S. Estimating a binary character’s effect on speciation and extinction. Syst. Biol. 56, 701–710 (2007).PubMed 

Google Scholar 
FitzJohn, R. G., Maddison, W. P. & Otto, S. P. Estimating trait-dependent speciation and extinction rates from incompletely resolved phylogenies. Syst. Biol. 58, 595–611 (2009).PubMed 

Google Scholar 
Beaulieu, J. M. & O’Meara, B. C. Detecting hidden diversification shifts in models of trait-dependent speciation and extinction. Syst. Biol. 65, 583–601 (2016).PubMed 

Google Scholar 
Herrera-Alsina, L., van Els, P. & Etienne, R. S. Detecting the dependence of diversification on multiple traits from phylogenetic trees and trait data. Syst. Biol. 68, 317–328 (2018).
Google Scholar 
Strathmann, R. R. Hypotheses on the origins of marine larvae. Annu. Rev. Ecol. Syst. 24, 89–117 (1993).
Google Scholar 
Wray, G. A. Evolution of larvae and developmental modes. In Ecology of Marine Invertebrate Larvae. (ed. McEdward, L.) 413–447 (CRC Press, 1995).Raff, R. A. Origins of the other metazoan body plans: the evolution of larval forms. Philos. T R. Soc. B 363, 1473–1479 (2008).
Google Scholar 
Collin, R. & Moran, A. in Evolutionary Ecology of Marine Invertebrate Larvae 50–66 (Oxford University Press, 2018).Collin, R. & Miglietta, M. P. Reversing opinions on Dollo’s Law. Trends Ecol. Evol. 23, 602–609 (2008).PubMed 

Google Scholar 
Wiens, J. J. Re-evolution of lost mandibular teeth in frogs after more than 200 million years, and re-evaluating Dollo’s law. Evolution 65, 1283–1296 (2011).PubMed 

Google Scholar 
Wiens, J. J., Kuczynski, C. A., Duellman, W. E. & Reeder, T. W. Loss and re-evolution of complex life cycles in marsupial frogs: does ancestral trait reconstruction mislead? Evolution 61, 1886–1899 (2007).CAS 
PubMed 

Google Scholar 
Chippindale, P. T., Bonett, R. M., Baldwin, A. S. & Wiens, J. J. Phylogenetic evidence for a major reversal of life-history evolution in plethodontid salamanders. Evolution 58, 2809–2815 (2004).CAS 
PubMed 

Google Scholar 
Castroviejo-Fisher, S. et al. Phylogenetic systematics of egg-brooding frogs (Anura: Hemiphractidae) and the evolution of direct development. Zootaxa 4004, 1–75 (2015).PubMed 

Google Scholar 
Naumann, B., Schweiger, S., Hammel, J. U. & Müller, H. Parallel evolution of direct development in frogs—skin and thyroid gland development in African Squeaker Frogs (Anura: Arthroleptidae: Arthroleptis). Dev. Dyn. 250, 584–600 (2021).PubMed 

Google Scholar 
Goldberg, J., Taucce, P. P. G., Quinzio, S. I., Haddad, C. F. B. & Candioti, F. V. Increasing our knowledge on direct-developing frogs: the ontogeny of Ischnocnema henselii (Anura: Brachycephalidae). Zool. Anz. 284, 78–87 (2020).
Google Scholar 
Wassersug, R. J. & Duellman, W. E. Oral structures and their development in egg-brooding hylid frog embryos and larvae: evolutionary and ecological implications. J. Morphol. 182, 1–37 (1984).PubMed 

Google Scholar 
Kerney, R. R., Blackburn, D. C., Müller, H. & Hanken, J. Do larval traits re-evolve? Evidence from the embryogenesis of a direct-developing salamander, Plethodon cinereus. Evolution 66, 252–262 (2011).PubMed 

Google Scholar 
Theska, T. Musculoskeletal development of the Central African caecilian Idiocranium russeli (Amphibia: Gymnophiona: Indotyphlidae) and its bearing on the re-evolution of larvae in caecilian amphibians. Zoomorphology 138, 137–158 (2019).
Google Scholar 
Laslo, M., Denver, R. J. & Hanken, J. Evolutionary conservation of thyroid hormone receptor and deiodinase expression dynamics in ovo in a direct-developing frog, Eleutherodactylus coqui. Front. Endocrinol. 10, 307 (2019).
Google Scholar 
Gao, W. et al. Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity. Proc. Natl Acad. Sci. USA 116, 3646–3655 (2019).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Altig, R. & Crother, B. I. The evolution of three deviations from the biphasic anuran life cycle: alternatives to selection. Herpetol. Rev. 37, 321–356 (2006).
Google Scholar 
Venturelli, D. P., da Silva, W. R. & Giaretta, A. A. Tadpoles’ resistance to desiccation in species of Leptodactylus (Anura, Leptodactylidae). J. Herpetol. 55, 265–270 (2021).
Google Scholar 
Seymour, R. S. Respiration of aquatic and terrestrial amphibian embryos. American Zoologist 39, 261–270 (1999).Blackburn, D. G. Convergent evolution of viviparity, matrotrophy, and specializations for fetal nutrition in reptiles and other vertebrates. Am. Zool. 32, 313–321 (1992).
Google Scholar 
Buckley, D., Alcobendas, M., García-París, M. & Wake, M. H. Heterochrony, cannibalism, and the evolution of viviparity in Salamandra salamandra. Evol. Dev. 9, 105–115 (2007).PubMed 

Google Scholar 
Kusrini, M. D., Rowley, J. J. L., Khairunnisa, L. R., Shea, G. M. & Altig, R. The reproductive biology and larvae of the first tadpole-bearing frog, Limnonectes larvaepartus. PLoS ONE 10, e116154–e116159 (2015).ADS 
PubMed 
PubMed Central 

Google Scholar 
Lanza, B. & Leo, P. Sul primo caso sicuro di riproduzione vivipara nel genere Speleomantes. 1–54 (2000).Lunghi, E. et al. Comparative reproductive biology of European cave salamanders (Genus Hydromantes): nesting selection and multiple annual breeding. Salamandra 54, 101–108 (2018).
Google Scholar 
Liedtke, H. C. et al. Terrestrial reproduction as an adaptation to steep terrain in African toads. Proc. R. Soc. Lond. B 284, 20162598–20162599 (2017).
Google Scholar 
Wake, M. H. The reproductive biology of Eleutherodactylus jasperi (Amphibia, Anura, Leptodactylidae), with comments on the evolution of live-bearing systems. J. Herpetol. 12, 121–133 (1978).
Google Scholar 
Jennings, D. H. & Hanken, J. Mechanistic basis of life history evolution in anuran amphibians: thyroid gland development in the direct-developing frog, Eleutherodactylus coqui. Gen. Comp. Endocr. 111, 225–232 (1998).CAS 
PubMed 

Google Scholar 
Callery, E. M. & Elinson, R. P. Thyroid hormone-dependent metamorphosis in a direct developing frog. Proc. Natl Acad. Sci. USA 97, 2615–2620 (2000).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Callery, E. M., Hung, F. & Elinson, R. P. Frogs without polliwogs: evolution of anuran direct development. BioEssays 23, 233–241 (2001).CAS 
PubMed 

Google Scholar 
Buckley, L. B. & Jetz, W. Environmental and historical constraints on global patterns of amphibian richness. Proc. R. Soc. Lond. B 274, 1167–1173 (2007).
Google Scholar 
Pyron, R. A. & Wiens, J. J. Large-scale phylogenetic analyses reveal the causes of high tropical amphibian diversity. Proc. R. Soc. Lond. B 280, 20131622 (2013).
Google Scholar 
Gómez-Rodríguez, C., Baselga, A. & Wiens, J. J. Is diversification rate related to climatic niche width? Glob. Ecol. Biogeogr. 24, 383–395 (2015).
Google Scholar 
Moen, D. S. & Wiens, J. J. Microhabitat and climatic niche change explain patterns of diversification among frog families. Am. Nat. 190, 29–44 (2017).PubMed 

Google Scholar 
Kozak, K. H. & Wiens, J. J. Accelerated rates of climatic-niche evolution underlie rapid species diversification. Ecol. Lett. 13, 1378–1389 (2010).PubMed 

Google Scholar 
Jaramillo, A. F. et al. Vastly underestimated species richness of Amazonian salamanders (Plethodontidae: Bolitoglossa) and implications about plethodontid diversification. Mol. Phylogenet. Evol. 149, 106841 (2020).PubMed 

Google Scholar 
Jetz, W. & Pyron, R. A. The interplay of past diversification and evolutionary isolation with present imperilment across the amphibian tree of life. Nat. Ecol. Evol. 2, 850–858 (2018).PubMed 

Google Scholar 
Bars-Closel, M., Kohlsdorf, T., Moen, D. S. & Wiens, J. J. Diversification rates are more strongly related to microhabitat than climate in squamate reptiles (lizards and snakes). Evolution 71, 2243–2261 (2017).PubMed 

Google Scholar 
Cyriac, V. P. & Kodandaramaiah, U. Digging their own macroevolutionary grave: fossoriality as an evolutionary dead end in snakes. J. Evolution. Biol. 31, 587–598 (2018).CAS 

Google Scholar 
Zamudio, K. R., Bell, R. C., Nali, R. C., Haddad, C. F. B. & Prado, C. P. A. Polyandry, predation, and the evolution of frog reproductive modes. Am. Nat. 188, S41–S61 (2016).PubMed 

Google Scholar 
Lion, M. B. et al. Global patterns of terrestriality in amphibian reproduction. Glob. Ecol. Biogeogr. 4, 679–13 (2019).
Google Scholar 
Müller, H. et al. Forests as promoters of terrestrial life-history strategies in East African amphibians. Biol. Lett. 9, 20121146–20121146 (2013).PubMed 
PubMed Central 

Google Scholar 
Velo-Antón, G., García-París, M., Galán, P. & Cordero Rivera, A. The evolution of viviparity in Holocene islands: ecological adaptation versus phylogenetic descent along the transition from aquatic to terrestrial environments. J. Zool. Syst. Evol. Res. 45, 345–352 (2007).
Google Scholar 
Liedtke, H. C. AmphiNom: an amphibian systematics tool. Syst. Biodivers. 17, 1–6 (2019).
Google Scholar 
R Core Team. R: a language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria, 2020). http://www.R-project.org/.IUCN. IUCN 2020. The IUCN Red List of Threatened Species. Version 2020-2. https://www.iucnredlist.org (2019).AmphibiaChina. The database of Chinese amphibians. Electronic Database (Kunming Institute of Zoology (CAS), Kunming, Yunnan, China, 2019) http://www.amphibiachina.org/.Ron, S. R., Yanez-Muñoz, M. H., Merino-Viteri, A. & Ortiz, D. A. Anfibios del Ecuador. Version 2019.0 (Museo de Zoología, Pontificia Universidad Católica del Ecuador, 2019). https://bioweb.bio/faunaweb/amphibiaweb.Greven, H. In Reproductive Biology and Phylogeny of Urodela 447–475 (Taylor & Francis, 2003).Marks, S. B. & Collazo, A. Direct development in Desmognathus aeneus (Caudata: Plethodontidae): a staging table. Copeia 1998, 637–648 (1998).
Google Scholar 
Müller, H., Loader, S. P., Ngalason, W., Howell, K. M. & Gower, D. J. Reproduction in brevicipitid frogs (Amphibia: Anura: Brevicipitidae)—evidence from Probreviceps m. macrodactylus. Copeia 2007, 726–733 (2007).
Google Scholar 
Velo-Antón, G., Santos, X., Sanmartín-Villar, I., Cordero-Rivera, A. & Buckley, D. Intraspecific variation in clutch size and maternal investment in pueriparous and larviparous Salamandra salamandra females. Evol. Ecol. 29, 185–204 (2014).
Google Scholar 
Beaulieu, J. M., O’Meara, B. C. & Donoghue, M. J. Identifying hidden rate changes in the evolution of a binary morphological character: the evolution of plant habit in campanulid angiosperms. Systematic biology 62, 725–737 (2013).Pupko, T., Pe’er, I., Shamir, R. & Graur, D. A fast algorithm for joint reconstruction of ancestral amino acid sequences. Mol. Biol. Evol. 17, 890–896 (2000).CAS 
PubMed 

Google Scholar 
Bollback, J. P. SIMMAP: stochastic character mapping of discrete traits on phylogenies. BMC Bioinforma. 7, 88 (2006).
Google Scholar 
Plummer, M., Best, N., Cowles, K. & Vines, K. CODA: convergence diagnosis and output analysis for MCMC. R. News 6, 7–11 (2006).
Google Scholar 
Tuffley, C. & Steel, M. Modeling the covarion hypothesis of nucleotide substitution. Math. Biosci. 147, 63–91 (1998).MathSciNet 
CAS 
PubMed 
MATH 

Google Scholar 
Xie, W., Lewis, P. O., Fan, Y., Kuo, L. & Chen, M. H. Improving marginal likelihood estimation for Bayesian phylogenetic model selection. Syst. Biol. 60, 150–160 (2011).PubMed 

Google Scholar 
Nakov, T., Beaulieu, J. M. & Alverson, A. J. Diatoms diversify and turn over faster in freshwater than marine environments. Evolution 73, 2497–2511 (2019).PubMed 

Google Scholar 
Etienne, R. S. et al. Diversity-dependence brings molecular phylogenies closer to agreement with the fossil record. Proc. R. Soc. Lond. B 279, 1300–1309 (2011).
Google Scholar  More

PTM concentrationsIn Table 1, the descriptive statistics of PTMs in 50 dust samples from Qom city are described. The background values were used based on the concentrations of metals in the Upper Continental Crust. The mean concentration of As, Cd, Cu, Mo, Pb, Sb, and, Zn exceeded the background value. Also, Cd, Cu, Mo, Pb, Sb, and, Zn had a coefficient of variation (C.V.) greater than 50%, indicating a severe variability in PTMs concentrations in the atmospheric dust of the studied area2. Metals with C.V.  Pb  > Zn  > As  > Cd  > Cu  > Mo  > Cr  > Mn  > Ni = Co. Antimony (38.55) and Pb (35.13) had the highest average EF values, which means they were enriched very high in the windowsill dust. Also, they had a wide range of EF values in the 50 stations: from 4.0 to 227.0 for Sb, and from 8.3 to 140.8 for Pb which might reflect the existence of discrete multiple sources in the studied area. The degree of enrichment for Pb and Sb in the industrial sector was extreme and in the commercial sector was very high; also, the other sectors were significantly enriched. Zinc and As had a more homogenous enrichment in the area. In all the functional sectors, 95% and 84% of stations were significantly enriched by As and Zn, respectively. Copper, Cd, and Mo were moderately enriched in all functional sectors, but the greenspace sector had minimal enrichment by these elements. Some areas in the industrial sector had significant to very high enrichment of Cd. The EF value indicated Co, Cr, Mn, and Ni were minimally enriched in all the stations.Figure 2Box plot of the (a) enrichment factor (EF), and (b) geo accumulation index (Igeo) for the dust samples in the studied area.Full size imageThe highest average values of Igeo were obtained in the order of Pb  > Sb  > As  > Zn. PTMs included Co, Cr, Ni, Fe, and Mn were categorized as unpolluted and Cd, Cu, and Mo were in the category of unpolluted to moderately polluted. In the industrial zone, the windowsill dust was extremely polluted with Sb and Pb. The sequence of contamination intensity with Pb, Zn and Sb according to land use was: industrial  > commercial  > residential  > greenspace. The highest concentration of arsenic in the study area belongs to the industrial area.To evaluate the pollution level based on land use, PLI and mCd indices were utilized (Fig. 3). These cumulative indices showed that the dust in Qom city is considerably contaminated with PTMs. According to the PLI index, all the stations were categorized as polluted sites. The PLIzone values were in the order of industrial (3.77)  > commercial (2.05)  > residential (1.67)  > green space (1.38). This pattern was also repeated with the mCd index. The mCd for the industrial sector ranged from 6.98 (high contamination level) to 39.60 (ultra-high contamination level). In the commercial sector, fifty percent of dust samples were classified as having a high degree of contamination. All the greenspace stations were in the moderate pollution category. This shows the possible effect of tree density in diminishing the risk of dust pollution to the receptors.Figure 3Pollution level indexes (a) mCd and (b) PLI, based on four functional areas.Full size imageSpatial distribution of PTMsThe As, Cd, Cu, Sb, Pb, and Zn content in 100% of the dust samples exceeded the background value. Spatial distribution maps were generated for the hotspot PTMs (As, Sb, Pb, Cd, Cu, Mo and Zn) by applying the inverse distance weighted (IDW) interpolation method (ArcGIS 10.3). Figure 4 demonstrates that PTMs dispersions were slightly influenced by the prevailing wind direction (from the west), suggesting they came from the point- or area- sources. On the other hand, the K–S test showed that the overall distribution of PTMs was not normal in the studied region. This might signify the influence of industrial activities and the presence of multiple sources of dust.Figure 4Spatial distribution maps of seven PTMs in windowsill dusts of Qom, Iran. This map was constructed using ArcGIS version 10.3. (https://www.esri.com/en-us/arcgis/products/arcgisdesktop/overview).Full size imageThe highest pollution load of PTMs belonged to the industrial section. The level of pollution gradually decreased from Shokouhieh to Mahmoudabad industrial zones. The reason is related to more active industries, a closed environment, and more construction existing in Shokouhieh industrial town than in Mahmoudabad industrial town.There is a clear decreasing trend from the central part to southern (downtown area) and southwestern (suburb area) parts of the city. In fact, these parts are diffusely populated and the southwestern part is almost new with lots of barren lands. Copper, Mo, and Cd show high concentrations toward the central part of the city. Educational, cultural and commercial activities are mainly located in the central part of the city. Also, historical and religious districts in the city center are accompanied by a huge influx of tourists throughout the year. For this reason, the central part of the city has various public transportations such as bus stands and taxi stations, and is dominated by a high load of motorcycles.In the eastern part of the city, some hotspots can be observed (Fig. 4). This part includes an important transportation system (like highways and a complex interchange) where exhaust traffic emissions might be a probable source of As, Sb, Pb, and Zn. Unlike Pb and Zn, several peaks of As are scattered in the western part, suggesting an area source might exist in the region. It is noteworthy that the western area is densely populated with lots of residential buildings. Bisht et al. (2022)35 also observed hotspots of As in the residential area of Dehradun, India.PTM potential sourcesTo evaluate the relationship between PTMs in dust samples, the Spearman correlation and PCA were developed (Fig. 5) and more details are given in Table S7. Statistical analysis can help to identify the potential source of contamination in urban dust. The Spearman correlation was significant at p  residential ≈ greenspace. The five PTMs with the highest overall HI are ranked as follows: Pb  > As  > Cr  > Mn  > Sb (Fig. 7). The HI values in all the sections were lower than the permissible level (1.00), except for Pb. In the industrial section, Pb recorded the highest HI value for children (HI = 1.73) which exceeded the acceptable value. The HI values were 10 times higher for children than adults indicating they are more susceptible to PTMs in the dust.The dominant pathway for noncancerous risk was ingestion followed by dermal contact and inhalation. The trend is in line with previous research25,51,52. However, for Co and Mn, the descending order was different as follows ingestion  > inhalation  > dermal contact. The highest contribution of HQinh and HQderm to HI was measured for Co (34.0%) and Cd (29.0%), respectively.In this study, the carcinogenic risk from windowsill dust was estimated for the carcinogens including Cd, Co, Cr and Ni, Pb, and As through the possible routes (Fig. 8, Table S9). The contribution of PTMs to CR decreased in the order of Cr (3.24E−05)  > As (2.05E−05)  > Pb (2.52E−06)  > Co (6.91E−09)  > Ni (1.72E−09)  > Cd (2.58E−10). The average CR values for target PTMs through inhalation ranged from 7.9E−10 to 1.7E−07, which remained in the safety zone (CR  inhalation  > dermal (Fig. 8). While the contribution of Cr to carcinogenic risk was higher through inhalation than ingestion. The reports concluded that the primary exposure route of Cr is inhalation54. Considering the predominant forms of Cr in the environment, CrVI is more toxic than CrIII. Exposure to CrVI can cause immunological diseases, dental effects and carcinogenic effects (lung cancer, nose and nasal sinus cancer, suspected laryngeal and stomach cancers)54,55.The result of health risk from target PTMs in windowsills of Qom indicates significant chronic exposure to Pb can take place for children in the industrial zone. The ingestion route is the most probable pathway for children due to their hand–to–mouth behavior56. Lead can bio-accumulate in the body without any obvious symptoms of toxicity56. The total CR values for Pb, Cr and As in different land-use types were in the range of tolerable carcinogenic risk (1 × 10−4  More

.readcube-buybox { display: none !important;}
Protecting an overfished lobster species helps the crustaceans to grow big, according to an analysis of European lobsters in marine sanctuaries1.

Access options

/* style specs start */
style{display:none!important}.LiveAreaSection-193358632 *{align-content:stretch;align-items:stretch;align-self:auto;animation-delay:0s;animation-direction:normal;animation-duration:0s;animation-fill-mode:none;animation-iteration-count:1;animation-name:none;animation-play-state:running;animation-timing-function:ease;azimuth:center;backface-visibility:visible;background-attachment:scroll;background-blend-mode:normal;background-clip:borderBox;background-color:transparent;background-image:none;background-origin:paddingBox;background-position:0 0;background-repeat:repeat;background-size:auto auto;block-size:auto;border-block-end-color:currentcolor;border-block-end-style:none;border-block-end-width:medium;border-block-start-color:currentcolor;border-block-start-style:none;border-block-start-width:medium;border-bottom-color:currentcolor;border-bottom-left-radius:0;border-bottom-right-radius:0;border-bottom-style:none;border-bottom-width:medium;border-collapse:separate;border-image-outset:0s;border-image-repeat:stretch;border-image-slice:100%;border-image-source:none;border-image-width:1;border-inline-end-color:currentcolor;border-inline-end-style:none;border-inline-end-width:medium;border-inline-start-color:currentcolor;border-inline-start-style:none;border-inline-start-width:medium;border-left-color:currentcolor;border-left-style:none;border-left-width:medium;border-right-color:currentcolor;border-right-style:none;border-right-width:medium;border-spacing:0;border-top-color:currentcolor;border-top-left-radius:0;border-top-right-radius:0;border-top-style:none;border-top-width:medium;bottom:auto;box-decoration-break:slice;box-shadow:none;box-sizing:border-box;break-after:auto;break-before:auto;break-inside:auto;caption-side:top;caret-color:auto;clear:none;clip:auto;clip-path:none;color:initial;column-count:auto;column-fill:balance;column-gap:normal;column-rule-color:currentcolor;column-rule-style:none;column-rule-width:medium;column-span:none;column-width:auto;content:normal;counter-increment:none;counter-reset:none;cursor:auto;display:inline;empty-cells:show;filter:none;flex-basis:auto;flex-direction:row;flex-grow:0;flex-shrink:1;flex-wrap:nowrap;float:none;font-family:initial;font-feature-settings:normal;font-kerning:auto;font-language-override:normal;font-size:medium;font-size-adjust:none;font-stretch:normal;font-style:normal;font-synthesis:weight style;font-variant:normal;font-variant-alternates:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-position:normal;font-weight:400;grid-auto-columns:auto;grid-auto-flow:row;grid-auto-rows:auto;grid-column-end:auto;grid-column-gap:0;grid-column-start:auto;grid-row-end:auto;grid-row-gap:0;grid-row-start:auto;grid-template-areas:none;grid-template-columns:none;grid-template-rows:none;height:auto;hyphens:manual;image-orientation:0deg;image-rendering:auto;image-resolution:1dppx;ime-mode:auto;inline-size:auto;isolation:auto;justify-content:flexStart;left:auto;letter-spacing:normal;line-break:auto;line-height:normal;list-style-image:none;list-style-position:outside;list-style-type:disc;margin-block-end:0;margin-block-start:0;margin-bottom:0;margin-inline-end:0;margin-inline-start:0;margin-left:0;margin-right:0;margin-top:0;mask-clip:borderBox;mask-composite:add;mask-image:none;mask-mode:matchSource;mask-origin:borderBox;mask-position:0 0;mask-repeat:repeat;mask-size:auto;mask-type:luminance;max-height:none;max-width:none;min-block-size:0;min-height:0;min-inline-size:0;min-width:0;mix-blend-mode:normal;object-fit:fill;object-position:50% 50%;offset-block-end:auto;offset-block-start:auto;offset-inline-end:auto;offset-inline-start:auto;opacity:1;order:0;orphans:2;outline-color:initial;outline-offset:0;outline-style:none;outline-width:medium;overflow:visible;overflow-wrap:normal;overflow-x:visible;overflow-y:visible;padding-block-end:0;padding-block-start:0;padding-bottom:0;padding-inline-end:0;padding-inline-start:0;padding-left:0;padding-right:0;padding-top:0;page-break-after:auto;page-break-before:auto;page-break-inside:auto;perspective:none;perspective-origin:50% 50%;pointer-events:auto;position:static;quotes:initial;resize:none;right:auto;ruby-align:spaceAround;ruby-merge:separate;ruby-position:over;scroll-behavior:auto;scroll-snap-coordinate:none;scroll-snap-destination:0 0;scroll-snap-points-x:none;scroll-snap-points-y:none;scroll-snap-type:none;shape-image-threshold:0;shape-margin:0;shape-outside:none;tab-size:8;table-layout:auto;text-align:initial;text-align-last:auto;text-combine-upright:none;text-decoration-color:currentcolor;text-decoration-line:none;text-decoration-style:solid;text-emphasis-color:currentcolor;text-emphasis-position:over right;text-emphasis-style:none;text-indent:0;text-justify:auto;text-orientation:mixed;text-overflow:clip;text-rendering:auto;text-shadow:none;text-transform:none;text-underline-position:auto;top:auto;touch-action:auto;transform:none;transform-box:borderBox;transform-origin:50% 50%0;transform-style:flat;transition-delay:0s;transition-duration:0s;transition-property:all;transition-timing-function:ease;vertical-align:baseline;visibility:visible;white-space:normal;widows:2;width:auto;will-change:auto;word-break:normal;word-spacing:normal;word-wrap:normal;writing-mode:horizontalTb;z-index:auto;-webkit-appearance:none;-moz-appearance:none;-ms-appearance:none;appearance:none;margin:0}.LiveAreaSection-193358632{width:100%}.LiveAreaSection-193358632 .login-option-buybox{display:block;width:100%;font-size:17px;line-height:30px;color:#222;padding-top:30px;font-family:Harding,Palatino,serif}.LiveAreaSection-193358632 .additional-access-options{display:block;font-weight:700;font-size:17px;line-height:30px;color:#222;font-family:Harding,Palatino,serif}.LiveAreaSection-193358632 .additional-login >li:not(:first-child)::before{transform:translateY(-50%);content:””;height:1rem;position:absolute;top:50%;left:0;border-left:2px solid #999}.LiveAreaSection-193358632 .additional-login >li:not(:first-child){padding-left:10px}.LiveAreaSection-193358632 .additional-login >li{display:inline-block;position:relative;vertical-align:middle;padding-right:10px}.BuyBoxSection-683559780{display:flex;flex-wrap:wrap;flex:1;flex-direction:row-reverse;margin:-30px -15px 0}.BuyBoxSection-683559780 .box-inner{width:100%;height:100%}.BuyBoxSection-683559780 .readcube-buybox{background-color:#f3f3f3;flex-shrink:1;flex-grow:1;flex-basis:255px;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .subscribe-buybox{background-color:#f3f3f3;flex-shrink:1;flex-grow:4;flex-basis:300px;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .subscribe-buybox-nature-plus{background-color:#f3f3f3;flex-shrink:1;flex-grow:4;flex-basis:100%;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .title-readcube{display:block;margin:0;margin-right:20%;margin-left:20%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .title-buybox{display:block;margin:0;margin-right:29%;margin-left:29%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .title-asia-buybox{display:block;margin:0;margin-right:5%;margin-left:5%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .asia-link{color:#069;cursor:pointer;text-decoration:none;font-size:1.05em;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:1.05em6}.BuyBoxSection-683559780 .access-readcube{display:block;margin:0;margin-right:10%;margin-left:10%;font-size:14px;color:#222;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .access-asia-buybox{display:block;margin:0;margin-right:5%;margin-left:5%;font-size:14px;color:#222;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .access-buybox{display:block;margin:0;margin-right:30%;margin-left:30%;font-size:14px;color:#222;opacity:.8px;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .usps-buybox{display:block;margin:0;margin-right:30%;margin-left:30%;font-size:14px;color:#222;opacity:.8px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .price-buybox{display:block;font-size:30px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;padding-top:30px;text-align:center}.BuyBoxSection-683559780 .price-from{font-size:14px;padding-right:10px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .issue-buybox{display:block;font-size:13px;text-align:center;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:19px}.BuyBoxSection-683559780 .no-price-buybox{display:block;font-size:13px;line-height:18px;text-align:center;padding-right:10%;padding-left:10%;padding-bottom:20px;padding-top:30px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif}.BuyBoxSection-683559780 .vat-buybox{display:block;margin-top:5px;margin-right:20%;margin-left:20%;font-size:11px;color:#222;padding-top:10px;padding-bottom:15px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:17px}.BuyBoxSection-683559780 .button-container{display:flex;padding-right:20px;padding-left:20px;justify-content:center}.BuyBoxSection-683559780 .button-container >*{flex:1px}.BuyBoxSection-683559780 .button-container >a:hover,.Button-505204839:hover,.Button-1078489254:hover,.Button-2808614501:hover{text-decoration:none}.BuyBoxSection-683559780 .readcube-button{background:#fff;margin-top:30px}.BuyBoxSection-683559780 .button-asia{background:#069;border:1px solid #069;border-radius:0;cursor:pointer;display:block;padding:9px;outline:0;text-align:center;text-decoration:none;min-width:80px;margin-top:75px}.BuyBoxSection-683559780 .button-label-asia,.ButtonLabel-3869432492,.ButtonLabel-3296148077,.ButtonLabel-1566022830{display:block;color:#fff;font-size:17px;line-height:20px;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;text-align:center;text-decoration:none;cursor:pointer}.Button-505204839,.Button-1078489254,.Button-2808614501{background:#069;border:1px solid #069;border-radius:0;cursor:pointer;display:block;padding:9px;outline:0;text-align:center;text-decoration:none;min-width:80px;max-width:320px;margin-top:10px}.Button-505204839 .readcube-label,.Button-1078489254 .readcube-label,.Button-2808614501 .readcube-label{color:#069}
/* style specs end */Subscribe to Nature+Get immediate online access to Nature and 55 other Nature journal$29.99monthlySubscribe to JournalGet full journal access for 1 year$199.00only $3.90 per issueAll prices are NET prices.VAT will be added later in the checkout.Tax calculation will be finalised during checkout.Buy articleGet time limited or full article access on ReadCube.$32.00All prices are NET prices.

Additional access options:

doi: https://doi.org/10.1038/d41586-022-03708-2

References

Subjects

Conservation biology More

Bardgett, R. D. et al. Combatting global grassland degradation. Nat. Rev. Earth Environ. 2, 720–735 (2021)Reichstein, M. et al. Climate extremes and the carbon cycle. Nature 500, 287–295 (2013).ADS 
CAS 
PubMed 

Google Scholar 
Seneviratne, S. I. et al. Weather and climate extreme events in a changing climate. In Climate change 2021: the physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2021).Pretty, J. et al. Global assessment of agricultural system redesign for sustainable intensification. Nat. Sustain 1, 441–446 (2018).
Google Scholar 
Allan, E. et al. Land use intensification alters ecosystem multifunctionality via loss of biodiversity and changes to functional composition. Ecol. Lett. 18, 834–843 (2015).PubMed 
PubMed Central 

Google Scholar 
Bardgett, R. D. & Cook, R. Functional aspects of soil animal diversity in agricultural grasslands. Appl. Soil Ecol. 10, 263–276 (1998).
Google Scholar 
Postma-Blaauw, M. B., de Goede, R. G. M., Bloem, J., Faber, J. H. & Brussaard, L. Soil biota community structure and abundance under agricultural intensification and extensification. Ecology 91, 460–473 (2010).PubMed 

Google Scholar 
Vályi, K., Rillig, M. C. & Hempel, S. Land-use intensity and host plant identity interactively shape communities of arbuscular mycorrhizal fungi in roots of grassland plants. N. Phytologist 205, 1577–1586 (2015).
Google Scholar 
de Vries, F. T., Hoffland, E., van Eekeren, N., Brussaard, L. & Bloem, J. Fungal/bacterial ratios in grasslands with contrasting nitrogen management. Soil Biol. Biochem. 38, 2092–2103 (2006).
Google Scholar 
de Vries, F. T. et al. Extensive Management Promotes Plant and Microbial Nitrogen Retention in Temperate Grassland. PLoS ONE 7, e51201 (2012).ADS 
PubMed 
PubMed Central 

Google Scholar 
de Vries, F. T., van Groenigen, J. W., Hoffland, E. & Bloem, J. Nitrogen losses from two grassland soils with different fungal biomass. Soil Biol. Biochem. 43, 997–1005 (2011).
Google Scholar 
Malik, A. A. et al. Soil fungal: bacterial ratios are linked to altered carbon cycling. Front. Microbiol. 7, 1247 (2016).Bardgett, R. D., Streeter, T. C. & Bol, R. Soil Microbes Compete Effectively with Plants for Organic-Nitrogen Inputs to Temperate Grasslands. Ecology 84, 1277–1287 (2003).
Google Scholar 
Bardgett, R. D. & McAlister, E. The measurement of soil fungal:bacterial biomass ratios as an indicator of ecosystem self-regulation in temperate meadow grasslands. Biol. Fertil. Soils 29, 282–290 (1999).
Google Scholar 
Gordon, H., Haygarth, P. M. & Bardgett, R. D. Drying and rewetting effects on soil microbial community composition and nutrient leaching. Soil Biol. Biochem. 40, 302–311 (2008).CAS 

Google Scholar 
Duffy, J. E. et al. The functional role of biodiversity in ecosystems: incorporating trophic complexity. Ecol. Lett. 10, 522–538 (2007).PubMed 

Google Scholar 
Wang, S. & Brose, U. Biodiversity and ecosystem functioning in food webs: the vertical diversity hypothesis. Ecol. Lett. 21, 9–20 (2018).ADS 
CAS 
PubMed 

Google Scholar 
Ruf, A., Kuzyakov, Y. & Lopatovskaya, O. Carbon fluxes in soil food webs of increasing complexity revealed by C-14 labelling and C-13 natural abundance. Soil Biol. Biochem. 38, 2390–2400 (2006).CAS 

Google Scholar 
Pollierer, M. M., Langel, R., Koerner, C., Maraun, M. & Scheu, S. The underestimated importance of belowground carbon input for forest soil animal food webs. Ecol. Lett. 10, 729–736 (2007).PubMed 

Google Scholar 
Eissfeller, V. et al. Incorporation of plant carbon and microbial nitrogen into the rhizosphere food web of beech and ash. Soil Biol. Biochem. 62, 76–81 (2013).CAS 

Google Scholar 
Gilbert, K. J. et al. Exploring carbon flow through the root channel in a temperate forest soil food web. Soil Biol. Biochem. 76, 45–52 (2014).CAS 

Google Scholar 
Goncharov, A. A., Tsurikov, S. M., Potapov, A. M. & Tiunov, A. V. Short-term incorporation of freshly fixed plant carbon into the soil animal food web: field study in a spruce forest. Ecol. Res. 31, 923–933 (2016).CAS 

Google Scholar 
Chomel, M. et al. Drought decreases incorporation of recent plant photosynthate into soil food webs regardless of their trophic complexity. Glob. Change Biol. 25, 3549–3561 (2019).ADS 

Google Scholar 
Moore, J. C., de Ruiter, P. C. & Hunt, H. W. Influence of productivity on the stability of real and model ecosystems. Science 261, 906–908 (1993).ADS 
CAS 
PubMed 
MATH 

Google Scholar 
de Ruiter, P. C., Neutel, A.-M. & Moore, J. C. Energetics, Patterns of Interaction Strengths, and Stability in Real Ecosystems. Science 269, 1257–1260 (1995).ADS 
PubMed 

Google Scholar 
Rooney, N., McCann, K., Gellner, G. & Moore, J. C. Structural asymmetry and the stability of diverse food webs. Nature 442, 265–269 (2006).ADS 
CAS 
PubMed 

Google Scholar 
Rooney, N. & McCann, K. S. Integrating food web diversity, structure and stability. Trends Ecol. Evolution 27, 40–46 (2012).
Google Scholar 
de Vries, F. T. et al. Land use alters the resistance and resilience of soil food webs to drought. Nat. Clim. Change 2, 276 (2012).ADS 

Google Scholar 
Ingrisch, J. et al. Land Use Alters the Drought Responses of Productivity and CO2 Fluxes in Mountain Grassland. Ecosystems 21, 689–703 (2018).PubMed 

Google Scholar 
Karlowsky, S. et al. Land use in mountain grasslands alters drought response and recovery of carbon allocation and plant‐microbial interactions. J. Ecol. 106, 1230–1243 (2017).PubMed 
PubMed Central 

Google Scholar 
Vilonen, L., Ross, M. & Smith, M. D. What happens after drought ends: synthesizing terms and definitions. N. Phytologist 235, 420–431 (2022).
Google Scholar 
Ingrisch, J., Karlowsky, S., Hasibeder, R., Gleixner, G. & Bahn, M. Drought and recovery effects on belowground respiration dynamics and the partitioning of recent carbon in managed and abandoned grassland. Glob. Change Biol. 26, 4366–4378 (2020).ADS 

Google Scholar 
Ward, S. E. et al. Legacy effects of grassland management on soil carbon to depth. Glob. Change Biol. 22, 2929–2938 (2016).ADS 

Google Scholar 
Henry, C. et al. A stomatal safety-efficiency trade-off constrains responses to leaf dehydration. Nat. Commun. 10, 3398 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Baptist, F. et al. 13C and 15N allocations of two alpine species from early and late snowmelt locations reflect their different growth strategies. J. Exp. Bot. 60, 2725–2735 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
Bergmann, J. et al. The fungal collaboration gradient dominates the root economics space in plants. Sci. Adv. 6, eaba3756 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Williams, A. et al. Root functional traits explain root exudation rate and composition across a range of grassland species. J. Ecol. 110, 21–33 (2022).
Google Scholar 
Deyn, G. B. D., Quirk, H., Oakley, S., Ostle, N. J. & Bartgett, R. D. Rapid transfer of photosynthetic carbon through the plant-soil system in differently managed species-rich grasslands. Biogeosciences 8, 1131–1139 (2011).Pausch, J. et al. Small but active – pool size does not matter for carbon incorporation in below‐ground food webs.Functional Ecol. 30, 479–489 (2016).
Google Scholar 
Morriën, E. et al. Soil networks become more connected and take up more carbon as nature restoration progresses. Nat. Commun. 8, 14349 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Li, Z. et al. The flux of root-derived carbon via fungi and bacteria into soil microarthropods (Collembola) differs markedly between cropping systems. Soil Biol. Biochem. 160, 108336 (2021).CAS 

Google Scholar 
Joergensen, R. Ergosterol and microbial biomass in the rhizosphere of grassland soils. Soil Biol. Biogeochemistry 32, 647–652 (2000).CAS 

Google Scholar 
Staddon, P. L., Ramsey, C. B., Ostle, N., Ineson, P. & Fitter, A. H. Rapid Turnover of Hyphae of Mycorrhizal Fungi Determined by AMS Microanalysis of 14C. Science 300, 1138–1140 (2003).ADS 
CAS 
PubMed 

Google Scholar 
Johnson, D., Leake, J. R., Ostle, N., Ineson, P. & Read, D. J. In situ 13CO2 pulse-labelling of upland grassland demonstrates a rapid pathway of carbon flux from arbuscular mycorrhizal mycelia to the soil. N. Phytologist 153, 327–334 (2002).CAS 

Google Scholar 
Johnson, D., Leake, J. R. & Read, D. J. Transfer of recent photosynthate into mycorrhizal mycelium of an upland grassland: short-term respiratory losses and accumulation of C-14. Soil Biol. Biochem. 34, 1521–1524 (2002).CAS 

Google Scholar 
Schimel, J., Balser, T. C. & Wallenstein, M. Microbial Stress-Response Physiology and Its Implications for Ecosystem Function. Ecology 88, 1386–1394 (2007).PubMed 

Google Scholar 
Strickland, M. S. & Rousk, J. Considering fungal:bacterial dominance in soils – Methods, controls, and ecosystem implications. Soil Biol. Biochem. 42, 1385–1395 (2010).CAS 

Google Scholar 
Manzoni, S., Schimel, J. P. & Porporato, A. Responses of soil microbial communities to water stress: results from a meta-analysis. Ecology 93, 930–938 (2012).PubMed 

Google Scholar 
Holden, S. R. & Treseder, K. K. A meta-analysis of soil microbial biomass responses to forest disturbances. Front Microbiol 4, 163 (2013).PubMed 
PubMed Central 

Google Scholar 
Guhr, A., Borken, W., Spohn, M. & Matzner, E. Redistribution of soil water by a saprotrophic fungus enhances carbon mineralization. PNAS 112, 14647–14651 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
de Vries, F. T. et al. Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9, 3033 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
Allen, M. F. Mycorrhizal Fungi: Highways for Water and Nutrients in Arid Soils. Vadose Zone J. 6, 291–297 (2007).
Google Scholar 
Kakouridis, A. et al. Routes to roots: direct evidence of water transport by arbuscular mycorrhizal fungi to host plants. bioRxiv https://doi.org/10.1101/2020.09.21.305409 (2020).Leake, J. R., Ostle, N. J., Rangel-Castro, J. I. & Johnson, D. Carbon fluxes from plants through soil organisms determined by field 13CO2 pulse-labelling in an upland grassland. Appl. Soil Ecol. 33, 152–175 (2006).
Google Scholar 
Maaß, S., Migliorini, M., Rillig, M. C. & Caruso, T. Disturbance, neutral theory, and patterns of beta diversity in soil communities. Ecol. Evolution 4, 4766–4774 (2014).
Google Scholar 
Barnard, R. L., Osborne, C. A. & Firestone, M. K. Changing precipitation pattern alters soil microbial community response to wet-up under a Mediterranean-type climate. ISME J. 9, 946–957 (2015).CAS 
PubMed 

Google Scholar 
Lennon, J. T. & Jones, S. E. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol 9, 119–130 (2011).CAS 
PubMed 

Google Scholar 
Meisner, A., Bååth, E. & Rousk, J. Microbial growth responses upon rewetting soil dried for four days or one year. Soil Biol. Biochem. 66, 188–192 (2013).CAS 

Google Scholar 
Meisner, A., Rousk, J. & Bååth, E. Prolonged drought changes the bacterial growth response to rewetting. Soil Biol. Biochem. 88, 314–322 (2015).CAS 

Google Scholar 
Blazewicz, S. J., Schwartz, E. & Firestone, M. K. Growth and death of bacteria and fungi underlie rainfall-induced carbon dioxide pulses from seasonally dried soil. Ecology 95, 1162–1172 (2014).PubMed 

Google Scholar 
Butterbach-Bahl, K., Baggs, E. M., Dannenmann, M., Kiese, R. & Zechmeister-Boltenstern, S. Nitrous oxide emissions from soils: how well do we understand the processes and their controls? Philos. Trans. R. Soc. B: Biol. Sci. 368, 20130122 (2013).
Google Scholar 
Baggs, E. M., Rees, R. M., Smith, K. A. & Vinten, A. J. A. Nitrous oxide emission from soils after incorporating crop residues. Soil Use & Manag. 16, 82–87 (2000).
Google Scholar 
Le Roux, X., Bardy, M., Loiseau, P. & Louault, F. Stimulation of soil nitrification and denitrification by grazing in grasslands: do changes in plant species composition matter? Oecologia 137, 417–425 (2003).ADS 
PubMed 

Google Scholar 
Morley, N. & Baggs, E. M. Carbon and oxygen controls on N2O and N2 production during nitrate reduction. Soil Biol. Biochem. 42, 1864–1871 (2010).CAS 

Google Scholar 
Davidson, E. A. & Kanter, D. Inventories and scenarios of nitrous oxide emissions. Environ. Res. Lett. 9, 105012 (2014).ADS 

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

Google Scholar 
Lehmann, J., Bossio, D. A., Kögel-Knabner, I. & Rillig, M. C. The concept and future prospects of soil health. Nat. Rev. Earth Environ. 1, 544–553 (2020).ADS 
PubMed 
PubMed Central 

Google Scholar 
Knapp, A. K. et al. Pushing precipitation to the extremes in distributed experiments: recommendations for simulating wet and dry years. Glob. Change Biol. 23, 1774–1782 (2017).ADS 

Google Scholar 
Cole, A. J. et al. Grassland biodiversity restoration increases resistance of carbon fluxes to drought. J. Appl. Ecol. 56, 1806–1816 (2019).CAS 

Google Scholar 
Fuchslueger, L., Bahn, M., Fritz, K., Hasibeder, R. & Richter, A. Experimental drought reduces the transfer of recently fixed plant carbon to soil microbes and alters the bacterial community composition in a mountain meadow. N. Phytologist 201, 916–927 (2014).CAS 

Google Scholar 
Buyer, J. S. & Sasser, M. High throughput phospholipid fatty acid analysis of soils. Appl. Soil Ecol. 61, 127–130 (2012).
Google Scholar 
Frostegård, Å., Bååth, E. & Tunlio, A. Shifts in the structure of soil microbial communities in limed forests as revealed by phospholipid fatty acid analysis. Soil Biol. Biochem. 25, 723–730 (1993).
Google Scholar 
Olsson, P. A., Thingstrup, I., Jakobsen, I. & Bååth, E. Estimation of the biomass of arbuscular mycorrhizal fungi in a linseed field. Soil Biol. Biochem. 31, 1879–1887 (1999).CAS 

Google Scholar 
Hopkin, S. P. A key to the Collembola (springtails) of Britain and Ireland (FSC, 2007).Krantz, G. W. & Walter, D. E. A manual of acarology (Texas Tech Universty Press, 2009).Caruso, T. & Migliorini, M. Euclidean geometry explains why lengths allow precise body mass estimates in terrestrial invertebrates: The case of oribatid mites. J. Theor. Biol. 256, 436–440 (2009).ADS 
MathSciNet 
CAS 
PubMed 
MATH 

Google Scholar 
Ganihar, S. R. Biomass estimates of terrestrial arthropods based on body length. J. Biosci. 22, 219–224 (1997).
Google Scholar 
Johnson, D., Vachon, J., Britton, A. J. & Helliwell, R. C. Drought alters carbon fluxes in alpine snowbed ecosystems through contrasting impacts on graminoids and forbs. N. Phytologist 190, 740–749 (2011).CAS 

Google Scholar 
Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280 (2001).ADS 
PubMed 

Google Scholar 
Anderson, M. J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46 (2001).
Google Scholar 
Anderson, M. J. Distance-Based Tests for Homogeneity of Multivariate Dispersions. Biometrics 62, 245–253 (2006).MathSciNet 
PubMed 
MATH 

Google Scholar 
Zuur, A., Ieno, E., Walker, N., Saveliev, A. & Smith, G. Mixed effects models and extensions in ecology with R. (Springer, 2009). More

Barbier, E. B. Marine ecosystem services. Curr. Biol. 27, R507–R510 (2017).CAS 
PubMed 

Google Scholar 
Liquete, C., Piroddi, C., Macías, D., Druon, J.-N. & Zulian, G. Ecosystem services sustainability in the Mediterranean Sea: Assessment of status and trends using multiple modelling approaches. Sci. Rep. 6, 34162 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Halpern, B. S. et al. Recent pace of change in human impact on the world’s ocean. Sci. Rep. 9, 1–8 (2019).CAS 

Google Scholar 
Duarte, C. M. et al. Rebuilding marine life. Nature 580, 39–51 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Long, R. D., Charles, A. & Stephenson, R. L. Key principles of marine ecosystem-based management. Mar. Policy 57, 53–60 (2015).
Google Scholar 
Link, J. S. & Browman, H. I. Operationalizing and implementing ecosystem-based management. ICES J. Mar. Sci. 74, 379–381 (2017).
Google Scholar 
EC. A Farm to Fork Strategy for a Fair, Healthy and Environmentally-Friendly Food System. Brussels: European Commission. (2020).EC. Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions, “Pathway to a healthy planet for all” with the sub-title “EU action Plan: ’Towards zero pollution for air, water and soil, COM (2021) 400. (2021).EC. The European Green Deal. Communication from the Commission to the European Parliament, the European Council, the European Economic and Social Committee and the Committee of the Regions, COM (2019) 640. (2019).Alexander, K. & Haward, M. The human side of marine ecosystem-based management (EBM): ‘Sectoral interplay’ as a challenge to implementing EBM. Mar. Policy 101, 33–38 (2019).
Google Scholar 
EC. The EU Blue Economy Report 2021. (2021).Ostlaender, N. et al. Modelling Inventory and Knowledge Man-agement System of the European Commission (MIDAS) (Publications Office of the European Union, 2019).
Google Scholar 
Friedland, R. et al. Effects of nutrient management scenarios on marine eutrophication indicators: A Pan-European, multi-model assessment in support of the Marine Strategy Framework Directive. Front. Mar. Sci. 8, 596126 (2021).
Google Scholar 
Piroddi, C. et al. Effects of nutrient management scenarios on marine food webs: A Pan-European Assessment in support of the Marine Strategy Framework Directive. Front. Mar. Sci. 8, 179 (2021).
Google Scholar 
Corrales, X. et al. Multi-zone marine protected areas: Assessment of ecosystem and fisheries benefits using multiple ecosystem models. Ocean Coast. Manag. 193, 105232 (2020).
Google Scholar 
Bentley, J. W. et al. Refining fisheries advice with stock-specific ecosystem information. Front. Mar. Sci. 8, 602072 (2021).
Google Scholar 
Steenbeek, J. et al. Making spatial-temporal marine ecosystem modelling better—A perspective. Environ. Model. Softw. 145, 105209 (2021).PubMed 
PubMed Central 

Google Scholar 
Heymans, J. J. et al. The ocean decade: A true ecosystem modelling challenge. Front. Mar. Sci. 7, 554573 (2020).
Google Scholar 
Hernvann, P.-Y. et al. The Celtic sea through time and space: Ecosystem modeling to unravel fishing and climate change impacts on food-web structure and dynamics. Front. Mar. Sci. 7, 1018 (2020).
Google Scholar 
Christensen, V. & Walters, C. J. Ecopath with Ecosim: Methods, capabilities and limitations. Ecol. Model. 172, 109–139 (2004).
Google Scholar 
Steenbeek, J. et al. Bridging the gap between ecosystem modeling tools and geographic information systems: Driving a food web model with external spatial–temporal data. Ecol. Model. 263, 139–151 (2013).
Google Scholar 
Christensen, V. et al. Representing variable habitat quality in a spatial food web model. Ecosystems 17, 1397–1412 (2014).CAS 

Google Scholar 
de Mutsert, K., Lewis, K., Milroy, S., Buszowski, J. & Steenbeek, J. Using ecosystem modeling to evaluate trade-offs in coastal management: Effects of large-scale river diversions on fish and fisheries. Ecol. Model. 360, 14–26 (2017).
Google Scholar 
Serpetti, N. et al. Modelling small scale impacts of Multi-Purpose Platforms: An ecosystem approach. Front. Mar. Sci. 8, 778 (2021).
Google Scholar 
DFO. Technical review of Roberts Bank Terminal 2 environmental assessment: section 10.3—assessing ecosystem productivity. DFO Can. Sci. Advis. Sec. Sci. Resp. 2016/050 (2016).Coll, M., Pennino, M. G., Steenbeek, J., Solé, J. & Bellido, J. M. Predicting marine species distributions: Complementarity of food-web and Bayesian hierarchical modelling approaches. Ecol. Model. 405, 86–101 (2019).
Google Scholar 
Coll, M. et al. The biodiversity of the Mediterranean Sea: Estimates, patterns, and threats. PLoS One 5, e11842 (2010).ADS 
PubMed 
PubMed Central 

Google Scholar 
Coll, M. et al. The Mediterranean Sea under siege: Spatial overlap between marine biodiversity, cumulative threats and marine reserves. Glob. Ecol. Biogeogr. 21, 465–480 (2012).
Google Scholar 
Micheli, F. et al. Cumulative human impacts on mediterranean and black sea marine ecosystems: Assessing current pressures and opportunities. PLoS One 8, e79889 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Piroddi, C., Colloca, F. & Tsikliras, A. C. The living marine resources in the Mediterranean Sea large marine ecosystem. Environ. Dev. 36, 100555 (2020).PubMed 
PubMed Central 

Google Scholar 
Barale, V. & Gade, M. Remote Sensing of the European Seas. (Springer, 2008).Siokou-Frangou, I. et al. Plankton in the open Mediterranean Sea: A review. Biogeosciences 7, 1543–1586 (2010).ADS 

Google Scholar 
Spalding, M. D. et al. Marine ecoregions of the world: A bioregionalization of coastal and shelf areas. Bioscience 57, 573–583 (2007).
Google Scholar 
Bianchi, C. N. et al. In Life in the Mediterranean Sea: A Look at Habitat Changes, vol. 1 55 (2012).Danovaro, R. et al. Deep-sea biodiversity in the Mediterranean Sea: The known, the unknown, and the unknowable. PLoS One 5, e11832 (2010).ADS 
PubMed 
PubMed Central 

Google Scholar 
Moullec, F. et al. Capturing the big picture of Mediterranean marine biodiversity with an end-to-end model of climate and fishing impacts. Prog. Oceanogr. 178, 102179 (2019).
Google Scholar 
Macias, D., Garcia-Gorriz, E., Piroddi, C. & Stips, A. Biogeochemical control of marine productivity in the Mediterranean Sea during the last 50 years. Glob. Biogeochem. Cycles 28, 897–907 (2014).ADS 
CAS 

Google Scholar 
Piroddi, C. et al. Historical changes of the Mediterranean Sea ecosystem: Modelling the role and impact of primary productivity and fisheries changes over time. Sci. Rep. 7, 44491 (2017).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lotze, H. K. et al. Depletion, degradation, and recovery potential of estuaries and coastal seas. Science 312, 1806–1809 (2006).ADS 
CAS 
PubMed 

Google Scholar 
Lotze, H. K., Coll, M. & Dunne, J. A. Historical changes in marine resources, food-web structure and ecosystem functioning in the Adriatic Sea, Mediterranean. Ecosystems 14, 198–222 (2011).
Google Scholar 
Macias, D., Huertas, I. E., Garcia-Gorriz, E. & Stips, A. Non-Redfieldian dynamics driven by phytoplankton phosphate frugality explain nutrient and chlorophyll patterns in model simulations for the Mediterranean Sea. Prog. Oceanogr. 173, 37–50 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Spedicato, M. T. et al. The MEDITS trawl survey specifications in an ecosystem approach to fishery management. Sci. Mar. 83, 9–20 (2019).
Google Scholar 
Corrales, X. et al. Hindcasting the dynamics of an Eastern Mediterranean marine ecosystem under the impacts of multiple stressors. Mar. Ecol. Prog. Ser. 580, 17–36 (2017).ADS 

Google Scholar 
FAO. The State of the Mediterranean and Black Sea fisheries (General Fisheries Commission for the Mediterranean (GFCM), 2020).Ferrà, C. et al. Mapping change in bottom trawling activity in the Mediterranean Sea through AIS data. Mar. Policy 94, 275–281 (2018).
Google Scholar 
Russo, T. et al. Trends in effort and yield of trawl fisheries: A case study from the Mediterranean Sea. Front. Mar. Sci. 6, 153 (2019).ADS 

Google Scholar 
Ramírez, F., Coll, M., Navarro, J., Bustamante, J. & Green, A. J. Spatial congruence between multiple stressors in the Mediterranean Sea may reduce its resilience to climate impacts. Sci. Rep. 8, 1–8 (2018).
Google Scholar 
Coll, M., Steenbeek, J., Ben Rais Lasram, F., Mouillot, D. & Cury, P. ‘Low-hanging fruit’ for conservation of marine vertebrate species at risk in the Mediterranean Sea. Glob. Ecol. Biogeogr. 24, 226–239 (2015).
Google Scholar 
Ruiz, J. et al. “Strengthening regional cooperation in the area of large pelagic fishery data collection (RECOLAPE)”, Annex III “Biological data collection for fisheries on highly migratory species” (2019).Boerder, K., Schiller, L. & Worm, B. Not all who wander are lost: Improving spatial protection for large pelagic fishes. Mar. Policy 105, 80–90 (2019).
Google Scholar 
Giakoumi, S. et al. Conserving European biodiversity across realms. Conserv. Lett. 12, e12586 (2019).
Google Scholar 
Gascuel, D. & Cheung, W. W. In Predicting Future Oceans 79–85 (Elsevier, 2019).Macias, D., Garcia-Gorriz, E. & Stips, A. Major fertilization sources and mechanisms for Mediterranean Sea coastal ecosystems. Limnol. Oceanogr. 63, 897–914 (2018).ADS 
CAS 

Google Scholar 
Alvarez-Berastegui, D., Tugores, M., Ottmann, D., Martín-Quetglas, M. & Reglero, P. Bluefin tuna larval indices in the Western Mediterranean, ecological and analytical sources of uncertainty. Collect. Vol. Sci. Pap. ICCAT. 77, 289–311 (2020).
Google Scholar 
ICCAT. 2020 SCRS Advice to the Commission (Madrid, Spain, 2020).Clavel-Henry, M., Piroddi, C., Quattrocchi, F., Macias, D. & Christensen, V. Spatial distribution and abundance of mesopelagic fish biomass in the Mediterranean Sea. Front. Mar. Sci. 7, 1136 (2020).
Google Scholar 
García-Ruiz, C. et al. Spatio-temporal patterns of macrourid fish species in the northern Mediterranean Sea. Sci. Mar. 83, 117–127 (2019).
Google Scholar 
Ainsworth, C. Quantifying species abundance trends in the Northern Gulf of California using local ecological knowledge. Mar. Coast. Fish. 3, 190–218 (2011).
Google Scholar 
Morris, E. K. et al. Choosing and using diversity indices: Insights for ecological applications from the German Biodiversity Exploratories. Ecol. Evol. 4, 3514–3524 (2014).PubMed 
PubMed Central 

Google Scholar 
Coll, M. et al. Ecological indicators to capture the effects of fishing on biodiversity and conservation status of marine ecosystems. Ecol. Ind. 60, 947–962 (2016).
Google Scholar 
Swartz, W., Sala, E., Tracey, S., Watson, R. & Pauly, D. The spatial expansion and ecological footprint of fisheries (1950 to present). PLoS One 5, e15143 (2010).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Damasio, L. M., Peninno, M. G. & Lopes, P. F. Small changes, big impacts: Geographic expansion in small-scale fisheries. Fish. Res. 226, 105533 (2020).
Google Scholar 
Coll, M. et al. Assessing fishing and marine biodiversity changes using fishers’ perceptions: The Spanish Mediterranean and Gulf of Cadiz case study. PLoS One 9, e85670 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Tsikliras, A. C., Dinouli, A., Tsiros, V.-Z. & Tsalkou, E. The Mediterranean and Black Sea fisheries at risk from overexploitation. PLoS ONE 10, e0121188 (2015).PubMed 
PubMed Central 

Google Scholar 
Pittman, S. et al. Seascape ecology: Identifying research priorities for an emerging ocean sustainability science. Mar. Ecol. Prog. Ser. 663, 1–29 (2021).ADS 

Google Scholar 
Kritzer, J. P. & Liu, O. R. In Stock Identification Methods 29–57 (Elsevier, 2014).Piroddi, C., Heymans, J. J., Macias, D., Gregoire, M. & Townsend, H. Editorial: Using ecological models to support and shape environmental policy decisions. Front. Mar. Sci. https://doi.org/10.3389/fmars.2021.815313 (2021).
Google Scholar 
Macias, D. et al. JRC marine modelling framework in support of the marine strategy framework directive: Inventory of models, basin configurations and datasets. Update 2018. (2018).Piante, C. & Ody, D. Blue Growth in the Mediterranean Sea: The Challenge of Good Environmental Status. 192 (France, 2015).Borja, A. et al. Past and future grand challenges in marine ecosystem ecology. Front. Mar. Sci. 7, 362 (2020).
Google Scholar 
Claudet, J., Loiseau, C., Sostres, M. & Zupan, M. Underprotected marine protected areas in a global biodiversity hotspot. One Earth 2, 380–384 (2020).ADS 

Google Scholar 
Piroddi, C., Coll, M., Steenbeek, J., Moy, D. M. & Christensen, V. Modelling the Mediterranean marine ecosystem as a whole: Addressing the challenge of complexity. Mar. Ecol. Prog. Ser. 533, 47–65 (2015).ADS 

Google Scholar 
Walters, C., Pauly, D. & Christensen, V. Ecospace: Prediction of mesoscale spatial patterns in trophic relationships of exploited ecosystems, with emphasis on the impacts of marine protected areas. Ecosystems 2, 539–554 (1999).
Google Scholar 
Christensen, V., Walters, C., Pauly, D. & Forrest, R. Ecopath with Ecosim 6: A User’s Guide (University of British Columbia, 2008).
Google Scholar 
Kaschner, K. et al. AquaMaps: Predicted range maps for aquatic species. In World Wide Web Electronic Publication, www.aquamaps.org, Version, vol. 8, 2016 (2016).De Mutsert, K., Lewis, K. A., White, E. D. & Buszowski, J. End-to-End modeling reveals species-specific effects of large-scale coastal restoration on living resources facing climate change. Front. Mar. Sci. 8, 104 (2021).
Google Scholar 
Coll, M. et al. Advancing global ecological modeling capabilities to simulate future trajectories of change in marine ecosystems. Front. Mar. Sci. 741, 567877 (2020).
Google Scholar 
Shannon, C. & Weaver, W. (Univ. Illinois Press, 1949).Ainsworth, C. H. & Pitcher, T. J. Modifying Kempton’s species diversity index for use with ecosystem simulation models. Ecol. Ind. 6, 623–630 (2006).
Google Scholar 
Coll, M. & Steenbeek, J. Standardized ecological indicators to assess aquatic food webs: The ECOIND software plug-in for Ecopath with Ecosim models. Environ. Model. Softw. 89, 120–130 (2017).
Google Scholar 
Taconet, M., Kroodsma, D. & Fernandes, J. Global Atlas of AIS-Based Fishing Activity—Challenges and Opportunities (2021). More

We investigated whether ricefield mosquito larval control and/or rice cultivation practices are associated with malaria vector densities through a systematic review and meta-analysis. Forty-seven experimental studies were eligible for inclusion in the qualitative analysis and thirty-three studies were eligible for the meta-analysis. It was demonstrated that the use of fish, chemical and biological larvicides in rice fields were effective in controlling larval malaria vector densities at all developmental stages. Intermittent irrigation, however, could only significantly reduce late-stage larvae. Based on a limited number of studies, meta-analyses on other forms of larval control such as monomolecular surface films (MSFs), neem, copepods and Azolla failed to demonstrate any consistent reduction in anopheline numbers. Similarly, rice cultivation practices such as plant variety and density, type of levelling and pesticide application were not generally associated with reduced malaria vectors. Nonetheless, in one study, minimal tillage was observed to reduce average numbers of larvae throughout a cropping season. In another study, herbicide application increased larval abundance over a 4-week period, as did one-time drainage in a third study.
Despite their different modes of action, the use of chemical and bacterial larvicides and MSFs were all relatively effective measures of larval control in rice fields, varying between a 57% to 76% reduction in vector abundance compared to no larviciding. Their effects were highest (often reaching 100% reduction) only shortly following application but did not persist for longer than two weeks. These larvicides mostly had short residual half-lives because they were applied to paddy water which was naturally not completely stagnant: there was a small but constant process of water loss (through drainage, evapotranspiration and percolation) and replacement through irrigation. Hence, even with a residual formulation, weekly re-application would be needed for sustained control47,40,41,50. This would be very labour- and cost-intensive to scale-up, to ensure that larvicides are evenly distributed across vast areas (even at plot/sub-plot level) throughout at least one 5-month long rice-growing season per year42,51. Aerial application (including unmanned aerial vehicles), although widely used in the US and Europe, is unlikely to be a feasible delivery system for smallholders in SSA, even in large irrigation schemes26,27,48,49. Furthermore, if synthetic organic chemicals were to be considered for riceland malaria vector control, their management in the current landscape of insecticide resistance across Africa must be considered.Biological control using fish was found to be, in general, slightly more effective than (chemical, bacterial and MSF) larviciding. The degree of effectiveness was dependent on the fish species and their feeding preferences: surface-feeding, larvivorous species provided better anopheline control than bottom-feeding selective feeders4,43. Selecting the most suitable fish for local rice fields is not straightforward; many criteria need to be considered4,52,53. Generally, fish were well-received by rice farmers, perceived to contribute to increased yield by reducing weeds and pests and providing fertiliser through excrement43,44. This was reportedly also observed in Guangxi, China, where a certain proportion of the field had to be deepened into a side-trench where the fish could take shelter when the fields were drained. Even with this reduction in rice production area, carp rearing still increased yields by 10% and farmer’s income per hectare by 70%53. Unfortunately, none of the eligible studies in this review had included yield or water use as an outcome. Future entomological studies need to measure these critical agronomic variables so that studies of vector control in rice can be understood by, and transferred to, agronomists. In SSA, irrigated rice-fish farming can be scaled up provided that an inventory of fish species suitable for specific locations is available and that water is consistently available in fields (an important limiting factor in African irrigation schemes)54. Lessons can be learnt from successful large-scale rice-fish systems in Asia, where they have served as win–win solutions for sustainable food production and malaria control16,55.Overall, there was only limited evidence that intermittent irrigation is effective at reducing late-instar anopheline larvae in rice fields. This finding contrasts with prior reviews, which found mixed results (regardless of larval stage) but emphasised that success was site-specific4,17,56. This contrast is presumably due to the inclusion criteria of our systematic review. These reviews excluded studies in various geographical settings and some older studies that reported successful anopheline control with intermittent irrigation but lacked either a contemporaneous control arm, adequate replication or adequate differentiation between culicines and anophelines16,57,50,51,52,61. It seems, from our review, that intermittent irrigation does not prevent the recruitment of early instars (and in one case, may have encouraged oviposition31) but tends to prevent their development into late-stage immatures. This important conclusion is, however, based only on four studies; more evidence is urgently needed where future trials should consider the basic principles of modern trials with adequate replication, controls and differentiation between larval instars and species.Generally, it is observed that drainage, passive or active, did not reliably reduce overall numbers of mosquito immatures. In India and Kenya, closer inspection revealed that soils were not drying sufficiently, so any stranded larvae were not killed31,46. Highlighted by van der Hoek et al.29 and Keiser et al.17, water management in rice fields is very dependent on the physical characteristics of the soil and the climate and is most suited to places that not only favour rapid drying, but also have a good control of water supply17,56. Moreover, repeated drainage, although directed against mosquitoes, can also kill their aquatic predators62. Since mosquitoes can re-establish themselves in a newly flooded rice field more quickly than their predators, intermittent irrigation with more than a week between successive drying periods can permit repeated cycles of mosquito breeding without any predation pressure. Its efficacy against malaria vectors is therefore highly reliant on the timing of the wetting and drying periods. Further site-specific research on timing, especially with regards to predator–prey interactions within the rice agroecosystem, is required to find the perfect balance.Another limitation in intermittent irrigation is that it cannot be applied during the first two to three weeks following transplanting, because rice plants must remain flooded to recover from transplanting shock. Unfortunately, this time coincides with peak vector breeding. Thus, other methods of larval control would be required to fill this gap. To agronomists, intermittent irrigation provides benefits to farmers, as it does not penalise yield but significantly reduces water consumption. Nonetheless, farmer compliance seems to be variable, especially in areas where water availability is inconsistent and intermittent irrigation would potentially require more labour31,32,39. Importantly, rice farmers doubted their ability to coordinate water distribution evenly amongst themselves, suggesting that there may be sharing issues, as in the “tragedy of the commons”63. Instead, they said that they preferred to have an agreed authority to regulate water46.No general conclusions could be made on the effect on malaria vectors of other rice cultivation practices (apart from water management) because only one study was eligible for each practice. Nevertheless, these experiments on pesticide application, tillage and weed control, as well as another study on plant spacing (not eligible since glass rods were used to simulate rice plants), do illustrate that small changes in agronomic inputs and conditions can have considerable effects on mosquito densities, not just rice yield36,38,64. Moreover, in partially- or shallowly-flooded plots, the larvae are often concentrated in depressions (usually footprints), suggesting that rice operations which leave or remove footprints (e.g. hand-weeding, drum seeders, levelling) will influence vector breeding4.Our study has some important limitations. First, in most trials, the units of intervention were replicate plots of rice, and success was measured as a reduction in larval densities within treated plots. This design focuses on the identification of effective and easy-to-implement ways of growing rice without growing mosquitoes, on the assumption that higher vector densities are harmful. However, from a public health perspective, the need for epidemiological outcomes is often, and reasonably, stressed22,65. Nonetheless, from a farmers’ perspective, it is also important to consider whether the vectors emerging from their rice fields significantly contribute to the local burden of malaria and to determine how this contribution can be minimised. There is evidence that riceland vectors do increase malaria transmission, since human biting rates are much higher in communities living next to rice schemes than their non-rice counterparts66 and that additional riceland vectors may intensify transmission and malaria prevalence in rice communities15. Hence, when investigating how rice-attributed malaria risk can be minimised, mosquito abundance as measured in the experimental rice trials is a useful indicator of potential impact on epidemiological outcomes.Second, larval density was not always separated into larval developmental stages. This can be misleading because some interventions work by reducing larval survival (but not by preventing oviposition) and development to late instars and pupae. Therefore, an intervention could completely eliminate late-stage larvae and pupae but have little effect on the total number of immatures. This was illustrated in our meta-analyses of intermittent irrigation in Table 3 and Supplementary Table 5, and could have been the case for some studies that failed to demonstrate consistent reductions in overall anopheline numbers but did not differentiate between larval instars34,45,67,60,69. We infer that when monitoring mosquito immatures in rice trials, it is important to distinguish between larval instars and pupae. Pupae should always be counted separately since its abundance is the most direct indicator of adult productivity70.Third, experimental trials rarely reported the timing of intervention application or accounted for different rice-growing phases, or “days after transplantation”, in the outcome. Both aspects are important to consider since an intervention may be suited to control larvae during certain growth phases but not others. This is illustrated by Djegbe et al.38, where, compared to deep tillage, minimal tillage could significantly reduce larvae during the early stages of rice cultivation but not during tillering and maturation38. In contrast, other interventions, such as Azolla and predatory copepods, took time to grow and accumulate, and were more effective during the later stages of a rice season45,67,71. This differentiation is important because it can identify components that could potentially form a complementary set of interventions against riceland malaria vectors, each component being effective at different parts of the season. Since rice fields, and hence the dynamics of riceland mosquito populations, vary from place to place, this set of interventions must also be robust. Special attention must be paid to the early stages of rice cultivation, particularly the first few weeks after transplanting (or sowing), since, with many vector species, a large proportion of adult mosquitoes are produced during this time.Fourth, the analysis of entomological counts is often inadequate. Many studies failed to provide the standard deviation (or any other measure of error) for larval counts and could not be included in the quantitative analysis. Often, due to the extreme (and not unexpected) variability of larval numbers, sample sizes were insufficient to calculate statistically significant differences between treatments. Fifth, a high risk of bias was found across both CTS and CITS studies, including high heterogeneity and some publication bias. Study quality was, in general, a shortcoming and limited the number of eligible studies for certain interventions, including intermittent irrigation. Moreover, there are conspicuous a priori reasons for bias in such experimental trials: trial locations are frequently chosen to maximise the probability of success.Finally, few studies were conducted in African countries, where the relationship between rice and malaria is most important because of the efficiency, and the “rice-philic” nature, of the vector An. gambiae s.l.15. In particular, there was a lack of studies on the effectiveness and scalability of biological control and rice cultivation practices. There is also very little information (particularly social science studies) on the views and perspectives of African rice farmers on mosquitoes in rice and interventions to control them72,73.In the future, as malaria declines (particularly across SSA), the contribution of rice production to increased malaria transmission is likely to become more conspicuous15. Unless this problem is addressed, rice growing will probably become an obstacle to malaria elimination. Current default methods of rice production provide near-perfect conditions for the larvae of African malaria vectors. Therefore, we need to develop modified rice-growing methods that are unfavourable to mosquitoes but still favourable for the rice. Although larviciding and biological control may be appropriate, their unsustainable costs remain the biggest barrier to uptake amongst smallholder farmers. Future investigations into riceland vector control should pay more attention to interventions that may be useful to farmers.Supported by medical entomologists, agronomists should lead the research task of identifying cultivation methods that achieve high rice productivity whilst suppressing vector productivity. Rice fields are a major global source of greenhouse gases, and agronomists have responded by successfully developing novel cultivation methods that minimise these emissions while maintaining yield. We need the same kind of response from agronomists, to achieve malaria control co-benefits within rice cultivation. At present, only a few aspects of rice cultivation have been investigated for their effects on mosquitoes, and the potential of many other practices for reducing anopheline numbers are awaiting study. Due to the spatial and temporal heterogeneity of rice agroecosystems, it is likely that no single control method can reduce mosquito numbers throughout an entire cropping season and in all soil types and irrigation methods. Thus, effective overall control is likely to come from a combination of local, site-specific set of complementary methods, each of which is active and effective during a different phase of the rice-growing season. More

Using a robust recent dataset, our analyses show that cocoa production is significantly affected by the maximum magnitude of ENSO phase during the current and previous purchase years (Fig. 2). The instantaneous effect is negative, followed by delayed positive effects in the two following years and negative in the third following year, combining to give a picture of multi-year fluctuations in cocoa production as a result of El Niño/La Niña events. Using a 70-year dataset, we show significant changes in these instantaneous and delayed ENSO-production relationships between recent and past time periods (Fig. 3). Using ERA5 data for the cocoa production area of Ghana, summarised at the same temporal resolution as the production data, we demonstrate significant relationships between ENSO phase and climate, with significant changes in mean climate and in ENSO-climate relationships (Fig. 4) between recent and past time periods. This agrees with prior work suggesting that ENSO may impact West Africa5,15, despite no current evidence of teleconnections between ENSO phase and West African climate17.Our 70-year production dataset represents a temporal extent unmatched by other research, however was aggregated to fewer replicates than the 21-year analysis (6 regions vs 68 districts). While this may represent reduced power, results from the overlapping time period of the two datasets strongly agree. The computation of yield, a more comparable metric between different-sized areas than total production, was not possible because data on area under production (AUP) were not available. However, the detrending process employed successfully eliminated variation between districts or regions (of which AUP is likely a substantial component) and long-term technological trends that would otherwise confound our ability to isolate the ENSO signal (Supplementary results).Perennial crops have multi-year growing patterns, with allocation of resources to growth, development and reproduction driven by climate in ways that are not fully understood29. ENSO generally peaks between October and December, also the busiest cocoa purchase period: thus we observe a relatively instantaneous apparent effect of ENSO phase on cocoa production. This reduction in cocoa production under El Niño inis consistent with results from farm monitoring8 and large-scale farm surveys30 evidencing production declines in from other regions (where teleconnections are understood), and with analyses of production data from West Africa31. During the main cocoa purchase period, coinciding with the minor wet and major dry seasons, we observe increases in water deficit during El Niño, leading to drought stress conditions. In small-scale cocoa studies, drought stress is correlated with reduction in pod production and increased tree mortality8,32, and in similar studies of other tree crops drought is directly linked to reduction in fruit or nut production33, although in all cases the mechanisms are unclear. Drought may generally create unfavourable conditions for growth and reproduction through reduced availability of water for vital processes, or more specifically by promoting disease incidence and pod rot8, increasing the chance of fire, increasing competition for soil moisture32, and/or reducing pollinator populations34. Alternatively, cocoa may respond to reduced water availability by reallocation of resources away from energetically expensive reproduction: rainfall exclusion experiments suggest that in the medium term, while bean production drops, vegetative growth is not significantly reduced during drought32.The significant increases in mean temperature and average drought stress we observed in some seasons over time is such that the climate experienced during El Niño events in recent decades represent novel extreme conditions for Ghana’s cocoa agriculture. This causes significant changes in the responses of cocoa production to ENSO phase over the same time period. One explanation for this may be that the warm, dry El Niño conditions in Ghana in the past were within the environmental tolerance of cocoa, leading to allocation of resources to reproduction in response to drought, increasing cocoa bean production and resulting in less severe instantaneous and delayed responses to ENSO phase (Fig. 3a–d) However, in recent decades this level or greater drought stress has become the norm (Fig. 4i–l), with El Niño conditions apparently triggering a different response mode, allocating resources away from reproduction in the short term and creating oscillating resource allocation over the following years.However, understanding the delayed responses of cocoa is challenging, especially as these represent a novel finding. There is little research that explores multi-annual physiological or ecological responses of cocoa to drought, and the explanation is likely to be a combination of both residual/delayed climatic responses to ENSO phase, and of life history strategies. The observed increase in production during the two years following El Niño may be explained by post-drought reallocation of resources to reproduction as remediation for lost reproductive output in the instantaneous response, or a shift to a ‘faster’ strategy by allocating resources to reproduction over the longer term, becoming evident in the data in subsequent years. Alternatively, this may be explained by favourable climatic conditions occurring during an El Niño event that impact the following years’ crop. March and April is a crucial time for cocoa pod development in Ghana and in recent years El Niño appears to bring greater rainfall during these months. Given the 6–9 months development of cocoa beans, the effects of this increased rainfall and reduced water deficit on cocoa production will be seen in the delayed response. We see evidence of this in the climate-change driven reversal of March–April rainfall patterns: while in the past El Niño has consistently resulted in drought stress, this reversal provides a respite from drought, buffering trees from reduced rainfall during the major wet season and giving sufficient resources for improved production in the following year.The robustness of our results provide evidence that may aid development of resilience strategies for ENSO-driven cocoa production variation in Ghana, but we may also consider whether these results can be generalised to the production of cocoa and/or perennial tree crops globally. The climatic impact of ENSO observed in Ghana is broadly consistent with many regions of the tropics2, the instantaneous cocoa production responses to El Niño are consistent with findings in these regions, and so we may expect these regions to see a similar pattern of multi-annual cocoa production variation in response to ENSO phase. However, there is considerable variation in ENSO responses among and within other perennial tree crops in regions where climatic responses to ENSO are similar to Ghana. Oil palm yields have been negatively associated with ENSO phase in Malaysia9, as have olive yields in Morocco (delayed by a year)33. Conversely, apple yields have been positively associated with ENSO phase in China10, as have coffee yields in Brazil35; however, no effect at all is seen in coffee in India over a 35-year time series7. Most of these analyses considered only a single ENSO phase (usually El Niño), and most considered only instantaneous impacts. However, it is clear that most of these crops do respond to ENSO, and given the shared biology it is reasonable to assume that delayed effects of ENSO phase are likely and should be considered to understand the full picture of ENSO impacts on perennial tree crops.The larger body of research into ENSO impacts on annual crops includes many studies using long time series, reporting high heterogeneity in space and among crops11,36,37. However, there appears to be little examination of changes in the direction and magnitude of ENSO responses over time; thus our findings are timely and signal that further research is needed to examine how changing climates may force novel extreme climatic conditions and shift response patterns to ENSO phase. Given that perennial tree crops are generally cash crops, and the utility of these crops to farmers are to a greater or lesser extent mediated by market forces, there is a need for improved forecasting of yield in response to changing climate and ENSO patterns to withstand production fluctuations. The low perishability of many perennial tree crops means that with accurate forecasting, supply may be managed or even exploited to ensure consistency of income both for farmers and those whose livelihoods depend on related food manufacturing industries.Our approach to understanding the responses of a perennial tree crop to ENSO phase and anthropogenic climate change exploited existing global, national and subnational datasets for climate and production with appropriate spatial and temporal resolution. We use freely available geographic and climate data, and employ highly replicable methods: a simple pipeline of climate data aggregation and summary computation, coupled with standard detrending and straightforward analytical methods with a relatively small computational requirement. This “big data” approach to agriculture-climate research demonstrates a relatively straightforward framework for understanding responses of agricultural productivity to climate and identifying temporal changes in these relationships. While small-scale studies examine the mechanisms of climate impacts through the interacting effects of agricultural practices, abiotic conditions, disease incidence and multi-trophic interactions, large-scale studies across regions and over time scales encompassing many ENSO oscillations are required to understand the global picture of perennial tree crop production security. Combined with local context-specific studies on governance arrangements16, such approaches could be crucial for reducing future vulnerability of these industries to increasing volatility under anthropogenic climate change. The main barrier to this research is the availability of production data from state or commercial entities. More

The effects of extreme drought on soil biochemical propertiesAs shown in Fig. 1A, the range of SOC during the early, midterm and late extreme drought experiments, were 73.53–251.44 g kg−1, 54.75–256.16 g kg−1, and 66.37–282.16 g kg−1, respectively. Concomitantly, DOC was 171.85–323.74 mg kg−1, 158.15 – 504.62 mg kg−1, and 166.63–418.43 mg kg−1, MBC was 247.80 – 461.69 mg kg−1, 257.90–450.98 mg kg−1, and 264.10–458.15 mg kg−1, respectively (Fig. 1B, C). The variation ranges of soil TN were 3.50–16.60 g kg−1, 4.70–34.5 g kg−1, and 6.70–32.50 g kg−1, respectively (Fig. 1D). Similarly, the variation ranges of NH4+ were 5.96–12.03 g kg−1, 5.39–12.59 g kg−1, and 5.74–13.03 g kg−1, NO3− were 2.27–8.79 mg kg−1, 5.07–9.62 mg kg−1, and 5.09–9.52 mg kg−1, respectively (Fig. 1E, F). The changes of SOC and NH4+ with soil depth were significantly different in different extreme drought periods and decreased significantly with the increase of soil depth (Table 1, P  More