Philippa Kaur

Polis, G. A., Anderson, W. B. & Holt, R. D. Toward an integration of landscape and food web ecology: the dynamics of spatially subsidized food webs. Annu. Rev. Ecol. Evol. Syst. 28, 289–316 (1997).Article 

Google Scholar 
Doughty, C. E. et al. Global nutrient transport in a world of giants. Proc. Natl Acad. Sci. USA 113, 868–873 (2016).Article 
CAS 
PubMed 

Google Scholar 
Burpee, B. T. & Saros, J. E. Cross-ecosystem nutrient subsidies in Arctic and alpine lakes: implications of global change for remote lakes. Environ. Sci. 22, 1166–1189 (2020).CAS 

Google Scholar 
Gallardo, B., Clavero, M., Sánchez, M. I. & Vilà, M. Global ecological impacts of invasive species in aquatic ecosystems. Glob. Change Biol. 22, 151–163 (2016).Article 

Google Scholar 
Justino, D. G., Maruyama, P. K. & Oliveira, P. E. Floral resource availability and hummingbird territorial behaviour on a Neotropical savanna shrub. J. Ornithol. 153, 189–197 (2012).Article 

Google Scholar 
Van Overveld, T. et al. Food predictability and social status drive individual resource specializations in a territorial vulture. Sci. Rep. 8, 15155 (2018).Gunn, R. L., Hartley, I. R., Algar, A. C., Nadiarti, N. & Keith, S. A. Variation in the behaviour of an obligate corallivore is influenced by resource availability. Behav. Ecol. Sociobiol. https://doi.org/10.1007/s00265-022-03132-6 (2022).Keith, S. A. et al. Synchronous behavioural shifts in reef fishes linked to mass coral bleaching. Nat. Clim. Change 8, 986–991 (2018).Article 

Google Scholar 
Davies, N. B. & Hartley, I. R. Food patchiness, territory overlap and social systems: an experiment with dunnocks Prunella modularis. J. Anim. Ecol. 65, 837–846 (1996).Article 

Google Scholar 
Cahill, A. E. et al. How does climate change cause extinction? Proc. R. Soc. B 280, 20121890 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Delarue, P. E. M., Kerr, S. E. & Rymer, T. L. Habitat complexity, environmental change and personality: a tropical perspective. Behav. Process. 120, 101–110 (2015).Stimson, J. The role of the territory in the ecology of the intertidal limpet Lottia gigantea (Gray). Ecology 54, 1020–1030 (1973).Article 

Google Scholar 
Sells, S. N. & Mitchell, M. S. The economics of territory selection. Ecol. Modell. 438, 109329 (2020).Article 

Google Scholar 
Graf, P. M., Mayer, M., Zedrosser, A., Hackländer, K. & Rosell, F. Territory size and age explain movement patterns in the Eurasian beaver. Mamm. Biol. 81, 587–594 (2016).Article 

Google Scholar 
Simon, C. The influence of food abundance on territory size in the Iguanid lizard Sceloporus jarrovi. Ecology 56, 993–998 (1975).Article 

Google Scholar 
Ippi, S., Cerón, G., Alvarez, L. M., Aráoz, R. & Blendinger, P. G. Relationships among territory size, body size, and food availability in a specialist river duck. Emu 118, 293–303 (2018).Article 

Google Scholar 
Berumen, M. L. & Pratchett, M. S. Effects of resource availability on the competitive behaviour of butterflyfishes (Chaetodontidae). In Proc. 10th International Coral Reef Symposium 644–650 (ReefBase, 2006); http://reefbase.org/resource_center/publication/icrs.aspx?icrs=ICRS10Brown, J. L. The evolution of diversity in avian territorial systems. Wilson Bull. 76, 160–169 (1964).
Google Scholar 
Peiman, K. S. & Robinson, B. W. Ecology and evolution of resource-related heterospecific aggression. Q. Rev. Biol. 85, 133–158 (2010).Article 
PubMed 

Google Scholar 
Grant, J. W. A., Girard, I. L., Breau, C. & Weir, L. K. Influence of food abundance on competitive aggression in juvenile convict cichlids. Anim. Behav. 63, 323–330 (2002).Article 

Google Scholar 
Duda, M. P. et al. Long-term changes in terrestrial vegetation linked to shifts in a colonial seabird population. Ecosystems 23, 1643–1656 (2020).Article 
CAS 

Google Scholar 
Graham, N. A. J. et al. Seabirds enhance coral reef productivity and functioning in the absence of invasive rats. Nature 559, 250–253 (2018).Article 
CAS 
PubMed 

Google Scholar 
Jones, H. P. et al. Severity of the effects of invasive rats on seabirds: a global review. Conserv. Biol. 22, 16–26 (2008).Article 
PubMed 

Google Scholar 
Honig, S. E. & Mahoney, B. Evidence of seabird guano enrichment on a coral reef in Oahu, Hawaii. Mar. Biol. 163, 22 (2016).Benkwitt, C. E., Gunn, R. L., Le Corre, M., Carr, P. & Graham, N. A. J. Rat eradication restores nutrient subsidies from seabirds across terrestrial and marine ecosystems. Curr. Biol. 31, 2704–2711.e4 (2021).Article 
CAS 
PubMed 

Google Scholar 
Savage, C. Seabird nutrients are assimilated by corals and enhance coral growth rates. Sci. Rep. 9, 4284 (2019).Benkwitt, C. E., Wilson, S. K. & Graham, N. A. J. Seabird nutrient subsidies alter patterns of algal abundance and fish biomass on coral reefs following a bleaching event. Glob. Change Biol. 25, 2619–2632 (2019).Article 

