Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sensible flux”

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 OSTI2026

Monthly Mean In Situ Surface Flux Observations Paired with Satellite-Derived and Reanalysis-Based Flux Data for the Great Lakes Region, 2001–2020

Surface radiative and turbulent heat fluxes over the Great Lakes strongly influence regional hydrological and meteorological processes, and their accurate representation is critical for numerical weather prediction and coupled atmosphere–lake modeling. However, direct flux observations are spatially sparse across the region, so gridded reanalysis and satellite-derived products are often used for climatological analyses and model evaluation despite differences in their flux representations. This dataset provides processed, quality-controlled, monthly mean surface flux observations from the Great Lakes Evaporation Network (GLEN), AmeriFlux, and the National Data Buoy Center, paired with spatiotemporally matched flux estimates from two reanalysis products, the fifth generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset (ERA5) and the Modern Era Reanalysis for Research and Applications, version 2 (MERRA-2), and two satellite-derived products, the Clouds and Earth's Radiant Energy Systems Energy Balanced and Filled (CERES-EBAF) and the Cloud, Albedo and Surface Radiation dataset from AVHRR data - Edition 3 (CLARA-A3). The dataset includes sixteen observational stations with variable temporal coverage within 2001–2020. For each station, a CSV file contains monthly time series of available flux variables, including surface downwelling shortwave radiation (SW), surface downwelling longwave radiation (LW), sensible heat (SH) flux, and latent heat flux (LH), alongside matched gridded product values where available. Columns in the CSV file correspond to different variables sourced from each dataset, with column titles structured as "{dataset}_{variable}". Columns with relevant metadata are also provided in each CSV file, including station latitude and longitude, monthly timestamps, and the name of the sourced observational data. These files are structured for direct use in common analysis tools, including Microsoft Excel, Python pandas, and Python matplotlib. This dataset supports climatological analysis of the Great Lakes regional surface energy budget, evaluation of satellite-derived and reanalysis-based flux products, and development or validation of flux representations in numerical weather prediction and coupled atmosphere–lake models.

Great Lakes↗
DOE OSTI2021

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part I: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LESs. We derive these quantities from measurements and reanalysis obtained for two cold-air outbreak (CAO) events during Phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) in February–March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from ERA5 data agree reasonably well with those derived from ACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time-evolving heat fluxes, wind, and advective tendencies profiles from ERA5 data to drive the LES. We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LESs. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path simulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol–cloud–meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Geospatial coherence of surface-atmosphere fluxes in the upper Great Lakes region

Surface-atmosphere fluxes are known to vary at multiple time scales, but uncertainty is high as to how fluxes change spatially within regions. With an increase in the number of eddy covariance towers, we are now able to examine the geospatial coherence of ecosystem fluxes, using time-series correlation. Eighteen sites from Michigan and Wisconsin were used in this study, ranging from 100 m to 600 km apart. Surface-atmosphere fluxes from a six-month period were used to quantify spatial coherence on a pair-wise basis. Using geospatial statistics, carbon and sensible heat (H) fluxes were found to be 95% correlated directly outside of their flux footprint and 56% correlated up to a distance of ~35 km. Latent (LE) and momentum (τ) fluxes were less correlated, 83% directly outside of their flux footprint and 40% at a distance of ~130 km albeit, at a much larger spatial distance than for the carbon and sensible heat fluxes. All fluxes showed strong spectral resonance at diel and seasonal timescales, with 1-, 2- and 3-month periods being common modes of variability among H, LE, and τ fluxes. Results based on Empirical Orthogonal Function show distinct transitions of net ecosystem exchange from fall to winter before photosynthesis or respiration while H and τ do not exhibit coherent trends. Here, this work demonstrates the potential of quantifying geospatial coherence of surface-atmosphere fluxes in the Midwestern United States, with the ability to predict fluxes beyond the spatial limit of a single flux tower footprint. Ultimately, expanding the flux measurements to larger scales would allow better spatial scaling of terrestrial surface-atmosphere fluxes between tower footprint and modeling or remote sensing scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Disentangling the Relative Drivers of Seasonal Evapotranspiration Across a Continental-Scale Aridity Gradient

