Predicting the ocean state in support of human activities, environmental monitoring, and policymaking across different regions worldwide is fundamental. To properly address physical, dynamical, ice, and biogeochemical processes, numerical strategies must be employed. The authors provide an outlook on the status of operational ocean forecasting systems in eight key regions including the global ocean: the West Pacific and Marginal Seas of South and East Asia, the Indian Seas, the African Seas, the Mediterranean and Black Sea, the North East Atlantic, South and Central America, North America (including the Canadian coastal region, the United States, and Mexico), and the Arctic. The authors initiate their discussion by addressing the specific regional challenges that must be addressed and proceed to discuss the numerical strategy and the available operational systems, ranging from regional to coastal scales. This compendium serves as a foundational reference for understanding the global offering, demonstrating how the diverse physical environment – ranging from waves to ice – and the biogeochemical features besides ocean dynamics can be systematically addressed through regular, coordinated prediction efforts.

In order to understand the contribution of the continental crust on arc magmas, we proceeded with a geochemical study of the Pleistocene eruptive products emitted by Sara Sara volcano in Southern Peru. This quaternary stratovolcano, and the surrounding monogenetic cones, are located at the northern arc termination of the Central Volcanic Zone (CVZ) of the Andes, an arc segment constructed on a thick continental crust. Based on major, trace elements and Sr-Nd-Pb isotopic ratios, we demonstrate that the observed geochemical variations are related to an assimilation-fractional crystallisation (AFC) process dominated by plagioclase and amphibole fractionation coupled with assimilation of the upper continental crust. To decipher the magmatic processes located in the lower continental crust, we extrapolated these geochemical trends towards more basic terms using the compositions of surrounding monogenic volcanoes. We achieve this by using trace element ratios vs. SiO₂diagrams to estimate the hypothetical composition of Sara Sara’s parental magmas before their differentiation by fractional crystallisation. These hypothetical compositions are similar to those of the cones in the AndahuaHuambo and Yura monogenic volcanic fields, representing the most primitive lavas in the Peruvian CVZ. By doing so, we demonstrate that the parental andesitic magmas of Sara Sara volcano exhibit a garnet-rutile signature, characterised by HREE depletion, and elevated Dy/Yb, Sr/Y and Nb/Ta ratios. We discuss the origin of this signature and conclude that it is acquired by partial melting of mafic lithologies leaving residues with garnet and rutile at the deep continental crust. These processes formed enriched basaltic melts and rhyo-dacitic melts which were mixed with mantle-derived magma, thereby generating hybrid andesitic parental magmas. In our scenario, arc andesites are generated in two steps including lower continental crust partial melting and subsequent hybridation with mantle-derived magmas, followed by an assimilation-fractional process in the upper crust. One sentence summary: In the context of the termination of a continental volcanic arc, the high-pressure melting process at the base of the continental crust allows the andesites to acquire a garnet-rutile signature, which is progressively modified by assimilation-crystallisation processes during differentiation at lower pressures.

In this paper, we outline the need for a coordinated international effort toward the building of an open-access Global Ocean Oxygen Database and ATlas (GO₂DAT) complying with the FAIR principles (Findable, Accessible, Interoperable, and Reusable). GO₂DAT will combine data from the coastal and open ocean, as measured by the chemical Winkler titration method or by sensors (e.g., optodes, electrodes) from Eulerian and Lagrangian platforms (e.g., ships, moorings, profiling floats, gliders, ships of opportunities, marine mammals, cabled observatories). GO₂DAT will further adopt a community-agreed, fully documented metadata format and a consistent quality control (QC) procedure and quality flagging (QF) system. GO₂DAT will serve to support the development of advanced data analysis and biogeochemical models for improving our mapping, understanding and forecasting capabilities for ocean O₂ changes and deoxygenation trends. It will offer the opportunity to develop quality-controlled data synthesis products with unprecedented spatial (vertical and horizontal) and temporal (sub-seasonal to multi-decadal) resolution. These products will support model assessment, improvement and evaluation as well as the development of climate and ocean health indicators. They will further support the decision-making processes associated with the emerging blue economy, the conservation of marine resources and their associated ecosystem services and the development of management tools required by a diverse community of users (e.g., environmental agencies, aquaculture, and fishing sectors). A better knowledge base of the spatial and temporal variations of marine O₂ will improve our understanding of the ocean O₂ budget, and allow better quantification of the Earth’s carbon and heat budgets. With the ever-increasing need to protect and sustainably manage ocean services, GO₂DAT will allow scientists to fully harness the increasing volumes of O₂ data already delivered by the expanding global ocean observing system and enable smooth incorporation of much higher quantities of data from autonomous platforms in the open ocean and coastal areas into comprehensive data products in the years to come. This paper aims at engaging the community (e.g., scientists, data managers, policy makers, service users) toward the development of GO₂DAT within the framework of the UN Global Ocean Oxygen Decade (GOOD) program recently endorsed by IOC-UNESCO. A roadmap toward GO₂DAT is proposed highlighting the efforts needed (e.g., in terms of human resources).

