Philippa Kaur

Editorial Committee of Vegetation Map China. Vegetation Map of China (1:1000 000) (Geology Press, 2007).Fu, B. et al. Current condition and protection strategies of Qinghai–Tibet Plateau ecological security barrier. Bull. Chin. Acad. Sci. 36, 1298–1306 (2021).
Google Scholar 
Zhang, Y. et al. Spatial and temporal variability in the net primary production (NPP) of alpine grassland on Tibetan Plateau from 1982 to 2009. Acta Geogr. Sin. 68, 1197–1211 (2013).
Google Scholar 
Bao, C. & Liu, R. Spatiotemporal evolution of the urban system in the Tibetan Plateau. J. Geoinf. Sci. 21, 1330–1340 (2019).
Google Scholar 
Miehe, G. et al. The Kobresia pygmaea ecosystem of the Tibetan highlands — origin, functioning and degradation of the world’s largest pastoral alpine ecosystem: Kobresia pastures of Tibet. Sci. Total Environ. 648, 754–771 (2019). This work describes features of K. pygmaea grassland and reveals that overstocking has caused pasture degradation and soil deterioration.Article 

Google Scholar 
Liu, Y. et al. Grassland dynamics in responses to climate variation and human activities in China from 2000 to 2013. Sci. Total Environ. 690, 27–39 (2019).Article 

Google Scholar 
Cao, J. et al. Grassland degradation on the Qinghai–Tibetan Plateau: reevaluation of causative factors. Rangel. Ecol. Manag. 72, 988–995 (2019).Article 

Google Scholar 
Zhao, X. Restoration and Sustainable Management of Degradaded Grassland in the Three Rivers Headwater (Science Press, 2011).Gao, Q. Exploration and Study on Eoclogical Revelization Fuatures in Qiangtang Plateau (China Agriculture Press, 2015).Gu, X. et al. Soil extractable organic C and N contents, methanotrophic activity under warming and degradation in a Tibetan alpine meadow. Agric. Ecosyst. Environ. 278, 6–14 (2019).Article 

Google Scholar 
Li, Y. et al. Changes of soil microbial community under different degraded gradients of alpine meadow. Agric. Ecosyst. Environ. 222, 213–222 (2016).Article 

Google Scholar 
Wang, W., Wang, Q. & Wang, H. The effect of land management on plant community composition, species diversity, and productivity of alpine Kobersia steppe meadow. Ecol. Res. 21, 181–187 (2005).Article 

Google Scholar 
Xu, H., Wang, X. & Zhang, X. Alpine grasslands response to climatic factors and anthropogenic activities on the Tibetan Plateau from 2000 to 2012. Ecol. Eng. 92, 251–259 (2016).Article 

Google Scholar 
Yu, L., Tang, L., Wei, D., Mei, M. & Zhou, H. Characteristics and causes of changes of alpine grassland productivity in the source region of Yellow River. Int. Conf. Geoinformatics https://doi.org/10.1109/GEOINFORMATICS.2010.5567879 (2010).Article 

Google Scholar 
Yang, Y. et al. Responses of the functional structure of soil microbial community to livestock grazing in the Tibetan alpine grassland. Glob. Chang. Biol. 19, 637–648 (2013). This work shows that soil microbial community functional structure is very sensitive to livestock grazing.Article 

Google Scholar 
Gao, Y. Z. et al. Belowground net primary productivity and biomass allocation of a grassland in Inner Mongolia is affected by grazing intensity. Plant Soil 307, 41–50 (2008).Article 

Google Scholar 
Dlamini, A. P. et al. Controlling factors of sheet erosion under degraded grasslands in the sloping lands of KwaZulu-Natal, South Africa. Agric. Water Manag. 98, 1711–1718 (2011).Article 

Google Scholar 
Niemandt, C. & Greve, M. Fragmentation metric proxies provide insights into historical biodiversity loss in critically endangered grassland. Agric. Ecosyst. Environ. 235, 172–181 (2016).Article 

Google Scholar 
Kang, S. C. et al. Review of climate and cryospheric change in the Tibetan Plateau. Environ. Res. Lett. 5, 15101–15101 (2010).Article 

Google Scholar 
Shen, H. et al. Effects of simulated N deposition on photosynthesis and productivity of key plants from different functional groups of alpine meadow on Qinghai–Tibetan Plateau. Environ. Pollut. 251, 731–737 (2019).Article 

Google Scholar 
Yu, G. R. et al. Stabilization of atmospheric nitrogen deposition in China over the past decade. Nat. Geosci. 12, 424 (2019).Article 

Google Scholar 
Lu, C. & Tian, H. Spatial and temporal patterns of nitrogen deposition in China: synthesis of observational data. J. Geophys. Res. 112, D22S05 (2007).
Google Scholar 
National Bureau of Statistics of China. China Statistics Yearbook (China Statistics Press, 2020).Mo, X. Sustainable livestock carring capacity and overgrazing rate of grassland over Qinghai–Tibet plateau since 1980. Natl Tibetan Plateau Data Center https://doi.org/10.11888/Socioeco.tpdc.270347 (2020).Article 

Google Scholar 
Tian, Y. Y., Jiang, G. H., Zhou, D. Y. & Li, G. Y. Systematically addressing the heterogeneity in the response of ecosystem services to agricultural modernization, industrialization and urbanization in the Qinghai–Tibetan Plateau from 2000 to 2018. J. Clean. Prod. https://doi.org/10.1016/j.jclepro.2020.125323 (2021).Article 

Google Scholar 
Yao, Y. et al. Spatiotemporal pattern of gross primary productivity and its covariation with climate in China over the last thirty years. Glob. Chang. Biol. 24, 184–196 (2018).Article 

Google Scholar 
Li, L. et al. Current challenges in distinguishing climatic and anthropogenic contributions to alpine grassland variation on the Tibetan Plateau. Ecol. Evol. 8, 5949–5963 (2018). This work finds that large inconsistencies exist in distinguishing the respective contribution of climatic and anthropogenic forces to grassland dynamics.Article 

Google Scholar 
Zhong, L., Ma, Y. M., Xue, Y. K. & Piao, S. L. Climate change trends and impacts on vegetation greening over the Tibetan Plateau. J. Geophys. Res. Atmos. 124, 7540–7552 (2019). This work demonstrates that the general increasing trends in vegetation density and greening of the QTP are mainly caused by climate factors, using satellite-derived climate and vegetation data from 1999 to 2014.Article 

Google Scholar 
Yang, K. & He, J. China meteorological forcing dataset (1979–2018). Natl Tibetan Plateau Data Center https://doi.org/10.11888/AtmosphericPhysics.tpe.249369.file (2019).Article 

Google Scholar 
Xiong, Q. et al. Monitoring the impact of climate change and human activities on grassland vegetation dynamics in the northeastern Qinghai–Tibet Plateau of China during 2000–2015. J. Arid. Land 11, 637–651 (2019).Article 

Google Scholar 
Pan, T., Zou, X. T., Liu, Y. J., Wu, S. H. & He, G. M. Contributions of climatic and non-climatic drivers to grassland variations on the Tibetan Plateau. Ecol. Eng. 108, 307–317 (2017).Article 

Google Scholar 
Hou, X. 1:1 Million vegetation map of China (National Tibetan Plateau Data Center, 2019).Peng, S. S. et al. Recent change of vegetation growth trend in China. Environ. Res. Lett. 6, 044027 (2011).Article 

Google Scholar 
Yuan, W. et al. Increased atmospheric vapor pressure deficit reduces global vegetation growth. Sci. Adv. 5, eaax1396 (2019).Article 

Google Scholar 
Zhu, Z. C. et al. Greening of the earth and its drivers. Nat. Clim. Chang. 6, 791 (2016).Article 

Google Scholar 
Vermote, E. et al. NOAA climate data record (CDR) of normalized difference vegetation index (NDVI), version 4. NOAA Natl Cent. Environ. Inf. https://doi.org/10.7289/V5PZ56R6 (2014).Article 

Google Scholar 
Chen, H. et al. Attribution analyses of changes in alpine grasslands on the Qinghai–Tibetan Plateau. Chin. Sci. Bull. 65, 2406–2418 (2020). This work demonstrates that human activities play an increasingly important role in the restoration of degraded grasslands.Article 

Google Scholar 
Shen, M. et al. Evaporative cooling over the Tibetan Plateau induced by vegetation growth. Proc. Natl Acad. Sci. USA 112, 9299–9304 (2015).Article 

Google Scholar 
Cai, D. et al. Vegetation dynamics on the Tibetan Plateau (1982–2006): an attribution by ecohydrological diagnostics. J. Clim. 28, 4576–4584 (2015).Article 

Google Scholar 
Zhou, W. et al. Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 83, 303–313 (2017).Article 

Google Scholar 
Liu, Z. et al. Patterns of plant species diversity along an altitudinal gradient and its effect on above-ground biomass in alpine meadows in Qinghai–Tibet Plateau. Biodivers. Sci. 23, 451–462 (2015).Article 

Google Scholar 
Harris, R. B. Rangeland degradation on the Qinghai–Tibetan Plateau: a review of the evidence of its magnitude and causes. J. Arid. Environ. 74, 1–12 (2010).Article 

Google Scholar 
Lu, S. et al. Basic characteristics of Stipa sareptana var. krylovii communities in China. Chin. J. Plant. Ecol. 44, 1087–1094 (2020).Article 

Google Scholar 
Qiao, X. et al. Distribution, community characteristics and classification of Stipa tianschanica var. klemenzii steppe in China. Chin. J. Plant. Ecol. 41, 231–237 (2017).Article 

Google Scholar 
Qiao, X., Guo, K., Zhao, L., Yang, Y. & Zhao, H. Distribution, community characteristics and classification of Stipa basiplumosa steppe on Tibetan Plateau. Geogr. Res. 36, 2432–2440 (2017).
Google Scholar 
Qiao, X., Guo, K., Zhao, L., Wang, Z. & Liu, C. Community characteristics of Stipa bungeana alliance in China. Chin. J. Plant Ecol. 44, 986–994 (2020).Article 

Google Scholar 
Zhu, Y., Qiao, X., Guo, K., Xu, R. & Zhao, L. Distribution, community characteristics and classification of Stipa tianschanica var. gobica steppe in China. Chin. J. Plant Ecol. 42, 785–792 (2018).Article 

Google Scholar 
Li, X. R., Jia, X. H. & Dong, G. R. Influence of desertification on vegetation pattern variations in the cold semi-arid grasslands of Qinghai–Tibet plateau, north-west China. J. Arid. Environ. 64, 505–522 (2006).Article 

Google Scholar 
Tang, L. et al. Changes in vegetation composition and plant diversity with rangeland degradation in the alpine region of Qinghai–Tibet Plateau. Rangel. J. 37, 107–115 (2015).Article 

Google Scholar 
Zhou, X. Chinese Kobresia Meadow (Science Press, 2001).Wang, B. Z. et al. Potential distribution patterns of Stipa bungeana in China and the major factors influencing distribution. Acta Prataculturae Sinica 28, 3–13 (2019).
Google Scholar 
Sun, H., Li, W., Zhang, M. & Han, Y. A comprehensive scientific expedition to the Qinghai–Tibet Plateau. Resour. Sci. 8, 22–30 (1986).
Google Scholar 
Zhu, F. X. et al. Spatiotemporal variations of annual shallow soil temperature on the Tibetan Plateau during 1983–2013. Clim. Dyn. 51, 2209–2227 (2018).Article 

Google Scholar 
Chen, L. T. et al. Changes of carbon stocks in alpine grassland soils from 2002 to 2011 on the Tibetan Plateau and their climatic causes. Geoderma 288, 166–174 (2017).Article 

Google Scholar 
Ding, J. et al. Decadal soil carbon accumulation across Tibetan permafrost regions. Nat. Geosci. 10, 420–424 (2017).Article 

Google Scholar 
Tian, L. M. et al. Variations in soil nutrient availability across Tibetan grassland from the 1980s to 2010s. Geoderma 338, 197–205 (2019).Article 

Google Scholar 
Pepin, N. et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 5, 424–430 (2015).Article 

Google Scholar 
Chen, H. et al. The impacts of climate change and human activities on biogeochemical cycles on the Qinghai–Tibetan Plateau. Glob. Chang. Biol. 19, 2940–2955 (2013). This work suggests that warming enhanced NPP and soil respiration but many uncertainties remain.Article 

Google Scholar 
Shen, M. G. et al. Plant phenological responses to climate change on the Tibetan Plateau: research status and challenges. Natl Sci. Rev. 2, 454–467 (2015).Article 

Google Scholar 
Liu, X. D., Yin, Z. Y., Shao, X. M. & Qin, N. S. Temporal trends and variability of daily maximum and minimum, extreme temperature events, and growing season length over the eastern and central Tibetan Plateau during 1961–2003. J. Geophys. Res. Atmos. https://doi.org/10.1029/2005jd006915 (2006).Article 

Google Scholar 
Yang, K. et al. Recent climate changes over the Tibetan Plateau and their impacts on energy and water cycle: a review. Glob. Planet. Change 112, 79–91 (2014). This work reviews the main spatio-temporal characteristics of climate change on the QTP.Article 

Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Experimental warming, not grazing, decreases rangeland quality on the Tibetan Plateau. Ecol. Appl. 17, 541–557 (2007).Article 

Google Scholar 
Li, C. et al. Productivity and quality of alpine grassland vary with soil water availability under experimental warming. Front. Plant. Sci. 9, 1790 (2018).Article 

Google Scholar 
Peng, A. H. et al. Plant community responses to warming modified by soil moisture in the Tibetan Plateau. Arct. Antarct. Alp. Res. 52, 60–69 (2020).Article 

Google Scholar 
Li, F. et al. Leaf area rather than photosynthetic rate determines the response of ecosystem productivity to experimental warming in an alpine steppe. J. Geophys. Res. Biogeosci. 124, 2277–2287 (2019).Article 

Google Scholar 
Zong, N. et al. Responses of ecosystem CO2 fluxes to short-term experimental warming and nitrogen enrichment in an alpine meadow, northern Tibet Plateau. Sci. World J. 2013, 415318 (2013).Article 

Google Scholar 
Chen, Q., Niu, B., Hu, Y., Luo, T. & Zhang, G. Warming and increased precipitation indirectly affect the composition and turnover of labile-fraction soil organic matter by directly affecting vegetation and microorganisms. Sci. Total Environ. 714, 136787 (2020).Article 

Google Scholar 
Wang, X. X. et al. Effects of short-term and long-term warming on soil nutrients, microbial biomass and enzyme activities in an alpine meadow on the Qinghai–Tibet Plateau of China. Soil. Biol. Biochem. 76, 140–142 (2014).Article 

