Search NASA⌕ Search

SEARCH · Search NASA

Results for “velocity deficit”

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 37 records · Page 2

Wake Statistics of Different-Scale Wind Turbines under Turbulent Boundary Layer Inflow

Subscale wind turbines can be installed in the field for the development of wind technologies, for which the blade aerodynamics can be designed in a way similar to that of a full-scale wind turbine. However, it is not clear whether the wake of a subscale turbine, which is located closer to the ground and faces different incoming turbulence, is also similar to that of a full-scale wind turbine. In this work we investigate the wakes from a full-scale wind turbine of rotor diameter 80 m and a subscale wind turbine of rotor diameter of 27 m using large-eddy simulation with the turbine blades and nacelle modeled using actuator surface models. The blade aerodynamics of the two turbines are the same. In the simulations, the two turbines also face the same turbulent boundary inflows. The computed results show differences between the two turbines for both velocity deficits and turbine-added turbulence kinetic energy. Such differences are further analyzed by examining the mean kinetic energy equation.

17 WIND ENERGY↗

Power Production, Inter- and Intra-Array Wake Losses from the U.S. East Coast Offshore Wind Energy Lease Areas

There is an urgent need to develop accurate predictions of power production, wake losses and array–array interactions from multi-GW offshore wind farms in order to enable developments that maximize power benefits, minimize levelized cost of energy and reduce investment uncertainty. New, climatologically representative simulations with the Weather Research and Forecasting (WRF) model are presented and analyzed to address these research needs with a specific focus on offshore wind energy lease areas along the U.S. east coast. These, uniquely detailed, simulations are designed to quantify important sources of wake-loss projection uncertainty. They sample across different wind turbine deployment scenarios and thus span the range of plausible installed capacity densities (ICDs) and also include two wind farm parameterizations (WFPs; Fitch and explicit wake parameterization (EWP)) and consider the precise WRF model release used. System-wide mean capacity factors for ICDs of 3.5 to 6.0 MWkm−2 range from 39 to 45% based on output from Fitch and 50 to 55% from EWP. Wake losses are 27–37% (Fitch) and 11–19% (EWP). The discrepancy in CF and wake losses from the two WFPs derives from two linked effects. First, EWP generates a weaker ‘deep array effect’ within the largest wind farm cluster (area of 3675 km2), though both parameterizations indicate substantial within-array wake losses. If 15 MW wind turbines are deployed at an ICD of 6 MWkm−2 the most heavily waked wind turbines generate an average of only 32–35% of the power of those that experience the freestream (undisturbed) flow. Nevertheless, there is no evidence for saturation of the resource. The wind power density (electrical power generation per unit of surface area) increases with ICD and lies between 2 and 3 Wm−2. Second, EWP also systematically generates smaller whole wind farm wakes. Sampling across all offshore wind energy lease areas and the range of ICD considered, the whole wind farm wake extent for a velocity deficit of 5% is 1.18 to 1.38 times larger in simulations with Fitch. Over three-quarters of the variability in normalized wake extents is attributable to variations in freestream wind speeds, turbulent kinetic energy and boundary layer depth. These dependencies on meteorological parameters allow for the development of computationally efficient emulators of wake extents from Fitch and EWP.

17 WIND ENERGY↗

Comparison of Coupled and Uncoupled Modeling of Floating Wind Farms with Shared Anchors

As design options for floating wind farms continue to be explored, shared (or multiline) anchors that secure mooring lines from multiple turbines remain a promising technology that can potentially reduce the number of anchors and overall mooring costs. This study evaluates two methods for analyzing the loads on shared anchors: one in which floating offshore wind turbines are simulated individually (using the software OpenFAST), and one in which an entire floating wind farm is simulated collectively (using the software FAST.Farm). A three-line shared anchor is evaluated for multiple loading scenarios in deep water, using the International Energy Agency 15 MW turbine on the VolturnUS-S semisubmersible platform. While the two methods produce broadly comparable results, the coupled wave loading on platforms within the farm results in wave force cancellations and amplifications that decrease multiline force directional ranges and increase multiline force extreme values (up to 7%) and standard deviations (up to 11%) for wave-driven load cases. The inclusion of wakes in FAST.Farm also reduces the net load on the shared anchor due to the velocity deficit, leading to larger differences between OpenFAST and FAST.Farm (up to 3% difference in mean loads) for load cases with operational turbines.

17 WIND ENERGY↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗

Differences in cluster and internal wake effects from mesoscale and large-eddy simulations off the US East Coast

