Results for “Correlation coefficient” — Page 14

DOE OSTI2022

A new method for inferring city emissions and lifetimes of nitrogen oxides from high-resolution nitrogen dioxide observations: a model study

We present a new method to infer emissions and lifetimes of nitrogen oxides (NO x ) based on tropospheric nitrogen dioxide (NO 2 ) observations together with reanalysis wind fields for cities located in polluted backgrounds. Since the accuracy of the method is difficult to assess due to lack of “true values” that can be used as a benchmark, we apply the method to synthetic NO 2 observations derived from the NASA-Unified Weather Research and Forecasting (NU-WRF) model at a high horizontal spatial resolution of 4 km × 4 km for cities over the continental United States. We compare the inferred emissions and lifetimes with the values given by the NU-WRF model to evaluate the method. The method is applicable to 26 US cities. The derived results are generally in good agreement with the values given by the model, with the relative differences of 2 % ± 17 % (mean ± standard deviation) and 15 % ± 25 % for lifetimes and emissions, respectively. Our investigation suggests that the use of wind data prior to the satellite overpass time improves the performance of the method. The correlation coefficients between inferred and NU-WRF lifetimes increase from 0.56 to 0.79 and for emissions increase from 0.88 to 0.96 when comparing results based on wind fields sampled simultaneously with satellite observations and averaged over 9 h data prior to satellite observations, respectively. We estimate that uncertainties in NO x lifetime and emissions arising from the method are approximately 15 % and 20 %, respectively, for typical (US) cities. The total uncertainties reach up to 43% (lifetimes) and 45% (emissions) by considering the additional uncertainties associated with satellite NO 2 observations and wind data. We expect this new method to be applicable to NO 2 observations from the TROPOspheric Monitoring Instrument (TROPOMI) and geostationary satellites, such as Geostationary Environment Monitoring Spectrometer (GEMS) or the Tropospheric Emissions: Monitoring Pollution (TEMPO) instrument, to estimate urban NO x emissions and lifetimes globally.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Summertime ozone pollution in China affected by stratospheric quasi-biennial oscillation

In recent years, the near-surface ozone (O 3 ) level has been rising fast in China, with increasing damage to human health and ecosystems. In this study, the impact of stratospheric quasi-biennial oscillation (QBO) on interannual variations in summertime tropospheric O 3 over China is investigated based on GEOS-Chem model simulations and satellite retrievals. QBO has a significant positive correlation with near-surface O 3 concentrations over central China (92.5–112.5° E, 26–38° N) when the sea surface temperature (SST) over the eastern tropical Pacific is warmer than normal, with a correlation coefficient of 0.53, but QBO has no significant effect on O 3 under the cold SST anomaly. Compared to the easterly phase of QBO, the near-surface O 3 concentrations have an increase of up to 3 ppb (5% relative to the average) over central China during its westerly phase under the warm SST anomaly. O 3 also increases above the surface and up to the upper troposphere, with a maximum increase of 2–3 ppb (3%–5%) in 850–500 hPa over central China when comparing westerly phase to easterly phase. Process-based analysis and sensitivity simulations suggest that the O 3 increase over central China is mainly attributed to the anomalous downward transport of O 3 during the westerly phase of QBO when a warm SST anomaly occurs in the eastern tropical Pacific, while the local chemical reactions and horizontal transport processes partly offset the O 3 increase. This work suggests a potentially important role of QBO and the related vertical transport process in affecting near-surface O 3 air quality, with an indication for O 3 pollution prediction and prevention.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Measurement report: The Fifth International Workshop on Ice Nucleation phase 1 (FIN-01): intercomparison of single-particle mass spectrometers

