ALBERT

Description: The geochemical monitoring of volcanic activity today relies largely on remote sensing, but the combination of this approach together with soil gas monitoring, using the appropriate parameters, is still not widely used. The main purpose of this study was to correlate data from crater gas emissions with flank emissions of soil gases at Mt. Etna volcano from June 2006 to December 2020. Crater SO2 fluxes were measured from fixed stations around the volcano using the DOAS technique and applying a modeled clear-sky spectrum. The SO2/HCl ratio in the crater plume was measured with the OP-FTIR technique from a transportable instrument, using the sun as an IR source. Soil CO2 efflux coupled with the 220Rn/222Rn activity ratio in soil gases (named SGDI) were measured at a fixed monitoring site on the east flank of Etna. All signals acquired were subject both to spectral analysis and to filtering of the periodic signals discovered. All filtered signals revealed changes that were nicely correlated both with other geophysical signals and with volcanic eruptions during the study period. Time lags between parameters were explained in terms of different modes of magma migration and storage inside the volcano before eruptions. A comprehensive dynamic degassing model is presented that allows for a better understanding of magma dynamics in an open-conduit volcano.

Description: Published

Description: 1122

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: In recent decades, the Campi Flegrei caldera (Italy) showed unrest characterized by in- creases in seismicity, ground uplift, and hydrothermal activity. Currently, the seismic and hydrother- mal phenomena are mostly concentrated in the Solfatara–Pisciarelli area, which presents a wide fumarolic field and mud emissions. The main fumarole in Pisciarelli is associated with a boiling mud pool. Recently, episodes of a sudden increase in hydrothermal activity and expansion of mud and gas emissions occurred in this area. During these episodes, which occurred in December 2018 and September 2020, Short Duration Events (SDEs), related to the intensity of mud pool boiling, were recorded in the fumarolic seismic tremor. We applied a Self-Organizing Map (SOM) neural network to recognize the occurrence of SDEs in the fumarolic tremor of Campi Flegrei, which provides important information on the state of activity of the hydrothermal system and about the possible phreatic activity. Our method, based on an ad hoc feature extraction procedure, effectively clustered the seismic signals containing SDEs and separated them from those representing the normal fumarolic tremor. This result is useful for improving the monitoring of the Solfatara–Pisciarelli hydrothermal area which is a high-risk zone in Campi Flegrei.

Description: Published

Description: 5505

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: This paper addresses the classification of images depicting the eruptive activity of Mount Etna, captured by a network of ground-based thermal cameras. The proposed approach utilizes Convolutional Neural Networks (CNNs), focusing on pretrained models. Eight popular pretrained neural networks underwent systematic evaluation, revealing their effectiveness in addressing the classification problem. The experimental results demonstrated that, following a retraining phase with a limited dataset, specific networks such as VGG-16 and AlexNet, achieved an impressive total accuracy of approximately 90%. Notably, VGG-16 and AlexNet emerged as practical choices, exhibiting individual class accuracies exceeding 90%. The case study emphasized the pivotal role of transfer learning, as attempts to solve the classification problem without pretrained networks resulted in unsatisfactory outcomes.

Description: Supported by Italian Research Center on High Performance Computing Big Data and Quantum Computing (ICSC), project funded by European Union—NextGenerationEU—and National Recovery and Resilience Plan (NRRP)—Mission 4 Component 2 within the activities of Spoke 3 (Astrophysics and Cosmos Observations). Sonia Calvari also acknowledges the financial support of the Project FIRST ForecastIng eRuptive activity at Stromboli volcano (Delibera n. 144/2020; Scientific Responsibility: S.C.) Vulcani 2019.

Description: Published

Description: 124-137

Description: OSV3: Sviluppo di nuovi sistemi osservazionali e di analisi ad alta sensibilità

Description: JCR Journal

Description: On 21 May 2023, a hidden eruption occurred at the Southeast Crater (SEC) of Etna (Italy); indeed, bad weather prevented its direct and remote observation. Tephra fell toward the southwest, and two lava flows propagated along the SEC’s southern and eastern flanks. The monitoring system of the Istituto Nazionale di Geofisica e Vulcanologia testified to its occurrence. We analyzed the seismic and infrasound signals to constrain the temporal evolution of the fountain, which lasted about 5 h. We finally reached Etna’s summit two weeks later and found an unexpected pyroclastic density current (PDC) deposit covering the southern lava flow at its middle portion. We performed unoccupied aerial system and field surveys to reconstruct in 3D the SEC, lava flows, and PDC deposits and to collect some samples. The data allowed for detailed mapping, quantification, and characterization of the products. The resulting lava flows and PDC deposit volumes were (1.54 ± 0.47) × 106 m3 and (1.30 ± 0.26) × 105 m3, respectively. We also analyzed ground-radar and satellite data to evaluate that the plume height ranges between 10 and 15 km. This work is a comprehensive analysis of the fieldwork, UAS, volcanic tremor, infrasound, radar, and satellite data. Our results increase awareness of the volcanic activity and potential dangers for visitors to Etna’s summit area.

Description: Published

Description: 1555

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: The Sciara del Fuoco (SdF) collapse scar at Stromboli is an active volcanic area affected by rapid morphological changes due to explosive/effusive eruptions and mass-wasting processes. The aim of this paper is to demonstrate the importance of an integrated analysis of multi-temporal remote sensing (photogrammetry, COSMO-SkyMed Synthetic Aperture Radar amplitude image) and marine geophysical data (multibeam and side scan sonar data) to characterize the main morphological, textural, and volumetric changes that occurred along the SdF slope in the 2020–2021 period. The analysis showed the marked erosive potential of the 19 May 2021 pyroclastic density current generated by a crater rim collapse, which mobilized a minimum volume of 44,000 m^3 in the upper Sciara del Fuoco slope and eroded 350,000–400,000 m^3 of material just considering the shallow-water setting. The analysis allowed us also to constrain the main factors controlling the emplacement of different lava flows and overflows during the monitored period. Despite the morphological continuity between the subaerial and submarine slope, textural variations in the SdF primarily depend on different processes and characteristics of the subaerial slope, the coastal area, the nearshore, and “deeper” marine areas.

Description: Published

Description: 4605

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(12),(2022): 3199-3219, https://doi.org/10.1175/jpo-d-22-0009.1.

Description: The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed water-mass transformations are dominated by rough topography “hotspots,” where the bottom enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger upwelling in a thin bottom boundary layer (BBL). These water-mass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downward diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model that approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the nonlocal coupling to the basin-scale circulation.

Description: We acknowledge funding support from National Science Foundation Awards 1536515, 1736109, and 2149080. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant 174530.

Description: 2023-05-18

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(6), (2022): 1091–1110, https://doi.org/10.1175/JPO-D-21-0068.1.

Description: Hundreds of full-depth temperature and salinity profiles collected by Deepglider autonomous underwater vehicles (AUVs) in the North Atlantic reveal robust signals in eddy isopycnal vertical displacement and horizontal current throughout the entire water column. In separate glider missions southeast of Bermuda, subsurface-intensified cold, fresh coherent vortices were observed with velocities exceeding 20 cm s−1 at depths greater than 1000 m. With vertical resolution on the order of 20 m or less, these full-depth glider slant profiles newly permit estimation of scaled vertical wavenumber spectra from the barotropic through the 40th baroclinic mode. Geostrophic turbulence theory predictions of spectral slopes associated with the forward enstrophy cascade and proportional to inverse wavenumber cubed generally agree with glider-derived quasi-universal spectra of potential and kinetic energy found at a variety of locations distinguished by a wide range of mean surface eddy kinetic energy. Water-column average spectral estimates merge at high vertical mode number to established descriptions of internal wave spectra. Among glider mission sites, geographic and seasonal variability implicate bottom drag as a mechanism for dissipation, but also the need for more persistent sampling of the deep ocean.

Description: This work was funded by NSF Grant 1736217 and would not have been possible without the help of Kirk O’Donnell, James Bennett, Noel Pelland, and all contributors to Deepglider development. We additionally thank the captain crew of the R/V Atlantic Explorer and the BATS team at the Bermuda Institute of Ocean Sciences, particularly Rod Johnson, as well as Seakeepers International for their professionalism, capability, and generous assistance in deploying and recovering gliders.

Description: During the last two decades, the Etna volcano has undergone several sequences of lava fountaining (LF) events that have had a major impact on road conditions, infrastructure and the local population. In this paper, we consider the LF episodes occurring between 2011 and 2022, calculating their erupted volumes using the images recorded by the monitoring thermal cameras and applying a manual procedure and a dedicated software to determine the lava fountain height over time, which is necessary to obtain the erupted volume. The comparison between the results indicates the two procedures match quite well, the main differences occurring when the visibility is poor and data are interpolated. With the aim of providing insights for hazard assessment, we have fitted some probabilistic models of both the LF inter-event times and the erupted volumes of pyroclastic material. In more detail, we have tested power-law distributions against log-normal, Weibull, generalised Pareto and log-logistic. Results show that the power-law distribution is the most likely among the alternatives. This implies the lack of characteristic scales for both the inter-event time and the pyroclastic volume, which means that we have no indication as to when a new episode of LF will occur and/or how much material will be erupted. What we can reasonably say is only that short inter-event times are more frequent than long inter-event times, and that LF characterised by small volumes are more frequent than LF with high volumes. However, if the hypothesis that magma accumulates on Etna at a rate of about 0.8 m3 s −1 holds, the material accumulated in the source region from the beginning of the observation period (2011) to the present (2022) has already been ejected. In simple terms, there is no accumulated magma in the shallow storage that is prone to be erupted in the near future.

Description: This research was funded by Project FIRST—ForecastIng eRuptive activity at Stromboli volcano: timing, eruptive style, size, intensity, and duration, INGV-Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020; Scientific Responsibility: S.C.).

Description: Published

Description: 6183

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Offshore data in the western Ionian Sea indicate that the NW–SE-trending dextral shear zone of the Alfeo–Etna Fault System turns to the N–S direction near the Ionian coastline, where the extensional Timpe Fault System is located. Morpho-structural data show that NW–SE-trending right-lateral strikeslip faults connect the Timpe Fault System with the upper slope of the volcano, where the eruptive activity mainly occurs along the N–S to NE–SW-trending fissures. Fault systems are related to the ~E–Wtrending extension and they are seismically active having given rise to shallow and low-moderate magnitude earthquakes in the last 150 years. As a whole, morpho-structural, geodetic and seismological data, seismic profiles and bathymetric maps suggest that similar geometric and kinematic features characterize the shear zone both on the eastern flank of the volcano and in the Ionian offshore. The Alfeo– Etna Fault System probably represents a major kinematic boundary in the western Ionian Sea associated with the Africa–Europe plate convergence since it accommodates, by right-lateral kinematics, the differential motion of adjacent western Ionian compartments. Along this major tectonic alignment, crustal structures such as releasing bends, pull-apart basins and extensional horsetails occur both offshore and on-land, where they probably represent the pathway for magma uprising from depth

Description: This research was funded by the Catania University PIA.CE.RI. Project (linea 2) “Interaction between volcanic activity and active tectonic processes in the Mt. Etna area (InvultEtna). The research has moreover benefited from funding provided by the agreement between Istituto Nazionale di Geofisica e Vulcanologia (INGV) and the Italian Presidenza del Consiglio dei Ministri, Dipartimento della Protezione Civile

Description: Published

Description: 128

Description: 2T. Deformazione crostale attiva

Description: JCR Journal

Description: The re-mobilization of volcaniclastic material poses a hazard factor which, although it decreases with time since the last eruption, remains present in the hydrographic basins of volcanic areas. Herein, we present the results of the numerical modelling of erosive phenomena of volcanic deposits, as well as of flooding in the volcanic area. The proposed approach includes runoff estimation, land use analysis, and the application of hydraulic and erosion modelling. It exploits the Iber software, a widely used and validated model for rainfall-runoff, river flooding, and erosion and sediment transport modelling. The methodology was applied to the Island of Vulcano (Italy), known for the erosion phenomena that affect the slopes of one of its volcanic cones (La Fossa cone). The rainfall excess was calculated using a 19-year dataset of hourly precipitations, and the curve number expressed by the information on soil cover in the area, derived from the land cover and land use analysis. The erosion and flow models were performed considering different rainfall scenarios. Results show a particularly strong erosion, with thicknesses greater than 0.4 m. This is consistent with field observations, in particular with some detailed data collected both after intense events and by long-term observation. Results of the hydraulic simulations show that moderate and torrential rainfall scenarios can lead to flood levels between 0.2 and 0.6 m, which mostly affect the harbours located in the island’s inhabited area.

Description: This project was partially funded by the “Fondi di Ateneo 2022 (ex 60%)” by the Università degli Studi di Firenze (project “VOLFLANK—Use of remote sensing data for the stability analysis of active volcanoes”; P.I.: F.D.T.). A.F. and M.F. carried out this work in the frame of INGV Progetti Ricerca Libera 2022 (project “VOLF—VOlcaniclastic debris flows at La Fossa cone (Volcano Island): evolution and hazard implication”).

Description: Published

Description: 16549

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(12), (2022): 3221–3240, https://doi.org/10.1175/jpo-d-22-0010.1.

Description: Small-scale mixing drives the diabatic upwelling that closes the abyssal ocean overturning circulation. Indirect microstructure measurements of in situ turbulence suggest that mixing is bottom enhanced over rough topography, implying downwelling in the interior and stronger upwelling in a sloping bottom boundary layer. Tracer release experiments (TREs), in which inert tracers are purposefully released and their dispersion is surveyed over time, have been used to independently infer turbulent diffusivities—but typically provide estimates in excess of microstructure ones. In an attempt to reconcile these differences, Ruan and Ferrari derived exact tracer-weighted buoyancy moment diagnostics, which we here apply to quasi-realistic simulations. A tracer’s diapycnal displacement rate is exactly twice the tracer-averaged buoyancy velocity, itself a convolution of an asymmetric upwelling/downwelling dipole. The tracer’s diapycnal spreading rate, however, involves both the expected positive contribution from the tracer-averaged in situ diffusion as well as an additional nonlinear diapycnal distortion term, which is caused by correlations between buoyancy and the buoyancy velocity, and can be of either sign. Distortion is generally positive (stretching) due to bottom-enhanced mixing in the stratified interior but negative (contraction) near the bottom. Our simulations suggest that these two effects coincidentally cancel for the Brazil Basin Tracer Release Experiment, resulting in negligible net distortion. By contrast, near-bottom tracers experience leading-order distortion that varies in time. Errors in tracer moments due to realistically sparse sampling are generally small (〈20%), especially compared to the O(1) structural errors due to the omission of distortion effects in inverse models. These results suggest that TREs, although indispensable, should not be treated as “unambiguous” constraints on diapycnal mixing.

Description: We acknowledge funding support from National Science Foundation Awards 1536515 and 1736109. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant 174530. This research is also supported by the NOAA Climate and Global Change Postdoctoral Fellowship Program, administered by UCAR’s Cooperative Programs for the Advancement of Earth System Science (CPAESS) under Award NA18NWS4620043B.

Description: 2023-05-18

Description: Numerical modelling of tephra fallout is a fast-developing research area in volcanology. Several models are currently available both to forecast the dispersion of volcanic particles in the atmosphere and to calculate the particles deposited at different locations on the ground. Data from these simulations can then be used both to manage volcanic crises (e.g., protect air traffic) or perform long-term hazard assessment studies (e.g., through hazard maps). Given the importance of these tasks, it is important that each model is thoroughly tested in order to assess advantages and limitations, and to provide useful information for quantifying the model uncertainty. In this study we tested the coupled PLUME-MoM/HYSPLIT models by applying them to the Puyehue–Cordon Caulle 2011 sub-Plinian eruption. More specifically, we tested new features recently introduced in these well-established models (ash aggregation, external water addition, and settling velocity models), we implemented a new inversion procedure, and we performed a parametric analysis. Our main results reaffirm the pivotal role played by mass eruption rate on the final deposit and show that some choices for the input parameters of the model can lead to the large overestimation in total deposited mass (which can be reduced with our inversion procedure). The parametric analysis suggests a most likely value of the mass eruption rate in the range 2.0–6.3 × 106 kg/s. More studies with a similar approach would be advisable in order to provide final users with useful indications about the parameters that should be carefully evaluated before being used as input for this kind of model.

Description: Published

Description: 784

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Between December 2020 and February 2022, the South East Crater of Etna has been the source of numerous eruptions, mostly characterized by the emission of lava fountains, pyroclastic material and short-lasting lava flows. Here we estimate the volume and distribution of the lava deposits by elaborating multi-source satellite imagery. SEVIRI data have been elaborated using CL-HOTSAT to estimate the lava volume emitted during each event and calculate the cumulative volume; Pléiades and WorldView-1 data have been used to derive Digital Surface Models, whose differences provide thickness distributions and hence volumes of the volcanic deposits. We find a good agreement, with the total average lava volume obtained by SEVIRI reaching 73.2 × 106 m3 and the one from optical data amounting to 67.7 × 106 m3. This proves the robustness of both techniques and the accuracy of the volume estimates, which provide important information on the lava flooding history and evolution of the volcano.

Description: This work was supported by the INGV project Pianeta Dinamico (CUP D53J19000170001) funded by MIUR (“Fondo finalizzato al rilancio degli investimenti delle amministrazioni centrali dello Stato e allo sviluppo del Paese,” legge 145/2018), Tema 8—PANACEA, Scientific Responsibility: A.C.). The research was also funded by “TUNE—Effusion rate estimates at Etna and Stromboli: constraints imposed by a variety of satellite remote sensing data” (Bando di Ricerca Libera 2019 of INGV; Scientific Responsibility: G.G.). This research was also supported by the Project FIRST—ForecastIng eRuptive activity at Stromboli volcano: timing, eruptive style, size, intensity, and duration, INGV-Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020; Scientific Responsibility: S.C.).

Description: Published

Description: 916

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: The volatiles released by the volcanic structures of the world contribute to natural environmental pollution both during the passive and active degassing stages. The Island of Vulcano is characterized by solfataric degassing mainly localized in the summit part (Fossa crater) and in the peripheral part in the Levante Bay. The normal solfataric degassing (high-temperature fumarolic area of the summit and boiling fluids emitted in the Levante Bay area), established after the last explosive eruption of 1888–90, is periodically interrupted by geochemical crises characterized by anomalous degassing that are attributable to increased volcanic inputs, which determine a sharp increase in the degassing rate. In this work, we have used the data acquired from the INGV (Istituto Nazionale di Geofisica e Vulcanologia) geochemical monitoring networks to identify, evaluate, and monitor the geochemical variations of the extensive parameters, such as the SO2 flux from the volcanic plume (solfataric cloud) and the CO2 flux from the soil in the summit area outside the fumaroles areas. The increase in the flux of volatiles started in June–July 2021 and reached its maximum in November of the same year. In particular, the mean monthly flux of SO2 plume of 22 tons day−1 (t d−1) and of CO2 from the soil of 1570 grams per square meter per day (g m2 d −1) increased during this event up to 89 t d−1 and 11,596 g m2 d −1, respectively, in November 2021. The average annual baseline value of SO2 output was estimated at 7700 t d−1 during normal solfataric activity. Instead, this outgassing increased to 18,000 and 24,000 t d−1 in 2021 and 2022, respectively, indicating that the system is still in an anomalous phase of outgassing and shows no signs of returning to the pre-crisis baseline values. In fact, in the first quarter of 2023, the SO2 output shows average values comparable to those emitted in 2022. Finally, the dispersion maps of SO2 on the island of Vulcano have been produced and have indicated that the areas close to the fumarolic source are characterized by concentrations of SO2 in the atmosphere higher than those permitted by European legislation (40 µg m−3 for 24 h of exposition) on human health.

