Philippa Kaur

McCallum, M. L. Vertebrate biodiversity losses point to a sixth mass extinction. Biodivers. Conserv. 24, 2497–2519 (2015).
Google Scholar 
Wake, D. B. & Vredenburg, V. T. Are we in the midst of the sixth mass extinction? A view from the world of amphibians. Proc. Natl. Acad. Sci. 105, 11466–11473. https://doi.org/10.1073/pnas.0801921105 (2008).Article 
PubMed 
PubMed Central 

Google Scholar 
Blehert, D. S. et al. Bat white-nose syndrome: An emerging fungal pathogen?. Science 323, 227. https://doi.org/10.1126/science.1163874 (2009).CAS 
Article 
PubMed 

Google Scholar 
Pautasso, M., Aas, G., Queloz, V. & Holdenrieder, O. European ash (Fraxinus excelsior) dieback—A conservation biology challenge. Biol. Cons. 158, 37–49 (2013).
Google Scholar 
Daszak, P., Cunningham, A. A. & Hyatt, A. D. Infectious disease and amphibian population declines. Divers. Distrib. 9, 141–150 (2003).
Google Scholar 
Fisher, M. C., Gow, N. A. R. & Gurr, S. J. Tackling emerging fungal threats to animal health, food security and ecosystem resilience. Philos. Trans. R. Soc. B Biol. Sci. https://doi.org/10.1098/rstb.2016.0332 (2016).Article 

Google Scholar 
Fisher, M. C. et al. Emerging fungal threats to animal, plant and ecosystem health. Nature 484, 186–194 (2012).CAS 
PubMed 

Google Scholar 
Lips, K. R., Reeve, J. D. & Witters, L. R. Ecological traits predicting amphibian population declines in Central America. Conserv. Biol. 17, 1078–1088 (2003).
Google Scholar 
Zipkin, E. F., DiRenzo, G. V., Ray, J. M., Rossman, S. & Lips, K. R. Tropical snake diversity collapses after widespread amphibian loss. Science 367, 814–816. https://doi.org/10.1126/science.aay5733 (2020).CAS 
Article 
PubMed 

Google Scholar 
Berger, L. et al. Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America. Proc. Natl. Acad. Sci. 95, 9031–9036 (1998).CAS 
PubMed 
PubMed Central 

Google Scholar 
Martel, A. et al. Recent introduction of a chytrid fungus endangers Western Palearctic salamanders. Science 346, 630–631. https://doi.org/10.1126/science.1258268 (2014).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Yap, T. A., Koo, M. S., Ambrose, R. F., Wake, D. B. & Vredenburg, V. T. Averting a North American biodiversity crisis. Science 349, 481–482 (2015).CAS 
PubMed 

Google Scholar 
Weldon, C., du Preez, L. H., Hyatt, A. D., Muller, R. & Speare, R. Origin of the amphibian chytrid fungus. Emerg. Infect. Dis. 10, 2100–2105 (2004).PubMed 
PubMed Central 

Google Scholar 
Talley, B. L., Muletz, C. R., Vredenburg, V. T., Fleischer, R. C. & Lips, K. R. A century of Batrachochytrium dendrobatidis in Illinois amphibians (1888–1989). Biol. Cons. 182, 254–261 (2015).
Google Scholar 
Rodriguez, D., Becker, C., Pupin, N., Haddad, C. & Zamudio, K. Long-term endemism of two highly divergent lineages of the amphibian-killing fungus in the Atlantic Forest of Brazil. Mol. Ecol. 23, 774–787 (2014).CAS 
PubMed 

Google Scholar 
Goka, K. et al. Amphibian chytridiomycosis in Japan: Distribution, haplotypes and possible route of entry into Japan. Mol. Ecol. 18, 4757–4774 (2009).CAS 
PubMed 

Google Scholar 
Bataille, A. et al. Genetic evidence for a high diversity and wide distribution of endemic strains of the pathogenic chytrid fungus Batrachochytrium dendrobatidis in wild Asian amphibians. Mol. Ecol. 23, 4196–4209. https://doi.org/10.1111/mec.12385 (2013).CAS 
Article 

Google Scholar 
O’Hanlon, S. J. et al. Recent Asian origin of chytrid fungi causing global amphibian declines. Science 360, 621–627. https://doi.org/10.1126/science.aar1965 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Swei, A. et al. Is chytridiomycosis an emerging infectious disease in Asia?. PLoS ONE 6, e23179 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Bai, C. M., Garner, T. W. J. & Li, Y. M. First evidence of Batrachochytrium dendrobatidis in China: Discovery of chytridiomycosis in introduced American bullfrogs and native amphibians in the Yunnan Province, China. EcoHealth 7, 127–134. https://doi.org/10.1007/s10393-010-0307-0 (2010).Article 
PubMed 

Google Scholar 
Yang, H. et al. First detection of the amphibian chytrid fungus Batrachochytrium dendrobatidis in free-ranging populations of amphibians on mainland Asia: Survey in South Korea. Dis. Aquat. Org. 86, 9–13 (2009).
Google Scholar 
Fong, J. J. et al. Early 1900s detection of Batrachochytrium dendrobatidis in Korean amphibians. PLoS ONE 10, e0115656 (2015).PubMed 
PubMed Central 

Google Scholar 
Kusrini, M., Skerratt, L., Garland, S., Berger, L. & Endarwin, W. Chytridiomycosis in frogs of Mount Gede Pangrango, Indonesia. Diseases Aquat. Organ. 82, 187–194 (2008).CAS 

Google Scholar 
Laking, A. E., Ngo, H. N., Pasmans, F., Martel, A. & Nguyen, T. T. Batrachochytrium salamandrivorans is the predominant chytrid fungus in Vietnamese salamanders. Sci. Rep. 7, 44443. https://doi.org/10.1038/srep44443 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Zhu, W. et al. A survey for Batrachochytrium salamandrivorans in Chinese amphibians. Curr. Zool. 60, 729–735 (2014).
Google Scholar 
Beukema, W. et al. Environmental context and differences between native and invasive observed niches of Batrachochytrium salamandrivorans affect invasion risk assessments in the Western Palaearctic. Divers. Distrib. 24, 1788–1801. https://doi.org/10.1111/ddi.12795 (2018).Article 

Google Scholar 
Auliya, M. et al. The global amphibian trade flows through Europe: The need for enforcing and improving legislation. Biodivers. Conserv. https://doi.org/10.1007/s10531-016-1193-8 (2016).Article 

Google Scholar 
Scheffers, B. R., Edwards, D. P., Diesmos, A., Williams, S. E. & Evans, T. A. Microhabitats reduce animal’s exposure to climate extremes. Glob. Change Biol. 20, 495–503 (2014).
Google Scholar 
Schmeller, D. S. et al. People, pollution and pathogens—Global change impacts in mountain freshwater ecosystems. Sci. Total Environ. 622–623, 756–763. https://doi.org/10.1016/j.scitotenv.2017.12.006 (2018).CAS 
Article 
PubMed 

Google Scholar 
Bernardo-Cravo, A., Schmeller, D. S., Chatzinotas, A., Vredenburg, V. T. & Loyau, A. Environmental factors and host microbiomes shape host-pathogen dynamics. Trends Parasitol. 36, 29–36 (2020).
Google Scholar 
Harris, R. N. et al. Skin microbes on frogs prevent morbidity and mortality caused by a lethal skin fungus. ISME J. 3, 818–824. https://doi.org/10.1038/ismej.2009.27 (2009).CAS 
Article 
PubMed 

Google Scholar 
Harris, R. N., James, T. Y., Lauer, A., Simon, M. A. & Patel, A. Amphibian pathogen Batrachochytrium dendrobatidis is inhibited by the cutaneous bacteria of amphibian species. EcoHealth 3, 53–56. https://doi.org/10.1007/s10393-10005-10009-10391 (2006).Article 

Google Scholar 
Piovia-Scott, J. et al. Greater species richness of bacterial skin symbionts better suppresses the amphibian fungal pathogen Batrachochytrium dendrobatidis. Microb. Ecol. 74, 217–226 (2017).PubMed 

Google Scholar 
Ellison, S., Knapp, R. A., Sparagon, W., Swei, A. & Vredenburg, V. T. Reduced skin bacterial diversity correlates with increased pathogen infection intensity in an endangered amphibian host. Mol. Ecol. 28, 127–140 (2019).PubMed 

Google Scholar 
Jani, A. J. & Briggs, C. J. The pathogen Batrachochytrium dendrobatidis disturbs the frog skin microbiome during a natural epidemic and experimental infection. Proc. Natl. Acad. Sci. USA 111, E5049-5058. https://doi.org/10.1073/pnas.1412752111 (2014).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Kueneman, J. G. et al. The amphibian skin-associated microbiome across species, space and life history stages. Mol. Ecol. 23, 1238–1250 (2014).PubMed 

Google Scholar 
Kueneman, J. G. Ecology of the Amphibian Skin-Associated Microbiome and Its Role in Pathogen Defense (University of Colorado at Boulder, 2015).
Google Scholar 
Kueneman, J. G. et al. Community richness of amphibian skin bacteria correlates with bioclimate at the global scale. Nat. Ecol. Evolut. 3, 381–389. https://doi.org/10.1038/s41559-019-0798-1 (2019).Article 

Google Scholar 
Jiménez, R. R. & Sommer, S. The amphibian microbiome: Natural range of variation, pathogenic dysbiosis, and role in conservation. Biodivers. Conserv. 26, 763–786. https://doi.org/10.1007/s10531-016-1272-x (2017).Article 

