Philippa Kaur

Jiang, Z. H. Bamboo and Rattan in the World (China Forest Publishing House, 2007).
Google Scholar 
Zhao, J., Wang, B., Li, Q., Yang, H. & Xu, K. Analysis of soil degradation causes in Phyllostachys edulis forests with different mulching years. Forests 9(3), 149 (2018).Article 

Google Scholar 
Su, W., Fan, S., Zhao, J. & Cai, C. Effects of various fertilization placements on the fate of urea-15N in moso bamboo forests. For. Ecol. Manag. 453, 117632 (2019).Article 

Google Scholar 
Zhao, J. et al. Ammonia volatilization and nitrogen runoff losses from moso bamboo forests under different fertilization practices. Can. J. For. Res. 49(3), 213–220 (2019).CAS 
Article 

Google Scholar 
Yin, J. et al. Abandonment lead to structural degradation and changes in carbon allocation patterns in Moso bamboo forests. For. Ecol. Manag. 449, 117449 (2019).Article 

Google Scholar 
Xu, Q. F. et al. Rapid bamboo invasion (expansion) and its effects on biodiversity and soil processes. Glob. Ecol. Conserv. 21, e00787 (2020).Article 

Google Scholar 
Prayogo, C., Sholehuddin, N., Putra, E. Z. H. S. & Rachmawati, R. Soil macrofauna diversity and structure under different management of pine-coffee agroforestry system. J. Degrade. Min. Land Manage. 6(3), 1727–1736 (2019).Article 

Google Scholar 
Coleman, B. R., Martin, A. R., Thevathasan, N. V., Gordon, A. M. & Isaac, M. E. Leaf trait variation and decomposition in short-rotation woody biomass crops under agroforestry management. Agric. Ecosyst. Environ. 298, 106971 (2020).CAS 
Article 

Google Scholar 
Cai, C. J., Fan, S. H., Liu, G. L., Wang, S. M. & Feng, Y. Research and development advance of compound management of bamboo forests. World Bamboo Rattan 16(5), 47–52 (2018) (in Chinese).
Google Scholar 
Song, Z. et al. Characteristics of Se-enriched mycelia by Stropharia rugoso-annulata and its antioxidant activities in vivo. Biol. Trace Elem. Res. 113(1), 81–89 (2009).Article 

Google Scholar 
Wang, Q. et al. Effects of drying on the structural characteristics and antioxidant activities of polysaccharides from Stropharia rugosoannulata. J. Food Sci. Technol. 58, 3622–3631 (2021).CAS 
Article 

Google Scholar 
Yan, P., Jiang, J. & Cui, W. Characterization of protoplasts prepared from the edible fungus, Stropharia rugoso-annulata. World J. Microbiol. Biotechnol. 20(2), 173–177 (2004).CAS 
Article 

Google Scholar 
Frouz, J. Effects of soil macro- and mesofauna on litter decomposition and soil organic matter stabilition. Geoderma 332, 161–172 (2018).ADS 
CAS 
Article 

Google Scholar 
Lin, D. et al. Soil fauna promote litter decomposition but do not alter the relationship between leaf economics spectrum and litter decomposability. Soil Biol. Biochem. 136, 107519 (2019).CAS 
Article 

Google Scholar 
Meehan, M. L. et al. Response of soil fauna to simulated global change factors depends on ambient climate conditions. Pedobiologia 83, 150672 (2020).Article 

Google Scholar 
Tan, B. et al. Soil fauna show different degradation patterns of lignin and cellulose along an elevational gradient. Appl. Soil Ecol. 155, 103673 (2020).Article 

Google Scholar 
John, K., Zaitsev, A. S. & Wolters, V. Soil fauna groups respond differentially to changes in crop rotation cycles in rice production systems. Pedobiologia 84, 150703 (2021).Article 

Google Scholar 
Qin, Z. et al. Changes in the soil meso- and micro-fauna community under the impacts of exotic Ambrosia artemisiifolia. Ecol. Res. 34(2), 265–276 (2019).Article 

Google Scholar 
Chauvat, M., Titsch, D., Zaytesev, A. S. & Wolters, V. Changes in soil faunal assemblages during conversion from pure to mixed forest stands. For. Ecol. Manag. 262(3), 317–324 (2011).Article 

Google Scholar 
Yan, S. et al. A soil fauna index for assessing soil quality. Soil Biol. Biochem. 47(2), 158–165 (2012).CAS 
Article 

Google Scholar 
Reeve, J. R. et al. Effects of soil type and farm management on soil ecological functional genes and microbial activities. ISME J. 4, 1099–1107 (2010).Article 

Google Scholar 
Lavelle, P., Bignell, D. & Lepage, M. Soil function in a changing world: The role of invertebrate engineers. Eur. J. Soil Biol. 33, 159–193 (1997).CAS 

Google Scholar 
Zhu, X. & Zhu, B. Diversity and abundance of soil fauna as influenced by long-term fertilization in cropland of purple soil, China. Soil Till. Res. 146, 39–46 (2015).Article 

Google Scholar 
Zhang, L., Wang, G. & Cao, F. The effect of ginkgo agroforestry patterns on soil fauna diversity. J. Nanjing For. Univ. 39(2), 27–32 (2015) (in Chinese).
Google Scholar 
Liu, P. et al. Impact of straw returning on cropland soil mesofauna community in the western part of black soil area. Chin. J. Ecol. 37(1), 139–146 (2018) (in Chinese).
Google Scholar 
Liu, M. Study on the model of interplanting edible fungi under bamboo (Phyllostachys edulis) forest and comprehensive benefit comparative. Master’s Thesis, Chinese Academy of Forestry (2021) (in Chinese).Wang, B., Shen, Q., Zhu, W., Shen, X. & Li, Q. Effects of interplanting Dictyophora echinovolvata on physicochemical properties, phospholipid fatty acids characters and enzyme activities in soil of Phyllostachy heterocycla cv. pubescens. For. Environ. Sci. 32(4), 28–32 (2016) (in Chinese).Article 

Google Scholar 
Ying, G. H. et al. Effect of cultivation of Dictyophora echinovolvata on shoot yield and soil under Phyllostachy heterocycla cv. pubescens stand. J. Zhejiang For. Sci. Technol. 34(6), 65–67 (2014) (in Chinese).
Google Scholar 
Sokol, N. W. et al. Life and death in the soil microbiome: How ecological processes influence biogeochemistry. Nat. Rev. Microbiol. 20, 415–430 (2022).CAS 
Article 

Google Scholar 
Fujii, K., Hayakawa, C., Inagaki, Y. & Kosaki, T. Effects of land use change on turnover and storage of soil organic matter in a tropical forest. Plant Soil 446(1), 425–439 (2020).CAS 
Article 

Google Scholar 
Fujii, K. & Toma, T. Comparison of soil acidification rates under different land uses in Indonesia. Plant Soil 465(1–2), 1–17 (2021).CAS 
Article 

Google Scholar 
Poss, R., Smith, C. J., Dunin, F. X. & Angus, J. F. Rate of soil acidification under wheat in a semi-arid environment. Plant Soil 177, 85–100 (1995).CAS 
Article 

Google Scholar 
Yin, X. et al. Distribution and diversity partterns of soil fauna in different salinization habitats of Songnen Grasslands, China. Appl. Soil Ecol. 123, 375–383 (2018).Article 

Google Scholar 
Luo, M. L. et al. Effects of different rice straw returning quantities on soil fauna community structure. J. Zhejiang A&F Univ. 37(1), 85–92 (2020) (in Chinese).
Google Scholar 
Peng, C. Y. et al. Community structure characteristics of medium- and small-sized soil faunas in typical artificial plantation in the upper reaches of Yangtze River. J. Zhejiang Univ. 45(5), 585–595 (2019) (in Chinese).
Google Scholar 
Carmen, M. U., Edmond, R. Z. & Michelle, M. W. Nematode indicators as integrative measures of soil condition in organic cropping systems. Soil Biol. Biochem. 64, 103–113 (2013).Article 

Google Scholar 
Kamau, S., Karanja, N. K., Ayuke, F. O. & Lehmann, J. Short-term influence of biochar and fertilizer-biochar blends on soil nutrients, fauna and maize growth. Biol. Fertil. Soils 55(7), 661–673 (2019).CAS 
Article 

Google Scholar 
Fu, X., Shao, M., Wei, X. & Horton, R. Soil organic carbon and total nitrogen as affected by vegetation types in Northern Loess Plateau of China. Geoderma 155(1–2), 31–35 (2010).ADS 
CAS 
Article 

Google Scholar 
Guan, F., Tang, X., Fan, S., Zhao, J. & Peng, C. Changes in soil carbon and nitrogen stocks followed the conversion from secondary forest to Chinese fir and Moso bamboo plantations. Catena 133, 455–460 (2015).CAS 
Article 

Google Scholar 
Liu, Y. et al. Higher soil fauna abundance accelerates litter carbon release across an alpine forest-tundra ecotone. Sci. Rep. 9, 10561 (2019).ADS 
CAS 
Article 

Google Scholar  More

