Ecology

Lynn, D. H. Ciliated Protozoa: Characterization, Classification, and Guide to the Literature 3rd edn, 1–455 (Springer, 2008).
Google Scholar 
Stoecker, D. K., Michaels, A. E. & Davis, L. H. Grazing by the jellyfish, Aurelia aurita, on microzooplankton. J. Plankton Res. 9, 901–915 (1987).Article 

Google Scholar 
Dolan, J. R., Vidussi, F. & Claustre, H. Planktonic ciliates in the Mediterranean Sea: Longitudinal trends. Deep-Sea Res. I(46), 2025–2039 (1999).Article 

Google Scholar 
Gómez, F. Trends on the distribution of ciliates in the open Pacific Ocean. Acta Oecol. 32, 188–202 (2007).Article 
ADS 

Google Scholar 
Azam, F. et al. The ecological role of water column microbes in the sea. Mar. Ecol. Prog. Ser. 10, 257–263 (1983).Article 
ADS 

Google Scholar 
Pierce, R. W. & Turner, J. T. Ecology of planktonic ciliates in marine food webs. Rev. Aquat. Sci. 6, 139–181 (1992).
Google Scholar 
Calbet, A. & Saiz, E. The ciliate-copepod link in marine ecosystems. Aquat. Microb. Ecol. 38, 157–167 (2005).Article 

Google Scholar 
Kim, Y. O. et al. Tintinnid species as biological indicators for monitoring intrusion of the warm oceanic waters into Korean coastal waters. Ocean Sci. J. 47, 161–172 (2012).Article 
ADS 

Google Scholar 
Wang, C. F. et al. Impact of the warm eddy on planktonic ciliate, with an emphasis on tintinnids as bioindicator species. Ecol. Indic. 133, 108441 (2021).Article 

Google Scholar 
Wang, C. F. et al. Planktonic tintinnid community structure variations in different water masses of the Arctic Basin. Front. Mar. Sci. 8, 775653 (2022).Article 

Google Scholar 
Haney, J. F. Diel patterns of zooplankton behavior. Bull. Mar. Sci. 43, 583–603 (1988).ADS 

Google Scholar 
Vaulot, D. & Marie, D. Diel variability of photosynthetic picoplankton in the equatorial Pacific. J. Geophys. Res-Oceans 104, 3297–3310 (1999).Article 
ADS 
CAS 

Google Scholar 
Hays, G. C., Webb, P. I. & Frears, S. L. Diet changes in the carbon and nitrogen content of the copepod Metridia lucens. J. Plankton Res. 4, 727–737 (1998).Article 

Google Scholar 
Hays, G. C., Harris, R. P. & Head, R. N. Diel changes in the near-surface biomass of zooplankton and the carbon content of vertical migrants. Deep-Sea Res. II(48), 1063–1068 (2001).ADS 

Google Scholar 
Anna, A., Enric, S. & Albert, C. Towards an understanding of diel feeding phythms in marine protists: Consequences of light manipulation. Microb. Ecol. 79, 64–72 (2020).Article 

Google Scholar 
Vaulot, D., Marie, D., Olson, R. J. & Chisholm, S. W. Growth of Prochlorococcus, a photosynthetic prokaryote, in the equatorial Pacific Ocean. Science 268, 1480–1482 (1995).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Binder, B. J. & DuRand, M. D. Diel cycles in surface waters of the equatorial Pacific. Deep-Sea Res. II(49), 2601–2617 (2002).ADS 

Google Scholar 
Li, C. L. et al. Quasi-antiphase diel patterns of abundance and cell size/biomass of picophytoplankton in the oligotrophic ocean. Geophys. Res. Lett. 49, e2022GL097753 (2022).ADS 

Google Scholar 
Ohman, M. D. The demographic benefits of diel vertical migration by zooplankton. Ecol. Monogr. 60, 257–281 (1990).Article 

Google Scholar 
Ringelberg, J. The photo behavior of Daphnia spp. as a model to explain diel vertical migration in zooplankton. Biol. Rev. 74, 397–423 (1999).Article 

Google Scholar 
Tarling, G. A., Jarvis, T., Emsley, S. M. & Matthews, J. B. L. Midnight sinking behaviour in Calanus finmarchicus: A response to satiation or krill predation?. Mar. Ecol. Prog. 240, 183–194 (2002).Article 

Google Scholar 
Cohen, J. H. & Forward, R. B. Diel vertical migration of the marine copepod Calanopia americana. I. Twilight DVM and its relationship to the diel light cycle. Mar. Biol. 147, 387–398 (2005).Article 

Google Scholar 
Cohen, J. H. & Forward, R. B. Diel vertical migration of the marine copepod Calanopia americana. II. Proximate role of exogenous light cues and endogenous rhythms. Mar. Biol. 147, 399–410 (2005).Article 

Google Scholar 
Ringelberg, J. Diel Vertical Migration of Zooplankton in Lakes and Oceans 1–347 (Springer, 2010).
Google Scholar 
Liu, H. J., Zhu, M. L., Guo, S. J., Zhao, X. H. & Sun, X. X. Effects of an anticyclonic eddy on the distribution and community structure of zooplankton in the South China Sea northern slope. J. Mar. Syst. 205, 103311 (2020).Article 

Google Scholar 
Tao, Z. C. et al. The diel vertical distribution and carbon biomass of the zooplankton community in the Caroline Seamount area of the western tropical Pacific Ocean. Sci. Rep. 12, 18908 (2022).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Dale, T. Diel vertical distribution of planktonic ciliates in Lindåspollene, Western Norway. Mar. Microb. Food Webs 2, 15–28 (1987).
Google Scholar 
Jonsson, P. R. Vertical distribution of planktonic ciliates–an experimental analysis of swimming behavior. Mar. Ecol. Prog. Ser. 52, 39–53 (1989).Article 
ADS 

Google Scholar 
Stocker, D. K., Taniguchi, A. & Michaels, A. E. Abundance of autotrophic, mixotrophic and heterotrophic ciliates in shelf and slope waters. Mar. Ecol. Prog. Ser. 50, 241–254 (1989).Article 
ADS 

Google Scholar 
Passow, U. Vertical migration of Gonyaulax catenata and Mesodinium rubrum. Mar. Biol. 110, 455–463 (1991).Article 

Google Scholar 
Suzuki, T. & Taniguchi, A. Temporal change of clustered distribution of planktonic ciliates in Toyama Bay in summers of 1989 and 1990. J. Oceanogr. 53, 35–40 (1997).Article 
CAS 

Google Scholar 
Olli, K. Diel vertical migration of phytoplankton and heterotrophic flagellates in the Gulf of Riga. J. Mar. Syst. 23, 145–163 (1999).Article 

Google Scholar 
Pérez, M. T., Dolan, J. R., Vidussi, F. & Fukai, E. Diel vertical distribution of planktonic ciliates within the surface layer of the NW Mediterrean (May 1995). Deep-Sea Res. I(47), 479–503 (2000).Article 

Google Scholar 
Rossberg, M. & Wickham, S. A. Ciliate vertical distribution and diel vertical migration in a eutrophic lake. Fund. Appl. Limnol. 171, 1–14 (2008).Article 

Google Scholar 
Gu, B. W. et al. High dynamics of ciliate community revealed via short-term, high-frequency sampling in a subtropical estuarine ecosystem. Front. Microbiol. 13, 797638 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Su, J. L. Overview of the South China Sea circulation and its influence on the coastal physical oceanography near the Pearl River Estuary. Cont. Shelf Res. 24, 1745–1760 (2004).Article 

Google Scholar 
Cravatte, S., Delcroix, T., Zhang, D., Mcphaden, M. & Leloup, J. Observed freshening and warming of the western pacific warm pool. Clim. Dyn. 33, 565–589 (2009).Article 

Google Scholar 
Feng, M. P., Zhang, W. C., Yu, Y., Xiao, T. & Sun, J. Horizontal distribution of tintinnids in the western South China Sea during summer 2007. J. Trop. Oceanogr. 32, 86–92 (2013).
Google Scholar 
Liu, H. X. et al. Composition and distribution of planktonic ciliates in the southern South China Sea during late summer: Comparison between surface and 75 m deep layer. J. Ocean Univ. China 15, 171–176 (2016).Article 
ADS 

Google Scholar 
Wang, C. F. et al. Vertical distribution of planktonic ciliates in the oceanic and slope areas of the western Pacific Ocean. Deep-Sea Res. II(167), 70–78 (2019).
Google Scholar 
Sun, P., Zhang, S. L., Wang, Y. & Huang, B. Q. Biogeographic role of the Kuroshio Current Intrusion in the microzooplankton community in the boundary zone of the northern South China Sea. Microorganisms 9, 1104 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Sohrin, R., Imazawa, M., Fukuda, H. & Suzuki, Y. Full-depth profiles of prokaryotes, heterotrophic nanoflagellates, and ciliates along a transect from the equatorial to the subarctic central Pacific Ocean. Deep-Sea Res. II(57), 1537–1550 (2010).ADS 

Google Scholar 
Wang, C. F. et al. Difference of planktonic ciliate communities of the tropical West Pacific, the Bering Sea and the Arctic Ocean. Acta Oceanol. Sin. 39, 9–17 (2020).Article 
CAS 

