The difficulty with the standard approach is well known. In general relativity, solving a system of interacting bodies means finding the spacetime metric while simultaneously computing how the bodies move through that metric. The two problems are coupled, and the coupling grows ferociously with each additional body. Beyond the two-body case, explicit solutions become impractical, and researchers typically resort to large-scale numerical simulations of the Einstein field equations, slicing time and continually updating the metric as the system evolves. Silverberg and Eischen take a different path. Rather than reorganizing the geometry of general relativity, as the ADM or teleparallel formulations do, they return to the classical tradition in which gravity is expressed through force relations between bodies, updated to match relativistic predictions.

Their formulation, which they call general mechanics, rests on what they term the parity hypothesis: that a force-based description of gravity can be mathematically faithful to the trajectories predicted by general relativity, at least in the absence of external radiation sources and gravitational waves. The centerpiece is a relativistic universal law of gravitation that modifies the familiar inverse-square law by the factor 1 + 3(h_R/cr)^2, where h_R is the specific relativistic angular momentum of the pair, c is the speed of light, and r is the separation of the bodies. This factor injects a rotational component of energy that becomes significant only when motion is relativistic. In the slow, weak-field limit it vanishes and Newton’s law reappears. When retained, the law reproduces the classic relativistic signatures, including perihelion precession, light deflection, and the photon sphere, in exact agreement with the Schwarzschild two-body solution.

The new contribution is the gravitational frame itself. For each interacting pair of bodies, the authors define the gravitational frame as the reference frame in which the pair’s total relativistic linear momentum vanishes, mathematically the familiar center-of-mass frame but reinterpreted as something far more consequential. In their view, this is not merely a convenient coordinate choice but the frame in which nature actually formulates gravitational interaction. The law of gravity, they argue, must first be evaluated in this pairwise frame, and only then Lorentz-transformed into whatever global frame the problem requires. The distinction is subtle but, they contend, physically essential.

The authors support the hypothesis with an energetic argument. The one-body formulation of the two-body problem, the setup Schwarzschild himself used, is independent of the velocity of the system’s mass center, so it corresponds to a whole family of possible two-body problems differing in their net linear momentum. Comparing the relativistic kinetic energies of the one-body and two-body descriptions, the two are equal only when the pair’s relativistic linear momentum is zero, that is, only in the gravitational frame. Among all admissible two-body solutions, only this one preserves energetic equivalence between the two formulations, which the authors take as strong evidence that nature selects this frame. They also note that the historic tests of general relativity, from Mercury’s perihelion to Eddington’s 1919 eclipse expedition, were implicitly performed in exactly such frames, where a dominant source and a light test body leave the pair with effectively zero net momentum.

To extend the idea to N bodies, the authors apply the gravitational law pairwise, computing each interaction in the gravitational frame of the relevant pair and then transforming the resulting accelerations into a common global frame using Lorentz transformations. The paper also defends the need for a global frame at all, using a matrix argument showing that frames sharing Lorentz transformability to a common global reference yield consistent physics, a point illustrated through the twin paradox. Evolution proceeds with respect to proper time in flat Minkowski spacetime rather than the coordinate-time slicing of curved-spacetime relativity, with radiation and gravitational-wave degrees of freedom set aside so that the conservation laws of mechanics survive.

The real test comes when the source is not a single point mass but a rotating, extended body. In general relativity, a spinning source drags spacetime around with it, the frame-dragging effect described by the Kerr metric, so that a body falling radially inward develops an azimuthal drift. A naive force law applied to the whole source as one object would miss this entirely, predicting purely radial infall. Silverberg and Eischen show that when the source is instead treated as an aggregate of moving constituents, each interacting pairwise in its own gravitational frame, the tangential deflections emerge naturally. No Kerr-like metric term, no extra transverse force, and no separately imposed frame-dragging term is added; the rotational behavior arises from the frame-selection rule itself.

The numerical experiments are striking. In one set, a test body starting at three Schwarzschild radii and moving inward at 0.8c approaches a source with an artificially pinned upward velocity. With a stationary source the body falls straight in; at 0.5c it deflects upward, and at 0.9c the deflection becomes pronounced. In the flagship example, the source is a rigid ring of 800 point sources, with total mass equal to the Sun’s, rotating at 0.95c at a radius of one-tenth the Schwarzschild radius. A test body aimed radially inward is swept sideways by roughly 0.044 Schwarzschild radii as it crosses the ring, deflected in the direction of rotation, exactly the qualitative signature of frame dragging. Near the photon sphere at 1.5 Schwarzschild radii, a light-speed test body orbiting a rotating ring shows increased attraction as the ring spins faster, with small but measurable differences between prograde and retrograde configurations, the retrograde orbit radius slightly larger than the prograde one, consistent in character with Kerr-type behavior.

The authors are careful about scope. They emphasize that their rotating-ring examples are behavioral tests, not quantitative comparisons with the full Kerr solution, which is a three-dimensional, axisymmetric black-hole exterior with an intrinsic spin parameter, horizon conditions, and a fixed multipole structure. In their formulation, rotation is not imposed through a spin parameter at all but constructed from the motion of constituent masses, much as rigid-body motion is built from moving point masses in classical mechanics. A decisive quantitative comparison with Kerr, they acknowledge, would require extending the planar formulation to fully three-dimensional rotating sources and carefully translating Kerr’s spin, horizon, and multipole assumptions into the mechanics framework. They also distinguish their exact two-body law from post-Newtonian approximations, which are systematically accurate only in the weak-field, slow-motion limit and become strained as velocities approach the speed of light.

