Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ecosystem model”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Boosts in leaf-level photosynthetic capacity aid Pinus ponderosa recovery from wildfire

Forests mitigate climate change by sequestering massive amounts of carbon, but recent increases in wildfire activity are threatening carbon storage. Currently, our understanding of wildfire impacts on forest resilience and the mechanisms controlling post-fire recovery remains unresolved due to a lack of empirical data on mature trees in natural settings. Here, we quantify the physiological mechanisms controlling carbon uptake immediately following wildfire in mature individuals of ponderosa pine (Pinus ponderosa), a wide-spread and canopy-dominant tree species in fire-prone forests. While photosynthetic capacity was lower in burned than unburned trees due to an overall depletion of resources, we show that within the burned trees, photosynthetic capacity increases with the severity of damage. Our data reveal that boosts in the efficiency of carbon uptake at the leaf-level may compensate for whole-tree damage, including the loss of leaf area and roots. We further show that heightened photosynthetic capacity in remaining needles on burned trees may be linked with reduced water stress and leaf nitrogen content, providing pivotal information about post-fire physiological processes. Our results have implications for Earth system modeling efforts because measurements of species-level physiological parameters are used in models to predict ecosystem and landscape-level carbon trajectories. Finally, current land management practices do not account for physiological resilience and recovery of severely burned trees. Our results suggest premature harvest may remove individuals that may otherwise survive, irrevocably altering forest carbon balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Preserving isohydricity: vertical environmental variability explains Amazon forest water-use strategies

Abstract Increases in hydrological extremes, including drought, are expected for Amazon forests. A fundamental challenge for predicting forest responses lies in identifying ecological strategies which underlie such responses. Characterization of species-specific hydraulic strategies for regulating water-use, thought to be arrayed along an ‘isohydric–anisohydric’ spectrum, is a widely used approach. However, recent studies have questioned the usefulness of this classification scheme, because its metrics are strongly influenced by environments, and hence can lead to divergent classifications even within the same species. Here, we propose an alternative approach positing that individual hydraulic regulation strategies emerge from the interaction of environments with traits. Specifically, we hypothesize that the vertical forest profile represents a key gradient in drought-related environments (atmospheric vapor pressure deficit, soil water availability) that drives divergent tree water-use strategies for coordinated regulation of stomatal conductance (gs) and leaf water potentials (ΨL) with tree rooting depth, a proxy for water availability. Testing this hypothesis in a seasonal eastern Amazon forest in Brazil, we found that hydraulic strategies indeed depend on height-associated environments. Upper canopy trees, experiencing high vapor pressure deficit (VPD), but stable soil water access through deep rooting, exhibited isohydric strategies, defined by little seasonal change in the diurnal pattern of gs and steady seasonal minimum ΨL. In contrast, understory trees, exposed to less variable VPD but highly variable soil water availability, exhibited anisohydric strategies, with fluctuations in diurnal gs that increased in the dry season along with increasing variation in ΨL. Our finding that canopy height structures the coordination between drought-related environmental stressors and hydraulic traits provides a basis for preserving the applicability of the isohydric-to-anisohydric spectrum, which we show here may consistently emerge from environmental context. Our work highlights the importance of understanding how environmental heterogeneity structures forest responses to climate change, providing a mechanistic basis for improving models of tropical ecosystems.

Forestry↗
DOE OSTI2024

Tropical forests and global change: biogeochemical responses and opportunities for cross-site comparisons, an organized INSPIRE session at the 108th Annual Meeting, Ecological Society of America, Portland, Oregon, USA, August 2023

In this study, tropical forests play a critical role in the global carbon (C) cycle. These ecosystems maintain the highest rates of net primary production (NPP) on Earth, contain c. 30% of terrestrial C stocks, and have some of the largest stores of fine-root biomass globally, as well as higher fine-root production and turnover rates compared with other biomes. Tropical forest responses to projected warming, altered rainfall regimes, and elevated CO 2 concentrations are likely to be different from other ecosystems because of their unique characteristics (Box 1), making targeted research and model development important for understanding tropical forest–climate feedbacks. There is now a critical mass of long-term global change field experiments and modeling efforts in tropical forests, yet thus far there has been little synthesis, cross-site comparison, or multi-site standardized experimentation among tropical forests to help us understand how these biomes are changing. An organized INSPIRE session at the 108th Annual Meeting of the Ecological Society of America set out to tackle just this. Speakers covered large-scale tropical forest field experiments and modeling efforts, with an emphasis on changes in ecosystem biogeochemistry under warming, drying, elevated atmospheric CO 2 , and changing nutrient status. In this Meeting report, we provide an overview of the large-scale global change experiments presented and highlight the main objectives and opportunities for tropical forest research that emerged, including cross-site comparisons and integration with ecosystem-scale models (Fig. 1).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ultrafast studies of elusive chemical reactions in the gas phase

