Ecology

1.Müller, A. et al. The Role of Biomass in the Sustainable Development Goals: A Reality Check and Governance Implications. (Institue for Advances Sustainability Studies, 2015).2.Keesstra, S. et al. Soil-related sustainable development goals: Four concepts to make land degradation neutrality and restoration work. Land 7, 133 (2018).Article 

Google Scholar 
3.Li, X., Rubæk, G. H. & Sørensen, P. High plant availability of phosphorus and low availability of cadmium in four biomass combustion ashes. Sci. Total Environ. 557–558, 851–860 (2016).ADS 
PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
4.Hannam, K. D., Deschamps, C., Kwiaton, M., Venier, L. & Hazlett, P. W. Regulations and Guidelines for the Use of Wood Ash as a Soil Amendment in Canadian Forests (Natural Resources Canada Canadian Forest Service, 2016).5.Nieminen, M., Laiho, R., Sarkkola, S. & Penttilä, T. Whole-tree, stem-only, and stump harvesting impacts on site nutrient capital of a Norway spruce-dominated peatland forest. Eur. J. For. Res. 135, 531–538 (2016).CAS 
Article 

Google Scholar 
6.Adams, M., Burger, J., Jenkins, A. & Zelazny, L. Impact of harvesting and atmospheric pollution on nutrient depletion of eastern US hardwood forests. For. Ecol. Manag. 138, 301–319 (2000).Article 

Google Scholar 
7.Ågren, A., Buffam, I., Bishop, K. & Laudon, H. Sensitivity of pH in a boreal stream network to a potential decrease in base cations caused by forest harvest. Can. J. Fish. Aquat. Sci. 67, 1116–1125 (2010).Article 
CAS 

Google Scholar 
8.Jerabkova, L., Prescott, C. E., Titus, B. D., Hope, G. D. & Walters, M. B. A meta-analysis of the effects of clearcut and variable-retention harvesting on soil nitrogen fluxes in boreal and temperate forests. Can. J. For. Res. 41, 1852–1870 (2011).CAS 
Article 

Google Scholar 
9.Cronan, C. S. & Grigal, D. F. Use of calcium/aluminum ratios as indicators of stress in forest ecosystems. J. Environ. Qual. 24, 209–226 (1995).CAS 
Article 

Google Scholar 
10.Augusto, L., Bakker, M. R. & Meredieu, C. Wood ash applications to temperate forest ecosystems—Potential benefits and drawbacks. Plant Soil 306, 181–198 (2008).CAS 
Article 

Google Scholar 
11.Pitman, R. M. Wood ash use in forestry—A review of the environmental impacts. Forestry 79, 563–588 (2006).Article 

Google Scholar 
12.Pugliese, S. et al. Wood ash as a forest soil amendment: The role of boiler and soil type on soil property response. Can. J. Soil Sci. 94, 621–634 (2014).CAS 
Article 

Google Scholar 
13.Saarsalmi, A., Mälkönen, E. & Kukkola, M. Effect of wood ash fertilization on soil chemical properties and stand nutrient status and growth of some coniferous stands in Finland. Scand. J. For. Res. 19, 217–233 (2004).Article 

Google Scholar 
14.Zimmermann, S. & Frey, B. Soil respiration and microbial properties in an acid forest soil: Effects of wood ash. Soil Biol. Biochem. 34, 1727–1737 (2002).CAS 
Article 

Google Scholar 
15.Perkiömäki, J. & Fritze, H. Short and long-term effects of wood ash on the boreal forest humus microbial community. Soil Biol. Biochem. 34, 1343–1353 (2002).Article 

Google Scholar 
16.Etiégni, L. & Campbell, A. G. Physical and chemical characteristics of wood ash. Bioresour. Technol. 37, 173–178 (1991).Article 

Google Scholar 
17.Saarsalmi, A., Smolander, A., Kukkola, M. & Arola, M. Effect of wood ash and nitrogen fertilization on soil chemical properties, soil microbial processes, and stand growth in two coniferous stands in Finland. Plant Soil 331, 329–340 (2010).CAS 
Article 

Google Scholar 
18.Saarsalmi, A. & Levula, T. Wood ash application and liming: Effects on soil chemical properties and growth of Scots pine transplants. Balt. For. 13, 149–157 (2007).
Google Scholar 
19.Demeyer, A., Voundi Nkana, J. & Verloo, M. Characteristics of wood ash and influence on soil properties and nutrient uptake: An overview. Bioresour. Technol. 77, 287–295 (2001).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
20.Zhou, H. Reducing, reusing and recycling solid wastes from wood fibre processing. In Towards Sustainable Management of the Boreal Forest (eds Burton, P. J. et al.) 759–798 (NRC Research Press, 2003) https://doi.org/10.1139/9780660187624.Chapter 

Google Scholar 
21.Rosenberg, O., Persson, T., Högbom, L. & Jacobson, S. Effects of wood-ash application on potential carbon and nitrogen mineralisation at two forest sites with different tree species, climate and N status. For. Ecol. Manag. 260, 511–518 (2010).Article 

Google Scholar 
22.Saarsalmi, A., Smolander, A., Kukkola, M., Moilanen, M. & Saramäki, J. 30-Year effects of wood ash and nitrogen fertilization on soil chemical properties, soil microbial processes and stand growth in a Scots pine stand. For. Ecol. Manag. 278, 63–70 (2012).Article 

Google Scholar 
23.McDonald, M. A., Hawkins, B. J., Prescott, C. E. & Kimmins, J. P. Growth and foliar nutrition of western red cedar fertilized with sewage sludge, pulp sludge, fish silage, and wood ash on northern Vancouver Island. Can. J. For. Res. 24, 297–301 (1994).Article 

Google Scholar 
24.Park, B. B., Yanai, R. D., Sahm, J. M., Lee, D. K. & Abrahamson, L. P. Wood ash effects on plant and soil in a willow bioenergy plantation. Biomass Bioenergy. 28, 355–365 (2005).CAS 
Article 

Google Scholar 
25.Batjes, N. H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 47, 151–163 (1996).CAS 
Article 

Google Scholar 
26.Nohrstedt, H. -Ö. Response of coniferous forest ecosystems on mineral soils to nutrient additions: A review of Swedish experiences. Scand. J. For. Res. 16, 555–573 (2001).Article 

Google Scholar 
27.Rennenberg, H. & Dannenmann, M. Nitrogen nutrition of trees in temperate forests—The significance of nitrogen availability in the pedosphere and atmosphere. Forests 6, 2820–2835 (2015).Article 

Google Scholar 
28.Solla-Gullón, F., Santalla, M., Pérez-Cruzado, C., Merino, A. & Rodríguez-Soalleiro, R. Response of Pinus radiata seedlings to application of mixed wood-bark ash at planting in a temperate region: Nutrition and growth. For. Ecol. Manag. 255, 3873–3884 (2008).Article 

Google Scholar 
29.Saarsalmi, A., Kukkola, M., Moilanen, M. & Arola, M. Long-term effects of ash and N fertilization on stand growth, tree nutrient status and soil chemistry in a Scots pine stand. For. Ecol. Manag. 235, 116–128 (2006).Article 

Google Scholar 
30.Ohno, T. Neutralization of soil acidity and release of phosphorus and potassium by wood ash. J. Environ. Qual. 21, 433–438 (1992).Article 

Google Scholar 
31.Bang-Andreasen, T. et al. Wood ash induced pH changes strongly affect soil bacterial numbers and community composition. Front. Microbiol. 8, 1400 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
32.Clapham, W. M. & Zibilske, L. M. Wood ash as a liming amendment. Commun. Soil Sci. Plant Anal. 23, 1209–1227 (1992).CAS 
Article 

