Multifractal Description of Streamflow and Suspended Sediment Concentration Data from Indian River Basins

Adarsh Sankaran; Anand Vishnu Babykuttan; Drisya Sasi Dharan; Marcin W\k{a}torek; Nandhu Ambika Raju; Vysakh Kunju Kunju

REVIEW 4 major objections 6 minor 11 references

Multifractal Description of Streamflow and Suspended Sediment Concentration Data from Indian River Basins

T0 review · 4 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Analysis of 192 daily streamflow records from 13 Indian river basins finds multifractality and long-term persistence with a mean Hurst exponent of about 0.585, and a joint flow–sediment persistence that is nearly the average of the…

desk verdict First India-wide multifractal streamflow baseline using standard tools; the quantitative exponents are likely biased by zero-stretch handling, but the descriptive compilation is still useful. read the letter →

arxiv 1909.01605 v1 pith:J3XROX5S submitted 2019-09-04 physics.geo-ph physics.ao-phphysics.data-anstat.AP

Adarsh Sankaran , Drisya Sasi Dharan , Anand Vishnu Babykuttan , Nandhu Ambika Raju , Vysakh Kunju Kunju , Marcin Wątorek This is my paper · ORCID

classification physics.geo-phphysics.ao-phphysics.data-anstat.AP

keywords streamflowmultifractalityHurstexponentlong-termpersistencesuspendedsedimentcross-correlationMFCCAIndianriverbasins

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper analyzes daily streamflow records from 192 stations across 13 Indian river basins using multifractal detrended fluctuation analysis (MF-DFA) and reports that the series are multifractal and long-term persistent, with a mean Hurst exponent of about 0.585. It then pairs streamflow with total suspended sediment data at 95 stations in five basins, using multifractal cross-correlation analysis (MFCCA). The central empirical finding is that the joint persistence of streamflow and suspended sediment is close to the arithmetic mean of the two individual persistence values. If these results hold, Indian river flows follow the long-memory, multifractal behavior observed in other regions, and sediment-transport memory at a station can be summarized from the flow and sediment scaling exponents alone.

What carries the argument

The method that carries the argument is multifractal detrended fluctuation analysis (MF-DFA), which detrends the cumulative profile in segments, computes $q$-th order fluctuation functions over a range of scales, and estimates a generalized Hurst exponent $h(q)$ from their log-log slope; q-dependence of $h(q)$ is the signal of multifractality. The second machine is multifractal cross-correlation analysis (MFCCA), which extends the same detrending to pairs of series and yields a cross-correlation exponent $\lambda(q)$ and a scale-dependent detrended cross-correlation coefficient. The identity that carries the coupling result is the arithmetic-mean rule $H_{xy} \approx (H_x + H_y)/2$, which the paper checks station by station.

What would settle it

Run a controlled test on synthetic series with a known Hurst exponent, insert zero stretches of the lengths found in peninsular Indian rivers, apply the paper's exact scale choices (minimum scale 10 for MF-DFA; minimum scale exceeding the longest zero run for MFCCA), and check whether the estimates recover the prescribed exponents and the mean-of-individuals relation. If they do not, the reported persistence values and basin rankings are artifacts of the zeros.

Watch

Extended reading notes

Core claim

The paper's core claim is that Indian streamflow is neither random nor single-scaling: the generalized Hurst exponent $h(q)$ varies with $q$, so the records are multifractal, and $h(q=2)$ sits mostly above 0.5, indicating long-term persistence; the basin-wide mean is 0.585, lower than a previously reported global runoff mean of 0.73. A second claim concerns coupling: using MFCCA, the joint persistence exponent $H_{xy}$ for streamflow and total suspended sediment is approximately $(H_x + H_y)/2$, the mean of the individual series exponents, across all five basins studied. A third comparative claim is that Krishna basin streamflow has the least persistence, Godavari the strongest multifractality and complexity (attributed to extensive flow regulation), and that streamflow persistence exceeds sediment persistence at about 60 percent of stations, with Godavari the reverse.

