⚡ electrical engineering

Mechanism-Based Screening of Seismic Instability in the Tarbela Reservoir Sediment Delta Using a Coupled QUAKE/W–SLOPE/W Analysis

This study employs a coupled QUAKE/W–SLOPE/W analysis to demonstrate that the Tarbela Reservoir's sediment delta is most vulnerable to localized seismic instability at the foreset toe and bottomset contact rather than generalized liquefaction, identifying a critical hazard pathway that supports targeted mitigation strategies like underwater dikes and dredging.

Original authors: Saad Hanif, Abdullah Arshad, Muhammad Adnan Hanif, Qudeer Hussain

Published 2026-07-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Saad Hanif, Abdullah Arshad, Muhammad Adnan Hanif, Qudeer Hussain

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Growing Mountain Underwater

Imagine the Tarbela Dam in Pakistan as a giant bathtub. Over the decades, the river flowing into it has been dumping sand, silt, and mud. Instead of just filling the bottom, this sediment has built a massive underwater mountain that is slowly growing toward the dam's water intake tunnels (the "drains" of the bathtub).

This underwater mountain isn't a solid block of rock; it's more like a layered cake made of loose, wet sand and mud. The researchers were worried: What happens if a big earthquake shakes this underwater mountain? Could it slide off and clog the dam's tunnels, stopping the water flow?

The Problem: It's Not Just "Liquefaction"

Usually, when people think of underwater sand and earthquakes, they think of liquefaction—where the ground turns into a thick soup (like quicksand) and everything floats away.

The researchers found that for Tarbela, it's not quite that simple. The main body of the underwater mountain is actually quite old and stiff (like a hardened crust). It's not likely to turn into soup. However, the very front edge of the mountain (where the new sediment is piling up) is loose, wet, and unstable.

The Analogy: Think of a sandcastle on the beach. The base has been packed down by the tide for hours (stiff and safe). But the very top tip of the turret, which you just built, is loose and wet. If you shake the beach, the whole castle doesn't collapse, but that loose tip might crumble and slide off.

How They Studied It: The "Digital Twin"

The team didn't go to the dam to dig holes. Instead, they built a digital twin of the sediment mountain using computer software (QUAKE/W and SLOPE/W).

The Shake: They simulated two types of earthquakes: a moderate one (0.28 g) and a severe, "worst-case" design earthquake (0.5 g).
The Stress Test: They watched how the pressure changed inside the sand. They looked for "hotspots" where the stress was highest and the water pressure was pushing the sand apart.
The Overlay: They took five different maps (showing stress, water pressure, susceptibility, movement, and sliding paths) and stacked them on top of each other like a stack of transparent sheets.

The Discovery: The "Sweet Spot" of Danger

When they looked at the stacked maps, a clear pattern emerged. All the bad signs pointed to one specific location: the toe of the front slope.

The Stress: The shaking concentrated force right at the front edge.
The Water: The water pressure was highest right under that front edge, acting like a lubricant.
The Slide: The computer calculated that if the mountain were to slide, it would slide along the weak layer right under that front edge.

The Result:

No Earthquake: The mountain is stable (Safety Score: 1.4).
Moderate Earthquake: It's still stable, but the safety margin shrinks (Safety Score: 1.2).
Severe Earthquake: The mountain is on the verge of sliding (Safety Score: 0.95). Note: A score below 1.0 means failure is imminent.

The researchers found that the mountain would likely start to slide at an earthquake intensity of about 0.45 g. This is a "yield point"—the moment the sediment starts to move permanently toward the dam.

Why This Matters: The "Choking" Risk

If that front tip of the underwater mountain slides, it won't just move a little bit. It will rush toward the dam's low-level tunnels.

The Analogy: Imagine trying to drink a milkshake through a straw, but someone pushes a glob of thick pudding right up against the tip of the straw. You can't get any liquid through. That is the risk: the sliding sediment could "choke" the tunnels, cutting off water for power and irrigation for months while workers try to clear it.

The Proposed Solutions

Since the danger is localized to that specific front tip, the researchers suggest targeted fixes rather than trying to fix the whole reservoir:

The Underwater Wall: Build a rock wall underwater in front of the sediment to hold it back, effectively removing the steep slope that wants to slide.
The Bypass: Create a way to flush the new sediment out of the reservoir before it builds up that dangerous front tip.
The Vacuum: Use dredges to suck out the loose sediment right at the dangerous front edge.

The Bottom Line

This study is a screening tool, not a final construction blueprint. It uses a "low-data" method to quickly identify where the danger is.

The main takeaway: The Tarbela Dam's underwater sediment isn't going to turn into a giant soup and collapse entirely. However, the loose, fresh front edge is a ticking time bomb during a major earthquake. If it slides, it could block the dam's tunnels. The solution is to reinforce or remove that specific front edge before a big quake hits.

