key: cord-0979639-20pzleok
authors: Li, Yaohan; Dong, You; Qian, Jing
title: Higher-Order Analysis of Probabilistic Long-Term Loss under Nonstationary Hazards
date: 2020-06-24
journal: Reliability engineering & system safety
DOI: 10.1016/j.ress.2020.107092
sha: 18e579d09abddbf1a9dff15f267d6cabcb669d97
doc_id: 979639
cord_uid: 20pzleok

Abstract Quantification of hazard-induced losses plays a significant role in risk assessment and management of civil infrastructure subjected to hazards in a life-cycle context. A rational approach to assess long-term loss is of vital importance. The loss assessment associated with stationary hazard models and low-order moments (i.e., expectation and variance) has been widely investigated in previous studies. This paper proposes a novel framework for the higher-order analysis of long-term loss under both stationary and nonstationary hazards. An analytical approach based on the moment generating function is developed to assess the first four statistical moments of long-term loss under different stochastic models (e.g., homogeneous Poisson process, non-homogeneous Poisson process, renewal process). Based on the law of total expectation, the developed approach expands the application scope of the moment generating function to nonstationary models and higher-order moments (i.e., skewness and kurtosis). Furthermore, by employing the convolution technique, the proposed approach effectively addresses the difficulty of assessing higher-order moments in a renewal process. Besides the loss analysis, the mixed Poisson process, a relatively new stochastic model, is introduced to consider uncertainty springing from the stochastic occurrence rate. Two illustrative examples are presented to demonstrate practical implementations of the developed approach. Ultimately, the proposed framework could aid the decision-maker to select the optimal option by incorporating higher-order moments of long-term loss within the decision-making process.

1 Highlights:

 An efficient and novel approach to derive higher-order statistical moments of the discounted long-term loss  Higher-order analysis to aid decision-making and optimal structural design of civil infrastructure  Large uncertainty associated with the stochastic occurrence and magnitude of hazards  An integrated framework for long-term loss assessment under stationary and nonstationary hazards

In recent decades, the devastating effects that hazards have on societies worldwide have intermittently raised the attention of governments and the public to hazard risk assessment and management. For civil infrastructure, various hazards (e.g., earthquakes, hurricanes, and progressive deterioration) may impair structural functionality, thus resulting in severe consequences. The hazard-induced consequences are commonly measured in terms of financial losses (e.g., repair cost), social losses (e.g., downtime, deaths), and environmental losses (e.g., carbon dioxide emissions). Due to various sources of uncertainty, the accumulated losses can be aggravated over the service life of civil infrastructure. Hence, the quantification of long-term loss of civil infrastructure subjected to hazards is of significant importance to aid the decision-maker to mitigate potential losses and enhance preparedness [1, 2] .

Uncertainty quantification plays a significant role in long-term loss assessment. There are large uncertainties associated with the frequency and magnitude of hazards [3, 4] . In previous studies, stationary models (e.g., homogeneous Poisson process) have been widely used to characterize the probabilistic behavior of hazards. Nowadays, studies show strong evidence that hazard arrivals may follow nonstationary behavior [5] . For instance, the renewal process is proposed to model the earthquake arrivals, in order to incorporate the time-varying energy accumulation of the fault [6, 7] . Such time-dependent trends are also identified in other hazards, such as more frequent hurricane landfalls under climate change [8] , increased wind speeds [9] , extreme precipitation [10] , and sea level rise [11] . In addition to natural hazards, the progressive deterioration of structural systems is also stochastic and time-variant [12] . Given the time-dependent trend of hazards, a general framework is needed to evaluate the long-term loss of civil infrastructure under both stationary and nonstationary models. This concern is discussed and addressed in this paper. 4 The analytical formulation of long-term loss is essential for risk assessment and management. Although numerical modeling is accessible, simulations are usually computationally expensive and time-consuming. Low-order moments (i.e., expectation and variance) of long-term loss based on stationary stochastic models have been investigated by previous studies [13, 14] . Based on the homogeneous Poisson process, an analytical formulation of the expected life-cycle cost of buildings under single and multiple hazards was presented by Wen and Kang [15] . Recently, several studies assessed the long-term loss (e.g., mean and variance) of civil infrastructure under nonstationary processes. For instance, Yeo and Cornell [16] proposed analytical expressions for the expected loss caused by earthquakes using homogeneous and non-homogeneous Poisson models. Wang et al. [17] computed the mean and variance of hurricane-induced damage loss using the non-homogeneous Poisson process. Lin and Shullman [18] assessed the risk of New York City being damaged by hurricanes and surge flooding in a nonstationary environment. To simplify the computational process, the nonstationary Poisson model was converted into a stationary one in these studies.

