Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ECOSYSTEM”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Ecosystem carbon transit versus turnover times in response to climate warming and rising atmospheric CO 2 concentration

Ecosystem carbon (C) transit time is a critical diagnostic parameter to characterize land C sequestration. This parameter has different variants in the literature, including a commonly used turnover time. However, we know little about how different transit time and turnover time are in representing carbon cycling through multiple compartments under a non-steady state. In this study, we estimate both C turnover time as defined by the conventional stock over flux and mean C transit time as defined by the mean age of C mass leaving the system. We incorporate them into the Community Atmosphere Biosphere Land Exchange (CABLE) model to estimate C turnover time and transit time in response to climate warming and rising atmospheric [CO 2 ]. Modelling analysis shows that both C turnover time and transit time increase with climate warming but decrease with rising atmospheric [CO 2 ]. Warming increases C turnover time by 2.4 years and transit time by 11.8 years in 2100 relative to that at steady state in 1901. During the same period, rising atmospheric [CO 2 ] decreases C turnover time by 3.8 years and transit time by 5.5 years. Our analysis shows that 65 % of the increase in global mean C transit time with climate warming results from the depletion of fast-turnover C pool. The remaining 35 % increase results from accompanied changes in compartment C age structures. Similarly, the decrease in mean C transit time with rising atmospheric [CO 2 ] results approximately equally from replenishment of C into fast-turnover C pool and subsequent decrease in compartment C age structure. Greatly different from the transit time, the turnover time, which does not account for changes in either C age structure or composition of respired C, underestimated impacts of warming and rising atmospheric [CO 2 ] on C diagnostic time and potentially led to deviations in estimating land C sequestration in multi-compartmental ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microbial dormancy and its impacts on northern temperate and boreal terrestrial ecosystem carbon budget

A large amount of soil carbon in northern temperate and boreal regions could be emitted as greenhouse gases in a warming future. However, lacking detailed microbial processes such as microbial dormancy in current biogeochemistry models might have biased the quantification of the regional carbon dynamics. Here the effect of microbial dormancy was incorporated into a biogeochemistry model to improve the quantification for the last century and this century. Compared with the previous model without considering the microbial dormancy, the new model estimated the regional soils stored 75.9 Pg more C in the terrestrial ecosystems during the last century and will store 50.4 and 125.2 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, in this century. This study highlights the importance of the representation of microbial dormancy in earth system models to adequately quantify the carbon dynamics in the northern temperate and boreal natural terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Biogeochemical and plant trait mechanisms drive enhanced methane emissions in response to whole-ecosystem warming

Abstract. Climate warming perturbs ecosystem carbon (C) cycling, causing both positive and negative feedbacks on greenhouse gas emissions. In 2016, we began a tidal marsh field experiment in two vegetation communities to investigate the mechanisms by which whole-ecosystem warming alters C gain, via plant-driven sequestration in soils, and C loss, primarily via methane (CH4) emissions. Here, we report the results from the first 4 years. As expected, warming of 5.1 ∘C more than doubled CH4 emissions in both plant communities. We propose this was caused by a combination of four mechanisms: (i) a decrease in the proportion of CH4 consumed by CH4 oxidation, (ii) more C substrates available for methanogenesis, (iii) reduced competition between methanogens and sulfate-reducing bacteria, and (iv) indirect effects of plant traits. Plots dominated by Spartina patens consistently emitted more CH4 than plots dominated by Schoenoplectus americanus, indicating key differences in the roles these common wetland plants play in affecting anaerobic soil biogeochemistry and suggesting that plant composition can modulate coastal wetland responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Extending a land-surface model with <i>Sphagnum</i> moss to simulate responses of a northern temperate bog to whole ecosystem warming and elevated CO<sub>2</sub>

Abstract. Mosses need to be incorporated into Earth system models to better simulate peatland functional dynamics under the changing environment. Sphagnum mosses are strong determinants of nutrient, carbon, and water cycling in peatland ecosystems. However, most land-surface models do not include Sphagnum or other mosses as represented plant functional types (PFTs), thereby limiting predictive assessment of peatland responses to environmental change. In this study, we introduce a moss PFT into the land model component (ELM) of the Energy Exascale Earth System Model (E3SM) by developing water content dynamics and nonvascular photosynthetic processes for moss. The model was parameterized and independently evaluated against observations from an ombrotrophic forested bog as part of the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project. The inclusion of a Sphagnum PFT with some Sphagnum-specific processes in ELM allows it to capture the observed seasonal dynamics of Sphagnum gross primary production (GPP) albeit with an underestimate of peak GPP. The model simulated a reasonable annual net primary production (NPP) for moss but with less interannual variation than observed, and it reproduced aboveground biomass for tree PFTs and stem biomass for shrubs. Different species showed highly variable warming responses under both ambient and elevated atmospheric CO2 concentrations, and elevated CO2 altered the warming response direction for the peatland ecosystem. Microtopography is critical: Sphagnum mosses on hummocks and hollows were simulated to show opposite warming responses (NPP decreasing with warming on hummocks but increasing in hollows), and hummock Sphagnum was modeled to have a strong dependence on water table height. The inclusion of this new moss PFT in global ELM simulations may provide a useful foundation for the investigation of northern peatland carbon exchange, enhancing the predictive capacity of carbon dynamics across the regional and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quantifying the role of moss in terrestrial ecosystem carbon dynamics in northern high latitudes

Mosses are ubiquitous in northern terrestrial ecosystems, and play an important role in regional carbon, water and energy cycling. Current global land surface models that do not consider mosses may bias the quantification of regional carbon dynamics. Here we incorporate mosses as a new plant functional type into the process-based Terrestrial Ecosystem Model (TEM 5.0), to develop a new model (TEM_Moss). The new model explicitly quantifies the interactions between vascular plants and mosses and their competition for energy, water, and nutrients. Compared to the estimates using TEM 5.0, the new model estimates that the regional terrestrial soils currently store 132.7 Pg more C and will store 157.5 and 179.1 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, by the end of the 21st century. Ensemble regional simulations forced with different parameters for the 21st century with TEM_Moss predict that the region will accumulate 161.1±142.1 Pg C under the RCP2.6 scenario and 186.7±166.1 Pg C under the RCP8.5 scenario over the century. Our study highlights the necessity of coupling moss into Earth system models to adequately quantify terrestrial carbon–climate feedbacks in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River basin

Wildfires impact vegetation mortality and productivity and are increasing in intensity, frequency, and spatial area in the western United States. The rates of vegetation recovery after fires play a major role in the reestablishment of biomass and ecosystem functioning (e.g., structure, resilience, and productivity), but such recovery rates are poorly understood. Here we use remotely sensed data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to quantify the resistance and resilience of leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET) to 138 wildfires of various burn severity across the Columbia River basin (CRB) of the Pacific Northwest in 2015. Increasing burn severity caused lower resistance and resilience for all three variables. Resistance and resilience are highest in grasslands, intermediate in savanna, and lowest in needleleaf evergreen forests, consistent with the adaptation of these vegetation types to fire. LAI has consistently lower resistance and resilience than GPP and ET, which is consistent with physical and physiological mechanisms that compensate for reduced LAI. Resilience is influenced by precipitation, vapor pressure deficit (VPD), and burn severity across all three vegetation types; however, burn severity plays a more minor role in grasslands. Increasing wildfire severity will reduce the resistance and resilience and lengthen the recovery time of vegetation structure and fluxes with climate change, with significant consequences for the provision of ecosystem functioning and implications for model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reviews and syntheses: Variable inundation across Earth's terrestrial ecosystems

Abstract. The structure, function, and dynamics of Earth's terrestrial ecosystems are profoundly influenced by how often (frequency) and how long (duration) they are inundated with water. A diverse array of natural and human-engineered systems experience temporally variable inundation whereby they fluctuate between inundated and non-inundated states. Variable inundation spans extreme events to predictable sub-daily cycles. Variably inundated ecosystems (VIEs) include hillslopes, non-perennial streams, wetlands, floodplains, temporary ponds, tidal systems, storm-impacted coastal zones, and human-engineered systems. VIEs are diverse in terms of inundation regimes, water chemistry and flow velocity, soil and sediment properties, vegetation, and many other properties. The spatial and temporal scales of variable inundation are vast, ranging from sub-meter to whole landscapes and from sub-hourly to multi-decadal. The broad range of system types and scales makes it challenging to predict the hydrology, biogeochemistry, ecology, and physical evolution of VIEs. Despite all experiencing the loss and gain of an overlying water column, VIEs are rarely considered together in conceptual, theoretical, modeling, or measurement frameworks and approaches. Studying VIEs together has the potential to generate mechanistic understanding that is transferable across a much broader range of environmental conditions, relative to knowledge generated by studying any one VIE type. We postulate that enhanced transferability will be important for predicting changes in VIE function in response to global change. Here we aim to catalyze cross-VIE science that studies drivers and impacts of variable inundation across Earth's VIEs. To this end, we complement expert mini-reviews of eight major VIE systems with overviews of VIE-relevant methods and challenges associated with scale. We conclude with perspectives on how cross-VIE science can derive transferable understanding via unifying conceptual models in which the impacts of variable inundation are studied across multi-dimensional environmental space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production

Terrestrial photosynthesis is the basis for vegetation growth and drives the land carbon cycle. Accurately simulating gross primary production (GPP, ecosystem-level apparent photosynthesis) is key for satellite monitoring and Earth system model predictions under climate change. While robust models exist for describing leaf-level photosynthesis, predictions diverge due to uncertain photosynthetic traits and parameters which vary on multiple spatial and temporal scales. Here, we describe and evaluate a GPP (photosynthesis per unit ground area) model, the P-model, that combines the Farquhar–von Caemmerer–Berry model for C 3 photosynthesis with an optimality principle for the carbon assimilation–transpiration trade-off, and predicts a multi-day average light use efficiency (LUE) for any climate and C 3 vegetation type. The model builds on the theory developed in Prentice et al. (2014) and Wang et al. (2017a) and is extended to include low temperature effects on the intrinsic quantum yield and an empirical soil moisture stress factor. The model is forced with site-level data of the fraction of absorbed photosynthetically active radiation (fAPAR) and meteorological data and is evaluated against GPP estimates from a globally distributed network of ecosystem flux measurements. Although the P-model requires relatively few inputs, the R 2 for predicted versus observed GPP based on the full model setup is 0.75 (8 d mean, 126 sites) – similar to comparable satellite-data-driven GPP models but without predefined vegetation-type-specific parameters. The R 2 is reduced to 0.70 when not accounting for the reduction in quantum yield at low temperatures and effects of low soil moisture on LUE. The R 2 for the P-model-predicted LUE is 0.32 (means by site) and 0.48 (means by vegetation type). Applying this model for global-scale simulations yields a total global GPP of 106–122 Pg C yr –1 (mean of 2001–2011), depending on the fAPAR forcing data. The P-model provides a simple but powerful method for predicting – rather than prescribing – light use efficiency and simulating terrestrial photosynthesis across a wide range of conditions. The model is available as an R package (rpmodel).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Implementation and evaluation of the unified stomatal optimization approach in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES)

Abstract. Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Global sensitivity analysis, probabilistic calibration, and predictive assessment for the data assimilation linked ecosystem carbon model

In this paper we propose a probabilistic framework for an uncertainty quantification (UQ) study of a carbon cycle model and focus on the comparison between steady-state and transient simulation setups. A global sensitivity analysis (GSA) study indicates the parameters and parameter couplings that are important at different times of the year for quantities of interest (QoIs) obtained with the data assimilation linked ecosystem carbon (DALEC) model. We then employ a Bayesian approach and a statistical model error term to calibrate the parameters of DALEC using net ecosystem exchange (NEE) observations at the Harvard Forest site. The calibration results are employed in the second part of the paper to assess the predictive skill of the model via posterior predictive checks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Sample Identifiers and Metadata to Support Data Management and Reuse in Multidisciplinary Ecosystem Sciences

Physical samples are foundational entities for research across biological, Earth, and environmental sciences. Data generated from sample-based analyses are not only the basis of individual studies, but can also be integrated with other data to answer new and broader-scale questions. Ecosystem studies increasingly rely on multidisciplinary team-science to study climate and environmental changes. While there are widely adopted conventions within certain domains to describe sample data, these have gaps when applied in a multidisciplinary context. In this study, we reviewed existing practices for identifying, characterizing, and linking related environmental samples. We then tested practicalities of assigning persistent identifiers to samples, with standardized metadata, in a pilot field test involving eight United States Department of Energy projects. Participants collected a variety of sample types, with analyses conducted across multiple facilities. We address terminology gaps for multidisciplinary research and make recommendations for assigning identifiers and metadata that supports sample tracking, integration, and reuse. Furthermore, our goal is to provide a practical approach to sample management, geared towards ecosystem scientists who contribute and reuse sample data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Active Layer Hydrology in an Arctic Tundra Ecosystem: Quantifying Water Sources and Cycling Using Water Stable Isotopes: Supporting Data

