人口流動(dòng)的影響機(jī)制研究
——基于氣候變化及區(qū)域差異

鄧宏乾 劉夢(mèng)曉 邵 勝

(1.華中師范大學(xué) 經(jīng)濟(jì)與工商管理學(xué)院， 湖北 武漢 430079； 2.華中師范大學(xué) 湖北房地產(chǎn)發(fā)展研究中心， 湖北 武漢 430079)

一、問(wèn)題提出

全球氣候變化近年來(lái)成為各界學(xué)者競(jìng)相關(guān)注的熱點(diǎn)，氣候條件也越來(lái)越被認(rèn)為是影響人類當(dāng)前活動(dòng)和未來(lái)決策的重要因素。人類對(duì)氣候變化的敏感程度可能超乎其自知水平，原生活環(huán)境細(xì)微的氣溫和降水量變化，會(huì)改變?nèi)藗兊男袨槠?，自?dòng)驅(qū)使人們“用腳投票”，流動(dòng)到其他氣候條件相對(duì)較好的地區(qū)。改革開放后，我國(guó)人口流動(dòng)規(guī)模急劇擴(kuò)張，據(jù)統(tǒng)計(jì)，截至2018年底，全國(guó)流動(dòng)人口數(shù)已達(dá)2.41億，這意味著全國(guó)17%①以上的人口都曾經(jīng)甚至依然在發(fā)生遷移。回顧歷史，中國(guó)人存在趨于安穩(wěn)、不愿意輕易背井離鄉(xiāng)的戀鄉(xiāng)情節(jié)，但目前我國(guó)人口流動(dòng)規(guī)模不斷擴(kuò)張，究竟是什么力量驅(qū)使人們離開自己的土地？以往學(xué)者對(duì)人口遷移或者人口流動(dòng)的因素分析，幾乎集中于城市經(jīng)濟(jì)發(fā)展水平和社會(huì)政治環(huán)境層面，鮮少探討氣候變化因素。鑒于全球性氣候問(wèn)題的日益突出，我國(guó)區(qū)域間氣候變化也日趨明顯，本文試圖探究氣候因素變化對(duì)中國(guó)人口流動(dòng)的影響機(jī)理，從氣候變化視角解釋人們“用腳投票”選擇區(qū)域流動(dòng)的原因。

隨著我國(guó)城鎮(zhèn)化加速發(fā)展，戶籍制度改革和農(nóng)業(yè)人口市民化政策的實(shí)施，破除了改革開放前農(nóng)業(yè)人口固化和“畫地為牢”體制，大量農(nóng)業(yè)人口逐步“市民化”，人口流動(dòng)呈現(xiàn)從農(nóng)村流向城市，從中西部向東部地區(qū)流動(dòng)的特點(diǎn)，人口流動(dòng)促進(jìn)了勞動(dòng)力資源區(qū)域間的優(yōu)化配置，促進(jìn)了經(jīng)濟(jì)發(fā)展。目前我國(guó)人口流動(dòng)以農(nóng)業(yè)人口為主，根據(jù)2016年全國(guó)人口衛(wèi)生計(jì)生動(dòng)態(tài)監(jiān)測(cè)調(diào)查數(shù)據(jù)統(tǒng)計(jì)，農(nóng)業(yè)流動(dòng)人口占總流動(dòng)人口比例高達(dá)83%，而非農(nóng)業(yè)人口占比僅為17%。人口流動(dòng)除了受經(jīng)濟(jì)、社會(huì)、政策的因素影響外，氣候自然條件成為人口遷移的一個(gè)重要原因，特別是破除城鄉(xiāng)二元戶籍制度，消除了人口跨區(qū)域、規(guī)模性流動(dòng)的體制障礙。氣溫和降水量是影響氣候的主要因素，氣候變化直接影響農(nóng)作物的產(chǎn)出，氣溫是土壤中有機(jī)碳變化的主要?dú)夂蛞蛩兀瑲鉁厣邥?huì)使得自然條件下土壤中有機(jī)碳含量減少，從而導(dǎo)致糧食作物和經(jīng)濟(jì)作物減產(chǎn)②。降水量是影響植被覆蓋變化的重要?dú)夂蛞蛩刂?，降水量的增加可以增?qiáng)植被修復(fù)生長(zhǎng)能力③，植被覆蓋面積的增加可以改善區(qū)域生態(tài)環(huán)境，增強(qiáng)地區(qū)吸引力和可持續(xù)發(fā)展能力，吸引人口的流入和聚集。氣溫升高和降水量減少將不利于地區(qū)農(nóng)作物的生長(zhǎng)④，糧食作物和經(jīng)濟(jì)作物產(chǎn)出的減少直接導(dǎo)致以農(nóng)業(yè)產(chǎn)出為重要經(jīng)濟(jì)來(lái)源的農(nóng)業(yè)人口經(jīng)濟(jì)收入下降，預(yù)期收入水平的下降成為農(nóng)業(yè)人口流出的重要誘因。農(nóng)業(yè)是我國(guó)經(jīng)濟(jì)發(fā)展的基礎(chǔ)產(chǎn)業(yè)，氣候條件的不利變化直接制約著地區(qū)經(jīng)濟(jì)發(fā)展水平和農(nóng)業(yè)人口的生存與發(fā)展，流動(dòng)到氣候條件更好、預(yù)期收入水平相對(duì)高的地區(qū)是農(nóng)業(yè)人口面對(duì)經(jīng)濟(jì)收入風(fēng)險(xiǎn)的普遍應(yīng)對(duì)方法?；谝陨?，本文提出兩個(gè)假設(shè)：

假設(shè)H1:氣溫變化與人口流動(dòng)間存在顯著線性負(fù)相關(guān)關(guān)系，且存在區(qū)域化差異。