Description: This research was funded by the INGV-DPCN (Italian National Institute of Geophysics and Volcanology-Italian National Department for Civil Protection) volcanic surveillance program of Vulcano island, ObFu 0304.010. Moreover, this investigation was partially funded by TORS project in the framework of institutional INGV projects (“Ricerca Libera”, ObFu 9999.549 and Pianeta Dinamico Task V2, ObFu 1020.010).

Description: Published

Description: 3086

Description: 4V. Processi pre-eruttivi

Description: JCR Journal

Description: The 2021 eruption at Tajogaite (Cumbre Vieja) volcano (La Palma, Spain) was characterized by Strombolian eruptions, Hawaiian fountaining, white gasdominated and grey ash-rich plumes, and lava effusion from multiple vents. The variety of eruptive styles displayed simultaneously and throughout the eruption presents an opportunity to explore controls on explosivity and the relationship between explosive and effusive activity. Explosive eruption dynamics were recorded using ground-based thermal photography and videography. We show results from the analysis of short (〈5 min) near-daily thermal videos taken throughout the eruption from multiple ground-based locations and continuous time-lapse thermal photos over the period November 16 to November 26. We measure the apparent radius, velocity, and volume flux of the high-temperature gas-and-ash jet and lava fountaining behaviors to investigate the evolution of the explosive activity over multiple time scales (seconds-minutes, hours, and daysweeks). We find fluctuations in volume flux of explosive material that correlate with changes in volcanic tremor and hours-long increases in explosive flux that are immediately preceded by increases in lava effusion rate. Correlated behavior at multiple vents suggests dynamic magma ascent pathways connected in the shallow (tens to hundreds of meters) sub-surface. We interpret the changes in explosivity and the relative amounts of effusive and explosivity to be the result of changes in gas flux and the degree of gas coupling.

Description: Authors from INVOLCAN and ITER acknowledge support under projects VOLRISKMAC (MAC/3.5b/124) and VOLRISKMAC II (MAC2/3.5b/328), financed by the Program INTERREG V-A Spain-Portugal MAC 2014–2020 of the European Commission; Cumbre Vieja Emergencia, financed by the Science and Innovation Ministry, Spanish Government; TFassistance, financed by the Cabildo Insular de Tenerife; and LPvolcano, financed by the Cabildo Insular de La Palma. We are grateful to Antoni Álvarez for his logistical support. Authors from LDEO acknowledge support from the Gordon and Betty Moore Foundation under grant GBMF8995 for the AVERT project, and from NSF under EAR-1654588.

Description: Published

Description: 1193436

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Stromboli is an open-conduit active volcano located in the southern Tyrrhenian Sea and is the easternmost island of the Aeolian Archipelago. It is known as “the lighthouse of the Mediterranean” for its continuous and mild Strombolian-type explosive activity, occurring at the summit craters. Sometimes the volcano undergoes more intense explosions, called “major explosions” if they affect just the summit above 500 m a.s.l. or “paroxysms” if the whole island is threatened. Effusive eruptions are less frequent, normally occurring every 3–5 years, and may be accompanied or preceded by landslides, crater collapses and tsunamis. Given the small size of the island (maximum diameter of 5 km, NE–SW) and the consequent proximity of the inhabited areas to the active craters (maximum distance 2.5 km), it is of paramount importance to use all available information to forecast the volcano’s eruptive activity. The availability of a detailed record of the volcano’s eruptive activity spanning some centuries has prompted evaluations on its possible short-term evolution. The aim of this paper is to present some statistical insights on the eruptive activity at Stromboli using a catalogue dating back to 1879 and reviewed for the events during the last two decades. Our results confirm the recent trend of a significant increase in major explosions, small lava flows and summit crater collapses at the volcano, and might help monitoring research institutions and stakeholders to evaluate volcanic hazards from eruptive activity at this and possibly other open-vent active basaltic volcanoes.

Description: This research was funded by Project FIRST—ForecastIng eRuptive activity at Stromboli volcano: timing, eruptive style, size, intensity, and duration, INGV–Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020; Scientific Responsibility: S.C.).

Description: Published

Description: 4822

Description: 1V. Storia eruttiva

Description: JCR Journal

Description: Constraining the magmatic 3He/4He signature of fluids degassed from a magmatic system is crucial for making inferences on its mantle source. This is especially important in arc volcanism, where variations in the composition of the wedge potentially induced by slab sediment fluids must be distinguished from the effects of magma differentiation, degassing, and crustal contamination. The study of fluid inclusions (FIs) trapped in minerals of volcanic rocks is becoming an increasingly used methodology in geochemical studies that integrates the classical study of volcanic and geothermal fluids. Here, we report on the first noble gas (He, Ne, Ar) concentrations and isotopic ratios of FI in olivine (Ol) and pyroxene (Px) crystals separated from eruptive products of the Telica and Baru volcanoes, belonging to the Nicaraguan and Panamanian arc-segments of Central America Volcanic arc (CAVA). FIs from Telica yield air corrected 3He/4He (Rc/Ra) of 7.2–7.4 Ra in Ol and 6.1–7.3 in Px, while those from Baru give 7.1–8.0 Ra in Ol and 4.2–5.8 Ra in Px. After a data quality check and a comparison with previous 3He/4He measurements carried out on the same volcanoes and along CAVA, we constrained a magmatic Rc/Ra signature of 7.5 Ra for Telica and of 8.0 Ra for Baru, both within the MORB range (8 1 Ra). These 3He/4He differences also reflect variations in the respective arc-segments, which cannot be explained by radiogenic 4He addition due to variable crust thickness, as the mantle beneath Nicaragua and Panama is at about 35 and 30 km, respectively. We instead highlight that the lowest 3He/4He signature observed in the Nicaraguan arc segment reflects a contamination of the underlying wedge by slab sediment fluids. Rc/Ra values up to 9.0 Ra are found at Pacaya volcano in Guatemala, where the crust is 45 km thick, while a 3He/4He signature of about 8.0 Ra was measured at Turrialba volcano in Costa Rica, which is similar to that of Baru, and reflects possible influence of slab melting, triggered by a change in subduction conditions and the contemporary subduction of the Galapagos hot-spot track below southern Costa Rica and western Panama.

Description: Published

Description: 2076-3417

Description: 1V. Storia eruttiva

Description: 2V. Struttura e sistema di alimentazione dei vulcani

Description: 3V. Proprietà chimico-fisiche dei magmi e dei prodotti vulcanici

Description: JCR Journal

Description: The identification of the mechanisms responsible for the deformation of calderas is of primary importance for our understanding of the dynamics of magmatic systems and the evaluation of volcanic hazards. We analyze twenty years (1997–2018) of geodetic measurements on Ischia Island (Italy), which include the Mt. Epomeo resurgent block, and is affected by hydrothermal manifestations and shallow seismicity. The data from the GPS Network and the leveling route show a constant subsidence with values up to 􀀀15 2.0 mm/yr and a centripetal displacement rate with the largest deformations on the southern flank of Mt. Epomeo. The joint inversion of GPS and levelling data is consistent with a 4 km deep source deflating by degassing and magma cooling below the southern flank of Mt. Epomeo. The depth of the source is supported by independent geophysical data. The Ischia deformation field is not related to the instability of the resurgent block or extensive gravity or tectonic processes. The seismicity reflects the dynamics of the shallow hydrothermal system being neither temporally nor spatially related to the deflation.

Description: Published

Description: 4648

Description: 1V. Storia eruttiva

Description: 2V. Struttura e sistema di alimentazione dei vulcani

Description: 4V. Processi pre-eruttivi

Description: 5V. Processi eruttivi e post-eruttivi

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: We propose a multi-temporal-scale analysis of ground deformation data using both high-rate tilt and GNSS measurements and the DInSAR and daily GNSS solutions in order to investigate a sequence of four paroxysmal episodes of the Voragine crater occurring in December 2015 at Mt. Etna (Italy). The analysis aimed at inferring the magma sources feeding a sequence of very violent eruptions, in order to understand the dynamics and to image the shallow feeding system of the volcano that enabled such a rapid magma accumulation and discharge. The high-rate data allowed us to constrain the sources responsible for the fast and violent dynamics of each paroxysm, while the cumulated deformation measured by DInSAR and daily GNSS solutions, over a period of 12 days encompassing the entire eruptive sequence, also showed the deeper part of the source involved in the considered period, where magma was stored. We defined the dynamics and rates of the magma transfer, with a middle-depth storage of gas-rich magma that charges, more or less continuously, a shallower level where magma stops temporarily, accumulating pressure due to the gas exsolution. This machine-gun-like mechanism could represent a general conceptual model for similar events at Etna and at all volcanoes.

Description: Published

Description: 4630

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Immersive virtual reality can potentially open up interesting geological sites to students, academics and others who may not have had the opportunity to visit such sites previously. We study how users perceive the usefulness of an immersive virtual reality approach applied to Earth Sciences teaching and communication. During nine immersive virtual reality-based events held in 2018 and 2019 in various locations (Vienna in Austria, Milan and Catania in Italy, Santorini in Greece), a large number of visitors had the opportunity to navigate, in immersive mode, across geological landscapes reconstructed by cutting-edge, unmanned aerial system-based photogrammetry techniques. The reconstructed virtual geological environments are specifically chosen virtual geosites, from Santorini (Greece), the North Volcanic Zone (Iceland), and Mt. Etna (Italy). Following the user experiences, we collected 459 questionnaires, with a large spread in participant age and cultural background. We find that the majority of respondents would be willing to repeat the immersive virtual reality experience, and importantly, most of the students and Earth Science academics who took part in the navigation confirmed the usefulness of this approach for geo-education purposes.

Description: This research has been provided in the framework of the following projects: (i) the MIUR project ACPR15T4_00098–Argo3D (http://argo3d.unimib.it/ (accessed on 26 November 2021)); (ii) 3DTeLC Erasmus + Project 2017-1-UK01-KA203-036719 (http://www.3dtelc.com (accessed on 26 November 2021)); (iii) EGU 2018 Public Engagement Grant (https://www.egu.eu/outreach/peg/ (accessed on 26 November 2021)). Agisoft Metashape is acknowledged for photogrammetric data processing. This article is also an outcome of Project MIUR–Dipartimenti di Eccellenza 2018–2022. Finally, this paper is an outcome of the Virtual Reality lab for Earth Sciences—GeoVires lab (https://geovires.unimib.it/ (accessed on 26 November 2021)). The work supports UNESCO IGCP 692 ‘Geoheritage for Resilience’.

Description: Published

Description: 9

Description: 1TM. Formazione

Description: JCR Journal

Description: Volcanology, seismology and Earth Sciences in general, like all quantitative sciences, are increasingly dependent on the quantity and quality of data acquired. In recent dec-ades, a marked evolution has characterized Earth sciences towards a greater use of ana-lytical and numerical approaches, shifting these fields from the natural to the physical sciences. Understanding the physical behavior of active volcanoes and faults is critical to as-sess the hazards affecting the population living close to active volcano and seismic areas, and thus to mitigate the risks posed by those threats [1,2]. The knowledge of a physical process requires the acquisition of a huge amount of information (data) on that particular phenomenon. Today, different kinds of data record the processes that operate in volcanic and tec-tonic systems and provide insights that can lead to improved predictions of potential hazards, both immediate and long term. The geoscience community has collected an enormous wealth of data that require further analysis. The diversity and quantity of these geoscience data and collections continue to expand [3]. The increasing amount of data and the availability of new technologies and instru-mentation at an ever-greater rate open new frontiers and challenges for acquiring, trans-mitting, archiving, processing and analyzing the newly available datasets. Guo [4] pre-dicted growth for the general digital universe size of factor 10 from 2016 to 2025. Among all digital data, scientific data are those relevant to the observation of natural phenomena and characterized by non-repeatability, high uncertainty, high dimensionality and a high degree of computational complexity [4]. This means that scientific data need to be well preserved, due to the non-repeatability, and implies a parallel growth of processing capa-bilities to be well exploited. Cheng et al. [5] highlighted the striking growth of Earth Sci-ence data from molecular to astronomical scales and the increasing use of supercompu-ting tools for supporting geoscience research. The authors evidence how, with the contin-uously increasing availability of digital data, Earth Sciences are also turning from the tra-ditional question-driven or problem-driven approach, where scientists seek to find an-swers to known questions, to the new data-driven one where scientists apply a data dis-covery process that might find answers to still unknown questions. In agreement with Cheng et al. [5], we believe that new integrated multi-disciplinary knowledge systems and new data discovery techniques for handling and mining big data for knowledge discovery would spur the integration of transdisciplinary and mul-ti-dimensional Earth science data. Furthermore, this will help the transition from a nar-row focus on separate disciplines to a holistic, comprehensive and integrative focus of the different disciplines linked to the Earth Sciences. With this aim, for this special issue titled “Data Processing and Modeling on Volcan-ic and Seismic Areas”, we invited articles on all aspects of solid Earth Science that made use of data to analyze and model processes related to volcanoes or earthquakes. Manuscripts with various types of analyses, including volcanic ground deformation modeling, seismic swarm characterization and volcanic gas measurement, have been proposed and published. The collection provides an insight into the enormous need for increasingly complex data analysis and modeling techniques to try to describe the natural phenomena here considered. This special issue was introduced to collect the latest research on the processing and modeling of Earth Sciences data, and to address challenging problems with all topics re-lated to volcanoes and seismic areas. Various subjects have been addressed in this collec-tion, mainly on data processing for volcanic studies (three papers), tectonics (two papers) and one paper on data analysis of a new instrument to measure gases. The first contribution to this collection [6] reports the results of the processing and combination of high-rate and low-rate geodetic data for revealing the dynamics underly-ing violent volcanic eruptions at Mount Etna. This study evidences the wide spectrum of ground deformation produced by these phenomena, to be investigated, processed and modeled in order to generate a picture of the feeding system of the volcano and better un-derstand its dynamics and rates of magma transfer in the upper crust. Another contribution focuses on volcanoes [7]: the authors exploit 20 years of high temporal resolution satellite Thermal Infra-Red (TIR) data collected over three active vol-canoes (Etna, Shishaldin and Shinmoedake). They present the results of an analysis of this dataset performed through a preliminary RST (Robust Satellite Techniques) algorithm implementation to TIR data from the Advanced Spaceborne Thermal Emission and Re-flection Radiometer (ASTER). This approach ensures efficient identification and mapping of volcanic thermal features even of a low intensity level, which is also useful in the per-spective of an operational multi-satellite observing system. The contribution by Woohyun Son et al. [8] proposes specific depth-domain data processing of migration velocity analysis (MVA) of seismic data collected during a survey on a saline aquifer sediment in the Southern Continental Shelf of Korea. This analysis al-lowed the authors to identify and determine the precise depth of a basalt flow that could act as a cap rock for CO2 storage beneath the aquifer. The investigation, through the geo-logical model obtained from both time- and depth-domain processing, provides suitable information for locating the best drilling sites for CO2 injection, maximizing the storage volume. In volcanic areas, gases represent important physical evidence of volcanic processes that need to be measured. Parracino et al. [9] have shown how novel range-resolved DI-AL-Lidar (Differential Absorption Light Detection and Ranging) could herald a new era in the observation of long-term volcanic CO2 gases. An accurate and integrated analysis of different types of data such as GNSS, seismic and MT-InSAR, has led, in the work by Gatsios et al. [10], to a first account of deformation processes and their temporal evolution over recent years for Methana (Greece), thus providing initial information to feed into a volcano baseline hazard assessment and mon-itoring system. Seismic data are among the most important data to understand the dynamics of the Earth’s interior. A consistent analysis of a seismic swarm allowed Kostoglou et al. [11] to shed more light on the regional geodynamics of the Kefalonia Transform Fault Zone (Greece), and to follow the temporal evolution of the b-value to distinguish between fore-shock and aftershock behaviors.

Description: Published

Description: 10759

Description: 6SR VULCANI – Servizi e ricerca per la società

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1705-1730, https://doi.org/10.1175/jpo-d-21-0243.1.

Description: Formation and evolution of barrier layers (BLs) and associated temperature inversions (TIs) were investigated using a 1-yr time series of oceanic and air–sea surface observations from three moorings deployed in the eastern Pacific fresh pool. BL thickness and TI amplitude showed a seasonality with maxima in boreal summer and autumn when BLs were persistently present. Mixed layer salinity (MLS) and mixed layer temperature (MLT) budgets were constructed to investigate the formation mechanism of BLs and TIs. The MLS budget showed that BLs were initially formed in response to horizontal advection of freshwater in boreal summer and then primarily maintained by precipitation. The MLT budget revealed that penetration of shortwave radiation through the mixed layer base is the dominant contributor to TI formation through subsurface warming. Geostrophic advection is a secondary contributor to TI formation through surface cooling. When the BL exists, the cooling effect from entrainment and the warming effect from detrainment are both significantly reduced. In addition, when the BL is associated with the presence of a TI, entrainment works to warm the mixed layer. The presence of BLs makes the shallower mixed layer more sensitive to surface heat and freshwater fluxes, acting to enhance the formation of TIs that increase the subsurface warming via shortwave penetration.

Description: SK is supported by JSPS Overseas Research Fellowships. JS and SK are supported by NASA Grant 80NSSC18K1500. JTF and the mooring deployment were funded by NASA Grants NNX15AG20G and 80NSSC18K1494. DZ is supported by NASA Grant 80NSSC18K1499. This publication is partially funded by the Cooperative Institute for Climate, Ocean, and Ecosystem Studies (CICOES) under NOAA Cooperative Agreement NA20OAR4320271, Contribution 2021-1152. This is PMEL Contribution 5268.

Description: 2023-01-27

Description: Two paroxysmal explosions occurred at Stromboli on July 3 and August 28, 2019, the first of which caused the death of a young tourist. After the first paroxysm an effusive activity began from the summit vents and affected the NW flank of the island for the entire period between the two paroxysms. We carried out an unsupervised analysis of seismic and infrasonic data of Strombolian explosions over 10 months (15 November 2018–15 September 2019) using a Self- Organizing Map (SOM) neural network to recognize changes in the eruptive patterns of Stromboli that preceded the paroxysms. We used a dataset of 14,289 events. The SOM analysis identified three main clusters that showed different occurrences with time indicating a clear change in Stromboli’s eruptive style before the paroxysm of 3 July 2019. We compared the main clusters with the recordings of the fixed monitoring cameras and with the Ground-Based Interferometric Synthetic Aperture Radar measurements, and found that the clusters are associated with different types of Strombolian explosions and different deformation patterns of the summit area. Our findings provide new insights into Strombolian eruptive mechanisms and new perspectives to improve the monitoring of Stromboli and other open conduit volcanoes.

Description: This work was supported by the project Progetto Strategico Dipartimentale INGV 2019 “Forecasting eruptive activity at Stromboli volcano: timing, eruptive style, size, intensity and duration” (FIRST). This work is also supported by a Marie Sklodowska-Curie Innovative Training Network Fellowship of the European Commission’s Horizon 2020 Programme under Contract Number 765710 INSIGHTS. This work benefited from the EU (DG ECHO) Project EVE n. 826292. This work has been partially supported by the “Presidenza del Consiglio dei Ministri–Dipartimento della Protezione Civile” (Presidency of the Council of Ministers–Department of Civil Protection; Scientific Responsibility: N.C.). However, this publication does not necessarily represent the official opinion and policies of the department.

Description: Published

Description: 1287

Description: 4V. Processi pre-eruttivi

Description: JCR Journal

Description: The Etna volcano is renowned worldwide for its extraordinary lava fountains that rise several kilometers above the vent and feed eruptive columns, then drift hundreds of kilometers away from the source. The Italian Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE) is responsible for the monitoring of Mt. Etna, and for this reason, has deployed a network of visible and thermal cameras around the volcano. From these cameras, INGV-OE keeps a keen eye, and is able to observe the eruptive activity, promptly advising the civil protection and aviation authorities of any changes, as well as quantifying the spread of lava flows and the extent of pyroclastic and ash plumes by using a careful analysis of the videos recorded by the monitoring cameras. However, most of the work involves analysis carried out by hand, which is necessarily approximate and time-consuming, thus limiting the usefulness of these results for a prompt hazard assessment. In addition, the start of lava fountains is often a gradual process, increasing in strength from Strombolian activity, to intermediate explosive activity, and eventually leading to sustained lava fountains. The thresholds between these different fields (Strombolian, Intermediate, and lava fountains) are not clear cut, and are often very difficult to distinguish by a manual analysis of the images. In this paper, we presented an automated routine that, when applied to thermal images and with good weather conditions, allowed us to detect (1) the starting and ending time of each lava fountain, (2) the area occupied by hot pyroclasts, (3) the elevation reached by the lava fountains over time, and (4) eventually, to calculate in real-time the erupted volume of pyroclasts, giving results close to the manual analysis but more focused on the sustained portion of the lava fountain, which is also the most dangerous. This routine can also be applied to other active volcanoes, allowing a prompt and uniform definition of the timing of the lava fountain eruptive activity, as well as the magnitude and intensity of the event.

