---
title: 'Spatial Ancestral Chains: Theory & Applications'
url: https://www.emergentmind.com/topics/spatial-ancestral-chains-sac
type: topic
---

# Spatial Ancestral Chains: Theory & Applications

Searching arXiv for the cited SAC-related papers and adjacent work to ground the article in current arXiv records.
I'll look up the arXiv entries directly to confirm the core papers and neighboring literature.
{"query":"1709.04416 Registering the evolutionary history in individual-based models of speciation", "max_results": 5}
Spatial Ancestral Chains (SAC) denote ancestor-tracing structures in which lineage is represented together with an explicit spatial embedding through time. In spatially explicit, individual-based simulations of speciation, a SAC for an extant individual is the ordered sequence of ancestor–descendant links, augmented by spatial coordinates at each generation, tracing the lineage back to a founder [1709.04416]. In spatial population models with local regulation, SAC refers to the genealogical path of a sampled individual through space-time, represented as a random walk in a dynamic random environment, with collections of such paths forming coalescing systems [1912.02558]. In monocular dynamic reconstruction, the same term is used in a different sense: a Spatial Ancestral Chain is the ancestor sequence of nodes in a temporal partition tree, queried recursively to compose spatial dynamics across hierarchy levels [2602.11845]. The shared structure is ancestry organized as a chain and coupled to spatial state, but the mathematical objects, observables, and inferential roles differ across these literatures.

## 1. Core definitions and formal structure

In the spatially explicit speciation setting, for individual $i$ extant at generation $t^*$, the ancestral chain is
$$
A_i = \{ i_{t^*}, p(i,t^*), p(p(i,t^*),t^*-1), \dots, \mathrm{founder} \},
$$
and the spatial SAC augments each node with coordinates $(x_{\cdot}(\tau), y_{\cdot}(\tau))$ for every generation $\tau$ on the path [1709.04416]. A single-lineage SAC is therefore a backward-time lineage record with explicit space-time localization. For two sampled individuals, SAC also supports a spatial coalescent representation: one traces both lineages back to the first common ancestor and extracts the two spatial paths from the extant samples to their MRCA [1709.04416].

In locally regulated spatial population models, SAC is operationalized as the spatial position of the ancestor $n$ generations back. The central construction is a process $X=(X_n)_{n\ge 0}$, where $X_n$ is the spatial location of the ancestor $n$ generations in the past of a sampled individual [1912.02558]. Conditional on the environment, $X$ is a time-inhomogeneous Markov chain. Multiple SACs then form a coalescing system of random walks in the same environment, and in one spatial dimension the diffusive scaling limit of the full collection is the Brownian web [1912.02558].

A precursor to this viewpoint appears in the branching Markov chain framework, where the lineage of a typical generation-$n$ individual is represented by a non-homogeneous Markov chain $Y^{(n)}$ with transition kernel
$$
Q_n(x,e,dy):=\frac{m_1(x,e,dy)\,m_{n-1}(y,Te)}{m_n(x,e)}.
$$
This is a Doob $h$-transform of the one-step mean kernel and functions as the ancestral chain of a typical descendant under a size-biased change of measure [1308.2835]. In homogeneous branching random walks, this chain reduces to the underlying spatial random walk [1308.2835]. This suggests that SAC in later spatial-population work can be viewed as a concrete, spatially embedded specialization of a broader ancestral-chain formalism.

In WorldTree, SAC is defined on tree nodes rather than biological individuals. For a temporal partition tree node $S_j(M_j,G_j)$, the chain
$$
C_j = \{S_{2^a}(M_{2^a},G_{2^a}) \mid k \in j/2^a,\ a = 1 \ldots \log_2 j\}
$$
collects ancestors along the root-to-node path, and the rendered representation at time $T$ is the union of the current node’s deformed Gaussians with those of all ancestors [2602.11845]. Here “ancestral” refers to hierarchy within a learned decomposition rather than genealogy in a population.

## 2. SAC in spatially explicit individual-based speciation models

