Search NASA⌕ Search

SEARCH · Search NASA

Results for “inversion and monitoring”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗
DOE OSTI2023

Microseismic Monitoring of Laboratory Hydraulic Fracturing Experiments in Granitic Rocks for Different Fracture Propagation Regimes

While hydraulic fracturing (HF) is a widely employed process, the underlying fracturing processes are still heavily contested. The attributes of the HF generated fracture network can exhibit substantial variation when dealing with specific HF propagation regimes encountered in the field. In this study, HF experiments were performed on true-triaxially loaded Barre granite cubes, with microseismic monitoring, to identify and characterize the fracturing mechanisms associated with different viscosity injection fluids. Utilizing fluids with high (oil/1450 cP) and low (water/1 cP) viscosity represented two key HF propagation regimes: viscosity- and toughness-dominated. The experiments conducted with oil involved higher breakdown pressures, larger fluid volumes, and slower fracture propagation speeds. Furthermore, the frequency–magnitude distribution (b value) for all experiments (1.9–2.3) was similar to those encountered for large-scale operations. Slightly larger b values were encountered during the initiation phase (2.4–2.7) relative to the fracture propagation and post-fracturing phases (1.9–2.2). Techniques such as polarity and moment-tensor inversion were utilized to characterize the source mechanisms. For the HF experiments with oil, tensile fractures were most dominant (92%) in the initiation phase compared to fracture propagation and post-fracturing phases (70–75%). Similar tensile fracturing dominancy was not observed with water, attributable to fluid permeation and leak-off. Regardless of the injection fluid or classification criteria employed, tensile fractures were the dominant type consistently, with fewer occurring in water experiments but the specific ratio of crack types varied with different source mechanism criteria employed.

58 GEOSCIENCES↗
DOE OSTI2020

Novel Perovskite Semiconductors for the Detection of Special Nuclear Material

Special nuclear material detection using portable, highly efficient detectors is of paramount importance to national security and environmental monitoring. Current neutron detectors, including gas proportional counters and scintillators, are not suitable for discrete or portable monitoring due to size, weight, and power limitations. For example, 3 He gas proportional detectors have significant drawbacks, including long absorption lengths, high voltage requirements, and dependence on the availability of 3 He. The goal of this project is to synthesize a novel, boron-based perovskite semiconductor crystals to create a more efficient and more portable neutron detector. This innovative material will lead to a significant decrease in power and size requirements with higher operating temperatures than other traditional solid-state detectors. Directly incorporating 10 B eliminates the need for a conversion layer and leads to high neutron absorption coefficients. In this seedling project, we used inverse temperature crystallization with boron precursors to create CH 3 NH 3 PbI 3-x (BF 4 ) x . The crystals were analyzed with X-ray diffraction, but the results were inconclusive. Future work will focus on tailoring the solution growth and selecting appropriate solvents and precursors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Dynamic changes in evapotranspiration, canopy photosynthesis and expansive growth in open field under rapid fluctuating radiation

Abstract In open fields, plants experience dynamic changes in environment, particularly radiation, temperature, wind, and humidity but their short-term responses have not been adequately characterized under natural conditions. In this study, we assessed causal effects of rapid radiation fluctuation on seven plant parameters in open fields of cotton and sweet corn canopies. The parameters are evapotranspiration ( ET ), canopy photosynthesis indicated by CO2 flux from air to canopy ( FCO2 ), sensible heat energy flux ( H ), canopy conductance for water vapor ( gcw ), canopy surface temperature ( Ts ), shoot and leaf elongation rate, and stem diameter change. The energy and CO2 fluxes were measured with Bowen ratio/energy balance/CO2 gradient (BREB+) technique, using averaging time of 5 min. Shoot + leaf elongation and stem diameter change were monitored with position transducers using averaging time of 1 min. All parameters were all found to respond to change in radiation and transpiration within minutes or sooner. While radiation effects on canopy gas exchanges are expected, illuminating are the indirect but immediate effects on shoot + leaf growth and stem diameter change through radiation effects on transpiration and plant water status. A novel finding is that gcw also responded within minutes or sooner to radiation fluctuations and that FCO2 was related almost linearly to gcw . Results are discussed in terms of soil-plant-atmosphere continuum, and interpreted in terms of dynamic interactions between transpiration and plant water status. The clear inverse relationship between ET and elongation rate or stem diameter changes provides additional evidence supporting the validity of 5-min averaging for the BREB + technique.

