Slow Earthquakes Research Articles

Slow periodic oscillation without radiation damping: new evolution laws for rate and state friction

SUMMARY The dynamics of sliding friction is mainly governed by the frictional force. Previous studies have shown that the laboratory-scale friction is well described by an empirical law stated in terms of the slip velocity and the state variable. The state variable represents the detailed physicochemical state of the sliding interface. Despite some theoretical attempts to derive this friction law, there has been no unique equation for time evolution of the state variable. Major equations known to date have their own merits and drawbacks. To shed light on this problem from a new aspect, here we investigate the feasibility of periodic motion without the help of radiation damping. Assuming a patch on which the slip velocity is perturbed from the rest of the sliding interface, we prove analytically that three major evolution laws fail to reproduce stable periodic motion without radiation damping. Furthermore, we propose two new evolution equations that can produce stable periodic motion without radiation damping. These two equations are scrutinized from the viewpoint of experimental validity and the relevance to slow earthquakes.

Empirical Low‐Frequency Earthquakes Synthesized From Tectonic Tremor Records

AbstractTectonic tremor and low‐frequency earthquakes (LFEs) are two different representations of the high‐frequency (>1 Hz) components of broadband slow earthquakes, which have been discovered in various tectonic regions. Although LFEs are considered building blocks of tremor, it is difficult to find constituent LFEs for some long‐duration tremor sequences. Here we introduce a new scheme to synthesize an impulsive response from a tremor source, or empirical LFE (eLFE), from complex tremor sequences. The eLFE is constructed by maximizing the absolute amplitude (L∞ norm) of the normalized and stacked waveforms of many tremor sequences at many stations. This problem includes a highly multimodal objective function, which is maximized via a combined optimization method that incorporates a genetic algorithm and iterative improvements. We first apply this method to several validation tests using synthetic tremor waveforms and real LFEs, and then tremor activities where few nearby LFEs have been detected. One example is the Okayama tremor, which has been known to be extraordinarily sensitive to tidal stress. The eLFE is helpful in determining the accurate hypocentral location and recovering the long‐term comprehensive tremor activity via a matched‐filter analysis. The well‐constrained depth of the Okayama tremor sources is less than 30 km, which suggests that the subducting Philippine Sea Plate is very shallow at this distant location from the trench axis. This eLFE scheme highlights the fact that tremor sequences may be composed of random fluctuations that are not necessarily impulsive and emphasizes the importance of characterizing such randomness to better understand slow earthquake behavior.

Thermal regime and slab dehydration beneath the Izu‐Bonin arc: Implications for fast and slow subduction earthquakes

AbstractThe transition from cold forearc to hot arc is considered fundamental to subduction dynamics. However, the range and depth of the cold nose changing significantly along the strike are difficult to constrain. By incorporating the three‐dimensional high‐resolution slab topography and MORVEL plate motion field to constrain thermomechanical modelling, we estimate the subduction thermal state and petrological dehydration in Izu‐Bonin. The multilayered hydrothermal regime gradually transitions from cold to less cold forearc, associated with a large‐scale subduction dewatering process varying along the strike. The dehydration embrittlement resulting from eclogitization and harzburgitization influences the occurrence of fast to slow earthquakes on megathrusts. Impressively, the cold nose is constrained to a depth of 60 km between the arc and Moho‐depth plate interface. The lowered temperature and delay of rock harzburgitization in the cold center of the subducted plate at depth contribute to the occurrence of deep earthquakes below the hot backarc interface.

Slow slip in subduction zones: Reconciling deformation fabrics with instrumental observations and laboratory results

