On the extinction of continuous-state branching processes in random environments

Abstract: This paper establishes a model of continuous-state branching processes with time inhomogeneous competition in Lévy random environments. Some results on extinction are presented, including the distribution of the extinction time, the limiting distribution conditioned on large extinction times and the asymptotic behavior near extinction. This paper also provides a new timespace transformation which can be used for further exploration in similar models. The results are applied to an epidemic model to describe the dynamics of infectious population and a virus model to describe the dynamics of viral load.


Introduction
Branching processes have wide applications and have raised the interest of many researchers, see [10]. The model was first established by Galton and Watson [8] in order to study the demographic phenomenon that a large proportion of families were continually dying out. Decades later, Feller [7] formalized the model to the modern version. He provided the recursion of the population of the n-th generation Z n as follows, where ξ n,i are random variables representing the number of offspring of the i-th individual in the n-th generation. In this model, time is discretized by generations. However, in real life, population develops continuously in time. Therefore, researchers improved the model to a continuous-time version. The first general formulation of the continuous-time discrete-state branching processes was given by Kolmogorov and Dmitrive, and has been one of the most popular demographic models, see [21]. In microbiology, the population is large and the events of birth and death occur very frequently. Scientists usually use terminologies like density or viral load to describe the amount, and the number is often not integers, see [20]. In this situation, a branching model of continuous-time and continuousstate is more appropriate and convenient to approximate the dynamics. Ba and Pardoux [1] proved that the continuous-state branching processes can be considered as the scaling limit of the continuoustime discrete-state branching processes. They considered the interactions between individuals and established the limiting stochastic equation, (1.2) where x → f (x) is a function on [0, ∞) representing the interactions, c is a positive constant and B(s) is a standard Brownian motion. This has the general form of a Feller branching diffusion.
There are plenty of biological articles focusing on how the population is influenced by the environment, see [5] for an example. But the investigation on the mathematical models of the continuous-state branching process in random environment, which is usually abbreviated as CBRE-process, is just in its infancy. Pioneering works on CBRE-processes include [2] and [9]. To describe the drastic reduction in population by catastrophes, Bansaye et al. [2] added a Lévy process with bounded variation into the continuous-state branching model. They used the following equation to model infected cells, where c > 0, b are constants, B(s) is a Brownian motion and N(ds, dθ) is Poisson random measure representing the catastrophes. They studied the behavior of P(Y t > 0) as t → ∞. This work was generalised by Li and Xu [17], who set the random environment to be a Lévy process {L(t) : t ≥ 0} with no jump less than −1. This jump restriction prevents the population from becoming negative. They constructed an equation for (α + 1)-stable branching processes in Lévy random environment (CBRE-processes) as follows, where c > 0, α ∈ (0, 1] are constants, {Z α (t) : t ≥ 0} is a spectrally positive (α + 1)-stable process. Motivated by the models above, we construct an (α + 1)-stable continuous-state branching process with both the effect of random environments and a time inhomogeneous competition, which will be abbreviated as the (α + 1)-stable CBCRE-process. The equation of our model is where s → f (s) is a locally bounded non-negative function on [0, ∞) describing the competition between individuals. When f ≡ 0, our model degenerates to the model of (α + 1)-stable CBREprocesses in [17]. Furthermore, [17] studied the asymptotic behavior of the survival probability.
But our results include the asymptotic behavior of the process itself. In other words, we study the asymptotic property of the paths instead of probability. When α = 1, f ≡ 0 and {L(t) : t ≥ 0} takes the special form of t 0 1 0 (1 − θ)N(ds, dθ), our model degenerates to the model of [2]. When α = 1 and L(t) ≡ 0, it becomes a model with only interactions but no environments, like [1].
In this paper, we define the extinction time of the process as the first time that the process reaches zero, which is denoted by T 0 . We provide the distribution of T 0 , and give a necessary condition for T 0 to be almost surely finite. The long-time behavior is an important topic for stochastic processes, see [11][12][13]. Our results also include the long-time behavior of the model. We present the limit distribution of CBCRE-processes conditioned on large extinction times. Furthermore, we find a determining function g to describe the asymptotic behavior near extinction of the (α+1)-stable CBCRE-process {X(t) : t ≥ 0}. In other words, This is a result more accurate than [2] and [17], since we describe the behavior of the paths rather than the extinction and survival probability. And this result is new even in CBRE-processes without competition. The paper is organized as follows, in Section 2, we establish the model, present the stochastic equation and state the main results. In Section 3, we give the proofs of the main results including a time-space transformation. In Section 4, we apply the results to two mathematical biological models. In the first part of Section 4, we improved the stochastic Susceptible-Infected-Recovered (SIR) epidemic model by considering the influence of environment. In the second part of Section 4, we apply our results to a virus elimination model and conclude that the differences between the immune systems of different hosts only affects the distribution of their recovery time. The asymptotic behavior of the viral load near the host's recovery time (equivalently, the extinction time of virus) depends only on the characteristics of the virus.