Still, the implications are tantalizing. If the gravitational-frame construction holds up under deeper scrutiny, relativistic N-body dynamics, from merging compact binaries to accretion disks threading spinning black holes, could one day be attacked with the same force-based, pairwise computational machinery that has served celestial mechanics since Newton, integrated here with a fourth-order Runge-Kutta scheme in proper time. The authors argue that this transition mirrors the shift toward modern N-body computational methods that transformed other fields of mechanics over the past half-century. Whether general mechanics can ultimately match the precision of full numerical relativity remains an open question, but the paper offers a provocative demonstration that frame dragging, strong-field precession, and spin-enhanced attraction can emerge from a force law and a well-chosen frame, without a single line of curved-spacetime geometry.

Subject of Research: A mechanics-based, force-law formulation of relativistic gravitation extended to N-body systems via pairwise gravitational frames

Article Title: The gravitational frame for solving relativistic N-body problems

Article References: Silverberg, L. M., & Eischen, J. W. (2026). The gravitational frame for solving relativistic N-body problems. Celestial Mechanics and Dynamical Astronomy, 138(5), Article 60. https://doi.org/10.1007/s10569-026-10326-x

Image Credits: AI Generated

DOI: 10.1007/s10569-026-10326-x

Keywords: general relativity, N-body problem, gravitational frame, frame dragging, Schwarzschild solution, Kerr metric, celestial mechanics, special relativity, Lorentz transformation, relativistic gravitation, numerical simulation, rotating black holes

A credit transfer system is, at its core, an accounting mechanism for learning. It records verified educational achievements so that a course completed at one institution can be recognized at another, supporting lifelong learning in an economy where skills rapidly become obsolete. Similar systems have been tried internationally, from the European Credit Transfer and Accumulation System to South Korea’s Academic Credit Bank and proposals for aligned Asian academic credits. Yet the Chinese version is distinctive in scale and in the complexity of its governance: it must coordinate individual learners deciding whether to participate, schools deciding how much to invest in compatible teaching resources, and a government deciding how aggressively to supervise and subsidize the enterprise. Each actor’s payoff depends on what the others do, which is precisely the situation game theory was built to describe.

He and Bu constructed a three-subject dynamic evolutionary game model. In an evolutionary game, agents are not treated as perfectly rational calculators who instantly find the best move; instead, populations of players adjust their strategies over time, with more successful strategies spreading through imitation and reinforcement. The approach, which traces back to Lewontin’s application of game theory to evolutionary biology and has since become standard in economics, is well suited to education policy, where participants learn, imitate, and adapt rather than optimize from the start. The researchers incorporated external constraints, strategic assumptions, and payment assumptions into the model, grounding the payoff structure in the benefit relationships among the stakeholders.

The model assigns each player a binary strategic choice. Learners can adopt strong participation, actively engaging with the system to improve their abilities, or hold back. Schools can participate actively by building high-quality teaching resources that feed the credit system, or invest minimally. The government can exert strong dominance through supervision and support, or take a laxer stance. Each choice carries costs and expected benefits: learners pay tuition and effort in exchange for portable credentials and improved skills; schools bear the expense of developing quality courses and administrative compatibility in exchange for enrollment, reputation, and possibly government incentives; the government funds oversight and subsidies in exchange for social returns such as a better-trained workforce.

To move beyond abstract equilibrium conditions, the authors ran numerical simulations of the model. Simulations of this kind allow researchers to trace how the probability of each strategy changes over successive rounds, revealing the decision-making mechanism of the three parties under different parameter settings. The exercise is particularly valuable because the equilibria of evolutionary games depend sensitively on cost allocations and benefit magnitudes; a system that looks stable on paper can collapse if, for example, the cost of building compatible teaching resources falls entirely on schools while the benefits accrue mostly to learners and the state.

The central finding is a matching rule. The optimal configuration for constructing a suitable credit transfer system combines strong participation by learners focused on ability improvement, active participation by schools in creating high-quality teaching resources, and strong dominance by the government in supervision and support. In other words, no single actor can carry the system: learner enthusiasm without quality courses produces credentials with little content, and quality resources without learners generate waste, while absent government supervision undermines trust in the credits themselves. The system functions as a three-legged stool, and the evolutionary dynamics determine whether all three legs stabilize simultaneously.

The analysis also illuminates why credit transfer efforts elsewhere have struggled. International experience, from pan-European grading scales to attempts at equivalency between the European credit systems for higher and vocational education, shows that technical compatibility is only part of the problem; incentives matter just as much. The Chinese study suggests that misalignment of cost and benefit among stakeholders, rather than flaws in credit accounting itself, may explain stalled adoption. When the government’s supervision and support are strong enough to shift expected payoffs, the simulations indicate that learners and schools can be tipped into the cooperative strategies that make the whole arrangement viable.

