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

    三江平原濕地植被—水位轉(zhuǎn)換函數(shù)研究

    2017-10-10 02:03:25王春玲郭晉平王國平婁彥景
    關(guān)鍵詞:樣方植被物種

    王春玲,郭晉平,王國平,婁彥景*

    (1.山西農(nóng)業(yè)大學 城鄉(xiāng)建設學院,山西 太谷 030801;2.中國科學院濕地生態(tài)與環(huán)境重點實驗室,中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所,吉林 長春 130102)

    三江平原濕地植被—水位轉(zhuǎn)換函數(shù)研究

    王春玲1,2,郭晉平1,王國平2,婁彥景2*

    (1.山西農(nóng)業(yè)大學 城鄉(xiāng)建設學院,山西 太谷 030801;2.中國科學院濕地生態(tài)與環(huán)境重點實驗室,中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所,吉林 長春 130102)

    [目的]植物殘體是沼澤濕地古生態(tài)學研究的重要生物指標,現(xiàn)代植被與環(huán)境因子的轉(zhuǎn)換函數(shù)為準確定量重建古環(huán)境和古氣候提供了有效途徑。[方法]本文利用R軟件“rioja”軟件包中的加權(quán)平均(Weighted Averaging,WA)和加權(quán)平均偏最小二乘法(Weighted Averaging Partial Least Squares,WAPLS)模型,以三江平原不同類型植物群落數(shù)據(jù)與相應的水位數(shù)據(jù)為樣本,建立現(xiàn)代植被-水位轉(zhuǎn)換函數(shù)。[結(jié)果]WAPLS模型的預測性能比WA模型好,WAPLS模型第3組分的預測性能最好,有最低的RMSEP值(3.69 m)和最高的R2值(0.52)。經(jīng)數(shù)據(jù)過濾后,兩個模型的預測性能有所提高,其中WAPLS模型第4組分的預測性能最好,RMSEP為2.83,降低了23.3 %,其精度達到了±2.83 cm,結(jié)果較理想。[結(jié)論]研究結(jié)果可為三江平原古水位定量重建提供方法支撐,提高其準確性。

    除趨勢對應分析(DCA);WA模型;WAPLS模型;預測能力

    Abstract:[Objective]Plant macrofossil is an important biological proxy for palaeoecology study. The transfer function of vegetation-water level can provide an effective way to reconstruct accurate quantitative palaeoenvironment and palaeoclimate. [Methods] In this paper, the WA(Weighted Averaging) and WAPLS(Weighted Averaging Partial Least Squares) models ofriojapackage in R were used to establish the vegetation-water level transfer function of Sanjiang Plain.[Results]The results showed that the prediction performance of the WAPLS model was better than that of the WA model. The third component of the WAPLS model had the best prediction performance, with the lowest RMSEP value (3.69 m) and the highest R2value (0.52). After data filtering, the prediction performance of the two models were improved, the prediction performance of the fourth component of the WAPLS model was the best, with the lowest RMSEP 2.83, which reduced by 23.3% compared to the previous value. The accuracy reached to ±2.83 cm, this result was relatively ideal.[Conclusion]This study provided a method and guarantee for the accurate quantitative reconstruction of the water level in Sanjiang Plain.

    Keywords:Detrended Correspondence Analysis(DCA), WA model, WAPLS model, Predication performance

    全球氣候變化及其對生態(tài)系統(tǒng)的影響是目前國內(nèi)外科學研究的熱點領域[1],定量重建古環(huán)境和古氣候是其中非常重要的組成部分。泥炭是記錄古植被、古環(huán)境和古氣候的天然地質(zhì)檔案,在古生態(tài)學中的地位越來越重要[2~8]。泥炭是沼澤濕地的特有產(chǎn)物,是以生物殘體為主的沉積物,植物殘體是其中最主要的組成成分,因此生物指標是泥炭古環(huán)境研究中最重要的代用指標[8~11]。定量重建古環(huán)境的重要途徑是沼澤濕地現(xiàn)代生物群落與環(huán)境因子間建立定量關(guān)系,目前研究人員普遍采用轉(zhuǎn)換函數(shù)的方法[12~21]。轉(zhuǎn)換函數(shù)在孢粉[13,14,16,20,22~29]、硅藻[21,30~32]、有孔蟲[17,33]、搖蚊[34~36]等生物指標的定量重建古環(huán)境研究中得到了廣泛的應用。沼澤濕地植物對環(huán)境、氣候變化響應敏感,很好的記錄了環(huán)境的變化過程[37~39],研究表明水位是沼澤濕地植物群落多樣性及分布格局的主要驅(qū)動因素[38,39]。泥炭地植物殘體是古生態(tài)學研究的一個非常重要的生物指標,主要用來重建古植被的演替過程[8~11,40,41]和古環(huán)境演化序列[42~46]。因此建立沼澤濕地現(xiàn)代植被與水位的轉(zhuǎn)換函數(shù)是定量重建沼澤濕地古水位的關(guān)鍵。國外已有一些學者利用植被—水位轉(zhuǎn)換函數(shù)定量重建古水位[28,47],而國內(nèi)相關(guān)研究僅僅是利用泥炭地植物殘體指標定性重建了泥炭地全新世以來的古環(huán)境[42,44~46]。本文利用三江平原典型沼澤濕地植被與水位間的關(guān)系建立轉(zhuǎn)換函數(shù),為進一步開展定量重建古水位研究奠定基礎。

    1 研究方法

    1.1 野外植被與水位調(diào)查

    2015年7-8月,在三江平原洪河保護區(qū)、試驗場、洪河塔、興凱湖等濕地進行植物樣方調(diào)查,同時記錄每個樣方所在位置的水位(cm)。共調(diào)查107個植物樣方,包含植物物種53種。每個植物樣方記錄每種植物的高度(cm)和蓋度(%)。調(diào)查的植物群落主要有毛果苔草群落,烏拉苔草群落,灰脈苔草群落,小葉章群落,甜茅群落,毛果—漂筏苔草群落,小葉章—烏拉苔草群落,這些群落是三江平原的典型濕地植物群落。

    1.2 轉(zhuǎn)換函數(shù)的建立

    本研究首先利用CANOCO4.5軟件對植物群落數(shù)據(jù)進行DCA排序分析,探討影響不同類型植物群落分布的主要環(huán)境因子。本文轉(zhuǎn)換函數(shù)的建立所需要的樣本數(shù)據(jù)集為物種的蓋度及相應的水位數(shù)據(jù)。本研究利用R軟件里的rioja軟件包的加權(quán)平均模型(Weighted Averaging, WA)和加權(quán)平均偏最小二乘法模型(Weighted Averaging-Partial Least Square,WA-PLS)模型進行植被—水位轉(zhuǎn)換函數(shù)的建立[14]。轉(zhuǎn)換函數(shù)的推導能力采用刀割法(leave-one-out或jack-knifing)進行交叉驗證,該驗證方法所用的測試集為每個單獨樣本,訓練集為剩余的n-1個樣本,以此來建立模型[33]。轉(zhuǎn)換函數(shù)的預測性能主要通過RMSEP(Root Mean Square Error of Prediction)和r2(squared correlation between inferred and observed values)兩個指標進行評價。RMSEP為預測值與實際觀測值之間的誤差,R2為預測值和觀測值的相關(guān)關(guān)系,RMSEP越小,R2越大,則模型的預測性能越好。選取目標函數(shù)的依據(jù)為:R2最大而RMSEP最小的函數(shù)。如果模型的預測模擬結(jié)果不理想,其可能的原因為殘差值過大的異常樣品使轉(zhuǎn)換函數(shù)模型的預測性能降低。因此如果想提高模型的預測能力,則需要將殘差值過大的樣本去掉,目前國際上常將殘差值超過樣本梯度20%的樣本剔除以優(yōu)化模型[14]。

    調(diào)查的植物物種共53種,選擇至少在6個樣方中出現(xiàn),且蓋度含量超過3%的相對優(yōu)勢物種作為建立轉(zhuǎn)換函數(shù)的物種,經(jīng)篩選共有21個物種可用于本研究中物種—水位轉(zhuǎn)換函數(shù)的建立。

    2 結(jié)果與分析

    2.1 DCA分析結(jié)果

    對調(diào)查的107個樣方進行DCA排序分析,由表1可見,第1軸的特征值明顯高于第2軸,表明植物的分布主要由第一軸的環(huán)境因子揭示;第一軸對物種的解釋量最大,為22.5%,前4軸可以解釋物種累積變異的48.7%。以第一軸和第二軸為坐標做群落排序軸的二維散點圖,如圖1和圖2所示。

    表1植被樣方DCA排序分析結(jié)果
    Table1 The results of DCA analysis of vegetation communities

    排序軸Axes1234特征值0.7210.5790.1610.099梯度長度3.0743.3812.3801.641累計物種變化百分比/%22.540.645.648.7總特征值3.206

    圖1a為植被調(diào)查物種的DCA排序圖,從圖中看出,物種排列方式揭示了一定的環(huán)境梯度。從DCA排序第1軸發(fā)現(xiàn),指示長期積水環(huán)境的物種如毛果苔草(Carexlasiocarpa,C.lasio)、睡菜(Menyanthestrifoliate,Menytrifo)、沼委陵菜(Comarumpalustre,Compalus)等位于最右端,而繡線菊(Spiraeasalicifolia,Spirsalic)、小葉章(Calamagrostisangustifolia,Calaangus)、地榆(Sanguisorbaofficinalis,Sangoffic)、沼柳(Salixrosmarinifoliavar.brachypoda,Salirosma)等喜干的植物則位于排序軸的左端,由此認為第一排序軸從左到右物種呈現(xiàn)干濕分布,即從左到右水分梯度呈現(xiàn)由干到濕的變化過程。物種DCA第2排序軸則相對比較復雜,沒有明確的環(huán)境梯度指示意義。從圖中還可以看出,每種植物在排序圖中有各自的分布中心和區(qū)域,這由每種植物都有各自所適宜的生長環(huán)境和氣候條件所決定。

    圖1 DCA排序圖,a為物種排序圖,b為樣方排序圖Fig.1 DCA ordination diagram of species (a) and samples (b)

    圖1b為樣方DCA排序圖,從圖中可以看出,樣方在排序圖上的分布有明顯的分區(qū)效應,相同區(qū)域內(nèi)的樣品有明顯的重疊現(xiàn)象,表明這些樣方具有相似的物種組成。相同區(qū)域內(nèi)的樣方群落具有相似的結(jié)構(gòu)組成和相似的生境條件。反之,在排序軸上距離越遠的樣方則它們的物種組成和生境條件等差異性越大。

    2.2 轉(zhuǎn)換函數(shù)建立

    首先使用全部樣方數(shù)據(jù)進行轉(zhuǎn)換函數(shù)模型的構(gòu)建。WA和WAPLS模型的預測性能見表2。從表中可以看出,WA.inv.tol模型的RMSE=3.94 cm,R2=0.47,相對于WA其它模型來說,RESM最小,R2最大,預測性能較好。WAPLS模型中,WAPLS模型組分3的RMSE最小,為3.69 cm,R2最大,為0.52,預測性能較其它組分要好。通過對WA和WAPLS兩模型相比較,發(fā)現(xiàn)WAPLS模型有較小的RMSE,較大的R2,因此WAPLS模型第3組分的預測性能最好,WAPLS模型是建立植被—水位轉(zhuǎn)換函數(shù)的最佳選擇。 很多研究證明,影響植被植被—水位轉(zhuǎn)換函數(shù)模型的預測性能最主要的因素為殘差值較大的樣品,因此為了為降低預測誤差,提高轉(zhuǎn)換函數(shù)模型的預測性能,通常需要從訓練樣本中將殘差值大于環(huán)境梯度寬度20%的樣品剔除掉。經(jīng)過這樣的數(shù)據(jù)處理, 8個樣品從訓練樣本中濾除,剩余99個樣品進行轉(zhuǎn)換函數(shù)預測性能分析,其結(jié)果見表3。

    表2 植被—水位轉(zhuǎn)換函數(shù)預測性能比較(全部樣方)Table 2 Predication performance of transfer functions for water level (all samples)

    表3 植被—水位轉(zhuǎn)換函數(shù)預測性能比較(優(yōu)化后)Table 3 Predication performance of transfer functions for water level (after data filtering)

    從表3中可以看出,經(jīng)過優(yōu)化后WA和WALPS模型的預測性能都有明顯的提高,RMSE都有所減少,R2都有所增加,所有的模型中,WAPLS模型第4分量的RMSEP 最小,為2.83,比全樣降低了23.3%;R2也最大,為0.70,比全樣提高了34.6%;另一個判斷模型預測性能的指標最大偏差也從原來的8.85降到3.35。從所有的評價模型預測性能的指標來看,WAPLS模型第4分量的預測性能最好,可以用來定量重建三江平原沼澤濕地古水位。

    通過繪制轉(zhuǎn)換函數(shù)模型的實測值與預測值散點圖、實測值與殘差值散點圖,可以直觀地展示模型的預測性能。圖2為水位—植被轉(zhuǎn)換函數(shù)預測性能圖,從下列散點圖中可以發(fā)現(xiàn),圖a和圖e中大部分樣點都比較集中,只有部分樣點比較偏離1∶1線,圖b圖f中大部分樣點都集中于0軸附近,少部分樣點距離0軸較遠,說明轉(zhuǎn)換函數(shù)模型預測性能相對較差。圖c,d,g,h為優(yōu)化后的預測性能圖,圖c和g中的樣點更集中線于1∶1,殘差值較之前變小,圖d和h中樣點更集中于0軸附近。這說明優(yōu)化后的模型具有更好的預測性能。

    3 討論

    定量重建古環(huán)境、古氣候是目前全球氣候變化的重要研究內(nèi)容之一。濕地植被對氣候、環(huán)境變化相應敏感,因此可以通過重建古植被來反演全新世以來的環(huán)境變化信息[37~39]。目前國外利用泥炭地植物殘體定性重建古環(huán)境的研究較多,而定量研究則相對較少,我國研究學者對植物殘體指標的研究相對比較薄弱[28,47]。以往的研究發(fā)現(xiàn),影響三江平原濕地植物群落分布的主要環(huán)境因素是水位[38,39],因此本研究直接探討水位與植物群落之間的關(guān)系,進而建立轉(zhuǎn)換函數(shù),探討其預測性能。

    轉(zhuǎn)換函數(shù)的建立首先要考慮的是模型的選擇,目前用來建立現(xiàn)代生物與環(huán)境因子之間關(guān)系的轉(zhuǎn)換函數(shù)模型有WA,PLS,WAPLS,ML,MAT,LWWA等[15]。從國內(nèi)外研究來看,WA和WAPLS模型是最常用的兩個模型,轉(zhuǎn)換函數(shù)的預測性能主要根據(jù)RMSEP和R2來判斷。通過分析國內(nèi)外目前相關(guān)研究現(xiàn)狀,本研究采用WA和WAPLS模型建立三江平原濕地植被—水位轉(zhuǎn)換函數(shù),從中選擇有較低RMSEP值和較高R2值的模型。

    影響轉(zhuǎn)換函數(shù)模型預測性能的另一個主要因素為訓練樣本集數(shù)據(jù)質(zhì)量。首先,需要將所調(diào)查的生物物種進行篩選,將調(diào)查中出現(xiàn)的稀有種及含量很少物種去除,然后進行轉(zhuǎn)換函數(shù)的建立。對于不同的生物指標,不同的學者篩選建立轉(zhuǎn)換函數(shù)物種的方法并不完全相同,如有的將未在所有樣本中出現(xiàn)或在所有樣本中出現(xiàn)但含量小于3%的物種過濾掉;有的選擇至少在3個樣品中超過2%的相對優(yōu)勢屬種作為研究對象[17]。本研究中,樣方調(diào)查的植物物種共53種,選擇至少在6個樣方中出現(xiàn),且其蓋度含量超過3%的相對優(yōu)勢物種作為建立轉(zhuǎn)換函數(shù)的物種,經(jīng)過篩選后共有21個物種用于建立物種—水位轉(zhuǎn)換函數(shù)。

    圖2 WA(a,b,c,d)和WAPLS(e,f,g,h)模型預測性能(a,b,e,f優(yōu)化前,c,d,g,h為優(yōu)化后)Fig.2 Performance of transfer function models for water level for WA(a,b,c,d)and WAPLS(e,f,g,h), pre-(a,b,e,f) and after-data filtering(c,d,g,h)

    其次,轉(zhuǎn)換函數(shù)建立以后,通常為了提高模型的預測性能,需要將殘差值較大的異常樣本過濾掉,然后重新構(gòu)建轉(zhuǎn)換函數(shù)。對于不同的生物指標,不同的訓練樣本集,研究學者過濾樣本的標準不一樣,目前國內(nèi)外普遍采納的標準是將殘差值大于環(huán)境梯度寬度20%的樣本過濾掉[14],本研究也采用此標準重新建立轉(zhuǎn)換函數(shù)。通過分析發(fā)現(xiàn),WAPLS模型比WA模型的預測性能更好,WAPLS模型的第4組分預測性能最好,有最小RMSEP(2.83),比全樣(3.69)降低了23.3%;最大的R2(0.70),比全樣(0.52)提高了34.6%,因此可以利用此模型來定量重建三江平原全新世以來的古水位。

    雖然通過數(shù)據(jù)的過濾處理可以提高模型的預測性能,但是這種數(shù)據(jù)過濾方法也存在一定的弊端,它使訓練樣本的環(huán)境梯度寬度縮小了,通常被過濾掉的樣品都是處于環(huán)境梯度寬度的兩端,這樣影響定量重建的結(jié)果。本研究中經(jīng)過數(shù)據(jù)過濾后,水位寬度由原來的25 cm變成20 cm,降幅較大。所以在過濾數(shù)據(jù)的時候需要慎重考慮,用過濾后的轉(zhuǎn)換函數(shù)模型定量重建古環(huán)境值會出現(xiàn)偏差。如果三江平原濕地泥炭柱芯中植物殘體組成與本研究訓練樣本集植物群落組成相似,水位可以用±2.83 cm的平均誤差進行重建。反之,如果組成有較大差異,則不能利用本文的研究結(jié)果來定量重建古水位。

    本研究調(diào)查的植物群落是三江平原典型的植物群落,如小葉章群落,毛果苔草群落,甜茅群落,烏拉苔草群落等,水位梯度范圍為25 cm,且只有一年的調(diào)查數(shù)據(jù),不能代表三江平原地區(qū)所有的植物群落分布情況,所以本研究建立的植被—水位轉(zhuǎn)換函數(shù)還不能運用到本地區(qū)的所有泥炭沉積柱芯中。因此要想建立一個適合于三江平原全部泥炭柱芯的植被—水位轉(zhuǎn)換函數(shù),就需要在目前植被群落調(diào)查的基礎上進行大量的數(shù)據(jù)補充,盡可能涵蓋所有的植物群落類型,進而可以使柱芯植物殘體組成和現(xiàn)代植物組成有良好的對比,兩者的相似度越大,則通過轉(zhuǎn)換函數(shù)定量重建的古水位越精確。同時可以通過定量重建的百年尺度的水位變化與現(xiàn)代觀測值進行對比分析,來檢驗轉(zhuǎn)換函數(shù)的預測性能,并進一步提高轉(zhuǎn)換函數(shù)的推導能力。目前國內(nèi)對濕地植物群落與環(huán)境因子轉(zhuǎn)換函數(shù)的研究幾乎空白,因此本文具有很重要的研究意義,不僅可以促進定量重建全新世以來三江平原濕地高分辨率的古水位研究,同時為濕地的管理保護及退化濕地的恢復重建工作提供科學依據(jù)。

    4 結(jié)論

    本研究采用WA和WAPLS模型構(gòu)建了三江平原濕地植被—水位轉(zhuǎn)換函數(shù),研究表明WAPLS模型的預測性能比WA模型好,WAPLS模型第3組分的預測性能最好,有最低的RMSEP(3.69)值和最高的R2(0.52)值。為了提高模型的預測性能,采取了數(shù)據(jù)過濾方法,通過數(shù)據(jù)過濾后的兩個模型預測性能都有所提高,其中WAPLS模型的第4組分的預測性能最好,其精度達到了±2.83 cm,結(jié)果較理想。

    [1]王紹武, 葛全勝, 王芳, 等. 全球氣候變暖爭議中的核心問題[J]. 地球科學進展, 2010, 25(6): 656-665.

    [2]Barber KE. Peatlands as scientific archives of past biodiversity[J]. Biodiversity and conservation, 1993, 2(5): 474-489.

    [3]Barber KE, Charman DJ. Holocene palaeoclimate records from peatlands[M]. Global change in the Holocene. Arnold, Hodder Headline Group, 2003: 210-226.

    [4]Blackford J. Palaeoclimatic records from peat bogs[J]. Trends in Ecology and Evolution, 2000, 15(5): 193-198.

    [5]Chambers FM, Charman DJ. Holocene environmental change: contributions from the peatland archive[J]. The Holocene, 2004, 14(1): 1-6.

    [6]Charman D. Peatlands and environmental change[M]. Chichester:John Wiley and Sons, 2002:3-23.

    [7]Jong RD, Blaauw M, Chambers FM, et al. Climate and peatlands[M]. Changing Climates, Earth Systems and Society. Netherlands:Springer Science Business Media, 2010: 85-121.

    [8]Mauquoy D, Van Geel B. Mire and Peat Macros[M]. Encyclopedia of Quaternary Science. Poland: Elsevier, 2007: 2315-2336.

    [9]Chambers FM, Booth RK, Vleeschouwer FD, et al. Development and refinement of proxy-climate indicators from peats[J]. Quaternary International, 2012, 268(9):21-33.

    [10]Mauquoy D, Hughes PDM, Van Geel B. A protocol for plant macrofossil analysis of peat deposits[J]. Mires and Peat, 2010, 7(6): 1-5.

    [11]郎惠卿, 金樹仁. 泥炭植物殘體的研究方法[J]. 植物學通報, 1984, 2(5): 49-51.

    [12]Turner TE, Swindles GT, Charman DJ, et al. Comparing regional and supra-regional transfer functions for palaeohydrological reconstruction from Holocene peatlands[J]. Palaeogeography Palaeoclimatology Palaeoecology, 2013, 369(1): 395-408.

    [13]焦文建. 小興安嶺泥炭蘚泥炭地有殼變形蟲—環(huán)境因子轉(zhuǎn)換函數(shù)模型研究[D]. 長春: 東北師范大學, 2012.

    [14]李鴻凱, 李微微, 蒲有寶, 等. 應用rioja軟件包建立有殼變形蟲-環(huán)境因子轉(zhuǎn)換函數(shù)[J]. 地理科學, 2013, 33(8): 1022-1028.

    [15]秦養(yǎng)民, 謝樹成, 顧延生, 等. 第四紀環(huán)境重建的良好代用指標-有殼變形蟲記錄與古生態(tài)學研究進展[J]. 地球科學進展, 2008, 23(8): 803-812.

    [16]宋理洪.三江平原泥炭沼澤有殼變形蟲多樣性及其在古水位重建中的作用[D]. 長春: 中國科學院研究生院東北地理與農(nóng)業(yè)生態(tài)研究所, 2015.

    [17]孫榮濤, 李鐵剛, 常鳳鳴, 等. 基于底棲有孔蟲轉(zhuǎn)換函數(shù)的南黃海近岸古環(huán)境演化過程定量重建[J]. 中國科學:地球科學, 2015, 45(9): 1383-1393.

    [18]張新榮, 胡克, 方石, 等. 東北泥炭表土沉積植硅體—氣候因子轉(zhuǎn)換函數(shù)建立及應用[J]. 沉積學報, 2008, 26(4): 676-682.

    [19]朱誠, 陳星, 馬春梅, 等. 神農(nóng)架大九湖孢粉氣候因子轉(zhuǎn)換函數(shù)與古氣候重建[J]. 科學通報, 2008, 53(S1): 38-44.

    [20]朱忠秋. 基于有殼變形蟲的大興安嶺泥炭地水位定量重建[D]. 長春: 東北師范大學, 2015.

    [21]莊陳程. 利用長江口現(xiàn)代潮灘硅藻建立海平面轉(zhuǎn)換函數(shù)的嘗試性研究[D]. 上海: 華東師范大, 2012.

    [22]李珍, 王開發(fā), 王永吉, 等.紅樹林孢粉-氣候因子轉(zhuǎn)換函數(shù)恢復古環(huán)境的可行性初探[J]. 海洋科學進展, 2002, 20(3): 73-78.

    [23]Bellen SV, Mauquoy D, Payne RJ, et al. An alternative approach to transfer functions? Testing the performance of a functional trait-based model for testate amoebae[J]. Palaeogeography Palaeoclimatology Palaeoecology, 2017, 468:173-183.

    [24]Mariusz G, Kazimierz T,ukasz L, et al. Unveiling exceptional Baltic bog ecohydrology, autogenic succession and climate change during the last 2000 years in CE Europe using replicate cores, multi-proxy data and functional traits of testate amoebae[J]. Quaternary Science Reviews, 2017, 156:90-106.

    [25]Panait A, Diaconu A, Galka M, et al. Hydrological conditions and carbon accumulation rates reconstructed from a mountain raised bog in the Carpathians: A multi-proxy approach[J]. Catena, 2017, 152:57-68.

    [26]Payne RJ, Kishaba K, Blackford JJ, et al. Ecology of testate amoebae (Protista) in South-Central Alaska peatlands: building transfer-function models for palaeoenvironmental studies[J]. The Holocene, 2006, 16(3): 403-414.

    [27]Payne RJ, Telford RJ, Blackford JJ, et al. Testing peatland testate amoeba transfer functions: Appropriate methods for clustered training-sets[J]. The Holocene, 2012, 22(7): 819-825.

    [28]V?liranta M, Blundell A, Charman DJ, et al. Reconstructing peatland water tables using transfer functions for plant macrofossils and testate amoebae: A methodological comparison[J]. Quaternary International, 2012, 268(9):34-43.

    [29]Wilmshurst JM, Wiser SK, Charman DJ. Reconstructing Holocene water tables in New Zealand using testate amoebae: differential preservation of tests and implications for the use of transfer functions[J]. The Holocene, 2003, 13(1): 61-72.

    [30]成小英,李世杰. 利用硅藻-pH轉(zhuǎn)換函數(shù)定量重建Llyn Hir湖pH的變化[J]. 中國科學院研究生院學報, 2006, 23(3): 357-363.

    [31]董旭輝, 羊向東, 王榮, 等.長江中下游地區(qū)湖泊硅藻-總磷轉(zhuǎn)換函數(shù)[J]. 湖泊科學, 2006, 18(1): 1-12.

    [32]莊陳程, 戰(zhàn)慶, 王張華.利用長江口現(xiàn)代潮灘沉積硅藻建立海平面轉(zhuǎn)換函數(shù)的嘗試[J]. 古地理學報, 2014, 16(4): 557-568.

    [33]向榮, 劉芳, 孫有斌, 等.沖繩海槽末次冰消期以來的浮游有孔蟲群落變化及5種古溫度轉(zhuǎn)換函數(shù)對比[J]. 地球科學(中國地質(zhì)大學學報), 2011, 36(4): 599-609

    [34]Chang JC, Shulmeister J, Woodward C. A chironomid based transfer function for reconstructing summer temperatures in southeastern Australia[J]. Palaeogeography Palaeoclimatology Palaeoecology, 2015, 423: 109-121.

    [35]Velle G, Telford RJ, Heiri O, et al. Testing intra-site transfer functions: an example using chironomids and water depth[J]. Journal of Paleolimnology, 2012, 48(3): 545-558.

    [36]胡玉, 陳建徽, 王海鵬, 等.基于搖蚊的古環(huán)境和古氣候國內(nèi)外研究進展與展望[J]. 地球科學進展, 2016, 31(8): 870-884.

    [37]Lou YJ, Lu XG, Wang GP, et al. Vegetative zonation patterns in depression and riparian wetlands of the Sanjiang Plain, Northeastern China[J]. African Journal of Biotechnology, 2012, 11(2): 355-365.

    [38]Lou YJ, Wang GP, Lu XG, et al. Zonation of plant cover and environmental factors in wetlands of the Sanjiang Plain, northeast China[J]. Nordic Journal of Botany, 2013, 31(6): 748-756.

    [39]Zhou DM, Luan ZQ, Guo XY, et al. Spatial distribution patterns of wetland plants in relation to environmental gradient in the Honghe National Nature Reserve, Northeast China[J]. Journal of Geographical Sciences, 2012, 22(1): 57-70.

    [40]Birks HH. Plant macrofossils[M]. Tracking environmental change using lake sediments. Netherlands: Kluwer Academic Publishers, 2002: 49-74.

    [41]Birks HH, Birks HJB. Future uses of pollen analysis must include plant macrofossils[J]. Journal of Biogeography, 2000, 27(1): 31-35.

    [42]Wang CL, Zhao HY, Wang GP. Vegetation development and water level changes in Shenjiadian peatland in Sanjiang Plain, Northeast China[J]. Chinese Geographical Science, 2015, 25(4): 451-461.

    [43]Wang CL, Zhao HY, Yu XY, et al. Palaeovegetation of Honghe wetland in Sanjiang Plain as a basis for conservation management and restoration[J]. Ecological Engineering, 2016, 96:79-85.

    [44]韓毅. 哈泥泥炭蘚泥炭記錄的沼澤地表濕度變化[D]. 長春: 東北師范大學, 2009.

    [45]馬云艷. 植物大化石和腐殖化度記錄的近500年來長白山西麓的干濕變化[D]. 長春: 東北師范大學, 2010.

    [46]趙紅艷, 李鴻凱, 韓毅, 等. 長白山西側(cè)哈泥沼澤表面濕度的多指標記錄及其可能的驅(qū)動因素[J]. 第四紀研究, 2014, 34(2): 434-442.

    [47]V?liranta M, Korhola A, Sepp? H, et al. High-resolution reconstruction of wetness dynamics in a southern boreal raised bog, Finland, during the late Holocene: a quantitative approach[J]. The Holocene, 2007, 17(8): 1093-1107.

    (編輯:韓志強)

    Studyonvegetation-waterleveltransferfunctionforwetlandsinSanjiangPlain.

    Wang Chunling1,2, Guo Jinping1, Wang Guoping2, Lou Yanjing2*

    (1.CollegeofUrbanandRuralConstruction,ShanxiAgriculturalUniversity,Taigu030801,China; 2.KeyLaboratoryofWetlandEcologyandEnvironment,NortheastInstituteofGeographyandAgroecology,ChineseAcademyofSciences,Changchun130102,China)

    X14

    A

    1671-8151(2017)10-0719-07

    2016-06-30

    2016-08-16

    王春玲(1984-),女(漢),山西忻州人,講師,博士,研究方向:濕地生態(tài)系統(tǒng)與環(huán)境

    *通信作者:婁彥景,副研究員,碩士生導師,Tel: 0431-85542272; E-mail: louyj@neigae.ac.cn

    國家自然科學基金(41401099);山西農(nóng)業(yè)大學科技創(chuàng)新基金項目(2015YJ21、2016015)

    猜你喜歡
    樣方植被物種
    藥王谷多種因素對漏蘆、野韭生長的影響
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    基于植被復綠技術(shù)的孔植試驗及應用
    昆明金殿國家森林公園林業(yè)有害生物調(diào)查研究
    典型野生刺梨(Rosa roxburghii Tratt.)灌草叢植物多樣性研究
    云南九龍山自然保護區(qū)水青岡群落物種多樣性及群落特征
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    綠色植被在溯溪旅游中的應用
    基于原生植被的長山群島植被退化分析
    极品教师在线视频| 男女边吃奶边做爰视频| 美女黄网站色视频| 最新中文字幕久久久久| 精品一区二区免费观看| 在线免费观看的www视频| 少妇的逼好多水| 十八禁网站网址无遮挡 | 成人国产麻豆网| 久久久久久久久久久免费av| 国产精品一区二区三区四区免费观看| 国产精品一区二区在线观看99 | 日韩大片免费观看网站| 亚洲怡红院男人天堂| 亚洲精品色激情综合| 国产高清国产精品国产三级 | 高清日韩中文字幕在线| 免费观看无遮挡的男女| 国产在线一区二区三区精| 黄色配什么色好看| 精品久久久久久久久亚洲| 成人漫画全彩无遮挡| 麻豆成人午夜福利视频| 看黄色毛片网站| 中文天堂在线官网| 精品久久久久久久久av| 好男人在线观看高清免费视频| 国产综合懂色| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 国产美女午夜福利| 亚洲三级黄色毛片| 久久久亚洲精品成人影院| 国产av不卡久久| 国产亚洲av嫩草精品影院| 中文在线观看免费www的网站| 伦理电影大哥的女人| 91av网一区二区| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产| 国产伦精品一区二区三区四那| 国内精品美女久久久久久| 一级毛片 在线播放| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 一级毛片黄色毛片免费观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人av在线播放网站| 日韩视频在线欧美| 麻豆久久精品国产亚洲av| 国产成人91sexporn| 一区二区三区免费毛片| 日日啪夜夜爽| 身体一侧抽搐| 国产亚洲最大av| 欧美xxxx性猛交bbbb| 大又大粗又爽又黄少妇毛片口| 大话2 男鬼变身卡| 久久精品综合一区二区三区| 欧美丝袜亚洲另类| 日韩,欧美,国产一区二区三区| 国产 一区精品| 综合色丁香网| 亚洲不卡免费看| 精品酒店卫生间| 亚洲av成人精品一二三区| 色播亚洲综合网| 蜜桃久久精品国产亚洲av| 大香蕉久久网| 亚洲一级一片aⅴ在线观看| 午夜免费激情av| 观看免费一级毛片| 一区二区三区高清视频在线| 深爱激情五月婷婷| 欧美成人精品欧美一级黄| 国产一区二区三区综合在线观看 | 2021天堂中文幕一二区在线观| 青青草视频在线视频观看| 女人被狂操c到高潮| 伊人久久精品亚洲午夜| 国产一区二区三区综合在线观看 | 2021天堂中文幕一二区在线观| 青青草视频在线视频观看| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 国产亚洲5aaaaa淫片| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 国产伦在线观看视频一区| 日本与韩国留学比较| 国产综合精华液| 欧美成人精品欧美一级黄| 精品酒店卫生间| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 可以在线观看毛片的网站| 久久热精品热| 最近中文字幕高清免费大全6| 99热这里只有精品一区| 永久网站在线| 国产国拍精品亚洲av在线观看| 国产精品久久久久久精品电影小说 | 亚洲av电影在线观看一区二区三区 | 91精品伊人久久大香线蕉| 大香蕉97超碰在线| 免费观看精品视频网站| 在线天堂最新版资源| 色吧在线观看| 一区二区三区免费毛片| 真实男女啪啪啪动态图| 亚洲第一区二区三区不卡| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽人人片av| av在线亚洲专区| 国产 一区精品| 毛片一级片免费看久久久久| 久久精品熟女亚洲av麻豆精品 | 中文字幕人妻熟人妻熟丝袜美| 欧美不卡视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 成人高潮视频无遮挡免费网站| 热99在线观看视频| 五月天丁香电影| 日韩电影二区| 黄片无遮挡物在线观看| 插阴视频在线观看视频| 久久99热6这里只有精品| 精品欧美国产一区二区三| 男的添女的下面高潮视频| 99久久中文字幕三级久久日本| 国产精品女同一区二区软件| 69人妻影院| 六月丁香七月| 色视频www国产| 十八禁国产超污无遮挡网站| 日日摸夜夜添夜夜添av毛片| 国内精品美女久久久久久| 久久99热这里只有精品18| 青青草视频在线视频观看| 欧美xxxx性猛交bbbb| 99热这里只有是精品50| 一个人看的www免费观看视频| 国产三级在线视频| 国产精品日韩av在线免费观看| 午夜老司机福利剧场| 91狼人影院| 国产成人福利小说| 欧美激情在线99| 如何舔出高潮| 伦理电影大哥的女人| 久久久久精品久久久久真实原创| 一级毛片久久久久久久久女| 一级黄片播放器| a级一级毛片免费在线观看| 午夜福利网站1000一区二区三区| 一级a做视频免费观看| 亚洲欧美一区二区三区国产| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区| 免费看美女性在线毛片视频| 18禁动态无遮挡网站| 免费无遮挡裸体视频| 国产伦精品一区二区三区视频9| 一本久久精品| 国产 一区 欧美 日韩| 少妇高潮的动态图| 亚洲欧美精品专区久久| 高清日韩中文字幕在线| 日日摸夜夜添夜夜爱| 成人二区视频| 中文资源天堂在线| 亚洲激情五月婷婷啪啪| 欧美精品国产亚洲| 99久国产av精品国产电影| 国产在视频线在精品| 免费大片黄手机在线观看| 亚洲熟女精品中文字幕| 最近视频中文字幕2019在线8| 久久精品夜色国产| 国产在线男女| 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| videos熟女内射| 久久精品人妻少妇| 久久99热这里只频精品6学生| 国产精品一二三区在线看| 国产精品久久久久久久电影| 一个人观看的视频www高清免费观看| 中文字幕免费在线视频6| 一夜夜www| 亚洲欧美一区二区三区国产| a级一级毛片免费在线观看| 亚洲精品日韩在线中文字幕| 国产男人的电影天堂91| 日韩av在线大香蕉| 免费看av在线观看网站| 中文资源天堂在线| 婷婷色麻豆天堂久久| 久久久久久久午夜电影| 亚洲av免费高清在线观看| 久久久a久久爽久久v久久| 你懂的网址亚洲精品在线观看| 久久久久久伊人网av| av国产久精品久网站免费入址| 内地一区二区视频在线| 精品酒店卫生间| 免费播放大片免费观看视频在线观看| 欧美zozozo另类| 亚洲天堂国产精品一区在线| 美女大奶头视频| 男人和女人高潮做爰伦理| 亚洲在久久综合| 欧美精品国产亚洲| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 赤兔流量卡办理| 亚洲av成人精品一二三区| 欧美xxxx性猛交bbbb| 欧美最新免费一区二区三区| 最近视频中文字幕2019在线8| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜| 在线免费十八禁| 欧美激情国产日韩精品一区| a级毛片免费高清观看在线播放| 成人亚洲精品av一区二区| 成人毛片60女人毛片免费| 亚洲aⅴ乱码一区二区在线播放| 欧美97在线视频| 欧美日韩精品成人综合77777| 久99久视频精品免费| 一个人看的www免费观看视频| 亚洲精品国产成人久久av| 又爽又黄a免费视频| 国产黄频视频在线观看| 国产精品一区www在线观看| 国产白丝娇喘喷水9色精品| av一本久久久久| 老师上课跳d突然被开到最大视频| 日产精品乱码卡一卡2卡三| 插阴视频在线观看视频| 免费大片黄手机在线观看| 国产精品日韩av在线免费观看| 日日啪夜夜爽| 成人国产麻豆网| 三级国产精品片| 永久网站在线| 卡戴珊不雅视频在线播放| 国产精品人妻久久久久久| 一级毛片 在线播放| 免费看日本二区| 边亲边吃奶的免费视频| av国产久精品久网站免费入址| 国产成人一区二区在线| 久久久久久国产a免费观看| 99久久人妻综合| 欧美日韩在线观看h| 国产精品一区www在线观看| 久久久久久久久久人人人人人人| 欧美激情在线99| 免费看不卡的av| 国产精品人妻久久久影院| 亚洲av不卡在线观看| 亚洲精品乱码久久久久久按摩| 一本久久精品| 国产成人精品福利久久| 一个人免费在线观看电影| 99久久精品一区二区三区| 美女国产视频在线观看| 高清毛片免费看| 亚洲自拍偷在线| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 久久午夜福利片| 亚洲性久久影院| 色哟哟·www| 丝袜喷水一区| 亚洲国产精品专区欧美| 欧美 日韩 精品 国产| 91午夜精品亚洲一区二区三区| 久久午夜福利片| 人妻制服诱惑在线中文字幕| 亚洲精品久久午夜乱码| 内射极品少妇av片p| 色哟哟·www| 观看免费一级毛片| 日韩大片免费观看网站| 国产综合懂色| 我的女老师完整版在线观看| 欧美日韩亚洲高清精品| 97人妻精品一区二区三区麻豆| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 2021少妇久久久久久久久久久| 亚洲成色77777| 嫩草影院新地址| 免费看不卡的av| 高清av免费在线| 国产亚洲av片在线观看秒播厂 | 欧美激情在线99| 亚洲国产日韩欧美精品在线观看| a级毛片免费高清观看在线播放| 免费少妇av软件| 精品久久久久久电影网| 久久久久久久久久久丰满| 婷婷色综合大香蕉| 麻豆精品久久久久久蜜桃| 国产精品综合久久久久久久免费| 国产老妇伦熟女老妇高清| 在线天堂最新版资源| 六月丁香七月| 老司机影院毛片| 69av精品久久久久久| 久久久久免费精品人妻一区二区| 我要看日韩黄色一级片| 80岁老熟妇乱子伦牲交| 亚洲av福利一区| 在线观看人妻少妇| 美女黄网站色视频| 日本一本二区三区精品| 日韩三级伦理在线观看| 日韩视频在线欧美| 亚洲内射少妇av| 天堂网av新在线| 国产成人免费观看mmmm| 天堂影院成人在线观看| 亚洲av男天堂| videos熟女内射| 精品久久久久久电影网| 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| 亚洲精品一区蜜桃| 久久久久久久久大av| 别揉我奶头 嗯啊视频| 在线观看一区二区三区| 中文字幕免费在线视频6| 日本午夜av视频| 久久99热6这里只有精品| 亚洲av国产av综合av卡| 久久99热6这里只有精品| 午夜日本视频在线| 一级黄片播放器| 午夜日本视频在线| 成年免费大片在线观看| 亚洲欧美精品专区久久| 免费看日本二区| 深夜a级毛片| 女人被狂操c到高潮| 黄片无遮挡物在线观看| 如何舔出高潮| 精品国产三级普通话版| 狠狠精品人妻久久久久久综合| 精华霜和精华液先用哪个| 嘟嘟电影网在线观看| 国产精品福利在线免费观看| 日本黄色片子视频| 精华霜和精华液先用哪个| 蜜桃亚洲精品一区二区三区| 久久久久精品久久久久真实原创| 亚洲精品一区蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 婷婷色麻豆天堂久久| 国产精品三级大全| 日韩欧美一区视频在线观看 | 国产成人freesex在线| 一级毛片久久久久久久久女| 成人高潮视频无遮挡免费网站| 69人妻影院| 久久鲁丝午夜福利片| 啦啦啦中文免费视频观看日本| 国产一区二区三区综合在线观看 | 一边亲一边摸免费视频| 日韩欧美一区视频在线观看 | 成人亚洲精品一区在线观看 | 亚洲一级一片aⅴ在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲人成网站在线播| 国内揄拍国产精品人妻在线| 久热久热在线精品观看| 久久精品久久久久久噜噜老黄| 18禁动态无遮挡网站| 日韩精品有码人妻一区| 精品亚洲乱码少妇综合久久| 亚洲精品乱久久久久久| 啦啦啦韩国在线观看视频| 麻豆国产97在线/欧美| 亚洲丝袜综合中文字幕| 舔av片在线| 在线播放无遮挡| 国产综合精华液| 国产一区二区三区综合在线观看 | 少妇裸体淫交视频免费看高清| 高清视频免费观看一区二区 | 国产有黄有色有爽视频| 麻豆成人午夜福利视频| 国产极品天堂在线| 如何舔出高潮| 亚洲av男天堂| 久久精品熟女亚洲av麻豆精品 | 国产成人精品一,二区| 国产一区亚洲一区在线观看| 免费人成在线观看视频色| 色综合站精品国产| 91狼人影院| 啦啦啦中文免费视频观看日本| 伦理电影大哥的女人| 国模一区二区三区四区视频| 日本-黄色视频高清免费观看| 国产片特级美女逼逼视频| 男的添女的下面高潮视频| 国产中年淑女户外野战色| av播播在线观看一区| 国产精品蜜桃在线观看| 久久久久久久久大av| 亚洲在线观看片| 97超视频在线观看视频| 日本三级黄在线观看| 日本色播在线视频| av国产久精品久网站免费入址| 搡老妇女老女人老熟妇| 国产日韩欧美在线精品| 自拍偷自拍亚洲精品老妇| 在现免费观看毛片| 三级经典国产精品| 男女啪啪激烈高潮av片| 99久国产av精品国产电影| 一级毛片黄色毛片免费观看视频| 日韩欧美精品免费久久| 大香蕉久久网| 国产黄片视频在线免费观看| 成人欧美大片| www.色视频.com| 99久国产av精品| 国产精品久久久久久久电影| 亚洲精品国产成人久久av| 精品国内亚洲2022精品成人| 国产精品国产三级专区第一集| 91在线精品国自产拍蜜月| 熟女电影av网| 嫩草影院新地址| 丝袜喷水一区| av国产免费在线观看| 春色校园在线视频观看| 黄色日韩在线| 欧美激情在线99| 不卡视频在线观看欧美| 国产成年人精品一区二区| 高清毛片免费看| 汤姆久久久久久久影院中文字幕 | 久久久久九九精品影院| 国产综合懂色| 纵有疾风起免费观看全集完整版 | 91精品国产九色| 国产av在哪里看| 国产精品一区www在线观看| 99久久九九国产精品国产免费| 久久久精品94久久精品| 中文字幕制服av| 三级国产精品欧美在线观看| 男女视频在线观看网站免费| 国产又色又爽无遮挡免| 久久人人爽人人爽人人片va| 亚洲成人一二三区av| 国产黄色免费在线视频| 亚洲精品影视一区二区三区av| 久久久精品免费免费高清| xxx大片免费视频| 99久久人妻综合| av播播在线观看一区| 精品欧美国产一区二区三| 黄色日韩在线| 性插视频无遮挡在线免费观看| 看免费成人av毛片| 搡老妇女老女人老熟妇| 赤兔流量卡办理| 在线免费观看的www视频| 老师上课跳d突然被开到最大视频| 日韩不卡一区二区三区视频在线| 一级黄片播放器| 久久热精品热| 欧美日韩国产mv在线观看视频 | 国产片特级美女逼逼视频| 亚洲电影在线观看av| 最近手机中文字幕大全| 精品国产一区二区三区久久久樱花 | 欧美高清性xxxxhd video| 黄片wwwwww| 免费黄网站久久成人精品| 九九久久精品国产亚洲av麻豆| 亚洲av在线观看美女高潮| 久久99热这里只有精品18| av黄色大香蕉| 久久久久久久久久黄片| 国产免费又黄又爽又色| 2022亚洲国产成人精品| 99热6这里只有精品| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 亚洲激情五月婷婷啪啪| 国产亚洲一区二区精品| 国产av在哪里看| 青春草国产在线视频| 蜜臀久久99精品久久宅男| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| 国产黄色视频一区二区在线观看| 国产永久视频网站| 能在线免费观看的黄片| 日韩人妻高清精品专区| 97精品久久久久久久久久精品| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影小说 | 99九九线精品视频在线观看视频| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 中国美白少妇内射xxxbb| 18禁在线播放成人免费| 亚洲,欧美,日韩| 校园人妻丝袜中文字幕| 亚洲成人av在线免费| 97精品久久久久久久久久精品| 最新中文字幕久久久久| 国产激情偷乱视频一区二区| 日韩av在线大香蕉| 国产乱人视频| 国产在线一区二区三区精| 91久久精品国产一区二区三区| 女人久久www免费人成看片| 午夜爱爱视频在线播放| 成人国产麻豆网| 大又大粗又爽又黄少妇毛片口| 91在线精品国自产拍蜜月| 国产熟女欧美一区二区| 一本一本综合久久| 成人毛片a级毛片在线播放| 成年人午夜在线观看视频 | 国产精品一及| 精品99又大又爽又粗少妇毛片| 国产三级在线视频| 午夜精品国产一区二区电影 | 国产不卡一卡二| 又爽又黄无遮挡网站| 亚洲av国产av综合av卡| 亚洲国产精品成人久久小说| 成年女人看的毛片在线观看| 美女内射精品一级片tv| 美女脱内裤让男人舔精品视频| 亚洲国产精品成人综合色| 啦啦啦韩国在线观看视频| www.色视频.com| 国产 一区 欧美 日韩| av免费在线看不卡| 亚洲人成网站高清观看| 国产精品一区二区三区四区久久| 成人综合一区亚洲| 少妇的逼好多水| 91精品一卡2卡3卡4卡| 日韩成人av中文字幕在线观看| 久久久久国产网址| 性插视频无遮挡在线免费观看| 三级国产精品欧美在线观看| 久久久久性生活片| 神马国产精品三级电影在线观看| 在线免费十八禁| 国产高潮美女av| 午夜精品在线福利| 插逼视频在线观看| 毛片一级片免费看久久久久| 网址你懂的国产日韩在线| 两个人视频免费观看高清| kizo精华| 亚洲欧洲日产国产| 国产精品不卡视频一区二区| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区| 一个人观看的视频www高清免费观看| 我的老师免费观看完整版| 美女高潮的动态| av免费在线看不卡| 97超碰精品成人国产| 在线观看免费高清a一片| 两个人的视频大全免费| 又爽又黄无遮挡网站| 99久久中文字幕三级久久日本| 午夜免费观看性视频| 一级毛片久久久久久久久女| 内射极品少妇av片p| 欧美成人a在线观看| 91久久精品电影网| 超碰97精品在线观看| 亚洲一级一片aⅴ在线观看| 在线播放无遮挡| 欧美xxxx黑人xx丫x性爽| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av蜜桃| 国产精品人妻久久久影院| 色5月婷婷丁香| 一级毛片我不卡| 婷婷色麻豆天堂久久| 亚洲图色成人| 丰满乱子伦码专区| h日本视频在线播放| 熟女人妻精品中文字幕| 久久精品久久久久久噜噜老黄| 天堂网av新在线| 激情五月婷婷亚洲| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 国产极品天堂在线| 欧美极品一区二区三区四区| 国产一区亚洲一区在线观看| 内射极品少妇av片p| 国产精品久久视频播放|