Ribeiro’s scientific foundation was laid in geostatistics, the branch of applied mathematics developed originally for mining estimation and later adopted across the environmental sciences. His 1992 doctoral thesis at Instituto Superior Técnico, written within the framework of mining engineering, was a substantial work of more than four hundred pages devoted to the geostatistical characterisation of hydrogeological systems, including contributions to structural model inference under nonstationary geostatistics. This was a technically demanding area: nonstationary models allow the statistical structure of a field, such as a groundwater level surface, to vary across space rather than assuming uniform behaviour, which is far closer to how real aquifers behave. As early as 1989, Ribeiro had co-authored work with F. Muge on a geostatistical approach to modelling a piezometric field, presented in the influential Geostatistics book series, demonstrating how kriging-based methods could be brought to bear on the spatial structure of groundwater heads.

That quantitative grounding soon produced practical tools that spread well beyond Portugal. Perhaps the best known is the Susceptibility Index, a method Ribeiro developed to assess the vulnerability of aquifers to diffuse agricultural pollution, first applied in the Alentejo region of southern Portugal and published in the proceedings of the 2001 Future Groundwater Resources at Risk conference. Vulnerability mapping itself has a distinguished lineage, traceable to the French work of Margat in 1968 and Albinet and Margat in 1970, and later systematised in the American DRASTIC scheme of 1987. Ribeiro’s contribution was to adapt the concept to the specific dynamics of agricultural contamination, where nitrate leaching from fertilisers and irrigation return flows threatens shallow aquifers. The method was subsequently applied in strikingly different settings, including the Daule aquifer in Ecuador, where Ribeiro and colleagues published a vulnerability assessment in the Science of the Total Environment in 2017, showing how a Portuguese-designed index could be transferred to a tropical coastal aquifer facing salinisation and agrochemical pressure.

The vulnerability work was inseparable from a broader research programme on groundwater quality and monitoring, much of it conducted with Stigter and Manuel João Carvalho Dill. In a pair of papers published in Hydrogeology Journal and the Journal of Hydrology in 2006, the team evaluated intrinsic and specific vulnerability methods against observed salinisation and nitrate contamination in two agricultural regions of southern Portugal, and then developed a groundwater quality index as a communication tool for agro-environmental policy. The point was not merely academic: indices translate complex multivariate chemistry into a form that regulators, farmers and the public can act upon. Later work pushed the statistical sophistication further, with factorial regression models built to explain and predict nitrate concentrations under agricultural land, published in the Journal of Hydrology in 2008, and a 2011 study in Environmental Science & Technology that examined the efficiency of monitoring programmes for nitrate-contaminated groundwater, asking hard questions about whether existing networks actually detect what they are supposed to detect.

Monitoring network design became a speciality in its own right. Together with Luís Nunes and Maria da Conceição Cunha, Ribeiro published a series of papers in the mid-2000s on optimising groundwater monitoring networks through redundancy reduction and on optimal space-time coverage under exploration cost constraints, appearing in the Journal of Water Resources Planning and Management and Environmental Monitoring and Assessment. The underlying problem is one that every environmental agency faces: sensors and sampling campaigns are expensive, yet removing the wrong wells destroys the information content of the entire network. By combining variance-reduction criteria from geostatistics with optimisation algorithms, the team provided a quantitative basis for deciding where each measurement earns its keep. A 2007 paper in Computer-Aided Civil and Infrastructure Engineering compared variance-reduction and space-filling design approaches, clarifying the trade-offs between statistical efficiency and spatial coverage.

Ribeiro’s geostatistical toolkit also reached back toward its origins in reservoir characterisation. A 1990 paper in Mathematical Geology, co-authored with colleagues including Amílcar Soares, applied geostatistical and multivariate data analysis techniques to improving reservoir description, connecting hydrogeology with petroleum geoscience. Multivariate statistics became a recurring thread throughout his career. In Colombia, working with Teresita Betancur-Vargas and her students at the Universidad de Antioquia, principal component analysis was used to untangle the hydrochemistry of the Urabá aquifer, published in the Journal of Geochemical Exploration in 2013, and to identify regional flow systems in the Gulf of Urabá. Similar approaches informed studies of pesticide exposure in the shallow groundwater of the Tagus vulnerable zone in Portugal, published in Environmental Science and Pollution Research in 2012, and of natural background levels assessment combined with indicator kriging for groundwater quality interpretation, published with Italian colleagues in the Science of the Total Environment in 2016.

