2.1. Spatial Risk Measures and Corresponding Axioms
First, we describe the setting required for a proper definition of spatial risk measures. Let
be the set of all compact subsets of
with a positive Lebesgue measure and
be the set of all convex elements of
. Denote by
the set of all real-valued and measurable
2 random fields on
having almost surely (a.s.)
3 locally integrable sample paths. Let
be the family of all possible distributions of random fields belonging to
. 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
. In the following,
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
be the set of all real-valued random variables defined on
. A risk measure typically will be some function
. 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
,
only depends on the distribution of
.
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 and , the normalized spatially aggregated loss is defined bywhere the random field belongs to and has distribution P. The quantity
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
, which is the loss per surface unit and can be interpreted, in a discrete setting
4 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 is measurable, and are well-defined random variables. Moreover, they are a.s. finite as A is compact and 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 and be a measurable random field. If the functionis locally integrable, then Q has a.s. locally integrable sample paths. Proof. Let
A be a compact subset of
. First, since
Q is measurable,
is a well-defined random variable. By Fubini’s theorem, we have
which necessarily implies that
Since this is true for all A being a compact subset of , 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 that assigns a real number to any region and distribution :where Π
is a classical and law-invariant risk measure and 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 , the quantity 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 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
, let
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 , we define the following axioms for the spatial risk measure associated with Π and induced by P, :
- 1.
Spatial invariance under translation:
for all and , where denotes the region A translated by the vector .
- 2.
Spatial sub-additivity:
for all .
- 3.
Asymptotic spatial homogeneity of order :
for all ,where is the area obtained by applying to A a homothety with center and ratio , and , are functions depending on P.
It is also reasonable to introduce the axiom of spatial anti-monotonicity: for all
,
. 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
(e.g., stationarity
5 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
and
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 which is much less risky than a region , it is very unlikely that the company reduces its risk by covering . 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
be such a sub-region (a subset of
) and
be the set of all compact subsets of
with a positive Lebesgue measure. The axiom of spatial invariance under translation becomes: for all
and
such that
,
; spatial sub-additivity is now written: for all
.
Of course, the fact that the axioms of Definition 3 are satisfied depends on both the classical risk measure and the cost field . 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 and be a measurable random field having a.s. locally integrable sample paths. Moreover, let A be a compact subset of with positive Lebesgue measure. Then the distribution ofonly 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
. For a fixed
, we denote by
the corresponding realization of
H on
and by
the realization of
H at location
. By definition, we have, for almost every
, that
Now, let
be a probability space different from the probability space
. Let
be a random vector defined on
and following the uniform distribution on
A, with density
. From (
3), it directly follows that, for almost every
,
Let us denote by
the expectation with respect to the probability measure
. We have
giving, using (
4),
Now, let
be independent replications of
(which are independent of the random field
H). The strong law of large numbers gives that, for almost every
,
Therefore, using Fubini’s theorem, we deduce that, for
-almost every
,
Now, we choose
such that the (non-random) sequence
satisfies (
6). We obtain
Equation (
7) says that the distribution of
is determined by the finite-dimensional distributions at the points belonging to the set
. 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 Let be a random field with distribution P. Although a particular realization of obviously depends on (through its corresponding realization), we know from Theorem 1 that its distribution is entirely characterized by A and P. This explains our notation instead of in Definition 1. More precisely, let be random fields having the same distribution P. Then, and have the same finite-dimensional distributions, which implies that .
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 that assigns a real number to any region and random field :where Π
is a classical risk measure. For a given classical and law-invariant risk measure
and a given region
, the value of the spatial risk measure of Definition 5 is completely determined by the distribution of
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
,
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 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 , we define the following axioms for the spatial risk measure associated with Π and induced by C, :
- 1.
Spatial invariance under translation:
for all and , where denotes the region A translated by the vector .
- 2.
Spatial sub-additivity:
for all .
- 3.
Asymptotic spatial homogeneity of order :
for all ,where is the area obtained by applying to A a homothety with center and ratio , and , 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
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
where
is the total loss,
N denotes the number of insurance policies and, for
is the claim related to the
i-th policy. The
’s are generally assumed to be independent but not necessarily identically distributed. In
, each location
corresponds to a specific insurance policy and thus each
is equivalent to a
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.,
, where
is a finite set of locations in
(e.g., part of a lattice in
). It is worth mentioning that the ideas of this paper can easily be applied to such a framework.
Even if dependence between the
,
in (
10) was allowed, considering
(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
denote either the claims reserves, revenues or any relevant related quantity
6 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
for a given hazard and potentially aims at covering also a region
disjoint of
. We assume that Ins properly hedges its risk on
, i.e.,
by positive homogeneity. Using again the same property,
Combined with
this yields
Hence, by translation invariance,
It follows from (
11) that
which gives
The combination of (
12) and (
13) yields that
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
, the risk is even better hedged on
. Exactly the same reasoning holds for
.
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 , . 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
,
. Assume that it is satisfied with
(e.g., we will see that for
being VaR or ES,
typically equals 1). It follows from Definition 6, Point 3, that
Since
, the dominant term as
is
. Assume that
and
. This is true under the conditions of
Section 3 for VaR and ES: we have
, which is positive as the cost field can be assumed to be non-negative and not a.s. equal to 0; regarding
, this is always true for ES and, provided that the confidence level
is greater than
, also for VaR. Consequently, for
large enough, the total risk of the company,
, is a decreasing function of
as soon as the revenue per surface unit (or claims reserves, …) satisfies
. Under the conditions of
Section 3, for VaR and ES,
for all
, 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
corresponds to the second highest power with respect to
. Provided that
and
(which is true for VaR and ES under the conditions of
Section 3), the corresponding term,
, 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
. For
large, the values of
,
and
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,
and
.
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 , , would become: for all ,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
where
is the exposure field,
D a damage function and
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
. 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 model
7 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
be a real-valued random field. For
a closed subset, we denote by
the
-field generated by the random variables
. Let
be disjoint closed subsets. The
-mixing coefficient (introduced by
Rosenblatt 1956) between the
-fields
and
is defined by
The
-mixing coefficient or absolute regularity coefficient (attributed to Kolmogorov in
Volkonskii and Rozanov 1959) between the
-fields
and
is given by
where the supremum is taken over all partitions
and
of
with the
’s in
and the
’s in
. These coefficients satisfy the useful inequality
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
of spatial risk measures associated with VaR (at a confidence level
) and induced by a cost field
is satisfied as soon as
C fulfills the CLT and has a constant expectation (see below). For
and
, we introduce
, where
stands for the Euclidean distance. Additionally, we denote by
the boundary of
V. A Van Hove sequence in
is a sequence
of bounded measurable subsets of
satisfying
,
, and
. 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
such that, for all
,
, satisfies the CLT if
and, for any Van Hove sequence
in
,
where
denotes the normal distribution with expectation
and variance
. In the case of a random field satisfying the CLT, we have the following result.
Theorem 2. Let . Assume moreover that C has a constant expectation (i.e., for all , ) and satisfies the CLT. Then, we have, for all , that 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
and any positive non-decreasing sequence
such that
, the sequence
is a Van Hove sequence. Therefore, since
C satisfies the CLT and has a constant expectation, we obtain
Second, we deduce (see the proof of Theorem 4, after (44), in
Koch 2017) that, for all
,
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
,
and
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:
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
, max-stable random fields are defined as follows.
Definition 7 (Max-stable random field)
. A real-valued random field is said to be max-stable if there exist sequences of functions and such that, for all ,where the are independent replications of Z. A max-stable random field is termed simple if it has standard Fréchet margins, i.e., for all , .
Now, let
be independent replications of a random field
. Let
and
be sequences of functions. If there exists a non-degenerate random field
such that
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
where the
are the points of a Poisson point process on
with intensity
and the
, are independent replications of a random field
such that, for all
,
. 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
is said to have stationary increments if the distribution of the random field
does not depend on
. Provided the increments of
W have a finite second moment, the variogram of
W,
, is defined by
where Var denotes the variance. The Brown–Resnick random field is specified as follows.
Definition 8 (Brown–Resnick random field)
. Let be a centred Gaussian random field with stationary increments and with variogram . Let us consider the random field Y defined byThen the simple max-stable random field defined by (
17)
with Y is referred to as the Brown–Resnick random field associated with the variogram . In the following, we will also call this field the Brown–Resnick random field built with W.8 Now, let
be the points of a Poisson point process on
with intensity function
. Independently, let
, be independent replicates of some non-negative random function
f on
satisfying
. Then, it is known that the mixed moving maxima (M3) random field
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 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
, 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
be a simple max-stable random field. In the case of two locations, the extremal coefficient function
is defined by
where
.