Next Article in Journal
The Population Accuracy Index: A New Measure of Population Stability for Model Monitoring
Next Article in Special Issue
Special Issue “Risk, Ruin and Survival: Decision Making in Insurance and Finance”
Previous Article in Journal
The Optimum Leverage Level of the Banking Sector
Previous Article in Special Issue
Practice Oriented and Monte Carlo Based Estimation of the Value-at-Risk for Operational Risk Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Risk Measures and Rate of Spatial Diversification

EPFL, Chair of Statistics STAT, EPFL-SB-MATH-STAT, MA B1 433 (Bâtiment MA), Station 8, 1015 Lausanne, Switzerland
Risks 2019, 7(2), 52; https://doi.org/10.3390/risks7020052
Submission received: 21 December 2018 / Revised: 25 March 2019 / Accepted: 4 April 2019 / Published: 2 May 2019
(This article belongs to the Special Issue Risk, Ruin and Survival: Decision Making in Insurance and Finance)

Abstract

:
An accurate assessment of the risk of extreme environmental events is of great importance for populations, authorities and the banking/insurance/reinsurance industry. Koch (2017) introduced a notion of spatial risk measure and a corresponding set of axioms which are well suited to analyze the risk due to events having a spatial extent, precisely such as environmental phenomena. The axiom of asymptotic spatial homogeneity is of particular interest since it allows one to quantify the rate of spatial diversification when the region under consideration becomes large. In this paper, we first investigate the general concepts of spatial risk measures and corresponding axioms further and thoroughly explain the usefulness of this theory for both actuarial science and practice. Second, in the case of a general cost field, we give sufficient conditions such that spatial risk measures associated with expectation, variance, value-at-risk as well as expected shortfall and induced by this cost field satisfy the axioms of asymptotic spatial homogeneity of order 0, −2, −1 and −1, respectively. Last but not least, in the case where the cost field is a function of a max-stable random field, we provide conditions on both the function and the max-stable field ensuring the latter properties. Max-stable random fields are relevant when assessing the risk of extreme events since they appear as a natural extension of multivariate extreme-value theory to the level of random fields. Overall, this paper improves our understanding of spatial risk measures as well as of their properties with respect to the space variable and generalizes many results obtained in Koch (2017).

1. Introduction

Hurricane Irma, which affected many Caribbean islands and parts of Florida in September 2017 caused at least 134 deaths and catastrophic damage exceeding 64.8 billion USD in value. Such an example shows the prime importance for civil authorities and for the insurance1 industry of the accurate assessment of the risk of natural disasters, particularly as, in a climate change context, certain types of extreme events become more and more frequent (e.g., Bevere and Mueller 2014).
Motivated by the spatial feature of natural disasters, Koch (2017) introduced a new notion of spatial risk measure, which makes explicit the contribution of the space and enables one to account for at least part of the spatial dependence in the risk measurement. He also introduced a set of axioms describing how the risk is expected to evolve with respect to the space variable, at least under some conditions. These notions constitute relevant tools for risk assessment. For instance, the knowledge of the order of asymptotic spatial homogeneity allows the quantification of the rate of spatial diversification. Hence, they may be appealing for the banking/insurance industry. It should be highlighted that the literature about risk measures in a spatial context is very limited. To the best of our knowledge, the paper by Koch (2017) constitutes the first attempt to establish a theory about risk measures in a spatial context where the risks spread over a continuous geographical region.
In the following, the spatial risk measure associated with a classical risk measure Π and induced by a cost random field C (e.g., modelling the cost due to damage caused by a natural disaster) consists in the function of space arising from the application of Π to the normalized integral of C on various geographical areas. The contribution of this paper is threefold. First, we further explore the notions of spatial risk measure and corresponding axioms introduced in Koch (2017). Among others, we show that, for a given region, the distribution of the normalized spatially aggregated loss is entirely determined by the finite-dimensional distributions of the cost field, and propose alternative definitions of the concepts developed in Koch (2017). Additionally, we deeply explain why this whole theory about spatial risk measures is fruitful for both actuarial science and practice; e.g., we show that considering the risk related to the normalized loss does not prevent our theory from being successful for the study of the risk related to the non-normalized loss. We also point out how it can be used by insurance companies to tackle concrete issues. Second, in the case of a general cost field, we give sufficient conditions such that spatial risk measures associated with expectation, variance, value-at-risk (VaR) as well as expected shortfall (ES) and induced by this cost field satisfy the axiom of asymptotic spatial homogeneity of order 0, −2, −1 and −1, respectively. Last but not least, we focus on the case where the cost field is a function of a max-stable random field. We provide sufficient conditions on both the function and the max-stable field such that spatial risk measures associated with expectation, variance, VaR as well as ES and induced by the resulting cost field satisfy the axiom of asymptotic spatial homogeneity of order 0, −2, −1 and −1, respectively. Max-stable random fields naturally appear when one is interested in extreme events having a spatial extent since they constitute an extension of multivariate extreme-value theory to the level of random fields (in the case of stochastic processes, see, e.g., de Haan 1984; de Haan and Ferreira 2006). They are particularly well suited to model the temporal maxima of a given variable (for instance a meteorological variable) at all points in space since they arise as the pointwise maxima taken over an infinite number of appropriately rescaled independent and identically distributed random fields. On the whole, this study improves our comprehension of spatial risk measures as well as of their properties with respect to the space variable and generalizes many results by Koch (2017).
The remainder of the paper is organized as follows. In Section 2, we recall and further study the notion of spatial risk measure and the corresponding set of axioms introduced in Koch (2017). Furthermore, we thoroughly demonstrate their usefulness for both actuarial science and practice. Then, we introduce some concepts about mixing and central limit theorems for random fields. Finally, we provide some insights about max-stable random fields. Then, Section 3 presents our results relating to the properties of some spatial risk measures. Finally, Section 4 contains a short summary as well as some perspectives.
Throughout the paper, ( Ω , F , P ) is an adequate probability space and = d and d designate equality and convergence in distribution, respectively. In the case of random fields, distribution has to be understood as the set of all finite-dimensional distributions. Finally, we denote by ν the Lebesgue measure.

2. Spatial Risk Measures and Other Concepts

2.1. Spatial Risk Measures and Corresponding Axioms

