Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data dependencies”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Accelerated Deterministic Phonon Transport With Consistent Material Temperature and Intensities

Abstract We present a method for deterministically solving the frequency and temperature dependent phonon radiative transport (PRT) equation in the single-mode relaxation time (SMRT) approximation in the self-adjoint angular flux (SAAF) form. To handle the nonlinear coupling between the phonon intensities and the material temperature, we apply a linearization approach that is similar to one in thermal radiative transport. This procedure leads to the PRT equation with pseudo-scattering. The method presented includes acceleration of both the inner pseudo-scattering source iterations and outer temperature iteration with a gray diffusion synthetic acceleration (DSA) and Anderson acceleration, respectively. We use the finite-element method to discretize the PRT equation in space and the method of discrete ordinates (SN) for angular discretization. The proposed method is verified by a gray method of manufactured solutions problem and demonstrated on a problem using temperature and direction dependent multigroup data from lithium aluminate (LiAlO2). The iterative performance of the acceleration method in each test is then compared to the unaccelerated method.

Engineering↗
DOE OSTI2025

Detection of local-scale changes in greenhouse gas emissions in urban environments using micrometeorological methods and comparison to a high-resolution inventory

We used the Monin–Obukhov similarity theory (MOST) flux-variance relationship to estimate greenhouse gas (GHG) fluxes from high-precision mole fraction measurements at 3 instrumented urban communication towers over several years, demonstrating the ability of this method to detect and quantify changes in emissions. Depending on data availability, we used carbon dioxide (CO 2 ) and carbon monoxide (CO) measurements and/or tracer ratios to estimate fluxes at 1 urban site (Site 3) and 1 suburban site (Site 7) in Indianapolis, IN, USA, and 1 urban site (COM) in Los Angeles, CA, USA. We also compared the estimated fluxes of CO 2 from fossil fuel sources (CO 2 ff) at Sites 3 and 7 and the total CO 2 fluxes at Site 3 to 20 m, hourly resolution subdomains of the high-resolution CO 2 emissions inventory, Hestia, for the year 2020, introducing a new way to evaluate emissions inventories at small spatial and temporal scales. Using the flux-variance relationship, we detected and quantified abrupt decreases in CO and CO 2 fluxes at Site 3 and COM in April 2020, coinciding with the stay-at-home order due to COVID-19 pandemic, as well as abrupt decreases in CO and CO 2 fluxes at Site 3 in July 2018 coinciding with a highway closure next to the site. The Hestia emissions inventory detected a decrease in emissions in April 2020 at Sites 3 and 7, but this decrease differed in magnitude from those detected in the atmospheric estimates. Seasonal trends in emissions are similar between Hestia and the atmospheric estimates at Site 7. We use differences in seasonal and spatial trends between the flux estimation methods to identify potential sources of uncertainty in both the atmospheric and inventory methods. The results from this study show that the flux-variance estimation method is a useful tool to monitor local-scale emissions and evaluate high-resolution emissions inventories.

Carbon dioxide↗
DOE OSTI2024

Updated Filter Leak Frequencies for Use in Risk Assessments

Quantitative risk assessment (QRA) is highly dependent on data, leading to more robust models as new and updated data is acquired. The Hydrogen Plus Other Alternative Fuels Risk Assessment (HyRAM+) QRA capabilities include calculations of individual risk from leaks in a gaseous hydrogen facility due to the potential effects of jet fires and explosions. Leak frequencies are acquired through statistical analysis of published data from a variety of sources and industries. The filter leak frequencies in previous versions of the HyRAM+ software are substantially greater than the leak frequencies of other components, leading to QRA results for gaseous hydrogen in which filters consistently dominate the overall risk. Data that were previously used to derive the filter leak frequencies were reevaluated for applicability and additional data points were added to update the filter leak frequencies. The new frequencies are more comparable to leak frequencies for other components.

08 HYDROGEN↗
DOE OSTI2025

MassIVE MSV000095714

Proteomic analysis reveals translation-related proteins are significantly altered in response to stabilized G4s. Identification of the importance of translation factors in overcoming G4s led to the question of how E. coli cells generally respond to chemicals that stabilize such structures. As a first step in addressing this question, a proteomic analysis was carried out to measure the quantitative effects of NMM on the levels of individual proteins in E. coli. Protein levels from early log-phase cultures of delta-tolC and delta-tolC tufA::kan strains grown in the presence or absence of NMM were measured to assess how reduced EF- Tu levels and NMM impacted expression.

Bottom-up proteomics↗
DOE OSTI2022

Towards Enhancing Coding Productivity for GPU Programming Using Static Graphs

