# The Stochastic Dynamics for Ecological Tourism System with Visitor Educational Intervention

Authors: Dongping Wei; Shouwen Wen; Yizeng Chen; Yong Xu; Bing Liang
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The ever-increasing visitation in parks and protected areas continues to present a considerable challenge for worldwide land managers with allowing recreational use while preserving natural conditions. In China, the fast expanding visitation in protected areas is quickly damaging the natural resources and precious culture without effective visitor education, while regulation and site management are also gaining very limited efficacy. We propose a differential equation to describe the ecological tourism system. Shown by the theoretical proof and numerical simulation, the ecological tourism system is unstable without any perturbed factors, especially visitor educational intervention, because the solution of the dynamic system explodes in a finite time given any initial value. Supposing that the intrinsic increasing rate of stakeholders in the systems stochastically perturbed by the visitor educational intervention, we discover that the stochastic dynamic model can effectively suppress the explosion of the solution. As such, we demonstrate that the tourism system can develop steadily and safely even under a large amount of visitors in public vacation, when employing continuous visitor education intervention programmes.


---

## Body

## 1. Introduction

Currently, the Chinese tourism market is the primary domestic tourism market owning about 1.7 billion visitors per year and the fourth inbound tourism market and the primary outbound tourism market in the Asian-Pacific region owning about 45.84 million visitors per year [1]. Many scholars proclaim in surprise that the era of recreation is approaching ordinary Chinese daily life. In the public holidays such as golden week of Labor’s day and National day, it is crowded and chaotic in every popular scene site such as Huangshan Mountain, Huashan Mountain, and the Great Wall, In these crowded parks or natural protected areas, the negative impacts, such as litter, tree damage, noise, and the rude behaviors of visitors, shock the environmentalists. In the United States, scholars concluded that the undesirable impacts on the tourism destination with fast expanding visitation include resource degradation and social impacts [2, 3]. The resource impacts associated with large amount of visitors include the trampled vegetation and soil damage, soil compaction and erosion, litter, human waste, and wildlife disturbance. The social impacts include the diminished satisfaction of visitors and recreation conflicts between visitors, visitors and park managers, visitors and neighbors of tourism destinations, and the culture conflicts [4]. Ferreira and Rosso [5] discovered that the Sao Paulo coast is highly depredated by the increased tourist activity. Manning and Aderson [6] also found that the dramatically grown visitors’ use of Arches National Park has had several important impacts on the park including trampling of fragile soils and vegetation (impacts on soil and vegetation) and crowding on trails and at attraction sites (crowding, impacts on trails and attraction sites). Left unmanaged, these impacts can lead to unacceptable changes in resource and social conditions such as the loss of sensitive or rare plants and animals, or declines in visitor satisfaction [7]. In China, many visitors feel awful when they arrive the tourism destination because of crowded people and the rude behaviors.Facing the increasingly serious negative tourism impacts, many scholars propose different strategies to reduce the impacts. Ferreira and Rosso [5] proposed that the management strategies should include isolation of sensitive areas, construction of boardwalks, visitor education, and monitoring programmes to reduce the degradation of rocky shore fauna on the Sao Paulo coast caused by the development of tourism. In United States, many park managers commonly employ regulations, site management, and visitor education to address the resource and social impacts. The regulations directly curtail tourism freedom and antagonize visitors by their enforcement. But the enforcements of regulations are often costly and lack efficacy. Peterson and Lime [8], Hendee and Dawson [3], and Stankey and Schreyer [9] noted that site management actions were less direct but they were also costly and permanently altered the natural setting and the nature of recreational experiences. The visitor education is a more appropriate, light-handed and indirect management response to reduce resource and social impacts [4, 10]. Docutte and Cole [11] presented that visitor education has been widely touted as the most appropriate approach to managing recreation in wilderness. Visitor education is also easier to win the support of visitors [3]. Hence, in the United States, the federal land management agencies employed a Leave No Trace (LNT) programme for visitor education in order to reduce the resource and social impacts. And many other international organizations such as Members of the World Tourism Organization developed a Global Code of Ethics for Tourism in 1999 to help minimize the negative impact of tourism on the environment and the cultural heritage (WTO, 2006). The Ecotourism society (1993) has also put forward “Ecotourism Guidelines for Nature Operators” to guide the visitor educational efforts to minimize the social and environment impacts. Hendee and Dawson [3] concluded that effective educational programmes can (1) reduce per capita and cumulative resource and experiential impacts; (2) aid in keeping resource and social conditions within acceptable limits; (3) reduce the need for site management on regulatory actions such as limits; (4) provide a key component of sustainable tourism programmes.There are a lot of references to study the visitor educational message content, message delivery approaches, audience characteristics and theoretical grounding including moral development theory and decisionmaking theory, reasoned actions and planned behavior theory. However, there is little research to study the relationship between the stakeholders and the human ecological system of tourist destinations. In fact, the stakeholders in ecotourism destinations cannot survive without the support from one another, and they must unite to a population in a certain space [12]. In 2012, Shou-wen  started to analyze the dynamic property of visitors and the ecological tourism system under the visitor educational intervention. Shou-wen et al. [13] set up the stochastic increasing model of stakeholders perturbed by the continuous visitor education. In this paper, we make efforts to the further study to better understand the relationship between the stakeholders and ecological tourism system intervened by visitor educational programmes. Based on stochastic differential equation theory [14], we explore whether the ecological tourism system can develop sustainably when perturbed by the continual visitor educational intervention.In order to clearly present the mathematical models in the next three sections, we introduce the key definitions and notations of variables in the following.x(t):
stakeholders in an ecological tourism system at timet.x0:
the initial value of differential equation at timet0.α:
the intrinsic increasing rate of stakeholders within the ecological tourism system.β:
the immigration rate of stakeholders from outside the system.γ:
the emigration rate of stakeholders from the inner system.K:
the maximum carrying capacity of the ecological tourism system.ω(t):
the Brownian movement.ε:
the intensity of visitor education intervention.τe:
the explosion time of stochastic differential equation.τk:
the stopping time of stochastic differential equation.k0:
a sufficiently large value.k:
any nonnegative integer.V(x):
aC2-function, V:R+→R+.M:
a sufficient large positive value.T:
a sufficient large positive constant.


