基于超越對數(shù)生產(chǎn)函數(shù)的氣象條件變化對京津冀地區(qū)行業(yè)產(chǎn)出的影響研究

王子忠 潘志華 宋玉 黃娜 龍步菊 倫飛

（1 中國農(nóng)業(yè)大學資源與環(huán)境學院，北京 100193；2 山西省氣象臺，太原 030006）

0 引言

氣象條件作為經(jīng)濟發(fā)展的重要環(huán)境要素，幾乎時刻都在制約或者促進經(jīng)濟的發(fā)展[1]。伴隨著現(xiàn)代社會生產(chǎn)力水平的突飛猛進和產(chǎn)業(yè)結構的日益完善，人類逐漸降低氣象條件變化帶來的影響，但是仍然無法避免氣象災害與氣候波動帶來的短期和長期影響。美國學者Dutton[2]曾指出，美國1/3的國內生產(chǎn)總值受到氣象條件的影響，約3萬億美元的年收入有一定程度的天氣和氣候風險。雖然依靠的是主觀定性分析，但是這個觀點被廣泛引用。迄今，國內外基于計量經(jīng)濟學方法圍繞氣象條件與經(jīng)濟發(fā)展之間開展了很多研究。Richard[3]利用多元回歸分析方法，分析了天氣變化對荷蘭國際旅游的影響，發(fā)現(xiàn)上一年適宜的夏季氣溫有利于下一年度假人數(shù)的增長。Dell等[4]分析了近半個世紀各個國家溫度沖擊和經(jīng)濟增長之間的關系，氣溫上升對窮國不利，對富裕國家影響不明顯。Naeem[5]在生產(chǎn)函數(shù)的基礎上進行拓展，以GDP增長率為因變量，以城市化率、人口增長率等為控制變量，添加氣溫和降水因子，建立了亞洲國家面板數(shù)據(jù)模型，結果表明氣溫、降水與GDP增長率之間存在負相關關系。Roy等[6]利用建立的時間序列回歸模型表明孟加拉國年平均氣溫每升高1℃，GDP增長率平均下降0.44%。Larsen[7]基于生產(chǎn)函數(shù)，在勞動力、資本以及能源生產(chǎn)要素基礎上，引入氣溫和降水要素研究了美國各州11個經(jīng)濟部門對氣候的敏感性，結果表明不同地區(qū)、不同行業(yè)的氣候敏感性存在差異，其中對西部地區(qū)影響最大。也有研究認為氣候與生產(chǎn)之間是非線性關系，Marshall等[8]認為幾乎所有國家的整體生產(chǎn)率與氣溫之間是非線性的。Zhang等[9]基于工業(yè)企業(yè)數(shù)據(jù)與氣溫建立的面板模型，證實氣溫與工業(yè)生產(chǎn)率之間存在倒“U”型關系。中國在這方面的研究雖然起步較晚，但仍有很多有意義的研究，并且對多種研究方法進行了探索。鞏在武[10]利用灰色關聯(lián)法，以蘇皖豫為研究區(qū)域,分析了氣溫、降水與三大產(chǎn)業(yè)存在相關性。吳息等[11]利用簡單的回歸統(tǒng)計方法，建立了氣象條件與工業(yè)、能源的相關關系，表明夏季升溫對工業(yè)生產(chǎn)不利。羅慧等[12]利用相關性分析方法分析了陜西省氣象要素與經(jīng)濟行業(yè)之間的相關關系，結果表明國民經(jīng)濟與氣象條件存在明顯相關性。丑潔明等[13]認為應該將生產(chǎn)活動與氣候變化研究結合起來，在Cobb-Douglass生產(chǎn)函數(shù)（簡稱C-D生產(chǎn)函數(shù)）中加入了氣象要素，通過評價氣候變化對糧食產(chǎn)量的影響，驗證了其假設合理性。羅慧等[14]基于C-D生產(chǎn)函數(shù)研究表明中國省域行業(yè)經(jīng)濟對氣象條件的敏感性總體上北方大于南方，西部大于東部，且農(nóng)業(yè)與建筑行業(yè)對氣象條件變化最為敏感。

