The research, published in Discover Geoscience by Dionisio Uendro Carlos, César Augusto Moreira, and Guilherme Libardi Maestrelo, introduces a moving-window ordinary kriging (MW-OK) algorithm designed specifically for airborne and ground radiometric data. Conventional interpolators such as inverse distance weighting and minimum curvature are fast and widely used, but they share a fundamental weakness: they do not explicitly account for the local spatial structure of the data. In practice, that means they can oversmooth short-wavelength variability, blurring exactly the sharp radiometric gradients that mark boundaries between rock units or zones of hydrothermal alteration.

Kriging, by contrast, is a geostatistical method built on a statistical model of how a property varies in space. Because radioactive decay is inherently a stochastic process, the authors argue, an interpolator that respects this randomness is desirable. Yet standard kriging typically relies on a single global variogram, the mathematical function describing spatial correlation, assumed to hold across the entire survey area. Radiometric variables rarely obey that assumption at regional scale. Lithological boundaries, weathering domains, and survey-related effects all conspire to make the spatial structure change from place to place, so a global model can over-smooth local features or produce biased estimates where the structure shifts.

MW-OK sidesteps this problem by abandoning the global model altogether. For every prediction node on the output grid, the algorithm defines a circular window of investigation, in this case with a radius of 750 meters, roughly 1.5 times the 500-meter flight-line spacing, so that each window intersects two adjacent flight lines. Within that window, the method assumes second-order stationarity holds only approximately and locally. An empirical semivariogram is computed from all unique pairs of samples in the neighborhood, binned into eight lag classes, with any bin containing fewer than three valid pairs rejected. An exponential model is then fitted to each local semivariogram by bounded non-linear least squares, yielding a nugget effect, partial sill, and range parameter tailored to that specific patch of ground.

The computational scale is striking. The study area, a 40 by 40 kilometer subset of the Sul-Riograndense Shield Aerogeophysical Project flown by the Geological Survey of Brazil in 2010, contains 50,374 measurement points. Because the algorithm estimates a local variogram for each measurement point, it effectively generates 50,374 variograms. Neighbor retrieval is accelerated with a k-dimensional tree using ball queries, and safeguards keep the local kriging systems well conditioned: a minimum of eight neighbors, a maximum of thirty, removal of duplicate coordinates to prevent singular matrices, and a tiny ridge term added to the kriging matrix diagonal for numerical robustness. If a non-linear variogram fit fails, the workflow falls back on physically sensible defaults derived from the local data variance and median nearest-neighbor spacing.

The fitted local variograms themselves carry geological information. For potassium, the nugget effect was essentially zero at the scale of the data support, meaning microscale variability and measurement noise are minor compared with the spatially structured component, and correlation decayed over a practical range of about 1.3 kilometers. Equivalent uranium showed an even longer continuity, around 1.7 kilometers, consistent with smoother spatial behavior, while equivalent thorium displayed a larger partial sill, indicating higher structured variance, with a moderate range near 1.5 kilometers. These parameters confirm that, within each window, the radiometric variability is predominantly spatially structured rather than random noise.

When the researchers applied MW-OK to the survey data, the resulting potassium, uranium, and thorium maps reproduced the broad patterns that geologists expect. Potassium enrichment above 4 percent traced the São Sepé Granite in the north, elevated values of roughly 4 percent marked the Ramada Granite in the southwest, and low concentrations coincided with the metatuff-dominated Bossoroca Complex and the sandstones of the Maricá Group. Thorium highs of about 21 parts per million again outlined the São Sepé Granite, with additional anomalies over the Cerro da Cria Granite and nearby rhyolites. But the kriged surfaces went further than the conventional minimum-curvature grids: they better represented local gradient variations near lithological contacts and revealed internal heterogeneity within the felsic bodies that the smoother deterministic algorithm muted. Notably, three east-west linear artifacts aligned with the survey tie lines, visible in the minimum-curvature uranium map, were suppressed in the kriged surface.

A key advantage of kriging is that it also delivers an estimate of its own uncertainty. The kriging variance maps for all three radioelements bore the clear imprint of the acquisition geometry: uncertainty was minimized along and near the north-south flight lines and increased in the gaps between them, telling interpreters exactly where the estimates are most reliable. To quantify how closely MW-OK matched the industry-standard approach, the team performed pixel-wise comparisons with the minimum-curvature grids. The agreement was remarkably strong, with coefficients of determination of 0.9938 for potassium, 0.9727 for equivalent uranium, and 0.9946 for equivalent thorium, and residual histograms tightly centered near zero. Ordinary least-squares fits revealed a subtle but systematic amplitude scaling, with slopes slightly below unity, suggesting MW-OK modestly attenuates the highest values while remaining close to the identity overall. The authors are careful to frame this as evidence of methodological consistency rather than broad superiority: the two methods agree at the domain scale, and their differences emerge primarily at local scales.