## 2. The Logical Differential Equation Model of Ecological Tourism System

In the literature, the ecological tourism system is defined as all the natural resources including air, rocks, mountains, soils, trees, wildlife, and all the stakeholders in visitation destination. And the stakeholders [12] in ecotourism include the seven factors: (1) the tourism resource suppliers’ group such as tourism resource managers or park managers; (2) the visitation product sellers group including tourism company, and tourism website managers; (3) the group of relative products or services suppliers who are providing transportation and hotels, catering, other tourism memories products; (4) the group of local communities including the staff of tourism companies and workers in other company and local people; (5) the group of governments including super government organizations and centeral and local governments; (6) the special interest group including nongovernment organizations like environmental organizations, human rights and labor rights organizations, or trust and environmental charity institutes, academics, and social medias; and (7) the core tourism products and services consumers like visitors.We supposed that there arex(t) stakeholders in an ecological tourism system at time t, considering that there are some stakeholders entering the system from outside the ecological system and some others leave the system when they finish their activities or business. We denote α as the intrinsic increasing rate of stakeholders within the ecological tourism system, β as the immigration rate of stakeholders from outside the system, and γ as the emigration rate of stakeholders from the inner system. K stands for the maximum carrying capacity of the ecological tourism system. And then 1/K represents the recreation impact factors per capita. So x(t)/K means the recreation impact factor at time t. So mapping the logistical differential equation on the ecological system, we have
(1)dxdt=x(α+β-γ-1Kx).
It is the original differential equation to describe the relationship between the stakeholders and the ecological system.Theorem 1.
For any given initial valuex0, the solution of (1) x(t) cannot avoid an explosion in a finite time.Proof.
Since the variablex(t) denotes the population of stakeholder,  x(t) must be a positive value at time t≥0. When the parameters of (1) α,β,γ∈(0,1), K∈R+, (1) only has a local solution as follows
(2)x(t)=α+β-γ-(1/K)+e-(α+β-γ)t(α+β-γ-(1/K)x0)/x0.
Given that x0>0, it is not difficult to concludethat  x(t) explodesto infinity when t→𝒯,
(3)𝒯=-1α+β-γlog-(1/K)x0(α+β-γ-(1/K)x0).
Hence, the solution of (1) x(t) will explode to infinity in a finite time.In other words, the original differential equation model just has the local solution. The simulation data in Table1 and Figure 1 show that the dynamic model (1) explodes in a finite time for any initial value x0=30,40,50,70,80,100, when α=0.25, β=0.12, γ=0.08. and K=1000. For instance, according to simulation data in Table 1 and Figure 1, we can easily find that the number of stakeholders x(10)=1081 is a little larger than the maximum capacity of the system K=1000, but at time t=11, x(11)=2990 is much larger than K=1000, when the initial value of stakeholder x0=40 without any control on the system. Hence, the system is chaotic and unsustainable after t=10. Eventually, the system is collapsed at time t=12, because x(12)=-5679<0. When x0=30,40,50,70,80,100, the system is collapsed at time t=14,12,10,7,6,4, see Figure 1; therefore, it can be concluded that the system is collapsed faster with larger initial value x0.Table 1
The simulation ofx(t) given the different initial value x0.
t
x
(
t
)
α
=
0.25,  β=0.12,  γ=0.08,  K=1000
x
0
=
30
x
0
=
40
x
0
=
50
x
0
=
70
x
0
=
80
x
0
=
100
1
53
79
112
215
301
688
2
63
95
138
280
415
1262
3
75
116
172
382
617
4545
4
91
143
219
555
1069
−3635
5
110
179
287
916
2907
−1422
6
135
229
392
2087
−5956
−931
7
169
302
575
−19651
−1638
−717
8
215
417
961
−1963
−1006
−598
9
281
622
2303
−1104
−755
−523
10
382
1081
−11137
−800
−621
−472
11
556
2990
−1841
−647
−538
−435
12
918
−5679
−1069
−555
−482
−408
13
2098
−1619
−785
−494
−443
−386
14
−18836
−1000
−638
−451
−413
−370
15
−1956
−752
−549
−420
−391
−357
16
−1102
−619
−490
−396
−373
−346
17
−800
−537
−448
−377
−359
−337
18
−646
−481
−417
−362
−348
−330
19
−554
−442
−394
−351
−339
−324Figure 1
The simulation ofx(t) given the different initial values of x0.Similarly, Figure2 demonstrates that the explosion of the dynamic model (1) cannot avoid any initial value x0=30,40,45,60,70,80, given that α=0.15, β=0.08, γ=0.06, and K=1000. For instance, Table 2 shows that the number of stakeholders x(16)=1698 is much larger than the maximum capacity of K=1000, given the initial value of stakeholder x0=30 without any control on the system. Hence, the system is chaotic and unsustainable after time t=16. Eventually, the system will be collapsed at time t=18. When x0=30,40,45,60,70,80, the system is collapsed at time t=18,13,12,10,6,3, see Figure 1. Therefore, it can be concluded that the system is collapsed faster with larger initial value x0. Comparing with Figure 1, Figure 2 indicates that the system is collapsed slower with smaller parameters α, β, and γ.Table 2
The simulation ofx(t) given the different initial value x0 and α,  β,  and  γ.
t
x
(
t
)
α
=
0.15,  β=0.08,  γ=0.06,  K=1000
x
0
=
30
x
0
=
40
x
0
=
45
x
0
=
60
x
0
=
70
x
0
=
80
1
58
98
127
166
311
828
2
65
112
148
198
404
1588
3
73
129
174
240
558
10110
4
83
151
208
299
856
−2574
5
94
179
254
385
1685
−1196
6
108
214
319
525
15181
−802
7
125
263
417
787
−2388
−615
8
145
332
581
1456
−1157
−507
9
171
437
908
6590
−786
−437
10
204
616
1878
−2942
−607
−387
11
249
989
102870
−1263
−502
−351
12
312
2224
−2126
−829
−433
−323
13
405
−15070
−1097
−630
−384
−301
14
559
−1851
−760
−516
−349
−283
15
860
−1025
−592
−443
−321
−268
16
1698
−727
−493
−392
−299
−256
17
16191
−574
−427
−354
−282
−246
18
−2366
−481
−380
−325
−268
−237
19
−1153
−418
−345
−303
−256
−230Figure 2
The simulation ofx(t) given the different initial values of x0 and α, β, and γ.The visitorsx0 are very large in some protected areas especially on the public vacations in China. So the solution of the differential dynamic equation will inevitably explode without any interventional programme. The tourism system of some famous protected areas or precious fragile natural history park is currently in the great danger of large visitations.


## 3. The Stochastically Dynamics Model Perturbed by Visitor Educational Intervention

As demonstrated in the simulations of Section2, with the high and disordered consumption, the ecological tourism system especially precious fragile protected areas or parks will eventually collapse without any visitor education or other visitor interventions. The booming China tourism market is extremely crowded in many famous tourism destinations on the public vacations. Although there are a lot of visitor regulations and rules in the tourism destinations, the lack of the effective visitors educational intervention coused the natural resource and culture resource to be badly degraded. The natural resource impacts and social impacts are extremely serious. So it is critical importance to develop effective and systematical visitor educational programme such as the Chinese version of LNT Programme to avoid systematical explosion of the ecological tourism system. Increasing tourist activity in protected regions demands management strategies including visitor education to reduce the impacts on the ecological tourism system and avoid the explosion of solution of the dynamic equation describing the system. In this paper, we suppose that the intrinsic increasing rate of stakeholders is intervened by visitor education, but the efficacy of visitor education is different from diverse visitor characteristics, the educational message content, message delivery approaches, and so forth. Hence, the intrinsic increasing rate of stakeholders behaves stochastically when perturbed by visitor intervention. Naturally, we can suppose that the intrinsic rate of increase is stochastically perturbed with
(4)α⟶α+εω(t),
where ω(t) and ε represent the Brownian movement and the intensity of visitor education, respectively. Substituting (4) into the original differential equation (1), we can obtain the stochastic dynamic model of ecological tourism system perturbed by visitor education intervention programme,
(5)dx=x(α+β-γ-1Kx)dt+εx2dω(t).
The critical purpose of the paper is to study whether the stochastic parameter of intrinsic increasing rate of stakeholders can effectively suppress the explosion of the solution of (2) in a finite time with any initial value x0 and stop the collapse of the ecological tourism system eventually. In this paper, we should prove that the solution of (5) can no longer explode in a finite time with a probability one when α+β<γ+(1/K)x0. And we will then see that for arbitrary parameters α,β,γ∈[0,1] and K∈R+, the dynamic system will not explode in a finite time with a probability one provided any intensity of the perturbed visitor educational intervention ε. The results demonstrate that effective visitor education programmes will support the sustainable development of ecological tourism system and suppress the explosion of the stochastic dynamic model.


