Search NASA⌕ Search

SEARCH · Search NASA

Results for “tree decomposition”

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

On-shell structures of MHV amplitudes beyond the planar limit

We initiate an exploration of on-shell functions in N=4 SYM beyond the planar limit by providing compact, combinatorial expressions for all leading singularities of MHV amplitudes and showing that they can always be expressed as a positive sum of differently ordered Parke-Taylor tree amplitudes. This is understood in terms of an extended notion of positivity in G(2, n), the Grassmannian of 2-planes in n dimensions: a single on-shell diagram can be associated with many different “positive” regions, of which the familiar G + (2, n) associated with planar diagrams is just one example. The decomposition into Parke-Taylor factors is simply a “triangulation” of these extended positive regions. The U(1) decoupling and Kleiss-Kuijf (KK) relations satisfied by the Parke-Taylor amplitudes also follow naturally from this geometric picture. Finally, these results suggest that non-planar MHV amplitudes in N=4 SYM at all loop orders can be expressed as a sum of polylogarithms weighted by color factors and (unordered) Parke-Taylor amplitudes.

1/N Expansion↗

Management issues in systems engineering

When applied to a system, the doctrine of successive refinement is a divide-and-conquer strategy. Complex systems are sucessively divided into pieces that are less complex, until they are simple enough to be conquered. This decomposition results in several structures for describing the product system and the producing system. These structures play important roles in systems engineering and project management. Many of the remaining sections in this chapter are devoted to describing some of these key structures. Structures that describe the product system include, but are not limited to, the requirements tree, system architecture and certain symbolic information such as system drawings, schematics, and data bases. The structures that describe the producing system include the project's work breakdown, schedules, cost accounts and organization.

Shishko, Robert↗

Quantifying Rates of Erosion Using Historical Aerial Photos - 19478

Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. In this work multi-temporal historical aerial images are analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. Scans of photographic negatives were obtained for multiple photographs and several negatives were scanned with a calibrated photogrammetric scanner by Geomni (formerly Blue Skies Consulting) in Belen, New Mexico. Other negatives were provided by USDA or scanned on a high-resolution, large format media scanner in the United States National Archives Cartographic Research Room, College Park, Maryland (1939). An additional set of aerial photographs was obtained from the USGS EROS Data Center web site. For other years, photographic prints were available for aerial photography. Several prints were scanned with a large-format media scanner at 1200 dpi. Ground control points for selected digital images used in the stereoscopic analyses were identified in ArcMap. Ground control points were collected by matching features which could be identified in the aerial photography as well as the orthophotography or shaded relief imagery from the Lidar dataset. Such features included fence line intersections, road intersections, railroad bridges, and structures, depending on the age of the photography. Many old fence lines and abandoned roads remained visible in the shaded relief imagery while they were completely obscured by tree canopy in modern photography. Ground control elevations were sampled from the Lidar-based DEM or directly from point clouds. The image coordinates and matching three-dimensional projected coordinates of the ground control points were exported to IDL where parameters for a three-dimensional direct linear transformation (DLT) were obtained. Built-in singular value decomposition routines were used to solve an overdetermined set of equations for the eleven parameters of the DLT for a particular historical photograph. A solution to the DLT parameters was obtained for each historical aerial photograph in which features of interest were identified using the stereoscope. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. There was inevitable uncertainty associated with the estimation of points from the historical aerial imagery. This includes estimation of changes in locations of features like gully heads through time. In this work a method to estimate the uncertainty associated with measurements of locations of features using the stereoscopic analysis is developed and applied. For features with locations in space that do not change (e.g., road intersections, etc.) the difference in measurements of locations from the Lidar dataset and that of the estimate from the historical aerial image provides an observation of the measurement error. These points used to estimate the measurement error are referred to as control points. For each of the historical aerial images, a different measurement error distribution is estimated using control points. These measurement errors are combined with the estimated distances from the stereoscopic analyses to produce distributions of gully retreat rates and hillslope erosion rates used along with output from Erosion Working Group (EWG) Landscape Evolution Models (LEMs) to provide statistical probability distributions for input into the PPA Model. Specifically, gully head retreat rate and gully widening rate were characterized using statistical probability distributions. The central limit theorem was then applied to develop rates that represent long-term averages that are used to inform a PPA Model built in GoldSim. The importance of these rates with respect to characterization of future risk can then assessed using sensitivity analyses. This provides an approach to determine if subsequent reduction of the uncertainty associated with the probability distributions of erosion rates will correspond to meaningful reduction in the estimation of future risk from the site. Ultimately this is a general method that utilizes historical aerial imagery to estimate rates of erosion with uncertainty. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Parallel Adaptive Mesh Refinement Library