Google Scholar 
33.Ulery, A. L., Graham, R. C. & Amrhein, C. Wood-ash composition and soil pH following intense burning. Soil Sci. 156, 358–364 (1993).ADS 
CAS 
Article 

Google Scholar 
34.Muse, J. K. & Mitchell, C. C. Paper mill boiler ash and lime by-products as soil liming materials. Agron. J. 87, 432–438 (1995).Article 

Google Scholar 
35.Bramryd, T. & Fransman, B. Silvicultural use of wood ashes – Effects on the nutrient and heavy metal balance in a pine (Pinus sylvestris, L) forest soil. Water Air Soil Pollut. 85, 1039–1044 (1995).ADS 
CAS 
Article 

Google Scholar 
36.Jones, D. L. & Quilliam, R. S. Metal contaminated biochar and wood ash negatively affect plant growth and soil quality after land application. J. Hazard. Mater. 276, 362–370 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
37.Kim, K.-D. Species alterations caused by nitrogen and carbon addition in nutrient-deficient municipal waste landfills. J. Ecol. Environ. 30, 161–170 (2007).Article 

Google Scholar 
38.Jacobson, S. Addition of stabilized wood ashes to Swedish coniferous stands on mineral soils—Effects on stem growth and needle nutrient concentrations. Silva Fenn. 37, 437–450 (2003).Article 

Google Scholar 
39.Hånell, B. & Magnusson, T. An evaluation of land suitability for forest fertilization with biofuel ash on organic soils in Sweden. For. Ecol. Manag. 209, 43–55 (2005).Article 

Google Scholar 
40.Hannam, K. D. et al. Wood ash as a soil amendment in Canadian forests: What are the barriers to utilization?. Can. J. For. Res. 48, 442–450 (2018).Article 

Google Scholar 
41.Woods, S. W. & Balfour, V. N. The effect of ash on runoff and erosion after a severe forest wildfire, Montana, USA. Int. J. Wildl. Fire 17, 535–548 (2008).Article 

Google Scholar 
42.Someshwar, A. Wood and combination wood-fired boiler ash characterization. J. Environ. Qual. 25, 962–972 (1996).CAS 
Article 

Google Scholar 
43.National Institute of Agricultural Science and Technology. Methods of Soil and Plant Analysis (ed. Im, J.-N.) 202 (National Institute of Agricultural Science and Technology, RDA, 2000).44.Sparks, D. L., Page, A. L., Helmke, P. A. & Leoppert, R. H. Methods of Soil Analysis Part3—Chemical Methods, SSSA Book Ser. 5.3. (Soil Science Society of America, American Society of Agronomy, 1996) https://doi.org/10.2136/sssabookser5.3.45.Walkley, A. & Black, I. A. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 37, 29–38 (1934).ADS 
CAS 
Article 

Google Scholar 
46.Nelson, D. W. & Sommers, L. E. Total carbon, organic carbon, and organic matter. In Methods of Soil Analysis Part3—Chemical Methods, SSSA Book Ser. 5.3 (eds Sparks, D. L. et al.) 961–1010 (Soil Science Society of America, American Society of Agronomy, 1996) https://doi.org/10.2136/sssabookser5.3.c34.Chapter 

Google Scholar 
47.Bremner, J. M. Nitrogen-total. In Methods of Soil Analysis Part 3—Chemical Methods, SSSA Book Ser. 5.3 (eds Sparks, D. L. et al.) 1085–1121 (Soil Science Society of America, American Society of Agronomy, 1996) https://doi.org/10.2136/sssabookser5.3.c37.Chapter 

Google Scholar 
48.Kuo, S. Phosphorus. In Methods of Soil Analysis Part 3—Chemical Methods, SSSA Book Ser. 5.3 (eds Sparks, D. L. et al.) 869–919 (Soil Science Society of America, American Society of Agronomy, 1996) https://doi.org/10.2136/sssabookser5.3.c32.Chapter 

Google Scholar 
49.Cox, M. S. The Lancaster soil test method as an alternative to the Mehlich 3 soil test method1. Soil Sci. 166, 484–489 (2001).ADS 
CAS 
Article 

Google Scholar 
50.Sumner, M. E. & Miller, W. P. Cation exchange capacity and exchange coefficients. In Methods of Soil Analysis Part 3—Chemical Methods, SSSA Book Ser. 5.3 (eds Sparks, D. L. et al.) 1201–1229 (Soil Science Society of America, American Society of Agronomy, 1996).
Google Scholar  More