Google Scholar 
Wang, C. F. et al. Planktonic ciliate trait structure variation over Yap, Mariana and Caroline seamounts in the tropical western Pacific Ocean. J. Oceanol. Limnol. 39, 1705–1717 (2021).Article 
ADS 

Google Scholar 
McLaren, I. A. Demographic strategy of vertical migration by a marine copepod. Amer. Nat. 108, 91–102 (1974).Article 

Google Scholar 
Loose, C. J., Von Elert, E. & Dawidowicz, P. Chemically-induced diel vertical migration in Daphnia: A new bioassay for kairomones exuded by fish. Arch. Hydrobiol. 126, 329–337 (1993).Article 

Google Scholar 
Bandara, K., Varpe, Ø., Wijewardene, L., Tverberg, V. & Eiane, K. Two hundred years of zooplankton vertical migration research. Biol. Rev. 96, 1–43 (2021).Article 

Google Scholar 
Oubelkheir, K. & Sciandra, A. Diel variations in particle stocks in the oligotrophic waters of the Ionian Sea (Mediterranean). J. Mar. Syst. 74, 364–371 (2008).Article 

Google Scholar 
Yang, E. J., Choi, J. K. & Hyun, J. H. Distribution and structure of heterotrophic protist communities in the northeast equatorial Pacific Ocean. Mar. Biol. 146, 1–15 (2004).Article 

Google Scholar 
Wang, C. F. et al. Planktonic ciliate community structure and its distribution in the oxygen minimum zones in the Bay of Bengal (Eastern Indian Ocean). J. Sea Res. 190, 102311 (2022).Article 

Google Scholar 
Daro, M. H. Migratory and grazing behavior of copepods and vertical distribution of phytoplankton. Bull. Mar. Sci. 43, 710–729 (1988).
Google Scholar 
Ursella, L., Cardin, V., Batistić, M., Garić, R. & Gačić, M. Evidence of zooplankton vertical migration from continuous Southern Adriatic buoy current-meter records. Prog. Oceanogr. 167, 78–96 (2018).Article 
ADS 

Google Scholar 
Roman, M. R., Dam, H. G., Le Borgne, R. & Zhang, X. Latitudinal comparisons of equatorial Pacific zooplankton. Deep-Sea Res. II(49), 2695–2711 (2002).ADS 

Google Scholar 
Steinberg, D. K., Cope, J. S., Wilson, S. E. & Kobari, T. A comparison of mesopelagic mesozooplankton community structure in the subtropical and subarctic North Pacific Ocean. Deep-Sea Res. II(55), 1615–1635 (2008).ADS 

Google Scholar 
Isla, A., Scharek, R. & Latasa, M. Zooplankton diel vertical migration and contribution to deep active carbon flux in the NW Mediterranean. J. Mar. Syst. 143, 86–97 (2015).Article 

Google Scholar 
Dolan, J. R. Morphology and ecology in tintinnid ciliates of the marine plankton: Correlates of lorica dimensions. Acta Protozoologica 49, 235–244 (2010).
Google Scholar 
Jacquet, S., Partensky, F., Lennon, J. F. & Vaulot, D. Diel patterns of growth and division in marine picoplankton in culture. J. Phycol. 37, 357–369 (2001).Article 

Google Scholar 
Pitta, P., Giannakourou, A. & Christaki, U. Planktonic ciliates in the oligotrophic Mediterranean Sea: Longitudinal trends of standing stocks, distributions and analysis of food vacuole contents. Aquat. Microb. Ecol. 24, 297–311 (2001).Article 

Google Scholar 
Weisse, T. & Montagnes, D. J. S. Ecology of planktonic ciliates in a changing world: Concepts, methods, and challenges. J. Eukaryot. Microbiol. 69, e12879 (2022).Article 
PubMed 

Google Scholar 
Heinbokel, J. F. Diel periodicities and rates of reproduction in natural populations of tintinnines in the oligotrophic waters off Hawaii, September 1982. Mar. Microb. Food Webs 2, 1–14 (1987).
Google Scholar 
Tsai, A. Y., Chiang, K. P., Chang, J. & Gong, G. C. Seasonal diel variations of picoplankton and nanoplankton in a subtropical western Pacific coastal ecosystem. Limnol. Oceanogr. 50, 1221–1231 (2005).Article 
ADS 
CAS 

Google Scholar 
Ribalet, F. et al. Light-driven synchrony of Prochlorococcus growth and mortality in the subtropical Pacific gyre. Proc. Natl. Acad. Sci. U. S. A. 112, 8008–8012 (2015).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Connell, P. E., Ribalet, F., Armbrust, E. V., White, A. & Caron, D. A. Diel oscillations in the feeding activity of heterotrophic and mixotrophic nanoplankton in the North Pacific Subtropical Gyre. Aquat. Microb. Ecol. 85, 167–181 (2020).Article 

Google Scholar 
Cheung, K. C., Poon, B., Lan, C. Y. & Wong, M. H. Assessment of metal and nutrient concentrations in river water and sediment collected from the cities in the Pearl River Delta, South China. Chemosphere 52, 1431–1440 (2003).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Huang, X. P., Huang, L. M. & Yue, W. Z. The characteristics of nutrients and eutrophication in the Pearl River estuary. South China. Mar. Pollut. Bull. 47, 30–36 (2003).Article 
CAS 
PubMed 

Google Scholar 
Liu, S. M. et al. Nutrient dynamics in the winter thermohaline frontal zone of the northern shelf region of the South China Sea. J. Geophys. Res. 115, C11020 (2010).Article 
ADS 

Google Scholar 
Shu, Y. Q., Wang, Q. & Zu, T. T. Progress on shelf and slope circulation in the northern South China Sea. Sci. China Earth Sci. 61, 560–571 (2018).Article 
ADS 

Google Scholar 
Dai, S. et al. The effects of a warm-core eddy on chlorophyll a distribution and phytoplankton community structure in the northern South China Sea in spring 2017. J. Mar. Syst. 210, 103396 (2020).Article 

Google Scholar 
He, X. Q. et al. Eddy-entrained Pearl River plume into the oligotrophic basin of the South China Sea. Cont. Shelf Res. 124, 117–124 (2016).Article 
ADS 

Google Scholar 
Pan, X. J. et al. Remote sensing of surface [nitrite + nitrate] in river-influenced shelf-seas: The northern South China Sea Shelf-sea. Remote Sens. Environ. 210, 1–11 (2018).Article 
ADS 

Google Scholar 
Xu, J. et al. Phosphorus limitation in the northern South China Sea during late summer: Influence of the Pearl River. Deep-Sea Res. I. 55, 1330–1342 (2008).Article 
CAS 

Google Scholar 
Caron, D. Inorganic nutrients, bacteria, and the microbial loop. Microb. Ecol. 28, 295–298 (1994).Article 
CAS 
PubMed 

Google Scholar 
Kirchman, D. The uptake of inorganic nutrients by heterotrophic bacteria. Microb. Ecol. 28, 255–271 (1994).Article 
CAS 
PubMed 

Google Scholar 
Song, J. M. Biogeochemical Processes of Biogenic Elements in China Marginal Seas 1–657 (Springer, 2011).
Google Scholar 
Zhang, W. C. et al. Review of nutrient (nitrogen and phosphorus) regeneration in the marine pelagic microbial food web. Mar. Sci. Bull. 35, 241–251 (2016).CAS 

Google Scholar 
Ma, J. et al. Effects of Y3 seamount on nutrients influencing the ecological environment in the Western Pacific Ocean. Earth Sci. Front. 27, 322–331 (2020).
Google Scholar 
Li, H. B. et al. Tintinnid diversity in the tropical West Pacific Ocean. Acta Oceanol. Sin. 37, 218–228 (2018).Article 
CAS 

Google Scholar 
Dolan, J. R., Ritchie, M. E. & Ras, J. The, “neutral” community structure of planktonic herbivores, tintinnid ciliates of the microzooplankton, across the SE Tropical Pacific Ocean. Biogeosciences 4, 297–310 (2007).Article 
ADS 
CAS 

Google Scholar 
Dolan, J. R., Ritchie, M. E., Tunin-Ley, A. & Pizay, M. Dynamics of core and occasional species in the marine plankton: Tintinnid ciliates in the north-west Mediterranean Sea. J. Biogeogr. 36, 887–895 (2009).Article 

Google Scholar 
Dolan, J. R. & Marrasé, C. Planktonic ciliate distribution relative to a deep chlorophyll maximum: Catalan Sea, NW Mediterranean, June 1993. Deep-Sea Res. I(42), 1965–1987 (1995).Article 

Google Scholar 
Suzuki, T. & Taniguchi, A. Standing crops and vertical distribution of four groups of marine planktonic ciliates in relation to phytoplankton chlorophyll a. Mar. Biol. 132, 375–382 (1998).Article 

Google Scholar 
Utermöhl, H. Zur vervollkommnung der quantitativen phytoplankton Methodik. Mit. Int. Ver. Theor. Angew. Limnol. 9, 1–38 (1958).
Google Scholar 
Lund, J. W. G., Kipling, C. & Cren, E. D. L. The inverted microscope method of estimating algal numbers and the statistical basis of estimations by counting. Hydrobiologia 11, 143–170 (1958).Article 