Parallel Adaptive Mesh Refinement Library (PARAMESH) is a package of Fortran 90 subroutines designed to provide a computer programmer with an easy route to extension of (1) a previously written serial code that uses a logically Cartesian structured mesh into (2) a parallel code with adaptive mesh refinement (AMR). Alternatively, in its simplest use, and with minimal effort, PARAMESH can operate as a domain-decomposition tool for users who want to parallelize their serial codes but who do not wish to utilize adaptivity. The package builds a hierarchy of sub-grids to cover the computational domain of a given application program, with spatial resolution varying to satisfy the demands of the application. The sub-grid blocks form the nodes of a tree data structure (a quad-tree in two or an oct-tree in three dimensions). Each grid block has a logically Cartesian mesh. The package supports one-, two- and three-dimensional models.

Mac-Neice, Peter↗

An MPI + $X$ implementation of contact global search using Kokkos

This paper describes an approach that seeks to parallelize the spatial search associated with computational contact mechanics. In contact mechanics, the purpose of the spatial search is to find “nearest neighbors,” which is the prelude to an imprinting search that resolves the interactions between the external surfaces of contacting bodies. In particular, we are interested in the contact global search portion of the spatial search associated with this operation on domain-decomposition-based meshes. Specifically, we describe an implementation that combines standard domain-decomposition-based MPI-parallel spatial search with thread-level parallelism (MPI-X) available on advanced computer architectures (those with GPU coprocessors). Our goal is to demonstrate the efficacy of the MPI-X paradigm in the overall contact search. Standard MPI-parallel implementations typically use a domain decomposition of the external surfaces of bodies within the domain in an attempt to efficiently distribute computational work. This decomposition may or may not be the same as the volume decomposition associated with the host physics. The parallel contact global search phase is then employed to find and distribute surface entities (nodes and faces) that are needed to compute contact constraints between entities owned by different MPI ranks without further inter-rank communication. Key steps of the contact global search include computing bounding boxes, building surface entity (node and face) search trees and finding and distributing entities required to complete on-rank (local) spatial searches. To enable source-code portability and performance across a variety of different computer architectures, we implemented the algorithm using the Kokkos hardware abstraction library. While we targeted development towards machines with a GPU accelerator per MPI rank, we also report performance results for OpenMP with a conventional multi-core compute node per rank. Results here demonstrate a 47 % decrease in the time spent within the global search algorithm, comparing the reference ACME algorithm with the GPU implementation, on an 18M face problem using four MPI ranks. As a result, while further work remains to maximize performance on the GPU, this result illustrates the potential of the proposed implementation.

97 MATHEMATICS AND COMPUTING↗

Multistage Stochastic optimization for mid-term integrated generation and maintenance scheduling of cascaded hydroelectric system with renewable energy uncertainty

