Research article

Assessing the stability of the reservoir rim in moraine deposits for a mega RCC dam

  • Received: 22 February 2024 Revised: 02 April 2024 Accepted: 15 April 2024 Published: 28 April 2024
  • Diamer Basha Dam is an under-construction, 272-meter-high, roller compacted concrete (RCC) dam on the Indus River in Pakistan. Once constructed, it will be the world's highest RCC gravity dam with a 105-kilometer-long reservoir. Most of the reservoir lies in unstable moraine deposits with steep slopes. Events like saturation during reservoir filling, alternate wetting, drawdown during reservoir operation, or a seismic event could trigger a large mass movement of these slopes into the reservoir to disrupt the dam functionality. This work identified the 15 most vulnerable slide areas using digital slope maps, elevation maps, and satellite imagery. Deterministic slope stability analysis was carried out on the identified sections under various stages of reservoir operation for static and seismic loading, using pseudo-static and dynamic analysis approaches. Probabilistic analysis was then performed using Monte Carlo simulation. The findings showed that most moraine deposits would collapse under reservoir filling, rapid drawdown, or seismic activity. Following the assessments, landslide susceptibility maps were generated, and an assessment of potential impacts, including the generation of dynamic waves, reservoir blockage, increased sediment loads, and reduced reservoir storage capacity, was also performed.

    Citation: Khalid Ahmad, Umair Ali, Khalid Farooq, Syed Kamran Hussain Shah, Muhammad Umar. Assessing the stability of the reservoir rim in moraine deposits for a mega RCC dam[J]. AIMS Geosciences, 2024, 10(2): 290-311. doi: 10.3934/geosci.2024017

    Related Papers:

    [1] Mengze Wu, Dejun Xie . The impact of ESG performance on the credit risk of listed companies in Shanghai and Shenzhen stock exchanges. Green Finance, 2024, 6(2): 199-218. doi: 10.3934/GF.2024008
    [2] Samuel Asante Gyamerah, Clement Asare . A critical review of the impact of uncertainties on green bonds. Green Finance, 2024, 6(1): 78-91. doi: 10.3934/GF.2024004
    [3] Aikaterini (Katerina) Tsoukala, Georgios Tsiotas . Assessing green bond risk: an empirical investigation. Green Finance, 2021, 3(2): 222-252. doi: 10.3934/GF.2021012
    [4] Rongwu Zhang, Huaqian Chen, Wenjia Zhang, Tong Lu . Does ESG performance inhibit or promote herding behavior of institutional investors?. Green Finance, 2025, 7(2): 363-380. doi: 10.3934/GF.2025013
    [5] Thabo J. Gopane . Policy implications of ESG-moderated credit risk on bank profitability. Green Finance, 2025, 7(3): 406-428. doi: 10.3934/GF.2025015
    [6] Raja Elyn Maryam Raja Ezuma, Nitanan Koshy Matthew . The perspectives of stakeholders on the effectiveness of green financing schemes in Malaysia. Green Finance, 2022, 4(4): 450-473. doi: 10.3934/GF.2022022
    [7] Amr ElAlfy, Adel Elgharbawy, Tia Rebecca Driver, Abdul-Jalil Ibrahim . Sustainability disclosure in the Gulf Cooperation Council (GCC) countries: Opportunities and Challenges. Green Finance, 2025, 7(1): 40-82. doi: 10.3934/GF.2025003
    [8] Zach Williams, Heather Apollonio . The causation dilemma in ESG research. Green Finance, 2024, 6(2): 265-286. doi: 10.3934/GF.2024011
    [9] María Cantero-Saiz, Sergio Sanfilippo-Azofra, Begoña Torre-Olmo, Violeta Bringas-Fernández . ESG and bank profitability: the moderating role of country sustainability in developing and developed economies. Green Finance, 2025, 7(2): 288-331. doi: 10.3934/GF.2025011
    [10] Tao Lin, Mingyue Du, Siyu Ren . How do green bonds affect green technology innovation? Firm evidence from China. Green Finance, 2022, 4(4): 492-511. doi: 10.3934/GF.2022024
  • Diamer Basha Dam is an under-construction, 272-meter-high, roller compacted concrete (RCC) dam on the Indus River in Pakistan. Once constructed, it will be the world's highest RCC gravity dam with a 105-kilometer-long reservoir. Most of the reservoir lies in unstable moraine deposits with steep slopes. Events like saturation during reservoir filling, alternate wetting, drawdown during reservoir operation, or a seismic event could trigger a large mass movement of these slopes into the reservoir to disrupt the dam functionality. This work identified the 15 most vulnerable slide areas using digital slope maps, elevation maps, and satellite imagery. Deterministic slope stability analysis was carried out on the identified sections under various stages of reservoir operation for static and seismic loading, using pseudo-static and dynamic analysis approaches. Probabilistic analysis was then performed using Monte Carlo simulation. The findings showed that most moraine deposits would collapse under reservoir filling, rapid drawdown, or seismic activity. Following the assessments, landslide susceptibility maps were generated, and an assessment of potential impacts, including the generation of dynamic waves, reservoir blockage, increased sediment loads, and reduced reservoir storage capacity, was also performed.



    1. Introduction

    Relatively large areas in the western region of the United States are classified as arid or semi-arid environments, which are characterized in part by their limited and variable precipitation. Semi-arid regions are expected to receive around 10 to 30 inches of average annual precipitation ( 254-762 mm year1) but climate change is predicted to increase the intensity and frequency of droughts globally [24]. For instance, in August of 2016, abnormally dry to moderate drought conditions were observed in several locations of the western region of the United States, with ranges between severe to extreme drought occurring in the northern portions, and severe to exceptional drought extending from California into Nevada [20]. The lack of precipitation puts indigenous species of plants and animals in semi-arid environments under unusual stress and the parallel habitat loss might pose a threat to local biodiversity [16]. Under these circumstances, having educated guesses of the potential vegetation biomass responses in semi-arid landscapes to long term changes in precipitation could serve to put forward the design of adaptation and conservation policies [9]. The estimation of the expected time of transition to a desert state (or bare-soil), as a conceivable measure of those responses, presents difficulties due to the complexities associated with specific water-vegetation systems. For instance, the inherent non-linear characteristics of semi-arid landscapes may trigger desertification in response to slow changes [23]. Results from simulations of simple mathematical models could still offer a hint of the relationship with the parameters that might be driving decline.

    Extensive mathematical modeling and analysis of semi-arid water-vegetation systems emerged for at least the past twenty years, especially since the appearance of deterministic ecohydrological models with broad developments focused on vegetation pattern formation, see for instance [5,12,19,22,25,26,27]. The effects of noise on dryland ecosystems that are usually described by deterministic models showing bistability have also been analysed, [6,21], suggesting the possibility of creation or disappearance of vegetated states in the form of noise-induced transitions, [11]. In contrast, the goals in this paper are to present a stochastic differential equation approximation for an idealized water-vegetation (non-spatial) discrete system and the estimation of mean transition times into a desert state. We emphasize that our model is not obtained by adding "noise" to a differential equation as previously done [6,21], but by the construction of a diffusion approximation.

    As a first step, we set up a Markov jump process that incorporates the interactions in an idealized water-vegetation system. Similar conceptual models have been used successfully in other biological contexts, see for instance [17] and [18,15]. The model involves only water and vegetation biomass, in an environment of limited capacity. When this capacity (or "system size") increases it gives rise to a deterministic system of differential equations for the mean densities. We deduce an intermediate mesoscale stochastic model between the jump process and the differential equations obtained for the means. Using estimated parameters for vegetation and precipitation in semi-arid landscapes from the literature, and data for state precipitation anomalies in California as baseline, we estimate the mean times for a system to reach desertification in a range of realistic precipitation anomalies, i.e. departures from long term mean. With these results we finally quantify, for this simple model, the dependence between changes in precipitation anomalies and mean transition times to the desert state.


    2. Theoretical framework for water-vegetation dynamics


    2.1. The stochastic and ODE models

    We start by defining a Markov jump process that represents a simplified version of the real interactions between water and vegetation at a small scale. We do this through the discretization of (alive) vegetation biomass and water volume in small units (individuals), for which a specific set of stochastic events can be explicitly characterized. Naturally, as we transition into larger scales discreteness is lost, and the continuous state space takes place.

    We first consider a patch with finite capacity, say N, containing three types of individuals: (ⅰ) vegetation biomass units; (ⅱ) water volume units; (ⅲ) empty locations. The dynamics of plant biomass and water interactions is driven by events involving a few processes: vegetation biomass loss, incoming water, water evaporation, and increase vegetation yield by taking up water. If we represent the state of the system, i.e. the number of biomass units, n, and the number of water volume units, m, by the pair (n,m) then these events correspond to the transitions detailed in Table 1.

    Table 1. Possible transition events with their associated jumps if the system is at state (n,m), where n and m represent units of biomass and water, respectively.
    Event Transition Jump Jump rate
    Vegetation biomass loss (n,m)(n1,m) (1,0) d
    Incoming water (n,m)(n,m+1) (0,1) s
    Water evaporation (n,m)(n,m1) (0,1) v
    Increase vegetation by (n,m)(n+1,m1) (1,1) b
    water take up
     | Show Table
    DownLoad: CSV

    It is then straightforward to find the probability rates of transition from a state (n,m) to (n,m), T(n,m|n,m):

    (a) T(n+1,m1|n,m)=bnNmN1

    (b) T(n1,m|n,m)=dnN

    (c) T(n,m+1|n,m)=sNnmN

    (d) T(n,m1|n,m)=vmN

    Using these rates we can write the associated Kolmogorov equation (see [7] for instance),

    dP(n,m,t)dt=T(n,m|n1,m+1)P(n1,m+1,t)+T(n,m|n+1,m)P(n+1,m,t)+T(n,m|n,m+1)P(n,m+1,t)+T(n,m|n,m1)P(n,m1,t)(T(n+1,m1|n,m)+T(n1,m|n,m)+T(n,m1|n,m)+T(n,m+1|n,m))P(n,m,t).

    where P(n,m,t) is the probability that the system is at the state (n,m) at time t. Imposing zero boundary conditions, multiplying the Kolmogorov equation by n, and summing over m and n gives the expression for the rate of change of the mean biomass,

    dndt=Nn,m=0[T(n+1,m1|n,m)T(n1,m|n,m)]P(n,m,t)=bnNmN1dnN, (1)

    where the correlations between the random variables are neglected under the assumption of a large N. Writing the mean density of vegetation as ρv=n/N finally gives

    dρvdt=˜bρvρw˜dρv, (2)

    where ˜b=b/(N1) and ˜d=d/N. Similarly,

    dmdt=Nn,m=0[T(n,m+1|n,m)T(n,m1|n,m)]P(n,m,t)        =ssn+mNbnNmN1vmN (3)

    Similarly, writing the mean density of water by ρw=m/N then

    dρwdt=˜s(1(ρv+ρw))˜bρvρw˜vρw, (4)

    where ˜s=s/N and ˜v=v/N. We remark at this point that water infiltration in the soil is generally improved by the presence of vegetation. As a consequence, the process of taking up water by plants gets more efficient. This fact can be incorporated into the model by letting the rate for taking up water be density dependent, that is, to the jump (1, -1) we associate a new rate bn/N. This leads to having T(n+1,m1|n,m)=bn2N2mN1, which replaced into the equation for dn/dt produces

    dρvdt=˜bρ2vρw˜dρv (5)

    instead of equation (2). A similar change happens in equation (4), which is now

    dρwdt=˜s(1(ρv+ρw))˜bρ2vρw˜vρw. (6)

    Equations (5) and (6) constitute a system of differential equations that serves as approximation to the dynamics of the mean densities for large values of N, also known as mean field equations in the Physics literature. We remark that for N relatively large the factor 1(ρv+ρw) is close to one, making our equations an approximation to the reaction part of Klausmeier's reaction-diffusion-advection equations for water-vegetation systems, [12].


    2.2. The diffusion approximation

    The diffusion approximation to our model (the mesoscale model) describes the system as an intermediate approximation that emerges between the Markov jump process model (microscale model) and the differential equations for the mean densities (macroscale model). For this approximation the state variables are continuous but include random fluctuations. We expect the new model to incorporate the differential equations and terms that express random fluctuations around the mean densities.

    To obtain a representation of the Markov jump process as a diffusion process one can follow either Kurtz's method [14], or find the same equations via the Fokker-Planck equation [18,7]; see [4] for a nice introduction. In Kurtz's approach, which is the one we use here, the jump process is represented by

    X(t)=X(0)+rrN(r)(t0NΦ(1NX(s);r)ds), (7)

    where r is a jump (see the third column in Table 1), N(r)(t) is a collection of independent rate 1 scalar Poisson processes and NΦ(1NX(s);r) is the rate at which a transition with jump r occurs. Kurtz showed that the process X(t)/N can be approximated by a process Y(t) that satisfies the stochastic differential equation

    dY(t)=rrΦ(Y(t);r)dt+1Nrr(Φ(Y(t);r)dW(r)(t), (8)

    where the W(r) are the Brownian motions associated with the jump types of the Markov jump process. The error introduced on bounded intervals of time by replacing X(t)/N with Y(t) is O(logN/N) as N. The first term in the sum of the right hand side of (8) corresponds to the vector field of mean densities, that is,

    A=[AvAw]=[˜bρ2vρw˜dρv,˜s(1ρvρw)˜bρ2vρw˜vρw,].

    If we denote by B the covariance matrix

    B=[BvvBvwBwvBww],

    where

    Bvv=˜bρ2vρw+˜dρv,Bvw=Bwv=˜bρ2vρw,Bww=˜bρ2vρw+˜s(1ρvρw)+˜vρw,

    and factor it as B=ggT for some g, then the stochastic differential equations system (8) can be written as

    dY(t)=Adt+1NgdW, (9)

    where dW is a two dimensional Brownian motion. Figure 1 compares several paths that correspond to the Markov jump process and the diffusion approximation, generated with the same set of (arbitrarily chosen) parameters. It is apparent that the stochastic differential equation provides a good approximation to the jump process.

    Figure 1. Simulations corresponding to the Markov jump process (left) and the diffusion approximation (right). For comparison purposes the paths in both panels were generated using the same parameters and the same scaled time.

    3. Simulations

    We use the stochastic differential equation (9) to simulate the water-vegetation system and obtain averages of the expected time to desertification (see [10] or [13] for a quick or an extensive introduction respectively to the numerical solution of stochastic differential equations).

    By identifying the parameters of the nondimensional deterministic (non-spatial) model in [12] with the mean field system obtained above (Av=Av(ρ) and Aw=Aw(ρ), where ρ=(ρv,ρw)) we obtain S=AR1/2J/L3/2 and ˜d=M/L, with the meaning and realistic values for these parameters listed in Table 3. Thus, for instance, Klausmeier's parameter ranges for S are [0.077,0.23] and ˜d=0.045 for trees, and [0.94,2.81] and ˜d=0.45 for grass. Also, the corresponding value for ˜b and ˜v in both cases is one. Regarding the average evaporation rate, we follow Klausmeier assuming that the equilibrium of water (in his deterministic model) is at w=75 mm, and then computing the associated evaporation rate given the averaged annual precipitation, [12]. For example, with A=300 mm year1 the evaporation rate is A/w=4 year1.

    Table 2. Parameters for semi-arid landscapes, taken from [12].
    Parameter Definition Estimated values Units
    R uptake rate of water 1.5(trees) - 100(grass) mm year1 (kg dry mass)2
    J yield of plant biomass 0.002(trees) - 0.003(grass) kg dry mass (mm)1
    M mortality rate 0.18(trees) - 1.8(grass) year1
    A precipitation 250 - 750 mm year1
    L evaporation rate 4 year1
     | Show Table
    DownLoad: CSV

    Precipitation anomalies records, i.e. records of the deviations from a long term precipitation mean, have a negative trend in specific geographic drought events. In the state of California for instance, which experienced unusually long drought conditions, the (state) average of the precipitation anomalies for the past 16 years is 2.07 (inches year1) (i.e. -52.58 mm year1), see Figure 2. Although there is a lot of variability across the state, we use this value for illustration purposes, and plot the data as if this negative deviation from the long term precipitation average were steady in time.

    Figure 2. State averages of precipitation anomalies for 2000-2016 in California (inches year1). The averaged anomaly (difference from long term average) during that period is -2.07 (inches year1) (-52.58 mm year1). The precipitation increase expected from El Niño for the winter 2015-2016 was scarcely above the long term state average. Data/image provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, from their Web site at http://www.esrl.noaa.gov/psd/.

    The results of the simulations are shown in Figure 3, portraying a roughly linear relationship between the time to desertification and the anomalies in precipitation in the range selected. The simulations were run using parameters for trees (top panel) and grass (low panel), with system capacity N=500. If we denote by T the average time to desertification and by P the average annual precipitation we can define the (dimensionless) sensitivity index S0 as the ratio S0=(ΔT/T)/(ΔP/P), where ΔT and ΔP represent the absolute change in the variables T and P, i.e. S0 measures the relative change in T with respect to the relative change in P, see for instance [2]. Direct computation from the averaged results gives S02. Similar results were obtained with the larger system capacity N=104.

    Figure 3. Examples of how average time to desertification might be affected by a reduction in average annual precipitation. Parameters for trees were used in panel (a) and for grass in panel (b). The average of negative anomalies similar to that observed for the last years in California is around 50 mm year1. The model suggests that the sensitivity index S02, i.e. relative changes in the mean time to desertification are roughly twice the relative changes in average annual precipitation. For the simulations, N=500 and the initial conditions were ρ(0)=(0.1,0.1) (squares), ρ(0)=(0.5,0.5) (stars) and ρ(0)=(0.9,0.1) (triangles). Each time average was obtained from 50000 simulations. Panel (c) shows the histograms corresponding to the simulated times to desertification with an average annual precipitation of 200 and 250 mm year1 (for grass) on the left and right, respectively. The simulations used the same initial conditions ρ(0)=(0.1,0.1).

    4. Conclusions and discussion

    A traditional approach for modeling interacting populations at the macroscopic level assumes that the terms in the equations that drive the dynamics represent the average effects of individual interactions in a general, all-inclusive way. Subsequent developments use those models as departing points for building theoretical extensions by incorporating further complexity to the equations, like the inclusion of spatial dependence by adding diffusion, or the introduction of "noise" terms. A different modeling approach is to start at the individual level, with explicit rules describing the interactions between individuals and their environment, [17]. This alternative implies the definition of a Markov jump process that constitutes the foundation for developing definitive model approximations that relate macro and microscopic dynamical levels. In this paper we have taken the latter approach for constructing a stochastic differential equation (continuous state space) that approximates the dynamics of an idealized water-vegetation system, initially conceived as a Markov jump process (with discrete states as a proxy for small scale). Thus, our work complements the existing literature on modeling noise in drylands, [21].

    The diffusion approximation obtained, together with parameter values for vegetation and precipitation for semi-arid landscapes extracted from the literature, and data on decreasing precipitation trends in California, were used to estimate average times for desertification. For a fixed system capacity (N=500) the simulations for trees and grass suggests that the sensitivity of the time to desertification from the annual precipitation is roughly similar, and approximately equal to 2, i.e. the relative change in the transition time into a desert state is equal to twice the relative change in precipitation. Repeating the simulations, for different initial conditions and for a larger capacity (N=104), provided the same approximate numerical relation. The simulations (see Figure 3) suggest, for instance, that a decrease of roughly 0.4 inches of precipitation (10 mm year1) might reduce times to desertification in some cases by more than 25 years for the case of trees, and around 5 years in the case of grass. Put another way, current trends of desertification could be significantly boosted if the patterns of increasing precipitation anomalies are maintained. However, looking at the basic transition mechanisms considered in formulating the Markov jump process, it is clear that the model should be used with care to draw any conclusions on specific vegetation-water systems.

    We remark that extended droughts may resemble desertification, but the return of seasonal precipitation events may recover the vegetation (see for instance [1] where desertification was limited to spatially localized areas). This suggests that the inclusion of patterns of precipitation anomalies restricted to relatively small areas would provide more reliable results. For systems with relatively small capacity we notice that the times to desertification may be reduced dramatically (see Figure 4).

    Figure 4. Left: time to desertification for A=250 (dashes) and A=200 (dot-dashes) as function of the system capacity N. The sensitivity of the time to desertification from the annual precipitation was computed for N=10000 showing to be the same as when N=500, i.e. 2. As N increases both times to desertification also increase, but they get reduced dramatically as N gets smaller. Right: Difference between the curves in the contiguous plot. Although the difference increases, the sensitivity of the time to desertification from the annual precipitation is apparently similar in relatively large systems.

    Further work should also include long term variations of other climate related parameters. For instance, it has been documented that higher temperatures increase evapotranspiration rates [3], which have been observed over most of the United States, [8]. Another aspect that deserves attention is the inclusion of changes in the vegetation dynamics during dry periods, where vegetation mortality could be exacerbated. As is clearly pointed out in [28], neglecting the effects of intermittent precipitation on vegetation dynamics may influence the results considerably. Finally, the understanding of desertification will demand treatment with insightful stochastic space-time models.


    Acknowledgments

    The authors are grateful to C. Kribs and three anonymous reviewers for comments and suggestions that led to significant improvement of the paper.




    [1] Wang L, Wu CZ, Tang LB, et al. (2020) Efficient reliability analysis of earth dam slope stability using extreme gradient boosting method. Acta Geotech 15: 3135-3150. https://doi.org/10.1007/s11440-020-00962-4 doi: 10.1007/s11440-020-00962-4
    [2] Genevois R, Ghirotti M (2005) The 1963 Vaiont Landslide. G di Geol Appl 1: 41-52.
    [3] Babar MS, Israr J, Zhang G, et al. (2022) Analysis of rock slope stability for Mohmand Dam—A comparative study. Cogent Eng 9: 2124939. https://doi.org/10.1080/23311916.2022.2124939 doi: 10.1080/23311916.2022.2124939
    [4] Mushtaq F, Rehman H, Ali U, et al. (2023) An Investigation of Recharging Groundwater Levels through River Ponding: New Strategy for Water Management in Sutlej River. Sustainability 15: 1047. https://doi.org/10.3390/SU15021047 doi: 10.3390/SU15021047
    [5] Basharat M, Ali SU, Azhar AH (2014) Spatial variation in irrigation demand and supply across canal commands in Punjab: a real integrated water resources management challenge. Water Policy 16: 397-421. https://doi.org/10.2166/WP.2013.060 doi: 10.2166/WP.2013.060
    [6] Ali U, Otsubo M, Ebizuka H, et al. (2020) Particle-scale insight into soil arching under trapdoor condition. Soils Found 60: 1171-1188. https://doi.org/10.1016/j.sandf.2020.06.011 doi: 10.1016/j.sandf.2020.06.011
    [7] Otsubo M, Kuwano R, Ali U, et al. (2018) Trapdoor model test and DEM simulation associated with arching. Physical Modelling in Geotechnics, CRC Press, 233-238.
    [8] Ali Naqvi U, Otsubo M (2020) A study on arching mechanism in trapdoor model test and equivalent discrete element simulations. 16th Asian Regional Conference on Soil Mechanics and Geotechnical Engineering.
    [9] Wang L, Wu C, Gu X, et al. (2020) Probabilistic stability analysis of earth dam slope under transient seepage using multivariate adaptive regression splines. Bull Eng Geol Environ 79: 2763-2775. https://doi.org/10.1007/s10064-020-01730-0 doi: 10.1007/s10064-020-01730-0
    [10] Santosh S, Khazanchi RN, Mishra A, et al. (2017) Reservoir Rim Stability: A Case Study of Punatsangchhu-I Hydroelectric Project Bhutan. INDOROCK 2017: 7th Indian Rock Conference, New Delhi, India.
    [11] MWH Consultants (2013) Preliminary Reservoir Slope Stability Assessment for Susitna-Watana Dam. Available from: https://www.arlis.org/docs/vol1/Susitna2/2/SuWa289/SuWa289sec4-5att3.pdf.
    [12] Wang Y, Liu X, Zhang X, et al. (2023) Dynamic response characteristics of moraine-soil slopes under the combined action of earthquakes and cryogenic freezing. Cold Reg Sci Technol 211: 103854. https://doi.org/10.1016/j.coldregions.2023.103854 doi: 10.1016/j.coldregions.2023.103854
    [13] Klimeš J, Novotný J, Novotná I, et al. (2016) Landslides in moraines as triggers of glacial lake outburst floods: example from Palcacocha Lake (Cordillera Blanca, Peru). Landslides 13: 1461-1477. https://doi.org/10.1007/s10346-016-0724-4 doi: 10.1007/s10346-016-0724-4
    [14] Ghimire A, Gajurel A (2020) Moraine dam stability analysis of Ngozumpa glacier in Gokyo area, Nepal. American Geophysical Union, Fall Meetin. Available from: https://www.agu.org/fall-meeting-2020.
    [15] Anbalagan A, Kumar A (2015) Reservoir induced landslides-a case study of reservoir rimregion of Tehri Dam. TIFAC-IDRiM Conference, New Delhi, India.
    [16] Chen XP, Huang JW (2011) Stability analysis of bank slope under conditions of reservoir impounding and rapid drawdown. J Rock Mech Geotech Eng 3: 429-437. https://doi.org/10.3724/SP.J.1235.2011.00429 doi: 10.3724/SP.J.1235.2011.00429
    [17] Yin YP, Huang BL, Wang WP, et al. (2016) Reservoir-induced landslides and risk control in Three Gorges Project on Yangtze River, China. J Rock Mech Geotech Eng 8: 577-595. https://doi.org/10.1016/J.JRMGE.2016.08.001 doi: 10.1016/J.JRMGE.2016.08.001
    [18] Chuaiwate P, Jaritngam S, Panedpojaman P, et al. (2022) Probabilistic Analysis of Slope against Uncertain Soil Parameters. Sustainability 14: 14530. https://doi.org/10.3390/su142114530 doi: 10.3390/su142114530
    [19] Kumar S, Choudhary SS, Burman A, et al. (2023) Probabilistic Slope Stability Analysis of Mount St. Helens Using Scoops3D and a Hybrid Intelligence Paradigm. Mathematics 11: 3809. https://doi.org/10.3390/math11183809 doi: 10.3390/math11183809
    [20] Petterson MG (2010) A Review of the geology and tectonics of the Kohistan island arc, north Pakistan. Geol Soc London Spec Publ 338: 287-327. https://doi.org/10.1144/SP338.14 doi: 10.1144/SP338.14
    [21] Khan SD, Walker DJ, Hall SA, et al. (2009) Did the Kohistan-Ladakh island arc collide first with India? GSA Bull 121: 366-384. https://doi.org/10.1130/B26348.1 doi: 10.1130/B26348.1
    [22] Ewing TA, Müntener O (2018) The mantle source of island arc magmatism during early subduction: Evidence from Hf isotopes in rutile from the Jijal Complex (Kohistan arc, Pakistan). Lithos 308-309: 262-277. https://doi.org/10.1016/J.LITHOS.2018.03.005 doi: 10.1016/J.LITHOS.2018.03.005
    [23] Bignold SM, Treloar PJ, Petford N (2006) Changing sources of magma generation beneath intra-oceanic island arcs: An insight from the juvenile Kohistan island arc, Pakistan Himalaya. Chem Geol 233: 46-74. https://doi.org/10.1016/J.CHEMGEO.2006.02.008 doi: 10.1016/J.CHEMGEO.2006.02.008
    [24] Jagoutz O, Müntener O, Burg JP, et al. (2006) Lower continental crust formation through focused flow in km-scale melt conduits: The zoned ultramafic bodies of the Chilas Complex in the Kohistan island arc (NW Pakistan). Earth Planet Sci Lett 242: 320-342. https://doi.org/10.1016/J.EPSL.2005.12.005 doi: 10.1016/J.EPSL.2005.12.005
    [25] Dhuime B, Bosch D, Garrido CJ, et al. (2009) Geochemical Architecture of the Lower- to Middle-crustal Section of a Paleo-island Arc (Kohistan Complex, Jijal-Kamila Area, Northern Pakistan): Implications for the Evolution of an Oceanic Subduction Zone. J Petrol 50: 531-569. https://doi.org/10.1093/PETROLOGY/EGP010 doi: 10.1093/PETROLOGY/EGP010
    [26] Iturrizaga L (2001) Lateroglacial valleys and landforms in the Karakoram Mountains (Pakistan). GeoJournal 54: 397-428. https://doi.org/10.1023/A:1021365416056 doi: 10.1023/A:1021365416056
    [27] GEOSLOPE (2012) User Manual, GEO-SLOPE. GEO-SLOPE International Ltd Calgary, Canada.
    [28] Lebourg T, Riss J, Pirard E (2004) Influence of morphological characteristics of heterogeneous moraine formations on their mechanical behaviour using image and statistical analysis. Eng Geol 73: 37-50. https://doi.org/10.1016/J.ENGGEO.2003.11.004 doi: 10.1016/J.ENGGEO.2003.11.004
    [29] Wang BL, Yao LK, Zhao HX, et al. (2018) The maximum height and attenuation of impulse waves generated by subaerial landslides. Shock Vib 2018.https://doi.org/10.1155/2018/1456579
  • geosci-10-02-017-s001.pdf
  • This article has been cited by:

    1. E. M. Starodubtseva, I. N. Tsymbalov, D. A. Gorlova, K. A. Ivanov, A. B. Savel’ev, Polarization state control for high-peak-power applications, 2025, 42, 1084-7529, 285, 10.1364/JOSAA.544487
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1822) PDF downloads(91) Cited by(0)

Article outline

Figures and Tables

Figures(16)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog