Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reflection coefficient”

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 181 records · Page 10

Analysis of First-Order Gratings in Silicon Photonic Waveguides

A simple thin film effective index analysis for first order gratings in Si photonic waveguides is shown to provide highly accurate results for reflected and transmitted power spectrums as long as the waveguide remains single mode and non-radiating. A cover layer can be added to the grating region of a Si photonic waveguide to increase the strength of the grating, modify transition losses from the input waveguide to the grating waveguide region, and/or modify the width of the reflectivity spectrum. For a given grating period, the peak reflection and spectral width of the reflectivity decrease as the duty cycle is decreased or increased from ~50%. For both radiating and multimode structures, the coupling between all modes, power radiated towards the superstrate (upwards), power radiated downwards (substrate) and transmitted power analyzed by Floquet-Bloch, Eigenmode Expansion and Finite Difference Time Domain methods show excellent agreement. Coupling coefficients calculated using analytic formulas are shown to be accurate only for shallow grating depths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Performance Comparison of Low-Cost Near-Infrared (NIR) Spectrometers to a Conventional Laboratory Spectrometer for Rapid Biomass Compositional Analysis

The performance of a conventional laboratory near-infrared (NIR) spectrometer and two NIR spectrometer prototypes (a Texas Instruments NIRSCAN Nano evaluation model (EVM) and an InnoSpectra NIR-M-R2 spectrometer) are compared by collecting reflectance spectra of 270 well-characterized herbaceous biomass samples, building calibration models using the partial least squares (PLS-2) algorithm to predict five constituents of the samples from the reflectance spectra, and comparing the resulting model statistics. The prediction models developed using spectra from the Foss XDS spectrometer were slightly better than the prediction models developed using spectra from either the TI NIRSCAN Nano EVM and the InnoSpectra NIR-M-R2 as measured by the root mean square error (RMSECV) and the correlation coefficient (R 2 _cv) for “leave-one-out” cross-validation (CV). The models built from the two prototype units were not statistically significantly different from each other (p = 0.05). The Foss spectrometer has a larger wavelength range (400–2500 nm) compared with the two prototypes (900–1700 nm). When the spectra from the Foss XDS spectrometer were truncated so their wavelength range matched the wavelength range of the two prototype units, the resulting model was not statistically significantly different from the models from either prototype.

47 OTHER INSTRUMENTATION↗

Stopping to reflect: Asymptotic static moving mirrors as quantum analogs of classical radiation

Radiation from an accelerating charge is a basic process that can serve as an intersection between classical and quantum physics. We present two exactly soluble electron trajectories that permit analysis of the radiation emitted, exploring its time evolution and spectrum by analogy with the moving mirror model of the dynamic Casimir effect. These classical solutions are finite energy, rectilinear (nonperiodic), asymptotically zero velocity worldlines with corresponding quantum analog beta Bogolyubov coefficients. One of them has an interesting connection to uniform acceleration and Leonardo da Vinci's water pitcher experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗

Impact of Water Ingress on Molybdenum Thin Films and Its Effect on Cu(In,Ga)Se 2 Solar Cells

Solar cell degradation can occur through many pathways and at the different stages of the fabrication process, notably because of the water condensation. In the case of Cu(In,Ga)Se 2 (CIGS) solar cells, the impact of water ingress on Mo back contact can be substantial for not only the film itself but also the device performance and reliability. Microstructural modification, change in film morphology, loss in reflectance, and increased resistivity because of the moisture ingress were observed via X-ray diffraction, transmission electron microscopy, spectroscopic ellipsometry, and four-point probe measurements, respectively. Secondary ion mass spectrometry measurements revealed drastic changes in the alkali (Na, K) profiles in both the CIGS and Mo layers, likely because of the modification of their diffusion coefficient through Mo. This, in turn, negatively impacts the solar cell efficiency by decreasing both fill factor and open-circuit voltage. Additional experiments, modifying the substrate and utilizing a NaF postdeposition treatment, highlight the mechanisms of degradation as being due to both a modification of the Mo/CIGS interface and the lack of alkali diffusion after water ingress.