Methodologically, the paper illustrates a growing trend in education research: importing formal tools from economics and biology to model policy systems with many interacting agents. The authors draw on stakeholder theory, which originated in management science with Donaldson and Preston’s influential formulation, to justify the three-party structure, and then embed it in a replicator-style dynamic framework. Numerical simulation serves as the bridge between analytic equilibrium conditions and practical policy levers, letting the researchers ask counterfactual questions about what happens when costs are reallocated or benefit demands change.

For policymakers, the practical implications are concrete. Investment in high-quality teaching resources should be encouraged so that schools find active participation genuinely rewarding, and government supervision should be robust enough to certify credit quality, since the value of any transferable credit rests on trust in its standards. Meanwhile, the design should make ability improvement, not merely credential collection, the visible payoff for learners, aligning individual motivation with the system’s educational purpose. The work was supported by China’s National Program for Funding Postdoctoral Researchers, and both authors are affiliated with Beijing Normal University.

As lifelong learning becomes an economic necessity in an era of rapid technological change, credit transfer systems are likely to spread globally, and the Chinese experience will be watched closely. He and Bu’s evolutionary game model provides a template for anticipating how such systems behave before they are built: identify the players, quantify their costs and benefits, simulate the dynamics, and design interventions that steer all parties toward the cooperative equilibrium. The finding that strong learner participation, active school engagement, and strong government dominance must arrive together is a warning against piecemeal reform, and a demonstration that the mathematics of strategic interaction can inform one of education’s most intricate institutional projects.

Subject of Research: Game-theoretic modeling of stakeholder participation strategies in China's credit transfer system for lifelong learning

Article Title: Participating Strategy of the Constructors in the Construction of Credit Transfer System in China

Article References: Participating Strategy of the Constructors in the Construction of Credit Transfer System in China. (n.d.). https://doi.org/10.1007/s44366-026-0077-z

Image Credits: AI Generated

DOI: 10.1007/s44366-026-0077-z

Keywords: credit transfer system, evolutionary game theory, lifelong learning, China, stakeholder theory, education policy, learners, schools, government supervision, numerical simulation, Beijing Normal University, Participating

The research focuses on what mining engineers call the goaf-side narrow coal pillar, or GSNCP, a thin ribbon of coal left standing between an active tunnel and a mined-out void. In extra-thick coal seams averaging eleven meters, these pillars are routinely riddled with legacy roadways from earlier, less regulated mining eras. When a new roadway, known as a goaf-side entry, is driven alongside such a compromised pillar, the results on site have been alarming: severe rib heaving, roof convergence, rib collapse, and even complete roof caving in zones where abandoned roadways cluster together. Traditional support methods such as bolting and grouting, the authors note, fail entirely under these conditions.

To understand why, the team built a large three-dimensional numerical model using FLAC3D 7.0, spanning 500 by 500 meters laterally and 300 meters in depth, discretized into half-meter cubic elements with more than 7.6 million cells in the critical study area alone. The model reproduced the actual excavation sequence of the mine, including dynamic compaction of the adjacent goaf using double-yield elements and stepwise extraction of the working face at one-meter cycles. The coal was governed by a strain-softening criterion while the overlying strata followed the Mohr-Coulomb failure criterion, allowing the researchers to track how stress redistributes as abandoned roadways slice the pillar into segments.

The simulations revealed a striking sensitivity to the spacing between abandoned roadways. When intervening coal ribs between two roadways were narrower than five meters, the segments entered a low-strength bearing state and simply shed their load onto neighboring pillars, triggering severe rib heaving. At a width of seven meters, the internal abutment pressure peaked at a dangerous 27.9 megapascals in a single sharp spike, a condition the authors describe as an ultra-high strength bearing state in which the pillar risks internal plasticization and fracturing. Only when spacing exceeded nine meters did the stress curve split into a healthier bimodal shape, with the peak dropping to 23.1 megapascals and the load distributed across a broad high-strength bearing zone.

The second half of the puzzle lies overhead. Using the classical OX break theory of main roof behavior, the team developed a theoretical model linking the arc-shaped triangular roof block that forms above the goaf edge to the abandoned roadways hollowing out the pillar beneath it. This triangular block, whose dimensions follow the periodic weighting span of the main roof, rotates and subsides around a fracture line after the roof breaks. The researchers identified nine possible relative positions among the solid coal rib, the triangular block, and the abandoned roadways, and derived moment-balance equations for each configuration to determine when the block remains in equilibrium and when it fails.

The analysis produced two instability coefficients with clear physical meaning. When the coefficient K2 exceeds one, the triangular block undergoes rotational deformation instability, crushing the rock at its corners. When the coefficient K1 falls between zero and one, the block slides instead. Crucially, the calculations showed that as more abandoned roadways accumulate beneath a single triangular block, both coefficients decrease simultaneously, meaning the weakened pillar can no longer supply the support force needed to prevent sliding. Because at most three abandoned roadways can fit beneath one block, the team showed that the failure mode shifts from rotation to sliding as roadway density increases, with the block’s gravitational load transferring through the roadway roof to the solid coal rib and producing uncontrollable deformation, roof fracturing, or even roof cutting.

