Comparison of Methods for Mapping Rhizosphere Processes in the Context of Their Surrounding Root and Soil Environments
Not Available
59 BASIC BIOLOGICAL SCIENCES↗
SEARCH · Search NASA
Results for “spatiotemporal analysis”
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.
Assessment of the Biomass Productivity Decline in the Lower Mekong Basin
This study aimed to delineate the geographic hotspots of negative trends in biomass productivity in the Lower Mekong Basin countries (Vietnam, Cambodia, Laos, and Thailand) and identify correlated regional environmental and anthropogenic factors. A long-term time-series (1982–2015) of Normalized Difference Vegetation Index at a resolution of approximately 9.16 km x 9.16 km was used to specify the areas with significant decline or increase in productivity. The relationships between vegetation changes and land attributes, such as climate, population density, soil/terrain conditions, and land-cover types, were examined. Rainfall time-series maps were used to identify areas that might have been affected by land degradation from those correlated with rainfall. Most of the detected potentially degraded areas were found in Cambodia, the Northwest and the Highland of Vietnam, the Northern Mountains of Thailand and Laos, and the mountainous border between Laos, Vietnam, and Cambodia. About 15% of the total land area of these four countries experienced a reduction in biomass productivity during the 34-year study period. The map of hotspots of changes in productivity can be used to direct further studies, including those at finer spatial resolution that may support policy makers and researchers in targeting the strategies for combating land degradation.
land degradation↗
Benchmarking the PCMCI Causal Discovery Algorithm for Spatiotemporal Systems
Causal discovery algorithms construct hypothesized causal graphs that depict causal dependencies among variables in observational data. While powerful, the accuracy of these algorithms is highly sensitive to the underlying dynamics of the system in ways that have not been fully characterized in the literature. In this report, we benchmark the PCMCI causal discovery algorithm in its application to gridded spatiotemporal systems. Effectively computing grid-level causal graphs on large grids will enable analysis of the causal impacts of transient and mobile spatial phenomena in large systems, such as the Earth’s climate. We evaluate the performance of PCMCI with a set of structural causal models, using simulated spatial vector autoregressive processes in one- and two-dimensions. We develop computational and analytical tools for characterizing these processes and their associated causal graphs. Our findings suggest that direct application of PCMCI is not suitable for the analysis of dynamical spatiotemporal gridded systems, such as climatological data, without significant preprocessing and downscaling of the data. PCMCI requires unrealistic sample sizes to achieve acceptable performance on even modestly sized problems and suffers from a notable curse of dimensionality. This work suggests that, even under generous structural assumptions, significant additional algorithmic improvements are needed before causal discovery algorithms can be reliably applied to grid-level outputs of earth system models.
54 ENVIRONMENTAL SCIENCES↗
Characterizing spatiotemporal patterns of crop phenology across North America during 2000–2016 using satellite imagery and agricultural survey data
Crop phenology represents an integrative indicator of climate change and plays a vital role in terrestrial carbon dynamics and sustainable agricultural development. However, spatiotemporal variations of crop phenology remain unclear at large scales. This knowledge gap has hindered our ability to realistically quantify the biogeochemical dynamics in agroecosystems, predict future climate, and make informed decisions for climate change mitigation and adaptation. In this study, we improved an EVI-curve-based approach and used it to detect spatiotemporal patterns in cropping intensity and five major phenological stages over North America during 2000–2016 using vegetation index in combination with agricultural survey data and other ancillary maps. Our predicted crop phenological stages showed strong linear relationships with the survey-based datasets, with R(exp 2), RMSEs, and MAEs in the ranges of 0.35 –0.99, three to ten days, and two to eight days, respectively. During the study period, the planting dates were advanced by 0.60 days/year (p < 0.01), and harvesting dates were delayed by 0.78 days/year (p < 0.01) over North America. A minimum temperature increase by 1 °C caused a 4.26-day planting advance (r = −0.50, p < 0. 01) or a 0.66-day harvest delay (r = 0.10, p < 0.01). While, a higher maximum temperature resulted in a planting advance by 4.48 days/°C (r = −0.62, p < 0.01) or a harvest advance by 2.22 days/°C (r = −0.40, p < 0.01). Our analysis illustrated evident spatiotemporal variations in crop phenology in response to climate change and management practices. The derived crop phenological datasets and cropping intensity maps can be used in regional climate assessments and in developing adaptation strategies.
North America↗
An evaluation of air quality in major urban areas of India
Rapid economic growth and burgeoning population have contributed to enhanced levels of PM 2.5 concentrations in urban regions of India. Evaluation of ambient air quality facilitates the assessment of effectiveness of emission control measures and early identification of new sources. This study provides a comprehensive statistical analysis of PM 2.5 concentrations in key urban areas across India, including Delhi, Kolkata, Mumbai, Chennai, Hyderabad, and several regional centers. Data from 2017 to 2023 was analyzed using trend analysis, cluster analysis, principal component analysis, and geostatistical interpolation to understand spatiotemporal variations and sources. The analysis reveals significant differences in spatial distribution of PM 2.5 concentrations with high annual averages in urban regions in Indo-Gangetic plain (82–123 μg m −3 ) and relatively lower concentrations (29–46 μg m −3 ) in southern urban areas of Kerala, Tamil Nadu and Andhra Pradesh. Delhi state had the highest 24-averaged PM 2.5 concentrations (112 μg m −3 ) followed by urban regions in Uttar Pradesh, Bihar and West Bengal (94 μg m −3 ). Trend analysis from 2017 to 2023 revealed an overall 2.5% decline in site-wide PM2.5 concentrations, with the exception of Ludhiana, which exhibited a consistent annual increase of 10%. Principal component analysis (PCA) attributes 30% of the variance to wintertime emissions, 13% to biomass burning, and 18% to the regional haze in the northern Indo-Gangetic Plain. Different analyses clearly demonstrates the contribution of biomass burning to pollution in Delhi and surrounding cities. Transboundary pollution to Kolkata is likely from the highly polluted region in Indo-Gangetic Plain. Coastal cities of Mumbai and Chennai has relatively lower pollution attributed to the influence of sea breeze dilution, with mostly local contribution and some potential transport from upwind industry clusters. Hyderabad also has local contribution due to high density of vehicular traffic and local small industries. This study shows that mitigation efforts targeting clusters of regions should be undertaken to curb the high PM2.5 pollution. Policy measures should be implemented both at local and the intra-state level to address shared sources and transport of pollution.
Hysplitbacktrajectories↗
Spatiotemporal Associations Between Social Vulnerability, Environmental Measurements, and COVID-19 in the Conterminous United States
This study summarizes the results from fitting a Bayesian hierarchical spatiotemporal model to coronavirus disease 2019 (COVID-19) cases and deaths at the county level in the United States for the year 2020. Two models were created, one for cases and one for deaths, utilizing a scaled Besag, York, Mollié model with Type I spatial-temporal interaction. Each model accounts for 16 social vulnerability and 7 environmental variables as fixed effects. The spatial pattern between COVID-19 cases and deaths is significantly different in many ways. The spatiotemporal trend of the pandemic in the United States illustrates a shift out of many of the major metropolitan areas into the United States Southeast and Southwest during the summer months and into the upper Midwest beginning in autumn. Analysis of the major social vulnerability predictors of COVID-19 infection and death found that counties with higher percentages of those not having a high school diploma, having non-White status and being Age 65 and over to be significant. Among the environmental variables, above ground level temperature had the strongest effect on relative risk to both cases and deaths. Hot and cold spots, areas of statistically significant high and low COVID-19 cases and deaths respectively, derived from the convolutional spatial effect show that areas with a high probability of above average relative risk have significantly higher Social Vulnerability Index composite scores. The same analysis utilizing the spatiotemporal interaction term exemplifies a more complex relationship between social vulnerability, environmental measurements, COVID-19 cases, and COVID-19 deaths.
spatial epidemiology↗
Applying Gaussian Process Machine Learning and Modern Probabilistic Programming to Satellite Data to Infer CO 2 Emissions
Satellite data provides essential insights into the spatiotemporal distribution of CO 2 concentrations. However, many atmospheric inverse models fail to adequately incorporate the spatial and temporal correlations inherent in satellite observations and often lack rigorous methods for estimating parameters like spatial length scales. We introduce an inference model that processes the spatiotemporal covariance in satellite data and estimates hyperparameters such as covariance length scales. Our approach uses the Gaussian process (GP) machine learning (ML) and modern probabilistic programming languages (PPLs) to perform atmospheric inversions of emissions from satellite data. We develop a GP ML inversion system based on modern PPLs and the GEOS-Chem chemical transport model, simulating atmospheric CO 2 concentrations corresponding to the Orbiting Carbon Observatory-2/3 (OCO-2/3) data for July 2020. In our supervised learning framework, we treat the GEOS-Chem simulated data set as the target, with predictors derived by scaling the target with sector-specific factors hidden from the GP machine. Our results show that the GP model, combined with GPU-enabled PPLs, effectively retrieves true emission scaling factors and infers noise levels concealed within the data. This suggests that our method could be applied over larger areas with more complex covariance structures, enabling comprehensive analysis of the spatiotemporal patterns observed in OCO-2/3 and similar satellite data sets.
54 ENVIRONMENTAL SCIENCES↗
Integrating Maximum Entropy Production Theory and Machine Learning to Improve Global Evapotranspiration Modeling
Accurate estimation of terrestrial evapotranspiration (ET) is vital for understanding global water and energy cycles. However, current global ET estimations are not well constrained. This study introduces an integrated framework combining the Maximum Entropy Production (MEP) theory with Random Forest (RF) model to improve global ET estimation. Specifically, in contrast to direct ET estimation by the RF model, the integrated framework (MEP‐RF) trains to predict error of MEP‐simulated ET. MEP‐RF outperforms RF in spatiotemporal extrapolation. Attribution analysis with in situ observations reveals that the inputs of MEP are the most critical variables for the ET process, including net radiation, vegetated area, soil moisture, and surface temperature. We further drive MEP‐RF with global reanalysis and satellite data sets of these four inputs, yielding a global mean terrestrial ET of 548 mm/year, with 77% attributed to transpiration. The global ET increased at a rate of 0.85 mm/year per year during 2003–2021, primarily due to vegetation greening rather than rising temperature, while decreasing soil moisture led to decreasing regional ET. The integrated framework provides a novel approach for the estimation of global ET without the need for hard‐to‐obtain and thus uncertain inputs, such as wind speed, surface roughness, aerodynamic and canopy stomatal resistance. Therefore, MEP‐RF offers an independent method on existing global ET products. It represents a promising physically based approach that can be incorporated into Earth System Models to enhance water and energy cycle simulations.
54 ENVIRONMENTAL SCIENCES↗
Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition
Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.
60 APPLIED LIFE SCIENCES↗
Absolute Radiation Measurements in Earth and Mars Entry Conditions
This paper reports on the measurement of radiative heating for shock heated flows which simulate conditions for Mars and Earth entries. Radiation measurements are made in NASA Ames' Electric Arc Shock Tube at velocities from 3-15 km/s in mixtures of N2/O2 and CO2/N2/Ar. The technique and limitations of the measurement are summarized in some detail. The absolute measurements will be discussed in regards to spectral features, radiative magnitude and spatiotemporal trends. Via analysis of spectra it is possible to extract properties such as electron density, and rotational, vibrational and electronic temperatures. Relaxation behind the shock is analyzed to determine how these properties relax to equilibrium and are used to validate and refine kinetic models. It is found that, for some conditions, some of these values diverge from non-equilibrium indicating a lack of similarity between the shock tube and free flight conditions. Possible reasons for this are discussed.
Radiative Heating↗
Stratospheric Instrusion Catalog: A 10-Year Compilation of Events Identified by using TRACK with NASA's MERRA-2 Reanalysis
Stratospheric intrusions "the introduction of ozone-rich stratospheric air into the troposphere" have been linked with surface ozone air quality exceedances, especially at the high elevations in the western USA in springtime. However, the impact of stratospheric intrusions in the remaining seasons and over the rest of the USA is less clear. A new approach to the study of stratospheric intrusions uses NASA's Goddard Earth Observing System Model (GEOS) model and assimilation products with an objective feature tracking algorithm to investigate the atmospheric dynamics that generate stratospheric intrusions and the different mechanisms through which stratospheric intrusions may influence tropospheric chemistry and surface air quality seasonally over both the western and the eastern USA. A catalog of stratospheric intrusions identified in the MERRA-2 reanalysis was produced for the period 2004-2015 and validated against surface ozone observations (focusing on those which exceed the national air quality standard) and a recent data set of stratospheric intrusion-influenced air quality exceedance flags from the US Environmental Protection Agency (EPA). Considering not all ozone exceedances have been flagged by the EPA, a collection of stratospheric intrusions can support air quality agencies for more rapid identification of the impact of stratospheric air on surface ozone and demonstrates that future operational analyses may aid in forecasting such events. An analysis of the spatiotemporal variability of stratospheric intrusions over the continental US was performed, and while the spring over the western USA does exhibit the largest number of stratospheric intrusions affecting the lower troposphere, the number of intrusions in the remaining seasons and over the eastern USA is sizable. By focusing on the major modes of variability that influence weather in the USA, such as the Pacific North American (PNA) teleconnection index, predicative meteorological patterns associated with stratospheric intrusions and their regional effects on tropospheric ozone were identified. Improved understanding of the connections between large-scale climate variability and local-scale dynamically-driven air quality events may support improved seasonal prediction of such events.
PN↗
Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter
Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.
Williford, Kenneth H.↗
Identification and correction of temporal and spatial distortions in scanning transmission electron microscopy
Scanning transmission electron microscopy (STEM) has become the technique of choice for quantitative characterization of atomic structure of materials, where the minute displacements of atomic columns from high-symmetry positions can be used to map strain, polarization, octahedra tilts, and other physical and chemical order parameter fields. The latter can be used as inputs into mesoscopic and atomistic models, providing insight into the correlative relationships and generative physics of materials on the atomic level. However, these quantitative applications of STEM necessitate understanding the microscope induced image distortions and developing the pathways to compensate them both as part of a rapid calibration procedure for in situ imaging, and the post-experimental data analysis stage. Here, we explore the spatiotemporal structure of the microscopic distortions in STEM using multivariate analysis of the atomic trajectories in the image stacks. Based on the behavior of principal component analysis (PCA), we develop the Gaussian process (GP)-based regression method for quantification of the distortion function. The limitations of such an approach and possible strategies for implementation as a part of in-line data acquisition in STEM are discussed. Here, the analysis workflow is summarized in a Jupyter notebook that can be used to retrace the analysis and analyze the reader's data.
36 MATERIALS SCIENCE↗
Spectrometer-free quantitative vapor sensing and classification via spatiotemporal imaging of porous silicon metasurfaces
Metasurfaces offer a compact platform for optical vapor sensing, but their practical deployment has been limited by weak evanescent light–matter interactions and reliance on spectrally resolved instrumentation. Here, we report porous silicon (pSi) metasurfaces for spectrometer-free quantitative detection of volatile organic compounds (VOCs) with strongly enhanced light–matter interaction. The engineered porosity increases sensitivity by >100× relative to non-porous dielectric metasurfaces, enabling limits of detection of 1.65 ppm for methanol and 9.1 ppm for ethanol across a broad dynamic range (<10 ppm to >103 ppm). Imaging-based readout provides a lightweight, spectrometer-free pathway for real-time quantitative sensing. Beyond quantitative detection, the mesoporous architecture introduces adsorption–desorption kinetics as an additional information channel. Analysis of the resulting spatiotemporal signatures enables kinetic fingerprinting without reliance on infrared spectral features or surface functionalization, and a lightweight machine-learning classifier differentiates acetone, methanol, ethanol, and isopropanol with 91.6% accuracy. These results establish porous metasurfaces as spatiotemporal sensing elements that couple quantitative vapor detection with kinetic fingerprinting through real-time dynamical responses, enabling low-cost, high-performance optical sensors.
Dash, Tomoshree [Clemson University]↗
Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep
Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.
58 GEOSCIENCES↗
Bias Correction and Statistical Downscaling of Future Solar Irradiance Projections Using the NSRDB
Assessing renewable energy resources under future climate scenarios has been highlighted to understand potential impacts of future climate change in renewable generation on the power sector. Climate model projection has been recognized by the renewable energy community as a useful data set to analyze the impacts of future climate change on renewable resources. However, future climate projections generated from general circulation models (GCMs) contain inherent biases that need to be corrected for accurate analysis of future projections of climate variables. In addition, the coarse spatiotemporal resolution of GCMs needs to be improved for regional climate studies. In this work, we develop statistical methods to downscale future projections of global horizontal irradiance (GHI) in a computationally efficient way. Our approach builds statistical downscaling models that correct bias of climate projection of GHI and downscale the future GHI projection from daily-scale to hourly-scale. The National Solar Radiation Database (NSRDB) is used to calibrate the statistical models and validate the downscaled GHI projections across the contiguous United State (CONUS). Preliminary results show that the statistical approach efficiently downscales climate projections of GHI with a nBIAS of 3%, nMAE of 34 % and nRMSE of 46% calculated against NSRDB for CONUS. This study describes the implemented methodology and initial results as well as future research to create high-resolution climate data sets for solar energy applications.
analytical models↗
Investigating the opioid epidemic across the United States: Associations between county-level characteristics and overdose mortality
The opioid crisis remains a critical public health challenge in the United States. Despite national efforts that reduced opioid prescribing by nearly 44% between 2011 and 2021, opioid overdose deaths more than tripled during the same period. This alarming trend reflects a major shift in the crisis, with illegal opioids now driving the majority of overdose deaths instead of prescription opioids. Although supply-side factors fueling this transition have been widely studied, the structural and community-level conditions that shape overdose mortality are less well understood. To help address this gap, this study has three primary objectives: (1) overcome structural gaps in national data to construct a complete nationwide county-level dataset from 2010 to 2022; (2) using data analysis, identify and investigate spatiotemporal anomalies in overdose mortality; and (3) using two machine-learning models, quantify the importance of thirteen social vulnerability variables in predicting overdose mortality. Our results identify unemployment and limited vehicle access as key county-level predictors of overdose mortality. Higher levels of these vulnerabilities are associated with elevated mortality, whereas lower levels are associated with reduced mortality. These findings highlight factors that may be relevant for public health planning and policy prioritization within the context of the opioid crisis.
Anomaly analysis↗
A New Study of Sea Spray Optical Properties from Multi-Sensor Spaceborne Observations
Retrievals of aerosol optical depth (AOD) from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) satellite sensor require the assumption of an extinction-to-backscatter ratio, also known as the lidar ratio. This paper evaluates a new method to calculate lidar ratio of sea spray aerosol using two independent sources: the AOD from Synergized Optical Depth of Aerosols (SODA) and the integrated attenuated backscatter from CALIOP. The method applied in this study allows particulate lidar ratio to be calculated for individual CALIOP retrievals of single aerosol layer columns over the ocean. Analyses are carried out using CALIOP level 2, 5km sea spray aerosol layer products and collocated SODA nighttime data from December 2007 to December 2009. The global mean lidar ratio for sea spray aerosols was found to be 26 sr, roughly 30 higher than the one prescribed by CALIOP. Data analysis also showed considerable spatiotemporal variability in calculated lidar ratio over different parts of the remote oceans. The calculated aerosol lidar ratios are shown to be inversely related to the mean ocean surface wind speed: increase in ocean surface wind speed (U10) from 0 to 15 ms-1 reduces the mean lidar ratios for sea spray particles from 32 sr (for 0 U10 4 ms-1) to 22 sr (for U10 15 ms-1). Such changes in the lidar ratio are expected to have a corresponding effect on sea spray AOD. The outcomes of this study are relevant for future improvements of the SODA and CALIOP operational product and could lead to more accurate retrievals of sea spray AOD.
aerosol optical depth↗