SEARCH · Search NASA

Results for “spatiotemporal data”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Evapotranspiration Partitioning Using Flux Tower Data in a Semi-Arid Ecosystem

Information about evapotranspiration (ET) and its components, that is, evaporation and transpiration, is crucial for a wide range of water and ecosystem management applications. However, partitioning ET into its two components is often challenging because of their spatiotemporal variabilities and lack of process understanding. This study developed a machine learning (ML) framework to shed light on ET processes and assess the relative importance of different drivers by incorporating hydrometeorology and biomass productivity variables. The Shapley Additive Explanations (SHAP) approach was applied to enhance explainability and rank the importance of ET drivers and their components. A total of 62 variables covering hydrometeorological and biomass productivity dimensions were considered from the Reynolds Creek Critical Zone Observatory (CZO) station in Idaho. The variable importance assessment identified the leading drivers individually for evaporation, transpiration and ET (soil water content for evaporation, vapour pressure deficit for transpiration and soil water content for ET). The results further highlighted the value of combining hydrometeorological and biomass productivity variables to achieve better predictability of ET processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Multichannel Doppler fiber-imager spectrometer for spatiotemporal characterization of high-intensity laser-driven plasmas

Here, we demonstrate a flexible multichannel fiber-based imaging Doppler spectrometer to characterize plasmas in high intensity (≥1 × 10 18 W/cm 2 ) laser-plasma experiments at high repetition rates. This instrument collects data from ×21 different plasma locations combining optical fibers and a single imaging spectrometer. This diagnostic maps the plasma velocity evolution as a function of time with sub-pico-second resolution. Experimental results showing 2D velocity measurements of plasma with 20 μm spatial resolution are presented. Intensities of the order of 10 18 W/cm 2 were used to generate a plasma, while a much less intense, frequency doubled (400 nm), probe beam (10 11 W/cm 2 ) was used to measure the Doppler shift from the plasma critical surface. The instrument can be scaled to a larger number of channels (e.g., 100) still using a single spectrometer.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Influence of soil heterogeneity on soybean plant development and crop yield evaluated using time-series of UAV and ground-based geophysical imagery

Understanding the interactions among agricultural processes, soil, and plants is necessary for optimizing crop yield and productivity. This study focuses on developing effective monitoring and analysis methodologies that estimate key soil and plant properties. These methodologies include data acquisition and processing approaches that use unmanned aerial vehicles (UAVs) and surface geophysical techniques. In particular, we applied these approaches to a soybean farm in Arkansas to characterize the soil–plant coupled spatial and temporal heterogeneity, as well as to identify key environmental factors that influence plant growth and yield. UAV-based multitemporal acquisition of high-resolution RGB (red–green–blue) imagery and direct measurements were used to monitor plant height and photosynthetic activity. We present an algorithm that efficiently exploits the high-resolution UAV images to estimate plant spatial abundance and plant vigor throughout the growing season. Such plant characterization is extremely important for the identification of anomalous areas, providing easily interpretable information that can be used to guide near-real-time farming decisions. Additionally, high-resolution multitemporal surface geophysical measurements of apparent soil electrical conductivity were used to estimate the spatial heterogeneity of soil texture. By integrating the multiscale multitype soil and plant datasets, we identified the spatiotemporal co-variance between soil properties and plant development and yield. Our novel approach for early season monitoring of plant spatial abundance identified areas of low productivity controlled by soil clay content, while temporal analysis of geophysical data showed the impact of soil moisture and irrigation practice (controlled by topography) on plant dynamics. Our study demonstrates the effective coupling of UAV data products with geophysical data to extract critical information for farm management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Ultra-Short-Term Spatiotemporal Forecasting of Renewable Resources: An Attention Temporal Convolutional Network Based Approach

