• <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)圖的建立
    久久人人爽人人片av| 亚洲欧美精品自产自拍| 亚洲欧美成人精品一区二区| 日韩一区二区视频免费看| 成人综合一区亚洲| 欧美变态另类bdsm刘玥| 国产伦理片在线播放av一区| 不卡视频在线观看欧美| 国产精品一区二区三区四区免费观看| 人妻制服诱惑在线中文字幕| kizo精华| 国产午夜精品一二区理论片| 亚洲av中文字字幕乱码综合| 国产伦精品一区二区三区四那| 久久久久久久久久成人| 一级毛片aaaaaa免费看小| 久久久久久久久中文| 国产亚洲精品av在线| 国产毛片a区久久久久| 免费av毛片视频| 夫妻性生交免费视频一级片| 人妻夜夜爽99麻豆av| 国产精品国产三级专区第一集| 国产真实伦视频高清在线观看| 久久国产乱子免费精品| 26uuu在线亚洲综合色| 免费不卡的大黄色大毛片视频在线观看 | 国产在视频线精品| 禁无遮挡网站| 老司机影院成人| 国产乱来视频区| 久久精品国产亚洲av天美| 高清在线视频一区二区三区 | 亚洲av中文av极速乱| 亚洲av成人精品一二三区| 又粗又硬又长又爽又黄的视频| 色哟哟·www| 日韩av在线大香蕉| 一本久久精品| 久久精品久久久久久久性| 青春草国产在线视频| 国产久久久一区二区三区| 又黄又爽又刺激的免费视频.| 国产探花在线观看一区二区| 中文资源天堂在线| 伊人久久精品亚洲午夜| 三级男女做爰猛烈吃奶摸视频| 久久久午夜欧美精品| 国产 一区精品| 97超视频在线观看视频| 亚洲av电影在线观看一区二区三区 | 天堂影院成人在线观看| 一级爰片在线观看| 日韩欧美在线乱码| 中文资源天堂在线| 国产高清不卡午夜福利| 国产三级中文精品| 国产精品久久久久久精品电影| 中国国产av一级| 精品久久久久久久末码| 亚洲成人久久爱视频| 国产精品久久电影中文字幕| 亚洲三级黄色毛片| 少妇的逼水好多| 久久久久久久久久久免费av| 男人舔奶头视频| 色噜噜av男人的天堂激情| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品自拍成人| 男人舔女人下体高潮全视频| 亚洲精品色激情综合| 在线观看美女被高潮喷水网站| 国产成人a∨麻豆精品| 国产老妇伦熟女老妇高清| 欧美成人免费av一区二区三区| 丝袜喷水一区| 伦精品一区二区三区| 九九在线视频观看精品| 国产三级中文精品| av国产免费在线观看| 久久久久久久久久黄片| 18禁在线播放成人免费| 国产免费视频播放在线视频 | 久热久热在线精品观看| 日韩欧美国产在线观看| 久久亚洲国产成人精品v| 精品国产一区二区三区久久久樱花 | 欧美潮喷喷水| 国产免费又黄又爽又色| 免费人成在线观看视频色| kizo精华| 黄色日韩在线| 国产高潮美女av| 精品久久久久久久久亚洲| 国产黄a三级三级三级人| 日本av手机在线免费观看| 在现免费观看毛片| 热99re8久久精品国产| 国产综合懂色| 91精品一卡2卡3卡4卡| 国产精品一区二区在线观看99 | 久久精品国产亚洲av涩爱| 国产黄片美女视频| 日产精品乱码卡一卡2卡三| 美女国产视频在线观看| 一级二级三级毛片免费看| 97热精品久久久久久| 精品无人区乱码1区二区| 欧美成人免费av一区二区三区| 99久久精品一区二区三区| 色综合色国产| 国产高清不卡午夜福利| 精品熟女少妇av免费看| 国产三级中文精品| 亚洲综合色惰| 青春草亚洲视频在线观看| 一级毛片我不卡| 久久久国产成人精品二区| 国产欧美日韩精品一区二区| 免费看日本二区| 在线播放无遮挡| 国产极品天堂在线| 综合色av麻豆| 国产中年淑女户外野战色| 国产精品电影一区二区三区| 精品一区二区三区人妻视频| 亚洲自偷自拍三级| 亚洲av免费高清在线观看| 国产大屁股一区二区在线视频| 亚洲经典国产精华液单| 精品国产三级普通话版| 丝袜美腿在线中文| 99热这里只有是精品在线观看| 国产成人a区在线观看| 看免费成人av毛片| 婷婷色麻豆天堂久久 | 青春草视频在线免费观看| 国产精品一区二区在线观看99 | 三级男女做爰猛烈吃奶摸视频| 神马国产精品三级电影在线观看| 我的女老师完整版在线观看| 亚洲av一区综合| 大香蕉97超碰在线| 老司机福利观看| 联通29元200g的流量卡| 国产毛片a区久久久久| 日本免费一区二区三区高清不卡| 久久久久久久久大av| 免费大片18禁| 男插女下体视频免费在线播放| av线在线观看网站| 搞女人的毛片| 国产真实伦视频高清在线观看| 97在线视频观看| 一级毛片我不卡| 老司机福利观看| 精品少妇黑人巨大在线播放 | 免费看av在线观看网站| 午夜免费男女啪啪视频观看| 国产熟女欧美一区二区| www日本黄色视频网| 69人妻影院| 超碰97精品在线观看| 日韩在线高清观看一区二区三区| 国产午夜福利久久久久久| 亚洲综合色惰| 亚洲精品影视一区二区三区av| 免费看av在线观看网站| 波多野结衣巨乳人妻| 天堂中文最新版在线下载 | 久久人人爽人人片av| 成人综合一区亚洲| 99久久九九国产精品国产免费| 亚洲欧美日韩无卡精品| 久久久久久大精品| 精品欧美国产一区二区三| 亚洲av不卡在线观看| 久久久久久久久中文| 在线观看av片永久免费下载| 亚洲国产欧洲综合997久久,| 午夜免费激情av| 又爽又黄a免费视频| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久久久精品电影| 免费搜索国产男女视频| 久久精品国产99精品国产亚洲性色| 免费一级毛片在线播放高清视频| 你懂的网址亚洲精品在线观看 | 一级爰片在线观看| 欧美+日韩+精品| 91久久精品电影网| 精品一区二区免费观看| 99国产精品一区二区蜜桃av| 久久久久性生活片| 一级毛片aaaaaa免费看小| 国产女主播在线喷水免费视频网站 | 熟女电影av网| 亚洲成人精品中文字幕电影| 综合色丁香网| 久久精品国产亚洲网站| 日本黄大片高清| 女人十人毛片免费观看3o分钟| 少妇熟女欧美另类| 一级爰片在线观看| 国产精品1区2区在线观看.| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 欧美人与善性xxx| 日本色播在线视频| 嫩草影院精品99| 中文字幕精品亚洲无线码一区| 啦啦啦韩国在线观看视频| 边亲边吃奶的免费视频| 少妇的逼好多水| 99久国产av精品国产电影| 国产乱人偷精品视频| 色网站视频免费| 久久99热这里只有精品18| 97热精品久久久久久| videos熟女内射| 国产精品一区二区在线观看99 | 中文字幕av在线有码专区| 特大巨黑吊av在线直播| 波多野结衣巨乳人妻| 亚洲精品日韩在线中文字幕| 色哟哟·www| 能在线免费观看的黄片| 国产色婷婷99| 日本黄大片高清| 国产一区二区在线av高清观看| 水蜜桃什么品种好| 欧美高清成人免费视频www| 亚洲欧美日韩无卡精品| av视频在线观看入口| 中文字幕亚洲精品专区| 久久久久久久久久成人| 精华霜和精华液先用哪个| 视频中文字幕在线观看| 久久久国产成人精品二区| 国产成人freesex在线| 久久这里只有精品中国| 午夜福利视频1000在线观看| 成人一区二区视频在线观看| 久久精品人妻少妇| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 久久亚洲国产成人精品v| 深夜a级毛片| 久久精品久久精品一区二区三区| 亚洲av电影不卡..在线观看| 伦理电影大哥的女人| 五月伊人婷婷丁香| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播| 在线a可以看的网站| 国产成年人精品一区二区| 成人性生交大片免费视频hd| 大又大粗又爽又黄少妇毛片口| 日韩人妻高清精品专区| 亚洲av电影在线观看一区二区三区 | 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 美女xxoo啪啪120秒动态图| 嘟嘟电影网在线观看| 一个人看的www免费观看视频| 免费电影在线观看免费观看| 最近手机中文字幕大全| av在线老鸭窝| 国产极品精品免费视频能看的| 亚洲18禁久久av| 亚洲最大成人手机在线| 国产成人精品一,二区| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 日本黄色片子视频| 一个人看的www免费观看视频| 国产亚洲91精品色在线| 成人性生交大片免费视频hd| 97超碰精品成人国产| 少妇丰满av| 丰满少妇做爰视频| 亚洲人成网站高清观看| 日本免费一区二区三区高清不卡| 视频中文字幕在线观看| 亚洲图色成人| 久久热精品热| 女人十人毛片免费观看3o分钟| 日本一本二区三区精品| 汤姆久久久久久久影院中文字幕 | 中国国产av一级| 内地一区二区视频在线| 在线天堂最新版资源| 国产精品日韩av在线免费观看| 91精品国产九色| 看十八女毛片水多多多| 爱豆传媒免费全集在线观看| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验| 日本五十路高清| 在线观看美女被高潮喷水网站| 久久精品久久久久久久性| 一个人看的www免费观看视频| 长腿黑丝高跟| 人妻系列 视频| 国产探花极品一区二区| 男人舔女人下体高潮全视频| 午夜老司机福利剧场| 欧美97在线视频| 亚洲国产精品sss在线观看| 日韩高清综合在线| 亚洲人成网站高清观看| 国产精品一区二区性色av| 乱人视频在线观看| 欧美成人免费av一区二区三区| 精品国产一区二区三区久久久樱花 | 三级男女做爰猛烈吃奶摸视频| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久久久免| 午夜福利在线观看免费完整高清在| 欧美性猛交黑人性爽| 如何舔出高潮| 你懂的网址亚洲精品在线观看 | 看黄色毛片网站| 久久国产乱子免费精品| 国产精华一区二区三区| 秋霞伦理黄片| 人人妻人人澡欧美一区二区| av卡一久久| 波多野结衣巨乳人妻| 九九久久精品国产亚洲av麻豆| 精品欧美国产一区二区三| 国产探花极品一区二区| 免费不卡的大黄色大毛片视频在线观看 | 一级毛片久久久久久久久女| 亚洲图色成人| 欧美精品一区二区大全| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 亚洲欧美精品自产自拍| 亚洲av一区综合| 中文字幕免费在线视频6| 91久久精品电影网| 日韩av在线大香蕉| 成人欧美大片| 97人妻精品一区二区三区麻豆| 亚洲最大成人手机在线| 内射极品少妇av片p| 爱豆传媒免费全集在线观看| 午夜激情福利司机影院| 大香蕉久久网| 一边亲一边摸免费视频| 蜜桃亚洲精品一区二区三区| 国产高潮美女av| 91久久精品国产一区二区成人| 99热精品在线国产| 日韩中字成人| 亚洲人与动物交配视频| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 黄片wwwwww| 一本一本综合久久| 黄色一级大片看看| 九色成人免费人妻av| 国产一级毛片在线| av在线蜜桃| 搡老妇女老女人老熟妇| 国产高清三级在线| 2021天堂中文幕一二区在线观| av国产久精品久网站免费入址| 看片在线看免费视频| 国产精品av视频在线免费观看| 久久人人爽人人片av| 老司机影院毛片| 舔av片在线| 国产成人精品一,二区| 看非洲黑人一级黄片| 一级黄片播放器| 国产午夜福利久久久久久| 久久99热这里只频精品6学生 | 日本免费一区二区三区高清不卡| 国产精品久久久久久av不卡| kizo精华| 国产高清视频在线观看网站| 极品教师在线视频| 我的老师免费观看完整版| 一区二区三区四区激情视频| 日韩一区二区视频免费看| 国产成人91sexporn| 国产真实乱freesex| 国产精品1区2区在线观看.| av在线天堂中文字幕| 日韩视频在线欧美| 美女大奶头视频| 51国产日韩欧美| 全区人妻精品视频| 狠狠狠狠99中文字幕| 99久久精品一区二区三区| 搡女人真爽免费视频火全软件| 久久久久性生活片| 日本-黄色视频高清免费观看| 噜噜噜噜噜久久久久久91| av线在线观看网站| 亚洲国产精品sss在线观看| 久久久国产成人免费| 伦理电影大哥的女人| 欧美性猛交黑人性爽| 日韩欧美国产在线观看| 男人舔女人下体高潮全视频| 国产真实乱freesex| 免费观看在线日韩| 欧美成人a在线观看| 欧美日韩国产亚洲二区| 国产亚洲91精品色在线| 嘟嘟电影网在线观看| 亚洲精品自拍成人| 99久国产av精品国产电影| www.av在线官网国产| 色哟哟·www| 观看免费一级毛片| 看非洲黑人一级黄片| 99热全是精品| 两个人的视频大全免费| 能在线免费看毛片的网站| 国产成人freesex在线| 国产精品一区二区三区四区久久| 高清毛片免费看| 亚洲欧美精品专区久久| 国产真实乱freesex| .国产精品久久| 最近中文字幕高清免费大全6| 尾随美女入室| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕| 26uuu在线亚洲综合色| 美女高潮的动态| 国产精品综合久久久久久久免费| 国模一区二区三区四区视频| 欧美+日韩+精品| 长腿黑丝高跟| 高清在线视频一区二区三区 | 日韩成人av中文字幕在线观看| 少妇的逼水好多| 亚洲av成人av| 国产精品麻豆人妻色哟哟久久 | 国产精品一区二区三区四区久久| 别揉我奶头 嗯啊视频| 美女大奶头视频| 久久99热6这里只有精品| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区 | 免费av观看视频| 国产亚洲午夜精品一区二区久久 | 久久99精品国语久久久| 99久久无色码亚洲精品果冻| 欧美日韩国产亚洲二区| 色综合站精品国产| 黄色欧美视频在线观看| 亚洲精品乱久久久久久| 亚洲精品一区蜜桃| 热99在线观看视频| 精品99又大又爽又粗少妇毛片| 亚洲国产精品sss在线观看| 少妇熟女欧美另类| 日本免费一区二区三区高清不卡| 在现免费观看毛片| 天堂中文最新版在线下载 | h日本视频在线播放| 我要搜黄色片| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 看片在线看免费视频| 久久精品影院6| 国产一区有黄有色的免费视频 | 精品人妻视频免费看| 少妇熟女欧美另类| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 一区二区三区高清视频在线| 亚洲欧美日韩高清专用| 成年女人看的毛片在线观看| 午夜福利在线观看吧| 女人久久www免费人成看片 | 岛国毛片在线播放| 校园人妻丝袜中文字幕| 特级一级黄色大片| 欧美一区二区亚洲| 国产伦一二天堂av在线观看| 成人漫画全彩无遮挡| 亚洲国产日韩欧美精品在线观看| 亚洲精品国产成人久久av| 亚洲av中文字字幕乱码综合| 麻豆国产97在线/欧美| 国产色婷婷99| 亚洲av成人av| 国产探花在线观看一区二区| 欧美精品国产亚洲| av在线老鸭窝| 亚洲av福利一区| 联通29元200g的流量卡| 一边亲一边摸免费视频| 国产淫语在线视频| 九色成人免费人妻av| 联通29元200g的流量卡| 床上黄色一级片| 日韩一区二区视频免费看| 亚洲精品国产成人久久av| 亚洲人成网站在线观看播放| or卡值多少钱| 国产亚洲精品av在线| 观看美女的网站| 日韩国内少妇激情av| 久久99热这里只频精品6学生 | 亚洲av.av天堂| 久久久a久久爽久久v久久| 亚洲av.av天堂| 最近最新中文字幕大全电影3| 亚洲一级一片aⅴ在线观看| 一级爰片在线观看| АⅤ资源中文在线天堂| 波多野结衣巨乳人妻| 一级二级三级毛片免费看| 国产成人freesex在线| 国产成人a∨麻豆精品| 美女脱内裤让男人舔精品视频| 性色avwww在线观看| 亚洲欧美精品自产自拍| 久久久国产成人免费| 天天躁夜夜躁狠狠久久av| 韩国av在线不卡| 亚洲精品久久久久久婷婷小说 | 内射极品少妇av片p| 校园人妻丝袜中文字幕| 精品酒店卫生间| 国产精品久久久久久av不卡| 免费av不卡在线播放| 国产亚洲精品久久久com| 国产女主播在线喷水免费视频网站 | 一级二级三级毛片免费看| 国产免费男女视频| 久久韩国三级中文字幕| 国产熟女欧美一区二区| 国产精品久久久久久精品电影小说 | 深夜a级毛片| 欧美成人午夜免费资源| 我的女老师完整版在线观看| 青春草视频在线免费观看| 久久久久久国产a免费观看| 尾随美女入室| 91在线精品国自产拍蜜月| 亚洲色图av天堂| 精品久久久久久久末码| 中文亚洲av片在线观看爽| 日韩亚洲欧美综合| ponron亚洲| 人妻系列 视频| www日本黄色视频网| 亚洲美女搞黄在线观看| 美女国产视频在线观看| 十八禁国产超污无遮挡网站| АⅤ资源中文在线天堂| 少妇的逼水好多| a级毛片免费高清观看在线播放| av在线亚洲专区| 99国产精品一区二区蜜桃av| 三级男女做爰猛烈吃奶摸视频| 亚洲最大成人手机在线| 亚洲精品自拍成人| 一级av片app| 久久99热这里只有精品18| 成人美女网站在线观看视频| 成年女人永久免费观看视频| 亚洲成av人片在线播放无| 精品免费久久久久久久清纯| 一个人看的www免费观看视频| 日本一二三区视频观看| 日韩av不卡免费在线播放| 男人狂女人下面高潮的视频| 搡女人真爽免费视频火全软件| 日本一本二区三区精品| 91久久精品国产一区二区三区| 精品少妇黑人巨大在线播放 | 国产高清不卡午夜福利| 国内精品宾馆在线| 蜜臀久久99精品久久宅男| 亚洲久久久久久中文字幕| 亚洲精品色激情综合| 免费av观看视频| 色播亚洲综合网| 国产精品一区二区性色av| 桃色一区二区三区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产免费又黄又爽又色| 欧美一区二区亚洲| 三级经典国产精品| 亚洲av免费高清在线观看| 日韩av在线免费看完整版不卡| 如何舔出高潮| 草草在线视频免费看| 国产91av在线免费观看| 免费观看a级毛片全部| 国产又黄又爽又无遮挡在线| 国产一区亚洲一区在线观看| 久久久久久久亚洲中文字幕| 国产精品国产三级国产av玫瑰| 久久久久久久久中文| 日韩av在线大香蕉| 亚洲第一区二区三区不卡| 乱系列少妇在线播放| 亚洲精品,欧美精品| 国产av一区在线观看免费| 国产亚洲午夜精品一区二区久久 | 韩国av在线不卡| 搡女人真爽免费视频火全软件| 中文字幕久久专区|