Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Calibrating the absolute magnitude of type Ia supernovae in nearby galaxies using [O II ] and implications for H 0

The present state of cosmology is facing a crisis where there is a fundamental disagreement in measurements of the Hubble constant ($H_{0}$), with significant tension between the early and late Universe methods. Type Ia supernovae (SNe Ia) are important to measuring $H_{0}$ through the astronomical distance ladder. However, there remains potential to better standardize SN Ia light curves by using known dependencies on host galaxy properties after the standard light curve width and colour corrections have been applied to the peak SN Ia luminosities. To explore this, we use the 5-yr photometrically identified SNe Ia sample obtained by the Dark Energy Survey, along with host galaxy spectra obtained by the Australian Dark Energy Survey. Using host galaxy spectroscopy, we find a significant trend with the equivalent width (EW) of the [O II ] $\lambda \lambda$ 3727, 29 doublet, a proxy for specific star formation rate, and Hubble residuals. We find that the correlation with [O II ] EW is a powerful alternative to the commonly used mass step after initial light-curve corrections. Applying this [O II ] EW correction to 20 SNe Ia in calibrator galaxies observed with WiFeS, we examined the impact on SN Ia absolute magnitudes and $H_{0}$. Our [O II ] EW corrections result in $H_{0}$ values ranging between 73.04 and 73.51 $km$ $s^{-1} \, Mpc^{-1}$, with a combined statistical and systematic uncertainty of $\sim 1.31$ $km$ $s^{-1} \, Mpc^{-1}$. However, even with this additional correction, the impact of host galaxy properties in standardizing SNe Ia appears limited in reducing the current tension ($\sim 5\sigma$) with the cosmic microwave background result for $H_{0}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evaluating Synthetic Smart Meter Locations For Communication System Modeling

An important part of good communication models for smart grid applications, particularly when wireless protocols are used, is the location information of the nodes. We’ve developed metrics to assist in evaluating whether a given model with such information is representative of feeders found in the world. These metrics are applied to a set of feeder models with known good location information for the smart meter locations and compared to a set of feeder models with similar location information that is not expected to be representative of real-world feeders. The comparison reveals that the suspect models do not pass statistical tests utilizing the developed metrics, providing initial validation of the analysis technique.

smart meter, feeder model, GridLAB-D↗
DOE OSTI2023

LPG Component Leak Frequency Estimation

Liquefied petroleum gas (LPG) is used in heating, cooking, and as a vehicle fuel (called autogas). A safety risk assessment may be needed to assess potential hazard scenarios and inform the regulations, codes, and standards that apply to LPG facilities, such as autogas refueling facilities. The frequency of unintended releases in an LPG system is an important aspect of a system quantitative risk assessment. This report documents estimation of leakage frequencies for individual components of LPG systems. These frequencies are described using uncertainty distributions obtained with Bayesian statistical methods, generic data, and LPG data which were publicly available. There was a lack of LPG data in the literature, so frequencies for most components were developed with generic data and should be used cautiously; without additional information about component leak frequencies in LPG systems, it is not known whether these generic frequencies may be conservative or non-conservative.

02 PETROLEUM↗
DOE OSTI2021

Visible-light Phase Curves from the Second Year of the TESS Primary Mission

We carried out a systematic study of full-orbit phase curves for known transiting systems in the northern ecliptic sky that were observed during Year 2 of the TESS primary mission. We applied the same methodology for target selection, data processing, and light-curve fitting as we did in our Year 1 study. Out of the 15 transiting systems selected for analysis, seven—HAT-P-7, KELT-1, KELT-9, KELT-16, KELT-20, Kepler-13A, and WASP-12—show statistically significant secondary eclipses and day–night atmospheric brightness modulations. Small eastward dayside hot-spot offsets were measured for KELT-9b and WASP-12b. KELT-1, Kepler-13A, and WASP-12 show additional phase-curve variability attributed to the tidal distortion of the host star; the amplitudes of these signals are consistent with theoretical predictions. We combined occultation measurements from TESS and Spitzer to compute dayside brightness temperatures, TESS-band geometric albedos, Bond albedos, and phase integrals for several systems. The new albedo values solidify the previously reported trend between dayside temperature and geometric albedo for planets with 1500 K < T {sub day} < 3000 K. For Kepler-13Ab, we carried out an atmospheric retrieval of the full secondary eclipse spectrum, which revealed a noninverted temperature–pressure profile, significant H{sub 2}O and K absorption in the near-infrared, evidence for strong optical atmospheric opacity due to sodium, and a confirmation of the high geometric albedo inferred from our simpler analysis. We explore the implications of the phase integrals (ratios of Bond to geometric albedos) for understanding exoplanet clouds. We also report updated transit ephemerides for all of the systems studied in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

