Search NASA⌕ Search

SEARCH · Search NASA

Results for “Land surface”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Assessing the Consistency of Estimated Ground Cover Fractions between the BLM AIM Method and Optical Remote Sensing Method

Since 2012, Argonne National Laboratory (Argonne) has supported the Bureau of Land Management (BLM) in developing remote sensing methodologies for long-term environmental monitoring of Palo Verde Mesa in eastern Riverside County, California, including methods for: detailed mapping of ephemeral streams, estimating fractional cover of desert-land surface components (e.g., trees, shrubs, litters, and bare ground), evaluating erosion risk or land stability, and characterizing vegetation alliances using spatial structure and geostatistical approaches. These studies showed the promise of remote sensing for monitoring changes in desert landscapes by providing information that would be difficult to obtain through field surveys. During this time the BLM has also worked to establish long-term monitoring protocols and compiled field-observation data collected using standardized protocols from the Assessment, Inventory, and Monitoring (AIM) strategy. The AIM data can be compared to data derived from remote sensing methods to evaluate their relative operational utility in monitoring landscape change. If ground cover estimated using remotely sensed imagery is comparable to AIM ground cover estimates, then remote sensing can be used to monitor whether any land cover change in desert landscapes may be related to solar energy development. Therefore, the goal of this study was to determine the consistency in ground cover estimates between AIM data and those derived from publicly- available remotely sensed imagery, such as that available through the U.S. Department of Agriculture, National Agricultural Imagery Program (NAIP), to examine feasibility of a remote sensing method for complementing AIM monitoring. Based on the image analysis in this study, we also provide recommendations for how small unmanned aerial system (sUAS) data may be used to complement BLM’s AIM data and NAIP imagery for future vegetation monitoring. The ground cover types we originally planned to investigate were trees, shrubs, and bare ground. However, the small sample size and a limited range of cover fraction of trees and shrubs in the AIM dataset (e.g., 33 samples with a maximum shrub cover of 18%, 14 samples with a maximum tree cover of 17%) did not allow for performing a meaningful evaluation for the remote sensing approach. Therefore, we conducted the study focusing on bare ground, foliar, and rock cover, all of which are indicators reported in the AIM remote sensing dataset.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

First Assessment of Cloud‐Land Coupling in LASSO Large‐Eddy Simulations

Abstract To enhance our understanding of cloud simulations over land, this study provides the first assessment of coupling between cloud and land surface in the Large‐Eddy Simulation (LES) Atmospheric Radiation Measurement Symbiotic Simulation and Observation (LASSO) activity for the shallow convection scenario. The analysis of observation data reveals a diurnal cycle of cloud‐land coupling, which co‐varies with surface fluxes. However, coupled (or decoupled) cumulus clouds are inadequately simulated, manifesting as a too‐high (or low) occurrence frequency during the afternoon. This discrepancy is mirrored by the overestimated cloud liquid water path and cloud‐top height. These overestimations are linked to the overpredicted boundary‐layer development and the easier trigger of shallow convection misrepresented in LES runs. Our study underscores the need to improve the representations of boundary‐layer processes and cloud‐land interactions within LES to better simulate shallow clouds in the future.

58 GEOSCIENCES↗
DOE OSTI2024

Assessment of U.S. Urban Surface Temperature Using GOES-16 and GOES-17 Data: Urban Heat Island and Temperature Inequality

Abstract This study utilizes hourly land surface temperature (LST) data from the Geostationary Operational Environmental Satellite (GOES) to analyze the seasonal and diurnal characteristics of surface urban heat island intensity (SUHII) across 120 largest U.S. cities and their surroundings. Distinct patterns emerge in the classification of seasonal daytime SUHII and nighttime SUHII. Specifically, the enhanced vegetation index (EVI) and albedo (ALB) play pivotal roles in influencing these temperature variations. The diurnal cycle of SUHII further reveals different trends, suggesting that climate conditions, urban and nonurban land covers, and anthropogenic activities during nighttime hours affect SUHII peaks. Exploring intracity LST dynamics, the study reveals a significant correlation between urban intensity (UI) and LST, with LST rising as UI increases. Notably, populations identified as more vulnerable by the social vulnerability index (SVI) are found in high UI regions. This results in discernible LST inequality, where the more vulnerable communities are under higher LST conditions, possibly leading to higher heat exposure. This comprehensive study accentuates the significance of tailoring city-specific climate change mitigation strategies, illuminating LST variations and their intertwined societal implications.

