Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “land surface model”

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 1,171 records · Page 65

A method for coupling a parameterization of the planetary boundary layer with a hydrologic model

Deardorff's parameterization of the planetary boundary layer is adapted to drive a hydrologic model. The method converts the atmospheric conditions measured at the anemometer height at one site to the mean values in the planetary boundary layer; it then uses the planetary boundary layer parameterization and the hydrologic variables to calculate the fluxes of momentum, heat and moisture at the atmosphere-land interface for a different site. A simplified hydrologic model is used for a simulation study of soil moisture and ground temperature on three different land surface covers. The results indicate that this method can be used to drive a spatially distributed hydrologic model by using observed data available at a meteorological station located on or nearby the site.

Lin, J. D.↗

Sarasota Climate: Monitoring Heat and Assessing Heat Vulnerability to Identify Locations for Heat Mitigation Efforts in Sarasota, Florida

The coastal county of Sarasota, Florida, is located within the humid subtropical climate region and experiences an average of 250 days of sunshine every year. The county has large, urbanized communities which are vulnerable to urban heat island (UHI) effects. These rapidly growing communities contribute to the increasing surface temperatures and placing more residents at risk of heat-related illness. In partnership with Sarasota County Sustainability, the Sarasota Climate team utilized Earth observation data from NASA Landsat 8 and Landsat 9’s Thermal Infrared Sensors (TIRS), and the International Space Station’s Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) to model UHI effects within the county during the Summer for the last five years, 2019 to 2023. Data analysis with the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) and Urban Heat Exposure Assessment Tempe 1.0 (UHEAT 1.0) models and within ModelBuilder produced maps that identified the land surface temperature (LST) variance within the county, the regions that are most susceptible to extreme heat, and areas least capable of mitigating the effects of UHI. The results revealed that heat intensity varies significantly across Sarasota County with the highest temperatures in the more developed western part of the county. Additionally, the team identified that there are at least three vulnerable communities that exist in high-heat regions, including North Sarasota, Venice, and North Port. These regions have an overlap between socioeconomic sensitivity and environmental hazard that indicate a high priority in future heat mitigation efforts.

Remote sensing↗

Desert Peak Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to the Desert Peak Geothermal Field. It includes all input and output files used in the project. The files include data categories of raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs including Radar, SWIR, Thermal, Geophysics, Geology, and Wells. The files for the Desert Peak Geothermal Site are used with the Geothermal Exploration Artificial Intelligence to identify indicators of blind geothermal systems. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Desert Peak Geothermal Field.

15 GEOTHERMAL ENERGY↗

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗

Global Simulation of Atmospheric Mercury Concentrations and Deposition Fluxes

Results from a numerical model of the global emissions, transport, chemistry, and deposition of mercury (Hg) in the atmosphere are presented. Hg (in the form of Hg(O) and Hg(II)) is emitted into the atmosphere from natural and anthropogenic sources (estimated to be 4000 and 2100 Mg/ yr, respectively). It is distributed between gaseous, aqueous and particulate phases. Removal of Hg from the atmosphere occurs via dry deposition and wet deposition, which are calculated by the model to be 3300 and 2800 Mg/ yr, respectively. Deposition on land surfaces accounts for 47% of total global deposition. The simulated Hg ambient surface concentrations and deposition fluxes to the Earth's surface are consistent with available observations. Observed spatial and seasonal trends are reproduced by the model, although larger spatial variations are observed in Hg(O) surface concentrations than are predicted by the model. The calculated atmospheric residence time of Hg is -1.7 years. Chemical transformations between Hg(O) and HG(II) have a strong influence on Hg deposition patterns because HG(II) is removed faster than Hg(O). Oxidation of Hg(O) to HG(II) occurs primarily in the gas phase, whereas HG(II) reduction to Hg(O) occurs solely in the aqueous phase. Our model results indicated that in the absence of the aqueous reactions the atmospheric residence time of Hg is reduced to 1.2 from 1.7 years and the Hg surface concentration is -25% lower because of the absence of the HG(II) reduction pathway. This result suggests that aqueous chemistry is an essential component of the atmospheric cycling of Hg.

Shia, Run-Lie↗

Analysis of Atmospheric Mesoscale Models for Entry, Descent and Landing

Each Mars Exploration Rover (MER) is sensitive to the martian winds encountered near the surface during the Entry, Descent and Landing (EDL) process. These winds are strongly influenced by local (mesoscale) conditions. In the absence of suitable wind observations, wind fields predicted by martian mesoscale atmospheric models have been analyzed to guide landing site selection. Two different models were used, the MRAMS model and the Mars MM5 model. In order to encompass both models and render their results useful to the EDL engineering team, a series of statistical techniques were applied to the model results. These analyses cover the high priority landing sites during the expected landing times (1200 to 1500 local time). The number of sites studied is limited by the computational and analysis cost of the mesoscale models.