Knowledge of the chemical composition and mixing state of aerosols at a single-particle level is critical for gaining insights into atmospheric processes. One common tool to make these measurements is single-particle mass spectrometry. There remains a need to compare the performance of different single-particle mass spectrometers (SPMSs). An intercomparison of SPMSs was conducted at the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) chamber at the Karlsruhe Institute of Technology (KIT) in November 2014, as part of the first phase of the Fifth International Workshop on Ice Nucleation (FIN-01). In this paper we compare size distributions and mass spectra of atmospherically relevant particle types measured by five SPMSs. These include different minerals, desert and soil dusts, soot, bioaerosol (Snomax; protein granule), secondary organic aerosol (SOA), and SOA-coated mineral particles. Most SPMSs reported similar vacuum aerodynamic diameter (d va ) within typical instrumental ranges from ~100–200 nm (lower limit) to ~2–3 µm (upper limit). In general, all SPMSs exhibited a wide dynamic range (up to ~10 3 ) and high signal-to-noise ratio (up to ~10 4 ) in mass spectra. Common spectral features with small diversities in mass spectra were found with high average Pearson's correlation coefficients, i.e., for average positive spectra r avg-pos =0.74 ± 0.12 and average negative spectra r avg-neg =0.67 ± 0.22. We found that instrument-specific detection efficiency (DE) was more dependent on particle size than particle type, and particle identification favored the use of bipolar, rather than monopolar, instruments. Particle classification from “blind experiments” showed that all instruments differentiated SOA, soot, and soil dust and detected subtle changes in the particle internal mixing but had difficulties differentiating among specific mineral types and dusts. This study helps to further understand the capabilities and limitations of the single-particle mass spectrometry technique in general and the specific performance of the instrument in characterizing atmospheric aerosol particles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Cloud condensation nuclei concentrations derived from the CAMS reanalysis

