![Page 1: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/1.jpg)
Wildfires and Elk and the US Forest Service
David R. Brillinger
Department of Statistics
University of California, Berkeley
www.stat.berkeley.edu/~brill
AIF/CIF Salt Lake City, October 9, 2014
![Page 2: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/2.jpg)
The USFS, USGS, and US Census can carry out the big
projects, ones that smaller institutions cannot. They also
have the resources to make the data conveniently
available..
About 15 years ago Ager, Brillinger, Preisler (ABP)
began produced a series of papers on animal movement
and wildfire risk.
Other USFS researchers involved were: Burgan, Finney,
McRoberts, Benoit, Wisdom, Kie, and Fujioka .
Canadian researchers involved too, a succession of
meetings at BIRS. Martell, McAlpine, Taylor, Woolford,
Wooton.
![Page 3: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/3.jpg)
“Communicating Science in an Uncertain World”
Scientists communicate by talks, papers, teaching and nowadays
even creating dotcoms.
Today focus on 4 AGP papers.
An ‘uncertain world”. Indeed and statisticians make it even more
uncertain. In part by formalizing and estimating aleatoric
(statistical) and epistemic (systematic)uncertainty.
Four papers of our papers will now discussed with graphics and
uncertainty highlighted, two concern elk trajectories/paths and
two concern wildfire risk,
![Page 4: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/4.jpg)
Overall approach involves exploratory data analysis followed by
formal statistical modelling, fitting and model assessment.
The trajectory papers proceed by working with potential functions
(i.e. formalizing attraction and repulsion of moving objects).
The risk papers make use of spline functions to bring in smooth
relationships.
![Page 5: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/5.jpg)
Elk trajectories.
Elk, deer, and cows,. live in the Starkey Reserve while hikers, riders, bikers and hunters visit.
Starkey is a major experimental field area set up by the USFS
- the largest research enclosure ever built to study usage and interactions.
Trajectories are defined by position estimates and times.
Source of data, (x(tj), y(tj), tj), j = 1,…,J
http://www.fs.fed.us/pnw/starkey/introduction.shtml
![Page 6: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/6.jpg)
Motivation: differential equation based on Newton’s
second law of motion
.
r = (x,y)
![Page 7: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/7.jpg)
Earliest paper
Brillinger, D. R., Preisler, H. K., Ager, A. A. and Kie, J. G. "The use of potential
functions in modelling animal movement". Data Analysis from Statistical
Foundations (Ed. A. K. M. E. Saleh). Nova Science, Huntington (2001).
LORAN-C for location; 1 < Δt < 4hr; restrict to Δt < 1.5hr
An SDE
A discrete approximation
LHS velocity μ estimated by R-gam
![Page 8: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/8.jpg)
Data
![Page 9: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/9.jpg)
Locations
![Page 10: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/10.jpg)
Simulation of trajectories used to estimate potential function, H.
Potential estimate
![Page 11: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/11.jpg)
But does a potential function exist?
∂Hx / ∂y = ∂Hy / ∂x ?
t-statistic
Jackknife employed, dropping 5 elk out of 50 ten
times
![Page 12: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/12.jpg)
Proportions exceeding 2.262
day 036
night .026
nominal ,050
![Page 13: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/13.jpg)
Brillinger, D. R., Preisler, H. K . and Wisdom, M. J. Modelling particles moving in a potential field with pairwise interactions and an application . Brazilian J. Prob. Statist. 25, 421-436 (2011).
radio collared Δt ≈ 5 min
Most recent paper
![Page 14: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/14.jpg)
Hierarchical
clustering
![Page 15: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/15.jpg)
398 wrt 395 V(x,y) cubic in x, y
![Page 16: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/16.jpg)
![Page 17: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/17.jpg)
Hunters and elk, Suzette Puente’s work
![Page 18: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/18.jpg)
bagplots
![Page 19: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/19.jpg)
Wildfire data.
Probabilistic risk assessment
insurance approach helpful
The data are fire locations (xj, y j, t j) , j = 1, …, J
voxel (x,y,t) N(dx,dy,dt) = 1 or 0 discrete Nx,y,t
ηx,y,t = g1(x,y) + g2(d) + g3(e) linear predictor
Pr{N = 1 | η} = exp{η}/(1 + exp{η}) logistic regression R-gam
![Page 20: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/20.jpg)
Earliest wildfire paper Brillinger, D.R., Preisler, H. K. and Benoit, J. W. "Risk assessment: a forest fire example”. Lecture Notes in Statistics 40, 177-196. IMS (2003)
Federal Lands in Oregon
1989 – 1996, 15786 fires > 0.1 acre
Fire locations
http://www.fs.usda.gov/rds/archive/Product/RDS-2013-
0009.
![Page 21: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/21.jpg)
spatial effect g3
![Page 22: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/22.jpg)
seasonal and elevation effects g2 , g3
![Page 23: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/23.jpg)
assessment of logistic link
![Page 24: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/24.jpg)
Umatilla Forest result
![Page 25: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/25.jpg)
Pr{ one or more fires
}
![Page 26: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/26.jpg)
Most recent paper Brillinger, D. R., Preisler, H. K., and Benoit, J. "Probabilistic risk assessment for wildfires.” Environmetrics 17, 623-633. (2006)
![Page 27: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/27.jpg)
![Page 28: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/28.jpg)
Random effects (empirical Bayes)
![Page 29: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/29.jpg)
![Page 30: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/30.jpg)
Model assessment
![Page 31: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/31.jpg)
Summary.
The USFS has advanced forest science by collecting data in
novel ecological experiments and making it publicall yavailable
and then encouraging its scientists collaborate with outside
researchers.
There has been an emphasis on communication. I communicate as
an academic:, consulting, papers, talks, teaching, thesis
supervision
There has been an emphasis on assessing uncertainty.
I take this opportunity to thank the many foresters particularly
USFS people who made this research possible and Abdel for
publishing much of it. Thank you all.
![Page 32: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/32.jpg)
![Page 33: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/33.jpg)
D. R. Brillinger, B. S. Autrey and M. D. Cattaneo, "Probabilistic Risk Modelling at the Wildland Urban Interface: the 2003 Cedar Fire". Environmetrics
The 2003 Cedar Fire.
Southern California near San Diego. Largest fire in California history. Large data set created for analytic use by San Diego community.
![Page 34: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/34.jpg)
![Page 35: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/35.jpg)
![Page 36: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/36.jpg)
![Page 37: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/37.jpg)
![Page 38: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/38.jpg)
Brillinger, D. R. and Finney, M. An exploratory data analysis of the temperature
fluctuations in a spreading fire Environmetrics DOI:10.1002/env.2279 (2014)
![Page 39: PowerPoint Presentationbrill/AIFCIFtalk18.pdfTitle: PowerPoint Presentation Author: david Created Date: 4/13/2015 2:32:44 PM](https://reader036.vdocument.in/reader036/viewer/2022071218/604eb7cb31217933ab2bc696/html5/thumbnails/39.jpg)