First, we describe the setting required for a proper definition of spatial risk measures. Let A be the set of all compact subsets of R 2 with a positive Lebesgue measure and A c be the set of all convex elements of A . Denote by C the set of all real-valued and measurable2 random fields on R 2 having almost surely (a.s.)3 locally integrable sample paths. Let P be the family of all possible distributions of random fields belonging to C . Each random field represents the economic or insured cost caused by the events belonging to specified classes and occurring during a given time period, say [ 0 , T L ] . In the following, T L is considered as fixed and does not appear anymore for the sake of notational parsimony. Each class of events (e.g., European windstorms or hurricanes) will be referred to as a hazard in the following. Let L be the set of all real-valued random variables defined on ( Ω , F , P ) . A risk measure typically will be some function Π : L R . This kind of risk measure will be called a classical risk measure in the following. A classical risk measure Π is termed law-invariant if, for all X ˜ L , Π ( X ˜ ) only depends on the distribution of X ˜ .
We first remind the reader of the definition of the normalized spatially aggregated loss, which enables one to disentangle the contribution of the space and the contribution of the hazards and underpins our definition of spatial risk measure.
Definition 1
(Normalized spatially aggregated loss as a function of the distribution of the cost field). For A A and P P , the normalized spatially aggregated loss is defined by
L N ( A , P ) = 1 ν ( A ) A C P ( x ) ν ( d x ) ,
where the random field { C P ( x ) } x R 2 belongs to C and has distribution P.
The quantity
L ( A , P ) = A C P ( x ) ν ( d x )
corresponds to the total economic or insured loss over region A due to specified hazards. For technical reasons and to favour a more intuitive understanding, we base our definition of spatial risk measures on L N ( A , P ) , which is the loss per surface unit and can be interpreted, in a discrete setting4 and in an insurance context, as the mean loss per insurance policy. Among other advantages, this normalization enables a fair comparison of the risks related to regions having different sizes.
Since the field C P is measurable, L ( A , P ) and L N ( A , P ) are well-defined random variables. Moreover, they are a.s. finite as A is compact and C P has a.s. locally integrable sample paths. The following proposition gives a sufficient condition for a random field to have a.s. locally integrable sample paths.
Proposition 1.
Let d 1 and { Q ( x ) } x R d be a measurable random field. If the function
E : R d R x E [ | Q ( x ) | ]
is locally integrable, then Q has a.s. locally integrable sample paths.
Proof. 
Let A be a compact subset of R d . First, since Q is measurable, A | Q ( x ) | ν ( d x ) is a well-defined random variable. By Fubini’s theorem, we have
E A | Q ( x ) | ν ( d x ) = A E [ | Q ( x ) | ] ν ( d x ) < ,
which necessarily implies that
A | Q ( x ) | ν ( d x ) < a . s .
Since this is true for all A being a compact subset of R d , we obtain the result.  □
We now recall the notion of spatial risk measure introduced by Koch (2017), which makes explicit the contribution of the space in the risk measurement.
Definition 2
(Spatial risk measure as a function of the distribution of the cost field). A spatial risk measure is a function R Π that assigns a real number to any region A A and distribution P P :
R Π : A × P R ( A , P ) Π ( L N ( A , P ) ) ,
where Π is a classical and law-invariant risk measure and L N ( A , P ) is defined in (1).
This extends the notion of classical risk measure to the spatial and infinite-dimensional setting as we now have a function of both the space and the distribution of a random field (or directly a random field, see below) instead of a function of a unique real-valued random variable. Note that law-invariance of Π is necessary for spatial risk measures to be defined in this way; see below for more details. For a given Π and a fixed P P , the quantity R Π ( · , P ) is referred to as the spatial risk measure associated with Π and induced by P. A nice feature is that, for many useful classical risk measures Π such as, e.g., variance, VaR and ES, this notion of spatial risk measure allows one to take (at least) part of the spatial dependence structure of the field C P into account. We could define spatial risk measures in the same way but using the non-normalized spatially aggregated loss; this is not what we do for reasons explained above and in Remark 2 below.
Now, we remind the reader of the set of axioms for spatial risk measures developed in Koch (2017). It concerns the spatial risk measures properties with respect to the space and not to the cost distribution, the latter being considered as given by the problem at hand. For any A A , let b A denote its barycenter.
Definition 3
(Set of axioms for spatial risk measures induced by a distribution). Let Π be a classical and law-invariant risk measure. For a fixed P P , we define the following axioms for the spatial risk measure associated with Π and induced by P, R Π ( · , P ) :
1. 
Spatial invariance under translation:
for all v R 2 and A A , R Π ( A + v , P ) = R Π ( A , P ) , where A + v denotes the region A translated by the vector v .
2. 
Spatial sub-additivity:
for all A 1 , A 2 A , R Π ( A 1 A 2 , P ) min { R Π ( A 1 , P ) , R Π ( A 2 , P ) } .
3. 
Asymptotic spatial homogeneity of order γ , γ 0 :
for all A A c ,
R Π ( λ A , P ) = λ K 1 ( A , P ) + K 2 ( A , P ) λ γ + o 1 λ γ ,
where λ A is the area obtained by applying to A a homothety with center b A and ratio λ > 0 , and K 1 ( · , P ) : A c R , K 2 ( · , P ) : A c R \ { 0 } are functions depending on P.
It is also reasonable to introduce the axiom of spatial anti-monotonicity: for all A 1 , A 2 A , A 1 A 2 R Π ( A 2 , P ) R Π ( A 1 , P ) . The latter is equivalent to the axiom of spatial sub-additivity. These axioms appear natural and make sense at least under some conditions on the cost field C P (e.g., stationarity5 in the case of spatial invariance under translation and spatial sub-additivity) and for some classical risk measures Π . The axiom of spatial sub-additivity indicates spatial diversification. If it is satisfied with strict inequality, an insurance company would be well advised to underwrite policies in both regions A 1 and A 2 instead of only one of them. This axiom involves the minimum operator because the concept of spatial risk measure is based on the normalized spatially aggregated loss; using the summation operator instead would not provide information about spatial diversification. On the other hand, if spatial risk measures were defined using the non-normalized loss, then summation would make sense; see Remark 2 below for more details. Originally, Koch (2017) used the term “sub-additivity”, among other reasons, by analogy with the axiom of sub-additivity by Artzner et al. (1999), which also conveys a diversification idea. The axiom of asymptotic spatial homogeneity of order γ quantifies the rate of spatial diversification when the region becomes large. Consequently, determining the value of γ is of interest for the insurance industry; see Section 2.2 for further details.
The axioms of spatial invariance under translation and spatial sub-additivity a priori make sense only if the cost field satisfies at least some kind of stationarity. If an insurance company covers a region A 1 which is much less risky than a region A 2 , it is very unlikely that the company reduces its risk by covering A 1 A 2 . For a given hazard (e.g., hurricanes), the cost resulting from a single specific event (e.g., a particular hurricane) generally varies across space, making any particular realization of the cost field spatially inhomogeneous. Nevertheless, the cost field (and not one realization of it) related to this hazard can be stationary or, at least, piecewise stationary; see immediately below.
In concrete actuarial applications, the cost field (for a given hazard) is often non-stationary over the entire region covered by the insurance company, unless it is a very small area. In many cases, however, it can reasonably be considered as locally stationary; see, e.g., Dahlhaus (2012) for an excellent review about locally stationary processes, and Eckley et al. (2010) as well as Anderes and Stein (2011) for papers dealing with local non-stationarity in the case of random fields. Locally stationary processes can be well approximated by piecewise stationary processes (e.g., Ombao et al. 2001, Section 2.2) and, assuming this to be also true for random fields, we can reasonably consider the cost field to be stationary over sub-regions, at least in most cases. In the latter, the axioms of spatial invariance under translation and spatial sub-additivity make sense separately on each sub-region over which the field is stationary. Let S ub be such a sub-region (a subset of R 2 ) and S ub be the set of all compact subsets of S ub with a positive Lebesgue measure. The axiom of spatial invariance under translation becomes: for all v R 2 and A S ub such that A + v S ub , R Π ( A + v , P ) = R Π ( A , P ) ; spatial sub-additivity is now written: for all A 1 , A 2 S ub , R Π ( A 1 A 2 , P ) min { R Π ( A 1 , P ) , R Π ( A 2 , P ) } .
Of course, the fact that the axioms of Definition 3 are satisfied depends on both the classical risk measure Π and the cost field C P . It may be interesting to determine for which classical risk measures the axioms are satisfied for the broadest class of cost fields. These classical risk measures could be considered as “adapted” to the spatial context.
Remark 1.
Although the concept of spatial risk measure and related axioms naturally apply in an insurance context (see Section 2.2 for further details), they can also be used in the banking industry and on financial markets. A potential application is the assessment of the risk related to event-linked securities such as catastrophe bonds. Furthermore, they can be used for a wider class of risks than those linked with damage due to environmental events. These concepts are actually insightful as soon as the risks spread over a geographical region. One might think, e.g., about the loss in value of real estate due to adverse economic conditions.
We close this section by deeply commenting on the previous concepts and giving slightly modified and more natural versions of previous definitions. First, we need the following useful result.
Theorem 1.
Let d 1 and { H ( x ) } x R d be a measurable random field having a.s. locally integrable sample paths. Moreover, let A be a compact subset of R d with positive Lebesgue measure. Then the distribution of
L N ( A , H ) = 1 ν ( A ) A H ( x ) ν ( d x )
only depends on A and the finite-dimensional distributions of H.
Proof. 
The proof is partly inspired from the proof of Theorem 11.4.1 in Samorodnitsky and Taqqu (1994). We assume that the random field H is defined on the probability space ( Ω , F , P ) . For a fixed ω Ω , we denote by H ω the corresponding realization of H on R d and by H ω ( x ) the realization of H at location x . By definition, we have, for almost every ω Ω , that
L N ( A , H ω ) = 1 ν ( A ) A H ω ( x ) ν ( d x ) .
Now, let ( Ω 1 , F 1 , P 1 ) be a probability space different from the probability space ( Ω , F , P ) . Let U be a random vector defined on ( Ω 1 , F 1 , P 1 ) and following the uniform distribution on A, with density f U ( x ) = I { x A } / ν ( A ) , x R d . From (3), it directly follows that, for almost every ω Ω ,
L N ( A , H ω ) = R d H ω ( x ) f U ( x ) ν ( d x ) .
Let us denote by E 1 the expectation with respect to the probability measure P 1 . We have
E 1 [ H ω ( U ) ] = R d H ω ( x ) f U ( x ) ν ( d x ) ,
giving, using (4),
L N ( A , H ω ) = E 1 [ H ω ( U ) ] .
Now, let U 1 , , U n be independent replications of U (which are independent of the random field H). The strong law of large numbers gives that, for almost every ω Ω ,
L N ( A , H ω ) = lim n 1 n i = 1 n H ω ( U i ) P 1 a . s .
Therefore, using Fubini’s theorem, we deduce that, for P 1 -almost every ω 1 Ω 1 ,
L N ( A , H ω ) = lim n 1 n i = 1 n H ω ( U i ( ω 1 ) ) P a . s .
Now, we choose ω 0 Ω 1 such that the (non-random) sequence ( U 1 ( ω 0 ) , U 2 ( ω 0 ) , ) satisfies (6). We obtain
L N ( A , H ω ) = lim n 1 n i = 1 n H ω ( U i ( ω 0 ) ) P a . s .
Equation (7) says that the distribution of L N ( A , H ) is determined by the finite-dimensional distributions at the points belonging to the set { U i ( ω 0 ) : i N } . This yields the result. □
It is more natural, especially in terms of interpretation, to introduce the normalized spatially aggregated loss as a function of the cost field instead of its distribution, as shown immediately below.
Definition 4
(Normalized spatially aggregated loss as a function of the cost field). The normalized spatially aggregated loss function is defined by
L N : A × C R ( A , C ) 1 ν ( A ) A C ( x ) ν ( d x ) .
Let C P C be a random field with distribution P. Although a particular realization of L N ( A , C P ) obviously depends on C P (through its corresponding realization), we know from Theorem 1 that its distribution is entirely characterized by A and P. This explains our notation L N ( A , P ) instead of L N ( A , C P ) in Definition 1. More precisely, let C P ( 1 ) , C P ( 2 ) C be random fields having the same distribution P. Then, C P ( 1 ) and C P ( 2 ) have the same finite-dimensional distributions, which implies that L N ( A , C P ( 1 ) ) = d L N ( A , C P ( 2 ) ) .
Similarly, it can appear more natural to define spatial risk measures as functions of the cost field instead of its distribution. Moreover, this allows spatial risk measures to be defined even when the classical risk measure Π is not law-invariant.
Definition 5
(Spatial risk measure as a function of the cost field). A spatial risk measure is a function R Π that assigns a real number to any region A A and random field C C :
R Π : A × C R ( A , C ) Π ( L N ( A , C ) ) ,
where Π is a classical risk measure.
For a given classical and law-invariant risk measure Π and a given region A A , the value of the spatial risk measure of Definition 5 is completely determined by the distribution of L N ( A , C ) by law-invariance of Π . Consequently, using Theorem 1, it is completely determined by A and the distribution of the cost field C. This explains why Koch (2017) has introduced the notion of spatial risk measure as a function of the distribution of C (see the reminder in Definition 2); if Π is law-invariant, the spatial risk measures described in Definitions 2 and 5 refer to the same notion. For a given Π and a fixed C C , R Π ( · , C ) is referred to as the spatial risk measure associated with Π and induced by C.
Of course, we can also express the axioms recalled in Definition 3 for the spatial risk measures induced by a cost field C C introduced immediately above. On top of being more natural, it enables one to leave out the assumption of law-invariance for the classical risk measure Π .
Definition 6
(Set of axioms for spatial risk measures induced by a cost field). Let Π be a classical risk measure. For a fixed C C , we define the following axioms for the spatial risk measure associated with Π and induced by C, R Π ( · , C ) :
1. 
Spatial invariance under translation:
for all v R 2 and A A , R Π ( A + v , C ) = R Π ( A , C ) , where A + v denotes the region A translated by the vector v .
2. 
Spatial sub-additivity:
for all A 1 , A 2 A , R Π ( A 1 A 2 , C ) min { R Π ( A 1 , C ) , R Π ( A 2 , C ) } .
3. 
Asymptotic spatial homogeneity of order γ , γ 0 :
for all A A c ,
R Π ( λ A , C ) = λ K 1 ( A , C ) + K 2 ( A , C ) λ γ + o 1 λ γ ,
where λ A is the area obtained by applying to A a homothety with center b A and ratio λ > 0 , and K 1 ( · , C ) : A c R , K 2 ( · , C ) : A c R \ { 0 } are functions depending on C.
All the explanations and interpretations given for Definitions 1–3 remain valid in the case of Definitions 4–6. For the reasons mentioned above, our opinion is that Definitions 4–6 rather than previous ones should be used. This is what is done in the following.

2.2. Concrete Applications to Insurance

This section is dedicated to the connections between the concepts described above and actuarial risk theory as well as real insurance practice. We especially show how they can be used for concrete purposes. In an insurance context, the quantity
L ( A , C ) = A C ( x ) ν ( d x )
appearing in Definition 4 (or equivalently in (2)) can be seen as a continuous and more complex version of the classical actuarial individual risk model. The latter can be formulated as
L ind = i = 1 N X i ,
where L ind is the total loss, N denotes the number of insurance policies and, for i = 1 , , N , X i is the claim related to the i-th policy. The X i ’s are generally assumed to be independent but not necessarily identically distributed. In L ( A , C ) , each location x corresponds to a specific insurance policy and thus each C ( x ) is equivalent to a X i in (10). By the way, by choosing ν to be a counting measure instead of the Lebesgue measure, the integral in (9) can be reduced to a sum, e.g., x A C ( x ) , where A is a finite set of locations in R 2 (e.g., part of a lattice in Z 2 ). It is worth mentioning that the ideas of this paper can easily be applied to such a framework.
Even if dependence between the X i , i = 1 , , N , in (10) was allowed, considering L ( A , C ) (see (9)) would appear more promising. Indeed, the geographical information of each risk (i.e., insurance policy) is explicitly taken into account and, consequently, the dependence between all risks can be modelled in a more realistic way than in (10). The dependence between the risks directly inherits from their respective associated geographical positions and, thus, ignoring their localizations as in (10) makes the modelling of their dependence more arbitrary and likely less reliable. In our approach, this dependence is fully characterized by the spatial dependence structure of the cost field C. Potential central limit theorems (see below) would have stronger implications because the dependence is more realistic. For these reasons, Models (8) and (9) allow a more accurate assessment of spatial diversification. The same remarks hold if we compare our loss models with the classical actuarial collective risk model.
Our risk models (8) and (9) and more generally our theory about spatial risk measures may be particularly relevant for an insurance company willing to adapt its policies portfolio. For example, the axioms of spatial sub-additivity and asymptotic spatial homogeneity can help it to assess the potential relevance of extending its activity to a new geographical region. Such an analysis requires the company to have an accurate view of the dependence between its risks (inter alia between the possible new risks and those already present in the portfolio), as allowed by Models (8) and (9) through the cost field C. Model (10) would not enable the insurer to precisely account for the dependence between the new risks and those already in the portfolio and hence to properly quantify the impact of a geographical expansion, i.e., of an increase of the number of contracts N.
At present, we show that, consistently with our intuition, considering the risk related to the normalized spatially aggregated loss is also insightful when the insurer is interested in the risk related to its non-normalized counterpart, which is often the case. Let Π be a positive homogeneous and translation invariant classical risk measure and p r denote either the claims reserves, revenues or any relevant related quantity6 per surface unit (possibly the mean premium per surface unit) of an insurance company Ins.
We first consider the axiom of spatial sub-additivity, which is assumed to be satisfied. Ins covers region A 1 for a given hazard and potentially aims at covering also a region A 2 disjoint of A 1 . We assume that Ins properly hedges its risk on A 1 , i.e.,
ν ( A 1 ) p r Π ( L ( A 1 , C ) ) , i . e . , p r Π ( L N ( A 1 , C ) ) ,
by positive homogeneity. Using again the same property,
Π ( L ( A 1 A 2 , C ) ) = ν ( A 1 A 2 ) Π ( L N ( A 1 A 2 , C ) ) .
Combined with
Π ( L N ( A 1 A 2 , C ) ) Π ( L N ( A 1 , C ) ) ,
this yields
Π ( L ( A 1 A 2 , C ) ) ν ( A 1 A 2 ) ν ( A 1 ) Π ( L ( A 1 , C ) ) .
Hence, by translation invariance,
Π ( L ( A 1 A 2 , C ) ν ( A 1 A 2 ) p r ) = Π ( L ( A 1 A 2 , C ) ) ν ( A 1 A 2 ) p r ν ( A 1 A 2 ) ν ( A 1 ) Π ( L ( A 1 , C ) ) ν ( A 1 A 2 ) p r .
It follows from (11) that
p r [ ν ( A 1 A 2 ) ν ( A 1 ) ] Π ( L ( A 1 , C ) ) ν ( A 1 ) [ ν ( A 1 A 2 ) ν ( A 1 ) ] ,
which gives
ν ( A 1 A 2 ) ν ( A 1 ) Π ( L ( A 1 , C ) ) ν ( A 1 A 2 ) p r Π ( L ( A 1 , C ) ) ν ( A 1 ) p r .
The combination of (12) and (13) yields that
Π ( L ( A 1 A 2 , C ) ν ( A 1 A 2 ) p r ) Π ( L ( A 1 , C ) ν ( A 1 ) p r ) .
The last inequality is strict if that in the axiom of spatial sub-additivity or in (11) is so. Thus, if Ins suitably hedges its risk on A 1 , the risk is even better hedged on A 1 A 2 . Exactly the same reasoning holds for A 2 .
Remark 2.
For spatial risk measures defined using the non-normalized spatially aggregated loss, we could propose the following axiom of spatial sub-additivity: for all disjoint A 1 , A 2 A , Π ( L ( A 1 A 2 , C ) ) Π ( L ( A 1 , C ) ) + Π ( L ( A 2 , C ) ) . Nevertheless, this property is trivially satisfied as soon as the classical risk measure Π is sub-additive and therefore its validity does not depend on the properties of the cost field C. Basing the axiom of spatial sub-additivity on the normalized spatially aggregated loss as we did is more appealing since it allows a diversification effect coming from C (and not only from Π ). This argument is in favour of defining spatial risk measures using the normalized spatially aggregated loss.
We now consider the axiom of asymptotic spatial homogeneity of order γ , γ 0 . Assume that it is satisfied with γ > 0 (e.g., we will see that for Π being VaR or ES, γ typically equals 1). It follows from Definition 6, Point 3, that
Π ( L ( λ A , C ) ν ( λ A ) p r ) = λ λ 2 ν ( A ) K 1 ( A , C ) + ν ( A ) K 2 ( A , C ) λ 2 γ + o λ 2 γ λ 2 ν ( A ) p r = λ λ 2 ν ( A ) ( K 1 ( A , C ) p r ) + ν ( A ) K 2 ( A , C ) λ 2 γ + o λ 2 γ .
Since γ > 0 , the dominant term as λ is λ 2 ν ( A ) ( K 1 ( A , C ) p r ) . Assume that K 1 ( A , C ) > 0 and K 2 ( A , C ) > 0 . This is true under the conditions of Section 3 for VaR and ES: we have K 1 ( A , C ) = E [ C ( 0 ) ] , which is positive as the cost field can be assumed to be non-negative and not a.s. equal to 0; regarding K 2 ( A , C ) , this is always true for ES and, provided that the confidence level α is greater than 1 / 2 , also for VaR. Consequently, for λ large enough, the total risk of the company, Π ( L ( λ A , C ) ν ( λ A ) p r ) , is a decreasing function of λ as soon as the revenue per surface unit (or claims reserves, …) satisfies p r > K 1 ( A , C ) . Under the conditions of Section 3, for VaR and ES, K 1 ( A , C ) = E [ C ( x ) ] for all x R 2 , and therefore the latter inequality entails that the revenue per surface unit (e.g., the mean premium) exceeds the expected cost at each location, which appears natural. The term 2 γ corresponds to the second highest power with respect to λ . Provided that K 2 ( A , C ) > 0 and 0 < γ < 2 (which is true for VaR and ES under the conditions of Section 3), the corresponding term, ν ( A ) K 2 ( A , C ) λ 2 γ , increases the total risk of the company as λ increases. However, the highest the value of γ , the fastest the decrease of the total risk as λ increases owing to the term in λ 2 . For λ large, the values of γ , K 1 ( A , C ) , K 2 ( A , C ) and p r allow one to determine the value of λ necessary to reach a targeted sufficiently low level of the total risk. Note that in the case of the variance, at least under the conditions of Section 3, K 1 ( A , C ) = 0 and γ = 2 .
Remark 3.
The axioms of spatial invariance under translation and asymptotic spatial homogeneity could also be defined for spatial risk measures based on the non-normalized spatially aggregated loss. Spatial invariance under translation would be unchanged and asymptotic spatial homogeneity of order γ , γ 0 , would become: for all A A c ,
Π ( L ( λ A , C ) ) = λ λ 2 ν ( A ) K 1 ( A , C ) + ν ( A ) K 2 ( A , C ) λ 2 γ + o λ 2 γ .
In this case, we would obtain the risk related to the non-normalized loss without assuming that Π is positive homogeneous.
Finally, we discuss a possible way for a company to develop an adequate model for the cost field C in regions where it is still inactive. The general model for the cost field introduced in Koch (2017, Section 2.3), is written
C ( x ) x R 2 = E ( x ) D Z ( x ) x R 2 ,
where { E ( x ) } x R 2 is the exposure field, D a damage function and { Z ( x ) } x R 2 the random field of the environmental variable generating risk. The cost is assumed to be only due to a unique class of events, i.e., to a unique natural hazard. The latter (e.g., heat waves or hurricanes) is described by the random field of an environmental variable (e.g., the temperature or the wind speed, respectively), Z. We assume that Z is representative of the risk during the whole period [ 0 , T L ] . The application of the damage function (also referred to as vulnerability curve in the literature) D to the natural hazard random field gives the destruction percentage at each location. Finally, multiplying the destruction percentage by the exposure gives the cost at each location. For more details, we refer the reader to Koch (2017, Section 2.3). In order to obtain an adequate model C in regions where it has no policies yet, the company can for instance consider crude estimates of the exposure field in the new region, develop a detailed statistical model7 for the environmental field Z responsible of the risk insured (e.g., wind speed in the case of hurricanes) using appropriate data and apply the same damage functions as in the region it already covers. The company can then simulate from this cost model, hence obtaining an empirical distribution of the loss appearing in (8) and (9). This makes it possible to check whether the axiom of spatial sub-additivity is satisfied or not. Furthermore, if the spatial domain is large (which is generally the case for reinsurance companies), considering potential central limit theorems and determining the order of asymptotic spatial homogeneity (by checking if the conditions of Section 3 are satisfied) is useful as it allows the company to quantify the rate of spatial diversification.
Remark 4.
Strictly speaking, the terms of the insurance policies should be accounted for in Model (14). By the way, the latter model can be interpreted differently from what is done here. For instance, we can imagine that Z represents the random field of the real cost and D accounts for the terms of the policies.

