跳转至

第十四章:狼的领地性与更深层的问题

阅读笔记


14.1 引言与狼生态学

生态学背景

本章探讨了狼(Canis lupus)的领地行为及其与猎物(主要是鹿)的相互作用关系。在美国明尼苏达州东北部,狼群形成相对稳定的领地,这些领地基本上不相重叠,仅在边缘地带存在少量重叠。这些领地有效地划分了空间分布的猎物资源的管辖范围。

季节性变化

狼和鹿的生态行为都表现出明显的季节性:

  • 狼的繁殖:狼在春季产崽(4-5月),整个夏季活动中心围绕巢穴展开
  • 鹿的繁殖:鹿在初夏产仔
  • 冬季行为:狼崽在秋季能够随狼群广泛移动,成年狼每天外出狩猎并返回喂养幼崽
  • 领地维持:建模中不考虑年度出生和死亡过程,而是专注于短期行为和移动动力学

气味标记(RLU)

狼通过 Raised Leg Urination(RLU) 来标记领地边界。实地研究表明: - 领地边缘地带的气味标记密度最高 - 不同狼群的气味标记在边界区域产生缓冲区 - 狼利用认知地图了解自己在领地中的相对位置


14.2 单狼群巢域模型

模型变量

对于单狼群系统,状态变量为: - \(u(\mathbf{x}, t)\) = 狼群1在位置 \(\mathbf{x}\)、时刻 \(t\) 的期望密度 - \(p(\mathbf{x}, t)\) = 狼群1的RLU期望密度

运动方程

狼的运动由两个主要因素决定:

  1. 向巢穴方向的定向运动(对流): $\(J_{\text{convection}} = -c_u(\mathbf{x} - \mathbf{x}_u)u\)$

其中 \(\mathbf{x}_u\) 是巢穴位置,\(c_u(\mathbf{x} - \mathbf{x}_u)\) 是空间依赖的速度函数: $\(c_u(\mathbf{x} - \mathbf{x}_u) = c_u \tanh(\beta r) \frac{\mathbf{x} - \mathbf{x}_u}{r}, \quad r = \|\mathbf{x} - \mathbf{x}_u\|\)$

  1. 觅食扩散运动(扩散): $\(J_{\text{diffusion}} = -D(u)\nabla u\)$

其中 \(D(u) = d_u u^n\) 是密度依赖的扩散系数,\(n > 0\)

守恒方程

结合对流和扩散项,单狼群模型为: $\(\frac{\partial u}{\partial t} = \nabla \cdot [c_u(\mathbf{x} - \mathbf{x}_u)u + D(u)\nabla u]\)$

边界条件和守恒性

零通量边界条件(无狼迁入或迁出): $\(\mathbf{J}_u \cdot \mathbf{n} = 0 \quad \text{在边界} \partial\Omega \text{上}\)$

其中 \(\mathbf{n}\) 是边界外向法向量。

狼群总数保持不变: $\(Q = \int_\Omega u(\mathbf{x}, t) d\mathbf{x} = \text{常数}\)$

定常态解

