Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wind Speed”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions

Wind profiling by Doppler lidar is common practice and highly useful in a wide range of applications. Airborne Doppler lidar can provide additional insights relative to ground-based systems by allowing for spatially distributed and targeted measurements. Providing a link between theory and measurement, a first large eddy simulation (LES)-based airborne Doppler lidar simulator (ADLS) has been developed. Simulated measurements are conducted based on LES wind fields, considering the coordinate and geometric transformations applicable to real-world measurements. The ADLS provides added value as the input truth used to create the measurements is known exactly, which is nearly impossible in real-world situations. Thus, valuable insight can be gained into measurement system characteristics as well as retrieval strategies. As an example application, airborne Doppler lidar wind profiling is investigated using the ADLS. For commonly used airborne velocity azimuth display (AVAD) techniques, flow homogeneity is assumed throughout the retrieval volume, a condition which is violated in turbulent boundary layer flow. Assuming an ideal measurement system, the ADLS allows to isolate and evaluate the error in wind profiling which occurs due to the violation of the flow homogeneity assumption. Overall, the ADLS demonstrates that wind profiling is possible in turbulent wind field conditions with reasonable errors (root mean squared error of 0.36 m s –1 for wind speed when using a commonly used system setup and retrieval strategy for the conditions investigated). Nevertheless, flow inhomogeneity, e.g., due to boundary layer turbulence, can cause an important contribution to wind profiling error and is non-negligible. Results suggest that airborne Doppler lidar wind profiling at low wind speeds (<5 m s –1 ) can be biased, if conducted in regions of inhomogeneous flow conditions.

17 WIND ENERGY↗
DOE OSTI2020

The Effects of Wind Veer During the Morning and Evening Transitions

Direction and speed wind shear modify turbine performance by changing inflow conditions on turbine blades. Using observations from the 2013 CWEX campaign, we found the daily atmospheric boundary layer transitions (morning and evening) match periods of high electricity demand for a wind farm in central Iowa. Power production during these periods was undermined for large direction shear and low speed shear scenarios. The morning transition displayed larger direction shear over the rotor layer for most wind speeds compared to the evening period. This resulted in lower turbine performance for the morning compared to the evening and whole day. This study shows that the combined effect of direction and speed shear is affecting turbine operation during high electricity demand times of day, predominantly during the morning at this wind farm.

17 WIND ENERGY↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗
DOE OSTI2023

Validation of turbulence intensity as simulated by the Weather Research and Forecasting model off the US northeast coast

We report turbulence intensity (TI) is often used to quantify the strength of turbulence in wind energy applications and serves as the basis of standards in wind turbine design. Thus, accurately characterizing the spatiotemporal variability in TI should lead to improved predictions of power production. Nevertheless, turbulence measurements over the ocean are far less prevalent than over land due to challenges in instrumental deployment, maintenance, and operation. Atmospheric models such as mesoscale (weather prediction) and large-eddy simulation (LES) models are commonly used in the wind energy industry to assess the spatial variability of a given site. However, the TI derivation from atmospheric models has not been well examined. An algorithm is proposed in this study to realize online calculation of TI in the Weather Research and Forecasting (WRF) model. Simulated TI is divided into two components depending on scale, including sub-grid (parameterized based on turbulence kinetic energy (TKE)) and grid resolved. The sensitivity of sea surface temperature (SST) on simulated TI is also tested. An assessment is performed by using observations collected during a field campaign conducted from February to June 2020 near the Woods Hole Oceanographic Institution Martha's Vineyard Coastal Observatory. Results show that while simulated TKE is generally smaller than the lidar-observed value, wind speed bias is usually small. Overall, this leads to a slight underestimation in sub-grid-scale estimated TI. Improved SST representation subsequently reduces model biases in atmospheric stability as well as wind speed and sub-grid TI near the hub height. Large TI events in conjunction with mesoscale weather systems observed during the studied period pose a challenge to accurately estimating TI from models. Due to notable uncertainty in accurately simulating those events, this suggests summing up sub-grid and resolved TI may not be an ideal solution. Efforts in further improving skills in simulating mesoscale flow and cloud systems are necessary as the next steps.

17 WIND ENERGY↗
DOE OSTI2019

An Analysis of the Processes Affecting Rapid Near-Surface Water Vapor Increases during the Afternoon to Evening Transition in Oklahoma

Abstract This study used 20 years of Oklahoma Mesonet data to investigate the changes of near-surface water vapor mixing ratio qυ during the afternoon to evening transition (AET). Similar to past studies, increases in qυ are found to occur near sunset. However, the location, magnitude, and timing of the qυ maximum occurring during the AET are shown to be dependent on the seasonal growth and harvest of vegetation across Oklahoma in the spring and summer months. Particularly, the late spring harvest of winter wheat grown in Oklahoma appears to modify the relative contribution of local and nonlocal processes on qυ. By analyzing time series of qυ during the AET, it is found that the likelihood of a presunset qυ maximum is strongly dependent upon vegetation, soil moisture, wind speed, and cloud cover. Analysis also reveals that the increase in qυ during the AET can increase the parcel conditional instability despite the surface cooling produced by loss of insolation. Next to known changes in low-level wind shear, these changes in instability and moisture demonstrate new ways the AET can modify the presence of the key ingredients relevant to explaining the climatological increase in severe convective storm hazards around sunset.

Blumberg, W. G.↗
DOE OSTI2024

Measuring the Effects of Wake Steering With Nacelle-Mounted Scanning Lidar

In large wind plants, wakes from upwind turbines affect downwind turbines by reducing wind speeds and increasing turbulence. Wake steering may mitigate this effect by deflecting the wakes of upwind turbines away from downwind units. Here, we characterize the impact of wakes from yawed turbines at a commercial-scale wind plant under varying atmospheric and turbine operating conditions. Six months of nacelle-based lidar measurements were collected as part of a field campaign in 2019-2020 in the northern US Great Plains to test the effectiveness of wake steering. We separate these lidar scans by atmospheric stability and turbine operating condition to summarize how yawed wind turbine wakes vary with these input parameters in the atmospheric boundary layer. We summarize the impact of wake steering on various wake characteristics including velocity deficit, wake width, and wake center as retrieved from these lidar data. Yawed wakes have significantly different centerline characteristics compared to unyawed wakes, with large regions of meander in the mid wake region. Yawed wakes are also deflected farther than unyawed wakes in less turbulent conditions with velocity deficits persisting further downstream in stable atmospheric conditions. Overall, yawed wind turbine wakes are larger and wider in lower wind speed environments, and yawed wakes are deflected farther in less turbulent conditions, suggesting that wake steering is most effective in stable atmospheric stratification.

17 WIND ENERGY↗
DOE OSTI2023

A Comparison of Adverse Meteorology Inputs for Emergency Management and Safety Basis Dispersion Modeling Applications at SRS

The MELCOR Accident Consequence Calculation System (MACCS), used for safety basis accident analyses at SRS, implements adjustment factors to consider the increased vertical turbulence associated with local surface roughness and the increased lateral turbulence associated with time-based plume meander. The two models used for SRS emergency planning radiological consequence assessment, i.e., Hotspot and Puff/Plume, currently do not have an option to input or apply diffusion adjustments for meander or surface roughness. Adverse meteorology conditions at SRS have been determined as E stability class and 1.3 meters per second wind speed for ground level releases. This is the input used for emergency planning models. Contrastingly safety basis analyses use a different input, representing adverse meteorology values as F stability class and wind speeds of 1.3 meters per second (for a ground level release). This results in an apparent discrepancy between the two approaches. In this report we document the differing modeling assumptions, discuss the logic behind each modeling approach used to represent local roughness and plume meander conditions, and numerically compare dispersion coefficients and concentration estimates obtained from each set of model inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

CONDENSATION-EVAPORATION MODELING OF DRIFT FROM A MECHANICAL COOLING TOWER

Numerical studyfor investigatingthe phase change behavior of a drift from a mechanical draft cooling tower is presented. A k-ε turbulence model, discrete phase model, species transport equations, and a heat-mass exchange model are solved simultaneouslyto capture the heat and mass transfer effects of a drift condensing and evaporatingeffects. The simulation result is compared to an available real scale test data and its limitation is also discussed. A parametric study is carried out by varying the ambient dry bulb temperature, ambient relative humidity, wind speed, liquid droplet diameter, and the drift amount. It is found that by combining the evaporation and condensation models, a drift phase change phenomena is accurately depicted. A parametric study reveals that wind speed impact is affected by the ambient temperature and the number of liquid droplets and their size play a critical role in the phase change phen

Han, Kai↗
DOE OSTI2020

A digital twin based on OpenFAST linearizations for real-time load and fatigue estimation of land-based turbines

