Morphological Defenses and Dynamic Resource Reallocation

This document provides the formal mathematical and biological formulation for plant morphological (constitutive) defenses and dynamic resource reallocation (apparent nutrition withdrawal) in the PHIDS ecosystem simulation model.

1. Biological and Ecological Context

Plants possess two distinct evolutionary defense paradigms against herbivory:

Constitutive (Morphological) Defenses: Permanent structural barriers such as trichomes, thorns, spines, and indigestible cell-wall compounds (lignin and silica). These defenses are metabolically fixed and directly impair herbivore feeding mechanics or reduce nutrient assimilation efficiency during every trophic attack.
Inducible Defenses & Dynamic Allocation: Physiological responses catalyzed by herbivore damage or volatile organic compound (VOC) warning signals. Upon perception, plants actively withdraw mobile nutrients and carbohydrates from vulnerable tissues (leaves and shoots) via vascular phloem transport into below-ground root sinks, temporarily lowering their nutritional attractiveness to foraging herbivores.

2. Mathematical Formulation

2.1 Constitutive Mechanical Attrition

Mechanical defenses (thorns and spines) inflict physical damage on grazing mouthparts during feeding events. Let \(m_{\text{bite}} \ge 0\) represent the species-specific mechanical damage coefficient per bite, and \(\rho_{\text{morph}} \in [0, 1]\) represent the herbivore's morphological resistance. The integer headcount reduction \(\Delta n\) for a swarm cohort of size \(n(t)\) feeding on plant tissue is given by:

\[\Delta n = \left\lfloor m_{\text{bite}} \cdot (1 - \rho_{\text{morph}}) \right\rfloor\]

\[n(t + \Delta t) = \max\left(0, n(t) - \Delta n\right)\]

Where \(\lfloor \cdot \rfloor\) denotes the floor function, enforcing discrete integer mortality within the herbivore population cohort.

Execution and Biological Rationale: Mechanical Attrition

Rather than artificially altering the core biochemical caloric density of the plant tissue, mechanical defenses such as trichomes, thorns, and rigid spines act exclusively to inflict localized physical trauma upon the grazing apparatus of the herbivore. This mandates a strict, immediate integer population penalty for every unit of biomass extracted, establishing an evolutionary gradient that heavily favors herbivore phenotypes possessing morphological counter-adaptations (\(\rho_{\text{morph}}\)).

Computationally, this dynamic is resolved as a non-allocating, \(O(1)\) integer subtraction evaluated inline within the synchronous trophic interaction loop. The strict enforcement of the mathematical floor function \(\lfloor \cdot \rfloor\) guarantees that fractional "ghost" casualties cannot propagate into the ECS state matrices, perfectly preserving integral population boundaries without triggering slow floating-point arithmetic traps.

2.2 Caloric Attenuation & Digestibility Discounting

Structural cell-wall barriers (lignin and silica) reduce the fraction of ingested plant biomass that herbivores can digest and metabolize into reproductive energy reserves. Let \(\mu_{\text{digest}} \in [0, 1]\) represent the plant's digestibility modifier and \(\delta_{\text{eff}} \ge 0\) represent the herbivore's digestive adaptation efficiency.

Given total plant energy consumed \(E_{\text{consumed}}\), the net metabolized energy \(E_{\text{metabolized}}\) added to the herbivore energy pool is:

\[\eta_{\text{net}} = \min\left(1.0, \max\left(0.0, \mu_{\text{digest}} \cdot \delta_{\text{eff}}\right)\right)\]

\[E_{\text{metabolized}} = E_{\text{consumed}} \cdot \eta_{\text{net}}\]

Execution and Biological Rationale: Caloric Attenuation

Integrating heavily cross-linked structural biopolymers, such as lignin and crystalline silica, directly into leaf cell walls drastically impairs the digestive mechanics of grazing swarms. While the raw biomass is successfully severed and ingested, the unyielding structural matrix prevents the herbivore's enzymatic pathways from extracting usable metabolic energy. This biological "Attrition Trap" effectively starves the herbivore cohort while their stomachs remain physically full, severely curtailing runaway reproductive cycles without requiring the plant to synthesize lethal, high-cost chemical toxins.

Within the simulation architecture, this attenuation is implemented as a mathematically bounded scalar multiplier \(\eta_{\text{net}} = \min(1.0, \max(0.0, \mu_{\text{digest}} \cdot \delta_{\text{eff}}))\). By intercepting the energy transfer pipeline immediately after gross biomass consumption is calculated, the engine decouples the physical removal of plant tissue from the actual caloric absorption of the swarm in a singular, vectorized arithmetic operation.

2.3 Rate-Limited Phloem Translocation Kinetics

When a plant activates a resource_withdrawal defense action, it commits to reducing its apparent nutritional factor \(N(t)\) toward a target factor \(N_{\text{target}} \in [0, 1]\) (where lower values represent lower nutritional attraction).

