---
title: "在活动恒星旁，读懂地球质量信号"
description: "对照固定版本原论文，读联合活动推断、时间窗假峰与 GJ 1002，追踪协方差、仪器零点和后验条件。"
canonical_url: "https://sandboxphysics.com/zh/stories/rv-activity/earth-mass-signals-and-active-stars/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/rv-activity/earth-mass-signals-and-active-stars/"
last_modified: "2026-10-28"
---

[恒星活动与径向速度](https://sandboxphysics.com/zh/stories/rv-activity/)

AC18 · L4

# 在活动恒星旁，读懂地球质量信号

Sandbox Physics · 2026-10-28

每秒一米量级的信号很小，对恒星、仪器或观测时间窗的一项选择，可能和轨道一样影响结论。带着这些选择读论文，图就不只是几条平滑曲线。

## 开始前需要什么

已读表面与联合后验两章。这里导读固定版本原文，没有声称复现每篇论文的全部分析。

## 先问谱线模型分辨了什么

先读 SOAP 2.0 的 1409.3594v1，第 2.1—2.3 节和图 1，分清拿走光量与抑制对流蓝移。作者使用的局部线型信息比这里相同的高斯暗槽丰富。现场模型计算加权一阶矩机制，并没有继承 SOAP 对线平分线或互相关函数的全部刻画。

再读 Haywood 等人的 1407.1044v1，第 2—3 节和附表 A1，注意同期测光如何帮助约束恒星变化。公开重跑保留其中 2012 年数值表，采用的速度—钙线过程却是另一项有明确边界的分析。振幅恰好接近论文，也不能算复现了原文的测光条件与完整推断。

对应阅读：Dumusque, Boisse & Santos — SOAP 2.0, arXiv:1409.3594v1 · Haywood et al. — CoRoT-7, arXiv:1407.1044v1

## 顺着协方差，读到结果

在 Rajpaul 等人的 1506.07304v1，沿第 3.1—3.4 节，从共享潜过程走到多个观测量。先推附录 A.2，再看第 4.1—4.2 节和图 3—4。问清哪个指标响应过程本身，哪个响应导数，交叉协方差的符号才有落脚处。

还要区分“用模型自身生成”的合成检查与“用更丰富的物理模拟器生成”的检查。前者检查能否恢复自己的假设，后者才触及部分模型失配；两者都不能单独证明真实恒星变化已被包办。我们的高斯、多峰重要抽样基准检查积分，而只有星斑的物理注入检查另一类失败。

对应阅读：Rajpaul et al. — joint activity processes, arXiv:1506.07304v1 · Rasmussen & Williams — Gaussian Processes for Machine Learning, Chapter 2

## 把观测时间窗也放进论证

读《Ghost in the time series》的 1510.05598v1，第 2—3 节及图 1—3，留意保持采样、改变合成信号内容的比较。选中的周期残差能在扣除干扰后留下来，原因可能就在扣除方法和时间窗。这里没有重跑半人马座阿尔法 B 的数据；教程是在 CoRoT-7 的日期上做相关诊断。

把峰值当证据前，先交代搜过哪些频率、用了什么干扰模型、怎样选峰。最强峰上的条件区间，没有自动计入搜索过程。同样，若活动过程吸收了注入的行星，也该跨周期与先验查看恢复情况，不能只因残差变小而满意。

对应阅读：Rajpaul, Aigrain & Roberts — Ghost in the time series, arXiv:1510.05598v1 · Rajpaul et al. — joint activity processes, arXiv:1506.07304v1

## 同一仪器名称，也可能有几个状态

采用 Suárez Mascareño 等人的 2212.07332v2，重点读第 2.1—2.2、4.1—4.4 节、表 2、表 B.1 和图 9—10。139 次观测合并 CARMENES 与 ESPRESSO，后者又按硬件变化分成三个状态。先追踪各组独立零点与抖动，再去看轨道振幅。

这项研究已发表于《Astronomy & Astrophysics》670卷A5（2023），这里固定阅读修订后的 arXiv 作者稿。旧摘要把第二周期印成20.202天，修订版和参数表则为21.202天。先把数字与表格对上，就不会把一处笔误误读成新的物理问题。

对照原时间序列、GP 分量和扣除其他分量后的相位折叠。折叠曲线已经经过模型分配，不是原始观测；表 B.1 的先验也是模型的一部分。两个最小质量都约为一个地球质量，中心恒星却轻得多，周期约 10.35 和 21.20 天。地球—太阳例子的周期为一年，振幅只有每秒 0.089 米，需要区分这两种尺度。

真实教程只有一段 HARPS 观测，不能从中实证恢复那些多仪器零点。另设的零点跳变实验明确标为合成。倾角若只是采用值，把最小质量换成的也只是条件质量，不能因此变成实测倾角。

对应阅读：Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2 · GJ 1002 — published A&A 670, A5 (2023) · IAU 2015 Resolution B3 — nominal conversion constants

## 写下什么新观测会改变判断

每读一篇，记下观测跨度、独立通道、外部先验、仪器分组和合成检查，再挑一个最能帮助区分解释的新增测量：更长的相位覆盖、误差明确的活动线指标、同时测光，或独立天空轨道。值得做的新观测应让不同解释给出不同预测，而不只是往双方都能解释的位置再添一个点。

现场 Lab 给出表面机制，冻结观测给出实际条件后验，原论文展示更丰富的观测决策。把三个范围说清，才能来回比较，也不把好看的动画误认成实测恒星表面。

对应阅读：Haywood et al. — CoRoT-7, arXiv:1407.1044v1 · Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2

## 练习与解释性答案

1. 论文把活动和另一颗伴星扣掉后，画出很干净的行星折叠曲线。称它为原始发现图之前，应先看什么？

   展开解释

   找原时间序列、被扣掉的分量与后验条件。折叠图是在给定分配下展示某个分量，图很清楚，不等于扣除项的不确定性也消失了。

2. 仪器名称没变，为什么硬件更新后仍可能需要新零点？

   展开解释

   光纤等光学部件改变后，谱线形状和速度标定都可能改变。同一个名称不保证同一种响应。独立零点与抖动表达这种可能性，它们的需要和先验应结合实际仪器历史说明。

## 动手检验预测

### [计划一次能分辨解释的测量](https://sandboxphysics.com/zh/labs/binary-star-mass/)

1. 比较 Lab 的纯光量活动与明亮磁活动区预设。
2. 选择一个对两种解释响应不同的新增观测量，在把它加入似然之前，先写明所需的误差模型。

[上一篇：把这 71 次 HARPS 观测重新算一遍](https://sandboxphysics.com/zh/stories/rv-activity/replay-the-harps-campaign/)

## 来源与计算范围

1. [Dumusque, Boisse & Santos — SOAP 2.0, arXiv:1409.3594v1](https://arxiv.org/pdf/1409.3594v1)

   2014年9月11日版本。第2.1—2.3节、式1—7及图1：表面线型积分、光量对比与对流蓝移。这里相同的局部高斯线是更简单的独立机制模型。

2. [Haywood et al. — CoRoT-7, arXiv:1407.1044v1](https://arxiv.org/pdf/1407.1044v1)

   2014年7月3日版本。第2—3及4.2节、式1—12、附表A1。保留2012年HARPS全部71行数值，教程没有复现完整同期测光分析。

3. [Rajpaul et al. — joint activity processes, arXiv:1506.07304v1](https://arxiv.org/html/1506.07304v1)

   2015年6月24日版本。第3.1—3.4、4.1—4.2节，图3—4及附录A.2。导数符号采用第一时刻减第二时刻的时差；核宽度必须连同定义使用。

4. [Rajpaul, Aigrain & Roberts — Ghost in the time series, arXiv:1510.05598v1](https://arxiv.org/pdf/1510.05598v1)

   2015年10月20日版本。第2—3节和图1—3：采样窗口与合成信号内容的比较。半人马座阿尔法B分析属于原文导读，没有被说成本站复算的数据。

5. [Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2](https://arxiv.org/pdf/2212.07332v2)

   arXiv于2023年1月25日提交的版本，稿件页眉为1月26日。第2.1—2.2、4.1—4.4节，表2、表B.1和图9—10。重点读仪器状态与最小质量，没有转载原图。

6. [IAU 2015 Resolution B3 — nominal conversion constants](https://www.iau.org/common/Uploaded%20files/IAUGA2015-Resolution-B3-recommended-nominal-conversion.pdf)

   正式通过的表：名义太阳、地球引力参数和名义太阳半径是精确换算常数，不是没有测量误差的真实太阳与地球。

7. [Rasmussen & Williams — Gaussian Processes for Machine Learning, Chapter 2](https://gaussianprocess.org/gpml/chapters/RW2.pdf)

   作者提供的2006年第2章，第2.1—2.2及2.7节：高斯条件分布、正规边缘密度与线性基函数。这里的重要抽样器与表面积分为独立编写。

8. [GJ 1002 — published A&A 670, A5 (2023)](https://www.aanda.org/articles/aa/pdf/2023/02/aa44991-22.pdf)

   出版方版本及DOI 10.1051/0004-6361/202244991可核对发表状态。旧摘要保留20.202天笔误，固定的arXiv v2与参数表为21.202天；导读采用修订作者稿并说明版本区别。