SEARCH · Engineering Papers

Results for “Population estimates”

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 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

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 OSTI2021

Automating Bayesian inference and design to quantify acoustic particle levitation

Self-propulsion of micro- and nanoparticles powered by ultrasound provides an attractive strategy for the remote manipulation of colloidal matter using biocompatible energy inputs. Quantitative understanding of particle motion and its dependence on size, shape, and composition requires accurate characterization of the acoustic field, which depends sensitively on the experimental setup. Here, we show how automated experiments based on Bayesian inference and design can accurately and efficiently characterize the acoustic field within resonant chambers used to propel acoustic nanomotors. Repeated cycles of observation, inference, and design (OID) are guided by a physical model that describes the rate at which levitating particles approach the nodal plane. Using video microscopy, we observe the relaxation of tracer particles to this plane following the application of the acoustic field. We use sequential Monte Carlo methods to infer model parameters such as the amplitude and frequency of the resonant chamber while accounting for particle-level measurement noise and population-level heterogeneity in the field. Guided by simulated outcomes, we select the optimal design for the next experiment as to maximize the information gain in the relevant parameters. We show how this iterative process serves to discriminate between competing hypotheses and efficiently converges to accurate parameter estimates using only few automated experiments. We discuss the need for model criticism to ensure the validity of the guiding model throughout automated cycles of observation, inference, and design. Furthermore, this work demonstrates how Bayesian methods can learn the parameters of nonlinear, hierarchical models used to describe video microscopy data of active colloids.

42 ENGINEERING↗

DOE OSTI2023

Investigating the vertical distribution of the disk as a function of radial action

As heating processes can broaden the distributions of radial actions and the vertical distributions of the Galactic disks, we investigate the vertical distribution of the Galactic disks as a function of radial action based on Apache Point Observatory Galactic Evolution Experiment (APOGEE) and Gaia data in order to deepen our understanding of the formation and heating history of the Galactic disks. The vertical distributions of the thin and thick Galactic disks defined in the chemical plane were fitted with a simple exponential function with a free parameter of scale height in different radial action ranges. Therefore, we were able to analyze the scale height as a function of radial action for different disk populations. We find that the distributions of radial action for both the thin and thick disks can be approximately described by pseudo-isothermal distributions, which give a statistical measurement for the temperature of the disk as indicated by the mean radial action of the star sample. Estimations of the scale heights in different radial action ranges for these pseudo-isothermal distributions of the disks seem to show fixed relationships between radial action J R and scale height h. We describe these relationships with a two-parameter function of h = √(J R /a)+b, where a and b are free parameters. When testing with a three-parameter function of h = α √(J R /a)+b, we find that this two-parameter function describes the thin disk well, but we note the function should be used with care for the thick disk. When comparing the best-fit relationships between the inner and outer disk for both of the thin and thick disks, we find that the relationships are nearly the same for the thin disks but are different for the thick disks. The inner thick disk shows a nearly flattened relationship, while the outer thick disk presents a gradually increasing relationship. This work highlights an alternative way to unveil the heating history of the disks by investigating the relationship between scale height and radial action, as these relationships encode the formation and heating history of the Galactic disks.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Forecasting fluid-injection induced seismicity to choose the best injection strategy for safety and efficiency

Induced seismicity poses a challenge to the development of Enhanced Geothermal Systems (EGS). Improving monitoring and forecasting techniques is essential to mitigate induced seismicity and thereby fostering a positive perception of EGS projects among local authorities and population. Induced seismicity is the result of complex and coupled thermo-hydro-mechanical-chemical mechanisms. Injection flux and pressure are crucial controlling parameters for both hydraulic stimulation and circulation protocols. We develop a methodology combining a hydro-mechanical model with a seismicity rate model to estimate the magnitude and frequency of mainshocks and aftershocks induced by fluid injection. We apply the methodology to the case of the Basel EGS (2006, Switzerland) to compare the effects of progressive, cyclic and constant injections on the mechanical response of discrete faults. Results from the coupled hydro-mechanical models show that the pore pressure diffusion and consequent enhancement of fault permeability are limited to the vicinity of the injection well during cyclic injection. Additionally, constant injection induces seismicity from the start of the injection but enhances the permeability of most of the faults within a shorter duration, inducing less post-injection seismicity. The methodology can be adapted to any numerical model and allows new projects to be developed by anticipating the safest injection protocol. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.