## 4. The Existence of Nonexplosion Solutions of the Ecological Tourism Dynamic System with Visitor Education of Intervention

In the paper, we denote(Ω,ℱ,{ℱt}t≥0,𝒫) to be a complete probability space with a filtration {ℱt}t≥0 satisfying the usual condition (i.e., right continuous, and ℱ0 contains all p-null sets) [14]. Let ω(t) be one-dimensional Brownian motion defined on the probability space.The definitionx(t) should be nonnegative because it represents the stakeholders. For any given initial value x0, if the stochastic differential equation (5) has a unique global solution, then the coefficient of the equation should be satisfy the linear growth condition and local Lipschilz condition [14–17]. Unfortunately, the coefficients of (5) only satisfy the locally Lipschitz continuous condition, so the solution of (5) only has local solution which may explode at a finite time [18].Now, we will show that the visitor education intervention can suppress the explosion of the solution as proved by the following theorem.Theorem 2.
Given any system parametersα,β,γ∈[0,1], K∈R+ and any initial value x0∈R+ if intensity of visitor education intervention ε≠0, there is a unique solution x(t) of (3) at t≥0, and the solution will be positive with probability one, namely, x(t)>0 when t≥0 a.s.Proof.
For the coefficients of (5) are the locally Lipschitz continuous, for any given initial value x0∈R+, there is a unique solution x(t) on the t∈(0,τe) where τe is the explosion time [18, 19]. To prove the solution of (5) x(t) is global and positive, we need to show that τe=∞ almost surely. Let k0>0 be sufficiently large value and the initial value of the equation x0 belong to the interval [1/k0,k0]. For any nonnegative integer k≥k0, we define the stopping time as follows
(6)τk=inf{t∈[0,τe):x(t)∉(1k,k)},
where setting inf∅=∞ throughout the whole paper, as usual ∅ denotes the empty set. Obviously, τk is increasing when k→∞. Set τ∞=limk→∞τk hence τ∞≤τe a.s. If it can be proved that τ∞=∞ a.s., then it is natural that τe=∞ a.s. and x(t)∈R+ for any t≥0. In other words, in order to show that x(t)∈R+ will not explode, we should show that τ∞=∞ a.s. first. Employing contradiction methodology, if this statement τ∞=∞ a.s. is false, then there is a pair of constant T≥0 and ε∈(0,1), such that
(7)P{τ∞≤T}>ε.
Therefore, an integer k1≥k0 can be found, such that
(8)P{τ∞≤T}>ε∀k≥k1.
Define a C2-function V:R+→R+,
(9)V(x)=x-1-logx0.5.
It is easily proved that C2-function V(x) is nonnegative onx>0.
Ifx(t)∈R+, it can be calculated from the Ito formula as follows
(10)dV(x)={0.5(x-0.5-x-1)x[(α+β-γ-1Kx)dt+εxdω(t)]}+0.5(0.5x-2-0.25x-1.5)ε2x4dt={0.5(x0.5-1)[(α+β-γ-1Kx)+(0.25-0.125x0.5)ε2x2(α+β-γ-1Kx)]}dt+0.5(x0.5-1)εxdω(t),
where we denote x(t)=x for simplification.
Then (10) can be simplified as follows
(11)0.5(x0.5-1)(α+β-γ-1Kx)=0.5(α+β-γ)x0.5-0.51Kx1.5-0.5(α+β-γ)+0.51Kx=-0.51Kx1.5+0.51Kx+0.5(α+β-γ)x0.5-0.5(α+β-γ),(0.25-0.125x0.5)ε2x2=0.25ε2x2-0.125ε2x2.5.
We note that
(12)F(x)=-0.125ε2x2.5+14ε2x2-0.51Kx1.5+0.51Kx+0.5(α+β-γ)x0.5-0.5(α+β-γ)≤-0.125ε2x2.5+14ε2x2+0.51Kx1.5+0.51Kx+0.5(α+β-γ)x0.5+0.5(α+β-γ).
Then we have dV(x)=F(x)dt+0.5(x0.5-1)εxdω(t).
Ifε≠0, then it is not difficult to yield that F(x) is bounded by sufficiently large positive value M, so that dV(x)≤M+0.5(x1.5-x)εdωT(t).
It is clear to obtain the following expression by integrating the both sides of the proceeding inequality from 0 toT∧τk as follows
(13)∫0T∧τkdV(x)≤∫0T∧τkMdt+∫0T∧τk0.5(x0.5-1)εxdω(t).
With (T∧τk)∈R+, taking expectations on both sides of inequality (8), we can yield
(14)EV(x(T∧τk))≤V(x0)+ME(T∧τk)≤V(x0)+MT.
Let Ω𝓀={τ𝓀≤T} for 𝓀≥k and the inequality (7) so that P(Ω𝓀)>ε. Taking any w∈Ω𝓀, with the definition of stopping time (6), x(τk,w) equals either 1/k ork, and then it is easy to conclude that V(x(τk,w)) is not less than either k-1-0.5logk or 1/k-1-0.5log(1/k). So x(t)≥k or x(t)≤1/k when t≥τk. Therefore, it is not difficult to get the following inequality:
(15)V(x(τk,w))≥[k-1-0.5logk]∧1k-1-0.5log(1k).
Taking expectation on both sides of (15), we get
(16)EV(x(τk,w))≥[k-1-0.5logk]∧1k-1-0.5log(1k).
It is easy to show that
(17)V(x0)+MT≥E[1Ω𝓀(w)V(x(T∧τk))]≥ε[k-1-0.5logk]∧1k-1-0.5log(1k),
where 1Ω𝓀(w) is the indicator function of Ω𝓀. Letting k→∞ in (17), then
(18)∞>V(x0)+MT=∞,
where T is a constant denoted in (7). Obviously, (18) is contradicted. So we can conclude that the statement of (7) is false. Hence, there are constants T>0 and ε∈(0,1) such that P{τ∞≤T}>ε. Consequently, the original statement is true in which P{τk≤τ∞}<ε. In other words, the statement τ∞=∞ is true almost surely.It is clear to derive from the Theorem2 that there is a positive solution x(t) of (7) and x(t) will never explode in a finite time. Following the condition of Theorem 2, we can deduce that the visitor educational intervention brings the stochastic intrinsic increasing rate of stakeholders into the tourism ecological system. When the intensity of intervention ε≠0, the stochastic intrinsic increasing rate can effectively suppress the explosion of the solution x(t) for any given x0 and avoid the collapse mentioned in Section 2. It means that the slight or tiny stochastical change of intrinsic increasing rate of stakeholder introduced by visitor educational intervention can keep the tourism ecological system developing safely and steadily. In fact, early in 1970s, Fishbein, Fazio [20], and Bradley [21] already noticed that the negative recreation impacts on wildlands and wildlife, and the ecological system or environment could not keep sustainable development without recreation planning and management. According to Theorem 2, the visitor educational intervention is a necessary condition for the safe and stable developing of the tourism ecological system. Especially in China, the booming recreation market is in urgent need of the visitor educational intervention to stop the chaotic tourism activities.


