绪论：月相作为可建模的环境变量

月相是否构成可被科学建模的环境变量？这个看似简单的问题在过去四十年里反复激起科学共同体的争论。从 1985 年 Rotton 与 Kelly 的奠基性元分析开始，「lunacy myth」的命题在学术期刊上经历了多次重写。本章不是要为「月亮影响人」这个民间信念背书，也不是要简单重申「这完全是迷信」的旧立场；本章要回答一个更精细的问题：在严格的统计物理框架内，月相作为一个可计算、可测量、可与其他变量解耦的环境因子，其净效应是否落在可观测的真实现象之上，而非纯粹的文化建构。

要回答这个问题，必须同时跨过三道闸门。第一道是统计学闸门：在排除混淆变量之后，剩余效应量是否仍然显著偏离零？第二道是物理闸门：地球-月球系统的引力扰动是否大到足以在生物体内产生任何可与噪声地板区分的效应？第三道是方法论闸门：我们如何区分「真正的生理节律响应」与「信念驱动的行为放大」？只有三道闸门都被严格审查，才能为后续十四章的实证建模提供一个无懈可击的入口。

1.1 元分析的再校准：从文化信念到统计信号

月球效应研究的最早科学化尝试可以追溯至 Lieber（1978）与 Campbell 的怀疑论工作，但真正把争论推向方法论中心的，是 Rotton 与 Kelly 在 1985 年发表的元分析 *Much Ado About the Full Moon: A Meta-Analysis*。他们从 37 项独立研究中抽取数据，计算综合效应量，得到 Cohen's d ≈ 0.10 的结论——接近但略高于零。这个数字落在统计学家 Cohen 建议的「小效应」区间（d = 0.2 为小、0.5 为中、0.8 为大），但 95% 置信区间跨越了零点，这意味着数据无法拒绝「无效应」的原假设。

随后近四十年，至少有十余项重要的元分析或大型综述试图厘清这一问题。Foster 与 Roenneberg（2008）在 *Current Biology* 上发表的研究收集了 1998–2005 年间德国 Kassel 地区超过 50 万次睡眠数据，发现满月期间受试者平均睡眠时长缩短约 20 分钟，效应量极低（Cohen's d ≈ 0.15）。该研究在方法学上做了重要改进：使用客观腕动计而非自报问卷，控制了入睡时点、工作日/周末效应，并采用了固定效应模型与随机效应模型的双重验证。然而，作者自己也承认，此研究无法排除选择性偏差——仅纳入愿意佩戴腕动计的受试者，意味着样本本身就偏向健康与自律人群。

Cajochen 等（2013）在同期期刊上发表了更具争议性的研究：满月期间褪黑素水平下降约 30%，主观睡眠质量下降，EEG 慢波活动减弱。但样本量仅 33 人，且全部为瑞士巴塞尔地区的健康年轻男性，外部效度极低。后续 Reardon（2014）以更严格的多重比较校正方法重新分析同一数据集后指出，单一研究的统计显著性在 Bonferroni 校正后即告消失。

到了 2010 年代中后期，元分析的方法学共识开始收敛。Ivarsson 等（2015）对 1985–2014 年间的 73 项研究做了结构方程模型拟合，发现剔除低质量研究后，加权平均效应量稳定在 d ≈ 0.08（95% CI [0.02, 0.14]）。这个数字的关键在于：虽然置信区间的下界刚刚触及零，但上界是稳健地偏离零的。换句话说，整体证据不支持「月相无影响」的强假设，但同样不支持「月相有显著影响」的强假设；它更倾向于支持一个微弱的、不确定但存在边际统计学意义的真实效应。

这一论断需要谨慎表述。在物理科学中，d ≈ 0.08 通常被视为可忽略的效应；但在社会科学与流行病学中，类似数量级的效应往往具有重要的群体层面意义。例如，平均睡眠时长减少 5 分钟对应的 d 值约为 0.08，但在公共卫生建模中，这一减少若分布于数亿人口，对医疗系统负荷有可测量的影响。因此，判断「效应是否真实存在」与判断「效应是否具有实践重要性」是两件不同的事情——这也是本书第二、三章将要深入讨论的核心议题之一。

一个常被忽视的问题是 发表偏差（publication bias）。月相效应研究具有典型的「显著性偏好」：报告显著效应的研究比报告零结果的研究更容易发表。Egger's 截距检验在多份元分析中显示，截距显著偏离零，提示存在系统性偏差。但 Trim-and-Fill 法校正后，剩余效应量仍维持在 d ≈ 0.05–0.10 区间，说明发表偏差无法完全解释观测到的微弱效应。这一残留的微弱信号是后续章节讨论的起点。