14 SOLAR ENERGY↗

Temperature-dependent spectral emittance of bauxite and silica particle beds

Bauxite and silica particles are candidate materials for solar thermal energy storage at high temperatures. The temperature-dependent emittance of packed beds with bauxite and silica particles was measured using a newly upgraded emissometer at wavelengths 2 μm ≤ λ ≤ 16 μm and temperatures up to ~730 K. The room-temperature emittance was obtained from the measured directional-hemispherical reflectance. A fused silica disc was used to test the emissometer by comparing the measured spectral emittance with the calculated emittance from a fitted Lorentz oscillator model. For the polycrystalline silica particles and the fused silica disc, the measured emittance increases with temperature in the mid-infrared region. The underlying mechanism is interpreted as the temperature-dependent damping coefficient in the Lorentz oscillator model. Two types of bauxite particles with different compositions and sizes were investigated. For λ > 10 μm, the measured emittance at elevated temperatures is higher than that at room temperature. In the region 2 μm < λ < 6 μm, the temperature dependence varies for different types of particles. The total emittance of bauxite particle beds was calculated by spectral integration using Planck’s distribution at the prescribed temperature. As a result, the calculated total emittance is between 0.89 and 0.96, but it does not change monotonically with temperature.

42 ENGINEERING↗

Pebble Tanker Model for Nuclear Criticality Safety Needs

This report documents a study performed to investigate the requirements for criticality safety benchmark experiments for high-assay, low-enriched uranium (HALEU) fuel in transportation applications. In this work, an exploratory application model, the “Pebble Tanker,” was developed to represent TRISO fuel in a transportation scenario for an analysis of the validation basis in industrial quantities. An aspect of the criticality validation process involves assessing the “similarity” between application and experimental benchmark systems through an integral index parameter evaluation. Here, this includes propagating nuclear data uncertainties and calculating a correlation coefficient (hereinafter referred to as “c k ”) to evaluate the similarity of benchmark experiments compared with the application Pebble Tanker model. Finding sufficient critical benchmark experiments allows for the evaluation of bias and bias uncertainty, thus determining the upper subcritical limit (USL) of the transportation package. A target k eff of ~0.94 was used in this work to establish appropriate modeling conditions, reflecting a reasonable estimate for a USL. Two container models were investigated: one with the Hermes-type pebble and one with the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering only fuel, containment structure, and either water or air. This allows a focus on the underlying physics of applications involving TRISO fuel pebbles using the Pebble Tanker model. A crucial consideration is the transport package's ability to safely hold pebbles while flooded, maintaining subcritical conditions. Tools available in the SCALE 6.3.1 suite—the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences—were employed for neutronics and sensitivity and uncertainty (S/U) analysis of the Pebble Tanker. Findings demonstrated sufficient available critical experiment benchmarks to perform a validation of the Pebble Tanker in the most reactive state, i.e., when the Tanker is flooded.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗

Optical Classification of Lower Amazon Waters Based on In Situ Data and Sentinel-3 Ocean and Land Color Instrument Imagery

Optical water types (OWTs) were identified from an in situ dataset of concomitant biogeochemical and optical parameters acquired in the Amazon River and its tributaries, in the Lower Amazon region, at different hydrological conditions from 2014 to 2017. A seasonal bio-optical characterization was performed. The k-means classification was applied to the in situ normalized reflectance spectra (rn(λ)), allowing the identification of four OWTs. An optical index method was also applied to the rn(λ) defining the thresholds of the OWTs. Next, level-3 Sentinel-3 Ocean and Land Color Instrument images representative of the seasonal discharge conditions were classified using the identified in situ OWTs as reference. The differences between Amazon River and clearwater tributary OWTs were dependent on the hydrological dynamics of the Amazon River, also showing a strong seasonal variability. Each OWT was associated with a specific bio-optical and biogeochemical environment assessed from the corresponding absorption coefficient values of colored dissolved organic matter (aCDOM) and particulate matter (ap), chlorophyll-a and suspended particulate matter (SPM) concentrations, and aCDOM/ap ratio. The rising water season presented a unique OWT with high SPM concentration and high relative contribution of ap to total absorption compared to the other OWTs. This bio-optical characterization of Lower Amazon River waters represents a first step for developing remote sensing inversion models adjusted to the optical complexity of this region.

