Entre septiembre de 2024 y septiembre de 2025, la actividad eruptiva del volcán Sabancaya evidenció una transición, desde una fase de disminución progresiva de la actividad hacia un evento explosivo violento que tuvo lugar el 13 de septiembre de 2025. Durante la mayor parte del periodo de análisis (septiembre de 2024 a julio de2025), se observó el descenso sostenido de actividad sísmica, el cese de las emisiones de ceniza (a fines de enero de 2025), deformación mínima del edificio volcánico (0.2 ± 0.4 cm/año), emisiones reducidas de dióxido de azufre (con un promedio de 129 toneladas diarias) y ausencia de anomalías térmicas significativas entre febrero y junio de2025. Sin embargo, el 13 de septiembre de 2025, el volcán registró la explosión más energética desde el inicio de su actual periodo eruptivo en 2016, con una energía sísmica estimada en 380 Megajoules (MJ)y una columna eruptiva que superó los 5000 metros sobre la cima del volcán. El evento generó flujos piroclásticos incandescentes que descendieron por las laderas hasta una distancia de 700 metros, proyectiles balísticos con alcances de hasta 2.3 kilómetros y una dispersión de ceniza superior a 50 kilómetros, en dirección oeste, sur y sureste. La explosión destruyó parcialmente el domo de lava, originando una cavidad de aproximadamente 105 metros de profundidad. Previo a dicho evento, los parámetros de monitoreo permanecieron estables, con excepción de un ligero incremento en la ocurrencia de sismos de tipo tornillo (44 eventos en agosto de 2025, de los 121 registrados en total durante el último año), lo que ha sido asociado a un proceso de presurización interna previo a la explosión. Finalmente, el análisis conjunto de todos los datos permite sugerir tres escenarios sobre la evolución futura del volcán Sabancaya: una actividad explosiva moderada acompañada por la formación de un nuevo domo de lava; una fase de relativa estabilidad caracterizada por actividad fumarólica sostenida; o una intensificación de los episodios explosivos, lo que implicaría un incremento significativo del riesgo eruptivo.

Se describe la actividad eruptiva del volcán Sabancaya desarrollada entre agosto de 2023 y agosto de 2024 a partir del análisis de la variación de datos multiparamétricos, incluyendo la sismicidad, deformación del terreno, emisiones de dióxido de azufre (SO₂), características de las emisiones de cenizas y anomalías térmicas. Durante este período, el volcán presentó una actividad explosiva constante y moderada, caracterizada por la ocurrencia frecuente de explosiones vulcanianas y la formación y evolución de domos de lava. Se registraron 37 explosiones diarias, con predominio de sismos de tipo Largo Período (LP) asociados a la circulación de fluidos magmáticos. En el sector norte del volcán Sabancaya se observó una tasa de inflación vertical de 1.5 ± 0.3 cm/año, inferior a la media observada desde 2013. Las emisiones de ceniza alcanzaron alturas entre 1500 y 3500 m sobre el borde del cráter. El promedio diario del flujo de SO₂ fue de 807 toneladas, con picos que superaron las 4000 toneladas diarias. Se identificaron anomalías térmicas con niveles de Potencia de Radiación Volcánica (VRP) entre 4 MW y 30 MW, y picos de 55 MW. Se identificaron tres fases para el crecimiento de un domo: 1) la primera fase en diciembre 2023 a marzo 2024; 2) una segunda fase de destrucción parcial entre abril y mayo 2024, y 3) un crecimiento significativo en julio 2024. Estos cambios coincidieron con oscilaciones en la deformación del terreno e incremento en la actividad sísmica y emisiones de ceniza. El análisis integrado de los datos indica la presencia de un sistema magmático caracterizado por continuos aportes de magma asociados a ciclos de presurización y despresurización desarrollados al interior del volcán. Según los datos de monitoreo, para un futuro próximo, es probable que la actividad eruptiva continúe, aunque no se puede descartar la posibilidad de un ligero incremento.