假設(shè)H2:降水量變化與人口流動(dòng)間存在顯著線性正相關(guān)關(guān)系，且存在區(qū)域化差異。

二、文獻(xiàn)綜述

氣候變化對(duì)人類行為產(chǎn)生影響，這一觀點(diǎn)已得到多方學(xué)者的認(rèn)證。早期關(guān)于氣候變化對(duì)人類活動(dòng)的研究主要是自然災(zāi)害和極端氣候變化引發(fā)的人口遷徙，如地震、海嘯、颶風(fēng)、山洪、旱災(zāi)、雹災(zāi)、火災(zāi)、核輻射等⑤。Hinnawi因此提出“環(huán)境難民”的概念，即由于自然原因?qū)ι鷳B(tài)環(huán)境造成嚴(yán)重破壞，嚴(yán)重危及人們的生存環(huán)境，從而人們被迫永久性離開原住地⑥；而后Swain提出“環(huán)境移民”的概念，弱化了“環(huán)境難民”的政治屬性，強(qiáng)化了氣候環(huán)境變化對(duì)人口流動(dòng)的影響結(jié)果⑦。國(guó)內(nèi)學(xué)者首先提出了“災(zāi)害移民”概念，即自然災(zāi)害引發(fā)的人口流動(dòng)⑧；繼而提出“生態(tài)移民”，即生態(tài)環(huán)境變化引發(fā)的人口流動(dòng)⑨；目前學(xué)界開始研究“氣候移民”問(wèn)題，即氣候變化引發(fā)的人口流動(dòng)⑩，相對(duì)于早期的極端天氣變化而引起人口遷移，緩慢、漸變式的氣候變化對(duì)人類行為的影響是不同的，因此氣候移民逐漸成為近來(lái)研究的重點(diǎn)，因?yàn)檫@類氣候變化更為普遍，如空氣質(zhì)量變化對(duì)來(lái)華留學(xué)生城市選擇的影響；氣溫變化對(duì)不同收入水平國(guó)家人口流動(dòng)的影響；降水量變化對(duì)國(guó)際人口流動(dòng)的影響等。

已有文獻(xiàn)表明氣溫變化與人口流動(dòng)存在緊密聯(lián)系，但學(xué)者們的研究結(jié)果存在差異。絕大部分學(xué)者認(rèn)為氣溫變化與人口流動(dòng)存在線性關(guān)系。董文麗、李王鳴指出，嚴(yán)寒地區(qū)的人們更有向城市地區(qū)流動(dòng)的趨勢(shì)，即氣溫條件好的城市地區(qū)吸引著更多的人口流動(dòng)。盧洪友等研究表明，氣溫變化與人口流動(dòng)成線性正相關(guān)關(guān)系，氣溫每上升1攝氏度，流動(dòng)人口比率平均上漲0.46個(gè)百分點(diǎn)。類似地，Mueller等在研究非洲東部人口流動(dòng)問(wèn)題時(shí)發(fā)現(xiàn)，氣溫偏離平均氣溫每增加1個(gè)標(biāo)準(zhǔn)單位，凈遷出率平均下降10個(gè)百分點(diǎn)。部分學(xué)者則認(rèn)為氣溫變化與人口流動(dòng)存在著線性負(fù)相關(guān)關(guān)系。王春超、王志文研究結(jié)果表明，氣溫變化與城市人口流動(dòng)的關(guān)系取決于城市是否為糧食主產(chǎn)區(qū)。其進(jìn)一步指出，糧食主產(chǎn)區(qū)的氣溫變化與流動(dòng)人口比率呈顯著負(fù)相關(guān)關(guān)系，而在非糧食主產(chǎn)區(qū)呈顯著正相關(guān)關(guān)系。這一結(jié)論與Cai等研究結(jié)果基本吻合，在農(nóng)業(yè)為主的國(guó)家，氣溫上升會(huì)導(dǎo)致人口流出，流動(dòng)人口比率下降。還有部分外國(guó)學(xué)者認(rèn)為氣溫變化對(duì)人口流動(dòng)的影響不是簡(jiǎn)單的線性關(guān)系，而是非線性的，如Bohra-Mishra等發(fā)現(xiàn)在印度尼西亞，氣溫變化與人口流動(dòng)間存在非線性負(fù)相關(guān)關(guān)系，當(dāng)氣溫超過(guò)25攝氏度時(shí)，氣溫的上升會(huì)導(dǎo)致人口流出的顯著增加。

降水量變化對(duì)人口流動(dòng)也存在重要影響。部分學(xué)者認(rèn)為降水量與人口流動(dòng)緊密相關(guān)，其中部分學(xué)者認(rèn)為降水量變化與人口流動(dòng)呈正相關(guān)關(guān)系，如溫小潔等指出，降水條件良好的地區(qū)植被生長(zhǎng)及修復(fù)能力較強(qiáng)，會(huì)促進(jìn)人口聚集；Abel等則認(rèn)為降水量的增加會(huì)減少社會(huì)沖突，減少城市人口的流出。但部分學(xué)者卻得出降水量變化與人口流動(dòng)呈負(fù)相關(guān)關(guān)系的結(jié)論，如Barrios等發(fā)現(xiàn)撒哈拉以南的非洲，降水量的增加將不利于人口向城市的流入；Mueller等同樣發(fā)現(xiàn)，在非洲東部降水量平均每增加1單位，城市凈遷出率平均上漲10個(gè)百分點(diǎn)。此外，降水量變化與人口流動(dòng)間非線性關(guān)系也被學(xué)者提出并印證。鄢姣和胡必亮在研究降雨量與城鎮(zhèn)化率的關(guān)系中發(fā)現(xiàn)，降雨量變化與人口流動(dòng)間并非是簡(jiǎn)單的線性關(guān)系，而是U型非線性關(guān)系：隨著降雨量的增加，城鎮(zhèn)人口流動(dòng)呈現(xiàn)先減后增的態(tài)勢(shì)。

