skip to main content


Title: Suyya’s Flood: Numerical Models of Kashmir’s Medieval Megaflood and Ancient Lake Kerewa Drainage Events
In the mid-ninth century, an earthquake triggered a landslide that blocked the narrow gorge of the Jhelum River where it exits the Kashmir Valley. The landslide impounded a lake that extended ≈100 km along the floor of the valley, implying an impounded volume of ≤21 km 3 , flooding the capital, Srinagar, and much agricultural land. An engineered breach of the landslide was contrived by a Medieval engineer resulting in the catastrophic release of flood waters. Using reasonable assumptions we calculate the probable minimum drainage time of this Medieval flood (<4 days) and maximum downstream surge velocities (≈12 m/s). These would have been sufficient to transport boulders in the bed of the Jhelum with dimensions of ≈6 m, consistent with those currently present in some reaches of the river. Given the morphology of the Jhelum gorge we consider that landslide outburst floods may have been common in Kashmir’s history. Ancient shorelines indicate that paleo-lake volumes in the Kashmir Valley may have exceeded 400 km 3 which, were they released in catastrophic floods, would have been associated with potential downstream outburst velocities >32 m/s, able to transport boulders with dimensions ≈40 m, far in excess of any found in the course of the Jhelum or in the Punjab plains. Their absence suggests that Kashmir’s ancient lakes were not lowered by outburst mechanisms much exceeding those associated with Suyya’s flood. Present-day floods have been many tens of meters shallower than those impounded by landslides in the Jhelum in the past several thousands of years. A challenge for future study will be to date Kashmir’s ancient shorelines to learn how often landslides and major impoundment events may have occurred in the valley.  more » « less
Award ID(s):
2028557
NSF-PAR ID:
10323672
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Earth Science, Systems and Society
Volume:
1
ISSN:
2634-730X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Profound effects of episodic megafloods (≥106 m3/s) have been observed on Earth and Mars. Quaternary megafloods sourced from valley‐blocking glaciers on the Tibetan Plateau likely play an important role in the geomorphic evolution of the Yarlung‐Tsangpo Gorge and mountain landscape of the eastern Himalaya. We use the first 2D numerical simulation of a megaflood sourced from a reconstructed 81 km3Tibetan lake to analyze flood hydraulics and examine the erosional and depositional potential of megafloods in mountain landscapes. The simulated flood has a duration >60 hr and a peak discharge of 3.1 × 106 m3/s. We find that the extent of inundated features like terraces, narrow valley sections, tight meander bends, and overtopped ridges influences locations of observed maximum depth (370 m), speed (76 m/s), and bed shear stress (>100 kPa), creating dynamic patterns of erosive potential. Consequently, it is difficult to predict local (≤1 km) patterns of megaflood erosional potential from either unit stream power or flood power from smaller magnitude outburst floods. However, both are useful when predicting regional (≥25 km) order‐of‐magnitude shifts in megaflood flood power. Portions of the flood domain downstream of the Gorge experience lower bed shear stresses and flood power <5 kW/m2, indicating potential for significant deposition. We suggest widespread deposition of boulders within the modern channel and fine‐grained particles on hillslopes during a megaflood likely impedes subsequent erosion and affects channel width and longitudinal form throughout the flood pathway. Our findings show the legacy of megaflooding in mountainous terrain includes both extensive erosion and deposition.

     
    more » « less
  2. Abstract

    Bedrock landsliding, including the formation of landslide dams, is a predominant geomorphic process in steep landscapes. Clarifying the importance of hydrologic and seismic mechanisms for triggering deep‐seated landslides remains an ongoing effort, and formulation of geomorphic metrics that predict dam preservation is crucial for quantifying secondary landslide hazards. Here, we identify >200 landslide‐dammed lakes in western Oregon and utilize dendrochronology and enhanced14C dating (“wiggle matching”) of “ghost forests” to establish slope failure timing at 20 sites. Our dated landslide dataset reveals bedrock landsliding has been common since the last Cascadia Subduction Zone earthquake in January 1700 AD. Our study does not reveal landslides that date to 1700 AD. Rather, we observe temporal clustering ofat leastfour landslides in the winter of 1889/1890 AD, coincident with a series of atmospheric rivers that generated one of the largest regionally recorded floods. We use topographic and field analyses to assess the relation between dam preservation and topographic characteristics of the impounded valleys. In contrast to previous studies, we do not observe systematic scaling between dam size and upstream drainage area, though dam stability indices for our sites correspond with “stable” dams elsewhere. Notably, we observe that dams are preferentially preserved at drainage areas of ∼1.5 to 13 km2and valley widths of ∼25 to 80 m, which may reflect the reduced downstream influence of debris flows and the accumulation of mature conifer trees upstream from landslide‐dammed lake outlets. We suggest that wood accumulation upstream of landslide dams tempers large stream discharges, thus inhibiting dam incision.

     
    more » « less
  3. Abstract

    In May 2012, a sediment‐laden flood along the Seti Khola (= river) caused 72 fatalities and widespread devastation for > 40 km in Pokhara, Nepal's second largest city. The flood was the terminal phase of a hazard cascade that likely began with a major rock‐slope collapse in the Annapurna Massif upstream, followed by intermittent ponding of meltwater and subsequent outburst flooding. Similar hazard cascades have been reported in other mountain belts, but peak discharges for these events have rarely been quantified. We use two hydrodynamic models to simulate the extent and geomorphic impacts of the 2012 flood and attempt to reconstruct the likely water discharge linked to even larger medieval sediment pulses. The latter are reported to have deposited several cubic kilometres of sediment in the Pokhara Valley. The process behind these sediment pulses is debated. We traced evidence of aggradation along the Seti Khola during field surveys and from RapidEye satellite images. We use two steady‐state flood models, HEC‐RAS and ANUGA, and high‐resolution topographic data, to constrain the initial flood discharge with the lowest mismatch between observed and predicted flood extents. We explore the physically plausible range of simplified flood scenarios, from meteorological (1000 m3 s−1) to cataclysmic outburst floods (600,000 m3 s−1). We find that the 2012 flood most likely had a peak discharge of 3700 m3 s−1in the upper Seti Khola and attenuated to 500 m3s−1when arriving in Pokhara city. Simulations of larger outburst floods produce extensive backwater effects in tributary valleys that match with the locations of upstream‐dipping medieval‐age slackwater sediments in several tributaries of the Seti Khola. Our findings are consistent with the notion that the medieval sediment pulses were linked to outburst floods with peak discharges of >50,000 m3 s−1, though discharge may have been an order of magnitude higher.

     
    more » « less
  4. Abstract

    Catastrophic drainage of glacial Lake Missoula through the Columbia River Gorge, USA, produced some of the largest floods ever known. However, erosion of the gorge during flooding has not been quantified, hindering discharge reconstructions and our understanding of landscape change by megafloods. Using a neural network and geomorphic observations, we reconstructed the gorge topography and found ∼7.4 km3of rock was eroded from gorge walls. Accounting for a narrower canyon and matching flood high‐water marks resulted in peak‐flood discharge reconstructions of 6 × 106–7 × 106 m3 s−1, which are 30%–40% lower than prior estimates based on the present‐day topography. Sediment transport modeling indicated that more frequent intermediate‐sized floods transported most of the eroded rock. Thus, similar to alluvial rivers, discharge magnitude‐frequency tradeoffs may also govern canyon formation by repeated megafloods.

     
    more » « less
  5. Abstract. To explore the sensitivity of rivers to blocking from landslidedebris, we exploit two similar geomorphic settings in California'sFranciscan mélange where slow-moving landslides, often referred to asearthflows, impinge on river channels with drainage areas that differ by afactor of 30. Analysis of valley widths and river long profiles over∼19 km of Alameda Creek (185 km2 drainage area) andArroyo Hondo (200 km2 drainage area) in central California shows avery consistent picture in which earthflows that intersect these channelsforce tens of meters of gravel aggradation for kilometers upstream, leadingto apparently long-lived sediment storage and channel burial at these sites.In contrast, over a ∼30 km section of the Eel River (5547 km2 drainage area), there are no knickpoints or aggradation upstreamof locations where earthflows impinge on its channel. Hydraulic andhydrologic data from United States Geological Survey (USGS) gages on Arroyo Hondo and the Eel River, combinedwith measured size distributions of boulders input by landslides for bothlocations, suggest that landslide derived boulders are not mobile at eithersite during the largest floods (>2-year recurrence) with field-measured flow depths. We therefore argue that boulder transport capacity isan unlikely explanation for the observed difference in sensitivity tolandslide inputs. At the same time, we find that earthflow fluxes per unitchannel width are nearly identical for Oak Ridge earthflow on Arroyo Hondo,where evidence for blocking is clear, and for the Boulder Creek earthflow onthe Eel River, where evidence for blocking is absent. These observationssuggest that boulder supply is also an unlikely explanation for the observedmorphological differences along the two rivers. Instead, we argue that thedramatically different sensitivity of the two locations to landslideblocking is related to differences in channel width relative to typicalseasonal displacements of earthflows. A synthesis of seasonal earthflowdisplacements in the Franciscan mélange shows that the channel width ofthe Eel River is ∼5 times larger than the largest annualseasonal displacement. In contrast, during wet winters, earthflows arecapable of crossing the entire channel width of Arroyo Hondo and AlamedaCreek. In support of this interpretation, satellite imagery shows thatimmobile earthflow-derived boulders are generally confined to the edges ofthe channel on the Eel River. By contrast, immobile earthflow-derivedboulders jam the entire channel on Arroyo Hondo. Our results imply that lower drainage area reaches of earthflow-dominated catchments may be particularly prone to blocking. By inhibiting the upstreampropagation of base-level signals, valley-blocking earthflows may thereforepromote the formation of so-called “relict topography”. 
    more » « less