2.3. Mixing and Central Limit Theorems for Random Fields

We first remind the reader of the definition of the α - and β -mixing coefficients which will be used in Section 3. Let { X ( x ) } x R d be a real-valued random field. For S R d a closed subset, we denote by F S X the σ -field generated by the random variables { X ( x ) : x S } . Let S 1 , S 2 R d be disjoint closed subsets. The α -mixing coefficient (introduced by Rosenblatt 1956) between the σ -fields F S 1 X and F S 2 X is defined by
α X ( S 1 , S 2 ) = sup | P ( A B ) P ( A ) P ( B ) | : A F S 1 X , B F S 2 X .
The β -mixing coefficient or absolute regularity coefficient (attributed to Kolmogorov in Volkonskii and Rozanov 1959) between the σ -fields F S 1 X and F S 2 X is given by
β X ( S 1 , S 2 ) = 1 2 sup i = 1 I j = 1 J | P ( A i B j ) P ( A i ) P ( B j ) | ,
where the supremum is taken over all partitions { A 1 , , A I } and { B 1 , , B J } of Ω with the A i ’s in F S 1 X and the B j ’s in F S 2 X . These coefficients satisfy the useful inequality
α X ( S 1 , S 2 ) 1 2 β X ( S 1 , S 2 ) , for   all   S 1 , S 2 R d .
Now, we recall the concepts of Van Hove sequence and central limit theorem (CLT) in the case of random fields. This will be useful, since, for instance, asymptotic spatial homogeneity of order 1 of spatial risk measures associated with VaR (at a confidence level α ( 0 , 1 ) \ { 1 / 2 } ) and induced by a cost field C C is satisfied as soon as C fulfills the CLT and has a constant expectation (see below). For V R d and r > 0 , we introduce V + r = { x R d : dist ( x , V ) r } , where dist stands for the Euclidean distance. Additionally, we denote by V the boundary of V. A Van Hove sequence in R d is a sequence ( V n ) n N of bounded measurable subsets of R d satisfying V n R d , lim n ν ( V n ) = , and lim n ν ( ( V n ) + r ) / ν ( V n ) = 0 for   all r > 0 . The assumption “bounded” does not always appear in the definition of a Van Hove sequence. Let Cov denote the covariance. In the following, we say that a random field { X ( x ) } x R d such that, for all x R d , E [ X ( x ) ] 2 < , satisfies the CLT if
R d | Cov ( X ( 0 ) , X ( x ) ) | ν ( d x ) < ,
σ X = R d Cov ( X ( 0 ) , X ( x ) ) ν ( d x ) 1 2 > 0 ,
and, for any Van Hove sequence ( V n ) n N in R d ,
1 [ ν ( V n ) ] 1 / 2 V n ( X ( x ) E [ X ( x ) ] ) ν ( d x ) d N ( 0 , σ X 2 ) , as   n ,
where N ( μ , σ 2 ) denotes the normal distribution with expectation μ R and variance σ 2 > 0 . In the case of a random field satisfying the CLT, we have the following result.
Theorem 2.
Let { C ( x ) } x R 2 C . Assume moreover that C has a constant expectation (i.e., for all x R 2 , E [ C ( x ) ] = E [ C ( 0 ) ] ) and satisfies the CLT. Then, we have, for all A A c , that
λ L N ( λ A , C ) E [ C ( 0 ) ] d N 0 , σ C 2 ν ( A ) , f o r   λ .
Proof. 
The result is essentially based on part of the proof of Theorem 4 in Koch (2017). We refer the reader to this proof for details and only provide the main ideas here. First, we show (see the third paragraph of the proof of Theorem 4 in Koch 2017) that, for any A A c and any positive non-decreasing sequence ( λ n ) n N R such that lim n λ n = , the sequence ( λ n A ) n N is a Van Hove sequence. Therefore, since C satisfies the CLT and has a constant expectation, we obtain
λ n L N ( λ n A , C ) E [ C ( 0 ) ] d N 0 , σ C 2 ν ( A ) , for   n .
Second, we deduce (see the proof of Theorem 4, after (44), in Koch 2017) that, for all A A c ,
λ L N ( λ A , C ) E [ C ( 0 ) ] d N 0 , σ C 2 ν ( A ) , for   λ .
This concludes the proof.  □
This theorem will be useful in the following since it will allow us to prove asymptotic spatial homogeneity of order respectively 2 , 1 and 1 for spatial risk measures associated with variance, VaR as well as ES and induced by a cost field satisfying the CLT and additional conditions. Moreover, if λ is large enough, it gives an approximation of the distribution of the normalized spatially aggregated loss:
L N ( λ A , C ) N E [ C ( 0 ) ] , σ C 2 λ 2 ν ( A ) ,
where ≈ means “approximately follows”. Such an approximation can be fruitful in practice, e.g., for an insurance company.

2.4. Max-Stable Random Fields

This concise introduction to max-stable fields is partly based on Koch et al. (2018, Section 2.2). Below, “⋁” denotes the supremum when the latter is taken over a countable set. In any dimension d 1 , max-stable random fields are defined as follows.
Definition 7
(Max-stable random field). A real-valued random field Z ( x ) x R d is said to be max-stable if there exist sequences of functions ( a T ( x ) , x R d ) T 1 > 0 and ( b T ( x ) , x R d ) T 1 R such that, for all T 1 ,
t = 1 T Z t ( x ) b T ( x ) a T ( x ) x R d = d { Z ( x ) } x R d ,
where the { Z t ( x ) } x R d , t = 1 , , T , are independent replications of Z.
A max-stable random field is termed simple if it has standard Fréchet margins, i.e., for all x R d , P ( Z ( x ) < z ) = exp 1 / z , z > 0 .
Now, let { T i ˜ ( x ) } x R d , i = 1 , , n , be independent replications of a random field { T ˜ ( x ) } x R d . Let ( c n ( x ) , x R d ) n 1 > 0 and ( d n ( x ) , x R d ) n 1 R be sequences of functions. If there exists a non-degenerate random field { G ( x ) } x R d such that
i = 1 n T i ˜ ( x ) d n ( x ) c n ( x ) x R d d G ( x ) x R d , for   n ,
then G is necessarily max-stable; see, e.g., de Haan (1984). This explains the relevance and significance of max-stable random fields in the modelling of spatial extremes.
Any simple max-stable random field Z can be written (see, e.g., de Haan 1984) as
Z ( x ) x R d = d i = 1 { U i Y i ( x ) } x R d ,
where the ( U i ) i 1 are the points of a Poisson point process on ( 0 , ) with intensity u 2 ν ( d u ) and the Y i , i 1 , are independent replications of a random field { Y ( x ) } x R d such that, for all x R d , E [ Y ( x ) ] = 1 . The field Y is not unique and is called a spectral random field of Z. Conversely, any random field of the form (17) is a simple max-stable random field. Hence, (17) enables the building up of models for max-stable fields. We now present one of the most famous among such models, the Brown–Resnick random field, which is defined in Kabluchko et al. (2009) as a generalization of the stochastic process introduced in Brown and Resnick (1977). We recall that a random field { W ( x ) } x R d is said to have stationary increments if the distribution of the random field { W ( x + x 0 ) W ( x 0 ) } x R d does not depend on x 0 R d . Provided the increments of W have a finite second moment, the variogram of W, γ W , is defined by
γ W ( x ) = Var ( W ( x ) W ( 0 ) ) , x R d ,
where Var denotes the variance. The Brown–Resnick random field is specified as follows.
Definition 8
(Brown–Resnick random field). Let { W ( x ) } x R d be a centred Gaussian random field with stationary increments and with variogram γ W . Let us consider the random field Y defined by
{ Y ( x ) } x R d = exp W ( x ) Var ( W ( x ) ) 2 x R d .
Then the simple max-stable random field defined by (17) with Y is referred to as the Brown–Resnick random field associated with the variogram γ W . In the following, we will also call this field the Brown–Resnick random field built with W.8
The Brown–Resnick field is stationary (Kabluchko et al. 2009, Theorem 2) and its distribution only depends on the variogram (Kabluchko et al. 2009, Proposition 11).
Now, let ( U i , C i ) i 1 be the points of a Poisson point process on ( 0 , ) × R d with intensity function u 2 ν ( d u ) × ν ( d c ) . Independently, let f i , i 1 , be independent replicates of some non-negative random function f on R d satisfying E R d f ( x ) ν ( d x ) = 1 . Then, it is known that the mixed moving maxima (M3) random field
{ Z ( x ) } x R d = i = 1 { U i f i ( x C i ) } x R d
is a stationary and simple max-stable field. The so-called Smith random field introduced by Smith (1990) is a specific case of M3 random field and is defined immediately below.
Definition 9
(Smith random field). Let Z be written as in (18) with f being the density of a d-variate Gaussian random vector with mean 0 and covariance matrix Σ . Then, the field Z is referred to as the Smith random field with covariance matrix Σ .
As the Brown–Resnick and Smith fields are defined using the random fields-based and M3 representations (17) and (18), respectively, it is usual in the spatial extremes literature to distinguish both models, although the Smith field with covariance matrix Σ corresponds to the Brown–Resnick field associated with the variogram γ W ( x ) = x Σ 1 x , x R d , where designates transposition; see, e.g., Huser and Davison (2013).
Finally, we briefly present the extremal coefficient (see, e.g., Schlather and Tawn 2003) which is a well-known measure of spatial dependence for max-stable random fields. Let { Z ( x ) } x R d be a simple max-stable random field. In the case of two locations, the extremal coefficient function θ is defined by
P Z ( x 1 ) u , Z ( x 2 ) u = exp θ ( x 1 , x 2 ) u , x 1 , x 2 R d ,
where u > 0 .

3. Properties of Some Induced Spatial Risk Measures

