pith. sign in

arxiv: 2604.13487 · v1 · submitted 2026-04-15 · 🌌 astro-ph.IM · astro-ph.EP

Site Quality Analysis for an Indian Submillimeter Telescope: A Reanalysis-Based Approach

Tanmay Singh , Mayuri Sathyanarayana Rao , Ritoban Basu Thakur This is my paper

Pith reviewed 2026-05-10 12:49 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.EP
keywords Ladakhsubmillimeter astronomyprecipitable water vaporsite selectionERA5 reanalysisatmospheric transmittance
0
0 comments X

The pith

Ladakh offers sites with low water vapor for 19-23 percent of the time, outperforming nearby Hanle and Merak.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses ERA5 reanalysis data to map precipitable water vapor across Ladakh over 184 months and identifies two grid cells where PWV stays at or below 1 mm for 19 and 23 percent of the period. These fractions exceed the 5 and 8 percent recorded at the Hanle and Merak cells. The authors then run the am radiative transfer model to calculate atmospheric transmittance and photon noise from 10 to 1000 GHz at the candidate sites. They show the resulting performance is competitive with established submillimeter facilities worldwide. The work concludes that Ladakh is favorable and recommends in-situ radiometer checks to confirm the reanalysis results before final site selection.

Core claim

Multiple ERA5 grid cells in Ladakh reach PWV ≤1 mm, with Site A near 34.25°N, 78.75°E satisfying the threshold for about 23 percent of the study duration and Site B near 32.50°N, 79.00°E for about 19 percent, compared with 5 percent at Hanle and 8 percent at Merak. Transmittance and brightness-temperature calculations indicate that these conditions support submillimeter observations at levels comparable to current or planned sites elsewhere.

What carries the argument

ERA5 reanalysis for monthly PWV statistics across Ladakh, followed by the am radiative transfer code to derive transmittance and atmospheric photon noise over 10-1000 GHz.

If this is right

  • Ladakh can support submillimeter observations with transmission windows competitive to global sites.
  • Site A and Site B outperform the existing Hanle and Merak locations within ERA5 resolution.
  • Atmospheric photon noise estimates allow direct comparison of expected data quality across frequency bands.
  • The analysis motivates targeted ground-based verification before observatory construction.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Confirmation would add a high-altitude Asian site to the global submillimeter network, potentially improving coverage for time-critical observations.
  • The reanalysis-plus-radiative-transfer method can be applied to other high-plateau regions with limited local weather data.
  • Lower PWV fractions at these sites could enable more efficient high-frequency observations than at the comparison Indian locations.

Load-bearing premise

ERA5 reanalysis at native grid resolution accurately represents actual PWV at the specific candidate locations without systematic bias from local topography or microclimates.

What would settle it

A radiometer campaign at Site A or Site B recording PWV below 1 mm for substantially less than 15 percent of the time over a full year.

Figures

Figures reproduced from arXiv: 2604.13487 by Mayuri Sathyanarayana Rao, Ritoban Basu Thakur, Tanmay Singh.