These theoretical predictions were tested in a 1:100 scale physical analogue model built from layered sand, lime, and gypsum mixtures matched to the site’s borehole geology, with hydraulic jacks applying 34.6 kilopascals to reproduce the overburden load at the mine’s 250-meter depth. The experiments confirmed the mechanism in all three fracture-line scenarios: with backfill in place, the pillar and roof stayed intact, but removing the backfill reproduced the full disaster pattern of pillar crushing and severe roof collapse observed in the field. The observed diagonal shear fractures, angled at roughly 55 to 60 degrees, matched the theoretical shear angle predicted by the Mohr-Coulomb criterion for coal with internal friction angles between 20 and 30 degrees.

Having diagnosed the disease, the researchers prescribed a cure that is as economical as it is effective: backfill the abandoned roadways with the very coal excavated during roadway driving. Laboratory tests determined the optimal recipe, mixing crushed coal, P425 Portland cement, and water in a ratio of 6 to 1 to 0.7. With a particle gradation of roughly 60 percent fines under five millimeters, the solidified mixture achieved a uniaxial compressive strength of 9.8 megapascals, about 80 percent of the raw coal strength of 12.3 megapascals, while still meeting the pumping requirements of mining concrete pumps. The mixture is pumped into mining-purpose backfill bags that prevent leakage into the goaf, with roof bolts fitted with end caps to avoid puncturing the bags.

Numerical simulations of the backfilled system showed a clear strength threshold. At just 20 percent of coal strength, the backfill carried little load and stress remained dangerously concentrated in the pillar. At 60 percent, the backfill began sharing roof load, though stress zones remained isolated. Only at 80 percent of coal strength did the stresses within the pillar and backfill merge into a single continuous high-strength bearing zone without peak eccentricity. The backfill also exerts lateral confinement on the flanking coal segments, creating what the authors call a strong-weak-strong configuration that boosts the pillar’s effective bearing capacity by approximately 27.9 percent when a single roadway among three is filled.

Field deployment at the 110203 longwall face validated the approach under production conditions. Borehole peeping through the pillar rib in the abandoned roadway concentration zone showed high internal integrity after treatment, and convergence monitoring recorded roughly 16 centimeters of rib-to-rib closure and 14 centimeters of roof-to-floor deformation at backfilled locations, figures notably smaller than in untreated areas. The study’s broader message for the industry is twofold: roadway density in protective pillars should be carefully controlled during mine planning, and where legacy roadways already compromise pillar integrity, solidified excavated coal offers a low-cost, locally sourced route to restoring load-bearing capacity and keeping goaf-side entries safe for mining.

Subject of Research: Instability of goaf-side narrow coal pillars beneath concentrated abandoned roadways in extra-thick coal seams and solidification-backfilling control

Article Title: Instability mechanism of narrow coal pillars in extra-thick coal seam within concentrated abandoned roadway areas and solidification-backfilling control of roadway excavation coal mass

Article References: Instability mechanism of narrow coal pillars in extra-thick coal seam within concentrated abandoned roadway areas and solidification-backfilling control of roadway excavation coal mass. (n.d.). https://doi.org/10.1016/j.rineng.2026.112978

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.112978

Keywords: coal mining, narrow coal pillar, abandoned roadways, goaf-side entry, roof stability, triangular block, numerical simulation, backfilling, Portland cement, strata control, FLAC3D, mine safety

The core of the work, titled Numerical Simulation of High-Temperature Thermochemical Non-equilibrium Flows Under Aerodynamic Heating, is a computational framework built on a conformal mesh node approach. In conventional hypersonic simulations, the fluid dynamics, the gas chemistry, and the thermal response of the vehicle structure are often computed separately and then loosely linked, if they are linked at all. Each hand-off between models introduces error, and errors compound precisely where engineers care most: at the vehicle surface, where heat flux, wall temperature, and near-wall chemistry jointly determine whether a thermal protection system survives re-entry or mission flight. The new method instead solves the wall aerodynamic heating, a two-temperature model of the gas, and finite-rate chemical reactions as one interacting system, so that heat entering the structure and the changing wall temperature feed directly back into the near-wall flow field at every step of the calculation.

The two-temperature model at the heart of the formulation reflects a crucial piece of high-enthalpy physics. Behind a strong shock wave, the translational and rotational modes of air molecules adjust almost instantly, while the vibrational modes lag behind, and chemical dissociation proceeds on its own finite time scale. Rather than forcing a single thermodynamic temperature onto the gas, the model carries a separate temperature for the translational-rotational degrees of freedom and another for vibrational excitation, allowing energy exchange between them through relaxation terms of the kind introduced in the classic Landau-Teller framework. Chemical source terms for the dissociation and exchange reactions of an eleven-species air model, with reaction rates of the Park type, are coupled to these temperatures, so the chemistry and the thermal nonequilibrium evolve together rather than in sequence.

Transport properties, often treated as an afterthought, are handled with equal care. Mixture viscosities follow established mixing rules, diffusion of species is represented through formulations rooted in Fick’s law, and the reaction-rate and thermodynamic data draw on widely used NASA reference compilations for high-temperature air. These choices matter because near a hot wall the composition of the gas changes rapidly: molecular oxygen and nitrogen dissociate, atoms accumulate, and the mixture’s viscosity, conductivity, and diffusivity all shift. A simulation that freezes these properties at freestream values can mispredict wall heat flux by a meaningful margin, and it is wall heat flux that sizes the thermal protection system.