The rapid increase in the penetration of renewable energy resources characterized by high variability and uncertainty is bringing new challenges to the power system operation. To ensure the efficient and reliable operation of electric grid, an accurate and general short-term forecasting algorithm with interpretability is desired. Moreover, the extensive off-site information provided by the proliferation of new renewable plants stimulates the interests in the spatiotemporal forecasting. In this paper, an attention temporal convolutional network, which is built on stacked dilated causal convolutional networks and attention mechanisms, is proposed to perform the ultra-short-term spatiotemporal forecasting of renewable resources. Compared with the existing spatiotemporal forecasting methods, the presented model needs no domain knowledge and can be applied to different forecasting tasks such as solar generation and wind speed forecasting. Here, the attention mechanism improves the interpretability. The algorithm can be used to produce both point and probabilistic forecasts. Numerical results on the data sets from National Renewable Energy Laboratory show superior performance over five baselines, in terms of skill scores. Compared with the baselines, the average improvements of accuracy introduced by the proposed method for the point and probabilistic forecasting are 15.08% and 15.85%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Planar laser induced fluorescence mapping of a carbon laser produced plasma

We present measurements of ion velocity distribution profiles obtained by laser induced fluorescence (LIF) on an explosive laser produced plasma. The spatiotemporal evolution of the resulting carbon ion velocity distribution was mapped by scanning through the Doppler-shifted absorption wavelengths using a tunable, diode-pumped laser. The acquisition of these data was facilitated by the high repetition rate capability of the ablation laser (1 Hz), which allowed for the accumulation of thousands of laser shots in short experimental times. By varying the intensity of the LIF beam, we were able to explore the effects of fluorescence power against the laser irradiance in the context of evaluating the saturation vs the non-saturation regime. The small size of the LIF beam led to high spatial resolution of the measurement compared to other ion velocity distribution measurement techniques, while the fast-gate operation mode of the camera detector enabled the measurement of the relevant electron transitions.

Instruments & Instrumentation↗

DOE OSTI2023

Assessing residential PM 2.5 concentrations and infiltration factors with high spatiotemporal resolution using crowdsourced sensors

Building conditions, outdoor climate, and human behavior influence residential concentrations of fine particulate matter (PM 2.5 ). To study PM 2.5 spatiotemporal variability in residences, we acquired paired indoor and outdoor PM 2.5 measurements at 3,977 residences across the United States totaling >10,000 monitor-years of time-resolved data (10-min resolution) from the PurpleAir network. Time-series analysis and statistical modeling apportioned residential PM 2.5 concentrations to outdoor sources (median residential contribution = 52% of total, coefficient of variation = 69%), episodic indoor emission events such as cooking (28%, CV = 210%) and persistent indoor sources (20%, CV = 112%). Residences in the temperate marine climate zone experienced higher infiltration factors, consistent with expectations for more time with open windows in milder climates. Likewise, for all climate zones, infiltration factors were highest in summer and lowest in winter, decreasing by approximately half in most climate zones. Large outdoor–indoor temperature differences were associated with lower infiltration factors, suggesting particle losses from active filtration occurred during heating and cooling. Absolute contributions from both outdoor and indoor sources increased during wildfire events. Infiltration factors decreased during periods of high outdoor PM 2.5 , such as during wildfires, reducing potential exposures from outdoor-origin particles but increasing potential exposures to indoor-origin particles. Time-of-day analysis reveals that episodic emission events are most frequent during mealtimes as well as on holidays (Thanksgiving and Christmas), indicating that cooking-related activities are a strong episodic emission source of indoor PM 2.5 in monitored residences.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Using Image Processing Techniques to Identify and Quantify Spatiotemporal Carbon Cycle Extremes

