Philippa Kaur

Hauer, M. E. et al. Sea-level rise and human migration. Nat. Rev. Earth Environ. 1, 28–39 (2020).ADS 
Article 

Google Scholar 
Duy, P., Chapman, L., Tight, M., Thuong, L. & Linh, P. Urban resilience to floods in coastal cities: Challenges and opportunities for Ho Chi Minh city and other emerging cities in southeast Asia. J. Urban Plan. Dev. 144, 05017018 (2018).Article 

Google Scholar 
Magno, R. et al. Semi-automatic operational service for drought monitoring and forecasting in the Tuscany region. Geosciences 8, 49 (2018).ADS 
Article 

Google Scholar 
Rico, A., Olcina, J., Baños, C., Garcia, X. & Sauri, D. Declining water consumption in the hotel industry of mass tourism resorts: Contrasting evidence for Benidorm, Spain. Curr. Issues Tour. 23, 770–783 (2020).Article 

Google Scholar 
Hasnat, G. T., Kabir, M. A. & Hossain, M. A. Major environmental issues and problems of South Asia, particularly Bangladesh. Handb. Environ. Mater. Manag., 1–40 (2018).Neumann, B., Vafeidis, A. T., Zimmermann, J. & Nicholls, R. J. Future coastal population growth and exposure to sea-level rise and coastal flooding—A global assessment. PLoS One 10, e0118571 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Cao, A. et al. Future of Asian Deltaic Megacities under sea level rise and land subsidence: Current adaptation pathways for Tokyo, Jakarta, Manila, and Ho Chi Minh City. Curr. Opin. Environ. Sustain. 50, 87–97 (2021).Article 

Google Scholar 
Rahmasary, A. N. et al. Overcoming the challenges of water, waste and climate change in Asian cities. Environ. Manag. 63, 520–535 (2019).ADS 
Article 

Google Scholar 
Smol, M., Adam, C. & Preisner, M. Circular economy model framework in the European water and wastewater sector. J. Mater. Cycles Waste Manag. 22, 682–697 (2020).Article 

Google Scholar 
Islam, M. F., Bhattacharya, B. & Popescu, I. Flood risk assessment due to cyclone-induced dike breaching in coastal areas of Bangladesh. Nat. Hazards Earth Syst. Sci. 19, 353–368 (2019).ADS 
Article 

Google Scholar 
Salim, M. A. & Siswanto, A. B. Kajian Penanganan Dampak Banjir Kabupaten Pekalongan. Rang Tek. J. 4, 295–303 (2021).Article 

Google Scholar 
Endo, A. et al. Describing and visualizing a water–energy–food nexus system. Water 10, 1245 (2018).Article 

Google Scholar 
Gurdak, J. J., Geyer, G. E., Nanus, L., Taniguchi, M. & Corona, C. R. Scale dependence of controls on groundwater vulnerability in the water–energy–food nexus, California Coastal Basin aquifer system. J. Hydrol. Reg. Stud. 11, 126–138 (2017).Article 

Google Scholar 
Lu, J., Lin, Y., Wu, J. & Zhang, C. Continental-scale spatial distribution, sources, and health risks of heavy metals in seafood: Challenge for the water-food-energy nexus sustainability in coastal regions?. Environ. Sci. Pollut. Res. 28, 63815–63828 (2021).CAS 
Article 

Google Scholar 
Miller-Robbie, L., Ramaswami, A. & Amerasinghe, P. Wastewater treatment and reuse in urban agriculture: Exploring the food, energy, water, and health nexus in Hyderabad, India. Environ. Res. Lett. 12, 075005 (2017).ADS 
Article 
CAS 

Google Scholar 
Taniguchi, M., Endo, A., Gurdak, J. J. & Swarzenski, P. Water-energy-food nexus in the Asia-Pacific region. J. Hydrol. 11, 1–8 (2017).
Google Scholar 
Bahri, M. Analysis of the water, energy, food and land nexus using the system archetypes: A case study in the Jatiluhur reservoir, West Java, Indonesia. Sci. Total Environ. 716, 137025 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lubis, R., Delinom, R., Martosuparno, S. & Bakti, H. Water-Food Nexus in Citarum Watershed, Indonesia Vol. 118, 012023 (IOP Publishing, 2018).
Google Scholar 
Pawitan, H., Delinom, R. & Taniguchi, M. The human–environment sustainability in Indonesia: The case of the Citarum basin Vol. 23 (UNESCO-IHP, 2015).Carmichael, L. et al. Urban planning as an enabler of urban health: Challenges and good practice in England following the 2012 planning and public health reforms. Land Use Policy 84, 154–162 (2019).Article 

Google Scholar 
World Health Organization. Addressing the Social Determinants of Health: The Urban Dimension and the Role of Local Government (World Health Organization, 2012).
Google Scholar 
Trencher, G. & Karvonen, A. Stretching, “smart”: Advancing health and well-being through the smart city agenda. Local Environ. 24, 610–627 (2019).Article 

Google Scholar 
Yang, L. et al. Can an island economy be more sustainable? A comparative study of Indonesia, Malaysia, and the Philippines. J. Clean. Prod. 242, 118572 (2020).Article 

Google Scholar 
Choirunisa, A. K. & Giyarsih, S. R. Kajian Kerentanan Fisik, Sosial, dan Ekonomi Pesisir Samas Kabupaten Bantul Terhadap Erosi Pantai. J. Bumi Indones. 5 (2016).Gumay, A. Validity and reliability maritime English seafarers proficiency test. INFERENCE J. Engl. Lang. Teach. 3, 64–69 (2021).Article 

Google Scholar 
Tarigan, M. S. Perubahan garis pantai di wilayah pesisir perairan Cisadane, Provinsi Banten. Makara J. Sci. (2010).Pruss-Ustun, A., Corvalán, C. F., World Health Organization. Preventing Disease Through Healthy Environments: Towards an Estimate of the Environmental Burden of Disease (World Health Organization, 2006).
Google Scholar 
Baasanjargal, T., Soon-Joo, A. & Mi-Jeong, K. Comparative analysis of Indonesian Batik traditional patterns: Focused on patterns of Yogyakarta and Pekalongan in Java Island. 한복문화 22, 75–91 (2019).Rismawati, S. D., Sofiani, T. & Rahmawati, D. R. Legal culture of religious capitalism on Batik business (a case study in Pekalongan Indonesia). JL Pol. Glob. 33, 107 (2015).
Google Scholar 
Pekalongan, B. K. Kota Pekalongan dalam Angka 2021 (2021).Google Maps. Pekalongan, Central Java (2022).Sunarjo, W. A., Ilmiani, A. & Ardianingsih, A. Analisis SWOT Sebagai Pengembangan UMKM Berbasis Ekonomi Kreatif Destinasi Pariwisata Batik Kota Pekalongan. Pena J. Ilmu Pengetah. Dan Teknol. 33, 34–43 (2019).Article 

Google Scholar 
Perpustakaan Provinsi Jawa Tengah. Museum Batik Pekalongan (2017).Brzezina, N. et al. Development of organic farming in Europe at the crossroads: Looking for the way forward through system archetypes lenses. Sustainability 9, 821 (2017).Article 

Google Scholar 
Gillies, A. & Maliapen, M. Using healthcare system archetypes to help hospitals become learning organisations. J. Model. Manag. (2008).Braun, W. The System Archetypes. The Systems Modeling Workbook, 1–26 (2002).Sterman, J. System Dynamics: Systems thinking and modeling for a complex world (2002).Islam, M. & Raja, D. R. Waterlogging risk assessment: An undervalued disaster risk in coastal urban community of Chattogram, Bangladesh. Earth 2, 151–173 (2021).Article 

Google Scholar 
Brzezina, N., Kopainsky, B. & Mathijs, E. Can organic farming reduce vulnerabilities and enhance the resilience of the European food system? A critical assessment using system dynamics structural thinking tools. Sustainability 8, 971 (2016).Article 

Google Scholar 
Nguyen, N. C. & Bosch, O. J. A systems thinking approach to identify leverage points for sustainability: A case study in the Cat Ba Biosphere Reserve, Vietnam. Syst. Res. Behav. Sci. 30, 104–115 (2013).Article 

Google Scholar 
Maani, K. E. & Cavana, R. Y. Systems Thinking, System Dynamics: Managing Change and Complexity (Pearson Prentice Hall, 2007).
Google Scholar 
Braun, W. The System Archetypes-the Systems Modeling Workbook. Available Wwwu Uniklu Ac Atgossimitpapsdwbsysarch Pdf (2002).Senge, P. M. The Fifth Discipline: The Art and Practice of the Learning Organization (Currency, 2006).
Google Scholar 
Bahri, M. et al. Deliverable 3.3: Integrated model with ad-hoc systems model of urban water supply (2018).Pekalongan, B. P. P. D. K. Pekalongan dalam Angka (2021).Fajar, M., Mediani, A. & Finesa, Y. Analisis Peranan IPAL dalam Strategi Penanganan Limbah Industri Batik di Kota Pekalongan. in Prosiding Seminar Nasional Geografi UMS X 2019 (2019).Kartika, F. D. S. & Helmi, M. Meta-analysis of Community’s Adaptation Pattern with Tidal Flood in Pekalongan City, Central Java, Indonesia Vol. 125, 09001 (EDP Sciences, 2019).
Google Scholar 
Kartika, F. D. S., Helmi, M. & Amirudin, A. Analisis Perubahan Penggunaan Lahan di Wilayah Pesisir Kota Pekalongan Menggunakan Citra Lansat 8, vol. 1 (2019).Damayanti, M. & Latifah, L. Strategi Kota Pekalongan dalam pengembangan wisata kreatif berbasis industri batik. J. Pengemb. Kota 3, 100–111 (2017).Article 

Google Scholar 
Pekalongan, B. K. Kota Pekalongan dalam Angka 2002 (2002).Andreas, H., Abidin, H. Z., Sarsito, D. A. & Pradipta, D. Adaptation of ‘early climate change disaster’ to the Northern coast of Java Island Indonesia. Eng. J. 22, 207–219 (2018).Article 

Google Scholar 
Marfai, M. A. et al. The impact of tidal flooding on a coastal community in Semarang, Indonesia. Environmentalist 28, 237–248 (2008).Article 

Google Scholar 
Chaussard, E., Amelung, F., Abidin, H. & Hong, S.-H. Sinking cities in Indonesia: ALOS PALSAR detects rapid subsidence due to groundwater and gas extraction. Remote Sens. Environ. 128, 150–161 (2013).ADS 
Article 

Google Scholar 
Andreas, H., Abidin, H. Z., Sarsito, D. A. & Pradipta, D. Remotes Sensing Capabilities on Land Subsidence and Coastal Water Hazard and Disaster Studies Vol. 500, 012036 (IOP Publishing, 2020).
Google Scholar 
Shofiana, R., Subardjo, P. & Pratikto, I. Analisis perubahan penggunaan lahan di wilayah pesisir Kota pekalongan menggunakan data landsat 7 etm+. J. Mar. Res. 2, 35–43 (2013).
Google Scholar 
Wijaya, A. Analisis Dinamika Pola Spasial Penggunaan Lahan Pada Wilayah Terdampak Kenaikan Muka Air Laut di Kota Pekalongan (2017).El-Fath, D. D. I., Atmodjo, W., Helmi, M., Widada, S. & Rochaddi, B. Analisis Spasial Area Genangan Banjir Rob Setelah Pembangunan Tanggul di Kabupaten Pekalongan, Jawa Tengah. Indones. J. Oceanogr. 4, 96–110 (2022).
Google Scholar 
Novita, M. G., Helmi, M., Widiaratih, R., Hariyadi, H. & Wirasatriya, A. Mengkaji Area Genangan Banjir Pasang Terhadap Penggunaan Lahan Pesisir Tahun 2020 Menggunakan Metode Geospasial di Kabupaten Pekalongan, Provinsi Ja. Indones. J. Oceanogr. 3, 14–26 (2021).Article 

Google Scholar 
Salim, M. A. Penanganan Banjir dan Rob di Wilayah Pekalongan. J. Tek. Sipil 11, 15–23 (2018).
Google Scholar 
Jumatiningrum, N. & Indrayati, A. Strategi Adaptasi Masyarakat Kelurahan Bandengan Kecamatan Pekalongan Utara dalam Menghadapi Banjir Pasang Air Laut (Rob). Edu Geogr. 9, 136–143 (2021).
Google Scholar 
BNPB. Data Kebencanaan Nasional (BNPB, 2021).
Google Scholar 
Giampietro, M., Aspinall, R. J., Ramos-Martin, J. & Bukkens, S. G. Resource Accounting for Sustainability Assessment: The Nexus Between Energy, Food, Water and Land Use (Routledge, 2014).Book 

Google Scholar 
Meadows, D. H., Randers, J. & Meadows, D. L. The Limits to Growth (1972) (Yale University Press, 2013).MATH 

Google Scholar 
Albrecht, T., Crootof, A. & Scott, C. Trends in the development of water–energy–food nexus methods (2017).Leck, H., Fitzpatrick, D. & Burchell, K. Energy, water and food: Towards a critical nexus approach. in Handbook on the Geographies of Energy (Edward Elgar Publishing, 2017).Scott, C. A., Kurian, M. & Wescoat, J. L. The water–energy–food nexus: Enhancing adaptive capacity to complex global challenges. in Governing the Nexus 15–38 (Springer, 2015).Wanty, E. E. Analisis Produksi Batik Cap Dari UKM Batik Kota Pekalongan (Studi Pada Sentra Batik Kota Pekalongan-Jawa Tengah, 2006).
Google Scholar 
Mankiw, N. G. Macroeconomics Vol. 41 (Worth Publishers, 2003).
Google Scholar 
Shen, J. & Kee, G. Development and Planning in Seven Major Coastal Cities in Southern and Eastern China (Springer, 2017).Book 

Google Scholar 
Xu, C., Haase, D., Su, M. & Yang, Z. The impact of urban compactness on energy-related greenhouse gas emissions across EU member states: Population density vs physical compactness. Appl. Energy 254, 113671 (2019).Article 

Google Scholar 
Marfai, M. A. & Cahyadi, A. Dampak bencana banjir pesisir dan adaptasi masyarakat terhadapnya di kabupaten Pekalongan (2017).Wartadesa.net. Tiga hari banjir rendam Pekalongan (2018).Google Maps. A dike in Pekalongan (n.d).Anindita, R. M., Susilowati, I. & Muhammad, F. Analisis Efektifitas Tanggul Laut di Pesisir Pekalongan Terhadap Penurunan Intensitas Banjir, vol. 2 80–88 (2020).Taniguchi, M. Groundwater and Subsurface Environments: Human Impacts in Asian Coastal Cities (Springer Science & Business Media, 2011).Book 

