Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Comparative framework”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating environmentally-sensitive tree recruitment in vegetation demographic models

Vegetation demographic models (VDMs) endeavor to predict how global forests will respond to climate change. This requires simulating which trees, if any, are able to recruit under changing environmental conditions. We present a new recruitment scheme for VDMs in which functional-type-specific recruitment rates are sensitive to light, soil moisture and the productivity of reproductive trees. We evaluate the scheme by predicting tree recruitment for four tropical tree functional types under varying meteorology and canopy structure at Barro Colorado Island, Panama. We compare predictions to those of a current VDM, quantitative observations and ecological expectations. We find that the scheme improves the magnitude and rank order of recruitment rates among functional types and captures recruitment limitations in response to variable understory light, soil moisture and precipitation regimes. Furthermore, our results indicate that adopting this framework will improve VDM capacity to predict functional-type-specific tree recruitment in response to climate change, thereby improving predictions of future forest distribution, composition and function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

HRMT-60- RaDIATE Material Studies

HRMT-60 experiment was performed at the CERN-HiRadMat facility in October 2022 to understand thermal shock response of conventional materials and novel materials to support the design and operation of future multi-MW accelerator beam windows and secondary particle-production targets. This experiment, organized within the framework of the RaDIATE collaboration, builds on the previous HRMT-43 (BeGrid2) experiment, where a variety of materials in both non-irradiated and previously proton-irradiated conditions were tested. The primary goal was to understand the failure mechanisms, limits and flow behavior of the various material specimens, as well as compare and contrast the thermal shock response of previously irradiated materials to their non-irradiated counterparts. A total of 120 samples were tested at different beam conditions. This poster will present the preliminary results of several materials tested during this experiments.

Ammigan, K.↗
DOE OSTI2022

Designing Porous Ion Emitters for Thermal Ionization Mass Spectrometry: Evaluating Metal-Organic Frameworks

In this work we describe the first exploration of metal-organic frameworks (MOFs) as ‘next generation’ ion emitters for thermal ionization mass spectrometry (TIMS). MOFs were identified as promising candidates for this application given the synthetic control over their desired structural properties. This tunability results in well-ordered, high surface area, high porosity frameworks with targeted sorption affinities. Here we explored an aluminum-based, bipyridine-containing MOF (MOF-253) with and without incorporating the high work function metal, rhenium (Re). After analysis of a Nd-bearing MOF, we hypothesized that the well-dispersed, sponge-like interconnected network of the degraded structure would enhance Nd ionization more than traditional TIMS loading activators (i.e., phosphoric acid). Compared to filaments loaded with phosphoric acid that require an additional benzene carburization step, the Nd ionization efficiencies (atoms detected relative to atoms loaded) for heated filaments loaded with MOF-253 were similar (~1%). Electron microscopy after TIMS analysis demonstrated that the MOF was retained on the filament and, as hypothesized, preserved sub-micron porosity post analysis. While these results are preliminary, they demonstrate that MOFs have potential to enhance ionization and exceed the performance of traditional loading techniques by forming nano-porous ion emitters (nano-PIEs). Thus, further experimentation is likely to exceed this performance through more specific selection of the base MOF structure and modifications to porosity and composition. This work represents a novel application of MOFs and a next step in the pursuit of advanced thermal ionization with potential to expand across the periodic table.

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

Quantifying Uncertainty in Ice Particle Velocity–Dimension Relationships Using MC3E Observations

