MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 731042 10.1155/2013/731042 731042 Research Article Nonsmooth Dynamic Behaviors Inherited from an Ecohydrological Model: Mutation, Bifurcation, and Chaos Lin Mu 1,2 Tian Fuqiang 2 Hu Heping 2 Liu Dengfeng 3 Li Yongping 1 School of Applied Mathematics Central University of Finance and Economics Beijing 100081 China cufe.edu.cn 2 Department of Hydraulic Engineering State Key Laboratory of Hydroscience and Engineering Tsinghua University Beijing 100084 China tsinghua.edu.cn 3 Key Laboratory of Northwest Water Resources and Environmental Ecology of MOE Xi'an University of Technology Xi'an 710048 China xaut.edu.cn 2013 19 3 2013 2013 17 12 2012 28 01 2013 2013 Copyright © 2013 Mu Lin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The existence of nontrivial dynamic behaviors in a hydrological system is intensively discussed in the literature. However, most of the work has been done from the nonlinear data analysis perspective, with only a few exceptions, due to the mathematical difficulties for theoretical analysis. In this study, a simple but comprehensive enough ecohydrological model with the pulsed atmospheric forcing was developed from the process analysis perspective. The model was then utilized to analyze the non-trivial dynamic behaviors in a coupled ecohydrological system qualitatively and numerically. Our results confirm the existence of multiple stationary states discussed by many researchers. Furthermore, parameter bifurcation was studied and the phenomenon of mutation is found to be rather common. Also, the chaotic characteristic of the system state is obtained under some specific parameters. Parts of these behaviors were seldom reported through the deterministic dynamic analysis done previously.

1. Introduction

As stated by Peterson et al. , the dynamic behavior of a hydrological system is always assumed to be primary in traditional hydrology; that is, at the cessation of a transient hydrological disturbance of any magnitude, the system will return to the same stable state and thus show an infinite resilience. This hypothesis is challenged by the dynamic theory, especially the nonlinear theory raised in meteorology . In fact, quite a few studies have been devoted to the “nontrivial” dynamic characters in hydrological systems in the past several decades . For example, chaos has been intensively studied in many hydrology-related processes including rainfall, river flow, rainfall-runoff, lake volume, temperature, pressure, and wind velocity (see the summary by Sivakumar ). These studies are, however, principally based on the nonlinear data analysis method alone, which can obtain the statistical features from the existing data but cannot indicate the causal relationship between physical processes and dynamic phenomena.

Another point of view, as this study mainly concerns, uses process-based model incorporating intertwined interactions and feedbacks to explore the nontrivial dynamic behaviors. For example, considering the coupled system including the upper soil layer and the atmosphere planetary boundary layer, D'Andrea et al.  found the bimodal distribution of soil moisture resulting from the existence of multiple equilibria in the continental water balance. Ridolfi et al.  developed a conceptual framework to study the interaction and positive feedback in wetland forests and riparian ecosystems, which shows that the feedbacks may lead to the emergence of multiple stable states. Based on the field observations from savanna ecosystems in Kalahari, D'Odorico et al.  found that the difference between the soil moisture in the canopy and intercanopy spaces is related to the aridity level, and the system may exhibit two stable states corresponding to the conditions with and without tree canopy cover. With positive feedbacks between vegetation and soil moisture in arid ecosystems, Borgogno et al.  found that random rainfall fluctuation may turn the bistable deterministic system into a stochastic system with only one statistically stable state in an analytical model. In order to challenge the potentially significant assumption of the single stable state in the ecohydrological models, Peterson et al.  developed a numerical distributed ecohydrological model and found that the multiple stable states may arise from the reduction of leaf area index as the saline water table approaches the surface. Also, they employed limit cycle continuation to locate the unstable stationary state in the phase space. Runyan and D'Odorico  developed a modeling framework relating vegetation-soil salinity feedbacks to the emergence of multiple stable states in the underlying dynamics and the results show that the presence of a strong feedback can result in bistable dynamics for a wide range of environmental conditions.

