Ecology

Carballido, J. L. et al. A new giant titanosaur sheds light on body mass evolution among sauropod dinosaurs. Proc. R. Soc. B Biol. Sci. 284, 20171219 (2017).
Google Scholar 
Hechenleitner, E. M. et al. Two Late Cretaceous sauropods reveal titanosaurian dispersal across South America. Commun. Biol. 3, 622 (2020).PubMed 
PubMed Central 

Google Scholar 
Otero, A., Carballido, J. L., Salgado, L., Canudo, I. & Garrido, A. C. Report of a giant titanosaur sauropod from the Upper Cretaceous of Neuquén Province, Argentina. Cretaceous Res. 122, 104754 (2021).
Google Scholar 
Sander, P. M. et al. Biology of the sauropod dinosaurs: The evolution of gigantism. Biol. Rev. 86, 117–155 (2011).PubMed 

Google Scholar 
Gorscak, E. & O’Connor, P. M. A new African titanosaurian sauropod dinosaur from the middle Cretaceous Galula Formation (Mtuka Member), Rukwa Rift Basin, southwestern Tanzania. PLoS One 14, e0211412 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Poropat, S. F. et al. Second specimen of the Late Cretaceous Australian sauropod dinosaur Diamantinasaurus matildae provides new anatomical information on the skull and neck of early titanosaurs. Zool. J. Linn. Soc. 192, 610–674 (2021).
Google Scholar 
González Riga, B. J. et al. An overview of the appendicular skeletal anatomy of South American titanosaurian sauropods, with definition of a newly recognized clade. An. Acad. Bras. Ciênc. 91, e20180374 (2019).
Google Scholar 
Gorscak, E. & O’connor, P. M. Time-calibrated models support congruency between Cretaceous continental rifting and titanosaurian evolutionary history. Biol. Lett. 12, 20151047 (2016).PubMed 
PubMed Central 

Google Scholar 
Pérez Moreno, A., Carballido, J. L., Otero, A., Salgado, L. & Calvo, J. O. The axial skeleton of Rinconsaurus caudamirus (Sauropoda: Titanosauria) from the Late Cretaceous of Patagonia, Argentina. Ameghiniana 59, 1–46 (2022).
Google Scholar 
Chiappe, L. M., Jackson, F., Coria, R. A. & Dingus, L. Nesting titanosaurs from Auca Mahuevo and adjacent sites. In The Sauropods: Evolution and Paleobiology (eds Curry Rogers, K. A. & Wilson, J. A.) 285–302 (University of California Press, 2005).
Google Scholar 
Hechenleitner, E. M., Grellet-Tinner, G. & Fiorelli, L. E. What do giant titanosaur dinosaurs and modern Australasian megapodes have in common?. PeerJ 3, e1341 (2015).PubMed 
PubMed Central 

Google Scholar 
Cojan, I., Renard, M. & Emmanuel, L. Palaeoenvironmental reconstruction of dinosaur nesting sites based on a geochemical approach to eggshells and associated palaeosols (Maastrichtian, Provence Basin, France). Palaeogeogr. Palaeoclimatol. Palaeoecol. 191, 111–138 (2003).
Google Scholar 
Grellet-Tinner, G., Codrea, V., Folie, A., Higa, A. & Smith, T. First evidence of reproductive adaptation to “island effect” of a dwarf Cretaceous Romanian titanosaur, with embryonic integument in ovo. PLoS One 7(3), e32051 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Khosla, A. & Lucas, S. Late Cretaceous dinosaur eggs and eggshells of peninsular India. Top. Geobiol. 51, 1–295 (2020).
Google Scholar 
Vila, B., Galobart, A., Oms, O., Poza, B. & Bravo, A. M. Assessing the nesting strategies of Late Cretaceous titanosaurs: 3-D clutch geometry from a new megaloolithid egg site. Lethaia 43, 197–208 (2009).
Google Scholar 
Chiappe, L. M. et al. Sauropod dinosaur embryos from the Late Cretaceous of Patagonia. Nature 396, 258–261 (1998).ADS 
CAS 

Google Scholar 
Salgado, L. et al. Upper Cretaceous dinosaur nesting sites of Río Negro (Salitral Ojo de Agua and Salinas de Trapalco-Salitral de Santa Rosa), Northern Patagonia, Argentina. Cretaceous Res. 28, 392–404 (2007).
Google Scholar 
Salgado, L., Magalhães Ribeiro, C. M., García, R. A. & Fernández, M. A. Late Cretaceous megaloolithid eggs from Salitral de Santa Rosa (Río Negro, Patagonia, Argentina): Inferences on the titanosaurian reproductive biology. Ameghiniana 46, 605–620 (2009).
Google Scholar 
Grellet-Tinner, G. & Fiorelli, L. E. A new Argentinean nesting site showing neosauropod dinosaur reproduction in a Cretaceous hydrothermal environment. Nat. Commun. 1, 32 (2010).ADS 
PubMed 

Google Scholar 
Hechenleitner, E. M. et al. A new Upper Cretaceous titanosaur nesting site from La Rioja (NW Argentina), with implications for titanosaur nesting strategies. Palaeontology 59, 433–446 (2016).
Google Scholar 
Kundrát, M. et al. Specialized craniofacial anatomy of a titanosaurian embryo from Argentina. Curr. Biol. 30, 4263–4269 (2020).PubMed 

Google Scholar 
Chiappe, L. M., Salgado, L. & Coria, R. A. Embryonic skulls of titanosaur sauropod dinosaurs. Science 293, 2444–2446 (2001).ADS 
CAS 
PubMed 

Google Scholar 
Grellet-Tinner, G., Chiappe, L. M., Norell, M. & Bottjer, D. Dinosaur eggs and nesting behaviors: A paleobiological investigation. Palaeogeogr. Palaeoclimatol. Palaeoecol. 232, 294–321 (2006).
Google Scholar 
Coria, R. A. & Chiappe, L. M. Embryonic skin from Late Cretaceous sauropods (Dinosauria) of Auca Mahuevo, Patagonia, Argentina. J. Paleontol. 81, 1528–1532 (2007).
Google Scholar 
Leuzinger, L. et al. Life and reproduction of titanosaurians: Isotopic hallmark of mid-palaeolatitude eggshells and its significance for body temperature, diet, and nesting. Chem. Geol. 583, 120452 (2021).ADS 
CAS 

Google Scholar 
Faccio, G. Dinosaurian eggs from the upper Cretaceous of Uruguay. In Dinosaur Eggs and Babies (eds Carpenter, K. et al.) 47–53 (Cambridge University Press, 1994).
Google Scholar 
Vianey-Liaud, M., Hirsch, K., Sahni, A. & Sigé, B. Late Cretaceous Peruvian eggshells and their relationships with Laurasian and Eastern Gondwanian material. Geobios 30, 75–90 (1997).
Google Scholar 
Magalhães Ribeiro, C. M. Microstructural analysis of dinosaur eggshells from Bauru Basin (Late Cretaceous), Minas Gerais, Brasil. In First International Symposium on Dinosaur Eggs and Babies (eds Bravo, A. M. & Reyes, T.) 117–122 (Isona I Conca Dellà Catalonia, 2000).
Google Scholar 
Magalhães Ribeiro, C. M. Ovo e fragmentos de cascas de ovos de dinossauros, provenientes de região de Peirópolis, Uberaba, Minas Gerais. Arq. Mus. Nac. 60, 223–228 (2002).
Google Scholar 
Grellet-Tinner, G. & Zaher, H. Taxonomic identification of the Megaloolithidae egg and eggshells from Cretaceous Bauru Basin (Minas Gerais, Brazil): Comparison with the Auca Mahuevo (Argentina) titanosaurid eggs. Papéis Avulsos Zool. 47, 105–112 (2007).
Google Scholar 
Martinelli, A. G. & Teixeira, V. P. A. The Late Cretaceous vertebrate record from the Bauru Group in the Triângulo Mineiro, southeastern Brazil. Bol. Geol. Minero 126, 129–158 (2015).
Google Scholar 
Soares, M. V. T. et al. Sedimentology of a distributive fluvial system: The Serra da Galga Formation, a new lithostratigraphic unit (Upper Cretaceous, Bauru Basin, Brazil). Geol. J. 56, 951–975 (2021).
Google Scholar 
Soares, M. V. T. et al. Landscape and depositional controls on palaeosols of a distributive fluvial system (Upper Cretaceous, Brazil). Sedim. Geol. 409, 105774 (2020).
Google Scholar 
Martinelli, A. G. et al. Palaeoecological implications of an Upper Cretaceous tetrapod burrow (Bauru Basin; Peirópolis, Minas Gerais, Brazil). Palaeogeogr. Palaeoclimatol. Palaeoecol. 528, 147–159 (2019).
Google Scholar 
Grellet-Tinner, G., Chiappe, L. M. & Coria, R. A. Eggs of titanosaurid sauropods from the Upper Cretaceous of Auca Mahuevo (Argentina). Can. J. Earth Sci. 41, 949–960 (2004).ADS 

Google Scholar 
Wilson, J. A., Mohabey, D. M., Peters, S. E. & Head, J. J. Predation upon hatchling dinosaurs by a new snake from the Late Cretaceous of India. PLoS Biol. 8, e1000322 (2010).PubMed 
PubMed Central 

