Search NASA⌕ Search

SEARCH · Search NASA

Results for “peak period”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sumner County Field Validation: An Assessment of Thermal Energy Storage in Municipal Buildings in Mixed-Humid Climates

The overall aim of the project was to conduct a full-scale implementation of TES in an operational administrative building in Sumner County, Kansas with a special focus of phase change material (PCM) passive implementation in the ceiling envelope. PCM was used with the intention of increasing the ability of public buildings to shift the peak HVAC energy demand and consumption. The study was conducted during the span of three years in 3 phases, (i) phase I: energy audit (pre-retrofit) (ii) phase II: ceiling tile installation (PCM retrofit), and (iii) phase III: energy evaluation (post-retrofit). Temperature and HVAC energy measurements as well as weather data collected is used for an extensive validation study and further parametric assessment to find (a) annual HVAC energy savings (b) peak period electricity and load savings focusing on cooling dominant months of the year.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Integrating AI Data Centers with the Power Grid

The rapid expansion of artificial intelligence (AI) has triggered an unprecedented surge in electricity demand, with US data center energy use projected to double or triple 2023 levels by 2028. This exponential growth places strain on grid infrastructure, which can hinder timely construction of desired computing capacity. To bridge this supply-demand gap, utilities and AI developers are increasingly turning to demand flexibility, a strategy that incentivizes shifting or reducing power use during peak periods of grid stress. Data centers are uniquely equipped for flexible operations due to their digital workloads, built-in redundancy, and onsite energy assets. This article outlines four primary mechanisms to enable data center flexibility: computational load flexibility (shifting tasks temporally or geographically), flexible use of core facility infrastructure adjustments, energy storage utilization, and onsite electricity generation. To encourage adoption, utilities are deploying new tariff designs, including voluntary interruptible service riders, mandated flexibility requirements, and streamlined interconnection processes for flexible loads. For the highly capitalized and rapidly growing AI industry, the primary motivators for embracing these strategies are expediting facility interconnection, satisfying emerging regulatory mandates, and mitigating community resistance. While demand flexibility cannot substitute the long-term need for new bulk power generation, it serves as an essential, immediate solution for enabling near-term deployment. By transforming data centers from grid stressors into stabilizing assets, flexible operations can ensure reliable grid integration, ease market pressures, and support a resilient power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Salt Hydrate Eutectic Thermal Energy Storage for Building Thermal Regulation (Final Technical Report)

Thermal energy storage is anticipated to play an important role in developing the power grid of the future - a power grid that meets increasing demands of users, is resistant to disruptions, but also allows for greater penetration of renewable resources. Specifically, thermal energy storage materials can be integrated into HVAC systems and building envelopes, where they can be used to shift power demands for building climate control from periods of peak demand to periods of low demand. Phase change materials (PCMs) are compelling as low-cost, high energy density thermal energy storage materials for building thermal management. However, there is a lack of high performance low-cost PCMs within the specific temperature ranges which would most effectively allow for power load shifting. Inorganic salt hydrates represent a promising class of PCMs, but their inherent limitations cause them to be currently unavailable for reliable building applications. The overarching goals of this research effort are to: 1) Discover low-cost, high volumetric density salt hydrate eutectic PCMs to store low-quality heat (10 to 40 °C); 2) Introduce a high thermal conductivity matrix to reduce the time constant for energy storage to ~0.1 to 1 hr, incorporate nucleation catalysts to decrease undercooling, and utilize microencapsulation and shape stabilization approaches, minimizing moisture loss/gain, mitigating phase separation, and maintaining stable melting behavior over the lifetime of the compounds; 3) Evaluate the impact these systems have on peak load shifting, and the potential for overall energy savings under different climatic scenarios and building configurations. These goals will be achieved by an integrated research program consisting of six cohesive research subtasks: 1) Materials discovery of eutectic salt hydrate PCMs by using computationally predicted thermodynamic equilibria, coupled with high-throughput experimental validation, 2) Rapid experimental screening of nucleation catalysts identified through robust computational databases, 3) Embedding salt hydrate PCM into a low cost and scalable high conductivity matrix, 4) Microencapsulation of salt hydrate microspheres using hybrid inorganic-polymer microencapsulation approach, 5) Shape stabilization by thermoreversible salt hydrate salogels, and 6) Analysis of end-use using thermal simulations, and characterization of mock-up energy storage finished components.