Figure 1
Figure 1. Figure 1: Comparison of precipitable water vapor (PWV) values derived from different methods with radiometer mea￾surements at Mauna Kea (top) and the ALMA site at Chaj￾nantor (bottom). Red triangles, magenta circles and blue squares correspond to the PWV values derived from the Methods M1, M2 and M3 respectively as described in this work. All other markers denote radiometer measurements compiled from L. Valeria et a… view at source ↗
Figure 2
Figure 2. Figure 2: The monthly average PWV across global sites over 2010 to April 2025 calculated using M2. Dome A has the best PWV, consistently remaining below 1 mm across all months over the full timescale. directly using the nearest or interpolated ERA5 grid￾point surface pressure in Equation 1 produces unrealis￾tically high PWV values (∼ 10 mm), inconsistent with in situ measurements (M. A. M. van Kooten & J. G. Izett 2… view at source ↗
Figure 3
Figure 3. Figure 3: Annual median PWV at the Indian Astronomical Observatory (IAO, Hanle) for 1998–2017. The plot compares radiometer measurements (IAO–Hanle; blue-white markers) with ERA5-derived PWV computed using the site-specific method based on geopotential-height–interpolated surface pressure (ERA5: IAO–Hanle, M2; orange-black markers) and the native-grid method using the nearest ERA5 grid cell (ERA5: Hanle Pixel, M3; g… view at source ↗
Figure 5
Figure 5. Figure 5: ERA5 gridded mean surface pressure map aver￾aged over the period from 2010 to April 2025, used for com￾puting precipitable water vapor (PWV) at each grid point within the region of interest. The standard deviation of sur￾face pressure at each grid is noted in each pixel. For completeness, the interpolated gridcell series (ERA5: Hanle Pixel) achieves smaller error scores against the radiometer (monthly RMSE… view at source ↗
Figure 7
Figure 7. Figure 7: Top: Monthly precipitable water vapor (PWV) time series for the candidate sites on the Ladakh plateau over 2010–2025. The lower sub-panel shows a zoomed view around the PWV ≈ 1 mm benchmark commonly used for submillimeter observing conditions. Bottom: PWV quan￾tile distributions for the same sites restricted to the winter months (NDJF) over 2010–2025. a static, plane–parallel, layered atmosphere with line … view at source ↗
Figure 8
Figure 8. Figure 8: Top: Zenith transmittance spectra, tx(f), from am for the Hanle pixel, Merak pixel, Site A, and Site B. For each site, four curves are shown: the overall median (solid blue), NDJF median (solid red), NDJF 25th percentile (dotted orange), and the APEX median profile (dashed gray) from L. Valeria et al. (2024), based on MERRA-2 reanalysis. A horizontal line at tx = 0.8 marks a reference transmission level. B… view at source ↗
Figure 9
Figure 9. Figure 9: Top: Digital elevation maps for the four candi￾date Ladakh site regions, each at ∼30 m resolution shown in a 0.5 ◦×0.5 ◦ window (±0.25◦ about the site). A white cross marks the region center and the white square corresponds to the grid size in ERA5 (0.25◦×0.25◦ ). Bottom: Probability density of surface elevation within 0.25◦×0.25◦ boxes centered on the four candidate sites. Each curve shows the normal￾ized… view at source ↗
read the original abstract

The Himalayan plateau region of Ladakh, India, is a potential host for a science-class submillimeter observatory, building on existing astronomical infrastructure near Hanle and Merak. Using the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5) data, we analyze precipitable water vapor (PWV) at monthly resolution over 184 months from January 2010 to April 2025, map PWV statistics across Ladakh, and identify candidate regions that reach PWV $\leq 1$ mm. For promising locations, we compute atmospheric transmittance and the corresponding atmospheric photon noise using the am (Atmospheric Model) radiative transfer code; we present transmittance and brightness temperature estimates over 10--1000 GHz and compare the inferred performance to sites hosting current or planned submillimeter facilities worldwide. We find Ladakh to be favorable for submillimeter observations, with multiple ERA5 grid cells reaching PWV $\leq 1$ mm. Within ERA5's spatial resolution, two regions emerge as particularly promising: Site A ($\approx 34.25^\circ$N, $78.75^\circ$E) and Site B ($\approx 32.50^\circ$N, $79.00^\circ$E), which satisfy PWV $\leq 1$ mm for about 23% and 19% of the study duration, respectively, compared to about 5% and 8% for the Hanle and Merak grid cells. These results motivate targeted in situ radiometer measurements for final site selection.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript uses ERA5 reanalysis data over 184 months (Jan 2010–Apr 2025) to map PWV across Ladakh at monthly resolution, identify grid cells with PWV ≤1 mm, and rank two candidate sites (A at ≈34.25°N, 78.75°E and B at ≈32.50°N, 79.00°E) that meet this threshold 23% and 19% of the time versus 5% and 8% for the Hanle and Merak cells. It further employs the am radiative transfer code to compute transmittance and brightness temperature over 10–1000 GHz at the candidate sites and compares performance metrics to existing or planned submillimeter facilities worldwide, concluding that Ladakh is favorable and that the two sites warrant in-situ radiometer follow-up.

Significance. If the ERA5-derived PWV statistics prove representative, the work supplies a transparent, reproducible preliminary ranking that could usefully guide site selection for a new Indian submillimeter telescope. The analysis draws on publicly documented reanalysis fields and the standard am code, includes direct comparisons to other observatories, and explicitly calls for targeted in-situ validation—strengths that increase its utility as a planning document even if the absolute fractions require later refinement.

major comments (2)
  1. [Results (PWV statistics and site identification)] The headline site-ranking claim (23% and 19% PWV ≤1 mm fractions for Sites A and B versus 5% and 8% for Hanle/Merak) is obtained directly from native-resolution ERA5 grid-cell statistics without downscaling, bias correction, or any comparison to local radiosonde/radiometer data. In the steep Ladakh terrain, sub-grid topographic and microclimate effects on scales <<31 km can alter PWV; this assumption is load-bearing for the claim that the two sites are “particularly promising.”
  2. [Methods and Results] No quantitative assessment of ERA5 PWV uncertainty or resolution limitations is provided for the high-elevation candidate cells, nor is there a sensitivity test showing how the reported percentages would change under plausible sub-grid corrections. This omission leaves the central comparative result only partially supported.
minor comments (2)
  1. [Abstract] The abstract states the PWV fractions without accompanying uncertainties or error bars; adding these (or at least noting the underlying monthly sampling) would improve clarity.
  2. [Figures and Results] Figure captions and text should explicitly state that all percentages and rankings are conditional on ERA5’s native grid resolution, consistent with the cautious phrasing already used in the abstract.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive review and for recognizing the manuscript's value as a preliminary planning document based on publicly available reanalysis data. We address each major comment below and describe the revisions we will make to improve the presentation of limitations and robustness of the results.

read point-by-point responses

  1. Referee: [Results (PWV statistics and site identification)] The headline site-ranking claim (23% and 19% PWV ≤1 mm fractions for Sites A and B versus 5% and 8% for Hanle/Merak) is obtained directly from native-resolution ERA5 grid-cell statistics without downscaling, bias correction, or any comparison to local radiosonde/radiometer data. In the steep Ladakh terrain, sub-grid topographic and microclimate effects on scales <<31 km can alter PWV; this assumption is load-bearing for the claim that the two sites are “particularly promising.”

    Authors: We agree that the reported fractions derive from native-resolution ERA5 cells and that sub-grid effects in Ladakh's terrain could influence local PWV. Our analysis is intentionally a reanalysis-based screening using publicly documented fields, as stated in the title, abstract, and conclusion; the percentages are therefore presented as indicative at the 31 km scale and are accompanied by an explicit call for in-situ radiometer validation. To address the referee's concern, we will revise the manuscript to add a dedicated paragraph in the Discussion section that (i) explicitly caveats the site-ranking claim as resolution-limited, (ii) references existing literature on ERA5 performance over complex high-altitude topography, and (iii) reiterates that the two sites are flagged for follow-up measurements rather than declared optimal. revision: yes

  2. Referee: [Methods and Results] No quantitative assessment of ERA5 PWV uncertainty or resolution limitations is provided for the high-elevation candidate cells, nor is there a sensitivity test showing how the reported percentages would change under plausible sub-grid corrections. This omission leaves the central comparative result only partially supported.

    Authors: We acknowledge that the submitted manuscript lacks a formal uncertainty quantification or sensitivity test for the PWV fractions. In the revised version we will add a short subsection (or appendix) that (i) summarizes published ERA5 PWV uncertainties at high-elevation arid sites and (ii) presents a simple sensitivity test in which the PWV time series for the candidate cells are perturbed by representative bias ranges (e.g., ±10 % and ±20 %) before recomputing the PWV ≤1 mm fractions. The resulting ranges will be reported alongside the nominal 19 % and 23 % values. We will also state clearly that a full bias correction or dynamical downscaling lies outside the scope of this reanalysis-focused study and is left for future work once in-situ data become available. revision: yes

standing simulated objections not resolved
  • Direct comparison of the ERA5-derived PWV statistics to local radiosonde or radiometer measurements at the specific candidate grid cells, as no such long-term in-situ observations exist for these remote locations.

Circularity Check

0 steps flagged

No circularity; direct statistical extraction from external ERA5 reanalysis and am code.

full rationale

The paper performs straightforward monthly PWV statistics and transmittance calculations on ERA5 grid data (184 months) plus external am radiative transfer modeling, with no equations, fitted parameters, self-referential predictions, or load-bearing self-citations. All reported fractions (e.g., 23% and 19% for Sites A/B) are direct counts from the input dataset at native resolution; the derivation chain contains no reductions to its own outputs or prior author work.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis rests on the domain assumption that ERA5 reanalysis faithfully represents local PWV at the stated grid cells and that the 1 mm threshold is the appropriate figure of merit; no free parameters are fitted and no new entities are postulated.

axioms (1)
  • domain assumption ERA5 reanalysis data accurately captures precipitable water vapor at the spatial and temporal resolution used for site selection.
    The entire site-ranking result depends on this assumption; the abstract does not report any local validation.

pith-pipeline@v0.9.0 · 5597 in / 1419 out tokens · 37555 ms · 2026-05-10T12:49:54.521305+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

    work page

  3. [3]

    arXiv.org , author =

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page doi:10.3847/1538-3881/aabc4f 2017