To test whether the coupled approach actually improves fidelity, the team first applied it to a deceptively simple geometry: laminar flow over a circular cylinder, one of the canonical validation cases in hypersonic aerothermodynamics, with experimental shock-layer data available from high-enthalpy ground tests. The quantity of greatest diagnostic value here is the shock stand-off distance, the gap between the bow shock and the body surface. That distance is governed by the density rise across the shock layer, which in turn is controlled by real-gas effects; if the simulation gets the thermochemistry wrong, the shock sits in the wrong place. The results showed that the shock stand-off distance computed with the fully coupled method agrees better with experimental data than predictions that ignore the structural heat transfer. The computed wall friction coefficient, meanwhile, increased slightly relative to solutions that neglect heat conduction into the solid, a sign that cooling of the near-wall gas through the wall thickens the boundary layer’s influence on the surface shear in ways that uncoupled models miss.

That modest increase in skin friction is more than a numerical curiosity. It signals that the energy exchanged between the fluid and the structure is large enough to reshape the flow itself, and any design tool that pretends the wall is adiabatic or held at a fixed temperature will inherit that blind spot. For a hypersonic vehicle flying a long trajectory, wall temperature rises over minutes, not seconds, and the evolving thermal state of the structure continuously modifies the chemistry and heat flux at the surface. A coupled method of the kind developed here allows engineers to simulate that feedback loop rather than assuming it away, which is precisely what the authors identify as the key challenge in the thermal protection design of hypersonic vehicles.

Encouraged by the cylinder validation, the researchers then scaled up to a case with genuine engineering relevance: a hypersonic wing flying under aerodynamic heating conditions. Here the emphasis fell on two quantities. The first was again the shock stand-off distance along the leading edge, where the sweep and curvature of a real wing produce shock layers that vary spanwise in ways a cylinder never can. The second was the near-wall chemically non-equilibrium flow field, the thin region where dissociated atoms recombine, vibrational temperatures relax toward translational values, and species gradients are steepest. According to the study, the results demonstrate that the coupled method can reasonably characterize how aerodynamic heating influences these near-wall nonequilibrium characteristics, giving designers a tool that links the structural thermal answer and the fluid thermochemical answer in a single, consistent solution rather than two partially reconciled ones.

The implications reach across the current wave of hypersonic development. Reusable launch systems, glide vehicles, and planetary entry capsules all spend critical portions of their trajectories in exactly the regime this method targets, where flight enthalpies are too high for the gas to behave as a calorically perfect ideal. Ground-test facilities can reproduce some of these conditions, but rarely all of them at once; vibrationally cold but chemically energetic flows in one facility, clean equilibrium flows in another. High-fidelity simulation that honestly couples the structure to the flow offers a way to bridge the gaps between sparse test data, and the authors note that the conformal mesh node formulation is what makes the fluid-solid coupling seamless at the shared boundary, avoiding the interpolation losses that plague loosely coupled schemes.

The team is explicit that the current framework is a foundation rather than a finished product. Future work, they write, may incorporate additional physical fields such as turbulence, radiation heat transfer, and wall catalytic effects. Each addition addresses a known gap: turbulence alters heat transfer dramatically along real vehicle surfaces, radiative heating becomes significant at entry speeds where shock layers glow, and wall catalycity, the tendency of a surface to promote recombination of dissociated atoms, can dump substantial additional energy into the wall. Turbulence-chemistry-radiation interactions coupled through a structural thermal solver represent one of the remaining grand challenges in hypersonics, and the architecture described in this study, funded by the National Natural Science Foundation of China under Grant No. 52275143, provides a credible platform on which those effects can be layered. For now, the message for the field is concrete: when the gas outside a hypersonic vehicle is hotter than the surface of a star’s atmosphere, the wall and the flow must be solved as one problem, and this work shows a validated way to do it.

Subject of Research: Fully coupled numerical simulation of high-temperature thermochemical non-equilibrium flows under aerodynamic heating for hypersonic vehicle thermal protection

Article Title: Numerical Simulation of High-Temperature Thermochemical Non-equilibrium Flows Under Aerodynamic Heating

Article References: Numerical Simulation of High-Temperature Thermochemical Non-equilibrium Flows Under Aerodynamic Heating. (n.d.). https://doi.org/10.1007/s42405-026-01256-x

Image Credits: AI Generated

DOI: 10.1007/s42405-026-01256-x

Keywords: hypersonic, aerodynamic heating, thermochemical non-equilibrium flow, numerical simulation, shock wave, boundary layer, two-temperature model, shock stand-off distance, thermal protection, computational fluid dynamics, high-enthalpy flow, chemical reactions

The research, led by Alejandro González and Jorge Macías of the University of Málaga together with Cléa Denamiel of the Ruđer Bošković Institute in Croatia, evaluates the performance of the Meteo-HySEA model against the state-of-the-art AdriSC-ADCIRC modeling suite. Meteo-HySEA, developed by the EDANYA research group in Málaga, belongs to the HySEA family of codes, which also includes the widely used Tsunami-HySEA and Landslide-HySEA models recognized as European flagship tools for simulating seismically and landslide-triggered tsunamis. What sets Meteo-HySEA apart is its combination of high-resolution finite volume numerical schemes with a native multi-GPU framework, allowing it to incorporate time-dependent atmospheric pressure fields as forcing while running orders of magnitude faster than conventional CPU-based models. Crucially, it can also simulate the onshore inundation that follows, extending the modeling chain from offshore wave generation all the way to flooded harbor quays.