The chemical composition of the interstellar medium and planetary atmospheres is constantly in flux as atoms and molecules collide and interact with high-energy particles such as electrons, protons, and photons. These transformative processes ultimately lead to the coalescence of molecules and eventually the birth of stars. Our understanding of these chemical ecosystems relies on models that synthesize data from gas-phase experiments, providing insights into reaction cross sections. This Review examines efforts to delve into the fundamental bond-forming and bond-breaking dynamics that occur during bimolecular and electron-initiated reactions. Furthermore, these experiments involve clever approaches to establish a time reference and the collision geometry necessary for tracking atomic motion with femtosecond time resolution. Findings from these efforts enhance present models and improve predictions for molecule-molecule and electron-molecule collisions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Integrated Simulation of PIP-II at Fermilab

We describe progress towards a community software ecosystem for efficient modeling of the Fermilab PIP-II complex for design validation and virtual test stands.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Research Development and Partnership Pilot (RDPP): Developing plans and partnerships for incorporating tree reproduction to understand Earth system change

The goals for this Research Development and Partnership Pilot (RDPP) proposal were to i) learn about Department of Energy (DOE) research and use the PI’s expertise on patterns and environmental drivers of tree reproduction as a basis to make connections with individuals and research groups at DOE National Laboratories, and ii) to develop plans and form partnerships that will both enhance the PI’s research. The objectives of the proj ect were for the PI to i) participate in the Department of Energy's Office of Science program in Biological and Environmental Research training and outreach activities, ii) conduct directed fact-finding on Earth and Environmental Systems Sciences Division (EESSD) research projects and new partnerships with individuals and groups at National Laboratories, ii) conduct meetings with potential collaborators at National Laboratories to discuss EESSD-relevant research ideas, and iv) to develop a plan for future research efforts. During the period of the award, the objectives of the proposal were met, as the PI attended Environmental System Sciences meetings, American Geophysical Union meetings, SPRUCE experiment meetings, and visited the Oak Ridge National Lab. These activities led to increased understanding of the research being conducted by DOE scientists and the PI had meetings with potential collaborators to develop future research plans that leverage expertise at DOE and the research interests of the PI to increase understanding of the role of tree reproduction in future carbon allocation and in tree regeneration models in boreal ecosystems.

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 OSTI2022

Modeling the Joint Effects of Vegetation Characteristics and Soil Properties on Ecosystem Dynamics in a Panama Tropical Forest

Abstract In tropical forests, both vegetation characteristics and soil properties are important not only for controlling energy, water, and gas exchanges directly but also determining the competition among species, successional dynamics, forest structure and composition. However, the joint effects of the two factors have received limited attention in Earth system model development. Here we use a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM‐FATES, to explore how plant traits and soil properties affect tropical forest growth and composition concurrently. A large ensemble of simulations with perturbed vegetation and soil hydrological parameters is conducted at the Barro Colorado Island, Panama. The simulations are compared against observed carbon, energy, and water fluxes. We find that soil hydrological parameters, particularly the scaling exponent of the soil retention curve ( B sw ), play crucial roles in controlling forest diversity, with higher B sw values (>7) favoring late successional species in competition, and lower B sw values (1 ∼ 7) promoting the coexistence of early and late successional plants. Considering the additional impact of soil properties resolves a systematic bias of FATES in simulating sensible/latent heat partitioning with repercussion on water budget and plant coexistence. A greater fraction of deeper tree roots can help maintain the dry‐season soil moisture and plant gas exchange. As soil properties are as important as vegetation parameters in predicting tropical forest dynamics, more efforts are needed to improve parameterizations of soil functions and belowground processes and their interactions with aboveground vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Valued peaks: Sustainable water allocation for small hydropower plants in an era of explicit ecological needs