Google Scholar 
Garcia, G. & Vianey-Liaud, M. Nouvelles données sur les coquilles d’œufs de dinosaures Megaloolithidae du Sud de la France: Systématique et variabilité intraspécifique. C. R. Acad. Sci. A. 332, 185–191 (2001).
Google Scholar 
Vianey-Liaud, M., Khosla, A. & Garcia, G. Relationships between European and Indian dinosaur eggs and eggshells of the oofamily Megaloolithidae. J. Vertebr. Paleontol. 23, 575–585 (2003).
Google Scholar 
Campos, D. A., Kellner, A. W. A., Bertini, R. J. & Santucci, R. M. On a titanosaurid (Dinosauria, Sauropoda) vertebral column from the Bauru Group, Late Cretaceous of Brazil. Arq. Mus. Nac. 63, 565–593 (2005).
Google Scholar 
Kellner, A. W. A., Campos, D. A. & Trotta, M. N. F. Description of a titanosaurid caudal series from the Bauru Group, Late Cretaceous of Brazil. Arq. Mus. Nac. 63, 529–564 (2005).
Google Scholar 
Machado, E. B., Avilla, L. S., Nava, W. R., Campos, D. A. & Kellner, A. W. A. A new titanosaur sauropod from the Late Cretaceous of Brazil. Zootaxa 3701, 301–321 (2013).PubMed 

Google Scholar 
Silva Junior, J. C. G. et al. Reassessment of Aeolosaurus maximus, a titanosaur dinosaur from the Late Cretaceous of Southeastern Brazil. Hist. Biol. 34, 402–411 (2022).
Google Scholar 
Bandeira, K. L. N. et al. A new giant Titanosauria (Dinosauria: Sauropoda) from the Late Cretaceous Bauru Group, Brazil. PLoS One 11, e0163373 (2016).PubMed 
PubMed Central 

Google Scholar 
Pearce, J. M. Philopatry: A return to origins. Auk 124, 1085–1087 (2007).
Google Scholar 
Kokko, H. & López-Sepulcre, A. From individual dispersal to species ranges: Perspectives for a changing world. Science 313, 789–791 (2006).ADS 
CAS 
PubMed 

Google Scholar 
Basilici, G., Hechenleitner, E. M., Fiorelli, L. E., Dal Bó, P. F. & Mountney, N. P. Preservation of titanosaur egg clutches in Upper Cretaceous cumulative palaeosols (Los Llanos Formation, La Rioja, Argentina). Palaeogeogr. Palaeoclimatol. Palaeoecol. 482, 83–102 (2017).
Google Scholar 
Mueller-Töwe, I. J., Sander, P. M., Schüller, H. & Thies, D. Hatching and infilling of dinosaur eggs as revealed by computed tomography. Palaeontogr. Abt. A 267, 119–168 (2002).
Google Scholar 
Paganelli, C. V. The physics of gas exchange across the avian eggshell. Am. Zool. 20, 329–338 (1980).
Google Scholar 
Grellet-Tinner, G., Fiorelli, L. E. & Salvador, R. B. Water vapor conductance of the Lower Cretaceous dinosaurian eggs from Sanagasta, La Rioja, Argentina: Paleobiological and paleoecological implications for South American faveoloolithid and megaloolithid eggs. Palaios 27, 35–47 (2012).ADS 

Google Scholar 
Grellet-Tinner, G. Phylogenetic interpretation of eggs and eggshells: Implications for phylogeny of Palaeognathae. Alcheringa 30, 141–182 (2006).
Google Scholar 
Mikhailov, K. E., Bray, E. S. & Hirsch, K. F. Parataxonomy of fossil egg remains (Veterovata): Principles and applications. J. Vertebr. Paleontol. 16, 763–769 (1996).
Google Scholar 
Nys, Y., Gautron, J., Garcia-Ruiz, J. M. & Hincke, M. T. Avian eggshell mineralization: Biochemical and functional characterization of matrix proteins. C. R. Palevol 3, 549–562 (2004).
Google Scholar 
Fernández, M. S., Passalacqua, K., Arias, J. I. & Arias, J. L. Partial biomimetic reconstitution of avian eggshell formation. J. Struct. Biol. 148, 1–10 (2004).PubMed 

Google Scholar 
Arias, J. L. & Fernández, M. S. Biomimetic processes through the study of mineralized shells. Mater. Charact. 50, 189–195 (2003).CAS 

Google Scholar 
Arias, J. L., Mann, K., Nys, Y., Garcia Ruiz, J. M. & Fernández, M. S. Eggshell growth and matrix macromolecules. In Handbook of Biomineralization (ed. Baeuerlein, E.) 309–328 (Wiley, 2007).
Google Scholar 
Chien, Y.-C., Hincke, M. T., Vali, H. & Mckee, M. D. Ultrastructural matrix–mineral relationships in avian eggshell, and effects of osteopontin on calcite growth in vitro. J. Struct. Biol. 163, 84–99 (2008).CAS 
PubMed 

Google Scholar 
Hincke, M. T. et al. The eggshell: Structure, composition and mineralization. Front. Biosci. 17, 1266–1280 (2012).CAS 

Google Scholar 
Stapane, L. et al. Avian eggshell formation reveals a new paradigm for vertebrate mineralization via vesicular amorphous calcium carbonate. J. Biol. Chem. 295, 15853–15869 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Norell, M. A. et al. A theropod dinosaur embryo and the affinities of the Flaming Cliffs dinosaur eggs. Science 266, 779–782 (1994).ADS 
CAS 
PubMed 

Google Scholar 
Araújo, R. et al. Filling the gaps of dinosaur eggshell phylogeny: Late Jurassic Theropod clutch with embryos from Portugal. Sci. Rep. 3, 1924 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Xing, L. et al. An exquisitely preserved in-ovo theropod dinosaur embryo sheds light on avian-like prehatching postures. iScience 103516, 24 (2021).
Google Scholar 
Sato, T., Cheng, Y.-N., Wu, X.-C., Zelenitsky, D. K. & Hsiao, Y.-F. A pair of shelled eggs inside a female dinosaur. Science 308, 375 (2005).CAS 
PubMed 

Google Scholar 
Norell, M. A., Clark, J. M. & Chiappe, L. M. A nesting dinosaur. Nature 378, 774–776 (1995).ADS 
CAS 

Google Scholar 
Dong, Z. M. & Currie, P. J. On the discovery of an oviraptorid skeleton on a nest of eggs at Bayan Mandahu, Inner Mongolia, People’s Republic of China. Can. J. Earth Sci. 33, 631–636 (1996).ADS 

Google Scholar  More