國(guó)內(nèi)外學(xué)者關(guān)于氣溫變化和降水量變化對(duì)人口流動(dòng)影響得出的結(jié)論之所以存在差異，原因在于學(xué)者們選取的研究對(duì)象所在國(guó)別和地區(qū)初始?xì)夂驐l件的差異較大以及氣候變化對(duì)人口流動(dòng)的影響尚缺乏理論依據(jù)。學(xué)界同時(shí)研究氣溫及降水量變化對(duì)人口流動(dòng)影響的文獻(xiàn)相對(duì)較少。本文的邊際貢獻(xiàn)是，利用全國(guó)省會(huì)城市和計(jì)劃單列市數(shù)據(jù)，同時(shí)考察氣溫變化及降水量變化對(duì)人口流動(dòng)的影響，采用氣溫和降水量的相對(duì)變化量而非絕對(duì)量作為核心解釋變量，更能客觀地描述氣溫變化和降水量變化對(duì)人口流動(dòng)的影響及其區(qū)域差異。

三、數(shù)據(jù)來(lái)源和模型設(shè)定

1.數(shù)據(jù)來(lái)源

本文研究2000-2018年省會(huì)城市和計(jì)劃單列市氣溫降水量變化對(duì)人口流動(dòng)的影響。數(shù)據(jù)來(lái)源于2000-2019年度的《中國(guó)城市統(tǒng)計(jì)年鑒》、各省會(huì)城市及計(jì)劃單列市統(tǒng)計(jì)年鑒、城市國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)、國(guó)家統(tǒng)計(jì)局和全國(guó)溫室數(shù)據(jù)系統(tǒng)。由于個(gè)別城市相關(guān)資料較多缺失，本文剔除了拉薩、蘭州、南寧這3個(gè)省會(huì)城市；此外對(duì)于個(gè)別城市年份數(shù)值的缺失，采用插值法補(bǔ)齊。

2.變量選取與設(shè)定

根據(jù)研究?jī)?nèi)容并參考以往文獻(xiàn)，本文選取的變量和測(cè)度方法如下(參見(jiàn)表1)：

測(cè)度人口流動(dòng)的被解釋變量為流動(dòng)人口比率變化值，即當(dāng)年流動(dòng)人口比率-上年流動(dòng)人口比率。由于流動(dòng)人口數(shù)據(jù)的缺乏和各年數(shù)據(jù)統(tǒng)計(jì)口徑的差異，本文將采取統(tǒng)一計(jì)算方法進(jìn)行量化。根據(jù)各城市相應(yīng)年份統(tǒng)計(jì)年鑒的人口數(shù)據(jù)，設(shè)計(jì)如下計(jì)算方法：流動(dòng)人口數(shù)=常住人口數(shù)-戶籍人口數(shù)，流動(dòng)人口比率=流動(dòng)人口數(shù)/常住人口數(shù)×100%，流動(dòng)人口比率變化值=當(dāng)年流動(dòng)人口比率-上年流動(dòng)人口比率。對(duì)于部分城市或年份年鑒中常住人口或戶籍人口數(shù)據(jù)缺失，但卻報(bào)告了流入人口和流出人口數(shù)的，采用如下測(cè)度方法：常住人口數(shù)=戶籍人口數(shù)+流入人口數(shù)-流出人口數(shù)；流動(dòng)人口數(shù)=流入人口數(shù)-流出人口數(shù)=凈流入人口數(shù)，流動(dòng)人口比率=凈流入人口數(shù)/(戶籍人口數(shù)+凈流入人口數(shù))。此外，在穩(wěn)健性檢驗(yàn)中，將進(jìn)一步引入被解釋變量的滯后一期數(shù)據(jù)。

測(cè)度氣候變化的核心解釋變量為氣溫變化和降水量變化。氣溫變化采用年平均氣溫變化值表示：年平均氣溫變化值=當(dāng)年年平均氣溫-上年年平均氣溫，年平均氣溫根據(jù)相應(yīng)年份日度平均氣溫計(jì)算得出。此外本文也采用月平均氣溫變化值做穩(wěn)健性檢驗(yàn)，月平均氣溫根據(jù)相應(yīng)年份日度平均氣溫計(jì)算得出。

降水量變化采用年總降水量變化值和年總降水量變化率表示：年總降水量變化值=當(dāng)年年總降水量-上年年總降水量，年總降水量變化率=(當(dāng)年年總降水量-上年年總降水量)/上年年總降水量，年總降水量根據(jù)相應(yīng)年份日度降水量計(jì)算得出。另外本文也采用月平均降水量變化值做穩(wěn)健性檢驗(yàn)，月平均降水量根據(jù)相應(yīng)年份日度降水量計(jì)算得出。

為了減少遺漏解釋變量帶來(lái)的估計(jì)偏誤，將引入如下與城市人口流動(dòng)相關(guān)的控制變量：

經(jīng)濟(jì)發(fā)展水平。城市的經(jīng)濟(jì)發(fā)展水平無(wú)疑是吸引人口流入的重要因素，本文采用人均地區(qū)生產(chǎn)總值來(lái)表示這一指標(biāo)，并設(shè)定該指標(biāo)的對(duì)數(shù)形式表示該變量。

城市居民收入水平。本文采用城鎮(zhèn)居民人均可支配收入來(lái)表示城市的平均收入水平，并設(shè)定城鎮(zhèn)居民人均可支配收入的對(duì)數(shù)形式來(lái)代表該變量。

表1 變量描述性統(tǒng)計(jì)

城市留居難易程度。人們考慮遷移居住地的前提，首先考慮的是住房可得性。能否負(fù)擔(dān)購(gòu)房壓力從而長(zhǎng)期定居，或者承擔(dān)租房租金從而長(zhǎng)期留居都會(huì)成為人們選擇是否流入該城市的考量指標(biāo)，而城市房?jī)r(jià)可以很好地衡量這一指標(biāo)。本文采取各城市住宅商品房平均銷售價(jià)格來(lái)表示這一指標(biāo)，并設(shè)定其對(duì)數(shù)形式來(lái)代表該變量。

城市公共服務(wù)水平。城市的教育資源、醫(yī)療水平、基礎(chǔ)交通設(shè)施是城市的競(jìng)爭(zhēng)優(yōu)勢(shì)，城市的公共服務(wù)水平通常與其財(cái)政收支緊密相關(guān)。城市公共服務(wù)水平的高低影響著人口流動(dòng)，一般來(lái)說(shuō)，公共服務(wù)水平較高的城市也吸引著更多的人口流入，本文采用人均地方一般公共預(yù)算收入來(lái)表示這一指標(biāo)。鑒于沒(méi)有統(tǒng)計(jì)年鑒公布這一指標(biāo)，本文采用如下方式計(jì)算得出：人均地方一般公共預(yù)算收入=地方一般公共預(yù)算收入/年末總?cè)丝跀?shù)。類似地，設(shè)定人均地方一般公共預(yù)算收入的對(duì)數(shù)形式代表該變量。

城市產(chǎn)業(yè)結(jié)構(gòu)。隨著城市化進(jìn)程的深入，產(chǎn)業(yè)結(jié)構(gòu)也發(fā)生著顯著的變化。第三產(chǎn)業(yè)范圍涵蓋廣泛，崗位豐富，可吸納更多的流動(dòng)人口。本文采用第三產(chǎn)業(yè)占GDP比重代表該指標(biāo)。

城市就業(yè)難易程度。就業(yè)問(wèn)題是關(guān)乎民生和社會(huì)穩(wěn)定的大問(wèn)題，城市潛在的勞動(dòng)崗位和就業(yè)的難易程度都成為人口流入的考慮因素。失業(yè)率偏高的城市，通常說(shuō)明潛在空余勞動(dòng)力需求不太旺盛，增大了新進(jìn)城市人口就業(yè)生存的難度。本文采用城鎮(zhèn)登記失業(yè)率來(lái)代表城市就業(yè)難易程度。

3.計(jì)量模型

基于上文中所提出的兩個(gè)非線性關(guān)系假設(shè)，區(qū)分氣溫和降水量設(shè)如下三個(gè)計(jì)量方程：

(1)

(2)

(3)

四、實(shí)證結(jié)果

Hausman檢驗(yàn)結(jié)果顯示，氣溫、降水量和聯(lián)合方程均適用固定效應(yīng)模型，本文對(duì)上述3個(gè)方程分別進(jìn)行基準(zhǔn)回歸分析。

1.氣溫基準(zhǔn)方程

基準(zhǔn)氣溫方程中的氣溫變化選取的測(cè)度指標(biāo)為年平均氣溫變化值(ctem)，固定效應(yīng)回歸結(jié)果如表2所示。列(1)(2)僅考慮氣溫變化一次項(xiàng)與人口流動(dòng)的關(guān)系，結(jié)果顯示年平均氣溫變化值一次項(xiàng)系數(shù)在5%水平上顯著為負(fù)。列(3)(4)中繼續(xù)引入年平均氣溫變化值的二次項(xiàng)以考察氣溫變化對(duì)人口流動(dòng)的非線性關(guān)系，結(jié)果顯示氣溫變化值的二次項(xiàng)系數(shù)并不顯著，而年平均氣溫變化值的一次項(xiàng)系數(shù)依然在5%水平上顯著為負(fù)。氣溫基準(zhǔn)方程的回歸結(jié)果表明，年平均氣溫變化值與流動(dòng)人口比率變化值間存在顯著線性負(fù)相關(guān)關(guān)系，即氣溫變化與人口流動(dòng)間存在顯著線性負(fù)相關(guān)關(guān)系，假設(shè)H1成立。且回歸結(jié)果顯示年平均氣溫變化值系數(shù)穩(wěn)定在-0.17左右，表明當(dāng)期年平均氣溫每升高1攝氏度，當(dāng)期流動(dòng)人口比率較上期平均會(huì)下降0.17個(gè)百分點(diǎn)。氣溫升高，流動(dòng)人口比率趨于下降，人口流出趨于增加，人口流入趨于減少。

表2 氣溫基準(zhǔn)方程回歸結(jié)果

續(xù)表

2.降水量基準(zhǔn)方程

降水量基準(zhǔn)方程中的降水量變化采用年總降水量變化值(ctpre)測(cè)度，固定效應(yīng)回歸結(jié)果如表3所示。列(5)(6)僅考慮降水量變化與人口流動(dòng)間的線性關(guān)系，結(jié)果顯示，年總降水量變化值的系數(shù)在5%水平上顯著為正，表明降水量變化與人口流動(dòng)間存在顯著線性正相關(guān)關(guān)系。列(7)(8)繼續(xù)加入年總降水量變化值的二次項(xiàng)，以探究降水量變化對(duì)人口流動(dòng)的非線性影響，結(jié)果顯示年總降水量變化值二次項(xiàng)系數(shù)均不顯著，而其一次項(xiàng)系數(shù)依然維持在5%水平上顯著為正，這表明降水量變化與人口流動(dòng)間僅存在線性正相關(guān)關(guān)系，假設(shè)H2成立。年總降水量變化值系數(shù)穩(wěn)定在0.0007左右，表明年總降水量每增加1000毫米，當(dāng)期流動(dòng)人口比率較上期平均上漲0.7個(gè)百分點(diǎn)。年總降水量增加，意味著該地區(qū)氣候條件有所改善，地區(qū)競(jìng)爭(zhēng)力吸引力提高，人口流入量預(yù)期會(huì)增加，當(dāng)期流動(dòng)人口比率也趨于上升。