Mesoscale simulations are increasingly used to estimate wake effects within and between large wind farms, despite limited validation for large-scale wake effects. This study evaluates the capabilities and limitations of mesoscale simulations in capturing wake-induced impacts on wind turbine power production through a direct comparison with large-domain large-eddy simulations (LESs) for three planned offshore wind farms under realistic atmospheric conditions and a range of atmospheric stabilities. We assess mesoscale performance in replicating wake characteristics behind single and multiple turbine clusters and quantify the resulting variability in mean turbine power. Results show that mesoscale Weather Research and Forecasting simulations with the Fitch wind farm parameterization capture key features of the velocity deficit downstream of both single and multiple wind farms, with mean root-mean-square errors near 5 % and good agreement with stability-driven wake behavior. However, in these simulations, the mesoscale Fitch parameterization underestimates power losses from internal wake effects, particularly when turbines align with the prevailing wind direction or under stable stratification. In these conditions, individual wakes persist and dominate downstream power deficits. The coarse resolution of the mesoscale simulations limits their ability to resolve individual wind turbine wakes that drive power fluctuations within wind farms. Nonetheless, mesoscale simulations can yield accurate estimates of combined wake losses from internal and cluster effects across some wind direction sectors, where errors in wake representation may cancel each other out. These findings underscore the strengths of mesoscale simulations for capturing broader wake patterns while highlighting their limitations for modeling turbine-level power losses. Future work should explore hybrid modeling approaches to capture both long-range cluster wake propagation and localized internal wake dynamics.

17 WIND ENERGY↗

Quantifying wind plant blockage under stable atmospheric conditions

Wind plant blockage reduces the wind velocity upstream undermining turbine performance for the first row of the plant. We assess how atmospheric stability modifies the induction zone of a wind plant in flat terrain. We also explore different approaches to quantifying the magnitude and extent of the induction zone from field-like observations. To investigate the influence from atmospheric stability, we compare simulations of two stable boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. We find a faster cooling rate at the surface, which produces a stronger stably stratified boundary layer, amplifies the induction zone of both an isolated turbine and of a large wind plant. A statistical analysis on the hub-height wind speed field shows wind slowdowns only extend far upstream (up to 15D) of a wind plant in strong stable boundary layers. To evaluate different ways of measuring wind plant blockage from field-like observations, we consider various ways of estimating the freestream velocity upstream of the plant. Sampling a large area upstream is the most accurate approach to estimating the freestream conditions, and thus of measuring the blockage effect. Also, the choice of sampling method may induce errors of the same order as the velocity deficit in the induction zone.

54 ENVIRONMENTAL SCIENCES↗

The revised FLORIDyn model: implementation of heterogeneous flow and the Gaussian wake

In this paper, a new version of the FLOw Redirection and Induction Dynamics (FLORIDyn) model is presented. The new model uses the three-dimensional parametric Gaussian FLORIS model and can provide dynamic wind farm simulations at a low computational cost under heterogeneous and changing wind conditions. Both FLORIS and FLORIDyn are parametric models which can be used to simulate wind farms, evaluate controller performance and can serve as a control-oriented model. One central element in which they differ is in their representation of flow dynamics: FLORIS neglects these and provides a computationally very cheap approximation of the mean wind farm flow. FLORIDyn defines a framework which utilizes this low computational cost of FLORIS to simulate basic wake dynamics. This is achieved by creating so-called observation points (OPs) at each time step at the rotor plane which inherit the turbine state. In this work, we develop the initial FLORIDyn framework further considering multiple aspects. The underlying FLORIS wake model is replaced by a Gaussian wake model. The distribution and characteristics of the OPs are adapted to account for the new parametric model but also to take complex flow conditions into account. To achieve this, a mathematical approach is developed to combine the parametric model and the changing, heterogeneous world conditions and link them with each OP. We also present a computationally lightweight wind field model to allow for a simulation environment in which heterogeneous flow conditions are possible. FLORIDyn is compared to Simulator for Offshore Wind Farm Applications (SOWFA) simulations in three- and nine-turbine cases under static and changing environmental conditions. The results show a good agreement with the timing of the impact of upstream state changes on downstream turbines. They also show a good agreement in terms of how wakes are displaced by wind direction changes and when the resulting velocity deficit is experienced by downstream turbines. A good fit of the mean generated power is ensured by the underlying FLORIS model. In the three-turbine case, FLORIDyn simulates 4 s simulation time in 24.49 ms computational time. The resulting new FLORIDyn model proves to be a computationally attractive and capable tool for model-based dynamic wind farm control.

17 WIND ENERGY↗

Investigating the interactions between wakes and floating wind turbines using FAST.Farm

As floating offshore wind progresses to commercial maturity, wake and array effects across a farm of floating offshore wind turbines (FOWTs) will become increasingly important. While wakes of land-based and bottom-fixed offshore wind turbines have been extensively studied, only recently has this topic become relevant for floating turbines. This work presents an investigation of the mutual interaction between the motions of floating wind turbines and wakes using FAST.Farm. While FAST.Farm has been extensively validated across a wide range of conditions, it has never been validated for FOWT applications. Hence, in the first part of this work, we validate FAST.Farm by comparing simulations of a single FOWT against high-fidelity results from large-eddy simulations available in the literature. The validation is based on wake meandering, mean wake deflection, and velocity deficit at different downstream locations. This validation showed that the original axisymmetric (polar) wake model of FAST.Farm overpredicts the vertical wake deflection induced by shaft tilt and floater pitch, while the new curled wake model is capable of properly capturing the vertical wake deflection. In the second part, we use FAST.Farm to analyze a small three-unit array of FOWTs with a spacing of 7 diameters across a wide range of environmental conditions. The same National Renewable Energy Laboratory 5 MW reference wind turbine atop the OC4-DeepCwind semisubmersible is adopted for the three FOWTs and for the validation against high-fidelity simulations. To assess the effect of the floating substructure, we compare the power production, tower-base moments, and blade-root moments obtained for the floating turbines with the results obtained in a fixed-bottom configuration. The main differences introduced by the floating substructure are the motions induced by the waves, the change in the natural frequencies of the tower caused by differences in the boundary condition at its base, and the larger vertical deflection of the wake deficit due to the mean pitch of the platform. The impact of these differences, as well as other minor effects, are analyzed in detail.

16 TIDAL AND WAVE POWER↗

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗

A Study of Wake Characteristics of Marine Turbine Arrays: Preprint

As the marine energy industry continues to develop and installations grow from single turbines to small arrays and larger farms, understanding and predicting wake behavior using modeling tools is necessary for array design and optimization. The wake characteristics of upstream turbines, including velocity deficit, wake swirling, tip vortices, wake shape and direction, and flow recovery, influence the inflow conditions experienced by downstream turbines in an array and, subsequently, affect the individual turbine loading and power generated by the array. These wake characteristics are influenced by many factors, including site geography, inflow shear profile, turbulence, waves, the strength of the current resource, gravitational forces, storm surges, and seasonal effects like temperature and wind conditions. Efficient array modeling tools that capture these wake effects and resulting turbine-to-turbine interactions need to be developed to optimize farm layouts and design individual turbines for operation within arrays. Unlike the wind industry, the marine energy research community is largely limited to the use of either numerical tools or tank/flume testing to understand marine turbine wake interactions. This paper will investigate recent advances in numerical and experimental marine turbine wake modeling, with a focus on recent publications and on both axial-flow and cross-flow marine turbines. A summary of recent research on these topics is provided. Additionally, a demonstration of turbine wake interactions using computational fluid dynamics is presented.

marine hydrokinetic turbine↗

Soil moisture thresholds explain a shift from light-limited to water-limited sap velocity in the Central Amazon during the 2015–16 El Niño drought

Transpiration is often considered to be light- but not water-limited in humid tropical rainforests due to abundant soil water, even during the dry seasons. The record-breaking 2015–16 El Niño drought provided a unique opportunity to examine whether transpiration is constrained by water under severe lack of rainfall. We measured sap velocity, soil water content, and meteorological variables in an old-growth upland forest in the Central Amazon throughout the 2015–16 drought. We found a rapid decline in sap velocity (-38 ± 21%, mean ± SD.) and in its temporal variability (-88%) during the drought compared to the wet season. Such changes were accompanied by a marked decline in soil moisture and an increase in temperature and vapor pressure deficit. Sap velocity was largely limited by net radiation during the wet and normal dry seasons; however, it shifted to be primarily limited by soil moisture during the drought. The threshold in which sap velocity became dominated by soil moisture was at 0.33 m 3 m -3 (around -150 kPa in soil matric potential), below which sap velocity dropped steeply. Our study provides evidence for a soil water threshold on transpiration in a moist tropical forest, suggesting a shift from light limitation to water limitation under future climate characterized by increased temperature and an increased frequency, intensity, duration and extent of extreme drought events.

54 ENVIRONMENTAL SCIENCES↗

Iron Bonding with Light Elements: Implications for Planetary Cores Beyond the Binary System

Light element alloying in iron is required to explain density deficit and seismic wave velocities in Earth’s core. However, the light element composition of the Earth’s core seems hard to constrain as nearly all light element alloying would reduce the density and sound velocity (elastic moduli). The alloying light elements include oxidizing elements like oxygen and sulfur and reducing elements like hydrogen and carbon, yet their chemical effects in the alloy system are less discussed. Moreover, Fe-X-ray Absorption Near Edge Structure (Fe-XANES) fingerprints have been studied for silicate materials with ferrous and ferric ions, while not many X-ray absorption spectroscopy (XAS) studies have focused on iron alloys, especially at high pressures. To investigate the bonding nature of iron alloys in planetary interiors, we presented X-ray absorption spectroscopy of iron–nitrogen and iron–carbon alloys at high pressures up to 50 GPa. Together with existing literature on iron–carbon, –hydrogen alloys, we analyzed their edge positions and found no significant difference in the degree of oxidation among these alloys. Pressure effects on edge positions were also found negligible. Our theoretical simulation of the valence state of iron, alloyed with S, C, O, N, and P also showed nearly unchanged behavior under pressures up to 300 GPa. This finding indicates that the high pressure bonding of iron alloyed with light elements closely resembles bonding at the ambient conditions. We suggest that the chemical properties of light elements constrain which ones can coexist within iron alloys.

Advanced Photon Source↗

FLOW Estimation and Rose Superposition (FLOWERS): an integral approach to engineering wake models

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each discrete wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions and speeds. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed 2 orders of magnitude faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a relatively smooth solution space for the gradient-based optimizer to excel in comparison to the existing weighted-averaging power calculation.

17 WIND ENERGY↗

Monitoring Terrestrial Water Storage, Drought and Seasonal Changes in Central Oklahoma With Ambient Seismic Noise

Abstract Significant imbalances in terrestrial water storage (TWS) and severe drought have been observed around the world as a consequence of climate changes. Improving our ability to monitor TWS and drought is critical for water‐resource management and water‐deficit estimation. We use continuous seismic ambient noise to monitor temporal evolution of near‐surface seismic velocity, dv / v , in central Oklahoma from 2013 to 2022. The derived dv / v is found to be negatively correlated with gravitational measurements and groundwater depths, showing the impact of groundwater storage on seismic velocities. The hydrological effects involving droughts and recharge of groundwater occur on a multi‐year time scale and dominate the overall derived velocity changes. The thermoelastic response to atmospheric temperature variations occurs primarily on a yearly timescale and dominates the superposed seasonal velocity changes in this study. The occurrences of droughts appear simultaneously with local peaks of dv / v , demonstrating the sensitivity of near‐surface seismic velocities to droughts.

58 GEOSCIENCES↗

Reduced Surface Recombination Velocity in Mixed-Cation Mixed-Halide Hybrid Perovskite Devices Leads to Improved Open-Circuit Voltage

Surface recombination plays an important role in perovskite solar cell performance.Despite rapid progress, the record power conversion efficiency (PCE, ~ 25.7%) for perovskite single-junction solar cells still lags behind its theoretical limit (~ 32%), largely due to the open circuit voltage (VOC) deficit that results from non-radiative recombination at the bulk of perovskite films and the interfaces between perovskite and transport layers. In this work, we show that (3-aminopropyl)trimethoxysilane (APTMS) surface passivation can reduce the surface recombination velocities (SRVs) of perovskite solar cells, thus improving device VOC and PCE.

perovskite solar cell, surface recombination, surf↗

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗

Thermoelastic Anomaly of Iron Carbonitride Across the Spin Transition and Implications for Planetary Cores

Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe-C-N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal-structured Fe 7 (N 0.75 C 0.25 ) 3 up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First-principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe 7 (N 0.75 C 0.25 ) 3 . Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe-C-N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)-rich planetary cores, thereby influencing the dynamics of such cores.

Geology↗

Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part II: Dissolving Stage

A minimalistic analytical model allowing analysis of the dissolving stage of nonprecipitating convective clouds is proposed. The model takes into account two mechanisms: turbulent mixing with a dry environment and cloud volume settling. The temporal changes in the spatial structure of a cloud and in its immediate environment in the course of cloud dissolving are analyzed. The comparison of the effects of a temperature increase in the course of cloud descent and mixing with dry surrounding air shows that the descent is a dominating factor determining a decrease in the liquid water content (LWC), while mixing has a stronger effect on the cloud shape. Narrowing/broadening of clouds due to lateral mixing with dry air during cloud dissolving is determined by the potential evaporation parameter proportional to the ratio of the saturation deficit in the cloud environment to LWC. An equation for cloud dissolving time is obtained. After a cloud totally dissolves, it leaves behind an area with humidity exceeding that of the environment. The main parameter determining the dissolving time is the downdraft velocity. It should exceed 50 cm s -1 in order to provide reasonable dissolving time. The turbulent intensity, LWC, and humidity of the environment air also have an impact on dissolving time: the lower the LWC and the humidity of environment air, the faster cloud dissolving is. The simple solution presented in this paper can be useful for evaluation of cloud characteristics at the dissolving stage and can be included in procedures of parameterization of cloud cover formed by nonprecipitating or slightly precipitating cumulus clouds (Cu). Values of the environment humidity and temperature, LWC at cloud top, cloud width, vertical velocity of downdraft, and the turbulent coefficient should be parameters of this parameterization.

54 ENVIRONMENTAL SCIENCES↗