Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial heterogeneity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Design Principles in Engineering of Multigrain Nanocatalysts via Multiscale Electronic Structure Characterization

Engineering grain boundary (GB) strain provides a promising pathway to tune the catalytic properties of nanocrystals. However, structural heterogeneity from random grain orientation and geometry has limited clear structure–property correlations. Here, we utilize a multigrain Co3O4/Mn3O4 core/shell nanocrystal platform as a model system to systematically investigate how geometric misfit strain at GBs serves as catalytically active sites for the oxygen reduction reaction. Through precise subnanometer-level control over grain morphology and by integrating multiscale electronic structure characterization, we identify the electronic structural signature of GB defects and establish a direct correlation between localized strain fields and modified electronic states. Strain modulation at GBs alters the eg orbital energy levels, with elongation along the z-axis combined with shear strain stabilizing the eg states, in contrast to the destabilization observed under pure shear strain. This stabilization mechanism enhances the electrocatalytic activity and selectivity of strained GBs compared with strain-relaxed grain surfaces. Furthermore, we reveal that GBs exhibit a radial strain gradient, producing a spatial energy shift that further modulates local electronic structures, as resolved through the classification of electron energy loss spectroscopy data. Together, these findings demonstrate that geometric misfit strain enables precise tuning of grain geometry and the resulting electronic structures, offering a robust strategy for engineering next-generation nanocatalysts.

Cho, Min Gee↗
DOE OSTI2024

The Green's Function Model Intercomparison Project (GFMIP) Protocol

The atmospheric Green's function method is a technique for modeling the response of the atmosphere to changes in the spatial field of surface temperature. While early studies applied this method to changes in atmospheric circulation, it has also become an important tool to understand changes in radiative feedbacks due to evolving patterns of warming, a phenomenon called the “pattern effect.” To better study this method, this paper presents a protocol for creating atmospheric Green's functions to serve as the basis for a model intercomparison project, GFMIP. The protocol has been developed using a series of sensitivity tests performed with the HadAM3 atmosphere-only general circulation model, along with existing and new simulations from other models. Our preliminary results have uncovered nonlinearities in the response of the atmosphere to surface temperature changes, including an asymmetrical response to warming versus cooling patch perturbations, and a change in the dependence of the response on the magnitude and size of the patches. These nonlinearities suggest that the pattern effect may depend on the heterogeneity of warming as well as its location. These experiments have also revealed tradeoffs in experimental design between patch size, perturbation strength, and the length of control and patch simulations. The protocol chosen on the basis of these experiments balances scientific utility with the simulation time and setup required by the Green's function approach. Running these simulations will further our understanding of many aspects of atmospheric response, from the pattern effect and radiative feedbacks to changes in circulation, cloudiness, and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Advances in imaging of chemically reacting flows

Many important chemically reacting systems are inherently multi-dimensional with spatial and temporal variations in the thermochemical state, which can be strongly coupled to interactions with transport processes. Fundamental insights into these systems require multi-dimensional measurements of the thermochemical state as well as fluid dynamics quantities. Laser-based imaging diagnostics provide spatially and temporally resolved measurements that help address this need. The state of the art in imaging diagnostics is continually progressing with the goal of attaining simultaneous multi-parameter measurements that capture transient processes, particularly those that lead to stochastic events, such as localized extinction in turbulent combustion. Development efforts in imaging diagnostics benefit from advances in laser and detector technology. This article provides a perspective on the progression of increasing dimensionality of laser-based imaging diagnostics and highlights the evolution from single-point measurements to 1D and 2D multi-parameter imaging and 3D high-speed imaging. This evolution is demonstrated using highlights of laser-based imaging techniques in combustion science research as an exemplar of a complex multi-dimensional chemically reacting system with chemistry–transport coupling. Imaging diagnostics impact basic research in other chemically reacting systems as well, such as measurements of near-surface gases in heterogeneous catalysis. The expanding dimensionality of imaging diagnostics leads to larger and more complex datasets that require increasingly demanding approaches to data analysis and provide opportunities for increased collaboration between experimental and computational researchers in tackling these challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗
DOE OSTI2022

Designing Unmanned Aerial Survey Monitoring Program to Assess Floating Litter Contamination