3.聯(lián)合基準(zhǔn)方程

聯(lián)合基準(zhǔn)方程旨在探究氣溫變化和降水量變化聯(lián)合對(duì)人口流動(dòng)的影響，氣溫變化和降水量變化選取的測(cè)度指標(biāo)分別為年平均氣溫變化值(ctem)和年總降水量變化值(ctpre)，回歸結(jié)果見(jiàn)表4。列(9)(10)顯示年平均氣溫變化值一次項(xiàng)系數(shù)在10%水平上顯著為負(fù)，年總降水量變化值一次項(xiàng)系數(shù)在5%水平上顯著為正。列(11)(12)加入了氣溫變化和降水量變化的二次項(xiàng)，結(jié)果顯示年平均氣溫變化值二次項(xiàng)系數(shù)不顯著，其一次項(xiàng)系數(shù)維持在10%水平上顯著為負(fù)；年總降水量變化值二次項(xiàng)系數(shù)亦不顯著，其一次項(xiàng)系數(shù)依然在5%水平上顯著為正。聯(lián)合基準(zhǔn)方程回歸結(jié)果中核心解釋變量的符號(hào)和顯著性水平與氣溫、降水量基準(zhǔn)方程所得結(jié)果基本一致，只是系數(shù)大小稍有變化，這表明氣溫變化與流動(dòng)人口間存在顯著線性負(fù)相關(guān)關(guān)系，降水量變化與人口流動(dòng)間存在顯著線性正相關(guān)關(guān)系，假設(shè)H1、H2依然成立。

表3 降水量基準(zhǔn)方程回歸結(jié)果

五、穩(wěn)健性檢驗(yàn)

1.氣溫方程穩(wěn)健性檢驗(yàn)

為了檢驗(yàn)氣溫變化對(duì)人口流動(dòng)影響結(jié)果的穩(wěn)健性，從以下幾方面進(jìn)行檢驗(yàn)討論：

(1)在氣溫基準(zhǔn)方程中引入被解釋變量滯后一期項(xiàng)(L.cmig)，采用差分GMM方法進(jìn)行動(dòng)態(tài)面板回歸(見(jiàn)表5)。列(13)(14)僅考慮一階自相關(guān)基礎(chǔ)上年平均氣溫變化值與人口流動(dòng)的線性關(guān)系，列(15)(16)在此基礎(chǔ)上加入年平均氣溫變化值的二次項(xiàng)，以探究氣溫變化與人口流動(dòng)的非線性關(guān)系。結(jié)果顯示年平均氣溫變化值的一次項(xiàng)系數(shù)穩(wěn)定在10%水平上顯著為負(fù)，而二次項(xiàng)系數(shù)均不顯著，表明氣溫變化與人口流動(dòng)間存在顯著線性負(fù)相關(guān)關(guān)系，這與氣溫基準(zhǔn)方程得到的結(jié)論基本一致，假設(shè)H1得到檢驗(yàn)。此外，表5結(jié)果顯示當(dāng)期流動(dòng)人口比率變化值與滯后一期流動(dòng)人口比率變化值均呈顯著負(fù)相關(guān)關(guān)系，若上期流動(dòng)人口比率增長(zhǎng)，那么當(dāng)期流動(dòng)人口比率則會(huì)出現(xiàn)小幅度下降；若上期流動(dòng)人口比率下降，那么當(dāng)期流動(dòng)人口比率則會(huì)出現(xiàn)小幅度上漲，這種變動(dòng)關(guān)系會(huì)有助于流動(dòng)人口比率維持在一個(gè)相對(duì)穩(wěn)定的增長(zhǎng)狀態(tài)，不會(huì)出現(xiàn)單方向的爆炸式增長(zhǎng)或下跌。

表5 氣溫方程穩(wěn)健性檢驗(yàn)結(jié)果

(2)將氣溫方程中的核心解釋變量替換為月平均氣溫變化值(cmtem)再次進(jìn)行穩(wěn)健性檢驗(yàn)(見(jiàn)表6)，前4列結(jié)果顯示月平均氣溫變化值一次項(xiàng)系數(shù)在5%水平上顯著為負(fù)，其二次項(xiàng)系數(shù)均不顯著；表6中后4列繼續(xù)引入被解釋變量滯后一期(L.cmig)進(jìn)行動(dòng)態(tài)面板回歸，結(jié)果顯示月平均氣溫變化值一次項(xiàng)系數(shù)依然在10%水平上顯著為負(fù)，系數(shù)大小略有降低，而其二次項(xiàng)系數(shù)依然不顯著。這表明月平均氣溫變化值與流動(dòng)人口比率變化值間存在顯著線性負(fù)相關(guān)關(guān)系，假設(shè)H1得到檢驗(yàn)。且月平均氣溫變化值系數(shù)約為-2.0，表示月平均氣溫每升高1攝氏度，當(dāng)期流動(dòng)人口比率較上期平均會(huì)下降2個(gè)百分點(diǎn)。這與氣溫基準(zhǔn)方程中年平均氣溫每升高1攝氏度，當(dāng)期流動(dòng)人口比率較上期平均會(huì)下降0.17個(gè)百分點(diǎn)的結(jié)論基本吻合。

