# Measuring Spatial Dependence among Maxima

Measuring Spatial Dependence among Maxima ... Davis and Mikosch (2007) ... Complex non-parametric structure Vh(u,v) = 2 Z1 0...

Measuring Spatial Dependence among Maxima • P. Naveau Laboratoire des Sciences du Climat et de l’Environnement IPSL, CNRS, France • Guillou, A.; Cooley, D.; Diebolt, J. http://amath.colorado.edu/faculty/naveau/

[email protected]

Outline of the Talk

1. Motivations 2. Maxima distribution 3. Geostatistics 4. Estimation

[email protected]

2

2 1

spatial dependence as

0

20

How to describe the

a function of the

−1

10

y

30

3

40

Spatial Statistics for Extremes

distance between two points?

10

20

30

40

x

Motivations Max

Geostat

Estimation

3

40

Spatial Statistics for Extremes

30

● ●

20

10 ●

for extreme events?

● ●

● ●

● ●

interpolation

spatial

●●

How to perform ●●

● ●

● ●

● ●

0

y

● ●

10

20

30

40

x

Motivations Max

Geostat

Estimation

4

Spatial Statistics for Extremes

A few Approaches for modeling spatial extremes Max-stable processes: Adapting asymptotic results for multivariate extremes Schlather & Tawn (2003), Naveau et al. (2007), de Haan & Pereira (2005) Bayesian or latent models: spatial structure indirectly modeled via the EVT parameters distribution Coles & Tawn (1996), Cooley et al. (2005) Linear filtering: Auto-Regressive spatio-temporal heavy tailed processes, Davis and Mikosch (2007) Gaussian anamorphosis: Transforming the field into a Gaussian one Wackernagel (2003)

Motivations Max

Geostat

Estimation

5

Univariate case for Maxima Convergence of sample maxima Normal density ⇒

⇐ Gumbel density

Uniform density ⇒

⇐ Weibull density

Cauchy density ⇒

⇐ Fr´ echet density

n = 50

Motivations

Max Geostat

n = 100

Estimation

6

Assumptions

Suppose we know the marginal distributions of maxima M (x) with M (x) = the maximum recorded at the location x from a stationary and isotropic field.

Without loss of generality, we first assume that the margins follow an unit Fr´ echet

F (u) = P[M (x) ≤ u] = exp(−1/u)

Motivations

Max Geostat

Estimation

7

A central question

For large n

P [Mn(x) < u, Mn(x + h) < v ] = ??

Motivations

Max Geostat

Estimation

8

Bivariate case for Maxima Asymptotic theory If one assumes that we have unit Fr´ echet margins then

"

lim P

n→∞

#

Mn(x) − an Mn(x + h) − an 6 u, 6 v = exp [−Vh(u, v)] bn bn

where Z 1

w 1−w Vh(u, v) = 2 max , u v 0 



dLh(w)

R1 with Lh(.) a distribution function on [0, 1] such that 0 w dLh(w) = 0.5.

Motivations

Max Geostat

Estimation

9

Bivariate case (M (x), M (x + h)) Complex non-parametric structure

Z 1

w 1−w max Vh(u, v) = 2 , u v 0 Special case u = v 



dLh(w)

Note Vh(u, u) = Vh(1, 1)/u Notations:

θ(h) := Vh(1, 1)

P [M (x) < u, M (x + h) < u] = exp(−θ(h)/u) = F (u)θ(h) because F (u) = exp(−1/u) ⇓

Motivations

Max Geostat

Estimation

10

θ(h) = Extremal coefficient

P [M (x) < u, M (x + h) < u] = F (u)θ(h)

Interpretation

Independence ⇒ θ(h) = 2 M (x) = M (x + h) ⇒ θ(h) = 1 Similar to correlation coefficients for Gaussian but ... No characterization of the full bivariate dependence ⇓

Motivations

Max Geostat

Estimation

11

An important question (1) How to estimate θ(h)?

  R1 w 1−w Vh(u, v) = 2 0 max u , v dLh(w)

12

Geostatistics: Variograms 2 γ(h) = 1 2 E|Z(x + h) − Z(x)|

1.0

● ●

Finite if light tails Capture

0.6

all

spatial

if

{Z(x)}

structure 0.4

Gaussian fields

0.2

semivariance

0.8

0.0

0.2

0.4

0.6

for extremes

0.8

distance

Motivations

Max

Geostat Estimation

13

A Different Variogram

|F (M (x + h)) − F (M (x))| with F (u) = exp(−1/u)

Motivations

Max

Geostat

Estimation

14

A Different Variogram

1 νh = E |F (M (x + h)) − F (M (x))| 2 with F (u) = exp(−1/u)

Defined for light & heavy tails Called a Madogram Nice links with extreme value theory

Motivations

Max

Geostat

Estimation

15

A Different Variogram

1 νh = E |F (M (x + h)) − F (M (x))| 2 Why does it work?

1 1 |a − b| = max(a, b) − (a + b) 2 2

a = F (M (x + h)) and b = F (M (x)) Ea = Eb = 1/2 θ(h) E max(a, b) = EF (max(M (x {z + h), M (x)})) = | 1 + θ(h) max-stable

Motivations

Max

Geostat

Estimation

16

Madogram νh ⇒ Extremal coeff θ(h)

1 + 2νh θ(h) = 1 − 2νh The madogram νh gives the extremal coefficient θ(h)

Motivations

Max

Geostat

Estimation

17

Comparisons with other estimators Gumbel (1960)   

  1  1 α P (X ≤ x, Y ≤ y ) = exp −  +  x 

  ! 1 α 1 α    y 

Four estimators - Pickands’ estimator (1975) - Deheuvels’ estimator (1991) - Hall and Tajvidi’s estimator (2000) - Cap´ era` a et al. (1997) estimator

Motivations

Max

Geostat

Estimation

18

Comparisons with other estimators

α = 0.3

Motivations

α = 0.7

Max

Geostat

Estimation

19

1 −1

10

0

20

y

2

30

3

40

Schlather’s models (2003)

10

20

30

40

x

s

θ(h) = 1 +

Motivations

1 1 − (exp(−h/40) + 1) 2 Max

Geostat

Estimation

20

Madogram νh ⇒ Extremal coeff θ(h) Schlather’s fields Extremal coeff

2.0

0.8

● ●

● ●

1.8

1.6

● ●

1.4

0.4

thetaHat

● ●

● ●

0.2

16

18

20

● ●

1.2

● ●

● ●

● ●

1.0

1

4

6

8

10

12

14

16

18

20

1

4

6

8

distance

0.0

0.6

Motivations

10

12

14

distance

Max

Geostat

Estimation

21

0

1

20 10

y

2

30

3

40

Smith’s models (2003)

10

20

30

40

x

θ(h) = 2Φ

Motivations

Max

p

hT Σ−1h/2

Geostat



Estimation

22

Madogram νh ⇒ Extremal coeff θ(h) Smith’s fields

2.0

Extremal coeff

1.0

● ●

1.8

0.8

● ●

thetaHat

● ● ●

0.4

● ●

1.6

1.4

0.6

● ●

● ●

1.2

0.2

● ●

● ●

1.0

0.0

1

4

6

8

10

12

14

16

18

20

1

4

6

8

distance

Motivations

10

12

14

16

18

20

distance

Max

Geostat

Estimation

23

Building valid Extremal coeff Proposition A Any extremal coefficient function θ(h) is such that 2 − θ(h) is positive semi-definite. Proposition B Any extremal coefficient function θ(h) satisfies the following inequalities θ(h + k) ≤ θ(h)θ(k), θ(h + k)τ ≤ θ(h)τ + θ(k)τ − 1, for all 0 ≤ τ ≤ 1, θ(h + k)τ ≥ θ(h)τ + θ(k)τ − 1, for all τ ≤ 0.

Motivations

Max

Geostat

Estimation

24

An important question (1) How to estimate θ(h) = Vh(1, 1)?

Done!!

(2) How to estimate Vh(u, v)? Note: Because Vh(u, v) = Vh(u/(u + v), v/(u + v))/(u + v) is sufficient to only estimate Vh(λ, 1 − λ) for λ ∈ [0, 1].

  R1 w 1−w Vh(u, v) = 2 0 max u , v dLh(w)

25

1 λ 1−λ νh(λ) = E F (M (x + h)) − F (M (x)) 2

Defined for light & heavy tails Called a λ-Madogram Nice links with extreme value theory νh(0) = νh(1) = 0.25

Motivations

Max

Geostat Estimation

26

λ 1 1−λ νh(λ) = 2 E F (M (x + h)) − F (M (x))

27

A fundamental relationship Vh(λ, 1 − λ) 3 νh(λ) = − c(λ), with c(λ) = 1 + Vh(λ, 1 − λ) 2(1 + λ)(2 − λ) Conversely, Vh(λ, 1 − λ) =

Motivations

Max

c(λ) + νh(λ) 1 − c(λ) − νh(λ)

Geostat Estimation

28

Estimation of Vh(u, v) Suppose that we have T iid years of daily annual maxima fields with unknown margins. λ 1 1−λ (M (x)) ? How to estimate νh(λ) = 2 E F (M (x + h)) − F

A naive estimator T 1 X λ 1−λ ν ˆh(λ) = Fn,T (Mn,t(x + h)) − Gn,T (Mn,t(x)) 2T t=1

with T T 1 X 1 X Fn,T (u) = 1l and Gn,T (u) = 1l T t=1 {Mn,t(x+h)≤u} T t=1 {Mn,t(x)≤u}

But, the conditions Eν ˆh(0) = Eν ˆh(1) = 0.25 are not satisfied

Motivations

Max

Geostat Estimation

29

Estimation of Vh(u, v) λ 1 1−λ How to estimate νh(λ) = 2 E F (M (x + h)) − F (M (x)) ?

A modified estimator T 1 X λ 1−λ ν ˆh(λ) = Fn,T (Mn,t(x + h)) − Gn,T (Mn,t(x)) 2T t=1 T   λ X λ 1 − Fn,T (Mn,t(x + h)) − 2T t=1 T   1−λ X 1−λ − 1 − Gn,T (Mn,t(x)) 2T t=1

1 1 − λ + λ2 + 2 (2 − λ)(1 + λ)

Motivations

Max

Geostat

Estimation

30

Simulations: 300 iid Schalther’s fields

Motivations

Max

Geostat

Estimation

31

Mean Square Error from simulations

300 iid Schalther’s fields

32

Notations for asymptotic results Margins of (X, Y ): unknown continuous margins: F, G Bivariate distribution H and copula: 

H(x, y) = C(F (x), G(y)) and C(u, v) = H F ←(u), G←(v) 

φ(H)(u, v) := H F ←(u), G←(v)





Bivariate empirical process ZT (u, v): √

ZT (u, v) :=

T φ(HT )(u, v) − φ(H)(u, v)

with T 1 X HT (u, v) = 1l{Xt≤u,Yt≤v} T t=1

Motivations

Max

Geostat

Estimation

33

Asymptotic properties

 √  ZT (u, v) = T φ(HT )(u, v) − φ(H)(u, v)

Proposition 1. Let (Xt, Yt)t=1,...,T be a sample of T bivariate random vectors with df H, continuous margins F and G, and with its associated copula C whose partial derivatives are continuous. Then, the process {ZT (u, v), 0 ≤ u, v ≤ 1} converges weakly to the Gaussian process {NC (u, v), 0 ≤ u, v ≤ 1} in `∞([0, 1]2) that is defined as ∂C ∂C NC (u, v) = BC (u, v) − BC (u, 1) (u, v) − BC (1, v) (u, v), ∂u ∂v where BC is a Brownian bridge on [0, 1]2 with covariance function 



E BC (u, v) · BC (u0, v 0) = C(u ∧ u0, v ∧ v 0) − C(u, v) · C(u0, v 0)

Motivations

Max

Geostat

Estimation

34

Convergence of the λ−madogram (X1, Y1), ..., (XT , YT ) T bivariate rv with unknown margins F and G λ  1−λ T  1 X F (Xt) − G (Yt) νbT (λ) := T T 2T t=1

Proposition 2. Under the assumptions of Proposition 1, let J be a function of bounded variation, continuous. Then, we have  T X 

 

1 √ J FT (Xt), GT (Yt) − EJ F (X), G(Y )   T t=1 

d

−→



Z [0,1]2

NC (u, v)dJ(u, v)

λ − y 1−λ|, provides the weak convergence of The special case, J(x, y) := 1 |x 2 the λ−madogram estimator

Motivations

Max

Geostat

Estimation

35

Madogram & EVT • (Z(x), Z(x + h))= precip. measurements at two nearby locations 





Mn(x), Mn(x + h) =



max Zi(x), max Zi(x + h)

i=1,...,n

i=1,...,n

n = recording unit, either hourly, daily or monthly

• Suppose that such bivariate vectors can be computed for a series of years and that these vectors are assumed to be iid in time T T 1 X 1 X Fn,T (u) = 1l{Mn,t(x+h)≤u} and Gn,T (u) = 1l{Mn,t(x)≤u} T t=1 T t=1 T   1 X λ 1−λ νbn,T (h, λ) = Fn,T Mn,t(x + h) − Gn,T (Mn,t(x)) 2T t=1

Motivations

Max

Geostat

Estimation

36

Madogram & EVT (cont’d) Proposition 3. Let (Mn,t(x), Mn,t(x + h)) be a sample of T bivariate vectors with that satisfies the assumptions of Proposition 1 and such that  Mn,t(x)−an Mn,t(x+h)−an , bn bn

converges in distribution to a bivariate EV distri-

bution with an extremal function defined by Vh(., .). Then, we have

 √  1 T νbn,T (h, λ) − E|F λ(M (x + h)) − F 1−λ(M (x))| 2 d

−→

Z [0,1]2

NC (u, v)dJ(u, v)

where n tends to ∞ as T goes to ∞

Motivations

Max

Geostat

Estimation

37

An application: in Bourgogne (Dijon)

Locations (in Lambert coordinates). Pre-processed 30-year maxima of daily precipitation

Motivations

Max

Geostat

Estimation

38

X 1 1 1 − λ + λ2 λ 1−λ (M (xi)) + F (M (xj )) − F 2|Nh| (x ,x )∈N 2 (2 − λ)(1 + λ) i

j

h

    X X λ 1−λ λ 1−λ − 1 − F (M (xi)) − 1−F (M (xi)) 2|Nh| (x ,x )∈N 2|Nh| (x ,x )∈N i

j

i

h

j

h

where Nh is the set of sample pairs lagged by the distance h.

Motivations

Max

Geostat

Estimation

39

Estimated λ−madogram for the field of maxima of daily precipitation over 1970-1999

Motivations

Max

Geostat

Estimation

40

Take-home messages Fields of maxima 6= Gaussian ones Spatial structure defined by the function Vh(u, v) λ−Madogram νh ⇒ dependence function Vh(u, v) We have proposed and study an estimator ν ˆh(λ)

Future research Develop spatial interpolation methods for maxima Derive statistical schemes for downscaling for maxima

Motivations

Max

Geostat

Estimation

41