## 5. The Boundedness of the Solution of the Dynamic Model

With the intensities of the visitor educational intervention, the solution of the dynamic model (5) cannot explode in a finite time following by the Theorem 2. However, it is not enough to show that the solution exists and does not explode. We consider how the solution of the dynamic model varies in a finite interval in this section. In other words, we will show that the solution of (5) is ultimately bounded.Denotex(t,x0) as the unique global solution of (5) for any given initial value x0>0. We define stochastically ultimate bounded property of the solution of the dynamic model as follows.Definition 3.
For any initial valuex0>0, if the solution of (3) x(t,x0) has the property
(19)limt→∞infx(t,x0)>0,
then (5) is defined to be stochastically ultimately bounded with probability one. The following theorem shows (5) is stochastically ultimately bounded.Theorem 4.
If the coefficient of (5) satisfies the following inequality:
(20)12ε2-1K<0,(α+β-γ+1K)2+4(12ε2-1K)(α+β-γ)≤0.
Then the solution of (5) x(t,x0) is stochastically ultimately bounded with a probability one. Given any initial value x0∈R+, the solution of (3) has the property
(21)limt→∞supx(t,x0)<∞.
So the solution of (5) x(t,x0) belongs to a finite interval
(22)[limt→∞infx(t),limt→∞supx(t)].Proof.
It is necessary to define aC2-function V:R+→R+ by
(23)V(x)=x-1-logx.
Applying the Ito formula, we derive from (23) that
(24)dV(x)=LV(x)dt+(1-x-1)x2εdω(t),
where
(25)LV(x)=(1-x-1)x(α+β-γ-1Kx)+12x-2ε2x4=(x-1)(α+β-γ-1Kx)+12ε2x2=(12ε2-1K)x2+(α+β-γ+1K)x-(α+β-γ).
It is obvious to see that LV(x(t)) is a quadratic function. Applying the property of quadratic function, when
(26)12ε2-1K<0,(27)Δ=(α+β-γ+1K)2+4(12ε2-1K)(α+β-γ)≤0.
If (20) is satisfied, the quadratic function LV(x(t)) must be less than zero. Therefore, we can get that
(28)dV(x(t))≤0dt+(x2-x)εdω(t).
Integrating both sides of inequality (28) from 0 to t, we can get
(29)∫0tdV(x(t))≤∫0t0dt+∫0t(x2-x)εdω(t).
Thus,
(30)V(x(t))≤V(x0)+∫0t(x2-x)εdω(t).
For convenience, we define
(31)X(t)≔V(x0)+∫0t(x2-x)εdω(t).
Since V(x) is a nonnegative function, applying Doob’s martingale convergence theorem, we can yield
(32)limt→∞x(t)<∞a.s.
Thus, it is natural to yield limt→∞supx(t)<∞ a.s.
Recalling the definition ofC2-function V:R+→R+, if and only if x→∞ or x→0, the limit of V(x)=x-1-logx is ∞. Hence limt→∞supV(x(t))=limt→∞sup(x(t))-1-log(x(t))<∞ a.s.
Then it is not difficult to deduce that(33)0<limt→∞infx(t)<limt→∞supx(t)<∞a.s.
In conclusion, given any initial value x0, the solution of (3) x(t) will vary in the interval of [limt→∞infx(t),limt→∞supx(t)].From Theorem4, the stakeholders of the ecological tourism system with visitor educational intervention can develop within an interval sustainably and steadily. Different kinds of stakeholders retained in a safety interval can immigrate from outside the system or emigrate from inside the system with regular rate. Hence, visitor education programmes are important elements of management efforts to ensure the protection of protected areas including Natural History Park and precious culture resource areas, and so forth.


