Jump-Diffusion Models for Valuing the Future: Discounting under Extreme Situations

We develop the process of discounting when underlying rates follow a jump-diffusion process, that is, when, in addition to diffusive behavior, rates suffer a series of finite discontinuities located at random Poissonian times. Jump amplitudes are also random and governed by an arbitrary density. Such a model may describe the economic evolution, specially when extreme situations occur (pandemics, global wars, etc.). When, between jumps, the dynamical evolution is governed by an Ornstein–Uhlenbeck diffusion process, we obtain exact and explicit expressions for the discount function and the long-run discount rate and show that the presence of discontinuities may drastically reduce the discount rate, a fact that has significant consequences for environmental planning. We also discuss as a specific example the case when rates are described by the continuous time random walk.


Introduction
The importance of discounting, particularly in the long run, does not exclusively refer to finance, but to many aspects of global economy. This is the case of long-term environmental planning, which is certainly acute in climate action. In essence, an environmental disaster that could cost X to repair at a future time t is worth exp(−rt)X today, where r is the interest rate assuming that is continuously compounded. This simple analysis assumes that interest rates remain constant between today and the distant future t, which may be decades ahead. The rate r thus becomes a key magnitude to decide whether it is more beneficial to take action today with a significant investment or whether the discount gives negligible value to today's investment.
No wonder that the estimation of the discount rate has enormous consequences and has been the object of intense work and controversy for a long time. While, for instance, the reputed and highly influential British economist Nicholas Stern [1] had been using a discounting rate of 1.4%, William Nordhaus [2] did propose a discount rate of 4% and even a higher rate (6%) [3]. The two estimates stand for very different perspectives on how to face climate change and other catastrophic events. Stern's rate entails for an immediate spending, and Nordhaus's figures say that action could be less urgent and not that strong. The choice of discount rate lays at the center of the debate on the urgency of climate change mitigation.
When we refer to climate, discount rates choice is based on ethical grounds [4,5] and on future economic growth assumptions. Economic arguments involve the maximization of utility functions [6]. Discounting not only includes economic growth; it also includes behavioral aspects such as impatience or the possibility of having a declining marginal utility. All of these aspects are covered in the Ramsey formula [7], which constitutes the basis of the more traditional approach to discounting [8]. It is, however, not realistic to represent discounting by deterministic functions of time. To consider decreasing exponentials with fixed rates is then too simplistic and it is necessary to consider an average overall interest rate path. In global warming, this issue becomes particularly sensitive as one must consider costs and benefits for long time horizons. Quantitative finance has indeed provided a robust framework (e.g., the so called Heath-Jarrow-Morton framework) [9], while several works have long recognized that interest rates must be modeled as random processes [10][11][12][13][14].
An economist working with statistics might decide to compute the discount rate as the average over empirical interest rates during the last 200 years (which is 2.7% in stable countries [15,16]), or estimate the average of Wall Street forward looking models with the 30 year bond prices. However, our recent work [17] shows that the rate is considerably below these averages, and this can be attributed to historical fluctuations. As a result, any proper analysis must take into account fluctuations in the real interest rate (obtained by subtracting nominal rates from inflation), which are fundamentally due to fluctuations in economic growth [18][19][20][21][22].
The function r(t) can, in principle, be described by any random process. Markovian processes are the simplest and most common hypothesis and they consist of continuous sample paths. In other words, real rates r(t) are modeled as diffusion processes. Throughout this paper, we will assume the so-called "Local Expectation Hypothesis", in which there is no market price of risk (investors are assumed to be risk neutral) and rates are based on the data generating measure [23][24][25].
We have been analyzing these issues by assuming three of the most popular stochastic models for the dynamics of interest rates [26]: Ornstein-Uhlenbeck [27], Feller [28], and lognormal [29] processes, which are also relevant in the context of statistical physics. However, we are interested in real rates which can be negative even during prolonged periods of time [15,17]: recall, real rates are nominal rates corrected by inflation, that is, r(t) = n(t) − i(t), where n(t) are nominal rates (usually positive, but not always) and computed from government bonds, and i(t) is the rate of inflation constructed out of consumer price indexes. The Ornstein-Uhlenbeck (OU) model is the only one that allows for negative rates while still considering simplest (and linear) mean reversion towards a normal interest rate. Moreover, discount asymptotic expression thus has an exponential decay with a long-run rate r ∞ that differs from historical average interest rates by being substantially smaller, zero or eventually negative [15][16][17]26,30].
We will go one step further and assume that, in addition to diffusive and continuous behavior, the sample paths of real rates r(t) also exhibit discontinuities. That is, we will model rates by a jump-diffusion process. During the last fifty years, jump-diffusion models have been extensively used not only in many branches of statistical mechanics and condensed matter physics [31], but also in economics and finance [32,33]. Specially relevant here is the work by Ahn and Thompson [13], which in 1988 already added jumps to a quite general diffusion model to investigate the effect of discontinuities on interest rates and highlighted the effects with the Feller process [28] (also called Cox-Ingersoll-Ross model [23]). Thus, the economic evolution is known to occasionally have sudden bursts that hardly adjust to continuous diffusion-like evolution. The fact that these discontinuities do not occur frequently sustains in a great measure the use of diffusion models for the economic evolution. However, many empirical observations of economic time series tend to show the appearance of many outliers in which changes of great magnitude occur during small intervals of time, in opposition to the basic diffusive hypothesis for which changes during short intervals of time are only by small amounts.
One very recent example is provided by the COVID-19 pandemics, where prices dropped worldwide approximately 40% in less than 3 weeks. Pandemic episodes are rather recurrent, thus during the 20th century, there have been reported several pandemic incidents from the Spanish flu of 1918 to the Hong Kong flu in 1968 or AIDS starting in 1981. Several similar episodes are reported during the 19th century (cholera, etc.). One could approximately quantify the appearance of 3 to 5 of such episodes per century. Other "nightmare scenarios" [34,35] of environmental disasters would include climate change, biotechnology, asteroid impacts, runaway computer systems, and nuclear proliferation, among others (see Refs. [36,37] for a thorough discussion).
The main objectives of this work are to provide a general framework and to elucidate the effect on the long-run discount-and, hence, on how we should value the future-of discontinuities that reflect the existence of high-risk events. For the sake of completeness, we consider not only sudden negative bursts, but also positive bursts as the two show distinctive behaviors. Section 2 has two parts. Section 2.1 provides the main and broad definitions being used for the analysis, and Section 2.2 presents the specificities of the Ornstein-Uhlenbeck jump-diffusion process to model interest rates. The first paragraphs in Section 3 provide the main result. Section 3.1 presents the asymptotic discount function jointly with the longrun discount rate. Two specific jump distributions are studied in Section 3.2: fixed jump amplitudes and Laplacian jump amplitudes. Before a final discussion and the conclusions, Section 3.3 finally provides the discount within the continuous time formalism, that is: a purely discontinuous process.