Abstract Ice particle terminal fall velocity ( V t ) is fundamental for determining microphysical processes, yet remains extremely challenging to measure. Current theoretical best estimates of V t are functions of Reynolds number. The Reynolds number is related to the Best number, which is a function of ice particle mass, area ratio ( A r ), and maximum dimension ( D max ). These estimates are not conducive for use in most models since model parameterizations often take the form , where ( α , β ) depend on habit and D max . A previously developed framework is used to determine surfaces of equally plausible ( α , β ) coefficients whereby ice particle size/shape distributions are combined with V t best estimates to determine mass- ( V M ) or reflectivity-weighted ( V Z ) velocities that closely match parameterized V M ,SD or V Z ,SD calculated using the ( α , β ) coefficients using two approaches. The first uses surfaces of equally plausible ( a , b ) coefficients describing mass ( M )–dimension relationships (i.e., ) to calculate mass- or reflectivity-weighted velocity from size/shape distributions that are then used to determine ( α , β ) coefficients. The second investigates how uncertainties in A r , D max , and size distribution N ( D ) affect V M or V Z . For seven of nine flight legs flown on 20 and 23 May 2011 during the Mesoscale Continental Convective Clouds Experiment (MC3E), uncertainty from natural parameter variability—namely, the variability in ice particle parameters in similar meteorological conditions—exceeds uncertainties arising from different A r assumptions or D max estimates. The combined uncertainty between A r , D max , and N ( D ) produced smaller variability in ( α , β ) compared to varying M ( D ), demonstrating M ( D ) must be accurately quantified for model fall velocities. Primary sources of uncertainty vary considerably depending on environmental conditions. Significance Statement Ice particle fall velocity is fundamental for numerous processes within clouds, and hence is a critical property that must be accurately represented in weather and climate models. Using aircraft observations of ice particle shapes and sizes obtained in clouds behind midlatitude thunderstorms, this work develops a new framework for estimating ice particle fall velocities and their uncertainty, including quantifying the importance of different uncertainty sources from cloud microphysics measurements. Natural parameter variability contributes the most uncertainty in ice particle fall velocity estimates, although other sources can also be important contributors to uncertainty in certain conditions. Additional work examining ice particle data is needed to further understand how dependent uncertainty in certain ice particle properties are to local environmental conditions.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Isomerization and Selective Hydrogenation of Propyne: Screening of Metal–Organic Frameworks Modified by Atomic Layer Deposition

Various metal oxide clusters upward of 8 atoms (Cu, Cd, Co, Fe, Ga, Mn, Mo, Ni, Sn, W, Zn, In, and Al) were incorporated into the pores of the metal–organic framework (MOF) NU-1000 via atomic layer deposition (ALD) and tested via high-throughput screening for catalytic isomerization and selective hydrogenation of propyne. Cu and Co were found to be the most active for propyne hydrogenation to propylene, and synergistic bimetallic combinations of Co and Zn, along with standalone Zn and Cd, were established as the most active for conversion to the isomerized product, propadiene. The combination of Co and Zn in NU-1000 diminished the propensity for full hydrogenation to propane as well as coking compared to its individual components. Finally, this study highlights the potential for high-throughput screening to survey monometallic and bimetallic cluster combinations that best affect the efficient transformation of small molecules, while discerning mechanistic differences in isomerization and hydrogenation by different metals.

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

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

SymbolFit: Automatic Parametric Modeling with Symbolic Regression

