Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “measurement errors”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CeNTREX: a new search for time-reversal symmetry violation in the 205 Tl nucleus

The Cold molecule Nuclear Time-Reversal EXperiment (CeNTREX) is a new effort aiming for a significant increase in sensitivity over the best present upper bounds on the strength of hadronic time reversal (T) violating fundamental interactions. The experimental signature will be shifts in nuclear magnetic resonance frequencies of 205 Tl in electrically-polarized thallium fluoride (TlF) molecules. Here we describe the motivation for studying these T-violating interactions and for using TlF to do so. To achieve higher sensitivity than earlier searches for T-violation in TlF, CeNTREX uses a cryogenic molecular beam source, optical state preparation and detection, and modern methods of coherent quantum state manipulation. Details of the measurement scheme and the current state of the apparatus are presented, with quantitative measurements of the TlF beam. Lastly, the estimated sensitivity and methods to control systematic errors are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Analysis of Three Types of Collocated Disdrometer Measurements at the ARM Southern Great Plains Observatory

To better provide U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s disdrometer deployment strategy and support ARM precipitation-related projects, this study analyzes 18 months of drop size distribution (DSD) observations from six collocated disdrometers deployed at the ARM Southern Great Plains (SGP) site. Emphasis is placed on quantifying the uncertainties related to DSD and rainfall properties using different types of disdrometers and different pairing concepts. The instrument sampling errors are also discussed, which improves our understanding of instrument-specific impacts on rainfall measurements. The key findings are as follows: 1) All the disdrometers show consistent behaviors overall regarding accumulated precipitation, mean rainfall rate, DSD parameters, and radar reflectivity values; and 2) Strong agreement is shown between three disdrometer types in terms of their DSDs for mid-size drop range (D = 1 mm-4 mm). The two-dimensional video disdrometer performance meets expectations as the most accurate unit among all. The paired disdrometer shows a reduced statistical sampling error in terms of drop counts and other DSD properties. A concept of deploying two or more collocated disdrometers side by side is recommended for future ARM field campaigns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Balancing Impedance and Controllability in Response Reconstruction with $\mathrm{TS-IMMAT}$

In a prior work the authors proposed a variant on the Impedance Matched Multi-Axis Test (IMMAT) in which a fixture is defined, called the Transmission Simulator (TS), and the desired environment is matched at a set of sensors on the TS. If the motion of the TS is matched then the response of the rest of the component will also match, provided that the attached component has the same dynamics as it did when the environment was measured. Hence, one would like the TS to be flexible so that it reproduces the boundary conditions that the component of interest experiences during flight, but the more flexible the TS, the more shakers might be needed to control its response. This work presents a derivation that gives expressions for these two potential error sources in TS-IMMAT. Then, various case studies are presented, both on simulated and real hardware, to understand the importance of each error term in practical testing. The theory explains the phenomena that were observed when using measurements from a component that flew on a sounding rocket. The environmental response was measured and then various fixtures were attached, each comprising more of the next assembly, or the hardware to which the component was attached in flight. MIMO testing was repeated with each fixture and the results were compared to seek to understand the role of the impedance match in this type of testing. Here the results show that the number of modes that are active in the transmission simulator is also very important, and so the best solution balances these two considerations. An improved method of simulating the MIMO test is then proposed, so simulations can be used to predict what fixture, or transmission simulator, will give the best results in a TS-IMMAT test.

42 ENGINEERING↗
DOE OSTI2021

Laboratory and field studies of ice-nucleating particles from open-lot livestock facilities in Texas

Abstract. In this work, an abundance of ice-nucleating particles (INPs) from livestock facilities was studied through laboratory measurements from cloud-simulation chamber experiments and field investigation in the Texas Panhandle. Surface materials from two livestock facilities, one in the Texas Panhandle and another from McGregor, Texas, were selected as dust proxies for laboratory analyses. These two samples possessed different chemical and biological properties. A combination of aerosol interaction and dynamics in the atmosphere (AIDA) measurements and offline ice spectrometry was used to assess the immersion freezing mode ice nucleation ability and efficiency of these proxy samples at temperatures above −29 ∘C. A dynamic filter processing chamber was also used to complement the freezing efficiencies of submicron and supermicron particles collected from the AIDA chamber. For the field survey, periodic ambient particle sampling took place at four commercial livestock facilities from July 2017 to July 2019. INP concentrations of collected particles were measured using an offline freezing test system, and the data were acquired for temperatures between −5 and −25 ∘C. Our AIDA laboratory results showed that the freezing spectra of two livestock dust proxies exhibited higher freezing efficiency than previously studied soil dust samples at temperatures below −25 ∘C. Despite their differences in composition, the freezing efficiencies of both proxy livestock dust samples were comparable to each other. Our dynamic filter processing chamber results showed on average approximately 50 % supermicron size dominance in the INPs of both dust proxies. Thus, our laboratory findings suggest the importance of particle size in immersion freezing for these samples and that the size might be a more important factor for immersion freezing of livestock dust than the composition. From a 3-year field survey, we measured a high concentration of ambient INPs of 1171.6 ± 691.6 L−1 (average ± standard error) at −25 ∘C for aerosol particles collected at the downwind edges of livestock facilities. An obvious seasonal variation in INP concentration, peaking in summer, was observed, with the maximum at the same temperature exceeding 10 000 L−1 on 23 July 2018. The observed high INP concentrations suggest that a livestock facility is a substantial source of INPs. The INP concentration values from our field survey showed a strong correlation with measured particulate matter mass concentration, which supports the importance of size in ice nucleation of particles from livestock facilities.

58 GEOSCIENCES↗
DOE OSTI2025

Impact of dynamics, entanglement and Markovian noise on the fidelity of few-qubit digital quantum simulation

Quantum algorithms have been proposed to accelerate the simulation of the chaotic dynamical systems that are ubiquitous in the physics of plasmas. Quantum computers without error correction might even use noise to their advantage to calculate the Lyapunov exponent by measuring the Loschmidt echo fidelity decay rate. For the first time, digital Hamiltonian simulations of the quantum sawtooth map, performed on the IBM-Q quantum hardware platform, show that the fidelity decay rate of a digital quantum simulation increases during the transition from dynamical localization to chaotic diffusion in the map. The observed error per CNOT gate increases by $1.5{\times }$ as the dynamics varies from localized to diffusive, while only changing the phases of virtual RZ gates and keeping the overall gate count constant. A gate-based Lindblad noise model that captures the effective change in relaxation and dephasing errors during gate operation qualitatively explains the effect of dynamics on fidelity as being due to the localization and entanglement of the states created. Specifically, highly delocalized states that are entangled with random phases show an increased sensitivity to dephasing and, on average, a similar sensitivity to relaxation as localized states. In contrast, delocalized unentangled states show an increased sensitivity to dephasing but a lower sensitivity to relaxation. This gate-based Lindblad model is shown to be a useful benchmarking tool by estimating the effective Lindblad coherence times during CNOT gates and finding a consistent $2\unicode{x2013}3{\times }$ shorter $T_2$ time than reported for idle qubits. Thus, the interplay of the dynamics of a simulation with the noise processes that are active can strongly influence the overall fidelity decay rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparing Outdoor to Indoor Performance for Bifacial Modules Affected by Polarization-Type Potential-Induced Degradation