Beyond the three radioelement grids, the workflow generates the interpretive products that exploration geologists actually use. It computes standard radiometric ratios (K/eTh, K/eU, eU/eTh) and the F-parameter, a widely used proxy for hydrothermal alteration that relates potassium abundance to the uranium-to-thorium ratio. It also produces ternary radioelement maps in three complementary forms: the classical K-eU-eTh display, normalized ternaries, and compositional-data transformations using additive, centered, and isometric log-ratios that respect the closure constraint and reduce spurious correlation when treating the three elements as parts of a composition. In the study area, a historically significant gold district in southern Brazil where mineralization is hydrothermal and structurally controlled in quartz veins hosted by schists and granites, these products delineated the main lithological contacts and showed spatial correspondence with known gold-mineralized zones and areas of gold enrichment identified by stream-sediment prospecting.

Perhaps most importantly for the wider community, the entire workflow is open source and exports standard deliverables, including GXF grids, GeoTIFF rasters, KMZ overlays, and PNG figures, with proper spatial metadata, making it directly interoperable with both open-source tools such as Python, QGIS, and Google Earth and commercial platforms like Oasis Montaj, Intrepid Geophysics, and ArcGIS. The authors acknowledge the method’s limitations: effective spatial resolution remains constrained by flight-line spacing and sampling interval, computational cost grows with finer grids and larger neighborhoods, and sparsely sampled sectors can yield under-populated windows with unstable variograms, which the workflow mitigates through minimum-neighbor thresholds and k-nearest-neighbor fallbacks. Looking ahead, they plan to incorporate co-regionalized cokriging constrained by lithology and structure, and deep-kriging approaches, within the same moving circular-window framework. For now, MW-OK stands as a robust, locally adaptive alternative that offers incremental rather than transformative gains over deterministic gridding, but gains that could make the difference between a blurred contact and a clearly resolved exploration target.

Subject of Research: A locally adaptive geostatistical interpolation method for processing airborne and ground gamma-ray spectrometric data

Article Title: A moving-window ordinary kriging approach for airborne and ground gamma-ray spectrometric data processing

Article References: Carlos, D. U., Moreira, C. A., & Maestrelo, G. L. (2026). A moving-window ordinary kriging approach for airborne and ground gamma-ray spectrometric data processing. Discover Geoscience, 4(1), Article 289. https://doi.org/10.1007/s44288-026-00674-4

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00674-4

Keywords: gamma-ray spectrometry, kriging, geostatistics, radiometric mapping, mineral exploration, semivariogram, interpolation, gold prospecting, geological mapping, Brazil, open-source software, compositional data

The research team, led by Carina Spengler and Simon Saurbier at the Karlsruhe Institute of Technology together with colleagues including Andreas Lindenmann, Dieter Krause, and Sven Matthiesen, tackled a stubborn gap in the science of hand-arm vibration. Prolonged exposure to tool vibration is a well-documented occupational hazard, linked to vascular, neurological, and musculoskeletal disorders, including the painful condition known as white finger syndrome. The international standard ISO 5349-1 governs how this exposure is measured, and a companion standard, ISO 10068:2012, provides reference values for the mechanical impedance of the hand-arm system—the ratio of dynamic force to vibration velocity at the point where the hand grips the handle. But those reference values were built on limited grip and push force ranges, and earlier studies typically tested only a handful of force combinations with small groups of subjects.

What makes the new work stand out is its systematic scope. Rather than testing two or three discrete force settings, the team measured the biodynamic response across a full ten-by-ten grid of grip and push forces, with grip forces spanning 10 to 135 newtons and push forces spanning 10 to 110 newtons. Five subjects—three men and two women aged 23 to 34—each completed 300 measurements over multiple sessions spread across roughly twelve days, with the order of force combinations randomized to guard against fatigue and learning effects. In total, 1500 measurements were recorded on a test rig equipped with two electromechanical shakers, one axial and one rotational, driving a standardized cylindrical measuring handle instrumented with one-axial and three-axial force sensors and accelerometers.

The excitation itself was engineered for efficiency and safety. Instead of sweeping through frequencies one at a time, the shaker delivered a multisine signal—a superposition of sine waves from 10 to 500 hertz in 5-hertz steps—which compressed each measurement to just 16.5 seconds. The unweighted root-mean-square acceleration was 30 meters per second squared, far higher than the 2.5 to 9 meters per second squared typical of earlier single-sine studies, and the calculated total vibration value of 6.5 meters per second squared per measurement remained well below the exposure limits that would trigger protective measures. Subjects watched a real-time graphical display of their applied grip and push forces against the targets, adjusting their effort as needed, and any measurement deviating more than 5 newtons or 14 percent from the target was discarded—150 of the 1500 recordings failed this test.