Google Scholar 
Benkwitt, C. E., Taylor, B. M., Meekan, M. G. & Graham, N. A. J. Natural nutrient subsidies alter demographic rates in a functionally important coral-reef fish. Sci. Rep. 11, 12575 (2021).Benkwitt, C. E., Wilson, S. K. & Graham, N. A. J. Biodiversity increases ecosystem functions despite multiple stressors on coral reefs. Nat. Ecol. Evol. 4, 919–926 (2020).Article 
PubMed 

Google Scholar 
Robles, H. & Martin, K. Resource quantity and quality determine the inter-specific associations between ecosystem engineers and resource users in a cavity-nest web. PLoS ONE 8, e74694 (2013).Catano, L. B., Gunn, B. K., Kelley, M. C. & Burkepile, D. E. Predation risk, resource quality, and reef structural complexity shape territoriality in a coral reef herbivore. PLoS ONE 10, e0118764 (2015).Wilcox, K. A., Wagner, M. A. & Reynolds, J. D. Salmon subsidies predict territory size and habitat selection of an avian insectivore. PLoS ONE 16, e0254314 (2021).Frost, S. K. & Frost, P. G. H. Territoriality and changes in resource use by sunbirds at Leonotis leonurus (Labiatae). Oecologia 45, 109–116 (1980).Maynard Smith, J. Evolution and the Theory of Games (Cambridge Univ. Press, 1982).Book 

Google Scholar 
Dochtermann, N. A., Schwab, T., Anderson Berdal, M., Dalos, J. & Royauté, R. The heritability of behavior: a meta-analysis. J. Hered. 110, 403–410 (2019).Article 
PubMed 

Google Scholar 
Sheppard, C. R. C. et al. Reefs and islands of the Chagos Archipelago, Indian Ocean: why it is the world’s largest no-take marine protected area. Aquat. Conserv. 22, 232–261 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Soeparno, Y. N., Shibuno, T. & Yamaoka, K. Relationship between pelagic larval duration and abundance of tropical fishes on temperate coasts of Japan. J. Fish. Biol. 80, 346–357 (2012).Article 
CAS 
PubMed 

Google Scholar 
Green, A. L. et al. Larval dispersal and movement patterns of coral reef fishes, and implications for marine reserve network design. Biol. Rev. 90, 1215–1247 (2015).Article 
PubMed 

Google Scholar 
Dall, S. R. X., Houston, A. I. & McNamara, J. M. The behavioural ecology of personality: consistent individual differences from an adaptive perspective. Ecol. Lett. 7, 734–739 (2004).Article 

Google Scholar 
Klumpp, D., McKinnon, D. & Daniel, P. Damselfish territories: zones of high productivity on coral reefs. Mar. Ecol. Prog. Ser. 40, 41–51 (1987).Article 

Google Scholar 
Carr, P. et al. Status and phenology of breeding seabirds and a review of important bird and biodiversity areas in the British Indian Ocean Territory. Bird Conserv. Int. 31, 14–34 (2020).Article 

Google Scholar 
Hoey, A. S. & Bellwood, D. R. Damselfish territories as a refuge for macroalgae on coral reefs. Coral Reefs 29, 107–118 (2010).Article 

Google Scholar 
Samways, M. J. Breakdown of butterflyfish (Chaetodontidae) territories associated with the onset of a mass coral bleaching event. Aquat. Conserv. 15, 101–107 (2005).Article 

Google Scholar 
Morgan, I. E. & Kramer, D. L. Determinants of social organization in a coral reef fish, the blue tang, Acanthurus coeruleus. Environ. Biol. Fishes 72, 443–453 (2005).Article 

Google Scholar 
Ceccarelli, D. M. Modification of benthic communities by territorial damselfish: a multi-species comparison. Coral Reefs 26, 853–866 (2007).Article 

Google Scholar 
Gochfeld, D. J. Territorial damselfishes facilitate survival of corals by providing an associational defense against predators. Mar. Ecol. Prog. Ser. 398, 137–148 (2010).Article 

Google Scholar 
Gordon, T. A. C., Cowburn, B. & Sluka, R. D. Defended territories of an aggressive damselfish contain lower juvenile coral density than adjacent non-defended areas on Kenyan lagoon patch reefs. Coral Reefs 34, 13–16 (2015).Article 

Google Scholar 
Hays, G. C. et al. A review of a decade of lessons from one of the world’s largest MPAs: conservation gains and key challenges. Mar. Biol. 167, 159–167 (2020).Article 

Google Scholar 
Nanninga, G. B., Côté, I. M., Beldade, R. & Mills, S. C. Behavioural acclimation to cameras and observers in coral reef fishes. Ethology 123, 705–711 (2017).Article 

Google Scholar 
Polunin, N. V. C. & Klumpp, D. W. Ecological correlates of foraging periodicity in herbivorous reef fishes of the Coral Sea. J. Exp. Mar. Biol. Ecol. 126, 1–20 (1989).Article 

Google Scholar 
Friard, O. & Gamba, M. BORIS: a free, versatile open-source event-logging software for video/audio coding and live observations. Methods Ecol. Evol. 7, 1325–1330 (2016).Article 

