民國時(shí)期四川省霍亂流行的時(shí)空規(guī)律及驅(qū)動(dòng)因素

閆晉博, 張 濤*, 龔勝生

(1.華中師范大學(xué)地理過程分析與模擬湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430079;2.華中師范大學(xué)可持續(xù)發(fā)展研究中心, 武漢 430079)

霍亂是由霍亂弧菌引起的急性消化道傳染病,是我國法定的兩種甲類傳染病之一.自19世紀(jì)以來全球暴發(fā)了7次霍亂大流行,是威脅全球人類健康的重大公共衛(wèi)生問題.據(jù)史料記載,霍亂于清朝嘉慶末年(1820年)傳入我國[1],在醫(yī)療水平落后、衛(wèi)生意識(shí)薄弱的晚清、民國時(shí)期,霍亂流行讓百姓的生活雪上加霜.歷史健康地理學(xué)關(guān)注歷史時(shí)期疾病的時(shí)空擴(kuò)散規(guī)律和人地相關(guān)性,針對(duì)歷史時(shí)期的傳染病、疫災(zāi)等的研究已取得顯著成果[2],包括3 000年來我國總體、朝代和區(qū)域疫災(zāi)時(shí)空規(guī)律和影響因素研究[3-5]及災(zāi)害鏈分析[6].晚清、民國時(shí)期的疫病研究更加細(xì)化,更多關(guān)注疫病的傳播與防控措施[7-10].具體到霍亂,李孜沫[11]認(rèn)為清代霍亂在中國的適生區(qū)為東部沿海地帶,疫病分布自東向西遞減、南方高于北方,具有明顯的流行分界;張萍分析了1932年陜西霍亂傳播與環(huán)境擾動(dòng)的關(guān)系[12],張萍[13]、孟文科等[14]、單麗[15]進(jìn)一步探究了霍亂及防控對(duì)地區(qū)近代化的影響;李玉尚等探究了商業(yè)型聚落[16]、戰(zhàn)爭起義[17]與霍亂疫情的關(guān)系;李春曉等探究了民國時(shí)期重慶霍亂的主要傳播因素以及政府的應(yīng)對(duì)和控制模式[18].霍亂與地理環(huán)境之間的關(guān)系錯(cuò)綜復(fù)雜,溫度、降水、海拔和到海岸線的距離對(duì)霍亂的分布有顯著影響[19-20].霍亂是民國時(shí)期四川省發(fā)生次數(shù)最多、影響最廣的疫病,常常朝患夕亡[21]、全家盡數(shù)死絕[22].從歷史健康地理學(xué)視角出發(fā),基于疫災(zāi)史料和環(huán)境資料,對(duì)民國時(shí)期四川省霍亂的時(shí)空規(guī)律及驅(qū)動(dòng)因素進(jìn)行分析,可為烈性傳染病的防控提供歷史經(jīng)驗(yàn).

1 研究區(qū)域、數(shù)據(jù)來源與研究方法

1.1 研究區(qū)域

四川省位于中國西南腹地,青藏高原和長江中下游平原間的過渡地帶,地勢(shì)西高東低,中東部的四川盆地為其主體區(qū)域,西部為川西高原.以亞熱帶濕潤氣候?yàn)橹?川西高原山地區(qū)為高原山地氣候.地處長江上游,主要河流包括雅礱江、岷江、大渡河、嘉陵江等.1925年四川省共轄5道146縣,包括西川道31縣、東川道36縣、建昌道28縣、永寧道25縣、嘉陵道26縣(圖1).自古以來四川省發(fā)生過大小戰(zhàn)爭數(shù)百次,隨著國民政府的西遷,重慶成為陪都,四川省成為當(dāng)時(shí)國家的戰(zhàn)略大后方.

注:參照《四川省明細(xì)全圖》(上海商務(wù)印書館1927年版)和《中國分省新地圖》(亞光輿地學(xué)社1947年版,第24、32-33頁)改繪.圖3～圖7同.圖1 民國時(shí)期四川省行政區(qū)劃圖Fig.1 Map of administrative divisions of Sichuan Province during the period of the Republic of China (1912-1949)

1.2 數(shù)據(jù)來源

