Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Soil aggregation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Remote spectral identification of surface aggregates by thermal imaging techniques - Progress report

The NASA Thermal Infrared Multispectral Scanner (TIMS) has been successfully used for the remote identification of a variety of soil and aggregate deposits in vegetated areas of two states. Over three million cubic meters of gravel deposits were identified from the imagery during a two year period. Verification was accomplished by ground reconnaissance using drilling machinery and by ground instrumentation. The method has been used to differentiate between fine and coarse grained soils, and gravel deposits. The deposits were found to have been naturally sorted according to grain size by depositional processes, providing each deposit with distinct spectral qualities. It was found that the masking effects of relatively dense vegetation were largely overcome by using imagery acquired at higher altitudes above terrain than 9000 meters, due to loss of resolution of the finer detail. The mechanics of image resolution are discussed, a method of data analysis used is described, and sample spectral signatures are illustrated.

Scholen, Douglas E.↗
NASA NTRS1994

Aspects of spatial and temporal aggregation in estimating regional carbon dioxide fluxes from temperate forest soils

We examine the influence of aggregation errors on developing estimates of regional soil-CO2 flux from temperate forests. We find daily soil-CO2 fluxes to be more sensitive to changes in soil temperatures (Q(sub 10) = 3.08) than air temperatures (Q(sub 10) = 1.99). The direct use of mean monthly air temperatures with a daily flux model underestimates regional fluxes by approximately 4%. Temporal aggregation error varies with spatial resolution. Overall, our calibrated modeling approach reduces spatial aggregation error by 9.3% and temporal aggregation error by 15.5%. After minimizing spatial and temporal aggregation errors, mature temperate forest soils are estimated to contribute 12.9 Pg C/yr to the atmosphere as carbon dioxide. Georeferenced model estimates agree well with annual soil-CO2 fluxes measured during chamber studies in mature temperate forest stands around the globe.

Kicklighter, David W.↗
DOE OSTI2025

WHONDRS laboratory time series moisture manipulative experiment from soil core layers across eastern contiguous US: time series aerobic respiration, geochemistry, and aggregates

