Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable conductance”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Implications of increased spatial and trophic overlap between juvenile Pacific salmon and Sablefish in the northern California Current

Abstract Objective The study was designed to assess long-term variability in the distribution of juvenile Pacific salmon Oncorhynchus spp. and Sablefish Anoplopoma fimbria. The study also evaluated whether Sablefish and Pacific salmon shared food resources and looked to characterize Sablefish during an understudied period of their life cycle. Methods To meet the objectives, the study used data from 26 years of surface trawls conducted in Oregon and Washington coastal waters (1998–2023). Spatial–temporal models were used to measure changes in abundance and distribution of Pacific salmon and Sablefish along with covariates of ocean temperature. The study evaluated trophic characteristics of Pacific salmon and Sablefish from 2020 for differences. The temporal variation in size and diets of Sablefish were also analyzed, along with energy density of fish caught in 2020. Result The spatial–temporal model demonstrated that there has been a nearshore expansion of juvenile Sablefish over the past 26 years that was correlated with increased ocean temperature. The nearshore expansion of Sablefish resulted in increased spatial and trophic overlap with juvenile Pacific salmon. While feeding in nearshore waters, juvenile Sablefish demonstrated competitive feeding advantages over juvenile Pacific salmon during a critical phase of salmonid early marine life history. Juvenile Sablefish exhibited significant ontogenetic diet and energetic shifts, and even the smallest (68–80 mm fork length) were piscivorous. Conclusions If juvenile Sablefish numbers continue to increase relative to Pacific salmon, they could exert more competitive pressure, especially if food resources become limited. Pacific salmon may experience adverse effects from competition, regardless of whether or not juvenile Sablefish, which have recently expanded into nearshore waters, successfully recruit to the adult population.

Daly, Elizabeth A. (ORCID:0000000195334457)↗
DOE OSTI2024

Modeling the Impacts of Material Properties on Oscillatory Neuron Behavior

In this study, neuromorphic computing, which mimics the functions of biological brains, offers improvements in both latency and energy efficiency over typical von Neumann computing architectures. Spiking neural networks can be especially power-efficient because they encode information temporally and can use more sparse electrical inputs. Here, we study the design of volatile memristors (variable resistors with memory) for neuronal devices, with particular consideration toward the feasibility of all-on-chip oscillation using built-in capacitance. We use circuit simulations to model the behavior of oscillator neurons with a range of realistic material properties. We find that energy inputs increase with insulating-phase resistivity, thermal conductivity, and device aspect ratio. However, we also find that the minimum capacitance needed for oscillation decreases with increasing insulating-phase resistivity, which opposes the constraints for power efficiency. Based on published data on NbO 2 , VO 2 , and EuNiO 3 , we find that existing materials can be engineered for all-on-chip spiking using their parasitic capacitance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Preliminary Design of Engineering-Scale Salt Accident Analysis Facility to Support Molten Salt Reactor Licensing

Systems-level nuclear accident analysis codes for reactor licensing must be validated using experimental data that represent behaviors expected during actual full-scale accidents. Some behaviors may arise from coupled processes and only manifest at large scales. This report presents the preliminary design of the Salt Accident Analysis Facility (SAAF, pronounced “safe”), which is an experimental test facility to be constructed at Argonne that can be used to conduct integrated salt accident tests at an engineering scale. The measurement capabilities of the SAAF are based on previously developed methods and will provide the representative datasets that are needed to support molten salt reactor (MSR) licensing. Details of the design, the processes to be quantified, the measurement techniques for quantifying the processes, the variables that can be adjusted to simulate different accident scenarios, and operational considerations are presented herein. This report provides stakeholders the opportunity to give feedback on the test facility capabilities and planned analyses before it is constructed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Laboratory Heat Flux Estimates of Seawater Foam for Low Wind Speeds

Laboratory experiments were conducted to measure the heat flux from seafoam continuously generated in natural seawater. Using a control volume technique, heat flux was calculated from foam and foam-free surfaces as a function of ambient humidity (ranged from 40% to 78%), air–water temperature difference (ranged from -9 °C to 0 °C), and wind speed (variable up to 3 m s-1). Water-surface skin temperature was imaged with a calibrated thermal infrared camera, and near-surface temperature profiles in the air, water, and foam were recorded. Net heat flux from foam surfaces increased with increasing wind speed and was shown to be up to four times greater than a foam-free surface. The fraction of the total heat flux due to the latent heat flux was observed for foam to be 0.75, with this value being relatively constant with wind speed. In contrast, for a foam-free surface the fraction of the total heat flux due to the latent heat flux decreased at higher wind speeds. Temperature profiles through foam are linear and have larger gradients, which increased with wind speed, while foam free surfaces show the expected logarithmic profile and show no variation with temperature. The radiometric surface temperatures show that foam is cooler and more variable than a foam-free surface, and bubble-resolving thermal images show that radiometrically transparent bubble caps and burst bubbles reveal warm foam below the cool surface layer, contributing to the enhanced variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A VLA Survey of Late-time Radio Emission from Superluminous Supernovae and the Host Galaxies

We present the results of 3 GHz radio continuum observations of 23 superluminous supernovae (SLSNe) and their host galaxies by using the Karl G. Jansky Very Large Array conducted 5–21 yr after the explosions. The sample consists of 15 Type I and 8 Type II SLSNe at z < 0.3, providing one of the largest samples of SLSNe with late-time radio data. We detected radio emission from one SLSN (PTF10hgi) and five hosts with a significance of >5σ. No time variability is found in late-time radio light curves of the radio-detected sources in a timescale of years except for PTF10hgi, whose variability is reported in a separate study. Comparison of star formation rates (SFRs) derived from the 3 GHz flux densities with those derived from SED modeling based on UV–NIR data shows that four hosts have an excess of radio SFRs, suggesting obscured star formation. Upper limits for undetected hosts and stacked results show that the majority of the SLSN hosts do not have a significant obscured star formation. By using the 3 GHz upper limits, we constrain the parameters for afterglows arising from interaction between initially off-axis jets and circumstellar medium (CSM). Here, we found that the models with higher energies ( E iso ≳ several × 10 53 erg) and CSM densities ( n ≳ 0.01 cm –3 ) are excluded, but lower energies or CSM densities are not excluded with the current data. We also constrained the models of pulsar wind nebulae powered by a newly born magnetar for a subsample of SLSNe with model predictions in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The 16th Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

The Sloan Digital Sky Surveys (SDSS) have been observing the skies from Apache Point Observatory (APO) since 1998 (using the 2.5m Sloan Foundation Telescope, Gunn et al. 2006) and from Las Campanas Observatory (LCO) since 2017 (using the du Pont 2.5m Telescope). Representing the fourth phase of SDSS, SDSS-IV (Blanton et al. 2017) consists of three main surveys; the Extended Baryon Oscillation Spectroscopic Survey (eBOSS; Dawson et al. 2016), Mapping Nearby Galaxies at APO (MaNGA; Bundy et al. 2015), and the APO Galactic Evolution Experiment 2 (APOGEE-2; Majewski et al. 2017). Within eBOSS, SDSS-IV has also conducted two smaller programs: the SPectroscopic IDentification of ERosita Sources (SPIDERS; Clerc et al. 2016; Dwelly et al. 2017) and the Time Domain Spectroscopic Survey (TDSS; Morganson et al. 2015). These programs have investigated a broad range of cosmological scales, including cosmology with large-scale structure in eBOSS, the population of quasars and variable or X-ray-emitting stars with TDSS and SPIDERS; nearby galaxies in MaNGA; and the Milky Way and its stars in APOGEE-2.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Development of an Unmanned Mobile Current Turbine Platform: Preprint

The design and development of a prototype unmanned mobile floating platform, equipped with a custom low-flow marine current turbine for autonomously seeking and harnessing tidal/coastal currents is described. The platform is an unmanned surface vehicle in the form of a catamaran with two electric outboard motors and capabilities for autonomous navigation. An undershot water wheel, aided by a custom flow concentrator, has been selected as the basic design for the marine current turbine, which is mounted on the stern of the unmanned surface vehicle platform. The concept of operation is that the platform would navigate to a designated marine current resource, autonomously anchor at the location, align itself in the current and deploy the current turbine using a custom cable-lift deployment mechanism. As the turbine harnesses the local current, an onboard power-take-off device converts the harnessed mechanical energy to electricity which is stored in onboard batteries. Considerations of deployment in tidal and coastal waters required obtaining the necessary environmental permits for conducting in-water testing; developing required mitigation measures in protecting local wildlife and their habitats; and identifying potential in-water test sites and surveying them for their suitability in terms of the current resource, bottom type, water depth, and local boat traffic. The design and development of the turbine and the results of initial in-water testing are discussed. The potential for scaling up the system for extended capacity is presented.