The uncertainties resulting from the escalating penetration of renewable energy resources pose severe challenges to the efficient operation of modern power systems. Hydroelectricity is characterized by its flexibility, controllability, and reliability, and thus becomes one of the most ideal energy resources to hedge against such uncertainties. This paper studies the mid-term integrated generation and maintenance scheduling of a cascaded hydroelectric system (CHS) consisting of multiple cascaded reservoirs and hydroelectric units. To precisely describe the mid-term water regulation policies, the hydraulic coupling relationship and water-energy nexus of CHS are incorporated into the proposed optimization model. The uncertainties of natural water inflow and the power outputs of wind/solar energy generation are taken into consideration and captured via a stochastic process modeled by a scenario tree. A multistage stochastic optimization (MSO) approach is developed to coordinate the complementary operations of multiple energy resources, by optimizing the mid-term water resource management, generation scheduling, and maintenance scheduling of CHS. The proposed MSO model is formulated as a large-scale mixed-integer linear program that presents significant computational intractability. To address this issue, a tailored Benders decomposition algorithm is developed. Two real-world case studies are conducted to demonstrate the capability and characteristics of the proposed model and algorithm. The computational results show that the proposed MSO model can exploit the flexibility of hydroelectricity to efficiently respond to variable wind and solar power, and reserve water resources for the generation in peak months to reduce the consumption of fossil fuel. Furthermore, the proposed solution approach also exhibits promising computational efficiency when handling large-scale models.

13 HYDRO ENERGY↗

Functional activity and functional gene diversity of a Cu-contaminated soil remediated by aided phytostabilization using compost, dolomitic limestone and a mixed tree stand

Trace elements (TEs) availability, biochemical activity and functional gene diversity was studied in a Cu-contaminated soil, revegetated for six years with a mixed stand of willow, black poplar, and false indigo-bush, and amended or not with compost plus dolomitic limestone (OMDL). The OMDL amendment significantly reduced Cu and As availability and soil toxicity, and increased the biochemical activity and microbial functional diversity assessed with the GEOCHIP technique, as compared to the unamended soil (Unt). The OMDL soil showed significantly higher abundance of 25 functional genes involved in decomposition organic compounds, and 11, 3 and 11 functional genes involved in the N, P and S biogeochemical cycles. Functional gene abundance was positively correlated with nutrient contents but negatively correlated with Cu availability and soil toxicity. The abundance of microbial functional genes encoding for resistance to various TEs also increased, possibly due to the microbial proliferation and lower Cu exposure in the presence of high total soil Cu concentration. Genes encoding for antibiotic resistance due to the co-occurrence of TEs and antibiotic resistant genes on genetic mobile elements. Finally, overall, phytomanagement confirmed its potential to restore the biological fertility and diversity of a severely Cu-contaminated soil, but the increase of TEs and antibiotic resistant gene abundances deserve attention in future studies.

54 ENVIRONMENTAL SCIENCES↗

Nitrogen and phosphorus cycling in an ombrotrophic peatland: a benchmark for assessing change

Aims Slow decomposition and isolation from groundwater mean that ombrotrophic peatlands store a large amount of soil carbon (C) but have low availability of nitrogen (N) and phosphorus (P). To better understand the role these limiting nutrients play in determining the C balance of peatland ecosystems, we compile comprehensive N and P budgets for a forested bog in northern Minnesota, USA. Methods N and P within plants, soils, and water are quantified based on field measurements. The resulting empirical dataset are then compared to modern-day, site-level simulations from the peatland land surface version of the Energy Exascale Earth System Model (ELM-SPRUCE).

Salmon, Verity G.↗

Sparse Gaussian elimination with controlled fill-in on a shared memory multiprocessor

It is shown that in sparse matrices arising from electronic circuits, it is possible to do computations on many diagonal elements simultaneously. A technique for obtaining an ordered compatible set directly from the ordered incompatible table is given. The ordering is based on the Markowitz number of the pivot candidates. This technique generates a set of compatible pivots with the property of generating few fills. A novel heuristic algorithm is presented that combines the idea of an order-compatible set with a limited binary tree search to generate several sets of compatible pivots in linear time. An elimination set for reducing the matrix is generated and selected on the basis of a minimum Markowitz sum number. The parallel pivoting technique presented is a stepwise algorithm and can be applied to any submatrix of the original matrix. Thus, it is not a preordering of the sparse matrix and is applied dynamically as the decomposition proceeds. Parameters are suggested to obtain a balance between parallelism and fill-ins. Results of applying the proposed algorithms on several large application matrices using the HEP multiprocessor (Kowalik, 1985) are presented and analyzed.