Additionally, the long-term loss is limited to Poisson models and the first two moments.

In addition to Poisson processes, recent studies proposed some new nonstationary models for the loss assessment. For instance, Pandey and Van Der Weide [19] used a stochastic renewal process with Brownian Passage Time distribution to formulate the expectation and variance of the discounted damage cost of a structure under earthquakes. The derivations were based on the renewal decomposition properties of renewal processes. The renewal model was also used to evaluate the lifetime resilience and cost of structural systems considering progressive deterioration [12] . Although the renewal approach provides an alternative option to assess the loss under nonstationary hazards, it cannot be applied to other stochastic models, such as the non-homogeneous Poisson process. Meanwhile, higher-order moments are not taken into account in previous studies.

Though the minimum expected loss has been widely used as a standard decision criterion, it is only suitable for risk-neutral decision makers. This criterion cannot cope with different attitudes [20]. Goda and Hong [21] indicated the structural design based on the expected life-cycle cost may not be optimal, and stated the need for statistical moments (e.g., variance, skewness, and kurtosis) of the cost. Furthermore, the mean-variance criterion is sufficient only when the utility function within the decision-making process is quadratic or the investment return (e.g., loss represents a negative return) follows a normal distribution [22] .

Higher-order moments, i.e., skewness and kurtosis, measure asymmetry and tail conditions of the distribution with respect to the long-term loss. In risk management, large skewness and kurtosis of loss imply heavy tail risks. Such undesired risks are associated with low-probability events with disastrous consequences, e.g., credit risk crisis [23] and COVID-19 pandemic [24] . Higher-order moments are required when risk preferences of decision makers are considered, e.g., in the stochastic dominance criteria [20, 25] . For instance, a decision maker with the absolute risk-averse attitude prefers positive skewness and small kurtosis of the investment return, as highly skewed data with large kurtosis indicate an increased likelihood of extreme losses [26] . Different decision results may be obtained due to the exclusion of these moments. Therefore, the assessment of higher-order moments of longterm loss is necessary. These moments can be used to aid the decision-making and optimal structural design of civil infrastructure by considering different attitudes.

To the best knowledge of the authors, mathematical expressions of higher-order moments of long-term loss have not been developed for Poisson and renewal models. In this paper, a novel approach based on the moment generating function is developed to formulate the higher-order moments of long-term loss of civil infrastructure subjected to hazards. In this context, the long-term loss refers to hazard-induced financial losses in a long time interval. In addition to the loss analysis, a new stochastic model of the mixed Poisson process is presented to take the uncertainty springing from stochastic occurrence rate into consideration.

This new model has been recently introduced in hurricane landfall simulations [27] and rainfall occurrence models [28] to consider environmental variability.

Overall, this paper develops an integrated framework to assess the probabilistic longterm loss under stationary and nonstationary hazards. Higher-order moments of the loss are formulated based on moment generating functions for the homogeneous Poisson, nonhomogeneous Poisson, mixed Poisson, and renewal processes. The merits, application scope, and limitations of the proposed method are explained and discussed. Two illustrative examples are presented to demonstrate implementations of the approach. The remainder of the paper is organized as follows. The subsequent two sections introduce the model setup for long-term loss and the moment generating function-based approach. Section 4 presents two case studies to interpret the applications. The last section summarizes major outcomes and conclusions. 6

