A satellite-based burned area dataset for the northern boreal region from 1982 to 2020

SPECIAL ISSUE |RESEARCH PAPER h ps://doi.o g/10.1071/WF22102 A sa elli e-based bu ned a ea da ase o he no he n bo eal egion om 1982 o 2020 José-And és Mo eno-Ruiz A , José-Ra ael Ga cía-Láza o A , Manuel A belo B,* and Ped o A. He nández-Leal B ABSTRACT Backg ound. Fi es in he bo eal o es occu wi h na u al equencies and pa e ns. Bu ned a ea (BA) is an essen ial a iable in assessing he impac o clima e change in bo eal egions. Aims. Spa ial wild i e occu ence da a since he 1950s a e a ailable o No h Ame ica. Howe e , he e a e no eliable da a o Eu asia, mainly o Sibe ia, du ing he 1980s and 1990s. Me hods. A Bayesian- ne wo k algo i hm was applied o he Long-Te m Da a Reco d (LTDR) Ve sion 5 o gene a e a BA Da aSe (BA-LTDR-DS) o he Bo eal egion om 1982 o 2020, alida ed using o icial e e ence da a and compa ed wi h he MODIS MCD64A1 p oduc . Key esul s. A high co ela ion (>93%) wi h all he e e ence BA da ase s was ound. BA-LTDR-DS da a g ouped by decades es ima ed a linea inc ease in BA o 4.47 million ha/decade. This end p o ides e idence o how global wa ming a ec s i e ac i i y in hese bo eal o es s. Conclusions. BA-LTDR-DS cons i u es a unique da a sou ce o he p e-MODIS e a, and becomes a eliable sou ce when o he p oduc s wi h highe spa ial/spec al esolu ion a e no a ailable. Implica ions. The BA-LTDR-DS da ase cons i u es he longes ime se ies de eloped o he bo eal egion a his spa ial esolu ion. BA-LTDR-DS could be used as inpu in global clima e models, helping imp o e wild i e p edic ion capabili ies and unde s and he in e ac ions be ween i e, clima e and ege a ion dynamics. Keywo ds: AVHRR, Bayesian ne wo k algo i hm, bo eal o es , bu ned a ea mapping, Eu asia, LTDR, MODIS, No h Ame ica, emo e sensing, Sibe ia, ime se ies analysis. In oduc ion The bo eal o es is he mos ex ensi e e es ial biome; i occupies ~14% o land on Ea h in a ci cumpola bel su ounding he suba c ic egions o he no he n hemisphe e. Two- hi ds o hese o es s a e in Eu asia (Scandina ia and Russia), he emaining hi d in No h Ame ica (Canada and Alaska). Bo eal o es s play a c i ical ole in egula ing clima e and he global ca bon cycle in he Ea h–a mosphe e sys em (Chapin e al. 2000; Kasischke e al. 2005). Bo eal o es s ha e been conside ed o yea s as a ca bon sink (Jobbágy and Jackson 2000; Ciais e al. 2010; Pan e al. 2011). In ac , mo e ca bon is s o ed in he bo eal o es egions han in any o he egion o he plane , possibly up o wice as much ca bon as is s o ed in opical o es s (B adshaw and Wa ken in 2015). Howe e , ecen s udies show ha bo eal o es ca bon sinks could e en be becoming a ne sou ce o emissions (Bonan 2008; Ku z e al. 2008; B adshaw and Wa ken in 2015; Po ie e al. 2019; Eckdahl e al. 2022). The apid wa ming ha he bo eal and A c ic egions ha e expe ienced o e he pas 30 yea s is modi ying he dynamics o na u al dis u bances ha we e his o ically domina ed by i e (Walsh 2014). This anomalous inc ease in empe a u e al e s he humidi y o he uels and, he e o e, modi ies he se e i y, he na u al egime o bo eal o es i es and he bu ned a ea, which leads, in u n, o possible eedback e ec s on clima e change (Goldamme and Fu yae 1996; Flannigan e al. 2005; Balshi e al. 2009a, 2009b; Tchebako a e al. 2009; Geo giadi e al. 2010; de G oo e al. 2013; Kelly e al. 2013; Co ield e al. 2019). The moni o ing o hese changes and he analysis o u u e scena ios a e o i al impo ance o implemen Fo ull lis o au ho a ilia ions and decla a ions see end o pape *Co espondence o: Manuel A belo Depa amen o de Física, Uni e sidad de La Laguna, 38200 San C is óbal de La Laguna, Spain. Email: [email p o ec ed] Recei ed: 22 June 2022 Accep ed: 15 Ap il 2023 Published: 4 May 2023 Ci e his: Mo eno-Ruiz J-A e al. (2023) In e na ional Jou nal o Wildland Fi e 32(6), 854–871. doi:10.1071/WF22102 © 2023 The Au ho (s) (o hei employe (s)). Published by CSIRO Publishing on behal o IAWF. This is an open access a icle dis ibu ed unde he C ea i e Commons A ibu ion- NonComme cial-NoDe i a i es 4.0 In e na ional License (CC BY-NC-ND) OPEN ACCESS managemen policies on clima e change ha p o ec bo h he bo eal o es and he ca bon i s o es (Bonan e al. 1992; Kasischke e al. 1995; Fuchs e al. 2009; Shuman e al. 2011; Loboda e al. 2012; K ylo e al. 2014; Ponoma e e al. 2016). Sa elli e emo e sensing has become an e ec i e ech- nique o he iden i ica ion and spa io- empo al cha ac e - isa ion o o es i es and bu ned a eas owing o i s co e age on bo h a global and egional scale (And eae 1991; Cooke e al. 1996; Dwye e al. 2000; Duncan 2003; Chu and Guo 2015). Cu en ly, di e en pla o ms p o ide sa elli e images ha a e used in algo i hms o map bu ned a ea (BA) (Vi cha 2011; Mouillo e al. 2014; Campagnolo e al. 2016; Chen e al. 2016a, 2016b). The Mode a e- Resolu ion Imaging Spec o adiome e (MODIS) on boa d NASA’s Te a and Aqua sa elli es and he Ad anced Ve y High Resolu ion Radiome e (AVHRR) on boa d he se ies o Na ional Oceanic and A mosphe ic Adminis a ion (NOAA) sa elli es a e examples o such senso s. They p o ide daily image y and highe -le el land and a mosphe e p oduc s o global mapping. Fo eliabili y and pe iod co e ed, he mos equen ly used BA p oduc s a e MODIS Collec ion 6 (MCD64A1) (Giglio e al. 2018; Bosche i e al. 2019) om he yea 2000 onwa ds, he Global Fi e Emissions Da abase e sion 4 (GFED-4) BA a ailable om mid-1995 h ough o he p esen (Giglio e al. 2013), and he mos ecen MODIS p oduc , Fi e_cci 5.1, a ailable o he yea s 2001–2020 (Chu ieco e al. 2018). As Chu ieco e al. (2008) a gue, ‘Longe ime se ies da a a e equi ed o acqui e a be e unde s anding o i e egimes, and hei mu ual ela ionships wi h global wa ming.’ In spi e o he ac ha cu en sa elli e emo e sensing sys ems (and hei de i ed i e p oduc s) ha e enhanced empo al, spa ial and spec al esolu ions, he a ailabili y o well-buil geospa ial ime se ies is sca ce o Eu asia, especially o Sibe ia (de G oo e al. 2013; Chen e al. 2016a; Ebe le e al. 2016; Ponoma e e al. 2016). In addi ion, measu es om sa elli e da a p esen s ong disc ep- ancies in BA es ima ions wi h ega d o epo ed da a in he o icial eco ds (Soja e al. 2004; Sukhinin e al. 2004; Vi cha 2011; Kuka skaya 2013; Chen e al. 2016b), unlike No h Ame ica, which has been well s udied (Kasischke and F ench 1995; Al-Saadi e al. 2008; Chu ieco e al. 2008; Soja e al. 2009; Kasischke e al. 2011; Mo eno Ruiz e al. 2012; Loboda e al. 2013; Mo eno-Ruiz e al. 2014a, 2014b, 2019). In his pape , we p esen he Bu ned A ea Long-Te m Da a Reco d Da aSe (BA-LTDR-DS), a unique long- e m BA p od- uc o he bo eal o es in he Clima e Modelling G id (CMG) esolu ion o ou decades (1982–2020). The CMG o ma , in a la i ude/longi ude geog aphic p ojec ion wi h a esolu- ion o 0.05°, allows he BA-LTDR-DS da ase o be used as inpu in global clima e models, helping o imp o e wild i e p edic ion capabili ies and unde s and he in e ac ions be ween i e, clima e and ege a ion dynamics in he no h- e n bo eal egion. Ma e ials and me hods S udy egion The s udy egion is geog aphically delimi ed by he pa allels 60°N and 72.5°N, di ided in u n in o wo sub- egions, No h Ame ica and Eu asia (Fig. 1). The i s sub- egion includes Alaska and he no he n pa o he Canadian bo eal egion and is bounded a he uppe igh co ne o he map a 72.5°N, 168.5°W and he bo om le co ne a 60°N 43.5°W. The No h Ame ican sub- egion con ains app oxima ely one- hi d o he en i e bo eal egion o his con inen based on he map desc ibed by B and (2009). The second egion s e ches om Scandina ia o he Paci ic coas s o Sibe ia and is limi ed a he uppe igh co ne a 72.5°N 5°E and he bo om le co ne a 60°N 180°E. E e g een coni e ous o es s (pine and sp uce) p edomina e in he wo sub- egions. Howe e , deciduous o es s, mainly bi ch and la ch, can also be ound, wi h a spa ial dis ibu ion in luenced by pos - i e dynamics (Roge s e al. 2015). The o es s o No h Ame ica end o ha e mo e black sp uce, whi e sp uce and pine species wi h b anches lowe o he g ound, hinne ba k and se o inous cones ha open a e being bu ned by i e. The Eu asian o es s ha e mo e i e- esis an species wi h hick ba k, we e needles and ewe low b anches. In addi- ion, be ween 25 and 30% o he landscape o he bo eal o es egion is pea land (o ganic soils) (Go ham 1991; Wiede e al. 2006; Beaulne e al. 2021; Nelson e al. 2021). The p esence o di e en species be ween No h Ame ican and Eu asian o es s ma ks a no able di e ence in i e egimes be ween he wo egions (de G oo e al. 2013). Fi e egimes can di e in he same biome, and he bo eal o es s o No h Ame ica and Eu asia a e an example o his (Haas e al. 2022). Fi es in No h Ame ica end o be la ge and mo e in ense, wi h highe uel consump ion (c own i es); in con as , i es in Eu asia end o be less in ense, wi h lowe uel consump ion (su ace i es) (Woos e and Zhang 2004; Wi h 2005; de G oo e al. 2013; Si no and Mokho 2018). Re e ence da a The only da abases ha include eliable in o ma ion on i es de ec ed in bo eal o es s o mo e han ou decades a e hose c ea ed and main ained by US and Canadian o es y agencies. Fo Eu asia, howe e , he e is no eliable BA e e ence se co e ing he s udy pe iod unde conside a ion (1980–2020), wi h only MODIS senso -de i ed BA p oduc s a ailable since 2000. The Alaska Fi e Se ice (AFS), a Fo Wainw igh , AK, USA, main ains a de ailed eco d o all de ec ed i e e en s since 1940 (h ps:// i e.ak.blm.go /). This da abase, in addi ion o p o iding he pe ime e o each BA, includes addi ional i e- ela ed in o ma ion such as he managemen o ice, i e name, geog aphical coo dina es, es ima ed a ea, cause and ele an commen s. The pe ime e s a e always delinea ed om he bes a ailable da a sou ce, which can www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 855 include ae ial and high spa ial esolu ion sa elli e image y (Landsa ype), as well as opog aphic maps. AFS ecognises di e ences in he scale and accu acy o he pe ime e s depending on he pe iod. Fi es la ge han 400 ha a e included o i es be o e 1987, hose la ge han 40 ha om 1987 o 1989, whe eas om 1990 onwa ds, all i es wi h a BA la ge han 4 ha a e conside ed. The Canadian Na ional Fi e Da abase (CNFDB) is com- piled and main ained by he Canadian Fo es Se ice (h ps:// cw is.c s.n can.gc.ca/ha/n db). CNFDB eco ds i e da a o all sizes since 1959 ha include i e loca ion and pe ime e da a supplied by Canadian e i o ial i e managemen agencies. The in o ma ion con ained in he CNFDB may no be comple e o e o - ee owing o he di e en mapping echniques used. In addi ion, he comple eness and quali y o he da a may a y be ween agencies and be ween yea s. The quali y o his da abase and i s use ulness ha e been demons a ed in se e al publica ions (Ami o e al. 2001; S ocks e al. 2003; Pa isien e al. 2006; Bu on e al. 2008; Hanes e al. 2019). Recen ly, he Canada Cen e o Mapping and Ea h Obse a ion and he Canadian Fo es Se ice de eloped a new da abase called he Na ional Bu ned A ea Composi e (NBAC) (Hall e al. 2020; Skakun e al. 2021). NBAC imp o es he BA de e mina ion o he CNFDB da abase using an au o- ma ic me hod based on ho spo s and he No malized Di e ence Vege a ion Index (NDVI) applied o high spa ial esolu ion (less han 30 m) sa elli e image y (Hall e al. 2020). NBAC is a ailable o i es eco ded om 1986 onwa ds (Skakun e al. 2022). The MODIS p oduc selec ed was he MCD64A1 Collec ion 6 (Giglio e al. 2018), based on da a om he Te a and Aqua sa elli es and dis ibu ed by he Land P ocesses Dis ibu ed Ac i e A chi e Cen e (LP DAAC). MCD64A1 is he NASA o icial BA p oduc . I is a global mon hly g idded p oduc o 500-m spa ial esolu ion. The algo i hm used o de ec bu ned pixels conside s a mosphe ically co ec ed su ace e lec ances o he sho wa e in a ed Bands 5 and 7 h ough a no malised ege a ion index in conjunc ion wi h ac i e i e da a a 1-km esolu ion. The p oduc con ains he es ima ed bu n da e, unbu ned, o e en unmapped a eas i he e we e no da a o es ablish bu ned/unbu ned s a us (Giglio e al. 2018). MCD64A1 C6 is cu en ly he p oduc wi h he high- es eliabili y compa ed wi h he o he BA p oduc s (Padilla e al. 2015; Mo eno-Ruiz e al. 2020), and i has he lowes commission and omission e o s in he bo eal o es egion (Bosche i e al. 2019). P e-p ocessing o he Long-Te m Da a Reco d (LTDR) The LTDR unded by he Clima e Da a Reco d P og am o he NOAA Na ional Clima ic Da a Cen e is a consis en long- e m da ase a a spa ial esolu ion o 0.05° (~5 km) based on daily da a om he AVHRR onboa d he NOAA sa elli es and daily da a acqui ed by MODIS onboa d NASA’s Te a and Aqua sa elli es (Pedel y e al. 2007). The daily global LTDR e sion 5 (1981–2021) used in he p esen s udy was downloaded om h ps://l d .nascom.nasa.go /cgi-bin/l d /l d Page.cgi. 130°W 50°W 40°W 30°W 20°W 10°W 0°E 10°E 20°E 30°E 40°E 50°E 140°W 150°W 160°W 170°W 180°W 170°E 160°E 150°E 140°E 130°E 120°E 110°E 100°E 90°E 80°E 70°E 60°E 120°W 110°W 100°W 90°W 80°W 70°W 60°W Fig. 1. The s udy egion ( ed pe ime e ) co e s he no he n bo eal egion. I is di ided in o wo sub- egions: No h Ame ica (72.5°N 168.5°W; 60°N 43.5°W) and Eu asia (72.5°N 5°E; 60°N 180°E). Bo eal o es (g een) is di e en ia ed om all he o he land co e s (b own), ice (whi e) and wa e bodies (blue). J-A Mo eno-Ruiz e al. In e na ional Jou nal o Wildland Fi e 856 The o iginal iles, in hie a chical da a o ma , co e he globe a a 0.05° esolu ion CMG wi h 7200 × 3600 cells. The o igi- nal iles we e ans o med in o a bina y sequen ial o ma (BSQ). Con e sion o physical alues (su ace e lec ance and b igh ness empe a u e) conside ed he quali y assessmen (QA) ields o he Daily Su ace Re lec ance p oduc (AVH09C1) o emo e possible snow-co e ed pixels and o il e o he p esence o clouds using he CLAVR-1 (Clouds om AVHRR-Phase I) algo i hm (S owe e al. 1999). Missing and in alid b igh ness empe a u e alues ound o he NOAA-16 and 18 sa elli es om 2000 o 2008 we e di ec ly eplaced by he equi alen alues om he MOD09CMG p od- uc o ha pe iod (h ps://ladsweb.modaps.eosdis.nasa.go / a chi e/allDa a/6/MOD09CMG). In addi ion, he signi ican ly decaying o bi o he NOAA-19 sa elli e om 2018 (Julien and Sob ino 2021; Giglio and Roy 2022) made i necessa y o use he MOD09CMG p oduc o his pe iod as well. Finally, o elimina e esidual clouds and cloud shadows ha could in e e e wi h he disc imina ion o bu ned pixels, he maximum b igh ness empe a u e (BT_CH3: 3.55–3.93 µm) c i e ion cons uc ed 10-day composi es om he LTDR-BSQ iles (Ba bosa e al. 1998). This composi ing c i e ion has p o ed o be e ec i e in disce ning bu ned om unbu ned a eas (Chu ieco e al. 2005). Tempe a u e alues abo e 350 K we e conside ed as e oneous alues. Nex , wo ege- a ion indices de i ed om he o iginal bands we e calcu- la ed, which we e use ul o BA disc imina ion in he bo eal egions: Global En i onmen al Moni o ing Index (GEMI) (Pin y and Ve s ae e 1992) and Bu ned Bo eal Fo es Index (BBFI) (Mo eno Ruiz e al. 2012). Table 1 desc ibes he band con igu a ions o each 10-day composi e ile in he BSQ o ma wi h loa ing da a ype. The bu ned a ea de ec ion algo i hm The me hodology de eloped by Mo eno Ruiz e al. (2012), based on a Bayesian ne wo k algo i hm (BA-LTDR), was applied o he en i e s udy a ea o ob ain he annual BA maps and hei co esponding empo al dis ibu ion (Mo eno Ruiz e al. 2012). This me hodology has been p e iously applied and success ully alida ed in di e en bo eal egions and o di e en pe iods (Núñez-Casillas e al. 2013; Mo eno- Ruiz e al. 2014a; Ga cía-Láza o e al. 2018). In he cu en wo k, he s udy egion was ex ended o he no he n bo eal egion (abo e 60°N) using a single algo i hm. This app oach allowed he cohe ence in he wo sub- egions (No h Ame ica and Eu asia) and he wo s udy pe iods (MODIS and p e-MODIS e as) o be assessed join ly, such ha he es ima es ob ained o one sub- egion and pe iod can be ex apola ed o o he bo eal sub- egions and pe iods. The di e en s eps ha cons i u e his me hodology a e sum- ma ised in Fig. 2 and a e desc ibed in g ea e dep h in p e ious s udies by he same au ho s (Mo eno-Ruiz e al. 2012, 2014b; Núñez-Casillas e al. 2013; Ga cía-Láza o e al. 2018; Guindos-Rojas e al. 2018). The algo i hm calcula es 12 s a is ical a iables based on he su ace e lec ance bands ρ 1 and ρ 2 , he b igh ness empe a u e T, and he BBFI and GEMI indices, o he 10-day composi e o po en- ial i e da es be o e and a e he i e o he yea o he i e e en , he yea be o e and he yea a e . In he no he n bo eal o es , ege a ion akes se e al yea s o eco e , e en mo e han a decade, which is why i was decided o also analyse he yea be o e and a e he i e, wi h a du a ion o 2 mon hs o he p e- i e and pos - i e pe iods o each yea , s a ing om he hypo he ical igni ion da e de e mined by he alue o he maximum o he BBFI (Mo eno Ruiz e al. 2012). The Bayesian ne wo k classi ie calcula ed he no - malised p obabili y o he unbu ned and bu ned classes using a aining se based on he pe ime e s o he BA ha was la ge han 1000 ha in he NE Sibe ia egion in 2010. These pe ime e s we e gene a ed om 53 pai s o Landsa - TM images ob ained om he Uni ed S a es Geological Su ey (USGS) conside ing p e- and pos - i e in o ma ion a 30-m spa ial esolu ion (Ga cía-Láza o e al. 2018). To imp o e he esul ing BA p obabili y maps, he spa ial cohe ence was analysed using a il e ing p ocess based on cellula au oma a heo y (Moja adi e al. 2004; Espinola e al. 2015). Finally, he BA-LTDR algo i hm de eloped was applied o he s udy egion in o de o gene a e a da ase o BA annual maps o he no he n bo eal o es Table 1. Bands con igu a ion o 10-days composi e iles. Band name Desc ip ion Equa ion ρ 1 Su ace e lec ance o ed channel SREFL_CH1 (0.5–0.7 µm) ρ 2 Su ace e lec ance o nea -in a ed channel SREFL_CH2 (0.7–1.0 µm) T Top o a mosphe e b igh ness empe a u e (K) BT_CH3 (3.55–3.93 µm) GEMI Global En i onmen al Moni o ing Index n n× (1 0.25 × ) 0.125 1 1 1 n=2 × ( ) + 1.5 × + 0.5 × + + 0.5 22122 1 2 1 BBFI Bu ned Bo eal Fo es Index +T1 2 2 QA Quali y assessmen ield www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 857 (abo e 60°N) o he 1982–2020 pe iod a 0.05° (~5 km) esolu ion. We e e o his da ase as he BA-LTDR-DS. Technical alida ion To assess he accu acy o he BA-LTDR-DS da ase , we spli he ime se ies o bo h sub- egions (Eu asia and No h Ame ica), in o wo – p e-MODIS and pos -MODIS. Since 2000 (MODIS e a), accu acy assessmen o he Eu asian sub- egion was accomplished using only he MCD64A1 C6 BA p oduc . Al hough MODIS da a a e no he mos sui able o quan i a i e alida ion o he BA-LTDR-DS p oduc , as hey a e a ec ed by e o s due o hei spa ial esolu ion o 500 m, hei use may be conside ed acco ding o he p o o- col o he Commi ee on Ea h Obse a ion Sa elli es (CEOS). CEOS ecommends ha an assessmen can be made o sys ema ic quali y con ol o a p oduc by s a is i- cal compa ison wi h independen ly ob ained BA da ase s o be e spa ial esolu ion when no o he o icial e e ence se exis s (Mo ise e e al. 2006; Bosche i e al. 2009). Be o e 2000 (p e-MODIS e a), accu acy assessmen s could no be e alua ed in his way because o he BA p oduc s wi h a highe spa ial esolu ion o his s udy sub- egion (Eu asia) we e no a ailable. Fo una ely, his is no he case o he No h Ame ican sub- egion, o which o icial e e ence da a a e a ailable o compa e wi h he BA-LTDR-DS da ase , as desc ibed in Re e ence da a sec ion. The AFS da abase and he CNFDB ha e al eady been used success ully as e e ence da a o assess he accu acy o sa elli e-de i ed BA p oduc s (Chu ieco e al. 2008; Chang and Song 2009; Giglio e al. 2009; Núñez-Casillas e al. 2013; Mo eno-Ruiz e al. 2019, 2020). The new e e ence da abase NBAC, a ailable o he pe iod 1986–2020, was used as he e e ence se o Canada. The i s 4 missing yea s (1982–1985) we e comple ed wi h da a om CNFDB. NBAC signi ican ly imp o es he CNFDB BA polygons by including small i es and some i es in emo e loca ions no p e iously conside ed in CNFDB, and by emo ing unbu ned islands and wa e bodies wi hin hose polygons (Hall e al. 2020; Skakun e al. 2021, 2022). The polygons o all he i es egis e ed in he AFS and CNFDB + NBAC be ween 1982 and 2020 we e used o p oduce he annual ec o laye s o g ound- u h e i ica- ion. Nex , hese ec o laye s we e ep ojec ed o a geo- g aphic p ojec ion wi h a pixel size o 0.005° (~500 m) o gene a e annual g ound- u h maps. To de e mine how he pixel was assigned a bu ned/non-bu ned alue, he me hod o maximum a ea wi hin he pixel was used (A none e al. 2016). Howe e , o compu e spa ial and empo al accu acy, he annual BA maps om he BA-LTDR-DS and he MCD64A1 C6 da ase s we e clipped o he No h Ame ican (72.5°N, 168.5°W, 60°N, 141°W) and Eu asian (72.5°N, 5°E, 60°N, 180°E) bo eal sub- egions. All e e ence da a maps we e esized o a geog aphic p ojec ion wi h a pixel size o 5 × 5 km by pixel agg ega ion (an agg ega ed pixel ep e- sen s he pe cen age o BA a he subpixel le el). Daily LTDR-BSQ (No he n bo eal egion) MOD09CMG 10-day composi es (ρ, T, GEMI, BBFI) p e- and pos - i e No h Ame ica Eu asia AFS + (CNFDB + NBAC) (1982–2020) MCD64A1 C6 (2001–2020) MCD64A1 C6 (2001–2020) Accu acy assessmen Bayesian ne wo k algo i hm Bu ned p obabili y maps Spa ial cohe ence analysis 160°W 70°N 65°N 60°N 150°W Daily global LTDR e sion 5 (1981–2021) Fig. 2. Flowcha o he p ocess o ob ain annual maps o he BA-LTDR-DS (Bu ned A ea Long-Te m Da a-Reco d Da ase ) in he no he n bo eal egion o 1982–2020 and o assess hei accu acy agains e e ence da a ( Mo eno Ruiz e al. 2012; Núñez- Casillas e al. 2013; Mo eno-Ruiz e al. 2014b; Ga cía-Láza o e al. 2018; Guindos-Rojas e al. 2018). J-A Mo eno-Ruiz e al. In e na ional Jou nal o Wildland Fi e 858 The empo al accu acy o he BA-LTDR-DS p oduc in each sub- egion was assessed conside ing he o al calcu- la ed annual BA. A iming dis ibu ion o he BA-LTDR-DS p oduc was ep esen ed on a cha oge he wi h he ime se ies o e e ence BAs, and a co ela ion analysis was ca - ied ou . The ela i e pe cen ages o he annual BA o he BA-LTDR-DS p oduc we e calcula ed wi h espec o he e e ence da a o he common yea s when a ailable. Fo he spa ial accu acy assessmen , sca e plo s o he annual BA p opo ions on 50 × 50 km g ids dis ibu ed uni- o mly o he BA-LTDR-DS da ase agains he e e ence da a in each sub- egion we e cons uc ed and a linea eg ession analysis was pe o med. Nex , a de ailed analysis o he spa ial accu acy o he BA-LTDR-DS da ase was made based on e o ma ixes e sus he e e ence maps a he pixel le el on an annual basis, calcula ing commission and omission e o s o he bu ned class (S ehman 1997). Omission e o s we e calcula ed as he a io o bu ned pixels classi ied as unbu ned o he o al bu ned pixels in he BA e e ence map, while commission e o s we e calcula ed as he a io o unbu ned pixels classi ied as bu ned o he o al bu ned pixels in each BA p oduc unde analysis. We con- side ed a pixel size o 50 km o p e en e o s de i ed om geo- e e encing o he images due o he di e ence in spa ial esolu ion (Mo eno-Ruiz e al. 2014a). Resul s Annual bu ned a ea maps o he no he n bo eal egion (LTDR-BA-DB) Fig. 3 p esen s he g ouping by decade o he annual maps o BAs o he wo bo eal sub- egions conside ed (No h Ame ica and Eu asia) ob ained om he LTDR-BA-DS p od- uc wi h he bu ned pixels ep esen ed in ed. Fig. 4 shows a composi ion wi h he ou decades o he BA de ec ed by BA- LTDR-DS o he no he n bo eal egion conside ed be ween he pa allels 60°N and 72.5°N. Annual dis ibu ion o he bu ned a ea es ima es Fig. 5 shows he es ima ed annual dis ibu ion o BA in he no he n bo eal egion o he pe iod 1982–2020 om he BA-LTDR-DS. A non-uni o m pa e n was obse ed wi h yea s whe e s ong i e ac i i y was de ec ed (BA > 3 million ha) and o he yea s whe e ba ely 0.5 million ha BA was de ec ed. On a e age, o he en i e egion and pe iod, ~1.84 million ha bu ned pe yea , bu wi h high a iabili y. Fo example, in 2014 ( he yea wi h he g ea es BA de ec ed), 5.41 million ha bu ned whe eas in 1992 ( he yea wi h he leas BA de ec ed), ~0.26 million ha bu ned. Fo he en i e pe iod analysed, No h Ame ica had 35.9% o he BA compa ed wi h 64.1% in Eu asia. Howe e , his a e age con ibu ion o each sub- egion seems o ha e no s a is ical signi icance owing o he la ge annual luc ua ions. I should be no ed, o example, ha he g ea es imbalances occu ed in he yea 2020, wi h a con ibu ion o 0.8% om No h Ame ica compa ed wi h 99.2% om Eu asia, o he opposi e case o 2004, whe e No h Ame ica con ibu ed 96.7% o he o al compa ed wi h 3.3% om he es . Tempo al accu acy No h Ame ica In he No h Ame ica bo eal sub- egion, he o al BA egis e ed by he AFS and CNFDB + NBAC da abases in he 1982–2020 pe iod was 41.23 million ha, wi h an i egula annual dis ibu ion. The highes i e ac i i y (4.79 million ha) occu ed in 2004, and he lowes in 1984, when only 0.10 million ha bu ned. Fig. 6 shows he annual dis ibu ion o BA om he e e ence AFS and CNFDB + NBAC da abases, he es ima ed BA o he BA-LTDR-DS p oduc o he pe iod 1982–2020 and es ima es o he MCD64A1 C6 p oduc o he pe iod 2000–2020. The BA-LTDR-DS de ec ed 62% o he e e ence BA, unde es ima ing he BA in all yea s in he ime se ies. The e is a s ong co ela ion (0.93) be ween he LTDR- BA-DS da a and he e e ence da a. I we di ide he ime se ies in o wo pa s, p e-MODIS (1982–1999) and MODIS (2000–2020), he BA-LTDR-DS p esen s almos homogeneous beha iou , bo h in he pe cen age o he es ima e o BA (66% s 60%) and in he co ela ion coe icien (0.97 s 0.93) wi h espec o he e e ence da a. Fo i s pa , he MCD64A1 C6 p oduc unde es ima ed BA in No h Ame ica by app oxi- ma ely 69% wi h a co ela ion coe icien o 0.99 wi h espec o he e e ence da a. Making an in e -compa ison o he BA-LTDR-DS wi h MCD64A1 C6 o he common pe iod (2000–2020), a co ela ion o 0.93 was ound. Eu asia Fo he ime accu acy assessmen in he Eu asian sub- egion, a co ela ion analysis be ween he ime se ies o he BA-LTDR-DS and he MCD64A1 C6 p oduc was conduc ed (Fig. 7). MCD64A1 C6 was used as he e e ence da a o he common yea s (2000–2020). MCD64A1 C6 es ima ed a BA o 47.10 million ha o hese yea s and he BA-LTDR-DS ~66% o ha alue. The BA-LTDR-DS unde es ima ed he BA in all common yea s excep o 2000. Howe e , i s empo al pa e n i s ema kably well wi h ha o he e e - ence, yielding a co ela ion coe icien be ween MCD64A1 C6 and BA-LTDR-DS o 0.95. Spa ial accu acy No h Ame ica Table 2 shows he esul s o he linea eg ession analysis o he BA pe cen ages o he BA-LTDR-DS and he MCD64A1 C6 p oduc s e sus he e e ence da a. The a e age de e mi- na ion coe icien (R 2 ) o e he 1982–2020 pe iod o he BA-LTDR-DS p oduc was 0.78, wi h a slope o 0.73, wi h no www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 859 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 70°N 65°N 60°N 70°N 65°N 60°N 60°W 50°W 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 60°W 50°W (a) No h Ame ica 1980s 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 70°N 65°N 60°N 70°N 65°N 60°N 60°W 50°W 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 60°W 50°W (b) No h Ame ica 1990s 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 70°N 65°N 60°N 70°N 65°N 60°N 60°W 50°W 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 60°W 50°W (c) No h Ame ica 2000s 160°W 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 70°N 65°N 60°N 70°N 65°N 60°N 70°N 65°N 60°N 70°N 65°N 60°N 60°W 50°W 160°W 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E 150°W 140°W 130°W 120°W 110°W 100°W 90°W 80°W 70°W 60°W 50°W (d) No h Ame ica 2010s (e) Eu asia 1980s 70°N 65°N 60°N 70°N 65°N 60°N 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E ( ) Eu asia 1990s 70°N 65°N 60°N 70°N 65°N 60°N 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E ( ) Eu asia 1990s 70°N 65°N 60°N 70°N 65°N 60°N 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E 10°E 20°E 30°E 40°E 50°E 60°E 70°E 80°E 90°E 100°E 110°E 120°E 130°E 140°E 150°E 160°E 170°E 180°E (h) Eu asia 2010s Fig. 3. Decadal bu ned a ea maps o No h Ame ica (a–d), and Eu asia (e–h) sub- egions om he Bu ned A ea Long-Te m Da a Reco d Da aSe (BA-LTDR-DS). Red, bu ned a ea; blue, wa e ; g een, non-bu ned. J-A Mo eno-Ruiz e al. In e na ional Jou nal o Wildland Fi e 860 signi ican di e ences be ween he wo ime sub-in e als (p e and MODIS e as). The a e age de e mina ion coe icien (R 2 ) o e he 2000–2020 pe iod o MCD64A1 C6 and he BA-LTDR-DS p oduc s was 0.88 and 0.78, wi h a slope o 0.69 and 0.69, espec i ely. A de ailed analysis o he spa ial accu acy o he No h Ame ican sub- egion on an annual basis o he BA p oduc s BA-LTDR-DS and MCD64A1 C6 is shown in Table 3. The a e age commission and omission e o s o he BA-LTDR- DS a e 0.14 and 0.47 espec i ely o he en i e s udy pe iod whe eas o he MODIS e a, he MCD64A1 C6 p oduc p es- en s a commission e o o 0.09 and omission e o o 0.37. Eu asia Table 4 shows he esul s o he accu acy o he BA es ima e ob ained om he linea eg ession analysis o he BA pe cen ages om he BA-LTDR-DS p oduc e sus MCD64A1 C6 using 50 × 50 km g ids. The a e age de e mi- na ion coe icien (R 2 ) o e he 2001–2020 pe iod o he BA-LTDR-DS p oduc was 0.78, wi h a slope o 0.81. Table 5 shows he commission and omission e o s de i ed om he e o ma ix o each yea , aking as e e ence he MCD64A1 C6 p oduc (only a ailable om he yea 2000). 140°W 160°W 160°E 140°E 130°E 110°E 90°E 70°E 50°E 130°W 110°W 90°W 70°W 50°W 40°W 1980s 1990s 2000s 2010s 20°W 20°E 40°E0° 0° Fig. 4. Mapping o bu ned a eas de ec ed by BA-LTDR-DS (Bu ned A ea Long-Te m Da a-Reco d Da ase ) om 1982 o 2020 o he no he n bo eal egion be ween 60°N and 72.5°N. Colou s co e- spond o bu ned a eas by decade. AFS + (CNFDB _NBAC) BA-LTDR-DS MCD64A1 C6 0 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 1 2 3 Bu ned a ea (million. ha) 4 5 Yea Fig. 6. Annual dis ibu ion o bu ned a ea es ima e (ha) in he No h Ame ican bo eal sub- egion om e e ence da a o AFS (Alaska Fi e Se ice) and CNFDB + NBAC (Canadian Na ional Fi e Da abase + Na ional Bu ned A ea Composi e), and he BA-LTDR-DS (Bu ned A ea Long-Te m Da a-Reco d Da ase ) and he MODIS Collec ion 6 MCD64A1 C6 bu ned a ea p oduc s. 0 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 1 2 3 Bu ned a ea (million. ha) 4 5Eu asia No h Ame ica Yea 6 Fig. 5. Annual dis ibu ion o bu ned a ea es ima e (ha) in he no he n bo eal egion om he BA-LTDR-DS (Bu ned A ea Long-Te m Da a-Reco d Da ase ). www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 861 Fo his pe iod, he a e age commission and omission e o s o he BA-LTDR p oduc we e 0.21 and 0.47, espec i ely. Discussion In his pape , we desc ibe he BA-LTDR-DS da ase , which p o ides annual BA a a spa ial esolu ion o 0.05° om 1982 o 2020 o he bo eal egion be ween 60°N and 72.5°N. We buil his da ase om he la es eleased Ve sion 5 o he LTDR da ase . The LTDR Ve sion 5 inco po a es imp o e- men s o he Bidi ec ional Re lec ance Dis ibu ion Func ion (BRDF) co ec ion, he calib a ion o AVHRR/3 da a on boa d pla o ms NOAA-16, 18 and 19, composi ing a mosphe ic co ec ions and he QAs (h ps://landweb.modaps.eosdis. nasa.go /cgi-bin/l d /l d /l d Page.cgi? ileName=LTDR_ upda e). Howe e , as has been discussed in p e ious s udies (O ón e al. 2019, 2021; Giglio and Roy 2022), hese e ec s canno be elimina ed. Indeed, in his new e sion, we de ec ed missing and w ong alues o he TOA b igh ness empe a u e bands om he NOAA-16 and 18 sa elli es ( om 2000 o 2008) in he no he n bo eal egion. To p ese e he cohe ence o he empe a u e bands in he LTDR da ase , we eplaced he Band T3 wi h he Band 20 B igh ness Tempe a u e (3.360–3.840 μm) o he MOD09CMG p oduc . In addi ion, he signi ican ly decaying o bi o he NOAA-19 sa elli e om 2018 (Julien and Sob ino 2021; Giglio and Roy 2022) made i necessa y o use he MOD09CMG p oduc o his pe iod. Finally, he QA bi was upda ed p ope ly. Using his new modi ied e sion o he LTDR, we ha e gene a ed he longes BA ime se ies ye buil a a spa ial esolu ion o 0.05° in he CMG o he no he n bo eal egion. To do his, a machine lea ning algo i hm based on a Bayesian ne wo k was used, de eloped speci ically o he de ec ion o BA in ha egion. The BA-LTDR-DS ex ends by mo e han 10 yea s he ime in e al o he Global Fi e Emissions Da abase (GFED4), which is om 1995 o he p esen (Giglio e al. 2013), and imp o es i s spa ial esolu ion by up o i e imes ( om 0.25° o 0.05°) as well as he be a long- e m BA da ase (Fi eCCILT1.0) de eloped by he Clima e Change Ini ia i e (CCI) p og am o he Eu opean Space Agency. In a ecen e sion o he CCI p og am, O ón e al. (2021) ob ained a new p oduc (Fi eCCILT11) a he same spa ial esolu ion as he BA-LTDR-DS. This p oduc uses a andom o es algo i hm ha calcula es he pe cen ages o BA o each pixel (so classi ica ion), unlike BA-LTDR-DS, which only de e mines whe he he pixel is comple ely bu ned o no (ha d classi i- ca ion). The main p oblem encoun e ed when compa ing BA-LTDR_DS wi h FIRECCILT11 ela es o he yea s compos- ing bo h ime se ies: Fi eCCILT11 is 2 yea s sho e (ending in 2018) and does no include he yea 1994, which conside a- bly dis o s he compa ison as 1994 was he yea wi h he la ges BA in No h Ame ica. E alua ing he spa ial and empo al accu acy o he en i e ime se ies o BA ob ained (BA-LTDR-DS) seems e y di icul gi en he non-exis ence o ano he se o e - e ence da a o he bo eal o es egion and he pe iod analysed ( om 1982 o 2020) compa ed wi h simila p od- uc s, ei he ob ained om i e eco ds pe ime e s o o icial agencies o om p oduc s de i ed om sa elli e images. Tha is why we we e o ced o e alua e ou p oduc o he wo sub- egions al eady desc ibed (Eu asia and No h Ame ica) and wo di e en pe iods ma ked by he yea 2000 when he MODIS senso was pu in o ope a ion. Fo he assessmen o he BA-LTDR-DS accu acy in he no he n sub- egion o No h Ame ica, he bes a ailable baseline da a (g ound- u h) we e used, i.e. he pe ime e s o BAs eco ded by he AFS and Canada (CFSFND + NBAC). Few coun ies ha e de ailed egis ies o BA pe ime e s a ailable o use o a ull assessmen o p oduc s esul ing om sa elli e images. The a ailabili y o he abo e in o ma- ion made i possible o conduc a de ailed s udy o he a ious p oduc s based on he o al amoun o i es a he han a simple sample (Mo eno-Ruiz e al. 2019). These da a- bases, al hough egula ly main ained and upda ed, may con ain e o s mainly due o he omission o small i es 0 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 1 2 3 Bu ned a ea (million. ha) 4 5BA-LTDR-DS MCD64A1 C6 Yea Fig. 7. Annual dis ibu ion o bu ned a ea es ima e (ha) in he Eu asian bo eal sub- egion om he BA-LTDR-DS (Bu ned A ea Long-Te m Da a-Reco d Da ase ) and he MODIS Collec ion 6 MCD64A1 C6 bu ned a ea p oduc s. J-A Mo eno-Ruiz e al. In e na ional Jou nal o Wildland Fi e 862 chemis y model. Jou nal o Geophysical Resea ch: A mosphe es 101, 21051–21065. doi:10.1029/96JD01835 de G oo WJ, Can in AS, Flannigan MD, Soja AJ, Gowman LM, Newbe y A (2013) A compa ison o Canadian and Russian bo eal o es i e egimes. Fo es Ecology and Managemen 294, 23–34. doi:10.1016/j. o eco.2012.07.033 Duncan BN (2003) In e annual and seasonal a iabili y o biomass bu ning emissions cons ained by sa elli e obse a ions. Jou nal o Geophysical Resea ch 108, 4040. doi:10.1029/2002JD002378 Dwye E, Pinnock S, G egoi e JM, Pe ei a JMC (2000) Global spa ial and empo al dis ibu ion o ege a ion i e as de e mined om sa elli e obse a ions. In e na ional Jou nal o Remo e Sensing 21, 1289–1302. doi:10.1080/014311600210182 Ebe le J, U ban M, Homolka A, Hü ich C, Schmullius C (2016) Mul i- Sou ce Da a In eg a ion and Analysis o Land Moni o ing in Sibe ia. In ‘No el Me hods o Moni o ing and Managing Land and Wa e Resou ces in Sibe ia’. (Eds L Muelle , AK Sheudshen, F Eulens ein) pp. 471–487. (Sp inge In e na ional Publishing: Cham) Eckdahl JA, K is ensen JA, Me cal e DB (2022) Clima ic a ia ion d i es loss and es uc u ing o ca bon and ni ogen in bo eal o es wild i e. Biogeosciences 19, 2487–2506. doi:10.5194/bg-19-2487-2022 Espinola M, Pied a-Fe nandez JA, Ayala R, I iba ne L, Wang JZ (2015) Con ex ual and hie a chical classi ica ion o sa elli e images based on cellula au oma a. IEEE T ansac ions on Geoscience and Remo e Sensing 53, 795–809. doi:10.1109/TGRS.2014.2328634 Flannigan MD, Logan KA, Ami o BD, Skinne WR, S ocks BJ (2005) Fu u e a ea bu ned in Canada. Clima ic Change 72, 1–16. doi:10.1007/ s10584-005-5935-y Flannigan MD, K awchuk MA, de G oo WJ, Wo on BM, Gowman LM (2009) Implica ions o changing clima e o global wildland i e. In e na ional Jou nal o Wildland Fi e 18, 483–507. doi:10.1071/ WF08187 F anquesa M, Rod iguez-Mon ellano AM, Chu ieco E, Aguado I (2022a) Re e ence da a accu acy impac s bu ned a ea p oduc alida ion: he ole o he expe analys . Remo e Sensing 14, 4354. doi:10.3390/ s14174354 F anquesa M, S ehman SV, Chu ieco E (2022b) Assessmen and cha - ac e iza ion o sou ces o e o impac ing he accu acy o global bu ned a ea p oduc s. Remo e Sensing o En i onmen 280, 113214. doi:10.1016/j. se.2022.113214 Fuchs H, Magdon P, Kleinn C, Flessa H (2009) Es ima ing abo eg ound ca bon in a ca chmen o he Sibe ian o es und a: Combining sa el- li e image y and ield in en o y. Remo e Sensing o En i onmen 113, 518–531. doi:10.1016/j. se.2008.07.017 Ga cía-Láza o JR, Mo eno-Ruiz JA, Riaño D, A belo M (2018) Es ima ion o bu ned a ea in he no heas e n Sibe ian bo eal o es om a Long-Te m Da a Reco d (LTDR) 1982–2015 ime se ies. Remo e Sensing 10, 940. doi:10.3390/ s10060940 Geo giadi AG, Milyuko a IP, Kashu ina EA (2010) En i onmen al change in Sibe ia. En i onmen al Change in Sibe ia: Ea h Obse a ion, Field S udies and Modelling 40, 157–169. doi:10.1007/ 978-90-481-8641-9 Giglio L, Roy DP (2022) Assessmen o sa elli e o bi -d i a i ac s in he long- e m AVHRR Fi eCCILT11 global bu ned a ea da a se . Science o Remo e Sensing 5, 100044. doi:10.1016/j.s s.2022.100044 Giglio L, Rande son JT, an de We GR, Kasibha la PS, Colla z GJ, Mo on DC, DeF ies RS (2009) Assessing a iabili y and long- e m ends in bu ned a ea by me ging mul iple sa elli e i e p oduc s. Biogeosciences Discussions 6, 11577–11622. doi:10.5194/bgd-6-11577-2009 Giglio L, Rande son JT, an de We GR (2013) Analysis o daily, mon hly, and annual bu ned a ea using he ou h-gene a ion global i e emissions da abase (GFED4). Jou nal o Geophysical Resea ch: Biogeosciences 118, 317–328. doi:10.1002/jg g.20042 Giglio L, Bosche i L, Roy DP, Humbe ML, Jus ice CO (2018) The Collec ion 6 MODIS bu ned a ea mapping algo i hm and p oduc . Remo e Sensing o En i onmen 217, 72–85. doi:10.1016/j. se.2018. 08.005 Goldamme JG, Fu yae V (1996) ‘Fi e in ecosys ems o bo eal Eu asia.’ (Sp inge Science+Business Media, B.V, Kluwe Academic Publishe s) Go ham E (1991) No he n pea lands: Role in he ca bon cycle and p obable esponses o clima ic wa ming. Ecological Applica ions 1, 182–195. doi:10.2307/1941811 Guindos-Rojas F, A belo M, Ga cía-Láza o JR, Mo eno-Ruiz JA, He nández-Leal PA (2018) E alua ion o a Bayesian algo i hm o de ec bu ned a eas in he Cana y Islands’ D y Woodlands and o es s eco egion using MODIS da a. Remo e Sensing 10, 789. doi:10.3390/ s10050789 Haas O, P en ice IC, Ha ison SP (2022) Global en i onmen al con ols on wild i e bu n a ea, size, and in ensi y. En i onmen al Resea ch Le e s 17, 065004. doi:10.1088/1748-9326/ac6a69 Hall RJ, Skakun RS, Me sa an a JM, Land y R, F ase RH, Raymond D, Ga ell M, Decke V, Li le J (2020) Gene a ing annual es ima es o o es i e dis u bance in Canada: he Na ional Bu ned A ea Composi e. In e na ional Jou nal o Wildland Fi e 29(10), 878–891. doi:10.1071/WF19201 Hanes CC, Wang X, Jain P, Pa isien M-A, Li le JM, Flannigan MD (2019) Fi e- egime changes in Canada o e he las hal cen u y. Canadian Jou nal o Fo es Resea ch 49, 256–269. doi:10.1139/cj - 2018-0293 Jobbágy EG, Jackson RB (2000) The e ical dis ibu ion o soil o ganic ca bon and i s ela ion o clima e and ege a ion. Ecological Applica ions 10, 423–436. doi:10.1890/1051-0761(2000)010[0423: TVDOSO]2.0.CO;2 Julien Y, Sob ino JA (2021) NOAA-AVHRR o bi al d i co ec ion: Valida ing me hods using msg-se i i da a as a benchma k da ase . Remo e Sensing 13, 925. doi:10.3390/ s13050925 Kasischke ES, F ench NHF (1995) Loca ing and es ima ing he a eal ex en o wild i es in alaskan bo eal o es s using mul iple-season AVHRR NDVI composi e da a. Remo e Sensing o En i onmen 51, 263–275. doi:10.1016/0034-4257(93)00074-J Kasischke ES, Ch is ensen J NL, S ocks BJ (1995) Fi e, global wa ming, and he ca bon balance o bo eal o es s. Ecological Applica ions 5, 437–451. doi:10.2307/1942034 Kasischke ES, Hye EJ, No elli PC, B uhwile LP, F ench NHF, Sukhinin AI, Hewson JH, S ocks BJ (2005) In luences o bo eal i e emissions on No he n Hemisphe e a mosphe ic ca bon and ca bon monoxide. Global Biogeochemical Cycles 19, GB1012. doi:10.1029/2004GB002300 Kasischke ES, Loboda T, Giglio L, F ench NHF, Hoy EE, de Jong B, Riano D (2011) Quan i ying bu ned a ea o No h Ame ican o es s: Implica ions o di ec educ ion o ca bon s ocks. Jou nal o Geophysical Resea ch: Biogeosciences 116, G04003. doi:10.1029/ 2011JG001707 Kelly R, Chipman ML, Higue a PE, S e ano a I, B ubake LB, Hu FS (2013) Recen bu ning o bo eal o es s exceeds i e egime limi s o he pas 10,000 yea s. P oceedings o he Na ional Academy o Sciences 110, 13055–13060. doi:10.1073/pnas.1305069110 K ylo A, McCa y JL, Po apo P, Loboda T, Tyuka ina A, Tu ubano a S, Hansen MC (2014) Remo e sensing es ima es o s and- eplacemen i es in Russia, 2002–2011. En i onmen al Resea ch Le e s 9, 105007. doi:10.1088/1748-9326/9/10/105007 Kuka skaya EA, Soja AJ, Pe ko AP, Ponoma e EI, I ano a GA, Cona d SG (2013) Fi e emissions es ima es in Sibe ia: e alua ion o unce ain ies in a ea bu ned, land co e , and uel consump ion. Canadian Jou nal o Fo es Resea ch 43, 493–506. doi:10.1139/cj - 2012-0367 Ku z WA, S inson G, Rampley GJ, Dymond CC, Neilson ET (2008) Risk o na u al dis u bances makes u u e con ibu ion o Canada’s o es s o he global ca bon cycle highly unce ain. P oceedings o he Na ional Academy o Sciences 105, 1551–1555. doi:10.1073/pnas. 0708133105 Loboda TV, Zhang Z, O’Neal KJ, Sun G, Csisza IA, Shuga HH, She man NJ (2012) Recons uc ing dis u bance his o y using sa elli e-based assessmen o he dis ibu ion o land co e in he Russian Fa Eas . Remo e Sensing o En i onmen 118, 241–248. doi:10.1016/j. se.2011.11.022 Loboda TV, F ench NHF, High -Ha C, Jenkins L, Mille ME (2013) Mapping i e ex en and bu n se e i y in Alaskan ussock und a: An analysis o he spec al esponse o und a ege a ion o wildland i e. Remo e Sensing o En i onmen 134, 194–209. doi:10.1016/j. se. 2013.03.003 Melchio e A, Bosche i L (2018) Global analysis o bu ned a ea pe sis - ence ime wi h MODIS Da a. Remo e Sensing 10, 750. doi:10.3390/ s10050750 Moja adi B, Lucas C, Va shosaz M (2004) Using lea ning cellula au oma a o pos classi ica ion sa elli e image y. In e na ional www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 869 A chi es o he Pho og amme y, Remo e Sensing and Spa ial In o ma ion Sciences - ISPRS A chi es 35, 991–995. Mo eno Ruiz JA, Riaño D, A belo M, F ench NHF, Us in SL, Whi ing ML (2012) Bu ned a ea mapping ime se ies in Canada (1984–1999) om NOAA-AVHRR LTDR: A compa ison wi h o he emo e sensing p oduc s and i e pe ime e s. Remo e Sensing o En i onmen 117, 407–414. doi:10.1016/j. se.2011.10.017 Mo eno-Ruiz JA, Ga cía-Láza o JR, del Águila Cano I, He nández-Leal P (2014a) Bu ned a ea mapping in he No h Ame ican bo eal o es using e a-MODIS LTDR (2001–2011): A compa ison wi h he MCD45A1, MCD64A1 and BA GEOLAND-2 p oduc s. Remo e Sensing 6, 815–840. doi:10.3390/ s6010815 Mo eno-Ruiz JA, Ga cia-Laza o JR, Riano D, Ke au e SC (2014b) The syne gy o he 0.05° (~5km) AVHRR long- e m da a eco d (LTDR) and landsa TM a chi e o map la ge i es in he No h Ame ican bo eal egion om 1984 o 1998. IEEE Jou nal o Selec ed Topics in Applied Ea h Obse a ions and Remo e Sensing 7, 1157–1166. doi:10.1109/JSTARS.2013.2292853 Mo eno-Ruiz JA, Ga cía-Láza o JR, A belo M, Riaño D (2019) A com- pa ison o bu ned a ea ime se ies in he Alaskan bo eal o es s om di e en emo e sensing p oduc s. Fo es s 10, 363. doi:10.3390/ 10050363 Mo eno-Ruiz JA, Ga cía-Láza o JR, A belo M, Can ón-Ga bín M (2020) MODIS senso capabili y o bu ned a ea mapping—assessmen o pe o mance and imp o emen s p o ided by he la es s anda d p od- uc s in bo eal egions. Senso s 20, 5423. doi:10.3390/s20185423 Mo ise e JT, Ba e F, P i e e JL, Myneni RB, Nickeson JE, Ga igues S, Shabano NV, Weiss M, Fe nandes RA, Leblanc SG, Kalacska M, Sanchez-Azo ei a GA, Chubey M, Ri a d B, S enbe g P, Rau iainen M, Voipio P, Manninen T, Pilan AN, Lewis TE, Iiames JS, Colombo R, Me oni M, Buse o L, Cohen WB, Tu ne DP, Wa ne ED, Pe e sen GW, Seu e G, Cook R (2006) Valida ion o Global Mode a e- Resolu ion LAI P oduc s: A F amewo k P oposed Wi hin he CEOS Land P oduc Valida ion Subg oup. IEEE T ansac ions on Geoscience and Remo e Sensing 44, 1804–1817. doi:10.1109/TGRS.2006.872529 Mouillo F, Schul z MG, Yue C, Cadule P, Tansey K, Ciais P, Chu ieco E (2014) Ten yea s o global bu ned a ea p oduc s om spacebo ne emo e sensing – A e iew: Analysis o use needs and ecommenda- ions o u u e de elopmen s. In e na ional Jou nal o Applied Ea h Obse a ion and Geoin o ma ion 26, 64–79. doi:10.1016/j.jag.2013. 05.014 Nelson K, Thompson D, Hopkinson C, Pe one R, Chasme L (2021) Pea land- i e in e ac ions: A e iew o wildland i e eedbacks and in e ac ions in Canadian bo eal pea lands. Science o The To al En i onmen 769, 145212. doi:10.1016/j.sci o en .2021.145212 Núñez-Casillas L, Ga cía Láza o JR, Mo eno-Ruiz JA, A belo M (2013) A compa a i e analysis o bu ned a ea da ase s in Canadian bo eal o es in 2000. The Scien i ic Wo ld Jou nal 2013, 289056. doi:10.1155/2013/289056 O ón G, Ramo R, Lizundia-Loiola J, Chu ieco E (2019) Global de ec ion o long- e m (1982–2017) bu ned a ea wi h AVHRR-LTDR da a. Remo e Sensing 11, 2079. doi:10.3390/ s11182079 O ón G, Lizundia-Loiola J, Pe ina i ML, Chu ieco E (2021) De elopmen o a consis en global long- e m bu ned a ea p oduc (1982–2018) based on AVHRR-LTDR da a. In e na ional Jou nal o Applied Ea h Obse a ion and Geoin o ma ion 103, 102473. doi:10.1016/j.jag.2021.102473 Padilla M, S ehman SV, Ramo R, Co i D, Han son S, Oli a P, Alonso- Canas I, B adley AV, Tansey K, Mo a B, Pe ei a JM, Chu ieco E (2015) Compa ing he accu acies o emo e sensing global bu ned a ea p od- uc s using s a i ied andom sampling and es ima ion. Remo e Sensing o En i onmen 160, 114–121. doi:10.1016/j. se.2015.01.005 Pan Y, Bi dsey RA, Fang J, Hough on R, Kauppi PE, Ku z WA, Phillips OL, Sh idenko A, Lewis SL, Canadell JG, Ciais P, Jackson RB, Pacala SW, McGui e AD, Piao S, Rau iainen A, Si ch S, Hayes D (2011) A la ge and pe sis en ca bon sink in he wo ld’s o es s. Science 333, 988–993. doi:10.1126/science.1201609 Pa isien MA, Pe e s VS, Wang Y, Li le JM, Bosch EM, S ocks BJ (2006) Spa ial pa e ns o o es i es in Canada 1980–1999. In e na ional Jou nal o Wildland Fi e 15, 361–374. doi:10.1071/WF06009 Pedel y J, De adiga S, Masuoka E, B own M, Pinzon J, Tucke C, Ve mo e E, P ince S, Nagol J, Jus ice C, Roy D, Ju J, Schaa C, Liu J, P i e e J, Pinhei o A (2007) Gene a ing a long- e m land da a eco d om he AVHRR and MODIS ins umen s. In ‘In e na ional Geoscience and Remo e Sensing Symposium (IGARSS)’. pp. 1021–1025. (IEEE) doi:10.1109/IGARSS.2007.4422974 Pin y B, Ve s ae e MM (1992) GEMI: a non-linea index o moni o global ege a ion om sa elli es. Vege a io 101, 15–20. doi:10.1007/ BF00031911 Ponoma e EI, Kha uk VI, Ranson KJ (2016) Wild i es dynamics in Sibe ian la ch o es s. Fo es s 7, 125. doi:10.3390/ 7060125 Po ie J, Gau hie S, Be ge on Y (2019) Spa ial dis ibu ion o mean i e size and occu ence in eas e n Canada: in luence o clima e, physical en i onmen and ligh ning s ike densi y. In e na ional Jou nal o Wildland Fi e 28, 927–940. doi:10.1071/WF18220 Roge s BM, Soja AJ, Goulden ML, Rande son JT (2015) In luence o ee species on con inen al di e ences in bo eal i es and clima e eed- backs. Na u e Geoscience 8, 228–234. doi:10.1038/ngeo2352 Shuman JK, Shuga HH, O’Hallo an TL (2011) Sensi i i y o Sibe ian la ch o es s o clima e change. Global Change Biology 17, 2370–2384. doi:10.1111/j.1365-2486.2011.02417.x Si no SA, Mokho II (2018) A Compa a i e Analysis o he Cha ac e is ics o Ac i e Fi es in he Bo eal Fo es s o Eu asia and No h Ame ica Based on Sa elli e Da a. Iz es iya, A mosphe ic and Oceanic Physics 54, 966–978. doi:10.1134/S0001433818090347 Skakun R, Whi man E, Li le JM, Pa isien MA (2021) A ea bu ned adjus men s o his o ical wildland i es in Canada. En i onmen al Resea ch Le e s 16, 064014. doi:10.1088/1748-9326/ab b2c Skakun R, Cas illa G, Me sa an a J, Whi man E, Rod igue S, Li le J, G oenewegen K, Coyle M (2022) Ex ending he Na ional Bu ned A ea Composi e Time Se ies o Wild i es in Canada. Remo e Sensing 14, 3050. doi:10.3390/ s14133050 Soja AJ, Co e WR, Shuga HH, Sukhinin AI, S ackhouse PW, McRae DJ, Cona d SG (2004) Es ima ing i e emissions and dispa i ies in bo eal Sibe ia (1998–2002). Jou nal o Geophysical Resea ch: A mosphe es 109, D14S06. doi:10.1029/2004JD004570 Soja AJ, Al-Saadi J, Giglio L, Randall D, Ki aka C, Poulio G, Ko dzi JJ, Ra use S, Pace TG, Pie ce TE, Moo e T, Roy B, Pie ce RB, Szykman JJ (2009) Assessing sa elli e-based i e da a o use in he Na ional Emissions In en o y. Jou nal o Applied Remo e Sensing 3, 031504. doi:10.1117/1.3148859 S ehman SV (1997) Selec ing and in e p e ing measu es o hema ic classi ica ion accu acy. Remo e Sensing o En i onmen 62, 77–89. doi:10.1016/S0034-4257(97)00083-7 S ocks BJ, Mason JA, Todd JB, Bosch EM, Wo on BM, Ami o BD, Flannigan MD, Hi sch KG, Logan KA, Ma ell DL, Skinne WR (2003) La ge o es i es in Canada, 1959–1997. Jou nal o Geophysical Resea ch 108(D1), 8149. doi:10.1029/2001JD000484 S owe LL, Da is PA, McClain EP (1999) Scien i ic Basis and Ini ial E alua ion o he CLAVR-1 Global Clea /Cloud Classi ica ion Algo i hm o he Ad anced Ve y High Resolu ion Radiome e . Jou nal o A mosphe ic and Oceanic Technology 16, 656–681. doi:10.1175/ 1520-0426(1999)016<0656:SBAIEO>2.0.CO;2 Sukhinin AI, F ench NHF, Kasischke ES, Hewson JH, Soja AJ, Csisza IA, Hye EJ, Loboda T, Con ad SG, Romasko VI, Pa lichenko EA, Miski SI, Slinkina OA (2004) AVHRR-based mapping o i es in Russia: New p oduc s o i e managemen and ca bon cycle s udies. Remo e Sensing o En i onmen 93, 546–564. doi:10.1016/j. se.2004.08.011 Tchebako a NM, Pa eno a E, Soja AJ (2009) The e ec s o clima e, pe ma os and i e on ege a ion change in Sibe ia in a changing clima e. En i onmen al Resea ch Le e s 4, 045013. doi:10.1088/ 1748-9326/4/4/045013 Tu e sky MR, Kane ES, Ha den JW, O ma RD, Manies KL, Hoy E, Kasischke ES (2011) Recen accele a ion o biomass bu ning and ca bon losses in Alaskan o es s and pea lands. Na u e Geoscience 4, 27–31. doi:10.1038/ngeo1027 Vande hoo MK, Fai aux N, Beal YJG, Hawbake TJ (2017) Valida ion o he USGS Landsa bu ned a ea essen ial clima e a iable (BAECV) ac oss he con e minous Uni ed S a es. Remo e Sensing o En i onmen 198, 393–406. doi:10.1016/j. se.2017.06.025 Vi cha A (2011) Wild i es in Russia in 2000–2008: Es ima es o bu n a eas using he sa elli e MODIS MCD45 da a. Remo e Sensing Le e s 2, 81–90. doi:10.1080/01431161.2010.499138 Walsh JE (2014) In ensi ied wa ming o he A c ic: Causes and impac s on middle la i udes. Global and Plane a y Change 117, 52–63. doi:10.1016/j.gloplacha.2014.03.003 J-A Mo eno-Ruiz e al. In e na ional Jou nal o Wildland Fi e 870 Wiede RK, Vi DH, Bensco e BW (2006) Pea lands and he bo eal o es . In ‘Bo eal pea land ecosys ems’. (Eds RK Wiede , DH Vi ) pp. 1–8. (Sp inge -Ve lag: Be lin, Heidelbe g) Wi h C (2005) Fi e egime and ee di e si y in bo eal o es s: implica ions o he ca bon cycle Fo es Di e si y and Func ion. In ‘Fo es Di e si y and Func ion, Ecological S udies. Vol. 176’. (Eds M Sche e -Lo enzen, C Kö ne , ED Schulze) pp. 309–44. (Sp inge : Be lin, Heidelbe g) Woos e MJ, Zhang YH (2004) Bo eal o es i es bu n less in ensely in Russia han in No h Ame ica. Geophysical Resea ch Le e s 31, L20505. doi:10.1029/2004GL020805 Da a a ailabili y. The da a ha suppo his s udy a e a ailable in he a icle and accompanying online supplemen a y ma e ial. Con lic s o in e es . The au ho s decla e no con lic o in e es . Decla a ion o unding. This esea ch was unded by he Minis e io de Ciencia, Inno ación y Uni e sidades (MCIU), he Agencia Es a al de In es igación (AEI) and he Fondo Eu opeo de Desa ollo Regional (FEDER) h ough p ojec RTI2018-099171-B-I00. Acknowledgemen s. The au ho s wish o hank he LTDR p ojec and i s eam o making he da a a ailable, and Alaska Fi e Se ice, he Canada Cen e o Mapping and Ea h Obse a ion o Na u al Resou ces Canada, he Canadian Fo es Se ice, NASA, NOAA and USGS o p ocessing and dis ibu ing he AFS, CFNDB, NBAC, MCD64A1 and LTDR da ase s. Au ho con ibu ions. J.A.M.-R. and J.R.G.-L. concei ed, designed and applied he me hodology. All au ho s ob ained, analysed and discussed he esul s; M.A. and P.A.H.-L. in collabo a ion wi h he es o he au ho s w o e and con ibu ed o he edi ing o he manusc ip . Au ho a ilia ions A Depa amen o de In o má ica, Uni e sidad de Alme ía, 04120 Alme ía, Spain. Email: [email p o ec ed]; [email p o ec ed] B Depa amen o de Física, Uni e sidad de La Laguna, 38200 San C is óbal de La Laguna, Spain. Email: [email p o ec ed]; [email p o ec ed] www.publish.csi o.au/w In e na ional Jou nal o Wildland Fi e 871