Description: This research was funded by the Project FIRST—ForecastIng eRuptive activity at Stromboli volcano: timing, eruptive style, size, intensity, and duration, INGV-Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020; Scientific Responsibility: S.C.).

Description: Published

Description: 2392

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Between 13 December 2020 and 21 February 2022, Etna volcano produced a sequence of 66 paroxysmal explosive eruptions, with Strombolian activity at the summit craters climaxing in lava fountains and eruption columns extending several kilometers above the craters, accompanied by minor and short-lasting lava flows from the crater rim. We selected three of these episodes that occurred within a short space of time, between 13 December 2020 and 12 March 2021, of different magnitude (i.e., erupted volume) and intensity (i.e., mass eruption rate or instantaneous eruption rate), and analyzed them from a multidisciplinary perspective. The aim was to gain insights into those parameters that mostly reveal the eruptive process for hazard assessment purposes. The multidisciplinary data consist of calibrated visible images, thermal images, seismic and infrasound data, ground deformation detected from the strainmeters, as well as satellite SEVIRI images. From these data, we obtained the timing of each paroxysmal event, the erupted volume in terms of tephra and lava flows, and the corresponding deflation of the source region, together with the development of the lava fountains and eruption columns with time. The results enabled determining that the smallest episode was that of 13 December 2020, which comprised three distinctive pulses but did not produce an eruptive column detectable from either monitoring cameras or satellites. The 28 February 2021 episode was remarkable for the short amount of time required to reach the climax, and was the most intense, whereas the 12 March 2021 event showed the longest duration but with an intensity between that of the previous two. Our results show that these three paroxysmal events display a typical trend, with the first event also being the smallest in terms of both erupted volume and intensity, the second being the most intense, and the third the one of greatest magnitude but less intense than the second. This is coherent with the end of the first paroxysmal phase on 1 April 2021, which was followed by 48 days of eruptive pause before starting again. In this context, the end of the paroxysmal phase was anticipated by a more effusive episode, thus heralding a temporary decline in the gas content within the feeding magma batch.

Description: This research was funded by the Project FIRST-ForecastIng eRuptive activity at Stromboli volcano: Timing, eruptive style, size, intensity, and duration; INGV-Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020). A.C. thanks the CHANCE project, II Edition, Università degli Studi di Catania (principal investigator A.C.) and the grant PIACERI, 2020-22 programme (PAROSSISMA project, code 22722132140; principal investigator Marco Viccaro). A.I. thanks the IMPACT project—A Multidisciplinary Insight on the Kinematics and Dynamics of Magmatic Processes at Mt. Etna Aimed at Identifying Rrecursor Phenomena and Developing Early Warning Systems, funded by INGV-Progetto Strategico Dipartimento Vulcani 2019 (Delibera n. 144/2020). S.S. thanks the ‘e-shape’ project, which receives funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement 820852.

Description: Published

Description: 4006

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(4), (2022): 597–616, https://doi.org/10.1175/jpo-d-21-0121.1.

Description: We provide a first-principles analysis of the energy fluxes in the oceanic internal wave field. The resulting formula is remarkably similar to the renowned phenomenological formula for the turbulent dissipation rate in the ocean, which is known as the finescale parameterization. The prediction is based on the wave turbulence theory of internal gravity waves and on a new methodology devised for the computation of the associated energy fluxes. In the standard spectral representation of the wave energy density, in the two-dimensional vertical wavenumber–frequency (m–ω) domain, the energy fluxes associated with the steady state are found to be directed downscale in both coordinates, closely matching the finescale parameterization formula in functional form and in magnitude. These energy transfers are composed of a “local” and a “scale-separated” contributions; while the former is quantified numerically, the latter is dominated by the induced diffusion process and is amenable to analytical treatment. Contrary to previous results indicating an inverse energy cascade from high frequency to low, at odds with observations, our analysis of all nonzero coefficients of the diffusion tensor predicts a direct energy cascade. Moreover, by the same analysis fundamental spectra that had been deemed “no-flux” solutions are reinstated to the status of “constant-downscale-flux” solutions. This is consequential for an understanding of energy fluxes, sources, and sinks that fits in the observational paradigm of the finescale parameterization, solving at once two long-standing paradoxes that had earned the name of “oceanic ultraviolet catastrophe.”

Description: The authors gratefully acknowledge support from the ONR Grant N00014-17-1-2852. YL gratefully acknowledges support from NSF DMS Award 2009418.

Description: 2022-09-25

Description: Editorial to a Special Issue

Description: The monitoring of active volcanoes is a complex task based on multidisciplinary and integrated analyses that use ground, drones, and satellite monitoring devices. Over time, and with the development of new technology and increasing frequency of acquisition, the use of remote sensing to accomplish this important task has grown enormously. This is especially so with the use of drones and satellites for classifying eruptive events, detecting the opening of new vents, the spreading of lava flows on the surface or ash plumes in the atmosphere, the fallout of tephra on the ground, the intrusion of new magma within the volcano edifice, and the deformation preceding impending eruptions, and others besides. The main challenge in using remote sensing techniques is to develop automated and reliable systems that may assist the decision-maker in volcano monitoring, hazard assessment, and risk reduction. The integration with ground-based techniques represents a valuable additional aspect that makes the proposed methods more robust and reinforces the results obtained. This collection of papers is focused on several active volcanoes, such as Stromboli, Etna, and Vulcano in Italy; the Long Valley caldera and Kilauea volcano in the USA; and Cotopaxi in Ecuador. The authors make use of several different methods to predict and forecast the volcanoes’ future behavior, using insights from the available data or from new automated routines applied to the analysis of existing data. The aim is to enable rapid assessments of the state of a volcano, discovering the connection between variables apparently not related to each other and to the state of the volcano. The development of new or automated routines is an important step forward in the process of forecasting eruptive activities, and this collection comprises several such examples. This Special Issue on the monitoring of active volcanoes using an integration of remote sensing and ground-based techniques comprises twelve papers. Three are focused on the results obtained for Stromboli volcano (Italy), where eruptive activity varies from moderate Strombolian, often accompanied by summit overflows, to highly explosive paroxysms, which are very dangerous both for the local population and for the many tourists who frequently visit the island. The first paper [1] presents the precursors of the paroxysmal and devastating explosive eruptions occurring in 2019. This paper applied an unsupervised analysis of seismic and infrasonic data, comprising a dataset of 14,289 Strombolian explosions occurring over 10 months, using a Self-Organizing Map (SOM) neural network to recognize changes in the eruptive patterns preceding the paroxysms. The SOM analysis identified three main clusters indicating a clear change in Stromboli’s eruptive style before the paroxysm of 3 July 2019. The main clusters were then compared with the recordings of the fixed monitoring cameras and with the Ground-Based Interferometric Synthetic Aperture Radar measurements, showing that they were associated with different types of Strombolian explosions and different deformation patterns of the summit area.

Description: Published

Description: 3626

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1593-1611, https://doi.org/10.1175/jpo-d-21-0180.1.

Description: This study presents novel observational estimates of turbulent dissipation and mixing in a standing meander between the Southeast Indian Ridge and the Macquarie Ridge in the Southern Ocean. By applying a finescale parameterization on the temperature, salinity, and velocity profiles collected from Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats in the upper 1600 m, we estimated the intensity and spatial distribution of dissipation rate and diapycnal mixing along the float tracks and investigated the sources. The indirect estimates indicate strong spatial and temporal variability of turbulent mixing varying from O(10−6) to O(10−3) m2 s−1 in the upper 1600 m. Elevated turbulent mixing is mostly associated with the Subantarctic Front (SAF) and mesoscale eddies. In the upper 500 m, enhanced mixing is associated with downward-propagating wind-generated near-inertial waves as well as the interaction between cyclonic eddies and upward-propagating internal waves. In the study region, the local topography does not play a role in turbulent mixing in the upper part of the water column, which has similar values in profiles over rough and smooth topography. However, both remotely generated internal tides and lee waves could contribute to the upward-propagating energy. Our results point strongly to the generation of turbulent mixing through the interaction of internal waves and the intense mesoscale eddy field.

Description: The observations were funded through grants from the Australian Research Council Discovery Project (DP170102162) and Australia’s Marine National Facility. Surface drifters were provided by Dr. Shaun Dolk of the Global Drifter Program. AC was supported by an Australian Research Council Postdoctoral Fellowship. AC, HEP, and NLB acknowledge support from the Australian Government Department of the Environment and Energy National Environmental Science Program and the ARC Centre of Excellence in Climate Extremes. KP acknowledges the support from the National Science Foundation.

Description: This work is devoted to the analysis of the background seismic noise acquired at the volcanoes (Campi Flegrei caldera, Ischia island, and Vesuvius) belonging to the Neapolitan volcanic district (Italy), and at the Colima volcano (Mexico). Continuous seismic acquisition is a complex mixture of volcanic transients and persistent volcanic and/or hydrothermal tremor, anthropogenic/ambient noise, oceanic loading, and meteo-marine contributions. The analysis of the background noise in a stationary volcanic phase could facilitate the identification of relevant waveforms often masked by microseisms and ambient noise. To address this issue, our approach proposes a machine learning (ML) modeling to recognize the “fingerprint” of a specific volcano by analyzing the background seismic noise from the continuous seismic acquisition. Specifically, two ML models, namely multi-layer perceptrons and convolutional neural network were trained to recognize one volcano from another based on the acquisition noise. Experimental results demonstrate the effectiveness of the two models in recognizing the noisy background signal, with promising performance in terms of accuracy, precision, recall, and F1 score. These results suggest that persistent volcanic signals share the same source information, as well as transient events, revealing a common generation mechanism but in different regimes. Moreover, assessing the dynamic state of a volcano through its background noise and promptly identifying any anomalies, which may indicate a change in its dynamics, can be a practical tool for real-time monitoring.

Description: Published

Description: 6835

Description: 4V. Processi pre-eruttivi

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(2), (2022): 223–235, https://doi.org/10.1175/JTECH-D-21-0110.1.

Description: Previous work with simulations of oceanographic high-frequency (HF) radars has identified possible improvements when using maximum likelihood estimation (MLE) for direction of arrival; however, methods for determining the number of emitters (here defined as spatially distinct patches of the ocean surface) have not realized these improvements. Here we describe and evaluate the use of the likelihood ratio (LR) for emitter detection, demonstrating its application to oceanographic HF radar data. The combined detection–estimation methods MLE-LR are compared with multiple signal classification method (MUSIC) and MUSIC parameters for SeaSonde HF radars, along with a method developed for 8-channel systems known as MUSIC-Highest. Results show that the use of MLE-LR produces similar accuracy, in terms of the RMS difference and correlation coefficients squared, as previous methods. We demonstrate that improved accuracy can be obtained for both methods, at the cost of fewer velocity observations and decreased spatial coverage. For SeaSondes, accuracy improvements are obtained with less commonly used parameter sets. The MLE-LR is shown to be able to resolve simultaneous closely spaced emitters, which has the potential to improve observations obtained by HF radars operating in complex current environments.

Description: This work was supported by the National Science Foundation (NSF) under Grant OCE-1658475. Computing resources were provided by the UCSB Center for Scientific Computing through an NSF MRSEC (DMR-1720256) and NSF CNS-1725797.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(1),(2022): 75–97, https://doi.org/10.1175/JPO-D-21-0099.1.

Description: Mesoscale eddies contain the bulk of the ocean’s kinetic energy (KE), but fundamental questions remain on the cross-scale KE transfers linking eddy generation and dissipation. The role of submesoscale flows represents the key point of discussion, with contrasting views of submesoscales as either a source or a sink of mesoscale KE. Here, the first observational assessment of the annual cycle of the KE transfer between mesoscale and submesoscale motions is performed in the upper layers of a typical open-ocean region. Although these diagnostics have marginal statistical significance and should be regarded cautiously, they are physically plausible and can provide a valuable benchmark for model evaluation. The cross-scale KE transfer exhibits two distinct stages, whereby submesoscales energize mesoscales in winter and drain mesoscales in spring. Despite this seasonal reversal, an inverse KE cascade operates throughout the year across much of the mesoscale range. Our results are not incompatible with recent modeling investigations that place the headwaters of the inverse KE cascade at the submesoscale, and that rationalize the seasonality of mesoscale KE as an inverse cascade-mediated response to the generation of submesoscales in winter. However, our findings may challenge those investigations by suggesting that, in spring, a downscale KE transfer could dampen the inverse KE cascade. An exploratory appraisal of the dynamics governing mesoscale–submesoscale KE exchanges suggests that the upscale KE transfer in winter is underpinned by mixed layer baroclinic instabilities, and that the downscale KE transfer in spring is associated with frontogenesis. Current submesoscale-permitting ocean models may substantially understate this downscale KE transfer, due to the models’ muted representation of frontogenesis.

Description: The OSMOSIS experiment was funded by the U.K. Natural Environment Research Council (NERC) through Grants NE/1019999/1 and NE/101993X/1. ACNG acknowledges the support of the Royal Society and the Wolfson Foundation, and XY that of a China Scholarship Council PhD studentship.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1677-1691, https://doi.org/10.1175/jpo-d-21-0269.1.

Description: Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed.

Description: This work was generously funded by NSF Grants OCE-1912302, OCE-1912125 (Drushka), and OCE-1912325 (Abernathey) as part of the Ocean Energy and Eddy Transport Climate Process Team.

Description: Author Posting. © American Meteorological Society , 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Zippel, S. F., Farrar, J. T., Zappa, C. J., Miller, U., St Laurent, L., Ijichi, T., Weller, R. A., McRaven, L., Nylund, S., & Le Bel, D. Moored turbulence measurements using pulse-coherent doppler sonar. Journal of Atmospheric and Oceanic Technology, 38(9), (2021): 1621–1639, https://doi.org/10.1175/JTECH-D-21-0005.1.

Description: Upper-ocean turbulence is central to the exchanges of heat, momentum, and gases across the air–sea interface and therefore plays a large role in weather and climate. Current understanding of upper-ocean mixing is lacking, often leading models to misrepresent mixed layer depths and sea surface temperature. In part, progress has been limited by the difficulty of measuring turbulence from fixed moorings that can simultaneously measure surface fluxes and upper-ocean stratification over long time periods. Here we introduce a direct wavenumber method for measuring turbulent kinetic energy (TKE) dissipation rates ϵ from long-enduring moorings using pulse-coherent ADCPs. We discuss optimal programming of the ADCPs, a robust mechanical design for use on a mooring to maximize data return, and data processing techniques including phase-ambiguity unwrapping, spectral analysis, and a correction for instrument response. The method was used in the Salinity Processes Upper-Ocean Regional Study (SPURS) to collect two year-long datasets. We find that the mooring-derived TKE dissipation rates compare favorably to estimates made nearby from a microstructure shear probe mounted to a glider during its two separate 2-week missions for O(10−8) ≤ ϵ ≤ O(10−5) m2 s−3. Periods of disagreement between turbulence estimates from the two platforms coincide with differences in vertical temperature profiles, which may indicate that barrier layers can substantially modulate upper-ocean turbulence over horizontal scales of 1–10 km. We also find that dissipation estimates from two different moorings at 12.5 and at 7 m are in agreement with the surface buoyancy flux during periods of strong nighttime convection, consistent with classic boundary layer theory.

Description: This work was funded by NASA as part of the Salinity Processes in the Upper Ocean Regional Study (SPURS), supporting field work for SPURS-1 (NASA Grant NNX11AE84G), for SPURS-2 (NASA Grant NNX15AG20G), and for analysis (NASA Grant 80NSSC18K1494). Funding for early iterations of this project associated with the VOCALS project and Stratus 9 mooring was provided by NSF (Awards 0745508 and 0745442). Additional funding was provided by ONR Grant N000141812431 and NSF Award 1756839. The Stratus Ocean Reference Station is funded by the Global Ocean Monitoring and Observing Program of the National Oceanic and Atmospheric Administration (CPO FundRef Number 100007298), through the Cooperative Institute for the North Atlantic Region (CINAR) under Cooperative Agreement NA14OAR4320158. Microstructure measurements made from the glider were supported by NSF (Award 1129646).

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(3), (2022): 363–382, https://doi.org/10.1175/jpo-d-21-0084.1.

Description: Meltwater from Greenland is an important freshwater source for the North Atlantic Ocean, released into the ocean at the head of fjords in the form of runoff, submarine melt, and icebergs. The meltwater release gives rise to complex in-fjord transformations that result in its dilution through mixing with other water masses. The transformed waters, which contain the meltwater, are exported from the fjords as a new water mass Glacially Modified Water (GMW). Here we use summer hydrographic data collected from 2013 to 2019 in Upernavik, a major glacial fjord in northwest Greenland, to describe the water masses that flow into the fjord from the shelf and the exported GMWs. Using an optimum multi-parameter technique across multiple years we then show that GMW is composed of 57.8% ± 8.1% Atlantic Water (AW), 41.0% ± 8.3% Polar Water (PW), 1.0% ± 0.1% subglacial discharge, and 0.2% ± 0.2% submarine meltwater. We show that the GMW fractional composition cannot be described by buoyant plume theory alone since it includes lateral mixing within the upper layers of the fjord not accounted for by buoyant plume dynamics. Consistent with its composition, we find that changes in GMW properties reflect changes in the AW and PW source waters. Using the obtained dilution ratios, this study suggests that the exchange across the fjord mouth during summer is on the order of 50 mSv (1 Sv ≡ 106 m3 s−1) (compared to a freshwater input of 0.5 mSv). This study provides a first-order parameterization for the exchange at the mouth of glacial fjords for large-scale ocean models.

