arxiv: 2604.14196 · v1 · submitted 2026-04-02 · ⚛️ physics.geo-ph

Recognition: 2 theorem links

· Lean Theorem

Toward buoyancy-driven flow at Campi Flegrei: coupled phase change and asymmetric geometry

Arturo Tozzi

Pith reviewed 2026-05-13 20:47 UTC · model grok-4.3

classification ⚛️ physics.geo-ph

keywords buoyancy-driven flowphase changeasymmetric geometryCampi FlegreiDarcy flowvolcanic unresthydrothermal pressurizationfluid dynamics

0 comments

The pith

Buoyancy from phase change and asymmetric geometry drives localized upward flow and uplift at Campi Flegrei.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models fluid dynamics in a shallow gas-rich reservoir where phase transitions, density gradients, and geometric asymmetry interact to generate buoyancy-driven flow along preferential pathways. This setup draws an analogy to self-propulsion in melting asymmetric ice blocks, with simulations incorporating time-dependent pressure forcing and hydraulic relaxation. The runs, calibrated to observed deformation rates and seismicity depths, reproduce temporal uplift variations and the persistence of spatially localized flow. The model yields predictions of nonlinear uplift acceleration, shallow localized seismicity, and flow velocities that scale with pressure and buoyancy forces.

Core claim

Coupling phase change with asymmetric geometry in a heterogeneous subsurface promotes channelized upward transport while phase transitions enhance buoyancy and redistribute pressure. Numerical simulations of buoyancy-enhanced Darcy flow in a gas-rich reservoir, constrained by Campi Flegrei deformation and seismicity data, reproduce observed temporal variations in uplift together with the persistence of spatially localized flow. Within this framework the same processes are shown to produce nonlinear uplift acceleration, shallow localized seismicity, and velocity scaling with pressure and buoyancy.

What carries the argument

buoyancy-enhanced Darcy flow along prescribed preferential pathways, coupled to phase transition and geometric asymmetry

If this is right

Asymmetric geometry channels upward fluid transport and sustains localized flow.
Phase change enhances buoyancy and contributes to pressure redistribution.
Uplift accelerates nonlinearly under continued time-dependent forcing.
Seismicity remains localized at shallow depths consistent with reported observations.
Flow velocity increases proportionally with pressure and buoyancy.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same buoyancy mechanism could be tested in other caldera systems that exhibit asymmetric deformation patterns.
Integration with multiphase reservoir models would allow quantitative assessment of how buoyancy alters long-term pressure evolution.
Time-series monitoring of localized flow rates versus surface deformation could provide a direct check on the predicted velocity scaling.
If the mechanism holds, asymmetric structural features mapped by tomography would become priority targets for unrest forecasting.

Load-bearing premise

Phase change and density gradients in subsurface volcanic fluids produce net flow in a manner directly analogous to the self-propulsion of melting asymmetric ice blocks.

What would settle it

Measurements showing no scaling of fluid velocity with buoyancy or pressure, or seismic tomography revealing uniform rather than channelized flow, would falsify the central mechanism.

read the original abstract

Bradyseism at Campi Flegrei is usually interpreted in terms of hydrothermal pressurization and magmatic degassing. Fluid flow, often treated as a passive response to pressure accumulation, is commonly modeled using simplified geometries and homogeneous permeability fields. We introduce a model in which phase transition, structural heterogeneity and geometric asymmetry jointly influence fluid flow and pressure distribution within a heterogeneous subsurface environment. We hypothesize that coupling among phase change, density gradients and flows may follow a mechanism similar to the self-propulsion observed in asymmetric floating bodies like melting ice blocks, where phase change generates buoyancy-driven currents along their inclined surfaces and net motion in the opposite direction. We simulate pressure evolution in a shallow gas-rich reservoir subject to time-dependent forcing and hydraulic relaxation, coupled to buoyancy-enhanced Darcy flow along prescribed preferential pathways. Our numerical simulations, grounded in reported deformation rates and seismicity depths at Campi Flegrei, reproduce temporal variations in uplift and the persistence of spatially localized flow. Within this framework, asymmetric geometry may promote channelized upward transport, while phase change may enhance buoyancy and contribute to pressure redistribution. Our model predicts nonlinear uplift acceleration, shallow localized seismicity and velocity scaling with pressure and buoyancy. Integration with existing multiphase models would enable the examination of how buoyancy-driven flows influence pressure evolution and deformation during volcanic unrest.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The buoyancy analogy from ice blocks is a reasonable extension for Campi Flegrei flows, but localization stays tied to prescribed pathways rather than emerging cleanly from the new coupling.