Load-bearing premise

The entire analysis assumes that the Hurst and cross-correlation exponents estimated by MF-DFA and MFCCA are unbiased for intermittent streamflow series with long runs of zero flow; the paper sets the minimum scale above the longest zero run but does not validate this on synthetic intermittent series.

Editorial extensions

If this is right

Daily streamflow in these basins can be treated as a long-memory process, so flood events are more likely to be followed by flood events rather than being independent across days.
Joint flow–sediment persistence at a station can be estimated from the two single-series Hurst exponents without running the cross-correlation analysis.
Annual-scale flow–sediment correlation is generally stronger than the overall correlation, so averaging over all scales hides the seasonal coupling that actually drives sediment transport.
Godavari's strongest multifractality and reversed flow–sediment persistence ordering indicate that dams and flow regulation alter the scaling signature of sediment transport.
The Indian mean persistence of 0.585 is below the global runoff value of 0.73, so extrapolating global persistence estimates to Indian basins would overestimate memory strength.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Inference: if the averaging rule for joint persistence generalizes, it suggests a one-parameter shortcut for paired flow–sediment records: estimate the two single-series exponents and take their mean as the coupled memory, which could be tested on new basins with paired data.
Inference: because the zero-run bias question is the weakest point, a synthetic intermittent-streamflow benchmark would either confirm the paper's scale choices or reveal which basin rankings change, making it the most direct next step.
Inference: the stronger annual than overall correlations imply that reservoir siltation and sediment-load models should use seasonal, scale-dependent coupling rather than a single annual aggregate.
Inference: comparing multifractal spectra of regulated versus unregulated rivers could separate human intervention effects from climatic controls on streamflow complexity, extending the paper's attribution of Godavari's behavior to regulation.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

Quick take: this is the first India-wide multifractal dataset for streamflow and streamflow-sediment links, built with standard tools. It's descriptive, not methodologically new, and it has a couple of real weaknesses. But as a baseline it's worth keeping.

The new part: 192 stations across 13 basins, six multifractal parameters, and MFCCA for five basins with sediment. That is a genuinely useful compilation for anyone working on Indian water resources. The authors are appropriately cautious in the discussion about dam effects and intermittency. The figures and tables are consistent with the claims at face value, and the citation pattern is solid, including the key prior results from Kantelhardt and Koscielny-Bunde.

The soft spots, in rough order of severity. First, the zero-stretch issue: southern peninsular rivers are intermittent, with long runs of zero or near-zero discharge. MF-DFA uses a fixed minimum scale of 10 days, while MFCCA uses a minimum scale exceeding the longest zero run. That asymmetry is unexplained, and no synthetic or surrogate test shows how much bias zero stretches produce in the estimated Hurst exponents. For a daily series with a long dry spell, a 10-day window can be entirely zeros; after profile construction that becomes a linear segment that the detrending polynomial removes, which can inflate the scaling exponent. Since the headline mean of 0.585 and the basin rankings all rest on these exponents, this is the load-bearing concern. It may not overturn the qualitative conclusion—Indian streamflow does look persistent and multifractal, consistent with global results—but the exact numbers need error bars and a sensitivity check. Second, the "joint persistence is the mean of individual persistences" result is not an empirical discovery. At q=2 the detrended covariance scales as the product of the individual variances, so Hxy ≈ (Hx+Hy)/2 is built into the method. The authors present this as a finding; a sentence acknowledging the mathematical origin would be honest. Third, details: the detrending polynomial order is never stated, and the abstract says mean 0.585 while the conclusions say 0.583. Minor, but sloppy. Fourth, no uncertainty quantification anywhere; basin comparison PDFs are shown without confidence intervals.

Who this is for: hydrologists wanting a quick India-wide reference on streamflow persistence and multifractality, or a comparison set for future method papers. With revisions—surrogate validation, explicit error bars, fixed text, and an honest note about the joint-persistence relation—it would be a solid regional contribution. I'd send it to review, not desk reject.

Would I cite it? Only if I needed the Indian statistics.

Referee Report