Google Scholar 
Baños, C. J., Hernández, M., Rico, A. M. & Olcina, J. The hydrosocial cycle in coastal tourist destinations in Alicante, Spain: Increasing resilience to drought. Sustainability 11, 4494 (2019).Article 

Google Scholar 
Sauda, R. H. & Nugraha, A. L. Kajian pemetaan kerentanan banjir rob di kabupaten pekalongan. J. Geod. Undip 8, 466–474 (2019).
Google Scholar 
Wartadesa.net. Ratusan warga Sragi masih mengungsi (2022).Buchori, I. et al. Adaptation to coastal flooding and inundation: Mitigations and migration pattern in Semarang City, Indonesia. Ocean Coast. Manag. 163, 445–455 (2018).Article 

Google Scholar 
Setiadi, R. & Nalau, J. Can urban regeneration improve health resilience in a changing climate? (2015).Isham, A., Mair, S. & Jackson, T. Wellbeing and productivity: A review of the literature (2020).Banson, K. E., Nguyen, N. C. & Bosch, O. J. Using system archetypes to identify drivers and barriers for sustainable agriculture in Africa: A case study in Ghana. Syst. Res. Behav. Sci. 33, 79–99 (2016).Article 

Google Scholar 
Lavrnić, S., Zapater-Pereyra, M. & Mancini, M. Water scarcity and wastewater reuse standards in Southern Europe: Focus on agriculture. Water. Air Soil Pollut. 228, 1–12 (2017).Article 
CAS 

Google Scholar 
Tortajada, C. & Nam Ong, C. Reused water policies for potable use (2016).Murali, R. M., Riyas, M., Reshma, K. & Kumar, S. S. Climate change impact and vulnerability assessment of Mumbai city, India. Nat. Hazards 102, 575–589 (2020).Article 

Google Scholar 
Abdullah, A. Y. M. et al. Spatio-temporal patterns of land use/land cover change in the heterogeneous coastal region of Bangladesh between 1990 and 2017. Remote Sens. 11, 790 (2019).ADS 
Article 

Google Scholar 
Ginanjar, A., Rezagama, A. & Handayani, D. S. Rencana Induk Sistem Penyediaan Air Minum Kota Pekalongan (2015).Reiblich, J., Hartge, E., Wedding, L., Killian, S. & Verutes, G. Bridging climate science, law, and policy to advance coastal adaptation planning. Mar. Policy 104, 125–134 (2019).Article 

Google Scholar 
Cook, B. I. et al. Revisiting the leading drivers of Pacific coastal drought variability in the contiguous United States. J. Clim. 31, 25–43 (2018).ADS 
Article 

Google Scholar 
Jodar-Abellan, A., Valdes-Abellan, J., Pla, C. & Gomariz-Castillo, F. Impact of land use changes on flash flood prediction using a sub-daily SWAT model in five Mediterranean ungauged watersheds (SE Spain). Sci. Total Environ. 657, 1578–1591 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Thanvisitthpon, N., Shrestha, S. & Pal, I. Urban flooding and climate change: A case study of Bangkok, Thailand. Environ. Urban. Asia 9, 86–100 (2018).Article 

Google Scholar 
Laksmi, G. S. Dampak Alih Fungsi Lahan dan Curah Hujan terhadap Banjir di Kota Pekalongan, Jawa Tengah, 382–391 (2020).Dhiman, R., VishnuRadhan, R., Eldho, T. & Inamdar, A. Flood risk and adaptation in Indian coastal cities: Recent scenarios. Appl. Water Sci. 9, 1–16 (2019).ADS 
Article 

Google Scholar 
Bahri, M. & Cremades, R. The Urban Drought Nexus Tool. Zenodo (2021). More

There are two phases involved in supporting the technical quality of the dataset: (i) the process used by Taronga when collecting data for each shark-bite incident, and (ii) the consistency modifications that we made during manuscript development.Phase oneFor each shark-bite incident, Taronga attempts to contact the most-relevant person involved in the event as possible — e.g., victim, victim’s family, witnesses. Those contacted are asked to complete a questionnaire with information relating to the shark bite. If applicable, questionnaires can also be completed by a fisheries officer in the relevant State and sent to Taronga. Taronga works closely with experts in each State’s fisheries department to validate information sourced in media reports. Each shark-bite case is unique, so the validation process varies depending on details specific to each incident. Forensic shark scientists within each State department are contacted after each shark bite to confirm details related to the incident. For example, validating species responsible for the bite often requires forensic analysis through expert examination of bite marks or artefacts21,22,23,24,25. When available, video footage is analysed for validation of information, such as confirmation of shark species and length.The database is cross-checked annually for the previous year with the International Shark-Attack File in Florida, USA, as well as with fisheries officers to ensure consistency. The database is cross-checked with the acknowledgement that there are discrepancies between versions due to differences in inclusion criteria. For example, the International Shark-Attack File includes bites where the victim was bitten aboard a boat, whereas the database we present here does not include bites aboard a boat.We acknowledge the limitations associated with this database, such as differences in reporting over time. For example, incidents might be reported more in recent decades due to technological advances making reporting more accessible or media publicising these events more widely18. There might also be reporting biases, for example, victims could be more likely to report a bite by a large, potentially dangerous shark (e.g., white, tiger, or bull shark) rather than a smaller, less-dangerous shark species (e.g., wobbegong shark). We also completed a quality assessment of the original database fields and redesigned the data acquisition and entry process (see Phase two) to allow exploration of shark-bite trends and patterns in Australia.Phase twoWe identified errors and inconsistencies in database fields. To avoid additional errors and inconsistencies, and to obtain a quality-controlled database, we redesigned the process for gaining and entering information into the database. This included creating a data descriptor (Supplementary File 2) used as a protocol to inform which questions to ask in the questionnaire. The data descriptor also directs the format of database entries by specifying information required in each field and by indicating the format of each entry (i.e., numeric, descriptive, or categorical). We manually inspected all previously entered data and adapted them to match the data descriptor. We checked each entry using the filter function in Microsoft Excel to identify any spelling and grammatical errors in the fields and ensure that all categories were grammatically identical. We standardised all metrics during this process (e.g., the data descriptor now stipulates that all length measurements should be recorded in metres).We validated and standardised the geographical locations of shark bites by converting all coordinates into decimal degrees using Microsoft Excel. We subsequently plotted all coordinates using the ggmap library27 in R (Version 4.0.2) (R Core Team 2020) (Fig. 2). We corrected any unusual coordinates (e.g., outside of Australian waters) and crosschecked them with site descriptions and states to ensure validity.Fig. 2Geographical locations of 1,196 shark bites in Australia. Each shark-bite incident is indicated by a red dot. (a) all shark-bite incidents; (b) bites most likely inflicted by bull sharks, (c) tiger sharks, and (d) white sharks. Two bite incidents that occurred at the Australian external territory, Cocos (Keeling) Islands, are not included on these maps. Background layers show elevation and major, perennial watercourses.Full size imageDuring the development of the data-descriptor protocol, we converted some previous descriptive columns into categorical columns. These columns included (but are not limited to) victim activity, attractant, injury location, injury severity, and weather condition. Converting these columns into categories facilitates analysis to investigate shark-bite patterns. For example, we converted victim activity into a categorical field to restrict answers to the following: snorkelling, motorised boating, unmotorised boating, boarding, swimming, standing, diving, fishing, or other, rather than allowing answers in any format. We used this information to create a time series to show the activity of shark-bite victims in Australia over time (Figs. 3 and 4). Shark bites have increased for boarders (including surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding) over time, particularly since 1960 (Figs. 3b, 4). This is likely due to the increase in popularity of board sports, particularly surfing, since the 1960s28. This trend is likely not reflected in Fig. 3a because shark bites are unlikely to be classed as ‘provoked’ during board riding.Fig. 3Number of shark bites (black, dashed line) and proportion of activity done by the victim at the time of shark-bite incidents in Australia from 1900 to 2022. Panels represent shark bites in Australia that are; (a) provoked or (b) unprovoked. Boarding includes surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding. Swimming includes snorkelling, spearfishing, freediving, body surfing, clinging to an object, falling into water, floating, or wading. Diving includes scuba-diving, hookah diving, or hard-hat diving. Fishing includes cleaning fish. No data for years 1908 and 1970.Full size imageFig. 4Number of shark bites (black, dashed line) and proportion of activity done by the victim at the time of shark-bite incidents in Australia from 1900 to 2022. Panels represent shark bites in Australia that are; (a) fatal or (b) non-fatal. Boarding includes surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding. Swimming includes snorkelling, spearfishing, freediving, body surfing, clinging to an object, falling into water, floating, or wading. Diving includes scuba-diving, hookah diving, or hard-hat diving. Fishing includes cleaning fish. No data for years 1908 and 1970.Full size imageThe column representing a shark-bite victim’s recovery status also requires a categorical response, restricting answers to: fatal, injured, or uninjured (Fig. 5). Since 1900, the proportion of shark-bite-related fatalities has decreased (Fig. 5). This trend is also true for the three species most attributed to shark-bite-related fatalities, white (Carcharodon carcharias), tiger (Galeocerdo cuvier), and bull (Carcharhinus leucas) sharks. The decrease in shark-bite-related deaths is likely due to advancements in medical responses to shark-bite victims over time14 and better understanding among surfers about using tourniquets to stem bleeding following increased certification as first responders in workplace occupational health and safety requirements. Bites resulting in an uninjured victim includes interactions where the shark might have bitten the victim’s equipment (i.e., surfboard, bodyboard, kayak) rather than biting the person.Fig. 5Proportion of victim-recovery status (fatal = grey; red = injured; blue = uninjured) resulting from all unprovoked shark bites in Australia between 1900–2022. Blank years represent years without any reported occurrences. (a) all shark-bite incidents, (b) bites most likely inflicted by bull sharks, (c) tiger sharks, or (d) white sharks. No data for years 1908 and 1970.Full size imagePreviously, entries in the injury location column were descriptive. We converted the column to a categorical field restricting answers to the following: arm, hand, lower arm, upper arm, shoulder, neck, head, torso, leg, foot, calf, thigh, pelvic region, or other. During the analyses process, we further categorised these injury locations into four body areas (head, arm, torso, and leg) to assess how injury location affects recovery status (fatal or injured) (Fig. 6). Fatality most often occurred following shark bites to the torso (Fig. 6). This is likely due to the injuries to organs and major arteries resulting in blood loss, which is a leading cause of shark-bite fatalities29. This is the first time that the location of a shark bite on the body has been assessed relative to recovery status.Fig. 6Proportion of Australian shark bites resulting in either fatality or injury categorised by injury location on the victim’s body (left panel; 250 bites resulting in fatal injury, 723 bites resulting in non-fatal injury) and by species (right panel; bites by 201 tiger, 170 bull, 258 white, and 303 bites other sharks).Full size imageUnderstanding how the location of a shark-bite wound relates to victim recovery has value in informing the development of shark-bite mitigations. For example, the development of shark-bite-resistant wetsuits could potentially result in higher survival rates of the user if the fabric is concentrated around the torso region30,31. Redesigning data acquisition and entry process to allow for categorical columns permits these types of analyses.Some detail of a shark-bite incident might be lost by converting previously descriptive columns into categorical columns. We addressed this by complementing categorical columns with accompanying fields to retain details of the incident. For example, injury location and injury severity columns are both categorical and allow for user-friendly data analysis, whereas the injury description column is descriptive and provides added detail about the victim’s injuries if applicable. Individually, all three columns address certain aspects of the victim’s injuries, and together, all three columns comprehensively summarise injuries to the shark-bite victim.Our analysis of the Australian Shark-Incident Database suggests that tiger sharks are proportionally responsible for the most fatalities of all shark species in Australia (38% of all tiger shark bites result in fatality), followed by bull sharks (32% of all bull shark bites result in fatality), and white sharks (25% of all white shark bites result in fatality) (Fig. 6). We emphasise that these figures represent the overall percentage of bites resulting in fatality since 1791 and do not account for possible changes over time. These figures are also proportional to the number of bites by each respective species. In Australia, white sharks are responsible for the largest number of bites on humans (361 total) compared to tiger (229 total) and bull sharks (197 total). At the time of publication, white and tiger sharks were each responsible for 91 and 86 total fatalities on humans in Australia, respectively.There were 540 incidents in which time of day was recorded. We used these data to assess whether particular shark-bite incidents are more likely to occur at specific times of day (Fig. 7). To standardise reported 24-hour times, we took the location (latitude, longitude) information and reported time and date of each incident using the getSunlightTimes function in the suncalc library in R32 to calculate whether the incident occurred in one of four light-availability categories: dawn, day, dusk, or night. Shark bites occur mostly during the day, which likely reflects time of day when there are more ocean users present. However, there is a slightly higher proportion of bites at dusk for bull sharks compared to tiger and white sharks (Fig. 7). Identifying these trends can assist authorities in developing data-driven educational messaging as a shark-bite mitigation measure. This is important considering that enhanced education is the preferred mitigation measure scored by ocean users in New South Wales33.Fig. 7Period of day (corrected for local time) distribution of provoked and unprovoked shark bites in Australia by shark species from 1791 to 2022 (n = 540 bites from all species, 70 from bull, 59 from tiger, and 217 from white sharks).Full size imageThese examples demonstrate that the Australian Shark-Incident Database will be useful for scientists to analyse environmental, social, and biological related shark-bite patterns in Australia. Use of the newly developed data descriptor to standardise future applications and account for quality assurance and control will aid in keeping the database consistent for ease of analysis and interpretation. Ultimately, the publishing of this database will improve our understanding of shark-bite incidents in Australia and will equip us with the knowledge to aim to avoid or predict these events in the future. More

Clemens, R. et al. Oats, more than just a whole grain: an introduction. Br. J. Nutr. 112, S1–S3 (2014).CAS 
PubMed 
Article 

Google Scholar 
Stewart, D. & Mcdougal, G. Oat agriculture, cultivation and breeding targets: implications for human nutrition and health. Br. J. Nutr. 2, 50–57 (2014).Article 
CAS 

Google Scholar 
Ren, C. Z. et al. “Twelfth Five-Year” Development Report of China’s Oat and Buckwheat Industry. Xi’an: Shaanxi Science and Technology Press, 2011–2015 (2016).Gleick, P. H. & Palaniappan, M. Peak water limits to freshwater withdrawal and use. Proc. Indian Natl. Sci. Acad. 107, 11155–11162 (2010).ADS 
CAS 

Google Scholar 
Yu, L., Zhao, X., Gao, X. & Siddique, K. H. M. Improving/maintaining water-use efficiency and yield of wheat by deficit irrigation: A global meta-analysis. Agric. Water Manag. 228, 105906 (2020).Article 

Google Scholar 
Bai, W., Zhang, H., Liu, B., Wu, Y. & Song, J. Effects of super-absorbent polymers on the physical and chemical properties of soil following different wetting and drying cycles. Soil Use Manag. 26, 253–260 (2010).Article 

