Journal of Geophysical Research: Solid Earth最新文献

Lahars are among the most frequent hazards associated with Volcán de Fuego, Guatemala. Despite their recurrence, early detection and automated alerts remain challenging since they often rely on manual monitoring and sparse visual confirmation. Yet, we can harness the high number of flows triggered every rainy season to characterize their seismic signatures and quantify their size and behavior. For this, we used seismic stations located along two active lahar channels on Fuego where this characterization describes a somewhat stable long-term flow behavior. This work revealed more varied short-term behavior characterized by increasing seismic activity in the time domain and a shift toward lower frequencies as these flows propagate downstream. Building on this characterization, we implemented K-nearest neighbor (KNN) based detectors using seismic signal attributes describing samples of the data in the time and frequency domains, as well as statistical functions of these samples. We trained generalized and station-specific detectors that achieved high accuracy for detecting moderate-to-large flows with lower performance for smaller or ambiguous events. We found that root mean square amplitude, a proxy for flow size, appears to control detector performance more than other signal features. The detector is computationally efficient and, in the case of Fuego, did not require additional instrumentation. This framework presents a portable solution for enhancing automated lahar detection while minimizing the use of location-specific parameters required by other methods.

Volcano deformation measured through Interferometric Synthetic Aperture Radar (InSAR) is ideal for volcano monitoring in many regions due to its global coverage, characteristic spatio-temporal patterns, and modeling insights. Routinely acquired and processed Sentinel-1 InSAR datacubes provide the first opportunity to systematically catalog, model and compare volcano deformation globally. Here, we present a framework (GBIS-BULK) to systematically pre-process and model volcano deformation signals, designed to be applied to routinely processed InSAR data sets. This requires a robust (semi-) automated approach to estimate signal locations and footprints for effective pre-processing and modeling. Our approach combines (a) filtering and clustering to locate the signal center; (b) noise reduction using Independent Component Analysis (ICA); and (c) image classification using Otsu thresholding to delimit the signal footprint. We invert for the best-fit point source model using constraints from existing global volcano deformation catalogs. First, we examine the influence of downsampling schemes, image noise and coherence using synthetic interferograms, showing nested-uniform downsampling is more suited to automated processing than quadtree methods which typically require manual tuning. Then, we validate the approach using Sentinel-1 deformation images from the East African Rift System (EARS). The pre-processing steps reasonably locate the signal at 15/16 of the EARS volcanoes, and the signal footprint at 14/16. ICA reduces or approximately maintains the image RMS in all cases. Our systematic point source estimates showed consistency when directly compared with previous (bespoke) modeling studies. This approach has the potential to be integrated with existing toolkits for routinely processing and analyzing Sentinel-1 InSAR data and hence applied globally.

The geology of New Zealand has been shaped by tectonic plate interactions driven by mantle convection over the past 60 million years, but the effects of these interactions on the transition to the lower mantle are not yet well understood. We analyze 10 years of teleseismic $��P�$-wave receiver functions using common conversion point stacking to investigate mantle transition zone discontinuities. The resulting topography of the 410-km and 660-km discontinuities, smoothed using kriging interpolation, reveals localized thermal anomalies. We observe more than 10 km of thinning in the mantle transition zone beneath central North Island coincident with the Hikurangi slab, accompanied by a significantly thickened and hydrous melt-rich layer above the 410-km discontinuity. In northern South Island, the 410-km discontinuity is uplifted by approximately 15 km, likely reflecting subducted materials from the oblique Hikurangi slab material reaching mantle transition zone depths. These observations are consistent with patterns of seismicity and seismic velocity anomalies from tomographic models. Additionally, we identify significant thinning of the mantle transition zone beneath both the Northland region ($��\sim �$10 km) and the Great South Basin ($��\sim �$15 km), indicating the influence of potential local thermal upwellings. The different slab polarities and transpressional boundary may drive lateral mantle flow, shaping the complex transition zone topography beneath New Zealand.