4 major / 6 minor

Summary. The paper applies Multifractal Detrended Fluctuation Analysis (MF-DFA) to daily streamflow records from 192 stations in 13 Indian river basins, and Multifractal Cross-Correlation Analysis (MFCCA) to streamflow–total suspended sediment (TSS) pairs at 95 stations in five basins. It reports that Indian streamflow is multifractal and exhibits long-term persistence with a mean Hurst exponent of about 0.58, that the Krishna basin is least persistent while the Godavari basin is most multifractal and complex, that the joint persistence of streamflow and TSS is approximately the mean of the individual persistence exponents, and that annual cross-correlations are generally higher than seasonal ones.

Significance. If the reported exponents were unbiased, the paper would provide a useful regional confirmation of the multifractality and long-term persistence seen in global runoff records, and it would extend MFCCA to a sediment-transport setting, which is comparatively rare. The compilation of 192 stations across 13 basins is a valuable empirical contribution. However, the headline quantitative claims—the mean Hurst exponent, the basin rankings, and the joint-persistence rule—rest on estimation choices and statistical hypotheses that are not validated in the manuscript, so the significance of the results as reported is presently uncertain.

major comments (4)

[Study area and Data; Results and Discussions (MF-DFA)] The paper acknowledges that peninsular river stations are intermittent with continuous zero or very low discharge, yet the MF-DFA analysis fixes the minimum scale at 10 days while the MFCCA analysis sets the minimum scale above the longest zero stretch. No synthetic or surrogate test is provided to show that the fixed 10-day minimum does not bias Hurst exponents for series containing multi-month zero runs. Because the mean exponent of 0.585, the 71.3% persistence percentage, and the basin rankings all derive from these estimates, this is a load-bearing gap. I ask for controlled experiments on synthetic intermittent series with known Hurst exponents (for example, fractionally integrated noise masked by dry periods calibrated to the observed zero-run lengths) and for a comparison of the two minimum-scale choices, together with bootstrap confidence intervals for H.
[Materials and Methods, Eqs. (13)-(14); Tables 2-6] The result that the joint persistence is approximately the mean of individual persistences is not an independent empirical law in the form presented. For q=2, F_2^{XY}(s) is the average over segments of products of detrended fluctuations of the two series, so if the DCCA correlation coefficient is approximately scale-invariant, the scaling exponent of F_2^{XY} is forced to be close to (H_x+H_y)/2 by the construction of the estimator. Presenting this as a discovery requires either a null-model comparison (for example, independent surrogates with the same individual H, showing that observed Hxy deviates from (Hx+Hy)/2 beyond what the estimator produces) or an explicit report of the deviations with their uncertainties. As written, Tables 2–6 cannot support the claim that this is a physical property of the streamflow–TSS system.
[Materials and Methods (scale selection); Results and Discussions (MF-DFA)] The Methods state that the maximum scale should be below 1/10 of the sample size, while the Results state that the maximum scale is N/2. This is an internal contradiction, and the actual scale range used for the reported exponents is ambiguous. Furthermore, the minimum-scale rule differs between MF-DFA (fixed at 10 days) and MFCCA (above the longest zero stretch) with no justification. The manuscript should specify the exact scale ranges used for each station and demonstrate that the fitted exponents are insensitive to the choice, or present results under both rules. Since all basin comparisons depend on these fitted exponents, the ambiguity is not merely cosmetic.
[Results and Discussions, Fig. 4] Comparative statements such as 'Krishna has least persistence' and 'Godavari has strongest multifractality' are based on overlapping PDF/CDF curves in Fig. 4, with no error bars or significance tests reported. Because the exponents are estimates with substantial sampling variability, a formal comparison—for example, bootstrap confidence intervals for H and spectral width, or a two-sample test on station-level exponents—is needed before these basin rankings can be accepted.

minor comments (6)