Monitoring a wind turbine requires intensive instrumentation, which would be too cost-prohibitive to deploy to an entire wind plant. This work presents a technique that uses readily available measurements to estimate signals that would otherwise require additional instrumentation. This study presents a digital twin concept with a focus on estimating wind speed, thrust, torque, tower-top position, and loads in the tower using supervisory control and data acquisition (SCADA) measurements. The model combines a linear state-space model obtained using OpenFAST linearizations, a wind speed estimator, and a Kalman filter algorithm that integrates measurements with the state model to perform state estimations. The measurements are: top acceleration, generator torque, pitch, and rotational speed. The article extends previous work that derived the linear state-space model using a different method. The new implementation, based on OpenFAST linearization capability, allows for a systematic extension of the method to more states, inputs, outputs, and to the offshore environment. Results from the two methods are compared, and the validation is made with additional measurements using the GE 1.5-MW turbine located at the National Renewable Energy Laboratory test site. Real-time damage equivalent loads of the tower bottom moment are estimated with an average accuracy of approximately 10%. Overall, the results from this proof of concept are encouraging, and further application of the model will be considered.

17 WIND ENERGY↗
DOE OSTI2024

Validation of Weather Forecasting Products

Meteorological modeling plays a pivotal role in operational safety and emergency response at the Savannah River Site (SRS). This study focuses on verifying the Regional Atmospheric Modelling System (RAMS) Version 4.3 through Mean Bias Error (MBE) and Root Mean Square Error (RMSE) analyses of temperature, dew point, and wind speed over a decade. Using observed data from SRS, we assessed RAMS' accuracy, revealing seasonal biases and error trends. Results indicate RAMS' strengths in mild weather conditions but challenges during seasonal extremes and wind speed predictions due to measurement disparities. Future research aims to expand verification to other models and parameters, advocating for enhanced forecasting accuracy crucial for safeguarding personnel and community well-being.

42 ENGINEERING↗
DOE OSTI2025

Evaluating the Trustworthiness of Explainable Artificial Intelligence (XAI) Methods Applied to Regression Predictions of Arctic Sea Ice Motion

Abstract Recent advances in explainable artificial intelligence (XAI) methods show promise for understanding predictions made by machine learning (ML) models. XAI explains how the input features are relevant or important for the model predictions. We train linear regression (LR) and convolutional neural network (CNN) models to make 1-day predictions of sea ice velocity in the Arctic from inputs of present-day wind velocity and previous-day ice velocity and concentration. We apply XAI methods to the CNN and compare explanations to variance explained by LR. We confirm the feasibility of using a novel XAI method [i.e., global layerwise relevance propagation (LRP)] to understand ML model predictions of sea ice motion by comparing it to established techniques. We investigate a suite of linear, perturbation-based, and propagation-based XAI methods in both local and global forms. Outputs from different explainability methods are generally consistent in showing that wind speed is the input feature with the highest contribution to ML predictions of ice motion, and we discuss inconsistencies in the spatial variability of the explanations. Additionally, we show that the CNN relies on both linear and nonlinear relationships between the inputs and uses nonlocal information to make predictions. LRP shows that wind speed over land is highly relevant for predicting ice motion offshore. This provides a framework to show how knowledge of environmental variables (i.e., wind) on land could be useful for predicting other properties (i.e., sea ice velocity) elsewhere. Significance Statement Explainable artificial intelligence (XAI) is useful for understanding predictions made by machine learning models. Our research establishes trustability in a novel implementation of an explainable AI method known as layerwise relevance propagation for Earth science applications. To do this, we provide a comparative evaluation of a suite of explainable AI methods applied to machine learning models that make 1-day predictions of Arctic sea ice velocity. We use explainable AI outputs to understand how the input features are used by the machine learning to predict ice motion. Additionally, we show that a convolutional neural network uses nonlinear and nonlocal information in making its predictions. We take advantage of the nonlocality to investigate the extent to which knowledge of wind on land is useful for predicting sea ice velocity elsewhere.

