Hydrology · Energy · Quantum computation

Water, energy and computation in one open-system view.

ShakyaQX is a research theme for Earth-system science. Its foundation is the SQX manifold, a peer-reviewed, three-dimensional generalisation of the Budyko hypothesis for open hydrological systems. The theme builds from there towards coupled water–energy behaviour, land-surface information and quantum-enhanced computation.

Developed by Chandra Mohan Shakya (Ph.D. Scholar) under the supervision of Prof. Rajib Kumar Bhattacharjya,
Department of Civil Engineering, Indian Institute of Technology Guwahati.
Published in Journal of Hydrology, Vol. 679 (2026), Article 136406.
Shakya QX HYDROLOGY ENERGY QUANTUM
The ShakyaQX theme links three areas around the SQX core: open-system hydrology, coupled water–energy behaviour and quantum-enhanced computation.
The ShakyaQX theme

More than one paper: a research programme

The published SQX manifold is the first layer. ShakyaQX is meant to grow from it into a connected framework for water, energy and computation in open Earth systems. Each layer below is labelled honestly by its current status.

01
Published foundation

Open-system hydrology

The SQX manifold makes basin openness an explicit coordinate of the Budyko framework. This covers glacier melt, groundwater depletion, reservoir regulation and inter-basin transfers, which classical theory assumes away.

02
Open preprint · developing

Coupled water–energy

A dual Budyko view that links how water is partitioned with how surface energy is partitioned, so a basin's hydrological state and thermal state can be read together.

03
Patent pending

Quantum-enhanced computation

Quantum–classical hybrid methods for calibrating and estimating parameters in Budyko-based models, developed under the filed patent application.

04
Future direction

Land surface & decision tools

Forward estimation of openness from land-use and land-cover information, together with spatial (GIS) tools that make the framework usable for water planning.

The framework

From a curve to a manifold

Classical Budyko-type equations describe how precipitation splits into evaporation and runoff in a closed catchment at steady state. Real basins are often open: groundwater is pumped or recharged, storage drains after the monsoon, and water moves across boundaries. SQX makes that openness an explicit state variable instead of treating it as error.

The SQX partition equation
E/P = x·[ 1 + (x/φ)γ ]−1/γ
with the generalised water balance  P + ΔS = E + Q
xAridity index, PET / P (the climatic coordinate) φOpenness coordinate, φ = 1 + ΔS/P  (φ > 0). φ = 1 for a closed basin γLandscape response parameter that indexes the family of surfaces E/PEvaporative index; runoff follows as Q/P = φ − E/P
1

Reduces to the classical form

At φ = 1 the surface becomes the Mezentsev–Choudhury–Yang (MCY) Budyko curve, so all classical results are kept as a special case.

2

Physically bounded

Evaporation is limited by energy (E/P → x) in humid climates and by available water, now including storage supply (E/P → φ), in arid ones.

3

Structurally stable

Strictly monotonic and globally concave in x. The input-rescaling singularity is moved to the boundary of the admissible domain (φ → 0) rather than appearing inside it. One continuous equation covers storage-subsidized and storage-depleted states, with no piecewise forms.

4

A stable climatic coordinate

Openness enters as a separate axis, so the aridity index keeps its climatic meaning instead of being rescaled by storage terms.

Sign convention. SQX treats ΔS as a net subsidy to the active water balance. ΔS > 0 is a net water gain (for example glacier melt or groundwater inflow), so φ > 1 is storage-subsidized. ΔS < 0 is depletion (for example groundwater extraction or drought), so φ < 1 is storage-depleted. This differs from the usual storage-anomaly convention, so φ should not be read as a direct restatement of GRACE storage change.

012340.00.40.81.2Aridity index x = PET / PE / Penergy limitφ = 0.70φ = 0.85φ = 1.00φ = 1.15φ = 1.30
One equation, many basins

The SQX family of curves

Each curve is a slice of the manifold at fixed openness φ, all with γ = 1.454. Storage-depleted states (φ < 1) sit below the classical Budyko curve, and storage-subsidized states (φ > 1) sit above it. The classical curve is the dashed line at φ = 1.

φ = 0.700.851.001.151.30
Interactive explorer

See how openness reshapes the partition

Move the sliders to change the openness coordinate φ and the landscape parameter γ. The curve, the water limit and the readouts update live from the published SQX equation.

φ < 1: storage-depleted (less water available to partition). φ > 1: storage-subsidized (for example glacier melt or groundwater inflow).
γ = 1.454 is the annual optimum calibrated at Goalpara on the Brahmaputra.
E / P–
Q / P–
0123450.00.40.81.21.6Aridity index x = PET / PEvaporative index E / Penergy limitwater limit E/P = φ
SQX curve (current φ, γ) Classical case, φ = 1 Energy and water limits
The manifold in 3D

