Philippa Kaur

Data collectionIn phase 1 (2017), we first circulated a comprehensive multi-language questionnaire and associated guidelines on protected area personnel numbers to major national protected area agencies, focusing on the 50 countries listed in the WDPA as having the most protected areas. The questionnaire requested information on personnel numbers, type of employers and management levels (from executive to skilled practical workers). Protected area personnel were defined as those spending at least 50% of their work time on protected area-related tasks. The questionnaire also requested information about job titles used for personnel equivalent to rangers. This phase produced usable data for 28 countries/territories.In phase 2 (2018 onwards), we conducted online searches for published data on protected area personnel numbers in the countries/territories not included in the questionnaire survey or where questionnaire responses were incomplete or unclear. The resulting information came from official organizational reports (10 countries/territories), published external studies, project documents and journal papers (35 countries/territories) and websites of protected area organizations or individual sites (9 countries/territories).In phase 3 (2018–2021), we directly requested personal contacts to locate or supply information from official sources both for the remaining countries/territories and to improve or verify data from phases 1 and 2. The minimum data requested were the overall number of protected area personnel, the number of those personnel that could be categorized as rangers, the terrestrial area of protected areas managed by the listed personnel and the source of the information. This phase contributed usable data for 68 countries and territories. Data for a further 17 countries/territories were assembled from multiple sources.The final dataset covered 176 countries/territories: 167 surveyed countries/territories and a further 9 countries/territories that have no WDPA-listed protected areas (Supplementary Table 1), with contributions from more than 150 individuals.Initial data processingTo assess and, where necessary, improve the reliability of data obtained in a wide range of formats and levels of detail and from multiple sources, we scored the data for each country/territory from 0 to 5 for each of four criteria—detail, accuracy, source and age of the data—with a maximum score of 20 (Supplementary Table 1 and Supplementary Fig. 1). For all low-scoring records (a score of less than 15), we sought more-reliable sources in later phases of the study, rejecting any final scores of less than 10.On reviewing the data, we excluded from the analysis protected areas identified in the WDPA as predominantly or entirely marine, Antarctica and countries/territories categorized in the WDPA as polar (Greenland, French Southern Territories, Bouvet Island, Heard Island and McDonald Islands, South Georgia and the South Sandwich Islands). These large, remote and/or largely uninhabited areas are likely to have quite different management models and scales of staffing from terrestrial protected areas (although marine protected areas are also widely understaffed11). For example, in 2012 the 972,000 km2 of Northeast Greenland Protected Area (categorized by the WDPA as polar) was only periodically visited by six two-person teams of naval personnel47, and the 2008 management plan of the 1.51 million km2 Papahānaumokuākea Marine National Monument (Hawai’i, USA) specifies just nine personnel, working in conjunction with several other agencies48. Data for one country were supplied by officials on the agreement that the country was not specifically identified in publications (the country is given the three-letter code ZZZ in relevant tables and figures).Because the format, completeness and level of detail of the data varied widely, from comprehensive personnel lists to single figures, we restricted our raw dataset to six variables that could be consistently extracted from data obtained for each country/territory:

1.

Total number of non-ranger personnel (if known)

2.

Total number of rangers (if known)

3.

Total number of protected area personnel (either the sum of 1 and 2 or provided as an undifferentiated total)

4.

Terrestrial area of protected areas covered by surveyed personnel (km2)

5.

Total terrestrial area of protected areas of the country/territory (km2)

6.

Year of the data