ENGINEERING,TIDAL AND WAVE POWER↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Data-assisted combustion simulations with dynamic submodel assignment using random forests

This investigation outlines a data-assisted approach that employs random forest classifiers for local and dynamic submodel assignment in turbulent-combustion simulations. This method is demonstrated in simulations of a single-element GOX/GCH4 rocket combustor; a priori as well as a posteriori assessments are conducted to (i) evaluate the accuracy and adjustability of the classifier for targeting different quantities of interest (QoIs), and (ii) assess improvements, resulting from the data-assisted combustion model assignment, in predicting target QoIs during simulation runtime. Results from the a priori study show that random forests, trained with local flow properties as input variables and combustion model errors as training labels, assign three different combustion models – finite-rate chemistry (FRC), flamelet progress variable (FPV) model, and inert mixing (IM) – with reasonable classification performance even when targeting multiple QoIs. Applications in a posteriori studies demonstrate improved predictions from data-assisted simulations, in temperature and CO mass fraction, when compared with monolithic FPV calculations. An additional a posteriori data-assisted simulation of a modified configuration demonstrates that the present approach can be successfully applied to different configurations, as long as thermophysical behavior can be represented by the training data. Furthermore, these results demonstrate that this data-driven framework holds promise for dynamic combustion submodel assignments in reacting flow simulations.

42 ENGINEERING↗
DOE OSTI2020

Mechanistic Insights into C–H Borylation of Arenes with Organoiridium Catalysts Embedded in a Microporous Metal–Organic Framework

Organometallic iridium catalysts can be used in conjunction with bispinacolatodiboron (B 2 Pin 2 ) to effect the borylation of a variety of substrates such as arenes, alkanes, heteroarenes, and oxygenates. Recently, efforts have also focused on integrating these catalysts into porous supports, such as metal–organic frameworks (MOFs). While the mechanism of homogeneous borylation systems has been thoroughly investigated experimentally and computationally, analogous studies in MOF-supported iridium catalysts have not been conducted. In this paper, we report the mechanistic investigation of a phenanthroline-iridium catalyst immobilized in the organic linker of Universitetet i Oslo (UiO)-67 (Zr 6 O 4 (OH) 4 (BPDC) 4 (PhenDC) 2 , BPDC = biphenyl-4,4'-dicarboxylate, PhenDC = 1,10-phenanthroline-4,4'-dicarboxylate). By using benzene as a model substrate, variable time normalization analysis (VTNA) of the kinetic data suggested a rate law consistent with zero-order in B 2 Pin 2 , and first-order in arene. A primary kinetic isotope effect (KIE) in the time course of benzene-d 6 borylation also provided complementary information about the role of the arene in the rate-determining step of the reaction. Characterization by techniques such as X-ray absorption spectroscopy (XAS) confirmed the presence of Ir(III), while pair distribution function (PDF) analysis suggested structures containing an Ir–Cl bond, further substantiated by X-ray photoelectron spectroscopy (XPS). Analysis of postcatalysis materials by inductively coupled plasma–optical emission spectroscopy (ICP-OES) revealed low boron accumulation, which may indicate an absence of boron in the resting state of the catalyst. Finally, in comparing borylation of benzene and toluene, a slight selectivity for benzene is observed, which is similar to the analogous homogeneous reaction, indicating the influence of substrate sterics on reactivity.

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

Drag Coefficient Analysis of Arbitrary Fragment Geometry with Emphasis on Meshing Sensitivity Studies.