The northwestern Peruvian Amazon (NWPA) basin (78.4–75.8° W, 7.9–5.4° S) is an important region for coffee and rice production in Peru. Currently, no prediction models are available for estimating rainfall in advance during the wet season (January–February–March, JFM). Hence, we developed multiple linear regression (MLR) models using predictors derived from sea surface temperature (SST) indices of the Pacific, Atlantic, and Indian Oceans, including central El Niño (C), eastern El Niño (E), tropical South Atlantic (tSATL), tropical North Atlantic (tNATL), extratropical North Atlantic (eNATL), and Indian Ocean basin-wide with E and C removed (IOBW*) indices. Additionally, we utilized large-scale convection indices, namely, the eastern Pacific intertropical convergence zone (ITCZe) and South American Monsoon System (SAMSi) indices, for the 1981–2018 period. Rainfall in the lowland NWPA exhibits a bimodal annual cycle, whereas rainfall in the highland NWPA exhibits a unimodal annual cycle. The MLR model can be used to accurately capture the interannual variability during the wet season in the highland NWPA by utilizing predictors derived from the C and SAMSi indices. In contrast, regarding rainfall in the lowland NWPA, the Pacific SST variability, SAMS and tropical North Atlantic index were relevant. For long lead times, the MLR model provided reliable forecasts of JFM rainfall anomalies in the highlands (R3, approximately 2700 m asl) as these regions are governed by Pacific variability. However, the MLR model exhibited limitations in accurately estimating the wettest JFM season in the highlands due to the absence of a predictor for the amplified effect of the Madden–Julian Oscillation on rainfall.

Water Stable Isotopes (δ18O, δ2H) are valuable tools for tracing sources and interactions in the water cycle, providing important information dedicated to understanding physical mechanisms related to global climate. Despite their significance, the topic of isotopic research in South America has been hindered by limited data. To address this gap, we launched a national-level water stable isotope dataset covering different water sources in Peru (WSI-PeruDB). The dataset contains curated in-house data and incorporates previously published records from various locations collected between 2000 and 2021. The WSIPeruDB dataset is composed of 489 water collection sites and allows a comprehensive use of the dataset by implementing standardized metadata templates containing essential geographical information such as latitude, longitude, and altitude (from sea level to 5000 m a.s.l), and sampling information such as sample type (e.g. groundwater, precipitation, river, spring, and others) and sampling frequency (e.g. biweekly, daily, monthly). The WSIPeruDB dataset is publicly available on Zenodo, facilitating access and use for the scientific community.

We present results of a dedicated global storm-time model for the reconstruction of ionospheric electron density in three-dimensions. Using the storm criterion of |Dst| ≥ 50 nT and Kp ≥ 4, the model is constructed using a combination of radio occultation and ionosonde data during the periods of 2006–2021 and 2000–2020, respectively. From the ionosonde data, only the bottomside electron density profiles up to the maximum height of the F2 layer (hmF2) are considered. In addition to the selection of storm-time data only for the model development, we have investigated the inclusion of time history for the geomagnetic storm indicator Kp at 9 and 33 h in an attempt to take into account the delay of physical processes related to atmospheric gravity waves or traveling ionospheric disturbances and thermospheric composition changes which drive varying ionospheric storm effects during storm conditions. Based on incoherent scatter radar data and in comparison with the IRI 2020 model, the developed storm-time model provides foF2 modelling improvement of above 50% during the storm main phase over Millstone Hill (42.6°N, 71.5°W) and Tromsø (69.6°N, 19.2°E) for the storm periods of 3–6 November 2021 and 23–25 March 2023, respectively. Modelled results for Jicamarca (11.8°S, 77.2°W) show that the storm-time model estimates foF2 by an improvement of over 20% during the main phase of the 07–10 September 2017 storm period. As the ionospheric conditions return to quiet time levels, the IRI 2020 model perform better than the constructed storm -time model.

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.