Water ecological security evaluation results of Yellow River BasinIndex weight analysisThis study selects the index weights in 2009, 2014, and 2019 for comparative analysis. As shown in Table 3, in terms of space, in the pressure layer, indicator A6 (Water area) has the most prominent weight, and indicator A3 (Natural population growth rate) has the most negligible weight; in the state layer, indicator B6 (Proportion of wetland area to total area) has the most prominent weight, and B1 (COD emissions per 10,000 yuan GDP) has the most negligible weight; in the response layer, indicator C3 (Green area rate of built-up area) has the most prominent weight, and indicator C2 (Centralized treatment rate of urban domestic sewage) has the most negligible weight. In summary, water area, wetland area, and built-up green space are the key indicators affecting the water ecology of the Yellow River Basin, including natural factors and economic and social factors.Table 3 Water ecological security index weight.Full size tableIn terms of time, indicators A6 and B6 have equal weights in three years and have always been in an important position. The weight of indicator C1 (the rate of stable compliance of wastewater discharge by industrial enterprises) has fallen for three consecutive years, from 0.38 to 0.09. It shows that after years of environmental management in various cities, the rate of compliance with wastewater discharge standards of industrial enterprises has been continuously increasing. It plays a positive role in the construction of water ecological security. The weight of indicator C3 has increased significantly in three years, from 0.31 in 2009 to 0.90 in 2019, indicating that with the continuous development of urbanization, the built-up area has become larger and larger, which has a massive impact on water ecological security. Therefore, the green area in the built-up area is vital, which is the key to ensuring the urban ecological environment. It is also a critical factor in maintaining the water ecological security.Trend analysis of water ecological securityThis study is based on Eq. (4) to calculate the WESI of the nine provinces in the past ten years, as shown in Fig. 3. From the perspective of the changes in WESI from 2009 to 2019, the overall trend is slowly increasing. Compared with 2009, WESI increased by 5.96% in 2019, but the average annual growth rate was only 0.59%. The sharp rise stage was in 2009–2012, with an average annual growth rate of 1.84%. Since 2009, there has been no inferior V water in the main stream of the Yellow River, and the water quality has been improving year by year. During this period, the nine provinces implemented the Yellow River Basin Flood Control Plan under the guidance of The State Council. The plan calls for strengthening infrastructure construction in the Yellow River Basin and conducting work such as river improvement and soil and water conservation. Therefore, we will promote the restoration of water ecology in the river basin and improve the safety of water ecology. From 2012 to 2019, WESI showed a trend of ups and downs. This is because the provinces have gradually shifted their development focus to the economy after achieving significant results in restoring water ecology in the river basin. The rapid economic development has brought more significant pressure to environmental governance and hindered water ecological safety improvement.Figure 3Trend map of water ecological security index (WESI) of nine provinces.Full size imageCriterion layer quantitative resultsTo further study and appraise the water ecological security of the study area, this paper quantifies the criteria layers (i.e., pressure, state, response) on account of the SMI-P method. It selects 2009, 2014, and 2019 for comparative analysis. As shown in Fig. 4, the criterion layer has undergone specific changes over time. First of all, the distribution of pressure in 62 cities has not changed much in three years. The areas with more tremendous pressure on water ecological security are mainly concentrated in eastern cities, including Shuozhou, Taiyuan, Jinzhou, Luliang, Linfen, Jincheng, and Changzhi, Anyang, Hebi, Jiaozuo, Puyang, Liaocheng, and other cities. Areas with less pressure are mainly concentrated in western and eastern cities, including Guoluo Tibetan Autonomous Prefecture, Hainan Tibetan Autonomous Prefecture, Haibei Tibetan Autonomous Prefecture, Ordos, Bayannaoer, Yulin, and other cities. In 2009, the precipitation in spring and winter in Lanzhou is less, the degree of drought is serious, and the flood disaster is more severe in flood season, which brings tremendous pressure to the water ecological security. After 2015, Lanzhou continued to implement the Action Plan for Prevention and Control of Water Pollution and then the river chief system was implemented. In 2019, The Work Plan of Lanzhou Municipal Water Pollution Prevention and Control Action in 2019 was issued and implemented. All these measures and actions have laid a foundation for water ecological security. On the contrary, with the rapid development of urbanization and economy and society, the pressure of water ecological security in Jinan has increased.Figure 4Quantitative spatial distribution map of the 62 cities in the Yellow River Basin. Note This was created by ArcMap-GIS, version 10.5. https://www.esri.com/.Full size imageThe larger the value of the status layer, the better the aquatic ecological status. On the contrary, the worse the aquatic ecological security. The overall spatial distribution of the status layer has not changed significantly in the past three years, and the changes are mainly concentrated in some cities. For example, the water ecological security status of Wuhan and Ulan Chab has gradually deteriorated in three years. The reason is that the urban population is becoming denser and sewage discharge is increasing, but related management and measures have not been fully implemented. In Dongying, the water ecological security status improved in 2014 and 2019. According to the Environmental Status Bulletin, in 2014, Dongying deepened its drainage basin pollution control system, continuously strengthened the restraint mechanism to improve river water quality, and carried out a pilot wetland ecological restoration.In the three years of 2009, 2014, and 2019, the response layer has changed more significantly than the pressure and status layers. It can be seen that the degree of response scarcity has gradually shifted from western cities to eastern cities. The reason can be understood as that due to their superior natural conditions, western cities have relatively weak awareness of water ecological protection and governance, and their ability to respond to emergencies is insufficient. However, with the increasingly prominent ecological and environmental problems, the awareness of maintaining water ecological safety is increasing, and the protection and governance measures are constantly improving. For example, Guoluo Tibetan Autonomous Prefecture, Hainan Tibetan Autonomous Prefecture, and Haibei Tibetan Autonomous Prefecture. Eastern cities are densely populated, urbanization development is faster than western cities, and environmental problems occur more frequently. Therefore, the awareness of ecological and environmental protection is more substantial, the governance system is relatively complete, and responsiveness is relatively good. However, as time progresses, some cities have somewhat slackened their ecological environment governance, and therefore their responsiveness has also weakened. For example, Shuozhou, Jinzhou, Lvliang, Linfen, and other places.Final quantitative resultsIn order to show the water ecological security status of 62 cities more intuitively, this paper shows the water ecological status level in Table 2 through the GIS spatial distribution map (Fig. 5).Figure 5Distribution map of water ecological security status in 62 cities of the Yellow River Basin. Note This was created by ArcMap-GIS, version 10.5. https://www.esri.com/.Full size imageLooking at the overall situation in the past three years, the water ecological security status is relatively stable, with little overall change. The reasons mainly include natural geographical location and economic and social development. In terms of physical geography, the safer areas are concentrated in the upper reaches of the Yellow River Basin, all of which have the characteristics of large land and sparsely populated areas and relatively superior natural conditions. They provide good conditions and foundations for the construction of water ecological security. The moderate warning cities are primarily located in the Loess Plateau and the North China Plain, where water resources are scarce, and the dense population, posing a threat to water ecological security. In terms of economic and social development, relatively safe areas are located in remote areas with inconvenient transportation. The region is dominated by agriculture and animal husbandry, with relatively backward economic development and a low level of urbanization. In addition, the threat to water ecological security is relatively tiny. Residents in the moderate warning area have a significant living demand, and the over-exploitation and utilization of natural resources have led to the destruction of the ecological environment. Therefore, it poses a more significant threat to water ecological security.Combining Fig. 5 and Table A.2 of appendix, it can be seen that in 2009, there were 8 safer cities, 22 with early warning level, and 32 with moderate warning. Relatively safe cities are concentrated in the southwest and north of the Yellow River Basin; cities with moderate warning level are distributed in the central and eastern areas. In 2014, the number of safer cities increased to 10, and the number of cities with moderate warning level decreased to 30. The means that water ecological security has received more and more attention, and cities have consciously strengthened the protection and governance of water ecology to maintain water ecological security. In 2019, there are 11 relatively safe cities, 21 cities with warning level, and 30 cities with moderate warning level. The overall situation has not changed much, and some cities have changed significantly. For example, Erdos had increased from an early warning status in 2009 to a safer status in 2014, and its safety index has risen from 0.57 to 0.65. Wuzhong has been upgraded from the warning level in 2009 (0.39) to the relatively safe in 2014 (0.44), and the safety index (0.47) in 2019 has also increased. Binzhou had improved from its early warning status (0.60) in 2009 to a relatively safe level (0.64) in 2014, and its safety index (0.66) has also increased in 2019, but the increase is not significant. On the contrary, Jinan has deteriorated from the early warning level in 2009 and 2014 to the moderate warning level in 2019, indicating that the water ecological security of Jinan has been seriously threatened in the process of rapid development.Spatial autocorrelation analysis of 62 cities in the Yellow River BasinGlobal spatial autocorrelation analysisThis paper selects 2009, 2014 and 2019, and analyzes the global spatial autocorrelation based on GeoDa. Combining Table 4 and Fig. 6, the Moran index for these three years was 0.298, 0.359, and 0.334 respectively, which were all in the [0,1] interval, indicating the water ecological security of 62 cities in the past three years showed significant spatial autocorrelation. Moreover, there is a positive spatial correlation, and the spatial autocorrelation is strong. The four quadrants of the scatter chart are high-high (i.e., first quadrant) aggregation area, low–high (i.e., second quadrant) aggregation area, low-low (i.e., third quadrant) aggregation area, and high-low (i.e., fourth quadrant) aggregation area. After testing, z-value  > 1.96, p-value  More

Detection of archival Phyllonorycter mines in historical herbariaOnly 1.5% (225 out of 15,009) of herbarium specimens of Tilia spp. examined from the Palearctic contained Ph. issikii leaf mines. These 225 herbarium specimens occurred in 185 geographical locations across the Palearctic, with the westernmost point in Germany (Hessen; the herbarium specimen dated by 2004) to the most eastern locations in Japan (on the island of Hokkaido; 1885–1974) (Fig. 1).Figure 1The localities where herbarium specimens of Tilia spp. carrying Phyllonorycter mines were collected in the Palearctic in the last 253 years. The dotted line divides Ph. issikii range to native (below the line) and invaded (above the line). The map was generated using ArcGIS 9.3 (Release 9.3. New York St., Redlands, CA. Environmental Systems Research Institute, http://www.esri.com/software/arcgis/eval-help/arcgis-93).Full size imageMost specimens with leaf mines (90%; 203/225) originated from Eastern Palearctic, in particular from the Russian Far East (RFE) (67.5%, 137/203) (Fig. 2a). In some cases, leaves were severely attacked, carrying up to 12 mines per leaf (as documented in the Russian Far East in 1930s–1960s). On the other hand, we found only 22 herbarium specimens with mines (10%; 22/225) from the putative invaded region in Western Palearctic, with the majority of herbarium specimens with mines (7% 15/225) from European Russia (Fig. 2b).Figure 2The presence of Phyllonorycter issikii mines in the herbarium specimens collected in the putative native (a) and invaded (b) ranges over the past 253 years (1764–2016). The number of herbarium specimens with and without mines and the percentage of the specimens with mines in each region or country from all herbarium specimens examined in a region or country (in brackets) are given next to each graph. The total number of herbarium specimens, including those with and without mines, is given for Eastern (a) and Western Palearctic (b) separately and altogether (a + b).Full size imageThe average number of leaf mines per herbarium specimen found in native (5.68 ± 0.77) and invaded regions (6.09 ± 1.70) was not significantly different (Mann–Whitney U-test: U = 20,145; Z = 0.43; p = 0.43). However, the infestation rate by Ph. issikii, i.e. percentage of leaves with mines per herbarium specimen was statistically higher in the West than in the East: 35% ± 8.19 versus 23% ± 1.94 (Mann–Whitney U-test: U = 1339; Z = 2.30; p = 0.02).Leaf mines from the East were significantly older than those from the West (Mann–Whitney U-test: U = 81; Z =  − 4.4; p  400 bp) were obtained for 71 archival specimens that were between 7 and 162 years old (Fig. 4, the points in dashed frame) (Table S4). Nine of these 71 specimens were over one century old (106–162-year-old): eight originated from the Palearctic and one from the Nearctic (Fig. 4, the points in gray cloud).In the Palearctic, the oldest successfully DNA barcoded Ph. issikii specimen (obtained sequence length 408 bp) was a 162-year-old larva dissected from the leaf mine on Tilia amurensis from the RFE (village Busse, Amur Oblast, the year 1859), sequence ID LMINH119-19 (Fig. 5, Table S5). In the Nearctic, the oldest sequenced specimen (obtained sequence length 658 bp) was 127-year-old larva of Ph. tiliacella on T. americana from USA, Pennsylvania (Fig. 5, Table S5).Figure 5A maximum likelihood tree of 81 COI sequences of Phyllonorycter spp. Overall, 71 archival sequenced specimens were dissected from herbaria collected in the Palearctic and the Nearctic in 1859–2014 and ten specimens (highlighted in blue) originated from the modern range20. The tree was generated with the K2P nucleotide substitution model and bootstrap method (2500 iterations), p  More

Graded lockdowns imposed by the South African government to manage the COVID-19 pandemic27,28,29 has afforded us a unique opportunity to quantify shorebird responses to increasing human density in Muizenberg Beach over 8 months in 2020, including a 2-month period of virtual human exclusion. In spite of our study being limited to one beach over 2 years, we were able to take advantage of data collected prior to- (2019) and during the 2020 COVID lockdowns, to better understand a pervasive feature of sandy beach ecosystems (human recreation) that is predicted to intensify in future10.Findings for the 2019–2020 component of our study generally conformed to hypotheses posed. Firstly, shorebird abundance was inversely associated with human abundance and was positively related to lockdown level in 2020. Secondly, shorebird abundance was generally greatest during lockdown levels 5 and 4, when humans were effectively absent from the beach. To contextualise, shorebird abundance was roughly six times greater at the start of lockdown level 5 (2020) than the equivalent period in 2019. Thirdly, lowest shorebird abundance occurred during lockdown level 1 when human abundance was greatest in 2020. Collectively, these findings indicate a strong inverse association between shorebird- and human abundance on Muizenberg Beach and align with results of other studies36,37,38,39. Cumulatively, our findings, allied with prior research highlight the potential for human recreational activity, particularly at high intensities, to impact shorebird utilisation of sandy beach ecosystems, which may in turn affect ecological functions they provide that contribute to ecosystem multifunctionality.The inverse relationship that we recorded between human- and shorebird abundance likely manifests through the diverse ways in which recreational activity impacts fundamental processes and ecosystem components, which in turn link ecologically to shorebirds10,36,37,38,39,40. Muizenberg Beach is popular for surfing, bait-harvesting and general recreational activities, and it is these activities that likely drive the human-shorebird relationship that we report, particularly in 2020. When carried out under high human densities, such activities can lead to a reduction in space available, rendering the ecosystem less suitable as a substrate for birds36. Noise pollution and the presence of dogs may further depress habitat suitability41. Repeated trampling of sediment can negatively impact macrofaunal populations, which together with altered sedimentary biogeochemistry (e.g. increased anoxia), can reduce trophic resource availability to shorebirds, with benthic bait-collecting compounding these effects42,43. At the start of our data collection in 2020, we were unable to identify shorebird species due to lockdown levels 5 and 4 prohibiting human presence on the beach27,28,29. It is probable though that shorebird assemblages during lockdown levels 5 and 4 were not the same as those we identified between lockdown level 3 to 1 (mainly gulls; Table 3). This is based on research showing that increasing environmental disturbances can induce switches in biotic assemblages to those that can tolerate human activities44. Thus, the shorebird assemblages we identified during lockdown levels 3 to 1 is potentially the end-result of the mechanisms highlighted above (space reduction, noise, reduced resource availability) acting on shorebird assemblages in the absence of humans (lockdown levels 5 and 4) following humans being permitted onto the beach.At an inter-annual level, our data revealed idiosyncratic patterns that raise interesting questions about human-shorebird relationships. In 2019, in the absence of any lockdowns, shorebird abundance rose over the winter period (May–August). Winter peaks in abundance have previously been recorded in the literature45,46,47, including for kelp gulls (Larus dominicanus), which were the dominant shorebird in Muizenberg Beach. Specifically, winter abundance peaks for this species have been recorded in sandy beaches in the Eastern Cape, the Swartkops Estuary and Algoa Bay in South Africa (southeast coast)45,46,47. However, the absence of a winter abundance peak in 2020 raises the possibility that the 2019 winter-peak was not seasonal but an opportunistic response to decreased human abundance (see Fig. 4A). In South Africa, coastal ecosystems generally experience greatest human numbers in summer, due to warmer conditions and long end-of-year-vacation periods, based on our observations and experiences.The second inter-annual trend worth noting in our findings is that shorebird abundance was greater in 2019 than 2020, despite lockdowns being implemented in 2020. This counterintuitive finding is likely due to lockdowns that excluded people from the beach in 2020 (levels 5 to 3) being too short in duration to facilitate increases in bird numbers in 2020 beyond the 2019 level. This is supported by our data showing that humans were excluded from the beach for a total of 2 months (April and May 2020; levels 5-4) out of the 8-month period during which photographs were analysed. It would have been expected at the onset of the study that humans would be excluded from the beach during lockdown level 329, which would have resulted in an additional two and a half months of human exclusion and potentially a higher mean shorebird abundance for 2020. However, it is clear from our data that humans were present on the beach during level 3. On closer inspection, it is evident that human numbers increased even prior to the end of lockdown level 4. In fact, human abundance was greater under lockdown level 3 in 2020 than in the same period in 2019. Such high numbers of humans on the beach despite prohibitions are likely due to a lack of compliance, confusion around regulations and/or ‘covid fatigue’, which describes the propensity of humans to grow tired of COVID-19 regulations48. An additional consideration is that human numbers on the beach increased dramatically during lockdown levels 2 and 1, being almost twice the level recorded in 2019 in the same period. The lower 2020 bird count that we recorded is thus likely a product of the short duration of human exclusions in 2020 (lockdown levels 4 and 5) and the magnitude and rate of increase in human numbers thereafter (levels 3-1). Separately, our findings additionally suggest that surrogates (lockdown levels in our case) are unreliable estimators of human presence or abundance and align with findings elsewhere24.The last noteworthy inter-annual trend in our data was the difference in strength of human-shorebird relationships. While the inverse relationship between human and shorebird numbers was evident in both years, it was only during 2020, when humans were excluded from Muizenberg Beach, that the extent of this relationship was revealed. Specifically, in 2020, human exclusion at the start of lockdown level 5 was accompanied by a six-fold increase in shorebird abundance relative to 2019 at the same period. Additional support for the difference in strength of the human-shorebird relationship is the (1) significant interaction recorded between human numbers and year in explaining shorebird abundance and (2) the almost twofold stronger negative relationship (based on regression slopes) between shorebird and human abundance in 2020 vs 2019. These findings suggest that were it not for the COVID lockdowns in 2020, the extent of increasing human numbers on shorebirds may have been masked. However, it must be borne in mind that inter-annual variation may have played some role in the difference in trends recorded for 2019 versus 2020, though we cannot quantify this, given that we only have data for 2 years. Nevertheless, we suggest that when making conservation/management recommendations, decision-makers need to be cognisant of the potential for human effects on sandy beach ecosystems to be underestimated in studies based on variation in human density, in which human exclusion at appropriate spatial and temporal scales is absent24. Concerns have been expressed in the past about the failure of studies to consistently detect large-scale changes in sandy beach ecosystems, including those induced by recreational activities19. We suggest that such deficiencies may relate in part to the scarcity of true human exclusions in disturbance studies at meaningful scales in space and time.Findings from the in situ component of our study suggested that shorebird assemblages were negligibly affected by the transition from lockdown level 3 to 1, but that spatial differences among zones were more prominent. The lack of cases in which lockdown levels interacted statistically with zones (Tables 2, 4) further reinforces our conclusion regarding lockdown effects. Shorebird assemblage structure did vary between lockdown levels 3 and 2, due mainly to increasing contributions of Chroicocephalus hartlaubii (Hartlaub’s Gull) from level 3 to 2 and the opposite for L. dominicanus. Contrary to our hypothesis, differences in assemblage (Shannon–Wiener diversity was the exception) and species metrics were not detected among lockdown levels. This was likely due to the gradient in human abundance being weak among lockdown levels 3 to 1, relative to levels 5 and 4, with there being no virtual exclusion of humans under level 3 lockdown, as would have been expected given government regulations29. It is also possible that under lockdown levels 3, 2 and 1, the shorebird assemblage was simplified and comprised species tolerant of human activities44. The increase in Shannon–Wiener diversity value from lockdown level 3 to 2 was counter expectation, but likely reflects increased evenness during lockdown level 2, brought on by the declining dominance of L. dominicanus and a greater contribution of C. hartlaubii.Taken in its entirety, our findings provide valuable perspectives on human-shorebird interactions in sandy beaches. Based on our 2020 data spanning lockdowns of decreasing severity, our findings suggest that shorebirds are likely to benefit from human-free periods. This benefit is in reality likely to extend across multiple-trophic levels and is unlikely to be shorebird-specific, based on prior research reporting positive organism metrics at lower trophic levels in low human and/or human-free conditions in beach ecosystems20. Broadly, our findings attest to the value of using current and future lockdowns associated with managing the global COVID-19 pandemic to provide data on responses of birds and other organism groups to human-free spaces and times25,26,49. These human-free conditions can additionally provide invaluable data on sensitivities of ecosystem components and processes to increasing human density25,26,49. Data collected during lockdowns can provide better approximations of baseline conditions in sandy beach ecosystems, thereby providing a more meaningful basis for (1) evaluating future ecosystem change in response to human and global change stressors and (2) developing ecosystem restoration programs. This would be central to preventing long-term ecosystem degradation through the shifting base-line syndrome, where successive generations of decision makers/scientists judge the magnitude of change experienced by ecosystem components against increasingly deteriorating conditions over generational time-scales50. We also advocate for data emanating from COVID lockdown studies to be used in public education initiatives, so that beach users are made aware of the ways in which recreational activities can influence beach ecosystems. Such initiatives can improve involvement of public stakeholders in management of sandy beach ecosystems, which has been shown to provide cost-effective and sound decision-making, while increasing support for conservation initiatives51,52,53.Lastly, our findings have shed light on the sensitivity of shorebirds to increasing human numbers, mainly for recreational purposes. By moving beyond binary contrasts of human presence/absence, our work has also shown the magnitude of increasing human numbers on shorebirds, by virtue of the 34.18% increase in human abundance in our study corresponding with a 79.63% decline in bird numbers during the transition from lockdown level 4 to 3 in 2020. This finding is highly relevant considering that our work was based on an urban ecosystem—such systems are thought to have avian communities that are more disturbance tolerant relative to rural or suburban ecosystems54. Broadly, our work emphasises the need for environmental managers and city planners to be cognisant of the sensitivity of shorebirds to human recreational activities, even in urban settings, and to develop appropriate management plans in conjunction with scientists and stakeholders51,52,53. It should be noted that bird responses that we recorded in 2020 are unlikely to be driven solely by changing human numbers in Muizenberg Beach. Processes influencing bird assemblages in beaches surrounding our focal study area, including changes in human numbers and behaviour, may also have been influential determinants of trends recorded. We lack the data to comment meaningfully on this, but is an area worth exploring in future studies.Concluding perspectivesThe global COVID-19 anthropause has been described as the greatest large-scale experiment in modern history. This period has afforded scientists a unique opportunity to refine understanding of the consequences of human activities on Earth’s natural environments25,26,49. This is particularly relevant for human-dominated ecosystems such as sandy beaches, which are arguably the most utilised of Earth’s ecosystems for recreational purposes. In the absence of the COVID-19 anthropause, it is doubtful whether human exclusions could be carried out at scales that would allow meaningful detection of responses to human recreational disturbance. Our findings broadly attest to the points raised thus far, illustrating not only the potential for conventional approaches to underestimate human effects in sandy beaches, but also the sensitivity of shorebirds to human recreation and the magnitude of human influence. We hope that our findings stimulate further research on human recreational effects on sandy beach ecosystems, particularly with a view towards quantifying disturbance sensitivities and response thresholds of fundamental processes that drive multifunctionality in these heavily utilised, yet highly significant coastal ecosystems. We suggest that this is an imperative, given the exponential human population growth expected in the future, particularly along the coast, and the increasing demand predicted on sandy beach ecosystems from recreation, tourism and commercial sectors10,18. At its broadest level, our work dovetails with prior calls for scientists to capitalise on current and future COVID lockdowns to refine our understanding of human-nature interactions25, so that ecosystems and socio-ecological services provided can be sustainably utilised in the future. More

Whittaker, R. H. Vegetation of the Siskiyou mountains, Oregon and California. Ecol. Monogr. 30, 279–338 (1960).
Google Scholar 
Wilson, R. J., Thomas, C. D., Fox, R., Roy, D. B. & Kunin, W. E. Spatial patterns in species distributions reveal biodiversity change. Nature 432, 393–396 (2004).CAS 
PubMed 

Google Scholar 
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).CAS 
PubMed 