This report characterizes Mach number effects in the hypersonic range from Mach 3 to Mach 7 for an arbitrary, computationally generated fragment. Force and moment coefficients in the x, y, z direction were successfully obtained and resulting trends were compared with theoretical expectations. Three orientations were tested in the form of unit vectors: [1,0,0], [0,1,0], and [0,0,1]; these represent the i, j, k principle axes. The final results ultimately showcased a trend very similar to that shown in previous literature — drag coefficients for a given orientation decreased with increasing Mach number by a very small amount in the hypersonic range. It was therefore concluded that the aerodynamic quantities used to obtain the fragment trajectory at low hypersonic Mach numbers (such as Mach 3) can be used to characterize the aerodynamic qualities at higher hypersonic Mach numbers (Mach 5, 7) with a reasonably small margin. A sensitivity study was also conducted which ultimately defined drag effects as a result of changing grid spacing and normal extrusion parameters for the unstructured tetrahedral mesh used in the simulation. Further work can be done in the form of experimental validation, specifically with regards to wind tunnels and ballistic range testing. Error can be reduced by testing more mesh variabilities and capturing a larger amount of fragment orientations. The results for force and moment characteristics dependent on meshing parameters and high Mach numbers are satisfactory and consistent with expected trends.

42 ENGINEERING↗
DOE OSTI2020

Monitoring Avian Productivity and Survivorship at the Sandia Wetlands at Los Alamos National Laboratory (2014-2019)

Biologists at Los Alamos National Laboratory completed six years (2014–2019) of demographic monitoring of passerines (songbirds) during the breeding season in Los Alamos, New Mexico. We captured songbirds at a mist netting station located in the Sandia Wetlands in Sandia Canyon (Technical Areas 60 and 61). We identified, measured, and then banded captured birds with federal migratory bird bands. Banding operations took place between May and August of each year, with 10 mist netting sessions annually. We conducted this project as part of the implementation of the Biological Resources Management Plan. The project complies with the 2013 Memorandum of Understanding between the United States Fish and Wildlife Service (USFWS) and the United States Department of Energy regarding the implementation of Executive Order 13186. The overall number of birds and species captured was variable, and we saw no decreases over the course of six years in bird species diversity, abundance, or the percentage of breeding birds. Our results indicate that the Sandia Wetlands support numerous species of breeding birds, including species of conservation concern. We did not see a declining trend for the top 10 most-captured species or for any of the species of conservation concern over time; however, we need more years of data to make robust conclusions about population trends through time. Our results indicated that adult recapture rates at the Sandia Wetlands are at least twice as low when compared with estimated adult survival in migratory passerines from peer-reviewed literature. We also did not recapture any birds that were banded as juveniles in subsequent years throughout the entire six-year period. Low recapture rates for adults and the non-existent recapture rates for juveniles should be explored further. To accurately assess individual species survival rates for adults and juveniles using mark-recapture models, we need a longer dataset, which would allow for a more robust analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Do Piperonyl Butoxide Long-Lasting Insecticide Treated Nets Provide Additional Protection Against Malaria Infections Compared with Conventional Nets in an Operational Setting in Western Kenya?

Malaria control in sub-Saharan Africa has stagnated despite widespread adoption of control measures such as long-lasting insecticidal nets (LLINs). Progress has stalled, in part, because of pyrethroid insecticide resistance, driving the need for retooling to increase the effectiveness of bed nets. Consequently, LLINs have been treated with the chemical synergist piperonyl butoxide (PBO). Piperonyl butoxide LLINs have been shown to be efficacious in controlled settings; however, their effectiveness in real-world settings warrants investigation. In Bungoma County, Western Kenya, a cohort of 768 participants was followed from June 2017 to December 2023 via active and passive surveillance. Household visits were conducted monthly, during which LLIN use for nets distributed in 2017 and 2021 was recorded, and symptomatic malaria cases were identified using rapid diagnostic tests (RDTs). The comparative effectiveness of PBO versus conventional LLINs was assessed in terms of malaria infections. A multilevel logistic regression model was fit with monthly RDT results as the dependent variable. The study results indicate that PBO LLINs provide greater protection against malaria at the individual level than conventional LLINs (odds ratio: 0.70; 95% CI: 0.47–1.03), although the findings were not statistically significant. The added protection against malaria infections provided by PBO LLINs compared with conventional LLINs observed in the current study aligns with findings from most previous studies, although this finding was not statistically significant. In areas with documented pyrethroid resistance, the use of LLINs with an added synergist, such as PBO, can provide additional protection against malaria infections (compared with pyrethroid-only LLINs) and should be considered for scaled-up scenarios despite the additional cost.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