The speciation framework is a spatially explicit, individual-based model with constant population size $N$, haploid genomes of length $B$ binary loci, sex labels in $\{\mathrm{male}, \mathrm{female}\}$, and coordinates $(x_i(t),y_i(t))$ on a periodic $L\times L$ lattice [1709.04416]. Distances are computed with the torus metric
$$
d_{ij}(t)=\sqrt{\min(|x_i-x_j|,L-|x_i-x_j|)^2+\min(|y_i-y_j|,L-|y_i-y_j|)^2}.
$$
Mating is local and compatibility-constrained: individual $i$ can mate only with partners $j$ satisfying $d_{ij}(t)\le S$ and genetic compatibility $h_{ij}(t)\le G$, where $h_{ij}(t)$ is the Hamming distance between binary genotypes [1709.04416]. Species at generation $t$ are defined as connected components of the genetic compatibility graph $G_t$, so conspecific individuals need not be directly compatible if they are connected through compatible intermediates [1709.04416].

Within this model, SAC sits between two bookkeeping schemes. The Most Recent Common Ancestor Time matrix (MRCAT) stores pairwise times to MRCA among extant individuals and is driven by parent pointers $p(\cdot,\cdot)$, while the Sequential Speciation and Extinction Events algorithm (SSEE) records forward-time species-level branching and extinction and yields the true species phylogeny [1709.04416]. SAC links these resolutions by using MRCAT to extract individual-level spatial lineages and coalescent paths, and by using SSEE to attach species-level spatial attributes such as ranges or centroids through time to branches traversed by those lineages [1709.04416].

For asexual models, the MRCAT recurrence is
$$
T_{t+1}(i,j)=T_t(P(i,t+1),P(j,t+1))+1,
$$
with $T_t(i,i)=0$ and $T_0(i,j)=1$ for $i\neq j$. For sexual models, maternal and paternal matrices evolve separately,
$$
T^M_{t+1}(i,j)=T^M_t(P_1(i,t+1),P_1(j,t+1))+1,\qquad
T^F_{t+1}(i,j)=T^F_t(P_2(i,t+1),P_2(j,t+1))+1,
$$
and a combined-lineage version may use
$$
T_{t+1}(i,j)=1+\min_{k,l\in\{1,2\}} T_t(P_k(i,t+1),P_l(j,t+1)).
$$
Coordinates $x(:,t),y(:,t)$ are stored per generation for SAC reconstruction [1709.04416].

The lineage reconstruction itself is direct. Starting from a focal individual at time $t^*$, one follows parent pointers backward, storing identifiers, generation stamps, and coordinates:
$$
\mathrm{SAC}_i=\{(A_i(g),\tau,x_{A_i(g)}(\tau),y_{A_i(g)}(\tau)):g=0,1,\dots\}.
$$
For two individuals, one can recover the MRCA from $T_{ij}(t^*)$, trace both parent chains until the first common ancestor, and extract the two coordinate sequences down to the coalescence time [1709.04416].

The model demonstrations use $N=1500$, $L=100$, $B=150$, $S=5$, $G=7$, $\mu=0.001$, and $Q=0.05$ [1709.04416]. Across 50 simulations at $T=500$ and $T=1000$ generations, trees from MRCAT were consistently closer to the true SSEE phylogeny than distance-based trees by Robinson–Foulds topology and $\gamma$ branch-length distributions, and the advantage increased at $T=1000$ [1709.04416]. Because SAC is reconstructed from parent pointers rather than from extant genetic distances alone, this result supports its use for fine-scale spatial genealogy in simulations.

## 3. SAC in locally regulated spatial population models

In spatial population models with local regulation, SAC is not primarily a storage device but a stochastic process induced by forward demographic dynamics. Two prototypical cases are emphasized: a discrete-time contact process interpreted through oriented percolation, and the logistic branching random walk (LBRW) [1912.02558].

For the contact process, the environment is an i.i.d. Bernoulli field $\omega(x,n)$ on $\mathbb{Z}^d\times\mathbb{Z}$. A site is open if $\omega(x,n)=1$ and closed otherwise. Occupancy evolves by open nearest-neighbor propagation, and local regulation arises because an inhabitable site chooses a single parent uniformly among occupied neighbors, so competition intensifies when many candidate parents are present [1912.02558]. In the stationary regime for $p>p_c$, the time-reversed indicator field $\xi_n(x)=1$ iff there is an infinite directed open path starting at $(x,n)$ serves as the dynamic environment for the ancestral walk [1912.02558]. On the event $B_0=\{\xi_0(0)=1\}$, the ancestral lineage from the origin moves by
$$
P(X_{n+1}=y\mid X_n=x,\xi)=\frac{1}{|U(x,n)\cap C|}
$$
for $(y,n+1)\in U(x,n)\cap C$, where $C=\{(x,n):\xi_n(x)=1\}$ is the oriented percolation backbone [1912.02558].

