Philippa Kaur

Sloan, R. E. in Essays on palaeontology in honour of Loris Shano Russell (ed C. S. Churcher) 134–155 (Royal Ontario Museum, 1976).Dodson, P. J. A faunal review of the Judith River (Oldman) Formation, Dinosaur Provincial Park, Alberta. Mosasaur 1, 89–118 (1983).
Google Scholar 
Clemens, W. A. in Dynamics of extinction (ed D. K. Elliott) 63–85 (John Wiley & Sons, 1986).Dodson, P. J. & Tatarinov, L. P. in The Dinosauria (eds D. B. Weishampel, P. J. Dodson, & H. Osmólska) 55–62 (University of California Press, 1990).Lehman, T. M. in Dinofest International (eds D. L. Wolberg, E. Stump, & G. D. Rosenberg) 223–240 (Philadelphia Academy of Natural Sciences, 1997).Lehman, T. M. in Mesozoic Vertebrate Life (eds D. H. Tanke & K. Carpenter) 310–328 (Indiana University Press, 2001).Sampson, S. D. et al. New horned dinosaurs from Utah provide evidence for intracontinental dinosaur endemism. PLoS ONE 5, e12292. https://doi.org/10.1371/journal.pone.0012292 (2010).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Mannion, P. D., Upchurch, P., Carrano, M. T. & Barrett, P. M. Testing the effect of the rock record on diversity: a multidisciplinary approach to elucidating the generic richness of sauropodomorph dinosaurs through time. Biol. Rev. 86, 157–181. https://doi.org/10.1111/j.1469-185X.2010.00139.x (2011).Article 
PubMed 

Google Scholar 
Upchurch, P., Mannion, P. D., Benson, R. B. J., Butler, R. J. & Carrano, M. T. Geological and anthropogenic controls on the sampling of the terrestrial fossil record: a case study from the Dinosauria. Geol. Soc. Spec. Publ 358, 209–240. https://doi.org/10.1144/SP358.14 (2011).Article 

Google Scholar 
Haq, B. U. Cretaceous eustasy revisited. Glob. Planet. Change 113, 44–58. https://doi.org/10.1016/j.gloplacha.2013.12.007 (2014).ADS 
Article 

Google Scholar 
Miller, K. G., Barrera, E., Olsson, R. K., Sugarman, P. J. & Savin, S. M. Does ice drive early Maastrichtian eustasy?. Geology 27, 783. https://doi.org/10.1130/0091-7613(1999)027%3c0783:dideme%3e2.3.co;2 (1999).ADS 
Article 

Google Scholar 
Catuneanu, O., Sweet, A. R. & Miall, A. D. Reciprocal stratigraphy of the Campanian-Paleocene Western Interior of North America. Sediment. Geol. 134, 235–255. https://doi.org/10.1016/S0037-0738(00)00045-2 (2000).ADS 
Article 

Google Scholar 
Smith, R. L. Ash flows. Geol. Soc. Am. Bull. 71, 795–841. https://doi.org/10.1130/0016-7606(1960)71[795:af]2.0.co;2 (1960).ADS 
Article 

Google Scholar 
Smedes, H. W. Geology and igneous petrology of the northern Elkhorn mountains. 116 (United States Geological Survey Professional Paper 510 1966).Rutland, C., Smedes, H. W., Tilling, R. I. & Greenwood, W. R. in Cordilleran volcanism, plutonism, and magma generation at various crustal levels, Montana and Idaho. 28th International Geological Congress, Field Trip Guidebook T337 (ed D. W. Hyndman) 16–31 (American Geophysical Union, 1989).Harlan, S. S. et al. 40Ar/39Ar and K-Ar Geochronology and Tectonic Significance of the Upper Cretaceous Adel Mountain Volcanics and Spatially Associated Tertiary Igneous Rocks, Northwestern Montana. 29 (United States Geological Survey Professional Paper 1696, 2005).Breyer, J. A. et al. Evidence for late cretaceous volcanism in Trans-Pecos Texas. J. Geol. 115, 243–251. https://doi.org/10.1086/510640 (2007).ADS 
Article 

Google Scholar 
Jennings, G. R., Lawton, T. E. & Clinkscales, C. A. Late cretaceous U-Pb tuff ages from the, Skunk Ranch Formation and their implications for age of Laramide deformation, Little Hatchet Mountains, southwestern New Mexico, USA. Cretac. Res. 43, 18–25. https://doi.org/10.1016/j.cretres.2013.02.001 (2013).Article 

Google Scholar 
Roberts, E. M. & Hendrix, M. S. Taphonomy of a petrified forest in the Two Medicine Formation (Campanian), northwest Montana: implications for palinspastic restoration of the Boulder batholith and Elkhorn Mountains Volcanics. Palaios 15, 476–482. https://doi.org/10.2307/3515516 (2000).ADS 
Article 

Google Scholar 
Sewall, J. O. et al. Climate model boundary conditions for four Cretaceous time slices. Clim. Past. 3, 647–657. https://doi.org/10.5194/cp-3-647-2007 (2007).Article 

Google Scholar 
Bertog, J. Stratigraphy of the lower Pierre Shale (Campanian): implications for the tectonic and eustatic controls on facies distributions. J. Geol. Res. 2010, 910243. https://doi.org/10.1155/2010/910243 (2010).ADS 
Article 

Google Scholar 
Fricke, H. C., Foreman, B. Z. & Sewall, J. O. Integrated climate model-oxygen isotope evidence for a North American monsoon during the Late Cretaceous. Earth Planet. Sci. Lett. 289, 11–21. https://doi.org/10.1016/j.epsl.2009.10.018 (2010).ADS 
CAS 
Article 

Google Scholar 
Obradovich, J. D. in Evolution of the Western Interior Basin (eds W. G. E. Caldwell & E. G. Kaufman) 379–396 (Geological Association of Canada Special Paper 39, 1993).Cobban, W. A., Walaszczyk, I., Obradovich, J. D. & McKinney, K. C. A USGS Zonal Table for the Upper Cretaceous Middle Cenomanian–Maastrichtian of the Western Interior of the United States Based on Ammonites, Inoceramids, and Radiometric Ages. (United States Geological Survey Open-File Report 2006–1250, 2006).Rogers, R. R., Swisher, C. C. & Horner, J. R. 40Ar/39Ar age and correlation of the nonmarine Two Medicine Formation (Upper Cretaceous), northwestern Montana, U.S.A. Can J Earth Sci 30, 1066–1075. https://doi.org/10.1139/e93-090 (1993).CAS 
Article 

Google Scholar 
Goodwin, M. B. & Deino, A. L. The first radiometric ages from the Judith River Formation (Upper Cretaceous), Hill County, Montana. Can. J. Earth Sci. 26, 1384–1391. https://doi.org/10.1139/e89-118 (1989).ADS 
CAS 
Article 

Google Scholar 
Thomas, R. G., Eberth, D. A., Deino, A. L. & Robinson, D. Composition, radioisotopic ages, and potential significance of an altered volcanic ash (bentonite) from the Upper Cretaceous Judith River Formation, Dinosaur Provincial Park, southern Alberta, Canada. Cretac. Res. 11, 125–162. https://doi.org/10.1016/s0195-6671(05)80030-8 (1990).CAS 
Article 

Google Scholar 
Roberts, E. M., Deino, A. L. & Chan, M. A. 40Ar/39Ar age of the Kaiparowits Formation, southern Utah, and correlation of contemporaneous Campanian strata and vertebrate faunas along the margin of the Western Interior Basin. Cretac. Res. 26, 307–318. https://doi.org/10.1016/j.cretres.2005.01.002 (2005).Article 

Google Scholar 
Fassett, J. E. & Steiner, M. B. in Mesozoic Geology and Paleontology of the Four Corners Region (eds O. Anderson, B. S. Kues, & S. G. Lucas) 239–247 (New Mexico Geological Society 48th Field Conference Guidebook, 1997).Sprain, C. J., Renne, P. R., Wilson, G. P. & Clemens, W. A. High-resolution chronostratigraphy of the terrestrial Cretaceous-Paleogene transition and recovery interval in the Hell Creek region, Montana. Geol. Soc. Am. Bull. 127, 393–409. https://doi.org/10.1130/B31076.1 (2015).ADS 
Article 

Google Scholar 
Clyde, W. C., Ramezani, J., Johnson, K. R., Bowring, S. A. & Jones, M. M. Direct high-precision U-Pb geochronology of the end-Cretaceous extinction and calibration of Paleocene astronomical timescales. Earth Planet. Sci. Lett. 452, 272–280. https://doi.org/10.1016/j.epsl.2016.07.041 (2016).ADS 
CAS 
Article 

Google Scholar 
Wang, T. T. et al. High-precision U-Pb geochronologic constraints on the Late Cretaceous terrestrial cyclostratigraphy and geomagnetic polarity from the Songliao Basin, Northeast China. Earth Planet. Sci. Lett. 446, 37–44. https://doi.org/10.1016/j.epsl.2016.04.007 (2016).ADS 
CAS 
Article 

Google Scholar 
Blakey, R. C. Paleogeography and Paleotectonics of the Western Interior Seaway, Jurassic-Cretaceous of North America. (American Association of Petroleum Geologists Search and Discovery Article 30392, 2014).Archibald, J. D. Dinosaur Extinction and the End of an Era: What the Fossils Say 240 (Columbia University Press, London, 1996).
Google Scholar 
Currie, P. J. & Russell, D. A. in Dinosaur Provincial Park: A Spectacular Ancient Ecosystem Revealed (eds P. J. Currie & E. B. Koppelhus) 537–569 (Indiana University Press, 2005).Eberth, D. A. & Hamblin, A. P. Tectonic, stratigraphic, and sedimentologic significance of a regional discontinuity in the upper Judith River Group (Belly River Wedge) of southern Alberta, Saskatchewan, and northern Montana. Can. J. Earth Sci. 30, 174–200. https://doi.org/10.1139/e93-016 (1993).ADS 
Article 

Google Scholar 
Eberth, D. A. in Dinosaur Provincial Park: A spectacular Ancient Ecosystem Revealed (eds P. J. Currie & E. B. Koppelhus) Ch. 3, 54–82 (Indiana University Press, 2005).Eberth, D. A. Origin and significance of mud-filled incised valleys (Upper Cretaceous) in southern Alberta, Canada. Sedimentology 43, 459–477. https://doi.org/10.1046/j.1365-3091.1996.d01-15.x (1996).ADS 
Article 

Google Scholar 
Russell, D. A. A new specimen of Stenonychosaurus from the Oldman Formation (Cretaceous) of Alberta. Can. J. Earth Sci. 6, 595–612. https://doi.org/10.1139/e69-059 (1969).ADS 
Article 

Google Scholar 
Dodson, P. Sedimentology and taphonomy of Oldman formation (Campanian), Dinosaur-Provincial-Park, Alberta (Canada). Palaeogeogr. Palaeocl. 10, 21–000. https://doi.org/10.1016/0031-0182(71)90044-7 (1971).Article 

Google Scholar 
Farlow, J. O. Consideration of trophic dynamics of a late cretaceous large dinosaur community (Oldman formation). Ecology 57, 841–857. https://doi.org/10.2307/1941052 (1976).Article 

Google Scholar 
Beland, P. & Russell, D. A. Paleoecology of Dinosaur-Provincial-Park (Cretaceous), Alberta, interpreted from distribution of articulated vertebrate remains. Can. J. Earth Sci. 15, 1012–1024. https://doi.org/10.1139/e78-109 (1978).ADS 
Article 

Google Scholar 
MacDonald, M., Currie, P. J. & Spencer, W. A. in Dinosaur Provincial Park: A Spectacular Ancient Ecosystem Revealed (eds P. J. Currie & E. B. Koppelhus) 478–485 (Indiana University Press, 2005).Eberth, D. A., Brinkman, D. B. & Barkas, V. in New Perspectives on Horned Dinosaurs: The Royal Tyrrell Museum Ceratopsian Symposium (eds M. J. Ryan, B. J. Chinnery-Allgeier, & D. A. Eberth) 495–508 (Indiana University Press, 2010).Mallon, J. C., Evans, D. C., Ryan, M. J. & Anderson, J. S. Megaherbivorous dinosaur turnover in the Dinosaur Park Formation (upper Campanian) of Alberta, Canada. Palaeogeogr. Palaeocl. 350, 124–138. https://doi.org/10.1016/j.palaeo.2012.06.024 (2012).Article 

Google Scholar 
Brown, C. M., Evans, D. C., Campione, N. E., O’Brien, L. J. & Eberth, D. A. Evidence for taphonomic size bias in the Dinosaur Park Formation (Campanian, Alberta), a model Mesozoic terrestrial alluvial-paralic system. Palaeogeogr Palaeocl 372, 108–122. https://doi.org/10.1016/j.palaeo.2012.06.027 (2013).Article 

Google Scholar 
Eberth, D. A. & Getty, M. A. in Dinosaur Provincial Park: A Spectacular Ancient Ecosystem Revealed (eds P. J. Currie & E. B. Koppelhus) 501–536 (Indiana University Press, 2005).Brown, C. M., Herridge-Berry, S., Chiba, K., Vitkus, A. & Eberth, D. A. High-resolution (centimetre-scale) GPS/GIS-based 3D mapping and spatial analysis of in situ fossils in two horned-dinosaur bonebeds in the Dinosaur Park Formation (Upper Cretaceous) at Dinosaur Provincial Park, Alberta, Canada. Can. J. Earth Sci. 58, 225–246. https://doi.org/10.1139/cjes-2019-0183 (2021).ADS 
Article 

Google Scholar 
Eberth, D. A., Braman, D. R. & Tokaryk, T. T. Stratigraphy, Sedimentology and vertebrate paleontology of the Judith River Formation (Campanian) near Muddy Lake, West-Central Saskatchewan. Bull. Can. Petrol. Geol. 38, 387–406 (1990).
Google Scholar 
Rogers, R. R. Sequence analysis of the Upper Cretaceous Two Medicine and Judith River formations, Montana; nonmarine response to the Claggett and Bearpaw marine cycles. J. Sediment. Res. 68, 615–631. https://doi.org/10.2110/jsr.68.604 (1998).ADS 
Article 

Google Scholar 
Rogers, R. R. Taphonomy of three dinosaur bone beds in the Upper Cretaceous Two Medicine Formation of Northwestern Montana: evidence for drought-related mortality. Palaios 5, 394–413. https://doi.org/10.2307/3514834 (1990).ADS 
Article 

Google Scholar 
Falcon-Lang, H. J. Growth interruptions in silicified conifer woods from the Upper Cretaceous Two Medicine Formation, Montana, USA: implications for palaeoclimate and dinosaur palaeoecology. Palaeogeogr. Palaeocl. 199, 299–314. https://doi.org/10.1016/S0031-0182(03)00539-X (2003).Article 

Google Scholar 
Horner, J. R. & Makela, R. Nest of juveniles provides evidence of family-structure among dinosaurs. Nature 282, 296–298. https://doi.org/10.1038/282296a0 (1979).ADS 
Article 

Google Scholar 
Horner, J. R., Varricchio, D. J. & Goodwin, M. B. Marine transgressions and the evolution of Cretaceous dinosaurs. Nature 358, 59–61. https://doi.org/10.1038/358059a0 (1992).ADS 
Article 

Google Scholar 
Sampson, S. D. Two new horned dinosaurs from the Upper Cretaceous Two Medicine Formation of Montana; With a phylogenetic analysis of the Centrosaurinae (Ornithischia:Ceratopsidae). J. Vertebr. Paleontol. 15, 743–760. https://doi.org/10.1080/02724634.1995.10011259 (1995).Article 

Google Scholar 
Carr, T. D., Varricchio, D. J., Sedlmayr, J. C., Roberts, E. M. & Moore, J. R. A new tyrannosaur with evidence for anagenesis and crocodile-like facial sensory system. Sci. Rep. 7, 44942. https://doi.org/10.1038/srep44942 (2017).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Wilson, J. P., Ryan, M. J. & Evans, D. C. A new, transitional centrosaurine ceratopsid from the Upper Cretaceous Two Medicine Formation of Montana and the evolution of the “Styracosaurus-line” dinosaurs. R. Soc. Open Sci. 7, 200284. https://doi.org/10.1098/rsos.200284 (2020).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Foreman, B. Z., Rogers, R. R., Deino, A. L., Wirth, K. R. & Thole, J. T. Geochemical characterization of bentonite beds in the Two Medicine Formation (Campanian, Montana), including a new 40Ar/39Ar age. Cretac. Res. 29, 373–385. https://doi.org/10.1016/j.cretres.2007.07.001 (2008).Article 

