Philippa Kaur

Bellot, S. et al. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01858-0 (2022).Article 

Google Scholar 
Eiserhardt, W. L. et al. Ann. Bot. 108, 1391–1416 (2011).Article 

Google Scholar 
Muscarella, R. et al. Glob. Ecol. Biogeogr. 29, 1495–1514 (2020).Article 

Google Scholar 
Cámara-Leret, R. et al. Nat. Plants 3, 16220 (2017).Article 

Google Scholar 
The IUCN Red List of Threatened Species (IUCN, 2018).BGCI ThreatSearch Online Database (BGCI, 2018).Carlos-Júnior, L. A. et al. Divers. Distrib. 25, 743–757 (2019).Article 

Google Scholar 
Pollock, L. J., Thuiller, W. & Jetz, W. Nature 546, 141–144 (2017).CAS 
Article 

Google Scholar 
Rice, J. et al. The IPBES Regional Assessment Report on Biodiversity and Ecosystem Services for the Americas (IPBES, 2018).Coelho de Souza, F. et al. Nat. Ecol. Evol. 3, 1754–1761 (2019).Article 

Google Scholar  More

Nansen, C., Phillips, T. W., Parajuleeb, M. N. & Franqui, R. A. Comparison of direct and indirect sampling procedures for Plodia interpunctella in a maize storage facility. J. Stored Prod. Res. 40, 151–168 (2004).Article 

Google Scholar 
Gerken, A. R. & Campbell, J. F. Using long-term capture data to predict Trogoderma variabile Ballion and Plodia interpunctella (Hübner) population patterns. Insects 10, 93. https://doi.org/10.3390/insects10040093 (2019).Article 
PubMed Central 

Google Scholar 
Athanassiou, C. G. & Buchelos, C. T. Grain properties and insect distribution trends in silos of wheat. J. Stored Prod Res. 88, 101632 (2020).Article 

Google Scholar 
Campbell, J., Mullen, M. & Dowdy, A. Monitoring stored-product pests in food processing plants with pheromone trapping, contour mapping, and mark-recapture. J. Econ. Entomol. 95, 1089–1101 (2002).CAS 
PubMed 
Article 

Google Scholar 
Arbogast, R. T., Weaver, D. K., Kendra, P. E. & Brenner, R. J. Implications of spatial distribution of insect populations in storage ecosystems. Environ. Entomol. 27, 202–216 (1998).Article 

Google Scholar 
Brenner, R. J., Focks, D. A., Arbogast, R. T., Weaver, D. K. & Shuman, D. Practical use of spatial analysis in precision targeting for integrated pest management. Am. Entomol. 44, 79–102 (1998).Article 

Google Scholar 
Arbogast, R. T., Kendra, P. E., Mankin, R. W. & McGovern, J. E. Monitoring insect pests in retail stores by trapping and spatial analysis. J. Econ. Entomol. 93, 1531–1542 (2000).CAS 
PubMed 
Article 

Google Scholar 
Arthur, F. & Phillips, T.W. Stored-product insect pest management and control. In Food Plant Sanitation; Hui, Y.H., Bruinsma, B.L., Gorham, J.R., Nip, W.-K., Tong, P.S., Ventresca, P., Eds.; Marcel Dekker, Inc, pp. 341–348(2003).Campbell, J. F., Toews, M. D., Arthur, F. H. & Arbogast, R. T. Long-term monitoring of Tribolium castaneum in two flour mills: Seasonal patterns and impact of fumigation. J. Econ. Entomol. 103, 991–1001 (2010).PubMed 
Article 

Google Scholar 
Doud, C. W. & Phillips, T. W. Activity of Plodia interpunctella (Lepidoptera: Pyralidae) in and around flour mills. J. Econ. Entomol. 93, 1842–1847 (2000).CAS 
PubMed 
Article 

Google Scholar 
Campbell, J. & Mullen, M. Distribution and dispersal behavior of Trogoderma variabile and Plodia interpunctella outside a food processing plant. J. Econ. Entomol. 97, 1455–1464 (2004).CAS 
PubMed 
Article 