Bifacial photovoltaic (PV) modules have the advantage of using light reflected off of the ground to contribute to power production. Predicting the energy gain is challenging and requires complex models to do so accurately. Often, module degradation over time is neglected in models for the sake of simplicity or is underestimated. Comparing outdoor and indoor current–voltage (I–V) performance for bifacial modules is more challenging than for monofacial modules, as there are additional variables to consider such as rear albedo non-uniformity, cell mismatch, and their effects on temperature. This challenge is compounded when heterogeneous degradation modes occur, such as polarization-type potential-induced degradation (PID-p). To examine the effects of PID-p on I–V predictions using an empirical data-driven approach, 16 bifacial PERC modules are installed outdoors on racks with different albedo conditions. A subset is exposed to high-voltage biases of −1500 V or +1500 V. Outdoor data are traced at irradiance ranges of 150–250 W/m 2 , 500–600 W/m 2 , and 900–1000 W/m 2 . These curves are corrected using control module temperature, wire resistivity, and module resistance measured indoors. We examine several methods to transform indoor I–V curves to accurately, and more simply than existing methods, approximate outdoor performance for bifacial modules without and with varying levels of PID-p degradation. This way, bifacial performance modeling can be more accessible and informed by fielded, degraded modules. Distributions of percent errors between indoor and outdoor performance parameters and Mean Absolute Percent Errors (MAPEs) are used to assess method quality. Results including low-irradiance data (150–250 W/m 2 ) are discussed but are filtered for quantifying method quality as these data introduce substantial errors. The method with the most optimal tradeoff between low MAPE and analysis simplicity involves measuring the front side of a module indoors at an irradiance equal to plane-of-array irradiance plus the product of module bifaciality and albedo irradiance. This method gives MAPE values of 1–6.5% for non-degraded and 1.6–5.9% for PID-p degraded module performance.

14 SOLAR ENERGY↗
DOE OSTI2023

A nanocryotron memory and logic family

The development of superconducting electronics based on nanocryotrons has been limited so far to few device circuits, in part due to the lack of standard and robust logic cells. Here, we introduce and experimentally demonstrate designs for a set of nanocryotron-based building blocks that can be configured and combined to implement memory and logic functions. The devices were fabricated by patterning a single superconducting layer of niobium nitride and measured in liquid helium on a wide range of operating points. The tests show 10 - 4 bit error rates with above ± 20 % margins up to 50 MHz and the possibility of operating under the effect of an out-of-plane 36 mT magnetic field, with ± 30 % margins at 10 MHz. Additionally, we designed and measured an equivalent delay-flip-flop made of two memory cells to show the possibility of combining multiple building blocks to make larger circuits. These blocks may constitute a solid foundation for the development of nanocryotron logic circuits and finite-state machines with potential applications in the integrated processing and control of superconducting nanowire single-photon detectors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Optimal Power Flow With State Estimation in the Loop for Distribution Networks

Here in this article, we propose a framework for running optimal control-estimation synthesis in distribution networks. Our approach combines a primal-dual gradient-based optimal power flow solver with a state estimation feedback loop based on a limited set of sensors for system monitoring, instead of assuming exact knowledge of all states. The estimation algorithm reduces uncertainty on unmeasured grid states based on certain online state measurements and noisy "pseudomeasurements." We analyze the convergence of the proposed algorithm and quantify the statistical estimation errors based on a weighted least-squares estimator. The numerical results on a 4521-node network demonstrate that this approach can scale to extremely large networks and provide robustness to both large pseudomeasurement variability and inherent sensor measurement noise.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Beam dynamics corrections in the measurement of the anomalous precession frequency at the Muon $g-2$ experiment at Fermilab