[Abstract vs Conclusions] The mean Hurst exponent is given as 0.585 in the Abstract and 0.583 in the Conclusions; the correct value should be reported consistently along with the exact number of stations used in its calculation.
[Table 1 vs Tables 2-6] The number of stations differs between Table 1 and the MFCCA tables (Godavari 23 vs 26, WFR T-K 28 vs 19, Krishna 31 vs 23, Mahanadi 19 vs 16). Please clarify whether these are the numbers of stations with available sediment data and, if so, state this explicitly in the text.
[Materials and Methods (MFCCA)] The statement that 'MFCCA is retrieved for the moment order q=2' contradicts the later statement that moment orders from -4 to +4 are used, and the figures show q-dependent H(q) for the joint series; the operative protocol should be clarified.
[Results and Discussions (MFCCA)] The meaning of 'seasonal correlation' (90-day or 3-month scale) and 'annual correlation' should be defined explicitly in the Methods, since these terms are used for all five basins.
[Results and Discussions, Fig. 5] The statement that drainage area has no effect on persistence should be supported by a correlation coefficient or regression slope rather than visual inspection of the scatter plot.
[Figure captions] The typo 'Assymetry' appears in the captions of Figs. 2 and 4 and should be corrected to 'Asymmetry'.

Circularity Check

1 steps flagged · score 6.0 of 10

MFCCA 'joint persistence equals mean of individual persistences' is built into the q=2 estimator, not an empirical discovery.

self definitional [Materials and Methods, MFCCA (around Eq. 13, q=2 note); Abstract; Results and Discussions, Tables 2-6]
"It is to be noted that in this study MFCCA is retrieved for the moment order q=2. ... At all stations of Cauvery basin, the joint persistence is found to be nearly the average of individual persistence of streamflow and TSS."

At q=2, Eq. (13) defines the MFCCA cross fluctuation function as F_2^XY(s) = [mean |f_XY^2(s,nu)|]^(1/2), where f_XY^2 is the segment detrended covariance. For two scaling series with detrended variances scaling as s^(2Hx) and s^(2Hy), the covariance term scales as s^(Hx+Hy), so the slope fitted from log F_2^XY vs log s is lambda(2) = (Hx+Hy)/2. The paper estimates Hx and Hy from the same MF-DFA/MFCCA machinery and then reports Hxy as an independent result; e.g., Biligundulu: (0.797+0.669)/2 = 0.733 equals the tabulated Hxy = 0.733. Thus the headline 'joint persistence is approximately the mean of the persistence of individual series' is a restatement of the q=2 estimator's built-in scaling, not a new property of Indian streamflow-sediment data.

full rationale

This is a partial circularity. The MF-DFA streamflow characterization (mean exponent 0.585, multifractal widths, basin ordering, insensitivity to drainage area) is a standard estimator applied to external data and does not reduce to its own inputs; the zero-stretch and scale-selection concerns are correctness risks, not circularity. The load-bearing step that reduces by construction is the cross-correlation conclusion Hxy approximately equals (Hx+Hy)/2. With q=2 the MFCCA fluctuation function is the square root of a mean detrended covariance, whose scaling exponent is the sum of the individual exponents under the standard scaling assumption, so the average relation is inherent to the estimator. Since this relation is advertised in the abstract and was treated as an empirical finding for every basin, the score is 6 rather than 0. No self-citation chain or imported uniqueness theorem is involved.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The analysis rests entirely on standard, cited algorithms and on the public WRIS data. The only ad hoc elements are the scale and moment-order choices, and the implicit assumption that intermittent zero-heavy series can be analyzed without bias. No new physical entities or parameters are introduced.

free parameters (4)

Minimum scale for MFCCA = greater than longest zero stretch (station-specific)
Hand-chosen threshold to avoid zero-dominated segments; not systematically varied or justified.
Scale range for MF-DFA = 10 to N/2
Ad hoc bounds; no sensitivity analysis is provided.
Moment order range q = -4 to +4
Chosen by hand following common practice; no convergence or sensitivity analysis.
Detrending polynomial order m = not reported (stated as 1-3 in text)
The paper says order 1-3 is sufficient but never states the value used in the reported results, leaving the analysis under-specified.