## 6. The Simulations of the Stochastic Dynamic Model of Ecological Tourism System

Following the Euler-Maruyama method [22] and using Matlab program, we can obtain the numerical solutions and their figures of the stochastic dynamic model of ecological tourism system (5) with different parameters α, β, γ and different initial values x0.Givenα=0.25, β=0.12, γ=0.08, and K=1000, it is easy to gain the simulation data of x(t) and draw Figure 3 when x0=30 and x0=100. Figure 3 demonstrates that the stochastic dynamic model is varying far away from the maximum carrying capacity of the ecological system K=1000. And the peak value of x(t) is 625 much smaller than K=1000. Comparing Figure 1 with Figure 3, it can be concluded that the stochastic dynamic model (5) can effectively suppress the explosion and avoid the chaotic and unsustainable situation with the same parameters α, β, γ.Figure 3
The simulation ofx(t) with α=0.25, β=0.12, γ=0.08, and K=1000 and different initial values x0.Givenα=0.15, β=0.08, γ=0.06, K=1000, it is not difficult to get the simulation data of x(t) and draw Figure 4 when x0=30 and x0=80. Figure 2 shows that the original logistical differential Equation (1) collapses in a finite time. But Figure 4 shows that the stochastic dynamic model (5) does not explode with the same parameters α, β, γ. The peak value of x(t) is 298 with smaller parameters α, β, γ, see Figure 4.Figure 4
The simulation ofx(t) with α=0.15, β=0.08, γ=0.06, K=1000 and different initial value x0.Shown by Figure5, the ecological tourism system perturbed by visitor educational intervention does not collapse even with larger initial value x0=900 which is very close to K=1000. When there is a large number of visitors in the ecological tourism system, the effective visitor education programme will also keep resource and social conditions within acceptable limits by employing low impact behaviors. Figure 5 demonstrates that the system restores to normal condition quickly. Following Figures 3 to 5, it can be concluded that the ecological tourism system with visitor educational intervention can develop safely and steadily.Figure 5
The simulation ofx(t) with α=0.25, β=0.12, γ=0.08, K=100 and x0=900.


