Ecology

Gruninger, R. J., Ribeiro, G. O., Cameron, A. & McAllister, T. A. Invited review: Application of meta-omics to understand the dynamic nature of the rumen microbiome and how it responds to diet in ruminants. Animal 13, 1843–1854 (2019).CAS 

Google Scholar 
Morgavi, D. P., Kelly, W. J., Janssen, P. H. & Attwood, G. T. Rumen microbial (meta)genomics and its application to ruminant production. Animal 7, 184–201 (2013).CAS 

Google Scholar 
Bergman, E. N. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol. Rev. 70, 567–590 (1990).CAS 

Google Scholar 
Flint, H. J. The rumen microbial ecosystem—Some recent developments. Trends Microbiol. 5, 483–488 (1997).CAS 

Google Scholar 
Hobson, P. N. & Stewart, C. S. The Rumen Microbial Ecosystem. (Springer, 2012).Moraïs, S. & Mizrahi, I. The road not taken: The rumen microbiome, functional groups, and community states. Trends Microbiol. 27, 538–549 (2019).
Google Scholar 
Cheng, K. J., Forsberg, C. W., Minato, H. & Costerton, J. W. in Physiological Aspects of Digestion and Metabolism in Ruminants (eds T. Tsuda, Y. Sasaki, & R. Kawashima) 595–624 (Academic Press, 1991).McSweeney, C. S., Palmer, B., McNeill, D. M. & Krause, D. O. Microbial interactions with tannins: Nutritional consequences for ruminants. Anim. Feed Sci. Technol. 91, 83–93 (2001).CAS 

Google Scholar 
Skene, I. K. & Brooker, J. D. Characterization of tannin acylhydrolase activity in the ruminal bacterium Selenomonas ruminantium. Anaerobe 1, 321–327 (1995).CAS 

Google Scholar 
Khanbabaee, K. & van Ree, T. Tannins: Classification and definition. Nat. Prod. Rep. 18, 641–649 (2001).CAS 

Google Scholar 
Makkar, H. P. S. & Becker, K. Isolation of tannins from leaves of some trees and shrubs and their properties. J. Agric. Food Chem. 42, 731–734 (1994).CAS 

Google Scholar 
Bhat, T. K., Kannan, A., Singh, B. & Sharma, O. P. Value addition of feed and fodder by alleviating the antinutritional effects of tannins. Agr. Res. 2, 189–206 (2013).CAS 

Google Scholar 
Shimada, T. Salivary proteins as a defense against dietary tannins. J. Chem. Ecol. 32, 1149–1163 (2006).CAS 

Google Scholar 
Zhu, J., Filippich, L. J. & Alsalami, M. T. Tannic acid intoxication in sheep and mice. Res. Vet. Sci. 53, 280–292 (1992).CAS 

Google Scholar 
Kohl, K. D., Stengel, A. & Dearing, M. D. Inoculation of tannin-degrading bacteria into novel hosts increases performance on tannin-rich diets. Environ. Microbiol. 18, 1720–1729 (2016).CAS 

Google Scholar 
Kumar, K., Chaudhary, L. C., Agarwal, N. & Kamra, D. N. Isolation and characterization of tannin-degrading bacteria from the rumen of goats fed oak (Quercus semicarpifolia) leaves. Agr. Res. 3, 377–385 (2014).
Google Scholar 
Odenyo, A. A. et al. Characterization of tannin-tolerant bacterial isolates from East African ruminants. Anaerobe 7, 5–15 (2001).CAS 

Google Scholar 
Grilli, D. J. et al. Analysis of the rumen bacterial diversity of goats during shift from forage to concentrate diet. Anaerobe 42, 17–26 (2016).
Google Scholar 
Tong, J. et al. Illumina sequencing analysis of the ruminal microbiota in high-yield and low-yield lactating dairy cows. PLoS ONE 13, e0198225 (2018).
Google Scholar 
Pope, P. B. et al. Metagenomics of the Svalbard reindeer rumen microbiome reveals abundance of polysaccharide utilization loci. PLoS ONE 7, e38571 (2012).ADS 
CAS 

Google Scholar 
Østbye, K., Wilson, R. & Rudi, K. Rumen microbiota for wild boreal cervids living in the same habitat. FEMS Microbiol. Lett. 363, fnw233 (2016).
Google Scholar 
Gruninger, R. J., Sensen, C. W., McAllister, T. A. & Forster, R. J. Diversity of rumen bacteria in Canadian cervids. PLoS ONE 9, e89682 (2014).ADS 

Google Scholar 
Henderson, G. et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5, 14567 (2015).CAS 

Google Scholar 
Reese, A. T. & Kearney, S. M. Incorporating functional trade-offs into studies of the gut microbiota. Curr. Opin. Microbiol. 50, 20–27 (2019).CAS 

Google Scholar 
Moeller, A. H. et al. Social behavior shapes the chimpanzee pan-microbiome. Sci. Adv. 2, e1500997 (2016).ADS 

Google Scholar 
Okano, T. & Matsuda, H. Biocultural diversity of Yakushima Island: Mountains, beaches, and sea. J. Mar. Isl. Cult. 2, 69–77 (2013).
Google Scholar 
Agetsuma, N., Agetsuma-Yanagihara, Y. & Takafumi, H. Food habits of Japanese deer in an evergreen forest: Litter-feeding deer. Mamm. Biol. 76, 201–207 (2011).
Google Scholar 
Higashi, Y., Hirota, S. K., Suyama, Y. & Yahara, T. Geographical and seasonal variation of plant taxa detected in faces of Cervus nippon yakushimae based on plant DNA analysis in Yakushima Island. Ecol. Res. 37, 582–597 (2022).CAS 

Google Scholar 
Kuroiwa, A. Nutritional ecology of the Yakushika (Cervus nippon yakushimae) population under high density Ph.D. thesis, Kyushu University, (2017).Koda, R., Agetsuma, N., Agetsuma-Yanagihara, Y., Tsujino, R. & Fujita, N. A proposal of the method of deer density estimate without fecal decomposition rate: A case study of fecal accumulation rate technique in Japan. Ecol. Res. 26, 227–231 (2011).
Google Scholar 
Yahara, T. in Deer eats world heritages: Ecology of deer and forets (eds T. Yumoto & H. Matsuda) 168–187 (Bunichi-Sogo-Shuppan, 2006).Onoda, Y. & Yahara, T. in Challenges for Conservation Ecology in Space and Time. (eds T. Miyashita & J. Nishihiro) 126–149 (University of Tokyo Press, 2015).Kagoshima Prefecture Nature Conservation Division. The current status of Yakusika in FY 2020, available at https://www.rinya.maff.go.jp/kyusyu/fukyu/shika/attach/pdf/yakushikaWG_R3_2-23.pdf (2020).Kuroiwa, A., Kuroe, M. & Yahara, T. Effects of density, season, and food intake on sika deer nutrition on Yakushima Island, Japan. Ecol. Res. 32, 369–378 (2017).
Google Scholar 
Hiura, T., Hashidoko, Y., Kobayashi, Y. & Tahara, S. Effective degradation of tannic acid by immobilized rumen microbes of a sika deer (Cervus nippon yesoensis) in winter. Anim. Feed Sci. Technol. 155, 1–8 (2010).CAS 

Google Scholar 
Kawarai, S. et al. Seasonal and geographical differences in the ruminal microbial and chloroplast composition of sika deer (Cervus nippon) in Japan. Sci. Rep. 12, 6356 (2022).ADS 
CAS 

Google Scholar 
Li, Z. et al. Response of the rumen microbiota of sika deer Cervus nippon fed different concentrations of tannin rich plants. PLoS ONE 10, e0123481 (2015).
Google Scholar 
McDonald, D. et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618 (2012).CAS 

Google Scholar 
Kim, M., Morrison, M. & Yu, Z. Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol. Ecol. 76, 49–63 (2011).CAS 

Google Scholar 
Weimer, P. J. Redundancy, resilience, and host specificity of the ruminal microbiota: Implications for engineering improved ruminal fermentations. Front. Microbiol. 6, 296 (2015).
Google Scholar 
Scott, K. P., Gratz, S. W., Sheridan, P. O., Flint, H. J. & Duncan, S. H. The influence of diet on the gut microbiota. Pharmacol. Res. 69, 52–60 (2013).CAS 

Google Scholar 
Tapio, I. et al. Taxon abundance, diversity, co-occurrence and network analysis of the ruminal microbiota in response to dietary changes in dairy cows. PLoS ONE 12, e0180260 (2017).
Google Scholar 
Ohara, M. in Agriculture in Hokkaido v2 (ed K. Iwama, Ohara, M., Araki, H., Yamada, T., Nakatsuji, H., Kataoka, T., Yamamoto, Y.) 1–18(Faculty of Agriculture, Hokkaido Univ., 2009).Igota, H., Sakuragi, M. & Uno, H. in Sika Deer: Biology and Management of Native and Introduced Populations (eds. Dale R. McCullough, Seiki Takatsuki, & Koichi Kaji) 251–272 (Springer Japan, 2009).Fernando, S. C. et al. Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ. Microbiol. 76, 7482–7490 (2010).ADS 
CAS 

Google Scholar 
Hu, X. et al. High-throughput analysis reveals seasonal variation of the gut microbiota composition within forest musk deer (Moschus berezovskii). Front. Microbiol. 9, (2018).Artzi, L., Morag, E., Shamshoum, M. & Bayer, E. A. Cellulosomal expansion: Functionality and incorporation into the complex. Biotechnol. Biofuels 9, 61 (2016).
Google Scholar 
Biddle, A., Stewart, L., Blanchard, J. & Leschine, S. Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in diverse gut communities. Diversity 5, (2013).Eisenhauer, N., Scheu, S. & Jousset, A. Bacterial diversity stabilizes community productivity. PLoS ONE 7, e34517 (2012).ADS 
CAS 

Google Scholar 
Miller, A. W., Oakeson, K. F., Dale, C. & Dearing, M. D. Effect of dietary oxalate on the gut microbiota of the mammalian herbivore Neotoma albigula. Appl. Environ. Microbiol. 82, 2669–2675 (2016).ADS 
CAS 

Google Scholar 
Adams, J. M., Rehill, B., Zhang, Y. & Gower, J. A test of the latitudinal defense hypothesis: Herbivory, tannins and total phenolics in four North American tree species. Ecol. Res. 24, 697–704 (2009).CAS 

Google Scholar 
Nabeshima, E., Murakami, M. & Hiura, T. Effects of herbivory and light conditions on induced defense in Quercus crispula. J. Plant Res. 114, 403–409 (2001).
Google Scholar 
Yang, C.-M., Yang, M.-M., Hsu, J.-M. & Jane, W.-N. Herbivorous insect causes deficiency of pigment–protein complexes in an oval-pointed cecidomyiid gall of Machilus thunbergii leaf. Bot. Bull. Acad. Sin. 44, 315–321 (2003).
Google Scholar 
Agetsuma, N., Agetsuma-Yanagihara, Y., Takafumi, H. & Nakaji, T. Plant constituents affecting food selection by sika deer. J. Wildl. Manag. 83, 669–678 (2019).
Google Scholar 
Couch, C. E. et al. Diet and gut microbiome enterotype are associated at the population level in African buffalo. Nat. Commun. 12, 2267 (2021).ADS 
CAS 

