The tsunami-generated magnetic field is a magnetic field that show up with the moving of tsunami. In the previous studies, researchers claimed that the tsunami-generated magnetic field arrives earlier than the tsunami sea level change based on analytical solutions and numerical simulations. In this paper, we used the world's first simultaneous data of sea level change and magnetic field in the 2009 Samoa and 2010 Chile tsunamis to study the relation between these two physical quantities. We found that the vertical component of tsunami magnetic field arrives earlier than the sea level change. Moreover, the horizontal component of tsunami magnetic field arrives even earlier than the vertical component. The tsunami magnetic field was also revealed that it can be used to estimate the tsunami wave height very accurately. We investigated the observed tsunami magnetic field by our 3-D time-domain simulation. However, the currently existing tsunami source models were unable to reproduce the observation in our research area. We confirmed that a better source model can improve the simulation. It follows that our high precision tsunami wave height data converted from the magnetic field should be used to construct a better tsunami source model.

Long-period (T > 10 s) shear-wave reverberations between the surface and reflecting boundaries below seismic stations are useful for studying the mantle transition zone (MTZ) but finite-frequency effects may complicate the interpretation of waveform stacks. Using waveform data from the USArray and spectral-element method synthetics for 3-D seismic models, we illustrate that a common-reflection point (CRP) modeling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes. Our CRP mapping of recorded waveforms places the 410-km and 660-km phase boundaries about 15 km deeper beneath the western US than beneath the central-eastern US if it is based on the 1-D PREM model. The apparent east-to-west deepening of the MTZ disappears in the CRP image if we account for shear-wave velocity variations in the mantle. We also find that ray theory overpredicts the traveltime delays of the reverberations if 3-D velocity variations in the mantle are prescribed by global models S40RTS, SEMUCB-WM1, and TX2015. Undulations of the 410-km and 660-km are underestimated in the analysis when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ.

Joint probabilistic inversions of magnetotelluric (MT) and seismic data has great potential for imaging the thermochemical structure of the lithosphere as well as mapping fluid/melt pathways and regions of mantle metasomatism. In this contribution we present a novel probabilistic (Bayesian) joint inversion scheme for 3D MT and surface-wave dispersion data particularly designed for large-scale lithospheric studies. The approach makes use of a recently developed strategy for fast solutions of the 3D MT forward problem (Manassero et al., 2020) and combines it with adaptive Markov chain Monte Carlo (MCMC) algorithms and parallel-in-parallel strategies to achieve extremely efficient simulations. To demonstrate the feasibility, benefits and performance of our joint inversion method to image the conductivity, temperature and velocity structures of the lithosphere, we apply it to two numerical examples of increasing complexity. The inversion approach presented here is timely and will be useful in the joint analysis of MT and surface wave data that are being collected in many parts of the world. This approach also opens up new avenues for the study of translithospheric and transcrustal magmatic systems, the detection of metasomatised mantle and the incorporation of MT into multi-observable inversions for the physical state of the Earth’s interior.