Alaghband, Gita↗

Alternative regularizations for Outer-Approximation algorithms for convex MINLP

In this work, we extend the regularization framework from Kronqvist et al. (Math Program 180(1):285–310, 2020) by incorporating several new regularization functions and develop a regularized single-tree search method for solving convex mixed-integer nonlinear programming (MINLP) problems. We propose a set of regularization functions based on distance metrics and Lagrangean approximations, used in the projection problem for finding new integer combinations to be used within the Outer-Approximation (OA) method. The new approach, called Regularized Outer-Approximation (ROA), has been implemented as part of the open-source Mixed-integer nonlinear decomposition toolbox for Pyomo—MindtPy. We compare the OA method with seven regularization function alternatives for ROA. Moreover, we extend the LP/NLP Branch and Bound method proposed by Quesada and Grossmann (Comput Chem Eng 16(10–11):937–947, 1992) to include regularization in an algorithm denoted RLP/NLP. We provide convergence guarantees for both ROA and RLP/NLP. Finally, we perform an extensive computational experiment considering all convex MINLP problems in the benchmark library MINLPLib. The computational results show clear advantages of using regularization combined with the OA method.

Convex Mixed-integer nonlinear programming↗

Toward a Predictive Understanding of the Response of Belowground Microbial Carbon Turnover to Climate Change Drivers in a Boreal Peatland

Peatlands are wetlands that act as giant "carbon banks", storing approximately one-third of Earth's soil carbon. Peatland carbon banks cover large areas in northern latitudes including the Arctic and evidence from previous work suggests that northern peatlands store a lot of carbon partly because they are cold and microbes don't work very fast to decompose organic matter in the cold. The concern is that with increasing temperatures, microbes will be become more active and release more carbon from the "bank" as greenhouse gases, carbon dioxide and natural gas (methane). However, there are many more factors that could determine why peatlands store so much carbon such as the fact that they are acidic, low in oxygen, or because they contain a peat moss carpet. All of these factors have been shown to slow down the activity of microbes that release greenhouse gases through the decomposition of organic matter. For example, peat mosses naturally produce organic chemical compounds that are thought to poison microbes, thereby limiting their activity. For this reason, peat moss was used in the past as a food preservative. As the amount of carbon dioxide increases in Earth's atmosphere, it also predicted that plants will become more active and produce more usable organic matter for microbes. Therefore, this project addressed the resilience of the peatland carbon bank to warming and the accumulation carbon dioxide in the atmosphere. The project was conducted at the Marcell Experimental Forest (MEF), northern Minnesota, where the Oak Ridge National Lab (ORNL) has established an experimental site known as Spruce and Peatland Response Under Climatic and Environmental Change (SPRUCE). The SPRUCE site simulates environmental change at the ecosystem scale, allowing us to understand how all organisms, from microbes to trees, will respond to higher temperatures and increased carbon dioxide in the atmosphere. From 2014 to 2017, the response of peatlands to environmental change was evaluated through multiple field campaigns each year along with laboratory research. Environmental changes were initiated in stages; soil heating began in 2014, air heating in 2015, and carbon dioxide enrichment in 2016. It was shown that after 1 year, soil heating alone resulted in an exponential increase in methane release to the atmosphere. However, this response was due solely to surface processes and not degradation of deep peat, which represents much older carbon stores. Studies in the laboratory showed that only the top 20–30 cm of peat from experimental plots had higher methane production rates at elevated temperatures. Further investigations using the latest advanced analytical chemistry and microbiology methods confirmed that greenhouse gas production primarily originated from decomposition of organic matter produced recently by plants and not from old peat. No differences in microbial abundances, dissolved organic matter concentrations or degradative enzyme activities were observed among treatments after 1 year of soil heating alone. The activity and not the abundance of microbial communities responded over this time scale. Subsequent field investigations from 2015 to 2017 supported these initial investigations that changes in chemistry and the release of gases from microbial activity was primarily occurring at the surface of the peatland. Thus, our initial results over the first few years of simulated environmental change indicated that the large store of deep peat carbon will be resistant to anaerobic degradation under future climatic warming. Peatland soils are poor in mineral nutrients, and it was thought previously that the decomposition of plant matter in the soil is controlled by processes leading primarily to methane production (methanogenesis), ultimately producing similar amounts of carbon dioxide and methane during organic matter degradation. However, we observed that much more carbon dioxide is produced than methane from the surface peatland. Using a combination of state of the art molecular techniques, this project showed that novel microbial processes involving the respiration or “breathing” of organic matter itself may be responsible for most of gas release from the peatland as carbon dioxide. This has opened up a whole new area of research in how wetlands respond to environmental change.