Environmental Sciences & Ecology↗
DOE OSTI2021

Land cover change-induced decline in terrestrial gross primary production over the conterminous United States from 2001 to 2016

As one of the most dynamic aspects of global environmental change, land cover change (LCC) has a profound impact on terrestrial carbon sequestration. However, LCC-induced carbon fluxes are still the most uncertain terms in global and regional carbon budgets. Ecosystem gross primary production (GPP) is the total carbon uptake by vegetation through photosynthesis, serving as a major control on ecosystem function and land carbon balance during and after the modification of the land surface. However, accurately capturing LCC-induced GPP changes requires both high-quality land cover data and controlling for variation driven by other environmental factors such as climate. In this study, we comprehensively examined the effects of LCC on annual GPP trends over the conterminous United States (CONUS) from 2001 to 2016 using the USGS National Land Cover Database, a remote sensing-driven ecosystem model, and the Google Earth Engine cloud computing platform. We designed a series of model experiments to identify LCC effects on GPP by controlling climate effects. During the study period, LCC exerted a strong negative effect on total GPP across the CONUS ([-2.2, -1.8] Tg C yr -2 ), while climate had smaller positive effects ([0.17, 0. 92] Tg C yr -2 ). The LCC-induced reduction of GPP was mainly caused by net forest loss ([-1.98, -1.39] Tg C yr -2 ) and urban expansion ([-2.03, -1.92] Tg C yr -2 ), but was partially offset by increases in crop area ([+0.66, +0.79] Tg C yr -2 ). Ensemble simulations from TRENDY did not capture the strong negative LCC influences on GPP, likely due to limitations of the adopted land use/cover data. Overall, our study provides a novel perspective on LCC-induced GPP changes, which could help to improve our understanding of ecosystem function changes and constrain the estimation of land carbon balance in the context of anthropogenic activity and climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Enlarged Nonclosure of Surface Energy Balance With Increasing Atmospheric Instabilities Linked to Changes in Coherent Structures

The nonclosure of surface energy balance remains an outstanding problem in eddy covariance (EC) measurements of land-surface fluxes of heat, water vapor, and CO 2 . Our data collected from an EC tower over a semiarid sagebrush ecosystem indicate that under unstable atmospheric conditions, the nonclosure becomes increased with increasing instability, consistent with many other studies. It is demonstrated that the increased nonclosure is not caused by the inadequate sampling of large-scale turbulent motions using a 30-min averaging interval, though the scales of turbulent motions dominating sensible and latent heat fluxes become enlarged with increasing instability. Quadrant analysis is then used to reveal that the flux contributions from ejections remain nearly constant with increasing instability, whereas the flux contributions from sweeps are reduced and their time fractions increase. Our results imply that the increased nonclosure of surface energy balance is associated with changes in turbulent structures including their dominant time scales and flux contributions of ejections and sweeps as the atmospheric instability increases, which require further studies using vertically distributed observations and/or large eddy simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Diversity and Composition of Methanotroph Communities in Caves

Recent observations have shown the atmospheric greenhouse gas methane (CH 4 ) is consumed by microorganisms (methanotrophs) in caves at rates comparable to CH 4 oxidation in surface soils. Caves are abundant in karst landscapes that comprise 14% of Earth’s land surface area, and therefore may represent a potentially important, but overlooked, CH 4 sink.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Evaluation and development of surface layer scheme representation of temperature inversions over boreal forests in Arctic wintertime conditions

Abstract. In this study, the Noah land surface model used in conjunction with the Mellor–Yamada–Janjić surface layer scheme (hereafter, Noah-MYJ) and the Noah multiphysics scheme (Noah-MP) from the Weather Research and Forecasting (WRF) 4.5.1 mesoscale model are evaluated with regard to their performance in reproducing positive temperature gradients over forested areas in the Arctic winter. First, simplified versions of the WRF schemes, recoded in Python, are compared with conceptual models of the surface layer in order to gain insight into the dependence of the temperature gradient on the wind speed at the top of the surface layer. It is shown that the WRF schemes place strong limits on the turbulent collapse, leading to lower surface temperature gradient at low wind speeds than in the conceptual models. We implemented modifications to the WRF schemes to correct this effect. The original and modified versions of Noah-MYJ and Noah-MP are then evaluated compared to long-term measurements at the Ameriflux Poker Flat Research Range, a forest site in interior Alaska. Noah-MP is found to perform better than Noah-MYJ because the former is a two-layer model which explicitly takes into account the effect of the forest canopy. Indeed, a non-negligible temperature gradient is maintained below the canopy at high wind speeds, leading to overall larger gradients than in the absence of vegetation. Furthermore, the modified versions are found to perform better than the original versions of each scheme because they better reproduce strong temperature gradients at low wind speeds.

Geology↗
DOE OSTI2022

Simultaneous Observations of Surface Layer Profiles of Humidity, Temperature, and Wind Using Scanning Lidar Instruments

We demonstrate the combination of three Doppler lidars (DLs), a water vapor differential absorption lidar and a temperature rotational Raman lidar for the investigation of the interactions between the land-surface and the atmospheric boundary layer. This combination of scanning lidars was operated for the first time during the Land-Atmosphere Feedback Experiment at the Atmospheric Radiation Measurement program's Southern Great Plain site, Oklahoma, USA, in August 2017, and provided simultaneous surface layer profiles of horizontal wind, humidity and temperature. The horizontal wind profiles were determined using the dual-Doppler method with two DLs. The scans were performed above four towers providing atmospheric variables and fluxes. These combined lidar data allowed for the estimation of the friction velocity as well as the surface latent and sensible heat fluxes. For this purpose, profiles calculated with Monin-Obukhov similarity theory were fitted with the Levenberg-Marquardt nonlinear least squares curve-fitting method to the measured surface-layer lidar profiles. We present case studies of three 50-min periods to illustrate the new method. The measurements resulted in friction velocities of 0.31, 0.29, and 0.38 ms –1 , sensible surface heat fluxes of 312.1, 234.1, and 183.3 Wm –2 and latent surface heat fluxes of 251.8, 227.3, and 274.4 Wm –2 . These values were compared with the in-situ measurements of the towers. Considering the sampling differences and the error analyses, the results agree within 25.1%, 42.5%, and 28.8%, respectively, demonstrating that this synergy of scanning active remote sensing systems can be used to derive surface fluxes with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Albedo of crops as a nature-based climate solution to global warming

Abstract Surface albedo can affect the energy budget and subsequently cause localized warming or cooling of the climate. When we convert a substantial portion of lands to agriculture, land surface properties are consequently altered, including albedo. Through crop selection and management, one can increase crop albedo to obtain higher levels of localized cooling effects to mitigate global warming. Still, there is little understanding about how distinctive features of a cropping system may be responsible for elevated albedo and consequently for the cooling potential of cultivated lands. To address this pressing issue, we conducted seasonal measurements of surface reflectivity during five growing seasons on annual crops of corn-soybean–winter wheat ( Zea mays L.- Glycine max L. Merrill - Triticum aestivum L. ; CSW) rotations at three agronomic intensities, a monoculture of perennial switchgrass ( Panicum virgatum L. ), and perennial polycultures of early successional and restored prairie grasslands. 
We found that crop-species, agronomic intensity, seasonality, and plant phenology had significant effects on albedo. The mean±SD albedo was highest in perennial crops of switchgrass (0.179±0.04), intermediate in early successional crops (0.170±0.04), and lowest in a reduced input corn systems with cover crops (0.154±0.02). The strongest cooling potentials were found in soybean (-0.450 kg CO 2 e m -2 yr -1 ) and switchgrass (-0.367 kg CO 2 e m -2 yr -1 ), with up to -0.265 kg CO 2 e m -2 yr -1 of localized climate cooling annually provided by different agroecosystems. We also demonstrated how diverse ecosystems, leaf canopy, and agronomic practices can affect surface reflectivity and provide another potential nature-based solution for reducing global warming at localized scales.