The benchmark for the new model was the AdriSC system, a CPU-based coupled atmosphere-ocean modeling framework that has become the reference tool for meteotsunami research in the Adriatic. AdriSC couples the Weather Research and Forecasting model, downscaled to 1.5-kilometer resolution over the Adriatic, with the two-dimensional depth-integrated ADCIRC ocean model running on an unstructured mesh refined to spatial resolutions of up to 10 meters in vulnerable zones. This minute-scale atmospheric forcing is essential to capture the speed and amplitude of tsunamigenic pressure disturbances that excite Proudman resonance across the basin and the fundamental oscillation modes of the bays, which typically fall in the 10-to-40-minute period band. However, the computational cost of such CPU-based frameworks creates serious bottlenecks for ensemble forecasting, probabilistic hazard assessment, and real-time early warning, where reducing simulation times from hours to minutes is essential.

The team focused on three well-documented events that struck the meteotsunami-prone harbors of Vela Luka on Korčula Island and Stari Grad and Vrboska on Hvar Island. The first, on 25–26 June 2014, was triggered by a train of atmospheric gravity waves propagating from the Tyrrhenian Sea across the Adriatic, producing rapid pressure perturbations of up to 2.4 hectopascals in five minutes and maximum sea levels reaching 1.5 meters inside the harbors. The second, from late June to early July 2017, was associated with a synoptic cyclone over the central Mediterranean and upper-level jet stream winds exceeding 55 meters per second; oscillations lasted nearly 24 hours, with amplitudes up to 0.69 meters at tide gauges and exceeding one meter as captured in videos from Vrboska. The third, a remarkable multi-day sequence between 11 and 19 May 2020, produced repeated waves of 0.6 to 0.8 meters driven by recurring high-frequency pressure disturbances of 2 to 4 hectopascals, flooding harbors and leaving boats ashore on multiple days.

A distinctive feature of the study is its deliberately strict evaluation protocol. No correction or tuning of the atmospheric forcing was applied, ensuring that both modeling systems were tested under realistic operational conditions rather than with the benefit of hindsight. This choice matters because retrospective adjustments of pressure fields could improve agreement with observations, but such approaches rely on prior knowledge of the event and are therefore useless in a real-time forecasting context. The evaluation drew on high-frequency observations from the MESSI observational network, which includes microbarographs measuring atmospheric pressure at one-minute intervals with a precision of plus or minus 0.01 hectopascals and radar tide gauges recording sea level with millimeter accuracy, supplemented for the earlier events by five-minute pressure records from the Crometeo network of amateur weather stations across the Croatian coast.

The results reveal a fundamental constraint that applies to any meteotsunami modeling system: simulation accuracy is ultimately limited by the quality of the atmospheric forcing. Comparing weather model simulations driven by the ERA-Interim and ERA5 reanalyses, the researchers found that the ERA-Interim forcing systematically generated more intense and widespread pressure anomalies, sometimes exceeding observed amplitudes and producing spurious fluctuations, while ERA5-driven simulations yielded more localized and lower-amplitude disturbances that aligned better with observations in timing and structure but often missed the strongest pressure jumps. Most strikingly, the ERA5-driven simulation completely failed to reproduce the sharp 2.5-hectopascal pressure jump recorded on 11 May 2020, the most intense disturbance of the entire study period. Because observational coverage over the central Adriatic remains sparse, the authors note, fully assessing the skill of atmospheric models in this region is itself a difficult task.

Against this backdrop, Meteo-HySEA performed encouragingly. The model successfully reproduced the timing and spatial variability of sea-level oscillations across all three events and generally yielded higher amplitudes than AdriSC-ADCIRC under the same forcing. During the June 2014 event, for example, Meteo-HySEA produced peak amplitudes of roughly one meter in Vela Luka compared with 0.6 meters from AdriSC-ADCIRC, suggesting it may better capture the potential extremes of meteotsunami events. However, a systematic shortcoming also emerged: Meteo-HySEA consistently overestimated the dominant wave periods, particularly in semi-enclosed basins. In Stari Grad, its median simulated wave periods exceeded observed values by up to 50 percent. The authors attribute this to the wet-dry technique used to simulate inundation, which causes the effective geomorphology of harbors to evolve over time, and to limited coastal bathymetric coverage. Both models, it should be noted, also overestimated observed periods, indicating that harbor resonance characteristics remain difficult to represent accurately in general.

To probe the differences between the two systems more deeply, the team ran controlled numerical experiments with synthetic pressure disturbances based on the analytical formulation developed for the AdriSC meteotsunami surrogate model. These idealized disturbances, defined by parameters such as origin, propagation direction, amplitude, translation speed, period, and spatial width, were chosen to produce the strongest plausible meteotsunamis in each harbor. The comparison exposed systematic differences in how the two models trap and dissipate energy within semi-enclosed basins. In Vela Luka, Meteo-HySEA predicted extreme elevations exceeding five meters under the most severe synthetic forcing, compared with peaks closer to four meters for AdriSC-ADCIRC, along with a slower decay of oscillations. The analysis also uncovered a genuine physical phenomenon: a massive, localized meteotsunami setup of up to 1.8 meters at the tip of Vela Luka harbor, generated by resonant mass transport, on which the high-frequency seiches ride. Both components, the authors stress, must be retained together to assess total coastal flooding hazard.

The study also candidly documents technical limitations. Spurious hotspots appeared in maximum sea-level maps for Stari Grad and Vrboska, driven by sparse raw bathymetric data near coastlines and numerical artifacts of projecting complex coastal boundaries onto high-resolution structured grids, where abrupt bathymetric steps cause the flow to behave as if hitting a vertical wall. In Vrboska, the narrow channel characteristics could not be resolved even at seven-meter grid spacing, shifting the modeled coastline relative to satellite imagery. The authors emphasize that numerical interpolation cannot artificially reconstruct missing physical topography, making high-quality, high-resolution topobathymetric surveys a critical requirement for nearshore areas if the inundation capabilities of the new model are to be fully exploited.

Nevertheless, the overall verdict is that GPU-based solvers like Meteo-HySEA represent a promising pathway toward next-generation meteotsunami forecasting and hazard assessment. Because the model runs orders of magnitude faster than CPU-based alternatives, it opens the door to real-time operational use, ensemble forecasting, and explicit treatment of atmospheric and boundary-condition uncertainties, all of which are essential given that the lead time between disturbance detection and coastal impact in the Adriatic is typically less than three hours. The authors outline future work on three fronts: implementing wind stress, tidal, and general circulation forcing in the next release of the code; systematic validation with dense observational networks to determine whether the longer persistence of oscillations in Meteo-HySEA reflects more realistic harbor seiches or whether ADCIRC’s stronger damping better represents physical energy dissipation; and testing the model’s operational potential through integration with real-time atmospheric forecasts and early warning protocols in collaboration with civil protection agencies. With atmospheric forcing remaining the dominant source of forecast uncertainty, advances in convection-permitting ensembles and machine-learning-based nowcasting are expected to further improve fidelity, strengthening coastal resilience in the Adriatic and potentially at meteotsunami hotspots worldwide.

Subject of Research: Evaluation of the GPU-based Meteo-HySEA model for simulating atmospherically driven meteotsunami events and harbor flooding in the Adriatic Sea

Article Title: Assessing Meteo-HySEA performance for Adriatic meteotsunami events

Article References: Assessing Meteo-HySEA performance for Adriatic meteotsunami events. (n.d.). https://doi.org/10.1007/s11069-026-08351-y

Image Credits: AI Generated

DOI: 10.1007/s11069-026-08351-y

Keywords: meteotsunami, Adriatic Sea, Meteo-HySEA, GPU modeling, harbor resonance, Proudman resonance, coastal flooding, early warning systems, atmospheric pressure disturbances, AdriSC-ADCIRC, sea-level oscillations, numerical simulation

Traditional mineral prospectivity modeling has long relied on what geologists call post-mineralization data: patterns recorded in rocks, soils, and geophysics long after the ore-forming event ended. These data-driven approaches treat mineralization as a static snapshot, effectively asking where gold is known to occur and searching for similar patterns elsewhere. The problem, the researchers argue, is that this strategy overlooks the dynamic controls that operated during the metallogenic evolution itself, the shifting stresses, migrating fluids, and thermal gradients that determined where gold was precipitated in the first place. Where a deposit sits today is the end product of a long, physically coupled process, and the fingerprints of that process are often subtle, deeply buried, and invisible to conventional exploration datasets. This limitation becomes especially severe in the search for concealed ore bodies at depth, where surface anomalies fade and deposit models extrapolated from shallow workings begin to lose their predictive power.

To overcome this, the team designed a three-stage workflow that moves from static geometry to dynamic process to intelligent integration. The first stage uses 3D spatial analysis to quantify the morphological features of the ore-controlling faults, the geological structures that channeled mineralizing fluids, and the primary geochemical halos, the chemical dispersal zones that surround ore bodies. Rather than simply drawing buffers around faults, the method extracts quantitative shape descriptors from the three-dimensional geometry of these features, capturing how fault orientations, curvature, and intersections create favorable sites for fluid focusing and gold deposition. This converts qualitative geological intuition, the sense that a fault bend or junction might be favorable, into explicit numerical predictor layers that a machine learning model can digest.

The second stage is the scientific heart of the approach: a coupled mechanical-thermal-hydrological, or MTH, numerical simulation of the ore-forming process itself. By building a three-dimensional computational model of the deposit’s structural framework and assigning rock properties drawn from established geomechanics and hydrogeology references, the researchers simulated how tectonic stresses deformed the rock mass, how heat redistributed through the system, and how hydrothermal fluids were driven through the permeable fault networks. Crucially, this simulation yields quantities that no drill core or geochemical survey can directly measure: the evolution of tectonic strain through time and the spatial distribution of fluid flux, two of the implicit geodynamic predictors that control whether dissolved gold is carried, concentrated, or dropped from solution. Where deformation localizes and fluids converge, epithermal gold systems like Guilaizhuang tend to deposit their metal, and the simulation pinpoints those zones in three dimensions.

The researchers describe the simulation output as effectively characterizing the spatiotemporal evolution of the metallogenic process at Guilaizhuang, providing crucial physical constraints that conventional prospectivity modeling lacks. Instead of inferring favorable conditions purely from where known ore is found, the model can identify where the physics of the system says ore formation was most likely, including in deep and lateral regions that have never been drilled. This coupling of process simulation with exploration targeting reflects a growing trend in computational geoscience, in which numerical experiments on coupled deformation, fluid flow, and heat transport serve as virtual laboratories for reconstructing mineral systems that humans can never observe directly.