One surface for every storage state

Here E/P is drawn over the aridity index x and the openness φ, using the γ value set above. The highlighted line is the curve for your current φ. The classical Budyko curve is the slice at φ = 1.

Drag to rotate. The surface turns slowly until you touch it.

low E/Phigh E/P
Empirical validation

Tested on the Brahmaputra at Goalpara

The framework was checked against observed discharge at the Goalpara gauging station, on the Brahmaputra for 2002–2021, using GRACE/GRACE-FO storage and TerraClimate climate data. It was then tested independently on seven basins in Myanmar. The results show both where SQX adds value and where it does not.

Monthly skill
0.650
Nash–Sutcliffe efficiency over the 63 admissible months, against −0.005 for the classical closed-basin form
Kling–Gupta
0.651
Kling–Gupta efficiency over the same admissible months (φ > 0)
Calibrated parameter
1.454
Annual optimum of the landscape parameter γ for the Goalpara catchment
Out-of-sample
0.476
Leave-one-out cross-validated NSE (out-of-sample check)

Where SQX adds value

  • Monthly scale, Goalpara: over the 63 admissible months, SQX reaches NSE 0.650, while the classical φ = 1 form reaches −0.005. The effective monthly γ is about 0.68.
  • Robustness: the result holds under leave-one-out cross-validation (NSE 0.476) and under GRACE Monte Carlo uncertainty propagation (2,000 realisations; mean NSE 0.603, 90% interval 0.16–0.90).
  • Diagnostics: the admissibility violations in 2010, 2013 and 2017 coincide with anomalously large positive ΔS. This is consistent with GRACE leakage from nearby glacierised Himalayan terrain, so SQX also helps flag data problems.

Where it does not, yet

  • Annual scale: at Goalpara, φ only ranges from 0.92 to 1.09, and SQX shows no predictive advantage (NSE −1.407 against −0.601 for the classical form).
  • Irrawaddy–Chindwin, Myanmar: across seven GRDC-gauged basins, SQX does not outperform the classical form in six. Diagnostic tests trace this to numerical conditioning of φ = 1 + ΔS/P under low dry-season precipitation, not to basin size, noise or a failure of the openness concept.
  • Domain: September–December at Goalpara and severely storage-depleted states (φ ≤ 0) lie outside the present equilibrium formulation.
Research programme

Where SQX stands, and where it is going

Each piece of work is labelled by its status, so peer-reviewed results stay separate from filed intellectual property and from work still in progress.

Peer-reviewed

Generalised 3D Budyko manifold

Analytical derivation, structural-stability proofs and Brahmaputra validation of the SQX surface.

Journal of Hydrology 679 (2026) 136406
Open preprint

Coupled water–energy extension

A dual Budyko view that links the hydrological partition to surface energy partitioning in open systems. It is being developed further.

Open-access preprint · ESS Open Archive
Patent pending

Quantum-enhanced Budyko modelling

Systems and methods that build on the SQX framework. Applicant: IIT Guwahati.

Indian patent application · complete specification filed
Doctoral thesis

Brahmaputra catchment behaviour

Climate and human drivers of hydrological change in the Brahmaputra basin, with SQX as the central methodological contribution.

Ph.D., Civil Engineering, IIT Guwahati
Peer-reviewed

Transferability across basins

An independent test on seven Irrawaddy–Chindwin basins in Myanmar defines an explicit boundary condition: the present φ = 1 + ΔS/P form is poorly conditioned when dry-season precipitation is low.

Reported in Journal of Hydrology 679 (2026) 136406
Future direction

Beyond ΔS/P: the ω–θ openness

Estimating openness as φ = 1 + ω − θ from independent land-surface information, so that it does not inherit the 1/P sensitivity. The plan also includes dynamic γ, routing-corrected φ and tests on more glacier-fed basins.

Priorities named in the JoH paper
Publications & preprints

Scholarly record

2026
Journal article

A generalized three-dimensional manifold for the Budyko hypothesis: analytical derivation and structural stability in open hydrological systems

Shakya, C. M., & Bhattacharjya, R. K.

Journal of Hydrology, 679, 136406. Elsevier.

2026
Preprint · open access

Coupled water–energy extension of the ShakyaQX framework for open hydrological systems

Shakya, C. M., & Bhattacharjya, R. K.

ESS Open Archive (ESSOAr) preprint. Free to read.

2025
Preprint

Bridging hydrological theory and machine learning: multi-scale Budyko–Fu diagnostics of runoff sensitivity in the Brahmaputra Basin