Google Scholar 
Döll, P. Impact of climate change and variability on irrigation requirements: A global perspective. Clim. Change 54, 269–293 (2002).ADS 
Article 

Google Scholar 
Harris, F. et al. The water use of Indian diets and socio- demographic factors related to dietary blue water footprint. Sci. Total Environ. 587, 128–136 (2017).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Unesco. Water and jobs: Facts and figures. Perugia, Italy: UNESCO, World Water Assessment Program. Retrieved from http://unesdoc.unesco.org/images/0024/002440/244041e.pdf (2016).Landi, A. et al. Land suitability evaluation for surface, sprinkle and drip irrigation methods in Fakkeh Plain. Iran. J. Appl. Anim. Sci. 8, 3646–3653 (2008).ADS 

Google Scholar 
Kang, S. et al. Improving agricultural water productivity to ensure food security in China under changing environment: From research to practice. Agric. Water Manag. 179, 5–17 (2017).Article 

Google Scholar 
Yang, D. et al. Effect of drip irrigation on wheat evapotranspiration, soil evaporation and transpiration in Northwest China. Agric. Water Manag. 232, 106001 (2020).Article 

Google Scholar 
Xu, S. T., Zhang, L., Neil, B. & McLaughlin, Mi. Effect of synthetic and natural water absorbing soil amendment soil physical properties under potato production in a semi-arid region. Soil Till. Res. 148, 31–39 (2015).Article 

Google Scholar 
Roper, M. M., Ward, P. R., Keulen, A. F. & Hill, J. R. Under no-tillage and stubble retention, soil water content and crop growth are poorly related to soil water repellency. Soil Till. Res. 126, 143–150 (2013).Article 

Google Scholar 
Zhao, H. et al. Ridge-furrow with full plasticfilm mulching improves water use efficiency and tuber yields of potato in a se miarid rainfed ecosystem. Field Crop Research. 161, 137–148 (2014).Article 

Google Scholar 
Li, J. et al. Effects of micro-sprinkling with different irrigation amount on grain yield and water use efficiency of winter wheat in the North China Plain. Agric. Water Manag. 224, 105736 (2019).Article 

Google Scholar 
Chouhan, S. S., Awasthi, M. K. & Nema, R. K. Studies on water productivity and yields responses of wheat based on drip irrigation systems in clay loam soil. Indian J. Sci. Technol. 8, 650 (2015).Article 

Google Scholar 
Liao, L., Zhang, L. & Bengtsson, L. Soil moisture variation and water consumption of spring wheat and their effects on crop yield under drip irrigation. Irrigat. Drainag. Syst. 22, 253–270 (2008).Article 

Google Scholar 
Jha, S. K. et al. Response of growth, yield and water use efficiency of winter wheat to different irrigation methods and scheduling in North China Plain. Agric. Water Manag. 217, 292–302 (2019).Article 

Google Scholar 
Yan, Z., Fengxin, W., Qi, Z., Kaijing, Y. & Youliang, Z. Effect of drip tape distance and irrigation amount on spring wheat yield and water use efficiency. Chin. Agric. Sci. Bull. 32, 194–199 (2016).
Google Scholar 
Chen, R. et al. Lateral spacing in drip-irrigated wheat: the effects on soil moisture, yield, and water use efficiency. Field Crop Res. 179, 52–62 (2015).Article 

Google Scholar 
Shock, C. C., Feibert, E.B.G., & Saunders, L. D. Water management for drip-irrigated spring wheat. Annual Rep. Med. Chem.. 2007 (2005).Bhardwaj, A. K., Shainberg, I., Goldstein, D., Warrington, D. N. & Levy, G. J. Water retention and hydraulic conductivity of cross-linked polyacrylamides in sandy soils. Soil Sci. Soc. Am. J. 71, 406–412 (2007).ADS 
CAS 
Article 

Google Scholar 
Demitri, C., Scalera, F., Madaghiele, M., Sannino, A. & Maffezzoli, A. Potential of cellulose-based superabsorbent hydrogels as water reservoir in agriculture. Int. J. Polym. Sci. 2013, 1–6 (2013).Article 
CAS 

Google Scholar 
Islam, M. R. et al. Effectiveness of a water-saving super-absorbent polymer in soil water conservation for corn (Zea mays L.) based on eco- physiological parameters. J. Agric. Food Sci. 91, 1998–2005 (2011).CAS 
Article 

Google Scholar 
Nazarli, H., Zardashti, M. R., Darvishzadeh, R. & Najafi, S. The effect of water stress and polymer on water use efficiency, yield, and several morphological traits of sunflower under greenhouse condition. Notulae Scientia Biologicae. 2, 53–58 (2010).Article 

Google Scholar 
Huettermann, A., Orikiriza, L. J. & Agaba, H. Application of superabsorbent polymers for improving the ecological chemistry of degraded or polluted lands. Clean: Soil, Air, Water 37, 517–526 (2009).CAS 

Google Scholar 
Jain, N. K., Meena, H. N. & Bhaduri, D. Improvement in productivity, water use efficiency, and soil nutrient dynamics of summer peanut (Arachis hypogaea L) through use of polythene mulch, hydrogel, and nutrient management. Commun. Soil Sci. Plant Anal. 48, 549–564 (2017).CAS 
Article 

Google Scholar 
Shekari, F., Javanmard, A. & Abbasi, A. Effects of super absorbent polymer application on yield and yield components of rapeseed. Notulae Scientia Biologicae. 7, 361–366 (2015).Article 

Google Scholar 
Wang, L. et al. Drip irrigation mode and water-retaining agent on growth regulation and water-saving effect of small coffee. Chin. J. Drainag. Irrigat. Mech. Eng. 33, 796–801 (2015).
Google Scholar 
Liu, P. et al. Effects of soil treatments on soil moisture and soybean yield under the condition of underground drip irrigation. Water Saving Irrigat. 25–28 (2019).Li, R. et al. Effects of water-retaining agent on soil water, fertilizer and corn yield under drip irrigation. J. Drainag. Irrigat. Mech. Eng. 36, 1337–1344 (2018).
Google Scholar 
Ma, B. L., Biswas, D. K., Zhou, Q. P. & Ren, C. Z. Comparisons among cultivars of wheat, hulled and hulless oats: Effects of N fertilization on growth and yield. Can. J. Plant Sci. 92, 1213–1222 (2012).Article 

Google Scholar 
He, W. Effects of different irrigation methods on photosynthesis and soil biological characteristics of oat. Inner Mongolia: Hohhot, Inner Mongolia Agricultural University Master’s Thesis (2013).Wu, N. et al. Effects of water-retaining agent dosage on the yield and quality of naked oats under two irrigation methods. J. Crops 35, 1552–1557 (2009).CAS 

Google Scholar 
Gee, G.W., Bauder, J.W.,. Particle-size analysis. In: Klute, A. (Ed.), Methods of Soil Analysis, Part 1. Soil Science Society of America, South Segoe Road, Madison, WI 53711 USA. 383–409 (1986).Lu, R. Soil Agricultural Chemical Analysis Method (China Agricultural Science and Technology Press, 2000).
Google Scholar 
Wang, D. Water use efficiency and optimal supplemental irrigation in a high yield wheat field. Field Crop Res. 213, 213–220 (2017).Article 

Google Scholar 
Chen, Y. et al. Straw strips mulch on furrows improves water use efficiency and yield of potato in a rainfed semiarid area. Agric. Water Manag. 211, 142–151 (2019).Article 

Google Scholar 
Finn, D. et al. Effect of added nitrogen on plant litter decomposition depends on initial soil carbon and nitrogen stoichiometry. Soil Biol. Biochem. 91, 160–168 (2015).CAS 
Article 

Google Scholar 
Mo, F., Wang, J. Y., Xiong, Y. C., Nguluu, S. N. & Li, F. M. Ridge-furrow mulching system in semiarid Kenya: A promising solution to improve soil water availability and maize productivity. Eur. J. Agron. 80, 124–136 (2016).Article 

Google Scholar 
Luo, C. L. et al. Dual plastic film and straw mulching boosts wheat productivity and soil quality under the El Nino in semiarid Kenya. Sci. Total Environ. 738, 139808 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bengough, A. G. Water dynamics of the root zone: Rhizosphere biophysics and its control on soil hydrology. Vadose Zone Journal. 11, 1–6 (2012).Article 

Google Scholar 
Zobel, R. W. Plant Roots: Rowth, Activity and Interaction with Soils. Crop Sci. 46, 2699 (2006).Article 

Google Scholar 
Scholl, P. et al. Root induced changes of effective 1D hydraulic properties in a soil column. Plant Soil 381, 193–213 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Williams, S. M. & Weil, R. R. Crop cover root channels may alleviate soil compaction effects on soybean crop. Soil Sci. Soc. Am. J. 68, 1403–1409 (2010).Article 

Google Scholar 
Farrell, C., Ang, X. Q. & Rayner, J. P. Water-retention additives increase plant available water in green roof substrates. Ecol. Eng. 52, 112–118 (2013).Article 

Google Scholar 
Agaba, H. et al. Effects of hydrogel amendment to different soils on plant available water and survival of trees under drought conditions. Clean: Soil, Air, Water 38, 328–335 (2010).CAS 

Google Scholar 
Wu, L., Liu, M. Z. & Liang, R. Preparation and properties of a double-coated slow release NPK compound fertilizer with superabsorbent and water-retention. Biores. Technol. 99, 547–554 (2008).CAS 
Article 

Google Scholar 
Afshar, R. K. et al. Interactive effect of deficit irrigation and soil organic amendments on seed yield and flavonolignan production of milk thistle (Silybum marianum L. Gaertn.). Ind. Crops Prod. 58, 166–172 (2014).CAS 
Article 

Google Scholar 
Wang, L. Effects of different sowing dates and fertilizer rates on the growth and yield of oats in Yinshan hilly area. Hohhot, Inner Mongolia Agricultural University Master’s Thesis (2020).Liu, Y. G. et al. Influence of planting density on the yield of naked oats and its constituent factors. J. Wheat Crops 28, 140–143 (2008).
Google Scholar 
Jia, Z. F. Effects of sowing rate and row spacing on grain quality of naked oat. Seed. 32, 67–69 (2013).
Google Scholar 
Lascano, R. J. & Van Bavel, C. H. M. Stimulation and measurement of evaporation from bare soil. Soil Sci. Soc. Am. J. 50, 1127–1132 (1986).ADS 
Article 

Google Scholar 
Lv, P. et al. Effects of descending distance under wide sowing conditions on wheat yield and dry matter accumulation and transport. J. Wheat Crops 40, 1–6 (2020).
Google Scholar 
Sun, H. Y. et al. Effects of different row spacing on evapotranspiration and yield of winter wheat in wheat fields. Chin. J. Agric. Eng. 1, 22–26 (2006).
Google Scholar 
Li, G. X. et al. Effects of sowing row spacing on yield and water use efficiency of dryland wheat in different years. Agric. Technol. Equipm. 1, 22–26 (2012).ADS 

Google Scholar 
Chen, S. Y. et al. Effects of planting row spacing on soil evaporation and water use in winter wheat fields. Chin. J. Ecol. Agric. 14, 86–89 (2006).
Google Scholar 
Yang, Y. H. et al. Effects of water-retaining agent on soil moisture and utilization of winter wheat at different growth stages. Chin. J. Agric. Eng. 26, 19–26 (2010).
Google Scholar 
Yang, Y. H. et al. Effects of different moisture conservation tillage measures on water consumption characteristics and annual water use of wheat and maize. North China Agric. J. 32, 103–110 (2017).
Google Scholar 
Du, S. N. et al. Effects of Water and PAM Application Modes on Soil Moisture and Maize Growth. Chin. J. Agric. Eng. 24, 30–35 (2008).
Google Scholar 
Tian, L. et al. Effects of combined application of water-retaining agent and microbial fertilizer on dry matter accumulation, distribution, transport and yield of dry oat. J. Ecol. 39, 2996–3003 (2020).
Google Scholar  More

Pachauri, R. K. et al. Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change. (IPCC, 2014).Thackeray, S. J. et al. Phenological sensitivity to climate across taxa and trophic levels. Nature 535, 241–245 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Walther, G.-R. et al. Ecological responses to recent climate change. Nature 416, 389–395 (2002).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Root, T. L., MacMynowski, D. P., Mastrandrea, M. D. & Schneider, S. H. Human-modified temperatures induce species changes: Joint attribution. Proc. Natl. Acad. Sci. USA 102, 7465–7469 (2005).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Menzel, A. et al. European phenological response to climate change matches the warming pattern. Glob. Change Biol. 12, 1969–1976 (2006).ADS 
Article 

Google Scholar 
Rosenzweig, C. et al. Attributing physical and biological impacts to anthropogenic climate change. Nature 453, 353–357 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Myneni, R. B., Keeling, C. D., Tucker, C. J., Asrar, G. & Nemani, R. R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 386, 698–702 (1997).ADS 
CAS 
Article 

Google Scholar 
Walkovszky, A. Changes in phenology of the locust tree (Robinia pseudoacacia L.) in Hungary. Int. J. Biometeorol. 41, 155–160 (1998).ADS 
Article 

Google Scholar 
Crick, H. Q. P. & Sparks, T. H. Climate change related to egg-laying trends [8]. Nature 399, 423–424 (1999).ADS 
CAS 
Article 

Google Scholar 
Both, C., Van Asch, M., Bijlsma, R. G., Van Den Burg, A. B. & Visser, M. E. Climate change and unequal phenological changes across four trophic levels: Constraints or adaptations?. J. Anim. Ecol. 78, 73–83 (2009).PubMed 
Article 

Google Scholar 
Brown, J. L., Li, S. H. & Bhagabati, N. Long-term trend toward earlier breeding in an American bird: A response to global warming?. Proc. Natl. Acad. Sci. USA. 96, 5565–5569 (1999).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Miller-Rushing, A. J., Lloyd-Evans, T. L., Primack, R. B. & Satzinger, P. Bird migration times, climate change, and changing population sizes. Glob. Chang. Biol. 14, 1959–1972 (2008).ADS 
Article 

Google Scholar 
Hughes, L. Biological consequences of global warming: Is the signal already apparent?. Trends Ecol. Evol. 15, 56–61 (2000).CAS 
PubMed 
Article 

Google Scholar 
Kudo, G., Nishikawa, Y., Kasagi, T. & Kosuge, S. Does seed production of spring ephemerals decrease when spring comes early?. Ecol. Res. 19, 255–259 (2004).Article 

Google Scholar 
Doi, H., Gordo, O. & Katano, I. Heterogeneous intra-annual climatic changes drive different phenological responses at two trophic levels. Clim. Res. 36, 181–190 (2008).Article 