In this section, we provide sufficient conditions on the cost field such that some induced spatial risk measures satisfy the axioms presented in Definition 6. First, we consider the case of a general cost field before investigating the relevant case of a cost field being a function of a max-stable random field. In the following, for α ( 0 , 1 ) , q α and ϕ denote the quantile at level α and the standard Gaussian density, respectively. We recall that for a random variable X ˜ with distribution function F, its value-at-risk at confidence level α ( 0 , 1 ) is written VaR α ( X ˜ ) = inf { x R : F ( x ) α } . Moreover, provided E [ | X ˜ | ] < , its expected shortfall at confidence level α ( 0 , 1 ) is defined as
ES α X ˜ = 1 1 α α 1 VaR u X ˜ ν ( d u ) .
Typical values for α are 0.95 and 0.99 . It should be noted that in the actuarial literature, ES is sometimes called tail value-at-risk (see, e.g., Definition 2.4.1 in Denuit et al. 2005).
In the following, we mainly consider the spatial risk measures
R 1 ( A , C ) = E [ L N ( A , C ) ] , A A , C C , R 2 ( A , C ) = Var ( L N ( A , C ) ) , A A , C C , R 3 , α ( A , C ) = VaR α ( L N ( A , C ) ) , A A , C C , R 4 , α ( A , C ) = ES α ( L N ( A , C ) ) , A A , C C .
As a classical risk measure, the expectation is not very satisfying since it does not provide any information about variability. Moreover, as will be seen, the associated spatial risk measures do not take into account the spatial dependence of the cost field. An advantage of variance, VaR and ES lies in the fact that their associated spatial risk measures all take into account (at least) part of this spatial dependence. Historically, the variance has been the dominating risk measure in finance, primarily due to the huge influence of the portfolio theory of Markowitz which uses variance as a measure of risk. However, using the variance is only possible when the normalized spatially aggregated loss has a finite second moment. Moreover, since it allocates the same weight to positive and negative deviations from the expectation, variance is a good risk measure only for distributions which are approximately symmetric around the expectation. Currently, VaR is probably the most widely used risk measure in the finance/insurance industry. Nevertheless, it does not provide any information about the severity of losses which occur with a probability lower than 1 α , which is obviously a serious shortcoming. Moreover, VaR is in general not sub-additive and hence not coherent in the sense of Artzner et al. (1999). ES overcomes these two drawbacks of VaR. Pertaining to the first one, it can be seen from the fact that, if a random variable X ˜ has a continuous distribution function, then
ES α X ˜ = E X ˜ | X ˜ > VaR α ( X ˜ ) .
Hence, the Basel Committee on Banking Supervision proposed the use of ES instead of VaR for the internal models-based approach (Basel Committee on Banking Supervision 2012, Section 3.2.1). However, contrary to VaR, ES is not elicitable (Gneiting 2011), implying that backtesting for ES is more difficult than for VaR.

3.1. General Cost Field

The next result provides sufficient conditions on the cost field C such that the induced spatial risk measure R 1 ( · , C ) satisfies the axioms presented in Definition 6.
Theorem 3.
Let { C ( x ) } x R 2 be a measurable random field having a constant expectation and such that, for all x R 2 , E [ | C ( x ) | ] = E [ | C ( 0 ) | ] < . Then, we have, for all A A , that R 1 ( A , C ) = E [ C ( 0 ) ] . Hence, the spatial risk measure induced by C R 1 ( · , C ) satisfies the axioms of spatial invariance under translation and spatial sub-additivity. If, moreover, E [ C ( 0 ) ] 0 , then R 1 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 0 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = E [ C ( 0 ) ] , A A c .
Proof. 
By assumption, the function x E [ | C ( x ) | ] is constant and hence obviously locally integrable. Consequently, as C is measurable, Proposition 1 gives that C has a.s. locally integrable sample paths. Using Fubini’s theorem and the fact that C has a constant expectation, we have, for all A A , that
R 1 ( A , C ) = 1 ν ( A ) A E [ C ( 0 ) ] ν ( d x ) = E [ C ( 0 ) ] .
Thus, for all v R 2 and A A , R 1 ( A + v , C ) = R 1 ( A , C ) . Moreover, for all A 1 , A 2 A ,
R 1 ( A 1 A 2 , C ) = R 1 ( A 1 , C ) = R 1 ( A 2 , C ) = min { R 1 ( A 1 , C ) , R 1 ( A 2 , C ) } .
Finally, for all A A c and λ > 0 , we have R 1 ( λ A , C ) = E [ C ( 0 ) ] . As | E [ C ( 0 ) ] | E [ | C ( 0 ) | ] < , we have | E [ C ( 0 ) ] | ( 0 , ) , which concludes the proof.  □
The next result is a generalization of Theorem 2 in Koch (2017) and will be useful in the following.
Theorem 4.
Let { C ( x ) } x R 2 C and such that, for all x R 2 , E [ C ( x ) ] 2 < . Moreover, assume that, for all A A ,
A A | E C ( x ) C ( y ) | ν ( d x ) ν ( d y ) < .
Then, for all A A and λ > 0 , we have
R 2 ( λ A , C ) = 1 λ 4 [ ν ( A ) ] 2 λ A λ A Cov ( C ( x ) , C ( y ) ) ν ( d x ) ν ( d y ) .
Condition (19) is satisfied for instance in the following cases:
1. 
For any A A ,
sup x A E [ C ( x ) ] 2 < .
2. 
For all x , y R 2 ,
Cov ( C ( x ) , C ( y ) ) = Cov ( C ( 0 ) , C ( x y ) ) ,
and
R 2 | Cov ( C ( 0 ) , C ( x ) ) | ν ( d x ) < .
Proof. 
For all A A , we consider L ( A , C ) = ν ( A ) L N ( A , C ) . Thus, using Fubini’s theorem and (19), we obtain
E [ L ( A , C ) ] 2 = E A C ( x ) ν ( d x ) 2 = E A C ( x ) ν ( d x ) A C ( y ) ν ( d y ) = A A E C ( x ) C ( y ) ν ( d x ) ν ( d y ) .
Moreover, it is clear that
E L ( A , C ) 2 = A A E C ( x ) E C ( y ) ν ( d x ) ν ( d y ) .
The combination of (23) and (24) gives that
E [ L ( A , C ) ] 2 E L ( A , C ) 2 = A A Cov ( C ( x ) , C ( y ) ) ν ( d x ) ν ( d y ) ,
which implies that
R 2 ( A , C ) = 1 [ ν ( A ) ] 2 A A Cov ( C ( x ) , C ( y ) ) ν ( d x ) ν ( d y ) .
The result is obtained by replacing A with λ A .
We now prove the second part of the theorem, concerning (19). Let A A . In the first case, we obtain, using Cauchy–Schwarz inequality and (20),
A A | E C ( x ) C ( y ) | ν ( d x ) ν ( d y ) A A E [ C ( x ) ] 2 1 2 E [ C ( y ) ] 2 1 2 ν ( d x ) ν ( d y ) A A sup x A E [ C ( x ) ] 2 1 2 sup x A E [ C ( y ) ] 2 1 2 ν ( d x ) ν ( d y ) < .
In the second case, it follows from (21) and (22) that
A A | Cov ( C ( x ) , C ( y ) ) | ν ( d x ) ν ( d y ) = A A | Cov ( C ( 0 ) , C ( x y ) ) | ν ( d x ) ν ( d y ) = A A y | Cov ( C ( 0 ) , C ( z ) ) | ν ( d z ) ν ( d y ) A R 2 | Cov ( C ( 0 ) , C ( z ) ) | ν ( d z ) ν ( d y ) = ν ( A ) σ ˜ C < ,
where
σ ˜ C = R 2 | Cov ( C ( 0 ) , C ( x ) ) | ν ( d x ) .
Thus, (19) is obviously satisfied.  □
We recall that for a random field { C ( x ) } x R 2 such that, for all x R 2 , E [ C ( x ) ] 2 < , we note
σ C = R 2 Cov ( C ( 0 ) , C ( x ) ) ν ( d x ) 1 2 .
The next theorem provides the main results of this subsection. In particular, it gives sufficient conditions on the cost field C such that the induced spatial risk measures R 2 ( · , C ) , R 3 , α ( · , C ) and R 4 , α ( · , C ) satisfy the axioms of asymptotic spatial homogeneity of order 2 , 1 and 1 , respectively.
Theorem 5.
Let { C ( x ) } x R 2 C .
1. 
Assume that C is stationary. Then, provided it exists, any spatial risk measure associated with a law-invariant classical risk measure Π and induced by C satisfies the axiom of spatial invariance under translation.
2. 
Assume that C is such that, for all x R 2 ,
E [ C ( x ) ] 2 < ,
and satisfies (21) and (22). Then, we have, for all A A c , that
R 2 ( λ A , C ) = λ σ C 2 λ 2 ν ( A ) + o 1 λ 2 .
Hence, if σ C > 0 , R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
3. 
Assume that C has a constant expectation and satisfies the CLT. Then, we have, for all A A c , that
R 3 , α ( λ A , C ) = λ E [ C ( 0 ) ] + σ C q α λ [ ν ( A ) ] 1 2 + o 1 λ .
Hence, if α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
4. 
Assume that C has a constant expectation, satisfies the CLT and is such that the random variables λ L N ( λ A , C ) E [ C ( 0 ) ] , λ > 0 , are uniformly integrable. Then, we have, for all A A c , that
R 4 , α ( λ A , C ) = λ E [ C ( 0 ) ] + σ C ϕ ( q α ) λ [ ν ( A ) ] 1 2 ( 1 α ) + o 1 λ .
Hence, R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
1. Let A A , v R 2 and Π be a law-invariant classical risk measure. Using the fact that ν ( A + v ) = ν ( A ) and a change of variable, we obtain
R Π ( A + v , C ) = Π 1 ν ( A + v ) A + v C ( x ) ν ( d x ) = Π 1 ν ( A ) A C ( y + v ) ν ( d y ) .
Due to the stationarity of C, we have, for all v R 2 , that { C ( x ) } x R 2 = d { C ( x + v ) } x R 2 , yielding, since Π is law-invariant,
Π 1 ν ( A ) A C ( y + v ) ν ( d y ) = Π 1 ν ( A ) A C ( x ) ν ( d x ) = R Π ( A , C ) .
The combination of (29) and (30) provides the result.
2. As (21) and (22) are satisfied, we know from Theorem 4 that, for all A A c and λ > 0 , R 2 ( λ A , C ) is well-defined. The result follows from an adapted version of the proof of Theorem 3, Point 3, in Koch (2017). We refer the reader to that proof for the technical parts. We only highlight some of the main steps as well as the main differences here.
The first part consists in showing that
lim λ λ 2 ν ( A ) R 2 ( λ A , C ) = σ C 2 .
Using Theorem 4 and (21), it follows that
R 2 ( λ A , C ) = 1 λ 4 [ ν ( A ) ] 2 λ A λ A Cov ( C ( 0 ) , C ( x y ) ) ν ( d x ) ν ( d y ) .
Let A λ = λ A , λ > 0 . Then, we consider the quantity
T λ = 1 λ 2 ν ( A ) A λ A λ k ( x y ) ν ( d x ) ν ( d y ) , λ > 0 ,
where
k ( x ) = Cov ( C ( 0 ) , C ( x ) ) , x R 2 .
The next step consists in showing that
lim λ T λ = σ C 2 .
For this purpose, we proceed similarly as in Koch (2017), proof of Theorem 3, Point 3. The only difference consists in the fact that here k is not necessarily non-negative. Hence, in order to bound | T 1 , λ | and | T 3 , λ | (quantities defined in Koch 2017) from above, k must be replaced with its absolute value in the corresponding integrals. This is where Condition (22) plays a role. Finally, since, for all λ > 0 ,
T λ = λ 2 ν ( A ) R 2 ( λ A , C ) ,
(31) follows from (32).
In a second part, we easily derive (26) from (31). Now, as a compact subset of R 2 , A is bounded, giving that ν ( A ) ( 0 , ) . Since, moreover, σ C 2 ( 0 , ) , σ C 2 / ν ( A ) ( 0 , ) . Hence, the second part of the result follows from (26).
3. Theorem 2 gives that, for all A A c ,
λ L N ( λ A , C ) E [ C ( 0 ) ] d N 0 , σ C 2 ν ( A ) , for   λ .
Hence, the fact that the quantile function of a normal random variable is continuous on ( 0 , 1 ) , Proposition 0.1 in Resnick (1987) and easy computations (see Koch 2017, proof of Theorem 5) yield (27). Since C satisfies the CLT, we have E [ C ( 0 ) ] 2 < and thus E C ( 0 ) < . Additionally, as α 1 / 2 , we have q α 0 . Moreover, as σ C > 0 (because C satisfies the CLT) and ν ( A ) > 0 , we obtain that σ C q α / [ ν ( A ) ] 1 2 0 . Finally, since α { 0 , 1 } , | q α | < . Furthermore, σ C < and ν ( A ) < , giving that | σ C q α / [ ν ( A ) ] 1 2 | < . The result follows by definition.
4. Since C satisfies the CLT, we have, for all x R 2 , E [ C ( x ) ] 2 < , which implies that, for all x R 2 , E | C ( x ) | < . We easily deduce, using Fubini’s theorem, that E [ | L N ( λ A , C ) | ] is finite, and, therefore, that R 4 , α ( λ A , C ) is well-defined for all A A c and λ > 0 . Theorem 2 gives that, for all A A c ,
λ L N ( λ A , C ) E [ C ( 0 ) ] d N 0 , σ C 2 ν ( A ) , for   λ .
Now, ES is known to be continuous with respect to convergence in distribution in the case of uniformly integrable random variables. For details, see, e.g., Wang et al. (2018), Theorem 3.2 and Example 2.2, Point (ii); the authors’ results concern bounded random variables but the mentioned result can be extended to the case of integrable random variables. Hence, it follows from the fact that the random variables λ L N ( λ , C ) E [ C ( 0 ) ] , λ > 0 , are uniformly integrable, and the expression of ES α for the Gaussian distribution, that
lim λ 1 1 α α 1 VaR u ( λ [ L N ( λ A , C ) E [ C ( 0 ) ] ] ) ν ( d u ) = σ C ϕ ( q α ) [ ν ( A ) ] 1 2 ( 1 α ) .
Moreover, we have
1 1 α α 1 VaR u ( λ [ L N ( λ A , C ) E [ C ( 0 ) ] ] ) ν ( d u ) = 1 1 α α 1 λ VaR u ( L N ( λ A , C ) ) E [ C ( 0 ) ] ν ( d u ) = λ R 4 , α ( λ A , C ) E [ C ( 0 ) ] .
Thus, (33) gives, for all A A c ,
λ R 4 , α ( λ A , C ) E [ C ( 0 ) ] = λ σ C ϕ ( q α ) [ ν ( A ) ] 1 2 ( 1 α ) + o ( 1 ) ,
which yields (28). Now, we have E C ( 0 ) < . Moreover, using the fact that, for all α ( 0 , 1 ) , ϕ ( q α ) ( 0 , ) , and arguments stated at the end of the proof of Point 3, we obtain | σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } | ( 0 , ) . Consequently, the result follows by definition.  □
Remark 5.
In order to establish Points 3 and 4, we take advantage of the fact that both VaR and ES are continuous with respect to convergence in distribution under appropriate assumptions. Hence, similar results may hold for other classical risk measures satisfying continuity with respect to convergence in distribution.
Theorem 5 entails the following important result.
Corollary 1.
Let { C ( x ) } x R 2 C . Moreover, assume that C satisfies (21) and the CLT. Then, we have, for all A A c , that
R 2 ( λ A , C ) = λ σ C 2 λ 2 ν ( A ) + o 1 λ 2 .
Hence, R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
Proof. 
Since C satisfies the CLT, it satisfies (22), (25) and σ C > 0 . Thus, the result follows from Theorem 5, Point 2.  □
The next result provides a convenient condition ensuring the uniform integrability required in Theorem 5, Point 4.
Proposition 2.
Let { C ( x ) } x R 2 C . Assume moreover that C has a constant expectation and satisfies the CLT. If C satisfies (21), then the random variables λ L N ( λ A , C ) E [ C ( 0 ) ] , λ > 0 , are uniformly integrable.
Proof. 
Let, for λ > 0 , M λ = λ L N ( λ A , C ) E [ C ( 0 ) ] . Theorem 2 gives that, for all A A c , M λ d M , for   λ , where M N 0 , σ C 2 / ν ( A ) . Therefore, by the continuous mapping theorem, we obtain
M λ 2 d M 2 , for   λ .
Now, it is clear that, for all λ > 0 , Var ( M λ ) = λ 2 R 2 ( λ A , C ) . Hence, it follows from (34) that Var ( M λ ) λ σ C 2 / ν ( A ) , which gives, since for all λ > 0 , E [ M λ ] = 0 , that E M λ 2 λ E M 2 . Additionally, M 2 is non-negative and integrable. Furthermore, the M λ 2 are non-negative and, for all λ > 0 , E M λ 2 = λ 2 R 2 ( λ A , C ) , which is finite according to Theorem 4 as (19) is satisfied. Therefore, the M λ 2 are integrable. Consequently, using (35) and Theorem 3.6 in Billingsley (1999), we know that the random variables M λ 2 , λ > 0 , are uniformly integrable. This directly yields that the random variables M λ , λ > 0 , are uniformly integrable.  □