The above process-based work demonstrates the existence of multiple stationary states in some specific ecohydrological systems, which depends on no particular data or conditions. But more sophisticated dynamic behaviors, such as mutation and chaos, have not been explored yet. Nevertheless, Rodriguez-Iturbe and the collaborators [12, 13] conducted a study with the similar purpose but for the hydrometeorological system operating at the climate scales (i.e., continental spatial scale and annual temporal scale). They analyzed the soil moisture balance equation and demonstrated that soil moisture states could present abundant nontrivial dynamic behaviors including bimodal states, fixed points, limit cycles, and chaos due to the feedback between soil moisture and precipitation (i.e., local recycling of moisture). For the ecohydrological system running at much smaller spatial and temporal scales interested in this study, the dominating feedback is not the local recycle of moisture, but the vegetation growth fed by soil moisture. The evolution of ecohydrological system has recently drawn a lot of attentions and its dynamic behavior deserves further study, which is the purpose of this study.

The remainder of this paper is organized as follows. In Section 2, a simple but comprehensive enough ecohydrological model is established under pulsed atmospheric forcing. Dynamic analysis follows in Section 3, and the so-called nontrivial dynamic behaviors are explored by the way of numerical simulation. The paper is concluded in Section 4.

2. Framework of the Ecohydrological Model 2.1. Functional Expression of Nonsmooth Rainfall Process

As mentioned above, Rodriguez-Iturbe et al. conducted the notable work in the dynamic analysis of hydrometeorologic process. However, the annual rainfall amount in their analysis is regarded as a simple function (in the mathematical sense) of advective precipitation, local evaporation, and local soil moisture, and hence the smaller scale rainfall variability and nonsmooth feature (see below) are ignored. Such description of rainfall process is inapplicable to ecohydrological system running at a smaller spatiotemporal scale. For the detail study of hydrological processes at the finer scale, the Poisson process with white noise is often utilized to describe the stochastic properties of rainfall [14, 15]. However, such stochastic manner will go against the qualitative analysis of the dynamic process due to the mathematical complexity.

A deterministic functional expression of rainfall with moderate complexity is necessary for dynamic analysis of hydrological process. From the mathematical perspective, the rainfall intensity is a nonsmooth function of time ; that is, a pulsing function which means the time derivative is not continuous everywhere. In general, nonsmooth modeling has been widely used in many different areas  for its intuitional advantage; although nonsmooth terms can introduce additional mathematical difficulties for dynamic analysis [20, 21]. In hydrology, a simplified nonsmooth function, that is, a rectangular pulsing function, is frequently used for theoretical analysis of ecohydrological or hydrometeorologic processes. For example, a simple step function was used to describe the pulsing rainfall forcing by Collins and Bras  to examine how the gradient of mean annual precipitation is weaved in the topography of water-limited ecosystems; an ideal quasiperiodic pulsing precipitation was utilized by Lu  to study the physical basis of daily flood and drought index; a periodic step function was adopted by Robinson and Sivapalan  to identify different hydrological regimes according to the interactions between rainfall variability and catchment response time within a theoretical framework. Following these studies, a quasiperiodical rectangular pulsing function was chosen to represent real stochastic rainfall process by keeping between-storm temporal structure while ignoring more variability for a preliminary dynamic analysis of ecohydrological system.

Remarkably, the chosen pulsing function is just like the water application manner in an automatic irrigation system under the equilibrium state, see also Section 3 and Figure 8(c). The so-called automatic irrigation  means the duration of irrigation is controlled by the real-time observation on soil moisture content, with the switching bounds predetermined. As the dynamic analysis in this study focuses on the equilibrium state of ecohydrological system, we will use automatic irrigation scheme to represent the quasiperiodical pulsing function. Mathematically, such substitution of automatic irrigation scheme for stochastic rainfall process can ensure that the model is deterministic and autonomous in the resultant ordinary differential equations (ODEs), and therefore the dynamic phenomena discovered in the model depends not on any particular selection of external rainfall time series but on the internal system configuration. Practically, the automation has been widely used in water-saving irrigation systems [25, 26] under dry as well as wet conditions , and therefore we can assume a no-rainfall area and focus on the automatic irrigation system only, for the irrigation amount is much larger than the rainfall in the arid especially hyperarid regions .

2.2. Ecohydrological Model

An ecohydrological model is then introduced as follows. The schematic diagram is shown as Figure 1. As mentioned above, it can be assumed to be in the hyperarid area where water supply comes from irrigation alone. Our work imitates the tank model  to describe the storage and release of aquifer.