Clinopyroxene and orthopyroxene are the two major repositories of rare earth elements (REE) in spinel peridotites. Most geochemical studies of REE in mantle samples focus on clinopyroxene. Recent advances in in situ trace element analysis has made it possible to measure REE abundance in orthopyroxene. The purpose of this study is to determine what additional information one can learn about mantle processes from REE abundances in orthopyroxene coexisting with clinopyroxene in residual spinel peridotites. To address this question, we select a group of spinel peridotite xenoliths (9 samples) and a group of abyssal peridotites (12 samples) that are considered residues of mantle melting and that have major element and REE compositions in the two pyroxenes reported in the literature. We use a disequilibrium double-porosity melting model and the Markov chain Monte Carlo method to invert melting parameters from REE abundance in the bulk sample. We then use a subsolidus reequilibration model to calculate REE redistribution between cpx and opx at the extent of melting inferred from the bulk REE data and at the closure temperature of REE in the two pyroxenes. We compare the calculated results with those observed in clinopyroxene and orthopyroxene in the selected peridotitic samples. Results from our two-step melting followed by subsolidus reequilibration modeling show that it is more reliable to deduce melting parameters from REE abundance in the bulk peridotite than in clinopyroxene. We do not recommend the use of REE in clinopyroxene alone to infer the degree of melting experienced by the mantle xenolith, as HREE in clinopyroxene in the xenolith are reset by subsolidus reequilibration. In general, LREE in orthopyroxene and HREE in clinopyroxene are more susceptible to subsolidus redistribution. The extent of redistribution depends on the modes of clinopyroxene and orthopyroxene in the sample and thermal history experienced by the peridotite. By modeling subsolidus redistribution of REE between orthopyroxene and clinopyroxene after melting, we show that it is possible to discriminate mineral mode of the starting mantle and cooling rate experienced by the peridotitic sample. We conclude that endmembers of the depleted MORB mantle and the primitive mantle are not homogeneous in mineral mode. A modally heterogeneous peridotitic starting mantle provides a simple explanation for the large variations of mineral mode observed in mantle xenoliths and abyssal peridotites. Finally, by using different starting mantle compositions in our simulations, we show that composition of the primitive mantle is more suitable for modeling REE depletion in cratonic mantle xenoliths than the composition of the depleted MORB mantle.

Connectivity of material constituents govern the transport, mechanical, chemical, thermal, and electromagnetic properties. Energy storage, recovery and conversion depends on connectivity of material constituents. High-resolution microscopy image of a material captures the microstructural aspects describing the distribution, topology and morphology of various material constituents. In this study, six metrics are developed and tested for quantifying the connectivity of material constituents as captured in the high-resolution microscopy images. The six metrics are as follows: geobody connectivity metric based on percolation theory, Euler number based on integral geometry, indicator variogram based on geostatistics, two-point cluster function, connectivity function, and travel-time histogram based on fast marching method. The performances of these metrics are tested on 3000 images representing six levels of connectivity. The metrics are also evaluated on the organic constituent captured in the scanning electron microscopy (SEM) images of organic-rich shale samples. The connectivity function and travel-time histogram based on fast-marching method are the most robust and reliable metrics. Material constituents exhibiting high connectivity result in large values of average travel time computed using fast-marching method and average connected distance computed using connectivity function. The proposed metrics will standardize and speed-up the analysis of connectivity to facilitate the characterization of properties and processes of energy-relevant materials.

Low-velocity accretionary wedges and sedimentary layers overlying continental plates widely exist in subduction zones. However, the two structures are commonly neglected in velocity models used in slip inversion, ground motion estimation, and dynamic rupture simulation, which may cause a biased estimation of coseismic slip and near-fault ground motions during subduction zone earthquakes. We use the 2011 Mw 9.0 Tohoku-Oki earthquake as an example and reproduce the observed seafloor deformation using 2-D dynamic rupture models with or without an accretionary wedge and a sedimentary layer. We find that the co-existence of the accretionary wedge and sedimentary layer significantly enhances the shallow coseismic slip and amplifies ground accelerations near the accretionary wedge. Hence, stress drop on the shallow fault estimated from the coseismic slip or surface deformation is overestimated when the two structures are neglected. We further simulate a suite of earthquakes where the up-dip rupture terminates at different depths. Results show that a sedimentary layer enhances coseismic slip in all cases, while an accretionary wedge can lead to a sharper decline in slip when negative dynamic stress drop exists on the shallow fault. However, a combination of the two structures tends to enhance fault slip, especially when rupture breaks through a trench. Thus, their combined effects are nonlinear and can be larger than the respective contribution of each structure. Our results emphasize that subduction zones featuring a co-existence of an accretionary wedge and a sedimentary layer may have inherently higher earthquake and tsunami hazards.