Google Scholar 
Varricchio, D. J. et al. in Large Meteorite Impacts and Planetary Evolution IV Vol. 465 (eds R. L. Gibson & W. U. Reimold) 269–299 (Geological Society of America Special Paper 465, 2010).Meek, F. B. & Hayden, F. V. Descriptions of new species of acephala and gasteropoda, from the tertiary formations of Nebraska Territory, with some general remarks on the geology of the country about the sources of the Missouri River. Ceratites Americanus. Proc. Acad. Nat. Sci. Phila. 8, 111–126 (1856).
Google Scholar 
Hayden, F. V. Notes explanatory of a map and section illustrating the geologic structure of the country bordering the Missouri River from the mouth of the Platte River to Fort Benton. Proc. Acad. Natl. Sci. Phila. 9, 109–148 (1857).
Google Scholar 
Hayden, F. V. in [Fourth Annual] Preliminary Report of the United States Geological Survey of Wyoming and portions of contiguous Territories 85–98 (U.S. Geological Survey, 1871).Dawson, G. M. in Report on the Geology and Resources of the Region in the Vicinity of the Forty-Ninth Parallel, from the Lake of the Woods to the Rocky Mountains 1–18 (British North American Boundary Commission, 1875).Stanton, T. W., Hatcher, J. B. & Knowlton, F. H. Geology and Paleontology of the Judith River Beds (United States Geological Survey Bulletin No. 257, 1905).Bowen, C. F. in Shorter Contributions to General Geology 1914 95–153 (United States Geological Survey Professional Paper 90, 1915).Waage, K. M. in The Cretaceous System in the Western Interior of North America: The Proceedings of an International Symposium Organized by the Geological Association of Canada, Saskatoon, Saskatchewan, May 23–26, 1973 (ed W. G. E. Caldwell) 55–81 (Geological Association of Canada Special paper 13, 1975).Leidy, J. Notice of remains of extinct reptiles and fishes, discovered by Dr. FV Hayden in the Bad Lands of the Judith River, Nebraska Territory. Proc. Acad. Nat. Sci. Phila. 8, 72–73. https://doi.org/10.5281/zenodo.1038128 (1856).Article 

Google Scholar 
Leidy, J. Extinct vertebrata from the Judith River and Great Lignite formations of Nebraska. Trans. Am. Philos. Soc. 11, 139–154. https://doi.org/10.2307/3231936 (1860).Article 

Google Scholar 
Cope, E. D. On some extinct reptiles and Batrachia from the Judith River and Fox Hills beds of Montana. Proc. Acad. Natl. Sci. Phila. 28, 340–359 (1876).
Google Scholar 
Sternberg, C. H. Notes on the fossil vertebrates collected on the Cope expedition to the Judith River and Cow Island beds, Montana, in 1876. Science 40, 134–135. https://doi.org/10.1126/science.40.1021.134 (1914).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Sahni, A. The vertebrate fauna of the Judith River Formation, Montana. Bull. Am. Mus. Nat. Hist. 147, 325–412 (1972).
Google Scholar 
Tschudy, B. D. Palynology of the upper Campanian (Cretaceous) Judith River Formation, north-central Montana. 42 (United States Geological Survey Professional Paper 770, 1973).Case, G. R. A new Selachian Fauna from the Judith River formation (Campanian) of Montana. Palaeontogr. Abt. A Band A 160, 176–205 (1978).
Google Scholar 
Horner, J. R. A new hadrosaur (Reptilia, Ornithischia) from the Upper Cretaceous Judith River Formation of Montana. J. Vertebr. Paleontol. 8, 314–321. https://doi.org/10.1080/02724634.1988.10011714 (1988).Article 

Google Scholar 
Fiorillo, A. R. & Currie, P. J. Theropod teeth from the Judith River formation (Upper Cretaceous) of south-central Montana. J. Vertebr. Paleontol. 14, 74–80. https://doi.org/10.1080/02724634.1994.10011539 (1994).Article 

Google Scholar 
Prieto-Marquez, A. New information on the cranium of Brachylophosaurus canadensis (Dinosauria, Hadrosauridae), with a revision of its phylogenetic position. J. Vertebr. Paleontol. 25, 144–156. https://doi.org/10.1671/0272-4634(2005)025[0144:Niotco]2.0.Co;2 (2005).Article 

Google Scholar 
Fricke, H. C., Rogers, R. R., Backlund, R., Dwyer, C. N. & Echt, S. Preservation of primary stable isotope signals in dinosaur remains, and environmental gradients of the Late Cretaceous of Montana and Alberta. Palaeogeogr. Palaeocl. 266, 13–27. https://doi.org/10.1016/j.palaeo.2008.03.030 (2008).Article 

Google Scholar 
Fricke, H. C., Rogers, R. R. & Gates, T. A. Hadrosaurid migration: inferences based on stable isotope comparisons among Late Cretaceous dinosaur localities. Paleobiology 35, 270–288. https://doi.org/10.1666/08025.1 (2009).Article 

Google Scholar 
Tweet, J. S., Chin, K., Braman, D. R. & Murphy, N. L. Probable gut contents within a specimen of Brachylophosaurus canadensis (Dinosauria: Hadrosauridae) from the Upper Cretaceous Judith River formation of Montana. Palaios 23, 624–635. https://doi.org/10.2110/palo.2007.p07-044r (2008).ADS 
Article 

Google Scholar 
Ryan, M. J., Evans, D. C., Currie, P. J. & Loewen, M. A. A new chasmosaurine from northern Laramidia expands frill disparity in ceratopsid dinosaurs. Naturwissenschaften 101, 505–512. https://doi.org/10.1007/s00114-014-1183-1 (2014).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Arbour, V. M. & Evans, D. C. A new ankylosaurine dinosaur from the Judith River formation of Montana, USA, based on an exceptional skeleton with soft tissue preservation. R. Soc. Open Sci. 4, 161086. https://doi.org/10.1098/rsos.161086 (2017).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Chiba, K., Ryan, M. J., Fanti, F., Loewen, M. A. & Evans, D. C. New material and systematic re-evaluation of Medusaceratops lokii (Dinosauria, Ceratopsidae) from the Judith River formation (Campanian, Montana). J. Paleontol. 92, 272–288. https://doi.org/10.1017/jpa.2017.62 (2017).Article 

Google Scholar 
Rogers, R. R. et al. Age, correlation, and lithostratigraphic revision of the Upper Cretaceous (Campanian) Judith River formation in its type area (north-central Montana), with a comparison of low- and high-accommodation alluvial records. J. Geol. 124, 99–135. https://doi.org/10.1086/684289 (2016).ADS 
CAS 
Article 

Google Scholar 
Lawton, T. F., Pollock, S. L. & Robinson, R. A. J. Integrating sandstone petrology and nonmarine sequence stratigraphy: application to the late cretaceous fluvial systems of southwestern Utah, USA. J. Sediment. Res. 73, 389–406. https://doi.org/10.1306/100702730389 (2003).ADS 
Article 

Google Scholar 
Jinnah, Z. A. et al. New 40Ar/39Ar and detrital zircon U-Pb ages for the Upper Cretaceous Wahweap and Kaiparowits formations on the Kaiparowits Plateau, Utah: implications for regional correlation, provenance, and biostratigraphy. Cretac. Res. 30, 287–299. https://doi.org/10.1016/j.cretres.2008.07.012 (2009).Article 

Google Scholar 
Beveridge, T. L. et al. Refined geochronology and revised stratigraphic nomenclature of the Upper Cretaceous Wahweap Formation, Utah, U.S.A. and the age of early Campanian vertebrates from southern Laramidia. Palaeogeogr. Palaeoclimatol. Palaeoecol. 591, 110876. https://doi.org/10.1016/j.palaeo.2022.110876 (2022).Article 

Google Scholar 
Jinnah, Z. A. & Roberts, E. M. Facies associations, paleoenvironment, and base-level changes in the Upper Cretaceous Wahweap Formation, Utah, USA. J. Sediment. Res. 81, 266–283. https://doi.org/10.2110/jsr.2011.22 (2011).ADS 
Article 

Google Scholar 
Gregory, H. E. & Moore, R. C. The Kaiparowits region, a geographic and geologic reconnaissance of parts of Utah and Arizona. Report No. 164, 161 (United States Geological Survey Professional Paper 164, 1931).Lohrengel, C. F. II. Palynology of Kaiparowits Formation, Garfield County, Utah. AAPG Bull. 53, 729–729. https://doi.org/10.1306/5d25c75f-16c1-11d7-8645000102c1865d (1969).Article 

Google Scholar 
Roberts, E. M. Facies architecture and depositional environments of the Upper Cretaceous Kaiparowits Formation, southern Utah. Sediment. Geol. 197, 207–233. https://doi.org/10.1016/j.sedgeo.2006.10.001 (2007).ADS 
Article 

Google Scholar 
Lawton, T. F. & Bradford, B. A. Correlation and provenance of Upper Cretaceous (Campanian) fluvial strata, Utah, USA, from Zircon U-Pb geochronology and petrography. J. Sediment. Res. 81, 495–512. https://doi.org/10.2110/jsr.2011.45 (2011).ADS 
Article 

Google Scholar 
Beveridge, T. L., Roberts, E. M. & Titus, A. L. Volcaniclastic member of the richly fossiliferous Kaiparowits Formation reveals new insights for regional correlation and tectonics in southern Utah during the latest Campanian. Cretac. Res. https://doi.org/10.1016/j.cretres.2020.104527 (2020).Article 

Google Scholar 
Titus, A. L. et al. in Interior Western United States (ed C. M. Dehler) 1–28 (Geological Society of America Field Guide 6, 2005).Titus, A. L. & Loewen, M. A. At the Top of the Grand Staircase: The Late Cretaceous of Southern Utah (Indiana University Press, 2013).Cifelli, R. L. Cretaceous mammals of southern Utah. I. Marsupials from the Kaiparowits Formation (Judithian). J. Vertebr. Paleontol. 10, 295–319. https://doi.org/10.1080/02724634.1990.10011816 (1990).Article 

Google Scholar 
Eaton, J., Cifelli, R., Hutchison, J. H., Kirkland, J. & Parrish, J. in Vertebrate Paleontology in Utah (ed D. D. Gillette) 345–353 (Utah Geological Survey Miscellaneous Publication 99–1, 1999).Zanno, L. E. & Sampson, S. D. A new oviraptorosaur (Theropoda, Maniraptora) from the Late Cretaceous (Campanian) of Utah. J. Vertebr. Paleontol. 25, 897–904. https://doi.org/10.1671/0272-4634(2005)025[0897:Anotmf]2.0.Co;2 (2005).Article 

Google Scholar 
Gates, T. A. & Sampson, S. D. A new species of Gryposaurus (Dinosauria : Hadrosauridae) from the late Campanian Kaiparowits Formation, southern Utah, USA. Zool J Linn Soc-Lond 151, 351–376. https://doi.org/10.1111/j.1096-3642.2007.00349.x (2007).Article 

Google Scholar 
Sampson, S. D., Lund, E. K., Loewen, M. A., Farke, A. A. & Clayton, K. E. A remarkable short-snouted horned dinosaur from the Late Cretaceous (late Campanian) of southern Laramidia. Proc. Biol. Sci. 280, 20131186. https://doi.org/10.1098/rspb.2013.1186 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Carr, T. D., Williamson, T. E., Britt, B. B. & Stadtman, K. Evidence for high taxonomic and morphologic tyrannosauroid diversity in the Late Cretaceous (Late Campanian) of the American Southwest and a new short-skulled tyrannosaurid from the Kaiparowits formation of Utah. Naturwissenschaften 98, 241–246. https://doi.org/10.1007/s00114-011-0762-7 (2011).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Zanno, L. E., Varricchio, D. J., O’Connor, P. M., Titus, A. L. & Knell, M. J. A new troodontid theropod, Talos sampsoni gen. et sp. Nov., from the Upper Cretaceous Western Interior Basin of North America. PLoS ONE 6, e24487. https://doi.org/10.1371/journal.pone.0024487 (2011).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Loewen, M. A., Irmis, R. B., Sertich, J. J., Currie, P. J. & Sampson, S. D. Tyrant dinosaur evolution tracks the rise and fall of Late Cretaceous oceans. PLoS ONE 8, e79420. https://doi.org/10.1371/journal.pone.0079420 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Wiersma, J. P. & Irmis, R. B. A new southern Laramidian ankylosaurid, Akainacephalus johnsoni gen. et sp. Nov., from the upper Campanian Kaiparowits Formation of southern Utah, USA. Peerj 6, e5016. https://doi.org/10.7717/peerj.5016 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Titus, A. L. et al. Geology and taphonomy of a unique tyrannosaurid bonebed from the upper Campanian Kaiparowits Formation of southern Utah: implications for tyrannosaurid gregariousness. PeerJ 9, e11013. https://doi.org/10.7717/peerj.11013 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Roberts, E., Sampson, S., Deino, A., Bowring, S. & Buchwaldt, R. in At the Top of the Grand Staircase: The Late Cretaceous of Southern Utah (eds A. L. Titus & M. A. Loewen) 85–106 (Indiana University Press, 2013).Fassett, J. E. & Hinds, J. S. Geology and fuel resources of the Fruitland Formation and Kirtland Shale of the San Juan Basin, New Mexico and Colorado. Report No. 676, 76 (United States Geological Survey Professional Paper 676, 1971).Fassett, J. E. in Geologic Assessment of Coal in the Colorado Plateau: Arizona, Colorado, New Mexico, and Utah (eds M. A. Kirschbaum, L. N. R. Roberts, & L. Biewick) Q1-Q132 (U.S. Geological Survey Professional Paper 1625–B, 2000).Flynn, A. G. et al. Early Paleocene magnetostratigraphy and revised biostratigraphy of the Ojo Alamo Sandstone and Lower Nacimiento Formation, San Juan Basin, New Mexico, USA. GSA Bull. 132, 2154–2174. https://doi.org/10.1130/b35481.1 (2020).Article 

Google Scholar 
Hay, O. P. On the habits and the pose of the Sauropodous dinosaurs, especially of Diplodocus. Am. Nat. 42, 672–681. https://doi.org/10.1086/278992 (1908).Article 

Google Scholar 
Gilmore, C. W. in Shorter Contributions to General Geology 1916 279–308 (United States Geological Survey Professional Paper 98-Q, 1916).Gilmore, C. W. On the Replilia of the Kirtland formation of New Mexico, with descriptions of new species of fossil turtles. Proc. U.S. Natl. Mus. 83, 159–188 (1935).Article 

Google Scholar 
Hunt, A. P. Integrated vertebrate, invertebrate and plant taphonomy of the Fossil Forest area (Fruitland and Kirtland formations: Late Cretaceous), San-Juan-County, New-Mexico, USA. Palaeogeogr. Palaeocl. 88, 85–107. https://doi.org/10.1016/0031-0182(91)90016-K (1991).Article 

Google Scholar 
Hunt, A. P. & Lucas, S. G. in New Mexico Geological Society 43rd Field Conference Guidebook Vol. 43 (eds S. G. Lucas, B. S. Kues, T. E. Williamson, & A. P. Hunt) 217–239 (New Mexico Geological Society, 1992).Fassett, J. E. & Heizler, M. T. in The Geology of the Ouray-Silverton Area (eds K. E. Karlstrom et al.) 115–121 (68th New Mexico Geological Society Field Conference Guidebook, 2017).Folinsbee, R., Lipson, J. & Baadsgaard, H. Potassium-argon dates of upper cretaceous ash falls, Alberta, Canada. Ann. N. Y. Acad. Sci. 91, 352. https://doi.org/10.1111/j.1749-6632.1961.tb35475.x (1961).ADS 
Article 

Google Scholar 
Lerbekmo, J. F. Petrology of the belly river formation, southern Alberta foothills. Sedimentology 2, 54–86. https://doi.org/10.1111/j.1365-3091.1963.tb01200.x (1963).ADS 
Article 

Google Scholar 
Min, K. W., Renne, P. R. & Huff, W. D. 40Ar/39Ar dating of Ordovician K-bentonites in Laurentia and Baltoscandia. Earth Planet. Sci. Lett. 185, 121–134. https://doi.org/10.1016/S0012-821x(00)00365-4 (2001).ADS 
CAS 
Article 

Google Scholar 
Steiger, R. H. & Jäger, E. Subcommission on geochronology: convention on the use of decay constants in geo- and cosmochronology. Earth Planet. Sci. Lett. 36, 359–362. https://doi.org/10.1016/0012-821x(77)90060-7 (1977).ADS 
CAS 
Article 

Google Scholar 
Samson, S. D. & Alexander, E. C. Calibration of the interlaboratory 40Ar-39Ar dating standard, Mmhb-1. Chem. Geol. 66, 27–34. https://doi.org/10.1016/0168-9622(87)90025-X (1987).CAS 
Article 

Google Scholar 
Deino, A. & Potts, R. Single-crystal 40Ar/39Ar dating of the Olorgesailie formation, Southern Kenya Rift. J. Geophys. Res. 95, 8453. https://doi.org/10.1029/JB095iB06p08453 (1990).ADS 
CAS 
Article 

Google Scholar 
Renne, P. R. et al. Intercalibration of standards, absolute ages and uncertainties in 40Ar/39Ar dating. Chem Geol 145, 117–152. https://doi.org/10.1016/s0009-2541(97)00159-9 (1998).ADS 
CAS 
Article 

Google Scholar 
Kuiper, K. F. et al. Synchronizing rock clocks of Earth history. Science 320, 500–504. https://doi.org/10.1126/science.1154339 (2008).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Fowler, D. W. Revised geochronology, correlation, and dinosaur stratigraphic ranges of the Santonian-Maastrichtian (Late Cretaceous) formations of the Western Interior of North America. PLoS ONE 12, e0188426. https://doi.org/10.1371/journal.pone.0188426 (2017).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Turrin, B. D. et al. in American Geophysical Union, Fall Meeting Vol. 2016 V23A–2969 (San Francisco, California, 2016).Phillips, D., Matchan, E. L., Dalton, H. & Kuiper, K. F. Revised astronomically calibrated 40Ar/39Ar ages for the Fish Canyon Tuff sanidine—closing the interlaboratory gap. Chem. Geol. 597, 120815. https://doi.org/10.1016/j.chemgeo.2022.120815 (2022).ADS 
CAS 
Article 