3.2. Cost Field Being a Function of a Max-Stable Random Field

We now consider a cost field model written as in (14), i.e.,
C ( x ) x R 2 = E ( x ) D Z ( x ) x R 2 ,
where Z is max-stable and the exposure is uniformly equal to unity. The relevance of using max-stable random fields has been previously highlighted.
In the following, all theorems and corollaries assume Z to be simple, although max-stable fields fitted to real data have generalized extreme-value (GEV) univariate marginal distributions with location, scale and shape parameters η R , τ > 0 and ξ R . However, this does not cause any loss of generality. If { Z ( x ) } x R 2 is a max-stable field with such GEV parameters, we can write
Z ( x ) = η + τ ( Z ˜ ( x ) ξ 1 ) / ξ if ξ 0 , η + τ log ( Z ˜ ( x ) ) if ξ = 0 , x R 2 ,
where { Z ˜ ( x ) } x R 2 is simple max-stable. Thus, there exists a function D 1 such that Z ( x ) = D 1 ( Z ˜ ( x ) ) and Model (36) can be written C ( x ) = D ˜ ( Z ˜ ( x ) ) , where Z ˜ is simple max-stable and D ˜ = D D 1 , with “ ” denoting function composition. On ( 0 , ) , for any ξ 0 the transformation z ˜ η + τ ( z ˜ ξ 1 ) / ξ is increasing and the same holds for the transformation z ˜ η + τ log ( z ˜ ) , implying that D 1 is increasing. Most often, the damage function D is also increasing (e.g., the higher the wind speed, temperature or rainfall amount, the higher the cost) and thus the same is true for D ˜ = D D 1 . Consequently, the requirement in Corollaries 3–5 (see below) on the function applied to the simple max-stable field to be non-decreasing and non-constant is generally satisfied in the applications motivating the present work.
For the sake of notational simplicity, in the following, we denote by Z (instead of Z ˜ ) the simple max-stable field and by D (instead of D ˜ ) the quantity D D 1 . Accordingly, the reader should pay attention to the fact that Z models the standardized environmental field (and not the real one) and D consists in the composition of the marginal transformation of Z and the damage function.
We first give sufficient conditions on the function D and the field Z such that the spatial risk measure R 1 ( · , D ( Z ) ) induced by the cost field D ( Z ) satisfies the axioms presented in Definition 6.
Corollary 2.
Let { Z ( x ) } x R 2 be a simple max-stable random field and D a measurable function such that { C ( x ) } x R 2 = { D ( Z ( x ) ) } x R 2 C and E [ | C ( 0 ) | ] < . Then, for all A A , R 1 ( A , C ) = E [ C ( 0 ) ] . Hence, R 1 ( · , C ) satisfies the axioms of spatial invariance under translation and spatial sub-additivity. If, moreover, E [ C ( 0 ) ] 0 , then R 1 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 0 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = E [ C ( 0 ) ] , A A c .
Proof. 
Since Z has identical margins, for all x R 2 , E | C ( x ) | = E | C ( 0 ) | . Therefore, the result directly follows from Theorem 3.  □
The result below gives sufficient conditions on D and Z such that the spatial risk measure R 2 ( · , D ( Z ) ) induced by the cost field D ( Z ) satisfies the axiom of asymptotic spatial homogeneity of order 2 .
Theorem 6.
Let { Z ( x ) } x R 2 be a simple and sample-continuous max-stable random field and D a measurable function such that { C ( x ) } x R 2 = { D ( Z ( x ) ) } x R 2 C and such that there exist p , q > 0 satisfying 2 / p + 1 / q = 1 such that
E | C ( 0 ) | p <
and
R 2 [ 2 θ ( 0 , x ) ] 1 q ν ( d x ) < ,
where θ is the extremal coefficient function of Z. Then, we have
R 2 | Cov ( C ( 0 ) , C ( x ) ) | ν ( d x ) < .
Additionally, assume that C satisfies (21) and σ C > 0 . Then R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
Proof. 
Since Z has identical margins, (38) yields that, for all x R 2 , E | C ( x ) | p < . Thus, using the fact that 2 / p + 1 / q = 1 , Davydov’s inequality (Davydov 1968, (2.2)) gives that
| Cov ( C ( 0 ) , C ( x ) ) | 12 α C ( { 0 } , { x } ) 1 q E | C ( 0 ) | p 1 p E | C ( x ) | p 1 p .
For all x R 2 , since D is measurable, C ( x ) = D ( Z ( x ) ) is F { x } Z -measurable. Hence, F { x } C F { x } Z , which gives by (15) that, for all x R 2 ,
α C { 0 } , { x } α Z { 0 } , { x } .
Now, using (16) and Corollary 2.2 in Dombry and Eyi-Minko (2012), we obtain that, for all x R 2 ,
α Z { 0 } , { x } 2 [ 2 θ ( 0 , x ) ] .
Thus, the combination of (41) and (42) gives that
α C { 0 } , { x } 2 [ 2 θ ( 0 , x ) ] .
Consequently, (40) gives that
| Cov ( C ( 0 ) , C ( x ) ) | 12 2 1 q E | C ( 0 ) | p E | C ( x ) | p 1 p [ 2 θ ( 0 , x ) ] 1 q .
Therefore, using (38) and (39), we obtain
R 2 | Cov ( C ( 0 ) , C ( x ) ) | ν ( d x ) < .
Since p , q > 0 and 2 / p + 1 / q = 1 , we have p > 2 . Consequently, for all x R 2 , E [ C ( x ) ] 2 < . Thus, Theorem 5, Point 2, gives the result.  □
Until the end, the following results provide sufficient conditions on D and Z such that the induced spatial risk measures R 2 ( · , D ( Z ) ) , R 3 , α ( · , D ( Z ) ) and R 4 , α ( · , D ( Z ) ) satisfy the axioms of asymptotic spatial homogeneity of order 2 , 1 and 1 , respectively. In order to establish them, we take advantage of the results in Koch et al. (2018) about the existence of a CLT for functions of stationary max-stable random fields. Let B ( R ) and B ( ( 0 , ) ) be the Borel σ -fields on R and ( 0 , ) , respectively. For h = ( h 1 , h 2 ) Z 2 , we adopt the notation [ h , h + 1 ] = [ h 1 , h 1 + 1 ] × [ h 2 , h 2 + 1 ] . Next theorem considers a general simple, stationary and sample-continuous max-stable random field.
Theorem 7.
Let { Z ( x ) } x R 2 be a simple, stationary and sample-continuous max-stable random field and D be a measurable function from ( ( 0 , ) , B ( ( 0 , ) ) ) to ( R , B ( R ) ) satisfying
E | D ( Z ( 0 ) ) | 2 + δ < ,
for some δ > 0 . Furthermore, assume that, for all h Z 2 ,
E min sup x [ 0 , 1 ] 2 { Y ( x ) } , sup x [ h , h + 1 ] { Y ( x ) } K h b ,
for some K > 0 , b > 2 max 2 , ( 2 + δ ) / δ and where { Y ( x ) } x R 2 is a spectral random field of Z (see (17)). Let { C ( x ) } x R 2 = { D ( Z ( x ) ) } x R 2 . Then, if σ C > 0 :
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
Since Z is sample-continuous, it is measurable. Thus, the function D being measurable from ( ( 0 , ) , B ( ( 0 , ) ) ) to ( R , B ( R ) ) , we obtain that C is measurable. Moreover, it follows from the stationarity of C (due to the stationarity of Z) and Condition (43) that, for all x R 2 , E | C ( x ) | = E | C ( 0 ) | < . Therefore, the function x E [ | C ( x ) | ] is constant and hence obviously locally integrable. Consequently, Proposition 1 gives that C has a.s. locally integrable sample paths. Therefore, C C .
Furthermore, the assumptions enable us to apply Theorem 2 in Koch et al. (2018). The latter yields that the random field C satisfies the CLT. Finally, since C is stationary, it satisfies (21) and has a constant expectation. Hence, Corollary 1 gives the first result. The second result follows from Theorem 5, Point 3. The combination of Proposition 2 and Point 4 in Theorem 5 yields the third result. □
Theorem 7 directly entails the following result.
Corollary 3.
Let Z, D and C be as in Theorem 7 (but without assuming that σ C > 0 ). Moreover, assume that D is non-decreasing and non-constant. Then:
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
Proposition 1 in Koch et al. (2018) gives that σ C > 0 . Therefore, Theorem 7 yields the result.  □
The next results concern the Brown–Resnick and Smith max-stable random fields. The Brown–Resnick model is of high practical interest since, owing to its flexibility, it appears as one of the best (if not the best) models among currently available max-stable models, at least for environmental data; see, e.g., Davison et al. (2012, Section 7.4), in the case of rainfall.
Theorem 8.
Let { Z ( x ) } x R 2 be the Brown–Resnick random field associated with the variogram γ W ( x ) = m x ψ , where m > 0 and ψ ( 0 , 2 ] , or the Smith random field with covariance matrix Σ, and D be as in Theorem 7. Let { C ( x ) } x R 2 = { D ( Z ( x ) ) } x R 2 . Then, if σ C > 0 :
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
We start with the proof in the case of the Brown–Resnick field. As previously mentioned, the Brown–Resnick random field is stationary. Thus, C is stationary and hence satisfies (21) and has a constant expectation. Moreover, we can see from the proof of Theorem 3 in Koch et al. (2018) that Z is sample-continuous. Consequently, the same arguments as in the proof of Theorem 7 yield that C C . Furthermore, Theorem 3 in Koch et al. (2018) gives that C satisfies the CLT. Therefore, Corollary 1 yields the first result. The second result follows from Theorem 5, Point 3. The combination of Proposition 2 and Point 4 in Theorem 5 gives the third result.
The Smith random field is stationary as an instance of M3 random field. Thus, C is stationary and consequently satisfies (21) and has a constant expectation. Moreover, as the Smith field is sample-continuous, the same arguments as in the proof of Theorem 7 yield that C C . Additionally, Theorem 4 in Koch et al. (2018) gives that C satisfies the CLT. Therefore, Corollary 1 yields the first result. The second result follows from Theorem 5, Point 3. The combination of Proposition 2 and Point 4 in Theorem 5 gives the third result.  □
Next corollary easily follows from Theorem 8.
Corollary 4.
Let Z, D and C be as in Theorem 8 (but without assuming that σ C > 0 ). Moreover, assume that D is non-decreasing and non-constant. Then:
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
As explained in the proof of Theorem 8, both such Brown–Resnick fields and the Smith field are stationary and sample-continuous. Furthermore, they are simple max-stable. Thus, Proposition 1 in Koch et al. (2018) gives that σ C > 0 . Hence, Theorem 8 yields the result.  □
Let . denote the Euclidean norm in R 2 . We introduce B 1 = x R 2 : x = 1 , the unit ball of R 2 . For two functions g 1 and g 2 from R 2 to R , the notation g 1 ( h ) = h o ( g 2 ( h ) ) means that lim h sup u B 1 | g 1 ( h u ) / g 2 ( h u ) | = 0 . Moreover, lim h g 1 ( h ) = must be understood as lim h inf u B 1 g 1 ( h u ) = .
Theorem 9.
Let { Z ( x ) } x R 2 be the Brown–Resnick random field built with a random field { W ( x ) } x R 2 which is sample-continuous and whose variogram satisfies
sup x [ 0 , 1 ] 2 { γ W ( h ) γ W ( x + h ) } = h o ( γ W ( h ) ) ,
and
lim h γ W ( h ) ln ( h ) = .
Moreover, let D be as in Theorem 7. Let { C ( x ) } x R 2 = { D ( Z ( x ) ) } x R 2 . Then, if σ C > 0 :
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
The same arguments as in the proof of Theorem 8 show that C satisfies (21) and has a constant expectation. As W is sample-continuous, Proposition 13 in Kabluchko et al. (2009) gives that Z is sample-continuous. Thus, the same arguments as in the proof of Theorem 7 show that C C . Moreover, Remark 3 in Koch et al. (2018) gives that C satisfies the CLT. Hence, Corollary 1 gives the first result. The second result follows from Theorem 5, Point 3. The combination of Proposition 2 and Point 4 in Theorem 5 yields the third result.  □
The following result is a direct consequence of Theorem 9.
Corollary 5.
Let Z, D and C be as in Theorem 9 (but without assuming that σ C > 0 ). Moreover, assume that D is non-decreasing and non-constant. Then:
1. 
R 2 ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 2 with K 1 ( A , C ) = 0 and K 2 ( A , C ) = σ C 2 / ν ( A ) , A A c .
2. 
For all α ( 0 , 1 ) \ { 1 / 2 } , R 3 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C q α / [ ν ( A ) ] 1 2 , A A c .
3. 
For all α ( 0 , 1 ) , R 4 , α ( · , C ) satisfies the axiom of asymptotic spatial homogeneity of order 1 with K 1 ( A , C ) = E [ C ( 0 ) ] and K 2 ( A , C ) = σ C ϕ ( q α ) / { [ ν ( A ) ] 1 2 ( 1 α ) } , A A c .
Proof. 
The random field Z is simple, stationary, sample-continuous (see the proof of Theorem 9) and max-stable. Thus, Proposition 1 in Koch et al. (2018) gives that σ C > 0 . Consequently, Theorem 9 yields the result.  □
We conclude this section by commenting on the damage function D ( z ) = I { z > u } , z > 0 , for u > 0 , which is considered in Koch (2017). This function is measurable from ( ( 0 , ) , B ( ( 0 , ) ) ) to ( R , B ( R ) ) . Moreover, it is bounded and hence obviously satisfies (43) for every random field Z. Additionally, this function is non-decreasing and non-constant. Consequently, the results of Theorem 3, Point 3 and Theorem 5, Point 2 in Koch (2017) concerning the Brown–Resnick random field associated with the variogram γ W ( x ) = m x ψ , where m > 0 and ψ ( 0 , 2 ] , and the Smith random field, are particular cases of Corollary 4.

