Skip to content

leonfoks/landslides-trigrs

 
 

Repository files navigation

README

Description

The Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Model (TRIGRS) is a Fortran program designed for modeling the timing and distribution of shallow, rainfall-induced landslides. The program computes transient pore-pressure changes, and attendant changes in the factor of safety, due to rainfall infiltration. The program models rainfall infiltration, resulting from storms that have durations ranging from hours to a few days, using analytical solutions for partial differential equations that represent one-dimensional, vertical flow in isotropic, homogeneous materials for either saturated or unsaturated conditions. Use of step-function series allows the program to represent variable rainfall input, and a simple runoff routing model allows the user to divert excess water from impervious areas onto more permeable downslope areas. The TRIGRS program uses a simple infinite-slope model to compute factor of safety on a cell-by-cell basis. An approximate formula for effective stress in unsaturated materials aids computation of the factor of safety in unsaturated soils. Horizontal heterogeneity is accounted for by allowing material properties, rainfall, and other input values to vary from cell to cell. This command-line program is used in conjunction with geographic information system (GIS) software to prepare input grids and visualize model results.

Purpose and limitations

TRIGRS is a tool to be used by investigators who have some knowledge and experience concerning landslide behavior. Selecting input parameters and interpreting results requires geologic and engineering judgment. The user should understand the theory and limitations behind TRIGRS, which are outlined in the documentation (Baum and others, 2002, 2008). The user must also be aware of the limitations of the digital elevation model, physical properties and hydrologic data that TRIGRS uses as input for analyses.

Typical Workflow

At the beginning of a project, the user would prepare a digital elevation model, slope grid, flow-direction grid, and physical properties zone grid files using Geographic Information System (GIS) software. Next, if needed, the user would run UnitConvert to put the grid files into a consistent system of units. Use GridMatch to ensure that the grid files are congruent. A common source of errors running TRIGRS is grids that contain the same number of cells, but the number or locations of no-data cells differ slightly because they have been generated from different sources. The user should then run TopoIndex using the digital elevation model and direction grid file as input to define the flow distribution pattern and compute weighting factors for distributing surface runoff. Routing and distribution of surface runoff is optional; however, TopoIndex will compute array sizes from the digital elevation model to be used by TRIGRS, and so the user should run TopoIndex once at the beginning of a project to determine array sizes. After preparing the topographic and physical properties data, converting the data to consistent units, computing array sizes and (if desired) preparing the flow distribution data, the project is ready for analysis using TRIGRS. The user may configure TRIGRS to establish a time-series response of shallow pore water and stability of shallow slope deposits to rainfall infiltration. Several initial runs of TRIGRS may be needed to verify calibration of the input parameters, the colluvial depth model and initial groundwater conditions.

What's included

This distribution includes source code files for the program TRIGRS and three companion utility programs, TopoIndex, GridMatch, and UnitConvert. It also includes sample data in the data folder and sample initialization files in the main folder. Empty folders for documentation, "doc" and executable binaries, "bin", are also included in the top-level directory.

Compiling the source code

We have chosen to use CMAKE to build the Fortran codes due to its easy cross platform nature. To do so, create a new folder in the root directory called "build". Type "cmake [other options] ../src" to generate the make files.

You will need a minimum cmake version of 3.12. You will also need a fortran compiler that can handle modules.

Then type "make" to compile the source.

This process will compile the serial version of the landslide_triggers code and, if MPI is detected, an MPI enabled version. Both versions need access to a HDF5 library with Fortran and high level capabilities. If during the compilation process the build fails to detect HDF5, set an environment variable called "HDF5_ROOT" that contains the directory of the HDF5 install you wish to link to.

For an MPI enabled build, specify the c, cxx, and fortran compilers during the cmake stage as follows, "cmake -DCMAKE_C_COMPILER=mpicc -DCMAKE_CXX_COMPILER=mpicxx -DCMAKE_Fortran_COMPILER=mpif90 ../src"

User Interface

TRIGRS and its companion utility programs, TopoIndex, GridMatch, and UnitConvert, run from the command line and have limited user interaction. Each program uses an initialization file that contains basic data needed to run the program as well as the names of other input files.

Latest version (August 2016)

This release, 2.1.0a, includes bug fixes, new output formats, and various optimizations to improve performance of the TRIGRS serial code. These improvements have also been incorporated into a new parallel implementation of the TRIGRS program (Alvioli and Baum, 2016). We have parallelized the four time-demanding execution modes of TRIGRS, namely both the saturated and unsaturated model with finite and infinite soil depth options, within the Message Passing Interface (MPI) framework. Performance gain with respect to the serial code was tested both on commercial hardware and on a high-performance multi-node machine. We also compared results of the parallel code against results of the serial code for a large study area in Colorado to verify accuracy of the results.

Testing

Throughout its development, the code has passed through various kinds of testing to verify that it (1) reproduces results of the basic formulas for transient infiltration and slope stability and correctly maps grid cells to their spatial locations (Baum and others, 2002, 2008, 2010, 2013), (2) that it produces results that are consistent with actual case studies (Savage and others, 2003, 2005, Salciarini and others, 2006; Godt and others, 2008a, 2008b; Baum and others, 2010, 2011; Raia and others 2014; Gioia and others, 2015).

References cited

  • Alvioli Massimiliano, and Baum, R.L., 2016, Parallelization of the TRIGRS model for rainfall-induced landslides using the message passing interface: Environmental Modeling & Software, v. 81, p. 122 - 135, doi: 10.1016/j.envsoft.2016.04.002.

  • Baum, R.L., Savage, W.Z., and Godt, J.W., 2002, TRIGRS--A FORTRAN Program for Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Analysis: U.S. Geological Survey Open-File Report 02-0424, 35 p., 2 appendices.

  • Baum, R.L., Savage, W.Z., and Godt, J.W., 2008, TRIGRS--A Fortran program for transient rainfall infiltration and grid-based regional slope-stability analysis, version 2.0: U.S. Geological Survey Open-File Report, 2008-1159, 75 p.

  • Baum, R. L., Godt, J.W., and Savage, W. Z., 2010, Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration: Journal of Geophysical Research, Earth Surface. v. 115, F03013, doi:10.1029/2009JF001321.

  • Baum, R.L., Godt, J.W., and Coe, J.A., 2011, Assessing susceptibility and timing of shallow landslide and debris flow initiation in the Oregon Coast Range, USA, In Genevois, R. Hamilton, D.L. and Prestininzi, A. (eds.) Proceedings of the Fifth International Conference on Debris Flow Hazards Mitigation--Mechanics, Prediction, and Assessment, Padua, Italy, June 7-11, 2011, p. 825-834. Rome: Casa Editrice Universitá La Sapienza (doi: 10.4408/IJEGE.2011-03.B-090).

  • Baum, R.L., and Godt, J.W., 2013, Correction to "Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration": Journal of Geophysical Research, Earth Surface, v. 118, 20100, doi:10.1002/jgrf.20100.

  • Gioia, Eleonora, Speranza, Gabriella, Ferretti, Maurizio, Godt, Jonathan W., Baum, Rex L., Marincioni, Fausto, 2015, Application of a process-based shallow landslide hazard model over a broad area in Central Italy: Landslides, doi: 10.1007/s10346-015-0670-6, p. 1-18.

  • Godt, J.W., Baum, R.L., Savage, W.Z., Salciarini, D., Schulz, W.H., and Harp, E.L., 2008a, Transient deterministic shallow landslide modeling--Requirements for susceptibility and hazard assessments in a GIS framework: Engineering Geology, v. 102, no. 3-4, Landslide Susceptibility, Hazard and Risk Zoning for Land Use Planning, p. 214-226, DOI: 10.1016/j.enggeo.2008.03.019.

  • Godt, J.W., Schulz, W.H., Baum, R.L., and Savage, W.Z., 2008b, Modeling rainfall conditions for shallow landsliding in Seattle, Washington, in Baum, R.L., Godt, J.W., and Highland, L.M., eds., Engineering geology and landslides of the Seattle, Washington, area: Geological Society of America Reviews in Engineering Geology, v. XX, p.137-152, doi: 10.1130/2008.4020(08)

  • Raia, S., Alvioli, M., Rossi, M., Baum, R. L., Godt, J. W., and Guzzetti, F., 2014, Improving predictive power of physically based rainfall-induced shallow landslide models: a probabilistic approach: Geoscientific Model Development, Vol. 7, p. 495-514, doi:10.5194/gmd-7-495-2014, 2014.

  • Salciarini, D., Godt, J.W., Savage, W.Z., Conversini, P., Baum, R.L., and Michael, J.A., 2006, Modeling regional initiation of rainfall-induced shallow landslides in the eastern Umbria Region of central Italy: Landslides, v. 3, no. 3, p. 181-194, DOI 10.1007/s10346-006-0037-0

  • Savage, W.Z., Godt, J.W., and Baum, R.L., 2003, A model for spatially and temporally distributed shallow landslide initiation by rainfall infiltration, in Rickenmann, D. and Chen, C., eds., Debris-Flow Hazards Mitigation--Mechanics, Prediction, and Assessment: Rotterdam, Millpress (Proceedings of the 3rd International conference on Debris Flow Hazards, Davos, Switzerland, September 10-13, 2003), p. 179-187.

  • Savage, W.Z., Godt, J.W., and Baum, R.L., 2004, Modeling time-dependent aerial slope stability, in Lacerda, W.A., Erlich, M., Fontoura, S.A.B., and Sayao, A.S.F., eds., Landslides--evaluation and stabilization, Proceedings of the 9th International Symposium on Landslides: London, A.A. Balkema Publishers, v. 1, p. 23-36.

About

No description, website, or topics provided.

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages

  • Fortran 95.2%
  • CMake 4.4%
  • Makefile 0.4%