Google Scholar 
Eberth, D. A. & Kamo, S. L. High-precision U-Pb CA-ID-TIMS dating and chronostratigraphy of the dinosaur-rich Horseshoe Canyon Formation (Upper Cretaceous, Campanian-Maastrichtian), Red Deer River valley, Alberta, Canada. Can. J. Earth Sci. 57, 1220–1237. https://doi.org/10.1139/cjes-2019-0019 (2020).ADS 
CAS 
Article 

Google Scholar 
Gale, A. S. et al. in Geologic Time Scale 2020 (eds F. M. Gradstein, J. G. Ogg, M. D. Schmitz, & G. M. Ogg) 1023–1086 (Elsevier, 2020).Condon, D. J., Schoene, B., McLean, N. M., Bowring, S. A. & Parrish, R. R. Metrology and traceability of U-Pb isotope dilution geochronology (EARTHTIME Tracer Calibration Part I). Geochim. Cosmochim. Acta 164, 464–480. https://doi.org/10.1016/j.gca.2015.05.026 (2015).ADS 
CAS 
Article 

Google Scholar 
Mattinson, J. M. Zircon U-Pb chemical abrasion (“CA-TIMS”) method: combined annealing and multi-step partial dissolution analysis for improved precision and accuracy of zircon ages. Chem. Geol. 220, 47–66. https://doi.org/10.1016/j.chemgeo.2005.03.011 (2005).ADS 
CAS 
Article 

Google Scholar 
McLean, N. M., Condon, D. J., Schoene, B. & Bowring, S. A. Evaluating uncertainties in the calibration of isotopic reference materials and multi-element isotopic tracers (EARTHTIME Tracer Calibration Part II). Geochim. Cosmochim. Acta 164, 481–501. https://doi.org/10.1016/j.gca.2015.02.040 (2015).ADS 
CAS 
Article 

Google Scholar 
Lu, J. et al. Volcanically driven lacustrine ecosystem changes during the Carnian Pluvial Episode (Late Triassic). Proc. Natl. Acad. Sci. U.S.A. 118, e2109895118. https://doi.org/10.1073/pnas.2109895118 (2021).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Jiang, B., Harlow, G. E., Wohletz, K., Zhou, Z. & Meng, J. New evidence suggests pyroclastic flows are responsible for the remarkable preservation of the Jehol biota. Nat. Commun. 5, 3151. https://doi.org/10.1038/ncomms4151 (2014).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Gates, T. A. et al. Biogeography of terrestrial and freshwater vertebrates from the late Cretaceous (Campanian) Western Interior of North America. Palaeogeogr. Palaeocl. 291, 371–387. https://doi.org/10.1016/j.palaeo.2010.03.008 (2010).Article 

Google Scholar 
Eaton, J. G. in Stratigraphy, depositional environments; and sedimentary tectonics of the western margin, Cretaceous Western Interior Seaway (eds J. Dale Nations & J. G. Eaton) 47–63 (Geological Society of America Special Paper 260, 1991).Sankey, J. T. Late Campanian southern dinosaurs, Aguja Formation, Big Bend, Texas. J. Paleontol. 75, 208–215. https://doi.org/10.1666/0022-3360(2001)075%3c0208:Lcsdaf%3e2.0.Co;2 (2001).Article 

Google Scholar 
Sullivan, R. & Lucas, S. G. Vertebrate faunal succession in the Upper Cretaceous, San Juan Basin, New Mexico, with implications for correlations within the north American western interior. J. Vertebr. Paleontol. 23, 102a–102a (2003).
Google Scholar 
Currie, P. J. in Dinosaur Provincial Park: A Spectacular Ancient Ecosystem Revealed (eds P. J. Currie & E. B. Koppelhus) 3–33 (Indiana University Press, 2005).Kirkland, J. I. & Deblieux, D. D. in New Perspectives on Horned Dinosaurs: The Royal Tyrrell Museum Ceratopsian Symposium (eds M. J. Ryan, B. J. Chinnery-Allgeier, & D. A. Eberth) 117–140 (Indiana University Press, 2010).Miller, I. M., Johnson, K., Kline, D. E., Nichols, D. J. & Barclay, R. in At the Top of the Grand Staircase: The Late Cretaceous of southern Utah (eds A. Titus & M. Loewen) 107–131 (Indiana University Press, 2013).Tapanila, L. & Roberts, E. in At the Top of the Grand Staircase: The Late Cretaceous of Southern Utah (eds A. L. Titus & M. A. Loewen) 132–152 (Indiana University Press, 2013).Schmitt, J. & Varricchio, D. J. Volcano-tectonic partitioning of Laramidia: Influence on Campanian terrestrial environments and ecosystems. Program and Abstracts. J. Vertebr. Paleontol. 31, 188. https://doi.org/10.1080/02724634.2011.10635174 (2011).Article 

Google Scholar 
Burgener, L. et al. An extreme climate gradient-induced ecological regionalization in the Upper Cretaceous Western Interior Basin of North America. GSA Bull. https://doi.org/10.1130/b35904.1 (2021).Article 

Google Scholar 
Sullivan, R. M. Revision of the dinosaur Stegoceras Lambe (Ornithischia, Pachycephalosauridae). J. Vertebr. Paleontol. 23, 181–207. https://doi.org/10.1671/0272-4634(2003)23[181:ROTDSL]2.0.CO;2 (2003).Article 

Google Scholar 
Sullivan, R. & Lucas, S. The Kirtlandian land-vertebrate “age”-faunal composition, temporal position and biostratigraphic correlation in the nonmarine Upper Cretaceous of western North America. N. M. Mus. Nat. Hist. Sci. Bull. 35, 7–29 (2006).
Google Scholar 
Lucas, S. G., Sullivan, R. M., Lichtig, A., Dalman, S. & Jasinski, S. E. in Cretaceous Period: Biotic Diversity and Biogeography Vol. New Mexico Museum of Natural History and Science Bulletin 71 (eds S. G. Lucas & A. Khosla) 195–213 (2016).Dean, C. D., Chiarenza, A. A. & Maidment, S. C. R. Formation binning: a new method for increased temporal resolution in regional studies, applied to the Late Cretaceous dinosaur fossil record of North America. Palaeontology 63, 881–901. https://doi.org/10.1111/pala.12492 (2020).Article 

Google Scholar 
Maidment, S. C. R., Dean, C. D., Mansergh, R. I. & Butler, R. J. Deep-time biodiversity patterns and the dinosaurian fossil record of the Late Cretaceous Western Interior, North America. Proc. Biol. Sci. 288, 20210692. https://doi.org/10.1098/rspb.2021.0692 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Loughney, K. M. & Badgley, C. The influence of depositional environment and basin history on the taphonomy of mammalian assemblages from the Barstow Formation (middle Miocene), California. Palaios 35, 175–190. https://doi.org/10.2110/palo.2019.067 (2020).ADS 
Article 

Google Scholar 
Sakamoto, M., Benton, M. J. & Venditti, C. Dinosaurs in decline tens of millions of years before their final extinction. Proc. Natl. Acad. Sci. 113, 5036–5040. https://doi.org/10.1073/pnas.1521478113 (2016).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Condamine, F. L., Guinot, G., Benton, M. J. & Currie, P. J. Dinosaur biodiversity declined well before the asteroid impact, influenced by ecological and environmental pressures. Nat. Commun. 12, 3833. https://doi.org/10.1038/s41467-021-23754-0 (2021).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Therrien, F. O. & Fastovsky, D. E. Paleoenvironments of early theropods, Chinle Formation (Late Triassic), Petrified Forest National Park, Arizona. Palaios 15, 194–211. https://doi.org/10.1669/0883-1351(2000)015%3c0194:poetcf%3e2.0.co;2 (2000).ADS 
Article 

Google Scholar 
Hoke, G. D., Schmitz, M. D. & Bowring, S. A. An ultrasonic method for isolating nonclay components from clay-rich material. Geochem. Geophys. Geosyst. 15, 492–498. https://doi.org/10.1002/2013GC005125 (2014).ADS 
Article 

Google Scholar 
Ramezani, J. et al. High-precision U-Pb zircon geochronology of the Late Triassic Chinle Formation, Petrified Forest National Park (Arizona, USA): temporal constraints on the early evolution of dinosaurs. Geol. Soc. Am. Bull. 123, 2142–2159. https://doi.org/10.1130/b30433.1 (2011).ADS 
CAS 
Article 

Google Scholar 
Widmann, P., Davies, J. H. F. L. & Schaltegger, U. Calibrating chemical abrasion: its effects on zircon crystal structure, chemical composition and U-Pb age. Chem. Geol. 511, 1–10. https://doi.org/10.1016/j.chemgeo.2019.02.026 (2019).ADS 
CAS 
Article 

Google Scholar 
Krogh, T. E. Low-contamination method for hydrothermal decomposition of zircon and extraction of U and Pb for isotopic age determinations. Geochim. Cosmochim. Acta 37, 485–494. https://doi.org/10.1016/0016-7037(73)90213-5 (1973).ADS 
CAS 
Article 

Google Scholar 
Gerstenberger, H. & Haase, G. A highly effective emitter substance for mass spectrometric Pb isotope ratio determinations. Chem. Geol. 136, 309–312. https://doi.org/10.1016/S0009-2541(96)00033-2 (1997).ADS 
CAS 
Article 

Google Scholar 
Bowring, J. F., McLean, N. M. & Bowring, S. A. Engineering cyber infrastructure for U-Pb geochronology: Tripoli and U-Pb_Redux. Geochem. Geophys. Geosyst. https://doi.org/10.1029/2010gc003479 (2011).Article 

Google Scholar 
McLean, N. M., Bowring, J. F. & Bowring, S. A. An algorithm for U-Pb isotope dilution data reduction and uncertainty propagation. Geochem. Geophys. Geosyst. https://doi.org/10.1029/2010gc003478 (2011).Article 

Google Scholar 
Machlus, M. L. et al. A strategy for cross-calibrating U-Pb chronology and astrochronology of sedimentary sequences: an example from the Green River Formation, Wyoming, USA. Earth Planet. Sci. Lett. 413, 70–78. https://doi.org/10.1016/j.epsl.2014.12.009 (2015).ADS 
CAS 
Article 

Google Scholar 
Hiess, J., Condon, D. J., McLean, N. & Noble, S. R. 238U/235U systematics in terrestrial uranium-bearing minerals. Science 335, 1610–1614. https://doi.org/10.1126/science.1215507 (2012).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Schoene, B., Crowley, J. L., Condon, D. J., Schmitz, M. D. & Bowring, S. A. Reassessing the uranium decay constants for geochronology using ID-TIMS U-Pb data. Geochim. Cosmochim. Acta 70, 426–445. https://doi.org/10.1016/j.gca.2005.09.007 (2006).ADS 
CAS 
Article 

Google Scholar 
Mattinson, J. M. Analysis of the relative decay constants of 235U and 238U by multi-step CA-TIMS measurements of closed-system natural zircon samples. Chem. Geol. 275, 186–198. https://doi.org/10.1016/j.chemgeo.2010.05.007 (2010).ADS 
CAS 
Article 

Google Scholar 
Jaffey, A. H., Flynn, K. F., Glendenin, L. E., Bentley, W. C. & Essling, A. M. Precision measurement of half-lives and specific activities of 235U and 238U. Phys. Rev. C 4, 1889–1906. https://doi.org/10.1103/PhysRevC.4.1889 (1971).ADS 
Article 

Google Scholar 
Nasdala, L. et al. GZ7 and GZ8—two zircon reference materials for SIMS U-Pb geochronology. Geostand. Geoanal. Res. 42, 431–457. https://doi.org/10.1111/ggr.12239 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Haslett, J. & Parnell, A. A simple monotone process with application to radiocarbon-dated depth chronologies. J. R. Stat. Soc. C Appl. Stat. 57, 399–418. https://doi.org/10.1111/j.1467-9876.2008.00623.x (2008).MathSciNet 
Article 
MATH 

Google Scholar 
Parnell, A. C., Haslett, J., Allen, J. R. M., Buck, C. E. & Huntley, B. A flexible approach to assessing synchroneity of past events using Bayesian reconstructions of sedimentation history. Quat. Sci. Rev. 27, 1872–1885. https://doi.org/10.1016/j.quascirev.2008.07.009 (2008).ADS 
Article 

Google Scholar  More

Jiang, D., Klaus, S., Zhang, Y.-P., Hillis, D. M. & Li, J.-T. Asymmetric biotic interchange across the Bering land bridge between Eurasia and North America. Natl Sci. Rev. 6, 739–745 (2019).Article 

Google Scholar 
Bailey, G. et al. in The Archaeology of Europe’s Drowned Landscapes (eds Bailey, G. et al.) 189–219 (Springer, 2020).Erlandson, J., Braje, T., Gill, K. & Graham, M. Ecology of the kelp highway: did marine resources facilitate human dispersal from northeast Asia to the Americas? J. Isl. Coast. Archaeol. 10, 392–411 (2015).Article 

Google Scholar 
McLaren, D. et al. Terminal Pleistocene epoch human footprints from the Pacific coast of Canada. PLoS ONE 13, e0193522 (2018).Article 

Google Scholar 
Woodward, J. C. The Ice Age: A Very Short Introduction (Oxford Univ. Press, 2014).Pettitt, P. & White, M. The British Palaeolithic: Human Societies at the Edge of the Pleistocene World (Routledge, 2012).Bell, M. Prehistoric Coastal Communities: The Mesolithic in Western Britain (Council for British Archaeology, 2007).Gaffney, V., Fitch, S. & Smith, D. Europe’s Lost World: The Rediscovery of Doggerland (Council for British Archaeology, 2009).Gutiérrez–Zugasti, I. et al. Shell midden research in Atlantic Europe: state of the art, research problems and perspectives for the future. Quat. Int. 239, 70–85 (2011).Article 

Google Scholar 
Milner, N., Craig, O.E. & Bailey, G. N. Shell Middens in Atlantic Europe (Oxbow Books, 2007).Neto de Carvalho, C. et al. First tracks of newborn straight-tusked elephants (Palaeoloxodon antiquus). Sci. Rep. 11, 17311 (2021).CAS 
Article 

Google Scholar 
Ashton, N. et al. Hominin footprints from Early Pleistocene deposits at Happisburgh, UK. PLoS ONE 9, e88329 (2014).Article 

Google Scholar 
Duveau, J., Berillon, G. & Verna, C. in Reading Prehistoric Human Tracks: Methods and Material (eds Pastoors, A. & Tilman, L.) 183–200 (Springer, 2021).Allen, J. R. L. Subfossil mammalian tracks (Flandrian) in the Severn Estuary, S.W. Britain: mechanics of formation, preservation and distribution. Philos. Trans. R. Soc. B 352, 481–518 (1997).Article 

Google Scholar 
Aldhouse-Green, S. H. R. et al. Prehistoric footprints from the Severn Estuary at Uskmouth and Magor Pill, Gwent, Wales. Archaeol. Cambrensis 141, 14–55 (1992).
Google Scholar 
Barr, K. & Bell, M. Neolithic and Bronze Age ungulate footprint-tracks of the Severn Estuary: species, age, identification and the interpretation of husbandry practices. Environ. Archaeol. 22, 1–14 (2016).Article 

Google Scholar 
Bennett, M. R. & Morse, S. A. Human Footprints: Fossilised Locomotion? (Springer, 2014).Polton, J. A., Palmer, M. R. & Howarth, M. J. Physical and dynamical oceanography of Liverpool Bay. Ocean Dynam. 91, 1421–1439 (2011).Article 

Google Scholar 
Roberts, G. Ephemeral, subfossil mammalian, avian and hominid footprints within Flandrian sediment exposures at Formby Point, Sefton Coast, North West England. Ichnos 16, 33–48 (2009).Article 

Google Scholar 
Burns, A. An 8000-Year Record of Prehistoric Footprints in a Dynamic Coastal Landscape, Formby Point, UK. PhD thesis, Univ. Manchester (2019).Tooley, M. J. Sea level changes in Northern England. Proc. Geol. Assoc. 93, 43–51 (1982).Article 

Google Scholar 
Pye, K. & Neal, A. in The Dynamics and Environmental context of Aeolian Sedimentary Systems (ed. Pye, K.) 201–217 (Geological Society Special Publication, 1993).Pye, K., Stokes, S. & Neal, A. Optical dating of aeolian sediments from the Sefton Coast, Northwest England. Proc. Geol. Assoc. 106, 281–292 (1995).Article 

Google Scholar 
Reimer, P. J. et al. The IntCal20 Northern Hemisphere radiocarbon calibration curve (0–55 cal kBP). Radiocarbon 62, 725–757 (2020).CAS 
Article 

Google Scholar 
Walker, M. et al. Subdividing the Holocene Series/Epoch: formalisation of stages/ages and subseries/subepochs, and designation of GSSPs and auxiliary stratotypes. J. Quat. Sci. 34, 173–186 (2019).Article 