58 GEOSCIENCES↗

DOE OSTI2025

Nonmagnetic fractional conductance in high mobility InAs quantum point contacts

In this paper, we report the magnetoelectronic properties of high mobility InAs quantum point contacts grown on InP substrates. The InAs quantum well is embedded between In 0.72 ⁢Ga 0.28 ⁢As cladding layers and In 0.81 ⁢Al 0.19 ⁢As barrier layers, and is populated via self-accumulation. The one-dimensional (1D) conductance reaches a maximum value of 17 plateaus, quantized in units of 2⁢𝑒 2 /ℎ, where 𝑒 is the fundamental unit of charge and ℎ is Planck's constant. The in-plane effective 𝑔 factor was estimated to be −10.9 ± 1.5 for subband 𝑁 = 1 and −10.8 ± 1.6 for subband 𝑁 = 2. Furthermore, a study of the nonmagnetic fractional conductance states at 0.2 (𝑒 2 /ℎ) and 0.1 (𝑒 2 /ℎ) is provided. While their origin remains under discussion, evidence suggests that they arise from strong electron-electron interactions and momentum-conserving backscattering between electrons in two distinct channels within the 1D region. This phenomenon may also be interpreted as an entanglement between the two channel directions facilitated by momentum-conserving backscattering.

Ballistic transport↗

DOE OSTI2023

Smoothing Lexis diagrams using kernel functions: A contemporary approach

Lexis diagrams are rectangular arrays of event rates indexed by age and period. Analysis of Lexis diagrams is a cornerstone of cancer surveillance research. Typically, population-based descriptive studies analyze multiple Lexis diagrams defined by sex, tumor characteristics, race/ethnicity, geographic region, etc. Inevitably the amount of information per Lexis diminishes with increasing stratification. Several methods have been proposed to smooth observed Lexis diagrams up front to clarify salient patterns and improve summary estimates of averages, gradients, and trends. In this article, we develop a novel bivariate kernel-based smoother that incorporates two key innovations. First, for any given kernel, we calculate its singular values decomposition, and select an optimal truncation point—the number of leading singular vectors to retain—based on the bias-corrected Akaike information criterion. Second, we model-average over a panel of candidate kernels with diverse shapes and bandwidths. The truncated model averaging approach is fast, automatic, has excellent performance, and provides a variance-covariance matrix that takes model selection into account. We present an in-depth case study (invasive estrogen receptor-negative breast cancer incidence among non-Hispanic white women in the United States) and simulate operating characteristics for 20 representative cancers. The truncated model averaging approach consistently outperforms any fixed kernel. Our results support the routine use of the truncated model averaging approach in descriptive studies of cancer.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

pygwb: A Python-based Library for Gravitational-wave Background Searches

The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the universe and the population of GW sources within it. We present a new, user-friendly, Python-based package for GW data analysis to search for an isotropic GWB in ground-based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one’s own needs. We describe the individual modules that make up pygwb, following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline that combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1979

Plasma property and performance prediction for mercury ion thrusters

The discharge chambers of mercury ion thrusters are modelled so the principal effects and processes which govern discharge plasma properties and thruster performance are described. The conservation relations for mass, charge and energy when applied to the Maxwellian electron population in the ion production region yield equations which may be made one-dimensional by the proper choice of coordinates. Solutions to these equations with the appropriate boundary conditions give electron density and temperature profiles which agree reasonably well with measurements. It is then possible to estimate plasma properties from thruster design data and those operating parameters which are directly controllable. By varying the operating parameter inputs to the computer code written to solve these equations, perfromance curves are obtained which agree quite well with measurements.