Seismic interferometry is a powerful tool to monitor the seismic velocity change associated with volcanic eruptions. For the monitoring, changes in seismic velocity with environmental origins (such as precipitation) are problematic. In order to model the environmental effects, we propose a new technique based on a state-space model. An extended Kalman filter estimates seismic velocity changes as state variables, with a first-order approximation of the stretching method. We apply this technique to three-component seismic records in order to detect the seismic velocity change associated with the Shinmoe-dake eruptions in 2011 and 2018. First, ambient noise cross-correlations were calculated from May 2010 to April 2018. We also modeled seismic velocity changes resulting from precipitation and the 2016 Kumamoto earthquake, with exponential type responses. Most of the results show no significant changes associated with the eruptions, although gradual inflation of the magma reservoir preceded the 2011 eruption by one year. The observed low sensitivity to static stress changes suggests that the fraction of geofluid and crack density at about 1 km depth is small, and the crack shapes could be circular. Only one station pair west of the crater shows the significant drop associated with the eruption in 2011. The gradual drop of seismic velocity up to 0.05% preceded the eruption by one month. When the gradual drop began, volcanic tremors were activated at about 2 km depth. These observations suggest that the drop could be caused by damage accumulation due to vertical magma migration beneath the summit.

Compositional heterogeneities within Europa’s ice shell likely impact the dynamics and habitability of the ice and subsurface ocean, but the total inventory and distribution of impurities within the shell is unknown. In sea ice on Earth, the thermochemical environment at the ice-ocean interface governs impurity entrainment into the ice. Here, we simulate Europa’s ice-ocean interface and bound the impurity load (1.053-14.72 g/kg (parts per thousand weight percent, or ppt) bulk ice shell salinity) and bulk salinity profile of the ice shell. We derive constitutive equations that predict ice composition as a function of the ice shell thermal gradient and ocean composition. We show that evolving solidification rates of the ocean and hydrologic features within the shell produce compositional variations (ice bulk salinities of 5-50% of the ocean salinity) that can affect the material properties of the ice. As the shell thickens, less salt is entrained at the ice-ocean interface, which implies Europa’s ice shell is compositionally homogeneous below ~ 1 km. Conversely, the solidification of water filled fractures or lenses introduces substantial compositional variations within the ice shell, creating gradients in mechanical and thermal properties within the ice shell that could help initiate and sustain geological activity. Our results suggest that ocean materials entrained within Europa’s ice shell affect the formation of geologic terrain and that these structures could be confirmed by planned spacecraft observations.

We conduct a high-resolution teleseismic receiver function investigation of the subducting plate interface within the Alaskan forearc beneath Kodiak Island using data collected as part of the Alaska Amphibious Community Seismic Experiment in 2019. The Kodiak node array consisted of 398 nodal geophones deployed at ~200-300 m spacing on northeastern Kodiak Island within the southern asperity of the 1964 Mw9.2 Great Alaska earthquake. Receiver function images at frequencies of 1.2 and 2.4 Hz show a coherent, slightly dipping velocity increase at ~30-40 km depth consistent with the expected slab Moho. In contrast to studies within the northern asperity of the 1964 rupture, we find no evidence for a prominent low-velocity layer above the slab Moho thick enough to be resolved by upgoing P-to-S conversions. These results support evidence from seismicity and geodetic strain suggesting that the 1964 rupture connected northern (Kenai) and southern (Kodiak) asperities with different plate interface properties.

Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1,640 source time functions (STFs) of Mw 5.5-8.0 earthquakes in the SCARDEC catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and therefore the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress drop estimates.

Characterizing the large M4.7+ seismic events during the 2018 Kīlauea eruption is important to understand the complex subsurface deformation at the Kīlauea summit. The first 12 events (May 17 - May 26) are associated with long-duration seismic signals and the remaining 50 events (May 29 - August 02) are accompanied by large-scale caldera collapses. Resolving the source location and mechanism is challenging because of the shallow source depth, significant non double-couple components, and complex velocity structure. We demonstrate that combining multiple geophysical data from broadband seismometers, accelerometers and infrasound is essential to resolve different aspects of the seismic source. Seismic moment tensor solutions using near-field summit stations show the early events are highly volumetric. Infrasound data and particle motion analysis identify the inflation source as the Halema’uma’u reservoir. For the later collapse events, two independent moment tensor inversions using local and global stations consistently show that asymmetric slips occur on inward-dipping normal faults along the northwest corner of the caldera. While the source mechanism from May 29 onwards is not fully resolvable seismically using far-field stations, infrasound records and simulations suggest there may be inflation during the collapse. The summit events are characterized by both inflation and asymmetric slip, which are consistent with geodetic data. Based on the location of the slip and microseismicity, the caldera may have failed in a ‘see-saw’ manner: small continuous slips in the form of microseismicity on the southeast corner of the caldera, compensated by large slips on the northwest during the large collapse events.