54 ENVIRONMENTAL SCIENCES↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

Simulation of the impact of particle recycling on the plasma in MPS-LD device based on the BOUT++ LPD module

A linear plasma device (LPD) module has been developed under the BOUT++ framework to simulate plasma transport in the MPS-LD. However, previously, the LPD module used a simplistic neutral particle model that only includes particle density and velocity, which prevents the full understanding of the plasma-neutrals interactions. In this work, we further optimize the neutral model by using a more complete neutral fluid model containing the continuity equation, momentum equation, and energy equation. The reactions such as charge exchange, excitation, and radiation collisions are included. Since the neutral particle source is mainly provided by particle recycling from the target, a particle recycling model is employed, which includes both fast reflection and slow thermal release. The upgraded LPD module is applied to simulate the argon (Ar) discharge experiment of MPS-LD, and the benchmark against experiment measurement and SOLPS-ITER simulation results are presented. Good agreements are obtained, showing the validation of the upgraded module. After that, the impact of particle recycling on Ar plasma is investigated. It is found that a higher recycling coefficient (R) promotes the achievement of high-density plasma at the target. In conclusion, the recycled Ar atoms change target plasma pressure as well as plasma-neutral collisions, which both contribute to plasma momentum loss, thus promoting the rollover of ion flux to the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the Impact of Temporal and Directional Traffic Distribution on Crash Frequencies

Safety Performance Functions (SPFs) are mathematical models that establish relationships between the frequency of various crash types and site-specific characteristics, serving as essential tools for traffic safety analysis and roadway design. Traditional SPFs, however, often overlook the temporal fluctuations in traffic flow (such as peak-hour surges) and directional imbalances between opposing traffic streams. These traffic patterns can exacerbate congestion, disrupt driver behavior, and create unexpected conflict points, potentially leading to increased crash frequencies and more severe accidents. In light of this gap, this study aims to explore the potential of incorporating K-factors (representing peak-hour traffic proportions) and D-factors (reflecting the imbalance of directional traffic) into the development of SPFs to assess whether these factors can effectively represent the impact of temporal and spatial traffic distribution on roadway safety. Using crash data from Pennsylvania urban-suburban collector roadways, it is found that the D-factor plays a significant role in predicting the frequency of total crashes, fatal + injury crashes, and angle crashes, with positive coefficient signs indicating that higher directional imbalances correspond to increased crash risks. Similarly, the K-factor emerges as a critical predictor for fatal + injury crashes and rear-end crashes, with negative coefficients suggesting that a more pronounced traffic peak is associated with a reduction in expected crash frequencies. These results highlight the importance of accounting for uneven traffic distribution in both time and direction when developing SPFs, offering deeper insights into crash patterns and supporting more effective safety interventions and roadway designs.

Xu, Guanhao [ORNL] (ORCID:0000000214326357)↗

Temperature evolution of the phonon dynamics in the Kitaev spin liquid