Longhurst, G. R.↗

NASA NTRS1992

The collisional lifetime of asteroid 951 Gaspra

We estimate the collisional lifetime of asteroid 951 Gaspra, the first target of an asteroid encounter by an interplanetary space probe. Plausible values range between 10 exp 8 and a few times 10 exp 9 yr, mainly depending on the impact strength of the body and the existing population of sub-km sized projectiles - both quantities for which only indirect observational evidence is currently available. Thus, Gaspra is most likely to be a fragment from a larger asteroid that was shattered in the past. Moreover, unless its material is very resistant to impact break-up (such as iron), Gaspra's surface is among the youngest extraterrestrial surfaces seen in the solar system (apart from active comet nuclei). However, better data on the abundance of very small asteroids are needed to determine its age in a more reliable way.

Farinella, P.↗

NASA NTRS2016

Assessing Woody Vegetation Trends in Sahelian Drylands Using MODIS Based Seasonal Metrics

Woody plants play a major role for the resilience of drylands and in peoples' livelihoods. However, due to their scattered distribution, quantifying and monitoring woody cover over space and time is challenging. We develop a phenology driven model and train/validate MODIS (MCD43A4, 500m) derived metrics with 178 ground observations from Niger, Senegal and Mali to estimate woody cover trends from 2000 to 2014 over the entire Sahel. The annual woody cover estimation at 500 m scale is fairly accurate with an RMSE of 4.3 (woody cover %) and r(exp 2) = 0.74. Over the 15 year period we observed an average increase of 1.7 (+/- 5.0) woody cover (%) with large spatial differences: No clear change can be observed in densely populated areas (0.2 +/- 4.2), whereas a positive change is seen in sparsely populated areas (2.1 +/- 5.2). Woody cover is generally stable in cropland areas (0.9 +/- 4.6), reflecting the protective management of parkland trees by the farmers. Positive changes are observed in savannas (2.5 +/- 5.4) and woodland areas (3.9 +/- 7.3). The major pattern of woody cover change reveals strong increases in the sparsely populated Sahel zones of eastern Senegal, western Mali and central Chad, but a decreasing trend is observed in the densely populated western parts of Senegal, northern Nigeria, Sudan and southwestern Niger. This decrease is often local and limited to woodlands, being an indication of ongoing expansion of cultivated areas and selective logging.We show that an overall positive trend is found in areas of low anthropogenic pressure demonstrating the potential of these ecosystems to provide services such as carbon storage, if not over-utilized. Taken together, our results provide an unprecedented synthesis of woody cover dynamics in theSahel, and point to land use and human population density as important drivers, however only partially and locally offsetting a general post-drought increase.

Brandt, Martin↗

NASA NTRSDate not supplied

A New Statistical Estimate of the Radar Coverage of the Low Earth Orbit Debris Environment