Google Scholar 
Larson, Z., Subramanyam, B. & Herrman, T. Stored-product insects associated with eight feed mills in the Midwestern United States. J. Econ. Entomol. 101, 998–1005 (2008).PubMed 
Article 

Google Scholar 
Trematerra, P., Paula, M. C., Sciarretta, A. & Lazzari, S. Spatio-temporal analysis of insect pests infesting a paddy rice storage facility. Neotrop. Entomol. 33, 469–479 (2004).Article 

Google Scholar 
Arthur, F. H., Campbell, J. F. & Toews, M. D. Distribution, abundance, and seasonal patterns of Plodia interpunctella (Hübner) in a commercial food storage facility. J. Stored Prod. Res. 53, 7–14 (2013).Article 

Google Scholar 
McKay, T., White, A. L., Starkus, L. A., Arthur, F. H. & Campbell, J. F. Seasonal patterns of stored-product insects at a rice mill. J. Econ. Entomol. 110, 1366–1376 (2017).PubMed 
Article 

Google Scholar 
Roesli, R., Subramanyam, B., Fairchild, F. J. & Behnke, K. C. Trap catches of stored-product insects before and after heat treatment in a pilot feed mill. J. Stored Prod. Res. 39, 521–540 (2003).Article 

Google Scholar 
Campbell, J., Ching’oma, G.M., Toews, M.D. & Ramaswamy, S. Spatial distribution and movement patterns of stored-product insects. In Proceedings of the 9th International Working Conference on Stored Product Protection, Campinas, Sao Paulo, Brazil, 15–18 October 2006; Lorini, I., Bacaltchuk, B., Beckel, H., Deckers, D., Sundfeld, E., Santos, J.P.D., Biagi, J.D., Celaro, J.C., Faroni, L.R.D., Bortolini, L.D.F., Eds.; Brazilian Post-harvest Association—ABRAPOS: Passo Fundo, RS, Brazil, p. 18 (2006).Trematerra, P., Gentile, P., Brunetti, A., Collins, L. & Chambers, J. Spatio-temporal analysis of trap catches of Tribolium confusum du Val in a semolina-mill, with a comparison of female and male distributions. J. Stored Prod. Res. 43, 315–322 (2007).Article 

Google Scholar 
Semeao, A. A., Campbell, J. F., Whitworth, R. J. & Sloderbeck, P. E. Influence of environmental and physical factors on capture of Tribolium castaneum (Coleoptera: Tenebrionidae) in a flour mill. J. Econ. Entomol. 105, 686–702 (2012).PubMed 
Article 

Google Scholar 
Campbell, J.F., Perez-Mendoza, J. &Weier, J. Insect Pest Management Decisions in Food Processing Facilities. In Stored Product Protection; Hagstrum, D.W., Phillips, T.W., Cuperus, G., Eds.; Kansas State University, pp. 219–232 (2012).Mohandass, S., Arthur, F. H., Zhu, K. & Throne, J. E. Biology and management of Plodia interpunctella (Lepidoptera:Pyralidae) in stored products. J. Stored Prod. Res. 43, 302–311 (2007).Article 

Google Scholar 
Hamlin, J.C., Reed, W.D. & Phillips, M.E. Biology of the Indianmeal Moth on Dried Fruits in California. USDA Technical Bulletin No. 242, (1931)Hagstrum, D.W. & Subramanyam, B. Review of Stored-Product Insect Resource. AACC International (2009).Soderstrom, T., Stoica, P. & Trulsson, E. Instrumental variable methods for closed loop systems. IFAC 10th Triennial World Congress, Munich, FRG. pp. 363–368(1987).Johnson, J. A., Valero, K. A., Hannel, M. M. & Gill, R. F. Seasonal occurrence of post harvest dried fruit insects and their parasitoids in a culled fig warehouse. J. Econ. Entomol. 93, 1380–1390 (2000).CAS 
PubMed 
Article 

Google Scholar 
Nansen, C., Subramanyam, B. & Roesli, R. Characterizing spatial distribution of trap captures of beetles in retail pet stores using SADIE® software. J. Stored Prod. Res. 40, 471–483 (2004).Article 