Google Scholar 
Walke, J. B. et al. Amphibian skin may select for rare environmental microbes. ISME J 8, 2207–2217. https://doi.org/10.1038/ismej.2014.77 (2014).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
McKenzie, V. J., Bowers, R. M., Fierer, N., Knight, R. & Lauber, C. L. Co-habiting amphibian species harbor unique skin bacterial communities in wild populations. ISME J 6, 588–596. https://doi.org/10.1038/ismej.2011.129 (2012).CAS 
Article 
PubMed 

Google Scholar 
Bates, K. A. et al. Amphibian chytridiomycosis outbreak dynamics are linked with host skin bacterial community structure. Nat. Commun. 9, 693. https://doi.org/10.1038/s41467-018-02967-w (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Ellison, S. et al. The influence of habitat and phylogeny on the skin microbiome of amphibians in Guatemala and Mexico. Microb. Ecol. 78, 257–267 (2019).PubMed 

Google Scholar 
Fisher, M. C., Pasmans, F. & Martel, A. Virulence and pathogenicity of chytrid fungi causing amphibian extinctions. Annu. Rev. Microbiol. https://doi.org/10.1146/annurev-micro-052621-124212 (2021).Article 
PubMed 

Google Scholar 
Haver, M. et al. The role of abiotic variables in an emerging global amphibian fungal disease in mountains. Sci. Total Environ. 815, 152735 (2021).PubMed 

Google Scholar 
Turner, A., Wassens, S., Heard, G. & Peters, A. Temperature as a driver of the pathogenicity and virulence of amphibian chytrid fungus Batrachochytrium dendrobatidis: A systematic review. J. Wildl. Dis. 57, 477–494 (2021).PubMed 

Google Scholar 
Woodhams, D., Alford, R., Briggs, C., Johnson, M. & Rollins-Smith, L. Life history trade-offs influence disease in changing climates: Strategies of an amphibian pathogen. Ecology 89, 1627–1639 (2008).PubMed 

Google Scholar 
Sonn, J. M., Berman, S. & Richards-Zawacki, C. L. The influence of temperature on chytridiomycosis in vivo. EcoHealth 14, 762–770. https://doi.org/10.1007/s10393-017-1269-2 (2017).Article 
PubMed 

Google Scholar 
Schmidt, B., Küpfer, E., Geiger, C., Wolf, S. & Schär, S. Elevated temperature clears chytrid fungus infections from tadpoles of the midwife toad, Alytes obstetricans. Amphibia-Reptilia 32, 276–280 (2011).
Google Scholar 
Bielby, J., Cooper, N., Cunningham, A. A., Garner, T. W. J. & Purvis, A. Predicting susceptibility to future declines in the world’s frogs. Conserv. Lett. 1, 82–90 (2008).
Google Scholar 
Gray, M. J., Miller, D. L. & Hoverman, J. T. Ecology and pathology of amphibian ranaviruses. Dis. Aquat. Org. 87, 243–266 (2009).
Google Scholar 
Murray, K., Skerratt, L., Speare, R. & McCallum, H. Impact and dynamics of disease in species threatened by the amphibian chytrid fungus, Batrachochytrium dendrobatidis. Conserv. Biol. 23, 1242–1252 (2009).PubMed 

Google Scholar 
Schmeller, D. S. et al. Microscopic aquatic predators strongly affect infection dynamics of a globally emerged pathogen. Curr. Biol. 24, 176–180. https://doi.org/10.1016/j.cub.2013.11.032 (2014).CAS 
Article 
PubMed 

Google Scholar 
Metzger, M. J. et al. Environmental stratifications as the basis for national, European and global ecological monitoring. Ecol. Ind. 33, 26–35. https://doi.org/10.1016/j.ecolind.2012.11.009 (2013).Article 

Google Scholar 
Metzger, M. J. et al. A high-resolution bioclimate map of the world: A unifying framework for global biodiversity research and monitoring. Glob. Ecol. Biogeogr. 22, 630–638. https://doi.org/10.1111/geb.12022 (2013).Article 

Google Scholar 
Clare, F., Daniel, O., Garner, T. & Fisher, M. Assessing the ability of swab data to determine the true burden of infection for the amphibian pathogen Batrachochytrium dendrobatidis. EcoHealth 13, 360–367. https://doi.org/10.1007/s10393-016-1114-z (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Cheng, T. L., Rovito, S. M., Wake, D. B. & Vredenburg, V. T. Coincident mass extirpation of neotropical amphibians with the emergence of the infectious fungal pathogen Batrachochytrium dendrobatidis. Proc. Natl. Acad. Sci. 108, 9502–9507 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Vredenburg, V. T. et al. Pathogen invasion history elucidates contemporary host pathogen dynamics. PLoS ONE 14, e0219981. https://doi.org/10.1371/journal.pone.0219981 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Hyatt, A. D. et al. Diagnostic assays and sampling protocols for the detection of Batrachochytrium dendrobatidis. Dis. Aquat. Org. 73, 175–192 (2007).CAS 

Google Scholar 
Blooi, M. et al. Duplex real-time PCR for rapid simultaneous detection of Batrachochytrium dendrobatidis and B. salamandrivorans in amphibian samples. J. Clin. Microbiol. 51, 4173–4177 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Boyle, D. G., Boyle, D. B., Olsen, V., Morgan, J. A. T. & Hyatt, A. D. Rapid quantitative detection of chytridiomycosis (Batrachochytrium dendrobatidis) in amphibian samples using real-time Taqman PCR assay. Dis. Aquat. Org. 60, 141–148 (2004).CAS 

Google Scholar 
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12 (2011).
Google Scholar 
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Pruesse, E., Peplies, J. & Glöckner, F. O. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823–1829 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2012).PubMed 
PubMed Central 

Google Scholar 
Bokulich, N. A. & Mills, D. A. Improved selection of internal transcribed spacer-specific primers enables quantitative, ultra-high-throughput profiling of fungal communities. Appl. Environ. Microbiol. https://doi.org/10.1128/aem.03870-12 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
McMurdie, P. J. & Holmes, S. Waste not, want not: Why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10, e1003531 (2014).PubMed 
PubMed Central 

Google Scholar 
Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191. https://doi.org/10.1038/sdata.2017.191 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Wells, N., Goddard, S. & Hayes, M. J. A self-calibrating Palmer Drought Severity Index. J. Clim. 17, 2335–2351 (2004).
Google Scholar 
Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60. https://doi.org/10.1186/gb-2011-12-6-r60 (2011).Article 
PubMed 
PubMed Central 

Google Scholar 
Fisher, M. C. et al. RACE: Risk assessment of chytridiomycosis to European Amphibian Biodiversity. Froglog 101, 45–47 (2012).
Google Scholar  More

Berghauser Pont, M. Y., Perg, P. G., Haupt, P. A. & Heyman, A. A systematic review of the scientifically demonstrated effects of densification. IOP Conf. Ser. Earth Environ. Sci. 588, 052031 (2020).
Google Scholar 
Cimburova, Z. & Berghauser Pont, M. Location matters: A systematic review of spatial contextual factors mediating ecosystem services of urban trees. Ecosyst. Serv. 50, 101296 (2021).
Google Scholar 
De Valck, J. et al. Valuing urban ecosystem services in sustainable brownfield redevelopment. Ecosyst. Serv. 35, 139–149 (2019).
Google Scholar 
Zuzolo, D. et al. Divide et disperda: Thirty years of fragmentation and impacts on the eco-mosaic in the case study of the metropolitan city of Naples. Land 10, 485 (2021).
Google Scholar 
Nelson, E. The Economics of Ecosystems and Biodiversity: Ecological and Economic Foundations , edited by Pushpam Kumar, London, Earthscan Publications, United Nations Environment Programme, 2010, xxxix + 410 pp., US$76.95 (hardback), ISBN 978-1-84971-212-5. J. Nat. Resour. Policy Res. 5, 68–70 (2013).
Google Scholar 
Duraiappah, A. K. et al. Millennium Ecosystem Assessment, 2005. Ecosystems and human well-being: Synthesis. World Resources Institute vol. 5 http://www.who.int/entity/globalchange/ecosystems/ecosys.pdf (2005).Cariñanos, P., Casares-Porcel, M. & Quesada-Rubio, J. M. Estimating the allergenic potential of urban green spaces: A case-study in Granada, Spain. Landsc. Urban Plan. 123, 134–144 (2014).
Google Scholar 
Haase, D. et al. A quantitative review of urban ecosystem service assessments: Concepts, models, and implementation. Ambio 43, 413–433 (2014).PubMed 
PubMed Central 

Google Scholar 
Mexia, T. et al. Ecosystem services: Urban parks under a magnifying glass. Environ. Res. 160, 469–478 (2018).CAS 
PubMed 

Google Scholar 
Brzoska, P., Grunewald, K. & Bastian, O. A multi-criteria analytical method to assess ecosystem services at urban site level, exemplified by two German city districts. Ecosyst. Serv. 49, 101268 (2021).
Google Scholar 
Zulian, G. et al. Practical application of spatial ecosystem service models to aid decision support. Ecosyst. Serv. 29, 465–480 (2018).PubMed 
PubMed Central 

Google Scholar 
Balmford, A. et al. Ecology: Economic reasons for conserving wild nature. Science (80-). 297, 950–953 (2002).ADS 
CAS 

Google Scholar 
Koulov, B., Ivanova, E., Borisova, B., Assenov, A. & Ravnachka, A. GIS-based valuation of ecosystem services in mountain regions: A case study of the Karlovo municipality in Bulgaria. One Ecosyst. 2, e14062 (2017).
Google Scholar 
Robertson, G. P. & Swinton, S. M. Reconciling agricultural productivity and environmental integrity: A grand challenge for agriculture. Front. Ecol. Environ. 3, 38–46 (2005).
Google Scholar 
Sandhu, H. S., Wratten, S. D., Cullen, R. & Case, B. The future of farming: The value of ecosystem services in conventional and organic arable land. An experimental approach. Ecol. Econ. 64, 835–848 (2008).
Google Scholar 
Berglihn, E. C. & Gómez-Baggethun, E. Ecosystem services from urban forests: The case of Oslomarka, Norway. Ecosyst. Serv. 51, 101358 (2021).
Google Scholar 
Nowak, D. J. Understanding i-Tree. (2020). https://doi.org/10.2737/NRS-GTR-200.Selvakumaran, S., Plank, S., Geiß, C., Rossi, C. & Middleton, C. Remote monitoring to predict bridge scour failure using Interferometric Synthetic Aperture Radar (InSAR) stacking techniques. Int. J. Appl. Earth Obs. Geoinf. 73, 463–470 (2018).ADS 

Google Scholar 
Gómez-Baggethun, E. & Barton, D. N. Classifying and valuing ecosystem services for urban planning. Ecol. Econ. 86, 235–245 (2013).
Google Scholar 
Gren, Å. & Andersson, E. Being efficient and green by rethinking the urban-rural divide—Combining urban expansion and food production by integrating an ecosystem service perspective into urban planning. Sustain. Cities Soc. 40, 75–82 (2018).
Google Scholar 
Grêt-Regamey, A., Celio, E., Klein, T. M. & Wissen Hayek, U. Understanding ecosystem services trade-offs with interactive procedural modeling for sustainable urban planning. Landsc. Urban Plan. 109, 107–116 (2013).
Google Scholar 
Bennett, E. M., Peterson, G. D. & Gordon, L. J. Understanding relationships among multiple ecosystem services. Ecol. Lett. 12, 1394–1404 (2009).PubMed 

Google Scholar 
Bradford, J. B. & D’Amato, A. W. Recognizing trade-offs in multi-objective land management. Front. Ecol. Environ. 10, 210–216 (2012).
Google Scholar 
Cueva, J. et al. Synergies and trade-offs in ecosystem services from urban and peri-urban forests and their implication to sustainable city design and planning. Sustain. Cities Soc. 82, 103903 (2022).
Google Scholar 
Allocca, V., Coda, S., Calcaterra, D. & De Vita, P. Groundwater rebound and flooding in the Naples’ periurban area (Italy). J. Flood Risk Manag. 15, e12775 (2022).
Google Scholar 
Padulano, R. et al. Using the present to estimate the future: A simplified approach for the quantification of climate change effects on urban flooding by scenario analysis. Hydrol. Process. 35, e14436 (2021).
Google Scholar 
D’Amato, G. et al. Allergenic pollen and pollen allergy in Europe. Allergy 62, 976–990 (2007).PubMed 

Google Scholar 
Prigioniero, A., Zuzolo, D., Sciarrillo, R. & Guarino, C. Assessing pollinosis risk in the Vesuvius National Park: A novel approach for Index of Urban Green Zones Allergenicity. Environ. Res. 197, 111063 (2021).CAS 
PubMed 

Google Scholar 
AgCult 2020 Classifica visitatori 2019: Capodimonte rientra nella classifica dei primi 30 musei d’Italia.La Valva, V., Guarino, C., De Natale, A., Cuozzo, V., Menale, B. La flora del Parco di Capodimonte di Napoli. in 33–34: 143–177. (Delpinoa, 1992).Stevens, P. F. Angiosperm Phylogeny Website. 2001. http://www.mobot.org/MOBOT/research/APweb/. (2017).James Barth, B., Ian FitzGibbon, S. & Stuart Wilson, R. New urban developments that retain more remnant trees have greater bird diversity. Landsc. Urban Plan. 136, 122–129 (2015).
Google Scholar 
Heckmann, K. E., Manley, P. N. & Schlesinger, M. D. Ecological integrity of remnant montane forests along an urban gradient in the Sierra Nevada. For. Ecol. Manage. 255, 2453–2466 (2008).
Google Scholar 
Prigioniero, A. et al. Role of historic gardens in biodiversity-conservation strategy: the example of the Giardino Inglese of Reggia di Caserta (UNESCO) (Italy). Plant Biosyst. 155, 983–993 (2021).
Google Scholar 
Song, Q., Wang, B., Wang, J. & Niu, X. Endangered and endemic species increase forest conservation values of species diversity based on the Shannon-Wiener index. IForest 9, 469–474 (2016).
Google Scholar 
Hess, M. C. M., Mesléard, F. & Buisson, E. Priority effects: Emerging principles for invasive plant species management. Ecol. Eng. 127, 48–57 (2019).
Google Scholar 
Carli, E. et al. Using vegetation dynamics to face the challenge of the conservation status assessment in semi-natural habitats. Rend. Lincei. Sci. Fis. e Nat. 29, 363–374 (2018).
Google Scholar 
Canedoli, C. et al. Evaluation of ecosystem services in a protected mountain area: Soil organic carbon stock and biodiversity in alpine forests and grasslands. Ecosyst. Serv. 44, 101135 (2020).
Google Scholar 
FAO. Global Forest Resources Assessment 2010. Main report. (2010).Lindén, L., Riikonen, A., Setälä, H. & Yli-Pelkonen, V. Quantifying carbon stocks in urban parks under cold climate conditions. Urban For. Urban Green. 49, 126633 (2020).
Google Scholar 
Nowak, D. J., Hirabayashi, S., Bodine, A. & Greenfield, E. Tree and forest effects on air quality and human health in the United States. Environ. Pollut. 193, 119–129 (2014).CAS 
PubMed 

Google Scholar 
Nowak, D. J., Crane, D. E. & Stevens, J. C. Air pollution removal by urban trees and shrubs in the United States. Urban For. Urban Green. 4, 115–123 (2006).
Google Scholar 
Nowak, D. J. & Crane, D. E. Carbon storage and sequestration by urban trees in the USA. Environ. Pollut. 116, 381–389 (2002).CAS 
PubMed 

Google Scholar 
Kocić, K., Spasić, T., Urošević, M. A. & Tomašević, M. Trees as natural barriers against heavy metal pollution and their role in the protection of cultural heritage. J. Cult. Herit. 15, 227–233 (2014).
Google Scholar 
Yang, J., McBride, J., Zhou, J. & Sun, Z. The urban forest in Beijing and its role in air pollution reduction. Urban For. Urban Green. 3, 65–78 (2005).
Google Scholar 
Zupancic, T., Westmacott, C., Bulthuis, M. The impact of green space on heat and air pollution in urban communities: A meta-narrative systematic review (2015).Cariñanos, P., Adinolfi, C., Díaz de la Guardia, C., De Linares, C. & Casares-Porcel, M. Characterization of Allergen Emission Sources in Urban Areas. J. Environ. Qual. 45, 244–252 (2016).PubMed 

Google Scholar 
D’Auria, A., De Toro, P., Fierro, N. & Montone, E. Integration between GIS and multi-criteria analysis for ecosystem services assessment: A methodological proposal for the National Park of Cilento, Vallo di Diano and Alburni (Italy). Sustain 10, 3329 (2018).
Google Scholar 
Prigioniero, A., Zuzolo, D., Niinemets, Ü. & Guarino, C. Nature-based solutions as tools for air phytoremediation: A review of the current knowledge and gaps. Environ. Pollut. 277, 116817 (2021).CAS 
PubMed 

Google Scholar 
Szkop, Z. Evaluating the sensitivity of the i-Tree Eco pollution model to different pollution data inputs: A case study from Warsaw, Poland. Urban For. Urban Green. 55, 126859 (2020).
Google Scholar 
Tao, J. et al. Elevation-dependent effects of growing season length on carbon sequestration in Xizang Plateau grassland. Ecol. Indic. 110, 105880 (2020).CAS 

Google Scholar 
Chen, Y. et al. Grassland carbon sequestration ability in China: A new perspective from Terrestrial Aridity Zones. Rangel. Ecol. Manag. 69, 84–94 (2016).
Google Scholar 
Gopalakrishnan, V., Hirabayashi, S., Ziv, G. & Bakshi, B. R. Air quality and human health impacts of grasslands and shrublands in the United States. Atmos. Environ. 182, 193–199 (2018).ADS 
CAS 

Google Scholar 
Pace, R. et al. Comparing i-Tree eco estimates of particulate matter deposition with leaf and canopy measurements in an urban mediterranean Holm Oak Forest. Environ. Sci. Technol. 55, 6613–6622 (2021).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Losos, J. B., Walton, B. M. & Bennett, A. F. Trade-offs between sprinting and clinging ability in Kenyan Chameleons. Funct. Ecol. 7, 281 (1993).
Google Scholar 
Pretzsch, H., Moser-Reischl, A., Rahman, M. A., Pauleit, S. & Rötzer, T. Towards sustainable management of the stock and ecosystem services of urban trees. From theory to model and application. Trees – Struct. Funct. (2021). https://doi.org/10.1007/s00468-021-02100-3.Grunewald, K. et al. Lessons learned from implementing the ecosystem services concept in urban planning. Ecosyst. Serv. 49, 101273 (2021).
Google Scholar 
Baldacchini, C., Sgrigna, G., Clarke, W., Tallis, M. & Calfapietra, C. An ultra-spatially resolved method to quali-quantitative monitor particulate matter in urban environment. Environ. Sci. Pollut. Res. 26, 18719–18729 (2019).CAS 