The Central Volcanic Zone of the Andes (CVZA) has been the focus of volcanological research for decades, becoming a very important site to understand a number of volcanic processes. Despite most of the research in the CVZA being carried out by foreign scientists, the last two decades have seen a significant increase in contributions by regional researchers. This surge has been facilitated by the creation of new volcanic observatories, improvement of the monitoring networks, creation of postgraduate programs where new local volcanologists are trained, creation of specialized research nuclei or groups, and increasing investment in research. This article presents a review of the evolution of the contributions of the regional volcanological community to the knowledge of the CVZA in the last 20 years (2000–2019), both from research and monitoring institutions in Peru, Bolivia, Argentina, and Chile. Based on updates made by the regional groups, a new list of active/potentially active volcanoes of the CVZA is presented, as is a complete database for article published on the CVZA. We find that a significant motivator has been regional volcanic unrest that has triggered new investment. Perú is the country with the highest investment in monitoring and research and is the best instrumented, Argentina is the country with the highest number of local participation in published papers in the domain of volcanology and magmatic systems, and Chilean volcanoes are the focus of the highest number of articles published. The current situation and general projections for the next decade (2020–2030) are also presented for each country, where we believe that the over the next 10 years, will be increased the monitoring and research capabilities, improved the scientific knowledge with more participation of regional institutions, and strengthen the collaboration and integrated work between CVZA countries, especially in border volcanoes.

Un sistema ciclónico, denominado Yaku, se configuró en medio de las condiciones anómalas de El Niño costero 2023 durante marzo, contribuyendo a un conjunto de cambios en el sistema acoplado océano-atmósfera frente a la costa peruana. Datos satelitales de temperatura superficial del mar (TSM), esfuerzo del viento y precipitación (PR) fueron analizados (anomalías) para marzo de 2023. Los resultados revelaron que además de las fuertes anomalías positivas para la TSM y PR, el afloramiento costero se vio fuertemente debilitado, especialmente entre las latitudes de los 3° S a los 12° S, del 7 al 13 de marzo, periodo en el que el ciclón Yaku estuvo más próximo a la costa peruana.

La gran preocupación de muchos peruanos en relación a El Niño en este momento es si el verano será como el de los años 1982-1983 o 1997 1998, es decir, eventos El Niño extraordinarios. El peor escenario concebible sería que ocurra un evento de esa magnitud y que el Perú no esté preparado. Ante la incertidumbre, el ENFEN dio el aviso de la posibilidad de tal evento y se están tomando medidas. De todas formas, es claro que es importante tener información más concreta sobre la probabilidad de estos intensos eventos, lo cual implica entender sus mecanismos y ser capaces de modelarlos correctamente. De gran urgencia el presente año, ¿son nuestra comprensión y los modelos científicos lo suficientemente buenos para la predicción de un El Niño extraordinario?

To better understand the recent internal structure of Misti volcano, we determined a 3D S-wave velocity model applying Ambient Noise Tomography (ANT). We used data from 23 broadband and short-period seismic stations temporarily installed at Misti volcano between March and December 2011. This dataset allowed us to obtain empirical Green's functions by cross-correlating seismic ambient noise signals. Then, we retrieved 104 dispersion curves using the frequency-time analysis (FTAN) and, through a non-linear multiscale inversion, we obtained nine 2-D Rayleigh waves group velocity maps for periods in the range 0.7 s - 2 s. Finally, we carried out the depth inversion through a Bayesian transdimensional inversion to obtain a 3-D S-wave velocity model down to 3 km depth. Our study highlights five relevant seismic velocity anomalies. We observed the presence of three high-velocity zones located in the west-northwest, southwest and southeast parts of the crater, that could be related to intrusive bodies possibly associated with the formation of Misti volcano. We also observed two low-velocity anomalies in the volcano's western and central parts, which coincide with previous studies' findings and are related to fractured and weakened materials associated with the external caldera collapse and recent eruption episodes.

Se presenta los resultados del procesamiento y tratamiento de imágenes satelitales Landsat (TM, ETM+ y OLI) del volcán Ubinas de los periodos eruptivos 2006-2009, 2013-2017 y 2019, para la detección y seguimiento de anomalías térmicas empleando el sistema VOLCANOMS. Se analizaron 53 escenas de imágenes Landsat 8 y 7 del periodo 2006-2019, identificándose 23 anomalías térmicas. El valor más alto de radiancia registrada ocurrió el 11 de mayo del 2007, con 277 W/ m²srµm para la banda SWIR1, 125 W/m²srµm para la banda SWIR2 y 23 W/ m²srµm para la banda TIR. Debido a las diferencias en la frecuencia de toma de datos y número limitado de anomalías identificadas, no se ha establecido una relación clara entre los valores de radiancia térmica con los datos de sismicidad, aunque estas anomalías se han estado observando durante periodos de intensa actividad sísmica de tipo largo periodo, asociadas a paso de fluidos. Los parámetros físicos estimados indican que las anomalías estarían relacionadas a presencia de cuerpos de lava en la superficie del cráter para los periodos eruptivos 2006- 2009 y en algunos registros del periodo 2013-2017, sin embargo, para el proceso eruptivo 2019 las anomalías estarían relacionado a procesos de desgasificación principalmente.