Hazard arrivals are commonly modeled as a stochastic process using historical observations. For instance, the homogeneous Poisson process is widely used to model the earthquake recurrence [29, 30] . The occurrence of hurricanes could also be represented by the homogeneous Poisson process based on historical records [31, 32] . In addition to using historical data, the future projection of hazard arrivals requires the consideration of variability in hazard characteristics [4, 33] . Various non-stationary models are developed to capture such variability. For instance, the renewal process is used to incorporate the time-dependent characteristics of ground motion activities for long-term seismic forecasting [34] . The nonhomogeneous and mixed Poisson processes are utilized to model the occurrence of hurricanes in a changing climate [35, 36] . In this section, theoretical descriptions and notations of four typical stochastic models (i.e., homogeneous Poisson process, renewal process, nonhomogeneous Poisson process, and mixed Poisson process) are presented. Based on stochastic models, the long-term loss is formulated.

As one of the most classical stationary models, the homogeneous Poisson process is adopted for stationary hazard arrivals and illustrated by using two different models: the homogenous Poisson and renewal processes.

A homogeneous Poisson process is an arrival process with stationary increments. It has a constant occurrence rate λ. The occurrence rate is defined as the number of hazard events within the time unit (e.g., per year), which is also known as the intensity function. A homogeneous Poisson process can be denoted by the number of occurrence {N(t int ), t int > 0} Poisson process. In this context, a process is considered as 'nonstationary' when the occurrence rate is not constant.

The occurrence of hazards could be time-dependent. For instance, after an earthquake, there can be a long period before the next earthquake, in which the accumulated elastic strain energy is released by the fault. By considering such time-dependent characteristics, a non-Poisson process can be used to model the occurrence of earthquakes over a long period of time [37, 34] . In the nonstationary renewal processes, different probabilistic models of interarrival times are used to quantify the time-dependent characteristics. For instance, several distributions of the inter-arrival time are provided in earthquake engineering, including lognormal [38] , gamma [39] , and Brownian Passage Time distribution [34] .

Another nonstationary model used for the long-term loss assessment is the non-homogeneous Poisson process, with a time-dependent occurrence rate. A renewal process is not necessarily 8 limited to the exponentially distributed inter-arrival times, while the Poisson process is not subjected to the time-independent occurrence rate. For instance, the non-homogeneous

Poisson process can be used to model hurricane landfalls in a changing climate [18, 35] . The frequency of hurricanes may increase due to the impact of the warming climate. For this process, the occurrence rate is time-varying, e.g., λ(t). The expected number of hazard arrivals over the time interval (0, t int ] can be computed as

The probability of having n number of arrivals within the time interval (0, t int ] is computed as

In the homogeneous and non-homogenous Poisson processes, the occurrence rate is restricted to a deterministic intensity function (i.e., either a constant λ or the time-dependent λ(t)).

However, the deterministic function may not be sufficient to capture the uncertainty in a long-term trend. The mixed Poisson process, which covers the uncertainty within intensity function, is receiving increased attention in the hazard model and should also be evaluated.

For instance, Xiao et al. [27] suggested using a stochastic intensity function to model hurricane occurrence when considering seasonal variability, based on the recorded hurricane activities along the U.S. Gulf and Atlantic coasts between 1900 and 2010.

The mixed Poisson process is known as a special case of the cox process. A cox process is a Poisson process, in which the occurrence rate is a stochastic process (denoted as {G(t int ), t int > 0}). A cox process becomes a mixed Poisson process when the rate is a random variable, e.g., {G(t int )} = Λ. Given the rate Λ, increments in the mixed Poisson process stay stationary [40] . However, these increments are no longer independent, as the number of arrivals relies on the distribution of Λ. Herein, the stochastic rate Λ is assigned to follow a gamma distribution with shape parameter α and rate parameter β. The rate Λ is continuous with a probability density function g(x). Hence, the probability density function of the gamma distributed rate Λ ~ Γ(α, β) can be assessed as