Google Scholar 
Kofoid, C. A. & Campbell, A. S. A Conspectus of the Marine and Fresh-Water Ciliata Belonging to the Suborder Tintinnoinea: With Descriptions of New Species Principally from the Agassiz Expedition to the Eastern Tropical Pacific 1904–1905 (University of California Press, 1929).
Google Scholar 
Kofoid, C. A., & Campbell, A. S. Reports on the scientific results of the expedition to the eastern tropical Pacific, in charge to Alexander Agassiz, by US Fish commission steamer “Albatross”, from October 1904 to March 1905, The Ciliata: The Tintinnoinea (Bulletin of the Museum of Comparative Zoology of Harvard College), vol. XXXVII. Cambridge University, Harvard (Lieut.-Commander LM Garrett, USN commanding) (1939).Zhang, W. C., Feng, M. P., Yu, Y., Zhang, C. X. & Xiao, T. An Illustrated Guide to Contemporary Tintinnids in the World 1–499 (Science Press, 2012).
Google Scholar 
Paranjape, M. A. & Gold, K. Cultivation of marine pelagic protozoa. Ann. Inst. Oceanogr. Paris 58, 143–150 (1982).
Google Scholar 
Alder, V. A. Tintinnoinea. In South Atlantic zooplankton (ed. Boltovskoy, D.) 321–384 (Backhuys, 1999).
Google Scholar 
Verity, P. G. & Langdon, C. Relationships between lorica volume, carbon, nitrogen, and ATP content of tintinnids in Narragansett Bay. J. Plankton R. 6, 859–868 (1984).Article 
CAS 

Google Scholar 
Putt, M. & Stoecker, D. K. An experimentally determined carbon: Volume ratio for marine “oligotrichous” ciliates from estuarine and coastal waters. Limnol. Oceanogr. 34, 1097–1103 (1989).Article 
ADS 

Google Scholar 
Yu, Y. et al. Basin-scale variation in planktonic ciliate distribution: A detailed temporal and spatial study of the Yellow Sea. Mar. Biol. Res. 10, 641–654 (2014).Article 

Google Scholar 
Wang, C. F. et al. Hydrographic feature variation caused pronounced differences of planktonic ciliate community in the Pacific Arctic Region in summer 2016 and 2019. Front. Microbiol. 13, 881048 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Margalef, R. Information theory in ecology. Gen. Syst. 3, 36–71 (1958).
Google Scholar 
Shannon, C. E. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).Article 
MathSciNet 
MATH 

Google Scholar 
Dolan, J. R. & Pierce, R. W. Diversity and distributions of tintinnid ciliates. In The Biology and Ecology of Tintinnid Ciliates: Models for Marine Plankton (eds Dolan, J. R. et al.) 214–243 (Wiley-Blackwell, 2013).
Google Scholar 
Xu, Z. L. & Chen, Y. Q. Aggregated intensity of dominant species of zooplankton in autumn in the East China Sea. J. Ecol. 8, 13–15 (1989).
Google Scholar 
Anderson, M. J., Gorley, R. N. & Clarke, K. R. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods (PRIMER-E, 2008).
Google Scholar 
Jiang, Y., Xu, G. & Xu, H. Use of multivariate dispersion to assess water quality based on species composition data. Environ. Sci. Pollut. Res. 23, 3267–3272 (2016).Article 
CAS 

Google Scholar  More

When studying windfarm developments at the European region scale, the high densities of windfarm developments on blanket bog in Galicia and Greater Manchester (north England) are influenced by the total extent of the recognised blanket bog which is lower in Spain (31.2 km2 total extent) and in Greater Manchester (40.8 km2 total extent) in comparison to other regions (Fig. 1), although no relationship between the total extent of recognised blanket bog and the windfarm developments (wind turbines, tracks and total affected area) was found. Although the rest of the European regions across Spain showed lower densities of windfarm infrastructures (Fig. 2), the total extent of recognised blanket bogs across those regions was under 1 km2 (Fig. 1) meaning that the majority of recognised Spanish blanket bogs could be under threat due to their small size and the potential impact of windfarm infrastructures, if installed. In addition to this, previously unrecognised Spanish blanket bogs that have now been reported17 that could also be under pressure as the lack of formal recognition and protection leaves this habitat exposed to a range of anthropogenic activities, including windfarm developments. In fact, some examples of blanket bogs with extensive damage have been identified and reported in Galicia25, and more recently in Cantabrian blanket bogs17.Spanish unmapped areas of blanket bog at the edge-of-range of this habitat in the south of Europe are, therefore, particularly at risk from windfarm developments, and may disappear before their extent and importance can be defined40. Currently, new renewable energy regulations have been developed as a result of the climate emergency, and several windfarm developments have been proposed in ecologically sensitive areas, where blanket bogs have been reported (e.g. Sierra del Escudo, Spain) increasing the pressure on this habitat. Spanish blanket bogs also have specific characteristics, such as their small size as a consequence of the topographical limitations (e.g. slope) for their development26, meaning that they usually only cover the hill summits, where wind energy potential is at its greatest. Since blanket bogs are small and the windfarm development may cover all of the hill summit for their installation, many blanket bogs will be irrevocably damaged40.Most of the Galician blanket bogs were protected in 1999, under the Natura 2000 network and were declared as Special Area of Conservation (SAC) in 2014. However, between 1999 and 2012, Galician blanket bogs underwent severe and significant alterations in the peatland surface as a consequence of the large number of windfarm developments41 that were established during the period (Table A—Supplementary information), even when the site was incorporated into the Natura 2000 network (Table B—Supplementary information). Despite available scientific evidence that showed the potential environmental risks for these vulnerable ecosystems, windfarms were installed in what this work found to be the most extensive windfarm infrastructures across recognised European blanket bogs (Fig. 2).The incomplete current understanding of the extent of Spanish blanket bogs highlights the need to improve the completeness and representativeness of their current records across the Spanish Atlantic biogeographical region to include, within Natura 2000, a sufficient cover of their occupied area, in proportion to the representation of this natural habitat type in the Member state, for which it could therefore be concluded that the network is complete. Due to the increasing evidence highlighting how important the transitional areas are within the blanket bog complex42, other peatland types and wet heaths should be also considered when recognising and protecting blanket bogs. Mapping unrecorded blanket bogs must be a priority to fully understand the geographical and climatic range of this habitat, and obligatory protection under the Habitats Directive (92/43/EEC) is key to protecting the southern edge-of-range of this habitat.In addition to the lack of protection and updated inventories, the priority status included in the Habitats Directive, key to promoting their protection and restoration, is only for active blanket bogs, excluding other degraded blanket bogs with the potential to be active (carbon sinks), if they are restored. An approach similar to that of Scotland, where degraded blanket bogs are included33,39, could promote blanket bog restoration across Europe and improve the protection of this natural carbon storage.Many countries have also misinterpreted the active status of the blanket bog meaning that it is difficult to define whether the recognised blanket bog habitat is classed as a priority or not. Some countries, such as the Republic of Ireland, have classified as 7130 only active blanket bogs36, meaning that degraded blanket bogs lack appropriate classification and incorrectly applying the Habitat Directive designation as not all blanket bogs are included. The priority status is given when the habitat is particularly vulnerable or unique to the EU and necessitates additional measures for their protection and surveillance; however, whilst some blanket bogs may not act currently as carbon sinks, they still contain large amounts of carbon, and when restored they can recover their carbon sink function1, and then act to mitigate climate change.The issue of windfarm developments across the Republic of Ireland has been previously reported using a peat map43. However, despite researchers highlighting the importance of excluding vulnerable peatland ecosystems in future developments44, new areas of windfarms have been built affecting further recognised blanket bogs. At least 79 wind turbines have been installed in the Republic of Ireland since 2008 on recognised blanket bogs (Table A—Supplementary information) representing the 9.8% of the total onshore turbines installed in the country (Table 3), highlighting the importance of this conflict. The contribution of wind energy production to electricity supply was predicted to be up to 30% by 202044. In 2020, wind energy consumed in the Republic of Ireland represented 36%45. This represented an average annual increase of wind energy consumption of 16.9%45 between 2005 and 2020, which may explain part of the increase of 42% in wind turbines since 2008 (Table A—Supplementary information).Table 3 Total % of turbines on blanket bog (recognised/national inventories) in relation with the total turbines installed by country.Full size tableAcross Europe, several governments have developed climate action plans that over the next decade promote renewable energies to reduce carbon emissions. The government of the Republic of Ireland is aiming to generate up to 80% of electricity from renewable energy by 2030, providing support for onshore windfarm developments with an increase of up to 32% of the renewable energy production by 2030, but with a favourable preference for offshore wind energy production (up to 52% of the renewable energy production)46. This may help to reduce the conflict between blanket bogs and windfarm developments. Currently, windfarm annual energy production on blanket bogs accounts for 263.4 MW, 6.1% of the total production of wind energy in the Republic of Ireland47.The promotion of onshore wind energy production46 and the lack of protection of the full extent of blanket bogs are also threats that need to be considered in the Republic of Ireland. In 2008, a peat map was published showing the distribution of blanket bogs and raised bogs across the Republic of Ireland43. However, the inventory of current recognised blanket bogs under the Habitats Directive does not cover the full extent reported in this research43. While the total extent of recognised blanket bogs under the Habitats Directive 92/43/ECC reported a total of 3621 km2 of blanket bogs36, the real extent of blanket bogs across the country could be up to 2.5 times more (9202 km2)43, highlighting the lack of protection and the potential further increase of the windfarms and peatlands conflict in the Republic of Ireland as it happens in Spain and Scotland.The lack of recognition of blanket bog habitat in combination with the promotion of wind energy production across the island of Ireland could affect further areas of blanket bog, increasing the degradation of blanket bogs. An urgent review of inventories needs to be promoted in both countries, the Republic of Ireland and Northern Ireland, to fully assess the impact of the extensive areas of windfarms across the whole island.In Scotland, the pressure of windfarm developments on blanket bogs is also evident, where the Scottish Planning Policy considers classes 1 and 2 as areas of significant protection; although, windfarm developments may be possible under some circumstances48 as is permitted under the Habitats Directive across the EU29. However, to assess the impacts of windfarms on peatlands in a consistent way and evaluate the environmental impact of potential new developments on carbon-rich soils, a carbon calculator has been developed by the Scottish Government49. The carbon calculator allows users to estimate the carbon savings of windfarms installed on peatlands, although they highlight the importance of long-term management in relation to the final net carbon calculation49. Nonetheless, installing windfarms on non-degraded peatlands has been reported as unlikely to reduce carbon emissions even when the management has been considered carefully and it should be avoided 30. Therefore, peatlands under classes 1 and 2 considered by the Scottish government as a priority should be excluded from any windfarm developments (currently representing over 16% of onshore turbines, Table 3); especially considering the current policy of increasing onshore windfarms in Scotland50. Long-term research is needed to fully assess the impacts before new windfarm developments are installed.The difference between the recognised blanket bogs included in the EU Habitats Directive and the Scottish national inventory highlights the importance of updating and defining the complete extent of blanket bogs to facilitate their protection and restoration.In this novel research, the extent of windfarm developments across all recognised European blanket bogs under the Habitats Directive have been assessed. Large extents of blanket bogs have already been damaged, concentrated in the edge-of-range of this habitat and directly affecting hundreds of hectares of blanket bog across the rest of Europe. The full potential long-term damage to the habitat functionality is still unclear, but scientific evidence supports the negative impacts of windfarm developments on this critical habitat. European blanket bogs need further scientific evidence to demonstrate the real benefit of incentivising the reduction of carbon emissions by installing onshore windfarm infrastructures on peatlands which are causing the degradation of the most important long-term natural carbon sink and storage ecosystems. A strategic restoration plan and appropriate relevant legislation would be beneficial to promote the safeguarding of blanket bogs in the UK after Brexit. An urgent revision and compliance of the legislation regarding the protection of blanket bogs needs to be implemented, especially under the current trend of promotion and increasing legislation on renewable energy to reduce carbon emissions. An improvement of the national inventories across the EU and UK protected area networks is critical to implement the recognition, protection, and restoration of this habitat, in order to guarantee its favourable conservation status and its function as a long-term carbon sink to mitigate climate change. More