Google Scholar 
Phillips, T.W., Berbert, R.C. &Cuperus, G.W. Post-harvest integrated pest management. In: Francis, F.J. (Ed.), Encyclopedia of Food Science and Technology. 2nd ed. Wiley Inc., pp. 2690–2701(2000).Phillips,T.W., Cogan, P.M. & Fadamiro, H.Y. Pheromones. In: Subramanyam, B., Hagstrum, D.W. (Eds.), Alternatives to Pesticides in Stored-product IPM. Kluwer Academic Publishers, pp. 273–302 (2000).Mullen, M. A. & Dowdy, A. K. A pheromone-baited trap for monitoring the Indian meal moth, Plodia interpunctella (Hübner) (Lepidoptera: Pyralidae). J. Stored Prod. Res. 37, 231–235 (2001).PubMed 
Article 

Google Scholar 
Nansen, C. & Phillips, T. W. Ovipositional responses of the Indian meal moth, Plodia interpunctella (Hübner) (Lepidoptera: Pyralidae) to oils. Ann. Entomol. Soc. Am. 96, 524–531 (2003).Article 

Google Scholar 
Hagstrum, D. W. Using five sampling methods to measure insect distribution and abundance in bins storing wheat. J. Stored Prod. Res. 36, 253–262 (2000).CAS 
PubMed 
Article 

Google Scholar 
Athanassiou, C. G., Kavallieratos, N. G., Sciarretta, A. & Trematerra, P. Mating disruption of Ephestia kuehniella (Zeller) (Lepidoptera: Pyralidae) in a storage facility: spatio-temporal distribution changed after long-term application. J. Stored Prod. Res. 67, 1–12 (2016).Article 

Google Scholar 
Lee, W. H., Jung, J. M., Kim, J., Lee, H. & Jung, S. Analysis of the spatial distribution and dispersion of Plodia interpunctella (Lepidoptera: Pyralidae) in South Korea. J. Stored Prod. Res. 86, 101577 (2020).Article 

Google Scholar 
Gerken, A.R. & Campbell, J.F. Spatial and temporal variation in stored-product insect pest distributions and implications for pest management in processing and storage facilities. Ann. Entomol. Soc. Am. saab049(2021).Athanassiou, C. G. & Buchelos, CTh. Detection of stored-wheat beetle species and estimation of population density using unbaited probe traps and grain trier samples. Ent. Exp. et Applic. 98, 67–78 (2001).Article 

Google Scholar 
Subramanyam, B. & Hagstrum, D.W. Sampling. In: Subramanyam B. & Hagstrum D.W. (eds), Integrated Management of Insects in Stored Products. Marcel Dekker Inc., pp. 135–193 (1995).Morrison, W. R. et al. Aeration to manage insects in wheat stored in the Balkan peninsula: Computer simulations using historical weather data. Agronomy 10, 1927 (2020).Article 

Google Scholar 
Toews, M. D., Campbell, J. F. & Arthur, F. H. Temporal dynamics and response to fogging or fumigation of stored-product Coleoptera in a grain processing facility. J. Stored Prod. Res. 42, 480–498 (2006).Article 

Google Scholar 
Buckman, K. A., Campbell, J. F. & Subramanyam, B. Tribolium castaneum (Coleoptera: Tenebrionidae) associated with rice mills: Fumigation efficacy and population rebound. J. Econ. Entomol. 106, 499–512 (2013).PubMed 
Article 

Google Scholar 
Campbell, J. F., Buckman, K. A., Fields, P. G. & Subramanyam, Bh. Evaluation of structural treatment efficacy against Tribolium castaneum and Tribolium confusum (Coleoptera: Tenebrionidae) using meta-analysis of multiple studies conducted in food facilities. J. Econ. Entomol. 108, 2125–2140 (2015).CAS 
PubMed 
Article 

Google Scholar 
Levene, H. Robust tests for equality of variances. In Ingram Olkin; Harold Hotelling; et al. (eds.). Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling. Stanford University Press. pp. 278–292(1960).SAS Institute. SAS/STAT 9.2 User’s guide. SAS Institute (2008).Taylor, L. R. Aggregation, variance and mean. Nature 189, 732–735 (1961).ADS 
Article 

Google Scholar 
Iwao, S. A new method of sequential sampling to classify populations according to a critical density. Res. Popln. Ecol. 16, 281–288 (1975).
Google Scholar 
Green, R. H. Measurement of non-randomness in spatial distribution. Res. Popln. Ecol. 8, 1–17 (1966).
Google Scholar 
Hillhouse, T. L. & Pitre, H. N. Comparison of sampling techniques to obtain measurements of insect populations on soybeans. J. Econ. Entomol. 67, 411–414 (1974).Article 

Google Scholar 
Cassie, R. M. Frequency distribution models in the ecology of plankton and other organisms. J. Anim. Ecol. 31, 65–92 (1962).Article 

Google Scholar 
Southwood, T. R. E. Ecological Methods, with Particular Reference to the Study of Insect Population (Chapman and Hall, 1995).
Google Scholar 
Costa, M. G., Barbosa, J. C., Yamamoto, P. T. & Leal, R. M. Spatial distribution of Diaphorina citri Kuwayama (Hemiptera: Psyllidae) in citrus orchards. Scientia Agric 67, 546–554 (2010).Article 

Google Scholar 
Patil, G. P. & Stiteler, W. M. Concepts of aggregation and their quantification: A critical review with some new result and applications. Pop. Ecol. 15, 238–254 (1974).Article 

Google Scholar 
David, F. N. & Moor, P. G. Notes on contagious distribution in plant populations. Ann. Bot. 18, 47–53 (1954).Article 

Google Scholar 
Lloyd, M. Mean crowding. J. Anim. Ecol. 36, 1–30 (1967).Article 

Google Scholar 
Southwood, T. R. E. & Henderson, P. A. Ecological Methods 3rd edn. (Blackwell Sciences, 2000).
Google Scholar 
Feng, M. G. & Nowierski, R. M. Spatial distribution and sampling plans for four species of cereal aphids (Homoptera: Aphididae) infesting spring wheat in southwestern Idaho. J. Econ. Entomol. 85, 830–837 (1992).Article 

Google Scholar  More