Google Scholar 
Goel, G., Puniya, A. K. & Singh, K. Tannic acid resistance in ruminal streptococcal isolates. J. Basic Microbiol. 45, 243–245 (2005).CAS 

Google Scholar 
Jiménez, N. et al. Genetic and biochemical approaches towards unravelling the degradation of gallotannins by Streptococcus gallolyticus. Microb. Cell Fact. 13, 154 (2014).
Google Scholar 
Nelson, K. E., Thonney, M. L., Woolston, T. K., Zinder, S. H. & Pell, A. N. Phenotypic and phylogenetic characterization of ruminal tannin-tolerant bacteria. Appl. Environ. Microbiol. 64, 3824–3830 (1998).ADS 
CAS 

Google Scholar 
Selwal, M. K. et al. Optimization of cultural conditions for tannase production by Pseudomonas aeruginosa IIIB 8914 under submerged fermentation. World J. Microbiol. Biotechnol. 26, 599–605 (2010).CAS 

Google Scholar 
Kohl, K. D., Weiss, R. B., Cox, J., Dale, C. & Denise Dearing, M. Gut microbes of mammalian herbivores facilitate intake of plant toxins. Ecol. Lett. 17, 1238–1246 (2014).
Google Scholar 
Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Method 7, 335–336 (2010).CAS 

Google Scholar 
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200 (2011).CAS 

Google Scholar 
Caporaso, J. G. et al. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26, 266–267 (2009).
Google Scholar 
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).ADS 

Google Scholar 
R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria, 2020).Osawa, R. Formation of a clear zone on tannin-treated brain heart infusion agar by a Streptococcus sp. isolated from feces of koalas. Appl. Environ. Microbiol. 56, 829–831 (1990).ADS 
CAS 

Google Scholar 
Hamamura, N., Olson, S. H., Ward, D. M. & Inskeep, W. P. Diversity and functional analysis of bacterial communities associated with natural hydrocarbon seeps in acidic soils at Rainbow Springs, Yellowstone National Park. Appl. Environ. Microbiol. 71, 5943–5950 (2005).ADS 
CAS 

Google Scholar 
Benson, D. A. et al. GenBank. Nucleic Acids Res. 41, D36–D42 (2012).ADS 

Google Scholar 
Chen, I.-M. A. et al. The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities. Nucleic Acids Res. 49, D751–D763 (2020)Suzuki, M. T., Taylor, L. T. & Delong, E. F. Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5 ’-nuclease assays. Appl. Environ. Microbiol. 66, 4605–4614 (2000).ADS 
CAS 

Google Scholar  More

Study area and soil collectionTwenty NRCS map units were selected across Hawaii Commercial & Sugar Company (HC&S) in central Maui that represented seven soil orders, 10 NRCS soil series, and approximately 77% of the total plantation area (Fig. 1). Soil heterogeneity across the landscape allowed for the comparison of a continuum of soil and soil properties that have experienced the same C4 grass inputs and agricultural treatment under sugarcane production for over 100 years. Conventional sugarcane production involved 2-year growth followed by harvest burn, collection of remaining stalks by mechanical ripper, deep tillage to 40 cm, no crop rotations, and little to no residue return. The sampled soils, collected from September-August 2015, thus represent a baseline of SOC after input-intensive tropical agriculture and long-term soil disturbance. Elemental analyses from this work show consistent agricultural disturbances led to degraded SOC content ranging from 0.23 to 2.91% SOC of soil mass with an average of only 1.16% SOC across all locations and depths.Figure 1Hawaiian Commercial and Sugar in central Maui with main Hawaiian Islands inset (left). Soil series identified by NRCS across HC&S fields (right) with black dots indicating 20 locations where soils were sampled to test landscape level differences in topical soil kinetics and associated soil properties under conventional sugarcane. Maps from Ref.19 created using ESRI ArcGIS with soil series data from: Soil Survey Staff, Natural Resources Conservation Service, United States Department of Agriculture, Web Soil Survey, Available online at http://websoilsurvey.nrcs.usda.gov/. Accessed [07/30/2016]19.Full size imageThe homogenized land use history allowed focused investigation of soil property effects on SOC storage across heterogenous soils (Table 1). Though soil inputs (e.g. water, nutrients, root inputs, residue removal) and disturbance regimes (e.g. burn, rip, till, compaction, no crop rotation) were consistent across the 20 field locations, average annual surface temperatures varied from 22.9 to 25.1 °C with a mean of 24.4 °C, average annual relative humidity varied from 70.4 to 79.2% with a mean of 73.4%, and average annual rainfall varied from 306 to 1493 mm with a mean of 575 mm. Gradients of rainfall, relative humidity, and elevation across the site generally increase in an east/north-east direction towards the prevailing winds and up the western slope of Haleakalā. In contrast, surface temperatures increase in the opposite direction towards Kihei and the southern tip of the West Maui Mountains.Table 1 NRCS soil classification and environmental conditions at 20 field sites.Full size tableaSoil descriptions26.bInterpolated estimates from Ref.25.Soil sampling and analysisPit locations were identified with a handheld GPS and were sampled using NRCS Rapid Carbon Assessment methods27. A total of 75 horizons were identified from the 20 selected map units to a depth of 1 m28,29. The central depth of each horizon was sampled using volumetric bulk density cores up to 50 cm. After 50 cm, a hand auger was used to check for any further horizon changes. The bulk density of horizons past 50 cm were estimated using collected soil mass and known auger size. Collected soils were air dried, processed through a 2 mm sieve, and analyzed for total C and nitrogen percent, SOC percent, soil texture, iron (Fe) and aluminum (Al) minerals, pH, cation and anion exchange capacity, extractable cations, wet and dry size classes, aggregate stability, and soil water potential at -15 kPa. Total C and nitrogen were measured by elemental analysis (Costech, ECS 4010, Valencia, CA), with SOC content determined by elemental analysis after hydrochloric acid digestion to remove carbonates. Soil texture was measured using sedimentary separation, while a 10:1 soil slurry in water was used to test soil pH. Soil pressure plates were used to measure soil water potential at -15 kPa.Fe and Al oxides were quantified in mineral phases using selective dissolutions of collected soils, including: (1) a 20:1 sodium citrate to sodium dithionite extraction, shaken 16 h, to quantify total free minerals30, (2) 0.25 M hydroxylamine hydrochloride and hydrochloric acid extraction, shaken 16 h, to quantify amorphous minerals31, and (3) 0.1 M sodium pyrophosphate (pH 10), shaken 16 h and centrifuged at 20,000g, to quantify organo-bound metals30. Extracted Fe, Al, and Si from al extractions were measured by inductively coupled plasma analysis (PerkinElmer, Optima ICP-OES, Norwalk, CT). Exploratory ratios of Fe/Al, Fe/Si, and Al/Si for the citrate/dithionite (c), hydroxylamine (h), and pyrophosphate (p) extractions were calculated. Crystalline Fe, operationally-defined as the difference between the citrate dithionite and hydroxylamine extraction, and Al + ½ Fe32 were calculated for each extraction.Plant-available phosphorus was extracted by the Olsen method using 0.5 M sodium bicarbonate adjusted to pH 8.5 and measured by continuous flow colorimetry (Hach, Lachat Quickchem 8500, Loveland, CO). Exchangeable cations (i.e. calcium, magnesium, potassium, and sodium), effective cation exchange capacity, and anion exchange capacity were measured by compulsive exchange using barium chloride and magnesium sulfate33. Cations were quantified by continuous flow colorimetry and flame-spectroscopy (Hach, Lachat Quickchem 8500, Loveland, CO). Field soils were air dried and initially passed through a 2 mm sieve before size classes of macroaggregate (2 mm – 250 µm) and microaggregate ( More