ADDITIVELY MANUFACTURED SURFACE HEAT TRANSFER ENHANCEMENTS FOR THE TRANSFORMATIONAL CHALLENGE REACTOR

The Transformational Challenge Reactor (TCR) is a high-temperature gas-cooled reactor design that uses additively manufactured fuel elements. TCR fuel elements have walls made of silicon carbide and are filled with tristructural-isotropic fuel particles. These fuel elements can have radically different shapes and integrated features than existing designs due to the reduced cost for complex structures in binder jet additive manufacturing. As such, binder jet additive manufacturing enables wall surface features to be embedded that can deliver superior heat transfer performance than smooth wall designs. In this work, the authors conducted a computational fluid dynamics study to evaluate selected wall features integrated into TCR fuel elements. Results show that surface features can outperform smooth wall designs; however, there are unique challenges for gas-cooled reactor core designs that have not been fully explored by previous research. For example, the rough surface finish and process variability of ceramics additive manufacturing make it challenging to predict surface roughness effects before fabrication. Additionally, the small hydraulic diameters of coolant channels in reactor cores make it difficult to engineer surface features that do not significantly increase the pressure drop. Engineers must carefully size surface features for heat transfer enhancement in additive fuel elements to operate above the base material's surface roughness effects and below the coolant channel size.

Weinmeister, Justin↗
DOE OSTI2024

Estimating irrigation water use from remotely sensed evapotranspiration data: Accuracy and uncertainties at field, water right, and regional scales

Irrigated agriculture is the dominant user of water globally, but most water withdrawals are not monitored or reported. As a result, it is largely unknown when, where, and how much water is used for irrigation. Here, we evaluated the ability of remotely sensed evapotranspiration (ET) data, integrated with other datasets, to calculate irrigation water withdrawals and applications in an intensively irrigated portion of the United States. We compared irrigation calculations based on an ensemble of satellite-driven ET models from OpenET with reported groundwater withdrawals from hundreds of farmer irrigation application records and a statewide flowmeter database at three spatial scales (field, water right group, and management area). At the field scale, we found that ET-based calculations of irrigation agreed best with reported irrigation when the OpenET ensemble mean was aggregated to the growing season timescale (bias = 1.6–4.9%, R 2 = 0.53–0.74), and agreement between calculated and reported irrigation was better for multi-year averages than for individual years. At the water right group scale, linking pumping wells to specific irrigated fields was the primary source of uncertainty. At the management area scale, calculated irrigation exhibited similar temporal patterns as flowmeter data but tended to be positively biased with more interannual variability. Disagreement between calculated and reported irrigation was strongly correlated with annual precipitation, and calculated and reported irrigation agreed more closely after statistically adjusting for annual precipitation. The selection of an ET model was also an important consideration, as variability across ET models was larger than the potential impacts of conservation measures employed in the region. From these results, we suggest key practices for working with ET-based irrigation data that include accurately accounting for changes in soil moisture, deep percolation, and runoff; careful verification of irrigated area and well-field linkages; and conducting application-specific evaluations of uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Testing and signal identification for two-sample high-dimensional covariances via multi-level thresholding

The paper considers testing and signal identification for covariance matrices from two populations of marginally sub-Gaussian distributed. A multi-level thresholding procedure is proposed for testing the equality of two high-dimensional covariance matrices, which is designed to detect sparse and faint differences between the covariances. A novel -statistic composition is developed to establish the asymptotic distribution of the thresholding statistics in conjunction with the matrix blocking and the coupling techniques. It is shown that the proposed test is more powerful than the existing tests in detecting sparse and weak signals in covariances. Multiple testing procedures are constructed to discover different covariances and the sub-groups of variables with different covariance structures between the two populations. The proposed procedures are based on the multi-level thresholding test, which are able to control the false discovery proportion () with high power. Simulation experiments and a case study on the returns of the S&P 500 stocks before and after the COVID-19 pandemic are conducted to demonstrate and compare the utilities of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Geostatistical interpolation of streambed hydrologic attributes with addition of left censored data and anisotropy

Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.

58 GEOSCIENCES↗