摘 要：地球的吸熱散熱及海平面溫度是氣候變化重要因素，為建立刻畫氣候變化的模型，首先采用地球表面平均溫度、海洋表面溫度、降水量以及CO2濃度作為模型的內(nèi)生變量。考慮影響氣候的因素之間也相互影響，選擇不事先約束條件的VAR模型來分析氣候變化，同時(shí)將月份作為外生變量用以去除時(shí)間上的趨勢(shì)。模型經(jīng)最優(yōu)滯后階數(shù)檢驗(yàn)后選定滯后5階，且模型通過了穩(wěn)定經(jīng)檢驗(yàn)。然后經(jīng)Granger因果檢驗(yàn)和方差分解后，得出各變量普遍能Granger引起其他變量和各結(jié)構(gòu)沖擊對(duì)內(nèi)生變量變化的貢獻(xiàn)度。最后在此模型基礎(chǔ)上，利用選取數(shù)據(jù)預(yù)測(cè)未來25年的氣候變化，經(jīng)預(yù)測(cè)在2045年全球表面溫度約將上漲0.24℃，全球海平面溫度約將上漲2.09℃，降水量約將上漲1.90mm，CO2濃度約將上漲4.27ppm。

關(guān)鍵詞：氣候變化預(yù)測(cè);HP濾波;向量自回歸模型

中圖分類號(hào)：TP319 ? ? ? ?文獻(xiàn)標(biāo)志碼：A ? ? ?文章編號(hào)：1673-291X（2021）15-0146-04

一、研究背景

對(duì)于全球變暖的問題，科學(xué)界已做出很多研究，但對(duì)于全球變暖是否主要由人類活動(dòng)造成，科學(xué)界尚有分歧，在非科學(xué)界更是爭(zhēng)論四起[1]。為使非專業(yè)人士更好地了解全球變暖現(xiàn)象及氣候變化的因素和趨勢(shì)，本文將專業(yè)復(fù)雜的氣候模型簡(jiǎn)化處理，建立更為直觀易懂的簡(jiǎn)化模型分析氣候變化的因素、趨勢(shì)，得出更易于理解直觀的結(jié)論。

首先預(yù)測(cè)模型需要考慮到地球的吸熱、散熱。地球吸收的熱量來源于太陽(yáng)的可見光輻射，即短波輻射，而地球的散熱屬于長(zhǎng)波輻射。當(dāng)太陽(yáng)光透過大氣層時(shí)，短波輻射透過率很高，但當(dāng)?shù)厍蛳蛲馍釙r(shí)，大氣層中的水汽、CO2等卻吸收了大部分的長(zhǎng)波輻射，由此導(dǎo)致了地球熱量的不斷積累?？梢?，影響地球吸熱散熱的主要因素正是大氣中這些能吸收長(zhǎng)波輻射的溫室氣體[2～5]。

另外，海洋也能吸收大量的熱量，水的比熱容比陸地的比熱容更大，其吸收的熱量對(duì)地球的表面溫度也具有一定影響[6～7]。所以我們將海洋表面溫度這一表征作為影響地球吸熱、散熱的影響因素或者被影響因素來看待，具體將在模型中進(jìn)一步分析。

對(duì)于建立刻畫氣候的模型，考慮到本文主要目的為幫助增強(qiáng)非專業(yè)人士對(duì)氣候的認(rèn)識(shí)，此題應(yīng)將復(fù)雜問題簡(jiǎn)單化，選取少量但關(guān)鍵的因素作為模型變量。同時(shí)考慮到引起氣候變化的滯后因素以及各個(gè)變量之間的相互影響關(guān)系，本文將采用向量自回歸模型（VAR模型），將每一個(gè)內(nèi)生變量作為整個(gè)模型每一個(gè)內(nèi)生變量的滯后因素，來探究考慮滯后及相互影響關(guān)系的變量對(duì)于氣候的影響[8～12]。

二、模型

（一）數(shù)據(jù)選取及處理

本文選取了1980年1月至2017年6月的地球表面平均溫度、海洋表面溫度、降水量以及CO2濃度數(shù)據(jù)，數(shù)據(jù)來源于NOAA ICOADS數(shù)據(jù)集。

為了簡(jiǎn)化計(jì)算、縮小量值，此題使用SPSS軟件對(duì)數(shù)據(jù)進(jìn)行了歸一化處理，數(shù)據(jù)樣本見表1，變量簡(jiǎn)寫及含義見下頁(yè)表2。

（二）建立模型

經(jīng)過初步分析，我們認(rèn)為地球表面溫度、海洋表面溫度、降水量以及CO2濃度這4個(gè)變量為相互聯(lián)系的時(shí)間序列變量，即它們之間互相影響，很難從實(shí)際角度直觀判斷各自的因果關(guān)系。據(jù)此，我們選擇建立不帶有事先約束條件的向量自回歸模型（下文簡(jiǎn)稱“VAR模型”）。該模型的一大特點(diǎn)是通過把系統(tǒng)中每一個(gè)內(nèi)生變量作為系統(tǒng)中所有內(nèi)生變量的滯后值的函數(shù)來構(gòu)造模型，從而回避了結(jié)構(gòu)化模型的要求，其常用于預(yù)測(cè)相互聯(lián)系的時(shí)間序列系統(tǒng)以及分析隨機(jī)擾動(dòng)對(duì)變量系統(tǒng)的動(dòng)態(tài)影響，能夠滿足此題的分析要求[13～16]。模型表達(dá)式為：

yt=A1yt-1+…+Apyt-p+BXt+εt

其中，yt是k維內(nèi)生變量向量，Xt是d維外生變量向量，p是滯后階數(shù)，樣本個(gè)數(shù)為T，k×k維矩陣A1，…Ap和k×d維矩陣B是要被估計(jì)的系數(shù)矩陣。εt是k維擾動(dòng)向量。

第一步，建立初步向量自回歸模型，將所有序列，即surf sst prec CO2時(shí)間序列作為內(nèi)生變量導(dǎo)入到Eviews軟件中。此外，將變化量月份作為模型的外生變量，建立初步的向量自回歸模型。

第二步，進(jìn)行滯后階數(shù)檢驗(yàn)，得出最優(yōu)模型的滯后階數(shù)為5階，按照最優(yōu)的滯后階數(shù)重新建立VAR模型，模型結(jié)果見表3。

首先，輸出的第一部分和第二部分是每個(gè)方程的標(biāo)準(zhǔn)OLS 回歸統(tǒng)計(jì)量，即以每一列變量各自作為因變量，計(jì)算方程的結(jié)果。從輸出結(jié)果可以看出，各因變量對(duì)滯后5階的變量普遍顯著，且第一個(gè)方程和第四個(gè)方程的調(diào)整擬合優(yōu)度都較高，說明擬合效果較好。雖然第二個(gè)和第三個(gè)方程的調(diào)整擬合優(yōu)度偏低，但因本模型的特點(diǎn)，對(duì)方程的擬合優(yōu)度沒有嚴(yán)格要求，為非重要因素，依然可以利用這個(gè)模型進(jìn)行下一步的預(yù)測(cè)分析。

其次，輸出的第三部分為整個(gè)VAR 模型的回歸統(tǒng)計(jì)量，整體AIC，SC都較小，Log likehood也較小，擬合效果均較好，且模型滿足穩(wěn)定性條件。

第三步，格蘭杰因果檢驗(yàn)，格蘭杰因果檢驗(yàn)?zāi)軌驒z驗(yàn)一個(gè)變量的滯后變量是否可以引入到其他變量方程中，如果其受到其他變量的滯后影響，則稱其具有格蘭杰因果關(guān)系。本文采取兩種方法進(jìn)行格蘭杰因果檢驗(yàn)，其結(jié)果見表4。

從表的數(shù)據(jù)不同可以看出，Eviews中不同的檢驗(yàn)處理方法，得出數(shù)據(jù)有略微差異，從表4數(shù)據(jù)中得出SST與CO2均不能Granger引起Prec，Prec與CO2也均不能Granger引起SST，而結(jié)果除Prec接受不能Granger引起Surf的假設(shè)，其余均顯著拒絕原假設(shè)，即表現(xiàn)為具有相互Granger關(guān)系。綜合兩種檢驗(yàn)結(jié)果，我們大體可以得出，在滯后5 階，各變量的滯后變量均對(duì)其他變量有影響。

第四步，方差分解，這一步檢驗(yàn)可得出每一個(gè)結(jié)構(gòu)沖擊對(duì)內(nèi)生變量變化（通常用方差來衡量）的貢獻(xiàn)度，將穩(wěn)定時(shí)期的貢獻(xiàn)度羅列在表 5中。

由表5可得出，在排除非外界因素影響的情況下，對(duì)Surf變量來說，對(duì)其影響效果最大的變量依次為Prec，CO2和SST;對(duì)SST影響最大的變量依次為Surf，CO2，Prec;對(duì)Prec影響最大的變量依次為Surf，SST，CO2;而對(duì)CO2影響最大的變量依次為Surf，Prec，自身及SST，其中有兩個(gè)變量的影響效果大于CO2對(duì)自身的影響效果。

（三）預(yù)測(cè)數(shù)據(jù)

根據(jù)上述模型，在無(wú)外界變量干擾下，預(yù)測(cè)未來25年變動(dòng)數(shù)據(jù)樣例見表 6，各變量的預(yù)測(cè)見圖 1。

由表6和圖1可以預(yù)測(cè)未來25年Surf上漲0.24℃左右，SST上漲2.09℃左右，Prec上漲1.90mm左右，CO2上漲4.27ppm左右。

結(jié)語(yǔ)

本文采用VAR模型將地球表面平均溫度（Surf）、海洋表面溫度（SST）、降水量（Prec）以及CO2濃度（CO2）作為模型內(nèi)生變量的滯后因素，探究其各個(gè)因素對(duì)氣候變化模型的影響。同時(shí)將月份作為外生變量用以去除時(shí)間上的趨勢(shì)，并依據(jù)所建立的模型預(yù)測(cè)未來25年的氣候變化趨勢(shì)，經(jīng)預(yù)測(cè)在2045年全球表面溫度約將上漲0.24℃，全球海平面溫度約將上漲2.09℃，降水量約將上漲1.90mm，CO2濃度約將上漲4.27ppm。

參考文獻(xiàn)：

[1] ?P.H.Gleick，R.M.Adams，R.M.Amasino.Climate Change and the Integrity of Science[J].Science，2010，（328）：689-690.

[2] ?李正.末次盛冰期以來巴倫支?！９藕Ｑ蟓h(huán)境及海冰研究進(jìn)展[C]//浙江省地理學(xué)會(huì).浙江省地理學(xué)會(huì)2020年學(xué)術(shù)年會(huì)暨環(huán)境與健康高峰論壇論文摘要集[C].浙江省地理學(xué)會(huì)，2020：1.

[3] ?Davidson Samuel G，Hesp Patrick A，Silva Graziela Miot da.Controls on dune scarping[J].Progress in Physical Geography，2020，（6）.

[4] ?Gulliver John S..Climate Change Adaptation：Civil and Environmental Engineering Challenge[J].Journal of Hydraulic Engineering，2020，（12）.

[5] ?武勝利，劉強(qiáng)吉.近50a新疆巴州地區(qū)氣溫與降水時(shí)空變化特征[J].干旱氣象，2016，（4）：610-616.

[6] ?董奎，劉兵，鄭旭榮.金溝河流域氣溫、降水與年徑流量演變規(guī)律研究[J].人民黃河，2010，（5）：32-33.

[7] ?李鄒，方雙喜，和春榮，等.香格里拉本底站大氣CO2濃度及變化特征初步研究[J].環(huán)境化學(xué)，2012，（12）.

[8] ?李藝珍，毛建剛，張明，岳春芳.基于VAR模型的新疆金溝河氣象要素與徑流關(guān)系分析[J].水資源與水工程學(xué)報(bào)，2020，（5）：80-86.

[9] ?趙娟，穆興民，高鵬，等.基于VAR模型分析徑流對(duì)降水和水土保持措施的響應(yīng)[J].泥沙研究，2017，（4）：37-44.

[10] ?劉盈曦，彭貴芬，陳先剛，楊宇明.香格里拉未來50a主要?dú)夂颦h(huán)境要素變化預(yù)估——基于小波分析和多元VAR回歸預(yù)估模型[J].資源科學(xué)，2016，（9）.

[11] ?王國(guó)慶，張建云，李巖，等.海河流域基于多模式的未來氣候變化趨勢(shì)[J].資源科學(xué)，2014，（5）.

[12] ?周兆媛，張時(shí)煌，高慶先，等.京津冀地區(qū)氣象要素對(duì)空氣質(zhì)量的影響及未來變化趨勢(shì)分析[J].資源科學(xué)，2014，（1）：191-199.

[13] ?趙娟.基于VAR模型的典型流域水沙變化及其對(duì)降水與水土保持措施的動(dòng)態(tài)響應(yīng)[D].咸陽(yáng)：西北農(nóng)林科技大學(xué)，2019.

[14] ?李雙，李哲，杜建括，王淑新，邢海虹.陜西省水資源利用與經(jīng)濟(jì)增長(zhǎng)動(dòng)態(tài)關(guān)系的VAR模型分析[J].生態(tài)經(jīng)濟(jì)，2020，（10）：146-154.

[15] ?Sims C A.Macroeconomics and reality[J].Econometrica，1980，（1）：1-48.

[16] ?朱慧明，劉智偉.時(shí)間序列向量自回歸模型的貝葉斯推斷理論[J].統(tǒng)計(jì)與決策，2004，（1）：11-12.

Simplified Climate Prediction Based on VAR Model

JIN Zi-qi

（School of Management，University of Shanghai for Science and Technology，Shanghai 200093，China）

Abstract：The earths heat absorption and heat dissipation and sea-level temperature are important factors of climate change.In order to establish a model to describe climate change，the average surface temperature of the earth，sea surface temperature，precipitation and carbon dioxide concentration are used as the endogenous variables of the model.Considering that the factors affecting climate also affect each other，the VAR model without prior constraints is selected to analyze climate change，and the month is used as exogenous variable to remove the trend of time.After the model is tested by the optimal lag order，5-order lag is selected，and the model has passed the stability test.Then，through the Granger causality test and variance decomposition，we can get the contribution of other variables caused by Granger and the structural impact to the changes of endogenous variables.Finally，based on this model，the climate change in the next 25 years is predicted by using selected data.It is predicted that the global surface temperature will rise by 0.24℃，the global sea level temperature will rise by 2.09℃，the precipitation will increase by 1.90 mm，and the carbon dioxide concentration will increase by 4.27 ppm.

Key words：climate change prediction;HP filter;vector autoregressive model

收稿日期：2020-09-07

作者簡(jiǎn)介：金子琪（1993-），女（朝鮮族），黑龍江佳木斯人，碩士研究生，從事數(shù)量經(jīng)濟(jì)研究。