綜合來看，衡量氣候和氣象條件變化對經(jīng)濟的影響，目前主要有兩種方法：一種是利用綜合評估模型，是現(xiàn)階段用于模擬不同氣候變化情景下造成的可能經(jīng)濟影響和計算減緩氣候變化成本的主流工具[15]，在IPCC第六次評估報告中得到了廣泛應用，最新的預估表明，減少排放不僅能降低平均預期損害，還能降低災難性事件的風險[16]，但綜合評估模型始終具有很大的不確定性[17]；另一種方法是本文主要討論的，即基于計量經(jīng)濟學的實證研究方法，比如C-D生產(chǎn)函數(shù)[18-23]、協(xié)整理論[24]、灰色關聯(lián)分析[10]以及相關分析方法[12]等都是度量氣象條件和氣候變化對社會經(jīng)濟影響的有益嘗試，從實證角度證明氣候波動與經(jīng)濟生產(chǎn)之間的關系，彌補了氣候變化綜合評估模型經(jīng)驗證據(jù)不足的缺陷[25]。其中，基于C-D生產(chǎn)函數(shù)的研究應用最為廣泛，在生產(chǎn)函數(shù)中加入氣象要素說明研究者已經(jīng)開始逐步重視氣象條件與經(jīng)濟產(chǎn)出的關系，但是卻忽略了氣象條件對生產(chǎn)要素的影響，實際上氣象條件不僅影響最后的產(chǎn)出，而且對于投入要素也會產(chǎn)生影響，并且生產(chǎn)要素與產(chǎn)出之間具有非線性關系。為此，本文拓展應用生產(chǎn)函數(shù)，以京津冀地區(qū)為例，建立超越對數(shù)生產(chǎn)函數(shù)定量研究不同行業(yè)對氣象條件變化的敏感性，以期為制定氣候變化適應對策提供科學依據(jù)。

1 數(shù)據(jù)與方法

1.1 生產(chǎn)函數(shù)

C-D生產(chǎn)函數(shù)的一般形式為：

式中，Q為行業(yè)產(chǎn)值；L為勞動投入量；K為資本投入量；A表示技術水平，且A＞0；α為勞動力產(chǎn)出的彈性系數(shù)，β為資本產(chǎn)出的彈性系數(shù)。

C-D生產(chǎn)函數(shù)具有經(jīng)濟學家所關心的許多性質，同時具有參數(shù)求取方便等特點。隨著經(jīng)濟理論的不斷完善，人們不斷嘗試將更多的生產(chǎn)要素納入到生產(chǎn)函數(shù)中[23]。其推廣形式可以將需要討論的諸多因子引入模型進行分析，這為改進模型提供了理論基礎[13]。將上式改寫為傳統(tǒng)的生產(chǎn)函數(shù)形式：

式中，βK、βL、βW分別代表各個投入要素的產(chǎn)出彈性；，借鑒Fraser的處理方法，在本文中將技術變化的影響表述為一個時間趨勢函數(shù)[7，26]；氣象要素W的選取，可以有多種指標，如溫度、降水或者極端氣象條件等，這些指標均可當成投入要素加入函數(shù)進行分析[13]。

1.2 超越對數(shù)生產(chǎn)函數(shù)

將單個或者多個氣象要素作為生產(chǎn)要素引入模型，這種方法在一定程度上固然可以揭示氣象條件變化對某部門的影響，但是卻忽略了各個因子之間的相互作用?？紤]到各個因子之間的相互作用，同時也考慮到經(jīng)濟生產(chǎn)與各個影響因子之間的非線性關系，超越對數(shù)生產(chǎn)函數(shù)引入了生產(chǎn)要素的交叉項和平方項，可以解釋要素之間的相互影響[23,27]。本文嘗試將氣象要素引入超越對數(shù)生產(chǎn)函數(shù)，來研究氣象要素對部門經(jīng)濟產(chǎn)出的影響。引入氣象要素后，其變化不僅影響最終產(chǎn)出，且對其他投入要素也會同時產(chǎn)生影響，比如某一年發(fā)生了較大的氣候波動，不僅會影響最終產(chǎn)出，而且會對投資和勞動力投入產(chǎn)生影響。

超越對數(shù)生產(chǎn)函數(shù)允許輸入不同變量之間的替代彈性，是一種易于估計且包容性很強的變彈性系數(shù)對數(shù)生產(chǎn)函數(shù)，可以很好地反映投入項的相互影響和平方項的非線性影響[23]。所以本文采用超越對數(shù)生產(chǎn)函數(shù)的重要假設前提是，氣象要素不僅影響經(jīng)濟部門的最終產(chǎn)值，即被解釋變量，還影響模型中的其他生產(chǎn)要素。本文中假設氣象條件作為生產(chǎn)要素，對勞動力投入和資本投入都會產(chǎn)生一定的影響。

本文擬在資本存量和勞動力的基礎上，引入采暖度日、降溫度日及降水距平三個氣象要素，采用與面板模型相結合的擬合方式。以建筑業(yè)為例，建立如下模型：

式中，H為采暖度日（℃·d）；C為降溫度日（℃·d）；P為降水距平百分率（%）；K為資本存量（億元）；L為勞動力數(shù)量（萬人）；Q為生產(chǎn)總值（億元）；δ、βk等為系數(shù)。在此模型中引入了時間趨勢變量t，利用其捕捉其他隨時間變化的變量，如技術變化趨勢、勞動者素質的提高等時間趨勢變量?？紤]到采暖度日和降溫度日沒有實質性的交叉影響，所以未引入其交叉項。

各要素的產(chǎn)出彈性不再簡單的等同于其系數(shù)，包含了交叉項和平方項的信息[7,23]?；谏鲜?，可以得出H對建筑業(yè)產(chǎn)出的邊際產(chǎn)出彈性為：

1.3 面板模型擬合方法的選擇

面板模型在氣候變化計量經(jīng)濟學領域得到了廣泛的應用[28]。在本研究中，利用京津冀地區(qū)三個?。ㄊ校└餍袠I(yè)經(jīng)濟數(shù)據(jù)和氣象數(shù)據(jù)組成面板數(shù)據(jù)模型，分析產(chǎn)出對氣象因子的敏感性。采用雙向固定效應模型，具有時間和個體兩個維度的信息，即對于不同的時點、不同的個體都有不同截距的模型，可以在一定程度上減少隨機誤差項帶來的異方差、自相關以及變量之間的共線性對估計的影響[29]。其次，本文中所構建的面板模型，時間序列大于截面數(shù)量，為長面板。

式中，Yit為被解釋變量；αi為個體（橫截面）影響；αt為時點（時序）影響；εit為殘差；i為個體，i＝1, 2, 3,…,N；t為時間，t＝1, 2, 3,…,T。

表1中給出了面板數(shù)據(jù)模型樣本，個體效應使用虛擬變量來控制，對于時間效應，為了避免損失較多的自由度，使用時間趨勢量控制[29]。

表1 面板數(shù)據(jù)模型樣本數(shù)據(jù)示例Table 1 Sample data example of panel data model

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

1.4.1 社會經(jīng)濟數(shù)據(jù)

社會經(jīng)濟統(tǒng)計資料主要來源于《中國統(tǒng)計年鑒》《北京統(tǒng)計年鑒》《天津統(tǒng)計年鑒》《河北經(jīng)濟年鑒》。經(jīng)濟部門選擇具有代表性的農(nóng)業(yè)、工業(yè)、建筑業(yè)、交通運輸倉儲郵政業(yè)、批發(fā)零售業(yè)。