表6 氣溫方程穩(wěn)健性檢驗(yàn)結(jié)果

(3)分區(qū)域回歸。根據(jù)城市所處地區(qū)將該樣本劃分為中、東、西和東北部四個(gè)地區(qū)進(jìn)行子樣本分類回歸，用費(fèi)舍爾組合檢驗(yàn)進(jìn)行組間系數(shù)差異檢驗(yàn)以考察氣溫變化對(duì)人口流動(dòng)的影響是否存在區(qū)域化差異，氣溫變化測(cè)度指標(biāo)為年平均氣溫變化值(ctem)。氣溫基準(zhǔn)方程分區(qū)域回歸結(jié)果(見(jiàn)表7)顯示年平均氣溫變化值二次項(xiàng)系數(shù)均不顯著，東部和中部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)分別在10%和5%水平上顯著為負(fù)，而西部和東北部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)均不顯著，這表明在東部和中部地區(qū)氣溫變化與人口流動(dòng)間的線性負(fù)相關(guān)關(guān)系顯著不同于西部和東北部，氣溫變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異，假設(shè)H1成立。東部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)約為-0.35，表明東部地區(qū)年平均氣溫每升高1攝氏度，當(dāng)期流動(dòng)人口比率較上期平均會(huì)下降0.35個(gè)百分點(diǎn)，下降幅度高于全國(guó)平均水平0.17個(gè)百分點(diǎn)。中部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)約為-0.28，表明中部地區(qū)年平均氣溫每升高1攝氏度，當(dāng)期流動(dòng)人口比率較上期平均會(huì)下降0.28個(gè)百分點(diǎn)，下降幅度高于全國(guó)平均水平0.17個(gè)百分點(diǎn)。東部和中部地區(qū)人口流動(dòng)對(duì)于氣溫變化較全國(guó)其他地區(qū)更為敏感，這是由于東部和中部地區(qū)的年平均氣溫相較于西部和東北部而言整體偏高，若年平均氣溫繼續(xù)升高會(huì)抑制人口流入，從而導(dǎo)致流動(dòng)人口比率下降。

表7 氣溫基準(zhǔn)方程分區(qū)域回歸結(jié)果

2.降水量方程穩(wěn)健性檢驗(yàn)

為了進(jìn)一步驗(yàn)證降水量變化與人口流動(dòng)間的線性正相關(guān)關(guān)系，采用如下方法對(duì)基準(zhǔn)方程結(jié)果進(jìn)行穩(wěn)健性檢驗(yàn)：

(1)將被解釋變量滯后一期(L.cmig)加入降水量基準(zhǔn)方程中，采用差分GMM法進(jìn)行動(dòng)態(tài)面板回歸，計(jì)量結(jié)果見(jiàn)表8。列(33)(34)結(jié)果顯示加入滯后一期被解釋變量后，年總降水量變化值系數(shù)依然在5%水平上顯著為正，表明降水量變化對(duì)人口流動(dòng)的顯著線性影響依然穩(wěn)健。列(35)(36)結(jié)果顯示加入滯后一期被解釋變量后，年總降水量變化值二次項(xiàng)系數(shù)均不顯著，而一次項(xiàng)系數(shù)依然在5%水平上顯著為正，這表明降水量變化與人口流動(dòng)間存在線性正相關(guān)關(guān)系，假設(shè)H2得到檢驗(yàn)。此外，加入之后一期被解釋變量后，年總降水量變化值系數(shù)穩(wěn)定在0.0007左右，這與降水基準(zhǔn)方程的系數(shù)基本保持一致。

(2)將降水量變化指標(biāo)更換為年總降水量變化率(rtpre)進(jìn)行檢驗(yàn)，列(37)(38)顯示年總降水量變化率的一次項(xiàng)系數(shù)在5%水平上顯著為正，表明年總降水量變化率會(huì)對(duì)流動(dòng)人口比率變化值產(chǎn)生顯著正向影響，即降水量變化會(huì)對(duì)人口流動(dòng)產(chǎn)生顯著線性正影響，假設(shè)H2得到檢驗(yàn)。列(39)(40)引入年總降水量變化率的二次項(xiàng)，結(jié)果顯示其系數(shù)均不顯著，而年總降水量變化率的一次項(xiàng)系數(shù)在10%水平上顯著為正，且年總降水量變化率系數(shù)穩(wěn)定在0.005左右，這表明當(dāng)期年總降水量較上期每增加1%，則當(dāng)期流動(dòng)人口比率較上期平均會(huì)上漲0.005個(gè)百分點(diǎn)。

表8 降水量方程穩(wěn)健性檢驗(yàn)結(jié)果

(3)將降水量方程中降水量變化指標(biāo)更換為月平均降水量變化值(cmpre)進(jìn)一步進(jìn)行穩(wěn)健性檢驗(yàn)，回歸結(jié)果見(jiàn)表9。列(41)(42)顯示月平均降水量變化值系數(shù)在5%水平上顯著為正；列(45)(46)表示加入滯后一期被解釋變量進(jìn)行動(dòng)態(tài)回歸，結(jié)果顯示月平均降水量變化值系數(shù)依然維持在5%水平上顯著為正。列(43)(44)顯示月平均降水量變化值二次項(xiàng)系數(shù)均不顯著，而一次項(xiàng)系數(shù)依然維持在5%顯著性水平；列(47)(48)顯示加入滯后一期被解釋變量后，月平均降水量變化值二次項(xiàng)系數(shù)依然不顯著，而其一次項(xiàng)系數(shù)始終維持在5%水平上顯著為正。這表明月平均降水量變化值與流動(dòng)人口比率變化值存在顯著線性正相關(guān)關(guān)系，假設(shè)H2再次得到檢驗(yàn)。且月平均降水量變化值一次項(xiàng)系數(shù)穩(wěn)定在0.008左右，表示月平均降水量每增加100毫米，當(dāng)期流動(dòng)人口比率較上期平均上漲0.8個(gè)百分點(diǎn)。這與降水基準(zhǔn)方程中的年總降水量每增加1000毫米，當(dāng)期流動(dòng)人口比率較上期平均上漲0.7個(gè)百分點(diǎn)的結(jié)論基本吻合。

表9 降水量方程穩(wěn)健性檢驗(yàn)結(jié)果

(4)分區(qū)域進(jìn)行回歸。類似地，將該樣本根據(jù)中、東、西和東北部地區(qū)分為四個(gè)子樣本，進(jìn)行組間系數(shù)差異檢驗(yàn)，降水量變化選取的指標(biāo)為年總降水量變化值(ctpre)。分區(qū)域回歸結(jié)果顯示東部地區(qū)年總降水量變化值一次項(xiàng)系數(shù)在1%水平上顯著為正，中、西、東北部地區(qū)的年總降水量變化值一次項(xiàng)系數(shù)均不顯著，而所有地區(qū)年總降水量變化值二次項(xiàng)系數(shù)均不顯著(見(jiàn)表10)，這表明在東部地區(qū)，年總降水量變化值對(duì)流動(dòng)人口比率變化值的線性影響程度顯著不同于中、西、東北部地區(qū)，降水量變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異，假設(shè)H2成立。東部地區(qū)年總降水量變化值系數(shù)約為0.0009，表示東部地區(qū)年總降水量每增加1000毫米，當(dāng)期流動(dòng)人口比率較上期平均會(huì)上漲0.9個(gè)百分點(diǎn)，略高全國(guó)其他地區(qū)平均上漲0.7個(gè)百分點(diǎn)的水平。

3.聯(lián)合方程穩(wěn)健性檢驗(yàn)

為了進(jìn)一步驗(yàn)證氣溫降水量的聯(lián)合變化對(duì)人口流動(dòng)的影響，采用如下方法進(jìn)行穩(wěn)健性檢驗(yàn)：

(1)更換核心解釋變量測(cè)度指標(biāo)，將聯(lián)合基準(zhǔn)方程中氣溫變化指標(biāo)替換為月平均氣溫變化值(cmtem)，降水量變化指標(biāo)替換為月平均降水量變化值(cmpre)。聯(lián)合方程穩(wěn)健性結(jié)果(見(jiàn)表11)顯示月平均氣溫變化值和月平均降水量變化值的二次項(xiàng)系數(shù)均不顯著，而月平均氣溫變化值一次項(xiàng)系數(shù)始終維持在10%水平上顯著為負(fù)，月平均降水量變化值一次項(xiàng)始終維持在5%水平上顯著為正，這表示氣溫變化與人口流動(dòng)間存在顯著線性負(fù)相關(guān)關(guān)系，降水量變化與人口流動(dòng)間存在顯著線性正相關(guān)關(guān)系，假設(shè)H1、H2得到檢驗(yàn)。

表10 降水量基準(zhǔn)方程分區(qū)域回歸結(jié)果

表11 聯(lián)合方程的穩(wěn)健性檢驗(yàn)結(jié)果

續(xù)表

(2)分區(qū)域進(jìn)行回歸。將該樣本根據(jù)中、東、西和東北部地區(qū)分為四個(gè)子樣本進(jìn)行組間系數(shù)差異檢驗(yàn)，選取的氣溫變化測(cè)度指標(biāo)為年平均氣溫變化值(ctem)，降水量變化測(cè)度指標(biāo)為年總降水量變化值(ctpre)。聯(lián)合方程分區(qū)域回歸結(jié)果(見(jiàn)表12)顯示東部和中部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)均在10%水平上顯著為負(fù)，西部和東北部地區(qū)系數(shù)則不顯著，這表示東部和中部地區(qū)的氣溫變化與人口流動(dòng)間線性負(fù)相關(guān)關(guān)系顯著不同于西部和東北部地區(qū)，氣溫變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異，假設(shè)H1得到檢驗(yàn)。且聯(lián)合方程中東部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)約為-0.25，中部地區(qū)年平均氣溫變化值一次項(xiàng)系數(shù)約為-0.27，與氣溫方程的分區(qū)域回歸結(jié)果基本接近。東部地區(qū)年總降水量變化值一次項(xiàng)系數(shù)在5%水平上顯著為正，中、西、東北部地區(qū)系數(shù)均不顯著，且所有地區(qū)年總降水量變化值的二次項(xiàng)系數(shù)均不顯著，這表明東部地區(qū)的降水量變化與人口流動(dòng)間的線性正相關(guān)關(guān)系顯著不同于其他地區(qū)，降水量變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異，假設(shè)H2得到檢驗(yàn)。且聯(lián)合方程中東部地區(qū)年總降水量變化值系數(shù)與降水量方程分區(qū)域回歸結(jié)果中的系數(shù)較接近。由此可知，聯(lián)合方程的分區(qū)域回歸結(jié)果與氣溫和降水量的分區(qū)域回歸結(jié)果基本吻合，假設(shè)H1、H2得到檢驗(yàn)。

表12 聯(lián)合基準(zhǔn)方程分區(qū)域回歸結(jié)果

六、結(jié)論與建議

通過(guò)構(gòu)建人口流動(dòng)的省級(jí)和計(jì)劃單列市面板模型，本文研究了氣溫變化和降水量變化對(duì)人口流動(dòng)的影響機(jī)制。結(jié)果表明，氣溫變化與人口流動(dòng)存在顯著線性負(fù)相關(guān)關(guān)系，氣溫升高，流動(dòng)人口比率趨于下降，氣溫升高對(duì)人口流動(dòng)有明顯的抑制作用；且東部和中部地區(qū)氣溫升高對(duì)人口流動(dòng)的抑制作用顯著高于西部和東北部地區(qū)，氣溫變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異。降水量變化與人口流動(dòng)存在顯著線性正相關(guān)關(guān)系，降水量增加，流動(dòng)人口比率趨于上升，降水量增加對(duì)人口流動(dòng)存在顯著促進(jìn)作用；且東部地區(qū)降水量增加對(duì)人口流動(dòng)的促進(jìn)作用顯著高于其他地區(qū)，降水量變化對(duì)人口流動(dòng)的影響存在區(qū)域化差異?；诒疚牡贸龅臍鉁刈兓徒邓孔兓瘜?duì)人口流動(dòng)影響的結(jié)論，本文提出以下三點(diǎn)建議：

第一，社會(huì)經(jīng)濟(jì)活動(dòng)排放溫室氣體導(dǎo)致大氣中溫室氣體濃度增加，引起了以氣溫升高為主要特征的氣候變化，城市人口和產(chǎn)業(yè)聚集加劇了溫室氣體的排放，導(dǎo)致城市局部溫度的明顯升高。為減緩氣溫升高對(duì)人口流動(dòng)的抑制作用，應(yīng)建立覆蓋所有地區(qū)的溫室氣體排放監(jiān)測(cè)和統(tǒng)計(jì)制度，因地制宜地建立溫室氣體排放管控和預(yù)警機(jī)制?；跂|部和中部地區(qū)氣溫升高對(duì)人口流動(dòng)抑制作用高于其他地區(qū)的結(jié)論，東部和中部地區(qū)城市更需要強(qiáng)化溫室氣體排放控制，實(shí)施低碳綠色經(jīng)濟(jì)發(fā)展戰(zhàn)略，平衡好經(jīng)濟(jì)發(fā)展速度與環(huán)境質(zhì)量改善，從而推動(dòng)經(jīng)濟(jì)高質(zhì)量發(fā)展。

第二，中國(guó)是嚴(yán)重缺水的國(guó)家之一，不少地區(qū)存在著資源性、工程性、水質(zhì)性等缺水問(wèn)題。中國(guó)西部地區(qū)荒漠化引發(fā)的缺水問(wèn)題較為突出，降水量的不足抑制了西部地區(qū)的人口流入，嚴(yán)重制約著西部地區(qū)經(jīng)濟(jì)發(fā)展。西部地區(qū)應(yīng)著力提高沙漠科學(xué)治理水平，選取耐寒、耐旱、耐鹽堿的植物種子開展人工造林種草，建設(shè)綠色沙漠屏障，改善土壤吸水蓄水能力，樹立綠色發(fā)展理念，增強(qiáng)西部地區(qū)人口流入意愿，探索西部地區(qū)“生態(tài)和扶貧”雙驅(qū)動(dòng)的可持續(xù)發(fā)展之路。同時(shí)東部地區(qū)也應(yīng)做好生態(tài)保護(hù)和修復(fù)工作，避免以犧牲生態(tài)環(huán)境為代價(jià)的經(jīng)濟(jì)發(fā)展。

第三，加大對(duì)西部地區(qū)生態(tài)環(huán)境治理的政策、資金支持，以減輕不利氣候變化對(duì)西部地區(qū)人口流動(dòng)的抑制作用，以防止大量“環(huán)境移民”現(xiàn)象，創(chuàng)造條件承接?xùn)|、中部地區(qū)產(chǎn)業(yè)轉(zhuǎn)移，大力發(fā)展勞動(dòng)密集型產(chǎn)業(yè)和資源密集型產(chǎn)業(yè)，實(shí)現(xiàn)勞動(dòng)力資源的合理配置與回流，促進(jìn)西部地區(qū)經(jīng)濟(jì)社會(huì)發(fā)展。

注釋

①根據(jù)國(guó)家統(tǒng)計(jì)局2018年中國(guó)總?cè)丝跀?shù)和流動(dòng)人口數(shù)計(jì)算得出。

②高魯鵬、梁文舉、趙軍：《氣候變化對(duì)黑土有機(jī)碳庫(kù)影響模擬研究》，《遼寧工程技術(shù)大學(xué)學(xué)報(bào)(自然科學(xué)版)》2005年第2期。

③韓貴鋒、徐建華：《人口與經(jīng)濟(jì)發(fā)展對(duì)植被的影響研究——以重慶市為例(英文)》，《長(zhǎng)江流域資源與環(huán)境》2008年第5期。

④地區(qū)氣溫升高，原適合該地區(qū)生長(zhǎng)的農(nóng)作物因土壤中有機(jī)碳含量減少，從而導(dǎo)致減產(chǎn)。

⑤C. Gray and V. Mueller， “Natural Disasters and Population Mobility in Bangladesh，”ProceedingsoftheNationalAcademyofSciences， vol.109， no.16， 2012， pp. 6000-6005．

⑥參見(jiàn)E. EI-Hinnawi，EnvironmentalRefugees， Nairobi: United Nations Environmental Programme， 1985．

⑦A. Swain， “Environmental Migration and Conflict Dynamics: Focus on Developing Regions，”ThirdWorldQuarterly， vol.17， no.5，1996， pp.959-973．

⑧余慶年、施國(guó)慶：《環(huán)境、氣候變化和人口遷移》，《中國(guó)人口·資源與環(huán)境》2010年第7期。

⑨稅偉、徐國(guó)偉、蘭肖雄、王雅文、馬菁：《生態(tài)移民國(guó)外研究進(jìn)展》，《世界地理研究》 2012年第1期。