Google Scholar 
De Luca, P., Guarino, C., Gullo, G., La Valva V., 1992. Il Parco di Capodimonte di Napoli: storia ed attualità. in 33–34: 143–177. (Delpinoa, 1992).Pignatti, S. Flora d’Italia vol.2. (2017).Braun-Blanquet, J. Plant Sociology (McGraw-Hill Book Company, 1932).
Google Scholar 
Catorci, A. et al. Reproductive traits variation in the herb layer of a submediterranean deciduous forest landscape. Plant Ecol. 214, 737–749 (2013).
Google Scholar 
Šumrada, T. et al. Are result-based schemes a superior approach to the conservation of High Nature Value grasslands? Evidence from Slovenia. Land Use Policy 111, 105749 (2021).
Google Scholar 
POWO. Plants of the World Online. Facilitated by the Royal Botanic Gardens, Kew. Board of Trustees of the Royal Botanic Gardens, Kew http://www.plantsoftheworldonline.org/ (2022).Bímová, K., Mandák, B. & Kašparová, I. How does Reynoutria invasion fit the various theories of invasibility?. J. Veg. Sci. 15, 495–504 (2004).
Google Scholar 
Wild, J., Neuhäuslová, Z. & Sofron, J. Changes of plant species composition in the Šumava spruce forests, SW Bohemia, since the 1970s. For. Ecol. Manag. 187, 117–132 (2004).
Google Scholar 
Damato, G. & Lobefalo, G. Allergenic pollens in the southern Mediterranean area. J. Allergy Clin. Immunol. 83, 116–122 (1989).CAS 

Google Scholar 
Cariñanos, P. et al. Assessing allergenicity in urban parks: A nature-based solution to reduce the impact on public health. Environ. Res. 155, 219–227 (2017).PubMed 

Google Scholar 
Cariñanos, P. et al. Estimation of the allergenic potential of urban trees and urban parks: Towards the healthy design of urban green spaces of the future. Int. J. Environ. Res. Public Health 16, 1357 (2019).PubMed Central 

Google Scholar  More

Systematic palaeontologyOrder Hemiptera Linnaeus, 1758Suborder Sternorrhyncha Amyot et Audinet-Serville, 1843Superfamily Protopsyllidioidea Carpenter, 1931Family Postopsyllidiidae Hakim, Azar et Huang, 2019Genus Megalophthallidion Drohojowska et Szwedo, gen. nov.LSID urn:lsid:zoobank.org:act:A6F71390-9B8E-4A19-8F30-C2A024B6EFB1Type speciesMegalophthallidion burmapateron Drohojowska et Szwedo, sp. nov.; by present designation and monotypy.EtymologyGeneric name is derived from Classic Greek megas (μέγας)—large, ophthalmos (ὀφθαλμός)—an eye and Greek form of generic name Psyllidium. Gender: masculine.Type localityNorthern Myanmar: state of Kachin, Noije bum 2001 Summit Site amber mine in the Hukawng Valley, SW of Maingkhwan.Type stratumLowermost Cenomanian, Upper Cretaceous (‘mid-Cretaceous’).DiagnosisHead capsule with 12 stiff setae on tubercles (18 setae in Postopsyllidium); fore wing without pterostigma (tiny pterostigma, widening of ScP + RA present in Postopsyllidium); vein CuP not thickened distally (distinctly thickened distally in Postopsyllidium); profemur with a row of ventral (ventrolateral) setae (two rows in Postopsyllidium).Megalophthallidion burmapateron Drohojowska et Szwedo, sp. nov.LSID urn:lsid:zoobank.org:act:F3F971F4-AE04-4F41-98B0-9A0A04470625.(Figs. 1A–F, 2A–I).Figure 1Megalophthallidion burmapteron gen. et sp. nov., holotype (MAIG 6687), imago. (A) Photo of body, ventral side; (B) photo of right antennae and (C) drawing of antenna; (D) drawing of body, dorsal side; (E) drawing of thorax structure with sclerites marked: red—pronotum; orange—mesopraescutum; yellow—mesoscutum; light green—mesoscutellum, dark green—mesopostnotum; light blue—metascutum; dark blue—metascutellum; violet—metapostnotum; (F) photo of thorax dorsal side. Scale bars: 0.5 mm (A), 0.2 mm (B–D), 0.1 mm (F).Full size imageFigure 2Megalophthallidion burmapteron gen. et sp. nov., holotype (MAIG 6687), imago. (A) Photo of right fore wing; (B) photo of right wings; (C) photo of antenna and proleg; (D) photo of proleg and mesoleg, and (E) photo of femur of proleg, and (F) photo of right metatarsus and left mesotarsus in the background, and (G) photo of right mesotarsus of mesoleg, and (H) Photo of tarsi; (I) photo of male genital block. Scale bars: 0.5 mm (A–D), 0.2 mm (B,E,F,H), 0.1 mm (G,I).Full size imageMaterialHolotype, number MAIG 6687 (BUB 96), deposited in Museum of Amber Inclusions (MAIG), University of Gdańsk, Poland. Imago, a complete and well-preserved male. Piece of amber 8 × 6 × 3 mm, cut from larger lump, polished flat on both sides.Type localityNorthern Myanmar: state of Kachin, Noije bum 2001 Summit Site amber mine in the Hukawng Valley, SW of Maingkhwan.Type stratumLowermost Cenomanian, Upper Cretaceous (‘mid-Cretaceous’).DiagnosisAs for the genus with the following additions: three ocelli distinct, antennomere IX the longest, about as long as pedicel, antennomeres III–VII and XI of similar length, antennomere XII the shortest, subconically tapered in apical portion. Paramere lobate, ventral margin with acute, small process, apical and dorsal margins rounded. Aedeagus geniculately bent at base, directed dorsally, tapered apicad.DescriptionMale (Figs. 1A–F, 2A–I). Head with compound eyes distinctly wider than pronotum (Fig. 1D–F). Compound eyes subglobular, protruding laterally. Vertex short in midline, about 2.5 times as wide as posterior margin and as long in middle; trapezoidal, anterior margin slightly arched, lateral margins diverging posteriad, posterior margin shallowly arched, disc of vertex with distinct setae on large tubercles: four setae at posterior margin, two at anterior angles of compound eyes, two medial, over the median ocellus. Three ocelli present, median ocellus distinct, visible from above, lateral ocelli near anterior angles of compound eyes. Frons about as wide as long in midline, two rows of setae on tubercles, upper row at level of median ocellus, lower one, below half of compound eye height. Clypeus, elongate, triangular, in lower portion roof-like; two setae on tubercles near upper margin. Genae very narrow. Rostrum reaching slightly beyond mesocoxae, apical segment slightly shorter than subapical one, darker. Antennae bases placed at lower margin of compound eyes; antennal fovea elevated; scapus shorter than pedicel, cylindrical; pedicel cylindrical; antennomeres IIIrd–VIIth and XIth of similar length, VIIIth slightly longer than VIIth, as long as Xth antennomere, IXth the longest, XIIth the shortest, tapered apically; rhinaria absent.Thorax (Fig. 1D–F): pronotum quadrangular, about as long as mesothorax; pronotum with anterior and posterior margins parallel, merely arcuate, disc with transverse groove in the median portion, lateral margins slightly arcuate, two distinct setae on tubercles in anterolateral angle, two setae on tubercles anterior margin at distance1/3 to median line, three distinct setae on tubercles in posterolateral angles. Mesopraescutum subtriangular, with apex widely rounded, about 0.4 times as wide as pronotum, about 0.4 times as long as wide, delicately separated from mesoscutum. Mesoscutum as wide as pronotum at widest point, distinctly narrowed medially, anterior angles rounded, anterolateral margin sigmoid, lateral angle acute, posterior angles wide, posterior margin V-shape incised, posterolateral areas of mesoscutum disc declivent posteriorly; disc with two setae on tubercles, at 1/3 of mesoscutum width. Mesoscutellum about as long as wide, diamond-shape, anterior and lateral angles acute, posterior angle rounded. Mesopostnotum in form of transverse band, slightly widened in median portion. Metascutum narrower than mesoscutum, anterior angles widely rounded, lateral angles acute, anterolateral margin concave, posterior margin arcuate, with deep median arcuate incision. The suture between metascutum and metascutellum weakly visible, metascutellum subtriangular, longer than wide at base.Parapteron with three distinct setae.Fore wing (Fig. 2A,B) membranous, narrow, elongate, about 3.5 times as long as wide, widest at 2/3 of length. Anterior margin merely arcuate, slightly bent at very base, anteroapical angle widely arcuate, apex rounded, posteroapical angle widely arcuate, tornus arcuate, claval margin straight, with incision between terminals of Pcu (claval apex) and A1. Stem ScP + R + MP + CuA slightly arcuate, very short stalk ScP + R + MP + CuA leaving basal cell, stem ScP + R oblique, straight, forked in basal half of fore wing length, branch ScP + RA, oblique, reaching anterior margin slightly distally of half of fore wing length, slightly distally of ending of CuA2 branch; branch RP slightly arcuate, a little more curved in basal section, reaching margin at anteroapical angle; stalk MP + CuA slightly shorter than basal cell; stem MP almost straight, forked in apical half of fore wing, at about 2/3 of fore wing length, with three terminals reaching margin between apex and posteroapical angle; stem CuA shorter than branches CuA1 and CuA2, about half as long as branch CuA1; claval vein CuP weak at base, not thickened distally; claval vein Pcu straight, claval vein A1 straight. Basal cell present, subtriangular, about twice as long as wide, basal veinlet cua-cup oblique, no other veinlets present; cell r (radial) very long, longer than half of fore wing length; cell m (medial) the shortest, shorter than cell cu (areola postica). Margins of fore wing with fringe of long setae, starting on costal margin near base of fore wing, ending at level of middle of cell cu; longitudinal veins with distinct, scarcely but evenly dispersed, movable setae; terminal section of CuP with two setae; costal margin with row of short, densely distributed setae, apical margin, tornus and claval margin with rows of scaly setae.Hind wing (Fig. 2B) membranous, shorter than fore wing, 3.23 times as long as wide. Costal margin bent at base, then almost straight up to the level of ScP + RA end and wing coupling lobe, then straight to anteroapical angle, anteroapical angle widely arcuate, apex arcuate, posteroapical angle arcuate, tornus straight, claval margin merely arcuate, posteroclaval angle angulate; stem ScP + R + MP bent at base, then straight, stem ScP + R short, branch ScP + RA short, about as long as stem ScP + R, branch RP arcuate basally than straight, reaching apex; stem MP arcuate, forked slightly distad CuA1 terminus level, branch MP1+2 slightly arcuate, reaching margin at posteroapical angle, branch MP3+4 straight, reaching tornus; stem CuA slightly bent at base, then straight, forked slightly distad ScP + R forking, branch CuA1 arcuate, branch CuA2 short, straight, slightly oblique, reaching tornus; claval vein CuP weak, visible only at base, claval vein Pcu slightly arcuate; wing coupling apparatus (fold) with a few short setae.Legs slender, relatively long, profemora armed (Fig. 2C–H). Procoxa as long as profemur, narrow, flattened. Protrochanter scaphoid, elongate, with long apical and subapical setae. Profemur flattened laterally, about as long as protibia, ventrally armed with four large setae on elevated plinths; dorsal margin with row of short, decumbent setae. Protibia narrow, rounded in cross section, covered with short setae, a few longer setae in distal portion. Protarsus—single, long tarsomere, plantar surface with row of semi-erect setae; tarsal claws long, straight, directed ventrally, no arolium nor empodium.Mesocoxa elongate, narrow, slightly flattened. Mesotrochanter scaphoid. Mesofemur slender, flattened laterally, dorsal margin with short setae. Mesotibia subequal to mesofemur, slender, covered with setae, two apical setae slightly thicker and longer. Mesotarsus with three tarsomeres, basimesotarsomere the longest, shorter than cumulative length of mid- and apical mesotarsomere, plantar margins with setae, two apical setae slightly longer and thicker; midmesotarsomere the shortest, 1/3 of basimesotarsomere length, a few setae on plantar surface; apical tarsomere shorter than basimesotarsomere, twice as long as midmesotarsomere, plantar surface with a few, scarcely dispersed setae, tarsal claws long, narrow, directed ventrally, no arolium nor empodium.Metacoxa conical, narrow. Metatrochanter scaphoid, elongate. Metafemur slender, laterally flattened, longer than mesofemur, dorsal margin with row of short setae. Metatibia, long, slender, 1.6 times as long as metafemur, with suberect setae of different size, two larger and longer and two shorter setae subapical setae. Metatarsus slightly less than half of metatibia length, with three tarsomeres, basimetatarsomere the longest, more than twice as long as apical metatarsomere, 1.5 times as long as combined length of mid- and apical metatarsomere, plantar surface with scarce decumbent setae; mid metatarsomere the shortest, 1/4 of basimetatarsomere length, plantar surface with a few setae, two apical ones slightly thicker; apical metatarsomere about 0.4 of basimetatarsomere length, with scarcely dispersed setae on along plantar surface; tarsal claws, long, slender, other pretarsal structures absent.Abdomen (Fig. 1F) narrowly attached to thorax, tergite segment shorter, 2nd tergite distinctly longer, 3rd to 8th tergites of similar length; pygofer narrowing apicad, ventral margin strongly elongated posteriorly; anal tube short, directed posterodorsad, anal style shorter than anal tube. Paramere lobate, ventral margin with acute, small process, apical and dorsal margins rounded. Aedeagus (Fig. 2I) geniculately bent at base, directed dorsad, tapered apicad.Female. Unknown.Megalophthallidion sp. (5th instar nymph)(Figs. 3A–D, 4A–F)Figure 3Megalophthallidion sp. (MAIG 6688), nymph. (A) Photo of body, dorsal side and (B) drawing of body dorsal side; (C) photo of body dorsal side and (D) drawing of body ventral side. Scale bars: 0.5 mm (A–D).Full size imageFigure 4Megalophthallidion sp. (MAIG 6688), nymph. Photo of clypeus and (B) drawing of clypeus; (C) photo of proleg, and (D) photo of mesoleg, and (E) photo of metaleg; (F) photo of posterior part of abdomen ventral side. Scale bars: 0.1 mm (A–F).Full size imageMaterialNymph, 5th instar, MAIG 6688 (BUB 1799), deposited in Museum of Amber Inclusions (MAIG), University of Gdańsk, Poland. Piece of amber 13 × 6 × 2 mm, cut from larger lump, polished flat on one side, more convex on the other.Diagnostic charactersThe nymph of Megalophthallidion gen. nov. is similar in general body shape to the only known fossil protopsyllidioidean nymph described from Lower Cretaceous Lebanese amber—Talaya batraba Drohojowska et Szwedo, 2013. The nymph of Talaya batraba is 2nd or 3rd instar, therefore some features are difficult to compare with this last instar nymph of Megalophthallidion gen. nov. The morphological states observed in those two specimens are: head covered with strongly expanded disc and expanded disc of pronotum, however shapes and ratios of these structures differ; compound eyes on ventral side of head, shifted laterad (ommatidia on cones in T. batraba, while ventroposterior expansions are present in Megalophthallidion gen. nov.); compound eyes visible from above as short, stout cones in fissure between posterior margin of disc (hypertrophied vertex) and anterior margin of pronotum (compound eyes (?) are visible on dorsal side of Permian Aleuronympha bibulla Riek, 1974); in Megalophthallidion gen. nov. rostrum reached mesocoxa, while in Talaya batraba distinctly exceeds length of the body; abdomen with 9 segments; tergites of abdominal segments 5th–9th expanded posterolaterad in form of fan-like expansion; 9th abdominal segment short, placed ventral; anal tube short, cylindrical, epiproct (?) globular.DescriptionNymph, 5th instar (Figs. 3A–D, 4A–F). Body oval shaped, dorso-ventrally flattened, 1.5 times longer than wide with segmentation visible; on the ventral side slightly concave. Length of body c. 1.56 mm long, outline, in dorsal view, maximum width of body 0.94 mm; length of head and pronotum (cephaloprothorax) c. 0.46 mm in midline, width c. 0.83 mm; cumulative length of mesonotum + metanotum c. 0.25 mm; abdomen c. 0.8 mm long. Dorsal side (Fig. 3A,B) with distinct median line (ecdysial line), not reaching anterior or posterior margin of the body, the line distinctly roof-like in abdominal portion. Anterior margin of head (cephaloprothorax) disc arcuate, lateral angles rounded; anterior margin of pronotum arcuate, lateral margins arcuately diverging posteriad, posterior margin distinctly arcuate, anterior angles widely rounded, posterior angles acutely rounded, disc elevated, convex, lateral portions declivitous; the fissure between posterior margin of head disc and anterior margin of pronotum narrow, widened medially, with stalked compound eyes popping out.Head partly separated from prothorax, wide in ventral view. Bases of antennae protruding anterolaterally, wide, anterior margin arcuate, with a small lump extending anteriorly connecting margin with vertex expansion. Suture separating anteclypeus and postclypeus visible in ventral aspect (Fig. 4A,B). Postclypeus about three times as long as wide, oval, slightly swollen, without any setae; weak traces of salivary pump muscle attachments visible. Anteclypeus about as long as postclypeus, widened in upper section below clypeal suture, convex, carinately elevated in lower section, with sides distinctly declivitous, clypellus long, carinately elevated. Lora (mandibulary plates) distinct, separated from anteclypeus by shallow suture, with upper angles at half of postclypeus length, lower angles at half of anteclypeus length, about as wide as half of postclypeus width. Maxillary plates narrow. Genal portion of head enlarged, medial portion arcuately convex; lateral sections narrowing laterally, terminally encircling bases of compound eyes. Antennae short (Fig. 3C,D), placed in front of genal portion. Antennal flagellum indistinctly subdivided into four segments. Rostrum (Fig. 4A,B) three-segmented, 0.2 mm long, with apex reaching apex of mesocoxae; apical segment about 2.5 times as long as subapical one.No lateral sclerites on meso- and metathorax, only one plus one large medial sclerite on both meso- and metathorax. Mesothoracic and metathoracic wing pads distinct, wide, subtriangular, with posterior apices directed posteriorly; lateral portions of mesothoracic wing pads arcuate. Fore wing pad 0.6 mm long, with small, straight humeral lobe, forming a right angle, not protruding anteriorly. Mesothoracic tergites slightly larger than metathoracic segments (respectively c. 0.14 mm and c. 0.12 mm long in midline, 0.26 mm and 0.27 mm in lateral lines); mesothoracic tergum with distinct median elevation (low double crest with ecdysial line in between), slightly wider than long in midline, anterior margin arcuate, lateral margins straight, subparallel, posterior margin concave. Metathoracic wing pad apex slightly exceeding mesothoracic wing pad. Metathoracic tergum wider than long, slightly shorter than mesothoracic tergum, with distinct elevation in the middle.Legs relatively long (Figs. 3C,D, 4C–E). Coxae of legs placed near the median axis of the body. Prolegs: procoxal pit with margins elevated, procoxa conical (c. 0.1 mm long), protrochanter scaphoid, about as long as procoxa, profemur c. 0.13 mm long, slightly flattened laterally, merely thickened, protibia longer than profemur, c. 0.23 mm long; tarsus shorter than protibia, basiprotarsomere about as long as apical protarsomere, the latter with distinct tarsal claws, and wide arolium. Mesoleg similar to proleg, mesocoxa conical (c. 0.1 mm long), mesotrochanter scaphoid, mesofemur (c. 0.13 mm) slightly flattened laterally, mesotibia slightly longer than mesofemur (c. 0.18 mm), mesotarsus slightly shorter than mesotibia, three-segmented, basimesotarsomere the longest (c. 0.07 mm), about as long as combined length of mid- and apical mesotarsomeres (c. 0.04 mm respectively), arolium wide, tarsal claws distinct. Metaleg: metacoxa conical (c. 0.1 mm), metatrochanter scaphoid, about as long as metacoxa (c. 0.12 mm). Metafemur (c. 0.17 mm) slightly more thickened than pro- and mesofemur, metatibia slightly longer (0.19 mm) than pro- and mesotibiae. Metatarsus three-segmented: basimetatarsomere about as long (0.08 mm) as combined length of mid- and apical metatarsomeres (0.04 mm respectively), arolium lobate, wide, tarsal claws distinct, widely spread.Abdomen (Fig. 3A–D) 9-segmented, narrow at base, widening fan-shape posteriorly, 1st segment visible from above, segmentation visible, abdominal terga 5th–9th expanded posterolaterally. Tergites carinately elevated in the middle, separated by ecdysial line. 1st sternite visible in ventral view, sternites 2nd–4th fused medially, sternites 5th–9th separated; 9th abdominal segment short (Fig. 4F), placed ventrally, under tergal expansion; anal tube short, cylindrical, epiproct (?) globular. More

