--- title: Parisian Ruin Probability url: https://www.emergentmind.com/topics/parisian-ruin-probability type: topic --- # Parisian Ruin Probability Parisian ruin probability is the probability that a surplus or reserve process exhibits a continuous excursion below a prescribed level for at least a specified delay, rather than merely crossing the level instantaneously. In the literature, the delay may be deterministic, exponentially distributed, or replaced by a cumulative occupation-time requirement, and the underlying risk process may be Gaussian, Brownian, spectrally negative Lévy, discrete-time, or multidimensional. The subject therefore sits at the intersection of fluctuation theory, Gaussian extremes, occupation-time asymptotics, and stochastic control [1610.00522] [1102.4055] [2104.12638] [1811.10110]. ## 1. Formal definition and principal variants A standard deterministic-delay formulation declares ruin when there exists a time interval of fixed length during which the reserve remains below zero. In one dimension, a representative event is \[ \mathbb{P}\Big\{\exists\, t\ge 0:\forall s\in[t,t+T]\ \ R(s)<0\Big\}, \] while on a finite horizon one often studies \[ \mathcal{P}_S(u,T_u)=\mathbb{P}\Big\{\inf_{t\in[0,S]}\sup_{s\in[t,t+T_u]}R_u(s)<0\Big\}. \] These formulations make the contrast with classical ruin explicit: classical ruin requires a single hit below zero, whereas Parisian ruin requires a sustained subzero excursion [2103.03213] [1610.00522]. A second major variant replaces the deterministic delay by an exponential clock. In the lifetime exponential Parisian formulation, if \(g_t:=\sup\{s\in[0,t]:W_s\ge 0\}\) is the last time wealth was nonnegative, then ruin occurs when the current negative excursion length exceeds an independent exponential clock; one form is \[ \kappa^\pi:=\inf\{t>0:t-g_t>e_\rho^{g_t}\}, \] with the clock reset at each return above zero. A closely related discounted formulation uses \[ \tau_P=\inf\{t>0:t-g_t>T_p\}, \] where \(T_p\) is exponential with hazard rate \(p\) [2104.12638] [2001.08344]. A third variant is cumulative Parisian ruin. Here the relevant object is not a single uninterrupted excursion but the total time spent in the ruin set. In the multi-line Brownian setting, \[ \tau_r(u)=\inf\Big\{t>0:\int_0^t \mathbf{1}\{U(s)<0\}\,ds>r\Big\}, \] with the inequality understood componentwise. In two-dimensional Brownian models over \([0,1]\), the cumulative condition is often scaled as \[ \int_0^1 \mathbf{1}\{R_1(t)<0,\ R_2(t)<0\}\,dt>\frac{L}{u^2}. \] This replaces excursion length by occupation time as the operative persistence variable [1811.10110] [2001.09302] [2109.12906]. The literature also includes a drawdown-based formulation for spectrally negative Lévy insurance risk processes. If \(Y_t=\bar X_t\vee y-X_t\) is the drawdown from the last record maximum, then Parisian drawdown ruin occurs when \(Y\) stays above a fixed drawdown level \(a\) for at least \(r\) consecutive units of time. This is distinct from fixed-level ruin because the barrier is measured relative to the running maximum rather than an absolute solvency threshold [1806.02083]. ## 2. Gaussian and Brownian asymptotics For finite-horizon Gaussian risk models, the principal asymptotic structure is governed by the local variance maximum and local correlation decay near the most likely ruin point. In the model \[ R_u(t)=u+ct-X(t), \] with \(X\) centered Gaussian, one assumes that the standard deviation \(\sigma(t)\) attains its unique maximum at the horizon endpoint \(S\), with local variance decay exponents \(\beta_1,\beta_2\), and local correlation roughness exponent \(\alpha\). The exact asymptotics of \[ \mathcal{P}_S(u,T_u)=\mathbb{P}\Big\{\inf_{t\in[0,S]}\sup_{s\in[t,t+T_u]}R_u(s)<0\Big\} \] then split into three regimes. If \(\alpha<\beta_1\) and \(T_u u^{2/\alpha}\to T\), generalized Pickands constants appear together with a polynomial factor in \(u\); if \(\alpha=\beta_1\), generalized Piterbarg-type constants arise; if \(\alpha>\beta_1\) and \(T_u\) vanishes sufficiently fast, then \[ \mathcal{P}_S(u,T_u)=\Psi\Big(\frac{u+cS}{\sigma(S)}\Big)(1+o(1)). \] Thus the Parisian delay may alter the leading constant and even the power prefactor, while leaving the Gaussian tail form intact [1504.07061]. The integrated Gaussian risk model introduces discounting directly into the reserve: \[ R_u(t)=u+c\int_0^t e^{-\delta(s)}\,ds-\int_0^t e^{-\delta(s)}Z(s)\,ds, \] with \(Z\) centered Gaussian and nonnegative covariance. In this setting, the variance function \[ \sigma^2(t)=2\int_0^t\int_0^v e^{-\delta(w)-\delta(v)}\operatorname{Cov}(Z(w),Z(v))\,dw\,dv \] is strictly increasing. For any bounded \(T_u\ge 0\), the Parisian and classical ruin probabilities are asymptotically identical on the log-scale: \[ \lim_{u\to\infty}\frac{\log \mathcal{P}_S(u,T_u)}{u^2} = \lim_{u\to\infty}\frac{\log \psi_S(u)}{u^2} = -\frac{1}{2\sigma^2(S)}. \] If \(T_u\to 0\), then the equivalence is sharper: \[ \mathcal{P}_S(u,T_u)\sim \psi_S(u)\sim \Psi\!\left(\frac{u+c\widetilde{\delta}(S)}{\sigma(S)}\right), \] and the conditional ruin time satisfies \[ u^2(S+T_u-\tau(u))\xrightarrow{d}\mathrm{Exp}\!\left(\frac{\sigma'(S)}{\sigma^3(S)}\right). \] This identifies a boundary-localized exponential limit law for the Parisian ruin time [1610.00522]. Self-similar Gaussian risk processes \[ R_u(t)=u+c\,t^\beta-X_H(t) \] with \(\operatorname{Var}(X_H(t))=t^{2H}\) exhibit a related but infinite-horizon picture. Under local stationarity of \(X_H(t)/t^H\) near the unique maximizer \[ t_0=\Big(\frac{H}{c(\beta-H)}\Big)^{1/\beta}, \] the Parisian ruin probability satisfies \[ \mathbb{P}\{\tau_u^*<\infty\} = \frac{\mathcal{F}_\alpha(D_0T)}{\mathcal{H}_\alpha}\, \mathbb{P}\{\tau(u)<\infty\}\,(1+o(1)), \] where \(T\) is the asymptotic delay scale, \(\mathcal{H}_\alpha\) is the Pickands constant, and \(\mathcal{F}_\alpha\) is its Parisian analogue. The Parisian and classical ruin times have the same Gaussian limit after centering and scaling, and their difference is negligible on that scale [1405.2958]. Brownian models with constant force of interest provide explicit finite- and infinite-horizon analogues. For \[ R_u^\delta(t)=e^{\delta t}\Big(u+c\int_0^t e^{-\delta s}\,ds-\sigma\int_0^t e^{-\delta s}\,dB(s)\Big), \] finite-horizon asymptotics take the form \[ \mathcal{K}_S^\delta(u,T_u)\sim \mathcal{P}(aT)\,\Psi(z_u), \] with \(u^2T_u\to T\), whereas the infinite-horizon model yields \[ \mathcal{K}^\delta(u,T_u)\sim \mathcal{P}_a^f[0,\infty)\, \Psi\!\left(\frac{1}{\sigma}\sqrt{2\delta u^2+4cu}\right). \] In both cases the Parisian factor is a generalized Parisian/Piterbarg constant, and the conditional ruin time is asymptotically exponential after \(u^2\)-scaling [1606.07339] [1702.06091]. Recent work extends this framework to locally self-similar Gaussian drivers with power-type trend. The asymptotic form remains \(c\,u^p\Psi(u)\), but the constants and exponents depend on the interaction between the local self-similarity index, the variance-drop exponent, the trend exponent, and the local roughness parameter. Parisian Pickands-type constants continue to govern the leading term, now for families of limiting Gaussian fields rather than a single limiting process [2604.00916]. ## 3. Spectrally negative Lévy and drawdown formulations For spectrally negative Lévy processes, Parisian ruin admits compact fluctuation-theoretic expressions in terms of scale functions. If \(X\) is a spectrally negative Lévy process with \(E[X_1]>0\), the deterministic-delay Parisian ruin time is \[ \kappa_r=\inf\{t>r:t-g_t>r\}, \qquad g_t=\sup\{0\le s\le t:X_s\ge 0\}, \] and the ruin probability satisfies \[ P_x(\kappa_r<\infty) = 1-E[X_1]\, \frac{\int_0^\infty W(x+z)\,z\,P(X_r\in dz)} {\int_0^\infty z\,P(X_r\in dz)}. \] This formula is valid for all \(x\in\mathbb{R}\), involves only the scale function \(W\) and the law of \(X_r\), and yields the correct limiting cases: \[ P_x(\kappa_r<\infty)\to P_x(\tau_0^-<\infty) \quad\text{as }r\downarrow 0, \qquad P_x(\kappa_r<\infty)\to 0 \quad\text{as }r\uparrow\infty. \] The monotonicity in both \(x\) and \(r\) is immediate from the same representation [1102.4055]. An earlier deterministic-delay analysis for spectrally negative Lévy insurance risk processes expresses Parisian ruin through classical ruin, the undershoot distribution, and the probability that a negative excursion survives longer than the delay. Its central decomposition is \[ P_x(\tau^\zeta<\infty) = P_x(\tau_0^-<\infty)\,P(\tau^\zeta<\infty) + \int_{(0,\infty)} P_{-z}(T_0^+>\zeta)\, P_x(\tau_0^-<\infty,-X_{\tau_0^-}\in dz). \] Within this framework the paper derives both Cramér-type and convolution-equivalent asymptotics. In the Cramér regime, \[ P_x(\tau^\zeta<\infty)\sim C(\zeta)e^{-Rx}, \] while in the convolution-equivalent regime the asymptotic scale is the integrated tail of the Lévy measure. This shows that the deterministic Parisian delay modifies constants but preserves the principal exponential or heavy-tail scale determined by the underlying Lévy fluctuations [1003.4299]. The drawdown formulation changes the geometry of ruin. For the drawdown process \[ Y_t=\bar X_t\vee y-X_t, \] Parisian drawdown ruin above level \(a\) over delay \(r\) is defined by the first excursion of \(Y\) above \(a\) that lasts at least \(r\). The resulting transform identities are semi-explicit in terms of the scale function and the one-time law of the Lévy process, and the joint Laplace transform of ruin time and position at ruin can be written in closed transform form. The most striking feature is qualitative rather than algebraic: \[ P_{y,x}(\tau_r^a<\infty)=1. \] This contrasts sharply with fixed-level Parisian ruin, where ruin probability can be strictly less than one under positive drift [1806.02083]. ## 4. Discrete-time and grid-based ruin In discrete time, Parisian ruin requires the surplus to remain nonpositive for a prescribed number of consecutive periods. For the classical discrete-time reserve \[ R_n=u+n-S_n, \qquad S_n=\sum_{i=1}^n Y_i, \] with \(\mu=E[Y_1]<1\), the Parisian ruin time with delay \(d\in\{1,2,\dots\}\) is \[ \tau^d = \inf\left\{ n\in\mathbb{N}: n-\sup\{s0\}>d,\ R_n\le 0 \right\}. \] Finite-time survival admits an exact decomposition through the first classical ruin epoch, the deficit at ruin, and Kendall’s identity for subsequent recovery. Infinite-time ruin is then expressed in terms of classical ruin probabilities and the probability of returning to \(+1\) within \(d\) steps after entering the nonpositive region. On the asymptotic side, light-tailed claims yield \[ \mathbb{P}_u(\tau^d<\infty)\sim C\Big[1-(1-\mathbb{P}_1(\tau^d<\infty))\,f(d)\Big]e^{-\gamma u}, \] so the Parisian delay changes the prefactor but not the Cramér exponent. In heavy-tailed regimes, the effect depends on the tail class: for \(F_I\in\mathcal{S}^{(0)}\), the leading asymptotic is unchanged from classical ruin, whereas for \(F\in\mathcal{S}^{(\alpha)}\), \(\alpha>0\), the delay modifies the leading constant [1403.7761]. The discrete-time dual risk model reverses the drift structure: \[ R_n^*=u-n+X_n, \qquad X_n=\sum_{i=1}^n Y_i, \] with \(E[Y_1]>1\). Here classical ruin is hitting \(0\), and Parisian ruin requires the process to remain strictly negative for \(r\) consecutive periods. The finite-time Parisian ruin probability is obtained recursively from classical dual ruin probabilities. In infinite time, the model has an especially simple structure: if \(A\in[0,1)\) is the unique solution of \[ \widetilde p(A)=A, \] where \(\widetilde p\) is the claim-size pgf, then the classical dual ruin probability is \(\psi^*(u)=A^u\), and the Parisian probability factorizes as \[ \psi_r^*(u)=D A^u \] for a constant \(D\) depending on the delay and the gain distribution but not on \(u\). The paper works out this factorization explicitly for the Binomial/Geometric model and for the Gambler’s ruin problem [1708.06785]. A third discrete framework retains continuous-time Brownian dynamics but allows ruin observation only on a uniform grid \(G(\delta)=\{0,\delta,2\delta,\dots\}\). In this model, \[ R_u(t)=u+ct-B(t), \] and grid-based Parisian ruin with delay \(T\) has the asymptotic form \[ \mathcal{P}_\delta(u,T)\sim \mathcal{H}_{2c^2\delta,\,2c^2T}\,e^{-2cu}. \] Thus the grid affects only the multiplicative constant, through a discrete Parisian Pickands-type constant, while the leading exponential scale remains \(e^{-2cu}\). The associated ruin time, conditionally on ruin, has the same Gaussian limit as in the continuous model after centering by \(u/c\) and scaling by \(\sqrt{u}/c^{3/2}\) [2001.10311]. ## 5. Multidimensional Brownian and reinsurance geometries Multidimensional Parisian ruin introduces both persistence and geometry. In the quota-share insurer–reinsurer model, both companies share a single Brownian or fractional Brownian loss stream, so simultaneous Parisian ruin reduces to crossing a one-dimensional piecewise-linear barrier: \[ \mathbb{P}\Big\{\exists t\ge 0:\forall s\in[t,t+T_u]\ B_H(s)>q_1u+c_1s,\ B_H(s)>q_2u+c_2s\Big\}. \] The asymptotic behavior depends critically on the intersection time \[ t_*=\frac{q_2-q_1}{c_1-c_2} \] relative to the variance-maximizing points \[ t_1=\frac{Hq_1}{(1-H)c_1}, \qquad t_2=\frac{Hq_2}{(1-H)c_2}, \] and on the scaling \[ T_u\,u^{1/H-2}\to T. \] In the outer regimes \(t_*\notin(t_1,t_2)\), the asymptotics are controlled by one side of the barrier; in the intersection regime \(t_*\in(t_1,t_2)\), two-sided Piterbarg-type constants enter. For \(H<1/2\), the growth condition on \(T_u\) is necessary; keeping \(T_u\equiv T>0\) produces substantially smaller probabilities [2103.03213]. In two-dimensional Brownian risk models with correlated marginals, simultaneous Parisian and cumulative Parisian ruin over a finite horizon typically live on the \(u^{-2}\) time scale. For \[ R_1(t)=u+c_1t-W_1(t),\qquad R_2(t)=au+c_2t-W_2(t), \] with correlation \(\rho\), simultaneous Parisian ruin with delay \(S/u^2\) satisfies \[ v_S(u,au)\sim C_{a,\rho}(S)\, \mathbb{P}\{W_1(1)>u+c_1,\ W_2(1)>au+c_2\}, \] and cumulative Parisian ruin with threshold \(L/u^2\) has \[ Y_L(u,au)\sim K_{a,\rho}(L)\, \mathbb{P}\{W_1(1)>u+c_1,\ W_2(1)>au+c_2\}. \] The conditional cumulative ruin time then has an exponential limit, with rate depending on \((a,\rho)\) in the genuinely two-dimensional regime and equal to \(1/2\) in the effectively one-dimensional regime [2001.09302]. The non-simultaneous Brownian Parisian model replaces a common time window by separate windows for each component. The asymptotics are then conditional on non-simultaneous classical ruin and exhibit a phase transition at \[ A_a=\frac{1-\sqrt{8a^2+1}}{4a}. \] When \(\rho>A_a\), the limiting conditional Parisian probability is a ratio of two-dimensional Parisian constants \(\mathcal{R}_{S_1,S_2}/\mathcal{R}_{0,0}\); when \(\rho=A_a\) or \(\rho0\), \(\pi^*(w)\) coincides with the classical lifetime-ruin minimizer and is independent of \(\rho\). For \(w<0\), \(\pi^*(w)\) is strictly larger than the classical rule, increases with the excursion hazard \(\rho\), and also increases with the mortality rate \(\lambda\). For small \(\rho\), \[ \psi(w;\rho)=\rho\,m(w)+O(\rho^2), \] so the minimal probability of lifetime exponential Parisian ruin is asymptotically proportional to the minimum expected occupation time below zero [2104.12638]. A reinsurance analogue is developed for the discounted probability of exponential Parisian ruin. In the classical risk model controlled by per-loss reinsurance \(R_t(y)\), with mean-variance reinsurance premium principle, the value function \[ V(x)=\inf_{R\in\mathcal R}E_x\big[e^{-\beta\tau_P}\mathbf 1_{\{\tau_P<\infty\}}\big] \] solves an integro-differential HJB equation of the form \[ 0= \beta V(x)+p(V(x)-1)\mathbf 1_{x<0} +\inf_{R\in\mathcal R} \Big\{[-\kappa+\lambda G(R)]V'(x)+\lambda E[V(x-R(Y))-V(x)]\Big\}, \] with boundary conditions \[ \lim_{x\to-\infty}V(x)=\frac{p}{p+\beta}, \qquad \lim_{x\to+\infty}V(x)=0. \] The analysis uses stochastic Perron’s method and proves that \(V\) is the unique continuous viscosity solution despite the discontinuity of the Hamiltonian at \(x=0\). The optimal control has a sharp structural split: full retention is optimal for \(x<0\), while for \(x\ge 0\) the optimal retention is characterized implicitly by the first-order condition \[ (1+\theta)+\eta y-\eta R(y;\gamma_1)-e^{\gamma_1 R(y;\gamma_1)}=0 \] whenever an interior solution exists. This yields excess-of-loss behavior under the expected-value principle and proportional behavior under the variance principle as special cases [2001.08344]. Taken together, these control results show that Parisian ruin has become a genuine optimization criterion. A plausible implication is that the delay mechanism is not merely a technical modification of classical ruin, but a state-dependent objective that can reverse comparative statics when the process is already below zero: in both investment and reinsurance formulations, the negative region induces qualitatively different optimal behavior than the solvent region [2104.12638] [2001.08344]. Source: https://www.emergentmind.com/topics/parisian-ruin-probability