Decline in organism size is seen as a major biological response to climate change, and can be particularly pronounced in aquatic ectotherms such as fish, with subsequent implications for fishery yield and food security. However, as well as being modulated by climate factors, the fish population size structure can also be impacted by biotic (competition, predation) and other human factors (harvesting). For migrating species such as salmon, while smaller size may represent reduced size at maturity, it may also indicate faster maturation. More

Shuang Gao from Bjerkens Center for Climate Research in Norway, and colleagues from Germany and the United States explored future changes in marine primary production and carbon uptake under climate scenarios using the Norwegian Earth-system model, with four river transport configurations incorporating established future economic development and nutrient-use efficiency pathways. The researchers find that riverine nutrient inputs lessen nutrient limitation under warmer conditions. In the future, the effect of increased riverine carbon may be larger than the effect of nutrient inputs on the projections of ocean carbon uptake. In the historical period, increased nutrient inputs are considered the most prominent driver of carbon uptake. The results of this study are subject to model limitations, and high-resolution models should be used to assess the future impact. More

Boudreau BP, Huettel M, Forster S, Jahnke RA, McLachlan A, Middelburg JJ, et al. Permeable marine sediments: Overturning an old paradigm. Eos, Trans Am Geophys Union. 2001;82:133–6.Article 

Google Scholar 
Cook PL, Wenzhöfer F, Rysgaard S, Galaktionov OS, Meysman FJ, Eyre BD, et al. Quantification of denitrification in permeable sediments: Insights from a two‐dimensional simulation analysis and experimental data. Limnol Oceanogr Methods. 2006;4:294–307.Article 
CAS 

Google Scholar 
Evrard V, Glud RN, Cook PL. The kinetics of denitrification in permeable sediments. Biogeochemistry. 2013;113:563–72.Article 

Google Scholar 
Huettel M, Berg P, Kostka JE. Benthic exchange and biogeochemical cycling in permeable sediments. Ann Rev Marine Sci. 2014;6:23–51.Article 

Google Scholar 
Rao AMF, McCarthy MJ, Gardner WS, Jahnke RA. Respiration and denitrification in permeable continental shelf deposits on the South Atlantic Bight: Rates of carbon and nitrogen cycling from sediment column experiments. Continental Shelf Res. 2007;27:1801–19.Article 

Google Scholar 
Billerbeck M, Werner U, Polerecky L, Walpersdorf E, de Beer D, Huettel M Surficial and deep pore water circulation governs spatial and temporal scales of nutrient recycling in intertidal sand flat sediment. Marine Ecol Progr Series. 2006;326:61–76.Jansen S, Walpersdorf E, Werner U, Billerbeck M, Böttcher ME, de Beer D. Functioning of intertidal flats inferred from temporal and spatial dynamics of O2, H2S and pH in their surface sediment. Ocean Dyn. 2009;59:317–32.Article 

Google Scholar 
de Beer D, Wenzhöfer F, Ferdelman TG, Boehme SE, Huettel M, van Beusekom JEE, et al. Transport and mineralization rates in North Sea sandy intertidal sediments, Sylt-Rømø Basin, Wadden Sea. Limnol Oceanogr. 2005;50:113–27.Article 

Google Scholar 
Gao H, Matyka M, Liu B, Khalili A, Kostka JE, Collins G, et al. Intensive and extensive nitrogen loss from intertidal permeable sediments of the Wadden Sea. Limnol Oceanogr. 2012;57:185–98.Article 
CAS 

Google Scholar 
Gao H, Schreiber F, Collins G, Jensen MM, Kostka JE, Lavik G, et al. Aerobic denitrification in permeable Wadden Sea sediments. ISME J. 2009;4:417.Article 
PubMed 

Google Scholar 
Elliott AH, Brooks NH. Transfer of nonsorbing solutes to a streambed with bed forms: Theory. Water Resour Res. 1997;33:123–36.Article 
CAS 

Google Scholar 
Precht E, Huettel M. Advective pore‐water exchange driven by surface gravity waves and its ecological implications. Limnol Oceanogr. 2003;48:1674–84.Article 

Google Scholar 
Ahmerkamp S, Marchant HK, Peng C, Probandt D, Littmann S, Kuypers MM, et al. The effect of sediment grain properties and porewater flow on microbial abundance and respiration in permeable sediments. Sci Rep. 2020;10:1–12.Article 

Google Scholar 
Ahmerkamp S, Winter C, Krämer K, Beer DD, Janssen F, Friedrich J, et al. Regulation of benthic oxygen fluxes in permeable sediments of the coastal ocean. Limnol Oceanogr. 2017;62:1935–54.Article 
CAS 

Google Scholar 
Cardenas MB, Wilson JL Dunes, turbulent eddies, and interfacial exchange with permeable sediments. Water Resour Res. 2007;43:W08412.Santos IR, Eyre BD, Huettel M. The driving forces of porewater and groundwater flow in permeable coastal sediments: a review. Estuarine, Coastal Shelf Sci. 2012;98:1–15.Article 

Google Scholar 
Probandt D, Eickhorst T, Ellrott A, Amann R, Knittel K. Microbial life on a sand grain: from bulk sediment to single grains. ISME J. 2018;12:623–33.Article 
PubMed 

Google Scholar 
Marchant HK, Ahmerkamp S, Lavik G, Tegetmeyer HE, Graf J, Klatt JM, et al. Denitrifying community in coastal sediments performs aerobic and anaerobic respiration simultaneously. ISME J. 2017;11:1799.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Marchant HK, Holtappels M, Lavik G, Ahmerkamp S, Winter C, Kuypers MMM. Coupled nitrification–denitrification leads to extensive N loss in subtidal permeable sediments. Limnol Oceanogr. 2016;61:1033–48.Article 

Google Scholar 
Marchant HK, Tegetmeyer HE, Ahmerkamp S, Holtappels M, Lavik G, Graf J, et al. Metabolic specialization of denitrifiers in permeable sediments controls N2O emissions. Environ Microbiol. 2018;20:4486–502.Article 
CAS 
PubMed 