AbstractCataclasites are a characteristic rock type found in drill cores from active faults as well as in exposed fossil subduction faults. Here, cataclasites are commonly associated with evidence for pervasive pressure solution and abundant hydrofracturing. They host the principal slip of regular earthquakes and the family of so-called slow earthquakes (episodic slip and tremor, low to very low frequency earthquakes, etc.). Slip velocities associated with the formation of the different types of cataclasites and conditions controlling slip are poorly constrained both from direct observations in nature as well as from experimental research. In this study, we explore exposed sections of subduction faults and their dominant microstructures. We use recently proposed constitutive laws to estimate deformation rates, and we compare predicted rates with instrumental observations from subduction zones. By identifying the maximum strain rates using fault scaling relations to constrain the fault core thickness, we find that the instrumental shear strain rates identified for the family of “slow earthquakes” features range from 10−3s−1 to 10−5s−1. These values agree with estimated rates for stress corrosion creep or brittle creep possibly controlling cataclastic deformation rates near the failure threshold. Typically, pore-fluid pressures are suggested to be high in subduction zones triggering brittle deformation and fault slip. However, seismic slip events causing local dilatancy may reduce fluid pressures promoting pressure-solution creep (yielding rates of <10−8 to 10−12s−1) during the interseismic period in agreement with dominant fabrics in plate interface zones. Our observations suggest that cataclasis is controlled by stress corrosion creep and driven by fluid pressure fluctuations at near-lithostatic effective pressure and shear stresses close to failure. We posit that cataclastic flow is the dominant physical mechanism governing transient creep episodes such as slow slip events (SSEs), accelerating preparatory slip before seismic events, and early afterslip in the seismogenic zone.

Autonomous extraction of millimeter-scale deformation in InSAR time series using deep learning

Systematically characterizing slip behaviours on active faults is key to unraveling the physics of tectonic faulting and the interplay between slow and fast earthquakes. Interferometric Synthetic Aperture Radar (InSAR), by enabling measurement of ground deformation at a global scale every few days, may hold the key to those interactions. However, atmospheric propagation delays often exceed ground deformation of interest despite state-of-the art processing, and thus InSAR analysis requires expert interpretation and a priori knowledge of fault systems, precluding global investigations of deformation dynamics. Here, we show that a deep auto-encoder architecture tailored to untangle ground deformation from noise in InSAR time series autonomously extracts deformation signals, without prior knowledge of a fault’s location or slip behaviour. Applied to InSAR data over the North Anatolian Fault, our method reaches 2 mm detection, revealing a slow earthquake twice as extensive as previously recognized. We further explore the generalization of our approach to inflation/deflation-induced deformation, applying the same methodology to the geothermal field of Coso, California.

The Role of Along‐Fault Dilatancy in Fault Slip Behavior

AbstractEarthquakes result from fast slip that occurs along a fault surface. Interestingly, numerous dense geodetic observations over the last two decades indicate that such dynamic slip may start by a gradual unlocking of the fault surface and related progressive slip acceleration. This first slow stage is of great interest, because it could define an early indicator of a devastating earthquake. However, not all slow slip turns into fast slip, and sometimes it may simply stop. In this study, we use a numerical model based on the discrete element method to simulate crustal strike‐slip faults of 50 km length that generate a wide variety of slip‐modes, from stable‐slip, to slow earthquakes, to fast earthquakes, all of which show similar characteristics to natural cases. The main goal of this work is to understand the conditions that allow slow events to turn into earthquakes, in contrast to those that cause slow earthquakes to stop. Our results suggest that fault surface geometry and related dilation/contraction patterns along strike play a key role. Slow earthquakes that initiate in large dilated regions bounded by neutral or low contracted domains, might turn into earthquakes. Slow events occurring in regions dominated by closely spaced, alternating, small magnitude dilational and contractional zones tend not to accelerate and may simply stop as isolated slow earthquakes.

Modelling the seismic potential of the Indo-Burman megathrust