Fluids released from subducting hydrated rocks influence volcanism, tectonics, and geochemical cycling, but the mechanisms of fluid escape in subduction zones remain poorly understood. We address this issue by investigating the Erro-Tobbio meta-serpentinites (ET-MS), Italy, exhumed serpentinite rocks that preserve extensive dehydration vein networks formed by the porosity-generating breakdown of antigorite and brucite. We characterized the structure and morphology of these self-organized vein networks and evaluated their hydrodynamic properties using a novel approach. Specifically, we combined X-ray tomography and drone imagery with generative machine learning, electron microscopy, and equilibrium thermodynamics to model and analyze fluid pathways in the ET-MS. In both natural and simulated samples, these dehydration vein networks act as efficient drainage systems, enabling rapid fluid percolation even at porosities below 1%. The maximum network permeability is $��\approx �1�0^{�-�15�}�{��m�}^2�$, several orders of magnitude higher than that of intact serpentinite. Fe-rich olivine and monticellite occur alongside relict brucite and magnetite in these veins. This assemblage indicates that the high permeability arises from porosity localized along brucite- and magnetite-rich veins, where infiltration of reducing fluids enhanced dehydration reactions. These findings demonstrate that serpentinite dehydration in subduction zones can produce flow-optimized vein structures that efficiently channel fluids at low porosity, potentially influencing fluid migration on local to regional scales before widespread dehydration occurs.

The Greater Alpine crust provides a natural laboratory for investigating tectonic and geodynamic processes owing to its strong structural and lithological heterogeneity. A key challenge is constraining its thermal and compositional properties, given the limited direct observations and the uncertainties of existing models. We analyze two years of seismic ambient noise recorded at about 700 broadband seismic stations, and eight years of teleseismic earthquakes recorded at approximately 400 broadband stations. Using ambient noise tomography and receiver functions, we derive the Vs of the crust, the average Vp/Vs ratio, and a new Moho depth map for the Greater Alpine region. Moho depths range from 15 to 55 km, and Vp/Vs ratios vary from <1.65 in the Variscan domain to >1.8 in the Apennines, Dinarides and sedimentary basins. Thermodynamic modeling translates these seismic results into quantitative estimates of silica content and linear thermal gradients. We find that the Variscan crust is relatively cold (<20°C/km) and silica-rich (>62 wt% $��S�i�O_2�$), in contrast with Alpine–Apenninic domains that show elevated thermal gradients (>25°C/km) and more mafic compositions (<60 wt% $��S�i�O_2�$). However, elastic velocities predicted by mineral physics systematically exceed observed seismic velocities, reflecting effects of sediments, porosity, anelastic relaxation, and non-equilibrium processes not captured in thermodynamic equilibrium models. We apply empirical corrections to address these discrepancies, emphasizing caution when interpreting seismic velocities from mineral physics models. These quantitative constraints refine the thermal and chemical distinction between stable Paleozoic lithosphere and actively deforming orogens in the Greater Alpine crust.

The study of silicate glasses is important to understand the physical and chemical properties of silicate melts under high-pressure conditions relevant to planetary interiors. We conducted in situ Brillouin spectroscopy measurements on two endmember, low-impurity CaSiO₃ glasses and one Fe, Al, Mg, Ti-bearing CaSiO₃ glass up to 23 GPa. We obtained pressure-dependent acoustic velocities and derived elastic moduli that exhibit discontinuities indicative of structural transitions in all compositions. The endmember CaSiO₃ glasses exhibit velocity softening below 2.6 GPa, consistent with earlier findings in other silicate glasses, while the Fe, Al, Mg, Ti-bearing CaSiO₃ glass displays a delayed onset of softening and densification. This softening is attributed to intermediate-range structural rearrangements. Above ∼8 GPa, both glasses show rapid increases in velocities and elastic moduli, reflecting densification associated with structural transitions. Comparison with other silicate glasses demonstrates that Ca acts as a strong network modifier, significantly reducing the stiffness of the glasses, while the addition of Fe, Al, Mg, and Ti collectively has a mixed effect on elasticity. CaSiO₃ glass crosses over in density with its counterpart crystal, wollastonite, at only ∼3 GPa—a pressure lower than any other crossover pressure observed in silicate glasses—suggesting that Ca-rich melts may become gravitationally stable at much shallower depths in planetary interiors than other silicate melts. These results provide new constraints on the structural evolution and elasticity of Ca-rich silicate glasses under compression and have implications for modeling the mobility of silicate melts in deep planetary environments.