Analysis of contents of heavy metals in wasteland soilThe test results show (Table 5) that the contents of Hg, Cd, As, Pb, Cr, Zn, Ni and Cu in the surface soil within Shigetai Coal Mine vary from 0.043 to 0.255, 0.44 to 2.23, 2.66 to 18.40, 11.80 to 42.80, 40.50 to 118.60, 18.90 to 70.10, 4.31 to 28.10, 4.96 to 46.25 mg/kg, respectively; the average contents of Hg, Cd, As, Pb, Cr, Zn, Ni and Cu are 0.128, 1.03, 4.73, 23.08, 76.22, 46.94, 16.11 and 12.10 mg/kg, respectively. The average contents of Hg, Cd, Pb and Cr in soil within the research area are 2.03, 1.36, 1.11 and 1.23 times of the soil background values in Shaanxi Province, respectively. The average contents of As, Zn and Cu are lower than the soil background value in Shaanxi Province, but the maximum contents of these three elements are 1.65, 1.01 and 2.16 times of the soil background values in Shaanxi Province, respectively. It is reported that the average concentration of lead in agricultural soil affected by coal mines is relatively high (433 mg kg−1)38. Lead is usually related to minerals in coal and occurs mainly in the form of sulfide such as PbS and PbSe39. In addition, aluminosilicate and carbonate also contain lead40. Chromium is a non-volatile element, which is related to aluminosilicate minerals41. In the mining process, chromium may be accumulated in coal, gangue or other tailings, and then enter the soil or water body through rain leaching42.Table 5 Statistics of contents of heavy metals in wasteland soil (n = 79).Full size tableThe coefficient of variation (CV) of Hg and Cd contents in soil within the research area is 0.050 and 0.37, respectively, with moderate variation, indicating that the content of these two heavy metals is less affected by the external factors; the coefficient of variation (CV) of As, Pb, Cr, Zn, Ni and Cu contents is 2.81, 7.46, 18.00, 13.51, 5.44 and 5.64, respectively, with strong variation (CV  > 0.50)43, indicating that the content of these eight heavy metals may be affected by some local pollution sources. The skewness coefficient (SK) ranges from − 3 to 3, and the larger its absolute value, the greater its skewness. When SK  > 0, it is positive skewness; when SK  More

ChemicalsErbium (III) oxide nanoparticles (Er2O3-NPs) were purchased from Sigma-Aldrich Chemical Company (Saint Louis, USA) with pink appearance and product number (203,238). Powders of Er2O3-NPs with 99.9 trace metals basis were suspended in deionized distilled water to prepare the required concentrations and ultra-sonicated prior use.Cell lineHuman hepatocellular carcinoma (Hep-G2) cells were obtained from Nawah Scientific Inc., (Mokatam, Cairo Egypt). Cells were maintained in Dulbecco’s Modified Eagle Medium (DMEM) media supplemented with streptomycin (100 mg/mL), penicillin (100 units/mL) and heat-inactivated fetal bovine serum (10) in humidified, 5% (v/v) CO2 atmosphere at 37 °C.Characterization of Er2O3-NPsThe purchased powders of Er2O3-NPs were characterized using a charge coupled device diffractometer (XPERT-PRO, PANalytical, Netherlands) to determine its X-ray diffraction (XRD) pattern. Zeta potential and particles’ size distribution of Er2O3-NPs were also detected using Malvern Instrument Zeta sizer Nano Series (Malvern Instruments, Westborough, MA) equipped with a He–Ne laser (λ = 633 nm, max 5mW). Moreover, transmission electron microscopy (TEM) imaging was done to detect the shape and average particles’ size of Er2O3-NPs suspension.Sulforhodamine B (SRB) cytotoxicity assaySulforhodamine B (SRB) assay was conducted to assess the influence of Er2O3-NPs on the proliferation of cancerous Hep-G2 cells12. Aliquots of 100 µl of Hep-G2 cells suspension containing 5 × 103 cells were separately cultured in 96-well plates and incubated for 24 h in complete media. Hep-G2 Cells were then treated with five different concentrations of Er2O3-NPs (0.01, 0.1, 1, 10 and 100 µg/ml) incubated for 24 h or (0.1, 1, 10, 100 and 1000 µg/ml) incubated for 72 h. After 24 or 72 h of Er2O3-NPs exposure, cultured cells were fixed by replacing media with 10% trichloroacetic acid (TCA) and incubated for one hour at 4 °C. Cells were then washed five times with distilled water, SRB solution (0.4% w/v) was added and incubated cells in a dark place at room temperature for 10 min. All plates were washed three times with 1% acetic acid and allowed to air-dry overnight. Then, protein-bound SRB stain was dissolved by adding TRIS (10 mM) and the absorbance was measured at 540 nm using a BMG LABTECH-FLUO star Omega microplate reader (Ortenberg, Germany).Cells treatmentCancerous Hep-G2 cells were cultured at the appropriate conditions and dived into control and treated cells. The control cells were treated with an equal volume of the vehicle (DMSO; final concentration, ≤ 0.1%), while the treated cells were treated with the IC50 of Er2O3-NPs. All cells were left for 72 h after nanoparticles treatment and were harvested by brief trypsinization and centrifugation. Each treatment was conducted in triplicate. Cells were washed twice with ice-cold PBS and used for different molecular assays.Estimation of genomic DNA integrityThe impact of Er2O3-NPs exposure on the integrity of genomic DNA in cancerous Hep-G2 cells was estimated using alkaline Comet assay13,14. Treated and control cells were mixed with low melting agarose and spread on clean slides pre-coated with normal melting agarose. After drying, slides were incubated in cold lysis buffer for 24 h in dark and then electrophoresed in alkaline electrophoresis buffer. Electrophoresed DNA was neutralized in Tris buffer and fixed in cold absolute ethanol. For analysis slides were stained with ethidium bromide, examined using epi-fluorescent microscope at magnification 200× and fifty comet nuclei were analyzed per sample using Comet Score software.Estimation of intracellular ROS generationThe effect of Er2O3-NPs exposure on intracellular ROS production in cancer Hep-G2 cells was studied using 2,7-dichlorofluorescein diacetate dye15. Cultured cells were washed with phosphate buffered saline (PBS) and then 2,7-dichlorofluorescein diacetate dye was added. Mixed cells and dye were left for 30 min in dark and spread on clean slides. The resultant fluorescent dichlorofluorescein complex from interaction of intracellular ROS with dichlorofluorescein diacetate dye was examined under epi-fluorescent at 20× magnification.Measuring the expression levels of apoptotic and anti-apoptotic genesQuantitative real time Polymerase chain reaction (RT-PCR) was conducted to measure the mRNA expression levels of apoptotic (p53 and Bax) and anti-apoptotic (Bcl2) genes in control and treated Hep-G2 cells. Whole cellular RNA was extracted according to the instructions listed by the GeneJET RNA Purification Kit (Thermo scientific, USA) (Thermo scientific, USA) and using Nanodrop device purity and concentration of the extracted RNAs were determined. These RNAs were then reverse transcribed into complementary DNA (cDNA) using the instructions of the Revert Aid First Strand cDNA Synthesis Kit (Thermo scientific, USA). For amplification, RT-PCR was performed using the previously designed primers shown in Table 116,17 by the 7500 Fast system (Applied Biosystem 7500, Clinilab, Egypt). A comparative Ct (DDCt) method was conducted to measure the expression levels of amplified genes and GAPDH gene was used as a housekeeping gene. Results were expressed as mean ± S.D.Table 1 Sequences of the used primers in qRT-PCR.Full size tableAnalysis of cell cycle distributionDistribution of cell cycle was analyzed using flow cytometry. Control and treated cancer Hep-G2 cells with IC50 of Er2O3-NPs for 72 h were harvested, washed with PBS and re-suspended in 1 mL of PBS containing RNAase A (50 µg/mL) and propidium iodide (10 µg/mL) (PI). Cells were incubated for 20 min in dark at 37 C and analyzed for DNA contents using FL2 (λex/em 535/617 nm) signal detector (ACEA Novocyte flow cytometer, ACEA Biosciences Inc., San Diego, CA, USA). For each sample, 12,000 events are acquired and cell cycle distribution is calculated using ACEA NovoExpress software (ACEA Biosciences Inc., San Diego, CA, USA).Estimation of apoptosis inductionApoptotic and necrotic cell populations were determined using Annexin V- Fluorescein isothiocyanate (FITC) apoptosis detection kit (Abcam Inc., Cambridge Science Park Cambridge, UK) coupled with two fluorescent channels flow cytometry. After treatment with Er2O3-NPs for 72 h and doxorubicin as a positive control, Hep-G2 cells were collected by trypsinization and washed twice with ice-cold PBS (pH 7.4). Harvested cells are incubated in dark with Annexin V-FITC/ propidium iodide (PI) solution for 30 min at room temperature, then injected via ACEA Novocyte flowcytometer (ACEA Biosciences Inc., San Diego, CA, USA) and analyzed for FITC and PI fluorescent signals using FL1 and FL2 signal detector, respectively (λex/em 488/530 nm for FITC and λex/em 535/617 nm for PI). For each sample, 12,000 events were acquired and positive FITC and/or PI cells are quantified by quadrant analysis and calculated using ACEA NovoExpress software (ACEA Biosciences Inc., San Diego, CA, USA).Statistical analysisResults of the current study are expressed as mean ± Standard Deviation (S.D) and were analyzed using the Statistical Package for the Social Sciences (SPSS) (version 20) at the significance level p  More

