The rain gauge is still the primary instrument, but the network is sparse, concentrated in developed countries, and carries measurement error and poor spatial representation.
IMERG and similar products face near-real-time aggregation and inter-calibration challenges across a heterogeneous satellite constellation, plus high operational cost.
SM2RAIN and SM2RAIN-NWF both require reference rainfall at spatial and temporal resolution comparable to the soil-moisture data, just to fix Z and a. That reliance constrains the method to gauged regions.
You need rainfall observations to build a model whose purpose is estimating rainfall where observations don't exist. This paper breaks that circle by transferring parameters from hydroclimatically analogous sites.
EvapotranspirationSoil moisturePrecipitationWindTemperatureRunoffSnowPressureVegetationLakeRadiationLat / Lon
Iterative Variance Inflation Factor screening, target VIF < 10, while keeping every process group represented. All variables passed except two, kept deliberately so their categories stayed in the feature space:
t2m · VIF 23.27str · VIF 13.79 → 18 variables retained for UMAP and clustering.
A static soil-properties map was not an explicit feature; its influence enters through land-surface response variables (long-term volumetric soil water, runoff, evaporation) that integrate soil, geology and vegetation effects on the regional water balance.
Hourly soil moisture and precipitation records; stations assigned to clusters by geographic location. Station counts differ by measurement depth:
Global coverage is dense over North America and Europe, with smaller groupings in West Africa, East Africa, eastern China and south-eastern Australia (Fig. 1 of the paper).
For POD and FAR, daily precipitation is first converted to a binary rain / no-rain series with a 1 mm day⁻¹ threshold — any day above it (in situ or SM2RAIN-NWF) counts as a rain event.
Two independent justifications converge: the Köppen–Geiger classification contains approximately 30 global climate subtypes, and GNG — which does not require k a priori — independently identified 30. The choice is grounded both climatologically and in the data.
K-means finds spherical clusters efficiently; GNG grows an adaptive topology capturing non-linear structure; GMM allows probabilistic, ellipsoidal assignment; AC reveals hierarchical nesting. Using all four tests whether the conclusion depends on the clustering method — it does not.
LOOCV — for a cluster with N stations, recompute the cluster mean N times, each time excluding one station, then apply to the held-out station.
20% hold-out — split each cluster 80/20, derive (ā, Z̄) from the training subset with B = 10,000 bootstrap resamples for 95% CIs, apply to the untouched test set.
Under the assumption that evapotranspiration is negligible during rainfall. A calibration procedure estimates all four parameters Z, a, b, c — and that procedure is exactly what this paper removes.
| Z | depth-related soil parameter |
| θ(t) | volumetric soil moisture |
| P, R, ET, G | rainfall, runoff, evapotranspiration, drainage |
Only Z and a. The fluctuation of c was found negligible — a stable value of 15.90 ± 0.005 across all stations — so c was held constant and the investigation restricted to the two sensitive parameters.
SM2RAIN-NWF replaces empirical drainage with the net water flux derived by inverting Warrick's (1975) solution of the linearized Richards equation. Reported to improve rainfall estimation accuracy by ~30% over SM2RAIN at basin (Lake Urmia) and national (Iran) scales.
| T | dimensionless time |
| U | dimensionless function from Warrick's solution |
| N, i | number of discrete time steps, and each step within the total N |
K-means sweeps all three — highest silhouette (0.43), lowest Davies–Bouldin (0.82), highest Calinski–Harabasz (197,736). GNG is the clear runner-up (0.40 / 0.86 / 190,865) and beats GMM and AC on Calinski–Harabasz, so both were carried into the calibration analysis. GMM is weakest on every measure.
The 3-D UMAP embedding separates the ERA5-Land feature space into 30 well-defined groups for both K-means and GNG (Fig. 2a, 2c). Mapped back onto the globe (Fig. 2b, 2d), those groups form coherent geographic regions of shared hydroclimatic behaviour — which is precisely what makes a shared calibration coefficient defensible.
GNG produces clusters that differ from partition-based methods, reflecting subtle non-linear patterns — some clusters (e.g. 21) show distinctly wider parameter spreads, others (e.g. 2) much narrower, indicating more homogeneous soil or climatic properties.
Only clusters with sufficient in-situ station coverage were retained. Clusters lacking adequate observational support were excluded to avoid unstable parameter estimation — so the framework needs a minimum of observations somewhere inside each zone, and is not yet directly transferable to regions with none.
This table is the product of the paper. Give a new, ungauged site its cluster and its sensor depth, read off (Z̄, ā), and SM2RAIN-NWF runs — no reference rainfall required. Z generally ranges higher in deeper layers (e.g. K-means cluster 23 at 10–20 cm, cluster 11 at 20–30 cm), whereas a may rise, fall, or turn bimodal depending on cluster composition.
Z is broadly distributed at low elevations; a shows a sharp, highly concentrated peak at low values in the same elevation band — strongly constrained behaviour near the surface.
Diversity of calibration coefficients is highest here: a broad, near-uniform distribution across clusters.
The distribution begins to concentrate at lower values, suggesting the parameter is becoming more constrained by underlying soil texture rather than diverse surface-level influences.
The low-elevation distribution of Z becomes distinctly bimodal. Evapotranspiration, infiltration and soil layering likely play a greater combined role in setting the coefficients.
The primary low-elevation mode consolidates toward the middle of the parameter range and a separate mode of lower Z appears at high elevations. The sharp surface peak in a has become broad and diffuse.
Many clusters show scattered patterns in both Z and a, possibly from significant changes in infiltration.
Wider violins and multi-peaked distributions (e.g. cluster 12 in some depth ranges, cluster 21 across multiple depths under GNG) indicate greater heterogeneity among stations within one cluster — the same zone can contain very different calibration behaviours. Narrow violins with tight boxes (cluster 4 at the surface; cluster 10 at shallow and mid depths; cluster 2 under GNG) mean more consistent station characteristics and more robust mean estimates.
The authors note a second, methodological contributor: Kim and Crow (2024) show that simplifying the water-balance physics and the limitations of soil-moisture data — noise and low sampling frequency — can themselves drive apparent variation in Z. The broadening seen in the deeper-layer clusters (notably 23 and 11) is probably partly this.
Some clusters show minor local degradation, but a significant number improve on the fully calibrated baseline. Most evident in correlation at 5–10 cm:
Why improvement is possible: the cluster-average parameter is a more robust and generalizable estimator — it is less prone to overfitting the noise and idiosyncrasies of a single site's calibration record.
Performance shifts are minimal across all depths and clusters. Changes in R are typically less than ±4%, and POD is almost perfectly preserved — 0.00% to 0.85% change in most cases.
Together with GNG's gains, this is the paper's central conclusion: cluster-based parameterization is a practical, robust pathway for scalable deployment in data-limited regions.
Calibrated runs show systematically slightly lower RMSE and reduced bias, but the differences are not large. FAR and POD show no significant difference between the calibrated (green) and cluster-parameterized (red) cases in Figs. 6 and 8 — event detection survives intact.
Both were carried through the same hold-out design (Tables A4, A5). They show the same qualitative result — cluster-derived coefficients hold up — but with larger negative excursions in R. Note the paper reports these tables without interpretive discussion; the reading below is drawn directly from their values.
The paper's own method-independence claim rests on K-means and GNG, whose results it calls uniform across both frameworks. Tables A4 and A5 point the same way, with thinner margins.
From point-wise calibration to regionalized calibration: each region can be assigned a fixed — or narrowly distributed — parameter set for a and Z, instead of calibrating the model separately at every station.
It opens SM2RAIN-NWF to remote and data-limited areas, and supports scalable climate monitoring, reanalysis and prediction.