Evapotranspiration (ET) is a significant ecosystem flux, governing the partitioning of energy at the land surface. Understanding the seasonal pattern and magnitude of ET is critical for anticipating a range of ecosystem impacts, including drought, heat-wave events, and plant mortality. Here, in this study, we identified the relative controls of seasonal variability in ET, and how these controls vary among ecosystems. We used overlapping AmeriFlux and PhenoCam time series at a daily timestep from 20 sites to explore these linkages (# site-years >100), and our study area covered a broad climatological aridity gradient in the U.S. and Canada. We focused on disentangling the most important controls of bulk surface conductance (G s ) and evaporative fraction (EF = LE/[H + LE]), where LE and H represent latent and sensible heat fluxes, respectively. Specifically, we investigated how vegetation phenology varied in importance relative to meteorological variables (vapor pressure deficit and antecedent precipitation) as a driver of G s and EF using path analysis, a framework for quantifying and comparing the causal linkages among multiple response and explanatory variables. Our results revealed that the drivers of G s and EF seasonality varied significantly between energy- and water-limited ecosystems. Specifically, precipitation had a much higher effect in water-limited ecosystems, while seasonal patterns in canopy greenness emerged as a stronger control in energy-limited ecosystems. Given that phenology is expected to shift under future climate, our findings provide key information for understanding and predicting how phenology may impact 21st-century hydroclimate regimes and the surface-energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Global mapping of urban thermal anisotropy reveals substantial potential biases for remotely sensed urban climates

Urban thermal anisotropy (UTA) drastically biases satellite-derived urban surface temperatures and fluxes, and consequently inhibits a more comprehensive understanding of global urban climates. However, UTA patterns and associated biases in observed urban climate variables have not been investigated across an adequate number of global cities with diverse contexts; nor is it known whether there are globally measurable factors that are closely related to the UTA intensity (UTAI, quantified as the maximum difference between the nadir and off-nadir urban thermal radiation). Here we investigate the UTAI over more than 5500 cities worldwide using the moderate resolution imaging spectroradiometer (MODIS) multi-angle land surface temperature (LST) observations from 2003 to 2021. The results show that the global mean UTAI can reach 5.1, 2.7, 2.4, and 1.7 K during summer daytime, winter daytime, summer nighttime, and winter nighttime, respectively. Using nadir LST observations as a reference, our analysis reveals that UTA can lead to an underestimation of satellite-based urban surface sensible heat fluxes ( H ) by 45.4% and surface urban heat island intensity ( I s ) by 43.0% when using LST observations obtained from sensor viewing zenith angles (VZAs) of ±60°. However, practitioners can limit the biases of H and I s within ±10% by using LSTs from sensor VZAs within ±30°. We find that UTAI is closely related to urban impervious surface percentage and surface air temperature across global cities. In conclusion, these findings have implications for angular normalization of satellite-retrieved instantaneous LST observations across cities worldwide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic changes in evapotranspiration, canopy photosynthesis and expansive growth in open field under rapid fluctuating radiation

Abstract In open fields, plants experience dynamic changes in environment, particularly radiation, temperature, wind, and humidity but their short-term responses have not been adequately characterized under natural conditions. In this study, we assessed causal effects of rapid radiation fluctuation on seven plant parameters in open fields of cotton and sweet corn canopies. The parameters are evapotranspiration ( ET ), canopy photosynthesis indicated by CO2 flux from air to canopy ( FCO2 ), sensible heat energy flux ( H ), canopy conductance for water vapor ( gcw ), canopy surface temperature ( Ts ), shoot and leaf elongation rate, and stem diameter change. The energy and CO2 fluxes were measured with Bowen ratio/energy balance/CO2 gradient (BREB+) technique, using averaging time of 5 min. Shoot + leaf elongation and stem diameter change were monitored with position transducers using averaging time of 1 min. All parameters were all found to respond to change in radiation and transpiration within minutes or sooner. While radiation effects on canopy gas exchanges are expected, illuminating are the indirect but immediate effects on shoot + leaf growth and stem diameter change through radiation effects on transpiration and plant water status. A novel finding is that gcw also responded within minutes or sooner to radiation fluctuations and that FCO2 was related almost linearly to gcw . Results are discussed in terms of soil-plant-atmosphere continuum, and interpreted in terms of dynamic interactions between transpiration and plant water status. The clear inverse relationship between ET and elongation rate or stem diameter changes provides additional evidence supporting the validity of 5-min averaging for the BREB + technique.

Hsiao, Theodore C.↗
DOE OSTI2021

Analysis of Coastal Fog from a Ship During the C-FOG Campaign

This work presents ship-based measurements of fog off St John’s Newfoundland on 13 September 2018 during the Coastal Fog (C-FOG) field campaign. The measurements included cloud-particle spectra, cloud-base height and aerosol backscatter, radiation, turbulence, visibility, and sea surface temperature. Radiosonde soundings were made at intervals less than 2 hours. Fog occurred in two episodes under the center of eastward-moving synoptic low-pressure system. The boundary layer structure during the first fog episode consisted of 3 layers, separated by 2 saturated air temperature inversions and capped by a subsidence inversion. The lowest layer was fog and the upper layers were cloud. The second fog episode consisted of one well-mixed fog layer capped by a subsidence inversion. Low wind speeds and stable stratification maintained low surface-layer turbulence during fog. Droplet size distributions had typical bimodal distributions. The visibility correlated with the droplet number concentration and liquid water content. The air temperature was higher than the sea surface temperature for the first 30 minutes of the first fog episode but was colder than the sea for the remainder of all fog. The sensible heat flux was upward, from sea to air, for the first 62 % of the first fog episode and then reversed to downward, from air to sea, for the remainder of the first fog episode and the second fog episode. Here, the counter-gradient heat fluxes observed (i.e., opposite to what is expected from the instantaneous air-sea temperature difference) appear to be related to turbulence, entrainment and stratification in the fog layer that overwhelmed the influence of air-sea temperature difference. While the synoptic-scale dynamics preconditioned the area for fog formation, the final step of fog appearance in this case was nuanced by stratification-turbulence interactions, local advective processes and microphysical environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Lower Tropospheric Processes: A Control on the Global Mean Precipitation Rate

The spread in global mean precipitation among climate models is explored in two ensembles using the complementary perspectives of surface evaporation and energy budgets. Models with higher global mean precipitation have stronger oceanic evaporation, driven by drier near-surface air. The drier surface conditions occur alongside increases in near-surface temperature and moisture at 925 hPa, which point to stronger boundary layer mixing. Correlations suggest that the degree of lower tropospheric mixing explains 18%–49% of the intermodel precipitation variance. To test this hypothesis, the degree of mixing is indirectly varied in a single-model experiment by adjusting the relative humidity threshold that controls low-cloud fraction. Indeed, increasing lower tropospheric mixing results in more global mean precipitation. Energetically, increased precipitation rates are associated with more downwelling longwave radiation to the surface and weaker sensible heat fluxes. These results highlight how lower-tropospheric processes must be better constrained to reduce the precipitation discrepancy among climate models.

58 GEOSCIENCES↗
DOE OSTI2020

Characterization of Surface Heterogeneity-Induced Convection Using Cluster Analysis

Land-atmosphere interactions and boundary layer processes often control the formation of shallow clouds and subsequently deep convective precipitation over the Southern Great Plains. In this study, we examine the impacts of large-scale advection on the cloud populations and land-atmospheric coupling observed during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign in 2016. In this work, we performed two Large Eddy Simulations (LES) using the Weather Research and Forecasting Model for a day in which the transition from clear-sky to shallow precipitating clouds and isolated deep convection was observed. The control simulation reproduced the overall distribution of cloud populations by using realistic soil conditions with an interactive land model. In the sensitivity simulation where large-scale advection is removed, a strong relationship between the land and boundary layer is found. To study the timing, location, and intensity of convective initiation and the relationship of clouds with land surface properties, a cluster analysis of equivalent potential temperature is performed for the simulation without large-scale advection. That analysis shows that convective clouds first form over regions with higher surface sensible heat flux. Precipitation from those convective clouds likely triggers new updrafts nearby about two hours later through the lifting associated with cold pools. Furthermore, the cluster analysis also shows that in addition to the spatial pattern of soil moisture, land use and soil texture in western Oklahoma also influence the location of convective initiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interpreting the Diurnal Cycle of Clouds and Precipitation in the ARM GoAmazon Observations: Shallow to Deep Convection Transition

The Green Ocean Amazon (GoAmazon) 2014/5 field campaign data are used to study the diurnal cycle of clouds and precipitation. Through a careful classification of days with shallow cumulus, congestus and deep convection, we investigate the major differences among locally generated convection regimes and the most important environmental factors governing the shallow-to-deep convection transition. On shallow cumulus days, a greater sensible heat flux drives deeper boundary layer growth, which entrains drier free-tropospheric air and lowers the relative humidity, thus leading to a significantly higher cloud base than those on days with deeper convection. Congestus and deep convection regimes exhibit distinct cloud top height distributions with noticeable differences in the vertical wind shear in the mid-troposphere, suggesting an important role of wind shear in limiting the vertical extent of convection. On deep convection days, with preexisting nocturnal convection or cold-pools from external disturbances, the timing of peak surface precipitation (12:00–13:00 LST) tends to be in-phase with the diurnal variation in surface fluxes. However, it takes longer for local deep convection to develop without these disturbances. A plume model with thermodynamic and dynamical constraints is developed to explore the relative importance of various convection-controlling factors. Initial cloud-base vertical velocity and buoyancy are important in helping parcels ascend to the level of free convection (LFC). After parcels reach the LFC, entrainment of environmental air and lower free troposphere humidity become crucial in determining cloud top. Entrainment rate differentiates among convection regimes, which may be tied to the cloud size distribution at cloud base.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impacts of Degradation on Water, Energy, and Carbon Cycling of the Amazon Tropical Forests

Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high-resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED-2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED-2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Winter Warming in North America Induced by Urbanization in China

Rapid urbanization in China in the past several decades has been shown to affect surface air temperature (SAT) locally, but its global impact is unclear. Using a global climate model, we have investigated the potential impact of urbanization in China on SAT at the global scale. The expansion of urban areas in China from 1985 to 2017 leads to winter surface air warming by as much as 0.8°C in North America (NA), largely caused by warm temperature advection due to global-scale atmospheric circulation changes. The East Asian jet stream moves northward in eastern and northeastern China where the expansion of urban areas has increased SAT and sensible heat flux, forcing change of the atmospheric stationary waves and so the structure of the jet streams.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A New Post‐Hoc Method to Reduce the Energy Imbalance in Eddy Covariance Measurements

Abstract Latent and sensible heat flux observations are essential for understanding land–atmosphere interactions. Measurements from the eddy covariance technique are widely used but suffer from systematic energy imbalance problems, partly due to missing large eddies from sub‐mesoscale transport. Because available energy drives the development of large eddies, we propose an available energy based correction method for turbulent flux measurements. We apply our method to 172 flux tower sites and show that we can reduce the energy imbalance from −14.99 to −0.65 W m −2 on average, together with improved consistency between turbulent fluxes and available energy and associated increases inr 2 at individual sites and across networks. Our results suggest that our method is conceptually and empirically preferable over the method implemented in the ONEFlux processing. This can contribute to the efforts in understanding and addressing the energy imbalance issue, which is crucial for the evaluation and calibration of land surface models.

Geology↗
DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulating Pyrocumulonimbus Clouds Using a Multiscale Wildfire Simulation Framework

Pyrocumulonimbus (pyroCb) clouds, driven by extreme fires under favorable meteorological conditions, can inject smoke into the stratosphere at magnitudes comparable to those of moderate volcanic eruptions, potentially altering the global radiative balance and atmospheric composition. However, simulating pyroCb is particularly challenging in Earth system models. Using the Energy Exascale Earth System Model (E3SM), we developed a novel global multiscale framework to model pyroCb events in California, which includes a high‐resolution fire radiative power time series, a one‐dimensional plume‐rise parameterization, a fire‐induced vertical water vapor transport scheme, and a surface wildfire sensible heat flux representation. Our simulation successfully reproduces many pyroCb features, including cloud height, spatiotemporal evolution, and convective intensity in comparison with satellite and ground‐based observations. Sensitivity experiments show that realistic pyroCb simulation depends on vertical water vapor transport. These advances provide a basis for future exploration of pyroCb impacts at regional and global scales within climate models.

E3SM↗
DOE OSTI2025

Influence of Soil Moisture on the Development of Organized Convective Systems in South America

Abstract This study investigates the role of soil moisture (SM) on the initiation and organization of convective systems using the convection‐permitting ICOsahedral Non‐hydrostatic (ICON) model. We conduct two sets of experiments: a Control experiment with interactive SM and a fixed SM experiment (FixedSM) with invariable SM conditions. We focus on two regions in South America: the Amazon and southeastern South America (SESA). Larger organized convective systems are associated with greater SM heterogeneity in both regions, though other large‐scale synoptic influences affect the robustness of this relationship in SESA. These results remain largely unaffected by disabling the effects of precipitation on SM in the FixedSM experiment, and complementary analyses using satellite‐based estimates of SM and precipitation support these findings. Spatial compositing of mesoscale environments in the Amazon shows the presence of well‐defined SM gradients, at a length scale of a few hundred kilometers, many hours before convective system detection. Larger SM gradients correspond to larger gradients in thermodynamic variables, particularly surface temperature and sensible heat flux, and are associated with larger convective systems. Overall, our findings suggest that surface heterogeneities such as SM gradients not only affect deep convection initiation, as previously suggested, but they can also encourage the growth and organization of convective systems into larger clusters, particularly in the absence of significant synoptic influences.

54 ENVIRONMENTAL SCIENCES↗