B18
城市噪声与鸟类鸣声最低频率的因果分离：以同一样点内「高噪声时段 vs 低噪声时段」的时间对照判定频率上移是否独立于栖息地结构
1 · 研究问题
在同一录音点位、同一物种上比较交通高峰时段与低噪声时段（清晨/节假日）的鸣声，最低频率（minimum frequency）是否随环境噪声即时上移？该时间尺度上的上移幅度，与传统「城市 vs 郊区」空间比较得到的幅度是否一致——若空间比较的幅度显著更大，差额可归因于栖息地结构等空间混杂因子而非噪声本身吗？
2 · 研究背景与空白
背景。鸟类在噪声环境中提高鸣声频率是行为生态学中被引用最多的现象之一，但其解释一直有争议：是适应性的频谱调整（主动避开低频交通噪声的掩蔽），还是唱得更响的物理副产物（Lombard 效应导致振幅升高、而振幅与频率在发声机制上耦合）。区分这两者需要控制振幅或控制空间混杂——而空间比较（城区 vs 郊区）天然混入了栖息地结构、植被密度、种群密度等一整套差异。
已有工作到哪一步。城市梯度上的频率上移已被反复证实并具体到物种：麻雀（Passer montanus）与长尾缝叶莺（Orthotomus sutorius）在城区的鸣唱频率均显著高于近郊；红耳鹎（Pycnonotus jocosus）在城区仅将鸣唱的第一个音节上移；欧乌鸫（Turdus merula）的城市梯度变异已有专文。机制侧也有明确进展：有工作论证高频鸣唱可能只是唱得更响的生理后果（Lombard 效应），但振幅升高并不能解释暗眼灯草鹀最低频率的升高。更重要的一条是，2021 年 Behavioral Ecology 的工作指出城市鸟类通信同时受非生物噪声与生物噪声约束，且在同时存在虫鸣时，麻雀与缝叶莺不再提高最低频率——说明这个响应是情境依赖的。栖息地结构与噪声共同塑造优势频率也已有专文。
空白在：绝大多数研究是跨样点的空间比较，把噪声与栖息地结构、种群密度绑在一起；在同一点位内利用时段差做时间对照、从而把空间混杂完全消掉的设计非常少，而这正是中学生用一部录音笔就能做的设计。这个空白适合一年期课题，但季节窗口是本题最大的死穴：鸣唱高峰集中在繁殖期，错过就要等一年。
3 · 可检验假设
- H1：同一点位同一物种，高噪声时段的鸣声最低频率显著高于低噪声时段，效应量对应每 10 dB(A) 噪声升高带来 ≥ 100 Hz 的上移，且混合模型中噪声固定效应的 95% CI 不含 0。
- H2（混杂量化备择）：跨样点空间比较得到的「每 10 dB 上移量」显著大于同点位时间对照得到的量（两者 95% CI 不重叠）；若成立，差额即为空间混杂因子（栖息地结构等）对表观噪声效应的贡献，可给出定量估计。
4 · 量化验收标准
- 方法学校验（硬门槛）：(a) 声压计与频率标定——用手机/录音笔录制已知频率的标准音（用另一设备播放 1、2、4 kHz 纯音），频谱峰值频率误差 ≤ 2%；用已知声压级参照（廉价校准声级计或播放定标信号）确认噪声测量的相对可比性，重复测量 CV ≤ 5%。(b) 最低频率测量的方法学陷阱必须显式处理——最低频率对信噪比与频谱图参数（窗长、动态范围阈值）高度敏感，须证明：在人工加噪的同一段录音上，所采用的测量流程给出的最低频率变化 ≤ 100 Hz；若变化超过阈值，必须改用对信噪比稳健的测量方式（如 −20 dB 相对阈值法或峰值频率），并在论文中把这一点作为方法学结论报告。此项不过关，后续全部频率数据无效——这是本课题最容易被评委击穿的地方。
- 采样规模写死：≥ 2 个物种 × ≥ 6 个固定点位 × 每点位 ≥ 2 个噪声时段 × 每时段 ≥ 5 个可用鸣唱片段；生物学重复 = 点位（近似个体/领域），n ≥ 6；技术重复 = 同一录音会话内的多个鸣唱片段。同一录音会话内的多个片段不得当作独立样本，必须以点位为随机效应。
- 每次录音同步记录：等效连续 A 计权声级 L_Aeq（1 min）、时间、气温、风速（> 3 级风不录）、与最近道路距离、点位周边 25 m 半径内的植被盖度（照片 + ImageJ 定量，判据写死）。
- 统计口径：以「点位」为随机截距的线性混合模型，固定效应为噪声级、时段、物种及其交互；报固定效应系数 + 95% CI 与边际/条件 R²，不得只报 p 值；同一模型再加入植被盖度作为协变量，比较噪声系数的变化幅度以量化混杂。空间比较（H2）用同一模型框架在跨点位数据上重跑，两个系数用自助法比较。
- 鸣唱片段的选取判据写死并盲法执行：信噪比阈值、片段完整性、不重复计入同一次连续鸣唱；由两名成员独立筛选片段，一致性 Cohen's κ ≥ 0.80；物种鉴定须双人独立确认，存疑片段剔除并记录剔除率。
- 全部原始音频、标注表、声学测量脚本（Praat 脚本或 Python/librosa）开源，第三方可重跑复现全部最低频率值（差异 ≤ 50 Hz）。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 录音 | 便携数字录音笔（线性 PCM，44.1 kHz/16 bit 以上）+ 防风罩（约 350 元）；或高质量手机录音 App（须核实是否有自动增益控制 AGC——AGC 必须可关闭，否则振幅信息全部作废） |
| 噪声测量 | 廉价数字声级计（A 计权，约 150 元）；须核实其与标准声级计的偏差，本课题只用相对值，绝对精度不作强要求但须在论文中说明 |
| 声学分析 | Praat（免费）或 Python + librosa/scipy（免费），纯 CPU；频谱图参数（窗长、重叠、动态范围）全程固定并写入方法 |
| 物种识别 | 现场目视 + 图鉴 + 鸣声比对；可用开源自动识别工具（如 BirdNET）作辅助，但每一条纳入分析的片段必须经人工确认，自动识别结果不得直接入库 |
| 环境协变量 | 手持风速计（约 80 元）、温湿度计；植被盖度用手机照片 + ImageJ |
| 道路/土地利用 | 公开地图服务的道路等级与距离量算；土地利用可用公开遥感产品（如 ESA WorldCover）需核实下载入口与分辨率 |
| 统计 | R（lme4、lmerTest、performance）或 Python（statsmodels），纯 CPU |
| 文献基准 | 麻雀/缝叶莺城市梯度频率数据、红耳鹎首音节上移、暗眼灯草鹀振幅–频率解耦、生物噪声调制效应等；仅用于对比与定位，不计入本项目数据贡献 |
预算合计约 500–700 元。
6 · 方法路径
- 建立录音与分析规程：固定录音增益（关闭 AGC）、固定录音距离范围、固定频谱图参数；完成标准音频率标定与噪声测量重复性校验。
- 完成最低频率的信噪比稳健性检验（第 1 条硬门槛 b）：对同一批录音人工加入不同水平的白噪声/交通噪声，量化测量流程的偏移，据此选定最终采用的频率测量口径并锁定不再更改。
- 踏查选点：在城市内选 ≥ 6 个点位，覆盖噪声梯度，且每个点位内部都存在明确的时段噪声差（例如临主干道公园入口）；记录点位环境参数。
- 采样：在繁殖期鸣唱高峰（当地通常 3–6 月，须先核实目标物种物候）按「同点位、同物种、不同噪声时段」配对录音，每点位每时段 ≥ 3 次会话分布于不同日期。
- 声学测量与筛选：盲法编号、双人筛选片段、双人确认物种，提取最低频率、峰值频率、频宽、持续时间。
- 拟合混合模型，估计时间尺度的噪声效应；再用同一批点位做跨样点空间比较，得到空间尺度的噪声效应，两者对比检验 H2。
- 独立交叉校验：在 2 个点位做一次播放实验的替代方案——被动噪声事件对照（利用突发交通噪声事件前后 30 s 的自然对照），检验个体在秒–分钟尺度上是否即时调整；若无法实施，改为在同一点位比较工作日与节假日的同一时段。
7 · 新颖性边界
本课题不声称发现噪声导致鸣声频率上移（这一现象已被反复证实），不做任何声音回放实验（playback）以免干扰繁殖期鸟类，不捕捉、标记或接触任何鸟类，因此无法追踪个体身份——「同一点位」只是个体的近似，这一局限必须在论文中显式声明并作为不确定性来源讨论。不声称区分了适应性调整与 Lombard 效应的最终机制。
已有工作完成了什么：城市梯度上的频率上移已在麻雀、长尾缝叶莺、红耳鹎（仅第一音节）、欧乌鸫等物种上定量报道；「高频可能只是唱得更响的生理后果」的 Lombard 解释与「振幅无法解释暗眼灯草鹀最低频率升高」的反证均已发表；2021 年 Behavioral Ecology 的工作证明生物噪声（虫鸣）在场时麻雀与缝叶莺不再提高最低频率；栖息地结构与噪声共同塑造优势频率已有专文。历届丘奖生物获奖论文中，2024 年 The Harker School 的《跨注意力多模态视听融合的蜜蜂健康评估系统》同样用声学做动物状态判断，但其目标是分类系统的性能，本课题的目标是效应量的因果分离，判据与产出物完全不同。
本项目的贡献（主结论）：用同一点位内的时间对照给出噪声效应的估计，并与传统空间对照的估计做定量比较，从而给出「空间比较高估了多少」的一个具体数字与置信区间。这是评价维度的转换——已有工作评估的是「有没有效应」，本项目评估的是「常用设计对该效应的偏置有多大」。
价值：如果时间对照的效应显著小于空间对照，说明大量文献中的效应量含有栖息地混杂；如果两者一致，则空间设计得到一次独立的方法学支持。两种结果都是有效结论，但必须用置信区间宽度证明本设计有能力分辨两者的差异。
8 · 决策门槛（go / no-go）
- 第 6 周末（2026 年 10 月中，季节性前置门）：必须完成录音设备标定与最低频率稳健性检验，并用秋季非繁殖期的零星鸣声/鸣叫（call）跑通一次完整流程。这一步必须在冬季之前完成，因为繁殖期窗口（次年 3–6 月）一旦开始就没有调试时间。若第 6 周未跑通，立即降级：把研究对象从「鸣唱（song）」改为全年可得的「鸣叫（call）」或改为城市昆虫（蟋蟀/蝉）的鸣声——后者季节窗口不同、密度高、易采样，分析框架与主结论完全保留。
- 第 10 周末：确认目标物种在本地的繁殖期物候与鸣唱高峰月份（查阅地方鸟类志与本地观鸟记录）。若目标物种在本地不常见，改选本地最常见的 2 个鸣唱物种，并在第 12 周前完成踏查定点。
- 第 12 周末：确认 ≥ 6 个点位同时满足「可安全到达 + 有明确时段噪声差 + 目标物种稳定出现」。若不足，降级：把设计从「时间对照 + 空间对照」缩减为纯时间对照（放弃 H2），把资源集中于把每点位的会话数提高到 ≥ 6。
- 第 28–40 周（2027 年 3–6 月，唯一的数据采集窗口）：这是不可延期的硬窗口，须在第 26 周前完成全部准备。窗口内每两周复盘一次可用片段数，若到第 34 周可用片段 < 目标量的 50%，立即增加点位或延长每次会话时长，不得等到窗口结束才发现数据不够。
- 第 44 周末：完成分析。若时间对照的噪声系数 CI 包含 0，以「在秒–小时尺度上未检出频率调整」为主结论成文，并把它与空间对照结果的对比作为核心论证——这恰恰支持了「空间比较的效应含混杂」这一命题，是强结论而非失败。
- 需提前核实（第 2 周内）：录音设备能否关闭 AGC（决定整个课题成立与否）；本地是否有观鸟/野生动物记录规范或保护区限制；学生在清晨时段野外活动的安全与监护安排。
- 预算裁剪顺序：先砍风速计（改用手机气象数据 + 现场定性判断，须写明）→ 再砍第二个物种（统计效力下降，主结论仍成立）。不可砍可关闭 AGC 的录音设备与声级计。
安全与伦理要点
- 全程为非侵入式被动观测：不捕捉、不标记、不接触、不投喂、不使用回放（playback）——回放会干扰繁殖期鸟类的领域行为，本方案明确排除。
- 录音距离保持在不引起警戒或飞离的范围（观察到警戒姿态或飞离即停止并后退），不接近巢区，不为拍摄或录音而搜寻鸟巢。
- 严格遵守《中华人民共和国野生动物保护法》与地方规定；不在自然保护区、禁入区域内作业；若目标物种属保护物种，仅做远距离被动录音，不做任何可能构成"干扰"的行为，必要时事先向当地林业/园林主管部门报备。
- 清晨野外作业须两人以上同行、告知教师与家长、避开无照明的偏僻路段与水边；不进入私人领地。
- 录音不可避免会采集到路人语音——须在方法中写明处理规则：不针对人声录音，含可识别对话的片段在分析后立即删除、不公开发布，公开的音频数据集须做人声段静音处理。
- 不记录车牌、不拍摄可识别人脸。