54 ENVIRONMENTAL SCIENCES↗

Unitarization from geometry

We study the perturbative unitarity of scattering amplitudes in general dimensional reductions of Yang–Mills theories and general relativity on closed internal manifolds. For the tree amplitudes of the dimensionally reduced theory to have the expected high-energy behavior of the higher-dimensional theory, the masses and cubic couplings of the Kaluza–Klein states must satisfy certain sum rules that ensure there are nontrivial cancellations between Feynman diagrams. These sum rules give constraints on the spectra and triple overlap integrals of eigenfunctions of Laplacian operators on the internal manifold and can be proven directly using Hodge and eigenfunction decompositions. One consequence of these constraints is that there is an upper bound on the ratio of consecutive eigenvalues of the scalar Laplacian on closed Ricci-flat manifolds with special holonomy. This gives a sharp bound on the allowed gaps between Kaluza–Klein excitations of the graviton that also applies to Calabi–Yau compactifications of string theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cosmohedra

It has been a long-standing challenge to find a geometric object underlying the cosmological wavefunction for Tr(ϕ 3 ) theory, generalizing associahedra and surfacehedra for scattering amplitudes. In this note, we describe a new class of polytopes — “cosmohedra” — that provide a natural solution to this problem. The faces of associahedra capture the combinatorics of non-overlapping chords of the momentum polygon, reflecting all partial factorizations of amplitudes. Cosmohedra are far richer — instead of non-overlapping chords, their faces capture the “russian doll” structure of non-overlapping subpolygons that determine the wavefunction. We show that cosmohedra are intimately related to associahedra, obtained by “blowing up” faces of the associahedron in a simple way. We give a full combinatorial description of cosmohedron faces and their factorization properties, and provide an explicit realization in terms of facet inequalities that further “shave” the facet inequalities of the associahedron. We also discuss a novel way for computing the wavefunction from cosmohedron geometry that extends the usual connection with polytope canonical forms. We illustrate cosmohedra with examples at tree-level and one loop; the close connection to surfacehedra suggests the generalization to all loop orders. Moving beyond the wavefunction, we briefly describe “cosmological correlahedra” for full correlators, which are one higher-dimensional polytopes, interpolating between associahedra and cosmohedra on opposite facets in an extra direction associated with the total energy. We speculate on how the existence of cosmohedra might suggest a “stringy” formulation for the cosmological wavefunction/correlators, generalizing the way in which the Minkowski sum decomposition of associahedra naturally extend particle to string amplitudes.

Scattering Amplitudes↗

NANO.PTML model for read-across prediction of nanosystems in neurosciences. computational model and experimental case of study

Abstract Neurodegenerative diseases involve progressive neuronal death. Traditional treatments often struggle due to solubility, bioavailability, and crossing the Blood-Brain Barrier (BBB). Nanoparticles (NPs) in biomedical field are garnering growing attention as neurodegenerative disease drugs (NDDs) carrier to the central nervous system. Here, we introduced computational and experimental analysis. In the computational study, a specific IFPTML technique was used, which combined Information Fusion (IF) + Perturbation Theory (PT) + Machine Learning (ML) to select the most promising Nanoparticle Neuronal Disease Drug Delivery (N2D3) systems. For the application of IFPTML model in the nanoscience, NANO.PTML is used. IF-process was carried out between 4403 NDDs assays and 260 cytotoxicity NP assays conducting a dataset of 500,000 cases. The optimal IFPTML was the Decision Tree (DT) algorithm which shown satisfactory performance with specificity values of 96.4% and 96.2%, and sensitivity values of 79.3% and 75.7% in the training (375k/75%) and validation (125k/25%) set. Moreover, the DT model obtained Area Under Receiver Operating Characteristic (AUROC) scores of 0.97 and 0.96 in the training and validation series, highlighting its effectiveness in classification tasks. In the experimental part, two samples of NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) were synthesized by thermal decomposition of an iron(III) oleate (FeOl) precursor and structurally characterized by different methods. Additionally, in order to make the as-synthesized hydrophobic NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) soluble in water the amphiphilic CTAB (Cetyl Trimethyl Ammonium Bromide) molecule was employed. Therefore, to conduct a study with a wider range of NP system variants, an experimental illustrative simulation experiment was performed using the IFPTML-DT model. For this, a set of 500,000 prediction dataset was created. The outcome of this experiment highlighted certain NANO.PTML systems as promising candidates for further investigation. The NANO.PTML approach holds potential to accelerate experimental investigations and offer initial insights into various NP and NDDs compounds, serving as an efficient alternative to time-consuming trial-and-error procedures.

60 APPLIED LIFE SCIENCES↗

Partitioning CO 2 fluxes with isotopologue measurements and modeling to understand mechanisms of forest carbon sequestration

1. Objectives This project combines automated in situ observations of the isotopologues of CO 2 with root observations, novel experimental manipulations of belowground processes, and isotope-enabled ecosystem modeling to investigate mechanisms of below- vs. aboveground carbon sequestration at the Harvard Forest Environmental Measurements Site (EMS). The proposed objectives, which have now been largely accomplished, include: A. Partitioning of net ecosystem CO 2 exchange (NEE) into photosynthesis and respiration using long-term continuous observations of the isotopic composition of NEE, and analysis of their dynamics ; B. Investigation of the influence of vegetation phenology on the timing and magnitude of carbon allocated belowground using measurements of root growth and indices of belowground autotrophic vs. heterotrophic respiration (via trenched plots and isotope measurements); C. Testing whether plant allocation of carbon belowground stimulates the microbial decomposition of soil organic matter, using in situ rhizosphere simulation experiments wherein realistic quantities of artificial isotopically-labeled exudates are released into the soil; and D. Synthesis and interpretation of the above data using the Ecosystem Demography Model 2 (ED2). 2. Highlights Accomplishments: • Our isotopic eddy flux record has completed its 5th full year and has been used to independently estimate ecosystem-scale respiration and photosynthesis. • Soil surface chamber isotopic flux measurements were carried out during three growing seasons, in conjunction with a trenching manipulation. Key findings to date (listed by objective): A. Partitioning of Net Ecosystem Exchange: 1. Ecosystem respiration is lower during the day than at night—the first robust evidence of the inhibition of leaf respiration by light (the “Kok effect”) at the ecosystem scale. 2. Because it neglects the Kok effect, the standard NEE partitioning approach overestimates ecosystem photosynthesis (by ~25%) and daytime respiration (by ~100%) in the first half of the growing season at our site, and portrays ecosystem photosynthetic light-use efficiency as declining when in fact it is stable until autumnal senescence. B. Vegetation Phenology and belowground allocation: Findings: 1. Autotrophic respiration (Ra) showed a seasonal pattern, peaking in mid-summer when trees were most active. 2. The effective age of the substrate for belowground respiration is less than 2 weeks. 3. Above and belowground phenology are more synchronous in deciduous hardwood stands than evergreen hemlock stands. 4. The decline in root respiration rates in the fall is related to temperature rather than acclimation of root respiration or substrate limitations. Methodological Issues: 5. The isotopic signatures of autotrophic and heterotrophic respiration are too similar for isotopic partitioning of belowground respiration into these two components at our site—in keeping with the recent findings of Bowling et al. (2015) in a subalpine conifer forest. 6. Artifacts of the trenching method, such as changes in soil moisture and increased carbon substrate from the newly severed roots, are significant and need to be quantified when determining daily to annual estimates of autotrophic and heterotrophic respiration. C. Effects of simulated exudates on priming of microbial decomposition: The stoichiometry of root exudates influences both the amount and the mechanism by which priming occurs. At low C:N, SOC loss is caused by an increase in microbial efficiency. At high C:N, SOC loss is caused by an increase in microbial biomass. D. Modeling with the Ecosystem Demography Model (ED2): 1. Incorporation of 13C tracking to create an isotopically-enabled Ecosystem Demography v2 model (ED2) 2. State-of-the-art parameter optimization methodology developed for improving ED2 model predictions and parameters. 3. Significantly improved model predictions of growth- and maintenance-related carbon fluxes and 13C fluxes