Google Scholar 
Both, C., Bouwhuis, S., Lessells, C. M. & Visser, M. E. Climate change and population declines in a long-distance migratory bird. Nature 441, 81–83 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Post, E. & Forchhammer, M. C. Climate change reduces reproductive success of an Arctic herbivore through trophic mismatch. Philos. Trans. R. Soc. B Biol. Sci. 363, 2369–2375 (2008).Article 

Google Scholar 
Thackeray, S. J. et al. Trophic level asynchrony in rates of phenological change for marine, freshwater and terrestrial environments. Glob. Chang. Biol. 16, 3304–3313 (2010).ADS 
Article 

Google Scholar 
Visser, M. E. & Both, C. Shifts in phenology due to global climate change: The need for a yardstick. Proc. R. Soc. B Biol. Sci. 272, 2561–2569 (2005).Article 

Google Scholar 
Hegland, S. J., Nielsen, A., Lázaro, A., Bjerknes, A. L. & Totland, Ø. How does climate warming affect plant-pollinator interactions?. Ecol. Lett. 12, 184–195 (2009).PubMed 
Article 

Google Scholar 
Brown, C. J. et al. Quantitative approaches in climate change ecology. Glob. Chang. Biol. 17, 3697–3713 (2011).ADS 
PubMed Central 
Article 

Google Scholar 
Parmesan, C., Duarte, C., Poloczanska, E., Richardson, A. J. & Singer, M. C. Overstretching attribution. Nat. Clim. Chang. 1, 2–4 (2011).ADS 
Article 

Google Scholar 
Damos, P. & Savopoulou-Soultani, M. Temperature-driven models for insect development and vital thermal requirements. Psyche (London) 2012, (2012).Osawa, T. et al. Climate-mediated population dynamics enhance distribution range expansion in a rice pest insect. Basic Appl. Ecol. 30, 41–51 (2018).Article 

Google Scholar 
Kiritani, K. Predicting impacts of global warming on population dynamics and distribution of arthropods in Japan. Popul. Ecol. 48, 5–12 (2006).Article 

Google Scholar 
Ozawa, A., Uchiyama, T. & Kasai, A. Estimating the day on which the numbers of adult tea spiny whiteflies (Aleurocanthus camelliae Kanmiya and Kasai) peaked and the number of generations produced in a year based on the effective cumulative temperature in tea fields. Annu. Rep. Kansai Plant Prot. Soc. 58, 57–64 (2016) (in Japanese).Article 

Google Scholar 
Ozawa, A., Saito, T. & Ikeda, F. Effect of host plant and temperature on reproduction of American Serpentine Leafminer, Liriomyza trifolii (Burgess). Japan. J. Appl. Entomol. Zool. 43, 41–48 (1999) (in Japanese).Article 

Google Scholar 
Baier, P., Pennerstorfer, J. & Schopf, A. PHENIPS—A comprehensive phenology model of Ips typographus (L.) (Col., Scolytinae) as a tool for hazard rating of bark beetle infestation. For. Ecol. Manage. 249, 171–186 (2007).Article 

Google Scholar 
Karuppaiah, V. & Sujayanad, G. K. Impact of climate change on population dynamics of insect pests. (2012).Lipper, L. et al. Climate-smart agriculture for food security. Nat. Clim. Chang. 4, 1068–1072 (2014).ADS 
Article 

Google Scholar 
Campbell, B. M. et al. Reducing risks to food security from climate change. Glob. Food Sec. 11, 34–43 (2016).Article 

Google Scholar 
Kiritani, K. The low development threshold temperature and the thermal constant in insects and mites in Japan (2nd edition). Bull. Nationai Instirute Agro-Environmental Sci. 31, 1–74 (2012) (in Japanese).
Google Scholar 
Nagasawa, A., Takahashi, A. & Higuchi, H. Host plant use for oviposition by Trigonotylus caelestialium (Hemiptera: Miridae) and Stenotus rubrovittatus (Hemiptera: Miridae). Appl. Entomol. Zool. 47, 331–339 (2012).Article 

Google Scholar 
Higuchi, H. Ecology and management of rice bugs causing pecky rice. Japan. J. Appl. Entomol. Zool. 54, 171–188 (2010) (in Japanese).Article 

Google Scholar 
Ohtomo, R. Occurrence and control of Stenotus rubrovittatus (Hemiptera: Miridae) in Touhoku area in Japan. Japan. J. Appl. Entomol. Zool. 57, 137–149 (2013) (in Japanese).Article 

Google Scholar 
Kiritani, K. The impact of global warming and land-use change on the pest status of rice and fruit bugs (Heteroptera) in Japan. Glob. Chang. Biol. 13, 1586–1595 (2007).ADS 
Article 

Google Scholar 
Nagasawa, A. & Higuchi, H. Suitability of poaceous plants for nymphal growth of the pecky rice bugs Trigonotylus caelestialium and Stenotus rubrovittatus (Hemiptera: Miridae) in Niigata, Japan. Appl. Entomol. Zool. 47, 421–427 (2012).Article 

Google Scholar 
Gordo, O. & Sanz, J. J. Temporal trends in phenology of the honey bee Apis mellifera (L.) and the small white Pieris rapae (L.) in the Iberian Peninsula (1952–2004). Ecol. Entomol. 31, 261–268 (2006).Article 

Google Scholar 
Sparks, T. H. & Yates, T. J. The effect of spring temperature on the appearance dates of British butterflies 1883–1993. Ecography (Cop.) 20, 368–374 (1997).Article 

Google Scholar 
Stefanescu, C., Penuelas, J. & Filella, I. Effects of climatic change on the phenology of butterflies in the northwest Mediterranean Basin. Glob. Chang. Biol. 9, 1494–1506 (2003).ADS 
Article 

Google Scholar 
Dell, D., Sparks, T. H. & Dennis, R. L. H. Climate change and the effect of increasing spring temperatures on emergence dates of the butterfly Apatura iris (Lepidoptera: Nymphalidae). Eur. J. Entomol. 102, 161–167 (2005).Article 

Google Scholar 
Tabuchi, K. et al. Rice bugs in the Tohoku Region: Their occurrence and damage from 2003 to 2013. Bull. Natl. Agric. Res. Cent. Tohoku Reg. 117, 63–115 (2015) (in Japanese).
Google Scholar 
Jolly, W. M. & Running, S. W. Effects of precipitation and soil water potential on drought deciduous phenology in the Kalahari. Glob. Chang. Biol. 10, 303–308 (2004).ADS 
Article 

Google Scholar 
Kriticos, D. J. et al. CliMond: Global high-resolution historical and future scenario climate surfaces for bioclimatic modelling. Methods Ecol. Evol. 3, 53–64 (2012).Article 

Google Scholar 
Takeda, A. & Shimizu, K. Characteristics of the damaged grain caused by the sorghum plant bug, Stenotus rubrovittatus (Hemiptera: Miridae) at different infection periods. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2009, 85–87 (2009) (in Japanese).
Google Scholar 
Ishimoto, M. Seasonal Prevalence of Occurrence of the Rice Leaf Bug, Trigonotylus caelestialium (Kirkaldy) (Heteroptera: Miridae) on Paddy Rice Plants. Jpn. J. Appl. Entomol. Zool. 48, 79–85 (2004) (in Japanese).
Katase, M., Shimizu, K., Siina, S., Hagiwara, K. & Iwai, H. Seasonal occurrence of rice bugs in the northern part of Chiba Prefecture. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2007, 99–104 (2007) (in Japanese).
Google Scholar 
Shintani, Y. Effect of seasonal variation in host–plant quality on the rice leaf bug, Trigonotylus caelestialium. Entomol. Exp. Appl. 133, 128–135 (2009).Article 

Google Scholar 
Takada, M. B., Yoshioka, A., Takagi, S., Iwabuchi, S. & Washitani, I. Multiple spatial scale factors affecting mirid bug abundance and damage level in organic rice paddies. Biol. Control 60, 169–174 (2012).Article 

Google Scholar 
Yoshioka, A., Takada, M. B. & Washitani, I. Landscape effects of a non-native grass facilitate source populations of a native generalist bug, Stenotus rubrovittatus, in a heterogeneous agricultural landscape. J. Insect Sci. 14, 1–14 (2014).Article 

Google Scholar 
Watanabe, T. & Higuchi, H. Recent occurrence and problem of rice bugs. Plant Prot. 60, 201–203 (2006) (in Japanese).
Google Scholar 
Seino, H. An estimation of distribution of meteorological elements using GIS and AMeDAS data. J. Agric. Meteorol. 48, 379–383 (1993) (in Japanese).Article 

Google Scholar 
Ishigooka, Y., Tsuneo, K., Nishimori, M., Hasegawa, T. & Ohno, H. Spatial characterization of recent hot summers in Japan with agro-climatic indices related to rice production. J. Agric. Meteorol. 67, 209–224 (2011).Article 

Google Scholar 
Yamasaki, K. et al. Intraspecific variations in life history traits of two pecky rice bug species from Japan: Mapping emergence dates and number of annual generations. Ecol. Evol. (2021).Sakagami, Y. & Korenaga, R. Triangle method—A simple method for the estimation of total effective temperature. Japan. J. Appl. Entomol. Zool. 25, 52–54 (1981) (in Japanese).Article 

Google Scholar 
Shigehisa, S. Seasonal changes in egg diapause induction and effects of photoperiod and temperature on egg diapause in the sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae). Japan. J. Appl. Entomol. Zool. 52, 229–232 (2008) (in Japanese).Article 

Google Scholar 
Ishimoto, M. Oviposition of sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae) on rice plants. Japan. J. Appl. Entomol. Zool. 55, 193–197 (2011) (in Japanese).Article 

Google Scholar 
Ogata, M. et al. Fecundity and longevity in adults of the sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae) under laboratory conditions. Japan. J. Entomol. 13, 129–132 (2010) (in Japanese).
Google Scholar 
Niiyama, T. & Iitomi, A. Optimal control timing of the rice leaf bug, Trigonotylus caelestialium (Heteroptera: Miridae). Annu. Rep. Soc. Plant Prot. North Japan. (in Japanese) (2003).Takeda, A., Shimizu, K., Shiina, S., Hagiwara, K. & Katase, M. Seasonal prevalence of Stenotus rubrovittatus (Hemiptera: Miridae) in a gramineous weed field and a rice field. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2008, 97–102 (2008) (in Japanese).
Google Scholar 
Niiyama, T. Studies on the ecology of Trigonotylus caelestialium and establishment of pesticide control techniques. Annu. Rep. Akita Prefect. Agric. Inst. 49, 147–180 (2009) (in Japanese).
Google Scholar 
Sato, M. Latest infestation period of Miridae causing pecky rice damage in Aomori Prefecture. Annu. Rep. Soc. Plant Prot. North Japan 2014, 129–134 (2014) (in Japanese).
Google Scholar 
Johnson, J. B. & Omland, K. S. Model selection in ecology and evolution. Trends Ecol. Evol. 19, 101–108 (2004).PubMed 
Article 

Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest package: tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).Article 

Google Scholar  More