Google Scholar 
Tooley, M. J. The peat beds of the southwest Lancashire coast. Nat. Lancs. 1, 19–21 (1970).
Google Scholar 
Gonzalez, S., Huddart, D. & Roberts, G. Holocene development of the Sefton coast: a multidisciplinary approach to understanding the archaeology. In Archaeological Sciences 1995 Proc. Conference on the Application of Scientific Techniques to the Study of Archaeology (eds Sinclair, A. et al.) 289–299 (Oxbow Books, 1997).Huddart, D., Roberts, G. & Gonzalez, S. Holocene human and animal footprints and their relationships with coastal environmental change, Formby Point, NW England. Quat. Int. 55, 29–41 (1999).Article 

Google Scholar 
Galbraith, H. et al. Global climate change and sea level rise: potential losses of intertidal habitat for shorebirds. Waterbirds 25, 173–183 (2002).Article 

Google Scholar 
Bellard, C., Leclerc, C. & Courchamp, F. Sea level rise and insular hotspots. Glob. Ecol. Biogeogr. 23, 203–212 (2014).Article 

Google Scholar 
Editorial Why biodiversity matters. Nat. Ecol. Evol. 1, 0042 (2017).Article 

Google Scholar 
Hall, J. G. A comparative analysis of the habitat of the extinct aurochs and other prehistoric mammals in Britain. Ecography 31, 187–190 (2008).Article 

Google Scholar 
Milner, N., Conneller, C. & Taylor, B. Star Carr: A Persistent Place in a Changing World Vol. 1 (White Rose Univ. Press, 2018).Conneller, C. The Mesolithic in Britain: Landscape and Society in Times of Change (Routledge, 2022).Hernandez, L. & Laundre, J. W. Foraging in the ‘landscape of fear’ and its implications for habitat use and diet quality of elk (Cervus elaphus) and bison (Bison bison). Wildl. Biol. 11, 215–220 (2005).Article 

Google Scholar 
Overton, N. J. in Multispecies Archaeology (ed. Pilaar-Birch, S.) 295–309 (Routledge, 2018).Moore, E. K., Britton, A. J., Iason, G., Pemberton, J. & Pakeman, R. J. Landscape-scale vegetation patterns influence small-scale grazing impacts. Biol. Conserv. 192, 218–225 (2015).Article 

Google Scholar 
Gonzalez, S. & Huddart, D. in The Quaternary of Northern England (eds Huddart, D. & Glasser, N. F.) 582–588 (Joint Nature Conservation Committee, 2002).Cowell, R. W. & Innes, J. The Wetlands of Merseyside (Lancaster Univ. Archaeological Unit, 1994).Leonard, P. B. et al. Landscape connectivity losses due to sea level rise and land use change. Anim. Conserv. 20, 80–90 (2017).Article 

Google Scholar 
Myers, N. et al. Biodiversity hotspots for conservation priorities. Nature 403, 853–858 (2000).CAS 
Article 

Google Scholar 
Stuart, A. J. in Island Britain: a Quaternary Perspective (ed. Preece, R. C.) 111–125 (Geological Society Special Publication, 1995).Maroo, S. & Yalden, D. W. The Mesolithic mammal fauna of Great Britain. Mammal. Rev. 30, 243–248 (2000).Article 

Google Scholar 
Yalden, D. W. The History of British Mammals (T. & A.D. Poyser, 1999).Barnett, R. The Missing Lynx: The Past and Future of Britain’s Lost Mammals (Bloomsbury, 2019).Crees, J. J., Carbone, C., Sommer, R. S., Benecke, N. & Turvey, S. T. Millennial-scale faunal record reveals differential resilience of European large mammals to human impacts across the Holocene. Proc. R. Soc. B 283, 20152152 (2016).Article 

Google Scholar 
Burns, A. in Reading Prehistoric Human Tracks (eds Pastoors, A. & Lenssen-Erz, T.) 295–315 (Springer, 2021).Brown, R., Lawrence, M. & Pope, J. Animals: Tracks, Trails and Signs (Octopus Publishing, 2004).Donovan, S. K. Animal and bird tracks. Ichnos 16, 238–238 (2009).Article 

Google Scholar 
Roberts, G., Gonzalez, S. & Huddart, D. Intertidal Holocene footprints and their archaeological significance. Antiquity 70, 647–651 (1996).Article 

Google Scholar 
Scales, R. in Prehistoric Coastal Communities: The Mesolithic in Western Britain (ed. Bell, M.) 139–159 (Council for British Archaeology, 2007).Robbins, L. M. Estimating height and weight from size of footprints. J. Forensic Sci. 31, 143–152 (1986).CAS 
Article 

Google Scholar 
Stuiver, M. & Polach, H. A. Reporting of C-14 data—discussion. Radiocarbon 19, 355–363 (1977).Article 

Google Scholar  More

Countries have yet to agree to protect at least 30% of land, a crucial target proposed in the global biodiversity deal.Credit: Roberto Schmidt/AFP via Getty

A beleaguered global deal to save the environment got a financial boost last week when Germany announced that it was upping its funding for international biodiversity conservation to €1.5 billion (US$1.49 billion) a year — an increase of €0.87 billion — making it the largest national financial pledge yet to save nature. The announcement came at a 20 September meeting in New York City, where political leaders, businesses and conservation and Indigenous-rights groups came together to rally momentum and support ahead of the United Nations biodiversity summit in Montreal, Canada, in December.Conservationists welcomed the extra funding, but warned that other wealthy countries must also reach deeper into their pockets to ensure that nations agree on a new biodiversity agreement, called the Post-2020 Global Biodiversity Framework. Estimates suggest that an additional US$700 billion annually is needed to protect the environment.Concerns over insufficient financing for global biodiversity conservation have stalled negotiations and threaten to derail attempts to finalize a deal in Montreal. The forthcoming summit will be the 15th meeting of the Conference of the Parties (COP15) to the UN’s Convention on Biological Diversity.Announcing the new funds, German Chancellor Olaf Scholz said: “With this contribution, we want to send a strong signal for an ambitious outcome of the biodiversity COP-15.”Claire Blanchard, head of global advocacy at WWF, a conservation group, told Nature that the extra funding “is highly significant” and sends an important signal that rich countries are prepared to step up.But she adds: “More signals of this kind will be needed to create the environment conducive to constructive dialogue in the negotiation room.”Andrew Deutz, a specialist in biodiversity law and finance at the Nature Conservancy, a conservation group in Arlington, Virginia, says he expects further funding announcements to come in the run up to and at the COP15.Other pledgesSeveral key political leaders, including Justin Trudeau, Canada’s prime minister, echoed calls for rich nations to make urgent progress to secure the biodiversity deal. Trudeau urged countries to agree on two crucial targets proposed in the biodiversity framework, both to be met by 2030: to halt and reverse biodiversity loss, and to protect at least 30% of land and seas.The new funding was bolstered by other pledges and developments, including a promise from a partnership of some of the world’s wealthiest private philanthropic foundations and charities to add to the $5 billion they have already committed to conservation, if other countries promise more funds.The partnership — which includes the Bezos Earth Fund, an environmental fund financed by entrepreneur Jeff Bezos — has already spent around $1 billion of its promised financing over the past two years, says Cristián Samper, head of the Wildlife Conservation Society, a not-for-profit group. Samper was speaking on behalf of the partnership at the meeting in New York City.Frans Timmermans, vice-president of the European Commission, reaffirmed that Europe would double its international biodiversity funding to $1.13 billion annually — a promise originally announced in September last year. Timmermans told the meeting that the European Union would set out more details about the funding soon.Funding shortfallAlso at the meeting, a group of four countries comprising Ecuador, Gabon, the Maldives and the United Kingdom launched a joint 10-point plan to bridge the biodiversity finance gap, which is estimated at $700 billion annually.The plan sets out the financial commitments and policy reforms needed to finance biodiversity on the required scale. For example, it encourages wealthy and lower-income nations to allocate new funds for biodiversity and to quickly deliver on their existing financial pledges. It requires donor countries to ensure that funds for overseas development do no harm to biodiversity. And it asks countries to dedicate a portion of their national funding for climate change to activities that also protect and conserve nature.The plan also commits countries to ensuring that public finance is invested in ways that benefit biodiversity, and to reviewing national subsidies and redirecting those that are harmful to nature. It calls on businesses to assess and disclose commercial risks associated with biodiversity decline, and to set quantitative targets to reduce their impact on the natural world. And it encourages multilateral development banks — such as the World Bank in Washington DC — and international financial institutions to ensure that their investments benefit biodiversity, and asks that they report on their biodiversity funding in time for COP-15.So far, 15 countries, including Canada, Germany and Norway, as well as the EU have endorsed the plan.“The plan provides a clear pathway for bridging the global biodiversity finance gap. Its significance lies in the political signal it sends,” says Blanchard.António Guterres, secretary-general of the UN, urged political leaders to “act now and at scale” to secure biodiversity financing and ensure agreement on the framework. “If negotiations continue at their slow pace, we are headed to failure,” he told the meeting. More

Gardner, T. A., Côté, I. M., Gill, J. A., Grant, A. & Watkinson, A. R. Long-term region-wide declines in Caribbean corals. Science 301, 958–960 (2003).ADS 
PubMed 
Article 

Google Scholar 
Bruno, J. F. & Selig, E. R. Regional decline of coral cover in the Indo-Pacific: Timing, extent, and subregional comparisons. PLoS ONE 2, e711 (2007).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
De’Ath, G., Fabricius, K. E., Sweatman, H. & Puotinen, M. The 27-year decline of coral cover on the Great Barrier Reef and its causes. Proc. Natl. Acad. Sci. U. S. A. 109, 17995–17999 (2012).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hughes, T. P., Graham, N. A. J., Jackson, J. B. C., Mumby, P. J. & Steneck, R. S. Rising to the challenge of sustaining coral reef resilience. Trends Ecol. Evol. 25, 633–642 (2010).PubMed 
Article 

Google Scholar 
Pandolfi, J. M., Connolly, S. R., Marshall, D. J. & Cohen, A. L. Projecting coral reef futures under global warming and ocean acidification. Science 333, 418–422 (2011).ADS 
PubMed 
Article 

Google Scholar 
Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377 (2017).ADS 
PubMed 
Article 

Google Scholar 
Bindoff, N. L. et al. Chapter 5: Changing Ocean, Marine Ecosystems, and Dependent Communities. IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (2019).Richmond, R. H. Reproduction and recruitment in corals: Critical links in the persistence of reefs. In Life and Death of Coral Reefs (ed. Birkeland, C. E.) 175–197 (Springer, 1997).Chapter 

Google Scholar 
Trapon, M. L., Pratchett, M. S., Hoey, A. S. & Baird, A. H. Influence of fish grazing and sedimentation on the early post-settlement survival of the tabular coral Acropora cytherea. Coral Reefs 32, 1051–1059 (2013).ADS 
Article 

Google Scholar 
Gallagher, C. & Doropoulos, C. Spatial refugia mediate juvenile coral survival during coral–predator interactions. Coral Reefs 36, 51–61 (2017).ADS 
Article 

Google Scholar 
Vermeij, M. J. A. & Sandin, S. A. Density-dependent settlement and mortality structure the earliest life phases of a coral population. Ecology 89, 1994–2004 (2008).PubMed 
Article 

Google Scholar 
Vermeij, M. J. A., Smith, J. E., Smith, C. M., Vega Thurber, R. & Sandin, S. A. Survival and settlement success of coral planulae: Independent and synergistic effects of macroalgae and microbes. Oecologia 159, 325–336 (2009).ADS 
PubMed 
Article 

Google Scholar 
Ricardo, G. F., Jones, R. J., Nordborg, M. & Negri, A. P. Settlement patterns of the coral Acropora millepora on sediment-laden surfaces. Sci. Total Environ. 609, 277–288 (2017).ADS 
PubMed 
Article 

Google Scholar 
Brunner, C. A., Uthicke, S., Ricardo, G. F., Hoogenboom, M. O. & Negri, A. P. Climate change doubles sedimentation-induced coral recruit mortality. Sci. Total Environ. 768, 143897 (2021).ADS 
PubMed 
Article 

Google Scholar 
Birrell, C. L., McCook, L. J., Willis, B. L. & Diaz-Pulido, G. A. Effects of benthic algae on the replenishment of corals and the implications for the resilience of coral reefs. In Oceanography and Marine Biology: An Annual Review 25–63 (CRC Press, 2008).Chapter 

Google Scholar 
Karcher, D. B. et al. Nitrogen eutrophication particularly promotes turf algae in coral reefs of the central Red Sea. PeerJ 2020, 1–25 (2020).
Google Scholar 
Kirschner, C. M. & Brennan, A. B. Bio-inspired antifouling strategies. Annu. Rev. Mater. Res. 42, 211–229 (2012).ADS 
Article 

Google Scholar 
Webster, N. S. et al. Metamorphosis of a scleractinian coral in response to microbial biofilms. Appl. Environ. Microbiol. 70, 1213–1221 (2004).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Heyward, A. J. & Negri, A. P. Natural inducers for coral larval metamorphosis. Coral Reefs 18, 273–279 (1999).Article 

Google Scholar 
Negri, A. P., Webster, N. S., Hill, R. T. & Heyward, A. J. Metamorphosis of broadcast spawning corals in response to bacteria isolated from crustose algae. Mar. Ecol. Prog. Ser. 223, 121–131 (2001).ADS 
Article 

Google Scholar 
Tebben, J. et al. Induction of larval metamorphosis of the coral Acropora millepora by tetrabromopyrrole isolated from a Pseudoalteromonas bacterium. PLoS ONE 6, 1–8 (2011).Article 

Google Scholar 
Sneed, J. M., Sharp, K. H., Ritchie, K. B. & Paul, V. J. The chemical cue tetrabromopyrrole from a biofilm bacterium induces settlement of multiple Caribbean corals. Proc. R. Soc. B Biol. Sci. 281, 1–9 (2014).
Google Scholar 
Tebben, J. et al. Chemical mediation of coral larval settlement by crustose coralline algae. Sci. Rep. 5, 1–11 (2015).Article 

Google Scholar 
Carpenter, R. C. & Edmunds, P. J. Local and regional scale recovery of Diadema promotes recruitment of scleractinian corals. Ecol. Lett. 9, 268–277 (2006).Article 

Google Scholar 
Box, S. J. & Mumby, P. J. Effect of macroalgal competition on growth and survival of juvenile Caribbean corals. Mar. Ecol. Prog. Ser. 342, 139–149 (2007).ADS 
Article 

Google Scholar 
Linares, C., Cebrian, E. & Coma, R. Effects of turf algae on recruitment and juvenile survival of gorgonian corals. Mar. Ecol. Prog. Ser. 452, 81–88 (2012).ADS 
Article 

Google Scholar 
McCook, L. J., Jompa, J. & Diaz-Pulido, G. Competition between corals and algae on coral reefs: A review of evidence and mechanisms. Coral Reefs 19, 400–417 (2001).ADS 
Article 

Google Scholar 
Nugues, M. M., Smith, G. W., Van Hooidonk, R. J., Seabra, M. I. & Bak, R. P. M. Algal contact as a trigger for coral disease. Ecol. Lett. 7, 919–923 (2004).Article 

Google Scholar 
Fong, J. et al. Allelopathic effects of macroalgae on Pocillopora acuta coral larvae. Mar. Environ. Res. 151, 104745. https://doi.org/10.1016/j.marenvres.2019.06.007 (2019).Article 
PubMed 

Google Scholar 
Hauri, C., Fabricius, K. E., Schaffelke, B. & Humphrey, C. Chemical and physical environmental conditions underneath mat- and canopy-forming macroalgae, and their effects on understorey corals. PLoS ONE 5, 1–9 (2010).Article 

Google Scholar 
Bay, L. K. et al. Reef Restoration and Adaptation Program : Intervention Technical Summary. A report provided to the Australian Government by the Reef Restoration and Adaptation Program. (2019).Anthony, K. R. N. et al. Interventions to help coral reefs under global change—A complex decision challenge. PLoS ONE 15, 1–14 (2020).Article 

Google Scholar 
Vardi, T. et al. Six priorities to advance the science and practice of coral reef restoration worldwide. Restor. Ecol. 29, 1–7 (2021).Article 

Google Scholar 
Heyward, A. J., Rees, M. & Smith, L. D. Coral spawning slicks harnessed for large-scale coral culture. Progr. Abstr. Int. Conf. Sci. Asp. Coral Reef Assess. Monit. Restor. 104, 188–189 (1999).
Google Scholar 
Harrison, P., Villanueva, R. & De la Cruz, D. Coral Reef Restoration using Mass Coral Larval Reseeding (Southern Cross University, 2016).
Google Scholar 
de la Cruz, D. W. & Harrison, P. L. Enhancing coral recruitment through assisted mass settlement of cultured coral larvae. PLoS ONE 15, e0242847. https://doi.org/10.1371/journal.pone.0242847 (2020).Article 

Google Scholar 
Chamberland, V. F., Snowden, S., Marhaver, K. L., Petersen, D. & Vermeij, M. J. A. The reproductive biology and early life ecology of a common Caribbean brain coral, Diploria labyrinthiformis (Scleractinia: Faviinae). Coral Reefs 36, 83–94 (2017).ADS 
Article 