本文選用年數(shù)據(jù)進行分析，時間為1985—2017年。具體計算時，使用GDP平減指數(shù)去除由于物價上漲或者通貨膨脹而增長的那一部分產(chǎn)值。固定資產(chǎn)折舊率為5%[30]；勞動力數(shù)據(jù)采用年末勞動力人數(shù)。

1.4.2 氣象數(shù)據(jù)

研究所用的氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)共享網(wǎng)（http://www.nmic.cn/）地面氣候數(shù)據(jù)資料集，包括京津冀地區(qū)27個氣象觀測點1961—2017年逐日氣象數(shù)據(jù)。研究選取采暖度日、降溫度日及年降水距平百分率作為氣象指標。度日是指日平均溫度與規(guī)定閾值溫度的離差，可以在一定程度上反映人生活最舒適的環(huán)境溫度[31]。結合中國《節(jié)能設計標準》和在多數(shù)研究中使用的閾值[32-33]，選取18 ℃和26 ℃作為基礎溫度，計算方法如下：

式中，Ti為日平均氣溫小于18 ℃的當天氣溫（℃）；Zi為日平均氣溫高于26 ℃的當天氣溫（℃）。

降水指標選取降水距平百分率，表示在一定時段內降水偏離多年平均狀態(tài)的離差程度，在一定程度可以反映降水量的年際差異。由于需要將數(shù)據(jù)對數(shù)化，降水距平難免有負值，所以本文對所有值統(tǒng)一加50，這并不會改變數(shù)據(jù)的波動特征，仍然是有效的。

式中，P為降水距平百分率（%），Pi為第i年年降水量（mm），ˉP為1961—1990年平均降水量（mm）。

2 結果與分析

2.1 京津冀地區(qū)1985—2017 年降溫度日和采暖度日變化特征

圖1為京津冀地區(qū)1985—2017年降溫度日和采暖度日變化特征?？梢钥闯?，1985年以來，京津冀地區(qū)采暖度日不斷減少，每10年減少68.4 ℃·d；降溫度日不斷增加，每10年增加13.2 ℃·d，且降溫度日變化的顯著性大于采暖度日。北京市采暖度日多年平均值為3175.43 ℃·d，氣候傾向率為－85.8 ℃·d/10 a；降溫度日波動上升趨勢明顯，氣候傾向率為15.4 ℃·d/10 a。天津市多年平均采暖度日為2819.68 ℃·d，呈波動下降趨勢；降溫度日上升趨勢更明顯，氣候傾向率為16.32 ℃·d/10 a，多年平均值為80.25 ℃·d。河北省采暖度日也呈明顯的下降趨勢，多年平均值為3407.53 ℃·d，降溫度日多年平均值為49.36 ℃·d，氣候傾向率為7.9 ℃·d/10 a。三個?。ㄊ校┎膳热铡⒔禍囟热毡憩F(xiàn)出相近的年際變化特征，京津冀地區(qū)升溫趨勢非常顯著，其中，北京市度日趨勢變化最為顯著，且三個?。ㄊ校┑慕禍囟热兆兓@著性大于采暖度日。

圖1 1985—2017年京津冀地區(qū)降溫度日和采暖度日變化特征（a）京津冀地區(qū)；（b）北京市；（c）河北?。唬╠）天津市Fig.1 Variation characteristics of cooling degree days and heating degree days in Beijing-Tianjin-Hebei region from 1985 to 2017(a) Beijing-Tianjin-Hebei region; (b) Beijing; (c) Hebei; (d) Tianjin

2.2 京津冀地區(qū)1985—2017年經(jīng)濟變化特征