read the letter

The paper's main move is to adapt the self-propulsion seen in melting asymmetric ice blocks to a volcanic setting, where phase change and density gradients could drive buoyancy-enhanced Darcy flow along asymmetric paths in a shallow reservoir. It couples this to time-dependent forcing and hydraulic relaxation, then runs simulations tuned to Campi Flegrei uplift rates and seismicity depths. The results match observed temporal uplift variations and keep flow localized, which is consistent with field reports. The suggestion that geometric asymmetry channels upward transport while phase change boosts buoyancy gives a concrete way to link multiphase effects to deformation patterns. That part is useful for people already working on hydrothermal pressurization models. The central limitation is that the preferential pathways are set in advance from seismicity data. This means the persistence of localized flow is not generated by the phase-change and buoyancy mechanism alone; it is helped by the fixed channels in the setup. Without sensitivity tests that remove or randomize those paths, it is hard to tell how much the ice-block analogy actually contributes versus the boundary conditions. The abstract also gives no equations, error bars, or independent validation checks, so the reproduction of uplift could partly reflect parameter choices. This work is aimed at volcanologists and fluid-dynamics modelers who study caldera unrest and want to add buoyancy considerations to existing Darcy and multiphase frameworks. It is worth a serious referee because the hypothesis is testable and relevant to hazard models, even though the current implementation needs clearer separation between imposed and emergent behavior.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes a numerical model for fluid dynamics at Campi Flegrei that couples phase transitions, structural heterogeneity, and geometric asymmetry. Drawing an analogy to self-propulsion in melting asymmetric ice blocks, it simulates pressure evolution in a shallow gas-rich reservoir under time-dependent forcing and hydraulic relaxation, using buoyancy-enhanced Darcy flow along prescribed preferential pathways. The simulations are claimed to reproduce observed temporal variations in uplift and the persistence of spatially localized flow, while predicting nonlinear uplift acceleration, shallow localized seismicity, and velocity scaling with pressure and buoyancy.

Significance. If the central mechanism could be shown to operate without imposed pathways and with independent validation, the work would provide a useful extension to existing hydrothermal models by highlighting how buoyancy and asymmetry might channelize flow and redistribute pressure in volcanic systems, potentially informing interpretations of deformation and seismicity during unrest.

major comments (2)

[Abstract / Model Description] Abstract and model setup: the preferential pathways for buoyancy-enhanced Darcy flow are prescribed a priori from reported seismicity depths rather than emerging from the hypothesized phase-change/density/asymmetry coupling. Consequently the reported persistence of localized flow and reproduction of uplift variations are built into the boundary conditions and forcing, not generated by the ice-block analogy; this directly undermines the claim that the simulations demonstrate the proposed mechanism.
[Abstract / Results] Abstract / Results: the manuscript states that simulations reproduce uplift variations and predict nonlinear acceleration but supplies no governing equations, constitutive relations for phase change, quantitative error metrics, or comparison against independent benchmarks. Without these, it is impossible to assess whether the match to reported deformation rates arises from the physics or from tuning of the time-dependent forcing and relaxation parameters.

minor comments (1)

[Methods] Notation for permeability heterogeneity and buoyancy terms should be defined explicitly in the methods section with reference to the Darcy formulation used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive review and the opportunity to address these points. We respond to each major comment below, indicating planned revisions where appropriate.

read point-by-point responses

Referee: [Abstract / Model Description] Abstract and model setup: the preferential pathways for buoyancy-enhanced Darcy flow are prescribed a priori from reported seismicity depths rather than emerging from the hypothesized phase-change/density/asymmetry coupling. Consequently the reported persistence of localized flow and reproduction of uplift variations are built into the boundary conditions and forcing, not generated by the ice-block analogy; this directly undermines the claim that the simulations demonstrate the proposed mechanism.

Authors: We agree that the preferential pathways are prescribed a priori using reported seismicity depths, as explicitly stated in the model description. This choice incorporates direct observational constraints from Campi Flegrei rather than allowing pathways to emerge self-consistently from the phase-change and asymmetry coupling in the present implementation. The ice-block analogy is invoked to motivate buoyancy enhancement along inclined surfaces within those pathways, but we acknowledge that the current setup does not demonstrate spontaneous channelization. In the revised manuscript we will clarify this distinction in the abstract, methods, and discussion, stating that the model illustrates the effects of buoyancy-driven flow along observationally informed pathways and will outline future work on self-consistent pathway emergence. revision: partial
Referee: [Abstract / Results] Abstract / Results: the manuscript states that simulations reproduce uplift variations and predict nonlinear acceleration but supplies no governing equations, constitutive relations for phase change, quantitative error metrics, or comparison against independent benchmarks. Without these, it is impossible to assess whether the match to reported deformation rates arises from the physics or from tuning of the time-dependent forcing and relaxation parameters.