On 15 January 2022, the submarine volcano on the southwest Pacific island of Tonga violently erupted. Thus far, the ionospheric oscillation features caused by the volcanic eruption have not been identified. Here, observations from the Super Dual Auroral Radar Network (SuperDARN) radars and digisondes \change{are}{were} employed to analyze ionospheric oscillations in the Northern Hemisphere caused by the volcanic eruption in Tonga. Due to the magnetic field conjugate effect, the ionospheric oscillations were observed much earlier than the arrival of surface air pressure waves, and the maximum negative line-of-sight (LOS) velocity of the ionospheric oscillations exceeded 100 m/s in the F layer. After the surface air pressure waves arrived, the maximum LOS velocity in the E layer approached 150 m/s. A maximum upward displacement of 100 km was observed in the ionosphere. This work provides a new perspective for understanding the strong ionospheric oscillation caused by geological hazards observed on Earth.

Open cracks and cavities play important roles in fluid transport. Underground water penetration induces microcrack activity, which can lead to rock failure and earthquake. Fluids in cracks can affect earthquake generation mechanisms through physical and physicochemical effects. Methods for characterizing the crack shape and water saturation of underground rock are needed for many scientific and industrial applications. The ability to estimate the status of cracks by using readily observable data such as elastic-wave velocities would be beneficial. We have demonstrated a laboratory method for estimating the crack status inside a cylindrical rock sample based on a vertically cracked transversely isotropic solid model by using measured P- and S-wave velocities and porosity derived from strain data. During injection of water to induce failure of a stressed rock sample, the crack aspect ratio changed from 1/400 to 1/160 and the degree of water saturation increased from 0 to 0.6. This laboratory-derived method can be applied to well-planned observations in field experiments. The in situ monitoring of cracks in rock is useful for industrial and scientific applications such as the sequestration of carbon dioxide and other waste, induced seismicity, and measuring the regional stress field.

Subsurface sequestration of carbon dioxide (CO2) requires long-term monitoring of the injected CO2 plume to prevent CO2 leakage along the wellbore or across the caprock. Accurate knowledge of the location and movement of the injected CO2 is crucial for risk management at a geological CO2-storage complex. Conventional methods for locating/assessing the injected CO2 plume in the subsurface assume a geophysical model, which is specific and may not be applicable to all types of CO2-injection reservoirs and scenarios. We developed an unsupervised-learning-based visualization of the subsurface CO2 plume that adapts and scales based on the data without requiring an assumption of the geophysical model. The data-processing workflow was applied to the cross-well tomography data from the SECARB Cranfield carbon geo-sequestration project. A multi-level clustering approach was developed to account for data imbalance due to the absence of CO2 in the large portion of the imaged reservoir. The first level of clustering differentiated CO2-bearing regions from the non-CO2 bearing regions and achieved a silhouette score of 0.85, a Calinski-Harabasz index of 160666, and a Davies-Bouldin index of 0.43, which are indicative of high quality, reliable clustering. The second level of clustering further differentiated the CO2-bearing regions into regions containing low, medium, and high CO2 content. Overall, the multi-level clustering achieved a silhouette score, Calinski-Harabasz index, and Davies-Bouldin index of 0.74, 59656, and 0.32, which confirm the high quality and reliability of the newly proposed unsupervised-learning-based visualization. Three distinct clustering techniques, namely k-means, mean-shift, and agglomerative, generated similar visualizations. In terms of the adjusted Rand index, the similarity of clusters identified by the three distinct clustering techniques is around 0.98, which indicates the robustness of the cluster labels assigned to various regions of the CO2-injection reservoir. Further, we find certain geophysical signatures, such as Fourier transform and wavelet transform, to be highly relevant and informative indicators of the spatial distribution of CO2 content.