The Indo-Burman arc is the boundary between the India and Burma plates, north of the Sumatra–Andaman subduction zone. The existence of active subduction in the Indo-Burman arc is a debatable issue because the Indian plate converges very obliquely beneath the Burma plate. Recent GPS measurements in Bangladesh, Myanmar, and northeast India indicate 13–17 mm/y of plate convergence along a shallow dipping megathrust while most of the strike-slip motion occurs on several steep faults, consistent with patterns of strain partitioning at subduction zones. A short period of instrumentally recorded seismicity and sparse historical records are insufficient to assess the possibility of great earthquakes at the Indo-Burman megathrust. Using the advantage of the Block-and-Fault Dynamics model allowing simultaneous simulation of slow tectonic motions and earthquakes, we test the hypothesis whether the India-Burma detachment is locked and able to produce great earthquakes, or it slips aseismically? We have shown that the model of locked detachment is preferred because it more adequately reproduces observed tectonic velocities. The integral characteristics of synthetic seismicity, the earthquake size distribution, and the rate of seismic activity are consistent with those derived from observations. Our results suggest that the megathrust is locked and can generate great M8+ earthquakes. The estimated average return period of great events exceeds one thousand years. Earthquakes of this size pose a great threat to NE India, Bangladesh and Myanmar, the most densely populated areas of the world.

P‐ and S‐Wave Velocities of Exhumed Metasediments From the Alaskan Subduction Zone: Implications for the In Situ Conditions Along the Megathrust

AbstractThe in situ state and properties of sediments entrained along subduction megathrusts exert key controls on their mechanics and slip behavior. Low seismic velocity and high Vp/Vs are hypothesized to indicate highly elevated fluid pressure, and are invoked as conditions in the source areas of slow earthquakes and tremor. We report on Vp and Vs measurements for exhumed metasediments from Kodiak Island, AK, representative of materials along the modern megathrust. Our data reveal anisotropy of ∼8–28% in Vp and ∼6.5–8% in Vs at effective stresses ranging from ∼1 to 90 MPa, with lower wavespeeds perpendicular to the dominant fabric. The fabric‐normal velocities are sufficiently low to explain observations from regional geophysical surveys and are consistent with rock physics‐based models that incorporate small (<1%) crack porosity. We suggest that low velocity at ∼8–20 km depth along megathrusts may arise simply from the presence of foliated metasediments, without requiring near‐lithostatic fluid pressure.

Heterogeneous Sediment Input at the Nankai Trough Subduction Zone: Implications for Shallow Slow Earthquake Localization

AbstractSubducted sediment plays a key role in modulating pore pressure and seismic behavior at subduction zones. We investigated the seismic character of incoming sediments to test how sediment and basement variations relate to the along‐strike changes within the accretionary prism and plate boundary conditions at the Nankai Trough. High‐resolution seismic data reveal for the first time the presence of countourite mounded drifts in the Shikoku Basin. These features have probably introduced permeability heterogeneities into an otherwise homogenous mud‐dominant unit. Additionally, we found that normal faults in this unit are more extensive than previously documented, which probably enhances along‐strike fluid transport. The wedge taper is more correlative with the thickness of the mud‐dominant facies than the turbidite thickness. This may be due to the permeability heterogeneities associated with contourite deposits and normal faults, or due to the absence of thick turbidite deposits. Turbidite deposits can either aid the drainage of the margin where they are not confined by basement topography, or contribute to high pore fluid pressures where they are confined by less permeable mudstone or basement topography. When confined turbidite deposits and contourite mounded drifts are subducted, they may contribute to localized compartments of excess pore pressure which provide the necessary conditions for slow slip behavior. We determined that along‐strike variations in seismic behavior are likely related to the subducting basement topographic and sediment characteristics, which vary on a local (<10 km) scale.

Shallow Slow Earthquake Episodes Near the Trench Axis off Costa Rica

