The satellites in question belong to the Gravity Recovery and Climate Experiment, known as GRACE, a joint NASA and German Aerospace Center mission that operated from 2002 to 2017 and was later continued by its successor, GRACE Follow-On. Rather than photographing the surface, GRACE tracks the distance between two orbiting spacecraft with extraordinary precision. When the mass of water stored in soils, aquifers and surface reservoirs shifts beneath them, the resulting gravitational tug subtly alters that distance. By measuring these minute changes, scientists can calculate how much water a region gains or loses each month, a quantity called the change in terrestrial water storage. The technique has transformed hydrology, but it comes with a fundamental catch: the gravity signal is smeared across enormous areas, so a single GRACE data point represents a block of land roughly 150,000 to 200,000 square kilometers in size. For a basin where a farmer’s tube well in one district may be draining an aquifer while a neighboring district recharges after monsoon rains, that blur hides exactly the local detail that water managers need most.

In a study published in Environmental Monitoring and Assessment, Chethan Varadaganahalli Anandagowda, Bhaskar R. Nikam and Suresh Kannaujiya tackled this resolution problem by teaching statistical models to translate the coarse satellite signal into fine-grained maps. Their approach, known as downscaling, exploits the fact that water storage changes are closely linked to variables that satellites and models can observe at much finer scales, including rainfall, soil moisture, evapotranspiration, temperature and terrain. The researchers built monthly climatology-based regression models that relate these high-resolution environmental datasets to the GRACE measurements, then used the learned relationships to estimate water storage change on a grid of 0.25 degrees by 0.25 degrees, roughly 25 to 28 kilometers per side. That is a dramatic sharpening of the effective resolution, bringing the satellite picture down to a scale where individual districts and aquifer systems begin to separate.

The heart of the method was a comparison of several statistical and machine learning approaches, including multiple linear regression, partial least squares regression and random forest modeling. The random forest technique, an ensemble method that builds hundreds of decision trees and averages their predictions, emerged as the clear winner. It achieved a Pearson correlation coefficient of 0.78 against validation data, meaning it captured more than three quarters of the variance in the observed storage signal, with a root mean square error of 12.3 centimeters of equivalent water height. In practical terms, the model could reproduce the month-to-month ebb and flow of the basin’s water storage with an accuracy that makes the downscaled maps genuinely useful rather than merely illustrative. The team validated their results using k-fold cross validation, a technique that repeatedly trains the model on most of the data and tests it on the withheld remainder, guarding against the model simply memorizing patterns rather than learning real physical relationships.

With the sharpened storage maps in hand, the researchers turned to the question that matters most for the Ganga Basin: what is happening to its groundwater specifically. GRACE measures all water stored on and below the land surface, so isolating the groundwater component requires subtracting the contributions of soil moisture, surface water and canopy storage. The team used products from the Global Land Data Assimilation System, a NASA modeling framework that simulates land surface water fluxes, together with satellite soil moisture records from the European Space Agency’s Climate Change Initiative, to strip away these shallower reservoirs. What remained was an estimate of groundwater storage change, which they then compared against observations from monitoring wells across the basin, converted to storage terms using specific yield values for the aquifer material.

The agreement between the space-based and ground-based estimates was striking, with both showing a consistent declining trend across the basin. But the fine-resolution data revealed something the coarse GRACE footprint had smoothed away: the decline is not uniform, and in some places it is far steeper than previous studies had suggested. Over the Delhi National Capital Region and Western Uttar Pradesh, regions where explosive population growth and intensive irrigation place enormous pressure on the aquifers, the downscaled estimates showed groundwater storage falling at rates of 4.69 to 7.16 centimeters of equivalent water height per year. Those figures are more severe than the trends reported in earlier basin-wide analyses, and they pinpoint the western Ganga plains as the epicenter of the basin’s groundwater crisis. The finding matters because it demonstrates that coarse satellite data, by averaging heavily depleted zones with less stressed ones, can systematically understate the intensity of depletion in the worst-affected areas.

The study went beyond storage trends to tackle drought, using multiple drought indices computed from the water storage record to identify major hydrometeorological drought events across the basin. Drought indices condense complex hydrological information into standardized measures of severity, allowing events to be compared across time and space. The analysis identified several significant drought episodes during the satellite era, but one stood out above the rest: the period from August 2015 to October 2017, which the researchers identified as the most intense drought in the record, persisting for a remarkable 27 consecutive months. That two-and-a-quarter-year event, which coincided with consecutive weak monsoon seasons, placed sustained stress on both surface reservoirs and the aquifers that farmers turned to when canals and rivers ran low. Because the downscaled data resolves the drought’s footprint at district scale, it shows which parts of the basin bore the brunt of that prolonged deficit rather than blurring the event into a basin-wide average.

The implications extend well beyond academic hydrology. India’s Central Ground Water Board assesses the country’s groundwater resources annually, but its estimates depend on the density and quality of well observations, which vary widely across states. Satellite-derived storage change offers an independent check on those assessments, and the downscaling approach demonstrated here brings that check down to a scale where it can inform district-level management decisions, from regulating well permits to planning managed aquifer recharge. The researchers note that their findings enhance understanding of local-scale groundwater abstractions and can contribute to sustainable water resource management in a basin where climate change is expected to make the monsoon more erratic even as demand continues to grow. Every centimeter of equivalent water height in these maps represents billions of cubic meters, so a decline of seven centimeters per year in a heavily farmed region translates into a withdrawal of water that took centuries or millennia to accumulate.

There are also lessons here for how the global scientific community monitors water in data-sparse regions. The Ganga Basin is hardly unique: aquifers from the North China Plain to California’s Central Valley to the Middle East are being drawn down faster than they refill, and in many of these places the observational infrastructure is thin. The success of the random forest approach, which outperformed classical regression methods by capturing nonlinear relationships between climate variables and storage change, suggests a transferable recipe. The datasets the team relied on, including GRACE mass grids processed at the University of Texas Center for Space Research, CHIRPS rainfall estimates, GLDAS model outputs, ESA soil moisture records and Indian well observations, are all freely available, meaning the workflow could be replicated for other stressed basins at little cost. As the GRACE Follow-On mission continues to deliver gravity data and machine learning tools grow more capable, the gap between what satellites can see and what water managers need to know keeps narrowing. For the Ganga Basin, that narrowing gap has just exposed an uncomfortable truth: the water is leaving faster, and in more places, than the blur of the old maps ever let on.

Subject of Research: Downscaling GRACE satellite gravity data with machine learning to assess groundwater depletion and drought in the Ganga Basin

Article Title: Downscaling GRACE storage change to analyze groundwater dynamics and drought evaluation in the Ganga Basin

Article References: Anandagowda, C. V., Nikam, B. R., & Kannaujiya, S. (2026). Downscaling GRACE storage change to analyze groundwater dynamics and drought evaluation in the Ganga Basin. Environmental Monitoring and Assessment, 198(11), Article 1146. https://doi.org/10.1007/s10661-026-16002-9

Image Credits: AI Generated

DOI: 10.1007/s10661-026-16002-9

Keywords: GRACE, groundwater, Ganga Basin, downscaling, random forest, machine learning, drought, terrestrial water storage, remote sensing, Delhi-NCR, Uttar Pradesh, water resources

The stakes could hardly be higher. Inorganic arsenate, or As(V), is a class-one carcinogen, and prolonged exposure through contaminated drinking water and food chains is linked to cancers and a range of other serious illnesses. Across India, the sand, silt and clay deposits laid down by the Ganga and its tributaries are a major reservoir of the element, which can shift from inert mineral phases into soil solution and groundwater whenever chemistry and redox conditions favor mobilization. Whether arsenic stays locked in place or leaks into aquifers depends largely on adsorption, the process by which dissolved ions attach themselves to reactive mineral surfaces. Soils that adsorb arsenate strongly act as natural buffers; soils that adsorb it weakly let the poison travel.

Most laboratory research on arsenic removal has focused on pure mineral phases, especially iron oxides, which are celebrated sorbents of arsenate. But natural soils are messy composites of clay minerals, oxides, organic matter and variable surface chemistry, and the team wanted to know how the whole package performs. They collected three alluvial soils from Varanasi spanning a textural gradient: a sandy loam, a sandy clay loam and a clay loam. Each was characterized in detail using a battery of techniques, including nitrogen gas adsorption to measure surface area, scanning electron microscopy with energy-dispersive X-ray spectroscopy to image surface morphology, Fourier-transform infrared spectroscopy to identify reactive functional groups, and X-ray diffraction to pin down the mineral assemblage.