The Muon $g-2$ experiment at Fermilab (E989) aims to measure the anomalous magnetic moment of the muon with an accuracy of 140 ppb (parts per billions). This accuracy, obtained by adding in quadrature a statistical and a systematic contribution of comparable value (100 ppb), will allow to reduce the experimental uncertainty (from the previous E821 experiment at BNL \cite{bnl}) of a factor of 4, and represents one of the most precise tests of the Standard Model (SM) theory of elementary particles. The first result on the Run-1 dataset \cite{prl} was released on April 7, 2021, showing a very good agreement with the previous result from BNL experiment, with a slightly better uncertainty. The corresponding experimental average increases the significance of the discrepancy between the measured and Standard Model prediction of 4.2$\sigma$ \cite{white_paper}.\\ \noindent The measured quantity is the muon magnetic anomaly a$_\mu$=$\frac{g_\mu-2}{2}$ where g$_\mu$ is the gyromagnetic factor of the muon. Dirac's equation predicts g$_\mu$ = 2, while radiative corrections, dominated by the QED contribution due to an exchange of a virtual photon, causes a per-mille correction on this quantity. By including all the SM contributions, a$_\mu$ is known at 370 ppb. The E989 experiment measures $a_\mu$ injecting positive muons with momentum of 3.1 GeV/$c$ polarized longitudinally in bunches (called $fills$) with an average rate of 12 Hz, in a storage ring of 14 meters diameter. Due to the parity violation in the weak muon decay, high energy positrons produced are emitted preferably in the muon's spin direction. By counting the number of positrons with energy greater then 1.7 GeV in function of the time, the frequency precession of the muon spin is measured, that together with the measurement of the magnetic field, allows to extract $a_\mu$. The positrons are detected with 24 electromagnetic calorimeters, that measure the energy and the arrival time of the positrons, each made of 54 crystals of lead fluoride (PbF$_2$) read by silicon photomultipliers (SiPM). Together with calorimeters, two tracking detectors are used to make non-destructive measurements of the muon beam characteristic by reconstructing the muon decay position extrapolating backward the decay positrons. The knowledge of the beam motion inside the ring plays a fundamental role in the analysis of $a_{\mu}$, where the measured anomalous precession frequency $\omega_a$ must be corrected for four main beam dynamics effects. A first correction is associated with the presence of an electric field responsible of the vertical focusing on the storage ring, where vertical direction is orthogonal to the orbit and horizontal direction is along the storage ring radius. Due to the oscillations in the vertical direction, the so called vertical betatron oscillations, a second correction is necessary to account for an average angle associated with the muon motion off the ideal orbit. A third correction is caused by lost muons in the ring which have a different spin phase at the injection respect to the decay ones. Finally, due to the correlation between the vertical and horizontal beam motion of the muons and the acceptance of the calorimeter, a correction (called ``phase-acceptance") arises. In Run1, due to the presence of two damaged resistors in one electrostatic quadrupole, this correction had a prominent role. The replacement of the damaged resistors before Run2 decreased this effect by one order of magnitude, and a further reduction in Run-3 was provided by an improved orbit.\\ \noindent The work of this Thesis focuses on the beam dynamics corrections on $\omega_a$. Due to the prominent role in Run-1 a special attention has been put to the phase-acceptance correction\footnote{The effect of this correction on the E821 BNL measurement of the $g-2$ has been evaluated to be within the quoted systematic error.}. Chapter \ref{ch:anomalous_magnetic_moment} introduces the anomalous magnetic moment of the muon. Chapter \ref{ch:early_experiment} describes the history of the Muon $g-2$ experiments. Chapter \ref{ch:theory} discusses the Standard Model prediction and possible new physics scenario. Chapter \ref{ch:muon_g2_experiment} describes the E989 experiment, whose experimental technique and the apparatus are discussed focusing on the improvements needed to reach the final goal on $a_\mu$ measurement. The original work of the Thesis is discussed in the last 5 chapters. Chapter \ref{ch:experiment_technique} presents the analysis technique to extract $\omega_a$, Chapter \ref{ch:beam_dynamics} describes the beam dynamics in the E989 experiment, Chapter \ref{ch:beam_dynamics_correction} discusses the beam dynamics correction to $\omega_a$, and Chapter \ref{sec:pa_corr} focuses on the phase acceptance correction. Finally Chapter \ref{ch:run23_analysis_improvements} contains the Run-2/3 improvements which are expected to allow for an increased precision ons the beam dynamics co...

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Analysis of limited coverage effects on areal density measurements in inertial confinement fusion implosions

Accurate diagnosis of areal density ($\textit{ρR}$) is critical for the inference of performance metrics in inertial confinement fusion implosions. One potential source of error in this diagnosis is the existence of low mode perturbations in the imploding target, which lead to asymmetries in the inference of the $\textit{ρR}$ from different lines of sight. Here, the error accrued as a result of limited coverage of the sphere due to a finite number of detectors is quantified, and the development of a forward scatter measurement from the OMEGA neutron time-of-flight detectors is motivated. In this work, a method by which the 1D-equivalent 4 π-averaged $\langleρR\rangle$ can be reconstructed, if accurate mode information can be diagnosed by other means, is validated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Validation of the DESI DR2 measurements of baryon acoustic oscillations from galaxies and quasars

The Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) galaxy and quasar clustering data represents a significant expansion of data from Data Release 1 (DR1), providing improved statistical precision in baryon acoustic oscillation (BAO) constraints across multiple tracers, including bright galaxies, luminous red galaxies, emission line galaxies, and quasars. In this paper, we validate the BAO analysis of DR2. We present the results of robustness tests on the blinded DR2 data and, after unblinding, consistency checks on the unblinded DR2 data. All results are compared with those obtained from a suite of mock catalogs that replicate the selection and clustering properties of the DR2 sample. We confirm the consistency of DR2 BAO measurements with DR1 while achieving a reduction in statistical uncertainties due to the increased survey volume and completeness. The combined BAO precision, including both statistical and systematic errors, improves from ∼0.52% in DR1 to 0.30% in DR2—a factor of 1.7 gain. We assess the impact of analysis choices, including different data vectors (correlation function vs power spectrum), modeling approaches and systematics treatments, and an assumption of the Gaussian likelihood, finding that our BAO constraints are stable across these variations and assumptions with a few minor refinements to the baseline setup of the DR1 BAO analysis. We summarize a series of pre-unblinding tests that confirmed the readiness of our analysis pipeline, the final systematic errors, and the DR2 BAO analysis baseline. The successful completion of these tests led to the unblinding of the DR2 BAO measurements, ultimately leading to the DESI DR2 cosmological analysis, with their implications for the expansion history of the Universe and the nature of dark energy presented in the DESI key paper (companion paper).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Spectroscopic Measurements and Models of Energy Deposition in the Substrate of Quantum Circuits by Natural Ionizing Radiation