Description: This work was partially supported by the Centre for Climate Dynamics (SKD) at the Bjerknes Centre for Climate Research. The authors thank NASA and the OMG consortium for making observational data freely available, and acknowledge M. Morlighem for good support in the early stages of this project. MM and LHS and would also like to thank Ø. Paasche, the ACER project, and the U.S. Norway Fulbright Foundation for the Norwegian Arctic Chair Grant 2019–20 that made the visit to Scripps Institution of Oceanography possible. FS acknowledges support from the DOE Office of Science Grant DE-SC0020073, Heising-Simons Foundation and from NSF and OCE-1756272. DAS acknowledges support from U.K. NERC Grants NE/P011365/1, NE/T011920/1, and NERC Independent Research Fellowship NE/T011920/1. MW was supported by an appointment to the NASA Postdoctoral Program at the Jet Propulsion Laboratory, California Institute of Technology, administered by the Universities Space Research Association under contract with NASA. CSA would like to acknowledge Geocenter Denmark for support to the project “Upernavik Glacier.”

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(12),(2021): 3663–3678, https://doi.org/10.1175/JPO-D-21-0058.1.

Description: The wind-driven exchange through complex ridges and islands between marginal seas and the open ocean is studied using both numerical and analytical models. The models are forced by a steady, spatially uniform northward wind stress intended to represent the large-scale, low-frequency wind patterns typical of the seasonal monsoons in the western Pacific Ocean. There is an eastward surface Ekman transport out of the marginal sea and westward geostrophic inflows into the marginal sea. The interaction between the Ekman transport and an island chain produces strong baroclinic flows along the island boundaries with a vertical depth that scales with the ratio of the inertial boundary layer thickness to the baroclinic deformation radius. The throughflows in the gaps are characterized by maximum transport in the center gap and decreasing transports toward the southern and northern tips of the island chain. An extended island rule theory demonstrates that throughflows are determined by the collective balance between viscosity on the meridional boundaries and the eastern side boundary of the islands. The outflowing transport is balanced primarily by a shallow current that enters the marginal sea along its equatorward boundary. The islands can block some direct exchange and result in a wind-driven overturning cell within the marginal sea, but this is compensated for by eastward zonal jets around the southern and northern tips of the island chain. Topography in the form of a deep slope, a ridge, or shallow shelves around the islands alters the current pathways but ultimately is unable to limit the total wind-driven exchange between the marginal sea and the open ocean.

Description: This research is supported in part by the China Scholarship Council (201906330102). H. G. is financially supported by the China Scholarship Council to study at WHOI for 2 years as a guest student. M. A. S. is supported by the National Science Foundation Grant OCE-1922538.

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Fine, E., MacKinnon, J., Alford, M., Middleton, L., Taylor, J., Mickett, J., Cole, S., Couto, N., Boyer, A., & Peacock, T. Double diffusion, shear instabilities, and heat impacts of a pacific summer water intrusion in the Beaufort Sea. Journal of Physical Oceanography, 52(2), (2022): 189–203, https://doi.org/10.1175/jpo-d-21-0074.1.

Description: Pacific Summer Water eddies and intrusions transport heat and salt from boundary regions into the western Arctic basin. Here we examine concurrent effects of lateral stirring and vertical mixing using microstructure data collected within a Pacific Summer Water intrusion with a length scale of ∼20 km. This intrusion was characterized by complex thermohaline structure in which warm Pacific Summer Water interleaved in alternating layers of O(1) m thickness with cooler water, due to lateral stirring and intrusive processes. Along interfaces between warm/salty and cold/freshwater masses, the density ratio was favorable to double-diffusive processes. The rate of dissipation of turbulent kinetic energy (ε) was elevated along the interleaving surfaces, with values up to 3 × 10−8 W kg−1 compared to background ε of less than 10−9 W kg−1. Based on the distribution of ε as a function of density ratio Rρ, we conclude that double-diffusive convection is largely responsible for the elevated ε observed over the survey. The lateral processes that created the layered thermohaline structure resulted in vertical thermohaline gradients susceptible to double-diffusive convection, resulting in upward vertical heat fluxes. Bulk vertical heat fluxes above the intrusion are estimated in the range of 0.2–1 W m−2, with the localized flux above the uppermost warm layer elevated to 2–10 W m−2. Lateral fluxes are much larger, estimated between 1000 and 5000 W m−2, and set an overall decay rate for the intrusion of 1–5 years.

Description: This work was supported by ONR Grant N00014-16-1-2378 and NSF Grants PLR 14-56705 and PLR-1303791, NSF Graduate Research Fellowship Grant DGE-1650112, as well as by the Postdoctoral Scholar Program at Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(4), (2022): 491–502, https://doi.org/10.1175/jtech-d-21-0046.1.

Description: The Air-Launched Autonomous Micro Observer (ALAMO) is a versatile profiling float that can be launched from an aircraft to make temperature and salinity observations of the upper ocean for over a year with high temporal sampling. Similar in dimensions and weight to an airborne expendable bathythermograph (AXBT), but with the same capability as Argo profiling floats, ALAMOs can be deployed from an A-sized (sonobuoy) launch tube, the stern ramp of a cargo plane, or the door of a small aircraft. Unlike an AXBT, however, the ALAMO float directly measures pressure, can incorporate additional sensors, and is capable of performing hundreds of ocean profiles compared to the single temperature profile provided by an AXBT. Upon deployment, the float parachutes to the ocean, releases the air-deployment package, and immediately begins profiling. Ocean profile data along with position and engineering information are transmitted via the Iridium satellite network, automatically processed, and then distributed by the Global Telecommunications System for use by the operational forecasting community. The ALAMO profiling mission can be modified using the two-way Iridium communications to change the profiling frequency and depth. Example observations are included to demonstrate the ALAMO’s utility.

Description: This work was supported by the National Oceanographic and Atmospheric Administration under Grants NA13OAR4830233 (as part of CINAR Sandy Supplemental funding from the Disaster Relief Appropriations Act of 2013) and NA14OAR4320158 and by Office of Naval Research under Grants N0001416WX01384, N0001416WX01262, and N000141512293. ALAMO floats are commercially available from MRV Systems, LLC (https://www.mrvsys.com).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of [publisher] for personal use, not for redistribution. The definitive version was published in Farrar, J. T., Durland, T., Jayne, S. R., & Price, J. F. Long-distance radiation of Rossby Waves from the equatorial current system. Journal of Physical Oceanography, 51(6), (2021): 1947–1966, https://doi.org/10.1175/JPO-D-20-0048.1.

Description: Measurements from satellite altimetry are used to show that sea surface height (SSH) variability throughout much of the North Pacific Ocean is coherent with the SSH signal of the tropical instability waves (TIWs) that result from instabilities of the equatorial currents. This variability has regular phase patterns consistent with freely propagating barotropic Rossby waves radiating energy away from the unstable equatorial currents, and the waves clearly propagate from the equatorial region to at least 30°N. The pattern of SSH variance at TIW frequencies exhibits remarkable patchiness on scales of hundreds of kilometers, which we interpret as being due to the combined effects of wave reflection, refraction, and interference. North of 40°N, more than 6000 km from the unstable equatorial currents, the SSH field remains coherent with the near-equatorial SSH variability, but it is not as clear whether the variability at the higher latitudes is a simple result of barotropic wave radiation from the tropical instability waves. Even more distant regions, as far north as the Aleutian Islands off of Alaska and the Kamchatka Peninsula of eastern Russia, have SSH variability that is significantly coherent with the near-equatorial instabilities. The variability is not well represented in the widely used gridded SSH data product commonly referred to as the AVISO or DUACS product, and this appears to be a result of spatial variations in the filtering properties of the objective mapping scheme.

Description: This work was supported by NASA Grants NNX13AE46G, NNX14AM71G, and NNX17AH54G.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(8), (2022): 1183-1198, https://doi.org/10.1175/jtech-d-21-0068.1.

Description: Horizontal kinematic properties, such as vorticity, divergence, and lateral strain rate, are estimated from drifter clusters using three approaches. At submesoscale horizontal length scales O(1–10)km, kinematic properties become as large as planetary vorticity f, but challenging to observe because they evolve on short time scales O(hourstodays). By simulating surface drifters in a model flow field, we quantify the sources of uncertainty in the kinematic property calculations due to the deformation of cluster shape. Uncertainties arise primarily due to (i) violation of the linear estimation methods and (ii) aliasing of unresolved scales. Systematic uncertainties (iii) due to GPS errors, are secondary but can become as large as (i) and (ii) when aspect ratios are small. Ideal cluster parameters (number of drifters, length scale, and aspect ratio) are determined and error functions estimated empirically and theoretically. The most robust method—a two-dimensional, linear least squares fit—is applied to the first few days of a drifter dataset from the Bay of Bengal. Application of the length scale and aspect-ratio criteria minimizes errors (i) and (ii), and reduces the total number of clusters and so computational cost. The drifter-estimated kinematic properties map out a cyclonic mesoscale eddy with a surface, submesoscale fronts at its perimeter. Our analyses suggest methodological guidance for computing the two-dimensional kinematic properties in submesoscale flows, given the recently increasing quantity and quality of drifter observations, while also highlighting challenges and limitations.

Description: This research was supported by the Office of Naval Research (ONR) Departmental Research Initiative ASIRI under Grant N00014-13-1-0451 (SE and AM) and Grant N00014-13-1-0477 (VH and LC). The authors thank the captain and crew of the R/V Roger Revelle, and Andrew Lucas with the Multiscale Ocean Dynamics group at the Scripps Institution for Oceanography for providing the FastCTD data collected in 2015, which was supported by ONR Grant N00014-13-1-0489, as well as Eric D’Asaro for helpful discussions and Lance Braasch for assistance with the drifter dataset. AM and SE further thank NSF (Grant OCE-I434788) and ONR (Grant N00014-16-1-2470) for support. VH and LC were additionally supported by ONR Grants N00014-15-1-2286, N00014-14-1-0183, N00014-19-1-26-91 and NOAA Global Drifter Program (GDP) Grant NA15OAR4320071.

Description: 2023-02-01

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 35(17), (2022): 5465-5482, https://doi.org/10.1175/jcli-d-21-0671.1.

Description: Understanding the contribution of ocean circulation to glacial–interglacial climate change is a major focus of paleoceanography. Specifically, many have tried to determine whether the volumes and depths of Antarctic- and North Atlantic–sourced waters in the deep ocean changed at the Last Glacial Maximum (LGM; ∼22–18 kyr BP) when atmospheric pCO2 concentrations were 100 ppm lower than the preindustrial. Measurements of sedimentary geochemical proxies are the primary way that these deep ocean structural changes have been reconstructed. However, the main proxies used to reconstruct LGM Atlantic water mass geometry provide conflicting results as to whether North Atlantic–sourced waters shoaled during the LGM. Despite this, a number of idealized modeling studies have been advanced to describe the physical processes resulting in shoaled North Atlantic waters. This paper aims to critically assess the approaches used to determine LGM Atlantic circulation geometry and lay out best practices for future work. We first compile existing proxy data and paleoclimate model output to deduce the processes responsible for setting the ocean distributions of geochemical proxies in the LGM Atlantic Ocean. We highlight how small-scale mixing processes in the ocean interior can decouple tracer distributions from the large-scale circulation, complicating the straightforward interpretation of geochemical tracers as proxies for water mass structure. Finally, we outline promising paths toward ascertaining the LGM circulation structure more clearly and deeply.

Description: S.K.H. was supported by the Investment in Science Fund at WHOI and the John E. and Anne W. Sawyer Endowed Fund in Support of Scientific Staff. F.J.P. was supported by a Stanback Postdoctoral Fellowship at Caltech.

Description: 2023-02-01

Description: Data availability: The source data (ground temperature measurements recorded -hourly- from 2009 to 2012) are submitted to PANGAEA — PDI31617.

Description: Mild thermal anomalies are sensitive to changes in the advection processes in a volcanic system. A mild thermal anomaly, near the top of the North-East Rift of Mt. Etna (Italy), has been monitored from January 2010 to September 2012 by means of four temperature sensors buried in the shallow ground. The pulses of the convective circulation have been tracked and the diffuse heat flux has been evaluated. The positive pulses of the convective front reflected the local increases of volcanic degassing; conversely, the negative pulses showed the contraction of the convective front emerging through the North-East Rift. The steam condensation depth fluctuated below the monitoring site, from depths of a couple of meters to more than 30 meters, while the New South-East crater was erupting. The data hourly recorded, relative to the 2012 eruptive period, were compared to the radiant energy released by the paroxysms. We registered a dramatic decrease in the diffuse heat flux several hours before the onset of the two most energetic paroxysms (12 and 23 April). Thereafter, the convective front (the steam condensation depth) showed many negative pulses, reaching the deepest recorded levels. Thermal transients could be one of the early signals, possibly heralding transitions in the dynamic equilibrium conditions.

Description: Published

Description: 4471

Description: 2V. Struttura e sistema di alimentazione dei vulcani

Description: 4V. Processi pre-eruttivi

Description: 5V. Processi eruttivi e post-eruttivi

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: 6A. Geochimica per l'ambiente e geologia medica

Description: 1IT. Reti di monitoraggio e sorveglianza

Description: 6IT. Osservatori non satellitari

Description: JCR Journal

Description: Following the 2004 seismic unrest at Tenerife and the 2011–2012 submarine eruption at El Hierro, the number of Global Navigation Satellite System (GNSS) observation sites in the Canary Islands (Spain) has increased, offering scientists a useful tool with which to infer the kinematics and present-day surface deformation of the Canary sector of the Atlantic Ocean. We take advantage of the common-mode component filtering technique to improve the signal-to-noise ratio of the velocities retrieved from the daily solutions of 18 permanent GNSS stations distributed in the Canaries. The analysis of GNSS time series spanning the period 2011–2017 enabled us to characterize major regions of deformation along the archipelago through the mapping of the 2D infinitesimal strain field. By applying the triangular segmentation approach to GNSS velocities, we unveil a variable kinematic behaviour within the islands. The retrieved extension pattern shows areas of maximum deformation west of Tenerife, Gran Canaria and Fuerteventura. For the submarine main seismogenic fault between Tenerife and Gran Canaria, we simulated the horizontal deformation and strain due to one of the strongest (mbLg 5.2) earthquakes of the region. The seismic areas between islands, mainly offshore Tenerife and Gran Canaria, seem mainly influenced by the regional tectonic stress, not the local volcanic activity. In addition, the analysis of the maximum shear strain confirms that the regional stress field influences the E–W and NE–SW tectonic lineaments, which, in accordance with the extensional and compressional tectonic regimes identified, might favour episodes of volcanism in the Canary Islands.

Description: Published

Description: 3297

Description: 2T. Deformazione crostale attiva

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: In the presented work, the spectral emissivity of basaltic melts at magmatic temperatures was retrieved in a laboratory-controlled experiment by measuring their spectral radiance. Granulated bombs of Etnean basalts were melted and the radiant energy from the melting surface was recorded by a portable spectroradiometer in the short wavelength infrared (SWIR) spectral range between 1500 and 2500 nm. The Draping algorithm, an improved algorithm for temperature and emissivity separation, was applied for the first time to SWIR hyperspectral data in order to take into account the non-uniform temperature distribution of the melt surface and, at the same time, solving the two temperatures and the spectral emissivity. The results have been validated by comparing our results with the emissivity measured at a "lava simulator". Basalt spectral emissivity does not vary significantly at magmatic temperature, but shows an absorption feature in the range 2180–2290 nm, an atmospheric window pivotal for the IR remote sensing of active volcanoes

Description: Published

Description: 2046

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Explosive basaltic eruptions eject a great amount of pyroclastic material into the atmosphere, forming columns rising to several kilometers above the eruptive vent and causing significant disruption to both proximal and distal communities. Here, we analyze data, collected by an X-band polarimetric weather radar and an L-band Doppler fixed-pointing radar, as well as by a thermal infrared (TIR) camera, in relation to lava fountain-fed tephra plumes at the Etna volcano in Italy. We clearly identify a jet, mainly composed of lapilli and bombs mixed with hot gas in the first portion of these volcanic plumes and here called the incandescent jet region (IJR). At Etna and due to the TIR camera configuration, the IJR typically corresponds to the region that saturates thermal images. We find that the IJR is correlated to a unique signature in polarimetric radar data as it represents a zone with a relatively high reflectivity and a low copolar correlation coe cient. Analyzing five recent Etna eruptions occurring in 2013 and 2015, we propose a jet region radar retrieval algorithm (JR3A), based on a decision-tree combining polarimetric X-band observables with L-band radar constraints, aiming at the IJR height detection during the explosive eruptions. The height of the IJR does not exactly correspond to the height of the lava fountain due to a di erent altitude, potentially reached by lapilli and blocks detected by the X-band weather radar. Nonetheless, it can be used as a proxy of the lava fountain height in order to obtain a first approximation of the exit velocity of the mixture and, therefore, of the mass eruption rate. The comparisons between the JR3A estimates of IJR heights with the corresponding values recovered from TIR imagery, show a fairly good agreement with di erences of less than 20% in clear air conditions, whereas the di erence between JR3A estimates of IJR height values and those derived from L-band radar data only are greater than 40%. The advantage of using an X-band polarimetric weather radar in an early warning system is that it provides information in all weather conditions. As a matter of fact, we show that JR3A retrievals can also be obtained in cloudy conditions when the TIR camera data cannot be processed.

Description: Published

Description: 3629

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: On the morning of 24 December 2018, an eruptive event occurred at Etna, which was followed the next day by a strong sequence of shallow earthquakes. The eruptive episode lasted until 30 December, ranging from moderate strombolian to lava fountain activity coupled with vigorous ash/gas emissions and a lava flow e usion toward the eastern volcano flank of Valle del Bove. In this work, the data collected from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) instruments on board the Meteosat Second Generation (MSG) geostationary satellite are used to characterize the Etna activity by estimating the proximal and distal eruption parameters in near real time. The inversion of data indicates the onset of eruption on 24 December at 11:15 UTC, a maximum Time Average Discharge Rate (TADR) of 8.3 m3/s, a cumulative lava volume emitted of 0.5 Mm3, and a Volcanic Plume Top Height (VPTH) that reached a maximum altitude of 8 km above sea level (asl). The volcanic cloud ash and SO2 result totally collocated, with an ash amount generally lower than SO2 except on 24 December during the climax phase. A total amount of about 100 and 35 kt of SO2 and ash respectively was emitted during the entire eruptive period, while the SO2 fluxes reached peaks of more than 600 kg/s, with a mean value of about 185 kg/s. The SEVIRI VPTH, ash/SO2 masses, and flux time series have been compared with the results obtained from the ground-based visible (VIS) cameras and FLux Automatic MEasurements (FLAME) networks, and the satellite images collected by the MODerate resolution Imaging Spectroradiometer (MODIS) instruments on board the Terra and Aqua- polar satellites. The analysis indicates good agreement between SEVIRI, VIS camera, and MODIS retrievals with VPTH, ash, and SO2 estimations all within measurement errors. The SEVIRI and FLAME SO2 flux retrievals show significant discrepancies due to the presence of volcanic ash and a gap of data on the FLAME network. The results obtained in this study show the ability of geostationary satellite systems to characterize eruptive events from the source to the atmosphere in near real time during the day and night, thus o ering a powerful tool to mitigate volcanic risk on both local population and airspace and to give insight on volcanic processes.