Hsiao, Theodore C.↗
DOE OSTI2021

A deep learning-accelerated data assimilation and forecasting workflow for commercial-scale geologic carbon storage

Fast assimilation of monitoring data to update forecasts of pressure buildup and carbon dioxide (CO 2 ) plume migration under geologic uncertainties is a challenging problem in geologic carbon storage. The high computational cost of data assimilation with a high-dimensional parameter space impedes fast decision-making for commercial-scale reservoir management. We propose to leverage physical understandings of porous medium flow behavior with deep learning techniques to develop a fast data assimilation-reservoir response forecasting workflow. Applying an Ensemble Smoother Multiple Data Assimilation (ES-MDA) framework, the workflow updates geologic properties and predicts reservoir performance with quantified uncertainty from pressure history and CO 2 plumes interpreted through seismic inversion. As the most computationally expensive component in such a workflow is reservoir simulation, we developed surrogate models to predict dynamic pressure and CO 2 plume extents under multi-well injection. The surrogate models employ deep convolutional neural networks, specifically, a wide residual network and a residual U-Net. The workflow is validated against a flat threedimensional reservoir model representative of a clastic shelf depositional environment. Intelligent treatments are applied to bridge between quantities in a true-3D reservoir model and those in a single-layer reservoir model. The workflow can complete history matching and reservoir forecasting with uncertainty quantification in less than one hour on a mainstream personal workstation.

25 ENERGY STORAGE↗
DOE OSTI2020

Investigating effects of soil chemicals on density of small mammal bioindicators using spatial capture-recapture models

Monitoring the ecological impacts of environmental pollution and the effectiveness of remediation efforts requires identifying relationships between contaminants and the disruption of biological processes in populations, communities, or ecosystems. Wildlife are useful bioindicators, but traditional comparative experimental approaches rely on a staunch and typically unverifiable assumption that, in the absence of contaminants, reference and contaminated sites would support the same densities of bioindicators, thereby inferring direct causation from indirect data. We demonstrate the utility of spatial capture-recapture (SCR) models for overcoming these issues, testing if community density of common small mammal bioindicators was directly influenced by soil chemical concentrations. By modeling density as an inhomogeneous Poisson point process, we found evidence for an inverse spatial relationship between Peromyscus density and soil mercury concentrations, but not other chemicals, such as polychlorinated biphenyls, at a site formerly occupied by a nuclear reactor. Although the coefficient point estimate supported Peromyscus density being lower where mercury concentrations were higher (β = –0.44), the 95% confidence interval overlapped zero, suggesting no effect was also compatible with our data. Estimated density from the most parsimonious model (2.88 mice/ha; 95% CI = 1.63–5.08), which did not support a density-chemical relationship, was within the range of reported densities for Peromyscus that did not inhabit contaminated sites elsewhere. Environmental pollution remains a global threat to biodiversity and ecosystem and human health, and our study provides an illustrative example of the utility of SCR models for investigating the effects that chemicals may have on wildlife bioindicator populations and communities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Utah FORGE: Phase 3 Magnetotelluric (MT) Data

New high-quality tensor MT data at 122 sites, including the vertical magnetic field and utilizing ultra-remote referencing, have been acquired over the Utah FORGE project area. The results will be used to delineate the densities of faults and fractures in crystalline basement rocks, to define the heat sources, and to derive baseline 3D resistivity structure for later MT monitoring of temporal changes in resistivity structure following well stimulation in the EGS reservoir. There are three files here related to Utah FORGE magnetotelluric (MT) data acquisition and processing. The FORGE MT EDIs zip file contains the observed MT responses in industry-standard EDI format. For each site, there is an EDI response file that utilized a local independent reference for noise cancellation, and a file that utilized a distant reference cancelling noise associated with the DC transmission line of the Delta IPP passing down the west side of Milford Valley. These two site files could be merged as appropriate. The FORGE Model Cell Center file contains the model volume of the 3D Forge MT inversion for characterizing the resistivity structure in the project area. It was derived using finite element inversion methodology described in Wannamaker et al, in the attached FORGE Phase3 Geophysics paper, from the MT observation EDI files. It is ASCII format (.dat) and entries are defined at top of the file in a simple x-y-z-Rho listing in UTM coordinates. The element layers drape the topography so the Rho value layers are not purely horizontal slices. This greatly simplifies the listing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Rapid Optimization of Total Variation with Applications in Imaging, Additive Manufacturing, and Qualification