Monitoring marine contamination by floating litter can be particularly challenging since debris are continuously moving over a large spatial extent pushed by currents, waves, and winds. Floating litter contamination have mostly relied on opportunistic surveys from vessels, modeling and, more recently, remote sensing with spectral analysis. This study explores how a low-cost commercial unmanned aircraft system equipped with a high-resolution RGB camera can be used as an alternative to conduct floating litter surveys in coastal waters or from vessels. The study compares different processing and analytical strategies and discusses operational constraints. Collected UAS images were analyzed using three different approaches: (i) manual counting (MC), using visual inspection and image annotation with object counts as a baseline; (ii) pixel-based detection, an automated color analysis process to assess overall contamination; and (iii) machine learning (ML), automated object detection and identification using state-of-the-art convolutional neural network (CNNs). Our findings illustrate that MC still remains the most precise method for classifying different floating objects. ML still has a heterogeneous performance in correctly identifying different classes of floating litter; however, it demonstrates promising results in detecting floating items, which can be leveraged to scale up monitoring efforts and be used in automated analysis of large sets of imagery to assess relative floating litter contamination.

Almeida, Sílvia (ORCID:0000000320636826)↗
DOE OSTI2021

Excess deaths reveal the true spatial, temporal and demographic impact of COVID-19 on mortality in Ecuador

Accepted Background In early 2020, Ecuador reported one of the highest surges of per capita deaths across the globe. Methods We collected a comprehensive dataset containing individual death records between 2015 and 2020, from the Ecuadorian National Institute of Statistics and Census and the Ecuadorian Ministry of Government. We computed the number of excess deaths across time, geographical locations and demographic groups using Poisson regression methods. Results Between 1 January and 23 September 2020, the number of excess deaths in Ecuador was 36 402 [95% confidence interval (CI): 35 762–36 827] or 208 per 100 000 people, which is 171% of the expected deaths in that period in a typical year. Only 20% of the excess deaths are attributable to confirmed COVID-19 deaths. Strikingly, in provinces that were most affected by COVID-19 such as Guayas and Santa Elena, the all-cause deaths are more than double the expected number of deaths that would have occurred in a normal year. The extent of excess deaths in men is higher than in women, and the number of excess deaths increases with age. Indigenous populations had the highest level of excess deaths among all ethnic groups. Conclusions Overall, the exceptionally high level of excess deaths in Ecuador highlights the enormous burden and heterogeneous impact of COVID-19 on mortality, especially in older age groups and Indigenous populations in Ecuador, which was not fully revealed by COVID-19 death counts. Together with the limited testing in Ecuador, our results suggest that the majority of the excess deaths were likely to be undocumented COVID-19 deaths.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Electronic and catalytic engineering in two-dimensional vdW metal–organic frameworks through alloying

Bimetallic metal-organic framework (MOF s ) alloys, in which heterogeneous metal clusters are incorporated into their backbone, are capable of highly selective separations and catalysis. Due to limitations in our fundamental understanding of their alloying, however, established methods result in phase-separated or amorphous two-dimensional (2D) MOF s or lack precise control over alloy ratios. Here, our results demonstrate 2D MOF alloys where metal cation ratios (M 1 and M 2 ) in M 1 xM 2 1-x BDC (M 1 or M 2 = Zn, Cu, Ni, Co, Fe, Mn) can be engineered on demand by controlling the metal salt dissociation constants. Resulting MOF alloys exhibit a highly 2D nature with excellent crystallinity and minute control over metal cation ratios. Overall, our experimental and theoretical results show that their electronic bandgaps and photoexcited carrier lifetimes can be engineered by metal cation alloying. Interestingly, 2D alloyed MOF s enable high-efficiency photo-catalytic water reduction performance in Co/Ni MOF alloys owing to the spatially separated metal clusters in 2D MOF alloys.

2D materials↗
DOE OSTI2026

Comparing Delay-, Distance-, and Cordon-Based Congestion Pricing Strategies Via Large-Scale Simulation