Lei, Cheyenne (ORCID:0000000272746309)↗
DOE OSTI2023

Measuring surface albedo for various crops

Surface albedo can affect the energy budget and subsequently cause localized warming or cooling of the climate. When we convert a substantial portion of lands to agriculture, land surface properties are consequently altered, including albedo. Through crop selection and management, one can increase crop albedo to obtain higher levels of localized cooling effects to mitigate global warming. Still, there is little understanding about how distinctive features of a cropping system may be responsible for elevated albedo and consequently for the cooling potential of cultivated lands. To address this pressing issue, we conducted seasonal measurements of surface reflectivity during five growing seasons on annual crops of corn-soybean–winter wheat

albedo↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

A Century of Reforestation Reduced Anthropogenic Warming in the Eastern United States

Abstract Restoring and preserving the world's forests are promising natural pathways to mitigate some aspects of climate change. In addition to regulating atmospheric carbon dioxide concentrations, forests modify surface and near‐surface air temperatures through biophysical processes. In the eastern United States (EUS), widespread reforestation during the 20th century coincided with an anomalous lack of warming, raising questions about reforestation's contribution to local cooling and climate mitigation. Using new cross‐scale approaches and multiple independent sources of data, we uncovered links between reforestation and the response of both surface and air temperature in the EUS. Ground‐ and satellite‐based observations showed that EUS forests cool the land surface by 1–2°C annually compared to nearby grasslands and croplands, with the strongest cooling effect during midday in the growing season, when cooling is 2–5°C. Young forests (20–40 years) have the strongest cooling effect on surface temperature. Surface cooling extends to the near‐surface air, with forests reducing midday air temperature by up to 1°C compared to nearby non‐forests. Analyses of historical land cover and air temperature trends showed that the cooling benefits of reforestation extend across the landscape. Locations surrounded by reforestation were up to 1°C cooler than neighboring locations that did not undergo land cover change, and areas dominated by regrowing forests were associated with cooling temperature trends in much of the EUS. Our work indicates reforestation contributed to the historically slow pace of warming in the EUS, underscoring reforestation's potential as a local climate adaptation strategy in temperate regions.

Environmental Sciences & Ecology↗
DOE OSTI2020

A new framework to map fine resolution cropping intensity across the globe: Algorithm, validation, and implication

We report accurate estimation of cropping intensity (CI), an indicator of food production, is well aligned with the ongoing efforts to achieve sustainable development goals (SDGs) under diminishing natural resources. The advancement in satellite remote sensing provides unprecedented opportunities for capturing CI information in a spatially continuous manner. However, challenges remain due to the lack of generalizable algorithms for accurately and efficiently mapping global CI with a fine spatial resolution. In this study, we developed a 30-m planetary-scale CI mapping framework with the reconstructed time series of Normalized Difference Vegetation Index (NDVI) from multiple satellite images. Using a binary crop phenophase profile indicating growing and non-growing periods, we estimated pixel-by-pixel CI by enumerating the total number of valid cropping cycles during the study years. Based on the Google Earth Engine cloud computing platform, we implemented the framework to estimate CI during 2016–2018 in eight geographic regions across continents that are representative of global cropping system diversity. Comparison with PhenoCam network data in four cropland sites suggests that the proposed framework is capable of capturing the seasonal dynamics of cropping practices. Spatially, overall accuracies based on validation samples range from 80.0% to 98.9% across different regions worldwide. Regarding the CI classes, single cropping systems are associated with more robust and less biased estimations than multiple cropping systems. Finally, our CI estimates reveal high agreement with two widely used land surface phenology products, including Vegetation Index and Phenology V004 (VIP4) and Moderate Resolution Imaging Spectroradiometer Land Cover Dynamics (MCD12Q2), meanwhile providing much more spatial details. Due to its robustness, the developed CI framework can be potentially generalized to produce global fine resolution CI products for food security and other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using machine learning to improve land use/cover characterization and projection for scenario-based global modeling

The characterization of the land surface in models is critical for robust estimation of the water cycle and associated extremes. Land use determines water demand, and land cover affects water availability. The interaction between water demand and availability is the key determinant of whether systems are resilient to variability in the overall water cycle. Furthermore, activities such as irrigation and deforestation or afforestation may influence regional to global precipitation patterns, depending on the extent of the activity. Thus this proposed work aims to address the following question: How does integrated land use/cover data and improved land use/cover projection, as informed by machine learning approaches, better resolve human-earth system resilience to water cycle variability and extremes?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗
DOE OSTI2020

Future Warming and Intensification of Precipitation Extremes: A “Double Whammy” Leading to Increasing Flood Risk in California

This study focuses on quantifying future anthropogenic changes in surface runoff associated with extreme precipitation in California's Sierra Nevada. The method involves driving a land surface model with output from a high resolution regional atmospheric simulation of the most extreme atmospheric rivers (ARs). AR events were selected from an ensemble of global climate model simulations of historical and late 21st century climate under the “high-emission” RCP8.5 scenario. Average precipitation during the future ARs increases by ~25% but a much lower proportion falls as snow. The resulting future runoff increase is dramatic—nearly 50%, reflecting both the precipitation increase and simultaneous conversion of snow to rain. The “double whammy” impact on runoff is largest in the 2,000–2,500 m elevation band, where the snowfall loss and precipitation increase are both especially large. Overall, this huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Diagnosing Near-Surface Model Errors with Candidate Physics Parameterization Schemes for the Multiphysics Rapid Refresh Forecast System (RRFS) Ensemble during Winter over the Northeastern United States and Southern Great Plains

Abstract During the winter of 2020/21 an ensemble of FV3-LAM forecasts was produced over the contiguous United States for the Winter Weather Experiment using five physics suites. These forecasts are evaluated with the goal of optimizing physics parameterizations within the future operational Rapid Refresh Forecast System (RRFS) in the Unified Forecast System (UFS) realm and for selecting suitable physics suites for a multiphysics RRFS ensemble. The five physics suites have different combinations of land surface models (LSMs), planetary boundary layer (PBL) parameterizations, and surface layer schemes, chosen from those used in current and possible future operational systems and likely to be supported in the operational UFS. Full-season evaluation reveals a persistent near-surface cold bias in the U.S. Northeast from one suite and a nighttime warm bias in the southern Great Plains in another suite, while other suites have smaller biases. A representative case is chosen to diagnose the cause for each of these biases using sensitivity simulations with different physics combinations or modified parameters and verified with additional mesonet observations. The cold bias in the Northeast is attributed to aspects of the Noah-MP LSM over snow cover, where Noah-MP simulates lower soil water content, and thus lower thermal conductivity than other LSMs, leading to less upward ground heat flux during nighttime and consequently lower surface temperature. The nighttime warm bias found in the southern Great Plains is attributed to overestimation of vertical mixing in the K -profile-based eddy-diffusivity mass-flux (K-EDMF) PBL scheme and insufficient land–atmospheric coupling from the GFS surface layer scheme over short vegetation. A few key parameters driving these systematic biases are identified.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Effective climate sensitivity distributions from a 1D model of global ocean and land temperature trends, 1970–2021

Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.

54 ENVIRONMENTAL SCIENCES↗