25 ENERGY STORAGE↗
DOE OSTI2020

Combined Heat and Power Technology Fact Sheet Series: Thermal Energy Storage

This fact sheet provides an overview of thermal energy storage (TES) technologies, which heat or cool a storage medium and, when needed, deliver the stored thermal energy to meet heating or cooling needs. TES systems are used in commercial buildings, industrial processes, and district energy installations to deliver stored thermal energy during peak demand periods, thereby reducing peak energy use. TES systems are often integrated with electric or absorption chillers to reduce peak electricity costs and, in the case of new construction, to reduce capital costs by optimizing chiller size. TES technologies can support sites that have either renewable or fossil power generation, including combined heat and power (CHP) installations. With CHP, TES can help optimize equipment size by reducing the required peak CHP thermal capacity and increasing annual CHP usage. TES can also provide turbine inlet cooling for gas turbines used in CHP applications, which increases power production in hot ambient conditions.

Combined Heat and Power, CHP, Thermal Energy, Tech↗
DOE OSTI2022

Deficient precipitation sensitivity to Sahel land surface forcings among CMIP5 models

Abstract The overall performance of the simulated seasonal precipitation response to local terrestrial forcings, namely vegetation abundance and soil moisture, in the Sahel among the Coupled Model Intercomparison Project Phase Five (CMIP5) Earth System Models (ESMs) is systematically investigated and compared with its observational counterpart using a multivariate statistical method. The observed seasonal precipitation response is evaluated against a large ensemble of observational, reanalysis, and satellite data sets to provide quantification of uncertainties. The behaviour of models with and without a Dynamic Global Vegetation Model (DGVM) component is also explored, along with the mechanisms responsible for terrestrial feedback on rainfall. In general, the CMIP5 models can reasonably capture the seasonal evolution of Sahel precipitation and soil moisture, albeit with wet biases during the pre‐monsoon period and dry biases during the peak monsoon period. The non‐DGVM ESMs simulate comparable leaf area indices (LAIs) with observations, while DGVM‐enabled ESMs simulate too much year‐round LAI. The variance of precipitation that is attributed to oceanic forcings in CMIP5 is comparable with observations; however, the variance of precipitation that is attributed to terrestrial forcings is smaller in CMIP5 models than observed, especially for non‐DGVM ESMs. CMIP5 models, especially those without DGVMs, undervalue precipitation's observed response strength to soil moisture anomalies. In both observations and CMIP5 models, none of the atmospheric variables show significant responses to direct vegetation forcing, except for the response in transpiration. Although vegetation has minimal direct effect on the atmospheric state, it can affect the atmosphere by modifying soil moisture and transpiration rate indirectly, which helps explain the more realistic simulation of rainfall in DGVM‐enabled ESMs than non‐DGVM ESMs. Coupling of an ESM to a DGVM is critical in generating reasonable land–atmosphere feedback and examining future ecological and climatic changes over the Sahel.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Deep Reinforcement Learning Based HVAC Control for Reducing Carbon Footprint of Buildings

In this paper, we present our work on deep reinforcement learning (DRL) based intelligent control of Heating, Ventilation, and Air Conditioning (HVAC) with the goal of reducing carbon emission. We performed this task using 1) Marginal Operating Emission Rates (MOER), where the objective was to shift the demand to the low emission period of the day and 2) Time-Of-Use (TOU) demand-response price where the objective was to shift the demand to low price period of the day. This was achieved by learning an optimal pre-cooing strategy. We found the carbon emission reduction in the range of 6%-16% depending on the opportunity presented by the MOER signal. Similarly, we observed the carbon emission reduction in the range of 23%-29% during the peak price period when TOU price was used. The results clearly demonstrated the applicability of our approach in reducing the carbon footprint of the building.

