REVIEW 3 major objections 4 minor 28 references
Atmospheric muon flux suppression at potential new low-radiation underground physics laboratory in Israel
T0 review · 3 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read A tunnel under Kokhav HaYarden suppresses atmospheric muons by a factor of 4,456, with an effective overburden of about 873 m.w.e.
desk verdict A useful first muon-flux measurement at a candidate Israeli underground site, with a plausible suppression factor, but the claimed precision rests on an unquantified above-ground dead-time correction. 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 central object is a portable muon hodoscope: four vertically stacked plastic scintillator plates ($144 \times 60 \times 1.2$ cm$^3$) with photomultiplier tubes at both ends of each plate, triggered by a coincidence on the top two plates and read out by two oscilloscopes. Track angles are reconstructed from the time delay between the two PMTs on each plate, giving a one-dimensional position along the plate and hence a zenith angle. The argument is carried by the ratio of trigger rates at the two sites: because the same detector, threshold, and geometry are used above ground and underground, the efficiency and acceptance factors cancel, so the rate ratio is the flux suppression factor. The measured suppression is then converted to an equivalent vertical depth using the differential muon intensity function $I_\mu(h_0) = 68 \times 10^{-6} e^{-h_0/285} + 2 \times 10^{-6} e^{-h_0/698}$ (with $h_0$ in m.w.e.), and cross-checked against a simulation that propagates a Gaisser-parametrized muon spectrum through a two-dimensional slant-depth map built from elevation data.
What would settle it
Run the same hodoscope above ground while lowering the photomultiplier high voltage in steps and record the trigger rate; if the rate does not scale linearly with the muon rate, then the efficiency is rate-dependent and the sea-level rate used in the suppression factor is biased. A second check is to compare the measured angular muon distribution at Kokhav HaYarden with a full muon transport simulation using the actual rock density; a systematic angular offset would implicate the slant-depth model rather than the flux measurement.
Extended reading notes
Core claim
The central claim is that the integrated atmospheric muon flux at the Kokhav HaYarden site is suppressed by a factor of $4456 \pm 77$ relative to sea level, corresponding to $3.75 \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ and an equivalent vertical depth of about 873 m.w.e. The authors establish this by measuring the trigger rate of a hodoscope with identical geometry at the above-ground Tel Aviv University site and at the underground site, in two orientations, and taking the ratio so that detector efficiency and acceptance cancel. They further show that the angular distribution of reconstructed muon tracks is biased away from the thickest rock overburden, revealing the muon shadow of the mountain, and that the local radon concentration is low. A simplified simulation that propagates a standard parametrization of the sea-level muon spectrum through the slant-depth map predicts about 860 m.w.e., in good agreement with the measured value. The conclusion is that the tunnel branch offers working conditions suitable for a new low-radiation underground laboratory, with the option of 25% more overburden deeper in the tunnel.
Load-bearing premise
The suppression factor assumes that the above-ground trigger rate of about $3.8 \times 10^1$ s$^{-1}$ at Tel Aviv University is not distorted by dead time or by two muons arriving close together, so that the detector's efficiency and acceptance cancel exactly between the above-ground and underground runs; the paper quotes only statistical errors and does not quantify rate-dependent effects.
Editorial extensions
If this is right
- At roughly 873 m.w.e., the Kokhav HaYarden site sits between shallow and deep underground laboratories, making it a candidate for material screening, high-voltage component testing, low-mass dark matter searches, and some neutrino experiments.
- The deeper, currently inaccessible part of the tunnel would add about 25% more vertical rock, which would further reduce the muon flux below the already measured suppression.
- The low radon activity of $28.3 \pm 14.0$ Bq m$^{-3}$ means radon is not a limiting background for the planned uses, although the high temperature and humidity will require attention.
- The measured angular muon shadow confirms that the slant-depth map built from elevation data captures the real overburden, validating the simulation-based estimate of about 860 m.w.e.
- The same portable hodoscope and ratio method can be carried to the deeper Manara pumped-storage site, which could reach roughly 2000 m.w.e. and enter the depth class of established deep laboratories.
Reading between the lines
- Inference: if the rate-ratio cancellation is robust, the same hodoscope could be redeployed at Manara with a longer sea-level control run to test for rate-dependent losses, since a small dead-time correction at the high above-ground rate would shift the quoted suppression factor.
- Inference: the agreement between the measured 873 m.w.e. and simulated 860 m.w.e. suggests the elevation-derived slant-depth map and a uniform rock density of 2.6 g cm$^{-3}$ are adequate for survey-level planning, but a direct density measurement would tighten the depth estimate.
- Inference: for superconducting qubit experiments, a factor-of-4,456 reduction in ionizing radiation could meaningfully lower radiation-induced quasiparticle poisoning, provided the site's humidity and temperature are controlled.
- Inference: the measured flux could be compared against depth-flux curves from other intermediate-depth labs; a systematic offset would point to local rock composition or the sea-level normalization rather than the detector.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a first exploratory muon-flux measurement at the proposed Kokhav HaYarden (KY) underground site in Israel, using a four-plate plastic scintillator hodoscope with a two-fold (top-plate) coincidence trigger. Data were taken above ground at Tel Aviv University (TAU) and underground at KY in two orientations, South-North and West-East. The central result is a suppression factor of 4456 ± 77 relative to TAU, corresponding to an integrated muon flux of 3.75 ± 0.06 × 10^-6 cm^-2 s^-1 and an effective vertical overburden of roughly 873 m.w.e. The paper also presents angular distributions that show the mountain's muon shadow, a simulation-based slant-depth map, and ancillary radon, temperature, and humidity measurements. The authors state that the quoted uncertainty is statistical only and that constant-efficiency systematics cancel in the ratio.
Significance. The direct ratio measurement is a simple and appropriate way to obtain a site-specific suppression factor without relying on detailed geological simulation, and the cancellation of detector-efficiency and acceptance terms that are constant between runs is a genuine strength. The angular distributions and the full-sky slant-depth map make the topographic effect directly visible, and the radon measurement is a useful auxiliary datum for site evaluation. If the rate-dependent systematic discussed below is quantified, the result would be a useful addition to the sparse literature on shallow and intermediate underground sites. At present, however, the headline precision is not demonstrated because the high-rate TAU run is not checked for dead time, pileup, or accidental coincidences.
major comments (3)
- [§4, Table 2] The central suppression factor is defined as the ratio of the TAU rate (3.80 ± 0.04 × 10^1 s^-1) to the average KY rate (8.535 × 10^-3 s^-1), but no dead-time, pileup, or accidental-coincidence correction is reported for the TAU run. Rate-dependent losses do not cancel in the TAU/KY ratio because the trigger rates differ by more than three orders of magnitude; a 2% trigger loss at TAU changes the suppression factor by about 90, which is outside the quoted ±77 statistical uncertainty. The text explicitly reports "only the statistical error here," so the claim of 4456 ± 77 is not yet demonstrated at the stated precision. Please use the stored TAU triggers or waveforms to measure dead time and check rate linearity, or add a quantitative conservative systematic uncertainty.
- [§4, correlated-muon sentence] The sentence claiming that correlated muons "can only increase the inferred shielding" is unclear and appears to contradict the preceding multiplicity statement. If two or more muons arriving together are counted as a single trigger, the measured TAU rate is an underestimate of the true muon arrival rate, which would reduce the measured suppression factor rather than increase it. Please clarify the logic and, if multi-muon events are intended to be a correction, provide a quantitative estimate of their fraction at TAU.
- [§4, Eq. (1)] The effective overburden of 873 m.w.e. is obtained by inverting Eq. (1), but the cited Mei-Hime relation is stated to be appropriate for 1–10 km.w.e., and 873 m.w.e. lies below that range. Extrapolating the formula outside its stated validity could bias the inferred depth; please justify the extrapolation or use a depth-intensity relation validated at shallow depths. Also clarify whether the quoted 873 m.w.e. depends on the rock-density assumption (2.6 g cm^-3) used for the m.w.e. conversion in Fig. 2.
minor comments (4)
- [Abstract and §4] The flux value is written as "3.75 ± 0.06 × 10^-6 cm^-2 s^-1" in both the abstract and Section 4; please parenthesize the uncertainty as (3.75 ± 0.06) × 10^-6 cm^-2 s^-1 to avoid ambiguity.
- [§4 and Fig. 10] The text says the muon shadow of the mountain is visible "in the western direction," while the Fig. 10 caption describes an eastward bias in the W-E orientation; please make this consistent.
- [§4, Fig. 11] Please state explicitly that the "faulty" events excluded from the angular distributions are nevertheless included in the trigger-rate counts used for the suppression factor, since the rate calculation is based on all triggers.
- [Throughout] There are several typographical and formatting issues, including "Tel A viv" in the author affiliations, "T able" in table captions, "auxillary" in Table 3, and "Galille" in Section 1; these should be corrected.
Circularity Check
No circularity: the suppression factor is a directly measured rate ratio, the m.w.e. conversion uses an external empirical function, and self-citations are not load-bearing.
full rationale
The paper's central claim (4456 +/- 77 suppression, 3.75e-6 cm^-2 s^-1) is a measured ratio of trigger rates at TAU and KY (Table 2), not the output of a fitted model. The conversion to an equivalent vertical depth uses the external Mei-Hime intensity function (Eq. 1, ref [19]), an independent empirical benchmark whose parameters are not fit to this dataset; agreement with the simplified simulation (approximately 860 m.w.e.) is a consistency check, not an input. The sea-level flux normalization (1.67e-2 cm^-2 s^-1) likewise comes from PDG [18]. No load-bearing argument reduces to a self-citation: references [16] and [17] describe the scintillator hardware, and [24] is an application example, not evidence for the flux result. The paper explicitly reports only statistical errors and leaves rate-dependent systematics (dead time, pileup) unquantified; that is a correctness or robustness limitation, not circularity, because the quantity quoted remains a direct measurement rather than a quantity constructed from its own assumptions. No fitted parameter is renamed as a prediction, and no uniqueness theorem or ansatz is imported from the authors' prior work.
Assumptions & free parameters
free parameters (1)
- Rock density for m.w.e. conversion =
2.6 g cm^-3
assumptions (5)
- domain assumption Sea-level muon flux of 1.67 x 10^-2 cm^-2 s^-1 from PDG is accurate for the TAU site, neglecting a small 44 m elevation deviation.
- domain assumption Detector efficiency, acceptance, and trigger behavior are identical at TAU and KY and therefore cancel in the rate ratio.
- domain assumption The underground trigger rate contains no significant accidental coincidence or radioactive background, so no background subtraction is needed.
- domain assumption The Mei-Hime intensity function in Eq. (1) is valid at the derived depth of roughly 873 m.w.e., even though the paper notes its stated range is 1 to 10 km.w.e.
- domain assumption The simulation using the Gaisser muon spectrum and continuous energy loss model adequately describes the site's topography and effective overburden.
Cite this review
Pith. "Pith review of Atmospheric muon flux suppression at potential new low-radiation underground physics laboratory in Israel." pith.science (2026). https://pith.science/paper/5IQNOCYS
@misc{pith2026250415102,
author = {Pith},
title = {Pith review of: Atmospheric muon flux suppression at potential new low-radiation underground physics laboratory in Israel},
year = {2026},
howpublished = {\url{https://pith.science/paper/5IQNOCYS}},
note = {Machine review of arXiv:2504.15102}
}
abstract
The residual atmospheric muon flux was measured at a candidate site for a new underground, low-radiation physics laboratory beneath the Kokhav HaYarden national park in Israel. Located inside the tunnels of a hydroelectric pumped-storage facility operating since 2024, the proposed site benefits from a vertical rock overburden of 361 meters, large potential floorspace, and easy access by road. A muon hodoscope of vertically stacked wide-area $144 \times 60 \times 1.2 $ cm$^3$ plastic scintillator plates was employed to measure the suppression in the integrated muon flux at the site as compared with above ground at sea level. Data-taking took place in mid-August of 2024 for several days and was split into South-North and West-East orientations to account for the directional acceptances due to the geometry of the detector. The suppression factor is reported at $4456 \pm 77$, expressed as $3.75 \pm 0.06 \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ in absolute terms, corresponding to an effective overburden of roughly $873$ m.w.e.. A deeper location at the site may also be available, but it could not be reached at this time. Furthermore, the asymmetric topography of the mountain above and its muon shadow are clearly visible in the angular data. Finally, auxiliary environmental measurements recorded low background radon activity at $28.3 \pm 14.0$ Bq m$^{-3}$. The experimental campaign thus succeeded in demonstrating the viability of the site's working conditions for future scientific research.
Reference graph
Works this paper leans on
-
[1]
Journal of Physics G: Nuclear and Particle Physics 43(1), 013001 (2015) https://doi.org/10
Undagoitia, T.M., Rauch, L.: Dark mat- ter direct-detection experiments. Journal of Physics G: Nuclear and Particle Physics 43(1), 013001 (2015) https://doi.org/10. 1088/0954-3899/43/1/013001
work page 2015
-
[2]
Cooley, J.: Dark Matter direct detection of classical WIMPs. SciPost Phys. Lect. Notes, 55 (2022) https://doi.org/10.21468/ SciPostPhysLectNotes.55
work page 2022
-
[3]
Annual Review of Nuclear and Particle Science69(Volume 69, 2019), 219–251 (2019) https://doi.org/10
Dolinski, M.J., Poon, A.W.P., Rodejo- hann, W.: Neutrinoless double-beta decay: Status and prospects. Annual Review of Nuclear and Particle Science69(Volume 69, 2019), 219–251 (2019) https://doi.org/10. 1146/annurev-nucl-101918-023407
2019
-
[4]
Kumaran, S., Ludhova, L., Penek, Ö., Set- tanta, G.: Borexino results on neutrinos from thesunandearth.Universe 7(7)(2021)https: //doi.org/10.3390/universe7070231
-
[5]
Joo, K.K.: New results from reno amp; prospects with reno-50. Journal of Physics: Conference Series 888(1), 012012 (2017) https://doi.org/10.1088/1742-6596/888/1/ 012012
-
[6]
: Reducing the impact of radioactivity on quantum circuits in a deep-underground facility
Cardani, L., et al. : Reducing the impact of radioactivity on quantum circuits in a deep-underground facility. Nature Commun. 12(1), 2733 (2021) https://doi.org/10.1038/ s41467-021-23032-z arXiv:2005.02286 [cond- mat.supr-con]
arXiv 2021
-
[7]
: Impact of ionizing radiation on superconducting qubit coher- ence
Vepsäläinen, A., et al. : Impact of ionizing radiation on superconducting qubit coher- ence. Nature 584(7822), 551–556 (2020) https://doi.org/10.1038/s41586-020-2619-8 arXiv:2001.09190 [quant-ph]
arXiv 2020
-
[8]
New Journal of Physics27(2), 023014 (2025) https://doi.org/10.1088/1367-2630/adaedc
Bertoldo, E., Pérez Sánchez, V., Martínez, M., Martínez, M., Khalife, H., Forn-Díaz, P.: Cosmic muon flux attenuation meth- ods for superconducting qubit experiments. New Journal of Physics27(2), 023014 (2025) https://doi.org/10.1088/1367-2630/adaedc
Show all 28 references
-
[9]
: First limit on wimp cross section with low background csi(tl) crystal detector
Lee, H.S., et al. : First limit on wimp cross section with low background csi(tl) crystal detector. Phys. Lett. B633, 201–208 (2006) https://doi.org/10.1016/j.physletb.2005.12. 035 arXiv:astro-ph/0509080
2006 arXiv
-
[10]
Astroparticle Physics 126, 102533 (2021) https://doi.org/10.1016/ j.astropartphys.2020.102533
Yoon, Y.S., Kim, J., Park, H.: Neutron background measurement for rare event search experiments in the yangyang under- ground laboratory. Astroparticle Physics 126, 102533 (2021) https://doi.org/10.1016/ j.astropartphys.2020.102533
2021
-
[11]
https://www.hutchisonwater.com/ projects/star-pumped-storage/
Hutchison Water: Star Pumped Storage Project. https://www.hutchisonwater.com/ projects/star-pumped-storage/. Accessed: 2025-04-17
2025
-
[12]
https://ellomay.com/ under-construction/
Ellomay Capital Ltd.: Under Construc- tion Projects. https://ellomay.com/ under-construction/. Accessed: 2025-04-17
2025
-
[13]
https://earth
Google Earth: Satellite imagery of Kokhav HaYarden pumped storage facility, Israel (32.5969° N, 35.5292 ° E). https://earth. google.com/
-
[14]
https://developers.google.com/ maps/documentation/elevation (2025)
Google Maps Platform: Google Maps Ele- vation API. https://developers.google.com/ maps/documentation/elevation (2025)
2025
-
[15]
: Review of particle physics
Workman, R.L., et al. : Review of particle physics. Prog. Theor. Exp. Phys.2022, 083– 01 (2022). (Particle Data Group)
2022
-
[16]
In: 2003 IEEE Nuclear Science Symposium
Etzion, E., al.: The cosmic ray hodoscopes for testing thin gap chambers at the tech- nion and tel aviv university. In: 2003 IEEE Nuclear Science Symposium. Confer- ence Record (IEEE Cat. No.03CH37515), vol. 5, pp. 3744–37495 (2003). https://doi.org/10. 1109/NSSMIC.2003.1352720 7
2003 arXiv
-
[17]
In: IEEE Symposium Conference Record Nuclear Science 2004., vol
Etzion, E., al.: The certification of atlas thin gap chambers produced in israel and china. In: IEEE Symposium Conference Record Nuclear Science 2004., vol. 1, pp. 236–2411 (2004). https://doi.org/10.1109/ NSSMIC.2004.1462188
2004 arXiv
-
[18]
Tanabashi, M.: Review of particle physics. Phys. Rev. D98, 030001 (2018) https://doi. org/10.1103/PhysRevD.98.030001
2018 doi
-
[19]
Mei, D.-M., Hime, A.: Muon-induced back- ground study for underground laboratories. Phys. Rev. D73, 053004 (2006) https://doi. org/10.1103/PhysRevD.73.053004
2006 doi
-
[20]
https://www.airthings
Airthings:AirthingsforBusiness–IndoorAir Quality Monitors. https://www.airthings. com/en/pro. Accessed: 2025-04-17
2025
-
[21]
: Gator: a low-background counting facility at the Gran Sasso Under- ground Laboratory
Baudis, L., et al. : Gator: a low-background counting facility at the Gran Sasso Under- ground Laboratory. JINST 6, 08010 (2011) https://doi.org/10.1088/1748-0221/6/08/ P08010 arXiv:1103.2125 [astro-ph.IM]
2011 arXiv
-
[22]
JINST 11(12), 12017 (2016) https: //doi.org/10.1088/1748-0221/11/12/P12017 arXiv:1606.03983 [physics.ins-det]
Sivers, M., Hofmann, B.A., Rosén, r., Schumann, M.: The GeMSE facility for low-background γ-ray spectrome- try. JINST 11(12), 12017 (2016) https: //doi.org/10.1088/1748-0221/11/12/P12017 arXiv:1606.03983 [physics.ins-det]
2016 arXiv
-
[23]
Heusser, G., Weber, M., Hakenmüller, J., Laubenstein, M., Lindner, M., Maneschg, W., Simgen, H., Stolzen- burg, D., Strecker, H.: GIOVE - A new detector setup for high sensitivity ger- manium spectroscopy at shallow depth. Eur. Phys. J. C 75(11), 531 (2015) https: //doi.org/10...
2015 arXiv
-
[24]
Budnik,R.,Chesnovsky,O.,Slone,O.,Volan- sky, T.: Direct Detection of Light Dark Mat- ter and Solar Neutrinos via Color Center Production in Crystals. Phys. Lett. B 782, 242–250 (2018) https://doi.org/10.1016/j. physletb.2018.04.063 arXiv:1705.03016 [hep- ph]
2018 arXiv
-
[25]
: First Direct-Detection Results on Sub-GeV Dark Matter Using the SENSEI Detector at SNOLAB
Adari, P., et al. : First Direct-Detection Results on Sub-GeV Dark Matter Using the SENSEI Detector at SNOLAB. Phys. Rev. Lett. 134(1), 011804 (2025) https: //doi.org/10.1103/PhysRevLett.134.011804 arXiv:2312.13342 [astro-ph.CO]
2025 arXiv
-
[26]
: First Constraints from DAMIC-M on Sub-GeV Dark-Matter Par- ticles Interacting with Electrons
Arnquist, I., et al. : First Constraints from DAMIC-M on Sub-GeV Dark-Matter Par- ticles Interacting with Electrons. Phys. Rev. Lett. 130(17), 171003 (2023) https: //doi.org/10.1103/PhysRevLett.130.171003 arXiv:2302.02372 [hep-ex]
2023 arXiv
-
[27]
Esch, E.-I., Bowles, T.J., Hime, A., Pichlmaier, A., Reifarth, R., Wollnik, H.: The Cosmic ray muon flux at WIPP. Nucl. Instrum. Meth. A 538, 516–525 (2005) https://doi.org/10.1016/j.nima.2004.09.005 arXiv:astro-ph/0408486
2005 arXiv
-
[28]
: The EXO-200 detector, part I: Detector design and construc- tion
Auger, M., et al. : The EXO-200 detector, part I: Detector design and construc- tion. JINST 7, 05010 (2012) https: //doi.org/10.1088/1748-0221/7/05/P05010 arXiv:1202.2192 [physics.ins-det] 8
2012 arXiv
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.