Google Scholar 
Ings, T. C. et al. Ecological networks—beyond food webs. J. Anim. Ecol. 78, 253–269 (2009).PubMed 

Google Scholar 
Dunne, J. A. & Williams, R. J. Cascading extinctions and community collapse in model food webs. Philos. Trans. R. Soc. Lond. B: Biol. Sci. 364, 1711–1723 (2009).
Google Scholar 
Leclère, D. et al. Bending the curve of terrestrial biodiversity needs an integrated strategy. Nature 585, 551–556 (2020).PubMed 

Google Scholar 
Vellend, M. The Theory of Ecological Communities Vol. 57 229 (Princeton University Press, 2016).Altermatt, F. Diversity in riverine metacommunities: a network perspective. Aquat. Ecol. 47, 365–377 (2013).
Google Scholar 
Peterson, E. E. et al. Modelling dendritic ecological networks in space: an integrated network perspective. Ecol. Lett. 16, 707–719 (2013).PubMed 

Google Scholar 
Tonkin, J. D. et al. The role of dispersal in river network metacommunities: patterns, processes, and pathways. Freshw. Biol. 63, 141–163 (2018).
Google Scholar 
Muneepeerakul, R. et al. Neutral metacommunity models predict fish diversity patterns in Mississippi-Missouri basin. Nature 453, 220–222 (2008).CAS 
PubMed 

Google Scholar 
Besemer, K. et al. Headwaters are critical reservoirs of microbial diversity for fluvial networks. Proc. Biol. Sci. 280, 20131760 (2013).PubMed 
PubMed Central 

Google Scholar 
Finn, D. S., Bonada, N., Múrria, C. & Hughes, J. M. Small but mighty: headwaters are vital to stream network biodiversity at two levels of organization. J. North Am. Benthol. Soc. 30, 963–980 (2011).
Google Scholar 
Altermatt, F., Seymour, M. & Martinez, N. River network properties shape α-diversity and community similarity patterns of aquatic insect communities across major drainage basins. J. Biogeogr. 40, 2249–2260 (2013).
Google Scholar 
Harvey, E., Gounand, I., Fronhofer, E. A. & Altermatt, F. Disturbance reverses classic biodiversity predictions in river-like landscapes. Proc. R. Soc. B: Biol. Sci. 285, 20182441 (2018).
Google Scholar 
Tylianakis, J. M., Laliberté, E., Nielsen, A. & Bascompte, J. Conservation of species interaction networks. Biol. Conserv. 143, 2270–2279 (2010).
Google Scholar 
Thompson, R. M. et al. Food webs: reconciling the structure and function of biodiversity. Trends Ecol. Evol. 27, 689–697 (2012).PubMed 

Google Scholar 
Woodward, G. & Hildrew, A. G. Food web structure in riverine landscapes. Freshw. Biol. 47, 777–798 (2002).
Google Scholar 
Williams, R. J. & Martinez, N. D. Limits to trophic levels and omnivory in complex food webs: theory and data. Am. Nat. 163, 458–468 (2004).PubMed 

Google Scholar 
Thompson, R. M. & Townsend, C. R. The effect of seasonal variation on the community structure and food-web attributes of two streams: implications for food-web science. Oikos 87, 75–88 (1999).
Google Scholar 
Wood, S. A., Russell, R., Hanson, D., Williams, R. J. & Dunne, J. A. Effects of spatial scale of sampling on food web structure. Ecol. Evol. 5, 3769–3782 (2015).PubMed 
PubMed Central 

Google Scholar 
Tylianakis, J. M. & Morris, R. J. Ecological networks across environmental gradients. Annu. Rev. Ecol., Evolution, Syst. 48, 25–48 (2017).
Google Scholar 
Romanuk, T. N. et al. The structure of food webs along river networks. Ecography 29, 3–10 (2006).
Google Scholar 
Olivier, P. et al. Exploring the temporal variability of a food web using long‐term biomonitoring data. Ecography 42, 2107–2121 (2019).
Google Scholar 
Poisot, T., Canard, E., Mouillot, D., Mouquet, N. & Gravel, D. The dissimilarity of species interaction networks. Ecol. Lett. 15, 1353–1361 (2012).PubMed 

Google Scholar 
Delmas, E. et al. Analysing ecological networks of species interactions. Biol. Rev. Camb. Philos. Soc. https://doi.org/10.1111/brv.12433 (2018).Article 
PubMed 

Google Scholar 
Tavares-Cromar, A. F. & Williams, D. D. The importance of temporal resolution in food web analysis: Evidence from a detritus-based stream. Ecol. Monogr. 66, 91–113 (1996).
Google Scholar 
Poisot, T., Stouffer, D. B. & Gravel, D. Beyond species: why ecological interaction networks vary through space and time. Oikos 124, 243–251 (2015).
Google Scholar 
Thomsen, P. F. & Willerslev, E. Environmental DNA—an emerging tool in conservation for monitoring past and present biodiversity. Biol. Conserv. 183, 4–18 (2015).
Google Scholar 
Deiner, K., Fronhofer, E. A., Mächler, E., Walser, J.-C. & Altermatt, F. Environmental DNA reveals that rivers are conveyer belts of biodiversity information. Nat. Commun. 7, 12544 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Dunne, J. A. In Ecological Networks: Linking Structure and Dynamics (eds. Pascual, J. A. & Dunne, J. A.) 27–86 (University Press, 2006).Neff, F. et al. Changes in plant-herbivore network structure and robustness along land-use intensity gradients in grasslands and forests. Sci Adv 7, eabf3985 (2021).O’Connor, M. J. et al. Unveiling the food webs of tetrapods across Europe through the prism of the Eltonian niche. J. Biogeogr. 47, 181–192 (2020).
Google Scholar 
Pellissier, L. et al. Comparing species interaction networks along environmental gradients. Biol. Rev. Camb. Philos. Soc. 93, 785–800 (2018).PubMed 

Google Scholar 
Saravia, L. A. et al. Ecological network assembly: how the regional metaweb influences local food webs. BioRxiv, https://doi.org/10.1101/340430 (2021).Blackman, R. C. et al. Mapping biodiversity hotspots of fish communities in subtropical streams through environmental DNA. Sci. Rep. 4, e65352 (2021).
Google Scholar 
Baselga, A. & Orme, C. D. L. betapart: an R package for the study of beta diversity: Betapart package. Methods Ecol. Evol. 3, 808–812 (2012).
Google Scholar 
Seymour, M. et al. Executing multi-taxa eDNA ecological assessment via traditional metrics and interactive networks. Sci. Total Environ. 729, 138801 (2020).CAS 
PubMed 

Google Scholar 
D’Alessandro, S. & Mariani, S. Sifting environmental DNA metabarcoding data sets for rapid reconstruction of marine food webs. Fish Fish 22, 822–833 (2021).
Google Scholar 
Zhang, Y. et al. Holistic pelagic biodiversity monitoring of the Black Sea via eDNA metabarcoding approach: From bacteria to marine mammals. Environ. Int. 135, 105307 (2020).PubMed 

Google Scholar 
Altermatt, F. et al. Uncovering the complete biodiversity structure in spatial networks: the example of riverine systems. Oikos 129, 607–618 (2020).
Google Scholar 
Widder, S. et al. Fluvial network organization imprints on microbial co-occurrence networks. Proc. Natl Acad. Sci. USA 111, 12799–12804 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Seymour, M. et al. Environmental DNA provides higher resolution assessment of riverine biodiversity and ecosystem function via spatio-temporal nestedness and turnover partitioning. Commun. Biol. 4, 512 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Mächler, E. et al. Assessing different components of diversity across a river network using eDNA. Environ. DNA 1, 290–301 (2019).
Google Scholar 
Peralta-Maraver, I., López-Rodríguez, M. J. & de Figueroa, J. M. T. Structure, dynamics and stability of a Mediterranean river food web. Mar. Freshw. Res. 68, 484–495 (2017).
Google Scholar 
Woodward, G. et al. Ecological networks in a changing climate. Ecol. Netw. 42, 71–138 (2010).
Google Scholar 
Kondoh, M., Kato, S. & Sakato, Y. Food webs are built up with nested subwebs. Ecology 91, 3123–3130 (2010).PubMed 

Google Scholar 
Vannote, R. L., Minshall, G. W., Cummins, K. W., Sedell, J. R. & Cushing, C. E. The River Continuum Concept. Can. J. Fish. Aquat. Sci. 37, 130–137 (1980).
Google Scholar 
Power, M. E. & Dietrich, W. E. Food webs in river networks. Ecol. Res. https://doi.org/10.1046/j.0912-3814.2002.00503.x (2002).Montoya, D., Yallop, M. L. & Memmott, J. Functional group diversity increases with modularity in complex food webs. Nat. Commun. 6, 7379 (2015).CAS 
PubMed 

Google Scholar 
Gravel, D., Albouy, C. & Thuiller, W. The meaning of functional trait composition of food webs for ecosystem functioning. Philos. Trans. R. Soc. Lond. B: Biol. Sci. 371, 20150268 (2016).Ruppert, K. M., Kline, R. J. & Rahman, M. S. Past, present, and future perspectives of environmental DNA (eDNA) metabarcoding: a systematic review in methods, monitoring, and applications of global eDNA. Glob. Ecol. Conserv. 17, e00547 (2019).
Google Scholar 
Carraro, L., Mächler, E., Wüthrich, R. & Altermatt, F. Environmental DNA allows upscaling spatial patterns of biodiversity in freshwater ecosystems. Nat. Commun. 11, 3585 (2020).PubMed 
PubMed Central 

Google Scholar 
Barnes, M. A. & Turner, C. R. The ecology of environmental DNA and implications for conservation genetics. Conserv. Genet. 17, 1–17 (2016).CAS 

Google Scholar 
Bista, I. et al. Annual time-series analysis of aqueous eDNA reveals ecologically relevant dynamics of lake ecosystem biodiversity. Nat. Commun. 8, 14087 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Erickson, R. A., Merkes, C. M., Jackson, C. A., Goforth, R. R. & Amberg, J. J. Seasonal trends in eDNA detection and occupancy of bigheaded carps. J. Gt. Lakes Res. 43, 762–770 (2017).
Google Scholar 
Troth, C. R., Sweet, M. J., Nightingale, J. & Burian, A. Seasonality, DNA degradation and spatial heterogeneity as drivers of eDNA detection dynamics. Sci. Total Environ. 768, 144466 (2021).CAS 
PubMed 

Google Scholar 
Thalinger, B. et al. The effect of activity, energy use, and species identity on environmental DNA shedding of freshwater fish. Front. Ecol. Evolution 9, 73 (2021).
Google Scholar 
Kelly, R. P., Port, J. A., Yamahara, K. M. & Crowder, L. B. Using environmental DNA to census marine fishes in a large mesocosm. PLoS ONE 9, e86175 (2014).PubMed 
PubMed Central 

Google Scholar 
Leray, M. et al. A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Front. Zool. 10, 34 (2013).PubMed 
PubMed Central 

Google Scholar 
Geller, J., Meyer, C., Parker, M. & Hawk, H. Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Mol. Ecol. Resour. 13, 851–861 (2013).CAS 
PubMed 

Google Scholar 
Liu, C. M. et al. BactQuant: An enhanced broad-coverage bacterial quantitative real-time PCR assay. BMC Microbiol. 12, 56 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Mansfeldt, C. et al. Microbial community shifts in streams receiving treated wastewater effluent. Sci. Total Environ. 709, 135727 (2020).CAS 
PubMed 

Google Scholar 
Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010).CAS 
PubMed 

Google Scholar 
Andrews, S. FASTQC A Quality Control tool for High Throughput Sequence Data (Babraham Institute, 2015).Magoč, T. & Salzberg, S. L. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963 (2011).PubMed 
PubMed Central 

Google Scholar 
Hänfling, B. et al. Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods. Mol. Ecol. 25, 3101–3119 (2016).PubMed 

Google Scholar 
Csárdi, G. & Nepusz, T. The igraph software package for complex network research. Int. J. Complex Syst. 1695, 1–9 (2006).
Google Scholar 
Oksanen, J. et al. vegan: Community Ecology Package 2.5-6. https://CRAN.Rproject.org/package=vegan (2019).Tachet, H., Bournaud, M., Richoux, P. & Usseglio-Polatera, P. Invertébrés d’eau douce—systématique, biologie, écologie (CNRS Editions, 2010).Schmidt-Kloiber, A. & Hering, D. www.freshwaterecology.info—an online tool that unifies, standardises and codifies more than 20,000 European freshwater organisms and their ecological preferences. Ecol. Indic. 53, 271–282 (2015).Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D. & Bertilsson, S. A guide to the natural history of freshwater lake bacteria. Microbiol. Mol. Biol. Rev. 75, 14–49 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Fortuna, M. A. et al. Nestedness versus modularity in ecological networks: two sides of the same coin? J. Anim. Ecol. 79, 811–817 (2010).PubMed 

Google Scholar 
Johnson, S., Domínguez-García, V., Donetti, L. & Muñoz, M. A. Trophic coherence determines food-web stability. Proc. Natl Acad. Sci. USA 111, 17923–17928 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wootton, K. L. Omnivory and stability in freshwater habitats: Does theory match reality? Freshw. Biol. 62, 821–832 (2017).
Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest package: tests in linear mixed effects models. J. Stat. Softw., Artic. 82, 1–26 (2017).
Google Scholar 
Lenth, R. V. Estimated Marginal Means, aka Least-Squares Means [R package emmeans version 1.6.1] (2021).RStudio Team RStudio: Integrated development for R. RStudio, PBC, Boston, MA. R version 4.0.4 Retrieved from http://www.rstudio.com/ (2021) More

Jeltsch, F. et al. Integrating movement ecology with biodiversity research—Exploring new avenues to address spatiotemporal biodiversity dynamics. Mov. Ecol. 1, 6 (2013).PubMed 
PubMed Central 

Google Scholar 
Dingle, H. Migration: The Biology of Life on the Move Migration (Oxford University Press, 2014).Joly, K. et al. Longest terrestrial migrations and movements around the world. Sci. Rep. 9, 15333 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Lundberg, J. & Moberg, F. Mobile link organisms and ecosystem functioning: Implications for ecosystem resilience and management. Ecosystems 6, 0087–0098 (2003).
Google Scholar 
Bauer, S. & Hoye, B. J. Migratory animals couple biodiversity and ecosystem functioning worldwide. Science 344, 1242552 (2014).CAS 
PubMed 

Google Scholar 
Nifong, J. C., Layman, C. A. & Silliman, B. R. Size, sex and individual-level behaviour drive intrapopulation variation in cross-ecosystem foraging of a top-predator. J. Anim. Ecol. 84, 35–48 (2015).PubMed 

Google Scholar 
Giroux, M.-A. et al. Benefiting from a migratory prey: Spatio-temporal patterns in allochthonous subsidization of an arctic predator. J. Anim. Ecol. 81, 533–542 (2012).PubMed 

Google Scholar 
Allen, A. M. & Singh, N. J. Linking movement ecology with wildlife management and conservation. Front. Ecol. Evol. 3, 155 (2016).
Google Scholar 
Bunnefeld, N. et al. A model-driven approach to quantify migration patterns: Individual, regional and yearly differences. J. Anim. Ecol. 80, 466–476 (2011).PubMed 

Google Scholar 
Teitelbaum, C. S. & Mueller, T. Beyond migration: Causes and consequences of nomadic animal movements. Trends Ecol. Evol. 34, 569–581 (2019).PubMed 

Google Scholar 
Berg, J. E., Hebblewhite, M., St. Clair, C. C. & Merrill, E. H. Prevalence and mechanisms of partial migration in ungulates. Front. Ecol. Evol. 7, 325 (2019).
Google Scholar 
Avgar, T., Street, G. & Fryxell, J. M. On the adaptive benefits of mammal migration. Can. J. Zool. 92, 481–490 (2014).
Google Scholar 
Barbour, M. G. & Billings, W. D. North American Terrestrial Vegetation (Cambridge University Press, 2000).
Google Scholar 
Smith, S. L., Throop, J. & Lewkowicz, A. G. Recent changes in climate and permafrost temperatures at forested and polar desert sites in northern Canada. Can. J. Earth Sci. 49, 914–924 (2012).ADS 

Google Scholar 
Lévesque, E. Plant Distribution and Colonization in Extreme Polar Deserts, Ellesmere Island, Canada (University of Toronto, 1997).
Google Scholar 
Bliss, L. C., Svoboda, J. & Bliss, D. I. Polar deserts, their plant cover and plant production in the Canadian High Arctic. Holarctic Ecol. 7, 305–324 (1984).
Google Scholar 
Berteaux, D. et al. Effects of changing permafrost and snow conditions on tundra wildlife: Critical places and times. Arctic Sci. 3, 65–90 (2017).
Google Scholar 
Duchesne, D., Gauthier, G. & Berteaux, D. Habitat selection, reproduction and predation of wintering lemmings in the Arctic. Oecologia 167, 967–980 (2011).ADS 
PubMed 

Google Scholar 
Fuglei, E., Blanchet, M.-A., Unander, S., Ims, R. A. & Pedersen, Å. Ø. Hidden in the darkness of the Polar night: A first glimpse into winter migration of the Svalbard rock ptarmigan. Wildl. Biol. 2017, SP1 (2017).
Google Scholar 
Schmidt, N. M. et al. Ungulate movement in an extreme seasonal environment: Year-round movement patterns of high-arctic muskoxen. Wildl. Biol. 22, 253–267 (2016).
Google Scholar 
Berteaux, D. & Lai, S. Walking on water: Terrestrial mammal migrations in the warming Arctic. Anim. Migr. 8, 65–73 (2021).
Google Scholar 
Gnanadesikan, G. E., Pearse, W. D. & Shaw, A. K. Evolution of mammalian migrations for refuge, breeding, and food. Ecol. Evol. 7, 5891–5900 (2017).PubMed 
PubMed Central 