En este estudio se analiza la actividad sísmica detectada en la zona del volcán Coropuna entre mayo de 2018 a abril de 2020. La clasificación sísmica permitió identificar principalmente eventos sismovolcánicos de tipo Volcano-Tectónico (VT), Largo Periodo (LP) y Tremor (TRE), los cuales representan hasta el 93 % del total de sismos, además de un grupo pequeño de sismos asociados a la dinámica del glaciar (DG). Respecto a los sismos de tipo VT, estos ocurrieron en tiempo y zonas distintas situadas entre 18 km (zona 1) y <8 km (zona 2) al noreste y oeste del Coropuna, con magnitudes de entre M2.0 a M3.8 y a menos de 15 km profundidad. De acuerdo con la distribución epicentral y el cálculo de mecanismos focales, la fuente de los sismos se debería a la activación de fallas normales paralelas al sistema de fallas dominantes en los Andes, de dirección NO-SE, NE-NO y O-E. Con relación a los sismos de tipo LP y TRE, estos fueron muy escasos, ya que, en promedio, en el periodo señalado, se registró 1 sismo por día. La localización de estos sismos fue superficial, debajo del volcán Coropuna (zona 2). Los tipos de sismos detectados indican que no hay indicios de ascenso de magma a la superficie; sin embargo, es probable que exista un “cuerpo de magma remanente” asociado a la última erupción del volcán Coropuna (~700 años AP), el cual generaría cambios de temperatura y esfuerzos del medio y el consecuente proceso de fracturamiento de rocas asociado a la ocurrencia de sismos de tipo VT. Asimismo, los cambios de temperatura del medio en contacto con el sistema hidrotermal del Coropuna generarían el movimiento de fluidos (gases magmáticos, vapor de agua, etc.) que ocasionarían sismos de tipo LP y TRE. Finalmente, se considera que el Coropuna presenta un sistema magmático latente, aunque su actividad es calificada como de bajo nivel. No se espera a corto y mediano plazo una reactivación de este volcán.

Se analiza la deformación que ocurre en los alrededores de los volcanes Sara Sara y Cerro Auquihuato, utilizando la técnica de interferometría diferencial de radar de apertura sintética (DInSAR) para el periodo 2014 a 2023. En el caso del volcán Sara Sara, en su flanco oeste y sobre el cono volcánico, se ha observado el levantamiento de la superficie que supera los 8 cm, posiblemente asociado a la acumulación de hielo y nieve. En el flanco noroeste, se han registrado desplazamientos opuestos, con un acercamiento de 15 cm en órbita ascendente y un alejamiento de 30 cm en órbita descendente, probablemente vinculados al desplazamiento lateral del flanco oeste del volcán Sara Sara, como resultado de la dinámica superficial, en especial, el movimiento de rocas originado por el derretimiento de la nieve en las laderas del flanco oeste del volcán. La ausencia de variaciones significativas en la sismicidad, respalda la hipótesis de que las deformaciones registradas no tienen origen en la dinámica interna del volcán. En el volcán Cerro Auquihuato, se han identificado cuatro zonas de deformación. 1) En el punto AUSE (sector sureste) se registró un alejamiento de 17 cm en órbita ascendente y un acercamiento de 7 cm en órbita descendente, indicando un desplazamiento lateral hacia el este, en dirección al río Huaccmi Mayoc. 2) En el punto AUSO (sector noroeste), se observó un alejamiento en ambas órbitas, con un máximo de 17 cm en órbita descendente, sugiriendo un posible hundimiento en esa área. 3) El punto AU MORALES (7 km al sureste del volcán) presenta la mayor deformación, superando los 30 cm, con acercamiento en órbita ascendente y alejamiento en descendente, interpretado como un desplazamiento lateral hacia el oeste, en dirección al río Huaccmi Mayoc. Inspecciones en campo han revelado grietas de 5 cm de ancho en esta zona, evidenciando su inestabilidad. 4) En el punto POMACOCHA (a 18 km al suroeste del volcán), se observa un alejamiento de 6 cm en órbita ascendente y un acercamiento de 17 cm en órbita descendente, indicando un desplazamiento lateral hacia el este, en dirección a la quebrada Alfahuaycco [...].