Flow velocity thresholdThere were 92 Chinese sturgeon signals from 2016 to 2019, which were identified with the DIDSON dual-frequency video sonar system. The distribution map of Chinese sturgeon signals was shown in Fig. 1. The number of monitored signals in 2016 was significantly higher than in 2017–2019. The latest wild reproduction of the Chinese sturgeon occurred in 2016. Overall, most Chinese sturgeon signals were distributed within 500 m downstream from Gezhouba Dam, and there were more in the right side(facing downstream) than in the left side. The flow field of each sturgeon signal was simulated by the model, and the velocity of each signal location was obtained. According to the statistical analysis of the flow velocity values, the frequency of the sturgeon signal at different flow velocity values was shown in Fig. 2. The results show that most signals were concentrated in areas with flow velocities of 0.6–1.5 m/s, which accounted for 88.1% of the signals; areas with flow velocities below 0.6 m/s accounted for 4.3% of the signals, and areas with flow velocities above 1.5 m/s accounted for 7.6%. Therefore, 0.6–1.5 m/s was selected as the preferred flow velocity range of the Chinese sturgeon for spawning activity. This result was approximately consistent with the ranges proposed by most other researchers. The low limit of the velocity range was lower than that of other researchers. There may be two reasons for this result: the first was that the bottom velocity we analysed was lower than the surface velocity and vertical average velocity under the same conditions; the second was that our research time was after 2016, and the discharge during the spawning period was relatively low, so the velocity of the Chinese sturgeon signal was also relatively low.Figure 1Distribution map of Chinese sturgeon signals, where ○ indicates Chinese sturgeon signals monitored in 2016, ∆ indicates those in 2017, □ indicates those in 2018, and ✩ indicates those in 2019. Map generated in ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/overview).Full size imageFigure 2Plots of the frequency for the different flow velocity ranges of Chinese sturgeon signals.Full size imageDifferent opening modes with identical dischargeThe discharge of 6150 m3/s on November 24, 2016, when the latest wild reproduction of Chinese sturgeon occurred, was used to study the flow velocity distribution with different opening modes. The specific opening mode cases are shown in Table 1. Case 1 was the actual situation, and the Dajiang Plant featured 7 open units: #8, #11, #13, #14, #16, #19, and #21. According to the amounts of electricity generated by Dajiang Plant and Erjiang Plant on that day, the proportion of the Dajiang River flow was 58.8%, and the average discharge of each unit was 516.6 m3/s. Case 2 and case 3 featured 7 open units with the same discharge, but in case 2, units #15–21 were continuously open on the right-side (facing downstream), and in case 3, units #8–14 were continuously open near the left side. Case 4 and case 5 were the most concentrated conditions with the discharge of 6150 m3/s because the maximum through-discharge for each unit in the Dajiang Plant is 825 m3/s19. In these cases, at least 5 units were open with an average discharge of 723 m3/s per unit. Case 4 involved continuously opening units #8–12 on the left side, and case 5 involved continuously opening units #17–21 on the right side. Case 6 involved simultaneously opening 14 units on Dajiang River, and the average discharge of each unit was 258.3 m3/s.Table 1 Calculation cases with different opening modes of units under the identical discharge.Full size tableFigure 3 shows the flow fields of the spawning ground under different opening modes with identical discharge. By comparing the areas with a velocity threshold range of 0.6–1.5 m/s in different cases, the most favourable opening mode was determined. In case 1, the velocity at the outlet of the units was higher than the 1.5 m/s velocity threshold, but the discharge of each unit was only 516.6 m3/s, so the high-velocity range was limited, and most areas were suitable. In case 2 and case 3, there was a large difference in proportions of suitable area. Because the left side was deeper than the right side, the flow velocity on the right side was higher under the same discharge, and case 3 more easily exceeded the flow threshold, which resulted in a larger unsuitable area. Case 2 was more suitable than case 1, which also demonstrated that opening the left-side units was more favourable. In case 4 and case 5, the proportions of suitable area were small. Because the units were concentrated, the discharge of each unit was too high, and the outlet velocity was more than 2 m/s, so a large area of high velocity appeared downstream of the units with backflow under the shut-down units. The proportion of suitable area in case 5 was larger than those in case 4 and case 3, which further indicates that opening the left-side units was more favourable than opening the right-side units. Case 6 was greater than that of any other case. Because the discharge of each unit was only 258.3 m3/s, the velocity of the unit outlet was less than 1.5 m/s, and almost all areas were suitable except for the small areas on both sides. The suitable-velocity area was the largest when all units of the Dajiang Plant of Gezhouba Dam were open; therefore, for a given discharge, it was best to open all units.Figure 3Flow field of the spawning ground in different opening modes with identical discharge, where the numbers at the top of each picture are the numbers of units to open, and the arrows indicate the direction of the water flow. Maps generated in Tecplot360 EX 2020 R1 (https://www.tecplot.com/products/tecplot-360/).Full size imageDifferent discharges under identical opening modeThe velocity distribution of the spawning field is affected by the opening mode of the units and discharge of Gezhouba Dam. To study the effect of different discharges, 14 cases were simulated, as shown in Table 2. All units of the Dajiang Plant were considered open because the proportion of suitable area was expected to be maximal under such circumstances. From 1982 to the present, the discharge during the spawning day of Chinese sturgeon under Gezhouba Dam has a wide range: the highest discharge was 27,290 m3/s in 1990, and the lowest discharge was 5590 m3/s in 2012. However, the highest design discharge of the Gezhouba Dam units is 17,930 m3/s20. Once the design discharge is exceeded, the spillway on Erjiang River discharges water, and the velocity distribution of the study area is not affected. Therefore, case 1 represents the lowest discharge of 5590 m3/s, and case 2 represents a discharge of 6000 m3/s. For each subsequent case, the discharge was increased by 1000 m3/s to case 13 with the highest flow of 17,930 m 3/s. In case 14, all units reached the design discharge, and the discharge of each unit was 825 m3/s19.Table 2 Calculation cases with the same opening mode under different discharges.Full size tableFigure 4 shows the proportion of suitable-velocity area with all units open under different discharges. According to the calculation results, the proportion of suitable area slightly fluctuated at approximately 96.2% for discharges of 5590–11,000 m3/s. Because the discharge of each unit was low, the velocity of the unit outlet was low, and most areas were within the velocity threshold. Therefore, it is advantageous to open all units when the discharge is low. After the discharge reached 12,000 m3/s, the proportion of suitable area rapidly decreased. Because the discharge of each unit was high, on the right side of Dajiang River, the velocity of the unit outlet exceeded the velocity threshold and increased with increases in discharge, and the range of effect gradually increased. In the last case, the proportion of suitable area was only 6% when the units reached the designed discharge of 825 m3/s. Because the discharge of each unit was too high, almost all areas exceeded the velocity threshold except for small areas on both sides. Therefore, at discharges below 12,000 m3/s, opening all units is favourable, and at discharge above 12,000 m3/s, a higher discharge corresponds to more unfavourable conditions.Figure 4Proportions of the suitable-velocity area with all units opened under different discharges.Full size imageOptimal scheme under high-flow conditionsHigh-flow conditions at Gezhouba Dam are considered those that exceed 12,000 m3/s because of the substantive decline in suitable habitat area at higher discharges. Because opening the units on the left side of the Dajiang Plant provides a more uniform, suitable habitat, we evaluated 20 cases with a left-side opening mode under different discharge, as shown in Table 3. Because the highest discharge of each unit in the Dajiang Plant is 825 m3/s, at least 9 units must be open when the discharge is 12,000 m3/s. Case 1 was designed to open 9 units on the left, i.e., units #13–21, and the discharge of each unit was 784 m3/s. Cases 2–5 increased by 1 unit from left to right until 13 units were opened. For discharges of 13,000 m3/s, 14,000 m3/s, 15,000 m3/s, and 16,000 m3/s, at least 10, 10, 11, and 12 units were opened. When the discharge was 17,000 m3/s and 17,930 m3/s, at least 13 units were open.Table 3 Calculation cases with different opening modes under high-flow conditions.Full size tableFigure 5 shows the proportions of suitable area for different opening modes under high-flow conditions. The calculation results show that when the discharge was 12,000 m3/s, 13,000 m3/s, and 14,000 m3/s, the proportion of suitable area showed a parabolic trend with the increase in number of units. When the discharge was 12,000 m3/s, the proportion of suitable area with 11 open units on the left was the largest, which was 8.7% larger than the value for all open units and 15% larger than the value for the lowest number of open units. When the discharge was 13,000 m3/s, 12 open units on the left had the largest proportion of suitable-flow-velocity area. When the discharge was 14,000 m3/s, the proportions of suitable area produced by opening 12 and 13 units on the left were the largest. The proportion of suitable area of the lowest number of open units was usually minimal because the discharge of each unit was too high, which resulted in a large area of high velocity that was not suitable for Chinese sturgeon to spawn. Because of the underwater topography, opening the left-side units was more favourable than opening the right-side units, so for all open units, the proportions of suitable area will be lower, and the number of units opened in the middle will be the most advantageous. For a discharge of 15,000 m3/s, with the increase in number of units, the proportion of suitable area increased, and there was no parabolic trend because the discharge of each unit exceeded 678 m3/s; thus, on the left side, there was a large area of high velocity, and the effect extended very far, which was not suitable for Chinese sturgeon.Figure 5Proportions of the suitable area for different opening modes under high-flow conditions, where 12,000–09 on the x-axis indicates that the discharge is 12,000 m3/s, and 9 units are open on the left.Full size image More

1.Maxwell, S. L., Fuller, R. A., Brooks, T. M. & Watson, J. E. M. Biodiversity: The ravages of guns, nets and bulldozers. Nature 536, 143–145 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
2.Ripple, W. J., Wolf, C., Newsome, T. M., Barnard, P. & Moomaw, W. R. World scientists’ warning of a climate emergency. Bioscience https://doi.org/10.1093/biosci/biz088 (2019).Article 

Google Scholar 
3.Seneviratne, S. I., Lüthi, D., Litschi, M. & Schär, C. Land–atmosphere coupling and climate change in Europe. Nature 443, 205–209 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
4.Liu, H. et al. Shifting plant species composition in response to climate change stabilizes grassland primary production. Proc. Natl. Acad. Sci. 115, 4051–4056 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
5.Schweiger, O., Settele, J., Kudrna, O., Klotz, S. & Kühn, I. Climate change can cause spatial mismatch of trophically interacting species. Ecology 89, 3472–3479 (2008).PubMed 
Article 