carbon emission↗
DOE OSTI2021

Open-Loop Control of Adjustable Tuned Mass Dampers for Floating Wind Turbine Platforms

In this study, we demonstrate the ability of adjustable tuned mass dampers (TMDs) to reduce the platform motion of floating offshore wind turbines (FOWTs). The TMDs are located in the hull and provide control authority by varying the amount of water ballast and compressed air in a reservoir over time. The optimal TMD settings depend on the wind and wave conditions of the FOWT. We present an open-loop control scheme that changes the TMD natural frequency based on sea state and a simple method for estimating the peak wave period and significant wave height. The performance of this open-loop control is compared to ideal control of the TMDs (set knowing the exact sea state), a TMD with a constant natural frequency targeting the design load case with the greatest platform motion, and a baseline set of simulations of the platform without TMDs. The constant natural frequency case performs nearly as effectively as the actively controlled case, suggesting that the TMDs can be simply designed for extreme load conditions and provide consistent control over a range of environments; however, the methods used to find the optimal TMD parameters and estimate sea state statistics could be used in other aspects of FOWT control and design.

artificial intelligence↗
DOE OSTI2020

Conceptual Design Study of a Thermal Energy Storage System Using Granular Solids for Coal-fired Boilers

<span style="font-size: 12pt;">Presentation for the 2020 International Pittsburgh Coal Conference held virtually from September 8-11, 2020.</span><span style="font-size: 12pt;">Abstract: As electricity from various renewable energy sources is increasingly supplied to the power grid, the grid’s stability will be challenged. Proposed grid-stability solutions will need to address variability of sources, such as solar and wind energy, in order to alleviate problems such as overgeneration and ramping. In regions where increases in solar energy are prominent, such as California, the net power demand curve often referred to as “the Duck Curve,” has caused problems for energy dispatchers, due to over-generation and fast ramping demands. In this project, we investigated a conceptual design for a Thermal Energy Storage System (TESS) integrated into a coal-fired power plant to support grid-scale storage. Implementing a TESS into power plants can allow for an increase in renewable sources, increase grid resilience, and allow conventional coal-fired plants to operate more efficiently. For this project, a TESS using granular solids as the energy-storage media is integrated into a conventional coal-fired power plant. The purpose of the TESS is to store energy from the flue gas during low demand periods and transfer the energy to steam during peak demand periods. The storage allows the stabilization of combustion and related boiler components at low loads and provides the opportunity to avoid some shutdowns. This study uses thermodynamic properties of several candidate thermal storage materials to estimate the size of storage system required for several operating scenarios. </span>

Huckaby, E. David↗
DOE OSTI2026

The past, present, and future of peaking thermal power plants in the United States

Today's power systems rely on "peaker plants" to reliably serve load during peak demand periods. In this study we consider the present competitiveness of different peaking options, how much and what kinds of plants have provided U.S. peaking capacity, and potential future peaking fleet compositions. We explore how capital intensity impacts the breakeven capacity factor between two potential resources: combustion turbines (CTs) and combined cycle plants (CCs). CTs outcompete CCs below 12%-17% annual capacity factor at today's prices, but can shift with changes to fuel or start costs. Historically, gas CTs and petroleum steam plants most closely fit the role of peakers. Peaking capacity could grow from approximately 280 GW today to 460-770 GW in 2050 composed of a wider range of resources. We conclude by discussing implications of this shift, with a focus on the potential planning considerations for shifting to a greater utilization of CCs for peaking needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗
DOE OSTI2021

The Evolving Role of Extreme Weather Events in the U.S. Power System with High Levels of Variable Renewable Energy