Google Scholar 
Jiang, L. L. et al. Plant organic N uptake maintains species dominance under long-term warming. Plant Soil 433, 243–255 (2018).Article 

Google Scholar 
Li, N. et al. Short-term effects of temperature enhancement on community structure and biomass of alpine meadow in the Qinghai–Tibet Plateau. Acta Ecol. Sin. 31, 895–905 (2011).
Google Scholar 
Jiang, Y., Fan, M. & Zhang, Y. Effect of short-term warming on plant community features of alpine meadow in northern Tibet. Chin. J. Ecol. 36, 616–622 (2017).
Google Scholar 
Wang, S. et al. Effects of warming and grazing on soil N availability, species composition, and ANPP in an alpine meadow. Ecology 93, 2365–2376 (2012). This work shows the effects of asymmetric warming and moderate grazing on plant composition, diversity, productivity and their relationships.Article 

Google Scholar 
Liu, P. et al. Ambient climate determines the directional trend of community stability under warming and grazing. Glob. Change Biol. 27, 5198–5210 (2021). This work finds that the negative effect of warming on plant diversity disappears with experimental duration, and ambient climate modulates the effects of warming and grazing on productivity stability.Article 

Google Scholar 
Zhang, B. et al. Responses of soil microbial communities to experimental warming in alpine grasslands on the Qinghai–Tibet Plateau. PLoS ONE 9, e103859 (2014).Article 

Google Scholar 
Chen, X. et al. Effects of warming and nitrogen fertilization on GHG flux in the permafrost region of an alpine meadow. Atmos. Environ. 157, 111–124 (2017).Article 

Google Scholar 
Zhang, Y. et al. Effects of grazing and climate warming on plant diversity, productivity and living state in the alpine rangelands and cultivated grasslands of the Qinghai–Tibetan Plateau. Rangel. J. 37, 57–65 (2015).Article 

Google Scholar 
Quan, Q. et al. High-level rather than low-level warming destabilizes plant community biomass production. J. Ecol. 109, 1607–1617 (2021).Article 

Google Scholar 
Wang, X. et al. Response of greenhouse gases emission fluxes to long-term warming in alpine meadow of northern Tibet. Chin. J. Agrometeorol. 39, 152–161 (2018).
Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Experimental warming causes large and rapid species loss, dampened by simulated grazing, on the Tibetan Plateau. Ecol. Lett. 7, 1170–1179 (2004).Article 

Google Scholar 
Zhang, C. H. et al. Recovery of plant species diversity during long-term experimental warming of a species-rich alpine meadow community on the Qinghai–Tibet Plateau. Biol. Conserv. 213, 218–224 (2017).Article 

Google Scholar 
Li, X. et al. Responses of biotic interactions of dominant and subordinate species to decadal warming and simulated rotational grazing in Tibetan alpine meadow. Sci. China Life Sci. 61, 849–859 (2018).Article 

Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Dynamic and complex microclimate responses to warming and grazing manipulations. Glob. Chang. Biol. 11, 1440–1451 (2005).Article 

Google Scholar 
Chen, J. et al. Warming effects on ecosystem carbon fluxes are modulated by plant functional types. Ecosystems 20, 515–526 (2017).Article 

Google Scholar 
Zhang, Y. Q. & Welker, J. M. Tibetan alpine tundra responses to simulated changes in climate: aboveground biomass and community responses. Arct. Alp. Res. 28, 203–209 (1996).Article 

Google Scholar 
Liu, H. et al. Shifting plant species composition in response to climate change stabilizes grassland primary production. Proc. Natl Acad. Sci. USA 115, 4051–4056 (2018). This work demonstrates that shifting plant species composition in response to climate change may have stabilized primary production in this high-elevation ecosystem, but also causes a shift from above-ground to below-ground productivity.Article 

Google Scholar 
Ganjurjav, H. et al. Differential response of alpine steppe and alpine meadow to climate warming in the central Qinghai–Tibetan Plateau. Agric. For. Meteorol. 223, 233–240 (2016).Article 

Google Scholar 
Jiang, F., Wei, X., Kang, B. & Shao, X. Effects of warming on alpine meadow diversity and primary productivity. Acta Agrestia Sin. 27, 298–305 (2019).
Google Scholar 
Zong, N. et al. Responses of plant community structure and species composition to warming and N addition in an alpine meadow, northern Tibetan Plateau, China. Chin. J. Appl. Ecol. 27, 3739–3748 (2016).
Google Scholar 
Peng, F. et al. Warming-induced shift towards forbs and grasses and its relation to the carbon sequestration in an alpine meadow. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/aa6508 (2017).Article 

Google Scholar 
Dorji, T. et al. Grazing and spring snow counteract the effects of warming on an alpine plant community in Tibet through effects on the dominant species. Agric. For. Meteorol. 263, 188–197 (2018).Article 

Google Scholar 
Xue, X., Peng, F., You, Q., Xu, M. & Dong, S. Belowground carbon responses to experimental warming regulated by soil moisture change in an alpine ecosystem of the Qinghai–Tibet Plateau. Ecol. Evol. 5, 4063–4078 (2015).Article 

Google Scholar 
Jing, X. et al. No temperature acclimation of soil extracellular enzymes to experimental warming in an alpine grassland ecosystem on the Tibetan Plateau. Biogeochemistry 117, 39–54 (2014).Article 

Google Scholar 
Yu, C. Q., Shen, Z. X., Zhang, X. Z., Sun, W. & Fu, G. Response of soil C and N, dissolved organic C and N, and inorganic N to short-term experimental warming in an alpine meadow on the Tibetan Plateau. Sci. World J. 2014, 152576 (2014).
Google Scholar 
Zhang, Y. et al. Simulated warming enhances the responses of microbial N transformations to reactive N input in a Tibetan alpine meadow. Environ. Int. 141, 105795 (2020).Article 

Google Scholar 
Jia, J. et al. Climate warming alters subsoil but not topsoil carbon dynamics in alpine grassland. Glob. Chang. Biol. 25, 4383–4393 (2019).Article 

Google Scholar 
Ding, X. L. et al. Warming increases microbial residue contribution to soil organic carbon in an alpine meadow. Soil. Biol. Biochem. 135, 13–19 (2019).Article 

Google Scholar 
Guan, S. et al. Climate warming impacts on soil organic carbon fractions and aggregate stability in a Tibetan alpine meadow. Soil. Biol. Biochem. 116, 224–236 (2018).Article 

Google Scholar 
Rui, Y. C. et al. Warming and grazing affect soil labile carbon and nitrogen pools differently in an alpine meadow of the Qinghai–Tibet Plateau in China. J. Soils Sediment. 11, 903–914 (2011).Article 

Google Scholar 
Heng, T., Wu, J., Xie, S. & Wu, M. The responses of soil C and N, microbial biomass C or N under alpine meadow of Qinghai–Tibet Plateau to the change of temperature and precipitation. Chin. Agric. Sci. Bull. 27, 425–430 (2011).
Google Scholar 
Li, N., Wang, G., Yang, Y., Gao, Y. & Liu, G. Plant production, and carbon and nitrogen source pools, are strongly intensified by experimental warming in alpine ecosystems in the Qinghai–Tibet Plateau. Soil. Biol. Biochem. 43, 942–953 (2011).Article 

Google Scholar 
Zhao, J. X. et al. Increased precipitation offsets the negative effect of warming on plant biomass and ecosystem respiration in a Tibetan alpine steppe. Agric. For. Meteorol. https://doi.org/10.1016/j.agrformet.2019.107761 (2019). This work shows that increased precipitation offsets the negative effect of warming on plant biomass and ecosystem respiration in a Tibetan alpine steppe.Article 

Google Scholar 
Wu, H. et al. Effects of increased precipitation combined with nitrogen addition and increased temperature on methane fluxes in alpine meadows of the Tibetan Plateau. Sci. Total Environ. 705, 135818 (2020).Article 

Google Scholar 
Shi, F. S., Chen, H., Chen, H. F., Wu, Y. & Wu, N. The combined effects of warming and drying suppress CO2 and N2O emission rates in an alpine meadow of the eastern Tibetan Plateau. Ecol. Res. 27, 725–733 (2012).Article 

Google Scholar 
Fu, G., Zhang, H. R. & Sun, W. Response of plant production to growing/non-growing season asymmetric warming in an alpine meadow of the northern Tibetan Plateau. Sci. Total Environ. 650, 2666–2673 (2019).Article 

Google Scholar 
Xiong, Q. L. et al. Warming and nitrogen deposition are interactive in shaping surface soil microbial communities near the alpine timberline zone on the eastern Qinghai–Tibet Plateau, southwestern China. Appl. Soil. Ecol. 101, 72–83 (2016).Article 

Google Scholar 
Wang, C. et al. Responses of plant leaf traits to simulated rainfall changes in alpine region. Acta Ecol. Sin. 41, 1–13 (2021).Article 

Google Scholar 
Zhang, K. et al. Effects of short-term warming and altered precipitation on soil microbial communities in alpine grassland of the Tibetan Plateau. Front. Microbiol. 7, 1032 (2016).
Google Scholar 
Evans, R. D. & Ehleringer, J. R. Water and nitrogen dynamics in an arid woodland. Oecologia 99, 233–242 (1994).Article 

Google Scholar 
Swap, R. J., Aranibar, J. N., Dowty, P. R., Gilhooly, W. P. III & Macko, S. A. Natural abundance of 13C and 15N in C3 and C4 vegetation of southern Africa: patterns and implications. Glob. Chang. Biol. 10, 350–358 (2004).Article 

Google Scholar 
Jia, Y. et al. Spatial and decadal variations in inorganic nitrogen wet deposition in China induced by human activity. Sci. Rep. 4, 3763–3763 (2014).Article 

Google Scholar 
Liu, Y. W., Xu, R., Wang, Y. S., Pan, Y. P. & Piao, S. L. Wet deposition of atmospheric inorganic nitrogen at five remote sites in the Tibetan Plateau. Atmos. Chem. Phys. 15, 11683–11700 (2015).Article 

Google Scholar 
Wang, W. et al. Atmospheric nitrogen deposition to a southeast Tibetan forest ecosystem. Atmosphere https://doi.org/10.3390/atmos11121331 (2020).Article 

Google Scholar 
Zou, X. et al. Ice-core based assessment of nitrogen deposition in the central Tibetan Plateau over the last millennium. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2021.152692 (2022).Article 

Google Scholar 
Aerts, R., Wallen, B. & Malmer, N. Growth-limiting nutrients in sphagnum-dominated bogs subject to low and high amospheric nitrogen supply. J. Ecol. 80, 131–140 (1992).Article 

Google Scholar 
Bai, Y. F. et al. Tradeoffs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning: evidence from Inner Mongolia grasslands. Glob. Chang. Biol. 16, 358–372 (2010).Article 

Google Scholar 
Du, Y. Population statistics of Qinghai–Tibet Plateau (1952–2016) (National Tibetan Plateau Data Center, 2019).Zhang, Y. J., Zhang, X. Q., Wang, X. Y., Liu, N. & Kan, H. M. Establishing the carrying capacity of the grasslands of China: a review. Rangel. J. 36, 1–9 (2014).Article 

Google Scholar 
Bardgett, R. D. et al. Combatting global grassland degradation. Nat. Rev. Earth Environ. 2, 720–735 (2021). This work shows that socio-ecological solutions are needed to combat degradation and promote restoration.Article 

Google Scholar 
Liu, M. et al. Effects of rotational and continuous overgrazing on newly assimilated C allocation. Biol. Fertil. Soils 57, 193–202 (2021).Article 

Google Scholar 
Yang, X. X. et al. Different responses of soil element contents and their stoichiometry (C:N:P) to yak grazing and Tibetan sheep grazing in an alpine grassland on the eastern Qinghai–Tibetan Plateau. Agric. Ecosyst. Environ. https://doi.org/10.1016/j.agee.2019.106628 (2019).Article 

Google Scholar 
Lin, B., Zhao, X. R., Zheng, Y., Qi, S. & Liu, X. Z. Effect of grazing intensity on protozoan community, microbial biomass, and enzyme activity in an alpine meadow on the Tibetan Plateau. J. Soils Sediment. 17, 2752–2762 (2017).Article 

Google Scholar 
Ma, W. M., Ding, K. Y. & Li, Z. W. Comparison of soil carbon and nitrogen stocks at grazing-excluded and yak grazed alpine meadow sites in Qinghai–Tibetan Plateau, China. Ecol. Eng. 87, 203–211 (2016).Article 

Google Scholar 
Li, W. et al. Effects of grazing regime on vegetation structure, productivity, soil quality, carbon and nitrogen storage of alpine meadow on the Qinghai–Tibetan Plateau. Ecol. Eng. 98, 123–133 (2017).Article 

Google Scholar 
Sun, J. et al. Effects of grazing regimes on plant traits and soil nutrients in an alpine steppe, northern Tibetan Plateau. PLoS ONE 9, e108821 (2014).Article 

Google Scholar 
Niu, K. C., He, J. S. & Lechowicz, M. J. Grazing-induced shifts in community functional composition and soil nutrient availability in Tibetan alpine meadows. J. Appl. Ecol. 53, 1554–1564 (2016).Article 

Google Scholar 
Luan, J. W. et al. Different grazing removal exclosures effects on soil C stocks among alpine ecosystems in east Qinghai–Tibet Plateau. Ecol. Eng. 64, 262–268 (2014).Article 

Google Scholar 
Wei, D. et al. Responses of CO2, CH4 and N2O fluxes to livestock exclosure in an alpine steppe on the Tibetan Plateau, China. Plant Soil 359, 45–55 (2012).Article 

Google Scholar 
Shen, H. et al. Grazing enhances plant photosynthetic capacity by altering soil nitrogen in alpine grasslands on the Qinghai–Tibetan Plateau. Agric. Ecosyst. Environ. 280, 161–168 (2019).Article 

Google Scholar 
Jiang, L. et al. Grazing modifies inorganic and organic nitrogen uptake by coexisting plant species in alpine grassland. Biol. Fertil. Soils 52, 211–221 (2016).Article 

Google Scholar 
Sun, Y., Schleuss, P. M., Pausch, J., Xu, X. L. & Kuzyakov, Y. Nitrogen pools and cycles in Tibetan Kobresia pastures depending on grazing. Biol. Fertil. Soils 54, 569–581 (2018).Article 

Google Scholar 
Chen, B. et al. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai–Tibet Plateau. Agric. For. Meteorol. 189-190, 11–18 (2014).Article 

Google Scholar 
Wang, Z. Q. et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai–Tibet Plateau, in China. Ecol. Inform. 33, 32–44 (2016).Article 