Model descriptionThe system dynamics model divided an urban dog population into the following subpopulations: (i) free-roaming dogs (both owned and unowned free-roaming, i.e. unrestricted dogs found on streets), (ii) shelter dogs (unowned restricted dogs living in shelters), and (iii) owned dogs (owned home-dwelling restricted dogs) (Fig. 1). The subpopulations change in size by individuals flowing between the different subpopulations or from flows extrinsically modelled (i.e. flows from subpopulations not included in the systems model; the acquisition of dogs from breeders and friends to the owned dog population, and the immigration/emigration of dogs from other neighbourhoods).Ordinary differential equations were used to describe the dog population dynamics. The models were written in R version 3.6.128, and numerically solved using the Runge–Kutta fourth order integration scheme with a 0.01 step sizes using the package “deSolve”29,30. For the baseline model, Eqs. (1–3) were used to describe the rates of change of dog subpopulations in the absence of management.Baseline free-roaming dog population (S):$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S$$
(1)
In the baseline model, the free-roaming dog population (Eq. 1) increases through the free-roaming dog intrinsic growth rate (rs), and the abandonment and roaming of dogs from the owned dog population (α) and decreases through adoption to the owned dog population (δ). The intrinsic growth rate is the sum of the effects of births, deaths, immigration, and emigration, which are not modelled separately. In this model, the growth rate of the free-roaming dog population is reduced depending on the population size in relation to the carrying capacity, through the logistic equation (rreal = rmax(1 − S/Ks))31. In the baseline simulation, the free-roaming dog population rises over time, until it stabilises at an equilibrium size.Baseline shelter dog population (H):$$frac{dH}{dt}=gamma times O-beta times H- {mu }_{h}times H$$
(2)
The shelter dog population (Eq. 2) increases through relinquishment of owned dogs (γ) and decreases through the adoption of shelter dogs to the owned dog population (β), and through the shelter dog death rate (µh). There is no carrying capacity for the shelter dog population as we assumed that more housing would be created as the population increases. This allowed calculation of the resources required to house shelter dogs.Baseline owned dog population (O),$$frac{dO}{dt}={r}_{o} times Otimes (1-frac{O}{{K}_{o}})+beta times H+delta times S-alpha times O-gamma times O$$
(3)
The owned dog population (Eq. 3) increases through the owned dog growth rate (ro), adoption of shelter dogs (β), and adoption of free-roaming dogs (δ); and decreases through abandonment/roaming (α) and relinquishment (γ) of owned dogs to the shelter dog population. The growth rate of the owned dog population (ro) combines the birth, death, and acquisition rates from sources other than the street or shelters (e.g. breeders, friends) and was modelled as density dependent by the limit to growth logistic formula (1 − O/Ko).Parameter estimatesDetailed descriptions of parameter estimates are provided in the supplementary information. The simulated environment was based on the city of Lviv, Ukraine. This city has an area of 182 km2 and a human population size of 717,803. Parameters were estimated from literature, where possible, and converted to monthly rates (Table 1). Initial sizes of the dog populations were estimated for the baseline simulation, based on our previous research in Lviv32. The carrying capacity depends on the availability of resources (i.e. food, shelter, water, and human attitudes and behaviour33) and is challenging to estimate. We assumed the initial free-roaming and owned dog populations were at carrying capacity. Initial population sizes for simulations including interventions were determined by the equilibrium population sizes from the baseline simulation (i.e. the stable population size, the points at which the populations were no longer increasing/decreasing).Table 1 Parameter description, parameter value, and minimum and maximum values used in the sensitivity analysis for the systems model.Full size tableEstimating the rate at which owned dogs are abandoned is difficult, as abandonment rates are often reported per dog-owning lifetime32,34 and owners are likely to under-report abandonment of dogs. Similarly, it is challenging to estimate the rate that owned dogs move from restricted to unrestricted (i.e. free-roaming). For simplicity, we modelled a combined abandonment/roaming rate (α) of 0.003 per month, estimated based on our previous research in Lviv and from literature34,35,36. We derive the owned dog relinquishment rate (γ) from New et al.37. We estimated shelter (β) and free-roaming adoption rates (δ) from shelter data in Lviv. We set the maximum intrinsic growth rate for the free-roaming dogs (rs) at 0.03 per month, similar to that reported in literature17,19,38. We assumed that demand for dogs was met quickly through a supply of dogs from births, breeders and friends and set a higher growth rate for the owned dog population (ro) at 0.07 per month.We assumed shelters operated with a “no-kill” policy (i.e. dogs were not killed in shelters as part of population management) and included a shelter dog death rate (µh) of 0.008 per month to incorporate deaths due to euthanasia for behavioural problems or health problems, or natural mortality. We modelled neutered free-roaming dog death rate (µn) explicitly for the CNR intervention at a minimum death rate of 0.02 per month38,39,40,41.InterventionsSix intervention scenarios were modelled (Table 2): sheltering; culling; CNR; responsible ownership; combined CNR and responsible ownership; and combined CNR and sheltering, representing interventions feasibly applied and often reported27. Table 2 outlines the equations describing each intervention. To simulate a sheltering intervention, a proportion of the free-roaming dog population was removed and added to the shelter dog population at sheltering rate (σ). In culling interventions, a proportion of the free-roaming dog population was removed through culling (χ).Table 2 Description of intervention parameters and coverages for simulations applied at continuous and annual periodicities.Full size tableFree-roaming dog population with sheltering intervention:$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S-sigma times S$$
(4)
Shelter dog population with sheltering intervention:$$frac{dH}{dt}=gamma times O-beta times H- {mu }_{h}times H+sigma times S$$
(5)
Free-roaming dog population with a culling intervention:$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S-chi times S$$
(6)
To simulate a CNR intervention, an additional subpopulation was added to the system (Eq. 7): (iv) the neutered free-roaming dog population (N; neutered, free-roaming). In this simulation, a proportion of the intact (I) free-roaming dog population was removed and added to the neutered free-roaming dog population. A neutering rate (φ) was added to the differential equations describing the intact free-roaming and the neutered free-roaming dog populations. Neutering was assumed to be lifelong (e.g. gonadectomy); a neutered free-roaming dog could not re-enter the intact free-roaming dog subpopulation. Neutered free-roaming dogs were removed from the population through the density dependent neutered dog death rate (µn); death rate increased when the population was closer to the carrying capacity. The death rate was a non-linear function of population size and carrying capacity modelled using a table lookup function (Fig. S1). Neutered free-roaming dogs were also removed through adoption to the owned dog population, and we assumed that adoption rates did not vary between neutered and intact free-roaming dogs.Neutered free-roaming dog population:$$frac{dN}{dt}=varphi times I-{mu }_{n}times N-delta times N$$
(7)
Intact free-roaming dog population with neutering intervention.$$frac{dI}{dt}={r}_{s}times Itimes left(1-frac{(I+N)}{{K}_{s}}right)+alpha times O-delta times I-varphi times I$$
(8)
To simulate a responsible ownership intervention, the baseline model was applied with decreased rate of abandonment/roaming (α) and increased rate of shelter adoption (β). To simulate combined CNR and responsible ownership, a proportion of the intact free-roaming dog population was removed through the neutering rate (φ), abandonments/roaming decreased (α) and shelter adoptions increased (β). In combined CNR and sheltering interventions, a proportion of the intact free-roaming dog population (I) was removed through neutering (φ) and added to the neutered free-roaming dog population (N), and a proportion was removed through sheltering (σ) and added to the shelter dog population (H).Intact free-roaming dog population with combined CNR and sheltering interventions:$$frac{dI}{dt}={r}_{s}times Stimes left(1-frac{(I+N)}{{K}_{s}}right)+alpha times O-delta times I-varphi times I- sigma times I$$
(9)
Intervention length, periodicity, and coverageAll simulations were run for 70 years to allow populations to reach equilibrium. It is important to note that this is a theoretical model; running the simulations for 70 years allows us to compare the interventions, but does not accurately predict the size of the dog subpopulations over this long time period. Interventions were applied for two lengths of time: (i) the full 70-year duration of the simulation; and (ii) a five-year period followed by no further intervention, to simulate a single period of investment in population management. In each of these simulations, we modelled the interventions as (i) continuous (i.e. a constant rate of e.g. neutering) and (ii) annual (i.e. intervention applied once per year). Interventions were run at low, medium, and high coverages (Table 2). As the processes are not equivalent, we apply different percentages for the intervention coverage (culling/neutering/sheltering) and the percent increase/decrease in parameter rates for the responsible ownership intervention. Intervention coverage refers to the proportion of dogs that are culled/neutered/sheltered per year (i.e. 20%, 40% and 70% annually) and, for responsible ownership interventions, the decrease in abandonment/roaming rate and increase in the adoption rate of shelter dogs (30%, 60% and 90% increase/decrease from baseline values). To model a low (20%), medium (40%) and high (70%) proportion of free-roaming dogs caught, but where half of the dogs were sheltered and half were neutered-and-returned, combined CNR and sheltering interventions were simulated at half-coverage (e.g. intervention rate of 0.7 was simulated by 0.35 neutered and 0.35 sheltered). For continuous interventions, sheltering (σ), culling (χ), and CNR (φ) were applied continuously during the length of the intervention. For annual interventions, σ, χ, and φ were applied to the ordinary differential equations using a forcing function applied at 12-month intervals. In simulations that included responsible ownership interventions, the decrease in owned dog abandonment/roaming (α) and the increase in shelter adoption (β) was assumed instantaneous and continuous (i.e. rates did not change throughout the intervention).Model outputsThe primary outcome of interest was the impact of interventions on free-roaming dog population size. For interventions applied for the duration of the simulation, we calculated: (i) equilibrium population size for each population; (ii) percent decrease in free-roaming dog population; (iii) costs of intervention in terms of staff-time; and (iv) an overall welfare score. For interventions applied for a five-year period, we also calculated: (v) minimum free-roaming dog population size and percent reduction from initial population size; and (vi) the length of time between the end of the intervention and time-point at which the free-roaming dog population reached above 20,000 dogs (the assumed initial free-roaming dog population size of Lviv, based on our previous research32, see Supplementary Information for detail).The costs of population management interventions vary by country (e.g. staff salaries vary between countries) and by the method of application (e.g. method of culling, or resources provided in a shelter). To enable a comparison of the resources required for each intervention, the staff time (staff working-months) required to achieve the intervention coverage was calculated. While this does not incorporate the full costs of an intervention, as equipment (e.g. surgical equipment), advertising campaigns, travel costs for the animal care team, and facilities (e.g. clinic or shelter costs) are not included, it can be used as a proxy for intervention cost. Using data provided from VIER PFOTEN International, we estimated the average number of staff required to catch and neuter the free-roaming dog population and to house the shelter dog population in each intervention, using this data as a proxy for catching and sheltering/culling. The number of dogs that can be cared for per shelter staff varies by shelter. To account for this, we estimated two staff-to-dog ratios (low and high). Table 3 describes the staff requirements for the different interventions.Table 3 Staff required for interventions and the number of dogs processed per staff per day.Full size tableUsing the projected population sizes, the staff time required for each staff type (e.g. number of veterinarian-months of work required) was calculated for each intervention. Relative salaries for the different staff types were estimated (Table 3). The relative salaries were used to calculate the cost of the interventions by:[Staff time required × relative salary ] × €20,000.Where €20,000 was the estimated annual salary of a European veterinarian, allowing relative staff-time costs to be compared between the different interventions. Average annual costs were reported.To provide overall welfare scores for each of the interventions, we apply the estimated welfare scores on a one to five scale, for each of the dog subpopulations, as determined by Hogasen et al. (2013)22. This scale is based on the Five Freedoms (freedom from hunger and thirst; freedom from discomfort; freedom from pain, injury, or disease; freedom to express normal behaviour; freedom from fear and distress42,43) and was calculated using expert opinions from 60 veterinarians in Italy22. The scores were weighted by the participants’ self-reported knowledge of different dog subpopulations, which resulted in the following scores: 2.8 for shelter dogs (WH); 3.5 for owned dogs (WO); 3.1 for neutered free-roaming dogs (WN); and 2.3 for intact free-roaming dogs (WI)22.Using these estimated welfare scores, we calculated an average welfare score for the total dog population based on the model’s projected population sizes for each subpopulation (Eq. 10). For interventions running for the duration of the simulation, the welfare score was calculated at the time point (t) when the population reached an equilibrium size. For interventions running for five years, the welfare score was calculated at the end of the five-year intervention. The percentage change in welfare scores from the baseline simulation were reported.$$Welfare score= frac{{H}_{t}times {W}_{H}+{O}_{t}times {W}_{O}+{N}_{t}times {W}_{N}+{I}_{t}times {W}_{I}}{{H}_{t}+{O}_{t}+{N}_{t}+{I}_{t}}$$
(10)
Model validation and sensitivity analysisA global sensitivity analysis was conducted on all parameters described in the baseline simulation and all interventions applied continuously, at high coverage, for the full duration of the simulation. A Latin square design algorithm was used in package “FME”44 to sample the parameters within their range of values (Table 1). For the global sensitivity analysis on interventions, all parameter values were varied, apart from the parameters involved in the intervention (e.g. culling, neutering, abandonment/roaming rates). The effects of altering individual parameters (local sensitivity analysis) on the population equilibrium was also examined for the baseline simulation using the Latin square design algorithm to sample each parameter, individually, within their range of values. Sensitivity analyses were run for 100 simulations over 50 years solved with 0.01 step sizes. More

For this initial compilation, we focused on gathering traits from field guides and species accounts rather than the primary research literature because each represents the culmination of a comprehensive effort to describe a regional flora/fauna by local experts25. Authors of these guides have already done the hard work of scouring the literature, corresponding with fellow naturalists, and compiling occurrence records to support range, phenology, and habitat associations26. We began by performing a comprehensive review of all the holdings in the Florida Museum of Natural History’s McGuire Center for Lepidoptera and Biodiversity library, at the University of Florida. This, and subsequent searches in online databases, allowed us to compile a list of references that currently has more than 800 relevant resources.We initially identified the categories of trait information available in each resource and its format to target volumes for trait extraction and processing. Given the unequal availability of resources among regions, we had the explicit goal of identifying a corpus that would maximize the number of extractable trait data from as many butterfly species as evenly across the globe as possible. This led to our choice of 117 volumes within several global regions (Fig. 2, Supplementary Material S1) and a focus on measurements (wingspan/forewing length), phenology (months of adult flight and total duration of flight in months) and voltinism (the number of adult flight periods per year), habitat affinities, and host plants as traits (Table 1, Supplementary Material S2).Table 1 The total number of species represented by each trait in LepTraits 1.0.Full size tableTo process these resources, we developed a protocol to scan each volume, extract verbatim natural language descriptions, provide quality control for extraction, and then resolve given taxonomic names to a standardized list27. This provided a database of trait information in which each “cell” included all text from a single resource relevant to one trait category of a single taxon. In order to “atomize” the raw text into standardized metrics or a controlled list of descriptive terms, we developed a methodology appropriate to each trait. This resulted in a more fine-grained dataset in which each “cell” included a single, standardized trait value. Since the values of these taxon-specific traits frequently differed among resources, we then calculated “consensus” traits for each species, for example, the average forewing length (Table 1). A graphical representation of this process with an example trait is illustrated in Fig. 1.Fig. 1A graphical illustration of the processing workflow used to compile, scan, digitize, extract, atomize, and compile species trait records from literature resources. (1) Literature resources were examined for potential trait data and compiled into a single library; (2) each literature resource was scanned into.pdf format so that text could be readily copy and pasted from species accounts; (3) each.pdf file was uploaded to an online database with associated metadata for each literature resource; (4) trait extractors utilized an online interface to extract verbatim, raw text from designated resources; (5) verbatim, raw text extracts were either automatically (via regular-expressions and keyword searches) or manually atomized to a controlled vocabulary; (6) species consensus traits were calculated by aggregating resource-level records by name-normalized taxonomy. Rulesets were used for consensus trait building and are detailed in the supplementary material. Both resource-level and species consensus traits are presented in the dataset.Full size imageResource compilation and ingestionText sources from the master list were digitized by multiple participating institutions. They scanned each page of the book and converted the images to editable text with Abbyy FineReader optical character recognition (OCR) software (abbyy.com). These PDFs with copy-and-pastable text were then uploaded to a secure, online database that included citation information about each resource. The geographic breadth covered by each resource was designated using the World Geographic Scheme (WGS)28; this information was used to assess geographic evenness of our trait compilation efforts. Resource metadata, including the WGS scheme, were kept with each resource in an online database where individuals could access scanned copies of the resource for trait extraction.Verbatim data extractionIndividual workers were assigned to a resource and instructed to copy verbatim trait information from the original source. They then pasted that text into the relevant data field in a standardized, electronic form on an online portal designed to facilitate extraction and processing. Most field guides and other book-length resources are organized within a taxonomic hierarchy to describe traits of a family with a contiguous block of text, for example, family, then genus, species, and finally subspecies within species. We call these text blocks describing a single taxon “accounts” (e.g., family account, species account), and we recorded data at the taxonomic resolution provided in the original source. These taxonomic ranks included family, subfamily, tribe, genus, species, and subspecies. When information for a taxon was encountered outside its own account, the “extractor” (project personnel trained to manually extract verbatim text) assigned to glean data from the book entered this text into a separate entry for the taxon. Trait information from figure captions and tables were also extracted from the resource. Graphical representations of phenology and voltinism were common, and these visual data were converted to text descriptions. Each resource was extracted in stages, and each stage was subjected to a quality assurance and control process (see Technical Validation). This process corrected mistakes and attempted to find unextracted data overlooked by the extractor. These problems were corrected before the extractor could proceed with further trait extraction from the resource and were also used for training purposes.AtomizationVerbatim text extracts were subjected to an “atomization” process in which raw text was standardized into disaggregated, readily computable data. This conversion into the final trait data format (numerical, categorical, etc.) was two-pronged and involved both manual editing and semi-automated atomization of verbatim text. Regular expressions were used for most semi-automated atomization, including extraction of wing measurements, which were converted into centimeters. Keyword searches were also performed in the semi-automated pipeline for phenology, voltinism, and oviposition traits. For example, “univoltine” or “uni*” was searched for across the voltinism raw text, along with other search terms. All semi-automated atomization outputs were subject to quality assurance and control detailed further in Technical Validation. Manual atomization tasks were performed by multiple team members for traits which presented higher complexity. For example, habitat affinities and host plant associations were atomized manually along with a quality control protocol based on predefined rule sets that are described further in the Supplementary Material S3.Normalization and consensus traitsTo provide consensus traits at the species (and sometimes genus) level, we standardized nomenclature through a process we called “name-normalization,” which harmonizes taxonomy across all of our resources29. This name-normalization procedure relied on a comprehensive catalog of valid names and synonyms27. Following taxonomic harmonization, we compiled consensus traits based on rule sets specified in the metadata of each trait. For example, species-level consensus of primary and secondary host plant families required that at least one-third of the records for a given taxon list a particular family of plants (when multiple records were available).Categorical traits such as voltinism list all known voltinism patterns for a species regardless of geographic context. To this end, it is important that users of these data are aware that not all traits may be applicable to their study region. For example, some species may be univoltine at higher latitudes or elevations, but bivoltine elsewhere. We therefore present both the resource-level records as well as the species consensus traits for use in analysis.For this initial synopsis of butterfly species traits, we extracted records from 117 literature/web-based resources, resulting in 75,103 individual trait extraction records across 12,448 unique species, out of the ca. 19,200 species described to date27. Figure 2 indicates the geographic regions covered by our 117 resources, mapped at the resolution level-two regions in the World Geographic Scheme28. A full list of resources can be found in the Supplemental Material S1 as a bibliography. Similarly, the geographic distribution of trait records is indicated in Fig. 3. Resource and consensus species trait records varied in number and in the scope of taxonomic coverage. Table 1 indicates the number of unique records and species level records for each trait. Table 2 indicates the number of species-level records by family. Measurement traits, including wingspan and forewing length, were the most comprehensive traits extracted from our resource set. This represents one of the largest trait datasets and the most comprehensive dataset for butterflies to date.Fig. 2Geographic breadth of our butterfly trait resources. Using a global map of level-two regions (World Geographic Scheme, Brummitt 2001), we have indicated the total number of resources available within each geographic area). Grey areas indicate that no resources were extracted from that region.Full size imageFig. 3Geographic breadth of our butterfly trait records. Using a global map of level-two regions(World Geographic Scheme, Brummitt 2001), we have indicated the total number of trait records from each geographic region). Grey areas indicate that trait records were not extracted from that region.Full size imageTable 2 The number of species represented within each family in LepTraits 1.0.Full size table More