Authors: The full manuscript contains the governing Darcy-flow equations, the constitutive relations for phase change (gas exsolution and associated density variations), and the numerical discretization. However, these elements were not presented with sufficient detail or quantitative validation in the abstract and results summary. In the revision we will add an explicit methods subsection with the complete equation set, include quantitative error metrics (e.g., normalized RMS misfit to observed uplift time series), and provide benchmark comparisons against simplified homogeneous-reservoir cases to demonstrate that the reported behavior is not solely the result of parameter tuning. revision: yes

Circularity Check

1 steps flagged

Prescribed preferential pathways and data-tuned forcing make localized flow and uplift 'reproduction' by construction

specific steps

fitted input called prediction [Abstract]
"We simulate pressure evolution in a shallow gas-rich reservoir subject to time-dependent forcing and hydraulic relaxation, coupled to buoyancy-enhanced Darcy flow along prescribed preferential pathways. Our numerical simulations, grounded in reported deformation rates and seismicity depths at Campi Flegrei, reproduce temporal variations in uplift and the persistence of spatially localized flow."

Preferential pathways are prescribed a priori (fixed locations from seismicity data) and forcing parameters are tuned to reported rates; the 'reproduction' of localized flow persistence and uplift variations is therefore imposed by the setup rather than generated by the hypothesized coupling of phase change, density gradients and asymmetry.

full rationale

The paper's derivation chain begins with a hypothesized analogy to ice-block self-propulsion but immediately implements buoyancy-enhanced Darcy flow only along explicitly prescribed preferential pathways whose locations are fixed from reported seismicity depths. Time-dependent forcing and hydraulic relaxation are likewise grounded in observed deformation rates. The claimed reproduction of temporally varying uplift and persistence of spatially localized flow therefore follows directly from these inputs rather than emerging from phase-change or geometric-asymmetry dynamics. The subsequent 'predictions' of nonlinear acceleration and velocity scaling inherit the same tuning. This matches the fitted-input-called-prediction pattern with no independent, unfitted benchmark shown.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The model rests on standard geophysical assumptions for Darcy flow and phase transitions without introducing new entities or fitted constants beyond those implied by matching reported deformation rates.

free parameters (1)

heterogeneous permeability fields
Prescribed preferential pathways and time-dependent forcing are set to match observed rates but not independently derived.

axioms (2)

standard math Darcy's law governs fluid flow in porous media
Invoked for buoyancy-enhanced flow along pathways in the shallow reservoir.
domain assumption Phase change generates density gradients that drive buoyancy
Core to the ice-block analogy applied to gas-rich reservoir.

pith-pipeline@v0.9.0 · 5526 in / 1262 out tokens · 45425 ms · 2026-05-13T20:47:00.754194+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/ArithmeticFromLogic.lean reality_from_one_distinction unclear
coupled to buoyancy-enhanced Darcy flow along prescribed preferential pathways... x_c(z) = -0.6 + 0.45(5-z) - 0.06(5-z)^2
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
u(t) = k/μ (P(t)×10^6/L + Δρ g)

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

[1]

Acocella, Valerio. 2021. Volcano-Tectonic Processes . Advances in V olcanology. Cham: Springer. https://doi.org/10.1007/978-3-030-65968-4

work page doi:10.1007/978-3-030-65968-4 2021
[2]

Tracking the 2007–2023 Magma-Driven Unrest at Campi Flegrei Caldera (Italy)

Astort, A., E. Trasatti, L. Caricchi, et al. 2024. “Tracking the 2007–2023 Magma-Driven Unrest at Campi Flegrei Caldera (Italy).” Communications Earth & Environment 5: 506. https://doi.org/10.1038/s43247-024-01665-4

work page doi:10.1038/s43247-024-01665-4 2024
[3]

Long-Term Thermo-Fluid Dynamic Modeling of Solfatara Hydrothermal System, Campi Flegrei Caldera