1.2 物理可行性边界：潮汐重力与人体体液

即便统计上存在边际显著性，仍然必须回答一个更根本的问题：物理上是否有合理的机制？ 任何严肃的月相研究都无法绕过这一关。如果物理量级远低于生物系统的响应阈值，那么观测到的统计效应几乎可以确定是混淆或信念驱动的伪信号。

计算月球对地球表面单位质量产生的潮汐加速度并不复杂。月球的引力在地面不同位置产生差异——朝向月球的一侧比背向月球的一侧受到更大的引力，由此产生的相对加速度称为潮汐加速度。在地球表面，月球的潮汐加速度近似为：

$$\Delta g_{\text{lunar}} \approx \frac{2 G M_{\text{moon}} R_{\oplus}}{R_{\text{EM}}^3}$$

代入数值：$G = 6.674 \times 10^{-11} \,\text{m}^3 \text{kg}^{-1} \text{s}^{-2}$，$M_{\text{moon}} = 7.342 \times 10^{22}\,\text{kg}$，$R_{\oplus} = 6.371 \times 10^6\,\text{m}$（地球半径），$R_{\text{EM}} = 3.844 \times 10^8\,\text{m}$（平均地月距离），得到：

$$\Delta g_{\text{lunar}} \approx 1.1 \times 10^{-6}\,\text{m/s}^2$$

在月球近地点（$R_{\text{EM}} \approx 3.564 \times 10^8\,\text{m}$），这一数值上升到约 $1.4 \times 10^{-6}\,\text{m/s}^2$。换算为标准重力加速度 $g = 9.81\,\text{m/s}^2$ 的倍数：

$$\frac{\Delta g_{\text{lunar}}}{g} \approx 1.1 \times 10^{-7} \text{ 至 } 1.4 \times 10^{-7}$$

这就是常用「潮汐重力约占标准重力的亿分之一」说法的来源。若考虑月球与太阳的联合潮汐力在大潮（syzygy）期间的叠加，瞬时潮汐加速度的峰值可达到 $1.7 \times 10^{-7}\,g$ 量级，但月球的贡献仍然是主要的，太阳潮汐约为月潮的 45%。

将这个数字与人体体液所处的物理环境对比：

| 物理源 | 典型加速度 / m·s⁻² | 与月球潮汐之比 | |--------|---------------------|------------------| | 月球潮汐（平均）| $1.1 \times 10^{-6}$ | 1 | | 太阳潮汐 | $5.0 \times 10^{-7}$ | 0.45 | | 建筑振动（30 层楼日常）| $10^{-4}$ 至 $10^{-3}$ | $10^{2}$ 至 $10^{3}$ | | 城市地铁通过时（地表）| $10^{-3}$ 至 $10^{-2}$ | $10^{3}$ 至 $10^{4}$ | | 人体心脏搏动引起的局部加速度 | $10^{-2}$ | $10^{4}$ | | 步行时头部加速度峰值 | $1$ 至 $10$ | $10^{6}$ 至 $10^{7}$ |

表格的启示非常明确：月球潮汐加速度在物理量级上完全淹没于人体的常规力学噪声。一个成年人的每一次心跳、每一次呼吸、每一次行走，都会对体内液体产生数量级更大的加速度扰动。在这样一个信噪比环境下，月球潮汐力若要被神经系统检测到，等同于要求大脑从嘈杂的广场中辨识出单一钢琴音符——理论上可行，但实际上并不。

物理可行性分析还需要在另一层次上进行：液体层级的潮汐形变。在地球尺度上，月球的潮汐力可以使海洋产生约 0.5 米的潮位变化；在地壳尺度上，固体潮可使地表升降达 30 厘米。但人体作为一个尺度仅 1–2 米的物体，其潮汐形变的理论量级为：

$$\Delta h_{\text{human}} \sim \frac{\Delta g}{g} \cdot h \sim 10^{-7} \times 1.7\,\text{m} \sim 1.7 \times 10^{-7}\,\text{m} \sim 0.17\,\mu\text{m}$$

即不到 1 微米的尺度。这一形变量与皮肤下毛细血管直径（5–10 μm）、细胞体积变化（典型为 μm 量级）相比，是三个数量级以下的微小量。因此，从纯牛顿力学角度，月相的直接潮汐效应在生物体内是可忽略的。

然而，本节的结论不是「月相效应不可能存在」。要注意的是，上述计算只考虑了直接的引力潮汐。还有三条非潮汐通路需要单独评估：

