Abstract
Volcanic eruptions necessitate precise monitoring of magma pressure and inflation for improved forecasting. Understanding deep magma storage is crucial for hazard assessment, yet imaging these systems is challenging due to complex heterogeneities that disrupt standard seismic migration techniques. Here we map the magmatic and hydrothermal system of the La Soufrière volcano in Guadeloupe by analyzing seismic noise data from a sparse geophone array under a matrix formalism. Seismic noise interferometry provides a reflection matrix containing the signature of echoes from deep heterogeneities. Using wave correlations resistant to disorder, matrix imaging successfully unscrambles wave distortions, revealing La Soufrière's internal structure down to 10 km with 100 m resolution. This method surpasses the diffraction limit imposed by geophone array aperture, providing crucial data for modeling and high-resolution monitoring. We see matrix imaging as a revolutionary tool for understanding volcanic systems and enhancing observatories' abilities to monitor dynamics and forecast eruptions.