Development of new equation to estimate the maximum soil depth by using the safety factor

Similar documents
(h+ ) = 0, (3.1) s = s 0, (3.2)

ANALYSIS AND MODELLING OF RAINFALL EVENTS

Table of Content. c 1 / 5

Population Dynamics Definition Model A model is defined as a physical representation of any natural phenomena Example: 1. A miniature building model.

A Mathematical Model for Unemployment-Taking an Action without Delay

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)

Comparing the Pre-image and Image of a Dilation

Activities. 4.1 Pythagoras' Theorem 4.2 Spirals 4.3 Clinometers 4.4 Radar 4.5 Posting Parcels 4.6 Interlocking Pipes 4.7 Sine Rule Notes and Solutions

Hyers-Ulam stability of Pielou logistic difference equation

COMPARISON BETWEEN TWO FRICTION MODEL PARAMETER ESTIMATION METHODS APPLIED TO CONTROL VALVES. Rodrigo Alvite Romano and Claudio Garcia

Review Topic 14: Relationships between two numerical variables

AP Calculus AB Unit 4 Assessment

99/105 Comparison of OrcaFlex with standard theoretical results

Chem Homework 11 due Monday, Apr. 28, 2014, 2 PM

, g. Exercise 1. Generator polynomials of a convolutional code, given in binary form, are g. Solution 1.

Electromagnetism Notes, NYU Spring 2018

Parabola and Catenary Equations for Conductor Height Calculation

Delay Variability at Signalized Intersections

T b a(f) [f ] +. P b a(f) = Conclude that if f is in AC then it is the difference of two monotone absolutely continuous functions.

16z z q. q( B) Max{2 z z z z B} r z r z r z r z B. John Riley 19 October Econ 401A: Microeconomic Theory. Homework 2 Answers

Part 4. Integration (with Proofs)

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx,

Iowa Training Systems Trial Snus Hill Winery Madrid, IA

Some Aspects of Non-Orthogonal Stagnation-Point Flow towards a Stretching Surface

Acceptance Sampling by Attributes

CHENG Chun Chor Litwin The Hong Kong Institute of Education

1.3 SCALARS AND VECTORS

Solutions to Assignment 1

New Expansion and Infinite Series

Modelling the Electrolyte Flow in a Full-scale Copper Electrorefining Tankhouse Cell

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

Engr354: Digital Logic Circuits

On the Scale factor of the Universe and Redshift.

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

The Emission-Absorption of Energy analyzed by Quantum-Relativity. Abstract

] dx (3) = [15x] 2 0

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Terminal Velocity and Raindrop Growth

Lecture 1 - Introduction and Basic Facts about PDEs

Lesson 2: The Pythagorean Theorem and Similar Triangles. A Brief Review of the Pythagorean Theorem.

Applications of Definite Integral

1 Bending of a beam with a rectangular section

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O

Applications of Definite Integral

8 THREE PHASE A.C. CIRCUITS

GM1 Consolidation Worksheet

CH 17: Flexible Mechanical Elements

Distance Measurement. Distance Measurement. Distance Measurement. Distance Measurement. Distance Measurement. Distance Measurement

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums

3.1 Review of Sine, Cosine and Tangent for Right Angles

CS 573 Automata Theory and Formal Languages

QUADRATIC EQUATION. Contents

Lecture Notes No. 10

Forces on curved surfaces Buoyant force Stability of floating and submerged bodies

Intermediate Math Circles Wednesday 17 October 2012 Geometry II: Side Lengths

ISOTHERMAL REACTOR DESIGN (4) Marcel Lacroix Université de Sherbrooke

a) Read over steps (1)- (4) below and sketch the path of the cycle on a P V plot on the graph below. Label all appropriate points.

8.2 ESTIMATING DETENTION VOLUMES

COMPARISON OF DIFFERENT PROCEDURES TO PREDICT UNSATURATED SOIL SHEAR STRENGTH. S.K. Vanapalli and D.G. Fredlund 1

Vadose Zone Hydrology

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4.

5.7 Improper Integrals

Global alignment. Genome Rearrangements Finding preserved genes. Lecture 18

Job No. Sheet 1 of 8 Rev B. Made by IR Date Aug Checked by FH/NB Date Oct Revised by MEB Date April 2006

Scientific notation is a way of expressing really big numbers or really small numbers.

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Intro to Nuclear and Particle Physics (5110)

MATH Final Review

Estimation of Global Solar Radiation in Onitsha and Calabar Using Empirical Models

Learning Objectives of Module 2 (Algebra and Calculus) Notes:

Defining Areas with Nil Landslide Hazard is a step toward a Comprehensive Landslide Loss Model

(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P.

Chemical Equilibrium

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

International Industrial Informatics and Computer Engineering Conference (IIICEC 2015)

Übungen zur Theoretischen Physik Fa WS 17/18

Section 3.6. Definite Integrals

Section 1.3 Triangles

Appendix C Partial discharges. 1. Relationship Between Measured and Actual Discharge Quantities

Metodologie di progetto HW Technology Mapping. Last update: 19/03/09

Environmental Science

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Songklanakarin Journal of Science and Technology SJST R1 Thongchan. A Modified Hyperbolic Secant Distribution

THE ASYMMETRY OF COASTAL WATER LEVEL RESPONSE TO LANDFALLING HURRICANES SIMULATED BY A THREE-DIMENSIONAL STORM SURGE MODEL

Electromagnetic-Power-based Modal Classification, Modal Expansion, and Modal Decomposition for Perfect Electric Conductors

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Introduction to Olympiad Inequalities

Example. Have precipitation and streamflow data, need to estimate losses

Novel Fiber-Optical Refractometric Sensor Employing Hemispherically-Shaped Detection Element

Maintaining Mathematical Proficiency

Non-Linear & Logistic Regression

Modeling of Catastrophic Failures in Power Systems

Tests for the Ratio of Two Poisson Rates

Application of the theory of compound cores for the assessment of stress pattern in the cross section of a strengthened beam column

Ellipses. The second type of conic is called an ellipse.

2 Fractions and ratios

Transcription:

Lndslides nd Engineered Slopes. Experiene, Theory nd Prtie Avers et l. (Eds) 2016 Assoiione Geoteni Itlin, Rome, Itly, ISBN 978-1-138-02988-0 Development of ne eqution to estimte the mximum soil depth by using the sfety ftor G.P. Mihel Snt Ctrin Stte University, Ibirm, Bril M. Kobiym Federl University of Rio Grnde do Sul, Porto Alegre, Bril ABSTRACT: Soil depth is n importnt prmeter hih is hrd to be mesured nd estimted in hydrogeomorphi nd pedologil studies. Furthermore, the soil depth is one of the most unknon vribles in the hillslope system. Therefore, the objetive of the present study s to develop ne eqution to estimte the mximum soil depth by using the ftor of sfety. For this development, uniform hydrologil model s ombined to n infinite slope model. The ne eqution is lled MEMPS (Modelo de Estimtiv d Máxim Profundidde do Solo). Differently from other equtions to lulte mximum soil depth, MEMPS does not require the ssumptions of omplete sturtion or totl bsene of ter in the soil. Thus, MEMPS represents ondition loser to the relity hen ompred to other equtions. Moreover, due to the inorportion of hydrogeomorphi prmeters in its formultion, MEMPS quires fitness to reflet the influene of hydrogeomorphi proesses on the soil depth estimtion. 1 INTRODUCTION Soil depth, defined here s the vertil distne beteen the soil surfe nd the bedrok, is n importnt prmeter hih is hrd to be mesured nd estimted in hydrogeomorphi nd pedologil studies. Furthermore, the soil depth is one of the most unknon vribles in the hillslope system (Ctni et l., 2010). Tht is hy there re vrious ttempts to estimte the soil depth. These ttempts n be lssified into to: field methods (Piegri et l., 2009, Bery, 2013) nd mthemtil methods (Sulnier et l., 1997, Segoni et l., 2012). Eh ttempt hs its on good nd bd points. Borg et l. (2002) performed sensitivity nlysis of slope stbility model nd demonstrted tht depth vrition n signifintly hnge the slope stbility results. It mens tht the suessful mpping of lndslide hrd res depends upon the orret determintion of the sptil distribution of soil depth in bsin. As the soil depth mesurement in hole bsin is quite diffiult nd time-onsuming (Dietrih et l., 1995), it is usully required to use some mthemtil formultions hih estimte the soil depth distribution. The proess-bsed methods re those tht use physilly-bsed equtions to desribe vrious mehnisms ffeting the soil formtion nd evolution. Though generting stisftory results, they require tht the physil proesses ting in the study re re ell knon nd tht input prmeters dt re urte. In this ontext, the objetive of the present study s to develop ne eqution to estimte the mximum soil depth by using the Ftor of Sfety (FS). This ne eqution resulted from the ombintion of n infinite slope stbility model ith uniform hydrologil model. The present study shos the eqution development nd omments its limittion use. 2 THEORETICAL DEVELOPMENT OF MEMPS Aording to Selby (1993), the vlue of FS n be formulted: FS s g g + ( ρ os θ ρ h os θ) tnφ ρ g osθ s 2 2 (1) here FS ftor of sfety; soil ohesion; ρ ter density; soil depth; h height of the ter olumn inside the soil lyer; ρ s soil density; g grvittionl elertion; θ slope ngle; nd f soil internl frition ngle (Figure 1). Iid (1999) exerised limit equilibrium nlysis for n infinite slope model, nd defined the ritil 1417

of ter in the soil olumn; nd (ii) its omplete sturtion. In ft, both onditions hrdly our. The present study, therefore, extended his theory in order to tret the more rel sitution, by ombining n infinite slope stbility model to uniform hydrologil model. Then, the present study modifies the Eqution 1 nd obtins: Figure 1. Representtion of n infinite slope model. P eight of the soil; τ sher stress; σ norml stress; u pore-pressure. (From Mihel et l., 2014). soil depth s the mximum vlue of the depth tht n settle on hillslope, by onsidering tht the FS vlue t prtiulr lotion does not reh loer vlues thn 1. In the se of ter bsene, here soil lyer is perfetly dry, the ritil depth ( 0 ) s desribed by: ρ g os θ tnθ tnφ 0 2 s ( ) (2) here 0 ritil depth for dry soil. The 0 vlue n be lulted, but it is never rehed, sine the soil depth sloly inreses over time nd periodi storms produe ertin hillslope sturtion (i.e. the ssumption of bsene of ter is not stisfied), using lndslides. Therefore, Iid (1999) proposed, for the se tht the soil is totlly sturted, nother eqution to determine the ritil depth ( 1 ): os θ g ( tnθ tnφ)+ ρ g tnφ 1 2 (3) here 1 ritil depth for the sturted soil. If the tul soil depth is smller thn 1, lndslide ill never our even ith very intensive rinfll. In this se, this lotion quires immunity sitution for ertin period. If the soil depth is greter thn 1, the lndslide n our hen the soil sturtion rehes threshold vlue. Iid (1999) demonstrted tht his equtions (Equtions 2 nd 3) hd good performne through omprison ith the vlues mesured in the field. As bove-mentioned, Iid (1999) determined the mximum soil depth, by onsidering to extreme onditions of soil moisture: (i) the totl bsene h s g ρ osθ 1 ρ θ tnφ os FS (4) The term h/ in Eqution 4 represents the sturtion degree of the hillslope soil. In stedy stte hydrologil model, hih is ommonly lled TOPOG, O Loughlin (1986) defined the etness index () ith the folloing eqution: h bt (5) here q uniform stedy stte rehrge; upstrem bsin re t one point; b ontour length; nd T soil trnsmissivity. Obviously, the vlue of is limited to 1, nd this ondition mens tht the height of the ter olumn inside the soil is equl to the soil depth. When the pplition of Eqution 5 results in vlue lrger thn one, is there overlnd flo. In Eqution 5 sturted hydruli ondutivity (Ks) is onsidered onstnt throughout the soil profile, nd then T n be lulted: T Ks osθ (6) By substituting the Eqution 6 in the Eqution 5, it is obtined: h or b Ks osθ h b Ks osθ (7) (8) As the ter blne of the soil lyer is onsidered permnent, the Eqution 7 is ble to estimte n verge behvior of the soil sturtion level over time. The mximum vlue of h is equl to. When the pplition of Eqution 8 results in h vlue lrger thn, the soil is ompletely sturted nd there is overlnd flo. Then, substituting the Eqution 7 in the Eqution 4, the folloing expression is obtined: 1418

q s g ρ os θ ρ osθ 1. b Ks osθ tn φ FS (9) The ssumption, tht the stbility threshold in hillslope ours hen FS is equl to 1, permits to estimte the mximum depth tht the soil n hieve ithout using instbility, only bsed on soil strength hrteristis, geomorphi hrteristis of the bsin nd hydrologil situtions tht tke ple. Hene, the Eqution 9 n be modified s follos: The present study lls the Eqution 13 Modelo de Estimtiv d Máxim Profundidde do Solo (MEMPS) hih mens model to estimte the mximum soil depth, in Portuguese. 1 ρ os θ. ρ g b Ks tnφ s osθ osθ 1 < (10) Isolting the term, the Eqution 10 beomes: ρ tnφ g osθ b Ks < osθ tnφ (11) The Eqution 11 performs n estimte of the mximum depth rehed by the soil under ertin hydrologil onditions. Then, by dividing the Eqution 11 by osθ nd treting the mximum soil depth, the vlue of the mximum soil depth ( ) n be obtined: tnφ ρ 2 g os θ b Ks osθ tnθ tnφ (12) The term in the brkets in the Eqution 12 represents the ter height of the soil lyer (h) of the Eqution 8. The mximum vlue of h must be equl or less thn the vlue of. Thus, it n be sid tht, under threshold blne ondition hen h beomes equl to, the ritil depth beomes 1 of the Eqution 3. Then, the Eqution 12 n be modified: tnφ ρ min s g s, 2 1 ρ os θ ρ b Ks osθ tnθ tnφ (13) 3 COMPARISON BETWEEN MEMPS AND IIDA (1999) The MEMPS (Eqution 13) hs intermedite hrteristis beteen to equtions (Equtions 2 nd 3) of Iid (1999) hih used n inversion of the FS eqution to limit the mximum soil depth physilly fesible on slope. The differene beteen the equtions of Iid (1999) nd the MEMPS is grphilly illustrted in Figure 2. For elborting these grphs, the lultion s done under the folloing onditions: 11.9 kp, q 0.005 m/dy, 300 m², b 5 m, f 30.5, Ks 0.38 m/dy, ρ s 1800 kg/m 3, nd ρ 1000 kg/m 3. Considering tht the im of the present study is just to demonstrte the performne of Eqution 13, nd not Figure 2. Vrition of mximum soil depth vlues estimted by Iid s equtions nd MEMPS. 1419

to modelling the tul mximum soil depth in slope, some of these prmeters ere defined from Mihel et l. (2014) nd other ones ere rbitrrily defined. The only re in the definition of these prmeters s to selet vlues of q,, b nd Ks tht ere not ble to generte omplete sturtion in the slope. It is lerly observed tht the Eqution 13 hs its pplition limit. When the stedy stte rehrge (q) used to define the soil depth distribution ith MEMPS results in omplete sturtion of the soil, the results of MEMPS re equl to those ith Eqution 3. When the vlue of q is set to ero (i.e. there is no ter in the bsin) the results of MEMPS re equl to those ith Eqution 2. Thus, MEMPS n desribe three different situtions (bsene of ter, omplete sturtion nd n intermedite ondition) ith only one eqution. Beuse it, the performne of the model depends on the orret estimton of q. Furthermore, this model is only pplible on slopes tht re steeper thn f, i.e., θ > f. In Figure 2, it is observed tht the etness ondition used by MEMPS is intermedite beteen omplete sturtion nd totl bsene of ter. For lrge vlues of soil depth, the etness ondition in MEMPS s more similr to bsene of ter. Insted, for smll vlues of soil depth, the etness ondition in MEMPS s more similr to sturtion. The prtilly sturtion ondition in the one bove the ter tble s not took into ount in this study. When nlying the mximum stble soil depth, only the onditions developed in possible rupture surfe re relevnt. In other ords, only the onditions in the interfe beteen soil nd bedrok re tking into ount, beuse this one is the one more suseptible to filure. In the pplition of MEMPS (Eqution 13), it is neessry to set uniform stedy stte rehrge rte, in order to mimi period of preipittion bove the regulr onditions, but ithout the pbility to trigger lndslides. The ondition ssumed genertes pttern of soil depth distribution in hih the hillslope tends to its stbility during norml or slightly more thn norml rinfll sesons. When hevy rinfll tht is more thn norml tkes ple, the soil grdully redues its FS vlue nd n reh unstble onditions. Thus, for exmple, n initil ondition of relted to the nnul men rinfll ould generte et ondition insuffiient to use the instbility. When hevy rinfll similr to tht used to define the soil depth distribution ours in the bsin, the sitution turns grdully from stble ondition to the threshold. If the rinfll persists nd the etness exeeds the ondition used to define the soil depth, the sitution evolves from limit equilibrium to unstble onditions. The hrteristi of MEMPS, in hih initilly there re stble onditions tht grdully turn in to unstble onditions ith the rinfll event, gives to this model the qulity to be quite dequte to define initil onditions to slope stbility modelling. Most of times, in the slope stbility modelling, the initil onditions re of instbility over lrge prt of the terrin (minly in steep onve res), even before the ourrene of the triggering event. In regions here the lndslides re triggered by hydrologil ftors, the initil onditions (before the rinfll or triggering event ourrene) of the terrin must be stble, ht is hieved ith pplition of MEMPS. Thus, ith this more relisti senrio, the slope stbility modelling ould be ble to desribe the speifi proess involved in the lndslide triggering. Furthermore, in lndslide rning systems, it is neessry to define the moment in hih the terrin ill turn from stble to unstble onditions. Thus, in these ses, n initil ondition of stbility is required. Therefore, MEMPS n be useful tool for soil depth nd lndslide hrd mpping, nd lndslide rning systems. 4 FINAL CONSIDERATIONS The present study developed ne eqution, lled MEMPS, to estimte the mximum soil depth through ombintion of stedy stte hydrologil model ith n infinite slope stbility model. This eqution does not need to hve ny of ssumptions of Iid (1999), representing ondition loser to the relity. In future, it should be verified ith field dt. Even if the MEMPS n be only used for res here the slope ngle is lrger thn the soil internl frition ngle, i.e., mountinous bsins, it is very useful. One of the gols of the development of this eqution is to rete more relisti initil senrios for slope stbility modelling. In these ses, MEMPS n be very useful both in to optimie lndslide modelling nd to improve lndslide rning systems. REFERENCES Bery, A.A. 2013. High resolution in seismi refrtion tomogrphy for environmentl study. Interntionl Journl of Geosienes 4: 792 796. Borg, M., Dll Fontn, G., Gregoretti, C., Mrhi, L. 2002. Assessment of shllo lndsliding by using physilly bsed model of hillslope stbility. Hydrologil Proesses 16: 2833 2851. Ctni, F., Segoni, S., Florni, G. 2010. An empiril geomorphology-bsed pproh to the sptil predition of soil thikness t thment sle. Wter Resoures Reserh 46: 1 15. 1420

Dietrih, W.E., Reis, R., Hsu, M.L., Montgomery, D.R. 1995. A proess-bsed model for olluvil soil depth nd shllo lndsliding using digitl elevtion dt. Hydrologil Proesses 9: 383 400. Iid, T. 1999. A stohsti hydro-geomorphologil model for shllo lndsliding due to rinstorm. Cten 34:293 313. Mihel, G.P., Kobiym, M., Goerl, R.F. 2014. Comprtive nlysis of SHALSTAB nd SINMAP for lndslide suseptibility mpping in the Cunh River bsin, southern Bril. Journl of Soils nd Sediments 14: 1266 1277. O Loughlin, E.M. 1986. Predition of surfe sturtion ones in nturl thments by topogrphi nlysis. Wter Resoures Reserh 22: 794 804. Piegri, E., Ctudell, V., Di Mrio, R., Milno. L., Niodemi, M., Soldovieri, M.G. 2009. Eletril resistivity tomogrphy nd sttistil nlysis in lndslide modelling: A oneptul pproh. Journl of Applied Geophysis 68(2): 151 158. Sulnier, G.M., Beven, K., Obled, C. 1997. Inluding sptilly vrible effetive soil depths in TOPMODEL. Journl of Hydrology 202: 158 172. Segoni, S., Rossi, G., Ctni, F. 2012. Improving bsin sle shllo lndslide modelling using relible soil thikness mps. Nturl Hrds 61: 85 101. Selby, M. 1993. Hillslope mterils nd proesses. 2nd ed. Oxford University Press, Oxford. 1421