Google Scholar 
Laverman AM, Pallud C, Abell J, Van, Cappellen P. Comparative survey of potential nitrate and sulfate reduction rates in aquatic sediments. Geochimica et Cosmochimica Acta. 2012;77:474–88.Article 
CAS 

Google Scholar 
Fenchel T, Jørgensen B. Detritus food chains of aquatic ecosystems: the role of bacteria. Adv Microb Ecol. 1977;1:1–58.Article 
CAS 

Google Scholar 
Canfield DE, Kristensen E, Thamdrup B Aquatic Geomicrobiology: Elsevier Science; 2005.Froelich PN, Klinkhammer G, Bender ML, Luedtke N, Heath GR, Cullen D, et al. Early oxidation of organic matter in pelagic sediments of the eastern equatorial Atlantic: suboxic diagenesis. Geochimica et cosmochimica Acta. 1979;43:1075–90.Article 
CAS 

Google Scholar 
Eckford RE, Fedorak PM. Chemical and microbiological changes in laboratory incubations of nitrate amendment “sour” produced waters from three western Canadian oil fields. J Ind Microbiol Biotechnol. 2002;29:243–54.Article 
CAS 
PubMed 

Google Scholar 
Grigoryan AA, Cornish SL, Buziak B, Lin S, Cavallaro A, Arensdorf JJ, et al. Competitive oxidation of volatile fatty acids by sulfate- and nitrate-reducing bacteria from an oil field in Argentina. Appl Environ Microbiol. 2008;74:4324.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hubert C, Nemati M, Jenneman G, Voordouw G. Containment of biogenic sulfide production in continuous up-flow packed-bed bioreactors with nitrate or nitrite. Biotechnol Progr. 2003;19:338–45.Article 
CAS 

Google Scholar 
Greene EA, Hubert C, Nemati M, Jenneman GE, Voordouw G. Nitrite reductase activity of sulphate-reducing bacteria prevents their inhibition by nitrate-reducing, sulphide-oxidizing bacteria. Environ Microbiol. 2003;5:607–17.Article 
CAS 
PubMed 

Google Scholar 
Wolfe BM, Lui SM, Cowan JA. Desulfoviridin, a multimeric-dissimilatory sulfite reductase from Desulfovibrio vulgaris (Hildenborough) Purification, characterization, kinetics and EPR studies. Eur J Biochem. 1994;223:79–89.Article 
CAS 
PubMed 

Google Scholar 
Fossing H, Gallardo VA, Jørgensen BB, Hüttel M, Nielsen LP, Schulz H, et al. Concentration and transport of nitrate by the mat-forming sulphur bacterium Thioploca. Nature. 1995;374:713–15.Article 
CAS 

Google Scholar 
Jørgensen BB. Big sulfur bacteria. ISME J. 2010;4:1083.Article 
PubMed 

Google Scholar 
Marzocchi U, Trojan D, Larsen S, Louise Meyer R, Peter Revsbech N, Schramm A, et al. Electric coupling between distant nitrate reduction and sulfide oxidation in marine sediment. ISME J. 2014;8:1682.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Londry KL, Suflita JM. Use of nitrate to control sulfide generation by sulfate-reducing bacteria associated with oily waste. J Ind Microbiol Biotechnol. 1999;22:582–9.Article 
CAS 
PubMed 

Google Scholar 
McInerney MJ, Bhupathiraju VK, Sublette KL. Evaluation of a microbial method to reduce hydrogen sulfide levels in a porous rock biofilm. J Ind Microbiol. 1992;11:53–8.Article 
CAS 

Google Scholar 
Schwermer CU, Ferdelman TG, Stief P, Gieseke A, Rezakhani N, Van Rijn J, et al. Effect of nitrate on sulfur transformations in sulfidogenic sludge of a marine aquaculture biofilter. FEMS Microbiol Ecol. 2010;72:476–84.Article 
CAS 
PubMed 

Google Scholar 
Thamdrup B, Fossing H, Jørgensen BB. Manganese, iron and sulfur cycling in a coastal marine sediment, Aarhus Bay, Denmark. Geochimica et Cosmochimica Acta. 1994;58:5115–29.Article 
CAS 

Google Scholar 
Al-Raei AM, Bosselmann K, Böttcher ME, Hespenheide B, Tauber F. Seasonal dynamics of microbial sulfate reduction in temperate intertidal surface sediments: controls by temperature and organic matter. Ocean Dyn. 2009;59:351–70.Article 

Google Scholar 
Dyksma S, Pjevac P, Ovanesov K, Mussmann M. Evidence for H2 consumption by uncultured Desulfobacterales in coastal sediments. Environ Microbiol. 2018;20:450–61.Article 
CAS 
PubMed 

Google Scholar 
Musat N, Werner U, Knittel K, Kolb S, Dodenhof T, van Beusekom JEE, et al. Microbial community structure of sandy intertidal sediments in the North Sea, Sylt-Rømø Basin, Wadden Sea. Syst Appl Microbiol. 2006;29:333–48.Article 
PubMed 

Google Scholar 
Mußmann M, Ishii K, Rabus R, Amann R. Diversity and vertical distribution of cultured and uncultured Deltaproteobacteria in an intertidal mud flat of the Wadden Sea. Environ Microbiol. 2005;7:405–18.Article 
PubMed 

Google Scholar 
Dyksma S, Lenk S, Sawicka JE, Mußmann M. Uncultured gammaproteobacteria and desulfobacteraceae account for major acetate assimilation in a coastal marine sediment. Front Microbiol. 2018;9:3124Article 
PubMed 
PubMed Central 

Google Scholar 
Chen J, Hanke A, Tegetmeyer HE, Kattelmann I, Sharma R, Hamann E, et al. Impacts of chemical gradients on microbial community structure. ISME J. 2017;11:920.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Saad S, Bhatnagar S, Tegetmeyer HE, Geelhoed JS, Strous M, Ruff SE. Transient exposure to oxygen or nitrate reveals ecophysiology of fermentative and sulfate‐reducing benthic microbial populations. Environ Microbiol. 2017;19:4866–81.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Brunet RC, Garcia-Gil LJ. Sulfide-induced dissimilatory nitrate reduction to ammonia in anaerobic freshwater sediments. FEMS Microbiol Ecol. 1996;21:131–8.Article 
CAS 

Google Scholar 
Murphy AE, Bulseco AN, Ackerman R, Vineis JH, Bowen JL. Sulphide addition favours respiratory ammonification (DNRA) over complete denitrification and alters the active microbial community in salt marsh sediments. Environ Microbiol. 2020;22:2124–39.Article 
CAS 
PubMed 

Google Scholar 
Krekeler D, Cypionka H. The preferred electron acceptor of Desulfovibrio desulfuricans CSN. FEMS Microbiol Ecol. 1995;17:271–7.Article 
CAS 

Google Scholar 
Seitz H-J, Cypionka H. Chemolithotrophic growth of Desulfovibrio desulfuricans with hydrogen coupled to ammonification of nitrate or nitrite. Arch Microbiol. 1986;146:63–7.Article 
CAS 

Google Scholar 
Dalsgaard T, Bak F. Nitrate reduction in a sulfate-reducing bacterium, Desulfovibrio desulfuricans, isolated from rice paddy soil: sulfide inhibition, kinetics, and regulation. Appl Environ Microbiol. 1994;60:291–7.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Marietou A. Nitrate reduction in sulfate-reducing bacteria. FEMS Microbiol Lett. 2016;363:fnw155.Article 
PubMed 

Google Scholar 
Marietou A, Griffiths L, Cole J. Preferential reduction of the thermodynamically less favorable electron acceptor, sulfate, by a nitrate-reducing strain of the sulfate-reducing bacterium Desulfovibrio desulfuricans 27774. J Bacteriol. 2009;191:882–889.Article 
CAS 
PubMed 

Google Scholar 
Korte HL, Saini A, Trotter VV, Butland GP, Arkin AP, Wall JD. Independence of nitrate and nitrite inhibition of Desulfovibrio vulgaris Hildenborough and use of nitrite as a substrate for growth. Environ Sci Technol. 2015;49:924–931.Article 
CAS 
PubMed 

Google Scholar 
Pereira IA, LeGall J, Xavier AV, Teixeira M. Characterization of a heme c nitrite reductase from a non-ammonifying microorganism, Desulfovibrio vulgaris Hildenborough. Biochimica et Biophysica Acta (BBA)-Protein Struct Mol Enzymol. 2000;1481:119–130.Article 
CAS 

Google Scholar 
Werner U, Billerbeck M, Polerecky L, Franke U, Huettel M, van Beusekom JEE, et al. Spatial and temporal patterns of mineralization rates and oxygen distribution in a permeable intertidal sand flat (Sylt, Germany). Limnol Oceanogr. 2006;51:2549–63.Article 
CAS 

Google Scholar 
Marchant HK, Lavik G, Holtappels M, Kuypers MMM. The fate of nitrate in intertidal permeable sediments. PLOS ONE. 2014;9:e104517.Article 
PubMed 
PubMed Central 