Climate change entered his research agenda early and stayed there. With Luís Veiga da Cunha, Ramiro Oliveira and João Nascimento, Ribeiro produced some of the first assessments of climate change impacts on groundwater resources in mainland Portugal, published in Portuguese water-sector literature in the mid-2000s and in IAHS proceedings in 2007. The work extended to the Madeira archipelago through the CLIMAAT II project and to the Mediterranean, where Stigter led a comparative assessment of climate change impacts on three coastal aquifers published in Regional Environmental Change in 2014. Further afield, Ribeiro contributed to evaluating piezometric trends using the seasonal Kendall test, a nonparametric method robust to seasonality, in the alluvial aquifers of the Elqui river basin in north-central Chile, published in Hydrological Sciences Journal in 2015, and to analysing how transmissivity uncertainty affects the reliability of models of the Northwestern Sahara Aquifer System, published in the Journal of African Earth Sciences in 2017 with Mounira Zammouri.

One of the most forward-looking strands of his later career concerned groundwater-dependent ecosystems and the living world hidden inside aquifers. Ribeiro co-edited the volume Groundwater and Ecosystems in the International Association of Hydrogeologists’ Selected Papers in Hydrogeology series, published by CRC Press in 2013, and worked with colleagues on groundwater ecosystems and bio-indicators. With Mafalda Shapouri, Luís Cancela da Fonseca and others, he studied the variation of stygofauna, the specialised animals inhabiting groundwater, along a gradient of salinisation risk in a coastal Mediterranean aquifer, published in Hydrology Research in 2016. Estuarine biodiversity was examined as an indicator of groundwater discharge in Estuarine, Coastal and Shelf Science in 2012. These interests anticipated by years the growing international recognition of groundwater as an ecosystem in its own right, a view crystallised in a 2023 Global Change Biology review describing groundwater as a hidden global keystone ecosystem. Ribeiro also contributed methodological work on identifying groundwater-dependent ecosystems in Portugal for the implementation of the European Water Framework Directive, presented at the 44th IAH Congress in Dubrovnik in 2017.

What truly set Ribeiro apart, and what gives the tribute its unusual warmth, is the cultural dimension of his legacy. He was a lifelong student of the relationship between water and music, tracing how composers from Handel to Toru Takemitsu drew inspiration from water in its many forms. He presented a radio programme on Antena 2 in Portugal devoted to Shostakovich, created the Musicágua project exploring water as a source of inspiration in musical creation, and in 2021 published a chapter on music inspired by groundwater and other components of the hydrological cycle in a Springer volume on geoethics and groundwater management. In the same volume he revisited ancestral groundwater techniques, arguing that traditional collection and distribution systems deserve renewed attention as nature-based solutions for modern water management, an argument he illustrated with a study of a public groundwater supply in the São Pedro do Sul municipality based on an ancestral capture system.

The tribute’s authors frame this breadth not as dilettantism but as a coherent humanistic vision: the conviction that groundwater science serves society best when it is quantitative yet humble, technically rigorous yet culturally literate. His career, spanning Portugal, Spain, Morocco, Chile, Colombia, Ecuador and beyond, and his mentorship of students and collaborators across three continents, built a community as much as a publication record. In an era when hydrogeology faces intensifying pressures from climate change, agricultural intensification and ecosystem degradation, the example of a scientist who could design a monitoring network, map an aquifer’s vulnerability, champion the unseen life within groundwater and then explain why a Shostakovich symphony and a spring share the same spirit, stands as a reminder of what the discipline, at its best, can be.

Subject of Research: The interdisciplinary career and legacy of hydrogeologist Luís Ribeiro in groundwater geostatistics, vulnerability assessment and water culture

Article Title: Luís Ribeiro: A visionary mind bridging groundwater science, culture and music through an interdisciplinary and humanistic legacy

Article References: Stigter, T. Y., Nascimento, J., de Melo, M. T. C., & Betancur-Vargas, T. (2026). Luís Ribeiro: A visionary mind bridging groundwater science, culture and music through an interdisciplinary and humanistic legacy. Hydrogeology Journal. https://doi.org/10.1007/s10040-026-03137-9

Image Credits: AI Generated

DOI: 10.1007/s10040-026-03137-9

Keywords: Luís Ribeiro, hydrogeology, geostatistics, groundwater vulnerability, Susceptibility Index, groundwater-dependent ecosystems, stygofauna, climate change, monitoring network design, nitrate contamination, water and music, Portugal

The logic behind a baseline study is deceptively simple but scientifically crucial. Every patch of ground on Earth emits ionizing radiation, driven by naturally occurring radionuclides such as uranium-238, thorium-232, and potassium-40 in soils and building materials, along with a contribution from cosmic rays. When a facility that handles radioactive material begins operating, regulators and the public will inevitably ask whether the environment has changed. Without a rigorous, statistically defensible picture of conditions before operations start, it becomes nearly impossible to distinguish pre-existing geological variability from any future operational influence. The new study addresses that problem head-on by characterizing the ambient gamma-radiation field at the South Cikarang site with unusual methodological care.