For the LBRW, $\eta_n(x)$ is the number of individuals at $x\in\mathbb{Z}^d$ in generation $n$, and each individual at $x$ produces offspring with Poisson mean
$$
\left(m-\sum_z \lambda_{z-x}\eta_n(z)\right)^+,
$$
followed by migration according to a symmetric, aperiodic finite-range kernel $p$ [1912.02558]. Local regulation is explicit in the competition term $\sum_z \lambda_{z-x}\eta_n(z)$, which depresses reproduction in crowded neighborhoods. In the stationary regime, the ancestral lineage of a uniformly sampled individual at $(0,0)$ has transition kernel
$$
P(X_{n+1}=y\mid X_n=x,\eta^{\mathrm{stat}})
=
\frac{
p_{x-y}\eta^{\mathrm{stat}}_{-n-1}(y)\left(m-\sum_z \lambda_{z-y}\eta^{\mathrm{stat}}_{-n-1}(z)\right)^+
}{
\sum_{y'} p_{x-y'}\eta^{\mathrm{stat}}_{-n-1}(y')\left(m-\sum_z \lambda_{z-y'}\eta^{\mathrm{stat}}_{-n-1}(z)\right)^+
}.
$$
The environment seen by the lineage is therefore the time-reversed stationary population field [1912.02558].

A central technical contribution is the regeneration construction. In the contact-process case, a local path exploration with additional randomness yields regeneration times
$$
T_0=0,\qquad
T_j=\inf\{k>T_{j-1}:\xi(\gamma_k^{(0,0)}(k))=1\},
$$
with increments $\tau_i=T_i-T_{i-1}$ and $Y_i=X_{T_i}-X_{T_{i-1}}$ such that $((Y_i,\tau_i))_{i\ge 1}$ are i.i.d., $Y_1$ is symmetrically distributed, and both $Y_1$ and $\tau_1$ have exponential tails [1912.02558]. In the LBRW, regeneration is obtained through coarse-grained good blocks, percolation comparison, strong coupling, and regeneration cones [1912.02558].

These renewal structures imply diffusive scaling. For the contact process, conditioning on $B_0$,
$$
P_\omega(X_n/n\to 0)=1
$$
for almost every environment, and both annealed and quenched central limit theorems hold with a centered, non-degenerate, isotropic Gaussian limit and variance
$$
\sigma^2(p)=\frac{E[Y_{1,1}^2]}{E[\tau_1]}\in(0,\infty).
$$
For the LBRW, under the small-competition regime ensuring a unique non-trivial stationary law, one has
$$
P(X_n/n\to 0\mid \eta_0^{\mathrm{stat}}(0)\neq 0)=1
$$
and an annealed CLT with a non-degenerate Gaussian limit [1912.02558]. The interpretation is that SAC behaves, at large scales, like a centered random walk whose effective variance encodes local regulation and environmental fluctuations.

Collections of SACs exhibit coalescence phenomena. In one spatial dimension, the entire set of ancestral paths on the oriented percolation backbone converges under diffusive scaling to the Brownian web, the continuum object of coalescing Brownian motions [1912.02558]. This is a system-level statement: the SAC ensemble, not only individual chains, has a universal scaling limit.

## 4. Pairwise coalescence, isolation by distance, and effective-size interpretations

The pairwise analysis of SAC is developed explicitly for two-dimensional LBRW lineages [2405.02090]. Two individuals are sampled at time $0$, one from site $0$ and one from site $x\in\mathbb{Z}^2\setminus\{0\}$, conditional on occupancy. Their ancestral locations are $X=(X_k)_{k\ge 0}$ and $X'=(X'_k)_{k\ge 0}$, and the coalescence time is
$$
\tau_{\mathrm{coal}}=\inf\{k\ge 0: X_k=X'_k \text{ and they choose the same parent}\}.
$$
Meeting time $\tau_{\mathrm{meet}}=\inf\{k\ge 0: X_k=X'_k\}$ satisfies $\tau_{\mathrm{meet}}\le \tau_{\mathrm{coal}}$ [2405.02090].

The main asymptotic statement is a two-dimensional tail law. Under $m\in(1,3)$, small competition parameters, and stationarity,
$$
\lim_{N\to\infty} P_{Nx}(\tau_{\mathrm{coal}}>N^{2\gamma})=\frac{1}{\gamma},\qquad \gamma\ge 1.
$$
An analogous result holds first for $\tau_{\mathrm{meet}}$, and finer analysis shows that $\tau_{\mathrm{coal}}-\tau_{\mathrm{meet}}$ is negligible at scale $N^{2\gamma}$ [2405.02090]. The result agrees with classical two-dimensional random-walk hitting asymptotics and with the suitably scaled stepping stone model [2405.02090].

The methodological core is a joint regeneration construction for the two lineages. Space-time is coarse-grained into boxes, regeneration times $R_i$ are defined for the pair, and the increments have polynomial tail control:
$$
\sup_{x,x'}P(R_i-R_{i-1}>t\mid X_0=x,X'_0=x')\le C t'^{-\beta}.
$$
A coupling to independent environments produces approximation bounds of order $\|x-x'\|^{-\beta}$ for the first regeneration increment, and a separation lemma gives exponential control on the time needed for the pair to separate to radius $n^\delta$ [2405.02090]. Hitting probabilities for spheres then reduce the large-distance genealogy to near-harmonic behavior of $\log \|x\|$ in two dimensions.

The same framework links pairwise SAC behavior to identity by descent. If the mutation rate is scaled as $\mu_N=mN^{-2\gamma}$, then
$$
\phi_N(Nx)=E_{Nx}[(1-\mu_N)^{2\tau_{\mathrm{coal}}}] \to 1-\frac{1}{\gamma},
$$
which matches the corresponding scaling of Malécot’s continuous-space approximation and the stepping stone model [2405.02090]. This supports the replacement of fluctuating local sizes by fixed effective sizes for large-scale inference in locally regulated populations [2405.02090]. In the terminology of SAC, pairwise links in the ancestral chain network have a predictable large-distance decay, and regeneration blocks provide the tractable decomposition that makes this asymptotic regime analyzable.

The spatially explicit speciation model gives related, but simulation-based, spatial diagnostics. Range overlap is
$$
O_{AB}(t)=\frac{|R_A(t)\cap R_B(t)|}{|R_A(t)\cup R_B(t)|},
$$
and isolation by distance can be quantified through a Mantel correlation between geographic and genetic distance matrices [1709.04416]. Centroid trajectories $c_s(t)$ align species movement with lineage chains and permit allopatric-versus-sympatric interpretation of branches traversed by SAC [1709.04416]. These are not asymptotic theorems, but they place SAC in a macroevolutionary spatial-signature analysis.

## 5. SAC in hierarchical dynamic scene reconstruction

In WorldTree, SAC is part of a monocular dynamic reconstruction framework that combines a Temporal Partition Tree (TPT) with a complementary spatial hierarchy [2602.11845]. The dynamic scene is represented by motion bases $M$ and canonical Gaussians $G$, with deformation driven by a KNN motion graph and blended $SE(3)$ motions via Dual Quaternion Blending [2602.11845]. TPT recursively partitions the time interval and assigns each node a subset of motion bases and Gaussians. A child focuses on a narrower temporal interval and inherits only a subset of its parent’s Gaussian primitives [2602.11845].

SAC supplies the missing spatial context created by this inheritance truncation. For node $S_j(M_j,G_j)$, the ancestral chain $C_j$ contains the node’s ancestors in TPT, and the deformed Gaussian set rendered at time $T$ is
$$
P_j=\{G_h(M_j,T)\mid G_h\in G_j\}
\cup
\left\{
\bigcup\left(\{G_h(M,T)\mid G_h\in G_k\}\mid S_k(M_k,G_k)\in C_j\right)
\right\}.
$$
Thus a node’s render is the union of its own deformed Gaussians and the deformed Gaussians from all ancestors on the chain [2602.11845]. SAC therefore performs composition rather than parameter fusion: ancestors remain separately optimized, and their contributions are aggregated only at render time [2602.11845].

The paper frames this as “hierarchical spatial composition and modeling specialization” and “hierarchical motion decoupling” [2602.11845]. Each ancestor has its own motion bases and Gaussians, optimized only on its own interval, while SAC recursively queries these ancestors to provide complementary spatial dynamics that were pruned during inheritance [2602.11845]. During breadth-first TPT construction, chains propagate by appending the current node to the parent chain:
$$
C_{2j}=\mathrm{Append}(C_j,S_j),\qquad C_{2j+1}=\mathrm{Append}(C_j,S_j).
$$
Nodes at the same depth have the same ancestor-chain initialization but are optimized independently, with frozen ancestors, enabling parallel depth-wise optimization [2602.11845].

WorldTree reports that SAC improves reconstruction quality beyond TPT alone. On NVIDIA-LS, enabling TPT on top of a baseline brings LPIPS from $0.121$ to $0.113$ and mAVGE from $0.109$ to $0.102$; adding SAC further improves LPIPS to $0.108$ and mAVGE to $0.096$ [2602.11845]. With BA and SW enabled, the final metrics are PSNR $24.06$, SSIM $0.804$, LPIPS $0.100$, AVGE $0.056$, mPSNR $18.55$, mSSIM $0.692$, and mAVGE $0.092$ [2602.11845]. On DyCheck, the reported mLPIPS is $0.240$ versus MoSca’s $0.264$, an improvement of $9.09\%$ relative to the second-best method [2602.11845].

Ablations emphasize chain length. With TPT depth fixed, the full chain yields mPSNR $18.69$, mSSIM $0.698$, and LPIPS $0.095$, whereas progressively reduced chains yield mPSNR $17.44/16.55/15.54$, mSSIM $0.627/0.581/0.521$, and LPIPS $0.117/0.128/0.143$ [2602.11845]. This indicates that the ancestor composition itself, rather than merely hierarchical splitting, is essential in that setting.

## 6. Algorithms, complexity, limitations, and cross-domain significance

The algorithmic profile of SAC depends sharply on domain. In the speciation model, MRCAT requires storing pairwise genealogical times over extant individuals, with $O(N^2)$ update time and $O(N^2)$ storage per generation; full parent-pointer storage for SAC reconstruction is $O(NT)$ in memory, while checkpointing reduces this to $O(NT/k)$ at an $O(k)$ extraction penalty per query [1709.04416]. Sparse representations are generally not useful for $T$ because most entries are nonzero, so scaling strategies rely instead on downsampling, windowing, and checkpointing [1709.04416]. SSEE preserves extinct species branches, but individual-level SACs for extinct individuals require archive storage if they are to remain queryable [1709.04416].

In the local-regulation literature, the main limitation is analytic rather than storage-related. The environments are non-elliptic and, in the contact-process case, not Markovian in a local sense, so general random-walk-in-dynamic-random-environment results do not directly apply [1912.02558]. Model-specific regeneration constructions are therefore essential. Open problems include characterizing $\sigma^2(p)$ as $p\downarrow p_c$ for the contact process, establishing a quenched CLT for the LBRW, and understanding higher-dimensional scaling limits of coalescing systems, since no Brownian-web analogue exists in $d\ge 2$ where Brownian motions do not meet [1912.02558]. In two-dimensional LBRW, the coalescence asymptotics require stationarity, small competition, and finite-range dispersal and competition; the results identify the decay profile but not the prefactors needed for direct numerical prediction [2405.02090].

In WorldTree, SAC adds a union over ancestors during rendering, so rasterization cost grows with the number of Gaussian primitives included from ancestral nodes [2602.11845]. Practical efficiency is maintained by limiting tree depth and inherited Gaussian counts. The experiments use TPT depth $2$, inherited Gaussian counts of $5000$ at layer $1$ and $2000$ at layer $2$, BA for $2000$ steps, Static Warm-up for $6000$ steps at the root, and layerwise optimization schedules of $4000/2000/2000$ steps for NVIDIA-LS and $4000/1000/1000$ for DyCheck [2602.11845]. The method depends on pretrained depth, tracking, and flow priors; failures are reported on small, fast-moving regions and under pervasive scene-wide non-rigid deformation [2602.11845].

Across domains, SAC consistently serves as an ancestry-conditioned spatial representation, but its formal role changes. In branching and population models, SAC is a stochastic lineage object used to derive many-to-one formulas, coalescence laws, CLTs, and scaling limits [1308.2835; 1912.02558; 2405.02090]. In speciation simulations, SAC is a reconstructible data structure that connects individual genealogy, space, and species-level phylogeny [1709.04416]. In dynamic reconstruction, SAC is a render-time composition mechanism over hierarchical ancestors that restores spatial support without re-entangling optimization [2602.11845]. This suggests that the unifying content of the term is not a single biological or computational ontology, but a common design principle: spatial information becomes tractable when ancestry is organized as an explicit chain and queried across scales.

Source: https://www.emergentmind.com/topics/spatial-ancestral-chains-sac