Technical Summary: Mechanism-Based Screening of Seismic Instability in the Tarbela Reservoir Sediment Delta

Problem Definition
The Tarbela Dam on the Indus River in Pakistan faces a persistent safety challenge driven by long-term reservoir sedimentation. Over decades, sediment accumulation has formed a prograding underwater delta that advances toward the dam's low-level tunnel intakes. This delta consists of a flat topset, a steep foreset, and a deeper bottomset. The primary concern is not merely the loss of storage capacity, but the risk of earthquake-triggered underwater slides. If the steep, loose, and saturated foreset fails during seismic shaking, the mobilized sediment could migrate toward the tunnel approaches, potentially choking power and irrigation releases for months. While the mechanics of such submarine mass movements are understood, there is a lack of mechanism-based interpretation in the literature that connects dynamic stress redistribution, pore-pressure evolution, and slope stability specifically for the Tarbela delta within a single framework to assess intake safety.

Methodology
This study employs a coupled numerical workflow using GeoStudio software to perform a mechanism-based seismic screening of a representative delta cross-section. The analysis proceeds in two stages:

Dynamic Response (QUAKE/W): A linear-elastic dynamic analysis is conducted to simulate the redistribution of effective stresses and pore-water pressures under horizontal peak ground accelerations (PGA) of 0.28 g (representing a moderate historical event) and 0.5 g (representing a design-level event). The model utilizes parameters derived from published literature for comparable loose, saturated deltaic sediments, as site-specific cyclic test data were unavailable.
Stability Assessment (SLOPE/W): The stress states generated by the dynamic analysis are imported into a limit-equilibrium stability analysis using the Morgenstern–Price method. This allows for the evaluation of the factor of safety (FoS) under non-uniform dynamic stress conditions rather than simplified pseudo-static assumptions.

The core innovation lies in an output-driven hotspot overlay framework. The authors integrate four distinct indicators—effective stress reduction, dynamic horizontal stress concentration, elevated pore-water pressure, and a relative cyclic-demand susceptibility index ( $I_{LS}$ )—with the critical slip surface and seismic deformation fields. This overlay identifies zones where multiple adverse indicators coincide, isolating the controlling failure mechanism.

Key Results

Localized Hazard Zone: The analysis consistently identifies the foreset toe and the underlying bottomset contact as the controlling zone of instability. While the bulk of the delta body has consolidated over decades, increasing its cyclic resistance, the advancing foreset toe remains loose and highly susceptible.
Subordinate Liquefaction: Contrary to the assumption of generalized liquefaction across the entire delta, the results show that liquefaction susceptibility is localized. The relative susceptibility index ( $I_{LS}$ ) reaches 2.2–2.6 at the fresh foreset toe but remains low (0.4–1.0) in the consolidated body. The governing mechanism is a stress-driven failure along the weak bottomset contact rather than a flow slide of the entire mass.
Stability Margins: The indicative minimum factor of safety decreases with increasing seismic demand:
- Static conditions: FoS ≈ 1.4
- 0.28 g shaking: FoS ≈ 1.2 (Stable, but reduced margin)
- 0.50 g shaking: FoS ≈ 0.95 (Imminent failure condition)
- The yield acceleration (where FoS = 1) is estimated via interpolation to be approximately 0.45 g.
Deformation Direction: Seismic deformation under 0.5 g loading is directed toward the Main Embankment Dam (MED), consistent with the historical advance of the delta. This movement aligns with the critical slip surface passing through the foreset-bottomset interface.
Field Corroboration: The model's prediction of sediment movement toward the dam is supported by hydrographic survey data from the 2005 earthquake, which recorded sediment shifts and deposition near the MED consistent with the modeled direction.

Significance and Claims
The paper positions itself as a screening framework rather than a final design prediction. Its primary contribution is demonstrating that the seismic hazard at Tarbela is a localized, intake-relevant pathway driven by the geometry and material properties of the foreset, rather than a basin-wide liquefaction event.

By integrating dynamic stress response with stability outcomes, the study provides a "low-data, transferable" method for identifying controlling mechanisms in sedimented reservoirs. The authors claim that this approach supports targeted mitigation strategies, specifically:

Construction of an underwater restraining dike to eliminate the steep foreset geometry.
Implementation of sediment bypass and flushing to reduce the rate of delta advance.
Targeted dredging at the intake approach and foreset toe.

The study concludes that while the absolute factors of safety are indicative and require site-specific calibration for design purposes, the identification of the foreset-toe/bottomset contact as the critical hazard zone is robust across parameter variations. This framework allows engineers to move beyond storage-loss assessments to address specific intake-safety risks in seismic regions.