Main Definitions
Suppose that r(t) is a random process representing the dynamical evolution of real rates. If we define the cumulative process the discount function is then defined as where the average is taken over all possible realizations of r(t) and market price of risk is assumed to be zero (see our previous work [16]). "Local Expectation Hypothesis" in which there is no market price of risk (investors are assumed to be risk neutral) is considered [23][24][25]. Closely related to the discount function, D(t) is the (average) discount rate defined as so that the discount function can be written in the standard exponential form D(t) = exp(−td(t)). Moreover, in terms of d(t), we can define the long-run discount rate, r ∞ , as the asymptotic value of the discount rate d(t) → r ∞ as t → ∞. That is, When introducing specific stochastic models, it then becomes particularly useful to consider the bidimensional process x(t), r(t) and denote by p(x, r, t|x 0 , r 0 , t 0 ) the probability density function (PDF) of such process (sometimes referred to as the data generating measure). This PDF is defined as p(x, r, t|x 0 , r 0 , t 0 )dxdr = Prob x ≤ x(t) < x + dx, r ≤ r(t) < r + dr|x(t 0 ) = x 0 , r(t 0 ) = r 0 , and the discount function defined as the average (2) can therefore be written as The joint characteristic function of the bidimensional process (x(t), r(t)) is defined as the Fourier transform of the joint PDF: Once we know the joint characteristic function obtaining discount is straightforward. Comparison of Equations (5) and (6) shows that and obtaining the discount function is equivalent to knowing the joint characteristic function of the bidimensional return process.