Rising atmospheric carbon dioxide due to human activities through fossil fuel emissions and land use changes have increased climate extremes such as heat waves and droughts that have led to and are expected to increase the occurrence of carbon cycle extremes. Carbon cycle extremes represent large anomalies in the carbon cycle that are associated with gains or losses in carbon uptake. Carbon cycle extremes could be continuous in space and time and cross political boundaries. Here, we present a methodology to identify large spatiotemporal extremes (STEs) in the terrestrial carbon cycle using image processing tools for feature detection. We characterized the STE events based on neighborhood structures that are three-dimensional adjacency matrices for the detection of spatiotemporal manifolds of carbon cycle extremes. We found that the area affected and carbon loss during negative carbon cycle extremes were consistent with continuous neighborhood structures. In the gross primary production data we used, 100 carbon cycle STEs accounted for more than 75% of all the negative carbon cycle extremes. This paper presents a comparative analysis of the magnitude of carbon cycle STEs and attribution of those STEs to climate drivers as a function of neighborhood structures for two observational datasets and an Earth system model simulation.

Sharma, Bharat↗

DOE OSTI2025

Brief communication: Decadal changes in topography, surface water and subsurface structure across an Arctic coastal tundra site

Abstract. In ice-rich polygonal tundra, spatiotemporal heterogeneity in ground-ice melt reshapes topography, impacting infrastructure, water and carbon cycles. This study evaluates changes in topography and subsurface structure at a coastal Arctic site by comparing data from two surveys conducted a decade apart. Each survey includes electrical resistivity tomography, active layer thickness, photogrammetry, and topographic data. Results reveal subsidence and decrease in permafrost table elevation with varying intensity and spatial distribution across polygons, alongside diverse thermal-hydrological responses, such as thermokarst pool formation in high-centered-polygons and more even subsidence in flat-centered-polygons. The study also underscores the value and limitations of sporadic surveys.

Bachman, Jonathan A↗

DOE OSTI2020

What can we learn from N 2 O isotope data? – Analytics, processes and modelling

The isotopic composition of nitrous oxide (N 2 O) provides useful information for evaluating N 2 O sources and budgets. Due to the co–occurrence of multiple N 2 O transformation pathways, it is, however, challenging to use isotopic information to quantify the contribution of distinct processes across variable spatiotemporal scales. Here, we present an overview of recent progress in N 2 O isotopic studies and provide suggestions for future research, mainly focusing on: analytical techniques; production and consumption processes; and interpretation and modelling approaches. Comparing isotope–ratio mass spectrometry (IRMS) with laser absorption spectroscopy (LAS), we conclude that IRMS is a precise technique for laboratory analysis of N 2 O isotopes, while LAS is more suitable for in situ/inline studies and offers advantages for site–specific analyses. When reviewing the link between the N 2 O isotopic composition and underlying mechanisms/processes, we find that, at the molecular scale, the specific enzymes and mechanisms involved determine isotopic fractionation effects. In contrast, at plot–to–global scales, mixing of N 2 O derived from different processes and their isotopic variability must be considered. We also find that dual isotope plots are effective for semi–quantitative attribution of co–occurring N 2 O production and reduction processes. More recently, process–based N 2 O isotopic models have been developed for natural abundance and 15 N–tracing studies, and have been shown to be effective, particularly for data with adequate temporal resolution. Here, despite the significant progress made over the last decade, there is still great need and potential for future work, including development of analytical techniques, reference materials and inter–laboratory comparisons, further exploration of N 2 O formation and destruction mechanisms, more observations across scales, and design and validation of interpretation and modelling approaches. Synthesizing all these efforts, we are confident that the N 2 O isotope community will continue to advance our understanding of N 2 O transformation processes in all spheres of the Earth, and in turn to gain improved constraints on regional and global budgets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗

DOE OSTI2020

Relative Roles of Low‐Level Wind Speed and Moisture in the Diurnal Cycle of Rainfall Over a Tropical Island Under Monsoonal Flows

