#### Abstract

Management of ecological tourism in protected areas faces many challenges, with visitation-related resource degradations and cultural impacts being two of them. To address those issues, several strategies including regulations, site managements, and visitor education programs have been commonly used in China and other countries. This paper presents a multiparameter stochastic differential equation model of an Ecological Tourism System to study how the populations of stakeholders vary in a finite time. The solution of Ordinary Differential Equation of Ecological Tourism System reveals that the system collapses when there is a lack of visitor educational intervention. Hence, the Stochastic Dynamic of Ecological Tourism System is introduced to suppress the explosion of the system. But the simulation results of the Stochastic Dynamic of Ecological Tourism System show that the system is still unstable and chaos in some small time interval. The Multiparameters Stochastic Dynamics of Ecological Tourism System is proposed to improve the performance in this paper. The Multiparameters Stochastic Dynamics of Ecological Tourism System not only suppresses the explosion of the system in a finite time, but also keeps the populations of stakeholders in an acceptable level. In conclusion, the Ecological Tourism System develops steadily and sustainably when land managers employ effective visitor education intervention programs to deal with recreation impacts.

#### 1. Introduction

According to the 2014 Yearbook of China Tourism Statistics, China enjoys the largest domestic tourism market in the world. Its inbound tourism market ranks the fourth and its outbound market ranks the first in the Asian-Pacific region. A large population gathers on tourism destinations in public vocational days including the Golden Weeks during the Labor Day (May 1) and National Day (October 1). Visitors crowd into the protected areas such as Huangshan National Natural Protected Park, Zhangjiajie National Forest Park, and Tai Mountain. These parks are so populated that it is hard to even walk. Due to the lack of visitor educational interventions, the parks suffer from quick unrecoverable degradation of resources. For example, litters are observed everywhere: precious trees are damaged, and wild animals are badly disturbed. The social impacts are no less: tensions among tourists due to crowding in the parks, conflicts between visitors and park managers, and, in an extreme case, tens of thousands of tourists trapped in Huashan Mountain, Shaanxi province, on October 2, 2012. In this case, a huge number of visitors challenged the capacity of the cable cars at Huashan Mountain, leaving tens of thousands stuck at the mountain peak late on Tuesday night (see http://english.cri.cn/6909/2012/10/04/2821s725417.htm). According to China Central Television, restless visitors demanded refunds from the tourism committee, and police was dispatched to help dealing with the crisis. Clearly, the negative resource and social impacts diminish the visitors’ satisfaction quickly.

Hence, the fast expanding visitation in some protected areas severely challenges the protection of the environment and precious natural resources, cultural resources, and social conditions. The United State faced the same situation in the 1970s. To address this issue, the land managers of protected areas, tourism providers, and other stakeholders commonly employed regulations and site management. As a direct management measure, the regulations curtailed visitors’ freedom and emphasized on the potential enforcement with punitive actions [1]. However, the enforcement was reported to be costly and lacks efficacy and antagonize visitors [2]. Site managements including developing and hardening a recreation site or trail are less direct, but they were also costly and could permanently alter the natural landscape and the nature of recreation experiences [3, 4]. In contrast, a visitor education program is designed to persuade visitors to adopt low impact practices. It is considered as a more appropriate, light-handed, and indirect management response to reduce resource impacts and improve visitor experiences [5]. Many visitor education programs such as Leave No Trace [6–8], Code of Conduct [9], and Environmental Guidelines for Tourists [10] have already been applied in many countries. But it is a pity that there has been no such programs being mandated in china where the tourism market is extremely booming and the resource and social impacts are serious. The advantages of visitor education programs have been recognized by more and more Chinese scholars and site managers in recent years. Some pilot studies of the visitor educational intervention to protect precious natural and cultural resource have been reported [11–13].

Many scholars from the U.S., Australia, and Canada have discussed the education messages contents (see, e.g., [14–16]), deliveries (see, e.g., [17–21]), audience characteristics (see, e.g., [22–24]) and theoretical grounding (see, e.g., [25–27]), the efficacy of educational efforts (see, e.g., [10, 28]), and so on. However, most of the studies only focus on the empirical results obtained by visitor education interventional experiment data [6]. Few researches were able to quantify the efficacy of the visitor education intervention in tourism management. The purpose of this investigation is to study the dynamics of Ecological Tourism System and the effects of visitor education programs on the population dynamics by using a stochastic dynamics model.

An Ecological Tourism System consists of the natural resources including air, rocks, mountains, soils, and trees and the stakeholders in visitation destination. The stakeholders [11] in ecotourism include the following seven groups:(1)the tourism resource supply group, including tourism resource managers and park managers;(2)the visitation product sellers group, including tourism companies and tourism website managers;(3)the group of related products or services suppliers, including sellers who are providing transportations, hotels, catering, and other tourism memories products;(4)the group of local communities, including the staff of tourism companies and workers in other companies and local people;(5)the government group, including center government, provincial governments, and municipalities;(6)the special interest group, which are often shown as nongovernment organizations (NGOs) including environmental organizations, human rights and labor rights organizations, Trust and environmental charity institutes, academics, and social media;(7)visitors, which are the core consumers of tourism products and services.

The stakeholders in ecotourism destinations cannot survive without the support from one another, and they must unite to a population in a certain space [11]. It is valuable to study the dynamic behavior of the Ecological Tourism System under the visitor educational intervention. Wei et al. (2013) proposed a stochastic dynamics to study the efficacy of visitor education interventions [29]. Applying stochastic differential equation theory [30], they proved that the Ecological Tourism System can develop sustainably when perturbed by the continual visitor educational intervention. But they only considered one parameter of the system that was perturbed by visitor education intervention. This study follows the same thread of research in Wei et al. (2013) and further assumes that three parameters (i.e., the intrinsic increasing rate of stakeholders, the immigration rate of stakeholders from outside of the system, and the emigration rate of stakeholders from the inner system) are perturbed by the visitor education intervention stochastically. The Multiparameter Stochastic Dynamics of Ecological Tourism System (MPSDETS) is proposed to study the efficacy of the visitor educational programs.

In order to clearly present Multiparameters Stochastic Dynamics of Ecological Tourism System in the following three sections, the key notations and definitions of variables are introduced as follows: : stakeholders in Ecological Tourism System at time ; : initial value of Ecological Tourism System at time ; : intrinsic increasing rate of stakeholders within Ecological Tourism System; : immigration rate of stakeholders from outside of Ecological Tourism System; : emigration rate of stakeholders from the inner Ecological Tourism System; : maximum carrying capacity of Ecological Tourism System; : Brownian motion; : intensity of visitor education interventions; : explosion time of stochastic differential equation; : stopping time of stochastic differential equation; : a sufficient large number; : any nonnegative integer; : a -function ; : a sufficient large positive number; : a sufficient large positive constant.

We define the system stability state, collapsed state, and chaos state as follows. When the number of stakeholders is smaller than the maximum capacity of system, then the system is under stability sate. When explodes to infinity or , then the system collapsed. When is a little larger than or near the maximum capacity of system, then the system is under chaos state. When the system is under stability state, then there are proper stakeholders in the system and the ETS are running safely. When the system is under chaos state, then there are too much stakeholders in the system and the ETS are running dangerously. When the system is under collapsed state, then the system cannot accept any visitors and the ETS are damaged seriously.

#### 2. The Multiparameter Stochastic Dynamics of Ecological Tourism System

There are stakeholders in the Ecological Tourism System (ETS) at time . Some stakeholders may enter the system from outside of ETS while some others leave the system after they finished their visitation activities or business affairs. is the maximum carrying capacity of the ETS, and then is the recreation impacts factor per capita. Then, means the recreation impact factor at time . The Ordinary Differential Equation of Ecological Tourism System (ODETS) without visitor educational intervention can be expressed as The solution of Ordinary Differential Equation (1) cannot avoid explosion when given a large initial value .

Since the variable denotes the population of stakeholders, must be positive value on time . When the parameters of (1) , , (1) only has a local solution as follows:Given that , it is not difficult to see that explodes to infinity when Hence, the solution of (1) explodes to infinity in a finite time . The initial quantity of visitors is very large in some protected area especially on the public vocation in China as mentioned in Section 1. So, the solution of the differential dynamic equation explodes inevitably. In other words, the ETS is unstable and easy to collapse. The Ecological Tourism System cannot develop sustainably in that case.

Now, consider the case that ETS is intervened by visitor education. In this case, the intrinsic rate of increase is stochastically perturbed with where . As a result, the stochastic dynamic model of Ecological Tourism System (SDETS) with one stochastic parameter can effectively suppress the explosion of the solution of (1) in a finite time with any initial value

Figure 1 shows that ODETS (1) without visitor educational intervention explodes at time when , , , , and . However, SDETS (5) does not explode at any time with the same parameters shown in Figure 1. Similar results can be drawn from Figure 2 for the case with parameters , , , , and . ODETS (1) explodes at time as shown in Figure 2, whereas the stochastic system is more tempered. However, Figure 2 also shows that the number of population of the stochastic system reaches a large number (1083), greater than the maximum capacity carrying of the system at . Although the value of stakeholders returns shortly to a normal level in a short time, the system is unstable or in chaos at the small interval time . It is obvious to see that SDETS with single-parameter intervention is not effective to describe the system at some time. Hence, it is valuable to improve SDETS (5) to avoid the chaos situation.

It is easy to find that other parameters , , in (1) can also affect the stability of the ETS. And these parameters are also stochastically perturbed by visitor educational intervention. It is important to study how the stochastic parameters affect Ecological Tourism System.

Because many stakeholders from outside competitors can immigrate into the system randomly, the immigration rate of the system may be larger than the acceptable level without visitor education intervention. If the immigration rate of the system is too large, then the stakeholders increase too fast to over the limit capacity in a short time leading to the chaos situation of the system. In fact, the immigration rate can be suppressed by the visitor education intervention. The effective visitor education programs can induce visitors and other stakeholders entering the system orderly, but the effects are changing stochastically. So, the immigration rate is perturbed by visitor education intervention stochastically as follows:where is a Brownian motion and is the intensity of visitor education.

Stakeholders such as visitors and other relative product suppliers can leave ETS when they have finished their activities and business. The emigration rate of stakeholders may be lower than the acceptable level of the system without visitor educational intervention. Then, the system may be too crowded to keep developing sustainably. The immigration rate of stakeholders can be restored to normal level because the effective visitor education programs can lead visitors and other stakeholders to leave ETS in time. But the effects of visitor education are changing stochastically in different situation. So, the immigration rate of stakeholders also fluctuates according to the effects of visitor education interventionswhere is a Brownian motion and is the intensity of visitor education to emigration rate of stakeholders.

The visitors and other stakeholders within ETS may adopt lower impact behaviors, keeping visitation activities or business affairs more orderly, smoothly, and eco-friendly when the visitor education interventions are highly effective. Hence, the effective visitor educational interventions may enhance the capacity of ETS with the same tourism resources and facilities. The different intensity of visitor education intervention results in different maximum capacity of Ecological Tourism System. is the maximum capacity of ETS with specific tourism resources and visitor education interventions. Hence, is stochastically perturbed by the visitor education as follows:where and represent a Brownian motion and the intensity of visitor education, respectively.

Then, the Multiparameter Stochastic Dynamics of Ecological Tourism System (MPSDETS) perturbed by visitor educational interventions may be described by Ito equation as follows:where and .

#### 3. The Solutions of MPSDETS-Model

With any specified definition, is denoted as a complete probability space with a filtration . satisfies the right continuous and contains all p-null sets [29, 30]. is denoted as a multidimensional standard Brownian motion defined on the proceeding probability space.

Because is the total number of stakeholders, it should be nonnegative. Furthermore, if the stochastic differential equation has a unique global solution for any given initial value , the coefficients of the equation should satisfy the linear growth condition and the local Lipschitz condition [31–34]. But the coefficients only satisfy the locally Lipschitz continuous condition, so the solution of (9) may explode at a finite time [31, 34].

The purpose of this section is to discover whether visitor education interventions can suppress the explosion of dynamic system or not. It is natural to impose the hypothesis on the visitor education intensities as follows:

The following theorem reveals that visitor education interventions can suppress the explosion of the solution as suggested by Wei et al. (2013). In other words, the theorem shows that the solution of (9) is positive and global under the hypothesis (H1).

Theorem 1. *If , , for any system parameters and any given initial value , there is a unique value to (9) on and the solution is positive with probability one; namely, when almost surely.*

*Proof. *Because the coefficients of (9) satisfy the locally Lipschitz continuous, there is a local solution on for any given initial value , when is the explosion time [29]. In order to prove the solution of (9) is global and positive, it should be shown that almost surely. The initial value of the equation varies within the interval for any sufficient large positive value . The stopping time is denoted aswhere and is denoted as the empty set. It is obvious that is increasing when . One can see that ; hence, a.s.

In order to show that a.s. and for any , we should prove that firstly. In other words, the equivalent statement of the theorem shows that a.s. If the statement a.s. is not true, then there is a pair of constant and such that Hence, there is an integer such that The -function is defined as The function is nonnegative on .

Since and , it can be derived from the first or second order condition that for any . When , it can be derived from the Ito formula that whereRecalling the hypothesis (H1), it is not difficult to see that is bounded by sufficient large positive value . Therefore, it can be obtained that Integrating both sides of proceeding inequality from 0 to , we can yieldAnd then take expectations on the both sides of inequality (18), and we can yield for and , so that . Take any , by the definition of stopping time (11), and there is which is equal to either or . Then, we can derive that is not less than either So, there is or when . Thus, we can getTaking expectation on both sides of (20), we can yieldThen, we can get the following inequality from (22):Moreover, is a constant which is denoted in (12), so that . When , we can obtainThen, we can get the contradiction by (24). So, we can draw a conclusion that the statement of (12) is false. Furthermore, the original statement is true which is . In other words, the statement is true almost surely.

It is easy to derive from Theorem 1 that the solution of (9) exists and nerve explodes in a finite time. Following hypothesis (11), the visitor educational interventions on ETS can introduce stochastic parameters including stochastic intrinsic increasing rate of stakeholders and stochastic emigration rate of stakeholders, stochastic immigration rate of stakeholders, and stochastic maximum capacity of tourism ecological system. When , , the stochastic parameters can suppress the explosion of the solution. It means that the slight or tiny stochastic change of any parameter introduced by the visitor educational intervention can maintain the stability and sustainability of the ETS. Therefore, the visitor educational intervention is necessary condition of safe and stable developing of ETS. Without the visitor education intervention, the system would collapse in a finite time as mentioned in Section 1. For instance, tens of thousands of tourists were trapped in the Huashan Mountain because of the conflict between visitors and park mangers.

#### 4. The Bounded Property of MPSDETS-Model

Under the hypothesis (H1), Theorem 1 shows that the solution of MPSDETS (9) does not explode in a finite time. The next step is to show that the solution varies in a finite interval. is denoted as the unique global solution of (9) for any given initial value . The stochastically ultimate bounded property of the solution of the MPSDETS is defined as follows.

*Definition 2. *For any initial value , if the solution of (9) has the property that then (9) is stochastically and ultimately bounded with probability one. It is useful to present the following theorem which shows that (9) is stochastically ultimate bounded.

Theorem 3. *If the coefficients of (9) satisfy the property as follows: then the solution of (9) is stochastically ultimate bounded with probability one. Meanwhile, given any initial value , the solution of (9) has the property that *

*Proof. *For convenience, we can define -function byBy Ito formula, we can get is denoted asObviously, is quadratic function. Applying the properties of quadratic function, the following inequalities can be obtained:So, when (26) are satisfied, the quadratic function must be less than zero. Therefore, we can obtain that Integrating both sides of inequality (33) from 0 to , we can yield is defined asApplying Doob’s martingale convergence theorem, we can yieldThen, a.s. is obtained.

Recalling the definition of , if and only if or , the limit of is . Hence, a.s.

Then, it is not difficult to see thatGiven any initial value , it is easy to obtain that the solution of (9) varies in the interval of as shown in Theorem 3. So, the stakeholders can develop within an interval steadily. The population of stakeholders can be retained in a safe interval. ETS under visitor educational interventions can keep sustainable opportunities for high quality visitor experiences while avoiding or minimizing associated negative impacts to protected area resources although there are very large tourists in a short time. In other words, the visitation negative impacts in protected areas caused by large amount of visitors in public vocations can be minimized by the high efficient visitor education programs.

#### 5. Simulations

Following the Euler-Maruyama method [35], it is easy to get the numerical solutions of MPSDTES (9) with different parameters by MATLAB program. And the figures can be drawn according to the simulation data.

Figure 3 shows that SDETS (2) can effectively suppress the explosion of Ordinary Differential Equation (1) when , , , , and . But SDETS (2) cannot avoid some chaos situations either. For example, the peak value of visitors and other stakeholders is larger than the maximum capacity . Although decreases to , the system is chaotic in the short interval time . The chaos situations appear two times in an interval time when , , , , and . Tables 1(a) and 1(b) show that chaos situations are getting larger when given larger initial value . The times of chaos situations are smaller when given smaller comparing Tables 1(a) and 1(b).

Figure 3 shows that MPSDETS (9) can also effectively suppress the explosion of ODETS (1) when , , , , and . MPSDETS (9) can keep the peak value of visitors and other stakeholders under the maximum capacity of Ecological Tourism System compared with SDETS (2). The chaos situations disappear in the interval with the same parameters. Table 1(a) shows that MPSDETS (9) can reduce the times of chaos situation obviously. For example, SDETS (2) has 5 times of chaos situation when , but MPSDETS (9) has none. Figure 4 also shows the same results with , , , , and .

It is interesting to find that MPSDETS (9) is much more stable than SDETS (2) from the simulations. The simulations also show that ETS can develop more sustainably when all parameters are perturbed by continuous visitor educational intervention. So, the stakeholders in ETS should be involved in the visitor education programs. The stakeholders such as the visitation product sellers group and the special interest group play necessary roles in visitor education programs.

#### 6. Conclusions

Section 1 concludes that the fast expanding visitation in protected areas causes some serious natural resource and culture impacts in China and other countries. The popular management strategies to deal with the negative impacts include regulations, site managements, and visitor education programs. The regulation and site management gain very limit efficacy. It is easy to conclude from literatures that the fast expanding visitation in parks and the protected area including some natural history park is quickly damaged without visitor education. And the visitor education intervention programs are the most effective way to deal with negative visitation impacts as shown in the literatures. Researchers only discussed visitor educational message content, delivery, and audience characteristics, and so on. However, there are few mathematical models to study the efficacy of visitor education interventions in views of dynamics of ETS. So, it is valuable to study the dynamics of ETS with one stochastic parameter under visitor education interventions. A similar study has been done by Wei et al. in 2013.

In Section 2, ODETS is established to describe Ecological Tourism System without visitor educational intervention. The model reveals that the ecological tourism collapses without any perturbed factors such as visitor educational intervention. SDETS is effective to suppress the explosion of ODETS as shown in Figures 1 and 2. But it is not effective to suppress the peak value of visitors as shown in Figure 3. The peak value of visitors is larger than the maximum capacity of the tourism system in a short time. Hence, it is necessary to assume that the intrinsic increasing rate of stakeholders, the immigration rate of stakeholders from outside of the system, and the emigration rate of stakeholders from the inner system are stochastically perturbed by the visitor education interventions. Then, the Multiparameter Stochastic Dynamics of ETS perturbed by visitor educational interventions can be described by Ito equation.

In Sections 3 and 4, we prove that MPSDETS does not explode in a finite time. Then, we prove that the solution of MPSDETS varies in a finite interval. The solution shows that the population of the system retains in an acceptable level which can sustain the opportunities for high quality visitor experiences while avoiding or minimizing the associated negative impacts to the protected area resources and neighbors communities. The educational interventions increase visitors’ knowledge and alter visitors’ undesired behaviors to low impact behaviors. In other words, the tourism system can develop steadily and safely even under a large amount of visitors in public vocations when employing visitor education intervention programs.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by the National Science Fund of China (Grant no. 71272177/G020902, Grant no. 31100958), Social Science Foundation of Guangdong (Grant no. GD14XGL03), Youth Innovation Fund of Shenzhen Polytechnic (Grant no. 601422K35018), Project of Higher Vocational Education Research Society of China (Grant no. GZYLX1213345), and Philosophical and Social Science Project of Shenzhen (Grant no. 125C09). The authors would like to thank the editor and an anonymous referee for their helpful comments on the paper.