avalanche crown-depth distributions edward (ned) bair 1, jeff dozier 1, and karl birkeland 2 photo...

25
Avalanche Crown-Depth Distributions Edward (Ned) Bair 1 , Jeff Dozier 1 , and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton, CO Photo courtesy of Mammoth Mountain Ski Patrol 1 Donald Bren School of Environmental Science and Management University of California – Santa Barbara 2 U.S.D.A. Forest Service National Avalanche Center

Upload: abigayle-knight

Post on 20-Jan-2016

214 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Avalanche Crown-Depth DistributionsEdward (Ned) Bair1, Jeff Dozier1, and Karl Birkeland2

Photo courtesy of Center for Snow and Avalanche Studies, Silverton, CO

Photo courtesy of Mammoth Mountain Ski Patrol

1 Donald Bren School of Environmental Science and ManagementUniversity of California – Santa Barbara2 U.S.D.A. Forest Service National Avalanche Center

Page 2: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

We are surrounded by high variability data

Earthquakes

Forest Fires

Stock Markets

Page 3: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

What is a power law?

Linear scale Log scale

Normal distribution Power law distribution

Log scale

Linear scale

Page 4: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Power Laws: More Normal Than Normal

Why? • Strong invariance properties specifically, for crown depths, maximization.

•Power laws are the most parsimonious model for high variability data (Willinger et al 2004)

W. Willinger et al., “More "normal" than normal: scaling distributions and complex systems,” Proceedings of the 2004 Simulation Conference, p. 141. doi: 10.1109/WSC.2004.1371310

Page 5: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Self Organized Criticality (SOC)

• Natural systems spontaneously organize into self-sustaining critical states.

Highly Optimized Tolerance (HOT)

• Systems are robust to common perturbations, but fragile to rare events.

Page 6: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Debate in snow science on power laws

• Do avalanches follow power law or lognormal distributions?

• What is the generating mechanism?

• Is universal?

Page 7: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Why is this important?•May answer why some avalanches are much deeper than others.

•Paths with low are stubborn!

•A universal exponent would mean all paths have the same proportion of large to small avalanches.

Page 8: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Mammoth Mountain Ski Patrol (1968-2008)3,106 crowns > 1/3 meter)

Westwide Avalanche Network (1968-1995)61,261 crowns > 1/3 meter from 29 avalanche areas

Data

Page 9: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

dMethods

•Maximum likelihood

•3 tests of significance: KS, 2,rank-sum

•Ranked by probability of fit

Page 10: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Distributions (excluding log normal)

Truncated power law:

Page 11: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Mammoth Patrol Crowns (N=3,106, =3.3-3.5)

Page 12: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,
Page 13: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Lex parsimoniae:All else being equal, the

simplest explanation is best.

• The generalized extreme value, and its special case, the Fréchet, provide the best fit.

• The simplest generating mechanism is a collection of maxima.

• The scaling exponent ( varies significantly with path and area.

Page 14: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Smoothing

We had to smooth the WAN data:

1) We assume WAN data are rounded uniformly ± 0.5 feet (i.e. a 2ft crown is between 1.5 ft and 2.5 ft).

2) We add a uniform random number on the unit interval (0 to 1).

3) We then subtract 0.5 ft, and convert to meters.

Page 15: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

WAN Alyeska Crowns (N=4,562, =3.5-3.7)

Page 16: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

WAN Snowbird Crowns (N=3,704, =3.6-3.8)

Page 17: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

WAN Squaw Valley Crowns (N=3,926, =3.9-4.1)

Page 18: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

WAN Alpine Meadows Crowns (N=3,435, =4.2-4.4)

Page 19: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,
Page 20: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Post Control Accident on a Stubborn Path

Mammoth Mountain, CA Path: Climax 1:52pm 4/17/06.

Page 21: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,
Page 22: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

4/17/2006 post-control video

Page 23: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Fréchet exponent () for selected Mammoth paths

Climax, =2.5, is the most stubborn of the 34 paths plotted here

Page 24: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

Fréchet MDA requires the parent distribution to be scaling

Page 25: Avalanche Crown-Depth Distributions Edward (Ned) Bair 1, Jeff Dozier 1, and Karl Birkeland 2 Photo courtesy of Center for Snow and Avalanche Studies, Silverton,

AcknowledgmentsWalter Rosenthal

National Science Foundation

Mammoth Mountain Ski Patrol

USFS and Know Williams