REVIEW 4 major objections 5 minor 1 cited by
Distributed Acoustic Sensing for Environmental Monitoring, and Newtonian Noise Mitigation:Comparable Sensitivity to Seismometers
T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Distributed acoustic sensing, using ordinary optical fiber, cancels a co-located seismometer's noise to the same residual factor as a geophone array (0.11 at 20 Hz), making it a scalable witness-sensor option for gravitational wave…
desk verdict First DAS-based cancellation of a colocated seismometer's noise is a real experimental step, but the headline residual is in-sample and internally inconsistent. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The mechanism is the equivalence between DAS strain-rate and ground velocity: $d\varepsilon/dt = [v(x+L_g/2)-v(x-L_g/2)]/L_g$, with the conversion $\varepsilon = du/dx = \pm (1/s)\,du/dt$, where $s$ is the local apparent slowness estimated by semblance analysis, a coherence-based scan over trial slowness values. Dividing the DAS strain-rate by $s$ yields ground acceleration comparable to a seismometer. On this sits the multichannel Wiener filter, whose residual is $R(\omega)=1-\vec{C}_{DS}^{\dagger}(\omega)\vec{C}_{DD}^{-1}(\omega)\vec{C}_{DS}(\omega)/C_{SS}(\omega)$; the square root $\sqrt{R}$ is the reported noise-reduction factor. The gauge length $L_g$ acts as a spatial averaging scale, and the dense channel spacing lets the array reconstruct wavelengths longer than a single gauge length.
What would settle it
Train the multichannel Wiener filter on the first half of the one-hour dataset and apply it to the second half, then recompute the residual factor at 20 Hz for the vertical seismometer component using six DAS channels; a held-out residual substantially above 0.11 (for example, above 0.3) would show the headline cancellation factor does not generalize.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that DAS is not just a detector of seismic events but a quantitatively equivalent seismic instrument. After converting DAS strain-rate to ground acceleration through the apparent slowness, the DAS waveforms, power spectra, and Bland–Altman agreement match a colocated broadband seismometer. Used as witness channels in a multichannel Wiener filter, six DAS sensors cancel a seismometer's vertical noise to a residual factor of 0.11 at 20 Hz, exactly matching six geophone channels, and outperform the geophones on the horizontal component (0.04 versus 0.15). The paper also reports that this high colocated correlation contradicts a prior theoretical expectation that strain meters must be placed far from the test mass to be useful for Newtonian-noise cancellation.
Load-bearing premise
The cancellation filter is tuned and evaluated on the same one-hour dataset, with no separate test period mentioned, so the reported residual factors may be optimistically biased by in-sample fitting.
Editorial extensions
If this is right
- Existing fiber infrastructure can be turned into dense seismic arrays, reducing the deployment cost of witness sensors for Newtonian-noise cancellation at gravitational wave observatories.
- DAS matches a geophone array for vertical-component cancellation (0.11 residual at 20 Hz with six channels) and exceeds it on the horizontal component (0.04 versus 0.15).
- DAS can reconstruct geophone signals with an average correlation of 0.97 and coherence above 0.7 above 3 Hz, allowing it to substitute for or supplement sparse geophone networks.
- DAS coherence lengths of about 11 m at 3–5 Hz and 23 m at 6–10 Hz resolve ground-motion structure finer than typical 25 m geophone spacing, which is relevant near detector test masses.
- Combining DAS with traditional seismometers and geophones is identified as a promising direction for further improving Newtonian-noise cancellation and environmental monitoring.
Reading between the lines
- A held-out evaluation is the natural next step: the paper's residual factors are in-sample, so a train/test split would give a realistic bound on DAS cancellation performance.
- Because DAS measures strain only along the fiber direction, the vertical-component improvement with more channels suggests that helical or multi-orientation cable layouts could push vertical cancellation below the quoted 0.11; that is a direct, testable extension.
- The reported coherence lengths of roughly 11 m and 23 m in the 3–10 Hz band imply an optimal witness-channel spacing for cancellation; arrays much denser than that may add little information while raising processing cost.
- If the colocated high correlation survives at quieter detector sites, existing telecommunication fiber around gravitational wave observatories could be repurposed for Newtonian-noise monitoring, with the open engineering question being real-time slowness estimation and adaptive filter updates.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper argues that distributed acoustic sensing (DAS) is a viable and scalable alternative to seismometers and geophones for environmental monitoring and Newtonian-noise mitigation in gravitational-wave detectors. Using field data from the DESY campus (vibrotruck, earthquake, thunderstorm, HVAC), the authors convert DAS strain-rate measurements to ground acceleration using a semblance-estimated slowness, compare waveforms and power spectra with co-located broadband seismometers, use multichannel Wiener filtering to cancel the vertical and east components of a seismometer with DAS and with geophones, and report residual noise factors down to 0.11 at 20 Hz. They also report a 0.97 correlation for Wiener-filter prediction of geophone signals from DAS, DAS coherence lengths of roughly 11 m and 23 m in the 3-10 Hz band, and several environmental monitoring examples.
Significance. If the central cancellation claim is established out of sample, the paper would provide a practically important result: DAS could serve as a scalable witness sensor for Newtonian-noise cancellation in current and next-generation gravitational-wave detectors, using existing fiber infrastructure and dense spatial sampling. The paper also contributes useful field demonstrations of DAS-to-seismometer conversion, coherence-length measurements, and environmental monitoring. However, the headline cancellation performance is currently supported only by in-sample Wiener-filter residuals, and the paper contains an internal inconsistency in the reported DAS vertical residual (0.11 vs 0.14). Those issues directly affect the central claim, so the significance of the result is not yet established.
major comments (4)
- [Section 4.3, Eq. (13); Section 5.2] The residual factors reported in Table 1 and the abstract are computed from cross-spectral matrices CDS, CDD, and CSS estimated with Daniell's method on the same 1-hour dataset used for evaluation; no train/test split or cross-validation is described. Since the Wiener filter coefficients are fitted to the same noise realization that is being cancelled, the in-sample residual sqrt(R(ω)) is optimistically biased, and the headline 'residual noise factor of 0.11 at 20 Hz' is not established as an out-of-sample cancellation factor. Please report out-of-sample residuals (or k-fold cross-validation) with uncertainties, and specify the number of frequency bins/taps used in the Wiener filter.
- [Abstract; Table 1; Section 5.2.1] The abstract and Table 1 report a DAS vertical-component residual of 0.11 at 20 Hz with six channels, while Section 5.2.1 states that 'DAS achieved 0.14' for the same configuration; Section 5.2 also gives 0.14 for six DAS channels at 20 Hz. This internal inconsistency affects the central comparison and must be resolved before the headline claim can be accepted.
- [Section 4.2; Section 5.3] The geophone-prediction correlation of 0.97 and the coherence 'above 0.7' are computed with Wiener filter coefficients estimated from the same data used for the prediction, so this metric is a fit-quality measure rather than a predictive validation. Please provide an out-of-sample evaluation and specify the Wiener filter length or number of taps, the frequency smoothing, and the dataset split.
- [Section 4.1; Section 3.1.1; Section 6] The DAS-to-acceleration conversion relies on an apparent slowness estimated from the DAS data itself via semblance (Eq. 8), and the selected DAS channels used for cancellation are not specified in Section 4.3 or Table 1. Because both the slowness and the channel selection are adapted to the record being analyzed, they can inflate the reported correlation and cancellation performance; please describe the channel-selection criteria, report sensitivity to slowness and channel choice, and ideally validate the conversion with an independent slowness estimate from the geophone array.
minor comments (5)
- [Section 6.0.1] The Gaussian correlation length Lc is a fitted parameter, and the reported coherence lengths (11 m and 23 m) come from a single vibrotruck event; please report the fit uncertainty and state whether the result is stable across events and frequency bands.
- [Equation (11)] The wavenumber vector is written as k = 2π s(f)·(sinθ, cosθ, 0), which appears dimensionally inconsistent with the slowness definition in Eq. (6); please define f and s explicitly and check the units.
- [References] Reference [20] is cited for LPSD, but the citation appears to point to a software package rather than the original LPSD algorithm; please cite the original algorithm and provide a version/commit for the software if used.
- [Throughout] There are several typographical and nomenclature issues, including 'Fracensca' in the author list, 'Deutches' in the acknowledgments, 'compliment' for 'complement' in Section 6, and the irregular spacing in 'W A VE'; these should be corrected.
- [Figure 13] The PDF plot of residuals should specify how the residuals are aggregated over frequency and channels; the statement that the maximum lies between 0.00 and 0.02 needs error bars or sample counts to be meaningful.
Circularity Check
In-sample Wiener filtering turns the headline residual 0.11 and geophone-prediction correlation 0.97 into fit statistics; no train/test split is described.
-
fitted input called prediction
[Section 4.2 (Wiener filter training), Section 5.3 and Figure 14 (reported correlation)]
"In this study, we used data from a 10-channel DAS array, down-sampled to 200 Hz, over a period of 1 hour. ... We assumed stationarity within the data set needed for Wiener filtering. ... The Wiener filter successfully estimated the geophone data from the DAS measurements, demonstrating the potential of DAS systems for seismic signal reconstruction. The correlation coefficient between the estimated and true geophone signals averaged 0.97, indicating a strong linear relationship between the two."
The Wiener filter coefficients h(t) in Eq. (12) are estimated from the same 1-hour DAS/geophone segment on which the 0.97 correlation is later reported. No train/test split or cross-validation is described. Since the optimal Wiener filter minimizes the mean-square error on the training segment by construction, the 'predicted' geophone trace is an in-sample reconstruction. The high correlation is therefore a fit statistic, not an out-of-sample prediction, despite being presented as evidence that DAS can 'predict' or 'estimate' geophone signals.
-
fitted input called prediction
[Section 4.3, Eq. (13); Section 5.2, 5.2.1; Table 1; abstract]
"The noise cancellation technique is based on the multichannel Wiener filter and utilizes cross-spectral matrices, expressed as [19]: ... . Residual p R(ω) represents the noise reduction factor that DAS is able to achieve. ... Cross-spectral and auto-spectral densities were estimated using Daniell's method [20]."
Eq. (13) defines the residual as the theoretical output of the optimal Wiener filter. The cross-spectral matrices CDD, CDS, and CSS are estimated on the same 1-hour dataset used to report the residuals in Section 5.2 and Table 1; no train/test split is mentioned. With up to six DAS channels and many frequency bins, the in-sample optimal residual is optimistically biased. Calling √R 'the noise reduction factor that DAS is able to achieve' and quoting 0.11 at 20 Hz presents this in-sample fit as achieved cancellation. The geophone comparison uses the same in-sample protocol, so the symmetry does not validate the absolute factor.
full rationale
The paper's two central quantitative claims—geophone prediction and noise-cancellation residuals—are computed with Wiener filters whose coefficients or cross-spectral matrices are estimated from the same 1-hour dataset on which the results are evaluated. Section 4.2 describes a 1-hour DAS/geophone dataset and assumes stationarity; Section 4.3 estimates cross-spectral densities with Daniell's method and then Eq. (13) reports the optimal residual. No train/test split, cross-validation, or independent validation segment is described. The 0.97 correlation and the 0.11 residual at 20 Hz are therefore in-sample fit statistics, making the headline 'prediction' partially reduce to a fit. There is independent non-circular content: the DAS-to-seismometer conversion uses slowness estimated from DAS moveout, not from the seismometer, and the waveform/PSD/Bland-Altman comparisons are external agreement checks. Self-citations [3] and [19] are not load-bearing because Eq. (13) is a standard multichannel Wiener formula and the geophone comparison is recomputed in this paper. An unrelated data-integrity concern is that the abstract and Table 1 give a DAS vertical residual of 0.11 at 20 Hz with six channels, while Sections 5.2 and 5.2.1 state 0.14; this inconsistency does not affect the circularity score but should be resolved. Overall, the central cancellation factor and geophone-prediction correlation are partly forced by the in-sample fitting procedure, so the score is 6.
Assumptions & free parameters
free parameters (4)
- Apparent slowness s =
estimated per time-frequency via semblance
- Wiener filter coefficients h(t) =
frequency-domain optimal filter
- Number and selection of DAS channels =
2-6 channels
- Gaussian correlation length Lc in coherence fits =
11.76 m and 23.55 m
assumptions (4)
- domain assumption Ground motion can be modeled as plane waves with a single apparent slowness s at each time/frequency.
- domain assumption DAS strain rate is related to seismometer velocity by a finite-difference approximation over the gauge length.
- domain assumption The noise field and signal are stationary over the 1-hour analysis window, so Wiener filtering is valid.
- domain assumption The vertical seismometer channel is a valid proxy for test-mass displacement noise in Newtonian noise cancellation.
Cite this review
Pith. "Pith review of Distributed Acoustic Sensing for Environmental Monitoring, and Newtonian Noise Mitigation:Comparable Sensitivity to Seismometers." pith.science (2026). https://pith.science/paper/U5BYWUAC
@misc{pith2026250713523,
author = {Pith},
title = {Pith review of: Distributed Acoustic Sensing for Environmental Monitoring, and Newtonian Noise Mitigation:Comparable Sensitivity to Seismometers},
year = {2026},
howpublished = {\url{https://pith.science/paper/U5BYWUAC}},
note = {Machine review of arXiv:2507.13523}
}
read the original abstract
Newtonian noise limits the low-frequency sensitivity of ground-based gravitational wave detectors. While seismometers and geophones are commonly employed to monitor ground motion for Newtonian noise cancellation, their limited spatial coverage and high deployment costs hinder scalability. In this study, we demonstrate that distributed acoustic sensing offers a viable and scalable alternative, providing performance comparable to that of conventional seismic instruments. Using data from acoustic sensing and colocated seismometers during both natural and controlled events, we observe a strong correlation, exceeding 0.8, between the two sensor types in the 3 to 20 Hz frequency band relevant for Newtonian noise. Moreover, when distributed acoustic sensing data are used to predict geophone signals, the correlation remains high, above 0.7, indicating that distributed acoustic sensing accurately captures both the spatial and spectral features of ground motion. As a case study, we apply distributed acoustic sensing data to cancel noise recorded by the vertical component of a seismometer and compare the results with those obtained using geophone data for the same task. Both distributed acoustic sensing and geophone-based cancellations yield a residual noise factor of 0.11 at 20 Hz. These findings confirm the feasibility of using distributed acoustic sensing for Newtonian noise mitigation and highlight its potential, in combination with traditional seismic sensors, to improve environmental monitoring and noise suppression in current and next-generation gravitational wave observatories.
Figures
Figures from the paper (15 more)
Forward citations
Cited by 1 Pith paper
-
LIGO A$^\sharp$: Detector Design and Science Prospects Beyond A+
LIGO A#, a room-temperature upgrade of the LIGO detectors proposed for the 2030s, would broaden sensitivity and increase projected compact-binary detection rates by factors of four to eight relative to A+.
Reference graph
Works this paper leans on
-
[1]
An Introduction to Distributed Optical Fibre Sensors
A. H. Hartog, “An Introduction to Distributed Optical Fibre Sensors”. CRC Press, 2017
work page 2017
-
[2]
Terrestrial gravity fluctuations
J. Harms, “Terrestrial gravity fluctuations”, Living Reviews in Relativity, 2019
work page 2019
-
[3]
Coherent and Incoherent Noise Cancellation using Distributed Optical Fiber Sensors
R. Rading, K.S. Isleif, “Coherent and Incoherent Noise Cancellation using Distributed Optical Fiber Sensors”, IEEE Photonic Society Summer Topicals, 2024
work page 2024
-
[4]
Y. Li, M. Karrenbach, and J. B. Ajo-Franklin, “A Literature Review: Distributed Acoustic Sensing (DAS) Geophysical Applications Over the Past 20 Years”, pg. 229–291, 2021
work page 2021
-
[5]
Seismic gravity-gradient noise in interferometric gravitational-wave detectors
S. Hughes and K. Thorne, “Seismic gravity-gradient noise in interferometric gravitational-wave detectors”, Physical Review D 58.12, 1998
work page 1998
-
[6]
Propagation in optical fibres,
A.H. Hartog, “Propagation in optical fibres,” in An Introduction to Distributed Optical Fibre Sensors, CRC Press, 2017
work page 2017
-
[7]
Photon probe — an optical-fiber time-domain reflectometer,
S.D. Personick, “Photon probe — an optical-fiber time-domain reflectometer,” The Bell System Technical Journal, vol. 56, no. 3, pp. 355–366, 1977. 24
work page 1977
-
[8]
DAS parameter definitions and tests
SEAFOM. “DAS parameter definitions and tests.”, [Online]. Available: http://seafom.com
Show all 20 references
-
[9]
Ground motion response to an ML 4.3 earthquake using co-located distributed acoustic sensing and seismometer arrays
H.F. Wang, X. Zeng, D.E. Miller, D. Fratta, K.L. Feigl, C.H. Thurber, and R.J. Mel- lors,“Ground motion response to an ML 4.3 earthquake using co-located distributed acoustic sensing and seismometer arrays”, Geophys. J. Int., 213, 2020–2036, 2018
2020
-
[10]
Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle
T.M. Daley, D.E. Miller, K. Dodds, P. Cook, and B.M. Freifeld, “Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle”, Alabama, Geophys. Prospect., 64(5), 1318–1334, 2015
2015
-
[11]
On the Broadband Instrument Response of Fiber-Optic DAS
N.J. Lindsey, H. Rademacher, and J.B. Ajo-Franklin, “On the Broadband Instrument Response of Fiber-Optic DAS” Arrays, J. Geophys. Res.-Sol. Ea., 125, 2020
2020
-
[12]
A Guide to the Seismic Wavefield as seen by DAS
L.N.B. Kennett. “A Guide to the Seismic Wavefield as seen by DAS”, Research School of Earth Sciences, Australian National University, 2024
2024
-
[13]
Complex seismic trace analysis
M.T. Taner, F. Koehler, and R.E. Sheriff, “Complex seismic trace analysis”, Geophysics, 44, 1041–1063, 1979
1979
-
[14]
Semblance and other Coherency Measures for Multichannel Data
N.S. Neidell and M.T. Taner, “Semblance and other Coherency Measures for Multichannel Data”, Geophysics, 36, 482–497, 1971
1971
-
[15]
Complex Semblance and Its Application
T. Shi and S. Huo,“Complex Semblance and Its Application”, J. Earth Sci., 30, 849–852, 2019
2019
-
[16]
Multiple emitter location and signal parameter estimation
R. Schmidt, “Multiple emitter location and signal parameter estimation”, IEEE Trans. An- tennas Propag., 34(3), 276–280, 1986
1986
-
[17]
Analysis of seismic waves crossing the Santa Clara valley using the three component MUSIQUE array algorithm
M. Hobiger, C. Cornou, P.-Y. Bard, N. Le Bihan, and W. Imperatori, “Analysis of seismic waves crossing the Santa Clara valley using the three component MUSIQUE array algorithm”, Geophys. J. Int., 207(1), 439–456, 2016
2016
-
[18]
W A VE Initiative: https://wave-hamburg.eu/
-
[19]
Newtonian Noise Studies in 2nd and 3rd Generation Gravitational Wave In- terferometric Detectors
F. Badaracco, “Newtonian Noise Studies in 2nd and 3rd Generation Gravitational Wave In- terferometric Detectors”, PhD thesis, Gran Sasso Science Institute, 2020
2020
- [20]
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.