ResStock Measure Documentation: Residential Variable-Speed Geothermal Heat Pump (4.4 COP, 30.9 EER)

The goal of this work is to develop energy efficiency, demand flexibility, and other retrofit end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to retrofits that can be applied to buildings during modeling. An "end-use savings shape" is the difference in energy consumption between a baseline building and a building with an energy efficiency, demand flexibility, or other retrofit measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ResStock (TM) is a highly granular, physics-based, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the residential building stock across the United States. The baseline model intends to represent the U.S. residential building stock as it existed in 2018. Technical documentation for the inputs and assumptions in the baseline building stock model is available in Reyna et al. (2025). Calibration and validation of the baseline model results are available in the final technical report of the End-Use Load Profiles project (Wilson et al. 2022). This documentation focuses on a single end-use savings shape measure: Residential Variable-Speed Geothermal Heat Pump (GHP). This document provides the relevant new modeling information for variable-speed systems not previously covered in either the single-stage or two-stage documents. Variable-speed GHPs represent the most efficient option available for this technology: They provide the most savings, with up to 46% for the applicable portion of the housing stock, compared to 31% for less efficient single-stage GHPs. Additional results shown here detail how the savings change for sections of the housing stock with different base heating fuels and in different climate zones, and they show the savings potential by state for both heating and cooling. Utility bills and electric panel impacts are also shown and discussed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Systematic identification of cis -regulatory variants that cause gene expression differences in a yeast cross

Sequence variation in regulatory DNA alters gene expression and shapes genetically complex traits. However, the identification of individual, causal regulatory variants is challenging. Here, we used a massively parallel reporter assay to measure the cis-regulatory consequences of 5832 natural DNA variants in the promoters of 2503 genes in the yeast Saccharomyces cerevisiae. We identified 451 causal variants, which underlie genetic loci known to affect gene expression. Several promoters harbored multiple causal variants. In five promoters, pairs of variants showed non-additive, epistatic interactions. Causal variants were enriched at conserved nucleotides, tended to have low derived allele frequency, and were depleted from promoters of essential genes, which is consistent with the action of negative selection. Causal variants were also enriched for alterations in transcription factor binding sites. Models integrating these features provided modest, but statistically significant, ability to predict causal variants. This work revealed a complex molecular basis for cis-acting regulatory variation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

The Argonne SuperGel for CBRN Decontamination - 20160

The Argonne SuperGel was developed between 2003 and 2015 to fill a gap in our nation's capability to quickly decontaminate important structures following a radiological or nuclear release event. Specifically, the decontamination technology was developed to minimize damage to monuments, high valued structures, and critical infrastructure while reducing environmental and health impacts. An important criterion during its development, common reagents were employed that could be easily acquired in order to minimize the timeline for its deployment. Its current formulation uses off-the-shelf super-absorbing hydrogels common to the food and agricultural industry and common salts. Over the years, two formulations of the Argonne SuperGel have been developed to specifically target radioactive cesium contaminations and then, more generically, actinide and fission product contaminations. A biodegradable derivative of phosphoric acid is used in small amounts to promote the removal of insoluble actinide species. The original order of unit operations for the SuperGel technology was 1) applying the wash solution to promote mobility of the contaminants on the surface, 2) applying the hydrogel to absorb the wash solution and contaminants, and 3) removing the hydrogel for disposal. During the course of our studies, we tested a new formulation of the gel technology that eliminated the need for a separate step to apply the wash solution without a compromise in decontamination factors. Greater than 70% and >95% of Cs-137 were removed from concrete and tile coupons, respectively, after two decontaminations, with the best results using wash solutions of 1.0 M KCl and 1.0 M NH{sub 4}Cl formulations. We found no statistical difference between results at 30 deg. C and 40 deg. C or at 70% and 90% relative humidity. We tested the gel under strong UV irradiation to simulate extreme environments. While a tropical noontime UV flux did not adversely affect the gel system decontamination of Cs-137 from concrete or tile coupons, the noontime flux resulted in significant dehydration of the gel after two hours of constant exposure. The rehydration of the gel by the humid air (at 90% humidity and 40 deg. C) was slow and may not be sufficient to offset dehydration by the sun. Americium decontamination from concrete was 70% after optimizing the phosphoric acid derivative and carbonate salt concentration in the wash solution. Since its development, we have had the opportunity to test the SuperGel in the removal of legacy contaminations in hot cell facilities and former glovebox facilities at Argonne. This has provided a unique opportunity to evaluate the SuperGel on a range of contaminants outside the original specifications for its use. We will report on the origin of the gel formulation, some highlighted experimental data including independent testing by the US EPA which has never been reported, and its efficacy for removing legacy alpha contaminations and its potential use for removal of chemical and biological hazardous agents. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Development of online pulse-finding algorithms for the Mu2e Stopping Target Monitor and an evaluation of the signal to background ratio for the Mu2e stopped-muon rate determination