圖2為京津冀地區(qū)分?。ㄊ校┲饕?jīng)濟變量隨時間的變化特征。GDP為實際GDP，即去除了通貨膨脹的影響。圖2a1—2a3為北京市主要經(jīng)濟變量隨時間變化趨勢圖?？梢钥闯觯本┦泄I(yè)和批發(fā)零售業(yè)產(chǎn)值增長速度最快，特別是工業(yè)產(chǎn)值。批發(fā)零售業(yè)增長迅速與北京市人口的迅速增長關系密切。建筑業(yè)和交通運輸倉儲郵政業(yè)產(chǎn)值增長水平近乎持平。農(nóng)業(yè)在1985—2017年，產(chǎn)值保持較為穩(wěn)定，增幅較小。工業(yè)和交通運輸倉儲郵政業(yè)資本存量最為龐大，其余三個部門資本存量逐年增長。工業(yè)部門勞動力數(shù)量呈下降趨勢，批發(fā)零售業(yè)上升趨勢明顯，其余部門基本保持在一定范圍，每年略有波動，與每年各部門對勞動力的需求有關。圖2b1—2b3為天津市五個行業(yè)主要經(jīng)濟變量隨時間變化趨勢圖?？梢钥闯觯鞑块T生產(chǎn)總值和資本存量在2001年后都迅速增長，部門生產(chǎn)總值在有的年份有小幅下降，但是基本上都呈現(xiàn)為逐年上升趨勢。其次批發(fā)零售業(yè)的上升趨勢也很顯著。工業(yè)行業(yè)的勞動力數(shù)量最多，在2001年，勞動力大幅度下降，之后穩(wěn)步回升，但在2013年后又呈下降趨勢。除農(nóng)業(yè)行業(yè)勞動力數(shù)量呈下降趨勢外，其他行業(yè)都呈穩(wěn)步上升趨勢。圖2c1—2c3為河北省五個行業(yè)經(jīng)濟變量的時間變化趨勢圖。可以看出，工業(yè)產(chǎn)值占的比重最大，1991年為214億，到2017年增長到了2300多億，增長了10倍左右。對應的資本存量也迅速積累，1991年資本存量為313億，到2017年達到接近17000億的資本存量，可見工業(yè)行業(yè)在1991年之后幾乎是爆發(fā)式增長，對應著工業(yè)勞動力數(shù)量也達到了1000萬。其他行業(yè)產(chǎn)值、資本存量以及勞動力數(shù)量也是穩(wěn)步上升。

圖2 京津冀地區(qū)1985—2017年社會經(jīng)濟變化特征Fig.2 Characteristics of socioeconomic changes in Beijing-Tianjin-Hebei region from 1985 to 2017

2.3 基于超越對數(shù)生產(chǎn)函數(shù)的行業(yè)敏感性分析

表2為各行業(yè)的超越對數(shù)生產(chǎn)函數(shù)系數(shù)。采用最小二乘虛擬變量法（LSDV）做回歸估計，最大程度減少異方差、自回歸及截面相關的影響。以面板數(shù)據(jù)為基礎，利用雙向固定效應模型（公式（5））得出各部門超越對數(shù)生產(chǎn)函數(shù)系數(shù)，再將系數(shù)以及各?。ㄊ校┑臍庀笠睾徒?jīng)濟變量序列數(shù)據(jù)代入公式（4），即可得出采暖度日、降溫度日以及降水距平對分行業(yè)產(chǎn)出的敏感性數(shù)值，通過取其多年平均值得到京津冀地區(qū)氣象要素的邊際產(chǎn)出彈性系數(shù)。

表2 各部門的超越對數(shù)生產(chǎn)函數(shù)系數(shù)Table 2 Coefficients of translog production function for each sector

圖3為氣象要素對部門產(chǎn)出的邊際產(chǎn)出彈性，可以看出，采暖度日對建筑業(yè)、交通運輸倉儲郵政業(yè)、工業(yè)、批發(fā)零售業(yè)以及農(nóng)業(yè)的邊際產(chǎn)出彈性分別為－0.1501%、0.6197%、0.9309%、0.4473%以及－0.7311%，除了建筑業(yè)和農(nóng)業(yè)之外，都表現(xiàn)為正向影響；降溫度日對五個經(jīng)濟行業(yè)的邊際產(chǎn)出彈性分別為0.0274%、－0.0208%、－0.0709%、－0.1028%以及0.0282%，對建筑業(yè)和農(nóng)業(yè)具有微弱的正向作用；降水距平對五個部門的邊際影響分別為－0.1643%、－0.13 65%、－0.0005%、－0.0817%以及0.0407%，除農(nóng)業(yè)之外，都表現(xiàn)為負向影響。

圖3 氣象要素對五個部門的邊際產(chǎn)出彈性的影響Fig.3 The impact of meteorological factors on the marginal output elasticity of five sectors

上述結果符合實際情況。采暖度日對建筑業(yè)表現(xiàn)為負向影響，即采暖度日每增加1%，則建筑業(yè)產(chǎn)出下降0.15%，建筑業(yè)主要包括建筑工程和建筑安裝工程，以戶外勞動為主。采暖度日的增加，意味著冷日數(shù)的增加，不利于戶外勞動作業(yè)，且混凝土易受凍破壞結構。北方冬季建筑業(yè)有冬歇期，變暖可延長施工期，但是目前采暖度日是減少的，所以采暖度日的減少對建筑業(yè)是有利的。采暖度日對交通運輸倉儲郵政業(yè)、工業(yè)以及批發(fā)零售業(yè)為正向影響。目前，工業(yè)生產(chǎn)大多為室內生產(chǎn)，氣溫降低對其影響并不是很明顯。從經(jīng)驗解釋角度看，冷日數(shù)的增加會增加對交通運輸?shù)男枨螅貏e是在北方地區(qū)，供暖季的燃料運輸是剛性需求。氣象要素對批發(fā)零售業(yè)的影響，一般是通過天氣變化影響消費者的行為，進而影響產(chǎn)值。

對比采暖度日邊際產(chǎn)出彈性和降溫度日邊際產(chǎn)出彈性，發(fā)現(xiàn)了截然相反的情況，降溫度日對建筑業(yè)為正向影響，雖然降溫度日的增加可能意味著高溫日數(shù)的上升，對戶外生產(chǎn)不利，但是目前采暖度日的大幅減少有利于戶外勞動時間的延長，對建筑業(yè)生產(chǎn)有利，所以很有可能采暖度日的有利影響抵消了降溫度日帶來的不利影響。降溫度日增加對工業(yè)生產(chǎn)不利，Zhang等[9]的研究表明，氣溫與工業(yè)生產(chǎn)率之間存在倒“U”型關系，并且表現(xiàn)出對高溫的敏感性；高溫天氣環(huán)境下會降低工人的生產(chǎn)積極性，也會導致工業(yè)機械故障率增加，從而降低生產(chǎn)效率，并且在高溫天氣下降溫需求增加，也會增加工業(yè)生產(chǎn)成本。降溫度日增加對交通運輸倉儲業(yè)和批發(fā)零售業(yè)不利，一方面是高溫天氣對交通設施和司機駕駛狀態(tài)都會造成負面影響，同時也影響消費者的消費欲望。在炎熱的天氣下，人們外出購物的想法明顯會受到影響。隨著全球氣候變暖，高溫熱浪天氣不斷增加，并且不斷刷新著記錄，所以務必要預防高溫天氣帶來的負面影響。

降水對多數(shù)部門都為負向影響（除了農(nóng)業(yè)），對建筑業(yè)影響最為明顯，這與孫鑒鋒等[21]的研究表明降水對大多數(shù)部門為負向產(chǎn)出效應的結果一致。通常認為降水對戶外勞動影響是顯著的，降水天氣下一般無法開展戶外勞動，以戶外勞動為主的建筑業(yè)一般都會停工。降水對交通運輸業(yè)和批發(fā)零售業(yè)也有著較大的負向影響，降水對交通不利，容易引發(fā)交通堵塞和交通事故，強降水甚至可能會導致交通癱瘓。北京市2012年7·21特大暴雨，導致全市道路以及橋梁多處損毀，多條鐵路線臨時停運，直接經(jīng)濟損失達116.4億元。與交通運輸倉儲業(yè)聯(lián)系緊密的批發(fā)零售業(yè)同樣也會受到影響。

總體來看，交通運輸倉儲郵政業(yè)、工業(yè)以及批發(fā)零售業(yè)對采暖度日的敏感性大于降溫度日，且采暖度日多為正向敏感性。但是分析1985—2017年的度日變化特征，可以發(fā)現(xiàn)，京津冀地區(qū)降溫度日的氣候傾向率為13.2 ℃·d/10 a，采暖度日的氣候傾向率為－68.4 ℃·d/10 a，但是采暖度日的多年平均值為3108.5 ℃·d，10 a的下降率只占到了2.2%，而降溫度日多年平均值為62.36 ℃·d，10 a的增加值占到了21.2%，即降溫度日的增加速率大于采暖度日的減少速率，所以綜合來看，降溫度日對這三個行業(yè)的產(chǎn)出影響迅速增大，且降溫度日每年波動比較大，其帶來的負面影響具有很大的不確定性。經(jīng)濟系統(tǒng)是復雜的，降溫度日的迅速增加在一定程度上反映了高溫熱浪天氣的增多，隨著全球氣候變暖趨勢加劇，對這三個行業(yè)的負面影響也不斷加強。對于建筑業(yè)和農(nóng)業(yè)，從目前來看，增暖對其是有利的，特別是建筑業(yè)，以戶外勞動為主，暖日數(shù)增加的同時，也增加了可作業(yè)時長，但是不能認為氣候變暖對其是持續(xù)有利的，很多研究認為氣候變暖與經(jīng)濟之間是非線性的關系，當氣候變暖到一定程度，負面效應將會占據(jù)主導作用。隨著變暖趨勢增強，對經(jīng)濟發(fā)展是否會產(chǎn)生抑制作用，仍需進一步研究。

氣象條件的變化和經(jīng)濟生產(chǎn)之間的關系非常復雜，可能簡單的線性關系和交互關系無法完整的定量分析它們之間的關系，但是，從生產(chǎn)函數(shù)的角度去探究，也是可行的，氣象要素的加入提高了方程的擬合能力。五個部門的R2都在90%以上，并且利用超越對數(shù)生產(chǎn)函數(shù)去擬合各行業(yè)最終產(chǎn)值，誤差位于1%～3%。

3 結論與討論

3.1 結論與建議

本文利用京津冀地區(qū)工業(yè)、建筑業(yè)等五個行業(yè)的生產(chǎn)總值、資本存量等經(jīng)濟數(shù)據(jù)和氣象數(shù)據(jù)，采用超越對數(shù)生產(chǎn)函數(shù)探究了氣象條件變化對京津冀地區(qū)行業(yè)的影響，主要結論如下。

1）采暖度日對建筑業(yè)、交通運輸倉儲郵政業(yè)、工業(yè)、批發(fā)零售業(yè)以及農(nóng)業(yè)的邊際產(chǎn)出彈性分別為－0.1501%、0.6197%、0.9309%、0.4473%以及－0.7311%，除了建筑業(yè)和農(nóng)業(yè)之外，都表現(xiàn)為正向影響。降溫度日對五個經(jīng)濟部門的邊際產(chǎn)出彈性分別為0.0274%、－0.0208%、－0.0709%、－0.1028%以及0.0282%，對建筑業(yè)和農(nóng)業(yè)具有微弱的正向作用。降水距平對五個部門的邊際影響分別為－0.1643%、－0.1365%、－0.0005%、－0.0407%以及0.0407%，除農(nóng)業(yè)之外，都表現(xiàn)為負向影響。

2）氣候變暖的負面效應逐漸增強。對于交通運輸倉儲郵政業(yè)、工業(yè)以及批發(fā)零售業(yè)，采暖度日對其影響為正向效應，但當前采暖度日為減少趨勢，所以在全球變暖的背景下，采暖度日對這三個行業(yè)的影響為負面影響；降溫度日對這三個行業(yè)的影響為負面影響。當前降溫度日的增加速率大于采暖度日的減少速率，所以降溫度日對這三個行業(yè)的產(chǎn)出影響增長迅速。盡管建筑業(yè)和農(nóng)業(yè)從目前來看，增暖對其有利的，但是隨著變暖趨勢增強，對其是否會產(chǎn)生抑制作用，仍需進一步關注。

3）降水對除農(nóng)業(yè)以外的其他部門都為負面影響。降水對建筑業(yè)和交通運輸業(yè)的影響最為明顯，這兩個行業(yè)都以戶外活動為主，在降水天氣下，對這兩個行業(yè)的抑制效應最為明顯。對批發(fā)零售業(yè)的負面影響是考慮到降水天氣下，對人們的外出購物影響非常大。

主動適應氣候變化是減少其帶來的風險的主要行動。在交通、建筑和工業(yè)行業(yè)，可以提前升級和改造基礎設施，提高防御高溫和極端降水帶來的破壞能力。在批發(fā)零售行業(yè)，可以提高對月、季及年的氣候預測能力，提前通過貿(mào)易等方式去規(guī)避可能帶來的氣候風險。在農(nóng)業(yè)行業(yè)，干旱和洪澇是造成農(nóng)業(yè)減產(chǎn)的重要因素，培育適合地區(qū)氣候特征的良種和興修水利設施，提高農(nóng)業(yè)生產(chǎn)者的氣候風險意識，可以有效地減少氣候變化帶來的風險。

3.2 討論

本文將氣象要素作為一種生產(chǎn)要素引入到超越對數(shù)生產(chǎn)函數(shù)中，充分考慮了氣象要素不僅影響產(chǎn)出而且還會影響其他投入要素，并在面板數(shù)據(jù)模型的基礎上進行擬合，盡量減少數(shù)據(jù)異方差、共線性等帶來的擬合偏差，得到了一些比較有意義的結果。但是本文受到統(tǒng)計數(shù)據(jù)較為粗糙的限制，改進的空間還非常大。而且經(jīng)濟生產(chǎn)活動具有很強的適應能力，很多生產(chǎn)活動會在下一個季度彌補上個季度由于天氣影響而帶來的損失，但由于部門產(chǎn)值等指標只在年度水平上統(tǒng)計，所以會對研究天氣與生產(chǎn)之間的關系帶來很多的不利影響。在當今時代，生產(chǎn)過程非常復雜，存在著很多的影響因素，很多因素很難定量化，本文將其他未知因素利用時間變量代替，可能欠妥，今后要對生產(chǎn)函數(shù)進一步完善。

由于部門數(shù)據(jù)在地級市難以收集到，只能采用省級數(shù)據(jù)，可能區(qū)域大小不對等，河北省相對北京、天津面積大很多，且經(jīng)濟活動過于分散，氣候差異也較大，所以本文盡量采用城市站點的氣象數(shù)據(jù)。經(jīng)濟活動大多數(shù)集中于城市，但是仍可能抹平了區(qū)域內部的差異，所以在今后研究中要縮小“點”研究區(qū)域和擴大整個“面”研究區(qū)域。

本文只是從研究框架和方法進行了一些探索，在未來的研究中要注重對高敏感性行業(yè)的深入研究，并結合未來氣候變化情景和社會經(jīng)濟情景，增強研究的預見性，從而能夠有效適應和減緩氣候變化帶來的影響。