Modeling of major cliff destabilizations and subsequent lahars in the Prêcheur catchment, Martinique.

Author(s):  
Marc Peruzzetto ◽  
Clara Levy ◽  
Yannick Thiery ◽  
Gilles Grandjean ◽  
Anne Mangeney ◽  
...  

<p>The Prêcheur river is located to the West of Montagne Pelée, in the Northern part of Martinique island. For several decades it has produced numerous lahars that directly threaten the Prêcheur village, at the mouth of the river. In recent years, the most important lahars have been correlated to massive collapses of the Samperre cliff, 9 km upstream from the sea, that create a reservoir of loose material at the bottom of the cliff. In 2010, a lahar started from this reservoir, destroyed a bridge and inundated part of the Prêcheur village. A new major period of collapses of the Samperre cliff started on 2 January 2018, involving more than 4 × 10<sup>6</sup> m<sup>3</sup> of material. In the following days, intense rainfalls triggered several lahars that reached the Prêcheur village but remained confined in the river bed. Since then, lahars and collapses have continued to occur, even though their frequency has decreased with time and their intensity is smaller compared to the onset of the crisis. One single lahar overflowed the river bed on 22 February 2018 without significant impact on infrastructures.</p><p>In this study, we test different possible scenarios to model the first and most important phase of the collapse of the Samperre cliff, that occurred in early January 2018, with the shallow-water model SHALTOP. We constrain the collapse geometry with photogrammetric 3D models and LIDAR topographic surveys, acquired in 2010 and in late January 2018. We also consider an intermediate volume to take into account a possible retreat of the cliff face between 2010 and 2018. The modeled traveled distances are compared to field observations. Finally, we use geomorphological and geological observations to identify potentially unstable structures within the cliff, and model the associated collapses.</p><p>These simulations provide insights on the possible geometry (extent and depth) of the debris reservoir at the bottom of the cliff, after a major collapse episode. This is of prior importance in order to estimate the location and volume of future lahars. In order to investigate their dynamics, we model the major 2010 lahar, for which the initial debris reservoir volume is known (about 2 Mm<sup>3</sup>). We first simulate the progressive remobilization of the reservoir by water with the r.avaflow numerical code. In a second test, we impose instead a constant flow discharge upstream until the same volume has been released. We test different parameters to identify which ones have the most significant influence on the lahar travel time, from its initiation until it reaches the Prêcheur village.</p>

2018 ◽  
Vol 850 ◽  
Author(s):  
James C. McWilliams ◽  
Cigdem Akan ◽  
Yusuke Uchiyama

Coherent vortices with horizontal swirl arise spontaneously in the wave-driven nearshore surf zone. Here, a demonstration is made of the much greater robustness of coherent barotropic dipole vortices on a sloping beach in a 2D shallow-water model compared with fully 3D models either without or with stable density stratification. The explanation is that active vortex tilting and stretching or instability in 3D disrupt an initially barotropic dipole vortex. Without stratification in 3D, the vorticity retains a dipole envelope structure but is internally fragmented. With stratification in 3D, the disrupted vortex reforms as a coherent but weaker surface-intensified baroclinic dipole vortex. An implication is that barotropic or depth-integrated dynamical models of the wave-driven surf zone misrepresent an important aspect of surf-eddy behaviour.


Water ◽  
2021 ◽  
Vol 13 (10) ◽  
pp. 1435
Author(s):  
Peng Hu ◽  
Junyu Tao ◽  
Aofei Ji ◽  
Wei Li ◽  
Zhiguo He

In this paper, a computationally efficient shallow water model is developed for sediment transport in the Yangtze estuary by considering mixed cohesive and non-cohesive sediment transport. It is firstly shown that the model is capable of reproducing tidal-hydrodynamics in the estuarine region. Secondly, it is demonstrated that the observed temporal variation of suspended sediment concentration (SSC) for mixed cohesive and non-cohesive sediments can be well-captured by the model with calibrated parameters (i.e., critical shear stresses for erosion/deposition, erosion coefficient). Numerical comparative studies indicate that: (1) consideration of multiple sediment fraction (both cohesive and non-cohesive sediments) is important for accurate modeling of SSC in the Yangtze Estuary; (2) the critical shear stress and the erosion coefficient is shown to be site-dependent, for which intensive calibration may be required; and (3) the Deepwater Navigation Channel (DNC) project may lead to enhanced current velocity and thus reduced sediment deposition in the North Passage of the Yangtze Estuary. Finally, the implementation of the hybrid local time step/global maximum time step (LTS/GMaTS) (using LTS to update the hydro-sediment module but using GMaTS to update the morphodynamic module) can lead to a reduction of as high as 90% in the computational cost for the Yangtze Estuary. This advantage, along with its well-demonstrated quantitative accuracy, indicates that the present model should find wide applications in estuarine regions.


2020 ◽  
Vol 32 (12) ◽  
pp. 124117
Author(s):  
M. W. Harris ◽  
F. J. Poulin ◽  
K. G. Lamb

Water ◽  
2021 ◽  
Vol 13 (16) ◽  
pp. 2152
Author(s):  
Gonzalo García-Alén ◽  
Olalla García-Fonte ◽  
Luis Cea ◽  
Luís Pena ◽  
Jerónimo Puertas

2D models based on the shallow water equations are widely used in river hydraulics. However, these models can present deficiencies in those cases in which their intrinsic hypotheses are not fulfilled. One of these cases is in the presence of weirs. In this work we present an experimental dataset including 194 experiments in nine different weirs. The experimental data are compared to the numerical results obtained with a 2D shallow water model in order to quantify the discrepancies that exist due to the non-fulfillment of the hydrostatic pressure hypotheses. The experimental dataset presented can be used for the validation of other modelling approaches.


Water ◽  
2021 ◽  
Vol 13 (15) ◽  
pp. 2054
Author(s):  
Naoki Kuroda ◽  
Katsuhide Yokoyama ◽  
Tadaharu Ishikawa

Our group has studied the spatiotemporal variation of soil and water salinity in an artificial salt marsh along the Arakawa River estuary and developed a practical model for predicting soil salinity. The salinity of the salt marsh and the water level of a nearby channel were measured once a month for 13 consecutive months. The vertical profile of the soil salinity in the salt marsh was measured once monthly over the same period. A numerical flow simulation adopting the shallow water model faithfully reproduced the salinity variation in the salt marsh. Further, we developed a soil salinity model to estimate the soil salinity in a salt marsh in Arakawa River. The vertical distribution of the soil salinity in the salt marsh was uniform and changed at almost the same time. The hydraulic conductivity of the soil, moreover, was high. The uniform distribution of salinity and high hydraulic conductivity could be explained by the vertical and horizontal transport of salinity through channels burrowed in the soil by organisms. By combining the shallow water model and the soil salinity model, the soil salinity of the salt marsh was well reproduced. The above results suggest that a stable brackish ecotone can be created in an artificial salt marsh using our numerical model as a design tool.


2009 ◽  
Vol 137 (10) ◽  
pp. 3339-3350 ◽  
Author(s):  
Ramachandran D. Nair

Abstract A second-order diffusion scheme is developed for the discontinuous Galerkin (DG) global shallow-water model. The shallow-water equations are discretized on the cubed sphere tiled with quadrilateral elements relying on a nonorthogonal curvilinear coordinate system. In the viscous shallow-water model the diffusion terms (viscous fluxes) are approximated with two different approaches: 1) the element-wise localized discretization without considering the interelement contributions and 2) the discretization based on the local discontinuous Galerkin (LDG) method. In the LDG formulation the advection–diffusion equation is solved as a first-order system. All of the curvature terms resulting from the cubed-sphere geometry are incorporated into the first-order system. The effectiveness of each diffusion scheme is studied using the standard shallow-water test cases. The approach of element-wise localized discretization of the diffusion term is easy to implement but found to be less effective, and with relatively high diffusion coefficients, it can adversely affect the solution. The shallow-water tests show that the LDG scheme converges monotonically and that the rate of convergence is dependent on the coefficient of diffusion. Also the LDG scheme successfully eliminates small-scale noise, and the simulated results are smooth and comparable to the reference solution.


Sign in / Sign up

Export Citation Format

Share Document