Because carbohydrate transport through vascular phloem sieve tubes is limited by hydrostatic pressure gradients and fluid viscosity, the transition is governed by a first-order rate-limited differential equation with translocation rate \(k_{\text{trans}} \in (0, 1]\):

\[\frac{d N(t)}{dt} = -k_{\text{trans}} \cdot \left(N(t) - N_{\text{target}}\right)\]

In discrete simulation time with tick step \(\Delta t = 1\):

Active Withdrawal Phase (\(\tau_{\text{withdrawal}} > 0\)):

\[N^{t+1} = N^t - k_{\text{trans}} \cdot \left(N^t - N_{\text{target}}\right)\]

Recovery Phase (\(\tau_{\text{withdrawal}} = 0\)):

\[N^{t+1} = N^t + k_{\text{trans}} \cdot \left(1.0 - N^t\right)\]

Execution and Biological Rationale: Rate-Limited Phloem Kinetics

Biological realities dictate that an autotroph cannot instantaneously evacuate complex carbohydrates and nitrogenous compounds from its vulnerable canopy tissues. Active vascular translocation is strictly rate-limited by internal hydrostatic pressure gradients and the inherent viscosity of the phloem sap (\(k_{\text{trans}}\)). This immutable physiological bottleneck enforces a temporal "window of vulnerability," during which incoming herbivores can continue to extract high-yield nutrition while the plant frantically attempts to bury its resources below ground.

To maintain continuous geometric stability within the environmental tensors, PHIDS approximates this fluid dynamic constraint via a discretized, first-order exponential relaxation recurrence: \(N^{t+1} = N^t + k (N_{\text{target}} - N^t)\). If the engine employed a naive, instantaneous step-function to modify nutritional values, it would introduce violent, discontinuous tears into the global flow-field potential matrix, inducing pathological routing errors and completely shattering the continuous biological tracking behaviors of the simulated swarms.

2.4 Attractant Field Gradient Scaling

The spatial attractiveness of a plant cell \((x, y)\) to foraging herbivores in the Numba-accelerated flow-field solver is scaled directly by the plant's current apparent nutrition factor \(N(x, y)\):

\[F(x, y) = \alpha \cdot \left(E_{\text{plant}}(x, y) \cdot N(x, y)\right) - \beta \cdot \sum_{k} T_k(x, y)\]

Where \(E_{\text{plant}}(x, y)\) is the total plant energy, \(T_k(x, y)\) represents localized repellent toxin concentrations, and \(\alpha, \beta\) are flow-field weighting coefficients. Down-regulating \(N(x, y)\) flattens the local attractant gradient, causing spatial flow-field navigation to steer approaching herbivore swarms away from the defended plant patch.

Dual Perspectives: Attractant Field Scaling

Biological Perspective: By translocating nutrients to root sinks, plants mask their nutritional value from sensory-seeking herbivores (trophic camouflage), forcing grazers to disperse toward other patches.
Computer Science & Mathematical Rationale: Multiplied directly into the 2D energy tensor \(E(x,y) \cdot N(x,y)\) prior to JIT flow-field generation. Modifies the global potential surface without requiring separate pathfinding re-computations or dynamic graph re-routing.

Appendix: Engine Implementation & Schema Mappings

For software engineers and data interface developers, this section maps the mathematical formulation above to concrete data schemas and system modules in the PHIDS codebase.

Schema Definitions (`src/phids/api/schemas/species.py` & `triggers.py`)

```python
class PassiveDefensesSchema(StrictBaseModel):
    """Constitutive morphological defenses of a flora species."""

    mechanical_damage_per_bite: float = Field(default=0.0, ge=0.0)
    digestibility_modifier: float = Field(default=1.0, ge=0.0, le=1.0)


class ResourceWithdrawalAction(StrictBaseModel):
    """Action configuring rate-limited phloem nutrient translocation."""

    type: Literal["resource_withdrawal"] = "resource_withdrawal"
    # Pydantic sentinel default: 1.0 = no-op (full apparent nutrition).
    # Explicit payloads must always supply the intended target factor.
    # The HTMX UI layer defaults to 0.2 for operator convenience.
    apparent_nutrition_factor: float = Field(default=1.0, ge=0.0, le=1.0)
    withdrawal_duration: int = Field(default=10, ge=1)
```

Component State (`src/phids/engine/components/plant.py`)

```python
@dataclass(slots=True)
class PlantComponent:
    apparent_nutrition_factor: float = 1.0
    target_nutrition_factor: float = 1.0
    translocation_rate: float = 0.2
    withdrawal_ticks_remaining: int = 0
```

System Execution Pipeline (`src/phids/engine/systems/lifecycle.py`)

```python
# Phloem translocation update during lifecycle tick
if plant.withdrawal_ticks_remaining > 0:
    plant.withdrawal_ticks_remaining -= 1
    plant.apparent_nutrition_factor += (
        plant.target_nutrition_factor - plant.apparent_nutrition_factor
    ) * plant.translocation_rate
else:
    plant.apparent_nutrition_factor += (
        1.0 - plant.apparent_nutrition_factor
    ) * plant.translocation_rate
```