Human-modified forests, such as selectively logged forests, are often characterized as degraded ecosystems because of their altered structure and low biomass. The concept of ecosystem degradation can be a double-edged sword. It rightly draws attention to the conservation value of old-growth systems and the importance of ecosystem restoration. However, it can also suggest that human-modified ecosystems are of low ecological value and therefore, in some cases, suitable for conversion to agriculture (such as oil palm plantations) and other land uses3,4,5.Selectively logged and other forms of structurally altered forests are becoming the prevailing vegetation cover in much of the tropical forest biome2. Such disturbance frequently leads to a decline in old-growth specialist species1, and also in non-specialist species in some contexts6,7,8. However, species-focused biodiversity metrics are only one measure of ecosystem vitality and functionality, and rarely consider the collective role that suites of species play in maintaining ecological functions9.An alternative approach is to focus on the energetics of key taxonomic groups, and the number and relative dominance of species contributing to each energetic pathway. Energetic approaches to examining ecosystem structure and function have a long history in ecosystem ecology10. Virtually all ecosystems are powered by a cascade of captured sunlight through an array of autotroph tissues and into hierarchical assemblages of herbivores, carnivores and detritivores. Energetic approaches shine light on the relative significance of energy flows among key taxa and provide insight into the processes that shape biodiversity and ecosystem function. The common currency of energy enables diverse guilds and taxa to be compared in a unified and physically meaningful manner: dominant energetic pathways can be identified, and the resilience of each pathway to the loss of individual species can be assessed. Quantitative links can then be made between animal communities and the plant-based ecosystem productivity on which they depend. The magnitude of energetic pathways in particular animal groups can often be indicators of key associated ecosystem processes, such as nutrient cycling, seed dispersal and pollination, or trophic factors such as intensity of predation pressure or availability of resource supply, all unified under the common metric of energy flux11,12.Energetics approaches have rarely been applied in biodiverse tropical ecosystems because of the range of observations they require11,12,13. Such analyses rely on: population density estimates for a very large number of species; understanding of the diet and feeding behaviour of the species; and reliable estimation of net primary productivity (NPP). Here we take advantage of uniquely rich datasets to apply an energetics lens to examine and quantify aspects of the ecological function and vitality of habitats in Sabah, Malaysia, that comprise old-growth forests, logged forest and oil palm plantation (Fig. 1 and Extended Data Fig. 1). Our approach is to calculate the short-term equilibrium production or consumption rates of food energy by specific species, guilds or taxonomic groups. We focus on three taxonomic groups (plants, birds and mammals) that are frequently used indicators of biodiversity and are relatively well understood ecologically.Fig. 1: Maps of the study sites in Sabah, Borneo.a–d, Maps showing locations of NPP plots and biodiversity surveys in old-growth forest, logged forest and oil palm plantations in the Stability of Altered Forest Ecosystems Project landscape (a), Maliau Basin (b), Danum Valley (c) and Sepilok (d). The inset in a shows the location of the four sites in Sabah. The shade of green indicates old-growth (dark green), twice-logged (intermediate green) or heavily logged (light green) forests. The camera and trap grid includes cameras and small mammal traps. White areas indicate oil palm plantations.Full size imageWe are interested in the fraction of primary productivity consumed by birds and mammals, and how it varies along the disturbance gradient, and how and why various food energetic pathways in mammals and birds, and the diversity of species contributing to those pathways, vary along the disturbance gradient. To estimate the density of 104 mammal and 144 bird species in each of the three habitat types, we aggregated data from 882 camera sampling locations (a total of 42,877 camera trap nights), 508 bird point count locations, 1,488 small terrestrial mammal trap locations (34,058 live-trap nights) and 336 bat trap locations (Fig. 1 and Extended Data Fig. 1). We then calculated daily energetic expenditure for each species based on their body mass, assigned each species to a dietary group and calculated total food consumption in energy units. For primary productivity, we relied on 34 plot-years (summation of plots multiplied by the number of years each plot is monitored) of measurements of the key components of NPP (canopy litterfall, woody growth, fine root production) using the protocols of the Global Ecosystem Monitoring Network14,15,16 across old-growth (n = 4), logged (n = 5) and oil palm (n = 1) plots. This dataset encompasses more than 14,000 measurements of litterfall, 20,000 tree diameter measurements and 2,700 fine root samples.Overall bird species diversity is maintained across the disturbance gradient and peaks in the logged forest; for mammals, there is also a slight increase in the logged forest, followed by rapid decline in the oil palm (Fig. 2b,c). Strikingly, both bird and mammal biomass increases substantially (144% and 231%, respectively) in the logged forest compared to the old-growth forest, with mammals contributing about 75% of total (bird plus mammal) biomass in both habitat types (Fig. 2b,c).Fig. 2: Variation of ecosystem energetics along the disturbance gradient from old-growth forest through logged forest to oil palm.a, Total NPP along the gradient (mean of intensive 1-ha plots; n = 4 for old growth (OG), n = 5 for logged and n = 1 for oil palm (OP); error bars are 95% confidence intervals derived from propagated uncertainty in the individually measured NPP components), with individual plot data points overlaid. b,c, Total body mass (bars, left axis) and number of species counted (blue dots and line, right axis) of birds (b) and mammals (c). d,e, Total direct energetic food intake by birds (d) and mammals (e). f,g, Percentage of NPP directly consumed by birds (f) and mammals (g). In b–e, body mass and energetics were estimated for individual bird and mammal species, with the bars showing the sum. Error bars denote 95% confidence intervals derived from 10,000 Monte Carlo simulation estimates incorporating uncertainty in body mass, population density, the daily energy expenditure equation, assimilation efficiency of the different food types, composition of the diet of each species and NPP. In f,g, the grey bars indicate direct consumption of NPP, white bars denote the percentage of NPP indirectly supporting bird and mammal food intake when the mean trophic level of consumed invertebrates is assumed to be 2.5, with the error bars denoting assumed mean trophic levels of 2.4 and 2.6. Note the log scale of the y axis in f,g. Numbers for d,e provided in Supplementary Data Tables 1, 2.Full size imageThe total flow of energy through consumption is amplified across all energetic pathways by a factor of 2.5 (2.2–3.0; all ranges reported are 95% confidence intervals) in logged forest relative to old-growth forest. In all three habitat types, total energy intake by birds is much greater than by mammals (Fig. 2d,e and Extended Data Table 1). Birds account for 67%, 68% and 90% of the total direct consumption by birds and mammals combined in old-growth forests, logged forests and oil palm, respectively. Although mammal biomass is higher than bird biomass in the old-growth and logged forests, the metabolism per unit mass is much higher in birds because of their small body size; hence, in terms of the energetics and consumption rates, the bird community dominates. The total energy intake by birds alone increases by a factor of 2.6 (2.1–3.2) in the logged forest relative to old-growth forest. This is mainly driven by a 2.5-fold (1.7–2.8) increase in foliage-gleaning insectivory (the dominant energetic pathway), and most other feeding guilds also show an even larger increase (Figs. 2d and 3). However, total bird energy intake in the oil palm drops back to levels similar to those in the old-growth forest, with a collapse in multiple guilds. For mammals, there is a similar 2.4-fold (1.9–3.2) increase in total consumption when going from old-growth to logged forest, but this declines sharply in oil palm plantation. Most notable is the 5.7-fold (3.2–10.2) increase in the importance of terrestrial mammal herbivores in the logged relative to old-growth forests. All four individual old-growth forest sites show consistently lower bird and mammal energetics than the logged forests (Extended Data Fig. 5).Fig. 3: Magnitude and species diversity of energetic pathways in old-growth forest, logged forest and oil palm.The size of the circles indicates the magnitude of energy flow, and the colour indicates birds or mammals. S, number of species; E, ESWI, an index of species redundancy and, therefore, resilience (high values indicate high redundancy; see main text). For clarity, guilds with small energetic flows are not shown, but are listed in Supplementary Data 4. Images created by J. Bentley.Full size imageThe fraction of NPP flowing through the bird and mammal communities increases by a factor of 2.1 (1.5–3.0) in logged forest relative to old-growth forest. There is very little increase in NPP in logged relative to old-growth forests (Fig. 2a) because increased NPP in patches of relatively intact logged forest is offset by very low productivity in more structurally degraded areas such as former logging platforms14,15. In oil palm plantations, oil palm fruits account for a large proportion of NPP, although a large fraction of these is harvested and removed from the ecosystem17. As a proportion of NPP, 1.62% (1.35–2.13%) is directly consumed by birds and mammals in the old-growth forest; this rises to 3.36% (2.57–5.07%) in the logged forest but drops to 0.89% (0.57–1.44%) in oil palm (Fig. 2f,g and Extended Data Table 2).If all invertebrates consumed are herbivores or detritivores (that is, at a trophic level of 2.0), and trophic efficiency is 10% (ref. 10), the total amount of NPP supporting the combined bird and mammal food intake would be 9%, 16% and 5% for old-growth forest, logged forest and oil palm, respectively. However, if the mean trophic level of consumed invertebrates is 2.5 (that is, a mix of herbivores and predators), the corresponding proportions would be 27%, 51% and 17% (Fig. 2f,g). As insectivory is the dominant feeding mode for the avian community, these numbers are dominated by bird diets. For birds in the old-growth forests, 0.35% of NPP supports direct herbivory and frugivory, but around 22% of NPP (assumed invertebrate trophic level 2.5) is indirectly required to support insectivory. The equivalent numbers for birds in logged forest are 0.83% and 46%. Hence, birds account for a much larger indirect consumption of NPP. Bird diet studies in old-growth and logged forest in the region suggest that consumed invertebrates have a mean trophic level of 2.5 (ref. 18; K. Sam, personal communication), indicating that the higher-end estimates of indirect NPP consumption (that is, around 50% in logged forests) are plausible.It is interesting to compare such high fractions of NPP to direct estimates of invertebrate herbivory. Scans of tree leaf litter from these forests suggest that just 7.0% of tree canopy leaf area (1–3% of total NPP) is removed by tree leaf herbivory14,16, but such estimates do not include other pathways available to invertebrates, including herbivory of the understorey, aboveground and belowground sap-sucking, leaf-mining, fruit- and wood-feeding, and canopy, litter and ground-layer detritivory. An increase in invertebrate biomass and herbivory in logged forest compared to old-growth forest has previously been reported in fogging studies in this landscape19. Such high levels of consumption of NPP by invertebrates could have implications on ecosystem vegetation biomass production, suggesting, first, that invertebrate herbivory has a substantial influence on recovery from logging and, second, that insectivorous bird densities may exert substantial indirect controls on ecosystem recovery.The distributions of energy flows among feeding guilds are remarkably stable among habitat types (Fig. 3), indicating that the amplified energy flows in the logged forests do not distort the overall trophic structure of vertebrate communities. Overall bird diet energetics are dominated by insectivory, which accounts for a strikingly invariant 66%, 63% and 66% of bird energetic consumption in old-growth forest, logged forest and oil palm, respectively. Foliage-gleaning dominates as a mode of invertebrate consumption in all three habitat types, with frugivory being the second most energetically important feeding mode (26%, 27% and 19%, respectively). Mammal diet is more evenly distributed across feeding guilds, but frugivory (31%, 30%, 30%) and folivory (24%, 38%, 26%) dominate. Small mammal insectivores are probably under-sampled (see Methods) so the contribution of mammal insectivory may be slightly greater than that estimated here. The apparent constancy of relative magnitude of feeding pathways across the intact and disturbed ecosystems is noteworthy and not sensitive to plausible shifts in feeding behaviour between habitat types (see Supplementary Discussion). There is no evidence of a substantial shift in dominant feeding guild: the principal feeding pathways present in the old-growth forest are maintained in the logged forest.When examining change at species level in the logged forests, the largest absolute increases in bird food consumption were in arboreal insectivores and omnivores (Fig. 4a and Extended Data Fig. 2a). In particular, this change was characterized by large increases in the abundance of bulbul species (Pycnonotus spp.). No bird species showed a significant or substantial reduction in overall energy consumption. In the oil palm plantation, total food consumption by birds was less than in logged forests, but similar to that in old-growth forests. However, this was driven by very high abundance of a handful of species, notably a single arboreal omnivore (yellow-vented bulbul Pycnonotus goiavier) and three arboreal insectivores (Mixornis bornensis, Rhipidura javanica, Copsychus saularis), whereas energy flows through most other bird species were greatly reduced (Fig. 4b and Extended Data Fig. 2b).Fig. 4: Changes in energy consumption by species in logged forest and oil palm relative to old-growth forest.a,b, Changes in energy consumption by species in logged forest relative to old-growth forest (a) and in oil palm relative to old-growth forest (b). The 20 species experiencing the largest increase (red) and decrease (blue) in both habitat types are shown. Bird species are shown in a lighter tone and mammal species are shown in a darker tone. The error bars denote 95% confidence intervals, derived from 10,000 Monte Carlo simulation estimates incorporating uncertainty in body mass, population density, the daily energy expenditure equation, assimilation efficiency of the different food types and composition of the diet of each species.Full size imageFor mammals, the increase in consumption in logged forests is dominated by consumption by large terrestrial herbivores increasing by a factor of 5.7 (3.2–10.2), particularly sambar deer (Rusa unicolor) and Asian elephant (Elephas maximus; Fig. 4a and Extended Data Figs. 2b and 3), along with that by small omnivores, predominantly rodents (native spiny rats, non-native black rat; Fig. 4). A few rainforest species show a strong decline (for example, greater mouse-deer Tragulus napu and brown spiny rat Maxomys rajah). In the oil palm, most mammal species collapse (Fig. 4b) and the limited consumption is dominated by a few disturbance-tolerant habitat generalists (for example, red muntjac Muntiacus muntjak, black rat Rattus rattus, civets), albeit these species are at lower densities than observed in old-growth forest (Extended Data Fig. 2).With very few exceptions, the amplified energy flows in logged forest seem to retain the same level of resilience as in old-growth forest. The diversity and dominance of species within any pathway can be a measure of the resilience of that pathway to loss of species. We assessed energetic dominance within individual pathways by defining an energetic Shannon–Wiener index (ESWI) to examine distribution of energy flow across species; low ESWI indicates a pathway with high dependence on a few species and hence potential vulnerability (Fig. 3). The overall ESWI across guilds does not differ between the old-growth and logged forest (t2,34 = −0.363, P = 0.930), but does decline substantially from old-growth forest to oil palm (t2,34 = −3.826, P = 0.0015), and from logged forest to oil palm (t2,34 = −3.639, P = 0.0025; linear mixed-effects models, with habitat type as fixed effect and guild as random effect; for model coefficients see Supplementary Table 3).Hence, for birds, the diversity of species contributing to dominant energetic pathways is maintained in the transition from old-growth to logged forests but declines substantially in oil palm. Mammals generally show lower diversity and ESWI than birds, but six out of ten feeding guilds maintain or increase ESWI in logged forest relative to the old-growth forests but collapse in oil palm (Fig. 3). Terrestrial herbivory is the largest mammal pathway in the logged forest but is dependent on only four species and is probably the most vulnerable of the larger pathways: a few large mammals (especially sambar deer) play a dominant terrestrial herbivory role in the logged forest. In parallel, bearded pigs (Sus barbatus), the only wild suid in Borneo, form an important and functionally unique component of the terrestrial omnivory pathway. These larger animals are particularly sensitive to anthropogenic pressures such as hunting, or associated pathogenic pressures as evidenced by the recent precipitous decline of the bearded pig in Sabah due to an outbreak of Asian swine fever (after our data were collected)20.Vertebrate populations across the tropics are particularly sensitive to hunting pressure21. Our study site has little hunting, but as a sensitivity analysis we explored the energetic consequences of 50% reduction in population density of those species potentially affected by targeted and/or indiscriminate hunting (Extended Data Fig. 4). Targeted hunted species include commercially valuable birds, and gun-hunted mammals (bearded pig, ungulates, banteng and mammals with medicinal value). Indiscriminately hunted species include birds and mammals likely to be trapped with nets and snares. Hunting in the logged forests lowers both bird and mammal energy flows but still leaves them at levels higher than in faunally intact old-growth forests. Such hunting brings bird energetics levels close to (but still above) those of old-growth forests. For mammals, however, even intensively hunted logged forests seem to maintain higher energetic flows than the old-growth forests. Hence, only very heavy hunting is likely to ‘offset’ the amplified energetics in the logged forest.The amplified energetic pathways in our logged forest probably arise as a result of bottom-up trophic factors including increased resource supply, palatability and accessibility. The more open forest structure in logged forest results in more vegetation being near ground level22,23 and hence more accessible to large generalist mammal herbivores, which show the most striking increase of the mammal guilds. The increased prioritization by plants of competition for light and therefore rapid vegetation growth strategies in logged forests results in higher leaf nutrient content and reduced leaf chemical defences against herbivory24,25, along with higher fruiting and flowering rates19 and greater clumping in resource supply9. This increased resource availability and palatability probably supports high invertebrate and vertebrate herbivore densities25. The act of disturbance displaces the ecosystem from a conservative chemically defended state to a more dynamic state with amplified energy and nutrient flow, but not to an extent that causes heavy disruption in animal community composition. Top-down trophic factors might also play a role in amplifying the energy flows in intermediate trophic levels, through mechanisms such as increased protection of ground-dwelling or nesting mammals and birds from aerial predators in the dense vegetation ground layer. This might partially explain the increased abundance of rodents, but there is little evidence of trophic release at this site because of the persisting high density of mammal carnivores26. Overall, the larger number of bottom-up mechanisms and surge in invertebrate consumption suggest that increased resource supply and palatability largely explains the amplification of consumption pathways in the logged forest. An alternative possibility is that the amplified vertebrate energetics do not indicate amplified overall animal energetics but rather a large diversion of energy from unmeasured invertebrate predation pathways (for example, parasitoids); this seems unlikely but warrants further exploration.Oil palm plantations show a large decline in the proportion of NPP consumed by mammals and birds compared to logged forests12. Mammal populations collapse because they are more vulnerable and avoid humans, and there is no suite of mammal generalists that can step in27,28. Birds show a more modest decline, to levels similar to those observed in old-growth forests, as there is a broad suite of generalist species that are able to adapt to and exploit the habitat types across the disturbance gradient, and because their small size and mobility render them less sensitive to human activity29. There is a consistent decline in the oil palm in ESWI for birds and especially for mammals, indicating a substantial increase in ecosystem vulnerability in many pathways.In conclusion, our analysis demonstrates the tremendously dynamic and ecologically vibrant nature of the studied logged forests, even heavily and repeatedly logged forests such as those found across Borneo. It is likely that the patterns, mechanisms and basic ecological energetics we describe are general to most tropical forests; amplification of multiple ecosystem processes after logging has also been reported for logged forests in Kenya9, but similar detailed analyses are needed for a range of tropical forests to elucidate the importance of biogeographic, climatic or other factors. We stress that our findings do not diminish the importance of protecting structurally intact old-growth forests, but rather question the meaning of degradation by shining a new light on the ecological value of logged and other structurally ‘degraded’ forests, reinforcing their significance to the conservation agenda30. We have shown that a wide diversity of species not only persist but thrive in the logged forest environment. Moreover, such ecological vibrancy probably enhances the prospects for ecosystem structural recovery. In terms of faunal intactness, our study landscape is close to a best-case scenario because hunting pressures were low. If logged forests can be protected from heavy defaunation, our analysis demonstrates that they can be vibrant ecosystems, providing many key ecosystem functions at levels much higher than in old-growth forests. Conservation of logged forest landscapes has an essential role to play in the in the protection of global biodiversity and biosphere function. More