Description: Published

Description: 1336

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Here we explore and review some of the latest ideas and applications of Raman spectroscopy to the volcanological sciences. Firstly, we provide a brief overview of how Raman spectral analysis works and how spectra from silicate glasses are interpreted. We then look at specific applications of Raman spectral analysis to the volcanological sciences based on measurements on and studies of natural materials in the laboratory. We conclude by examining the potential for Raman spectral analysis to be used as a field based aid to volcano monitoring via in situ studies of proximal deposits and; perhaps; in remote sensing campaigns.

Description: Published

Description: 805

Description: 3V. Proprietà chimico-fisiche dei magmi e dei prodotti vulcanici

Description: JCR Journal

Description: This work has been included in the discussion presented by the author at the 7th International Conference on Time Series and Forecasting, Gran Canaria, Spain,19–21 July 2021. Data Availability Statement: https://doi.org/10.26022/IEDA/112021, accessed on 11 June 2021.

Description: On the Island of Vulcano (Aeolian Archipelago, Italy) the temperatures of fumarole emissions, have ranged from about 700 ◦C to the boiling point. Since the end of the last eruption (1890 A.D.), many periods of increased heating of hydrothermal systems, underlying the La Fossa area have been identified, but an eruptive condition has not yet been reached. The time variation of the high temperature fumaroles has been tracked by the network of sensors located at a few discrete sites on the summit area of La Fossa cone. The same continuous monitoring network has been working for more than 30 years. The time series shows that a natural cyclic modulation has repeated after almost 20 years, and its periodicity yet has to be discussed and interpreted. The statistical approach and the spectral analysis could provide an objective evaluation to reveal the timing, intensity, and general significance of the thermodynamic perturbations that occurred in the hydrothermal circuits of La Fossa caldera, during the study period. The continuous monitoring data series avoid unrealistic interpolations and allow promptly recognizing changes, which perturb the hydrothermal circuits, highlighting—possibly in near real time—the transient phases of energy release from the different sources (hydrologic/magmatic).

Description: agreement between Istituto Nazionale di Geofisica e Vulcanologia and the Italian Presidenza del Consiglio dei Ministri, Dipartimento della Protezione Civile (DPC).

Description: Published

Description: 47

Description: 4V. Processi pre-eruttivi

Description: N/A or not JCR

Description: Volcanic emissions are a well-known hazard that can have serious impacts on local populations and aviation operations. Whereas several remote sensing observations detect high-intensity explosive eruptions, few studies focus on low intensity and long-lasting volcanic emissions. In this work, we have managed to fully characterize those events by analyzing the volcanic plume produced on the last day of the 2018 Christmas eruption at Mt. Etna, in Italy. We combined data from a visible calibrated camera, a multi-wavelength elastic/Raman Lidar system, from SEVIRI (EUMETSAT-MSG) and MODIS (NASA-Terra/Aqua) satellites and, for the first time, data from an automatic sun-photometer of the aerosol robotic network (AERONET). Results show that the volcanic plume height, ranging between 4.5 and 6 km at the source, decreased by about 0.5 km after 25 km. Moreover, the volcanic plume was detectable by the satellites up to a distance of about 400 km and contained very fine particles with a mean e ective radius of about 7 m. In some time intervals, volcanic ash mass concentration values were around the aviation safety thresholds of 2 10􀀀3 g m􀀀3. Of note, Lidar observations show two main stratifications of about 0.25 km, which were not observed at the volcanic source. The presence of the double stratification could have important implications on satellite retrievals, which usually consider only one plume layer. This work gives new details on the main features of volcanic plumes produced during low intensity and long-lasting volcanic plume emissions.

Description: Published

Description: 3866

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Volcanism is the most widespread expression of cyclic processes of formation and/or destruction that shape the Earth’s surface. Calderas are morphological depressions resulting from the collapse of a magma chamber following large eruptions and are commonly found in subduction-related tectono-magmatic regimes, such as arc and back-arc settings. Some of the most impressive examples of seafloor hydrothermal venting occur within submarine calderas. Here, we show the results of magnetic investigations at two hydrothermally active submarine calderas, i.e., Palinuro Seamount in the Southern Tyrrhenian Sea, Italy, and Brothers volcano of the Kermadec arc, New Zealand. These volcanoes occur in different geodynamic settings but show similarities in the development of their hydrothermal systems, both of which are hosted within calderas. We present a new integrated model based on morphological, geological and magnetic data for the Palinuro caldera, and we compare this with the well-established model of Brothers caldera, highlighting the differences and common features in the geophysical expressions of both hydrothermal systems. For consistency with the results at Brothers volcano, we build a model of demagnetised areas associated with hydrothermal alteration derived from 3D inversion of magnetic data. Both these models for Brothers and Palinuro show that hydrothermal up-flow zones are strongly controlled by caldera structures which provide large-scale permeability pathways, favouring circulation of the hydrothermal fluids at depth.

Description: Published

Description: 504

Description: 2V. Struttura e sistema di alimentazione dei vulcani

Description: 1A. Geomagnetismo e Paleomagnetismo

Description: 3A. Geofisica marina e osservazioni multiparametriche a fondo mare

Description: JCR Journal

Description: In July and August 2019, Stromboli volcano underwent two dangerous paroxysms previously considered “unexpected” because of the absence of significant changes in usually monitored parameters. We applied a multidisciplinary approach to search for signals able to indicate the possibility of larger explosive activity and to devise a model to explain the observed variations. We analysed geodetic data, satellite thermal data, images from remote cameras and seismic data in a timespan crossing the eruptive period of 2019 to identify precursors of the two paroxysms on a medium-term time span (months) and to perform an in-depth analysis of the signals recorded on a short time scale (hours, minutes) before the paroxysm. We developed a model that explains the observations. We call the model “push and go” where the uppermost feeding system of Stromboli is made up of a lower section occupied by a low viscosity, low density magma that is largely composed of gases and a shallower section occupied by the accumulated melt. We hypothesize that the paroxysms are triggered when an overpressure in the lower section is built up; the explosion will occur at the very moment such overpressure overcomes the confining pressure of the highly viscous magma above it.

Description: Published

Description: 4064

Description: 4V. Processi pre-eruttivi

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(9), (2020): 2491-2506, doi:10.1175/JPO-D-20-0056.1.

Description: An idealized two-layer shallow water model is applied to the study of the dynamics of the Arctic Ocean halocline. The model is forced by a surface stress distribution reflective of the observed wind stress pattern and ice motion and by an inflow representing the flow of Pacific Water through Bering Strait. The model reproduces the main elements of the halocline circulation: an anticyclonic Beaufort Gyre in the western basin (representing the Canada Basin), a cyclonic circulation in the eastern basin (representing the Eurasian Basin), and a Transpolar Drift between the two gyres directed from the upwind side of the basin to the downwind side of the basin. Analysis of the potential vorticity budget shows a basin-averaged balance primarily between potential vorticity input at the surface and dissipation at the lateral boundaries. However, advection is a leading-order term not only within the anticyclonic and cyclonic gyres but also between the gyres. This means that the eastern and western basins are dynamically connected through the advection of potential vorticity. Both eddy and mean fluxes play a role in connecting the regions of potential vorticity input at the surface with the opposite gyre and with the viscous boundary layers. These conclusions are based on a series of model runs in which forcing, topography, straits, and the Coriolis parameter were varied.

Description: This study was supported by National Science Foundation Grant OPP-1822334. Comments and suggestions from two anonymous referees greatly helped to improve the paper.

Description: 2021-02-17

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Blevins, M. G., Allen, H. L., Colson, B. C., Cook, A.-M., Greenbaum, A. Z., Hemami, S. S., Hollmann, J., Kim, E., LaRocca, A. A., Markoski, K. A., Miraglia, P., Mott, V. L., Robberson, W. M., Santos, J. A., Sprachman, M. M., Swierk, P., Tate, S., Witinski, M. F., Kratchman, L. B., & Michel, A. P. M. Field-portable microplastic sensing in aqueous environments: a perspective on emerging techniques. Sensors, 21(10), (2021): 3532, https://doi.org/10.3390/s21103532.

Description: Microplastics (MPs) have been found in aqueous environments ranging from rural ponds and lakes to the deep ocean. Despite the ubiquity of MPs, our ability to characterize MPs in the environment is limited by the lack of technologies for rapidly and accurately identifying and quantifying MPs. Although standards exist for MP sample collection and preparation, methods of MP analysis vary considerably and produce data with a broad range of data content and quality. The need for extensive analysis-specific sample preparation in current technology approaches has hindered the emergence of a single technique which can operate on aqueous samples in the field, rather than on dried laboratory preparations. In this perspective, we consider MP measurement technologies with a focus on both their eventual field-deployability and their respective data products (e.g., MP particle count, size, and/or polymer type). We present preliminary demonstrations of several prospective MP measurement techniques, with an eye towards developing a solution or solutions that can transition from the laboratory to the field. Specifically, experimental results are presented from multiple prototype systems that measure various physical properties of MPs: pyrolysis-differential mobility spectroscopy, short-wave infrared imaging, aqueous Nile Red labeling and counting, acoustophoresis, ultrasound, impedance spectroscopy, and dielectrophoresis.

Description: We greatly thank our funding agencies: Gerstner Philanthropies (to A.P.M.M.), the Richard Saltonstall Charitable Foundation (to A.P.M.M.), and the Wallace Research Foundation (to A.P.M.M. and S.S.H.). Funding for M.G.B. was provided by a Draper Fellowship and to B.C.C. by an MIT Martin Fellowship. Draper thanks EPA region 9 for their partnership and support through a Cooperative Research and Development Agreement, an industry/government agreement regarding funding and personnel contributions of time and expertise.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Meunier, T., Pallas Sanz, E., de Marez, C., Perez, J., Tenreiro, M., Ruiz Angulo, A., & Bower, A. The dynamical structure of a warm core ring as I\inferred from glider observations and along-track altimetry. Remote Sensing, 13(13), (2021): 2456, https://doi.org/10.3390/rs13132456.

Description: This study investigates the vertical structure of the dynamical properties of a warm-core ring in the Gulf of Mexico (Loop Current ring) using glider observations. We introduce a new method to correct the glider’s along-track coordinate, which is, in general, biased by the unsteady relative movements of the glider and the eddy, yielding large errors on horizontal derivatives. Here, we take advantage of the synopticity of satellite along-track altimetry to apply corrections on the glider’s position by matching in situ steric height with satellite-measured sea surface height. This relocation method allows recovering the eddy’s azimuthal symmetry, precisely estimating the rotation axis position, and computing reliable horizontal derivatives. It is shown to be particularly appropriate to compute the eddy’s cyclo-geostrophic velocity, relative vorticity, and shear strain, which are otherwise out of reach when using the glider’s raw traveled distance as a horizontal coordinate. The Ertel potential vorticity (PV) structure of the warm core ring is studied in details, and we show that the PV anomaly is entirely controlled by vortex stretching. Sign reversal of the PV gradient across the water column suggests that the ring might be baroclinically unstable. The PV gradient is also largely controlled by gradients of the vortex stretching term. We also show that the ring’s total energy partition is strongly skewed, with available potential energy being 3 times larger than kinetic energy. The possible impact of this energy partition on the Loop Current rings longevity is also discussed.

Description: This research was funded by a grant of the National Council of Science and Technology of Mexico—Secretariat of Energy Hydrocarbons Trust, project 201441. This is a contribution of the Gulf of Mexico Research Consortium (CIGoM).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1), (2021): 19-35, https://doi.org/10.1175/JPO-D-19-0233.1.

Description: In the Beaufort Sea in September of 2015, concurrent mooring and microstructure observations were used to assess dissipation rates in the vicinity of 72°35′N, 145°1′W. Microstructure measurements from a free-falling profiler survey showed very low [O(10−10) W kg−1] turbulent kinetic energy dissipation rates ε. A finescale parameterization based on both shear and strain measurements was applied to estimate the ratio of shear to strain Rω and ε at the mooring location, and a strain-based parameterization was applied to the microstructure survey (which occurred approximately 100 km away from the mooring site) for direct comparison with microstructure results. The finescale parameterization worked well, with discrepancies ranging from a factor of 1–2.5 depending on depth. The largest discrepancies occurred at depths with high shear. Mean Rω was 17, and Rω showed high variability with values ranging from 3 to 50 over 8 days. Observed ε was slightly elevated (factor of 2–3 compared with a later survey of 11 profiles taken over 3 h) from 25 to 125 m following a wind event which occurred at the beginning of the mooring deployment, reaching a maximum of ε= 6 × 10−10 W kg−1 at 30-m depth. Velocity signals associated with near-inertial waves (NIWs) were observed at depths greater than 200 m, where the Atlantic Water mass represents a reservoir of oceanic heat. However, no evidence of elevated ε or heat fluxes was observed in association with NIWs at these depths in either the microstructure survey or the finescale parameterization estimates.

Description: This work was supported by NSF Grants PLR 14-56705 and PLR-1303791 and by NSF Graduate Research Fellowship Grant DGE-1650112.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1),(2021): 3-17, https://doi.org/10.1175/JPO-D-20-0064.1.

Description: The strong El Niño of 2014–16 was observed west of the Galápagos Islands through sustained deployment of underwater gliders. Three years of observations began in October 2013 and ended in October 2016, with observations at longitudes 93° and 95°W between latitudes 2°N and 2°S. In total, there were over 3000 glider-days of data, covering over 50 000 km with over 12 000 profiles. Coverage was superior closer to the Galápagos on 93°W, where gliders were equipped with sensors to measure velocity as well as temperature, salinity, and pressure. The repeated glider transects are analyzed to produce highly resolved mean sections and maps of observed variables as functions of time, latitude, and depth. The mean sections reveal the structure of the Equatorial Undercurrent (EUC), the South Equatorial Current, and the equatorial front. The mean fields are used to calculate potential vorticity Q and Richardson number Ri. Gradients in the mean are strong enough to make the sign of Q opposite to that of planetary vorticity and to have Ri near unity, suggestive of mixing. Temporal variability is dominated by the 2014–16 El Niño, with the arrival of depressed isopycnals documented in 2014 and 2015. Increases in eastward velocity advect anomalously salty water and are uncorrelated with warm temperatures and deep isopycnals. Thus, vertical advection is important to changes in heat, and horizontal advection is relevant to changes in salt. Implications of this work include possibilities for future research, model assessment and improvement, and sustained observations across the equatorial Pacific.

Description: We gratefully acknowledge the support of the National Science Foundation (OCE-1232971, OCE-1233282) and the Ocean Observing and Monitoring Division of the National Oceanographic and Atmospheric Administration (NA13OAR4830216).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 38(1), (2021): 3-16, https://doi.org/10.1175/JTECH-D-20-0110.1.

Description: Airborne expendable bathythermographs (AXBTs) are air-launched, single-use temperature–depth probes that telemeter temperature observations as VHF-modulated frequencies. This study describes the AXBT Real-Time Editing System (ARES), which is composed of two components: the ARES Data Acquisition System, which receives telemetered temperature–depth profiles with no external hardware other than a VHF radio receiver, and the ARES Profile Editing System, which quality controls AXBT temperature–depth profiles. The ARES Data Acquisition System performs fast Fourier transforms on windowed segments of the demodulated signal transmitted from the AXBT. For each segment, temperature is determined from peak frequency and depth from elapsed time since profile start. Valid signals are distinguished from noise by comparing peak signal levels and signal-to-noise ratios to predetermined thresholds. When evaluated using 387 profiles, the ARES Data Acquisition System produced temperature–depth profiles nearly identical to those generated using a Sippican MK-21 processor, while reducing the amount of noise from VHF interference included in those profiles. The ARES Profile Editor applies a series of automated checks to identify and correct common profile discrepancies before displaying the profile on an editing interface that provides simple user controls to make additional corrections. When evaluated against 1177 tropical Atlantic and Pacific AXBT profiles, the ARES automated quality control system successfully corrected 87% of the profiles without any required manual intervention. Necessary future work includes improvements to the automated quality control algorithm and algorithm evaluation against a broader dataset of temperature–depth profiles from around the world across all seasons.

Description: This work was sponsored by the Office of Naval Research (Grants N000141812819 and N0001420WX00345) and the U.S. Navy’s Civilian Institution Office with the MIT–WHOI Joint Program.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(11), (2020): 3267–3294, https://doi.org/10.1175/JPO-D-19-0310.1.

Description: As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure p and round-trip acoustic travel time τ records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.

Description: Support for this research was provided by Office of Naval Research Grants N00014-16-1-2668, N00014-18-1-2406, N00014-15-1-2488, and N00014-15-1-2622. R.C.M. was additionally supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(11),(2020): 3331–3351, https://doi.org/10.1175/JPO-D-20-0035.1.

Description: This study examines the generation of warm spiral structures (referred to as spiral streamers here) over Gulf Stream warm-core rings. Satellite sea surface temperature imagery shows spiral streamers forming after warmer water from the Gulf Stream or newly formed warm-core rings impinges onto old warm-core rings and then intrudes into the old rings. Field measurements in April 2018 capture the vertical structure of a warm spiral streamer as a shallow lens of low-density water winding over an old ring. Observations also show subduction on both sides of the spiral streamer, which carries surface waters downward. Idealized numerical model simulations initialized with observed water-mass densities reproduce spiral streamers over warm-core rings and reveal that their formation is a nonlinear submesoscale process forced by mesoscale dynamics. The negative density anomaly of the intruding water causes a density front at the interface between the intruding water and surface ring water, which, through thermal wind balance, drives a local anticyclonic flow. The pressure gradient and momentum advection of the local interfacial flow push the intruding water toward the ring center. The large-scale anticyclonic flow of the ring and the radial motion of the intruding water together form the spiral streamer. The observed subduction on both sides of the spiral streamer is part of the secondary cross-streamer circulation resulting from frontogenesis on the stretching streamer edges. The surface divergence of the secondary circulation pushes the side edges of the streamer away from each other, widens the warm spiral on the surface, and thus enhances its surface signal.

Description: Authors W. G. Zhang and D. J. McGillicuddy are both supported by the National Science Foundation through Grant OCE 1657803.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1), (2021): 247-266, https://doi.org/10.1175/JPO-D-20-0098.1.