Hoffman, Lauren [Scripps Institution of Oceanograp↗
DOE OSTI2019

Improving Wind Energy Forecasting through Numerical Weather Prediction Model Development

The primary goal of the Second Wind Forecast Improvement Project (WFIP2) is to advance the state-of-the-art of wind energy forecasting in complex terrain. To achieve this goal, a comprehensive 18-month field measurement campaign was conducted in the region of the Columbia River basin. The observations were used to diagnose and quantify systematic forecast errors in the operational High-Resolution Rapid Refresh (HRRR) model during weather events of particular concern to wind energy forecasting. Examples of such events are cold pools, gap flows, thermal troughs/marine pushes, mountain waves, and topographic wakes. WFIP2 model development has focused on the boundary layer and surface-layer schemes, cloud–radiation interaction, the representation of drag associated with subgrid-scale topography, and the representation of wind farms in the HRRR. Additionally, refinements to numerical methods have helped to improve some of the common forecast error modes, especially the high wind speed biases associated with early erosion of mountain–valley cold pools. This study describes the model development and testing undertaken during WFIP2 and demonstrates forecast improvements. Specifically, WFIP2 found that mean absolute errors in rotor-layer wind speed forecasts could be reduced by 5%–20% in winter by improving the turbulent mixing lengths, horizontal diffusion, and gravity wave drag. The model improvements made in WFIP2 are also shown to be applicable to regions outside of complex terrain. Ongoing and future challenges in model development will also be discussed.

17 WIND ENERGY↗
DOE OSTI2021

Implementation of an Orographic Drag Scheme Considering Orographic Anisotropy in All Flow Directions in the Earth System Model CAS‐ESM 2.0

Abstract A reasonable representation of orographic anisotropy in earth system models is vital for improving weather and climate modeling. In this study, we implement the orographic drag scheme, including 3‐D orographic anisotropy (3D‐AFD), into the Chinese Academy of Sciences Earth System Model version 2 (CAS‐ESM 2.0). Three groups of simulations named sensitivity run, medium‐range forecast, and seasonal forecast respectively are conducted using the updated CAS‐ESM model and validated against station observation and reanalysis data. These simulations are run using the three schemes—3D‐AFD scheme, the 3D orographic anisotropy scheme for the eight‐direction (3D‐8x), and the 2D isotropic scheme (2D)—to compare their performance in CAS‐ESM 2.0. Sensitivity runs show that the 3D‐AFD provide more continuous transition of calculated drag as function of wind direction than the 3D‐8x, while the drag calculated using the 2D scheme show no variation with change of wind direction. Enhanced drag in the medium range and seasonal forecast using the updated CAS‐ESM alleviates part of the winter wind speed bias over the Tibetan Plateau (TP) and the cold bias over TP and the Siberian polar region. It is shown that the 3D‐AFD scheme alleviates more bias than that of the 3D‐8x scheme (by wind speed reduction of 1 ∼ 2 m/s and temperature of 1 ∼ 2 K) especially in the seasonal forecast. The results suggest that reasonable representation of the orographic anisotropy is important in climate modeling.

Xie, Jinbo↗
DOE OSTI2023

Simulations for Planning of Liquid Hydrogen Spill Test

In order to better understand the complex pooling and vaporization of a liquid hydrogen spill, Sandia National Laboratories is conducting a highly instrumented, controlled experiment inside their Shock Tube Facility. Simulations were run before the experiment to help with the planning of experimental conditions, including sensor placement and cross wind velocity. This paper describes the modeling used in this planning process and its main conclusions. Sierra Suite’s Fuego, an in-house computational fluid dynamics code, was used to simulate a RANS model of a liquid hydrogen spill with five crosswind velocities: 0.45, 0.89, 1.34, 1.79, and 2.24 m/s. Two pool sizes were considered: a diameter of 0.85 m and a diameter of 1.7. A grid resolution study was completed on the smaller pool size with a 1.34 m/s crosswind. A comparison of the length and height of the plume of flammable hydrogen vaporizing from the pool shows that the plume becomes longer and remains closer to the ground with increasing wind speed. The plume reaches the top of the facility only in the 0.45 m/s case. From these results, we concluded that it will be best for the spacing and location of the concentration sensors to be reconfigured for each wind speed during the experiment.

08 HYDROGEN↗
DOE OSTI2019

National Wind Technology Center-Characterization of Atmospheric Conditions

Data provided by the instrumentation located at the National Wind Technology Center (NWTC) is compiled below for the purposes of characterizing the atmospheric operating conditions at the National Renewable Energy Laboratories (NREL) wind site. Conditions are presented conditioned on month of the year to show seasonal variation, and based on characteristic quantities seen in the atmospheric flow including, bulk wind speed, direction, and atmospheric stability. All data reported and analyzed in the following document were collected from instrumentation mounted on the 135 m-tall masts erected at the NWTC near Boulder, Colorado. Each of the two met-masts (named M4 and M5) were sited specifically to measure the inflow wind to several of the utility-scale turbines at the NWTC. The construction, operation, and maintenance of the masts has been supported since their commissioning in 2011 by the Department of Energy, Office of Energy Efficiency and Renewable Energy, Wind and Water Power Program. Data from the masts are publicly available for download or basic visualization, hosted on NREL's website (https://nwtc.nrel.gov/135mData). Data from the masts have been used in the past for a wide range of projects including wind turbine power performance measurements, high fidelity model development and validation, and atmospheric science research. An interactive visualization platform has been developed on NREL's Open Energy Information platform and can be used to reproduce many of the graphics included in the current work. This report describes in detail the atmospheric conditions recorded by M5, the inflow mast for the DOE-owned GE 1.5 MW wind turbine. The report focuses on data recorded from August 2012 to May 2018 and includes: an introduction to the wind climate at the NWTC; a physical description of the met masts; a review of installed instrumentation; and data collection, processing, and quality control methods. In addition to a concise description of the instrumentation from the met masts and the data channels they provide, this report characterizes the behavior of the atmospheric resource first in a cumulative sense, discussing mainly the distributions wind speeds, flow directions and thermodynamic properties at several heights along the towers. Periodicity in atmospheric conditions is then detailed through analysis of diurnal and annual cycles of the data streams. Periods of both low turbulence and extreme atmospheric events that are of interest in terms of wind turbine operation and wind energy research are detailed in the final chapter of this report. Low turbulence events are somewhat rare at the NWTC and are of interest from the perspective of validating the operation of new wind turbine and blade designs. Similarly, the NTWC is known for a relatively higher rate of occurrence of extreme weather events.

17 WIND ENERGY↗
DOE OSTI2023

From Hydrometeor Size Distribution Measurements to Projections of Wind Turbine Blade Leading-Edge Erosion

Wind turbine blade leading-edge erosion (LEE) is a cause of increased operation and maintenance costs and decreased annual energy production. Thus, detailed, site-specific quantification of likely erosion conditions are critically needed to inform wind plant owner/operator decisions regarding mitigation strategies. Estimating the damage potential at a wind plant site requires accurate measurement of precipitation intensity, phase, droplet size distributions, wind speeds and their joint statistics. The current work quantifies the effect of disdrometer type on the characterization of LEE potential at a site in the US Southern Great Plains. using observations from three co-located disdrometers (an optical, an impact and a video disdrometer), along with hub-height wind-speed observations from a Doppler lidar and two LEE models: a kinetic energy model and the Springer model. Estimates of total kinetic energy of hydrometeor impacts over the four-year study period vary by as much as 38%, and coating lifetime derived from accumulated distance-to-failure estimates from the Springer model differ by an even greater amount, depending on disdrometer type. Damage potential at this site is concentrated in time, with 50% of impact kinetic energy occurring in 6–12 h per year, depending on which set of disdrometer observations is used. Rotor-speed curtailment during the most erosive 0.1–0.2% of 10 min periods is found to increase blade lifetimes and lead to the lowest levelized cost of energy.

17 WIND ENERGY↗
DOE OSTI2025

Sea spray aerosol production flux retrieval based on Doppler lidar measurements

Supermicron sea spray aerosols (SSA) play a crucial role in stratocumulus cloud microphysics and drizzle formation. However, our understanding of SSA's spatiotemporal distribution, production mechanisms, and fluxes in the marine boundary layer remains limited. This study introduces a novel approach to determine supermicron SSA production flux and number concentration at various heights above the ocean surface. Our method leverages Doppler lidar data, specifically attenuated backscatter and vertical velocity, collected at a 105 m range gate during the Eastern Pacific Clouds and Aerosols for Precipitation Experiment (EPCAPE) at Scripps Pier, La Jolla, California. To focus on nascent SSA, we analyzed data from periods with wind directions between 225° and 315° and wind speeds exceeding 4 m s −1 . Cloud and precipitation events were excluded using relative humidity, backscatter thresholds, and 2D Video Disdrometer observations. Results show that calculated supermicron SSA number concentrations at 105 m ranged from 0.26 to 2.93 cm −3 for surface wind speeds between 4 and 6 ms −1 . When exceeding the method's Limit of Detection, estimated production fluxes ranged from 0.20 to 1.53 cm −2 s −1 . This methodology has potential for application in various oceanic locations, likely enhancing our understanding of SSA's role in atmospheric processes and climate interactions.

54 ENVIRONMENTAL SCIENCES↗