The main contribution of this work is to increase the coding productivity of GPU programming by using the concept of Static Graphs. GPU capabilities have been increasing significantly in terms of performance and memory capacity. However, there are still some problems in terms of scalability and limitations to the amount of work that a GPU can perform at a time. To minimize the overhead associated with the launch of GPU kernels, as well as to maximize the use of GPU capacity, we have combined the new CUDA Graph API with the CUDA programming model (including CUDA math libraries) and the OpenACC programming model. We use as test cases two different, well-known and widely used problems in HPC and AI: the Conjugate Gradient method and the Particle Swarm Optimization. In the first test case (Conjugate Gradient) we focus on the integration of Static Graphs with CUDA. In this case, we are able to significantly outperform the NVIDIA reference code, reaching an acceleration of up to 11x thanks to a better implementation, which can benefit from the new CUDA Graph capabilities. In the second test case (Particle Swarm Optimization), we complement the OpenACC functionality with the use of CUDA Graph, achieving again accelerations of up to one order of magnitude, with average speedups ranging from 2x to 4x, and performance very close to a reference and optimized CUDA code. Our main target is to achieve a higher coding productivity model for GPU programming by using Static Graphs, which provides, in a very transparent way, a better exploitation of the GPU capacity. The combination of using Static Graphs with two of the current most important GPU programming models (CUDA and OpenACC) is able to reduce considerably the execution time w.r.t. the use of CUDA and OpenACC only, achieving accelerations of up to more than one order of magnitude. Finally, we propose an interface to incorporate the concept of Static Graphs into the OpenACC Specifications.

58 GEOSCIENCES↗
DOE OSTI2026

Added value of site load measurements in probabilistic lifetime extension: a Lillgrund case study

Site-specific fatigue estimation is an essential part of wind turbine lifetime extension, with various methods depending on data availability. The present study compares probabilistic lifetime extension assessment results for rotor blades with and without load measurements. It also addresses two key questions in such assessments: the applicability of the Frandsen model for estimating waked turbulence under complex and mixed wake conditions and the extrapolation of mid-term data over longer time periods. The case study wind turbine is SWT-2.3-93, located at the edge of the Lillgrund wind farm, situated in the Øresund Strait between Denmark and Sweden. The turbine is extensively instrumented, with 5 years of data available from its supervisory control and data acquisition (SCADA) system. Although the Frandsen turbulence estimates deviate in a different manner from measurements at below- and above-rated mean wind speeds, the model remains a conservative approach for fatigue load prediction and reliability. In the current case study, the site-specific assessment using strain gauge measurements yields a 33 % higher annual fatigue reliability index after 35 years compared to a scenario based on the Frandsen estimation combined with ambient environmental data and a generic aeroelastic model. The results also demonstrate that the sensitivity of fatigue reliability to load uncertainty is negligible when load measurements are used directly but relatively high when relying on the Frandsen model in combination with a generic aeroelastic model. Overall, the high variability of the lifetime extension in different scenarios of data availability and accuracy shows the importance and added value of high-quality measurements combined with wind-farm-level SCADA and a model updated in real time (digital twins).

17 WIND ENERGY↗
DOE OSTI2020

Synthesis, Crystal Structure, and Cooperative 3d–5d Magnetism in Rock Salt Type Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6

Two new transition metal oxides with the nominal chemical compositions of Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6 were successfully synthesized. Both compounds crystallize in an ordered rock salt structure type in the monoclinic C2/m space group. The crystal structures were determined using both synchrotron X-ray and time-of-flight neutron, powder diffraction data. In both phases, Ni2+ ions are present while oxidation states of osmium are +6 and +5 in Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6 , respectively. Ni2+ ions in the hypothetical fully ordered phase form a honeycomb arrangement in the ab crystallographic plane and these hexagons are centered by osmium ions. Furthermore, the magnetic layers are separated along the c axis by the octahedra, which are centered by Li+ (or Li+/Ni 2+ , depending on the chemical compositions). Crystal structure refinements reveal that there is some degree of mixed occupancy in cationic positions. Temperature dependent magnetic susceptibility data for both phases show ferrimagnetic transitions with predominant antiferromagnetic (AFM) interactions among 3d electrons of nickel and Sd electrons of osmium. Iso-thermal magnetization loops as a function of the applied magnetic field below the transition temperatures confirm the ferrimagnetic nature in magnetic transitions. Temperature dependent heat capacity data, however, did not exhibit any anomaly in either phase, indicating the absence of long-range magnetic ordering. The lack of long-range order for both Os 5+ and Os 6+ -based compounds was also confirmed by low temperature neutron diffraction data down to 10 K. Temperature dependent AC magnetic susceptibility data in various frequencies for both samples indicate that Li 4 NiOsO 6 exhibits spin-glass-like behavior, while the transition temperature for Li 3 Ni 2 OsO 6 is nearly frequency independent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A search for time-dependent astrophysical neutrino emission with icecube data from 2012 to 2017