Google Scholar 
6.Parmesan, C. et al. Poleward shifts in geographical ranges of butterfly species associated with regional warming. Nature 399, 579–583 (1999).ADS 
CAS 
Article 

Google Scholar 
7.Dieker, P., Drees, C. & Assmann, T. Two high-mountain burnet moth species (Lepidoptera, Zygaenidae) react differently to the global change drivers climate and land-use. Biol. Conserv. 144, 2810–2818 (2011).Article 

Google Scholar 
8.Habel, J. C., Rödder, D., Schmitt, T. & Nève, G. Global warming will affect the genetic diversity and uniqueness of Lycaena helle populations. Glob. Change Biol. 17, 194–205 (2011).ADS 
Article 

Google Scholar 
9.Grabherr, G., Gottfried, M. & Pauli, H. Climate change impacts in alpine environments. Geogr. Compass 4, 1133–1153 (2010).Article 

Google Scholar 
10.Alexander, J. M. et al. Lags in the response of mountain plant communities to climate change. Glob. Change Biol. 24, 563–579 (2018).11.Renner, S. S. & Zohner, C. M. Climate change and phenological mismatch in trophic interactions among plants, insects, and vertebrates. Annu. Rev. Ecol. Evol. Syst. 49, 165–182 (2018).Article 

Google Scholar 
12.Fleishman, E. & Murphy, D. D. A realistic assessment of the indicator potential of butterflies and other charismatic taxonomic groups. Conserv. Biol. 23, 1109–1116 (2009).PubMed 
Article 

Google Scholar 
13.Sexton, J. P., Montiel, J., Shay, J. E., Stephens, M. R. & Slatyer, R. A. Evolution of ecological niche breadth. Annu. Rev. Ecol. Evol. Syst. 48, 183–206 (2017).Article 

Google Scholar 
14.Herrera, J. M., Ploquin, E. F., Rasmont, P. & Obeso, J. R. Climatic niche breadth determines the response of bumblebees (Bombus spp.) to climate warming in mountain areas of the Northern Iberian Peninsula. J. Insect Conserv. 22, 771–779 (2018).Article 

Google Scholar 
15.Habel, J. C. et al. Butterfly community shifts over two centuries. Conserv. Biol. 30, 754–762 (2016).PubMed 
Article 

Google Scholar 
16.Descombes, P., Pradervand, J. N., Golay, J., Guisan, A. & Pellissier, L. Simulated shifts in trophic niche breadth modulate range loss of alpine butterflies under climate change. Ecography 39, 796–804 (2016).Article 

Google Scholar 
17.Kerr, J. T. Racing against change: Understanding dispersal and persistence to improve species’ conservation prospects. Proc. R. Soc. B 287, 20202061 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
18.Dapporto, L., Cini, A., Voda, R., Dinca, V., Wiemers, M., Menchetti, M., Magini, G., Talavera, G., Shreeve, T., Bonelli, S., Casacci, L. P., Balletto, E., Scalercio, S. & Vila, R. Data from: Integrating three comprehensive datasets shows that mitochondrial DNA variation is linked to species traits and paleogeographic events in European butterflies. (Version 2, p. 4647103 bytes). Dryad (2019).19.Wiemers, M. et al. An updated checklist of the European butterflies (Lepidoptera, Papilionoidea). ZooKeys 811, 9–45 (2018).Article 

Google Scholar 
20.Wiemers, M., Chazot, N., Wheat, C., Schweiger, O. & Wahlberg, N. A complete time-calibrated multi-gene phylogeny of the European butterflies. ZooKeys 938, 97–124 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
21.Middleton-Welling, J. et al. A new comprehensive trait database of European and Maghreb butterflies, Papilionoidea. Sci. Data 7, 351 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
22.Weckström, K. et al. Impacts of climate warming on alpine lake biota over the past decade. Arct. Antarct. Alp. Res. 48, 361–376 (2016).Article 

Google Scholar 
23.Steinbauer, K., Lamprecht, A., Winkler, M., Bardy-Curchhalter, M., Kreiner, D., Suen, M. & Pauli, H. Shifting composition and functioning in alpine plant communities—Evidence of climate warming effects from 14 years biodiversity observation in the Northeastern Alps. In Conference Vol. 621–622 (2017).24.Bräu, M., Arbeitsgemeinschaft Bayerischer Entomologen & Bayerisches Landesamt für Umwelt (Eds.). Tagfalter in Bayern: 26 Tabellen. (Ulmer, 2013).25.Weidemann, H.-J. Tagfalter Vol. 1 (Neumann-Neudamm, 1986).
Google Scholar 
26.Weidemann, H.-J. Tagfalter: Biologie-Ökologie-Biotopschutz Vol. 2 (Neumann-Neudamm, 1988).
Google Scholar 
27.Konvicka, M., Maradova, M., Benes, J., Fric, Z. & Kepka, P. Uphill shifts in distribution of butterflies in the Czech Republic: Effects of changing climate detected on a regional scale. Glob. Ecol. Biogeogr. 12, 403–410 (2003).Article 

Google Scholar 
28.Wilson, R. J., Gutiérrez, D., Gutiérrez, J. & Monserrat, V. J. An elevational shift in butterfly species richness and composition accompanying recent climate change. Glob. Change Biol. 13, 1873–1887 (2007).ADS 
Article 

Google Scholar 
29.Wilson, R. J. et al. Changes to the elevational limits and extent of species ranges associated with climate change: Elevational shifts accompany climate change. Ecol. Lett. 8, 1138–1146 (2005).PubMed 
Article 
PubMed Central 

Google Scholar 
30.Forister, M. L. et al. Compounded effects of climate change and habitat alteration shift patterns of butterfly diversity. Proc. Natl. Acad. Sci. 107, 2088–2092 (2010).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
31.Warren, M. S. et al. Rapid responses of British butterflies to opposing forces of climate and habitat change. Nature 414, 65–69 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
32.Hill, J. K. et al. Responses of butterflies to twentieth century climate warming: Implications for future ranges. Proc. R. Soc. Lond. Ser. B Biol. Sci. 269, 2163–2171 (2002).CAS 
Article 

Google Scholar 
33.Essens, T., van Langevelde, F., Vos, R. A., Van Swaay, C. A. M. & WallisDeVries, M. F. Ecological determinants of butterfly vulnerability across the European continent. J. Insect Conserv. 21, 439–450 (2017).Article 

Google Scholar 
34.van Swaay, C., Warren, M. & Loïs, G. Biotope use and trends of European butterflies. J. Insect Conserv. 10, 189–209 (2006).Article 

Google Scholar 
35.Pyke, G. H., Thomson, J. D., Inouye, D. W. & Miller, T. J. Effects of climate change on phenologies and distributions of bumble bees and the plants they visit. Ecosphere 7, e01267 (2016).Article 

Google Scholar 
36.Biella, P. et al. Distribution patterns of the cold adapted bumblebee Bombus alpinus in the Alps and hints of an uphill shift (Insecta: Hymenoptera: Apidae). J. Insect Conserv. 21, 357–366 (2017).Article 

Google Scholar 
37.Parolo, G. & Rossi, G. Upward migration of vascular plants following a climate warming trend in the Alps. Basic Appl. Ecol. 9, 100–107 (2008).Article 

Google Scholar 
38.Filazzola, A., Matter, S. F. & Roland, J. Inclusion of trophic interactions increases the vulnerability of an alpine butterfly species to climate change. Glob. Change Biol. 26, 2867–2877 (2020).ADS 
Article 

Google Scholar 
39.Schweiger, O. et al. Multiple stressors on biotic interactions: How climate change and alien species interact to affect pollination. Biol. Rev. 85, 777–795 (2010).PubMed 

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

Google Scholar 
41.Birkhofer, K. et al. Land-use type and intensity differentially filter traits in above- and below-ground arthropod communities. J. Anim. Ecol. 86, 511–520 (2017).PubMed 
Article 

Google Scholar 
42.Dapporto, L. & Dennis, R. L. H. The generalist–specialist continuum: Testing predictions for distribution and trends in British butterflies. Biol. Conserv. 157, 229–236 (2013).Article 

Google Scholar 
43.Bartoňová, A., Benes, J. & Konvicka, M. Generalist–specialist continuum and life history traits of Central European butterflies (Lepidoptera)—Are we missing a part of the picture?. Eur. J. Entomol. 111, 543–553 (2014).Article 

Google Scholar 
44.Bartoňová, A. et al. Isolated Asian steppe element in the Balkans: Habitats of Proterebia afra (Lepidoptera: Nymphalidae: Satyrinae) and associated butterfly communities. J. Insect Conserv. 21, 559–571 (2017).Article 

Google Scholar 
45.Hodkinson, I. D. Terrestrial insects along elevation gradients: Species and community responses to altitude. Biol. Rev. 80, 489 (2005).PubMed 
Article 

Google Scholar 
46.Roth, T., Plattner, M. & Amrhein, V. Plants, birds and butterflies: Short-term responses of species communities to climate warming vary by taxon and with altitude. PLoS ONE 9, e82490 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
47.Biesmeijer, J. C. et al. Parallel declines in pollinators and insect-pollinated plants in Britain and the Netherlands. Science 313, 351–354 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
48.Filz, K. J., Engler, J. O., Stoffels, J., Weitzel, M. & Schmitt, T. Missing the target? A critical view on butterfly conservation efforts on calcareous grasslands in south-western Germany. Biodivers. Conserv. 22, 2223–2241 (2013).Article 

Google Scholar 
49.Hiebl, J. & Frei, C. Daily temperature grids for Austria since 1961—Concept, creation and applicability. Theor. Appl. Climatol. 124, 161–178 (2016).ADS 
Article 

Google Scholar 
50.Hiebl, J. & Frei, C. Daily precipitation grids for Austria since 1961—Development and evaluation of a spatial dataset for hydroclimatic monitoring and modelling. Theor. Appl. Climatol. 132, 327–345 (2018).ADS 
Article 

Google Scholar 
51.Bivand, R. & Yu, D. spgwr: Geographically Weighted Regression (R Package Version 0.6-34) [Computer Software]. https://CRAN.R-project.org/package=spgwr (2019).52.Hijmans, R. J. raster: Geographic Data Analysis and Modeling (R Package Version 3.3-13) [Computer Software]. https://CRAN.R-project.org/package=raster (2019).53.Hijmans, R. J., Phillips, S., Leathwick, J. & Elith, J. dismo: Species Distribution Modeling (R Package Version 1.1-4) [Computer Software]. https://CRAN.R-project.org/package=dismo (2017)54.Höttinger, H. & Pennerstorfer, J. Rote Liste der Tagschmetterlinge Österreichs (Lepidoptera: Papilionoidea & Hesperioidea). In Rote Listen gefährdeter Tiere Österreichs. Checklisten, Gefährdungsanalysen, Handlungsbedarf. Teil 1: Säugetiere, Vögel, Heuschrecken, Wasserkäfer, Netzflügler, Schnabelfliegen, Tagfalter. Grüne Reihe des Bundesministeriums für Land- und Forstwirtschaft, Umwelt und Wasserwirtschaft (Gesamtherausgeberin Ruth Wallner) Band 14/1 (ed. Zulka, K. P.) 313–354 (Böhlau, 2005).55.Blonder, B. & Harris, D. J. hypervolume: High Dimensional Geometry and Set Operations Using Kernel Density Estimation, Support Vector Machines, and Convex Hulls (R Package Version 2.0.12) [Computer Software]. https://CRAN.R-project.org/package=hypervolume (2019).56.Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259 (2006).Article 

Google Scholar 
57.Phillips, S. J., Anderson, R. P., Dudík, M., Schapire, R. E. & Blair, M. E. Opening the black box: An open-source release of Maxent. Ecography 40, 887–893 (2017).Article 

Google Scholar 
58.Phillips, S. J., Dudík, M. & Schapire, R. E. Maxent Software for Modeling Species Niches and Distributions (Version 3.4.1) [Computer Software]. http://biodiversityinformatics.amnh.org/open_source/maxent/ (2017).59.Swets, J. Measuring the accuracy of diagnostic systems. Science 240, 1285–1293 (1988).ADS 
MathSciNet 
CAS 
PubMed 
MATH 
Article 

Google Scholar 
60.Weiss, M. & Banko, G. Ecosystem Type Map v3.1—Terrestrial and Marine Ecosystems. ETC/BD report to the EEA (2018). More

The gobi region ecosystem has low stability because of its single species composition and simple structure (Fig. 8a). Large-scale shrub planting destroyed the original stable state (Fig. 8b) and resulted in another stable state via self-adjustment. In this process, the planted shrubs deteriorated the original ecosystem by competing for water and a chain reaction may ensue, leading to greater ecological problems. The original intention of the large-scale planting of shrubs was to maintain regional ecological balance, protect biodiversity, and fix sand, thus improving the environment (Fig. 8c). However, given the poor choice of the planting location, the expected results were not achieved. In fact, the opposite results of the original good intentions were achieved (Fig. 8d).Figure 8Diagram of different development stages of large-scale afforestation in the gobi region (a: original ground surface; b: holes dug for afforestation; c the living trees planted; d: ground surface when the trees are dead).Full size imageChina has a large expanse of arid areas, and has suffered from droughts for a long time. Land afforestation has been at the forefront of China’s policy principles, and there are government departments specializing in this field. In recent years, the Chinese Government has recommended a series of major strategies, for example, the “construction of ecological civilization” and “lucid waters and lush mountains are invaluable assets”, and also promoted greening projects, including “Three North Shelterbelt Project”, “Beijing-Tianjin Sandstorm Source Control Project”, and the “Natural Forest Protection Project”. More recently, desert greening has been conducted by people and enterprises, for example, the Ant Forest and Society of Entrepreneurs & Ecology (SEE). As a result of these projects and initiatives, China’s greening has contributed to global greening totals15,16. For afforestation, China’s policy departments have recommended the principles of “sticking to local conditions, suitable land for green, suitable trees for trees, suitable shrub for shrub, suitable grass for grass” and promoting the overall protection of “Mountain-River-Forest-Farmland-Lake-Grass-Desert system”, with particular references to desert. Their goal is to scientifically promote afforestation of the land and to clarify “where to afforest, what to afforest, how to afforest, how to manage”. However, problems arise very easily when grassroots executors are involved.The total area of the gobi region in China is approximately 56.95 × 104 km2, accounting for 13.36% of the national area, and is primarily distributed in the northwest extreme arid regions17. As mentioned above, gobi refers to a special arid landform that has a notably low water supply and is unsuitable for growing trees and shrubs. As an important natural landform, the gobi plays a key role in ecological protection; hence, its reference as “black vegetation”. However, there is a lack of understanding of the gobi, and it is often regarded as an area that needs to be greened or reformed. However, gobi, as an extremely arid region, is a fragile ecosystem. Once the gravel on the gobi surface is destroyed, it could lead to a series of ecological and environmental problems. Therefore, afforestation in arid areas is both a scientific and technical issue which must be conducted according to different regional characteristics, rather than by blindly planting trees in unsuitable areas. This study aims to attract more attention from the government forestry department and implementation personnel involved in afforestation activities so as to revise relevant policies. In response to the findings of this study, we have several recommendations: (1) it is necessary to popularize the understanding of scientific greening within the general public; (2) scientific understanding of the gobi needs to be increased, and awareness must be raised to promote its protection; (3) afforestation projects and management must be scientifically and systematically improved to ensure long-term effectiveness, and; (4) restoration and protection measures should be taken immediately in the gobi regions that have been afforested or destroyed.One of the most important causes of all these problems is the implementation of national policies on subsidies for greening and planting trees in desert areas. According to our survey, personnel who specifically plant trees and engage in afforestation are businessmen, farmers, or others, with most of them being businessmen from abroad, and only a few being local people. All the personnel are more concerned about the subsidies than greening and planting trees itself. According to the policy, they will receive majority of the subsidy if the planted trees live for three years, irrespective of whether the trees survive after that. Therefore, to guarantee the survival of the planted trees for three years, they even use water tankers to carry water to the trees from a great distance. However, after three years, the people stop watering the trees planted in the Gobi region, thereby leading to the death of trees after a few years as they cannot survive only on natural precipitation and groundwater. In pursuit of maximum profits, these businessmen will pursue larger areas for planting trees, which will cause further damage to the ecological environment in the Gobi region. Based on the current situation, we propose the following suggestions: (1) Trees that are planted must be monitored over a long time period, which will greatly reduce the short-term profit motive of the people engaged in planting trees. (2) We must plan greening and planting trees according to local conditions, respecting the laws of nature. Not all areas should be greened; moreover, we should not plant trees, especially in the gobi region, where planting trees can possibly destroy the gobi ecological environment, which is a very fragile desert ecosystem. (4) Personnel responsible for the destruction of the gobi ecological environment by unscientific greening and planting of trees must be obligated to restore the surface conditions of the gobi to prevent the aggravation of wind erosion and desertification, which will increase their awareness of environmental protection and receive punishment for environmental damage. More

