• <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)分析
    日韩精品青青久久久久久| 午夜精品久久久久久毛片777| 午夜成年电影在线免费观看| 成人永久免费在线观看视频| 黄色 视频免费看| 香蕉av资源在线| 两个人视频免费观看高清| 久久久精品国产亚洲av高清涩受| 国产成人av激情在线播放| av福利片在线观看| 国产又色又爽无遮挡免费看| 岛国在线免费视频观看| 一卡2卡三卡四卡精品乱码亚洲| 成人欧美大片| 不卡一级毛片| 黄色视频不卡| 后天国语完整版免费观看| 国产片内射在线| 老熟妇乱子伦视频在线观看| 亚洲国产看品久久| 亚洲精华国产精华精| 精品欧美一区二区三区在线| 又黄又爽又免费观看的视频| 夜夜看夜夜爽夜夜摸| www.自偷自拍.com| 最近视频中文字幕2019在线8| 在线观看日韩欧美| 久久久久久久久中文| 少妇粗大呻吟视频| 午夜免费观看网址| 久久香蕉国产精品| 一级片免费观看大全| 欧美一级毛片孕妇| 曰老女人黄片| 日韩欧美一区二区三区在线观看| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| av福利片在线观看| 搡老岳熟女国产| 草草在线视频免费看| 日韩欧美免费精品| 午夜免费激情av| 熟女少妇亚洲综合色aaa.| 超碰成人久久| av在线播放免费不卡| 性色av乱码一区二区三区2| 97超级碰碰碰精品色视频在线观看| 久久 成人 亚洲| 每晚都被弄得嗷嗷叫到高潮| 九色成人免费人妻av| 午夜精品在线福利| 日韩大尺度精品在线看网址| 国内少妇人妻偷人精品xxx网站 | 搡老妇女老女人老熟妇| 性欧美人与动物交配| 一本综合久久免费| 中文字幕人成人乱码亚洲影| 久久精品成人免费网站| 亚洲在线自拍视频| 精品国产亚洲在线| 亚洲av日韩精品久久久久久密| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添小说| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕人妻熟女| 首页视频小说图片口味搜索| 国产精品一区二区三区四区免费观看 | 午夜免费激情av| 99久久无色码亚洲精品果冻| 两性夫妻黄色片| 国产69精品久久久久777片 | 波多野结衣高清无吗| 九九热线精品视视频播放| 精品一区二区三区四区五区乱码| 看黄色毛片网站| 亚洲av成人不卡在线观看播放网| 久久人妻福利社区极品人妻图片| 99热这里只有精品一区 | 亚洲成人久久性| 亚洲国产欧美人成| 又爽又黄无遮挡网站| 757午夜福利合集在线观看| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 成年女人毛片免费观看观看9| 欧美绝顶高潮抽搐喷水| 国产高清激情床上av| 国产精品 国内视频| 国产人伦9x9x在线观看| 久久婷婷成人综合色麻豆| 国产av不卡久久| 女警被强在线播放| 久久精品国产综合久久久| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 国产精品亚洲av一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 欧美在线黄色| 欧美zozozo另类| 亚洲男人的天堂狠狠| 在线观看66精品国产| 91麻豆精品激情在线观看国产| 男人的好看免费观看在线视频 | 国产精品精品国产色婷婷| 欧美黄色淫秽网站| 女同久久另类99精品国产91| 久久久久久大精品| 国产精品久久电影中文字幕| 三级国产精品欧美在线观看 | 久久这里只有精品19| 国产亚洲精品综合一区在线观看 | 亚洲天堂国产精品一区在线| 18禁观看日本| 国产成人精品久久二区二区91| 欧美中文综合在线视频| 99久久久亚洲精品蜜臀av| 久久久久久久午夜电影| 校园春色视频在线观看| 亚洲熟女毛片儿| 亚洲精品美女久久av网站| 18禁黄网站禁片午夜丰满| 亚洲av片天天在线观看| av欧美777| 男人的好看免费观看在线视频 | av视频在线观看入口| 欧美国产日韩亚洲一区| 身体一侧抽搐| 18禁美女被吸乳视频| www.精华液| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 正在播放国产对白刺激| 每晚都被弄得嗷嗷叫到高潮| 成人午夜高清在线视频| 色尼玛亚洲综合影院| 动漫黄色视频在线观看| 99国产综合亚洲精品| 婷婷丁香在线五月| 黄色成人免费大全| 久久久国产精品麻豆| 国产单亲对白刺激| 亚洲人成网站在线播放欧美日韩| 亚洲人成77777在线视频| 最新在线观看一区二区三区| 国产精品亚洲一级av第二区| 俄罗斯特黄特色一大片| 精品久久久久久成人av| 两个人的视频大全免费| 欧美极品一区二区三区四区| 在线a可以看的网站| 天堂动漫精品| 欧美成人免费av一区二区三区| 一二三四在线观看免费中文在| 亚洲av美国av| 中出人妻视频一区二区| 我的老师免费观看完整版| 露出奶头的视频| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| www.熟女人妻精品国产| 婷婷精品国产亚洲av在线| 国产成人精品无人区| 制服人妻中文乱码| 免费看a级黄色片| 中亚洲国语对白在线视频| 窝窝影院91人妻| 国产成人系列免费观看| 麻豆国产av国片精品| 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 精品午夜福利视频在线观看一区| 欧美三级亚洲精品| 性色av乱码一区二区三区2| 久久午夜综合久久蜜桃| 午夜福利18| 久久天躁狠狠躁夜夜2o2o| 丁香六月欧美| 久久婷婷人人爽人人干人人爱| 好男人在线观看高清免费视频| 亚洲国产欧美网| 国产不卡一卡二| 欧美成狂野欧美在线观看| 国产亚洲精品综合一区在线观看 | 国产激情欧美一区二区| 他把我摸到了高潮在线观看| 桃色一区二区三区在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲性夜色夜夜综合| 午夜成年电影在线免费观看| 国产三级在线视频| 一级作爱视频免费观看| 日本成人三级电影网站| 性色av乱码一区二区三区2| xxx96com| 午夜福利欧美成人| 精品一区二区三区四区五区乱码| 国产精品一区二区三区四区免费观看 | 三级国产精品欧美在线观看 | 午夜福利成人在线免费观看| 岛国在线观看网站| 久久精品人妻少妇| 国产午夜精品论理片| 欧美成人午夜精品| 村上凉子中文字幕在线| 天堂动漫精品| 久久午夜亚洲精品久久| 女同久久另类99精品国产91| av国产免费在线观看| 51午夜福利影视在线观看| 91麻豆av在线| 天天一区二区日本电影三级| 国产又色又爽无遮挡免费看| 亚洲精华国产精华精| 精品第一国产精品| 欧美三级亚洲精品| 日本精品一区二区三区蜜桃| av视频在线观看入口| 亚洲欧洲精品一区二区精品久久久| 日韩国内少妇激情av| 亚洲美女视频黄频| 久久久久九九精品影院| 99精品在免费线老司机午夜| 国产三级中文精品| 精品免费久久久久久久清纯| 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 精品第一国产精品| www.熟女人妻精品国产| 亚洲欧美日韩东京热| 午夜福利欧美成人| 制服丝袜大香蕉在线| 可以免费在线观看a视频的电影网站| 欧美日韩国产亚洲二区| 成人国产一区最新在线观看| 一个人免费在线观看的高清视频| 黄色片一级片一级黄色片| 国产亚洲精品综合一区在线观看 | 国产精品久久视频播放| 欧美中文综合在线视频| 国产精品国产高清国产av| 亚洲精品久久成人aⅴ小说| 亚洲七黄色美女视频| 可以免费在线观看a视频的电影网站| 超碰成人久久| 五月伊人婷婷丁香| 99国产综合亚洲精品| e午夜精品久久久久久久| 亚洲欧美日韩高清在线视频| 国产亚洲av高清不卡| 国产欧美日韩一区二区三| 国产成人精品无人区| 首页视频小说图片口味搜索| 亚洲av片天天在线观看| 国产一区二区三区视频了| 久久久久久久久久黄片| 一个人观看的视频www高清免费观看 | 国产99白浆流出| 日本五十路高清| 亚洲自拍偷在线| 久久久国产精品麻豆| 最近最新免费中文字幕在线| 少妇的丰满在线观看| 日日爽夜夜爽网站| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 久久久久久久久久黄片| 国产成人精品无人区| 很黄的视频免费| 国产精品野战在线观看| 搞女人的毛片| 国内少妇人妻偷人精品xxx网站 | 我要搜黄色片| 午夜精品在线福利| www国产在线视频色| 国产精品电影一区二区三区| 久久欧美精品欧美久久欧美| 国产精品久久久久久精品电影| 色综合站精品国产| 久久亚洲真实| 两人在一起打扑克的视频| 日本熟妇午夜| 国产aⅴ精品一区二区三区波| 欧美性长视频在线观看| 国产精品,欧美在线| 日本一区二区免费在线视频| 超碰成人久久| 久久精品人妻少妇| 长腿黑丝高跟| 男女之事视频高清在线观看| 国产精品美女特级片免费视频播放器 | 少妇被粗大的猛进出69影院| 国产成人精品无人区| 精品一区二区三区四区五区乱码| 国产1区2区3区精品| 国产亚洲精品久久久久久毛片| 欧美av亚洲av综合av国产av| 亚洲av电影在线进入| 后天国语完整版免费观看| 欧美大码av| 国产精品1区2区在线观看.| 国产精品香港三级国产av潘金莲| 男女之事视频高清在线观看| 欧美人与性动交α欧美精品济南到| 人妻丰满熟妇av一区二区三区| 久久国产精品影院| 久久99热这里只有精品18| 日本 欧美在线| 欧美成人免费av一区二区三区| 国产成+人综合+亚洲专区| 久久久久久人人人人人| 黄片小视频在线播放| 日本a在线网址| 国产av又大| 黄色a级毛片大全视频| 精品不卡国产一区二区三区| 亚洲欧美日韩无卡精品| 久久人妻福利社区极品人妻图片| 精品欧美一区二区三区在线| 日日摸夜夜添夜夜添小说| 国产一区二区三区视频了| 午夜免费观看网址| 日本a在线网址| 中出人妻视频一区二区| 99国产精品一区二区三区| 男女午夜视频在线观看| 久久久久久久午夜电影| 在线永久观看黄色视频| 婷婷精品国产亚洲av在线| 久久中文看片网| 1024手机看黄色片| 最新美女视频免费是黄的| 亚洲熟女毛片儿| 哪里可以看免费的av片| 国产成人影院久久av| 亚洲欧洲精品一区二区精品久久久| 成人特级黄色片久久久久久久| 欧美黑人欧美精品刺激| 亚洲色图 男人天堂 中文字幕| 久久国产精品人妻蜜桃| 一区二区三区激情视频| 亚洲熟妇中文字幕五十中出| 别揉我奶头~嗯~啊~动态视频| 一级毛片女人18水好多| 国内少妇人妻偷人精品xxx网站 | 日韩精品青青久久久久久| 99国产综合亚洲精品| 91麻豆av在线| 亚洲欧美日韩东京热| 少妇被粗大的猛进出69影院| 一边摸一边抽搐一进一小说| 日本一本二区三区精品| 人妻夜夜爽99麻豆av| 色综合站精品国产| av片东京热男人的天堂| 国产亚洲精品一区二区www| av欧美777| 男人的好看免费观看在线视频 | 色综合欧美亚洲国产小说| 特级一级黄色大片| 欧美黑人精品巨大| 免费看十八禁软件| 级片在线观看| 伊人久久大香线蕉亚洲五| 99国产精品一区二区蜜桃av| 精品欧美国产一区二区三| 欧美绝顶高潮抽搐喷水| 97人妻精品一区二区三区麻豆| 亚洲国产精品成人综合色| 精品第一国产精品| 嫩草影视91久久| 在线观看美女被高潮喷水网站 | 毛片女人毛片| 一本大道久久a久久精品| av欧美777| 国产欧美日韩一区二区精品| 搡老熟女国产l中国老女人| 国产一区在线观看成人免费| 两个人视频免费观看高清| www.精华液| 老汉色∧v一级毛片| 国产精品爽爽va在线观看网站| 老司机福利观看| 亚洲中文字幕日韩| 一个人免费在线观看的高清视频| 国产激情久久老熟女| 色哟哟哟哟哟哟| 精品国产乱码久久久久久男人| 亚洲天堂国产精品一区在线| 男人舔奶头视频| 国产精品久久久久久人妻精品电影| 50天的宝宝边吃奶边哭怎么回事| 国产男靠女视频免费网站| 亚洲欧美一区二区三区黑人| xxx96com| 亚洲成人精品中文字幕电影| a级毛片在线看网站| 国产精品永久免费网站| 一进一出抽搐gif免费好疼| 国产亚洲精品久久久久5区| √禁漫天堂资源中文www| 国产午夜福利久久久久久| 亚洲av成人av| 可以免费在线观看a视频的电影网站| 精品国产亚洲在线| 亚洲人成77777在线视频| √禁漫天堂资源中文www| 日本免费一区二区三区高清不卡| av视频在线观看入口| 亚洲在线自拍视频| 成人国产一区最新在线观看| 久久中文字幕人妻熟女| 国产久久久一区二区三区| 国产精品久久久久久久电影 | 亚洲专区字幕在线| 国产精品久久久人人做人人爽| 美女 人体艺术 gogo| 成年版毛片免费区| 长腿黑丝高跟| 日本撒尿小便嘘嘘汇集6| 无限看片的www在线观看| 精品乱码久久久久久99久播| 黄色片一级片一级黄色片| 国产精品乱码一区二三区的特点| 久久天堂一区二区三区四区| 99久久国产精品久久久| 曰老女人黄片| 久久精品影院6| 中出人妻视频一区二区| 夜夜看夜夜爽夜夜摸| 中文字幕熟女人妻在线| 法律面前人人平等表现在哪些方面| 老汉色∧v一级毛片| 欧美日韩乱码在线| 最近最新中文字幕大全电影3| 久久久久亚洲av毛片大全| 精品高清国产在线一区| 又粗又爽又猛毛片免费看| 51午夜福利影视在线观看| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线在线| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 午夜a级毛片| 一进一出抽搐动态| 欧美日韩精品网址| 床上黄色一级片| 欧美极品一区二区三区四区| 色哟哟哟哟哟哟| 亚洲一码二码三码区别大吗| 黄色毛片三级朝国网站| 日韩 欧美 亚洲 中文字幕| 欧美日韩一级在线毛片| 777久久人妻少妇嫩草av网站| 琪琪午夜伦伦电影理论片6080| 欧美成人免费av一区二区三区| 色哟哟哟哟哟哟| 久久久久久久精品吃奶| 波多野结衣巨乳人妻| 国产精品av视频在线免费观看| 在线观看免费午夜福利视频| 国产精品免费视频内射| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久电影中文字幕| 日本免费一区二区三区高清不卡| 久久久久久久精品吃奶| 亚洲av中文字字幕乱码综合| 99久久精品热视频| 老司机午夜十八禁免费视频| 小说图片视频综合网站| 日韩欧美免费精品| 91大片在线观看| 美女大奶头视频| 给我免费播放毛片高清在线观看| 久久精品国产综合久久久| 日韩欧美一区二区三区在线观看| 无人区码免费观看不卡| 国产精品98久久久久久宅男小说| 一个人免费在线观看的高清视频| 99热这里只有是精品50| 国产男靠女视频免费网站| 亚洲精品久久国产高清桃花| 免费在线观看视频国产中文字幕亚洲| 国产蜜桃级精品一区二区三区| 在线观看免费午夜福利视频| 国产黄片美女视频| 久久久久久久午夜电影| 久久午夜综合久久蜜桃| 美女大奶头视频| 午夜激情福利司机影院| 午夜两性在线视频| 亚洲精品国产精品久久久不卡| 丰满人妻熟妇乱又伦精品不卡| av天堂在线播放| 色精品久久人妻99蜜桃| 在线观看舔阴道视频| 免费在线观看完整版高清| 国产亚洲欧美98| 国产一区二区三区视频了| 两性夫妻黄色片| 成人特级黄色片久久久久久久| 日本 欧美在线| 国产午夜精品论理片| 久久久水蜜桃国产精品网| 亚洲av电影不卡..在线观看| 亚洲真实伦在线观看| 午夜精品在线福利| 色老头精品视频在线观看| 九九热线精品视视频播放| 久久热在线av| 国产又色又爽无遮挡免费看| 一进一出抽搐动态| 亚洲熟女毛片儿| 免费人成视频x8x8入口观看| 91在线观看av| АⅤ资源中文在线天堂| 色尼玛亚洲综合影院| 五月伊人婷婷丁香| 亚洲成a人片在线一区二区| 老汉色∧v一级毛片| 1024视频免费在线观看| 超碰成人久久| 九九热线精品视视频播放| 九色国产91popny在线| 久久亚洲精品不卡| 国产97色在线日韩免费| 久久精品人妻少妇| 亚洲成av人片免费观看| av福利片在线观看| 岛国在线免费视频观看| 国产亚洲精品综合一区在线观看 | 国产日本99.免费观看| 好男人电影高清在线观看| 成人手机av| 国产精品1区2区在线观看.| 国产黄色小视频在线观看| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 国产久久久一区二区三区| 亚洲成人中文字幕在线播放| 国产区一区二久久| 午夜福利在线在线| 国产成人aa在线观看| 在线观看免费视频日本深夜| 狠狠狠狠99中文字幕| 91九色精品人成在线观看| 国产午夜精品久久久久久| 最新美女视频免费是黄的| 久久精品国产亚洲av香蕉五月| 欧美成人免费av一区二区三区| 此物有八面人人有两片| 18禁黄网站禁片午夜丰满| 欧美日韩黄片免| 十八禁人妻一区二区| 亚洲av五月六月丁香网| 亚洲一码二码三码区别大吗| 搞女人的毛片| 在线观看美女被高潮喷水网站 | 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲一级av第二区| 国产精品 国内视频| 曰老女人黄片| 老鸭窝网址在线观看| 757午夜福利合集在线观看| 日本免费a在线| 欧美高清成人免费视频www| 精品久久久久久久毛片微露脸| 窝窝影院91人妻| 国产黄a三级三级三级人| 日本一本二区三区精品| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 丰满人妻一区二区三区视频av | 国产亚洲av高清不卡| 欧美黑人巨大hd| 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 又黄又粗又硬又大视频| 极品教师在线免费播放| 成人永久免费在线观看视频| 俄罗斯特黄特色一大片| 不卡av一区二区三区| 国产亚洲精品久久久久5区| 亚洲自偷自拍图片 自拍| 午夜福利免费观看在线| 国产亚洲精品综合一区在线观看 | 亚洲在线自拍视频| 国产精品免费一区二区三区在线| 日韩欧美国产在线观看| 亚洲,欧美精品.| 欧美一级毛片孕妇| ponron亚洲| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| or卡值多少钱| 成人一区二区视频在线观看| 变态另类丝袜制服| 精品高清国产在线一区| 九九热线精品视视频播放| 91九色精品人成在线观看| 叶爱在线成人免费视频播放| 日韩国内少妇激情av| 亚洲av熟女| 欧美精品亚洲一区二区| 这个男人来自地球电影免费观看| 国产精品亚洲一级av第二区| 久久九九热精品免费| 麻豆久久精品国产亚洲av| av福利片在线| 午夜免费成人在线视频| aaaaa片日本免费| 欧美在线黄色| 伦理电影免费视频| 99精品欧美一区二区三区四区| 高潮久久久久久久久久久不卡|