低轨卫星轨道衰减速率与太阳活动指数 F10.7 的相关性：基于公开 TLE 的第 25 周检验

来源说明：本则出自独立撰写的第二批方案。它与前 20 则同样遵循八块结构与硬门槛要求，但撰写时未做英文文献检索，新颖性边界依据的是历届获奖图谱与既有知识，而非当轮查新。因此其「需核实」条目更多，启动前须自行补一轮英文检索。

1 · 研究问题

由公开 TLE(两行轨道根数)序列导出的低轨卫星半长轴衰减速率，与太阳 10.7 cm 射电流量指数 F10.7 的相关性有多强、滞后多少天？第 25 周极大期(2024 前后)衰减率相对静期的增幅，能否用标准热层密度模型半定量解释？

2 · 研究背景与空白

400–600 km 高度的卫星受稀薄热层大气拖曳持续降轨；热层密度对太阳 EUV 辐照极端敏感，而 F10.7 是 EUV 的经典代理指数。太阳活动上升会使同一卫星的衰减速率成倍变化——2022 年一批刚发射的 Starlink 卫星因地磁暴密度激增而再入，是该效应的著名实例。TLE 是北美防空司令部对全部编目天体持续发布的轨道根数，配合 SGP4 模型可还原平均半长轴的时间序列：公开、免费、更新至今。

热层密度-太阳活动关系有成熟的经验模型(NRLMSIS 系列)与大量文献；用 TLE 反演大气密度也是航天动力学的成熟技术。

空白在于时间覆盖：第 25 周极大期(2024–2025)是 20 年来最强的一段热层加热，用公开 TLE 对"多颗同类卫星 × 完整极大期"做同口径相关与滞后分析、并与密度模型对照的公开学生级工作很少(需核实：搜索关键词 "TLE orbital decay F10.7 solar cycle 25 thermosphere")。全部为文本数据，纯 CPU。

3 · 可检验假设

• H1：去趋势后的衰减速率与 F10.7(81 天平滑)的互相关峰值 ≥ 0.6，滞后在 0–5 天量级(热层响应时标)。
• H2：极大期(2024)平均衰减速率相对静期(2020–2021)的比值，与 NRLMSIS 模型在相应 F10.7 下的密度比一致到 2 倍以内；若系统偏离，偏离方向与地磁活动(Ap 指数)残余相关(标注为待检验推测)。

4 · 量化验收标准

1. 方法学校验(硬门槛):SGP4 传播链校验——用 ISS 公开 TLE 导出的高度序列与 NASA 公开的 ISS 轨道再推升(reboost)事件时间线对齐，全部推升事件(高度突增)识别率 100%，推升间隙的自然衰减段高度与公开报道值偏差 ≤ 2 km。此条不过关，后续全部结论无效。
2. 统计口径预先写死：时间序列一律按时间划分做任何拟合的训练/验证；相关性在双序列去趋势(21 天高通)后计算，显著性用块自助法(block bootstrap，块长 27 天≈太阳自转周期)≥1000 次；滞后由互相关峰 + 自助法置信区间。
3. 卫星样本 ≥ 10 颗：选无主动轨道维持的低轨目标(如已停止维持的退役卫星/火箭体,450–600 km)，筛选规则冻结；含机动的目标用突变检测剔除机动段(判据写死)。
4. 半长轴提取规则冻结:TLE 平均运动 → 平均半长轴，历元不均匀采样用逐历元差分 + 30 天滑动中位数。
5. H2 模型对照:NRLMSIS(Python pymsis 包，开源)在卫星高度、对应 F10.7/Ap 下的密度比；模型参数全部列表。
6. 可复现性:TLE 快照存档(下载即存，避免源变动)+ 全脚本开源。

5 · 数据与工具

用途	来源 / 工具
TLE 当前与近期	CelesTrak(celestrak.org，公开无需注册,GP 数据接口)
TLE 历史序列	Space-Track(space-track.org，免费注册；批量历史下载受 API 条款约束，需核实当前限速与允许用量；不可用时见 go/no-go 降级)
轨道传播	Python sgp4 包(官方算法移植，开源)；全部计算 CPU 秒-分钟级
太阳/地磁指数	F10.7：加拿大自然资源部/NOAA SWPC 公开日值;Ap：德国 GFZ 公开(kp.gfz-potsdam.de)
大气模型	pymsis(NRLMSIS 官方 Python 封装，开源),仅用于校验与对比，不计入本项目的数据贡献

6 · 方法路径

1. 装 sgp4/pymsis，取 ISS 历史 TLE 完成推升事件校验(验收 1)。
2. 冻结目标筛选规则，建立 ≥ 10 颗无维持目标清单并核实其状态(公开编目信息)。
3. 批量拉取 2020–2026 TLE，生成半长轴与衰减速率序列，机动剔除。
4. 与 F10.7 做去趋势互相关 + 块自助法显著性 + 滞后估计(H1)。
5. NRLMSIS 密度比对照极大/静期衰减比(H2),Ap 残差分析。
6. 核实工具边界:TLE 平均根数的类型(Brouwer 平均)与直接差分的系统差——读 sgp4 文档并用短弧数值传播对照一次。
7. 独立交叉校验：换一组不同高度(~500 vs ~800 km)目标，验证相关强度随高度降低而增强的预期。

7 · 新颖性边界

• 本课题不声称发现热层-太阳活动耦合(教科书级已知)，不声称改进大气模型;NRLMSIS 与 TLE 反演技术均为他人已发表工作。
• 已有工作:TLE 反演密度的方法学文献、极大期卫星再入的工程报告。
• 本项目贡献(主结论)：第 25 周完整极大期 × 多目标同口径的相关-滞后定量分析 + 与 NRLMSIS 的偏差表。属"时间覆盖空白 + 新样本独立检验"。
• 价值：相关强、滞后短且模型符合 → 公开数据即可监测热层响应的可复现范例；模型系统偏离 → 有文献价值的极大期偏差记录。两个结果都成立，块自助法保证显著性可信。
• 核对获奖清单：本赛道无轨道力学/空间天气课题；与 T07/T08 同属太阳影响链但观测端不同(X 射线事件 vs 大气拖曳)，三题互不重叠。

8 · 决策门槛(go / no-go)

⚠️ 第 4 周末：数据入口判定。阈值:Space-Track 历史 TLE 批量获取跑通(或确认受限)。若受限,降级路径一：用 CelesTrak 近实时数据前瞻式采集(第 4–40 周每周存档，累计 36 周仍覆盖一段活动变化)+ 公开的历史衰减案例(如 2022 Starlink 事件公开数据，需核实可得性);降级路径二：改用 ISS 等 NASA 直接公布轨道数据的少数目标做深入单星分析。两条路径保留"衰减-F10.7 相关"主框架。 ⚠️ 第 8 周末：验收 1 达标；不达标先查 TLE 历元解析与时区。 ⚠️ 第 20 周检查点：若互相关 < 0.3(H1 大幅失败)，检查目标是否有未公开机动；确认后把"TLE 序列中机动/异常检测方法"升级为并列结论，框架不散。 ⚠️ 选择前提：需接受"数据获取端存在条款不确定性"这一风险；不接受者选 T07/T10。

---