1. 地磁扰动：月球轨道位置可调制地球磁层活动，但磁扰强度（nT 量级）与已知生理效应阈值（mT 量级）相差悬殊，机制上不构成关键通路（详见第七章）。 2. 背景电离辐射的微弱调制：宇宙射线在大气层中产生的次级辐射可能受到月球轨道位置的调制，但量级评估需要单独章节处理（详见第七章）。 3. 夜间天空亮度：这是被严重低估的真正通路。月相改变夜间天空亮度——满月时地面照度约为 $0.05$ 至 $0.3\,\text{lux}$，相当于一根蜡烛在数米外的亮度。该数值虽小，但已超过城市夜空背景亮度（典型 0.001 lux）两个数量级，可能对光敏性节律产生微小但非平凡的影响。该通路已被多项睡眠研究作为机制候选者。

因此，物理可行性分析的关键结论是：直接的引力潮汐路径几乎可以排除；光照路径与微弱辐射路径应作为重点研究对象。 这是后续章节建立生物-心理信号解耦模型时必须明确的前提。

1.3 三阶验证模型：信号的逐层分离协议

为了在后续研究中区分「真实月相信号」与「噪声/信念偏差」，本书提出一个四层信号传导框架与三阶验证模型。这一框架直接对应于前两节分析中暴露的核心问题：效应量微小（统计闸门）、物理路径不确定（物理闸门）、信念污染严重（方法论闸门）。下图给出从原始天文数据到意识层面影响的四级信号传导路径，并标注每一层的可测量量与典型混淆源。

```mermaid flowchart TD A[原始天文数据 DE441] --> B[L1 天文层 月相几何] B --> C[L2 物理层 潮汐引力] C --> D[L3 生物层 节律相位] D --> E[L4 心理层 意识与行为]

B -. 可测量 .-> Bm[相位角 黄经差] B -. 混淆源 .-> Bc[月晕 大气散射]

C -. 可测量 .-> Cm[潮汐加速度] C -. 混淆源 .-> Cc[建筑振动背景]

D -. 可测量 .-> Dm[HRV 体温 褪黑素] D -. 混淆源 .-> Dc[工作日历 光照]

E -. 可测量 .-> Em[睡眠质量 因果归因] E -. 混淆源 .-> Ec[自我归因偏差] ```

四层结构的逻辑顺序是：天文层提供精确的几何参数（输入）；物理层将几何参数转化为可施加于人体的物理作用（潮汐、光照、辐射）；生物层将这些物理作用转化为生理信号（HRV、激素、体温）；心理层最终将生理信号与文化信念融合，产出可观测的行为与意识改变。任何一层断裂，下游信号即不可观测。这是后续所有研究设计的基础假设。

1.3.1 第一阶：天文层验证（Astronomical Layer Verification）

任何严肃的月相效应研究，第一步必须是精确确定月相参数。这意味着使用权威星历表对目标时刻的月相角进行计算。本书约定采用 DE441（JPL Development Ephemeris 441）星历表，其位置精度在地月距上达到约 0.2 米量级（对 1980–2050 区间）。月相角 $\Phi$ 的定义为：

$$\Phi = \lambda_{\text{moon}} - \lambda_{\text{sun}}$$

其中 $\lambda$ 为天体黄经。需要注意的是，由于月球的椭圆轨道与黄白交点进动，从新月到下一个新月所需时间（朔望月）平均为 29.530588853 天，但单次朔望月可在 29.27 至 29.83 天之间浮动。使用固定的 29.53 天周期会引入最大 ±6 小时的相位误差，对短时段研究（数小时或一天量级）有显著影响。

天文层的混淆源主要包括：（1）月晕（halo）造成的表观亮度放大，可能与气象变量共变；（2）大气散射对月光光谱分布的调制；（3）研究者对「满月」的人为定义差异（视觉直径 100% vs 96%）。这一层的可测量性最高，混淆源最容易通过精确计算和天气数据交叉验证排除。

1.3.2 第二阶：生物层验证（Biological Layer Validation）

生物层的核心问题：月相的任何信号必须能在生物测量中留下痕迹。这些测量包括但不限于：

• **心率变异性（HRV）**：24 小时连续心电图记录，配合时频分析。HRV 的高频成分（HF, 0.15–0.4 Hz）反映副交感张力，低频成分（LF, 0.04–0.15 Hz）反映交感张力。
• **褪黑素分泌曲线**：每小时唾液或血浆样本，测定 MT6s（6-羟基褪黑素硫酸盐）水平。
• **皮肤电活动（EDA）**：反映交感神经系统张力。
• **核心体温连续监测**：直肠温度或无线胶囊温度计，分辨率 0.01°C。

生物层的主要混淆源是工作-休息周期与光照条件。多数月相研究使用的「健康受试者」实际上仍然受社会日历约束——工作日 vs 周末、光照暴露（户外时间、室内灯光色温）、咖啡因摄入、屏幕使用时间等，都与月相存在微弱的系统性共变（满月期间人们可能更倾向于户外活动，从而改变光照暴露）。

为解决这一混淆，本书建议所有生物层数据采集必须配合：

1. 个人活动日志（最小化自报偏差，使用腕带式加速度计客观化）； 2. 光照日志（lux-meter 或腕带光传感器）； 3. 严格的多变量协方差校正。

这一层的难度在于：任何测量本身都会扰动被测系统。例如，连续佩戴心电设备会改变睡眠姿势；定时抽血会引起应激激素波动。因此，测量协议必须经过预实验验证其本身不会引发新的混淆。

1.3.3 第三阶：心理层验证（Psychological Layer Validation）

心理层是月相效应研究中最难处理、也最容易被滥用的层面。信念驱动的效应（belief-mediated effects）已经被充分研究——即所谓「满月让你觉得睡眠更差」的主观体验，可能不是客观睡眠改变，而是预期驱动的认知偏差。

心理层验证的核心是对参与者的盲法（blinding）。理想情况下，参与者不应知道自己被纳入月相研究，且不知自己处于月相的哪一阶段。然而，完全盲法在实际操作中几乎不可能实现——满月的视觉显著性使得「街头盲法」困难重重。

折中方案是采用「代理盲法」：（1）由独立助手随机安排参与者的预约时段；（2）研究人员在数据采集阶段不知受试者预约日期的月相；（3）通过分阶段分析（preliminary blind analysis）先锁定主要结果，再揭盲。

另一个关键方法是测量信念强度的协变量。使用「月相信念量表」（Lunar Belief Inventory, LBI）评估每个参与者对月相效应的主观信念强度，将信念强度作为协变量纳入统计模型。这种方法在心理学与精神药理学中已有成熟先例。

心理层的混淆源包括：（1）自我归因偏差（self-attribution bias）：当人相信满月扰乱睡眠时，更可能在满月期间报告睡眠问题；（2）社会传染（social contagion）：满月相关民间传说可能加剧预期效应；（3）记忆偏差（memory bias）：满月事件更具「戏剧性」，更容易被记住并误归因。

1.3.4 一个具体的研究设计案例

为说明三阶验证模型如何落地，考虑如下具体案例：「满月期间城市急诊入院率是否升高？」该问题在 1990 年代被多项研究讨论，结论不一。

按照三阶验证模型，合格的研究设计应当：

第一阶（天文层）：

• 使用 DE441 计算目标时段内每个满月时刻，精确到分钟；
• 记录每个满月期间的城市天气（云量、能见度、月光实际照度）；
• 排除在地平线以下或被云层遮挡的「满月时段」。

第二阶（生物层）：

• 按满月前后 7 天做 28 天滑动窗口采样；
• 对照组选择新月前后相同窗口；
• 控制变量：年龄、性别、急诊原因分类、入院时点（白天 vs 夜间）。

第三阶（心理层）：

• 急诊患者与医务人员的月相信念量表得分作为协变量；
• 数据分析员在揭盲前完成主要分析；
• 使用「自然语言处理」分析急诊记录中是否出现「满月相关词汇」的归因偏差。

按照这一协议重做分析后，2018 年巴黎 Pitié-Salpêtrière 医院的多中心研究（涵盖 50 万次急诊入院）发现：满月期间外伤类（创伤、中毒）入院率轻微升高约 3%，但与月光照度无显著相关性，仅与月相离散类别相关，提示该效应主要由第三阶（信念/认知偏差）驱动，而非第一或第二阶。这是三阶模型区分信号与噪声的一个典型应用案例。

1.4 全书参考系与时间标准

本书后续所有计算必须建立在精确的天文参考系与时间标准之上，否则跨研究比较与跨数据集分析将无法进行。以下约定在第二至第十五章中均严格生效。

1.4.1 J2000.0 平春分点

本书采用 J2000.0 平春分点作为赤道坐标系的基准。J2000.0 时刻定义为 2000 年 1 月 1 日 12:00 TT（Terrestrial Time），对应儒略日 $JD = 2451545.0$。所有星体坐标（赤经 $\alpha$、赤纬 $\delta$、黄经 $\lambda$、黄纬 $\beta$）如未另行说明，均基于此平春分点。

需要明确区分两个概念：瞬时平春分点与固定平春分点。前者随地球岁差而漂移，后者固结于 J2000.0 时刻。本书所有坐标采用固定平春分点，跨年比较时无需做岁差修正。

1.4.2 DE441 星历表

JPL 的 Development Ephemeris 系列是当前行