Abstract Using satellite observations and reanalysis data, this study investigates the relative roles of low‐level wind speed and moisture content in the rainfall diurnal cycle over a tropical island influenced by monsoonal winds. Results show that the spatiotemporal characteristics of rainfall over the island are evidently influenced by the two factors. Specifically, the spatial distribution of rainfall is primarily determined by the low‐level wind speeds, while the rainfall diurnal cycle is influenced by both factors. Weaker low‐level wind speed leads to stronger and wider inland convergence and afternoon precipitation, which are mainly produced by the inland penetrating sea breeze fronts. Nighttime rainfall over the island is closely related to the inertial oscillation. Stronger nocturnal rainfall can be found with higher low‐level moisture content and stronger monsoonal winds. The multiscale interactions between the large‐scale circulation, land–sea breezes, and inertial oscillation play essential roles in modulating the rainfall over this tropical island.

Zhu, Lei↗

DOE OSTI2025

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

DOE OSTI2021

Bragg Coherent Diffraction Imaging for In Situ Studies in Electrocatalysis

Electrocatalysis is at the heart of a broad range of physicochemical applications that play an important role in the present and future of a sustainable economy. Among the myriad of different electrocatalysts used in this field, nanomaterials are of ubiquitous importance. An increased surface area/volume ratio compared to bulk makes nanoscale catalysts the preferred choice to perform electrocatalytic reactions. Bragg coherent diffraction imaging (BCDI) was introduced in 2006 and since has been applied to obtain 3D images of crystalline nanomaterials. BCDI provides information about the displacement field, which is directly related to strain. Lattice strain in the catalysts impacts their electronic configuration and, consequently, their binding energy with reaction intermediates. Even though there have been significant improvements since its birth, the fact that the experiments can only be performed at synchrotron facilities and its relatively low resolution to date (~10 nm spatial resolution) have prevented the popularization of this technique. Herein, we will briefly describe the fundamentals of the technique, including the electrocatalysis relevant information that we can extract from it. Subsequently, we review some of the computational experiments that complement the BCDI data for enhanced information extraction and improved understanding of the underlying nanoscale electrocatalytic processes. We next highlight success stories of BCDI applied to different electrochemical systems and in heterogeneous catalysis to show how the technique can contribute to future studies in electrocatalysis. Finally, we outline current challenges in spatiotemporal resolution limits of BCDI and provide our perspectives on recent developments in synchrotron facilities as well as the role of machine learning and artificial intelligence in addressing them.

bragg coherent diffraction imaging↗

DOE OSTI2021

Leveraging observed soil heterotrophic respiration fluxes as a novel constraint on global-scale models

Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.

Jian, Jinshi↗

DOE OSTI2024

Observed Variability in Convective Cell Characteristics and Near-Storm Environments across the Sea- and Bay-Breeze Fronts in Southeast Texas

Abstract During the DOE Atmospheric Radiation Measurement (ARM) Tracking Aerosol Convection Interactions Experiment (TRACER) IOP spanning June–September 2022, two fixed ARM sites and a mobile team concurrently sampled the airmass heterogeneity across sea- and bay-breeze fronts around the greater Houston metropolitan region. Here, we quantify the spatiotemporal variability between maritime (coastal/bay side of breeze fronts) and continental (inland side of breeze fronts) air masses over 15 IOP days characterized by strong sea-breeze forcing. We analyze environmental profile data from 177 radiosondes and use S- and C-band radar data to track and quantify the variability in attributes of more than 2300 shallow and transitioning cells across different air masses. The composite analysis of environmental profiles indicates that during the early afternoon, the sea-breeze maritime air mass exhibits lower convective available potential energy (CAPE) than the bay-breeze maritime air mass. As the sea breeze advances inland with time, CAPE within the maritime air mass exceeds that of the continental air mass to the north of the breeze fronts. In general, maritime cells have a larger mean composite reflectivity and cell widths than continental cells; however, the response varies between shallow and transitioning cells. Mean composite 20-dB Z echo-top heights, however, are similar across air masses for both shallow and transitioning cells. The continental and maritime inflow air mass for transitioning cells has significantly different mean values for mixed-layer entrainment CAPE, lifted condensation level, level of free condensation, boundary layer depth, and diluted equilibrium level. For shallow cells, only total precipitable water shows a significant difference. Significance Statement The greater Houston metropolitan area is a natural laboratory for understanding the individual impacts of background meteorology and aerosols on convective clouds. Due to its proximity to the Gulf Coast and Galveston Bay, the Houston region experiences a diurnal precipitation cycle in the summer, driven by convection triggered from sea- and bay-breeze fronts. These fronts act as a boundary between air masses with distinct thermodynamic and environmental characteristics. Convergence along these fronts and interactions between storm outflow and the fronts facilitate convection initiation in different mesoscale air masses. This study quantifies the heterogeneity among these air masses while investigating their influence on cloud microphysics. We find that the effect of airmass heterogeneity is more pronounced for the bulk microphysical properties in shallow clouds.