The characterization told a coherent story. The clay loam, the finest-textured of the three, had the largest BET surface area at 11.08 square meters per gram, the highest cation exchange capacity at 12.4 centimoles of charge per kilogram, and the most organic carbon at 0.42 percent. X-ray diffraction revealed quartz, feldspar and aluminosilicate phases in all three soils, along with clay minerals such as illite and kaolinite. Infrared spectroscopy detected broad hydroxyl bands and metal-oxygen vibrations, the very groups known to participate in ligand-exchange reactions with arsenate. Electron microscopy of the clay loam showed a heterogeneous, highly aggregated mineral matrix with platy particles and abundant inter-particle voids, exactly the kind of architecture that exposes reactive surfaces to dissolved contaminants.

With the soils characterized, the team ran batch adsorption experiments at 25 degrees Celsius and neutral pH. Kinetics experiments showed rapid initial uptake, with the clay loam adsorbing 0.296 milligrams of arsenate per gram within the first thirty minutes, rising to 0.494 milligrams per gram by twenty-four hours, after which adsorption plateaued. The researchers fitted the data to pseudo-first-order, pseudo-second-order, intra-particle diffusion and Elovich models. The pseudo-second-order model performed best, pointing to chemisorption, meaning arsenate forms actual chemical bonds with surface functional groups rather than merely sticking by weak physical forces. The strong performance of the Elovich model, meanwhile, indicated that the adsorption sites are energetically heterogeneous, a signature of natural soils with mixed mineralogy.

The pH dependence of adsorption proved dramatic. All three soils grabbed the most arsenate in the acidic range between pH 4 and pH 6, and adsorption fell sharply as conditions turned alkaline. At pH 7, the sandy loam, sandy clay loam and clay loam adsorbed 0.350, 0.420 and 0.490 milligrams per gram respectively, but at pH 10 those figures collapsed to just 0.100, 0.178 and 0.248. The explanation lies in surface charge. Under acidic conditions, protonated soil surfaces carry positive charge that electrostatically attracts the negatively charged arsenate species H2AsO4− and HAsO4²−. As pH rises, hydroxide ions compete for binding sites and clay surfaces become increasingly negative, repelling arsenate and driving it back into solution. This is bad news for alkaline soils and groundwater systems, where arsenate mobility is inherently higher.

Equilibrium isotherm experiments, spanning arsenate concentrations from 2.5 to 50 milligrams per liter, quantified the capacity of each soil. The Langmuir model estimated maximum monolayer adsorption capacities of 0.531, 0.612 and 0.662 milligrams per gram for the sandy loam, sandy clay loam and clay loam respectively, a clear texture-driven hierarchy. Yet the Langmuir model, which assumes uniform binding sites, actually fit the data relatively poorly. The Freundlich and Sips models, which explicitly account for surface heterogeneity, performed better, confirming that arsenate retention in these soils is distributed across a mosaic of sites with different affinities: clay mineral edges, hydroxyl-bearing surfaces, and reactive Al-OH groups. The Temkin model further revealed a higher heat of adsorption in the clay loam, signaling stronger overall interactions.

The most consequential part of the study came when the researchers introduced competition. Groundwater is never a pure arsenate solution; it carries phosphate from fertilizers, silicate from mineral weathering, and sulphate from a host of natural and anthropogenic sources. All three are oxyanions that can occupy the same inner-sphere binding sites arsenate uses. When the team added 10, 20 and 30 milligrams per liter of each competitor alongside 20 milligrams per liter of arsenate, adsorption dropped significantly in every soil. The competitive strength followed a clear order: phosphate was the fiercest rival, silicate came second, and sulphate was the weakest. At 30 milligrams per liter, phosphate slashed arsenate adsorption to 0.092, 0.082 and 0.130 milligrams per gram in the three soils, cutting retention by more than half in the sandy loam and sandy clay loam.

Phosphate’s dominance is chemically logical. Arsenate and phosphate are structural near-twins, with almost identical acid dissociation constants, so they compete for precisely the same hydroxyl-bearing mineral sites through ligand exchange and inner-sphere complexation. Silicate binds similarly and can even polymerize on surfaces, gradually clogging adsorption sites, while sulphate’s weaker affinity leaves more room for arsenate. Intriguingly, the inhibitory effect of every competitor weakened progressively from the sandy loam to the clay loam. The clay-rich soil, with its larger and more diverse pool of reactive sites, simply had more binding capacity to spare, buffering it against competition. The authors conclude that competitive adsorption is governed not only by the chemical affinity of each anion for mineral surfaces but also by the abundance and heterogeneity of reactive sites within the soil itself.