General description of the experimental areaThe experiment was performed for two years at the National Key Irrigation Experimental Station located on the Songnen Plain in Heping town, Qing’an County, Suihua, Heilongjiang, China, with a geographical location of 45° 63′ N and 125° 44′ E at an elevation of 450 m above sea level (Fig. 1). This region consists of plain topography and has a semiarid cold temperate continental monsoon climate, i.e., a typical cold region with a black soil distribution area. The average annual temperature is 2.5 °C, the average annual precipitation is 550 mm, the precipitation is concentrated from June to September of each year, and the average annual surface evaporation is 750 mm. The growth period of crops is 156–171 days, and there is a frost-free period of approximately 128 days year−122. The soil at the study site is albic paddy soil with a mean bulk density of 1.01 g/cm3 and a porosity of 61.8% prevails. The basic physicochemical properties of the soil were as follows: the mass ratio of organic matter was 41.8 g/kg, pH value was 6.45, total nitrogen mass ratio was 15.06 g/kg, total phosphorus mass ratio was 15.23 g/kg, total potassium mass ratio was 20.11 g/kg, mass ratio of alkaline hydrolysis nitrogen was 198.29 mg/kg, available phosphorus mass ratio was 36.22 mg/kg and available potassium mass ratio was 112.06 mg/kg.Figure 1Location of the study area. The map and inset map in this image were drawn by the authors using ArcGIS software. The software version used was ArcGIS software v.10.2, and its URL is http://www.esri.com/.Full size imageHumic acid fertilizerHumic acid fertilizer was produced by Yunnan Kunming Grey Environmental Protection Engineering Co., Ltd., China (Fig. 2). The organic matter was ≥ 61.4%, and the total nutrients (nitrogen, phosphorus and potassium) were ≥ 18.23%, of which N ≥ 3.63%, P2O5 ≥ 2.03%, and K2O ≥ 12.57%. The moisture content was ≤ 2.51%, the pH value was 5.7, the worm egg mortality rate was ≥ 95%, and the amount of faecal colibacillosis was ≤ 3%. The fertilizer contained numerous elements necessary for plants. The contents of harmful elements, including arsenic, mercury, lead, cadmium and chromium, were ≤ 2.8%, 0.01%, 7.6%, 0.1% and 4.7%, respectively; these were lower than the test standard.Figure 2Humic acid fertilizer in powder form.Full size imageExperimental design and observation methodsIrrigationIn this experiment, three irrigation practices, namely, control irrigation (C), wet irrigation (W) and flood irrigation (F), were designed (Table 1).Table 1 Different irrigation methods.Full size tableControl irrigation (C) of rice had no water layer in the rest of the growing stages, except for the shallow water layer at the regreen stage of rice, which was maintained at 0–30 mm, and the natural dryness in the yellow stage. The irrigation time and irrigation quota were determined by the root soil moisture content as the control index. The upper limit of irrigation was the saturated moisture content of the soil, the lower limit of soil moisture at each growth stage was the percentage of saturated moisture content, and the TPIME-PICO64/32 soil moisture analyser was used to determine the soil moisture content at 7:00 a.m. and 18:00 p.m., respectively. When the soil moisture content was close to or lower than the lower limit of irrigation, artificial irrigation occurred until the upper irrigation limit was reached. The soil moisture content was maintained between the upper irrigation limit and the lower irrigation limit of the corresponding fertility stage. Under the wet irrigation (W) and flood irrigation (F) conditions, it was necessary to read the depth of the water layer through bricks and a vertical ruler embedded in the field before and after 8:00 am every day to determine if irrigation was needed. If irrigation was needed, then the water metre was recorded before and after each irrigation. The difference between before and after was the amount of irrigation23.FertilizationIn our research, five fertilization methods were applied, as shown in Table 2. In this experiment, the rice cultivar “Suijing No. 18” was selected. Urea and humic acid fertilizer were applied according to the proportion of base fertilizer:tillering fertilizer:heading fertilizer (5:3:2). The amounts of phosphorus and potassium fertilizers were the same for all treatments, and P2O5 (45 kg ha−1) and K2O (80 kg ha−1) were used. Phosphorus was applied once as a basal application. Potassium fertilizer was applied twice: once as a basal fertilizer and at 8.5 leaf age (panicle primordium differentiation stage) at a 1:1 ratio22.Table 2 The fertilizer methods.Full size tableThis study was performed with a randomized complete block design with three replications. Three irrigation practices and five fertilizer methods were applied, for a total of 15 treatments as follows: CT1, CT2, CT3, CT4, CT5; WT1, WT2, WT3, WT4, WT5; FT1, FT2, FT3, FT4, and FT5 (C, W, and F represent control irrigation, wet irrigation, and flood irrigation; T represents fertilizer treatment).Measurements of the samplesA soil temperature sensor (HZTJ1-1) was buried in each experimental plot to monitor the temperature of each soil layer (5 cm, 10 cm, 15 cm, 20 cm and 25 cm depth). The transmission of photosynthetically active radiation was measured from 11:00 to 13:00 by using a SunScan Canopy Analysis System (Delta T Devices, Ltd., Cambridge, UK), and data during the crop-growing season were recorded every day24.Plant measurements were taken during the periods of tillering to ripening on days with no wind and good light. The fluorescence parameters were measured by a portable fluorescence measurement system (Li-6400XT, America). The detection light intensity was 1500 μmol m−2 s−1, and the saturated pulsed light intensity was 7200 μmolm−2 s−1. The functional leaves were dark adapted for 30 min, and then the maximum photosynthetic efficiency of PSII (Fv/Fm) was measured. Photochemical quenching (QP) and nonphotochemical quenching (NPQ) were measured with natural light. Simultaneously, the leaf chlorophyll relative content (SPAD) was monitored using SPAD 502 (Konica Minolta, Inc., Tokyo, Japan). For plant agronomic characteristics, the distance from the stem base to the stem tip was measured with a straight ruler to quantify plant height24.Statistical analysisExperimental data obtained for different parameters were analysed statistically using the analysis of variance technique as applicable to randomized complete block design. Duncan’s multiple range test was employed to assess differences between the treatment means at a 5% probability level. All statistical analyses were performed using SPSS 22.0 for Windows24.
Ethics approvalExperimental research and field studies on plants, including the collection of plant material, comply with relevant institutional, national, and international guidelines and legislation. We had appropriate permissions/licences to perform the experiment in the study area. More

The overall magnitude of changes in SOC, TN, and C:N in response to chemical nitrogen fertilizers reductionThe results showed that chemical nitrogen fertilizers reduction significantly decreased SOC and TN by 2.76% and 4.19% respectively, while increased C:N by 6.11% across all database (Fig. 1). SOC mainly derives from crop residues and secretions which closely related to crops growths, and crops growths were affected by fertilization, especially nitrogen fertilization20,21. The reduction of chemical nitrogen fertilizer led to poor crop growth, which reduced the amount of crop residues return, and then decreased SOC. Similarly, TN from crops was reduced due to poor crop growth. In addition, the reduction of chemical nitrogen fertilizers directly reduced the input of soil nitrogen. The increase of C:N was the result of the decrease of TN being greater than that of SOC. The responses of C:N to chemical nitrogen fertilizers reduction enhanced the comprehension of the couple relationship between SOC and TN, which was beneficial to the evolution of the C-N coupling models. Moreover, the accuracy of C-N coupling models depends on the precise quantification of the responses of SOC and TN to nitrogen fertilization. And our results accurately quantified the difference responses of SOC and TN to different nitrogen fertilization regimes, thus optimizing the C-N coupling models.Figure 1The weighted response ratio (RR++) for the responses to chemical nitrogen fertilizers of soil organic carbon (SOC, a), total nitrogen (TN, b), and their ratios (C:N, c). Bars denote the overall mean response ratio RR++ and 95% confidence intervals (CI). The star (*) indicates significance when the 95% CI that do not go across the zero line. The vertical lines are drawn at lnRR = 0. The value represents independent sample size.Full size imageResponses of SOC, TN and C:N to chemical nitrogen fertilizers reduction magnitudeWhen grouped by chemical nitrogen fertilizers reduction magnitude, SOC showed a significant increase by 6.9% in medium magnitude, while SOC was significantly decreased by 3.10% and 7.26% in high and total magnitude respectively (Fig. 1a). Liu and Greaver22 also stated the reduction of medium nitrogen fertilizer increased the average microbial biomass from 15 to 20%, thereby increasing the SOC content. Previous studies had reported that there were strong positive correlations between soil organic matter and soil microbial biomass in both the agricultural ecosystem and natural ecosystem23,24. Numerous researchers have demonstrated the significance of nitrogen availability in soil for the plant biomass across most ecosystems25,26. Moreover, nitrogen deficient would inhibit the activity of extracellular enzymes and root activities27. Generally, soil degradation caused by continuous rising chemical nitrogen fertilizers application may inhibit the growth of crops and ultimately reduce the SOC28.TN showed no significant difference in low and medium chemical nitrogen fertilizers reduction magnitude (p  > 0.05), while TN in high magnitude and total chemical nitrogen fertilizers reduction magnitude exhibited a decrease with 3.10% and 9.37% respectively (Fig. 1b). Numerous studies described that the amount of nitrogen fertilizers used in China was higher than the demand of N for crop, which caused serious N leaching and runoff29,30. Chemical nitrogen fertilizers in low and medium magnitude would not decrease the TN of soil by reducing N leaching and runoff. However, the residual nitrogen in soil cannot meet the requirement for the sustainable growth of plant with litter or without exogenous nitrogen supplement, which resulted in the decrease of TN in high and total chemical nitrogen fertilizers magnitude. Consequently, optimal nitrogen fertilizers application rates will take into account crops yield and environment friendliness.Additionally, C:N had a significant increase with ranging from 1.82% to 8.98% under the four chemical nitrogen fertilizers reduction magnitude (Fig. 1c), suggesting the relative increase of SOC compared to TN. Previous studies have revealed that C:N had significantly influence on the soil bacterial community structures31. And there were also considerable studies indicated that chemical nitrogen fertilizers have impact on the soil bacterial communities32,33. We may speculate that the change of C:N would bring about the variations of soil bacteria communities under the chemical nitrogen fertilizers regimes.Responses of SOC, TN, and C:N to chemical nitrogen fertilizers reduction durationNegative response of SOC to short-term chemical nitrogen fertilizers reduction was observed in our study, which was consistent with the study of Gong, et al.34 that chemical nitrogen fertilizers reduction decreased SOC by reducing crop-derived carbon by one year. However, SOC was significantly increased by 1.06% and 4.65% at mid-term and long-term chemical nitrogen fertilizers reduction respectively, which was similar with the findings of Ning, et al.11 that SOC was significantly increased under more than 5 years of chemical nitrogen fertilizers reduction treatment. TN was significantly decreased by 1.96% at short-term chemical nitrogen fertilizers reduction duration, while the results converted at mid-term chemical nitrogen fertilizers reduction duration. The effect of long-term chemical nitrogen fertilizers reduction on TN was not significant (p  > 0.05). The divergent response of TN to different chemical nitrogen fertilizers duration was mainly caused by the various treatments. In terms of C:N, a greater positive response was observed at short-term chemical nitrogen fertilizers duration (9.06%) than mid-term and long-term duration (1.99%). Moreover, with the prolongation of the chemical reduction time of nitrogen, the response ratio tends to zero, suggesting that the effect of chemical fertilizers gradually decrease. This may be ascribed to the buffer capacity of soil to resist the changes from external environment, including nutrients, pollutants, and redox substances35.Responses of SOC, TN, and C:N to different chemical nitrogen fertilizers reduction patternsUnder the pattern of chemical nitrogen fertilizers reduction without organic fertilizers supplement, SOC and TN significantly decreased by 3.83% and 11.46% respectively, however, chemical nitrogen fertilizers reduction with organic fertilizers supplement significantly increased SOC and TN by 4.92% and 8.33% respectively. Moreover, C:N significantly increased under the two chemical nitrogen fertilizers patterns (p  0.05), but there was a negative effect on SOC in high and total magnitude (p  0.05). The no significant decrease at mid-term duration might result from the limited information reported in original studies of this meta-analysis36. TN showed no significant response to chemical nitrogen fertilizers without organic fertilizers supplement in the low and medium magnitude (p  > 0.05). However, TN was significantly decreased by 8.62% and 16.7% respectively in the high and total magnitude. When regarding to chemical nitrogen fertilizers reduction duration, TN was significantly reduced at all of the categories, ranging from 3.13% to 13.4% (Fig. 2c). In the pattern of chemical nitrogen fertilizers reduction with organic fertilizers supplement, chemical nitrogen fertilizers reduction at medium, high, and total magnitudes significantly increased SOC by 13.85%, 13.03%, and 5.46%respectively, however, the response of SOC in the low chemical nitrogen fertilizers magnitude was not significant. Chemical nitrogen fertilizers reduction duration significantly increased SOC by 7.01%, 1.71%, and 22.02% in the short-term, mid-term, and long-term respectively. Comparatively, TN showed a significantly increase in most chemical nitrogen fertilizers categories expect for the long-term chemical nitrogen fertilizers duration, with an increasing from 4.90% to 14.69% (Fig. 2d).Figure 2The weighted response ratio (RR++) for the responses to chemical nitrogen fertilizers of soil organic carbon (SOC, a), total nitrogen (TN, b), and their ratios (C:N, c) under the two patterns (with organic fertilizers ; without organic fertilizers). Bars denote the overall mean response ratio RR++ and 95% confidence intervals (CI). The star (*) indicates significance when the 95% CI that do not go across the zero line. The vertical lines are drawn at lnRR = 0. The values represent independent sample size.Full size imageOrganic fertilizers were mainly derived from animal manure or crops straws, which contained large amount of organic matter and nitrogen elements37,38. The application of organic fertilizers increased the input of SOC and TN directly. Moreover, organic fertilizer could promote the growth of crops by releasing phenols, vitamins, enzymes, auxins and other substances during the decomposition process, thus the SOC derived from crops would be increased37,39. In addition, organic fertilizers provide various nutrients for microbial reproduction, which increase the microbial population and organic carbon and total nitrogen content37. More importantly, the application of organic fertilizers could improve organic carbon sequestration and maintain its stability in aggregates, thereby reducing losses of SOC and TN40.C:N showed an increase under all of the chemical nitrogen fertilizers reduction with organic fertilizer supplement. The positive response of C:N to organic fertilizer supplement may be related to the higher C:N of organic fertilizer than soil. The average values of C:N of the commonly used organic fertilizers including animal manure, crop straws and biochar were 14, 60 and 30 respectively, while the C:N of soil was lower than 10 in average according to extensive literature researches41. Therefore, organic fertilizers would be a favorable alternative of chemical fertilizers for the sustainable development of agriculture.The correlation between the response of SOC, TN, and C:N and environmental variablesThe analysis of linear regression was conducted to analyze the environmental variables including mean annual temperature (MAT), mean annual precipitation (MAP), accumulated temperature above 10 °C (MATA), which may exert influence on SOC, TN and C:N. No significant correlation among the lnRR of SOC, TN, C:N and environmental variables were observed among the whole database (p  > 0.05; Fig. S1). Rule out the interference of organic fertilizers supplement, we analyzed the relationship between lnRR of SOC, TN, C:N and environmental variables as the Figures showed in Figs. 3 and 4 respectively. Under chemical nitrogen fertilizers without organic fertilizers supplement, there was a significant negative correlation between lnRR of SOC and MAT (p  More

Irigoien, X. et al. Nat. Commun. 5, 3271 (2014).Article 

Google Scholar 
Ariza, A. et al. Nat. Clim. Change https://doi.org/10.1038/s41558-022-01479-2 (2022).Article 

Google Scholar 
Klevjer, T. A. et al. Sci. Rep. 6, 19873 (2016).CAS 
Article 

Google Scholar 
Braun, C. D. et al. Annu. Rev. Mar. Sci. 14, 129–159 (2022).Article 

Google Scholar 
Heneghan, R. F. et al. Prog. Oceanogr. 198, 102659 (2021).Article 

Google Scholar 
Polovina, J. J., Dunne, J. P., Woodworth, P. A. & Howell, E. A. ICES J. Mar. Sci. 68, 986–995 (2011).Article 

Google Scholar 
Cheung, W. W. L. et al. Fish Fish. 10, 235–251 (2009).Article 

Google Scholar 
Hazen, E. L. et al. Nat. Clim. Change 3, 234–238 (2013).Article 

Google Scholar 
Powers, R. P. & Jetz, W. Nat. Clim. Change 9, 323–329 (2019).Article 

Google Scholar 
Purves, D. et al. Nature 493, 295–297 (2013).CAS 
Article 

Google Scholar 
Hobday, A. J., Spillman, C. M., Paige Eveson, J. & Hartog, J. R. Fish. Oceanogr. 25, 45–56 (2016).Article 

Google Scholar 
Pons, M. et al. Proc. Natl Acad. Sci. USA 119, e2114508119 (2022).Article 

Google Scholar  More