Systematic palaeontologyOrder Hemiptera Linnaeus, 1758Suborder Sternorrhyncha Amyot et Audinet-Serville, 1843Superfamily Protopsyllidioidea Carpenter, 1931Family Postopsyllidiidae Hakim, Azar et Huang, 2019Genus Megalophthallidion Drohojowska et Szwedo, gen. nov.LSID urn:lsid:zoobank.org:act:A6F71390-9B8E-4A19-8F30-C2A024B6EFB1Type speciesMegalophthallidion burmapateron Drohojowska et Szwedo, sp. nov.; by present designation and monotypy.EtymologyGeneric name is derived from Classic Greek megas (μέγας)—large, ophthalmos (ὀφθαλμός)—an eye and Greek form of generic name Psyllidium. Gender: masculine.Type localityNorthern Myanmar: state of Kachin, Noije bum 2001 Summit Site amber mine in the Hukawng Valley, SW of Maingkhwan.Type stratumLowermost Cenomanian, Upper Cretaceous (‘mid-Cretaceous’).DiagnosisHead capsule with 12 stiff setae on tubercles (18 setae in Postopsyllidium); fore wing without pterostigma (tiny pterostigma, widening of ScP + RA present in Postopsyllidium); vein CuP not thickened distally (distinctly thickened distally in Postopsyllidium); profemur with a row of ventral (ventrolateral) setae (two rows in Postopsyllidium).Megalophthallidion burmapateron Drohojowska et Szwedo, sp. nov.LSID urn:lsid:zoobank.org:act:F3F971F4-AE04-4F41-98B0-9A0A04470625.(Figs. 1A–F, 2A–I).Figure 1Megalophthallidion burmapteron gen. et sp. nov., holotype (MAIG 6687), imago. (A) Photo of body, ventral side; (B) photo of right antennae and (C) drawing of antenna; (D) drawing of body, dorsal side; (E) drawing of thorax structure with sclerites marked: red—pronotum; orange—mesopraescutum; yellow—mesoscutum; light green—mesoscutellum, dark green—mesopostnotum; light blue—metascutum; dark blue—metascutellum; violet—metapostnotum; (F) photo of thorax dorsal side. Scale bars: 0.5 mm (A), 0.2 mm (B–D), 0.1 mm (F).Full size imageFigure 2Megalophthallidion burmapteron gen. et sp. nov., holotype (MAIG 6687), imago. (A) Photo of right fore wing; (B) photo of right wings; (C) photo of antenna and proleg; (D) photo of proleg and mesoleg, and (E) photo of femur of proleg, and (F) photo of right metatarsus and left mesotarsus in the background, and (G) photo of right mesotarsus of mesoleg, and (H) Photo of tarsi; (I) photo of male genital block. Scale bars: 0.5 mm (A–D), 0.2 mm (B,E,F,H), 0.1 mm (G,I).Full size imageMaterialHolotype, number MAIG 6687 (BUB 96), deposited in Museum of Amber Inclusions (MAIG), University of Gdańsk, Poland. Imago, a complete and well-preserved male. Piece of amber 8 × 6 × 3 mm, cut from larger lump, polished flat on both sides.Type localityNorthern Myanmar: state of Kachin, Noije bum 2001 Summit Site amber mine in the Hukawng Valley, SW of Maingkhwan.Type stratumLowermost Cenomanian, Upper Cretaceous (‘mid-Cretaceous’).DiagnosisAs for the genus with the following additions: three ocelli distinct, antennomere IX the longest, about as long as pedicel, antennomeres III–VII and XI of similar length, antennomere XII the shortest, subconically tapered in apical portion. Paramere lobate, ventral margin with acute, small process, apical and dorsal margins rounded. Aedeagus geniculately bent at base, directed dorsally, tapered apicad.DescriptionMale (Figs. 1A–F, 2A–I). Head with compound eyes distinctly wider than pronotum (Fig. 1D–F). Compound eyes subglobular, protruding laterally. Vertex short in midline, about 2.5 times as wide as posterior margin and as long in middle; trapezoidal, anterior margin slightly arched, lateral margins diverging posteriad, posterior margin shallowly arched, disc of vertex with distinct setae on large tubercles: four setae at posterior margin, two at anterior angles of compound eyes, two medial, over the median ocellus. Three ocelli present, median ocellus distinct, visible from above, lateral ocelli near anterior angles of compound eyes. Frons about as wide as long in midline, two rows of setae on tubercles, upper row at level of median ocellus, lower one, below half of compound eye height. Clypeus, elongate, triangular, in lower portion roof-like; two setae on tubercles near upper margin. Genae very narrow. Rostrum reaching slightly beyond mesocoxae, apical segment slightly shorter than subapical one, darker. Antennae bases placed at lower margin of compound eyes; antennal fovea elevated; scapus shorter than pedicel, cylindrical; pedicel cylindrical; antennomeres IIIrd–VIIth and XIth of similar length, VIIIth slightly longer than VIIth, as long as Xth antennomere, IXth the longest, XIIth the shortest, tapered apically; rhinaria absent.Thorax (Fig. 1D–F): pronotum quadrangular, about as long as mesothorax; pronotum with anterior and posterior margins parallel, merely arcuate, disc with transverse groove in the median portion, lateral margins slightly arcuate, two distinct setae on tubercles in anterolateral angle, two setae on tubercles anterior margin at distance1/3 to median line, three distinct setae on tubercles in posterolateral angles. Mesopraescutum subtriangular, with apex widely rounded, about 0.4 times as wide as pronotum, about 0.4 times as long as wide, delicately separated from mesoscutum. Mesoscutum as wide as pronotum at widest point, distinctly narrowed medially, anterior angles rounded, anterolateral margin sigmoid, lateral angle acute, posterior angles wide, posterior margin V-shape incised, posterolateral areas of mesoscutum disc declivent posteriorly; disc with two setae on tubercles, at 1/3 of mesoscutum width. Mesoscutellum about as long as wide, diamond-shape, anterior and lateral angles acute, posterior angle rounded. Mesopostnotum in form of transverse band, slightly widened in median portion. Metascutum narrower than mesoscutum, anterior angles widely rounded, lateral angles acute, anterolateral margin concave, posterior margin arcuate, with deep median arcuate incision. The suture between metascutum and metascutellum weakly visible, metascutellum subtriangular, longer than wide at base.Parapteron with three distinct setae.Fore wing (Fig. 2A,B) membranous, narrow, elongate, about 3.5 times as long as wide, widest at 2/3 of length. Anterior margin merely arcuate, slightly bent at very base, anteroapical angle widely arcuate, apex rounded, posteroapical angle widely arcuate, tornus arcuate, claval margin straight, with incision between terminals of Pcu (claval apex) and A1. Stem ScP + R + MP + CuA slightly arcuate, very short stalk ScP + R + MP + CuA leaving basal cell, stem ScP + R oblique, straight, forked in basal half of fore wing length, branch ScP + RA, oblique, reaching anterior margin slightly distally of half of fore wing length, slightly distally of ending of CuA2 branch; branch RP slightly arcuate, a little more curved in basal section, reaching margin at anteroapical angle; stalk MP + CuA slightly shorter than basal cell; stem MP almost straight, forked in apical half of fore wing, at about 2/3 of fore wing length, with three terminals reaching margin between apex and posteroapical angle; stem CuA shorter than branches CuA1 and CuA2, about half as long as branch CuA1; claval vein CuP weak at base, not thickened distally; claval vein Pcu straight, claval vein A1 straight. Basal cell present, subtriangular, about twice as long as wide, basal veinlet cua-cup oblique, no other veinlets present; cell r (radial) very long, longer than half of fore wing length; cell m (medial) the shortest, shorter than cell cu (areola postica). Margins of fore wing with fringe of long setae, starting on costal margin near base of fore wing, ending at level of middle of cell cu; longitudinal veins with distinct, scarcely but evenly dispersed, movable setae; terminal section of CuP with two setae; costal margin with row of short, densely distributed setae, apical margin, tornus and claval margin with rows of scaly setae.Hind wing (Fig. 2B) membranous, shorter than fore wing, 3.23 times as long as wide. Costal margin bent at base, then almost straight up to the level of ScP + RA end and wing coupling lobe, then straight to anteroapical angle, anteroapical angle widely arcuate, apex arcuate, posteroapical angle arcuate, tornus straight, claval margin merely arcuate, posteroclaval angle angulate; stem ScP + R + MP bent at base, then straight, stem ScP + R short, branch ScP + RA short, about as long as stem ScP + R, branch RP arcuate basally than straight, reaching apex; stem MP arcuate, forked slightly distad CuA1 terminus level, branch MP1+2 slightly arcuate, reaching margin at posteroapical angle, branch MP3+4 straight, reaching tornus; stem CuA slightly bent at base, then straight, forked slightly distad ScP + R forking, branch CuA1 arcuate, branch CuA2 short, straight, slightly oblique, reaching tornus; claval vein CuP weak, visible only at base, claval vein Pcu slightly arcuate; wing coupling apparatus (fold) with a few short setae.Legs slender, relatively long, profemora armed (Fig. 2C–H). Procoxa as long as profemur, narrow, flattened. Protrochanter scaphoid, elongate, with long apical and subapical setae. Profemur flattened laterally, about as long as protibia, ventrally armed with four large setae on elevated plinths; dorsal margin with row of short, decumbent setae. Protibia narrow, rounded in cross section, covered with short setae, a few longer setae in distal portion. Protarsus—single, long tarsomere, plantar surface with row of semi-erect setae; tarsal claws long, straight, directed ventrally, no arolium nor empodium.Mesocoxa elongate, narrow, slightly flattened. Mesotrochanter scaphoid. Mesofemur slender, flattened laterally, dorsal margin with short setae. Mesotibia subequal to mesofemur, slender, covered with setae, two apical setae slightly thicker and longer. Mesotarsus with three tarsomeres, basimesotarsomere the longest, shorter than cumulative length of mid- and apical mesotarsomere, plantar margins with setae, two apical setae slightly longer and thicker; midmesotarsomere the shortest, 1/3 of basimesotarsomere length, a few setae on plantar surface; apical tarsomere shorter than basimesotarsomere, twice as long as midmesotarsomere, plantar surface with a few, scarcely dispersed setae, tarsal claws long, narrow, directed ventrally, no arolium nor empodium.Metacoxa conical, narrow. Metatrochanter scaphoid, elongate. Metafemur slender, laterally flattened, longer than mesofemur, dorsal margin with row of short setae. Metatibia, long, slender, 1.6 times as long as metafemur, with suberect setae of different size, two larger and longer and two shorter setae subapical setae. Metatarsus slightly less than half of metatibia length, with three tarsomeres, basimetatarsomere the longest, more than twice as long as apical metatarsomere, 1.5 times as long as combined length of mid- and apical metatarsomere, plantar surface with scarce decumbent setae; mid metatarsomere the shortest, 1/4 of basimetatarsomere length, plantar surface with a few setae, two apical ones slightly thicker; apical metatarsomere about 0.4 of basimetatarsomere length, with scarcely dispersed setae on along plantar surface; tarsal claws, long, slender, other pretarsal structures absent.Abdomen (Fig. 1F) narrowly attached to thorax, tergite segment shorter, 2nd tergite distinctly longer, 3rd to 8th tergites of similar length; pygofer narrowing apicad, ventral margin strongly elongated posteriorly; anal tube short, directed posterodorsad, anal style shorter than anal tube. Paramere lobate, ventral margin with acute, small process, apical and dorsal margins rounded. Aedeagus (Fig. 2I) geniculately bent at base, directed dorsad, tapered apicad.Female. Unknown.Megalophthallidion sp. (5th instar nymph)(Figs. 3A–D, 4A–F)Figure 3Megalophthallidion sp. (MAIG 6688), nymph. (A) Photo of body, dorsal side and (B) drawing of body dorsal side; (C) photo of body dorsal side and (D) drawing of body ventral side. Scale bars: 0.5 mm (A–D).Full size imageFigure 4Megalophthallidion sp. (MAIG 6688), nymph. Photo of clypeus and (B) drawing of clypeus; (C) photo of proleg, and (D) photo of mesoleg, and (E) photo of metaleg; (F) photo of posterior part of abdomen ventral side. Scale bars: 0.1 mm (A–F).Full size imageMaterialNymph, 5th instar, MAIG 6688 (BUB 1799), deposited in Museum of Amber Inclusions (MAIG), University of Gdańsk, Poland. Piece of amber 13 × 6 × 2 mm, cut from larger lump, polished flat on one side, more convex on the other.Diagnostic charactersThe nymph of Megalophthallidion gen. nov. is similar in general body shape to the only known fossil protopsyllidioidean nymph described from Lower Cretaceous Lebanese amber—Talaya batraba Drohojowska et Szwedo, 2013. The nymph of Talaya batraba is 2nd or 3rd instar, therefore some features are difficult to compare with this last instar nymph of Megalophthallidion gen. nov. The morphological states observed in those two specimens are: head covered with strongly expanded disc and expanded disc of pronotum, however shapes and ratios of these structures differ; compound eyes on ventral side of head, shifted laterad (ommatidia on cones in T. batraba, while ventroposterior expansions are present in Megalophthallidion gen. nov.); compound eyes visible from above as short, stout cones in fissure between posterior margin of disc (hypertrophied vertex) and anterior margin of pronotum (compound eyes (?) are visible on dorsal side of Permian Aleuronympha bibulla Riek, 1974); in Megalophthallidion gen. nov. rostrum reached mesocoxa, while in Talaya batraba distinctly exceeds length of the body; abdomen with 9 segments; tergites of abdominal segments 5th–9th expanded posterolaterad in form of fan-like expansion; 9th abdominal segment short, placed ventral; anal tube short, cylindrical, epiproct (?) globular.DescriptionNymph, 5th instar (Figs. 3A–D, 4A–F). Body oval shaped, dorso-ventrally flattened, 1.5 times longer than wide with segmentation visible; on the ventral side slightly concave. Length of body c. 1.56 mm long, outline, in dorsal view, maximum width of body 0.94 mm; length of head and pronotum (cephaloprothorax) c. 0.46 mm in midline, width c. 0.83 mm; cumulative length of mesonotum + metanotum c. 0.25 mm; abdomen c. 0.8 mm long. Dorsal side (Fig. 3A,B) with distinct median line (ecdysial line), not reaching anterior or posterior margin of the body, the line distinctly roof-like in abdominal portion. Anterior margin of head (cephaloprothorax) disc arcuate, lateral angles rounded; anterior margin of pronotum arcuate, lateral margins arcuately diverging posteriad, posterior margin distinctly arcuate, anterior angles widely rounded, posterior angles acutely rounded, disc elevated, convex, lateral portions declivitous; the fissure between posterior margin of head disc and anterior margin of pronotum narrow, widened medially, with stalked compound eyes popping out.Head partly separated from prothorax, wide in ventral view. Bases of antennae protruding anterolaterally, wide, anterior margin arcuate, with a small lump extending anteriorly connecting margin with vertex expansion. Suture separating anteclypeus and postclypeus visible in ventral aspect (Fig. 4A,B). Postclypeus about three times as long as wide, oval, slightly swollen, without any setae; weak traces of salivary pump muscle attachments visible. Anteclypeus about as long as postclypeus, widened in upper section below clypeal suture, convex, carinately elevated in lower section, with sides distinctly declivitous, clypellus long, carinately elevated. Lora (mandibulary plates) distinct, separated from anteclypeus by shallow suture, with upper angles at half of postclypeus length, lower angles at half of anteclypeus length, about as wide as half of postclypeus width. Maxillary plates narrow. Genal portion of head enlarged, medial portion arcuately convex; lateral sections narrowing laterally, terminally encircling bases of compound eyes. Antennae short (Fig. 3C,D), placed in front of genal portion. Antennal flagellum indistinctly subdivided into four segments. Rostrum (Fig. 4A,B) three-segmented, 0.2 mm long, with apex reaching apex of mesocoxae; apical segment about 2.5 times as long as subapical one.No lateral sclerites on meso- and metathorax, only one plus one large medial sclerite on both meso- and metathorax. Mesothoracic and metathoracic wing pads distinct, wide, subtriangular, with posterior apices directed posteriorly; lateral portions of mesothoracic wing pads arcuate. Fore wing pad 0.6 mm long, with small, straight humeral lobe, forming a right angle, not protruding anteriorly. Mesothoracic tergites slightly larger than metathoracic segments (respectively c. 0.14 mm and c. 0.12 mm long in midline, 0.26 mm and 0.27 mm in lateral lines); mesothoracic tergum with distinct median elevation (low double crest with ecdysial line in between), slightly wider than long in midline, anterior margin arcuate, lateral margins straight, subparallel, posterior margin concave. Metathoracic wing pad apex slightly exceeding mesothoracic wing pad. Metathoracic tergum wider than long, slightly shorter than mesothoracic tergum, with distinct elevation in the middle.Legs relatively long (Figs. 3C,D, 4C–E). Coxae of legs placed near the median axis of the body. Prolegs: procoxal pit with margins elevated, procoxa conical (c. 0.1 mm long), protrochanter scaphoid, about as long as procoxa, profemur c. 0.13 mm long, slightly flattened laterally, merely thickened, protibia longer than profemur, c. 0.23 mm long; tarsus shorter than protibia, basiprotarsomere about as long as apical protarsomere, the latter with distinct tarsal claws, and wide arolium. Mesoleg similar to proleg, mesocoxa conical (c. 0.1 mm long), mesotrochanter scaphoid, mesofemur (c. 0.13 mm) slightly flattened laterally, mesotibia slightly longer than mesofemur (c. 0.18 mm), mesotarsus slightly shorter than mesotibia, three-segmented, basimesotarsomere the longest (c. 0.07 mm), about as long as combined length of mid- and apical mesotarsomeres (c. 0.04 mm respectively), arolium wide, tarsal claws distinct. Metaleg: metacoxa conical (c. 0.1 mm), metatrochanter scaphoid, about as long as metacoxa (c. 0.12 mm). Metafemur (c. 0.17 mm) slightly more thickened than pro- and mesofemur, metatibia slightly longer (0.19 mm) than pro- and mesotibiae. Metatarsus three-segmented: basimetatarsomere about as long (0.08 mm) as combined length of mid- and apical metatarsomeres (0.04 mm respectively), arolium lobate, wide, tarsal claws distinct, widely spread.Abdomen (Fig. 3A–D) 9-segmented, narrow at base, widening fan-shape posteriorly, 1st segment visible from above, segmentation visible, abdominal terga 5th–9th expanded posterolaterally. Tergites carinately elevated in the middle, separated by ecdysial line. 1st sternite visible in ventral view, sternites 2nd–4th fused medially, sternites 5th–9th separated; 9th abdominal segment short (Fig. 4F), placed ventrally, under tergal expansion; anal tube short, cylindrical, epiproct (?) globular. More