Google Scholar 
Paola, V. D., Vullioud, P., Demarta, L., Alwany, M. A. & Ros, A. F. H. Factors affecting interspecific aggression in a year-round territorial species, the jewel damselfish. Ethology 118, 721–732 (2012).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).Bürkner, P. brms: an R package for Bayesian multilevel models using Stan. J. Stat. Soft. 80, 1–28 (2017).Article 

Google Scholar 
RStan: the R interface to Stan. R package version 2.21.5 (Stan Development Team, 2022).Hadfield, J. D. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Soft. 33, 1–22 (2010).Article 

Google Scholar 
Gelman, A. & Rubin, D. B. Inference from iterative simulation using multiple sequences. Stat. Sci. 7, 457–472 (1992).Article 

Google Scholar 
Vehtari, A., Gelman, A., Simpson, D., Carpenter, B. & Burkner, P. C. Rank-normalization, folding, and localization: an improved (formula presented) for assessing convergence of MCMC (with Discussion). Bayesian Anal. 16, 667–718 (2021).Vehtari, A., Gelman, A. & Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 27, 1413–1432 (2017).Article 

Google Scholar  More

Potapov, P. et al. Global maps of cropland extent and change show accelerated cropland expansion in the twenty-first century. Nat. Food 3, 19–28 (2022).Article 

Google Scholar 
Land Use Statistics and Indicators. Global, Regional and Country Trends 2000–2020 FAOSTAT Analytical Brief Series No 48 https://www.fao.org/food-agriculture-statistics/data-release/data-release-detail/en/c/1599856/ (FAO, 2022).FAO. Land Statistics. Global, Regional and Country Trends, 1990–2018 FAOSTAT Analytical Brief Series No. 15 https://www.fao.org/3/cb2860en/cb2860en.pdf (FAO, 2021).Summary for policymakers in: Special Report on Climate Change and Land (eds Shukla, P. R. et al.) https://www.ipcc.ch/site/assets/uploads/sites/4/2020/02/SPM_Updated-Jan20.pdf (WMO, in the press).Sustainable Development Goals Indicator 2.4.1 (FAO, accessed); https://www.fao.org/sustainable-development-goals/indicators/241/en/Eggleston, H. S., Buendia, L., Miwa, K., Ngara, T. & Tanabe, K. 2006 IPCC Guidelines for National Greenhouse Gas Inventories (IGES, 2006).Grassi, G. et al. Carbon fluxes from land 2000–2020: bringing clarity on countries’ reporting. Earth Syst. Sci. Data 14, 4643–4666 (2022).Article 
ADS 

Google Scholar 
Tubiello, F. N. et al. Measuring Progress Towards Sustainable Agriculture FAO Statistical Working Papers Series No. 21–24 https://www.fao.org/3/cb4549en/cb4549en.pdf (FAO, 2021).Conchedda, G. & Tubiello, F. N. Drainage of organic soils and GHG emissions: validation with country data. Earth Syst. Sci. Data 12, 3113–3137 (2020).Article 
ADS 

Google Scholar 
Hanson, C., Mazur, E., Stolle, F., Davis, C. & Searchinger, T. 5 takeaways on cropland expansion and what it means for people and the planet. WRI Insights https://www.wri.org/insights/cropland-expansion-impacts-people-planet (2022).Potapov, P. et al. The Global 2000–-2020 land cover and land use change dataset derived from the Landsat archive: first results. Front. Remote Sens. 3, 856903 (2022).Article 

Google Scholar 
Hansen, M. C. et al. Global land use extent and dispersion within natural land cover using Landsat data. Environ. Res. Lett. 17, 034050 (2022).Article 
ADS 

Google Scholar 
Tubiello, F. N. et al. Carbon emissions and removals from forests: new estimates, 1990–2020. Earth Syst. Sci. Data. 13, 1681–1691 (2021).Article 
ADS 

Google Scholar  More

Berner, L. T. et al. Summer warming explains widespread but not uniform greening in the Arctic tundra biome. Nat. Commun. 11, 4621 (2020).Article 
ADS 
CAS 

Google Scholar 
Elmendorf, S. C. et al. Plot-scale evidence of tundra vegetation change and links to recent summer warming. Nat. Clim. Change 2, 453–457 (2012).Article 
ADS 

Google Scholar 
Overland, J. E. et al. Surface air temperature. In Arctic Report Card: Update for 2019 (eds Richter-Menge, J. et al.) (U.S. National Park Service, 2020).
Google Scholar 
Post, E., Steinman, B. A. & Mann, M. E. Acceleration of phenological advance and warming with latitude over the past century. Sci. Rep. 8, 3927 (2018).Article 
ADS 

Google Scholar 
Diepstraten, R. A. E., Jessen, T. D., Fauvelle, C. M. D. & Musiani, M. M. Does climate change and plant phenology research neglect the Arctic tundra?. Ecosphere 9, e02362 (2018).Article 

Google Scholar 
Flynn, D. F. B. & Wolkovich, E. M. Temperature and photoperiod drive spring phenology across all species in a temperate forest community. New Phytol. 219, 1353–1362 (2018).Article 
CAS 

Google Scholar 
Billings, W. D. & Bliss, L. C. An alpine snowbank environment and its effects on vegetation, plant development, and productivity. Ecology 40, 388–397 (1959).Article 

Google Scholar 
Billings, W. D. & Mooney, H. A. The ecology of arctic and alpine plants. Biol. Rev. 43, 481–529 (1968).Article 