Description: This study focuses on mechanisms of shelf valley bathymetry affecting the spread of riverine freshwater in the nearshore region. In the context of Changjiang River, a numerical model is used with different no-tide idealized configurations to simulate development of unforced river plumes over a sloping bottom, with and without a shelf valley off the estuary mouth. All simulated freshwater plumes are surface-trapped with continuously growing bulges near the estuary mouth and narrow coastal currents downstream. The simulations indicate that a shelf valley tends to compress the bulge along the direction of the valley long axis and modify the incident angle of the bulge flow impinging toward the coast, which then affects the strength of the coastal current. The bulge compression results from geostrophic adjustment and isobath-following tendency of the depth-averaged flow in the bulge region. Generally, the resulting change in the direction of the bulge impinging flow enhances down-shelf momentum advection and freshwater delivery into the coastal current. Sensitivity simulations with altered river discharges Q, Coriolis parameter, shelf bottom slope, valley geometry, and ambient stratification show that enhancement of down-shelf freshwater transport in the coastal current, ΔQc, increases with increasing valley depth within the bulge region and decreasing slope Burger number of the ambient shelf. Assuming potential vorticity conservation, a scaling formula of ΔQc/Q is developed, and it agrees well with results of the sensitivity simulations. Mechanisms of valley influences on unforced river plumes revealed here will help future studies of topographic influence on river plumes under more realistic conditions.

Description: This work is conducted by Canbo Xiao and Weifeng (Gordon) Zhang during CX’s one-year visit at Woods Hole Oceanographic Institution (WHOI) in 2018–19. CX was supported by China Scholarship Council.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Spingys, C. P., Garabato, A. C. N., Legg, S., Polzin, K. L., Abrahamsen, E. P., Buckingham, C. E., Forryan, A., & Frajka-Williams, E. E. Mixing and transformation in a deep western boundary current: a case study. Journal of Physical Oceanography, 51(4), (2021): 1205-1222, https://doi.org/10.1175/JPO-D-20-0132.1

Description: Water-mass transformation by turbulent mixing is a key part of the deep-ocean overturning, as it drives the upwelling of dense waters formed at high latitudes. Here, we quantify this transformation and its underpinning processes in a small Southern Ocean basin: the Orkney Deep. Observations reveal a focusing of the transport in density space as a deep western boundary current (DWBC) flows through the region, associated with lightening and densification of the current’s denser and lighter layers, respectively. These transformations are driven by vigorous turbulent mixing. Comparing this transformation with measurements of the rate of turbulent kinetic energy dissipation indicates that, within the DWBC, turbulence operates with a high mixing efficiency, characterized by a dissipation ratio of 0.6 to 1 that exceeds the common value of 0.2. This result is corroborated by estimates of the dissipation ratio from microstructure observations. The causes of the transformation are unraveled through a decomposition into contributions dependent on the gradients in density space of the: dianeutral mixing rate, isoneutral area, and stratification. The transformation is found to be primarily driven by strong turbulence acting on an abrupt transition from the weakly stratified bottom boundary layer to well-stratified off-boundary waters. The reduced boundary layer stratification is generated by a downslope Ekman flow associated with the DWBC’s flow along sloping topography, and is further regulated by submesoscale instabilities acting to restratify near-boundary waters. Our results provide observational evidence endorsing the importance of near-boundary mixing processes to deep-ocean overturning, and highlight the role of DWBCs as hot spots of dianeutral upwelling.

Description: CS, ACNG, AF, and EFW were supported by the U.K. Natural Environment Research Council (NERC) Grant NE/K013181/1. ACNG was supported by the Royal Society and Wolfson Foundation. EPA and CEB were supported by NERC Grant NE/K012843/1. CEB was funded by an MSCA grant (No. 798319) from the European Union’s Horizon 2020 program. EPA was supported by NERC Grant NE/N018095/1. SL and KP were supported by U.S. National Science Foundation Grants OCE-1536453 and OCE-1536779. SL acknowledges support of Award NA18OAR4320123 from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce. The statements, findings, conclusions, and recommendations are those of the authors, and do not necessarily reflect the views of the National Oceanic and Atmospheric Administration, or the U.S. Department of Commerce.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(7), (2021): 2087–2102, https://doi.org/10.1175/JPO-D-20-0255.1.

Description: The boundary current system in the Labrador Sea plays an integral role in modulating convection in the interior basin. Four years of mooring data from the eastern Labrador Sea reveal persistent mesoscale variability in the West Greenland boundary current. Between 2014 and 2018, 197 middepth intensified cyclones were identified that passed the array near the 2000-m isobath. In this study, we quantify these features and show that they are the downstream manifestation of Denmark Strait Overflow Water (DSOW) cyclones. A composite cyclone is constructed revealing an average radius of 9 km, maximum azimuthal speed of 24 cm s−1, and a core propagation velocity of 27 cm s−1. The core propagation velocity is significantly smaller than upstream near Denmark Strait, allowing them to trap more water. The cyclones transport a 200-m-thick lens of dense water at the bottom of the water column and increase the transport of DSOW in the West Greenland boundary current by 17% relative to the background flow. Only a portion of the features generated at Denmark Strait make it to the Labrador Sea, implying that the remainder are shed into the interior Irminger Sea, are retroflected at Cape Farewell, or dissipate. A synoptic shipboard survey east of Cape Farewell, conducted in summer 2020, captured two of these features that shed further light on their structure and timing. This is the first time DSOW cyclones have been observed in the Labrador Sea—a discovery that could have important implications for interior stratification.

Description: A. P. and R. S. P. were funded by National Science Foundation Grants OCE-1259618 and OCE-1756361. I. L. B. and F. S. were funded by National Science Foundation Grants OCE-1258823 and OCE-1756272. N. P. H. was supported by the Natural Environment Research Council U.K. OSNAP program (NE/K010875/1 and NE/K010700/1). M. A. S. was supported by NSF Grants OCE-1558742 and OPP-1822334.

Description: 2021-12-08

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(8),(2020): 2315-2321, doi:10.1175/JPO-D-19-0327.1.

Description: Low-frequency currents and eddies transport sediment, pathogens, larvae, and heat along the coast and between the shoreline and deeper water. Here, low-frequency currents (between 0.1 and 4.0 mHz) observed in shallow surfzone waters for 120 days during a wide range of wave conditions are compared with theories for generation by instabilities of alongshore currents, by ocean-wave-induced sea surface modulations, and by a nonlinear transfer of energy from breaking waves to low-frequency motions via a two-dimensional inverse energy cascade. For these data, the low-frequency currents are not strongly correlated with shear of the alongshore current, with the strength of the alongshore current, or with wave-group statistics. In contrast, on many occasions, the low-frequency currents are consistent with an inverse energy cascade from breaking waves. The energy of the low-frequency surfzone currents increases with the directional spread of the wave field, consistent with vorticity injection by short-crested breaking waves, and structure functions increase with spatial lags, consistent with a cascade of energy from few-meter-scale vortices to larger-scale motions. These results include the first field evidence for the inverse energy cascade in the surfzone and suggest that breaking waves and nonlinear energy transfers should be considered when estimating nearshore transport processes across and along the coast.

Description: Funding was provided by a Vannevar Bush Faculty Fellowship [from OUSD(R&E)] and NSF.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Cusack, J. M., Brearley, J. A., Garabato, A. C. N., Smeed, D. A., Polzin, K. L., Velzeboer, N., & Shakespeare, C. J. Observed eddy-internal wave interactions in the Southern Ocean. Journal of Physical Oceanography, 50(10), (2020): 3042-3062, doi:10.1175/JPO-D-20-0001.1.

Description: The physical mechanisms that remove energy from the Southern Ocean’s vigorous mesoscale eddy field are not well understood. One proposed mechanism is direct energy transfer to the internal wave field in the ocean interior, via eddy-induced straining and shearing of preexisting internal waves. The magnitude, vertical structure, and temporal variability of the rate of energy transfer between eddies and internal waves is quantified from a 14-month deployment of a mooring cluster in the Scotia Sea. Velocity and buoyancy observations are decomposed into wave and eddy components, and the energy transfer is estimated using the Reynolds-averaged energy equation. We find that eddies gain energy from the internal wave field at a rate of −2.2 ± 0.6 mW m−2, integrated from the bottom to 566 m below the surface. This result can be decomposed into a positive (eddy to wave) component, equal to 0.2 ± 0.1 mW m−2, driven by horizontal straining of internal waves, and a negative (wave to eddy) component, equal to −2.5 ± 0.6 mW m−2, driven by vertical shearing of the wave spectrum. Temporal variability of the transfer rate is much greater than the mean value. Close to topography, large energy transfers are associated with low-frequency buoyancy fluxes, the underpinning physics of which do not conform to linear wave dynamics and are thereby in need of further research. Our work suggests that eddy–internal wave interactions may play a significant role in the energy balance of the Southern Ocean mesoscale eddy and internal wave fields.

Description: Funding for DIMES was provided by U.K. Natural Environment Research Council (NERC) Grants NE/E007058/1 and NE/E005667/1. JMC acknowledges the support of a NERC PhD studentship, and ACNG that of the Royal Society and the Wolfson Foundation. NV acknowledges support from the ARC Centre of Excellence for Climate Extremes (CLEX) Honours Scholarship and the ANU PBSA Partnership - Spotless Scholarship. CJS acknowledges support from an ARC Discovery Early Career Researcher Award DE180100087 and an Australian National University Futures Scheme award. Numerical simulations were conducted on the National Computational Infrastructure (NCI) facility, Canberra, Australia. This study has been conducted using E.U. Copernicus Marine Service Information. We thank two anonymous reviewers for their comments which helped to improve the manuscript significantly. Codes and output files are available online at the project repository (https://github.com/jessecusack/DIMES_eddy_wave_interactions).

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(9),(2020): 2797-2814, https://doi.org/10.1175/JPO-D-19-0326.1.

Description: Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10−7 W kg−1 and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10−3 to 10.1 × 10−3 m2 s−1. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 106 m3 s−1) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10−3 ± 5 × 10−4 m2 s−1, with an associated downward turbulent heat flux of 9 W m−2.

Description: This study is supported by NSFC (91858204), the CAS Strategic Priority Research Program (XDB42000000), NSFC(41720104008, 41421005, 41876025), QMSNL (2018SDKJ0104-02), and the Shandong Provincial projects (U1606402). L. Pratt was supported by the U.S. NSF Grant OCE-1657870.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(12), (2020): 3379-3402, https://do.org/10.1175/JPO-D-20-0028.1.

Description: One of the largest and most persistent features in the Alboran Sea is the Western Alboran Gyre (WAG), an anticyclonic recirculation bounded by the Atlantic Jet (AJ) to the north and the Moroccan coast to the south. Eulerian budgets from several months of a high-resolution model run are used to examine the exchange of water across the Eulerian WAG’s boundary and the processes affecting the salinity, temperature, and vorticity of the WAG. The volume transport across the sides of the WAG is found to be related to vertical isopycnal movement at the base of the gyre. Advection is found to drive a decay in the salinity minimum and anticyclonic vorticity of the Eulerian WAG. Given the large contributions of advection, a Lagrangian analysis is performed, revealing geometric aspects of the exchange that are hidden in an Eulerian view. In particular, stable and unstable manifolds identify a stirring region around the outer reaches of the gyre where water is exchanged with the WAG on a time scale of weeks. Its complement is an inner core that expands with depth and exchanges water with its surroundings on much longer time scales. The 3D evolution of one parcel, or lobe, of water as it enters the WAG is also described, identifying a general Lagrangian subduction pathway.

Description: This work was supported on DOD (MURI) Grant N000141110087, ONR Grant N000141812417, and NSF Grant OCE-1558806.

Description: 2021-05-18

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(7), (2020): 1839-1852, https://doi.org/10.1175/JPO-D-19-0209.1.

Description: The Lagrangian characteristics of the surface flow field arising when an idealized, anticyclonic, mesoscale, isolated deep-ocean eddy collides with continental slope and shelf topography are explored. In addition to fluid parcel trajectories, we consider the trajectories of biological organisms that are able to navigate and swim, and for which shallow water is a destination. Of particular interest is the movement of organisms initially located in the offshore eddy, the manner in which the eddy influences the ability of the organisms to reach the shelf break, and the spatial and temporal distributions of organisms that do so. For nonswimmers or very slow swimmers, the organisms arrive at the shelf break in distinct pulses, with different pulses occurring at different locations along the shelf break. This phenomenon is closely related to the episodic formation of trailing vortices that are formed after the eddy collides with the continental slope, turns, and travels parallel to the coast. Analysis based on finite-time Lyapunov exponents reveals initial locations of all successful trajectories reaching the shoreline, and provides maps of the transport pathways showing that much of the cross-shelf-break transport occurs in the lee of the eddy as it moves parallel to the shore. The same analysis shows that the onshore transport is interrupted after a trailing vortex detaches. As the swimming speeds are increased, the organisms are influenced less by the eddy and tend to show up en mass and in a single pulse.

Description: IR and LP were supported by National Science Foundation (NSF) Grant OCE-1558806. DC was supported by NSF U.S. National Science Foundation’s Physical Oceanography program through Grants OCE-1059632 and OCE-1433953 as well as the Academic Programs Office, Woods Hole Oceanographic Institution. We acknowledge high-performance computing support from Yellowstone (http://n2t.net/ark:/85065/d7wd3xhc) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 34(5), (2021): 1767-1788, https://doi.org/10.1175/JCLI-D-19-1020.1.

Description: Marine heatwaves along the coast of Western Australia, referred to as Ningaloo Niño, have had dramatic impacts on the ecosystem in the recent decade. A number of local and remote forcing mechanisms have been put forward; however, little is known about the depth structure of such temperature extremes. Utilizing an eddy-active global ocean general circulation model, Ningaloo Niño and the corresponding cold Ningaloo Niña events are investigated between 1958 and 2016, with a focus on their depth structure. The relative roles of buoyancy and wind forcing are inferred from sensitivity experiments. Composites reveal a strong symmetry between cold and warm events in their vertical structure and associated large-scale spatial patterns. Temperature anomalies are largest at the surface, where buoyancy forcing is dominant, and extend down to 300-m depth (or deeper), with wind forcing being the main driver. Large-scale subsurface anomalies arise from a vertical modulation of the thermocline, extending from the western Pacific into the tropical eastern Indian Ocean. The strongest Ningaloo Niños in 2000 and 2011 are unprecedented compound events, where long-lasting high temperatures are accompanied by extreme freshening, which emerges in association with La Niñas, that is more common and persistent during the negative phase of the interdecadal Pacific oscillation. It is shown that Ningaloo Niños during La Niña phases have a distinctively deeper reach and are associated with a strengthening of the Leeuwin Current, while events during El Niño are limited to the surface layer temperatures, likely driven by local atmosphere–ocean feedbacks, without a clear imprint on salinity and velocity.

Description: The following support is gratefully acknowledged: the Feodor-Lynen Fellowship by the Alexander von Humboldt Foundation and the WHOI Postdoctoral Scholar program (to SR), the Office of Naval Research under project number N-00014-19-12646 (to GG), the James E. and Barbara V. Moltz Fellowship for Climate-Related Research (to CCU), and IndoArchipel from the Deutsche Forschungsgemeinschaft (DFG) as part of the Special Priority Program (SPP)-1889 “Regional Sea Level Change and Society” (SeaLevel) (for PW).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(8),(2021): 2425–2441, https://doi.org/10.1175/JPO-D-20-0317.1.

Description: The frequency and latitudinal dependence of the midlatitude wind-driven meridional overturning circulation (MOC) is studied using theory and linear and nonlinear applications of a quasigeostrophic numerical model. Wind forcing is varied either by changing the strength of the wind or by shifting the meridional location of the wind stress curl pattern. At forcing periods of less than the first-mode baroclinic Rossby wave basin crossing time scale, the linear response in the middepth and deep ocean is in phase and opposite to the Ekman transport. For forcing periods that are close to the Rossby wave basin crossing time scale, the upper and deep MOC are enhanced, and the middepth MOC becomes phase shifted, relative to the Ekman transport. At longer forcing periods the deep MOC weakens and the middepth MOC increases, but eventually for long enough forcing periods (decadal) the entire wind-driven MOC spins down. Nonlinearities and mesoscale eddies are found to be important in two ways. First, baroclinic instability causes the middepth MOC to weaken, lose correlation with the Ekman transport, and lose correlation with the MOC in the opposite gyre. Second, eddy thickness fluxes extend the MOC beyond the latitudes of direct wind forcing. These results are consistent with several recent studies describing the four-dimensional structure of the MOC in the North Atlantic Ocean.

Description: This study was supported by National Science Foundation Grant OCE-1947290.

Description: 2022-01-13

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(6), (2020): 1557-1582, doi:10.1175/JPO-D-19-0239.1.

Description: We examine various contributions to the vertical velocity field within large mesoscale eddies by analyzing multiple solutions to an idealized numerical model of a representative anticyclonic warm core Gulf Stream ring. Initial conditions are constructed to reproduce the observed density and nutrient profiles collected during the Warm Core Rings Program. The contributions to vertical fluxes diagnosed from the numerical simulations are compared against a divergence-based, semidiagnostic equation and a generalized omega equation to better understand the dynamics of the vertical velocity field. Frictional decay alone is found to be ineffective in raising isopycnals and transporting nutrients to the upper ocean. With representative wind forcing, the magnitude of vorticity gradient–induced Ekman pumping is not necessarily larger than the current-induced counterpart on a time scale relevant to ecosystem response. Under realistic forcing conditions, strain deformation can perturb the ring to be noncircular and induce vertical velocities much larger than the Ekman vertical velocities. Nutrient budget diagnosis, together with analysis of the relative magnitudes of the various types of vertical fluxes, allows us to describe the time-scale dependence of nutrient delivery. At time scales that are relevant to individual phytoplankton (from hours to days), the magnitudes of nutrient flux by Ekman velocities and deformation-induced velocities are comparable. Over the life span of a typical warm core ring, which can span multiple seasons, surface current–induced Ekman pumping is the most effective mechanism in upper-ocean nutrient enrichment because of its persistence in the center of anticyclones regardless of the direction of the wind forcing.

Description: This work was supported by the National Science Foundation Ocean Science Division under Grant OCE-1558960. PG also acknowledges support of the NASA Physical Oceanography Program Grant NNX16H59G. KC would like to thank D. McGillicuddy Jr. for inspiring discussions and suggestions during the course of this study. Constructive comments from two anonymous reviewers are appreciated.

Description: Recent volcanic gas compilations have urged the need to expand in-situ plume measurements to poorly studied, remote volcanic regions. Despite being recognized as one of the main volcanic epicenters on the planet, the Vanuatu arc remains poorly characterized for its subaerial emissions and their chemical imprints. Here, we report on the first plume chemistry data for Mount Garet, on the island of Gaua, one of the few persistent volatile emitters along the Vanuatu arc. Data were collected with a multi-component gas analyzer system (multi-GAS) during a field campaign in December 2018. The average volcanic gas chemistry is characterized by mean molar CO2/SO2, H2O/SO2, H2S/SO2 and H2/SO2 ratios of 0.87, 47.2, 0.13 and 0.01, respectively. Molar proportions in the gas plume are estimated at 95.9 11.6, 1.8 0.5, 2.0 0.01, 0.26 0.02 and 0.06 0.01, for H2O, CO2, SO2, H2S and H2. Using the satellite-based 10-year (2005–2015) averaged SO2 flux of ~434 t d􀀀1 for Mt. Garet, we estimate a total volatile output of about 6482 t d􀀀1 (CO2 ~259 t d􀀀1; H2O ~5758 t d􀀀1; H2S ~30 t d􀀀1; H2 ~0.5 t d􀀀1). This may be representative of a quiescent, yet persistent degassing period at Mt. Garet; whilst, as indicated by SO2 flux reports for the 2009–2010 unrest, emissions can be much higher during eruptive episodes. Our estimated emission rates and gas composition for Mount Garet provide insightful information on volcanic gas signatures in the northernmost part of the Vanuatu Arc Segment. The apparent CO2-poor signature of high-temperature plume degassing at Mount Garet raises questions on the nature of sediments being subducted in this region of the arc and the possible role of the slab as the source of subaerial CO2. In order to better address the dynamics of along-arc volatile recycling, more volcanic gas surveys are needed focusing on northern Vanuatu volcanoes.

Description: This research was conducted as part of the Trail by Fire II—Closing the Ring Project (PI: Y. Moussallam) funded by the National Geographic Society (grant number CP-122R-17), the Rolex Awards for Enterprise and the French national Research Institute for Development (IRD). J.L. also acknowledges travel funding support from Ministero dell’istruzione, dell’università e della ricerca (MIUR;) under grant n. PRIN2017-2017LMNLAW).