Analysis of contents of heavy metals in wasteland soilThe test results show (Table 5) that the contents of Hg, Cd, As, Pb, Cr, Zn, Ni and Cu in the surface soil within Shigetai Coal Mine vary from 0.043 to 0.255, 0.44 to 2.23, 2.66 to 18.40, 11.80 to 42.80, 40.50 to 118.60, 18.90 to 70.10, 4.31 to 28.10, 4.96 to 46.25 mg/kg, respectively; the average contents of Hg, Cd, As, Pb, Cr, Zn, Ni and Cu are 0.128, 1.03, 4.73, 23.08, 76.22, 46.94, 16.11 and 12.10 mg/kg, respectively. The average contents of Hg, Cd, Pb and Cr in soil within the research area are 2.03, 1.36, 1.11 and 1.23 times of the soil background values in Shaanxi Province, respectively. The average contents of As, Zn and Cu are lower than the soil background value in Shaanxi Province, but the maximum contents of these three elements are 1.65, 1.01 and 2.16 times of the soil background values in Shaanxi Province, respectively. It is reported that the average concentration of lead in agricultural soil affected by coal mines is relatively high (433 mg kg−1)38. Lead is usually related to minerals in coal and occurs mainly in the form of sulfide such as PbS and PbSe39. In addition, aluminosilicate and carbonate also contain lead40. Chromium is a non-volatile element, which is related to aluminosilicate minerals41. In the mining process, chromium may be accumulated in coal, gangue or other tailings, and then enter the soil or water body through rain leaching42.Table 5 Statistics of contents of heavy metals in wasteland soil (n = 79).Full size tableThe coefficient of variation (CV) of Hg and Cd contents in soil within the research area is 0.050 and 0.37, respectively, with moderate variation, indicating that the content of these two heavy metals is less affected by the external factors; the coefficient of variation (CV) of As, Pb, Cr, Zn, Ni and Cu contents is 2.81, 7.46, 18.00, 13.51, 5.44 and 5.64, respectively, with strong variation (CV  > 0.50)43, indicating that the content of these eight heavy metals may be affected by some local pollution sources. The skewness coefficient (SK) ranges from − 3 to 3, and the larger its absolute value, the greater its skewness. When SK  > 0, it is positive skewness; when SK  More

Filgueira, R. R., Fournier, L. L., Cerisola, C. I., Gelati, P. & Garcia, M. G. Particle-size distribution in soils: A critical study of the fractal model validation. Geoderma 134, 327–334 (2006).ADS 
Article 

Google Scholar 
Deng, J. F., Li, J. H., Deng, G., Zhu, H. Y. & Zhang, R. H. Fractal scaling of particle-size distribution and associations with soil properties of Mongolian pine plantations in the Mu Us Desert, China. Sci. Rep. 7, 6742 (2018).ADS 
Article 

Google Scholar 
Gao, Y. J. et al. “Fertile islands” beneath three desert vegetation on soil phosphorus fractions, enzymatic activities, and microbial biomass in the desert-oasis transition zone. CATENA 212, 106090 (2022).CAS 
Article 

Google Scholar 
Zeraatpisheh, M., Ayoubi, S., Mirbagheri, Z., Mosaddeghi, M. R. & Xu, M. Spatial prediction of soil aggregate stability and soil organic carbon in aggregate fractions using machine learning algorithms and environmental variables. Geoderma Regioanl. 27, e00440 (2021).Article 

Google Scholar 
Zha, C., Shao, M., Jia, X. & Zhang, C. Particle size distribution of soils (0–500 cm) in the Loess Plateau, China. Geoderma 7, 251–258 (2016).Article 

Google Scholar 
Callesen, I., Keck, H. & Andersen, T. J. Particle size distribution in soils and marine sediments by laser diffraction using Malvern Mastersizer 2000-method uncertainty including the effect of hydrogen peroxide pretreatment. J. Soils Sediments 18, 2500–2510 (2018).CAS 
Article 

Google Scholar 
He, Y. J. & Lv, D. Y. Fractal expression of soil particle-size distribution at the basin scale. Open Geosci. 14, 70–78 (2022).Article 

Google Scholar 
Besalatpour, A. A., Ayoubi, S., Hajabbasi, M. A., Mosaddeghi, M. R. & Schulin, R. Estimating wet soil aggregate stability from easily available properties in a highly mountainous watershed. CATENA 111, 72–79 (2013).Article 

Google Scholar 
Besalatpour, A. A., Ayoubi, S., Hajabbasi, M. A., Yousefian, J. A. & Gharipour, A. Feature selection using parallel genetic algorithm for the prediction of geometric mean diameter of soil aggregates by machine learning methods. Arid Land Res. Manag. 28, 383–394 (2014).Article 

Google Scholar 
Xu, G. C., Li, Z. B. & Li, P. Fractal features of soil particle-size distribution and total soil nitrogen distribution in a typical watershed in the source area of the middle Dan River, China. CATENA 101, 17–23 (2013).CAS 
Article 

Google Scholar 
Jia, W. R. et al. Grain size distribution at four developmental stages of crescent dunes in the hinterland of the Taklimakan Desert, China. J. Arid Land 8, 722–733 (2016).Article 

Google Scholar 
Rabot, E., Wiesmeier, M., Schlüter, S. & Vogel, H. J. Soil structure as an indicator of soil functions: A review. Geoderma 314, 122–137 (2018).ADS 
Article 

Google Scholar 
Ghanbarian, B. & Daigle, H. Fractal dimension of soil fragment mass-size distribution: A critical analysis. Geoderma 245–246, 98–103 (2015).ADS 
Article 

Google Scholar 
Deng, Y., Cai, C., Xia, D., Ding, S. & Chen, J. Fractal features of soil particle size distribution under different land-use patterns in the alluvial fans of collapsing gullies in the hilly granitic region of southern China. PLoS ONE 12, e0173555 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Zhai, J. Y. et al. Change in soil particle size distribution and erodibility with latitude and vegetation restoration chronosequence on the Loess Plateau, China. Int. J Environ. Res. Public Health 17, 822 (2020).CAS 
PubMed Central 
Article 

Google Scholar 
Gao, Z. Y., Niu, F. J., Lin, Z. J. & Luo, J. Fractal and multifractal analysis of soil particle-size distribution and correlation with soil hydrological properties in active layer of Qinghai-Tibet Plateau, China. CATENA 203, 105373 (2021).Article 

Google Scholar 
Chen, T. L. et al. Multifractal characteristics and spatial variability of soil particle-size distribution in different land use patterns in a small catchment of the Three Gorges Reservoir Region, China. J. Mt. Sci. 18, 111–125 (2021).Article 

Google Scholar 
Gui, D. W. et al. Characterizing variations in soil particle size distribution in oasis farmlands-a case study of the Cele Oasis. Math. Comput. Model. 51, 1306–1311 (2010).Article 

Google Scholar 
Millán, H., Gonzalez-Posada, M., Aguilar, M., Domınguez, J. & Céspedes, L. On the fractal scaling of soil data. Particle-size distributions. Geoderma 117, 117–128 (2003).ADS 
Article 

Google Scholar 
Qi, F. et al. Soil particle size distribution characteristics of different land-use types in the Funiu mountainous region. Soil Till. Res. 184, 45–51 (2018).Article 

Google Scholar 
Muhtar, Q., Hiroki, T. & Mijit, H. Formation and internal structure of Tamarix cones in the Taklimakan Desert. J. Arid Environ. 50, 81–97 (2002).Article 