D M Kass↗

Atmosphere-biosphere exchange flux of carbon dioxide in a tallgrass prairie modeled with satellite spectral data

The estimation of the rate of net CO2 uptake of vegetated land surfaces is essential for studies of global carbon cycle. The present paper demonstrates the use of spectral reflectance data from satellite remote sensing to model net CO2 flux (NCF) of a tallgrass canopy at the Konza prairie, Kansas. A bidirectional reflectance canopy model was used to estimate seasonal changes in canopy leaf area index (LAI) from surface reflectances remotely sensed by SPOT 1 and Landsat 5 satellites. The radiation model was also coupled with leaf conductance-photosynthesis models to scale up stomatal conductance and NCF from individual leaves to canopy level according to radiation distribution inside the canopy. The satellite-data-driven model was able to closely simulate the seasonal change in LAI as well as the short-term variation of canopy LAI caused by the dry period during late July and early August in the area. Modeled canopy stomatal conductance (g(sub c)) and NCF agree with measurements within 0.16 cm/s and 0.28 mg m(exp -2)/s, respectively, during the growth season from late May to late August. In October both measured and modeled NCF turned to small negative values as canopy photosynthesis diminished and predicted LAI approached zero. In addition to data scatter, some of the differences between modeled and measured g(sub c) and NCF may be attributed to uncertainties in seasonal changes of plant physiological status that were not detected by satellite data; some of the differences were caused by inadequate description of the dependence of nighttime CO2 flux of soil respiration on near-surface turbulent mixing.

Gao, W.↗

Representing Soil Microbial Dynamics and Organo‐Mineral Interactions in the E3SM Land Model (ELM‐ReSOM)

Explicit representation of soil microbial processes and interactions with biotic and abiotic processes in Earth System Models (ESMs) remains limited, despite their importance in biogeochemical cycles. To address this gap, which hinders prediction of global biogeochemial cycling and responses to atmospheric conditions, we integrated a microbe- and mineral-surface-explicit model, the Reaction-network-based model of soil organic matter and Microbes (ReSOM), into the Energy Exascale ESM (E3SM) land model (ELM). Here, we describe ELM-ReSOM and show a case study at a conifer forest in California. ELM-ReSOM accurately simulated surface CO 2 fluxes and SOM stocks, demonstrating improved representations of microbial and mineral interactions compared to the default ELM. We examined ELM-ReSOM sensitivity to microbial traits, enzyme properties, and organo-mineral interactions. Microbial traits such as the maximum mortality rate, transporter-density scaling factor, and maximum monomer assimilation rate were strong controllers of heterotrophic respiration, while these microbial traits and enzyme-related properties collectively influenced SOM stocks. Mineral surfaces primarily affected SOM stocks by adsorbing enzymes, thereby limiting depolymerization. Synergies among processes led to stronger impacts of parameters when evaluated together versus separately (i.e., most parameters had greater indirect than direct effects). For example, due to interactions of microbial necromass with mineral surface adsorption, the indirect effect of the maximum microbial mortality rate was 33% larger than its direct effect on SOM stock. Thus, microbial and enzyme dynamics and their interactions with mineral surfaces play critical roles in SOM cycling. Tackling the challenges of microbe-explicit models will advance understanding and modeling of SOM dynamics.

Tao, Jing [Lawrence Berkeley National Laboratory (↗

Current and Future Plans of the NASA Data Assimilation Office (DAO)

The mission of the Data Assimilation Office (DAO) is to advance the state of the art of data assimilation and produce research-quality assimilated data sets which make optimal use of space-based observations. Development efforts over the last few years have focused on delivering a production data assimilation system in support of NASA's Terra platform. That system, called the Goddard Earth Observing System - version 2 or GEOS-2, represents a major upgrade to the baseline GEOS-1 system employed in NASA's first reanalysis effort. GEOS-2 includes a physical-space three dimensional variational analysis algorithm (the Physical-space Statistical Analysis System or PSAS) and numerous improvements to the general circulation model. The latter include a Soil-Vegetation-Atmosphere Transfer (SVAT) land surface scheme, a level 2.5 moist turbulence scheme and new Short Wave (SW) and Long Wave (LW) radiation code. The system also includes an off-line ozone assimilation system, and the capability to assimilate scatterometer surface winds, and TIROS Operational Vertical Sounder (TOVS) and Special Sensor Microwave Imager (SSM/I) moisture data. GEOS-2 is currently run at 1 degree horizontal resolution and 48 levels extending to O.Olmb. Experimental versions of GEOS-2 are run with a global stretched grid allowing enhanced (e.g. 1/4 deg) regional resolution. Other capabilities being developed include, the assimilation of Tropical Rainfall Measuring Mission (TRMM) precipitation and Global Positioning System (GPS) data, an off-line land surface assimilation system, and a retrospective analysis scheme. The DAO is also engaged in a number of collaborative efforts to help accelerate the development of the next generation data assimilation system. These include, a joint modeling effort between the DAO and NCAR/CGDD to develop a new Global Circulation Model (GCM), and a Department of Energy Lawrence Livermore National Laboratory (DOE/LLNL) collaboration on model parallelization. Plans for the next reanalysis will be discussed in the context of current and near term system quality and computing capabilities, and the need for multiple reanalysis products.

Atlas, Robert↗

Brady Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to Brady's Geothermal Field. It includes all input and output files for the Geothermal Exploration Artificial Intelligence. Input and output files are sorted into three categories: raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs which are titled Radar, SWIR, Thermal, Geophysics, Geology, and Wells. These inputs and outputs were used with the Geothermal Exploration Artificial Intelligence to identify indicators of blind geothermal systems at the Brady Hot Springs Geothermal Site. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Brady Hot Springs Geothermal Site.

15 GEOTHERMAL ENERGY↗

Salton Sea Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to Salton Sea Geothermal Field. It includes all input and output files used with the Geothermal Exploration Artificial Intelligence. Input and output files are sorted into three categories: raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs which are titled Radar, SWIR, Thermal, Geophysics, Geology, and Wells. The files are used with the Geothermal Exploration Artificial Intelligence for the Salton Sea Geothermal Site to identify indicators of blind geothermal systems. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Salton Sea Geothermal Site.

15 GEOTHERMAL ENERGY↗

A Realization of Bias Correction Method in the GMAO Coupled System

Over the past several decades, a tremendous effort has been made to improve model performance in the simulation of the climate system. The cold or warm sea surface temperature (SST) bias in the tropics is still a problem common to most coupled ocean atmosphere general circulation models (CGCMs). The precipitation biases in CGCMs are also accompanied by SST and surface wind biases. The deficiencies and biases over the equatorial oceans through their influence on the Walker circulation likely contribute the precipitation biases over land surfaces. In this study, we introduce an approach in the CGCM modeling to correct model biases. This approach utilizes the history of the model's short-term forecasting errors and their seasonal dependence to modify model's tendency term and to minimize its climate drift. The study shows that such an approach removes most of model climate biases. A number of other aspects of the model simulation (e.g. extratropical transient activities) are also improved considerably due to the imposed pre-processed initial 3-hour model drift corrections. Because many regional biases in the GEOS-5 CGCM are common amongst other current models, our approaches and findings are applicable to these other models as well.

Chang, Yehui↗

A Comparison of Methods for a Priori Bias Correction in Soil Moisture Data Assimilation

Data assimilation is being increasingly used to merge remotely sensed land surface variables such as soil moisture, snow and skin temperature with estimates from land models. Its success, however, depends on unbiased model predictions and unbiased observations. Here, a suite of continental-scale, synthetic soil moisture assimilation experiments is used to compare two approaches that address typical biases in soil moisture prior to data assimilation: (i) parameter estimation to calibrate the land model to the climatology of the soil moisture observations, and (ii) scaling of the observations to the model s soil moisture climatology. To enable this research, an optimization infrastructure was added to the NASA Land Information System (LIS) that includes gradient-based optimization methods and global, heuristic search algorithms. The land model calibration eliminates the bias but does not necessarily result in more realistic model parameters. Nevertheless, the experiments confirm that model calibration yields assimilation estimates of surface and root zone soil moisture that are as skillful as those obtained through scaling of the observations to the model s climatology. Analysis of innovation diagnostics underlines the importance of addressing bias in soil moisture assimilation and confirms that both approaches adequately address the issue.

Kumar, Sujay V.↗

The Application of Satellite-Derived, High-Resolution Land Use/Land Cover Data to Improve Urban Air Quality Model Forecasts

Local and state agencies are responsible for developing state implementation plans to meet National Ambient Air Quality Standards. Numerical models used for this purpose simulate the transport and transformation of criteria pollutants and their precursors. The specification of land use/land cover (LULC) plays an important role in controlling modeled surface meteorology and emissions. NASA researchers have worked with partners and Atlanta stakeholders to incorporate an improved high-resolution LULC dataset for the Atlanta area within their modeling system and to assess meteorological and air quality impacts of Urban Heat Island (UHI) mitigation strategies. The new LULC dataset provides a more accurate representation of land use, has the potential to improve model accuracy, and facilitates prediction of LULC changes. Use of the new LULC dataset for two summertime episodes improved meteorological forecasts, with an existing daytime cold bias of approx. equal to 3 C reduced by 30%. Model performance for ozone prediction did not show improvement. In addition, LULC changes due to Atlanta area urbanization were predicted through 2030, for which model simulations predict higher urban air temperatures. The incorporation of UHI mitigation strategies partially offset this warming trend. The data and modeling methods used are generally applicable to other U.S. cities.

Quattrochi, D. A.↗

GLACE: The Global Land-Atmosphere Coupling Experiment: Overview - Part 1

GLACE is a model intercomparison study focusing on a typically neglected yet critical element of numerical weather and climate modeling: land-atmosphere coupling strength, or the degree to which anomalies in land surface state (e.g., soil moisture) can affect rainfall generation and other atmospheric processes. The twelve AGCM groups participating in GLACE performed a series of simple numerical experiments that allow the objective quantification of this element. The derived coupling strengths vary widely. Some similarity, however, is found in the spatial patterns generated by the models, enough similarity to pinpoint multi-model "hot spots" of land-atmosphere coupling. For boreal summer, such hot spots for precipitation and temperature are found over large regions of Africa, central North America and India; a hot spot for temperature is also found over eastern China. The design of the GLACE simulations are described in full detail so that any interested modeling group can repeat them easily and thereby place their model s coupling strength within the broad range of those documented here.

Koster, Randal D.↗

Drone flight visible images, digital elevation maps, and geotiffs, Council, Seward Peninsula, Alaska, USA, July 2017

Remote sensing data collected from Lawrence Berkeley National Laboratory’s (LBNL) unmanned aerial system (UAS) quadcopter platform – Inspire-1 – operated by the Ameriflux Management Project (ameriflux.lbl.gov) Technical Team. The Insprire-1 is a visible imagery UAS platform that simultaneously measures very high spatial resolution optical red/green/blue (RGB) at regular intervals along each flight path. This package includes data from four flight paths over the NGEE Arctic Council Road Site near mile marker 71 in the Seward Peninsula, Alaska in July 2017. Derived image products for each flight include ortho-mosaiced RGB, an RGB-based digital surface model (DSM) using the structure from motion (SfM) technique, and a digital terrain model (DTM). Ancillary aircraft data, flight mission parameters, and general flight conditions are also included. Unprocessed and processed data products are included in this package (processing levels 0-2). This metadata document contains flight campaign, instrument and file metadata, along with a description of data processing levels, data products and file naming scheme.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Drone flight visible images, digital elevation maps, and geotiffs, Teller, Seward Peninsula, Alaska, USA, July 2017

Remote sensing data collected from Lawrence Berkeley National Laboratory’s (LBNL) unmanned aerial system (UAS) quadcopter platform – Inspire-1 – operated by the Ameriflux Management Project (ameriflux.lbl.gov) Technical Team. The Insprire-1 is a visible imagery UAS platform that simultaneously measures very high spatial resolution optical red/green/blue (RGB) at regular intervals along the flight path. This package includes data from a flight path over the NGEE Arctic Teller Site, west of mile marker 27 Bob Blodgett Highway in the Seward Peninsula, Alaska in July 2017. Derived image products include ortho-mosaiced RGB, an RGB-based digital surface model (DSM) using the structure from motion (SfM) technique, and a digital terrain model (DTM). Ancillary aircraft data, flight mission parameters, and general flight conditions are also included. Unprocessed and processed data products are included in this package (processing levels 0-2). This metadata document contains flight campaign, instrument and file metadata, along with a description of data processing levels, data products and file naming scheme.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Drone flight visible images, digital elevation maps, and geotiffs, Kougarok, Seward Peninsula, Alaska, USA, July 2017

Remote sensing data collected from Lawrence Berkeley National Laboratory’s (LBNL) unmanned aerial system (UAS) quadcopter platform – Inspire-1 – operated by the Ameriflux Management Project (ameriflux.lbl.gov) Technical Team. The Insprire-1 is a visible imagery UAS platform that simultaneously measures very high spatial resolution optical red/green/blue (RGB) at regular intervals along each flight path. This package includes data from four flight paths over the NGEE Arctic Kougarok Site, west of mile marker 64 on Kuzitrin Road in the Seward Peninsula, Alaska in July 2017. Derived image products include ortho-mosaiced RGB, an RGB-based digital surface model (DSM) using the structure from motion (SfM) technique, and a digital terrain model (DTM). Ancillary aircraft data, flight mission parameters, and general flight conditions are also included. Unprocessed and processed data products are included in this package (processing levels 0-2). This metadata document contains flight campaign, instrument and file metadata, along with a description of data processing levels, data products and file naming scheme.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