This study compares the impacts of delay-, distance-, and cordon-based congestion pricing strategies for Austin, Texas, using the POLARIS agent-based activity-based travel demand simulation model. This approach enables agent-level heterogeneity and realistic choice options (including destination, mode, and activity scheduling) for dynamic traffic assignment and congestion feedbacks across a major metro region, which are features lacking in past work. To ensure comparability, distance-based tolls were set to generate the same revenue as delay-based tolling of $3.5 M/day, averaging $1.17/resident/day or $0.42/vehicle-trip. Delay-based pricing delivers 44% lower network delay and 13% lower VHT compared to the no-toll baseline, levels unmatched by other pricing strategies. At the height of the AM peak, drivers pay up to $0.13/mile on average, though most links in the network remain untolled. Distance-based pricing is the most effective at reducing VMT (by 4%), but VHT reductions (of 6%) primarily stem from drivers selecting closer destinations, achieving only one-fourth the delay reduction of delay-based pricing. Across various implementations of delay- and distance-based pricing, the results suggest that spatial variations of tolls are far more important than temporal variations. Cordon tolls produce minimal impacts at the network-wide level, but offer substantial delay reductions inside the cordon. Other major findings include: 1) delay-based pricing increases trip-making during the PM peak period due to backward shifts in discretionary-activity start times by higher-income residents; and 2) tolls’ spatial impacts, including changes in network flows and tolls paid by residents, vary substantially between delay- and distance-based pricing strategies.

Agent-based modeling↗
DOE OSTI2023

Kinetic Control via Binding Sites within the Confined Space of Metal Metalloporphyrin‐Frameworks for Enhanced Shape‐Selectivity Catalysis

Abstract One striking feature of enzyme is its controllable ability to trap substrates via synergistic or cooperative binding in the enzymatic pocket, which renders the shape‐selectivity of product by the confined spatial environment. The success of shape‐selective catalysis relies on the ability of enzyme to tune the thermodynamics and kinetics for chemical reactions. In emulation of enzyme's ability, we showcase herein a targeting strategy with the substrate being anchored on the internal pore wall of metal‐organic frameworks (MOFs), taking full advantage of the sterically kinetic control to achieve shape‐selectivity for the reactions. For this purpose, a series of binding site‐accessible metal metalloporphyrin‐frameworks (MMPFs) have been investigated to shed light on the nature of enzyme‐mimic catalysis. They exhibit a different density of binding sites that are well arranged into the nanospace with corresponding distances of opposite binding sites. Such a structural specificity results in a facile switch in selectivity from an exclusive formation of the thermodynamically stable product to the kinetic product. Thus, the proposed targeting strategy, based on the combination of porous materials and binding events, paves a new way to develop highly efficient heterogeneous catalysts for shifting selectivity.

Zhang, Weijie↗
DOE OSTI2026

Multiscale Cryo Electron Microscopy Reveals Interfacial Degradation and Stabilization in Battery Electrodes

Electrochemical interfaces are dynamic systems, evolving based on their local environment and reactant surface structures. The electrode-electrolyte interface in Li-ion batteries can be protective, limiting parasitic reactions with the electrolyte to passivate the surface [1]. Additionally, this interphase has an impact on the Li-ion transport through that layer based on its composition, bonding environment, and thickness. These parameters are challenging to collect and may vary depending on the electrode surface site investigated relative to its spatial position in a coin cell. This study will detail a multiscale cryogenic electron microscopy approach where millimeter-scale cross-sections through the coin cell batteries were made using a cryogenic stage within a fs-laser plasma focused ion beam (laser PFIB) with complementary energy dispersive X-ray spectroscopy able to detect variations in the composition at electrode interfaces [2]. Microscale cross-sectioning and lamella sample preparation of battery electrodes was conducted at the Center for Integrated Nanotechnologies using a Ga-ion focused ion beam (FIB) with air-free and cryo-transfer [3], followed by nanoscale mapping of composition and bonding within the CEI through cryo-scanning transmission electron microscopy (cryo-STEM) electron energy loss spectroscopy [4]. This multiscale approach enabled identification of millimeter-scale features of a battery stack with visualization of degradation in electrodes such as cracks in cathode particles, gas evolution, and SEI evolution; microscale interfacial characteristics, such as heterogeneity in the SEI or barrier layer and identification of electrolyte networks to the electrode surfaces; and nanoscale measurement of the CEI thickness, mapping of transition metal bonding within the cathode particles to identify loss of active materials, and identification of beneficial electrolyte additives incorporated into the CEI structure. This multiscale approach allows for a statistical understanding of the primary mechanisms and parasitic degradation pathways that impact performance by limiting the ion transport pathways within Li+ batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of a high fidelity model of the CROCUS experimental reactor