A volcanic eruption is usually preceded by seismic precursors, but their interpretation and use for forecasting the eruption onset time remain a challenge. A part of the eruptive processes in open conduits of volcanoes may be similar to those encountered in geysers. Since geysers erupt more often, they are useful sites for testing new forecasting methods. We tested the application of Permutation Entropy (PE) as a robust method to assess the complexity in seismic recordings of the Strokkur geyser, Iceland. Strokkur features several minute-long eruptive cycles, enabling us to verify in 63 recorded cycles whether PE behaves consistently from one eruption to the next one. We performed synthetic tests to understand the effect of different parameter settings in the PE calculation. Our application to Strokkur shows a distinct, repeating PE pattern consistent with previously identified phases in the eruptive cycle. We find a systematic increase in PE within the last 15s before the eruption, indicating that an eruption will occur. We quantified the predictive power of PE, showing that PE performs better than seismic signal strength or quiescence when it comes to forecasting eruptions.

Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid-ocean ridge dynamics. Yet, interpreting anisotropy in the context of grain-scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co-located compressional- and shear-wavespeed anisotropy observations at the NoMelt experiment located on ~70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil-spreading direction, consistent with corner-flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250–400% in the shallow mantle. We find evidence for D-type olivine lattice-preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D-type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during mid-ocean ridge spreading; we hypothesize that grain-boundary sliding during dislocation creep relaxes strain compatibility, allowing D-type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation-accommodated grain-boundary sliding (disGBS) has implications for in situ stress and grain size during mid-ocean ridge spreading and implies grain-size dependent deformation, in contrast to pure dislocation creep.

The 2019 Ridgecrest, California, earthquake sequence represents a complex pattern of seismicity that is characterized by the occurrence of a well defined foreshock sequence followed by a mainshock and subsequent aftershocks. In this work, a detailed statistical analysis of the sequence is performed. Particularly, the parametric modelling of the frequency-magnitude statistics and the earthquake occurrence rate is carried out. It is shown that the clustering of earthquakes plays an important role during the evolution of this sequence. In addition, the problem of constraining the magnitude of the largest expected aftershocks to occur during the evolution of the sequence is addressed. In order to do this, two approaches are considered. The first one is based on the extreme value theory, whereas the second one uses the Bayesian predictive framework. The latter approach has allowed to incorporate the complex earthquake clustering through the Epidemic Type Aftershock Sequence (ETAS) process and the uncertainties associated with the model parameters into the computation of the corresponding probabilities. The results indicate that the inclusion of the foreshock sequence into the analysis produces higher probabilities for the occurrence of the largest expected aftershocks after the M7.1 mainshock compared to the approach based on the extreme value distribution combined with the Omori-Utsu formula for the earthquake rate. Several statistical tests are applied to verify the forecast.

We investigate the spatiotemporal patterns of strain accommodation during multiphase rift evolution in the Shire Rift Zone (SRZ), East Africa. The NW-trending SRZ records a transition from magma-rich rifting phases (Permian-Early Jurassic: Rift-Phase 1 (RP1), and Late Jurassic-Cretaceous: Rift-Phase 2 (RP2)) to a magma-poor phase in the Cenozoic (ongoing: Rift-Phase 3 (RP3)). Our observations show that although the rift border faults largely mimic the pre-rift basement metamorphic fabrics, the rift termination zones occur near crustal-scale rift-orthogonal basement shear zones (Sanangoe (SSZ) and the Lurio shear zones) during RP1-RP2. In RP3, the RP1-RP2 sub-basins were largely abandoned, and the rift axes migrated northeastward (rift-orthogonally) into the RP1-RP2 basin margin, and northwestward (strike-parallel) ahead of the RP2 rift-tip. The northwestern RP3 rift-axis side-steps across the SSZ, with a rotation of border faults across the shear zone and terminates further northwest at another regional-scale shear zone. We suggest that over the multiple pulses of tectonic extension and strain migration in the SRZ, pre-rift basement fabrics acted as: 1) zones of mechanical strength contrast that localized the large rift faults, and 2) mechanical ‘barriers’ that refracted and possibly, temporarily halted the propagation of the rift zone. Further, the cooled RP1-RP2 mafic dikes facilitated later-phase deformation in the form of border fault hard-linking transverse faults that exploited mechanical anisotropies within the dike clusters and served as mechanically-strong zones that arrested some of the RP3 fault-tips. Overall, we argue that during pulsed rift propagation, inherited strength anisotropies can serve as both strain-localizing, refracting, and transient strain-inhibiting tectonic structures.