Google Scholar 
Huang, K. et al. The influences of climate change and human activities on vegetation dynamics in the Qinghai–Tibet Plateau. Remote Sens. 8, 876 (2016).Article 

Google Scholar 
Li, L. et al. Increasing sensitivity of alpine grasslands to climate variability along an elevational gradient on the Qinghai–Tibet Plateau. Sci. Total Environ. 678, 21–29 (2019).Article 

Google Scholar 
Wang, Z. et al. Vegetation expansion on the Tibetan Plateau and its relationship with climate change. Remote. Sens. https://doi.org/10.3390/rs12244150 (2020).Article 

Google Scholar 
Wu, J. et al. Disentangling climatic and anthropogenic contributions to nonlinear dynamics of alpine grassland productivity on the Qinghai–Tibetan Plateau. J. Environ. Manag. 281, 111875 (2021).Article 

Google Scholar 
Fu, G., Shen, Z. X. & Zhang, X. Z. Increased precipitation has stronger effects on plant production of an alpine meadow than does experimental warming in the northern Tibetan Plateau. Agric. For. Meteorol. 249, 11–21 (2018).Article 

Google Scholar 
Hu, Y. et al. Effect of increasing precipitation and warming on microbial community in Tibetan alpine steppe. Environ. Res. 189, 109917 (2020).Article 

Google Scholar 
Ma, Z. et al. Climate warming reduces the temporal stability of plant community biomass production. Nat. Commun. 8, 15378 (2017).Article 

Google Scholar 
Bai, W., Xi, J. & Wang, G. Effects of short-term warming and nitrogen addition on CO2 emission during growing season in an alpine swamp meadow ecosystem of Qinghai–Tibetan Plateau. Chin. J. Ecol. 38, 927–936 (2019).
Google Scholar 
Bai, W., Wang, G. X., Xi, J. Y., Liu, Y. W. & Yin, P. S. Short-term responses of ecosystem respiration to warming and nitrogen addition in an alpine swamp meadow. Eur. J. Soil Biol. 92, 16–23 (2019).Article 

Google Scholar 
Ge, Y. et al. Effects of warming and nitrogen addition on plant community structure and species diversity of alpine meadow in northern Tibet. Ecol. Environ. Sci. 28, 2185–2191 (2019).
Google Scholar 
Zong, N. et al. Effects of warming and nitrogen addition on community production and biomass allocation in an alpine meadow. Chin. J. Appl. Ecol. 29, 59–67 (2018).
Google Scholar 
Zhu, X. X. et al. Effects of warming, grazing/cutting and nitrogen fertilization on greenhouse gas fluxes during growing seasons in an alpine meadow on the Tibetan Plateau. Agric. For. Meteorol. 214, 506–514 (2015).Article 

Google Scholar 
Fu, G. et al. Clipping alters the response of biomass production to experimental warming: a case study in an alpine meadow on the Tibetan Plateau, China. J. Mt. Sci. 12, 935–942 (2015).Article 

Google Scholar 
Chen, S. P. et al. Plant diversity enhances productivity and soil carbon storage. Proc. Natl Acad. Sci. USA. 115, 4027–4032 (2018).Article 

Google Scholar 
Wu, J. S. et al. Effects of livestock exclusion and climate change on aboveground biomass accumulation in alpine pastures across the northern Tibetan Plateau. Chin. Sci. Bull. 59, 4332–4340 (2014).Article 

Google Scholar 
Sun, J., Cheng, G. W., Li, W. P., Sha, Y. K. & Yang, Y. C. On the variation of NDVI with the principal climatic elements in the Tibetan Plateau. Remote Sens. 5, 1894–1911 (2013).Article 

Google Scholar 
Sun, J. et al. Reconsidering the efficiency of grazing exclusion using fences on the Tibetan Plateau. Sci. Bull. 65, 1405–1414 (2020). This work finds that fencing enclosures lead to some negative impacts, such as hindering wildlife movement.Article 

Google Scholar 
Yu, C. et al. Grazing exclusion to recover degraded alpine pastures needs scientific assessments across the northern Tibetan Plateau. Sustainability https://doi.org/10.3390/su8111162 (2016).Article 

Google Scholar 
Wu, J. & Wang, X. Effect of enclosure ages on community characters and biomas of the degraded alpine steppe at the northern Tibet. Acta Agrestia Sin. 25, 261–266 (2017).
Google Scholar 
Zhao, J. X., Luo, T. X., Li, R. C., Li, X. & Tian, L. H. Grazing effect on growing season ecosystem respiration and its temperature sensitivity in alpine grasslands along a large altitudinal gradient on the central Tibetan Plateau. Agric. For. Meteorol. 218, 114–121 (2016).Article 

Google Scholar 
Deng, L., Sweeney, S. & Shangguan, Z. P. Grassland responses to grazing disturbance: plant diversity changes with grazing intensity in a desert steppe. Grass Forage Sci. 69, 524–533 (2014).Article 

Google Scholar 
Yuan, Z., Epstein, H. & Li, G. Grazing exclusion did not affect soil properties in alpine meadows in the Tibetan permafrost region. Ecol. Eng. https://doi.org/10.1016/j.ecoleng.2019.105657 (2020).Article 

Google Scholar 
Zhang, W. et al. Effects of banning grazing and delaying grazing on species diversity and biomass of alpine meadow in northern Tibet. J. Agric. Sci. Technol. 15, 143–149 (2013).
Google Scholar 
Miao, F., Guo, Y., Miao, P., Guo, Z. & Shen, Y. The northeast edge of Qinghai–Tibet Plateau area of alpine meadow community characteristics respond to nurture. Acta Prataculture Sin. 21, 11–16 (2012).
Google Scholar 
Lu, X. et al. Short-term grazing exclusion has no impact on soil properties and nutrients of degraded alpine grassland in Tibet, China. Solid Earth 6, 1195–1205 (2015).Article 

Google Scholar 
Gao, Y. H., Zeng, X. Y., Schumann, M. & Chen, H. Effectiveness of exclosures on restoration of degraded alpine meadow in the eastern Tibetan Plateau. Arid. Land. Res. Manag. 25, 164–175 (2011).Article 

Google Scholar 
Yao, X. X. et al. Effects of long term fencing on biomass, coverage, density, biodiversity and nutritional values of vegetation community in an alpine meadow of the Qinghai–Tibet Plateau. Ecol. Eng. 130, 80–93 (2019).Article 

Google Scholar 
Chen, W., Chang, H. & Liu, R. Fractal features of soil particle size distributions and their implications for indicating enclosure management in a semiarid grassland ecosystem. Pol. J. Ecol. 68, 132–144 (2020).
Google Scholar 
Smith, D., King, R. & Allen, B. L. Impacts of exclusion fencing on target and non-target fauna: a global review. Biol. Rev. 95, 1590–1606 (2020).Article 

Google Scholar 
Zhang, Y. et al. Assessment of effectiveness of nature reserves on the Tibetan Plateau based on net primary production and the large sample comparison method. J. Geogr. Sci. 26, 27–44 (2016).Article 

Google Scholar 
Hu, J. Research on the status quo and problems of natural reserve construction in Qinghai–Tibet Plateau. Environ. Dev. 32, 204–206 (2020).
Google Scholar 
Shao, Q., Fan, J., Liu, J., Cao, W. & Liu, L. 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. 32, 35–44 (2017).
Google Scholar 
Liu, F. & Zeng, Y. N. Spatial–temporal change in vegetation net primary productivity and its response to climate and human activities in Qinghai Plateau in the past 16 years. Acta Ecol. Sin. 39, 1528–1540 (2019).
Google Scholar 
Zhang, Y., Wu, X., Qi, W., Li, S. & Bai, W. Characteristics and protection effectiveness of nature reserves on the Tibetan Plateau, China. Resources. Science 37, 1455–1464 (2015).
Google Scholar 
Buckley, M. C. & Crone, E. E. Negative off-site impacts of ecological restoration: understanding and addressing the conflict. Conserv. Biol. 22, 1118–1124 (2008).Article 

Google Scholar 
Cao, S. X. & Zhang, J. Political risks arising from the impacts of large-scale afforestation on water resources of the Tibetan Plateau. Gondwana Res. 28, 898–903 (2015).Article 

Google Scholar 
Li, Y. & Li, W. Why “Balance of Forage and Livestock” system failed to reach sustainable grassland utilization. China Agric. Univ. J. Soc. Sci. Ed. 29, 124–131 (2012).
Google Scholar 
Du, S. A Study on the Satisfaction Degree of Herdsmen’s Income and Grassland Ecological Compensation Policy. Master thesis, Lanzhou Univ. (2019).Yu, H., Wang, G., Yang, Y. & Lü, Y. Concept of grassland green carrying capacity and its application framework in national park. Acta Ecol. Sin. 40, 7248–7254 (2020).
Google Scholar 
Deng, Y. & Li, C. The investigation and research about the Farmland Retirement and Environment Project in the Yangtze River headwaters area. Ecol. Econ. 2, 77–80 (2006).
Google Scholar 
Zhou, Q. et al. Analysis on the relationship between grassland area and forage-livestock balance in Qinghai–Tibet Plateau. Chin. J. Grassl. 41, 110–117 (2019).
Google Scholar 
Li, Y. et al. Awareness and reaction of herdsmen to the policy of returning grazing land to grasslands in the Changtang Plateau,Tibet. Pratacultural Sci. 30, 788–794 (2013).
Google Scholar 
Fan, J. et al. Third pole national park group construction is scientific choice for implementing strategy of major function zoning and green development in Tibet, China. Bull. Chin. Acad. Sci. 32, 932–944 (2017).
Google Scholar 
Xu, Z., Cheng, S. & Gao, L. Impacts of herders sedentarization on regional spatial heterogeneity and grassland ecosystem change in pastoral area. J. Arid. Land. Resour. Environ. 31, 8–13 (2017).
Google Scholar 
Ptackova, J. Sedentarisation of Tibetan nomads in China: implementation of the Nomadic settlement project in the Tibetan Amdo area; Qinghai and Sichuan Provinces. Pastoralism https://doi.org/10.1186/2041-7136-1-4 (2011).Article 

Google Scholar 
Weber, K. T. & Horst, S. Desertification and livestock grazing: the roles of sedentarization, mobility and rest. Pastoralism https://doi.org/10.1186/2041-7136-1-19 (2011).Article 

Google Scholar 
Zhang, J. et al. Ecological consequence of nomad settlement policy in the pasture area of Qinghai–Tibetan Plateau: from plant and soil perspectives. J. Environ. Manage. https://doi.org/10.1016/j.jenvman.2020.110114 (2020).Article 

Google Scholar 
Li, C. X., de Jong, R., Schmid, B., Wulf, H. & Schaepman, M. E. Spatial variation of human influences on grassland biomass on the Qinghai–Tibetan Plateau. Sci. Total Environ. 665, 678–689 (2019).Article 

Google Scholar 
Kuang, W. Dataset of Urban Distribution, Urban Population and Built-up Area in Tibetan Plateau (2000–2015) (National Tibetan Plateau Data Center, 2021).Tian, L. & Chen, J. Urban expansion inferenced by ecosystem production on the Qinghai–Tibet plateau. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/ac3178 (2022).Article 

Google Scholar 
Liu, Y. & Lu, C. Quantifying grass coverage trends to identify the hot plots of grassland degradation in the Tibetan Plateau during 2000–2019. Int. J. Environ. Res. Public Health https://doi.org/10.3390/ijerph18020416 (2021).Article 

Google Scholar 
Tang, L. et al. Warming counteracts grazing effects on the functional structure of the soil microbial community in a Tibetan grassland. Soil. Biol. Biochem. 134, 113–121 (2019).Article 

Google Scholar 
Zhong, L. Tourism Planning Case in Tibetan Plateau (China Tourism Press, 2018).La, M. Discussion of the coordinated development of tourism development and ecological Environment Protection in Tibetan. Soc. Sci. Res. 6, 118–120, (2013).
Google Scholar 
Zhuang, M. et al. Opportunities for household energy on the Qinghai–Tibet Plateau in line with United Nations’ Sustainable Development Goals. Renew. Sustain. Energy Rev. https://doi.org/10.1016/j.rser.2021.110982 (2021).Article 

Google Scholar 
Ruess, R. W. & Mcnaughton, S. J. Grazing and the dynamics of nutrient and energy regulated microbial processes in the serengeti grasslands. Oikos 49, 101–110 (1987).Article 

Google Scholar 
Li, M. et al. Changes in plant species richness distribution in Tibetan alpine grasslands under different precipitation scenarios. Glob. Ecol. Conserv. 21, 13 (2020).
Google Scholar 
Wang, Z. et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai–Tibet Plateau, in China. Ecol. Inform. 33, 32–44 (2016).Article 

Google Scholar 
Muñoz Sabater, J. ERA5-Land monthly averaged data from 1981 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS), https://doi.org/10.24381/cds.68d2bb30 (2019).Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. Terraclimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191 (2018).Article 

Google Scholar  More

Zaneveld, J. R., McMinds, R. & Thurber, R. V. Stress and stability: Applying the Anna Karenina principle to animal microbiomes. Nat. Microbiol. 2, 1–8 (2017).Article 
CAS 

Google Scholar 
Trevelline, B. K., Fontaine, S. S., Hartup, B. K. & Kohl, K. D. Conservation biology needs a microbial renaissance: A call for the consideration of host-associated microbiota in wildlife management practices. Proc. R. Soc. B 286, 2018–2448 (2019).Article 

Google Scholar 
Bennett, A. G. On the occurrence of diatoms on the skin of whales. Proc. R. Soc. Lond. B 91, 352–357 (1920).ADS 
Article 

Google Scholar 
Denys, L. Morphology and taxonomy of epizoic diatoms (Epiphalaina and Tursiocola) on a sperm whale (Physeter macrocephalus) stranded on the coast of Belgium. Diatom. Res. 12, 1–18 (1997).Article 

Google Scholar 
Majewska, R. Tursiocola neliana sp. nov (Bacillariophyceae) epizoic on South African leatherback sea turtles (Dermochelys coriacea) and new observations on the genus Tursiocola. Phytotaxa 453, 1–15 (2020).Article 

Google Scholar 
Majewska, R. et al. Chelonicola and Poulinea, two new gomphonemoid genera living on marine turtles from Costa Rica. Phytotaxa 233, 236–250 (2015).Article 

Google Scholar 
Majewska, R. et al. Shared epizoic taxa and differences in diatom community structure between green turtles (Chelonia mydas) from distant habitats. Microb Ecol. 74, 969–978 (2017).PubMed 
Article 

Google Scholar 
Majewska, R. et al. Two new epizoic Achnanthes species (Bacillariophyta) living on marine turtles from Costa Rica. Bot. Mar. 60, 303–318 (2017).Article 

Google Scholar 
Majewska, R., De Stefano, M. & Van de Vijver, B. Labellicula lecohuiana, a new epizoic diatom species living on green turtles in Costa Rica. Nova Hedwig Beih. 146, 23–31 (2018).Article 

Google Scholar 
Majewska, R. et al. Craspedostauros alatus sp. nov., a new diatom (Bacillariophyta) species found on museum sea turtle specimens. Diatom Res. 33, 229–240 (2018).Article 

Google Scholar 
Majewska, R. et al. Six new epibiotic Proschkinia (Bacillariophyta) species and new insights into the genus phylogeny. Eur. J. Phycol. 54, 609–631 (2019).Article 

Google Scholar 
Majewska, R., Robert, K., Van de Vijver, B. & Nel, R. A new species of Lucanicum (Cyclophorales, Bacillariophyta) associated with loggerhead sea turtles from South Africa. Bot. Lett. 167, 7–14 (2020).Article 

Google Scholar 
Frankovich, T. A., Sullivan, M. J. & Stacy, N. I. Tursiocola denysii sp. Nov. (Bacillariophyta) from the neck skin of Loggerhead sea turtles (Caretta caretta). Phytotaxa 234, 227–236 (2015).Article 

Google Scholar 
Frankovich, T. A., Ashworth, M. P., Sullivan, M. J., Vesela, J. & Stacy, N. I. Medlinella amphoroidea gen. et sp. Nov. (Bacillariophyta) from the neck skin of Loggerhead sea turtles (Caretta caretta). Phytotaxa 272, 101–114 (2016).Article 

Google Scholar 
Riaux-Gobin, C. et al. New epizoic diatom (Bacillariophyta) species from sea turtles in the Eastern Caribbean and South Pacific. Diatom Res. 32, 109–125 (2017).Article 

Google Scholar 
Riaux-Gobin, C., Witkowski, A., Chevallier, D. & Daniszewska-Kowalczyk, G. Two new Tursiocola species (Bacillariophyta) epizoic on green turtles (Chelonia mydas) in French Guiana and Eastern Caribbean. Fottea Olomouc 17, 150–163 (2017).Article 

Google Scholar 
Riaux-Gobin, C., Witkowski, A., Kociolek, J. P. & Chevallier, D. Navicula dermochelycola sp. Nov., presumably an exclusively epizoic diatom on sea turtles Dermochelys coriacea and Lepidochelys olivacea from French Guiana. Oceanol. Hydrobiol. Stud. 49, 132–139 (2020).CAS 
Article 

Google Scholar 
Robert, K., Bosak, S. & Van de Vijver, B. Catenula exigua sp. nov., a new marine diatom (Bacillariophyta) species from the Adriatic Sea. Phytotaxa 414, 113–118 (2019).Article 

Google Scholar 
Van de Vijver, B. & Bosak, S. Planothidium kaetherobertianum, a new marine diatom (Bacillariophyta) species from the Adriatic Sea. Phytotaxa 425, 105–112 (2019).Article 

Google Scholar 
Robinson, N. J. et al. Epibiotic diatoms are universally present on all sea turtle species. PLoS ONE 11, e0157011 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Van de Vijver, B. et al. Diversity of diatom communities (Bacillariophyta) associated with loggerhead sea turtles. PLoS ONE 15, e0236513 (2020).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Van de Vijver, B., Robert, K., Witkowski, A. & Bosak, S. Majewskaea gen. nov. (Bacillariophyta), a new marine benthic diatom genus from the Adriatic Sea. Fottea 20, 112–120 (2020).Article 

Google Scholar 
Majewska, R. Nagumoea hydrophicola sp. Nov. (Bacillariophyta), the first diatom species described from sea snakes. Diatom Res. 36, 49–59 (2021).Article 

Google Scholar 
Frankovich, T. A., Sullivan, M. J. & Stacey, N. I. Three new species of Tursiocola (Bacillariophyta) from the skin of the West Indian manatee (Trichechus manatus). Phytotaxa 204, 33–48 (2015).Article 

Google Scholar 
Frankovich, T. A., Ashworth, M. P., Sullivan, M. J., Theriot, E. C. & Stacy, N. I. Epizoic and apochlorotic Tursiocola species (Bacillariophyta) from the skin of Florida manatees (Trichechus manatus latirostris). Protist 169, 539–568 (2018).PubMed 
Article 

Google Scholar 
Azari, M. et al. Diatoms on sea turtles and floating debris in the Persian Gulf (Western Asia). Phycologia 59, 292–304 (2020).Article 

Google Scholar 
Majewska, R. & Goosen, W. E. For better, for worse: Manatee-associated Tursiocola (Bacillariophyta) remain faithful to their host. J. Phycol. 56, 1019–1027 (2020).CAS 
PubMed 
Article 

Google Scholar 
Smol, J. P. & Stoermer, E. F. The Diatoms: Applications for the Environmental and Earth Sciences (Cambridge University Press, 2010).Book 

Google Scholar 
Rivera, S. F. et al. DNA metabarcoding and microscopic analyses of sea turtles biofilms: Complementary to understand turtle behavior. PLoS ONE 13, e0195770 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Majewska, R. et al. On sea turtle-associated Craspedostauros with description of three novel species. J Phycol. 57, 199–208 (2021).CAS 
PubMed 
Article 

Google Scholar 
Holmes, R. W. The morphology of diatoms epizoic on cetaceans and their transfer from Cocconeis to two new genera, Bennettella and Epipellis. Br. Phycol. J. 20, 43–57 (1985).Article 

Google Scholar 
Woodworth, K. A., Frankovich, T. A. & Freshwater, D. W. Melanothamnus maniticola (Ceramiales, Rhodophyta): An epizoic species evolved for life on the West Indian Manatee. J. Phycol. 55, 1239–1245 (2019).CAS 
PubMed 
Article 

Google Scholar 
Vitt, L. J. & Caldwell, J. P. Herpetology: An Introductory Biology of Amphibians and Reptiles (Academic Press, 2013).
Google Scholar 
Pitman, L. R. et al. Skin in the game: Epidermal molt as a driver of long-distance migration in whales. Mar. Mamm. Sci. 36, 565–594 (2020).Article 

Google Scholar 
Pope, D. H. & Berger, L. R. Algal photosynthesis at increased hydrostatic pressure and constant pO2. Arch. Microbiol. 89, 321–325 (1973).CAS 

Google Scholar 
Calcagno, V., Jarne, P., Loreau, M., Mouquet, N. & David, P. Diversity spurs diversification in ecological communities. Nat. Commun. 8, 15810 (2017).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Robinson, N. J. & Pfaller, J. B. Sea turtle epibiosis: Global patterns and knowledge gaps. Trends Evol. Ecol. 10, 844021 (2021).
Google Scholar 
Conant, T. A., Dutton, P. H., Eguchi, T., Epperly, S. P., Fahy, C. C., Godfrey, M. H., MacPherson, S. L., Possardt, E. E., Schroeder, B. A., Seminoff, J. A., Snover, M. L. Loggerhead sea turtle (Caretta caretta) 2009 status review under the US Endangered Species Act. In Report of the loggerhead biological review Team to the National Marine Fisheries Service. 222, 1–230 (2009).Evans, K. M., Wortley, A. H. & Mann, D. G. An assessment of potential diatom ‘“barcode”’ genes (cox1, rbcL, 18S and ITS rDNA) and their effectiveness in determining relationships in Sellaphora (Bacillariophyta). Protist 158, 349–364 (2007).CAS 
PubMed 
Article 

Google Scholar 
Hamsher, S. E., Evans, K. M., Mann, D. G., Poulíčková, A. & Saunders, G. W. Barcoding diatoms: Exploring alternatives to COI-5P. Protist 162, 405–422 (2011).CAS 
PubMed 
Article 

Google Scholar 
Bowen, B. W. & Karl, S. A. Population genetics and phylogeography of sea turtles. Mol Ecol. 16, 4886–4907 (2007).CAS 
PubMed 
Article 

Google Scholar 
Shanker, K., Ramadevi, J., Choudhury, B. C., Singh, L. & Aggarwal, R. K. Phylogeography of olive ridley turtles (Lepidochelys olivacea) on the east coast of India: implications for conservation theory. Mol. Ecol. 13, 1899–1909 (2004).CAS 
PubMed 
Article 

Google Scholar 
Pinou, T. et al. Standardizing sea turtle epibiont sampling: Outcomes of the epibiont workshop at the 37th International Sea Turtle Symposium. Mar. Turt. Newsl. 157, 22–32 (2019).
Google Scholar 
Ehrhert L., Ogren L. H. Studies in foraging habitats: capturing and handling turtles. In Research and management techniques for the conservation of sea turtles (eds. Eckert, K. L., Bjorndal, K. A., Abreu-Grobois, F. A., Donnelly, M.). IUCN/SSC Marine Turtle Specialist Group. Publication No. 4. (1999).Guillard, R. R. Culture of phytoplankton for feeding marine invertebrates. In Culture of Marine Invertebrate Animals 29–60 (Springer, 1975).Theriot, E. C., Ashworth, M. P., Nakov, T., Ruck, E. & Jansen, R. K. Dissecting signal and noise in diatom chloroplast protein encoding genes with phylogenetic information profiling. Mol. Phylogenet. Evol. 89, 28–36 (2015).CAS 
PubMed 
Article 

Google Scholar 
Lobban, C. S., Ashworth, M. P., Calaor, J. J. & Theriot, E. C. Extreme diversity in fine-grained morphology reveals fourteen new species of conopeate Nitzschia (Bacillariophyta: Bacillariales). Phytotaxa. 401, 199–238 (2019).Article 

Google Scholar 
Lanfear, R., Frandsen, P. B., Wright, A. M., Senfeld, T. & Calcott, B. PartitionFinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 34, 772–773 (2017).CAS 
PubMed 

Google Scholar 
Lanfear, R., Calcott, B., Kainer, D., Mayer, C. & Stamatakis, A. Selecting optimal partitioning schemes for phylogenomic datasets. BMC Evol. Biol. 14, 1–14 (2014).Article 

Google Scholar 
Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nguyen, L.-T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
PubMed 
Article 