Total Variation optimization penalizes the gradient of a control variable or state. While this work focuses on image processing in particular, it has also found applications in inverse problems and topology optimization. In image processing, the goal is to maintain faithfulness to the original image while denoising and/or deblurring. Additionally, bilevel optimization over the spatially varying regularization weights can illuminate interfaces such as damage regions and other anomalies. We will address two fundamental challenges with TV-optimization: (i) the typical slow convergence of existing TV-optimization methods, and (ii) the selection of spatially varying TV parameters to promote interface detection. Additionally, we will apply such techniques to image data collected in additive manufacturing. In said context, stochasticity in build events induces flaws in the manufactured piece, compromising the integrity of said part. There is a critical need for in-situ monitoring to spot anomalies once they form, and in this setting we apply our total variation and hyperparameter solvers. We will develop a customized algorithm based on for extreme-scale TV-optimization that achieves super-linear or quadratic-convergence, a critical property for real-time, image-by-image analysis. A worst-case outcome is a preprocessing step that enhances image quality in-situ, specifically for out-of-focus and noisy images.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Acoustoelastic characterization of plates using zero group velocity Lamb modes

Acoustoelasticity, a characteristic of material anharmonicity, gives rise to a link between wave propagation velocity and the stress state in materials. Ultrasonic techniques to monitor this coupling, particularly with high sensitivity and in a noncontact manner, can have widespread application both in the quantification of applied and residual stress and in the characterization of nonlinear material behavior through measurement of higher order elastic constants. Here, we use a laser ultrasonic technique to excite and detect zero group velocity (ZGV) Lamb wave resonances in aluminum plates under uniaxial loading. A laser line source is used to excite these resonances at different orientations with respect to the applied load, and the signals are detected using an interferometer. The effects of stress and source orientation on ZGV resonance frequencies are validated using the theory of acoustoelastic Lamb wave propagation. In addition, a model-based inversion technique is used to extract Murnaghan's third-order elastic constants from measurements of the stress dependence of the first two ZGV modes generated parallel and perpendicular to the applied load. Laser generation and detection of ZGV resonances is shown to be an effective and powerful approach for the noncontact and nondestructive acoustoelastic characterization of elastic waveguides.

36 MATERIALS SCIENCE↗
DOE OSTI2024

AE33 Aethalometer Records at RMBL [AMF2] from April 2022 to October 2023

The dataset presents atmospheric mass loadings of black carbon (BC) and brown carbon (BrC), during the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign that took place from 2021 – 2023 in Gothic, CO. A Magee Scientific aethalometer (model AE33) was deployed at the Rocky Mountain Biological Laboratory, 200 m north of AMF2 to monitor mass loadings of optically absorbing aerosols. The instrument inlet was equipped with a BGI SCC 1.829 ambient cyclone with a particle cutoff of 2.5 μm and placed 5 m above ground level. The dataset reports hour-resolved BC mass loadings at 370, 470, 520, 590, 660, 880, and 950 (BC370, BC470, BC 520, BC590, BC660, BC880, BC950 in ng/m3) reported by the instrument built-in algorithm, which calculates mass loadings from the rate of change of the attenuation of light transmitted through the aerosol-laden filter. BC880 is the operationally defined reference of the BC concentration. Hour-resolved absorbance at each wavelength were calculated from their respective mass absorption cross-section in inverse megameters (Mm-1). The absorption Angstrom exponent (AAE) for BrC and BC were calculated for each hour. BrC AAE (AAE_BrC) was calculated from the slope of Log(wavelength) versus Log(absorbance) based on the measurements at 370, 470, and 520 nm. BC AAE (AAE_BrC) was calculated from the slope of Log(wavelength) versus Log(absorbance) based on the measurements at 590, 660, 880, and 950 nm. The aethalometer record spans time periods of 4/5/2022 – 10/3/2022 and 12/8/2022 – 10/9/2023. The ambient CO2 concentration recorded by a Vaisala CARBOCAP carbon dioxide probe GMP343 is also reported in ppm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Toward the implementation of self-powered, wireless, real-time reactor power sensing