Lithospheric discontinuities, including the lithosphere-asthenosphere boundary (LAB) and the enigmatic mid-lithospheric discontinuities (MLDs), hold important clues about the structure and evolution of tectonic plates. However, P- and S-receiver-function techniques (PRF and SRF), two traditional techniques to image Earth’s deep discontinuities, have some shortcomings in imaging lithosphere discontinuities. Here, we propose a new method using reflections generated by teleseismic S waves (hereafter S reflections) to image lithospheric discontinuities, which is less affected by multiple phases than PRFs and has better depth resolution than SRFs. We apply this method to data collected by the Transportable Array and other regional seismic networks and obtain new high-resolution images of the lithosphere below the contiguous US. Beneath the tectonically active Western US, we observe a negative polarity reflector (NPR) in the depth range of 60–110 km, with greatly varying amplitude and depth, which correlates with active tectonic processes. We interpret this feature as the lithosphere-asthenosphere boundary below the Western US. Beneath the tectonically stable Central and Eastern US, we observe two NPRs in the depth ranges of 60–100 km and 100–150 km, whose amplitude and depth also vary significantly, and which appear to correlate with past tectonic processes. We interpret these features as mid-lithospheric discontinuities below the Central and Eastern US. Our results show reasonable agreement with results from PRFs, which have similar depth resolution, suggesting the possibility of joint inversion of S reflections and PRFs to constrain the properties of lithospheric discontinuities.

Understanding the behavior of the shallow portion of the subduction zone, which generates the largest earthquakes and devastating tsunamis, is a vital step forward in earthquake geoscience. Monitoring only a fraction of a single megathrust earthquake cycle and the offshore location of the source of these earthquakes are the foremost reasons for the insufficient understanding. The frictional-elastoplastic interaction between the interface and its overlying wedge causes variable surface strain signals such that the wedge strain patterns may reveal the mechanical state of the interface. We employ Seismotectonic Scale Modeling and simplify elastoplastic megathrust subduction, generate hundreds of analog seismic cycles at laboratory scale, and monitor the surface strain signals over the model’s forearc over high to low temporal resolutions. We establish two coseismically compressional and extensional wedge configurations to explore the mechanical and kinematic interaction between the shallow wedge and the interface. Our results demonstrate that this interaction can partition the wedge into different segments such that the anlastic extensional segment overlays the seismogenic zone at depth. Moreover, the different segments of the wedge may switch their state from compression/extension to extension/compression domains. We highlight that a more segmented upper plate represents megathrust subduction that generates more characteristic and periodic events. Additionally, the strain time series reveals that the strain state may remain quasi-stable over a few seismic cycles in the coastal zone and then switch to the opposite mode. These observations are crucial for evaluating earthquake-related morphotectonic markers (i.e., marine terraces) and short-term interseismic GPS time-series onshore (coastal region).