We used the WDPA, official publications and websites to determine (or verify) the area of terrestrial protected areas covered by the personnel listed for each country/territory, using WDPA data if there were discrepancies. Total national terrestrial protected area coverage was taken from the WDPA, with the exception of Turkey, where the area officially reported to the WDPA is significantly less than the nationally published area.The raw data from the survey are shown in Supplementary Table 1.Candidate predictorsTo predict the number of rangers and non-rangers in countries and territories for which we had no data (Statistical analysis), we collected information on the following set of variables, hereafter referred to as candidate predictors:Location dataThe WGS84 latitude and longitude of the centroid of the largest land mass associated with each country/ territory (to obtain the polygons defining the land masses, we used the R package rnaturalearth version 0.1.0; https://github.com/ropensci/rnaturalearth)2020 data from the World Bank (https://data.worldbank.org/indicator)

Area of the country/territory

Population density: the mid-year population divided by land area

Gross domestic product (GDP) in US dollars

GDP per capita in US dollars (GDP divided by mid-year population)

Growth rate of GDP

The proportion of rural inhabitants

The proportion of unemployed inhabitants

The forested proportion of the country/territory

2020 data for each country/territory from the WDPA (https://www.protectedplanet.net/)

The total terrestrial area of WDPA-listed protected areas

The proportion of the terrestrial area of all IUCN-categorized protected areas (Categories I–VI) that falls within protected areas in Category I or II

The proportion of the terrestrial area of all IUCN-categorized protected areas (Categories I–VI) that falls within protected areas in Categories I–IV

2020 data from the Yale Center for Environmental Law and Policy Environmental Performance Index (https://epi.yale.edu/)

Environmental Performance Index (EPI): a composite index using 32 performance indicators across 11 categories

Ecosystem Vitality Index (EVI): an indicator of how well countries preserve, protect and enhance ecosystems and the services they provide

Species Protection Index (SPI): an indicator of the species-level ecological representativeness of each country’s/territory’s protected area network

Not all this information was available for all countries/territories. Most of the missing data were for small territories that account for only a very small proportion of the total area of protected areas worldwide (Supplementary Table 2c).Statistical analysisOur primary objective was to estimate the total number of all personnel engaged in managing all the world’s WDPA-listed terrestrial protected areas and the number categorized as rangers. Our raw data collection yielded full, partial or no information on total personnel and ranger numbers for each country/territory (Supplementary Table 1 shows the completeness of all the data collected). Our first task, therefore, was (1) to impute the information for unsurveyed protected areas on the basis of information from surveyed protected areas within the same countries/territories and (2) to predict those numbers for countries/territories where no information was available on overall personnel numbers and/or ranger numbers on the basis of relationships we could establish between available information and candidate predictors in other countries/territories (Supplementary Table 7). A brief description of these two approaches follows, and full details on the analysis are provided in Supplementary Information.Data imputationFor countries/territories where we had obtained information about numbers of personnel and/or rangers for only some protected areas, our strategy was to populate the unsurveyed protected areas in proportion to the densities of personnel or rangers from the surveyed protected areas of the same countries/territories. For example, for Spain we obtained evidence that there are 619 rangers responsible for protected areas covering 44,328 km2, out of a national total protected area system covering 142,573 km2. To impute the number of rangers for the remaining 98,245 km2, we used the density of rangers in the surveyed area (one ranger per 44,328/619 = 71.6 km2) and applied that to the unsurveyed area, giving a total of 1,991 rangers (619 + (98,245/71.6)). This imputation assumes that unsurveyed areas are staffed at the same density as surveyed areas, whereas in reality the relative densities are likely to vary in unknown ways within different countries/territories. To study the sensitivity of our results to the assumed proportion, we repeated our analysis using the following proportions of the observed densities: 0, 0.25, 0.50, 0.75 and 1.00. This provided a range of personnel numbers from a minimum (based on a proportion of 0) to a presumed maximum (based on a proportion of 1.00). From the data obtained, it was not possible to calculate the actual proportions, but based on the experience of the practitioners in the author team, the unsurveyed areas are highly unlikely to be staffed at higher densities than surveyed areas and, on average, are very likely to be staffed at lower densities. After all, most survey respondents were national or subnational agencies responsible for protected areas subject to stronger formal requirements for protection and management and therefore likely to have larger workforces. Unsurveyed protected areas are more likely to be managed by local entities, with fewer resources, less-stringent management obligations and therefore fewer personnel. The range of proportions we considered to populate unsurveyed areas should therefore yield predictions encompassing the actual (unknown) numbers of rangers and non-rangers with a conservative margin of error. In the main text, we have reported the results of imputation assuming a proportion of 1, which is probably the most optimistic assessment of the current workforce in protected areas within the proportions of the observed densities considered. Results using lower proportions are shown in Extended Data Fig. 2 and Supplementary Tables 4 and 5.Data predictionOur imputation approach was not possible for countries/territories where (1) zero ranger or personnel data had been obtained and (2) specific data had not been obtained that allowed imputation either for rangers or for total personnel (where only total personnel numbers or only ranger numbers had been obtained). To predict the missing information, we used two different statistical approaches: linear mixed models (LMMs)49 and a general implementation of random forests, which we term RF/ETs because it encompasses both random forests sensu stricto (RFs)50 and a variant called extremely randomized trees (ETs)51. LMMs and RFs have been extensively discussed and reviewed in the literature49,52,53. We adopted these approaches because both have proved successful in producing accurate predictions for a wide range of applications and because both are well suited to our data since they both produce predictions from a set of predictors and allow for the consideration of spatial effects54,55. Furthermore, comparing predictions generated through very different methods informs us about the robustness of our results with respect to key statistical assumptions. LMMs come from the ‘data modelling culture’56 and belong to parametric statistics; RF/ETs come from the ‘algorithmic modelling culture’ and belong to non-parametric statistics.We followed the same workflow for both statistical approaches, comprising eight steps: (1) general data preparation; (2) preparation of initial training datasets; (3) selection of predictor variables and of the method used for handling spatial autocorrelation; (4) preparation of final training datasets; (5) fine tuning; (6) final training; (7) preparation of datasets for predictions and simulations; and (8) predictions and simulations (see Supplementary Information for details).Both approaches yielded very similar results with our data. We chose to present the LMM results in the main text, but we provide and compare the results obtained by both approaches in Supplementary Information.SoftwareWe performed all the data analyses using the free open-source statistical software R version 4.157. We used the R package spaMM version 3.9.13 to implement LMMs58 and the R package ranger version 0.13.1 to implement RF/ETs59. To reformat and plot the data, we used the Tidyverse suite of packages60. Details are provided in an R package we specifically developed so that findings presented in this paper can readily be reproduced (see Code availability). Using a workstation with an AMD Ryzen Threadripper 3990 × 64-core processor and 256 GB of RAM, our complete workflow ran in ~3,000 CPU hours.Estimation of required numbers and densities of personnelTo estimate the numbers of personnel and rangers required for effective management of existing protected areas, we referred to ref. 25. This estimates that the minimum budget needed to adequately manage the existing protected area system is US$67.6 billion per year and that current annual expenditure is US$24.3 billion. From these figures, we can calculate that resources invested in the current global system of protected areas are approximately 36% of what is required. We consulted data from https://ourworldindata.org to determine that the proportion of global public expenditure on employee compensation has remained between 21.01% and 23.33% in the years from 2006 to 2019. We obtained these figures from the ‘Government Spending’ section of the site, consulting the chart ‘Share of employee compensation in public spending, 2002 to 2019’ and selecting data for ‘World’. On the basis of this broadly constant proportion and the assumption that total employee compensation is an indicator of total employee numbers, we inferred that current numbers of protected area employees are also around 36% of what is required. We therefore multiplied our estimations of personnel and ranger numbers by 1/0.36 and recalculated the densities on this basis (current requirement = 1/0.36 × current estimate).To estimate staffing requirements for 30% global coverage of protected areas—the global target intended to be reached by 2030—we used the mean personnel and ranger densities calculated as being required at present to ‘populate’ a global area of terrestrial protected areas if increased from the percentage at the time of our study (15.7%) to 30% (current requirement × (0.300/0.157)).Economic calculationsWe based our calculations on published data from 202025, which estimate that expanding the protected areas to 30% would generate higher overall output (revenues) than non-expansion (an extra US$64–454 billion per year by 2050). This figure is only an indicative, partial estimate, generated for the purposes of comparison and to illustrate the substantial return on investment that protected area staff investments imply. Using these figures and our estimates of personnel requirements to ensure effective management of 30% coverage, we calculated the range of sums that each additional protected area staff member has the potential to generate (Supplementary Table 8). For clarity, we rounded these figures to the nearest hundred US dollars in the main text.Our estimates of the gross value added per worker in forestry and agriculture (sectors responsible for similar proportions of the world as protected areas) are included to provide a point of comparison for the figures showing the economic benefit generated per protected area personnel member (see the preceding). The data for the gross annual value of world agricultural production (US$3,550,231,736,000) and the number of workers employed in agriculture (343,527,711) come from the Food and Agriculture Organization of the United Nations30, providing an average gross value of annual agricultural production per worker of US$10,335. We adjusted these 2018 data to 2020 price levels using a deflator based on the US consumer price index (CPI) from the World Economic Outlook database61 (Supplementary Table 9). This ensures that all the economic value data we present are directly comparable for protected area, agricultural and forestry workers. We calculated the gross value of forest production per worker on the basis of direct contribution of forestry of more than US$539 billion to world GDP in 201162 and total forest-sector employment of 11.881 million full-time-equivalent jobs in 201032. These were the most up-to-date global estimates we could locate from credible sources that presented comparable estimates of forest-sector employment and contribution to GDP. This gives an average gross value of forest production per worker of US$45,367 per year. We used the same method as for agriculture to bring these figures to 2020 price levels (Supplementary Table 9). These figures are rounded to the nearest hundred US dollars in the main text. More

We hypothesized that functionally distinct forest types can be mapped at moderate spatial resolutions, using a combination of canopy foliar traits and canopy structure information. Our analysis of LiDAR and imaging spectroscopy data at spatial resolutions ranging from 4 to 200 m (16 m2–40,000 m2), with an emphasis on the 30 m (900 m2) spaceborne hyperspectral spatial resolution, reveals that few remotely sensed canopy properties are needed to successfully identify ecologically distinct forest types at two diverse tropical forest sites in Malaysian Borneo. In testing our second hypothesis that mapped forest types exhibit distinct ecosystem function, we found that forest types identified using remotely sensed leaf P, LMA, Max H, and canopy cover at 20 m height (Cover20) closely align with forest types defined from field-based floristic surveys29,30,31,32,33 and inventory plot-based measurements of growth and mortality rates (Fig. 4b). Our approach, however, enables mapping of their entire spatial extent (Fig. 1) and reveals important structural and functional variation within areas characterized as a single forest type in previous studies (Fig. 3). Current and forthcoming satellite hyperspectral platforms, including PRISMA (30 m), CHIME (20–30 m), and SBG (30 m), have or will have comparable spectral resolution, higher temporal revisits, and much greater geographic coverage. The ability to conduct this type of analysis using remote sensing measurements at 30 m resolution suggests that our method can be applied to these emerging spaceborne imaging spectroscopy data to reveal important differences in structure and function across the world’s tropical forests.Nested functional forest types revealedTo test our first hypothesis, rather than making an a priori decision about the number of k-means clusters (k), we explored the capacity of remotely sensed data to reveal ecologically relevant variation in forest types. Baldeck and Asner took a similar unsupervised approach to estimating beta diversity in South Africa34. Because the choice of k directly influences analysis outcomes, careful selection of k is required. Different approaches for identifying the number of clusters, using the Gapk and Wk elbow metrics35, yielded varying optimal numbers of clusters for the Sepilok and Danum landscapes (Fig. 1, Supplementary Figs. 4 and 5). However, at both sites, a comparison of results based on different values of k revealed ecologically meaningful structural and functional differences and graduated transitions between forest types (Fig. 2, Supplementary Figs. 7 and 8), indicating that the exploration of traits that aggregate or separate forest types as k changes is a valuable exercise. Overlap between the remotely sensed forest type boundaries and inventory plots within distinct forest types indicate that the series of clustered forests align closely with forest types defined based on in situ data on species composition and ecosystem structure. In part, this type of analysis requires careful selection of the number of clusters. Additionally, however, we gained valuable insights via the exploration of varying numbers of clusters as it relates to biologically meaningful categorization of forest types. Extending this method to other parts of the tropics will require similar decision-making, which will either require user input, or the development of robust automated algorithms for selecting k.Forest types capture differences in ecosystem dynamicsWe further evaluated the canopy traits and structural attributes that were most critical for mapping distinct forest types, hypothesizing that mapped forest types exhibit distinct ecosystem function. Forest types revealed by the cluster analyses were distributed along the leaf economic spectrum, where the leaf economic spectrum characterizes a tradeoff in plant growth strategies36. LMA, which can covary strongly with leaf N and P, is a key indicator of plant growth strategies along the spectrum37. At the slow-return end of the leaf economics spectrum, plants in nutrient-poor conditions with low leaf nutrient concentrations invest in leaf structure and defense, expressed as high LMA, strategizing longer-lived, tougher leaves with slower decomposition rates. This strategy comes at the cost of slower growth. At the quick-return end of the spectrum, plants in nutrient-rich environments with higher leaf nutrient concentrations invest less in structure and defense, enabling faster growth and more rapid leaf turnover, i.e., shorter leaf lifespans. This quick-return growth strategy supports higher photosynthetic rates and more rapid carbon gain36.In this study, the principal components and clustering results yielded forest types that are indicative of community level differences associated with leaf economic spectrum differences. The nutrient rich sites (Danum1 and Danum2, Supplementary Fig. 8) show high canopy N and P and low LMA compared to the nutrient poor and acidic sites (Sandstone and Kerangas), which contributes to lower leaf photosynthetic capacity (Vcmax) and growth (Fig. 4b). Foliar N:P also increased with site fertility, confirming that tropical forests are primarily limited by phosphorus, and not nitrogen38,39, with large implications for carbon sequestration in these forests. Orthogonal differences in canopy structure and architecture between Danum forest types and Sepilok Sandstone and Alluvial forests could be indicative of ecosystem scale differences in the sensitivity of these forests to endogenous disturbance processes40.The significant differences in aboveground carbon stocks and growth and mortality rates between forest types further suggests strong differences in ecosystem dynamics. In general, growth rates varied inversely to aboveground carbon, and higher aboveground carbon corresponded to lower mortality rates. As an example, the Sepilok sandstone forests, which are largely comprised of slow-growing dipterocarp species29,33, had the highest median aboveground carbon (236 Mg C ha−1), with higher canopy P and N, and lower LMA. The taller canopy and low canopy leaf nutrient concentrations are consistent with the low growth and mortality rates found in the sandstone forest, indicating a slow-growth strategy yielding larger trees and higher aboveground carbon stocks. In contrast, alluvial forests exhibit high turnover with mortality and growth rates higher relative to Sandstone forests corresponding to lower aboveground carbon on average. Kerangas forests exhibited low aboveground carbon despite an intermediate plot-level growth rate, and mortality rates that were significantly lower than the Danum or alluvial forest types. Kerangas forests, which were characterized by the highest LMA, lowest foliar P and N (Fig. 2a), and the lowest plot-level aboveground carbon density (186 Mg C ha−1; Fig. 4a), are known to have higher stem densities, lower canopy heights, and long-lived leaves5,32,41, suggesting well-developed strategies for nutrient retention42. Interestingly, despite significantly different aboveground carbon and demography, the kerangas and sandstone forests did not differ in LAI or canopy architecture (P:H); although maximum height, Cover20, and Hpeak LAI were significantly higher in the sandstone forest, highlighting the need to account for differences beyond LAI when scaling processes from leaves to ecosystems.In addition, when three forest types were distinguished at Sepilok, the alluvial inventory plot had significantly higher aboveground carbon than the remote sensing-derived alluvial forest extent (Fig. 4a, p  More

Sakschewski, B. et al. Leaf and stem economics spectra drive diversity of functional plant traits in a dynamic global vegetation model. Glob. Change Biol. 21, 2711–2725 (2015).Article 

Google Scholar 
Berzaghi, F. et al. Towards a new generation of trait-flexible vegetation models. Trends Ecol. Evol. 35, 191–205 (2020).Article 
PubMed 

Google Scholar 
Bruelheide, H. et al. Global trait–environment relationships of plant communities. Nat. Ecol. Evol. 2, 1906–1917 (2018).Article 
PubMed 

Google Scholar 
Joswig, J. S. et al. Climatic and soil factors explain the two-dimensional spectrum of global plant trait variation. Nat. Ecol. Evol. 6, 36–50 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
van Bodegom, P. M., Douma, J. C. & Verheijen, L. M. A fully traits-based approach to modeling global vegetation distribution. Proc. Natl Acad. Sci. USA 111, 13733–13738 (2014).PubMed Central 

Google Scholar 
Moreno Martínez, A. et al. A methodology to derive global maps of leaf traits using remote sensing and climate data. Remote Sens. Environ. 218, 69–88 (2018).Article 

Google Scholar 
Pérez-Harguindeguy, N. et al. New handbook for standardized measurment of plant functional traits worldwide. Aust. J. Bot. 23, 167–234 (2013).Article 

Google Scholar 
Kattge, J. et al. TRY—a global database of plant traits. Glob. Change Biol. 17, 2905–2935 (2011).Article 

Google Scholar 
Kattge, J. et al. TRY plant trait database-enhanced coverage and open access. Glob. Change Biol. 26, 119–188 (2020).Article 

Google Scholar 
Jetz, W. et al. Monitoring plant functional diversity from space. Nat. Plants 2, 16024 (2016).Article 
PubMed 

Google Scholar 
Butler, E. E. et al. Mapping local and global variability in plant trait distributions. Proc. Natl Acad. Sci. USA 114, E10937–E10946 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Boonman, C. C. et al. Assessing the reliability of predicted plant trait distributions at the global scale. Glob. Ecol. Biogeogr. 29, 1034–1051 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Madani, N. et al. Future global productivity will be affected by plant trait response to climate. Sci. Rep. 8, 2870 (2018).Vallicrosa, H. et al. Global distribution and drivers of forest biome foliar nitrogen to phosphorus ratios (N:P). Glob. Ecol. Biogeogr. 31, 861–871 (2022).Article 

Google Scholar 
Meyer, H. & Pebesma, E. Predicting into unknown space? Estimating the area of applicability of spatial prediction models. Methods Ecol. Evol. 12, 1620–1633 (2021).Article 

Google Scholar 
Schiller, C. et al. Deep learning and citizen science enable automated plant trait predictions from photographs. Sci. Rep. 11, 16395 (2021).Aguirre-Gutiérrez, J. et al. Pantropical modelling of canopy functional traits using sentinel-2 remote sensing data. Remote Sens. Environ. 252, 112–122 (2021).Article 

Google Scholar 
Homolova, L. et al. Review of optical-based remote sensing for plant trait mapping. Ecol. Complex. 15, 1–16 (2013).Article 

Google Scholar 
Van Cleemput, E. et al. The functional characterization of grass-and-shrubland ecosystems using hyperspectral remote sensing: trends, accuracy and moderating variables. Remote Sens. Environ. 209, 747–763 (2018).Article 

Google Scholar 
Kattenborn, T., Fassnacht, F. E. & Schmidtlein, S. Differentiating plant functional types using reflectance: which traits make the difference? Remote Sens. Ecol. Conserv. 5, 5–19 (2019).Article 

Google Scholar 
Hauser, L. T. et al. Explaining discrepancies between spectral and in-situ plant diversity in multispectral satellite earth observation. Remote Sens. Environ. 265, 112684 (2021).Article 

Google Scholar 
Wäldchen, J. & Mäder, P. Plant species identification using computer vision techniques: a systematic literature review. Arch. Comput. Methods Eng. 25, 507–543 (2018).Article 
PubMed 

Google Scholar 
Jones, H. G. What plant is that? Tests of automated image recognition apps for plant identification on plants from the British flora. AoB Plants 12, plaa052 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Hampton, S. E. et al. Big data and the future of ecology. Front. Ecol. Environ. 11, 156–162 (2013).Article 

Google Scholar 
WÜest, R. O. et al. Macroecology in the age of big data—where to go from here? J. Biogeogr. 47, 1–12 (2020).Article 

Google Scholar 
Mäder, P. et al. The Flora Incognita app—interactive plant species identification. Methods Ecol. Evol. 12, 1335–1342 (2021).Article 

Google Scholar 
Di Cecco, G. J. et al. Observing the observers: how participants contribute data to iNaturalist and implications for biodiversity science. BioScience 71, 1179–1188 (2021).Article 

Google Scholar 
Mahecha, M. D. et al. Crowd-sourced plant occurrence data provide a reliable description of macroecological gradients. Ecography 44, 1131–1142 (2021).Article 

Google Scholar 
Botella, C. et al. Jointly estimating spatial sampling effort and habitat suitability for multiple species from opportunistic presence-only data. Methods Ecol. Evol. 12, 933–945 (2021).Article 

Google Scholar 
iNaturalist Research-Grade Observations (GBIF, accessed 5 January 2022); https://www.gbif.org/dataset/50c9509d-22c7-4a22-a47d-8c48425ef4a7Callaghan, C. T. et al. Three frontiers for the future of biodiversity research using citizen science data. BioScience 71, 55–63 (2020).
Google Scholar 
Dickinson, J. L., Zuckerberg, B. & Bonter, D. N. Citizen science as an ecological research tool: challenges and benefits. Ann. Rev. Ecol. Evol. Syst. 41, 149–172 (2010).Article 

Google Scholar 
Kosmala, M. et al. Assessing data quality in citizen science. Front. Ecol. Environ. 14, 551–560 (2016).Article 

Google Scholar 
Boakes, E. H. et al. Patterns of contribution to citizen science biodiversity projects increase understanding of volunteers’ recording behaviour. Sci. Rep. 6, 33051 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bowler, D.E. et al. Temporal trends in the spatial bias of species occurrence records. Ecography 2022, e06219 (2022). https://doi.org/10.1111/ecog.06219GBIF Occurrence Download (GBIF, 4 January 2022); https://doi.org/10.15468/dl.34tjreBruelheide, H. et al. sPlot—a new tool for global vegetation analyses. journal of vegetation science. J. Veg. Sci. 30, 161–186 (2019).Article 

Google Scholar 
Sabatini, F. et al. sPlotOpen—an environmentally balanced, open access, global dataset of vegetation plots. Glob. Ecol. Biogeogr. 30, 1740–1764 (2021).Article 

Google Scholar 
Whittaker, R.H. et al. Communities and Ecosystems (Macmillan/Collier Macmillan, 1970).Olson, D. M. et al. Terrestrial ecoregions of the world: a new map of life on Earth: a new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. BioScience 51, 933–938 (2001).Article 

Google Scholar 
Joswig, J., Wirth, C. & Schuman, M. Climatic and soil factors explain the two-dimensional spectrum of global plant trait variation. Nat. Ecol. Evol. 6, 36–50 (2022).Article 
PubMed 

Google Scholar 
Díaz, S. et al. The global spectrum of plant form and function. Nature 529, 167–171 (2016).Article 
PubMed 

Google Scholar 
Ploton, P. et al. Spatial validation reveals poor predictive performance of large-scale ecological mapping models. Nat. Commun. 11, 4540 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Meyer, H. & Pebesma, E. Spatial validation reveals poor predictive performance of large-scale ecological mapping models. Methods Ecol. Evol. 12, 1620–1633 (2021).Article 

Google Scholar 
Schrodt, F. et al. Bhpmf—a hierarchical Bayesian approach to gap filling and trait prediction for macroecology and functional biogeography. Glob. Ecol. Biogeogr. 24, 1510–1521 (2015).Article 

Google Scholar 
Kuppler, J. et al. Global gradients in intraspecific variation in vegetative and floral traits are partially associated with climate and species richness. Glob. Ecol. Biogeogr. 29, 992–1007 (2020).Article 

Google Scholar 
Scheiter, S., Langan, L. & Higgins, S. I. Next-generation dynamic global vegetation models: learning from community ecology. New Phytol. 198, 957–969 (2013).Article 
PubMed 

Google Scholar 
Taubert, F. et al. Confronting an individual-based simulation model with empirical community patterns of grasslands. PLoS ONE 15, e0236546 (2020).Roger, E. & Klistorner, S. (2016) Bioblitzes help science communicators engage local communities in environmental research. J. Sci. Commun. https://doi.org/10.22323/2.15030206 (2016).Legendre, P. & Legendre, L. Numerical Ecology 3rd edn (Elsevier, 2012).Warton, D. I. et al. Smatr 3—an R package for estimation and inference about allometric lines. Methods Ecol Evol 3, 257–259 (2012).Article 

Google Scholar 
Wolf, S. et al. iNaturalist_traits: iNaturalist trait maps version 1 (January 5, 2022) Zenodo https://doi.org/10.5281/zenodo.6671891 (2022). More

Field samplingLake Masoko fish were chased into fixed gill nets and SCUBA by a team of professional divers at different target depths determined by diver depth gauge (12× male benthic, 12× male littoral). Riverine fish (11× Mbaka River and 1× Itupi river) were collected by local fishermen. On collection, all fish were euthanized using clove oil. Collection of wild fish was done in accordance with local regulations and permits in 2015, 2016, 2018 and 2019. On collection, fish were immediately photographed with color and metric scales, and tissues were dissected and stored in RNAlater (Sigma-Aldrich); some samples were first stored in ethanol. Only male specimens (showing bright nuptial coloration) were used in this study for the practical reason of avoiding any misassignment of individuals to the wrong population (only male individuals show clear differences in phenotypes and could therefore be reliably assigned to a population). Furthermore, we assumed that any epigenetic divergence relevant to speciation should be contributing to between-population differences in traits possessed by both sexes (habitat occupancy, diet). To investigate the role of epigenetics in phenotypic diversification and adaptation to different diets, homogenized liver tissue – a largely homogenous and key organ involved in dietary metabolism, hormone production and hematopoiesis – was used for all RNA-seq and WGBS experiments.Common-garden experimentCommon-garden fish were bred from wild-caught fish specimens, collected and imported at the same time by a team of professional aquarium fish collectors according to approved veterinary regulations of the University of Bangor, UK. Wild-caught fish were acclimatized to laboratory tanks and reared to produce first-generation (G1) common-garden fish, which were reared under the same controlled laboratory conditions in separate tanks (light–dark cycles, diet: algae flakes daily, 2–3 times weekly frozen diet) for approximately 6 months (post hatching). G1 adult males showing bright nuptial colors were culled at the same biological stages (6 months post hatching) using MS222 in accordance with the veterinary regulations of the University of Bangor, UK. Immediately on culling, fish were photographed and tissues collected and snap-frozen in tubes.Stable isotopesTo assess dietary/nutritional profiles in the three ecomorph populations, carbon (δ13C) and nitrogen (δ15N) isotope analysis of muscle samples (for the same individuals as RRBS; 12, 12 and 9 samples for benthic, littoral and riverine populations, respectively) was undertaken by elemental analyzer isotope ratio mass spectrometry by Iso-Analytical Limited. It is important to note that stable isotope analysis does not depend on the use of the same tissue as the ones used for the RRBS/WGBS samples45. Normality tests (Shapiro–Wilk, using the R package rstatix v.0.7.0), robust for small sample sizes, were performed to assess sample deviation from a Gaussian distribution. Levene’s test for homogeneity of variance was then performed (R package carData v.3.0-5) to test for homogeneity of variance across groups. Finally, Welch’s ANOVA was performed followed by Games–Howell all-pairs comparison tests with adjusted P value using Tukey’s method (rstatix v.0.7.0). Mean differences in isotope measurements and 95% CI mean differences were calculated using Dabestr v.0.3.0 with 5,000 bootstrapped resampling.Throughout this manuscript, all box plots are defined as follows: centre line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; points, outliers.RNA-seqNext-generation sequencing library preparationTotal RNA from liver tissues stored in RNAlater was extracted using a phenol/chloroform approach (TRIzol reagent; Sigma-Aldrich). Of note, when tissues for bisulphite sequencing samples were not available, additional wild-caught samples were used (Supplementary Table 3). The quality and quantity of RNA extraction were assessed using TapeStation (Agilent Technologies), Qubit and NanoDrop (Thermo Fisher Scientific). Next-generation sequencing (NGS) libraries were prepared using poly(A) tail-isolated RNA fraction and sequenced on a NovaSeq system (S4; paired-end 100/150 bp; Supplementary Table 3), yielding on average 32.9 ± 3.9 Mio reads.Read alignment and differential gene expression analysisAdaptor sequence in reads, low-quality bases (Phred score  More

Díaz, S. et al. Pervasive human-driven decline of life on Earth points to the need for transformative change. Science 366, eaax3100 (2019).Chase, J. M. et al. Species richness change across spatial scales. Oikos 128, 1079–1091 (2019).Article 

Google Scholar 
Vellend, M. et al. Global meta-analysis reveals no net change in local-scale plant biodiversity over time. Proceedings of the National Academy of Sciences 110, 19456–19459 (2013).Article 
ADS 
CAS 

Google Scholar 
Blowes, S. A. et al. The geography of biodiversity change in marine and terrestrial assemblages. Science 366, 339–345 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Riedel, T., Polley, H. & Klatt, S. Germany. in National Forest Inventories (eds. Vidal, C., Alberdi, I. A., Hernández Mateo, L. & Redmond, J. J.) 405–421, https://doi.org/10.1007/978-3-319-44015-6 (Springer International Publishing, 2016).Braun-Blanquet, J. Pflanzensoziologie. Grundzüge der Vegetationskunde. vol. Seite: (Julius Springer, 1928).Bernhardt-Römermann, M. et al. Drivers of temporal changes in temperate forest plant diversity vary across spatial scales. Glob Change Biol 21, 3726–3737 (2015).Article 
ADS 

Google Scholar 
Ahrns, C. & Hofmann, G. Vegetationsdynamik und Florenwandel im ehemaligen mitteldeutschen Waldschutzgebiet ‘Hainich’ im Intervall 1963–1995. Hercynia N.F. 31, 33–64 (1998).
Google Scholar 
Dittmann, T., Heinken, T. & Schmidt, M. Die Wälder von Magdeburgerforth (Fläming, Sachsen-Anhalt) – eine Wiederholungsuntersuchung nach sechs Jahrzehnten, https://doi.org/10.14471/2018.38.009 (2018).Günther, K., Schmidt, M., Quitt, H. & Heinken, T. Veränderungen der Waldvegetation im Elbe-Havelwinkel von 1960 bis 2015. Tuexenia 41, 53–85 (2021).
Google Scholar 
Janiesch, P. Vegetationsökologische Untersuchungen in einem Erlenbruchwald im nördlichen Münsterland. 25 Jahre im Vergleich. Abhandlungen aus dem Westfälischen Museum für Naturkunde 71–80 (2003).Naaf, T. & Wulf, M. Habitat specialists and generalists drive homogenization and differentiation of temperate forest plant communities at the regional scale. Biological Conservation 143, 848–855 (2010).Article 

Google Scholar 
Reinecke, J., Klemm, G. & Heinken, T. Vegetation change and homogenization of species composition in temperate nutrient deficient Scots pine forests after 45 yr. J Veg Sci 25, 113–121 (2014).Article 

Google Scholar 
Mölder, A., Streit, M. & Schmidt, W. When beech strikes back: How strict nature conservation reduces herb-layer diversity and productivity in Central European deciduous forests. Forest Ecology and Management 319, 51–61 (2014).Article 

Google Scholar 
Fischer, C., Parth, A. & Schmidt, W. Vegetationsdynamik in Buchen-Naturwäldern. Ein Vergleich aus Süd-Niedersachsen. Hercynia N.F. 45–68 (2009).Schmidt, W. Die Naturschutzgebiete Hainholz und Staufenberg am Harzrand – Sukzessionsforschung in Buchenwäldern ohne Bewirtschaftung (Exkursion E). Tuexenia 22, 151–213 (2002).
Google Scholar 
Strubelt, I., Diekmann, M. & Zacharias, D. Changes in species composition and richness in an alluvial hardwood forest over 52 yrs. J Veg Sci 28, 401–412 (2017).Article 

Google Scholar 
Strubelt, I., Diekmann, M., Peppler-Lisbach, C., Gerken, A. & Zacharias, D. Vegetation changes in the Hasbruch forest nature reserve (NW Germany) depend on management and habitat type. Forest Ecology and Management 444, 78–88 (2019).Article 

Google Scholar 
Wilmanns, O. & Bogenrieder, A. Veränderungen der Buchenwälder des Kaiserstuhls im Laufe von vier Jahrzehnten und ihre Interpretation – pflanzensoziologische Tabellen als Dokumente. Abh. Landesmus. Naturk. Münster Westfalen 48, 55–79 (1986).
Google Scholar 
Huwer, A. & Wittig, R. Changes in the species composition of hedgerows in the Westphalian Basin over a thirty-five-year period. Tuexenia 32, 31–53 (2012).
Google Scholar 
Immoor, A., Zacharias, D., Müller, J. & Diekmann, M. A re-visitation study (1948–2015) of wet grassland vegetation in the Stedinger Land near Bremen, North-western Germany, https://doi.org/10.14471/2017.37.013 (2017).Rosenthal, G. Erhaltung und Regeneration von Feuchtwiesen. Vegetationsökologische Untersuchungen auf Dauerflächen. Diss. Bot. 182, 1–283 (1992).
Google Scholar 
Poptcheva, K., Schwartze, P., Vogel, A., Kleinebecker, T. & Hölzel, N. Changes in wet meadow vegetation after 20 years of different management in a field experiment (North-West Germany). Agriculture, Ecosystems & Environment 134, 108–114 (2009).Article 

Google Scholar 
Diekmann, M. et al. Patterns of long‐term vegetation change vary between different types of semi‐natural grasslands in Western and Central Europe. J Veg Sci 30, 187–202 (2019).Article 

Google Scholar 
Hundt, R. Ökologisch‐geobotanische Untersuchungen an den mitteldeutschen Wiesengesellschaften unter besonderer Berücksichtigung ihres Wasserhaushaltes und ihrer Veränderung durch die Intensivbewirtschaftung. (Wehry-Druck OHG, 2001).Kuhn, G., Heinz, S. & Mayer, F. Grünlandmonitoring Bayern. Ersterhebung der Vegetation 2002–2008. Schriftenreihe LfL Bayerische Landesanstalt für Landwirtschaft 3, 1–161 (2011).
Google Scholar 
Raehse, S. Veränderungen der hessischen Grünlandvegetation seit Beginn der 50er Jahre am Beispiel ausgewählter Tal- und Bergregionen Nord- und Mittelhessens. (University Press GmbH, 2001).Scheidel, U. & Bruelheide, H. Versuche zur Beweidung von Bergwiesen im Harz. Hercynia N.F 37, 87–101 (2004).
Google Scholar 
Sommer, S. & Hachmöller, B. Auswertung der Vegetationsaufnahmen von Dauerbeobachtungenflächen auf Bergwiesen im NSG Oelsen bei variierter Mahd im Vergleich zur Brache. Ber. Arbeitsgem. Sächs. Bot. N.F. 18, 99–135 (2001).
Google Scholar 
Wegener, U. Vegetationswandel des Berggrünlands nach Untersuchungen von 1954 bis 2016 – Wege zur Erhaltung der Bergwiesen. Mountain grasslands vegetation change after research from 1954 to 2016 – ways to preserve mountain meadows. Abhandlungen und Berichte aus dem Museum Heineanum 11, 35–101 (2018).
Google Scholar 
Wittig, B., Müller, J. & Mahnke-Ritoff, A. Talauen-Glatthaferwiesen im Verdener Wesertal (Niedersachsen). Tuexenia 39, 249–265 (2019).
Google Scholar 
Heinrich, W., Marstaller, R. & Voigt, W. Eine Langzeitstudie zur Sukzession in Halbtrockenrasen – Strukturwandlungen in einer Dauerbeobachtungsfläche im Naturschutzgebiet “Leutratal und Cospoth” bei Jena (Thüringen). Artenschutzreport Jena 30, 1–80 (2012).
Google Scholar 
Hüllbusch, E., Brand, L. M., Ende, P. & Dengler, J. Little vegetation change during two decades in a dry grassland complex in the Biosphere Reserve Schorfheide-Chorin (NE Germany). Tuexenia 36, 395–412 (2016).
Google Scholar 
Knapp, R. Dauerflächen-Untersuchungen über die Einwirkung von Haustieren und Wild während trockener und feuchter Zeiten in Mesobromion-Halbtrockenrasen in Hessen. Mitt. Florist.-Soziol. Arbeitsgem. N.F. 19/20, 269–274 (1977).
Google Scholar 
Matesanz, S., Brooker, R. W., Valladares, F. & Klotz, S. Temporal dynamics of marginal steppic vegetation over a 26-year period of substantial environmental change: Temporal dynamics of marginal steppic vegetation over a 26-year period. Journal of Vegetation Science 20, 299–310 (2009).Article 

Google Scholar 
Schwabe, A., Zehm, A., Nobis, M., Storm, C. & Süß, K. Auswirkungen von Schaf-Erstbeweidung auf die Vegetation primär basenreicher Sand-Ökosysteme. NNA Berichte 1/2004, 39–54 (2004).
Google Scholar 
Schwabe, A., Süss, K. & Storm, C. What are the long-term effects of livestock grazing in steppic sandy grassland with high conservation value? Results from a 12-year field study. Tuexenia 33, 189–212 (2013).
Google Scholar 
Peppler‐Lisbach, C., Stanik, N., Könitz, N. & Rosenthal, G. Long‐term vegetation changes in Nardus grasslands indicate eutrophication, recovery from acidification, and management change as the main drivers. Appl Veg Sci 23, 508–521 (2020).Article 

Google Scholar 
Peppler-Lisbach, C. & Könitz, N. Vegetationsveränderungen in Borstgrasrasen des Werra-Meißner-Gebietes (Hessen, Niedersachsen) nach 25 Jahren. Tuexenia 37, 201–228 (2017).
Google Scholar 
Wittig, B., Müller, J., Quast, R. & Miehlich, H. Arnica montana in Calluna-Heiden auf dem Schießplatz Unterlüß (Niedersachsen). Tuexenia 40, 131–146 (2020).
Google Scholar 
Rumpf, S. B. et al. Range dynamics of mountain plants decrease with elevation. Proc Natl Acad Sci USA 115, 1848–1853 (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kudernatsch, T. et al. Vegetationsveränderungen alpiner Kalk-Magerrasen im Nationalpark Berchtesgaden während der letzten drei Jahrzehnte. Tuexenia 36, 205–221 (2016).
Google Scholar 
Poschlod, P. et al. Long‐term monitoring in rivers of south Germany since the 1970ies. Macrophytes as indicators for the assessment of water quality. in Long‐term ecological research. Between Theory and Application (eds. Müller, F., Baessler, C., Schubert, H. & Klotz, S.) 189–199 (Springer, 2006).Dierschke, H. Dynamik und Konstanz an naturnahen Flussufern. 27 Jahre Dauerflächenuntersuchungen am Oderufer (Harzvorland). Braunschweiger Geobotanische Arbeiten 9, 119–138 (2008).
Google Scholar 
Kreyling, J. et al. Rewetting does not return drained fen peatlands to their old selves. Nat Commun 12, 5693 (2021).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bohn, U. & Schniotalle, S. Hochmoor-, Grünland- und Waldrenaturierung im Naturschutzgebiet ‘Rotes Moor’/Hohe Rhön 1981–2001. (Landwirtschaftsverlag, 2008).Koch, M. & Jurasinski, G. Four decades of vegetation development in a percolation mire complex following intensive drainage and abandonment. Plant Ecology & Diversity 8, 49–60 (2015).Article 

Google Scholar 
Walther, K. Die Vegetation des Maujahn 1984. Wiederholung der vegetationskundlichen Untersuchung eines wendländischen Moores. Tuexenia 6, 145–193 (1986).
Google Scholar 
Berg, C. & Mahn, E.-G. Anthropogene Vegetationsveränderungen der Straßenrandvegetation in den letzten 30 Jahren – die Glatthaferwiesen des Raumes Halle/Saale. Tuexenia 10, 185–195 (1990).
Google Scholar 
Meyer, S., Wesche, K., Krause, B. & Leuschner, C. Dramatic losses of specialist arable plants in Central Germany since the 1950s/60s – a cross-regional analysis. Diversity Distribution 19, 1175–1187 (2013).Article 

Google Scholar 
Meyer, S., Wesche, K., Krause, B. & Leuschner, C. Veränderungen in der Segetalflora in den letzten Jahrzehnten und mögliche Konsequenzen für Agrarvögel. Julius-Kühn-Archiv 442, 64–78 (2013).
Google Scholar 
Kutzelnigg, H. Veränderungen der Ackerwildkrautflora im Gebiet um Moers/Niederrhein seit 1950 und ihre Ursachen. Tuexenia 4, 81–102 (1984).
Google Scholar 
Milligan, G., Rose, R. J. & Marrs, R. H. Winners and losers in a long-term study of vegetation change at Moor House NNR: Effects of sheep-grazing and its removal on British upland vegetation. Ecological Indicators 68, 89–101 (2016).Article 

Google Scholar 
Wittig, B., Waldman, T. & Diekmann, M. Veränderungen der Grünlandvegetation im Holtumer Moor über vier Jahrzehnte. Hercynia N.F 40, 285–300 (2007).
Google Scholar 
Henning, K., Lorenz, A., von Oheimb, G., Härdtle, W. & Tischew, S. Year-round cattle and horse grazing supports the restoration of abandoned, dry sandy grassland and heathland communities by supressing Calamagrostis epigejos and enhancing species richness. Journal for Nature Conservation 40, 120–130 (2017).Article 

Google Scholar 
Blüml, V. Langfristige Veränderungen von Flora und Vegetation des Grünlandes in der Dümmerniederung (Niedersachsen) unter dem Einfluss von Naturschutzmaßnahmen. (Bremen, 2011).Von Oheimb, G. et al. Halboffene Weidelandschaft Höltigbaum. Perspektiven für den Erhalt und die naturverträgliche Nutzung von Offenlandlebensräumen. (Landwirschaftsverlag, 2006).Dornelas, M. et al. BioTIME: A database of biodiversity time series for the Anthropocene. Global Ecol Biogeogr 27, 760–786 (2018).Article 

Google Scholar 
Pimm, S. L. et al. The biodiversity of species and their rates of extinction, distribution, and protection. Science 344, 1246752–1246752 (2014).Article 
CAS 
PubMed 

Google Scholar 
Barnosky, A. D. et al. Has the Earth’s sixth mass extinction already arrived? Nature 471, 51–57 (2011).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Vellend, M. The Biodiversity Conservation Paradox. Am. Sci. 105, 94 (2017).Article 

Google Scholar 
Cardinale, B. J., Gonzalez, A., Allington, G. R. H. & Loreau, M. Is local biodiversity declining or not? A summary of the debate over analysis of species richness time trends. Biological Conservation 219, 175–183 (2018).Article 

Google Scholar 
Perring, M. P. et al. Understanding context dependency in the response of forest understorey plant communities to nitrogen deposition. Environmental Pollution 242, 1787–1799 (2018).Article 
CAS 
PubMed 

Google Scholar 
Steinbauer, M. J. et al. Accelerated increase in plant species richness on mountain summits is linked to warming. Nature 556, 231–234 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Braun-Blanquet, J. Prinzipien einer Systematik der Pflanzengesellschaften auf floristischer Grundlage. Jahrb. St. Gallischen Naturwiss. Ges. 57, 305–351 (1921).
Google Scholar 
Becking, R. W. The Zürich-Montpellier school of phytosociology. Bot. Rev. 23, 411–488 (1957).Article 

Google Scholar 
Bruelheide, H. et al. sPlot – A new tool for global vegetation analyses. J Veg Sci 30, 161–186 (2019).Article 

Google Scholar 
O L Pescott, T A Humphrey & K J Walker. A short guide to using British and Irish plant occurrence data for research, https://doi.org/10.13140/RG.2.2.33746.86720 (2018).Eichenberg, D. et al. Widespread decline in Central European plant diversity across six decades. Global Change Biology 27, 1097–1110 (2021).Article 
ADS 
CAS 

Google Scholar 
Chytrý, M. et al. European Vegetation Archive (EVA): an integrated database of European vegetation plots. Appl Veg Sci 19, 173–180 (2016).Article 

Google Scholar 
Van der Maarel, E. Transformation of cover-abundance values in phytosociology and its effects on community similarity. Vegetatio 39, 97–114 (1979).Article 

Google Scholar 
Tichý, L. et al. Optimal transformation of species cover for vegetation classification. Appl Veg Sci 23, 710–717 (2020).Article 

Google Scholar 
Podani, J. Braun-Blanquet’s legacy and data analysis in vegetation science. Journal of Vegetation Science 17, 113–117 (2006).Article 

Google Scholar 
Londo, G. Dezimalskala für die vegetationskundliche Aufnahme von Dauerquadraten. in Sukzessionsforschung (ed. Schmidt, W.). Ber. Int. Symp. Int. Vereinig. Vegetationsk. Rinteln vol. 1973, 613–617 (Cramer, 1975).Bruelheide, H. & Luginbühl, U. Peeking at ecosystem stability: making use of a natural disturbance experiment to analyze resistance and resilience. Ecology 90, 1314–1325 (2009).Article 
PubMed 

Google Scholar 
Hennekens, S. M. & Schaminée, J. H. J. TURBOVEG, a comprehensive data base management system for vegetation data. J. Veg. Sc. 12, 589–591 (2001).Article 

Google Scholar 
Gaston, K. J. & Curnutt, J. L. The dynamics of abundance-range size relationships. Oikos 81, 38 (1998).Article 

Google Scholar 
Gaston, K. J. et al. Abundance-occupancy relationships. J Appl Ecology 37, 39–59 (2000).Article 

Google Scholar 
Sporbert, M. et al. Testing macroecological abundance patterns: The relationship between local abundance and range size, range position and climatic suitability among European vascular plants. J Biogeogr jbi.13926, https://doi.org/10.1111/jbi.13926 (2020).European Commission. Report on the Conservation Status of Habitat Types and Species as required under Article 17 of the Habitats Directive. https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:52009DC0358 (2009).Poschlod, P. Geschichte der Kulturlandschaft. (Ulmer, 2017).Mcgill, B., Enquist, B., Weiher, E. & Westoby, M. Rebuilding community ecology from functional traits. Trends in Ecology & Evolution 21, 178–185 (2006).Article 

Google Scholar 
Jandt, U. et al. More losses than gains during one century of plant biodiversity change in Germany. Nature https://doi.org/10.1038/s41586-022-05320-w (2022).Schaminée, J. H. J., Hennekens, S. M., Chytrý, M. & Rodwell, J. S. Vegetation-plot data and databases in Europe: an overview. Preslia 81, 173–185 (2009).
Google Scholar 
ESA. Land Cover CCI product user guide ver. 2. Tech. Rep. https://maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdf (2017).Kadmon, R., Farber, O. & Danin, A. Effect of roadside bias on the accuracy of predictive maps produced by bioclimatic models. Ecological Applications 14, 401–413 (2004).Article 

Google Scholar 
Davies, C. E., Moss, D. & Hill, M. O. EUNIS Habitat Classification Revised 2004. 310 https://www.eea.europa.eu/data-and-maps/data/eunis-habitat-classification/documentation/eunis-2004-report.pdf/download (2004).Chytrý, M. et al. EUNIS Habitat Classification: Expert system, characteristic species combinations and distribution maps of European habitats. Appl Veg Sci 23, 648–675 (2020).Article 

Google Scholar 
Bruelheide, H., Tichý, L., Chytrý, M. & Jansen, F. Implementing the formal language of the vegetation classification expert systems (ESy) in the statistical computing environment R. Appl Veg Sci, https://doi.org/10.1111/avsc.12562 (2021).Jandt, U., Bruelheide, H. & ReSurveyGermany Consortium. ReSurvey Germany: vegetation-plot resurvey data from Germany. German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig https://doi.org/10.25829/idiv.3514-0qsq70 (2022).Jansen, F. & Dengler, J. GermanSL – eine universelle taxonomische Referenzliste für Vegetationsdatenbanken. Tuexenia 28, 239–253 (2008).
Google Scholar 
Wisskirchen, R. & Haeupler, H. Standardliste der Farn-und Blütenpflanzen Deutschlands. (Ulmer, 1998).Jansen, F. & Dengler, J. Plant names in vegetation databases–a neglected source of bias. Journal of Vegetation Science 21, 1179–1186 (2010).Article 

Google Scholar 
Fischer, H. S. On the combination of species cover values from different vegetation layers. Applied Vegetation Science 18, 169–170 (2015).Article 

Google Scholar 
Schwabe, A. & Kratochwil, A. Pflanzensoziologische Dauerflächen-Untersuchungen im Bannwald ‘Flüh’ (Südschwarzwald) unter besonderer Berücksichtigung der Weidfeld-Sukzession. Standort.Wald 49, 5–49 (2015).
Google Scholar 
Poschlod, P., Schreiber, K.-F., Mitlacher, K., Römermann, C. & Bernhardt-Römermann, M. Entwicklung der Vegetation und ihre naturschutzfachliche Bewertung. in Landschaftspflege und Naturschutz im Extensivgrünland. 30 Jahre Offenhaltungsversuche Baden-Württemberg (eds. Schreiber, K.-F., Brauckmann, H.-J., Broll, G., Krebs, S. & Poschlod, P.) vol. 97 243–288 (2009).Gonzalez, A. et al. Estimating local biodiversity change: a critique of papers claiming no net loss of local diversity. Ecology 97, 1949–1960 (2016).Article 
PubMed 

Google Scholar  More

Adolescents losing sleepEpidemiological studies in US school students aged 14–18 have shown that declines in mental health mirror reductions in the amount of sleep they are getting. Although it is hard to show a causal link between these changes, the authors of this Comment article argue that ensuring that young people get enough sleep is crucial for them to thrive. Various factors could be contributing to this drop-off in sleep, they say, including the use of digital media before bed, schoolwork pressures and extracurricular activities late in the evening or early in the morning.

Sources: J. M. Twenge et al. Sleep Med. 39, 47–53 (2017)/US CDC YRBSS

Monkeypox trajectoryAlmost six months after the monkeypox virus started to spread globally, vaccination efforts and behavioural changes seem to be containing the current strain — at least in the United States and Europe. The number of cases in these regions peaked in August and is now falling. But the situation could still play out in several ways, as this News story reports. At best, the outbreak might fizzle out over the next few months or years. At worst, the virus could become endemic outside Africa.

Source: WHO

The most valuable soilsThis map shows the regions of the world where the conservation of soil should be prioritized. Soils contain a wealth of biodiversity, such as bacteria, fungi, nematode worms and earthworms. These organisms have important roles in ecosystem processes, such as carbon and nutrient cycling, water storage and supporting plant growth. The authors of a paper in Nature set out to identify global hotspots for conservation by surveying soil biodiversity and ecosystem functions at 615 sites around the world. They found hotspots of biodiversity in temperate and Mediterranean regions and in alpine tundra, whereas hotspots of species uniqueness occurred in the tropics and drylands. More than 70% of the hotspots are not adequately covered by protected areas. More

Manjeri, G., Muhamad, R. & Tan, S. G. Oryctes rhinoceros beetles, an oil palm pest in Malaysia. Annu. Res. Rev. Biol. 4, 3429–3439 (2014).Article 

Google Scholar 
Allou, K., Morin, J. P., Kouassi, P., Nklo, F. H. & Rochat, D. Oryctes monoceros trapping with synthetic pheromone and palm material in Ivory Coast. J. Chem. Ecol. 32, 1743–1754 (2006).Article 
CAS 
PubMed 

Google Scholar 
Alibert, H. Study on the insect pests of oil palm in Dahomey. Rev. Botan. Appl. 18, 745–773 (1936).
Google Scholar 
Catley, A. The coconut rhinoceros beetle Oryctes rhinoceros (L) [Coleoptera: Scarabaeidae: Dynastinae]. PANS Pest Articles News Summar. 15, 18–30 (1969).Article 

Google Scholar 
Fauzana, H., Sutikno, A. & Salbiah, D. Population fluctuations Oryctes rhinoceros L. beetle in plant oil palm (Elaeis guineensis Jacq.) given mulching oil palm empty bunch. Cropsaver Int. J. Trop. Insect Sci. 1, 42–47 (2018).
Google Scholar 
Paudel, S., Mansfield, S., Villamizar, L. F., Jackson, T. A. & Marshall, S. D. Can biological control overcome the threat from newly invasive coconut rhinoceros beetle populations (Coleoptera: Scarabaeidae)? A review. Ann. Entomol. Soc. Am. 114, 247–256 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Molet, T. In CPHST Pest Datasheet for Oryctes rhinoceros. USDA-APHIS-PPQCPHST. Revised July 2014 (2013).Hinckley, A. D. Ecology of the coconut rhinoceros beetle, Oryctes rhinoceros (L.) (Coleoptera: Dynastidae). Biotropica 1973, 111–116 (1973).Article 

Google Scholar 
Sitepu, D., Kharie, S., Waroka, JS & Motulo, HFJ. Methods for the production and use of Marhizium anisopliae against Oryctes rhinoceros. In Integrated Coconut Pest Control Project—Annual report of Coconut Research Institute—Manado, North Sulawesi, Indonesia 104–111 (1988).Philippe, R. & Dery, S. K. Coconut research and development. CORD 20, 43–51 (2004).
Google Scholar 
Purrini, K. Baculovirus oryctes release into Oryctes monoceros population in Tanzania, with special reference to the interaction of virus isolates used in our laboratory infection experiments. J. Invertebr. Pathol. 53, 285–300 (1989).Article 

Google Scholar 
Ukeh, D. A., Usua, E. J. & Umoetok, S. B. A. Notes on the biology of Oryctes monoceros (OLIV.) A pest of palms in Nigeria. World J. Agric. Res. 2, 33–36 (2003).
Google Scholar 
Dry, F. W. Notes on the coconut beetle (Oryctes monoceros, Ol.) in Kenya Colony. Bull. Entomol. Res. 13, 103–107 (1922).Article 

Google Scholar 
Bedford, G. O. Biology, ecology, and control of palm rhinoceros beetles. Annu. Rev. Entomol. 25, 309–339 (1980).Article 

Google Scholar 
Khoo, K. C., Yusoff, M. N. M. & Lee, T. W. Pulp and paper of oil palm trunk. In Research Pamphlet No.107: Oil Palm Stem Utilisation, Kuala Lumpur, Malaysia, FRIM 51–65 (1991).Giblin-Davis, R. M. Borers of palms. In Insects on Palms (eds Moore, D. et al.) (CABI Publishing, Wallingford, 2001).
Google Scholar 
Drumoni, A. & Ponchel, Y. Première capture au Yémen d’ Oryctes (Rykanoryctes) monoceros (Olivier, 1789) et confirmation de la présence de cette espèce africaine dans la Péninsule Arabique (Coleoptera, Dynastidae). Entomol. Afr. 15, 25–29 (2010).
Google Scholar 
Lever, R. J. A. W. Pests of the Coconut Palm (Food and Agriculture Organization of the United Nations, Rome, 1969).Moore, A. Rhinoceros beetle pest found in Guam and Saipan. In Pest Alert, Suva, Fiji: Plant Protection Service, Secretariat of the Pacific Community (2007).Zhang, K., Yao, L., Meng, J. & Tao, J. Maxent modeling for predicting the potential geographical distribution of two peony species under climate change. Sci. Total Environ. Sci. 634, 1326–1334 (2018).Article 
ADS 
CAS 

Google Scholar 
Ding, F., Fu, J., Jiang, D., Hao, M. & Lin, G. Mapping the spatial distribution of Aedes aegypti and Aedes albopictus. Acta Trop. 178, 155–162 (2018).Article 
PubMed 

Google Scholar 
Valencia-Rodríguez, D., Jiménez-Segura, L., Rogéliz, C. A. & Parra, J. L. Ecological niche modeling as an effective tool to predict the distribution of freshwater organisms: The case of the Sabaleta Brycon henni (Eigenmann, 1913). PLoS ONE 16, e0247876 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Escobar, L. E., Qiao, H., Cabello, J. & Peterson, A. T. Ecological niche modeling re-examined: A case study with the Darwin’s fox. Ecol. Evol. 8, 4757–4770 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Warren, D. L. & Seifert, S. N. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecol. Appl. 21, 335–342 (2011).Article 
PubMed 

Google Scholar 
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 
Phillips, S. J. Transferability, sample selection bias and background data in presence-only modelling: A response to Peterson et al. (2007). Ecography 31, 272–278 (2008).Article 

Google Scholar 
Elith, J. et al. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17, 43–57 (2011).Article 

Google Scholar 
Phillips, S. J. & Dudík, M. Modeling of species distributions with MaxEnt: New extensions and a comprehensive evaluation. Ecography 31, 161–175 (2008).Article 

Google Scholar 
Arnold, J. D., Brewer, S. C. & Dennison, P. E. Modeling climate-fire connections within the Great basin and Upper Colorado River Basin. Fire Ecol. 10, 64–75 (2014).Article 

Google Scholar 
Phillips, J. S. & Elith, J. On estimating probability of presence from use-availability or presence-background data. Ecology 94, 1409–1419 (2013).Article 
PubMed 

Google Scholar 
Santana, P. A. Jr., Kumar, L., Da Silva, R. S., Pereira, J. L. & Picanço, M. C. Assessing the impact of climate change on the worldwide distribution of Dalbulus maidis (DeLong) using MaxEnt. Pest. Manag. Sci. 75, 2706–2715 (2019).Article 
CAS 
PubMed 

Google Scholar 
Li, et al. Predicting the current and future distributions of Brontispa longissima (Coleoptera: Chrysomelidae) under climate change in China. Glob. Ecol. Conserv. 25, e01444 (2021).Article 

Google Scholar 
Li, T. et al. Direct and indirect effects of environmental factors, spatial constraints, and functional traits on shaping the plant diversity of montane forests. Ecol. Evol. 10, 557–568 (2020).Article 
PubMed 

Google Scholar 
Namgung, H., Kim, M. J., Baek, S., Lee, J. H. & Kim, H. Predicting potential current distribution of Lycorma delicatula (Hemiptera: Fulgoridae) using MaxEnt model in South Korea. J. Asia Pac. Entomol. 23, 291–297 (2020).Article 

Google Scholar 
Ji, W., Gao, G. & Wei, J. Potential global distribution of Daktulosphaira vitifoliae under climate change based on MaxEnt. Insects. 12, 347 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Ji, W., Han, K., Lu, Y. & Wei, J. Predicting the potential distribution of the vine mealybug, Planococcus ficus under climate change by MaxEnt. J. Crop. Prot. 137, 105268 (2020).Article 

Google Scholar 
Sharma, HC & Prabhakar, CS. Impact of climate change on pest management and food security. In Integrated Pest Management 23–36 (Academic Press, Cambridge, 2014).Skendžić, S., Zovko, M., Živković, I. P., Lešic, V. & Lemić, D. The impact of climate change on agricultural insect pests. Insects. 12, 440 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Ward, N. L. & Masters, G. J. Linking climate change and species invasion: An illustration using insect herbivores. Glob. Change Biol. 13, 1605–1615 (2007).Article 
ADS 

Google Scholar 
De Queiroz, D. L., Burckhardt, D. & Majer, J. Integrated pest management of eucalypt psyllids (Insecta, Hemiptera, Psylloidea). In Integrated pest management and pest control-current and future tactics. INTECH 2012, 385–412 (2012).
Google Scholar 
Hochberg, M. E. & Waage, J. K. A model for the biological control of Oryctes rhinoceros (Coleoptera: Scarabaeidae) by means of pathogens. J. Appl. Ecol. 28, 514–531 (1991).Article 

Google Scholar 
Liu, Y. et al. MaxEnt modelling for predicting the potential distribution of a near threatened rosewood species (Dalbergia cultrata Graham ex Benth). Ecol. Eng. 141, 105612 (2019).Article 

Google Scholar 
Wang, R. et al. Predictions of potential geographical distribution of Diaphorina citri (Kuwayama) in China under climate change scenarios. Sci. Rep. 10, 1–9 (2020).CAS 

Google Scholar 
Wood, B. J. Studies on the effect of ground vegetation on infestations of Oryctes rhinoceros (L.) (Col., Dynastidae) in young oil palm replantings in Malaysia. Bull Entomol. Res. 59, 85–96 (1969).Article 

Google Scholar 
Mittal, I. C. Survey of scarabaeid (Coleoptera) fauna of Himachal Pradesh (India). J. Entomol. Res. 24, 259–269 (2000).
Google Scholar 
Zheng, C., Jiang, D., Ding, F., Fu, J. & Hao, M. Spatiotemporal patterns and risk factors for scrub typhus from 2007 to 2017 in southern China. Clin. Infect. Dis. 69, 1205–1211 (2019).Article 
PubMed 

Google Scholar 
Chen, S., Ding, F., Hao, M. & Jiang, D. Mapping the potential global distribution of red imported fire ant (Solenopsis invicta Buren) based on a machine learning method. Sustainability. 12, 10182 (2020).Article 

Google Scholar 
Ding, F. et al. Infection and risk factors of human and avian influenza in pigs in south China. Prev. Vet. Med. 190, 105317 (2021).Article 
PubMed 

Google Scholar 
Jiang, D. et al. Spatiotemporal patterns and spatial risk factors for Visceral leishmaniasis from 2007 to 2017 in Western and Central China: A modelling analysis. Sci. Total Environ Sci. 764, 144275 (2021).Article 
ADS 
CAS 

Google Scholar 
Méndez-Rojas, D. M., Cultid-Medina, C. & Escobar, F. Influence of land use change on rove beetle diversity: A systematic review and global meta-analysis of a mega-diverse insect group. Ecol. Indic. 122, 107239 (2021).Article 

Google Scholar 
Oke, T. R. City size and the urban heat island. Atmos. Environ. 7, 769–779 (1973).Article 
ADS 

Google Scholar 
Briere, J. F., Pracros, P., Le Roux, A. Y. & Pierre, J. S. A novel rate model of temperature-dependent development for arthropods. Environ. Entomol. 28, 22–29 (1999).Article 

Google Scholar 
Zeng, Y., Low, B. W. & Yeo, D. C. Novel methods to select environmental variables in MaxEnt: A case study using invasive crayfish. Eco. Model. 341, 5–13 (2016).Article 

Google Scholar 
Fand, B. B. et al. Invasion risk of the South American tomato pinworm Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae) in India: Predictions based on MaxEnt ecological niche modelling. Int. J. Trop. Insect Sci. 40, 1–11 (2020).Article 

Google Scholar 
Li, W. J. et al. Potential distribution prediction of natural Pseudotsuga sinensis forest in Guizhou based on Maxent model. J. For. Res. 48, 47–52 (2019).
Google Scholar 
McIntyre, S., Rangel, E. F., Ready, P. D. & Carvalho, B. M. Species-specific ecological niche modelling predicts different range contractions for Lutzomyia intermedia and a related vector of Leishmania braziliensis following climate change in South America. Parasit. Vectors 10, 1–15 (2017).Article 

Google Scholar 
Hao, M. et al. Global potential distribution of Oryctes rhinoceros, as predicted by boosted regression tree model. Glob. Ecol. Conserv. 37, e02175 (2022).Article 

Google Scholar 
Aidoo, O. F. et al. The impact of climate change on potential invasion risk of Oryctes monoceros worldwide. Front. Ecol. Evol. 10, 633 (2022).Article 

Google Scholar 
Aidoo, O. F. et al. Lethal yellowing disease: Insights from predicting potential distribution under different climate change scenarios. J. Plant Dis. Prot. 2021, 1–13 (2021).
Google Scholar 
Ruheili, A. M. A., Boluwade, A. & Subhi, A. M. A. Assessing the Impact of Climate Change on the Distribution of Lime (16srii-B) and Alfalfa (16srii-D) Phytoplasma Disease Using MaxEnt. Plants. 10, 460 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Wang, R. et al. Predicting the potential distribution of the Asian citrus psyllid, Diaphorina citri (Kuwayama), in China using the MaxEnt model. PeerJ 7, e7323 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
He, S. T. & Jing, P. F. Prediction of potential distribution areas of Salvia bowleyana Dunn. in China based on MaxEnt and suitability analysis. J Anhui Agri. Sci. 8, 2311–2314 (2014).
Google Scholar 
Chahouki, M. A. Z. & Sahragard, H. P. Maxent modelling for distribution of plant species habitats of rangelands (Iran). Pol. J. Ecol. 64, 453–467 (2016).
Google Scholar 
Shabani, F., Kumar, L. & Ahmadi, M. Assessing accuracy methods of species distribution models: AUC, specificity, sensitivity and the true skill statistic. Glob. Int. J. Hum. Soc. Sci. 18, 6–18 (2018).
Google Scholar 
Baloch, M. N., Fan, J., Haseeb, M. & Zhang, R. Mapping potential distribution of Spodoptera frugiperda (Lepidoptera: Noctuidae) in central Asia. Insects. 11, 172 (2020).Article 
PubMed Central 

Google Scholar 
Wang, N., Li, Z., Wu, J., Rajotte, E. G., Wan, F & Wang, Z. The potential geographical distribution of Bactrocera dorsalis (Diptera: Tephrididae) in China based on emergence rate model and ArcGIS. In International Conference on Computer and Computing Technologies in Agriculture 399–411. (Springer, Boston, 2008).Manrique, V., Cuda, J. P., Overholt, W. A. & Diaz, R. Temperature-dependent development and potential distribution of Episimus utilis (Lepidoptera: Tortricidae), a candidate biological control agent of Brazilian peppertree (Sapindales: Anacardiaceae) in Florida. Environ. Entomol. 37, 862–870 (2008).Article 
PubMed 

Google Scholar 
Das, D. K., Singh, J. & Vennila, S. Emerging crop pest scenario under the impact of climate change–a brief review. AgroPhysics. 11, 13–20 (2011).CAS 

Google Scholar 
Porter, J. H., Parry, M. L. & Carter, T. R. The potential effects of climatic change on agricultural insect pests. Agric. For. Meteorol. 57, 221–240 (1991).Article 
ADS 

Google Scholar 
Trenberth, K. E. Climate change caused by human activities is happening and it already has major consequences. J. Energy Nat. Resour. Law. 36, 463–481 (2018).Article 

Google Scholar 
Xu, D., Zhuo, Z., Li, X. & Wang, R. Distribution and invasion risk assessment of Oryctes rhinoceros (L.) in China under changing climate. J. Appl. Entomol. 146, 385–395 (2022).Article 

Google Scholar 
Sushil, K. & Mukhtar, A. Effect of temperature and humidity on biology of rhinoceros beetle, Oryctes rhinoceros Linn. on oil palm. J. Appl. Anim. Res. 18, 108–112 (2007).
Google Scholar 
Sabidin, N. N. E. The effect of climate change to the population of rhinoceros beetle (Oryctes rhinoceros) at selected oil palm plantation. In Bachelor of Science Thesis Dissertation. Universiti Teknologi MARA. https://ir.uitm.edu.my/id/eprint/22754. (2018).Yadav, R. & Chang, N. T. Effects of temperature on the development and population growth of the melon thrips, Thrips palmi, on eggplant, Solanum melongena. J. Insect Sci. 14, 78 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Ju, R. T., Wang, F. & Li, B. Effects of temperature on the development and population growth of the sycamore lace bug, Corythucha ciliata. J. Insect Sci. 11, 1–12 (2011).Article 

Google Scholar 
Zheng, F. S., Du, Y. Z., Wang, Z. J. & Xu, J. J. Effect of temperature on the demography of Galerucella birmanica (Coleoptera: Chrysomelidae). Insect Sci. 15, 375–380 (2008).Article 

Google Scholar 
Azrag, A. G. et al. Modelling the effect of temperature on the biology and demographic parameters of the African coffee white stem borer, Monochamus leuconotus (Pascoe) (Coleoptera: Cerambycidae). J. Therm. Biol. 89, 102534 (2020).Article 
CAS 
PubMed 

Google Scholar 
Aidoo, O. F. et al. The African citrus triozid Trioza erytreae Del Guercio (Hemiptera: Triozidae): Temporal dynamics and susceptibility to entomopathogenic fungi in East Africa. Int. J. Trop. Insect Sci. 41, 563–573 (2021).Article 

Google Scholar 
Leonard, A. et al. Predicting the current and future distribution of the edible long-horned grasshopper Ruspolia differens (Serville) using temperature-dependent phenology models. J. Therm. Biol. 95, 102786 (2021).Article 
PubMed 

Google Scholar 
Roy, B. A. et al. Increasing forest loss worldwide from invasive pests requires new trade regulations. Front. Ecol. Environ. 12, 457–465 (2014).Article 

Google Scholar 
Shabani, F., Kumar, L. & Ahmadi, M. A comparison of absolute performance of different correlative and mechanistic species distribution models in an independent area. Ecol. Evol. 6, 5973–5986 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Cianci, D., Hartemink, N. & Ibáñez-Justicia, A. Modelling the potential spatial distribution of mosquito species using three different techniques. Int. J. Health Geogr. 14, 10 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Zelazny, B. & Alfiler, A. Oryctes rhinoceros (Coleoptera: Scarabaeidae) larva abundance and mortality factors in the Philippines. Environ. Entomol. 15, 84–87 (1986).Article 

Google Scholar 
Wood, B.J. Studies on the effect of ground vegetation on infestations of Oryctes rhinoceros (L.)(Col., Dynastidae) in young oil palm replantings in Malaysia. Bull. Entomol. Res. 59, 85–96 (1969). More

Asian elephants spend most of their time outside protected areas because they prefer the food they find there, an international team of scientists reports. But this behaviour is putting the animals and people in harm’s way, say researchers.The finding has important implications for the long-term survival of the animals because protected areas are a cornerstone of global conservation strategies to protect threatened species, say researchers.If protected areas do not contain animals’ preferred habitats, they will wander out, says Ahimsa Campos-Arceiz, who studies Asian elephants (Elephas maximus) at the Chinese Academy of Sciences’ Xishuangbanna Tropical Botanical Garden in Menglun, China. “It’s a good intention, but doesn’t always work out that way.”Human–elephant conflict is the biggest threat for Asian elephants. Over the past few decades, animals in protected areas have increasingly wandered into villages. They often cause destruction, damaging crops and infrastructure and injuring and even killing people.Wandering elephantsTo understand how effective protected areas are for conserving Asian elephants, Campos-Arceiz and his colleagues set out to get a precise picture of Asian-elephant movements. They collared 102 individuals in Peninsular Malaysia and Borneo, recording 600,000 GPS locations over a decade. They found that most elephants spent most of their time in habitats outside the protected areas, at the forest edge and in areas of regrowth. The findings were published in the Journal of Applied Ecology1 on 18 October.The researchers suspect that the elephants venture out because they like to eat grasses, bamboo, palms and fast-growing trees, which are common in disturbed forests and relatively scarce under the canopy of old-growth forests.Philip Nyhus, a conservation biologist who specializes in human–wildlife conflict at Colby College in Waterville, Maine, says Asian elephants live deep in dense forest and so are much more difficult to study than African elephants, which roam open savannahs. “The sample size is impressive,” he says.The finding is not unexpected given past anecdotal observations of elephant behaviour, says Nyhus. But now the data show that this is a common strategy for the survival of these animals, and not just something seen in a subset of the population. The research provides strong evidence for how to set up suitable protected areas that reduce the risk of elephants wandering out, he says.‘There will be conflict’The results do not diminish the importance of protected areas, which provide long-term safety for the animals, says Campos-Arceiz, who did the field work while at the University of Nottingham Malaysia in Selangor. “But they are clearly not enough.”The study suggests that “there will be conflict between humans and elephants”, says Guo Xianming, director of the Research Institute of Xishuangbanna National Nature Reserve in Jinghong.Asian elephants wander into villages owing to a combination of reasons: an increase in elephant populations, forests in many reserves have grown denser and have become unsuitable for the animals, and increasing habitat loss and degradation outside.Last year, two herds of elephants made global headlines as they wandered out of the Xishuangbanna National Nature Reserve and travelled for hundreds of kilometers, wreaking havoc along the way. One herd spent five weeks at the botanical garden where Campos-Arceiz works. “It was intense,” he says.There is an urgent need to understand how people and elephants can better share the landscape, says Guo. And the first step is by better protecting people’s lives and livelihoods. “It’s the only way of peaceful co-existence.”
The reporting of the story was supported by International Women’s Media Foundation’s Howard G. Buffett Fund for Women Journalists. More