Google Scholar 
Sørensen, T. Temperature relations and phenology of the northeast Greenland flowering plants. Meddr Gronland 1–305 (1941).Barrett, R. T. & Hollister, R. D. Arctic plants are capable of sustained responses to long-term warming. Polar Res. 35, 25405 (2016).Article 

Google Scholar 
Julitta, T. et al. Using digital camera images to analyse snowmelt and phenology of a subalpine grassland. Agric. For. Meteorol. 198–199, 116–125 (2014).Article 
ADS 

Google Scholar 
Petraglia, A. et al. Responses of flowering phenology of snowbed plants to an experimentally imposed extreme advanced snowmelt. Plant Ecol. 215, 759–768 (2014).Article 

Google Scholar 
Semenchuk, P. R. et al. High Arctic plant phenology is determined by snowmelt patterns but duration of phenological periods is fixed: An example of periodicity. Environ. Res. Lett. 11, 125006 (2016).Article 
ADS 

Google Scholar 
Hollister, R. D., Webber, P. J. & Bay, C. Plant response to temperature in northern Alaska: Implications for predicting vegetation change. Ecology 86, 1562–1570 (2005).Article 

Google Scholar 
Oberbauer, S. et al. Phenological response of tundra plants to background climate variation tested using the International Tundra Experiment. Philos. Trans. R. Soc. B Biol. Sci. 368, 20120481 (2013).Article 
CAS 

Google Scholar 
Tieszen, L. L. Photosynthesis in the principal Barrow, Alaska, species: A summary of field and laboratory responses. In Vegetation and Production Ecology of an Alaskan Arctic Tundra (ed. Tieszen, L. L.) 241–268 (Springer, 1978).Chapter 

Google Scholar 
Körner, Ch. CO2 exchange in the alpine sedge Carex curvula as influenced by canopy structure, light and temperature. Oecologia 53, 98–104 (1982).Article 
ADS 

Google Scholar 
Tieszen, L. L. Photosynthesis and respiration in arctic tundra grasses: Field light intensity and temperature responses. Arct. Alp. Res. 5, 239–251 (1973).Article 
CAS 

Google Scholar 
Huang, M. et al. Air temperature optima of vegetation productivity across global biomes. Nat. Ecol. Evol. 3, 772–779 (2019).Article 

Google Scholar 
Marchand, F. L., Mertens, S., Kockelbergh, F., Beyens, L. & Nijs, I. Performance of high arctic tundra plants improved during but deteriorated after exposure to a simulated extreme temperature event. Glob. Change Biol. 11, 2078–2089 (2005).Article 
ADS 

Google Scholar 
Yan, W. An equation for modelling the temperature response of plants using only the cardinal temperatures. Ann. Bot. 84, 607–614 (1999).Article 

Google Scholar 
Zhou, G. & Wang, Q. A new nonlinear method for calculating growing degree days. Sci. Rep. 8, 10149 (2018).Article 
ADS 

Google Scholar 
Kramer, K. Selecting a model to predict the onset of growth of Fagus sylvatica. J. Appl. Ecol. 31, 172 (1994).Article 

Google Scholar 
Nakano, Y., Higuchi, Y., Sumitomo, K. & Hisamatsu, T. Flowering retardation by high temperature in chrysanthemums: Involvement of FLOWERING LOCUS T-like 3 gene repression. J. Exp. Bot. 64, 909–920 (2013).Article 
CAS 

Google Scholar 
del Olmo, I., Poza-Viejo, L., Piñeiro, M., Jarillo, J. A. & Crevillén, P. High ambient temperature leads to reduced FT expression and delayed flowering in Brassica rapa via a mechanism associated with H2A.Z dynamics. Plant J. 100, 343–356 (2019).Article 

Google Scholar 
Wolkovich, E. M. et al. Warming experiments underpredict plant phenological responses to climate change. Nature 485, 494 (2012).Article 
ADS 
CAS 

Google Scholar 
Hollister, R. D. et al. A review of open top chamber (OTC) performance across the ITEX Network. Arct. Sci. https://doi.org/10.1139/AS-2022-0030 (2022).Article 

Google Scholar 
Bütikofer, L. et al. The problem of scale in predicting biological responses to climate. Glob. Change Biol. 26, 6657–6666 (2020).Article 
ADS 

Google Scholar 
Gu, S. Growing degree hours—A simple, accurate, and precise protocol to approximate growing heat summation for grapevines. Int. J. Biometeorol. 60, 1123–1134 (2016).Article 
ADS 
CAS 

Google Scholar 
Roltsch, W. J., Zalom, F. G., Strawn, A. J., Strand, J. F. & Pitcairn, M. J. Evaluation of several degree-day estimation methods in California climates. Int. J. Biometeorol. 42, 169–176 (1999).Article 
ADS 

Google Scholar 
Richardson, A. D. et al. Ecosystem warming extends vegetation activity but heightens vulnerability to cold temperatures. Nature 560, 368–371 (2018).Article 
ADS 
CAS 

Google Scholar 
Ettinger, A. K., Buonaiuto, D. M., Chamberlain, C. J., Morales-Castilla, I. & Wolkovich, E. M. Spatial and temporal shifts in photoperiod with climate change. New Phytol. 230, 462–474 (2021).Article 
CAS 

Google Scholar 
Seyednasrollah, B., Swenson, J. J., Domec, J.-C. & Clark, J. S. Leaf phenology paradox: Why warming matters most where it is already warm. Remote Sens. Environ. 209, 446–455 (2018).Article 
ADS 