High-energy neutrinos are unique messengers of the high-energy universe, tracing the processes of cosmic ray acceleration. This work presents analyses focusing on time-dependent neutrino point-source searches. A scan of the whole sky, making no prior assumption about source candidates, is performed, looking for a space and time clustering of high-energy neutrinos in data collected by the IceCube Neutrino Observatory between 2012 and 2017. No statistically significant evidence for a time-dependent neutrino signal is found with this search during this period, as all results are consistent with the background expectation. Within this study period, the blazar 3C 279, showed strong variability, inducing a very prominent gamma-ray flare observed in 2015 June. This event motivated a dedicated study of the blazar, which consists of searching for a time-dependent neutrino signal correlated with the gamma-ray emission. No evidence for a time-dependent signal is found. Hence, an upper limit on the neutrino fluence is derived, allowing us to constrain a hadronic emission model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine-Learning-Based Mapping and Modeling of Solar Energy with Ultra-High Spatiotemporal Granularity

Despite the rapid growth of solar energy, we still lack a dynamic, high-fidelity database that tracks the spatiotemporal variations of solar PVs and their associated infrastructures across different places at a spatially resolved scale. The absence of such data presents a barrier to various applications such as solar PV growth projection, solar energy integration, solar incentive design, and climate risk assessment. In this project, we aim to bridge this gap by developing AI-based algorithms to extract granular information about solar PV installations and their associated infrastructures (i.e., distribution grids) from widely available unstructured data like remote sensing images and street views. As a result, we have built the Solar Energy Atlas, a fine-grained, large-scale geospatial overlay of distributed solar PVs and distribution grids. On top of it, we have advanced the understanding of solar adoption and distribution grid vulnerability to climate-induced extremes. Our major contributions can be summarized as follow: (1) By developing new AI algorithms, we have built the most comprehensive solar PV spatiotemporal database covering the entire US. This is the first time we obtained the exact GPS locations, size, subtype, and installation year information for rooftop solar PVs across the US. This database can be used for solar PV growth projection, solar energy integration, solar energy policy analysis and design, and spatially-resolved climate risk assessment. (2) Leveraging this database, we have uncovered the socioeconomic driving factors that are correlated with earlier onset of solar adoption and higher saturated adoption levels. We have identified the heterogeneity in the effects of different types of financial incentives on solar adoption and provided implications for tailoring incentive design based on local income levels to promote equitable solar adoption. (3) We have developed a distribution grid GIS mapping algorithm which can obtain granular geospatial and topology information about distribution grids using multi-modal open data, reducing the dependency on hard-to-obtain smart meter data of conventional approaches. It shows effectiveness in both the U.S. and Sub-Saharan Africa. Using this algorithm, we have uncovered the non-uniform vulnerability of distribution grids to wildfires in California in the aspects of undergrounding protection and Distributed Energy Resources (DER) preparedness. This has provided important implications for improving the affordability and equity of grid adaptation approaches. (3) We have made our produced database publicly available and provided user-friendly interface to enable various stakeholders and the general public to interact with the data. We have also integrated the produced data into the Data Commons platform to enable the public to access the data and correlate it with other location-specific characteristics simply using natural language as queries. The impact of our project is three-fold: (1) New algorithms for mapping solar PVs and distribution grids across space and time, which are open source to facilitate researchers and industry; (2) New databases of solar PVs and distribution grids that have been made publicly available for engineering, social, and policy applications; (3) New understandings and actionable insights on the potential approaches to promoting solar adoption and reducing energy infrastructure vulnerabilities. In this report, we start by discussing the project background and motivation (section 5), followed by the overview of project objectives (section 6). Results and discussion for each task are presented in section 7. Significant accomplishments are summarized in section 8. This report will be concluded by discussing the paths forwards (section 9), products (section 10), and team roles (section 11).

14 SOLAR ENERGY↗
DOE OSTI2022

Seismic Waveform Inversion Capability on Resource-Constrained Edge Devices

Seismic full wave inversion (FWI) is a widely used non-linear seismic imaging method used to reconstruct subsurface velocity images, however it is time consuming, has high computational cost and depend heavily on human interaction. Recently, deep learning has accelerated it’s use in several data-driven techniques, however most deep learning techniques suffer from overfitting and stability issues. In this work, we propose an edge computing-based data-driven inversion technique based on supervised deep convolutional neural network to accurately reconstruct the subsurface velocities. Deep learning based data-driven technique depends mostly on bulk data training. In this work, we train our deep convolutional neural network (DCN) (UNet and InversionNet) on the raw seismic data and their corresponding velocity models during the training phase to learn the non-linear mapping between the seismic data and velocity models. The trained network is then used to estimate the velocity models from new input seismic data during the prediction phase. The prediction phase is performed on a resource-constrained edge device such as Raspberry Pi. Raspberry Pi provides real-time and on-device computational power to execute the inference process. In addition, we demonstrate robustness of our models to perform inversion in the presence on noise by performing both noise-aware and no-noise training and feeding the resulting trained models with noise at different signal-to-noise (SNR) ratio values. We make great efforts to achieve very feasible inference times on the Raspberry Pi for both models. Specifically, the inference times per prediction for UNet and InversionNet models on Raspberry Pi were 22 and 4 s respectively whilst inference times for both models on the GPU were 2 and 18 s which are very comparable. Finally, we have designed a user-friendly interactive graphical user interface (GUI) to automate the model execution and inversion process on the Raspberry Pi.

