skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Relating the Strength of Density Dependence and the Spatial Distribution of Individuals
Spatial patterns in ecology contain useful information about underlying mechanisms and processes. Although there are many summary statistics used to quantify these spatial patterns, there are far fewer models that directly link explicit ecological mechanisms to observed patterns easily derived from available data. We present a model of intraspecific spatial aggregation that quantitatively relates static spatial patterning to negative density dependence. Individuals are placed according to the colonization rule consistent with the Maximum Entropy Theory of Ecology (METE), and die with probability proportional to their abundance raised to a power α, a parameter indicating the degree of density dependence. This model can therefore be interpreted as a hybridization of MaxEnt and mechanism. Our model shows quantitatively and generally that increasing density dependence randomizes spatial patterning. α = 1 recovers the strongly aggregated METE distribution that is consistent with many ecosystems empirically, and as α → 2 our prediction approaches the binomial distribution consistent with random placement. For 1 < α < 2, our model predicts more aggregation than random placement but less than METE. We additionally relate our mechanistic parameter α to the statistical aggregation parameter k in the negative binomial distribution, giving it an ecological interpretation in the context of density dependence. We use our model to analyze two contrasting datasets, a 50 ha tropical forest and a 64 m 2 serpentine grassland plot. For each dataset, we infer α for individual species as well as a community α parameter. We find that α is generally larger in the tightly packed forest than the sparse grassland, and the degree of density dependence increases at smaller scales. These results are consistent with current understanding in both ecosystems, and we infer this underlying density dependence using only empirical spatial patterns. Our model can easily be applied to other datasets where spatially explicit data are available.  more » « less
Award ID(s):
1751380
PAR ID:
10301156
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Frontiers in Ecology and Evolution
Volume:
9
ISSN:
2296-701X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Maximum Entropy Theory of Ecology (METE) predicts the shapes of macroecological metrics in relatively static ecosystems, across spatial scales, taxonomic categories and habitats, using constraints imposed by static state variables. In disturbed ecosystems, however, with time‐varying state variables, its predictions often fail. We extend macroecological theory from static to dynamic by combining the MaxEnt inference procedure with explicit mechanisms governing disturbance. In the static limit, the resulting theory, DynaMETE, reduces to METE but also predicts a new scaling relationship among static state variables. Under disturbances, expressed as shifts in demographic, ontogenic growth or migration rates, DynaMETE predicts the time trajectories of the state variables as well as the time‐varying shapes of macroecological metrics such as the species abundance distribution and the distribution of metabolic rates over individuals. An iterative procedure for solving the dynamic theory is presented. Characteristic signatures of the deviation from static predictions of macroecological patterns are shown to result from different kinds of disturbance. By combining MaxEnt inference with explicit dynamical mechanisms of disturbance, DynaMETE is a candidate theory of macroecology for ecosystems responding to anthropogenic or natural disturbances. 
    more » « less
  2. null (Ed.)
    ABSTRACT We present a new catalogue of Damped Lyman-α absorbers from SDSS DR16Q, as well as new estimates of their statistical properties. Our estimates are computed with the Gaussian process models presented in Garnett et al., Ho, Bird & Garnett with an improved model for marginalizing uncertainty in the mean optical depth of each quasar. We compute the column density distribution function (CDDF) at 2 < z < 5, the line density (dN/dX), and the neutral hydrogen density (ΩDLA). Our Gaussian process model provides a posterior probability distribution of the number of DLAs per spectrum, thus allowing unbiased probabilistic predictions of the statistics of DLA populations even with the noisiest data. We measure a non-zero column density distribution function for $$N_{\rm {HI}}\lt 3 \times 10^{22} \, \rm {cm}^{-2}$$ with $$95{{\ \rm per\ cent}}$$ confidence limits, and $$N_{\rm {HI}}\lesssim 10^{22} \, \rm {cm}^{-2}$$ for spectra with signal-to-noise ratios >4. Our results for DLA line density and total hydrogen density are consistent with previous measurements. Despite a small bias due to the poorly measured blue edges of the spectra, we demonstrate that our new model can measure the DLA population statistics when the DLA is in the Lyman-β forest region. We verify our results are not sensitive to the signal-to-noise ratios and redshifts of the background quasars although a residual correlation remains for detections from zQSO < 2.5, indicating some residual systematics when applying our models on very short spectra, where the SDSS spectral observing window only covers part of the Lyman-α forest. 
    more » « less
  3. Abstract Spatial aggregation of environmental or trophically transmitted parasites has the potential to influence host–parasite interactions. The distribution of parasites on hosts is one result of those interactions, and the role of spatial aggregation is unclear. We use a spatially explicit agent‐based model to determine how spatial aggregation of parasites influences the distribution of parasite burdens across a range of parasite densities and host recovery rates. Our model simulates the random movement of hosts across landscapes with varying spatial configurations of areas occupied by environmental parasites, allowing hosts to acquire parasites they encounter and subsequently lose them. When parasites are more spatially aggregated in the environment, the aggregation of parasite burdens on hosts is higher (i.e., more hosts with few parasites, fewer hosts with many parasites), but the effect is less pronounced at high parasite density and fast host recovery rates. In addition, the correlation between individual hosts' final parasite burdens and their cumulative parasite burdens (including lost parasites) is greater at higher levels of spatial parasite aggregation. Our work suggests that fine‐scale spatial patterns of parasites can play a strong role in shaping how hosts are parasitized, particularly when parasite density is low‐to‐moderate and recovery rates are slow. 
    more » « less
  4. Recent technology breakthroughs in spatially resolved transcriptomics (SRT) have enabled the comprehensive molecular characterization of cells whilst preserving their spatial and gene expression contexts. One of the fundamental questions in analyzing SRT data is the identification of spatially variable genes whose expressions display spatially correlated patterns. Existing approaches are built upon either the Gaussian process-based model, which relies onad hockernels, or the energy-based Ising model, which requires gene expression to be measured on a lattice grid. To overcome these potential limitations, we developed a generalized energy-based framework to model gene expression measured from imaging-based SRT platforms, accommodating the irregular spatial distribution of measured cells. Our Bayesian model applies a zero-inflated negative binomial mixture model to dichotomize the raw count data, reducing noise. Additionally, we incorporate a geostatistical mark interaction model with a generalized energy function, where the interaction parameter is used to identify the spatial pattern. Auxiliary variable MCMC algorithms were employed to sample from the posterior distribution with an intractable normalizing constant. We demonstrated the strength of our method on both simulated and real data. Our simulation study showed that our method captured various spatial patterns with high accuracy; moreover, analysis of a seqFISH dataset and a STARmap dataset established that our proposed method is able to identify genes with novel and strong spatial patterns. 
    more » « less
  5. Abstract Active galactic nuclei (AGNs) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) andSimbacosmological hydrodynamic simulations to show that the low-redshift Lyαforest can provide constraints on the impact of AGN feedback. We show that TNG overpredicts the number density of absorbers at column densitiesNHI< 1014cm−2compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducialSimbamodel, which includes AGN jet feedback, is the preferred fit to the observed CDD of thez= 0.1 Lyαforest across 5 orders of magnitude in column density. We show that theSimbaresults with jets produce a quantitatively better fit to the observational data than theSimbaresults without jets, even when the ultraviolet background is left as a free parameter. AGN jets inSimbaare high speed, collimated, weakly interacting with the interstellar medium (via brief hydrodynamic decoupling), and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the intergalactic medium when compared to TNG’s kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low-redshift Lyαforest provides plausible evidence for long-range AGN jet feedback. 
    more » « less