Diffusion Process in the Presence of Random Jumps
The process r(t) can be any random process, although the simplest and most usual assumption consists in modeling r(t) as a diffusion process, that is, a Markovian process with continuous sample paths. We here take the further step of assuming that r(t) is a compound process that combines an ordinary diffusion with random jumps. The diffusion process trajectory thus exhibits discontinuities at random instants of time. Discontinuities will here be assumed to be finite. The resulting bidimensional process (x(t), r(t)) is described by the following pair of stochastic differential equations (all stochastic differential equations are interpreted in the sense of Itô): where f (r) and g(r) are given functions (the drift and noise intensity, respectively), ξ(t) is a zero-mean Gaussian white noise (ξ(t)dt = dW(t) where W(t) is the standard Wiener process with unit variance), and n(t) is a white shot noise. The white shot noise can be written as [38,39] where γ j and t j (j = 1, 2, 3 . . . ) are independent and identically distributed random variables. The random quantities γ i characterize the size of jumps and are described by a given PDF, which we denote by h(u). For simplicity, the size of these discontinuities are taken to be identically distributed and independent of each other as well as independent of the instants of time at which they occur. We further assume that these random times form a Poisson set of events. In such a case, the time interval τ between two consecutive jumps {t j , t j+1 } is governed by the PDF [40] where λ > 0 is the rate of the Poisson process and λ −1 is the average time interval between two consecutive jumps.
To obtain the discount we need to look at the joint PDF p(x, r, t|x 0 , r 0 , t 0 ) of the bidimensional process (x(t), r(t)). It is, however, convenient to first consider the jump PDF characterizing the discontinuities of the return process. This density is defined as [31] W(x, r|x 0 , r 0 ) = lim p(x, r, t 0 + ∆t|x 0 , r 0 , t 0 ) , A standard reasoning-based on the Chapman-Kolmogorov equation and detailed in Gardiner's monograph [31] (see also [40]) shows that the PDF of the bidimensional jump-diffusion process defined in (8) obeys the integro-differential equation with the initial condition The assumptions made above about discontinuities allow us to obtain a a more explicit expression for both the transition density W and the integro-differential Equation (11). Let us recall that the magnitude of the discontinuities, expressed by the random variables γ i , is independent of the times t i where jumps occur. We see from the model expressed by Equation (8) that the instantaneous jumps only affect r(t), but not x(t). These considerations, along with the Poisson character of jump times, allow us to take as transition density the following expression [40]: Substituting this simpler expression for W into Equation (11) and taking into account the homogeneity of both x and t (which amounts to take x 0 = 0 and t 0 = 0) as well as the normalization condition on the PDF h(u), we see that the integro-differential equation for p(x, r, t|r 0 ), Equation (11) reads ∂p ∂t and the initial condition is Equation (14) is the most general formulation of the discount problem of timehomogeneous diffusion with independent and Poissonian random jumps. In order to proceed further we need to further specify a particular diffusion process for the continuous part of the return.

The Ornstein-Uhlenbeck Process and Poissonian Jumps
In the modeling of financial interest rates, the Ornstein-Uhlenbeck (OU) diffusion process was proposed by Oldrich Vasicek during the late nineteen seventies [10]. The model allows for both positive and negative rates and is, therefore, suitable for describing the so-called real interest rates. We have extensively used this process in the study of long-run discounting [15][16][17]26]. For the OU process, the drift is linear and the noise intensity constant: The parameter m (usually referred to as "normal level") is the mean value to which the process reverts in the long run, α > 0 is the strength of the reversion to the mean, and k > 0 is the amplitude of the fluctuations. In the stationary regime when t α −1 rates are explicitly given by [26] where ξ(t) is the Gaussian white noise defined above. We, thus, easily see that the normal level m is the stationary mean value of the return, while the stationary autocorrelation function C(τ) is given by [26] showing that α −1 is the autocorrelation time and σ 2 = k 2 /2α is the variance. For this continuous model, we have been able to obtain a closed expression for the discount function D(t), which in the long run, as t → ∞ (cf. Equation (4)), reads [26] Let us now assume that the rate process r(t) is governed by an OU process with random discontinuities described by Poissonian jumps. The integro-differential equation for the joint PDF p(x, r, t|r 0 ) will be given by (cf. Equations (14) and (16)) with the initial condition Fourier transforming Equations (18) and (19) results in a much simpler problem for the characteristic function, wherep =p(ω 1 , ω 2 , t|r 0 ) is the joint Fourier transform defined in Equation (6) and is the characteristic function of the jump PDF h(u). The initial condition is now given bỹ Equation (20) is a partial differential equation of first order whose solution can be obtained by the method of characteristics [41]. In the Appendix A, we show that the exact solution to the initial-value problem (20) and (22) is given bỹ and Looking at Equation (23), we see that when there are no discontinuities (i.e., λ = 0), the PDF (23) reduces to a Gaussian density as we had obtained in previous works [26]. Denoting this density by p (0) and setting λ = 0 in Equation (23), we get We can thus write Equation (23) as Let us finally recall that knowing the joint PDF of the two-dimensional process (x(t), r(t)), the distribution of the return r(t) is given by the marginal density, and the characteristic function of return,p(ω 2 , t|r 0 ), is simply obtained by setting ω 1 = 0 in the joint characteristic function. From Equations (23)-(27), we get In terms of the characteristic functionp(ω 2 , t|r 0 ), the moments of the return, r n (t) , are given by the derivatives i np(n) (ω 2 = 0, t|r 0 ) (n = 1, 2, . . . ). Thus, for example, the variance of the return, that is the volatility, is given by where µ is the second moment of the jump density, µ = −h (0). In the long-range (t → ∞) the volatility reaches the stationary value

