地表气温长程相关性 DFA 标度指数的新时间窗独立检验与"站点 vs ERA5 再分析"同点偏置量化

1 · 研究问题

用 1991–2025 年新 35 年时间窗与自选全球站网重新估计日气温距平的去趋势涨落分析（detrended fluctuation analysis, DFA）标度指数 α，"内陆站 α ≈ 0.65、海岸/岛屿更高"的经典分区结论是否依然成立？在同一地点上，ERA5 再分析格点序列的 α 是否相对实测站存在系统偏移？

2 · 研究背景与空白

气温距平的长程相关性表现为 DFA 涨落函数 F(s) ~ s^α 且 α > 0.5，被视为气候系统慢分量（海洋、土壤记忆）的统计指纹。Koscielny-Bunde et al. (1998, PRL) 提出内陆站 α ≈ 0.65 的"普适律"；Fraedrich & Blender 及后续工作（如中国 557 站研究、澳大利亚地理依赖研究、希腊 ERA-Interim 的 MF-DFA）确认了海陆差异与纬度依赖，也对普适性提出修正。方法与数据全公开。

空白在时间覆盖与评价维度：其一，经典结论的数据多止于 2010 年前后，1991–2025 这一增暖加速期的独立时间窗检验（趋势去除方法对 α 的影响也需重估）是标准的换样本轴；其二，"同一位置上站点观测与 ERA5 的 α 是否系统不同"直接量化再分析同化平滑对标度统计的影响——本次检索未找到完全相同的系统同点对照工作（未找到不等于不存在，第 4 周文献复查后冻结表述）。计算量几乎为零，正负结论都有效。

3 · 可检验假设

• H1：内陆站 α 分布中位数落在 0.60–0.70，海岸/岛屿站显著更高（Mann–Whitney U 检验 p < 0.01）。
• H2：ERA5 同点序列的 α 系统性高于配对站点值（再分析平滑抬升相关性），配对差中位数 > 0.02 且 Wilcoxon 符号秩 p < 0.01；若无系统差，则"再分析可替代站点做标度分析"获得一次定量支持——两种结果都有效。

4 · 量化验收标准

1. 方法学校验（硬门槛）：自写 DFA 对合成分数高斯噪声（H = 0.6/0.7/0.8 各 100 条、长度与真实序列相同）恢复 α 平均偏差 ≤ 0.02；并对 3 个文献明确给出 α 值的站点复现至偏差 ≤ 0.05。此步不过关，后续全部结论无效。
2. 以 DFA2 为主口径，DFA1/DFA3 做敏感性；拟合区间预先写死为 10 天–4 年并报告区间端点敏感性。
3. 站点质控口径预先写死：记录长度 ≥ 30 年、缺测 < 10%、GHCN 质量标志过滤规则全部写入脚本。
4. 季节循环去除双方案对照（逐日气候态扣除 vs 前 4 阶傅里叶谐波），两方案结论分别报告。
5. 所有 α 的误差棒用滑动窗 + bootstrap 给出；站网选择判据（分层随机抽样按大洲/海陆/纬度带）预注册。
6. 全部脚本与站点清单开源。

5 · 数据与工具

用途	来源 / 工具
站点日气温	GHCN-Daily（NOAA NCEI，ncei.noaa.gov 按站号直接下载 CSV，单站约 1–5 MB；覆盖全球十万余站、逐日、1880s 至今），目标站网 100–200 站
再分析	ERA5 单层 2 m 气温（Copernicus CDS，cds.climate.copernicus.eu，需免费注册；用 CDS API 按点/按年切片下载，单点 35 年小时序列 ~ 数十 MB，先日均化；0.25° 分辨率、1940 至今）。CDS 新版接口与配额需第 6 周核实
合成校验信号	自写分数高斯噪声生成器（numpy FFT 法）——仅作校验
文献 α 值	Koscielny-Bunde 1998 等表列值——仅作对照，不计入本项目的数据贡献
分析	numpy/scipy；DFA 自写（约 60 行，不依赖黑箱包，便于审查）

算力：全项目任何笔记本数小时内完成，无算力风险。降规模版本：站网 200 → 60 站、ERA5 只取配对子集 → H1/H2 框架完整保留，仅统计力下降。

6 · 方法路径

1. 自写 DFA 与 fGn 生成器，完成硬门槛校验（合成信号 + 文献站点复现）。
2. 按预注册判据抽取全球站网，写死质控脚本并输出站点清单。
3. 计算全部站点 α（DFA1/2/3、双去季节方案），做海陆/纬度分组统计检验（H1）。
4. 核实并跑通 CDS API，按站点坐标下载 ERA5 同点序列并日均化。
5. 做配对 α 差分析（H2），并按纬度带/地形分层看偏置结构。
6. 敏感性总表：拟合区间、DFA 阶数、去季节方案三因素对主结论的影响矩阵。
7. 独立交叉校验：对 20% 抽样站点用功率谱指数 β（α = (β+1)/2 关系）复核 DFA 结果。

7 · 新颖性边界

• 本课题不声称发现气温长程相关性（1998 年已发表）也不声称提出新统计方法。
• 已有工作：Koscielny-Bunde et al. (1998) 与后续区域研究（中国、澳大利亚、希腊等）已给出 α 的地理分布图景；MF-DFA 应用于再分析数据也已存在。
• 本项目贡献（主结论）：1991–2025 新时间窗的全球分层站网独立检验 + 站点/ERA5 同点配对偏置的系统量化。前者是最稳的换样本轴；后者是评价维度补充，对所有用再分析数据做标度分析的研究是直接有用的定标信息。
• 若分区差异不显著，"经典分区在新窗口不显著"同样是有效结论，需用误差棒证明可分辨 0.05 量级的 α 差。
• 历届获奖图谱中无大气统计类先例；2022 有一项微塑料–沉积物与气候相关的提名奖，主题无重叠。

8 · 决策门槛（go / no-go）

• 第 4 周末：硬门槛（合成信号 ≤ 0.02）通过；不过关 → 排查 DFA 实现。
• 第 6 周末：核实 CDS API 注册、配额与下载速度（不要边做边发现）。若 ERA5 获取受阻 → 具名降级：H2 改为"α 对海拔与离海距离的回归结构"（纯 GHCN 即可），H1 与全部统计框架保留。
• 第 16 周末：站网 α 全表完成。若海陆差异效应量 < 预注册最小可检测差 → 扩站网至 300 站（一周预算）或正式接受"不显著"结论并转入敏感性归因。
• 预算裁剪顺序：ERA5 配对范围 → DFA 阶数敏感性 → 站网规模下限 60 站不可再裁。
• 适合学生：统计兴趣强、编程入门即可的队伍；本题是十题中"每周都有产出"最有保障的一条。

---