For over three decades, the NASA Orbital Debris Program Office (ODPO) has used the Goldstone Orbital Debris Radar, Haystack Ultrawideband Satellite Imaging Radar (HUSIR), and Haystack Auxiliary (HAX) radar assets to collect data on the low Earth orbit (LEO) debris environment. Each radar, with its unique beamwidth, altitude and inclination coverage, and limiting size threshold, operates in a beam park mode to statistically sample the orbital debris population in LEO. Provided that these assets are shared with other users, the orbital debris data collection is not continuous; rather, intermittent data collects are acquired and sent to the NASA ODPO. To understand the sampling process conducted by each radar over time and any related observational biases, a Statistical Confirmation of Radar Uniformity or Bias (SCRUB) code has been implemented to model the coverage of these assets for informing future operations, as well as usage of the data collected from these ground-based sensors. For this analysis, Right Ascension of the Ascending Node (RAAN) is used as a metric to measure statistical coverage. A complete survey of the LEO environment is understood to be measurements that sufficiently sample all values of RAAN for each altitude-inclination pair visible from the radar asset. Regions of incompleteness or statistical bias can help inform future observation campaigns. The SCRUB tool evaluates the coverage of the LEO environment, not by examining individual objects that may pass through a sensor’s field of view (FOV), but by determining which orbit planes pass through the FOV. Once a pointing geometry for a radar site, observation time, and range extent are configured by the user, SCRUB computes the inclinations and altitudes that are visible by the sensor. For each inclination-altitude pair, there is a distinct pair of possible RAAN values, corresponding to the ascending and descending orbit passing through that point in the sensor cone. Repeating this process for all points in the beam, and for multiple time periods during an observation window, creates a matrix of all inclination-altitude-RAAN combinations that are visible during a sensor run. This process can then be repeated for all observations within a year (for an annual survey), propagating all the RAAN values to a common epoch, typically the start of a calendar year, and combining the observations to assemble a full estimate of the RAAN coverage of the LEO environment. This RAAN coverage can then be analyzed for uniformity of sampling, within a certain inclination-altitude pair, or between larger regions of the space environment. This paper provides a general overview of the radar assets utilized by ODPO, typical analysis data products assuming circular orbits, and a discussion of the algorithms that feed between modeling and measurement operations. Following the description of the algorithm, estimates of coverage using HUSIR radar data collected from multiple years are developed and compared.

Chris L Ostrom↗

DOE OSTI2020

Predicting the proximity to macroscopic failure using local strain populations from dynamic in situ X-ray tomography triaxial compression experiments on rocks

Predicting the proximity of large-scale dynamic failure is a critical concern in the engineering and geophysical sciences. Here we use evolving contractive, dilatational, and shear strain deformation preceding failure in dynamic X-ray tomography experiments to examine which strain components best predict the proximity to failure. We develop machine learning models to predict the proximity to failure using time series of three-dimensional local incremental strain tensor fields acquired in rock deformation experiments under stress conditions of the upper crust. Three-dimensional scans acquired in situ throughout triaxial compression experiments provide a distribution of density contrasts from which we estimate the three-dimensional incremental strain that accumulates between each scan acquisition. Training machine learning models on multiple experiments of six rock types provides suites of feature importance that indicate the predictive power of each feature. Comparing the average importance of groups of features that include information about each strain component quantifies the ability of the contractive, dilatational and shear strain to predict the proximity of macroscopic failure. A total of 24 models of four machine learning algorithms with six rock types indicate that 1) the dilatational strain provides the best predictive power of the strain components, and 2) the intermediate values (25th-75th percentile) of the strain population provide the best predictive power of the statistics of the strain populations. In addition, the success of the predictions of models trained on one rock type and tested on other rock types quantifies the similarities and differences of the precursory strain accumulation process in the six rock types. These similarities suggest the potential existence of a unified theory of brittle rock deformation for a range of rock types.

58 GEOSCIENCES↗

DOE OSTI2021

Hurricanes alter 10 Be concentrations in tropical river sediment but do not change regional erosion rate estimates

Tropical islands, including many in island arcs, are subjected to recurring disturbances from extreme storms such as tropical cyclones. To test whether such storms influence cosmogenic nuclide concentrations such that they do not reflect long-term rates of erosion, we measured meteoric and in situ 10Be in river sediment samples from Dominica, an andesitic island in the Caribbean, before and after category five Hurricane Maria (2017). Populations of before- and after-storm concentrations are statistically indistinguishable (n=7 pairs for in situ 10 Be, n=11 pairs for meteoric 10 Be). 10 Be concentrations vary from -138% to +73% within before-after sample pairs relative to the mean of the pair. These new data suggest that the effects of extreme storms on the depth and amount of near-surface erosion on Dominica vary spatially. Our data support calculations of Niemi, et al. (2005) and Yanites, et al. (2009) suggesting that basin-by-basin comparisons of erosion rates based on cosmogenic nuclides should be approached with caution in small (<~100 km 2 ) watersheds affected by mass movements and extreme storms. Erosion rates determined from in situ 10Be in Dominica (geometric mean=0.102 mm y-1, n=12) are low compared to similarly steep and wet areas globally and correlate positively with the spatial density of mass movements.