Consequently, the probability of having n number of arrivals within the time interval (0,

Based on the stochastic models, the long-term loss of civil infrastructure subjected to hazards can be formulated. The time interval (0, t int ] is defined as the service life of civil infrastructure. During the time interval, hazard arrivals are modeled as a stochastic process and the total number of hazard events is N(t int ) = n with N(0) = 0. The arriving time and interarrival time of the kth event can be defined as T k and W k , respectively. The arriving time T k equals the sum of inter-arrival times, i.e., T k = W 1 + W 2 + … + W k , as described previously. L is the loss severity. L k refers to the financial loss due to structural damage/failure under the kth hazard event. The loss severity L k and the inter-arrival time W k are assumed to be independent. The long-term loss, denoted as LTL(t int ), is the hazard-induced financial loss of civil infrastructure subjected to hazards within the investigated period t int . The loss is discounted to the present value using a monetary discount rate r, as indicated in Figure 1 . The

The long-term loss LTL(t int ) is related to the stochastic models of the hazard. For instance, if a homogeneous Poisson process is implemented for loss assessment, the inter-arriving time W k follows an exponential distribution and the number of arrivals N(t int ) has a Poisson distribution. These parameters are different under different models. 

This section introduces the theoretical fundamentals of the proposed moment generating function-based approach. Analytical derivations of moment generating functions and statistical moments of long-term loss under the homogeneous Poisson process, nonhomogeneous Poisson process, mixed Poisson process, and renewal process are presented.

The merit, capability, and limitations of the developed approach are also discussed.

The moment generating function uniquely dictates the probability distribution of a random variable. This property can be adopted to formulate probability distributions of random variables. For instance, random variables will have the same probability distributions if they have identical moment generating functions. Another property of the moment generating function is that raw moments (e.g., mean, second moment) can be obtained by taking derivatives [40, 41] . This property is utilized to formulate statistical moments of long-term loss in this paper.

For a random variable X, its moment generating function about η (η  ℝ) is defined as

The first two raw moments of X can be obtained by taking the first and second derivatives of the moment generating function at zero

Similarly, the mth-order moment can be assessed by taking the mth derivative at zero

Based on this concept, the key to derive higher-order moments of long-term loss LTL(t int )

is to compute its moment generating function int () LTL t  . (12) in which Φ L refers to the moment generating function of loss severity L. The derivation of Φ L requires the information with respect to the probabilistic distribution of loss severity.

Previous studies indicated that the probabilistic loss severity can be modeled by the exponential distribution [42, 43] . Herein, the loss severity L is assumed to follow an exponential distribution L ~ EXP(θ), with the mean E[L] = 1/θ. By substituting Φ L into Eq. (13) By taking the first and second derivatives at zero, the expectation and variance of longterm loss can be obtained

Likewise, the mth order moment can be assessed using Eq. (11) . The expressions of skewness and kurtosis are associated with the third and fourth-order raw moments. For instance, the skewness u 3 and kurtosis u 4 are (17) in which μ refers to the expected long-term loss E[LTL(t int )] and σ is the standard deviation of LTL(t int ). Consequently, skewness and kurtosis of long-term loss can be obtained int int 3 1/2

From Eqs. (18) and (19), it shows the skewness and kurtosis of long-term loss are not affected by θ. Hence, the skewness and kurtosis under the homogeneous Poisson process are independent of the exponentially distributed loss severity. When other distributions are used for the loss severity, the higher-order moments may be affected.

The same technique can be used for the non-homogeneous model by applying the law of total expectation. For the non-homogeneous Poisson process, the expected number of arrivals

Given the number of arrivals, the moment generating function of long-term loss can be derived as

The detailed derivations of Eqs. (12) and (21) are shown in Appendix A. If a linear function is used for the increasing occurrence rate, e.g., λ(t int ) = λ 0 (1 + ct int ), the moment generating function becomes (22) in which λ 0 is the initial stationary occurrence rate and c refers to an annual increase rate of hazard occurrence. For the given linear rate function, the moment generating function may not be differentiable at zero. Under this circumstance, the moments can be computed by taking derivatives and finding limits by approaching zero. Accordingly, the raw moments can be computed as

As mentioned before, the mixed Poisson process has a stochastic occurrence rate Λ. Though the random variable Λ affects the probability of the number of arrivals, the inter-arrival times 

For the illustrative purpose, the random variable Λ follows a gamma distribution Λ ~ Γ(α, 

According to the renewal theorem [40] , Eq. (27) (28) where F W indicates the cumulative distribution function of the inter-arrival time. Consequently, the moments of long-term loss under the renewal process can be obtained by taking derivatives of Eq. (28) at zero. The mth-order derivative of the moment generating function gives 

When η equals zero, the mth-order moments of long-term loss is 

The first four moments of long-term loss can be effectively obtained by using this 

The third and fourth-order moments obtained from the renewal function are also validated by comparing with the values computed using the moment generating function, as shown in Eqs (14) and (15).

In summary, the developed moment generating function-based analytical approach can effectively evaluate the higher-order moments of long-term loss under different stochastic models. Based on the law of total expectation, the proposed approach expands the application scope of the moment approach, which was formerly used for the homogeneous Poisson Apart from the loss assessment, higher-order moments can be used to compute the longterm reliability [45] . For instance, skewness was involved in the third-order moment method to assess the long-term reliability of reinforced concrete structures under chloride-induced corrosion [46] . Reliability analysis involving skewness and kurtosis was also conducted in Lu et al. [47] . Another application of the moment generating function is that statistical moments can be used if there is insufficient information. For instance, Zhao and Lu [48] used statistical moments to describe probabilistic characteristics of random variables.

Two illustrative examples are presented to demonstrate the feasibility and applicability of the developed framework. The first example focuses on the computation of loss severity. The long-term seismic loss of the investigated bridge under seismic hazard is computed using the renewal model. The second example aims to identify the impact of climate variability on the hurricane-induced loss by using different stochastic processes. The homogeneous, nonhomogeneous, and mixed Poisson processes are used to model the hurricane arrivals.

According to the mathematical derivations provided in Section 3, there need several inputs for the long-term loss assessment, including the occurrence rate of hazard, monetary discount rate, time interval, and the statistical moments of loss severity under the investigated hazard scenario. Generally, the loss assessment of civil infrastructure under seismic hazards consists of four components, as shown in Figure 2 , including hazard analysis, vulnerability analysis, damage loss estimation, and long-term loss assessment. The loss severity is an essential element. This example aims to compute the loss severity of a reinforced concrete bridge under the seismic hazard. The recurrence of earthquakes is modeled as a renewal process. Five damage states of the bridge are specified as no damage, slight, moderate, major, and complete, according to HAZUS [53] . The seismic demand is the displacement ductility of the 

where Z  refers to the standard normal cumulative distribution function; S D and S C represent the median of seismic demand and capacity, respectively; and β D and β C refer to the standard deviation of demand and capacity, respectively. From the regression relationship, the value of A is 2.8869 and the value of B is 1.0702. Accordingly, the seismic fragility curves of four different damage states considering displacement ductility are obtained, as shown in Figure 4 . Given the structural vulnerability, the loss severity is the product of the repair cost of the damaged bridge and the conditional probability that the bridge is in a specified damage state under the given intensity measure [57] . The cost of rebuilding equals the unit price times the length and width of the bridge, in which the unit price of bridge rebuilding is taken as 2,306 USD/m 2 [58, 59] . The repair cost of the bridge being in different damage states can be evaluated as a fraction of the rebuilding cost through repair ratios. Repair ratios of 0, 0.03, 0.25, 0.75, and 1 are used for none, slight, moderate, major, and complete damage, respectively [60] . Subsequently, the loss severity equals the summation of loss associated with all possible damage states. The expectation, standard deviation, the third, and the fourth- 

The occurrence of hurricanes can be modeled as a stochastic process using the observation data [61, 32] . For instance, a homogeneous Poisson process can be used to model hurricane arrivals in a stationary environment [61] . However, stochastic models using historical observations only may not be sufficient to project future scenarios, as the variability in characteristics of hazards is not considered [33] . In recent decades, hurricane arrivals in a changing environment considering the effects of climate change and variability have been modeled by the non-homogeneous and mixed Poisson models [18, 31] . This example aims to assess the impact of climate change and variability on the hurricane-induced loss, from a long-term perspective.

The homogeneous Poisson process is widely used to model hurricane arrivals in a stationary environment, which assumes a constant occurrence rate λ based on historical observations [61, 32] . This rate is typically determined by dividing the total number of hurricane landfalls by the observation period [61] . Recent studies observe the increasing trend on hurricane frequency in the warming climate [62] . For such scenarios, the nonhomogeneous Poisson process with an increasing occurrence rate, i.e., λ(t int ), could be used to predict the increase in the number of hurricane arrivals [35, 18] .

In a changing environment, in addition to the potential increasing trend, the occurrence of hurricanes can be significantly influenced by climate variability. Climate variability refers to variations in the mean state and characteristics of climate [63] . Previous studies stated the importance of considering the occurrence rate as a random variable in the Poisson process (i.e., the mixed Poisson process) for future hurricane predictions. For instance, Elsner and Bossak [31] projected the occurrence rate of the U.S. hurricane landfalls using the mixed Poisson model. Villarini et al. [36] assessed changes in hurricane frequency using the mixed Poisson process, by modeling the dependence of hurricane occurrence on different climate indices. In this example, a gamma distributed stochastic rate is utilized [31] .

The long-term loss analysis is performed on a multi-span simply supported girder bridge [64] . The bridge has six spans equally distributed with a length of 146 m. This type of bridge 22 is most susceptible to the deck unseating damage [65] . The given annual occurrence rate of hurricane λ 0 is 0.245 for the investigated area in the stationary environment (i.e., with a homogeneous Poisson process). The rate is determined by counting the total number of 27 hurricane landfalls in the investigated region from 1900 to 2100, i.e., λ 0 = 27/110 per year [17] . The expected loss severity E[L] is 1.283 million USD [64] . The detailed computation of λ 0 and loss severity under hurricane hazards can be performed according to the process shown in the earthquake example. The monetary discount rate r is 2% for the long-term evaluation [66] . The loss severity follows an exponential distribution.

In a stationary environment, the occurrence of hurricanes is modeled as a homogeneous

Poisson process, with a rate of λ = λ 0 throughout the lifetime. In a changing climate, the nonhomogeneous Poisson process is adopted and the occurrence rate is assumed to follow a linear increasing relationship λ(t int ) = λ 0 (1 + ct int ), in which c refers to an annual increase rate of hurricane landfalls [35] . The annual increase rate of 0.2% indicates that the number of hurricanes is increased by 20% in the next century [8, 35] . The mixed Poisson process is also adopted to compute the long-term loss of the bridge under hurricanes considering climate variability. Herein, the parameters within the stochastic occurrence rate Λ ~ Γ(α, β) are based on the information presented in Elsner and Bossak [31] . The rate has a mean of E[Λ] = λ 0 (the same value as the rate in the homogeneous model) and the gamma parameters are assumed as Λ ~ Γ(0.49, 2). Given more information (e.g., climate information), the parameters used in the non-homogenous and mixed Poisson processes can be upgraded and the relevant results would be computed.

Using Eqs. (13) , (21) , and (26), the expectation, standard deviation, coefficient of variation, skewness, and kurtosis of the long-term loss under the homogeneous, nonhomogeneous, and mixed Poisson processes are obtained for the investigated bridge, as shown in Table 1 . When the service life reaches 150 years, the expected long-term loss of bridge under the non-homogeneous Poisson model is approximately increased by 5.0% in the changing climate. As a relatively small increase in the occurrence rate is assumed in this example, the loss results with respect to the non-homogeneous model do not show large differences from those of the homogeneous model.

According to Table 1 , all skewness values are positive (right skewed), which indicates a longer tail on the right side of the distribution of long-term loss. Homogeneous and nonhomogeneous models are moderately skewed (skewness between 0.5 and 1), while the mixed model is highly skewed (skewness greater than 1). Meanwhile, all kurtosis values are greater than 3, indicating that all the tails are heavier and longer than a normal distribution. For the case using the mixed Poisson process, the expected loss is the same as that using the homogeneous Poisson process, but the standard deviation, skewness, and kurtosis are much larger, which indicated potential heavy tail risks. Special attention should be paid to this aspect within the decision-making process.

The proposed analytical method is validated by the Monte Carlo simulation. In the simulation approach, the homogeneous model is generated using the exponentially distributed inter-arrival times, while a non-homogeneous process is simulated using the thinning method [67] . The mixed Poisson process is simulated based on the stochastic occurrence rate [68] .

The results of the mean and standard deviation of the long-term loss under the homogenous and non-homogeneous models are indicated in Figure 5 . Figure A parametric study is developed to measure the sensitivity of long-term loss to the monetary discount rate and variability of loss severity. A series of incremental monetary discount rates from 1% to 5% is considered for the homogeneous Poisson model. Figure 7(a) illustrates that the expected long-term loss decreases rapidly with the increase of the monetary discount rate. The same trend is also observed for the standard deviation of long- 25 term loss. However, the tendency is opposite with respect to the skewness and kurtosis. 26 assessed. The long-term loss with the initial standard deviation SD[L] is taken as the reference value. The long-term loss under the other four cases with triple, double, one-half, and one-third SD[L] is assessed. The associated standard deviations are presented in Figure   8 (a). It shows that the dispersion of long-term loss is significantly enlarged with the increase of loss severity. However, the skewness of long-term loss is not influenced, as shown in Figure 8 (b). From Eqs. (18) and (19), the skewness and kurtosis are independent of the exponentially distributed loss severity.

In the illustrative examples, structures are assumed to be fully restored to the pre-hazard stage. This assumption is commonly used in the loss assessment [14, 15, 69] . The level of restoration usually depends on many factors, e.g., the acceptable level of structural performance, investment, and tradeoffs between appropriate performance levels and investment. Further studies are needed to incorporate different restoration models.

This paper develops a novel moment generating function-based analytical approach to assess the long-term loss under both stationary and nonstationary hazards. Higher-order moments of long-term loss under four stochastic models are assessed, including the homogeneous Poisson process, non-homogeneous Poisson process, mixed Poisson process, and renewal process.

The proposed approach is applied to two illustrative examples to assess the loss severity and the long-term loss incorporating parametric studies. The analytical approach is validated by the Monte Carlo simulations. With a satisfactory agreement of the results, the effectiveness and accuracy of the proposed approach are validated. Following conclusions are drawn:

1. The moment generating function-based approach is proposed for the higher-order analysis of long-term loss. By using the law of total expectation, the developed approach successfully expands the application scope of the moment generating function. with limited information, e.g., with a few observations. The application of higher-order moments in the decision-making process needs to be investigated in future studies.

The derivation of the moment generating function of long-term loss under the homogeneous and non-homogeneous Poisson models in Eqs. (12) and (21) Eq. (A2) requires the joint probability density function of the arriving times T 1 , T 2 , …, T k . The loss severity is independent of the number of arrivals N(t int ). Hence, for the homogeneous Poisson process, the conditional joint probability density function of the arriving times T 1 , in which 0 < s 1 < s 2 < … < s n < t int . The result shown in Eq. (A3) can be alternatively explained by the order statistics [40, 70] . Given N(t int ) = n, the arriving time can be expressed by a sequence of independently identically uniformly distributed random variables {U 1 , 

Bridge life-cycle performance and cost: analysis, prediction, optimisation and decision-making

Probabilistic assessment of an interdependent healthcare: bridge network system under seismic hazard

Modelling seismic hazard in earthquake loss models with spatially distributed exposure

Probabilistic seismic loss assessment of a Vancouver high-rise building

A decision model for intergenerational lifecycle risk assessment of civil infrastructure exposed to hurricanes under climate change. Reliability Engineering and System Safety

A physically-based earthquake recurrence model for estimation of long-term earthquake probabilities

Longterm time-dependent probabilities for the third Uniform California Earthquake Rupture Forecast (UCERF3)

Modeled impact of anthropogenic warming on the frequency of intense Atlantic hurricanes

Climatology models for extreme hurricane winds near the United States

Non-stationary frequency analysis of heavy rainfall events in southern France

Scenario-based projection of extreme sea levels

Life-cycle management of deteriorating civil infrastructure considering resilience to lifetime hazards: A general approach based on renewal-reward processes

Probabilistic time-dependent multihazard lifecycle assessment and resilience of bridges considering climate change

Building life-cycle cost analysis due to mainshock and aftershock occurrences

Minimum building life-cycle cost design criteria. I: Methodology

Stochastic characterization and decision bases under time-dependent aftershock risk in performance-based earthquake engineering

Assessing hurricane damage costs in the presence of vulnerability model uncertainty

Dealing with hurricane surge flooding in a changing environment: part I. Risk assessment considering storm climatology change, sea level rise, and coastal development

Stochastic renewal process models for estimation of damage cost over the life-cycle of a structure

Stochastic dominance: Investment decision making under uncertainty

Optimal seismic design considering risk attitude, societal tolerable risk level, and life quality criterion

Mean-variance analysis in portfolio choice and capital markets

Tail risk and asset prices

6 Finance in the times of coronavirus

A polynomial goal programming model for portfolio optimization based on entropy and higher moments

Global optimization of higher order moments in portfolio selection

Modeling for seasonal marked point processes: An analysis of evolving hurricane occurrences

Global multivariate point pattern models for rain type occurrence

Long-term earthquake clustering

Optimization and risk acceptability based on the life quality index

Bayesian analysis of US hurricane climate

Stochastic modeling of hurricane damage

The economics of climate change impacts and policy benefits at city scale: a conceptual framework

A Brownian model for recurrent earthquakes

Managing risks to civil infrastructure due to natural hazards: communicating long-term risks due to climate change

Modeling the dependence of tropical storm counts in the North Atlantic basin on climate indices

Temporal and magnitude dependence in earthquake recurrence models

Viscoelasticity, postseismic slip, fault interactions, and the recurrence of large earthquakes

Estimating background activity based on interevent-time distribution

Introduction to probability models

Stochastic calculus for finance II: Continuous-time models

Statistics of extremes, with applications in environment, insurance, and finance. Extreme Values in Finance, Telecommunications, and the Environment

Hazard function analysis for flood planning under nonstationarity

Long-term resilience and loss assessment of highway bridges under multiple natural hazards

Moment methods for structural reliability

Time-dependent probability assessment for chloride induced corrosion of RC structures using the third-moment method

Fast integration algorithms for time-dependent structural reliability analysis considering correlated random variables

Fourth-moment standardization for structural reliability assessment

Frequency of earthquakes in California

United States Geological Survey (USGS)

NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes

Open System for Earthquake Engineering Simulation (OpenSees

Earthquake loss estimation methodology: technical and user manuals. Federal Emergency Management Agency

Analytical fragility curves for highway bridges in moderate seismic zones (Doctoral dissertation

New ground motion selection procedures and selected motions for the PEER transportation research program

Hybrid multi-criteria decision making for seismic intensity measure selection considering uncertainty

Seismic loss and resilience assessment of singlecolumn rocking bridges

Resilience and life-cycle performance of smart bridges with shape memory alloy (SMA)-cable-based bearings

Performance-based assessment of bridges with steel-SMA reinforced piers in a life-cycle context by numerical approach

Redars 2 methodology and software for seismic risk analysis of Highway Systems

Secular changes to the ENSO-US hurricane relationship

Increasing destructiveness of tropical cyclones over the past 30 years

A globally consistent reanalysis of hurricane variability and trends

Risk-informed hazard loss of bridges in a life-cycle context

Experimental and 3D Numerical Investigation of Solitary Wave Forces on Coastal Bridges

Ethical discounting for civil infrastructure decisions extending over multiple generations

Transforming renewal processes for simulation of nonstationary arrival processes

Modeling of the Risk Process

Risk-based seismic life-cycle costbenefit (LCC-B) analysis for bridge retrofit assessment

Non-life insurance mathematics: an introduction with the Poisson process

Conditional, non-homogeneous and doubly stochastic compound Poisson processes with stochastic discounted claims

The study has been supported by The Hong Kong Polytechnic University under Start-Up