The successful development of real-time, self-powered, wireless reactor power sensing technology would enable key operating parameters of every fuel assembly in a commercial reactor core to be continuously monitored, and potentially be utilized for applications such as closed-loop reactor control. This paper discusses the expansion of successful work on the application of Self Powered Neutron Detectors (SPND) for reactor power monitoring and the wireless transmission of these detector signals using advanced vacuum micro-electronic (VME) transmitter technology developed by Westinghouse in collaboration with the Penn State University Radiation Science and Engineering Center. Results demonstrate that this technology may be successfully applied via long-term testing in the Penn State TRIGA reactor. SPND are inherently slower than, for example, ion chambers due to their dependence on radioactive decay to generate a signal. Additional research developed delay compensation algorithms to accelerate the response of these detectors for use in near-real-time closed-loop reactor control. Initial results indicate that practical signal acceleration can be obtained digitally or with an equivalent analog circuit. Future research will apply an inverse model of the sensor’s balance differential equations to more-effectively accelerate the SPND signal. Use of inverse models will also prove to be adaptive in the sense that time-dependent burn-up of these detectors, due to neutron activation of the detector material, can be modeled and applied in the acceleration algorithm.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Machine learning-driven design and self-sensing capabilities of automotive bumper lattices for adaptive impact response

We present a novel approach to design an automotive bumper energy absorber using carbon fiber reinforced polymer composites, optimized to meet conflicting performance requirements for two distinct impact scenarios. The design must satisfy both a low-speed (2.5 mph) pendulum intrusion test, simulating vehicle-to-vehicle collisions, and a high-speed (25 mph) leg flexion test, replicating pedestrian impacts. These tests demand opposing deformation characteristics: high flexibility (deformation < 85 mm) for the former and high stiffness (deformation < 22 mm) for the latter. To address these contradictory requirements, we developed a machine learning (ML) framework for inverse optimization of lattice designs and material selection. Unlike traditional iterative design processes, our ML model directly outputs optimal design parameters and material choices based on target performance inputs. The energy absorber was fabricated using advanced additive manufacturing techniques, including extrusion deposition and digital light processing. The integration of carbon fibers provides multifunctionality to the bumper structure, enabling self-sensing capabilities through changes in electrical resistivity under compression. This electrical response demonstrates high repeatability under multiple cycles at 2% compression and exhibits distinct signatures during crack formation under high deformation. This research offers adaptive performance through innovative design methodologies and smart material integration. The approach has potential applications in various fields requiring adaptive energy absorption and real-time structural health monitoring.

Chawla, Komal [ORNL] (ORCID:0000000190327565)↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2021

EXPRES. II. Searching for Planets around Active Stars: A Case Study of HD 101501

By controlling instrumental errors to below 10 cm s{sup −1}, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian processes (GP) have become a standard tool for modeling correlated noise in radial velocity data sets. While GPs are constrained and motivated by physical properties of the star, in some cases they are still flexible enough to absorb unresolved Keplerian signals. We apply GP regression to EXPRES radial velocity measurements of the 3.5 Gyr old chromospherically active Sun-like star, HD 101501. We obtain tight constraints on the stellar rotation period and the evolution of spot distributions using 28 seasons of ground-based photometry, as well as recent Transiting Exoplanet Survey Satellite data. Light-curve inversion was carried out on both photometry data sets to reveal the spot distribution and spot evolution timescales on the star. We find that the >5 m s{sup −1} rms radial velocity variations in HD 101501 are well modeled with a GP stellar activity model without planets, yielding a residual rms scatter of 45 cm s{sup −1}. We carry out simulations, injecting and recovering signals with the GP framework, to demonstrate that high-cadence observations are required to use GPs most efficiently to detect low-mass planets around active stars like HD 101501. Sparse sampling prevents GPs from learning the correlated noise structure and can allow it to absorb prospective Keplerian signals. We quantify the moderate to high-cadence monitoring that provides the necessary information to disentangle photospheric features using GPs and to detect planets around active stars.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Dynamic concurrent learning method to neutralize cyber attacks and faults for industrial asset monitoring nodes

Input signals may be received from monitoring nodes of the industrial asset, each input signal comprising time series data representing current operation. A neutralization engine may transform the input signals into feature vectors in feature space, each feature vector being associated with one of a plurality of overlapping batches of received input signals. A dynamic decision boundary may be generated based on the set of feature vectors, and an abnormal state of the asset may be detected based on the set of feature vectors and a predetermined static decision boundary. An estimated neutralized value for each abnormal feature value may be calculated based on the dynamic decision boundary and the static decision boundary such that a future set of feature vectors will be moved with respect to the static decision boundary. An inverse transform of each estimated neutralized value may be performed to generate neutralized signals comprising time series data that are output.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantitatively Monitoring Bubble-Flow at a Seep Site Offshore Oregon: Field Trials and Methodological Advances for Parallel Optical and Hydroacoustical Measurements