We introduce SymbolFit (API: https://github.com/hftsoi/symbolfit), a framework that automates parametric modeling by using symbolic regression to perform a machine-search for functions that fit the data while simultaneously providing uncertainty estimates in a single run. Traditionally, constructing a parametric model to accurately describe binned data has been a manual and iterative process, requiring an adequate functional form to be determined before the fit can be performed. The main challenge arises when the appropriate functional forms cannot be derived from first principles, especially when there is no underlying true closed-form function for the distribution. In this work, we develop a framework that automates and streamlines the process by utilizing symbolic regression, a machine learning technique that explores a vast space of candidate functions without requiring a predefined functional form because the functional form itself is treated as a trainable parameter, making the process far more efficient and effortless than traditional regression methods. We demonstrate the framework in high-energy physics experiments at the CERN Large Hadron Collider (LHC) using five real proton-proton collision datasets from new physics searches, including background modeling in resonance searches for high-mass dijet, trijet, paired-dijet, diphoton, and dimuon events. We show that our framework can flexibly and efficiently generate a wide range of candidate functions that fit a nontrivial distribution well using a simple fit configuration that varies only by random seed, and that the same fit configuration, which defines a vast function space, can also be applied to distributions of different shapes, whereas achieving a comparable result with traditional methods would have required extensive manual effort.

Tsoi, Ho Fung [Univ. of Pennsylvania, Philadelphia↗
DOE OSTI2026

Multiphysics simulation of TRISO fuel compacts and the effects of homogenization on silver release predictions

This work studies the impact of explicit and homogenized modeling approaches on the multiphysics simulation of TRistructural ISOtropic (TRISO) fuel compacts in prismatic High Temperature Gas Reactors (HTGRs) and the silver release predictions. TRISO fuel compacts exhibit complex double heterogeneity that significantly affects heat conduction, neutron transport, and fission products release. In this work, we use Cardinal, a multiphysics tool based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to couple neutron transport and heat conduction. OpenMC is used for neutron transport and the MOOSE heat transfer module is used for the heat conduction. Then, we use BISON for the silver release predictions. Two different modeling approaches—explicit modeling of individual TRISO particles and homogenized representation using effective thermal properties—are compared at high TRISO packing fractions (20% and 40%) across varying power densities. Results demonstrate that homogenization significantly under-predicts peak temperatures, for the case of high power/TRISO, there is a difference of 83.73 K in the maximum temperature. Additionally, homogenization underestimates the silver release predictions compared to explicit modeling. Specifically, the temperature and power density differences lead to significant differences in silver release predictions. This work demonstrates the importance of accurately modeling heterogeneity of the TRISO particles to reliably predict fission product release and assess reactor safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Hierarchical Truncations for Many-Body Expansion Potentials

In this work, a new strategy to truncate high-order terms in the many-body expansion (MBE) is proposed. This new approach, which we call a hierarchical many-body expansion (HMBE), is based on a hierarchical partition of the system into multitier fragments and can in principle be applied to any large molecular system. Numerical tests on a series of (H 2 O) 64 structures are presented, demonstrating satisfactory relative energies between the clusters and binding energies of individual clusters compared with full-cluster calculations, with significantly fewer high-order terms computed than conventional MBE. The hierarchical truncation can be augmented by certain many-body terms for fragments at the interface between the partitions (called “Schengen terms”) to further improve accuracy. This work establishes the HMBE scheme as a promising framework to model very large systems (e.g., proteins), which are naturally built on a hierarchical structure.

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

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A quantitative comparison of the fingerprint of twinned microstructures through surface and three-dimensional techniques

Assessing the fingerprint of a material’s microstructure is key for supporting materials design. With the emergence of a wide range of 3D characterization techniques, it is critical to understand the main differences in fingerprints reconstructed from 2D and 3D datasets. To this end, we introduce a graph-based microstructure reconstruction framework that enables structural comparisons of twin domain networks in high purity Ti using 3D and 2D electron backscatter diffraction. Insights into the structure of the twin networks are facilitated by combining statistical analysis of twin crystallography with visual and graphical analysis of the novel graph abstractions of the twins. We demonstrate that compared to 3D reconstructions, conventional 2D views of twinning miss key aspects of the microstructure including the high interconnectivity of domains into networks that span the full reconstruction volume. The reduced cross-grain and in-grain twin connectivity typically observed in 2D has notable implications on our understanding of how twinning mediates the plastic response of microstructures and how twin networks evolve. It is thus clear that 3D characterization is critical for accurately inferring both twin network morphologies as well as the key unit processes facilitating network formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Modular and Transferable Reinforcement Learning Framework for the Fleet Rebalancing Problem

Mobility on demand (MoD) systems show great promise in realizing flexible and efficient urban transportation. However, significant technical challenges arise from operational decision making associated with MoD vehicle dispatch and fleet rebalancing. For this reason, operators tend to employ simplified algorithms that have been demonstrated to work well in a particular setting. To help bridge the gap between novel and existing methods, we propose a modular framework for fleet rebalancing based on model-free reinforcement learning (RL) that can leverage an existing dispatch method to minimize system cost. In particular, by treating dispatch as part of the environment dynamics, a centralized agent can learn to intermittently direct the dispatcher to reposition free vehicles and mitigate against fleet imbalance. We formulate RL state and action spaces as distributions over a grid partitioning of the operating area, making the framework scalable and avoiding the complexities associated with multiagent RL. Numerical experiments, using real-world trip and network data, demonstrate that RL reduces waiting time by 28% to 38% for the same-day evaluation, 17% to 44% for cross-day evaluation, and 22% to 25% for cross-season evaluation compared with no rebalancing scenarios. This approach has several distinct advantages over baseline methods including: improved system cost; high degree of adaptability to the selected dispatch method; and the ability to perform scale-invariant transfer learning between problem instances with similar vehicle and request distributions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Convergence and divergence in science and practice of urban and rural forest restoration

Forest restoration has never been higher on policymakers' agendas. Complex and multi-dimensional arrangements across the urban–rural continuum challenge restorationists and require integrative approaches to strengthen environmental protection and increase restoration outcomes. It remains unclear if urban and rural forest restoration are moving towards or away from each other in practice and research, and whether comparing research outcomes can help stakeholders to gain a clearer understanding of the interconnectedness between the two fields. This study aims to identify the challenges and opportunities for enhancing forest restoration in both urban and rural systems by reviewing the scientific evidence, engaging with key stakeholders and using an urban–rural forest restoration framework. Using the Society for Ecological Restoration's International Principles as discussion topics, we highlight aspects of convergence and divergence between the two fields to broaden our understanding of forest restoration and promote integrative management approaches to address future forest conditions. Our findings reveal that urban and rural forest restoration have convergent and divergent aspects. We emphasise the importance of tailoring goals and objectives to specific contexts and the need to design different institutions and incentives based on the social and ecological needs and goals of stakeholders in different regions. Here we discuss the challenges of achieving high levels of ecological restoration and the need to go beyond traditional ecology to plan, implement, monitor, and adaptively manage restored forests. We suggest that rivers and watersheds could serve as a common ground linking rural and urban landscapes and that forest restoration could interact with other environmental protection measures. We note the potential for expanding the creative vision associated with increasing tree-containing environments in cities to generate more diverse and resilient forest restoration outcomes in rural settings. This study underscores the value of integrative management approaches in addressing future forest conditions across the urban–rural continuum. Our framework provides valuable insights for policymakers, researchers, and decision-makers to advance the field of forest restoration and address the challenges of restoration across the urban–rural continuum. The rural–urban interface serves as a convergence point for forest restoration, and both urban and rural fields can benefit from each other's expertise.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Substitution or Shared Utilization? Intrahousehold Vehicle Use in Mixed-Powertrain Households

While previous research has focused heavily on understanding the factors deriving alternative fuel vehicle adoption rates, there remains a significant gap in understanding how households distribute mileage across different powertrains. This study utilizes data from the 2022 Next Generation National Household Travel Survey to investigate vehicle miles traveled within a sample of 150 plug-in electric vehicle (PEV)-owning households (in which at least one battery electric vehicle is present), characterizing how different powertrains are integrated into daily mobility. Leveraging a Seemingly Unrelated Regression (SUR) framework the study jointly models the utilization of PEVs, hybrid electric vehicles (HEV), and internal combustion engine vehicles (ICEVs) while accounting for household-level substitution effects. The results provide evidence of an asymmetric substitution effect. In households with mixed-powertrain configurations, the ICEV captures a substantially higher share of household miles (compared with the PEV), acting as a utility sponge. Conversely, the model identifies specific socioeconomic and geographic cohorts that prioritize PEV as the primary household workhorse, indicating a systematic sorting effect. Although the sample size limits broader generalizability, these findings suggest that PEVs are used for frequent, specific routine-intensive roles, whereas the ICEV remains a specialized utility vehicle. These insights highlight distinct intrahousehold vehicle use behaviors that are often obscured by aggregate fleetwide statistics.

25 ENERGY STORAGE↗
DOE OSTI2022

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces – part II

The modeling of the near-field interaction in the scattering-type scanning near-field optical microscope (s-SNOM) is rapidly advancing, although an accurate yet versatile modeling framework that can be easily adapted to various complex situations is still lacking. In this work, we propose a time-efficient numerical scheme in the quasi-electrostatic limit to capture the tip-sample interaction in the near field. This method considers an extended tip geometry, which is a significant advantage compared to the previously reported method based on the point-dipole approximation. Using this formalism, we investigate, among others, nontrivial questions such as uniaxial and biaxial anisotropy in the near-field interaction, the relationship between various experimental parameters (e.g. tip radius, tapping amplitude, etc.), and the tip-dependent spatial resolution. The demonstrated method further sheds light on the understanding of the contrast mechanism in s-SNOM imaging and spectroscopy, while also representing a valuable platform for future quantitative analysis of the experimental observations.

Chen, Xinzhong (ORCID:0000000261033848)↗
DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bayesian Spatial Models for Projecting Corn Yields

Climate change is predicted to impact corn yields. Previous studies analyzing these impacts differ in data and modeling approaches and, consequently, corn yield projections. We analyze the impacts of climate change on corn yields using two statistical models with different approaches for dealing with county-level effects. The first model, which is novel to modeling corn yields, uses a computationally efficient spatial basis function approach. We use a Bayesian framework to incorporate both parametric and climate model structural uncertainty. We find that the statistical models have similar predictive abilities, but the spatial basis function model is faster and hence potentially a useful tool for crop yield projections. We also explore how different gridded temperature datasets affect the statistical model fit and performance. Compared to the dataset with only weather station data, we find that the dataset composed of satellite and weather station data results in a model with a magnified relationship between temperature and corn yields. For all statistical models, we observe a relationship between temperature and corn yields that is broadly similar to previous studies. We use downscaled and bias-corrected CMIP5 climate model projections to obtain detrended corn yield projections for 2020–2049 and 2069–2098. In both periods, we project a decrease in the mean corn yield production, reinforcing the findings of other studies. However, the magnitude of the decrease and the associated uncertainties we obtain differ from previous studies.

54 ENVIRONMENTAL SCIENCES↗