
Citation: Peter Warren. Business and policy models to incentivise utilities to engage with demand-side management[J]. Green Finance, 2019, 1(1): 4-29. doi: 10.3934/GF.2019.1.4
[1] | Jianxin Yang, Zhipeng Qiu, Xue-Zhi Li . Global stability of an age-structured cholera model. Mathematical Biosciences and Engineering, 2014, 11(3): 641-665. doi: 10.3934/mbe.2014.11.641 |
[2] | Mingtao Li, Xin Pei, Juan Zhang, Li Li . Asymptotic analysis of endemic equilibrium to a brucellosis model. Mathematical Biosciences and Engineering, 2019, 16(5): 5836-5850. doi: 10.3934/mbe.2019291 |
[3] | Linhua Zhou, Meng Fan, Qiang Hou, Zhen Jin, Xiangdong Sun . Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China. Mathematical Biosciences and Engineering, 2018, 15(2): 543-567. doi: 10.3934/mbe.2018025 |
[4] | Zongmin Yue, Yuanhua Mu, Kekui Yu . Dynamic analysis of sheep Brucellosis model with environmental infection pathways. Mathematical Biosciences and Engineering, 2023, 20(7): 11688-11712. doi: 10.3934/mbe.2023520 |
[5] | Mingtao Li, Guiquan Sun, Juan Zhang, Zhen Jin, Xiangdong Sun, Youming Wang, Baoxu Huang, Yaohui Zheng . Transmission dynamics and control for a brucellosis model in Hinggan League of Inner Mongolia, China. Mathematical Biosciences and Engineering, 2014, 11(5): 1115-1137. doi: 10.3934/mbe.2014.11.1115 |
[6] | Zhisheng Shuai, P. van den Driessche . Impact of heterogeneity on the dynamics of an SEIR epidemic model. Mathematical Biosciences and Engineering, 2012, 9(2): 393-411. doi: 10.3934/mbe.2012.9.393 |
[7] | Yaoyao Qin, Xin Pei, Mingtao Li, Yuzhen Chai . Transmission dynamics of brucellosis with patch model: Shanxi and Hebei Provinces as cases. Mathematical Biosciences and Engineering, 2022, 19(6): 6396-6414. doi: 10.3934/mbe.2022300 |
[8] | Qiuyi Su, Jianhong Wu . Impact of variability of reproductive ageing and rate on childhood infectious disease prevention and control: insights from stage-structured population models. Mathematical Biosciences and Engineering, 2020, 17(6): 7671-7691. doi: 10.3934/mbe.2020390 |
[9] | Pengyan Liu, Hong-Xu Li . Global behavior of a multi-group SEIR epidemic model with age structure and spatial diffusion. Mathematical Biosciences and Engineering, 2020, 17(6): 7248-7273. doi: 10.3934/mbe.2020372 |
[10] | Tsanou Berge, Samuel Bowong, Jean Lubuma, Martin Luther Mann Manyombe . Modeling Ebola Virus Disease transmissions with reservoir in a complex virus life ecology. Mathematical Biosciences and Engineering, 2018, 15(1): 21-56. doi: 10.3934/mbe.2018002 |
Brucellosis, a zoonotic disease, is a natural epidemic disease that is not only prevalent among livestock and humans but also widely spread in wild animals [1,2]. There are human and animal brucellosis in most countries of the world, and high-risk areas of epidemics are mainly distributed in developing countries, such as Syria, Jordan, Zambia, Mongolia [3,4,5,6]. It has brought huge losses to the livestock industry worldwide and serious health problems to livestock-related practitioners, there are hundreds of thousands of new cases reported annually [6,7,8]. Infected domestic and wild animals and their excreta are the main source of infection, it is contact and pathogen infection which are the main modes of brucellosis transmission [9]. Human brucellosis is rarely transmitted to susceptible animals, and there is no infection between people reported [10,11]. Therefore, it is the eradication of animal brucellosis that is the only way to solving human health problem, and understanding the mechanisms and risk factors of the spread of brucellosis is one of the first problems that must be solved.
Brucellosis is transmitted to susceptible individuals mainly by contact with infectious individuals or by sucking pathogens from the environment [8]. In animals, the transmission of brucellosis mainly occurs between sexually mature animals, young animals may also be infected, but generally do not have any clinical manifestations and serological tests are usually non-positive, and the infection between young animals is so little as to be almost invisible [12]. In other words, sexually mature animals are very susceptible to brucellosis, latent infections may be found in young animals, but they are generally resistant [9]. It is important that many infected animals have a longer incubation period and that these animals remain serologically non-positive during this period. That is to say, latent animals may not be infectious and it is almost impossible to be found through detection, which is an important risk factor in the elimination of brucellosis [9]. Therefore, the impact of these mechanisms on brucellosis is worthy to be studied using mathematical models.
Statistical methods such as descriptive statistics, correlation analysis and time series analysis have been widely used for quantitative assessment of risk management measures for brucellosis, among which is very worthy of concern about the study of bison, elk and livestock brucellosis (see [13,14,15,16,17]). Theoretical studies on the impacts of the transmission mechanisms of brucellosis and applied studies on the assessments of risk management measures have also been studied using the kinetic model (see [18,19,20,21,22]). Especially in recent years, many kinetic models with indirect transmission have been established to analyze brucellosis transmission (see [23,24,25,26,27]). Although there have been many studies on brucellosis, there are still many transmission mechanisms and risk factors that are not considered in existing models. For example, infected individuals may have no infectivity in the early stage, during which time is different; and animals that are not sexually mature are hardly infected with direct and indirect modes [9]. Therefore, in the present work, a multi-stage dynamic model with distributed time delay is proposed involving the above risk factors and general nonlinear incidences. The existence and uniqueness of the endemic equilibrium is analyzed, and the local and global asymptotic stability of equilibria is proved.
The rest of this manuscript is constructed as follows. The dynamic model with distributed time delay and some preliminary results are given in Section 2. In Section 3, the content of research is the global dynamics of the disease-free equilibrium. The global stability of the equilibrium of persistent infection are analyzed in Section 4. In Section 5, the stability results are further explained by numerical simulation. A summary and further discussion is proposed in Section 6.
According to the transmission mechanisms of animal brucellosis, we classify animal population into three compartments: The young susceptible compartment S1(t), the adult susceptible compartment S2(t) (sexually mature), the infected compartment I(t). Similar to the definition in literature [25], B(t) can be defined as the concentration or number of pathogens in the environment. Some explanations and assumptions about the kinetic model are listed as follows. (Ⅰ) Since brucellosis can cause abortion and significantly reduce the survival rate of young animals, we suppose that infected animals have no birth rate. (Ⅱ) The supplementary rate of young susceptible animal population S1(t) is mainly derived from birth and import, it is assumed that the supplementary rate is A+bS2(t). (Ⅲ) Usually, there are two different stages for infected animals, which are no infection force in the early stage and are infectious carriers when infected animals begin to shed brucella. Therefore, we make use of a delay τ to describe the time from a susceptible individual to a infectious individual, and since the length of this time is varied, thus the delay τ is a distributed parameter in the interval [0,h], where h is the maximum value of the delay. As a result, the susceptible class S2(t) is reduced due to infection at the rate S2g(B)+S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ. (Ⅳ) Since the infected animal sheds the pathogen after the time τ, then the pathogen shedding rate of the infected animal is ∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ. (Ⅴ) In animal breeding, only the basic ewes breed newborns, and adult animals are widely used for slaughter or trade, then the birth rate b of adult animals is assumed to be less than the elimination rate μ2 of adult animals. Therefore, the modelling for animal brucellosis is given through the following distributed time-delay system:
{dS1dt=A+bS2−μ1S1−σ1S1,dS2dt=σ1S1−S2g(B)−S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ−μ2S2,dIdt=S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+S2g(B)−(μ2+c)I,dBdt=∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ−dB. | (2.1) |
Here, b is the birth rate of adult animals. μ1 and μ2 are the elimination rates of young and adult animals, respectively. σ1 is the transfer rate from young susceptible individuals to sexually mature individuals. The culling rate is defined by c and d is the decay rate of the pathogen. f(I) and g(B) are contact and indirect infection rates, respectively. ρ(τ) is a distributed function, it is non-negative, continuous and ∫h0ρ(τ)dτ=1.
The initial conditions for system (2.1) are given as follows:
{S1(x)=ϕ1(x),S2(x)=ϕ2(x),I(x)=ϕ3(x),B(x)=ϕ4(x),x∈[−h,0],h>0,ϕ=(ϕ1,ϕ2,ϕ3,ϕ4)∈C+⊂C. | (2.2) |
Here, C denotes the Banach space C([−h,0],ℜ4) of continuous functions mapping the interval [−h,0] into ℜ4 with the sup-norm ‖ϕ‖=supx∈[−h,0]|ϕi(x)|,i=1,2,3,4 for ϕ∈C. The nonnegative cone of C is defined as C+=C([−h,0],ℜ4+).
In order to make epidemiological significance for system (2.1), f,g and h are assumed to be second-order continuous differentiable functions and satisfy the following hypotheses:
(H1) f(0)=g(0)=h(0)=0 and f(I),g(B),h(I)>0 for I,B>0;
(H2) f′(I),g′(B)>0 and f″(I),g″(B)≤0 for B,I≥0;
(H3) h′(I)>0 and h″(I)≤0 for I≥0.
The function f may be βIp or saturation incidences kln(1+λIk) and βIp1+kIp with constants β,p,λ,k>0 [28,29]. The pathogen infection rate g may be λB1+TB or η(1−e−αB) with constants λ,η,α>0 and T≥0 [30]. The function h may be kI with constants k>0 [31].
It is easy to verify that the system (2.1) always has a disease-free equilibrium E0=(S01,S02,0,0), where
S01=Aμ2μ2(μ1+σ1)−bσ1,S02=Aσ1μ2(μ1+σ1)−bσ1. |
The reproduction number of sytem (2.1) is given by the following expression:
R0=nS02fI(0)μ2+c+nS02gB(0)hI(0)d(μ2+c)=R01+R02. | (2.3) |
where n=∫h0ρ(τ)e−(μ2+c)τdτ. e−(μ2+c)τ denotes the probability of survival from a newly infected susceptible individual to an infectious individual. n is the total survival rate of infected individuals. fI(0) is the infection rate by an infectious individual. gB(0) is the infection rate of brucella. 1μ2+c is the average life span of infected individual. hI(0)d represents the total number of brucella shed by infectious individuals. Therefore, based on the explanation in the literature [26], R0 relies on direct and indirect infection and can be divided into two parts R01 and R02.
In the next section, we establish that the solution of system (2.1) is non-negative and bounded and analyze the uniqueness of the equilibrium point of persistent infection. The following qualitative results indicate the nonnegativity and boundedness of the solution of system (2.1) with initial conditions (2.2).
Theorem 2.1. (S1(t),S2(t),I(t),B(t)) is the solution of system (2.1) with the initial conditions (2.2), then S1(t),S2(t),I(t),B(t) are nonnegative and ultimately bounded.
Proof. Since (˙S1(t)+˙S2(t))=A>0 for t∈[0,+∞) when S1(t)=S2(t)=0, it implies that S1(t)+S2(t)≥0 for all t∈[0,+∞). Therefore, if S1(t)=0, then S2(t)≥0 and ˙S1(t)=A+bS2>0, so that S1(t)≥0. We denote
a1(t)=∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+g(B)+μ2,a2(t)=S2(t)(∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+g(B)), |
and
a3(t)=∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ. |
From the last three equations in (2.1), it can conclude that
S2(t)=S2(0)e−∫t0a1(η)dη+σ1∫t0e−∫tηa1(ξ)dξS1(η)dη≥0,I(t)=I(0)e−(μ2+c)t+∫t0e−(μ2+c)(t−η)a2(η)dη≥0, |
and
B(t)=B(0)e−dt+∫t0e−d(t−η)a3(η)dη≥0 |
for all t≥0. Thus, S1(t),S2(t),I(t),B(t)≥0 for t≥0.
We now analyze the boundedness of the solution of system (2.1). In fact,
(˙S1(t)+˙S2(t)+˙I(t))≤A−μ1S1−(μ2−b)S2−(μ2+c)I |
≤A−μ(S1(t)+S2(t)+I(t)), |
where μ=min{μ1,μ2−b}.
Hence, lim supt→∞(S1(t)+S2(t)+I(t))≤Aμ, it implies that lim supt→∞I(t)≤Aμ and lim supt→∞B(t)≤ndh(Aμ). Therefore, S1(t),S2(t),I(t),B(t) are ultimately bounded.
Since the solution of system (2.1) is nonnegative and ultimately bounded, the set
Ω={(S1(⋅),S2(⋅),I(⋅),B(⋅))∈C+:‖S1+S2+I‖≤Aμ,‖B‖≤ndh(Aμ)}, |
is positively invariant for system (2.1).
In order to analyze the uniqueness of the positive solution of system (2.1), the following lemma is given:
Lemma 2.1. Assume that conditions (H1)−(H3) are satisfied, the functions f(I)I, h(I)I, g(B)B and g(h(I))I are monotonic decreasing for I,B>0.
Proof. Since f″(I)≤0, it shows that f′(I) is monotonic decreasing, it follows that
f(I)I=f(I)−f(0)I−0=f′(ξ1)≥f′(I),ξ1∈(0,I), |
and
(f(I)I)′=f′(I)I−f(I)I2≤0. |
That is to say, f(I)I is a monotonically decreasing function. On the basis of the above method, it can show that h(I)I and g(B)B are also monotonic decreasing. Noting that
g(h(I))I=g(h(I))−g(h(0))h(I)−h(0)h(I)I=g′(h(ξ2))h′(ξ3)≥g′(h(I))h′(I),ξ2,ξ3∈(0,I), |
It can deduce that
(g(h(I))I)′=g′(h(I))h′(I)I−g(h(I))I2≤0. |
Therefore, g(h(I))I is also monotonic decreasing.
For system (2.1), the endemic equilibrium E∗=(S∗1,S∗2,I∗,B∗) can be derived from the following algebraic equations:
{A+bS2=μ1S1+σ1S1,σ1S1=S2(nf(I)+g(B))+μ2S2,S2(nf(I)+g(B))=(μ2+c)I,nh(I)=dB. | (2.4) |
By direct calculation, it can be written as
{S1=A+bS2μ1+σ1,m(S2)=S2(nf(I)+g(B)),S2(nf(I)+g(B))=(μ2+c)I,B=H(I). |
where m(S2)=Aσ1−(μ2(μ1+σ1)−bσ1)S2μ1+σ1 and H(I)=ndh(I).
Let us define
F1(S2,I)≜ |
Similar to the method of analysis in literature [26], using Lemma 2.1, the following result can be summarized:
Theorem 2.2. Assume that conditions (H_1)-(H_3) hold. Then there is a unique positive solution E_* = (S_1^*, S_2^*, I^*, B^*) of system (2.1) if and only if R_0 > 1 .
In this section, we show that the global stability of the equilibrium point E_0 of system (2.1) is independent of the initial value. The following conclusions are first obtained.
Lemma 3.1. Assume that conditions (H_1)-(H_{3}) hold. The disease-free equilibrium E_0 is locally asymptotically stable if R_0\leq1 and is unstable if R_0 > 1 .
Proof. The characteristic equation at E_0 is
\begin{equation} \left|\begin{array}{cccc} \lambda+\mu_1+\sigma_1 & -b & 0 & 0\\[0.1cm] -\sigma_1 & \lambda+\mu_2 & S_2^0f_I(0)\Gamma(\lambda) & S_2^0g_B(0)\\[0.1cm] 0 & 0 & \lambda+\mu_2+c-\Gamma(\lambda)S_2^0f_I(0) & -S_2^0g_B(0)\\[0.1cm] 0 & 0 & -h_I(0)\Gamma(\lambda) & \lambda+d \end{array}\right| = 0, \end{equation} | (3.1) |
where
\Gamma(\lambda) = \int_0^{h}\rho(\tau)e^{-(\mu_2+c+\lambda)\tau}d\tau. |
It follows from (3.1) that
\begin{equation*} \begin{split} &(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2) ((\lambda+d)(\lambda+\mu_2+c-\Gamma(\lambda)S_2^0f_I(0))\\ &-\Gamma(\lambda)S_2^0g_B(0)h_I(0))\\ = &b\sigma_1((\lambda+d)(\lambda+\mu_2+c-\Gamma(\lambda)S_2^0f_I(0))\\ &-\Gamma(\lambda)S_2^0g_B(0)h_I(0)).\\ \end{split} \end{equation*} |
That is,
(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)(\lambda+d)(\lambda+\mu_2+c) -b\sigma_1(\lambda+d)(\lambda+\mu_2+c) |
\quad\hskip 0.1cm = ((\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1) ((\lambda+d)\Gamma(\lambda)S_2^0f_I(0)+\Gamma(\lambda)S_2^0g_B(0)h_I(0)). |
It follows that
H_1(\lambda) ((\lambda+d)(\lambda+\mu_2+c)-(\lambda+d)\Gamma(\lambda)S_2^0f_I(0)-\Gamma(\lambda)S_2^0g_B(0)h_I(0)) = 0, |
where H_1(\lambda) = (\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1 . Since \mu_1+\sigma_1+\mu_2 > 0 and \mu_2(\mu_1+\sigma_1)-b\sigma_1 > 0 , then H_1(\lambda) consists of two roots which are negative real parts. Therefore, we only analyze the distribution of the roots of the following equation:
\begin{equation} H_2(\lambda) = (\lambda+d)(\lambda+\mu_2+c)-(\lambda+d)\Gamma(\lambda)S_2^0f_I(0) -\Gamma(\lambda)S_2^0g_B(0)h_I(0) = 0. \end{equation} | (3.2) |
Assume now that R_0 > 1 , then
H_2(0) = d(\mu_2+c)-dnS_2^0f_I(0)-nS_2^0g_B(0)h_I(0)) \lt 0, \hskip 0.5cm H_2(+\infty) = +\infty. |
Hence H_2(\lambda) has at least one positive root in [0, +\infty) , then E_0 is unstable if R_0 > 1 .
From (3.2), we have
\begin{equation*} \begin{split} (\lambda+d)(\lambda+\mu_2+c) = &\Gamma(\lambda)(\lambda S_2^0f_I(0)+dS_2^0f_I(0)+S_2^0g_B(0)h_I(0))\\ = &\frac{(\mu_2+c)F(\lambda)}{n}(\lambda R_{01}+dR_0),\\ \end{split} \end{equation*} |
or
\begin{equation} (\lambda+d)(\frac{\lambda}{\mu_2+c}+1) = \frac{R_0F(\lambda)}{n}(\lambda \frac{R_{01}}{R_0}+d). \end{equation} | (3.3) |
Next, considering the case R_0\leq1 . If \lambda = x+yi is a solution of (3.3), one shows that x < 0 . Otherwise, x\geq 0 implies
\left| \lambda+d\right| \gt \left| \lambda\frac{R_{01}}{R_0}+d\right| , \hskip 0.2cm \left| \frac{\lambda}{\mu_2+c}+1\right| \gt 1, \hskip 0.2cm \left| \frac{R_0F(\lambda)}{n}\right| \leq 1, |
and thus
\left|(\lambda+d)(\frac{\lambda}{\mu_2+c}+1)\right| \gt \left|\frac{R_0F(\lambda)}{n}(\lambda \frac{R_{01}}{R_0}+d)\right|, |
this is a contradiction to (3.3). Therefore, all roots of Eq (3.3) have no zero and positive real parts when R_0\leq1 , this shows that E_0 is locally asymptotically stable.
Theorem 3.1. Assume that conditions (H_1)-(H_{3}) are established. If R_0\leq1 , the disease-free equilibrium E_0 = (S_1^0, S_2^0, 0, 0) of system (2.1) is globally asymptotically stable.
Proof. Since the functions \frac{f(I)}{I} , \frac{g(B)}{B} and \frac{h(I)}{I} are decreasing, then we have
\begin{split} \frac{nS_2f(I)}{(\mu_2+c)I}\leq &\lim\limits_{I\rightarrow0^+}\frac{nS_2^0f(I)}{(\mu_2+c)I} = \frac{nS_2^0f_I(0)}{\mu_2+c}\triangleq b_1,\\ \frac{S_2g(B)}{dB}\leq&\lim\limits_{B\rightarrow0^+}\frac{S_2^0g(B)}{dB} = \frac{S_2^0g_B(0)}{d}\triangleq b_2,\\ \frac{nh(I)}{(\mu_2+c)I}\leq&\lim\limits_{I\rightarrow0^+}\frac{nh(I)}{(\mu_2+c)I} = \frac{nh_I(0)}{\mu_2+c}. \end{split} |
Define
J = \left(\begin{array}{cc} 1 & 0 \\ -\frac{nh_I(0)}{\mu_2+c} & 1\\ \end{array} \right), \hskip 0.5cm (a_1, a_2) = (b_1, b_2)J^{-1}. |
We find a_1 = R_0 and define a Lyapunov functional L as follows:
L(t) = L_1(t)+L_2(t)+L_3(t), |
where
\begin{split} L_1 = &R_0(S_1-S_1^0-S_1^0\ln\frac{S_1}{S_1^0}+S_2-S_2^0-S_2^0\ln\frac{S_2}{S_2^0}+I)+a_2B,\\ L_2 = &R_0S_2^0\int_0^{h}\Psi(\tau)f(I(t-\tau))d\tau,\\ L_3 = &a_2\int_0^{h}\Psi(\tau)h(I(t-\tau))d\tau, \end{split} |
and
\Psi(\tau) = \int_\tau^{h}\varphi(s)ds, \hskip 0.5cm \varphi(s) = \rho(s)e^{-(\mu_2+c)s}. |
Then the derivative of L_1 along the positive solutions of system (2.1) is
\begin{split} \frac{dL_1}{dt} = &R_0(1-\frac{S_1^0}{S_1})\frac{dS_1}{dt} +R_0(1-\frac{S_2^0}{S_2})\frac{dS_2}{dt} +R_0\frac{dI}{dt}+a_2\frac{dB}{dt}\\ = &R_0(2A+bS_2^0+\sigma_1S_1^0-\mu_1S_1-(\mu_2-b)S_2-S_1^0\frac{A}{S_1} -bS_1^0\frac{S_2}{S_1}-\sigma_1S_2^0\frac{S_1}{S_2})\\ &+R_0(S_2^0g(B)+S_2^0\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau-(\mu_2+c)I)\\ &+a_2(\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau-dB).\\ \end{split} | (3.4) |
Calculating the derivative of L_2(t) along the solutions of system (2.1), one obtains
\begin{split} \frac{dL_2}{dt} = &R_0S_2^0\int_0^{h}\Psi(\tau)\frac{df(I(t-\tau))}{dt}d\tau\\ = &-R_0S_2^0\int_0^{h}\Psi(\tau)\frac{df(I(t-\tau))}{d\tau}d\tau\\ = &-R_0S_2^0(\Psi(\tau)f(I(t-\tau)))|_0^h+R_0S_2^0\int_0^{h}\frac{d\Psi(\tau)}{d\tau}f(I(t-\tau))d\tau\\ = &R_0S_2^0nf(I(t))-R_0S_2^0\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau.\\ \end{split} | (3.5) |
Similar to the above-used method, it can obtain that
\frac{dL_3}{dt} = a_2nh(I(t))-a_2\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau. | (3.6) |
Combining the Eqs (3.4), (3.5) and (3.6), it follows that
\begin{split} \frac{dL}{dt} = &\frac{dL_1}{dt}+\frac{dL_2}{dt}+\frac{dL_3}{dt}\\ = &R_0\mu_1S_1^0(2-\frac{S_1}{S_1^0}-\frac{S_1^0}{S_1}) +R_0bS_2^0(2-\frac{S_1S_2^0}{S_1^0S_2}-\frac{S_1^0S_2}{S_1S_2^0})\\ &+R_0(\mu_2-b)S_2^0(3-\frac{S_2}{S_2^0}-\frac{S_1S_2^0}{S_1^0S_2}-\frac{S_1^0}{S_1})\\ &+R_0(\frac{nS_2^0f(I)}{(\mu_2+c)I}, \frac{S_2^0g(B)}{dB})((\mu_2+c)I, dB)^T\\ &-(R_0, a_2) \left(\begin{array}{cc} 1 & 0\\ -\frac{nh(I)}{(\mu_2+c)I} & 1\\ \end{array} \right) ((\mu_2+c)I, dB)^T\\ \leq &R_0(\frac{nS_2^0f_I(0)}{(\mu_2+c)}, \frac{S_2^0g_B(0)}{d}) ((\mu_2+c)I, dB)^T\\ &-(R_0, a_2) \left(\begin{array}{ccc} 1 & 0\\ -\frac{nh_I(0)}{\mu_2+c} & 1\\ \end{array} \right) ((\mu_2+c)I, dB)^T\\ = &(R_0-1)(b_1, b_2)((\mu_2+c)I, dB)^T\\ \leq&0.\\ \end{split} |
Similar to the analysis of Theorem 1 in the literature [32], it follows from Lemma 3.1 that the disease-free equilibrium E_0 is globally asymptotically stable by LaSalle's Invariance Principle [33].
According to Lemma 3.1, the disease-free steady state solution E_0 is unstable when R_0 > 1 . Using Theorem 4.2 in [34], the uniform permanence of system (2.1) can be proven, the process is ignored here. In this following section, by constructing a Lyapunov functional, the global stability of the equilibrium of persistent infection is proved. We first analyze its local stability.
Lemma 4.1. Assume that conditions (H_1)-(H_{3}) are ture, If R_0 > 1 , the endemic equilibrium E_* of system (2.1) exists and is locally asymptotically stable.
Proof. For system (2.1), the characteristic equation at E_* is
\left|\begin{array}{cccc} \lambda+\mu_1+\sigma_1 & -b & 0 & 0\\[0.1cm] -\sigma_1 & H_3(\lambda) & \Gamma(\lambda)S_2^*f_I(I^*) & S_2^*g_B(B^*)\\[0.1cm] 0 & -nf(I^*)-g(B^*) & H_4(\lambda) & -S_2^*g_B(B^*)\\[0.1cm] 0 & 0 & -\Gamma(\lambda)h_I(I^*) & \lambda+d \end{array}\right| = 0, |
where H_3(\lambda) = \lambda+nf(I^*)+g(B^*)+\mu_2 , H_4(\lambda) = \lambda+\mu_2+c-\Gamma(\lambda)S_2^*f_I(I^*) . By simple calculation, one obtains
\begin{split} &(\lambda+d)(\lambda+\mu_2+c)\Phi_1(\lambda)\\ = &\Phi_2(\lambda)\Gamma(\lambda)((\lambda+d)S_2^*f_I(I^*)+h_I(I^*)S_2^*g_B(B^*))\\ = &\Phi_2(\lambda)\frac{\Gamma(\lambda)}{n}\Lambda n (\lambda \frac{S_2^*f_I(I^*)}{\Lambda}+d), \end{split} |
or
\begin{equation} (\lambda+d)(\frac{\lambda}{\mu_2+c}+1)\Phi_1(\lambda) = \Phi_2(\lambda)\frac{\Gamma(\lambda)}{n}\frac{\Lambda n}{\mu_2+c} (\lambda \frac{S_2^*f_I(I^*)}{\Lambda}+d), \end{equation} | (4.1) |
where
\begin{split} \Phi_1(\lambda) = &(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2+nf(I^*)+g(B^*))-b\sigma_1,\\ \Phi_2(\lambda) = &(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1,\\ \Lambda = &S_2^*f_I(I^*)+\frac{1}{d}S_2^*h_I(I^*)g_B(B^*)\\ \end{split} |
Noting that from Lemma 1
\begin{split} \frac{\Lambda n}{\mu_2+c} = &\frac{S_2^*f_I(I^*)n+\frac{1}{d}S_2^*h_I(I^*)g_B(B^*) n}{\mu_2+c}\\ \leq &\frac{S_2^*(\frac{f(I^*)n}{I^*}+\frac{n}{d}\frac{h(I^*)}{I^*}\frac{g(B^*)}{B^*})}{\mu_2+c}\\ = &\frac{S_2^*(f(I^*)n+g(B^*))}{I^*(\mu_2+c)}\\ = &1. \end{split} |
Assume \lambda = x+yi is a solution of (4.1). If x\geq0 , we then have
\left| \lambda+d\right| \gt \left| \lambda\frac{S_2^*f_I(I^*)}{\Lambda}+d\right| , \hskip 0.1cm \left| \frac{\lambda}{\mu_2+c}+1\right| \gt 1, \hskip 0.1cm \left| \frac{F(\lambda)}{n}\right| \leq 1, |
and
\begin{split} \left| \frac{\Phi_1(\lambda)}{\Phi_2(\lambda)}\right| = &\left| \frac{(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2+nf(I^*)+g(B^*))-b\sigma_1} {(\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1}\right|\\ = &\left| \frac{\lambda+\mu_2+nf(I^*)+g(B^*)-\frac{b\sigma_1}{\lambda+\mu_1+\sigma_1}} {\lambda+\mu_2-\frac{b\sigma_1}{\lambda+\mu_1+\sigma_1}}\right|\\ = &\left| \frac{M_1+nf(I^*)+g(B^*)+M_2i}{M_1+M_2i}\right|\\ \gt &1, \end{split} |
where M_1 = x+\mu_2-\frac{b\sigma_1(x+\mu_1+\sigma_1)}{(x+\mu_1+\sigma_1)^2+y^2} > 0 , M_2 = y+\frac{b\sigma_1y}{(x+\mu_1+\sigma_1)^2+y^2} . So it concludes that
\left|(\lambda+d)(\frac{\lambda}{\mu_2+c}+1)\Phi_1(\lambda)\right| \gt \left| \Phi_2(\lambda)\frac{F(\lambda)}{n}\frac{\Lambda n}{\mu_2+c} (\lambda \frac{S_2^*f_I(I^*)}{\Lambda}+d)\right|, |
which is a contradiction to (4.1). Therefore, the Eq (4.1) can not have any roots with a nonnegative real part, this implies that E_* is locally asymptotically stable if R_0 > 1 .
Theorem 4.1. Assume that conditions (H_1)-(H_{3}) hold. If R_0 > 1 , the endemic equilibrium E_* = (S_1^*, S_2^*, I^*, B^*) of system (2.1) is globally asymptotically stable.
Proof. Define
\begin{split} L_1 = &S_1-S_1^*-S_1^*\ln\frac{S_1}{S_1^*}+S_2-S_2^*-S_2^*\ln\frac{S_2}{S_2^*}\\ &+I-I^*-I^*\ln\frac{I}{I^*}+\frac{S_2^*g(B^*)}{nh(I^*)}(B-B^*-B^*\ln\frac{B}{B^*}).\end{split} |
Finding the time derivative of L_1 along the positive solutions of system (2.1) gives
\begin{split} \frac{dL_1}{dt} = &(1-\frac{S_1^*}{S_1})\frac{dS_1}{dt}+ (1-\frac{S_2^*}{S_2})\frac{dS_2}{dt}+(1-\frac{I^*}{I})\frac{dI}{dt}+ \frac{S_2^*g(B^*)}{nh(I^*)}(1-\frac{B^*}{B})\frac{dB}{dt}\\ = &2A+bS_2^*+\sigma_1S_1^*-\mu_1S_1-(\mu_2-b)S_2-S_1^*\frac{A}{S_1}-bS_1^*\frac{S_2}{S_1} -\sigma_1S_2^*\frac{S_1}{S_2}\\ &+S_2^*g(B)+S_2^*\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau-(\mu_2+c)I\\ &-\frac{I^*}{I}S_2g(B)-\frac{I^*}{I}S_2\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau\\ &+\frac{S_2^*g(B^*)}{nh(I^*)}(\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau-dB+dB^*)\\ &-\frac{S_2^*g(B^*)}{nh(I^*)}\frac{B^*}{B}\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau.\\ \end{split} |
By the Eq (2.4), we get
\begin{split} \frac{dL_1}{dt} = &\mu_1S_1^*(2-\frac{S_1}{S_1^*}-\frac{S_1^*}{S_1}) +(\mu_2-b)S_2^*(3-\frac{S_2}{S_2^*}-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*}{S_1})\\ &+bS_2^*(2-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*S_2}{S_1S_2^*}) +nS_2^*f(I^*)(3-\frac{S_1^*}{S_1}-\frac{S_1S_2^*}{S_2S_1^*}-\frac{I}{I^*})\\ &+S_2^*g(B^*)(3+\frac{g(B)}{g(B^*)}-\frac{S_1^*}{S_1}-\frac{S_1S_2^*}{S_2S_1^*}-\frac{I}{I^*} -\frac{S_2g(B)I^*}{S_2^*g(B^*)I})\\ &+S_2^*\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau -\frac{I^*}{I}S_2\int_0^{h}\varphi(\tau)f(I(t-\tau))d\tau\\ &+S_2^*g(B^*)(\frac{h(I)}{h(I^*)}-\frac{B}{B^*}-\frac{B^*h(I)}{Bh(I^*)}+1)\\ &+\frac{S_2^*g(B^*)}{nh(I^*)}(\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau-nh(I)+\frac{B^*}{B}nh(I))\\ &-\frac{S_2^*g(B^*)}{nh(I^*)}\frac{B^*}{B}\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau.\\ \end{split} | (4.2) |
Since the function \nu(x) = 1-x+\ln x is nonpositive for x > 0 and \nu(x) = 0 if and only if x = 1 , so we define
L_2 = -S_2^*f(I^*)\int_0^{h}\Psi(\tau)\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau. |
A direct calculation shows that
\begin{split} \frac{dL_2}{dt} = &-S_2^*f(I^*)\int_0^{h}\Psi(\tau)\frac{d}{dt}\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau\\ = &S_2^*f(I^*)\int_0^{h}\Psi(\tau)\frac{d}{d\tau}\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau\\ = &S_2^*f(I^*)(\Psi(\tau)\nu(\frac{f(I(t-\tau))}{f(I^*)})\mid_0^h +\int_0^{h}\varphi(\tau)\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau)\\ = &-S_2^*f(I^*)\int_0^{h}\varphi(\tau)\nu(\frac{f(I(t))}{f(I^*)})d\tau +S_2^*f(I^*)\int_0^{h}\varphi(\tau)\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau\\ = &S_2^*f(I^*)\int_0^{h}\varphi(\tau)(\frac{f(I)}{f(I^*)} -\frac{f(I(t-\tau))}{f(I^*)}+\ln\frac{f(I(t-\tau))}{f(I)})d\tau\\ = &nS_2^*f(I^*)\frac{f(I)}{f(I^*)}+S_2^*f(I^*)\int_0^{h}\varphi(\tau) (-\frac{f(I(t-\tau))}{f(I^*)}+\ln\frac{f(I(t-\tau))}{f(I)})d\tau.\\ \end{split} | (4.3) |
Define
L_3 = -\frac{S_2^*g(B^*)}{n}\int_0^{h}\Psi(\tau)\nu(\frac{h(I(t-\tau))}{h(I^*)})d\tau. |
Calculating the time derivative of L_3(t) , one obtains
\begin{split} \frac{dL_3}{dt} = &-\frac{S_2^*g(B^*)}{n}\int_0^{h}\Psi(\tau)\frac{d}{dt}\nu(\frac{h(I(t-\tau))}{h(I^*)})d\tau\\ = &\frac{S_2^*g(B^*)}{n}\int_0^{h}\Psi(\tau)\frac{d}{d\tau}\nu(\frac{h(I(t-\tau))}{h(I^*)})d\tau\\ = &\frac{S_2^*g(B^*)}{n}\int_0^{h}\varphi(\tau)(\frac{h(I)}{h(I^*)} -\frac{h(I(t-\tau))}{h(I^*)}+\ln\frac{h(I(t-\tau))}{h(I)})d\tau\\ = &S_2^*g(B^*)\frac{h(I)}{h(I^*)}+\frac{S_2^*g(B^*)}{n}\int_0^{h}\varphi(\tau) (-\frac{h(I(t-\tau))}{h(I^*)}+\ln\frac{h(I(t-\tau))}{h(I)})d\tau. \end{split} | (4.4) |
For system (2.1), the following Lyapunov functional is considered:
L = L_1+L_2+L_3. |
From (4.2), (4.3) and (4.4), we can get
\begin{split} \frac{dL}{dt} = &\mu_1S_1^*(2-\frac{S_1}{S_1^*}-\frac{S_1^*}{S_1}) +(\mu_2-b)S_2^*(3-\frac{S_2}{S_2^*}-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*}{S_1})\\ &+bS_2^*(2-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*S_2}{S_1S_2^*})\\ &+nS_2^*f(I^*)(3+\frac{f(I)}{f(I*)}-\frac{S_1^*}{S_1}-\frac{S_1S_2^*}{S_2S_1^*} -\frac{I}{I^*}-\frac{S_2f(I)I^*}{S_2^*f(I^*)I})\\ &+S_2^*g(B^*)(3+\frac{g(B)}{g(B^*)}-\frac{S_1^*}{S_1}-\frac{S_1S_2^*}{S_2S_1^*}-\frac{I}{I^*} -\frac{S_2g(B)I^*}{S_2^*g(B^*)I})\\ &+S_2^*g(B^*)(\frac{h(I)}{h(I^*)}-\frac{B}{B^*}-\frac{B^*h(I)}{Bh(I^*)}+1)\\ &+S_2^*f(I^*)\int_0^{h}\varphi(\tau)\mathcal{F}(\tau)d\tau +\frac{S_2^*g(B^*)}{n}\int_0^{h}\varphi(\tau)\mathcal{H}(\tau)d\tau, \end{split} | (4.5) |
where
\begin{split} \int_0^{h}\varphi(\tau)\mathcal{F}(\tau)d\tau = &\int_0^{h}\varphi(\tau)(\frac{S_2f(I)I^*}{S_2^*f(I^*)I} -\frac{S_2I^*f(I(t-\tau))}{S_2^*f(I^*)I}+\ln\frac{f(I(t-\tau))}{f(I)})d\tau\\ = &\int_0^{h}\varphi(\tau)(\nu(\frac{S_2I^*f(I(t-\tau))}{S_2^*f(I^*)I}) -\nu(\frac{S_2f(I)I^*}{S_2^*f(I^*)I}))d\tau\\ = &\int_0^{h}\varphi(\tau)\nu(\frac{S_2I^*f(I(t-\tau))}{S_2^*f(I^*)I})d\tau -n\nu(\frac{S_2f(I)I^*}{S_2^*f(I^*)I}), \end{split} | (4.6) |
and
\begin{split} \int_0^{h}\varphi(\tau)\mathcal{H}(\tau)d\tau = &\int_0^{h}\varphi(\tau)(\frac{h(I)B^*}{h(I^*)B} -\frac{B^*h(I(t-\tau))}{h(I^*)B}+\ln\frac{h(I(t-\tau))}{h(I)})d\tau\\ = &\int_0^{h}\varphi(\tau)\nu(\frac{B^*h(I(t-\tau))}{h(I^*)B})d\tau -n\nu(\frac{h(I)B^*}{h(I^*)B}). \end{split} | (4.7) |
it follows from (4.5), (4.6) and (4.7) that
\begin{split} \frac{dL}{dt} = &\mu_1S_1^*(2-\frac{S_1}{S_1^*}-\frac{S_1^*}{S_1}) +(\mu_2-b)S_2^*(3-\frac{S_2}{S_2^*}-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*}{S_1})\\ &+bS_2^*(2-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*S_2}{S_1S_2^*}) +nS_2^*f(I^*)(\frac{f(I)}{f(I^*)}-1)(1-\frac{f(I^*)I}{f(I)I^*})\\ &+nS_2^*f(I^*)(\nu(\frac{S_1^*}{S_1})+\nu(\frac{S_1S_2^*}{S_1^*S_2}) +\nu(\frac{S_2f(I)I^*}{S_2^*f(I^*)I})+\nu(\frac{f(I^*)I}{f(I)I^*}))\\ &+S_2^*g(B^*)((\frac{g(B)}{g(B^*)}-1)(1-\frac{g(B^*)B}{g(B)B^*}) +(\frac{h(I)}{h(I^*)}-1)(1-\frac{h(I^*)I}{h(I)I^*}))\\ &+S_2^*g(B^*)(\nu(\frac{S_1^*}{S_1})+\nu(\frac{S_1S_2^*}{S_1^*S_2}) +\nu(\frac{S_2g(B)I^*}{S_2^*g(B^*)I})+\nu(\frac{g(B^*)B}{g(B)B^*}))\\ &+S_2^*g(B^*)(\nu(\frac{h(I)B^*}{h(I^*)B})+\nu(\frac{h(I^*)I}{h(I)I^*}))\\ &+S_2^*f(I^*)(\int_0^{h}\varphi(\tau)\nu(\frac{S_2I^*f(I(t-\tau))}{S_2^*f(I^*)I})d\tau -n\nu(\frac{S_2f(I)I^*}{S_2^*f(I^*)I}))\\ &+\frac{S_2^*g(B^*)}{n}(\int_0^{h}\varphi(\tau)\nu(\frac{B^*h(I(t-\tau))}{h(I^*)B})d\tau -n\nu(\frac{h(I)B^*}{h(I^*)B}))\\ \leq &\mu_1S_1^*(2-\frac{S_1}{S_1^*}-\frac{S_1^*}{S_1}) +(\mu_2-b)S_2^*(3-\frac{S_2}{S_2^*}-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*}{S_1})\\ &+bS_2^*(2-\frac{S_1S_2^*}{S_1^*S_2}-\frac{S_1^*S_2}{S_1S_2^*}) +nS_2^*f(I^*)(\frac{f(I)}{f(I^*)}-1)(1-\frac{f(I^*)I}{f(I)I^*})\\ &+S_2^*g(B^*)((\frac{g(B)}{g(B^*)}-1)(1-\frac{g(B^*)B}{g(B)B^*}) +(\frac{h(I)}{h(I^*)}-1)(1-\frac{h(I^*)I}{h(I)I^*}))\\ \end{split} |
By Lemma 2.1, it can conclude that
\frac{dL}{dt} = \frac{dL_1}{dt}+\frac{dL_2}{dt}+\frac{dL_3}{dt}\leq 0. |
The equality \frac{dL}{dt} = 0 suggests that \frac{S_1^*}{S_1} = 1 , \frac{S_1S_2^*}{S_1^*S_2} = 1 , \frac{f(I)}{f(I^*)} = 1 and \frac{g(B)}{g(B^*)} = 1 , it implies that E_* is the maximum invariant set of system (2.1) in the set \{\frac{dL}{dt} = 0\} . Using Lemma 4.1, the endemic steady state E_* is globally asymptotically stable.
By choosing a specific kernel function and some infection functions, system (2.1) can be evolved into different dynamic models with time delays. In this section, some such examples are used to further illustrate theoretical results.
Example: Consider the kernel function \rho(\tau) = \delta(\tau-\tau_0) , where \delta is the Dirac delta function. Then system (2.1) can be rewritten as
\begin{equation} \begin{cases} \frac{dS_1}{dt} = A+bS_2-\mu_1S_1-\sigma_1S_1,\\ \frac{dS_2}{dt} = \sigma_1S_1-\alpha S_2 B-\beta S_2f(I(t-\tau_0))e^{-(\mu_2+c)\tau_0}-\mu_2S_2, \\ \frac{dI}{dt} = \alpha S_2 B+\beta S_2f(I(t-\tau_0))e^{-(\mu_2+c)\tau_0}-(\mu_2+c)I,\\ \frac{dB}{dt} = kI(t-\tau_0)e^{-(\mu_2+c)\tau_0}-dB. \end{cases} \end{equation} | (5.1) |
Where A, b, \mu_1, \mu_2, \sigma_1, \alpha, \beta, c, k and d > 0 , g(B) = \alpha B and h(I) = kI . If f(I) = \beta I , it is easy to verify that f(I), g(B) and h(I) satisfy the assumptions (H_1)-(H_3) . Therefore, the basic reproduction number of system (5.1) can be defined as
\begin{equation*} R_0 = \frac{\beta S_2^0e^{-(\mu_2+c)\tau_0}}{\mu_2+c}+\frac{k\alpha S_2^0e^{-(\mu_2+c)\tau_0}}{d(\mu_2+c)}. \end{equation*} |
In this case, if R_0\leq1 , the disease-free equilibrium is globally asymptotically stable and the endemic equilibrium is also globally asymptotically stable if R_0 > 1 (see Figures 1 and 2). As is shown in Figure 2, initial values and time delays can not change dynamic properties of system (2.1) when R_0 is greater than one or less than or equal to one. However, from Figure 2b, it is easy to find that time delay \tau_0 has a significant impact on the positive equilibrium of system (5.1).
If f(I) = \beta I^2 [28], then f(I) does not satisfy hypothesis (H_2) , that is, \frac{f(I)}{I} is a monotonically increasing function. From Figure 3, it is easy to see that the system (5.1) appears periodic oscillation behavior under certain conditions.
In this article, a general S_1S_2IB dynamics model with distributed time delay for animal brucellosis is formulated. Under the assumptions of general biological significance, the non-negative and boundedness of the solution of system (2.1) is first proved. And then the global dynamics of the steady-state solution of system (2.1) are analyzed by constructing Lyapunov functional, it is found that the dynamic properties of equilibria depend on the basic reproduction number R_0 : If R_0\leq 1 , animal brucellosis will eventually die out regardless of the initial value; and if R_0 > 1 , the spread of animal brucellosis is persistent and it eventually reaches the endemic steady state. These results imply that distributed time delay does not change the dynamic properties of system (2.1) when R_0 is greater than one or less than or equal to one. Finally, similar to numerical methods in [35,36], the stability results and other dynamic behaviors are further illustrated through numerical simulation, it turns out that the system experiences periodic oscillations if the assumption (H_2) is not satisfied. In other words, the system (2.1) may exhibit more complex dynamic behaviors if the function \frac{f}{I}, \frac{g}{B} or \frac{h}{I} is monotonically increasing. In these cases, the impact of distributed time delay on other dynamical behaviors of system (2.1) is not completely clear, then we leave these for future work.
This research is partially supported by the National Youth Science Foundation of China (11501528), the National Sciences Foundation of China(11571324) and the Fund for Shanxi "1331KIRT".
The authors declare that they have no competing interests.
[1] | Albadi El-Saadany (2008) A summary of demand response in electricity markets. Electr Power Syst Res 78: 1989–1996. |
[2] | American Council for an Energy-Efficient Economy (ACEEE). Available from: http://aceee.org. |
[3] | Baxter LW (1995) Net Lost Revenue Adjustment (NLRA) Mechanisms for Utility DSM Programs. Energy 20: 1215–1223. |
[4] | Behrangrad M (2015) A review of demand side management business models in the electricity market. Renewable Sustainable Energy Rev 47: 270–283. |
[5] | Berg W, Nowak S, Kelly M, et al. (2016) The 2016 State Energy Efficiency Scorecard. Research Report U1606, American Council for an Energy-Efficient Economy (ACEEE). |
[6] | Bertoldi P, Labanca N, Rezessy S, et al. (2013) Where to place the saving obligation: Energy end-users or suppliers? Energy Policy 63: 328–337. |
[7] | Bertoldi P, Rezessy S, Anable J, et al. (2011) Energy Saving Obligations and White Certificates: Ideas and Considerations for the Transport Sector, Int J Sustainable Transp 5: 345–374. |
[8] | Bertoldi P, Rezessy S, Lees E, et al. (2010) Energy supplier obligations and white certificate schemes: Comparative analysis of experiences in the European Union. Energy Policy 38: 1455–1469. |
[9] | Bertrand A, Mian A, Kantor I, et al. (2019) Regional waste heat valorisation: A mixed integer linear programming method for energy service companies. Energy, 167: 454–468. |
[10] | Bhattacharyya SC (2011) Energy Economics: Concepts, Issues, Markets and Governance, Springer. |
[11] | Bigliani R, Eastman R, Segalotto J, et al. (2015) Designing the new utility business models. White paper, Energy Insights IDC. |
[12] | Boait P (2009) Energy services and ESCos—their benefits and implications for regulation and the consumer, research report, De Montfort University, UK. |
[13] | Bundgaard SS, Dyhr-Mikkelsen K, Larsen AE, et al. (2013a) Energy Efficiency Obligation Schemes in the EU-Lessons Learned from Denmark. First quarter 2013, Int Assoc Energy Econ (IAEE). |
[14] | Bundgaard SS, Togeby M, Dyhr-Mikkelsen K, et al. (2013b) Spending to save: Evaluation of the energy efficiency obligation in Denmark. ECEEE Summer Study 2013 Proceedings, European Council for an Energy-Efficient Economy (ECEEE). |
[15] | Bureau of Energy Efficiency (BEE) (2015) National Mission for Enhanced Energy Efficiency, Ministry of Power, Government of India, Available from: https://www.beeindia.gov.in/content/nmeee-1. |
[16] | CAISO (California Independent System Operator). Available from: http://www.caiso.com/Pages/default.aspx. |
[17] | Capgemini (2008) Demand Response: A decisive breakthrough for Europe, in collaboration with Vaasaett and Enerdata. |
[18] | Chapman AC, Verbic G (2017) An iterative on-line auction mechanism for aggregated demand-side participation. IEEE Trans Smart Grid 8: 158–168. |
[19] | Chen Y, Nie P, Yang Y (2017) Energy management contract with subsidy. J Renewable Sustainable Energy 9: 055903. |
[20] | Cheng C (2005) Electricity demand-side management for an energy efficient future in China: Technology options and policy priorities, PhD Thesis, Massachusetts Institute of Technology, USA, 43. |
[21] | Chinnow J, Bsufka K, Schmidt AD, et al. (2011) A simulation framework for smart meter security evaluation. 2011 IEEE International Conference on Smart Measurements for Future Grids (SMFG), 1–9, 14–16th November 2011. |
[22] | Cooke D (2011) Empowering Customer Choice in Electricity Markets, Information Paper, October 2011, International Energy Agency (IEA). |
[23] | Crossley D (2010) International best practice in using energy efficiency and demand management to support electricity networks. Report 4-Scaling the Peaks: Demand Management and Electricity Networks, Australian Alliance to Save Energy, Australia. |
[24] | Crossley D (2011) Tempo Electricity Tariff—France. 11th case study article of Task XV—Network-Driven DSM, IEA DSM Programme. |
[25] | Crossley D (2013) Effective mechanisms to increase the use of demand-side resources. Global Power Best Practice Series, Regulatory Assistance Project (RAP). |
[26] | Crossley D, Maloney M, Watt G (2000) Developing mechanisms for promoting demand-side management and energy efficiency in changing electricity businesses. IEA DSM Programme, Task Ⅵ Research Report No. 3, Energy Futures Australia, August 2000. |
[27] | Darby S (2006) The effectiveness of feedback on energy consumption–a review for DEFRA, Environmental Change Institute, University of Oxford. |
[28] | Data Communications Company (DCC). Available from: https: //www.smartdcc.co.uk. |
[29] | Department for Business, Energy & Industrial Strategy (BEIS) (2017a), Electricity Demand Reduction Pilot—Interim Evaluation Findings, February 2017, UK Government. |
[30] | Department for Business, Energy & Industrial Strategy (BEIS) (2017b) Realising the potential of demand-side response to 2025: A focus on small energy users, UK Government. |
[31] | Department for Business, Energy & Industrial Strategy (BEIS) (2016) Evidence Gathering: Thermal Energy Storage (TES) Technologies, UK Government. |
[32] | Department of Energy (DoE) (2017) Federal Energy Management Program. Department of Energy, Washington, USA. Available from: https://energy.gov/eere/femp/about-federal-energy-management-program. |
[33] | Depuru SSSR, Wang L, Devabhaktuni V (2011) Smart meters for power grid: Challenges, issues, advantages and status. Renewable Sustainable Energy Rev 15: 2736–2742. |
[34] | Devine-Wright H, Devine-Wright P (2004) From demand side management to demand side participation: Tracing an environmental psychology of sustainable electricity system evolution. J Appl Psychol 6: 167–177. |
[35] | Didden MH, D'haeseleer WD (2003) Demand Side Management in a competitive European market: Who should be responsible for its implementation?. Energy Policy 31: 1307–1314. |
[36] | Dreessen T (2003) Advantages and disadvantages of the two world ESCO models: Shared savings and guaranteed savings, in: Proceedings of the First Pan-European Conference on Energy Service Companies. |
[37] | Efficiency Vermont (2016a) Available from: https://www.efficiencyvermont.com/about/history. |
[38] | Efficiency Vermont (2016b) Available from: Annual Report 2015, Efficiency Vermont, Burlington: Vermont, USA, October 2016. |
[39] | EfficiencyOne (2015). Available from: https://www.efficiencyone.ca/who-we-are/. |
[40] | El Bakari K, Myrzik JMA, Kling WL (2009) Prospects of a virtual power plant to control a cluster of distributed generation and renewable energy sources, 2009 Proceedings of the 44th International Universities Power Engineering Conference (UPEC), 1–5, 1–4th September 2009, Glasgow. |
[41] | Energy Efficiency Services Limited (EESL) (2015) Available from: http://www.eeslindia.org/User_Panel/UserView.aspx?TypeID=1025&p=About%20Us. |
[42] | ERCOT (Electric Reliability Council of Texas). Available from: http://www.ercot.com. |
[43] | Eto J, Stoft S, Belden T (1997) The theory and practice of decoupling utility revenues from sales. Util Policy 6: 43–55. |
[44] | Eto JH (1994) The theory and practice of decoupling. eScholarship, Lawrence Berkeley Laboratory, USA. |
[45] | European Commission (2015) Energy Efficiency-Energy Performance Contracting. Available from: http://iet.jrc.ec.europa.eu/energyefficiency/european-energy-service-companies/energy-performance-contracting. |
[46] | European Council Directive 2009/72/EC on smart meter rollout. |
[47] | European Union Directive 2012/27/EU on energy efficiency. |
[48] | Eyre N, Pavan M, Bodineau L (2009) Energy company obligations to save energy in Italy, the UK and France: What have we learnt?. ECEEE Summer Study 2013 Proceedings, European Council for an Energy-Efficient Economy (ECEEE). |
[49] | Fattahi A, Nahavandi A, Jokarzadeh M (2018) A comprehensive reserve allocation method in a micro-grid considering renewable generation intermittency and demand side participation. Energy 155: 678–689. |
[50] | Federal Energy Management Program (FEMP) (2016) Overview—Federal Energy Management Program. Department of Energy, Washington, USA, Available from: https://energy.gov/sites/prod/files/2016/08/f33/femp_overview.pdf. |
[51] | Fetter J, Thomas S, Potes A, et al. (2012) Energy Efficiency in the Forward Capacity Market: Evaluating the Business Case for Building Energy Efficiency as a Resource for the Electric Grid. 2012 ACEEE Summer Study on Energy Efficiency in Buildings, American Council for an Energy-Efficient Economy (ACEEE), USA. |
[52] | Gillingham K, Newell RG, Palmer K (2009) Energy efficiency economics and policy. Annual Review of Resource Economics, NBER Working Paper Series, National Bureau of Economic Research. |
[53] | Granovetter MS (1973) The strength of weak ties. Am J Sociol 78: 1360–1380. |
[54] | Haas R, Nakicenovic N, Ajanovic A, et al. (2008) Towards sustainability of energy systems: A primer on how to apply the concept of energy services to identify necessary trends and policies. Energy Policy 36: 4012–4021. |
[55] | Hamidi V, Li F, Robinson F (2009) Demand Response in the UK's domestic sector. Electr Power Syst Res 79: 1722–1726. |
[56] | Hamilton B, Plunkett J, Wickenden M (2002) Gauging success of the nation's first efficiency utility: Efficiency Vermont's first two years. 2002 ACEEE Summer Study on Energy Efficiency in Buildings—Proceedings, American Council for an Energy-Efficient Economy (ACEEE), USA. |
[57] | Hannon M, Foxon TJ, Gale WF (2013) The co-evolutionary relationship between Energy Service Companies and the UK energy system: Implications for a low-carbon transition. Energy Policy 61: 1031–1045. |
[58] | Hargreaves T, Nye M, Burgess J (2010) Making energy visible: A qualitative field study of how householders interact with feedback from smart energy monitors. Energy Policy 38: 6111–6119. |
[59] | Hayes S, Nadel S, Kushler M, et al. (2011) Carrots for Utilities: Providing Financial Returns for Utility Investments in Energy Efficiency. Report Number U111, American Council for an Energy-Efficient Economy (ACEEE), Washington DC, USA. |
[60] | Hienerth C, Keinz P, Lettl C (2011) Exploring the nature and implementation of user-centred business models. Long Range Plann 44: 344–374. |
[61] | Hills J (2012) Getting the measure of fuel poverty—final report of the Fuel Poverty Review. CASE report 72, Centre for Analysis of Social Exclusion (CASE), March 2012. |
[62] | Hirst E (2002) The financial and physical insurance benefits of price-responsive demand. Electr J 4: 66–73. |
[63] | Hirst E, Blank E (1994) Solutions to regulatory disincentives for utility DSM programmes. Util Policy 4: 105–112. |
[64] | Home Energy Conservation Act of 1996, London: HMSO, UK. |
[65] | ICF International (2007) Introduction to energy performance contracting, report prepared for the US Environmental Protection Agency (EPA), USA. |
[66] | IEA (2015) Energy Efficiency Market Report 2015, International Energy Agency, Paris: France. |
[67] | IEA DSM Programme (2008) Evaluation and Acquisition of Network-driven DSM Resources. Task 15-Network driven DSM, International Energy Agency (IEA). |
[68] | International Energy Agency Demand-Side Management Technology Collaboration Programme (IEA DSM TCP) Task 25—Business models for a more effective market update of DSM energy services. Available from: http://www.ieadsm.org/task/task-25-business-models-for-a-more-effective-uptake/. |
[69] | International Finance Corporation (IFC) (2011) IFC Energy Service Company Market Analysis, Final Report (Revised), 23rd June 2011. Available from: http://www.ifc.org/wps/wcm/connect/dbaaf8804aabab1c978dd79e0dc67fc6/IFC+EE+ESCOS+Market+Analysis.pdf?MOD=AJPERES. |
[70] | ISO-NE (ISO-New England). Available from: https://www.iso-ne.com. |
[71] | Johnson-Brown J, Reingen PH (1987) Social ties and word-of-mouth referral behaviour. J Consum Res 14. |
[72] | Kangas HL, Lazarevic D, Kivimaa P (2018) Technical skills, disinterest and non-functional regulation: Barriers to building energy efficiency in Finland viewed by energy service companies. Energy Policy 114: 63–76. |
[73] | Kim JH, Shcherbakova A (2011) Common failures of demand response. Energy 36: 873–880. |
[74] | Kushler M, York D, Witte P (2006) Aligning Utility Interests with Energy Efficiency Objectives: A Review of Recent Efforts at Decoupling and Performance Incentives. Report Number U061, American Council for an Energy-Efficient Economy (ACEEE), Washington DC, USA. |
[75] | Lise W, Kruseman G (2008) Long-term price and environmental effects in a liberalized electricity market. Energy Econ 30: 230–248. |
[76] | Liu Y (2017) Demand response and energy efficiency in the capacity resource procurement: Case studies of forward capacity markets in ISO New England, PJM and Great Britain. Energy Policy 100: 271–282. |
[77] | Lovins AB (1990) The negawatts revolution. Conference Board Magasine 27: 18–23. |
[78] | Marino A, Bertoldi P, Rezessy S (2010) Energy Service Companies Market in Europe: Status Report 2010. Joint Research Centre (JRC) European Commission, EUR 24516 EN-2010, Italy. |
[79] | Marino A, Bertoldi P, Rezessy S, et al. (2011) A snapshot of the European energy service market in 2010 and policy recommendations to foster a further market development. Energy Policy 39: 6190–6198. |
[80] | Martínez Ceseña EA, Good N, Mancarella P (2015) Electrical network capacity support from demand side response: Techno-economic assessment of potential business cases for small commercial and residential end-users. Energy Policy 82: 222–232. |
[81] | Mau P, Eyzaguirre J, Jaccard M, et al. (2008) The 'neighbour effect': Simulating dynamics in consumer preferences for new vehicle technologies. Ecological Economics 68: 504–516. |
[82] | MISO (Midcontinent Independent System Operator). Available from: https://www.misoenergy.org/Pages/Home.aspx. |
[83] | Molina M, Kushler M (2015) Policies Matter: Creating a Foundation for an Energy-Efficient Utility of the Future. ACEEE White Paper, American Council for an Energy-Efficient Economy (ACEEE), Washington DC, USA. |
[84] | Mourik R, Bouwknegt R (2017) Mind your business: Entrepreneurs, their dynamic capabilities, context and new business models for energy efficiency services. ECEEE Summer Study 2017 Proceedings, European Council for an Energy-Efficient Economy (ECEEE). |
[85] | Murray J (2015) First government auction pays companies to power down. businessGreen, 4th February 2015. Available from: http://www.businessgreen.com/bg/news/2393653/first-government-auction-pays-companies-to-power-down-at-peak-times. |
[86] | Narayanamurthy R, Robinson J (2015) Bring Your Own Device Program Approaches. Palo Alto, Electric Power Research Institute (EPRI), California, USA. |
[87] | Natural Resources Defense Council (NRDC) (2013) Decoupling Across the United States. Map, August 2013. Available from: http://www.nrdc.org/energy/decoupling. |
[88] | Natural Resources Defense Council (NRDC) (2016) The Evidence Is In: Decoupling Spurs Energy Efficiency Investment. Available from: https://www.nrdc.org/experts/samantha-williams/evidence-decoupling-spurs-energy-efficiency-investment. |
[89] | Natural Resources Defense Council (NRDC). Available from: https://www.nrdc.org. |
[90] | Newborough M, Augood P (1999) Demand-side management opportunities for the UK domestic sector. IEE Proceedings—C Generation, Transmission and Distribution 146: 283–293. |
[91] | Nilsson H, Bangens L, Goven B, et al. (2012) We are lost if we don't develop new business models. ECEEE Summer Study 2012 Proceedings, European Council for an Energy-Efficient Economy (ECEEE). |
[92] | Nowak S, Baatz B, Gilleo A, et al. (2015) Beyond carrots for utilities: A national review of performance incentives for energy efficiency. Research Report U1504, American Council for an Energy-Efficient Economy (ACEEE), June 2015. |
[93] | NYISO (New York Independent System Operator). Available from: http://www.nyiso.com/public/index.jsp. |
[94] | Oikonomou V, Di Giacomo M, Russolillo D, et al. (2012) White Certificates in the Italian Energy Oligopoly Market. Energy Sources Part B-Economics Plann Policy 7: 104–111. |
[95] | Oikonomou V, Jepma C, Becchis F, et al. (2008) White Certificates for energy efficiency improvement with energy taxes: A theoretical economic model. Energy Econ 30: 3044–3062. |
[96] | Oikonomou V, Patel MK, van der Gaast W, et al. (2009) Voluntary agreements with white certificates for energy efficiency improvement as a hybrid policy instrument. Energy Policy 37: 1970–1982. |
[97] | Otashu JI, Baldea M (2018) Grid-level "battery" operation of chemical processes and demand-side participation in short-term electricity markets. Appl Energy 220: 562–575. |
[98] | Owen G, Pooley M, Ward J (2012) What demand side services could household customers offer?. Paper 3 of the GB Electricity Demand—realising the resource project, Sustainability First, April 2012, UK. |
[99] | Owens S, Driffil L (2008) How to change attitudes and behaviours in the context of energy. Energy Policy 36: 4412–4418. |
[100] | Pachauri S, Ürge-Vorsatz D, LaBelle M (2012) Synergies between energy efficiency and energy access policies and strategies. Global Policy 3: 187–197. |
[101] | Parrish B, Heptonstall P, Gross R (2016) The potential for UK residential demand-side participation. System Architecture Challenges: Supergen + for HubNet, HubNet. |
[102] | Paterakis NG, Erdinç O, Catalão PS (2017) An overview of demand response: Key elements and international experience. Renewable Sustainable Energy Rev 69: 871–891. |
[103] | PJM (2017) Demand Response Fact Sheet. Available from: https://www.pjm.com/~/media/about-pjm/newsroom/fact-sheets/demand-response-fact-sheet.ashx. |
[104] | PJM, Available from: http://www.pjm.com. |
[105] | Poudineh R, Jamasb T (2014) Distributed generation, storage, demand response and energy efficiency as alternatives to grid capacity enhancement. Energy Policy 67: 222–231. |
[106] | Prüggler N (2003) Economic potential of demand response at household level—are Central-European market conditions sufficient?. Energy Policy 60: 487–498. |
[107] | Regulatory Assistance Project (RAP), Available from: http://www.raponline.org. |
[108] | Rosenow J (2012) Energy savings obligations in the UK-A history of change. Energy Policy 49: 373–382. |
[109] | Rosenow J, Galvin R (2013) Evaluating the evaluations: Evidence from energy efficiency programmes in Germany and the UK. Energy Build 62: 450–458. |
[110] | Rosenow J, Platt R, Flanagan B (2013) Fuel poverty and energy efficiency obligations—A critical assessment of the supplier obligation in the UK. Energy Policy 62: 1194–1203. |
[111] | Saker D, Coker P, Vahdati M, et al. (2013) Unlocking the demand response potential from domestic hot water tanks. 4th Annual TSBE (Technologies for Sustainable Built Environments) EngD Conference Proceedings, University of Reading, UK. |
[112] | Satchwell A, Cappers P, Goldman C (2011) Carrots and sticks: A comprehensive business model for the successful achievement of energy efficiency resource standards. Util Policy 19: 218–225. |
[113] | Schultz D, Eto J (1990) Carrots and sticks: Shared-savings incentive programs for energy efficiency. Electr J 3: 32–46. |
[114] | Spees K, Lave LB (2007) Demand response and electricity market efficiency. Electr J 20: 69–85. |
[115] | Steinberger JK, van Niel J, Bourg D (2009) Profiting from negawatts: Reducing absolute consumption and emissions through a performance-based energy economy. Energy Policy 37: 361–370. |
[116] | Stoker L (2017) Electricity Demand Reduction pilot reports positive interim results. Clean Energy News, 10th February 2017. Available from: http://www.cleanenergynews.co.uk/news/efficiency/electricity-demand-reduction-pilot-reports-positive-interim-results. |
[117] | Strachan N, Warren P (2011) Incorporating behavioural complexity in energy-economic models. UK Energy Research Centre Conference on: Energy and People: Futures, Complexity and Challenges, 20–21 September 2011, Environmental Change Institute, Oxford, UK. |
[118] | Strbac G (2008) Demand side management: Benefits and challenges. Energy Policy 36: 4419–4426. |
[119] | Stromback J, Dromacque C, Yassin MH (2011) The potential of smart meter enabled programs to increase energy and systems efficiency: A mass pilot comparison, Vaasaett, 13. |
[120] | The Electricity Journal. Available from: https://www.journals.elsevier.com/the-electricity-journal. |
[121] | THINK (2012) How to refurbish all buildings by 2050. Topic 7, Final report, THINK, June 2012. |
[122] | Titus E, Michals J, Hurley D, et al. (2009) Energy efficiency as a resource in the PJM Capacity Market. Int Energy Program Eval Conf 2009. |
[123] | Torriti J, Hassan MG, Leach M (2010) Demand response experience in Europe: Policies, programmes and implementation. Energy 35: 1575–1583. |
[124] | United Nations/Framework Convention on Climate Change (2015) Adoption of the Paris Agreement, 21st Conference of the Parties, Paris: United Nations. |
[125] | Ürge-Vorsatz D, Hauff J (2001) Drivers of market transformation: Analysis of the Hungarian lighting success story. Energy Policy 29: 801–810. |
[126] | Utilities Policy: https://www.journals.elsevier.com/utilities-policy. |
[127] | Vermont Public Service Department (2014) Energy Efficiency. Available from: http://publicservice.vermont.gov/topics/energy_efficiency. |
[128] | Walawalkar R, Fernands S, Thakur N, et al. (2010) Evolution and current status of demand response (DR) in electricity markets: Insights from PJM and NYISO. Energy 35: 1553–1560. |
[129] | Ward J, Pooley M, Owen G (2012) What demand-side services can provide value to the electricity sector?. Paper 4, GB Electricity Demand—Realising the Resource, Sustainability First, June 2012. |
[130] | Warren P (2014) The use of systematic reviews to analyse demand-side management policy. Energy Effic 7: 417–427. |
[131] | Warren P (2014b) A review of demand-side management policy in the UK. Renewable Sustainable Energy Rev 29: 941–951. |
[132] | Warren P (2015) Demand-Side Management Policy: Mechanisms for Success and Failure, PhD Thesis, University College London (UCL), UK. |
[133] | Warren P (2017a) The potential of smart technologies and micro-generation in UK SMEs. Energies 10. |
[134] | Warren P (2017b) Transferability of demand-side policies between countries. Energy Policy 109: 757–766. |
[135] | Warren P (2018) Demand-side policy: Global evidence base and implementation patterns. Energy Environ 29: 706–731. |
[136] | Warren P (2019) Demand-side policy: Mechanisms for success and failure. Econ Energy Environ Policy. |