This thesis presents a determination of the acceptance of the Mu2e Stopping Target Monitor (STM) detector for both signal and background events. The rate of X-rays of a specific energy in the STM is used to determine the rate of stopped muons in the aluminium stopping target, which is used as the normalisation for all Mu2e measurements, particularly the branching ratio for a charged lepton flavour violating (cLFV) process. The acceptance for the 347\,keV $2p-1s$ aluminium signal X-ray is determined to be $(8.03 \pm 0.05) \times 10^{-10}$\,cm$^{-2}$ and the acceptance for the background photons (at 347\,keV) to be $(3.19 \pm 0.21) \times 10^{-8}$\,cm$^{-2}$. This results in a signal-to-background ratio of approximately two which is significantly higher than had previously been estimated. These signal-to-background rates are used to determine the time required to achieve a statistical uncertainty of 10\% on the number of 347\,keV X-rays as a function of the detector resolution and back ground rate. For a resolution of 2\,keV, and the background rate determined in this thesis, approximately one minute of data is required. The three corrections needed to be applied to the measured 347\,keV rate to determine the rate of stopped muons are also determined and potential sources of bias arising from the background shape, background fit range and energy calibration are investigated and found to be insignificant. The raw energy of the particles in the STM detector comes from a 16-bit ADC digitising at $\sim 300$\,MHz. It is not possible to store all the raw digitisations and online algorithms have been developed to rapidly select pulses in the detector and reconstruct their energy. The performance of these algorithms has been characterised with simulated data as a function of detector rate and resolution and optimum operating points for the algorithms have been defined. The integrity of the algorithms has also been demonstrated and characterised on test-beam data at high r ates and with radioactive sources at low rates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Detection Limits of Low-mass, Long-period Exoplanets Using Gaussian Processes Applied to HARPS-N Solar Radial Velocities

Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s{sup −1}. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s{sup −1}) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s{sup −1}) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Intrinsic and extinction colour components in SNe Ia and the determination of R V

Context.The colour fluctuations of type Ia supernovae (SNe Ia) include intrinsic and extrinsic components, which both contribute to the observed variability. Previous works proposed a statistical separation of these two contributions, but the individual intrinsic colour contributions of each SN Ia were not extracted. In addition, a large uncertainty remains on the value of the parameterR V , which characterises the dust extinction formula. Aims.Leveraging the known parameterisation of the extinction formula for dust in our Galaxy, and applying it to the host galaxy of SNe Ia, we propose a new method of separation –valid for each SN– using the correlations between colour fluctuations. This also allows us to derive a well-constrained value of the extinction parameterR V with different, possibly smaller systematic errors. We also define a three-dimensional space of intrinsic colour fluctuations. Methods.The key ingredients in this attempt at separating the intrinsic and extinction colour components for each SN –and subsequently measuringR V – are the assumption of a linearized dependence of magnitude on the extinction component of colour, a one-dimensional extra-intrinsic colour space (in addition to CaII H&Kλ3945 and Si IIλ4131 contributions) over four independent colours, and the absence of correlation between the intrinsic and extrinsic variabilities. Results.We show that a consistent solution is found under the previous assumptions, but the observed systematic trends point to a (small) inadequacy of the extinction formula. Once corrected, all systematic extinction effects can be cancelled by choosing a single scaling of the extinction colour component as well as an appropriate value ofR V = 2.181 ± 0.117. The observed colours are described within an accuracy of 0.025 mag. The resulting magnitude variability is 0.13 over allUBVRIbandpasses, and this fluctuation is shown to be independent of the bandpass to within 0.02 mag.

Astronomy & Astrophysics↗
DOE OSTI2022

A deep, multi-epoch Chandra HETG study of the ionized outflow from NGC 4051

Actively accreting supermassive black holes significantly impact the evolution of their host galaxies, truncating further star formation by expelling large fractions of gas with wide-angle outflows. The X-ray band is key to understanding how these black hole winds affect their environment, as the outflows have high temperatures ~10 5–8 K). We have developed a Bayesian framework for characterizing active galactic nucleus outflows with an improved ability to explore parameter space and perform robust model selection. We applied this framework to a new 700 ks and an archival 315 ks Chandra High Energy Transmission Gratings observation of the Seyfert galaxy NGC 4051. We have detected six absorbers intrinsic to NGC 4051. These wind components span velocities from 400 to 30000 km s -1 . We have determined that the most statistically significant wind component is purely collisionally ionized, which is the first detection of such an absorber. This wind has T ≈ 10 7 K and v ≈ 880 km s -1 and remains remarkably stable between the two epochs. Other slow components also remain stable across time. Fast outflow components change their properties between 2008 and 2016, suggesting either physical changes or clouds moving in and out of the line of sight. For one of the fast components, we obtain one of the tightest wind density measurements to date, log n/(cm -3 ) = 13.0$^{+0.01}_{-0.02}$, and determine that it is located at ~240 gravitational radii. The estimated total outflow power surpasses 5 percent of the bolometric luminosity (albeit with large uncertainties) making it important in the context of galaxy–black hole interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dark Energy Survey Year 6 results: Galaxy-galaxy lensing

Here, we present galaxy-galaxy lensing measurements from the full six years of data from the Dark Energy Survey (DES Y6), covering 4031 deg 2 and used in the DES Y6 3×2 pt cosmological analysis. We use the Maglim++ lens sample, containing ∼9 million galaxies divided into six redshift bins, and the Metadetection source catalog, including ∼140 million galaxies divided into four redshift bins. The mean tangential shear signal achieves a total signal-to-noise ratio (S/N) of 173, corresponding to a 17% improvement over DES Y3. After applying the scale cuts used in the cosmological analysis, with 𝑅 min =6 Mpc/ℎ (4 Mpc/ℎ) for the linear (nonlinear) galaxy-bias model, the S/N is reduced to 75 (90). A comprehensive suite of validation tests demonstrates that the measurement is robust against observational and astrophysical systematics at the statistical precision required for the DES Y6 analysis. Although not used in the main cosmological analysis, we extract high signal-to-noise geometric shear-ratio measurements from the galaxy-galaxy lensing signal on small angular scales. These measurements provide an internal consistency check on the photometric redshift distributions and shear calibration used in the 3 ×2 pt analysis.

cosmological parameters↗
DOE OSTI2023

Lattice quantum chromodynamics at large isospin density

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of n-pion correlation functions for n ϵ {1, … , N} using a single N × N matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 charged pions using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass m π = 170 MeV and spacetime volumes of (4.4 3 × 8.8) fm 4 and (5.8 3 × 11.6) fm 4 . We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well-approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential μ I ~ 1.5m π , signaling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound, c s , in the medium is seen to exceed the ideal-gas (conformal) limit ($c^{2}_{s} ≤ 1/3)$ over a wide range of chemical potential before falling towards the asymptotic expectation at μ I ~ 15m π . These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Comprehensive Margin Assessment of the ASME Section III, Division 5, Class A Primary Load Design Rules