Description: Published

Description: id 7293

Description: 4V. Processi pre-eruttivi

Description: JCR Journal

Description: Between 28 March and 1 April 2020, Stromboli volcano erupted, with overflows from the NE crater rim spreading along the barren Sciara del Fuoco slope and reaching the sea along the NW coast of the island. Poor weather conditions did not allow a detailed observation of the crater zone through the cameras monitoring network, but a clear view of the lower slope and the flows expanding in the area allowed us to characterize the flow features. This evidence was integrated with satellite, GBInSAR, and seismic data, thus enabling a reconstruction of the whole volcanic event, which involved several small collapses of the summit cone and the generation of pyroclastic density currents (PDCs) spreading along the slope and on the sea surface. Satellite monitoring allowed for the mapping of the lava flow field and the quantification of the erupted volume, and GBInSAR continuous measurements detected the crater widening and the deflation of the summit cone caused by the last overflow. The characterization of the seismicity made it possible to identify the signals that are associated with the propagation of PDCs along the volcano flank and, for the first time, to recognize the signal that is produced by the impact of the PDCs on the coast.

Description: This work has been financially supported by the “Presidenza del Consiglio dei Ministri—Dipartimento della Protezione Civile” (Presidency of the Council of Ministers–Department of Civil Protection) (DPC-UNIFI Agreement 2019–2021; Scientific Responsibility: N.C.); this publication, however, does not necessarily reflect the position and the official policies of the Department. Additional funds for paper publication have been provided by INGV-OE.

Description: Published

Description: 3010

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Accurate tracking and forecasting of ash dispersal in the atmosphere and quantification of its uncertainty are of fundamental importance for volcanic risk mitigation. Numerical models and satellite sensors offer two complementary ways to monitor ash clouds in real time, but limits and uncertainties affect both techniques. Numerical forecasts of volcanic clouds can be improved by assimilating satellite observations of atmospheric ash mass load. In this paper, we present a data assimilation procedure aimed at improving the monitoring and forecasting of volcanic ash clouds produced by explosive eruptions. In particular, we applied the Local Ensemble Transform Kalman Filter (LETKF) to the results of the Volcanic Ash Transport and Dispersion model HYSPLIT. To properly simulate the release and atmospheric transport of volcanic ash particles, HYSPLIT has been initialized with the results of the eruptive column model PLUME-MoM. The assimilation procedure has been tested against SEVIRI measurements of the volcanic cloud produced during the explosive eruption occurred at Mt. Etna on 24 December 2018. The results show how the assimilation procedure significantly improves the representation of the current ash dispersal and its forecast. In addition, the numerical tests show that the use of the sequential Ensemble Kalman Filter does not require a precise initialization of the numerical model, being able to improve the forecasts as the assimilation cycles are performed.

Description: Published

Description: id 359

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Volcanic plume height is a key parameter in retrieving plume ascent and dispersal dynamics, as well as eruption intensity; all of which are crucial for assessing hazards to aircraft operations. One way to retrieve cloud height is the shadow technique. This uses shadows cast on the ground and the sun geometry to calculate cloud height. This technique has, however, not been frequently used, especially not with high-spatial resolution (30 m pixel) satellite data. On 26 October 2013, Mt Etna (Sicily, Italy) produced a lava fountain feeding an ash plume that drifted SW and through the approach routes to Catania international airport. We compared the proximal plume height time-series obtained from fixed monitoring cameras with data retrieved from a Landsat-8 Operational Land Imager image, with results being in good agreement. The application of the shadow technique to a single high-spatial resolution image allowed us to fully document the ascent and dispersion history of the plume–cloud system. We managed to do this over a distance of 60 km and a time period of 50 min, with a precision of a few seconds and vertical error on plume altitude of ±200 m. We converted height with distance to height with time using the plume dispersion velocity, defining a bent-over plume that settled to a neutral buoyancy level with distance. Potentially, the shadow technique defined here allows downwind plume height profiles and mass discharge rate time series to be built over distances of up to 260 km and periods of 24 h, depending on vent location in the image, wind speed, and direction.

Description: This research was funded by CNES-TOSCA (Terre Solide), grant number 10 3703 “Integration of sample return data and remote sensing for advanced understanding of volcanic ash formation and dispersion” (PI: Lucia Gurioli).

Description: Published

Description: id 3951

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: The Etna flank eruption that started on 24 December 2018 lasted a few days and involved the opening of an eruptive fissure, accompanied by a seismic swarm and shallow earthquakes, significant SO2 flux release, and by large and widespread ground deformation, especially on the eastern flank of the volcano. Lava fountains and ash plumes from the uppermost eruptive fissure accompanied the opening stage, causing disruption to Catania International Airport, and were followed by a quiet lava effusion within the barren Valle del Bove depression until 27 December. This was the first flank eruption to occur at Etna in the last decade, during which eruptive activity was confined to the summit craters and resulted in lava fountains and lava flow output from the crater rims. In this paper, we used ground and satellite remote sensing techniques to describe the sequence of events, quantify the erupted volumes of lava, gas, and tephra, and assess volcanic hazards.

Description: Published

Description: id 905

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Hahn, L. C., Storelvmo, T., Hofer, S., Parfitt, R., & Ummenhofer, C. C. Importance of Orography for Greenland cloud and melt response to atmospheric blocking. Journal of Climate, 33(10), (2020): 4187-4206, doi:10.1175/JCLI-D-19-0527.1.

Description: More frequent high pressure conditions associated with atmospheric blocking episodes over Greenland in recent decades have been suggested to enhance melt through large-scale subsidence and cloud dissipation, which allows more solar radiation to reach the ice sheet surface. Here we investigate mechanisms linking high pressure circulation anomalies to Greenland cloud changes and resulting cloud radiative effects, with a focus on the previously neglected role of topography. Using reanalysis and satellite data in addition to a regional climate model, we show that anticyclonic circulation anomalies over Greenland during recent extreme blocking summers produce cloud changes dependent on orographic lift and descent. The resulting increased cloud cover over northern Greenland promotes surface longwave warming, while reduced cloud cover in southern and marginal Greenland favors surface shortwave warming. Comparison with an idealized model simulation with flattened topography reveals that orographic effects were necessary to produce area-averaged decreasing cloud cover since the mid-1990s and the extreme melt observed in the summer of 2012. This demonstrates a key role for Greenland topography in mediating the cloud and melt response to large-scale circulation variability. These results suggest that future melt will depend on the pattern of circulation anomalies as well as the shape of the Greenland Ice Sheet.

Description: This research was supported by the Woods Hole Oceanographic Institution Summer Student Fellow program, by the U.S. National Science Foundation under AGS-1355339 to C.C.U., and by the European Research Council through Grant 758005.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(2), (2020): 415-437, doi:10.1175/JPO-D-19-0019.1.

Description: Results are presented from two dye release experiments conducted in the seasonal thermocline of the Sargasso Sea, one in a region of low horizontal strain rate (~10−6 s−1), the second in a region of intermediate horizontal strain rate (~10−5 s−1). Both experiments lasted ~6 days, covering spatial scales of 1–10 and 1–50 km for the low and intermediate strain rate regimes, respectively. Diapycnal diffusivities estimated from the two experiments were κz = (2–5) × 10−6 m2 s−1, while isopycnal diffusivities were κH = (0.2–3) m2 s−1, with the range in κH being less a reflection of site-to-site variability, and more due to uncertainties in the background strain rate acting on the patch combined with uncertain time dependence. The Site I (low strain) experiment exhibited minimal stretching, elongating to approximately 10 km over 6 days while maintaining a width of ~5 km, and with a notable vertical tilt in the meridional direction. By contrast, the Site II (intermediate strain) experiment exhibited significant stretching, elongating to more than 50 km in length and advecting more than 150 km while still maintaining a width of order 3–5 km. Early surveys from both experiments showed patchy distributions indicative of small-scale stirring at scales of order a few hundred meters. Later surveys show relatively smooth, coherent distributions with only occasional patchiness, suggestive of a diffusive rather than stirring process at the scales of the now larger patches. Together the two experiments provide important clues as to the rates and underlying processes driving diapycnal and isopycnal mixing at these scales.

Description: Results are presented from two dye release experiments conducted in the seasonal thermocline of the Sargasso Sea, one in a region of low horizontal strain rate (~10−6 s−1), the second in a region of intermediate horizontal strain rate (~10−5 s−1). Both experiments lasted ~6 days, covering spatial scales of 1–10 and 1–50 km for the low and intermediate strain rate regimes, respectively. Diapycnal diffusivities estimated from the two experiments were κz = (2–5) × 10−6 m2 s−1, while isopycnal diffusivities were κH = (0.2–3) m2 s−1, with the range in κH being less a reflection of site-to-site variability, and more due to uncertainties in the background strain rate acting on the patch combined with uncertain time dependence. The Site I (low strain) experiment exhibited minimal stretching, elongating to approximately 10 km over 6 days while maintaining a width of ~5 km, and with a notable vertical tilt in the meridional direction. By contrast, the Site II (intermediate strain) experiment exhibited significant stretching, elongating to more than 50 km in length and advecting more than 150 km while still maintaining a width of order 3–5 km. Early surveys from both experiments showed patchy distributions indicative of small-scale stirring at scales of order a few hundred meters. Later surveys show relatively smooth, coherent distributions with only occasional patchiness, suggestive of a diffusive rather than stirring process at the scales of the now larger patches. Together the two experiments provide important clues as to the rates and underlying processes driving diapycnal and isopycnal mixing at these scales.

Description: 2020-08-06

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 37(5), (2020): 789-806, doi:10.1175/JTECH-D-18-0244.1.

Description: Realistic ocean state prediction and its validation rely on the availability of high quality in situ observations. To detect data errors, adequate quality check procedures must be designed. This paper presents procedures that take advantage of the ever-growing observation databases that provide climatological knowledge of the ocean variability in the neighborhood of an observation location. Local validity intervals are used to estimate binarily whether the observed values are considered as good or erroneous. Whereas a classical approach estimates validity bounds from first- and second-order moments of the climatological parameter distribution, that is, mean and variance, this work proposes to infer them directly from minimum and maximum observed values. Such an approach avoids any assumption of the parameter distribution such as unimodality, symmetry around the mean, peakedness, or homogeneous distribution tail height relative to distribution peak. To reach adequate statistical robustness, an extensive manual quality control of the reference dataset is critical. Once the data have been quality checked, the local minima and maxima reference fields are derived and the method is compared with the classical mean/variance-based approach. Performance is assessed in terms of statistics of good and bad detections. It is shown that the present size of the reference datasets allows the parameter estimates to reach a satisfactory robustness level to always make the method more efficient than the classical one. As expected, insufficient robustness persists in areas with an especially low number of samples and high variability.