4. Conclusions

In this paper, we first explore the notions of spatial risk measure and corresponding axioms introduced in Koch (2017) further as well as describe their utility for both actuarial science and practice. Second, in the case of a general cost field, we provide sufficient conditions such that spatial risk measures associated with expectation, variance, VaR as well as ES and induced by this cost field satisfy the axiom of asymptotic spatial homogeneity of order 0, −2, −1 and −1, respectively. Finally, in the case where the cost field is a function of a max-stable random field, we give sufficient conditions on both the function and the max-stable field such that spatial risk measures associated with expectation, variance, VaR as well as ES and induced by the resulting cost field satisfy the axiom of asymptotic spatial homogeneity of order 0, −2, −1 and −1, respectively. Hence, these conditions allow one to know the rate of spatial diversification when the region under study becomes large, which is valuable for the banking/insurance industry. Overall, this paper improves our comprehension of the concept of spatial risk measure as well as of their properties with respect to the space variable and, among others, generalizes several results to be found in Koch (2017).
Ongoing work consists of the study of concrete examples of spatial risk measures involving max-stable fields and relevant damage functions. Inter alia, we apply our theory to winter storm risk over a specific European region. As previously mentioned, max-stable fields have GEV univariate marginal distributions with three parameters. The first step involves jointly fitting the latter and the dependence parameters of different max-stable models (Smith, Brown–Resnick, …) to wind speed maxima using, e.g., composite likelihood methods (see, e.g., Padoan et al. 2010). Model selection has then to be performed employing, for instance, the composite likelihood information criterion. The second step consists in choosing an appropriate damage function and exposure field and leads, in combination with the first one, to the cost field model. If the sufficient conditions mentioned in Section 3.2 are met, then we can draw conclusions about the asymptotic rate of spatial diversification, and less importantly spatial invariance under translation. Simulating from the cost field, we obtain realizations of the normalized spatially aggregated loss on chosen sub-regions, which allow the estimation of the spatial risk measures of interest. This enables one to check, e.g., whether the axiom of spatial sub-additivity is satisfied.
Future work will include the study of spatial risk measures associated with other classical risk measures (e.g., more general distorsion risk measures than VaR or ES and expectile risk measures) and/or induced by cost fields involving other kinds of random fields than max-stable fields. For instance, it would be worthwhile to investigate whether spatial risk measures associated with VaR and ES can still satisfy the axiom of asymptotic spatial homogeneity of order −1 in the case where the cost field does not satisfy the CLT.

