Journal of Geophysical Research. Solid Earth

A series of discrete elements simulations is presented for the study of fault gouges' frictional response. The gouge is considered to have previously undergone ultra-cataclastic flow and long-time consolidation loading. We explore the effect of different particle characteristics such as size, polydispersity, and also shearing velocities on gouge's response under the conditions met in the seismogenic zone. Monte-Carlo analyses suggest that the local stick-slip events disappear when averaging over a large number of numerical samples...

Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high-frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring...

Our understanding of geomagnetic field intensity prior to the era of direct instrumental measurements relies on paleointensity analysis of rocks and archaeological materials that serve as magnetic recorders. Only in rare cases are absolute paleointensity data sets continuous over millennial timescales, in sub-centennial resolution, and directly dated using radiocarbon. As a result, fundamental properties of the geomagnetic field, such as its maximum intensity and rate of change have remained a subject of lively discussion...

Accurate earthquake location and magnitude estimation play critical roles in seismology. Recent deep learning frameworks have produced encouraging results on various seismological tasks (e.g., earthquake detection, phase picking, seismic classification, and earthquake early warning). Many existing machine learning earthquake location methods utilize waveform information from a single station. However, multiple stations contain more complete information for earthquake source characterization. Inspired by recent successes in applying graph neural networks (GNNs) in graph-structured data, we develop a Spatiotemporal Graph Neural Network (STGNN) for estimating earthquake locations and magnitudes...

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016-2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution...

FeO represents an important end-member for planetary interiors mineralogy. However, its properties in the liquid state under high pressure are poorly constrained. Here, in situ high-pressure and high-temperature X-ray diffraction experiments, ab initio simulations, and thermodynamic calculations are combined to study the local structure and density evolution of liquid FeO under extreme conditions. Our results highlight a strong shortening of the Fe-Fe distance, particularly pronounced between ambient pressure and ∼40 GPa, possibly related with the insulator to metal transition occurring in solid FeO over a similar pressure range...

Understanding the stress evolution of extinct volcanoes can improve efforts to forecast flank eruptions on active systems. Field, petrographic, and seismic data are combined with numerical modeling to investigate the paleo-stress field of New Zealand's Akaroa Volcano, or Akaroa Volcanic Complex. Field mapping identifies 86 radially oriented dikes and seven lava domes found only within a narrow elevation range along Akaroa's erosional crater rim. These observations suggest that crater rim dike emplacement resulted from lateral deflection of vertically ascending intrusions from a centralized magma source, which in turn may have facilitated formation of the lava domes, as well as two scoria cones...

Fault surfaces are rough at all scales, and this significantly affects fault-slip behavior. However, roughness is only occasionally considered experimentally and then often in experiments imposing a low-slip velocity, corresponding to the initiation stage of the earthquake cycle. Here, the effect of roughness on earthquake nucleation up to runaway slip is investigated through a series of dry load-stepping biaxial experiments performed on bare rock surfaces with a variety of roughnesses. These laboratory faults reached slip velocities of at least 100 mm/s...

Many described subduction complexes (or mélanges) exhumed from seismogenic depths comprise thick, turbidite-dominated sequences with deformed zones containing clasts or boudins of more competent sandstone and/or basalt. In contrast, many active subduction zones have a relatively small thickness of sedimentary inputs (<2 km), turbidite sequences are commonly accreted rather than subducted, and the role of pelagic sediments and basalt (lavas and hyaloclastites) in the deforming zone near the plate interface at <20 km depth is poorly understood...

High-resolution computer simulations of earthquake sequences in three or even two dimensions pose great demands on time and energy, making lower-cost simplifications a competitive alternative. We systematically study the advantages and limitations of simplifications that eliminate spatial dimensions in quasi-dynamic earthquake sequence models, from 3D models with a 2D fault plane down to 0D or 1D models with a 0D fault point. We demonstrate that, when 2D or 3D models produce quasi-periodic characteristic earthquakes, their behavior is qualitatively similar to lower-dimension models...