Two quantities anchored the analysis. The mechanical impedance, expressed in newton-seconds per meter, describes how the hand-arm system dynamically resists motion: its real part reflects energy dissipated per unit time, while its imaginary part captures energy that is reversibly stored and returned. The vibration power absorption, measured in watts, is derived from the real part of the impedance and the squared vibration velocity, and quantifies the actual mechanical power flowing into the body. Both were evaluated along the translational z-axis of the hand-arm system, the direction defined in ISO 5349-1 as the primary axis of exposure assessment.

The results reveal a striking frequency-dependent choreography. Within a single subject, varying the grip and push forces broadened the median impedance magnitude from values between roughly 100 and 430 newton-seconds per meter into a range spanning 50 to 950 newton-seconds per meter across the measured frequencies; for one subject the spread reached a 15.8-fold increase at 50 hertz. Grip force proved the dominant parameter, strongly influencing impedance between 30 and 70 hertz and again above 160 hertz, with a maximum difference of 318 newton-seconds per meter at 40 hertz and a difference of 290 newton-seconds per meter persisting even at 500 hertz. Push force mattered most in the mid-range between 50 and 110 hertz, peaking at 177 newton-seconds per meter at 55 hertz, but its influence faded above 160 hertz—consistent with the fact that grip is generated by the hand and forearm while push originates from the shoulder, and high-frequency vibrations barely propagate beyond the hand itself.

The power absorption data told an equally nuanced story. Most of the absorbed energy concentrated below 60 hertz, precisely the range in which hand-held tools transmit the most vibration to the body, and the curves peaked between 35 and 40 hertz before collapsing toward zero above 100 hertz. Higher forces generally increased absorption, but in a narrow low-frequency window—roughly 10 to 30 hertz—the trend reversed, with stronger grips actually reducing the energy absorbed. The researchers attribute this to stiffening: as muscle contraction increases, the hand-arm system becomes a harder spring, reflecting more energy rather than dissipating it, while simultaneously improving mechanical coupling that channels more vibration into the arm at higher frequencies. Notably, the two female subjects showed lower impedance values above 100 hertz and different resonance features than the male subjects, echoing earlier reports that anthropometry and sex-related differences shape the biodynamic response.

Perhaps the most practically consequential finding concerns prediction. Because the force-dependent curves proved approximately linear at nearly every frequency, the team tested whether intermediate force levels could be interpolated from a small number of measured reference curves. The answer was emphatic: with just three carefully spaced grid curves, linear interpolation predicted both impedance magnitude and power absorption at intermediate forces with a mean absolute percentage error below 5 percent, corresponding to roughly 95 percent accuracy. Adding a fourth curve improved matters further—dropping the mean error for impedance to 3.29 percent—but with diminishing returns, and the largest accuracy gain came from moving from two to three curves. For future experiments, this means the grueling full-grid protocol may be unnecessary: a handful of reference force levels, distributed across the range of interest, could capture the essential behavior while dramatically reducing the burden on human subjects.

The study also delivers a pointed critique of the existing standard. When the measured impedance curves were compared against the ISO 10068:2012 reference range, they frequently fell outside it, especially between 20 and 100 hertz and above 250 hertz—exactly the frequency bands where grip and push forces exert their strongest influence. The authors argue that interaction forces should therefore not be treated as a fixed boundary condition in vibration analysis but as an explicit modeling parameter, one that varies continuously during real tool use. Their data, published as an open dataset, and their interpolation framework pave the way for physical and simulation-based hand-arm models that incorporate user forces directly—enabling virtual testing of tools before prototypes exist, robotic replication of human dynamics, and ultimately the design of machines that shake less hard on the people who hold them.

Subject of Research: Effects of grip and push force variations on the mechanical impedance and vibration power absorption of the human hand-arm system

Article Title: Modeling the influence of grip and push force variations on the mechanical impedance and vibration power absorption of the hand-arm system

Article References: Spengler, C., Saurbier, S., Lindenmann, A., Doellken, M., Bunk, L., Krause, D., & Matthiesen, S. (2026). Modeling the influence of grip and push force variations on the mechanical impedance and vibration power absorption of the hand-arm system. Results in Engineering, 32, Article 113151. https://doi.org/10.1016/j.rineng.2026.113151

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113151

Keywords: hand-arm vibration, mechanical impedance, vibration power absorption, grip force, push force, biodynamics, power tools, occupational health, ISO 10068, vibration exposure, interpolation, human-machine systems