Neufeld JD, Wagner M, Murrell JC. Who eats what, where and when? Isotope-labelling experiments are coming of age. ISME J. 2007;1:103–10.CAS 
PubMed 
Article 

Google Scholar 
Boschker HTS, Nold SC, Wellsbury P, Bos D, de Graaf W, Pel R, et al. Direct linking of microbial populations to specific biogeochemical processes by 13C-labelling of biomarkers. Nature. 1998;392:801–5Jehmlich N, Schmidt F, von Bergen M, Richnow H-H, Vogt C. Protein-based stable isotope probing (Protein-SIP) reveals active species within anoxic mixed cultures. ISME J. 2008;2:1122–33.CAS 
PubMed 
Article 

Google Scholar 
Radajewski S, Ineson P, Parekh NR, Colin Murrell J. Stable-isotope probing as a tool in microbial ecology. Nature. 2000;403:646–9.CAS 
PubMed 
Article 

Google Scholar 
Manefield M, Whiteley AS, Griffiths RI, Bailey MJ. RNA stable isotope probing, a novel means of linking microbial community function to phylogeny. Appl Environ Microbiol. 2002;68:5367–73.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Berry D, Mader E, Lee TK, Woebken D, Wang Y, Zhu D, et al. Tracking heavy water (D2O) incorporation for identifying and sorting active microbial cells. Proc Natl Acad Sci USA. 2015;112:E194–203.CAS 
PubMed 

Google Scholar 
Jehmlich N, Vogt C, Lünsmann V, Richnow HH, von Bergen M. Protein-SIP in environmental studies. Curr Opin Biotechnol. 2016;41:26–33.CAS 
PubMed 
Article 

Google Scholar 
Haichar, FEZ, Achouak W, Christen R, Heulin T, et al. Identification of cellulolytic bacteria in soil by stable isotope probing. Environ Microbiol. 2007;9:625–34Rangel-Castro JI, Ignacio Rangel-Castro J, Killham K, Ostle N, Nicol GW, Anderson IC, et al. Stable isotope probing analysis of the influence of liming on root exudate utilization by soil microorganisms. Environ Microbiol. 2005;7:828–38.CAS 
PubMed 
Article 

Google Scholar 
Wang Y, Song Y, Tao Y, Muhamadali H, Goodacre R, Zhou N-Y, et al. Reverse and multiple stable isotope probing to study bacterial metabolism and interactions at the single cell level. Anal Chem. 2016;88:9443–50.CAS 
PubMed 
Article 

Google Scholar 
Sharma K, Palatinszky M, Nikolov G, Berry D, Shank EA. Transparent soil microcosms for live-cell imaging and non-destructive stable isotope probing of soil microorganisms. Elife. 2020;9:e56275.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lee KS, Landry Z, Pereira FC, Wagner M, Berry D, Huang WE, et al. Raman microspectroscopy for microbiology. Nat. Rev. Methods Primers. 2021;1:80.CAS 
Article 

Google Scholar 
Hatzenpichler R, Krukenberg V, Spietz RL, Jay ZJ. Next-generation physiology approaches to study microbiome function at single cell level. Nat Rev Microbiol. 2020;18:241–56.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wagner M. Single-cell ecophysiology of microbes as revealed by raman microspectroscopy or secondary ion mass spectrometry imaging. Ann Rev Microbiol. 2009;63:411–29Lennon JT, Jones SE. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol. 2011;9:119–30.CAS 
PubMed 
Article 

Google Scholar 
Lewis K. Persister cells, dormancy and infectious disease. Nat Rev Microbiol. 2006;5:48–56.PubMed 
Article 
CAS 

Google Scholar 
Nielsen KM, Johnsen PJ, Bensasson D, Daffonchio D. Release and persistence of extracellular DNA in the environment. Environ Biosafety Res. 2007;6:37–53.CAS 
PubMed 
Article 

Google Scholar 
Blazewicz SJ, Barnard RL, Daly RA, Firestone MK. Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 2013;7:2061–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nocker A, Sossa-Fernandez P, Burr MD, Camper AK. Use of propidium monoazide for live/dead distinction in microbial ecology. Appl Environ Microbiol. 2007;73:5111–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tawakoli PN, Al-Ahmad A, Hoth-Hannig W, Hannig M, Hannig C. Comparison of different live/dead stainings for detection and quantification of adherent microorganisms in the initial oral biofilm. Clin Oral Investig. 2013;17:841–50.CAS 
PubMed 
Article 

Google Scholar 
Netuschil L, Auschill TM, Sculean A, Arweiler NB. Confusion over live/dead stainings for the detection of vital microorganisms in oral biofilms-which stain is suitable? BMC Oral Health. 2014;14:2.PubMed 
PubMed Central 
Article 

Google Scholar 
Hatzenpichler R, Connon SA, Goudeau D, Malmstrom RR, Woyke T, Orphan VJ. Visualizing in situ translational activity for identifying and sorting slow-growing archaeal−bacterial consortia. Proc Natl Acad Sci USA. 2016;113:E4069–78.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kuru E, Hughes HV, Brown PJ, Hall E, Tekkam S, Cava F, et al. In Situ probing of newly synthesized peptidoglycan in live bacteria with fluorescent D-amino acids. Angew Chem Int Ed Engl. 2012;51:12519–23.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kopf SH, McGlynn SE, Green-Saxena A, Guan Y, Newman DK, Orphan VJ. Heavy water and15N labelling with NanoSIMS analysis reveals growth rate-dependent metabolic heterogeneity in chemostats. Environ Microbiol. 2015;17:2542–56Kopf SH, Sessions AL, Cowley ES, Reyes C, Van Sambeek L, Hu Y, et al. Trace incorporation of heavy water reveals slow and heterogeneous pathogen growth rates in cystic fibrosis sputum. Proc Natl Acad Sci USA. 2016;113:E110–6.CAS 
PubMed 
Article 

Google Scholar 
Neubauer C, Kasi AS, Grahl N, Sessions AL, Kopf SH, Kato R, et al. Refining the Application of Microbial Lipids as Tracers of Staphylococcus aureus Growth Rates in Cystic Fibrosis Sputum. J Bacteriol. 2018;200:e00365–18.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Haider S, Wagner M, Schmid MC, Sixt BS, Christian JG, Häcker G, et al. Raman microspectroscopy reveals long-term extracellular activity of Chlamydiae. Mol Microbiol. 2010;77:687–700.CAS 
PubMed 
Article 

Google Scholar 
Kloehn J, Boughton BA, Saunders EC, O’Callaghan S, Binger KJ, McConville MJ. Identification of Metabolically Quiescent Leishmania mexicana Parasites in Peripheral and Cured Dermal Granulomas Using Stable Isotope Tracing Imaging Mass Spectrometry. mBio. 2021;12:e00129–21.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kong L, Setlow P, Li Y-Q. Direct analysis of water content and movement in single dormant bacterial spores using confocal Raman microspectroscopy and Raman imaging. Anal Chem. 2013;85:7094–101.CAS 
PubMed 
Article 

Google Scholar 
Knudsen SM, Cermak N, Delgado FF, Setlow B, Setlow P, Manalis SR. Water and small-molecule permeation of dormant Bacillus subtilis spores. J Bacteriol. 2016;198:168–77.CAS 
PubMed 
Article 

Google Scholar 
Chen D, Huang S-S, Li Y-Q. Real-time detection of kinetic germination and heterogeneity of single Bacillus spores by laser tweezers Raman spectroscopy. Anal Chem. 2006;78:6936–41.CAS 
PubMed 
Article 

Google Scholar 
Devictor V, Clavel J, Julliard R, Lavergne S, Mouillot D, Thuiller W, et al. Defining and measuring ecological specialization. J Appl Ecol. 2010;47:15–25.Article 

Google Scholar 
Pereira FC, Berry D. Microbial nutrient niches in the gut. Environ Microbiol. 2017;19:1366–78.PubMed 
PubMed Central 
Article 

Google Scholar 
Shakya M, Lo C-C, Chain PSG. Advances and challenges in metatranscriptomic analysis. Front Genet. 2019;10:904.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Berry D, Loy A. Stable-Isotope probing of human and animal microbiome function. Trends Microbiol. 2018;26:999–1007.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Terrado R, Pasulka AL, Lie AA-Y, Orphan VJ, Heidelberg KB, Caron DA. Autotrophic and heterotrophic acquisition of carbon and nitrogen by a mixotrophic chrysophyte established through stable isotope analysis. ISME J. 2017;11:2022–34.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dekas AE, Parada AE, Mayali X, Fuhrman JA, Wollard J, Weber PK, et al. Characterizing Chemoautotrophy and Heterotrophy in Marine Archaea and Bacteria With Single-Cell Multi-isotope NanoSIP. Front Microbiol. 2019;10:2682.PubMed 
PubMed Central 
Article 

Google Scholar 
Wegener G, Bausch M, Holler T, Thang NM, Mollar XP, Kellermann MY, et al. Assessing sub-seafloor microbial activity by combined stable isotope probing with deuterated water and 13C-bicarbonate. Environ Microbiol. 2019;14:1517–27Jing X, Gou H, Gong Y, Su X, Xu L, Ji Y, et al. Raman-activated cell sorting and metagenomic sequencing revealing carbon-fixing bacteria in the ocean. Environ Microbiol. 2018;20:2241–55.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Xu J, Zhu D, Ibrahim AD, Allen CCR, Gibson CM, Fowler PW, et al. Raman deuterium isotope probing reveals microbial metabolism at the single-cell level. Anal Chem. 2017;89:13305–12.CAS 
PubMed 
Article 

Google Scholar 
Zhang M, Hong W, Abutaleb NS, Li J, Dong P-T, Zong C, et al. Rapid determination of antimicrobial susceptibility by stimulated Raman scattering imaging of D2O metabolic incorporation in a single bacterium. Adv Chem Microsc Life Sci Transl Med. 2021.Lima C, Muhamadali H, Xu Y, Kansiz M, Goodacre R. Imaging Isotopically Labeled Bacteria at the Single-Cell Level Using High-Resolution Optical Infrared Photothermal Spectroscopy. Anal Chem. 2021;93:3082–8.CAS 
PubMed 
Article 

Google Scholar 
Ackermann M. A functional perspective on phenotypic heterogeneity in microorganisms. Nat Rev Microbiol. 2015;13:497–508.CAS 
PubMed 
Article 

Google Scholar 
Balaban NQ, Merrin J, Chait R, Kowalik L, Leibler S. Bacterial persistence as a phenotypic switch. Science. 2004;305:1622–5.CAS 
PubMed 
Article 

Google Scholar 
Maamar H, Raj A, Dubnau D. Noise in gene expression determines cell fate in Bacillus subtilis. Science. 2007;317:526–9.CAS 
PubMed 
Article 

Google Scholar 
Emonet T, Cluzel P. Relationship between cellular response and behavioral variability in bacterial chemotaxis. Proc Natl Acad Sci USA. 2008;105:3304–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ozbudak EM, Thattai M, Lim HN, Shraiman BI, Van Oudenaarden A. Multistability in the lactose utilization network of Escherichia coli. Nature. 2004;427:737–40.CAS 
PubMed 
Article 

Google Scholar 
Kiviet DJ, Nghe P, Walker N, Boulineau S, Sunderlikova V, Tans SJ. Stochasticity of metabolism and growth at the single-cell level. Nature. 2014;514:376–9.CAS 
PubMed 
Article 

Google Scholar 
Kotte O, Volkmer B, Radzikowski JL, Heinemann M. Phenotypic bistability in Escherichia coli’s central carbon metabolism. Mol Syst Biol. 2014;10:736.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
New AM, Cerulus B, Govers SK, Perez-Samper G, Zhu B, Boogmans S, et al. Different levels of catabolite repression optimize growth in stable and variable environments. PLoS Biol. 2014;12:e1001764.PubMed 
PubMed Central 
Article 

Google Scholar 
Solopova A, van Gestel J, Weissing FJ, Bachmann H, Teusink B, Kok J, et al. Bet-hedging during bacterial diauxic shift. Proc Natl Acad Sci USA. 2014;111:7427–32.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Schreiber F, Littmann S, Lavik G, Escrig S, Meibom A, Kuypers MMM, et al. Phenotypic heterogeneity driven by nutrient limitation promotes growth in fluctuating environments. Nat Microbiol. 2016;1:16055.CAS 
PubMed 
Article 

Google Scholar 
Nikolic N, Schreiber F, Dal Co A, Kiviet DJ, Bergmiller T, Littmann S, et al. Cell-to-cell variation and specialization in sugar metabolism in clonal bacterial populations. PLoS Genet. 2017;13:e1007122.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Takhaveev V, Heinemann M. Metabolic heterogeneity in clonal microbial populations. Curr Opin Microbiol. 2018;45:30–8.CAS 
PubMed 
Article 

Google Scholar 
Altschuler SJ, Wu LF. Cellular heterogeneity: do differences make a difference? Cell. 2010;141:559–63.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Beaumont HJE, Gallie J, Kost C, Ferguson GC, Rainey PB. Experimental evolution of bet hedging. Nature. 2009;462:90–3.CAS 
PubMed 
Article 