Naturally occurring background radiation is a potential source of correlated decoherence events in superconducting qubits that will challenge error-correction schemes. In order to characterize the radiation environment in an unshielded laboratory representative of superconducting qubits’ environments, we performed broadband, spectroscopic measurements of background radiation events inside a millikelvin refrigerator. The spectrometer was designed to mimic the size and composition of a quantum circuit. Specifically, we measured the background radiation spectra in silicon substrates of two thicknesses, 500 and 1500 µm, and one area, 25 mm 2 . The observed spectra span energies from a few kilo-electron-volts up to nearly 10 MeV, are nearly featureless, and decrease in intensity by a factor of 40 000 between 100 keV and 3 MeV for the 500-µm substrate. We integrate the spectra to obtain the average event rates and deposited power levels. These quantities correspond to a rate of 0.023 events per second and a power of 4.9 keV s -1 , when counting events that deposit at least 40 keV for the 500-µm-thick substrate. We find that the cryogenic measurements are in good agreement with predictions based on simple measurements of the terrestrial gamma-ray flux outside the refrigerator, published models of cosmic-ray fluxes, a crude model of the cryostat, and radiation-transport simulations. This model requires no free parameters to predict the background radiation spectra in the silicon substrates. The agreement between measurements and predictions demonstrates that the model we present can be used to assess the relative contributions of terrestrial and cosmic-ray sources to background radiation interactions in silicon substrates of varying thickness. These spectroscopic measurements are performed with a novel combination of superconducting microresonators located on micromachined silicon islands that define the interaction volume with background radiation. The resonators transduce deposited energy to a readily detectable electrical signal. Microresonator readout closely resembles dispersive superconducting qubit readout, so similar devices—with or without micromachined islands—are suitable for integration with superconducting quantum circuits as detectors for background radiation events. For our specific laboratory conditions, we find that gamma-ray emissions from radioisotopes are responsible for the majority of events that deposit E < 1 ⁢Me⁢V. We present results demonstrating that the background radiation spectrum contains relevant contributions from cosmic-ray particles other than muons, particularly a tail of multi-mega-electron-volt events due to protons and neutrons. These observations suggest several paths to reducing the impact of background radiation on quantum circuits, supported by an empirically validated model for generating reliable predictions of radiation interactions with silicon substrates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

