• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    碳排放增加與氣溫變化統(tǒng)計因果關系的多重檢驗

    2013-04-29 01:27:44白萬平楊廣仁張學敏
    貴州財經(jīng)大學學報 2013年5期
    關鍵詞:同質性

    白萬平 楊廣仁 張學敏

    收稿日期:2013-05-12

    基金項目:

    作者簡介:白萬平(1963-),男,重慶人,博士,貴州財經(jīng)大學教授,研究方向為計量經(jīng)濟模型應用研究;楊廣仁(1978-),男,山東人,博士研究生,副教授,研究方向為非參數(shù)統(tǒng)計;張學敏(1975-),女,四川資陽人,經(jīng)濟師,人力資源管理師,國家二級企業(yè)培訓師,研究方向為人力資源管理。

    摘要: 在氣候變化問題上,“災難論”和“陰謀論”是由碳排放量增加與氣溫升高“因果論”和“非因果論”衍生的兩種主要觀點,但未見學者從統(tǒng)計因果關系檢驗角度深入探究。選取全球161年的時間序列數(shù)據(jù),運用線性和非線性等方法檢驗兩者之間的統(tǒng)計因果關系,結果表明碳排放增加是氣溫升高的原因,同時發(fā)現(xiàn)氣溫變化存在一個滯后期為29年的周期性變動;在面板數(shù)據(jù)因果檢驗中,從全球抽取20個代表性國家101年的數(shù)據(jù),檢驗結果仍然是碳排放增加是氣溫變化的原因。因此,有證據(jù)表明,碳排放增加是全球氣溫變化的統(tǒng)計原因。節(jié)能減排是共同責任,“陰謀論”不足為信。

    關鍵詞:碳排放量;氣溫變化;統(tǒng)計因果關系;PVAR模型;同質性;Hisao檢驗

    文章編號:2095-5960(2013)05-0046-06

    ;中圖分類號:C81

    ;文獻標識碼:A

    在碳排放增加與地球溫度升高之間關系問題上,學術界一直存在“因果論”和“非因果論”之爭,近年來又被演變成一種政治工具,對應為“災難論”和“陰謀論”兩種觀點,特別是美國、加拿大兩個發(fā)達國家相繼退出減排協(xié)定后,使問題變得更加撲朔迷離。誠然,地球氣溫變化是一個復雜過程,氣象科學家通過科學觀測、模擬等手段證實了變化的復雜性。但換一個角度思考,兩者所反映的是事物(現(xiàn)象)之間的因果關系,在獲取大樣本的觀察數(shù)據(jù)后,可以利用統(tǒng)計因果檢驗的方法,探討兩者之間的統(tǒng)計因果關系。遺憾的是,對于這樣一個重大問題,未見國內外學者的相關著述。通過運用多種現(xiàn)代統(tǒng)計因果關系檢驗方法探索兩者的關系,其一能填補兩者之間統(tǒng)計因果關系方面的缺失,為爭論提供統(tǒng)計證據(jù),其二拓展統(tǒng)計因果關系的應用范圍,這正是本文試圖作出的貢獻。

    本文安排如下:第一部分,統(tǒng)計因果關系檢驗進展回顧;第二部分,選擇161年的時間序列,檢驗全球碳排放增加與氣溫升高之間的統(tǒng)計因果關系;第三部分,在全球選擇20個樣本國家,從面板數(shù)據(jù)角度檢驗兩者之間的因果關系檢驗;最后,為本文結論。

    一、統(tǒng)計因果關系檢驗進展

    Granger(1969)[9]開創(chuàng)性地提出統(tǒng)計非因果關系檢驗方法后,經(jīng)過四十余年的發(fā)展,用統(tǒng)計方法檢驗變量因果關系的研究,已經(jīng)在時間序列向面板數(shù)據(jù)、線性向非線性、參數(shù)向非參數(shù)檢驗等多個維度上發(fā)展,成為科學檢驗變量間統(tǒng)計因果關系的必要手段和工具。

    從統(tǒng)計上講,假設序列Y和X都是平穩(wěn)的,我們說X不引起Y(即X不是Y的Granger原因)是指在對Y作預測時,假如由包含X過去值所得到的預測值并不優(yōu)于未含X過去值的預測?;蛘哒f,如果X不引起Y, X無助于對Y的預測。反之,則說X引起Y,即X是Y的Granger原因。對非平穩(wěn)時間序列,假設兩個序列之間有協(xié)整關系,可由Engle and Granger (1987)[7]單方程檢驗完成,或按Granger(1988)[10]兩步法得出誤差修正模型(ECM)處理,其中涉及滯后階數(shù)的敏感性問題,可由Hsiao (1981,1982)[12] [13]用Akaike (1970)[1]的最終預報誤差 (FPE) 結合Granger方法處理。Hsiao的方法還可以用于處理多變量線性因果關系及其最優(yōu)滯后階數(shù)的選擇上。但發(fā)展至今,非線性因果關系檢驗問題在實證中仍然未獲得滿意結果,盡管已經(jīng)出現(xiàn)了BDS(1992)[2]、HJ(1994) [11]

    理論上,時間序列中的雙變量線性與非線性因果檢驗方法較為系統(tǒng),包括多因一果和一因多果在內的多變量因果關系,以及面板數(shù)據(jù)變量因果關系的檢驗已經(jīng)獲得了一些解決的方法,但研究還不充分。如在檢驗多變量因果關系時,目前僅有的TY檢驗因只對少數(shù)具有特殊單整情形給出證明,尚不具備普適性;BWZ檢驗將BDS方法所做的推廣存在過度拒絕的傾向;同時多元非線性Granger因果檢驗的成果目前大多處在理論層面,在應用中有一定難度。在面板數(shù)據(jù)因果關系檢驗中,Hurlin等人的研究成果對于檢驗某些場合的因果關系奏效,但該方法還存在分類不完備,可能出現(xiàn)難以解釋的結果;Furkan的方法 在T較小或內生工具變量不易確定時用固定Granger因果估計存在固有的偏差且不可比較等問題,表明在統(tǒng)計因果關系理論有待發(fā)展。

    在實證研究中,統(tǒng)計因果關系在研究經(jīng)濟系統(tǒng)變量間關系時應用較多,但對其他系統(tǒng)內變量間關系,特別是爭論了一個多世紀的碳排放與氣候變化之間的統(tǒng)計因果關系,未發(fā)現(xiàn)有學者采用上述統(tǒng)計方法研究。因此,無論是擴展統(tǒng)計因果關系方法的應用范圍還是為一個重大關系提供證據(jù),都有必要從時間序列和面板數(shù)據(jù)兩個層面對此進行檢驗,為解釋兩者之間存在的關系提供佐證。

    二、全球氣溫變化與碳排放統(tǒng)計因果關系檢驗

    在以下的分析中,我們視序列的屬性和滿足的條件,采用適合的方法,由標準的Granger統(tǒng)計因果檢驗開始,層層推進,從線性、非線性時間序列到面板數(shù)據(jù)等多個維度對兩者之間的統(tǒng)計因果關系進行全面檢驗,獲得相應證據(jù)。

    變量說明和樣本選取。根據(jù)變量的代表性和數(shù)據(jù)的可得性,碳排放量和氣溫分別以xt,yt表示①①以每年消費化石能源的碳排放量表示年碳排放量,全球年均氣溫表示氣溫。數(shù)據(jù)來源:英國國家氣象臺、美國橡樹嶺實驗室。,樣本時間段為1850—2010年。在跨越三個世紀的161年中,地球年均氣溫最高的是1998年,達到0.529℃,最低的是1911年,為零下0.573℃,溫度的極差為1.112℃,超過1℃,但161年氣溫升高不到1℃。同期全球碳排放量快速增加,由1850年的54mtc增加到2010年的9167mtc,人均碳排放量也由0.05tc增長到2010年的1.34tc,增長28倍。

    單位根檢驗。對原序列采用6種方法檢驗后不能排除存在單位根的可能,進一步對兩個差分序列進行檢驗,結果如表1。從檢驗結果看,兩個差分序列平穩(wěn),不存在單位根。

    為了觀察不同滯后階數(shù)的反應,將滯后階數(shù)從1取到50,得出檢驗統(tǒng)計量對應的P值,在檢驗結果中,沒有一致證據(jù)表明氣溫升高與碳排放增加的因果關系。在碳排放增加與氣溫升高的檢驗中,滯后階數(shù)大于1后,沒有證據(jù)表明碳排放增加是氣溫升高的Granger原因。而滯后1階時,F(xiàn)檢驗值達到4.30431,P值為0.0397。顯然,在5%的顯著性水平下,拒絕原假設,即有統(tǒng)計證據(jù)表明,碳排放增加是氣溫升高的Granger原因。因果檢驗對滯后階數(shù)的敏感性,無法得出一致的證據(jù),需要進一步檢驗。兩個變量是一階單整的,接下來檢驗是否存在協(xié)整關系。利用Johansen檢驗,記Yt為碳排放、氣溫構成的二維列向量,Xt是確定性向量,分別在5種不同的設定進行檢驗。整理后的檢驗結果如表2所示。

    表 2Johansen協(xié)整檢驗

    數(shù)據(jù)趨勢無無線性拋物線

    檢驗類型無截距、無趨勢有截距、無趨勢有截距、無趨勢有截距、有趨勢有趨勢、有截距

    跡12210

    最大特征根數(shù)12210

    單純觀察檢驗結果,無法得出是否存在協(xié)整關系的結論。利用變量回歸,不能拒絕趨勢項冗余的假設,而趨勢項對判斷兩者之間是否存在協(xié)整關系至關重要。因此,用Johansen協(xié)整檢驗并不能確定兩者之間的長期穩(wěn)定因果關系,采用Engle-Granger和Phillips-Ouliaris檢驗,所得結果整理如表3:

    由此可知,E-G和PP兩種單方程協(xié)整檢驗都拒絕碳排放量與氣溫之間不存在協(xié)整關系的假設,可以認為在樣本期內,兩者之間存在長期穩(wěn)定的關系。

    由于Granger檢驗不能確定最佳滯后階數(shù),為提高精度,進一步采用Hsiao(1981)的方法:按1/5樣本選擇滯后階的原則,選取的最大滯后階數(shù)為32 ,計算出各滯后階數(shù)的FPE值,比較篩選的最小FPE=0.00838,相應的滯后階數(shù)m*=29,該滯后階數(shù)可解釋為地球氣溫大約在29年左右出現(xiàn)周期性波動的特性;固定m*=29,引入碳排放量,計算到32階的FPE值,得出最小FPE=0.008277,相應n*=1,由此確定的最佳滯后階數(shù)為1期,即上期碳排放量對本期氣溫影響最為顯著,該結果與通過Granger因果檢驗所得滯后階數(shù)一致。所以,采用線性Granger因果關系檢驗后的結論是:碳排放增加是氣溫升高的統(tǒng)計原因。

    上述檢驗的基本假設是兩者之間存在線性關系,如果兩者之間存在的關系不是線性的,則證據(jù)不足,需要進一步采用非線性統(tǒng)計檢驗方法檢驗。

    由兩個變量的特征,首先選擇非線性因果關系中的BDS檢驗方法。用VAR模型剔除主要線性關系,檢驗余下部分的獨立同分布特征。本問題中,所關注的問題是xt是否引起yt的變動,經(jīng)過甄別,無約束、滯后兩階向量自回歸模型是適宜模型,檢驗y關于x向量自回歸模型中的殘差,按照選取滯后階數(shù)的原則,選擇直到6階后的檢驗結果見表4。

    檢驗結果表明經(jīng)過線性關系過濾后的序列中已經(jīng)不含其他非線性因果關系,無需再采用其他非線性因果方法檢驗,支持碳排放增加是氣溫升高的統(tǒng)計原因這一結果。

    利用線性和非線性統(tǒng)計因果關系的多種方法檢驗,所得結論是碳排放增加是氣溫升高的統(tǒng)計原因,全球碳排放增加在一定程度上引起了全球氣溫升高。

    三、碳排放與氣溫面板數(shù)據(jù)因果關系檢驗

    對碳排放與氣溫變化時間序列統(tǒng)計因果關系的研究,只是從時間進程上為我們提供了碳排放增加引起氣溫升高的證據(jù)。由于各地區(qū)工業(yè)化程度不一、區(qū)位各異,導致碳排放分布極不均衡,各地緯度等不同,氣溫差異也較大,僅僅由時序獲取的證據(jù)略顯單薄。為此,我們在全球選擇樣本點,采用面板數(shù)據(jù)因果關系檢驗方法,做出進一步的論證。

    樣本選擇。綜合考慮地區(qū)工業(yè)化程度、區(qū)域分布均衡性和數(shù)據(jù)可得性,從全球選取20個國家①①各國氣溫以主要城市代替,德國、俄羅斯等有分合的國家數(shù)據(jù)已按統(tǒng)一口徑處理,數(shù)據(jù)來源同前。,分別為:美國、俄羅斯、英國、法國、德國、日本、波蘭、挪威、瑞典、奧地利、愛爾蘭、南非、墨西哥、印尼、印度、巴西、新西蘭、阿根廷、澳大利亞、中國。其中包含了早期工業(yè)化國家、主要大國和發(fā)展中國家,分布在不同緯度,樣本國家的總碳排放量占全球排放的近70%,樣本有較強的代表性。受觀測數(shù)據(jù)限制,為了兼顧各大洲皆有國家入選,最后確定樣本區(qū)間為1910年至2010年,采集了101年的面板數(shù)據(jù)。

    檢驗說明,對1、2、3階的滯后期數(shù),檢驗值都遠大于臨界值(顯著性水平為1%),拒絕原假設。因此,不能排除碳排放增加是氣溫變化的統(tǒng)計原因。

    至此,通過利用時間序列和面板數(shù)據(jù)的多種因果關系檢驗,結果都拒絕碳排放增加不是氣溫變化統(tǒng)計原因的假設。

    在現(xiàn)有研究成果中,基于各地區(qū)不同區(qū)位、氣溫和碳排放量,氣溫升高所產(chǎn)生的影響有別,“災難論”似乎證據(jù)更充分,人們總體上更傾向于氣溫升高對人類社會經(jīng)濟發(fā)展“弊大于利”的結論。假設氣溫持續(xù)升高,可以預見的危害包括海平面上升、病蟲害增加、氣候反常頻率增加和沙漠化加劇等,將直接挑戰(zhàn)可持續(xù)發(fā)展的基本條件。如果不能有效控制碳排放,阻止氣溫升高,對于經(jīng)濟發(fā)展主要依靠沿海地區(qū)的海岸國家,氣溫升高導致的經(jīng)濟損失后果將更嚴重??刂苹茉聪M、尋找清潔替代能源、有效減少碳排放從而實現(xiàn)社會經(jīng)濟可持續(xù)發(fā)展,理應成為人類的共同選擇和發(fā)展路徑。

    四、結論

    碳排放增加與氣候變化是復雜系統(tǒng)中的兩個方面,選取碳排放和氣溫兩個變量,在總結統(tǒng)計因果關系檢驗方法進展后,采用線性Granger因果檢驗捕捉到在碳排放增量滯后一期時,兩者存在因果關系。進一步采用Hisao檢驗法,得出最優(yōu)滯后階數(shù),表明碳排放增加是氣溫升高的原因,同時還發(fā)現(xiàn)地球氣溫在29年左右出現(xiàn)一個周期的結論??紤]兩者之間可能存在非線性因果關系,用非線性檢驗方法,同樣支持碳排放增加引起氣溫升高的結論。有證據(jù)表明,在全球范圍內,碳排放增加是氣溫升高的直接原因。為提供更充分的證據(jù),從全球選取20個樣本國家,收集從1910年到2010年共101年的數(shù)據(jù),采用面板數(shù)據(jù)因果檢驗法,在檢驗出這些國家的碳排放和氣溫具有同質性后,進一步檢驗兩者的因果關系,結果同樣支持碳排放增加是氣溫變化的統(tǒng)計原因。因此,從統(tǒng)計意義上講,不能拒絕碳排放增加是氣溫變化的原因。

    研究結論表明,按照無悔策略原則,中國政府在應對氣候變化問題上的一系列決策是合理的選擇。在資源趨緊和氣溫變化雙重約束下,各地區(qū)需要從實際出發(fā),選擇科學發(fā)展路徑,節(jié)能減排、保護環(huán)境,保持經(jīng)濟又好又快增長,這不僅是自身發(fā)展的需要,也是一個對世界負責的政府應有的態(tài)度。

    參考文獻

    [1]Akaike, H. (1970). Statistical predicator identification. Annals of the Institute of Statistical Mathematics[J], 21, 203–207.

    [2]Baek, E.G., Brock, W.A., 1992. A general test for nonlinear Granger causality: bivariate model[R]. Working Paper. Korea Development Institute, University of Wisconsin-Madison.

    [3]Cees Diks, Valentyn Panchenko, A new statistic and practical guidelines for nonparametric Granger causality testing[J], Journal of Economic Dynamics & Control 30 (2006) 1647–1669.

    [4]Douglas Holtz-Eakin, Whitney Newey and Harvey S. Rosen, Estimating vector autoregressions with panel data[J], Econometrica, 1988, 56( 6) :1371- 1395.

    [5]Dufour, J.-M., Pelletier, D., Renault, E, 2006. Short run and long run causality in time series: Inference[J]. Journal of Econometrics 132 (2), 337–362.

    [6]Dufour, Abderrahim Taamoutic, Short and long run causality measures: Theory and inference[J], Journal of Econometrics 154 (2010) 42–58.

    [7]Engle, R.F, Granger, C.W.J., 1987. Cointegration and error-correction: representation, estimation and testing[J]. Econometrica 55, 251–276.

    [8]Furkan Emirmahmutoglu ,Nezir KoseTesting for Granger causality in heterogeneous mixed panels[J],Economic Modelling,Volume 28, Issue 3, May 2011, Pages 870-876.

    [9]Granger, C.W.J, 1969. Investigating causal relations by econometric models and crossspectral methods[J]. Econometrica37, 424–438.

    [10]Granger, C.W.J., 1988. Causality, cointegration and control[J]. Journal of Economic Dynamics and Control 12, 551–559.

    [11]Hiemstra, C., & Jones, J., D. (1994). Testing for Linear and Non-linear Granger Causality in the Stock Price Volume Relation[J]. Journal of Finance, 49, 1639-1664.

    [12]Hsiao, C. (1981). Autoregressive modelling and money-income causality detection[J]. Journal of Monetary Economics, 7,85–106.

    [13]Hsiao, C. (1982) Autoregressive Modelling and Causal Ordering of Economic Variables[J].Journal of Economic Dynamics and Control 4: 243-259.

    [14]Hurlin Christophe, Venet Baptiste, Granger Causality Tests in Panel Data Models Fixed Coefficients[R] ,Working Paper, Eurisco, Université Paris IX Dauphine ,2001.

    [15]Hurlin Christophe, Venet Baptiste, Granger Causality Tests in Panel Data Models with Fixed Coefficients[R], Working Paper ,2003.

    [16]Hurlin C. Testing Granger causality in heterogenous panel data models with fixed coefficients[R]. Paris: Université Paris IX Dauphine,2004.

    [17]Hurlin, C., 2008. Testing for Granger non-causality in heterogeneous panels[R]. Working Papers halshs-00224434, version 1.

    [18]Jean-Marie Dufour, Abderrahim Taamoutic, Short and long run causality measures: Theory and inference[J], Journal of Econometrics 154 (2010) 42-58.

    [19]Johansen, S., Estimation and hypothesis testing of cointegration vectors in Gaussian vector autoregressive models[J]. Econometrica 59(1991) 1551–1580.

    [20]Konya, L., 2006. Exports and growth: Granger causality analysis on OECD countries with a panel data approach[J]. Economic Modeling 23, 978–992.

    [21]Michael Eichler,Granger causality and path diagrams for multivariate time series[J],Journal of Econometrics 137 (2007) 334–353.

    [22]P. Grassberger and I. Procaccia. Measuring the strangeness of strange attractors[J]. Physica D, 9:189--208, 1983.

    [23]Toda, H.Y., Yamamoto, T., 1995. Statistical inference in vector autoregression with possibly integrated processes[J]. Journal of Econometrics 66, 225–250.

    [24]Weinhold D, A dynamic‘fixed effectsmodel for heterogeneous panel data[R],London: London School of Economics, April 1999.

    [25]Zhidong Bai, Heng Li, Wing-Keung Wong, Bingzhi Zhang, Multivariate causality tests with simulation and application[J], Statistics and Probability Letters 81 (2011) 1063–1071.

    [26]Zhidong Bai, Wing-Keung Wong, Bingzhi Zhang,Multivariate linear and nonlinear causality tests[J],Mathematics and Computers in Simulation 81 (2010) 5–17.

    猜你喜歡
    同質性
    新媒體社群門檻設置的特點與功能
    新媒體研究(2019年2期)2019-03-29 12:05:58
    關于發(fā)揮高校學生群體同質性的積極作用探究
    卷宗(2019年4期)2019-03-18 11:44:28
    鄉(xiāng)土田園牧歌的同質性與異質性
    青年時代(2018年22期)2018-10-16 10:44:26
    汽車廣告的同質性與個性化
    高管團隊特征對企業(yè)績效的影響研究
    軟科學(2017年6期)2017-06-30 22:59:13
    媒介文化同質性和獨創(chuàng)性之間的矛盾審視
    魅力中國(2016年22期)2017-05-24 12:19:19
    同質性和社會影響對混合型社交網(wǎng)絡形成的仿真分析
    基于同質性審視的高職應用型本科工程教育研究
    職教論壇(2017年4期)2017-03-13 16:43:19
    理性程度的異質性:基于理論與實踐的考察
    高等工程教育與高等職業(yè)教育的同質性
    免费大片18禁| 人人妻人人澡欧美一区二区| 国产乱人视频| 精品一区二区三区av网在线观看| 在线免费观看不下载黄p国产 | 18禁裸乳无遮挡免费网站照片| 亚洲,欧美精品.| 婷婷精品国产亚洲av在线| 国产伦一二天堂av在线观看| 黄色丝袜av网址大全| 真人一进一出gif抽搐免费| 国产精品香港三级国产av潘金莲| 又黄又粗又硬又大视频| 亚洲五月天丁香| 十八禁网站免费在线| 搡女人真爽免费视频火全软件 | bbb黄色大片| 亚洲av熟女| 欧美一区二区国产精品久久精品| 啦啦啦免费观看视频1| 久久国产精品影院| 可以在线观看毛片的网站| av片东京热男人的天堂| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人av| 国产真人三级小视频在线观看| 久久久久久人人人人人| 欧美不卡视频在线免费观看| 99热精品在线国产| 精品日产1卡2卡| 一级作爱视频免费观看| 欧美激情在线99| 国产av不卡久久| 久久亚洲精品不卡| 999久久久精品免费观看国产| 国产熟女xx| 日本在线视频免费播放| 国产精品美女特级片免费视频播放器| 国产一区二区三区在线臀色熟女| 色吧在线观看| 亚洲精品在线观看二区| 国产精品98久久久久久宅男小说| 亚洲国产日韩欧美精品在线观看 | 伊人久久大香线蕉亚洲五| 亚洲精品一卡2卡三卡4卡5卡| 日日夜夜操网爽| 欧美国产日韩亚洲一区| 午夜激情欧美在线| 亚洲欧美日韩卡通动漫| 欧美性猛交黑人性爽| 日本精品一区二区三区蜜桃| 精品国产三级普通话版| 成人特级黄色片久久久久久久| 亚洲午夜理论影院| 黑人欧美特级aaaaaa片| 少妇人妻精品综合一区二区 | 精品无人区乱码1区二区| 免费看a级黄色片| 国产精品影院久久| 久久久久久久亚洲中文字幕 | av国产免费在线观看| 动漫黄色视频在线观看| 网址你懂的国产日韩在线| 国产午夜精品论理片| 亚洲在线自拍视频| www国产在线视频色| 真人一进一出gif抽搐免费| 18美女黄网站色大片免费观看| 精品久久久久久久毛片微露脸| 久久久久久久久久黄片| 精品国产三级普通话版| 午夜福利欧美成人| 免费看日本二区| 日韩av在线大香蕉| 一二三四社区在线视频社区8| 久久久久久久久大av| 亚洲成人精品中文字幕电影| 日韩免费av在线播放| 在线十欧美十亚洲十日本专区| 欧美绝顶高潮抽搐喷水| 搡老熟女国产l中国老女人| av专区在线播放| 亚洲精品乱码久久久v下载方式 | 欧洲精品卡2卡3卡4卡5卡区| 欧美色视频一区免费| 91在线精品国自产拍蜜月 | 精品一区二区三区av网在线观看| 婷婷丁香在线五月| 少妇人妻精品综合一区二区 | 中亚洲国语对白在线视频| 欧美黑人欧美精品刺激| 搡女人真爽免费视频火全软件 | 午夜精品久久久久久毛片777| 舔av片在线| 久久婷婷人人爽人人干人人爱| av在线蜜桃| 特级一级黄色大片| 国产色婷婷99| 亚洲第一电影网av| 亚洲一区高清亚洲精品| 天天躁日日操中文字幕| 90打野战视频偷拍视频| 亚洲乱码一区二区免费版| 日韩欧美 国产精品| 国产视频一区二区在线看| 亚洲国产日韩欧美精品在线观看 | 亚洲av成人不卡在线观看播放网| 国产精品乱码一区二三区的特点| 18+在线观看网站| 国产伦精品一区二区三区四那| 日韩欧美国产在线观看| 高清在线国产一区| 欧美日韩亚洲国产一区二区在线观看| 又爽又黄无遮挡网站| 亚洲人成伊人成综合网2020| 麻豆国产av国片精品| 午夜福利免费观看在线| 亚洲国产欧美网| 小蜜桃在线观看免费完整版高清| 亚洲欧美日韩东京热| 久久人人精品亚洲av| 亚洲精品在线美女| 国产午夜精品论理片| 99精品在免费线老司机午夜| 内地一区二区视频在线| 欧美乱妇无乱码| 18禁黄网站禁片午夜丰满| 少妇人妻精品综合一区二区 | 久久久久性生活片| 久久久国产成人免费| 好男人在线观看高清免费视频| 国产麻豆成人av免费视频| 成年人黄色毛片网站| 午夜老司机福利剧场| 夜夜夜夜夜久久久久| 麻豆国产av国片精品| 午夜福利视频1000在线观看| 观看美女的网站| 夜夜爽天天搞| 99riav亚洲国产免费| 国产精华一区二区三区| 日本一二三区视频观看| 老鸭窝网址在线观看| 99久久无色码亚洲精品果冻| 免费搜索国产男女视频| 久久久久九九精品影院| 51国产日韩欧美| 狂野欧美激情性xxxx| 欧美一级毛片孕妇| av黄色大香蕉| 中文字幕人妻熟人妻熟丝袜美 | 男人舔奶头视频| 少妇人妻一区二区三区视频| 日本 av在线| 国产欧美日韩一区二区三| 午夜免费男女啪啪视频观看 | 在线播放无遮挡| 亚洲人成网站高清观看| 午夜视频国产福利| 久久精品国产综合久久久| 特级一级黄色大片| 黄色成人免费大全| 久久精品91蜜桃| 伊人久久大香线蕉亚洲五| 少妇丰满av| 久久这里只有精品中国| 男女下面进入的视频免费午夜| 他把我摸到了高潮在线观看| 天堂√8在线中文| 综合色av麻豆| 午夜福利免费观看在线| а√天堂www在线а√下载| 性色avwww在线观看| 最近在线观看免费完整版| 成人av在线播放网站| 欧美丝袜亚洲另类 | 亚洲中文字幕一区二区三区有码在线看| 成人av一区二区三区在线看| 色哟哟哟哟哟哟| 在线视频色国产色| 久久中文看片网| 露出奶头的视频| 亚洲精华国产精华精| 一边摸一边抽搐一进一小说| 日韩欧美精品免费久久 | 深爱激情五月婷婷| 久久久久性生活片| 欧美xxxx黑人xx丫x性爽| 免费看十八禁软件| 亚洲国产色片| 97超级碰碰碰精品色视频在线观看| 非洲黑人性xxxx精品又粗又长| 美女 人体艺术 gogo| 老司机午夜福利在线观看视频| 精品人妻一区二区三区麻豆 | 香蕉丝袜av| 国产精品三级大全| www.999成人在线观看| 精品久久久久久久毛片微露脸| 欧美性感艳星| 欧美成人a在线观看| 两个人视频免费观看高清| 黄片大片在线免费观看| 51国产日韩欧美| 老司机福利观看| 99久久成人亚洲精品观看| 丁香六月欧美| 精品福利观看| 麻豆久久精品国产亚洲av| 国产野战对白在线观看| 国产欧美日韩一区二区三| 免费高清视频大片| 床上黄色一级片| 叶爱在线成人免费视频播放| 观看美女的网站| 欧美黑人巨大hd| 国产欧美日韩一区二区精品| 成年女人毛片免费观看观看9| 成熟少妇高潮喷水视频| 窝窝影院91人妻| 欧美最新免费一区二区三区 | 真实男女啪啪啪动态图| 国产免费男女视频| 色哟哟哟哟哟哟| 最新中文字幕久久久久| a级毛片a级免费在线| 日本一二三区视频观看| 成年女人毛片免费观看观看9| 怎么达到女性高潮| 国产精品98久久久久久宅男小说| av国产免费在线观看| 岛国视频午夜一区免费看| 久久久久久大精品| 18禁裸乳无遮挡免费网站照片| 精品日产1卡2卡| 亚洲中文日韩欧美视频| 国产黄片美女视频| 在线视频色国产色| 国产亚洲欧美在线一区二区| 日韩亚洲欧美综合| 蜜桃久久精品国产亚洲av| 可以在线观看毛片的网站| 网址你懂的国产日韩在线| 男人舔奶头视频| 国产欧美日韩精品亚洲av| 中国美女看黄片| 欧美中文日本在线观看视频| 午夜免费男女啪啪视频观看 | 国产视频一区二区在线看| 一级a爱片免费观看的视频| 99在线人妻在线中文字幕| 午夜福利在线观看吧| 91av网一区二区| 黄色女人牲交| 亚洲av一区综合| 一进一出抽搐动态| 国产高清三级在线| 久久欧美精品欧美久久欧美| 午夜福利欧美成人| 熟妇人妻久久中文字幕3abv| 成人性生交大片免费视频hd| 国产黄a三级三级三级人| 亚洲中文字幕日韩| www.熟女人妻精品国产| 别揉我奶头~嗯~啊~动态视频| 综合色av麻豆| 岛国视频午夜一区免费看| 国产精品一区二区三区四区久久| 国产高清激情床上av| 日日夜夜操网爽| 成人av在线播放网站| 免费看十八禁软件| 欧美乱妇无乱码| 午夜精品在线福利| 蜜桃亚洲精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 观看美女的网站| 青草久久国产| 久久精品国产99精品国产亚洲性色| 久久久久性生活片| 精品99又大又爽又粗少妇毛片 | 99在线视频只有这里精品首页| 91麻豆精品激情在线观看国产| 久久精品国产清高在天天线| 亚洲av日韩精品久久久久久密| 欧美成人a在线观看| 99精品久久久久人妻精品| 亚洲av电影不卡..在线观看| 深爱激情五月婷婷| 丰满的人妻完整版| 久久精品国产自在天天线| 日韩欧美三级三区| 亚洲熟妇熟女久久| 女同久久另类99精品国产91| 午夜两性在线视频| 亚洲精品亚洲一区二区| 午夜a级毛片| 丰满的人妻完整版| 精品国产三级普通话版| 天堂√8在线中文| 一边摸一边抽搐一进一小说| 黄色视频,在线免费观看| 亚洲精品国产精品久久久不卡| 俄罗斯特黄特色一大片| 久久九九热精品免费| 女警被强在线播放| 午夜激情欧美在线| 国产久久久一区二区三区| 国产激情欧美一区二区| 日本黄大片高清| 高清日韩中文字幕在线| 亚洲国产精品sss在线观看| 亚洲国产欧洲综合997久久,| 色吧在线观看| 国产成人系列免费观看| www日本在线高清视频| 国产在视频线在精品| 亚洲欧美日韩高清在线视频| 麻豆一二三区av精品| 变态另类成人亚洲欧美熟女| 久久精品国产自在天天线| 97超视频在线观看视频| 免费看美女性在线毛片视频| 99在线视频只有这里精品首页| 精品人妻一区二区三区麻豆 | 亚洲av免费高清在线观看| 久久伊人香网站| 搞女人的毛片| 麻豆久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 桃红色精品国产亚洲av| 午夜福利成人在线免费观看| 看片在线看免费视频| 国产精品久久久久久久电影 | 欧美日韩一级在线毛片| 日本在线视频免费播放| 亚洲欧美日韩东京热| 精品人妻偷拍中文字幕| 欧美黑人欧美精品刺激| 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 久久久精品大字幕| 婷婷精品国产亚洲av| av国产免费在线观看| 淫秽高清视频在线观看| 成人性生交大片免费视频hd| 黑人欧美特级aaaaaa片| 黄色日韩在线| 五月伊人婷婷丁香| 久久久精品欧美日韩精品| 精品一区二区三区视频在线 | av国产免费在线观看| 一个人免费在线观看电影| 真人一进一出gif抽搐免费| 欧洲精品卡2卡3卡4卡5卡区| 伊人久久精品亚洲午夜| 久久亚洲精品不卡| 日韩av在线大香蕉| 日本三级黄在线观看| 欧美日本视频| 91字幕亚洲| 长腿黑丝高跟| 俺也久久电影网| 观看免费一级毛片| 国内精品美女久久久久久| 丰满乱子伦码专区| 狂野欧美激情性xxxx| 两性午夜刺激爽爽歪歪视频在线观看| 国产av在哪里看| 搡老妇女老女人老熟妇| 日韩欧美在线二视频| 国产一区二区三区视频了| 午夜亚洲福利在线播放| 日韩欧美国产在线观看| 精品人妻偷拍中文字幕| 丰满乱子伦码专区| 日本一二三区视频观看| 非洲黑人性xxxx精品又粗又长| 亚洲男人的天堂狠狠| 国产亚洲精品综合一区在线观看| 久久人人精品亚洲av| 亚洲五月天丁香| 日日摸夜夜添夜夜添小说| 久久国产乱子伦精品免费另类| 日韩高清综合在线| 国产亚洲精品av在线| 美女大奶头视频| 18禁美女被吸乳视频| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 免费搜索国产男女视频| 亚洲中文日韩欧美视频| 国产精品香港三级国产av潘金莲| 亚洲最大成人手机在线| 国产精品一区二区三区四区免费观看 | 欧美高清成人免费视频www| 久久精品国产清高在天天线| 亚洲第一欧美日韩一区二区三区| 麻豆久久精品国产亚洲av| 男女之事视频高清在线观看| 欧美大码av| 又粗又爽又猛毛片免费看| 91在线精品国自产拍蜜月 | 99国产极品粉嫩在线观看| 偷拍熟女少妇极品色| 国产v大片淫在线免费观看| 亚洲av二区三区四区| 国产私拍福利视频在线观看| 久久久国产成人免费| 小蜜桃在线观看免费完整版高清| 日韩欧美国产一区二区入口| 免费搜索国产男女视频| 丰满人妻一区二区三区视频av | 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| av在线天堂中文字幕| 首页视频小说图片口味搜索| 熟女电影av网| 欧美性猛交╳xxx乱大交人| 国产 一区 欧美 日韩| 在线观看av片永久免费下载| 最近视频中文字幕2019在线8| 久久精品国产99精品国产亚洲性色| 在线观看免费视频日本深夜| 国产一区二区三区在线臀色熟女| 精品人妻一区二区三区麻豆 | 免费av不卡在线播放| 一二三四社区在线视频社区8| www.www免费av| 久久人妻av系列| 最好的美女福利视频网| 一区二区三区免费毛片| 91九色精品人成在线观看| 亚洲精品色激情综合| 99久久综合精品五月天人人| av在线蜜桃| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av| 国产免费av片在线观看野外av| 午夜免费激情av| 国产欧美日韩精品一区二区| 亚洲欧美日韩无卡精品| 亚洲精品影视一区二区三区av| 免费在线观看影片大全网站| 亚洲国产精品999在线| 免费人成视频x8x8入口观看| 蜜桃亚洲精品一区二区三区| 99视频精品全部免费 在线| a级毛片a级免费在线| 黄色丝袜av网址大全| 国产高潮美女av| 欧美zozozo另类| 人人妻人人看人人澡| 欧美在线黄色| 国产亚洲av嫩草精品影院| 人人妻人人看人人澡| 美女黄网站色视频| 国产v大片淫在线免费观看| 91久久精品电影网| 日韩欧美精品免费久久 | 制服人妻中文乱码| 在线观看66精品国产| 18禁黄网站禁片免费观看直播| 午夜福利免费观看在线| 夜夜夜夜夜久久久久| 身体一侧抽搐| 91久久精品国产一区二区成人 | 每晚都被弄得嗷嗷叫到高潮| 欧美在线一区亚洲| 制服人妻中文乱码| 国产三级中文精品| 最近最新免费中文字幕在线| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 日本撒尿小便嘘嘘汇集6| 一个人看视频在线观看www免费 | 黄色丝袜av网址大全| 少妇丰满av| 欧美中文日本在线观看视频| 此物有八面人人有两片| 久久久国产成人免费| 国产极品精品免费视频能看的| 午夜精品在线福利| 1000部很黄的大片| 毛片女人毛片| 免费无遮挡裸体视频| 成人午夜高清在线视频| 国产色爽女视频免费观看| 亚洲电影在线观看av| 欧美zozozo另类| 国产伦一二天堂av在线观看| 久久久久久国产a免费观看| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇仑乱视频hdxx| 啪啪无遮挡十八禁网站| 亚洲最大成人中文| 国产精品女同一区二区软件 | 99国产综合亚洲精品| 窝窝影院91人妻| 九九久久精品国产亚洲av麻豆| 午夜激情福利司机影院| 少妇的逼水好多| 色综合亚洲欧美另类图片| 中出人妻视频一区二区| 中文字幕av成人在线电影| 丝袜美腿在线中文| 天堂av国产一区二区熟女人妻| 国产欧美日韩一区二区三| 国产精品久久电影中文字幕| 97人妻精品一区二区三区麻豆| 免费在线观看亚洲国产| 伊人久久大香线蕉亚洲五| 一级作爱视频免费观看| 欧美不卡视频在线免费观看| 中文字幕熟女人妻在线| 国产激情偷乱视频一区二区| 长腿黑丝高跟| 观看免费一级毛片| 欧美午夜高清在线| 日韩欧美在线二视频| 91av网一区二区| 国产亚洲欧美98| 神马国产精品三级电影在线观看| 精品熟女少妇八av免费久了| 亚洲专区国产一区二区| 琪琪午夜伦伦电影理论片6080| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 国产美女午夜福利| 亚洲精品在线美女| 亚洲五月天丁香| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品日韩av在线免费观看| 真实男女啪啪啪动态图| 欧美av亚洲av综合av国产av| 成人鲁丝片一二三区免费| 999久久久精品免费观看国产| 偷拍熟女少妇极品色| 级片在线观看| 最好的美女福利视频网| 一级作爱视频免费观看| 国产精品99久久99久久久不卡| 亚洲av一区综合| 久久久久久久久中文| 99精品欧美一区二区三区四区| 午夜福利成人在线免费观看| 人人妻,人人澡人人爽秒播| 男人和女人高潮做爰伦理| 国产乱人伦免费视频| 51午夜福利影视在线观看| 午夜免费男女啪啪视频观看 | 国产三级中文精品| 午夜精品久久久久久毛片777| 一级黄色大片毛片| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 少妇的逼好多水| 国产在线精品亚洲第一网站| 国产精品久久久久久精品电影| 成人亚洲精品av一区二区| 日韩成人在线观看一区二区三区| 可以在线观看的亚洲视频| 成人特级黄色片久久久久久久| 男女之事视频高清在线观看| 在线播放无遮挡| 国产国拍精品亚洲av在线观看 | 国产在视频线在精品| 亚洲精品一区av在线观看| 国产黄色小视频在线观看| 在线国产一区二区在线| 亚洲精品色激情综合| 最新美女视频免费是黄的| 久久这里只有精品中国| 欧美另类亚洲清纯唯美| 欧美av亚洲av综合av国产av| 国产97色在线日韩免费| 国产男靠女视频免费网站| 成人鲁丝片一二三区免费| 日韩精品青青久久久久久| 欧美色视频一区免费| 亚洲一区二区三区色噜噜| 久99久视频精品免费| 亚洲精品在线观看二区| 久久久国产成人精品二区| 国产不卡一卡二| 色尼玛亚洲综合影院| 午夜免费观看网址| 国产成人a区在线观看| 亚洲国产欧洲综合997久久,| 国内精品久久久久精免费| 18禁国产床啪视频网站| 亚洲天堂国产精品一区在线| 国产私拍福利视频在线观看| 国产久久久一区二区三区| 丰满的人妻完整版| 久久久久久久亚洲中文字幕 | 欧美黑人欧美精品刺激| av专区在线播放| 亚洲狠狠婷婷综合久久图片| 99国产精品一区二区三区| 久久亚洲真实| 九九在线视频观看精品| 十八禁人妻一区二区| 不卡一级毛片| 操出白浆在线播放| 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 国产成人系列免费观看| 亚洲不卡免费看| 黄片小视频在线播放| 免费一级毛片在线播放高清视频| 舔av片在线| 国产免费男女视频| 免费搜索国产男女视频| 在线观看66精品国产|