Sharma, Milind↗

DOE OSTI2022

Progress Towards a Predictive Eagle Behavior and Risk Modeling Framework: Overview and Recent Validation Efforts

This presentation summarizes progress to date of the U.S. Department of Energy project, "Development of a computational framework for modeling golden eagles (Aquila chrysaetos) near wind farms," which focuses on stochastic behavioral modeling of soaring raptors across landscape, facility, and turbine spatiotemporal scales. This publicly available, open-source modeling framework includes behavioral models based on three different underlying principles: energy minimization at landscape scale, behavioral heuristics at landscape-facility scale, and data-driven behaviors at the facility-micro-scale. We will briefly overview the key advancements in the behavioral modeling state of the art, which leverages multiple high-resolution telemetry data sources combined with high-fidelity atmospheric flow modeling insights. We then present preliminary results from a validation study in Altamont, California. This new study involves a novel application of the Stochastic Soaring Raptor Simulator (SSRS), in a new geographic locale, to understand facility scale eagle movement patterns over time scales representative of a wind project's lifetime. For this desktop analysis (that does not depend on any high-performance computing resources), SSRS simultaneously considers a variety of wind conditions and eagle approach vectors toward a project site of interest. This work demonstrates the integration of publicly available landscape-scale atmospheric datasets, our recently improved engineering updraft models (see presentation from Thedin et al.), and our energy minimization behavioral models within the SSRS framework. While we only present results from a single behavioral model, the integration of these three modeling components forms the foundation for our more sophisticated behavioral models (see presentations from Brandes et al., Sandhu et al.) that are under active development. Results are presented in the form of presence maps, which may be applied to estimate risk to wildlife, augment ground survey data, inform wind-plant operations, or incorporated into wind-plant designs.

agent-based modeling↗

DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗

DOE OSTI2025

Identifying spatiotemporal patterns in opioid vulnerability: investigating the links between disability, prescription opioids and opioid-related mortality

Background: The opioid crisis remains one of the most daunting and complex public health problems in the United States. This study investigates the national epidemic by analyzing vulnerability profiles of three key factors: opioid-related mortality rates, opioid prescription dispensing rates, and disability rank ordered rates. Methods: This study utilizes county level data, spanning the years 2014 through 2020, on the rates of opioid-related mortality, opioid prescription dispensing, and disability. To successfully estimate and predict trends in these opioid-related factors, we augment the Kalman Filter with a novel spatial component. To define opioid vulnerability profiles, we create heat maps of our filter’s predicted rates across the nation’s counties and identify the hotspots. In this context, hotspots are defined on a year-by-year basis as counties with rates in the top 5% nationally. Results: Our spatial Kalman filter demonstrates strong predictive performance. From 2014 to 2018, these predictions highlight consistent spatiotemporal patterns across all three factors, with Appalachia distinguished as the nation’s most vulnerable region. Starting in 2019 however, the dispensing rate profiles undergo a dramatic and chaotic shift. Conclusions: The initial primary drivers of opioid abuse in the Appalachian region were likely prescription opioids; however, it now appears that abuse is sustained by illegal drugs. Additionally, we find that the disabled subpopulation may be more at risk of opioid-related mortality than the general population. Public health initiatives must extend beyond controlling prescription practices to address the transition to and impact of illicit drug use.

60 APPLIED LIFE SCIENCES↗