Description: This study has been conducted using EU Copernicus Marine Service Information and was supported by the European Union within the EU Copernicus Marine Service In Situ phase-I and phase-II contracts led by Ifremer. The publication was also supported by SOERE CTDO2 in France. The Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it (see http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System (http://doi.org/10.17882/42182). The marine mammal data were collected and made freely available by the International MEOP Consortium and the national programs that contribute to it (see http://www.meop.net; https://doi.org/10.17882/45461). Aleix Gelabert and Dídac Costa were the skippers of the OPOO, sponsored by the Intergovernmental Oceanographic Commission (UNESCO) and Pharmaton. The BWR is a periodic oceanic race organized by the Fundació Navegació Oceànica de Barcelona (FNOB). Reviewer D. Briand provided some useful comments on the final version of the draft paper before submission.

Description: 2020-11-04

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(5), (2020): 1245-1263, doi:10.1175/JPO-D-19-0213.1.

Description: We use laboratory experiments and theoretical modeling to investigate the surface expression of a subglacial discharge plume, as occurs at many fjords around Greenland. The experiments consider a fountain that is released vertically into a homogeneous fluid, adjacent either to a vertical or a sloping wall, that then spreads horizontally at the free surface before sinking back to the bottom. We present a model that separates the fountain into two separate regions: a vertical fountain and a horizontal, negatively buoyant jet. The model is compared to laboratory experiments that are conducted over a range of volume fluxes, density differences, and ambient fluid depths. It is shown that the nondimensionalized length, width, and aspect ratio of the surface expression are dependent on the Froude number, calculated at the start of the negatively buoyant jet. The model is applied to observations of the surface expression from a Greenland subglacial discharge plume. In the case where the discharge plume reaches the surface with negative buoyancy the model can be used to estimate the discharge properties at the base of the glacier.

Description: We gratefully acknowledge technical assistance from Anders Jensen and thank anonymous reviewers for improving the clarity of the manuscript. CM thanks the Weston Howard Jr. Scholarship for funding. Support to CC was given by NSF project OCE-1434041 and OCE-1658079.

Description: 2020-10-27

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(4), (2020): 887-905, doi:10.1175/JPO-D-19-0110.1.

Description: The Equatorial Undercurrent (EUC) encounters the Galápagos Archipelago on the equator as it flows eastward across the Pacific. The impact of the Galápagos Archipelago on the EUC in the eastern equatorial Pacific remains largely unknown. In this study, the path of the EUC as it reaches the Galápagos Archipelago is measured directly using high-resolution observations obtained by autonomous underwater gliders. Gliders were deployed along three lines that define a closed region with the Galápagos Archipelago as the eastern boundary and 93°W from 2°S to 2°N as the western boundary. Twelve transects were simultaneously occupied along the three lines during 52 days in April–May 2016. Analysis of individual glider transects and average sections along each line show that the EUC splits around the Galápagos Archipelago. Velocity normal to the transects is used to estimate net horizontal volume transport into the volume. Downward integration of the net horizontal transport profile provides an estimate of the time- and areal-averaged vertical velocity profile over the 52-day time period. Local maxima in vertical velocity occur at depths of 25 and 280 m with magnitudes of (1.7 ± 0.6) × 10−5 m s−1 and (8.0 ± 1.6) × 10−5 m s−1, respectively. Volume transport as a function of salinity indicates that water crossing 93°W south (north) of 0.4°S tends to flow around the south (north) side of the Galápagos Archipelago. Comparisons are made between previous observational and modeling studies with differences attributed to effects of the strong 2015/16 El Niño event, the annual cycle of local winds, and varying longitudes between studies of the equatorial Pacific.

Description: This work was supported by National Science Foundation (Grants OCE-1232971 and OCE-1233282) and the NASA Earth and Space Science Fellowship Program (Grant 80NSSC17K0443).

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(4), (2020): 1045-1064, doi:10.1175/JPO-D-19-0137.1.

Description: Three simulations of the circulation in the Gulf of Mexico (the “Gulf”) using different numerical general circulation models are compared with results of recent large-scale observational campaigns conducted throughout the deep (〉1500 m) Gulf. Analyses of these observations have provided new understanding of large-scale mean circulation features and variability throughout the deep Gulf. Important features include cyclonic flow along the continental slope, deep cyclonic circulation in the western Gulf, a counterrotating pair of cells under the Loop Current region, and a cyclonic cell to the south of this pair. These dominant circulation features are represented in each of the ocean model simulations, although with some obvious differences. A striking difference between all the models and the observations is that the simulated deep eddy kinetic energy under the Loop Current region is generally less than one-half of that computed from observations. A multidecadal integration of one of these numerical simulations is used to evaluate the uncertainty of estimates of velocity statistics in the deep Gulf computed from limited-length (4 years) observational or model records. This analysis shows that the main deep circulation features identified from the observational studies appear to be robust and are not substantially impacted by variability on time scales longer than the observational records. Differences in strengths and structures of the circulation features are identified, however, and quantified through standard error analysis of the statistical estimates using the model solutions.

Description: This work was supported by the Gulf Research Program of the National Academy of Sciences under Awards 2000006422 and 2000009966. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Gulf Research Program or the National Academy of Sciences. The authors acknowledge the GLORYS project for providing the ocean reanalysis data used in the ROMS simulation. GLORYS is jointly conducted by MERCATOR OCEAN, CORIOLIS, and CNRS/INSU. Installation, recovery, data acquisition, and processing of the CANEK group current-meter moorings were possible because of CICESE-PetróleosMexicanos Grant PEP-CICESE 428229851 and the dedicated work of the crew of the B/O Justo Sierra and scientists of the CANEK group. The authors thank Dr. Aljaz Maslo, CICESE, for assistance with analysis of model data. The Bureau of Ocean Energy Management (BOEM), U.S. Dept. of the Interior, provided funding for the Lagrangian Study of the Deep Circulation in the Gulf of Mexico and the Observations and Dynamics of the Loop Current study. HYCOM simulation data are available from the HYCOM data server (https://www.hycom.org/data/goml0pt04/expt-02pt2), MITgcm data are available from the ECCO data server (http://ecco.ucsd.edu/gom_results2.html), and the ROMS simulation data are available from GRIIDC (NA.x837.000:0001).

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(5),(2020): 1227-1244, doi:10.1175/JPO-D-19-0280.1.

Description: The Nordic seas are commonly described as a single basin to investigate their dynamics and sensitivity to environmental changes when using a theoretical framework. Here, we introduce a conceptual model for a two-basin marginal sea that better represents the Nordic seas geometry. In our conceptual model, the marginal sea is characterized by both a cyclonic boundary current and a front current as a result of different hydrographic properties east and west of the midocean ridge. The theory is compared to idealized model simulations and shows good agreement over a wide range of parameter settings, indicating that the physics in the two-basin marginal sea is well captured by the conceptual model. The balances between the atmospheric buoyancy forcing and the lateral eddy heat fluxes from the boundary current and the front current differ between the Lofoten and the Greenland Basins, since the Lofoten Basin is more strongly eddy dominated. Results show that this asymmetric sensitivity leads to opposing responses depending on the strength of the atmospheric buoyancy forcing. Additionally, the front current plays an essential role for the heat and volume budget of the two basins, by providing an additional pathway for heat toward the interior of both basins via lateral eddy heat fluxes. The variability of the temperature difference between east and west influences the strength of the different flow branches through the marginal sea and provides a dynamical explanation for the observed correlation between the front current and the slope current of the Norwegian Atlantic Current in the Nordic seas.

Description: We thank Ilker Fer and two anonymous reviewers whose comments improved this paper. S. L. Ypma and S. Georgiou were supported by NWO (Netherlands Organisation for Scientific Research) VIDI Grant 864.13.011 awarded to C. A. Katsman. M. A. Spall was supported by National Science Foundation Grants OCE-1558742 and OPP-1822334. E. Lambert is funded by the ERA4CS project INSeaPTION. The model data analyzed in this study are available on request from the corresponding author. This study has been conducted using E.U. Copernicus Marine Service Information. The altimeter products were produced by Ssalto/Duacs and distributed by Aviso+, with support from CNES (https://www.aviso.altimetry.fr).

Description: 2020-10-27

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of [publisher] for personal use, not for redistribution. The definitive version was published in Schlundt, M., Farrar, J. T., Bigorre, S. P., Plueddemann, A. J., & Weller, R. A. (2020). Accuracy of wind observations from open-ocean buoys: correction for flow distortion. Journal of Atmospheric and Oceanic Technology, 37(4), 687-703, doi:10.1175/JTECH-D-19-0132.1.

Description: The comparison of equivalent neutral winds obtained from (i) four WHOI buoys in the subtropics and (ii) scatterometer estimates at those locations reveals a root-mean-square (RMS) difference of 0.56–0.76 m s−1. To investigate this RMS difference, different buoy wind error sources were examined. These buoys are particularly well suited to examine two important sources of buoy wind errors because 1) redundant anemometers and a comparison with numerical flow simulations allow us to quantitatively assess flow distortion errors, and 2) 1-min sampling at the buoys allows us to examine the sensitivity of buoy temporal sampling/averaging in the buoy–scatterometer comparisons. The interanemometer difference varies as a function of wind direction relative to the buoy wind vane and is consistent with the effects of flow distortion expected based on numerical flow simulations. Comparison between the anemometers and scatterometer winds supports the interpretation that the interanemometer disagreement, which can be up to 5% of the wind speed, is due to flow distortion. These insights motivate an empirical correction to the individual anemometer records and subsequent comparison with scatterometer estimates show good agreement.

Description: We gratefully acknowledge the help of three anonymous reviewers, whose input greatly improved the paper. In particular, one reviewer pointed out a mistake in our initial interpretation of scatterometer stability, which was corrected in the final manuscript. JTF and MS were supported by NASA Grant NNX14AM71G (International Ocean Vector Winds Science Team). The SPURS observations were supported by NASA (Grants NNX11AE84G, NNX15AG20G, and 80NSSC18K1494). The Stratus, NTAS, and WHOTS ocean reference stations (ORS) are long-term surface moorings deployed as part of the OceanSITES (http://www.oceansites.org) component of the Global Ocean Observing System, and are supported by NOAA’s Climate Program Office’s Ocean Observing and Monitoring Division, as are RAW, AJP, and SPB through the Cooperative Institute for the North Atlantic Region (CINAR) under Cooperative Agreement NA14OAR4320158 with NOAA Climate Program Office (CPO) (FundRef No. 100007298). The technical staff of the UOP Group at WHOI and the crews of NOAA and UNOLS vessels have been essential to the successful long-term maintenance of the ORS.

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(11), (2019): 2867-2881, doi: 10.1175/JPO-D-19-0072.1.

Description: The Antarctic Circumpolar Current plays a central role in the ventilation of heat and carbon in the global ocean. In particular, the isopycnal slopes determine where each water mass outcrops and thus how the ocean interacts with the atmosphere. The region-integrated isopycnal slopes have been suggested to be eddy saturated, that is, stay relatively constant as the wind forcing changes, but whether or not the flow is saturated in realistic present day and future parameter regimes is unknown. This study analyzes an idealized two-layer quasigeostrophic channel model forced by a wind stress and a residual overturning generated by a mass flux across the interface between the two layers, with and without a blocking ridge. The sign and strength of the residual overturning set which way the isopycnal slopes change with the wind forcing, leading to an increase in slope with an increase in wind forcing for a positive overturning and a decrease in slope for a negative overturning, following the usual conventions; this behavior is caused by the dominant standing meander weakening as the wind stress weakens causing the isopycnal slopes to become more sensitive to changes in the wind stress and converge with the slopes of a flat-bottomed simulation. Eddy saturation only appears once the wind forcing passes a critical level. These results show that theories for saturation must have both topography and residual overturning in order to be complete and provide a framework for understanding how the isopycnal slopes in the Southern Ocean may change in response to future changes in wind forcing.

Description: MKY and RF acknowledge support through NSF Awards OCE-1536515 and AGS-1835576. MKY acknowledges funding from NDSEG. GRF was supported by NSF OCE-1459702. We are very grateful for conversations with David Marshall, Andrew Stewart, and two anonymous reviewers that greatly improved the manuscript. The code for running the model is found at https://github.com/mkyoungs/JPO-QG-Channel.

Description: 2020-04-30

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(12), (2019): 3061-3068, doi: 10.1175/JPO-D-18-0172.1.

Description: The calculation of energy flux in coastal trapped wave modes is reviewed in the context of tidal energy pathways near the coast. The significant barotropic pressures and currents associated with coastal trapped wave modes mean that large errors in estimating the wave flux are incurred if only the baroclinic component is considered. A specific example is given showing that baroclinic flux constitutes only 10% of the flux in a mode-1 wave for a reasonable choice of stratification and bathymetry. The interpretation of baroclinic energy flux and barotropic-to-baroclinic conversion at the coast is discussed: in contrast to the open ocean, estimates of baroclinic energy flux do not represent a wave energy flux; neither does conversion represent the scattering of energy from the tidal Kelvin wave to higher modes.

Description: This work was supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship, and by NSF under Grant OCE-1756781. I am grateful to K. Brink for the many useful conversations that contributed to this work and to J. Toole for providing detailed comments on an early version of this paper. The comments of three anonymous reviewers were very helpful in improving this paper.

Description: 2020-06-03

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 37(5), (2020): 807-824, doi:10.1175/JTECH-D-19-0054.1.

Description: Marine mammals are under growing pressure as anthropogenic use of the ocean increases. Ship strikes of large whales and loud underwater sound sources including air guns for marine geophysical prospecting and naval midfrequency sonar are criticized for their possible negative effects on marine mammals. Competent authorities regularly require the implementation of mitigation measures, including vessel speed reductions or shutdown of acoustic sources if marine mammals are sighted in sensitive areas or in predefined exclusion zones around a vessel. To ensure successful mitigation, reliable at-sea detection of animals is crucial. To date, ship-based marine mammal observers are the most commonly implemented detection method; however, thermal (IR) imaging–based automatic detection systems have been used in recent years. This study evaluates thermal imaging–based automatic whale detection technology for its use across different oceans. The performance of this technology is characterized with respect to environmental conditions, and an automatic detection algorithm for whale blows is presented. The technology can detect whales in polar, temperate, and subtropical ocean regimes over distances of up to several kilometers and outperforms marine mammal observers in the number of whales detected. These results show that thermal imaging technology can be used to assist in providing protection for marine mammals against ship strike and acoustic impact across the world’s oceans.

Description: This work was funded by the Office of Naval Research (ONR) under Award N000141310856, by the Environmental Studies Research Fund (ESRF; esrfunds.org) under Award 2014-03S and by the Alfred-Wegener-Institute Helmholtz Zentrum für Polar- und Meeresforschung. DPZ and OB declare competing financial interests: 1) Patent US8941728B2, DE102011114084B4: A method for automatic real-time marine mammal detection. The patent describes the ideas basic to the automatic whale detection software as used to acquire and process the data presented in this paper. 2) Licensing of the Tashtego automatic whale detection software to the manufacturer of IR sensor. The authors confirm that these competing financial interests did not alter their adherence good scientific practice. We thank P. Abgrall, J. Coffey, K. Keats, B. Mactavish, V. Moulton, and S. Penney-Belbin for data collection or IR image review. We thank S. Besaw, J. Christian, A. Coombs, P. Coombs, W. Costello, T. Elliott, E. Evans, I. Goudie, C. Jones, K. Knowles, R. Martin, A. Murphy, D. and J. Shepherd; and the staffs at the Irish Loop Express, the Myrick Wireless Interpretive Centre, the Mistaken Point Ecological Reserve, and the lighthouse keepers for logistical assistance at our remote field site. We thank D. Boutilier and B. McDonald (DFO) for assisting us in obtaining license to occupy permits for Cape Race. We thank D. Taylor (ESRF Research Manager) for his support.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(3), (2020): 715-726, doi:10.1175/JPO-D-19-0021.1.

Description: Closing the overturning circulation of bottom water requires abyssal transformation to lighter densities and upwelling. Where and how buoyancy is gained and water is transported upward remain topics of debate, not least because the available observations generally show downward-increasing turbulence levels in the abyss, apparently implying mean vertical turbulent buoyancy-flux divergence (densification). Here, we synthesize available observations indicating that bottom water is made less dense and upwelled in fracture zone valleys on the flanks of slow-spreading midocean ridges, which cover more than one-half of the seafloor area in some regions. The fracture zones are filled almost completely with water flowing up-valley and gaining buoyancy. Locally, valley water is transformed to lighter densities both in thin boundary layers that are in contact with the seafloor, where the buoyancy flux must vanish to match the no-flux boundary condition, and in thicker layers associated with downward-decreasing turbulence levels below interior maxima associated with hydraulic overflows and critical-layer interactions. Integrated across the valley, the turbulent buoyancy fluxes show maxima near the sidewall crests, consistent with net convergence below, with little sensitivity of this pattern to the vertical structure of the turbulence profiles, which implies that buoyancy flux convergence in the layers with downward-decreasing turbulence levels dominates over the divergence elsewhere, accounting for the net transformation to lighter densities in fracture zone valleys. We conclude that fracture zone topography likely exerts a controlling influence on the transformation and upwelling of bottom water in many areas of the global ocean.

Description: The data used in this study were collected in the context of several projects funded by the U.S. National Science Foundation (NSF), in particular BBTRE (OCE-9415589 and OCE-9415598) and DoMORE (OCE-1235094). Funding for the analysis was provided as part of the NSF DoMORE and DECIMAL (OCE-1735618) projects. Author Ijichi is a Japan Society for the Promotion of Science (JSPS) Overseas Research Fellow. Comments on an early draft of this paper by Jim Ledwell and Bryan Kaiser, as well as topical discussions with Jörn Callies and Trevor McDougall, are gratefully acknowledged. The paper was greatly improved during the review process, in particular because of the critical comments from one of the two anonymous reviewers.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(1), (2020): 255-268, doi:10.1175/JPO-D-19-0166.1.

Description: Regional connectivity is important to the global climate salinity response, particularly because salinity anomalies do not have a damping feedback with atmospheric freshwater fluxes and may therefore be advected over long distances by ocean circulation, resulting in nonlocal influences. Climate model intercomparison experiments such as CMIP5 exhibit large uncertainty in some aspects of the salinity response, hypothesized here to be a result of ocean dynamics. We use two types of Lagrangian particle tracking experiments to investigate pathways of exchange for salinity anomalies. The first uses forward trajectories to estimate average transport time scales between water cycle regimes. The second uses reverse trajectories and a freshwater accumulation method to quantitatively identify remote influences in the salinity response. Additionally, we compare velocity fields with both resolved and parameterized eddies to understand the impact of eddy stirring on intergyre exchange. These experiments show that surface anomalies are readily exchanged within the ocean gyres by the mean circulation, but intergyre exchange is slower and largely eddy driven. These dynamics are used to analyze the North Atlantic salinity response to climate warming and water cycle intensification, where the system is broadly forced with fresh surface anomalies in the subpolar gyre and salty surface anomalies in the subtropical gyres. Under these competing forcings, strong intergyre eddy fluxes carry anomalously salty subtropical water into the subpolar gyre which balances out much of the local freshwater input.

Description: We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. We also thank the creators of the SODA and ECCO reanalysis products. This work was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Award 80NSSC17K0372, and by National Science Foundation Award OCE-1433132. The SODA outputs used here can be accessed at http://www.atmos.umd.edu/~ocean/, and the ECCO outputs at https://ecco.jpl.nasa.gov/. Data from the CMIP5 ensemble is available at https://esgf-node.llnl.gov/projects/esgf-llnl/. The particle tracking code used for these experiments can be found at https://github.com/slevang/particle-tracking.

Description: 2020-07-20

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Martini, K. I., Murphy, D. J., Schmitt, R. W., & Larson, N. G. Reply to "comments on 'corrections for pumped SBE 41CP CTDs determined from stratified tank experiments'". Journal of Atmospheric and Oceanic Technology, 37(2), (2020): 357-363, doi:10.1175/JTECH-D-19-0171.1.

Description: The response in Johnson (2020) that the method used to determine cell thermal mass correction coefficients for SBE 41CP CTD data from Argo floats is biased as determined by Martini et al. (2019) is valid. However, the recommendation for correction coefficients should not be followed due to these three errors in Johnson (2020): Alignment is as large a source of dynamic error as cell thermal mass in the SBE 41CP CTD. Order of operations was overlooked, so that cell thermal mass is used to correct for alignment errors caused by the temporal mismatch of temperature and conductivity. The cell thermal mass corrections determined in Johnson et al. (2007) and Johnson (2020) also bias salinity. In this response we will do the following: Detail how the corrections in Johnson (2020) are biased because the optimization procedure does not accurately model physics in the tank and conductivity cell. Verify using in situ data from Argo floats deployed in the ocean that alignment is a significant source of error for the SBE 41CP as shown in Martini et al. (2019). Determine cell thermal mass correction coefficients from the stratified tank experiment merging the methods of Johnson (2020) and Martini et al. (2019) to optimize against a model that better represents the physics in the tank and conductivity cell. Compare the corrections using in situ data using the coefficients determined in Johnson et al. (2007), Martini et al. (2019), Johnson (2020), and this manuscript.

Description: Thanks to Pelle Robbins for finding the in situ profiles used for this analysis in the vast database of Argo floats, John Gilson showing me how to access that high-resolution data, Ray Schmitt for use of the stratified tank, Susan Wijffels, Breck Owens, and Annie Wong for intellectual support, and Diego Sorrentino and Vlad Simontov for validating the sampling scheme in the SBE 41CP.

Description: 2020-08-24

Description: The first in situ measurements of seawater density that referred to a geographical position at sea and time of the year were carried out by Count Luigi Ferdinando Marsili between 1679 and 1680 in the Adriatic Sea, Aegean Sea, Marmara Sea, and the Bosporus. Not only was this the first investigation with documented oceanographic measurements carried out at stations, but themeasurements were described in such an accurateway that the authorswere able to reconstruct the observations in modern units. These first measurements concern the ‘‘specific gravity’’ of seawaters (i.e., the ratio between fluid densities). The data reported in the historical oceanographic treatise Osservazioni intorno al Bosforo Tracio (Marsili) allowed the reconstruction of the seawater density at different geographic locations between 1679 and 1680. Marsili’s experimental methodology included the collection of surface and deep water samples, the analysis of the samples with a hydrostatic ampoule, and the use of a reference water to standardize the measurements.Acomparison of reconstructed densities with present-day values shows an agreement within 10%–20% uncertainty, owing to various aspects of the measurement methodology that are difficult to reconstruct from the documentary evidence. Marsili also measured the current speed and the depth of the current inversion in the Bosporus, which are consistent with the present-day knowledge. The experimental data collected in the Bosporus enabledMarsili to enunciate a theory on the cause of the two-layer flow at the strait, demonstrated by his laboratory experiment and later confirmed by many analytical and numerical studies.

Description: American Meteorological Society.

Description: Published

Description: 845 - 860

Description: 4A. Oceanografia e clima

Description: JCR Journal

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(7), (2019): 1889-1904, doi:10.1175/JPO-D-19-0053.1.

Description: A high-resolution numerical model, together with in situ and satellite observations, is used to explore the nature and dynamics of the dominant high-frequency (from one day to one week) variability in Denmark Strait. Mooring measurements in the center of the strait reveal that warm water “flooding events” occur, whereby the North Icelandic Irminger Current (NIIC) propagates offshore and advects subtropical-origin water northward through the deepest part of the sill. Two other types of mesoscale processes in Denmark Strait have been described previously in the literature, known as “boluses” and “pulses,” associated with a raising and lowering of the overflow water interface. Our measurements reveal that flooding events occur in conjunction with especially pronounced pulses. The model indicates that the NIIC hydrographic front is maintained by a balance between frontogenesis by the large-scale flow and frontolysis by baroclinic instability. Specifically, the temperature and salinity tendency equations demonstrate that the eddies act to relax the front, while the mean flow acts to sharpen it. Furthermore, the model reveals that the two dense water processes—boluses and pulses (and hence flooding events)—are dynamically related to each other and tied to the meandering of the hydrographic front in the strait. Our study thus provides a general framework for interpreting the short-time-scale variability of Denmark Strait Overflow Water entering the Irminger Sea.

Description: MAS was supported by the National Science Foundation (NSF) under Grants OCE-1558742 and OCE-1534618. RSP, PL, and DM were supported by NSF under Grants OCE-1558742 and OCE-1259618. WJvA was supported by the Helmholtz Infrastructure Initiative FRAM. TWNH and MA were supported by NSF under Grants OCE-1633124 and OCE-118123.

Description: 2020-07-01

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography, 49 (2), (2019): 607-630, doi:10.1175/JPO-D-18-0166.1.

Description: The Lagrangian motion in the eddy field produced from an unstable retrograde jet along the shelf break is studied from idealized numerical experiments with a primitive equation model. The jet is initially in thermal wind balance with a cross-isobath density gradient and is not subjected to any atmospheric forcing. Over the course of the model integration, the jet becomes unstable and produces a quasi-stationary eddy field over a 2-month period. During this period, the cross-slope flow at the shelf break is characterized by along-slope correlation scales of O(10) km and temporal correlation scales of a few days. The relative dispersion of parcels across isobaths is found to increase with time as tb, where 1 〈 b 〈 2. This mixed diffusive–ballistic regime appears to reflect the combined effects of (i) the short length scales of velocity correlation at the shelf break and (ii) the seaward excursion of monopolar and dipolar vortices. Cross-slope dispersion is greater offshore of the front than inshore of the front, as offshore parcels are both subducted onshore below density surfaces and translated offshore with eddies. Nonetheless, the exchange of parcels across the jet remains very limited on the monthly time scale. Particles originating from the bottom experience upward displacements of a few tens of meters and seaward displacements of O(100) km, suggesting that the eddy activity engendered by an unstable along-slope jet provides another mechanism for bottom boundary layer detachment near the shelf edge.

Description: The author expresses his gratitude to the researchers who contributed to the development and public dissemination of POM [for a list of contributors, see Mellor (2002) and comments in the source code]. Discussions with Kenneth Brink, Hyodae Seo, and Weifeng Zhang have been helpful. Comments provided by Kenneth Brink on a draft are gratefully acknowledged. The criticism from two anonymous reviewers allowed us to better focus the manuscript and to significantly improve its clarity. This work has been supported by Grant OCE-1556400 from the U.S. National Science Foundation.

Description: 2020-02-18