Google Scholar 
Randall, C. J. et al. Sexual production of corals for reef restoration in the Anthropocene. Mar. Ecol. Prog. Ser. 635, 203–232 (2020).ADS 
Article 

Google Scholar 
Miller, M. W. et al. Settlement yields in large-scale in situ culture of Caribbean coral larvae for restoration. Restor. Ecol. https://doi.org/10.1111/rec.13512 (2021).Article 

Google Scholar 
Baria-Rodriguez, M. V., de la Cruz, D. W., Dizon, R. M., Yap, H. T. & Villanueva, R. D. Performance and cost-effectiveness of sexually produced Acropora granulosa juveniles compared with asexually generated coral fragments in restoring degraded reef areas. Aquat. Conserv. Mar. Freshw. Ecosyst. 29, 891–900 (2019).Article 

Google Scholar 
Doropoulos, C., Elzinga, J., ter Hofstede, R., van Koningsveld, M. & Babcock, R. C. Optimizing industrial-scale coral reef restoration: Comparing harvesting wild coral spawn slicks and transplanting gravid adult colonies. Restor. Ecol. 27, 758–767 (2019).Article 

Google Scholar 
Kuffner, I. B., Andersson, A. J., Jokiel, P. L., Rodgers, K. S. & MacKenzie, F. T. Decreased abundance of crustose coralline algae due to ocean acidification. Nat. Geosci. 1, 114–117 (2008).ADS 
Article 

Google Scholar 
Webster, N. S., Uthicke, S., Botté, E. S., Flores, F. & Negri, A. P. Ocean acidification reduces induction of coral settlement by crustose coralline algae. Glob. Change Biol. 19, 303–315 (2013).ADS 
Article 

Google Scholar 
Randall, C. J., Giuliano, C., Heyward, A. J. & Negri, A. P. Enhancing coral survival on deployment devices with microrefugia. Front. Mar. Sci. 8, 662263. https://doi.org/10.3389/fmars.2021.662263 (2021).Article 

Google Scholar 
Kuffner, I. B. et al. Inhibition of coral recruitment by macroalgae and cyanobacteria. Mar. Ecol. Prog. Ser. 323, 107–117 (2006).ADS 
Article 

Google Scholar 
Arnold, S. N., Steneck, R. S. & Mumby, P. J. Running the gauntlet: Inhibitory effects of algal turfs on the processes of coral recruitment. Mar. Ecol. Prog. Ser. 414, 91–105 (2010).ADS 
Article 

Google Scholar 
Speare, K. E., Duran, A., Miller, M. W. & Burkepile, D. E. Sediment associated with algal turfs inhibits the settlement of two endangered coral species. Mar. Pollut. Bull. 144, 189–195 (2019).PubMed 
Article 

Google Scholar 
Tebben, J., Guest, J. R., Sin, T. M., Steinberg, P. D. & Harder, T. Corals like it waxed: Paraffin-based antifouling technology enhances coral spat survival. PLoS ONE 9, 1–8 (2014).Article 

Google Scholar 
Almeida, E., Diamantino, T. C. & de Sousa, O. Marine paints: The particular case of antifouling paints. Prog. Org. Coat. 59, 2–20 (2007).Article 

Google Scholar 
Negri, A. P., Smith, L. D., Webster, N. S. & Heyward, A. J. Understanding ship-grounding impacts on a coral reef: Potential effects of anti-foulant paint contamination on coral recruitment. Mar. Pollut. Bull. 44, 111–117 (2002).PubMed 
Article 

Google Scholar 
Smith, L. D., Negri, A. P., Philipp, E., Webster, N. S. & Heyward, A. J. The effects of antifoulant-paint-contaminated sediments on coral recruits and branchlets. Mar. Biol. 143, 651–657 (2003).Article 

Google Scholar 
Jacobson, A. H. & Willingham, G. L. Sea-nine antifoulant: An environmentally acceptable alternative to organotin antifoulants. Sci. Total Environ. 258, 103–110 (2000).ADS 
PubMed 
Article 

Google Scholar 
Silva, V. et al. Isothiazolinone biocides: Chemistry, biological, and toxicity profiles. Molecules 25, 991. https://doi.org/10.3390/molecules25040991 (2020).Article 
PubMed Central 

Google Scholar 
da Silva, A. R., da Guerreiro, A. S., Martins, S. E. & Sandrini, J. Z. DCOIT unbalances the antioxidant defense system in juvenile and adults of the marine bivalve Amarilladesma mactroides (Mollusca: Bivalvia). Comp. Biochem. Physiol. Part C 250, 109169 (2021).
Google Scholar 
Cima, F. et al. Preliminary evaluation of the toxic effects of the antifouling biocide Sea-Nine 211TM in the soft coral Sarcophyton cf. glaucum (Octocorallia, Alcyonacea) based on PAM fluorometry and biomarkers. Mar. Environ. Res. 83, 16–22 (2013).PubMed 
Article 

Google Scholar 
Wendt, I., Backhaus, T., Blanck, H. & Arrhenius, Å. The toxicity of the three antifouling biocides DCOIT, TPBP and medetomidine to the marine pelagic copepod Acartia tonsa. Ecotoxicology 25, 871–879 (2016).PubMed 
Article 

Google Scholar 
Chen, L. et al. Identification of molecular targets for 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) in teleosts: New insight into mechanism of toxicity. Environ. Sci. Technol. 51, 1840–1847 (2017).ADS 
PubMed 
Article 

Google Scholar 
Martins, S. E., Fillmann, G., Lillicrap, A. & Thomas, K. V. Review: Ecotoxicity of organic and organo-metallic antifouling co-biocides and implications for environmental hazard and risk assessments in aquatic ecosystems. Biofouling 34, 34–52 (2018).PubMed 
Article 

Google Scholar 
Moon, Y. S., Kim, M., Hong, C. P., Kang, J. H. & Jung, J. H. Overlapping and unique toxic effects of three alternative antifouling biocides (Diuron, Irgarol 1051 ®, Sea-Nine 211 ® ) on non-target marine fish. Ecotoxicol. Environ. Saf. 180, 23–32 (2019).PubMed 
Article 

Google Scholar 
Su, Y. et al. Toxicity of 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) in the marine decapod Litopenaeus vannamei. Environ. Pollut. 251, 708–716 (2019).PubMed 
Article 

Google Scholar 
Fonseca, V. B., da Guerreiro, A. S., Vargas, M. A. & Sandrini, J. Z. Effects of DCOIT (4,5-dichloro-2-octyl-4-isothiazolin-3-one) to the haemocytes of mussels Perna perna. Comp. Biochem. Physiol Part C 232, 108737. https://doi.org/10.1016/j.cbpc.2020.108737 (2020).Article 

Google Scholar 
Ferreira, V. et al. Effects of nanostructure antifouling biocides towards a coral species in the context of global changes. Sci. Total Environ. 799, 149324 (2021).ADS 
PubMed 
Article 

Google Scholar 
de Campos, B. G. et al. A preliminary study on multi-level biomarkers response of the tropical oyster Crassostrea brasiliana to exposure to the antifouling biocide DCOIT. Mar. Pollut. Bull. 174, 112141 (2022).Article 

Google Scholar 
Maia, F. et al. Incorporation of biocides in nanocapsules for protective coatings used in maritime applications. Chem. Eng. J. 270, 150–157 (2015).Article 

Google Scholar 
Santos, J. V. N. et al. Can encapsulation of the biocide DCOIT affect the anti-fouling efficacy and toxicity on tropical bivalves?. Appl. Sci. 10, 1–12 (2020).Article 

Google Scholar 
Detty, M. R., Ciriminna, R., Bright, F. V. & Pagliaro, M. Environmentally benign sol-gel antifouling and foul-releasing coatings. Acc. Chem. Res. 47, 678–687 (2014).PubMed 
Article 

Google Scholar 
Korschelt, K., Tahir, M. N. & Tremel, W. A Step into the future: Applications of nanoparticle enzyme mimics. Chemistry 24, 9703–9713 (2018).PubMed 
Article 

Google Scholar 
Herget, K. et al. Haloperoxidase mimicry by CeO2-x nanorods combats biofouling. Adv. Mater. 29, 1–8 (2017).Article 

Google Scholar 
Korschelt, K. et al. CeO2-: X nanorods with intrinsic urease-like activity. Nanoscale 10, 13074–13082 (2018).PubMed 
Article 

Google Scholar 
Herget, K., Frerichs, H., Pfitzner, F., Tahir, M. N. & Tremel, W. Functional enzyme mimics for oxidative halogenation reactions that combat biofilm formation. Adv. Mater. 30, 1–28 (2018).Article 

Google Scholar 
Doropoulos, C., Ward, S., Marshell, A., Diaz-Pulido, G. & Mumby, P. J. Interactions among chronic and acute impacts on coral recruits: The importance of size-escape thresholds. Ecology 93, 2131–2138 (2012).PubMed 
Article 

Google Scholar 
Ji, Z. et al. Designed synthesis of CeO2 nanorods and nanowires for studying toxicological effects of high aspect ratio nanomaterials. ACS Nano 6, 5366–5380 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Herget, K. et al. Supporting Information: Haloperoxidase mimicry by CeO2-x nanorods combats biofouling. Adv. Mater. 29, 1603823 (2017).Article 

Google Scholar 
Sokolova, A. et al. Spontaneous multiscale phase separation within fluorinated xerogel coatings for fouling-release surfaces. Biofouling 28, 143–157 (2012).PubMed 
Article 

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

Google Scholar 
ImageJ Release Notes. https://imagej.nih.gov/ij/notes.html.Arganda-Carreras, I. et al. Trainable Weka Segmentation: A machine learning tool for microscopy pixel classification. Bioinformatics 33, 2424–2426 (2017).PubMed 
Article 

Google Scholar 
Arganda-Carreras, I. et al. Supplementary Data: Trainable Weka Segmentation: A Machine Learning Tool for Microscopy Pixel Classification: Trainable Weka Segmentation User Manualhttps://doi.org/10.1093/bioinformatics/btx180 (2017).Vyas, N., Sammons, R. L., Addison, O., Dehghani, H. & Walmsley, A. D. A quantitative method to measure biofilm removal efficiency from complex biomaterial surfaces using SEM and image analysis. Sci. Rep. 6, 2–11 (2016).Article 

Google Scholar 
Carbone, D. A., Gargano, I., Pinto, G., De Natale, A. & Pollio, A. Evaluating microalgae attachment to surfaces: A first approach towards a laboratory integrated assessment. Chem. Eng. Trans. 57, 73–78 (2017).
Google Scholar 
Moreno Osorio, J. H. et al. Early colonization stages of fabric carriers by two Chlorella strains. J. Appl. Phycol. 32, 3631–3644 (2020).Article 

Google Scholar 
Ricardo, G. F. et al. Impacts of water quality on Acropora coral settlement: The relative importance of substrate quality and light. Sci. Total Environ. 777, 146079. https://doi.org/10.1016/j.scitotenv.2021.146079 (2021).ADS 
Article 
PubMed 

Google Scholar 
Macadam, A., Nowell, C. J. & Quigley, K. Machine learning for the fast and accurate assessment of fitness in coral early life history. Remote Sens. 13, 1–17 (2021).Article 

Google Scholar 
Negri, A. P. & Heyward, A. J. Inhibition of Fertilization and Larval Metamorphosis of the Coral Acropora millepora (Ehrenberg, 1834) by Petroleum Products. Mar. Pollut. Bull. 41, 420–427 (2000).Article 

Google Scholar 
Nordborg, F. M., Flores, F., Brinkman, D. L., Agustí, S. & Negri, A. P. Phototoxic effects of two common marine fuels on the settlement success of the coral Acropora tenuis. Sci. Rep. 8, 1–12 (2018).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. (2021).Wickham, H. et al. Welcome to the Tidyverse. J. Open Source Softw. 4, 1686. https://doi.org/10.21105/joss.01686 (2019).ADS 
Article 

Google Scholar 
Pinheiro, J., Bates, D., Debroy, S., Sarkar, D. & R Core Team. Linear and nonlinear mixed effects models contact. Linear nonlinear Mix. Eff. Model. 3, 103–135 (2021).
Google Scholar 
Fox, J. & Weisberg, S. An R Companion to Applied Regression (Sage Publications, 2019).
Google Scholar 
Lenth, R. V. Emmeans: Estimated Marginal Means. https://cran.r-project.org/package=emmeans (2021).Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378–400 (2017).Article 

Google Scholar 
Dafforn, K. A., Lewis, J. A. & Johnston, E. L. Antifouling strategies: History and regulation, ecological impacts and mitigation. Mar. Pollut. Bull. 62, 453–465 (2011).PubMed 
Article 

Google Scholar 
Wu, R. et al. Room temperature synthesis of defective cerium oxide for efficient marine anti-biofouling. Adv. Compos. Hybrid Mater. https://doi.org/10.1007/s42114-021-00256-7 (2021).Article 

Google Scholar 
Hu, M. et al. Nanozymes in nanofibrous mats with haloperoxidase-like activity to combat biofouling. ACS Appl. Mater. Interfaces 10, 44722–44730 (2018).PubMed 
Article 

Google Scholar 
He, X. et al. Haloperoxidase mimicry by CeO2-x nanorods of different aspect ratios for antibacterial performance. ACS Sustain. Chem. Eng. 8, 6744–6752 (2020).Article 

Google Scholar 
Saxena, P. & Harish,. Nanoecotoxicological reports of engineered metal oxide nanoparticles on algae. Curr. Pollut. Rep. 4, 128–142 (2018).Article 

Google Scholar 
Xu, Y. et al. Effects of cerium oxide nanoparticles on bacterial growth and behaviors: Induction of biofilm formation and stress response. Environ. Sci. Pollut. Res. 26, 9293–9304 (2019).Article 

Google Scholar 
Xu, Y. et al. Mechanistic understanding of cerium oxide nanoparticle-mediated biofilm formation in Pseudomonas aeruginosa. Environ. Sci. Pollut. Res. 25, 34765–34776 (2018).Article 

Google Scholar 
Tang, Y. et al. Hybrid xerogel films as novel coatings for antifouling and fouling release. Biofouling 21, 59–71 (2005).PubMed 
Article 

Google Scholar 
Gunari, N. et al. The control of marine biofouling on xerogel surfaces with nanometer-scale topography. Biofouling 27, 137–149 (2011).PubMed 
Article 

Google Scholar 
Maia, F. et al. Silica nanocontainers for active corrosion protection. Nanoscale 4, 1287–1298 (2012).ADS 
PubMed 
Article 

Google Scholar 
Martins, R. et al. Effects of a novel anticorrosion engineered nanomaterial on the bivalve: Ruditapes philippinarum. Environ. Sci. Nano 4, 1064–1076 (2017).Article 

Google Scholar 
Gutner-Hoch, E. et al. Antimacrofouling efficacy of innovative inorganic nanomaterials loaded with booster biocides. J. Mar. Sci. Eng. 6, 15. https://doi.org/10.3390/jmse6010006 (2018).Article 

Google Scholar 
Negri, A. P. & Heyward, A. J. Inhibition of coral fertilisation and larval metamorphosis by tributyltin and copper. Mar. Environ. Res. 51, 17–27 (2001).PubMed 
Article 

Google Scholar 
Morse, D. E., Hooker, N., Morse, A. N. C. & Jensen, R. A. Control of larval metamorphosis and recruitment in sympatric agariciid corals. J. Exp. Mar. Biol. Ecol. 116, 193–217 (1988).Article 

Google Scholar 
Harrington, L., Fabricius, K., De’ath, G. & Negri, A. Recognition and selection of settlement substrata determine post-settlement survival in corals. Ecology 85, 3428–3437 (2004).Article 

Google Scholar 
Jorissen, H., Baumgartner, C., Steneck, R. S. & Nugues, M. M. Contrasting effects of crustose coralline algae from exposed and subcryptic habitats on coral recruits. Coral Reefs 39, 1767–1778 (2020).Article 

Google Scholar 
Figueiredo, J. et al. Toxicity of innovative anti-fouling nano-based solutions to marine species. Environ. Sci. Nano 6, 1418–1429 (2019).Article 

Google Scholar 
Shafir, S., Abady, S. & Rinkevich, B. Improved sustainable maintenance for mid-water coral nursery by the application of an anti-fouling agent. J. Exp. Mar. Biol. Ecol. 368, 124–128 (2009).Article 

Google Scholar  More