Google Scholar 
Chernomor, O., Von Haeseler, A. & Minh, B. Q. Terrace aware data structure for phylogenomic inference from supermatrices. Syst. Biol. 65, 997–1008 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Aberer, A. J., Kobert, K. & Stamatakis, A. ExaBayes: Massively parallel bayesian tree inference for the whole-genome Era. Mol. Biol. Evol. 31, 2553–2556 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Sample collection and metagenomic sequencingWritten informed consent was obtained prior to participation in the project. The study protocol for the Japanese (Disease, Drug, Diet, Daily life) microbiome project was approved by the medical ethics committees of the Tokyo Medical University (Approval No: T2019-0119), National Center for Global Health and Medicine (Approval No: 1690), the University of Tokyo (Approval No: 2019185NI), Waseda University (Approval No: 2018-318), and the RIKEN Center for Integrative Medical Sciences (Approval No: H30-7). We conducted a prospective cross-sectional study of 4198 individuals participating in the Japanese 4D microbiome project, which commenced in January 2015 and is ongoing20.Participants registered in the project were those who visited hospitals in the area for disease diagnosis or a health checkup. Faecal samples are collected from both healthy and diseased participants. The eligibility criteria for participants are as follows: (1) born and raised in Japan; (2) age >15 years; (3) written informed consent provided; and (4) having an endoscopic diagnosis on colonoscopy; either having undergone a colonoscopy within the last 3 years or planning to undergo colonoscopy for colorectal cancer screening, surveillance, and diagnosis of various gastrointestinal symptoms. The exclusion criteria were as follows: (1) suspected acute infectious disease based on clinical findings (e.g., acute enterocolitis, pneumonia, tuberculosis etc.); (2) acute bleeding; (3) hearing loss; (4) unable to understand written documents; (5) unable to write and (6) limited ability to perform activities of daily living. No compensation was paid to participants.Participants collected faecal samples using a Cary–Blair medium-containing tube60 at home, and the samples were refrigerated for up to 2 days before the hospital visit. Immediately after participants arrived at the hospital, their faecal samples were frozen at −80 °C until DNA extraction. We avoided collecting samples within 1 month of administering bowel preparation for colonoscopy because it has a profound effect on the gut microbiome and metabolome61. Health professionals checked that the amount of stool was sufficient for analysis. Shotgun metagenomic sequencing was performed for 4241 faecal samples and quality controls were conducted20, from which 43 samples were excluded from further analyses due to the low number of high-quality reads (130 bp. Encoded genes in the contigs were predicted by MetaGeneMark (3.38)70. Assembled contigs were defined as phages if they passed all of the following six criteria.

1.

A genome size threshold was applied, and contigs less than 10 Kb were excluded, as typical dsDNA phages have genomes larger than >10 Kb71.

2.

Viral-specific k-mer patterns were checked by DeepVirFinder (v1.0)22. Contigs with p-values >0.05 were excluded from further analysis.

3.

To detect viral hallmark genes (VHGs) and plasmid hallmark genes, we performed a highly sensitive HMM-HMM search against the Pfam database72. First, the encoded genes were aligned to the viral protein database, collected from complete (circular) viral genomes (n = 13,628) in the IMG/VR v2 database30 using JackHMMER. The obtained HMM profiles were searched against the Pfam database using hhblits73 with a  >95% probability cut-off. These procedures were performed using the pipeline_for_high_sensitive_domain_search script (https://github.com/yosuken/pipeline_for_high_sensitive_domain_search)74,75. Contigs with plasmid hallmark genes or those without VHGs were excluded. The hallmark genes used in this analysis are summarised in Supplementary Data 3.

4.

The presence of housekeeping marker genes of prokaryotic species was checked by fetchMG (v1.0)76, and ribosomal RNA genes (5 S, 16 S and 23 S) were identified by barrnap (0.9) (https://github.com/tseemann/barrnap). Contigs with the marker genes and ribosomal RNA genes were excluded from further analysis.

5.

The encoded genes of each contig were aligned to the viral protein database and a plasmid protein database constructed from the reference plasmids in RefSeq (n = 16,136, in April 2020) using DIAMOND (v0.9.29.130)77 with the more-sensitive option. The number of genes aligned to each database was compared, and contigs with more genes aligned to the plasmid protein database were excluded from further analysis.

6.

The proportion of provirus regions was assessed by CheckV (v0.7)24, and contigs estimated with 70% and 10% contamination.To evaluate the performance of this custom pipeline, we applied the pipeline to reference phage genomes (n = 2609, as positive data) and plasmid sequences (n = 16,136, as negative data) in Refseq. The true positive rate was defined as the number of phages detected as phages by the pipeline divided by the number of reference phages. The false positive rate was defined as the number of plasmids detected as phages by the pipeline divided by the number of reference plasmids. DeepVirFinder22, VirSorter (v1.0.3)23 Virsorter2 (2.2.3)25, VIBRANT (v1.2.1)26, Seeker (v1.0.3)27 and ViralVerify (v1.1)28 were also applied to the same datasets with the default parameters, and the performance was compared among them.Analysis of phage genomesViral operational taxonomic units (vOTUs) were constructed by clustering phage genomes with a  > 95% identity29 using dRep (v2.2.3)78 with the default options. Representative sequences of each vOTU selected by dRep were further clustered with reference sequences in RefSeq, IMG/VR30, gut virome database (GVD)15, gut phage database (GPD)9, and metagenomic gut virus (MGV) database31 with >95% identity and >85% length coverage using aniclust.py script in the CheckV package to identify common sequences among the databases.To further construct broader viral clusters (VC), proportions of protein clusters shared between phages were assessed. First, to define protein clusters, similarity searches of all protein sequences from all the phages identified in this study were performed using DIAMOND with the more-sensitive option (e-value 20% of clusters were grouped as a VC, which corresponds approximately to family- or subfamily-level clusters7,37. Rarefaction curves of the vOTUs and VCs were estimated with the iNEXT function in the iNEXT package (v2.0.20)80. The similarity matrix of the phages based on the percentage of shared protein clusters was further projected by tSNE using the tsne function in the Rtsne package (v0.16).Taxonomy annotation of phages was performed with a voting approach described previously16 with minor modifications. First, the protein sequences of each phage were aligned to viral proteins detected from phage genomes in RefSeq (n = 2609, in April 2020) using DIAMOND with the more-sensitive option. Then, the best-hit taxonomy of each protein (family levels) was counted, and the most common taxonomy was assigned to the phage if >20% of proteins in the phage were aligned to the same taxonomy.Phage lifestyles (i.e. virulent or temperate) were predicted by BACPHLIP40 and alignments to reference bacterial genomes in the RefSeq. Phages were defined as temperate if the BACPHLIP score was >0.8 or the phage genome was aligned to any reference genomes with >1000 bp alignment length with >95% identity.Host predictionBacterial and archaeal genomes were downloaded from the RefSeq database (in April 2019). To reduce the redundancy of genomes from closely related strains in the same species (e.g. Escherichia coli), 10 genomes were selected randomly for species with more than 10 genomes, and other genomes were excluded from the dataset. The reference dataset consisted of 33,215 bacterial and 822 archaeal genomes.Host prediction of the identified phages was performed using CRISPR spacers81. CRISPR spacers were predicted from the reference microbial genomes and assembled contigs ( >10,000 bp) from the 4198 metagenomic datasets using PILER-CR (1.06)82. Short (100 bp) spacers were discarded. In total, 679,323 and 283,619 spacers were identified from the reference microbial genomes and assembled contigs, respectively. Taxonomy information was assigned to the assembled contigs if they were aligned to the microbial reference genomes with >90% identity and >70% length coverage thresholds using MiniMap283. The CRISPR spacers were mapped to the phage genomes using BLASTN with the option for short sequences: -a20 -m9 -e1 -G10 -E2 -q1 -W7 -F F81. CRISPR spacers, which were mapped with 100% identity or 1 mismatch/indel with >95% sequence alignment, were used for host assignment at the genus level. Assignments of host species were checked manually, and if any of the following non-human intestinal species were assigned, the host was excluded: Dickeya, Anaerobutyricum, Rubellimicrobium, Eisenbergiella, Harryflintia, Leucothrix, Photorhabdus, Spirosoma, Syntrophobotulus, Thermincola, Algoriphagus, Franconibacter, Kandleria, Lawsonibacter, Methylomonas, Provencibacterium, Pseudoruminoccoccus, Rhodanobacter, Romboutsia, Sharpea, Varibaculum and Thioalkalivibrio.Quantification of viral abundance and analysis of the virome profileTo quantify the viral abundances in each sample, metagenomic reads were mapped to the gene set of VHGs (Supplementary Data 3) of each representative vOTU using Bowtie2 with a  > 95% identity threshold, and reads per kilobase million (RPKM) were calculated for each vOTU. The reason for using only VHGs in the analysis was to avoid over-counting of viral reads, which could be caused by spurious mapping of reads from horizontally transferred genes of other phages or bacterial species. The α-diversity (Shannon diversity) of the vOTU-level viral profile was calculated using the diversity function in the vegan package. The β-diversity (Bray-Curtis distance) between individuals was assessed using the vegdist function, and the average distance against other individuals was calculated for each individual. The VC-level viral profile was obtained by summing all the RPKM of vOTUs for each VC.Phylogenetic analysis of novel VCsTo construct phylogenetic trees for the vOTUs and reference genomes, protein sequences of large terminases, portal proteins, and major capsid proteins (Supplementary Data 3), which are often used to construct phage phylogenetic trees7,9, were extracted from the vOTUs in the 10 most abundant VCs (VC_19, 1, 2, 24, 12, 15, 3, 44, 18, 6), and their homologues were searched for in the reference phage genomes in RefSeq using DIAMOND with the more-sensitive option (e-value 0.01% (n = 865) and genera with average relative abundance >0.5% (n = 32) were included in the analysis.Analysis of VLPs and whole metagenomes from 24 faecal samplesQuality filtering of sequenced reads from the 24 VLPs and whole metagenomes was performed using fastp (version 0.20.1)92 with the default parameters. Contamination with human (hg38) or phiX genomes was excluded by mapping the reads to the genomes using Bowtie2.To exclude bacterial DNA contamination in the VLP dataset, we performed further filtering. First, the VLP reads were assembled into contigs using MEGAHIT and the contigs were checked for virus or not. Contigs were defined as viral contigs if they were predicted as viruses by DeepVirFinder (P-value More

Moore, J., Pope, J., Woods, M. & Ellis, A. 2018 Aerial Survey Results: California (USDA, 2018).Stephens, S. L. et al. Drought, tree mortality, and wildfire in forests adapted to frequent fire. Bioscience 68, 77–88 (2018).Article 

Google Scholar 
Li, S. & Banerjee, T. Spatial and temporal pattern of wildfires in California from 2000 to 2019. Sci. Rep. 11, 8779 (2021).Article 

Google Scholar 
Wang, D. et al. Economic footprint of California wildfires in 2018. Nat. Sustain. 4, 252–260 (2020).Article 

Google Scholar 
Asner, G. P. et al. Progressive forest canopy water loss during the 2012–2015 California drought. Proc. Natl Acad. Sci. USA 113, E249–E255 (2016).
Google Scholar 
Brodrick, P. G., Anderegg, L. D. L. & Asner, G. P. Forest drought resistance at large geographic scales. Geophys. Res. Lett. 46, 2752–2760 (2019).Article 

Google Scholar 
Jump, A. S. et al. Structural overshoot of tree growth with climate variability and the global spectrum of drought-induced forest dieback. Glob. Change Biol. 23, 3742–3757 (2017).Article 

Google Scholar 
Goulden, M. L. & Bales, R. C. California forest die-off linked to multi-year deep soil drying in 2012–2015 drought. Nat. Geosci. 12, 632–637 (2019).Article 

Google Scholar 
Paz-Kagan, T. et al. What mediates tree mortality during drought in the southern Sierra Nevada? Ecol. Appl. 27, 2443–2457 (2017).Article 

Google Scholar 
Trugman, A. T., Anderegg, L. D. L., Anderegg, W. R. L., Das, A. J. & Stephenson, N. L. Why is Tree Drought Mortality so Hard to Predict? Trends Ecol. Evol. 36, 520–532.(2021).Goodfellow, B. W. et al. The chemical, mechanical, and hydrological evolution of weathering granitoid. J. Geophys. Res. Earth Surf. 121, 1410–1435 (2016).Article 

Google Scholar 
Shen, X., Arson, C., Ferrier, K. L., West, N. & Dai, S. Mineral weathering and bedrock weakening: modeling microscale bedrock damage under biotite weathering. J. Geophys. Res. Earth Surf. 124, 2623–2646 (2019).Article 

Google Scholar 
McLaughlin, B. C. et al. Weather underground: subsurface hydrologic processes mediate tree vulnerability to extreme climatic drought. Glob. Change Biol. 26, 3091–3107 (2020).Article 

Google Scholar 
Hahm, W. J. et al. Low subsurface water storage capacity relative to annual rainfall decouples Mediterranean plant productivity and water use from rainfall variability. Geophys. Res. Lett. 46, 6544–6553 (2019).Article 

Google Scholar 
Zhang, Y., Keenan, T. F. & Zhou, S. Exacerbated drought impacts on global ecosystems due to structural overshoot. Nat. Ecol. Evol. 5, 1490–1498 (2021).Article 

Google Scholar 
Tague, C. & Peng, H. The sensitivity of forest water use to the timing of precipitation and snowmelt recharge in the California Sierra: implications for a warming climate. J. Geophys. Res. Biogeosci. 118, 875–887 (2013).Article 

Google Scholar 
Hahm, W. J., Riebe, C. S., Lukens, C. E. & Araki, S. Bedrock composition regulates mountain ecosystems and landscape evolution. Proc. Natl Acad. Sci. USA 111, 3338–3343 (2014).Article 

Google Scholar 
Uhlig, D., Schuessler, J. A., Bouchez, J., Dixon, J. L. & von Blanckenburg, F. Quantifying nutrient uptake as driver of rock weathering in forest ecosystems by magnesium stable isotopes. Biogeosciences 14, 3111–3128 (2017).Article 

Google Scholar 
Stone, E. C. Dew as an ecological factor: II. The effect of artificial dew on the survival of Pinus ponderosa and associated species. Ecology 38, 414–422 (1957).Article 

Google Scholar 
Wald, J. A., Graham, R. C. & Schoeneberger, P. J. Distribution and properties of soft weathered bedrock at ≤1 m depth in the contiguous United States. Earth Surf. Process. Landf. 38, 614–626 (2013).Article 

Google Scholar 
Klos, P. Z. et al. Subsurface plant-accessible water in mountain ecosystems with a Mediterranean climate. WIREs Water 5, e1277 (2018).Article 

Google Scholar 
Dawson, T. E., Hahm, W. J. & Crutchfield-Peters, K. Digging deeper: what the critical zone perspective adds to the study of plant ecophysiology. N. Phytol. 226, 666–671 (2020).Article 

Google Scholar 
Rempe, D. M. & Dietrich, W. E. Direct observations of rock moisture, a hidden component of the hydrologic cycle. Proc. Natl Acad. Sci. USA 115, 2664–2669 (2018).Article 

Google Scholar 
Holbrook, W. S. et al. Links between physical and chemical weathering inferred from a 65-m-deep borehole through Earth’s critical zone. Sci. Rep. 9, 4495 (2019).Article 

Google Scholar 
Krone, L. V. et al. Deep weathering in the semi-arid Coastal Cordillera, Chile. Sci. Rep. 11, 13057 (2021).Article 

Google Scholar 
Callahan, R. P. et al. Subsurface weathering revealed in hillslope‐integrated porosity distributions. Geophys. Res. Lett. 47, e2020GL088322 (2020).Holbrook, W. S. et al. Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory. Earth Surf. Process. Landf. 39, 366–380 (2014).Article 

Google Scholar 
Hayes, J. L., Riebe, C. S., Holbrook, W. S., Flinchum, B. A. & Hartsough, P. C. Porosity production in weathered rock: where volumetric strain dominates over chemical mass loss. Sci. Adv. 5, eaao0834 (2019).Article 

Google Scholar 
Riebe, C. S. et al. Anisovolumetric weathering in granitic saprolite controlled by climate and erosion rate. Geology 49, 551–555 (2021).Article 

Google Scholar 
McCormick, E. L. et al. Widespread woody plant use of water stored in bedrock. Nature 597, 225–229 (2021).Article 

Google Scholar 
Vitousek, P. M., Porder, S. & Houlton, B. Z. Terrestrial phosphorus limitation: mechanisms, implications, and nitrogen–phosphorus interactions. Ecol. Appl. 20, 5–15 (2010).Article 

Google Scholar 
Bateman, P. C., Dodge, F. C. W. & Bruggman, P. E. Major Oxide Analyses, CPIW Norms, Modes, and Bulk Specific Gravities of Plutonic Rocks from the Mariposa 1° × 2° Sheet, Central Sierra Nevada, California Open-File Report 84–162 (USGS, 1984).Amundson, R., Richter, D. D., Humphreys, G. S., Jobbagy, E. G. & Gaillardet, J. Coupling between biota and earth materials in the critical zone. Elements 3, 327–332 (2007).Article 

Google Scholar 
Tune, A. K., Druhan, J. L., Wang, J., Bennett, P. C. & Rempe, D. M. Carbon dioxide production in bedrock beneath soils substantially contributes to forest carbon cycling. J. Geophys. Res. Biogeosci. 125, e2020JG005795 (2020).Gabet, E. J. & Mudd, S. M. Bedrock erosion by root fracture and tree throw: a coupled biogeomorphic model to explore the humped soil production function and the persistence of hillslope soils. J. Geophys. Res. 115, F04005 (2010).Bateman, P. C. Plutonism in the Central Part of the Sierra Nevada Batholith, California (USGS, 1992); http://pubs.er.usgs.gov/publication/pp1483Callahan, R. P. et al. Arrested development: erosional equilibrium in the southern Sierra Nevada, California, maintained by feedbacks between channel incision and hillslope sediment production. GSA Bull. 131, 1179–1202 (2019).Article 

Google Scholar 
Flinchum, B. A. et al. Estimating the water holding capacity of the critical zone using near-surface geophysics. Hydrol. Process. 32, 3308–3326 (2018).Article 

Google Scholar 
St. Clair, J. Geophysical Investigations of Underplating at the Middle American Trench, Weathering in the Critical Zone, and Snow Water Equivalent in Seasonal Snow. PhD thesis, Univ. Wyoming (2015).Dvorkin, J. & Nur, A. Elasticity of high‐porosity sandstones: theory for two North Sea data sets. Geophysics 61, 1363–1370 (1996).Article 

Google Scholar 
Gu, X. et al. Seismic refraction tracks porosity generation and possible CO2 production at depth under a headwater catchment. Proc. Natl Acad. Sci. USA 117, 18991–18997 (2020).Article 

Google Scholar 
Pasquet, S., Holbrook, W. S., Carr, B. J. & Sims, K. W. W. Geophysical imaging of shallow degassing in a Yellowstone hydrothermal system. Geophys. Res. Lett. 43, 12,027–12,035 (2016).Article 

Google Scholar 
Dahlgren, R. A., Boettinger, J. L., Huntington, G. L. & Amundson, R. G. Soil development along an elevational transect in the western Sierra Nevada, California. Geoderma 78, 207–236 (1997).Article 

Google Scholar 
Stone, E. L. & Kalisz, P. J. On the maximum extent of tree roots. For. Ecol. Manage. 46, 59–102 (1991).Article 

Google Scholar 
Carlson, T. N. & Ripley, D. A. On the relation between NDVI, fractional vegetation cover, and leaf area index. Remote Sens. Environ. 62, 241–252 (1997).Article 

Google Scholar 
Goulden, M. L. et al. Evapotranspiration along an elevation gradient in California’s Sierra Nevada. J. Geophys. Res. 117, G03028 (2012).Ma, Q. et al. Wildfire controls on evapotranspiration in California’s Sierra Nevada. J. Hydrol. 590, 125364 (2020).Article 

Google Scholar 
Roche, J. W., Goulden, M. L. & Bales, R. C. Estimating evapotranspiration change due to forest treatment and fire at the basin scale in the Sierra Nevada, California. Ecohydrology 11, e1978 (2018).Bales, R. C. et al. Mechanisms controlling the impact of multi-year drought on mountain hydrology. Sci. Rep. 8, 690 (2018).Article 

Google Scholar 
Roy, D. P. et al. Characterization of Landsat-7 to Landsat-8 reflective wavelength and normalized difference vegetation index continuity. Remote Sens. Environ. 185, 57–70 (2016).Article 

Google Scholar 
Su, Y. et al. Emerging stress and relative resiliency of giant sequoia groves experiencing multiyear dry periods in a warming climate. J. Geophys. Res. Biogeosci. 122, 3063–3075 (2017).Article 

Google Scholar 
Moore, J., McAfee, L. & Iaccarino, J. 2016 Aerial Survey Results: California (USDA, 2017).Budyko, M. I., Miller, D. H. & Miller, D. H. Climate and Life (Academic Press, 1974).Hargreaves, G. H. & Samani, Z. A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1, 96–99 (1985).Article 

Google Scholar 
PRISM Climate Group PRISM Climate Data (Oregon State Univ., 2019).Bales, R. et al. Spatially distributed water-balance and meteorological data from the rain–snow transition, southern Sierra Nevada, California. Earth Syst. Sci. Data 10, 1795–1805 (2018).Article 

Google Scholar 
Callahan, R. P. Supplement for “Forest vulnerability to drought controlled by bedrock composition”. Hydroshare https://doi.org/10.4211/hs.edbb6ebfbc744186b5800932cd00b507 (2022).Earth Resources Observation and Science (EROS) Center USGS EROS Archive—Aerial Phorography—National Agriculture Imagery Program (NAIP) (USGS, 2017); https://doi.org/10.5066/F7QN651G More

Bedrock composition can play a critical role in determining the structure and water demand of forests, influencing their vulnerability to drought. The properties of bedrock can help explain within-region patterns of tree mortality in the 2011–2017 California drought.Montane forests are iconic natural resources that provide habitat, carbon sequestration, regulation of water, and, for many cultures, profound meaning. A warming climate and prolonged droughts threaten these forests, as shown by the 2011–2017 drought in California, USA, which killed over 140 million trees. However, the vulnerability of forests to climate-driven risks is not evenly distributed across these landscapes. In the 2011–2017 drought, some contiguous forested areas (or forest stands) suffered more than 70% mortality while forests in other locations experienced few or no losses1. Understanding these spatial patterns is critical for the projection of future risks and for targeted forest management. Writing in Nature Geoscience, Callahan and colleagues look beneath the surface at the composition of bedrock and find a link to these patterns of drought mortality in the California Sierra2. More

Richardson, A. D. et al. Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis. Glob. Change Biol. 18, 566–584 (2012).Article 

Google Scholar 
Keenan, T. F. et al. Net carbon uptake has increased through warming-induced changes in temperate forest phenology. Nat. Clim. Change 4, 598–604 (2014).CAS 
Article 

Google Scholar 
Piao, S. L. et al. Leaf onset in the northern hemisphere triggered by daytime temperature. Nat. Commun. 6, 6911 (2015).CAS 
Article 

Google Scholar 
Penuelas, J., Rutishauser, T. & Filella, I. Phenology feedbacks on climate change. Science 324, 887–888 (2009).CAS 
Article 

Google Scholar 
Garonna, I. et al. Strong contribution of autumn phenology to changes in satellite-derived growing season length estimates across Europe (1982–2011). Glob. Change Biol. 20, 3457–3470 (2014).Article 

Google Scholar 
Piao, S. L. et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451, 49–52 (2008).CAS 
Article 

Google Scholar 
Zhao, Y. et al. ABA receptor PYL9 promotes drought resistance and leaf senescence. Proc. Natl Acad. Sci. USA 113, 1949–1954 (2016).CAS 
Article 

Google Scholar 
Keskitalo, J., Bergquist, G., Gardestrom, P. & Jansson, S. A cellular timetable of autumn senescence. Plant Physiol. 139, 1635–1648 (2005).CAS 
Article 

Google Scholar 
Liu, Q. et al. Delayed autumn phenology in the Northern Hemisphere is related to change in both climate and spring phenology. Glob. Change Biol. 22, 3702–3711 (2016).Article 

Google Scholar 
Wu, C. Y. et al. Contrasting responses of autumn-leaf senescence to daytime and night-time warming. Nat. Clim. Change 8, 1092–1096 (2018).CAS 
Article 

Google Scholar 
Zani, D., Crowther, T. W., Mo, L., Renner, S. S. & Zohner, C. M. Increased growing-season productivity drives earlier autumn leaf senescence in temperate trees. Science 370, 1066–1071 (2020).CAS 
Article 

Google Scholar 
Zhang, Y., Parazoo, N. C., Williams, A. P., Zhou, S. & Gentine, P. Large and projected strengthening moisture limitation on end-of-season photosynthesis. Proc. Natl Acad. Sci. USA 117, 9216–9222 (2020).CAS 
Article 

Google Scholar 
Grossiord, C. et al. Plant responses to rising vapor pressure deficit. New Phytol. 226, 1550–1566 (2020).Article 

Google Scholar 
Ciais, P. et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 437, 529–533 (2005).CAS 
Article 

Google Scholar 
Keenan, T. F. & Richardson, A. D. The timing of autumn senescence is affected by the timing of spring phenology: implications for predictive models. Glob. Change Biol. 21, 2634–2641 (2015).Article 

Google Scholar 
Liu, L. B. et al. Soil moisture dominates dryness stress on ecosystem production globally. Nat. Commun. 11, 4892 (2020).CAS 
Article 

Google Scholar 
Delpierre, N. et al. Modelling interannual and spatial variability of leaf senescence for three deciduous tree species in France. Agric. For. Meteorol. 149, 938–948 (2009).Article 

Google Scholar 
Piao, S. L. et al. Weakening temperature control on the interannual variations of spring carbon uptake across northern lands. Nat. Clim. Change 7, 359–363 (2017).CAS 
Article 

Google Scholar 
Fu, Y. S. H. et al. Declining global warming effects on the phenology of spring leaf unfolding. Nature 526, 104–107 (2015).CAS 
Article 

Google Scholar 
Seastedt, T. R. & Knapp, A. K. Consequences of nonequilibrium resource availability across multiple time scales: the transient maxima hypothesis. Am. Nat. 141, 621–633 (1993).CAS 
Article 

Google Scholar 
Korner, C. Paradigm shift in plant growth control. Curr. Opin. Plant Biol. 25, 107–114 (2015).CAS 
Article 

Google Scholar 
Huxman, T. E. et al. Convergence across biomes to a common rain-use efficiency. Nature 429, 651–654 (2004).CAS 
Article 

Google Scholar 
McDowell, N. et al. Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought. New Phytol. 178, 719–739 (2008).Article 

Google Scholar 
Nolan, R. H. et al. Differences in osmotic adjustment, foliar abscisic acid dynamics, and stomatal regulation between an isohydric and anisohydric woody angiosperm during drought. Plant Cell Environ. 40, 3122–3134 (2017).CAS 
Article 

Google Scholar 
Fan, Y., Miguez-Macho, G., Jobbágy, E. G., Jackson, R. B. & Otero-Casal, C. Hydrologic regulation of plant rooting depth. Proc. Natl Acad. Sci. USA 114, 10572–10577 (2017).CAS 
Article 

Google Scholar 
Choat, B. et al. Triggers of tree mortality under drought. Nature 558, 531–539 (2018).CAS 
Article 

Google Scholar 
Giardina, F. et al. Tall Amazonian forests are less sensitive to precipitation variability. Nat. Geosci. 11, 405–409 (2018).CAS 
Article 

Google Scholar 
Kannenberg, S. A., Driscoll, A. W., Szejner, P., Anderegg, W. R. L. & Ehleringer, J. R. Rapid increases in shrubland and forest intrinsic water-use efficiency during an ongoing megadrought. Proc. Natl Acad. Sci. USA https://doi.org/10.1073/pnas.2118052118 (2021).Liu, Q. et al. Extension of the growing season increases vegetation exposure to frost. Nat. Commun. https://doi.org/10.1038/s41467-017-02690-y (2018).Schuur, E. A. G. et al. Climate change and the permafrost carbon feedback. Nature 520, 171–179 (2015).CAS 
Article 

Google Scholar 
Samaniego, L. et al. Anthropogenic warming exacerbates European soil moisture droughts. Nat. Clim. Change 8, 421–426 (2018).Article 

Google Scholar 
Templ, B. et al. Pan European Phenological database (PEP725): a single point of access for European data. Int. J. Biometeorol. 62, 1109–1113 (2018).Article 

Google Scholar 
Shen, M. et al. Increasing altitudinal gradient of spring vegetation phenology during the last decade on the Qinghai-Tibetan Plateau. Agric. For. Meteorol. 189, 71–80 (2014).Article 

Google Scholar 
Zhang, X. Y. Reconstruction of a complete global time series of daily vegetation index trajectory from long-term AVHRR data. Remote Sens. Environ. 156, 457–472 (2015).Article 

Google Scholar 
Chen, J. et al. A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky–Golay filter. Remote Sens. Environ. 91, 332–344 (2004).Article 

Google Scholar 
White, M. A. et al. Intercomparison, interpretation, and assessment of spring phenology in North America estimated from remote sensing for 1982–2006. Glob. Change Biol. 15, 2335–2359 (2009).Article 

Google Scholar 
Zhang, X. et al. Monitoring vegetation phenology using MODIS. Remote Sens. Environ. 84, 471–475 (2003).Article 

Google Scholar 
Gonsamo, A., Chen, J. M., Price, D. T., Kurz, W. A. & Wu, C. Y. Land surface phenology from optical satellite measurement and CO2 eddy covariance technique. J. Geophys. Res. 117, G03032 (2012).
Google Scholar 
Muñoz, S. ERA5-Land Monthly Averaged Data from 1981 to Present (C3S CDS, date accessed:10-8-2021); https://doi.org/10.24381/cds.68d2bb30Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191 (2018).Article 

Google Scholar 
Müller, W. A. et al. A Higher-resolution version of the Max Planck Institute Earth System Model (MPI-ESM1.2-HR). J. Adv. Model. Earth Syst. 10, 1383–1413 (2018).Article 

Google Scholar 
Vicente-Serrano, S. M. et al. Response of vegetation to drought time-scales across global land biomes. Proc. Natl Acad. Sci. USA 110, 52–57 (2013).CAS 
Article 

Google Scholar 
Allen, R. G., Smith, M., Pereira, L. S. & Perrier, A. An update for the calculation of reference evapotranspiration. ICID Bull. 43, 64–92 (1994).
Google Scholar 
Gampe, D. et al. Increasing impact of warm droughts on northern ecosystem productivity over recent decades. Nat. Clim. Change https://doi.org/10.1038/s41558-021-01112-8 (2021).Sheffield, J., Wood, E. F. & Roderick, M. L. Little change in global drought over the past 60 years. Nature 491, 435–438 (2012).CAS 
Article 

Google Scholar 
Peng, J., Wu, C. Y., Zhang, X. Y., Wang, X. Y. & Gonsamo, A. Satellite detection of cumulative and lagged effects of drought on autumn leaf senescence over the Northern Hemisphere. Glob. Change Biol. 25, 2174–2188 (2019).Article 

Google Scholar 
Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations—the CRU TS3.10 Dataset. Int. J. Climatol. 34, 623–642 (2014).Article 

Google Scholar 
Beaudoing, H., Rodell, M. & NASA/GSFC/HSL. GLDAS Noah Land Surface Model L4 3 Hourly 0.25 × 0.25 Degree Version 2.0 (GES DISC, 2015); https://doi.org/10.5067/342OHQM9AK6QBeaudoing, H., Rodell, M. & NASA/GSFC/HSL. GLDAS Noah Land Surface Model L4 3 Hourly 0.25 ×0.25 Degree Version 2.1 (GES DISC, 2016); https://doi.org/10.5067/E7TYRXPJKWOQZheng, Y. et al. Improved estimate of global gross primary production for reproducing its long-term variation, 1982–2017. Earth Syst. Sci. Data 12, 2725–2746 (2020).Article 

Google Scholar 
Zhang, K. et al. Vegetation greening and climate change promote multidecadal rises of global land evapotranspiration. Sci. Rep. https://doi.org/10.1038/srep15956 (2015).Li, Y. et al. Estimating global ecosystem isohydry/anisohydry using active and passive microwave satellite data. J. Geophys. Res. 122, 3306–3321 (2017).Article 

Google Scholar 
Moesinger, L. et al. The global long-term microwave Vegetation Optical Depth Climate Archive (VODCA). Earth Syst. Sci. Data 12, 177–196 (2020).Gupta, H. V., Kling, H., Yilmaz, K. K. & Martinez, G. F. Decomposition of the mean squared error and NSE performance criteria: implications for improving hydrological modelling. J. Hydrol. 377, 80–91 (2009).Article 

Google Scholar 
Botta, A., Viovy, N., Ciais, P., Friedlingstein, P. & Monfray, P. A global prognostic scheme of leaf onset using satellite data. Glob. Change Biol. 6, 709–725 (2000).Article 

Google Scholar  More

Izdebski, A. et al. Palaeoecological data indicates land-use changes across Europe linked to spatial heterogeneity in mortality during the Black Death pandemic. Nat. Ecol. Evol. 6, 297–306 (2022).CAS 
Article 

Google Scholar 
Benedictow, O. J. The Complete History of the Black Death (The Boydell Press, 2021).Palermo, L. Mercati del Grano a Roma tra Medioevo e Rinascimento. Il Mercato Distrettuale del Grano in Età Comunale (Istituto Nazionale di Studi Romani, 1990).Cortonesi, A. I cereali nell’Italia del tardo medioevo. Note sugli aspetti qualitativi del consumo. Riv. Stor. Agricol. 37, 3–30 (1997).
Google Scholar 
Nanni, P. in The Crisis of the 14th Century. Teleconnections Between Environmental and Societal Change? (eds Bauch M. & Schenk G. J.) 169–189 (De Gruyter, 2020).Lagerås, P. Environment, Society and the Black Death: An Interdisciplinary Approach to the Late-Medieval Crisis in Sweden (Oxbow Books, 2016).Roosen, J. & Curtis, D. The ‘light touch’ of the Black Death in the southern Netherlands: an urban trick? Econ. Hist. Rev. 72, 32–56 (2019).Article 

Google Scholar 
Preiser-Kapeller, J. Der Lange Sommer und die Kleine Eiszeit: Klima, Pandemien und der Wandel der Alten Welt 500–1500 n. Chr. (Mandelbaum, 2021).Sadori, L. The Lateglacial and Holocene vegetation and climate history of Lago di Mezzano (central Italy). Quat. Sci. Rev. 202, 30–44 (2018).Article 

Google Scholar 
Cortonesi, A. Ruralia. Economie e Paesaggi del Medioevo Italiano (Il Calamo, 1995).Cortonesi, A. L’olivo nell’Italia medievale. Reti Medievali Riv. 6, 1–29 (2005).
Google Scholar 
Mensing, S. A. et al. Historical ecology reveals landscape transformation coincident with cultural development in central Italy since the Roman Period. Sci. Rep. 8, 2138 (2018).Article 

Google Scholar 
Cortonesi, A. in Il Paesaggio Agrario Italiano Medievale: Storia e Didattica, 113–120 (Istituto Alcide Cervi, 2011). More

l-cry mutants show higher spawning synchrony than wild-type animals under non-natural light conditionsIn order to test for a functional involvement of L-Cry in monthly oscillator function, we generated two l-cry mutant alleles (Δ34 and Δ11bp) (Fig. 1a) using TALENs28. In parallel, we generated a monoclonal antibody against Platynereis L-Cry. By testing mutant versus wildtype worms with the anti-L-Cry antibody in Western blots (Fig. 1b) and immunohistochemistry (Fig. 1e–j), we verified the absence of L-Cry protein in mutants. Furthermore, we confirmed that the staining of the antibody in wildtype worms (Fig. 1e–h) matches the regions where l-cry mRNA is expressed (Fig. 1d). These tests confirmed that the engineered l-cry mutations result in loss-of-function alleles. In turn, they validate the specificity of the raised anti-L-Cry antibody.Fig. 1: l-cry–/– mutants are loss-of-function alleles.a Overview of the l-cry genomic locus for wt and mutants. Both mutant alleles result in an early frameshift and premature stop codons. The Δ34 allele has an additional 9 bp deletion in exon 3. b Western Blots of P. dumerilii heads probed with anti-L-Cry antibody. In the context of further investigations such Western blots of mutant versus wild types have been performed more than 10 times with highly consistent results. Also see further analyses in this manuscript and ref. 36. c overview of P. dumerilii. d whole mount in situ hybridization against l-cry mRNA on worm head. ae, anterior eye; pe, posterior eye. e–j Immunohistochemistry of premature wild-type (e–h) and mutant (i, j) worm heads sampled at zt19/20 using anti-L-Cry antibody (green) and Hoechst staining (magenta), dorsal views, anterior up. e, f: z-stack images (maximal projections of 50 layers, 1.28 µm each) in the area highlighted by the rectangle in (d), whereas (g–j) are single layer images of the area highlighted by the white rectangles in (e, f). In the context of further investigations such stainings of mutant versus wild types have been performed more than 10 times with highly consistent results. Also see further analyses in this manuscript and ref. 36.Full size imageWe next assessed the circalunar maturation timing of wild types and l-cry mutant populations in conventional culture conditions, i.e. worms grown under typical indoor room lighting (named here artificial sun- and moonlight, Supplementary Fig. 1b).We expected either no phenotype (if L-Cry was not involved in circalunar clock entrainment) or a decreased spawning precision (if L-Cry was functioning as moonlight receptor in circalunar clock entrainment). Instead we observed an increased precision of the entrained worm population:We analyzed the maturation data using two statistical approaches, linear and circular statistics. We used the classical linear plots5 and statistics to compare the monthly spawning data distribution (Fig. 2a–c, i). This revealed a clear difference between mutant animals, which exhibited a stronger spawning peak at the beginning of the NM phase, compared to their wildtype and heterozygous counterparts (Fig. 2a–c, Kolmogorov–Smirnov test on overall data distribution, Fig. 2i).Fig. 2: L-Cry shields the circalunar clock from light that is not naturalistic moonlight.a–d, j Spawning of l-cry +/+ (a), l-cry +/– (Δ34) (b) and l-cry −/−(Δ34/ Δ34) (c) animals over the lunar month in the lab with 8 nights of artificial moonlight (a–c), under natural conditions in the sea (d, replotted from ref. 34,50,) and in the lab using naturalistic sun- and moonlight (j, 8 nights moonlight). e–h, k Data as in (a–d, j) as circular plot. 360° correspond to 30 days of the lunar month. The arrow represents the mean vector, characterized by the direction angle µ and r (length of µ). r indicates phase coherence (measure of population synchrony). p-values inside the plots: result of Rayleigh Tests. Significance indicates non-random distribution of data points. The inner circle represents the Rayleigh critical value (p = 0.05). i–l Results of two-sided multisample statistics on spawning data shown in (a–h, j, k). The phase differences in days can be calculated from the angle between the two mean vectors (i.e. 12°= 1 day).Full size imageWe then analyzed the same data using circular statistics (as the monthly cycle is repeating, see details in Methods section), which allowed us to describe the data with the mean vector (defined by the direction angle µ and its length r, shown as arrows in Fig. 2e–g). The phase coherence r (ranging from 0 to 1) serves as a measure for synchrony of the population data. As expected for entrained populations, all genotypes distributed their spawning across a lunar month significantly different from random (Fig. 2e–g, p values in circles, Rayleigh’s Uniformity test29). In line with the observed higher spawning peak of the l-cry−/− mutants in the linear plots, the circular analysis revealed a significant difference in spawning distribution (Mardia–Watson-Wheeler test, for details see Methods section) and higher spawning synchrony of mutants (r = 0.614) than in wild types and heterozygotes (r = 0.295 and r = 0.222) (Fig. 2i). The specificity of this phenotype of higher spawning precision for l-cry homozygous mutants was confirmed by analyses on trans-heterozygous l-cry (Δ34/Δ11) mutants (Supplementary Fig. 2), and by the fact that such a phenotype is not detectable in any other light receptor mutant available in Platynereis (r-opsin130: Supplementary Fig. 3a, b, e, f, i; c-opsin131: Supplementary Fig. 3c, d, g, h, i, Go-opsin: refs. 32, 33).The higher spawning synchrony of l-cry mutants under artificial light mimics the spawning precision of wild-type at its natural habitatThis increased spawning precision of l-cry mutants under artificial (but conventional indoor) laboratory light conditions let us wonder about the actual population synchrony of the worms under truly natural conditions. The lunar spawning synchrony of P. dumerilii at the Bay of Naples (the origin of our lab culture) has been worked on for more than 100 y. This allowed us to re-investigate very detailed spawning data records from the worms’ natural habitat published prior to environmental/light pollution. For better accessibility and comparability we combined all months and replotted the data published in 192934 (Fig. 2d, h, I; see details in Methods section; r = 0.631). This analysis revealed that the higher spawning synchrony in l-cry–/– worms mimics the actual spawning synchrony of P. dumerillii populations in their natural habitat34 (compare Fig. 2c, g with 2d, h.)Given that recent, non-inbred isolates from the same habitat as our lab inbred strains (which is the same habitat as the data collected in ref. 34) exhibit a broad spawning distribution under standard worm culture light conditions (which includes the bright artificial moonlight)35, we hypothesized that the difference in spawning synchrony between wildtype laboratory cultures and populations in their natural habitat is caused by the rather bright nocturnal light stimulus typically used for the standard laboratory culture (Supplementary Fig. 1a vs. b).Lunar spawning precision of wild-type animals depends on naturalistic moonlight conditionsWe next tested the resulting prediction that naturalistic moonlight should increase the spawning precision of the wildtype population, using naturalistic sun- and moonlight devices we specifically designed based on light measurements at the natural habitat of P. dumerilii31 (Supplementary Fig. 1a, c). We assessed the impact of the naturalistic sun- and moonlight (Supplementary Fig. 1a, c) on wildtype animals, maintaining the temporal aspects of the lab light regime (i.e. 8 nights of “full moon”). Indeed, merely adjusting the light intensity to naturalistic conditions increased the precision and phase coherence of population-wide reproduction: After several months under naturalistic sun- and moonlight, wildtype worms spawned with a major peak highly comparable to the wildtype precision reported at its natural habitat (Fig. 2d, h vs. j, k), and also exhibited an increased population synchrony (r = 0.398 compared to r = 0.295 under standard worm room light conditions). This increased similarity to the spawning distribution at the natural habitat (“Sea”) is confirmed by statistical analyses (Fig. 2l): The phase difference (angle between the two mean vectors) is only one day (corresponding to 12°). In contrast, the spawning distribution of wild types under standard worm room light versus naturalistic light conditions is highly significantly different in linear and circular statistical tests and has a phase difference of 7.7 days (Fig. 2l).These findings show that it is the naturalistic light that is critical for a highly precise entrainment of the monthly clock of wild-type worms. Given that l-cry–/– animals reach this high precision with the artificial light (i.e. standard lab light) implies that in wildtype L-Cry blocks artificial, but not naturalistic full-moonlight from efficiently synchronizing the circalunar clock. This block is removed in l-cry–/– animals, leading to a better synchronization of the l-cry–/– population. This finding suggests that L-Cry’s major role could be that of a gatekeeper controlling which ambient light is interpreted as full-moonlight stimulus for circalunar clock entrainment.
l-cry functions as a light signal gatekeeper for circalunar clock entrainmentA prediction of this hypothesis is that mutants should entrain better to an artificial full-moonlight stimulus provided out-of-phase than their wild type counterparts (in which L-Cry should block the “wrong” moonlight at least partially from re-entraining the circalunar oscillator).We thus compared the spawning rhythms of l-cry+/+ and l-cry–/– worms under a re-entrainment paradigm, where we provided our bright artificial culture full-moonlight at the time of the subjective new moon phase (Fig. 3a). In order to compare the spawning data distribution relative to the initial full moon (FM) stimulus, as well as to the new full moon stimulus (i.e. new FM), we used two nomenclatures for the months: months with numbers are analyzed relative to the initial nocturnal light stimulus (i.e. FM), whereas months with letters are analyzed relative to the new (phase-shifted) nocturnal light stimulus (i.e. new FM, Fig. 3a). When the nocturnal light stimulus is omitted (to test for the oscillator function) we then refer to ‘free-running FM’ (FR-FM) or ‘new free-running FM’ (new FR-FM), respectively (Fig. 3a). Using these definitions, the efficiency of circalunar clock re-entrainment will be reflected in the similarity of spawning data distributions between month 1 and month D, i.e. the more similar the distribution, the more the population has shifted to the new phase.Fig. 3: l-cry−/− mutants entrain the circalunar clock faster than wt to a high-intensity artificial moonlight stimulus.a Nocturnal moonlight exposure protocol of lunar phase shift (entrained by 8 nights, phased shifted by 6 nights of artificial culture moon, light green). b, c Number of mature animals (percent per month, rolling mean with a window of 3 days) of l-cry wild-type (b) and homozygous mutant (c) animals. p-values indicate results of Kolomogorov–Smirnov tests. Dark blue arrowheads- old FM phase: wt show a spawning minimum, indicative that the worms are not properly phase shifted. Mutants spawn in high numbers, but don’t spawn at the old NM indicated by light blue arrowhead. Also compare to initial FM and NM in months 1,2. d, e Circular plots of the data shown in (b) and (c). Each circle represents one lunar month. Each dot represents one mature worm. The arrow represents the mean vector characterized by the direction angle µ and r. r (length of µ) indicates phase coherence (measure of population synchrony). The inner circle represents the Rayleigh critical value (p = 0.05). f, g Results of two-sided multisample statistics of data in (d, e). Phase differences in days can be calculated from the angle between the two mean vectors (i.e. 12°= 1 day).Full size imageWhen using the artificial nocturnal light conditions, the re-entrainment of l-cry–/– animals was both faster and more complete than for their wildtype relatives, as predicted from our gate keeper hypothesis. This is evident from the linear data analysis and Kolmogorov–Smirnov tests when comparing the month before the entrainment (month 1) with two months that should be shifted after the entrainment (months C,D, Fig. 3b, c, f, g).Most notably, while l-cry−/− worms were fully shifted in month D (Fig. 3c: compare boxes and see complete lack of spawning at the light blue arrowhead indicating the old NM/new FR-FM phase versus massive spawning at new NM phase around dark blue arrowhead), wildtype animals were still mostly spawning according to the initial lunar phase (Fig. 3b: compare boxes and see spawning at the light blue arrowhead versus almost lack of spawning at dark blue arrowhead). The faster re-entrainment of l-cry–/–, compared to l-cry+/+ animals is also confirmed by the Mardia–Watson-Wheeler test (see Methods section for details). For l-cry+/+ animals, the comparisons of the spawning distributions before and after re-entrainment show a 1000-fold (months 1 versus C) and tenfold (months 1 versus D) higher statistical significance difference than the corresponding comparisons for l-cry−/− worms (Fig. 3f, g). Consistently, the phase differences in days calculated from the angle between the two mean vectors from the circular analysis is smaller in the mutants than in the wild types when comparing the phase of the month before the entrainment (month 1) with two months after the entrainment (months C, D) (Fig. 3d–g). The fact that there are still differences in the mutant population before and after entrainment is likely due to the fact that even the mutants are not fully re-entrained. However, they have shifted more robustly in response to an artificial nocturnal light stimulus than the wild types. This provides further evidence that in wildtype worms L-Cry indeed blocks the “wrong” light from entering into the circalunar clock and thus functions as a light gatekeeper.L-Cry functions mainly as light interpreter, while its contribution as direct moonlight entraining photoreceptor is (at best) minorWe next tested to which extent L-Cry is itself a sensor for the re-entrainment signal under naturalistic light conditions. Based on the finding that l-cry−/− worms can still re-entrain the circalunar oscillator (see above), it is clear that even if L-Cry also directly contributed to the entrainment, it cannot be the only moonlight receptor mediating entrainment. With the experiments below, we aimed to test if L-Cry has any role as an entraining photoreceptor to the monthly oscillator.Thus, we tested how the circalunar clock is shifted in response to a re-entrainment with naturalistic moonlight in Platynereis wt versus l-cry−/− worms. For this, animals initially raised and entrained under standard worm room light conditions of artificial sun- and moonlight (Supplementary Fig. 1b, e) were challenged by a deviating FM stimulus of 8 nights of naturalistic moonlight (Fig. 4a, Supplementary Fig. 1c, e). This re-entraining stimulus was repeated for three consecutive months (Fig. 4a).Fig. 4: l-cry has a minor contribution as entraining photoreceptor to circalunar clock entrainment.a Nocturnal moonlight exposure protocol of lunar phase shift with 8 nights of naturalistic moonlight (dark green). Number of mature animals (percent per month, rolling mean with a window of 3 days) of l-cry wild-type (b) and mutant (c) animals. p-values: Kolomogorov–Smirnov tests. Black arrowheads indicate spawning-free intervals of the wildtype, which shifted to the position of the new FM (under free-running conditions: FR-FM). d, e Data as in (b, c) plotted as circular data. 360° correspond to 30 days of the lunar month. The arrow represents the mean vector characterized by the direction angle µ and r. r (length of µ) indicates phase coherence (measure of population synchrony). p values are results of Rayleigh Tests: Significance indicates non-random distribution of data points. The inner circle represents the Rayleigh critical value (p = 0.05). f, g Results of two-sided multisample statistics on spawning data shown in (a–e). Phase differences in days can be calculated from the angle between the two mean vectors (i.e. 12°= 1 day).Full size imageThe resulting spawning distribution was analyzed for the efficacy of the naturalistic moonlight to phase-shift the circalunar oscillator. In order to test if the animals had shifted their spawning to the new phase, we again compared the spawning pattern before the exposure to the new full moon stimulus (months with numbers: data distribution analyzed relative to the initial/old FM, see Fig. 4a for an overview) to the spawning pattern after the exposure to the new full moon stimulus (months with letters: data distribution analyzed relative to the new FM, Fig. 4a). The more similar the data distributions of month 1 is to the months C, D, the more the population was shifted to the new phase.The first re-entraining full moon stimulus (Fig. 4b, first dark green box) is given in the middle of the main spawning period. The nocturnal light itself does not cause immediate effects on the number of spawning worms (Fig. 4b, see also Fig. 2b, c), but the repeated exposure resulted in a noticeable shift of the spawning distribution indicating a phase shift of the monthly oscillator in wildtype. Already at the third re-entraining full moon stimulus, wildtype animals exhibited a completely shifted spawning pattern (Fig. 4b, d-d″, month 1, 2 vs. month C). This is supported by statistical analyses: When comparing the months 1 and 2 (relative to the old FM before the shift) to the month C (relative to the new FM after the shift), both the Kolmogorov–Smirnov test (Fig. 4b: gray rectangles, 4f) and the Mardia–Watson–Wheeler test of the same data were non-significant (Fig. 4f), indicative of the population shifting to the new phase. Consistently, the direction angle (µ) of the mean vectors before and after the shift was highly similar, resulting in a phase difference of only 0.2 days between months 1 and C and 0.5 days between month 2 and month C (Fig. 4f, for details see methods). The month under circalunar free-running conditions (month D) supports this observation, albeit with lower statistical support (Fig. 4b, d″, f).Of note, wild-type worms would eventually reach the high spawning precision found under naturalistic moonlight only after several more months based on independent experiments (Fig. 2j, k).When we analyzed the spawning distribution of l-cry mutants in the same way as the wild types, we found that the data distribution exhibited significant differences in the linear Kolmogorov–Smirnov test when comparing months 1 and 2 before the shift to the months C and D after the shift (Fig. 4c: gray rectangles, Fig. 4g); as well as in the phase distribution in the circular analyses when comparing the months before the shift (months 1 and 2) with the last months of the shift (months C,D) (Fig. 4e, e′ versus e″, e‴, g). The populations also exhibited a noticeable phase difference of ≥3.5 days (Fig. 4g).Based on the statistical significant difference in the re-entrainment of l-cry–/–, but not wild-type populations under a naturalistic sun- and moonlight regime, we conclude that L-Cry also likely contributes to circalunar entrainment as a photoreceptor. However, as these differences are rather minor, compared to the much stronger differences seen under artificial light regime, we conclude that its major role is the light gatekeeping function.In an independent study that focused on the impact of moonlight on daily timing, we identified r-Opsin1 as a lunar light receptor that mediates moonlight effects on the worms’ ~24 h clock36. We tested if r-opsin1 is similarly important for mediating the moonlight effects on the monthly oscillator of the worm, analyzed here. This is not the case. r-opsin1–/– animals re-entrain as well as wildtype worms under naturalistic light conditions (Supplementary Fig. 4). This adds to and is also consistent with our above observation that the spawning distribution is un-altered between r-opsin1–/– and wildtype animals under artificial light conditions (Supplementary Fig. 3a, b, e, f). This finding also further enforces the notion that monthly and daily oscillators use distinct mechanisms, but both require L-Cry as light interpreter.L-Cry discriminates between naturalistic sun- and moonlight by forming differently photoreduced statesGiven that the phenotype of l-cry–/– animals suggests a role of L-Cry as light gatekeeper, i.e. only allowing the ‘right’ light to most efficiently impact on the circalunar oscillator, we next investigated how this could function on the biochemical and cell biological level.While we have previously shown that Pdu-L-Cry is degraded upon light exposure in S2 cell culture15, it has remained unclear if L-Cry has the spectral properties and sensitivity to sense moonlight and whether this would differ from sunlight sensation. To test this, we purified full length L-Cry from insect cells (Supplementary Fig. 5a–c). Multi-angle light scattering (SEC-MALS) analyses of purified dark-state L-Cry revealed a molar mass of about 130 kDa, consistent with the formation of an L-Cry homodimer (theoretical molar mass of L-Cry monomer is 65.6 kDa) (Fig. 5a). Furthermore, purified L-Cry binds Flavin Adenine Dinucleotide (FAD) as its chromophore (Supplementary Fig. 5d, e). We then used UV/Vis absorption spectroscopy to analyze the FAD photoreaction of purified L-Cry in presence of 1 mM TCEP to prevent protein oxidation. The absorption spectrum of dark-state L-Cry showed maxima at 450 nm and 475 nm, consistent with the presence of oxidized FAD (Supplementary Fig. 5f, black line). As basic starting point to analyze its photocycle, L-Cry was photoreduced using a LED (PerkinElmer ACULED Dyo) with a blue-light dominated spectrum and spectral peak at 450 nm (Supplementary Fig. 1d, d′, henceforth referred to as “blue-light”) for 110 s37. The light-activated spectrum showed that blue-light irradiation of L-Cry leads to the complete conversion of FADox into an anionic FAD radical (FADo-) with characteristic FADo- absorption maxima at 370 nm and 404 nm and reduced absorbance at 450 nm (Supplementary Fig. 5f, blue spectrum, black arrows). In darkness, L-Cry reverted back to the dark-state with time constants of 2 min (18 °C), 4 min (6 °C) and 4.7 min (ice) (Supplementary Fig. 5g–k).Fig. 5: L-Cry forms differently photoreduced sunlight- and moonlight states.a Multi-Angle Light Scattering (MALS) analyses of dark-state L-Cry fractionated by size exclusion chromatography (SEC). Black dashed line: normalized UV absorbance, solid line: normalized scattering signal. The molar mass of about 130 kDa derived from MALS (mass signal shown in red) corresponds to an L-Cry homodimer. b Absorption spectrum of L-Cry in darkness (black) and after sunlight exposure (orange). Additional timepoints: Supplementary Fig. 6a. c Dark recovery of L-Cry after 20 min sunlight on ice. Absorbance at 450 nm in Supplementary Fig. 6b. d, e Absorption spectra of L-Cry after exposure to naturalistic moonlight for different durations. f Full spectra of dark recovery after 6 h moonlight. Absorbance at 450 nm: Supplementary Fig. 6d. g Absorption spectrum of L-Cry after 6 h of moonlight followed by 20 min of sunlight. h Absorption spectrum of L-Cry after 20 min sunlight followed by moonlight first results in dark-state recovery. Absorbance at 450 nm: Supplementary Fig. 6e. i Absorption spectrum of L-Cry after 20 min sunlight followed by 4 h and 6 h moonlight builds up the moonlight state. j Model of L-Cry responses to sunlight (orange), moonlight (green) and darkness (black). Only transitions between stably accumulating states are shown. Absorbances in (b–i) were normalized when a shift in the baseline occurred between different measurements of the same measurement set, which is then indicated on the Y-axis as “normalized absorbance”.Full size imageWe then investigated the response of L-Cry to ecologically relevant light, i.e. sun- and moonlight using naturalistic sun- and moonlight devices that we designed based on light measurements at the natural habitat of P. dumerilii31 (Supplementary Fig. 1a, c, e). Upon naturalistic sunlight illumination, FAD was photoreduced to FADo-, but with slower kinetics than under the stronger blue-light source, likely due to the intensity differences between the two lights (Supplementary Fig. 1c–e).While blue-light illumination led to a complete photoreduction within 110 s (Supplementary Fig. 5f), sunlight induced photoreduction to FADo- was completed after 5–20 min (Fig. 5b) and did not further increase upon continued illumation for up to 2 h (Supplementary Fig. 6a). Dark recovery kinetics had time constants of 3.2 min (18 °C) and 5 min (ice) (Fig. 1c, Supplementary Fig. 6b, c).As the absorbance spectrum of L-Cry overlaps with that of moonlight at the Platynereis natural habitat (Supplementary Fig. 1a), L-Cry has the principle spectral prerequisite to sense moonlight. However, the most striking characteristic of moonlight is its very low intensity (5.8 × 1010 photons/cm2/s at −5m, Supplementary Fig. 1a–e). To test if Pdu-L-Cry is sensitive enough for moonlight, we illuminated purified L-Cry with our custom-built naturalistic moonlight, closely resembling full-moonlight intensity and spectrum at the Platynereis natural habitat (Supplementary Fig. 1a, c, e). Naturalistic moonlight exposure up to 2.75 h did not markedly photoreduce FAD, notably there was no difference between 1 h and 2.75 h (Fig. 5d). However, further continuous naturalistic moonlight illumination of 4 h and longer resulted in significant changes (Fig. 5d), whereby the spectrum transitioned towards the light activated state of FADo- (note peak changes at 404 nm and at 450 nm). This photoreduction progressed further until 6 h naturalistic moonlight exposure (Fig. 5d). No additional photoreduction could be observed after 9 h and 12 h of naturalistic moonlight exposure (Fig. 5e), indicating a distinct state induced by naturalistic moonlight that reaches its maximum after ~6 h, when about half of the L-Cry molecules are photoreduced. This time of ~6 h is remarkably consistent with classical work showing that a minimum of ~6 h of continuous nocturnal light is important for circalunar clock entrainment, irrespective of the preceding photoperiod5. The dark recovery of L-Cry after 6 h moonlight exposure occurred with a time constant of 6.7 min at 18 °C (Fig. 5f, Supplementary Fig. 6d). Given that both sunlight and moonlight cause FAD photoreduction, but with different kinetics and different final FADo- product/FADox educt ratios, we wondered how purified L-Cry would react to transitions between naturalistic sun- and moonlight (i.e. during “sunrise” and “sunset”).Mimicking the sunrise scenario, L-Cry was first illuminated with naturalistic moonlight for 6 h followed by 20 min of sunlight exposure. This resulted in an immediate enrichment of the FADo- state (Fig. 5g). Hence, naturalistic sunlight immediately photoreduces remaining oxidized flavin molecules, that are characteristic of moonlight activated L-Cry, to FADo-, to reach a distinct fully reduced sunlight state.In contrast, when we next mimicked the day-night transition (“sunset”) by first photoreducing with naturalistic sunlight (or strong blue-light) and subsequently exposed L-Cry to moonlight, L-Cry first returned to its full dark-state within about 30 min (naturalistic sunlight: τ = 7 min (ice), Fig. 5h, Supplementary Fig. 6e; blue-light: τ = 9 min (ice), Supplementary Fig. 6f–h), despite the continuous naturalistic moonlight illumination. Prolonged moonlight illumination then led to the conversion of dark-state L-Cry to the moonlight state (Fig. 5i, Supplementary Fig. 6f). Hence, fully photoreduced sunlight-state L-Cry first has to return to the dark-state before accumulating the moonlight state characterized by the stable presence of the partial FADo- product/FADox educt. In contrast to sunlight-state L-Cry, moonlight-state L-Cry does not return to the oxidized (dark) state under naturalistic moonlight (Fig. 5e), i.e. moonlight maintains the moonlight state, but not the sunlight state. We note, that a partially photoreduced L-Cry state may be formed transiently during dark-state recovery of the sunlight state under moonlight. However, this transiently occurring partially photoreduced L-Cry state would differ from the “true” moonlight state (e.g. by an allosteric change) preventing its accumulation (see discussion and Supplementary Fig. 6i).Given that L-Cry forms a homodimer and moonlight photoreduces about half of the FAD molecules, we propose that the moonlight state corresponds to a half-reduced FADo- FADox dimer, where FAD is only photoreduced in one L-Cry monomer, whereas in the sunlight state both monomers are photoreduced (FADo- FADo-) (Fig. 5j). This implies that the quantum yield for FADox to FADo- photoreduction differs between the two L-Cry monomers. One monomer (referred to as “A” in Fig. 5j) acts as “very low intensity light sensor” with a high quantum yield ΦA. Hence, the very low photon number provided after 6 h of moonlight illumination is sufficient to photoreduce its flavin co-factor, resulting in the partially photoreduced FADo- FADox moonlight state (Fig. 5j).For direct comparison, our naturalistic moonlight’s emission (in the main absorbance range of L-Cry: 330 nm–510 nm) is 5.4 × 1010 photons/cm2/s (Supplementary Fig. 1e), which accumulates to ~1.2 × 1015 photons/cm2 in the 6 h required to reach the half-reduced moonlight state (Fig. 5d, e). For naturalistic sunlight, emitting ~7.5 × 1014 photons/cm2/s (330–510 nm), at least 5 min of sunlight illumination (i.e. > ~1.8 × 1017 photons/cm2) are required to photoreduce the flavin in both L-Cry monomers in order to reach the fully photoreduced FADo- FADo- sunlight state (Fig. 5b, j). Thus, the second L-Cry monomer (monomer “B” in Fig. 5j) has a significantly lower quantum yield ΦB for FAD photoreduction (ΦB  More