Results
We know that in terms of the characteristic function of the bidimensional process (x(t), r(t)) the discount function D(t) is given by (cf. Equation (7)) Then from Equation (29) we see that where D (0) (t) is the discount function for the continuous process in the absence of jumps and φ(t) ≡ φ(−i, t). The explicit form of these quantities is, respectively, given by (cf. Equations (27) and (24)) and Equation (33) constitutes the main result of this work and expresses the discount function of the jump-diffusion process in terms of the discount function D (0) (t) of the continuous OU process and the function φ(t) related to discontinuities.

Asymptotic Discount Function
We next analyze the asymptotic behavior as t → ∞ of the discount function (33). Firstly from Equation (34) we easily see that On the other hand, in order to get the asymptotic behavior of φ(t) we expand the jump characteristic functionh(θ) around the value θ = −i/α, and plugging it into Equation (35) we have Thus and in the long-time limit we have Let us note that if the sum on the right hand side of this expression is convergent and h(−i/α) is finite, we then have Substituting Equations (36) and (38) into Equation (33) yields where term discarded first order approximation (other remaining terms are exponentially small) reads The asymptotic discount function can thus be written as where the long-run discount rate defined in Equation (4) reads and is the long-run discount rate in the absence of jumps [26]. From Equation (42) we see that discontinuities will reduce the long-run discount rate as long as Let us remark that all expressions involving the long-run rate are meaningful as long ash(−i/α) exists. We will, however, see below some cases in whichh(−i/α) is infinite and r ∞ is meaningless.

Bounded and Symmetric Jump Density
We next develop condition (44) when the jump density h(u) is bounded and symmetric around u = 0. In other words, when sudden ups and downs of r(t) are finite and equally likely. From the definition ofh(θ), and bearing in mind the symmetry of h(u) around u = 0 (implying that h(u) = h(−u)), we may writeh that is,h Since cosh(u/α) > 1 and recalling that the normalization and symmetry of h(u) imply henceh(−i/α) > 1 and condition (44) holds. Therefore, for finite and symmetric jumps where ups and downs in return are equally likely, discontinuities always reduce the long-run discount rate. Let us recall that this conclusion remains valid as long as the integral in (47) is finite.
This result coincides with the Cox-Ingersoll-Ross process (rates are always positive) with negative and fixed jumps [13]. It is also possible to write an explicit expression of the long-run discount rate in terms of the jump PDF h(u) instead of the expression given by Equation (42), which gives r ∞ in terms of the characteristic functionh(−i/α). Indeed, taking into account Equation (47) and recalling the normalization of h(u), we write and substituting this result into Equation (42), we obtain clearly showing that r ∞ < r (0) ∞ for symmetrical jumps. See also the final part of Appendix B for an alternative approach. Appendix B is, however, fundamentally focused in providing the discount rate when the jump distribution is asymmetric, that is, when sudden ups and downs are not equally likely.

Some Specific Jump Distributions
We now study two particular examples of the jump distribution h(u). For these examples, we obtain the long-run discount rate and elucidate the meaning of condition (44) assuring that discontinuities reduce the long-run rate.