In this work, we present a study of the phonon dynamics in the honeycomb Kitaev spin model at finite temperatures. We show that the fractionalized spin excitations of the Kitaev spin liquid, the itinerant Majorana fermions, and static Z 2 fluxes have distinct effects on the phonon dynamics, which makes the phonon dynamics a promising tool for exploring the Kitaev spin liquid candidate materials. In particular, we focus on the signature of the fractionalized excitations in the thermodynamic behavior of the sound attenuation and the phonon Hall viscosity: the former describes the phonon decay into the fractionalized excitations, and the latter is the leading order time-reversal symmetry breaking effect on the acoustic phonon. We find that the angular dependence of the attenuation coefficient and its magnitude is modified by the thermal excitation of the Z 2 fluxes. The strength of this effect strongly depends on the relative magnitude of the sound velocity and the Fermi velocity characterizing the low-energy Majorana fermions. We also show that the Hall viscosity is strongly suppressed by the increase of the density of the Z 2 fluxes at finite temperatures. All these observations reflect the effects of the emergent disorder on the Majorana fermions introduced by the Z 2 fluxes. Our analysis is based on the complementary analytical calculations in the low-temperature zero-flux sector, and numerical calculations in the inhomogeneous flux sectors at intermediate and high temperatures with stratified Monte Carlo (strMC) method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sulfur isotope fractionation between aqueous and carbonate-associated sulfate in abiotic calcite and aragonite

Sulfate (SO 4 2– ) incorporated into calcium carbonate minerals enables measurements of sulfur (S) isotope ratios in carbonate rocks. This Carbonate Associated Sulfate (CAS) in marine carbonate minerals is thought to faithfully represent the S isotope composition of the seawater sulfate incorporated into the mineral, with little or no S isotope fractionation in the process. However, comparison between different calcifying species reveals both positive and negative S isotope fractionation between CAS and seawater sulfate, and a large range of S isotope ratios can be found within a single rock sample, depending on the component measured. To better understand the isotopic effects associated with sulfate incorporation into carbonate minerals, we precipitated inorganic calcite and aragonite over a range covering more than two orders of magnitude of sulfate concentration and precipitation rate. Coupled measurements of CAS concentration, S isotope composition and X-ray absorption near-edge spectra (XANES) permit characterization and explanation of the observed dependence of S isotope fractionation between CAS and aqueous sulfate (CAS-SO 4 2– isotope fractionation) on sulfate concentration and precipitation rate. In aragonite, the CAS-SO 4 2– isotope fractionation is 1.0 ± 0.3‰ and independent of the sulfate (and CAS) concentration. In contrast, we note the CAS-SO 4 2– isotope fractionation in calcite covaries strongly with the sulfate concentration and weakly with the precipitation rate, between values of 1.3 ± 0.1 and 3.1 ± 0.6‰. We suggest that the correlation between aqueous sulfate concentration and CAS-SO 4 2– isotope fractionation in calcite reflects a dependence of the equilibrium S isotope fractionation on the concentration of CAS, through the effect of the sulfate impurity on the carbonate mineral’s energetic state.

58 GEOSCIENCES↗

The Cosmological Bootstrap: Spinning correlators from symmetries and factorization

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗

Evaluation of WSR-88D Level III and MRMS Rainfall Estimates Against Rain Gauge Observations at the Savannah River Site

Rainfall data for the Savannah River Site (SRS) has been historically measured by rain gauges. These instruments serve as ground truth for most climatological and weather applications; however, gauge measurements are prone to errors or biases under certain weather conditions. Rainfall estimates from radar reflectivity values have been developed and improved over the years and serve as an alternative or supplement for gauge measurements. This study compares measurements from tipping bucket rain gauges with rainfall estimates from the NOAA NWS WSR-88D Level III hourly rainfall product and the Multi-Radar Multi-Sensor (MRMS) gauge-corrected precipitation estimate. Results show good agreement between radar derived amounts and ground measurements, with MRMS values showing correlation coefficients of 0.60-0.95 and RMSE values of less than 1.0 cm (0.4 in). The WSR-88D Level III estimates result in correlation coefficients of 0.46-0.86 and RMSE values of less than 1.3 cm (0.5 in). Few outliers are observed for each data pair and are evaluated against precipitation classification products (WSR-88D Hybrid Hydrometeor Classification product and MRMS Precipitation Flag product). The rain gauges used in this study are not part of the Hydrometeorological Automated Data System (HADS) network used to correct the MRMS estimates and therefore, results of this work provide an independent validation of the MRMS gauge correction scheme.

54 ENVIRONMENTAL SCIENCES↗