A preprint posted to arXiv (arXiv:2604.14196) titled 'Toward buoyancy-driven flow at Campi Flegrei: coupled phase change and asymmetric geometry' uses numerical fluid-dynamics simulations to examine how buoyancy, rather than simple overpressure, can drive magma movement beneath the caldera west of Naples. The authors couple thermal convection with phase transitions — melting of host rock and crystallization of magma — inside a deliberately non-symmetric computational domain that approximates the real caldera's collapsed and faulted structure. As a purely computational study there is no sample size in the experimental sense; the authors instead run suites of simulations varying geometry, viscosity contrasts, and phase-change parameters. They openly note limitations: idealized two-dimensional slices, uncertain rock-property inputs, and the absence of full three-dimensional fault networks or chemical heterogeneity.

This modeling goes well beyond the generic 'magma is moving' statements common in mainstream coverage of Campi Flegrei's recent unrest. Since 2005 the Pozzuoli area has experienced approximately 1.2 m of uplift and repeated earthquake swarms, including more than 600 events in a single month in 2024. Most news reports treat these observations as isolated signals of possible eruption without linking them to specific geophysical drivers. The preprint demonstrates that asymmetric geometry creates preferential flow corridors; buoyancy-driven upwelling is focused along one flank while downwelling of cooler, denser material occurs on the other. Phase changes act as an energy buffer: latent heat absorbed during melting or released during crystallization can sustain slow ground inflation for years to decades without producing the rapid pressure spikes needed for eruption.

Synthesizing this with two peer-reviewed studies sharpens the insight. A 2022 paper in Geophysical Research Letters (Trasatti et al., DOI:10.1029/2022GL098180) inverted InSAR and GPS data to infer a shallow sill-like intrusion at 2–3 km depth that matches the buoyancy-driven emplacement levels explored in the new model. Separately, a 2023 Nature Communications study (Chiodini et al., DOI:10.1038/s41467-023-37877-5) documented increased CO₂ flux consistent with degassing from crystallizing magma — exactly the phase-change process highlighted in the preprint. Together these sources paint a caldera that is episodically intruded by small buoyant batches rather than a single giant chamber approaching critical overpressure.

Historical patterns reinforce the analysis. The 1538 Monte Nuovo eruption followed decades of uplift and seismicity similar to today, yet many earlier unrest episodes (including the 1982–84 crisis with 1.8 m uplift) ended without eruption. The new model offers a mechanistic explanation: asymmetric geometry can divert magma laterally along faults, while phase changes dissipate the thermal energy that would otherwise drive runaway ascent. Mainstream reporting has largely missed this nuance, defaulting to binary 'will it erupt soon?' framing that both alarms and under-informs the roughly three million people living in the greater Naples metropolitan area.

The preprint is not yet peer-reviewed, so its conclusions remain provisional. Nonetheless it supplies a missing bridge between detailed geophysical mechanisms and practical hazard assessment. Future work should embed these coupled equations in data-assimilation frameworks that ingest real-time seismic and deformation streams from Italy's INGV monitoring network. For now the clearest takeaway is that volcanic risk at Campi Flegrei is not uniformly distributed; certain districts sit above modeled upwelling zones while others overlie downwelling limbs. Updated hazard maps that incorporate geometric and phase-change controls could therefore sharpen civil-protection planning far more effectively than generic caldera-wide alerts.