As weather-dependent renewable generation grows, it is important for power system planning to understand the broad trends and correlations between weather, renewable resources, and load. The traditional planning, performed by utilities and system operators, includes the study of system resource adequacy during peak load periods in the summer and winter to ensure the generation and transmission system is appropriate to meet load. But in a power grid with a high penetration of variable renewable energy (i.e., wind and solar), periods of high risk to system resource adequacy may no longer correspond only to hours of peak load. In particular, high shares of variable renewable energy, even when well-forecasted to inform system operations, can further complicate the stress extreme weather events already place on the grid. They also may lead to changes to the types of weather conditions that are most problematic to system operations and resource adequacy due to widespread and extended deficits of wind and solar generation. Accordingly, the focus of reliability assessments in long-term planning studies may need to evolve in the coming years to more fully incorporate weather events that lead to these deficits. This report seeks to identify these new weather events and understand the characteristics of the events that lead to system risk of future systems with higher penetrations variable renewable energy.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

How the U.S. Power Grid Kept the Lights on in Summer 2024

Maintaining the reliability of the bulk power system, which supplies and transmits electricity, is a critical priority of electric grid planners, operators, and regulators. As demand for electricity increases and the U.S. resource mix changes, how grid operators meet peak demand is changing. In summer 2024, grid operators in all regions maintained enough capacity to keep the lights and air conditioners on during periods of peak demand, even as older generators have been retired. And, an increasing number of regions used more solar and storage to meet peak demand. In this publication, we describe grid operations on the highest demand day in ERCOT and a few other regions and how solar and storage in particular worked together to help meet peak demand.

14 SOLAR ENERGY↗
DOE OSTI2024

Advancing process-based flood frequency analysis for assessing flood hazard and population flood exposure

Recent studies have showcased the use of process-based hydrological models with Stochastic Storm Transposition (SST) techniques to conduct Flood Frequency Analysis (FFA). This framework, referred hereby FFA-SST, has proved to be a robust strategy to estimate peak flows of specific annual exceedance probability (e.g., 100-year peak flow) that can reflect natural and anthropogenic disturbances, including changes in land use and meteorological patterns. With the objective of advancing the FFA-SST framework, this study presents for the first time the use of an Integrated Surface-Subsurface Hydrological Model (ISSHM) to conduct FFA-SST by extending the analysis from peak flow responses to flood extent, enabling a unique view and analysis of flood hazard and population flood exposure at the basin scale. As a proof-of-concept, we used the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct FFA-SST by simulating the flood response to 5,000 annual synthetic storm events in a 2,227 $km^2$ Southeast Texas watershed. We demonstrate that ATS, without site-specific calibration, provides a robust process-based representation of peak flows, flood extent, streamflow, evapotranspiration, soil moisture content, and water storage changes. Our results and analyses, covering frequency curves up to a 500-year return period for peak flows, basin inundation fractions, and the number of people exposed to flooding, offer a unique perspective to analyze flood impacts across spatial scales. Overall, this study provides critical insights for flood risk management by extending the FFA-SST framework to include both flood hazard and population flood exposure analyses at the basin scale. Such an approach will empower stakeholders and disaster emergency agencies with a more comprehensive understanding of flood impacts across the entire basin domain, facilitating informed decision-making for flood risk assessment and management.

58 GEOSCIENCES↗
DOE OSTI2026

Smart building HVAC control challenge: experience and solutions from the ADRENALIN project

A smart building HVAC control competition crowdsourced and compared algorithms on fair and equal ground using the standardized BOPTEST framework. The competition attracted 138 participants, but only 9% submitted valid solutions for the final stage, highlighting the complexity of advanced HVAC control design. The winning solutions showed significant potential to reduce energy use and cost by shifting demand, without compromising occupant comfort. Across scenarios, thermal energy cost reductions of 36–76% relative to a baseline, were achieved. In peak heat periods, the cost reduction leveraged limited energy use reduction (0–15%), but more significant energy price reduction (34–62%). This shows smart controls' ability to avoid as much as possible consumption during the morning peak hours, when spot prices are tendentially the highest. Hosting the competition has highlighted challenges in creating competitions that both are fair and promotes solutions that are transferable to real life implementation.