Preliminaries and main results
In this section, we establish the model of (α + 1)-stable CBCRE-processes on a filtered probability space (Ω, F , F t , P). Suppose that constants b ∈ R, c > 0, α ∈ (0, 1]. A stable branching mechanism is a continuous function φ on [0, ∞) with the representation is an (F t )-Lévy process with Lévy-Itô decomposition: where {W(t) : t ≥ 0} is an (F t )-Brownian motion, and N(ds, dz) is an (F t )-Poisson random measure on (0, ∞)×R with intensity dsn(dz), whose compensated measure is denoted byÑ(ds, dz). Let s → f (s) be a locally bounded nonnegative function on [0, ∞). For any positive constant x, consider the following stochastic integral equation, where {Z α (t) : t ≥ 0} is a spectrally positive (α + 1)-stable process with α ∈ (0, 1]. When α = 1, this can be regarded as a Brownian motion. When α ∈ (0, 1), its Lévy measure is In this equation, t 0 f (s)X(s)ds describes the influence of competition on the population dynamics. We may consider f (s) as a competition rate. It changes with time and it stays non-negative, since at a fixed time point, the competition becomes fiercer as the population becomes larger. In Eq (2.3), there are two kinds of noise: (1 + α)cX(s−)dZ α (s) represents the randomness of branching, which is the same as the noise term in the equation of traditional stable branching processes, see P.55 of [16]; and t 0 X(s−)dL(s) represents the effect of random environments, see [9] and [17].
For convenience, define another (F t )-Lévy process {ξ(t) : t ≥ 0} by In the situation where f ≡ 0, the unique strong solution of (2.3) is the CBRE-process with branching mechanism φ and random environment {ξ(t) : t ≥ 0} defined by [9]. Denotē Proof. According to Theorem 1.3 of [6], has a unique strong solution. Following by similar arguments in [19], the existence and uniqueness of Eq (2.3) can be proved. The remaining procedures can be obtained by repeating the calculations in the proof of Theorem 3.6 in [9], with ξ in there replaced byξ.
The unique strong solution {X(t) : t ≥ 0} of Eq (2.3) is a CBCRE-process with initial value x, branching mechanism φ, and competition f in random environment ξ, which is abbreviated as (α + 1)stable CBCRE-process. By the definition in He et al. [9], the extinction time T 0 := inf{t ≥ 0 : X(t) = 0} is obviously an F t -stopping time. Again, by replacing parameters as in the proof of Theorem 3.3, we have the following conclusion paralleling to that in Section 5 of [17].
Proposition 2.2. The distribution of the extinction time T 0 is provided as Therefore, lim inf t→∞ξ (t) = −∞ is the necessary condition for the stable CBRE-process {X(t) : t ≥ 0} to extinct in a finite time almost surely. The main result of this paper is the extinction behavior of the process, including the limit distribution conditioned on large extinction times and the asymptotic behavior near extinction.

Theorem 2.3. [Conditional Limit
Theorem] Let x > 0 and t ≥ 0. {X(t) : t ≥ 0} is a stable CBCREprocess with initial value x and extinction time T 0 . Then for any F t -measurable bounded random variable F,P whereP is a probability measure such thatP The above result provides the distribution of {X(t) : t ≥ 0} conditioned on large extinction times. Under this new probability measure, {X(t) : t ≥ 0} is still a Markov process in (0, ∞). See [15] for the similar conditional limit theorems of continuous-state branching processes with neither environment nor competition.

Proofs of main results
In this section, we prove the main results and provide an interesting time-space transformation.
is the strong solution of (2.3). By Itô's formula, It is easy to see that {M t : t ≥ 0} is a positive martingale with PM t = x. The conclusion is obvious.
Proof of Theorem 2.3. By Proposition 3.1 and Girsanov's theorem, there is a unique probability measureP on (Ω, F ) such thatP for any F t -measurable bounded random variable F. Now it suffices to show that P[e −λX(t) ] = lim r→∞ P[e −λX(t) |T 0 > r + t]. Define P¯ξ to be the quenched law conditioned on the random environment, It is not difficult to get v¯ξ 0,r+t = v¯ξ 0,t (v¯ξ t,t+r ) and d dλ v¯ξ 0,t (λ) λ=0 = e¯ξ t . By the above and (2.13) we have, lim r→∞ v¯ξ t,t+r = 0 almost surely. Then By bounded convergence theorem, it is easy to see that That completes the proof. Consider the stochastic integral equation The unique strong solution of (3.2) is a continuous-state branching process with (α+1)-stable branching mechanism. We denote the process by {Y(t) : t ≥ 0}. There is a nice relationship between the process {Y(t) : t ≥ 0} and {X(t) : t ≥ 0}.   A(s)). ThenZ is also an (α + 1)-stable process with Lévy measure m(dz) defined by (2.4). Thus, We denote the extinction time of {Y(t) : t ≥ 0} by T 0 := inf{t ≥ 0 : Y(t) = 0}. There is a result on the asymptotic behavior near extinction of Y(t), which we state as the Lemma below. [14]] Suppose that {Y(t) : t ≥ 0} is a continuous-state branching process with branching mechanism φ given by (2.1). Then
Notice that the limitation is in the sense of P-a.s., and αξ(T 0 ) is pathwisely constant. Thus, when t → 0+, the influence ofξ(T 0 ) can be ignored compared to log 1 t . Hence we get the conclusion. The time-space transformation that we use in the proof builds a strong bridge between CBCREprocesses and traditional continuous-state branching processes. This novel approach contributes to the study of continuous-state branching processes in random environments and can be used for further exploration in similar models.