Putz, F. E. et al. Conserv. Lett. 5, 296–303 (2012).Article 

Google Scholar 
Blaser, J., Sarre, A., Poore, D. & Johnson, S. Status of Tropical Forest Management 2011. ITTO Tech. Ser. No. 38 (International Tropical Timber Organization, 2011); available at https://go.nature.com/3usq2an
Google Scholar 
Malhi, Y. et al. Nature https://doi.org/10.1038/s41586-022-05523-1 (2022).Article 

Google Scholar 
Zwerts, J. A. et al. Conserv. Sci. Pract. 3, e568 (2021).Article 

Google Scholar 
Wilkie, D. S., Bennett, E. L., Peres, C. A. & Cunningham, A. A. Ann. N. Y. Acad. Sci. 1223, 120–128 (2011).Article 
PubMed 

Google Scholar 
Putz, F. E., Blate, G. M., Redford, K. H., Fimbel, R. & Robinson, J. Conserv. Biol. 15, 7–20 (2001).Article 

Google Scholar  More

Doak, D. F. et al. The statistical inevitability of stability‐diversity relationships in community ecology. Am. Nat. 151, 264–276 (1998).Article 
CAS 

Google Scholar 
Schläpfer, F. & Schmid, B. Ecosystem effects of biodiversity: a classification hypotheses and exploration of empirical results. Ecol. Appl. 9, 893–912 (1999).Article 

Google Scholar 
Lehman, C. L. & Tilman, D. Biodiversity, stability, and productivity in competitive communities. Am. Nat. 156, 534–552 (2000).Article 

Google Scholar 
Allan, E. et al. More diverse plant communities have higher functioning over time due to turnover in complementary dominant species. Proc. Natl Acad. Sci. USA 108, 17034–17039 (2011).Article 
ADS 
CAS 

Google Scholar 
Isbell, F. et al. High plant diversity is needed to maintain ecosystem services. Nature 477, 199–202 (2011).Article 
ADS 
CAS 

Google Scholar 
Wagg, C. et al. Plant diversity maintains long-term ecosystem productivity under frequent drought by increasing short-term variation. Ecology 98, 2952–2961 (2017).Article 

Google Scholar 
Isbell, F. et al. Biodiversity increases the resistance of ecosystem productivity to climate extremes. Nature 526, 574–577 (2015).Article 
ADS 
CAS 

Google Scholar 
Reich, P. B. et al. Impacts of biodiversity loss escalate through time as redundancy fades. Science 336, 589–592 (2012).Article 
ADS 
CAS 

Google Scholar 
Guerrero-Ramírez, N. R. et al. Diversity-dependent temporal divergence of ecosystem functioning in experimental ecosystems. Nat. Ecol. Evol. 1, 1639–1642 (2017).Article 

Google Scholar 
Meyer, S. T. et al. Effects of biodiversity strengthen over time as ecosystem functioning declines at low and increases at high biodiversity. Ecosphere 7, e01619 (2016).Article 

Google Scholar 
Huang, Y. et al. Impacts of species richness on productivity in a large-scale subtropical forest experiment. Science 362, 80–83 (2018).Article 
ADS 
CAS 

Google Scholar 
Bongers, F. J. et al. Functional diversity effects on productivity increase with age in a forest biodiversity experiment. Nat. Ecol. Evol. 5, 1594–1603 (2021).Article 

Google Scholar 
Weisser, W. W. et al. Biodiversity effects on ecosystem functioning in a 15-year grassland experiment: Patterns, mechanisms, and open questions. Basic Appl. Ecol. 23, 1–73 (2017).Article 

Google Scholar 
Guerrero-Ramírez, N. R., Reich, P. B., Wagg, C., Ciobanu, M. & Eisenhauer, N. Diversity-dependent plant–soil feedbacks underlie long-term plant diversity effects on primary productivity. Ecosphere 10, e02704 (2019).Article 

Google Scholar 
Eisenhauer, N. The shape that matters: how important is biodiversity for ecosystem functioning. Sci. China Life Sci. 65, 651–653 (2022).Article 

Google Scholar 
Cardinale, B. J. et al. Impacts of plant diversity on biomass production increase through time because of species complementarity. Proc. Natl Acad. Sci. USA 104, 18123–18128 (2007).Article 
ADS 
CAS 

Google Scholar 
Marquard, E. et al. Plant species richness and functional composition drive overyielding in a six-year grassland experiment. Ecology 90, 3290–3302 (2009).Article 

Google Scholar 
Zuppinger-Dingley, D. et al. Selection for niche differentiation in plant communities increases biodiversity effects. Nature 515, 108–111 (2014).Article 
ADS 
CAS 

Google Scholar 
Loreau, M. & Hector, A. Partitioning selection and complementarity in biodiversity experiments. Nature 412, 72–76 (2001).Article 
ADS 
CAS 

Google Scholar 
Wang, S. et al. How complementarity and selection affect the relationship between ecosystem functioning and stability. Ecology 102, e03347 (2021).Article 

Google Scholar 
Yan, Y. et al. Mechanistic links between biodiversity effects on ecosystem functioning and stability in a multi-site grassland experiment. J. Ecol. 109, 3370–3378 (2021).Article 

Google Scholar 
Barry, K. E. et al. The future of complementarity: disentangling causes from consequences. Trends Ecol. Evol. 34, 167–180 (2019).Article 

Google Scholar 
Yachi, S. & Loreau, M. Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. Proc. Natl Acad. Sci. USA 96, 1463–1468 (1999).Article 
ADS 
CAS 

Google Scholar 
Gonzalez, A. & Loreau, M. The causes and consequences of compensatory dynamics in ecological communities. Annu. Rev. Ecol. Evol. Syst. 40, 393–414 (2009).Article 

Google Scholar 
Thibaut, L. M. & Connolly, S. R. Understanding diversity–stability relationships: towards a unified model of portfolio effects. Ecol. Lett. 16, 140–150 (2013).Article 

Google Scholar 
Craven, D. et al. Multiple facets of biodiversity drive the diversity–stability relationship. Nat. Ecol. Evol. 2, 1579–1587 (2018).Article 

Google Scholar 
Tilman, D., Reich, P. B. & Knops, J. M. H. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441, 629–632 (2006).Article 
ADS 
CAS 

Google Scholar 
Loreau, M. Biodiversity and Ecosystem Functioning (Princeton Univ. Press,2010).Loreau, M. & de Mazancourt, C. Biodiversity and ecosystem stability: a synthesis of underlying mechanisms. Ecol. Lett. 16, 106–115 (2013).Article 

Google Scholar 
Isbell, F. et al. Quantifying effects of biodiversity on ecosystem functioning across times and places. Ecol. Lett. 21, 763–778 (2018).Article 

Google Scholar 
Maron, J. L., Marler, M., Klironomos, J. N. & Cleveland, C. C. Soil fungal pathogens and the relationship between plant diversity and productivity. Ecol. Lett. 14, 36–41 (2011).Article 

Google Scholar 
Schnitzer, S. A. et al. Soil microbes drive the classic plant diversity–productivity pattern. Ecology 92, 296–303 (2011).Article 

Google Scholar 
Marquard, E. et al. Changes in the abundance of grassland species in monocultures versus mixtures and their relation to biodiversity effects. PLoS ONE 8, e75599 (2013).Article 
ADS 
CAS 

Google Scholar 
Roscher, C. et al. Identifying population- and community-level mechanisms of diversity–stability relationships in experimental grasslands. J. Ecol. 99, 1460–1469 (2011).Article 

Google Scholar 
Civitello, D. J. et al. Biodiversity inhibits parasites: broad evidence for the dilution effect. Proc. Natl Acad. Sci. USA 112, 8667–8671 (2015).Article 
ADS 
CAS 

Google Scholar 
Kulmatiski, A., Beard, K. H. & Heavilin, J. Plant–soil feedbacks provide an additional explanation for diversity–productivity relationships. Proc. R. Soc. B 279, 3020–3026 (2012).Article 

Google Scholar 
van Moorsel, S. J. et al. Co-occurrence history increases ecosystem stability and resilience in experimental plant communities. Ecology 102, e03205 (2021).Article 

Google Scholar 
Schöb, C., Brooker, R. W. & Zuppinger-Dingley, D. Evolution of facilitation requires diverse communities. Nat. Ecol. Evol. 2, 1381–1385 (2018).Article 

Google Scholar 
Temperton, V. M., Mwangi, P. N., Scherer-Lorenzen, M., Schmid, B. & Buchmann, N. Positive interactions between nitrogen-fixing legumes and four different neighbouring species in a biodiversity experiment. Oecologia 151, 190–205 (2007).Article 
ADS 

Google Scholar 
Furey, G. N. & Tilman, D. Plant biodiversity and the regeneration of soil fertility. Proc. Natl Acad. Sci. USA 118, e2111321118 (2021).Article 
CAS 

Google Scholar 
Gubsch, M. et al. Foliar and soil δ15N values reveal increased nitrogen partitioning among species in diverse grassland communities. Plant Cell Environ. 34, 895–908 (2011).Article 
CAS 

Google Scholar 
Roscher, C., Schmid, B., Buchmann, N., Weigelt, A. & Schulze, E.-D. Legume species differ in the responses of their functional traits to plant diversity. Oecologia 165, 437–452 (2011).Article 
ADS 

Google Scholar 
Eisenhauer, N. et al. Plant diversity effects on soil microorganisms support the singular hypothesis. Ecology 91, 485–496 (2010).Article 
CAS 

Google Scholar 
Fornara, D. A. & Tilman, D. Plant functional composition influences rates of soil carbon and nitrogen accumulation. J. Ecol. 96, 314–322 (2008).Article 
CAS 

Google Scholar 
Lange, M. et al. Plant diversity increases soil microbial activity and soil carbon storage. Nat. Commun. 6, 6707 (2015).Article 
ADS 
CAS 

Google Scholar 
Xu, S. et al. Species richness promotes ecosystem carbon storage: evidence from biodiversity-ecosystem functioning experiments. Proc. R. Soc. B 287, 20202063 (2020).Article 
CAS 

Google Scholar 
Cong, W.-F. et al. Plant species richness promotes soil carbon and nitrogen stocks in grasslands without legumes. J. Ecol. 102, 1163–1170 (2014).Article 
CAS 

Google Scholar 
Leimer, S. et al. Mechanisms behind plant diversity effects on inorganic and organic N leaching from temperate grassland. Biogeochemistry 131, 339–353 (2016).Article 
CAS 

Google Scholar 
Xu, Q. et al. Consistently positive effect of species diversity on ecosystem, but not population, temporal stability. Ecol. Lett. 24, 2256–2266 (2021).Article 

Google Scholar 
Hector, A. et al. General stabilizing effects of plant diversity on grassland productivity through population asynchrony and overyielding. Ecology 91, 2213–2220 (2010).Article 
CAS 

Google Scholar 
Turnbull, L. A., Levine, J. M., Loreau, M. & Hector, A. Coexistence, niches and biodiversity effects on ecosystem functioning. Ecol. Lett. 16, 116–127 (2013).Article 

Google Scholar 
Wright, A. J. et al. Flooding disturbances increase resource availability and productivity but reduce stability in diverse plant communities. Nat. Commun. 6, 6092 (2015).Article 
ADS 
CAS 

Google Scholar 
Fischer, F. M. et al. Plant species richness and functional traits affect community stability after a flood event. Philos. Trans. R. Soc. B 371, 20150276 (2016).Article 

Google Scholar 
Roscher, C. et al. A functional trait-based approach to understand community assembly and diversity–productivity relationships over 7 years in experimental grasslands. Perspect. Plant Ecol. Evol. Syst. 15, 139–149 (2013).Article 

Google Scholar 
Eisenhauer, N. et al. Biotic interactions, community assembly, and eco-evolutionary dynamics as drivers of long-term biodiversity–ecosystem functioning relationships. Res. Ideas Outcomes 5, e47042 (2019).Article 

Google Scholar 
van Moorsel, S. J., Schmid, M. W., Hahl, T., Zuppinger-Dingley, D. & Schmid, B. Selection in response to community diversity alters plant performance and functional traits. Perspect. Plant Ecol. Evol. Syst. 33, 51–61 (2018).Article 

Google Scholar 
van Moorsel, S. J. et al. Community evolution increases plant productivity at low diversity. Ecol. Lett. 21, 128–137 (2018).Article 

Google Scholar 
Roeder, A. et al. Plant diversity effects on plant longevity and their relationships to population stability in experimental grasslands. J. Ecol. 109, 2566–2579 (2021).Article 

Google Scholar 
Cadotte, M. W., Dinnage, R. & Tilman, D. Phylogenetic diversity promotes ecosystem stability. Ecology 93, S223–S233 (2012).Article 

Google Scholar 
Pu, Z., Daya, P., Tan, J. & Jiang, L. Phylogenetic diversity stabilizes community biomass. J. Plant Ecol. 7, 176–187 (2014).Article 

Google Scholar 
Carrara, F., Giometto, A., Seymour, M., Rinaldo, A. & Altermatt, F. Experimental evidence for strong stabilizing forces at high functional diversity of aquatic microbial communities. Ecology 96, 1340–1350 (2015).Article 

Google Scholar 
Hooper, D. U. et al. Effects of biodiversity on ecosystem functioning: a conceensus of current knowledge. Ecol. Monogr. 75, 3–35 (2005).Article 

Google Scholar 
Ruijven, J. V. & Berendse, F. Contrasting effects of diversity on the temporal stability of plant populations. Oikos 116, 1323–1330 (2007).Article 

Google Scholar 
Proulx, R. et al. Diversity promotes temporal stability across levels of ecosystem organization in experimental grasslands. PLoS ONE 5, e13382 (2010).Article 
ADS 

Google Scholar 
Hoaglin, D. C., Iglewicz, B. & Tukey, J. W. Performance of some resistant rules for outlier labeling. JASA 81, 991–999 (1986).Article 
MathSciNet 

Google Scholar 
Loreau, M. & de Mazancourt, C. Species synchrony and its drivers: neutral and nonneutral community dynamics in fluctuating environments. Am. Nat. 172, E48–E66 (2008).Article 

Google Scholar 
Gross, K. et al. Species richness and the temporal stability of biomass production: a new analysis of recent biodiversity experiments. Am. Nat. 183, 1–12 (2014).Article 

Google Scholar 
Schmid, B., Baruffol, M., Wang, Z. & Niklaus, P. A. A guide to analyzing biodiversity experiments. J. Plant Ecol. 10, 91–110 (2017).Article 

Google Scholar 
Rosseel, Y. lavaan: An R package for structural equation modeling. J. Stat. Softw. 48, 1–36 (2012).Article 

Google Scholar  More

Gao, A. L. & Zhou, K. Y. Growth and reproduction of three populations of finless porpoise, Neophocaena phocaenoides, in Chinese waters. Aquat Mamm 19, 3–12 (1993).
Google Scholar 
Jefferson, T. A. Preliminary analysis of geographic variation in cranial morphometrics of the finless porpoise (Neophocaena phocaenoides). Raffles Bull Zool 10, 3–14 (2002).
Google Scholar 
Pilleri, G. & Gihr, M. Contribution to the knowledge of the cetaceans of Pakistan with particular reference to the genera Neomeris, Sousa, Delphinus and Tursiops and description of a new Chinese porpoise (Neomeris asiaeorientalis). Investig Cetacea 4, 107–162 (1972).
Google Scholar 
Pilleri, G. & Gihr, M. On the taxonomy and ecology of the finless black porpoise, Neophocaena (Cetacea, Delphinidae). Mammalia 39, 657–673 (1975).Article 

Google Scholar 
Wang, P. L. The morphological characters and the problem of subspecies identifications of the finless porpoise. Fish Sci 11, 4–8 (1992).
Google Scholar 
Wang, P. L. On the taxonomy of the finless porpoise in China. Fish Sci 6, 10–14 (1992).
Google Scholar 
Gao, A. L. & Zhou, K. Y. Geographical variation of external measurements and three subspecies of Neophocaena phocaenoides in Chinese waters. Acta Theriol Sin 15, 81–92 (1995).
Google Scholar 
Wang, J. Y., Frasier, T. R., Yang, S. C. & White, B. N. Detecting recent speciation events: the case of the finless porpoise (genus Neophocaena). Heredity 101, 145–155 (2008).Article 