The schematic diagram of the model.

As shown in Figure 1, the height h(t)(cm) stands for the storage per unit area at time t(d). Roughly speaking, it is equivalent to water content in the soil. The tank has only one side outlet with height H(cm). The side runoff from tank is assumed to be linear with respect to tank water storage excess relative to initial height H, defined by a positive coefficient c. The surface runoff is ignored. The term v(t) (quantity of grasses per square meter) stands for biomass, that is, the vegetation density per unit surface area at time t. Generally, the evaporation and transpiration are related to biomass by α(h,v)v. Assume that α is constant for a preliminary study here. Then, the dynamical equation of water storage can be written as follows: (1)h˙=p-αv-{c(h-H),  ifhH,0,  ifh<H,ifi(t)=1,h˙=-αv-{c(h-H),  ifhH,0,  ifh<H,ifi(t)=-1. Here we divide the whole process into two phases and introduce a variable  i(t) with value domain {-1,+1} to depict them; that is, i(t)=1 means irrigation at time t, with constant intensity p(cm/d), while i(t)=-1 means no irrigation at time t, that is, irrigation intensity is 0.

Another important dynamical equation is the evolution equation of the vegetation, which can be described as (2)v˙=β(h,v)v, where β(h,v) indicates the varying rate per unit vegetation. It is a complicated problem to describe β function in ecology. However, a qualitative function can be described approximatively as follows: (3)with  β(h,v)={-β0,ifhv[0,μ1-δ1),β0+β1δ1+δ2(hv-μ1),ifhv[μ1-δ1,μ1+δ2),β1,ifhv[μ1+δ2,μ2-δ3),-β0+β1δ3+δ4(hv-μ2),ifhv[μ2-δ3,μ2+δ4),-β0,ifhv[μ2+δ4,+). The framework is similar to (the linearization of) the work in . It means that the plants would die down when h/v is too large or too small, that is, v˙<0 . A sketch map of relation between “breed” rate and available water amount (for unit plant) is shown as Figure 2. Here β1 (d-1) is the maximum growth rate per unit grass, while -β0 (d-1) is the rate of natural death. Also, there are two equilibrium states h/v=μ1(cm/n), and h/v=μ2(cm/n) during which the grass density keeps constant. Positive coefficients δk(cm/n), k=1,2,3,4, are used for separating the suitable growth parameter region.

The schematic diagram of an ideal β function.

2.3. Automatic Irrigation Scheme

To summarize, the model has three critical variables: h(t), v(t), and i(t). The value sets of h(t) and v(t) are continuous with their dynamic functions (1)–(3), while the value set of i(t) has only two discrete elements 1 and −1. We only have to define a dynamic function of i(t) to complete our model. As mentioned in Section 2.1, we use a scheme of automatic irrigation switching rule based on the soil moisture content indicated by h: (4)i(t+)={-i(t),ifi(t)=-1,h(t)w0,ori(t)=1,h(t)w1,i(t),else, where i(t+)=limΔt0+i(t+Δt). Bounds of water content w0(cm) and w1(cm) are introduced, w0<w1. In agricultural water management practice, w0 could be chosen near the wilting point, and w1 could be chosen near the field capacity . Automatic control can be realized by combining the equipments of soil moisture monitoring, pipe irrigation, and switching . Finally, an executable conceptual ecohydrological model is obtained by combining (1)–(4) together.

2.4. Qualitative Analysis on the Model

Plotting a trajectory {(h(t),v(t),t)t0} starting from (h(0),v(0)) in h-v space, we get a phase diagram; for example, see Figure 3. Thereinto, the state point is the pair (h(t),v(t)), and the limit stationary state (the state when t+) appears as a period “circle” motion between h[w0,w1]. Generally, there are two stable stationary states. v(t)0 is always a stationary state, which means that given plants cannot be generated without progenitor. Note the motion {(w0,0)(w1,0)(w0,0)} as the zero circle. However in most cases there should be another periodic orbit with v0 for some parameters and initial values. It means that plants can be grown steadily under certain conditions. Note this stable stationary state as the nonzero circle. Because the vegetation density changes periodically, we concern the average vegetation density in a period of nonzero circle, noted as v¯. Furthermore, note the average water supply and the average outflow in a nonzero circle as p¯ and Q-. They are two important efficiency indices.

Typical multiple stationary states.

Efforts for theoretical analysis can be done after sufficient simplification of (1)–(3). For example, the Fourier series of the irrigation intensity of the limit cycle can be computed as follows (assume that i(0-)=-1 and i(0)=1): (5)p(t)=pη+2Σn=1sin(nπη)cos(nπη-2nπηtT), where T is the time length of a period, and η is the proportion of irrigating time length in a period. Also restricting h<H,v<μ1, for all h,v, the first-order derivative of the so-called Poincaré map  at (w0,v) is (6)f(v)|(w0,v)=w1α(w1-w0)β0+β1δ1+δ2+w1α(p-αv)2(β0+β1δ1+δ2)2(αw0v-μ1p). This map indicates the iterative mapping of vegetation density after a loop of irrigation-nonirrigation. It transforms the different equations into the iterative equations and may be convenient to solve some problems; for example, the stationary state of the differential equations corresponds to the equilibrium of f(v). But generally, solving f(v) without any constraint is quite difficult and beyond the scope of this study. Numerical experiment is then carried out to explore the potential nontrivial dynamic behaviors.

3. Analysis of Nontrivial Dynamic Behaviors through Numerical Experiments

Generally speaking, for a given research area and particular vegetation, c,H,α,β,δi,i=1,2,3,4,μ1,μ2 are constants, and p,w0,w1 are artificial parameters. Furthermore, the adjustment range of w0 and w1 is quite small for they should be near to the wilting point and field capacity, respectively, in the irrigation application. Of course, they can also be extended when we move beyond the irrigation problem, see also the discussion in the conclusion section. During the following discussion, we choose α=1(cm·n-1·d-1),c=0.5(d-1),  H=5(cm),β0=1(d-1),β1=0.5(d-1),δ1=δ3=0.5(cm·n-1),δ2=δ4=1(cm·n-1),μ1=2(cm·n-1), and μ2=5(cm·n-1) as a numerical testing environment, tune w0(cm) and w1(cm) coarsely, and consider the dynamic behaviors as depending on fine tuning of p(cm/d). (Here n stands for the quantity of grasses.)

There are three stationary states in Figure 3. Here (w0,w1,p)=(1,7,4). The two thick “circles” are the two stable stationary states, and the thin line divide the first quadrant into two domains of attraction. Points starting from the upper region converge to the upper thick circle, while points starting from the lower region converge to the zero circle. Furthermore, the thin line between [w0,w1] is also a half part of a stationary state (from h=w0 to h=w1), and it is unstable, that is, any tiny disturbance would lead the trajectory to leave this stationary state.

Figure 4 shows some other types of the nonzero stable states for different parameters. The hysteresis is similar with the famous LotkaVolterra predator-prey models raised in the First World War ; that is, during different regions, vegetation intensity and soil moisture content can either be of positive correlation or negative correlation, similar to the relationship between cartilaginous fishes and food fishes in the Mediterranean.

Different types of multiple stationary states.

( w 0 , w 1 , p ) = ( 1,4 , 4 )

( w 0 , w 1 , p ) = ( 1,7 , 3 )

Also for some arbitrary parameters w0,w1,p, the nonzero stable circle vanishes, which means all plants will die after sufficiently long time. The alteration of the number of stationary states is identified as typical bifurcation. As other parameters fixed, the bifurcation map with p and v¯ can be drawn as in Figure 5. Bifurcations with other variables are omitted for conciseness. From the bifurcation map, we can see, as p increased, that the steady vegetation coverage vanishes at p4.4 (Figure 5(a)) and p8.3 (Figure 5(b)), respectively. That means overlarge irrigation intensity will lead plants to die out.

Bifurcation map with v¯ and p. Here (w0,w1)=(1,7) in (a) and (w0,w1)=(1,4) in ((b), (c)). (c) is a locally magnified (b).

Another important dynamic phenomenon presented in Figure 5 is the mutation, or so-called jump. The long-term behavior of our ecohydrological system is continuous with the parameter p for most of the time but jumps at some individual points. For example, in Figure 5(a), v¯ varies continually with p for most time and jumps near p=3.7,3.9,4.4. So along with the continuous increase of irrigation intensity, the limit vegetation density may suddenly go up or drop down discretely. A rough interpretation of this phenomenon is the mutation parameters collide the border of continuity intervals of the vegetation growth Equations (2) and (3) .

Further experiments are worthy to be done to verify the reliability of numerical analysis. In nature, the independent evolution events with the same parameters will never be repeated, and therefore these dynamical phenomena could not be directly verified by field observations. Nevertheless, the independent repeated artificial experiments could be conducted for the demonstration of these phenomena including the automatic irrigation experiments.

Analysis on the nontrivial dynamic behaviors can be used to guide human practice. For example, under the text condition, we plot the bifurcation maps v¯,p¯,Q¯ (stand for the average vegetation coverage, irrigation quantity, and lateral flux, resp.) with p in a same figure, see Figure 6. From the computation we can see that the maximal stable vegetation coverage is at p3.8, and the unit water consumption is also the least economical at this time.

Bifurcation map of v¯,p¯,  and  Q¯ with p.

Another typical example in which dynamic research plays a key role in practical problems is to qualitative study the influence of catastrophic mutation. For example, a fire can push the system from the attracting basin of the nonzero circle to the attracting basin of the zero circle, such as Figure 7.

Influence of catastrophic mutation—a fire case.

(a) The chaotic phase diagram; (b) The Lyapunov exponent. Here (w0,w1,p)=(1,4,8.2); (c) the irrigation time series.

Moreover, deterministic chaos also occurs in this model, for example, when p[7.8,8.3], as shown in Figures 5(b) and 5(c). The chaotic limit set in h-v space when p=8.2 is shown in Figure 8(a). In the scenario, with the constant irrigation strategy, the vegetation would not converge to a unique circle and the vegetation is “fluctuating” forever. Meanwhile, some chaotic characteristic can be computed. Consider the Poincaré map at the point (w0,v), the Lyapunov exponent  of the Poincaré map is larger than zero, as shown in Figure 8(b), which confirms the existence of chaos. The time series of irrigation can be drawn as in Figure 8(c). To be noted, the graph of irrigation time series indicates a quasiperiodic rectangular pulsing function (see Section 2.1). As mentioned above, planar chaos is a specific behavior of nonsmooth dynamic system, which maybe more natural to depict hydrology process. The possible cause of this behavior is that the observation of soil water content of automatic irrigation system is discrete. The scenario should be designed (or avoided) to depend on the purpose of the irrigation system in the practice.

4. Conclusion

An ideal ecohydrological model convenient for dynamic analysis was established in this paper. Compared with traditional models, the nonsmooth property is emphasized. Pulsing function was used to describe the rainfall intensity with a practical interpretation of automatic irrigation. Under this approach, the deterministic process-based ecohydrological model was analyzed theoretically. Based on the model, some deterministic nontrivial dynamic phenomena are revealed. For example, the existence of multiple stationary states, mutation, bifurcation, and chaos was discussed under numerical simulation, with an overview on their possible application in optimization, avoiding the influence of catastrophic mutation, and so on.

This work can be seen as a kind of demonstration to show the existence of some complicated dynamic phenomena in ecohydrological process. Further observational study should be done to determine the scope of parameters in the real world which has nontrivial dynamic phenomena, probably starting from hyperarid area just as Liu et al.  have done to validate ecohydrological feedbacks in hyperarid Tarim River. Noticing that for the convenience of theoretical study, we used a step function to indicate the atmospheric forcing term, instead of considering the stochastic features of rainfall. Actually, one can change the switching policy (4) to a more complicated one, for example, through fixing the time length of i(t)=1 corresponding to a timing irrigation, or randomizing the time length of i(t)=1 corresponding to a stochastic irrigation. A bold hypothesis is that if we understand the rainfall mechanism thoroughly, a proper switching policy can be invented to simulate the fully stochastic features of real rainfall process by introducing additional oceanic and climatic variables, as implicated by the simple formulation of local recycling rainfall in [12, 13]. Also, the model could be further improved by refining the existing feedbacks between hydrological and ecological processes and by incorporating more feedbacks among hydrological, atmospheric, and other related processes.

Nomenclature c :

Side outlet coefficient (d-1)

h :

Storage of aquifer (cm)

H :

Height of the side outlet of the tank (cm)

i :

Irrigation state i{±1}

p :

Irrigation intensity (cm/d)

Q :

Side outlet intensity (cm/d)

t :

Time (d)

T :

Time length of a period (d)

v :

Biomass (n) ((n) is short for quantity of grasses)

w 0 , w 1 :

Bounds of storage of aquifer (cm)

α :

Transpiration coefficient (cm/(d·n))

β 0 , β 1 :

Bounds of evolution rate (d-1)

δ 1 ~ 4 :

Distance between kickpoints and equilibria in (3) (n)

μ 1 , μ 2 :

Equilibria of evolution (n)

η :

The proportion of irrigating time length in a period (100%)

x ˙ :

d x / d t , x stands for any variable

x ¯ :

Average x in a nonzero circle, x stands for any variable.

Acknowledgments

This research was funded by the National Science Foundation of China (NSFC51179084, 51190092, 51222901, and 11072274), the foundation of State Key Laboratory of Hydroscience and Engineering of Tsinghua University (2012-KY-03), and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.

Peterson T. J. Argent R. M. Western A. W. Chiew F. H. S. Multiple stable states in hydrological models: an ecohydrological investigation Water Resources Research 2009 45 3 W03406 10.1029/2008WR006886 Schuster H. G. Just W. Deterministic Chaos, An Introduction 2005 Wiley-VCH 10.1002/3527604804 MR2190040 Kumar P. Variability, feedback, and cooperative process dynamics: elements of a unifying hydrologic theory Geophysical Compass 2007 1 6 1338 1360 Laio F. Porporato A. Ridolfi L. Rodriguez-Iturbe I. Plants in water-controlled ecosystems: active role in hydrologie processes and responce to water stress II. Probabilistic soil moisture dynamics Advances in Water Resources 2001 24 7 707 723 2-s2.0-0035398604 10.1016/S0309-1708(01)00005-7 Sivakumar B. Chaos theory in hydrology: important issues and interpretations Journal of Hydrology 2000 227 1–4 1 20 2-s2.0-0034737012 10.1016/S0022-1694(99)00186-9 Sivakumar B. Chaos theory in geophysics: past, present and future Chaos, Solitons and Fractals 2004 19 2 441 462 10.1016/S0960-0779(03)00055-9 MR2007457 ZBL1080.37610 D’Andrea F. Provenzale A. Vautard R. De Noblet-Decoudre N. Hot and cool summers: multiple equilibria of the continental water cycle Geophysical Research Letters 2006 33 24 L24807 10.1029/2006GL027972 Ridolfi L. D’Odorico P. Laio F. Tamea S. Rodriguez-Iturbe I. Effect of vegetation-water table feedbacks on the stability and resilience of plant ecosystems Water Resources Research 2006 42 1 W01201 10.1029/2005WR004444 D’Odorico P. Caylor K. Okin G. S. Scanlon T. M. On soil moisture-vegetation feedbacks and their possible effects on the dynamics of dryland ecosystems Journal of Geophysical Research 2007 112 G4 G04010 10.1029/2006JG000379 Borgogno F. D'Odorico P. Laio F. Ridolfi L. Effect of rainfall interannual variability on the stability and resilience of dryland plant ecosystems Water Resources Research 2007 43 6 2-s2.0-34547227218 10.1029/2006WR005314 W06411 Runyan C. W. D’Odorico P. Ecohydrological feedbacks between salt accumulation and vegetation dynamics: the role of vegetation-groundwater interactions Water Resources Research 2010 46 11 W11561 10.1029/2010WR009464 Rodriguez-Iturbe I. Entekhabi D. Bras R. L. Nonlinear dynamics of soil moisture at climate scales, 1. Stochastic analysis Water Resources Research 1991 27 8 1899 1906 2-s2.0-0026281797 10.1029/91WR01035 Rodriguez-Iturbe I. Entekhabi D. Lee J.-S. Bras R. L. Nonlinear dynamics of soil moisture at climate scales, 2. Chaotic analysis Water Resources Research 1991 27 8 1907 1915 2-s2.0-0026304150 10.1029/91WR01036 Laio F. Tamea S. Ridolfi L. D'Odorico P. Rodriguez-Iturbe I. Ecohydrology of groundwater-dependent ecosystems: 1. Stochastic water table dynamics Water Resources Research 2009 45 5 2-s2.0-67650280449 10.1029/2008WR007292 W05419 Robinson J. S. Sivapalan M. Temporal scales and hydrological regimes: implications for flood frequency scaling Water Resources Research 1997 33 12 2981 2999 2-s2.0-0030659491 Ghorbani M. A. Khatibi R. Sivakumar B. Cobb L. Study of discontinuities in hydrological data using catastrophe theory Hydrological Sciences Journal 2010 55 7 1137 1151 2-s2.0-78650471763 10.1080/02626667.2010.513477 van der Schaft A. Schumacher H. An Introduction to Hybrid Dynamical Systems 2000 251 Berlin, Germany Springer MR1734638 ZBL0988.93060 Bernardo M. Chris J. Alan R. Piotr K. Arne B. Gerard O. Petri T. Numerical simulation of finite dimensional multibody nonsmooth mechanical systems ASME Applied Mechanics Reviews 2002 55 107 150 Matveev S. A. Savkin V. A. Qualitative Theory of Hybrid Dynamical Systems 2000 Berlin, Germany di Bernardo M. Budd C. J. Champneys A. R. Kowalczyk P. Nordmark A. B. Olivar Tost G. Piiroinen P. T. Bifurcations in nonsmooth dynamical systems SIAM Review 2008 50 4 629 701 10.1137/050625060 MR2460800 ZBL1168.34006 Kunze M. Non-Smooth Dynamical Systems 2000 1744 Berlin, Germany Springer 10.1007/BFb0103843 MR1789550 Collins D. B. G. Bras R. L. Climatic and ecological controls of equilibrium drainage density, relief, and channel concavity in dry lands Water Resources Research 2010 46 4 W04508 10.1029/2009WR008615 Lu E. Determining the start, duration, and strength of flood and drought with daily precipitation: rationale Geophysical Research Letters 2009 36 12 L12707 10.1029/2009GL038817 Singh B. Boivin J. Kirkpatrick G. Hum B. Automatic irrigation scheduling system (AISSUM): principles and applications Journal of Irrigation & Drainage Engineering—ASCE 1995 121 1 43 56 2-s2.0-0028972093 Abraham N. Hema P. S. Saritha E. K. Subramannian S. Irrigation automation based on soil electrical conductivity and leaf temperature Agricultural Water Management 2000 45 2 145 157 2-s2.0-0034064441 10.1016/S0378-3774(99)00078-5 Yildirim M. The performance of irrigation automation system based on soil moisture level for pepper (Capsicum annuum L.) growth Journal of Food, Agriculture and Environment 2010 8 3-4 629 634 2-s2.0-78751624955 Cardenas-Lailhacar B. Dukes M. D. Miller G. L. Sensor-Based automation of irrigation on bermudagrass, during wet weather conditions Journal of Irrigation and Drainage Engineering 2008 134 2 120 128 2-s2.0-41749085776 10.1061/(ASCE)0733-9437(2008)134:2(120) Ruisen Z. Xinguang D. Yingjie M. Sustainable water saving: new concept of modern agricultural water saving, starting from development of Xinjiang's agricultural irrigation over the last 50 years Irrigation and Drainage 2009 58 4 383 392 2-s2.0-70349983281 10.1002/ird.414 Sugawara M. Automatic calibration of the tank model Hydrological Sciences Bulletin 1979 24 3 375 388 2-s2.0-0018517222 Chui T. F. M. Low S. Y. Liong S.-Y. An ecohydrological model for studying groundwater—vegetation interactions in wetlands Journal of Hydrology 2011 409 1-2 291 304 10.1016/j.jhydrol.2011.08.039 Allen R. G. Pereira L. S. Raes D. Smith M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements 1998 56 Irrigation and Drainage Burghes D. N. Borrie M. S. Modelling with Differential Equations 1981 Chichester, UK Horwood MR616822 Liu D. Tian F. Hu H. Lin M. Cong Z. Ecohydrological evolution model on riparian vegetation in hyper-arid regions and its validation in the lower reach of Tarim River Hydrological Processes 2012 26 13 2049 2060 10.1002/hyp.8313