New Lightning-Derived Vertical Total Electron Content Data Provide Unique Global Ionospheric Measurements

A newly released, novel ionospheric data set of global gridded vertical total electron content (VTEC) is introduced in this paper. This VTEC data set, provided by Los Alamos National Laboratory (LANL), is derived from very high frequency (VHF; defined as 30–300 MHz) broadband radio-frequency (RF) measurements of lightning made by U.S. Department of Defense sensing systems on board Global Positioning System (GPS) satellites. This paper presents the new data set (LANL VTEC), discusses the errors inherent in VHF TEC estimation due to ionospheric dispersion, and compares the LANL VTEC to two community standard VTEC gridded products: Jet Propulsion Laboratory's Global Ionospheric Model (JPL GIM) and the CEDAR community's Open Madrigal VTEC gridded measurements of L-band GNSS (global navigation satellite systems) TEC. We find that the LANL VTEC data have an offset of 3 TECU from CEDAR Madrigal GNSS VTEC and a full-width-half-maximum (FWHM) of 6 TECU. In comparison, the offset between LANL VTEC and the JPL GIM model is –3 TECU, but with a FWHM of 5 TECU. We also compare to Jason-3 VTEC measurements over the ocean, finding an offset of less than 0.5 TECU and an FWHM of <5 TECU. Because this technique uses a completely different methodology to determine TEC, the sources of errors are distinct from the typical ground-based GNSS L-band (GHz) TEC measurements. Also, because it is derived from RF lightning signals, this data set provides measurements in regions that are not well covered by ground-based GPS measurements, such as over oceans and over central Africa.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Correction for Neutral Pressure-Driven Signal in Radiated Power Measurements on Proto-MPEX

A pinhole camera using four-channel resistive bolometers was designed and implemented on the Prototype Material Plasma Exposure eXperiment (Proto-MPEX) to measure the plasma radiated power. Quantification of the radiated power is required for comprehensive power balance assessments but faces unique challenges on linear plasma devices due to weak absolute signal, compared to typical use in toroidal fusion plasmas, and high neutral pressure changes. Around 1 mTorr, the pressure-driven signal contribution was shown to overwhelm the radiated power signal in Proto-MPEX, found to be ~3 μW . Since the pressure-driven signals could not be predicted from first principle arguments, the signals were corrected using empirical models that take two experimental pressure measurements and the raw bolometer signal to reveal the radiated power signals. Four “gas-only” Proto-MPEX shots were taken, without the radio frequency (rf) plasma discharge, to train, optimize, and validate empirical models. The data included two pressure measurements at different locations along with the bolometer signal. Multiple models were tested; partial least-squares regression with a cross correlation method to account for the time offset in the pressure measurements provided the best combination of accuracy and stability. The model was validated to have a modeling error of 0.0036 V, which was lower than 10% of the expected radiated power signal. The model was tested on a different “gas-only” shot with a complete time-series data. For plasma operation shots, the corrected bolometer measurements were compared with an absolute extreme ultraviolet photodiode measurement and similarities and differences are discussed.

Lee, Seungsup↗
DOE OSTI2021

Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis

The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey Year 3 BAO sample over its full photometric redshift range 0.6 < zphoto < 1.1. The mocks are based upon 488 ICE-COLA fast N-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution – halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function w(θ), angular power spectrum Cℓ, and projected clustering ξp(r⊥) to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.Key words: catalogs / large-scale structure of Universe / galaxies: distances and redshifts / Galaxy: halo / methods: numerical

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On negative induced polarization in frequency domain measurements

SUMMARY Induced polarization (IP) has been widely used to non-invasively characterize electrical conduction and polarization in the subsurface resulting from an applied electric field. Earth materials exhibit a lossy capacitance defined by a negative intrinsic phase in frequency-domain IP (FDIP) or a positive intrinsic chargeability in time-domain IP (TDIP). However, error-free positive apparent phase or negative apparent chargeability (i.e. negative IP effects) can occur in IP measurements over heterogeneous media. While negative IP effects in TDIP data sets have been discussed, no studies have addressed this topic in detail for FDIP measurements. We describe theory and numerical modelling to explain the origin of negative IP effects in FDIP measurements. A positive apparent phase may occur when a relatively high polarizability feature falls into negative sensitivity zones of complex resistivity measurements. The polarity of the apparent phase is determined by the distribution of subsurface intrinsic phase and resistivity, with the resistivity impacting the apparent phase polarity via its control on the sensitivity distribution. A physical explanation for the occurrence of positive apparent phase data is provided by an electric circuit model representing a four-electrode measurement. We also show that the apparent phase polarity will be frequency dependent when resistivity changes significantly with frequency (i.e. in the presence of significant IP effects). Consequently, negative IP effects manifest themselves in the shape of apparent phase spectra recorded with multifrequency (spectral IP) data sets. Our results imply that positive apparent phase measurements should be anticipated and should be retained during inversion and interpretation of single frequency and spectral IP data sets.

Geochemistry & Geophysics↗
DOE OSTI2023

Warming of the lower Columbia River , 1853 to 2018

Abstract Water temperature is a critical ecological indicator; however, few studies have statistically modeled century‐scale trends in riverine or estuarine water temperature, or their cause. Here, we recover, digitize, and analyze archival temperature measurements from the 1850s onward to investigate how and why water temperatures in the lower Columbia River are changing. To infill data gaps and explore changes, we develop regression models of daily historical Columbia River water temperature using time‐lagged river flow and air temperature as the independent variables. Models were developed for three time periods (mid‐19th, mid‐20th, and early 21st century), using archival and modern measurements (1854–1876; 1938–present). Daily and monthly averaged root‐mean‐square errors overall are 0.89°C and 0.77°C, respectively for the 1938–2018 period. Results suggest that annual averaged water temperature increased by 2.2°C ± 0.2°C since the 1850s, a rate of 1.3°C ± 0.1°C/century. Increased water temperatures are seasonally dependent. An increase of approximately 2.0°C ± 0.2°C/century occurs in the July–Dec time‐frame, while springtime trends are statistically insignificant. Rising temperatures change the probability of exceeding ecologically important thresholds; since the 1850s, the number of days with water temperatures over 20°C increased from ~5 to 60 per year, while the number below 2°C decreased from ~10 to 0 days/per year. Overall, the modern system is warmer, but exhibits less temperature variability. The reservoir system reduces sensitivity to short‐term atmospheric forcing. Statistical experiments within our modeling framework suggest that increased water temperature is driven by warming air temperatures (~29%), altered river flow (~14%), and water resources management (~57%).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An augmented Lagrangian filter method

Here, we introduce a filter mechanism to enforce convergence for augmented Lagrangian methods for nonlinear programming. In contrast to traditional augmented Lagrangian methods, our approach does not require the use of forcing sequences that drive the first-order error to zero. Instead, we employ a filter to drive the optimality measures to zero. Our algorithm is flexible in the sense that it allows for equality-constrained quadratic programming steps to accelerate local convergence. We also include a feasibility restoration phase that allows fast detection of infeasible problems. We provide a convergence proof that shows that our algorithm converges to first-order stationary points. We provide preliminary numerical results that demonstrate the effectiveness of our proposed method.

97 MATHEMATICS AND COMPUTING↗