This dataset supports a broader study examining the effects of wetting and drying on soil layers across the eastern contiguous United States (CONUS). The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata. Samples were collected as part of a collaboration between WHONDRS (Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems; https://whondrs.pnnl.gov) and MONet (Molecular Observation Network; https://www.emsl.pnnl.gov/monet). The field samples (soil cores) were labeled as MEL_##_COR and subsequent subsamples begin with MEL_##. Additional subsamples were taken for the laboratory experiment and were labeled as EL_##. The labels from the MEL field samples and the EL subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EL_01 is a subsample from MEL_01). See the critical details section below for more details on sample naming and experimental design.For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions.This dataset is comprised of (1) a folder containing environmental context photos; (2) file-level metadata; (3) data dictionary; (4) field metadata; (5) readme; (6) international generic sample number (IGSN) mapping file; and (7) a subfolder with soil sample data from field samples and the incubation experiment. The sample data subfolder contains (1) effect size; (2) gravimetric moisture from field samples and incubation experiment; (3) respiration rates, raw dissolved oxygen values, and plots; (4) specific conductance, pH, and temperature from the incubation; (5) soil aggregates; (6) a summary containing median values of each data type for each treatment (wet and dry) in the incubation; (7) a summary containing averages for each data type of each soil layer; and (8) methods codes. All files are .csv, .pdf, .jpeg, or .jpg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Land management effects on wet aggregate stability and carbon content

Land management affects soil structure and many other soil properties and processes. Our objectives were to evaluate soil organic C (SOC), aggregate size distribution, aggregate-associated C, and soil structure as affected by long-term land management and slope. A chronosequence of 38 on-farm sites with low to high (5–18%) slopes was selected to evaluate 5–40 yr of management. The sites were classified as business as usual (BAU) cropland (BAU-Crop), BAU pasture (BAU-Past), newly established conservation reserve program (CRP) areas (CRP-New), and established CRP (CRP-Old). Soil samples were collected from the 0-to-5- and 5-to-15-cm depth increments and processed for soil property measurements including fractionation by wet sieving into five aggregate size classes (>2,000, 1,000–2,000, 500–1,000, 250–500, and 53–250 μm). Within the surface 5 cm, mean weight diameter (MWD) and geometric mean diameter (GMD) were used to characterize soil structural stability. The BAU-Past and CRP-Old sites had 79% more macroaggregates (>2,000, 1,000–2,000, and 500–1,000 μm), 123% higher MWD, 38% higher GMD, and 47% higher SOC than BAU-Crop or CRP-New sites. The 5-to-15-cm depth increment showed a similar but lower magnitude response. Aggregate-associated C was quantified using a constant soil mass that reflected aggregate size distribution to prevent overestimating C content. Lower-slope locations had more SOC, more macroaggregates, more C associated with macroaggregates, and higher GMD and MWD compared with high-slope locations across all management classifications and soil depths. The results support our hypothesis that the high-slop soils may benefits from specific management decisions than the lower-sloping soils as a function of landscape property. We recommend reestablishing grassland on sloping land that is susceptible to excessive soil erosion, although those practices will likely take a long time to restore soil structural stability and SOC content to precultivation levels.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

Spatial and Temporal Dust Source Variability in Northern China Identified Using Advanced Remote Sensing Analysis

The aim of this research is to provide a detailed characterization of spatial patterns and temporal trends in the regional and local dust source areas within the desert of the Alashan Prefecture (Inner Mongolia, China). This problem was approached through multi-scale remote sensing analysis of vegetation changes. The primary requirements for this regional analysis are high spatial and spectral resolution data, accurate spectral calibration and good temporal resolution with a suitable temporal baseline. Landsat analysis and field validation along with the low spatial resolution classifications from MODIS and AVHRR are combined to provide a reliable characterization of the different potential dust-producing sources. The representation of intra-annual and inter-annual Normalized Difference Vegetation Index (NDVI) trend to assess land cover discrimination for mapping potential dust source using MODIS and AVHRR at larger scale is enhanced by Landsat Spectral Mixing Analysis (SMA). The combined methodology is to determine the extent to which Landsat can distinguish important soils types in order to better understand how soil reflectance behaves at seasonal and inter-annual timescales. As a final result mapping soil surface properties using SMA is representative of responses of different land and soil cover previously identified by NDVI trend. The results could be used in dust emission models even if they are not reflecting aggregate formation, soil stability or particle coatings showing to be critical for accurately represent dust source over different regional and local emitting areas.

spectral mixing analysis↗
DOE OSTI2026

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗
NASA NTRS1973

The distribution in lunar soil of carbon released by pyrolysis

The carbon contents of various lunar soil particle types and sieve fractions of Apollo 15 and 16 samples have been determined by the pyrolysis method. The mineral, glass, and high-grade breccia fragments in the soils examined contain relatively low amounts of carbon (approximately 8, 25, and 25 microg C/g sample respectively in 149-250 micron grains). Most low-grade breccias and all agglutinates examined have high carbon contents (approximately 52 and 80 microg C/g sample respectively), and agglutinate abundance is indicative of the carbon content and maturity of a soil. The distribution of carbon with respect to particle size in mature soils generally reveals a minimum in carbon content at about 100 micron particle diameter. At smaller particle diameters, carbon content is directly proportional to particle surface area and therefore increases with the ratio (surface area)/(particle mass). A model relating the cycle of comminution and aggregation of soil particles to the redistribution of surface implanted carbon is developed.

Desmarais, D. J.↗
NASA NTRS1971

Lunar-Surface Closeup Stereoscopic Photography on the Sea of Tranquility (Apollo 11 Landing Site)

Analysis of returned lunar samples provides limited information about lunar geology. To obtain information about in-place lunar material, a closeup stereoscopic camera capable of photographing small-scale surface features was built and was used at the Apollo 11 landing site. Stereoscopic photographs were taken of surface areas relative to the lunar module, and the surfaces photographed were analyzed. The photographs are classified into five groups: soil disturbed by astronaut activities, generally undisturbed soil, loose aggregate surface material, crater bottoms with prominent glass deposits, and hard rock deposits. Glass deposits in the returned samples are described for comparison with the features observed in the photographs. The stereoscopic photographs were of outstanding quality and show the nature of lunar-surface material in detail. Lunar topography was reconstructed from the photographs with an analytical plotter. The photography results indicate that the closeup composition and genesis of lunar soil at the Apollo 11 landing site.

Greenwood, W. R.↗
DOE OSTI2025

A framework for testing soil carbon dynamics post land-use transition in a multisector dynamics model

Soil carbon plays a crucial role in the global carbon cycle. Changes in land use can determine whether carbon is stored or is emitted into the atmosphere as carbon dioxide, which has broad implications for the human and Earth systems. These feedbacks to the carbon cycle and their socio-economic drivers are modelled by many global multisector dynamics models to project future possibilities for the human-Earth system. One notable model of this class is the Global Change Analysis Model (GCAM), which uses a simplified process to model soil organic carbon (SOC) content after land-use transition across 384 land units. While the current GCAM soil carbon framework is based on scientific principles, it has not been tested against experimental data. This work examines rates of SOC change from GCAM input data. Specifically, first order rate constants derived from model inputs were compared to values from two syntheses to assess GCAM’s accuracy. Welch’s t-tests and linear models were used to determine if rate constants were consistent across all tested geographical areas and land-use transition types. While we found that there was general agreement on the direction and magnitude (i.e., rate) of SOC change, the rate constant derived from GCAM and empirical values differed strongly in a subset of specific instances. These results indicate that GCAM’s current SOC dynamics during land use transition successfully capture broad patterns of change in this critical carbon pool, but should be interpreted with caution at finer spatial scales. One potential cause of these discrepancies is our highly aggregated variable, soil timescale, which could be made more granular to improve accuracy. When using economically rooted multisector dynamics models, such as GCAM, it is critical to understand such model limitations for representing specific Earth system processes.

carbon↗
NASA NTRS1975

Chemistry of agglutinate fractions in lunar soils

Agglutinates are aggregates of crystalline grains and lithic fragments bonded together by glass. It is thought that glassy agglutinates are formed at the upper surface of the lunar regolith by impact-related melting and welding of soil particles, in response to meteoroid and micrometeoroid bombardment. A description is presented of an investigation in which bulk soils were separated into 'agglutinate' and 'nonagglutinate' fractions. The obtained fractions were analyzed for major, minor, and trace elements. The obtained chemical data for agglutinate and nonagglutinate fractions of lunar soils indicate that agglutinitic glass is enriched in mafic and most lithophile elements relative to the bulk soils. A model involving preferential melting and assimilation of mesostasis material and mafic soil components is proposed to account for the observed chemical data. It is suggested that glassy agglutinates may form more readily in mafic soils than in more feldspathic ones. Such selectivity should be most effective between mare and highland soils, but may possibly operate on a more local scale.

Rhodes, J. M.↗
NASA NTRS1993

Spatial and temporal variability of soil temperature, moisture and surface soil properties

The overall objectives of this research were to: (l) Relate in-situ measured soil-water content and temperature profiles to remotely sensed surface soil-water and temperature conditions; to model simultaneous heat and water movement for spatially and temporally changing soil conditions; (2) Determine the spatial and temporal variability of surface soil properties affecting emissivity, reflectance, and material and energy flux across the soil surface. This will include physical, chemical, and mineralogical characteristics of primary soil components and aggregate systems; and (3) Develop surface soil classes of naturally occurring and distributed soil property assemblages and group classes to be tested with respect to water content, emissivity and reflectivity. This document is a report of studies conducted during the period funded by NASA grants. The project was designed to be conducted over a five year period. Since funding was discontinued after three years, some of the research started was not completed. Additional publications are planned whenever funding can be obtained to finalize data analysis for both the arid and humid locations.

Hajek, B. F.↗
NASA NTRS1998

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Aggregation - 2

This is a follow up on the preceding presentation by Crosson. The grid size for remote microwave measurements is much coarser than the hydrological model computational grids. To validate the hydrological models with measurements we propose mechanisms to aggregate the hydrological model outputs for soil moisture to allow comparison with measurements. Weighted neighborhood averaging methods are proposed to facilitate the comparison. We will also discuss such complications as misalignment, rotation and other distortions introduced by a generalized sensor image.

Schamschula, Marius↗
NASA NTRS1975

Dry heat effects on survival of indigenous soil particle microflora and particle viability studies of Kennedy Space Center soil

Research efforts were concentrated on attempts to obtain data concerning the dry heat resistance of particle microflora in Kennedy Space Center soil samples. The in situ dry heat resistance profiles at selected temperatures for the aggregate microflora on soil particles of certain size ranges were determined. Viability profiles of older soil samples were compared with more recently stored soil samples. The effect of increased particle numbers on viability profiles after dry heat treatment was investigated. These soil particle viability data for various temperatures and times provide information on the soil microflora response to heat treatment and are useful in making selections for spacecraft sterilization cycles.

Ruschmeyer, O. R.↗
NASA NTRS2015

How the Assumed Size Distribution of Dust Minerals Affects the Predicted Ice Forming Nuclei

The formation of ice in clouds depends on the availability of ice forming nuclei (IFN). Dust aerosol particles are considered the most important source of IFN at a global scale. Recent laboratory studies have demonstrated that the mineral feldspar provides the most efficient dust IFN for immersion freezing and together with kaolinite for deposition ice nucleation, and that the phyllosilicates illite and montmorillonite (a member of the smectite group) are of secondary importance.A few studies have applied global models that simulate mineral specific dust to predict the number and geographical distribution of IFN. These studies have been based on the simple assumption that the mineral composition of soil as provided in data sets from the literature translates directly into the mineral composition of the dust aerosols. However, these tables are based on measurements of wet-sieved soil where dust aggregates are destroyed to a large degree. In consequence, the size distribution of dust is shifted to smaller sizes, and phyllosilicates like illite, kaolinite, and smectite are only found in the size range 2 m. In contrast, in measurements of the mineral composition of dust aerosols, the largest mass fraction of these phyllosilicates is found in the size range 2 m as part of dust aggregates. Conversely, the mass fraction of feldspar is smaller in this size range, varying with the geographical location. This may have a significant effect on the predicted IFN number and its geographical distribution.An improved mineral specific dust aerosol module has been recently implemented in the NASA GISS Earth System ModelE2. The dust module takes into consideration the disaggregated state of wet-sieved soil, on which the tables of soil mineral fractions are based. To simulate the atmospheric cycle of the minerals, the mass size distribution of each mineral in aggregates that are emitted from undispersed parent soil is reconstructed. In the current study, we test the null-hypothesis that simulating the presence of a large mass fraction of phyllosilicates in dust aerosols in the size range 2 m, in comparison to a simple model assumption where this is neglected, does not yield a significant effect on the magnitude and geographical distribution of the predicted IFN number. Results from sensitivity experiments are presented as well.

Dust↗
DOE OSTI2023

Evaluation of Multi-Fidelity Soil Moisture Products Across the Continental United States

We have aggregated the most recent soil moisture datasets from a diverse range of sources, encompassing the Continental United States (CONUS). These sources encompass gridded data from remote sensing products, reanalysis products, machine learning-based projects, and land surface modeling products. Additionally, we have obtained and processed in-situ soil moisture observations from the International Soil Moisture Network. The collected datasets exhibit variations in both temporal and spatial resolutions. Among the 20 datasets, six are available at a spatial resolution of 0.25 degrees, while three are at a coarser spatial resolution of 25 km. To minimize spatial interpolation, we conducted data uncertainty evaluations at the 0.25-degree spatial resolution. For our data evaluations, we maintained a monthly temporal resolution, which effectively captures soil moisture seasonality and interannual variability. Our data processing strategy preserves the raw data and interpolated data at their original temporal resolutions. Datasets with higher temporal resolutions, including daily, three-hourly, and hourly datasets, are set aside for subsequent analyses. These analyses will delve into topics such as soil moisture changes and recovery during extreme weather events. Furthermore, we have processed auxiliary data to enhance our evaluation, leveraging tools such as Google Earth Engine. This includes incorporating topography data, land use land cover data, Köppen-Geiger climate classification, and more to provide a comprehensive assessment from multiple sources.

Li, Lingcheng↗
DOE OSTI2021

Microbiome Aggregated Traits and Assembly Are More Sensitive to Soil Management than Diversity

How soil is managed, particularly for agriculture, exerts stresses upon soil microbiomes, resulting in altered community structures and functional states. Understanding how soil microbiomes respond to combined stresses is important for predicting system performance under different land use scenarios, aids in identification of the most environmentally benign managements, and provides insight into how system function can be recovered in degraded soils. We use a long-established field experiment to study the effects of combined chronic (press) disturbance of the magnitude of organic carbon inputs with acute (pulse) effects of physical disturbance by tillage and chemical disturbance due to inorganic fertilization and pesticide application. We show that because of the variety of ways it can be assessed, biodiversity—here based on microbial small subunit rRNA gene phylotypes—does not provide a consistent view of community change. In contrast, aggregated traits associated with soil microbiomes indicate general loss of function, measured as a reduction of average genome lengths, associated with chronic reduction of organic inputs in arable or bare fallow soils and altered growth strategies associated with rRNA operon copy number in prokaryotes, as well as a switch to pathogenicity in fungal communities. In addition, pulse disturbance by soil tillage is associated with an increased influence of stochastic processes upon prokaryote community assembly, but fungicide used in arable soils results in niche assembly of fungal communities compared to untilled grassland. Overall, bacteria, archaea, and fungi do not share a common response to land management change, and estimates of biodiversity do not capture important facets of community adaptation to stresses adequately.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Effects of soil particles and convective transport on dispersion and aggregation of nanoplastics via small-angle neutron scattering (SANS) and ultra SANS (USANS)

Terrestrial nanoplastics (NPs) pose a serious threat to agricultural food production systems due to the potential harm of soil-born micro- and macroorganisms that promote soil fertility and ability of NPs to adsorb onto and penetrate into vegetables and other crops. Very little is known about the dispersion, fate and transport of NPs in soils. This is because of the challenges of analyzing terrestrial NPs by conventional microscopic techniques due to the low concentrations of NPs and absence of optical transparency in these systems. Herein, we investigate the potential utility of small-angle neutron scattering (SANS) and Ultra SANS (USANS) to probe the agglomeration behavior of NPs prepared from polybutyrate adipate terephthalate, a prominent biodegradable plastic used in agricultural mulching, in the presence of vermiculite, an artificial soil. SANS with the contrast matching technique was used to study the aggregation of NPs co-dispersed with vermiculite in aqueous media. We determined the contrast match point for vermiculite was 66% D 2 O / 33% H 2 O. At this condition, the signal for vermiculite was ~50–100%-fold lower that obtained using neat H 2 O or D 2 O as solvent. According to SANS and USANS, smaller-sized NPs (50 nm) remained dispersed in water and did not undergo size reduction or self-agglomeration, nor formed agglomerates with vermiculite. Larger-sized NPs (300–1000 nm) formed self-agglomerates and agglomerates with vermiculite, demonstrating their significant adhesion with soil. However, employment of convective transport (simulated by ex situ stirring of the slurries prior to SANS and USANS analyses) reduced the self-agglomeration, demonstrating weak NP-NP interactions. Convective transport also led to size reduction of the larger-sized NPs. Therefore, this study demonstrates the potential utility of SANS and USANS with contrast matching technique for investigating behavior of terrestrial NPs in complex soil systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Aligning theoretical and empirical representations of soil carbon-to-nitrogen stoichiometry with process-based terrestrial biogeochemistry models

Soil carbon-nitrogen (C:N) stoichiometry acts as a control over decomposition and soil organic matter formation and loss, making it a key soil property for understanding ecosystem dynamics and projected ecosystems responses to global environmental change. However, the controls of soil C:N and how they respond to increasing pressures from global change agents are not fully understood. The “foundational” controls on soil C:N, namely plant and microbial C:N, have been used to predict soil C:N, but fail to accurately simulate all ecosystems and may be insufficient for predictions under global environmental change. Here we present an “emerging” representation of controls of soil C:N that includes plant-microbe-mineral feedbacks that have been shown to regulate soil C:N. We argue that including representation of these emerging drivers in process-based terrestrial biogeochemistry models, which include biological N fixation, mycorrhizae, priming, root exudation of organic acids, and mineralogy (including soil texture, mineral composition, and aggregation), will improve mechanistic representation of soil C:N and associated processes. Such improvements will produce models that will better simulate a variety of ecological states and predict soil C:N when global changes modify plant-microbe-mineral interactions. Here, we align our empirical understanding of controls of soil C:N with those controls represented in models, identifying contexts where emerging drivers might be particularly important to represent (e.g., priming and root exudation in nutrient-limited conditions) and areas of future work. Additionally, we show that implementing emerging drivers of soil C:N results in different simulated outcomes at steady state and in response to elevated atmospheric CO2. Our review and preliminary simulations support the need to incorporate emerging drivers of soil C:N into process-based terrestrial biogeochemistry models, allowing for both theoretical exploration of mechanisms and potentially more accurate predictions of land biogeochemical responses to global change.

54 ENVIRONMENTAL SCIENCES↗