Optimizing hydropower operations to balance economic profitability and support functioning ecosystem services is integral to river management policy. In this article, we propose a dynamic, constrained optimization framework for small hydropower plants (SHPs) to evaluate trade-offs between economic profitability and socio-ecological requirements. Specifically, we examine the balance between short-term losses in hydropower generation and the potential for compensatory benefits in the form of revenue from recreational ecosystem services, irrespective of the direct beneficiary. Our framework integrates a fish habitat model, a hydropower optimization model, and a recreational ecosystem service estimate to evaluate different environmental flow scenarios. The optimization process gives three outflow release scenarios, informed by previous streamflow realisations (dam inflow), and designed environmental flow constraints. The framework is applied and tested for the river Kuusinkijoki in North-eastern Finland, which is a habitat for migratory brown trout and grayling populations. We show that the revenue loss due to the environmental flow constraints arises through a reduction in revenue per generated energy unit and through a reduction in turbine efficiency. Additionally, the simulation results reveal that all the designed environmental flow constraints cannot be met simultaneously. Under the environmental flow scenario with both minimum flow and flow ramping rate constraints, the annual hydropower revenue decreases by 16.5 %. An annual increase of 8 % in recreational fishing visits offsets the revenue loss. In conclusion, the developed framework provides knowledge of the costs and benefits of hydropower environmental flow constraints and guides the prioritizing process of environmental measures.

13 HYDRO ENERGY↗
DOE OSTI2025

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Machine learning models inaccurately predict current and future high-latitude C balances

The high-latitude carbon (C) cycle is a key feedback to the global climate system, yet because of system complexity and data limitations, there is currently disagreement over whether the region is a source or sink of C. Recent advances in big data analytics and computing power have popularized the use of machine learning (ML) algorithms to upscale site measurements of ecosystem processes, and in some cases forecast the response of these processes to climate change. Due to data limitations, however, ML model predictions of these processes are almost never validated with independent datasets. To better understand and characterize the limitations of these methods, we develop an approach to independently evaluate ML upscaling and forecasting. We mimic data-driven upscaling and forecasting efforts by applying ML algorithms to different subsets of regional process-model simulation gridcells, and then test ML performance using the remaining gridcells. In this study, we simulate C fluxes and environmental data across Alaska using ecosys, a process-rich terrestrial ecosystem model, and then apply boosted regression tree ML algorithms to training data configurations that mirror and expand upon existing AmeriFLUX eddy-covariance data availability. We first show that a ML model trained using ecosys outputs from currently-available Alaska AmeriFLUX sites incorrectly predicts that Alaska is presently a modeled net C source. Increased spatial coverage of the training dataset improves ML predictions, halving the bias when 240 modeled sites are used instead of 15. However, even this more accurate ML model incorrectly predicts Alaska C fluxes under 21st century climate change because of changes in atmospheric CO 2 , litter inputs, and vegetation composition that have impacts on C fluxes which cannot be inferred from the training data. Our results provide key insights to future C flux upscaling efforts and expose the potential for inaccurate ML upscaling and forecasting of high-latitude C cycle dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A multi-objective optimization model for cropland design considering profit, biodiversity, and ecosystem services

More sustainable agricultural methods are needed to alleviate the decreases in biodiversity and ecosystem services that have occurred because of industrial agriculture. One such method is the inclusion of alternative crops into croplands that can support biodiversity, reduce erosion and chemical runoff, and sequester carbon in the soil. However, the question of where such crops should be planted to balance competing economic and environmental objectives remains open. To this end, we develop a mixed-integer quadratically constrained program to optimize the layout of a cropland considering economic, biodiversity, greenhouse gas emissions, and water quality objectives. We include spatially varying fertilization as a decision variable in addition to crop establishment location. We further include the effect of core area and edges between different crops on biodiversity. To demonstrate the applicability of the model, we apply it to an example field, showing how the optimal cropland design changes as a decision-maker prioritizes different objectives and as edges have different impacts on biodiversity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Diffusion‐Based Uncertainty Quantification Method to Advance E3SM Land Model Calibration

Abstract Calibrating land surface models and accurately quantifying their uncertainty are crucial for improving the reliability of simulations of complex environmental processes. This, in turn, advances our predictive understanding of ecosystems and supports climate‐resilient decision‐making. Traditional calibration methods, however, face challenges of high computational costs and difficulties in accurately quantifying parameter uncertainties. To address these issues, we develop a diffusion‐based uncertainty quantification (DBUQ) method. Unlike conventional generative diffusion methods, which are computationally expensive and memory‐intensive, DBUQ innovates by formulating a parameterized generative model and approximates this model through supervised learning, which enables quick generation of parameter posterior samples to quantify its uncertainty. DBUQ is effective, efficient, and general‐purpose, making it suitable for site‐specific ecosystem model calibration and broadly applicable for parameter uncertainty quantification across various earth system models. In this study, we applied DBUQ to calibrate the Energy Exascale Earth System Model land model at the Missouri Ozark AmeriFlux forest site. Results indicated that DBUQ produced accurate parameter posterior distributions similar to those from Markov Chain Monte Carlo sampling but with 30 times less computing time. This significant improvement in efficiency suggests that DBUQ can enable rapid, site‐level model calibration at a global scale, enhancing our predictive understanding of climate impacts on terrestrial ecosystems.

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 OSTI2020

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