基于《中國三千年疫災(zāi)史料匯編》(清代卷、民國卷(上)、民國卷(下))[23],以縣為單位梳理統(tǒng)計(jì)民國時(shí)期四川霍亂疫情史料.由于記載偏差,民國時(shí)期四川省的部分痢疾疫情實(shí)質(zhì)上也是霍亂,將郫縣(1928年)[24]、平武縣(1935年)[25]兩次明確為霍亂的痢疾疫情納入統(tǒng)計(jì)分析.霍亂流行的驅(qū)動(dòng)因素包括自然環(huán)境和人文環(huán)境兩方面,自然環(huán)境考慮氣候、地形、水源、洪澇災(zāi)害等因素,人文環(huán)境考慮人口、交通、戰(zhàn)亂等因素.影響因素?cái)?shù)據(jù)中,使用當(dāng)前自然環(huán)境數(shù)據(jù)模擬[26]民國時(shí)期四川省的自然地理環(huán)境,氣象數(shù)據(jù)來源于國家氣象信息中心,海拔、河道數(shù)據(jù)來自四川省及重慶市DEM數(shù)據(jù)提取和全國五級(jí)河道數(shù)據(jù)提取,洪澇災(zāi)害數(shù)據(jù)來自《四川省兩千年洪災(zāi)史料匯編》[27].人口數(shù)據(jù)來自《民國人口戶籍史料匯編》[28],交通數(shù)據(jù)來自《中國分省精圖》[29],戰(zhàn)亂資料來自《中國三千年疫災(zāi)史料匯編》,在霍亂疫情史料中進(jìn)一步梳理兵疫霍亂信息參與分析,具體見下述表1.

1.3 研究方法

時(shí)間分析中,基于霍亂疫情史料建立流行年表,參考前期研究[30-33],使用霍亂縣數(shù)、面積、頻度、廣度、厚度、趨勢(shì)、周期、季發(fā)率等指標(biāo),分析霍亂的總體流行特征、年際變化和季節(jié)變化趨勢(shì).霍亂頻度指霍亂流行年數(shù)與歷經(jīng)年數(shù)的百分比.霍亂廣度指霍亂流行的縣數(shù)或土地面積,也可用霍亂縣數(shù)與區(qū)域總縣數(shù)的百分比來表示.霍亂厚度指霍亂多年累計(jì)波及的土地面積與該區(qū)域土地總面積的比值.霍亂季發(fā)率指發(fā)生霍亂的季節(jié)與該時(shí)段總的霍亂季節(jié)數(shù)的百分比.空間分析中,使用空間自相關(guān)分析方法分析疫情空間集聚性特征.驅(qū)動(dòng)因素分析中,以斯皮爾曼相關(guān)性分析篩選出具有顯著相關(guān)性的霍亂流行驅(qū)動(dòng)因素,采用地理加權(quán)回歸分析明確影響因素的空間異質(zhì)性,采用地理探測(cè)器探測(cè)影響因素之間的交互作用.以上步驟在Excel 2019、ArcGIS 10.7、SPSS 22.0中完成.

地理加權(quán)回歸模型(GWR)可以考慮變量的空間異質(zhì)性,反映處于不同地理位置時(shí),變量對(duì)該區(qū)域的影響程度.其計(jì)算公式為:

(1)

式中,yi和xij(j=1,2,…,k)分別為點(diǎn)i處的因變量和自變量;(ui,vi)指點(diǎn)i的空間坐標(biāo),β0(ui,vi)為截距項(xiàng),βj(ui,vi)指回歸系數(shù);系數(shù)符號(hào)表示正負(fù)相關(guān)性質(zhì),大小表示相關(guān)程度的強(qiáng)弱.

地理探測(cè)器可用于分析地理事物背后的驅(qū)動(dòng)因素,因素間的交互作用探測(cè)是其特色功能[34].通過q統(tǒng)計(jì)量比較,如下式:

(2)

2 民國時(shí)期四川省霍亂流行的時(shí)空規(guī)律

民國時(shí)期,四川省共發(fā)生大小疫災(zāi)1 000余次,其中發(fā)生在夏季、秋季以及貫穿夏秋兩季的疫災(zāi)最多,占近80%;各種疫災(zāi)累計(jì)波及全域880縣次,其中排在前三的疫災(zāi)依次是霍亂(434縣次)、天花(182縣次)、傷寒(108縣次).

2.1 四川省霍亂流行的時(shí)間分布特征