Google Scholar 
Canfield DE. Reactive iron in marine sediments. Geochimica et cosmochimica acta. 1989;53:619–632.Article 
CAS 
PubMed 

Google Scholar 
Billerbeck M, Werner U, Bosselmann K, Walpersdorf E, Huettel M. Nutrient release from an exposed intertidal sand flat. Marine Ecol Progr Series. 2006;316:35–51.Article 
CAS 

Google Scholar 
Canfield DE, Stewart FJ, Thamdrup B, De Brabandere L, Dalsgaard T, Delong EF, et al. A cryptic sulfur cycle in oxygen-minimum–zone waters off the Chilean Coast. Science. 2010;330:1375.Article 
CAS 
PubMed 

Google Scholar 
Jørgensen BB, Findlay AJ, Pellerin A. The biogeochemical sulfur cycle of marine sediments. Front Microbiol. 2019;10:849.Article 
PubMed 
PubMed Central 

Google Scholar 
Nemati M, Mazutinec TJ, Jenneman GE, Voordouw G. Control of biogenic H2S production with nitrite and molybdate. J Ind Microbiol Biotechnol. 2001;26:350–355.Article 
CAS 
PubMed 

Google Scholar 
Haveman SA, Greene EA, Stilwell CP, Voordouw JK, Voordouw G. Physiological and Gene Expression Analysis of Inhibition of Desulfovibrio vulgaris Hildenborough by Nitrite. J Bacteriol. 2004;186:7944–7950.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Behrendt A, de Beer D, Stief P. Vertical activity distribution of dissimilatory nitrate reduction in coastal marine sediments. Biogeosciences. 2013;10:7509–23.Article 

Google Scholar 
Findlay AJ, Pellerin A, Laufer K, Jørgensen BB. Quantification of sulphide oxidation rates in marine sediment. Geochimica et Cosmochimica Acta. 2020;280:441–52.Article 
CAS 

Google Scholar 
Waite DW, Chuvochina M, Pelikan C, Parks DH, Yilmaz P, Wagner M, et al. Proposal to reclassify the proteobacterial classes Deltaproteobacteria and Oligoflexia, and the phylum Thermodesulfobacteria into four phyla reflecting major functional capabilities. Int J Syst Evolut Microbiol. 2020;70:5972–6016.Article 
CAS 

Google Scholar 
Dyksma S, Bischof K, Fuchs BM, Hoffmann K, Meier D, Meyerdierks A, et al. Ubiquitous Gammaproteobacteria dominate dark carbon fixation in coastal sediments. ISME J. 2016;10:1939–1953.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lenk S, Arnds J, Zerjatke K, Musat N, Amann R, Mußmann M. Novel groups of Gammaproteobacteria catalyse sulfur oxidation and carbon fixation in a coastal, intertidal sediment. Environ Microbiol. 2011;13:758–774.Article 
CAS 
PubMed 

Google Scholar 
An S, Gardner W S. Dissimilatory nitrate reduction to ammonium (DNRA) as a nitrogen link, versus denitrification as a sink in a shallow estuary (Laguna Madre/Baffin Bay, Texas). Marine Ecol Progr Series. 2002;237:41–50.Article 
CAS 

Google Scholar 
Wankel SD, Ziebis W, Buchwald C, Charoenpong C, de Beer D, Dentinger J, et al. Evidence for fungal and chemodenitrification based N2O flux from nitrogen impacted coastal sediments. Nat Commun. 2017;8:15595.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Moura I, Bursakov S, Costa C, Moura JJ. Nitrate and nitrite utilization in sulfate-reducing bacteria. Anaerobe. 1997;3:279–290.Article 
CAS 
PubMed 

Google Scholar 
Song G, Liu S, Zhang J, Zhu Z, Zhang G, Marchant HK, et al. Response of benthic nitrogen cycling to estuarine hypoxia. Limnol Oceanogr. 2021;66:652–66.Article 
CAS 

Google Scholar 
Tiedje JM, Sexstone AJ, Myrold DD, Robinson JA. Denitrification: ecological niches, competition and survival. Antonie van Leeuwenhoek. 1983;48:569–583.Article 

Google Scholar 
Strohm TO, Griffin B, Zumft WG, Schink B. Growth yields in bacterial denitrification and nitrate ammonification. Appl Environ Microbiol. 2007;73:1420–1424.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rütting T, Boeckx P, Müller C, Klemedtsson L. Assessment of the importance of dissimilatory nitrate reduction to ammonium for the terrestrial nitrogen cycle. Biogeosciences. 2011;8:1779–91.Article 

Google Scholar 
Røy H, Lee JS, Jansen S, de Beer D. Tide-driven deep pore-water flow in intertidal sand flats. Limnol Oceanogr. 2008;53:1521–30.Article 

Google Scholar 
Cline JD. Spectrophotometric determination of hydrogen sulfide in natural waters 1. Limnol Oceanogr. 1969;14:454–458.Article 
CAS 

Google Scholar 
Viollier E, Inglett P, Hunter K, Roychoudhury A, Van, Cappellen P. The ferrozine method revisited: Fe (II)/Fe (III) determination in natural waters. Appl Geochem. 2000;15:785–90.Article 
CAS 

Google Scholar 
Røy H, Weber HS, Tarpgaard IH, Ferdelman TG, Jørgensen BB. Determination of dissimilatory sulfate reduction rates in marine sediment via radioactive 35S tracer. Limnol Oceanogr Methods. 2014;12:196–211.Article 

Google Scholar 
García-Robledo E, Corzo A, Papaspyrou S. A fast and direct spectrophotometric method for the sequential determination of nitrate and nitrite at low concentrations in small volumes. Marine Chem. 2014;162:30–36.Article 

Google Scholar 
Miranda KM, Espey MG, Wink DA. A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide. 2001;5:62–71.Article 
CAS 
PubMed 

Google Scholar 
Holtappels M, Lavik G, Jensen MM, Kuypers MMM Chapter ten – 15N-Labeling Experiments to Dissect the Contributions of Heterotrophic Denitrification and Anammox to Nitrogen Removal in the OMZ Waters of the Ocean. In: Klotz MG, editor. Methods in Enzymology. 486: Academic Press; 2011. p. 223-251.Preisler A, De Beer D, Lichtschlag A, Lavik G, Boetius A, Jørgensen BB. Biological and chemical sulfide oxidation in a Beggiatoa inhabited marine sediment.ISME J. 2007;1:341–353.Article 
CAS 
PubMed 

Google Scholar 
Warembourg FR 5 – Nitrogen fixation in soil and plant systems. In: Knowles R, Blackburn TH, editors. Nitrogen Isotope Techniques. San Diego: Academic Press; 1993. p. 127-156.Orellana LH, Rodriguez-R LM, Konstantinidis KT. ROCker: accurate detection and quantification of target genes in short-read metagenomic data sets by modeling sliding-window bitscores. Nucleic Acids Res. 2016;45:e14–e14.PubMed Central 

Google Scholar 
Menzel P, Ng KL, Krogh A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7:11257.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Anantharaman K, Hausmann B, Jungbluth SP, Kantor RS, Lavy A, Warren LA, et al. Expanded diversity of microbial groups that shape the dissimilatory sulfur cycle. ISME J. 2018;12:1715–28.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Watanabe T, Kojima H, Fukui M. Identity of major sulfur-cycle prokaryotes in freshwater lake ecosystems revealed by a comprehensive phylogenetic study of the dissimilatory adenylylsulfate reductase. Sci Rep. 2016;6:36262.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar  More

