-
Notifications
You must be signed in to change notification settings - Fork 0
/
bushfire_conversation.tex
236 lines (185 loc) · 13.6 KB
/
bushfire_conversation.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
% Options for packages loaded elsewhere
\PassOptionsToPackage{unicode}{hyperref}
\PassOptionsToPackage{hyphens}{url}
%
\documentclass[
11pt,
a4paper,
]{article}
\usepackage{lmodern}
\usepackage{amssymb,amsmath}
\usepackage{ifxetex,ifluatex}
\ifnum 0\ifxetex 1\fi\ifluatex 1\fi=0 % if pdftex
\usepackage[T1]{fontenc}
\usepackage[utf8]{inputenc}
\usepackage{textcomp} % provide euro and other symbols
\else % if luatex or xetex
\usepackage{unicode-math}
\defaultfontfeatures{Scale=MatchLowercase}
\defaultfontfeatures[\rmfamily]{Ligatures=TeX,Scale=1}
\fi
% Use upquote if available, for straight quotes in verbatim environments
\IfFileExists{upquote.sty}{\usepackage{upquote}}{}
\IfFileExists{microtype.sty}{% use microtype if available
\usepackage[]{microtype}
\UseMicrotypeSet[protrusion]{basicmath} % disable protrusion for tt fonts
}{}
\makeatletter
\@ifundefined{KOMAClassName}{% if non-KOMA class
\IfFileExists{parskip.sty}{%
\usepackage{parskip}
}{% else
\setlength{\parindent}{0pt}
\setlength{\parskip}{6pt plus 2pt minus 1pt}}
}{% if KOMA class
\KOMAoptions{parskip=half}}
\makeatother
\usepackage{xcolor}
\IfFileExists{xurl.sty}{\usepackage{xurl}}{} % add URL line breaks if available
\IfFileExists{bookmark.sty}{\usepackage{bookmark}}{\usepackage{hyperref}}
\hypersetup{
pdftitle={Yes, 2019-2020 Victorian bushfires were most likely caused by lightning and there's we can learn from open data},
pdfauthor={Di Cook, Weihao Li, Emily Dodwell},
hidelinks,
pdfcreator={LaTeX via pandoc}}
\urlstyle{same} % disable monospaced font for URLs
\usepackage[margin=1in]{geometry}
\usepackage{longtable,booktabs}
% Correct order of tables after \paragraph or \subparagraph
\usepackage{etoolbox}
\makeatletter
\patchcmd\longtable{\par}{\if@noskipsec\mbox{}\fi\par}{}{}
\makeatother
% Allow footnotes in longtable head/foot
\IfFileExists{footnotehyper.sty}{\usepackage{footnotehyper}}{\usepackage{footnote}}
\makesavenoteenv{longtable}
\usepackage{graphicx,grffile}
\makeatletter
\def\maxwidth{\ifdim\Gin@nat@width>\linewidth\linewidth\else\Gin@nat@width\fi}
\def\maxheight{\ifdim\Gin@nat@height>\textheight\textheight\else\Gin@nat@height\fi}
\makeatother
% Scale images if necessary, so that they will not overflow the page
% margins by default, and it is still possible to overwrite the defaults
% using explicit options in \includegraphics[width, height, ...]{}
\setkeys{Gin}{width=\maxwidth,height=\maxheight,keepaspectratio}
% Set default figure placement to htbp
\makeatletter
\def\fps@figure{htbp}
\makeatother
\setlength{\emergencystretch}{3em} % prevent overfull lines
\providecommand{\tightlist}{%
\setlength{\itemsep}{0pt}\setlength{\parskip}{0pt}}
\setcounter{secnumdepth}{5}
%% Any special functions or other packages can be loaded here.
\usepackage[no-weekday]{eukdate}
\usepackage{sourceserifpro}
\usepackage[scaled=0.86]{DejaVuSansMono}
\usepackage{float,bm,setspace}
\setstretch{1.2}
%% CAPTIONS
\RequirePackage{caption}
\DeclareCaptionStyle{italic}[justification=centering]
{labelfont={bf},textfont={it},labelsep=colon}
\captionsetup[figure]{style=italic,format=hang,singlelinecheck=true}
\captionsetup[table]{style=italic,format=hang,singlelinecheck=true}
%% GRAPHICS
\RequirePackage{graphicx}
\setcounter{topnumber}{2}
\setcounter{bottomnumber}{2}
\setcounter{totalnumber}{4}
\renewcommand{\topfraction}{0.85}
\renewcommand{\bottomfraction}{0.85}
\renewcommand{\textfraction}{0.15}
\renewcommand{\floatpagefraction}{0.8}
%% BIBLIOGRAPHY
\makeatletter
\@ifpackageloaded{biblatex}{}{\usepackage[style=authoryear-comp, backend=biber, natbib=true]{biblatex}}
\makeatother
\ExecuteBibliographyOptions{bibencoding=utf8,minnames=1,maxnames=3, maxbibnames=99,dashed=false,terseinits=true,giveninits=true,uniquename=false,uniquelist=false,doi=false, isbn=false,url=true,sortcites=false}
\DeclareFieldFormat{url}{\texttt{\url{#1}}}
\DeclareFieldFormat[article]{pages}{#1}
\DeclareFieldFormat[inproceedings]{pages}{\lowercase{pp.}#1}
\DeclareFieldFormat[incollection]{pages}{\lowercase{pp.}#1}
\DeclareFieldFormat[article]{volume}{\mkbibbold{#1}}
\DeclareFieldFormat[article]{number}{\mkbibparens{#1}}
\DeclareFieldFormat[article]{title}{\MakeCapital{#1}}
\DeclareFieldFormat[inproceedings]{title}{#1}
\DeclareFieldFormat{shorthandwidth}{#1}
% No dot before number of articles
\usepackage{xpatch}
\xpatchbibmacro{volume+number+eid}{\setunit*{\adddot}}{}{}{}
% Remove In: for an article.
\renewbibmacro{in:}{%
\ifentrytype{article}{}{%
\printtext{\bibstring{in}\intitlepunct}}}
\makeatletter
\DeclareDelimFormat[cbx@textcite]{nameyeardelim}{\addspace}
\makeatother
\renewcommand*{\finalnamedelim}{%
%\ifnumgreater{\value{liststop}}{2}{\finalandcomma}{}% there really should be no funny Oxford comma business here
\addspace\&\space}
\usepackage[style=authoryear-comp,]{biblatex}
\addbibresource{references.bib}
\title{Yes, 2019-2020 Victorian bushfires were most likely caused by lightning and there's we can learn from open data}
\author{Di Cook, Weihao Li, Emily Dodwell}
\date{}
\begin{document}
\maketitle
\hypertarget{intro}{%
\section{Introduction}\label{intro}}
The 2019-2020 Australia bushfires, compared with other major bushfires in history, had a more devastating impact on the environment and properties: 3094 houses were destroyed and over 17M hectares of land burned. Fortunately, fewer lives were lost.
A Jan 11, 2020 ABC article ``The truth about Australia's fires --- arsonists aren't responsible for many this season'' reports about Victoria:
\begin{quote}
\emph{The Country Fire Authority (CFA) said the majority of fires were not arson-related. `Most of the fires have been caused by lightning', said Brett Mitchell, the CFA incident controller in Bairnsdale, in East Gippsland. `Our intelligence suggests there are no deliberate lightings that we are aware of.' Victoria Police had no arson figures available for this bushfire season, but said in the 12 months to September 2019, a dozen people had been arrested for causing bushfires.}
\end{quote}
These words were spoken in response to a fierce war of hashtags on Twitter: ``\#ClimateEmergency'' vs ``\#ArsonEmergency''. In the ABC article, most other states quoted numbers for arson and lightning, but Victoria's response was more vague. This prompted us to dig into open data resources to learn as much as we could about the causes of the bushfires in Victoria.
Historical data collected manually on Victorian fires from 2000-2019 attributes 41\% of fires to lightning, 17\% to arson, 34\% from accidents and 7\% from burning off. When this data is combined with information on weather, distance to roads, camp sites, CFA stations, and fuel load, as we might expect, we learn that proximity to human activity, and winds, tends to be the distinguishing factors between the causes.
Although the cause of the 2019-2020 season fires is not officially known, our model predictions, using satellite hotspots data, weather, distance to roads, camp sites, CFA stations, and fuel load, suggest that 82\% of the fires in Victoria can be attributed to lightning. Only 4\% due to arson, 14\% through accidents, and 1\% from burning off. A benefit of the hot spot data is that it provides eyes in locations that are hard for people to get to. Our methodology utilises open data and open source software that could be used for the rest of Australia, too.
\hypertarget{data}{%
\section{Data sources and collation}\label{data}}
Historical fire origins data from the Victorian Department of Environment, Land, Water and Planning (DELWP), spanning 2000-2019 summers, suggests that most fires have been caused by lightning. The number of fires has been increasing in recent years, and most of these have been attributed to accidental causes. (This latter observation is possibly disturbing, and should be explained by the relevant authorities.) Arson ignitions appear to be decreasing.
\begin{figure}[H]
{\centering \includegraphics[width=1\linewidth]{bushfire_conversation_files/figure-latex/fireorigins-1}
}
\caption{Causes of fires from 2000-2019. Lightning is most common cause. The number of fires is increasing, and this is mostly due to accidents.}\label{fig:fireorigins}
\end{figure}
To track bushfires in Australia remotely with high temporal and spatial resolution, we used hotspot data taken from the Himawari-8 satellite. This data was filtered using the firepower information to keep only hotspots that were most likely fires. We clustered the hotspots in time and space to estimate the ignition point and time of fires, and to label individual fires. This data was supplemented with data from other sources: temperature, rainfall, wind, solar exposure, fuel load, as well as distance to camp sites, roads and CFA stations. Rainfall was aggregated at different time intervals to provide cumulative moisture for locations.
A model was trained to predict one of four causes (lightning, accident, arson, burning off) on the historical data, which was supplemented with the same weather, fuel layer, and distances data as the hot spot data. The performance of the model was very good: 75\% overall accuracy, 90\% accurate on lightning, 78\% for accidents, and 54\% for arson, which was mostly confused with accident, as would make sense. The most important contributors to distinguishing between lightning and arson (or accident) ignition were distances to CFA stations, roads and camp sites, and average wind speed over the past 12-24 months. Smaller distances were most likely arson or accident, as might be expected.
\hypertarget{models}{%
\section{Model predictions}\label{models}}
The model predicts the majority of fires in the 2019-2020 bushfires to be due to lightning. Most of the fires were located in densely vegetated areas of the state, inaccessible by road, which is similar to the historical locations.
\begin{figure}[H]
{\centering \includegraphics[width=1\linewidth]{bushfire_conversation_files/figure-latex/maps-1}
}
\caption{Spatial distribution of causes of fires from 2000-2019, and predictions for 2019-2020 season.}\label{fig:maps}
\end{figure}
When the fire ignitions are examined by month during the season, it can be seen that lightning ignition has occurred throughout the summer. All of the fires in February were predicted to be due to lightning. Accident and arson were commonly predicted causes in March, and early in the season. Reassuringly, ignition due to burning off was predicted primarily in October, prior to the fire restrictions.
\begin{figure}[H]
{\centering \includegraphics[width=1\linewidth]{bushfire_conversation_files/figure-latex/monthlymaps-1}
}
\caption{Spatio-temporal distribution of cause predictions for 2019-2020 season. Reassuringly, fires due to burning off primarily occured in October, prior to fire restrictions. February fires were all predicted to be due to lightning.}\label{fig:monthlymaps}
\end{figure}
\hypertarget{validation-of-predictions}{%
\section*{Validation of predictions}\label{validation-of-predictions}}
\addcontentsline{toc}{section}{Validation of predictions}
In order to validate predictions, we have searched media releases and descriptions. Information is sparse to non-existent. Wikipedia's page \href{https://en.wikipedia.org/wiki/2019-20_Australian_bushfire_season}{2019--20 Australian bushfire season} reports
\begin{quote}
\emph{On 21 November 2019, lightning strikes ignited a series of fires in East Gippsland}
\end{quote}
which is consistent with our predictions.
Close to Melbourne there was disturbing news of a bushfire in a high density koala habitat on French Island on Jan 18. \href{http://www.mpnews.com.au/2020/01/21/all-clear-after-island-fire-fright/}{Mornington Peninsula News} reports:
\begin{quote}
\emph{Dry lightning is being blamed for the blaze which was thought to have started on the previous Wednesday and smouldered for three days before an east-wind change brought it surging to life.}
\end{quote}
Contrary to this, our model predicts this to be arson or accident with 72\% confidence. The other data about this ignition, used in our model prediction, was that it is non-forest, 3m to a road, quite a warm day (21-37\(^o\)C), moderate winds and no rain. If we zoom into a map of the area to the first detected hotspot it can be seen to be right near a road, making plausible evidence for this to actually have been a human-caused fire.
\includegraphics{french_island.png}
\hypertarget{summary}{%
\section{Summary}\label{summary}}
In summary, the ability to detect fire starts from remotely sensed hot spots makes it possible to more rapidly predict fire cause, particularly when combined with weather, distance to human activity, and vegetation information. Our analysis using open data supports that lightning is the most common cause of bushfires in Victoria.
\hypertarget{acknowledgements}{%
\section{Acknowledgements}\label{acknowledgements}}
This analysis is based on thesis research conducted by Weihao Li, a Monash University Honours student, under the supervision of Professor Di Cook from the Department of Econometrics and Business Statistics at the Monash Business School, and Emily Dodwell, former Principal Inventive Scientist at AT\&T Labs Research in New York City. The full analysis is available at \url{https://github.com/TengMCing/bushfire-conversation}. The Australian Centre of Excellence for Mathematical and Statistical Frontiers supported Emily's travel to Australia to start this project.
The code used for the analysis can be found in the github repository \href{https://github.com/TengMCing/thesis}{github.com/TengMCing/thesis}.
You can explore the historical fire data, predictions for 2019-2020 fires, and a fire risk map for Victoria using the shiny app at \href{https://ebsmonash.shinyapps.io/VICfire/}{ebsmonash.shinyapps.io/VICfire/}
\printbibliography
\end{document}