assumptions (4)

standard math MF-DFA and MFCCA algorithms, as introduced by Kantelhardt et al. (2002) and Oświęcimka et al. (2014), are valid for the analyzed series.
The paper applies these methods without modification or formal proof, relying on the cited literature.
domain assumption The WRIS streamflow and TSS data are accurate, homogeneous, and long enough for multifractal analysis.
No quality control or homogeneity testing is reported; the analysis assumes the public dataset is reliable.
ad hoc to paper Intermittent series with long zero stretches still yield reliable Hurst exponents after setting the minimum scale above the longest zero run (MFCCA) or at 10 (MF-DFA).
The scale choice is introduced specifically for this analysis and is not validated with surrogate or synthetic tests.
domain assumption The detrending polynomial order (DFA1-3) does not materially affect the reported exponents.
The paper discusses polynomial order as a key issue but does not state which order was used in the reported results.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Multifractal Description of Streamflow and Suspended Sediment Concentration Data from Indian River Basins." pith.science (2026). https://pith.science/paper/J3XROX5S

@misc{pith2026190901605,
  author       = {Pith},
  title        = {Pith review of: Multifractal Description of Streamflow and Suspended Sediment Concentration Data from Indian River Basins},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/J3XROX5S}},
  note         = {Machine review of arXiv:1909.01605}
}

read the original abstract

This study investigates the multifractality of streamflow data of 192 stations located in 13 river basins in India using the Multifractal Detrended Fluctuation Analysis (MF-DFA). The streamflow datasets of different river basins displayed multifractality and long term persistence with a mean exponent of 0.585. The streamflow records of Krishna basin displayed least persistence and that of Godavari basin displayed strongest multifractality and complexity. Subsequently, the streamflow-sediment links of five major river basins are evaluated using the novel Multifractal Cross Correlation Analysis (MFCCA) method of cross correlation studies. The results showed that the joint persistence of streamflow and total suspended sediments (TSS) is approximately the mean of the persistence of individual series. The streamflow displayed higher persistence than TSS in 60 % of the stations while in majority of stations of Godavari basin the trend was opposite. The annual cross correlation is higher than seasonal cross correlation in majority of stations but at these time scales strength of their association differs with river basin.

Figures

Figures reproduced from arXiv: 1909.01605 by the authors.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

[1]

Krishna 31 1850 251360 1095 18615

work page
[2]

Brahmani-Baitarani 9 830 33955 4015 15330

work page
[3]

Sabarmati 6 1421 19636 5840 9490

work page
[4]

Mahi 7 1510 32510 3285 13805

work page
[5]

Mahanadi 19 1100 124450 4015 15695

work page
[6]

Subarnarekha 5 1330 12649 6205 14235

work page
[7]

Tapi 5 8487 58400 3285 5110

work page
[8]

Cauvery 31 258 66243 2555 16425

work page

Show all 11 references

[9]

WFR Tadri-Kanyakumari 28 238 5755 1460 16425
[10]

EFR Pennar-Kanyakumari 13 850 16230 4015 16060
[11]

Godavari 23 2500 307800 1019 13111 12 Pennar 7 2486 37981 1245 10606 13 WFR-Kutch- Saurashtra-Luni 8 345 6960 6865 15111 22 Table 2 Hurst exponents of streamflow and TSS data of Cauvery basin along with the cross correlation Station Hx (Streamflow) Hy (TSS ) Scaling Exponent (...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

Krishna 31 1850 251360 1095 18615

work page

[2] [2]

Brahmani-Baitarani 9 830 33955 4015 15330

work page

[3] [3]

Sabarmati 6 1421 19636 5840 9490

work page

[4] [4]

Mahi 7 1510 32510 3285 13805

work page

[5] [5]

Mahanadi 19 1100 124450 4015 15695

work page

[6] [6]

Subarnarekha 5 1330 12649 6205 14235

work page

[7] [7]

Tapi 5 8487 58400 3285 5110

work page

[8] [8]

Cauvery 31 258 66243 2555 16425

work page