54 ENVIRONMENTAL SCIENCES↗

Partitioning CO2 fluxes with isotopologue measurements and modeling to understand mechanisms of forest carbon sequestration

1. Objectives This project combines automated in situ observations of the isotopologues of CO2 with root observations, novel experimental manipulations of belowground processes, and isotope-enabled ecosystem modeling to investigate mechanisms of below- vs. aboveground carbon sequestration at the Harvard Forest Environmental Measurements Site (EMS). The proposed objectives, which have now been largely accomplished, include: A. Partitioning of net ecosystem CO2 exchange (NEE) into photosynthesis and respiration using long-term continuous observations of the isotopic composition of NEE, and analysis of their dynamics ; B. Investigation of the influence of vegetation phenology on the timing and magnitude of carbon allocated belowground using measurements of root growth and indices of belowground autotrophic vs. heterotrophic respiration (via trenched plots and isotope measurements); C. Testing whether plant allocation of carbon belowground stimulates the microbial decomposition of soil organic matter, using in situ rhizosphere simulation experiments wherein realistic quantities of artificial isotopically-labeled exudates are released into the soil; and D. Synthesis and interpretation of the above data using the Ecosystem Demography Model 2 (ED2). 2. Highlights Accomplishments: • Our isotopic eddy flux record has completed its 5th full year and has been used to independently estimate ecosystem-scale respiration and photosynthesis. • Soil surface chamber isotopic flux measurements were carried out during three growing seasons, in conjunction with a trenching manipulation. Key findings to date (listed by objective): A. Partitioning of Net Ecosystem Exchange: 1. Ecosystem respiration is lower during the day than at night—the first robust evidence of the inhibition of leaf respiration by light (the “Kok effect”) at the ecosystem scale. 2. Because it neglects the Kok effect, the standard NEE partitioning approach overestimates ecosystem photosynthesis (by ~25%) and daytime respiration (by ~100%) in the first half of the growing season at our site, and portrays ecosystem photosynthetic light-use efficiency as declining when in fact it is stable until autumnal senescence. B. Vegetation Phenology and belowground allocation: Findings: 1. Autotrophic respiration (Ra) showed a seasonal pattern, peaking in mid-summer when trees were most active. 2. The effective age of the substrate for belowground respiration is less than 2 weeks. 3. Above and belowground phenology are more synchronous in deciduous hardwood stands than evergreen hemlock stands. 4. The decline in root respiration rates in the fall is related to temperature rather than acclimation of root respiration or substrate limitations. Methodological Issues: 5. The isotopic signatures of autotrophic and heterotrophic respiration are too similar for isotopic partitioning of belowground respiration into these two components at our site—in keeping with the recent findings of Bowling et al. (2015) in a subalpine conifer forest. 6. Artifacts of the trenching method, such as changes in soil moisture and increased carbon substrate from the newly severed roots, are significant and need to be quantified when determining daily to annual estimates of autotrophic and heterotrophic respiration. C. Effects of simulated exudates on priming of microbial decomposition: The stoichiometry of root exudates influences both the amount and the mechanism by which priming occurs. At low C:N, SOC loss is caused by an increase in microbial efficiency. At high C:N, SOC loss is caused by an increase in microbial biomass. D. Modeling with the Ecosystem Demography Model (ED2): 1. Incorporation of 13C tracking to create an isotopically-enabled Ecosystem Demography v2 model (ED2) 2. State-of-the-art parameter optimization methodology developed for improving ED2 model predictions and parameters. 3. Significantly improved model predictions of growth- and maintenance-related carbon fluxes and 13C fluxes

54 ENVIRONMENTAL SCIENCES↗

Low-intensity frequent fires in coniferous forests transform soil organic matter in ways that may offset ecosystem carbon losses

The impact of shifting disturbance regimes on soil carbon (C) storage is a key uncertainty in global change research. Wildfires in coniferous forests are becoming more frequent in many regions, potentially causing large C emissions. Repeated low-intensity prescribed fires can mitigate wildfire severity, but repeated combustion may decrease soil C unless compensatory responses stabilize soil organic matter. Here, we tested how 30 years of decadal prescribed burning affected C and nitrogen (N) in plants, detritus, and soils in coniferous forests in the Sierra Nevada mountains, USA. Tree basal area and litter stocks were resilient to fire, but fire reduced forest floor C by 77% (-36.4 Mg C/ha). In mineral soils, fire reduced C that was free from minerals by 41% (-4.4 Mg C/ha) but not C associated with minerals, and only in depths ≤ 5 cm. Fire also transformed the properties of remaining mineral soil organic matter by increasing the proportion of C in a pyrogenic form (from 3.2% to 7.5%) and associated with minerals (from 46% to 58%), suggesting the remaining soil C is more resistant to decomposition. Laboratory assays illustrated that fire reduced microbial CO2 respiration rates by 55% and the activity of eight extracellular enzymes that degrade cellulosic and aromatic compounds by 40–66%. Lower decomposition was correlated with lower inorganic N (-49%), especially ammonium, suggesting N availability is coupled with decomposition. The relative increase in forms of soil organic matter that are resistant to decay or stabilized onto mineral surfaces, and the associated decline in decomposition suggest that low-intensity fires may promote mineral soil C storage in pools with long mean residence times in coniferous forests.

54 ENVIRONMENTAL SCIENCES↗

A Method to Access Absolute fIPAR fo Vegetation in Spatially Complex Ecosystems

Arid and semi-arid lands compose a large fraction of the earth's terrestrial vegetation, and thereby contribute significantly to global atmospheric-biospheric interactions. The thorny shrubs and small trees in these semi-arid shrub lands have counterparts throughout much of the world's tropical and subtropical zones and have captured substantial areas of the world's former grasslands. The objective of our field and remotely sensed measurements in the semi-arid shrublands of Texas is to monitor interannual variability and directional change in landscape structure, ecosystem processes and atmosphere-biosphere exchanges. To understand the role ecosystems play in controlling the composition of the atmosphere, it is necessary to quantify processes such as photosynthesis and primary production, decomposition and soil carbon storage, and trace gas exchanges. Photosynthesis is the link whereby surface-atmosphere exchanges such as the radiation balance and exchange of heat, moisture, and gas can be inferred. It also describes the efficiency of carbon dioxide exchange and is directly related to the primary production of vegetation. Our efforts in this paper focus on the indirect, quantification of photosynthesis, and thereby carbon flux and net primary production, via remote sensing and direct measurements of intercepted photosynthetically active radiation (IPAR).

Wessman, Carol A.↗