Google Scholar 
Calabrese F, Voloshynovska I, Musat F, Thullner M, Schlömann M, Richnow HH, et al. Quantitation and comparison of phenotypic heterogeneity among single cells of monoclonal microbial populations. Front Microbiol. 2019;10:2814.PubMed 
PubMed Central 
Article 

Google Scholar 
Zimmermann M, Escrig S, Hübschmann T, Kirf MK, Brand A, Inglis RF, et al. Phenotypic heterogeneity in metabolic traits among single cells of a rare bacterial species in its natural environment quantified with a combination of flow cell sorting and NanoSIMS. Front Microbiol. 2015;6:243.PubMed 
PubMed Central 
Article 

Google Scholar 
Zimmermann M, Escrig S, Lavik G, Kuypers MMM, Meibom A, Ackermann M, et al. Substrate and electron donor limitation induce phenotypic heterogeneity in different metabolic activities in a green sulphur bacterium. Environ Microbiol Rep. 2018;10:179–83.CAS 
PubMed 
Article 

Google Scholar 
Sheik AR, Muller EE, Audinot J-N, Lebrun LA, Grysan P, Guignard C, et al. In situ phenotypic heterogeneity among single cells of the filamentous bacterium Candidatus Microthrix parvicella. ISME J. 2016;10:1274–9.CAS 
PubMed 
Article 

Google Scholar 
Seyedsayamdost MR, Case RJ, Kolter R, Clardy J. The Jekyll-and-Hyde chemistry of Phaeobacter gallaeciensis. Nat Chem. 2011;3:331–5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ferrier-Pagès C, Leal MC. Stable isotopes as tracers of trophic interactions in marine mutualistic symbioses. Ecol Evol. 2019;9:723–40.PubMed 
Article 

Google Scholar 
Pasulka AL, Thamatrakoln K, Kopf SH, Guan Y, Poulos B, Moradian A, et al. Interrogating marine virus-host interactions and elemental transfer with BONCAT and nanoSIMS-based methods. Environ Microbiol. 2018;20:671–92.CAS 
PubMed 
Article 

Google Scholar 
Kopp C, Domart-Coulon I, Escrig S, Humbel BM, Hignette M, Meibom A. Subcellular investigation of photosynthesis-driven carbon assimilation in the symbiotic reef coral Pocillopora damicornis. mBio. 2015;6:e02299–14.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rädecker N, Pogoreutz C, Gegner HM, Cárdenas A, Roth F, Bougoure J, et al. Heat stress destabilizes symbiotic nutrient cycling in corals. Proc Natl Acad Sci U S A. 2021;118:e2022653118.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Krueger T, Bodin J, Horwitz N, Loussert-Fonta C, Sakr A, Escrig S, et al. Temperature and feeding induce tissue level changes in autotrophic and heterotrophic nutrient allocation in the coral symbiosis – a NanoSIMS study. Sci Rep. 2018;8:12710.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Gibbin E, Gavish A, Krueger T, Kramarsky-Winter E, Shapiro O, Guiet R, et al. Vibrio coralliilyticus infection triggers a behavioural response and perturbs nutritional exchange and tissue integrity in a symbiotic coral. ISME J. 2019;13:989–1003.CAS 
PubMed 
Article 

Google Scholar 
Rix L, Ribes M, Coma R, Jahn MT, de Goeij JM, van Oevelen D, et al. Heterotrophy in the earliest gut: a single-cell view of heterotrophic carbon and nitrogen assimilation in sponge-microbe symbioses. ISME J. 2020;14:2554–67.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Thomas T, Moitinho-Silva L, Lurgi M, Björk JR, Easson C, Astudillo-García C, et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun. 2016;7:11870.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mills MM, Turk-Kubo KA, van Dijken GL, Henke BA, Harding K, Wilson ST, et al. Unusual marine cyanobacteria/haptophyte symbiosis relies on N2 fixation even in N-rich environments. ISME J. 2020;14:2395–406.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Turk-Kubo KA, Mills MM, Arrigo KR, van Dijken G, Henke BA, Stewart B, et al. UCYN-A/haptophyte symbioses dominate N2 fixation in the Southern California Current System. ISME Commun. 2021;1:1–13.Article 

Google Scholar 
Moore CM, Mills MM, Arrigo KR, Berman-Frank I, Bopp L, Boyd PW, et al. Processes and patterns of oceanic nutrient limitation. Nat Geosci. 2013;6:701–10.CAS 
Article 

Google Scholar 
Scheller S, Yu H, Chadwick GL, McGlynn SE, Orphan VJ. Artificial electron acceptors decouple archaeal methane oxidation from sulfate reduction. Science. 2016;351:703–7.CAS 
PubMed 
Article 

Google Scholar 
Pereira FC, Wasmund K, Cobankovic I, Jehmlich N, Herbold CW, Lee KS, et al. Rational design of a microbial consortium of mucosal sugar utilizers reduces Clostridiodes difficile colonization. Nat Commun. 2020;11:5104.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mooshammer M, Kitzinger K, Schintlmeister A, Ahmerkamp S, Nielsen JL, Nielsen PH, et al. Flow-through stable isotope probing (Flow-SIP) minimizes cross-feeding in complex microbial communities. ISME J. 2021;15:348–53.CAS 
PubMed 
Article 

Google Scholar 
Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL. Solutions to the public goods dilemma in bacterial biofilms. Curr Biol. 2014;24:50–5.CAS 
PubMed 
Article 

Google Scholar 
Słomka J, Alcolombri U, Secchi E, Stocker R, Fernandez VI. Encounter rates between bacteria and small sinking particles. New J Phys. 2020;22:043016.Article 

Google Scholar 
Alcolombri U, Peaudecerf FJ, Fernandez VI, Behrendt L, Lee KS, Stocker R. Sinking enhances the degradation of organic particles by marine bacteria. Nat Geosci. 2021;14:775–80.CAS 
Article 

Google Scholar 
University of Massachusetts Amherst Massachusetts Lynn Margulis, Margulis L, Fester R. Symbiosis as a source of evolutionary innovation: speciation and morphogenesis. MIT Press; 1991. 454 p.Legin AA, Schintlmeister A, Sommerfeld NS, Eckhard M, Theiner S, Reipert S, et al. Nano-scale imaging of dual stable isotope labeled oxaliplatin in human colon cancer cells reveals the nucleolus as a putative node for therapeutic effect. Nanoscale Adv. 2021;3:249–62.CAS 
Article 

Google Scholar 
Schaible GA, et al. Correlative SIP-FISH-Raman-SEM-NanoSIMS links identity, morphology, biochemistry, and physiology of environmental microbes. ISME COMMUN. 2022;2:52.Article 

Google Scholar 
Yu G-H, Chi Z-L, Kappler A, Sun F-S, Liu C-Q, Teng HH, et al. Fungal nanophase particles catalyze iron transformation for oxidative stress removal and iron acquisition. Curr Biol. 2020;30:2943–50.e4.CAS 
PubMed 
Article 

Google Scholar 
Subirana MA, Riemschneider S, Hause G, Dobritzsch D, Schaumlöffel D, Herzberg M. High spatial resolution imaging of subcellular macro and trace element distribution during phagocytosis. Metallomics. 2022;14:mfac011.PubMed 
Article 

Google Scholar 
Bonnin EA, Fornasiero EF, Lange F, Turck CW, Rizzoli SO. NanoSIMS observations of mouse retinal cells reveal strict metabolic controls on nitrogen turnover. BMC Mol Cell Biol. 2021;22:5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Jo MC, Liu W, Gu L, Dang W, Qin L. High-throughput analysis of yeast replicative aging using a microfluidic system. Proc Natl Acad Sci U S A. 2015;112:9364–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Anggraini D, Ota N, Shen Y, Tang T, Tanaka Y, Hosokawa Y, et al. Recent advances in microfluidic devices for single-cell cultivation: methods and applications. Lab Chip. 2022;22:1438–68.CAS 
PubMed 
Article 

Google Scholar 
Eriksen R, Daria V, Gluckstad J. Fully dynamic multiple-beam optical tweezers. Opt Express. 2002;10:597–602.PubMed 
Article 

Google Scholar 
Dai X, Fu W, Chi H, Mesias VSD, Zhu H, Leung CW, et al. Optical tweezers-controlled hotspot for sensitive and reproducible surface-enhanced Raman spectroscopy characterization of native protein structures. Nat Commun. 2021;12:1292.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Collins DJ, Morahan B, Garcia-Bustos J, Doerig C, Plebanski M, Neild A. Two-dimensional single-cell patterning with one cell per well driven by surface acoustic waves. Nat Commun. 2015;6:8686.CAS 
PubMed 
Article 

Google Scholar 
Hu F, Shi L, Min W. Biological imaging of chemical bonds by stimulated Raman scattering microscopy. Nat Methods. 2019;16:830–42.CAS 
PubMed 
Article 

Google Scholar 
Ge X, Pereira FC, Mitteregger M, Berry D, Zhang M, Hausmann B, et al. SRS-FISH: A high-throughput platform linking microbiome metabolism to identity at the single-cell level. Proc Natl Acad Sci U S A. 2022;119:e2203519119.PubMed 
PubMed Central 
Article 

Google Scholar 
Vandergrift GW, Kew W, Lukowski JK, Bhattacharjee A, Liyu AV, Shank EA, et al. Imaging and direct sampling capabilities of nanospray desorption electrospray ionization with absorption-mode 21 Tesla Fourier transform ion cyclotron resonance mass spectrometry. Anal Chem. 2022;94:3629–36.CAS 
PubMed 
Article 

Google Scholar 
Harrison JP, Berry D. Vibrational spectroscopy for imaging single microbial cells in complex biological samples. Front Microbiol. 2017;8:675.PubMed 
PubMed Central 
Article 

Google Scholar 
Mayali X. NanoSIMS: microscale quantification of biogeochemical activity with large-scale impacts. Ann Rev Mar Sci. 2020;12:449–67.PubMed 
Article 

Google Scholar 
Alexandrov T. Spatial metabolomics and imaging mass spectrometry in the age of artificial intelligence. Annu Rev Biomed Data Sci. 2020;3:61–87.PubMed 
PubMed Central 
Article 

Google Scholar 
Boschker HTS, Middelburg JJ. Stable isotopes and biomarkers in microbial ecology. FEMS Microbiol Ecol. 2002;40:85–95.CAS 
PubMed 
Article 

Google Scholar 
Mayali X, Weber PK, Nuccio E, Lietard J, Somoza M, Blazewicz SJ, et al. Chip-SIP: Stable Isotope Probing analyzed with rRNA-targeted microarrays and nanoSIMS. Methods Mol Biol. 2019;2046:71–87.PubMed 
Article 

Google Scholar 
Chokkathukalam A, Kim D-H, Barrett MP, Breitling R, Creek DJ. Stable isotope-labeling studies in metabolomics: new insights into structure and dynamics of metabolic networks. Bioanalysis. 2014;6:511–24.CAS 
PubMed 
Article 

Google Scholar 
Hiller K, Metallo CM, Kelleher JK, Stephanopoulos G. Nontargeted elucidation of metabolic pathways using stable-isotope tracers and mass spectrometry. Anal Chem. 2010;82:6621–8.CAS 
PubMed 
Article 

Google Scholar 
Rusconi R, Garren M, Stocker R. Microfluidics expanding the frontiers of microbial ecology. Annu Rev Biophys. 2014;43:65–91.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lee KS, Pereira FC, Palatinszky M, Behrendt L, Alcolombri U, Berry D, et al. Optofluidic Raman-activated cell sorting for targeted genome retrieval or cultivation of microbial cells with specific functions. Nat Protoc. 2021;16:634–76.CAS 
PubMed 
Article 

Google Scholar 
Wagner M, Haider S. New trends in fluorescence in situ hybridization for identification and functional analyses of microbes. Curr Opin Biotechnol. 2012;23:96–102.CAS 
PubMed 
Article 

Google Scholar  More

USDA, N. Census of Horticultural Specialties (USDA, 2014).
Google Scholar 
USDA, N. Census of Horticultural Specialties (USDA, 2019).
Google Scholar 
Soliman, A. S. & Shanan, N. T. The role of natural exogenous foliar applications in alleviating salinity stress in Lagerstroemia indica L. seedlings. J. Appl. Hortic. 19, 35–45 (2017).Article 

Google Scholar 
Chappell, M. R., Braman, S. K., Williams-Woodward, J. & Knox, G. J. J. o. E. H. Optimizing plant health and pest management of Lagerstroemia spp. in commercial production and landscape situations in the southeastern United States: A review. 30, 161–172 (2012).Gu, M., Merchant, M., Robbins, J. & Hopkins, J. Crape Myrtle Bark Scale: A New Exotic Pest. Texas A&M AgriLife Ext. Service. EHT 49 (2014).Kondo, T., Gullan, P. J. & Williams, D. J. Coccidology. The study of scale insects (Hemiptera: Sternorrhyncha: Coccoidea). Ciencia y Tecnología Agropecuaria 9, 55–61 (2008).Article 

Google Scholar 
Jiang, N. & Xu, H. Observertion on Eriococcus lagerostroemiae Kuwana. J. Anhui Agric. Coll. 25, 142–144 (1998).
Google Scholar 
He, D., Cheng, J., Zhao, H. & Chen, S. Biological characteristic and control efficacy of Eriococcus lagerstroemiae. Chin. Bull. Entomol. 45, 812–814 (2008).
Google Scholar 
Harcourt, D. The development and use of life tables in the study of natural insect populations. Annu. Rev. Entomol. 14, 175–196 (1969).Article 

Google Scholar 
Leslie, P. H. On the use of matrices in certain population mathematics. Biometrika 33, 183–212 (1945).MathSciNet 
CAS 
PubMed 
MATH 
Article 

Google Scholar 
Birch, L. The intrinsic rate of natural increase of an insect population. J. Anim. Ecol., 15–26 (1948).Chi, H. Life-table analysis incorporating both sexes and variable development rates among individuals. Environ. Entomol. 17, 26–34 (1988).Article 

Google Scholar 
Chi, H. & Liu, H. Two new methods for the study of insect population ecology. Bull. Inst. Zool. Acad. Sin 24, 225–240 (1985).
Google Scholar 
Fathipour, Y. & Maleknia, B. in Ecofriendly Pest Management for Food Security (ed Omkar) 329–366 (Academic Press, 2016).Auad, A. et al. The impact of temperature on biological aspects and life table of Rhopalosiphum padi (Hemiptera: Aphididae) fed with signal grass. Fla. Entomol. 569–577 (2009).Qu, Y. et al. Sublethal and hormesis effects of beta-cypermethrin on the biology, life table parameters and reproductive potential of soybean aphid Aphis glycines. Ecotoxicology 26, 1002–1009 (2017).CAS 
PubMed 
Article 

Google Scholar 
Araujo, E. S., Benatto, A., Mogor, A. F., Penteado, S. C. & Zawadneak, M. A. Biological parameters and fertility life table of Aphis forbesi Weed, 1889 (Hemiptera: Aphididae) on strawberry. Braz. J. Biol. 76, 937–941. https://doi.org/10.1590/1519-6984.04715 (2016).CAS 
Article 
PubMed 

Google Scholar 
Krishnamoorthy, S. V. & Mahadevan, N. R. Life table studies of sugarcane scale, Melanaspis glomerata G. J. Entomol. Res. 27, 203–212 (2003).
Google Scholar 
Uematsu, H. Studies on life table for an armored scale insect, Aonidiella taxus Leonardi (Homoptera: Diaspididae). J. Fac. Agric. Kyushu Univ. (1979).Hill, M. G., Mauchline, N. A., Hall, A. J. & Stannard, K. A. Life table parameters of two armoured scale insect (Hemiptera: Diaspididae) species on resistant and susceptible kiwifruit (Actinidia spp.) germplasm. N. Z. J. Crop Hortic. Sci. 37, 335–343 (2009).Article 

Google Scholar 
Yong, C. X. W. Z. C. & Shaoyun, Z. J. Y. S. W. Age-specific life table of chinese white wax scale (Ericerus pela) natural population and analysis of death key factors. Scientia Silvae Sinica 9 (2008).Rosado, J. F. et al. Natural biological control of green scale (Hemiptera: Coccidae): a field life-table study. Biocontrol. Sci. Technol. 24, 190–202 (2014).Article 

Google Scholar 
Fand, B. B., Gautam, R. D., Chander, S. & Suroshe, S. S. Life table analysis of the mealybug, Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae) under laboratory conditions. J. Entomol. Res. 34, 175–179 (2010).
Google Scholar 
Vargas-Madríz, H. et al. Life and fertility table of Bactericera cockerelli (Hemiptera: Triozidae), under different fertilization treatments in the 7705 tomato hybrid. Rev. Chil. entomol. 39 (2014).Huang, Y. B. & Chi, H. Age-stage, two-sex life tables of Bactrocera cucurbitae (Coquillett)(Diptera: Tephritidae) with a discussion on the problem of applying female age-specific life tables to insect populations. Insect Sci. 19, 263–273 (2012).Article 

Google Scholar 
Saska, P. et al. Leaf structural traits rather than drought resistance determine aphid performance on spring wheat. J. Pest. Sci. 94, 423–434 (2021).Article 

Google Scholar 
Ma, K., Tang, Q., Xia, J., Lv, N. & Gao, X. Fitness costs of sulfoxaflor resistance in the cotton aphid, Aphis gossypii Glover. Pestic. Biochem. Physiol. 158, 40–46 (2019).CAS 
PubMed 
Article 

Google Scholar 
Ullah, F. et al. Fitness costs in clothianidin-resistant population of the melon aphid, Aphis gossypii. PLoS ONE 15, e0238707 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Güncan, A. & Gümüş, E. Influence of different hazelnut cultivars on some demographic characteristics of the filbert aphid (Hemiptera: Aphididae). J. Econ. Entomol. 110, 1856–1862 (2017).PubMed 
Article 

Google Scholar 
Bailey, R., Chang, N.-T., Lai, P.-Y. & Hsu, T.-C. Life table of cycad scale, Aulacaspis yasumatsui (Hemiptera: Diaspididae), reared on Cycas in Taiwan. J. Asia Pac. Entomol. 13, 183–187 (2010).Article 

Google Scholar 
Wang, Z., Chen, Y. & Diaz, R. Temperature-dependent development and host range of crapemyrtle bark scale, Acanthococcus lagerstroemiae (Kuwana)(Hemiptera: Eriococcidae). Fla. Entomol. 102, 181–186 (2019).Article 

Google Scholar 
Zhang, Z.-J. et al. A determining factor for insect feeding preference in the silkworm, Bombyx mori. PLoS Biol. 17, e3000162 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Wang, Z., Chen, Y., Diaz, R. & Laine, R. A. Physiology of crapemyrtle bark scale, Acanthococcus lagerstroemiae (Kuwana), associated with seasonally altered cold tolerance. J. Insect Physiol. 112, 1–8 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Suh, S.-J. Notes on some parasitoids (Hymenoptera: Chalcidoidea) associated with Acanthococcus lagerstroemiae (Kuwana)(Hemiptera: Eriococcidae) in the Republic of Korea. Insecta mundi 0690, 1–5 (2019).
Google Scholar 
Meindl, G. A., Bain, D. J. & Ashman, T.-L. Edaphic factors and plant–insect interactions: Direct and indirect effects of serpentine soil on florivores and pollinators. Oecologia 173, 1355–1366 (2013).ADS 
PubMed 
Article 

Google Scholar 
Wielgolaski, F. E. Phenological modifications in plants by various edaphic factors. Int. J. Biometeorol. 45, 196–202 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Uchida, R. in Plant nutrient management in Hawaii’s soils (ed Raymond S. Uchida James A. Silva) 31–55 (University of Hawaii at Manoa, College of Agriculture & Tropical Resources, 2000).Flanders, S. E. Observations on host plant induced behavior of scale insects and their endoparasites. Can. Entomol. 102, 913–926 (1970).Article 

Google Scholar 
Yang, T.-C. & Chi, H. Life tables and development of Bemisia argentifolii (Homoptera: Aleyrodidae) at different temperatures. J. Econ. Entomol. 99, 691–698 (2006).PubMed 
Article 

Google Scholar 
Tuan, S. J., Lee, C. C. & Chi, H. Population and damage projection of Spodoptera litura (F.) on peanuts (Arachis hypogaea L.) under different conditions using the age-stage, two-sex life table. Pest Manag. Sci. 70, 805–813 (2014).CAS 
PubMed 
Article 

Google Scholar 
Vafaie, E. et al. Seasonal population patterns of a new scale pest, Acanthococcus lagerstroemiae Kuwana (Hemiptera: Sternorrhynca: Eriococcidae), of Crapemyrtles in Texas, Louisiana, and Arkansas. J. Environ. Hortic. 38, 8–14 (2020).Article 

Google Scholar 
Vafaie, E. K. Bark and systemic insecticidal control of Acanthococcus (= Eriococcus) lagerstroemiae (Hemiptera: Eriococcidae) on Potted Crapemyrtles, 2017. Arthropod manag. tests 44, tsy109 (2019).Vafaie, E. K. & Knight, C. M. J. A. M. T. Bark and systemic insecticidal control of Acanthococcus (= Eriococcus) lagerstroemiae (Crapemyrtle Bark Scale) on Landscape Crapemyrtles, 2016. 42, tsx130 (2017).Vafaie, E. & Gu, M. Insecticidal control of crapemyrtle bark scale on potted crapemyrtles, Fall 2018. Arthropod. Manag. Tests 44, tsz061 (2019).Article 

Google Scholar 
Aktar, M. W., Sengupta, D. & Chowdhury, A. J. I. t. Impact of pesticides use in agriculture: their benefits and hazards. 2, 1 (2009).Grafton-Cardwell, E. & Vehrs, S. Monitoring for organophosphate-and carbamate-resistant armored scale (Homoptera: Diaspididae) in San Joaquin valley citrus. J. Econ. Entomol. 88, 495–504 (1995).CAS 
Article 

Google Scholar 
Almarinez, B. J. M. et al. Biological control: A major component of the pest management program for the invasive coconut scale insect, Aspidiotus rigidus Reyne, in the Philippines. Insects 11, 745 (2020).PubMed Central 
Article 

Google Scholar 
Grout, T. & Richards, G. Value of pheromone traps for predicting infestations of red scale, Aonidiella aurantii (Maskell)(Hom., Diaspididae), limited by natural enemy activity and insecticides used to control citrus thrips, Scirtothrips aurantii Faure (Thys., Thripidae). J. Appl. Entomol. 111, 20–27 (1991).Article 

Google Scholar 
Grafton-Cardwell, E., Millar, J., O’Connell, N. & Hanks, L. Sex pheromone of yellow scale, Aonidiella citrina (Homoptera: Diaspididae): Evaluation as an IPM tactic. J. Agric. Urban. Entomol. 17, 75–88 (2000).CAS 

Google Scholar 
Jactel, H., Menassieu, P., Lettere, M., Mori, K. & Einhorn, J. Field response of maritime pine scale, Matsucoccus feytaudi Duc. (Homoptera: Margarodidae), to synthetic sex pheromone stereoisomers. J. Chem. Ecol. 20, 2159–2170 (1994).CAS 
PubMed 
Article 

Google Scholar 
Mendel, Z. et al. Outdoor attractancy of males of Matsucoccus josephi (Homoptera: Matsucoccidae) and Elatophilus hebraicus (Hemiptera: Anthocoridae) to synthetic female sex pheromone of Matsucoccus josephi. J. Chem. Ecol. 21, 331–341 (1995).CAS 
PubMed 
Article 

Google Scholar 
Zada, A. et al. Sex pheromone of the citrus mealybug Planococcus citri: Synthesis and optimization of trap parameters. J. Econ. Entomol. 97, 361–368 (2004).CAS 
PubMed 
Article 

Google Scholar 
Zhang, Z. & Shi, Y. Studies on the Morphology and Biology of Eriococcus Lagerstroemiae Kuwana. J. Shandong Agri. Univ. 2 (1986).Savopoulou-Soultani, M., Papadopoulos, N. T., Milonas, P. & Moyal, P. Abiotic factors and insect abundance. PSYCHE 2012 (2012).Vandegehuchte, M. L., de la Pena, E. & Bonte, D. Relative importance of biotic and abiotic soil components to plant growth and insect herbivore population dynamics. PLoS ONE 5, e12937 (2010).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Clavijo McCormick, A. Can plant–natural enemy communication withstand disruption by biotic and abiotic factors?. Ecol. Evol. 6, 8569–8582 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Nebapure, S. M. & Sagar, D. Insect-plant interaction: A road map from knowledge to novel technology. Karnataka J. Agric. Sci. 28, 1–7 (2015).
Google Scholar 
Murashige, T. & Skoog, F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant. 15, 473–497 (1962).CAS 
Article 

Google Scholar 
Hogendorp, B. K., Cloyd, R. A. & Swiader, J. M. Effect of nitrogen fertility on reproduction and development of citrus mealybug, Planococcus citri Risso (Homoptera: Pseudococcidae), feeding on two colors of coleus Solenostemon scutellarioides L. Codd. Environ. Entomol. 35, 201–211 (2006).Article 

Google Scholar 
Lema, K. & Mahungu, N. in Tropical root crops: Production and uses in Africa: proceedings of the Second Triennial Symposium of the International Society for Tropical Root Crops-Africa Branch held in Douala, Cameroon, 14-19 Aug. 1983. (IDRC, Ottawa, ON, CA).McClure, M. S. Dispersal of the scale Fiorinia externa (Homoptera: Diaspididae) and effects of edaphic factors on its establishment on hemlock. Environ. Entomol. 6, 539–544 (1977).Article 

Google Scholar 
Salama, H., Amin, A. & Hawash, M. Effect of nutrients supplied to citrus seedlings on their susceptibility to infestation with the scale insects Aonidiella aurantii (Maskell) and Lepidosaphes beckii (Newman)(Coccoidea). Zeitschrift für Angewandte Entomologie 71, 395–405 (1972).Article 

Google Scholar 
Rasmann, S. & Pellissier, L. in Climate Change and Insect Pests Vol. 8 (ed P. Niemelä C. Björkman) 38–53 (Wallingford, UK: CAB Int., 2015).Wang, Z. & Li, S. Effects of nitrogen and phosphorus fertilization on plant growth and nitrate accumulation in vegetables. J. Plant Nutr. 27, 539–556 (2004).CAS 
Article 

Google Scholar 
Da Costa, P. B. et al. The effects of different fertilization conditions on bacterial plant growth promoting traits: Guidelines for directed bacterial prospection and testing. Plant Soil. 368, 267–280 (2013).Article 

Google Scholar 
Dong, H., Kong, X., Li, W., Tang, W. & Zhang, D. Effects of plant density and nitrogen and potassium fertilization on cotton yield and uptake of major nutrients in two fields with varying fertility. Field Crops Res. 119, 106–113 (2010).Article 

Google Scholar 
Aulakh, M., Dev, G. & Arora, B. Effect of sulphur fertilization on the nitrogen–sulphur relationships in alfalfa (Medicago sativa L. Pers.). Plant Soil. 45, 75–80 (1976).CAS 
Article 

Google Scholar 
Powell, G., Tosh, C. R. & Hardie, J. Host plant selection by aphids: Behavioral, evolutionary, and applied perspectives. Annu. Rev. Entomol. 51, 309–330 (2006).CAS 
PubMed 
Article 

Google Scholar 
Sauge, M. H., Grechi, I. & Poëssel, J. L. Nitrogen fertilization effects on Myzus persicae aphid dynamics on peach: Vegetative growth allocation or chemical defence?. Entomol. Exp. Appl. 136, 123–133 (2010).CAS 
Article 

Google Scholar 
Chen, Y., Serteyn, L., Wang, Z., He, K. & Francis, F. Reduction of plant suitability for corn leaf aphid (Hemiptera: Aphididae) under elevated carbon dioxide condition. Environ. Entomol. (2019).Miller, D. R. & Kosztarab, M. Recent advances in the study of scale insects. Annu. Rev. Entomol. 24, 1–27 (1979).CAS 
Article 

Google Scholar 
Hardy, N. B., Peterson, D. A. & Normark, B. B. Scale insect host ranges are broader in the tropics. Biol. Lett. 11, 20150924 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Chen, Q. et al. Age-stage, two-sex life table of Parapoynx crisonalis (Lepidoptera: Pyralidae) at different temperatures. PLoS ONE 12, e0173380 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Li, X. et al. Density-dependent demography and mass-rearing of Carposina sasakii (Lepidoptera: Carposinidae) incorporating life table variability. J. Econ. Entomol. 112, 255–265 (2019).PubMed 
Article 

Google Scholar 
Ning, S., Zhang, W., Sun, Y. & Feng, J. Development of insect life tables: comparison of two demographic methods of Delia antiqua (Diptera: Anthomyiidae) on different hosts. Sci. Rep. 7, 1–10 (2017).ADS 
Article 

Google Scholar 
TWOSEX-MSChart: A computer program for the age-stage, two-sex life table analysis (2020).Goodman, D. Optimal life histories, optimal notation, and the value of reproductive value. Am. Nat. 119, 803–823 (1982).MathSciNet 
Article 

Google Scholar 
Efron, B. & Tibshirani, R. J. An Introduction to the Bootstrap (CRC Press, 1994).MATH 
Book 

Google Scholar  More