Data include results from water isotope analyses (one *.csv file) for samples collected in Utqiagvik (Barrow), Alaska during August and September 2012. Samples were from surface and soil pore waters from 17 drainages that could be interlake (basins with polygonal terrain), different-aged drain thaw lake basins (young, medium, old, or ancient), or a combination of different aged basins. Samples taken in different drainage flow types at three different depths at each location in and around the Barrow Environmental Observatory. Precipitation stable isotope data are also included (added in October 2019 with no changes to previously released data). This dataset used in Throckmorton, et.al. 2016.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Radiocarbon in Ecosystem Respiration and Soil Pore-Space CO2 with Surface Gas Flux, Air Temperature, and Soil Temperature and Moisture, Utqiagvik (Barrow), Alaska, 2012-2014

Dataset includes Delta14C measurements made from CO2 that was collected and purified in 2012-2014 from surface soil chambers, soil pore space, and background atmosphere. In addition to 14CO2 data, dataset includes co-located measurements of CO2 and CH4 flux, soil and air temperature, and soil moisture. Measurements and field samples were taken from intensive study site 1 areas A, B, and C, and the site 0 and AB transects, from specified positions in high-centered, flat-centered, and low centered polygons. The column "sample_size" was added to "radiocarbon_field_Barrow_2012_2013_2014.csv" September 2018 with no changes to other data. This dataset contains one pdf and four csv files.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Baltimore Ecosystem Study: Stream chemistry for core sites in Gwynns Falls

In the Baltimore urban long-term ecological research (LTER) project, (Baltimore Ecosystem Study, BES) we use the watershed approach to evaluate integrated ecosystem function. The LTER research is centered on the Gwynns Falls watershed, a 17,150 ha catchment that traverses a gradient from the urban core of Baltimore, through older urban residential (1900 - 1950) and suburban (1950- 1980) zones, rapidly suburbanizing areas and a rural/suburban fringe. Our long-term sampling network includes four longitudinal sampling sites along the Gwynns Falls as well as several small (40 - 100 ha) watersheds located within or near to the Gwynns Falls. The longitudinal sites provide data on water and nutrient fluxes in the different land use zones of the watershed (rural/suburban, rapidly suburbanizing, old suburban, urban core) and the small watersheds provide more focused data on specific land use areas (forest, agriculture, rural/suburban, urban). Each of the gaging sites is continuously monitored for discharge and is sampled weekly for chemistry. Additional chemical sampling is carried out in a supplemental set of sites to provide a greater range of land use. Weekly analyses includes nitrate, phosphate, total nitrogen, total phosphorus, chloride and sulfate, turbidity, fecal coliforms, temperature, dissolved oxygen and pH. Cations, dissolved organic carbon and nitrogen and metals are measured on selected samples. Streamflow data for this site are posted at: http://waterdata.usgs.gov/md/nwis/nwisman?site_no=015835701

Groffman, Peter M (ORCID:0000000183716255)↗
NASA NTRS1974

A Simulation Model for Studying Effects of Pollution and Freshwater Inflow on Secondary Productivity in an Ecosystem

A mathematical model of an ecosystem is developed. Secondary productivity is evaluated in terms of man related and controllable factors. Information from an existing physical parameters model is used as well as pertinent biological measurements. Predictive information of value to estuarine management is presented. Biological, chemical, and physical parameters measured in order to develop models of ecosystems are identified.

Johnson, R. W.↗
NASA NTRS1976

Assessment of the impact of increased solar ultraviolet radiation upon marine ecosystems

Data was provided to assess the potential impact upon marine ecosystems if space shuttle operations contribute to a reduction of the stratospheric ozone layer. The potential for irreversible damage to the productivity, structure and/or functioning of a model estuarine ecosystem by increased UV-B radiation was established. The sensitivity of key community components (the primary producers) to increased UV-B radiation was delineated.

Vandyke, H.↗
NASA NTRS1977

Estimation of old field ecosystem biomass using low altitude imagery

Color-infrared photography was used to evaluate the biomass of experimental plots in an old-field ecosystem that was treated with different levels of waste water from a sewage treatment facility. Cibachrome prints at a scale of approximately 1:1,600 produced from 35 mm color infrared slides were used to analyze density patterns using prepared tonal density scales and multicell grids registered to ground panels shown on the photograph. Correlation analyses between tonal density and vegetation biomass obtained from ground samples and harvests were carried out. Correlations between mean tonal density and harvest biomass data gave consistently high coefficients ranging from 0.530 to 0.896 at the 0.001 significance level. Corresponding multiple regression analysis resulted in higher correlation coefficients. The results of this study indicate that aerial infrared photography can be used to estimate standing crop biomass on waste water irrigated old field ecosystems. Combined with minimal ground truth data, this technique could enable managers of wastewater irrigation projects to precisely time harvest of such systems for maximal removal of nutrients in harvested biomass.

Nor, S. M.↗