We propose a theory for preventing instabilities and inducing controlled, slow-slip in frictionally unstable systems, such as the Generalized-Burridge-Knopoff (GBK) model and seismic fault models. We exploit the dependence of friction on pressure and use it as a backdoor for altering the dynamics of the underlying dynamical system. We use the mathematical Theory of Control and, for the first time, we manage to (a) stabilize and restrict chaos in this kind of systems, (b) guarantee slow frictional dissipation and (c) tune the system toward desirable global asymptotic equilibria of lower energy...

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick-slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa...

The permeability of volcanic rock controls the distribution of pore fluids and pore fluid pressure within a volcanic edifice, and is therefore considered to influence eruptive style and volcano deformation. We measured the porosity and permeability of a porous volcanic rock during deformation in the brittle and ductile regimes. In the brittle regime, permeability decreases by a factor of 2-6 up to the peak stress due the closure of narrow pore throats but, following shear fracture formation, remains approximately constant as strain is accommodated by sliding on the fracture...

Micromagnetic tomography (MMT) is a technique that combines X-ray micro computed tomography and scanning magnetometry data to obtain information about the magnetic potential of individual grains embedded in a sample. Recovering magnetic signals of individual grains in natural and synthetic samples provides a new pathway to study the remanent magnetization that carries information about the ancient geomagnetic field and is the basis of all paleomagnetic studies. MMT infers the magnetic potential of individual grains by numerical inversion of surface magnetic measurements using spherical harmonic expansions...

We present the first three-dimensional (3D) anisotropic teleseismic P -wave tomography model of the upper mantle covering the entire Central Mediterranean. Compared to isotropic tomography, it is found that including the magnitude, azimuth, and, importantly, dip of seismic anisotropy in our inversions simplifies isotropic heterogeneity by reducing the magnitude of slow anomalies while yielding anisotropy patterns that are consistent with regional tectonics. The isotropic component of our preferred tomography model is dominated by numerous fast anomalies associated with retreating, stagnant, and detached slab segments...

Aseismic afterslip is postseismic fault sliding that may significantly redistribute crustal stresses and drive aftershock sequences. Afterslip is typically modeled through geodetic observations of surface deformation on a case-by-case basis, thus questions of how and why the afterslip moment varies between earthquakes remain largely unaddressed. We compile 148 afterslip studies following 53 M w 6.0-9.1 earthquakes, and formally analyze a subset of 88 well-constrained kinematic models. Afterslip and coseismic moments scale near-linearly, with a median Spearman's rank correlation coefficient (CC) of 0...

Thermobarometry is a fundamental tool to quantitatively interrogate magma plumbing systems and broaden our appreciation of volcanic processes. Developments in random forest-based machine learning lend themselves to a data-driven approach to clinopyroxene thermobarometry, allowing users to access large experimental data sets that can be tailored to individual applications in Earth Sciences. We present a methodological assessment of random forest thermobarometry using the R freeware package extraTrees. We investigate the model performance, the effect of hyperparameter tuning, and assess different methods for calculating uncertainties...

Ab initio molecular dynamics calculations on a carbonate-silicate-metal melt were performed to study speciation and coordination changes as a function of pressure and temperature. We examine in detail the bond abundances of specific element pairs and the distribution of coordination environments over conditions spanning Earth's present-day mantle. Average coordination numbers increase continuously from 4 to 8 for Fe and Mg, from 4 to 6 for Si, and from 2 to 4 for C from 1 to 148 GPa (4,000 K). Speciation across all pressure and temperature conditions is complex due to the unusual bonding of carbon...

Global Navigation Satellite System (GNSS) vertical displacements measuring the elastic response of Earth's crust to changes in hydrologic mass have been used to produce terrestrial water storage change (∆TWS) estimates for studying both annual ∆TWS as well as multi-year trends. However, these estimates require a high observation station density and minimal contamination by nonhydrologic deformation sources. The Gravity Recovery and Climate Experiment (GRACE) is another satellite-based measurement system that can be used to measure regional TWS fluctuations...

Continuous seismograms contain a wealth of information with a large variety of signals with different origin. Identifying these signals is a crucial step in understanding physical geological objects. We propose a strategy to identify classes of signals in continuous single-station seismograms in an unsupervised fashion. Our strategy relies on extracting meaningful waveform features based on a deep scattering network combined with an independent component analysis. Based on the extracted features, agglomerative clustering then groups these waveforms in a hierarchical fashion and reveals the process of clustering in a dendrogram...