A Psammodromus algirus lizard in Spain, where wildfires can confer long-lasting relief from parasites. Credit: Philippe Clement/Nature Picture Library

Ecology
13 July 2021
Destructive fires serve as pest control for lizards

Mediterranean lizards in burnt areas are less likely to be afflicted by mites than their neighbours in unburnt woodlands.

Share on Twitter
Share on Twitter

Share on Facebook
Share on Facebook

Share via E-Mail
Share via E-Mail

Occasional wildfires can help lizards to keep a clean house: the blazes cleanse natural areas of mites that can infest the reptiles’ skin.High-intensity fires in Mediterranean shrublands and woodlands renew vegetation, shoo away seed eaters and keep tree diseases in check. Lola Álvarez-Ruiz at the Desertification Research Centre in Valencia, Spain, and her colleagues were curious to know whether the flames could also be beneficial to animals.Between 2016 and 2018, the researchers sampled Psammodromus algirus, a species of ground-dwelling lizard, in eight burnt and adjacent unburnt areas in Spain. They then counted either how many mites were attached to the creatures’ skin or how many raised scales the lizards had — an indication of previous infection with the parasite.Lizards that lived in unburnt areas were four times more likely to carry mites than were those in recently scorched environments, and were also more likely to have raised scales. The results suggest that there was a lower incidence of parasitism even several years after a fire had occurred.

Proc. R. Soc. B (2021)

Ecology More

Oceanographic context of the sampling locationAll Trichodesmium colonies used in this study were collected from the same phytoplankton net which sampled a surface-ocean Southern Caribbean Sea community (Fig. 1a). At the sampling station the phosphate concentration was low (0.13 μM at 100 m) as is typical in an oligotrophic environment, while the surface dissolved iron concentration was relatively high (2.02 nM at 100 m), consistent with coastal or atmospheric inputs being mobilized in this region (Fig. 1a). By far the most abundant Trichodesmium species at this location was an uncharacterized Trichodesmium thiebautii species, as determined by Trichodesmium-specific metagenome-assembled-genome recruiting (see Table S1).Thirty individual colonies of mixed morphology were separated by hand-picking, immediately examined, and photographed by fluorescent microscopy (385 excitation, >420 nm emission), then frozen individually for particle characterization and metaproteomic analysis (Fig. S1). All colonies used in this study presented as healthy with reddish-orange pigmentation and well-defined shape. When the particles were present they auto-fluoresced in the visual light range, appearing as yellow, red, or blue dots. In general, the particles were concentrated in the center of puff-type colonies, though they were also present in tufts but in smaller numbers. Strikingly, colonies either had many such particles or none at all. Based on prior experimental evidence demonstrating that Trichodesmium colonies can capture mineral particles and access iron from them [14,15,16, 21], we hypothesized that these particles were terrestrially derived minerals (Fig. 1c–h). Therefore, we embarked to understand the morphological heterogeneity by characterizing the particles and the colony’s molecular response to them.Mineralogical characterization of the colony-associated particlesTo find out whether these natural colonies of Trichodesmium had captured iron-rich mineral particles, we performed synchrotron-based micro-X-ray fluorescence (μ-XRF) element mapping of representative colonies with the observed particle associations. Prior evidence of Trichodesmium–particle associations has been based mainly on experimental “feeding” of dust to cultured or captured colonies [15, 16, 20, 22,23,24,25], and it was therefore important to establish these specific Trichodesmium–particle relationships, which developed in nature. We examined one tuft- and two puff-type colonies, all of which had particles associated with them. The element maps were consistent with the hypothesis that there were mineral particles enriched in iron (Fe), copper (Cu), zinc (Zn), titanium/barium (Ti/Ba, which cannot be distinguished by this method), manganese (Mn) and cobalt (Co), though the concentrations approached the limit of detection for the latter two elements (Fig. 2, Figs. S2 and S3). Iron concentrations were particularly high in the particles. Micro-X-ray absorption near-edge structure (μ-XANES) spectra for iron were collected on six particles—three each from the two puffs (Fig. 2 and Fig. S4). The particles contained mineral-bound iron with average oxidation states of 2.6, 2.7, two of oxidation state 2.9, and two of oxidation state 3.0 (Table S2, Fig. S5). While the mineralogy of these particles could not be definitively resolved using μ-XANES, the structure of the absorption edge and post-edge region provided insight into broad mineral groups. Both Fe(III) (oxy/hydro)oxides and mixed-valence iron-bearing minerals consistent with iron silicates were present, suggesting heterogeneous mineral character. While we could not positively identify the silicate mineral phases based on XANES, the spectroscopic similarity of some samples to iron-smectite and the geologic context suggest iron-bearing clays were present (Fig. S5). In this geographic region, iron oxides and clays could be sourced from atmospheric dust deposition, which is common in this region [27, 28] and/or from riverine sources such as the Orinoco and/or Amazon rivers [29, 30].Fig. 2: μ-XRF-based element maps of a Trichodesmium tuft (left) and puff (right) colony (beamsize 3 ×3 μm).White/gray contours, based on the sulfur panel, which is indicative of biomass, have been provided (white = high [S] threshold, gray = lower [S] threshold). The color scale is the same for each image, with the maximum concentration for each element indicated in parentheses; iron is displayed using two scales. Iron oxidation states were determined via μ XANES for three particles in the puff colony, and these are annotated in yellow. The corresponding XANES spectra are shown in Fig. S4 and tabulated data in Table S2.Full size imageThese colony-associated mineral particles likely serve as a simultaneous source of nutritional (Fe, Ni, Co, Mn) and toxic (Cu) metals to the colonies. The elemental composition of the particles is similar to a recent characterization of Trichodesmium-particle associations in the South Atlantic [30]. Release of metals from the particles likely vary over time, with copper, nickel, zinc, and cobalt continually leaching and iron leaching initially, then re-adsorbing back onto particles unless organic chelates assist in solubilization [31].Proteome composition is altered by particle presenceTo understand the impact of the particles on colony diversity and function, we performed comparative metaproteomic analysis of the individual Trichodesmium colonies and their microbiota. Seven puffs without particles, 14 puffs with particles, and 4 tufts with particles were analyzed by a new single-colony metaproteomic method. This approach allowed for the first time the molecular profiles of heterogeneous Trichodesmium colonies to be examined individually. Compared to bulk population-level metaproteomes from this location, which achieved deeper resolution of low-abundance proteins by integrating biomass from 50 to 100 colonies (4478 proteins identified) [32], proteome coverage for the low-biomass single colonies was lower yet sufficient for characterizing colony function (2078 proteins identified, Fig. S6) [32]. In total, 1591 Trichodesmium and 487 epibiont proteins were identified across the 25 single-colony metaproteomes versus 2944 Trichodesmium and 1534 epibiont proteins across triplicate population-level metaproteomes (Tables S3 and S4). Phylogenetic exclusivity was checked such that peptides used to identify epibiont proteins were not present in the Trichodesmium genome (Table S5 and Fig. S7) [33, 34].Trichodesmium’s epibiont community plays crucial roles in colony health and physiology, and together the single-colony proteomes demonstrated a diverse and functionally active microbiome associated with the colonies (Fig. 1b and Fig. S8). The proteomic analysis generally reflected the more abundant, “core” members of the epibiont community as was expected given their low-biomass proportion relative to Trichodesmium cells. Many commonly identified epibiont groups were present including Alphaproteobacteria, Microscilla, and non-Trichodesmium cyanobacteria [12, 35, 36]. In general, epibiont abundance was unaffected by particle presence, with one exception: Firmicute proteins were more abundant in tufts and puffs with particles, suggesting enhanced, possibly anaerobic, metabolism. Greater differences were identified between the puff and tuft morphologies, independent of particle presence and consistent with prior characterizations finding that puffs and tufts harbor distinct epibiont communities [12]. Specifically, eukaryotic proteins were more abundant in puffs compared to tufts. These proteins likely represent copepods due to sequence similarity to the model organism Calanus finmarchicus, and this result is consistent with observed associations between copepods and puffs at this location (Fig. S8B). Notably, proteins from the PVC superphylum, particularly an uncharacterized eukaryote pathogen species related to Chlamydia, were also more abundant in puffs. Eukaryotes are often observed in association with Trichodesmium colonies, but are not always identified due to differences in sampling protocols that could wash them away [12], as well as due to biases in analytical methods, for instance in studies with a focus on bacterial 16S or metagenomic analyses. Overall, the differences in the epibiont community were small, suggesting that these do not explain the observed morphological heterogeneity. We therefore turn our attention to describing how the particles impacted the proteome of Trichodesmium specifically.Mineral presence was associated with significant differences in the Trichodesmium proteome. In total, 131 proteins were differently abundant in puffs with particles versus without particles (p  More

1.Valles, S. M. et al. A picorna-like virus from the red imported fire ant, Solenopsis invicta: Initial discovery, genome sequence, and characterization. Virology 328(1), 151–157. https://doi.org/10.1016/j.virol.2004.07.016 (2004).CAS 
Article 
PubMed 

Google Scholar 
2.Russell, J. A. The ants (Hymenoptera: Formicidae) are unique and enigmatic hosts of prevalent Wolbachia (Alphaproteobacteria) symbionts. Myrmecol. News 16, 7–23 (2012).
Google Scholar 
3.Mongkolsamrit, S. et al. Life cycle, host range and temporal variation of Ophiocordyceps unilateralis/Hirsutella formicarum on Formicine ants. J. Invertebr. Pathol. 111(3), 217–224. https://doi.org/10.1016/j.jip.2012.08.007 (2012).CAS 
Article 
PubMed 

Google Scholar 
4.Espadaler, X. & Santamaria, S. Ecto-and endoparasitic fungi on ants from the Holarctic region. Psyche https://doi.org/10.1155/2012/168478 (2012).Article 

Google Scholar 
5.Boomsma, J. J., Schmid-Hempel, P. & Hughes, W. O. H. Life histories and parasite pressure across the major groups of social insects. in Insect Evolutionary Ecology: Proceedings of the Royal Entomological Society’s 22nd Symposium, Reading, UK, 2003. (CABI Publishing, 2005).6.Lachaud, J. P., Lenoir, A. & Hughes, D. P. Ants and their parasites. Psyche https://doi.org/10.1155/2013/264279 (2013).Article 

Google Scholar 
7.de Bekker, C., Will, I., Das, B. & Adams, R. M. The ants (Hymenoptera: Formicidae) and their parasites: Effects of parasitic manipulations and host responses on ant behavioral ecology. Myrmecol. News 28, 1–24. https://doi.org/10.25849/myrmecol.news_028:001 (2018).Article 

Google Scholar 
8.Witek, M., Barbero, F. & Markó, B. Myrmica ants host highly diverse parasitic communities: From social parasites to microbes. Insect. Soc. 61(4), 307–323. https://doi.org/10.1007/s00040-014-0362-6 (2014).Article 

Google Scholar 
9.Tartally, A., Szücs, B. & Ebsen, J. R. The first records of Rickia wasmannii Cavara, 1899, a myrmecophilous fungus, and its Myrmica Latreille, 1804 host ants in Hungary and Romania (Ascomycetes: Laboulbeniales, Hymenoptera: Formicidae). Myrmecol. News 10, 123 (2007).
Google Scholar 
10.García, F., Espadaler, X., Echave, P. & Vila, R. Hormigas (Hymenoptera, Formicidae) de los acantilados de l’Avenc de Tavertet (Barcelona, Península Ibérica). Boletín de la Sociedad entomológica Aragonesa 47, 363–367 (2010).
Google Scholar 
11.Bezdĕčková, K. & Bezdĕčka, P. First records of the myrmecophilous fungus Rickia wasmannii (Ascomycetes: Laboulbeniales) in the Czech Republic. Acta Musei Moraviae Scientiae Biologicae 96(1), 193–197 (2011).
Google Scholar 
12.Cavara, F. Di una nuova Laboulbeniacea: Rickia wasmannii, nov. gen. e nov.spec. Malpighia 13(1-2), 173–188. (1899).13.Hulden, L. Floristic notes on palaearctic Laboulbeniales (Ascomycetes). Karstenia 25, 1–16 (1985).Article 

Google Scholar 
14.Csata, E., Erős, K. & Markó, B. Effects of the ectoparasitic fungus Rickia wasmannii on its ant host Myrmica scabrinodis: Changes in host mortality and behavior. Insect. Soc. 61(3), 247–252. https://doi.org/10.1007/s00040-014-0349-3 (2014).Article 

Google Scholar 
15.Báthori, F., Csata, E. & Tartally, A. Rickia wasmannii increases the need for water in Myrmica scabrinodis (Ascomycota: Laboulbeniales; Hymenoptera: Formicidae). J. Invertebr. Pathol. 126, 78–82. https://doi.org/10.1016/j.jip.2015.01.005 (2015).Article 
PubMed 

Google Scholar 
16.Báthori, F., Rádai, Z. & Tartally, A. The effect of Rickia wasmannii (Ascomycota, Laboulbeniales) on the aggression and boldness of Myrmica scabrinodis (Hymenoptera, Formicidae). J. Hymenopt. Res. 58, 41. https://doi.org/10.3897/jhr.58.13253 (2017).Article 

Google Scholar 
17.Csata, E. et al. Lock-picks: Fungal infection facilitates the intrusion of strangers into ant colonies. Sci. Rep. UK 7(1), 1–14. https://doi.org/10.1038/srep46323 (2017).CAS 
Article 

Google Scholar 
18.Csata, E., Billen, J., Bernadou, A., Heinze, J. & Markó, B. Infection-related variation in cuticle thickness in the ant Myrmica scabrinodis (Hymenoptera: Formicidae). Insect. Soc. 65(3), 503–506. https://doi.org/10.1007/s00040-018-0628-5 (2018).Article 

Google Scholar 
19.Tartally, A. et al. Patterns of host use by brood parasitic Maculinea butterflies across Europe. Philos. T. Roy. Soc. B. 374(1769), 20180202. https://doi.org/10.1098/rstb.2018.0202 (2019).Article 

Google Scholar 
20.Fauser-Misslin, A., Sadd, B. M., Neumann, P. & Sandrock, C. Influence of combined pesticide and parasite exposure on bumblebee colony traits in the laboratory. J. Appl. Ecol. 51(2), 450–459. https://doi.org/10.1111/1365-2664.12188 (2014).Article 

Google Scholar 
21.Müller, C. B. & Schmid-Hempel, P. Correlates of reproductive success among field colonies of Bombus lucorum: The importance of growth and parasites. Ecol. Entomol. 17, 343–353. https://doi.org/10.1111/j.1365-2311.1992.tb01068.x (1992).Article 

Google Scholar 
22.Porter, S. D. Impact of temperature on colony growth and developmental rates of the ant, Solenopsis invicta. J. Insect Physiol. 34(12), 1127–1133. https://doi.org/10.1016/0022-1910(88)90215-6 (1988).Article 

Google Scholar 
23.Nooten, S. S. & Rehan, S. M. Historical changes in bumble bee body size and range shift of declining species. Biodivers. Conserv. 29, 451–467. https://doi.org/10.1007/s10531-019-01893-7 (2020).Article 

Google Scholar 
24.Schmid-Hempel, P. On the evolutionary ecology of host–parasite interactions: Addressing the question with regard to bumblebees and their parasites. Naturwissenschaften 88, 147–158. https://doi.org/10.1007/s001140100222 (2001).ADS 
CAS 
Article 
PubMed 

Google Scholar 
25.Molet, M., Péronnet, R., Couette, S., Canovas, C. & Doums, C. Effect of temperature and social environment on worker size in the ant Temnothorax nylanderi. J. Therm. Biol. 67, 22–29. https://doi.org/10.1016/j.jtherbio.2017.04.013 (2017).Article 
PubMed 

Google Scholar 
26.Haelewaters, D., Boer, P., Gort, G. & Noordijk, J. Studies of Laboulbeniales (Fungi, Ascomycota) on Myrmica ants (II): Variation of infection by Rickia wasmannii over habitats and time. Anim. Biol. 65(3–4), 219–231. https://doi.org/10.1163/15707563-00002472 (2015).Article 

Google Scholar 
27.Báthori, F., Pfliegler, W. P., Rádai, Z. & Tartally, A. Host age determines parasite load of Laboulbeniales fungi infecting ants: Implications for host-parasite relationship and fungal life history. Mycoscience 30, 1–6. https://doi.org/10.1016/j.myc.2017.09.004 (2017).Article 

Google Scholar 
28.Bezdĕčka, P. & Bezdĕčková, K. First record of the myrmecophilous fungus Rickia wasmannii (Ascomycetes: Laboulbeniales) in Slovakia. Folia Faunistica Slovaca 16(2), 71–72 (2011).
Google Scholar 
29.Csősz, S. & Majoros, G. Ontogenetic origin of mermithogenic Myrmica phenotypes (Hymenoptera, Formicidae). Insect. Soc. 56, 70–76. https://doi.org/10.1007/s00040-008-1040-3 (2009).Article 

Google Scholar 
30.Di Salvo, M. et al. The microbiome of the Maculinea-Myrmica host-parasite interaction. Sci. Rep. UK 9(1), 1–10. https://doi.org/10.1038/s41598-019-44514-7 (2019).ADS 
CAS 
Article 

Google Scholar 
31.Tragust, S., Tartally, A., Espadaler, X. & Billen, J. Histopathology of Laboulbeniales (Ascomycota: Laboulbeniales): Ectoparasitic fungi on ants (Hymenoptera: Formicidae). Myrmecol. News 23, 81–89. https://doi.org/10.25849/myrmecol.news_023:081 (2016).Article 

Google Scholar 
32.Konrad, M., Grasse, A. V., Tragust, S. & Cremer, S. Anti-pathogen protection versus survival costs mediated by an ectosymbiont in an ant host. P. Roy. Soc. B-Biol. Sci. 282(1799), 20141976. https://doi.org/10.1098/rspb.2014.1976 (2015).CAS 
Article 

Google Scholar 
33.Peeters, C., Molet, M., Lin, C. C. & Billen, J. Evolution of cheaper workers in ants: A comparative study of exoskeleton thickness. Biol. J. Linn. Soc. 121, 556–563. https://doi.org/10.1093/biolinnean/blx011 (2017).Article 

Google Scholar 
34.Cammaerts-Tricot, M. C. Production and perception of attractive pheromones by differently aged workers of Myrmica rubra (Hymenoptera Formicidae). Insect. Soc. 21(3), 235–247. https://doi.org/10.1007/BF02226916 (1974).Article 

Google Scholar 
35.Csősz, S. et al. Insect morphometry is reproducible under average investigation standards. Ecol. Evol. 11(1), 547–559. https://doi.org/10.1002/ece3.7075 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
36.Wolak, M. E., Fairbairn, D. J. & Paulsen, Y. R. Guidelines for estimating repeatability. Methods Ecol. Evol. 3(1), 129–137. https://doi.org/10.1111/j.2041-210X.2011.00125.x (2012).Article 

Google Scholar 
37.R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. https://www.R-project.org/. R version 4.0.2 (2020-06-22) (Vienna, Austria, 2020).38.Wright, K. nipals: Principal Components Analysis using NIPALS or Weighted EMPCA, with Gram-Schmidt Orthogonalization. R package version 0.7., published 2020-01-24. https://CRAN.R-project.org/package=nipals/index.html (2020).39.Bates, D., Machler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67(1), 1–48. https://doi.org/10.18637/jss.v067.i01 (2015).40.Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. lmerTest package: Tests in linear mixed effects models. J. Stat. Soft. 82(13), 1–26. https://doi.org/10.18637/JSS.V082.I13 (2017).Article 

Google Scholar 
41.Fox, J. & Weisberg, S. An R Companion to Applied Regression 3rd edn. (Sage, 2019).
Google Scholar 
42.Jennrich, R. I. An asymptotic χ2 test for the equality of two correlation matrices. J. Am. Stat. Assoc. 65(330), 904–912. https://doi.org/10.1080/01621459.1970.10481133 (1970).MathSciNet 
Article 
MATH 

Google Scholar  More

The authors acknowledge valuable discussion with other members of the International Synthesis Consortium, including M. de Araujo Mamede, M. Palmer, J. Kramer, R. Beilinson, J. Arnott and P. Kille. J.B. thanks USGS internal reviewers K. Bagstad and W. Sanford. We also acknowledge all synthesis group members who have been incredibly creative in continuing their research dynamic despite going 100% virtual. M.W. acknowledges funding by the DFG (via iDiv: FZT 118, 02548816). Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government. More