Google Scholar 
Best, T. L. & Henry, T. H. Lepus arcticus. Mamm. Species 1–9 (1994).Dalerum, F. et al. Exploring the diet of arctic wolves (Canis lupus arctos) at their northern range limit. Can. J. Zool. 96, 277–281 (2018).
Google Scholar 
Mech, L. D. Annual arctic wolf pack size related to arctic hare numbers. Arctic 60, 309–311 (2007).
Google Scholar 
Small, R. J., Keith, L. B. & Barta, R. M. Demographic responses of Arctic hares Lepus arcticus placed on two predominantly forested islands in Newfoundland. Ecography 15, 161–165 (1992).
Google Scholar 
Small, R. J., Keith, L. B. & Barta, R. M. Dispersion of introduced arctic hares (Lepus arcticus) on islands off Newfoundland’s south coast. Can. J. Zool. 69, 2618–2623 (1991).
Google Scholar 
Hearn, B. J., Keith, L. B. & Rongstad, O. J. Demography and ecology of the arctic hare (Lepus arcticus) in southwestern Newfoundland. Can. J. Zool. 65, 852–861 (1987).
Google Scholar 
Harper, F. The Mammals of Keewatin Vol. 12 (Miscellaneaous Publications, Museum of Natural History, University of Kansas, 1956).
Google Scholar 
Dalerum, F. et al. Spatial variation in Arctic hare (Lepus arcticus) populations around the Hall Basin. Polar Biol. 40, 2113–2118 (2017).
Google Scholar 
Fraser, K. C. et al. Tracking the conservation promise of movement ecology. Front. Ecol. Evol. 6, 150 (2018).
Google Scholar 
CAFF. Arctic Biodiversity Assessment. Status and trends in Arctic biodiversity. Conservation of Arctic Flora and Fauna, Akureyri (2013).Desjardins, É. et al. Survey of the vascular plants of Alert (Ellesmere Island, Canada), a polar desert at the northern tip of the Americas. CheckList 17, 181–225 (2021).
Google Scholar 
Keith, L. B., Meslow, E. C. & Rongstad, O. J. Techniques for snowshoe hare population studies. J. Wildl. Manag. 32, 801–812 (1968).
Google Scholar 
Davidson, S. C. et al. Ecological insights from three decades of animal movement tracking across a changing Arctic. Science 370, 712–715 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Wikelski, M., Davidson, S. C. & Kays, R. Movebank: Archive, analysis and sharing of animal movement data. Hosted by the Max Planck Institute of Animal Behavior. http://www.movebank.org (2021).Berteaux, D. Data from: Study ‘Arctic hare Alert—Argos tracking’. MoveBank Data Repository https://doi.org/10.5441/001/1.d5d912c4 (2021).Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2021).Christin, S., St-Laurent, M.-H. & Berteaux, D. Evaluation of Argos telemetry accuracy in the High-Arctic and implications for the estimation of home-range size. PLoS One 10, e0141999 (2015).PubMed 
PubMed Central 

Google Scholar 
QGIS Association. QGIS Geographic Information System (2021).Harris, S. et al. Home-range analysis using radio-tracking data? A review of problems and techniques particularly as applied to the study of mammals. Mamm. Rev. 20, 97–123 (1990).
Google Scholar 
Le Corre, M., Dussault, C. & Côté, S. D. Detecting changes in the annual movements of terrestrial migratory species: Using the first-passage time to document the spring migration of caribou. Mov. Ecol. 2, 19 (2014).PubMed 
PubMed Central 

Google Scholar 
Nicholson, K. L., Arthur, S. M., Horne, J. S., Garton, E. O. & Vecchio, P. A. D. Modeling caribou movements: Seasonal ranges and migration routes of the central Arctic herd. PLoS One 11, e0150333 (2016).PubMed 
PubMed Central 

Google Scholar 
Nelson, M. E., Mech, L. D. & Frame, P. F. Tracking of white-tailed deer migration by global positioning system. J. Mammal. 85, 505–510 (2004).
Google Scholar 
Singh, N. J. & Ericsson, G. Changing motivations during migration: Linking movement speed to reproductive status in a migratory large mammal. Biol. Lett. 10, 20140379 (2014).PubMed 
PubMed Central 

Google Scholar 
Jakes, A. F. et al. Classifying the migration behaviors of pronghorn on their northern range. J. Wildl. Manag. 82, 1229–1242 (2018).
Google Scholar 
Bates, D., Maechler, M., Bolker, B. & Walker, S. lme4: Linear mixed-effects models using Eigen and S4. (2015).Duong, T. ks: Kernel density estimation and kernel discriminant analysis for multivariate data in R. J. Stat. Softw. 21, 1–16 (2007).
Google Scholar 
Gitzen, R. A., Millspaugh, J. J. & Kernohan, B. J. Bandwidth selection for fixed-kernel analysis of animal utilization distributions. J. Wildl. Manag. 70, 1334–1344 (2006).
Google Scholar 
Austin, R. E. et al. Patterns of at-sea behaviour at a hybrid zone between two threatened seabirds. Sci. Rep. 9, 14720 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Gillis, E. A. & Krebs, C. J. Natal dispersal of snowshoe hares during a cyclic population increase. J. Mammal. 80, 933–939 (1999).
Google Scholar 
Dahl, F. & Willebrand, T. Natal dispersal, adult home ranges and site fidelity of mountain hares (Lepus timidus) in the boreal forest of Sweden. Wildl. Biol. 11, 309–317 (2005).
Google Scholar 
Angerbjörn, A. & Flux, J. E. C. Lepus timidus. Mamm. Species 495, 1–11 (1995).
Google Scholar 
Smith, G. W., Stoddart, L. C. & Knowlton, F. F. Long-distance movements of black-tailed jackrabbits. J. Wildl. Manag. 66, 463 (2002).
Google Scholar 
Cote, J. et al. Behavioural synchronization of large-scale animal movements—Disperse alone, but migrate together?. Biol. Rev. 92, 1275–1296 (2017).PubMed 

Google Scholar 
Bauer, S., McNamara, J. M. & Barta, Z. Environmental variability, reliability of information and the timing of migration. Proc. R. Soc. B 287, 20200622 (2020).PubMed 
PubMed Central 

Google Scholar 
Couzin, I. D. Collective animal migration. Curr. Biol. 28, R976–R980 (2018).CAS 
PubMed 

Google Scholar 
Lai, S. et al. Unsuspected mobility of Arctic hares revealed by longest journey ever recorded in a lagomorph. Ecology 103(3), e3620 https://doi.org/10.1002/ecy.3620 (2022).PubMed 

Google Scholar 
Abrahms, B. et al. Suite of simple metrics reveals common movement syndromes across vertebrate taxa. Mov. Ecol. 5, 12 (2017).PubMed 
PubMed Central 

Google Scholar 
Chapman, B. B., Brönmark, C., Nilsson, J. -Å. & Hansson, L.-A. The ecology and evolution of partial migration. Oikos 120, 1764–1775 (2011).
Google Scholar 
Singh, N. J., Börger, L., Dettki, H., Bunnefeld, N. & Ericsson, G. From migration to nomadism: Movement variability in a northern ungulate across its latitudinal range. Ecol. Appl. 22, 2007–2020 (2012).PubMed 

Google Scholar 
Bastille-Rousseau, G. et al. Flexible characterization of animal movement pattern using net squared displacement and a latent state model. Mov. Ecol. 4, 15 (2016).PubMed 
PubMed Central 

Google Scholar 
Mueller, T. & Fagan, W. F. Search and navigation in dynamic environments—From individual behaviors to population distributions. Oikos 117, 654–664 (2008).
Google Scholar 
Krebs, C. J., Boonstra, R. & Boutin, S. Using experimentation to understand the 10-year snowshoe hare cycle in the boreal forest of North America. J. Anim. Ecol. 87, 87–100 (2018).PubMed 

Google Scholar 
Reid, N. & Harrison, A. Post-release GPS tracking of hand-reared Irish hare Lepus timidus hibernicus leverets, Slemish, Co. Antrim, Northern Ireland. J. Wildl. Rehabil. 31, 25 (2011).
Google Scholar 
Weterings, M. J. A. et al. Strong reactive movement response of the medium-sized European hare to elevated predation risk in short vegetation. Anim. Behav. 115, 107–114 (2016).
Google Scholar 
Krebs, C. J., Boutin, S. & Boonstra, R. Ecosystem Dynamics of the Boreal Forest: The Kluane Project (Oxford University Press, 2001).
Google Scholar 
Feierabend, D. & Kielland, K. Movements, activity patterns, and habitat use of snowshoe hares (Lepus americanus) in interior Alaska. J. Mammal. 95, 525–533 (2014).
Google Scholar 
Levänen, R., Pohjoismäki, J. L. O. & Kunnasranta, M. Home ranges of semi-urban brown hares (Lepus europaeus) and mountain hares (Lepus timidus) at northern latitudes. Ann. Zool. Fenn. 56, 107–120 (2019).
Google Scholar 
Nathan, R. et al. A movement ecology paradigm for unifying organismal movement research. Proc. Natl. Acad. Sci. U.S.A. 105, 19052–19059 (2008).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Abrahms, B. et al. Emerging perspectives on resource tracking and animal movement ecology. Trends Ecol. Evol. 36, 308–320 (2021).PubMed 

Google Scholar 
France, R. L. The Lake Hazen trough: A late winter oasis in a polar desert. Biol. Conserv. 63, 149–151 (1993).
Google Scholar 
Jenkins, D. A., Campbell, M., Hope, G., Goorts, J. & McLoughlin, P. Recent trends in abundance of Peary caribou (Rangifer tarandus pearyi) and muskoxen (Ovibos moschatus) in the Canadian Arctic Archipelago, Nunavut 233.Mech, L. Proportion of calves and adult muskoxen, Ovibos moschatus killed by gray wolves, Canis lupus, in July on Ellesmere Island (USGS Northern Prairie Wildlife Research Center, 2010).
Google Scholar 
Gunn, A., Miller, F., Barry, S. & Buchan, A. A near-total decline in caribou on Prince of Wales, Somerset, and Russell Islands, Canadian Arctic. Arctic 59, 1–13 (2006).
Google Scholar 
Edwards, J. Diet shifts in moose due to predator avoidance. Oecologia 60, 185–189 (1983).ADS 
PubMed 

