Fast and slow earthquakes are predominantly generated along faults constituting active plate boundaries. Characterized by repeated devastating earthquakes and frequent slow slip events and tremors, the Alaska megathrust presents a chance to understand the complicated dynamics of a subduction system changing from steep to shallow dips associated with enigmatically abundant fast and slow seismic events. Based on three-dimensional thermal modeling, we find that the downgoing metamorphosed oceanic crust containing bound water releases a large amount of fluid and causes the recurrence of fast and slow earthquakes by elevated pore fluid pressure and hydrofracturing. The seismogenic interface and the slow slip events (SSEs) identified beneath the Upper Cook Inlet coincide well with the slab metamorphic dehydration regions. The observed slow earthquakes with quasi-stable fault slips preferentially occur, accompanied by high dehydration and temperature downdip along the transition zone.

Single-stationed geomagnetic vertical intensity polarization (GVIP) anomalies have demonstrated good predictions of the occurrence of large earthquakes in Japan. Nonetheless, due to the lack of a previously densified geomagnetic network, how the multistationary GVIP anomaly (MGVIPA) corresponds to impending earthquakes remains poorly understood. Based on the newly constructed geomagnetic network from 2014 in Qinghai, China, which is composed of 23 electromagnetic stations, we suggested an MGVIPA method to analyze the correlation with large earthquakes since 2015. The results show that (1) the occurrence of MGVIPA is characterized by clusters in time that continue in a short period; (2) the spatial distribution of MGVIPA usually occurs with high values synchronously at several places over the same period; and (3) the Mw ≥ 6 earthquakes occurred in the regions indicated by MGVIPA within a period ranging from 3 months to 1 year from 2015 to 2021 in Qinghai, China.

Himalayan orogenesis remains enigmatic in terms of Tibetan Plateau geodynamics originating from the Cenozoic India–Eurasian continental collision. India underthrusts below Tibet to the Yarlung–Tsangpo suture, which has been identified as the northernmost boundary for underplating. However, the way in which the historical evolution of continental subduction induces plateau uplift and the way it controls the variation in uplift between outboard and inboard areas is still unclear. To interpret the evolutionary mechanisms involved in the Himalayan growth history, we constructed different 3-D dynamic models at important stages to address these questions related to the formation of the Himalayas on the basis of paleoenthalpy evidence encoded in fossil leaves from recently documented assemblages in southern Tibet. The results show that (1) the effect of crustal thickening was the predominant factor in the early evolution from the Paleocene to the early Eocene, which resulted in a moderate growth rate. (2) The consecutive slab break-off eastward from the western syntaxis and the associated slab rebound significantly accelerated orogenesis from the late Eocene to the Oligocene. The upwelling asthenospheric flow was a key control of increasing crustal buoyancy, which resulted in the fastest growth of the Himalayas during the early Miocene. (3) Thereafter, the gradually enhanced monsoon and surface erosion during accompanying the increasing mountain height resulted in a slowdown of the orogenic rate, which counterbalanced the buoyant force produced by asthenospheric flow driving continuous Himalayan growth.

Long-term stress accumulation influenced by coseismic stress changes and postseismic viscoelastic relaxation is considered critical to triggering giant earthquakes. Nevertheless, how the stress increase is interrupted by aftershocks and how it influences the megaseismic cycle remain enigmatic. In this study, based on the Mohr–Coulomb failure criterion at the nucleated segments of the 2008 great Sichuan earthquake, the stress variation associated with four M > 6 aftershocks was calculated for the period from 2010 to 2017. The results show that (1) the spatial distribution of coseismic stress change is correlated with the rupture pattern of large events and has a fundamental impact on triggering subsequent earthquakes and (2) postseismic viscoelastic relaxation leads to increased Coulomb stress accumulation at the northern and southern edges of the seismogenic Longmenshan fault, which results in enhanced fault instability and the potential for future large events.

Compelling evidence has shown that geomagnetic disturbances in vertical intensity polarization before great earthquakes are promising precursors across diverse rupture conditions. However, the geomagnetic vertical intensity polarization method uses the spectrum of smooth signals, and the anomalous waveforms of seismic electromagnetic radiation, which are basically nonstationary, have not been adequately considered. By combining pulse amplitude analysis and an experimental study of the cumulative frequency of anomalies, we found that the pulse amplitudes before the 2022 Luding M6.8 earthquake show characteristics of multiple synchronous anomalies, with the highest (or higher) values occurring during the analyzed period. Similar synchronous anomalies were observed before the 2021 Yangbi M6.4 earthquake, the 2022 Lushan M6.1 earthquake and the 2022 Malcolm M6.0 earthquake, and these anomalies indicate migration from the periphery toward the epicenters over time. The synchronous changes are in line with the recognition of previous geomagnetic anomalies with characteristics of high values before an earthquake and gradual recovery after the earthquake. Our study suggests that the pulse amplitude is effective for extracting anomalies in geomagnetic vertical intensity polarization, especially in the presence of nonstationary signals when utilizing observations from multiple station arrays. Our findings highlight the importance of incorporating pulse amplitude analysis into earthquake prediction research on geomagnetic disturbances.

Geomagnetic vertical intensity polarization is a method with a clear mechanism, mature processing methods, and a strong ability to extract anomalous information in the quantitative analysis of seismogenic geomagnetic disturbances. The existing analyses of geomagnetic vertical intensity polarization are all based on the 5~100 s frequency band without refinement of the partitioning process. Although many successful results have been obtained, there are still two problems in the process of extracting anomalies: the geomagnetic anomalies that satisfy the determination criteria are still high in occurrence frequency; and the anomalies are distributed over too large an area in space, which leads to difficulties in determining the location of the epicenter. In this study, based on observations from western China, where fluxgate observation points are positioned in areas with frequent, densely distributed medium-strength earthquakes, we refined the frequency bands of geomagnetic vertical intensity polarization, recalculated the spatial and temporal evolution characteristics of geomagnetic disturbances before earthquakes, and improved the crossover frequency anomaly prediction index while promoting the application of the method in earthquake forecasting.