The implications ripple outward in two directions. For arsenic mitigation, the study suggests that clay-rich alluvial soils are naturally better at immobilizing arsenate, and that soil texture and mineralogy deserve more attention when mapping contamination risk across the Indo-Gangetic plain. For agriculture, the phosphate result is a caution: heavy phosphorus fertilization on arsenic-affected soils could displace arsenate from grain surfaces and push it into soil solution and crops, a trade-off farmers and land managers rarely see. The authors are careful to note that their batch experiments, conducted at fixed pH and single-anion conditions, simplify the far messier reality of field soils, where redox fluctuations, dissolved organic matter, seasonal water flow and multiple competing ions interact. Still, by connecting spectroscopic mineral characterization to quantitative adsorption behavior under realistic competitive pressure, the work delivers a practical message: the ground beneath arsenic-affected regions is not just a passive contaminant source but an active, texture-dependent filter, and understanding its chemistry may be one of the cheapest tools available for protecting water and food supplies.

Subject of Research: Competitive adsorption of arsenate on alluvial soils under phosphate, silicate and sulphate conditions

Article Title: Assessment of adsorption potential of alluvial soils for arsenic removal under competing anion conditions

Article References: Chattopadhyay, A., Singh, A. P., Rakshit, A., Barman, A., Chourasia, A., & Mondal, T. (2026). Assessment of adsorption potential of alluvial soils for arsenic removal under competing anion conditions. Discover Soil, 3(1), Article 177. https://doi.org/10.1007/s44378-026-00332-8

Image Credits: AI Generated

DOI: 10.1007/s44378-026-00332-8

Keywords: arsenic, arsenate, alluvial soil, adsorption, phosphate, silicate, sulphate, clay minerals, groundwater contamination, soil chemistry, Ganga basin, chemisorption

The stakes could hardly be higher. The glaciers of the Indian Central Himalaya feed the headwaters of the Indus, Ganga and Brahmaputra river systems, which sustain hundreds of millions of people across South Asia. Glacier meltwater acts as a seasonal buffer for river flow, releasing stored ice during dry summer months when downstream demand peaks. As the ice reservoirs shrink, that buffering capacity erodes, threatening both water security and agricultural stability. At the same time, deglaciation is spawning new hazards: the review documents the proliferation of glacial lakes, including roughly 47 lakes that have formed in front of the retreating Milam Glacier, raising the specter of glacial lake outburst floods, debris flows and ice avalanches in densely populated valleys below.

To build their synthesis, the researchers applied the PRISMA framework, a structured protocol for systematic reviews. Starting from 3,412 records identified in the ScienceDirect and Springer Nature databases, they screened out duplicates, non-English publications and studies lacking quantitative recession data, ultimately retaining 50 peer-reviewed studies published between 1990 and 2024. The team standardized the extracted measurements — frontal retreat rates in meters per year, mass balance in meters of water equivalent, and equilibrium line altitudes in meters above sea level — so that figures from different studies, sensors and observation periods could be compared directly. They also catalogued the sources of uncertainty, from differences in satellite spatial resolution and georeferencing errors to inconsistencies in how glacier boundaries were delineated in the field versus from orbit.

The headline numbers are striking. Across 23 glaciers with compiled retreat records, the mean frontal retreat rate is 23.23 meters per year, but the spread is enormous: from roughly 3.7 meters per year at the relatively stable Bhagirathi Kharak Glacier to a staggering 88 meters per year at Khatling Glacier, which has lost more than 4,000 meters of length since 1965. Approximately 56 percent of the glaciers examined show moderate retreat of 10 to 30 meters per year, while 26 percent are retreating faster than 30 meters per year. Lake-terminating glaciers retreat significantly faster, averaging about 34.9 meters per year, compared with 22.5 meters per year for land-terminating ones — clear evidence that calving and water-ice thermal erosion accelerate ice loss where glaciers end in proglacial lakes.

Individual glacier case studies illustrate the complexity. The Milam Glacier in the Kumaon sector receded by an average of 480 meters in a single decade between 2004 and 2014, its snout retreating so rapidly that the frontal area shrank and the glacier became detached from a tributary. The Kafni Glacier lost 1,057 meters of length between 1972 and 2018, an average of nearly 23 meters per year. The Gangotri Glacier, one of the most closely watched ice masses in the region, shows surface velocities ranging from about 25 to 29 meters per year near the snout to as much as 70 meters per year in its interior, indicating that the lower tongue is being drawn downslope and consumed faster than the upper reaches. Meanwhile, the Satopanth Glacier has retreated at a comparatively sedate 2.9 to 8.5 meters per year, depending on the period and method of measurement.