Google Scholar 
Zhao, Y. J. & Xia, X. C. Research on the Relationship Between Tamarix Cone and Environmental Change in Lop Nur Region of Xinjiang 38–142 (Sci. Press, 2011) (in Chinese).
Google Scholar 
Yin, C. H., Shi, Q. M., Liang, F. & Tian, C. Y. Distribution pattern of soil salinity in Tamarix Nebkhas in Tarim Basin. Bull Soil Water Conserv. 33, 287–293 (2013) (in Chinese).
Google Scholar 
Zheng, T., Li, J. G., Li, W. H. & Wan, J. H. Soil heterogeneity and its effects on plant community in oasis desert transition zone in the lower peaches of Tarim River. J. Desert Res. 30, 128–134 (2010) (in Chinese).
Google Scholar 
Liu, J. H., Wang, X. Q., Ma, Y. & Tan, F. Z. Spatial variation of soil salinity on Tamarix ramosissima nebkhas and interdune in oasis-desert ecotone. J. Desert Res. 36, 181–189 (2016) (in Chinese).CAS 

Google Scholar 
Dong, Z. W. et al. Stoichiometric features of C, N, and P in soil and litter of Tamarix cones and their relationship with environmental factors in the Taklimakan Desert, China. J. Soils Sediments 20, 690–704 (2020).CAS 
Article 

Google Scholar 
Dong, Z. W., Li, S. Y., Mao, D. L. & Lei, J. Q. Distribution pattern of soil grain size in Tamarix sand dune in the southwest of Gurbantunggut Desert. J. Soil Water Conserv. 35, 64-72/79 (2021) (in Chinese).
Google Scholar 
Dong, Z. W., Zhao, Y., Lei, J. Q. & Xi, Y. Q. Distribution pattern and influencing factors of soil salinity at Tamarix cones in the Taklimakan Desert. Chin. J. Plant Eco. 42, 873–884 (2018) (in Chinese).Article 

Google Scholar 
Xu, L. S. et al. Oasis microclimate effect on the dust deposition in Cele Oasis at southern Tarim Basin, China. Arab J. Geosci. 9, 294 (2016).Article 

Google Scholar 
Liu, J. H., Wang, X. Q., Ma, Y. & Tan, F. Z. Spatial heterogeneity of soil grain size on Tamarix ramosissima nebkhas and interdune in desert-oasis ecotone. J. Beijing For. Univ. 37, 89–99 (2015) (in Chinese).CAS 

Google Scholar 
Mao, D. L. et al. Fractal characteristics of grain size of sand and dust in aeolian sand movement in Cele oasis-desert ecotone in Xinjiang, China. Acta Pedol. Sinica 55, 88–99 (2018) (in Chinese).
Google Scholar 
Li, J. R. & Ravib, S. Interactions among hydrological-aeolian processes and vegetation determine grain-size distribution of sediments in a semi-arid coppice dune (nebkha) system. J. Arid Environ. 154, 24–33 (2018).ADS 
Article 

Google Scholar 
Ayoubi, S., Karchegani, P. M., Mosaddeghi, M. R. & Honarjoo, N. Soil aggregation and organic carbon as affected by topography and land use change in western Iran. Soil Till. Res. 121, 18–26 (2012).Article 

Google Scholar 
Wang, X. M., Dong, Z. B., Zhang, J. W. & Chen, G. T. Geomorphology of sand dunes in the Northeast Taklimakan Desert. Geomorphology 42, 183–195 (2002).ADS 
Article 

Google Scholar 
Liu, W. G. et al. Onset of permanent Taklimakan Desert linked to the mid-Pleistocene transition. Geology 48, 782–786 (2020).ADS 
CAS 
Article 

Google Scholar 
Yang, X. H. et al. Characteristics of soil particle size distribution and its effect on dust emission in Taklimakan Desert. Trans. CSAE 36, 167–174 (2020) (in Chinese).
Google Scholar 
Bao, S. D. Soil agricultural chemistry analysis 152–200 (China Agr. Press, 2000) (in Chinese).
Google Scholar 
Folk, R. L. & Ward, W. C. Brazos Riverbar: A study in the significance of grain size parameters. J. Sediment Petrol. 27, 3–26 (1957).ADS 
Article 

Google Scholar 
Weil, R. R. & Brady, N. C. The Nature and Properties of Soils 15th edn. (PrenticeHall Press, 2017).
Google Scholar 
Churchman, G. J. Game changer in soil science. Functional role of clay minerals in soil. J. Plant Nutr. Soil Sci. 181, 99–103 (2018).CAS 
Article 

Google Scholar 
Tyler, S. W. & Wheatcraft, S. W. Fractal scaling of soil particle size distributions: Analysis and limitations. Soil Sci. Soc. Am. J. 56, 362–369 (1992).ADS 
Article 

Google Scholar 
Lin, Y. C., Mu, G. J., Xu, L. S. & Zhao, X. The origin of bimodal grain-size distribution for aeolian deposits. Aeolian Res. 20, 80–88 (2016).ADS 
Article 

Google Scholar 
Sha, G. L., Wei, T. X., Chen, Y. X., Fu, Y. C. & Ren, K. Characteristics of soil particle size distribution of typical plantcommunities on the hilly areas of Loess Plateau. Arid Land Geogr. https://doi.org/10.12118/j.issn.1000-6060.2021.487 (2022) (in Chinese).Article 

Google Scholar 
Yang, J. D., Li, G. J., Dai, Y., Rao, W. B. & Ji, J. F. Isotopic evidences for provenances of loess of the Chinese Loess Plateau. Earth Sci. Front. 16, 195–206 (2009) (in Chinese).
Google Scholar 
Wu, L. & Zhang, Y. M. Precipitation and soil particle size co-determine spatial distribution of biological soil crusts in the Gurbantunggut Desert, China. J. Arid Land 10, 701–711 (2018).Article 

Google Scholar 
Li, X. B. et al. Relationship between soil particle size distribution and soil nutrient distribution characteristics in typical communities of desert grassland. Actabot. Boreal-Occident Sin. 37, 1635–1644 (2017) (in Chinese).
Google Scholar 
Gao, G. L. et al. Fractal scaling of particle size distribution and relationships with topsoil properties affected by biological soil crusts. PLoS ONE 9, e88559 (2014).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhao, Y., Feng, Q. & Yang, H. Soil salinity distribution and its relationship with soil particle size in the lower reaches of Heihe River, Northwestern China. Environ. Earth Sci. 75, 1–18 (2016).ADS 
Article 

Google Scholar 
Zhang, X. Y. et al. Sources of Asian dust and role of climate change versus desertification in Asian dust emission. Geophys. Res. Lett. 30, 2272 (2003).ADS 
Article 

Google Scholar  More