Scanning electron microscopy and energy spectrum analysisThe SEM images show the morphological structures of DA and Al-DA before and after adsorption (Fig. 1). DA and Al-DA have disk-like microstructures29 with sur-faces containing both large and small pores, that is, DA and Al-DA have unique multi-level pore structures. The main component of DA and Al-DA is silica, which has a large specific surface area, good thermal stability, and is a natural green material for use as a water treatment agent with a porous structure31. The micrographs show that before adsorption, the DA surface is smooth with a distinct pore structure, whereas modification with aluminium hydroxide makes DA coarse and loose because of the formation of amorphous aluminium hydroxide colloids32. After adsorption, the surface pore structure is covered over for DA and completely covered over for Al-DA, which indicates that F− reacts with Al3+ to form nanoscale precipitates22. The results of the EDS analysis (Fig. 2) show that the content of elemental Al increased from 3.96 to 12.74% after DA was modified with aluminium hydroxide, indicating that Al adhered effectively to the modified DA surface. After adsorption, the content of elemental Al decreased from 3.96 to 1.36% for DA and from 12.74 to 2.03% for Al-DA, which fully confirmed that fluorine preferentially combined with Al to form aluminium precipitates during adsorption, thereby decreasing the Al content.Figure 1SEM images of DA and Al-DA before and after adsorption. (A) Before DA adsorption. (B) After DA adsorption. (C) Before Al-DA adsorption. (D) After Al-DA adsorption.Full size imageFigure 2EDS graphs of DA and Al-DA before and after adsorption. (A) Before DA adsorption. (B) After DA adsorption. (C) Before Al-DA adsorption. (D) After Al-DA adsorption.Full size imageXRD analysisThe surface mineral composition and crystallinity of the materials before and after adsorption were analyzed by XRD (Fig. 3). In the DA and Al-DA patterns, the wide diffraction peaks at approximately 22.0°, 26.0°, and 50.0° mainly correspond to amorphous SiO2, and the diffraction peak at approximately 35° mainly corresponds to amorphous Al2O3, indicating that the material is polycrystalline29. It has been re-ported that amorphous materials may be good adsorbents because of a large specific surface area and numerous active sites33. Many Al(OH)3 peaks and NaCl peaks appear in the XRD pattern of Al-DA, indicating the successful modification of DA by aluminium hydroxide. After adsorption, Na3AlF6 peaks appear in the DA pattern, and Na3AlF6 and AlF3 peaks appear in the Al-DA pattern, whereas the characteristic peaks of NaCl are absent in the Al-DA pattern, which indicates the participation of NaCl in the adsorption process. It has been demonstrated that in the presence of excess sodium fluoride in the reaction solution, the generated aluminium fluoride combines with sodium fluoride to form a NaAlF4 intermediate, which is subsequently converted to cryolite complexes by further adsorption of sodium fluoride34. This result confirms the XRD mapping results.Figure 3XRD patterns of DA and Al-DA before and after adsorption.Full size imageInfrared analysisFigure 4 shows the FTIR spectra of DA and Al-DA before and after adsorption: peaks at 3418, 1635, 1096, 791, and 538 cm−1 appear in the spectrum of DA spectrum before adsorption, and peaks at 3630, 3449, 1637, 1094, 913, 793, and 538 cm−1, appear in the Al-DA spectrum before adsorption. The strong and broad band centered at 3418 cm−1 is due to the stretching vibration of the adsorbed water hydroxyl group (O–H) and the surface hydroxyl group, the vibrational peak at approximately 1635 cm−1 is probably from bound water or the surface hydroxyl group; the peaks at 1096 cm−1 and 538 cm−1 correspond to siloxane groups (Si–O–Si–) and an Al–O absorption band, respectively; and the strong oscillations at 791 cm−1 may be attributed to inorganic Al salts35,36,37. The original absorption peak in the DA spectrum is shifted in the spectrum of DA modified with aluminium hydroxide, confirming the successful modification of DA. The shift of the band at 3418 cm−1 in the DA spectrum to a higher frequency at 3623 cm−1 in the DA spectrum after fluoride absorption is caused by fluoride bonding and has been previously reported38. Another noticeable change in the spectra of DA and Al-DA before and after adsorption is the increase or decrease in the intensity of bending vibrations of specific peaks because the highly electronegative fluoride may have an inductive effect on the respective groups that leads to a blueshift, and the formation of hydrogen bonds leads to a redshift and broadening of the spectral band. The shifts and changes of these peaks indicate the interaction of fluoride with the respective groups29. The new peak at approximately 1170 cm−1 in the spectra of DA and Al-DA with adsorbed fluoride may be due to the formation of Al-F bonds6. The IR spectra show that the formation of a new bonding electronic structure by surface complexation with F− is one of the main mechanisms for the adsorption of F−.Figure 4FTIR spectra of DA and Al-DA before and after adsorption.Full size imageZeta potential analysisThe zeta potential of the material surface plays a very key role in the adsorption process, which reflects the surface charge properties of the material under different pH conditions, and also reflects the surface properties of the material. To obtain the zero charge point of the material, we studied the potential change of the material under different pH values. The results are shown in Fig. 5. In the range of pH 3–11, the zeta potential of the two adsorbents decreased linearly with the increase in pH, and the pHzpc of DA and Al-DA were 9.84 and 10.61, respectively. When pH  More