Applications
In this section we apply our model to mathematical biology. We consider a particle system with X(0) particles at time t = 0. After each unit of time, each particle in the system splits into a random number of new particles following a fixed reproduction law. This procedure can approximate the reproduction of a virus at the micro level, as well as the increase in the infectious population at the macro level. If we denote the population of the n-th generation by Z n , we can get where ξ n,i are random variables representing the number of offspring of the i-th individual in the nth generation. This is the definition of classic branching processes, or GW-processes. However, the discrete setting is idealized. This is problematic since the reproduction of the virus and the speed of contagion are always quite fast, and the quantities we discuss are always quite large. Thus, it is more appropriate and convenient to compress the original process defined by (4.1) both in time and state. It was provided in 36 and 37 pages of [4] that the scaling limits of a series of GW-processes with a fixed reproduction law are continuous-time and continuous-state stable branching processes. Moreover, when we consider the random environment and competition, we can have stable CBCRE-processes as the scaling limits, see [3]. Thus, it is reasonable to apply our model to the study of epidemic and viral load.

Application to an epidemic model
In this part, we apply our results to an improved stochastic Susceptible-Infected-Recovered (SIR) epidemic model. The stochastic SIR-model was introduced in [18] as the nonlinear stochastic differential equations: In this model, the whole population N(t) is divided into three parts: the susceptible part S (t), the infectious part I(t) and the recovered part R(t). Parameters α, β and γ denote the disease mortality rate, the transmission rate and the recovery rate, respectively. B 1 (t), B 2 (t) and B 3 (t) are mutually independent standard Brownian motions.
To better approximate the actual conditions, we consider environmental influences on the epidemic model. The environment may affect the dynamics of an infectious population via fast traffic tools like airplanes and high-speed railways, and policies like quarantine and lock-down, among others. We describe the environment by the Lévy processes ξ(t) and L(t) defined in Section 2 and build the following equations: When the infectious population is small compared to the large whole population, the infectious population can be approximated by a branching process defined by where B(s) is a standard Brownian motion and L(t) is defined by (2.2). According to the results in previous sections, if the condition lim inf t→∞ [ξ(t) + (β − α − γ)t] = −∞ is satisfied, the epidemic will be eliminated in a finite amount of time T. The distribution of T is given by Moreover, the infectious population will tend toward zero as t → T in the behavior described by the function (α + β + γ)(T − t) 2 log log 1 .

Application to a virus elimination model
The continuous-state branching process is a natural model for population which die and multiply fast and randomly, see Bansaye et al. [2]. Assume that at time 0 a person is infected with initial viral load x. Then the viruses in the person replicate with the branching mechanism φ(λ) = cλ α+1 , λ > 0, where c > 0, α ∈ (0, 1]. The immune cells of the host are produced according to a Poisson distribution with parameter r > 0, and kill a random fraction θ ∈ (0, 1) of the viruses. The distribution of the killing fraction θ is denoted by ν, where ν satisfies In this equation, {Z α (t) : t ≥ 0} is a spectrally positive (α + 1)-stable process with α ∈ (0, 1], as defined in Section 2. M(ds, dθ) is a Poisson random measure with intensity rdsν(dθ). A host recovers when the viral load reaches zero. And the time it takes to recover is denoted by T, which amounts to the extinction time of {X(t) : t ≥ 0}. According to the conclusions in Section 2, T is a random time with distribution P(T ≤ t) = P(e −xη 0,t ), Moreover, the viral load tends toward zero near the recovery moment T in a behavior that can be approximated by function (cα(T − t)) α −1 (log log 1 (T − t) ) −α .
These results demonstrate • The fewer viruses with which a host was infected at beginning, the faster the host is expected to recover; • Differences in the sensitivity and effectiveness between the immune systems of different hosts (corresponding to r and ν in the model) result in different lengths of recovery time; • The proliferative ability of the virus (corresponding to c, α in the model) affects both the recovery time and the behavior of the viruses being eliminated. Indeed, it is the only factor determining this behavior.