曇 娜，阿拉騰圖婭，包玉龍，高彥哲，敖日格樂

（內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院, 內(nèi)蒙古 呼和浩特 010022）

火作為一種自然生態(tài)因子，對草地生態(tài)系統(tǒng)具有積極或消極的影響，不僅影響著草地生態(tài)系統(tǒng)發(fā)展、演替和消亡過程[1]，具有偶然發(fā)生、迅速發(fā)展、難以控制的特點，對防范能力弱的草原牧區(qū)生命與財產(chǎn)帶來非常嚴(yán)重的威脅，給畜牧業(yè)為主的草原地區(qū)經(jīng)濟(jì)發(fā)展、社會安定帶來巨大影響。蒙古高原草原區(qū)是歐亞大陸溫帶草原的主要組成部分，也是火行為極為活躍的地區(qū)[2-3]，特別是東部區(qū)。據(jù)統(tǒng)計，1994 年春，錫林郭勒盟東烏珠穆沁旗發(fā)生3 起特大火災(zāi)，燒死大小牲畜共0.93 萬頭，造成直接經(jīng)濟(jì)損失5 330萬元，間接經(jīng)濟(jì)損失 82.5 萬元[4]。1996 年和1997 年的2 月—6 月，蒙古國發(fā)生大面積的森林草原火災(zāi)，共計29 人死亡，82 人受傷，1.17 萬頭牲畜死亡，多處基礎(chǔ)設(shè)施、以及0.23 萬hm2牧場和森林被燒毀，財產(chǎn)損失總計172.29 萬元[5]。因此，探究蒙古高原東部草原火的時空演變模式，以此為據(jù)制定草原火災(zāi)管理對策、減災(zāi)規(guī)劃，部署防災(zāi)、抗災(zāi)、救災(zāi)工作和草原發(fā)展規(guī)劃是學(xué)術(shù)界和管理部門的迫切任務(wù)。

關(guān)于蒙古高原草原火分布研究已有較多成果，分布在新巴爾虎地區(qū)[6]、中蒙邊境地區(qū)[7]、蒙古高原[8]、呼倫貝爾地區(qū)[9-10]、錫林郭勒盟地區(qū)[11]、內(nèi)蒙古自治區(qū)東部和牧區(qū)[12-14]。上述草原火的時空分布研究均采用時間序列和空間統(tǒng)計分析法，即使用過程線表達(dá)某一區(qū)域草原火發(fā)生時間變化特征，用地理信息系統(tǒng)空間統(tǒng)計法研究某個時間斷面上的草原火發(fā)生區(qū)域差異特征。對草原火的空間特征與時間因素有效融合考慮方面不夠充分，將草原火的空間分布與隨時間變化特征割裂開來，不能有效地挖掘草原火的時空分布演變趨勢。然而，草原火是在特定條件下才發(fā)生。滿足草原火發(fā)生的各種條件不僅不同地區(qū)有差異，同一地區(qū)不同時間段也有很大差異，隨機性、復(fù)雜性非常突出。因此，要想精確了解各點位的草原火發(fā)生模式則需要將空間點位和時間序列一一對應(yīng)尋找規(guī)律。同時，其他突發(fā)性災(zāi)難事件也警示著需要用網(wǎng)格化管理方式精準(zhǔn)把握突發(fā)事件的發(fā)生發(fā)展規(guī)律，采取精準(zhǔn)的防范措施。Purwanto 等[15]和Mo 等[16]分別以印度尼西亞東爪哇和中國為研究區(qū)，基于時空立方體分析技術(shù)，識別新型冠狀病毒感染的時空聚類格局，確定冷、熱點趨勢，并展示疫情可視化輸出結(jié)果，以直觀觀察疫情高發(fā)區(qū)和低風(fēng)險區(qū)，為下一步防控策略提供依據(jù)。這一技術(shù)可將整個區(qū)域劃分成多個網(wǎng)格，大區(qū)域火事件碎化成小網(wǎng)格火事件，同時發(fā)現(xiàn)每個網(wǎng)格點位時間上的歷史規(guī)律性，強調(diào)各網(wǎng)格不同警種間的聯(lián)合、以及空間上相鄰點位之間的關(guān)聯(lián)，實現(xiàn)草原火研究更進(jìn)一步精細(xì)化，竭力化危為零、化危為機提供新的思路。