AbstractSlow earthquakes are mainly distributed in regions surrounding seismogenic zones along the plate boundaries of subduction zones. In the Central American subduction zone, large regular interplate earthquakes with magnitudes of 7–8 occur repeatedly around the Nicoya Peninsula, in Costa Rica, and a tsunami earthquake occurred off Nicaragua, just north of Costa Rica, in 1992. To clarify the spatial distribution of various slip behaviors at the plate boundary, we detected and located very low frequency earthquakes (VLFEs) around the Nicoya Peninsula using a grid‐search matched‐filter technique with synthetic templates based on a regional three‐dimensional model. VLFEs were active in September 2004 and August 2005, mainly near the trench axis, updip of the seismogenic zone. The distribution of VLFEs overlapped with large slip areas of slow slip events. Low frequency tremor signals were also found in high‐frequency seismogram envelopes within the same time windows as detected VLFEs; thus, we also investigated the energy rates of tremors accompanied by VLFEs. The range of scaled energy, which is the ratio of the seismic energy rate of a tremor to the seismic moment rate of accompanying VLFE and related to the rupture process of seismic phenomena, was 10−9–10−8. The along‐dip separation of shallow slow and large earthquakes and the range of the scaled energy off Costa Rica are similar to those in shallow slow earthquakes in Nankai, which shares a similar thermal structure along the shallow plate boundary.

Constraints on Element Mobility During Deformation Within the Seismogenic Zone, Shimanto Belt, Japan

AbstractSubduction interfaces exhibit a variety of slip behaviors, including megathrust and slow earthquakes. Field observations are consistent with crack‐seal deformation, in which tensile cracks are sealed by fluid‐transported solute. However, there are few constraints on the mass fluxes and length and time scales of such deformation and attendant increases in cohesion within the seismogenic zone. Here, we present a systematic geochemical investigation of mass transport associated with development of crack‐seal veins in the Shimanto Belt—an accretionary complex that preserves a record of plate boundary slip behavior at temperatures relevant to the seismogenic zone: 150–350°C. These mélanges show evidence for shear across decameter‐scale zones of deformation dominated by anastomosing scaly fabrics and pervasive veins. We use meso‐ and microstructural observations with geochemical observations of scaly fabrics and crack‐seal veins to evaluate the role of silica redistribution in healing fracture porosity along the plate interface and modulating slip behavior. Crack‐seal veins contain primarily quartz with albite and calcite, and vein textures provide evidence for partial sealing. Bulk rock analyses determined that the amount of phyllosilicates, specifically illite and chlorite, increases with temperature. A simple mass‐balance model based on the immobile chemical component TiO2 shows a systematic trend in mobility for all three mélanges and increased element mobility as a function of temperature. Scaly fabrics and veins show compositional evidence for locally sourced mineral redistribution. This study supports a model where development of scaly slip surfaces and fracture healing through temperature‐dependent mineral redistribution can impact slip behavior and fluid flow along the subduction plate interface.

Development of a Compact Broadband Ocean-Bottom Seismometer

Abstract Studies of very-low-frequency earthquakes and low-frequency tremors (slow earthquakes) in the shallow region of plate boundaries need seafloor broadband seismic observations. Because it is expected that seafloor spatially high-density monitoring requires numerous broadband sensors for slow earthquakes near trenches, we have developed a long-term compact broadband ocean-bottom seismometer (CBBOBS) by upgrading the long-term short-period ocean-bottom seismometer that has seismic sensors with a natural frequency of 1 Hz and is being mainly used for observation of microearthquakes. Because many long-term ocean-bottom seismometers with short-period sensors are available, we can increase the number of broadband seafloor sensors at a low cost. A short-period seismometer is exchanged for a compact broadband seismometer with a period of 20 or 120 s. Because the ocean-bottom seismometers are installed by free fall, we have no attitude control during an installation. Therefore, we have developed a new leveling system for compact broadband seismic sensors. This new leveling system keeps the same dimensions as the conventional leveling system for 1 Hz seismometers so that the broadband seismic sensor can be installed conveniently. Tolerance for leveling is less than 1°. A tilt of up to 20° is allowed for the leveling operation. A microprocessor controls the leveling procedure. Some of the newly developed ocean-bottom seismometers were deployed in the western Nankai trough, where slow earthquakes frequently occur. The data from the ocean-bottom seismometers on the seafloor were evaluated, and we confirmed that the long-term CBBOBS is suitable for observation of slow earthquakes. The developed ocean-bottom seismometer is also available for submarine volcanic observation and broadband seafloor observation to estimate deep seismic structures.

Shallow slow earthquakes to decipher future catastrophic earthquakes in the Guerrero seismic gap

The Guerrero seismic gap is presumed to be a major source of seismic and tsunami hazard along the Mexican subduction zone. Until recently, there were limited observations at the shallow portion of the plate interface offshore Guerrero, so we deployed instruments there to better characterize the extent of the seismogenic zone. Here we report the discovery of episodic shallow tremors and potential slow slip events in Guerrero offshore. Their distribution, together with that of repeating earthquakes, seismicity, residual gravity and bathymetry, suggest that a portion of the shallow plate interface in the gap undergoes stable slip. This mechanical condition may not only explain the long return period of large earthquakes inside the gap, but also reveals why the rupture from past M < 8 earthquakes on adjacent megathrust segments did not propagate into the gap to result in much larger events. However, dynamic rupture effects could drive one of these nearby earthquakes to break through the entire Guerrero seismic gap.

Seismogenic and tremorgenic slow slip near the stability transition of frictional sliding

Slow-slip events and tremors occur below the seismogenic zone of major plate boundaries. While the physics of aseismic slow-slip events is relatively well understood, the mechanics of seismogenic slow slip remains elusive because the conditions leading to slow or fast ruptures are thought to be mutually exclusive. Here, we explore fault dynamics in the parametric space of frictional conditions to show that seismogenic slow-slip events are the natural behavior of homogeneous faults, as long as the velocity dependence approaches velocity neutral with a small characteristic nucleation size. Tremors can originate from rapid bursts of slow earthquakes that are triggered as the slow-slip rupture spreads over small-scale asperities. The near velocity-neutral conditions explain the underlying mechanics of collocated slow and fast slip of seismogenic slow-slip events commonly found below the seismogenic zone. The presence of material heterogeneity may explain the spatio-temporal clustering and migration features of tremor activity.

Data‐Driven Clustering Reveals More Than 900 Small Magnitude Slow Earthquakes and Their Characteristics

AbstractSmall magnitude slow earthquakes remain largely undetected in geodetic data due to noise levels. However, tremor and low‐frequency earthquakes (LFE) may manifest slowly slipping fault motion as a cluster of events, i.e., a slow earthquake. Here, we identify &gt;900 slow earthquakes in southwest Japan via data‐driven clustering of tremor and LFE catalogs. We establish a more complete database for slow earthquakes in southwest Japan and demonstrate their characteristics and long‐term behavior. While sometimes sub‐episodes of well‐known episodic tremor and slip, the small slow earthquake clusters share similar scaling properties–energy, duration, and rupture area–with larger magnitude fast and slow earthquakes. The small slow earthquake clusters tend to rupture faster when migrating in their preferred rupture direction, but inter‐event rupture speeds are on average similar to those of larger slow earthquakes. This suggests that rupture speed does not necessarily control slow earthquake spatial extent or eventual magnitude.

Crustal and Upper‐Mantle Structure Below Central and Southern Mexico

AbstractThe Mexican subduction zone is a complex convergent margin characterized by a slab with an unusual morphology, an abnormal location of the volcanic arc, and the existence of slow earthquakes (slow slip events and tectonic tremors). The number of seismic imaging studies of the Mexican subduction zone has increased over the past decade. These studies have revealed the seismic wave velocity structure beneath Central and Southern Mexico. However, the existing tomographic models are either confined to specific areas or have coarse resolution. As a consequence, they do not capture the lateral and finer‐scale variations in the mantle and the crust structure. Here, we fit frequency‐dependent traveltime differences between observed and synthetic seismograms in a three‐dimensional model of Central and Southern Mexico to constrain crustal and upper mantle seismic velocity structures jointly. We use ∼3,300 seismic records, filtered between 5 and 50 s, from 74 regional earthquakes. We perform point spread function tests to assess the resolution and trade‐offs between model parameters. Our tomographic model reveals several seismic wave velocity features. These features include the geometry of the Cocos slab, the partial melting zone beneath the Trans‐Mexican Volcanic Belt, and the Yucatan slab diving below Los Tuxtlas Volcanic Field that are consistent with previous studies. We also identify horizontal structures in the lower crust and the shallow upper mantle that were not previously identified. These structures include slow seismic wave velocity anomalies that may be correlated with ultraslow seismic velocity zones and high conductivity regions, where slow earthquakes have been identified.

High Fluid‐Pressure Patches Beneath the Décollement: A Potential Source of Slow Earthquakes in the Nankai Trough off Cape Muroto

AbstractPore pressure plays a key role in the generation of earthquakes in subduction zones. However, quantitative constraints for its determination are quite limited. Here, we estimate the subsurface pore pressure by analyzing the transient upwelling flow of drilling mud from borehole C0023A of the International Ocean Discovery Program (IODP) Expedition 370, in the Nankai Trough off Cape Muroto. This upward flow provided the first direct evidence of an overpressured aquifer in the underthrust sediments off Cape Muroto. To estimate the pre‐drilling pore pressure in the overpressured aquifer around a depth of 950–1,050 m below sea floor, we examined the measured porosities of core samples retrieved from nearby IODP wells; we then proceeded to explain the observed time evolution of the flow rate of the upwelling flow by modeling various sized aquifers through solving a radial diffusion equation. It was observed that for a permeability of 10−13 m2, the aquifer possessed an initial excess pore pressure of ∼5–10 MPa above the hydrostatic pressure, with a lateral dimension of several hundred meters and thickness of several tens of meters. The overpressure estimates from the porosity‐depth profile at Site C0023 differ from those at other drill sites in the region, suggesting the possible existence of multiple overpressured aquifers with a patchy distribution in the underthrust sediments of the Nankai Trough. As pore pressure is relevant in maintaining fault stability, the overpressured aquifers may be the source of slow earthquakes that have been observed around the drilling site.

3-D thermal regime and dehydration processes around the regions of slow earthquakes along the Ryukyu Trench

Several interplate seismic events, such as short-term slow slip events (S-SSEs) and low-frequency earthquakes (LFEs), have been identified in the Ryukyu Trench, southwestern Japan. As one of the specific characteristics of this seismicity, the depths at which S-SSEs occur at the plate interface beneath Okinawa Island are approximately 5–10 km shallower than those beneath the Yaeyama Islands. To elucidate the cause of this difference in depth, we constructed a three-dimensional, Cartesian thermomechanical subduction model and applied the subduction history of the Philippine Sea (PHS) plate in the model region. As a result, the interplate temperatures at which S-SSEs take place were estimated to range from 350 to 450 °C beneath Okinawa Island and from 500 to 600 °C beneath the Yaeyama Islands. The former temperature range is consistent with previous thermal modelling studies for the occurrence of slow earthquakes, but the latter temperature range is by approximately 150 °C higher than the former. Therefore, explaining how the depth difference in S-SSEs could be caused from the aspect of only the thermal regime is difficult. Using phase diagrams for hydrous minerals in the oceanic crust and mantle wedge, we also estimated the water content distribution on and above the plate interface of the PHS plate. Near the S-SSE fault planes, almost the same amount of dehydration associated with phase transformations of hydrous minerals from blueschist to amphibolite and from amphibolite to amphibole eclogite within the oceanic crust were inferred along Okinawa Island and the Yaeyama Islands, respectively. On the other hand, the phase transformations within the mantle wedge were inferred only beneath the Yaeyama Islands, whereas no specific phase transformation was inferred beneath Okinawa Island around the S-SSE occurrence region. Therefore, we conclude that dehydrated fluid derived from the oceanic crust at the plate interface would play a key role in the occurrence of S-SSEs.

A Decade of Lessons Learned from the 2011 Tohoku‐Oki Earthquake

AbstractThe 2011 Mw 9.0 Tohoku‐oki earthquake is one of the world's best‐recorded ruptures. In the aftermath of this devastating event, it is important to learn from the complete record. We describe the state of knowledge of the megathrust earthquake generation process before the earthquake, and what has been learned in the decade since the historic event. Prior to 2011, there were a number of studies suggesting the potential of a great megathrust earthquake in NE Japan from geodesy, geology, seismology, geomorphology, and paleoseismology, but results from each field were not enough to enable a consensus assessment of the hazard. A transient unfastening of interplate coupling and increased seismicity were recognized before the earthquake, but did not lead to alerts. Since the mainshock, follow‐up studies have (1) documented that the rupture occurred in an area with a large interplate slip deficit, (2) established large near‐trench coseismic slip, (3) examined structural anomalies and fault‐zone materials correlated with the coseismic slip, (4) clarified the historical and paleoseismic recurrence of M∼9 earthquakes, and (5) identified various kinds of possible precursors. The studies have also illuminated the heterogeneous distribution of coseismic rupture, aftershocks, slow earthquakes and aseismic afterslip, and the enduring viscoelastic response, which together make up the complex megathrust earthquake cycle. Given these scientific advances, the enhanced seismic hazard of an impending great earthquake can now be more accurately established, although we do not believe such an event could be predicted with confidence.

Transient motion of the largest landslide on earth, modulated by hydrological forces

Sea-level rise of the Caspian Sea (CS) during the early Khvalynian (approximately 40–25 ka BP) generated hundreds of giant landslides along the sea’s ancient coastlines in western Kazakhstan, which extended hundreds of kilometers. Although similar landslides have been observed along the present-day coastlines of the CS in the area of a prominent high escarpment, it remains unclear whether some of these ancient landslides are still active and whether the movement is slow or catastrophic, as previously suggested. The present study is the first to show evidence proving that the geomorphic responses to sea-level changes of the CS that were triggered in the Pleistocene are currently active. Using interferometric synthetic aperture radar (InSAR) data, we show that one of these giant landslides occurring along the western shore of the Kara-Bogaz-Gol (KBG) lagoon of the CS presents active transient motion, which makes it the world’s largest active landslide reported thus far. Extending more than 25 km along the eastern coast of the inundated KBG depression in a N–S direction with maximum landward expansion of 5 km from the shoreline to the flat Ustyurt Plateau, this landslide conveys ~ 10 × 109 m3 rocks toward the lagoon at a rate of ~ 2.5 cm/year. This event releases a nearly episodic aseismic moment of 6.0 × 1010 Nm annually, which is equivalent to the response of an Mw 5.1 earthquake. We analyze the present-day evolution of this giant coastal landslide at high temporal and spatial resolutions using Sentinel-1 radar images acquired on descending and ascending modes every 12 days between 2014 and 2020. Modelling with elastic dislocations suggests that the KBG landslide was accommodated mostly by a shallow basal décollement with a nearly horizontal listric slip plane. Moreover, our analysis reveals week-long accelerating slip events at changing amplitudes that occur seasonally with slow, lateral spreading rather than sudden catastrophic motion. A strong correlation between the episodic slip events and seasonal water-level changes in the KBG lagoon suggests a causative mechanism for the transient accelerating slip events. Although water-level changes are widely acknowledged to trigger transient motion on a land mass, such movement, which is similar to a silent earthquake, has not been observed thus far at this mega scale; on an extremely low-angle detachment planes at < 5° with modulation by sea-level changes. This study suggests that present-day sea-level changes can reactivate giant landslides that originated 40–25 ka.