Measurements of scalar flux distributions with fine spatial and energy resolutions are needed to remedy one of the validation shortcomings of the novel neutronics full core solvers, such as MPACT and nTRACER. While a very detailed resolution of the flux can be calculated with such codes, only a limited experimental data set is available to check their accuracy. Such type of measurements are on-going at the zero power reactor CROCUS, operated at the Laboratory for Reactor Physics and System Behaviours of the EPFL, thanks to the development of advanced miniature neutron detection systems. This kind of experimental data would provide the community with a suitable benchmark for the validation of high fidelity neutronics solvers. In parallel, a multi-physics solver for steady-state and transient analysis of nuclear reactors, named GeN-Foam, has been developed. Based on the finite-volume OpenFOAM library, GeN-Foam provides us with enough flexibility to analyze non-conventional reactor geometries such as that of CROCUS. While CROCUS heterogeneities cannot be modeled by MPACT and nTRACER for the moment, GeN-Foam offers a unique opportunity to build a high fidelity model which mimics these codes' method to reach sub-pin simulation resolution. This document aims at describing the work achieved to get from the existing GeN-Foam model of the CROCUS reactor based on a structured mesh and using the neutron diffusion, the first high-fidelity model using discrete ordinates method as an approximation to neutron transport and an unstructured mesh for inter lattice water gap description and sub-pin heterogeneous modeling. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Nanotomography for Quantitative 3D Particle Reconstruction

Particulates are ubiquitous across fuel cycle operations and carry critical information about particle formation, processing, and potential proliferation-related activities. Traditional analytical techniques, including micro-Raman spectroscopy and standard electron microscopy, are often limited in spatial resolution or dimensionality, particularly when used to examine metallic or submicron-scale features. Understanding particle morphology, phase distribution, and internal porosity is essential for constraining formation conditions, thermodynamic environments, and material transport behavior. In this report, we demonstrate the application of plasma focused ion beam nanotomography to reconstruct micron-scale particulates at nanoscale resolution. Using high-resolution backscattered electron imaging and Avizo software, we obtained 3D reconstructions that enabled quantitative analysis of particle morphology, phase composition, and internal voids. Representative examples include a Ta particle with a large central void and a composite particle with embedded tetrahedral crystalline structures. These reconstructions reveal structural and compositional details that are inaccessible through conventional 2D imaging. The results demonstrate that nanotomography provides both qualitative and quantitative insights into particle formation and behavior. Using nanotomography, porosity and phase distributions can be quantified to inform models of particle density, transport, and solidification conditions. Beyond technical insights, the workflow developed here establishes a transferable capability for analyzing heterogeneous particles and has potential applications in bulk materials studies via x-ray computed tomography or other volumetric imaging modalities. Ongoing efforts are focused on optimizing the workflow to process multiple particles simultaneously, increasing throughput and statistical robustness. Overall, this work illustrates the power of nanotomography as a tool for connecting particulate morphology to formation mechanisms, composition, and transport, thereby strengthening analytical capabilities for nuclear forensics, fuel cycle analysis, and related scientific investigations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Characterization of Surface Heterogeneity-Induced Convection Using Cluster Analysis

Land-atmosphere interactions and boundary layer processes often control the formation of shallow clouds and subsequently deep convective precipitation over the Southern Great Plains. In this study, we examine the impacts of large-scale advection on the cloud populations and land-atmospheric coupling observed during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign in 2016. In this work, we performed two Large Eddy Simulations (LES) using the Weather Research and Forecasting Model for a day in which the transition from clear-sky to shallow precipitating clouds and isolated deep convection was observed. The control simulation reproduced the overall distribution of cloud populations by using realistic soil conditions with an interactive land model. In the sensitivity simulation where large-scale advection is removed, a strong relationship between the land and boundary layer is found. To study the timing, location, and intensity of convective initiation and the relationship of clouds with land surface properties, a cluster analysis of equivalent potential temperature is performed for the simulation without large-scale advection. That analysis shows that convective clouds first form over regions with higher surface sensible heat flux. Precipitation from those convective clouds likely triggers new updrafts nearby about two hours later through the lifting associated with cold pools. Furthermore, the cluster analysis also shows that in addition to the spatial pattern of soil moisture, land use and soil texture in western Oklahoma also influence the location of convective initiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Kinetic Control via Binding Sites within the Confined Space of Metal Metalloporphyrin-Frameworks for Enhanced Shape-Selectivity Catalysis