Shakya, C. M., & Bhattacharjya, R. K.

SSRN preprint.

Intellectual property

Patent application

Patent pending · India

Quantum-Enhanced Systems and Methods for Budyko-Based Hydrological Modeling

The application covers systems and methods that build on the SQX open-system partition framework. Only the published equation and results are described on this website. The claimed technical details are in the official patent record.

Applicant
Indian Institute of Technology Guwahati
Inventors
Chandra Mohan Shakya
Prof. Rajib Kumar Bhattacharjya
Application No.
202531055230 · Indian Patent Office, Kolkata
Status
Complete specification filed (2026)
Priority date: 7 June 2025
Milestones

How ShakyaQX developed

2020

Doctoral research begins

Ph.D. in Civil Engineering at IIT Guwahati, studying hydrological change in the Brahmaputra basin.

Groundwork

Limits of closed-basin theory

Budyko–Fu and machine-learning diagnostics across the basin show where classical assumptions break down.

Jun 2025

Provisional patent

The priority date is established for the ShakyaQX systems and methods.

2026

SQX manifold published

Journal of Hydrology 679, 136406: derivation, stability and validation.

2026

Complete specification filed

Indian patent application by IIT Guwahati, with both researchers as inventors.

Next

Water–energy & beyond

Coupled water–energy theory, forward-computed openness and quantum-enhanced calibration.

People

Research team

Subsurface Hydrology Laboratory, Department of Civil Engineering, Indian Institute of Technology Guwahati, Assam 781039, India.

Chandra Mohan Shakya

Ph.D. Scholar, Department of Civil Engineering

Indian Institute of Technology Guwahati

Developer of ShakyaQXFirst author · Inventor

Works on hydrological modelling and water resources engineering, combining theoretical hydrology with machine-learning analysis across the Brahmaputra and central Indian basins. He has about 3.5 years of teaching experience, including at NIT Raipur and NIT Hamirpur, and was previously a Project Associate/JRF at IIT Mandi on a SERB-funded project.

Interests: Budyko-type frameworks · open-system water balance · ML in hydrology · optimisation of water resources systems

chandra.mohan@iitg.ac.inLinkedIn ↗

Prof. Rajib Kumar Bhattacharjya

Professor, Department of Civil Engineering

Indian Institute of Technology Guwahati

Ph.D. SupervisorCo-author · Co-inventor

Prof. Bhattacharjya has been a Professor at IIT Guwahati since 2015 and received his Ph.D. in Civil Engineering from IIT Kanpur. He has been a visiting faculty member at Dalhousie University, Canada. His research covers optimisation methods, computational hydraulics, the impact of climate change, and surface water–groundwater interaction.

Faculty profile, IIT Guwahati ↗
How to cite

Citing ShakyaQX

If you use the SQX framework or equation, please cite the original publication:

Shakya, C. M., & Bhattacharjya, R. K. (2026). A generalized three-dimensional manifold for the Budyko hypothesis: Analytical derivation and structural stability in open hydrological systems. Journal of Hydrology, 679, 136406. https://doi.org/10.1016/j.jhydrol.2026.136406
@article{Shakya2026SQX,
  author  = {Shakya, Chandra Mohan and Bhattacharjya, Rajib Kumar},
  title   = {A generalized three-dimensional manifold for the
             Budyko hypothesis: analytical derivation and structural
             stability in open hydrological systems},
  journal = {Journal of Hydrology},
  volume  = {679},
  pages   = {136406},
  year    = {2026},
  doi     = {10.1016/j.jhydrol.2026.136406}
}
Acknowledgements

Credits

  • This work was carried out in the Department of Civil Engineering, Indian Institute of Technology Guwahati, under the supervision of Prof. Rajib Kumar Bhattacharjya.
  • The authors thank the members of the Doctoral Committee, chaired by Prof. Suresh A. Kartha, with members Prof. Arup Kumar Sarma and Prof. Archana M. Nair (Department of Civil Engineering, IIT Guwahati), for their guidance and review.
  • Observed discharge at Goalpara: Central Water Commission (CWC), Government of India. Myanmar discharge: Global Runoff Data Centre (GRDC).
  • Terrestrial water storage: GRACE / GRACE-FO CSR RL06 mascon solutions (Tapley et al., 2019; Save et al., 2020).
  • Precipitation and PET: TerraClimate v1.0 (Abatzoglou et al., 2018). Elevation: SRTM DEM (Farr et al., 2007).
  • The framework builds on the work of Budyko, Mezentsev, Choudhury, Yang and Fu on long-term water–energy balance.