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

    利用馬尾松解析木對林地立地指數(shù)和立地形的研究

    2021-08-16 01:50:54李清順王得軍孫景梅
    西北林學(xué)院學(xué)報 2021年4期
    關(guān)鍵詞:馬尾松胸徑生長量

    李清順,王得軍,孫景梅,劉 勇

    (1.國家林業(yè)和草原局 西北調(diào)查規(guī)劃設(shè)計院,陜西 西安 710048;2.旱區(qū)生態(tài)水文與災(zāi)害防治國家林業(yè)和草原局 局重點實驗室,陜西 西安 710048)

    立地質(zhì)量是指在某一立地上既定森林或者其他植被類型的生產(chǎn)潛力,對林地立地質(zhì)量的有效科學(xué)評價是林地管理和經(jīng)營的基礎(chǔ)。目前,最為流行的立地質(zhì)量評價方法是采用一定年齡下的優(yōu)勢高,即立地指數(shù)[1-4]。為了獲取立地指數(shù),需要測得立木的年齡和優(yōu)勢高數(shù)據(jù)進而建立年齡-優(yōu)勢高方程。但是,由于生長環(huán)境和過程等因素的影響,只有年齡大、同齡、保存良好且未受破壞的純林或單一優(yōu)勢木的林分才能達到建立模型的要求[5-7]。正是由于這些限制,導(dǎo)致立地指數(shù)的應(yīng)用受到影響。在此情況下,立地形(或立地產(chǎn)量)作為一種替代立地指數(shù)的方法被提出[8]。立地形定義為指定胸徑下的優(yōu)勢高或亞優(yōu)勢高平均值,立地形的提出基于2個假設(shè):(1)立地形隨著削度(胸徑樹高比)的降低而升高;(2)在異齡或混交林分中,林分密度并不影響優(yōu)勢木或亞優(yōu)勢木的胸徑-樹高關(guān)系[8-9]。立地形從被提出就一直備受爭議,很多學(xué)者依據(jù)立地形對一定地區(qū)的立地質(zhì)量作了評價,并認為該方法對于異齡、混交林有效[10-14],而另外的學(xué)者則得出了相反的結(jié)論[15-16]。

    對于立地質(zhì)量的評價一般采用固定樣地、臨時樣地或者解析木數(shù)據(jù)。不管采用哪種方法,年齡的獲取和精度保證都是評價的關(guān)鍵。因此,采用解析木數(shù)據(jù)進行評價是比較準確和客觀的。本研究采用重慶市南川區(qū)馬尾松解析木128株、圓盤記錄1 282對,建立了馬尾松立地指數(shù)和立地形表,同時,對立地指數(shù)和立地形作了對比驗證,并對使用立地形對馬尾松立地質(zhì)量評價的可行性進行了評價。在此基礎(chǔ)上給出了使用立地形對立地質(zhì)量進行評價的條件,從而避免了對立地形進行立地質(zhì)量進行評價的直接否定或肯定。

    1 材料與方法

    1.1 研究區(qū)概況

    南川區(qū)位于重慶市南部,介于28°46′-29°30′N,106°54′-107°27′E,南與貴州省道真縣、正安縣、桐梓縣接壤,東北與武隆縣為鄰,北接涪陵區(qū),西連巴南區(qū)、綦江區(qū)、萬盛經(jīng)濟開發(fā)區(qū)。南川區(qū)地處四川盆地與云貴高原的過渡地帶,地質(zhì)構(gòu)造復(fù)雜,平均海拔340~2 251 m;屬亞熱帶濕潤性季風(fēng)氣候區(qū),氣候溫和,全區(qū)年平均氣溫16.6℃,常年平均降雨量1 185 mm。南川區(qū)的喬木林樹種以馬尾松(Pinusmassoniana)、柏木(Cupressusfunebris)、杉木(Cunninghamialanceolata)為主,其中馬尾松面積最大,占整個喬木林面積的67.7%。

    1.2 數(shù)據(jù)來源

    研究區(qū)數(shù)據(jù)來源為2017年南川區(qū)森林資源二類調(diào)查樣地補充調(diào)查數(shù)據(jù)。樣地采用大小為1畝(邊長為25.82 m的方形樣地)的臨時樣地,臨時樣地參考前期確認的馬尾松小班分布圖,按立地類型全覆蓋進行設(shè)置,保證每一種立地類型有3個以上的樣地。共設(shè)置臨時樣地35個,樣地海拔472~1 583 m,坡度2°~37°,坡向以陽坡為主(近70%為陽坡,陽坡包括南坡、東南坡和西南坡)。每個樣地選擇生長正常、無病蟲害、無斷梢的優(yōu)勢木3~5株。優(yōu)勢木在伐倒前,應(yīng)先準確確定根頸位置和實測胸徑,并在樹干上標明胸高直徑的位置和正南正北方向。伐倒后,先測定由根頸至第1個死枝和活枝在樹干上的高度,然后打去枝椏,在全樹干上用紅色標明南北方向。并在全樹高1/4、1/2、3/4處做上紅色標記,然后測量帶皮和去皮直徑(去皮直徑結(jié)合圓盤截取進行)。同時,按2 m進行造材,用藍色在每2 m一段處做上標記,測量帶皮和去皮直徑(去皮直徑結(jié)合圓盤截取進行)。結(jié)合樹干全長測定,采用中央斷面區(qū)分求積法,按1 m(樹干全長3~12 m)或2 m(樹干長度在12 m以上)對樹干進行區(qū)分,并做好標記,標記處是圓盤段面高處。不足一個區(qū)分段長度的樹干為梢頭木。需要注意的是:胸徑與梢頭底直徑位置處必須截取圓盤,做好標記。最終獲取馬尾松解析木128株,圓盤記錄1 282對。通過散點圖和三倍標準差進行控制,刪除雙河林場187號樣地4號解析木10個圓盤記錄對,然后進行描述性統(tǒng)計(表1)。

    表1 優(yōu)勢木特征值

    1.3 立地指數(shù)和立地形模型的建立

    1.3.1 數(shù)據(jù)劃分 根據(jù)研究需要,依據(jù)鄉(xiāng)鎮(zhèn)、村全覆蓋原則,進行劃分,4/5用以建模,1/5用以驗證。劃分結(jié)果102株解析木(包含圓盤記錄1 021對)用以建模,25株解析木(包含圓盤記錄251對)用于驗證。

    1.3.2 導(dǎo)向曲線的擬合 導(dǎo)向曲線是指中等立地條件下,林分優(yōu)勢高隨年齡變化的平均生長曲線。導(dǎo)向曲線的選擇直接影響立地質(zhì)量評價的精度,這就需要導(dǎo)向曲線不僅能反映生物的生長規(guī)律,同時要能對生長過程做出最佳擬合。本研究采用了5個常用的理論生長方程Logistic、Mitscherlich、Gompertz、Korf和 Richard[17](表2)。另外,為了保證立地指數(shù)和立地形對比的一致性,年齡-樹高和直徑-樹高采用了相同的生長模型。表2中AD表示年齡或直徑(胸徑),a、b、c為模型參數(shù)。模型根據(jù)決定系數(shù)R2,均方根誤差RMES、AIC及方程的參數(shù)變動系數(shù)進行選擇。

    表2 樹高生長備選模型

    1.3.3 基準年齡、基準胸徑和指數(shù)級距的確定 確定基準年齡和基準胸徑就是尋找樹高生長趨于穩(wěn)定且能夠靈敏反映立地差異的年齡?;鶞誓挲g和基準胸徑對于立地指數(shù)表或立地形表的準確編制影響十分明顯。其一般考慮以下3點:(1)樹高生長趨于穩(wěn)定后的一個齡(徑)階;(2)采伐年齡(胸徑);(3)自然成熟年齡(胸徑)的1/2[18-19]。樹高的連續(xù)生長量和平均生長量曲線在樹高生長趨于穩(wěn)定后,將會出現(xiàn)相交,相交后的后一個齡階的年齡(胸徑)可以作為基準年齡(胸徑)。樹高的連續(xù)生長量定義為樹木某一年(或1 cm)的實際生長速度,連續(xù)生長量是樹木年齡t(或胸徑D)的函數(shù)。其生長方程為(方程中t換為D,即為以胸徑為單位的連續(xù)生長量):

    (1)

    樹高平均生長量說明樹木在某一時刻t(或某一胸徑D)的平均生長速度,其方程為(同樣t換成D,即為以胸徑為單位的平均生長量):

    (2)

    式中,H(t)為年齡為t時的樹高,H(t0)為年齡為0時的理論樹高。由于所選用生長方程在年齡或胸徑為0時,樹高并不一定為0,因此方程中要減去理論樹高。

    指數(shù)級距的確定根據(jù)基準年齡(胸徑)時的樹高絕對變幅和指數(shù)級個數(shù)確定。公式如下:

    CV=Si/HI

    (3)

    ΔCV=CVi+1/CVi

    (4)

    C=ΔH/k

    (5)

    式中,CV為變異系數(shù);ΔCV為變異系數(shù)變化幅度;Si為第i齡(徑)階樹高標準差;HI為第i齡(徑)階主曲線樹高值;C為指數(shù)級距;ΔH為基準年齡(胸徑)時的樹高絕對變幅;k為指數(shù)級個數(shù)。

    1.3.4 立地指數(shù)表和立地形表的編制 立地指數(shù)表的編制方法主要有標準差調(diào)整法、變動系數(shù)調(diào)整法、相對優(yōu)勢高法。很多學(xué)者對3種編制方法的優(yōu)劣性進行了對比,并對其適用性作了說明[21-22]。本研究采用比較流行的標準差調(diào)整法進行編表,所用公式如下:

    (6)

    式中,Hij為第i齡階第j指數(shù)級調(diào)整后的樹高;Hik為第i齡階的導(dǎo)向曲線樹高;H0j為基準年齡時第j指數(shù)級的樹高;H0k為基準年齡時的導(dǎo)向曲線樹高;SA0為基準年齡所在齡階樹高標準差理論值;SAi為第i齡階樹高標準差理論值。

    (7)

    (8)

    (9)

    2 結(jié)果與分析

    2.1 導(dǎo)向曲線的擬合結(jié)果

    采用SPSS26軟件對所選5個備選方程進擬合,擬合結(jié)果見表3。

    表3 導(dǎo)向曲線擬合結(jié)果

    在同一坐標系下分別作出年齡-樹高和胸徑-樹高曲線(圖1、圖2)[23]。

    圖1 年齡-樹高擬合曲線

    從表3可以看出,無論從決定系數(shù)、均方根誤差還是從AIC結(jié)果來看,年齡-樹高的方程M4(Korf)和M5(Richard)的擬合效果都優(yōu)于其他3個模型。查看方程M4和M5的3個參數(shù)的變動系數(shù),方程M4參數(shù)變動系數(shù)分別為34.74%、4.24%、16.22%,方程M5的參數(shù)變動系數(shù)分別為7.44%、15.16%、7.53。從參數(shù)變動系數(shù)來看,方程5的穩(wěn)定性要優(yōu)于方程M4。再參照圖1的方程擬合曲線,方程M1(Logistic)和方程M3(Gompertz)在幼齡階段過高擬合,而在成、過熟林階段存在過低擬合;而方程M2(Mitscherlich)擬合結(jié)果整體偏高;方程M4在幼齡林階段擬合偏低,在成、過熟林階段擬合偏高。

    從表3來看,對胸徑-樹高模型方程M3(Gompertz)從各項指標來看都優(yōu)于其他4個方程。再參照圖2的擬合曲線,方程M1(Logistic)在幼齡階段過高擬合,而在成、過熟林階段存在過低擬合;方程M2(Mitscherlich)整體擬合偏高;方程M4(Korf)、M5(Richard)在幼齡林階段擬合偏低,在成、過熟林階段擬合偏高。

    綜上所述,年齡-樹高選用方程M5(Richard),胸徑-樹高選用方程M3(Gompertz),具體方程如下:

    H=22.083(1-e-0.041t)1.417

    (10)

    H=17.577e-4.565e-0.171D

    (11)

    圖2 胸徑-樹高擬合曲線

    2.2 基準年齡、胸徑與指數(shù)級距

    根據(jù)年齡-樹高模型、胸徑-樹高模型,采用方程(1)、(2)分別求出馬尾松生長量曲線方程如下:

    Z(t)=1.283e-0.041t(1-e-0.041t)0.417

    (12)

    (13)

    Z(D)=13.721e(-4.565e-0.171D-0.171D)

    (14)

    (15)

    根據(jù)方程(12)、(13)和方程(14)、(15)利用Matlab軟件分別繪出馬尾松生長量曲線圖(圖3、圖4)。從圖4可以看出,馬尾松在8 a時連年生長量達到最大,在16 a時樹高生長趨于穩(wěn)定,平均生長量達到最大,并且與連年生長量出現(xiàn)相交,以后平均生長量和連年生長量都下降。16 a后的下一個齡階為20 a,因此,將馬尾松的基準年齡定為20 a。從圖2可以看出,馬尾松在胸徑為9 cm時連續(xù)生長量達到最大,在胸徑接近14 cm時樹高生長趨于穩(wěn)定,平均生長量達到最大,并且與連續(xù)生長量出現(xiàn)相交,以后平均生長量和連續(xù)生長量都下降。因此選擇平均胸徑14 cm作為基準胸徑。

    基準年齡為20 a時,樹高變動范圍為5.7~15.6 m,確定級數(shù)距為2 m,劃分6~16共6個指數(shù)級。基準胸徑為14 cm時,樹高變動范圍為5.6~19.8 m,級數(shù)距仍為2 m,劃分出6~20共8個指數(shù)級。

    2.3 立地指數(shù)表和立地形表的編制

    利用標準差調(diào)整法編制地位指數(shù)表和地位形表(表4、表5),同時作出其曲線簇(圖5、圖6)。

    2.4 模型的評價與檢驗結(jié)果

    根據(jù)編表齡(徑)階范圍,刪除圓盤記錄中低于最小編表齡(徑)階的圓盤記錄,剩余圓盤記錄179對。經(jīng)成對T檢驗算出立地指數(shù)和立地形評價指標結(jié)果(表5)。從表6可以看出,立地指數(shù)95%顯著水平下的T值為1.528,平均相對誤差1.309%,相對總誤差0.03%,預(yù)估精度99.18%;立地形95%顯著水平下的T值為0.903,平均總相對誤差0.181%,相對總誤差0.02%,預(yù)估精度99.22%。查T檢驗表T0.05(178)=1.96,立地指數(shù)和立地形的T值均小于該值,說明立地指數(shù)和立地形在95%水平下顯著。平均總相對誤差和相對總誤差都<5%,而預(yù)估精度分別達到了99.18%和99.22%,說明所編立地指數(shù)表和立地形表精度較高,能夠滿足編表需要。

    圖3 馬尾松生長量曲線(胸徑)

    圖4 馬尾松生長量曲線(年齡)

    圖5 立地指數(shù)曲線

    表4 馬尾松立地指數(shù)表(基準年齡:20 a)

    表5 馬尾松立地形表(基準胸徑:14 cm)

    圖6 立地形曲線

    表6 立地指數(shù)和立地形檢驗結(jié)果

    3 討論

    3.1 基準年齡和基準胸徑的選擇

    基準年齡和基準胸徑的選擇對于立地指數(shù)和立地形模型的影響顯著,如果選擇不當(dāng)或選擇錯誤將嚴重影響立地指數(shù)和立地形的編表質(zhì)量。目前常用的比較流行的方法中,采用樹種生活史(主伐年齡或胸徑)的1/2或成熟年齡(胸徑)的年數(shù),不僅缺乏對當(dāng)?shù)貙嶋H立地情況的考慮,而且缺乏理論基礎(chǔ)。相比較而言,利用某一樹種連續(xù)生長量和平均生長量趨于穩(wěn)定后的年齡(或胸徑),具有堅實的理論基礎(chǔ),而且比較客觀的反應(yīng)了所研究地區(qū)的立地質(zhì)量情況。但是,在采用該方法時,一般都是采用實測值的連續(xù)生長量和變異系數(shù)的變化趨于穩(wěn)定時的年齡(胸徑)[18-20],這樣對于實測值的要求較高,當(dāng)樣本量不夠時很難準確的估計出實際年齡(胸徑)。本研究根據(jù)所確定的導(dǎo)向曲線,直接求出該樹種的連續(xù)生長量和平均生長量方程,然后利用Matlab繪出其生長量曲線,根據(jù)連續(xù)生長量和平均生長量曲線次相交時的年齡(胸徑),采用其下一個齡階或徑階作為基準年齡(胸徑)。該方法具有嚴謹?shù)臄?shù)學(xué)理論基礎(chǔ),而且繪出的曲線圖直觀明了,可以方便準確地確定基準年齡(胸徑)。

    3.2 指數(shù)級的選擇對編表精度的影響

    本研究中立地指數(shù)和立地形的指數(shù)級分別為6級和8級,指數(shù)級別的不同是由于齡階和徑階的不同劃分造成的結(jié)果。齡階采用的是5 a一個齡階,基準年齡為20 a,即20齡階,實際采用的年齡范圍為15~20 a。胸徑采用的是2 cm一個徑階,基準胸徑為14 cm,即14徑階,實際采用的胸徑范圍13~15 cm。在基準年齡和基準胸徑范圍內(nèi)樹高變化范圍為5.7~15.6 m和5.6~19.8 m,因此造成了指數(shù)級別的不同。指數(shù)級別的不同也在一定程度上影響了立地指數(shù)表和立地形表的精度。由于立地形表的級別較多,每一級別所對應(yīng)的變動范圍較窄,編制精度也就越高,這與模型的評價與檢驗結(jié)果相一致。

    3.3 立地形和立地質(zhì)量評價結(jié)果的差異

    由于在評價立地質(zhì)量時采用的是立地指數(shù)(立地形)基于基準年齡(胸徑)時的實際值,而不是其調(diào)整值,因此,為了檢驗2種方法的差異,我們首先采用差分方程法重新計算立地指數(shù)和立地形。所用方程如下:

    (16)

    (17)

    式中,SI為立地指數(shù),H為樹高,t為樹齡,t0為基準年齡;SF為立地形,D為胸徑,D0為基準胸徑。

    根據(jù)以上計算結(jié)果對立地指數(shù)和立地形進行線性擬合,得到的線性方程如下:

    SF=0.631SI+4.215

    (18)

    從方程(16)可以看出,在立地指數(shù)<11.4時,立地形大于立地指數(shù),當(dāng)立地指數(shù)>11.4時,立地形小于立地指數(shù)。

    進一步,用以上5個經(jīng)驗生長模型擬合年齡和胸徑,得到胸徑年齡方程如下:

    D=61.03exp (-4.829t-0.372)

    (19)

    當(dāng)基準年齡為20 a時,胸徑為12.5 cm,明顯小于所選基準胸徑14 cm,此時對應(yīng)的立地指數(shù)為10,而基準胸徑對應(yīng)的立地形為12。

    由此可以看出,當(dāng)樹高生長趨于穩(wěn)定時,由于胸徑生長滯后于樹高,導(dǎo)致基準胸徑的選擇大于對應(yīng)基準年齡時的胸徑,從而導(dǎo)致立地形在生長的不同階段與立地指數(shù)反映出的立地質(zhì)量的差異性。也就是說在樹木生長的早期,由于胸徑的生長率低于樹高生長率,導(dǎo)致立地形偏低,而在樹木生長的中后期由于胸徑的生長率高于樹高的生長率,導(dǎo)致立地形偏高。立地形和立地指數(shù)在評價立地質(zhì)量的這種差異性是在對立地質(zhì)量進行客觀評價時需要注意的。

    3.4 立地形的適用性分析

    考查削度與立地指數(shù)和立地形的相關(guān)性,削度與立地指數(shù)的相關(guān)系數(shù)為-0.569,與立地形的相關(guān)系數(shù)為-0.309。由此可見立地指數(shù)與立地形都與削度呈現(xiàn)負相關(guān)性,即立地質(zhì)量和立地形都隨著削度的增加而降低,只是立地指數(shù)的這種負相關(guān)性更加明顯。由此可以看出,立地形不管是在同齡純林還是異齡混交林種,只要胸徑的生長能夠具有相對獨立性而不受周圍樹木競爭的影響,就可以相對比較客觀的反映立地質(zhì)量,其評價結(jié)果也是可靠的。但是,一旦胸徑的生長受到林分密度等因素的影響,立地形的評價結(jié)果就會出現(xiàn)偏差,從而失去其評價的客觀性。另外,立地質(zhì)量與立地形一個明顯的區(qū)別是一個采用了無量綱參數(shù)t(年齡),一個采用了有量綱參數(shù)D(胸徑)。而年齡是不受外界環(huán)境影響的,而胸徑明顯受到周圍各種因素的影響,因此在采用立地形進行立地質(zhì)量評價時必須首先排除這種影響才能比較客觀的對立地質(zhì)量做出科學(xué)的評價。也就是說當(dāng)胸徑與樹高的相關(guān)性越強,立地形與立地指數(shù)反映出的林地質(zhì)量結(jié)果越接近,立地形也就可以替代立地指數(shù)作為評價立地質(zhì)量的指標。由此給出如下推斷:不管是在同齡純林還是異齡交林,優(yōu)勢木的選擇至關(guān)重要,只要能夠選擇出符合條件的優(yōu)勢木,用立地形代替立地指數(shù)是完全可行的。當(dāng)然,這個推斷需要相關(guān)實驗來驗證,也是未來一個研究的方向。

    4 結(jié)論

    本研究采用重慶市南川區(qū)馬尾松解析木數(shù)據(jù),選擇合適的導(dǎo)向曲線方程,建立了基于導(dǎo)向曲線的立地指數(shù)和立地形模型。結(jié)果顯示立地指數(shù)和立地形表達精度在95%水平下顯著,平均總相對誤差和相對總誤差都<5%,而預(yù)估精度分別達到了99.18%和99.22%,說明所編立地指數(shù)表和立地形表精度較高,能夠滿足編表需要。在此基礎(chǔ)上,對研究過程和結(jié)果進行了分析討論,得出了如下結(jié)論:(1)根據(jù)導(dǎo)向曲線繪制生長量曲線方程的方法確定基準年齡和基準胸徑的方法準確可行;(2)由于齡階和徑階的選擇差異,導(dǎo)致立地指數(shù)級不同,同時影響了立地形與立地指數(shù)的編表質(zhì)量;(3)由于樹木的生長特性,立地形在表示立地質(zhì)量時,在林木的早期將會低于立地指數(shù),在樹木生長的中后期將會高于立地指數(shù);(4)立地形的適用性取決于胸徑與樹高的相關(guān)性,相關(guān)性越強,立地形就越準確,在選擇合適的優(yōu)勢木的情況下,立地形可以取代立地指數(shù)來表示林地質(zhì)量。結(jié)論(3)和(4)是在分析本研究的結(jié)果上提出的,需要相關(guān)的現(xiàn)地研究進行驗證,也是下一步研究的一個方向。

    猜你喜歡
    馬尾松胸徑生長量
    白城山新1號楊育苗密度研究
    武漢5種常見園林綠化樹種胸徑與樹高的相關(guān)性研究
    日本落葉松人工林生長規(guī)律分析
    綠色科技(2021年21期)2021-11-26 09:12:06
    北大河林業(yè)局森林生長量、枯損量調(diào)查與分析
    馬尾松栽培技術(shù)及撫育管理
    綠色科技(2019年5期)2019-11-29 13:17:37
    宜春區(qū)域南方紅豆杉生境及其生長量分析
    福建省森林資源監(jiān)測體系抽樣調(diào)查中胸徑測量精度范圍的精準確定
    馬尾松果糖-1,6-二磷酸酶基因克隆及表達模式分析
    馬尾松初級種子園復(fù)壯技術(shù)
    24年生馬尾松種子園自由授粉子代測定及家系選擇
    午夜免费成人在线视频| 看黄色毛片网站| av一本久久久久| 制服人妻中文乱码| 久久人妻av系列| xxx96com| 亚洲色图综合在线观看| 国产成+人综合+亚洲专区| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费电影在线观看| 啦啦啦免费观看视频1| av不卡在线播放| 国产日韩欧美亚洲二区| 国产在线一区二区三区精| 99re在线观看精品视频| 人妻 亚洲 视频| av有码第一页| 免费人成视频x8x8入口观看| 麻豆国产av国片精品| 成年女人毛片免费观看观看9 | 免费在线观看亚洲国产| 中文字幕色久视频| 久久午夜综合久久蜜桃| 午夜免费鲁丝| 日日摸夜夜添夜夜添小说| 两性午夜刺激爽爽歪歪视频在线观看 | 无限看片的www在线观看| 窝窝影院91人妻| 欧美日韩国产mv在线观看视频| 99久久99久久久精品蜜桃| 精品国内亚洲2022精品成人 | 黄色 视频免费看| 精品一区二区三区视频在线观看免费 | 51午夜福利影视在线观看| 俄罗斯特黄特色一大片| 日韩视频一区二区在线观看| 日本a在线网址| 最新的欧美精品一区二区| 一边摸一边抽搐一进一小说 | 国产精品98久久久久久宅男小说| 丰满饥渴人妻一区二区三| videos熟女内射| www.自偷自拍.com| 国产av一区二区精品久久| 日韩三级视频一区二区三区| 国产精品一区二区在线观看99| 岛国毛片在线播放| 国产一区二区三区在线臀色熟女 | 涩涩av久久男人的天堂| 国产成人精品久久二区二区免费| 一级毛片高清免费大全| 无人区码免费观看不卡| 老司机午夜十八禁免费视频| 成人国语在线视频| 正在播放国产对白刺激| 啦啦啦在线免费观看视频4| 好看av亚洲va欧美ⅴa在| 欧美精品啪啪一区二区三区| 中文亚洲av片在线观看爽 | 天堂动漫精品| 久久久水蜜桃国产精品网| 亚洲精品一二三| 久久国产精品男人的天堂亚洲| 亚洲av成人不卡在线观看播放网| 国产日韩欧美亚洲二区| 村上凉子中文字幕在线| 无遮挡黄片免费观看| 亚洲少妇的诱惑av| 国产不卡av网站在线观看| 亚洲中文av在线| 大香蕉久久网| 一本一本久久a久久精品综合妖精| 久久中文字幕人妻熟女| 中文字幕人妻丝袜一区二区| 人成视频在线观看免费观看| 午夜亚洲福利在线播放| 一级作爱视频免费观看| 老司机午夜福利在线观看视频| 国产日韩一区二区三区精品不卡| 黄色毛片三级朝国网站| 免费日韩欧美在线观看| 丰满饥渴人妻一区二区三| 一级a爱片免费观看的视频| 国产免费av片在线观看野外av| 亚洲精品粉嫩美女一区| 久久香蕉国产精品| videosex国产| av超薄肉色丝袜交足视频| 黄色片一级片一级黄色片| 啦啦啦在线免费观看视频4| 成人av一区二区三区在线看| 午夜福利免费观看在线| 精品国产一区二区久久| 男男h啪啪无遮挡| 久久人人97超碰香蕉20202| 亚洲精品粉嫩美女一区| 成人国语在线视频| 757午夜福利合集在线观看| 国产又色又爽无遮挡免费看| 日韩大码丰满熟妇| 脱女人内裤的视频| 91成人精品电影| av有码第一页| 婷婷丁香在线五月| 国产精品亚洲一级av第二区| 久久香蕉国产精品| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 亚洲精品中文字幕一二三四区| 女人爽到高潮嗷嗷叫在线视频| 免费看a级黄色片| 久久香蕉国产精品| 国产亚洲欧美98| 国产激情欧美一区二区| 乱人伦中国视频| 黑丝袜美女国产一区| 99热网站在线观看| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 精品亚洲成国产av| 狠狠狠狠99中文字幕| 久久久国产一区二区| 欧美日韩精品网址| 中文字幕人妻丝袜一区二区| 国产成人一区二区三区免费视频网站| av电影中文网址| 亚洲国产毛片av蜜桃av| 亚洲五月天丁香| 亚洲精品国产色婷婷电影| 黄网站色视频无遮挡免费观看| 精品人妻熟女毛片av久久网站| 后天国语完整版免费观看| 亚洲 欧美一区二区三区| 亚洲全国av大片| 国产男女超爽视频在线观看| 国产精品 欧美亚洲| 亚洲三区欧美一区| 女性生殖器流出的白浆| 欧美成人免费av一区二区三区 | 国产高清激情床上av| 又大又爽又粗| 午夜日韩欧美国产| 国产主播在线观看一区二区| 一区二区三区精品91| 我的亚洲天堂| 天天躁狠狠躁夜夜躁狠狠躁| 嫁个100分男人电影在线观看| 淫妇啪啪啪对白视频| ponron亚洲| 免费高清在线观看日韩| svipshipincom国产片| 麻豆乱淫一区二区| 欧美黑人欧美精品刺激| 欧美黄色淫秽网站| 丰满饥渴人妻一区二区三| 国产精品九九99| 精品久久久久久电影网| 欧美不卡视频在线免费观看 | 国产成人精品无人区| 夫妻午夜视频| 国产一区二区三区视频了| 90打野战视频偷拍视频| 色尼玛亚洲综合影院| 成人特级黄色片久久久久久久| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| 色老头精品视频在线观看| 1024香蕉在线观看| 夜夜夜夜夜久久久久| 亚洲国产精品合色在线| 日本黄色日本黄色录像| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品1区2区在线观看. | 国产日韩一区二区三区精品不卡| 91字幕亚洲| 午夜亚洲福利在线播放| 亚洲精品国产色婷婷电影| 女人被躁到高潮嗷嗷叫费观| 天堂中文最新版在线下载| 黑人猛操日本美女一级片| 精品一区二区三区视频在线观看免费 | 一级毛片精品| 美女国产高潮福利片在线看| 一进一出抽搐gif免费好疼 | 99国产精品免费福利视频| 天天添夜夜摸| 男人的好看免费观看在线视频 | 国产精品一区二区在线观看99| 国产精品久久视频播放| 国产成人免费观看mmmm| 亚洲情色 制服丝袜| 国产精品免费视频内射| 国产成人精品在线电影| 亚洲成国产人片在线观看| 国产亚洲av高清不卡| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 少妇粗大呻吟视频| 夜夜爽天天搞| 9191精品国产免费久久| 欧美老熟妇乱子伦牲交| 99久久综合精品五月天人人| 丝袜在线中文字幕| 国产成人精品无人区| 人妻丰满熟妇av一区二区三区 | tocl精华| 欧美乱色亚洲激情| 91麻豆精品激情在线观看国产 | 午夜日韩欧美国产| 美女高潮喷水抽搐中文字幕| 很黄的视频免费| 一级a爱视频在线免费观看| 久久久久精品人妻al黑| 亚洲avbb在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产精品偷伦视频观看了| 9191精品国产免费久久| 十八禁网站免费在线| 久久久国产成人精品二区 | 亚洲视频免费观看视频| 女人精品久久久久毛片| 少妇粗大呻吟视频| 又大又爽又粗| 日韩免费高清中文字幕av| 亚洲第一青青草原| 精品午夜福利视频在线观看一区| 18禁观看日本| 国产黄色免费在线视频| 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 国产精品 国内视频| 亚洲中文日韩欧美视频| 免费在线观看完整版高清| 亚洲一区高清亚洲精品| 中文亚洲av片在线观看爽 | 久热爱精品视频在线9| 色94色欧美一区二区| 男女床上黄色一级片免费看| 在线av久久热| 午夜视频精品福利| 久久99一区二区三区| 国精品久久久久久国模美| 99re在线观看精品视频| 两人在一起打扑克的视频| 国产亚洲一区二区精品| 亚洲aⅴ乱码一区二区在线播放 | 热re99久久精品国产66热6| 国产一区二区激情短视频| 女警被强在线播放| 精品国产一区二区三区久久久樱花| 日韩三级视频一区二区三区| 久久午夜亚洲精品久久| 99热网站在线观看| 高清欧美精品videossex| 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 9191精品国产免费久久| 亚洲片人在线观看| 成年人黄色毛片网站| 91在线观看av| 国产淫语在线视频| 大型av网站在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲午夜理论影院| 婷婷精品国产亚洲av在线 | av网站在线播放免费| 下体分泌物呈黄色| 午夜福利乱码中文字幕| 亚洲国产欧美一区二区综合| 日韩欧美在线二视频 | 麻豆av在线久日| 国产男女内射视频| 国产精品久久久久成人av| tube8黄色片| 满18在线观看网站| svipshipincom国产片| 欧美成人午夜精品| 一夜夜www| 久久国产亚洲av麻豆专区| 日韩欧美一区视频在线观看| 欧美成人免费av一区二区三区 | 在线观看免费日韩欧美大片| 黄色成人免费大全| 亚洲国产精品一区二区三区在线| 欧美乱色亚洲激情| 国产精品av久久久久免费| svipshipincom国产片| 一级黄色大片毛片| 国产精品自产拍在线观看55亚洲 | 搡老岳熟女国产| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 国产精品九九99| 久久午夜综合久久蜜桃| 丰满的人妻完整版| 欧美另类亚洲清纯唯美| 国产一区有黄有色的免费视频| 日本一区二区免费在线视频| 免费在线观看影片大全网站| 国产成人影院久久av| 757午夜福利合集在线观看| 露出奶头的视频| 人妻丰满熟妇av一区二区三区 | 久久久久视频综合| 国产野战对白在线观看| 亚洲人成伊人成综合网2020| 999久久久国产精品视频| 在线观看免费日韩欧美大片| 麻豆乱淫一区二区| 在线观看午夜福利视频| 变态另类成人亚洲欧美熟女 | 五月开心婷婷网| 国产精品一区二区免费欧美| 欧美午夜高清在线| 欧美精品av麻豆av| 亚洲欧美色中文字幕在线| 最新的欧美精品一区二区| 亚洲色图综合在线观看| 热re99久久国产66热| 看免费av毛片| 黄色 视频免费看| 露出奶头的视频| 一进一出抽搐动态| 色婷婷久久久亚洲欧美| 免费女性裸体啪啪无遮挡网站| 精品高清国产在线一区| 国产免费现黄频在线看| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻丝袜一区二区| 黄色丝袜av网址大全| 精品国产一区二区三区四区第35| 亚洲精品乱久久久久久| 日韩成人在线观看一区二区三区| 欧美黄色淫秽网站| 在线视频色国产色| 在线av久久热| tube8黄色片| 极品人妻少妇av视频| 日韩免费av在线播放| 国产色视频综合| 国产精品亚洲一级av第二区| 18禁裸乳无遮挡免费网站照片 | 国产野战对白在线观看| 国产亚洲欧美98| 黑人巨大精品欧美一区二区蜜桃| 精品一区二区三区四区五区乱码| 乱人伦中国视频| 免费观看a级毛片全部| 色老头精品视频在线观看| 人妻一区二区av| 狂野欧美激情性xxxx| 国产99白浆流出| 99国产综合亚洲精品| 日本vs欧美在线观看视频| 亚洲av日韩在线播放| www.精华液| 三上悠亚av全集在线观看| 又大又爽又粗| 成人永久免费在线观看视频| 丁香欧美五月| 在线观看午夜福利视频| 亚洲欧美日韩高清在线视频| 少妇粗大呻吟视频| 中文字幕另类日韩欧美亚洲嫩草| 久久久国产一区二区| 婷婷成人精品国产| 国产精品一区二区在线观看99| 免费在线观看日本一区| 嫩草影视91久久| 69av精品久久久久久| 亚洲情色 制服丝袜| 国产男靠女视频免费网站| 母亲3免费完整高清在线观看| 国产成人免费无遮挡视频| 午夜福利乱码中文字幕| 别揉我奶头~嗯~啊~动态视频| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| av电影中文网址| 国产有黄有色有爽视频| 老司机亚洲免费影院| 国产伦人伦偷精品视频| 手机成人av网站| 一夜夜www| 成人精品一区二区免费| 国产又色又爽无遮挡免费看| 视频区图区小说| 亚洲成人免费电影在线观看| 亚洲七黄色美女视频| 欧美 亚洲 国产 日韩一| 中文字幕另类日韩欧美亚洲嫩草| www日本在线高清视频| 男女之事视频高清在线观看| 欧美精品啪啪一区二区三区| 中亚洲国语对白在线视频| 国产乱人伦免费视频| 天天操日日干夜夜撸| 精品福利观看| 在线观看午夜福利视频| 不卡一级毛片| 三级毛片av免费| 人人妻,人人澡人人爽秒播| 两人在一起打扑克的视频| 午夜免费成人在线视频| 国产单亲对白刺激| 久久精品aⅴ一区二区三区四区| 国产精品 欧美亚洲| 国产欧美日韩一区二区三| 国产aⅴ精品一区二区三区波| 另类亚洲欧美激情| 亚洲情色 制服丝袜| 一二三四社区在线视频社区8| 夜夜躁狠狠躁天天躁| 在线天堂中文资源库| 一级毛片高清免费大全| 香蕉丝袜av| 欧美人与性动交α欧美软件| 国产精品久久久久久人妻精品电影| 一边摸一边抽搐一进一出视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲 国产 在线| 99久久精品国产亚洲精品| 色在线成人网| 校园春色视频在线观看| 一进一出好大好爽视频| 亚洲在线自拍视频| 女人久久www免费人成看片| 黑人操中国人逼视频| 五月开心婷婷网| 咕卡用的链子| 男女免费视频国产| 一级黄色大片毛片| 欧美人与性动交α欧美精品济南到| 少妇粗大呻吟视频| 最新在线观看一区二区三区| 午夜影院日韩av| 久久久国产一区二区| 另类亚洲欧美激情| av国产精品久久久久影院| 视频区欧美日本亚洲| netflix在线观看网站| 国产有黄有色有爽视频| 欧美精品亚洲一区二区| 久久九九热精品免费| 999久久久国产精品视频| 性少妇av在线| 日韩精品免费视频一区二区三区| 久久热在线av| 91成年电影在线观看| 亚洲专区字幕在线| 麻豆av在线久日| 身体一侧抽搐| 亚洲成a人片在线一区二区| 精品一区二区三卡| 午夜免费观看网址| 亚洲成国产人片在线观看| 大型黄色视频在线免费观看| 夜夜爽天天搞| 亚洲精华国产精华精| 日日夜夜操网爽| 黄片小视频在线播放| 欧美一级毛片孕妇| 身体一侧抽搐| 一边摸一边抽搐一进一小说 | 亚洲成人国产一区在线观看| 亚洲第一欧美日韩一区二区三区| 国产在线一区二区三区精| 欧美成狂野欧美在线观看| av免费在线观看网站| 亚洲国产精品一区二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 女人久久www免费人成看片| 人妻丰满熟妇av一区二区三区 | 中文欧美无线码| 九色亚洲精品在线播放| 亚洲av日韩在线播放| 国产在视频线精品| 成人亚洲精品一区在线观看| 国产99久久九九免费精品| 国产91精品成人一区二区三区| 午夜成年电影在线免费观看| 久久国产乱子伦精品免费另类| 99热网站在线观看| 午夜免费成人在线视频| 国产精品秋霞免费鲁丝片| 无遮挡黄片免费观看| av一本久久久久| 亚洲熟妇熟女久久| 久久精品亚洲熟妇少妇任你| 18禁裸乳无遮挡免费网站照片 | 日韩人妻精品一区2区三区| 久久久精品区二区三区| 亚洲欧美激情综合另类| 久久久精品免费免费高清| 脱女人内裤的视频| 亚洲精品中文字幕一二三四区| 久久香蕉精品热| 国产欧美日韩综合在线一区二区| 国产成人系列免费观看| 制服人妻中文乱码| 国产男女超爽视频在线观看| 国产精品免费视频内射| 国产国语露脸激情在线看| 欧美精品高潮呻吟av久久| 国产高清videossex| 国产精品.久久久| 又黄又爽又免费观看的视频| 香蕉国产在线看| 国产在线一区二区三区精| www.精华液| 亚洲精品国产精品久久久不卡| 欧美性长视频在线观看| 制服诱惑二区| 丁香欧美五月| 人妻丰满熟妇av一区二区三区 | 日日摸夜夜添夜夜添小说| 18禁美女被吸乳视频| 国产视频一区二区在线看| 久久这里只有精品19| 国产又色又爽无遮挡免费看| av国产精品久久久久影院| 日本a在线网址| 久久久久久免费高清国产稀缺| e午夜精品久久久久久久| 精品国产美女av久久久久小说| 一区二区三区国产精品乱码| 日韩人妻精品一区2区三区| 久久久久久久国产电影| 午夜久久久在线观看| 欧美激情极品国产一区二区三区| 亚洲一区中文字幕在线| 中文字幕人妻熟女乱码| 成年动漫av网址| 丁香六月欧美| 亚洲黑人精品在线| 久久人人97超碰香蕉20202| www.精华液| 久久久久久久午夜电影 | 国产又色又爽无遮挡免费看| 欧美最黄视频在线播放免费 | 91国产中文字幕| 亚洲成人免费av在线播放| xxx96com| 精品一区二区三区视频在线观看免费 | 男女免费视频国产| 免费在线观看影片大全网站| 大片电影免费在线观看免费| 十八禁网站免费在线| 热99国产精品久久久久久7| 亚洲色图 男人天堂 中文字幕| 国产一区有黄有色的免费视频| 淫妇啪啪啪对白视频| 51午夜福利影视在线观看| 老汉色∧v一级毛片| 亚洲综合色网址| 国产欧美亚洲国产| 亚洲午夜精品一区,二区,三区| 一区在线观看完整版| 精品一区二区三卡| 国产成人免费无遮挡视频| av超薄肉色丝袜交足视频| av线在线观看网站| 最新美女视频免费是黄的| 中文字幕另类日韩欧美亚洲嫩草| 超碰97精品在线观看| 精品午夜福利视频在线观看一区| 欧美激情极品国产一区二区三区| 午夜激情av网站| 精品久久久久久电影网| 欧美午夜高清在线| 91九色精品人成在线观看| 国产色视频综合| 国产真人三级小视频在线观看| 一级毛片女人18水好多| 在线看a的网站| 18禁观看日本| 伦理电影免费视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美亚洲日本最大视频资源| 国产精品免费大片| 午夜福利在线免费观看网站| 麻豆成人av在线观看| 欧美黄色淫秽网站| 久久久久精品人妻al黑| 国产精品国产高清国产av | 精品国产一区二区三区四区第35| 久久国产亚洲av麻豆专区| 久久精品成人免费网站| 极品少妇高潮喷水抽搐| 欧美久久黑人一区二区| 最近最新中文字幕大全免费视频| 亚洲 国产 在线| 久久影院123| 免费女性裸体啪啪无遮挡网站| 女人被躁到高潮嗷嗷叫费观| 男女高潮啪啪啪动态图| 91av网站免费观看| 老熟妇仑乱视频hdxx| 欧美国产精品一级二级三级| 精品国产一区二区三区久久久樱花| 国产一区二区三区视频了| 欧美中文综合在线视频| 精品国产一区二区三区久久久樱花| 亚洲精品中文字幕在线视频| 午夜久久久在线观看| 俄罗斯特黄特色一大片| 每晚都被弄得嗷嗷叫到高潮| 国内久久婷婷六月综合欲色啪| 中文字幕另类日韩欧美亚洲嫩草| 一进一出抽搐gif免费好疼 |