Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed estimation”

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 487 records · Page 27

The PAU Survey: narrow-band photometric redshifts using Gaussian processes

Here, we study the performance of the hybrid template machine learning photometric redshift (photo- z ) algorithm delight , which uses Gaussian processes, on a subset of the early data release of the Physics of the Accelerating Universe Survey (PAUS). We calibrate the fluxes of the 40 PAUS narrow bands with six broad-band fluxes ( uBVriz ) in the Cosmic Evolution Survey (COSMOS) field using three different methods, including a new method that utilizes the correlation between the apparent size and overall flux of the galaxy. We use a rich set of empirically derived galaxy spectral templates as guides to train the Gaussian process, and we show that our results are competitive with other standard photometric redshift algorithms. delight achieves a photo- z 68th percentile error of σ 68 = 0.0081(1 + z ) without any quality cut for galaxies with i auto < 22.5 as compared to 0.0089(1 + z ) and 0.0202(1 + z ) for the bpz and annz 2 codes, respectively. delight is also shown to produce more accurate probability distribution functions for individual redshift estimates than bpz and annz 2. Common photo- z outliers of delight and bcnz 2 (previously applied to PAUS) are found to be primarily caused by outliers in the narrow-band fluxes, with a small number of cases potentially indicating spectroscopic redshift failures in the reference sample. In the process, we introduce performance metrics derived from the results of bcnz 2 and delight , allowing us to achieve a photo- z quality of σ 68 < 0.0035(1 + z ) at a magnitude of i auto < 22.5 while keeping 50 per cent objects of the galaxy sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Stimulated emission tomography for efficient characterization of spatial entanglement

Stimulated-emission tomography (SET) is an excellent tool for characterizing the process of spontaneous parametric down conversion (SPDC), which is commonly used to create pairs of entangled photons for use in quantum information protocols. The use of stimulated emission increases the average number of detected photons by several orders of magnitude compared to the spontaneous process. In a SET measurement, the parametric down conversion is seeded by an intense signal field prepared with specified mode properties rather than by broadband multimodal vacuum fluctuations, as is the case for the spontaneous process. The SET process generates an intense idler field in a mode that is the complex conjugate to the signal mode. In this work we use SET to estimate the joint spatial mode distribution (JSMD) in the Laguerre-Gaussian (LG) basis of the two photons of an entangled photon pair. The pair is produced by parametric down conversion in a beta barium borate (BBO) crystal with type-II phase matching pumped at a wavelength of 405 nm along with a 780 nm seed signal beam prepared in a variety of LG modes to generate an 842 nm idler beam of which the spatial mode distribution is measured. We observe strong idler production and good agreement with the theoretical prediction of its spatial mode distribution. Our experimental procedure should enable the efficient determination of the photon-pair wavefunctions produced by low-brightness SPDC sources and the characterization of high-dimensional entangled-photon pairs. Published by the American Physical Society 2024

Xu, Yang (ORCID:0009000954547320)↗

Doubling protected land area may be inefficient at preserving the extent of undeveloped land and could cause substantial regional shifts in land use

Abstract Projection of land use and land‐cover change is highly uncertain yet drives critical estimates of carbon emissions, climate change, and food and bioenergy production. We use new, spatially explicit land availability data in conjunction with a model sensitivity analysis to estimate the effects of additional land protection on land use and land cover. The land availability data include protected land and agricultural suitability and is incorporated into the Moirai land data system for initializing the Global Change Analysis Model. Overall, decreasing land availability is relatively inefficient at preserving undeveloped land while having considerable regional land‐use impacts. Current amounts of protected area have little effect on land and crop production estimates, but including the spatial distribution of unsuitable (i.e., unavailable) land dramatically shifts bioenergy production from high northern latitudes to the rest of the world, compared with uniform availability. This highlights the importance of spatial heterogeneity in understanding and managing land change. Approximately doubling the current protected area to emulate a 30% protected area target may avoid land conversion by 2050 of less than half the newly protected extent while reducing bioenergy feedstock land by 10.4% and cropland and grazed pasture by over 3%. Regional bioenergy land may be reduced (increased) by up to 46% (36%), cropland reduced by up to 61%, pasture reduced by up to 100%, and harvested forest reduced by up to 35%. Only a few regions show notable gains in some undeveloped land types of up to 36%. Half of the regions can reach the target using only unsuitable land, which would minimize impacts on agriculture but may not meet conservation goals. Rather than focusing on an area target, a more robust approach may be to carefully select newly protected land to meet well‐defined conservation goals while minimizing impacts to agriculture.

54 ENVIRONMENTAL SCIENCES↗

A Review of Doses for Dental Imaging in 2010–2020 and Development of a Web Dose Calculator