Study siteOur study site is in an intensive agricultural landscape in NW Spain known as “Tierra de Campos”, which occupies part of three out of nine provinces of Castilla-y-León region (Palencia, Valladolid, and Zamora). This area is considered the main “hot-spot” of tularemia in Spain and Southern Europe16 and is characterized by higher-than-average vole abundances during outbreaks17.Iberian hare abundance indexYearly occurrence of vole outbreaks in NW Spain between 1996 and 2020 (i.e., 1997, 1998, 2007, 2008, 2014, 2019) were identified based on reports in the news (historical reconstruction18) and more recently (from 2009 onward) using common vole abundance indices obtained from live-trapping monitoring (i.e.4,19).To study the Iberian hare population trends we used regional hunting statistics available from the regional government (Junta de Castilla-y-León, CAZDATA Project, https://medioambiente.jcyl.es/web/es/caza-pesca/cazdata-banco-datos-actividad.html [Cited 2022 Sep 23]), which included hunting records as well as the number of hunting licences from 1974 to 2020. We used the number of hunted hares divided by the number of hunting licences each year as an abundance index for hares in “Tierra de Campos” (compiling data from the provinces of Palencia, Zamora and Valladolid). CAZDATA Project is an initiative proposed by the Hunting Federation of Castilla y León, which has the support of the regional government and, more importantly, the commitment of almost 60% of the hunting societies in the community to implement a system for monitoring hunting activity. Since this information is gathered by hunters for the benefit of the hunting activity, we are confidence on its reliability to carry out the present study.
Francisella tularensis prevalence in Iberian haresWe compiled data on F. tularensis prevalence in Iberian hares from 2007 to 2016 using previously published information from a passive surveillance program carried out by the “Regional Network of Epidemiological Surveillance” (Red de Vigilancia Epidemiológica de la Dirección General de Salud Pública) of Castilla-y-León region20. This provided us with information on hare tularemia prevalence (amount of positives/number of screened individual) each year within the three provinces from “Tierra de Campos”.Statistical analysesTo study Iberian hare population trends, we calculated an index of yearly hare population instantaneous growth rate (PGR) using the hunting bag data (hare abundance index) from 1996 to 2020. Hare PGR was calculated as follows:$$PGR= lnleft(frac{{AI}_{t}}{{AI}_{t+1}}right)$$where ln stands for natural logarithm, AIt is Iberian hare abundance index on year t. and AIt+1 is the Iberian hare abundance index on year t + 1. PGRs were estimated yearly from 1996 to 2019. This dependent variable was fitted to a Generalized Linear Mixed Model using the glmmTMB function (GLMMTMB, package glmmTMB21) and a gaussian family distribution and identity link function. The categorical variable vole outbreak year (i.e., with two levels: years with (1) or without vole outbreak (0), hereafter “Vole”) and “Province” (i.e., with three levels: Palencia, Valladolid and Zamora), and their interaction were used as explanatory variables. “Year” of sampling was included as a random factor (i.e., 1996–2019). Significance of the fixed effects in the models was calculated with Type II tests using the function Anova in the car package22. We previously checked the model for overdispersion and distribution fitting using function simulateResiduals (package DHARMa23, simulations = 999). The variable PGR expresses the change between year t and t + 1. We included AI at t as a covariate in the model, in order to take into account density-dependence in hare PGR (the extent to which the abundance changes in between year t and t + 1 depends on the abundance during year t). For this to make biological sense, we rescaled the covariable AI so that it has mean equal to zero. Thus, the effect of the other predictor variables in the model (i.e., “Vole” and “Province”) was interpreted as the effect that these variables have on PGR when the abundance value is at 0. Thus, the effect of “Vole” and “Province” on PGR will be obtained by the mean value of abundance.We assessed the effect of vole outbreak years on the Iberian hare’s population PGR by running a multiple Pearson correlation (function ggscatter) between PGR and AI, considering both, PGR for all the years of the study period (i.e., 1996–2019) and only those years where vole outbreaks were detected (i.e., 1997, 1998, 2007, 2008, 2014, 2019).Finally, we tested for difference in the prevalence of F. tularensis on Iberian hare’s during years with or without vole outbreaks using a GLMMTMB21 with a binomial family distribution and a logit link function, where prevalence of F. tularensis in hares was the dependent variable, and “Vole” outbreak years and “Province” (i.e. Palencia, Valladolid and Zamora) were the responses variables. In this case, the variable “Vole” outbreak years included three levels (i.e. 0 = no vole outbreak, 1 = vole outbreak year, 2 = one year after vole outbreak), to assess if F. tularensis prevalence in hare also persist one year after a vole outbreak. “Year” of sampling was included as a random factor (i.e., 2007–2016). Due to the limited sample size, we did not include the interaction between “Vole” and “Province” to not overfit the model. We also previously checked the model for overdispersion and distribution fitting using function simulateResiduals (package DHARMa23, simulations = 999). All analysis were carried out using the R statistical computing environment24. More

Visser, M. E. & Holleman, L. J. Warmer springs disrupt the synchrony of oak and winter moth phenology. Proc. Biol. Sci. 268, 289–294 (2001).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Duchenne, F. et al. Phenological shifts alter the seasonal structure of pollinator assemblages in Europe. Nat. Ecol. Evol. 4, 115–121 (2020).Article 
CAS 
PubMed 

Google Scholar 
Bartomeus, I. et al. Climate-associated phenological advances in bee pollinators and bee-pollinated plants. Proc. Natl. Acad. Sci. U. S. A. 108, 20645–20649 (2011).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lehmann, P. et al. Complex responses of global insect pests to climate warming. Front. Ecol. Environ. 18, 141–150 (2020).Article 

Google Scholar 
Organization, W. H. et al. Global vector control response 2017–2030. In Global Vector Control Response 2017–2030 (2017).Reisen, W. K. et al. Effects of warm winter temperature on the abundance and gonotrophic activity of Culex (Diptera: Culicidae) in California. J. Med. Entomol. 47, 230–237 (2010).Article 
PubMed 

Google Scholar 
Denlinger, D. L. & Armbruster, P. A. Mosquito diapause. Annu. Rev. Entomol. 59, 73–93 (2014).Article 
CAS 
PubMed 

Google Scholar 
Mullen, G. R. & Durden, L. A. Medical and Veterinary Entomology. (Academic Press, 2009).Messenger, P. S. Bioclimatic studies with insects. Annu. Rev. Entomol. 4, 183–206 (1959).Article 

Google Scholar 
Hard, J. J., Bradshaw, W. E. & Holzapfel, C. M. The genetic basis of photoperiodism and its evolutionary divergence among populations of the pitcher-plant mosquito, Wyeomyia smithii. Am. Nat. 142, 457–473 (1993).Article 
CAS 
PubMed 

Google Scholar 
Danilevskii, A. S. et al. Photoperiodism and seasonal development of insects. In Photoperiodism and Seasonal Development of Insects (1965).Nietschke, B. S., Magarey, R. D., Borchert, D. M., Calvin, D. D. & Jones, E. A developmental database to support insect phenology models. Crop Prot. 26, 1444–1448 (2007).Article 

Google Scholar 
Vinogradova, E. B. Diapause in aquatic insects, with emphasis on mosquitoes. In Diapause in Aquatic Invertebrates Theory and Human Use (eds. Alekseev, V. R., de Stasio, B. T. & Gilbert, J. J.) 83–113 (Springer Netherlands, 2007).Chown, S. L. & Terblanche, J. S. Physiological diversity in insects: Ecological and evolutionary contexts. Adv. Insect Phys. 33, 50–152 (2006).Article 

Google Scholar 
Belitz, M. W. et al. Climate drivers of adult insect activity are conditioned by life history traits. Ecol. Lett. 24, 2687–2699 (2021).Article 
PubMed 

Google Scholar 
Townroe, S. & Callaghan, A. British container breeding mosquitoes: The impact of urbanisation and climate change on community composition and phenology. PLoS ONE 9, e95325 (2014).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Westby, K. M., Adalsteinsson, S. A., Biro, E. G., Beckermann, A. J. & Medley, K. A. Aedes albopictus populations and larval habitat characteristics across the landscape: Significant differences exist between urban and rural land use types. Insects 12, 196 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Li, D., Stucky, B. J., Deck, J., Baiser, B. & Guralnick, R. P. The effect of urbanization on plant phenology depends on regional temperature. Nat. Ecol. Evol. 3, 1661–1667 (2019).Article 
PubMed 

Google Scholar 
Diniz, D. F. A., de Albuquerque, C. M. R., Oliva, L. O., de Melo-Santos, M. A. V. & Ayres, C. F. J. Diapause and quiescence: Dormancy mechanisms that contribute to the geographical expansion of mosquitoes and their evolutionary success. Parasit. Vectors 10, 310 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Rund, S. S. C., Moise, I. K., Beier, J. C. & Martinez, M. E. Rescuing troves of hidden ecological data to tackle emerging mosquito-borne diseases. J. Am. Mosq. Control Assoc. 35, 75–83 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Dewitz, J., U.S. Geological Survey. National Land Cover Database (NLCD) 2019 Products (ver. 2. 0, June 2021): U.S. Geological Survey data release (ver. 2. 0, June 2021): U.S. Geological Survey data release.Woodring, J. et al. Diapause, transovarial transmission, and filial infection rates in geographic strains of La Crosse virus-infected Aedes triseriatus. Am. J. Trop. Med. Hyg. 58, 587–588 (1998).Article 
CAS 
PubMed 

Google Scholar 
Ellwood, E. R. et al. Disentangling the paradox of insect phenology: Are temporal trends reflecting the response to warming?. Oecologia 168, 1161–1171 (2012).Article 
ADS 
PubMed 

Google Scholar 
Degaetano, A. T. Meteorological effects on adult mosquito (Culex) populations in metropolitan New Jersey. Int. J. Biometeorol. 49, 345–353 (2005).Article 
ADS 
PubMed 

Google Scholar 
Su, T., Webb, J. P., Meyer, R. P. & Mulla, M. S. Spatial and temporal distribution of mosquitoes in underground storm drain systems in Orange County, California. J. Vector Ecol. 28, 79–89 (2003).PubMed 

Google Scholar 
Harbison, J. E., Henry, M., Xamplas, C. & Dugas, L. R. Evaluation of Culex pipiens populations in a residential area with a high density of catch basins in a suburb of Chicago, Illinois. J. Am. Mosq. Control Assoc. 30, 228–230 (2014).Article 
PubMed 

Google Scholar 
Wang, X. et al. Impact of underground storm drain systems on larval ecology of Culex and Aedes species in urban environments of Southern California. Sci. Rep. 11, 12667 (2021).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Geery, P. R. & Holub, R. E. Seasonal abundance and control of Culex spp. in catch basins in Illinois. J. Am. Mosq. Control Assoc. 5, 537–540 (1989).CAS 
PubMed 

Google Scholar 
Nelms, B. M., Macedo, P. A., Kothera, L., Savage, H. M. & Reisen, W. K. Overwintering biology of Culex (Diptera: Culicidae) mosquitoes in the Sacramento Valley of California. J. Med. Entomol. 50, 773–790 (2013).Article 
PubMed 

Google Scholar 
Becker, N. et al. Mosquitoes and Their Control (Springer, 2010).Book 

Google Scholar 
Bartlett-Healy, K., Crans, W. & Gaugler, R. Temporal and spatial synchrony of Culex territans (Diptera: Culicidae) with their amphibian hosts. J. Med. Entomol. 45, 1031–1038 (2008).Article 
PubMed 

Google Scholar 
Reeves, L. E. et al. Identification of Uranotaenia sapphirina as a specialist of annelids broadens known mosquito host use patterns. Commun. Biol. 1, 92 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Burkett-Cadena, N. D. et al. Host reproductive phenology drives seasonal patterns of host use in mosquitoes. PLoS One 6, e17681 (2011).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
van Strien, A. J., Plantenga, W. F., Soldaat, L. L., van Swaay, C. A. M. & Wallisdevries, M. F. Bias in phenology assessments based on first appearance data of butterflies. Oecologia 156, 227–235 (2008).Article 
ADS 
PubMed 

Google Scholar 
Beck, H. E. et al. Present and future Köppen–Geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Giraldo-Calderón, G. I. et al. Vectorbase.org updates: Bioinformatic resources for invertebrate vectors of human pathogens and related organisms. Curr. Opin. Insect Sci. 50, 100860 (2022).Article 
PubMed 

Google Scholar 
Wickham, François, Henry & Müller. dplyr: A grammar of data manipulation. R package version 0.4.Hesselbarth, M. H. K., Sciaini, M., With, K. A., Wiegand, K. & Nowosad, J. landscapemetrics: An open-source R tool to calculate landscape metrics. Ecography 42, 1648–1657 (2019).Article 

Google Scholar 
Abatzoglou, J. T. Development of gridded surface meteorological data for ecological applications and modelling. Int. J. Climatol. 33, 121–131 (2013).Article 

Google Scholar 
Johnson, M. climateR: climateR. R package version 0.1.0. https://github.com/mikejohnson51/climateR.Wickham, H. et al. Welcome to the tidyverse. J. Open Source Softw. 4, 1686 (2019).Article 
ADS 

Google Scholar 
Schmucki, R., Harrower, C. A. & Dennis, E. B.. rbms: Computing generalised abundance indices for butterfly monitoring count data. R package version.Belitz, M. W., Larsen, E. A., Shirey, V., Li, D. & Guralnick, R. P. Phenological research based on natural history collections: Practical guidelines and a lepidopteran case study. Funct. Ecol. https://doi.org/10.1111/1365-2435.14173 (2022).Article 

Google Scholar 
Larsen, E. A., Belitz, M. W., Guralnick, R. P. & Ries, L. Consistent trait-temperature interactions drive butterfly phenology in both incidental and survey data. Sci. Rep. 12, 13370 (2022).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models Usinglme4. J. Stat. Softw. https://doi.org/10.18637/jss.v067.i01 (2015).Article 

Google Scholar 
Fox, J. & Weisberg, S. An R Companion to Applied Regression. (SAGE Publications, 2018).Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).Article 

Google Scholar 
Lüdecke, D., Ben-Shachar, M., Patil, I., Waggoner, P. & Makowski, D. Performance: An R package for assessment, comparison and testing of statistical models. J. Open Source Softw. 6, 3139 (2021).Article 
ADS 

Google Scholar 
Nakagawa, S., Johnson, P. C. D. & Schielzeth, H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R. Soc. Interface 14, 20170213 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Lüdecke, D. sjPlot: Data visualization for statistics in social science. R package version. More