The Yellowstone magmatic system is one of the largest magmatic systems on Earth, and thus an ideal location to study magmatic processes. Whereas previous seismic tomography results could only image a shallow magma reservoir, a recent study using more seismometers showed that a second and massive partially molten mush reservoir exists above the Moho \citep{huang2015yellowstone}. To understand the measurable surface response of this system to visco-elasto-plastic deformation, it is thus important to take the whole system from the mantle plume up to the shallow magma reservoirs into account. Here, we employ lithospheric-scale 3D visco-elasto-plastic geodynamic models to test the influence of parameters such as the connectivity of the reservoirs and rheology of the lithosphere on the dynamics of the system. A gravity inversion is used to constrain the effective density of the magma reservoirs, and an adjoint modelling approach reveals the key model parameters affecting the surface velocity. Model results show that a combination of connected reservoirs with plastic rheology can explain the recorded slow vertical surface uplift rates of around 1.2 cm/yr, as representing a long term background signal. A geodynamic inversion to fit the model to observed GPS surface velocities reveals that the magnitude of surface uplift varies strongly with the viscosity difference between the reservoirs and the crust. Even though stress directions have not been used as inversion parameters, modelled stress orientations are consistent with observations. However, phases of larger uplift velocities can also result from magma reservoir inflation which is a short term effect. We consider two approaches: 1) overpressure in the magma reservoir in the asthenosphere and 2) inflation of the uppermost reservoir prescribed by an internal kinematic boundary condition. We demonstrate that the asthenosphere inflation has a smaller effect on the surface velocities in comparison with the uppermost reservoir inflation. We show that the pure buoyant uplift of magma bodies in combination with magma reservoir inflation can explain (varying) observed uplift rates at the example of the Yellowstone volcanic system.

Europa Clipper will arrive at Jupiter at the end of this decade and will explore Europa through a series of flybys. One of its many goals is to characterize Europa’s topography and global shape using the EIS and REASON instruments. In addition, Europa Clipper’s UV Spectrograph will observe stars pass behind (be occulted by) Europa. The spectrograph has sufficiently precise timing, corresponding to a topographic precision of order meters, that these occultations can also serve as altimetric measurements. Because of gaps in the REASON radar altimeter coverage imposed by the flyby geometries, the addition of ~100 occultations results in a substantial improvement in the recovery of Europa’s long-wavelength shape. Typically five extra spherical harmonic degrees of topography can be recovered by combining occultations with radar altimetry.

This paper analyzes the effects of pore pressure rate for a spring - block system that is a simple model of a laboratory experiment. Pore pressure is increased at a constant rate in a remote reservoir and slip is governed by rate and state friction. The frequency of rapid slip events increases with the increase of a nondimensional pressure rate that is the ratio of the time scale of frictional sliding to that for pressure increase. Rate and state and pressure rate effects interact in a limited range of pressure rate and diffusivity. Above a critical value of the pressure rate there is transition to a significant downward linear trend of the stress, reflecting the increase of pore fluid pressure in the reservoir. This trend leads to Coulomb failure due to the decrease of the frictional resistance and the effective stress principle.

Declaration (Stauning, 2022) on “invalid PCS index” is based on the following arguments: PCS index is calculated with use of incorrect procedure (Janzhura and Troshichev, 2011), which is a part of the unified PC derivation method (Troshichev et al., 2006); PCS index used in analyses is a preliminary index, which was not approved by IAGA and, therefore, it cannot be regarded a correct index; PCN and PCS indices demonstrate, from time to time, large difference in value, which should be treated as evidence of the PCS index invalidity. This paper presents comments to these arguments. (1) Procedure (Janzhura and Troshichev, 2011) was never included in the unified PC derivation method (Troshichev et al., 2006), which was used for calculation of “correct” PCN index and “invalid” PCS index”. (2) IAGA approved the PC index and method of the index derivation in 2013, when the definitive PCN and PCS indices have not been produced. Definitive PCN and PCS indices were approved by IAGA in 2021. (3) Difference between PCN and PCS indices reached the maximum in solstice periods is natural phenomenon related to the PC seasonal variation (PCsummer > PCwinter) explained by peculiarities of the field-aligned currents closure in the winter and summer polar caps. Conclusion is made that criticism of the PCS index, presented in Stauning (2022), is based on groundless arguments.