With static geological predictors and dynamic simulation outputs in hand, the team faced a final challenge: how to fuse these multi-source, heterogeneous layers into a single coherent prospectivity map. Their answer is a purpose-built deep learning architecture called 3D CBAM-ResCNN, an attention-enhanced three-dimensional convolutional neural network that combines residual structures with the convolutional block attention module, or CBAM. Residual connections, popularized in computer vision, allow very deep networks to train stably by letting information bypass layers, while CBAM teaches the network to selectively emphasize the most informative channels and spatial locations in the data. In practical terms, the attention mechanism lets the model decide, voxel by voxel and feature by feature, which predictors genuinely matter for gold mineralization and which are redundant noise, a critical capability when combining dozens of overlapping geological, geochemical, and geodynamic layers.

The results show that the 3D CBAM-ResCNN achieves the best performance among the configurations tested, excelling at identifying the spatial dependencies that link mineralization to its controlling features while suppressing the feature redundancy that degrades simpler models. Standard three-dimensional convolutional networks without attention tend to treat all input layers equally, allowing noisy or correlated predictors to dilute the signal; the attention-enhanced architecture concentrates its learning capacity on the fault morphology descriptors and simulation-derived strain and fluid flux fields that carry the real predictive weight. The prospectivity volumes the network produces score highest in accuracy and reliability, correctly reproducing the spatial distribution of known mineralization while extending meaningful predictions into unexplored territory.

Perhaps the most consequential output for explorers is the delineation of two exploration targets, zones where the model’s probability estimates rise sharply despite lying beyond the currently well-understood footprint of the deposit. These targets provide a scientific basis for future deep drilling at Guilaizhuang, offering the kind of quantitative, physically grounded justification that exploration managers need before committing expensive drill campaigns. Given that the Guilaizhuang system is recognized as having significant deep-seated metallogenic potential, finding the next ore body at depth could meaningfully extend the life and economics of the mining district. The study was supported by China’s National Science and Technology Major Project, the National Natural Science Foundation of China, and the Key Research and Development Plan of Shandong Province, with exploration data supplied by the Shandong Provincial Lunan Geology and Exploration Institute.

Beyond one gold deposit in Shandong, the framework points toward a broader transformation in how hidden mineral resources are found worldwide. As shallow discoveries become rarer and exploration moves deeper, the industry increasingly needs methods that combine mechanistic understanding with machine learning rather than relying on correlation alone. The Guilaizhuang study demonstrates that simulated strain evolution and fluid flux can serve as first-class predictors alongside conventional geology and geochemistry, and that attention-based 3D neural networks can orchestrate this diverse evidence with measurable gains in accuracy. For a discipline racing to supply the metals of the energy transition, the message is clear: the fastest route to buried treasure may run through supercomputers first, with drills following the physics where the algorithms say to look. The datasets generated in the study are not publicly available due to a confidentiality agreement, but the published methodology offers a replicable blueprint for three-dimensional targeting wherever structurally controlled hydrothermal systems remain hidden in the deep subsurface.

Guilaizhuang belongs to a distinctive family of gold deposits in which gold occurs with telluride minerals, and earlier studies of the Pingyi area have documented telluride-bearing Au mineralization linked to fluid boiling, a process that can trigger rapid gold precipitation when pressure drops in rising hydrothermal fluids. That geological character helps explain why the fault-focused fluid pathways reconstructed by the MTH simulation carry such predictive weight: boiling and fluid mixing in epithermal systems are tightly controlled by where deformation localizes and where flow converges.

The deposit also sits on the southeastern margin of the North China Craton, a region whose lithospheric thinning and repeated magmatic pulses have long been linked to gold metallogeny in western Shandong. Pyrite chemistry and in situ sulfur isotope work on Guilaizhuang ores has further constrained gold enrichment mechanisms, giving later modelers a well-studied natural laboratory. Against that backdrop, coupling process simulation with attention-based deep learning offers a way to translate decades of deposit-scale research into quantitative, three-dimensional exploration guidance.

Subject of Research: Three-dimensional gold prospectivity modeling using coupled numerical simulation and attention-enhanced deep learning at the Guilaizhuang deposit, China

Article Title: Combining 3D Numerical Simulation and Attention-Enhanced 3D CNN for Mineral Prospectivity Modeling: A Case Study of the Guilaizhuang Gold Deposit, Western Shandong, China

Article References: Zou, Y., Chen, G., Li, J., & Mao, X. (2026). Combining 3D Numerical Simulation and Attention-Enhanced 3D CNN for Mineral Prospectivity Modeling: A Case Study of the Guilaizhuang Gold Deposit, Western Shandong, China. Natural Resources Research. https://doi.org/10.1007/s11053-026-10768-y

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10768-y

Keywords: 3D mineral prospectivity modeling, numerical simulation, 3D CBAM-ResCNN, gold exploration, Guilaizhuang gold deposit, epithermal gold, fluid flux, tectonic strain evolution, deep learning, 3D geological modeling, exploration targeting, Western Shandong