More stories

  • in

    Genome-wide identification and expression profile of Elovl genes in threadfin fish Eleutheronema

    Identification of Elovl genes from E. tetradactylum and E. rhadinumTotally, we successfully identified 9 Elovl genes, including elovl1a, elovl1b, elovl4a, elovl4b, elovl5, elovl6, elovl6l, elovl7a, and elovl8b, both from E. tetradactylum and E. rhadinum genome (Table 2). In E. rhadinum, the shortest and the longest putative CDS length among all Elovl genes was 810 bp and 2019 bp, respectively. Their encoded protein size ranged from 269 amino acids to 672 amino acids. The theoretical molecular weight of Elovl proteins varied from 31061.48 to 75051.42 Da, with the theoretical isoelectric points (pI) ranging from 7.86 to 9.59. Most of the Elovl proteins were characterized as stable and hydrophilic proteins. Signal peptide prediction analysis showed that the elovl1b, elovl5, and elovl6 contained signal peptide sequences. In addition to elovl8b, all Elovl proteins contained transmembrane domains ranging from 5 to 7. Almost all Elovl proteins were predicted to be endoplasmic reticulum-located except elovl8b, predominantly localized in the nucleus.Table 2 Basic information for the Elovl gene family members.Full size tableIn E. tetradactylum, the putative CDS length of Elovl genes ranged from 810 to 1824 bp, and their encoded protein size ranged from 269 amino acids to 409 amino acids. The molecular weight of Elovl proteins varied from 31049.42 to 68750.14 Da, with the pI ranging from 8.72 to 9.64. Like Elovl proteins in E. rhadinum, most elovl proteins were predicted to be stable and hydrophilic. Signal peptide prediction analysis revealed that elovl1a, elovl5, and elovl6 had signal peptide sequence, which was different from E. rhadinum that elovl1b contained signal peptide sequence, but elovl1a did not. In addition, seven members showed the same number of transmembrane structures with E. rhadinum, while the elovl8b contained three and elovl4b contained seven transmembrane structures in contrast to E. rhadinum. The elovl8b was predicted to be localized in nuclear, while other members were localized in the endoplasmic reticulum, similar to E. rhadinum.Evolution of divergence and conservation of Elovl genesDivergence and conservation accompany the process of species evolution. To elucidate the phylogenetic relationship of Elovl genes among different species, a maximum like-hood tree was constructed on the basis of 18 Elovl genes in E. tetradactylum and E. rhadinum and 106 publicly available Elovl protein sequences. As shown in Fig. 1, these Elovl genes can be divided into eight subfamilies, including elovl1a/1b, elovl2, elovl3, elovl4a, elovl5, elovl6/6 l, elovl7a/7b, elovl8a/8b. However, 6 subfamilies were presented in the Eleutheronema genus, and there was only one subtype for elovl7 (elovl7a) and elovl8 (elovl8b) in E. tetradactylum and E. rhadinum. The elovl3 was mainly identified in mammalians such as Homo sapiens and Mus musculus, while a recent study reported a full repertoire of Elovl genes in the Colossoma macropomum genome, including elovl330. The loss of elovl2 occurred in the vast majority of marine fish lineages, which was only presented in a few fish species, such as C. carpio, D. rerio, S. salar, and S. grahami.Figure 1Phylogenetic tree for 18 Elovl proteins from E. tetradactylum and E. rhadinum, and 106 publicly available Elovl proteins from other species. All these proteins were aligned using ClustalW and then subjected to MEGAX for phylogenetic tree construction using the maximum like-hood method with 1000 replicates.Full size imageWe further performed the gene structure analysis to visualize the exon–intron structure of each gene, and the results revealed that the elovl8b had the largest intron number, while the elovl6/6 l subfamily genes contained three introns (Fig. 2a). Except for elovl8, Elovl genes belonging to the same subfamily shared a similar gene structure. Additionally, we identified ten motifs in Elovl genes, and the conversed motif types, numbers, and distributions in Elovl proteins were much more similar except for the elovl8b (Fig. 2b, TableS1). Two conserved motifs were found in the Elovl gene family except for elovl8b in E. rhadinum, which were related to the ELO domain via SMART evaluation analysis (Fig. 2c and d). Gene structural variation is important for gene evolution. In E. tetradactylum and E. rhadinum, Elovl genes showed similar gene structure, and the proteins shared similar motif compositions, indicating that the Elovl genes were highly conserved in the Eleutheronema genus.Figure 2Gene structure and conserved motifs diagram of Elovl genes. (a) Gene structure of Elovl genes. Exons were represented by pink boxes and introns by black lines; (b) Conserved motifs of Elovl proteins; (c and d) Logo representations of the ELO domains, motifs 1 and 2, respectively.Full size imageIn the process of evolution via natural selection, adaptation to certain environmental conditions likely drove the changes in endogenous capacity for LC-PUFA biosynthesis between marine and freshwater fishes31. The Elovl gene family has been functionally studied and characterized in a variety of fish species, and the member of the Elovl gene family of each species varied greatly. In the present study, for a comprehensive analysis of Elovl genes in the Eleutheronema genus, the Elovl gene ortholog clusters of mammals and various teleosts with different ecological niches and habitats were collected. The results showed that only seven Elovl genes (one gene for each subtype) were observed in mammals; however, more members were variably presented in teleosts, which might be related to the teleost-specific duplication. A previous study revealed that Sinocyclocheilus graham and C. carpio possessed the highest number of Elovl genes, containing 21 members of subtypes, resulting from an extra independent 4th whole-genome duplication event32, 33. Interestingly, only 9 Elovl genes were observed in Eleutheronema genus, the same as T. rubripes, possibly due to gene loss and the asymmetric acceleration of the evolutionary rate in one of the paralogs following the whole-genome duplication in some teleost fishes34. Additionally, the elovl2 and elovl3 were absent, but a novel subtype, elovl8, was present in most marine fishes. The elovl8, the most recently identified and novel active member of the Elovl protein family member, has been proposed to be a fish-specific elongase with two gene paralogs (elovl8a and elovl8b) described in teleost35. In Eleutheronema, we also found that the elovl8b was presented in E. tetradactylum and E. rhadinum, indicating the important roles in the LC-PUFAs biosynthesis of Eleutheronema fish. Similar results were also observed in rabbitfish and zebrafish20. The Elovl gene family member number in Eleutheronema genus is the same as T. rubripes, but less than I. punctatus (10), Gadus morhua (10), D. rerio (14), S. salar (18), and C. carpio (21), which might be due to the differential expansion events during the evolutions of fish species.Predicting the protein structure is a fundamental prerequisite for understanding the function and possible interactions of a protein. In the present study, the secondary structures as well as three-dimensional structures of Elovl proteins in both E. tetradactylum and E. rhadinum were predicted using the SOPMA and Phyre2 programs, respectively. The protein structures of all the candidate Elovl proteins were modeled at  > 90% confidence. The secondary structures of these proteins in E. tetradactylum revealed 40.86–50.30% alpha helixes, 28.10–28.10% random coil, 13.75–20.67% extended strand and 2.38–4.47% beta turn, while these ratios were predicted to be 47.55–53.27, 30.00–36.01, 6.99–18.12 and 2.38–4.75%, respectively, in E. rhadinum (Table 3). High ratio of alpha helixes and random coil in the Elovl protein structure might play important roles in fatty acids biosynthesis in fish, in accordance with the literature for the order Perciformes in Perca fluviatilis36. Additionally, the secondary structure pattern of Elovl proteins in the candidate E. tetradactylum and E. rhadinum species were highly similar (Fig. 3), indicating the probable similar biological functions as well as highly evolutionarily conserved Elovl genes in Eleutheronema species.Table 3 Properties of the secondary structures of Elovl proteins.Full size tableFigure 3The secondary structure pattern, including alpha helix (blue color), random coil (purple color), extended strand (red color), and beta turn (green color), of Elovl proteins in E. tetradactylum and E. rhadinum.Full size imageThe 3D model results showed that all predicted Elovl proteins had complex 3D structures, composing of multiple secondary structures including alpha-helices, random coils, and others (Fig. 4). The Elovl proteins of different subfamilies showed different 3D configurations. The 3D structures of Elovl proteins also revealed the presence of the conserved domain in each Elovl protein, which showed a typical three-dimensional frame comprising of various parallel alpha-helixes. To assay the quality and accuracy of the predicted 3D model for the candidate Elovl proteins, the Ramachandran plot analysis was employed (Figure S1). In model validation, the qualities of the Elovl proteins model varied from 90 to 98% based on the Ramachandran plot analysis, suggesting the reasonably good quality and reliability of the predicted 3D models. These results indicated that the predicted 3D model of Elovl proteins could provide valuable information for the further comprehensive studies of molecular function in the fatty acids biosynthesis in Eleutheronema species. Additionally, the comparisons between these structures in E. tetradactylum and E. rhadinum suggested that the Elovl proteins encompassed the conserved structures. In addition, gene duplication resulted in obvious 3D structural variation in the duplicated genes, such as Elovl4 (elovl4a and elovl4b), Elovl6 (elovl6 and elovl6l). The ascertained variations were revealed in duplicated Elovl proteins, and the diversities in these proteins structure may reflect their different obligations in the fatty acid biosynthesis and other biological processes.Figure 4Three-dimensional modeling of Elovl proteins in E. tetradactylum and E. rhadinum. All models have confidence levels above 90%.Full size imageTo explore the functional selection pressures acting on Elovl gene family, Ka, Ks, and Ka/Ks ratios were calculated for each gene. Generally, Ka/Ks  1 indicates positive selection. In this study, we found that all the Ka/Ks ratios for each gene were less than 0.5, suggesting that they were subjected to strong purifying selection during evolution, and their functions might be evolutionarily conserved (Fig. 5). Therefore, theoretically, the Elovl genes in the Eleutheronema genus had eliminated deleterious mutations in the population through purification selection. Similar results were also observed in Elovl gene family of Gymnocypris przewalskii that no positive selection trace was detected in most members except elovl211. Moreover, elovl6l and elovl8b showed a higher average Ka/Ks ratio than the other seven members, indicating that the evolution of elovl6l and elovl8b might be much less conservative and thereby could provide more variants for natural selection in Eleutheronema species.Figure 5The evolutionary rates of the Elovl genes in (a) E. tetradactylum and (b) E. rhadinum. The Ka, Ks, and Ka/Ks values were demonstrated in boxplots with error lines.Full size imageChromosomal location, collinearity, and protein–protein interaction network analysis of Elovl genesAs shown in Fig. 6a and b, Elovl genes were randomly and unevenly distributed on seven chromosomes in both E. tetradactylum and E. rhadinum, including Chr5, Chr6, Chr8, Chr10, Chr11, Chr13, and Chr25. The Chr5 and Chr6 harbored two Elovl genes (elovl1b and elovl8b in Chr5, elovl5 and elovl6l in Chr6), while other chromosomes each carried a single Elovl gene. Collinearity relationship analysis was performed to further investigate the gene duplication events within the Elovl gene family. The results revealed that a pair of segmental duplication genes (elovl4a/4b) showed collinear relationships. A chromosome-wide collinearity analysis also showed that the chromosomes were highly homologous between E. tetradactylum and E. rhadinum, including the Elovl gene family (Figure S2). To infer the protein interaction within Elovl gene family, we constructed the protein–protein interaction (PPI) network of the Elovl proteins based on the interaction relationship of the homologous Elovl proteins in zebrafish. The results showed that Elovl genes had close interaction with other members except for the elovl4a/4b and elovl8b (Fig. 6c), which suggested that they might participate in diverse functions by interacting with other proteins. Thus far, elovl4a and elovl4b were widely identified in most fish, which could effectively elongate PUFA substrates37. In addition, the elovl4a/4b were identified to be homologous proteins of zebrafish, indicating that the elovl4 subtype was highly conserved during evolution and played important roles in the biosynthesis of LC-PUFA in Eleutheronema.Figure 6Chromosomal location and collinearity analysis of Elovl gene family in (a) E. tetradactylum and (b) E. rhadinum. Colored boxes represented chromosomes. Segmental duplication genes are connected with grey lines; (c) a protein–protein interaction network for Elovl genes based on their orthologs in zebrafish.Full size imageExpression patterns of ELOVL genes in different tissuesIn the present study, we aimed to determine the expression patterns and gained insights into the potential functions of Elovl genes in the brain, eye, gill, heart, kidney, liver, muscle, stomach, and intestine. The expression patterns of Elovl genes in different tissues and species were distinct, suggesting the diverse roles during fish development (Fig. 7a and b). In our present study, the elovl1a and elvovl1b were expressed in a relatively narrow range of tissues, including the liver, stomach, and intestine. Some Elovl genes had much higher relative expression rates, e.g., elovl1a and elovl7a. The elovl4a was primarily distributed in the brain and eye, slightly expressed in gills while hardly detectable in other tissues, consistent with previous studies37, 38, which might play an important role in endogenous biosynthesis of LC-PUFA in the neural system of fish. In contrast to elovl4a, elovl4b was ubiquitously, instead of tissue-specific, expressed in most tissues while hardly examined in the heart and kidney. The elovl4a and elovl4b were two commonly paralogues in evolutionarily diverged fish species, and the striking difference in expression patterns between elovl4a and elovl4b might be due to the potential functional divergence of these two paralogues. In addition, elovl8b, the novel active member of the Elovl protein family, was expressed in several tissues, suggesting the essential roles in LC-PUFAs biosynthesis of teleost as indicated by a previous study20. Moreover, the differences in expression patterns among different Elovl genes indicated that these genes might possibly undergo functional divergence during evolution in the Eleutheronema genus. Overall, our present study firstly provided the preliminary organ-specific expression data of the Elovl gene family in E. tetradactylum and E. rhadinum, which could provide the foundation for further clarifying the function of these genes in the evolutionary development of the Eleutheronema genus.Figure 7qPCR assessment of tissue distribution of elovl1a, elovl1b, elovl4a, elovl4b, elovl5, elovl6, elovl6l, elovl7a, and elovl8b gene expression in (a) E. tetradactylum and (b) E. rhadinum for various tissues including the brain, eye, gill, heart, kidney, liver, muscle, stomach, and intestine.Full size image More

  • in

    Localized coevolution between microbial predator and prey alters community-wide gene expression and ecosystem function

    Ehrlich PR, Raven PH. Butterflies and plants: a study in coevolution. Evolution. 1964;18:586–608.Article 

    Google Scholar 
    Marston MF, Pierciey FJ Jr, Shepard A, Gearin G, Qi J, Yandava C, et al. Rapid diversification of coevolving marine Synechococcus and a virus. Proc Natl Acad Sci USA. 2012;109:4544–9.Article 
    CAS 

    Google Scholar 
    Hall AR, Scanlan PD, Buckling A. Bacteria-phage coevolution and the emergence of generalist pathogens. Am Nat. 2011;177:44–53.Article 

    Google Scholar 
    Badger MR, Andrews TJ, Whitney SM, Ludwig M, Yellowlees DC, Leggat W, et al. The diversity and coevolution of Rubisco, plastids, pyrenoids, and chloroplast-based CO2-concentrating mechanisms in algae. Can J Bot. 1998;76:1052–71.CAS 

    Google Scholar 
    Schluter D. The ecology of adaptive radiation. Oxford, UK: University Press; 2000.Buckling A, Maclean CR, Brockhurst MA, Colegrave N. The Beagle in a bottle. Nature. 2009;457:824–9.Article 
    CAS 

    Google Scholar 
    Thompson JN. The coevolutionary process. Chicago, USA: University of Chicago Press; 1994.Vallina SM, Follows MJ, Dutkiewicz S, Montoya JM, Cermeno P, Loreau M. Global relationship between phytoplankton diversity and productivity in the ocean. Nat Commun. 2014;5:4299.Article 
    CAS 

    Google Scholar 
    Jürgens K, Matz C. Predation as a shaping force for the phenotypic and genotypic composition of planktonic bacteria. Antonie Van Leeuwenhoek. 2002;81:413–34.Article 

    Google Scholar 
    Wildschutte H, Wolfe DM, Tamewitz A, Lawrence JG. Protozoan predation, diversifying selection, and the evolution of antigenic diversity in Salmonella. Proc Natl Acad Sci USA. 2004;101:10644–9.Article 
    CAS 

    Google Scholar 
    Thompson JN. The geographic mosaic of coevolution. Chicago, USA: University of Chicago Press; 2005.Hahn MW, Höfle MG. Grazing of protozoa and its effect on populations of aquatic bacteria. FEMS Microbiol Ecol. 2001;35:113–21.Article 
    CAS 

    Google Scholar 
    Fuhrman JA, Noble RT. Viruses and protists cause similar bacterial mortality in coastal seawater. Limnol Oceanogr. 1995;40:1236–42.Article 

    Google Scholar 
    Azam F, Fenchel T, Field JG, Gray JS, Meyer-Reil LA, Thingstad F. The ecological role of water-column microbes in the sea. Mar Ecol Prog Ser. 1983;10:257–63.Article 

    Google Scholar 
    Lankau RA, Strauss SY. Mutual feedbacks maintain both genetic and species diversity in a plant community. Science. 2007;317:1561–3.Article 
    CAS 

    Google Scholar 
    Hogle SL, Hepolehto I, Ruokolainen L, Cairns J, Hiltunen T. Effects of phenotypic variation on consumer coexistence and prey community structure. Ecol Lett. 2022:25;307–19.Yoshida T, Jones LE, Ellner SP, Fussmann GF, Hairston NG Jr. Rapid evolution drives ecological dynamics in a predator-prey system. Nature. 2003;424:303–6.Article 
    CAS 

    Google Scholar 
    McClean D, McNally L, Salzberg LI, Devine KM, Brown SP, Donohue I. Single gene locus changes perturb complex microbial communities as much as apex predator loss. Nat Commun. 2015;6:8235.Article 

    Google Scholar 
    Gómez P, Paterson S, De Meester L, Liu X, Lenzi L, Sharma MD, et al. Local adaptation of a bacterium is as important as its presence in structuring a natural microbial community. Nat Commun. 2016;7:12453.Article 

    Google Scholar 
    Middelboe M, Holmfeldt K, Riemann L, Nybroe O, Haaber J. Bacteriophages drive strain diversification in a marine Flavobacterium: Implications for phage resistance and physiological properties. Environ Microbiol. 2009;11:1971–82.Article 
    CAS 

    Google Scholar 
    Lennon JT, Martiny JBH. Rapid evolution buffers ecosystem impacts of viruses in a microbial food web. Ecol Lett. 2008;11:1178–88.Article 

    Google Scholar 
    Cairns J, Jokela R, Hultman J, Tamminen M, Virta M, Hiltunen T. Construction and characterization of synthetic bacterial community for experimental ecology and evolution. Front Genet. 2018;9:312.Article 

    Google Scholar 
    Pascual-García A, Bell T. Community-level signatures of ecological succession in natural bacterial communities. Nat Commun. 2020;11:2386.Article 

    Google Scholar 
    Rivett DW, Bell T. Abundance determines the functional role of bacterial phylotypes in complex communities. Nat Microbiol. 2018;3:767–72.Article 
    CAS 

    Google Scholar 
    Cairns J, Moerman F, Fronhofer EA, Altermatt F, Hiltunen T. Evolution in interacting species alters predator life-history traits, behaviour and morphology in experimental microbial communities. Proc Biol Sci. 2020;287:20200652.
    Google Scholar 
    Cooke DP, Wedge DC, Lunter G. A unified haplotype-based method for accurate and comprehensive variant calling. Nat Biotechnol. 2021;39:885–92.Article 
    CAS 

    Google Scholar 
    Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012;6:80–92.Article 
    CAS 

    Google Scholar 
    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.Article 
    CAS 

    Google Scholar 
    Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, et al. Fast genome-wide functional annotation through orthology assignment by eggNOG-Mapper. Mol Biol Evol. 2017;34:2115–22.Article 
    CAS 

    Google Scholar 
    Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47:D309–14.Article 
    CAS 

    Google Scholar 
    Blin K, Shaw S, Steinke K, Villebro R, Ziemert N, Lee SY, et al. antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019;47:W81–7.Article 
    CAS 

    Google Scholar 
    Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. The dynamics of molecular evolution over 60,000 generations. Nature. 2017;551:45–50.Article 

    Google Scholar 
    Timonen J, Mannerström H, Vehtari A, Lähdesmäki H. lgpr: an interpretable nonparametric method for inferring covariate effects from longitudinal data. Bioinformatics. 2021;37:1860–7.Article 
    CAS 

    Google Scholar 
    Willis AD, Martin BD. Estimating diversity in networked ecological communities. Biostatistics. 2022;23:207–22.Article 

    Google Scholar 
    Willis A, Bunge J, Whitman T. Improved detection of changes in species richness in high diversity microbial communities. J R Stat Soc C. 2017;66:963–77.Article 

    Google Scholar 
    Anderson MJ. A new method for non‐parametric multivariate analysis of variance. Austral Ecol. 2001:26;32–46.Anderson MJ. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62:245–53.Article 

    Google Scholar 
    Martin BD, Witten D, Willis AD. Modeling microbial abundances and dysbiosis with beta-binomial regression. Ann Appl Stat. 2020;14:94–115.Article 

    Google Scholar 
    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.Article 

    Google Scholar 
    Zhang Y, Thompson KN, Huttenhower C, Franzosa EA. Statistical approaches for differential expression analysis in metatranscriptomics. Bioinformatics. 2021;37:i34–41.Article 
    CAS 

    Google Scholar 
    Abdi H, Williams LJ, Valentin D, Bennani-Dosse M. STATIS and DISTATIS: optimum multitable principal component analysis and three way metric multidimensional scaling. WIREs Comp Stat. 2012;4:124–67.Charrad M, Ghazzali N, Boiteau V, Niknafs A. NbClust: an R package for determining the relevant number of clusters in a data set. J Stat Softw. 2014;61:1–36.Article 

    Google Scholar 
    Zhu A, Ibrahim JG, Love MI. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. 2019;35:2084–92.Article 
    CAS 

    Google Scholar 
    Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.Article 
    CAS 

    Google Scholar 
    Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.Article 
    CAS 

    Google Scholar 
    Wichman HA, Badgett MR, Scott LA, Boulianne CM, Bull JJ. Different trajectories of parallel evolution during viral adaptation. Science. 1999;285:422–4.Article 
    CAS 

    Google Scholar 
    Lieberman TD, Michel J-B, Aingaran M, Potter-Bynoe G, Roux D, Davis MR Jr, et al. Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes. Nat Genet. 2011;43:1275–80.Article 
    CAS 

    Google Scholar 
    Lebeuf-Taylor E, McCloskey N, Bailey SF, Hinz A, Kassen R. The distribution of fitness effects among synonymous mutations in a gene under directional selection. Elife. 2019;8:e45952.Article 

    Google Scholar 
    Bailey SF, Hinz A, Kassen R. Adaptive synonymous mutations in an experimentally evolved Pseudomonas fluorescens population. Nat Commun. 2014;5:4076.Article 
    CAS 

    Google Scholar 
    Mukherjee S, Jemielita M, Stergioula V, Tikhonov M, Bassler BL. Photosensing and quorum sensing are integrated to control Pseudomonas aeruginosa collective behaviors. PLoS Biol. 2019;17:e3000579.Article 
    CAS 

    Google Scholar 
    Segura A, Hurtado A, Duque E, Ramos JL. Transcriptional phase variation at the flhB gene of Pseudomonas putida DOT-T1E is involved in response to environmental changes and suggests the participation of the flagellar export system in solvent tolerance. J Bacteriol. 2004;186:1905–9.Article 
    CAS 

    Google Scholar 
    Lee X, Reimmann C, Greub G, Sufrin J, Croxatto A. The Pseudomonas aeruginosa toxin L-2-amino-4-methoxy-trans-3-butenoic acid inhibits growth and induces encystment in Acanthamoeba castellanii. Microbes Infect. 2012;14:268–72.Article 
    CAS 

    Google Scholar 
    Montagnes DJS, Barbosa AB, Boenigk J, Davidson K, Jürgens K, Macek M, et al. Selective feeding behaviour of key free-living protists: avenues for continued study. Aquat Micro Ecol. 2008;53:83–98.Article 

    Google Scholar 
    Collins K, editor. Tetrahymena thermophila. New York: Academic Press, Elsevier; 2012.Ruehle MD, Orias E, Pearson CG. Tetrahymena as a unicellular model eukaryote: genetic and genomic tools. Genetics. 2016;203:649–65.Article 
    CAS 

    Google Scholar 
    Plum K, Tarkington J, Zufall RA. Experimental evolution in Tetrahymena. Microorganisms. 2022;10:1–11.Heard E, Martienssen RA. Transgenerational epigenetic inheritance: myths and mechanisms. Cell. 2014;157:95–109.Article 
    CAS 

    Google Scholar 
    Jones ML, Rivett DW, Pascual-García A, Bell T. Relationships between community composition, productivity and invasion resistance in semi-natural bacterial microcosms. Elife. 2021;10:1–25.Kertesz MA. Riding the sulfur cycle-metabolism of sulfonates and sulfate esters in gram-negative bacteria. FEMS Microbiol Rev. 2000;24:135–75.CAS 

    Google Scholar 
    Park C, Shin B, Park W. Protective role of bacterial alkanesulfonate monooxygenase under oxidative stress. Appl Environ Microbiol. 2020;86:1–14.Shatalin K, Shatalina E, Mironov A, Nudler E. H2S: a universal defense against antibiotics in bacteria. Science. 2011;334:986–90.Article 
    CAS 

    Google Scholar 
    Ong C-LY, Beatson SA, Totsika M, Forestier C, McEwan AG, Schembri MA. Molecular analysis of type 3 fimbrial genes from Escherichia coli, Klebsiella and Citrobacter species. BMC Microbiol. 2010;10:183.Article 

    Google Scholar 
    McNally L, Brown SP. Building the microbiome in health and disease: niche construction and social conflict in bacteria. Philos Trans R Soc Lond B Biol Sci. 2015;370:1–8.Scheuerl T, Cairns J, Becks L, Hiltunen T. Predator coevolution and prey trait variability determine species coexistence. Proc Biol Sci. 2019;286:20190245.
    Google Scholar 
    Koskella B, Brockhurst MA. Bacteria-phage coevolution as a driver of ecological and evolutionary processes in microbial communities. FEMS Microbiol Rev. 2014;38:916–31.Article 
    CAS 

    Google Scholar 
    Wilhelm Scherer H. Sulfur in soils. J Plant Nutr Soil Sci. 2009;172:326–35.Article 

    Google Scholar 
    Kaya K. Chemistry and biochemistry of taurolipids. Prog Lipid Res. 1992;31:87–108.Article 
    CAS 

    Google Scholar 
    Fukami T. Historical contingency in community assembly: integrating niches, species pools, and priority effects. Annu Rev Ecol Evol Syst. 2015;46:1–23.Article 

    Google Scholar 
    Debray R, Herbert RA, Jaffe AL, Crits-Christoph A, Power ME, Koskella B. Priority effects in microbiome assembly. Nat Rev Microbiol. 2022;20:109–21.Article 
    CAS 

    Google Scholar 
    Price TD, Qvarnström A, Irwin DE. The role of phenotypic plasticity in driving genetic evolution. Proc Biol Sci. 2003;270:1433–40.Article 

    Google Scholar  More

  • in

    Human-mediated dispersal drives the spread of the spotted lanternfly (Lycorma delicatula)

    Simberloff, D. et al. (eds) Invasive Species in a Globalized World (University of Chicago Press, 2015).
    Google Scholar 
    Gippet, J. M., Liebhold, A. M., Fenn-Moltu, G. & Bertelsmeier, C. Human-mediated dispersal in insects. Curr. Opin. Insect Sci. 35, 96–102 (2019).Article 

    Google Scholar 
    Hall, C. M. Biological invasion, biosecurity, tourism, and globalisation. In Handbook of Globalisation and Tourism (Edward Elgar Publishing, 2019).
    Google Scholar 
    Bertelsmeier, C. Globalization and the anthropogenic spread of invasive social insects. Curr. Opin. Insect Sci. https://doi.org/10.1016/j.cois.2021.01.006 (2021).Article 

    Google Scholar 
    Simberloff, D. How common are invasion-induced ecosystem impacts?. Biol. Invasions 13, 1255–1268 (2011).Article 

    Google Scholar 
    Hayes, K. R. & Barry, S. C. Are there any consistent predictors of invasion success?. Biol. Invasions 10, 483–506 (2008).Article 

    Google Scholar 
    Catford, J. A., Vesk, P. A., Richardson, D. M. & Pyšek, P. Quantifying levels of biological invasion: Towards the objective classification of invaded and invasible ecosystems. Glob. Change Biol. 18, 44–62 (2012).Article 
    ADS 

    Google Scholar 
    Arim, M., Abades, S. R., Neill, P. E., Lima, M. & Marquet, P. A. Spread dynamics of invasive species. Proc. Natl. Acad. Sci. USA 103, 374–378 (2006).Article 
    ADS 
    CAS 

    Google Scholar 
    Kamenova, S. et al. Invasions toolkit: Current methods for tracking the spread and impact of invasive species. Adv. Ecol. Res. 56, 85–182 (2017).Article 

    Google Scholar 
    Hulme, P. E. Trade, transport and trouble: Managing invasive species pathways in an era of globalization. J. Appl. Ecol. 46, 10–18 (2009).Article 

    Google Scholar 
    Banks, N. C., Paini, D. R., Bayliss, K. L. & Hodda, M. The role of global trade and transport network topology in the human-mediated dispersal of alien species. Ecol. Lett. 18, 188–199 (2015).Article 

    Google Scholar 
    Crooks, J. A. & Rilov, G. The establishment of invasive species. In Biological Invasions in Marine Ecosystems 173–175 (Springer, 2009).Chapter 

    Google Scholar 
    Lockwood, J. L., Cassey, P. & Blackburn, T. M. The more you introduce the more you get: The role of colonization pressure and propagule pressure in invasion ecology. Divers. Distrib. 15, 904–910 (2009).Article 

    Google Scholar 
    Sakai, A. K. et al. The population biology of invasive species. Annu. Rev. Ecol. Syst. 32, 305–332 (2001).Article 

    Google Scholar 
    O’Reilly-Nugent, A. et al. Landscape effects on the spread of invasive species. Curr. Landsc. Ecol. Rep. 1, 107–114 (2016).Article 

    Google Scholar 
    Simberloff, D. We can eliminate invasions or live with them. Successful management projects. In Ecological Impacts of Non-native Invertebrates and Fungi on Terrestrial Ecosystems 149–157 (Springer, 2008).
    Google Scholar 
    Gutierrez, A. P. & Ponti, L. Eradication of invasive species: Why the biology matters. Environ. Entomol. 42, 395–411 (2013).Article 

    Google Scholar 
    McLaughlin, G. M. & Dearden, P. K. Invasive insects: Management methods explored. J. Insect Sci. 19, 17 (2019).Article 

    Google Scholar 
    Han, J. M. et al. Lycorma delicatula (hemiptera: Auchenorrhyncha: Fulgoridae: Aphaeninae) finally, but suddenly arrived in Korea. Entomol. Res. 38, 281–286 (2008).Article 

    Google Scholar 
    Park, J.-D. et al. Biological characteristics of lycorma delicatula and the control effects of some insecticides. Korean J. Appl. Entomol. 48, 53–57 (2009).Article 

    Google Scholar 
    Shin, Y.-H., Moon, S.-R., Yoon, C.-M., Ahn, K.-S. & Kim, G.-H. Insecticidal activity of 26 insectcides against eggs and nymphs of Lycorma delicatula (hemiptera: Fulgoridae). Korean J. Pestic. Sci. 14, 157–163 (2010).
    Google Scholar 
    Dara, S. K., Barringer, L. & Arthurs, S. P. Lycorma delicatula (hemiptera: Fulgoridae): A new invasive pest in the United States. J. Integr. Pest Manag. 6, 20 (2015).Article 

    Google Scholar 
    Urban, J. M. Perspective: Shedding light on spotted lanternfly impacts in the USA. Pest Manag. Sci. 76, 10–17 (2020).Article 
    CAS 

    Google Scholar 
    Liu, G. Some extracts from the history of entomology in china. Psyche 46, 23–28 (1939).Article 

    Google Scholar 
    Barringer, L. E., Donovall, L. R., Spichiger, S.-E., Lynch, D. & Henry, D. The first new world record of Lycorma delicatula (insecta: Hemiptera: Fulgoridae). Entomol. News 125, 20–23 (2015).Article 

    Google Scholar 
    Parra, G., Moylett, H. & Bulluck, R. Technical Working Group Summary Report: Spotted Lanternfly, Lycorma Delicatula (White, 1845). (2018).Harper, J. K., Stone, W., Kelsey, T. W. & Kime, L. F. Potential Economic Impact of the Spotted Lanternfly on Agriculture and Forestry in Pennsylvania 1–84 (The Center for Rural Pennsylvania, 2019).
    Google Scholar 
    Kim, J. G., Lee, E.-H., Seo, Y.-M. & Kim, N.-Y. Cyclic behavior of Lycorma delicatula (insecta: Hemiptera: Fulgoridae) on host plants. J. Insect Behav. 24, 423–435 (2011).Article 

    Google Scholar 
    Albright, T. A. et al. Pennsylvania forests 2014. Resour. Bull. 111, 1–140 (2017).
    Google Scholar 
    Liu, H. Oviposition substrate selection, egg mass characteristics, host preference, and life history of the spotted lanternfly (hemiptera: Fulgoridae) in North America. Environ. Entomol. 48, 1452–1468 (2019).
    Google Scholar 
    Barringer, L. & Ciafré, C. M. Worldwide feeding host plants of spotted lanternfly, with significant additions from North America. Environ. Entomol. 49, 999–1011 (2020).Article 

    Google Scholar 
    Murman, K. et al. Distribution, survival, and development of spotted lanternfly on host plants found in North America. Environ. Entomol. 49, 1270–1281 (2020).Article 

    Google Scholar 
    Huron, N. A., Behm, J. E. & Helmus, M. R. Paninvasion severity assessment of a us grape pest to disrupt the global wine market. bioRxiv (2021).Dara, S. K. Update on the Spotted Lanternfly.Jung, J.-M., Jung, S., Byeon, D.-H. & Lee, W.-H. Model-based prediction of potential distribution of the invasive insect pest, spotted lanternfly Lycorma delicatula (hemiptera: Fulgoridae), by using climex. J. Asia-Pac. Biodivers. 10, 532–538 (2017).Article 

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

    Google Scholar 
    Wakie, T. T., Neven, L. G., Yee, W. L. & Lu, Z. The establishment risk of Lycorma delicatula (hemiptera: Fulgoridae) in the United States and globally. J. Econ. Entomol. 113, 306–314 (2020).
    Google Scholar 
    Grimm, V. et al. A standard protocol for describing individual-based and agent-based models. Ecol. Model. 198, 115–126 (2006).Article 

    Google Scholar 
    DeAngelis, D. L. Individual-Based Models and Approaches in Ecology: Populations, Communities and Ecosystems (CRC Press, 2018).Book 

    Google Scholar 
    Łomnicki, A. Individual-based models and the individual-based approach to population ecology. Ecol. Model. 115, 191–198 (1999).Article 

    Google Scholar 
    Grimm, V. & Railsback, S. F. A conceptual framework for designing individual-based models. In Individual-Based Modeling and Ecology 71–121 (Princeton University Press, 2005).Chapter 
    MATH 

    Google Scholar 
    Smith, N. R. et al. Agent-based models of malaria transmission: A systematic review. Malar. J. 17, 1–16 (2018).Article 
    CAS 

    Google Scholar 
    Venkatramanan, S. et al. Using data-driven agent-based models for forecasting emerging infectious diseases. Epidemics 22, 43–49 (2018).Article 

    Google Scholar 
    Harris, C. M., Park, K. J., Atkinson, R., Edwards, C. & Travis, J. Invasive species control: Incorporating demographic data and seed dispersal into a management model for rhododendron ponticum. Ecol. Inform. 4, 226–233 (2009).Article 

    Google Scholar 
    Gallien, L., Münkemüller, T., Albert, C. H., Boulangeat, I. & Thuiller, W. Predicting potential distributions of invasive species: Where to go from here?. Divers. Distrib. 16, 331–342 (2010).Article 

    Google Scholar 
    Rebaudo, F., Crespo-Pérez, V., Silvain, J.-F. & Dangles, O. Agent-based modeling of human-induced spread of invasive species in agricultural landscapes: Insights from the potato moth in ecuador. J. Artif. Soc. Soc. Simul. 14, 7 (2011).Article 

    Google Scholar 
    Day, C. C., Landguth, E. L., Bearlin, A., Holden, Z. A. & Whiteley, A. R. Using simulation modeling to inform management of invasive species: A case study of eastern brook trout suppression and eradication. Biol. Conserv. 221, 10–22 (2018).Article 

    Google Scholar 
    Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259 (2006).Article 

    Google Scholar 
    Phillips, S. J., Dudı’k, M. & Schapire, R. E. A maximum entropy approach to species distribution modeling. In Proceedings of the Twenty-first International Conference on Machine Learning 83 (2004).Phillips, S. J. et al. A brief tutorial on maxent. AT&T Res. 190, 231–259 (2005).
    Google Scholar 
    R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2020).
    Google Scholar 
    Urbanek, S. RJava: Low-Level R to Java Interface. (2020).Hijmans, R. J., Phillips, S., Leathwick, J., Elith, J. & Hijmans, M. R. J. Package ‘dismo’. Circles 9, 1–68 (2017).
    Google Scholar 
    Elith, J. et al. A statistical explanation of maxent for ecologists. Divers. Distrib. 17, 43–57 (2011).Article 

    Google Scholar 
    Lane, M. A. & Edwards, J. L. The global biodiversity information facility (gbif). Syst. Assoc. Spec. 73, 1 (2007).
    Google Scholar 
    O’Donnell, M. S. & Ignizio, D. A. Bioclimatic predictors for supporting ecological applications in the conterminous united states. US Geol. Surv. Data Ser. 691, 4–9 (2012).
    Google Scholar 
    Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

    Google Scholar 
    Hijmans, R. J. Raster: Geographic Data Analysis and Modeling. (2020).Venter, O. et al. Last of the wild project, version 3 (lwp-3): 2009 human footprint, 2018 release. NASA Socioeconomic Data and Applications Center (SEDAC) 10, H46T40JQ44 (2018).Park, M. Overwintering ecology and population genetics of Lycorma delicatula (hemiptera: Fulgoridae) in Korea. Seoul National University, Seoul, Korea Doctoral Thesis (2015).Pearson, K. I. Mathematical contributions to the theory of evolution. VII. On the correlation of characters not quantitatively measurable. Philos. Trans. R. Soc. Lond. Ser. A 195, 1–47 (1900).ADS 
    MATH 

    Google Scholar 
    Warmerdam, F. The geospatial data abstraction library. In Open Source Approaches in Spatial Data Handling 87–104 (Springer, 2008).Chapter 

    Google Scholar 
    Greenberg, J. A., Mattiuzzi, M. & SystemRequirements, G. Package ‘gdalUtils’. (2020).Domingue, M. J. & Baker, T. C. Orientation of flight for physically disturbed spotted lanternflies, Lycorma delicatula, (Hemiptera, fulgoridae). J. Asia-Pac. Entomol. 22, 117–120 (2019).Article 

    Google Scholar 
    Myrick, A. J. & Baker, T. C. Analysis of anemotactic flight tendencies of the spotted lanternfly (Lycorma delicatula) during the 2017 mass dispersal flights in pennsylvania. J. Insect Behav. 32, 11–23 (2019).Article 

    Google Scholar 
    Wolfin, M. S., Myrick, A. J. & Baker, T. C. Flight duration capabilities of dispersing adult spotted lanternflies, Lycorma delicatula. J. Insect Behav. 33, 125–137 (2020).Article 

    Google Scholar 
    Strömbom, D. & Pandey, S. Modeling the life cycle of the spotted lanternfly (Lycorma delicatula) with management implications. Math. Biosci. 340, 108670 (2021).Article 
    MATH 

    Google Scholar 
    Wellington, W. G. Conditions governing the distribution of insects in the free atmosphere. Can. Entomol. 77, 7–15 (1945).Article 

    Google Scholar 
    DeLong, D. M. The bionomics of leafhoppers. Annu. Rev. Entomol. 16, 179–210 (1971).Article 

    Google Scholar 
    Baker, T. et al. Progression of seasonal activities of adults of the spotted lanternfly, Lycorma delicatula, during the 2017 season of mass flight dispersal behavior in eastern Pennsylvania. J. Asia-Pac. Entomol. 22, 705–713 (2019).Article 

    Google Scholar 
    Leach, H. & Leach, A. Seasonal phenology and activity of spotted lanternfly (Lycorma delicatula) in eastern us vineyards. J. Pest Sci. 93, 1215–1224 (2020).Article 

    Google Scholar 
    Goutte, C. & Gaussier, E. A probabilistic interpretation of precision, recall and f-score, with implication for evaluation. In European Conference on Information Retrieval 345–359 (Springer, 2005).
    Google Scholar 
    Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).Article 

    Google Scholar 
    Tukey, J. Multiple comparisons. J. Am. Stat. Assoc. 48, 624–625 (1953).
    Google Scholar 
    Mendiburu, F. de & Mendiburu, M. F. de. Package ‘agricolae’. R Package, Version 1-2 (2019).McAvoy, T. J., Snyder, A. L., Johnson, N., Salom, S. M. & Kok, L. T. Road survey of the invasive tree-of-heaven (Ailanthus altissima) in Virginia. Invasive Plant Sci. Manag. 5, 506–512 (2012).Article 

    Google Scholar 
    Casella, F. & Vurro, M. Ailanthus altissima (tree of heaven): Spread and harmfulness in a case-study urban area. Arboricult. J. 35, 172–181 (2013).Article 

    Google Scholar 
    Takahashi, D. & Park, Y.-S. Spatial heterogeneities of human-mediated dispersal vectors accelerate the range expansion of invaders with source–destination-mediated dispersal. Sci. Rep. 10, 1–9 (2020).Article 

    Google Scholar 
    Meijer, J. R., Huijbregts, M. A., Schotten, K. C. & Schipper, A. M. Global patterns of current and future road infrastructure. Environ. Res. Lett. 13, 064006 (2018).Article 
    ADS 

    Google Scholar 
    Turner, R. M. et al. Worldwide border interceptions provide a window into human-mediated global insect movement. Ecol. Appl. 31, e02412 (2021).Article 

    Google Scholar 
    Ricciardi, A. Are modern biological invasions an unprecedented form of global change?. Conserv. Biol. 21, 329–336 (2007).Article 

    Google Scholar 
    Wilson, J. R., Dormontt, E. E., Prentis, P. J., Lowe, A. J. & Richardson, D. M. Something in the way you move: Dispersal pathways affect invasion success. Trends Ecol. Evol. 24, 136–144 (2009).Article 

    Google Scholar 
    Auffret, A. G., Berg, J. & Cousins, S. A. The geography of human-mediated dispersal. Divers. Distrib. 20, 1450–1456 (2014).Article 

    Google Scholar 
    Koch, F. H., Yemshanov, D., Magarey, R. D. & Smith, W. D. Dispersal of invasive forest insects via recreational firewood: A quantitative analysis. J. Econ. Entomol. 105, 438–450 (2012).Article 

    Google Scholar 
    Eyer, P.-A. et al. Extensive human-mediated jump dispersal within and across the native and introduced ranges of the invasive termite Reticulitermes flavipes. Mol. Ecol. 30, 3948–3964 (2020).Article 

    Google Scholar 
    Petrice, T. R. & Haack, R. A. Effects of cutting date, outdoor storage conditions, and splitting on survival of Agrilus planipennis (coleoptera: Buprestidae) in firewood logs. J. Econ. Entomol. 99, 790–796 (2006).Article 

    Google Scholar 
    Petrice, T. R. & Haack, R. A. Can emerald ash borer, Agrilus planipennis (coleoptera: Buprestidae), emerge from logs two summers after infested trees are cut?. Great Lakes Entomol. 40, 92–95 (2007).
    Google Scholar 
    Muirhead, J. R. et al. Modelling local and long-distance dispersal of invasive emerald ash borer Agrilus planipennis (coleoptera) in North America. Divers. Distrib. 12, 71–79 (2006).Article 

    Google Scholar 
    Güneralp, B., Reba, M., Hales, B. U., Wentz, E. A. & Seto, K. C. Trends in urban land expansion, density, and land transitions from 1970 to 2010: A global synthesis. Environ. Res. Lett. 15, 044015 (2020).Article 
    ADS 

    Google Scholar 
    Hulme, P. E. Unwelcome exchange: International trade as a direct and indirect driver of biological invasions worldwide. One Earth 4, 666–679 (2021).Article 
    ADS 

    Google Scholar  More

  • in

    Human fingerprint on structural density of forests globally

    Watson, J. E. M. et al. The exceptional value of intact forest ecosystems. Nat. Ecol. Evol. 2, 599–610 (2018).Article 

    Google Scholar 
    Potapov, P. et al. The last frontiers of wilderness: tracking loss of intact forest landscapes from 2000 to 2013. Sci. Adv. https://doi.org/10.1126/sciadv.1600821 (2017).Matricardi, E. A. T. et al. Long-term forest degradation surpasses deforestation in the Brazilian Amazon. Science 369, 1378–1382 (2020).Article 
    CAS 

    Google Scholar 
    Venter, O. et al. Sixteen years of change in the global terrestrial human footprint and implications for biodiversity conservation. Nat. Commun. 7, 12558 (2016).Article 
    CAS 

    Google Scholar 
    Grantham, H. S. et al. The emerging threat of extractives sector to intact forest landscapes. Front. For. Glob. Change https://doi.org/10.3389/ffgc.2021.692338 (2021).IPBES: Summary for Policymakers. In The Global Assessment Report on Biodiversity and Ecosystem Services (eds Díaz, S. et al.) (IPBES, 2019).Qin, Y. et al. Carbon loss from forest degradation exceeds that from deforestation in the Brazilian Amazon. Nat. Clim. Change https://doi.org/10.1038/s41558-021-01026-5 (2021).Article 

    Google Scholar 
    Maxwell, S. L. et al. Degradation and forgone removals increase the carbon impact of intact forest loss by 626%. Sci. Adv. 5, eaax2546 (2019).Article 
    CAS 

    Google Scholar 
    Betts, M. G. et al. Global forest loss disproportionately erodes biodiversity in intact landscapes. Nature 547, 441–444 (2017).Article 
    CAS 

    Google Scholar 
    Venter, O. et al. Targeting global protected area expansion for imperiled biodiversity. PLoS Biol. 12, e1001891 (2014).Article 

    Google Scholar 
    Laurance, W. F. et al. Averting biodiversity collapse in tropical forest protected areas. Nature 489, 290–294 (2012).Article 
    CAS 

    Google Scholar 
    Coad, L. et al. Measuring impact of protected area management interventions: current and future use of the global database of protected area management effectiveness. Phil. Trans. R. Soc. B 370, 20140281 (2015).Article 

    Google Scholar 
    Bonan, G. B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).Article 
    CAS 

    Google Scholar 
    Ehbrecht, M. et al. Global patterns and climatic controls of forest structural complexity. Nat. Commun. 12, 519 (2021).Article 
    CAS 

    Google Scholar 
    Zhang, J., Nielsen, S. E., Mao, L., Chen, S. & Svenning, J. C. Regional and historical factors supplement current climate in shaping global forest canopy height. J. Ecol. 104, 469–478 (2016).Article 

    Google Scholar 
    Ellis, E. C. et al. People have shaped most of terrestrial nature for at least 12,000 years. Proc. Natl Acad. Sci. USA 118, e2023483118 (2021).Article 
    CAS 

    Google Scholar 
    Knight, C. A. et al. Land management explains major trends in forest structure and composition over the last millennium in California’s Klamath Mountains. Proc. Natl Acad. Sci. USA 119, e2116264119 (2022).Article 
    CAS 

    Google Scholar 
    Stephens, L. et al. Archaeological assessment reveals Earth’s early transformation through land use. Science 365, 897–902 (2019).Article 
    CAS 

    Google Scholar 
    Asner, G. P., Llactayo, W., Tupayachi, R. & Luna, E. R. Elevated rates of gold mining in the Amazon revealed through high-resolution monitoring. Proc. Natl Acad. Sci. USA 110, 18454–18459 (2013).Article 
    CAS 

    Google Scholar 
    Hoang, N. T. & Kanemoto, K. Mapping the deforestation footprint of nations reveals growing threat to tropical forests. Nat. Ecol. Evol. 5, 845–853 (2021).Article 

    Google Scholar 
    Lim, C. L., Prescott, G. W., De Alban, J. D. T., Ziegler, A. D. & Webb, E. L. Untangling the proximate causes and underlying drivers of deforestation and forest degradation in Myanmar. Conserv. Biol. 31, 1362–1372 (2017).Article 

    Google Scholar 
    Sandel, B. & Svenning, J. C. Human impacts drive a global topographic signature in tree cover. Nat Commun. https://doi.org/10.1038/ncomms3474 (2013).Potapov, P. et al. Mapping the world’s intact forest landscapes by remote sensing. Ecol. Soc. 13, 51 (2008).Article 

    Google Scholar 
    Geldmann, J., Manica, A., Burgess, N. D., Coad, L. & Balmford, A. A global-level assessment of the effectiveness of protected areas at resisting anthropogenic pressures. Proc. Natl Acad. Sci. USA 116, 23209–23215 (2019).Article 
    CAS 

    Google Scholar 
    Yang, H. et al. A global assessment of the impact of individual protected areas on preventing forest loss. Sci. Total Environ. 777, 145995 (2021).Article 
    CAS 

    Google Scholar 
    Jones, K. R. et al. One-third of global protected land is under intense human pressure. Science 360, 788–791 (2018).Article 
    CAS 

    Google Scholar 
    Clerici, N. et al. Deforestation in Colombian protected areas increased during post-conflict periods. Sci. Rep. 10, 4971 (2020).Article 
    CAS 

    Google Scholar 
    Heino, M. et al. Forest loss in protected areas and intact forest landscapes: a global analysis. PLoS ONE 10, e0138918 (2015).Article 

    Google Scholar 
    Leberger, R., Rosa, I. M. D., Guerra, C. A., Wolf, F. & Pereira, H. M. Global patterns of forest loss across IUCN categories of protected areas. Biol. Conserv. 241, 108299 (2020).Article 

    Google Scholar 
    Wade, C. M. et al. What is threatening forests in protected areas? A global assessment of deforestation in protected areas, 2001–2018. Forests 11, 539 (2020).Article 

    Google Scholar 
    Transforming Our World: The 2030 Agenda for Sustainable Development (UN DESA, 2016).Burleson, E. Paris Agreement and consensus to address climate challenge. ASIL Insight 20, 8 (2016).
    Google Scholar 
    Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).Article 
    CAS 

    Google Scholar 
    Quegan, S. et al. The European Space Agency BIOMASS mission: measuring forest above-ground biomass from space. Remote Sens. Environ. 227, 44–60 (2019).Article 

    Google Scholar 
    Simard, M., Pinto, N., Fisher, J. B. & Baccini, A. Mapping forest canopy height globally with spaceborne lidar. J. Geophys. Res. Biogeosci. https://doi.org/10.1029/2011JG001708 (2011).Potapov, P. et al. Mapping global forest canopy height through integration of GEDI and Landsat data. Remote Sens. Environ. 253, 112165 (2021).Article 

    Google Scholar 
    Atkins, J. W., Fahey, R. T., Hardiman, B. S. & Gough, C. M. Forest canopy structural complexity and light absorption relationships at the subcontinental scale. J. Geophys. Res. Biogeosci. 123, 1387–1405 (2018).Article 

    Google Scholar 
    Scarth, P., Armston, J., Lucas, R. & Bunting, P. A structural classification of Australian vegetation using ICESat/GLAS, ALOS PALSAR, and Landsat sensor data. Remote Sens. 11, 147 (2019).Article 

    Google Scholar 
    Dubayah, R. et al. The global ecosystem dynamics investigation: high-resolution laser ranging of the Earth’s forests and topography. Sci. Remote Sens. 1, 100002 (2020).Article 

    Google Scholar 
    Lang, N. et al. Global canopy height regression and uncertainty estimation from GEDI LIDAR waveforms with deep ensembles. Remote Sens. Environ. 268, 112760 (2022).Article 

    Google Scholar 
    Marselis, S. M., Keil, P., Chase, J. M. & Dubayah, R. The use of GEDI canopy structure for explaining variation in tree species richness in natural forests. Environ. Res. Lett. 17, 045003 (2022).Article 

    Google Scholar 
    MacArthur, R. H. & MacArthur, J. W. On bird species diversity. Ecology 42, 594–598 (1961).Article 

    Google Scholar 
    Walter, J. A., Stovall, A. E. L. & Atkins, J. W. Vegetation structural complexity and biodiversity in the Great Smoky Mountains. Ecosphere 12, e03390 (2021).Article 

    Google Scholar 
    Camps-Valls, G. et al. A unified vegetation index for quantifying the terrestrial biosphere. Sci. Adv. 7, eabc7447 (2021).Article 
    CAS 

    Google Scholar 
    Kennedy, C. M., Oakleaf, J. R., Theobald, D. M., Baruch-Mordo, S. & Kiesecker, J. Managing the middle: a shift in conservation priorities based on the global human modification gradient. Glob. Change Biol. 25, 811–826 (2019).Article 

    Google Scholar 
    Weiss, D. J. et al. A global map of travel time to cities to assess inequalities in accessibility in 2015. Nature 553, 333–336 (2018).Article 
    CAS 

    Google Scholar 
    Chazdon, R. L. et al. A policy‐driven knowledge agenda for global forest and landscape restoration. Conserv. Lett. 10, 125–132 (2017).Article 

    Google Scholar 
    Skidmore, A. K. et al. Priority list of biodiversity metrics to observe from space. Nat. Ecol. Evol. 5, 896–906 (2021).Article 

    Google Scholar 
    Schneider, F. D. et al. Mapping functional diversity from remotely sensed morphological and physiological forest traits. Nat. Commun. 8, 1441 (2017).Article 

    Google Scholar 
    Grantham, H. S. et al. Anthropogenic modification of forests means only 40% of remaining forests have high ecosystem integrity. Nat. Commun. 11, 5978 (2020).Article 
    CAS 

    Google Scholar 
    Ponta, N. et al. Drivers of transgression: what pushes people to enter protected areas. Biol. Conserv. 257, 109121 (2021).Article 

    Google Scholar 
    Pack, S. M. et al. Protected area downgrading, downsizing, and degazettement (PADDD) in the Amazon. Biol. Conserv. 197, 32–39 (2016).Article 

    Google Scholar 
    Tollefson, J. Illegal mining in the Amazon hits record high amid Indigenous protests. Nature 598, 15–16 (2021).Article 
    CAS 

    Google Scholar 
    Thies, C., Rosoman, G., Cotter, J. & Meaden, S. Intact Forest Landscapes. Why It Is Crucial to Protect Them from Industrial Exploitation Technical Note Bd 5 (Greenpeace, 2011).Chazdon, R. L. Beyond deforestation: restoring forests and ecosystem services on degraded lands. Science 320, 1458–1460 (2008).Article 
    CAS 

    Google Scholar 
    Lindenmayer, D. B. et al. New policies for old trees: averting a global crisis in a keystone ecological structure. Conserv. Lett. 7, 61–69 (2014).Article 

    Google Scholar 
    Dave, R. et al. Second Bonn Challenge Progress Report: Application of the Barometer in 2018 (IUCN, 2018).Tang, H. & Armston, J. Algorithm Theoretical Basis Document (ATBD) for GEDI L2B Footprint Canopy Cover and Vertical Profile Metrics (Goddard Space Flight Center, 2019).Adam, M., Urbazaev, M., Dubois, C. & Schmullius, C. Accuracy assessment of GEDI terrain elevation and canopy height estimates in European temperate forests: influence of environmental and acquisition parameters. Remote Sens. 12, 3948 (2020).Article 

    Google Scholar 
    Dorado-Roda, I. et al. Assessing the accuracy of GEDI data for canopy height and aboveground biomass estimates in Mediterranean forests. Remote Sens. 13, 2279 (2021).Article 

    Google Scholar 
    Duncanson, L. et al. Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission. Remote Sens. Environ. 270, 112845 (2022).Article 

    Google Scholar 
    Hofton, M., Blair, J. B., Story, S. & Yi, D. Algorithm Theoretical Basis Document (ATBD) (NASA, 2020).Dubayah, R. et al. GEDI L3 Gridded Land Surface Metrics v.2 (ORNL DAAC, 2021).Roy, D. P., Kashongwe, H. B. & Armston, J. The impact of geolocation uncertainty on GEDI tropical forest canopy height estimation and change monitoring. Sci. Remote Sens. 4, 100024 (2021).Article 

    Google Scholar 
    Potapov, P., Hansen, M. C., Stehman, S. V., Loveland, T. R. & Pittman, K. Combining MODIS and Landsat imagery to estimate and map boreal forest cover loss. Remote Sens. Environ. 112, 3708–3719 (2008).Article 

    Google Scholar 
    Dinerstein, E. et al. An ecoregion-based approach to protecting half the terrestrial realm. Bioscience 67, 534–545 (2017).Article 

    Google Scholar 
    Silva, C. A. et al. rGEDI: NASA’s global ecosystem ynamics investigation (GEDI) data visualization and processing. R package version 0.1.2. (2020).The R Project for Statistical Computing (The R Foundation, 2014); https://www.R-project.org/Fischer, B., Smith, M., Pau, G., Morgan, M. & van Twisk, D. rhdf5: R interface to HDF5. R package version 2.40.0 (2022).Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191 (2018).Article 

    Google Scholar 
    Giglio, L., Loboda, T., Roy, D. P., Quayle, B. & Justice, C. O. An active-fire based burned area mapping algorithm for the MODIS sensor. Remote Sens. Environ. 113, 408–420 (2009).Article 

    Google Scholar 
    Hengl, T. & Wheeler, I. Soil organic carbon content in x 5 g/kg at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution. Zenodo https://doi.org/10.5281/zenodo.1475458 (2018).Farr, T. The shuttle radar topography mission. Rev. Geophys. https://doi.org/10.1029/2005RG000183 (2007).James, G., Witten, D., Hastie, T. & Tibshirani, R. An Introduction to Statistical Learning Vol. 112 (Springer, 2013).Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).Article 

    Google Scholar 
    Bivand, R. et al. Package ‘spdep’: spatial dependence: weighting schemes, statistics version 1.2-7 (The Comprehensive R Archive Network, 2015).Bivand, R., Yu, D., Nakaya, T., Garcia-Lopez, M.-A. & Bivand, M. R. Package ‘spgwr’: geographically eighted regression. R package version 0.6-35 (2020).Fotheringham, A. S., Brunsdon, C. & Charlton, M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships (Wiley, 2003). More

  • in

    Revisiting Mt Fuji’s groundwater origins with helium, vanadium and environmental DNA tracers

    Chakraborty, A. & Jones, T. E. in Natural Heritage of Japan Geoheritage, Geoparks and Geotourism (Conservation and Management Series) (eds Chakraborty, A. et al.) Ch. 16 (Springer, 2018).Nakamura, K. Possible nascent trench along the eastern Japan Sea as the convergent boundary between Eurasian and North American plates (in Japanese). Bull. Earthq. Res. Inst. 58, 711–722 (1983).
    Google Scholar 
    Seno, T. Is northern Honshu a microplate? Tectonophysics 115, 177–196 (1985).Article 

    Google Scholar 
    Ogawa, Y., Takami, Y. & Takazawa, S. in Formation and Applications of the Sedimentary Record in Arc Collision Zones Vol. 436 (eds Draut, A. E. at al.) 155–170 (Geological Society of America, 2008).Tsuya, H. & Morimoto, R. Types of volcanic eruptions in Japan (in Japanese). Bull. Volcanol. 26, 209–222 (1963).Article 
    CAS 

    Google Scholar 
    Aoki, Y., Tsunematsu, K. & Yoshimoto, M. Recent progress of geophysical and geological studies of Mt. Fuji Volcano, Japan. Earth Sci. Rev. 194, 264–282 (2019).Article 

    Google Scholar 
    Tsuchi, R. Geology and groundwater of Mt. Fuji, Japan (in Japanese). J. Geogr. 126, 33–42 (2017).Article 

    Google Scholar 
    Vittecoq, B., Reninger, P.-A., Lacquement, F., Martelet, G. & Violette, S. Hydrogeological conceptual model of andesitic watersheds revealed by high-resolution heliborne geophysics. Hydrol. Earth Sys. Sci. 23, 2321–2338 (2019).Article 
    CAS 

    Google Scholar 
    Yamamoto, S. Hydrologic study of volcano Fuji and its adjacent areas (in Japanese). Geogr. Rev. Jpn 43, 267–184 (1970).Article 

    Google Scholar 
    Yamamoto, T. & Nakada, S. in Volcanic Hazards, Risks, and Disasters (eds Shroder, J. F. & Papale, P.) 355–376 (Elsevier, 2015).Hasegawa, A. et al. Plate subduction, and generation of earthquakes and magmas in Japan as inferred from seismic observations: an overview. Gondwana Res. 16, 370–400 (2009).Article 

    Google Scholar 
    Kashiwagi, H. & Nakajima, J. Three‐dimensional seismic attenuation structure of central Japan and deep sources of arc magmatism. Geophys. Res. Lett. 46, 13746–13755 (2019).Article 

    Google Scholar 
    Obrochta, S. P. et al. Mt. Fuji Holocene eruption history reconstructed from proximal lake sediments and high-density radiocarbon dating. Quat. Sci. Rev. 200, 395–405 (2018).Article 

    Google Scholar 
    Tosaki, Y. & Asai, K. Groundwater ages in Mt. Fuji (in Japanese). J. Geogr. 126, 89–104 (2017).Article 

    Google Scholar 
    Imtiaz, M. et al. Vanadium, recent advancements and research prospects: a review. Environ. Int. 80, 79–88 (2015).Article 
    CAS 

    Google Scholar 
    Koshimizu, S., & Tomura, K. (2000). Geochemical behavior of trace vanadium in the spring, groundwater and lake water at the foot of Mt. Fuji, Central Japan. In K. Sato & Y. Iwasa (Eds.), Groundwater Updates. Springer, Tokyo. 171-176. https://doi.org/10.1007/978-4-431-68442-8_29Ono, M. et al. Regional groundwater flow system in a stratovolcano adjacent to a coastal area: a case study of Mt. Fuji and Suruga Bay, Japan. Hydrogeol. J. 27, 717–730 (2019).Article 

    Google Scholar 
    UNESCO Fujisan, Sacred Place and Source of Artistic Inspiration (World Heritage Convention, 2013); https://whc.unesco.org/en/list/1418Nationally Designated Cultural Properties Database (in Japanese) (Agency of Cultural Affairs Japan, 2020); https://kunishitei.bunka.go.jp/bsys/indexShowa’s 100 Famous Waters of Japan (Ministry of the Environment Japan (MOEJ), 1985); https://www.env.go.jp/water/meisui/Heisei’s 100 Famous Waters of Japan (MOEJ, 2009): https://www.env.go.jp/water/meisui/An Overview of the Bottled Water Market in Japan (Frost & Sullivan, 2016).Fujiyoshida Mineral Water Conservation Association FMWCA Regulations (in Japanese) (Mt. Fuji Springs Inc., 2016); http://fujiyoshida-hozen.org/aboutwater/Adachi, Y. et al. The physiological effects of the undercurrent water from Mt. Fuji on type 2 diabetic KK-Ay mice. Biomed. Res. Trace Elem. 15, 76–78 (2004).CAS 

    Google Scholar 
    Isogai, A., Kanada, R., Iawata, H. & Sudo, S. The influence of vanadium on the components of hineka (in Japanese). J. Brew. Soc. Jpn 107, 443–450 (2012).Article 

    Google Scholar 
    Tamada, Y., Tokui, M., Yamashita, N., Kubodera, T. & Akashi, T. Analyzing the relationship between the inorganic element profile of sake dilution water and dimethyl trisulfide formation using multi-element profiling. J. Biosci. Bioeng. 127, 710–713 (2019).Article 
    CAS 

    Google Scholar 
    London Sake Challenge 2018: Awarded Sake (Sake Somelier Association (SSA), 2018); https://londonsakechallenge.com/awarded-sake-2019/London Sake Challenge 2019: Awarded Sake (SSA, 2019); https://londonsakechallenge.com/awarded-sake-2019/Yasuhara, M., Hayashi, T. & Asai, K. Overview of the special issue “Groundwater in Mt. Fuji”. J. Geogr. 126, 25–27 (2017).Article 

    Google Scholar 
    Yasuhara, M., Hayashi, T., Asai, K., Uchiyama, M. & Nakamura, T. Overview of the special issue “Groundwater in Mt. Fuji (Part 2)”. J. Geogr. 129, 657–660 (2020).Article 

    Google Scholar 
    Gmati, S., Tase, N., Tsujimura, M. & Tosaki, Y. Aquifers interaction in the southwestern foot of Mt. Fuji, Japan, examined through hydrochemistry and statistical analyses. Hydrol. Res. Lett. 5, 58–63 (2011).Article 

    Google Scholar 
    Ikeda, K. Water-sediments interaction of salinized groundwater, and its chemical compositions in coastal areas (in Japanese). Jpn. J. Limnol. 46, 303–314 (1985).Article 
    CAS 

    Google Scholar 
    Kato, K. et al. Unveiled groundwater flushing from the deep seafloor in Suruga Bay. Limnology https://doi.org/10.1007/s10201-014-0445-0 (2015).Segawa, T. et al. Microbes in groundwater of a volcanic mountain, Mt. Fuji; 16S rDNA phylogenetic analysis as a possible indicator for the transport routes of groundwater. Geomicrobiol. J. 32, 677–688 (2015).Article 

    Google Scholar 
    Sugiyama, A., Masuda, S., Nagaosa, K., Tsujimura, M. & Kato, K. Tracking the direct impact of rainfall on groundwater at Mt. Fuji by multiple analyses including microbial DNA. Biogeosciences 15, 721–732 (2018).Article 
    CAS 

    Google Scholar 
    Yasuhara, M., Kazahaya, K. & Marui, A. in Fuji Volcano (eds Aramaki, S. et al.) 389–405 (Yamanashi Institute of Environmental Sciences, 2007).Tsuchi, R. in Fuji Volcano (eds Aramaki, S. et al.) 375–387 (Yamanashi Institute of Environmental Sciences, 2007).Takada, A., Yamamoto, T., Ishizuka, Y. & Nakano, S. in Miscellaneous Map Series No. 12, 56 (Geological Survey of Japan (GSJ), National Institute of Advanced Industrial Science and Technology (AIST), 2016).Uchiyama, T. Hydrogeological structure and hydrological characterization in the northern foot area of Fuji volcano, central Japan (in Japanese). J. Geogr. 129, 697–724 (2020).Article 

    Google Scholar 
    Ikawa, R. et al. in S-5: Seamless Geoinformation of Coastal Zone “Northern Coastal Zone of Suruga Bay” (GSJ, AIST, 2016).AIST 2014 Marine Geological and Environmental Survey Confirmation Technology Development Results Report (in Japanese) (AIST, 2015).AIST 2015 Marine Geological and Environmental Survey Confirmation Technology Development Results Report (in Japanese) (AIST, 2016).Lin, A., Iida, K. & Tanaka, H. On-land active thrust faults of the Nankai–Suruga subduction zone: the Fujikawa-kako Fault Zone, central Japan. Tectonophysics 601, 1–19 (2013).Article 

    Google Scholar 
    Fujita, E. et al. Stress field change around the Mount Fuji volcano magma system caused by the Tohoku megathrust earthquake, Japan. Bull. Volcanol. 75, 679 (2013).Article 

    Google Scholar 
    Kano, K.-I., Odawara, K., Yamamoto, G. & Ito, T. Tectonics of the Fujikawa-kako Fault Zone around the Hoshiyama Hills, central Japan, since 1 Ma. Geosci. Rep. Shizuoka Univ. 46, 19–49 (2019).
    Google Scholar 
    Schilling, O. S., Cook, P. G. & Brunner, P. Beyond classical observations in hydrogeology: the advantages of including exchange flux, temperature, tracer concentration, residence time and soil moisture observations in groundwater model calibration. Rev. Geophys. 57, 146–182 (2019).Article 

    Google Scholar 
    Schilling, O. S. et al. Quantifying groundwater recharge dynamics and unsaturated zone processes in snow-dominated catchments via on-site dissolved gas analysis. Water Resour. Res. 57, e2020WR028479 (2021).Article 

    Google Scholar 
    National Hydrological Environment Database of Japan (GSJ, AIST, 2020).Hayashi, T. Understanding the groundwater flow system at the northern part of Mt. Fuji: current issues and prospects (in Japanese). J. Geogr. 129, 677–695 (2020).Article 

    Google Scholar 
    Yasuhara, M., Marui, A., & Kazahaya, K. (1997). Stable isotopic composition of groundwater from Mt. Yatsugatake and Mt. Fuji, Japan. Proceedings of the Rabat Symposium. Rabat Symposium, April 1997, Wallingford, UK.Jasechko, S. Global isotope hydrogeology—review. Rev. Geophys. https://doi.org/10.1029/2018RG000627 (2019).Yaguchi, M., Muramatsu, Y., Chiba, H., Okumura, F. & Ohba, T. The origin and hydrochemistry of deep well waters from the northern foot of Mt. Fuji, central Japan. Geochem. J. 50, 227–239 (2016).Article 
    CAS 

    Google Scholar 
    Aizawa, K. et al. Gas pathways and remotely triggered earthquakes beneath Mount Fuji, Japan. Geology 44, 127–130 (2016).Article 
    CAS 

    Google Scholar 
    Kipfer, R. et al. Injection of mantle type helium into Lake Van (Turkey): the clue for quantifying deep water renewal. Earth Planet. Sci. Lett. 125, 357–370 (1994).Article 
    CAS 

    Google Scholar 
    Kipfer, R., Aeschbach-Hertig, W., Peeters, F. & Stute, M. in Noble Gases in Geochemistry and Cosmochemistry Reviews in Mineralogy and Geochemistry Vol. 47 (eds Porcelli, D. et al.) Ch. 14 (De Gruyter, 2002).Sano, Y. & Fischer, T. P. in The Noble Gases as Geochemical Tracers: Advances in isotope geochemistry (ed. Burnard, O.) Ch. 10 (Springer, 2013).Sano, Y. & Wakita, H. Distribution of 3He/4He ratios and its implications for geotectonic structure of the Japanese Islands. J. Geophys. Res. 90, 8729–8741 (1985).Article 
    CAS 

    Google Scholar 
    Tomonaga, Y. et al. Fluid dynamics along the Nankai Trough: He isotopes reveal direct seafloor mantle-fluid emission in the Kumano Basin (Southwest Japan). ACS Earth Space Chem. 4, 2015–2112 (2020).Article 

    Google Scholar 
    Chen, A. et al. Mantle fluids associated with crustal-scale faulting in a continental subduction setting, Taiwan. Sci Rep. 9, 10805 (2019).Article 

    Google Scholar 
    Crossey, L. J. et al. Continental smokers couple mantle degassing and distinctive microbiology within continents. Earth Planet. Sci. Lett. 435, 22–30 (2016).Article 
    CAS 

    Google Scholar 
    Crossey, L. J. et al. Degassing of mantle-derived CO2 and He from springs in the southern Colorado Plateau region—neotectonic connections and implications for groundwater systems. Geol. Soc. Am. Bull. 121, 1034–1053 (2009).Article 
    CAS 

    Google Scholar 
    Kusuda, C., Iwamori, H., Nakamura, H., Kazahaya, K. & Morikawa, N. Arima hot spring waters as a deep-seated brine from subducting slab. Earth Planets Space 66, 119 (2014).Article 

    Google Scholar 
    Sano, Y., Kameda, A., Takahata, N., Yamamoto, J. & Nakajima, J. Tracing extinct spreading center in SW Japan by helium-3 emanation. Chem. Geol. 266, 50–56 (2009).Article 
    CAS 

    Google Scholar 
    Sano, Y. et al. Groundwater helium anomaly reflects strain change during the 2016 Kumamoto earthquake in Southwest Japan. Sci. Rep. 6, 37939 (2016).Article 
    CAS 

    Google Scholar 
    Peeters, F. et al. Improving noble gas based paleoclimate reconstruction and groundwater dating using 20Ne/22Ne ratios. Geochim. Cosmochim. Acta 67, 587–600 (2002).Article 

    Google Scholar 
    Reimann, C. & de Caritat, P. Chemical Elements in the Environment 398 (Springer, 1998).Hamada, T. in Vanadium in the Environment. Part 1: Chemistry and Biochemistry Advances in Environmental Sciences and Technology Vol. 10 (ed. Nriagu, J. O.) 97–123 (Wiley & Sons, 1998).Koshimizu, S. & Kyotani, T. Geochemical behaviors of multi-elements in water samples from the Fuji and Sagami Rivers, Central Japan, using vanadium as an effective indicator. Jpn J. Limnol. 63, 113–124 (2002).Article 
    CAS 

    Google Scholar 
    Sohrin, R. in Green Science and Technology (eds Park, E. Y. et al.) Ch. 7 (CRC, 2019).Wehrli, B. & Stumm, W. Oxygenation of vanadyl(IV). Effect of coordinated surface hydroxyl groups and hydroxide ion. Langmuir 4, 753–758 (1988).Article 
    CAS 

    Google Scholar 
    Wright, M. T. & Belitz, K. Factors controlling the regional distribution of vanadium in groundwater. Ground Water 48, 515–525 (2010).Article 
    CAS 

    Google Scholar 
    Deverel, S. J., Goldberg, S. & Fujii, R. in Agricultural salinity assessment and management (eds W.W. Wallender & K.K. Tanji) 89–137 (American Society of Civil Engineers, 2012).Wehrli, B. & Stumm, W. Vanadyl in natural waters: adsorption and hydrolysis promote oxygenation. Geochim. Cosmochim. Acta 53, 69–77 (1989).Article 
    CAS 

    Google Scholar 
    Chen, G. & Liu, H. Understanding the reduction kinetics of aqueous vanadium(V) and transformation products using rotating ring-disk electrodes. Environ. Sci. Technol. 51, 11643–11651 (2017).Article 
    CAS 

    Google Scholar 
    Telfeyan, K., Johannesson, K. H., Mohajerin, T. J. & Palmore, C. D. Vanadium geochemistry along groundwater flow paths in contrasting aquifers of the United States: Carrizo Sand (Texas) and Oasis Valley (Nevada) aquifers. Chem. Geol. 410, 63–78 (2015).Article 
    CAS 

    Google Scholar 
    Kan, K. et al. Archaea in Yellowstone Lake. ISME J. 5, 1784–1795 (2011).Article 
    CAS 

    Google Scholar 
    Wong, H. L. et al. Dynamics of archaea at fine spatial scales in Shark Bay mat microbiomes. Sci. Rep. 7, 46160 (2017).Article 
    CAS 

    Google Scholar 
    Ikeda, K. A study on chemical characteristics of ground water in Fuji area (in Japanese). J. Groundw. Hydrol. 24, 77–93 (1982).
    Google Scholar 
    Aizawa, K. et al. Hydrothermal system beneath Mt. Fuji volcano inferred from magnetotellurics and electric self-potential. Earth Planet. Sci. Lett. 235, 343–355 (2005).Article 
    CAS 

    Google Scholar 
    Yamamoto, T., Takada, A., Ishizuka, Y., Miyaji, N. & Tajima, Y. Basaltic pyroclastic flows of Fuji volcano, Japan: characteristics of the deposits and their origin. Bull. Volcanol. 67, 622–633 (2005).Article 

    Google Scholar 
    Yamamoto, T., Takada, A., Ishizuka, Y. & Nakano, S. Chronology of the products of Fuji volcano based on new radiometoric carbon ages (in Japanese). Bull. Volcanol. 50, 53–70 (2005).CAS 

    Google Scholar 
    Aizawa, K., Yoshimura, R. & Oshiman, N. Splitting of the Philippine Sea Plate and a magma chamber beneath Mt. Fuji. Geophys. Res. Lett. 31, L09603 (2004).Article 

    Google Scholar 
    Nakamura, H., Iwamori, H. & Kimura, J.-I. Geochemical evidence for enhanced fluid flux due to overlapping subducting plates. Nat. Geosci. 1, 380–384 (2008).Article 
    CAS 

    Google Scholar 
    Kaneko, T., Yasuda, A., Fujii, T. & Yoshimoto, M. Crypto-magma chambers beneath Mt. Fuji. J. Volcanol. Geotherm. Res. 193, 161–170 (2010).Article 
    CAS 

    Google Scholar 
    Tsuya, H., Machida, H., & Shimozuru, D. (1988). Geology of volcano Mt. Fuji. Explanatory text of the geologic map of Mt. Fuji (scale 1:50,000; second printing). Geological Survey of Japan (GSJ), Tsukuba, Japan.Yoshimoto, M. et al. Evolution of Mount Fuji, Japan: inference from drilling into the subaerial oldest volcano, pre-Komitake. Isl. Arc. 19, 470–488 (2010).Article 

    Google Scholar 
    Shikazono, N., Arakawa, T. & Nakano, T. Groundwater quality, flow, and nitrogen pollution at the southern foot of Mt. Fuji (in Japanese). J. Geogr. 123, 323–342 (2014).Article 

    Google Scholar 
    Tosaki, Y., Tase, N., Sasa, K., Takahashi, T. & Nagashima, Y. Estimation of groundwater residence time using the 36Cl bomb pulse. Groundwater 49, 891–902 (2011).Article 
    CAS 

    Google Scholar 
    Yamamoto, T. Geology of the Southwestern Part of Fuji Volcano (in Japanese) 27 (GSJ, AIST, 2014).Tsuya, H. Geology of volcano Mt. Fuji. Explanatory text of the geologic map of Mt. Fuji (scale 1:50,000). Geological Survey of Japan, Tsukuba, Japan. (1968).Tomiyama, S., Ii, H., Miyaike, S., Hattori, R. & Ito, Y. Estimation of the sources and flow system of groundwater in Fuji-Gotenba area by stable isotopic analysis and groundwater flow simulation (in Japanese). Bunseki Kagaku 58, 865–872 (2009).Article 
    CAS 

    Google Scholar 
    Oguchi, T. & Oguchi, C. T. in Geomorphological Landscapes of the World (ed. Migoń, P.) Ch. 31 (Springer, 2010).Mean Annual Precipitation from 1981-2010 Recorded at the Four Mt. Fuji Observatories (Mishima, Fuji, Furuseki, Yamanaka) (Japan Meteorological Agency, 2015).Schilling, O. S., Park, Y.-J., Therrien, R. & Nagare, R. M. Integrated surface and subsurface hydrological modeling with snowmelt and pore water freeze-thaw. Groundwater 57, 63–74 (2018).Article 

    Google Scholar 
    Sakio, H. & Masuzawa, T. Advancing timberline on Mt. Fuji between 1978 and 2018. Plants 9, 1537 (2020).Article 

    Google Scholar 
    Asai, K. & Koshimizu, S. 3H/3He-based groundwater ages for springs located at the foot of Mt. Fuji (in Japanese). J. Groundw. Hydrol. 61, 291–298 (2019).Article 

    Google Scholar 
    Sakai, Y., Shita, K., Koshimizu, S. & Tomura, K. Geochemical study of trace vanadium in water by preconcentrational neutron activation analysis. J. Radioanal. Nucl. Chem. 216, 203–212 (1997).Article 
    CAS 

    Google Scholar 
    Nahar, S. & Zhang, J. Concentration and distribution of organic and inorganic water pollutants in eastern Shizuoka, Japan. Toxicol. Environ. Chem. https://doi.org/10.1080/02772248.2011.610498 (2011).Kamitani, T., Watanabe, M., Muranaka, Y., Shin, K.-C. & Nakano, T. Geographical characteristics and sources of dissolved ions in groundwater at the southern part of Mt. Fuji (in Japanese). J. Geogr. 126, 43–71 (2017).Article 

    Google Scholar 
    Kawagucci, S. et al. Disturbance of deep-sea environments induced by the M9.0 Tohoku earthquake. Sci Rep. 2, 270 (2012).Article 

    Google Scholar 
    Uchida, N. & Bürgmann, R. A decade of lessons learned from the 2011 Tohoku-Oki earthquake. Rev. Geophys. 59, e2020RG000713 (2021).Article 

    Google Scholar 
    Mahara, Y., Igarashi, T. & Tanaka, Y. Groundwater ages of confined aquifer in Mishima lava flow, Shizuoka (in Japanese). J. Groundw. Hydrol. 35, 201–215 (1993).Article 

    Google Scholar 
    Nakamura, T. et al. Sources of water and nitrate in springs at the northern foot of Mt. Fuji and nitrate loading in the Katsuragawa River (in Japanese). J. Geogr. 126, 73–88 (2017).Article 

    Google Scholar 
    Notsu, K., Mori, T., Sumino, H. & Ohno, M. in Fuji Volcano (eds Aramaki, S. et al.) 173–182 (Yamanashi Institute of Environmental Sciences, 2007).Ogata, M. & Kobayashi, H. Hydrologic Science Research for the Management and Utilization of Ground Water Resources in the Northern Piedmont Area of Mt. Fuji: Fluorine Ion and Vanadium Contained in Ground Water at the Northern Foot of Mt. Fuji (Yamanashi Industrial Technology Center, 2015).Ogata, M., Kobayashi, H. & Koshimizu, S. Concentration of fluorine in groundwater and groundwater table at the northern foot of Mt. Fuji (in Japanese). J. Groundw. Hydrol. 56, 35–51 (2014).Article 

    Google Scholar 
    Ohno, M., Sumino, H., Hernandez, P. A., Sato, T. & Nagao, K. Helium isotopes in the Izu Peninsula, Japan: relation of magma and crustal activity. J. Volcanol. Geotherm. Res. 199, 118–126 (2011).Article 
    CAS 

    Google Scholar 
    Okabe, S., Shibasaki, M., Oikawa, T., Kawaguchi, Y. & Nihongi, H. Geochemical studies of spring and lake waters on and around Mt. Fuji (in Japanese). J. Sch. Mar. Sci. Technol. Tokai Univ. 14, 81–105 (1981).CAS 

    Google Scholar 
    Ono, M., Ikawa, R., Machida, H. & Marui, A. Distribution of radon concentration in groundwater at the southwestern foot of Mt. Fuji (in Japanese). Radioisotopes 65, 431–439 (2016).Article 
    CAS 

    Google Scholar 
    Tosaki, Y. Estimation of Groundwater Residence Time Using Bomb-Produced Chlorine-36. PhD thesis, Univ. Tsukuba (2008).Umeda, K., Asamori, K. & Kusano, T. Release of mantle and crustal helium from a fault following an inland earthquake. Appl. Geochem. 37, 134–141 (2013).Article 
    CAS 

    Google Scholar 
    Yamamoto, C. Estimation of Groundwater Flow System Using Multi-tracer Techniques in Mt. Fuji, Japan. (in Japanese) PhD thesis, Univ. Tsukuba (2016).Yamamoto, S. & Nakamura, T. Visit to valuable water springs (129) valuable water at the northern foot of Mount Fuji (Fuji-Kawaguchiko Town) (in Japanese). J. Groundw. Hydrol. 62, 329–336 (2020).Article 

    Google Scholar 
    Yamamoto, S. et al. Water sources of lake bottom springs in Lake Kawaguchi, northern foot of Mount Fuji, Japan (in Japanese). J. Geogr. 129, 665–676 (2020).Article 

    Google Scholar 
    Yamamoto, S., Nakamura, T. & Uchiyama, T. Newly discovered lake bottom springs from Lake Kawaguchi, the northern foot of Mount Fuji, Japan (in Japanese). J. Jpn Assoc. Hydrol. Sci. 47, 49–59 (2017).
    Google Scholar 
    Yamamoto, S., Nakamura, T., Koishikawa, H. & Uchiyama, T. Water quality of shallow groundwater in the southern coast area of Lake Kawaguchi at the northern foot of Mt. Fuji, Yamanashi, Japan (in Japanese). Mt Fuji Res. 11, 1–9 (2017).
    Google Scholar 
    Coplen, T. B. Reporting of stable hydrogen, carbon, and oxygen isotopic abundances. Geothermics 66, 273–276 (1994).CAS 

    Google Scholar 
    Nimz, G. J. in Isotope Tracers in Catchment Hydrology (eds Kendall, C. & McDonnell, J. J.) Ch. 8 (Elsevier, 1998).Bullen, T. D. & Kendall, C. in Isotope Tracers in Catchment Hydrology (eds Kendall, C. & McDonnell, J. J.) Ch. 18 (Elsevier, 1998).Vanadium Pentoxide and Other Inorganic Vanadium Compounds Vol. 29 (WHO, 2001).Nagai, T., Takahashi, M., Hirahara, Y. & Shuto, K. Sr-Nd isotopic compositions of volcanic rocks from Fuji, Komitake and Ashitaka Volcanoes, Central Japan (in Japanese). Proc. Inst. Nat. Sci. Nihon Univ. 39, 205–215 (2004).CAS 

    Google Scholar 
    Hogan, J. F. & Blum, J. D. Tracing hydrologic flow paths in a small forested watershed using variations in 87Sr/86Sr, [Ca]/[Sr], [Ba]/[Sr] and δ18O. Water Resour. Res. 39, 1282 (2003).Article 

    Google Scholar 
    Koshikawa, M. K. et al. Using isotopes to determine the contribution of volcanic ash to Sr and Ca in stream waters and plants in a granite watershed, Mt. Tsukuba, central Japan. Environ. Earth Sci. 75, 501 (2016).Article 

    Google Scholar 
    Graustein, W. C. in Stable Isotopes in Ecological Research Ecological Studies (Analysis and Synthesis) (eds Rundel, JP.W. et al.) Ch. 28 (Springer, 1989).Cook, P. G. & Böhlke, J.-K. in Environmental Tracers in Subsurface Hydrology (eds Cook, P. G. & Herczeg, A. L.) Ch. 1 (Springer, 2000).Aeschbach-Hertig, W. & Solomon, D. K. in The Noble Gases as Geochemical Tracers (ed. Burnard, P.) Ch. 5 (Springer, 2013).Popp, A. L. et al. A framework for untangling transient groundwater mixing and travel times. Water Resour. Res. 57, e2020WR028362 (2021).Article 

    Google Scholar 
    Schilling, O. S. et al. Advancing physically-based flow simulations of alluvial systems through observations of 222Rn, 3H/3He, atmospheric noble gases and the novel 37Ar tracer method. Water Resour. Res. 53, 10465–10490 (2017).Article 

    Google Scholar 
    Tomonaga, Y. et al. Using noble-gas and stable-isotope data to determine groundwater origin and flow regimes: applicatoin to the Ceneri Base Tunnel (Switzerland). J. Hydrol. 545, 395–409 (2017).Article 
    CAS 

    Google Scholar 
    Niu, Y. et al. Noble gas signatures in the island of Maui, Hawaii – characterizing groundwater sources in fractured systems. Water Resour. Res. 53, 3599–3614 (2017).Article 

    Google Scholar 
    Warrier, R. B., Castro, M. C. & Hall, C. M. Recharge and source-water insights from the Galapagos Islands using noble gases and stable isotopes. Water Resour. Res. https://doi.org/10.1029/2011WR010954 (2012).Schilling, O. S. et al. Buried paleo-channel detection with a groundwater model, tracer-based observations, and spatially varying, preferred anisotropy pilot point calibration. Geophys. Res. Lett. 49, e2022GL098944 (2022).Article 

    Google Scholar 
    Brennwald, M. S., Schmidt, M., Oser, J. & Kipfer, R. A portable and autonomous mass spectrometric system for on-site environmental gas analysis. Environ. Sci. Technol. 50, 13455–12463 (2016).Article 
    CAS 

    Google Scholar 
    Tomonaga, Y. et al. On-line monitoring of the gas composition in the full-scale emplacement experiment at Mont Terri (Switzerland). Appl. Geochem. 100, 234–243 (2019).Article 
    CAS 

    Google Scholar 
    Brennwald, M. S., Tomonaga, Y. & Kipfer, R. Deconvolution and compensation of mass spectrometric overlap interferences with the miniRUEDI portable mass spectrometer. MethodsX 7, 101038 (2020).Article 
    CAS 

    Google Scholar 
    Van Rossum, G. & Drake, F. L. Python 3 Reference Manual (CreateSpace, 2009).Beyerle, U. et al. A mass spectrometric system for the analysis of noble gases and tritium from water samples. Environ. Sci. Technol. 34, 2042–2050 (2000).Article 
    CAS 

    Google Scholar 
    Clarke, W. B., Jenkins, W. J. & Top, Z. Determination of tritium by mass spectrometric measurement of 3He. Int. J. Appl. Radiat. Isotopes 27, 515–522 (1976).Article 
    CAS 

    Google Scholar 
    Bucci, A., Petrella, E., Celivo, F. & Naclerio, G. Use of molecular approaches in hydrogeological studies: the case of carbonate aquifers in southern Italy. Hydrogeol. J. 25, 1017–1031 (2017).Article 
    CAS 

    Google Scholar 
    Proctor, C. R. et al. Phylogenetic clustering of small low nucleic acid-content bacteria across diverse freshwater ecosystems. ISME J. 12, 1344–1359 (2018).Article 
    CAS 

    Google Scholar 
    Pronk, M., Goldscheider, N. & Zopfi, J. Microbial communities in karst groundwater and their potential use for biomonitoring. Hydrogeol. J. 17, 37–48 (2009).Article 

    Google Scholar 
    Miller, J. B., Frisbee, M. D., Hamilton, T. L. & Murugapiran, S. K. Recharge from glacial meltwater is critical for alpine springs and their microbiomes. Environ. Res. Lett. 16, 064012 (2021).Article 
    CAS 

    Google Scholar 
    Ginn, T. R. et al. in Encyclopedia of Hydrological Sciences (ed. Anderson, M.G.) Ch. 105 (John Wiley & Sons, 2005).Tufenkji, N. & Emelko, M. B. in Encyclopedia of Environmental Health (ed. Nriagu, J.O.) Vol. 2, 715–726 (Elsevier, 2011).Nevecherya, I. K., Shestakov, V. M., Mazaev, V. T. & Shlepnina, T. G. Survival rate of pathogenic bacteria and viruses in groundwater. Water Res. 32, 209–214 (2005).Article 
    CAS 

    Google Scholar 
    Franzosa, E. A. et al. Sequencing and beyond: integrating molecular ‘omics’ for microbial community profiling. Nature Rev. Microbiol. 13, 360–372 (2015).Article 
    CAS 

    Google Scholar 
    Kimura, H., Ishibashi, J. I., Masuda, H., Kato, K. & Hanada, S. Selective phylogenetic analysis targeting 16S rRNA genes of hyperthermophilic archaea in the deep-subsurface hot biosphere. Appl. Environ. Microbiol. 73, 2110–2117 (2007).Article 
    CAS 

    Google Scholar 
    Somerville, C. C., Knight, I. T., Straube, W. L. & Colwell, R. R. Simple, rapid method for direct isolation of nucleic-acids from aquatic environments. Appl. Environ. Microbiol. 55, 548–554 (1989).Article 
    CAS 

    Google Scholar 
    Takahashi, S., Tomita, J., Nishioka, K., Hisada, T. & Nishijima, M. Development of a prokaryotic universal primer for simultaneous analysis of bacteria and archaea using next-generation sequencing. PLoS ONE https://doi.org/10.1371/journal.pone.0105592 (2014).Wasimuddin et al. Evaluation of primer pairs for microbiome profiling from soils to humans within the One Health framework. Mol. Ecol. Resour. 20, 1558–1571 (2020).Article 
    CAS 

    Google Scholar 
    Suzuki, Y., Shimizu, H., Kuroda, T., Takada, Y. & Nukazawa, K. Plant debris are hotbeds for pathogenic bacteria on recreational sandy beaches. Sci Rep. 11, 11496 (2021).Article 
    CAS 

    Google Scholar 
    Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).Article 
    CAS 

    Google Scholar 
    Caporaso, J. G. et al. QIIME allows analysis of high- throughput community sequencing data. Nat. Methods 7, 335–336 (2010).Article 
    CAS 

    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).Article 
    CAS 

    Google Scholar 
    DeSantis, T. Z. et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072 (2006).Article 
    CAS 

    Google Scholar 
    McMurdie, P. J. & Holmes, S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217 (2013).Article 
    CAS 

    Google Scholar 
    R: A Language and Environment for Statistical Computing v.3.6.2 (R Foundation for Statistical Computing, 2019).Porter, K. G. & Feig, Y. S. The use of DAPI for identifying and counting aquatic microflora. Limnol. Oceanogr. 25, 943–948 (1980).Article 

    Google Scholar 
    Schilling, O. S. et al. Mt. Fuji hydrogeochemical and microbiological dataset. HydroShare https://doi.org/10.4211/hs.4eac370d12e142b5aa718e5deb57da39 (2022).Gotelli, N. J. & Chao, A. in Encyclopedia of Biodiversity Vol. 5 (ed. Levin, S. A.) 195–211 (Academic, 2013).World Imagery (Esri, 2021); https://www.arcgis.com/home/item.html?id=10df2279f9684e4a9f6a7f08febac2a9Elevation Tile Map of Japan (DEM5A; Resolution: 5m) (Geospatial Information Authority of Japan (GSI), 2021).Chiba, T., Kaneta, S. & Suzuki, Y. in The International Archives of the Photogrammetry Vol. XXXVII Ch. B2 (Remote Sensing and Spatial Information Sciences, 2008).Air Asia Survey Co. Ltd Red Relief Image Map of Japan (RRIM 10_2016) (GSI, 2016).Active Fault Database of Japan April 26 2019 edn Disclosure database DB095 (AIST, 2019).Bird, P. An updated digital model of plate boundaries. Geochem. Geophys. Geosyst. https://doi.org/10.1029/2001GC000252 (2003).Van Horne, A., Sato, H. & Ishiyama, T. Evolution of the Sea of Japan back-arc and some unsolved issues. Tectonophysics 710–711, 6–20 (2017).Article 

    Google Scholar 
    Shannon, C. E. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).Article 

    Google Scholar 
    2019 Coastal Disposal System Evaluation Confirmation Technology Results Report (in Japanese) (AIST, 2019). More

  • in

    Large variations in afforestation-related climate cooling and warming effects across short distances

    Study site descriptionThe study was carried out at the edge of an arid region in mature plantations dominated by P. halepensis, of age of 40-50 years (Pinus halepensis), and their adjacent non-forested ecosystems. The sites were distributed across a climatic gradient from arid and semi-arid to dry sub-humid (Fig. 1, Supplementary Table S1). Three selected paired sites included: (1) An arid site at Yatir forest (annual precipitation, P: 280 mm; aridity index, AI: 0.18; elevation: 650 m; light brown Rendzina soil, and forest density: 300 trees ha−1), where a permanent flux tower has been operating since the year 2000 (http://fluxnet.ornl.gov). Note that an AI of 0.2 marks the boundary between arid and semi-arid regions. Yatir, with AI = 0.18, is formally within an arid zone, but on the edge of a semi-arid zone. (2) An intermediate semi-arid site in Eshtaol forest (P: 480 mm; AI: 0.37; elevation: 350 m; light brown Rendzina soil, and forest density: 450 trees ha−1). (3) A dry-subhumid site in northern Israel at the Birya forest (P: 770 mm; AI: 0.64; elevation: 755 m; dark brown Terra-Rossa and Rendzina soil, and forest density: 600 trees ha−1). Non-forest ecosystems were sparse dwarf shrublands, dominated by Sarcopoterium spinosum in a patchy distribution with a wide variety of herbaceous species, mostly annuals, that grew in between the shrubs during winter to early spring, and then dried out. All non-forested sites had been subjected to livestock grazing (exposing soils). Finally, an additional site that was characterized as Oak-forest vegetation was added. The site was dominated by two oak species, Quercus calliprinos and Quercus ithaburensis (P: 540 mm; AI: 0.4; see previous publication for more details on the oak site53). All sites were under high solar radiation load, with an annual average of approximately 240 Wm−2 in the arid region and only 3% lower in the northern site in the dry-subhumid region (Table 1).Mobile laboratoryMeasurements were conducted on a campaign basis using a mobile lab with a flux measurement system at all sites except the Yatir forest, where the permanent flux tower was used (http://fluxnet.ornl.gov;54). Repeated campaigns of approximately two weeks at each site, along the seasonal cycle, were undertaken during 4 years of measurements, 2012–2015 (a total of 6-7 campaigns per site, evenly distributed between the seasons) the 4 years of measurements were found to be representative of previous 70 years of precipitation record (Supplementary Figure S1 and Table S2). The mobile lab was housed on a 12-ton 4 × 4 truck with a pneumatic mast with an eddy-covariance system and provided the facility for any auxiliary and related measurements. Non-radiative flux measurements were undertaken using an eddy-covariance system to quantify CO2, sensible heat (H), and latent heat (LE) fluxes using a 3D sonic anemometer (R3, Gill Instruments, Hampshire, UK) and an enclosed-path CO2/H2O infrared gas analyzer (IRGA; LI-7200, LI-COR). Non-radiative flux measurements were accompanied by meteorological sensors, including air temperature (Ta), relative humidity (RH), and pressure (Campbell Scientific Inc., Logan, UT, USA), radiation fluxes of net solar- and net long-wave radiation (SWnet and LWnet, respectively), and photosynthetic radiation sensors (Kipp & Zonen, Delft, Holland). Raw EC data and the data from the meteorological sensors were collected using a computer and a CR3000 logger (Campbell Sci., Logan, UT, USA), respectively. The EC system was positioned at the center of each field site with the location and height aimed at providing sufficient ‘fetch’ of relatively homogeneous terrains. For detailed information on the use of the mobile lab and the following data processing of short and long-term fluxes see previous publications29,55,56.Data processingMean 30-min fluxes (CO2, LE, and H) were computed using Eddy-pro 5.1.1 software (LiCor, Lincoln, Nebraska, USA). Quality control of the data included a spike removal procedure. A linear fit was used for filling short gaps (below three hours) of missing values due to technical failure. Information about background meteorological parameters, including P, Ta, RH, and global radiation (Rg), was collected from meteorological stations (standard met stations maintained by the Israel Meteorological Service, https://ims.gov.il/en). The data were obtained at half-hourly time resolution and for a continuous period of 15 years since 2000.Estimating continuous fluxes using the flux meteorological algorithmEstimation of the flux-based annual carbon and radiation budgets was undertaken using the short campaign measurements as a basis to produce a continuous, seasonal, annual, and inter-annual scale dataset of ecosystem fluxes (flux meteorological algorithm). The flux meteorological algorithm method was undertaken based on the relationships between measured fluxes (CO2, LE, H, SWnet, and LWnet) and meteorological parameters (Ta, RH, Rg, VPD, and transpiration deficit, a parameter that correlated well with soil moisture, see main text and supplementary material of previous publication29. A two-step multiple stepwise regression was established, first between the measured fluxes (H, LE, and the ecosystem net carbon exchange) and the meteorological parameters measured by the mobile lab devices, and then between the two meteorological datasets (i.e., the variables measured by the Israel meteorological stations) for the same measurement times. Annual fluxes were calculated for the combined dataset of all campaigns at each site using the following generic linear equation:$$y={{{{{rm{a}}}}}}+{Sigma }_{i}{b}_{i}{x}_{i}$$
    (1)
    where, y is the ecosystem flux of interest, the daily average for radiative fluxes (LWnet and SWnet), non-radiative fluxes (H and LE), and daily sum for net ecosystem exchange (NEE), a and ({b}_{i}) are parameters, and ({x}_{i}) is Ta, RH, Rg, vapor pressure deficit, or transpiration deficit. The meteorological variables (({x}_{i})) were selected by stepwise regression, with ({b}_{i}=0) when a specific ({x}_{i}) was excluded.Based on this methodology, ecosystem flux data were extrapolated to the previous 7–15 years (since 2000 in the dry-subhumid and arid sites, since 2004 in the semi-arid sites, and since 2008 in the Oak-forest site) using all the available continuous meteorological parameters from the meteorological stations associated with our field sites. The long-term annual sums and means of extrapolated ecosystem fluxes were averaged for multi-year means of each site for the period of available extrapolated data. In the first, previously published phase of this study29, the extrapolation method was extensively tested, including simulation experiments at the arid forest site, where continuous flux data from the 20 years old permanent flux tower were available. Five percent of the daily data were selected by bootstrap (with 20 repetitions), a stepwise regression was performed for this sample, and then, the prediction of fluxes using Eq. 1 above was performed for the entire observation period. In the current phase of the study additional flux measurements are included with the R2 coefficients for the additional measurements ranging between 0.4-0.9, see Supplementary Table S3.The aridity index of the Oak-forest was in between those of the semi-arid and dry-subhumid Pine-forests (0.4 compared to 0.37 and 0.67, respectively). Therefore, to compare the Oak-forest with Pine-forest and non-forest sites, the average results from the semi-arid and dry-subhumid paired sites were used.Radiative forcing and carbon equivalence equationsTo compare the changes in the carbon and radiation budgets caused by forestation, we adopted the approach of Myhre et al. 30, and used Eq. 2:$${{RF}}_{triangle C}=5.35{{{{mathrm{ln}}}}}left(1+frac{triangle C}{{C}_{0}}right),left[W,{m}^{-2}right]$$
    (2)
    where land-use changes in radiative forcing (RFΔC) are calculated based on the CO2 reference concentration, C0 (400 ppm for the measured period of study), and ΔC, which is the change in atmospheric CO2 in ppm, with a constant radiative efficiency (RE) value of 5.35. Here, ΔC is calculated based on the annual net ecosystem productivity (NEP; positive carbon gain by the forest, which is identical to net ecosystem exchange (NEE), the negative carbon removal from the atmosphere) as the difference between forested and non-forested ecosystems (ΔNEP) multiplied by a unit conversion constant:$$triangle C={left[{overline{{NEP}}}_{F}-{overline{{NEP}}}_{{NF}}right]}_{[gC{m}^{-2}{y}^{-1}]}cdot k ,[{ppm}]$$
    (3)
    where, k is a unit conversion factor, from ppm to g C (k = 2.13 × 109), calculated as the ratio between the air molar mass (Ma = 28.95; g mol−1), to carbon molar mass (Mc = 12.0107; g mol−1), and total air mass (ma = 5.15 × 109; g).Etminan et al. 57 introduced an updated approach to calculate the RE as a co-dependent of the change in CO2 concentration and atmospheric N2O:$${RE}={a}_{1}{left(triangle Cright)}^{2}+{b}_{1}left|triangle Cright|+{c}_{1}bar{N}+5.36,left[W,{m}^{-2}right]$$
    (4)
    where, (triangle {{{{{rm{C}}}}}}) is the change in atmospheric CO2 in ppm resulting from the forestation, as calculated in Eq. 3, (bar{{{{{{rm{N}}}}}}}) is the atmospheric N2O concentration in ppb (323), and the coefficients a1, b1, and c1 are −2.4 × 10−7 Wm−2ppm−1, 7.2 × 10−4 Wm−2 ppm−1, and −2.1 × 10−4 Wm−2ppb−1, respectively.Combining Eqs. 2 and 4 with an airborne fraction of (zeta =0.44)58, we obtain Eq. 5:$${{RF}}_{triangle C}={RE}cdot {{{{mathrm{ln}}}}}left(1+frac{zeta cdot triangle C}{{C}_{0}}right),left[W,{m}^{-2}right]$$
    (5)
    Next, the annual average radiative forcing due to differences in radiation flux was calculated as follows:$${{RF}}_{triangle R}=frac{triangle Rcdot {A}_{F}}{{A}_{E}},left[W,{m}^{-2}right]$$
    (6)
    where, ΔR is the difference between forest and non-forest reflected short-wave or emitted long-wave radiation (ΔSWnet and ΔLWnet, respectively), assuming that the atmospheric incoming solar and thermal radiation fluxes are identical for the two, normalized by the ratio of the forest area (({A}_{F})) to the Earth’s area (({A}_{E}=5.1times {10}^{14},{m}^{2})).As forest conversion mostly has a lower albedo, the number of years needed to balance (‘Break even time’) the warming effect of changes in radiation budget by the cooling effect of carbon sequestration is calculated by combining Eqs. 5 and 6:$$^prime{Break},{even},{time}^prime=frac{{{RF}}_{triangle alpha }}{{{RF}}_{triangle C}},[{{{{{rm{years}}}}}}]$$
    (7)
    The multiyear averages of NEP for each of the three paired sites (forest and non-forest) were then modeled over a forest life span of 80 years. This was done based on a logarithmic model, modified for dryland, which takes as an input the long-term averages of NEP ((overline{{NEP}})) as in Eq. 3:$${{NEP}}_{t}=overline{{NEP}}(1-{exp }^{bcdot t}),left[g{{{{{rm{C}}}}}},{{{{{{rm{m}}}}}}}^{-2}{{yr}}^{-1}right]$$
    (8)
    where annual carbon gain at time t (NEPt) is a function of the multiyear average carbon gain ((overline{{{{{{rm{NEP}}}}}}})), forest age (t), and growth rate (b). Parameter b is a constant (b = −0.17) and is calculated based on the global analysis of Besnard et al. 59, limiting the data to only dryland flux sites60. Note that this analysis indicates NEP reaching a steady state and was used here to describe the initial forest growth phase, while growth analyses indicate that carbon sequestration peaks after about 80 years, followed by a steep decline50. This is consistent with the time scale for forest carbon sequestration considered here.In contrast to the one-year differences presented in Table 1 (ΔNEP), the net sequestration potential ((triangle)SP) for each of the paired sites was calculated as the accumulated ecosystem ΔNEPt along with forest age ((triangle) is the difference between forest and non-forest sites):$$triangle {SP}=mathop{sum }limits_{t=0}^{{age}}{triangle {NEP}}_{t}/100,left[{{{{{rm{tC}}}}}},{{{{{{rm{ha}}}}}}}^{-1}{{age}}^{-1}right]$$
    (9)
    The ΔSP growth model was compared with previously published data of long-term carbon stock changes in arid forests (i.e., cumulative NEP over 50 years since forest establishment, t = 50), demonstrating agreement within ± 10%27.For comparison with previous studies, the carbon emission equivalent of shortwave forcing (EESF) was calculated using an inverse version of Eqs. 5 and 6 based on Betts1:$${EESF}={C}_{0}left({e}^{frac{{{RF}}_{triangle R}}{zeta cdot {RE}}}-1right)cdot k/100,left[{{{{{rm{tC}}}}}},{{{{{{rm{ha}}}}}}}^{-1}{{age}}^{-1}right]$$
    (10)
    where, C0 is the reference atmospheric CO2 concentration (400 ppm, the average atmospheric concentration for the past decade), RFΔR is the multiyear average change in radiative forcing as a result of the change in surface albedo (Eq. 6 W m−2), RE is the radiative efficiency (Eq. 4, W m−2), ζ is the airborne fraction (0.44 as in Eq. 5), and k is a conversion factor, from ppm to g C (2.13 × 109 as in Eq. 3). Equation 10 was also used to calculate the emission equivalent of longwave forcing (EELF) with the RFΔR as the multiyear average change in radiative forcing as a result of the change in net long-wave radiation (ΔLWnet).Finally, the net equivalent change in carbon stock due to both the cooling effect of carbon sequestration and the warming effect due to albedo change (net equivalent stock change; NESC), was calculated by a simple subtraction:$${NESC}=triangle {SP}-{EESF},left[{{{{{rm{tC}}}}}},{{{{{{rm{ha}}}}}}}^{-1}{{age}}^{-1}right]$$
    (11)
    A comparison of the ΔSP (Eq. 9), the EESF (Eq. 10), and NESC (Eq. 11) with the same metrics as those used in other studies1,14,61 was done when appropriate. An exception was made for Arora & Montenegro (2011), where only carbon stock changes (ΔSP) were available in carbon units, and biogeophysical (BGP) and biogeochemical (BGC) effects were expressed as temperature changes. To overcome this metric difference, we converted the biogeophysical to carbon equivalent units (EESF + EELF) by multiplying the carbon stock changes (ΔSP) by the ratio between the BGP and BGC effects on temperature (EESF + EELF = ΔSP × BGP/BGC).Statistical and data analysesThe paired t-test was used to compare multi-annual averages of all variables between forested and adjacent non-forested sites and between sites across the climatic gradient. The variables of interest that were detected for their significant differences were albedo, net radiation and its longwave and shortwave components, latent heat fluxes, sensible heat fluxes, and net ecosystem productivity. All statistical and data analyses were performed using R 3.6.0 (R Core Team, 2020)62.Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

  • in

    Author Correction: Measuring the world’s cropland area

    Authors and AffiliationsStatistics Division, Food and Agriculture Organization of the United Nations, Rome, ItalyFrancesco N. Tubiello, Giulia Conchedda, Leon Casse & Giorgia De SantisDigitization and Informatics Division, Food and Agriculture Organization of the United Nations, Rome, ItalyHao Pengyu & Chen ZhongxinInternational Institute for Applied Systems Analysis, Laxenburg, AustriaSteffen FritzGeospatial Unit, Land and Water Division, Food and Agriculture Organization of the United Nations, Rome, ItalyDouglas MuchoneyAuthorsFrancesco N. TubielloGiulia ConcheddaLeon CasseHao PengyuChen ZhongxinGiorgia De SantisSteffen FritzDouglas MuchoneyCorresponding authorCorrespondence to
    Francesco N. Tubiello. More

  • in

    Spring phenology alters vegetation drought recovery

    Mishra, A. K. & Singh, V. P. J. Hydrol. 391, 202–216 (2010).Article 

    Google Scholar 
    Jiao, W. et al. Nat. Commun. 12, 3777 (2021).Article 
    CAS 

    Google Scholar 
    Gampe, D. et al. Nat. Clim. Change 11, 772–779 (2021).Article 

    Google Scholar 
    IPCC Climate Change 2021: The Physical Science Basis (eds Masson-Delmotte, V. et al.) (Cambridge Univ. Press, 2021).Schwalm, C. R. et al. Nature 548, 202–205 (2017).Article 
    CAS 

    Google Scholar 
    Li, Y. et al. Nat. Clim. Change https://doi.org/10.1038/s41558-022-01584-2 (2023).Fourth National Climate Assessment: Volume II—Impacts, Risks, and Adaptation in the United States (US Global Change Research Program, 2018).Daryanto, S., Wang, L. & Jacinthe, P. A. PLoS ONE 11, e0156362 (2016).Article 

    Google Scholar 
    Jiao, W. et al. J. Geophys. Res. Biogeosci. 127, e2021JG006431 (2022).Augspurger, C. K. Oecologia 156, 281–286 (2008).Article 

    Google Scholar 
    Lian, X. et al. Nat. Commun. 12, 983 (2021).Article 
    CAS 

    Google Scholar 
    Buermann, W. et al. Nature 562, 110–114 (2018).Article 
    CAS 

    Google Scholar 
    Lian, X. et al. Sci. Adv. 6, eaax0255 (2020).Article 

    Google Scholar 
    Jiao, W., Wang, L. & McCabe, M. F. Rem. Sens. Environ. 256, 112313 (2021).Article 

    Google Scholar  More