Google Scholar 
Jefferson, T. A. & Wang, J. Y. Revision of the taxonomy of finless porpoises (genus Neophocaena): the existence of two species. J Mar Anim Ecol 4, 3–16 (2011).
Google Scholar 
Zhou, X. M. et al. Population genomics of finless porpoises reveal an incipient cetacean species adapted to freshwater. Nat Commun 9, 1276 (2018).Article 
ADS 

Google Scholar 
Wang, D., Turvey, S.T., Zhao, X. & Mei, Z. Neophocaena asiaeorientalis ssp. asiaeorientalis. The IUCN Red List of Threatened Species https://www.iucnredlist.org/species/43205774/45893487 (2013).Wang, J. Y. & Reeves, R. Neophocaena Asiaeorientalis. The IUCN Red List of Threatened Species https://www.iucnredlist.org/species/41754/50381766 (2017).Kasuya, T. Japanese whaling and other cetacean fisheries. Environ Sci Pollut Res Int 14, 39–48 (2007).Article 

Google Scholar 
Yoshida, H., Shirakihara, K., Kishino, H. & Shirakihara, M. A population size estimate of the finless porpoise, Neophocaena phocaenoides, from aerial sighting surveys in Ariake Sound and Tachibana Bay, Japan. Popul Ecol 39, 239–247 (1997).Article 

Google Scholar 
Amano, M., Nakahara, F., Hayano, A. & Shirakihara, K. Abundance estimate of finless porpoises off the Pacific coast of eastern Japan based on aerial surveys. Mamm Study 28, 103–110 (2003).Article 

Google Scholar 
Shirakihara, K., Shirakihara, M. & Yamamoto, Y. Distribution and abundance of finless porpoise in the Inland Sea of Japan. Mar Biol 150, 1025–1032 (2007).Article 

Google Scholar 
Zuo, T., Sun, J. Q., Shi, Y. Q. & Wang, J. Primary survey of finless porpoise population in the Bohai Sea. Acta Theriol Sin 38, 551–561 (2018).
Google Scholar 
Ruan, R., Guo, A. H., Hao, Y. J., Zheng, J. S. & Wang, D. De novo assembly and characterization of narrow-ridged finless porpoise renal transcriptome and identification of candidate genes involved in osmoregulation. Int J Mol Sci 16, 2220–2238 (2015).Article 

Google Scholar 
Li, S. H. et al. Echolocation click sounds from wild inshore finless porpoise (Neophocaena phocaenoides sunameri) with comparisons to the sonar of riverine N. p. asiaeorientalis. J Acoust Soc Am 121, 3938–3946 (2007).Article 
ADS 

Google Scholar 
Dong, J. H., Wang, G. J. & Xiao, Z. Z. Migration and population difference of the finless porpoise in China. Mar Sci 5, 42–45 (1993).
Google Scholar 
Lu, Z. C. et al. Analysis of the diet of finless porpoise (Neophocaena asiaeorientalis sunameri) based on prey morphological characters and DNA barcoding. Conserv Genet Resour 8, 523–531 (2016).Article 

Google Scholar 
Chen, B. et al. Finless porpoises (Neophocaena asiaeorientalis) in the East China Sea: insights into feeding habits using morphological, molecular, and stable isotopic techniques. Can J Fish Aquat Sci 74, 1628–1645 (2017).Article 

Google Scholar 
Nurk, S. et al. The complete sequence of a human genome. Science 376, 44–53 (2022).Article 
ADS 

Google Scholar 
Chen, Y. X. et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 7, 1–6 (2018).Article 
ADS 

Google Scholar 
Chikhi, R. & Medvedev, P. Informed and automated k-mer size selection for genome assembly. Bioinformatics 30, 31–37 (2014).Article 

Google Scholar 
Chin, C. S. et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 10, 563–569 (2013).Article 

Google Scholar 
Cheng, H. Y., Concepcion, G. T., Feng, X. W., Zhang, H. W. & Li, H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods 18, 170–175 (2021).Article 

Google Scholar 
Roach, M. J., Schmidt, S. A. & Borneman, A. R. Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinformatics 19, 1–10 (2018).Article 

Google Scholar 
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).Article 

Google Scholar 
Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst 3, 95–98 (2016).Article 

Google Scholar 
Dudchenko, O. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92–95 (2017).Article 
ADS 

Google Scholar 
Xiong, Y., Brandley, M. C., Xu, S. X., Zhou, K. Y. & Yang, G. Seven new dolphin mitochondrial genomes and a time-calibrated phylogeny of whales. BMC Evol Biol 9, 1–13 (2009).Article 

Google Scholar 
Alonge, M. et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol 20, 1–17 (2019).Article 

Google Scholar 
Mayer, A., Lahr, G., Swaab, D. F., Pilgrim, C. & Reisert, I. The Y-chromosomal genes SRY and ZFY are transcribed in adult human brain. Neurogenetics 1, 281–288 (1998).Article 

Google Scholar 
Sinclair, A. H. et al. A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif. Nature 346, 240–244 (1990).Article 
ADS 

Google Scholar 
Koopman, P., Gubbay, J., Vivian, N., Goodfellow, P. & Lovell-Badge, R. Male development of chromosomally female mice transgenic for Sry. Nature 351, 117–121 (1991).Article 
ADS 

Google Scholar 
Salo, P. et al. Molecular mapping of the putative gonadoblastoma locus on the Y chromosome. Genes Chromosomes Cancer 14, 210–214 (1995).Article 

Google Scholar 
Tsuchiya, K., Reijo, R., Page, D. C. & Disteche, C. M. Gonadoblastoma: molecular definition of the susceptibility region on the Y chromosome. Am J Hum Genet 57, 1400–1407 (1995).
Google Scholar 
Gegenschatz-Schmid, K., Verkauskas, G., Stadler, M. B. & Hadziselimovic, F. Genes located in Y-chromosomal regions important for male fertility show altered transcript levels in cryptorchidism and respond to curative hormone treatment. Basic Clin Androl 29, 1–8 (2019).Article 

Google Scholar 
Chen, N. Using Repeat Masker to identify repetitive elements in genomic sequences. Curr protoc Bioinf 5, 4–10 (2004).Article 

Google Scholar 
Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res 35, W265–W268 (2007).Article 

Google Scholar 
Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, i351–i358 (2005).Article 

Google Scholar 
Bao, W. D., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob DNA 6, 1–6 (2015).Article 

Google Scholar 
Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res 27, 573–580 (1999).Article 

Google Scholar 
Liu, W. et al. Blood Transcriptome Analysis Reveals Gene Expression Differences between Yangtze Finless Porpoises from Two Habitats: Natural and Ex Situ Protected Waters. Fishes 7, 96 (2022).Article 

Google Scholar 
Yin, D. H. et al. Integrated analysis of blood mRNAs and microRNAs reveals immune changes with age in the Yangtze finless porpoise (Neophocaena asiaeorientalis). Comp Biochem Physiol B Biochem Mol Biol 256, 110635 (2021).Article 

Google Scholar 
Kim, D., Paggi, J. M., Park, C., Bennett, C. & Salzberg, S. L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37, 907–915 (2019).Article 

Google Scholar 
Kovaka, S. et al. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol 20, 1–13 (2019).Article 

Google Scholar 
Stanke, M., Diekhans, M., Baertsch, R. & Haussler, D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24, 637–644 (2008).Article 

Google Scholar 
Keane, M. et al. Insights into the evolution of longevity from the bowhead whale genome. Cell Rep 10, 112–122 (2015).Article 

Google Scholar 
Yim, H. S. et al. Minke whale genome and aquatic adaptation in cetaceans. Nat Genet 46, 88–92 (2014).Article 

Google Scholar 
Jones, S. J. et al. The genome of the beluga whale (Delphinapterus leucas). Genes 8, 378 (2017).Article 
ADS 

Google Scholar 
Zhou, X. M. et al. Baiji genomes reveal low genetic variability and new insights into secondary aquatic adaptations. Nat Commun 4, 1–6 (2013).Article 
ADS 

Google Scholar 
Foote, A. D. et al. Convergent evolution of the genomes of marine mammals. Nat Genet 47, 272–275 (2015).Article 

Google Scholar 
Keilwagen, J., Hartung, F. & Grau, J. GeMoMa: homology-based gene prediction utilizing intron position conservation and RNA-seq data. Methods Mol Biol 1962, 161–177 (2019).Article 

Google Scholar 
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44, D457–D462 (2016).Article 

Google Scholar 
Bairoch, A. & Apweiler, R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res 28, 45–48 (2000).Article 

Google Scholar 
Korf, I. Gene finding in novel genomes. BMC bioinformatics 5, 1–9 (2004).Article 

Google Scholar 
Finn, R. D. et al. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res 45, D190–D199 (2017).Article 

Google Scholar 
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J Mol Biol 215, 403–410 (1990).Article 

Google Scholar 
Mulder, N. J. & Apweiler, R. InterPro and InterProScan: tools for protein sequence classification and comparison. Methods Mol Biol 396, 59–70 (2007).Article 

Google Scholar 
Ashburner, M. et al. Gene ontology: tool for the unification of biology. Nat Genet 25, 25–29 (2000).Article 

Google Scholar 
NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR21047154 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20760935 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20760936 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20997931 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20997932 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20997933 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20997934 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRR20997935 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRP389529 (2022).Yin, D. H. et al. Neophocaena asiaeorientalis sunameri isolate NAS202207, whole genome shotgun sequencing project. GenBank https://identifiers.org/insdc.gca:GCA_026225855.1 (2022).Yin, D. H. et al. Gapless genome assembly of East Asian finless porpoise, Neophocaena asiaeorientalis sunameri. figshare https://doi.org/10.6084/m9.figshare.20381274.v2 (2022).Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212 (2015).Article 

Google Scholar 
Marçais, G. et al. MUMmer4: A fast and versatile genome alignment system. PLoS Comput Biol 14, e1005944 (2018).Article 

Google Scholar  More

Tilman, D., May, R. M., Lehman, C. L. & Nowak, M. A. Habitat destruction and the extinction debt. Nature 371, 65–66 (1994).Article 

Google Scholar 
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).Article 
CAS 
PubMed 

Google Scholar 
Urban, M. C. Accelerating extinction risk from climate change. Science 348, 571–573 (2015).Article 
CAS 
PubMed 

Google Scholar 
Fonseca, C. R. et al. Conservation biology: four decades of problem- and solution-based research. Perspect. Ecol. Conserv. 19, 121–130 (2021).
Google Scholar 
Smits, P. & Finnegan, S. How predictable is extinction? Forecasting species survival at million-year timescales. Philos. Trans. R. Soc. B Biol. Sci. 374, 20190392 (2019).Article 

Google Scholar 
Hanski, I. & Ovaskainen, O. Extinction debt at extinction threshold. Conserv. Biol. 16, 666–673 (2002).Article 

Google Scholar 
Kuussaari, M. et al. Extinction debt: a challenge for biodiversity conservation. Trends Ecol. Evol. 24, 564–571 (2009).Article 
PubMed 

Google Scholar 
Ridding, L. E. et al. Inconsistent detection of extinction debts using different methods. Ecography 44, 33–43 (2021).Article 