時空立方體作為一種時空數(shù)據(jù)模型，能夠融合處理地理現(xiàn)象的時間、空間與屬性信息，實現(xiàn)數(shù)據(jù)的歷史狀態(tài)重建、時空變化跟蹤及發(fā)展態(tài)勢預(yù)測[17]，可以保證時空數(shù)據(jù)的連續(xù)性[18]?；诳臻g自相關(guān)理論和時間序列分析方法可以有效挖掘地理現(xiàn)象的空間依賴關(guān)系及時間依賴關(guān)系[19]。綜上，本研究輸入蒙古高原草原火高頻區(qū)2001—2020 年每次發(fā)生的草原火的空間數(shù)據(jù)，獲得時空立方體，利用相鄰空間上火點發(fā)生次數(shù)的關(guān)聯(lián)獲得草原火異常值與集聚強度，借助趨勢分析方法對立方體每個網(wǎng)格時間條柱上的火點集聚強度獲得草原火時空分布的熱點、冷點趨勢評估可為具有科學(xué)的理論邏輯、深厚的歷史邏輯與直觀的實踐邏輯草原火災(zāi)網(wǎng)格化管理和集中力量精準(zhǔn)防范提供輔助決策參考。

1 材料與方法

1.1 研究區(qū)概況

蒙古高原地處歐亞大陸腹地，位于37°46′～53°08′ N，87°40′～122°15′ E。東起大興安嶺、西至阿爾泰山，北界為薩彥嶺、雅布洛諾夫山脈，南界為陰山山脈[20]。因東西跨度大，自然條件差異懸殊而草原火發(fā)生面積和頻率也不同。按蒙古國的各省和內(nèi)蒙古的盟市統(tǒng)計結(jié)果，近20 年各行政單位平均每年過火次數(shù)32 次，其中蒙古國3 個省、內(nèi)蒙古2 個市(盟)的平均每年過火次數(shù)高達(dá)184 次，明顯高于其他地區(qū)的頻次。因此，本研究選中國內(nèi)蒙古自治區(qū)呼倫貝爾市、錫林郭勒盟和蒙古國東方省、蘇赫巴托爾省和肯特省為研究區(qū)(圖1)，總面積73.92 萬km2。研究區(qū)整體地勢自西向東逐漸降低，平均海拔1 200～1 500 m。植被類型以草原為主，從東北向西南依次按草甸草原、典型草原和荒漠草原演替，呈明顯的地帶性。屬溫帶大陸性氣候，干燥、風(fēng)大，降水量少且集中于夏季。處于東亞季風(fēng)邊緣區(qū)，同時受西風(fēng)帶環(huán)流的影響，春季干旱多風(fēng)，夏季高溫多雨[21]。每年10 月左右植被開始枯萎凋落，可燃物量增多，冬季被雪覆蓋，次年的5 月初草原植被返青期前枯草存量較豐富而秋季和春季易引發(fā)草原火災(zāi)[22]。

圖1 研究區(qū)地理位置圖Figure 1 Geographical location of the study area

1.2 數(shù)據(jù)來源及處理