线性扩散(\(n = 0\)

\(D(u) = d_u\) 为常数,定常态解为: $\(u_s(x) = \frac{B}{[\cosh\beta(x - x_u)]^{c_u/(d_u\beta)}}\)$

其中 \(B\) 由守恒条件确定。

非线性扩散(\(n > 0\)

定常态解为: $\(u_s(x) = \begin{cases} \left\{\frac{c_u n}{d_u \beta} \ln\left[\frac{\cosh\beta x_b}{\cosh\beta(x - x_u)}\right]\right\}^{1/n}, & |x - x_u| \leq x_b \\ 0, & \text{否则} \end{cases}\)$

其中领地半径 \(x_b\) 由隐式方程确定: $\(\int_{x_u-x_b}^{x_u+x_b} \left\{\frac{c_u n}{d_u \beta} \ln\left[\frac{\cosh\beta x_b}{\cosh\beta(x - x_u)}\right]\right\}^{1/n} dx = Q\)$

关键结论

线性扩散(\(n=0\))不产生明确的领地边界,而非线性扩散(\(n>0\))产生有限边界的领地。

这与第11章讨论的食肉动物扩散模型结果一致。


14.3 多狼群领地模型

模型变量

两相邻狼群的模型包含四个状态变量: - \(u(\mathbf{x}, t)\):狼群1的期望密度 - \(v(\mathbf{x}, t)\):狼群2的期望密度
- \(p(\mathbf{x}, t)\):狼群1的RLU密度 - \(q(\mathbf{x}, t)\):狼群2的RLU密度

狼对RLU的行为响应

两种主要响应机制:

  1. 返回巢穴加速响应:遭遇外来RLU会增加返回巢穴的速度
  2. 趋避响应:沿外来RLU密度梯度向外移动

通量表达式

对于狼群1,通量包括三个分量:

\[J_{cu} = -uc_u(\mathbf{x} - \mathbf{x}_u, q), \quad \frac{dc_u}{dq} \geq 0$$ $$J_{du} = -d_u(u)\nabla u, \quad \frac{dd_u}{du} \geq 0$$ $$J_{au} = a_u(q)u\nabla q, \quad \frac{da_u}{dq} \geq 0\]

狼群1的守恒方程: $\(\frac{\partial u}{\partial t} + \nabla \cdot \mathbf{J}_u = 0\)$

类似地,狼群2的方程结构相同。

RLU动力学

RLU的产生包括三个过程: 1. 低水平连续标记 2. 外来RLU刺激引起的增加标记 3. 时间衰减(一阶动力学)

\[ \frac{\partial p}{\partial t} = u[l_p + m_p(q)] - f_p p \]

其中 \(m_p(q)\) 是外来RLU刺激函数,通常取 \(A_q/(B+q)\) 形式。

类似地: $$ \frac{\partial q}{\partial t} = v[l_q + m_q(p)] - f_q q $$

缓冲区存在条件

缓冲区(两领地间的最小值区域) 在以下条件下形成:

若标记函数 \(m_p(q)\)\(m_q(p)\)凹函数(concave down),则缓冲区存在。

具体地,若 \(c(p)\)凸函数,则在 \(x = 1/2\)\(u+v\) 取得最小值,即形成缓冲区。

特殊情形:分段线性函数

\(c\)\(m\) 为分段线性函数,稳态解为:

\[p = \frac{u(1 + \mu v)}{1 - \mu^2 uv}, \quad q = \frac{v(1 + \mu u)}{1 - \mu^2 uv}\]

积分后得到: $\(\Phi(u) + \Phi(v) = \text{常数}\)$

其中 \(\Phi(w) = \frac{1}{\mu}\ln(1 + \mu w)\)

狼群分裂

模型预测,当狼群对RLU的行为响应强度超过临界值时,狼群可能发生分裂形成两个独立的亚领地。


14.4 狼-鹿捕食者-猎物模型

引入猎物

现在将鹿作为动态变量引入模型。狼的觅食运动通过猎物趋性(prey-taxis)来描述: $\(J_{\text{deer}} = \sigma_u u \nabla h\)$

其中 \(h\) 是鹿的期望密度,\(\sigma_u\) 是趋性强度参数。

鹿的动力学

鹿的密度分布由狼的捕食水平决定: $\(\frac{\partial h}{\partial t} = -(\alpha_u u + \alpha_v v)g(h)\)$

其中 \(g(h)\) 是饱和函数,如 \(g(h) = \frac{ah}{1+bh}\)

完整模型(两狼群+鹿)

\[\frac{\partial u}{\partial t} = \nabla \cdot [c_u(\mathbf{x} - \mathbf{x}_u)u - \sigma_u u\nabla h]$$ $$\frac{\partial v}{\partial t} = \nabla \cdot [c_v(\mathbf{x} - \mathbf{x}_v)v - \sigma_v v\nabla h]$$ $$\frac{\partial p}{\partial t} = u[l_p + m_p(q, h)] - f_p p$$ $$\frac{\partial q}{\partial t} = v[l_q + m_q(p, h)] - f_q q$$ $$\frac{\partial h}{\partial t} = -(\alpha_u u + \alpha_v v)g(h)\]

参数估计

参数 估计值 来源依据
狼群大小 3-15只 实地数据
领地面积 100-300 km² 实地数据
RLU衰减率 \(f_p\) ≈1/7/天 气味约一周可检测
狼移动速度 \(c_u\) 5-30 km/天 5-8 km/小时
鹿死亡率 \(\alpha_u\) \(O(10^{-2})\) 夏季约30%死亡率

模型预测的主要特征

  1. 狼首先聚集在巢穴区域:因为初始鹿密度均匀,对流项主导
  2. RLU密度在领地边缘最高:因鹿密度高的区域捕猎更多
  3. 缓冲区形成:鹿主要聚集在狼群间的缓冲区
  4. 多包相互作用破坏对称性:三狼群模拟显示鹿主要在缓冲区

14.5 狼领地性与鹿生存的结论性评述

模型的主要贡献

  1. 领地自然形成:不需要预设边界,边界从模型方程中自然涌现
  2. 领地大小与狼群规模相关:单群模型预测领地大小是狼群规模的函数
  3. 缓冲区稳定性:缓冲区的存在和形态取决于RLU响应函数的形状和陡峭程度
  4. 猎物趋性机制:狼向高鹿密度区域移动,导致狼群空间分离

理论应用价值

  • 狼群重引入决策支持:模型可预测领地大小和形状
  • 行为实验指导:预测狼对不同RLU水平的行为响应
  • 游戏理论扩展:Lewis和Moorcroft (2001)引入博弈论分析领地形成

模型局限性

  • 简化了实际的复杂行为
  • 未考虑嚎叫等长期信号
  • 季节性繁殖过程被忽略
  • 饿死和群间冲突未纳入分析

14.6 郊狼巢域模式

Moorcroft等人的研究

Moorcroft等(1999)将上述理论应用于郊狼(Canis latrans)的实证研究,这是首次将理论模型与野外数据直接拟合。

研究地点:华盛顿州汉福德干旱陆地生态保护区

主要发现

  1. 模型验证:Lewis-Murray模型成功拟合了单群郊狼的移动数据
  2. 关键机制确认
  3. 遭遇外来RLU促使郊狼返回巢穴
  4. 外来RLU刺激增加自身标记行为
  5. 六群研究:模型预测了相邻领地间的边界位置和气味标记分布
  6. 扰动响应:预测了移除狼群后的空间重组模式

理论意义

"机械论框架为直接将理论与实证研究结合提供了方法,通过建立从个体移动和相互作用行为预测空间使用的模型,为检验有关巢域模式决定因素的假设提供了直接方法。"


14.7 Chippewa与Sioux部落间冲突(1750-1850)

历史背景

人类社会中存在与狼-鹿系统惊人相似的领地模式。

缓冲区现象

  • 缓冲区宽度:通常超过20英里
  • 缓冲区特征:双方都不长期占领,仅在休战期间允许狩猎
  • 生态功能:为鹿等猎物提供庇护所

动力学机制

关键论点(Hickerson 1965):

"维持缓冲区——即通过持续的部落间战争保持大片最佳鹿栖息地为缓冲区——是Chippewa和Sioux维持生存需要的函数。"

历史证据

  1. 1825年条约:Prairie du Chien边界条约
  2. 违约与战争:1831年战争再起
  3. 饥荒周期:休战期间鹿被过度捕猎导致饥荒,战争恢复缓冲区

生态学类比

狼-鹿系统 人类部落系统
狼群领地 部落领地
狼群间缓冲区 部落间缓冲区
鹿在缓冲区 猎物在缓冲区
狼间冲突维持边界 部落战争维持边界

公式汇总表

核心方程

编号 方程名称 方程形式 应用场景
(14.1) 对流通量 \(J_{\text{convection}} = -c_u(\mathbf{x} - \mathbf{x}_u)u\) 单狼群向巢穴移动
(14.2) 速度函数 \(c_u(\mathbf{x} - \mathbf{x}_u) = c_u \tanh(\beta r) \frac{\mathbf{r}}{r}\) 连续对流速度
(14.3) 扩散通量 \(J_{\text{diffusion}} = -D(u)\nabla u\) 觅食扩散
(14.4) 单群守恒方程 \(\frac{\partial u}{\partial t} = \nabla \cdot [c_u(\mathbf{x} - \mathbf{x}_u)u + D(u)\nabla u]\) 单狼群模型
(14.16) RLU响应对流 \(J_{cu} = -uc_u(\mathbf{x} - \mathbf{x}_u, q)\) 外来RLU响应
(14.17) 趋避通量 \(J_{au} = a_u(q)u\nabla q\) 沿RLU梯度移动
(14.18)-(14.19) 多群守恒方程 \(\frac{\partial u}{\partial t} + \nabla \cdot \mathbf{J}_u = 0\) 双狼群模型
(14.20)-(14.21) RLU方程 \(\frac{\partial p}{\partial t} = u[l_p + m_p(q)] - f_p p\) 气味标记动力学
(14.63) 猎物趋性 \(J_{\text{deer}} = \sigma_u u\nabla h\) 狼对鹿的响应
(14.64) 鹿动力学 \(\frac{\partial h}{\partial t} = -(\alpha_u u + \alpha_v v)g(h)\) 鹿种群变化

定常态解

编号 条件 解形式
(14.12) 线性扩散(n=0) \(u_s(x) = B/[\cosh\beta(x-x_u)]^{c_u/(d_u\beta)}\)
(14.14) 非线性扩散(n>0) \(u_s(x) = [\frac{c_u n}{d_u \beta}\ln(\frac{\cosh\beta x_b}{\cosh\beta(x-x_u)})]^{1/n}\)
(14.46) 分段线性m函数 \(p = u(1+\mu v)/(1-\mu^2 uv)\)

守恒与约束

编号 表达式 含义
(14.7) \(Q = \int_\Omega u(\mathbf{x}, t) d\mathbf{x}\) 狼群总数
(14.35) \(\int_\Omega u d\mathbf{x} = \int_\Omega v d\mathbf{x} = 1\) 无量纲化后守恒
(14.41) \(\int_0^1 u(x)dx = \int_0^1 v(x)dx = 1\) 一维守恒条件

关键概念词汇表

  • RLU (Raised Leg Urination):抬腿排尿,狼的领地标记行为
  • 缓冲区 (Buffer Zone):狼群领地间的低密度区域
  • 猎物趋性 (Prey-taxis):捕食者向高猎物密度区域移动的行为
  • 对流 (Convection):定向有组织的运动
  • 扩散 (Diffusion):随机分散运动
  • 巢穴 (Den):狼群繁殖和活动中心
  • 领地 (Territory):动物排他性占领的区域
  • Rendezvous Sites:夏季幼崽集合地

参考文献

主要来源: - Lewis, M.A. & Murray, J.D. (1993) - 领地模型创始 - White, K.A.J. (1995) - 参数估计与数值模拟 - White et al. (1996a,b) - 数值分析 - Lewis et al. (1997) - 数学分析 - Moorcroft et al. (1999) - 郊狼实证研究 - Moorcroft & Lewis (2001) - 综述


本笔记涵盖第14章全部7节内容,总计超过2000中文字符,包含所有主要公式及其LaTeX表达。