El presente avance de investigación analiza la distribución de masas de agua durante enero-marzo de 2017 (moderado cálido), 2018 (frío débil), 2019 (cálido débil) y 2020 (neutro) frente a la costa norte del Perú. Para construir un diagrama T-S se utilizan los datos colectados por un flotador (3901231) del programa ARGO que fue lanzado en el año 2016 durante El Niño 2015-2016. Los resultados preliminares muestran una variabilidad en la presencia de las masas agua en la costa norte de Perú.

A fin de evaluar la ocurrencia de peligros naturales en el Valle del Colca (distrito de Chivay), se realiza la revisión bibliográfica disponible y el análisis de la actividad sísmica ocurrida en la zona en los últimos 5 años. Los peligros más recurrentes en el tiempo y que afectan al Valle del Colca y a todas las localidades que se distribuyen en su entorno son los movimientos en masa, erupciones volcánicas y sismos. Al ser los peligros naturales cíclicos, es de esperarse que ellos se repitan en el futuro con la misma o mayor intensidad y con la posibilidad de que los daños y efectos se incrementen de acuerdo al crecimiento y permanencia de las áreas urbanas en zonas reconocidas, históricamente, como de alto riesgo. Ante esta realidad, el riesgo por exposición es alta y más aún si las acciones de Gestión del Riesgo de Desastres no se articulan correctamente entre la población y la propia naturaleza. El Valle del Colca es reconocido como una zona geodinámica joven y por lo tanto, muy activa, puesta en evidencia con la ocurrencia frecuente movimientos en masa, erupciones volcánicas y sismos.

El objetivo del presente avance de investigación es realizar el análisis espectral de ondeletas cruzadas (TOX) entre la serie temporal de precipitación comparada con las series temporales de humedad específica, temperatura superficial del aire e intensidad del viento horizontal antes, durante y después del evento de precipitación intensa ocurrido sobre el observatorio de Huancayo (OHY) el día 6 de octubre de 2019 a las 15:00 hora local. Este evento registró un máximo de precipitación acumulada horaria igual a 6.2 mm h−1. Los resultados del TOX muestran que existe una similitud entre los patrones de periodicidad de todas las series temporales analizadas con un valor mínimo de 0.3 horas (18 minutos) y una periodicidad máxima de 6 horas para la humedad específica y la temperatura del aire y de 8 horas para la intensidad del viento horizontal. En relación a los desfases o retardos temporales, los resultados muestran un desfase mínimo de 180° (9 minutos) y máximo de 90° (1.5 horas) para la humedad específica, mínimo de 180° (9 minutos) y máximo de 230° (3.8 horas) para la temperatura superficial del aire y mínimo de 270° (14 minutos) y máximo de 45° (1 hora) para la intensidad del viento horizontal. Estas diferencias ponen en evidencia el hecho de que los máximos de estas variables ocurren en momentos diferentes dentro de su ciclo diurno y, por lo tanto, tienen un desfase diferente en relación al momento de la ocurrencia de la máxima intensidad de precipitación.

En este estudio se muestran los resultados preliminares del análisis de sensibilidad del sistema acoplado regional CROCO-OASIS-WRF (COW), implementado para el sistema océano–atmósfera del Pacífico sudeste, haciendo uso de la infraestructura computacional del IGP denominada HPC-Linux-Cluster. Para este objetivo, se configuran dos simulaciones acopladas que son casi idénticas, ya que difieren en la versión del modelo atmosférico que están utilizando. Los resultados preliminares muestran que las configuraciones del modelo acoplado simulan con alta correlación el Índice Costero El Niño (ICEN) y, por lo tanto, son útiles para estudiar los procesos de interacción océano-atmósfera en el Pacífico Sudeste.

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.

Dentro de las fases cálida y fría del El Niño-Oscilación del Sur (ENOS), las ondas ecuatoriales juegan un rol importante tanto en su activación como terminación. Basándose en la teoría clásica de ondas, la aproximación lineal de agua somera en un plano ecuatorial beta presenta soluciones de onda que siguen una relación de dispersión, la cual depende de su longitud de onda y frecuencia. El presente trabajo tiene como objetivo identificar las ondas observadas en los datos de anomalía del nivel mar mediante las curvas de dispersión y compararlas con las curvas teóricas. Para esto se calculan los coeficientes de proyección meridional en base a la data observada y se identifican las bandas principales de actividad tanto en la frecuencia como en el número de onda. Los resultados preliminares muestran actividad de ondas de Kelvin concentradas en la banda de 30-100 días y 8880-11100 km (80°-100°), mientras que para las ondas Rossby se encuentran entre 27-37 días y 1110-2220 km (10°-20°).