Fixed Jump Amplitudes
Let us first assume that the amplitudes of the discontinuities consist of a series of N fixed values, γ 1 , . . . , γ N (N = 1, 2, 3, . . . ). If these values are equally likely, the jump distribution function is In this case, the function φ(t) defined in (35) can be written as where Ei(·) is the Exponential integral defined as [42] Ei Expanding the integrand we get: and where The discount function D(t), given in Equation (33), now reads with the jump-free discount D (0) (t) given by Equation (34). Figure 1 shows the effect on the average discount rate d(t), defined in (3), of the presence of jumps with the simplest possible case when there is only one jump amplitude γ < 0 (N = 1). In this case, the average discount rate reads (cf. Equations (53) and (50)) In Figure 1, we assume a jump frequency λ = 0.02 1/year (1 jump every 50 years) and show how the discount d(t) changes as a function of time when considering no jumps and negative jumps of size |γ|/α = 0.25 and 0.5, where α is the reversion to the mean of the OU process (cf. Equation (16)). When considering the case of United States of America [17], it can be shown that small changes in γ < 0 parameter can lead to very sensitive effects to the discount rate, lowering the rate to values close to 1% and lower, even if jumps size are small or very small.  In the long-run, as t → ∞, we see from (51) and (52) that and the discount function can be approximated by Equation (41): where the long-run discount rate given in Equation (42) now reads and r ∞ is the jump-free discount rate, Equation (43). Let us also note that r ∞ < r The simplest case, when all jumps have the same amplitude, γ j = γ, reads (cf. Equation (55)) In this case, if there is a sudden decrease in return (γ < 0) then r ∞ < r ∞ . Otherwise, an increase of return (γ > 0) implies the increase of the long-run discount ratio, r ∞ > r (0) ∞ . These results are consistent with condition (44). Figure 2 shows how r ∞ change as a function of the scaled dimensionless jump length magnitude |γ|/α where α is the reversion to the mean of the OU process (cf. Equation (16)). Figure 2 shows the opposite effects in the long run discount function for positive and negative jumps when r (0) ∞ = 1.81%, which corresponds to the long-run discount rate estimated in Ref. [17] with United States of America (USA) real interest rate ratio datasets. For positive jumps γ > 0, the long-run rate in the case of USA can increase up to 3% with jump amplitudes of size 5% taking place once every 50 years. For negative jumps γ < 0, the long-run discount rates can become negative for jumps amplitudes of only 4%. These results confirm the high sensitivity of the long-run discount rate r ∞ when considering jumps.  Figure 2. The long-run discount rate r ∞ as a function of the the scaled jump size (|γ|/α) for jumps with fixed and bounded amplitudes (γ > 0 and γ < 0, cf. Equation (57)), with two fixed and bounded symmetric jumps ±γ (cf. Equation (58)) and for Laplacian jumps with average absolute value equals to γ (cf. Equation (68)). In all cases, jumps frequency λ is 1/50 years and r (0) Another particular case consists in assuming that discontinuities have only two possible amplitudes, which are equal but of opposite signs, that is, γ 1 = γ and γ 2 = −γ. Now (cf. (49))h (θ) = cosh(γθ/α) and the long-run discount rate (42) is Note that since cosh(γ/α) > 1 for all values of γ/α, then r ∞ < r (0) ∞ . In this example in which return suddenly decreases or increases equally likely by a fixed quantity, discontinuities always reduce the long-run rate, as we have already proved in a previous section. Figure 2 shows that this decrease can be quite sensitive even when jumps amplitudes are relatively small (r ∞ can already be negative when jumps sizes are of the order of 5% when considering USA datasets).

