Mark Twain National Forest Ecological Forecasting: Using NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest
Explore the source record for details and available documents.
Kaitlyn Bretz↗
SEARCH · Search NASA
Results for “forest ecology”
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.
Mark Twain National Forest Ecological Forecasting - Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest
The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.
Kaitlyn Bretz↗
Mark Twain National Forest Ecological Forecasting: Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest
The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.
DEVELOP Tech Paper↗
Forest Attributes from Radar Interferometric Structure and its Fusion with Optical Remote Sensing
The possibility of global, three-dimensional remote sensing of forest structure with interferometric synthetic aperture radar (InSAR) bears on important forest ecological processes, particularly the carbon cycle. InSAR supplements two-dimensional remote sensing with information in the vertical dimension. Its strengths in potential for global coverage complement those of lidar (light detecting and ranging), which has the potential for high-accuracy vertical profiles over small areas. InSAR derives its sensitivity to forest vertical structure from the differences in signals received by two, spatially separate radar receivers. Estimation of parameters describing vertical structure requires multiple-polarization, multiple-frequency, or multiple-baseline InSAR. Combining InSAR with complementary remote sensing techniques, such as hyperspectral optical imaging and lidar, can enhance vertical-structure estimates and consequent biophysical quantities of importance to ecologists, such as biomass. Future InSAR experiments will supplement recent airborne and spaceborne demonstrations, and together with inputs from ecologists regarding structure, they will suggest designs for future spaceborne strategies for measuring global vegetation structure.
InSAR↗
Landsat TM as a Tool for Locating Habitat for Cerulean Warblers
I believe that I made significant strides in three areas between fall of 1997 and fall of 2000 when I concluded my participation in the JOVE program. First, I acquired skill in digital remote sensing. This was significant to me because it had been 20 years since I had done any work utilizing remote sensing. I used my new skills in two classroom settings (forest ecology and GIS). In addition, I will participate as an instructor of digital remote sensing in a workshop for secondary educators this coming spring. Second, I received funding from the Arkansas Game and Fish Commission and the U.S. Forest Service to supplement JOVE funds. Third, and most importantly, a students and I developed a technique using LandSAT TM for identifying habitat for cerulean warblers. We developed a habitat model using logistic regression to discriminate between pixels that had a high probability of representing good cerulean warbler habitat and pixels that had a low probability of representing cerulean warbler habitat. Using this model, we located five significant populations of cerulean warblers in the Ozark National Forest of Arkansas. These populations were unknown before the initiation of this research and further represent a significant proportion of the known cerulean warblers in Arkansas. Preliminary findings were presented at the Ornithological Societies of America meeting in August of 1999. I also presented findings at the Arkansas Game and Fish Commission Research Symposium held in June of 2000. Finally, one paper is in press: James, D. A., C.J. Kellner, J. Self, and J. Davis., 'Breeding season distribution of cerulean warblers in Arkansas in the 1990's'. In addition, one paper is under construction: 'Population fluctuation and habitat selection by cerulean warblers in upland forests of Arkansas,' and one paper is under consideration: 'LandSAT TM and Logistic regression for identification of cerulean warbler habitat in upland forests of Arkansas.'
Kellner, Chris↗
Postfire Biogeochemical Processes: Implications to Source Water Quality in Fire-Influenced Watersheds
Forested watersheds are instrumental in providing purified and reliable water to millions of people worldwide. The changing climate has increased the frequency and severity of global fire events. Forested watersheds and their ecosystem functions are greatly disrupted during fire activity. Postfire concerns in forested watersheds include unpredictable and potentially simultaneous alterations in source water quality and hydro-biogeochemical processes. Here, the degree of fire severity can complexly modify water quality through the production of fire-transformed constituents on the burned forest floor (i.e., nutrients, metal(loid)s, dissolved organic matter, and the formation of disinfection byproducts). Correspondingly, fire severity and postfire rainfall patterns can refine hydro-biogeochemical processes that influence the transport of the fire-transformed constituents (i.e., vegetation function, soil structure, hydrological pathways, and microbial communities). Postfire alterations to water quality and hydro-biogeochemical processes introduce further complexity with varying temporal influence, which ranges from months to decades. As postfire water quality and watershed response research progresses, it is essential to homogenize interdisciplinary expertise to bridge knowledge gaps between fields ranging from forest ecology, hydrology, microbiology, and geochemistry. A multidisciplinary approach in wildfire research will facilitate a comprehensive perception of the diverse water quality risks associated with fire activity and mitigate fire concerns on a global level.
Disinfection Byproducts↗
Carbon cycling across ecosystem succession in a north temperate forest: Controls and management implications
Despite decades of progress, much remains unknown about successional trajectories of carbon (C) cycling in north temperate forests. Drivers and mechanisms of these changes, including the role of different types of disturbances, are particularly elusive. To address this gap, we synthesized decades of data from experimental chronosequences and long-term monitoring at a well-studied, regionally representative field site in northern Michigan, USA. Our study provides a comprehensive assessment of changes in above- and belowground ecosystem components over two centuries of succession, links temporal dynamics in C pools and fluxes with underlying drivers, and offers several conceptual insights to the field of forest ecology. Our first advance shows how temporal dynamics in some ecosystem components are consistent across severe disturbances that reset succession and partial disturbances that slightly modify it: both of these disturbance types increase soil N availability, alter fungal community composition, and alter growth and competitive interactions between short-lived pioneer and longer-lived tree taxa. Further, these changes in turn affect soil C stocks, respiratory emissions, and other belowground processes. Second, we show that some other ecosystem components have effects on C cycling that are not consistent over the course of succession. For example, canopy structure does not influence C uptake early in succession but becomes important as stands develop, and the importance of individual structural properties changes over the course of two centuries of stand development. Third, we show that in recent decades, climate change is masking or overriding the influence of community composition on C uptake, while respiratory emissions are sensitive to both climatic and compositional change. In synthesis, we emphasize that time is not a driver of C cycling; it is a dimension within which ecosystem drivers such as canopy structure, tree and microbial community composition change. Changes in those drivers, not in forest age, are what control forest C trajectories, and those changes can happen quickly or slowly, through natural processes or deliberate intervention. Stemming from this view and a whole-ecosystem perspective on forest succession, we offer management applications from this work and assess its broader relevance to understanding long-term change in other north temperate forest ecosystems.
54 ENVIRONMENTAL SCIENCES↗
The use of four band multispectral photography to identify forest cover types
Four-band multispectral aerial photography and a color additive viewer were employed to identify forest cover types in Northern Alabama. The multispectral photography utilized the blue, green, red and near-infrared spectral regions and was made with black and white infrared film. On the basis of color differences alone, a differentiation between conifers and hardwoods was possible; however, supplementary information related to forest ecology proved necessary for the differentiation of various species of pines and hardwoods.
Downs, S. W., Jr.↗
Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems
Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. Here, we showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.
59 BASIC BIOLOGICAL SCIENCES↗
Integrating Forest Structural Diversity Measurement Into Ecological Research
The measurement of forest structure has evolved steadily due to advances in technology, methodology, and theory. Such advances have greatly increased our capacity to describe key forest structural elements and resulted in a range of measurement approaches from traditional analog tools such as measurement tapes to highly derived and computationally intensive methods such as advanced remote sensing tools (e.g., lidar, radar). This assortment of measurement approaches results in structural metrics unique to each method, with the caveat that metrics may be biased or constrained by the measurement approach taken. While forest structural diversity (FSD) metrics foster novel research opportunities, understanding how they are measured or derived, limitations of the measurement approach taken, as well as their biological interpretation is crucial for proper application. We review the measurement of forest structure and structural diversity—an umbrella term that includes quantification of the distribution of functional and biotic components of forests. We consider how and where these approaches can be used, the role of technology in measuring structure, how measurement impacts extend beyond research, and current limitations and potential opportunities for future research.
Jeff W. Atkins↗
Maine Ecological Forecasting II: Identifying Forest Cover and Assessing Federally Endangered Atlantic Salmon Habitat in Maine Using Earth Observations
Major declines in Atlantic salmon (Salmo salar) populations have occurred alongside dam construction, shifting temperatures, and changing land use and land cover (LULC), restricting the last remaining wild population in the United States to Maine. In collaboration with the Maine Department of Marine Resources and the Downeast Salmon Federation, the NASA DEVELOP team utilized Earth observations to assess changes within critical salmon habitat. The team used the Landsat 5 Thematic Mapper (TM), Landsat 8 Operation Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), and National Land Cover Database (NLCD) to refine historical LULC maps from 1985 to 2021. The team also used elevation data from the Shuttle Radar Topography Mission (SRTM) and IDRISI TerrSet Land Change Modeler to forecast LULC change to 2040. From Terra Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) data, the team derived a time series, summer averages, and anomaly maps between 2000 and 2021. The team analyzed LST in relation to LULC classes, specifically forest cover type, and in-situ stream temperature from the Spatial Hydro-Ecological Decision System (EcoSHEDS). LULC trends from 1985 to 2021 revealed a net transition of coniferous forest to other classes, such as deciduous forest and developed land. The team found an association between warmer temperatures and a greater presence of developed and mixed forest classes per 10,000 acres across Maine. These visualizations and analyses will aid partners in identifying riparian locations with ideal or unfavorable habitat conditions to inform Atlantic salmon population recovery efforts.
Tony Bowman↗
A primer on forest structure measurement with lidar for ecologists
Light detection and ranging (lidar) technology has fundamentally advanced the way we measure forest structure, facilitating new insights into ecological processes. Lidar for forest ecology applications is deployed on multiple types of platforms that operate from the ground, air, or space, and each has associated strengths and limitations. Ideally, the choice of what kind of lidar to use in a particular study should be guided by the ecological question of interest; however, practical considerations of cost, data availability, and processing tools can be equally important. This synthesis is a practical introduction to how different lidar platforms characterize forest structure (e.g., tree size/location, wood volume, branching structure, aboveground biomass, leaf properties), designed for a general audience of ecologists (not remote sensing scientists) seeking an accessible introduction to the use of lidar. We also provide examples of novel ecological insights from recent lidar research and describe current limitations and areas of expected improvement. Last, we include an appendix of data collected from terrestrial, mobile, unoccupied aerial system, airplane, and satellite lidar platforms within a common temperate forest area, with associated code to allow new lidar users to visualize and manipulate data in R.
Cushman, KC [ORNL] (ORCID:0000000234641151)↗
Mangrove Canopy Height Globally Related to Precipitation and Cyclone Frequency
Mangrove wetlands are among the most productive and carbon-dense ecosystems in the world. Their structural attributes vary considerably across spatial scales, yielding large uncertainties in regional and global estimates of carbon stocks. Here, we present a global analysis of mangrove canopy height gradients and aboveground carbon stocks based on remotely sensed measurements and field data. Our study highlights that precipitation, temperature and cyclone frequency explain 74% of the global trends in maximum canopy height, with other geophysical factors influencing the observed variability at local and regional scales. We find the tallest mangrove forests in Gabon, equatorial Africa, where stands attain 62.8 m. The total global mangrove carbon stock (above- and belowground biomass, and soil) is estimated at 5.03 Pg, with a quarter of this value stored in Indonesia. Our analysis implies sensitivity of mangrove structure to climate change, and offers a baseline to monitor national and regional trends in mangrove carbon stocks.
Carbon cycle↗
Coastal groundwater salinization impairs tree carbon–water balance
How coastal forest productivity varies with local nutrient availability and water supply remains a knowledge gap under climate change. In a two-decade field experiment manipulating fertilization and density in coastal pine forests, we show that a decade of growth enhancement by simulated sedimentary nutrient inputs has resulted in a striking reversal in growth and increased mortality risk as drought and sea level rise progressed. Recent groundwater salinization has further triggered a shift from nutrient to water limitation, causing severe stomatal closure and decoupling of tree carbon–water balance, which induces a negative intrinsic water use efficiency (iWUE)–growth relationship. Higher tree iWUE predicted sharper tree growth declines, and both nutrient enrichment and high stand density amplified this feedback, increasing the risk of hydraulic failure and mortality. These results suggest that a transient nutrient-stimulatory effect could drive further coastal forest degradation due to heightened belowground saltwater stress under sea-level rise.
Climate-change ecology↗
Hydraulic Diversity of Forests Regulates Ecosystem Resilience During Drought
Plants influence the atmosphere through fluxes of carbon, water and energy, and can intensify drought through land–atmosphere feedback effects. The diversity of plant functional traits in forests, especially physiological traits related to water (hydraulic) transport, may have a critical role in land–atmosphere feedback, particularly during drought. Here we combine 352 site-years of eddy covariance measurements from 40 forest sites, remote-sensing observations of plant water content and plant functional-trait data to test whether the diversity in plant traits affects the response of the ecosystem to drought. We find evidence that higher hydraulic diversity buffers variation in ecosystem flux during dry periods across temperate and boreal forests. Hydraulic traits were the predominant significant predictors of cross-site patterns in drought response. By contrast, standard leaf and wood traits, such as specific leaf area and wood density, had little explanatory power. Our results demonstrate that diversity in the hydraulic traits of trees mediates ecosystem resilience to drought and is likely to have an important role in future ecosystem–atmosphere feedback effects in a changing climate.
Climate and Earth system modelling↗
Data from TropiRoot 1.0 database: tropical root characteristics across environments
TropiRoot 1.0 is a new tropical root database with root characteristics across environment gradients. It has data extracted from 104 new sources, resulting in more than 8000 rows of data (either species or community data). Most of the data in TropiRoot 1.0 includes root characteristics such as root biomass, morphology, root dynamics, mass fraction, architecture, anatomy, physiology and root chemistry. This initiative represents an approximately 30% increase in the currently available data for tropical roots in the Fine Root Ecology Database (FRED). TropiRoot 1.0, contains root characteristics from 25 different countries where seven are located in Asia, six in South America, five in Central America and the Caribbean, four in Africa, two in North America, and 1 in Oceania. Due to the volume of data, when ancillary data was available, including soil data, these data was either extracted and included in the database or their availability was recorded in an additional column. Multiple contributors checked the entries for outliers during the collation process to ensure data quality. For text-based observations, we examined all cells to ensure that their content relates to their specific categories. For numerical observations, we ordered each numerical value from least to greatest and plotted the values, checking apparent outliers against the data in their respective sources and correcting or removing incorrect or impossible values. Some data (soil and aboveground) have different columns for the same variable presented in different units, including originally published units, but root characteristics data had units converted to match the ones reported in FRED. By filling a gap from global databases, TropiRoot 1.0 expands our knowledge of otherwise so far underrepresented regions, and our ability to assess global trends. This advancement can be used to improve tropical forest representation in vegetation models.
54 ENVIRONMENTAL SCIENCES↗
Maine Ecological Forecasting II: Identifying Forest Cover and Assessing Federally Endangered Atlantic Salmon Habitat in Maine Using Earth Observations
Explore the source record for details and available documents.
Tony Bowman↗
CORPSE model with litter decomposition parameters derived from the LIDET dataset
This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436
Juice, Stephanie↗