Google Scholar 
Breshears, D. D. et al. Underappreciated plant vulnerabilities to heat waves. New Phytol. 231, 32–39 (2021).Article 

Google Scholar 
Chaudhry, S. & Sidhu, G. P. S. Climate change regulated abiotic stress mechanisms in plants: A comprehensive review. Plant Cell Rep. 41, 1–31 (2022).Article 
CAS 

Google Scholar 
Sun, X. et al. Global diurnal temperature range (DTR) changes since 1901. Clim. Dyn. 52, 3343–3356 (2019).Article 

Google Scholar 
Ballinger, T. J. NOAA Arctic Report Card 2021: Surface Air Temperature. https://doi.org/10.25923/53XD-9K68 (2021).Jagadish, S. V. K., Way, D. A. & Sharkey, T. D. Plant heat stress: Concepts directing future research. Plant Cell Environ. 44, 1992–2005 (2021).Article 
CAS 

Google Scholar 
Gilmore, E. C. Jr. & Rogers, J. S. Heat units as a method of measuring maturity in corn. Agron. J. 50, 611–615 (1958).Article 

Google Scholar 
Sánchez, B., Rasmussen, A. & Porter, J. R. Temperatures and the growth and development of maize and rice: A review. Glob. Change Biol. 20, 408–417 (2014).Article 
ADS 

Google Scholar 
Molitor, D., Junk, J., Evers, D., Hoffmann, L. & Beyer, M. A high-resolution cumulative degree day-based model to simulate phenological development of grapevine. Am. J. Enol. Vitic. 65, 72–80 (2014).Article 

Google Scholar 
CaraDonna, P. J., Iler, A. M. & Inouye, D. W. Shifts in flowering phenology reshape a subalpine plant community. Proc. Natl. Acad. Sci. 111, 4916–4921 (2014).Article 
ADS 
CAS 

Google Scholar 
Inouye, B. D., Ehrlén, J. & Underwood, N. Phenology as a process rather than an event: From individual reaction norms to community metrics. Ecol. Monogr. 89, e01352 (2019).Article 

Google Scholar 
Miles, W. T. S. et al. Quantifying full phenological event distributions reveals simultaneous advances, temporal stability and delays in spring and autumn migration timing in long-distance migratory birds. Glob. Change Biol. 23, 1400–1414 (2017).Article 
ADS 

Google Scholar 
Moussus, J.-P., Julliard, R. & Jiguet, F. Featuring 10 phenological estimators using simulated data. Methods Ecol. Evol. 1, 140–150 (2010).Article 

Google Scholar 
Dowle, M. & Srinivasan, A. data.table: Extension of ‘data.frame’ (2019).Auguie, B. egg: Extensions for ‘ggplot2’: Custom Geom, Custom Themes, Plot Alignment, Labelled Panels, Symmetric Scales, and Fixed Panel Size (2019).Wood, S. & Scheipl, F. gamm4: Generalized Additive Mixed Models using ‘mgcv’ and ‘lme4’ (2020).Auguie, B. gridExtra: Miscellaneous Functions for ‘Grid’ Graphics (2017).Hamner, B. & Frasco, M. Metrics: Evaluation Metrics for Machine Learning (2018).Gilli, M., Maringer, D. & Schumann, E. Numerical Methods and Optimization in Finance (Elsevier/Academic Press, 2019).MATH 

Google Scholar 
Garnier, S. viridis: Default Color Maps from ‘matplotlib’ (2018).Wickham, H. et al. Welcome to the Tidyverse. J. Open Source Softw. 4, 1686 (2019).Article 
ADS 

Google Scholar  More

Study populationThe study is based on data from mothers and children participating in a longitudinal Southwest Finland cohort, Steps to Healthy development of Children (the STEPS Study) (described in detail in Lagström et al.31). The STEPS study is an ongoing population-based and multidisciplinary study that investigates children’s physical, psychological and social development, starting from pregnancy and continuing until adolescence. All Finnish- and Swedish-speaking mothers delivering a child between 1 January, 2008 and 31 March, 2010 in the Hospital District of Southwest Finland formed the cohort population (in total, 9811 mothers and their 9936 children). Altogether, 1797 mothers with 1805 neonates volunteered as participants for the intensive follow-up group of the STEPS Study. Mothers were recruited by midwives either during the first trimester of pregnancy while visiting maternity health care clinics, or after delivery on the maternity wards of Turku University Hospital or Salo Regional Hospital, or by a letter mailed to the mothers. The participating mothers differ slightly from the whole cohort population in some background characteristics (being older, with first-born child and higher socioeconomic status)31. The ethics committee of the Hospital District of Southwest Finland has approved the STEPS Study (2/2007) and all methods were performed in accordance with relevant guidelines and regulations. Written informed consent was obtained from all the participants and, for children, from one parent for study participation. Subjects have been and are free to withdraw from the study at any time without any specific reason. The STEPS Study have the appropriate government authorization to the use of the National birth register (THL/974/5.05.00/2017).Breastmilk collection and HMO analysisMothers from the STEPS Study were asked to collect breastmilk samples when the infant was approximately 3 months old. In total, 812 of the 1797 mothers (45%) provided a breastmilk sample. There were only slight differences in maternal and child characteristics between the participants providing breastmilk samples and the total STEPS Study cohort40. Altogether, 795 breastmilk samples were included in this study (excluding the duplicate observations and the 2nd born twins, samples with technical unclarity or insufficient sample quantity, one breastmilk sample collected notably later than the other samples, at infant age of 14.5 months (range for the other breastmilk samples: 0.6–6.07 months), one sample with missing information on the date of collection and six mothers missing data on residential green environment) (Supplementary Fig. 2). Mothers received written instructions for the collection of breastmilk samples: samples were collected by manual expression in the morning from one single breast, first milking a few drops to waste before collecting the actual sample (~ 10 ml) into a plastic container (pre-feed sample). The samples were stored in the fridge and the mothers brought the samples to the research center or the samples were collected from their homes on the day of sampling. All samples were frozen and stored at − 70 °C until analysis.High Performance Liquid Chromatography (HPLC) was used to identify HMOs in breastmilk as previously described40,57,58 at the University of California, San Diego (methods described in detail in Berger et al.58). Milk samples were spiked with raffinose (a non-HMO carbohydrate) as an internal standard to allow absolute quantification. HMOs were extracted by high-throughput solid-phase extraction, fluorescently labelled, and measured using HPLC with fluorescent detection (HPLC-FLD)58. Absolute concentrations for the 19 HMOs were calculated based on standard retention times and corrected for internal standard recovery. Quantified HMOs included: 2′-fucosyllactose (2′FL), 3-fucosyllactose (3FL), lacto-N-neotetraose (LNnT), 3′-sialyllactose (3′SL), difucosyllactose (DFlac), 6′-sialyllactose (6′SL), lacto-N-tetraose (LNT), lacto-Nfucopentaose (LNFP) I, LNFP II, LNFP III, sialyl-LNT (LST) b, LSTc, difucosyllacto-LNT (DFLNT), lacto-N-hexaose (LNH), disialyllacto-N-tetraose (DSLNT), fucosyllacto-Nhexaose (FLNH), difucosyllacto-N-hexaose (DFLNH), fucodisialyllacto-lacto-N-hexaose (FDSLNH) and disialyllacto-N-hexaose (DSLNH). HMOs were also summed up to seven groups based on structural features: small HMOs (2′FL, 3FL, 3′SL, 6′SL, and DFLac), type 1 HMOs (LNT, LNFP I, LNFP II, LSTb, DSLNT), type 2 HMOs (LNnT, LNFP III, LSTc), α-1-2-fucosylated HMOs (2’FL, LNFP I), terminal α-2-6-sialylated HMOs (6′SL, LSTc), internal α-2-6-sialylated HMOs (DSLNT, LSTb), terminal α-2-3-sialylated HMOs (3′SL, DSLNT). The total concentration of HMOs was calculated as the sum of the 19 oligosaccharides. HMO-bound fucose and HMO-bound sialic acid were calculated on a molar basis. The proportion of each HMO comprising the total HMO concentration was also calculated. HMO Simpson’s diversity was calculated as Simpson’s Reciprocal Index 1/D, which is the reciprocal sum of the square of the relative abundance of each of the measured 19 HMOs57,59. The higher the diversity value, the more heterogenous is the HMO composition in the sample.Properties of the residential green environmentThe selected residential green environment variables measure the properties of the green environments surrounding the homes of the participants and do not include any measures of the house characteristics, indoor environment or the actual use of green spaces by the participants. The residential green environment variables were selected due to their previously observed associations with residential microbiota and health33,34,35. The variables of the residential green environments were derived from multispectral satellite images series, with a 30 m × 30 m of spatial resolution (Landsat TM 5, National Aeronautics and Space Administration—NASA) and land cover data (CORINE). We used Landsat TM images obtained over the summertime (June–August, greenest months in Finland), to minimize the seasonal variation of living vegetation and cloud cover as well as to better identify vegetation areas and maximise the contrast in our estimated exposure. In each selected Landsat TM 5 images, the cloud was masked out, and the Normalized Difference Vegetation Index (NDVI)36 was calculated. The final NDVI map was the mean of NDVI images collected over three consecutive years (2008–2010), to make an NDVI map with non-missing values due to cloud cover for the study area. NDVI map measures the vegetation cover, vitality and density. The NDVI can get values ranging from − 1 to 1 where values below zero represent water surfaces, values close to zero indicate areas with low intensity of living vegetation and values close to one indicate high abundance of living vegetation. For the analyses, areas covered by water were removed and the value ranged from 0 to 1, to prevent negative values for underestimating the greenness values of the residential area like in some prior studies60. We assumed that summertime NDVI identified the green space and vegetation density well, but greenness intensity might vary seasonally.Second, we used calculated indicators related to the diversity and naturalness of the land cover from CORINE Land Cover data sets of the year 201261. The 12 land cover types include: (1) Residential area, (2) Industrial/commercial area, (3) Transport network, (4) Sport/leisure, (5) Agriculture, (6) Broad-leave forest, (7) Coniferous forest, (8) Mixed forest, (9) Shrub/grassland, (10) Bare surface, (11) Wetland, and (12) Water bodies. From this information, we calculated two vegetation cover indexes. The Vegetation Cover Diversity Index (Simpson’s Diversity Index of Vegetation Cover, VCDI)37, only includes vegetation classes from CORINE land cover types (categories 5–9 and 11). VCDI approaches 1 as the number of different vegetation classes increases and the proportional distribution of area among the land use classes becomes more equitable. Furthermore, because we were particularly interested in the natural vegetation cover in the residential area, we calculated the area-weighted Naturalness Index (NI)38. This is an integrated indicator used to measure the human impact and degree of all human interventions on ecological components. The index is based on CORINE Land Cover data but reclassified to 15 classes. Residential areas have been divided to two classes: Continuous residential area (High density buildings) and Discontinuous residential area (Low density, mostly individual houses area). Agricultural area has also been divided to two classes: Agricultural area (Cropland) and Pasture as well as class 9 (Shrub/grassland) has been separated to Woodland and Natural grassland. Assignment of CORINE Land Cover classes to degrees of naturalness has been made based on Walz and Stein 201438. The area-weighted NI range from 1 to 7, where low values represent low human impact (≤ 3 = Natural), medium values moderate human impact (4–5 = Semi-natural) and high values strong human impact (6–7 = Non-Natural). To ease the interpretation of results and to correspond to the same direction than the other residential green environment variables, we have reverse-scaled the NI values, so that higher values illustrate more natural residential areas.Background factorsAs genetics is strongly linked to HMO composition, maternal secretor status was determined by high abundance (secretor) or near absence (non-secretor) of the HMO 2’FL in the breastmilk samples. Mothers with active secretor (Se) genes and FUT2 enzyme produce high amounts of α-1-2-fucosylated HMOs such as 2′-fucosyllactose (2′FL), whereas in the breastmilk of non-secretor mothers these HMOs are almost absent. Beyond genetics, other maternal and infant characteristics may influence HMO composition. So far, several associations have been reported, including lactation stage, maternal pre-pregnancy BMI, maternal age, parity, maternal diet, mode of delivery, infant gestational age and infant sex22,40. Information on the potential confounding factors, child sex, birth weight, maternal age at birth, number of previous births, marital status, maternal occupational status, smoking during pregnancy (before and during pregnancy), maternal pre-pregnancy BMI, mode of delivery, duration of pregnancy and maternal diseases [including both maternal disorders predominantly related to pregnancy such as pre-eclampsia and gestational diabetes and chronic diseases (diseases of the nervous, circulatory, respiratory, digestive, musculoskeletal and genitourinary systems, cancer and mental and behavioral disorders, according to ICD-10 codes, i.e. WHO International Classification of Diseases Tenth Revision)], were obtained from Medical Birth Registers. Self-administered questionnaires upon recruitment provided information on family net income and maternal education level. Those who had no professional training or a maximum of an intermediate level of vocational training (secondary education) were classified as “basic”. Those who had studied at a University of Applied Sciences or higher (tertiary education) were classified as “advanced”. The advanced level included any academic degree (bachelor’s, master’s, licentiate or doctoral degree). Maternal diet quality was assessed in late pregnancy using the Index of Diet Quality (IDQ62) which measures adherence to health promoting diet and nutrition recommendations. The IDQ score was used in its original form by setting the statistically defined cut-off value at 10, with scores below 10 points indicating unhealthy diets and non-adherence and scores of 10–15 points indicating a health-promoting diet and adherence dietary guidelines. Lactation time postpartum (child age) and season were received from the recorded breastmilk collection dates. Lactation status (exclusive/partial/unknown breastfeeding) at the time of breastmilk collection were gathered from follow-up diaries. From partially breastfeeding mothers (n = 277) 253 had started formula feeding and 28 solids at the time of milk collection. Last, a summary z score representing socio-economic disadvantage in the residential area was obtained from Statistics Finland grid database for the year 2009 and is based on the proportion of adults with low level of education, the unemployment rate, and proportion of people living in rented housing at each participant’s residential area55.Statistical analysesTo harmonize the residential green environment variables we calculated the mean values for NDVI, VCDI and NI in 750 × 750 m squares (and 250 × 250 m) around participant homes in a Geographical Information System (QGIS, www.qgis.org). The same grid sizes were used to calculate residential socioeconomic disadvantage in the residential area55 at the time of child birth. The geographical coordinates (latitude/longitude) of the cohort participants’ home address (795 mothers) were obtained from the Population Register Centre at the time of their child birth.The background characteristics of the mothers and children are given as means and standard deviations (SD) for continuous variables and percentages for categorical variables. Due to non-normal distribution, natural logarithmic transformation was performed for all HMO variables (19 individual components, sum of HMOs, HMO-bound sialic acid, HMO-bound fucose and HMO groups (all in nmol/mL)) except for HMO diversity. Associations between each background factor and HMO diversity and 19 individual HMO components were analysed with univariate generalized linear models to identify factors independently associated with HMO composition. All factors demonstrating a significant association (p  More

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.This is a summary of: Yin, J. et al. Future socio-ecosystem productivity threatened by compound drought– heatwave events. Nat. Sustain. https://doi.org/10.1038/10.1038/s41893-022-01024-1 (2023). More

Coleman, D. J. et al. Limnol. Oceanogr. Lett. 7, 140–149 (2022).Article 

Google Scholar 
Noyce, G. L. et al. https://doi.org/10.1038/s41561-022-01070-6 (2022).Kirwan, M. L. & Megonigal, J. P. Nature 504, 53–60 (2013).Article 

Google Scholar 
Morris, J. T., Sundareshwar, P. V., Nietch, C. T., Kjerve, B. & Cahoon, D. R. Ecology 83, 2869–2877 (2002).Article 

Google Scholar 
Noyce, G. L., Kirwan, M. L., Rich, R. L. & Megonigal, J. P. Proc. Natl Acad. Sci. 116, 21623–21628 (2019).Article 

Google Scholar 
Langley, J. A., McKee, K. L., Cahoon, D. R., Cherry, J. A. & Megonigal, J. P. Proc. Natl Acad. Sci. 106, 182–6186 (2009).Article 

Google Scholar 
Dean, J. F. et al. Rev. Geophys. 56, 207–250 (2018).Article 

Google Scholar 
IPCC Climate Change 2021: The Physical Science Basis (eds Masson-Delmotte, V. et al.) (Cambridge University Press, 2021).Lin, Y. et al. Water Res. 205, 117682 (2021).Article 

Google Scholar 
Zakharova, L., Meyer, K. M. & Seifan, M. Ecol. Modell. 407, 108703 (2019).Article 

Google Scholar  More

Sample collection and postmortem evaluationBald eagle carcasses, and/or oropharyngeal and cloacal swabs were collected in the field and submitted to the Southeastern Cooperative Wildlife Disease Study Research and Diagnostic Service. In some cases, live bald eagles were found moribund and transported to wildlife rehabilitation clinics and either died in transit or soon after arrival. Carcasses underwent postmortem evaluation, including gross and histopathology. Tissue samples [heart, brain, kidney, spleen, lung, adrenal gland, pancreas, liver, small and large intestine, and cloacal bursa (if present)] were fixed in 10% neutral buffered formalin and routinely processed for histopathology23 at the Athens Veterinary Diagnostic Laboratory. Histopathology was assessed by a board-certified veterinary pathologist.Additional bald eagle and waterfowl species mortality dataData on wild bird deaths attributed to highly pathogenic influenza A viruses were retrieved from the U.S. Department of Agriculture, Animal and Plant Health Inspection Service website, at: https://www.aphis.usda.gov/aphis/ourfocus/animalhealth/animal-disease-information/avian/avian-influenza/hpai-2022/2022-hpai-wild-birds. These data are publicly available and include state, county, date detected, and species of individual birds that tested positive for HP IAV.ImmunohistochemistryImmunohistochemistry (IHC) for avian influenza virus was performed in select cases on brain, pancreas, spleen, liver, and/or adrenal gland at the Athens Veterinary Diagnostic Laboratory. IHC was performed on an automated stainer (Nemesis 3600, Biocare Medical). Polyclonal antiserum against influenza A virus was used as the primary antibody (ab155877, Abcam), diluted 1:3000, and incubated for 60 min at 37 °C with agent-positive control. Antigen retrieval was with Target Retrieval Solution (S2367, Dako) pH (10x) at 110 °C for 15 min. Enzyme blockage was via 3% H2O2 for 20 min (H324-500, Fisher Scientific); protein blockage was with Universal Blocking Reagent (10x) Power Block diluted at 1:10 for 5 min (HK085-5 K, BioGenex); link was by biotinylated rabbit anti-goat (BA-5000, Vector) at a 1:100 dilution for 10 min with 4 + streptavidin alkaline phosphatase label for 10 min (AP605H, BioCare Medical). Staining was with warp red chromogen kit for 5 min (WR8065, BioCare Medical). Known influenza A-virus positive control tissues were tested alongside each case.Polymerase chain reactionOropharyngeal and cloacal swabs from bald eagle carcasses were pooled for each individual eagle and tested by real-time reverse transcription polymerase chain reaction (rRT-PCR). Briefly, swabs samples were extracted with the KingFisher magnetic particle processer using the MagMAX-96 AI/ND Viral RNA isolation Kit (Ambion/Applied Biosystems, Foster City, CA) following a modified MagMAX-S protocol24. Resultant nucleic acids were screened against primers specific for H5 IAV in rRT-PCR; samples that yielded a cycle threshold value  More

FAO. Handbook on crop statistics: improving methods for measuring crop area, production and yield. (FAO, Rome, Italy, 2018).FAO. Land use statistics and indicators: global, regional and county trends 1990-2019. FAOSTAT Anal. Brief Ser. No 28 (2021).Potapov, P. et al. Global maps of cropland extent and change show accelerated cropland expansion in the twenty-first century. Nat. Food 1–10 (2021) https://doi.org/10.1038/s43016-021-00429-zFAO. A system of integrated agricultural censuses and surveys. (FAO, 2005).FAO. Land use statistics and indicators. Global, regional and country trends, 2000–2020. (FAO, Rome, Italy, 2022).Loveland, T. R. et al. Development of a global land cover characteristics database and IGBP DISCover from 1 km AVHRR data. Int. J. Remote Sens. 21, 1303–1330 (2000).Article 

Google Scholar 
Zanaga, D. et al. ESA WorldCover 10 m 2020 v100. (2021) https://doi.org/10.5281/zenodo.5571936Cochran, W. G. Sampling techniques. (Wiley, 1977).Stehman, S. V. Estimating area and map accuracy for stratified random sampling when the strata are different from the map classes. Int. J. Remote Sens. 35, 4923–4939 (2014).Article 

Google Scholar 
Tsujino, R., Kaijisa, T. & Yumoto, T. Causes and history of forest loss in Cambodia. Int. For. Rev. 21, 372–384 (2019).
Google Scholar 
Hu, Q. et al. Global cropland intensification surpassed expansion between 2000 and 2010: A spatio-temporal analysis based on GlobeLand30. Sci. Total Environ. 746, 141035 (2020).Grainger, A. Difficulties in tracking the long-term global trend in tropical forest area. Proc. Natl Acad. Sci. 105, 818–823 (2008).Article 
ADS 
CAS 

Google Scholar 
FAO. FAOSTAT. https://www.fao.org/faostat/en/#home (2021). More