Google Scholar 
Berglund, H. & Jonsson, B. G. Verifying an extinction debt among lichens and fungi in northern Swedish boreal forests. Conserv. Biol. 19, 338–348 (2005).Article 

Google Scholar 
Jones, I. L., Bunnefeld, N., Jump, A. S., Peres, C. A. & Dent, D. H. Extinction debt on reservoir land-bridge islands. Biol. Conserv. 199, 75–83 (2016).Article 

Google Scholar 
Triantis, K. et al. Extinction debt on oceanic islands. Ecography 33, 285–294 (2010).
Google Scholar 
Wearn, O. R., Reuman, D. C. & Ewers, R. M. Extinction debt and windows of conservation opportunity in the Brazilian Amazon. Science 337, 228–232 (2012).Article 
CAS 
PubMed 

Google Scholar 
Pan, Y. et al. Spatial and temporal scales of landscape structure affect the biodiversity-landscape relationship across ecologically distinct species groups. Landsc. Ecol. 37, 2311–2325 (2022).Article 

Google Scholar 
Soga, M. & Koike, S. Mapping the potential extinction debt of butterflies in a modern city: Implications for conservation priorities in urban landscapes. Anim. Conserv. 16, 1–11 (2013).Article 

Google Scholar 
Knapp, S., Winter, M. & Klotz, S. Increasing species richness but decreasing phylogenetic richness and divergence over a 320-year period of urbanization. J. Appl. Ecol. 54, 1152–1160 (2017).Article 

Google Scholar 
McGill, B. J., Dornelas, M., Gotelli, N. J. & Magurran, A. E. Fifteen forms of biodiversity trend in the anthropocene. Trends Ecol. Evol. 30, 104–113 (2015).Article 
PubMed 

Google Scholar 
Chen, Y. & Peng, S. Evidence and mapping of extinction debts for global forest-dwelling reptiles, amphibians and mammals. Sci. Rep. 7, 1–10 (2017).
Google Scholar 
Krauss, J. et al. Habitat fragmentation causes immediate and time-delayed biodiversity loss at different trophic levels. Ecol. Lett. 13, 597–605 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Cowlishaw, G. Predicting the pattern of decline of African primate diversity: An extinction debt from historical deforestation. Conserv. Biol. 13, 1183–1193 (1999).Article 

Google Scholar 
Figueiredo, L., Krauss, J., Steffan-Dewenter, I. & Sarmento Cabral, J. Understanding extinction debts: spatio–temporal scales, mechanisms and a roadmap for future research. Ecography 42, 1973–1990 (2019).Article 

Google Scholar 
Aerts, R. & Honnay, O. Forest restoration, biodiversity and ecosystem functioning. BMC Ecol. 11, 1–21 (2011).Article 

Google Scholar 
Haddad, N. M. et al. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci. Adv. 1, e1500052 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Maxwell, S. L. et al. Area-based conservation in the twenty-first century. Nature 586, 217–227 (2020).Article 
CAS 
PubMed 

Google Scholar 
IUCN. The IUCN Red List of Threatened Species, Version 2019-1. https://www.iucnredlist.org. Downloaded on 23 February 2022. (2019).Brown, J. L. et al. Spatial biodiversity patterns of Madagascar’s amphibians and reptiles. PLoS ONE 11, e0144076 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Powney, G. D., Grenyer, R., Orme, C. D. L., Owens, I. P. F. & Meiri, S. Hot, dry and different: Australian lizard richness is unlike that of mammals, amphibians and birds. Glob. Ecol. Biogeogr. 19, 386–396 (2010).Article 

Google Scholar 
Pianka, E. R. Desert lizard diversity: additional comments and some data. Am. Nat. 134, 344–364 (1989).Article 

Google Scholar 
Chen, Y. H. Combining the species-area-habitat relationship and environmental cluster analysis to set conservation priorities: A study in the Zhoushan Archipelago, China. Conserv. Biol. 23, 537–545 (2009).Article 
PubMed 

Google Scholar 
Ricklefs, R. E. & Lovette, I. J. The roles of island area per se and habitat diversity in the species-area relationships of four Lesser Antillean faunal groups. J. Anim. Ecol. 68, 1142–1160 (1999).Article 

Google Scholar 
Souza, F. L., Martins, F. I. & Raizer, J. Habitat heterogeneity and anuran community of an agroecosystem in the Pantanal of Brazil. Phyllomedusa 13, 41–50 (2014).Article 

Google Scholar 
Kelt, D. A. & Van Vuren, D. H. The ecology and macroecology of mammalian home range area. Am. Nat. 157, 637–645 (2001).Article 
CAS 
PubMed 

Google Scholar 
McNab, B. K. Bioenergetics and the determination of home range size. Am. Nat. 97, 133–140 (1963).Article 

Google Scholar 
Powell, R. A. & Mitchell, M. S. What is a home range? J. Mammal. 93, 948–958 (2012).Article 

Google Scholar 
Hoffmann, S., Irl, S. D. H. & Beierkuhnlein, C. Predicted climate shifts within terrestrial protected areas worldwide. Nat. Commun. 10, 1–10 (2019).Article 

Google Scholar 
Giam, X. et al. Reservoirs of richness: least disturbed tropical forests are centres of undescribed species diversity. Proc. R. Soc. B 279, 67–76 (2012).Article 
PubMed 

Google Scholar 
Pillay, R. et al. Tropical forests are home to over half of the world’s vertebrate species. Front. Ecol. Environ. 20, 10–15 (2022).Article 
PubMed 

Google Scholar 
Li, H. et al. Large numbers of vertebrates began rapid population decline in the late 19th century. Proc. Natl Acad. Sci. USA 113, 14079–14084 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pringle, R. M. Upgrading protected areas to conserve wild biodiversity. Nature 546, 91–99 (2017).Article 
CAS 
PubMed 

Google Scholar 
Forzieri, G., Dakos, V., McDowell, N. G., Ramdane, A. & Cescatti, A. Emerging signals of declining forest resilience under climate change. Nature 608, 534–539 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Diamond, J. M. Biogeographic kinetics: estimation of relaxation times for Avifaunas of southwest Pacific islands. Proc. Natl Acad. Sci. USA 69, 3199–3203 (1972).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Jackson, S. T. & Sax, D. F. Balancing biodiversity in a changing environment: extinction debt, immigration credit and species turnover. Trends Ecol. Evol. 25, 153–160 (2010).Article 
PubMed 

Google Scholar 
Foley, J. A. et al. Amazonia revealed: forest degradation and loss of ecosystem goods and services in the Amazon Basin. Front. Ecol. Environ. 5, 25–32 (2007).Article 

Google Scholar 
Asamoah, E. F., Beaumont, L. J. & Maina, J. M. Climate and land-use changes reduce the benefits of terrestrial protected areas. Nat. Clim. Chang. 11, 1105–1110 (2021).Article 

Google Scholar 
Hurtt, G. C. et al. Harmonization of land-use scenarios for the period 1500–2100: 600 years of global gridded annual land-use transitions, wood harvest, and resulting secondary lands. Clim. Change 109, 117–161 (2011).Article 

Google Scholar 
Peng, S. et al. Sensitivity of land use change emission estimates to historical land use and land cover mapping. Glob. Biogeochem. Cycles 31, 626–643 (2017).Article 
CAS 

Google Scholar 
Jain, A. K., Meiyappan, P., Song, Y. & House, J. I. CO2 emissions from land-use change affected more by nitrogen cycle, than by the choice of land-cover data. Glob. Chang. Biol. 19, 2893–2906 (2013).Article 
PubMed 

Google Scholar 
Poulter, B. et al. Plant functional type classification for earth system models: results from the European Space Agency’s Land Cover Climate Change Initiative. Geosci. Model Dev. 8, 2315–2328 (2015).Article 

Google Scholar 
Pongratz, J., Reick, C., Raddatz, T. & Claussen, M. A reconstruction of global agricultural areas and land cover for the last millennium. Global Biogeochem. Cycles 22, (2008).Dietz, F. C. The industrial revolution. In the Hands of a Child (1970).Gütschow, J., Jeffery, L. & Gieseke, R. The PRIMAP-hist national historical emissions time series (1850-2016). V. 2.0. GFZ Data Services (2019).Dinerstein, E. et al. An ecoregion-based approach to protecting half the terrestrial realm. Bioscience 67, 534–545 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Protected Planet: The World Database on Protected Areas (UNEP-WCMC and IUCN, accessed 9 January 2022); www.protectedplanet.net.Butchart, S. H. M. et al. Shortfalls and solutions for meeting national and global conservation area targets. Conserv. Lett. 8, 329–337 (2015).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing Version 4.0.2 (2020). More

Whale entanglements in fishing gear threaten whale populations, seafood production and long-term sustainability of commercial fisheries. While multiple mitigation strategies to reduce entanglements exist, there has been minimal consideration of the economic impact of these strategies. Here, we estimated retrospective losses to ex-vessel revenues for one of California’s most lucrative fisheries. Overall, we found fishery closures decreased ex-vessel revenue, with results showing some uncertainty due to large model prediction error. Regional differences in losses revealed interesting trends in the capacity for the fishery to recoup costs. For example, in the NMA, relatively small losses at the fishery level were predicted ($0.3 million in total) for the 2019 season despite an early closure to the season due to whale entanglement risk.NMA fishers collectively were able to meet predicted revenue for the season despite a shortening of the fishing 2019 season. In the 2020 season however, the NMA did not experience disturbances due to whale entanglements but larger ex-vessel losses (of $3.9 million) were predicted. This suggests that other disturbances such as a delay to the season due to crab meat quality, lost fishing opportunity related to the COVID-19 pandemic, or other unknown factors, had an influence on ex-vessel revenue during the 2020 season. While most of the 2020 season landings in the NMA occurred before COVID-19 arrived in the US, there is evidence that prices in latter part of the season may have been depressed due to loss of export markets for live crab47.In the CMA however, despite landing the majority of crab available during the 2019 season (see Fig. 2c), losses of $9.4 million were experienced across the fishery. While total fishery catch was not greatly reduced, closure to the fishery in the spring may be responsible for revenue losses through other mechanisms (e.g. price). In the 2020 season, whale entanglement risk substantially shortened the fishing season in the CMA, through a delay at the beginning of the season and an early closure in the spring. Estimated losses were largest ($14.4 million) during this season. It is likely that the COVID-19 pandemic was also responsible for some of this estimated loss in the CMA in the 2020 season47. Our model did not control for impacts of the COVID-19 pandemic. However, price trends suggest that that price of Dungeness Crab in California was not affected until mid-March 2020, at which point the fishery had caught 92% of the seasons catch (see Supplementary File S2). Prices then returned to normal levels in mid-May. If we apply extrapolated prices between mid-March and mid-May by replacing observed prices with linearly increasing prices by week, revenues would have been $753,754 higher in total across the fishery. This rough estimate suggests we can attribute 4.1% of overall estimated revenue losses during the 2020 season to COVID-19 impacts, with the caveat that we do not know what prices would have been in the absence of the pandemic. A counterfactual approach has been used to disentangle multiple stressors to infer causal impacts of management interventions elsewhere48, however as these closures, and the COVID pandemic, potentially impacted all fishers in the California Dungeness crab fishery, there are no control groups available for comparison and therefore this approach would not be appropriate.Closures and other disturbances appear to have been less impactful in the NMA and high price for Dungeness crab may have contributed to the ability of vessels operating in the NMA to withstand disturbances (Supplementary Fig. S2). Prices were particularly high during the summer portion of the season in 2020 during which time the CMA was closed to Dungeness crab fishing (Supplementary Fig. S2). The NMA did not experience closures due to whale entanglement during 2020 and was predicted to have lower than average pre-season abundance (lower catch potential) during 2020 (see Fig. 2.b), while the CMA was predicted to have high catch potential for 2020 (Fig. 2.c), therefore differences in management measures implemented, and seasons’ catch potential, also contributed to differences in losses estimated.The CMA also experienced high prices, including decadal high prices for crab during the November–December of the 2019 fishing season (Supplementary Fig. S2). However, losses observed overall across the two seasons suggest the fishery, unlike the NMA, did not get much overall benefit from the high price in 2019 or the high pre-season abundance of crab (i.e. catch potential) estimated for the 2020 season in the CMA. A number of factors may have contributed to a poor season in the CMA including catchability or biology of Dungeness crab as well as external factors such as the COVID-19 pandemic behavioral choice factors, for example deciding not to fish45. Temporally shifting or reducing the opportunity for participation through closed periods due to whale entanglement risk may have exacerbated other impacts on revenues in the CMA which were not as impactful on revenues in the NMA.The high variability in estimated economic impacts per vessel reported here demonstrates that closures did not affect all vessels equally, similarly to impacts observed following a climate related harmful algal bloom in the 2016 season which were variable by vessel size and between communities45. The estimated losses we present at the fishery level in the NMA and CMA may therefore be underestimated, or overestimated, for particular groups of vessels within those management areas. This reflects the diverse nature of the Dungeness Crab fishery in behaviour and fishing strategy and highlights the importance of capturing impacts at finer scales than the fishery level alone.Limitations to the estimation of closure impactsA limitation of the hurdle model is that there are other latent factors influencing fishery participation and revenues that our model does not incorporate, particularly those determining fisher behavior such as fuel price, shipyard backlogs and market demand. A behavioral choice model, for example one that incorporates location or fishing alternative choice given a closure50,51,52 would be a potential method to better understand how spatial management strategies affect fisher behavior and is recommended as a future analysis to assess trade-offs involving socio-economic risk. Our results, reporting losses from Dungeness crab fishing revenue only, also do not account for the ability of some fishers to mitigate revenue losses by participating in other fisheries. Dungeness crab fishing is highly connected within west coast fishery participation networks44,45. Thus, it is important to note that our results for the 2019 and 2020 seasons present only losses from Dungeness crab fishing and may overestimate total annual revenue losses by some vessels that are able to mitigate impacts with participation in other fisheries.The model, predicting out-of-sample, over-estimated revenues in recent years suggesting that our predictions of revenues may also be over-predicted. An improved estimation at the vessel level, given some over-estimation of vessels that did not fish, could be investigated through a selection model approach rather than a two-part model approach54. However, two-part models are most appropriate for estimation of conditional (actual) outcomes as was intended here rather than unconditional (potential) outcomes and they do not require separate drivers for the selection and estimation model, which we did not have available54. When the impacts of policy interventions are difficult to disentangle from other impacts, approaches such as a counterfactual synthetic control48 approach could be used to separate the impacts of the policy alone. In this context, however, it is useful to report the cumulative impact of disturbances given that these disturbances (e.g., delays due to crab quality, harmful algal blooms) happen frequently and therefore the closures will rarely happen in isolation.Whilst there are limitations to our approach, revenue predictions presented here offer more insight compared to predicting revenues based only on a 5-year average of total fishery revenues (Supplementary Table S3) as is commonly conducted to calculate disaster assistance requirement, as our analysis includes an estimation of crab abundance as well as historical vessel level data in its estimation. Accounting for the influence of crab abundance is critical in this fishery given abundance is highly variable and the majority of fishable biomass is taken each year. Estimation of revenue at the individual vessel level allows for consideration of fishery heterogeneity (e.g., by vessel size). Revenues calculated on a 5-year average would suggest total California Commercial Dungeness crab fishery revenues would have been $10.62 million higher than observed in 2019 and $12.73 million higher than observed in 2020 (Supplementary Table S3). Thus, revenues estimated on the 5-year average suggest that losses would have been $0.97 million higher than our model prediction across the fishery for 2019 and $5.56 million lower than our model prediction for 2020. Our predictions suggest that delays and closures due to whale entanglement mitigation and other disturbances in to the 2019 and 2020 seasons were similar to the impact of closures due to the HAB in the 2016 season, which were estimated at $13.6 million in losses from Dungeness Crab revenues across the fishery38.Economic cost of mitigationMany strategies that prevent fishery interactions with marine mammals exist, including gear reductions or modifications, depth limitations and dynamic or seasonal time-area closures13,14,22,23,24,25,26,55. Whilst the fishery does implement pro-active gear modification measures set out in the best practices guide34, only two management intervention options were enacted in the 2019 and 2020 seasons to mitigate against entanglements of marine life with Dungeness crab gear; delays to the start of the crab season in the winter and early closures in spring due to overlap with whale distribution in fishing grounds. These delays and closures can have differential impacts on the fishery as the fishing season is not heterogeneously prosperous. An example is that closures during the holiday season (Nov–Dec) when Dungeness crab is traditionally consumed can cause substantial lost revenue opportunity for fishers at a time when price and demand are highest35,49. The fishery operates as a derby in which the majority of revenues are made in the first month of the fishery being open. The strong seasonal dynamics of the Dungeness crab fishery, largely driven by rapid depletion of legal sized crab, mean that the timing of management actions can have important impacts on fishing revenues. Across the fishery, based on observed vessel level revenues during the 2011–2018 baseline period, vessels earned an average of 62.33% (SD 24.04) of annual ex-vessel revenue during the first month of the season (15th Nov–15th Dec for the CMA/1st Dec–31st Dec for the NMA). After April 1st, vessels on average earn 10.54% (SD 18.98) of annual ex-vessel revenue. This average, based only on vessels that historically have actively participate past April 1st, (283 vessels in the NMA, 346 vessels in the CMA) rises to 20.36% (SD 13.37) of ex-vessel revenue. Thus, while the majority of the overall fisheries revenue is taken at the start of the season, an April 1st closure could still have a substantial impact on the revenues of active fishing vessels in the spring. Determination of economic risk for the fishery, at a minimum, should consider timing of closures in addition to total revenue losses, in order to quantify losses that will be felt at the individual vessel level. We suggest further research to investigate how closures affect different groups of fishers through stakeholder participation.Socio-economic impacts from whale mitigation measures could permeate into communities further than our analysis (based on ex-vessel revenue only) conveys35,36,37,49, and further investigation into these community level impacts is necessary to understand and sustain an equitable fishery supply chain even where there is no absolute revenue loss. Some of the communities influenced by whale entanglement mitigation in California rely heavily on ocean resources for employment, through fishing occupations but also through hospitality and tourism. Managing this issue in a way that minimizes the burden on resource dependent communities is strongly in line with the objectives set out in the UN Sustainable Development Goals (SDG’s), especially SDG 14 (life below water) but also related goals such as human well-being, reducing inequality and reducing the impacts of climate change56.Management ImplicationsBalancing socio-economic impacts against whale entanglement risk is challenging given the legally protected status of whale populations. However, potential economic losses reported here should motivate the development of mitigation measures (through cooperative innovation between industry, researchers and managers) that allow fishery production to be optimized whilst ensuring successful whale protection. At present, entire management areas, which constitute large regions of the coast, are closed in response to whale entanglement risk in California. Investigating how to minimize the spatiotemporal footprint of closures, such as by defining high risk zones dynamically based on fine-scale information of whale density and fishing effort, could provide an alternative mitigation structure. This could better consider the economic and conservation trade-offs while still being sensitive to changing environmental conditions. The introduction of dynamic zone closures, often broadly referred to as dynamic ocean management, has been demonstrated to reduce risk whilst minimizing lost fishing opportunities12,26,57,58, especially when environmental variability is high or species have a dynamic distribution59. Moreover, analysis of policy instruments to reduce whale entanglements with the American lobster fishery on the US Northeast coast found that economic costs of risk reduction could be 20% lower when mitigation decisions considered fishing opportunity costs alongside non-monetary benefits (biological risk), compared to non-monetary benefits alone12. This is promising for the implementation of such strategies in the California Current System.The caveat of this strategy is that dynamic zone closures require spatially and temporally explicit information on whale density and fishing effort which can be costly to attain. The use of ropeless gear has also been suggested as an alternative whale entanglement mitigation measure that requires further research and development before being initiated as an alternative regulatory tool60. The costs of monitoring or technical advancements however may outweigh the financial and societal cost of fishery closures. Revenue losses for Dungeness crab estimated here for the 2019 and 2020 seasons are on par with losses experienced during the HAB period. During the delays to the 2016 fishing season an estimated $26.1 million was lost from ex-vessel revenues from all species that crab fishers target, including $13.6 million from Dungeness crab alone38, requiring $25 million in government aid. Whale mitigation under the RAMP regulation will potentially delay or close the fishery year after year with uncertain economic impact that cannot be sustainably resolved with government aid. Development of tools to mitigate against economic loss while achieving whale protection will be necessary to come to a sustainable solution. This can only be achieved by first including economic loss in risk assessments. Doing so may also provide balance to partnerships between fishery managers and fishers.Regulators are obligated to protect Humpback whales, blue whales and Leatherback turtles using the best available science33. In this fishery, current triggers to open and close are based on a range of factors, but thus ultimately depend on the number of whales present within a management region33. Regulators have a number of alternative regulatory options available to them, which include depth restrictions, gear restrictions or modifications and fleet advisories, if they can offer the same level of whale protection33. Yet, the RAMP process lacks the socio-economic information needed to consider the socio-economic risk of regulatory actions, and that of the alternatives, to the fishing community. Results presented here highlight that the economic effects and that risk to fishing communities should be considered when designing whale entanglement mitigation programs33. Having this economic information will facilitate the ability of managers, as set out in the RAMP regulation (subsection d4)33, to consider the socio-economic impact if deciding between management measures that equivalently reduce entanglement risk.We have used two fishing seasons as an example of the economic impacts of these new whale entanglement regulations which will be implemented each year going forward. Synthesis of ex-vessel revenues is not a complete picture of the socio-economic impacts of regulations, but it provides a starting point for protecting both whales and fishing communities. While reported whale entanglements remain higher than pre-2014 totals, reported whale entanglements in California have declined markedly in the years following the 2014–2016 large marine heatwave (Fig. 1b). This is a success for this fishery and attributed to increased awareness, development of best practices for fishing gear and the mitigation program to protect whales. We now need to be successful at protecting and mitigating the socio-economic impacts to fishery participants and the fishing communities they support. More