2.1.1 霍亂流行的總體特征 統(tǒng)計(jì)結(jié)果顯示,在研究時(shí)間范圍內(nèi),四川省共有39年發(fā)生過霍亂,霍亂的流行頻度高達(dá)97.50%;四川省146個(gè)縣中共有121個(gè)縣遭受霍亂波及,霍亂的流行廣度為82.88%;在四川省41.86萬km2國土中,霍亂流行區(qū)累計(jì)面積達(dá)68.36萬km2,霍亂流行厚度為1.63層.從總體流行特征而言,民國時(shí)期四川省霍亂流行十分嚴(yán)重.

2.1.2 霍亂流行的年際變化 從霍亂流行的縣數(shù)與流行面積來看,民國時(shí)期四川省霍亂疫情縣數(shù)呈現(xiàn)十分明顯的上升趨勢(shì),說明了霍亂流行的范圍不斷擴(kuò)大,霍亂疫情越來越嚴(yán)重(圖2).霍亂的趨勢(shì)變化與所有疫災(zāi)的趨勢(shì)變化相吻合.自20世紀(jì)30年代開始,霍亂縣數(shù)較民國前期明顯上升,所有疫災(zāi)的區(qū)縣數(shù)也呈現(xiàn)明顯增長,這與自20世紀(jì)30年代開始的日本侵華戰(zhàn)爭造成社會(huì)動(dòng)蕩、大西遷密切相關(guān).從流行周期來看,霍亂爆發(fā)的時(shí)間周期為3～5年,前后共出現(xiàn)四個(gè)明顯的霍亂流行高峰,即1916年、1920年、1939年和1945年,也與所有疫災(zāi)峰值年份相符,進(jìn)一步驗(yàn)證霍亂在四川省疫災(zāi)中占主導(dǎo)地位.

圖2 1910-1949年四川省疫災(zāi)及霍亂變化趨勢(shì)圖Fig.2 Trend of epidemic disasters and cholera in Sichuan Province in 1910-1949

2.1.3 霍亂流行的季節(jié)分布 霍亂疫情往往來源于食物和飲用水的污染變質(zhì),對(duì)于季節(jié)環(huán)境變化十分敏感.據(jù)明確記載季節(jié)的史料統(tǒng)計(jì),研究時(shí)段159個(gè)季節(jié)中霍亂流行的季節(jié)有67個(gè),霍亂流行的季發(fā)率為42.14%,其中:流行于春季4次,占5.97%;夏季36次,占53.73%;秋季26次,占38.80%;冬季1次,占1.49%;夏、秋兩季總共占比92.54%.從季發(fā)率來看,霍亂高發(fā)于夏、秋兩季.從霍亂流行縣次的季節(jié)狀況來看,具有明確記載的329個(gè)霍亂縣次里,有320縣次分布在夏季、秋季或是貫穿夏秋兩季數(shù)月,占比高達(dá)97.26%.霍亂流行于夏、秋季節(jié),或與四川省夏、秋季節(jié)高溫多雨的氣候條件相關(guān).

2.2 四川省霍亂流行的空間分布規(guī)律

2.2.1 霍亂流行的地域演變 以1910—1919年、1920—1929年、1930—1939年、1940—1949年每十年為一段,統(tǒng)計(jì)每十年各縣霍亂流行年數(shù)和全省累計(jì)霍亂縣數(shù),分析民國時(shí)期四川省的霍亂災(zāi)情的空間分布狀況.結(jié)果顯示,四川省全域霍亂的流行頻度和廣度都呈現(xiàn)出波動(dòng)上升的特點(diǎn).

1) 1910—1919年,霍亂流行區(qū)域相對(duì)分散,呈斑塊狀散布(圖3a):以川西北、川東和川中部永寧道內(nèi)部轄縣為主要疫災(zāi)縣.此時(shí)全域霍亂流行年數(shù)維持低位,以遂寧縣流行頻度最高,十年內(nèi)有3次流行.

2) 1920—1929年,霍亂流行的區(qū)域有所收束,出現(xiàn)條帶斑塊狀分布(圖3b):川西邊陲的懋功、蘆山、天全、雅安、滎經(jīng)、峨眉、樂山7縣連成一片,集中在大渡河及其支流青衣江流域;川北由南江、蒼溪開始,沿嘉陵道中部的閬中、鹽亭、蓬溪,向川內(nèi)腹地形成帶狀霍亂分布,與川中地區(qū)的樂至、安岳、富順、瀘縣等縣連片;東川道仍呈散布態(tài)勢(shì),達(dá)縣霍亂流行頻度有所上升.