The bee health modelThe conceptual model of the interactions of various stressors with honey bee health is described by the following system of ordinary differential equations (ODEs)$${{tau }_{{HB}}dot{x}}_{{HB}}= {-{delta }_{{HB}}x}_{{HB}}+{g}_{{TC}}left({x}_{{TC}}right)+{g}_{{VA}}left({x}_{{VA}}right)+{g}_{{VI}}left({x}_{{VI}}right) \ +{bar{f}}_{S,C}left({u}_{S},{u}_{C},{x}_{{TC}},{x}_{{VA}}right)+{bar{f}}_{P}left({u}_{P},{x}_{{TC}}right)+{underline{f}}_{{HB}}left({u}_{T}right)$$
(1)
$${{tau }_{{TC}}dot{x}}_{{TC}}={-{delta }_{{TC}}x}_{{TC}}+{g}_{{HB}}left({x}_{{HB}}right)$$
(2)
$${{tau }_{{VA}}dot{x}}_{{VA}}={-{delta }_{{VA}}x}_{{VA}}+{h}_{{VA}}left({{x}_{{HB}},x}_{{TC}},varepsilon {x}_{{VI}}right)+{underline{f}}_{{VA}}left({u}_{T}right)$$
(3)
$${{tau }_{{VI}}dot{x}}_{{VI}}={-{delta }_{{VI}}x}_{{VI}}+{h}_{{VI}}left({{x}_{{HB}},x}_{{TC}},{varepsilon x}_{{VI}}right)$$
(4)
for the state variables ({x}_{{HB}}) representing honey bee health, ({x}_{{TC}}) the stress due to toxic compounds (e.g., neonicotinoid insecticides), ({x}_{{VA}}) the stress due to parasites (e.g., V. destructor) and ({x}_{{VI}}) the stress due to pathogens (e.g., DWV). The system includes the effects of external inputs as sugar ({u}_{S}), pollen ({u}_{P}), absolute deviation from desired temperature ({u}_{T}) and sub-optimal temperature ({u}_{C}). All the inputs and possible parameters are non-negative; the coefficients (tau) denote the time constants; the coefficients (delta) denote the self-regulation parameters; (varepsilon) in the last two equations allows to account for pathogens that can ((varepsilon , > , 0)) or cannot ((varepsilon=0)) impair the immune system (through link m in Fig. 1). We assume that the functions (g) are smooth, bounded, positive, convex and decreasing to 0; the functions (bar{f}) are smooth, bounded, non-negative, concave and increasing with respect to (w.r.t.) (u) arguments (vanishing only when the first (u) argument vanishes) while convex and decreasing to 0 w.r.t. (x) arguments; the functions ({underline{f}}) are smooth, bounded, non-positive and decreasing (vanishing only when (u=0)); the functions (h) are smooth, bounded, positive, convex and decreasing to 0 w.r.t. the first argument while concave and increasing w.r.t. all the other arguments. For a detailed description of the various functions, together with a summary of the biological effects they account for and a reference to the conceptual model in Fig. 1, see Supplementary Table 3.Structural analysis of the bee health modelWe describe here the structural considerations and computations that yield the structural influence matrix for the honey bee health system.The structural influence matrix (M) is defined as follows. (M) is a symbolic matrix with entries ({M}_{{ij}}) chosen among: +,−,0,?, according to the criteria described below. Consider an equilibrium point (bar{x}) and a constant perturbation (u) applied on the (j)-th system variable (small enough not to compromise the stability of the equilibrium). The equilibrium value will be modified as (bar{x}+delta bar{x}). Consider the sign of the perturbation of the (i)-th variable, (delta bar{{x}_{i}}). Then ({M}_{{ij}}) = + if (delta bar{{x}_{i}}) always has the same sign as (u); ({M}_{{ij}}=) − if (delta bar{{x}_{i}}) always has the opposite sign as (u); ({M}_{{ij}}) = 0 if always (delta bar{{x}_{i}}=0); regardless of the system parameters. Conversely, if the sign does depend on the system parameters, we set ({M}_{{ij}}) = ?.In this section we prove that the influence matrix of the honey bee health system is structurally determined, i.e., there are no “?”‘ entries in (M).We start with the following proposition.
Proposition 1
Assume that a matrix
(J)
is Hurwitz stable (i.e., all its eigenvalues have negative real part) and has the sign pattern
$${sign}left(Jright)=left[begin{array}{cccc}- & – & – & -\ – & – & 0 & 0\ – &+& – &+\ – &+& 0 & -end{array}right]$$
Then, the sign pattern of
({adj}left(-Jright))
, the adjoint of
(-J)
, is
$${sign}left({adj}left(-Jright)right)=left[begin{array}{cccc}+& – & – & -\ – &+&+&+\ – &+&+&+\ – &+&+&+end{array}right]$$
Proof To prove the statement, we just change the sign of the first variable, hence we change sign to the first row and column of matrix (J). The resulting matrix (M) is such that$${sign}left(Mright)=left[begin{array}{cccc}- &+&+&+\+& – & 0 & 0\+&+& – &+\+&+& 0 & -end{array}right]$$We observe that (M) is a Metzler matrix, namely, all its off-diagonal entries are non-negative. Moreover, the matrix is Hurwitz stable. Then, we can proceed as in the proof of Proposition 4 in a previous report16. Given a Metzler matrix that is Hurwitz stable, its inverse has non-positive entries; hence, the inverse of (-M) has non-negative entries: ({left(-Mright)}^{-1}ge 0) elementwise. Moreover, we observe that(,M) is an irreducible matrix, i.e., there is no variable permutation that brings the matrix in a block (either upper or lower) triangular form. This implies that the inverse of (-M) has strictly positive entries: ({left(-Mright)}^{-1} , > , 0) elementwise. Also, stability implies that the determinant of (-M) is positive: ({det }left(-Mright) , > , 0). Then, ({adj}left(-Mright)={left(-Mright)}^{-1}{det }left(-Mright) > 0), hence the adjoint of (-M) is also positive elementwise. To consider again the original sign of the variables, we change sign to the first row and column of ({adj}left(-Mright)), and we get the signature above for ({adj}left(-Jright)).The next step is the characterization of the structural influence matrix, which corresponds to the sign pattern of the adjoint of the negative Jacobian matrix in Proposition 1.To this aim, we first consider the linearized system and write it in a matrix-vector form$$dot{x}left(tright)={Jx}left(tright)+{e}_{j}u$$where (dot{x}left(tright)) is the time derivative of the four-dimensional vector (xleft(tright)) and ({e}_{k}), (k={{{{mathrm{1,2,3,4}}}}}), is an input vector, constant in time, with a single non-zero component, the (k)-th, equal to 1, while the scalar (u , > , 0) is the magnitude of the input. We wish to assess the (i)-th component of (xleft(tright)), ({x}_{i}left(tright)={e}_{i}^{T}xleft(tright)). If (J) is Hurwitz, as assumed, the steady-state value of variable ({x}_{i}left(tright)) due to the input perturbation ({e}_{k}) applied to the equation of variable ({x}_{k}left(tright)) is achieved for$$0=Jbar{x}+{e}_{k}u,$$namely$${x}_{i}=-{e}_{i}^{T}{J}^{-1}{e}_{k}u,$$which implies that the sign of the steady-state value ({bar{x}}_{i}) of variable ({x}_{i}) due to a persistent positive input acting on the (k)-th equation has the same sign as ({(-{J}^{-1})}_{{ik}}), the (left(i,kright)) entry of matrix ({left(-Jright)}^{-1}). Since we assume Hurwitz stability, we have that ({det }left(-Jright)) is positive, hence the sign pattern of the inverse ({left(-Jright)}^{-1}) corresponds to the sign pattern of the adjoint, ({adj}left(-Jright)). In fact, ({adj}left(-Jright)={left(-Jright)}^{-1}{det }left(-Jright)).We next consider the nonlinear system under investigation, which we write in the form$$dot{x}left(tright)=fleft(xleft(tright)right)$$and without restriction we assume that the zero vector is an equilibrium point: (0=fleft(0right)). This condition can be always achieved, without loss of generality, by a translation of coordinates. We also consider a stable equilibrium: we assume that the linearized system at the equilibrium is asymptotically stable, namely its Jacobian (J), which has the sign pattern considered in Proposition 1 above, is Hurwitz. We also assume that a constant input perturbation of magnitude (u) is applied to the system, affecting the (k)-th equation, i.e.,$$dot{x}left(tright)=fleft(xleft(tright)right)+{e}_{k}u,$$and that the perturbation is small enough to keep the state in the domain of attraction of the considered equilibrium. Due to this perturbation, a new steady state (bar{x}left(uright)) is reached that satisfies the condition$$0=fleft(bar{x}left(uright)right)+{e}_{k}u$$To determine the sign of the new equilibrium components (bar{x}left(uright)), we consider this new equilibrium vector as a function of (u) in a small interval (left[0,{x}_{{MAX}}right]). Adopting the implicit function theorem yields$$frac{d}{{dx}}bar{x}left(uright)=-J{left(uright)}^{-1}{e}_{k}u,$$where we have denoted by (Jleft(uright)) the Jacobian matrix computed at the perturbed equilibrium (bar{x}left(uright)). Hence, for (u) small enough, the sign of the derivatives of the entries of the new, perturbed equilibrium are, structurally, the same as those in the (k)-th column of matrix (-{J}^{-1}). Since, by construction, (xleft(0right)=0), this is also the sign of the elements of vector (bar{x}left(uright)), for (u) in the interval (left[0,{x}_{{MAX}}right]).We have therefore proved that the original nonlinear system describing honey bee health admits the following structural influence matrix:$$left[begin{array}{cccc}+& – & – & -\ – &+&+&+\ – &+&+&+\ – &+&+&+end{array}right]$$System equilibriaThe results concerning the system equilibria were obtained through a standard analytical treatment of the nonlinear equations describing the equilibrium conditions of the system of differential Eqs. (1), (2), (3), (4). A detailed description of methods is reported in Supplementary Methods.Laboratory experiments using honey beesTo confirm the bistability of the system representing honey bee health as affected by multiple stressors, we used data from several survival experiments, carried out in a laboratory environment according to the same standardized method, over a 6-year period (Source data file).All experiments involved Apis mellifera worker bees, sampled at the larval stage or before eclosion, from the hives of the experimental apiary of the University of Udine (46°04′54.2″N, 13°12′34.2″E). Previous studies indicated that the local bee population consists of hybrids between A. mellifera ligustica and A.m. carnica62,63. Ethical approval was not required for this study.We considered experiments on the effect of the following stressors: infection with 1000 DWV genome copies administered through the diet before pupation, feeding with a 50 ppm nicotine in a sugar solution at the adult stage, exposition to a sub-optimal temperature of 32 °C at the adult stage. All experiments were replicated 3 to 13 times, using, in total, the number of bees reported in Table 1.For the artificial infection with DWV, we collected with soft forceps individual L4 larvae from the brood cells of several combs. Groups of 20–30 of such larvae were placed in Petri dishes with an artificial diet made of 50% royal jelly, 37% distilled water, 6% glucose, 6% fructose, and 1% yeast. 25 DWV copies per mg of diet were added or not to the diet according to the experimental group (note that a bee larva at this stage consumes about 40 mg of larval food per day, thus the viral infection per bee was 1000 viral copies). After 24 h larvae were transferred onto a piece of filter paper to remove the residues of the diet and then into a clean Petri dish, where they were maintained until eclosion. At the day of emergence, bees were transferred to plastic cages in a thermostatic cabinet, where they were kept until death. The DWV extract was prepared according to previously described protocols64 and quantified according to standard methods.For the treatment with nicotine, 10 µL of pure nicotine were added to 200 g of the sugar solution used for the feeding of the caged bees, to reach the concentration of 50 ppm.Finally, to expose bees to a 32 °C temperature, the plastic cages with the adult bees were kept in a thermostatic cabinet whose temperature was set accordingly.To monitor the survival of the adult bees treated as above, they were maintained from eclosion until death in plastic cages in a dark incubator at 34.5 °C (or 32 °C, according to the experiment), 75% R.H.; two syringes were used to supply a sugar solution made of 2.4 mol/L of glucose and fructose (61% and 31%, respectively) and water, respectively; dead bees were counted daily.All the results of these experiments are reported in Source data file.All experiments were carried out during the summer months, from June to September for 6 consecutive years. Previous data indicated that, in this region, virus prevalence increases along the active season starting from very low levels in spring and reaching 100% of virus-infected honey bees by the end of the summer; virus abundance in infected honey bees follows a similar trend28. For this reason, it can be assumed that bees sampled early in the season are either uninfected or they bear only a very low viral infection level, whereas bees sampled later in the season are likely to be virus-infected, bearing moderate to high viral infections. To confirm this assumption and identify a method for filtering our data according to viral infection, we assessed viral infection in a sample of bees from the untreated control group of each experiment, by means of qRT-PCR. According to standard practice, we assumed that Ct values below 30 are indicative of an effective viral infection, whereas Ct above that threshold are more likely in virus negative bees. As expected, we found that virus prevalence increases from June to September (Supplementary Figure 1a), in such a way that up to mid July only the minority of bees can be considered as viral infected (Supplementary Figure 1b). Therefore, we classified as “early” all the samples collected up to mid July and assumed that viral infection in those samples was low; on the other hand, samples collected from mid July till September were classified as “late” and we assumed that viral infection in those samples was high.qRT-PCR analysis of viral infection was carried out as follows. At the beginning of every experiment (i.e., at day 0), two to five bees for each replication were sampled in liquid nitrogen and transferred in a −80 °C refrigerator. After defrosting of samples in RNA later, the gut of each honey bee was eliminated to avoid the clogging of the mini spin column used after. The whole body of sampled bees was homogenized using a TissueLyser (Qiagen®, Germany). Total RNA was extracted from each bee according to the procedure provided with the RNeasy Plus mini kit (Qiagen®, Germany). The amount of RNA in each sample was quantified with a NanoDrop® spectrophotomer (ThermoFisher™, USA). cDNA was synthetized starting from 500 ng of RNA following the manufacturer specifications (PROMEGA, Italy). Additional negative control samples containing no RT enzyme were included. DWV presence was verified by qRT-PCR considering as positive all samples with a Ct value lower than 30. The following primers were adopted: DWV (F: GGTAAGCGATGGTTGTTTG, R: CCGTGAATATAGTGTGAGG65). 10 ng of cDNA from each sample were analyzed using SYBR®green dye (Ambion®) according to the manufacturer specifications, on a BioRad CFX96 Touch™ Real time PCR Detector. Primer efficiency was calculated according to the formula (E={10}^{left(-1/{{{{{{rm{slope}}}}}}}-1right)*100}). The following thermal cycling profiles were adopted: one cycle at 95 °C for 10 min, 40 cycles at 95 °C for 15 s and 60 °C for 1 min, and one cycle at 68 °C for 7 min.Individual survival and colony stabilityTo investigate how the death rate of forager bees affects colony growth, a compartment model of honey bee colony population dynamics was proposed50. This model showed that death rates over a critical threshold led to colony failure. Here we modified this model to include premature death of bees at younger age, as predicted by our model of individual bee health in the presence of an immuno-suppressive virus. We show that the critical threshold found in the previously published model50 becomes a decreasing function of the death rate of the younger individuals, so that premature death (and, in turn, immune-suppression) favors colony collapse.In more details, we first summarize the results of the previously published model50 where two populations (F) (forager) and (H) (hive) of bees are considered and where conditions are provided on the mortality (m) of (F) under which the whole population collapses: namely, mathematically stated, the system admits the zero equilibrium only. Here we extend the model partitioning (H) in two categories, (Y) (younger hive bees) and (O) (older hive bees), asintroducing an early mortality factor (n) for the young population, showing how such a factor worsens the collapsing condition.The previously published model50 concerns the interaction between hive bees (H) and forager bees (F) and is described by the ODEs$$dot{H}=Lfrac{H+F}{w+H+F}-Hleft(alpha -sigma frac{F}{H+F}right)$$$$dot{F}=Hleft(alpha -sigma frac{F}{H+F}right)-{mF}.$$Above, (L) is the queen’s eggs laying rate, (w) is the rate at which (L) is reached as the total population (H+F) gets large, (alpha) is the maximum rate at which hive bees become forager bees in the absence of the latter, (sigma) measures the reduction of recruitment of hive bees in the presence of forager bees and, finally, (m) is the death rate of forager bees (while the death rate of hive bees is assumed to be negligible).We first summarize the main results in terms of a threshold value for (m) in view of colony collapse, as our further analysis will follow a similar approach. All the parameters are assumed to be positive.The search for the equilibria of the above ODEs leads to the unique nontrivial equilibrium (beyond the trivial one)$$bar{H}=frac{L}{{mJ}}-frac{w}{1+J}$$$$bar{F}=Jbar{H}$$for$$J=Jleft(mright):=frac{alpha -sigma -m+sqrt{{left(alpha -sigma -mright)}^{2}+4malpha }}{2m}.$$Note that (J) is alway positive (and, moreover, it is independent of (L) and (w)). It follows that (bar{F}) and (bar{H}) have the same sign, so that the existence of the nontrivial equilibrium is equivalent to (bar{F}+bar{H} , > , 0). It is not difficult to recover that$$bar{F}+bar{H}=frac{w}{m}left(lfrac{1+J}{J}-mright)$$where (l:=L/w) is introduced for brevity. Then if (alpha le l) we get$$bar{F}+bar{H}=frac{w}{m}left(lfrac{1+J}{J}-mright)ge frac{w}{m}left(alpha frac{1+J}{J}-mright)=frac{w}{m}left(sigma+{mJ}right) , > , 0,$$with the last equality following from$$alpha -sigma frac{J}{1+J}-{mJ}=0,$$which in turn comes from annihilating the right-hand side of the second ODE and from using (J=bar{F}/bar{H}) while searching for equilibria. We conclude that, independently of (m), the colony never collapses if the recruitment rate (alpha) of forager bees is sufficiently low.Hence, we assume (alpha , > , l). Observe that$$bar{F}+bar{H}iff l , > , Jleft(m-lright)$$guarantees existence whenever (m) is sufficiently small, viz. (mle l). Assume then (m , > , l), so that the above condition reads$$J , < , frac{l}{m-l}$$leading to the threshold condition$$m , < , bar{m}:=frac{l}{2}frac{alpha+sigma+sqrt{{left(alpha -sigma right)}^{2}+4sigma l}}{alpha -l}$$by using the definition of (J), see Eq. (2) the previously published model50.A standard stability analysis shows that, assuming (alpha,m , > , l), the nontrivial equilibrium is (globally) asymptotically stable whenever it exists (positive), i.e., whenever (m , < , bar{m}). Otherwise, the only (globally) attracting equilibrium is the trivial one, corresponding to colony collapse (see Fig. 5 for the previously published model50 or Fig. 4 for (n=0)). In the mathematical jargon, the disappearance of the positive equilibrium, for (m) exceeding (bar{m}), is referred to as a transcritical bifurcation43.Now, in view of the outcome of the analysis of our model of individual bee health, we introduce a mortality term for the younger bees. As forager bees are recruited from adult hive bees, we divide the class of hive bees (H) in younger (Y) and older (O), assuming that the former die at a rate (n), while the death rate of the latter remains negligible according to the previously published model50. Obviously, (H=Y+O). The original ODEs are consequently modified as$$dot{Y}=Lfrac{H+F}{w+H+F}-Y$$$$dot{O}=left(1-nright)Y-Hleft(alpha -sigma frac{F}{H+F}right)$$$$dot{F}=Hleft(alpha -sigma frac{F}{H+F}right)-{mF}.$$Note that the sum of the first two equations above gives$$dot{H}=Lfrac{H+F}{w+H+F}-Hleft(alpha -sigma frac{F}{H+F}right)-{nY}.$$The new negative mortality term for younger hive bees, (-{nY}), models the fact that only the younger hive bees die prematurely while the rest of the dynamics is unchanged with respect to the original model.The search for equilibria soon gives$$bar{Y}=Lfrac{bar{H}+bar{F}}{w+bar{H}+bar{F}}$$from the first ODE above, so that the remaining two equilibrium conditions lead to$$bar{H}=frac{{L}_{n}}{{mJ}}-frac{w}{1+J}$$$$bar{F}=Jbar{H}$$for the same (J) originally defined and ({L}_{n}:=Lleft(1-nright)) (note that (nin left({{{{mathrm{0,1}}}}}right)), and the case (n=0) brings us back to the original model). From this point on the analysis is the same as that previously summarized for the original model, but for replacing (L) with ({L}_{n}) and (l) with (l:=lleft(1-nright)). Consequently, by assuming (alpha,m , > , {l}_{n}) (which is less restrictive when (n , > , 0)), the threshold condition (m < bar{m}) becomes$$m , < , bar{m}left(nright):=frac{{l}_{n}}{2}frac{alpha+sigma+sqrt{{left(alpha -sigma right)}^{2}+4sigma {l}_{n}}}{alpha -{l}_{n}},$$which clearly returns the original threshold condition when (n=0). Since$$frac{dbar{m}}{{dn}}left(nright) , < , 0$$as it can be immediately verified, it follows that the critical value for (m), (bar{m}left(nright)), beyond which the colony system admits only the zero equilibrium, i.e., the transcritical bifurcation value, decreases with (n) (Fig. 4). We thus conclude that colony collapse is favored by the premature death of younger hive bees, possibly caused by a virus impairing the immune system as shown by the analysis of our model of individual bee health.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