One striking feature of enzyme is its controllable ability to trap substrates via synergistic or cooperative binding in the enzymatic pocket, which renders the shape-selectivity of product by the confined spatial environment. The success of shape-selective catalysis relies on the ability of enzyme to tune the thermodynamics and kinetics for chemical reactions. In emulation of enzyme's ability, we showcase herein a targeting strategy with the substrate being anchored on the internal pore wall of metal-organic frameworks (MOFs), taking full advantage of the sterically kinetic control to achieve shape-selectivity for the reactions. For this purpose, a series of binding site-accessible metal metalloporphyrin-frameworks (MMPFs) have been investigated to shed light on the nature of enzyme-mimic catalysis. They exhibit a different density of binding sites that are well arranged into the nanospace with corresponding distances of opposite binding sites. Furthermore, such a structural specificity results in a facile switch in selectivity from an exclusive formation of the thermodynamically stable product to the kinetic product. Thus, the proposed targeting strategy, based on the combination of porous materials and binding events, paves a new way to develop highly efficient heterogeneous catalysts for shifting selectivity.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Formation and coarsening of epitaxially-supported metal nanoclusters

This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Soil respiration spatial and temporal variability in China between 1961 and 2014

Soil respiration (Rs) plays an important role in terrestrial-atmospheric carbon exchange but remains one of the least studied components of the carbon cycle. How environmental changes influence Rs, and in turn, how Rs influences terrestrial carbon storage in China is unclear. Here, we estimated spatial patterns and temporal trends in Rs from 1961 to 2014 to determine the influence of the recent warming hiatus on the Rs temporal variability in China. We evaluated the relationship between Rs and a set of environmental factors and developed biome-specific Rs models that were parameterized based on 2,862 Rs records and 18 continuous variables (e.g., air temperature, precipitation and leaf area index). Based on the biome-specific Rs models and environmental information, we estimated the annual Rs as well as the change trend of Rs in China during the study period. The mean annual Rs from 1961 to 2014 was 3.50 (confidence interval +/- 0.87) Pg center dot C center dot year(-1). Rs in China increased at a rate of 0.008 Pg center dot C center dot year(-2) between 1982 and 1998 (warming period), but was static between 1999 and 2014 (warming hiatus). This reduction in Rs increase rate may have an important role in enhancing terrestrial carbon sequestration. These findings confirm that temporal heterogeneity in Rs could be influenced by regional environmental changes, and further help in improving the understanding of soil carbon dynamics in China.

Jian, Jinshi↗
DOE OSTI2020

Temporal Resource Continuity Increases Predator Abundance in a Metapopulation Model: Insights for Conservation and Biocontrol

The amount of habitat in a landscape is an important metric for evaluating the effects of land cover on biodiversity, yet it fails to capture complex temporal dimensions of resource availability that could be consequential for species population dynamics. Here, we use a spatially-explicit predator–prey metapopulation model to test the effect of different spatiotemporal resource patterns on insect predators and their prey. We examined population responses in model landscapes that varied in both the amount and temporal variability of basal vegetation. Further, we examined cases where prey comprised either a single generalist species or two specialist species that use different resources available either early or late in the growing season. We found that predators and generalist prey benefitted from lower temporal variance of basal resources, which increased landscape-scale abundances. However, increasing the amount of basal resources also increased the variability of generalist prey populations. Specialist prey, on the other hand, did not benefit from less temporally variable basal resources, as they were restricted by habitat type, while also suffering greater predation. Predators achieved greater prey suppression in landscapes with less temporally variable resources, but the overall effects on prey abundance depended on prey habitat specialization. Our simulations demonstrate the joint importance of both the amount and temporal variability of resources for understanding how landscape heterogeneity influences biodiversity and ecosystem services such as the biological control of agricultural pests.

Spiesman, Brian↗