The researchers designated nineteen monitoring locations across terrain that is anything but uniform, mixing engineered fill, natural clay, and vegetated ground. Because one designated point, labeled Point B, was physically inaccessible, direct measurements were ultimately obtained at eighteen locations. At each accessible point, the team measured the ambient dose equivalent rate, written H*(10), which is the standard operational quantity for estimating external exposure to penetrating radiation. Crucially, measurements were taken at two heights: at the ground surface and at one meter above ground, the height at which a standing human body would receive most of its external dose. Each detector-and-height configuration was repeated ten times, allowing the team to quantify measurement uncertainty rather than relying on single spot readings.

What makes the study technically interesting is its dual-detector strategy. The team used two fundamentally different radiation detection technologies in parallel: a Geiger-Müller tube, a gas-filled detector prized for ruggedness and simplicity, and a sodium iodide scintillation detector, NaI(Tl), which converts gamma-ray energy into flashes of light and offers higher sensitivity. The two detector types respond differently to the energy spectrum of environmental gamma rays, so agreement between them is a powerful internal check. Across all detector systems and both measurement heights, the ambient dose equivalent rates ranged from 0.030 to 0.120 microsieverts per hour, values comfortably within the range typical of natural background worldwide. The median reduction in dose rate from the surface to one meter height was 11.8 percent for the GM detector and 14.2 percent for the NaI(Tl) system, a pattern consistent with gamma rays originating primarily from radionuclides in the soil itself.

The statistical agreement between the two instruments was striking. Pearson correlation coefficients reached 0.958 at ground level and 0.919 at one meter, indicating that despite their different physical principles, the two detector systems told essentially the same story across the site. In radiation metrology, this kind of cross-validation matters enormously. A single detector can drift, suffer from energy-response artifacts, or be fooled by site-specific conditions. When two independent technologies converge on the same spatial pattern with nearly identical rankings of hot and cold spots, confidence in the resulting baseline rises dramatically. The authors also subjected their data to formal uncertainty analysis, following international guidance from the International Atomic Energy Agency on quantifying uncertainty in nuclear analytical measurements, which strengthens the study’s defensibility in a regulatory context.

Perhaps the most visually compelling part of the work is its use of geostatistics. Rather than treating the eighteen measurement points as isolated numbers, the team applied Ordinary Kriging, a spatial interpolation technique borrowed from mining geology and soil science. Kriging does not simply draw smooth contours between points; it models the underlying spatial structure of the data, using a variogram to describe how quickly radiation values become uncorrelated with distance. The analysis identified meaningful spatial dependence over a range of approximately 100 meters, meaning that measurements taken closer together than that distance carry information about one another, while points farther apart are effectively independent. That length scale is itself a finding: it tells future monitoring programs how densely they need to sample to capture the site’s true variability.

The geostatistical mapping also revealed a clear physical explanation for the radiation pattern. Higher dose rates clustered over areas of engineered sand and crushed-stone fill, materials that typically contain mineral grains enriched in naturally occurring radionuclides, while lower values were generally associated with natural clay-dominated and vegetated ground. This association between construction fill and elevated gamma readings echoes findings from other parts of Indonesia, including previous work by overlapping research teams in the Mamuju region of West Sulawesi, which is known for unusually high natural radiation. In a rapidly industrializing landscape like Cikarang, where ground is routinely reshaped and imported fill is common, understanding that human-modified ground can shift the local radiation field is essential context for any future facility monitoring.

To translate the measurements into human-health terms, the team performed screening-level dose and risk assessments. Annual external dose estimates for a hypothetical person at the site ranged from 0.052 to 0.172 millisieverts per year, and excess lifetime cancer risk estimates ranged from 1.76 times ten to the minus four to 5.86 times ten to the minus four. These figures should be read as conservative screening values rather than predictions of harm; they sit well within the range of doses people receive from natural background radiation everywhere on the planet, which averages a few millisieverts per year globally according to United Nations Scientific Committee assessments. The point of the exercise is not alarm but calibration: by quantifying the risk landscape before the cyclotron is built, the researchers have created a reference line against which any future change can be judged.

Medical cyclotrons occupy a distinctive niche in the radiological world. Unlike nuclear power plants, they accelerate charged particles to bombard targets and produce short-lived isotopes such as fluorine-18 and technetium precursors used in diagnostic imaging. Their routine emissions are modest, but neutron activation, activated components, and gaseous releases during target irradiation mean that regulators, following standards set by Indonesia’s nuclear regulator BAPETEN and international guidance from the IAEA and the International Commission on Radiological Protection, still require careful environmental surveillance. A credible pre-operational baseline is the foundation of that surveillance, and the South Cikarang study demonstrates how to build one with limited resources: replicate measurements, dual detectors, height-resolved sampling, and spatial statistics.