Laplacian Jump Amplitudes
As a second example, we suppose that jump amplitudes are not fixed but distributed according to the Laplace ("tent shape") density: (γ > 0). In this model, increasing and decreasing jumps are equally likely with zero average discontinuity, ∆r = 0, and σ = ∆r 2 = γ. Thus, the parameter γ represents the average of absolute values of the amplitude of discontinuities. Bearing in mind that Equation (59) represents a symmetric distribution around u = 0 we see that the Fourier transform of h(u) can be written as When θ ∈ R is real direct integration [42] yields Suppose, however, that θ = iξ (ξ ∈ R) is an imaginary number, in such a case since cos iθ = cosh ξ the integral in Equation (60) diverges when |ξ| ≥ √ 2/γ. We thus havẽ For the special case ξ = −1/α, we havẽ where is a dimensionless parameter that combines the average absolute jump γ and the strength α of the reversion to the mean of the OU process.
Recall that in terms of the jump-free discount, the discount function is given by (cf. Equation (33)) where D (0) (t) and φ(t) are, respectively, given by Equations (34) and (35). In the Appendix C, we show that for the Laplace density (59), the form and behavior of the discount function depends on the value of the parameter c defined in Equation (64). We have three cases: 1. When c < 1 (i.e., γ < α √ 2) we prove in the Appendix C that the function φ(t) is given by In this case, the discount function is finite and follows from Equation (65) after substituting Equation (66). Figure 3 illustrate this result considering the OU parameters estimated in Ref. [17] while considering different jumps frequencies λ and different jumps amplitudes in terms of c. For large values of t, we have (36) and (43)), we finally obtain and the expression for the long-run rate reads with r (0) ∞ given in Equation (43). In this model, discontinuities reduce the long-run discount rate if c < 1, which implies that r ∞ is finite. That is, when γ < α √ 2 and the average of the absolute value of discontinuities is smaller than the strength of the reversion to the mean represented by α √ 2. The behavior of the long-run discount rate r ∞ for Laplacian jumps can be compared to the fixed and bounded jumps amplitude case provided in Section 3.2.1. As shown in Figure 2, the behavior is qualitatively similar to the case of two symmetric jumps. The differences among both examples become relevant when c → 1 − , being the curve consistent with the critical behavior described in Equation (63).

2.
When c > 1 (i.e., γ > α √ 2), we prove in the Appendix C that the discount becomes infinite for times greater than a critical time, where while for t < t * the discount function is given as in case (i) above (even though now it has no sense asking for the asymptotic behavior of discount as t → ∞). Figure 4 shows the sensitivity of the critical time with respect to jumps amplitude. We there show that critical time t * can become shorter than a year or be strongly reduced as c increases.

3.
For the threshold case c = 1 (i.e., γ = α √ 2), the discount function grows exponentially. Thus, in the Appendix C we show that Note that this behavior is not contradictory with our previous results since, as r ∞ becomes negative, and discount turns into an increasing function for t large enough.   (66)) as a function of time (in years) and for different jumps time frequency λ. We take the Ornstein-Unlenbeck parameters estimated somewhere else with initial interest rate r 0 = 1% [17] with United States of America (USA) and Sweden (SWE) dates, which are considered to be stable countries. From this case, we conclude that if jump amplitudes are on average smaller than the restoring force toward the normal level (γ < α √ 2) jumps reduce the long-run discount rate. However, when jump amplitudes are larger (γ > α √ 2), the discount function becomes infinite at a finite time.

Discount in the Continuous Time Random Walk Formalism
Up to this point, we have dealt with a rate process described by a diffusion process in which there are superimposed finite discontinuities. At one extreme of the model, we find continuous diffusion processes with no discontinuities, which have been developed in our previous works [15,17,26]. At the other extreme, we find a purely discontinuous process where the return starting at some initial value keeps this value during a random interval of time and makes a sudden jump with a random amplitude to a new value, keeps the new value another random time interval, makes another jump, and so on. This is precisely the Continuous Time Random Walk (CTRW) with countless applications in many branches of natural sciences, engineering and economics and social sciences [43,44]. Let us observe that if waiting times between jumps are Poissonian with rate λ and independent of the jump amplitudes, then the PDF of the bidimensional process (x(t), r(t))-which we denote by p 0 (x, r, t|r 0 )-is described by the integro-differential Equation (14) with f (r) = g(r) = 0 with the initial condition By Fourier transforming this problem, we easily obtain the expression for the characteristic functionp(ω 1 , ω 2 , t|r 0 ). It reads (this solution can be also obtained from Equation (23) with the values α = 0 and k = 0): The marginal distribution of the return,p 0 (ω 2 , t|r 0 ), is obtained from (74) after setting and the return variance reads σ 2 0 (t) = λµt, where µ = −h (0) is the jump second moment. The discount function is obtained by setting ω 1 = −i and ω 2 = 0 in the joint characteristic function (74). For the CTRW model, this yields where Equation (76) is the expression of the discount function when rates are modeled by a Markovian CTRW with a general jump density h(u).
In the special case of Laplacian jumpsh(−iξ) (cf. Equation (62)) For times such that t < √ 2/γ, we see from Equation (77) In this case, the discount function reads On the other hand, if t > √ 2/γ, we can write However, due to Equation (78), the second integral is infinite, hence φ 0 (t) = −∞ and Therefore, for the CTRW model with Laplacian jumps, the discount function becomes infinite in the finite time t * = √ 2/γ.