This report provides a comprehensive margin assessment of the ASME Section III, Division 5, Subsection HB, Subpart B rules for the design of high temperature reactor components against load controlled stress limits. These rules, often called the Code primary load design rules, provide protection against creep rupture and plastic collapse under steady conditions. The method adopted here is to compare the deterministic ASME design life, according to the primary load design rules, against the expected, statistical service life of the component, considering creep rupture as the relevant failure mode. The report provides the design margin of the ASME rules in terms of the probability of premature failure -- the probability that the component fails in service before reaching the ASME design life. To complete this assessment, the report describes the development of a statistical creep life assessment procedure accounting for variations in the component loading, material creep rate, and the material rupture time. This includes a novel method for correlating creep rupture data using a Gaussian process accounting for heat-to-heat variation in the Larson-Miller time-temperature parameter. The report applies the complete margin assessment process to quantify the design margin inherent in the ASME primary load design rules as a function of temperature for relevant component geometries and materials. The overall conclusion is that the ASME rules are conservative, providing designs with very small probability of premature failure, but some optimization of the Code design rules is possible to provide a more uniform design margin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Evaluation of Four Ground-based Retrievals of Cloud Droplet Number Concentration in Marine Stratocumulus with Aircraft In Situ Measurements

Cloud droplet number concentration (N d ) is crucial for understanding aerosol-cloud interactions (ACI) and associated radiative effects. We present evaluations of four ground-based N d retrievals based on comprehensive datasets from the Atmospheric Radiation Measurements (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign. The N d retrieval methods use ARM ENA observatory ground-based remote sensing observations from a Micropulse lidar, Raman lidar, cloud radar, and the ARM NDROP Value-added Product (VAP), all of which also retrieve cloud effective radius (r e ). The retrievals are compared against aircraft measurements from the Fast-Cloud Droplet Probe (FCDP) and the Cloud and Aerosol Spectrometer (CAS) obtained from low-level marine boundary layer clouds on 12 flight days during summer and winter seasons. Additionally, the in situ measurements are used to validate the assumptions and characterizations used in the retrieval algorithms. Statistical comparisons of the probability distribution function (PDF) of the N d and cloud r e retrievals with aircraft measurements demonstrate that these retrievals align well with in situ measurements for overcast clouds, but they may substantially differ for broken clouds or clouds with low liquid water path (LWP). The retrievals are applied to four years of ground-based remote sensing measurements of overcast marine boundary layer clouds at the ARM ENA observatory to find that N d (r e ) values exhibit seasonal variations, with higher (lower) values during the summer season and lower (higher) values during the winter season. The ensemble of various retrievals using different measurements and retrieval algorithms such as those in this paper can help to quantify N d retrieval uncertainties and identify reliable N d retrieval scenarios. Of the retrieval methods, we recommend using the using the Micropulse lidar-based method given its good agreement with in situ measurements, it has less sensitivity to issues arising from precipitation and low cloud LWP/optical depth, and it has broad applicability by functioning for both day and nighttime conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A modeling pipeline to relate municipal wastewater surveillance and regional public health data

As COVID-19 becomes endemic, public health departments require improved passive indicators, which are independent of voluntary testing data, to estimate the prevalence of COVID-19 in local communities. Quantification of SARS-CoV-2 RNA from wastewater has the potential to be a powerful passive indicator. However, connecting measured SARS-CoV-2 RNA to community prevalence is challenging. We have developed a generalized methodology to improve the predictive power of wastewater measurements and applied it to data collected from treatment plants in the Chicago area. We built and compared a set of multi-linear regression models, which incorporate pepper mild mottle virus (PMMoV) as a population biomarker, Bovine coronavirus (BCoV) as a recovery control, and wastewater system flow rate into a corrected estimate for SARS-CoV-2 RNA concentration. For our data, models with BCoV performed better than those with PMMoV, but all model estimates of prevalence significantly improved correlation compared to doing no correction. We also investigate the utility of RNA measurements in wastewater as a leading indicator of COVID-19 trends. We do this in a rolling manner for corrected wastewater data and for other prevalence indicators, and statistically compare the temporal relationship between new increases in the wastewater data and those in other prevalence indicators. We find that wastewater trends often lead other COVID-19 indicators in predicting new surges.

60 APPLIED LIFE SCIENCES↗