The broader significance of the study lies in its transferability. The authors describe their approach as a framework integrating measurement validation, uncertainty analysis, and spatial modeling for environmental surveillance around planned radiological facilities, and the recipe is deliberately generic. Any country planning a cyclotron, a radioisotope production plant, or even a research reactor could adapt the same workflow, adjusting the sampling density to the spatial correlation scale revealed by its own variogram. As nuclear medicine expands across Southeast Asia and the global south, demand for radiopharmaceuticals is growing faster than the infrastructure to produce them, and with that growth comes a parallel need for environmental science that earns public trust. This study, supported by BRIN and an ITB research grant, with site access provided by the state pharmaceutical company PT Bio Farma, shows that meticulous baseline work can be done before the first beam is ever switched on, ensuring that when questions arise about a facility’s environmental footprint, the answer will rest on data rather than doubt.

Subject of Research: Pre-operational environmental gamma-radiation baseline characterization at a proposed medical cyclotron facility in Indonesia

Article Title: Environmental radiological baseline characterization of a proposed medical cyclotron facility in Indonesia using dual-detector measurements and geostatistical analysis

Article References: Purnama, D. S., Muliawan, D., Santoso, M., Fitriana, R., Nugraha, E. D., Abdullah, A. R. A., & Permana, S. (2026). Environmental radiological baseline characterization of a proposed medical cyclotron facility in Indonesia using dual-detector measurements and geostatistical analysis. Environmental Geochemistry and Health, 48(15), Article 617. https://doi.org/10.1007/s10653-026-03505-0

Image Credits: AI Generated

DOI: 10.1007/s10653-026-03505-0

Keywords: medical cyclotron, radiological baseline, ambient dose equivalent, Geiger-Müller detector, NaI(Tl) scintillator, Ordinary Kriging, geostatistics, environmental monitoring, radiation protection, nuclear medicine, Indonesia, West Java

The stakes are high. Arsenic is a naturally occurring metalloid that becomes toxic to humans at relatively low exposures, and the World Health Organization recognizes it as a major public health concern. In Guizhou, the problem has a notorious history: chronic arsenic poisoning has previously been documented in villages where residents burned coal with exceptionally high arsenic content indoors. When arsenic sits in the topsoil of farmland, it can enter the food chain through crops, be inhaled as dust, or leach into the groundwater that threads through karst aquifers. Mapping its spatial distribution is therefore not an academic exercise but a prerequisite for protecting food safety in one of China’s most geologically distinctive agricultural regions.

The research team, led by Zhizhuo Liu and Lang Zhang with colleagues from Beijing Normal University, the Institute of Geophysical and Geochemical Exploration, and Tianjin Chengjian University, collected agricultural topsoil samples across Zhijin County and measured arsenic concentrations in the laboratory. The results revealed a striking range: values spanned from as low as 1.79 milligrams per kilogram of soil to more than 40 milligrams per kilogram, with a mean of 17.7 milligrams per kilogram. That spread matters, because it means some fields are relatively clean while others approach or exceed thresholds of concern, and the difference between them is not random. Understanding what drives that heterogeneity was the central question of the study.

Before any modeling could begin, the researchers had to confront a statistical challenge common in soil geochemistry: arsenic concentration data are typically skewed, with a long tail of high values that can distort conventional analyses. The team applied a Yeo-Johnson transformation, a flexible mathematical procedure that reshapes the distribution to reduce skewness. The transformation helped, but formal tests still rejected normality of the transformed data, with a p-value of 0.003. This detail is more than statistical housekeeping. It signals that arsenic in these soils is genuinely patchy and structured by underlying processes, rather than varying smoothly and randomly, which shaped the choice of methods that followed.

To characterize the spatial structure, the researchers turned to semivariograms, the workhorse tool of geostatistics. A semivariogram describes how similar soil values are as a function of the distance between sampling points, and fitting a mathematical model to it reveals the scale over which the variable behaves predictably. Among the candidate models, an exponential semivariogram fit best, achieving a coefficient of determination of 0.820. Two parameters stood out. The nugget-to-sill ratio of 0.627 indicated that a substantial fraction of the variation occurs at very short distances or within measurement error, a signature of strong local heterogeneity. The range of 37.02 kilometers showed that arsenic values remain spatially correlated over tens of kilometers, implying that broad regional forces, not just field-scale quirks, shape the pattern.

Using that fitted model, the team produced maps of arsenic across the county with ordinary kriging, a geostatistical interpolation technique that weights nearby observations according to the modeled spatial structure. After back-transforming the predictions to the original concentration scale, the maps revealed relatively high arsenic values concentrated mainly in the northern, northeastern, and central-eastern parts of Zhijin County. These hotspots provide exactly the kind of actionable intelligence that environmental agencies need: instead of monitoring uniformly, regulators can focus verification sampling and agricultural inspections on the zones where the geostatistical model suggests arsenic is most likely to be elevated.

Mapping where arsenic is high is only half the story; the other half is explaining why. For this, the researchers employed Geodetector, a statistical framework designed specifically to quantify how much of the spatial variation of a variable can be explained by a categorical environmental factor. The method computes a q statistic that measures explanatory power, and it can also test whether pairs of factors interact to explain more variation together than either does alone. In the Zhijin analysis, several factors showed statistically significant associations with arsenic patterns, with nominal p-values at or below 0.003: mercury concentration, soil organic carbon, distance to mining sites, elevation, distance to rivers, and lithology, the underlying rock type from which the soils developed.

The single most powerful factor was mercury, with a q statistic of 0.4241, the largest of any individual variable tested. The pairing of mercury and arsenic is geologically meaningful, because both elements are often enriched together by the same mineralization and coal-related geological processes that characterize parts of Guizhou. Even more intriguing was the interaction analysis: the combination of mercury and lithology produced the largest joint q statistic of the study, 0.5030, meaning that rock type and mercury together explained half of the spatial variation in soil arsenic. This suggests that arsenic accumulation in Zhijin is not driven by a single cause but by the interplay of geological substrate and geochemical processes that concentrate multiple potentially toxic elements simultaneously.

To push the explanatory analysis further, the team compared five candidate machine learning regressors for predicting arsenic concentrations from environmental covariates, including satellite-derived and terrain-based variables such as Sentinel-2 imagery, elevation from the Shuttle Radar Topography Mission, and river networks from OpenStreetMap. The winner was XGBRegressor, an implementation of extreme gradient boosting, a tree-based ensemble method that builds many sequential decision trees to capture nonlinear relationships. It achieved a pooled out-of-fold coefficient of determination of 0.4197, a root mean square error of 6.7793 milligrams per kilogram, and a mean absolute error of 4.5831 milligrams per kilogram, outperforming the other four models on every metric.

Crucially, the researchers did not treat the model as a black box. They applied SHAP, or SHapley Additive exPlanations, a technique borrowed from game theory that assigns each input variable a signed contribution to every individual prediction. At the global level, SHAP ranked mercury, elevation, and lithology as the three most influential predictors, and mercury and elevation remained among the top three across all five cross-validation folds, indicating that the ranking was stable rather than an artifact of a particular data split. The team then went a step further and produced sample-level SHAP maps, which visualize how each factor pushes arsenic predictions up or down at specific locations. These maps revealed heterogeneous signed contributions across the county, showing that the same factor can raise predicted arsenic in one area and lower it in another, a nuance that global averages completely obscure.

The authors are careful, and appropriately so, about what these findings do and do not prove. Predictive performance was moderate, with the best model explaining roughly 42 percent of the variance, and spatial transferability to other counties remains unverified. More fundamentally, the study distinguishes between spatial explanatory power, predictive contribution, and causality: the fact that mercury and lithology statistically explain arsenic patterns does not by itself establish a causal mechanism. What the results do support is targeted verification and zoned monitoring. For a region where karst hydrology can rapidly transport contaminants and where millions of people depend on local agriculture, that is a pragmatic and valuable outcome. The study also demonstrates a methodological template, pairing geostatistics, Geodetector, gradient boosting, and SHAP-based interpretability, that other regions with complex geology and legacy mining can adapt to trace the hidden geography of soil contamination.

Subject of Research: Spatial distribution and environmental drivers of arsenic in karst agricultural soils

Article Title: Spatial structure and multiscale environmental associations of arsenic in karst agricultural soils of Zhijin County, southwestern China

Article References: Spatial structure and multiscale environmental associations of arsenic in karst agricultural soils of Zhijin County, southwestern China. (n.d.). https://doi.org/10.1007/s10661-026-15960-4

Image Credits: AI Generated

DOI: 10.1007/s10661-026-15960-4

Keywords: arsenic, karst soils, soil contamination, geostatistics, ordinary kriging, Geodetector, XGBoost, SHAP, heavy metals, Guizhou, agricultural soil, environmental monitoring

The research, led by Aaron Jerome Rampersad of the University of Canterbury with Christianne Marie-Claire Faith Zakour of the Loss and Damage Youth Coalition, addresses a problem that has haunted Caribbean hydrology for decades. Rainfall networks across the region are sparse, records are riddled with gaps, and conventional methods for building design rainfall curves quietly assume data that simply do not exist. The stakes are not abstract. Hurricane Ivan in 2004 damaged or destroyed approximately 89 percent of Grenada’s housing stock, inflicting losses near 900 million US dollars, roughly twice the national GDP. Hurricane Beryl in 2024 caused an estimated 218 million dollars in damage and triggered parametric insurance payouts of 55.6 million dollars. Designing infrastructure against the wrong rainfall statistics has direct, measurable consequences.

The team assembled an archive of 28 rain gauges drawing on records from Grenada’s National Water and Sewerage Authority and the Grenada Airports Authority. Before any analysis, the raw material was daunting: 7,671 observed station-days from the water authority network alongside 14,885 daily values from the airport, with most non-airport stations suffering gaps ranging from isolated days to entire missing years. Only Point Salines, with an approximately continuous 40-year record, approached the completeness that standard frequency analysis assumes. For many stations, the eventual curated dataset guaranteed a minimum of 12 years of usable records, a thin foundation for estimating rainfall quantities associated with 50- or 100-year return periods.

The workflow begins with a diagnostic innovation. Rather than relying on the classical semivariogram, the geostatistical workhorse that bins station pairs by separation distance, the authors computed site-specific Pearson correlations directly between every pair of stations. This approach, adapted from techniques used to study non-stationary spatial correlation in earthquake ground motions, sidesteps a known weakness: with so few stations, lag-bin averaging obscures the very structure the analysis is meant to reveal. The verdict was clear. Daily rainfall dependence in Grenada is governed primarily by how far apart two stations are, with elevation acting as a secondary influence whose effect shifts with the seasons. March, one of the driest months, showed strikingly coherent spatial rainfall, while June, as the Intertropical Convergence Zone migrates northward, produced far more scattered behavior.

Those diagnostics fed directly into the gap-filling stage. The team fitted spatial correlation models to the daily rainfall field using Bayesian inference, implemented in Python with the NumPyro library and the No-U-Turn Sampler, treating model parameters as probability distributions rather than fixed numbers. Ordinary kriging driven by these Bayesian-inferred models then reconstructed missing daily values, and it outperformed a full bench of competitors including inverse distance weighting, radial basis functions, and Gaussian process regression. Adding elevation dissimilarity as a covariate delivered marginal but consistent gains. Strict quality control followed: imputations were only retained when at least six donor stations contributed and the kriging prediction variance stayed below 25 percent of the marginal daily variance, lifting every station’s completeness above 70 percent.

With the reconstructed dataset in hand, the researchers turned to extreme value theory. Instead of the Gumbel distribution used in earlier Grenadian studies, which effectively fixes the shape parameter at zero and can underestimate rare rainfall quantiles, they fitted the full generalized extreme value distribution to annual maxima and the generalized Pareto distribution to peaks over threshold, the latter using declustering to ensure that exceedances separated by less than 24 hours counted only once. Because most stations contributed just 12 years of annual maxima, weakly informative priors were specified through an empirical Bayes-style strategy informed by preliminary analytical fits, stabilizing inference on the tail-shape parameter that controls how heavy the rainfall distribution’s upper end truly is. Markov chain Monte Carlo sampling, with four chains and extended warm-up, produced full posterior distributions for every parameter.

The physical signals that emerged are plausible for a mountainous Caribbean island. The shape parameter correlated moderately with elevation, at 0.25 for the generalized extreme value model and 0.34 for the generalized Pareto model, consistent with orographic enhancement steepening the upper tail of the rainfall distribution where moisture-laden trade winds are forced over the central highlands. Interpolating the posterior return levels across the island required choosing among kriging variants, and intrinsic collocated cokriging with a residual correlogram, which exploits elevation as a secondary variable, won on leave-one-out cross-validation. Five poorly constrained stations with only about four years of reliable data were excluded after sensitivity testing showed they destabilized the fitted spatial structures. The final maps show the highest predicted extremes in the island’s northeast, broadly matching Grenada’s known climatic zoning, though accompanied by appropriately large uncertainty estimates.

The study is unusually candid about its own limitations. A sensitivity analysis traced every annual maximum and threshold exceedance back to its source, classifying each as observed or imputed, and then refitted the models using observed extremes only. Where imputed values dominated a station’s extreme sample, return levels shifted dramatically, with differences approaching 80 percent for some generalized Pareto estimates and exceeding 50 percent for some extreme value estimates; at the 25-year return period, most stations stayed within roughly 25 percent. The authors attribute this partly to the smoothing inherent in kriging, which produces conditional-mean predictions rather than stochastic realizations and can dampen localized extremes. They suggest that future work propagate imputation uncertainty directly, through empirical Bayesian kriging, multiple conditional realizations, or hierarchical models treating missing rainfall as latent quantities.

What elevates the work beyond a single-country case study is its transferability. The complete workflow, from site-specific correlation diagnostics through Bayesian model fitting to island-wide interpolation, is publicly available through GitHub and Zenodo repositories, and the underlying Grenada Daily Rainfall Database has been released on Zenodo. The authors also outline an engineering validation path, proposing two-dimensional flood simulations in flood-prone catchments such as St. John’s and Great River to test whether the estimated rainfall fields produce physically reasonable inundation. For small island developing states facing intensifying hurricanes and rising adaptation costs, the message is straightforward: with Bayesian methods, even a fragmented, decades-old network of rain gauges can yield defensible, spatially explicit design rainfall, provided the uncertainties are confronted rather than hidden.

Subject of Research: Bayesian estimation of multi-station rainfall intensity-duration-frequency curves and extreme daily rainfall mapping in data-limited island settings, applied to Grenada

Article Title: A Bayesian workflow for multi-station IDF curve development in data-limited island settings: application to Grenada

Article References: Rampersad, A. J., & Zakour, C. M.-C. F. (2026). A Bayesian workflow for multi-station IDF curve development in data-limited island settings: application to Grenada. Theoretical and Applied Climatology, 157(10), Article 633. https://doi.org/10.1007/s00704-026-06558-4

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06558-4

Keywords: IDF curves, extreme rainfall, Bayesian inference, Grenada, kriging imputation, generalized extreme value distribution, generalized Pareto distribution, spatial correlation, small island states, geostatistics, return levels, climate risk

The study area lies between roughly 2°57′ and 3°26′ N and 10°14′ and 10°41′ E on the southern Cameroon Plateau, at an average elevation of about 550 metres, where the Lokoundjie River and its tributaries drain a landscape of Archean to Paleoproterozoic basement. Geologically, the region belongs to the Nyong Group, a reactivated segment of the Congo Craton’s northwestern margin that was transformed during a Paleoproterozoic tectono-metamorphic event around 2050 million years ago. The Nyong Group hosts greenstone-related lithologies including pyroxenites, amphibolites, peridotites, talc schists and banded iron formations, alongside foliated tonalite-trondhjemite-granodiorite suites, orthogneisses, granodiorites and syenites. This cratonic terrane has long attracted prospectors: gold occurrences are documented across southern Cameroon around Bipindi, Lolodorf and Akom II, and artisanal miners routinely work alluvial gravels and altered quartz veins along an established gold corridor.

Fieldwork centred on the tributaries of the Tyango River, where the researchers collected six fresh outcrop samples to characterise the basement and ten alluvial sediment samples, labelled BIP-01 to BIP-10, from hand-dug pits in active riverbed deposits at depths of 50 to 100 centimetres, targeting gravel-rich horizons. Thin-section petrography revealed three principal basement lithologies: dark grey, weakly foliated pyroxene-epidote gneisses with heterogranular granoblastic textures and abundant pyroxene and epidote; massive, fine- to medium-grained amphibolites dominated by amphibole with secondary epidote replacing it; and whitish to grey-black quartzites exposed near the confluence of the Nyaba’ah and Tyango rivers, composed mainly of quartz and feldspar with muscovite, rare pyroxene relics and opaque minerals. Crucially, optical microscopy did not confirm any discrete gold- or platinum-bearing grains in these rocks, so the lithologies serve as provenance indicators rather than proven ore sources.

The sedimentological analysis painted a picture of a proximal, texturally immature system. Granulometric sieving showed that most samples are dominated by fine to medium fractions between 0.5 and 0.063 millimetres, with poor to moderate sorting and cumulative curves whose slopes range from steep to gentle. Steeply declining curves in samples such as BIP-01 and BIP-07 record high-energy deposition in fast-flowing water, while fine-dominated curves in BIP-04 and BIP-10 point to quiet, lake- or floodplain-like settings. Histograms revealed bimodal distributions in six samples, with coarse particles concentrated at pit bottoms beneath fines, a pattern consistent with density-driven sorting. Quartz grain morphoscopy proved especially telling: very angular to angular grains make up the overwhelming majority of all samples, with some samples containing up to 96 percent very angular grains and low sphericity throughout, indicating that the sediment travelled only short distances from nearby metamorphic sources with negligible mechanical wear.

Heavy-mineral concentrates extracted from the sediments were dominated by opaque minerals, which account for about 56.67 percent of the assemblage, and pink, prismatic to pyramidal zircon at roughly 30.67 percent, with subordinate garnet, epidote, hornblende, diopside, kyanite, sillimanite, andalusite and mica. This mix points to short transport from heterogeneous metamorphic source rocks. The researchers are careful to stress, however, that without reflected-light microscopy, scanning electron microscopy with energy-dispersive spectroscopy, or electron microprobe data, the opaque grains cannot yet be classified as platinum minerals or gold-bearing phases, and the heavy minerals should be read as provenance and hydraulic concentration indicators rather than established pathfinders for gold in Bipindi.