Discussion
In a series of recent works [15][16][17]26,30], we have analyzed the process of discounting using mostly methods borrowed from non-equilibrium statistical physics and stochastic processes. In these works, we have considered three of the most popular stochastic models for the dynamics of interest rates: Ornstein-Uhlenbeck [27], Feller [28], and log-normal [29] processes, which are also very relevant in statistical physics. However, we are interested in real rates (that is, nominal rates corrected by inflation), which can be negative even during prolonged periods of time [15,17] and, since Feller and log-normal models deal exclusively with positive quantities, this leads to the Ornstein-Uhlenbeck (OU) process as the only model with mean reversion and allowing for negative rates. The Ornstein-Uhlenbeck (OU) model is the only one that allows for negative rates while still considering simplest (and linear) mean reversion towards a normal interest rate. Mean-reversion assumes that the interest rate follows a stationary process and this can be considered as a limitation when the model is contrasted with empirical data.
The work presented here continues with such an undertaking, but we go one step further and assume that, in addition to diffusive and continuous behavior, the sample paths of real rates r(t) also exhibit discontinuities. That is, we will model rates by a jumpdiffusion process as the economic evolution is known to occasionally have sudden bursts that hardly adjust to continuous diffusion-like evolution. We have thus wanted to elucidate the effect on the long-run discount of discontinuities that reflect the high-risk events that might occur in the future.
We have obtained a very general formula of the discount function for processes that combine an Ornstein-Uhlenbeck (OU) dynamics with the presence of Poissonian jumps with frequency λ (see also the Feller process with jumps in [13]). Equation (33) shows this key result. Two almost immediate questions are: how this general formula behaves for very long times? Additionally, which is the resulting long-run discount rate r ∞ ? Thus, ifh(θ) is the jump characteristic function, we have proved that as long ash(−i/α) exists, one obtains a value for r ∞ . Otherwise, the discount becomes infinite for a finite time horizon. An infinite discount is indeed catastrophic for the economy because it implies that any future value is zero. Furthermore, in case thath(−i/α) > 1, the addition of jumps to the model results into a decrease of r ∞ (cf. Equations (43) and (44)). This latter case entails a call for a more immediate action to climate change and it applies to those processes with symmetric jump amplitude. The obtention of lower long-run discount rates r ∞ for symmetric jumps amplitude is of particular importance, as it deepens in the idea already suggested by the OU process with no sudden jumps. That is to say, the fact that bounded unbiased uncertainty, no matter whether continuous or not, increases the urgency for immediate action.
To go deeper in evaluating the effect of discontinuities, we have gone through three different scenarios. The simplest case refers to the existence of a fixed jumps amplitude. We have been able to obtain the exact formula for the discount, which is represented in Figure 1 by considering the parameters of the OU process for United States of America [17] for a single negative jump γ. We there show that even by considering a frequency λ of one jump every 50 years and small jumps (γ = 0.5α) the effect is more than evident. As expected, sudden negative returns represent a quicker decrease in the discount function. Negative jumps of 5% return size (due to catastrophic news such as the COVID-19 outbreak) can already lead r ∞ to negative values, which in practice is telling us that immediate actions with strong investment are unavoidable to face, even if they are unknown, climate effects in the future (see Figure 2). In contrast, the assumption of having future sudden positive returns-due, for instance, to positive news or a major technological breakthroughincrements r ∞ , thus releasing pressure for taking action rapidly. It is, however, important to mention that the increment is lower than the decrement observed for negative jumps of the same size (see Figure 2).
A more sophisticated scenario is to consider the possibility of having two jumps of the same amplitude, but of opposite signs. This symmetric scenario also allows us to obtain the discount function, and this case always lowers the rate because of its symmetry. In the case of USA, as can also be shown in Figure 2, this drives r ∞ from 1.8% to 1% with jump amplitude γ = 5 % and frequency λ = (1/50) year −1 .
While keeping the same number of parameters (so it does not provide further complexity in data calibration), but considering a symmetric and continuous distribution of jumps amplitude (Laplace density), it is possible to provide a third scenario. This case interestingly contrasts to the case of two fixed and symmetric jumps. Laplacian jumps brings up two different discount solutions depending on the value of c = γ/α √ 2. The jump size average γ must be compared to the mean-reversion intensity α (cf. Equation (64)). If c < 1, it is possible to obtain an analytical formula of the discount function, which is carefully explored in Figure 3 with the OU estimated parameters with real interest rate datasets from USA and Sweden [17]. We there extend the analysis by exploring the effect of different values for jumps frequency λ and scaled jump amplitude c. However, when c > 1, the discount function becomes infinite for a critical time t * . Figure 4 shows that the critical time t * is quite sensitive to jumps amplitude size in a wide variety of cases and t * can become rather small, even for small jump amplitudes (For instance, few years, less than a decade, with jumps amplitude size about 7%, if we consider OU parameter from the case of USA).
Finally, it is also possible to obtain the discount function if we disregard the diffusion contribution. This scenario corresponds to the Continuous Time Random Walk and the results when assuming Laplacian jumps amplitudes are different depending on the critical time √ 2/γ. At short times, when t < √ 2/γ, the discount function decreases exponentially (cf. Equation (79)). However, as t > √ 2/γ, the discount function becomes infinite.