The impacts of climate change on host–parasite dynamics are particularly complex to predict, as they involve an interplay of both physiological and behavioural factors, from both host and parasite. For example, while warming may increase parasite developmental rates and thus increase transmission, excessive heat may instead exceed thermal limits, leading to higher parasite mortality. Transmission also relates to both the distribution and abundance of host species, which may also shift under changing climates. More

Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: The role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010).PubMed 
Article 

Google Scholar 
Keuling, O., Lauterbach, K., Stier, N. & Roth, M. Hunter feedback of individually marked wild boar Sus scrofa L.: Dispersal and efficiency of hunting in northeastern Germany. Eur. J. Wildl. Res. 56, 159–167 (2010).Article 

Google Scholar 
Trondrud, L. M. et al. Fat storage influences fasting endurance more than body size in an ungulate. Funct. Ecol. 35, 1470–1480 (2021).CAS 
Article 

Google Scholar 
Wilmers, C. C. et al. The golden age of bio-logging: How animal-borne sensors are advancing the frontiers of ecology. Ecology 96, 1741–1753 (2015).PubMed 
Article 

Google Scholar 
Kukalová, M., Gazárková, A. & Adamík, P. Should i stay or should i go? The influence of handling by researchers on den use in an arboreal nocturnal rodent. Ethology 119, 848–859 (2013).Article 

Google Scholar 
Holt, R. D. et al. Estimating duration of short-term acute effects of capture handling and radiomarking. J. Wildl. Manag. 73, 989–995 (2009).Article 

Google Scholar 
Marco, I., Viñas, L., Velarde, R., Pastor, J. & Lavin, S. Effects of capture and transport on blood parameters in free-ranging mouflon (Ovis ammon). J. Zoo Wildl. Med. 28, 428–433 (1997).CAS 
PubMed 

Google Scholar 
Cattet, M., Boulanger, J., Stenhouse, G., Powell, R. A. & Reynolds-Hogland, M. J. An evaluation of long-term capture effects in ursids: Implications for wildlife welfare and research. J. Mammal. 89, 973–990 (2008).Article 

Google Scholar 
Mortensen, R. M. & Rosell, F. Long-term capture and handling effects on body condition, reproduction and survival in a semi-aquatic mammal. Sci. Rep. 10, 1–16 (2020).Article 

Google Scholar 
Soulsbury, C. D. et al. The welfare and ethics of research involving wild animals: A primer. Methods Ecol. Evol. 11, 1164–1181 (2020).Article 

Google Scholar 
Herman, J. P. et al. Regulation of the hypothalamic-pituitary- adrenocortical stress response. Compr. Physiol. 6, 603–621 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Sapolsky, R. M., Romero, L. M. & Munck, A. U. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr. Rev. 21, 55–89 (2000).CAS 
PubMed 

Google Scholar 
Sjaastad, V. Ø., Hove, K. & Sand, O. Physiology of Domestic Animals (Scandinavian Veterinary Press, 2016).
Google Scholar 
Omsjø, E. H. et al. Evaluating capture stress and its effects on reproductive success in Svalbard reindeer. Can. J. Zool. 87, 73–85 (2009).Article 

Google Scholar 
Marco, I., Viñas, L., Velarde, R., Pastor, J. & Lavin, S. The stress response to repeated capture in mouflon (Ovis ammon): Physiological, haematological and biochemical parameters. J. Vet. Med. Ser. A Physiol. Pathol. Clin. Med. 45, 243–253 (1998).CAS 
Article 

Google Scholar 
Hattingh, J., Pitts, N. I. & Ganhao, M. F. Immediate response to repeated capture and handling of wild impala. J. Exp. Zool. 248, 109–112 (1988).CAS 
PubMed 
Article 

Google Scholar 
Ortega, A. C. et al. Effectiveness of partial sedation to reduce stress in captured mule deer. J. Wildl. Manag. 84, 1445–1456 (2020).Article 

Google Scholar 
Arnemo, J. M. & Caulkett, N. Stress. In Zoo Animal and Wildlife Anesthesia and Immobilization (eds West, G. et al.) 103–109 (Blackwell Publications, 2007).
Google Scholar 
Sinclair, M. D. A review of the physiological effects of α2-agonists related to the clinical use of medetomidine in small animal practice. Can. Vet. J. 44, 885–897 (2003).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ranheim, B. et al. The effects of medetomidine and its reversal with atipamezole on plasma glucose, cortisol and noradrenaline in cattle and sheep. J. Vet. Pharmacol. Ther. 23, 379–387 (2000).CAS 
PubMed 
Article 

Google Scholar 
Carroll, G. L. et al. Effect of medetomidine and its antagonism with atipamezole on stress-related hormones, metabolites, physiologic responses, sedation, and mechanical threshold in goats. Vet. Anaesth. Analg. 32, 147–157 (2005).CAS 
PubMed 
Article 

Google Scholar 
Rode, K. D. et al. Effects of capturing and collaring on polar bears: finDings from long-term research on the southern Beaufort Sea population. Wildl. Res. 41, 311–322 (2014).Article 

Google Scholar 
Sakamoto, H., Misumi, K., Nakama, M. & Aoki, Y. The effects of xylazine on intrauterine pressure, uterine blood flow, maternal and fetal cardiovascular and pulmonary function in pregnant goats. J. Vet. Med. Sci. 58, 211–217 (1996).CAS 
PubMed 
Article 

Google Scholar 
Katila, T. & Oijala, M. The effect of detomidine (Domosedan) on the maintenance of equine pregnancy and foetal development: ten cases. Equine Vet. J. 20, 323–326 (1988).CAS 
PubMed 
Article 

Google Scholar 
Larsen, D. G. & Gauthier, D. A. Effects of capturing pregnant moose and calves on calf survivorship. J. Wildl. Manag. 53, 564 (1989).Article 

Google Scholar 
Côté, S. D., Festa-Bianchet, M. & Fournier, F. Life-history effects of chemical immobilization and radiocollars on mountain goats. J. Wildl. Manage. 62, 745–752 (1998).Article 

Google Scholar 
DelGiudice, G. D., Mech, L. D., Paul, W. J. & Karns, P. D. Effects on fawn survival of multiple immobilizations of captive pregnant white-tailed deer. J. Wildl. Dis. 22, 245–248 (1986).CAS 
PubMed 
Article 

Google Scholar 
Brivio, F., Grignolio, S., Sica, N., Cerise, S. & Bassano, B. Assessing the impact of capture on wild animals: The case study of chemical immobilisation on alpine ibex. PLoS ONE 10, e0130957 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Wingfield, J. C. et al. Ecological bases of hormone-behavior interactions: The ‘emergency life history stage’. Am. Zool. 38, 191–206 (1998).CAS 
Article 

Google Scholar 
Huber, S., Palme, R. & Arnold, W. Effects of season, sex, and sample collection on concentrations of fecal cortisol metabolites in red deer (Cervus elaphus). Gen. Comp. Endocrinol. 130, 48–54 (2003).CAS 
PubMed 
Article 

Google Scholar 
Morellet, N. et al. The effect of capture on ranging behaviour and activity of the European roe deer Capreolus capreolus. Wildlife Biol. 15, 278–287 (2009).Article 

Google Scholar 
Tarlow, E. M. & Blumstein, D. T. Evaluating methods to quantify anthropogenic stressors on wild animals. Appl. Anim. Behav. Sci. 102, 429–451 (2007).Article 

Google Scholar 
Hik, D. S. Does risk of predation influence the cyclic decline of snowshoe hares. Wildl. Res. 22, 115–129 (1995).Article 

Google Scholar 
Ordiz, A. et al. Lasting behavioural responses of brown bears to experimental encounters with humans. J. Appl. Ecol. 50, 306–314 (2013).Article 

Google Scholar 
Dechen Quinn, A. C., Williams, D. M. & Porter, W. F. Postcapture movement rates can inform data-censoring protocols for GPS-collared animals. J. Mammal. 93, 456–463 (2012).Article 

Google Scholar 
Cattet, M. R. L. Falling through the cracks: Shortcomings in the collaboration between biologists and veterinarians and their consequences for wildlife. ILAR J. 54, 33–40 (2013).CAS 
PubMed 
Article 

Google Scholar 
Albon, S. D. et al. Contrasting effects of summer and winter warming on body mass explain population dynamics in a food-limited Arctic herbivore. Glob. Change Biol. 23, 1374–1389 (2017).ADS 
Article 

Google Scholar 
Ovejero, R. et al. Do cortisol and corticosterone play the same role in coping with stressors? Measuring glucocorticoid serum in free-ranging guanacos (Lama guanicoe). J. Exp. Zool. Part A Ecol. Genet. Physiol. 319, 539–547 (2013).CAS 
Article 

Google Scholar 
Bonacic, C., Feber, R. E. & Macdonald, D. W. Capture of the vicuña (Vicugna vicugna) for sustainable use: Animal welfare implications. Biol. Conserv. 129, 543–550 (2006).Article 

Google Scholar 
Romero, L. M. & Beattie, U. K. Common myths of glucocorticoid function in ecology and conservation. J. Exp. Zool. Part A Ecol. Integr. Physiol. 337, 7–14 (2022).CAS 
Article 

Google Scholar 
Sire, J. E. et al. The effect of blood sampling on plasma cortisol in female reindeer (Rangifer tarandus tarandus L). Acta Vet. Scand. 36, 583–587 (1995).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Harlow, H. J., Thorne, E. T., Williams, E. S., Belden, E. L. & Gern, W. A. Adrenal responsiveness in domestic sheep ( Ovis aries ) to acute and chronic stressors as predicted by remote monitoring of cardiac frequency. Can. J. Zool. 65, 2021–2027 (1987).Article 

Google Scholar 
Pottinger, T. G. & Moran, T. A. Differences in plasma cortisol and cortisone dynamics during stress in two strains of rainbow trout (Oncorhynchus mykiss). J. Fish Biol. 43, 121–130 (1993).CAS 
Article 

Google Scholar 
Arnemo, J. M. & Ranheim, B. Effects of medetomidine and atipamezole on serum glucose and cortisol levels in captive reindeer (Rangifer tarandus tarandus). Rangifer 19, 85–89 (1999).Article 

Google Scholar 
Mentaberre, G. et al. Effects of azaperone and haloperidol on the stress response of drive-net captured Iberian ibexes (Capra pyrenaica). Eur. J. Wildl. Res. 56, 757–764 (2010).Article 

Google Scholar 
Northrup, J. M., Anderson, C. R. & Wittemyer, G. Effects of helicopter capture and handling on movement behavior of mule deer. J. Wildl. Manag. 78, 731–738 (2014).Article 

Google Scholar 
Jung, T. S. et al. Short-term effect of helicopter-based capture on movements of a social ungulate. J. Wildl. Manag. 83, 830–837 (2019).Article 

Google Scholar 
Nurmi, H., Laaksonen, S., Raekallio, M. & Hänninen, L. Wintertime pharmacokinetics of intravenously and orally administered meloxicam in semi-domesticated reindeer (Rangifer tarandus tarandus). Vet. Anaesth. Analg. 49, 423–428 (2022).CAS 
PubMed 
Article 

Google Scholar 
Chapple, R. S., English, A. W., Mulley, R. C. & Lepherd, E. E. Haematology and serum biochemistry of captive unsedated chital deer (Axis axis) in Australia. J. Wildl. Dis. 27, 396–406 (1991).CAS 
PubMed 
Article 

Google Scholar 
Brosh, A. Heart rate measurements as an index of energy expenditure and energy balance in ruminants: A review1. J. Anim. Sci. 85, 1213–1227 (2007).CAS 
PubMed 
Article 

Google Scholar 
Suazo, A. A., Delong, A. T., Bard, A. A. & Oddy, D. M. Repeated capture of beach mice (Peromyscus polionotus phasma and P. P. niveiventris) reduces body mass. J. Mammal. 86, 520–523 (2005).Article 

Google Scholar 
Hoyle, S. D., Horsup, A. B., Johnson, C. N., Crossman, D. G. & McCallum, H. Live-trapping of the northern hairy-nosed wombat (Lasiorhinus krefftii): Population-size estimates and effects on individuals. Wildl. Res. 22, 741–755 (1995).Article 

Google Scholar 
Estruelas, N. F. Short- and long-term physiological effects of capture and handling on free-ranging brown bears (Ursus arctos). PhD Thesis. (Inland Norway University of Applied Sciences, 2017).Veiberg, V. et al. Maternal winter body mass and not spring phenology determine annual calf production in an Arctic herbivore. Oikos 126, 980–987 (2017).Article 

Google Scholar 
Loe, L. E. et al. The neglected season: Warmer autumns counteract harsher winters and promote population growth in Arctic reindeer. Glob. Change Biol. 27, 993–1002 (2021).ADS 
CAS 
Article 

Google Scholar 
Larsen, T. S., Nilsson, N. & Blix, A. S. Seasonal changes in lipogenesis and lipolysis in isolated adipocytes from Svalbard and Norwegian reindeer. Acta Physiol. Scand. 123, 97–104 (1985).CAS 
PubMed 
Article 

Google Scholar 
Colman, J. E., Jacobsen, B. W. & Reimers, E. Summer response distances of Svalbard reindeer (Rangifer tarandus platyrhynchus) to provocation by humans on foot. Wildlife Biol. 7, 275–283 (2001).Article 

Google Scholar 
Trondrud, L. M. et al. Determinants of heart rate in Svalbard reindeer reveal mechanisms of seasonal energy management. Philos. Trans. R. Soc. B Biol. Sci. 376, 20200215 (2021).Article 

Google Scholar 
Pigeon, G. et al. Context-dependent fitness costs of reproduction despite stable body mass costs in an Arctic herbivore. J. Anim. Ecol. 91, 61–73 (2022).PubMed 
Article 

Google Scholar 
Peeters, B., Pedersen, Å., Veiberg, V. & Hansen, B. Hunting quotas, selectivity and stochastic population dynamics challenge the management of wild reindeer. Clim. Res. https://doi.org/10.3354/cr01668 (2021).Article 

Google Scholar 
Loe, L. E. et al. Activity pattern of arctic reindeer in a predator-free environment: No need to keep a daily rhythm. Oecologia 152, 617–624 (2007).ADS 
PubMed 
Article 

Google Scholar 
Dahl, S. R. et al. Assay of steroids by liquid chromatography–tandem mass spectrometry in monitoring 21-hydroxylase deficiency. Endocr. Connect. 7, 1542–1550 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Loe, L. E. et al. Testing five hypotheses of sexual segregation in an arctic ungulate. J. Anim. Ecol. 75, 485–496 (2006).PubMed 
Article 

Google Scholar 
Reimers, E., Lund, S. & Ergon, T. Vigilance and fright behaviour in the insular Svalbard reindeer (Rangifer tarandus platyrhynchus). Can. J. Zool. 89, 753–764 (2011).Article 

Google Scholar 
The R Core Team. R: A language and environment for statistical computing (2021).Burnham, K. P. & Anderson, D. R. in Model selection and multimodel inference. A Practical Information-Theoretic Approach. Ecological Modelling (Springer, 2002).Blanchet, F. G., Tikhonov, G. & Norberg, A. HMSC: Hierarchical modelling of species community. R package version 2.2-0 (2019).Ovaskainen, O. et al. How to make more out of community data? A conceptual framework and its implementation as models and software. Ecol. Lett. 20, 561–576 (2017).PubMed 
Article 

Google Scholar 
Legendre, P. & Legendre, L. Numerical Ecology (Elsevier Science BV, 2012).MATH 

Google Scholar 
Diggle, P. J., Heagerty, P., Liang, K.-Y. & Zeger, S. L. Analysis of Longitudinal Data (Oxford University Press, 2013).MATH 

Google Scholar  More

Barnett, T. P. et al. Penetration of human-induced warming into the world’s oceans. Science 309, 284–287 (2005).CAS 
Article 

Google Scholar 
Levitus, S. et al. Global ocean heat content 1955–2008 in light of recently revealed instrumentation problems. Geophys. Res. Lett. 36, L07608 (2009).
Google Scholar 
Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Change 3, 919–925 (2013).Article 

Google Scholar 
García Molinos, J. et al. Climate velocity and the future global redistribution of marine biodiversity. Nat. Clim. Change 6, 83–88 (2016).Article 

Google Scholar 
Free, C. M. et al. Impacts of historical warming on marine fisheries production. Science 363, 979–983 (2019).CAS 
Article 

Google Scholar 
Hughes, N. F. & Grand, T. C. Physiological ecology meets the ideal-free distribution: predicting the distribution of size-structured fish populations across temperature gradients. Environ. Biol. Fishes 59, 285–298 (2000).Article 

Google Scholar 
Tittensor, D. P. et al. Global patterns and predictors of marine biodiversity across taxa. Nature 466, 1098–1101 (2010).CAS 
Article 

Google Scholar 
Sunday, J. M., Bates, A. E. & Dulvy, N. K. Global analysis of thermal tolerance and latitude in ectotherms. Proc. R. Soc. B 278, 1823–1830 (2011).Article 

Google Scholar 
Waldock, C., Stuart‐Smith, R. D., Edgar, G. J., Bird, T. J. & Bates, A. E. The shape of abundance distributions across temperature gradients in reef fishes. Ecol. Lett. 22, 685–696 (2019).Article 

Google Scholar 
Stuart-Smith, R. D., Edgar, G. J. & Bates, A. E. Thermal limits to the geographic distributions of shallow-water marine species. Nat. Ecol. Evol. 1, 1846–1852 (2017).Article 

Google Scholar 
Pinsky, M. L., Worm, B., Fogarty, M. J., Sarmiento, J. L. & Levin, S. A. Marine taxa track local climate velocities. Science 341, 1239–1242 (2013).CAS 
Article 

Google Scholar 
Beaugrand, G., Edwards, M., Raybaud, V., Goberville, E. & Kirby, R. R. Future vulnerability of marine biodiversity compared with contemporary and past changes. Nat. Clim. Change 5, 695–701 (2015).Article 

Google Scholar 
Trisos, C. H., Merow, C. & Pigot, A. L. The projected timing of abrupt ecological disruption from climate change. Nature 580, 496–501 (2020).CAS 
Article 

Google Scholar 
Levin, L. A. & Le Bris, N. The deep ocean under climate change. Science 350, 766–768 (2015).CAS 
Article 

Google Scholar 
Deutsch, C. A. et al. Impacts of climate warming on terrestrial ectotherms across latitude. Proc. Natl Acad. Sci. USA 105, 6668–6672 (2008).CAS 
Article 

Google Scholar 
Sunday, J. M., Bates, A. E. & Dulvy, N. K. Thermal tolerance and the global redistribution of animals. Nat. Clim. Change 2, 686–690 (2012).Article 

Google Scholar 
Radeloff, V. C. et al. The rise of novelty in ecosystems. Ecol. Appl. 25, 2051–2068 (2015).Article 

Google Scholar 
Lotterhos, K. E., Láruson, Á. J. & Jiang, L.-Q. Novel and disappearing climates in the global surface ocean from 1800 to 2100. Sci. Rep. 11, 15535 (2021).CAS 
Article 

Google Scholar 
Mora, C. et al. The projected timing of climate departure from recent variability. Nature 502, 183–187 (2013).CAS 
Article 

Google Scholar 
Henson, S. A. et al. Rapid emergence of climate change in environmental drivers of marine ecosystems. Nat. Commun. 8, 14682 (2017).Article 

Google Scholar 
Séférian, R. et al. Evaluation of CNRM Earth System Model, CNRM‐ESM2‐1: role of Earth system processes in present‐day and future climate. J. Adv. Model. Earth Syst. 11, 4182–4227 (2019).Article 

Google Scholar 
Gidden, M. J. et al. Global emissions pathways under different socioeconomic scenarios for use in CMIP6: a dataset of harmonized emissions trajectories through the end of the century. Geosci. Model Dev. 12, 1443–1475 (2019).CAS 
Article 

Google Scholar 
Eyring, V. et al. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev. 9, 1937–1958 (2016).Article 

Google Scholar 
Beszczynska-Möller, A., Fahrbach, E., Schauer, U. & Hansen, E. Variability in Atlantic water temperature and transport at the entrance to the Arctic Ocean, 1997–2010. ICES J. Mar. Sci. 69, 852–863 (2012).Article 

Google Scholar 
Sutton, T. T. Vertical ecology of the pelagic ocean: classical patterns and new perspectives. J. Fish. Biol. 83, 1508–1527 (2013).CAS 
Article 

Google Scholar 
Richter, I. Climate model biases in the eastern tropical oceans: causes, impacts and ways forward. WIREs Clim. Change 6, 345–358 (2015).Article 

Google Scholar 
Pozo Buil, M. et al. A dynamically downscaled ensemble of future projections for the California Current System. Front. Mar. Sci. 8, 612874 (2021).Article 

Google Scholar 
Leonard, M. et al. A compound event framework for understanding extreme impacts. WIREs Clim. Change 5, 113–128 (2014).Article 

Google Scholar 
Kwiatkowski, L. et al. Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections. Biogeosciences 17, 3439–3470 (2020).CAS 
Article 

Google Scholar 
Bopp, L. et al. Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences 10, 6225–6245 (2013).Article 

Google Scholar 
Cheng, L., Abraham, J., Hausfather, Z. & Trenberth, K. E. How fast are the oceans warming? Science 363, 128–129 (2019).CAS 
Article 

Google Scholar 
Hawkins, E. & Sutton, R. Time of emergence of climate signals. Geophys. Res. Lett. 39, L01702 (2012).Article 

Google Scholar 
Stuart-Smith, R. D., Edgar, G. J., Barrett, N. S., Kininmonth, S. J. & Bates, A. E. Thermal biases and vulnerability to warming in the world’s marine fauna. Nature 528, 88–92 (2015).CAS 
Article 

Google Scholar 
Filbee-Dexter, K. et al. Marine heatwaves and the collapse of marginal North Atlantic kelp forests. Sci. Rep. 10, 13388 (2020).CAS 
Article 

Google Scholar 
Román-Palacios, C. & Wiens, J. J. Recent responses to climate change reveal the drivers of species extinction and survival. Proc. Natl Acad. Sci. USA 117, 4211–4217 (2020).Article 
CAS 

Google Scholar 
Silvy, Y., Guilyardi, E., Sallée, J.-B. & Durack, P. J. Human-induced changes to the global ocean water masses and their time of emergence. Nat. Clim. Change 10, 1030–1036 (2020).CAS 
Article 

Google Scholar 
Cheng, L., Zheng, F. & Zhu, J. Distinctive ocean interior changes during the recent warming slowdown. Sci. Rep. 5, 14346 (2015).CAS 
Article 

Google Scholar 
Brito-Morales, I. et al. Climate velocity reveals increasing exposure of deep-ocean biodiversity to future warming. Nat. Clim. Change 10, 576–581 (2020).CAS 
Article 

Google Scholar 
Frölicher, T. L. & Laufkötter, C. Emerging risks from marine heat waves. Nat. Commun. 9, 650 (2018).Article 
CAS 

Google Scholar 
Oliver, E. C. J. et al. Marine Heatwaves. Ann. Rev. Mar. Sci. 13, 313–342 (2021).Article 

Google Scholar 
Perry, A. L., Low, P. J., Ellis, J. R. & Reynolds, J. D. Climate change and distribution shifts in marine fishes. Science 308, 1912–1915 (2005).CAS 
Article 

Google Scholar 
Chaudhary, C., Richardson, A. J., Schoeman, D. S. & Costello, M. J. Global warming is causing a more pronounced dip in marine species richness around the equator. Proc. Natl Acad. Sci. USA 118, e2015094118 (2021).CAS 
Article 

Google Scholar 
Burrows, M. T. et al. Ocean community warming responses explained by thermal affinities and temperature gradients. Nat. Clim. Change 9, 959–963 (2019).Article 

Google Scholar 
IPCC Climate Change 2022: Impacts, Adaptation, and Vulnerability (eds Pörtner, H.-O. et al.) (Cambridge Univ. Press, 2022).Cahill, A. E. et al. How does climate change cause extinction? Proc. R. Soc. B280, 20121890 (2013).Article 

Google Scholar 
Hastings, R. A. et al. Climate change drives poleward increases and equatorward declines in marine species. Curr. Biol. 30, 1572–1577.e2 (2020).CAS 
Article 

Google Scholar 
Jorda, G. et al. Ocean warming compresses the three-dimensional habitat of marine life. Nat. Ecol. Evol. 4, 109–114 (2020).Article 

Google Scholar 
Dulvy, N. K. et al. Climate change and deepening of the North Sea fish assemblage: a biotic indicator of warming seas. J. Appl. Ecol. 45, 1029–1039 (2008).Article 

Google Scholar 
Thatje, S. Climate warming affects the depth distribution of marine ectotherms. Mar. Ecol. Prog. Ser. 660, 233–240 (2021).Article 

Google Scholar 
Manuel, S. A., Coates, K. A., Kenworthy, W. J. & Fourqurean, J. W. Tropical species at the northern limit of their range: composition and distribution in Bermuda’s benthic habitats in relation to depth and light availability. Mar. Environ. Res. 89, 63–75 (2013).CAS 
Article 

Google Scholar 
Peck, L. S., Webb, K. E. & Bailey, D. M. Extreme sensitivity of biological function to temperature in Antarctic marine species. Funct. Ecol. 18, 625–630 (2004).Article 

Google Scholar 
Peck, L. S., Morley, S. A., Richard, J. & Clark, M. S. Acclimation and thermal tolerance in Antarctic marine ectotherms. J. Exp. Biol. 217, 16–22 (2014).Article 

Google Scholar 
Walsh, J. E. Climate of the Arctic marine environment. Ecol. Appl. 18, S3–S22 (2008).Article 

Google Scholar 
Storch, D., Menzel, L., Frickenhaus, S. & Pörtner, H.-O. Climate sensitivity across marine domains of life: limits to evolutionary adaptation shape species interactions. Glob. Change Biol. 20, 3059–3067 (2014).Article 

Google Scholar 
Araújo, M. B. et al. Heat freezes niche evolution. Ecol. Lett. 16, 1206–1219 (2013).Article 

Google Scholar 
Pörtner, H. O., Peck, L. & Somero, G. Thermal limits and adaptation in marine Antarctic ectotherms: an integrative view. Philos. Trans. R. Soc. B 362, 2233–2258 (2007).Article 
CAS 

Google Scholar 
Qu, Y.-F. & Wiens, J. J. Higher temperatures lower rates of physiological and niche evolution. Proc. R. Soc. B 287, 20200823 (2020).Article 

Google Scholar 
Cohen, D.M., Inada, T., Iwamoto, T. and Scialabba, N. FAO Species Catalogue, Vol. 10. Gadiform Fishes of the World (Order Gadiformes) (FAO, 1990).Strand, E. & Huse, G. Vertical migration in adult Atlantic cod (Gadus morhua). Can. J. Fish. Aquat. Sci. 64, 1747–1760 (2007).Article 

Google Scholar 
Frölicher, T. L., Fischer, E. M. & Gruber, N. Marine heatwaves under global warming. Nature 560, 360–364 (2018).Article 
CAS 

Google Scholar 
Wernberg, T. et al. Climate-driven regime shift of a temperate marine ecosystem. Science 353, 169–172 (2016).CAS 
Article 

Google Scholar 
Smale, D. A. et al. Marine heatwaves threaten global biodiversity and the provision of ecosystem services. Nat. Clim. Change 9, 306–312 (2019).Article 

Google Scholar 
Cheung, W. W. L. & Frölicher, T. L. Marine heatwaves exacerbate climate change impacts for fisheries in the northeast Pacific. Sci. Rep. 10, 6678 (2020).CAS 
Article 

Google Scholar 
Brierley, A. S. & Kingsford, M. J. Impacts of climate change on marine organisms and ecosystems. Curr. Biol. 19, R602–R614 (2009).CAS 
Article 

Google Scholar 
Bijma, J., Pörtner, H.-O., Yesson, C. & Rogers, A. D. Climate change and the oceans—what does the future hold? Mar. Pollut. Bull. 74, 495–505 (2013).CAS 
Article 

Google Scholar 
Jackson, J. B. C. et al. Historical overfishing and the recent collapse of coastal ecosystems. Science 293, 629–637 (2001).CAS 
Article 

Google Scholar 
Duarte, C. M. et al. The soundscape of the Anthropocene ocean. Science 371, eaba4658 (2021).CAS 
Article 

Google Scholar 
Rochman, C. M. & Hoellein, T. The global odyssey of plastic pollution. Science 368, 1184–1185 (2020).CAS 
Article 

Google Scholar 
Gruber, N., Boyd, P. W., Frölicher, T. L. & Vogt, M. Biogeochemical extremes and compound events in the ocean. Nature 600, 395–407 (2021).CAS 
Article 

Google Scholar 
Madec, G. et al. NEMO ocean engine. Zenodo https://www.earth-prints.org/handle/2122/13309 (2017).Mathiot, P., Jenkins, A., Harris, C. & Madec, G. Explicit representation and parametrised impacts of under ice shelf seas in the z∗- coordinate ocean model NEMO 3.6. Geosci. Model Dev. 10, 2849–2874 (2017).Article 

Google Scholar 
Dai, A. & Bloecker, C. E. Impacts of internal variability on temperature and precipitation trends in large ensemble simulations by two climate models. Clim. Dyn. 52, 289–306 (2019).Article 

Google Scholar 
Deser, C., Phillips, A., Bourdette, V. & Teng, H. Uncertainty in climate change projections: the role of internal variability. Clim. Dyn. 38, 527–546 (2012).Article 

Google Scholar 
Middag, R. et al. Intercomparison of dissolved trace elements at the Bermuda Atlantic Time Series station. Mar. Chem. 177, 476–489 (2015).CAS 
Article 

Google Scholar 
Welch, B. L. The generalization of Student’s’ problem when several different population variances are involved. Biometrika 34, 28 (1947).CAS 

Google Scholar 
Lenoir, J. et al. Species better track climate warming in the oceans than on land. Nat. Ecol. Evol. 4, 1044–1059 (2020).Article 

Google Scholar 
Janzen, D. H. Why mountain passes are higher in the Tropics. Am. Nat. 101, 233–249 (1967).Article 

Google Scholar 
Seebacher, F., White, C. R. & Franklin, C. E. Physiological plasticity increases resilience of ectothermic animals to climate change. Nat. Clim. Change 5, 61–66 (2015).Article 

Google Scholar 
Hoffmann, A. A. & Sgrò, C. M. Climate change and evolutionary adaptation. Nature 470, 479–485 (2011).CAS 
Article 

Google Scholar 
Sandblom, E. et al. Physiological constraints to climate warming in fish follow principles of plastic floors and concrete ceilings. Nat. Commun. 7, 11447 (2016).CAS 
Article 

Google Scholar 
Tewksbury, J. J., Huey, R. B. & Deutsch, C. A. Putting the heat on tropical animals. Science 320, 1296–1297 (2008).CAS 
Article 

Google Scholar 
Dahlke, F. T., Wohlrab, S., Butzin, M. & Pörtner, H.-O. Thermal bottlenecks in the life cycle define climate vulnerability of fish. Science 369, 65–70 (2020).CAS 
Article 

Google Scholar  More