Bulk-sediment geochemistry, performed at ALS Global in Vancouver using aqua regia digestion and inductively coupled plasma mass spectrometry with certified reference materials, added a chemical dimension. Aluminium oxide contents are low, below 2.31 percent, and titanium oxide ranges from 0.05 to 0.19 percent, while iron oxide is markedly enriched upstream, reaching 25.16 percent, and declines downstream, a trend the authors attribute to alteration and transport. Chemical index of alteration values mostly exceed 70 percent and climb as high as nearly 96 percent, indicating moderate to intense chemical weathering under the humid tropical climate. Upstream samples show aluminium-to-sodium ratios reaching 231, evidence of severe sodium leaching, while high thorium-to-uranium ratios above the upper continental crust average of about 3.8 confirm uranium loss during weathering. Provenance discrimination based on aluminium-to-titanium ratios, thorium-versus-scandium plots and lanthanum-versus-thorium plots indicates a mixed mafic to felsic source, consistent with derivation from the gneisses, amphibolites, quartzites and tonalite-trondhjemite-granodiorite lithologies of the Nyong Group, with only minimal sediment recycling.

The precious-metal results were striking in their asymmetry. Gold concentrations range from 0.0001 to 0.243 parts per million, but platinum remains at or below 0.001 parts per million and palladium at or below 0.003 parts per million, effectively at detection limits. The authors interpret the gold distribution cautiously as a local alluvial anomaly rather than evidence of substantial mineralization, and they explicitly decline to claim platinum-group-element mineralization without direct mineralogical confirmation. A Pearson correlation matrix reinforced this restraint: aluminium and iron oxides correlate strongly, as do zinc and copper, but gold, palladium and platinum show no strong positive relationships with the main lithogenic elements, indicating that precious-metal contents are low, discontinuous and weakly coupled to bulk-sediment chemistry.

To convert these point measurements into a spatial picture, the team built a geographic information system database in ArcGIS 10.8 and produced interpolated gold distribution maps and three-dimensional visualisations in Surfer 16, applying ordinary kriging guided by directional semi-variograms. The statistics revealed a strongly positively skewed distribution with a mean of 0.0248 parts per million, in which 90 percent of samples fall in a low-grade class below 0.0608 parts per million while the remaining 10 percent, averaging 0.2127 parts per million, occupy a high-grade class between 0.1823 and 0.243 parts per million. The fitted spherical semi-variogram model combines a nugget effect of 0.0018 with a sill variance of 0.0045 and a range of about 3.83, oriented 21.62 degrees toward the south, with spatial correlation fading beyond roughly 13.76 degrees in the southern direction. The strong nugget component reflects short-scale variability, sparse sampling and analytical noise, which is why the kriged maps are presented as exploratory guides rather than resource models. Nevertheless, the interpolation consistently shows higher gold values toward the southern part of the study area and at lower elevations, a pattern consistent with alluvial concentration and aligning with previously documented gold showings in altered rocks around Akom II.

The authors are candid that ten sediment samples cannot establish structural control or prove a mineralized body, and they warn against over-reading the apparent north-south trend until structural measurements, lineament analysis and bedrock lithogeochemistry are integrated. What the study delivers instead is a disciplined exploration framework: the southern sector toward Akom II and Grand Zambi emerges as the clear priority for denser sediment sampling, seasonal monitoring, structural mapping and mineralogical confirmation of opaque grains by scanning electron microscopy or electron microprobe. In a region where artisanal miners have long worked the rivers on intuition, this fusion of microscopic petrography, weathering geochemistry and geostatistics offers something more valuable than a quick strike, a transparent, testable map of where the next phase of exploration should dig.

Subject of Research: Integrated geochemical, sedimentological and geostatistical assessment of gold dispersion in alluvial sediments at Bipindi, southern Cameroon

Article Title: Integrated geochemistry, geostatistics, and sedimentology to identify potential gold-bearing zones at Bipindi, southern Cameroon

Article References: Gake Belle, R., Mbanga Nyobe, J., Mbabi Bitchong, A., Nga Essomba Tsoungui, P. E., Mimba, M. E., & Ndip Ojong, E. (2026). Integrated geochemistry, geostatistics, and sedimentology to identify potential gold-bearing zones at Bipindi, southern Cameroon. Discover Geoscience, 4(1), Article 350. https://doi.org/10.1007/s44288-026-00710-3

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00710-3

Keywords: gold exploration, Bipindi, southern Cameroon, Congo Craton, stream sediment geochemistry, geostatistics, kriging, heavy minerals, sedimentology, Nyong Group, platinum-group elements, alluvial gold