Manu, Daniel (ORCID:0000000154982677)↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spin Seebeck effect near the antiferromagnetic spin-flop transition

In this work, we develop a low-temperature, long-wavelength theory for the interfacial spin Seebeck effect (SSE) in easy-axis antiferromagnets. The field-induced spin-flop (SF) transition of Néel order is associated with a qualitative change in SSE behavior: Below SF, there are two spin carriers with opposite magnetic moments, with the carriers polarized along the field forming a majority magnon band. Above SF, the low-energy, ferromagneticlike mode has magnetic moment opposite the field. This results in a sign change of the SSE across SF, which agrees with recent measurements on Cr 2 O 3 /Pt and Cr 2 O 3 /Ta devices. In our theory, SSE is due to a Néel spin current below SF and a magnetic spin current above SF. Using the ratio of the associated Néel to magnetic spin-mixing conductances as a single constant fitting parameter, we reproduce the field dependence of the experimental data and partially the temperature dependence of the relative SSE jump across SF.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Using data-science approaches to unravel insights for enhanced transport of lithium ions in single-ion conducting polymer electrolyte

Solid polymer electrolytes have yet to achieve the an ionic conductivity > 1 mS/cm at room temperature for realistic applications. This target implies the need to reduce the effective energy barriers of ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer-Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules, 56, 15, 6051(2023)), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify a series of anions leading to the energy barriers less than 30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers less than 20 kJ/mol, we focused on boron-containing SICs, and performed density functional theory (DFT) based calculations to connect the chemical structures via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion transport. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity $> 1$ mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes. The published datasets include all the temperature-dependent ionic conductivity collected from the literature with literature DOIs, DFT calculated binding energies, and python scripts to analyze data, construct statistical models, and generate plots.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using Data-Science Approaches to Unravel Insights for Enhanced Transport of Lithium Ions in Single-Ion Conducting Polymer Electrolytes

Solid polymer electrolytes have yet to achieve the desired ionic conductivity (>1 mS/cm) near room temperature required for many applications. This target implies the need to reduce the effective energy barriers for ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer–Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules 2023, 56, 15, 6051), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify anions leading to the energy barriers <30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers <20 kJ/mol, we have performed density functional theory (DFT) based calculations to connect the chemical structures of boron-based anions via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion hopping. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity >1 mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Trion Formation Resolves Observed Peak Shifts in the Optical Spectra of Transition-Metal Dichalcogenides

Monolayer transition-metal dichalcogenides (ML-TMDs) have the potential to unlock novel photonic and chemical technologies if their optoelectronic properties can be understood and controlled. Yet, recent work has offered contradictory explanations for how TMD absorption spectra change with carrier concentration, fluence, and time. Here, we test our hypothesis that the large broadening and shifting of the strong band-edge features observed in optical spectra arise from the formation of negative trions. We do this by fitting an ab initio based, many-body model to our experimental electrochemical data. Our approach provides an excellent, global description of the potential-dependent linear absorption data. We further leverage our model to demonstrate that trion formation explains the nonmonotonic potential dependence of the transient absorption spectra, including through photoinduced derivative line shapes for the trion peak. Our results motivate the continued development of theoretical methods to describe cutting-edge experiments in a physically transparent way.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data and code for: "Soil respiration response to simulated precipitation change depends on ecosystem type and study duration"

Data, code, and figures from our meta-analysis of the effect of precipitation change on soil respiration, with an emphasis on ecosystem variability and the effect of study duration. We combined the results of 81 separate studies to determine the effect of altered rainfall on the amount of carbon that leaves soil in the form of carbon dioxide, a major component of global carbon cycling. In addition we looked at how long the changes lasted, and how different soil properties and the intensity of precipitation changes affected the study results. We found that more precipitation resulted in greater amounts of carbon dioxide leaving the soil, and less precipitation resulted in less. However, the changes weakened over time in ecosystems that typically receive plenty of rainfall (e.g., forests), in contrast to ecosystems that typically receive little rainfall (e.g., deserts) where changes strengthened over time. Changes in the amount of carbon dioxide that left the soil were also affected by the amount of carbon in the soil, which impacts how much water soil can hold.

54 ENVIRONMENTAL SCIENCES↗