Dental imaging is one of the most common types of diagnostic radiological procedures in modern medicine. We introduce a comprehensive table of organ doses received by patients in dental imaging procedures extracted from literature and a new web application to visualize the summarized dose information. We analyzed articles, published after 2010, from PubMed on organ and effective doses delivered by dental imaging procedures, including intraoral radiography, panoramic radiography, and cone-beam computed tomography (CBCT), and summarized doses by dosimetry method, machine model, patient age, and technical parameters. Mean effective doses delivered by intraoral, 1.32 (0.60–2.56) μSv, and panoramic, 17.93 (3.47–75.00) μSv, procedures were found to be about1% and 15% of that delivered by CBCT, 121.09 (17.10–392.20) μSv, respectively. In CBCT imaging, child phantoms received about 29% more effective dose than the adult phantoms received. The effective dose of a large field of view (FOV) (>150 cm2) was about 1.6 times greater than that of a small FOV (<50 cm2). The maximum CBCT effective dose with a large FOV for children, 392.2 μSv, was about 13% of theeffective dose that a person receives on average every year from natural radiation, 3110 μSv. Monte Carlo simulations of representative cases of the three dental imaging procedures were then conducted to estimate and visualize the dose distribution within the head. The user-friendly interactive web application (available at http://dentaldose.org) receives user input, such as the number of intraoral radiographs taken, and displays total organ and effective doses, dose distribution maps, and a comparison with other medical and natural sources of radiation. The web dose calculator provides a practical resource for patients interested in understanding the radiation doses delivered by dental imaging procedures.

Lee, Hawon↗

Single-Logarithmic Corrections to Small-$x$ Helicity Evolution

The small- x x quark helicity evolution equations at double-logarithmic order, with the kernel \sim\alpha_s\ln^2(1/x) ∼ α s ln 2 ( 1 / x ) , have been derived previously. In this work, we derive the single-logarithmic corrections to the equations, to order \alpha_s\ln(1/x) α s ln ( 1 / x ) of the evolution kernel. The new equations include the effects of the running coupling and the unpolarized small- x x evolution, both of which are parametrically significant at single-logarithmic order. The large- N_c N c and large- N_c\ N_f N c N f approximations to the equation are computed. (Here, N_c N c and N_f N f are the numbers of quark colors and flavors, respectively.) Their solutions will provide more precise estimates of the quark helicity distribution at small x x , contributing to the resolution of the proton spin puzzle.

Tawabutr, Yossathorn↗

2019 Cost of Wind Energy Review

This report uses representative utility-scale and distributed wind energy projects to estimate the levelized cost of energy (LCOE) for land-based and offshore wind power plants in the United States. Data and results detailed here are derived from 2019 commissioned plants and representative industry data as well as state-of-the-art modeling capabilities. Modeling is conducted to provide more granular detail on specific cost categories. This report represents the ninth annual installment and is intended to provide insight into current component-level costs as well as a basis for understanding variability in wind energy LCOE across the country.

17 WIND ENERGY↗

Thoughts on Analyzing Neutron Multiplicity Data

Bayesian methods offer many advantages for analyzing neutron multiplicity data. Ideally one would make nonparametric estimates of the the posterior distributions, the probabilities - given the data - for unknowns to take on particular values. Once one has the posterior, uncertainty quantification is very natural. We describe how different ways of formulating the statistical theory of fission chains are mathematically equivalent. When confronted with actual data and the burden of finite statistics, different methods of expressing the theory offer advantages for understanding the underlying source parameters from an information theory standpoint.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impacts of Regional Air Mobility and Electrified Aircraft on Local Grid Infrastructure and Airport Electricity Demand

The goal of this work is to analyze the electrical infrastructure that may be necessary to serve electrified aviation at a subset of airports where potential electrified flight demand has been provided. Additionally, this work will estimate the amount of onsite distributed energy resources (i.e., solar photovoltaics and battery energy storage systems) that could be used to serve EA in cost effective scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

2022 Cost of Wind Energy Review [Slides]

The 12th annual Cost of Wind Energy Review, now presented in slide deck format, uses representative utility-scale and distributed wind energy projects to estimate the levelized cost of energy (LCOE) for land-based and offshore wind power plants in the United States.

17 WIND ENERGY↗

PR100: Estimated Medium- and Heavy-Duty Electric Vehicle Adoption and Load Estimation in Puerto Rico through 2050

The 2-year Puerto Rico Grid Resilience and Transition to 100% Renewable Energy Study analyzed stakeholder-driven pathways to Puerto Rico’s clean energy future. Outputs relating to electricity demand modeling were partially informed by estimates of electric vehicle adoption across all classes of medium- and heavy-duty vehicles (MHDVs), and the ensuing charging loads. To create these estimates, the team developed a transportation model for MHDVs in Puerto Rico to estimate the amount and geospatial distribution of energy used. Charging schedules for the different end uses of MHDVs were then used to construct electric load shapes assuming a portion of those vehicles would be replaced by battery electric counterparts. Study results showed that, by 2050, electric vehicles may constitute roughly 50% of the MHDV population in Puerto Rico. The resulting electrical demand curve attributable to MHDV charging showed that, for solar energy-based electrical systems with limited energy storage, this demand may create challenges unless appropriately managed either on the demand or supply side.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Tuning Shinkarev’s Bicycle: Separating the Parallel Cycles of Photosystem II Using Empirical Wavelet Transform

The oxygen-evolving complex (OEC) of Photosystem II (PSII) catalyzes light-driven water oxidation, a process necessary to sustain Earth’s atmospheric oxygen. Oxygen yields measured during single-turnover flash sequences exhibit period-four oscillations, which form the basis of the Joliot–Kok (S-state) model. However, when the oscillations of other processes contribute to the measured oxygen yield, fitting methods can conflate these signals and distort estimates of inefficiencies and initial S-state populations. To address this, we applied the empirical wavelet transform (EWT) as a model-independent method to separate overlapping oscillators and capture damping dynamics that are not well represented in Fourier analysis. We tested this framework on polarographic flash-oxygen traces from both our Synechocystis sp. PCC 6803 thylakoid membrane preparations and archival datasets on Chlorella and isolated chloroplasts. EWT consistently resolves the expected period-four component alongside a distinct binary oscillation. Simulations suggest that fitting this isolated period-four signal recovers VZAD parameters more accurately than analysis of raw traces, yielding different estimates for S-state distributions and transition probabilities. Notably, this binary oscillation aligns closely with semiquinone dynamics predicted solely from period-four fit parameters. These findings indicate that EWT can effectively distinguish complex signals in oxygen evolution, offering a framework potentially applicable to other spectroscopic probes of the S-state cycle.

Ferrari, Nicholas [Louisiana State Univ., Baton Ro↗

Clumping index estimation with 30°-tilted cameras in row crops: Evaluation of methods and segment size effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Furthermore, applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

30°-tilted camera↗

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

Model ensembles are effective tools for estimating prediction uncertainty in deep learning atomistic force fields. However, their widespread adoption is hindered by high computational costs and overconfident error estimates. In this work, we address these challenges by leveraging distributions of per-sample errors obtained during training and employing a distance-based similarity search in the model latent space. Our method, which we call LTAU (Loss Trajectory Analysis for Uncertainty), efficiently estimates the full probability distribution function of errors for any test point using the logged training errors, achieving speeds that are 2–3 orders of magnitudes faster than typical ensemble methods and allowing it to be used for tasks where training or evaluating multiple models would be infeasible. We apply LTAU towards estimating parametric uncertainty in atomistic force fields (LTAU-FF), demonstrating that it produces well-calibrated confidence intervals and predicts errors that correlate strongly with the true errors for data near the training domain. Furthermore, we show that the errors predicted by LTAU-FF can be used in practical applications for detecting out-of-domain data, tuning model performance, and predicting failure during simulations. We believe that LTAU will be a valuable tool for uncertainty quantification in atomistic force fields and is a promising method that should be further explored in other domains of machine learning.

97 MATHEMATICS AND COMPUTING↗

Ensuring Solution Uniqueness in Three-Phase Power System State Estimation

This paper is concerned with the issue of potential non-unique solutions in three-phase state estimation. Theory of observability analysis for positive sequence power system state estimation is based on certain assumptions that avoid possibility of multiple solutions. Also, it is shown that observability of a positive sequence network remains independent of the network parameters or the operating state. When extending single-phase observability analysis directly to the three-phase case, this paper considers the possibility of converging to multiple solutions, i.e. solution non-uniqueness, even for cases where state estimator successfully converges. The study illustrates via numerical examples the likelihood of converging to entirely different solutions for certain network parameters. It also examines how the operating state, particularly under unbalanced loading, leads to solution non-uniqueness. The paper then describes an alternative approach to ensure a unique solution in three-phase state estimation. This method aims to accurately and uniquely estimate the state of any unbalanced three-phase system, irrespective of load imbalance, network configuration, existence of synchronous generators or transformers.

Power System State Estimation, Three-Phase, Distri↗

Probabilistic Analysis of Long-Term Degradation of Microwave Cavity Flow Sensor

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. We have performed a proof-of-concept experiment of flow sensing with the transducer in liquid sodium at 340°C in impinging liquid jet geometry. The transducer remained in liquid sodium for 70 days. After removal, no structural damage was observed, and the expected transducer response was verified in a water test. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance. Since a yellowish residue was observed on the transducer surface after 70 days of immersion in liquid sodium, we have investigated possible evidence of corrosion. Chromium depletion is a typical indicator of the corrosion process in stainless steel. Scraping off a residue from the transducer and performing scanning electron microscopy (SEM) with energy dispersive analysis (EDS) did not find any chromium in the residue. Approximately 60% of the residue consisted of copper, which can be attributed to contamination of sodium due to powder residue from machining of copper and brass components of the transducer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Final Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Chemistry imaging and distribution analysis of rare earth elements in coal using LIBS and LA-ICP-MS instruments

Currently, demand for rare earth elements (REEs) increased significantly. Coal is actively evaluated as potential economic sources for extraction of REEs. Here, in this work, laser-induced breakdown spectroscopy (LIBS) was evaluated for rapid estimation of REEs content and their distribution in the natural coal samples. The results were compared with similar laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) measurements. Thirteen coal samples (nine standard samples and five natural samples) were used in this study. Powder samples were pressed into pellets while coal chunks were directly ablated for data recording. Pellets of the powder standard samples were used to optimize the data acquisition system and then data recorded with this optimized system was used to identify the proper data acquisition and analysis models. After establishing the proper data acquisition system and analysis model using the standard samples, natural coal samples in powder form and their chunks were utilized to record LIBS and LA-ICP-MS spectra. Multivariate calibration models were developed using four of the natural samples, which were evaluated by predicting the REE content in the fifth sample. Principal component analysis was performed on the LIBS data obtained from the natural samples and it classified all the samples with high accuracy. Two-dimensional (2D) elemental mapping on coal chunk samples was also performed using both LIBS and LA-ICP-MS to study the distribution of REEs in the samples. The resulting elemental images and their correlations can be used to infer mineral distributions.

01 COAL, LIGNITE, AND PEAT↗

Application of Point Precipitation Frequency Estimates to Watersheds

This report documents work sponsored by the U.S. Nuclear Regulatory Commission (NRC) at the Oak Ridge National Laboratory (ORNL) as part of the RES project, “Application of Point Precipitation Frequency Estimates to Watersheds.” This project was implemented as part of the Probabilistic Flood Hazard Assessment (PFHA) Research Program. The objective of the PFHA Research Program is to develop tools and guidance on the use of PFHA methods to risk-inform NRC’s licensing of new facilities as well as licensing and oversight of currently operating facilities as they relate to flooding hazards. Many nuclear power plants (NPPs) are located on or near rivers so riverine flooding hazards need to be considered in their design and operation. Probabilistic riverine flood models are important tools for realistic assessment of flooding risks. However, these models require areal estimates of the depth, duration, and frequency of rainfall distributed over the watershed, which are not often available. Point precipitation frequency estimates are more widely available. For example, the National Oceanic and Atmospheric Administration (NOAA) has published NOAA Atlas 14, which provides point precipitation frequency estimates for 5-minute through 60-day durations at average recurrence intervals of 1-year through 1,000-year. The research documented in this report addresses areal reduction factors (ARFs), which can be used to convert the widely available point precipitation frequency estimates, to estimates of areal precipitation frequency over a watershed. The most widely used ARF source is Technical Paper 29 (TP-29) published by the then U.S. Weather Bureau in 1958. However, both the methods and the underlying precipitation data used to produce TP-29 are seriously out of date. For example, due to the small gauge network available at the time of TP-29’s compilation, ARF estimates developed are only for watersheds smaller than about 400 square miles. Due to the relatively short record lengths of precipitation data available, frequency considerations could not be accurately determined. Other factors such as regional climate and seasonality were not addressed. Several newer methods have been published since TP-29 was developed and both the type and quantity of precipitation data have increased significantly, along with computational resources and analytical tools such as geographic information systems. This report reviewed and assessed the available precipitation products and methods for conducting ARF analysis. The work applied up-to-date precipitation data products and analysis methods with a novel watershed-based approach to investigate how ARF estimates vary across different methods, data sources, geographical locations, return periods, and seasons. The overall findings reported here regarding basic ARF trends are in line with other recent studies showing that ARFs decrease with increasing area, increase with increasing duration, and decrease with increasing return period. This study found significant differences among the available ARF methods. This work also found a strong geographical variability across different US hydrologic regions, suggesting that the ARF are specific to regional climate patterns and geographical characteristics and should not be applied arbitrarily to other locations. The results also reveal the importance of data record length, especially for high return level ARFs. The work reported in NUREG/CR-7271 will assist NRC staff in assessing different classes of ARF methods in conjunction with available rainfall data sets. It will also support the development of guidance for application of point precipitation data in PFHAs. It should be noted that the ARF values presented in this report for any location or region were developed for the purposes of comparing methods and investigating the factors that influence ARFs. They should not be considered official and should not be used in leu of a site-specific analysis.

54 ENVIRONMENTAL SCIENCES↗