Determining number concentrations of cloud condensation nuclei (CCN) is one of the first steps in the chain in analysis of cloud droplet formation, the direct microphysical link between aerosols and cloud droplets, and a process key for aerosol–cloud interactions (ACI). However, due to sparse coverage of in situ measurements and difficulties associated with retrievals from satellites, a global exploration of their magnitude, source as well as temporal and spatial distribution cannot be easily obtained. Thus, a better representation of CCN numbers is one of the goals for quantifying ACI processes and achieving uncertainty-reduced estimates of their associated radiative forcing. Here, we introduce a new CCN dataset which is derived based on aerosol mass mixing ratios from the latest Copernicus Atmosphere Monitoring Service reanalysis (CAMSRA) in a diagnostic model that uses CAMSRA aerosol properties and a simplified kappa-Köhler framework suitable for global models. The emitted aerosols in CAMSRA are not only based on input from emission inventories using aerosol observations, they also have a strong tie to satellite-retrieved aerosol optical depth (AOD) as this is assimilated as a constraining factor in the reanalysis. Furthermore, the reanalysis interpolates for cases of poor or missing retrievals and thus allows for a full spatiotemporal quantification of CCN numbers. The derived CCN dataset captures the general trend and spatial and temporal distribution of total CCN number concentrations and CCN from different aerosol species. A brief evaluation with ground-based in situ measurements demonstrates the improvement of the modelled CCN over the sole use of AOD as a proxy for CCN as the overall correlation coefficient improved from 0.37 to 0.71. However, we find the modelled CCN from CAMSRA to be generally high biased and find a particular erroneous overestimation at one heavily polluted site which emphasises the need for further validation. The CCN dataset (https://doi.org/10.26050/WDCC/QUAERERE_CCNCAMS_v1, Block, 2023), which is now freely available to users, features 3-D CCN number concentrations of global coverage for various supersaturations and aerosol species covering the years 2003–2021 with daily frequency. This dataset is one of its kind as it offers lots of opportunities to be used for evaluation in models and in ACI studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗

DOE OSTI2023

Evaluation of E3SM land model snow simulations over the western United States

Abstract. Seasonal snow has crucial impacts on climate, ecosystems, and humans, but it is vulnerable to global warming. The land component (ELM) of the Energy Exascale Earth System Model (E3SM) mechanistically simulates snow processes from accumulation, canopy interception, compaction, and snow aging to melt. Although high-quality field measurements, remote sensing snow products, and data assimilation products with high spatio-temporal resolution are available, there has been no systematic evaluation of the snow properties and phenology in ELM. This study comprehensively evaluates ELM snow simulations over the western United States at 0.125∘ resolution during 2001–2019 using the Snow Telemetry (SNOTEL) in situ networks, MODIS remote sensing products (i.e., MCD43 surface albedo product), the spatially and temporally complete (STC) snow-covered area and grain size (MODSCAG) and MODIS dust and radiative forcing in snow (MODDRFS) products (STC-MODSCAG/STC-MODDRFS), and the snow property inversion from remote sensing (SPIReS) product and two data assimilation products of snow water equivalent and snow depth – i.e., University of Arizona (UA) and SNOw Data Assimilation System (SNODAS). Overall the ELM simulations are consistent with the benchmarking datasets and reproduce the spatio-temporal patterns, interannual variability, and elevation gradients for different snow properties including snow cover fraction (fsno), surface albedo (αsur) over snow cover regions, snow water equivalent (SWE), and snow depth (Dsno). However, there are large biases of fsno with dense forest cover and αsur in the Rocky Mountains and Sierra Nevada in winter, compared to the MODIS products. There are large discrepancies of snow albedo, snow grain size, and light-absorbing particle-induced snow albedo reduction between ELM and the MODIS products, attributed to uncertainties in the aerosol forcing data, snow aging processes in ELM, and remote sensing retrievals. Against UA and SNODAS, ELM has a mean bias of −20.7 mm (−35.9 %) and −20.4 mm (−35.5 %), respectively, for spring, and −13.8 mm (−27.8 %) and −10.2 mm (−22.2 %), respectively, for winter. ELM shows a relatively high correlation with SNOTEL SWE, with mean correlation coefficients of 0.69 but negative mean biases of −122.7 mm. Compared to the snow phenology of STC-MODSCAG and SPIReS, ELM shows delayed snow accumulation onset dates by 17.3 and 12.4 d, earlier snow end dates by 35.5 and 26.8 d, and shorter snow durations by 52.9 and 39.5 d, respectively. This study underscores the need for diagnosing model biases and improving ELM representations of snow properties and snow phenology in mountainous areas for more credible simulation and future projection of mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Can reanalysis products outperform mesoscale numerical weather prediction models in modeling the wind resource in simple terrain?

Mesoscale numerical weather prediction (NWP) models are generally considered more accurate than reanalysis products in characterizing the wind resource at heights of interest for wind energy, given their finer spatial resolution and more comprehensive physics. However, advancements in the latest ERA-5 reanalysis product motivate an assessment on whether ERA-5 can model wind speeds as well as a state-of-the-art NWP model – the Weather Research and Forecasting (WRF) Model. We consider this research question for both simple terrain and offshore applications. Specifically, we compare wind profiles from ERA-5 and the preliminary WRF runs of the Wind Integration National Dataset (WIND) Toolkit Long-term Ensemble Dataset (WTK-LED) to those observed by lidars at a site in Oklahoma, United States, and in a United States Atlantic offshore wind energy area. We find that ERA-5 shows a significant negative bias (~-1ms-1) at both locations, with a larger bias at the land-based site. WTK-LED-predicted wind speed profiles show a limited negative bias (~-0.5ms-1) offshore and a slight positive bias (~+0.5ms-1) at the land-based site. On the other hand, we find that ERA-5 outperforms WTK-LED in terms of the centered root-mean-square error (cRMSE) and correlation coefficient, for both the land-based and offshore cases, in all atmospheric stability conditions. We find that WTK-LED's higher cRMSE is caused by its tendency to overpredict the amplitude of the wind speed diurnal cycle. At the land-based site, this is partially caused by wind plant wake effects not being accurately captured by WTK-LED.

17 WIND ENERGY↗

DOE OSTI2024

Quantitative comparison of power production and power quality onshore and offshore: a case study from the eastern United States

Abstract. A major issue in quantifying potential power generation from prospective wind energy sites is the lack of observations from heights relevant to modern wind turbines, particularly for offshore where blade tip heights are projected to increase beyond 250 m. We present analyses of uniquely detailed data sets from lidar (light detection and ranging) deployments in New York State and on two buoys in the adjacent New York Bight to examine the relative power generation potential and power quality at these on- and offshore locations. Time series of 10 min wind power production are computed from these wind speeds using the power curve from the International Energy Agency 15 MW reference wind turbine. Given the relatively close proximity of these lidar deployments, they share a common synoptic-scale meteorology and seasonal variability with lowest wind speeds in July and August. Time series of power production from the on- and offshore location are highly spatially correlated with the Spearman rank correlation coefficient dropping below 0.4 for separation distances of approximately 350 km. Hence careful planning of on- and offshore wind farms (i.e., separation of major plants by > 350 km) can be used reduce the system-wide probability of low wind energy power production. Energy density at 150 m height at the offshore buoys is more than 40 % higher, and the Weibull scale parameter is 2 m s−1 higher than at all but one of the land sites. Analyses of power production time series indicate annual energy production is almost twice as high for the two offshore locations. Further, electrical power production quality is higher from the offshore sites that exhibit a lower amplitude of diurnal variability, plus a lower probability of wind speeds below the cut-in and of ramp events of any magnitude. Despite this and the higher resource, the estimated levelized cost of energy (LCoE) is higher from the offshore sites mainly due to the higher infrastructure costs. Nonetheless, the projected LCoE is highly competitive from all sites considered.

Foody, Rebecca↗

DOE OSTI2021

xsacrcfrspccopol (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

xsacrcfrspcxpol (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

xsacrcfrspccross (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

X-Band Scanning ARM Precipitation Radar, CF/Radial Formatted, Quality Controlled

The X-band Scanning ARM Precipitation Radar (X-SAPR) is an X-band dual-polarization Doppler weather radar. The X-SAPR operates in a simultaneous transmit and receive (STAR) mode, meaning that the transmit signal is split so that power is transmitted on both horizontal and vertical polarizations at the same time. The X-SAPR transmitter is a 200 kW (peak power) magnetron. The receiver is based upon the Vaisala RVP-900 and runs Vaisala’s IRIS software. In addition to the first three Doppler moments (reflectivity, radial velocity, and spectra width), the X-SAPR also provides differential reflectivity, correlation coefficient, and specific differential phases. The dual-polarization variables enable estimates of rainfall rates and identification of precipitation types. At the Southern Great Plains site, three X-SAPRs surround the Central Facility, allowing the use of multi-Doppler velocity retrievals to estimate wind fields.

xsaprcfrqc↗

DOE OSTI2023

Performance of methods for SARS-CoV-2 variant detection and abundance estimation within mixed population samples

The accurate identification of SARS-CoV-2 (SC2) variants and estimation of their abundance in mixed population samples (e.g., air or wastewater) is imperative for successful surveillance of community level trends. Assessing the performance of SC2 variant composition estimators (VCEs) should improve our confidence in public health decision making. Here, we introduce a linear regression based VCE and compare its performance to four other VCEs: two re-purposed DNA sequence read classifiers (Kallisto and Kraken2), a maximum-likelihood based method (Lineage deComposition for Sars-Cov-2 pooled samples (LCS)), and a regression based method (Freyja). We simulated DNA sequence datasets of known variant composition from both Illumina and Oxford Nanopore Technologies (ONT) platforms and assessed the performance of each VCE. We also evaluated VCEs performance using publicly available empirical wastewater samples collected for SC2 surveillance efforts. Bioinformatic analyses were performed with a custom NextFlow workflow (C-WAP, CFSAN Wastewater Analysis Pipeline). Relative root mean squared error (RRMSE) was used as a measure of performance with respect to the known abundance and concordance correlation coefficient (CCC) was used to measure agreement between pairs of estimators. Based on our results from simulated data, Kallisto was the most accurate estimator as it had the lowest RRMSE, followed by Freyja. Kallisto and Freyja had the most similar predictions, reflected by the highest CCC metrics. We also found that accuracy was platform and amplicon panel dependent. For example, the accuracy of Freyja was significantly higher with Illumina data compared to ONT data; performance of Kallisto was best with ARTICv4. However, when analyzing empirical data there was poor agreement among methods and variations in the number of variants detected (e.g., Freyja ARTICv4 had a mean of 2.2 variants while Kallisto ARTICv4 had a mean of 10.1 variants). This work provides an understanding of the differences in performance of a number of VCEs and how accurate they are in capturing the relative abundance of SC2 variants within a mixed sample (e.g., wastewater). Such information should help officials gauge the confidence they can have in such data for informing public health decisions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Analysis of Correlation between Cold Weather Meteorological Variables and Electricity Outages

The significance of the impact of weather on the electric grid has grown as climate change continues to increase the frequency and intensity of extreme weather events. In recent years (2021-2022) in particular, extreme winter weather has affected the grid in locations in the US rarely exposed to extreme low temperatures, snow and icing conditions. Here we analyze the correlation between cold weather meteorological variables and electricity outages during two large winter storm events, Uri (February 2021) and Landon (February 2022) using Random Forest machine learning and Pearson’s correlation coefficient. Our geographical focus across the two storms is the state of Texas. Extrapolation of the method to winter weather impacts over other years and additional locations is proposed.

Dumas, Melissa↗

DOE OSTI2024

Development and Validation of a Scoring System for Abnormalities in the Gopher Frog (Rana capito)

Headstarting efforts are thought to be critical in supplementing populations of the at-risk Gopher Frog (Rana capito); however, recent efforts have occasionally resulted in juveniles with developmental abnormalities. In response, we developed a scoring system to collect quantitative data on the presence and severity of these developmental abnormalities. Our objective was to describe and validate the abnormality scoring system so that it can be used by all Gopher Frog headstarting facilities. The scoring system covers five primary conditions encompassing commonly observed abnormalities. Two groups of participants with different levels of prior experience working with Gopher Frogs assigned scores to a set of images presented to them for each condition. We used intra-class correlation coefficients (ICC) to test the scoring system for inter- and intra-rater agreement as well as agreement with the benchmark standard (established by the authors). We found high ICC values for inter-rater agreement, intra-rater agreement, and agreement to the benchmark standard indicating either excellent or good reliability for all five conditions and for all raters when grouped together. These findings support the reliability and validity of the proposed developmental abnormality scoring system. Gopher Frog headstarting facilities can implement this scoring system to assist in tracking the frequency and severity of abnormalities observed in future headstarting efforts. We hope that by creating a reliable scoring system for Gopher Frogs, it can provide an overall framework and serve as a valuable resource to evaluate abnormalities across any amphibian species.

abnormality↗

DOE OSTI2025

Data-driven modeling of dynamic occupant thermostat override behavior for demand response applications

Buildings consume nearly 40% of global energy and produce similar emissions. Whiletechnological advances address efficiency, occupant behavior causes energy use variations up to 300% between identical buildings. This gap between predicted and actual building performance impacts building design, operations, and grid demand management programs. Through analyses of smart thermostat data from 1,400 single-occupant homes, the researchdemonstrates that occupants respond to 8°F thermostat setpoint changes within a median of 15 minutes, while 2°F changes trigger responses within a median of 30 minutes. This highlights an understudied temporal relationship between thermostat setbacks and response time of occupant behaviors. Models of such behavior dynamics are required to incorporate occupant impacts into building performance simulation. A key contribution of this dissertation is the Thermal Frustration Theory (TFT), which positsthat thermal discomfort driven behaviors are caused by the time-accumulation of discomfort, not simply a temperature deviation threshold or a delay from an initiating event. Using a dataset of 634 thermostats, each with 25+ manual setpoint changes, a comparative analysis of TFT and comfort zone and a delayed response theories demonstrated that personalized TFT models better predict when manual setpoint change occur. This was measured by the area under the curve statistical measure (AUC); all three models perform similarly by a Matthews Correlation Coefficient measure. Higher AUC performance is especially important for modeling occupant behavior in demand response programs where false negatives of rare occupant interactions could adversely affect grid stability. EnergyPlus based simulations were conducted with TFT-derived occupant models, demonstrating the ability to identify parameters of known TFT models from only data observable with smart thermostats, even under the presence of noise from routine overrides. Overall, the dissertation highlights that thermostat interactions are neither static,instantaneous, nor driven solely by the environment. Instead, temporal accumulation of discomfort and routine-based behavior play important roles. The methodology and results offer a pathway towards more accurate modeling of human-building interactions for policy assessment, building design, and demand response programs.

Sharma, Kunind [Northeastern University] (ORCID:00↗

DOE OSTI2025

Stream Temperature Responses to Summer Urban Rain Events Along the Savannah River

The hydrological urban heat island (HUHI) is a recent facet of the urban heat island (UHI), describing thermal enhancement of bodies of water in response to urbanization. Although the forefront of this work has been studied for large Metropolitan areas, this effect on developing cities is currently unknown. To locally measure and quantify HUHI effects in a developing city, we utilized Multi-Radar Multi Sensor (MRMS) radar-derived rainfall data to estimate rainfall in the Augusta Metropolitan Area (AMA) and the Savannah Metropolitan Area (SMA), two developing communities adjacent to the Savannah River. We analyzed mean temperature trends of two stream gauges measuring temperature at, and downstream of the central AMA and the SMA. Results show overall peak stream temperatures within AMA are reached within a 2hr timeframe with surge temperatures between 1-2 K, perhaps assisted by the presence of the Augusta Shoals further upstream. The SMA appears to have little HUHI effects due to urban greening and impacts to coastal meteorology. It is suggested that Moderate rainfall events (10-25mm of rainfall) have a significant correlation (~0.1 correlation coefficient with a significance level of 0.05) between stream temperature and rainfall, indicating potential warmer runoff input into the Savannah River. Although results seem promising, further research into moderate rainfall events are needed to determine the extent of the HUHI for AMA.

Wermter, Joseph E. [Savannah River National Labora↗

DOE OSTI2026

A Physics-Based Digital Twin for Wave Elevation and Seabed Moment Estimation of Offshore Monopiles: Preprint

In this work, we present a proof of concept of a physics-based digital twin for a monopile structure (with overhead inertia) subjected to wave loading. The digital twin is formulated using reduced-order models derived from first principles and combined with a Kalman filter for state estimation. The proposed framework estimates the monopile top motion, the wave elevation, and the section forces and moments along the pile using primarily acceleration measurements at the monopile top. Key innovations include the use of a hydrodynamic shape function to represent distributed wave loading in a compact and computationally efficient manner, and the introduction of a shaping filter to augment the state-space with wave kinematics. Synthetic measurement data are generated using OpenFAST and used as a reference to assess the performance of the digital twin. Results demonstrate that the wave elevation can be accurately reconstructed without direct sea-state measurements as long as the wave regime is inertia-dominated. Under the ideal tested conditions, the total hydrodynamic force and sea-bed bending moment are estimated with relative errors on the order of 1% and correlation coefficients exceeding 96%. Future work will evaluate the estimator's performance under operational uncertainties and more complex loading conditions.

17 WIND ENERGY↗