Two lander-based devices, the Bubble-Box and GasQuant-II, were used to investigate the spatial and temporal variability and total gas flow rates of a seep area offshore Oregon, United States. The Bubble-Box is a stereo camera–equipped lander that records bubbles inside a rising corridor with 80 Hz, allowing for automated image analyses of bubble size distributions and rising speeds. GasQuant is a hydroacoustic lander using a horizontally oriented multibeam swath to record the backscatter intensity of bubble streams passing the swath plain. The experimental set up at the Astoria Canyon site at a water depth of about 500 m aimed at calibrating the hydroacoustic GasQuant data with the visual Bubble-Box data for a spatial and temporal flow rate quantification of the site. For about 90 h in total, both systems were deployed simultaneously and pressure and temperature data were recorded using a CTD as well. Detailed image analyses show a Gaussian-like bubble size distribution of bubbles with a radius of 0.6–6 mm (mean 2.5 mm, std. dev. 0.25 mm); this is very similar to other measurements reported in the literature. Rising speeds ranged from 15 to 37 cm/s between 1- and 5-mm bubble sizes and are thus, in parts, slightly faster than reported elsewhere. Bubble sizes and calculated flow rates are rather constant over time at the two monitored bubble streams. Flow rates of these individual bubble streams are in the range of 544–1,278 mm 3 /s. One Bubble-Box data set was used to calibrate the acoustic backscatter response of the GasQuant data, enabling us to calculate a flow rate of the ensonified seep area (~1,700 m 2 ) that ranged from 4.98 to 8.33 L/min (5.38 × 10 6 to 9.01 × 10 6 CH 4 mol/year). Such flow rates are common for seep areas of similar size, and as such, this location is classified as a normally active seep area. For deriving these acoustically based flow rates, the detailed data pre-processing considered echogram gridding methods of the swath data and bubble responses at the respective water depth. The described method uses the inverse gas flow quantification approach and gives an in-depth example of the benefits of using acoustic and optical methods in tandem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interannual variability and trends of summertime PM2.5-based air quality in the Intermountain West

Abstract Summertime air quality is a growing public health concern in the populated region of Northern Utah. Whereas winter air pollution is highly linked with local atmospheric temperature inversions associated with upper atmospheric high-pressure and radiational cooling in valleys, the relationship between climate factors and the frequency of poor air quality during summer is still unknown. Analyzing the last 20 years of data, we demonstrated that summertime unhealthy days (as defined by PM2.5 air quality index level) in Northern Utah highly correlate with the number of dry-hot days, wildfire size, and an upper atmospheric ridge over the Northwestern United States. The persistent atmospheric ridge enhances lightning-caused fire burned areas in northwestern states and then transports the wildfire smoke toward Northern Utah. Similarly, climate model simulations confirm observational findings, such as an increasing trend of the upper atmospheric ridge and summertime dry days in the northwestern states. Such metrics developed in this study could be used to establish longer-term monitoring and seasonal forecasting for air quality and its compounding factors, which is currently limited to forecasting products for only several days.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A scalable variational method for estimating the latent infection-rate field of an outbreak

In this paper, we explore whether the infection-rate of a disease can serve as a robust monitoring variable in epidemiological surveillance algorithms. The infection-rate is dependent on population mixing patterns that do not vary erratically day-to-day; in contrast, daily case-counts used in contemporary surveillance algorithms are corrupted by reporting errors. The technical challenge lies in estimating the latent infection-rate from case-counts. Here we devise a Bayesian method to estimate the infection-rate across multiple adjoining areal units, and then use it, via an anomaly detector, to discern a change in epidemiological dynamics. We extend an existing model for estimating the infection-rate in an areal unit by incorporating a Markov random field model, so that we may estimate infection-rates across multiple areal units, while preserving spatial correlations observed in the epidemiological dynamics. To carry out the high-dimensional Bayesian inverse problem, we develop an implementation of mean-field variational inference specific to the infection model and integrate it with the random field model to incorporate correlations across counties. The method is tested on estimating the COVID-19 infection-rates across all 33 counties in New Mexico using data from the summer of 2020, and then employing them to detect the arrival of the Fall 2020 COVID-19 wave. We perform the detection using a temporal algorithm that is applied county-by-county. We also show how the infection-rate field can be used to cluster counties with similar epidemiological dynamics.

60 APPLIED LIFE SCIENCES↗