Conclusions
The main results of this paper are summarized in Table 1. Let us recall that the chief objective of the present work is to contribute to the mathematical modeling of discounting, within the context of environmental economics and climate action, by considering extreme situations or outliers. Although Martin Weitzman has already introduced the effect of fat tails in the economic evolution [34,35], our approach here is more dynamical and it is based not on fat tails, but on the addition of discontinuities in the economic time evolution. As we have shown such discontinuities may result, under certain conditions, into an infinite discount, confirming what Weitzman foresaw a few years ago. Let us emphasize that our approach (being complementary to that of Weitzman) provides a very general mathematical framework which allows for direct computation which can lead to future research in for instance option pricing by means of the martingale analysis and the characteristic function of the process. Table 1. Summary of the results proved in Section 3. The Ornstein-Uhlenbeck diffusion process already shows that the presence of noise (k = 0) reduces the long-run discount rate r (0) ∞ . The inclusion of Poissonian jumps with specific scenarios leads to several discount functions D(t) and several long-run discount rate r ∞ .

Main definitions
Discount function: Discount rate: ln D(t)/t Long-run discount rate: r ∞ = lim t→∞ ln D(t)/t Ornstein-Uhlenbeck (OU) dr/dt = −α(r − m)dt + kξ(t) r (0) ∞ = m − k 2 /(2α 2 ) OU and Poissonian jumps dr/dt = −α(r − m)dt + kξ(t) + n(t) Ifh(−i/α) is finite and jumps are symmetric Not defined Critical explosive time If jumps have one-fixed increasing amplitude We have also explored specific scenarios that show that discontinuities due to unexpected shocks (even if they represent downside shocks such as epidemics or a climate disaster or upside shocks due, for instance, to a new technological breakthrough) can severely affect current estimates of economic variables such as the long-run discount rate. Changes on this economic variables are shown to be strong enough to influence current decision mechanisms on whether and in which degree we shall take action today to face climate change.
We finally would like to stress the fact that that we here have wanted to include extreme and infrequent events to the rate process. The most natural choice to capture these events is to consider the addition of a Poisson jump process to an underlying diffusion process. However, there are other ways to incorporate a more generic class of jumps, which allows both frequent and infrequent jumps [45,46]. This generalization could lead to new and interesting theoretical results for the discount function that definitely would deserve careful attention and can drive new promising research avenues. Funding: This research was funded by MINEICO (Spain), Agencia Estatal de Investigación (AEI) and Fondo Europeo de Desarrollo Regional grant number PID2019-106811GB-C33; by Generalitat de Catalunya grant number 2017 SGR 608.

Conflicts of Interest:
The authors declare no conflict of interest.The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Observe that the two integrals are positive definite which implies that we will obtain a diminution of the long-run rate with respect to the jump free case, r ∞ < r (0) ∞ , if and only if the quantity within curly brackets is positive, that is Let us incidentally note that the numerator on the left-hand side of this inequality must be bounded, otherwise the expression of the long-run rate given in Equation (A17) is no longer valid.
For symmetrical jumps around the origin we have and Equation (A17) reduces to that is, which corresponds to the expression (48) of the main text.
Expanding this function in Taylor series around c = 1, we obtain which substituting into Equation (A26) and after taking the limit c → 1 yields and as t → ∞ we have and discount will be eventually dominated by this term, leading to which shows that the discount function increases in an explosive way.