Study site and animalsWe studied n = 64 flocks of Japanese tits in mixed deciduous-coniferous forests in Nagano and Gumma (36°17-31’N, 138°26-39’E), Japan. Although most of the birds had not been individually colour-ringed, all the experimental trials were conducted at least 400 m apart; previous observations on colour-ringed individuals showed that this distance was enough to ensure the collection of data from different individuals30. In this site, one of the major predators of small birds is the bull-headed shrike, which is often mobbed by small birds including Japanese tits.Playback stimulusTo test whether Japanese tits recognize an alert-recruitment call sequence as a single unit, we prepared four treatments: (i) one-speaker playback of alert-recruitment call sequences, (ii) two-speaker playback of alert-recruitment call sequences with alert and recruitment calls played from different speakers, (iii) one-speaker playback of recruitment-alert call sequences, (iv) two-speaker playback of recruitment-alert call sequences with recruitment and alert calls played from different speakers (Fig. 3). We created sound files for these treatments using the software program Audacity 2.1.3 (http://www.audacityteam.org). For one-speaker treatments, we composed mono sound files where call sequences were repeated onto a single channel, whereas for two-speaker treatments, we composed stereo sound files where either alert or recruitment calls were repeated onto the right or left channels, respectively. All the files contained an equal number of alert calls (30 calls) and recruitment calls (30 calls) at the same rate (one call every 3 s), resulting in 90-s of stimuli (Fig. 3), which corresponds to the range of the natural calling rate of alert-recruitment sequences during mobbing by Japanese tits10. For all stimuli, within-call-sequence intervals between alert and recruitment calls were constant (0.1 s), which is within the range of intervals of these calls in natural call sequences11,17. In contrast, between-call-sequence intervals varied from 1.50 to 1.81 (median = 1.68) due to the difference in call length, but were constant across playback stimuli within the same “block” where the four treatments were created using the same call exemplars (see below). While alert calls are composed of three distinct note types, recruitment calls are strings of the same note type that vary in repetition number. Since the repetition number can vary depending on predator type10, we conducted predator exposure experiments to Japanese tit flocks (n = 12) and recorded call sequences towards a bull-headed shrike life-like specimen. In response to a shrike specimen, tits produced alert-recruitment call sequences with a recruitment note repetition number ranging from 5 to 15. Since the interquartile range of repetition number was 6.75 to 10, we used recruitment calls with 7–10 notes as playback stimuli in this study. In consideration for the possible influence of sound editing procedure, we created all the stimuli in the same manner; we copied alert and recruitment call parts separately from recording files, and pasted them onto background noise files to produce all four types of stimuli. Playback amplitudes were constant across treatments, 70 dB at 1.0 m measured using a sound level meter (SM-325, AS ONE Corporation). Therefore, the differences between treatments only depend on whether these calls are produced as sequences from the same source and how the calls are ordered.We carefully designed experiments to control for the possibility that individual-based acoustic features in alert and recruitment calls might influence tits’ responses. First, we prepared 16 unique sets of alert and recruitment calls using either calls from the same bird (n = 8 source individuals, n = 8 unique call sets) or from two different birds (n = 16 source individuals, n = 8 unique call sets). Then, we created the four types of treatments (i.e., alert-recruitment call sequences from the same speaker, from different speakers, and in reversed order from the same speaker and from different speakers) from each of the alert-recruitment call sets, resulting in 16 blocks of playback stimuli (Supplementary Table 3). This allows us to test the possible influence of individual-based acoustic variation on receivers’ responses.We were also careful to avoid the possible influence of population-level signatures of acoustic features: we only used Japanese tits’ call sequences that had been previously recorded from the same study population. We saved the sound files in .wav format (16-bit accuracy, 48-kHz sampling rate) onto a playback device (iPhone 8, Apple Inc.). We used the default Music app (Apple Inc.) to playback the sound files.ExperimentWe (TNS and YKM) conducted experimental trials from 26 October to 4 December 2020 and during the period of 0800 and 1600 h (Japan Standard Time). We did not conduct trials under wet and windy weather conditions, since these may influence behavioural patterns of forest birds31. First, we searched for and located a flock of Japanese tits. Upon finding a flock, we fixed a taxidermic specimen of bull-headed shrike in a perching posture on the branch at 1.8 ± 0.2 m (mean ± s.d., n = 64) above the ground. Then, we placed either one or two Bluetooth speakers (SoundLink Micro, BOSE) on tree branches at 1.6 ± 0.2 m (mean ± s.d., n = 96) above the ground, and oriented them upwards to control for the possible influence of directionality. We set the distance between the shrike specimen and the speaker(s) at 5 m. For trials with two speakers, we set the distance between speakers at 10 m, mimicking the situation in which two birds are calling (Fig. 3). The shrike specimen was first covered with a black cloth and was exposed by removing the cloth just before each trial.We began playbacks when at least two Japanese tits were present within 15 m from the shrike specimen. During 90-s of playbacks, we recorded (i) whether birds approached within 2-m of the shrike specimen during the playback and (ii) whether birds exhibited wing flicking displays12,13. We counted the number of individuals within 15 m from the shrike during 90-s of playbacks and considered it as flock size. During trials, we sat on the ground at ca. 10 m from the shrike specimen to decrease the influence of the observers’ presence on bird behaviour. To account for the inter-observer reliability32, we calculated intra-class correlation coefficient (ICC; icc function in the R package irr) between us. The lowest ICC was 0.998, indicating high degree of inter-observer reliability for the two behavioural measurements. We also video-recorded the responses of tits using a digital video camera (FDR-AX60, SONY). After completion of each trial, we checked the video recording and made an on-the-spot confirmation of the exact location at which each bird made the closest approach to the shrike specimen during the 90-s of playbacks. Then, using a tape measure, we recorded the minimum approach distance of birds to the shrike specimen. Thus, our final data set consisted of the most reliable observations confirmed by two experimenters and video evidence.The order of trials was randomized within each block (n = 16 blocks), each of which is composed of a unique alert-recruitment call set but includes four treatments differing in the number of speakers and call order. Therefore, responses to all four treatments were observed under largely similar conditions. In a few trials, the first bird to approach the shrike specimen was from a heterospecific species, such as a varied tit (n = 1) or a long-tailed tit (n = 1). To account for the possibility that these birds evoke mobbing behaviour in Japanese tits, we only used the data from instances where the first individual to approach the shrike was a Japanese tit. Otherwise, we repeated the same treatment at a different site.We used 64 unique playbacks created from 16 unique sets of alert-recruitment calls for 64 trials in order to avoid pseudoreplication33. We prepared two specimens of male bull-headed shrikes and used each of them for the equal number of trials. We did not use specimens of female shrikes since females migrate from the study site in late summer and only males were observed during the study period.Statistical analysisWe analyzed the effect of playback treatments on the mobbing behaviours of Japanese tits using generalized linear mixed models in R34,35. We used the proportions of Japanese tits in flocks that (i) approached within 2-m of the shrike specimen and (ii) exhibited wing flicking displays. For the analysis of predator approach, we prepared two vectors (i.e., the number of Japanese tits that approached the shrike specimen and the number of Japanese tits that did not approach the shrike specimen). Then, we created a single response variable by binding together these two vectors using cbind function. Similarly, for the analysis of wing flicking displays, we created a single response variable by binding two vectors (i.e., the number of tits that exhibit wing flicking and the number of tits that did not exhibit wing flicking). We fitted playback treatments as a fixed term, and flock size (maximum number of Japanese tits observed during 90-s of playback) and the way of creating playback stimuli (whether the two call types were recorded from a single individual or two individuals) as covariates. We also included identity of alert-recruitment call sets that were used for creating playback stimuli (i.e., call sets from either one or two source individuals) and identity of shrike specimens as random terms. We used a binomial error distribution and logit-link function (glmer in the R package lme4) for these models. Statistical significance was calculated by log-likelihood ratio tests using anova in the R package stats. We further conducted post-hoc pairwise comparisons between treatments by using estimated marginal means (emmeans in the R package emmeans). When making pairwise comparisons, we adjusted p-values by applying a false discovery rate control for multiple testing36. All tests were two-sided and the significance level was set at α = 0.05. Exact p-values are reported when p ≥ 0.0001.Ethics statementAll protocols were approved by the ethics committee of Kyoto University, the Ministry of the Environment, and the Forestry Agency of Japan, and adhered to Guidelines for the Use of Animals of the Association for the Study of Animal Behaviour/Animal Behavior Society37.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Experimental designA 2-year field experiment was conducted at the Modern Agricultural Research and Development Base of Henan Province (113° 35′–114° 15′ E, 34° 53′–35° 11′ N). In order to enhance the diversity of LAI data, a split-plot design with a variety of field management measures and three replications was selected for the experiment (Fig. 1). The size of each experiment plot was 40 m2, the soil texture was predominantly sandy loam and sandy clay loam, as determined by textural analysis of soil samples collected before planting. Maize cultivar Dedan-5 was used in the experiment, which was planted on June 12, 2019, and June 20, 2020, with a row spacing of 42 cm and a planting density of 7 seedlings·m−2. The soil and cultivar in field experiments were representatives of those in the region. The irrigation, pesticide, and herbicide control practices followed local management for maize production.Figure 1The experimental design.Full size imageLAI measurements and UAV-based image acquisitionThe measurements of LAI were conducted at four growth stages including the tasseling stage (TS), flowering stage (FS), grain-filling stage (GS), and milk-ripe stage (MS) of maize in 2019 and 2020, a total of 264 LAI data of maize were collected during the 2-year field trial (Table 1). In order to reduce the impact of plant variability, the random sampling method was used to collect LAI samples. For each plot, three plants were randomly selected to measure the total green leaf area with the non-destructive portable leaf area meter (Laser Area Meter CI-203; CID Inc.). And the average leaf area of selected plants represented the single plant leaf area in each experiment plot. The LAI of each plot wasTable 1 Description of samplings.Full size table$$mathrm{LAI}=mathrm{LA}*mathrm{D}$$
(1)
where (mathrm{LA}) is the leaf area of a single plant in each plot; (mathrm{D}) is the planting density in one square meter.PHANTOM 4 PRO (DJI-Innovations Inc., Shenzhen, China) is a multi-rotor UAV equipped with a 20-megapixel visible-light camera that was employed to capture digital images. Aerial observations were conducted on the same dates as the LAI measurements, which was between 10:30 a.m. and 2:00 p.m. local time when the solar zenith angle was minimal. The UAV was flown automatically based on preset flight parameters and waypoints, with a forward overlap of 80% and a side overlap of 60%. A three-axis gimbal integrated with the inertial navigation system stabilized the camera, the automatic camera mode with fixed ISO (100) and a fixed exposure was used during the flight. Altogether, 4192 images were taken in eight flights from a flight height of 29.36 m above ground, with a spatial resolution of 0.008 m.The measurements of maize LAI were carried out with permission from the Modern Agricultural Research and Development Base of Henan Province. All experiments were carried out in accordance with relevant institutional, national, and international guidelines and legislation.Image pre-processingDJI Terra (version 2.3.3) was used to generate ortho-rectified images based on the structure from motion algorithms and a mosaic blending model. The main procedures are as follows: (1) extract feature points and match features according to the longitude, latitude, elevation, roll angle, pitch angle, and heading angle of each image; (2) build dense 3D point clouds by using dense multi-view stereo matching algorithm; (3) build a 3D polygonal mesh based on the vector relationship between each point in the dense cloud; (4) establish a 3D model with both external image and internal structure by merging the mosaic image into the 3D model; (5) generate digital orthophoto map (DOM).Vegetation indices (VIs) derived from the UAV-based digital imageryDigital imagery records the intensity of visible red (R), green (G), and blue (B) bands in individual pixels24. In order to enhance the vegetation parameters contained in the digital image, fourteen commonly used RGB-based VIs were collected, and their correlation with the LAI of maize at different growth stages was evaluated. Table 2 shows the detailed information of the selected RGB-based VIs.Table 2 RGB-based VIs for LAI estimation.Full size tableCentered on the point where LAI was measured, regions of interests (ROIs) with a size of 100*100 were clipped from the digital image. Python 3.7.3 was used for extracting the R, G, B information of maize and computing the RGB-based VIs from ROIs. In order to reduce the effects of light and shadow, the R, G, B color space of the image was normalized according to the followings:$$mathrm{r}=frac{R}{R+G+B}$$
(2)
$$g=frac{G}{R+G+B}$$
(3)
$$b=frac{B}{R+G+B}$$
(4)
where r, g, and b are the normalized values. R, G, B are the pixel values from the digital images based on each band.Pearson correlation analysisBefore regression analysis, the Pearson correlation analysis was performed to determine the relationship between maize LAI and different RGB-based VIs extracted from the digital image. Pearson correlation coefficient ((mathrm{r})) reflects the degree of linear correlation between two variables, which is between − 1 and 1. The calculation formula of Pearson correlation coefficient was expressed as follows:$$mathrm{r}= frac{sum_{i=1}^{n}left({X}_{i}-overline{X }right)left({Y}_{i}-overline{Y }right)}{sqrt{sum_{i=1}^{n}{left({X}_{i}-overline{X }right)}^{2}}sqrt{sum_{i=1}^{n}{left({Y}_{i}-overline{Y }right)}^{2}}}$$
(5)
where (X), (mathrm{Y}) are variables, (n) is the number of variables.Regression methodsLinear regression (LR)Linear regression is an approach for modelling the relationship between dependent and independent variables. The case of one independent variable is called unary linear regression (ULR), the expressions can be expressed as follows:$$mathrm{y}={beta }_{0}+{beta }_{1}x+varepsilon $$
(6)
where (varepsilon ) is deviation, which satisfies the normal distribution. (x), (mathrm{y}) are variables. ({beta }_{0}), ({beta }_{1}) are the intercept and slope of the regression line, respectively.For more than one independent variable, the regression process is called multiple linear regression (MLR), the expressions can be expressed as:$$mathrm{y}={beta }_{0}+{beta }_{1}{x}_{1}+{beta }_{2}{x}_{2}+dots +{beta }_{n}{x}_{n}$$
(7)
where ({x}_{1}),( {x}_{2}), …, ({x}_{n}), (mathrm{y}) are variables, ({beta }_{0}), ({beta }_{1}), ({beta }_{2}), …, ({beta }_{n}) are coefficients that determined by least square method and gradient descent method38.The RGB-based VIs with the highest Pearson correlation coefficient was used to establish the ULR model, and VIs with a correlation coefficient higher than 0.7 were used to establish the MLR model. In each growth stage, 70% of observation data were randomly selected for establishing models, and the remaining 30% of data were used as the testing dataset to assess the model performance.Back propagation neural networks (BPNN)In this study, a three-layer BPNN model was established for LAI estimation (Fig. 2). RGB-based VIs with a correlation coefficient higher than 0.7 were selected as the input variables. Tan-Sigmoid activation function was used in the hidden layer, and the Levenberg–Marquardt algorithm was selected as the training function. The maximum epoch of BPNN training was set to 1000, the learning rate was set to 0.005, and the MSE was set to 0.001. The observation data set was split into the training set and the testing dataset with a ratio of 7:3. The training dataset was used to fit the weights and bias of the BPNN model, the testing dataset was used to evaluate the model performance. Before training, data normalization was conducted for the input and output variables, and the denormalization was required to convent the output variable back into the original units after training.Figure 2Three-layer BPNN model.Full size imageRandom forest (RF)RF is a non-parametric ensemble ML method that operates by constructing a multitude of decision trees at training time and outputting the average prediction of the individual trees (Fig. 3). The bootstrapping approach was used to collect different sub-training data from the input training dataset to construct individual decision trees.Figure 3Random forest model.Full size imageThe construction process of RF regression model is as follows:

(1)

The value of (mathrm{n}_mathrm{estimators}) was tested from 50 to 1000 in increments of 50, and the value of 500 was finally selected according to higher R2 and lower RMSE.

(2)

At each node per tree, (mathrm{m}_mathrm{try}) RGB-based VIs was randomly selected from all 14 vegetation indices, and the best split was chosen according the lowest Gini Index. (mathrm{m}_mathrm{try}) was tested from 3 to 10, and the final value was 6.

(3)

The other parameters in the RF model were kept as default values according to the (mathrm{RandomForestRegressor}) function in (mathrm{Scikit}-mathrm{learn library}).

(4)

For each tree, the data splitting process in each internal node was repeated from the root node until a pre-defined stop condition was reached.

(5)

Similar with LR and BPNN model, the RGB-based VIs with a correlation coefficient higher than 0.7 were selected as the input variables, and the output variable is LAI.

Data analysis and performance evaluationThe repeated random sampling validation method was used to evaluate the generalization performance of different models. The training and testing dataset were randomly split 500 times. For each split, the LR, BPNN, and RF models were fitted to the training dataset, and the estimation accuracy was evaluated using the testing dataset. The coefficient of determination (R2), root mean square error (RMSE), and Akaike information criterion (AIC) of the training dataset were used for the assessment of models39, and the estimation accuracy was evaluated by R2 and RMSE of the testing dataset. Mathematically, a higher R2 corresponds to a smaller RMSE, and thus represents better model performance. The procedures of LAI inversion using UAV-based digital imagery and ML methods were shown in Fig. 4.Figure 4Flowchart of LAI inversion using UAV-based remote sensing and ML methods.Full size imageThe construction and evaluation of models was performed using Python 3.7.3 in Windows 10 operating system with Intel Core i7-9700 processor, 3.00 GHz CPU, and 32 GB RAM. The processing software is Spyder. The statistical analysis and figure plotting were performed in R × 64 4.0.3. More

Ethical statementThis research project was ethically reviewed and approved by the Ethics Committee of the Instituto Gulbenkian de Ciência (license reference: A009.2018), and by the Portuguese National Entity that regulates the use of laboratory animals (DGAV – Direção Geral de Alimentação e Veterinária (license reference: 008958). All experiments conducted on animals followed the Portuguese (Decreto-Lei n° 113/2013) and European (Directive 2010/63/EU) legislations, concerning housing, husbandry and animal welfare.Escherichia coli clonesThe ancestral invader E. coli strain expresses a Yellow Fluorescent Protein (YFP), and carries streptomycin and ampicillin resistance markers for easiness of isolation from the mouse feces [galK::amp (pZ12)::PLlacO−1-YFP, strR (rpsl150), ΔlacIZYA::scar]. An E. coli strain used for the in vivo competition experiments is isogenic to the ancestral invader but expresses a Cyan Fluorescent Protein (CFP) and carries streptomycin and chloramphenicol resistance markers [galK::chlor (pZ12)::PLlacO−1-CFP, strR (rpsl150), ΔlacIZYA::scar]. The resident E. coli lineage was isolated from the feces along time using McConkey + 0.4% lactose medium, as previously described9. All the resident clones sampled from each mouse belong to E.coli phylogenetic group B9. The invader E. coli strains (YFP and CFP) derive from the K-12 MG1655 strain (DM08) and exhibit a gat negative phenotype, gatZ::IS112. The resident E. coli clone used for the competition experiments in the mouse gut expresses a mCherry fluorescent protein and a chloramphenicol resistance marker, allowing to distinguish the invader and resident strains in the mice feces.E. coli clones were grown at 37 °C under aeration in liquid media Luria broth (LB) from SIGMA — or McConkey and LB agar plates. Media were supplemented with antibiotics streptomycin (100 µg/mL), ampicillin (100 µg/mL) or chloramphenicol (30 µg/mL) when specified.Serial plating of 1X PBS dilutions of feces in LB agar plates supplemented with the appropriate antibiotics were incubated overnight and YFP, CFP or mCherry-labeled bacterial numbers were assessed by counting the fluorescent colonies using a fluorescent stereoscope (SteREO Lumar, Carl Zeiss). The detection limit for bacterial plating was ~300 CFU/g of feces9.In vivo evolution and competition experimentsAll mice (Mus musculus) used in this study were supplied by the Rodent Facility at Instituto Gulbenkian de Ciência (IGC) and were given ad libitum access to food (Rat and Mouse No.3 Breeding (Special Diets Services) and water. Mice were kept at 20-24 °C and 40-60% humidity with a 12-h light-dark cycle. For the in vivo evolution experiment we used the gut colonization model previously established9. Briefly, mice drank water with streptomycin (5 g/L) only for 24 h before a 4 h starvation period of food and water. The animals were then inoculated by gavage with 100 µL of an E. coli bacterial suspension of ~108 colony-forming units (CFUs). Mice A2, B2, D2, E2, G2, H2 and I2 were successfully colonized with the invader E. coli, while mice C2 and F2 failed to be colonized. Six- to eight-week-old C57BL/6 J non-littermate female mice were kept in individually ventilated cages under specified pathogen-free (SPF) barrier conditions at the IGC animal facility. Fecal pellets were collected during more than one year ( >400 days) and stored in 15% glycerol at −80 °C for later analysis. In the competition experiments between the invader ancestral E. coli and evolved populations, we colonized the mice using a 1:1 ratio of each genotype, with bacterial loads being assessed and frozen on a daily basis after gavage.In vivo competition experiments in which the two modes of selection (directional and diversifying) were acting for a longer time period were performed using evolved invader E. coli populations colonizing mice D2, B2 and A2, H2. Here we used both male (n = 8) and female (n = 8) C57BL/6 J mice aged six- to eight-week-old treated with streptomycin during 3 days before gavage. E. coli populations evolving for short time periods do not allow for strong conclusions on which mode of selection is taking place. Evolved invader populations such as I2 or G2 were therefore not used for in vivo fitness assays. To assess the impact of the mouse resident E. coli in the competitive fitness of dgoR we performed one-to-one competitions between the invader ancestral and dgoR KO clones. We first homogenized the mice microbiotas by co-housing the animals during seven days. The animals (n = 6, female C57BL/6 J mice aged six- to eight-week-old) were then maintained under co-housing and given streptomycin-supplemented (5 g/L) water during seven days to break colonization resistance and eradicate their resident E. coli. At this point, the co-housed mice were removed from the antibiotic-supplemented water for two days. The following day, one group of mice was gavaged with an mCherry-expressing resident E. coli (n = 3 mice) while the other group (n = 3) was not, with all animals being individually caged from this point on and receiving normal water without antibiotic. The day after gavage, all mice were colonized with a mix (1:1) of the invader ancestral and the dgoR KO clones, and the bacterial loads were assessed and frozen on a daily basis.Microbiota analysisFecal DNA was extracted with a QIAamp DNA Stool MiniKit (Qiagen), according to the manufacturer’s instructions and with an additional step of mechanical disruption32. 16 S rRNA gene amplification and sequencing was carried out at the Gene Expression Unit from Instituto Gulbenkian de Ciência, following the service protocol. For each sample, the V4 region of the 16 S rRNA gene was amplified in triplicate, using the primer pair F515/R806, under the following PCR cycling conditions: 94 °C for 3 min, 35 cycles of 94 °C for 60 s, 50 °C for 60 s, and 72 °C for 105 s, with an extension step of 72 °C for 10 min. Samples were then pair-end sequenced on an Illumina MiSeq Benchtop Sequencer, following Illumina recommendations. Sampling for microbiota analysis was performed until the microbiota composition stabilized (~1 year after the antibiotic perturbation).QIIME2 version 2017.1133 was used to analyze the 16 S rRNA sequences by following the authors’ online tutorials (https://docs.qiime2.org/2017.11/tutorials/). Briefly, the demultiplexed sequences were filtered using the “denoise-single” command of DADA2 version 1.1434, and forward and reverse sequences were trimmed in the position in which the 25th percentile’s quality score got below 20. Diversity analysis was performed following the QIIME2 tutorial35. Beta diversity distances were calculated through Unweighted Unifrac36. For taxonomic analysis, OTU were picked by assigning operational taxonomic units at 97% similarity against the Greengenes database version 13 (Greengenes 13_8 99% OTUs (250 bp, V4 region 515 F/806 R))37.Whole-genome sequencing and analysis pipelineDNA was extracted38 from E. coli populations (mixture of  > 1000 clones) or a single clone growing in LB plates supplemented with antibiotic to avoid contamination. DNA concentration and purity were quantified using Qubit and NanoDrop, respectively. The DNA library construction and sequencing were carried out by the IGC genomics facility using the Illumina Miseq platform. Processing of raw reads and variants analysis was based on the previous work39. Briefly, sequencing adapters were removed using fastp version 0.20.040 and raw reads were trimmed bidirectionally by 4 bp window sizes across which an average base quality of 20 was required to be retained. Further retention of reads required a minimum length of 100 bps per read containing at least 50% base pairs with phred scores at or above 20. BBsplit (part of BBMap version 38.9)41 was used to remove likely contaminating reads as explained previously39. Separate reference genomes were used for the alignment of invader (K-12 (substrain MG1655; Accession Number: NC_000913.2)) and resident (Accession Number: SAMN15163749) E. coli genomes. Alignments were performed via three alignment approaches: BWA-sampe version 0.7.1742, MOSAIK version 2.743, and Breseq version 0.35.144,45. Final average alignment depths for invader and resident populations across time points equalled 302 (median = 236) and 253 (median = 235), respectively. While Breseq provides variant analysis in addition to alignment, other variant calling approaches were used to identify putative variation in the sequenced genomes, and to verify data from Breseq. A naïve pipeline39 using the mpileup utility within SAMtools version 1.946 and a custom script written in python was employed. Only reads with a minimum mapping quality of 20 were considered for analysis, and variant calling was limited to bases with call qualities of at least 30. At these positions, a minimum of 5 quality reads had to support a putative variant on both strands (with strand bias, pos. strand / neg. strand, above 0.2 or below 5) for further consideration. Finally, mutations were retained if detected in more than one of the alignment approaches, and if they reached a minimum frequency of 5% at a minimum of one time point sampled. Further simple and complex small variants were considered from freebayes version 0.9.2147 with similar thresholds, while insertion sequence movements and other mobile element activity was inferred via is mapper version 248 and panISa version 0.1.649, as well as Breseq, as previously described39. All putative variants were verified manually in IGV version 2.750,51. Raw sequencing reads were deposited in the sequence read archive under bioproject PRJNA666769. Population dynamics of lineage-specific dynamics and the resulting Muller plots were inferred manually and are meant strictly as a means of presenting the data. In order to generate these plots, mutations were sorted by frequency (descending for each time point at which the population was sampled). The largest frequency mutations were considered major lineages within which minor frequency mutations occurred. Assuming that a mutation, which arises subsequent to a preexisting mutation (an already differentiated lineage) cannot exceed the frequency of that preexisting mutation at any point, and will fluctuate in frequency with the preexisting one, we assigned mutations to the lineages within each population. While this resolved the majority of high frequency and medium frequency mutations, low-frequency mutations within the Muller plots cannot be placed with high confidence, and are only included for completeness.Prophage induction rateTo calculate the maximum prophage induction rate we grew E. coli lysogenic clones, starting with the same initial OD600 values: ~0.1 (Bioscreen C system, Oy Growth Curves Ab Ltd), with agitation at 37 °C in LB medium in the presence or absence of mitomycin C along time (5 µg/mL)9. The OD600 values were normalized by dividing the ones in the presence of mitomycin C by the ones in the absence of mitomycin C (sampling interval: 30 min). The LN of this ratios along time originates a lysis curve, where the maximum slope corresponds to the maximal prophage induction rate for each clone analyzed. We tested evolved clones from mouse A2, H2 and G2 against the ancestral clone which only carries the Nef and the KingRac prophages. We also tested clones of the resident strain that had evolved in the presence of the invader for more than 400 days (these clones were sampled from mouse A2).
E. coli growth rate, growth curves, cell aggregation, biofilm and motility capacityTo calculate the maximum bacterial growth rate, we grew E. coli lysogenic clones, starting with the same initial OD600 values: ~0.1 (Bioscreen C system, Oy Growth Curves Ab Ltd), with agitation at 37 °C in LB medium along time using reading intervals of 30 min. The LN of the OD600 values along time originates a growth curve, where the maximum slope corresponds to the maximum bacterial growth rate for each clone analyzed.To test for metabolic differences of the psuK/fruA mutation, growth curves of evolved lysogenic E. coli clones, bearing the Nef and KingRac prophages, with or without the psuK/fruA mutation were performed with the same initial OD600 value (~0.03) for each clone. The clones were grown in glucose (0.4%) minimal medium (MM9-SIGMA) with or without pseudouridine (80 μM) and absorbance values were obtained using the Bioscreen C apparatus during 12 h.Frozen stocks of E. coli clones were used to seed tubes with 5 mL of liquid LB. These were incubated overnight at 37 °C under static conditions to assess the formation of cell flocks/clumps, observable to the naked eye, in order to evaluate the formation of cell aggregates. Biofilm was tested according a previously published protocol52 and to evaluate the motility capacity we adapted the protocol from Croze and colleagues53. Briefly, overnight E. coli clonal cultures grown with agitation at 37 °C in 5 mL LB medium supplemented with streptomycin (100 ug/mL) were adjusted to the same absorbance and a 3uL volume was dropped on top of soft agar (0.25%). Plates were incubated at 37 °C and photos were taken at day 1, 2 and 5 post-inoculation to assess swarming motility phenotype.Number of E. coli generations during mouse gut colonizationTo estimate the number of generations of E. coli in the mouse gut, we used a previously described protocol to measure the fluorescent intensity of a probe specific to E. coli 23 S rRNA (as a measure of ribosomal content) that correlates with the growth rate of the bacterial cells54. We measured the number of generations of the ancestral E. coli clone while colonizing the gut of 2 mice, treated during 24 h with streptomycin (5 g/L) before gavage, during 25 days.Plasmid DNA extraction and PCR detection of ~69Kb (repA) and ~109Kb (repB) plasmidsPlasmid DNA was extracted from overnight cultures using a Plasmid Mini Kit (Qiagen), according to the manufacturer’s guidelines. Specific primers for the amplification of repA and repB genes, were used to determine the frequency of the 68935 bp (~69 Kb) and 108557 bp (~109 Kb) plasmids, respectively, in the invader E. coli population.The primers used for repA gene were:repA-Forward: 5’-CAGTCCCCTAAAGAATCGCCCC-3’ and repA-Reverse: 5’-TGACCAGGAGCGGCACAATCGC-3’.For repB the primer sequences were:repB-Forward: 5’-GTGGATAAGTCGTCCGGTGAGC-3’ and repB-Reverse: 5’-GTTCAAACAGGCGGGGATCGGC3’.PCR amplification of plasmid-specific genes was performed in 12 isolated random clones from mouse A2 at days 104 and 493. PCR reactions were performed in a total volume of 25 μL, containing 1 μL of plasmid DNA, 1X Taq polymerase buffer, 200 μM dNTPs, 0.2 μM of each primer and 1.25 U Taq polymerase. PCR reaction conditions: 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 65 °C for 30 s and 72 °C for 30 s, finalizing with 5 min at 72 °C. DNA was visualized on a 2% agarose gel stained with GelRed and run at 160 V for 60 min.Construction of the dgoR KO mutantP1 transduction was used to construct a ΔdgoR mutant (dgoR KO). This KO strain was created by replacing the wild-type dgoR in the invader ancestral YFP-expressing genetic background by the respective knock-out from the KEIO collection, strain JW562755, in which the dgoR sequence is replaced by a kanamycin resistance cassette. The presence of the cassette was confirmed by PCR using primers dgoK-F: GCGATGTAGCGAGCTGTC, and yidX-R: GGGAATAAACCGGCAGCC. PCR reactions were performed in a total volume of 25 μL, containing 1 μL of DNA, 1X Taq polymerase buffer, 200 μM dNTPs, 0.2 μM of each primer and 1.25 U Taq polymerase. PCR reaction conditions: 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 65 °C for 30 s and 72 °C for 30 s, finalizing with 5 min at 72 °C. DNA was visualized in a 2% agarose gel stained with GelRed and run at 160 V for 60 min.RNA extraction, DNAse treatment, RT-PCR and qPCRThe Qiagen RNeasy Mini Kit was used for RNA extraction. RNA concentration and quality were evaluated in the Nanodrop 2000 and by gel-electrophoresis. DNase treatment was performed with the RQ1 DNase (Promega) by adding 0.5 μl of DNase to 1 μg of RNA and 1 μl buffer in a final volume of 15 ul, followed by incubation 30 min at 37 °C. Afterwards, 1 ul of stop solution was added and incubation for 15 min at 65 °C was performed to inactivate the DNase. As a control for complete DNA digest a PCR was performed on the reactions including positive controls. Reverse transcription was performed with M-MLV RT[-H] (Promega) by mixing 1 μg of RNA with 0.5 μl random primers (Promega) and nuclease free water to a volume of 15 μl, incubation at 70 °C for 5 min and a quick cool down on ice. Afterwards the reverse transcription was accomplished by adding 5 μl of RT buffer, 0.5 μl RT enzyme and 2 μl dNTP mix, followed by incubation for 10 min at 25 °C, 50 min at 50 °C and 10 min at 70 °C. The resulting cDNA was diluted 100-fold in nuclease free water before changes in gene expression were detected using the The QuantStudio 7Flex (Applied Biosystems) with iTaq Universal SYBR Green Supermix (BioRad) and the following cycling protocol: Hold stage: 2 min at 50 °C, 10 min at 95 °C. PCR stage (40 cycles): 15 s at 95 °C, 30 s at 58 °C, 30 s at 60 °C. Melt curve stage: 15 s at 95 °C, 1 min at 50 °C then increments of 0.05 °C/s until 95 °C. Melt curve analysis was performed to verify product homogeneity. All reactions included six biological and three technical replicates for each sample. A relative quantification method of analysis with normalization against the endogenous control rrsA and employing the primer specific efficiencies was used according to the Pfaffl method (add reference). The primers used were designed with PrimerQuest (idt). The used primer sequences were: psuK – TGCGTTAGCAGCGATTGA, AATTTACGCCTGGTGGAGTAG; arcA – GATTCATGGTACGGGACAGTAG, CCGTGACAACGAAGTCGATAA; yjtD – CGCACATGGATCTGGTGATA, GGCGTGGCGTAGTAATGATA and rrsR – GTCAGCTCGTGTTGTGAAATG, CCCACCTTCCTCCAGTTTATC.Statistics and reproducibilityCorrelation between microbiota diversity measures and E. coli loads (CFU) or persistence (1-presence or 0-absence) was performed in R using the statistical package rmcorr (version 0.5.2)56 and lme4 (version 1.1-10)57, respectively. The rate of accumulation of new ISs in vivo was compared using Wilcoxon paired signed ranked test for expected and observed insertions, while the rate of selective sweeps correlation was performed using the Spearman Correlation test. Selective sweeps were taken to be mutations or HGT events that reached  > 95% frequency in the population and kept high frequency until the end of the colonization. Statistical analysis of prophage induction as well as biofilm levels was performed using the Mann-Whitney test in GraphPad Prism (version 8.4.3). A single sample T-Test was used test if the growth rate of evolved invader clones deviates from the mean of the ancestral. A Wilcoxon rank sum test with continuity correction was used to compare the relative expression levels of the evolved clones with the ancestral. P values of x0 are of order of their inverse selection coefficient (up to logarithmic corrections):$${{{{{rm{G}}}}}}left(xright)approx 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{s}.$$
(1)

Clonal interference under uniform directional selection. This mode occurs in asexual populations when adaptive mutations become frequent enough to interfere with one another59,60,61. Only a fraction of the established adaptive mutations reaches fixation; sojourn times to intermediate frequencies are set by a global coalescence rate (widetilde{sigma }) that is higher than the typical selection coefficient of individual mutations62:

$${{{{{rm{G}}}}}}left(xright) < 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{widetilde{sigma}}.$$ (2) Details of these dynamics depend on the spectrum of selection coefficients and on the overall mutation rate, which set the strength of clonal interference. For moderate interference, where a few concurrent beneficial mutations compete for fixation, we expect a roughly exponential drop of the frequency propagator, (Gleft(xright)sim {{exp }}left(-lambda xright)), reflecting the probability that a trajectory reaches frequency x without interference by a stronger competing clade. Moderate interference generates an effective neutrality for weaker beneficial mutations and at higher frequencies63. This regime has been mapped for influenza64. In the asymptotic regime of a travelling fitness wave, where many beneficial mutations are simultaneously present, the fate of a mutation is settled in the range of small frequencies; that is, at the tip of the wave65. In this regime, emergent neutrality affects the vast majority of beneficial mutations and most of the frequency regime66. Hence, the frequency propagator rapidly drops to its asymptotic value (Gleft(x=1right)ll 1.) Adaptation under diversifying selection. More complex selection scenarios involve selection within and between ecotypes, i.e., subpopulations occupying distinct ecological niches67,68. An important factor generating niches and ecotypes is the differential use of food and other environmental resources. In this mode, ecotype-specific, conditionally beneficial mutations reach intermediate frequencies after a time given by their within-ecotype selection coefficients, but fixation can be slowed down or suppressed by diversifying (negative frequency-dependent) cross-ecotype selection18, $${{{{{rm{G}}}}}}left(xright)approx 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{s},left(xlesssim ,frac{1}{2}right)$$ (3) $${{{{{rm{G}}}}}}left(xright) < 1,{{{{{rm{T}}}}}}left(xright)gg frac{1}{s}left(xto 1right).$$ (4) The details depend on the details of the eco-evolutionary model (synergistic vs. antagonistic interactions, carrying capacities, amount of resource competition vs. explicitly frequency-dependent selection). In a model with directional selection within ecotypes, conditionally beneficial mutations rapidly fix within ecotypes, but lead only to finite shifts of the ecotype frequencies. In the simplest case, the resulting dynamics of ecotype frequencies is diffusive, resulting in an effectively neutral turnover of ecotypes18. Given negative frequency-dependent selection between ecotypes, fixations become even rarer and can be completely suppressed; that is, ecotypes can become stable on the time scales of observation. The separation of time and selection scales between intra- and cross-ecotype frequency changes is expected to be a robust feature of ecotype-dependent selection: sojourn of adaptive alleles to intermediate frequencies is fast, fixation is slower and rarer. In other words, ecotype-dependent selection is characterized by two regimes of coalescence times T(x).Frequency propagators and the coalescence time spectra expected under these evolutionary modes are qualitatively sketched in Supplementary Fig. 11. For periodic sweeps under directional selection (dark green, left column), G(x) depends weakly on x and T(x) is set by rapid sweeps for all x. For clonal interference under directional selection (green, center column), G(x) decreases substantially with increasing x and T(x) becomes uniformly shorter. Under negative frequency-dependent selection (brown, right column), G(x) decreases substantially with increasing x, while T(x) substantially increases for large x and diverges in case of strong frequency-dependent selection generating stable ecotypes (dashed lines). (see Supplementary Fig. 11 for the results of simulations assuming a model of direction selection or assuming a resource competition model where ecotype formation occurs31.The ({{{{{boldsymbol{p}}}}}})-({{{{{boldsymbol{tau }}}}}}) selection testThis test is based on qualitative characteristics of the functions G(x), T(x) and does not depend on details of the evolutionary process. We evaluate G(x) and T(x) for host-specific families of frequency trajectories; sojourn times are counted from an initial frequency x0=0.01. Origination times at this frequency are inferred by backward extrapolation of the first observed trajectory segment; the reported results are robust under variations of the threshold x0 and the extrapolation procedure. We then compute two summary statistics: the probability (p) that a mutation established at an intermediate frequency xm reaches near-fixation at a frequency xf,$$p=frac{{{{{{rm{G}}}}}}({x}_{f})}{{{{{{rm{G}}}}}}({x}_{m})},$$ (5) and the corresponding fraction of sojourn times,$$tau=,frac{{{{{{rm{T}}}}}}({x}_{f})}{{{{{{rm{T}}}}}}({x}_{m})}.$$ (6) Here we use xm=0.3 and xf=0.95 to limit the uncertainties of empirical trajectories at low and high frequency; however, the selection test is robust under variation of these frequencies. We find evidence for different modes of evolution: The long-term frequency trajectories of mice B2, D2 and E2 are consistent with predominantly frequency-dependent selection (Fig. 2, Fig. 4a–c). The propagator G(x) is a strongly decreasing function of x, resulting in fixation probabilities (p) 0.6, as measured by time ratios τ  > 3.

The trajectories of mice A2, G2, and I2 show a signature of recurrent selective sweeps and clonal interference under uniform directional selection (Fig. 4a–c). The propagator G(x) is a decreasing function of x, resulting in fixation probabilities (p=0.2-0.8), depending on the strength of clonal interference. Fixation times are short, giving time ratios (tau lesssim 2).

The shorter trajectory of mouse H2 signals periodic sweeps under uniform directional selection (Fig. 3, Fig. 4a–c). The origination rate of mutations is lower than in the longer trajectories, and G(x) shows a weak decrease with (p=1.) Sojourn times T(x) are short and grow uniformly with x, resulting in a time ratio τ=2.25. This pattern is expected under directional selection in the low mutation regime: (Tleft(xright)={{log }}left[x/(1-x)right]/{s}) for individual mutations with a uniform selection coefficient s, leading to τ=2.0 for xm=0.3 and xf=0.95 (this value is marked as a dashed line in Fig. 4c).

The trajectories of non-mutator lines in the long-term in vitro evolution experiment of Good et al1, evaluated over the first 7500 generations, show an overall signal of clonal interference under uniform directional selection (Fig. 4c, Supplementary Fig. 12). The frequency propagators G(x) are strongly decreasing functions of x and sojourn times T(x) grow uniformly with x. We find (p=0.2-0.8) and (tau lesssim 2), similar to the pattern in mice A2, G2, and I2.

Code for Selection testsThe code for selection tests from the mutation frequency trajectories can be found in the Supplementary Information file.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Woodward SE, Krekhno Z, Finlay BB. Here, There, and Everywhere: How Pathogenic Escherichia coli Sense and Respond to Gastrointestinal Biogeography. Cell Microbiol. 2019;21:e13107.CAS 
Article 

Google Scholar 
Collins JW, Keeney KM, Crepin VF, Rathinam VAK, Fitzgerald KA, Finlay BB, et al. Citrobacter rodentium: Infection, Inflammation and the Microbiota. Nat Rev Microbiol. 2014;12:612–23.CAS 
Article 

Google Scholar 
Mowat AM, Agace WW. Regional Specialization Within the Intestinal Immune System. Nat Rev Immunol. 2014;14:nri3738.Article 

Google Scholar 
Population Dynamics – Latest Research and News | Nature [Internet]. [cited 2017 Nov 5]. Available from: https://www.nature.com/subjects/population-dynamicsAbel S, Abel zur Wiesch P, Davis BM, Waldor MK. Analysis of Bottlenecks in Experimental Models of Infection. PLoS Pathog. 2015;11:e1004823.Global Burden of Disease Study. 2013 Collaborators. Global, Regional, and National Incidence, Prevalence, and Years Lived with Disability for 301 Acute and Chronic Diseases and Injuries in 188 Countries, 1990-2013: A Systematic Analysis for the Global Burden of Disease Study 2013. Lancet Lond Engl 2015;386:743–800.Article 

Google Scholar 
WHO | Diarrhoeal Disease [Internet]. WHO. [cited 2017 Oct 22]. Available from: http://www.who.int/mediacentre/factsheets/fs330/en/.Crepin VF, Collins JW, Habibzay M, Frankel G. Citrobacter rodentium Mouse Model of Bacterial Infection. Nat Protoc. 2016;11:1851–76.CAS 
Article 

Google Scholar 
Bhinder G, Sham HP, Chan JM, Morampudi V, Jacobson K, Vallance BA. The Citrobacter rodentium Mouse Model: Studying Pathogen and Host Contributions to Infectious Colitis. J Vis Exp. 2013;19:e50222.
Google Scholar 
Barthold SW, Coleman GL, Bhatt PN, Osbaldiston GW, Jonas AM. The Etiology of Transmissible Murine Colonic Hyperplasia. Lab Anim Sci. 1976;26:889–94.CAS 
PubMed 

Google Scholar 
Deng W, Vallance BA, Li Y, Puente JL, Finlay BB. Citrobacter rodentium Translocated Intimin Receptor (Tir) is an Essential Virulence Factor Needed for Actin Condensation, Intestinal Colonization and Colonic Hyperplasia in Mice. Mol Microbiol. 2003;48:95–115.CAS 
Article 

Google Scholar 
Schauer DB, Falkow S. Attaching and Effacing Locus of a Citrobacter freundii Biotype That Causes Transmissible Murine Colonic Hyperplasia. Infect Immun. 1993;61:2486–92.CAS 
Article 

Google Scholar 
Elliott SJ, Yu J, Kaper JB. The Cloned Locus of Enterocyte Effacement From Enterohemorrhagic Escherichia coli O157:H7 is Unable to Confer the Attaching and Effacing Phenotype Upon E. coli K-12. Infect Immun. 1999;67:4260–3.CAS 
Article 

Google Scholar 
Mellies JL, Elliott SJ, Sperandio V, Donnenberg MS, Kaper JB. The Per Regulon of Enteropathogenic Escherichia coli: Identification of a Regulatory Cascade and a Novel Transcriptional Activator, the Locus of Enterocyte Effacement (LEE)-Encoded Regulator (Ler). Mol Microbiol. 1999;33:296–306.CAS 
Article 

Google Scholar 
Jarvis KG, Girón JA, Jerse AE, McDaniel TK, Donnenberg MS, Kaper JB. Enteropathogenic Escherichia coli Contains a Putative Type III Secretion System Necessary for the Export of Proteins Involved in Attaching and Effacing Lesion Formation. Proc Natl Acad Sci. 1995;92:7996–8000.CAS 
Article 

Google Scholar 
Moon HW, Whipp SC, Argenzio RA, Levine MM, Giannella RA. Attaching and Effacing Activities of Rabbit and Human Enteropathogenic Escherichia coli in Pig and Rabbit Intestines. Infect Immun. 1983;41:1340–51.CAS 
Article 

Google Scholar 
Abel S, Abel zur Wiesch P, Chang HH, Davis BM, Lipsitch M, Waldor MK. Sequence Tag–based Analysis of Microbial Population Dynamics. Nat Methods. 2015;12:223–6.CAS 
Article 

Google Scholar 
Cavalli-Sforza LL, Edwards AW. Phylogenetic analysis. Models and Estimation Procedures. Am J Hum Genet. 1967;19:233–57. 3.CAS 
PubMed 
PubMed Central 

Google Scholar 
Zhang T, Abel S, Abel Zur Wiesch P, Sasabe J, Davis BM, Higgins DE, et al. Deciphering the Landscape of Host Barriers to Listeria monocytogenes Infection. Proc Natl Acad Sci . 2017;114:6334–9.CAS 
Article 

Google Scholar 
Hullahalli K, Waldor MK. Pathogen Clonal Expansion Underlies Multiorgan Dissemination and Organ-Specific Outcomes During Murine Systemic Infection. eLife 2021;10:e70910.CAS 
Article 

Google Scholar 
Petty NK, Feltwell T, Pickard D, Clare S, Toribio AL, Fookes M, et al. Citrobacter rodentium is an Unstable Pathogen Showing Evidence of Significant Genomic Flux. PLOS Pathog. 2011;7:e1002018.CAS 
Article 

Google Scholar 
Yasutomi E, Hoshi N, Adachi S, Otsuka T, Kong L, Ku Y, et al. Proton Pump Inhibitors Increase the Susceptibility of Mice to Oral Infection with Enteropathogenic Bacteria. Dig Dis Sci. 2018;63:881–9.CAS 
Article 

Google Scholar 
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-Resolution Sample Inference From Illumina Amplicon Data. Nat Methods. 2016;13:581–3.CAS 
Article 

Google Scholar 
Krimbas CB, Tsakas S. The Genetics of Dacus Oleae. v. Changes of Esterase Polymorphism in a Natural Population Following Insecticide Control-Selection or Drift? Evol Int J Org Evol. 1971;25:454–60.Article 

Google Scholar 
Peter BM, Slatkin M. Detecting Range Expansions From Genetic Data. Evol Int J Org Evol. 2013;67:3274–89.Article 

Google Scholar 
Schwarz R, Kaspar A, Seelig J, Künnecke B. Gastrointestinal Transit Times in Mice and Humans Measured with 27Al and 19F Nuclear Magnetic Resonance. Magn Reson Med. 2002;48:255–61.Article 

Google Scholar 
McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PloS One. 2013;8:e61217.CAS 
Article 

Google Scholar 
Petty NK, Bulgin R, Crepin VF, Cerdeño-Tárraga AM, Schroeder GN, Quail MA, et al. The Citrobacter rodentium Genome Sequence Reveals Convergent Evolution with Human Pathogenic Escherichia coli. J Bacteriol. 2010;192:525–38.CAS 
Article 

Google Scholar 
Mundy R, Pickard D, Wilson RK, Simmons CP, Dougan G, Frankel G. Identification of a Novel Type IV Pilus Gene Cluster Required for Gastrointestinal Colonization of Citrobacter rodentium. Mol Microbiol. 2003;48:795–809.CAS 
Article 

Google Scholar 
Darwin AJ, Miller VL. Identification of Yersinia enterocolitica Genes Affecting Survival in an Animal Host Using Signature-Tagged Transposon Mutagenesis. Mol Microbiol. 1999;32:51–62.CAS 
Article 

Google Scholar 
Maroncle N, Balestrino D, Rich C, Forestier C. Identification of Klebsiella pneumoniae Genes Involved in Intestinal Colonization and Adhesion Using Signature-Tagged Mutagenesis. Infect Immun. 2002;70:4729–34.CAS 
Article 

Google Scholar 
Wiles S, Dougan G, Frankel G. Emergence of a ‘Hyperinfectious’ Bacterial State After Passage of Citrobacter rodentium Through the Host Gastrointestinal Tract. Cell Microbiol. 2005;7:1163–72.CAS 
Article 

Google Scholar 
Wiles S, Clare S, Harker J, Huett A, Young D, Dougan G, et al. Organ Specificity, Colonization and Clearance Dynamics In Vivo Following Oral Challenges with the Murine Pathogen Citrobacter rodentium. Cell Microbiol. 2004;6:963–72.CAS 
Article 

Google Scholar 
Kitamoto S, Nagao-Kitamoto H, Kuffa P, Kamada N. Regulation of Virulence: The Rise and Fall of Gastrointestinal Pathogens. J Gastroenterol. 2016;51:195–205.Article 

Google Scholar 
Takumi K, de Jonge R, Havelaar A. Modelling Inactivation of Escherichia coli by Low pH: Application to Passage Through the Stomach of Young and Elderly People. J Appl Microbiol. 2000;89:935–43.CAS 
Article 

Google Scholar 
Pienaar JA, Singh A, Barnard TG. Acid-Happy: Survival and Recovery of Enteropathogenic Escherichia coli (EPEC) in Simulated Gastric Fluid. Micro Pathog. 2019;128:396–404.CAS 
Article 

Google Scholar 
McConnell EL, Basit AW, Murdan S. Measurements of Rat and Mouse Gastrointestinal pH, Fluid and Lymphoid Tissue, and Implications for In-Vivo Experiments. J Pharm Pharm. 2008;60:63–70.CAS 
Article 

Google Scholar 
Tan A, Petty NK, Hocking D, Bennett-Wood V, Wakefield M, Praszkier J, et al. Evolutionary Adaptation of an AraC-Like Regulatory Protein in Citrobacter rodentium and Escherichia Species. Infect Immun. 2015;83:1384–95.CAS 
Article 

Google Scholar 
Sanchez KK, Chen GY, Schieber AMP, Redford SE, Shokhirev MN, Leblanc M, et al. Cooperative Metabolic Adaptations in the Host Can Favor Asymptomatic Infection and Select for Attenuated Virulence in an Enteric Pathogen. Cell 2018;175:146–.e15.CAS 
Article 

Google Scholar 
Cunningham R, Jones L, Enki DG, Tischhauser R. Proton Pump Inhibitor Use as a Risk Factor for Enterobacteriaceal Infection: A Case-Control Study. J Hosp Infect. 2018;100:60–4.CAS 
Article 

Google Scholar 
Kelly OB, Dillane C, Patchett SE, Harewood GC, Murray FE. The Inappropriate Prescription of Oral Proton Pump Inhibitors in the Hospital Setting: A Prospective Cross-Sectional Study. Dig Dis Sci. 2015;60:2280–6.CAS 
Article 

Google Scholar 
Federici S, Suez J, Elinav E. Our Microbiome: On the Challenges, Promises, and Hype. Results Probl Cell Differ. 2020;69:539–57.CAS 
Article 

Google Scholar  More