58 GEOSCIENCES↗

DOE OSTI2024

Coordinated JWST Imaging of Three Distance Indicators in a Supernova Host Galaxy and an Estimate of the Tip of the Red Giant Branch Color Dependence

Boasting a 6.5 m mirror in space, JWST can increase by several times the number of supernovae (SNe) to which a redshift-independent distance has been measured with a precision distance indicator (e.g., tip of the red giant branch (TRGB) or Cepheids); the limited number of such SN calibrators currently dominates the uncertainty budget in distance ladder Hubble constant (H 0 ) experiments. JWST/NIRCAM imaging of the Virgo Cluster galaxy NGC 4536 is used here to preview JWST program GO-1995, which aims to measure H 0 using three stellar distance indicators (Cepheids, TRGB, and J-branch asymptotic giant branch/carbon stars). Each population of distance indicator was here successfully detected—with sufficiently large number statistics, well-measured fluxes, and characteristic distributions consistent with ingoing expectations—so as to confirm that we can acquire distances from each method precise to about 0.05 mag (statistical uncertainty only). We leverage overlapping Hubble Space Telescope imaging to identify TRGB stars, crossmatch them with the JWST photometry, and present a preliminary constraint on the slope of the TRGB's F115W versus (F115W – F444W) relation equal to -0.99 ± 0.16 mag mag -1 . This slope is consistent with prior slope measurements in the similar Two Micron All-Sky Survey J band, as well as with predictions from the BaSTI isochrone suite. We use the new TRGB slope estimate to flatten the 2D TRGB feature and measure a (blinded) TRGB distance relative to a set of fiducial TRGB colors, intended to represent the absolute fiducial calibrations expected from geometric anchors such as NGC 4258 and the Magellanic Clouds. In doing so, we empirically demonstrate that the TRGB can be used as a standardizable candle at the IR wavelengths accessible with JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Estimating the Acute Health Impacts of Fire-Originated PM 2.5 Exposure During the 2017 California Wildfires: Sensitivity to Choices of Inputs

Exposure to wildfire smoke increases the risk of respiratory and cardiovascular hospital admissions. Health impact assessments, used to inform decision-making processes, characterize the health impacts of environmental exposures by combining preexisting epidemiological concentration–response functions (CRFs) with estimates of exposure. These two key inputs influence the magnitude and uncertainty of the health impacts estimated, but for wildfire-related impact assessments the extent of their impact is largely unknown. We first estimated the number of respiratory, cardiovascular, and asthma hospital admissions attributable to fire-originated PM 2.5 exposure in central California during the October 2017 wildfires, using Monte Carlo simulations to quantify uncertainty with respect to the exposure and epidemiological inputs. We next conducted sensitivity analyses, comparing four estimates of fire-originated PM 2.5 and two CRFs, wildfire and nonwildfire specific, to understand their impact on the estimation of excess admissions and sources of uncertainty. We estimate the fires accounted for an excess 240 (95% CI: 114, 404) respiratory, 68 (95% CI: -10, 159) cardiovascular, and 45 (95% CI: 18, 81) asthma hospital admissions, with 56% of admissions occurring in the Bay Area. Although differences between impact assessment methods are not statistically significant, the admissions estimates' magnitude is particularly sensitive to the CRF specified while the uncertainty is most sensitive to estimates of fire-originated PM 2.5 . Not accounting for the exposure surface's uncertainty leads to an underestimation of the uncertainty of the health impacts estimated. Employing context-specific CRFs and using accurate exposure estimates that combine multiple data sets generates more certain estimates of the acute health impacts of wildfires.