Se evalúa la atenuación de las ondas de coda Qc-1 en el Complejo Volcánico Ampato-Sabancaya (CVAS) y volcán Hualca Hualca (VHH), a partir del análisis de 4636 sismos Volcano-Tectónicos (VT) ocurridos durante el periodo 2020-2022, con magnitudes de M1.3 a M5.2 y profundidades menores a 17 km con el objetivo de identificar cambios asociados al proceso eruptivo del volcán Sabancaya. Se aplicó la metodología de Retrodispersión Simple y frecuencias de filtrado de 2, 4, 8 y 16 Hz. Se analizó la atenuación Qc-1 promedio, la variación temporal para la frecuencia de filtrado de 2 Hz, así como la distribución de la atenuación en profundidad mediante la pseudo-tomografía de Qc-1. Los resultados muestran que, la dependencia frecuencial promedio estimada para el CVAS y VHH tiene la forma Qc⁻¹= 0.021±0.016f ⁰.⁷⁴±⁰.¹⁸. sugiriendo que la zona de estudio presenta alta atenuación, comparada con su entorno. También, la variación temporal de la atenuación, analizada para la frecuencia de 2 Hz, evidencia cambios importantes los primeros meses del 2020, entre febrero y junio de 2021, y entre julio y agosto de 2022, asociadas con periodos de mayor y menor actividad eruptiva del volcán Sabancaya. Mientras tanto, el análisis de la pseudo-tomografía de Qc⁻¹ permitió identificar anomalías de alta atenuación asociadas al sistema hidrotermal del Complejo Volcánico Ampato Sabancaya y su interconexión con el posible reservorio magmático ubicado por debajo del volcán Hualca Hualca, entre 8 y 12 km; así como con el alto grado de fracturamiento del sistema de fallas Huambo-Cabanaconde, posiblemente en la etapa de cierre de fracturas.

La aplicación del método geofísico de magnetotelúrica ha permitido caracterizar la estructura interna del estratovolcán Misti. Se ha identificado en superficie el sistema hidrotermal, el cual se extiende lateralmente bajo todo el edificio volcánico, entre 1 km y 2 km sobre el nivel del mar. En la dirección perpendicular al cráter, este sistema es alimentado con fluidos magmáticos (principalmente gases volcánicos), a través de un conducto que se prolonga en profundidad. Asimismo, se identificó un dominio resistivo que abarca casi toda la estructura profunda del volcán, vinculado a rocas intrusivas que conforman el basamento. Otro hallazgo significativo es el cuerpo conductor ubicado ligeramente al este del cráter, a unos 10 km bajo el nivel del mar, interpretado como la cámara magmática de composición andesítica del sistema volcánico del Misti. La escasa actividad sísmica concentrada en la parte superior del edificio volcánico y los valores de resistividad ligeramente altos de la cámara magmática, indican una recarga de magma reducida o limitada. Los resultados obtenidos muestran la primera imagen en profundidad del sistema magmático del Misti y evidencian que se mantiene activo.

The region of Paucarani zone, Tacna (Peru) is formed mainly by mountains, hills, and volcanoes, in addition to wetlands. It is mainly composed of dasitic and andesitic rocks from the Holocene. The regional structural geology is characterized by fault planes in the NW–SE direction, which correspond to the fault system of Apurimac-Caylloma-Maure. This research work presents models of the geoelectrical structure after the analysis, processing, inversion, and interpretation of magnetotelluric data collected in the geothermal area of Paucarani, in order to determine areas of high electrical conductivity that are generally associated with anomalous temperatures, which are of great interest for the possible generation of geothermal energy. The data were acquired by the Instituto Geológico Minero Metalúrgico of Perú in the 2017 campaign covering an area of 120 km2. Forty-three longband MT soundings with frequencies ranging from 10−3 Hz to 103 Hz were processed and distributed in 9 NE–SW profiles. The distribution of soundings was conditioned by the varied topography. Data processing required several steps, such as the determination of apparent impedances, dimensionality analysis using the WALDIM method, strike analysis from geological (regional faults) and geophysical (magnetic prospecting) information. After processing the raw data, the 2D inversion was performed using the WinGlink software. This program builds a geoelectric model from data inversion by iterative processes and smoothness constraints from an initial model that is composed of a mesh with a given resistivity. As a product of all these steps, geoelectrical sections were generated capable of identifying resistivity variations throughout the study area. The results indicate that the most conductive zones, of greatest interest for the study of geothermal sources, are located in the western part of the study area.