• <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)分析
    18禁国产床啪视频网站| 国产片内射在线| 国产成人欧美| 18禁裸乳无遮挡免费网站照片 | 午夜福利免费观看在线| 国产一卡二卡三卡精品| 午夜福利一区二区在线看| 色在线成人网| 人人妻人人澡欧美一区二区 | √禁漫天堂资源中文www| 韩国av一区二区三区四区| 一边摸一边抽搐一进一出视频| 国产视频一区二区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 色综合亚洲欧美另类图片| 国产黄a三级三级三级人| 少妇裸体淫交视频免费看高清 | 777久久人妻少妇嫩草av网站| 亚洲成人久久性| 中文字幕人成人乱码亚洲影| 亚洲av电影在线进入| www.999成人在线观看| 国产精品爽爽va在线观看网站 | 欧美日韩瑟瑟在线播放| 一区二区三区激情视频| 亚洲美女黄片视频| 深夜精品福利| 久久天躁狠狠躁夜夜2o2o| 桃色一区二区三区在线观看| 日日夜夜操网爽| 亚洲国产精品合色在线| 神马国产精品三级电影在线观看 | 在线观看免费视频日本深夜| 精品电影一区二区在线| 免费在线观看完整版高清| 男女之事视频高清在线观看| 日本一区二区免费在线视频| 两个人看的免费小视频| 免费观看精品视频网站| 色综合婷婷激情| 黑人巨大精品欧美一区二区mp4| 成人三级黄色视频| 国产亚洲欧美在线一区二区| 男人操女人黄网站| 午夜成年电影在线免费观看| 久久久久久久久中文| 亚洲精品国产区一区二| 人人妻,人人澡人人爽秒播| 黄色 视频免费看| 九色亚洲精品在线播放| 国产免费男女视频| 国产99白浆流出| 精品久久久久久,| 日韩欧美三级三区| 日韩视频一区二区在线观看| 国产一区二区激情短视频| 日本a在线网址| 嫁个100分男人电影在线观看| 黄色片一级片一级黄色片| 国产片内射在线| 国产1区2区3区精品| 少妇熟女aⅴ在线视频| 亚洲三区欧美一区| 欧美av亚洲av综合av国产av| 女人高潮潮喷娇喘18禁视频| 欧美黄色淫秽网站| 成人手机av| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美另类亚洲清纯唯美| 成人永久免费在线观看视频| 老司机靠b影院| 丰满的人妻完整版| 免费看美女性在线毛片视频| 大型av网站在线播放| 免费看美女性在线毛片视频| 十八禁人妻一区二区| 女人精品久久久久毛片| 亚洲情色 制服丝袜| 久久久久久久久中文| 亚洲精品一卡2卡三卡4卡5卡| 自线自在国产av| av欧美777| 亚洲熟妇中文字幕五十中出| 黑丝袜美女国产一区| 亚洲人成77777在线视频| 麻豆av在线久日| 人人妻人人澡人人看| 免费久久久久久久精品成人欧美视频| 看免费av毛片| 欧美激情久久久久久爽电影 | 国产精品av久久久久免费| 亚洲五月天丁香| 69精品国产乱码久久久| 首页视频小说图片口味搜索| 级片在线观看| 久久精品亚洲精品国产色婷小说| 亚洲精品在线美女| 亚洲av成人不卡在线观看播放网| 欧美在线黄色| 国产精品九九99| av网站免费在线观看视频| 国产精华一区二区三区| 日本撒尿小便嘘嘘汇集6| 法律面前人人平等表现在哪些方面| 麻豆一二三区av精品| 无限看片的www在线观看| 精品一区二区三区视频在线观看免费| 午夜精品国产一区二区电影| 日韩 欧美 亚洲 中文字幕| 欧美黄色片欧美黄色片| 免费高清视频大片| 国产亚洲精品久久久久5区| 男女之事视频高清在线观看| tocl精华| 欧美一级毛片孕妇| 一级毛片精品| 变态另类成人亚洲欧美熟女 | 俄罗斯特黄特色一大片| 日日干狠狠操夜夜爽| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产一区二区精华液| 午夜免费鲁丝| 国语自产精品视频在线第100页| 国产精品综合久久久久久久免费 | 高潮久久久久久久久久久不卡| 91老司机精品| 欧美性长视频在线观看| 欧美日韩黄片免| 男人舔女人下体高潮全视频| 免费看美女性在线毛片视频| 国产精品,欧美在线| 国产欧美日韩一区二区三| 宅男免费午夜| 日日干狠狠操夜夜爽| 人人妻人人澡人人看| 国产成人欧美在线观看| 99精品在免费线老司机午夜| 在线十欧美十亚洲十日本专区| 国产视频一区二区在线看| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品999在线| 国产私拍福利视频在线观看| 欧美日本中文国产一区发布| 午夜福利免费观看在线| 国产成人精品久久二区二区91| 成年人黄色毛片网站| 亚洲第一电影网av| 丝袜美足系列| 一本大道久久a久久精品| 老熟妇乱子伦视频在线观看| 少妇粗大呻吟视频| 两个人免费观看高清视频| 99国产精品免费福利视频| 丝袜美腿诱惑在线| 欧美色视频一区免费| 久久伊人香网站| 色综合婷婷激情| 色哟哟哟哟哟哟| 久久天堂一区二区三区四区| 香蕉丝袜av| 美女 人体艺术 gogo| 国产在线精品亚洲第一网站| 日韩精品中文字幕看吧| 国产99白浆流出| 日韩 欧美 亚洲 中文字幕| 日韩欧美免费精品| 久99久视频精品免费| 窝窝影院91人妻| 激情视频va一区二区三区| 午夜免费观看网址| 十八禁人妻一区二区| 亚洲情色 制服丝袜| 欧美日韩瑟瑟在线播放| 精品人妻1区二区| 日韩三级视频一区二区三区| 波多野结衣一区麻豆| 免费在线观看日本一区| 国产精品免费视频内射| 午夜a级毛片| 精品久久久久久久毛片微露脸| 一本久久中文字幕| 老司机福利观看| 好看av亚洲va欧美ⅴa在| 97碰自拍视频| 久久久水蜜桃国产精品网| 18禁裸乳无遮挡免费网站照片 | 国产精品久久视频播放| 大码成人一级视频| svipshipincom国产片| 久久精品亚洲熟妇少妇任你| 久久香蕉精品热| 久久精品国产清高在天天线| 满18在线观看网站| 十分钟在线观看高清视频www| 咕卡用的链子| 久久欧美精品欧美久久欧美| 国产精品二区激情视频| 亚洲五月婷婷丁香| 香蕉国产在线看| 精品第一国产精品| 少妇 在线观看| 两个人免费观看高清视频| 国产伦人伦偷精品视频| av片东京热男人的天堂| 真人做人爱边吃奶动态| 国产一级毛片七仙女欲春2 | 亚洲av成人av| 男人的好看免费观看在线视频 | 老熟妇仑乱视频hdxx| 国产成人系列免费观看| 叶爱在线成人免费视频播放| 女人精品久久久久毛片| 成人永久免费在线观看视频| 亚洲色图av天堂| 两个人免费观看高清视频| 中亚洲国语对白在线视频| 色尼玛亚洲综合影院| 亚洲电影在线观看av| 88av欧美| 真人一进一出gif抽搐免费| 国产激情欧美一区二区| 大型黄色视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久亚洲精品不卡| 久久精品国产综合久久久| 欧美一区二区精品小视频在线| 日韩三级视频一区二区三区| 亚洲人成电影免费在线| 中文字幕av电影在线播放| 乱人伦中国视频| 亚洲无线在线观看| 日本vs欧美在线观看视频| 午夜日韩欧美国产| 人人妻人人澡欧美一区二区 | 亚洲欧洲精品一区二区精品久久久| 国语自产精品视频在线第100页| 黄色女人牲交| 在线观看午夜福利视频| 正在播放国产对白刺激| 亚洲三区欧美一区| 亚洲激情在线av| 日韩大码丰满熟妇| 国产av一区在线观看免费| 在线观看66精品国产| 亚洲,欧美精品.| 18禁裸乳无遮挡免费网站照片 | 精品欧美一区二区三区在线| 天天一区二区日本电影三级 | 色尼玛亚洲综合影院| 日韩中文字幕欧美一区二区| 国产私拍福利视频在线观看| 黄片大片在线免费观看| 国产精品电影一区二区三区| 悠悠久久av| 欧美日韩中文字幕国产精品一区二区三区 | 日韩一卡2卡3卡4卡2021年| 国产熟女xx| 国产不卡一卡二| 啪啪无遮挡十八禁网站| 久久久久国产一级毛片高清牌| 色婷婷久久久亚洲欧美| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 久久影院123| 国产欧美日韩一区二区三区在线| 俄罗斯特黄特色一大片| 成人三级做爰电影| 变态另类成人亚洲欧美熟女 | 国产1区2区3区精品| 亚洲自拍偷在线| 男女之事视频高清在线观看| 国产激情久久老熟女| 国产精品1区2区在线观看.| 国产精品免费一区二区三区在线| 嫁个100分男人电影在线观看| 日日爽夜夜爽网站| 国产男靠女视频免费网站| 久久久久久久久久久久大奶| 国产精品一区二区在线不卡| 欧美亚洲日本最大视频资源| 午夜福利免费观看在线| 亚洲第一欧美日韩一区二区三区| 又紧又爽又黄一区二区| 免费高清视频大片| 国产97色在线日韩免费| 免费无遮挡裸体视频| 天堂√8在线中文| 国产麻豆成人av免费视频| 精品久久久精品久久久| av在线播放免费不卡| 国产av又大| 亚洲国产日韩欧美精品在线观看 | 成人特级黄色片久久久久久久| 久久久久久大精品| 女人被狂操c到高潮| 变态另类成人亚洲欧美熟女 | 窝窝影院91人妻| 天天躁夜夜躁狠狠躁躁| 亚洲精品美女久久久久99蜜臀| 夜夜夜夜夜久久久久| 少妇 在线观看| 成人亚洲精品一区在线观看| 中出人妻视频一区二区| 少妇 在线观看| 黑丝袜美女国产一区| 色综合欧美亚洲国产小说| 中文亚洲av片在线观看爽| 黑丝袜美女国产一区| 久久午夜综合久久蜜桃| 女人高潮潮喷娇喘18禁视频| 又大又爽又粗| 久久国产亚洲av麻豆专区| 亚洲美女黄片视频| 在线观看午夜福利视频| 麻豆一二三区av精品| 国产真人三级小视频在线观看| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看 | 美女大奶头视频| 欧美日韩亚洲综合一区二区三区_| 91麻豆av在线| 国产午夜福利久久久久久| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 国产一区二区三区视频了| 国产伦人伦偷精品视频| 中文字幕最新亚洲高清| 久久欧美精品欧美久久欧美| 1024视频免费在线观看| 亚洲欧美一区二区三区黑人| 在线播放国产精品三级| 老司机靠b影院| 日本免费a在线| 国产日韩一区二区三区精品不卡| 午夜福利高清视频| 亚洲aⅴ乱码一区二区在线播放 | 日韩视频一区二区在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲成人免费电影在线观看| 免费高清在线观看日韩| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 亚洲精华国产精华精| 国产激情欧美一区二区| 欧美在线黄色| 久久久国产成人免费| 99国产精品一区二区三区| 91精品三级在线观看| 精品久久蜜臀av无| 丝袜在线中文字幕| 两个人视频免费观看高清| 91大片在线观看| 久久久久亚洲av毛片大全| cao死你这个sao货| 国产1区2区3区精品| 久久狼人影院| 午夜精品在线福利| 亚洲国产精品合色在线| 亚洲中文av在线| 久久人妻av系列| 亚洲国产精品999在线| 香蕉国产在线看| 曰老女人黄片| 亚洲成av片中文字幕在线观看| 国产成人欧美在线观看| 妹子高潮喷水视频| 午夜亚洲福利在线播放| 两个人看的免费小视频| 午夜福利18| 桃红色精品国产亚洲av| 国产熟女xx| 亚洲男人天堂网一区| 欧美一级毛片孕妇| 女人被躁到高潮嗷嗷叫费观| 亚洲七黄色美女视频| 我的亚洲天堂| 成人三级黄色视频| 久久久国产成人免费| 精品一区二区三区av网在线观看| 涩涩av久久男人的天堂| 亚洲精品久久国产高清桃花| 99久久99久久久精品蜜桃| 免费av毛片视频| 99久久精品国产亚洲精品| 成人特级黄色片久久久久久久| 黄频高清免费视频| 国产99久久九九免费精品| 母亲3免费完整高清在线观看| 精品久久久久久久人妻蜜臀av | 色播在线永久视频| 一区二区三区高清视频在线| 国产精品爽爽va在线观看网站 | 1024视频免费在线观看| 一区二区三区国产精品乱码| 搡老岳熟女国产| 一本久久中文字幕| 久久热在线av| 久久人人精品亚洲av| 免费看十八禁软件| 99精品在免费线老司机午夜| 国产精品久久久av美女十八| 如日韩欧美国产精品一区二区三区| 亚洲国产日韩欧美精品在线观看 | 搞女人的毛片| 在线免费观看的www视频| 免费不卡黄色视频| 亚洲电影在线观看av| 亚洲熟妇中文字幕五十中出| 精品久久久久久成人av| av天堂在线播放| 一边摸一边抽搐一进一小说| 韩国av一区二区三区四区| 国产精品久久久久久精品电影 | 十分钟在线观看高清视频www| 国产精品亚洲一级av第二区| 精品久久久久久久毛片微露脸| 久久久久精品国产欧美久久久| 此物有八面人人有两片| 中亚洲国语对白在线视频| 搞女人的毛片| av有码第一页| 伦理电影免费视频| 波多野结衣高清无吗| 午夜免费激情av| 韩国av一区二区三区四区| 天堂动漫精品| svipshipincom国产片| 久久久久亚洲av毛片大全| 一边摸一边抽搐一进一小说| 波多野结衣巨乳人妻| 成人三级黄色视频| 少妇熟女aⅴ在线视频| 免费观看人在逋| 在线播放国产精品三级| 最好的美女福利视频网| 欧美色视频一区免费| 精品不卡国产一区二区三区| 国产精品自产拍在线观看55亚洲| 我的亚洲天堂| 亚洲第一av免费看| 美女 人体艺术 gogo| 亚洲五月色婷婷综合| 亚洲九九香蕉| 精品国内亚洲2022精品成人| 在线永久观看黄色视频| 欧美一级毛片孕妇| 亚洲伊人色综图| 香蕉国产在线看| 国产成+人综合+亚洲专区| 日韩欧美三级三区| svipshipincom国产片| 人人妻,人人澡人人爽秒播| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 纯流量卡能插随身wifi吗| 中国美女看黄片| 亚洲一区二区三区不卡视频| 天堂√8在线中文| 欧美中文综合在线视频| 亚洲精品一区av在线观看| 午夜精品久久久久久毛片777| 国产精品1区2区在线观看.| 亚洲九九香蕉| 日韩大尺度精品在线看网址 | 欧美黄色淫秽网站| 久久久久久人人人人人| 1024视频免费在线观看| 无限看片的www在线观看| 亚洲中文av在线| 最近最新中文字幕大全电影3 | 亚洲国产高清在线一区二区三 | 精品乱码久久久久久99久播| 老司机福利观看| 97碰自拍视频| 波多野结衣一区麻豆| 一级a爱视频在线免费观看| 日韩欧美免费精品| 一区二区三区国产精品乱码| 欧美激情 高清一区二区三区| 国产精品久久电影中文字幕| 久久久国产欧美日韩av| 在线观看免费午夜福利视频| 精品一区二区三区视频在线观看免费| 母亲3免费完整高清在线观看| 超碰成人久久| 欧美激情高清一区二区三区| 一区二区三区国产精品乱码| 国产精品一区二区三区四区久久 | 精品熟女少妇八av免费久了| 亚洲三区欧美一区| 国产精品久久久久久人妻精品电影| 亚洲午夜理论影院| 啦啦啦 在线观看视频| 亚洲欧美日韩高清在线视频| 丰满的人妻完整版| 欧美中文综合在线视频| 国产一区在线观看成人免费| av片东京热男人的天堂| 亚洲全国av大片| 国产欧美日韩综合在线一区二区| 亚洲中文字幕一区二区三区有码在线看 | av视频免费观看在线观看| 精品午夜福利视频在线观看一区| 亚洲av成人不卡在线观看播放网| 一区在线观看完整版| 亚洲人成网站在线播放欧美日韩| 久久久久久久久久久久大奶| 色婷婷久久久亚洲欧美| 91av网站免费观看| 韩国av一区二区三区四区| 亚洲人成77777在线视频| 精品高清国产在线一区| 欧美日本中文国产一区发布| 9色porny在线观看| 免费av毛片视频| 久久伊人香网站| 国产又色又爽无遮挡免费看| 国产精品久久电影中文字幕| 亚洲第一av免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产又爽黄色视频| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 国产97色在线日韩免费| 人人妻,人人澡人人爽秒播| 极品人妻少妇av视频| 亚洲av成人不卡在线观看播放网| 免费在线观看黄色视频的| 嫩草影视91久久| 日本vs欧美在线观看视频| 亚洲欧美日韩无卡精品| 精品久久蜜臀av无| 亚洲片人在线观看| 国产精品一区二区在线不卡| 国产精品免费一区二区三区在线| av天堂久久9| 深夜精品福利| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 黄色毛片三级朝国网站| 亚洲中文av在线| 久久久久久久久中文| 精品国产乱码久久久久久男人| 欧美国产精品va在线观看不卡| 91av网站免费观看| 757午夜福利合集在线观看| 亚洲精华国产精华精| 亚洲欧美激情在线| 欧美色视频一区免费| 亚洲专区中文字幕在线| 亚洲一区高清亚洲精品| 免费少妇av软件| 国内久久婷婷六月综合欲色啪| 他把我摸到了高潮在线观看| 91大片在线观看| 欧美日本中文国产一区发布| 真人做人爱边吃奶动态| 丝袜美腿诱惑在线| 在线十欧美十亚洲十日本专区| 国产精品久久久久久人妻精品电影| 一二三四在线观看免费中文在| 一区二区三区精品91| 黄色 视频免费看| 美女 人体艺术 gogo| 亚洲国产毛片av蜜桃av| 欧美一级毛片孕妇| 国产欧美日韩一区二区三| 最新美女视频免费是黄的| 变态另类丝袜制服| 一本综合久久免费| 免费人成视频x8x8入口观看| 麻豆成人av在线观看| 别揉我奶头~嗯~啊~动态视频| www.自偷自拍.com| 久久精品国产99精品国产亚洲性色 | 久久国产精品男人的天堂亚洲| 丰满人妻熟妇乱又伦精品不卡| 一个人免费在线观看的高清视频| 精品久久久久久久人妻蜜臀av | 国产激情久久老熟女| 久久中文看片网| 日日爽夜夜爽网站| av在线天堂中文字幕| 老司机深夜福利视频在线观看| 电影成人av| 一区二区三区激情视频| 欧美色视频一区免费| 国产av在哪里看| 成熟少妇高潮喷水视频| 久久精品国产清高在天天线| 人妻丰满熟妇av一区二区三区| 久久久精品欧美日韩精品| 女性被躁到高潮视频| 视频在线观看一区二区三区| 久久精品91蜜桃| 麻豆av在线久日| 国产99白浆流出| 国产精品美女特级片免费视频播放器 | 999久久久国产精品视频| 久久久精品国产亚洲av高清涩受| 国产xxxxx性猛交| av天堂久久9| 久久久精品欧美日韩精品| 日韩精品中文字幕看吧| 在线观看免费日韩欧美大片| 久久欧美精品欧美久久欧美| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 丝袜美足系列| 99国产精品一区二区蜜桃av| 侵犯人妻中文字幕一二三四区| 久久午夜亚洲精品久久| www.自偷自拍.com| 亚洲男人的天堂狠狠| 老汉色av国产亚洲站长工具|