58 GEOSCIENCES↗

DOE OSTI2025

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗

DOE OSTI2023

A complete catalogue of broad-line AGNs and double-peaked emission lines from MaNGA integral-field spectroscopy of 10K galaxies: stellar population of AGNs, supermassive black holes, and dual AGNs

We analyse the integral field spectroscopy data for the ≈10 000 galaxies in final data release of the MaNGA survey. We identify 188 galaxies for which the emission lines cannot be described by single Gaussian components. These galaxies can be classified into (1) 38 galaxies with broad $\rm H\alpha$ and [O III] $\rm \lambda$5007 lines, (2) 101 galaxies with broad $\rm H\alpha$ lines but no broad [O III] $\rm \lambda$5007 lines, and (3) 49 galaxies with double-peaked narrow emission lines. Most of the broad-line galaxies are classified as active galactic nuclei (AGNs) from their line ratios. The catalogue helps us further understand the AGN-galaxy coevolution through the stellar population of broad-line region host galaxies and the relation between broad lines’ properties and the host galaxies’ dynamical properties. The stellar population properties (including mass, age, and metallicity) of broad-line host galaxies suggest there is no significant difference between narrow-line Seyfert-2 galaxies and Type-1 AGNs with broad $\rm H\alpha$ lines. We use the broad-$\rm H\alpha$ line width and luminosity to estimate masses of black hole in these galaxies, and test the MBH–σe relation in Type-1 AGN host galaxies. Furthermore, we find three dual AGN candidates supported by radio images from the VLA FIRST survey. This sample may be useful for further studies on AGN activities and feedback processes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Joint Analysis of the Epidemic Evolution and Human Mobility During the First Wave of COVID-19 in Spain: Retrospective Study

Background The initial wave of the COVID-19 pandemic placed a tremendous strain on health care systems worldwide. To mitigate the spread of the virus, many countries implemented stringent nonpharmaceutical interventions (NPIs), which significantly altered human behavior both before and after their enactment. Despite these efforts, a precise assessment of the impact and efficacy of these NPIs, as well as the extent of human behavioral changes, remained elusive. Objective In this study, we conducted a retrospective analysis of the initial wave of COVID-19 in Spain to better comprehend the influence of NPIs and their interaction with human behavior. Such investigations are vital for devising future mitigation strategies to combat COVID-19 and enhance epidemic preparedness more broadly. Methods We used a combination of national and regional retrospective analyses of pandemic incidence alongside large-scale mobility data to assess the impact and timing of government-implemented NPIs in combating COVID-19. Additionally, we compared these findings with a model-based inference of hospitalizations and fatalities. This model-based approach enabled us to construct counterfactual scenarios that gauged the consequences of delayed initiation of epidemic response measures. Results Our analysis demonstrated that the pre–national lockdown epidemic response, encompassing regional measures and heightened individual awareness, significantly contributed to reducing the disease burden in Spain. The mobility data indicated that people adjusted their behavior in response to the regional epidemiological situation before the nationwide lockdown was implemented. Counterfactual scenarios suggested that without this early epidemic response, there would have been an estimated 45,400 (95% CI 37,400-58,000) fatalities and 182,600 (95% CI 150,400-233,800) hospitalizations compared to the reported figures of 27,800 fatalities and 107,600 hospitalizations, respectively. Conclusions Our findings underscore the significance of self-implemented prevention measures by the population and regional NPIs before the national lockdown in Spain. The study also emphasizes the necessity for prompt and precise data quantification prior to enacting enforced measures. This highlights the critical interplay between NPIs, epidemic progression, and human behavior. This interdependence presents a challenge in predicting the impact of NPIs before they are implemented.

Public, Environmental & Occupational Health↗