We focused on the Northern High Latitudes (NHL, latitude > 50°N, excluding Greenland) due to their importance for carbon (CO2-C, the same hereafter)-climate feedbacks in the Earth system. To minimize the potential human influence on the CO2 cycle, we excluded areas under agricultural management (croplands, cropland/natural vegetation mosaic, and urban types), and considered only pixels of natural vegetation defined from the MODIS MCD12Q1 (v006) based IGBP land cover classification. Our main focus was the NHL permafrost region because permafrost plays a critical role in the ecology, environment, and society in the NHL. Permafrost, or permanently frozen ground, is defined as ground (soil, sediment, or rock) that remains at or below 0 °C for at least two consecutive years. The occurrence of permafrost is primarily controlled by temperature and has a strong effect on hydrology, soils, and vegetation composition and structure. Based on the categorical permafrost map from the International Permafrost Association58, the permafrost region (excluding permanent snow/ice and barren land), including sporadic (10–50%), discontinuous (50–90%), and continuous ( >90%) permafrost, encompasses about 15.7 × 106 km2, accounts for 57% of the NHL study dominion, and is dominated by tundra (shrubland and grass) and deciduous needleleaf (i.e., larch) forest that is regionally abundant in Siberia. The NHL non-permafrost region covers about 11.9 × 106 km2 and is dominated by mixed and evergreen needleleaf boreal forests (Fig. S1).Atmospheric CO2 inversions (ACIs)ACIs provide regionally-integrated estimates of surface-to-atmosphere net ecosystem CO2 exchange (NEEACI) fluxes by utilizing atmospheric CO2 concentration measurements and atmospheric transport models59. ACIs differ from each other mainly in their underlying atmospheric observations, transport models, spatial and temporal flux resolutions, land surface models used to predict prior fluxes, observation uncertainty and prior error assignment, and inversion methods. We used an ensemble mean of six different ACI products, each providing monthly gridded NEEACI at 1-degree spatial resolution, including Carbon‐Tracker 2019B (2000-2019, CT2019)60, Carbon‐Tracker Europe 2020 (2000–2019, CTE2020)61, Copernicus Atmosphere Monitoring Service (1979–2019, CAMS)62, Jena CarboScope (versions s76_v4.2 1976–2017, and s85_v4.2 1985-2017)63,64, and JAMSTEC (1996–2017)65. The monthly gridded ensemble mean NEEACI at 1-degree spatial resolution was calculated using the available ACIs from 1980-2017. Monthly ACI ensemble mean NEEACI data were summed to seasonal and annual values, and used to calculate the spatial and temporal trends of net CO2 uptake, and to investigate its relationship to climate and environmental controls.Productivity datasetDirect observations of vegetation productivity do not exist at a circumpolar scale. We therefore used two long-term gridded satellite-based estimates of vegetation productivity, including gross primary production (GPP) derived using a light use efficiency (LUE) approach (LUE GPP, 1982–1985)21,66 and satellite observations of Normalized Difference Vegetation Index (NDVI) from the Global Inventory Modeling and Mapping Studies (GIMMS NDVI, 1982–1985)67. LUE GPP (monthly, 0.5° spatial resolution, 1982–2015) is calculated from satellite observations of NDVI from the Advanced Very High-Resolution Radiometer (AVHRR; 1982 to 2015) combined with meteorological data, using the MOD17 LUE approach. LUE GPP has been extensively validated with a global array of eddy-flux tower sites68,69,70 and tends to provide better estimates in ecosystems with greater seasonal variability at high latitudes. Following66,71, we used the ensemble mean of GPP estimates from three of the most commonly used meteorological data sets: National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis; NASA Global Modeling and Assimilation Office (GMAO) Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2); and European Center for Medium-Range Weather Forecasting (ECMWF). GIMMS NDVI (bimonthly, 1/12 spatial resolution, 1982–2015) provides the longest satellite observations of vegetation “greenness”, and is widely used in studies of phenology, productivity, biomass, and disturbance monitoring as it has proven to be an effective surrogate of vegetation photosynthetic activity72.The gridded GPP data were resampled to 1-degree resolution at monthly time scales, to be consistent with NEEACI, and used to test (H1) whether greater temperature sensitivity of vegetation productivity explains the different trends in net CO2 uptake across the NHL. LUE GPP was also used to calculate monthly total ecosystem respiration (TER) as the difference between GPP and NEEACI (i.e., TERresidual =  GPP– NEEACI) from 1982-2015, as global observations of respiration do not exist. The NEEACI, GPP and TERresidual were used as observation-constrained top-down CO2 fluxes to investigate mechanisms underlying the seasonal CO2 dynamics in the structural equation modeling and additional decision tree-based analysis.Eddy Covariance (EC) measurements of bottom-up CO2 fluxesA total of 48 sites with at least three years of data representing the major NHL ecosystems were obtained from the FLUXNET2015 database (Table S1 and Fig. S1). EC measurements provide direct observations of net ecosystem CO2 exchange (NEE) and estimate the GPP and TER flux components of NEE using other climate variables. Daily GPP and TER were estimated as the mean value from both the nighttime partitioning method73 and the light response curve method74. More details on the flux partitioning and gap-filling methods used are provided by75. Daily fluxes were summed into seasonal and annual values and used to compare with trends from ACIs (Fig. S7), to estimate the climate and environmental controls on the CO2 cycle in the pathway analysis (Fig. 5), and to calculate the net CO2 uptake sensitivity to spring temperature (Fig. S14).Ensemble of dynamic global vegetation models (TRENDY simulations)The TRENDY intercomparison project compiles simulations from state-of-the-art dynamic global vegetation models (DGVMs) to evaluate terrestrial energy, water, and net CO2 exchanges76. The DGVMs provide a bottom-up approach to evaluate terrestrial CO2 fluxes (e.g., net biome production [NBP]) and allow deeper insight into the mechanisms driving changes in carbon stocks and fluxes. We used monthly NBP, GPP, and TER (autotrophic + heterotrophic respiration; Ra + Rh) from ten TRENDY v7 DGVMs76, including CABLE-POP, CLM5.0, OCN, ORCHIDEE, ORCHIDEE-CNP, VISIT, DLEM, LPJ, LPJ-GUESS, and LPX. We analyzed the “S3” simulations that include time-varying atmospheric CO2 concentrations, climate, and land use. All simulations were based on climate forcing from the CRU-NCEPv4 climate variables at 6-hour resolution. CO2 flux outputs were summarized monthly at 1-degree spatial resolution from 1980 to 2017. Monthly ensemble mean NBP, GPP, and TER were summed to seasonal and annual values, and then used to compare with observation-constrained ACI top-down CO2 fluxes (Figs. 4 and 5).Satellite data-driven carbon flux estimates (SMAP L4C)We also used a much finer spatio-temporal simulation of carbon fluxes from the NASA Soil Moisture Active Passive (SMAP) mission Level 4 Carbon product (L4C) to quantify the temperature and moisture sensitivity of NHL CO2 exchange77. The SMAP L4C provides global operational daily estimates of NEE and component CO2 fluxes for GPP and TER at 9 km resolution since 2015; whereas, an offline version of the L4C model provides a similar Nature Run (NR) carbon flux record over a longer period (2000-present), but without the influence of SMAP observational inputs. The L4C model has been calibrated against FLUXNET tower CO2 flux measurements and shows favorable performance and accuracy in high latitude regions4,77. In this analysis, daily gridded CO2 fluxes at 9-km resolution from the L4C NR record were summed to seasonal and annual values, and used to calculate the sensitivity of net C uptake in response to spring temperature (Fig. S14).CO2 fluxes in this analysis are defined with respect to the biosphere so that a positive value indicates the biosphere is a net sink of CO2 absorbed from the atmosphere. The different data products described above use different terminology (e.g., NEE, NBP) with slightly different meanings; however, they all provide estimates of net land-atmosphere CO2 exchange78.Climate, tree cover, permafrost, and soil moisture dataMonthly gridded air temperatures at 0.5-degree spatial resolution from 1980 to 2017 were obtained from the Climate Research Unit (CRU TS v4.02) at the University of East Anglia79. Air temperature was summarized at seasonal and annual scales to calculate temperature sensitivities of net CO2 uptake and to investigate the mechanism underlying the seasonal CO2 dynamics.Percent tree cover (%TC) at 0.05-degree spatial resolution was averaged over a 35-year (1982-2016) period using annual %TC layers derived from the Advanced Very High-Resolution Radiometer (AVHRR) (Fig. 1a)42. %TC was binned using 5% TC intervals to assess its relation to net CO2 uptake, or aggregated at a regional scale (e.g., TC  > 50% or TC  90%), discontinuous permafrost (DisconP, 10% < P  90%), discontinuous (DisconP, 10% < P  0.05 indicate a good fitting model), Bentler’s comparative fit index (CFI, where CFI ≈ 1 indicates a good fitting model), and the root mean square error of approximation (RMSEA; where RMSEA ≤ 0.05 and p  > 0.1 indicate a good fitting model). The standardized regression coefficient can be interpreted as the relative influences of exogenous (independent) variables. The R2 indicates the total variation in an endogenous (dependent) variable explained by all exogenous (independent) variables.Direct and legacy effects of temperature on seasonal net CO2 uptakeBecause landscape thawing and snow conditions regulate the onset of vegetation growth and influence the seasonal and annual CO2 cycles in the NHL24,84, we also analyzed the legacy effects of spring (May–Jun) temperature on seasonal net CO2 uptake. We regressed seasonal and annual net CO2 uptake from the site-level EC observations, regional-level ACI ensemble, and the TRENDY NBP ensemble against spring (May-June) air temperature. For EC observations, net CO2 uptake (i.e., NEE) and air temperature were summarized from site-level measurements. For the ACIs and TRENDY ensemble, net CO2 uptake (i.e., NEEACI and NBP) was summarized as regional means from the ACIs and TRENDY ensemble outputs, and air temperature was summarized as regional means from CRU temperature. The slope of the regression line was interpreted as the spring temperature sensitivity of the CO2 cycle. Simple linear regression was used here mainly due to the strong influence of spring temperature on the seasonal and annual CO2 cycle in NHL ecosystems30. Temperature sensitivity (γ: g C m−2 day−1 K−1) is the change in net CO2 flux (g C m−2 day−1) in response to a 1-degree temperature change. The sensitivity of net CO2 uptake to warm spring anomalies was calculated for different seasons (EGS, LGS, and annual) and regions (i.e., permafrost and non-permafrost), and the T-test was used to test for the difference in γ among different regions, seasons, and datasets. Similarly, direct effects of temperature on net CO2 uptake were calculated using the same season data (Fig. S14).Observationally-constrained estimates (EC and ACIs) showed that the sensitivity of net CO2 uptake in the EGS to spring temperature is positive (γ  > 0) and not statistically different (p  > 0.05) between permafrost and non-permafrost regions (({gamma }_{{ACI}}^{{np}})=0.125 ± 0.020 gC m−2 d−1 K−1; ({gamma }_{{EC}}^{{np}}) = 0.052 ± 0.013 gC m−2 d−1 K−1). In contrast, the sensitivity of net CO2 uptake in LGS to spring temperature is negative (γ  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

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 areasThe study took place in the cities of Basel, Lausanne, Lugano and Zurich, in Switzerland. Basel, Lausanne and Zurich are located north of the Alps, in the geographical region of the Central Plateau (Supplementary Fig. S1). This region stretches from Lake Geneva in the southwest to Lake Constance in the northeast and is the most densely populated region in Switzerland. Zurich is the largest city of Switzerland and encompasses 88 km2 with a total human resident population of 420,21741. Lausanne and Basel are smaller than Zurich, with a surface of 41 and 24 km2 and a total population of 139,408 and 173,232, respectively41. The climate in these three cities is moderately continental, with cold winters often reaching freezing temperatures in January, and warm summers. Lugano is located in Ticino, south of the Alps (Supplementary Fig. S1), where the climate is strongly affected by the Mediterranean Sea, with mild winters and summers warm and humid, sometimes hot. Lugano is the smallest of the four cities with 50,603 residents in 26 km241.Aedes albopictus is well established in Lugano since 2009 and an integrated vector management is constantly implemented to contain the numbers of the mosquito at a manageable level. This consists of an intensive surveillance, with oviposition traps distributed according to a grid system, several control interventions, such as the removal of breeding sites and the systematic application of larvicides in public areas, mainly in catch basins, and extensive public information campaigns24,26. In Basel, two populations of Ae. albopictus are established since 2018: a first population in an area adjacent to the motorway toll on the border with France and a second population in an area near the border with Germany27. The mosquito has also been recorded repeatedly at various locations in the city of Basel and the surveillance indicates that the mosquito is spreading42. Control actions are taken exclusively within the perimeter of repeated detections of the mosquito and include regular treatment of catch basins with larvicides, distribution of flyers and door-to-door information campaigns42. In Zurich, Ae. albopictus was first detected in 2016 in a bus station for international coach services located in the centre of the city, near the main train station. Thanks to immediate surveillance and control actions (i.e., treatment of catch basins in the area with larvicides), to date there is no established population within the perimeter of the bus station despite continuous repeated introductions40. A small population was also detected in 2018 in a suburban neighbourhood in the Wollishofen district of Zurich, approximately 5 km southwest from the international bus station. Also in this case, immediate surveillance and control actions, including larval control and door-to-door information, were taken with success and no adults, eggs or aquatic stages have been found in 2020 and 202140. In Lausanne, no tiger mosquito has been reported to date (Swiss Mosquito Network, http://www.mosquitoes-switzerland.ch (accessed on 17 February 2022)).Microclimate dataBased on a previous investigation we conducted in Ticino, Basel and Zurich20, we focused the microclimate monitoring on ordinary stormwater catch basins positioned on the side of public roads. In each city, we monitored ten catch basins located either in urban context (defined as areas with high-density development, consisting of apartment blocks, commercial or industrial units) or in residential areas consisting mainly of houses with private gardens located in peri-urban area (Supplementary Table S1, Supplementary Fig. S2). The catch basins were usually homogeneous in dimension, in the same city, although we recorded variations in depth. In Basel, we included catch basins located in the urban area near the border with France, in which Ae. albopictus is established. In Zurich, we included catch basins located in the international bus station, where Ae. albopictus was recorded in summer, and in the residential area of Wollishofen, where a small population of Ae. albopictus was detected and then likely eradicated. In Lausanne, some catch basins were selected in potential points of introduction of the mosquito (e.g., near a campsite, the main train station, etc.). In Lugano, Ae. albopictus was established in all the locations selected.A sensor device was installed in each selected catch basin. The sensor devices were built in house. The development of the devices and the Wireless Sensor Network (WSN) has been described in detail by Strigaro et al.29. Briefly, the device consisted of a waterproof plastic box containing a LoPy Micro-Controller Unit (Pycom, Guildford, United Kingdom), a waterproof temperature probe (accuracy of ± 0.5 °C), a light sensor (measuring illuminance arriving at the sensor device, in lux), an SD card, the rechargeable batteries and other parts. The main box, with the light sensor, was hung on the inside wall of the catch basin. The temperature probe was attached to the wall at a depth ranging from 0.3 to 0.5 m, depending on the depth of the catch basin and the level of the water in the catch basin. The probe was placed in direct contact with the inside wall of the catch basin, in order to measure the microclimatic conditions where the mosquito eggs are potentially laid. The data collected was transmitted to a data warehouse based on istSOS, an open-source Python based implementation of the Sensor Observation Service standard (SOS) of the Open Geospatial Consortium (OGC)43. The data was transmitted through the Swisscom Low Power Network (LPN) LoRaWAN (Swisscom Ltd, Ittigen, Switzerland): the data sent by the sensor devices was received by a Swisscom Gateway and then sent to the data warehouse29.In addition to the sensor devices installed in the catch basins, four devices were installed outside four catch basins in each city, except in Lugano, where three devices were installed. These external devices were placed in vegetation representing potential resting habitats for Ae. albopictus adults in the reproductive season, at 1–2 m above the ground and analyzed to confirm the close similarity between measured external temperatures and MeteoSwiss gridded temperature data. However, since the main goal of the data collection was to model the differences between MeteoSwiss gridded temperature data and catch basins’ temperatures, only a small number of external sensors were deployed. Microclimate data were collected from beginning of December 2019 to end of February 2020, a period defined as cold season, with acquisition interval set at one hour. In Lugano, data collection started on the 12th or 13th of December 2019.Local climate dataWe used two types of local climate data. The first type is the momentary hourly free-air temperatures recorded at 2 m above ground level by permanent weather stations. The weather stations belong to SwissMetNet, the automatic monitoring network of MeteoSwiss. For each city, we selected the weather station closest to the study area (Supplementary Table S1, Supplementary Fig. S2) and temperature data were retrieved from https://gate.meteoswiss.ch/idaweb (source: MeteoSwiss, Zurich-Airport, Switzerland; accessed on 12 August 2021).The second type of local climate data is the MeteoSwiss spatial climate daily datasets (source: MeteoSwiss). These temperature datasets are constructed through interpolation of daily minimum, maximum, and mean temperatures from a network of approximately 90 SwissMetNet permanent weather stations to a 1 km resolution grid in the Swiss coordinate system CH190344,45. This results in three temperature datasets describing the km-scale distribution of day-to-day temperature variations in Switzerland. We referred to them as gridded temperature data. Each monitored catch basin and external device was assigned, based on its geographical position, to the corresponding 1 km × 1 km cell of the climate grid. Each cell was identified with its MeteoSwiss (MS) number (Supplementary Table S1).Data analysisThe hourly temperatures were used to compute daily mean, maximum and minimum temperatures and daily temperature ranges, which were calculated as the difference between the maximum and minimum daily temperature. Temperatures of catch basins and external habitats were compared to temperatures of permanent weather stations and to the gridded temperatures both graphically and using the nonparametric Mann–Whitney U-test, for which a P value of  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

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