采用的遙感衛(wèi)星數(shù)據(jù)為2001—2020 年MCD64A1火燒跡地產(chǎn)品數(shù)據(jù)集，該數(shù)據(jù)集由NASA 網(wǎng)站提供(http://ladsweb.nascom.nasa.gov)。MCD64A1 是MODIS月尺度火燒跡地數(shù)據(jù)，其空間分辨率為500 m，Sinusoidal 投影，數(shù)據(jù)格式為HDF。數(shù)據(jù)由燃燒日期、不確定的燃燒日期、質(zhì)量評估、開始日期和結(jié)束日期等信息組成。本研究中利用MRT 工具，對MCD64A1 數(shù)據(jù)進(jìn)行預(yù)處理(數(shù)據(jù)拼接、投影轉(zhuǎn)換、數(shù)據(jù)格式轉(zhuǎn)換)，再利用ArcGIS 軟件按研究區(qū)范圍進(jìn)行裁剪。

1.3 研究方法

1.3.1 時空立方體模型

Hagerstrand[23]將時間概念引入空間數(shù)據(jù)中，并首先提出使用時空立方體模型來表示地理現(xiàn)象隨時間的變化。時空立方體模型在使用二維空間表達(dá)地理對象的傳統(tǒng)方法基礎(chǔ)上，增加了一維時間軸，用于表達(dá)地理實體或現(xiàn)象隨著時間的變化(圖2)。本研究中的草原火時空立方體模型采用二維坐標(biāo)軸來表示草原火的平面位置，再加一維縱向的時間軸表示平面位置上草原火點數(shù)量隨時間的變化，從而構(gòu)建一個由平面空間和時間軸組成的時空融合草原火立方體。通過統(tǒng)計各時空立方體中包含的草原火點數(shù)目，將草原火點聚合到每個空間網(wǎng)格的時空條柱中。這樣則可運用時間維的幾何特性來表達(dá)研究區(qū)內(nèi)各空間點位上的草原火隨時間的變化過程。

圖2 時空立方體模型Figure 2 Space-time cube model

利用 ArcGIS Pro 軟件提供的時空立方體創(chuàng)建工具，將草原火點聚合到時空條柱中，眾多的時空條柱組成一個三維的時空立方體，使用網(wǎng)絡(luò)公用數(shù)據(jù)格式(network common data format, NetCDF)存儲。構(gòu)建時空立方體時需要確定時間步長和空間步長。本研究中采用控制變量法對空間步長5、10、15、20 km和時間步長1、3、6、9、12 個月進(jìn)行試驗探究，發(fā)現(xiàn)以較小的空間步長，較大的時間步長聚合火點，會挖掘出更多的冷熱點(表1)。并通過對火燒跡地面積的統(tǒng)計發(fā)現(xiàn)平均面積為27.90 km2；以周期長度為時間間隔創(chuàng)建時空立方體，能夠排除周期性變化對時空演變模式中趨勢分析的干擾。從按月統(tǒng)計的蒙古高原草原火發(fā)生時間序列發(fā)現(xiàn)以一年為周期的規(guī)律性變化很明顯(圖3)。因此，將空間步長設(shè)置為5 km × 5 km，時間步長設(shè)置為1 年(12 個月)。

表1 蒙古高原高頻區(qū)不同時間步長、空間步長的探索性結(jié)果Table 1 Exploratory analysis results of different time and space steps in the high-frequency area of Mongolian Plateau

圖3 蒙古高原高頻區(qū)草原火2001—2020 年月平均頻率圖Figure 3 Monthly average frequency of grassland fires in high-frequency areas of Mongolian Plateau from 2001 to 2020

1.3.2 局部異常值分析

局部異常值分析可確定數(shù)據(jù)中的顯著聚類和異常值。在時空立方體中查找顯著聚類和異常值時通過局部空間自相關(guān)統(tǒng)計[24]來實現(xiàn)每個時空條柱時間與空間上與其鄰域間存在的統(tǒng)計差異分。即選擇時空模式挖掘功能模塊按鄰域距離和鄰域時間步長參數(shù)估計各立方體條柱的局部空間自相關(guān)統(tǒng)計量，檢測局部異常值，以無顯著性模式、高值—高值聚類模式、高值—低值異常模式、低值—低值聚類模式、低值—高值異常模式、多種類型模式6 種模式展示結(jié)果。

1.3.3 時空熱點分析

時空熱點分析方法以空間熱點分析和時間序列分析為基礎(chǔ)，利用時空立方體模型將時空關(guān)系納入到空間熱點分析中，探測時空中存在的熱點和冷點。基于時間序列分析的方法評估這些熱點和冷點隨時間的變化趨勢，分辨出新增、連續(xù)、加強、持續(xù)、逐漸減少、分散、震蕩、歷史以及無顯著性等不同變化趨勢的熱點或冷點。時空熱點分析使用Getis-Ord G*i統(tǒng)計[25]探測要素空間分布的熱點和冷點，檢測熱點或冷點的時間序列變化趨勢則是基于Mann-Kendall[26]趨勢分析法來實現(xiàn)。

2 結(jié)果與分析

2.1 草原火時空分布特征

2.1.1 時間分布特征

在2001—2020 年間，蒙古高原草原火高頻區(qū)過火次數(shù)及過火面積存在明顯的年際變化，表現(xiàn)為峰谷交替，每隔 3～4 年出現(xiàn)一次草原火發(fā)生的高峰，但總體呈現(xiàn)出波動中下降趨勢(圖4)。2001—2020年蒙古高原草原火高頻區(qū)共發(fā)生草原火18 138 起，累積過火面積達(dá)20.00 萬 km2。草原火過火次數(shù)較高的年份為2003、2006、2008、2014、2015 年。其中，2006 年過火次數(shù)最多，達(dá)0.18 萬次；過火次數(shù)最少的年份為2001 年，為417 次。2003、2007、2011、2012、2015 年是草原火過火面積較大的年份，其中，2003 年的草原火過火面積最大，達(dá)2.13 萬 km2；2018年的過火面積最小，僅為0.11 萬 km2。

圖4 2001—2020 年蒙古高原高頻區(qū)草原火年際變化Figure 4 Interannual variation of grassland fires in highfrequency areas of Mongolian Plateau from 2001 to 2020

2.1.2 空間分布特征

從2001—2020 年蒙古高原草原火高頻區(qū)燃燒面積統(tǒng)計結(jié)果看(圖5)，蒙古國草原火燒跡地明顯多于中國內(nèi)蒙古自治區(qū)。在蒙古國，東方省是研究區(qū)內(nèi)草原火過火面積最大的省，在研究時段內(nèi)過火面積達(dá)13.81 萬 km2，占整個研究區(qū)總過火面積的69%；蘇和巴托爾省和肯特省過火面積分別達(dá)2.53 萬和1.31 萬 km2。在中國內(nèi)蒙古自治區(qū)，火燒跡地主要分布在呼倫貝爾市、錫林郭勒盟東北部等地區(qū)。呼倫貝爾市是中國內(nèi)蒙古自治區(qū)內(nèi)草原火過火面積最大的市，過火面積達(dá)1.97 萬 km2。錫林郭勒盟過火面積達(dá)3.21 萬 km2。

圖5 2001—2020 年蒙古高原高頻區(qū)草原火燒跡地分布圖Figure 5 Distribution of burned grassland in high-frequency areas of Mongolian Plateau from 2001 to 2020

2.2 草原火時空立方體創(chuàng)建

以12 個月為時間步長、5 km 為空間步長創(chuàng)建的草原火時空立方體模型如圖6 所示，時空立方體在20 個時間步長間隔內(nèi)將 160.53 萬點聚合到 6.65萬格網(wǎng)位置。每個位置均是5 km × 5 km 的方形。整個時空立方體由西向東跨越0.12 萬 km ，由南向北跨越0.13 萬 km 的區(qū)域。時間軸上每個時段步長間隔為12 個月，因此時空立方體覆蓋的整個時間段為240 個月。

圖6 蒙古高原高頻區(qū)草原火時空立方體Figure 6 Space-time cube of grassland fires in highfrequency areas of Mongolian Plateau

2.3 草原火局部異常值分析

由于目前對時空立方體模型中的最優(yōu)鄰域時間和鄰域距離尚未形成明確的標(biāo)準(zhǔn)和依據(jù)，因此本研究采用控制變量法對不同的鄰域時間和鄰域距離進(jìn)行試驗。通過分析生成時空演變模式的數(shù)量和種類，并以生成時空演變模式的占比和數(shù)量，對比選擇出鄰域時間和鄰域距離[27]。對鄰域時間24～60個月和鄰域距離5～20 km 進(jìn)行試驗，發(fā)現(xiàn)當(dāng)鄰域時間為24 個月且鄰域距離為20 km 時，時空演變模式占比最高，且生成的時空聚類模式數(shù)量最多(圖7)。這一時空鄰域也符合當(dāng)年的草原火通過可燃物的消耗和植物的生長來影響次年的火發(fā)生概率，再往后滯后作用不大以及草原火以其發(fā)生網(wǎng)格為中心向各方向影響周圍最鄰近網(wǎng)格位置點的可能性最大的實際情況。因此，選擇鄰域時間為24 個月(鄰域時間步長為1)，鄰域距離為20 km。從草原火局部異常值格網(wǎng)統(tǒng)計結(jié)果來看(圖8)，主要為多種類型模式，其次為低—低聚類模式和無顯著性模式以及少量的高—高聚類模式、高—低異常模式和低—高異常模式。

圖7 蒙古高原高頻區(qū)不同鄰域時間、鄰域距離局部異常值分布圖Figure 7 Distribution of local outliers with different neighborhood times and neighborhood distances in high-frequency areas of Mongolian Plateau

圖8 蒙古高原高頻區(qū)草原火局部異常值格網(wǎng)統(tǒng)計圖Figure 8 Grid statistics of local outliers of grassland fires in high-frequency areas of Mongolian Plateau

從草原火聚類與異常值的空間分布來看(圖7)，呼倫貝爾市莫力達(dá)瓦達(dá)斡爾族自治旗紅彥鎮(zhèn)、哈達(dá)陽鎮(zhèn)、額爾和鄉(xiāng)、騰克鎮(zhèn)和登特科鎮(zhèn)出現(xiàn)高—高聚類模式，說明周圍鄰近地區(qū)是草原火高發(fā)區(qū)，其中該網(wǎng)點是高發(fā)區(qū)中的高發(fā)點。呼倫貝爾市中西部、東北部及東南部，錫林郭勒盟東北部、東部及東南角，東方省中部，肯特省西北部、中部及東南部，蘇赫巴托爾省北部部分地區(qū)出現(xiàn)低—低聚類模式，說明草原火發(fā)生一直維持在相對較低的水平，沒有出現(xiàn)異常點。東方省巴彥東縣和色楞格縣，肯特省賓代爾縣、德勒格爾汗縣、達(dá)爾罕縣和巴彥塔格縣，呼倫貝爾市鄂倫春自治旗克一河鎮(zhèn)、額爾古納市恩和哈達(dá)鎮(zhèn)及恩和俄羅斯族民族鄉(xiāng)、鄂溫克族自治旗錫尼河西蘇木、新巴爾虎左旗烏布爾寶力格蘇木和新巴爾虎右旗呼倫鎮(zhèn)，錫林郭勒盟東烏珠穆沁旗烏里雅斯太鎮(zhèn)和呼熱圖淖爾蘇木、西烏珠穆沁旗吉仁高勒鎮(zhèn)、錫林浩特市寶力根蘇木出現(xiàn)高—低異常模式，出現(xiàn)的高頻發(fā)生的異常點主要是沿公路人口集聚地，與周圍地區(qū)相比草原火發(fā)生頻率明顯高。東方省哈馬爾達(dá)板縣，蘇赫巴托爾省蒙赫汗縣和阿斯嘎特縣，呼倫貝爾市莫力達(dá)瓦達(dá)斡爾族自治旗紅彥鎮(zhèn)、奎勒河鎮(zhèn)、塔溫敖寶鎮(zhèn)、登特科鎮(zhèn)、巴彥鄂溫克民族鄉(xiāng)和額爾和鄉(xiāng)，錫林郭勒盟東烏珠穆沁旗滿都胡寶拉格鎮(zhèn)和呼熱圖淖爾蘇木、阿巴嘎旗吉爾嘎朗圖蘇木出現(xiàn)低—高異常模式，說明這些點為防火系統(tǒng)較完善的人口聚集區(qū)，而周圍地區(qū)多為森林草原火高發(fā)區(qū)。除此之外，20 年間草原火發(fā)生面積最大、最集中連片的區(qū)域呈現(xiàn)出多種類型模式的特征，主要包括呼倫貝爾市中西部和東部，東方省大部分地區(qū)、肯特省北部、蘇赫巴托省北部和東北部，錫林郭勒盟東烏珠穆沁旗等地區(qū)，說明這些網(wǎng)格小區(qū)域和周圍相鄰地區(qū)出現(xiàn)過草原火發(fā)生頻率均高或均低的聚類模式，也出現(xiàn)過網(wǎng)格小區(qū)域草原火發(fā)生頻率高于或低于相鄰地區(qū)的個別異常點，因此該點在總體上會表現(xiàn)出多種類型模式。

2.4 草原火時空熱點分析

與局部異常值分析相同，采用控制變量法對鄰域距離5～20 km 和鄰域時間24～60 個月進(jìn)行試驗，發(fā)現(xiàn)當(dāng)鄰域距離為20 km 且鄰域時間為48 個月時，時空演變模式占比最高且生成的時空熱點模式數(shù)量最多(圖9)。這一鄰域時間也符合每隔 3～4 年出現(xiàn)一次草原火發(fā)生高峰的時間周期。因此，選擇鄰域時間為48 個月(鄰域時間步長為3)，鄰域距離為20 km。從草原火時空熱點格網(wǎng)統(tǒng)計結(jié)果來看(圖10)，熱點演變模式的類別主要是振蕩熱點、新增熱點和分散熱點，還有個別的連續(xù)熱點。主要的冷點演變模式類別是振蕩冷點和分散冷點，此外還有連續(xù)冷點以及少量的新增冷點、加強冷點和持續(xù)冷點。無加強熱點、持續(xù)熱點、逐漸減少冷/熱點和歷史冷/熱點等模式類別。

圖9 蒙古高原高頻區(qū)不同鄰域時間、鄰域距離時空熱點分布圖Figure 9 Spatiotemporal distribution of hot spots with different neighborhood times and neighborhood distances in high-frequency areas of Mongolian Plateau

圖10 蒙古高原高頻區(qū)草原火時空熱點格網(wǎng)統(tǒng)計圖Figure 10 Grid statistics of spatiotemporal hot spots of grassland fires in high-frequency areas of Mongolian Plateau

蒙古高原高頻區(qū)草原火時空冷熱點演變模式中(圖9)，草原火的熱點主要分布在呼倫貝爾市東部、陳巴爾虎旗和新巴爾虎左旗國境線附近，肯特省北部，東方省西南角、東北部、東南部以及與呼倫貝爾市新巴爾虎左旗接壤處。其中，在呼倫貝爾市鄂倫春自治旗北部和阿榮旗東部與莫力達(dá)瓦達(dá)斡爾族自治旗南部交界處出現(xiàn)新增熱點，說明研究時段前期草原火發(fā)生較少但在最近草原火發(fā)生頻率較高。連續(xù)熱點主要分布在莫力達(dá)瓦達(dá)斡爾族自治旗額爾和鄉(xiāng)西南部和新巴爾虎左旗北部，說明草原火發(fā)生頻率較高且連續(xù)發(fā)生。分散熱點主要分布在莫力達(dá)瓦達(dá)斡爾族自治旗北部、陳巴爾虎旗國境線附近和東方省西南部國境線附近，說明草原火頻繁發(fā)生，特點為斷斷續(xù)續(xù)不連續(xù)出現(xiàn)。振蕩熱點的覆蓋范圍較廣，主要分布在呼倫貝爾市鄂倫春自治旗西南部，陳巴爾虎旗和新巴爾虎左旗國境線附近，東方省喬巴山縣東部、巴彥烏拉縣與巴彥東縣交界處、哈馬爾達(dá)板縣西南角、南部與錫林郭勒盟東烏珠穆沁旗交界處和東北部與呼倫貝爾市新巴爾虎左旗交界處，說明在研究時段內(nèi)交替出現(xiàn)高頻危險時期和低頻相對安全時期，研究期間內(nèi)近期草原火發(fā)生頻率較高。

草原火發(fā)生頻率低的冷點模式分布范圍較廣泛。其中，新增冷點主要分布在東方省哈馬爾達(dá)板縣中東部和東部與呼倫貝爾市新巴爾虎左旗接壤的國境線附近，說明草原火發(fā)生頻率在最近時間段內(nèi)顯著下降。在呼倫貝爾市鄂倫春自治旗北部國境線附近、額爾古納市北部、陳巴爾虎旗中北部、新巴爾虎左旗東南部與鄂溫克自治旗南部交界處，錫林郭勒盟正藍(lán)旗東南部與多倫縣西部交界處、太仆寺旗東南部呈連續(xù)冷點趨勢，說明草原火發(fā)生頻率較低，未曾間斷。加強冷點主要分布在呼倫貝爾市阿榮旗東南部，說明草原火發(fā)生頻率一直較低，近期進(jìn)一步降低呈強化趨勢。在呼倫貝爾市陳巴爾虎旗巴彥哈達(dá)蘇木東部、阿榮旗南部，錫林郭勒盟多倫縣蔡木山鄉(xiāng)西部呈持續(xù)冷點趨勢，說明在90%的研究時段內(nèi)草原火發(fā)生頻率較低。分散冷點主要分布在呼倫貝爾市中西部和東部部分地區(qū)，肯特省西北角和東部與蘇赫巴托爾省西北部交界處，蘇赫巴托爾省東部部分地區(qū)，錫林郭勒盟東南角，說明以上地區(qū)的草原火斷斷續(xù)續(xù)的偶爾發(fā)生，時序上冷點出現(xiàn)的規(guī)律性較差。在東方省大部分地區(qū)、肯特省北部、蘇赫巴托省東部和北部、呼倫貝爾市東部和錫林郭勒盟東北部部分地區(qū)呈振蕩冷點趨勢，說明以上地區(qū)草原火低頻發(fā)生時期具有周期性特點。

3 討論

已有關(guān)于蒙古高原草原火的研究成果較多，如曲炤鵬等[8]基于地理信息系統(tǒng)空間分析方法研究了2000—2007 年蒙古高原草原火行為的時空分布規(guī)律，研究者認(rèn)為蒙古國草原火多發(fā)地帶主要分布在東部與內(nèi)蒙古自治區(qū)交界處的草甸草原和蘇赫巴托爾省內(nèi)；內(nèi)蒙古自治區(qū)草原火多發(fā)區(qū)主要分布在呼倫貝爾地區(qū)和錫林郭勒盟中部和東北部地區(qū)。包剛[28]等采用時間序列和空間統(tǒng)計分析方法對蒙古高原火行為時空變化趨勢進(jìn)行了研究，發(fā)現(xiàn)2001—2012 年蒙古國多發(fā)地帶主要分布在東方省、蘇赫巴托爾省和肯特省等地區(qū)，在內(nèi)蒙古自治區(qū)呼倫貝爾草原區(qū)具有較集中的過火跡地，并且在錫林郭勒草原也有較頻繁發(fā)生的過火跡地。徐書興[29]等采用空間統(tǒng)計分析方法研究了2001—2017 年蒙古東部野火時空動態(tài)，結(jié)果顯示17 年間蒙古東部野火主要發(fā)生在蒙古東北部地區(qū)，特別是在蒙古國—中國沿邊地區(qū)。野火發(fā)生區(qū)域主要集中在東方省、肯特省和蘇赫巴托爾省。本研究中也得到了與采用傳統(tǒng)時空分布方法研究結(jié)果相同的蒙古高原草原火高頻區(qū)空間分布和隨時間變化的區(qū)域特征，但傳統(tǒng)時空分布方法僅考慮宏觀特征，未考慮時空積聚性和時空演變動態(tài)模式，對精準(zhǔn)管理還不能提供有效信息。

傳統(tǒng)的空間統(tǒng)計分析是以區(qū)域化變量為信息，通過空間統(tǒng)計手段，研究地理現(xiàn)象的空間變異與空間結(jié)構(gòu)[30]，輸入分析的是某個時間斷面或平均值數(shù)據(jù)。而時空立方體模型在傳統(tǒng)的統(tǒng)計方法上加入時間維度的考量，整合時空信息和時空關(guān)系，進(jìn)一步將“鄰近度”概念從“空間臨近”擴(kuò)展為“時空鄰近”[31]，研究整個地理現(xiàn)象隨時空的推移所表現(xiàn)出的復(fù)雜趨勢變化。與傳統(tǒng)的空間統(tǒng)計方法相比局部異常值分析和新興時空熱點分析類別增多，并且根據(jù)設(shè)定的鄰域距離和鄰域時間步長參數(shù)值在每個條柱上進(jìn)行分析，其精確度提高。為此，本研究基于時空立方體模型進(jìn)一步深入分析草原火點的時空分布特點，將時空數(shù)據(jù)挖掘的方法應(yīng)用到草原火點變化分析研究中，彌補了傳統(tǒng)空間統(tǒng)計將空間、時間或者屬性分開分析的不足，實現(xiàn)了空間、時間和屬性信息的完整表達(dá)，有效揭示了草原火現(xiàn)象的時空特征和規(guī)律。這可為準(zhǔn)確、快速、精細(xì)地識別草原火多發(fā)點位及其時空演變特性，網(wǎng)格化管理草原火，管護(hù)員分配巡查調(diào)度、草原火預(yù)測、降低火災(zāi)發(fā)生頻率等提供輔助決策，也可推廣應(yīng)用到其他領(lǐng)域。

草原火行為是一種受自然條件和人文因素共同影響的危害性事件。中國內(nèi)蒙古自治區(qū)和蒙古國草原火發(fā)生面積的差異懸殊特點表明，植被狀況可能是主要原因。人類活動也已成為影響草原火行為的重要因素。研究證明，放牧對草原火發(fā)生具有一定抑制作用[32]。中國內(nèi)蒙古自治區(qū)牧民生活方式已實現(xiàn)由游牧向定居轉(zhuǎn)變，牲畜數(shù)量增加，使得草地實際承載力加大，枯落物減少；而蒙古國仍保持著傳統(tǒng)的游牧方式，草原植被休養(yǎng)時間長，草原退化相對較輕[33]，地表覆蓋的枯落物豐富，為草原火的發(fā)生提供了豐富的燃料。另外，中國內(nèi)蒙古自治區(qū)嚴(yán)格的管制措施和撲救力度，使得該地區(qū)草原火行為受到強烈抑制；而蒙古國人口密度低，人力和財力有限，火災(zāi)撲救和管理措施相對落后[32]，當(dāng)?shù)厝嗣駥Σ菰鸬陌l(fā)生多采取放任的態(tài)度，導(dǎo)致該地區(qū)草原火蔓延迅速，過火面積大且頻率高。本研究雖然利用時空立方體、局部異常值分析和時空冷熱點分析挖掘蒙古高原草原火高頻區(qū)時空演變模式，也發(fā)現(xiàn)了些人口聚集相關(guān)的個別熱點和冷點，但未能全面揭示草原火時空分布的影響機制。從網(wǎng)格化精細(xì)管理視角，重視歷史邏輯性和空間自相關(guān)性的基礎(chǔ)上，有待進(jìn)一步深入研究各因素對于草原火發(fā)生點和波及區(qū)的具體影響程度。

4 結(jié)論

利用蒙古高原草原火高頻區(qū)2001—2020 年MCD64A1 火燒跡地數(shù)據(jù)，挖掘分析研究區(qū)草原火的時空分布模式，得出以下結(jié)論：

1) 蒙古高原草原火高頻區(qū)整體上呈從西南向東北草原火發(fā)生頻率與面積逐漸增多的趨勢，蒙古國和中國內(nèi)蒙古自治區(qū)境內(nèi)的草原火發(fā)生面積差異懸殊，草原火主要分布在蒙古國境內(nèi)。時間上，年際波動較大，整體呈波動中下降的趨勢。

2) 研究區(qū)草原火發(fā)生面積最大區(qū)的空間和時間自相關(guān)關(guān)系比較復(fù)雜，規(guī)律性差，高—高聚類、高—低異常、低—低聚類、低—高異常各模式均出現(xiàn)過，防火減災(zāi)難度大。而草原火發(fā)生頻率顯著高或低的個別點主要分布在人口聚集的鄉(xiāng)鎮(zhèn)所在地，這表明著草原火發(fā)生和防范與人口密度和行為有關(guān)。

3) 草原火發(fā)生面積最大區(qū)各網(wǎng)格上雖然草原火發(fā)生頻率低，但不連續(xù)、不定期地發(fā)生火，波及面廣。呼倫貝爾市東部、西部國境線附近，肯特省北部，東方省西南角、東北部、東南部以及與呼倫貝爾市新巴爾虎左旗交界處草原火雖也不連續(xù)或不定期發(fā)生，但近期發(fā)生頻率較高，需引起關(guān)注。呼倫貝爾市東部是草原火頻繁發(fā)生區(qū)域，需要高度警惕。