3) 1930—1939年,霍亂流行區(qū)域呈現(xiàn)集中連片態(tài)勢(shì),出現(xiàn)明顯的流行頻度高值區(qū)(圖3c):十年內(nèi)分別形成了川中以內(nèi)江、資中兩縣為核心以及川東北以宣漢、達(dá)縣兩縣為核心的兩個(gè)西北-東南連片延展霍亂流行板塊.總體而言,霍亂發(fā)生地更加向四川省中部集中.首次出現(xiàn)霍亂年數(shù)達(dá)5的縣(內(nèi)江、江北),即平均每隔一年便會(huì)遭受一次霍亂侵襲.1930年代相對(duì)前兩個(gè)十年段而言,累計(jì)霍亂縣數(shù)激增近30個(gè),說明霍亂波及范圍越來越廣.

4) 1940—1949年,霍亂流行區(qū)域在四川省中部集中連片,霍亂流行頻度和廣度都達(dá)到空前規(guī)模(圖3d):霍亂流行頻度高值區(qū)以資中、內(nèi)江、富順、榮縣等地為核心,形成沿沱江至長江分布的高頻度帶.富順和江北的流行頻度高達(dá)6次.

圖3 民國時(shí)期四川省霍亂流行空間變化Fig.3 Spatial change of cholera epidemic in Sichuan Province during the period of the Republic of China (1912-1949)

四個(gè)時(shí)段,四川省累計(jì)霍亂年數(shù)分別為57、105、103、169,累計(jì)霍亂縣數(shù)分別為50、53、80、92.從霍亂年數(shù)看,十年間發(fā)生4次及以上的疫災(zāi)縣由0個(gè)增加到10個(gè),高頻度區(qū)域范圍不斷擴(kuò)大;從霍亂縣數(shù)可知,民國時(shí)期四川省霍亂流行的范圍總體呈現(xiàn)擴(kuò)大趨勢(shì),由早期的斑塊分散狀到中后期逐漸集中連片擴(kuò)展.

2.2.2 霍亂流行的空間集聚 民國時(shí)期,霍亂在四川省的絕大部分縣皆有流行(圖4a).146縣只有鹽源、會(huì)理、峨邊、秀山、巫山等25縣未出現(xiàn)過霍亂流行,絕大部分位于川西南、川南、川東邊陲地區(qū).從霍亂流行頻度來看,以永寧道內(nèi)部各縣最為嚴(yán)重,建昌道則相對(duì)最輕.

對(duì)民國四川省霍亂進(jìn)行空間自相關(guān)分析.全局空間自相關(guān)以莫蘭指數(shù)衡量地理事物的空間集聚態(tài)勢(shì),結(jié)果顯示Moran’sI值為0.06(z=3.84,p=0),由此可知民國時(shí)期四川省的霍亂在空間上呈現(xiàn)出顯著集聚態(tài)勢(shì).通過局部空間自相關(guān)分析明確霍亂在空間上的集聚區(qū),從結(jié)果來看(圖4b),霍亂疫情的高值集聚區(qū)處于川中丘陵及成都平原的部分地區(qū),即永寧道、嘉陵道、東川道交界地帶.疫情的低值集聚區(qū)處于川西南的鹽邊、昭覺及川東南的酉陽等縣.

圖4 民國時(shí)期四川省霍亂總體空間分布與集聚區(qū)域圖Fig.4 Overall spatial distribution and agglomeration area of cholera in Sichuan Province during the period of the Republic of China (1912-1949)

3 民國時(shí)期四川省霍亂流行的驅(qū)動(dòng)因素

3.1 霍亂流行與環(huán)境因素的相關(guān)性

3.1.1 自然環(huán)境驅(qū)動(dòng)因素

1) 氣候因素.氣候決定的水熱條件是自然環(huán)境中各類生物活動(dòng)的基礎(chǔ).氣候?qū)θ祟惢顒?dòng)產(chǎn)生多階復(fù)合、連環(huán)的作用.氣候的冷暖干濕變化對(duì)疫病的流行成災(zāi)具有重要影響.在霍亂流行年數(shù)與年均溫、年降水量的相關(guān)性分析中,雖然兩者未通過顯著性檢驗(yàn),但從霍亂流行年數(shù)與年平均溫度、年平均降水量的耦合圖(圖5a)來看,霍亂的主要流行區(qū)位于年均溫17 ℃的四川盆地中南部,同時(shí)霍亂流行的重災(zāi)區(qū)主要位于年雨量1 000～1 200 mm的地區(qū).一般而言,濕熱環(huán)境十分有利于霍亂弧菌的繁殖,進(jìn)而易爆發(fā)大規(guī)模的霍亂疫情.通過民國時(shí)期四川省霍亂流行的季節(jié)性分析可知,霍亂在夏秋兩季的季發(fā)率非常高,一定程度上說明了濕熱氣候條件能夠驅(qū)動(dòng)霍亂疫情的暴發(fā).

2) 地形因素.地形起伏影響水熱條件、水系特征,進(jìn)一步影響其它自然及人文環(huán)境.從高程與霍亂流行年數(shù)的相關(guān)性分析結(jié)果來看,高程與霍亂流行年數(shù)呈現(xiàn)顯著的負(fù)相關(guān)(r=-0.306,p≤0.01).將霍亂流行年數(shù)與四川省地形圖疊合(圖5b),霍亂的流行區(qū)多為廣袤平原或低矮丘陵,而在高海拔的山地則較少發(fā)生.平原內(nèi)部地勢(shì)平坦,河網(wǎng)密集,河水流速相對(duì)較緩,流動(dòng)性差,污染物容易積聚,極易滋生細(xì)菌,誘發(fā)霍亂.此外,平原或丘陵地區(qū)優(yōu)越的水熱條件使得人類活動(dòng)豐富多樣,極易促成疫情傳播擴(kuò)散.《什邡縣志》載:“縣治方亭鎮(zhèn)地處平壩,但環(huán)境衛(wèi)生差,人口密集,流動(dòng)頻繁,都是導(dǎo)致霍亂流行因素.永興、靈杰、雙盛、馬井、禾豐等沿河一帶,當(dāng)年因遭水災(zāi),亦易霍亂蔓延”[35].

3) 水源因素.霍亂的水體傳播性,使得受污染的井水、河水往往是其擴(kuò)散的重要媒介.《威遠(yuǎn)縣志》載:“(1939年)6月,城北鋪?zhàn)訛嘲l(fā)生霍亂,病源順河而下,迅速蔓延,縣城及沿河居民先后死亡數(shù)百人”[36].選取河道長度、河網(wǎng)密度、縣治距河道的最近距離(包含五級(jí)河道)、縣治距主要河道的最近距離(僅取四級(jí)及以上河道)四組數(shù)據(jù)與霍亂流行年數(shù)進(jìn)行分析,結(jié)果顯示縣治距河道的最近距離(r=-0.230,p≤0.01)、縣治距主要河道的最近距離(r=-0.183,p≤0.05)與霍亂流行年數(shù)呈負(fù)相關(guān),河道長度、河網(wǎng)密度相關(guān)性未通過顯著性檢驗(yàn).將四級(jí)及以上河道分別作3 km和10 km緩沖區(qū)再與霍亂流行年數(shù)疊合(圖5b),霍亂流行年數(shù)大于9的高頻度縣位于3 km緩沖區(qū)內(nèi)的占66.7%,所有霍亂發(fā)生縣中有53.4%位于10 km緩沖區(qū)范圍內(nèi).

4) 洪澇災(zāi)害.遭受洪水侵襲的地區(qū)往往生靈涂炭,傷亡慘重,病菌也極易四散傳播.《射洪縣志》載:“(1945年)8月30日至31日,涪江及各支流洪水同時(shí)暴漲,沿河城鄉(xiāng)淹死2 000余人,災(zāi)后疫病(霍亂)流行”[37],又如《彭山縣志》載:“(1947年)秋,洪水泛濫,沿江一帶霍亂流行”[38].從洪澇災(zāi)害數(shù)與霍亂流行年數(shù)的相關(guān)性分析結(jié)果來看,洪澇災(zāi)害數(shù)與霍亂流行年數(shù)呈現(xiàn)顯著的正相關(guān)(r=0.303,p≤0.01).將霍亂流行年數(shù)與四川省地形圖疊合(圖5c),岷江和沱江流域洪澇災(zāi)害頻發(fā),民國時(shí)期四川省霍亂的主要流行區(qū)也多位于洪災(zāi)頻度較高的地帶.

圖5 民國時(shí)期四川省霍亂流行與地理環(huán)境的耦合Fig.5 Coupling of cholera epidemic and geographical environment in Sichuan Province during the period of the Republic of China (1912-1949)

3.1.2 人文環(huán)境驅(qū)動(dòng)因素

1) 人口分布.人口作為人地系統(tǒng)中最為活躍的要素,人口集聚對(duì)于任何傳染病而言都是疫情蔓延的催化劑,“人口流動(dòng)是疫災(zāi)流行的加速器”[39].人口數(shù)量(r=0.409,p≤0.01)、人口密度(r=0.355,p≤0.01)與霍亂年數(shù)呈現(xiàn)顯著的正相關(guān).將人口密度做等值線再與民國全歷史時(shí)期四川省的霍亂疫情頻數(shù)分布重疊(圖5d),可以發(fā)現(xiàn)人口集聚的地區(qū),霍亂疫情相對(duì)更為嚴(yán)重.

2) 交通因素.交通發(fā)展在帶來便利的同時(shí)很大程度上促進(jìn)了各類傳染病的蔓延.尤其近代以后,伴隨鐵路交通的發(fā)展,各類傳染病實(shí)現(xiàn)沿鐵路線長距離快速擴(kuò)散,其中1932年沿隴海鐵路深入關(guān)中一帶,造成西北地區(qū)暴發(fā)的大規(guī)模霍亂疫情便是很好的印證[40].《資中縣志》載:“(1947年秋)霍亂流行,交通沿線尤烈”[41].從四川省各縣交通線總里程、各縣治與主要交通線的距離與霍亂流行年數(shù)的相關(guān)性來看,交通線總里程與霍亂流行年數(shù)呈顯著正相關(guān)(r=0.323,p≤0.01),縣治與主要交通線的距離與霍亂流行年數(shù)呈顯著負(fù)相關(guān)(r=-0.293,p≤0.01),結(jié)合交通線分布狀況來看(圖5d),霍亂流行度較高的地區(qū)普遍為主要交通線沿線地區(qū),交通通達(dá)性較高.

3) 戰(zhàn)亂因素.霍亂疫情經(jīng)常伴隨行軍、駐軍、戰(zhàn)亂出現(xiàn).通過梳理史料發(fā)現(xiàn)(表1),四川省的霍亂疫情常從軍隊(duì)內(nèi)部爆發(fā)并向外傳染,這與行軍的艱苦條件以及駐軍內(nèi)部衛(wèi)生意識(shí)和條件較差密切相關(guān).1916—1920年的霍亂高發(fā)期或與當(dāng)時(shí)的軍閥混戰(zhàn)相關(guān),如《四川軍閥史料》記載:“(1917年)2月起,土匪數(shù)踞古宋城,軍閥互用匪勢(shì)傾軋……戰(zhàn)后,疫癘(霍亂)大作,傳遍數(shù)區(qū),死亡3 000余人,全家滅絕者200余戶,加上流亡轉(zhuǎn)徙,宋人殆減少大半”[42].當(dāng)時(shí)軍閥行軍、交戰(zhàn)主要位于閬中、梓潼、成都、簡陽、古宋等縣,這些縣也是當(dāng)時(shí)霍亂疫情較為流行的地區(qū).此外,戰(zhàn)亂帶來的社會(huì)動(dòng)蕩與難民潮也會(huì)對(duì)霍亂疫情蔓延起到推波助瀾的作用.《巫山縣志》載:“是年,日軍占領(lǐng)宜昌,大量難民及國民黨軍隊(duì)涌入巫山,一時(shí)瘟疫(疑似霍亂)大發(fā)”[43].

表1 民國四川省兵疫霍亂史料

3.2 霍亂驅(qū)動(dòng)因素的交互作用及其影響的空間異質(zhì)性

3.2.1 霍亂驅(qū)動(dòng)因素的交互作用 地理環(huán)境作為一個(gè)有機(jī)整體,自然環(huán)境因素和人文環(huán)境因素相互疊加、相互作用形成疫病的適生條件[33],各種社會(huì)、經(jīng)濟(jì)、政治、氣候、技術(shù)和環(huán)境因素可以塑造并影響疾病在人群中的出現(xiàn)模式[44].在環(huán)境驅(qū)動(dòng)因素相關(guān)性分析的基礎(chǔ)上,剔除未通過顯著性檢驗(yàn)的影響因子,運(yùn)用地理探測(cè)器的因子探測(cè)和交互探測(cè)探究各影響因子對(duì)霍亂疫災(zāi)空間異質(zhì)性的作用強(qiáng)度,所有影響因子均在95%或99%水平上顯著.

就單因子作用而言,影響因子作用強(qiáng)度從大到小排列依次為:人口數(shù)量>人口密度>交通線里程>高程>洪澇災(zāi)害>縣治與主要交通線距離>縣治與最近河道距離.單因子的作用中,人文社會(huì)環(huán)境因子的影響普遍強(qiáng)于自然環(huán)境的作用,對(duì)霍亂流行影響貢獻(xiàn)最大的因子都來源于人口的集聚分布,再次驗(yàn)證人口集聚之地更容易疫病流行.交通線里程及交通線距離的影響,也印證了交通要素對(duì)傳染病的擴(kuò)散作用,與史料及相關(guān)研究結(jié)論一致.洪災(zāi)與河道距離也對(duì)霍亂流行具有重要影響,霍亂作為經(jīng)水傳播的疫病,民國時(shí)期四川省居民的飲水主要來源于井水或河水,在洪澇災(zāi)害之后,污穢腐物進(jìn)入水體,病菌在水中大量滋生,極易通過飲水等途徑進(jìn)入人體導(dǎo)致疫災(zāi)爆發(fā).

從影響因子的交互作用的結(jié)果來看,任意兩個(gè)影響因子之間的交互作用均為非線性增強(qiáng)或雙因子增強(qiáng).自然環(huán)境因子與人文環(huán)境因子的交互增強(qiáng)效果尤為明顯,其中以洪澇災(zāi)害次數(shù)與人口數(shù)量的交互作用最為強(qiáng)勁,p值由分別作用的0.105和0.270上升至交互作用的0.560.鄰水地帶依托水源發(fā)展耕作和交通運(yùn)輸,人口集聚、經(jīng)濟(jì)發(fā)展,鄰水聚落迅速興旺,交通條件進(jìn)一步改善,人口流動(dòng)進(jìn)一步增強(qiáng).但與此同時(shí),霍亂等疫災(zāi)爆發(fā)和傳播的風(fēng)險(xiǎn)隨之增大.大災(zāi)之后必有大疫,當(dāng)經(jīng)歷洪澇之災(zāi)后出現(xiàn)霍亂疫情,在多種環(huán)境要素的疊加作用下,上述地區(qū)更容易成為霍亂重災(zāi)區(qū).民國時(shí)期的瀘縣、富順縣和江北縣等霍亂疫災(zāi)高發(fā)區(qū)都存在以上特點(diǎn).地理探測(cè)器結(jié)果進(jìn)一步驗(yàn)證了自然和人文環(huán)境的交互作用對(duì)于霍亂流行具有強(qiáng)大的推動(dòng)作用.

表2 地理探測(cè)器運(yùn)算結(jié)果

3.2.2 霍亂驅(qū)動(dòng)因素的空間異質(zhì)性 通過相關(guān)性分析和逐步回歸分析,篩選出通過0.01顯著性檢驗(yàn)和共線性診斷的驅(qū)動(dòng)因素,即人口數(shù)量、洪澇災(zāi)害次數(shù),進(jìn)一步與霍亂流行年數(shù)在ArcGIS軟件中進(jìn)行地理加權(quán)回歸分析.從結(jié)果來看,人口數(shù)量對(duì)霍亂流行的影響呈正相關(guān),相對(duì)呈現(xiàn)緯向遞變,影響強(qiáng)度由四川省內(nèi)東南向西北遞減(圖6a).人口數(shù)量對(duì)霍亂流行影響作用較大的地區(qū)主要是川中及川東南地區(qū),這與民國時(shí)期四川省的人口分布大體吻合,再次驗(yàn)證人口集聚之地更容易促發(fā)疫病流行.川南的古宋、興文、筠連一帶,川西南的鹽邊、昭覺、雷波等地人口數(shù)量回歸系數(shù)較高,但其人口數(shù)量相對(duì)較少,可能是因?yàn)檫@些地區(qū)較為偏遠(yuǎn),民眾衛(wèi)生防疫觀念較為落后,一定程度上導(dǎo)致對(duì)霍亂流行更加敏感.洪澇災(zāi)害的影響同樣呈正相關(guān),強(qiáng)度具有從西北向東南遞增.洪澇災(zāi)害對(duì)霍亂流行的推動(dòng)作用在四川省東南部較為強(qiáng)烈(圖6b).結(jié)合主要河流的分布,川內(nèi)的渠江、嘉陵江、沱江、岷江等偏南北向河流皆向東南注入長江,在河流下游及交匯處容易爆發(fā)洪澇災(zāi)害,成為大規(guī)模霍亂流行的導(dǎo)火索.

圖6 民國時(shí)期四川省霍亂流行驅(qū)動(dòng)因素的空間分異情況Fig.6 Spatial differentiation of driving factors of cholera epidemic in Sichuan Province during the period of the Republic of China (1912-1949)

4 結(jié)論與討論

4.1 結(jié)論

1) 從時(shí)間規(guī)律來看,民國時(shí)期四川省的霍亂疫情在波動(dòng)中發(fā)展,霍亂流行頻度達(dá)97.50%,廣度達(dá)82.88%,出現(xiàn)1916年、1920年、1939年、1945年四個(gè)明顯高峰,20世紀(jì)30年代開始疫情加重,霍亂季發(fā)率為42.14%,高發(fā)于夏秋兩季.

2) 從空間規(guī)律來看,霍亂的疫災(zāi)面積總體呈現(xiàn)擴(kuò)大趨勢(shì),由早期的斑塊分散狀到中后期逐漸集中連片擴(kuò)展,高發(fā)地區(qū)不斷增多,四川盆地內(nèi)部的川中丘陵及成都平原的部分地區(qū)為霍亂主要流行區(qū),川西南及東南的部分縣則較少遭受霍亂侵襲.

3) 霍亂與自然環(huán)境和人文環(huán)境都具有較強(qiáng)的相關(guān)性.部分驅(qū)動(dòng)因素存在顯著的交互作用,人口數(shù)量、洪澇災(zāi)害因素具有明顯的空間異質(zhì)性.氣候、地形、水文、洪災(zāi)等自然環(huán)境因素和人口、交通、戰(zhàn)亂等人文環(huán)境因素相互交織,對(duì)霍亂的產(chǎn)生和發(fā)展具有強(qiáng)大推動(dòng).

4.2 討論

研究表明,霍亂流行的頻度與河道的分布間具有明顯的距離衰減性,體現(xiàn)出霍亂的水體傳播特性,與霍亂適生度分析結(jié)果相吻合[11],但河道長度與河網(wǎng)密度因素與霍亂流行的相關(guān)性并不顯著,還需進(jìn)一步驗(yàn)證.從相關(guān)性分析及空間耦合狀況來看,民國時(shí)期四川省霍亂多流行于低海拔及濕熱地帶,相關(guān)研究也發(fā)現(xiàn)低海拔地區(qū)往往有更多的霍亂病例,霍亂發(fā)病率隨降水和氣溫的升高而上升[19].洪澇災(zāi)害對(duì)霍亂流行的重要推動(dòng)作用,但旱災(zāi)與霍亂的關(guān)系還需進(jìn)一步研究,關(guān)于1932年陜西霍亂的研究[12]及相關(guān)史料表明極端干旱可能也與霍亂爆發(fā)有關(guān),如《樂至縣志》載:“(1920年)6—7月,天旱,駐軍遍駐民房,各街馬糞馬尿橫流,穢氣薰蒸,引起疫病,麻腳瘟流行.秋涼下雨,疫情緩解”[45],又如《江油縣志》載:“(1935年)冬,久旱無雨,井塘干涸,人畜多飲用臟水,霍亂等疾病流行,大堰、東興一帶發(fā)病甚多”[46].后續(xù)研究可綜合多源歷史環(huán)境數(shù)據(jù),深入探究霍亂流行的驅(qū)動(dòng)因素及機(jī)制.