the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Global dynamic precipitation isoscapes over three-quarters of a century
Abstract. Stable precipitation isotopes are widely used as tracers in water cycling and material transport. However, observational networks (e.g., the Global Network of Isotopes in Precipitation, or GNIP) are sparse in space and time, constraining analyses and applications in data-poor regions. To bridge this gap, static isoscapes, which interpolate monthly observations to long-term annual or climatological monthly means, have been developed and applied across disciplines. In recent years, growing interest in event-scale processes and isotope-enabled hydrological modeling has increased the need for higher-frequency information.
In this study, a global precipitation isoscape using an offline isotope circulation model (ICM) forced by JRA-3Q on a 1.25° grid covering −80° to 80° of latitude is produced. The deliverables include daily values, precipitation-weighted monthly values, and climatological monthly means for 1948–2023; the daily and monthly time series cover September 1947 to March 2024. GNIP-based validation demonstrates high skill for climatology. When all the station–month pairs are analyzed in a single regression, R2 equals 0.86 for δ18O and 0.87 for δ2H, with root-mean-square errors (RMSEs) of 2.02 ‰ and 16.0 ‰, respectively. At the monthly scale, unweighted averaging of station-level metrics yields R2 values of 0.47 and 0.48, with RMSEs of 2.62 ‰ and 19.9 ‰, respectively. Daily performance was quantified from same-day regional means in East Asia and from 21 GNIP stations not used in the correction; R² is typically between 0.3 and 0.6, and the RMSE is between 2 and 4 ‰ for δ18O. By comparison, d-excess exhibits lower global skill, likely reflecting sensitivity to uncertainties in the forcing and to model simplifications.
In contrast to observation-interpolated static isoscapes, this reanalysis-driven dataset delivers seamless coverage, including data-sparse regions, with explicitly characterized error properties. The data are distributed as NetCDF/CSV (daily, monthly) and GeoTIFF (climatology) and support applications from isotope-enabled hydrology and source attribution to climate-impact assessments and global water-resource analyses.
- Preprint
(5797 KB) - Metadata XML
-
Supplement
(2624 KB) - BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on essd-2026-52', Anonymous Referee #1, 05 Apr 2026
-
RC2: 'Review of "Global dynamic precipitation isoscapes over three-quarters of a century" by Ikuya Adachi and Tsutomu Yamanaka', Anonymous Referee #2, 09 Jul 2026
The authors present a global product of isotopes in precipitation using an isotopic model driven by reanalysis data.
Consistent fields of isotopes in precipitation can be very valuable for a large variety of users from climatology, hydrology, ecology, to plant physiology. With this in mind, I would suggest a variety of amendments to the dataset.
- Isotopes in precipitation alone are not so useful without other information such as other environmental variables. They come from the reanalysis product JRA-3Q. I do not know how to access this data and I could not find out easily. It looks like it is only available to collaborative institutions. A solution would be to include a subset of variables in this data product. One could follow the subset that ECMWF chose for the timeseries downloads (https://cds.climate.copernicus.eu/datasets/reanalysis-era5-single-levels-timeseries?tab=overview).
- It is also important that I can do my own averages as a user. For example, I want to have averages per season. With the current product, this is only possible with the netcdf product because the csv product does not include precipitation. But I also do not know what is precipitation in the product: mean, sum, ...?
- In plant physiology, for example, isotopes in precipitation is not enough. I would need also isotopes in atmospheric vapour. This should be available from the isotope model as well so one could include it in the product.
- For quite a few applications, sub-daily information is important. Leaf water isotopes are, for example, very different in midday and at night. Other people started using the isotopic composition of single events such as storms, which a reanalysis-driven product could probably provide. The hourly product would hence be very valuable. I would rather have the hourly than the daily product.
- An important point of the current dataset is for me the physical consistency, which is more difficult with pure data-driven approaches. But in this way, the current bias correction poses problems. Especially once vapour isotopes are included, one cannot simply do a bias correction of the precipitation isotopes. A solution would be that corrected and uncorrected data is provided.
- I do not think that the current bias-correction is appropriate. By doing an inverse distance-weighted correction, there is a correction if a GNIP station is nearby but no correction if there is no GNIP station nearby. But the latter has probably the same bias than the former. The current bias correction does also not take into account known dependencies of water isotopes such as the altitude effect and continental effects. I would imagine rather a bias correction with an external drift such as kriging.
- I was also a bit surprised that IsoGSM was not mentioned, which provides a similar data product. It would be interesting to know the differences, the strength and weaknesses of both products so that I can choose the appropriate product for my needs.
- The author compared to SWING2 results but the latest results of this intercomparison are called WisoMIP now (Bong et al., Journal of Geophysical Research: Atmospheres 2026, https://doi.org/10.1029/2025JD044985)
- Units are missing in all products.
- I find the csv product not very useful. It would of course be fantastic if I only have to download the grid cell I am interested in but currently I have to download a zip file for a large area from which I then have to extract the grid cell.
- The description of the data files is insufficient in the manuscript and also on the Zenodo sites. For example, I had to find out by trial and error which zip file I have to download for my specific station.
Citation: https://doi.org/10.5194/essd-2026-52-RC2 - AC1: 'Comment on essd-2026-52', Ikuya Adachi, 02 Aug 2026
Data sets
Global climatology, monthly, daily precipitation isotope dataset, 1947–2024 Ikuya Adachi and Tsutomu Yamanaka https://doi.org/10.5281/zenodo.18309766
Globa daily precipitation isotope dataset, 1947–2024 Ikuya Adachi and Tsutomu Yamanaka https://doi.org/10.5281/zenodo.17337279
Global monthly precipitation isotope dataset, 1947–2024 Ikuya Adachi and Tsutomu Yamanaka https://doi.org/10.5281/zenodo.17301045
Model code and software
ICM.pro Ikuya Adachi, Tsutomu Yamanaka https://doi.org/10.5281/zenodo.18309766
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 293 | 143 | 35 | 471 | 70 | 23 | 42 |
- HTML: 293
- PDF: 143
- XML: 35
- Total: 471
- Supplement: 70
- BibTeX: 23
- EndNote: 42
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This manuscript presents a valuable global dynamic precipitation isotope dataset with a long temporal coverage, multiple temporal products, and a generally well-structured validation framework. Overall, the manuscript is of high quality, and the dataset has clear scientific value and broad potential applications, especially for isotope-related studies in data-sparse regions. In my view, the paper is already close to publishable in its current form.
I only suggest a few minor revisions: