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

    地極移動(dòng)與河川徑流的關(guān)系研究

    2014-03-02 02:47:37劉蘇峽王盛王月玲李明新黃玫彭公炳肖子牛
    關(guān)鍵詞:格蘭杰因果關(guān)系徑流

    劉蘇峽王盛,王月玲李明新黃玫彭公炳肖子牛

    (1 中國(guó)科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室,北京 100101;2 中國(guó)科學(xué)院大學(xué),北京 100049;3 長(zhǎng)江水利委員會(huì)水文局,武漢 430010;4 中國(guó)氣象局氣象干部培訓(xùn)學(xué)院,北京 100081)

    地極移動(dòng)與河川徑流的關(guān)系研究

    劉蘇峽1王盛1,2王月玲1李明新3黃玫1彭公炳1肖子牛4

    (1 中國(guó)科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室,北京 100101;2 中國(guó)科學(xué)院大學(xué),北京 100049;3 長(zhǎng)江水利委員會(huì)水文局,武漢 430010;4 中國(guó)氣象局氣象干部培訓(xùn)學(xué)院,北京 100081)

    在總結(jié)地極移動(dòng)(以下簡(jiǎn)稱極移)和徑流之間可能存在的相互作用機(jī)制的基礎(chǔ)上,利用受人類活動(dòng)影響較少的青藏高原雅魯藏布江的月徑流資料和極移資料以及格蘭杰因果關(guān)系檢驗(yàn)方法,從統(tǒng)計(jì)學(xué)角度探索了在月、季和年尺度上,極移變化與徑流變化之間可能存在的聯(lián)系。結(jié)果顯示,在月尺度上,極移X分量變化量在滯后1~8個(gè)月和10~21個(gè)月后是徑流變化量的格蘭杰原因;極移Y分量變化量在滯后1~9個(gè)月和17~24個(gè)月后是徑流變化量的格蘭杰原因。在季尺度上,極移X分量變化量只有在滯后2個(gè)季的情況下是徑流變化量的格蘭杰原因;而極移Y分量變化量在滯后4~6個(gè)季的情況下是徑流變化量的格蘭杰原因。在年尺度上,未檢測(cè)出格蘭杰因果關(guān)系。從徑流到極移的檢測(cè)中發(fā)現(xiàn),月尺度上,徑流變化量在滯后3~25個(gè)月后為極移X分量變化量的格蘭杰原因;徑流變化量在滯后1個(gè)月和3~25個(gè)月后為極移Y分量變化量的格蘭杰原因。在季尺度上,徑流變化量在滯后2~8個(gè)季的情況下是極移X分量變化量的格蘭杰原因;徑流變化量在滯后1~8個(gè)季后是極移Y分量變化量的格蘭杰原因。在年尺度上,未發(fā)現(xiàn)徑流變化量和極移變化量存在格蘭杰因果關(guān)系。在月、季和年尺度上,極移變化量和徑流變化量的不同格蘭杰因果關(guān)系表明,利用極移資料可能在月和季尺度上提高資料稀缺區(qū)域的水文預(yù)測(cè)精度。

    地極移動(dòng)(極移),徑流,格蘭杰因果關(guān)系,雅魯藏布江

    1 引言

    地極移動(dòng)(polar motion),簡(jiǎn)稱極移,是指地球自轉(zhuǎn)軸相對(duì)于地球本體的位置變化。影響地球自轉(zhuǎn)軸變化的力包括外力和內(nèi)力。因作用機(jī)制的不同,地軸運(yùn)動(dòng)的表現(xiàn)形式也不同,可分為兩類。第一類是被稱為歲差(precession)和章動(dòng)(nutation)的地軸的方向變動(dòng)。它們由不涉及機(jī)械能消長(zhǎng)(如日月潮)的外力所驅(qū)動(dòng)。展現(xiàn)在太空,歲差相當(dāng)于轉(zhuǎn)軸方向在地球頭上畫圓的運(yùn)動(dòng),周期約2.5萬年(常被稱為柏拉圖年,或者大年)。章動(dòng)相當(dāng)于轉(zhuǎn)軸方向沿著歲差的軌跡的微小顫動(dòng),周期約為18.6年。另一類即為極移,主要由內(nèi)力驅(qū)動(dòng),是三維物體轉(zhuǎn)動(dòng)時(shí)自發(fā)的擺動(dòng)(wobble)。目前已觀測(cè)到的1年和14個(gè)月的錢德勒(Chandler)周期。驅(qū)動(dòng)極移的內(nèi)力可能是大氣環(huán)流、雨雪消長(zhǎng)、洋流和地下水等,但一直還沒有定論[1]。 不過Chen等[2]最新的研究成果展示了北極海冰的消融對(duì)近10年的地軸向東移動(dòng)的驅(qū)動(dòng)聯(lián)系。

    從機(jī)理上講,極移和徑流之間可能存在兩種內(nèi)在聯(lián)系。一方面如圖1所示,極移可以通過變形力的作用影響地球表層的大氣運(yùn)動(dòng),進(jìn)而影響陸地水文循環(huán)過程[3-4]。彭公炳等[3]基于1900—1977年的極移資料和我國(guó)一些地區(qū)的氣溫、降水、氣壓、少數(shù)水位和流量資料以及北半球內(nèi)的一些大氣環(huán)流和大氣活動(dòng)中心指標(biāo),采用滑動(dòng)分析、方差分析和能譜分析,發(fā)現(xiàn)極移的12個(gè)月、14個(gè)月、6~7年和35年的周期在氣候上都有明顯反映。年周期一致不言而喻,而北京、沈陽、長(zhǎng)春、哈爾濱、鄭州、重慶、昆明、汕頭8個(gè)站10年逐月平均氣溫、500hPa副高指數(shù)、亞洲經(jīng)向環(huán)流指數(shù)以及武漢、九江、岳陽逐月平均降水都具有11.4和13.3個(gè)月的周期,即接近極移的12和14個(gè)月周期。而東北3站(哈爾濱、長(zhǎng)春、沈陽)1—2月平均氣溫和年絕對(duì)最低氣溫,長(zhǎng)江中下游5站(上海、南京、九江、蕪湖、武漢)1—2月絕對(duì)最低氣溫、6—8月絕對(duì)最低氣溫和年絕對(duì)氣溫,華南4站(廣州、南寧、貴陽、昆明)年絕對(duì)最低氣溫,長(zhǎng)江中下游5站5—8月降水,華北5站(北京、天津、保定、石家莊、營(yíng)口)7—8月降水,漢口站6—9月平均流量,太平洋副高西界和南界,赤道低壓位置,都有6~7年的周期。長(zhǎng)江中下游5站5—8月降水和年降水、華北5站年降水、東北3站年絕對(duì)最低氣溫、華南5站1—2月絕對(duì)最低氣溫、長(zhǎng)江5站6—8月平均氣溫、赤道低壓東界和西伯利亞高壓北界、7月印度低壓東界位置都有約30~36年的周期。

    另一方面,如圖2所示,徑流變化反映了區(qū)域水儲(chǔ)量信息,而水儲(chǔ)量的變化引起地球質(zhì)量再分配,對(duì)極移具有重要的激發(fā)作用[4-10]。

    在極移與徑流之間存在的上述可能的緊密聯(lián)系基礎(chǔ)上,本文運(yùn)用格蘭杰(Granger)因果關(guān)系檢驗(yàn)方法,從統(tǒng)計(jì)學(xué)角度進(jìn)一步探索了徑流和極移之間的關(guān)系。

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

    2.1 格蘭杰方法

    格蘭杰因果關(guān)系[12]是從統(tǒng)計(jì)學(xué)角度探索兩組時(shí)間序列之間的因果關(guān)系的方法,起初被用于經(jīng)濟(jì)學(xué),目前已開始被用于氣象和水文科學(xué)研究中[13-16]。對(duì)于變量Y的預(yù)測(cè),如果在加入變量X的信息后能夠顯著地提高單獨(dú)運(yùn)用Y的信息進(jìn)行預(yù)測(cè)的效果,則稱變量X是變量Y的格蘭杰原因。定義k階自回歸模型(AR)和k階

    的向量自回歸模型(VAR)分別為:

    圖1 極移對(duì)徑流影響的可能物理機(jī)制(圖中問號(hào)表示從極移變性力到大氣環(huán)流的機(jī)制迄今還不甚明了,極移變化曲線引自文獻(xiàn)[11])

    圖2 徑流對(duì)極移影響的可能物理機(jī)制

    式中,C1和C2為常量,α、β、γ是模型的系數(shù),εt0和εt1是模型的殘差。格蘭杰檢驗(yàn)的原假設(shè)(二者之間不存在因果關(guān)系):

    通過最小二乘法計(jì)算AR和VAR模型的參數(shù),然后用F檢驗(yàn)判定兩個(gè)模型之間是否存在顯著差異。

    式中,RSSAR是AR模型的殘差平方和,RSSVAR是VAR模型的殘差平方和,q是AR模型的系數(shù)的個(gè)數(shù)(q=k),m是VAR模型的系數(shù)的個(gè)數(shù)(m=2k+1),n是觀測(cè)數(shù)據(jù)樣本數(shù)。

    格蘭杰因果關(guān)系檢驗(yàn)的前提是序列必須是平穩(wěn)的,否則將出現(xiàn)假的因果關(guān)系[17-18]。如果序列不平穩(wěn),可通過進(jìn)行一階差分甚至是二階差分計(jì)算,把序列轉(zhuǎn)化為平穩(wěn)序列。單位根檢驗(yàn)(Augmented Dickey Fuller,ADF)[19]常被用來進(jìn)行序列平穩(wěn)性檢驗(yàn)[15]。

    2.2 數(shù)據(jù)

    極移資料來自國(guó)際地球自轉(zhuǎn)和參考系服務(wù)(International Earth Rotation and Reference Systems Service)的地球定位參數(shù)(The Earth Orientation Parameters)計(jì)劃。相對(duì)于國(guó)際協(xié)議原點(diǎn)(Conventional International Origin,CIO),地極用坐標(biāo)X和Y來記錄,單位為毫弧秒(mas)。X方向?yàn)镚reenwich子午線方向,從北極向南為正;Y方向?yàn)槲鹘?jīng)90°方向,從北極向南為正。資料起始年限為1846年,原始資料時(shí)間間隔包括0.1年和0.5年。

    鑒于全球大部分河流受人類活動(dòng)影響,因此選取位于青藏高原的受人類活動(dòng)影響小、區(qū)域水循環(huán)變化主要較多地受控于自然因素的雅魯藏布江(圖3),進(jìn)行極移與徑流之間的關(guān)系探索。徑流序列是1978年1月—2006年12月的奴下水文站的月徑流資料(Q),單位為m3·s-1。奴下水文站流域控制面積為191235km2。在徑流的組成成分里,冰川積雪融水占38%[20],徑流變化中包含了大量的冰川積雪變化和土壤水分變化信息。

    2.3 分析方法

    首先對(duì)時(shí)間序列X、Y和Q進(jìn)行平穩(wěn)性檢驗(yàn),以p=0.05為顯著性水平。然后從月、季和年3個(gè)不同尺度進(jìn)行格蘭杰因果檢驗(yàn),探索不同時(shí)間尺度上,極移和徑流之間可能存在的關(guān)系。為了跟徑流資料匹配,檢驗(yàn)所用資料年限統(tǒng)一為 1978年1月—2006年12月,極移資料統(tǒng)一處理為月尺度間隔。為了探索極移和徑流二者之間的關(guān)系,分別進(jìn)行從極移到徑流(徑流預(yù)測(cè)中加極移)和從徑流到極移(極移預(yù)測(cè)中加徑流)的格蘭杰因果關(guān)系檢驗(yàn)。為了兼顧模型的預(yù)測(cè)效率和自由度,在選取滯后項(xiàng)的時(shí)候,月尺度上進(jìn)行了1~36個(gè)滯后項(xiàng)檢驗(yàn),季尺度上進(jìn)行了1~12個(gè)滯后項(xiàng)檢驗(yàn),年尺度上進(jìn)行了1~3個(gè)滯后項(xiàng)檢驗(yàn)。

    圖3 雅雅魯藏布江及奴下水文站示意圖

    3 結(jié)果

    3.1 極移和徑流的時(shí)間變化趨勢(shì)

    圖4為極移X分量、Y分量和雅魯藏布江奴下站徑流在1978—2006年的變化趨勢(shì)??梢钥闯?在29年間,極移X分量沿Greenwich子午線方向相對(duì)于CIO南北來回?cái)[動(dòng),變幅在―250~300mas之間,1mas約為3cm,故變程大約為16m。極移Y分量沿西經(jīng)90°向南方向擺動(dòng)式漂移,變幅在 0~600mas之間,變程合約18m。

    雅魯藏布江奴下站徑流29年間變化表現(xiàn)為前20年平穩(wěn),后9年水量增大。在明顯的年內(nèi)季節(jié)變化的基礎(chǔ)上,每年的月徑流峰值相差8000m3·s-1,折合成徑流深,峰值相差0.108m。

    29年間,相對(duì)于徑流的29個(gè)峰值,極移X分量有24個(gè)峰值,極移Y分量有25個(gè)峰值。極移X和Y分量為正弦式脈沖,徑流Q是蛙跳式脈沖。

    從上述直觀分析可見,極移和徑流之間存在一定的聯(lián)系,下面詳細(xì)探討存在的聯(lián)系。

    3.2 平穩(wěn)性檢驗(yàn)

    對(duì)3個(gè)時(shí)間尺度序列進(jìn)行ADF單位根平穩(wěn)性檢驗(yàn)(表1—3)。在進(jìn)行ADF檢驗(yàn)時(shí),需要確定所構(gòu)建的統(tǒng)計(jì)量的一階滯后差分和的項(xiàng)數(shù)p*。根據(jù)AIC準(zhǔn)則(Akaike information criteria)[21],月、季和年尺度上

    p*的最大值分別取30、12和6。

    由表1—3可知,在置信度為0.05的水平上,系列X、Y和Q均是不平穩(wěn)的(概率大于5%)。對(duì)X、Y和Q分別求一階差分,得到序列ΔX、ΔY和ΔQ,即極移X、Y分量和徑流的變化量,再進(jìn)行平穩(wěn)性檢驗(yàn)(表1—3),3個(gè)序列均轉(zhuǎn)化為了平穩(wěn)序列,因此可以對(duì)它們進(jìn)行格蘭杰因果關(guān)系檢驗(yàn)。

    圖4 1978年1月—2006年12月的極移X分量、Y分量和雅魯藏布江奴下站徑流

    表1 月尺度序列ADF檢驗(yàn)結(jié)果(p*的最大值取為30,顯著性水平為0.05)

    表2 季尺度序列ADF檢驗(yàn)結(jié)果(p*的最大值取為12,顯著性水平為0.05)

    表3 年尺度序列ADF檢驗(yàn)結(jié)果(p*的最大值取為6,顯著性水平為0.05)

    3.3 徑流預(yù)測(cè)中加極移因子對(duì)其預(yù)測(cè)精度提高的效果檢驗(yàn)

    圖5顯示了徑流變化預(yù)測(cè)中加極移因子的變化量對(duì)預(yù)測(cè)精度提高的效果。根據(jù)格蘭杰檢驗(yàn)原理,若拒絕原假設(shè)的概率比顯著性水平(0.05)還小,則拒絕原假設(shè)。具體到圖5a,也就是說,若建立滯后分別為1~8個(gè)月、10~21個(gè)月的徑流變化量的AR和VAR模型,在加入ΔX(極移X分量的變化量)后,兩模型之間存在顯著差異,徑流預(yù)測(cè)精度可能提高。同理,若建立滯后為1~9個(gè)月、17~24個(gè)月的徑流變化的AR和VAR模型,在加入極移Y分量的變化量后,兩模型之間存在顯著差異,徑流預(yù)測(cè)精度可能提高(圖5b)。

    在季尺度上,僅僅在建立滯后為2個(gè)季度的徑流變化的AR和VAR模型中,在加入極移X分量的變化量后,兩模型之間存在顯著差異,徑流預(yù)測(cè)精度可能提高(圖5c)。而在建立滯后為1、4到6個(gè)季度的徑

    流變化的AR和VAR模型,在加入極移Y分量的變化量后,兩模型之間存在顯著差異,徑流預(yù)測(cè)精度可能提高(圖5d)。

    圖5 極移影響徑流的因果關(guān)系檢驗(yàn)結(jié)果(a,c,e對(duì)應(yīng)加入極移X分量的變化量結(jié)果;b,d,f對(duì)應(yīng)加入極移Y分量的變化量結(jié)果;每個(gè)圖的y軸是拒絕原假設(shè)的概率,虛線對(duì)應(yīng)顯著性水平p=0.05)

    在年尺度上,不論是滯后多少年,加入極移X或者Y分量的變化量,徑流變化的AR和VAR模型之間不存在顯著差異,徑流預(yù)測(cè)精度不可能提高。

    3.4 極移預(yù)測(cè)中加入徑流因子對(duì)其預(yù)測(cè)精度提高的效果檢驗(yàn)

    圖6顯示了極移變化預(yù)測(cè)中加徑流變化量對(duì)預(yù)測(cè)精度提高的效果。對(duì)于月尺度,若建立滯后為3~25個(gè)月的極移X分量變化量的AR和VAR模型,在加入徑流變化量后,兩模型之間存在顯著差異,極移預(yù)測(cè)精度可能提高(圖6a)。

    若建立滯后分別為1個(gè)月和3~25個(gè)月的極移Y分量變化量的AR和VAR模型,在加入徑流變化量后,兩模型之間存在顯著差異,極移變化量預(yù)測(cè)精度可能提高(圖6b)。

    在季尺度上,在建立滯后為2~8個(gè)季度的極移X分量變化量的AR和VAR模型中,在加入徑流變化量后,兩模型之間存在顯著差異,極移變化量預(yù)測(cè)精度可能提高(圖6c)。在建立滯后為1~8個(gè)季度的極移Y分量變化量的AR和VAR模型中,在加入徑流變化量后,兩模型之間存在顯著差異,極移預(yù)測(cè)精度可能提高(圖6d)。

    在年尺度上,不論滯后多少年,加入徑流變化量,不論是對(duì)極移X還是Y分量的變化量,其AR和VAR模型之間不存在顯著差異,極移預(yù)測(cè)精度不可能提高。

    圖6 徑流影響極移的因果關(guān)系檢驗(yàn)結(jié)果(a,c,e對(duì)應(yīng)極移X分量變化量的結(jié)果;b,d,f對(duì)應(yīng)極移Y分量變化量的結(jié)果;每個(gè)圖的y軸是拒絕原假設(shè)的概率,虛線對(duì)應(yīng)顯著性水平p=0.05)

    4 結(jié)果討論

    4.1 徑流和極移究竟在哪個(gè)尺度更可能存在聯(lián)系

    從以上的月、季和年尺度上分別進(jìn)行的極移與徑流之間的格蘭杰因果關(guān)系探索中可知,在月和季節(jié)尺度上,徑流和極移在0.05顯著性水平上顯著地互為格蘭杰因果。也就是說,總能找到一個(gè)滯后時(shí)間并在該時(shí)間上建立AR和VAR模型,使得兩模型呈顯著差異,從而能提供改進(jìn)預(yù)測(cè)精度的契機(jī)。而在年尺度上,尚未檢測(cè)出因果關(guān)系。這驗(yàn)證了本文總結(jié)的極移和徑流之間可能存在聯(lián)系,并和Chen等[6]認(rèn)為陸地水儲(chǔ)量變化對(duì)極移具有很強(qiáng)的季節(jié)激發(fā)結(jié)論相一致。而這暗示了利用長(zhǎng)序列的極移,可能在月和季尺度上提高水文預(yù)報(bào)。

    4.2 是極移影響徑流還是徑流影響極移

    如前所述,極移和徑流誰是因,誰是果,尚不明確。從上述的月、季和年尺度上分別進(jìn)行的極移與徑流之間的格蘭杰因果關(guān)系探索中也可看出,徑流變化和極移變化互為影響。有意思的是不同的滯后尺度影響不一樣。在相同的滯后項(xiàng)選取的前提下,統(tǒng)計(jì)了徑流變化預(yù)測(cè)加入極移變化量和極移變化量預(yù)測(cè)加入徑流變化量,發(fā)生拒絕原假設(shè)的滯后模型的個(gè)數(shù)(表4)。發(fā)現(xiàn),單純從這個(gè)意義上,似乎是極移變化的加入對(duì)徑流變化預(yù)測(cè)精度的影響范圍要比徑流變化的加入對(duì)極移變化預(yù)測(cè)精度的影響范圍大(滯后模型個(gè)數(shù)多)。這暗示,極移對(duì)徑流的影響作用大于徑流對(duì)極移的影響作用。如果把徑流作為地球氣候系統(tǒng)的一個(gè)組成單元,也許這是地球運(yùn)動(dòng)因子影響氣候變化的一個(gè)潛在信號(hào)。

    如圖5—6以及表4所示,極移X和Y分量的變化量與徑流變化之間的關(guān)系在時(shí)間上并不完全同步:在同一滯后時(shí)間下拒絕原假設(shè)的統(tǒng)計(jì)概率有所差異,拒絕原假設(shè)的滯后模型的個(gè)數(shù)也不相同??傮w上來說,徑流變化與極移Y分量變化量的統(tǒng)計(jì)關(guān)系要強(qiáng)于徑流變化與極移X分量變化量的關(guān)系。這可能與雅魯藏布江所在的位置和該區(qū)域的水循環(huán)特點(diǎn)有著緊密聯(lián)系。如圖7所示,雅魯藏布江所在區(qū)域位于90°E附近,即極移Y分量的負(fù)方向,該區(qū)域的水文循環(huán)過程主要是印度洋季風(fēng)沿著經(jīng)向北上形成降水,降水再主要通過徑

    流的形式又回到印度洋,這個(gè)水文循環(huán)過程主要是經(jīng)向上的,所以在統(tǒng)計(jì)關(guān)系上該區(qū)域的徑流變化與極移Y分量的變化統(tǒng)計(jì)關(guān)系更顯著。

    表4 拒絕原假設(shè)的滯后模型的個(gè)數(shù)

    4.3 雅魯藏布江徑流與極移

    和整個(gè)地球相比,雅魯藏布江流域面積不是很大,為何本文的統(tǒng)計(jì)結(jié)果顯示其徑流的變化在一定程度上竟能對(duì)極移可能有影響?朱琳等[22]的研究表明,在從1948—2011年的64年間,全球山脈力矩變化最為顯著的地區(qū)集中在歐亞大陸的中南部和南美洲安第斯山脈,尤其以青藏高原區(qū)域變化最為劇烈。本文統(tǒng)計(jì)得到的雅魯藏布江徑流對(duì)極移的影響是因?yàn)槠湮恢锰幱诿舾袇^(qū)的特殊性?還是一種統(tǒng)計(jì)的必然(即抑或任何河川徑流都存在關(guān)系)?下一步我們將選擇世界大河進(jìn)行深入研究,進(jìn)一步揭示徑流對(duì)極移的影響機(jī)制。

    4.4 研究的引申意義

    陸地水儲(chǔ)量變化是極移激發(fā)的一個(gè)重要源,而極移反過來又可能通過離心力的變化來改變陸地水文循環(huán)過程。鑒于極移具有較長(zhǎng)時(shí)間序列,弄清楚極移與徑流之間的關(guān)系,可為借用極移資料提高資料稀少地區(qū)的水文預(yù)測(cè)精度提供寶貴的資料支持。資料稀缺區(qū)域的水文預(yù)報(bào)(Predictions in Ungauged Basins, PUB)是國(guó)際水文科學(xué)協(xié)會(huì)(IAHS)于2002年啟動(dòng)的一個(gè)十年水文計(jì)劃。通過類比,劉蘇峽等[23-24]總結(jié)出PUB研究的“借”、“替”、“種”方法三元論。循此方法論,若能厘清極移和徑流的關(guān)系,就可“借”用有較長(zhǎng)時(shí)間系列的極移資料的優(yōu)勢(shì),展延水文資料,加強(qiáng)寒區(qū)水文研究[25],提高水文預(yù)測(cè)的精度。

    圖7 雅魯藏布江水文循環(huán)過程示意圖

    4.5 研究的局限

    格蘭杰檢驗(yàn)的一個(gè)局限是它只能檢測(cè)出加入某個(gè)因子后是否與不加因子的預(yù)測(cè)效果是否一樣,但它尚不能展示,加入某個(gè)因子后,是否就會(huì)提高預(yù)測(cè)精度。很顯然,不一樣的后果也可能是降低預(yù)測(cè)精度。所以徑流和極移關(guān)系的探索任重而道遠(yuǎn)。

    在空間尺度上,極移運(yùn)動(dòng)和激發(fā)涉及到整個(gè)地球系統(tǒng),所以嚴(yán)格意義上,完全的研究應(yīng)選取全世界大型河流進(jìn)行研究。這將是下一步工作目標(biāo)。

    另外,我們只做了有限的滯后時(shí)間的比較,對(duì)于其他滯后時(shí)間,是否會(huì)帶來不同的結(jié)論,還需要更多探索。

    最后要強(qiáng)調(diào)的是,即使統(tǒng)計(jì)方法能做得盡善盡美,欲真正弄清極移與徑流的關(guān)系還需要走“統(tǒng)計(jì)——統(tǒng)計(jì)推斷——?jiǎng)恿C(jī)制”相結(jié)合的路徑。彭公炳等[3]在統(tǒng)計(jì)事實(shí)基礎(chǔ)上推斷,極移振幅增大時(shí),極移與氣候的聯(lián)系的可能性增強(qiáng)。極移振幅的高值年,亞歐中緯地區(qū)經(jīng)向環(huán)流指數(shù)增強(qiáng),緯向環(huán)流指數(shù)減弱,副高偏南,中緯度地區(qū)海洋向大陸輸送的水汽減少,因?yàn)榻邓疁p少,冬半年氣溫也降低,而地處副熱帶的我國(guó)長(zhǎng)江中下游降水有所增加,反之亦然。更深入的統(tǒng)計(jì)推斷和動(dòng)力機(jī)制將是下一步工作的方向。

    5 結(jié)論

    本文采用格蘭杰因果檢驗(yàn)方法,基于國(guó)際地球自轉(zhuǎn)和參考系服務(wù)地球定位參數(shù)計(jì)劃的極移數(shù)據(jù)和位于青藏高原的受人類活動(dòng)影響小的雅魯藏布江奴下站1978—2006年的月徑流資料,在月、季和年尺度上分別進(jìn)行了從極移到徑流(徑流變化預(yù)測(cè)中加極移變化)和從徑流到極移(極移變化預(yù)測(cè)中加徑流變化)的格蘭杰因果關(guān)系檢驗(yàn)。

    在月和季節(jié)尺度上,徑流和極移在0.05顯著性水平上顯著地互為格蘭杰因果。也就是說總能找到一個(gè)滯后時(shí)間并在該時(shí)間上建立自回歸(AR)和向量自回歸(VAR)模型,使得兩模型呈顯著差異,從而能提供改進(jìn)預(yù)測(cè)精度的契機(jī)。而在年尺度上,尚未檢測(cè)出因果關(guān)系。

    比較發(fā)生拒絕原假設(shè)的滯后模型的個(gè)數(shù),發(fā)現(xiàn)極移加入徑流對(duì)徑流預(yù)測(cè)精度的影響范圍要比徑流加入極移對(duì)極移預(yù)測(cè)精度的影響范圍大,也即滯后模型個(gè)數(shù)多,暗示極移對(duì)徑流的影響大于徑流對(duì)極移的影響。

    研究成果將可能用極移資料提高資料稀缺地區(qū)的

    水文預(yù)測(cè)精度。

    極移和徑流關(guān)系的深入的研究,將是對(duì)格蘭杰方法的改進(jìn)、采用全世界大河進(jìn)行全面研究、討論更多的滯后時(shí)間,特別是,采用“統(tǒng)計(jì)——統(tǒng)計(jì)推斷——?jiǎng)恿C(jī)制”相結(jié)合的方法進(jìn)行。

    [1] 趙豐. 見微知著——話地球自轉(zhuǎn). 科學(xué)月刊, 1990, 21(11): 866-876.

    [2]Chen J L, Wilson C R, Ries J C, et al. Rapid ice melting drives Earth’s pole to the east. Geophysical Research Letters, 2013, 40(11): 2625-2630.

    [3]彭公炳, 陸巍, 殷延珍. 地極移動(dòng)與氣候的幾個(gè)問題. 大氣科學(xué), 1980, 4(4): 369-378.

    [4]Wahr J M. Deformation induced by polar motion. Journal of Geophysical Research: Solid Earth, 1985, 90(B11): 9363-9368.

    [5]閆昊明, 鐘敏, 朱耀仲, 等. 極移半年振蕩的年際變化與北大西洋濤動(dòng). 自然科學(xué)進(jìn)展, 2002, 12(1): 104-107.

    [6]Chen J L, Wilson C R. Hydrological excitations of polar motion, 1993-2002. Geophysical Journal International, 2005, 160(3): 833-839.

    [7]Nastula J, Ponte R M, Salstein D A. Comparison of polar motion excitation series derived from GRACE and from analyses of geophysical fluids. Geophys Res Lett, 2007, 34, L11306, doi:10.1029/2006GL028983.

    [8]Jin S, Chambers D P, Tapley B D. Hydrological and oceanic e ff ects on polar motion from GRACE and models. J Geophys Res, 2010,115, B02.

    [9]Jin S G, Hassana A A, Feng G P. Assessment of terrestrial water contributions to polar motion from GRACE and hydrological models. Journal of Geodynamics, 2012, 62: 40-48.

    [10]Chen J L, Wilson C R, Zhou Y H. Seasonal excitation of polar motion. Journal of Geodynamics, 2012, 62: 8-15.

    [11]Schuh H, B?hm S. Earth rotation - In: Gupta H K (Ed), Encyclopedia of solid Earth geophysics (Encyclopedia of Earth Sciences Series 1). Springer, 2011.

    [12]Granger C W J. Investigating causal relations by econometric methods and cross-spectral methods. Econometrica, 1969, 37: 424-438.

    [13]Elsner J B. Granger causality and Atlantic hurricanes. Tellus, 2007, 59: 476-485.

    [14]He D, Ren J, Fu K, et al. Sediment change under climate changes and human activities in the Yuanjiang-Red River Basin. Chinese Science Bulletin, 2007, 52: 164-171.

    [15]Attanasio A. Testing for linear Granger causality from natural/ anthropogenic forcings to global temperature anomalies.Theoretical Appl Climatol, 2012, 110: 281-289.

    [16]Attanasio A, Pasini A, Triacca U. A contribution to attribution of recent global warming by out-of-sample Granger causality analysis. Atmos Sci Lett, 2012, 13: 67-72.

    [17]Stock J H, Watson M W. Interpreting the evidence on moneyincome causality. Journal of Econometrics, 1989, 40: 161-181.

    [18]Sims C A, Stock J H , Watson M W. Inference in linear time series models with some unit roots. Econometrica, 1990, 58: 113-144.

    [19]Dickey D, Fuller W A. Likelihood ratio statistic for autoregressive time series with a unit root. Econometrica, 1981, 49: 427-431.

    [20]劉天仇. 雅魯藏布江水文特征. 地理學(xué)報(bào), 1999, 54(s): 157-164.

    [21]Akaike H. A new look at the statistical model identi fi cation. IEEE Transactions on Automatic Control, 1974, 19(6): 716-723.

    [22]朱琳, 黃玫, 鞏賀, 等. 全球山脈力矩時(shí)空變化及其與地球自轉(zhuǎn)的關(guān)系. 氣象科技進(jìn)展, 2014, 4(3): 32-35.

    [23]劉蘇峽, 劉昌明, 趙衛(wèi)民. 無測(cè)站流域水文預(yù)測(cè) (PUB) 的研究方法. 地理科學(xué)進(jìn)展, 2010, 29(11): 1333-1339.

    [24]Liu S, Mo X, Liu C, et al, On Pub methodology from Chinese lessons. Hydrological Science Journal, 2014, doi: 10.1080/02626667.2014.899702.

    [25]Wang S, Liu S. Exploring the relationship between polar motion and a primitive river’s runoff based on Granger causality. IAHS Publ. 360, 2013.

    Exploring the Relationship between Polar Motion and Runoff

    Liu Suxia1, Wang Sheng1,2, Wang Yueling1, Li Mingxin3, Huang Mei1, Peng Gongbing1, Xiao Ziniu4
    (1 Key Laboratory of Water Cycle and Related Land Surface Processes , Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences , Beijing 100101 2 University of Chinese Academy of Sciences, Beijing 100049 3 Bureau of Hydrology, Changjiang Water Resources Commission, Wuhan 430010 4 China Meteorological Administration (CMA) Training Centre, Beijing 100081)

    The evidence of the relationship between polar motion and runoff is reviewed. The relationship between the runoff of Yarlung Zangbo, a primitive river in Tibet, China and polar motion with data records from about 1846 to present is explored. First, the action path framework was structured based on geophysical principles. In order to make time series stationary, the initial time series of polar motionXcomponent,Ycomponent and runoff were transformed into their corresponding fi rst order difference time series ΔX, ΔY, and ΔQ. The Granger causality test between ΔX, ΔYand ΔQwas conducted on monthly, seasonal and annual time scales. It is found that on a monthly scale ΔXinf l uences ΔQat the lag being from the 1stto the 21stmonth, with the 9thmonth being an exception. ΔYinf l uences ΔQat the lag from the 1stto the 9thmonth and from the 17thto the 24thmonth. On a seasonal scale, the inf l uence of ΔXon ΔQcan be seen in the 2ndseason (i.e. from the 4thto the 6thmonth). The inf l uence of ΔYon ΔQcan be seen on a seasonal scale from the 4thto the 6thseason (i.e. from the 10thto 18thmonth). We cannot see evident Granger causality from ΔX, ΔYto ΔQ, on annual scales. For the Granger inf l uence of ΔQto ΔX, ΔY, it is found that on a monthly scale the inf l uences are prominent at the lag being from the 3rdto the 25thmonths for the ΔX, for the ΔYit is the 1stand the 3rdto the 25thmonths. On a seasonal scale these inf l uences can be seen at the lag from the 2ndto the 8thseason (corresponding to the 4thto the 24thmonth) for ΔXand at the lag from the 1stto the 8thseason (corresponding to from the 1stto the 24thmonth) for ΔY. Again, on a annual scale no evident granger causality can be found from runoff to polar motion. The different behaviours on monthly, seasonal and annual scales suggest that using the monthly data of polar motion to obtain the monthly runoff data is more practicable than borrowing the data from polar motion for river runoff on seasonal and annual scales.

    polar motion, runoff, Granger causality, Yarlung Zangbo River

    10.3969/j.issn.2095-1973.2014.03.001

    2013年11月28日;

    2014年4月7日

    劉蘇峽(1965—),Email: liusx@igsnrr.ac.cn

    資助信息:國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2012CB957802);國(guó)家自然科學(xué)基金資助項(xiàng)目(41071024);長(zhǎng)江水利委員會(huì)水文局委托項(xiàng)目“長(zhǎng)江流域水中長(zhǎng)期供水規(guī)劃——典型流域冰川徑流計(jì)算及徑流組成分析”

    猜你喜歡
    格蘭杰因果關(guān)系徑流
    玩忽職守型瀆職罪中嚴(yán)重不負(fù)責(zé)任與重大損害后果的因果關(guān)系
    做完形填空題,需考慮的邏輯關(guān)系
    幫助犯因果關(guān)系芻議
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    格蘭杰因果關(guān)系在神經(jīng)科學(xué)領(lǐng)域的發(fā)展及缺陷
    電子科技(2015年8期)2015-12-18 13:17:56
    探秘“大徑流”
    攻克“大徑流”
    介入因素對(duì)因果關(guān)系認(rèn)定的影響
    榜單
    江埡水庫降雨徑流相關(guān)圖的建立
    午夜福利在线在线| 亚洲五月天丁香| 久久天躁狠狠躁夜夜2o2o| 欧美最黄视频在线播放免费| 日韩在线高清观看一区二区三区 | 国产精品,欧美在线| 国产黄色小视频在线观看| 国产欧美日韩一区二区精品| 大又大粗又爽又黄少妇毛片口| av天堂中文字幕网| 男人舔女人下体高潮全视频| 最近最新免费中文字幕在线| 久久精品国产鲁丝片午夜精品 | 嫁个100分男人电影在线观看| 女的被弄到高潮叫床怎么办 | 亚洲精华国产精华精| 国产v大片淫在线免费观看| 国产精品不卡视频一区二区| 人人妻人人澡欧美一区二区| 日韩欧美免费精品| 国产精品伦人一区二区| 欧美激情久久久久久爽电影| 久久这里只有精品中国| 成人无遮挡网站| 校园春色视频在线观看| 婷婷亚洲欧美| 国产蜜桃级精品一区二区三区| 日本免费a在线| 99热6这里只有精品| 日日摸夜夜添夜夜添小说| 国产极品精品免费视频能看的| 少妇丰满av| 天美传媒精品一区二区| 色av中文字幕| 白带黄色成豆腐渣| 亚洲综合色惰| 一本久久中文字幕| 成人美女网站在线观看视频| 久久久久久九九精品二区国产| 国产精品人妻久久久影院| www.www免费av| 日本在线视频免费播放| 国产精品久久久久久久电影| 天堂动漫精品| 国产精品美女特级片免费视频播放器| netflix在线观看网站| 男女啪啪激烈高潮av片| 免费大片18禁| 在线免费观看不下载黄p国产 | 在线天堂最新版资源| 欧美日韩国产亚洲二区| 91久久精品电影网| 亚洲性夜色夜夜综合| 很黄的视频免费| 中文字幕高清在线视频| 国产一区二区三区av在线 | 久久人人爽人人爽人人片va| 亚洲成人久久爱视频| 国产av麻豆久久久久久久| 在线天堂最新版资源| 亚洲一级一片aⅴ在线观看| a在线观看视频网站| 亚洲天堂国产精品一区在线| 不卡一级毛片| 国产一区二区在线av高清观看| 国产淫片久久久久久久久| 99热精品在线国产| 午夜影院日韩av| 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 婷婷精品国产亚洲av| 久久精品国产亚洲av香蕉五月| 久久午夜福利片| 亚洲第一电影网av| 日本-黄色视频高清免费观看| 美女高潮喷水抽搐中文字幕| 三级毛片av免费| 午夜视频国产福利| 亚洲真实伦在线观看| 黄色视频,在线免费观看| 日本黄色片子视频| 少妇熟女aⅴ在线视频| 精品人妻视频免费看| 精品免费久久久久久久清纯| 别揉我奶头 嗯啊视频| h日本视频在线播放| 婷婷六月久久综合丁香| 国产欧美日韩一区二区精品| 国产成人一区二区在线| 国产免费av片在线观看野外av| 国产精品久久视频播放| 欧美潮喷喷水| 国产男靠女视频免费网站| 俄罗斯特黄特色一大片| 国产v大片淫在线免费观看| 国产白丝娇喘喷水9色精品| 久久这里只有精品中国| 精品欧美国产一区二区三| 免费观看精品视频网站| 2021天堂中文幕一二区在线观| 国产精品女同一区二区软件 | 免费不卡的大黄色大毛片视频在线观看 | 男插女下体视频免费在线播放| 熟女电影av网| 亚洲精品国产成人久久av| 日本黄色片子视频| 久久香蕉精品热| 亚洲av成人av| 男人狂女人下面高潮的视频| h日本视频在线播放| 一进一出好大好爽视频| 能在线免费观看的黄片| 国产伦精品一区二区三区四那| 免费看光身美女| 男人舔女人下体高潮全视频| 日韩欧美免费精品| 成年女人毛片免费观看观看9| 白带黄色成豆腐渣| 精品久久国产蜜桃| 国产日本99.免费观看| 国产免费男女视频| 亚洲一级一片aⅴ在线观看| 麻豆成人午夜福利视频| 亚洲无线观看免费| 一a级毛片在线观看| 亚洲精华国产精华精| 国产成人a区在线观看| 日韩人妻高清精品专区| 一区二区三区四区激情视频 | 欧美成人免费av一区二区三区| 亚洲狠狠婷婷综合久久图片| 日本三级黄在线观看| 熟女电影av网| 露出奶头的视频| 免费看av在线观看网站| 熟女电影av网| 亚洲性夜色夜夜综合| 色av中文字幕| 午夜福利高清视频| 免费搜索国产男女视频| 可以在线观看毛片的网站| 亚洲精品在线观看二区| 俺也久久电影网| 成人欧美大片| 欧美bdsm另类| 1024手机看黄色片| 日本熟妇午夜| 色吧在线观看| 日本黄色片子视频| 亚洲一区二区三区色噜噜| 91久久精品国产一区二区三区| 久久久国产成人精品二区| 亚洲色图av天堂| 亚洲自偷自拍三级| 亚洲人成伊人成综合网2020| 国产高清不卡午夜福利| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| 国产女主播在线喷水免费视频网站 | 最近中文字幕高清免费大全6 | 美女大奶头视频| 欧美黑人巨大hd| 亚洲欧美日韩高清在线视频| 日韩精品青青久久久久久| 最近中文字幕高清免费大全6 | 精品国产三级普通话版| 日韩 亚洲 欧美在线| 久久久久免费精品人妻一区二区| 99久久精品热视频| 麻豆av噜噜一区二区三区| 精品福利观看| 国产成人福利小说| 欧洲精品卡2卡3卡4卡5卡区| av在线蜜桃| 亚洲av一区综合| 久久久久久久久久久丰满 | 色精品久久人妻99蜜桃| av在线亚洲专区| 国产高清激情床上av| 亚洲黑人精品在线| 免费看av在线观看网站| 国产精品美女特级片免费视频播放器| 日韩欧美 国产精品| 草草在线视频免费看| 日本爱情动作片www.在线观看 | 精品日产1卡2卡| 久久久久久久久久成人| 日韩欧美国产一区二区入口| 国国产精品蜜臀av免费| 一本一本综合久久| 最近最新中文字幕大全电影3| 国产精品1区2区在线观看.| 男人舔女人下体高潮全视频| а√天堂www在线а√下载| 啦啦啦啦在线视频资源| 国产在线男女| 久久久久久伊人网av| 亚洲成人久久爱视频| 亚洲欧美日韩东京热| 亚洲av美国av| 亚洲av不卡在线观看| 麻豆久久精品国产亚洲av| 听说在线观看完整版免费高清| 在线播放国产精品三级| 亚洲自拍偷在线| 一本久久中文字幕| 国产精品一区二区三区四区免费观看 | 亚洲精品粉嫩美女一区| 别揉我奶头 嗯啊视频| 91午夜精品亚洲一区二区三区 | 老熟妇仑乱视频hdxx| 日本a在线网址| 成人综合一区亚洲| 黄色女人牲交| 特大巨黑吊av在线直播| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 欧美日韩国产亚洲二区| 亚洲精品色激情综合| 欧美中文日本在线观看视频| 国内少妇人妻偷人精品xxx网站| 99国产精品一区二区蜜桃av| 美女被艹到高潮喷水动态| 欧美日韩精品成人综合77777| 久久久国产成人免费| 亚洲最大成人中文| 最近视频中文字幕2019在线8| 尤物成人国产欧美一区二区三区| 日本色播在线视频| 日本撒尿小便嘘嘘汇集6| 日本成人三级电影网站| 日韩av在线大香蕉| 国产精品亚洲美女久久久| 美女高潮喷水抽搐中文字幕| 国产精品永久免费网站| av中文乱码字幕在线| 性色avwww在线观看| 国产伦一二天堂av在线观看| 成人特级av手机在线观看| 日本在线视频免费播放| 亚洲精品成人久久久久久| 国产人妻一区二区三区在| 网址你懂的国产日韩在线| 丰满乱子伦码专区| 免费人成在线观看视频色| 国产主播在线观看一区二区| 午夜激情福利司机影院| 1000部很黄的大片| 午夜老司机福利剧场| 国产精品亚洲美女久久久| 久久人人爽人人爽人人片va| 久久国产精品人妻蜜桃| 精品久久久久久,| 一本一本综合久久| 亚洲成人免费电影在线观看| 免费人成视频x8x8入口观看| 欧美日韩国产亚洲二区| 琪琪午夜伦伦电影理论片6080| 国产在视频线在精品| 国产精品福利在线免费观看| 熟女人妻精品中文字幕| 不卡一级毛片| 春色校园在线视频观看| 一个人看的www免费观看视频| 噜噜噜噜噜久久久久久91| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 18禁黄网站禁片午夜丰满| 99视频精品全部免费 在线| 嫩草影院新地址| 日日摸夜夜添夜夜添小说| 国产私拍福利视频在线观看| 国产精品久久电影中文字幕| 深夜a级毛片| 午夜精品久久久久久毛片777| 国产黄片美女视频| 亚洲七黄色美女视频| av专区在线播放| 一进一出好大好爽视频| 国产91精品成人一区二区三区| 99久国产av精品| 国产乱人视频| 99在线视频只有这里精品首页| 国产乱人伦免费视频| 一级毛片久久久久久久久女| 韩国av一区二区三区四区| 国产亚洲91精品色在线| av黄色大香蕉| 国产蜜桃级精品一区二区三区| 国产精品伦人一区二区| 日韩欧美精品v在线| 国产人妻一区二区三区在| 中文亚洲av片在线观看爽| 日本a在线网址| netflix在线观看网站| 亚洲内射少妇av| av黄色大香蕉| 久久久久久伊人网av| 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 热99re8久久精品国产| 欧美三级亚洲精品| 中国美女看黄片| 少妇高潮的动态图| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 欧美zozozo另类| 亚洲国产日韩欧美精品在线观看| 淫秽高清视频在线观看| 成人精品一区二区免费| 亚洲国产欧美人成| 成人精品一区二区免费| 综合色av麻豆| 亚洲午夜理论影院| 嫩草影院精品99| 午夜精品一区二区三区免费看| 99热网站在线观看| 久久精品国产鲁丝片午夜精品 | 日韩大尺度精品在线看网址| 亚洲av成人av| 联通29元200g的流量卡| 哪里可以看免费的av片| 99热精品在线国产| 女生性感内裤真人,穿戴方法视频| 99热这里只有是精品50| 国产男人的电影天堂91| 国产色婷婷99| 欧美精品国产亚洲| 99精品在免费线老司机午夜| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 麻豆av噜噜一区二区三区| 国内精品久久久久久久电影| 国产免费男女视频| 我的女老师完整版在线观看| 最近在线观看免费完整版| 特级一级黄色大片| 我的女老师完整版在线观看| 欧美色欧美亚洲另类二区| 网址你懂的国产日韩在线| 亚洲精品粉嫩美女一区| 亚洲人成网站在线播| 22中文网久久字幕| 18禁在线播放成人免费| 91麻豆av在线| 亚洲电影在线观看av| 俄罗斯特黄特色一大片| 级片在线观看| 日韩人妻高清精品专区| 免费人成视频x8x8入口观看| 免费观看精品视频网站| 亚洲精品成人久久久久久| 午夜精品久久久久久毛片777| 伊人久久精品亚洲午夜| 丰满人妻一区二区三区视频av| 精品无人区乱码1区二区| 中国美白少妇内射xxxbb| 精品乱码久久久久久99久播| 国产亚洲欧美98| 欧美区成人在线视频| 最近视频中文字幕2019在线8| 亚洲图色成人| 国产男靠女视频免费网站| 18禁裸乳无遮挡免费网站照片| 91在线观看av| 亚洲av.av天堂| 国产视频内射| 国产熟女欧美一区二区| 国产精品精品国产色婷婷| 国产高清三级在线| 亚洲专区中文字幕在线| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 久久精品夜夜夜夜夜久久蜜豆| 真人一进一出gif抽搐免费| 亚洲国产精品合色在线| 欧美性猛交黑人性爽| 国内久久婷婷六月综合欲色啪| 91午夜精品亚洲一区二区三区 | 99热精品在线国产| 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 国产伦一二天堂av在线观看| 国产私拍福利视频在线观看| 日本黄色片子视频| 国产一区二区三区视频了| 两个人视频免费观看高清| 免费在线观看日本一区| 黄片wwwwww| 国产精品美女特级片免费视频播放器| 国产黄片美女视频| 狂野欧美激情性xxxx在线观看| 国产精品永久免费网站| 此物有八面人人有两片| 中国美女看黄片| 欧美高清成人免费视频www| 色综合亚洲欧美另类图片| 在线看三级毛片| 国产极品精品免费视频能看的| 免费看美女性在线毛片视频| 成人美女网站在线观看视频| 午夜福利欧美成人| 可以在线观看的亚洲视频| 免费av不卡在线播放| 舔av片在线| 男女那种视频在线观看| 婷婷色综合大香蕉| 99热6这里只有精品| 免费观看在线日韩| 女人十人毛片免费观看3o分钟| 精品久久久久久久人妻蜜臀av| 99热网站在线观看| 极品教师在线视频| 此物有八面人人有两片| 国产亚洲av嫩草精品影院| 久久久久国产精品人妻aⅴ院| 丰满乱子伦码专区| 日本一本二区三区精品| 国产一区二区在线观看日韩| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 天美传媒精品一区二区| 天堂av国产一区二区熟女人妻| 免费av毛片视频| 亚洲av不卡在线观看| 91午夜精品亚洲一区二区三区 | 国产精品不卡视频一区二区| 88av欧美| 日本三级黄在线观看| 少妇人妻一区二区三区视频| 国产人妻一区二区三区在| 九九在线视频观看精品| 午夜亚洲福利在线播放| 小蜜桃在线观看免费完整版高清| 高清日韩中文字幕在线| 亚洲avbb在线观看| 人人妻,人人澡人人爽秒播| 亚洲最大成人手机在线| 日韩精品中文字幕看吧| 嫩草影院精品99| 联通29元200g的流量卡| 高清日韩中文字幕在线| 亚洲avbb在线观看| 国产 一区精品| 韩国av在线不卡| 午夜福利高清视频| 好男人在线观看高清免费视频| av福利片在线观看| 精品午夜福利视频在线观看一区| av在线观看视频网站免费| 简卡轻食公司| 亚洲第一区二区三区不卡| 久久亚洲真实| 国产探花极品一区二区| 亚洲美女视频黄频| 免费在线观看成人毛片| 麻豆久久精品国产亚洲av| 在线观看午夜福利视频| 久久精品国产亚洲av涩爱 | 国产极品精品免费视频能看的| 欧美成人a在线观看| 人妻制服诱惑在线中文字幕| 中文字幕熟女人妻在线| 真人做人爱边吃奶动态| 欧美精品国产亚洲| 国产色爽女视频免费观看| 国产蜜桃级精品一区二区三区| 精品午夜福利视频在线观看一区| 免费看a级黄色片| 在线免费观看不下载黄p国产 | 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 免费观看的影片在线观看| av福利片在线观看| 毛片一级片免费看久久久久 | 国产免费av片在线观看野外av| 九九爱精品视频在线观看| 精品人妻视频免费看| 97超级碰碰碰精品色视频在线观看| 91午夜精品亚洲一区二区三区 | 日本免费一区二区三区高清不卡| 99riav亚洲国产免费| 少妇高潮的动态图| 成人午夜高清在线视频| 国产精品亚洲美女久久久| 91久久精品国产一区二区三区| 成人欧美大片| 国产亚洲91精品色在线| 久久久精品大字幕| 在线免费观看的www视频| 日本一本二区三区精品| 亚洲av.av天堂| 天堂影院成人在线观看| 久久精品国产亚洲av香蕉五月| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 亚洲成人中文字幕在线播放| 久久久国产成人精品二区| 午夜福利18| 性欧美人与动物交配| 少妇人妻一区二区三区视频| 久久精品国产99精品国产亚洲性色| 精品人妻1区二区| x7x7x7水蜜桃| 国产真实伦视频高清在线观看 | 淫秽高清视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 乱系列少妇在线播放| 少妇被粗大猛烈的视频| 老师上课跳d突然被开到最大视频| 丰满的人妻完整版| 国产一区二区三区在线臀色熟女| 欧美丝袜亚洲另类 | 99久久九九国产精品国产免费| 日日干狠狠操夜夜爽| 尾随美女入室| 少妇裸体淫交视频免费看高清| 国产午夜精品论理片| 日本免费a在线| 看免费成人av毛片| 亚洲aⅴ乱码一区二区在线播放| 亚洲最大成人中文| 一个人观看的视频www高清免费观看| 天天一区二区日本电影三级| 日本爱情动作片www.在线观看 | av.在线天堂| 国产亚洲精品久久久com| 精品一区二区三区人妻视频| av女优亚洲男人天堂| 男人舔奶头视频| 午夜激情欧美在线| 日本 av在线| 成人欧美大片| 久久亚洲真实| 日韩欧美国产一区二区入口| 国产av一区在线观看免费| 综合色av麻豆| 国产亚洲精品久久久久久毛片| 久久久久性生活片| 免费电影在线观看免费观看| 尾随美女入室| 给我免费播放毛片高清在线观看| 美女大奶头视频| 免费在线观看日本一区| 久久6这里有精品| 午夜福利高清视频| 亚洲男人的天堂狠狠| 日韩欧美三级三区| 老司机深夜福利视频在线观看| 男人和女人高潮做爰伦理| www日本黄色视频网| 亚洲人成网站高清观看| АⅤ资源中文在线天堂| 亚洲国产日韩欧美精品在线观看| 69人妻影院| 一本一本综合久久| 老司机深夜福利视频在线观看| 国产黄色小视频在线观看| 国产伦在线观看视频一区| а√天堂www在线а√下载| 国产精品1区2区在线观看.| 久久久久国内视频| 国内精品美女久久久久久| 好男人在线观看高清免费视频| 亚洲av中文av极速乱 | 成人av一区二区三区在线看| 国产麻豆成人av免费视频| 精品久久久久久,| 成熟少妇高潮喷水视频| 全区人妻精品视频| 给我免费播放毛片高清在线观看| 22中文网久久字幕| 午夜精品在线福利| 老熟妇仑乱视频hdxx| 国产精品一区www在线观看 | 男人的好看免费观看在线视频| 成人特级黄色片久久久久久久| 成人av在线播放网站| 亚洲av中文av极速乱 | 深夜a级毛片| 极品教师在线视频| 日本黄色视频三级网站网址| 97热精品久久久久久| 97碰自拍视频| 内地一区二区视频在线| 99热网站在线观看| 淫妇啪啪啪对白视频| 极品教师在线免费播放| 婷婷亚洲欧美| 18禁裸乳无遮挡免费网站照片| 色哟哟哟哟哟哟| 欧美3d第一页| 国产91精品成人一区二区三区| 亚洲av成人精品一区久久| 91久久精品国产一区二区三区| 久久精品国产亚洲av涩爱 | 999久久久精品免费观看国产| 国产真实乱freesex| 国产三级在线视频| 91麻豆av在线| 国产私拍福利视频在线观看| 免费一级毛片在线播放高清视频| 乱人视频在线观看| 亚洲av五月六月丁香网| 午夜精品在线福利| 免费看a级黄色片| 日本三级黄在线观看| 精品久久久久久久久av| 大又大粗又爽又黄少妇毛片口| 我要搜黄色片| 亚洲中文字幕日韩| 18禁裸乳无遮挡免费网站照片| 99久久久亚洲精品蜜臀av|