## 7. Conclusions

The ever-increasing demand for outdoor recreation in parks and protected areas worldwide continues to present a considerable challenge for managers [22]. Allowing recreational use while preserving natural conditions, it also presents a challenge for park managers in China because of fast expanding visitation. The protected areas including some natural history parks will be quickly damaged without visitor education intervention, while regulation and site management are also gaining very limited efficacy.In Section2, supposing that the different kinds of stakeholders of ecological tourism system grow with the specified intrinsic increasing rate and the stakeholders from outside competitors can immigrate into the system and some insider of the system can also go out freely when they finish their activities or business, we set up a differential equation to describe the ecological tourism system. Without any perturbed factors especially visitor educational intervention shown by the theoretical proof and numerical simulation, the solution of the dynamic system will explode in a finite time given any initial value.In order to avoid the explosion of the solution of the dynamic model and the unstable situation of the ecological tourism system, we add visitor educational intervention into the system in Section3. Supposing that the intrinsic increasing rate of stakeholders is stochastically perturbed by the visitor educational intervention, we explore a stochastic dynamic model which can effectively suppress the explosion of the solution shown in Section 4. And in Section 5, we also prove that the solution of the stochastic dynamic model will vary in a finite interval. In Section 6, following the Euler-Maruyama method, we obtain the numerical solutions and their figures of the stochastic dynamic model of ecological tourism system. So we demonstrate that the tourism system will develop steadily and safely even under a large amount of visitors in public vocation, when employing continuous visitor education intervention programmes.Generally, the findings of the stochastic dynamic model support that the continuous visitor education intervention including ecological knowledge, minimum impact behavior or technique, and regulations message is the necessary response to the challenge of recreation impact issues which try to manage resources sustainably and provide quality visitor experiences.



---

*Source: 617608-2013-06-10.xml*