Search NASA⌕ Search

SEARCH · Search NASA

Results for “litter”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Invasions eliminate the legacy effects of substrate history on microbial nitrogen cycling

Abstract Changes in substrate quality driven by climate, land use, or other forms of global change may represent a strong selective force on microbial communities. Invasion of new taxa into a community through dispersal, evolution, or recolonization could impact the outcome of this environmental selection. Here, we simulated substrate change with a trait‐based model of microbial litter decomposition (DEMENTpy) to assess the legacy effects of past substrate quality and the impact of selection by a new substrate on community decomposition activity. Simulations were run with different levels of invasion, including invasion from communities long‐adapted to the new substrate. Legacy effects were evident with substrate change for native communities differing in composition. Protein was the only substrate that exerted a strong enough selective force to affect community composition. Legacy effects disappeared when invaders came from substrates similar to the new substrate. Together, our simulations demonstrate that substrate quality changes associated with global change can lead to legacy effects on substrate degradation. In decomposing plant litter, such legacy effects can occur if substrate inputs shift to higher protein content and if invasion is low.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old‐growth to second‐growth forests in headwater ecosystems

Abstract Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear‐cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old‐growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear‐cut harvesting, from mixed‐species old‐growth to lodgepole pine‐dominated second‐growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3‐year period (2016–2018) at the Fraser Experimental Forest in north‐central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old‐ versus second‐growth forest. Fluorescence spectroscopy showed that the DOM exported from the old‐growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second‐growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old‐growth forest relative to second growth. Old‐growth and second‐growth forests are co‐mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

Fegel, Timothy S.↗
DOE OSTI2022

Respiration regimes in rivers: Partitioning source–specific respiration from metabolism time series

Respiration in streams is controlled by the timing, magnitude, and quality of organic matter (OM) inputs from internal primary production and external fluxes. Here, we estimated the contribution of different OM sources to seasonal, annual, and event-driven characteristics of whole-stream ecosystem respiration (ER) using an inverse modeling framework that accounts for possible time-lags between OM inputs and respiration. We modeled site-specific, dynamic OM stocks contributing to ER: autochthonous OM from gross primary production (GPP); allochthonous OM delivered during flow events; and seasonal pulses of leaf litter. OM stored in the sediment and dissolved organic matter (DOM) transported during baseflow were modeled as a stable stock contributing to baseline respiration. We applied this modeling framework to five streams with different catchment size, climate, and canopy cover, where multi-year time series of ER and environmental variables were available. Overall, the model explained between 53% and 74% of observed ER dynamics. Respiration of autochthonous OM tracked seasonal peaks in GPP in spring or summer. Increases in ER were often associated with high-flow events. Respiration associated with litter inputs was larger in smaller streams. Time lags between leaf inputs and respiration were longer than for other OM sources, likely due to lower biological reactivity. Model estimates of source-specific ER and OM stocks compared well with existing measures of OM stocks, inputs, and respiration or decomposition. Our modeling approach has the potential to expand the scale of comparative analyses of OM dynamics within and among freshwater ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Salinity exposure affects lower-canopy specific leaf area of upland trees in a coastal deciduous forest

Sea level rise and increasing storm surges are likely to affect the canopy physiology, ecology, and structure of coastal forests, even well in advance of tree mortality. Laboratory and greenhouse studies have documented that saltwater exposure can trigger changes in leaf-level physiology and morphology, but few in situ studies have examined how tree-specific leaf area (SLA), the ratio of leaf area to mass and a crucial trait and model parameter, is affected by saline soils. We conducted an observational study of SLA in a mid-Atlantic (USA) coastal deciduous forest, taking advantage of a natural gradient in salinity along a tidal creek. Measured SLA of the 239 trees and seven species sampled ranged from Carya glabra (N = 6 trees, mean SLA = 277.9±36.3cm 2 /g) to Fagus grandifolia (N=60, 321.9±62.9cm 2 g); as expected, trees species and canopy position (sun versus shade) significantly affected SLA. For trees (N=100) directly exposed to the tidal creek, salinity was highly significant after accounting for species (P<0.001), with trees in the lower reaches of the creek having lower SLA. Leaf area index (LAI), computed from SLA and litter traps, ranged from 4.8 to 15.8 and was inversely related to salinity exposure; the spatial variability in leaf litter production contributed much more to LAI uncertainty than did SLA variability. These in situ results are correlative but consistent with the hypothesis, based on previous greenhouse studies, that the stress of chronic salinity exposure changes species’ leaf morphology. Our findings are useful for understanding the growing effects of saltwater intrusion into upland forests, as well as parameterizing and testing ecosystem-scale models simulating forest stressors and disturbances at the terrestrial-aquatic interface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Soil organic matter is principally root derived in an Ultisol under oak forest

To a large degree, the sources and stability of soil organic carbon remain poorly constrained. A clear understanding of links among the components of the soil C cycle is hampered by the complexity of the system as well as challenges associated with partitioning bulk soil C into meaningful fractions. A large accidental 14 CO 2 release at the Oak Ridge Reservation in Tennessee, USA provided a strong label pulse into adjacent, well-studied oak forests, resulting in highly elevated Δ 14 C values in leaf litter (~1000‰) and roots (~260–450‰). A four-year manipulative study was conducted to determine the relative contribution of litter versus roots to the bulk mineral soil C pool, as well as to free light, occluded light and heavy fractions. The heavy fraction was further split into fractions with densities of 1.7–2.4 g cm –3 and >2.4 g cm –3 to test the homogeneity of the mineral-associated fraction of C. Substantial concentrations of label were detected in all soil fractions within a year of the 14 CO 2 release, indicating rapid incorporation of newly fixed photosynthates in all fractions of soil organic C, regardless of differences in stability inferred by previous work. This rapid incorporation of label occurred only in treatments where roots were labeled, indicating that roots are the major source of inputs to mineral soil C stocks at these sites. Separation of the heavy fraction into subfractions of intermediate (1.7–2.4 g cm –3 ) and high (>2.4 g cm –3 ) density indicated that both subfractions incorporated label at similar rates, despite significant differences in degree of microbial processing. In general, the rate of label incorporation suggested a much faster turnover for all fractions than indicated by natural radiocarbon abundance values. This study suggests that within each soil fraction there are portions of slow-cycling and fast-cycling materials, and the determination of an average turnover time or mean age is dependent on experimental approach. The rapid incorporation of label into all fractions within a year regardless of inferred stability implies a high degree of heterogeneity in all fractions regardless of how finely the soils are partitioned. Further refinement of the nature and drivers of this heterogeneity could yield important insights into the soil C cycle.

58 GEOSCIENCES↗
DOE OSTI2020

Shale as a Source of Organic Carbon in Floodplain Sediments of a Mountainous Watershed

Abstract Shales contain high levels of organic carbon (OC) and represent a large fraction of the Earth's reduced carbon stocks. While recent evidence suggests that shale‐derived OC may be actively cycled in riverine systems, this process is poorly understood and not currently considered in global C models. Through the use of sediment density fractionations, extractions, radiocarbon measurements, and chemical characterization, we provide information on the abundance, chemistry, and mobility of shale‐derived OC in floodplain sediments of a shale‐rich mountainous watershed. The heavy fraction of the sediment, representing mineral‐associated OC, is the largest (84 ± 6% of TOC) and oldest (Δ 14 C values −224 to −853‰) OC pool. Evidence of shale‐derived OC is observed in all sediment C pools (i.e., occluded light fraction, water‐soluble, and pyrophosphate‐extractable) except the free light fraction, which is entirely modern. Relatively consistent chemistry was observed across samples for extracted and density‐separated OC, despite wide ranges of Δ 14 C values. Carbon spectroscopy revealed that floodplain sediments had a higher degree of functionalized aromatic groups and lower carbonate content compared to shale collected nearby, consistent with chemical alteration and mixing with other C sources in the floodplain. We estimate that approximately 23–34% of sediment OC is derived from shale, with implications for other shale‐derived elements (e.g., N). This study demonstrates the important contribution of shale‐OC, particularly in environments with low litter inputs. The large impact of radiocarbon‐dead shale‐OC, which has a thermally altered chemical structure distinct from plant litter, on Δ 14 C values and reactivity of sediment‐OC must be considered.

Fox, Patricia M.↗
DOE OSTI2026

Hydrological Control on Soil Redox Condition and Carbon Loss of Coastal Wetland Under Sea-Level Rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Chen, Kewei [ORNL] (ORCID:000000032580514X)↗
DOE OSTI2023

Contrasting geochemical and fungal controls on decomposition of lignin and soil carbon at continental scale

Lignin is an abundant and complex plant polymer that may limit litter decomposition, yet lignin is sometimes a minor constituent of soil organic carbon (SOC). Accounting for diversity in soil characteristics might reconcile this apparent contradiction. Tracking decomposition of a lignin/litter mixture and SOC across different North American mineral soils using lab and field incubations, here we show that cumulative lignin decomposition varies 18-fold among soils and is strongly correlated with bulk litter decomposition, but not SOC decomposition. Climate legacy predicts decomposition in the lab, and impacts of nitrogen availability are minor compared with geochemical and microbial properties. Lignin decomposition increases with some metals and fungal taxa, whereas SOC decomposition decreases with metals and is weakly related with fungi. Decoupling of lignin and SOC decomposition and their contrasting biogeochemical drivers indicate that lignin is not necessarily a bottleneck for SOC decomposition and can explain variable contributions of lignin to SOC among ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A global dataset of terrestrial biological nitrogen fixation

Biological nitrogen fixation (BNF) is the main natural source of new nitrogen inputs in terrestrial ecosystems, supporting terrestrial productivity, carbon uptake, and other Earth system processes. We assembled a comprehensive global dataset of field measurements of BNF in all major N-fixing niches across natural terrestrial biomes derived from the analysis of 376 BNF studies. The dataset comprises 32 variables, including site location, biome type, N-fixing niche, sampling year, quantification method, BNF rate (kg N ha −1 y −1 ), the percentage of nitrogen derived from the atmosphere (%N dfa ), N fixer or N-fixing substrate abundance, BNF rate per unit of N fixer abundance, and species identity. Overall, the dataset combines 1,207 BNF rates for trees, shrubs, herbs, soil, leaf litter, woody litter, dead wood, mosses, lichens, and biocrusts, 152 herb %N dfa values, 1,005 measurements of N fixer or N-fixing substrate abundance, and 762 BNF rates per unit of N fixer abundance for a total of 424 species across 66 countries. This dataset facilitates synthesis, meta-analysis, upscaling, and model benchmarking of BNF fluxes at multiple spatial scales.

Reis Ely, Carla R. [Oregon State Univ., Corvallis,↗
DOE OSTI2020

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity

Abstract As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates—namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose–methanol–choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases—we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Soil resistance and recovery during neotropical forest succession

The recovery of soil conditions is crucial for successful ecosystem restoration and, hence, for achieving the goals of the UN Decade on Ecosystem Restoration. Here, we assess how soils resist forest conversion and agricultural land use, and how soils recover during subsequent tropical forest succession on abandoned agricultural fields. Our overarching question is how soil resistance and recovery depend on local conditions such as climate, soil type and land-use history. For 300 plots in 21 sites across the Neotropics, we used a chronosequence approach in which we sampled soils from two depths in old-growth forests, agricultural fields (i.e. crop fields and pastures), and secondary forests that differ in age (1–95 years) since abandonment. We measured six soil properties using a standardized sampling design and laboratory analyses. Soil resistance strongly depended on local conditions. Croplands and sites on high-activity clay (i.e. high fertility) show strong increases in bulk density and decreases in pH, carbon (C) and nitrogen (N) during deforestation and subsequent agricultural use. Resistance is lower in such sites probably because of a sharp decline in fine root biomass in croplands in the upper soil layers, and a decline in litter input from formerly productive old-growth forest (on high-activity clays). Soil recovery also strongly depended on local conditions. During forest succession, high-activity clays and croplands decreased most strongly in bulk density and increased in C and N, possibly because of strongly compacted soils with low C and N after cropland abandonment, and because of rapid vegetation recovery in high-activity clays leading to greater fine root growth and litter input. Furthermore, sites at low precipitation decreased in pH, whereas sites at high precipitation increased in N and decreased in C : N ratio. Extractable phosphorus (P) did not recover during succession, suggesting increased P limitation as forests age. These results indicate that no single solution exists for effective soil restoration and that local site conditions should determine the restoration strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Human activities shape global patterns of decomposition rates in rivers

Rivers and streams contribute to global carbon cycling by decomposing immense quantities of terrestrial plant matter. However, decomposition rates are highly variable and large-scale patterns and drivers of this process remain poorly understood. Using a cellulose-based assay to reflect the primary constituent of plant detritus, we generated a predictive model (81% variance explained) for cellulose decomposition rates across 514 globally distributed streams. A large number of variables were important for predicting decomposition, highlighting the complexity of this process at the global scale. Predicted cellulose decomposition rates, when combined with genus-level litter quality attributes, explain published leaf litter decomposition rates with high accuracy (70% variance explained). Finally, our global map provides estimates of rates across vast understudied areas of Earth and reveals rapid decomposition across continental-scale areas dominated by human activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Testing the contribution of dispersal to microbial succession following a wildfire

ABSTRACT Given increased wildfire activity, there is growing interest in understanding the drivers of microbial succession after fire. Dispersal may be especially important to post-fire succession as biotic communities can be more susceptible to invasion following a disturbance. Here, we experimentally manipulated dispersal into disturbed leaf litter communities collected following a wildfire and tracked bacterial and fungal dispersal assemblages over time. We show that the identity and source of microbes immigrating into the soil surface post-fire change across time with seasonal shifts and the reemergence of aboveground vegetation. Further, dispersal significantly contributed to the reassembly of leaf litter microbial communities after the fire, producing an increasingly distinct assembly trajectory. The effect of dispersal on α-diversity and β-diversity was ecosystem dependent but, unexpectedly, influenced bacterial and fungal communities in a similar manner within ecosystems. Collectively, these results demonstrate that dispersal explicitly alters the course of microbial community succession following a wildfire and may impact bacteria and fungi in parallel ways, despite differing in traits expected to alter dispersal patterns. IMPORTANCE Identifying the mechanisms underlying microbial community succession is necessary for predicting how microbial communities, and their functioning, will respond to future environmental change. Dispersal is one mechanism expected to affect microbial succession, yet the difficult nature of manipulating microorganisms in the environment has limited our understanding of its contribution. Using a dispersal exclusion experiment, this study isolates the specific effect of environmental dispersal on bacterial and fungal community assembly over time following a wildfire. The work demonstrates the potential to quantify dispersal impacts on soil microbial communities over time and test how dispersal might further interact with other assembly processes in response to environmental change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

DayCent MUVP model calibrated for the simulation of bioenergy crops grown on set-aside land in the US

DayCent MUVP version (Methanogenesis, UV litter degradation and Photosynthesis). DAYCENT is the daily time-step version of the CENTURY biogeochemical model (Parton et al., 1994). DAYCENT simulates fluxes of C and N among the atmosphere, vegetation, and soil (Del Grosso et al., 2001a; Parton et al., 1998). Key submodels include soil water content and temperature by layer, plant production and allocation of net primary production (NPP), decomposition of litter and soil organic matter, mineralization of nutrients, N gas emissions from nitrification and denitrification, and CH4 oxidation in non-saturated soils.

biogeochemical model↗
DOE OSTI2021

Microbial community response to a decade of simulated global changes depends on the plant community

Global changes such as increased drought and atmospheric nitrogen deposition perturb both the microbial and plant communities that mediate terrestrial ecosystem functioning. However, few studies consider how microbial responses to global changes may be influenced by interactions with plant communities. To begin to address the role of microbial–plant interactions, we tested the hypothesis that the response of microbial communities to global change depends on the plant community. We characterized bacterial and fungal communities from 395 plant litter samples taken from the Loma Ridge Global Change Experiment, a decade-long global change experiment in Southern California that manipulates rainfall and nitrogen levels across two adjacent ecosystems, a grassland and a coastal sage scrubland. The differences in bacterial and fungal composition between ecosystems paralleled distinctions in plant community composition. In addition to the direct main effects, the global change treatments altered microbial composition in an ecosystem-dependent manner, in support of our hypothesis. The interaction between the drought treatment and ecosystem explained nearly 5% of the variation in bacterial community composition, similar to the variation explained by the ecosystem-independent effects of drought. Unexpectedly, we found that the main effect of drought was approximately four times as strong on bacterial composition as that of nitrogen addition, which did not alter fungal or plant composition. Overall, the findings underscore the importance of considering plant–microbe interactions when considering the transferability of the results of global change experiments across ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Soil and soil fraction radiocarbon, C and N from the Enriched Background Isotope Study, Oak Ridge National Laboratory, Tennessee, USA

A large accidental 14CO2 release at the Oak Ridge Reservation in Tennessee, USA provided a strong label pulse into adjacent, well-studied oak forests, resulting in highly elevated Δ14C values in leaf litter (~1000‰) and roots (~260–450‰). A four-year manipulative study was conducted to provide insight into mechanisms underlying soil Organic matter dynamics, specifically to determine the relative contribution of litter versus roots to the bulk mineral soil C pool, as well as to free light, occluded light and heavy fractions. The heavy fraction was further split into fractions with densities of 1.7–2.4 g cm-3 and >2.4 g cm-3 to test the homogeneity of the mineral-associated fraction of C. Substantial concentrations of label were detected in all soil fractions within a year of the 14CO2 release, indicating rapid incorporation of newly fixed photosynthates in all fractions of soil organic C, regardless of differences in stability inferred by previous work. This rapid incorporation of label occurred only in treatments where roots were labeled, indicating that roots are the major source of inputs to mineral soil C stocks at these sites. Separation of the heavy fraction into subfractions of intermediate (1.7–2.4 g cm -3) and high (>2.4 g cm-3) density indicated that both subfractions incorporated label at similar rates, despite significant differences in degree of microbial processing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Soil organic matter, tree communities, and fungal communities across mycorrhizal gradients in the Eastern United States

We collected these data to investigate how mycorrhizal associations are related to soil C and N across four sites in the eastern U.S. broadleaf forest biome, which capture broad variability in climate and tree species. Our sites represented each of the four ecoregions in the eastern United States temperate forest—warm continental (New Hampshire), hot continental (Wisconsin), Prairie (Illinios), and subtropical (Georgia). This breadth naturally provided tree species diversity and allowed us to also investigate how soil C and N are related to canopy tree and EcM fungal community composition. Specifically, we analyzed the effects of mycorrhizal association, canopy tree family, and EcM fungal taxa on the proportion of C and N associated with soil minerals (MAOM) and bulk soil C:N. These are plot level data from forests in New Hampshire, Georgia, Wisconsin, and Illinois. The plots are 10 m in radius and vary across gradients of ectomycorrhizal tree basal area. Within each plot, data include the basal area of all trees >2 cm in diameter identified to species, the abundance of different fungal taxa based on ITS sequences, and the %C, %N, C:N of bulk soil and of soil density fractions (free light, occluded light, and particulate) as well as the proportion of soil C and N within each fraction. The soil data (both organic matter and fungal communities) are from the top 10 cm of the mineral soil. This also includes a meta-analysis of leaf litter k and C:N values for the tree species present in the forests.The "Tree_BA_ALL.csv" files include the basal area of all trees >2 cm in diameter at breast height within our 10 cm radii plots. The "Soil_Data_Final.csv" file contains the soil carbon and nitrogen in bulk and density fractions (concentrations as well as proportion for free light fraction (fLF), occluded light fraction (oLF), and heavy or dense fraction (HF), oxalate extractable iron and aluminum concentrations, as well as the geographic coordinates of each plot's location. The "FungiNonRelative.csv" file contains the abundances of fungal taxa identified via ITS sequencing for each plot (with the exception of 4 plots from GA where there was not enough quality DNA to extract).The "R_ReadyK.csv" and "R_Ready_CN.csv" files contain leaf litter decomposition constant (k) and C:N values from the literature for the tree species found at our sites.

54 ENVIRONMENTAL SCIENCES↗