Barone, Andrea, Gianluca Gola, Stefano Caliro, Giovanni Chiodini, Pietro Tizzani, and Raffaele Castaldo. 2025. “Long-Term Thermo-Fluid Dynamic Modeling of Solfatara Hydrothermal System, Campi Flegrei Caldera.” Journal of Volcanology and Geothermal Research 459: 108277. https://doi.org/10.1016/j.jvolgeores.2025.108277

work page doi:10.1016/j.jvolgeores.2025.108277 2025
[4]

Self -Propulsion of Floating Ice Blocks Caused by Melting in Water

Berhanu, Michael, Amit Dawadi, Martin Chaigne, Jérôme Jovet, and Arshad Kudrolli. 2026. “Self -Propulsion of Floating Ice Blocks Caused by Melting in Water.” Physical Review Fluids 11: 033802. https://doi.org/10.1103/k6gl-1zqm

work page doi:10.1103/k6gl-1zqm 2026
[5]

Seismological Evidence for a Multifault Network at the Subduction Interface

Chalumeau, C., H. Agurto-Detzel, A. Rietbrock, et al. 2024. “Seismological Evidence for a Multifault Network at the Subduction Interface.” Nature 628: 558–562. https://doi.org/10.1038/s41586-024-07245-y

work page doi:10.1038/s41586-024-07245-y 2024
[6]

Evolution in Unrest Processes at Campi Flegrei Caldera as Inferred from Seismicity and Ground Deformation

Danesi, Sara, Christopher R. J. Kilburn, Giuseppe De Natale, et al. 2024. “Evolution in Unrest Processes at Campi Flegrei Caldera as Inferred from Seismicity and Ground Deformation.” Earth and Planetary Science Letters 630: 118600

work page 2024
[7]

3D Structure and Dynamics of Campi Flegrei Enhance Multi-Hazard Assessment

De Landro, G., T. Vanorio, T. Muzellec, et al. 2025. “3D Structure and Dynamics of Campi Flegrei Enhance Multi-Hazard Assessment.” Nature Communications 16: 4814. https://doi.org/10.1038/s41467-025-59821-z

work page doi:10.1038/s41467-025-59821-z 2025
[8]

Reconstructing Hydrothermal Fluid Pathways and Storage at the Solfatara Crater (Campi Flegrei, Italy) Using Seismic Scattering and Absorption

Di Martino, Maria Del Pilar, Luca De Siena, Vincenzo Serlenga, and Grazia De Landro. 2022. “Reconstructing Hydrothermal Fluid Pathways and Storage at the Solfatara Crater (Campi Flegrei, Italy) Using Seismic Scattering and Absorption.” Frontiers in Earth Science 10. https://doi.org/10.3389/feart.2022.852510

work page doi:10.3389/feart.2022.852510 2022
[9]

Earthquake Rupture Velocity and Stress Drop Interaction in the Campi Flegrei V olcanic Caldera

Di Toro, G. 2025. “Earthquake Rupture Velocity and Stress Drop Interaction in the Campi Flegrei V olcanic Caldera.” Communications Earth & Environment. https://doi.org/10.1038/s43247-025-02808-x

work page doi:10.1038/s43247-025-02808-x 2025
[10]

Tracking Transient Changes in the Plumbing System at Campi Flegrei Caldera

Giacomuzzi, G., C. Chiarabba, F. Bianco, P. De Gori, and N. Piana Agostinetti. 2024. “Tracking Transient Changes in the Plumbing System at Campi Flegrei Caldera.” Earth and Planetary Science Letters 637: 118744. https://doi.org/10.1016/j.epsl.2024.118744

work page doi:10.1016/j.epsl.2024.118744 2024
[11]

Burst-Like Swarms in the Campi Flegrei Caldera Accelerating Unrest from 2021 to 2024

Giudicepietro, Flora, Rosario Avino, Eliana Bellucci Sessa, Andrea Bevilacqua, Manuela Bonano, Stefano Caliro, et al. 2025. “Burst-Like Swarms in the Campi Flegrei Caldera Accelerating Unrest from 2021 to 2024.” Nature Communications 16 (1): 1548. https://doi.org/10.1038/s41467-025-56723-y

work page doi:10.1038/s41467-025-56723-y 2025
[12]

The Campi Flegrei Caldera: Historical Revision and New Data on Seismic Crises, Bradyseisms, the Monte Nuovo Eruption and Ensuing Earthquakes (Twelfth Century–1582 AD)

Guidoboni, Emanuela, and Caterina Ciuccarelli. 2011. “The Campi Flegrei Caldera: Historical Revision and New Data on Seismic Crises, Bradyseisms, the Monte Nuovo Eruption and Ensuing Earthquakes (Twelfth Century–1582 AD).” Bulletin of Volcanology 73: 655–677. https://doi.org/10.1007/s00445-010-0430-3

work page doi:10.1007/s00445-010-0430-3 2011
[13]

A Deformation-Driven Earthquake Interaction Model for Seismicity at Campi Flegrei

Hainzl, S., T. Dahm, and A. Tramelli. 2026. “A Deformation-Driven Earthquake Interaction Model for Seismicity at Campi Flegrei.” Communications Earth & Environment 7: 244

work page 2026
[14]

Campi Flegrei: The Relationship between Ground Uplift and Seismic Activity Has Been Quantified

INGV . 2024. “Campi Flegrei: The Relationship between Ground Uplift and Seismic Activity Has Been Quantified.” Press release, December 3, 2024

work page 2024
[15]

Campi Flegrei Monthly and Weekly Surveillance Bulletins, February 2026

INGV-Osservatorio Vesuviano. 2026. “Campi Flegrei Monthly and Weekly Surveillance Bulletins, February 2026”

work page 2026
[16]

Seismological Monitoring of Campi Flegrei Caldera

Martini, M., F. Giudicepietro, L. D’Auria, M. Orazi, G. Borriello, C. Buonocunto, et al. 2008. “Seismological Monitoring of Campi Flegrei Caldera.” Presented at the European Geosciences Union General Assembly 2008, Vienna, Austria. https://www.earth-prints.org/handle/2122/4168

work page 2008
[17]

Forecasting the Preparatory Phase of Induced Earthquakes by Recurrent Neural Network

Picozzi, Matteo, and Antonio Giovanni Iaccarino. 2021. “Forecasting the Preparatory Phase of Induced Earthquakes by Recurrent Neural Network.” Forecasting 3 (1): 17–36. https://doi.org/10.3390/forecast3010002

work page doi:10.3390/forecast3010002 2021
[18]

Coupled Earthquakes and Resonance Processes during the Uplift of Campi Flegrei Caldera

Rapagnani, Giacomo, Simone Cesca, Gilberto Saccorotti, Gesa Petersen, Torsten Dahm, Francesca Bianco, and Francesco Grigoli. 2025. “Coupled Earthquakes and Resonance Processes during the Uplift of Campi Flegrei Caldera.” Communications Earth & Environment 6 (1): 607. https://doi.org/10.1038/s43247-025-02604-7

work page doi:10.1038/s43247-025-02604-7 2025
[19]

A novel single-crystal & single-pass source for polarisation- and colour-entangled photon pairs,

Sabbarese, C., F. Ambrosino, G. Chiodini, et al. 2020. “Continuous Radon Monitoring during Seven Years of V olcanic Unrest at Campi Flegrei Caldera (Italy).” Scientific Reports 10: 9551. https://doi.org/10.1038/s41598- 020-66590-w

work page doi:10.1038/s41598- 2020
[20]

The Phlegrean Fields V olcanological Evolution

Sbrana, A., P. Marianelli, and G. Pasquini. 2021. “The Phlegrean Fields V olcanological Evolution.” Journal of Maps 17 (2): 557–570. https://doi.org/10.1080/17445647.2021.1982033

work page doi:10.1080/17445647.2021.1982033 2021
[21]

Beroza, William L

Tan, Xing, Anna Tramelli, Sergio Gammaldi, Gregory C. Beroza, William L. Ellsworth, and Warner Marzocchi

work page
[22]

A Clearer View of the Current Phase of Unrest at Campi Flegrei Caldera

“A Clearer View of the Current Phase of Unrest at Campi Flegrei Caldera.” Science 390 (6768): 70–75. https://doi.org/10.1126/science.adw9038

work page doi:10.1126/science.adw9038
[23]

b Value Enlightens Different Rheological Behaviour in Campi Flegrei Caldera

Tramelli, A., V . Convertito, and C. Godano. 2024. “b Value Enlightens Different Rheological Behaviour in Campi Flegrei Caldera.” Communications Earth & Environment 5: 275. https://doi.org/10.1038/s43247-024- 01447-y

work page doi:10.1038/s43247-024- 2024
[24]

Noise -Based Seismic Monitoring of the Campi Flegrei Caldera

Zaccarelli, Lucia, and Francesca Bianco. 2017. “Noise -Based Seismic Monitoring of the Campi Flegrei Caldera.” Geophysical Research Letters 44. https://doi.org/10.1002/2016GL072477

work page doi:10.1002/2016gl072477 2017