BOPTEST↗
DOE OSTI2026

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2023

Leveraging Groundwater Dynamics to Improve Predictions of Summer Low-Flow Discharges

Summer streamflow predictions are critical for managing water resources; however, warming-induced shifts from snow to rain regimes impact low-flow predictive models. Additionally, reductions in snowpack drive earlier peak flows and lower summer flows across the western United States increasing reliance on groundwater for maintaining summer streamflow. However, it remains poorly understood how groundwater contributions vary interannually. Here, we quantify recession limb groundwater (RLGW), defined as the proportional groundwater contribution to the stream during the period between peak stream flow and low flow, to predict summer low flows across three diverse western US watersheds. We ask (a) how do snow and rain dynamics influence interannual variations of RLGW contributions and summer low flows?; (b) which watershed attributes impact the effectiveness of RLGW as a predictor of summer low flows? Linear models reveal that RLGW is a strong predictor of low flows across all sites and drastically improves low-flow prediction compared to snow metrics at a rain-dominated site. Results suggest that strength of RLGW control on summer low flows may be mediated by subsurface storage. Subsurface storage can be divided into dynamic (i.e., variability saturated) and deep (i.e., permanently saturated) components, and we hypothesize that interannual variability in dynamic storage contribution to streamflow drives RLGW variability. In systems with a higher proportion of dynamic storage, RLGW is a better predictor of summer low flow because the stream is more responsive to dynamic storage contributions compared to deep-storage-dominated systems. Overall, including RLGW improved low-flow prediction across diverse watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Annual carbon dioxide flux over seasonal sea ice in the Canadian Arctic

Continuous measurements of carbon dioxide (CO 2 ) flux were collected from a 10 m eddy covariance tower in a coastal-marine environment in the Canadian Arctic Archipelago over the course of a 17-month period. The extended length of data collection resulted in a unique dataset that includes measurements from two spring melt and summer seasons and one autumn freeze-up. These field observations were used to verify findings from previous theoretical and laboratory experiments investigating air-sea gas exchange in connection with sea ice. The results corroborated previous findings showing that thick ice cover under winter conditions acts as a barrier to gas exchange. In the spring, CO 2 fluxes were downward (uptake) in both the presence of melt ponds and during ice break-up. However, diurnal cycles were present throughout the early spring melt period, corresponding to the opposing influences of freezing and melting at the ice surface. Fluxes measured during melt periods confirmed previous laboratory tank measurements that showed a gas transfer coefficient of melting ice of 0.4 mol m −2 d −1 atm −1 . Open water CO 2 fluxes showed outgassing in early summer and uptake in mid-to-late summer, tied closely to trends in surface water temperature and its effect on the partial pressure of CO 2 in the water. The autumn period of the field campaign represents the first eddy covariance CO 2 fluxes measured over naturally forming sea ice. Our measurements showed mean upward fluxes (outgassing) of 1.1±1.5 mmol m −2 d −1 associated with the freezing of ice – the same order of magnitude found by previous laboratory tank experiments. However, peak flux periods during ice formation had measured fluxes that were a factor of 3 higher than the tank experiments, suggesting the importance of natural conditions (e.g., wind) on air-ice gas exchange. Conducting an Arctic-wide extrapolation we estimate CO 2 outgassing from the freezing period to be a counterbalance equivalent to 5 to 15 % of the magnitude of the estimated Arctic CO 2 sink. Overall, there was no evidence of dramatically enhanced gas exchange in marginal ice conditions as proposed by previous studies. Although the different seasons showed active CO 2 exchange, there was a balance between upward and downward fluxes at this specific location, resulting in a small net CO 2 uptake over the annual cycle of −0.3 g C m −2 .

54 ENVIRONMENTAL SCIENCES↗