Google Scholar 
Gustine, D. D., Parker, K. L., Lay, R. J., Gillingham, M. P. & Heard, D. C. Calf survival of woodland caribou in a multi-predator ecosystem. Wildl. Monogr. 165, 1–32 (2006).
Google Scholar 
Klein, D. & Bay, C. Diet selection by vertebrate herbivores in the High Arctic of Greenland. Ecography 14, 152–155 (1991).
Google Scholar 
Parks Canada. Resource Description and Analysis—Ellesmere Island National Park Reserve Vol. 1 (Natural Resource Conservation Section, Parks Canada, Department of Canadian Heritage, 1994).
Google Scholar 
Parks Canada. Quttinirpaaq National Park of Canada: Management plan 76. https://www.pc.gc.ca/en/pn-np/nu/quttinirpaaq/info/index/gestion-management-2009 (2009).Winkler, D. W. et al. Cues, strategies, and outcomes: How migrating vertebrates track environmental change. Mov. Ecol. 2, 10 (2014).
Google Scholar 
Robinson, R. et al. Travelling through a warming world: Climate change and migratory species. Endang. Species Res. 7, 87–99 (2009).ADS 

Google Scholar  More

Site descriptionIn 2012, a field experiment was conducted in the Aksu National Experimental Station of Oasis Farmland Ecosystem27 (40°37′ N, 80°45′ E, altitude 1028 m) (Fig. 1), located in the west of Tarim River Basin in Xinjiang Province, China. The experimental area had a typical temperate arid climate. During the study period (May to October), the average minimum and maximum temperatures varied between 16.7 and 34.8 ℃ respectively.Figure 1Location of the Aksu National Experimental Station of Oasis Farmland Ecosystem (the map was created by software: QGIS Version 3.16.15 LTR: URL, https://www.qgis.org/en/site/).Full size imageThe cotton fields where the experiment conducted were public land, belong to Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, China. With the permissions of Xinjiang Institute of Ecology and Geography, we conducted experiments in the cotton field of the Aksu National Experimental Station of Oasis Farmland Ecosystem.Experimental designTwo treatments, each 10 m × 10 m in size, were established on one of cotton fields at the Aksu National Experimental Station of Oasis Farmland Ecosystem on April 5, 2012.One treatment planting cotton with TC method, the other with MC method. For the MC method, a high-density and air-tight transparent polythene film (0.01–0.02 mm thick, 1.25 m wide) was placed over the soil surface before sowing. Small holes (0.02 m × 0.02 m, at 0.1 m intervals within a row) in the plastic film were made to place cotton seeds. Four rows were sown on each strip of plastic film. For the TC treatment, the plants were sown as that for the MC treatment. The planting density (266 667plant ha−1) and irrigation pattern (frequency and volume of irrigation) for the TC method were entirely consistent with those for the MC method.Half-hourly measurements of soil CO2 efflux, soil temperature and moisture were made on 6 June 2012. The whole experiment was completed on 4 November 2012. According to irrigation, the whole experiment can be divided into three stages: stage before irrigation (from 6 to 24 June), during irrigation (from 25 June to 10 October) and irrigation stop stage (from 11October to 4 November). During the irrigation period, we conducted seven times of irrigation (once in two week). The water-soluble compound fertilizer (N + P2O5 + K2O ≥ 51%) was used for fertilization in the experimental field, and the application rate was 30 g m−2. We dissolved water-soluble compound fertilizer in water and sprayed into the field by sprayer. During the irrigation period, the fertilizer was applied for 5 times.The cottonseeds we used in this study comply with the provisions of the regulations of the People’s Republic of China on Seed Administration and the detailed rules for the implementation of crop seeds. The fertilization we used in this study comply with the provisions of the People’s Republic of China on Chemical fertilizer standard. All the experiments we conducted in the cotton field of Aksu oasis farmland ecosystem National Experimental Station met the provisions of the agricultural law of the People’s Republic of China. We also carried out the experiment of this study under the guidance of the provisions of the measures for the administration of national field scientific observation and research stations.Field measurement of soil CO2 concentrationSolid-state CO2 sensors (GMM221 and GMM222, Vaisala, Finland) were installed in the midpoint of each treatment to measure soil CO2 concentration. A cable connected each soil probe with a transimitter body placed on the ground. The transimitter sent output signals from the probe to a data logger (CR1000, Campbell Scientific Inc., Logan, UT, USA) and to an optional LCC display on the transmitter.In each treatment, four CO2 concentration sensors were buried at depths of 0 cm, 5 cm, 10 cm and 15 cm. Soil CO2 concentrations were recorded once in 30 min. The measurement of soil CO2 concentrations were conducted from 6 June 2012 to 4 November 2012.On 8 November, these sensors were excavated and recalibrated in the laboratory. We found no change in the slope or offset.Environmental and soil CO2 efflux measurementsThe soil water content and temperature at the same soil depth with solid-state CO2 sensors were measured on the cotton fields at the Aksu National Experimental Station of Oasis Farmland Ecosystem27,28, respectively. Soil volumetric water content and soil temperature were measured using soil moisture probes (pF-Meter, EcoTech GmbH, Bonn, Germany)26 and temperature probes (PT100,Heraeus Sensor Technology, Kleinostheim, Germany)26, respectively.Bulk density was determined by core method29. Briefly, a cylindrical metal sampler (volume of 100cm3) was inserted into the soil and carefully removed to preserve the sample. The sample was oven-dried at 105 °C and weighed. The ratio between dry weight of the soil sample and the cylinder volume was applied to provide the bulk density.Half-hourly soil CO2 efflux measurements were conducted using a closed dynamic chamber method26 (CIRAS-1 PP Systems, Hitchin, UK) on the TC treatment, beginning on 6 June 2012. A chamber, with a diameter of 9.96 cm and a volume of 1, 170 cm3 was inserted into the soil at depth of 3 cm. Soil CO2 concentrations were measured by infrared gas analyzer. The collecting of CO2 from each sampling point took 120 s to get reliable estimates of soil CO2 efflux.Data analysisIn order to calculate CO2 efflux in soil, Fick’s first law of diffusion was used:$$F_{i} = – D_{s} frac{dc}{{dz}}$$
(1)
where Fi is the CO2 efflux at depth zi, Ds the CO2 diffusion coefficient in the soil, and dc / dz the vertical soil CO2 gradient. In this study, the vertical CO2 gradient (dC/dz) was approximately a constant at different depths of soil in our site for the field conditions experienced in the TC treatment during study period. However, a quadratic function of depth to concentrations fitted to soil CO2 concentration gradients in the MC treatment.Ds can be estimated as$$D_{s} = xi D_{a}$$
(2)
where ξ is the gas tortuosity factor and Da is the CO2 diffusion coefficient in free air. The effect of temperature and pressure on Da is given by$$D_{a} = D_{a} 0left( {frac{T}{293.15}} right)^{1.75} left( {frac{P}{101.3}} right)$$
(3)
where T is the temperature (K), P the air pressure (kPa), Dao a reference value of Da at 20 °C (293.15 K) and 101.3 kPa, and is given as 14.7 mm2 s–130 .There are several empirical models in the literature for computing ξ31. We used the Millington–Quirk model32:$$xi = frac{{alpha^{10/3} }}{{phi^{2} }}$$
(4)
where a is the volumetric air content (air-filled porosity), Φ is the porosity. Note,$$phi = alpha + theta = 1 – frac{{rho_{b} }}{{rho_{m} }}$$
(5)
where ρb is the bulk density, and ρm is the particle density for the mineral soil.Soil surface CO2 efflux was calculated using the CO2 gradient flux method based on CO2 concentrations within the soil profile1. Briefly, the flux of CO2 between any two layers in the soil profile was calculated using the Moldrup model33.In order to determine soil CO2 storage, the equation for CO2 was performed.$${S}_{C{O}_{2}}=frac{partial (aC)}{partial t}$$
(6)
where C (ppm) is the concentration of CO2 within the soil pores, (a) is the aerial porosity of the soil layer, D is the molecular diffusivity of CO2 with the soil, and S(µmol m−3 s−1)is the source strength in the soil layer at depth.We determined temperature responses for soil CO2 efflux using the van’t Hoff equation34 (Eq. 7);$$R = R0e^{BT}$$
(7)
where R is soil CO2 efflux, T is soil temperature (°C) at 10 cm depth, and R0 is the soil respiration rate at a reference temperature of 0 °C (µmol m−3 s−1).The Q10 value for Eq. (8) was calculated according to definition as:$$Q_{{{1}0}} = R_{{{text{T}} + {1}0}} /R_{{text{T}}} = {text{ e}}^{{{1}0{text{B}}}}$$
(8)
where RT and RT+10 are Rr or Rd rates at temperature T and T + 10, respectively. The Q10 value is independent of temperature in Eq. (8). More

Global patterns of ARG distributionWe used a set of 4572 metagenomic samples to illustrate the global patterns of ARG distribution (Supplementary Data 1). These samples were collected from six types of habitats: air, aquatic, terrestrial, engineered, humans and other hosts (Fig. 1a and Supplementary Data 1). From these samples, we identified a total of 2561 ARGs that conferred resistance to 24 drug classes of antibiotics based on the Comprehensive Antibiotic Research Database (CARD). Of these, 2401 were genes conferring resistance to only one drug class, and 160 conferred resistances to multiple drug classes (Supplementary Data 2). Twenty-five ARGs were found in more than 75% samples, however, the frequency of most ARGs (2313/2561) were More