Funding

This research was partly funded by the Swiss National Science Foundation grant number 200021_178824.

Acknowledgments

The author would like to thank Anthony C. Davison and Christian Y. Robert for some interesting comments. He also acknowledges Paul Embrechts and Ruodu Wang for providing some references about robustness of risk measures as well as Gennady Samorodnitsky for a fruitful exchange about integrals of random fields. Finally, he is also grateful to the associate editor and two anonymous referees for insightful suggestions.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the writing of the manuscript and in the decision to publish the results.

References

  1. Anderes, Ethan B., and Michael L. Stein. 2011. Local likelihood estimation for nonstationary random fields. Journal of Multivariate Analysis 102: 506–20. [Google Scholar] [CrossRef]
  2. Artzner, Philippe, Freddy Delbaen, Jean-Marc Eber, and David Heath. 1999. Coherent measures of risk. Mathematical Finance 9: 203–28. [Google Scholar] [CrossRef]
  3. Basel Committee on Banking Supervision. 2012. Fundamental Review of the Trading Book. Available online: https://www.bis.org/publ/bcbs219.htm (accessed on 25 April 2019).
  4. Bevere, Lucia, and Lea Mueller. 2014. Natural catastrophes and man-made disasters in 2013: Large losses from floods and hail; Haiyan hits the Philippines. Sigma Swiss Re 1. [Google Scholar]
  5. Billingsley, Patrick. 1999. Convergence of Probability Measures. New York: John Wiley & Sons. [Google Scholar]
  6. Brown, Bruce M., and Sidney I. Resnick. 1977. Extreme values of independent stochastic processes. Journal of Applied Probability 14: 732–39. [Google Scholar] [CrossRef]
  7. Dahlhaus, Rainer. 2012. Locally stationary processes. In Handbook of Statistics. Kidlington: Elsevier, pp. 351–413. [Google Scholar] [CrossRef]
  8. Davison, Anthony C., Simone A. Padoan, and Mathieu Ribatet. 2012. Statistical modeling of spatial extremes. Statistical Science 27: 161–86. [Google Scholar] [CrossRef]
  9. Davydov, Yu A. 1968. Convergence of distributions generated by stationary stochastic processes. Theory of Probability and Its Applications 13: 691–96. [Google Scholar] [CrossRef]
  10. De Haan, Laurens. 1984. A spectral representation for max-stable processes. The Annals of Probability 12: 1194–204. [Google Scholar] [CrossRef]
  11. De Haan, Laurens, and Ana Ferreira. 2006. Extreme Value Theory: An Introduction. New York: Springer. [Google Scholar] [CrossRef]
  12. Denuit, Michel, Jan Dhaene, Marc Goovaerts, and Rob Kaas. 2005. Actuarial Theory for Dependent Risks: Measures, Orders and Models. Chichester: John Wiley & Sons. [Google Scholar]
  13. Dombry, Clément, and Frédéric Eyi-Minko. 2012. Strong mixing properties of max-infinitely divisible random fields. Stochastic Processes and their Applications 122: 3790–811. [Google Scholar] [CrossRef]
  14. Eckley, Idris A., Guy P. Nason, and Robert L. Treloar. 2010. Locally stationary wavelet fields with application to the modelling and analysis of image texture. Journal of the Royal Statistical Society: Series C (Applied Statistics) 59: 595–616. [Google Scholar] [CrossRef]
  15. Gneiting, Tilmann. 2011. Making and evaluating point forecasts. Journal of the American Statistical Association 106: 746–62. [Google Scholar] [CrossRef]
  16. Huser, Raphaël, and Anthony C. Davison. 2013. Composite likelihood estimation for the Brown–Resnick process. Biometrika 100: 511–18. [Google Scholar] [CrossRef]
  17. Kabluchko, Zakhar, Martin Schlather, and Laurens de Haan. 2009. Stationary max-stable fields associated to negative definite functions. The Annals of Probability 37: 2042–65. [Google Scholar] [CrossRef]
  18. Koch, Erwan. 2017. Spatial risk measures and applications to max-stable processes. Extremes 20: 635–70. [Google Scholar] [CrossRef]
  19. Koch, Erwan, Clément Dombry, and Christian Y. Robert. 2018. A central limit theorem for functions of stationary mixing max-stable random fields on R d . Stochastic Processes and Their Applications. [Google Scholar] [CrossRef]
  20. Ombao, Hernando C., Jonathan A. Raz, Rainer von Sachs, and Beth A. Malow. 2001. Automatic statistical analysis of bivariate nonstationary time series. Journal of the American Statistical Association 96: 543–60. [Google Scholar] [CrossRef]
  21. Padoan, Simone A., Mathieu Ribatet, and Scott A. Sisson. 2010. Likelihood-based inference for max-stable processes. Journal of the American Statistical Association 105: 263–77. [Google Scholar] [CrossRef]
  22. Resnick, Sidney I. 1987. Extreme Values, Regular Variation, and Point Processes. New York: Springer. [Google Scholar] [CrossRef]
  23. Rosenblatt, Murray. 1956. A central limit theorem and a strong mixing condition. Proceedings of the National Academy of Sciences of the United States of America 42: 43–47. [Google Scholar] [CrossRef] [PubMed]
  24. Samorodnitsky, Gennady, and Murad S. Taqqu. 1994. Stable Non-Gaussian Random Processes: Stochastic Models with Infinite Variance. London: Chapman and Hall/CRC. [Google Scholar]
  25. Schlather, Martin, and Jonathan A. Tawn. 2003. A dependence measure for multivariate and spatial extreme values: Properties and inference. Biometrika 90: 139–56. [Google Scholar] [CrossRef]
  26. Smith, Richard L. 1990. Max-Stable Processes and Spatial Extremes. Unpublished manuscript. Available online: https://pdfs.semanticscholar.org/033b/fb040b67dbe584a62ae900ad537f5f0eef0b.pdf (accessed on 25 April 2019).
  27. Volkonskii, V. A., and Yu A. Rozanov. 1959. Some limit theorems for random functions. I. Theory of Probability and Its Applications 4: 178–97. [Google Scholar] [CrossRef]
  28. Wang, Ruodu, Yunran Wei, and Gordon Willmot. 2018. Characterization, Robustness and Aggregation of Signed Choquet Integrals. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=2956962 (accessed on 25 April 2019).
1
Throughout the paper, insurance also refers to reinsurance.
2
Throughout the paper, when applied to random fields, the adjective “measurable” means “jointly measurable”.
3
Unless otherwise stated, by a.s., we mean P -a.s.
4
5
Throughout the paper, stationarity refers to strict stationarity.
6
It is out of the scope of this paper to enter into accounting details.
7
Potentially different from those developed in the natural catastrophes industry: e.g., a max-stable model.
8
In the following, when W is sample-continuous, what we refer to as the Brown–Resnick random field built with W is obtained by taking replications of W (see (17)) which are also sample-continuous.

Share and Cite

MDPI and ACS Style

Koch, E. Spatial Risk Measures and Rate of Spatial Diversification. Risks 2019, 7, 52. https://doi.org/10.3390/risks7020052

AMA Style

Koch E. Spatial Risk Measures and Rate of Spatial Diversification. Risks. 2019; 7(2):52. https://doi.org/10.3390/risks7020052

Chicago/Turabian Style

Koch, Erwan. 2019. "Spatial Risk Measures and Rate of Spatial Diversification" Risks 7, no. 2: 52. https://doi.org/10.3390/risks7020052

APA Style

Koch, E. (2019). Spatial Risk Measures and Rate of Spatial Diversification. Risks, 7(2), 52. https://doi.org/10.3390/risks7020052

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop