Ecology

Plankton are microscopic organisms at the base of aquatic food webs and therefore essential to all life on Earth. In our view, international legal protection of plankton is urgently needed because of their high susceptibility to the effects of climate change, including ocean warming and acidification.
Competing Interests
The authors declare no competing interests. More

Commensal rodents serve as important reservoirs of rodent-borne pathogens. Efforts to reduce the risk of pathogen transmission include decimating rodent populations, altering access pathways, upholding good waste management practices and denying easy access to food sources. In our study, we examined the incidence of rodent activity in waste collection premises in public residential estates in Singapore and examined the factors associated with rodent activity to inform the priority of rodent control measures of resource limited municipal estate managers.Of the three types of waste collection premises, rodent activity had the highest incidence in refuse bin centres followed by CRCs and IRC bin chambers. Refuse bin centres are prone to refuse spillage because refuse is manually transferred from refuse collection carts into bulk bins and refuse compactors located within the centres. Bin centres tend to be larger than CRCs and IRC bin chambers and the storage of bulky waste that provide additional areas of rodent harbourage are a common sight in Singapore. IRC chambers and refuse bin centres in combination far outnumber CRCs, and the former two are a distinct characteristic of older public housing estates in Singapore. This suggests that older public housing estates have a higher propensity for rodent infestation compared to newer ones. Aging infrastructure can also provide a greater number of harbourage areas and alternate access pathways for rodent travel that increase their ability to obtain food sources. Our study findings were in support of previous studies which found that older infrastructure was associated with a greater likelihood of rodent activity22,23.We also found that the number of IRC bin chutes was positively associated with rodent infestation. Fluids from food waste in IRC bin chambers are drained directly into a sanitary line that is common to all other bin chambers within the same building. A possible explanation therefore is that rodents which find their way into the sanitary line can easily access all bin chambers in the same building. This suggests that preventing individual bin chamber access may reduce food availability to rodents which traverse the sanitary line in search of food sources.In the present study, we observed that rodent sightings were relatively higher in some months in the first half of the calendar year compared to the second half. Even though our estimates were positive, those for some months were not statistically significant. In Singapore, end-December, January to February are usually associated with increased food production due to the year-end (Christmas and New Year celebrations) and early-year (Chinese New Year) festivities. A proportionate increase in food waste over that period could improve survivability of rodents that leads to increased mating and reproduction. We therefore postulate that the higher seasonal rodent activity is plausible but recommend that future studies be conducted with sufficient longitude to examine the differences in the seasonal pattern across the three categories of premises more closely. A previous study in Harbin, China27 reported a seasonal pattern in the age composition of R. norvegicus while an ecological study on R. norvegicus in Salvador Brazil did not find any difference in the number of rats trapped between the dry and rainy seasons28. The inconsistent seasonal findings between studies could be due to the differences in the climate, degree of urbanization and environmental conditions of study locations.The relative rise in rodent activity in the first half of the year coupled with older estates being at greater risk of rodent activity suggest that municipal town councils which prioritize regular infrastructural repairs and improvements in older estates and complete them in the second half of each calendar year would help mitigate the anticipated rise in the first half of the new calendar year.In our study, we examined the relations between visual cues and rodent activity to help estate managers prioritise their control efforts. We found that rodent droppings were a common positive predictor of rodent activity across all three categories of waste collection premises. In particular, the odds of droppings in IRC bin chambers were the highest among the three categories of premises. We hypothesize that the probability of rodent dropping sightings was in part related to the accessibility of food waste and thus time spent by rodents within the respective waste collection premises. Each IRC chamber contains an open top bin that receives waste that is disposed down the IRC chamber chute. Food waste in IRC bin chambers are thus more easily accessed by rodents compared to in CRCs where waste is stored in a compactor and in bin centres where bulk bins are covered until the waste is compacted or collected.In Salvador, Brazil, the presence of Rattus norvegicus droppings were independently associated with an increased risk of Leptospira infection in humans29. Further research on site-specific Leptospira infection risks in Singapore are required to affirm the utility of droppings as an indicator for Leptospira infection risk. In addition, rub marks and gnaw marks were also positive predictors of rodent activity in CRCs and IRC bin chambers. A study in Chile reported that gnaw marks and holes, as well as grease or rub marks left behind by rodent travel were indicators of rodent activity30. A previous study carried out in an urban city in Taiwan reported that rodent droppings and rub marks were well correlated with rodent infestation31. Our findings, which were in support of these previous studies, suggest that estate managers can maximise the cost effectiveness of their resources by focusing their control efforts based on visual cues without relying solely on trapping activities for surveillance.We found a positive relationship between the number of rodent burrows and rodent activity in all three waste collection premises, though this was only significant for refuse bin centres. That the direction of effect for burrows was consistent these three premises, was a reassuring observation. It is possible that we did not have enough study power to establish the observed positive relations in CRCs and IRC bin chambers. Therefore future studies should seek to confirm our findings. R. norvegicus excavate extensive burrow systems that are able to house a large number of rats32. They exhibit a strong preference for creating burrows in loose soil and on sloping terrain33 and construct shallow burrows in close proximity to water bodies and food sources34. As rodent burrows are primarily used for nesting, food storage and harbourage purposes35, burrows can provide important information about the extent of rodent activity in an area and may be used as an indicator for estate managers to focus their investigations.A previous study in New York, United States found that the presence of numerous restaurants, or having older infrastructure were associated with increased levels of R. norvegicus22. Unexpectedly, we did not find any evidence that the number of dining establishments was associated with rodent activity. However, instances of rodent activity have been reported in food establishments in Singapore36,37,38. We hypothesize that rodent movement is restricted to the surrounding area of the food establishments due to the plethora of food available, with little reason for rodents to venture into waste collection premises. Future studies examining the relationship between rodent activity in food establishments and waste collection premises are required to confirm this.In our study, the presence of gnaw marks (aOR: 5.61), rub marks (aOR: 5.04) in CRCs and rodent droppings in CRCs (aOR: 6.20), IRC bin chambers (aOR: 90.84) and bin centres (aOR: 3.61) had the largest strengths of association with rodent activity. Comparatively, in a study in Johannesburg, South Africa, predictors such as dampness (aOR: 2.54) and cracks (aOR: 1.92) in homes had relatively smaller effects on rodent activity20, while a study in Salvador, Brazil found relatively larger effects of homes with dilapidated fences and walls (aOR: 8.95) and those built on earthen slopes (aOR: 4.95)21. This suggests that rodent activity can be strongly influenced by site- and setting-specific factors, and supports the body of evidence on the strong adaptability of rodents in our urban environment”.Urban environments have the capacity to alter the biology of the pathogens, hosts and vectors, which can influence disease transmission39. The proximate setting of dense urban environments allows for close contact between humans and synanthropic rodents, thereby increasing the transmission risk of zoonotic diseases4. In addition to causing diseases in human populations, urban rats are also known to compromise food safety, damage infrastructure and cause mental health distress25,40. The responsibility of rodent control in residential estates is important but may be one among many other competing public health and estate management responsibilities that municipal town councils have to undertake. Consequently, estate managers have to prioritize their limited resources in order to maximise the cost effectiveness of their resource allocation choices. Based on our study findings, we recommend that estate managers adopt a risk-based approach in vector control resource allocation in waste collection premises according to infrastructural age and visual cues for rodent activity.IRC bin chambers which are a distinct feature of the oldest residential buildings, were observed with a substantially higher odds of rodent activity compared to the other categories of waste collection premises. This suggests that rodent control resource allocation should be prioritized in older residential estates. The clear seasonal pattern of rodent activity in CRCs suggests that estate managers can increase their rodent control activities thereat in the first half of the year. Finally, easy access to food waste directly increases the probability of survival and consequently the rodent population size. Future research should examine the quality of municipal solid waste management and the waste processing flow in residential estates to determine how rodent access to food waste can be further minimized to reduce the population of rodents.Study strengths and limitationsWe analysed data from all public residential estates in Singapore; our findings are thus generalizable at the national level. The use of outcomes and independent measures from individual waste collection premises over multiple cycles of inspection provided stronger evidence for causal inferences. We analysed data over 12 months to account for within-year variations that could influence the outcome measure. Rodents were visually identified without molecular speciation because no trapping was carried out. Though the majority of rodents were observed to be Rattus norvegicus, which is the most common species of rodents in public housing estates in Singapore, we could not rule out misclassification of rodents. However, our findings remain relevant for municipal authorities seeking to prioritize resources for vector control in waste collection premises under their care. More

Falkowski, P. The power of plankton. Nature 483, 17–20 (2012).Article 

Google Scholar 
Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 281, 237 (1998).Article 
CAS 
PubMed 

Google Scholar 
Grossart, H. P. & Simon, M. Significance of limnetic organic aggregates (lake snow) for the sinking flux of particulate organic matter in a large lake. Aquat. Microb. Ecol. 15, 115–125 (1998).Article 

Google Scholar 
Weyhenmeyer, G. A. & Bloesch, J. The pattern of particle flux variability in Swedish and Swiss lakes. Sci. Total Environ. 266, 69–78 (2001).Article 
CAS 
PubMed 

Google Scholar 
Fender, C. K. et al. Investigating particle size-flux relationships and the biological pump across a range of plankton ecosystem states from coastal to oligotrophic. Front. Marine Sci. 6, https://doi.org/10.3389/fmars.2019.00603 (2019).Iversen, M. H., Nowald, N., Ploug, H., Jackson, G. A. & Fischer, G. High resolution profiles of vertical particulate organic matter export off Cape Blanc, Mauritania: Degradation processes and ballasting effects. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 57, 771–784 (2010).Article 
CAS 

Google Scholar 
Griffiths, J. R. et al. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Change Biol. 23, 2179–2196 (2017).Article 

Google Scholar 
Diaz, R. J. & Rosenberg, R. Spreading dead zones and consequences for marine ecosystems. Science 321, 926–929 (2008).Article 
CAS 
PubMed 

Google Scholar 
Jenny, J. P. et al. Global spread of hypoxia in freshwater ecosystems during the last three centuries is caused by rising local human pressure. Glob. Chang Biol. 22, 1481–1489 (2016).Article 
PubMed 

Google Scholar 
Carstensen, J., Andersen, J. H., Gustafsson, B. G. & Conley, D. J. Deoxygenation of the Baltic Sea during the last century. Proc. Natl Acad. Sci. 111, 5628–5633 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Simon, M., Grossart, H. P., Schweitzer, B. & Ploug, H. Microbial ecology of organic aggregates in aquatic ecosystems. Aquat. Microb. Ecol. 28, 175–211 (2002).Article 

Google Scholar 
Burd, A. B. & Jackson, G. A. Particle aggregation. Ann. Rev. Mar. Sci. 1, 65–90 (2009).Article 
PubMed 

Google Scholar 
Kiørboe, T., Lundsgaard, C., Olesen, M. & Hansen, J. L. S. Aggregation and sedimentation processes during a spring phytoplankton bloom: A field experiment to test coagulation theory. J. Mar. Res. 52, 297–323 (1994).Article 

Google Scholar 
Boyd, P. W. & Trull, T. W. Understanding the export of biogenic particles in oceanic waters: Is there consensus? Prog. Oceanogr. 72, 276–312 (2007).Article 

Google Scholar 
Legendre, L. & Rivkin, R. B. Fluxes of carbon in the upper ocean: regulation by food-web control nodes. Mar. Ecol. Prog. Ser. 242, 95–109 (2002).Article 

Google Scholar 
Kaneko, H. et al. Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean. iScience 24, 102002 (2021).Article 
CAS 
PubMed 

Google Scholar 
Guidi, L. et al. Plankton networks driving carbon export in the oligotrophic ocean. Nature 532, 465 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Grossart, H.-P. et al. Fungi in aquatic ecosystems. Nat. Rev. Microbiol. 17, 339–354 (2019).Article 
CAS 
PubMed 

Google Scholar 
Amend, A. et al. Fungi in the marine environment: Open questions and unsolved problems. mBio 10, e01189–01118 (2019).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ortiz-Álvarez, R., Triadó-Margarit, X., Camarero, L., Casamayor, E. O. & Catalan, J. High planktonic diversity in mountain lakes contains similar contributions of autotrophic, heterotrophic and parasitic eukaryotic life forms. Sci. Rep. 8, 4457 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Gutiérrez, M. H., Pantoja, S., Tejos, E. & Quiñones, R. A. The role of fungi in processing marine organic matter in the upwelling ecosystem off Chile. Mar. Biol. 158, 205–219 (2011).Article 

Google Scholar 
Edgcomb, V. P., Beaudoin, D., Gast, R., Biddle, J. F. & Teske, A. Marine subsurface eukaryotes: The fungal majority. Environ. Microbiol. 13, 172–183 (2011).Article 
CAS 
PubMed 

Google Scholar 
Frenken, T. et al. Integrating chytrid fungal parasites into plankton ecology: research gaps and needs. Environ. Microbiol. 19, 3802–3822 (2017).Article 
PubMed 

Google Scholar 
Van den Wyngaert, S. et al. Seasonality of parasitic and saprotrophic zoosporic fungi: linking sequence data to ecological traits. ISME J. 16, 2242–2254 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Gsell, A. S. et al. Long-term trends and seasonal variation in host density, temperature, and nutrients differentially affect chytrid fungi parasitising lake phytoplankton. Freshwat. Biol. https://doi.org/10.1111/fwb.13958 (2022).Gutiérrez, M. H., Jara, A. M. & Pantoja, S. Fungal parasites infect marine diatoms in the upwelling ecosystem of the Humboldt current system off central Chile. Environ. Microbiol. 18, 1646–1653 (2016).Article 
PubMed 

Google Scholar 
Kilias, E. S. et al. Chytrid fungi distribution and co-occurrence with diatoms correlate with sea ice melt in the Arctic Ocean. Commun. Biol. 3, 183 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hassett, B. T., Ducluzeau, A. L. L., Collins, R. E. & Gradinger, R. Spatial distribution of aquatic marine fungi across the western Arctic and sub-arctic. Environ. Microbiol. 19, 475–484 (2017).Article 
CAS 
PubMed 

Google Scholar 
Lepelletier, F. et al. Dinomyces arenysensis gen. et sp. nov. (Rhizophydiales, Dinomycetaceae fam. nov.), a chytrid infecting marine dinoflagellates. Protist 165, 230–244 (2014).Article 
PubMed 

Google Scholar 
Hassett, B. T. & Gradinger, R. Chytrids dominate arctic marine fungal communities. Environ. Microbiol. 18, 2001–2009 (2016).Article 
CAS 
PubMed 

Google Scholar 
Le Calvez, T., Burgaud, G., Mahé, S., Barbier, G. & Vandenkoornhuyse, P. Fungal diversity in deep-sea hydrothermal ecosystems. Appl. Environ. Microbiol. 75, 6415–6421 (2009).Article 
PubMed 
PubMed Central 

Google Scholar 
Richards, T. A. et al. Molecular diversity and distribution of marine fungi across 130 european environmental samples. Proc. R. Soc. B Biol. Sci. 282, 20152243 (2015).Article 

Google Scholar 
Taylor, J. D. & Cunliffe, M. Multi-year assessment of coastal planktonic fungi reveals environmental drivers of diversity and abundance. ISME J. 10, 2118–2128 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Comeau, A. M., Vincent, W. F., Bernier, L. & Lovejoy, C. Novel chytrid lineages dominate fungal sequences in diverse marine and freshwater habitats. Sci. Rep. 6, 30120 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wang, Y., Sen, B., He, Y., Xie, N. & Wang, G. Spatiotemporal distribution and assemblages of planktonic fungi in the coastal waters of the Bohai Sea. Front. Microbiol. 9, 584 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Gao, Z., Johnson, Z. I. & Wang, G. Molecular characterization of the spatial diversity and novel lineages of mycoplankton in Hawaiian coastal waters. ISME J. 4, 111–120 (2009).Article 
PubMed 

Google Scholar 
Duan, Y. et al. A high-resolution time series reveals distinct seasonal patterns of planktonic fungi at a temperate coastal ocean site (Beaufort, North Carolina, USA). Appl. Environ. Microbiol. 84, e00967–00918 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cleary, A. C., Søreide, J. E., Freese, D., Niehoff, B. & Gabrielsen, T. M. Feeding by Calanus glacialis in a high arctic fjord: Potential seasonal importance of alternative prey. ICES J. Mar. Sci. 74, 1937–1946 (2017).Article 

Google Scholar 
Renaud, P. E., Morata, N., Carroll, M. L., Denisenko, S. G. & Reigstad, M. Pelagic–benthic coupling in the western Barents Sea: Processes and time scales. Deep Sea Res. Part II: Topical Stud. Oceanogr. 55, 2372–2380 (2008).Article 
CAS 

Google Scholar 
Lepère, C., Ostrowski, M., Hartmann, M., Zubkov, M. V. & Scanlan, D. J. In situ associations between marine photosynthetic picoeukaryotes and potential parasites – a role for fungi? Environ. Microbiol. Rep. 8, 445–451 (2016).Article 
PubMed 

Google Scholar 
Kagami, M., Gurung, T. B., Yoshida, T. & Urabe, J. To sink or to be lysed? Contrasting fate of two large phytoplankton species in Lake Biwa. Limnol. Oceanogr. 51, 2775–2786 (2006).Article 

Google Scholar 
Gerphagnon, M., Colombet, J., Latour, D. & Sime-Ngando, T. Spatial and temporal changes of parasitic chytrids of cyanobacteria. Sci. Rep. 7, 6056 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Ibelings, B. W. et al. Chytrid infections and diatom spring blooms: Paradoxical effects of climate warming on fungal epidemics in lakes. Freshwat. Biol. 56, 754–766 (2011).Article 

Google Scholar 
Gsell, A. S. et al. Spatiotemporal variation in the distribution of chytrid parasites in diatom host populations. Freshwat. Biol. 58, 523–537 (2013).Article 

Google Scholar 
Grami, B. et al. Functional effects of parasites on food web properties during the spring diatom bloom in Lake Pavin: A linear inverse modeling analysis. PLOS ONE. 6, e23273 (2011).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Klawonn, I. et al. Characterizing the “fungal shunt”: Parasitic fungi on diatoms affect carbon flow and bacterial communities in aquatic microbial food webs. Proc. Natl Acad. Sci. 118, e2102225118 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kagami, M., Miki, T. & Takimoto, G. Mycoloop: Chytrids in aquatic food webs. Front. Microbiol. 5, 166 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Laundon, D. & Cunliffe, M. A call for a better understanding of aquatic chytrid biology. Front. Fungal Biol. 2, https://doi.org/10.3389/ffunb.2021.708813 (2021).Ploug, H., Iversen, M. H. & Fischer, G. Ballast, sinking velocity, and apparent diffusivity within marine snow and zooplankton fecal pellets: Implications for substrate turnover by attached bacteria. Limnol. Oceanogr. 53, 1878–1886 (2008).Article 

Google Scholar 
Laurenceau-Cornec, E. C., Trull, T. W., Davies, D. M., De La Rocha, C. L. & Blain, S. Phytoplankton morphology controls on marine snow sinking velocity. Mar. Ecol. Prog. Ser. 520, 35–56 (2015).Article 

Google Scholar 
Tréguer, P. et al. Influence of diatom diversity on the ocean biological carbon pump. Nat. Geosci. 11, 27–37 (2018).Article 

Google Scholar 
Alldredge, A. L., Gotschalk, C., Passow, U. & Riebesell, U. Mass aggregation of diatom blooms: Insights from a mesocosm study. Deep Sea Res. Part II: Topical Stud. Oceanogr. 42, 9–27 (1995).Article 
CAS 

Google Scholar 
Seto, K., Van den Wyngaert, S., Degawa, Y. & Kagami, M. Taxonomic revision of the genus Zygorhizidium: Zygorhizidiales and Zygophlyctidales ord. nov. (Chytridiomycetes, Chytridiomycota). Fungal Syst. Evol. 5, 17–38 (2020).Article 
CAS 
PubMed 

Google Scholar 
Engel, A. in Practical Guidelines for the Analysis of Seawater (eds Wurl O & Raton B) (CRC Press, 2009).Cisternas-Novoa, C., Lee, C. & Engel, A. A semi-quantitative spectrophotometric, dye-binding assay for determination of Coomassie Blue stainable particles. Limnol. Oceanogr. Methods. 12, 604–616 (2014).Article 

Google Scholar 
Passow, U. & Alldredge, A. L. A dye-binding assay for the spectrophotometric measurement of transparent exopolymer particles (TEP). Limnol. Oceanogr. 40, 1326–1335 (1995).Article 
CAS 

Google Scholar 
Iversen, M. H. & Ploug, H. Temperature effects on carbon-specific respiration rate and sinking velocity of diatom aggregates – potential implications for deep ocean export processes. Biogeosciences 10, 4073–4085 (2013).Article 

Google Scholar 
van der Jagt, H., Friese, C., Stuut, J.-B. W., Fischer, G. & Iversen, M. H. The ballasting effect of Saharan dust deposition on aggregate dynamics and carbon export: Aggregation, settling, and scavenging potential of marine snow. Limnol. Oceanogr. 63, 1386–1394 (2018).Article 

Google Scholar 
Grossart, H. P. & Ploug, H. Microbial degradation of organic carbon and nitrogen on diatom aggregates. Limnol. Oceanogr. 46, 267–277 (2001).Article 
CAS 

Google Scholar 
Iversen, M. H. & Ploug, H. Ballast minerals and the sinking carbon flux in the ocean: Carbon-specific respiration rates and sinking velocity of marine snow aggregates. Biogeosciences 7, 2613–2624 (2010).Article 
CAS 

Google Scholar 
Ploug, H. & Grossart, H. P. Bacterial growth and grazing on diatom aggregates: Respiratory carbon turnover as a function of aggregate size and sinking velocity. Limnol. Oceanogr. 45, 1467–1475 (2000).Article 
CAS 

Google Scholar 
Belcher, A. et al. Depth-resolved particle-associated microbial respiration in the northeast Atlantic. Biogeosciences 13, 4927–4943 (2016).Article 

Google Scholar 
Ploug, H., Grossart, H. P., Azam, F. & Jørgensen, B. B. Photosynthesis, respiration, and carbon turnover in sinking marine snow from surface waters of Southern California Bight: Implications for the carbon cycle in the ocean. Mar. Ecol. Prog. Ser. 179, 1–11 (1999).Article 
CAS 

Google Scholar 
Turner, J. T. Zooplankton fecal pellets, marine snow, phytodetritus and the ocean’s biological pump. Prog. Oceanogr. 130, 205–248 (2015).Article 

Google Scholar 
Nguyen, T. T. H. et al. Microbes contribute to setting the ocean carbon flux by altering the fate of sinking particulates. Nat. Commun. 13, 1657 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Zakem, E. J., Cael, B. B. & Levine, N. M. A unified theory for organic matter accumulation. Proc. Natl Acad. Sci. 118, e2016896118 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Alcolombri, U. et al. Sinking enhances the degradation of organic particles by marine bacteria. Nat. Geosci. 14, 775–780 (2021).Article 
CAS 

Google Scholar 
Buesseler, K. O. & Boyd, P. W. Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean. Limnol. Oceanogr. 54, 1210–1232 (2009).Article 
CAS 

Google Scholar 
Henson, S., Le Moigne, F. & Giering, S. Drivers of carbon export efficiency in the global ocean. Glob. Biogeochem. Cycles. 33, 891–903 (2019).Article 
CAS 

Google Scholar 
Gsell, A. S., De Senerpont Domis, L. N., Verhoeven, K. J. F., Van Donk, E. & Ibelings, B. W. Chytrid epidemics may increase genetic diversity of a diatom spring-bloom. ISME J. 7, 2057–2059 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Agha, R., Saebelfeld, M., Manthey, C., Rohrlack, T. & Wolinska, J. Chytrid parasitism facilitates trophic transfer between bloom-forming cyanobacteria and zooplankton (Daphnia). Sci. Rep. 6, 35039 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rasconi, S. et al. Parasitic chytrids upgrade and convey primary produced carbon during inedible algae proliferation. Protist 171, 125768 (2020).Article 
CAS 
PubMed 

Google Scholar 
Guidi, L. et al. Effects of phytoplankton community on production, size, and export of large aggregates: A world-ocean analysis. Limnol. Oceanogr. 54, 1951–1963 (2009).Article 

Google Scholar 
Boyd, P. W. & Newton, P. P. Does planktonic community structure determine downward particulate organic carbon flux in different oceanic provinces?. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 46, 63–91 (1999).Article 
CAS 

Google Scholar 
van der Jagt, H., Wiedmann, I., Hildebrandt, N., Niehoff, B. & Iversen, M. H. Aggregate feeding by the copepods Calanus and Pseudocalanus controls carbon flux attenuation in the arctic shelf sea during the productive period. Front. Mar. Sci. 7, 543124 (2020).Article 

Google Scholar 
Steinberg, D. K. et al. Bacterial vs. zooplankton control of sinking particle flux in the ocean’s twilight zone. Limnol. Oceanogr. 53, 1327–1338 (2008).Article 

Google Scholar 
Cavan, E. L., Henson, S. A., Belcher, A. & Sanders, R. Role of zooplankton in determining the efficiency of the biological carbon pump. Biogeosciences 14, 177–186 (2017).Article 
CAS 

Google Scholar 
Gachon, C. M. M., Küpper, H., Küpper, F. C. & Šetlík, I. Single-cell chlorophyll fluorescence kinetic microscopy of Pylaiella littoralis (Phaeophyceae) infected by Chytridium polysiphoniae (Chytridiomycota). Eur. J. Phycol. 41, 395–403 (2006).Article 

Google Scholar 
Senga, Y., Yabe, S., Nakamura, T. & Kagami, M. Influence of parasitic chytrids on the quantity and quality of algal dissolved organic matter (AOM). Water Res. 145, 346––353 (2018).Article 
PubMed 

Google Scholar 
Roberts, C., Allen, R., Bird, K. E. & Cunliffe, M. Chytrid fungi shape bacterial communities on model particulate organic matter. Biol. Lett. 16, 20200368 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Blackburn, N., Fenchel, T. & Mitchell, J. Microscale nutrient patches in planktonic habitats shown by chemotactic bacteria. Science 282, 2254–2256 (1998).Article 
CAS 
PubMed 

Google Scholar 
Smriga, S., Fernandez, V. I., Mitchell, J. G. & Stocker, R. Chemotaxis toward phytoplankton drives organic matter partitioning among marine bacteria. Proc. Natl Acad. Sci. 113, 1576–1581 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Buchan, A., LeCleir, G. R., Gulvik, C. A. & González, J. M. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698 (2014).Article 
CAS 
PubMed 

Google Scholar 
Shibl, A. A. et al. Diatom modulation of select bacteria through use of two unique secondary metabolites. Proc. Natl Acad. Sci. 117, 27445–27455 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Guidi, L. et al. Relationship between particle size distribution and flux in the mesopelagic zone. Deep-Sea Res. Part I Oceanogr. Res. Papers. 55, 1364–1374 (2008).Article 
CAS 

Google Scholar 
Jackson, G. A. et al. Particle size spectra between 1 μm and 1 cm at Monterey Bay determined using multiple instruments. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 44, 1739–1767 (1997).Article 

Google Scholar 
Frenken, T. et al. Warming accelerates termination of a phytoplankton spring bloom by fungal parasites. Glob. Change Biol. 22, 299–309 (2016).Article 

Google Scholar 
Mari, X., Passow, U., Migon, C., Burd, A. B. & Legendre, L. Transparent exopolymer particles: Effects on carbon cycling in the ocean. Prog. Oceanogr. 151, 13–37 (2017).Article 

Google Scholar 
Passow, U. Transparent exopolymer particles (TEP) in aquatic environments. Prog. Oceanogr. 55, 287–333 (2002).Article 

Google Scholar 
Prieto, L. et al. Scales and processes in the aggregation of diatom blooms: high time resolution and wide size range records in a mesocosm study. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 49, 1233–1253 (2002).Article 

Google Scholar 
Kiørboe, T., Andersen, K. P. & Dam, H. G. Coagulation efficiency and aggregate formation in marine phytoplankton. Mar. Biol. 107, 235–245 (1990).Article 

Google Scholar 
Vidal-Melgosa, S. et al. Diatom fucan polysaccharide precipitates carbon during algal blooms. Nat. Commun. 12, 1150 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gärdes, A., Iversen, M. H., Grossart, H. P., Passow, U. & Ullrich, M. S. Diatom-associated bacteria are required for aggregation of Thalassiosira weissflogii. ISME J. 5, 436–445 (2011).Article 
PubMed 

Google Scholar 
Grossart, H. P. & Simon, M. Interactions of planktonic algae and bacteria: Effects on algal growth and organic matter dynamics. Aquat. Microb. Ecol. 47, 163–176 (2007).Article 

Google Scholar 
Short, S. M. The ecology of viruses that infect eukaryotic algae. Environ. Microbiol. 14, 2253–2271 (2012).Article 
PubMed 

Google Scholar 
Carlström, D. The crystal structure of α-chitin (Poly-N-acetyl-d-glucosamine). J. Biophysical Biochemical Cytol. 3, 669–683 (1957).Article 

Google Scholar 
Miklasz, K. A. & Denny, M. W. Diatom sinkings speeds: Improved predictions and insight from a modified Stokes’ law. Limnol. Oceanogr. 55, 2513–2525 (2010).Article 

Google Scholar 
Bidle, K. D. & Azam, F. Accelerated dissolution of diatom silica by marine bacterial assemblages. Nature 397, 508 (1999).Article 
CAS 

Google Scholar 
Gerphagnon, M. et al. Comparison of sterol and fatty acid profiles of chytrids and their hosts reveals trophic upgrading of nutritionally inadequate phytoplankton by fungal parasites. Environ. Microbiol. 21, 949–958 (2019).Article 
CAS 
PubMed 

Google Scholar 
Kagami, M., Von Elert, E., Ibelings, B. W., De Bruin, A. & Van Donk, E. The parasitic chytrid, Zygorhizidium, facilitates the growth of the cladoceran zooplankter, Daphnia, in cultures of the inedible alga, Asterionella. Proc. R. Soc. B Biol. Sci. 274, 1561–1566 (2007).Article 

Google Scholar 
Carney, L. T. & Lane, T. W. Parasites in algae mass culture. Front. Microbiol. 5, 1–8 (2014).Article 

Google Scholar 
Williams, D. M. Synedra, Ulnaria: definitions and descriptions – a partial resolution. Diatom Res. 26, 149–153 (2011).Article 

Google Scholar 
Arar, E. J. & Collins, G. B. Method 445.0: In vitro determination of chlorophyll and phaeophytin a in marine and freshwater algae by fluorescence. U.S. Environemental Protection Agency, Cinncinnati, Ohio Revision 1.2, 1–22 (1997).Klawonn, I., Dunker, S., Kagami, M., Grossart, H.-P., Van den Wyngaert, S. Intercomparison of two fluorescent dyes to visualize parasitic fungi (Chytridiomycota) on phytoplankton. Microb. Ecol. 85, 9–23 (2023).Article 
CAS 
PubMed 

Google Scholar 
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Alldredge, A. L. & Gotschalk, C. In situ settling behavior of marine snow. Limnol. Oceanogr. 33, 351 (1988).Article 

Google Scholar 
Jackson, G. A. Coagulation in a rotating cylinder. Limnol. Oceanogr. Methods. 13, e10018 (2015).Article 

Google Scholar 
Shanks, A. L. & Edmondson, E. W. Laboratory-made artificial marine snow: a biological model of the real thing. Mar. Biol. 101, 463–470 (1989).Article 

Google Scholar 
Cowen, R. K. & Guigand, C. M. In situ ichthyoplankton imaging system (ISIIS): system design and preliminary results. Limnol. Oceanogr. Methods. 6, 126–132 (2008).Article 

Google Scholar 
Jackson, G. A. & Burd, A. B. Simulating aggregate dynamics in ocean biogeochemical models. Prog. Oceanogr. 133, 55–65 (2015).Article 

Google Scholar 
Petrik, C. M., Jackson, G. A. & Checkley, D. M. Aggregates and their distributions determined from LOPC observations made using an autonomous profiling float. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 74, 64–81 (2013).Article 

Google Scholar 
Johnson, C. P., Li, X. & Logan, B. E. Settling velocities of fractal aggregates. Environ. Sci. Technol. 30, 1911–1918 (1996).Article 
CAS 

Google Scholar 
Laurenceau-Cornec, E. C. et al. New guidelines for the application of Stokes’ models to the sinking velocity of marine aggregates. Limnol. Oceanogr. 65, 1264–1285 (2020).Article 
CAS 

Google Scholar 
Ploug, H. & Grossart, H. P. Bacterial production and respiration in suspended aggregates – A matter of the incubation method. Aquat. Microb. Ecol. 20, 21–29 (1999).Article 

Google Scholar 
Berggren, M., Lapierre, J.-F. & del Giorgio, P. A. Magnitude and regulation of bacterioplankton respiratory quotient across freshwater environmental gradients. ISME J. 6, 984–993 (2012).Article 
CAS 
PubMed 

Google Scholar 
R.CoreTeam. R: A language and environment for statistical computing. Vienna, Austria. Retrieved from https://www.R-project.org/ (2016). More

Fan, M. & Yu, X. Impacts of the grain for green program on the spatial pattern of land uses and ecosystem services in mountainous settlements in southwest china. Glob. Ecol. Conserv. 21, 806. https://doi.org/10.1016/j.gecco.2019.e00806 (2020).Article 

Google Scholar 
Huang, L., Shao, Q. & Liu, J. The spatial and temporal patterns of carbon sequestration by forestation in Jiangxi Province. Acta Ecol. Sin. 35, 2105–2118 (2015).
Google Scholar 
Wang, B., Gao, P., Niu, X. & Sun, J. Policy-driven china’s grain to green program: Implications for ecosystem services. Ecosyst. Serv. 27, 38–47. https://doi.org/10.1016/j.gecco.2019.e00806 (2017).Article 

Google Scholar 
Deng, L., Zhou-ping, S. & Li, R. Effects of the grain-for-green program on soil erosion in China. Int. J. Sediment Res. 27, 120–127. https://doi.org/10.1016/s1001-6279(12)60021-3 (2012).Article 

Google Scholar 
Cragun, D. et al. Qualitative comparative analysis: A hybrid method for identifying factors associated with program effectiveness. J. Mix. Methods Res. 10, 251–272. https://doi.org/10.1177/1558689815572023 (2016).Article 
PubMed 

Google Scholar 
Biancardo, S. A. et al. An innovative framework for integrating cost-benefit analysis (cba) within building information modeling (bim). Socio-Econ. Plan. Sci. 85, 1014795. https://doi.org/10.1016/j.seps.2022.101495 (2022).Article 

Google Scholar 
Miller, C. J., Smith, S. N. & Pugatch, M. Experimental and quasi-experimental designs in implementation research. Psychiatry Res. 283, 112452. https://doi.org/10.1016/j.psychres.2019.06.027 (2020).Article 
PubMed 

Google Scholar 
Park, S.-G. et al. Characteristics of the flow field inside and around a square fish cage considering the circular swimming pattern of a farmed fish school: Laboratory experiments and field observations. Ocean. Eng. 261, 112097. https://doi.org/10.1016/j.oceaneng.2022.112097 (2022).Article 

Google Scholar 
Zhou, F. & Wang, X. The carbon emissions trading scheme and green technology innovation in china: A new structural economics perspective. Econ. Anal. Policy 74, 365–381. https://doi.org/10.1016/j.eap.2022.03.007 (2022).Article 

Google Scholar 
Gharehgozli, O. An empirical comparison between a regression framework and the synthetic control method. Q. Rev. Econ. Financ. 81, 70–81. https://doi.org/10.1016/j.qref.2021.05.002 (2021).Article 

Google Scholar 
Salman, M., Long, X., Wang, G. & Zha, D. Paris climate agreement and global environmental efficiency: New evidence from fuzzy regression discontinuity design. Energy Policy 168, 113128. https://doi.org/10.1016/j.enpol.2022.113128 (2022).Article 

Google Scholar 
Stanford, B., Zavaleta, E. & Millard-Ball, A. Where and why does restoration happen? Ecological and sociopolitical influences on stream restoration in coastal california. Biol. Conserv. 221, 219–227. https://doi.org/10.1016/j.biocon.2018.03.016 (2018).Article 

Google Scholar 
Milchakova, L., Bondareva, N. A. & Alexandrov, V. Core areas in the structure of the regional ecological framework of Sevastopol City. South Russia-Ecol. Dev. 17, 102–114. https://doi.org/10.18470/1992-1098-2022-2-102-114 (2022).Article 

Google Scholar 
Miller, R., Nielsen, E. & Huang, C.-H. Ecosystem service valuation through wildfire risk mitigation: Design, governance, and outcomes of the flagstaff watershed protection project (fwpp). Forests 8, 142. https://doi.org/10.3390/f8050142 (2017).Article 

Google Scholar 
Zhao, H. E. A. Spatiotemporal patterns of vegetation conversion under the grain for green program in southwest China. Conserv. Sci. Pract. 4, e604. https://doi.org/10.1111/csp2.604 (2022).Article 
MathSciNet 

Google Scholar 
Guo, B., Xie, T. & Subrahmanyam, M. The impact of china’s grain for green program on rural economy and precipitation: A case study of Yan river basin in the loess plateau. Sustainability 11, 5336. https://doi.org/10.3390/su11195336 (2019).Article 

Google Scholar 
Zuo, Y., Cheng, J. & Fu, M. Analysis of land use change and the role of policy dimensions in ecologically complex areas: A case study in Chongqing. Land 11, 627. https://doi.org/10.3390/land11050627 (2022).Article 

Google Scholar 
Delang, C. O. The effects of china’s grain for green program on migration and remittance. Econ. Agrar. Recursos. Nat. 18, 117–132. https://doi.org/10.2004/ag.econ.281239 (2019).Article 

Google Scholar 
Treacy, P. E. A. Impacts of china’s grain for green program on migration and household income. Environ. Manage. 62, 489–499. https://doi.org/10.1007/s00267-018-1047-0 (2018).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Han, R., Guo, L., Xu, N. & Wang, D. The effect of the grain for green program on ecosystem health in the upper reaches of the Yangtze river basin: A case study of eastern Sichuan, China. Int. J. Environ. Res. Public Heal. 16, 2112. https://doi.org/10.3390/ijerph16122112 (2019).Article 

Google Scholar 
Zhang, X., Liu, K., Li, X., Wang, S. & Wang, J. Vulnerability assessment and its driving forces in terms of ndvi and gpp over the loess plateau, China. Phys. Chem. Earth Parts A/B/C 125, 103106. https://doi.org/10.1016/j.pce.2022.103106 (2022).Article 

Google Scholar 
Liu, L., Yan, J. & Li, S. Spatial-temporal characteristics of vegetation restoration in Qinghai Province from 2000 to 2009. Bull. Soil Water Conserv. 34, 263–267 (2014).ADS 
CAS 

Google Scholar 
Shao, E. A. Q. Target-based assessment on effects of first-stage ecological conservation and restoration project in three-river source region, china and policy recommendations. Bull. Chin. Acad. Sci. (Chin. Vers.) 32, 35–44 (2017).
Google Scholar 
Guo, E. A. J. The dynamic evolution of the ecological footprint and ecological capacity of Qinghai Province. Sustainability 12, 3065. https://doi.org/10.3390/su12073065 (2020).Article 

Google Scholar 
Jiang, W., Yihe, L., Yuanxin, L. & Wenwen, G. Ecosystem service value of the Qinghai-Tibet plateau significantly increased during 25 years. Ecosyst. Serv. 44, 101146. https://doi.org/10.1016/j.ecoser.2020.101146 (2020).Article 

Google Scholar 
Xie, S. W. W. Research on the coupling coordination between economic development and ecological environment—a case study of ecological civilization construction of Qinghai Province. Plateau Sci. Res. 4, 36–45 (2020).
Google Scholar 
Jielan, L., Xingpeng, C., Yu, W. & Zilong, Z. Research on the sustainable development in Qinghai province based on systemdynamics. Resour. Sci. 31, 1624–1631 (2009).
Google Scholar 
Tabutin, E. T. D. The relationships between population growth and environment: From doctrinal to empirical. Revuet.-monde 33, 273–294 (1992).Article 
CAS 

Google Scholar 
Fu, E. A. G. Impact of the grain for green program on forest cover in China. J. Environ. Econ. Policy 8, 231–249. https://doi.org/10.1080/21606544.2018.1552626 (2019).Article 

Google Scholar 
Jia, E. A. X. The tradeoff and synergy between ecosystem services in the grain-for-green areas in northern Shaanxi, China. Ecol. Indic. 43, 103–113. https://doi.org/10.1016/j.ecolind.2014.02.028 (2014).Article 

Google Scholar 
Deng, L., Guo-bin, L. & Zhou-ping, S. Land-use conversion and changing soil carbon stocks in China’s ‘grain-for-green’program: A synthesis. Glob. Chang. Biol. 20, 3544–3556. https://doi.org/10.1111/gcb.12508 (2014).Article 
ADS 
PubMed 

Google Scholar 
Deng, E. A. L. Past and future carbon sequestration benefits of China’s grain for green program. Glob. Environ. Chang. 47, 13–20. https://doi.org/10.1016/j.gloenvcha.2017.09.006 (2017).Article 

Google Scholar 
Cao, E. A. S. Payoff of the grain for green policy. J. Appl. Ecol. 57, 1180–1188. https://doi.org/10.1111/1365-2664.13608 (2020).Article 

Google Scholar 
Guo, H., Li, B., Hou, Y., Lu, S. & Nan, B. Rural households’ willingness to participate in the grain for green program again: A case study of Zhungeer, China. For. Policy Econ. 44, 42–49. https://doi.org/10.1016/j.forpol.2014.05.002 (2014).Article 

Google Scholar 
Li, E. A. Y. Coupling between the grain for green program and a household level-based agricultural eco-economic system in Ansai, Shaanxi Province of China. J. Arid Land 12, 199–214. https://doi.org/10.1007/s40333-020-0060-3 (2020).Article 

Google Scholar 
Li, Y., Feng, Y., Guo, X. & Peng, F. Changes in coastal city ecosystem service values based on land use—a case study of Yingkou, China. Land Use Policy 65, 287–293. https://doi.org/10.1016/j.landusepol.2017.04.021 (2017).Article 

Google Scholar 
Peng, E. A. J. Ecosystem services response to urbanization in metropolitan areas: Thresholds identification. Sci. Total. Environ. 607, 706–714. https://doi.org/10.1016/j.scitotenv.2017.06.218 (2017).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Costanza, E. A. R. The value of the world’s ecosystem services and natural capital. Nature 387, 253–260. https://doi.org/10.1038/387253a0 (1997).Article 
ADS 
CAS 

Google Scholar 
Wu, X., Wang, S., Fu, B., Liu, Y. & Zhu, Y. Land use optimization based on ecosystem service assessment: A case study in the yanhe watershed. Land Use Policy 72, 303–312. https://doi.org/10.1016/j.landusepol.2018.01.003 (2018).Article 

Google Scholar 
Zhou, E. A. Y. Land use-driven changes in ecosystem service values and simulation of future scenarios: A case study of the Qinghai–Tibet plateau. Sustainability 13, 4079. https://doi.org/10.3390/su13074079 (2021).Article 

Google Scholar 
Han, W. S. X. Z. Responses of ecosystem service to land use change in Qinghai province. Energies 9, 303. https://doi.org/10.1016/j.landusepol.2018.01.003 (2016).Article 

Google Scholar 
Shooshtari, S. J., Shayesteh, K., Gholamalifard, M., Azari, M. & López-Moreno, J. I. Land cover change modelling in hyrcanian forests, northern Iran: A landscape pattern and transformation analysis perspective. Cuader. De Investig. Geogr. 44, 743–761. https://doi.org/10.18172/cig.3279 (2018).Article 

Google Scholar 
Boerema, A., Rebelo, A. J., Bodi, M. B., Esler, K. J. & Meire, P. Are ecosystem services adequately quantified?. J. Appl. Ecol. 54, 358–370. https://doi.org/10.1111/1365-2664.12696 (2017).Article 

Google Scholar 
Costanza, E. A. R. Changes in the global value of ecosystem services. Glob. Environ. Change 26, 152–158. https://doi.org/10.1016/j.gloenvcha.2014.04.002 (2014).Article 

Google Scholar 
Newton, E. A. A. Assessing, quantifying and valuing the ecosystem services of coastal lagoons. J. For. Nat. Conserv. 44, 50–65. https://doi.org/10.1016/j.jnc.2018.02.009 (2018).Article 

Google Scholar 
Wang, L.-J., Ma, S., Zhao, Y.-G. & Zhang, J.-C. Ecological restoration projects did not increase the value of all ecosystem services in northeast china. For. Ecol. Manag. 495, 119340. https://doi.org/10.1016/j.foreco.2021.119340 (2021).Article 

Google Scholar 
Sun, X. & Li, F. Spatiotemporal assessment and trade-offs of multiple ecosystem services based on land use changes in zengcheng, china. Sci. Total Env. 609, 1569–1581. https://doi.org/10.1016/j.scitotenv.2017.07.221 (2017).Article 
CAS 

Google Scholar 
Aulia, A., Sandhu, H. & Millington, A. Quantifying the economic value of ecosystem services in oil palm dominated landscapes in Riau Province in Sumatra, Indonesia. Land 9, 194. https://doi.org/10.3390/land9060194 (2020).Article 

Google Scholar 
Peng, E. A. J. Simulating the impact of grain-for-green programme on ecosystem services trade-offs in northwestern Yunnan, China. Ecosyst. Serv. 39, 100998 (2019).Article 

Google Scholar 
Shi, Y., Shi, D., Zhou, L. & Fang, R. Identification of ecosystem services supply and demand areas and simulation of ecosystem service flows in Shanghai. Ecol. Indic. 115, 106418. https://doi.org/10.1016/j.ecolind.2020.106418 (2020).Article 
CAS 

Google Scholar 
Zoderer, B. M., Tasser, E., Carver, S. & Tappeiner, U. Stakeholder perspectives on ecosystem service supply and ecosystem service demand bundles. Ecosyst. Serv. 37, 100938. https://doi.org/10.1016/j.ecoser.2019.100938 (2019).Article 

Google Scholar 
Bryan, B. A., Ye, Y., Zhang, J. & Connor, J. D. Land-use change impacts on ecosystem services value: Incorporating the scarcity effects of supply and demand dynamics. Ecosyst. services 32, 144–157. https://doi.org/10.1016/j.ecoser.2018.07.002 (2018).Article 

Google Scholar 
Xiaojuan, Q. & Yufen, T. Coordinative development between population, economy, resources and environment in north-west area of china. China Popul. Resour. Environ. 18, P110-114 (2008).
Google Scholar 
Haiyang, Z., Zhang, Z. & Zhang, P. Rs-and gis-based evaluation and dynamic monitoring of land desertification in Qinghai Province. Arid Zone Res. 24, 153–158 (2007).
Google Scholar 
Kang, B. et al. Research on grassland ecosystem service value in china under climate change based on meta-analysis: A case study of Qinghai Province. Int. J. Clim. Chang. Strateg. Manag. https://doi.org/10.1108/ijccsm-06-2020-0073 (2020).Article 

Google Scholar 
Wang, X. & Zang, Y. Carbon footprint of the agricultural sector in Qinghai Province, China. Appl. Sci. 9, 2047. https://doi.org/10.3390/app9102047 (2019).Article 
CAS 

Google Scholar 
Wei, E. A. W. The dynamic analysis and comparison of emergy ecological footprint for the qinghai–tibet plateau: A case study of Qinghai Province and Tibet. Sustainability 11, 5587. https://doi.org/10.3390/su11205587 (2019).Article 

Google Scholar 
Chen, E. A. W. An emergy accounting based regional sustainability evaluation: A case of Qinghai in China. Ecol. Indic. 88, 152–160. https://doi.org/10.1016/j.ecolind.2017.12.069 (2018).Article 

Google Scholar 
National Bureau of Statistics of China. Chinese Statistical Yearbook 2021 (China Statistics Press, 2021).
Google Scholar 
Statistics Bureau of Qinghai Province. Qinghai Statistical Yearbook 2021 (China Statistics Press, 2021).
Google Scholar 
Zhang, J. & Ren, Z. Spatiotemporal pattern and terrain gradient effect of land use change in Qinling-Bashan mountains. Trans. Chin. Soc. Agric. Eng. 32, 250–257 (2016).
Google Scholar 
Liao, K. The discussion and prospect for geo-informatic Tupu. Geo-inf. Sci. 1, 14–20 (2002).ADS 

Google Scholar 
Lu, X., Shi, Y., Chen, C. & Yu, M. Monitoring cropland transition and its impact on ecosystem services value in developed regions of china: A case study of Jiangsu Province. Land Use Policy 69, 25–40. https://doi.org/10.1016/j.landusepol.2017.08.035 (2017).Article 

Google Scholar 
Lyu, X., Shi, Y. Y., Huang, X. J., Sun, X. F. & Miao, Z. W. Geo-spectrum characteristics of land use change in Jiangsu Province, China. The J. Appl. Ecol. 27, 1077–1084. https://doi.org/10.13287/j.1001-9332.201604.006 (2016).Article 

Google Scholar 
Dong, S., Zhao, Y. & Li, X. Spatiotemporal patterns of land use change in plateau region based on the terrain gradient—a case study in Panxian county, Guizhou Province. Res. Soil Water Conserv. 24, 213–222 (2017).
Google Scholar 
Oneill, R. V., Riitters, K. H., Wickham, J. D. & Jones, K. B. Landscape pattern metrics and regional assessment. Ecosyst. Health 5, 225–233. https://doi.org/10.1046/j.1526-0992.1999.09942.x (1999).Article 

Google Scholar 
Xie, G. D., Zhang, C. X., Zhang, L. M., Chen, W. H. & Li, S. M Improvement of the evaluation method for ecosystem service value based on per unit area. J. Nat. Resour. 30, 1243 (2015).
Google Scholar 
Wang, Y., Erfu, D., Le, Y. & Liang, M. Land use/land cover change and the effects on ecosystem services in the Hengduan mountain region, China. Ecosyst. Serv. 34, 55–67. https://doi.org/10.1016/j.ecoser.2018.09.008 (2018).Article 

Google Scholar 
Li, R., Shi, Y., Feng, C.-C. & Guo, L. The spatial relationship between ecosystem service scarcity value and urbanization from the perspective of heterogeneity in typical arid and semiarid regions of china. Ecol. Indic. 132, 108299 (2021).Article 

Google Scholar 
Shi, Y., Feng, C.-C., Yu, Q. & Guo, L. Integrating supply and demand factors for estimating ecosystem services scarcity value and its response to urbanization in typical mountainous and hilly regions of south china. Sci. Total. Environ. 796, 149032 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Liu, C. F., Li, J. Z., Li, X. M., He, X. Y. & Chen, W. Selection of landscape metrics for urban forest based on simulated landscapes. J. Appl. Ecol. 20, 1125–1131 (2009).
Google Scholar 
Zhao, L., Fan, X., Lin, H. & Hong, T. W Hong Impact of expressways on land use changes, landscape patterns, and ecosystem services value in Nanping city, China. Pol. J. Environ. Stud. 30, 2935–2946. https://doi.org/10.15244/pjoes/128584 (2021).Article 

Google Scholar 
Xu, W., Dong, X. & Zhang, Z. Spatiotemporal scale effect of vegetation landscape pattern in Saihanba area. J. North-East For. Univ. 49, 106 (2021).
Google Scholar 
Schmidt, K., Sachse, R. & Walz, A. Current role of social benefits in ecosystem service assessments. Landsc. Urban Plan. 149, 49–64. https://doi.org/10.1016/j.landurbplan.2016.01.005 (2016).Article 

Google Scholar 
Akber, M., Khan, M., Islam-M, R., Munsur, R. & Mohammad, A. Impact of land use change on ecosystem services of southwest coastal bangladesh. J. land Use Science 13, 238–250. https://doi.org/10.1080/1747423x.2018.1529832 (2018).Article 

Google Scholar 
Wang, E. A. & Xiaobin, C. Linking land use change, ecosystem services and human well-being: A case study of the manas river basin of Xinjiang, China. Ecosyst. Serv. 27, 113–123. https://doi.org/10.1016/j.ecoser.2017.08.013 (2017).Article 

Google Scholar 
Ouyang, L. T., Xiao, W. & Yongh, L. X. Spatial interaction between urbanization and ecosystem services in chinese urban agglomerations. Land Use Policy 109, 105587. https://doi.org/10.1016/j.landusepol.2021.105587 (2021).Article 

Google Scholar 
Chen, G., Chi, J. & Li, W. X. The spatial association of ecosystem services with land use and land cover change at the county level in China, 1995–2015. Sci. Total. Environ. 669, 459–470. https://doi.org/10.1016/j.scitotenv.2019.03.139 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Singh, E. A. D. Distinct influences of land cover and land management on seasonal climate. J. Geophys. Res. Atmos. 123, 12017. https://doi.org/10.1029/2018JD028874 (2018).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Diffenbaugh, N. S. Influence of modern land cover on the climate of the united states. Clim. Dyn. 33, 945–958. https://doi.org/10.1007/s00382-009-0566-z (2009).Article 

Google Scholar  More

Hydrochemistry and methane oxidation ratesThe water column of the deep North Basin is considered meromictic (i.e., permanently stratified). At the time of sampling for methane oxidation rate measurements in November 2016, the redox transition zone extended from 79 to 105 m depth; as defined here, with an upper boundary set at O2  More

Herbert, E. R. et al. A global perspective on wetland salinization: Ecological consequences of a growing threat to freshwater wetlands. Ecosphere. https://doi.org/10.1890/es14-00534.1 (2015).Article 

Google Scholar 
Kelly, V. R. et al. Long-term sodium chloride retention in a rural watershed: Legacy effects of road salt on streamater concentration. Environ. Sci. Tech. 42, 410–415. https://doi.org/10.1021/es071391l (2008).Article 
CAS 

Google Scholar 
Tiwari, A. & Rachlin, J. W. A review of road salt ecological impacts. Northeast. Nat. 25, 123–142. https://doi.org/10.1656/045.025.0110 (2018).Article 

Google Scholar 
Hintz, W. D. & Relyea, R. A. A review of the species, community, and ecosystem impacts of road salt salinisation in fresh waters. Freshwater Biol. 64, 1081–1097. https://doi.org/10.1111/fwb.13286 (2019).Article 

Google Scholar 
Dugan, H. A. et al. Salting our freshwater lakes. Proc. Natl. Acad. of Sci. U.S.A 114, 4453–4458. https://doi.org/10.1073/pnas.1620211114 (2017).Article 
ADS 
CAS 

Google Scholar 
Kaushal, S. S. et al. Increased salinization of fresh water in the northeastern United States. Proc. Natl. Acad. of Sci. U.S.A. 102, 13517–13520. https://doi.org/10.1073/pnas.0506414102 (2005).Article 
ADS 
CAS 

Google Scholar 
Kaushal, S. S. et al. Freshwater salinization syndrome: from emerging global problem to managing risks. Biogeochemistry 154, 255–292. https://doi.org/10.1007/s10533-021-00784-w (2021).Article 

Google Scholar 
Kaushal, S. S. et al. Freshwater salinization syndrome on a continental scale. Proc. Natl. Acad. of Sci. U.S.A. 115, E574–E583. https://doi.org/10.1073/pnas.1711234115 (2018).Article 
CAS 

Google Scholar 
Hintz, W. D., Fay, L. & Relyea, R. A. Road salts, human safety, and the rising salinity of our fresh waters. Front. Ecol. Environ. 9, 22–30. https://doi.org/10.1002/fee.2433 (2022).Article 

Google Scholar 
Petranka, J. W. & Doyle, E. J. Effects of road salts on the composition of seasonal pond communities: Can the use of road salts enhance mosquito recruitment?. Aquat. Ecol. 44, 155–166. https://doi.org/10.1007/s10452-009-9286-z (2010).Article 
CAS 

Google Scholar 
Petranka, J. W. & Francis, R. A. Effects of road salts on seasonal wetlands: Poor prey performance may compromise growth of predatory salamanders. Wetlands 33, 707–715. https://doi.org/10.1007/s13157-013-0428-7 (2013).Article 

Google Scholar 
Searle, C. L., Shaw, C. L., Hunsberger, K. K., Prado, M. & Duffy, M. A. Salinization decreases population densities of the freshwater crustacean Daphnia dentifera. Hydrobiologia 770, 165–172. https://doi.org/10.1007/s10750-015-2579-4 (2016).Article 
CAS 

Google Scholar 
Hebert, M. P. et al. Lake salinization drives consistent losses of zooplankton abundance and diversity across coordinated mesocosm experiments. Limnol. Oceanogr. Let. https://doi.org/10.1002/lol2.10239 (2022).Article 

Google Scholar 
Collins, S. J. & Russell, R. W. Toxicity of road salt to nova scotia amphibians. Environ. Pollut. 157, 320–324. https://doi.org/10.1016/j.envpol.2008.06.032 (2009).Article 
CAS 
PubMed 

Google Scholar 
Milotic, D., Milotic, M. & Koprivnikar, J. Effects of road salt on larval amphibian susceptibility to parasitism through behavior and immunocompetence. Aquat. Toxicol. 189, 42–49. https://doi.org/10.1016/j.aquatox.2017.05.015 (2017).Article 
CAS 
PubMed 

Google Scholar 
Sanzo, D. & Hecnar, S. J. Effects of road de-icing salt (NaCl) on larval wood frogs (Rana sylvatica). Environ. Pollut. 140, 247–256. https://doi.org/10.1016/j.envpol.2005.07.013 (2006).Article 
CAS 
PubMed 

Google Scholar 
Arnott, S. E. et al. Road salt impacts freshwater zooplankton at concentrations below current water quality guidelines. Envir. Sci. Tech. 54, 9398–9407. https://doi.org/10.1021/acs.est.0c02396 (2020).Article 
CAS 

Google Scholar 
Elphick, J. R. F., Bergh, K. D. & Bailey, H. C. Chronic toxicity of chloride to freshwater species effects of hardness and implications for water quality guidelines. Environ. Toxicol. Chem. 30, 239–246. https://doi.org/10.1002/etc.365 (2011).Article 
CAS 
PubMed 

Google Scholar 
Mount, D. R. et al. The acute toxicity of major ion salts to Ceriodaphnia dubia: I. Influence of background water chemistry. Environ. Toxicol. Chem. 35, 3039–3057. https://doi.org/10.1002/etc.3487 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Soucek, D. J. Comparison of hardness- and chloride-regulated acute effects of sodium sulfate on two freshwater crustaceans. Environ. Toxicol. Chem. 26, 773–779. https://doi.org/10.1897/06-229r.1 (2007).Article 
CAS 
PubMed 

Google Scholar 
Bhateria, R. & Jain, D. Water quality assessment of lake water: A review. Sustain. Wat. Res. Manag. 2, 161–173. https://doi.org/10.1007/s40899-015-0014-7 (2016).Article 

Google Scholar 
USGS. Hardness of Water. https://www.usgs.gov/special-topics/water-science-school/science/hardness-water#overview, Accessed: 1 August 2022 (2018).Brown, A. H. & Yan, N. D. Food quantity affects the sensitivity of Daphnia to Road Salt. Environ. Sci. Tech. 49, 4673–4680. https://doi.org/10.1021/es5061534 (2015).Article 
CAS 

Google Scholar 
Smith, D. W. & Cooper, S. D. Competition among cladocera. Ecology 63, 1004–1015. https://doi.org/10.2307/1937240 (1982).Article 

Google Scholar 
Soucek, D. J. et al. Influence of water hardness and sulfate on the acute toxicity of chloride to sensitive freshwater invertebrates. Environ. Toxicol. Chem. 30, 930–938. https://doi.org/10.1002/etc.454 (2011).Article 
CAS 
PubMed 

Google Scholar 
Gust, K. A. et al. Daphnia magna’s sense of competition: Intra-specific interactions (ISI) alter life history strategies and increase metals toxicity. Ecotoxicology 25, 1126–1135. https://doi.org/10.1007/s10646-016-1667-1 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Liu, X. & Steiner, C. F. Ecotoxicology of salinity tolerance in Daphnia pulex: Interactive effects of clonal variation, salinity stress and predation. J. Plankton Res. 39, 687–697. https://doi.org/10.1093/plankt/fbx027 (2017).Article 
CAS 

Google Scholar 
Evans, M. & Frick, C. The effects of road salts on aquatic ecosystems. Report No. 02-308, (Environment Canada – Water Science and Technology Directorate, 2001).USEPA. (U.S. Environmental Protection Agency, 1988).Schuler, M. S. et al. Regulations are needed to protect freshwater ecosystems from salinization. Phil. Trans. R. Soc. B. https://doi.org/10.1098/rstb.2018.0019 (2019).Article 

Google Scholar 
Canadian Council of Ministers for the Environment. Candadian water Quality Guidelines for the Protection of Aquatic Life: Chloride. (Environment Canada, Gatineau, Canada, 2011).Valleau, R. E., Paterson, A. M. & Smol, J. P. Effects of road-salt application on Cladocera assemblages in shallow precambrian shield lakes in south-central Ontario, Canada. Freshwat. Sci. 39, 824–836. https://doi.org/10.1086/711666 (2020).Article 

Google Scholar 
Hintz, W. D. et al. Current water quality guidelines across North America and Europe do not protect lakes from salinization. Proc. Natl. Acad. of Sci. U.S.A. https://doi.org/10.1073/pnas.2115033119 (2022).Article 

Google Scholar 
Valleau, R. E., Celis-Salgado, M. P., Arnott, S. E., Paterson, A. M. & Smol, J. P. Assessing the effect of salinization (NaCl) on the survival and reproduction of two ubiquitous cladocera species (Bosmina longirostris and Chydorus brevilabris). Wat. Air Soil Pollut. 233, 135. https://doi.org/10.1007/s11270-021-05482-9 (2022).Article 
ADS 
CAS 

Google Scholar 
Celis-Salgado, M. P., Cairns, A., Kim, N. & Yan, N. D. The FLAMES medium: A new, soft-water culture and bioassay medium for Cladocera. SIL Proc. 1922–2010(30), 265–271. https://doi.org/10.1080/03680770.2008.11902123 (2008).Article 

Google Scholar 
USEPA. Methods for measuring the acute toxicity of effluents and receiving waters to freshwater and marine organisms, 5th ed. Office of Water Washington, DC (2002).Hintz, W. D. et al. Concurrent improvement and deterioration of epilimnetic water quality in an oligotrophic lake over 37 years. Limnol. Oceanogr. 65, 927–938. https://doi.org/10.1002/lno.11359 (2020).Article 
ADS 
CAS 

Google Scholar 
Winner, R. W. Interactive effects of water hardness and humic acid on the chronic toxicity of cadmium to Daphnia pulex. Aquat. Toxicol. 8, 281–293. https://doi.org/10.1016/0166-445X(86)90080-9 (1986).Article 
CAS 

Google Scholar 
Kaushal, S. S. et al. Novel “chemical cocktails” in inland waters are a consequence of the freshwater salinization syndrome. Phil. Trans. R. Soc. B. https://doi.org/10.1098/rstb.2018.0017 (2019).Article 

Google Scholar 
Kaushal, S. S. et al. Making “chemical cocktails”: Evolution of urban geochemical processes across the periodic table of elements. Appl. Geochem. https://doi.org/10.1016/j.apgeochem.2020.104632 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Cremona, F. et al. How warming and other stressors affect zooplankton abundance, biomass and community composition in shallow eutrophic lakes. Clim. Change 159, 565–580. https://doi.org/10.1007/s10584-020-02698-2 (2020).Article 
ADS 
CAS 

Google Scholar 
Lind, L. et al. Salty fertile lakes: How salinization and eutrophication alter the structure of freshwater communities. Ecosphere. https://doi.org/10.1002/ecs2.2383 (2018).Article 

Google Scholar 
Stoler, A. B. et al. Effects of a common insecticide on wetland communities with varying quality of leaf litter inputs. Environ. Pollut. 226, 452–462. https://doi.org/10.1016/j.envpol.2017.04.019 (2017).Article 
CAS 
PubMed 

Google Scholar 
Riessen, H. P. & Sprules, W. G. Demographic costs of antipredator defenses in Daphnia pulex. Ecology 71, 1536–1546. https://doi.org/10.2307/1938290 (1990).Article 

Google Scholar  More

Concentrations of PTEs in agricultural soil of oasesThe descriptive statistical results of 11 PTE concentrations in agricultural soils of the four oases are shown in Table 1. The concentration order of different PTEs in the four oases was different. In terms of mean values, the concentrations of Co, Cu, Ni, Pb, Sb, V and Zn in Kashgar Oasis agricultural soil were higher than those in the other three oases, and the concentrations of Sn and Tl in Hotan Oasis agricultural soil were the highest. The Cv value of As in the Hotan Oasis was the highest (0.31), indicating that compared with that in other oases and other PTEs, the As in Hotan Oasis soil was more likely to have external inputs other than natural sources59. According to the “Soil Environmental Quality · Agricultural Land Soil Pollution Risk Control Standard (Trial)” (GB15618-2018)60, the corresponding elements of oasis agricultural soil in the source region of the Tarim River did not exceed the pollution risk control standard. Compared with the soil environmental quality provisions in the “Environmental Quality Assessment Standard for Producing Areas of Edible Agricultural Products” (HJ332-2006)61, the corresponding elements in the soil of the four oases were also below the limits.Table 1 Descriptive statistics of PTEs of oasis agricultural soil in the source region of the Tarim River.Full size tableCompared with the PTE concentrations in the agricultural soil of oases in other arid areas, the concentrations of Cd, Cu and V in the soil of the study area were lower than those of the Bortala River Basin and Ili River Basin3 in the northern Tianshan Mountains, while the concentrations of Co, Ni, Pb and Zn were similar to those of corresponding elements in agricultural soil of the Bortala River Basin and Ili River Basin3. Compared with the concentrations of PTEs in agricultural soil of Wuwei and Jiuquan cities in the Hexi Corridor62,63, the concentrations of Cu, Ni, Pb and V elements in agricultural soil of the four oases in the study area were lower. Similarly, the concentrations of Cu, V and Zn in the study area were much lower than the corresponding concentrations in the agricultural soil of the Bahariya Oasis, Egypt64.Geochemical baseline values of PTEs in oasis agricultural soilsThe PTE concentration of the subsoil (60–80 cm) after logarithmic and reciprocal conversion passed the KS normal distribution test (p  > 0.05). The GBVs of PTEs in agricultural soils of the four oases in the source region of the Tarim River were calculated according to the cumulative frequency method, and the results are shown in Supplementary Figs. S1–4 and Table 2. The GBVs of PTEs in agricultural soils were very different among the four oases. For example, except for those of Ni, Sn, Tl and V, the GBVs of PTEs in agricultural soils of the Hotan Oasis were much lower than those of the other three oases (Table 2). There was also a large gap between the GBVs of PTEs in agricultural soil and the soil background values in China65, Xinjiang66,67 and Xinjiang agricultural land66. In particular, the GBV of Zn obtained in this study was 3.6–5 times the soil background value of Xinjiang agricultural land. Under the same land use type, there are differences in the abundance of environmental PTEs due to heterogeneity in the disturbance degree of human activities and the parent material and soil formation processes68. Therefore, there were some differences in the GBVs of each PTE among the four oases in this study when the land use types were all agricultural soils.Table 2 GBVs of PTEs in oasis agricultural soils of the source region of the Tarim River.Full size tablePollution characteristics of PTEs in agricultural soil of oasesAccording to the GBVs of the above elements, the PF, PLI and PN of PTEs in agricultural soil of the four oases were calculated, and the results are shown in Fig. 2. The results of PF showed that (Fig. 2a), except for As and Cd in the Hotan Oasis, the PTEs of agricultural soils in the four oases were at uncontaminated and slightly contaminated levels, indicating that the oasis agricultural soils in the source region of the Tarim River were affected by human activities to a certain extent. From the perspective of individual elements, the PF of Cd in agricultural soils of the four oases was higher than that of other PTEs. Among them, Cd in most soil samples from the Hotan Oasis (81.25%) showed mild contamination, while Cd in the remaining samples (18.75%) showed moderate contamination. Since Cd in soil is generally regarded as an indicator of agricultural activities involving fertilizer use, the results suggest that soil in oases (especially the Hotan Oasis) in the source region of the Tarim River may be seriously affected by agricultural activities.Figure 2Box-whisker plots of PF (a), PLI (b) and PN (c) of PTEs in oasis agricultural soils.Full size imageThe results of the pollution load index of multiple PTEs showed that the agricultural soils in the Yarkant River Oasis, Kashgar Oasis and Aksu Oasis were between uncontaminated and moderate contamination levels, and most soil samples in the Hotan Oasis were moderately contaminated. PN showed similar results. The reason for this result was the high PF values of As and Cd in the agricultural soil of the Hotan Oasis (Fig. 2a), which caused the high comprehensive pollution index in this area. It is worth noting that compared with those in the other three areas, most PTEs (As, Cd, Pb, Sb, V and Zn) in the agricultural soil of the Hotan Oasis had lower mean values, but their PLI and PN values were higher. This was closely related to the lower GBVs of PTEs in the Hotan area.Source apportionment and source specific risk assessment of PTEs in oasis agricultural soilSource apportionment of soil PTEsFactor analysis and correlation analysis were used to identify the main sources of PTEs in the main oasis agricultural soil in the source region of the Tarim River. The results are shown in Table 3 and Supplementary Fig. S5. Table 3 shows that PTEs in the agricultural soil of the Yarkant River Oasis had two principal factors. The first principal component F1, which explained 61.6% of the total variance, mainly described Co, V, Cu and Sb and moderately described Cd, Ni and Tl. The second principal component, F2, accounted for 21.8% of the variance and mainly described As, Pb and Sn. There were three main factors for PTEs in Kashgar Oasis agricultural soil. Factor 1, which explained 58.9% of the total variance, was mainly composed of As, Co, Cu, Ni, Pb, Sb, V and Zn; factor 2, which accounted for 24.8% of the variance, was mainly composed of Sn and Tl; and factor 3, which accounted for 7.9% of the variance, was composed of Cd. There were three main factors for the PTEs of Aksu Oasis agricultural soil. Factor 1 was mainly composed of Cd, Co, Cu, Ni, Pb, V and Zn; factor 2 was mainly composed of Sn and Tl; and factor 3 was composed of As and Sb. PTEs in agricultural soil of the Hotan Oasis had only one main factor, which was composed of all (11) elements. To verify the above results of PTE extraction from oasis agricultural soil by PCA, Spearman correlation analysis was performed on PTEs from the soils of the four oases. The results showed that the correlations between PTEs in agricultural soils of the four oases were consistent with the results of PCA, indicating the reliability of the PCA results.Table 3 Principal component analysis results of PTEs of oasis agricultural soil in the source region of the Tarim River.Full size tableTo further determine the specific sources and contribution rates of PTEs in oasis agricultural soil in the source region of the Tarim River, the PMF model was adopted for analysis, and the results are shown in Fig. 3. Overall, the classification results of PTEs in agricultural soils of the four oases by the PMF model were consistent with the results of PCA. According to the above results, eight environmental factors related to the source of soil PTEs were selected: distance from factory (DF), distance from road (DR), pH, soil type (ST), total nitrogen (TN), soil fine silty particle size percentage (Fine silty), soil silty particle size percentage (Silty), and soil coarse silty particle size percentage (Coarse silty). A Geodetector model of PTEs and environmental factors in soil was constructed, and the results are shown in Supplementary Table S5. Among them, the TN content in soil was considered to be an index related to the intensity of agricultural activities69. Since the particle composition of atmospheric dust in the study area is mainly silt70, the silty particle size content of soil was selected as the atmospheric dust index in this study.Figure 3PMF analysis results of PTEs of agricultural soils in the Yarkant River Oasis (a), Kashgar Oasis (b) and Aksu Oasis (c).Full size imageIn the Yarkant River Oasis, Co, Cu and Zn showed no pollution at most of the sampling points for the first source factor of PTEs in agricultural soil (Fig. 2), indicating that the levels of these elements in soil were less affected by human activities. Meanwhile, the results of the factor detector showed that ST was one of the factors explaining the spatial distribution of Co, Cu, Sb, V and Zn (Supplementary Table S5). The Geodetector results revealed TN as one of the strong drivers of Cd and Tl, the former of which is found in pesticides and fertilizers and the latter of which may also enter agricultural soil through sewage irrigation and fertilizers contaminated by industrial wastewater71,72. Therefore, it is inferred that F1 in the agricultural soil of the Yarkant River Oasis had a mixed source of natural resources and agricultural activities. The three elements in F2 (As, Pb and Sn) were closely related to fuel combustion and traffic factors, among which Pb is usually selected as an identifying element for traffic sources (including leaded exhaust gas, vehicle tires and brake pads)73. In the results of the Geodetector model, the best explanatory factors for Pb and Sn were Silt and DR, while the best explanatory factor for As was Slit. Therefore, F2 was inferred to be the source of road dust/atmospheric dust.The weights of F1 in Kashgar Oasis agricultural soil PTEs were mainly included As, Co, Cu, Ni, Pb, Sb, V and Zn; the weights of F2 were mainly included Sn and Tl; and F3 were mainly included Cd (Fig. 3b). Similar to the results observed for the Yarkant River Oasis, the factor detector results showed that F1 might have a mixed anthropogenic source, including transportation, industrial and other sources. Both Sn and Tl in F2 are widely used in industrial manufacturing, and Sn can be used as an additive to enhance the properties of steel or alloys74. Tl is used in many different industrial manufacturing and medical fields, and metal smelting, sulfuric acid production, coal burning, cement manufacturing and other industrial activities involving the use of Tl minerals are the main pathways by which this element enters the environment71. At the same time, the strongest explanatory factor for Sn and Tl in the results of the geographic detector was DF, so F2 was inferred to be an industrial source. F3 contained only Cd (related to agricultural activities), so F3 was inferred to be the source of agricultural activities (Supplementary Table S5).In the PMF results, the weights of F1 for the PTEs of Aksu Oasis agricultural soil were mainly included Cu, Ni, Cd, Zn, Co, V and Pb; the weights of F2 were mainly included Sn and Tl; and the weights of F3 were mainly included As and Sb (Fig. 3c). Similarly, from the factor detector calculation results, it was inferred that F1 in the Aksu Oasis was a mixed source. According to the results of Geodetector analysis, the best explanatory factors for Sn and Tl in the agricultural soil of the Aksu Oasis were TN and ST. Previous studies have also shown that Sn in soil may also come from agricultural practices (pesticides)75, so F2 was inferred to represent agricultural activities and natural sources (rock mineralization). Similarly, according to Supplementary Table S5, F3 was inferred to be the source of agricultural activities. Previous studies have also concluded that agricultural activities such as the application of phosphate fertilizer are also the main source of As and Sb in soil76,77, which was consistent with the results of the Geodetector model.The explanatory factors for the PTEs of agricultural soil in the Hotan Oasis were the total nitrogen content and the silty particle size percentage, indicating that the PTEs in agricultural soil in the Hotan Oasis mainly came from agricultural activities and atmospheric dust fall. The economy in the Hotan Oasis is dominated by irrigated agriculture, and since most oases border large deserts, wind and sand disasters are extremely serious (annual average dust days exceeding 220 days)78. The sources of PTEs inferred by the Geodetector model in this study were consistent with the reality.Source-specific ecological risk assessmentBased on the results of source analysis, the source-specific ecological risks posed by agricultural soil PTEs were evaluated in this study, and the results are shown in Fig. 4. The total ecological risks caused by PTEs in the agricultural soils in the Yarkant River Oasis, Kashgar Oasis and Aksu Oasis were null, while those at of all sampling sites in the Hotan Oasis were moderate risk. According to the results of the source-specific ecological risk of PTEs, there was no direct relationship between the contribution degree of source-specific risks to the total ecological risk and its contribution to the existence of PTEs in soil. In particular, the agricultural activity source (F3), which accounted for only 7.9% of the PTEs in Kashgar Oasis agricultural soil, contributed the most to the total risk, the mixed source (F1), which contributed the most to the existence of PTEs in soil, contributed the second most to the total risk, and the industrial source (F2) contributed the least. The reason for this result was that the toxicity factor of Cd (30) introduced by agricultural activities was significantly higher than that of PTEs released from other sources, which was consistent with the conclusion of other studies16,19 that the main source of highly toxic elements was more likely to cause ecological risks than that of low-toxicity elements.Figure 4Source-specific ecological risk boxplots of PTEs of agricultural soils in oases in the source region of the Tarim River.Full size imageTo explore the spatial distribution of ecological risks generated by PTEs from different sources, a distribution map of ecological risks posed by specific sources at each sampling point was generated (Fig. 5). The total ecological risk of the sampling sites near Awat County in the Aksu Oasis was higher than 40, indicating moderate risk. In addition, the total ecological risk at the other sampling sites was at a low level. However, the ecological risks caused by mixed sources (F1) and agricultural activities and natural sources (F2) were high in the southwestern Aksu Oasis and the sampling sites near Aksu city, resulting in high total ecological risks. The ecological risks at all the sampling sites in the two oases in Kashgar (Kashgar Oasis and Yarkant River Oasis) were low. Notably, the sampling sites near Kashgar city had high total ecological risks due to the high ecological risks caused by PTEs from mixed sources (F1) and agricultural activity sources (F3). The total ecological risk at all sampling sites in the Hotan Oasis was moderate. In particular, the sampling sites near Hotan city had the highest total ecological risk generated by PTEs from agricultural activities and atmospheric dust sources (73.71).Figure 5Spatial distribution of the ecological risks of PTEs from different sources in the four oases. The graphs were generated by QGIS 3.26.3 (https://www.qgis.org) and the land use data are from the ESA global 30 m resolution land use cover dataset (https://viewer.esa-worldcover.org/worldcover). The combination of graphs (a–d) was accomplished with linkscape 1.2.1 (https://inkscape.org).Full size imageSource-specific human health risk assessmentAccording to the proportions of different sources of each PTE obtained from the PMF model, the source-specific health risks of PTEs in agricultural soil in the four oases to the human body were calculated, and the results are shown in Tables 4 and 5 and Fig. 6. Although the results of noncarcinogenic risk caused by different sources of PTEs in agricultural soil in each oasis showed that adults and children were not at risk, the total THI values of PTEs for children in the four regions were all greater than 1, indicating that soil PTEs in the study area posed significant noncarcinogenic risks for children. In terms of carcinogenic risk, the total TCR values of PTEs for adults and children were on the order of 1E−05, within the range of 1E−06 and 1E−04, indicating that the carcinogenic risk of soil PTEs for the human body is acceptable in the study area. In general, the noncarcinogenic and carcinogenic risks of PTEs from different sources in children were higher than those in adults, which can be explained by children having more opportunities PTE contamination through hand-to-mouth ingestion and dermal contact than adults due to the areas where they play and unhealthy eating habits (e.g., children are more likely to suck their fingers)79. Therefore, different parameters were set when employing the health risk assessment model.Table 4 Specific noncarcinogenic risks of PTEs from different sources.Full size tableTable 5 Specific carcinogenic risks of PTEs from different sources.Full size tableFigure 6Human health risk proportions of PTEs from different sources in the Aksu Oasis (a), Kashgar Oasis (b) and Yarkant River Oasis (c).Full size imageIn terms of the proportion of human health risks caused by PTEs from different sources (Fig. 6), those from mixed sources, which were considered to have the largest contribution to PTEs in the Aksu Oasis, Kashgar Oasis and Yarkant River Oasis, were not the largest. This might be explained by the presence of more toxic elements, such as As and Tl, in other sources of the three oases contributing the most (Fig. 3), while their low RfD may explain their greater noncarcinogenic risk. PTEs from agricultural activities and natural sources accounted for the largest proportion of carcinogenic risk in the Aksu Oasis, while mixed sources accounted for the largest proportion of carcinogenic risk in the Kashgar Oasis and Yarkant River Oasis. This was mainly because the presence of As, Cd, Co, Ni and Pb elements with carcinogenic risk in the above two sources accounted for a large proportion (Fig. 3), thus posing a large carcinogenic risk. The results of the proportion of health risks attributed to different sources of PTEs to the human body showed that the control of human health risks of PTEs cannot be determined based on their concentration alone. Therefore, in view of the obvious noncarcinogenic risk of soil PTEs in the four oases to children, from the results of the contribution of different sources of PTEs to the noncarcinogenic risk to children, agricultural activities and natural sources, industrial sources and atmospheric dust fall were the priority control sources in the Aksu Oasis, Kashgar Oasis and Yarkant River Oasis, respectively. More

These authors contributed equally: Arjun Srivathsa, Divya Vasudev, Tanaya Nair, Jagdish Krishnaswamy, Uma Ramakrishnan.National Centre for Biological Science, TIFR, Bengaluru, IndiaArjun Srivathsa, Tanaya Nair, Mahesh Sankaran & Uma RamakrishnanWildlife Conservation Society-India, Bengaluru, IndiaArjun SrivathsaConservation Initiatives, Guwahati, IndiaDivya Vasudev & Varun R. GoswamiDivision of Biosciences, University College London, London, UKTanaya NairDepartments of Biology and Environmental Studies, Macalester College, Saint Paul, MN, USAStotra ChakrabartiWorld Wildlife Fund, Delhi, IndiaPranav Chanchani, Arpit Deomurari, Dipankar Ghose & Prachi ThatteDepartment of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, USARuth DeFriesAmity Institute of Forestry and Wildlife, Amity University, Noida, IndiaArpit DeomurariWildlife Institute of India, Dehradun, IndiaSutirtha DuttaFoundation for Ecological Research, Advocacy and Learning, Bengaluru, IndiaRajat Nayak & Srinivas VaidyanathanNetwork for Conserving Central India, Gurgaon, IndiaAmrita NeelakantanWorld Resources Institute, New Delhi, IndiaMadhu VermaSchool of Environment and Sustainability, Indian Institute for Human Settlements, Bengaluru, IndiaJagdish KrishnaswamyAshoka Trust for Research in Ecology and the Environment, Bengaluru, IndiaJagdish KrishnaswamyBiodiversity Collaborative, Bengaluru, IndiaJagdish Krishnaswamy, Mahesh Sankaran & Uma Ramakrishnan More