Recognition: unknown
A General Framework for Radial Velocity Calibration in Low-Resolution Spectroscopic Surveys: Correcting Wavelength-Dependent and Global Systematics with Application to LAMOST DR9
Pith reviewed 2026-05-10 02:27 UTC · model grok-4.3
The pith
A calibration framework fixes wavelength-dependent shifts and global offsets in LAMOST radial velocities, halving repeat scatter to 1.8 km/s.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors show that LAMOST low-resolution spectra exhibit both wavelength-dependent RV inconsistencies between spectral segments and global zero-point offsets. They correct the wavelength dependence by segmenting each spectrum into eight parts, measuring segment-wise offsets relative to the full spectrum at spectrograph and fiber levels, and fitting low-order polynomials to those offsets. Zero-point corrections are applied hierarchically: a joint chi-squared minimization at the spectrograph level constrained by repeat observations and cross-matches to APOGEE and Gaia RVS, followed by fiber-level averaging of seasonal offsets. The resulting catalog for 5.7 million spectra achieves a factor-
What carries the argument
Division of each spectrum into eight wavelength segments, measurement of segment-wise RV offsets relative to the full spectrum, low-order polynomial fits at spectrograph and fiber levels for wavelength dependence, and hierarchical chi-squared minimization plus seasonal averaging for global zero-point corrections.
If this is right
- Cross-night repeat RV differences at high signal-to-noise drop from 3.6 km/s to 1.8 km/s standard deviation, implying 1.3 km/s single-measurement precision.
- Dispersions against APOGEE and Gaia data fall from about 4.0 km/s to 2.0 km/s.
- A homogeneous value-added catalog of corrected radial velocities for approximately 5.7 million spectra is released.
- The same segment-polynomial and hierarchical zero-point procedure applies directly to radial-velocity calibration in other large-scale low-resolution spectroscopic surveys.
Where Pith is reading between the lines
- The cleaned LAMOST sample could support finer mapping of stellar motions across the Milky Way disk and halo than was previously possible with the uncorrected data.
- Adopting the framework in upcoming surveys with similar resolution and fiber-fed designs may allow their velocity catalogs to be combined with LAMOST without introducing large systematic mismatches.
- If low-order polynomials leave higher-frequency residual structure in some fibers, adding targeted higher-order terms or per-fiber spline corrections could be tested on repeat observations to check for further gains.
- The released catalog enables direct tests of whether the achieved precision is limited by photon noise or by uncorrected instrument effects at the current resolution.
Load-bearing premise
Reference radial velocities from external surveys are free of their own wavelength-dependent or zero-point systematics that could be absorbed into the LAMOST corrections.
What would settle it
A comparison of the corrected LAMOST velocities against radial velocities from a high-resolution survey that was not used in any part of the calibration process, checking whether the dispersion remains near 2 km/s.
Figures
read the original abstract
Radial velocity (RV) is crucial for stellar kinematics and Galactic archaeology. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) has obtained over ten million low-resolution spectra ($R \sim 1800$), yielding RVs for millions of stars, but these suffer from (1) wavelength-dependent inconsistencies (relative shifts between spectral segments) and (2) global zero-point offsets (uniform shifts of entire spectra). In this work, we comprehensively characterize and correct both. Each spectrum is first divided into eight segments of about 500 Angstrom. We organize the data at the spectrograph and fiber levels, measure segment-wise RV offsets relative to the full spectrum at each level, and then fit these offsets with low-order polynomials to correct wavelength-dependent systematics. We then correct zero-points hierarchically: at the spectrograph level by minimizing a joint chi-squared constrained by repeat observations and cross-matches with APOGEE and Gaia RVS, and at the fiber level by averaging seasonal offsets. After correction, RV precision improves significantly: for cross-night repeats, the standard deviation of RV differences at high signal-to-noise ratios drops by a factor of two from about 3.6 to about 1.8 km s$^{-1}$, implying a single-measurement precision of about 1.3 km s$^{-1}$. External checks with APOGEE and Gaia show dispersions drop from about 4.0 to about 2.0 km s$^{-1}$. The precision approaches, though slightly below, the theoretical limit at $R \sim 1800$. We release a value-added RV catalog with corrected velocities for about 5.7 million spectra, providing a homogeneous and systematically corrected dataset. The framework established in this work is also applicable to RV calibration in other large-scale spectroscopic surveys.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a general framework for correcting wavelength-dependent and global zero-point systematics in radial velocities from low-resolution LAMOST spectra. Spectra are divided into eight ~500 Å segments; segment-wise RV offsets are measured and fit with low-order polynomials at the spectrograph and fiber levels to address wavelength dependence. Global zero-points are then corrected hierarchically: spectrograph-level offsets via joint χ² minimization constrained by repeat observations and APOGEE/Gaia cross-matches, followed by fiber-level seasonal averaging. Post-correction, cross-night repeat RV difference SD at high S/N drops from ~3.6 to ~1.8 km s⁻¹ (implying ~1.3 km s⁻¹ single-measurement precision), and external dispersions with APOGEE/Gaia drop from ~4.0 to ~2.0 km s⁻¹. A value-added catalog of corrected RVs for ~5.7 million spectra is released.
Significance. If the corrections are robust and independently validated, the work delivers a substantially more homogeneous and precise RV dataset for Galactic archaeology and stellar kinematics, approaching the theoretical limit at R~1800. The hierarchical, multi-level approach and release of the catalog are practical strengths; the framework is also positioned as reusable for other surveys.
major comments (2)
- [§3.3] §3.3 (zero-point correction): The spectrograph-level zero-point offsets are obtained by minimizing a joint χ² explicitly constrained by the same repeat observations and APOGEE/Gaia cross-matches that are later used to quantify the post-correction scatter reductions (3.6→1.8 km s⁻¹ for repeats; 4.0→2.0 km s⁻¹ externally). This creates a risk that the headline precision gains are partly by construction. A held-out validation subset, temporal split, or explicit cross-validation procedure must be demonstrated to confirm that the reported improvements reflect removal of systematics rather than fitting to the validation data.
- [§3.2] §3.2 (wavelength-dependent correction): The selection of exactly eight segments and the specific polynomial orders are presented without quantitative justification or residual-structure diagnostics. It is unclear how these choices were optimized to avoid under- or over-fitting while fully capturing wavelength-dependent shifts; additional tests (e.g., segment-number sensitivity or post-fit residual RV maps) are needed to support that the polynomial model is sufficient and generalizable.
minor comments (2)
- [Abstract] The abstract and methods should explicitly state the exact wavelength boundaries of the eight segments and confirm that the ~500 Å figure is uniform across the LAMOST range.
- Figure captions for before/after RV comparison plots should include the exact S/N cuts and number of objects used in each panel to allow direct reproduction of the quoted SD values.
Simulated Author's Rebuttal
We thank the referee for the thoughtful and constructive comments. We address each major point below and will revise the manuscript accordingly to strengthen the validation and justification of our methods.
read point-by-point responses
-
Referee: [§3.3] §3.3 (zero-point correction): The spectrograph-level zero-point offsets are obtained by minimizing a joint χ² explicitly constrained by the same repeat observations and APOGEE/Gaia cross-matches that are later used to quantify the post-correction scatter reductions (3.6→1.8 km s⁻¹ for repeats; 4.0→2.0 km s⁻¹ externally). This creates a risk that the headline precision gains are partly by construction. A held-out validation subset, temporal split, or explicit cross-validation procedure must be demonstrated to confirm that the reported improvements reflect removal of systematics rather than fitting to the validation data.
Authors: We acknowledge the validity of this concern about potential circularity. The joint χ² minimization determines the offsets that best reconcile the repeat observations and external cross-matches, after which the scatter reductions are measured on the corrected values. While this is a common approach for deriving global corrections, we agree that independent validation is required to confirm the gains arise from systematic removal. In the revised manuscript, we will add a held-out validation analysis in §3.3: the repeat observations will be partitioned (e.g., 70/30 training/validation split or temporal split for cross-night data), the minimization performed only on the training subset, and the scatter reduction then evaluated on the held-out data. Analogous checks will be shown for the APOGEE/Gaia cross-matches. Results will be presented quantitatively to demonstrate that the improvements persist. revision: yes
-
Referee: [§3.2] §3.2 (wavelength-dependent correction): The selection of exactly eight segments and the specific polynomial orders are presented without quantitative justification or residual-structure diagnostics. It is unclear how these choices were optimized to avoid under- or over-fitting while fully capturing wavelength-dependent shifts; additional tests (e.g., segment-number sensitivity or post-fit residual RV maps) are needed to support that the polynomial model is sufficient and generalizable.
Authors: We agree that the manuscript would benefit from explicit justification and diagnostics for these choices. The eight ~500 Å segments were selected to partition the LAMOST wavelength range while ensuring sufficient spectral features and S/N per segment for stable RV measurements; low-order polynomials were adopted after inspecting the observed offset trends to capture smooth wavelength dependence. To address the referee's request, we will expand §3.2 with: (i) a sensitivity study varying segment number (e.g., 4, 8, 16) and reporting effects on final RV precision and residual structure; (ii) post-fit residual RV maps or binned diagnostics showing the reduction of wavelength-dependent trends; and (iii) a brief description of how polynomial orders were chosen (e.g., via residual minimization on a subset). These additions will support the robustness and generalizability of the model. revision: yes
Circularity Check
Zero-point corrections fitted via chi-squared to repeats and APOGEE/Gaia matches, with precision gains then reported on the same data
specific steps
-
fitted input called prediction
[Abstract]
"We then correct zero-points hierarchically: at the spectrograph level by minimizing a joint chi-squared constrained by repeat observations and cross-matches with APOGEE and Gaia RVS, and at the fiber level by averaging seasonal offsets. After correction, RV precision improves significantly: for cross-night repeats, the standard deviation of RV differences at high signal-to-noise ratios drops by a factor of two from about 3.6 to about 1.8 km s^{-1}, implying a single-measurement precision of about 1.3 km s^{-1}. External checks with APOGEE and Gaia show dispersions drop from about 4.0 to about "
The zero-point corrections are determined by minimizing chi-squared that incorporates the repeat observations and APOGEE/Gaia cross-matches as constraints. The reported post-correction reductions in RV scatter for cross-night repeats and dispersions versus APOGEE/Gaia are measured on these identical datasets, so the improvement is statistically enforced by the fitting process rather than an independent test.
full rationale
The wavelength-dependent polynomial corrections are derived from internal segment-to-full-spectrum offsets at spectrograph and fiber levels and appear independent of the target RV values. However, the global zero-point corrections minimize a joint chi-squared explicitly using the repeat observations and APOGEE/Gaia cross-matches. The headline precision claims (cross-night repeat SD dropping 3.6→1.8 km s^{-1}; external dispersions 4.0→2.0 km s^{-1}) are then evaluated on these same fitted datasets without held-out validation. This reduces the reported improvements partly by construction, matching the fitted-input-called-prediction pattern for the central validation step. No self-citations or other circular patterns are present.
Axiom & Free-Parameter Ledger
free parameters (2)
- number of spectral segments
- polynomial order
axioms (2)
- domain assumption Radial velocity offsets between spectral segments vary smoothly with wavelength and can be captured by low-order polynomials.
- domain assumption APOGEE and Gaia RVS radial velocities provide unbiased external anchors for zero-point calibration.
Reference graph
Works this paper leans on
-
[1]
2025, arXiv preprint arXiv:2503.14745 Bai, Z.-R., Zhang, H.-T., Yuan, H.-L., et al
Abdul-Karim, M., Adame, A., Aguado, D., et al. 2025, arXiv preprint arXiv:2503.14745 Bai, Z.-R., Zhang, H.-T., Yuan, H.-L., et al. 2017, Research in Astronomy and Astrophysics, 17, 091 Blomme, R., Fr´ emat, Y., Sartoretti, P., et al. 2023, A&A, 674, A7, doi: 10.1051/0004-6361/202243685 Bohlin, R. C., M´ esz´ aros, S., Fleming, S. W., et al. 2017, The Astr...
-
[2]
A3, but for Spectrograph
26 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=358) 9 11 1 3 5 2012–2013 (N=796) 9 11 1 3 5 2013–2014 (N=704) 9 11 1 3 5 2014–2015 (N=707) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[3]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=339) 9 11 1 3 5 2012–2013 (N=795) 9 11 1 3 5 2013–2014 (N=694) 9 11 1 3 5 2014–2015 (N=723) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[4]
A3, but for Spectrograph
27 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=341) 9 11 1 3 5 2012–2013 (N=786) 9 11 1 3 5 2013–2014 (N=714) 9 11 1 3 5 2014–2015 (N=723) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[5]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=338) 9 11 1 3 5 2012–2013 (N=776) 9 11 1 3 5 2013–2014 (N=698) 9 11 1 3 5 2014–2015 (N=702) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[6]
A3, but for Spectrograph
28 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=340) 9 11 1 3 5 2012–2013 (N=786) 9 11 1 3 5 2013–2014 (N=701) 9 11 1 3 5 2014–2015 (N=703) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[7]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=333) 9 11 1 3 5 2012–2013 (N=784) 9 11 1 3 5 2013–2014 (N=706) 9 11 1 3 5 2014–2015 (N=708) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[8]
A3, but for Spectrograph
29 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=349) 9 11 1 3 5 2012–2013 (N=777) 9 11 1 3 5 2013–2014 (N=705) 9 11 1 3 5 2014–2015 (N=709) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[9]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=336) 9 11 1 3 5 2012–2013 (N=768) 9 11 1 3 5 2013–2014 (N=688) 9 11 1 3 5 2014–2015 (N=686) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[10]
A3, but for Spectrograph
30 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=299) 9 11 1 3 5 2012–2013 (N=773) 9 11 1 3 5 2013–2014 (N=702) 9 11 1 3 5 2014–2015 (N=687) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[11]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=354) 9 11 1 3 5 2012–2013 (N=715) 9 11 1 3 5 2013–2014 (N=702) 9 11 1 3 5 2014–2015 (N=696) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[12]
A3, but for Spectrograph
31 1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=368) 9 11 1 3 5 2012–2013 (N=764) 9 11 1 3 5 2013–2014 (N=692) 9 11 1 3 5 2014–2015 (N=693) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 20...
2011
-
[13]
A3, but for Spectrograph
1012 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=359) 9 11 1 3 5 2012–2013 (N=788) 9 11 1 3 5 2013–2014 (N=712) 9 11 1 3 5 2014–2015 (N=716) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2015–...
2011
-
[14]
A3, but for Spectrograph
32 10 12 2 46 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2011–2012 (N=355) 9 11 1 3 5 2012–2013 (N=763) 9 11 1 3 5 2013–2014 (N=702) 9 11 1 3 5 2014–2015 (N=681) 9 11 1 3 5 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments 2...
2011
-
[15]
50 100 150 200 250 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–4400 ÅSpectral Segments Spectrograph 1 50 100 150 200 250 Spectrograph 2 50 100 150 200 250 Spectrograph 3 50 100 150 200 250 Spectrograph 4 50 100 150 200 250 8500–9000 Å 8000–8500 Å 6500–7000 Å 6000–6500 Å 5300–5800 Å 4900–5400 Å 4400–4900 Å 3900–...
2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.