Why such divergence under a shared regional climate? The answer, the review argues, lies in what glaciologists call morphometric controls — the physical characteristics of each glacier itself. Supraglacial debris cover emerges as the single most important protective factor: a thick blanket of rock and sediment insulates the ice beneath, slowing surface melt. The Satopanth Glacier, with roughly 32 percent of its surface mantled in debris and a favorable northeast orientation, shows a mass balance deficit of only about minus 0.28 meters of water equivalent per year, whereas the steep, nearly debris-free Khatling Glacier loses about minus 0.62 meters per year. Quantitatively, recent geodetic work suggests that a 10 percent increase in debris cover can buffer mass loss by approximately 0.36 meters of water equivalent annually, while a 10 percent increase in slope corresponds to roughly 0.86 meters of additional loss.

Glacier geometry matters too. The review highlights a ‘shape index’ — the ratio of glacier length to width — as a predictor of resilience. Elongated valley glaciers such as Gangotri, with a shape index of 4.8, maintain stronger ice flux from accumulation zones to ablation zones and greater structural integrity, while small, roughly circular glaciers like Chorabari, with a shape index of 1.9, have high perimeter-to-area ratios and are far more vulnerable to rapid thermal breakdown. Aspect plays a role as well: south-facing glaciers, which receive more solar radiation, have shrunk by about 18 percent, considerably more than their north-facing counterparts. Local topography, morphology and ice velocity together explain up to 46 percent of the variation in glacier response — a reminder that climate sets the direction of change but terrain sets its pace.

The climatic backdrop itself has shifted. Analysis of meteorological data from 1980 to 2023, drawing on CHIRPS precipitation and reanalysis temperature products, reveals a transition around the year 2000. The preceding two decades were characterized by a ‘warming and drying’ trend, with precipitation declining at about 21 millimeters per year and temperatures rising at 0.047 degrees Celsius per year. After 2000, precipitation increased at roughly 32 millimeters per year while temperature rise nearly flattened. This shift toward a warmer-but-wetter regime helps explain why retreat rates of some eastern glaciers, such as Pindari and Milam, slowed somewhat after 2010, as increased snowfall partially buffered mass losses. Yet the same humidity and meltwater ponding raises the risk of glacial lake expansion and outburst floods, and the snowline has climbed roughly 100 meters since 1994, fragmenting large stable ice bodies into smaller, more fragile pieces — the glacier count in the Alaknanda Basin has jumped about 26 percent through this breakup process.

The mass balance ledger is grimly consistent. In the Upper Ganges basin, eight glaciers showed cumulative mass losses between minus 0.15 and minus 5.04 gigatonnes between 2000 and 2014, a weighted average of minus 0.61 meters of water equivalent per year, equivalent to 16 gigatonnes of lost freshwater storage. Thinning rates in the Garhwal Himalaya quadrupled, from 0.07 meters per year between 2000 and 2015 to 0.31 meters per year between 2015 and 2020, and the equilibrium line altitude at Nanda Devi rose by nearly 500 meters between 1980 and 2017. The review’s authors warn that cumulative glacier length in the region could shrink by around 40 percent by 2050, a prospect with profound implications for the roughly 800 million people dependent on the Ganga basin. They call for expanded mass balance monitoring, better estimates of ice thickness and volume, and tighter integration of field surveys with optical and radar satellite data, unmanned aerial vehicles and digital elevation models. Only about eleven of the roughly 10,000 glaciers in the Indian Himalaya have been studied in detail for mass balance — a gap that, as this synthesis makes clear, the region can no longer afford.

Subject of Research: Multi-decadal glacier recession, mass balance and equilibrium line altitude changes in the Indian Central Himalaya

Article Title: Time dependent assessment of glacier recession in the Indian Central Himalaya

Article References: Kumari, J., Kothyari, G. C., Patidar, A. K., & Mehta, M. (2026). Time dependent assessment of glacier recession in the Indian Central Himalaya. Discover Geoscience, 4(1), Article 357. https://doi.org/10.1007/s44288-026-00720-1

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00720-1

Keywords: Himalayan glaciers, glacier retreat, mass balance, equilibrium line altitude, Uttarakhand, climate change, remote sensing, supraglacial debris, glacial lakes, Ganga basin, cryosphere, satellite imagery