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

    大興安嶺北部阿里河樟子松年輪氣候響應(yīng)及冬季降水重建

    2014-06-27 05:50:38呂姍娜王曉春
    關(guān)鍵詞:樹(shù)輪年表樟子松

    呂姍娜,王曉春

    (東北林業(yè)大學(xué)生態(tài)研究中心,黑龍江哈爾濱150040)

    大興安嶺北部阿里河樟子松年輪氣候響應(yīng)及冬季降水重建

    呂姍娜,王曉春

    (東北林業(yè)大學(xué)生態(tài)研究中心,黑龍江哈爾濱150040)

    對(duì)大興安嶺阿里河地區(qū)樟子松(Pinus sylvestris var.mongolica)樹(shù)輪年表與氣候要素的響應(yīng)關(guān)系進(jìn)行了分析研究.結(jié)果表明:年輪指數(shù)與該地區(qū)前一年12月至當(dāng)年1月平均降水存在顯著的正相關(guān)關(guān)系.基于樹(shù)輪寬度重建了阿里河1809年以來(lái)前一年12月至當(dāng)年1月平均降水量回歸方程,該方程解釋了47%的冬季降水變化,且該序列與黑龍江省干濕變化存在較好的一致性.對(duì)重建序列的分析表明,該地區(qū)冬季降水變化存在4個(gè)偏干階段和4個(gè)偏濕階段.阿里河地區(qū)過(guò)去200a冬季降水存在干濕交替時(shí)期,即19世紀(jì)初期較干旱,中期較濕潤(rùn),后期至20世紀(jì)初期較干旱,而20世紀(jì)中期較濕潤(rùn),后期則相對(duì)偏干.同時(shí),阿里河地區(qū)冬季降水存在85~38,5.0~4.5和2.2a顯著的變化周期.

    阿里河;樟子松;冬季降水;樹(shù)木年輪

    樹(shù)木作為生態(tài)系統(tǒng)的主要組成部分,在生長(zhǎng)過(guò)程中每一年的年輪形成都受當(dāng)年及生長(zhǎng)前期的許多氣候因子的影響,這種影響在樹(shù)木生長(zhǎng)和年輪結(jié)構(gòu)中尤為重要[1].而在干旱、半干旱地區(qū),降水對(duì)年輪寬窄變化的影響比較大[2].大興安嶺是我國(guó)東北最大的原始林區(qū),該地區(qū)降水十分豐富,會(huì)對(duì)河川徑流、農(nóng)牧業(yè)生產(chǎn)產(chǎn)生影響,甚至有時(shí)會(huì)導(dǎo)致洪澇災(zāi)害的發(fā)生.阿里河位于大興安嶺北部地區(qū),了解該地區(qū)的多年降雪情況對(duì)大興安嶺林區(qū)的森林培育管理以及林區(qū)人們的生產(chǎn)生活十分重要.

    我國(guó)北方的大部分地區(qū),氣象站的降雪記錄幾乎均不超過(guò)50年,對(duì)理解和探索多年降雪變化特征顯得太短.而樹(shù)木年輪以其準(zhǔn)確性、連續(xù)性強(qiáng)、分辨率高、地域分布范圍廣泛和復(fù)本量好等特點(diǎn),已成為在全球氣候變化研究中獲取過(guò)去氣候環(huán)境演變數(shù)據(jù)的重要方法之一[3-5].目前,已有許多學(xué)者在不同地區(qū)利用樹(shù)木年輪數(shù)據(jù)開(kāi)展了降水重建的研究[6-8],但是有關(guān)冬季降水(降雪)的重建研究相對(duì)較少[9-12],而在冬季降雪非常豐富的大興安嶺地區(qū)尚未開(kāi)展此項(xiàng)研究.本文利用大興安嶺北部地區(qū)的主要樹(shù)種樟子松的樹(shù)木年輪資料,將降雪量轉(zhuǎn)換為冬季降水量(采用國(guó)家標(biāo)準(zhǔn)轉(zhuǎn)換方法),重建了1742年以來(lái)前一年12月至當(dāng)年1月的平均降水序列,進(jìn)而討論了冬季降水變化的特征,以為大興安嶺北部地區(qū)的氣候重建提供一些數(shù)據(jù)資料,同時(shí)也對(duì)完善大興安嶺地區(qū)樹(shù)輪資料庫(kù)提供一些理論依據(jù).

    1 材料與方法

    1.1 研究區(qū)域概況

    本文的研究地點(diǎn)為阿里河林業(yè)局的興阿林場(chǎng)(50°38′37.3″N、124°28′28.7″E,海拔376m),位于大興安嶺的北部地區(qū)(見(jiàn)圖1).該地區(qū)屬寒溫帶大陸性季風(fēng)氣候,冬寒夏涼,晝夜溫差較大.年平均氣溫為-3.5℃,極端最低氣溫達(dá)到-52.3℃;多年平均降水量400~600mm,主要集中在6—8月(見(jiàn)圖2),約占全年降水量的65%,全年無(wú)霜期90~110d[13].

    圖1 大興安嶺阿里河樟子松樹(shù)輪采樣點(diǎn)和氣象站位置圖

    本研究所用樹(shù)種為樟子松(Pinus sylvestris var.mongolica),它耐寒抗旱、適應(yīng)性強(qiáng),是我國(guó)松屬中最耐寒的樹(shù)種之一,主要分布在大興安嶺北部海拔450~980m的地區(qū).樟子松多為純林,間或混有少量的興安落葉松.林下灌木主要包括興安杜鵑(Rhododendron dauricum)、山刺玫(Rosa davurica)、柳葉繡線菊(Spiraea salicifolia)、篤斯越桔(Vaccinium uliginosum)和小葉杜香(Ledum palustre)等;草本植物發(fā)育良好,多由耐旱種組成,如興安野青茅(Calamagrostic turczaninowii)、矮山黧豆(Lathyrus humilis)、廣布野豌豆(Vicia cracca)等.樟子松林一般面積不大,分散鑲嵌在興安落葉松林間[13].

    1.2 樣品的采集與處理

    野外取樣的時(shí)間為2010年7月份;取樣樹(shù)種為樟子松,所取樹(shù)木一般為胸徑36~58cm的健康活立木.基于國(guó)際樹(shù)木年輪庫(kù)的標(biāo)準(zhǔn),用內(nèi)徑為5.3mm的生長(zhǎng)錐采取了32棵樹(shù)木的年輪樣芯,每株樹(shù)木采集1個(gè),共采集32根樣芯.

    將采集的樣芯帶回實(shí)驗(yàn)室,根據(jù)樹(shù)木年輪樣本處理的基本程序[14],對(duì)采集的樣芯進(jìn)行預(yù)處理(粘貼、固定、打磨).在雙筒顯微鏡下用骨架示意圖法進(jìn)行目視交叉定年,然后使用LINTAB寬度測(cè)量?jī)x對(duì)樹(shù)木年輪寬度進(jìn)行逐年測(cè)量,該測(cè)量系統(tǒng)精度可達(dá)0.001mm.最后利用COFECHA程序?qū)Χ旰蜏y(cè)量結(jié)果進(jìn)行檢驗(yàn)[15],32根樹(shù)芯成功進(jìn)行了交叉定年,除3個(gè)樣芯(ALH15、ALH19、ALH22)局部片段上稍有問(wèn)題外,其他樣芯的COFECHA檢驗(yàn)質(zhì)量均較高,樣芯中沒(méi)有出現(xiàn)缺輪現(xiàn)象.

    1.3 年表建立

    為減少因年齡等因素帶來(lái)的影響,對(duì)經(jīng)過(guò)交叉定年的年輪序列,利用ARSTAN程序進(jìn)行去趨勢(shì)和標(biāo)準(zhǔn)化處理[16].去趨勢(shì)方法采用負(fù)指數(shù)或線性函數(shù)進(jìn)行,以便保留更多的低頻變化信號(hào);采用雙權(quán)重平均法進(jìn)行年輪曲線的標(biāo)準(zhǔn)化,最后得到標(biāo)準(zhǔn)年表、差值年表和自回歸年表(見(jiàn)圖3).綜合分析年表特征,最后選擇標(biāo)準(zhǔn)年表用于隨后的分析.

    圖2 大興安嶺小二溝氣象站月平均氣溫和月總降水量(1957—2009年)

    1.4 氣候資料

    本文采用的氣候數(shù)據(jù)取自離采樣地較近的小二溝氣象站(49°12′N(xiāo)、123°43′E,海拔286m,記錄時(shí)段1957—2009年),該地區(qū)年平均降水量為496.42mm,10—3月平均降水占全年的8.5%,全年平均溫度為-3.61℃.將月總降水量(此處的降水為液態(tài)和固態(tài)的降水量,降雪到降水的轉(zhuǎn)換采用國(guó)家標(biāo)準(zhǔn)進(jìn)行)、月平均溫度與樹(shù)輪寬度指數(shù)進(jìn)行相關(guān)分析.

    1.5 統(tǒng)計(jì)分析

    為分析樹(shù)木徑向生長(zhǎng)的影響因素,將每年各個(gè)單月及不同月份組合的平均降水量與樹(shù)輪寬度指數(shù)序列分別進(jìn)行相關(guān)分析.經(jīng)分析發(fā)現(xiàn)前一年12月至當(dāng)年1月平均降水量與年表相關(guān)性最好,因此用多元線性回歸來(lái)獲得重建方程.重建方程的可靠性用誤差縮減值(reduction of error)、符號(hào)檢驗(yàn)(sign test)和乘積平均值檢驗(yàn)(product means test)進(jìn)行檢驗(yàn),這幾個(gè)統(tǒng)計(jì)值的具體檢驗(yàn)方法參照Cook等文獻(xiàn)[17].用多帶譜分析(multi-taper methods analysis)檢查重建降水的周期性.文中所有的統(tǒng)計(jì)分析用SPSS 16.0完成.

    2 結(jié)果與討論

    2.1 年表統(tǒng)計(jì)特征

    從樟子松3個(gè)年表的主要特征參數(shù)及共同區(qū)間分析結(jié)果(見(jiàn)表1)來(lái)看,年輪指數(shù)的一階自相關(guān)較高,說(shuō)明樹(shù)木前一年的生長(zhǎng)狀況對(duì)當(dāng)年生長(zhǎng)有一定的影響.年表的平均敏感度、信噪比和第一特征解釋量都比較高,表明樟子松年輪中含有的高頻變化信號(hào)較為理想.并且標(biāo)準(zhǔn)年表(STD)和差值年表(RES)的樣本總體代表性也達(dá)到了0.87和0.89,說(shuō)明采集的樣本量所含有的信號(hào)是足以代表總體特征的.以上各項(xiàng)年表統(tǒng)計(jì)特征值表明,樣本之間的樹(shù)輪寬度變化有很好的一致性,可靠記錄了當(dāng)?shù)氐臍夂蛐畔?,所以可以用?lái)進(jìn)行年輪氣候?qū)W的分析.

    表1 阿里河樟子松標(biāo)準(zhǔn)、差值和自回歸年表的主要特征參數(shù)

    圖3 阿里河樟子松標(biāo)準(zhǔn)、殘差和自回歸年表及樣本量

    2.2 樟子松年輪-氣候關(guān)系

    通過(guò)對(duì)樹(shù)木生長(zhǎng)與氣候因子的相關(guān)分析和響應(yīng)分析,確定所要重建的氣候要素.利用SPSS軟件對(duì)樹(shù)輪寬度指數(shù)序列與1957—2009年的每個(gè)月平均溫度和月降水量以及各個(gè)季節(jié)的平均溫度和降水量進(jìn)行相關(guān)分析,結(jié)果見(jiàn)圖4,最終選擇最佳的重建時(shí)期.

    圖4 標(biāo)準(zhǔn)年表(STD)樹(shù)輪寬度指數(shù)序列(A)和差值年表(RES)樹(shù)輪寬度指數(shù)序列(B)與氣候要素的相關(guān)分析

    從表1和圖4(A)可以看出,標(biāo)準(zhǔn)年表(STD)包含了當(dāng)?shù)卦S多氣候信息,降水量除了前一年11月、當(dāng)年2月和5月外,年輪指數(shù)與其他月份的降水量都呈正相關(guān),而且與前一年12月至當(dāng)年1月平均降水量的相關(guān)性達(dá)到了顯著水平,適合進(jìn)行氣候重建.前一年冬季降雪增加,有利于來(lái)年春季土壤水分的存留,從而有利于第二年樹(shù)木的生長(zhǎng).

    從表1和圖4(B)可以看出,雖然差值年表(RES)也包含了許多氣候信息,并且年輪指數(shù)與當(dāng)年1月、6月和前一年至當(dāng)年1月平均降水量以及前一年11月的溫度通過(guò)了95%的置信水平,但是達(dá)不到重建的條件,所以不適合用差值年表(RES)進(jìn)行氣候的重建.

    溫度響應(yīng)和降水有一定的差異,除6月份以外,其他月份的溫度與樹(shù)輪寬度指數(shù)呈正相關(guān),并且2,5,8,9月以及前一年10月和11月的溫度也通過(guò)了95%的置信區(qū)間.隨著6,7月溫度的逐步升高,雨季來(lái)臨,水分條件可滿足樹(shù)木生長(zhǎng)的需要,此時(shí),溫度、降水對(duì)樹(shù)木生長(zhǎng)過(guò)程的限制比前期要弱得多,年輪寬度與氣候的相關(guān)性偏低[18].

    2.3 冬季降水重建

    通過(guò)樟子松年輪-氣候關(guān)系的分析,以及對(duì)每年各個(gè)單月及不同月份組合的平均降水量與樹(shù)輪寬度指數(shù)序列分別進(jìn)行相關(guān)分析,最終選取前一年12月和當(dāng)年1月平均降水量與標(biāo)準(zhǔn)年表進(jìn)行回歸方程的重建,其相關(guān)系數(shù)達(dá)0.68,超過(guò)0.001的顯著性水平,說(shuō)明樹(shù)輪指數(shù)與當(dāng)?shù)厍耙荒?2月至當(dāng)年1月份的平均降水相關(guān).

    建立的線性回歸方程為:Log Y=1.14+2.34Log XSTD(N=50,r=0.64,R2=0.41=0.40,F(xiàn)=34.466,P<0.000 1).式中:Y代表前一年12月和當(dāng)年1月的平均降水量;XSTD代表樹(shù)輪標(biāo)準(zhǔn)年表寬度指數(shù).重建方程的方差解釋量為41%,調(diào)整后的方差解釋量為40%,經(jīng)F檢驗(yàn)達(dá)到的99%的顯著水平.

    對(duì)氣象記錄時(shí)段的觀測(cè)值和重建值進(jìn)行比較(見(jiàn)圖5)可以看出,除個(gè)別時(shí)段有些差距外,其他時(shí)段這兩條曲線的變化趨勢(shì)比較一致.利用逐一剔除法檢測(cè)發(fā)現(xiàn),1992年和1997年的氣象數(shù)據(jù)使回歸方程欠穩(wěn)定,氣象記錄表明,1992年和1997年前一年12月至當(dāng)年1月的降水偏多(1992年為8.75mm,1997年為9.15mm;52a的平均值為3.91mm),一般來(lái)說(shuō)樹(shù)木徑向生長(zhǎng)對(duì)缺水響應(yīng)較好,也就是缺水時(shí)生長(zhǎng)減緩直至停止,而當(dāng)水分增加時(shí)年輪生長(zhǎng)并不能一直隨水分增加而加寬,所以這兩年降水雖多,但重建值并沒(méi)有那么高.如果剔除這兩年,回歸模型的質(zhì)量有所提高(N=48,r=0.68,R2=0.47,=0.46,F(xiàn)=41.975,P<0.000 1).由于1992年和1997年為氣候異常年,從研究長(zhǎng)期趨勢(shì)變化及預(yù)測(cè)的需要等方面考慮,我們將這兩年從模型中剔除.最終建立的回歸方程為:Log Y=0.478+2.376Log XSTD.

    2.4 重建方程的檢驗(yàn)

    為檢驗(yàn)重建方程的可靠性,采用分段方法,用誤差縮減值、符號(hào)檢驗(yàn)和乘積平均值檢驗(yàn)等幾個(gè)統(tǒng)計(jì)量對(duì)重建方程進(jìn)行檢驗(yàn)(見(jiàn)表2).誤差縮減值(RE)是普遍應(yīng)用的精確檢驗(yàn)估計(jì)氣候要素重建值可靠性的統(tǒng)計(jì)量,一般認(rèn)為RE值≥0就算較好地通過(guò)了該項(xiàng)檢驗(yàn),表明該重建方程是穩(wěn)定的,重建值是真實(shí)可信的.本文的結(jié)果(0.50和0.65)表明重建方程是可信的.符號(hào)檢驗(yàn)和乘積平均值檢驗(yàn)均達(dá)到顯著水平.這些統(tǒng)計(jì)檢驗(yàn)結(jié)果表明這個(gè)重建方程是可靠的,可以用于重建過(guò)去200多年來(lái)阿里河地區(qū)前一年12月至當(dāng)年1月的平均降水量.

    圖5 1957—2009年重建和實(shí)測(cè)冬季降水比較

    表2 重建方程的統(tǒng)計(jì)檢驗(yàn)

    2.5 近200年冬季降水變化分析

    根據(jù)上述回歸方程重建阿里河地區(qū)1809年以來(lái)(1742—1808年樣芯序列的復(fù)本量小于5,因此取1809—2009年來(lái)進(jìn)行重建)前一年12月至當(dāng)年1月份平均降水量序列,并對(duì)該序列進(jìn)行了11a的滑動(dòng)平均以獲取低頻變化的信息,結(jié)果見(jiàn)圖6A.

    圖6 阿里河地區(qū)12月—1月平均降水量重建序列(A)與海拉爾地區(qū)1865年以來(lái)的降水變化比較(B)

    從重建的阿里河地區(qū)冬季降水變化可以看出(見(jiàn)圖6A),重建序列表現(xiàn)出4個(gè)降水較多的時(shí)期(高于多年平均值):1837—1848,1867—1869,1932—1965和1988—2002年;同時(shí)也存在4個(gè)降水較少的時(shí)期(低于多年平均降水量):1810—1825,1850—1866,1901—1930和1970—1982年.我們把重建的阿里河前一年12月至當(dāng)年1月的平均降水變化與Liu等[19]用樟子松重建的大興安嶺南段海拉爾地區(qū)1865年以來(lái)前一年7月到當(dāng)年6月的降水變化進(jìn)行了比較(見(jiàn)圖6A和B),結(jié)果表明,二者的豐枯水期吻合的比較好.我們還注意到后2個(gè)干旱時(shí)段(1901—1930年和1970—1982年),以及1867—1887年,1932—1965年和1988—2002年3個(gè)濕潤(rùn)時(shí)期,與黑龍江近200年旱澇變化的干濕時(shí)期有一定的對(duì)應(yīng)[20],但是在起止時(shí)間和持續(xù)時(shí)間上有些差異,并且1872,1927,1949,1954,1980,1982,1988和1998年等極端旱澇年份在樹(shù)輪上也有體現(xiàn)[20-22];通過(guò)查閱歷史資料也發(fā)現(xiàn),1876—1877年和1928—1930年發(fā)生的全國(guó)大范圍的嚴(yán)重干旱事件在樹(shù)輪上也有很好的體現(xiàn)[23],這說(shuō)明阿里河地區(qū)降水的干濕時(shí)期是可信的.

    利用多帶譜分析方法對(duì)阿里河地區(qū)前一年12月到當(dāng)年1月平均降水重建序列進(jìn)行周期分析,以了解過(guò)去200年來(lái)降水變化的周期性(見(jiàn)圖7).通過(guò)分析可以看出,主要的振蕩準(zhǔn)周期有85~38,5.0~4.5和2.2a,均超過(guò)了0.01的顯著水平.氣象學(xué)上的“準(zhǔn)兩年震蕩”(QBO)則與分析得出的2.2a顯著準(zhǔn)周期相一致.

    3 結(jié)論

    通過(guò)比較樹(shù)輪寬度指數(shù)年表之間的參數(shù)特征值,發(fā)現(xiàn)樹(shù)輪對(duì)氣候變化敏感,適合進(jìn)行年輪氣候分析.樟子松標(biāo)準(zhǔn)年表與該地區(qū)前一年12至當(dāng)年1月份平均降水量存在顯著的正相關(guān)關(guān)系,相關(guān)系數(shù)可達(dá)0.68.

    采用線性回歸方程對(duì)阿里河地區(qū)200年來(lái)(1809—2009年)前一年12至當(dāng)年1月平均降水量進(jìn)行了重建,重建方程的方差解釋量為47%(調(diào)整值為46%,F(xiàn)=41.975,P<0.000 1),經(jīng)縮減誤差、符號(hào)檢驗(yàn)和乘積平均值檢驗(yàn),表明重建的降水序列是可信的.

    從整個(gè)重建序列來(lái)看,大興安嶺阿里河地區(qū)近200年冬季降水存在4個(gè)豐水期:1837—1848,1867—1887,1930—1965和1988—2002年;同時(shí)也存在4個(gè)降水較少的時(shí)期:1810—1825,1850—1866,1901—1929和1970—1982年.同時(shí),通過(guò)比對(duì)歷史資料,表明重建的冬季降水序列與黑龍江省近200年來(lái)干濕變化情況比較吻合,并且全國(guó)大范圍的嚴(yán)重干旱事件在樹(shù)輪上也有很好的體現(xiàn).

    多帶譜分析發(fā)現(xiàn),重建的200年來(lái)阿里河地區(qū)前一年12至當(dāng)年1月平均降水量存在85~38,5.0~4.5和2.2a的顯著變化周期.

    圖7 重建降水的多帶譜分析

    [1] 吳祥定.樹(shù)木年輪與氣候變化[M].北京:氣象出版社,1990.

    [2] 李江風(fēng),袁玉江,田希光.樹(shù)木年輪水文學(xué)研究與應(yīng)用[M].北京:科學(xué)出版社,2000.

    [3] GRAUMLICH L J.A 1000-year record of temperature and precipitation in the Sierra Nevada[J].Quaternary Research,1993,39(2):249-255.

    [4] LARA A,VILLALBA R.A 3620-year temperature record from Fitzroya cupressoides tree rings in southern South America[J].Science,1993,260:1104-1106.

    [5] LUCKMAN B H.Tree-ring based reconstruction of summer temperature at the Columbia Icefield Alberta,Canada.A.D.1073-1983[J].Holocene,1997(7):375-389.

    [6] 王亞軍,高尚玉,馬玉貞,等.基于樹(shù)輪重建的寧夏河?xùn)|沙地西部公元1899年來(lái)年降水量的變化[J].干旱區(qū)地理,2010,33(3):377-384.

    [7] 江善虎,任立良,雍斌,等.利用樹(shù)木年輪重建馬鞍山地區(qū)7—9月降水量序列[J].海河大學(xué)學(xué)報(bào):自然科學(xué)版,2010,38(1):31-35.

    [8] 張同文,王麗麗,袁玉江,等.利用樹(shù)輪寬度資料重建天山中段南坡巴侖臺(tái)地區(qū)過(guò)去645年來(lái)的降水變化[J].地理科學(xué),2011,31(2):251-256.

    [9] THOMSEN G.Response to winter precipitation in ring-width chronologies of Pinus sylvestris L.from the northwestern Siberian plain,Russia[J].Tree-Ring Research,2001,57(1):15-29.

    [10] PEDERSON G T,GRAY S T,WOODHOUSE C A,et al.The unusual nature of recent snowpack declines in the North American Cordillera[J].Science,2011,333:332-335.

    [11] LIU X,ZHAO L,CHEN T,et al.Combined tree-ring width andδ13C to reconstruct snowpack depth:apilot study in theGongga Mountain,west China[J].Theoretical Applied Climatology,2011,103:133-144.

    [12] HOLZK?MPER S,KUHRY P,KULTTI S,et al.Stable isotopes in tree rings as proxies for winter precipitation changes in the Russian arctic over the past 150years[J].Geochronometria,2008,32:37-46.

    [13] 徐化成.中國(guó)大興安嶺森林[M].北京:科學(xué)出版社,1998.

    [14] FRITTS H C.Tree rings and climate[M].London:Academic Press,1976.

    [15] HOLMES R L.Computer-assisted quality control in tree-ring dating and measurement[J].Tree-ring Bulletin,1983,43:69-78.

    [16] COOK E R.A time series analysis approach to tree-ring standardization[D].Tucson:University of Arizona,1985.

    [17] COOK E R,MEKO D M,STAHLE D W,et al.Drought reconstruction for the continental United States[J].Journal of Climate,1999,12:1145-1162.

    [18] 袁玉江,李江風(fēng).天山烏魯木齊河源450a冬季溫度序列的重建與分析[J].冰川凍土,1999,21(1):64-70.

    [19] LIU Y,BAO G,SONG H,et al.Precipitation reconstruction from Hailar pine(Pinus sylvestris var.mongolica)tree rings in the Hailar region,Inner Mongolia,China back to 1865AD[J].Palaeogeography Palaeoclimatology Palaeoecology,2009,282:81-87.

    [20] 徐南平,袁美英.近200年來(lái)旱澇變化的探討[J].黑龍江氣象,1994(2):18-21.

    [21] 韓曉君.黑龍江省旱澇特征分析[J].東北水利水電,1999(1):25-26.

    [22] 王艷秋,邢俊江,張麗娟,等.黑龍江旱澇低溫指標(biāo)及其氣候分析應(yīng)用[J].自然災(zāi)害學(xué)報(bào),2008,17(5):142-146.

    [23] 張德二,梁有葉.1876—1878年中國(guó)大范圍持續(xù)干旱事件[J].氣候變化研究進(jìn)展,2010,6(2):106-112.

    Growth-climate response and winter precipitation reconstruction of Pinus sylvestris var.mongolicain A'li River of Greater Khingan Range

    LV Shan-na,WANG Xiao-chun
    (Center for Ecological Research,Northeast Forestry University,Harbin 150040,China)

    Greater Khingan Range is one of most important forest area in northeast China.It has critical significance for forestry management and production to understand long-term climate conditions in this region.Tree-ring chronology of Pinus sylvestris var.mongolica in A'li River of Greater Khingan Range was significantly correlated with the precipitation in current January and previous December by using response function analysis.Therefore,winter precipitations based on tree-ring records in A'li River since 1809were reconstructed,which explained 47%of the variance in winter precipitation.Reconstructed winter precipitation series was coincided with wet-dry changes in Heilongjiang Province.Winter precipitation in A'li River for last 200years exhibited four dry and wet periods.In addition,the winter precipitation also showed that climate in this area experienced alternating dry and wet periods,including dry periods in the early 19thcentury,wet periods in the mid-19thcentury and dry periods in the late 19thcentury to the early 20thcentury.While,climate was relatively dry in the early 20thcentury,wet in the middle 20thcentury and dry at the end of 20thcentury.Reconstructed winter precipitations exhibited significant 85~38,5.0~4.5and 2.2years periodicities by using multi-taper methods.

    A'li River;Pinus sylvestris var.mongolica;winter precipitation;tree rings

    S 718.45 [學(xué)科代碼] 220·1060

    A

    (責(zé)任編輯:方 林)

    1000-1832(2014)02-00110-07

    10.11672/dbsdzk2014-01-022

    2013-04-22

    國(guó)家自然科學(xué)基金資助項(xiàng)目(30970481);中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)基金資助項(xiàng)目(DL13EA05-02);長(zhǎng)江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃資助項(xiàng)目(IRT1054);黑龍江省留學(xué)歸國(guó)基金資助項(xiàng)目(LC2012C09).

    呂姍娜(1984—),女,碩士研究生;通訊作者:王曉春(1975—),男,博士,教授,主要從事樹(shù)木年輪學(xué)和全球生態(tài)學(xué)研究.

    猜你喜歡
    樹(shù)輪年表樟子松
    年表
    世界建筑(2022年10期)2022-10-25 08:26:48
    樹(shù)輪研究方法與進(jìn)展
    天山南北坡樹(shù)輪穩(wěn)定碳同位素對(duì)氣候的響應(yīng)差異
    王錫良年表
    石魯年表
    塞罕壩樟子松幼林撫育與管理
    初探北方樟子松栽培關(guān)鍵技術(shù)
    我眼中的樟子松
    北極光(2018年12期)2018-03-07 01:01:58
    樹(shù)木年輪對(duì)生態(tài)環(huán)境變化的響應(yīng)及應(yīng)用研究
    阿爾泰山中部樹(shù)輪寬度年表特征及其氣候響應(yīng)分析
    国产亚洲精品久久久久5区| 亚洲男人天堂网一区| 国内毛片毛片毛片毛片毛片| 午夜精品在线福利| 成人18禁高潮啪啪吃奶动态图| 亚洲专区字幕在线| 国产99白浆流出| 亚洲国产看品久久| 日韩视频一区二区在线观看| 午夜精品久久久久久毛片777| 国产精品永久免费网站| 变态另类成人亚洲欧美熟女 | 母亲3免费完整高清在线观看| 最新在线观看一区二区三区| 精品久久久精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区在线臀色熟女| 国产成人av教育| 精品国产乱子伦一区二区三区| 黑人巨大精品欧美一区二区mp4| 日韩精品免费视频一区二区三区| 母亲3免费完整高清在线观看| 欧美成人免费av一区二区三区| 搡老熟女国产l中国老女人| 亚洲全国av大片| 欧美黄色片欧美黄色片| 久久午夜亚洲精品久久| 亚洲欧美日韩高清在线视频| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看| 国产区一区二久久| 搞女人的毛片| av在线播放免费不卡| 精品免费久久久久久久清纯| 成熟少妇高潮喷水视频| 亚洲午夜精品一区,二区,三区| 丁香欧美五月| 久久亚洲真实| xxx96com| 日韩免费av在线播放| 免费高清视频大片| 人人妻人人澡人人看| 久久青草综合色| 身体一侧抽搐| 亚洲av日韩精品久久久久久密| 99精品久久久久人妻精品| 我的亚洲天堂| 精品国产超薄肉色丝袜足j| 午夜免费观看网址| 国产熟女午夜一区二区三区| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av| 精品欧美一区二区三区在线| 国产精品98久久久久久宅男小说| 欧美黄色片欧美黄色片| 人人妻人人澡人人看| 免费一级毛片在线播放高清视频 | e午夜精品久久久久久久| 亚洲午夜理论影院| 国产精品久久久av美女十八| 午夜福利在线观看吧| 日本在线视频免费播放| 免费在线观看视频国产中文字幕亚洲| 757午夜福利合集在线观看| 999久久久精品免费观看国产| 欧美在线黄色| 亚洲熟妇中文字幕五十中出| 日本a在线网址| 国产国语露脸激情在线看| 国产精华一区二区三区| 狠狠狠狠99中文字幕| av网站免费在线观看视频| 中出人妻视频一区二区| 一级a爱视频在线免费观看| avwww免费| 老鸭窝网址在线观看| 法律面前人人平等表现在哪些方面| 成年版毛片免费区| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器 | 精品人妻在线不人妻| 69av精品久久久久久| 女人被狂操c到高潮| 色婷婷久久久亚洲欧美| 国产91精品成人一区二区三区| 一本久久中文字幕| 麻豆一二三区av精品| 日韩欧美国产一区二区入口| 亚洲九九香蕉| 国产精品 欧美亚洲| 在线观看免费视频日本深夜| 国产精品九九99| 免费在线观看亚洲国产| 亚洲色图av天堂| 午夜福利一区二区在线看| 久久精品成人免费网站| 看片在线看免费视频| 欧美+亚洲+日韩+国产| 国产av精品麻豆| 欧美丝袜亚洲另类 | 日本 欧美在线| 亚洲天堂国产精品一区在线| 91字幕亚洲| 成在线人永久免费视频| 人妻久久中文字幕网| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 少妇的丰满在线观看| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| www.www免费av| 我的亚洲天堂| 村上凉子中文字幕在线| 国产精品久久视频播放| 亚洲天堂国产精品一区在线| 韩国精品一区二区三区| 免费观看人在逋| 人妻久久中文字幕网| 一级黄色大片毛片| 成人免费观看视频高清| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区| 亚洲人成77777在线视频| 欧美黄色淫秽网站| 亚洲熟女毛片儿| 在线天堂中文资源库| aaaaa片日本免费| 欧美绝顶高潮抽搐喷水| 在线免费观看的www视频| 久久香蕉激情| 亚洲av成人一区二区三| 巨乳人妻的诱惑在线观看| 午夜福利在线观看吧| 日韩大尺度精品在线看网址 | 成熟少妇高潮喷水视频| 久久久久久人人人人人| 国产伦人伦偷精品视频| 淫妇啪啪啪对白视频| 97碰自拍视频| 欧美在线一区亚洲| 国产激情欧美一区二区| 精品国产一区二区三区四区第35| 波多野结衣一区麻豆| 午夜日韩欧美国产| 欧美午夜高清在线| 两人在一起打扑克的视频| 又黄又爽又免费观看的视频| 亚洲,欧美精品.| 操出白浆在线播放| 欧美日韩瑟瑟在线播放| 19禁男女啪啪无遮挡网站| 夜夜夜夜夜久久久久| 激情在线观看视频在线高清| 国产亚洲欧美在线一区二区| 成人18禁在线播放| 国产男靠女视频免费网站| 欧美激情 高清一区二区三区| 又大又爽又粗| 人妻丰满熟妇av一区二区三区| 日韩精品中文字幕看吧| 国产一区二区在线av高清观看| 国产亚洲欧美98| 九色亚洲精品在线播放| 国产一区二区三区在线臀色熟女| 欧美久久黑人一区二区| 亚洲国产精品999在线| 午夜福利在线观看吧| 日日夜夜操网爽| 午夜福利18| 老鸭窝网址在线观看| 中文字幕久久专区| 久久精品国产亚洲av香蕉五月| 欧美成狂野欧美在线观看| 99久久99久久久精品蜜桃| 国产精品香港三级国产av潘金莲| 非洲黑人性xxxx精品又粗又长| 欧美国产精品va在线观看不卡| 久久精品国产清高在天天线| 首页视频小说图片口味搜索| 国产高清视频在线播放一区| 欧美成人一区二区免费高清观看 | 18禁美女被吸乳视频| 正在播放国产对白刺激| 久久人人精品亚洲av| 这个男人来自地球电影免费观看| 咕卡用的链子| 他把我摸到了高潮在线观看| 亚洲免费av在线视频| 在线观看66精品国产| 黄色视频不卡| 正在播放国产对白刺激| 亚洲国产精品999在线| 91国产中文字幕| 99久久99久久久精品蜜桃| 男女做爰动态图高潮gif福利片 | 欧美一区二区精品小视频在线| 亚洲精品中文字幕一二三四区| 黑人欧美特级aaaaaa片| 国产成人系列免费观看| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| or卡值多少钱| 国产精品一区二区在线不卡| 欧美av亚洲av综合av国产av| 别揉我奶头~嗯~啊~动态视频| 国产精品乱码一区二三区的特点 | 搡老熟女国产l中国老女人| 精品国产美女av久久久久小说| 99久久国产精品久久久| 国产午夜福利久久久久久| av在线天堂中文字幕| 精品久久久久久成人av| 免费在线观看视频国产中文字幕亚洲| 欧美精品亚洲一区二区| 在线观看一区二区三区| 老鸭窝网址在线观看| 操美女的视频在线观看| 国产成人精品久久二区二区免费| 悠悠久久av| 91麻豆精品激情在线观看国产| 男人舔女人的私密视频| 亚洲国产看品久久| 国产欧美日韩一区二区精品| av免费在线观看网站| 视频区欧美日本亚洲| 欧美日韩亚洲国产一区二区在线观看| 免费看a级黄色片| 日韩高清综合在线| 久久久久国产精品人妻aⅴ院| 岛国视频午夜一区免费看| 少妇的丰满在线观看| 麻豆久久精品国产亚洲av| 日本a在线网址| 国产精品电影一区二区三区| 国产91精品成人一区二区三区| 成人国语在线视频| 亚洲全国av大片| 视频在线观看一区二区三区| www日本在线高清视频| 十八禁网站免费在线| 成熟少妇高潮喷水视频| 国产成人欧美| 又大又爽又粗| 又黄又粗又硬又大视频| 久久久久久久午夜电影| 精品国产乱子伦一区二区三区| 中文字幕人成人乱码亚洲影| 亚洲一区高清亚洲精品| 亚洲专区中文字幕在线| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 亚洲无线在线观看| av天堂久久9| 国产野战对白在线观看| 久久这里只有精品19| 在线观看免费视频日本深夜| 久久久久久久久久久久大奶| 亚洲国产日韩欧美精品在线观看 | 日韩高清综合在线| 国产精品一区二区精品视频观看| 黄频高清免费视频| 91在线观看av| 久久久久久大精品| 脱女人内裤的视频| 亚洲,欧美精品.| 亚洲av电影不卡..在线观看| 色在线成人网| 97超级碰碰碰精品色视频在线观看| 精品一区二区三区视频在线观看免费| 亚洲五月天丁香| 好看av亚洲va欧美ⅴa在| www.999成人在线观看| 亚洲av电影不卡..在线观看| 国产在线观看jvid| 成人国语在线视频| 国产亚洲精品第一综合不卡| 无人区码免费观看不卡| ponron亚洲| 丝袜人妻中文字幕| 国产精品99久久99久久久不卡| 久久人妻福利社区极品人妻图片| 日本在线视频免费播放| 国产一区二区三区视频了| 免费看十八禁软件| 级片在线观看| 欧美日韩黄片免| 国产精品自产拍在线观看55亚洲| 很黄的视频免费| 九色亚洲精品在线播放| 宅男免费午夜| 午夜精品在线福利| 国产成人啪精品午夜网站| 亚洲电影在线观看av| 99国产精品免费福利视频| 精品国产一区二区久久| 国产精品乱码一区二三区的特点 | 少妇的丰满在线观看| 亚洲成人免费电影在线观看| 久久久久久久午夜电影| 中文字幕色久视频| 欧美性长视频在线观看| 美女午夜性视频免费| 丝袜人妻中文字幕| 脱女人内裤的视频| 波多野结衣av一区二区av| 国产高清有码在线观看视频 | 不卡av一区二区三区| 嫩草影视91久久| 国产伦人伦偷精品视频| 国产区一区二久久| 又紧又爽又黄一区二区| 99国产精品免费福利视频| 国产精品秋霞免费鲁丝片| 国产精品亚洲美女久久久| 日韩成人在线观看一区二区三区| 国产伦人伦偷精品视频| 午夜成年电影在线免费观看| 999久久久精品免费观看国产| 久久精品国产亚洲av高清一级| 欧美最黄视频在线播放免费| 日本欧美视频一区| 91国产中文字幕| 韩国av一区二区三区四区| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 在线国产一区二区在线| 国产主播在线观看一区二区| 久久久国产成人免费| 欧美乱色亚洲激情| 男人舔女人下体高潮全视频| 成年人黄色毛片网站| 看黄色毛片网站| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 99久久99久久久精品蜜桃| 国产xxxxx性猛交| 午夜福利视频1000在线观看 | 亚洲av日韩精品久久久久久密| 久久久久久亚洲精品国产蜜桃av| 女同久久另类99精品国产91| 国产精品久久久久久精品电影 | av福利片在线| 国产免费男女视频| 国产欧美日韩一区二区三| 人人妻人人澡人人看| 色老头精品视频在线观看| 深夜精品福利| 在线观看免费午夜福利视频| 亚洲专区字幕在线| 国产男靠女视频免费网站| 最近最新免费中文字幕在线| 欧美最黄视频在线播放免费| 免费少妇av软件| 亚洲精品久久国产高清桃花| 麻豆一二三区av精品| 国产野战对白在线观看| 日日干狠狠操夜夜爽| 国产高清videossex| 国产一级毛片七仙女欲春2 | 老鸭窝网址在线观看| 99精品在免费线老司机午夜| 天天躁夜夜躁狠狠躁躁| 狠狠狠狠99中文字幕| 国产精品永久免费网站| 黄色 视频免费看| 国产不卡一卡二| 亚洲少妇的诱惑av| 日韩欧美免费精品| 国产成人啪精品午夜网站| 久久精品91蜜桃| 久久精品影院6| 99久久国产精品久久久| 一边摸一边抽搐一进一出视频| 亚洲av美国av| 久99久视频精品免费| 99国产精品99久久久久| 亚洲成人精品中文字幕电影| 高清黄色对白视频在线免费看| 老司机午夜福利在线观看视频| 欧美精品亚洲一区二区| 午夜免费观看网址| 欧美色欧美亚洲另类二区 | 性色av乱码一区二区三区2| 校园春色视频在线观看| 一边摸一边做爽爽视频免费| 免费av毛片视频| 国产精品秋霞免费鲁丝片| 美国免费a级毛片| 亚洲片人在线观看| a在线观看视频网站| 99精品在免费线老司机午夜| 黑人巨大精品欧美一区二区蜜桃| 性色av乱码一区二区三区2| 欧美一级毛片孕妇| 亚洲人成网站在线播放欧美日韩| 亚洲av成人av| 99在线视频只有这里精品首页| 亚洲中文日韩欧美视频| 国产又爽黄色视频| 一级,二级,三级黄色视频| 亚洲情色 制服丝袜| 91字幕亚洲| cao死你这个sao货| 精品国产乱子伦一区二区三区| 欧美日本亚洲视频在线播放| 日本免费a在线| 精品久久久久久成人av| 一边摸一边抽搐一进一出视频| 欧美午夜高清在线| 黄色片一级片一级黄色片| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| 禁无遮挡网站| 国产成+人综合+亚洲专区| 欧美av亚洲av综合av国产av| av视频在线观看入口| 亚洲av电影不卡..在线观看| 成人特级黄色片久久久久久久| 男女床上黄色一级片免费看| 极品人妻少妇av视频| 亚洲午夜精品一区,二区,三区| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 亚洲av成人一区二区三| 日韩 欧美 亚洲 中文字幕| 国产亚洲av嫩草精品影院| 亚洲av成人一区二区三| 婷婷六月久久综合丁香| 视频在线观看一区二区三区| 黄色a级毛片大全视频| 最好的美女福利视频网| 国产精品久久视频播放| 色综合婷婷激情| 免费少妇av软件| 国产精品久久久久久人妻精品电影| 级片在线观看| 亚洲国产欧美日韩在线播放| 久久香蕉国产精品| 国产精品久久电影中文字幕| 一边摸一边抽搐一进一出视频| 国产精品综合久久久久久久免费 | 满18在线观看网站| 日本在线视频免费播放| 在线观看www视频免费| 夜夜爽天天搞| 欧美日本视频| 国产精品免费一区二区三区在线| 久久午夜亚洲精品久久| 国产精品亚洲美女久久久| 人人妻人人爽人人添夜夜欢视频| 久久久久久亚洲精品国产蜜桃av| 一本大道久久a久久精品| 人人妻人人澡欧美一区二区 | 又黄又粗又硬又大视频| av视频免费观看在线观看| 国产97色在线日韩免费| 亚洲精品国产区一区二| 亚洲九九香蕉| 黑人欧美特级aaaaaa片| 天天躁狠狠躁夜夜躁狠狠躁| 久久性视频一级片| 熟妇人妻久久中文字幕3abv| 亚洲成人国产一区在线观看| 操美女的视频在线观看| 黄网站色视频无遮挡免费观看| 女性生殖器流出的白浆| 色综合欧美亚洲国产小说| 身体一侧抽搐| 夜夜爽天天搞| 看黄色毛片网站| 亚洲第一青青草原| 久久久久久亚洲精品国产蜜桃av| 亚洲中文字幕日韩| 欧美成人性av电影在线观看| 免费观看精品视频网站| 少妇被粗大的猛进出69影院| bbb黄色大片| 少妇 在线观看| 成人免费观看视频高清| 欧美中文综合在线视频| 自线自在国产av| 亚洲五月天丁香| 女性被躁到高潮视频| 国产三级黄色录像| 国产精品一区二区免费欧美| 一区二区日韩欧美中文字幕| 99精品久久久久人妻精品| 久久久久国产精品人妻aⅴ院| 黄色视频,在线免费观看| 看免费av毛片| 宅男免费午夜| www国产在线视频色| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕大全免费视频| 一级毛片女人18水好多| 一进一出好大好爽视频| 精品日产1卡2卡| 身体一侧抽搐| 激情视频va一区二区三区| 亚洲七黄色美女视频| 脱女人内裤的视频| 一进一出好大好爽视频| 老司机午夜福利在线观看视频| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看 | 国产在线观看jvid| 后天国语完整版免费观看| 亚洲五月色婷婷综合| 日韩av在线大香蕉| 日韩国内少妇激情av| 国产精品精品国产色婷婷| 中文字幕高清在线视频| www.999成人在线观看| 亚洲欧洲精品一区二区精品久久久| 国产又色又爽无遮挡免费看| 欧美色欧美亚洲另类二区 | 日本 欧美在线| 免费在线观看视频国产中文字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 国产单亲对白刺激| 久久人人精品亚洲av| 欧美人与性动交α欧美精品济南到| 国产精品1区2区在线观看.| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| 欧美激情 高清一区二区三区| 久久午夜综合久久蜜桃| 男女做爰动态图高潮gif福利片 | 亚洲欧美日韩高清在线视频| 国产免费av片在线观看野外av| 亚洲男人的天堂狠狠| 99re在线观看精品视频| 欧美成人免费av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 极品人妻少妇av视频| 大型黄色视频在线免费观看| 一区在线观看完整版| 国产99白浆流出| 别揉我奶头~嗯~啊~动态视频| 亚洲一区二区三区色噜噜| 中国美女看黄片| 欧美不卡视频在线免费观看 | 欧美精品亚洲一区二区| 极品教师在线免费播放| 热re99久久国产66热| АⅤ资源中文在线天堂| 久久国产精品人妻蜜桃| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 日韩成人在线观看一区二区三区| 人人妻人人澡人人看| 久久婷婷成人综合色麻豆| 精品欧美国产一区二区三| 中文亚洲av片在线观看爽| 十八禁人妻一区二区| 一二三四在线观看免费中文在| 午夜福利在线观看吧| 国产精品免费视频内射| 男人的好看免费观看在线视频 | 99在线人妻在线中文字幕| 日本三级黄在线观看| 巨乳人妻的诱惑在线观看| www国产在线视频色| av天堂在线播放| 欧美色视频一区免费| 国产精品综合久久久久久久免费 | a级毛片在线看网站| 久久热在线av| 国产91精品成人一区二区三区| 日本三级黄在线观看| 久久狼人影院| 亚洲精品国产一区二区精华液| 国产不卡一卡二| 极品人妻少妇av视频| 999久久久精品免费观看国产| 欧美人与性动交α欧美精品济南到| 亚洲熟女毛片儿| 高清毛片免费观看视频网站| 老熟妇乱子伦视频在线观看| 日本 欧美在线| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸| 免费观看精品视频网站| 国产精品爽爽va在线观看网站 | 亚洲精品久久国产高清桃花| 精品不卡国产一区二区三区| 日本 欧美在线| 9热在线视频观看99| av超薄肉色丝袜交足视频| 在线观看66精品国产| 大型av网站在线播放| 丁香六月欧美| 久久精品影院6| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 在线国产一区二区在线| 国产成人精品在线电影| 精品久久久久久久久久免费视频| 国产欧美日韩一区二区三| 天堂影院成人在线观看| 一二三四社区在线视频社区8| 欧美大码av| avwww免费| 欧美另类亚洲清纯唯美| 日韩欧美一区视频在线观看| 精品欧美一区二区三区在线| 国产真人三级小视频在线观看| 日韩欧美在线二视频| 在线观看免费视频网站a站| 久久草成人影院| 久久婷婷人人爽人人干人人爱 | 日韩成人在线观看一区二区三区| 婷婷丁香在线五月| 亚洲欧美日韩无卡精品|