Coupled Long-Term Experiment and Model Investigation of the Differential Response of Plants and Soil Microbes in a Changing Permafrost Tundra Ecosystem

The major research goal of this project was to understand and quantify the fate of carbon stored in permafrost ecosystems using a combination of field and laboratory experiments to measure and model isotope ratios and carbon fluxes in a tundra ecosystem exposed to experimental warming. Field measurements centered on a two-factor experimental warming to increase air and soil temperatures alone, and in combination, at a tundra field site in Alaska. A second manipulation of water table was embedded within the warming treatment such that both major environmental factors (temperature, moisture) controlling ecosystem carbon dynamics were experimentally altered. The experiment was interfaced with modeling activities using both data assimilation and forward modeling approaches. Models were used to make forecasts of ecosystem carbon dynamics beyond the time frame and environmental space of the experiment itself. This project was a follow-on to previous work that made use of the same field manipulation; here we continued the experimental manipulation in order to test new hypotheses regarding long-term effects of warming and permafrost thaw and subsequent changes in moisture availability on ecosystem carbon dynamics. These results were linked with ongoing synthesis and model intercomparison activities through the Permafrost Carbon Network, which led to additional synthesis publications as a result of collaboration and data sharing from this project. As key outcomes, this project has produced new papers published in the peer-reviewed literature, archived datasets used for model-data intercomparisons, trained graduate students and postdoctoral researchers, and has raised awareness about the vulnerability of permafrost carbon to a wider audience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microbial contribution to post-fire tundra ecosystem recovery over the 21st century

Abstract Tundra ecosystems have experienced an increased frequency of fire, and this trend is predicted to continue throughout the 21 st Century. Post-fire recovery is underpinned by complex interactions between microbial functional groups that drive nutrient cycling. Here we use a mechanistic model to demonstrate an acceleration of the nitrogen cycle post-fire driven by changes in niche space and microbial competitive dynamics. We show that over the first 5-years post-fire, fast-growing bacterial heterotrophs colonize regions of the soil previously occupied by slower-growing saprotrophic fungi. The bacterial heterotrophs mineralize organic matter, releasing nutrients into the soil. This pathway outweighs new sources of nitrogen and facilitates the recovery of plant productivity. We broadly show here that while consideration of distinct microbial metabolisms related to carbon and nutrient cycling remains rare in terrestrial ecosystem models, they are important when considering the rate of ecosystem recovery post-disturbance and the feedback to soil nutrient cycles on centennial timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scaling Global Warming Impacts on Ocean Ecosystems: Lessons From a Suite of Earth System Models

An important technique used by climate modelers to isolate the impacts of increasing greenhouse gasses on Earth System processes is to simulate the impact of an abrupt increase in carbon dioxide. The spatial pattern of change provides a “fingerprint” that is generally much larger than natural variability. Insofar as the response to radiative forcing is linear (the impact of quadrupling CO 2 is twice the impact of doubling CO 2 ) this fingerprint can then be used to estimate the impact of historical greenhouse gas forcing. However, the degree to which biogeochemical cycles respond linearly to radiative forcing has rarely been tested. In this paper, we evaluate which ocean biogeochemical fields are likely to respond linearly to changing radiative forcing, which ones do not, and where linearity breaks down. We also demonstrate that the representation of lateral mixing by mesoscale eddies, which varies significantly across climate models, plays an important role in modulating the breakdown of linearity. Globally integrated surface rates of biogeochemical cycling (primary productivity, particulate export) respond in a relatively linear fashion and are only moderately sensitive to mixing. By contrast, the habitability of the interior ocean (as determined by hypoxia and calcite supersaturation) behaves non-linearly and is very sensitive to mixing. This is because the deep ocean, as well as certain regions in the surface ocean, are very sensitive to the magnitude of deep wintertime convection. The cessation of convection under global warming is strongly modulated by the representation of eddy mixing.

54 ENVIRONMENTAL SCIENCES↗