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

    長白山闊葉紅松林中影響喬木幼苗存活的關(guān)鍵因素

    2015-01-19 07:03:36王均偉趙秀海唐景毅
    生態(tài)學(xué)報(bào) 2015年19期
    關(guān)鍵詞:紅松林樣方負(fù)相關(guān)

    肖 翠, 劉 帥, 黃 珍, 樊 瑩, 王均偉, 趙秀海,*, 唐景毅

    1 北京林業(yè)大學(xué)省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100083 2 吉林省白河林業(yè)局, 吉林 133600

    長白山闊葉紅松林中影響喬木幼苗存活的關(guān)鍵因素

    肖 翠1, 劉 帥1, 黃 珍1, 樊 瑩1, 王均偉1, 趙秀海1,*, 唐景毅2

    1 北京林業(yè)大學(xué)省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100083 2 吉林省白河林業(yè)局, 吉林 133600

    應(yīng)用廣義線性混合模型,對長白山闊葉紅松林中影響幼苗存活的生物因素和非生物因素進(jìn)行分析。結(jié)果表明:(1)在群落水平上,幼苗存活率和生物因素中同種幼苗鄰居顯著負(fù)相關(guān),說明在闊葉紅松林群落中存在負(fù)密度制約效應(yīng)。(2)生物因子和非生物因子對不同年齡階段的幼苗存活率影響不同。對于1年生的幼苗,幼苗存活率與異種大樹鄰居呈顯著負(fù)相關(guān),與同種大樹鄰居呈顯著正相關(guān);對于2—3年生的幼苗,其存活率和同種幼苗鄰居、同種大樹鄰居均呈顯著負(fù)相關(guān),和非生物因子相關(guān)不顯著;對于4年生以上的幼苗,其存活率和土壤主成分分析的PC1(低的有機(jī)質(zhì)、速效鉀、速效氮等比較貧瘠的土壤)顯著負(fù)相關(guān)。(3)種子的傳播方式不同,幼苗存活率的影響因子也不同。對于風(fēng)傳播的物種,存活率與同種幼苗鄰居密度顯著正相關(guān)。對于重力傳播的物種,幼苗存活率與土壤PC3(高的全氮和速效氮,含氮較高的土壤)、異種大樹鄰居、草本密度呈負(fù)相關(guān),與林冠開闊度和草本蓋度成正相關(guān)。(4)對于不同的物種,影響幼苗存活的因素也不同。紫椴的幼苗存活率與土壤PC3、異種大樹鄰居、草本密度呈顯著負(fù)相關(guān),與土壤主成分PC2(高的有機(jī)質(zhì)和全氮等養(yǎng)分比較好的土壤)呈顯著正相關(guān)。

    長白山闊葉紅松林; 廣義線性混合模型; 生物因素; 非生物因素; 幼苗存活

    在樹木生活史周期中,幼苗階段是個(gè)體生長最為脆弱、對環(huán)境變化最為敏感的時(shí)期,也是個(gè)體數(shù)量變化最大的時(shí)期[1]。樹木幼苗的存活狀況間接影響森林群落的物種組成、結(jié)構(gòu)以及穩(wěn)定性。影響幼苗存活的因素很多,如光照、土壤水分、小尺度環(huán)境干擾、坡向、坡度、林冠蓋度、枯枝落葉層等各種環(huán)境因子[2-3],以及林分條件、物種組成、草食動(dòng)物啃食、菌類感染、幼苗與母樹之間的距離等生物因子[4-6,1]。然而,以往的研究多集中在對單個(gè)種群幼苗存活問題的探討上[7-8],對整個(gè)森林群落幼苗存活問題的研究還相對較少。

    近年來,隨著森林固定監(jiān)測大樣地的發(fā)展[9],長期連續(xù)定位監(jiān)測森林群落的動(dòng)態(tài)變化得以實(shí)現(xiàn)[10]。幼苗監(jiān)測作為森林樣地監(jiān)測的一個(gè)重要方面,越來越受到重視[11-14]。充分結(jié)合和綜合利用大樣地中的監(jiān)測數(shù)據(jù),能夠幫助了解幼苗的動(dòng)態(tài)變化規(guī)律以及影響其存活的原因[15-17]。長白山闊葉紅松林是典型的溫帶森林群落,具有豐富的物種多樣性,是北半球同緯度地區(qū)樹種最為豐富的森林類型。針對長白山闊葉紅松林的研究一直是我國溫帶森林生態(tài)系統(tǒng)研究的熱點(diǎn)[18-19],對群落幼苗存活規(guī)律的研究將進(jìn)一步完善我們對闊葉紅松林物種共存的認(rèn)識(shí)。有關(guān)闊葉紅松林中幼苗的研究已有一些[20-21],而綜合探討生物和非生物等多種因素對幼苗存活影響的研究還比較有限。

    綜上所述,本文以長白山闊葉紅松林永久監(jiān)測樣地為研究平臺(tái),以喬木幼苗為研究對象,研究影響幼苗存活的生物和非生物因素,擬回答以下幾個(gè)問題:(1)哪些因素影響闊葉紅松林幼苗的存活?(2)負(fù)密度制約現(xiàn)象在溫帶森林群落中是否普遍存在?(3)生物因素和非生物因素對不同年齡階段、不同傳播方式、不同物種的幼苗存活率影響是否一樣?(4)生物因素中的草本植物對幼苗存活是否有顯著影響?

    1 研究地概況與研究方法

    1.1 研究地概況

    研究區(qū)域位于長白山國家級(jí)自然保護(hù)區(qū)闊葉紅松林內(nèi)。本區(qū)地帶性氣候?yàn)闇貛駶櫟募撅L(fēng)氣候,冬季寒冷干燥,夏季溫暖濕潤。研究區(qū)域立地類型為平緩的坡地,屬典型的溫帶大陸性山地氣候[22]。該地區(qū)年均氣溫0.6 ℃,年均降水量707 mm,土壤類型為山地暗棕色森林土[23]。

    1.2 研究方法

    1.2.1 樣地設(shè)置

    在2007年,參照CTFS (Center for Tropical Forest Science)樣地建設(shè)標(biāo)準(zhǔn),選擇典型的原始闊葉紅松林,應(yīng)用全站儀建立了面積1 hm2(100 m×100 m)的固定監(jiān)測樣地。樣地中心坐標(biāo)為N42°20.211′,E128°05.705′,海拔為784 m。

    1.2.2 喬木調(diào)查

    記錄樣地內(nèi)所有胸徑(Diameter of Breast Height)大于等于1 cm的木本植物(包括喬木,灌木),記錄其種類,測量其胸徑和坐標(biāo),并對其定位及掛牌標(biāo)記[24]。

    1.2.3 幼苗及草本植物調(diào)查

    在樣地內(nèi)均勻設(shè)置了50個(gè)種子收集器,在每個(gè)種子收集器旁邊分別設(shè)置2個(gè)面積1 m×1 m的幼苗調(diào)查樣方,共計(jì)100個(gè)幼苗樣方,具體位置如圖1所示。將樣方內(nèi)胸徑小于1 cm的木本植物作為幼苗,調(diào)查其種類,并掛牌標(biāo)記。同時(shí)測量樣方內(nèi)草本植物種類、高度、密度和蓋度。其中,蓋度用投影法測量。即1 m×1 m樣方內(nèi)草本植物投影面積占整個(gè)樣方面積的比值。每年7—8月復(fù)查1次,迄今為止共完成了7次調(diào)查。

    圖1 幼苗樣方示意圖Fig.1 The sketch mapofseed trap and seedling grid空心方框代表草本和幼苗調(diào)查樣方,黑色方框代表種子收集器

    1.2.4 幼苗年齡的確定

    喬木幼苗的年齡的確定是根據(jù)長期定位監(jiān)測和主莖上的芽痕等因素綜合判斷的。1年生幼苗有子葉,無脫落芽痕;2—3年生幼苗存在1—2圈芽痕,每個(gè)芽痕之間有生長量。4年生以上的幼苗存在3圈以上芽痕。

    1.2.5 環(huán)境因子測量

    以樣地內(nèi)20 m×20 m為尺度測量環(huán)境因子,包括土壤有機(jī)質(zhì)、全氮、有效氮、有效磷、有效鉀、土壤pH值、林冠開闊度。具體調(diào)查方法是:在每個(gè)20 m×20 m的樣方內(nèi),隨機(jī)各取2個(gè)500 g土壤表層(0—20 cm)的土壤樣品,帶回分析。室內(nèi)分析依據(jù)《土壤農(nóng)業(yè)化學(xué)分析方法》的步驟進(jìn)行。兩個(gè)樣方的平均值為每個(gè)樣方測量結(jié)果的最終值。林冠開闊度指林地被樹冠垂直投影所覆蓋的比例。其測量采用林冠分析儀(WinSCANOPY,Canada)在每個(gè)樣方中距離地面1m處隨機(jī)拍攝2張半球狀照片,從半球狀照片中獲得。林冠開闊度也取平均值。在每個(gè)樣地所處地帶地勢較為平緩,沒有明顯的溝、山脊和谷地等特殊地形,故不考慮地形的影響。影響每個(gè)20 m×20 m樣方內(nèi)的幼苗存活的環(huán)境因子一樣。

    為了減少多元共線性現(xiàn)象和降低表述土壤因素的變量,參照Bai等對環(huán)境因子的處理方法,我們對土壤的有機(jī)質(zhì)、全氮、有效氮、有效磷、有效鉀5個(gè)代表土壤養(yǎng)分的化學(xué)變量進(jìn)行主成分分析。PC1代表低的有機(jī)質(zhì)、速效鉀、速效氮等比較貧瘠的土壤,PC2代表高的有機(jī)質(zhì)和全氮等養(yǎng)分比較好的土壤,PC3代表高的全氮和速效氮,含氮較高的土壤。具體結(jié)果見表1。

    表1 土壤養(yǎng)分的主成分分析

    1.2.6 數(shù)據(jù)分析方法

    選取2012年和2013年的幼苗調(diào)查數(shù)據(jù),應(yīng)用廣義線性混合模型分析幼苗存活與生物因子和非生物因子之間的關(guān)聯(lián)。廣義線性混合模型中包含固定效應(yīng)和隨機(jī)效應(yīng)兩部分自變量,將生物因子和非生物因子作為固定效應(yīng)并都經(jīng)過標(biāo)準(zhǔn)化處理,小樣方號(hào)和樹種被用作隨機(jī)效應(yīng)[10]。為了比較生物因子和非生物因子對幼苗存活的重要性,共建立了以下3個(gè)備選模型:(1)生物模型:模型中僅包含大樹效應(yīng)、幼苗效應(yīng)和草本植物效應(yīng)。大樹效應(yīng)是指以幼苗樣方為中心直徑10 m范圍內(nèi)的同種大樹胸高斷面積和異種大樹胸高斷面積[12]。不同研究選取的尺度不同,大樹效應(yīng)的影響也不同。長白山地區(qū)的研究大多數(shù)選擇20 m范圍,這主要是考慮林冠的影響。而本文試驗(yàn)地面積較小,面積僅有1 hm2,同時(shí)兼顧大樹的地上林冠和地下根系兩方面對幼苗的影響,故選取10 m范圍[12]。幼苗樣方10 m范圍內(nèi)的大樹均位于樣地內(nèi),即不存在邊緣效應(yīng)。幼苗效應(yīng)是指幼苗所在幼苗調(diào)查樣方內(nèi)的同種幼苗密度和異種幼苗密度,草本植物效應(yīng)是指幼苗調(diào)查樣方草本植物密度和蓋度;(2)非生物模型:模型中僅包含環(huán)境因子,即土壤PC1、土壤PC2、土壤PC3、土壤pH值和林冠開闊度。(3)生物+非生物模型。用Akaike′s Information Criterion (AIC)比較模型,不同的模型有不同的AIC值,AIC差值小于2的看作無差別[25]。將幼苗分成以下幾個(gè)類別進(jìn)行上述模型檢驗(yàn),(1)群落水平:包括所有喬木幼苗;(2)不同的年齡:1a、2—3a和大于等于4a的幼苗;(3)不同的種子擴(kuò)散方式:重力傳播和風(fēng)傳播。長白山地區(qū)靠動(dòng)物傳播種子的物種主要是紅松,但因?yàn)檎{(diào)查中2012年至2013年存活下來的紅松幼苗僅有5株,沒有達(dá)到統(tǒng)計(jì)分析的數(shù)量,故本文沒有分析動(dòng)物傳播方式;(4)物種水平:2012—2013年調(diào)查數(shù)據(jù)中存活數(shù)量最多的前5個(gè)樹種。為了比較各個(gè)參數(shù)對幼苗存活的貢獻(xiàn),計(jì)算每個(gè)參數(shù)的優(yōu)勢比率(Odds ratio),優(yōu)勢比率大于1(95%的置信區(qū)間)表明與幼苗存活正相關(guān),優(yōu)勢比率小于1(95%的置信區(qū)間)表明與幼苗存活負(fù)相關(guān)[20]。

    采用R軟件的‘lme4’[26]資源包進(jìn)行了幼苗存活關(guān)鍵因素的分析。

    2 結(jié)果

    2.1 群落水平的幼苗存活

    2012年共調(diào)查到1288株喬木幼苗,2013年調(diào)查到965棵喬木幼苗,從2012年至2013年仍然存活的幼苗有549棵,占2012年幼苗總量的42.6%。幼苗存活模型中用到的參數(shù)見表2。

    表2 幼苗存活模型中用到的參數(shù)

    文章中共用了3個(gè)模型:生物模型,非生物模型以及生物+非生物模型。通過表3中的ΔAIC值(ΔAIC=AIC-min(AIC))可以判斷模型的好壞。ΔAIC值越小的模型越好。

    表3 個(gè)體水平幼苗存活模型中AIC and ΔAIC值

    圖2 群落水平上幼苗存活的優(yōu)勢比率分析 Fig.2 Odds ratios of a tree seedling survival for community level圓圈代表每個(gè)參數(shù)的優(yōu)勢比率,實(shí)心點(diǎn)表示顯著相關(guān),空心點(diǎn)表示相關(guān)不顯著

    群落水平上,最好的模型為生物模型(表3),生物鄰居對幼苗的存活率影響更大。由圖2的優(yōu)勢比率可以看出,在群落水平上,同種幼苗密度和幼苗存活率呈負(fù)相關(guān)(odds ratio=-0.28,P=0.013),即同種幼苗密度越大的地方,幼苗存活率越低。

    2.2 不同年齡的幼苗存活

    不同年齡階段影響幼苗存活的最優(yōu)模型不同(圖3)。2—3年生幼苗最好的模型是生物模型,1年生幼苗和4年生以上幼苗最優(yōu)模型是生物+非生物模型。1年生幼苗存活率與同種大樹胸高斷面積呈顯著正相關(guān)(odds ratio=0.52,P=0.013),與異種大樹胸高斷面積呈顯著負(fù)相關(guān)(odds ratio=-0.37,P=0.010)。2—3年生的幼苗與同種大樹胸高斷面積呈顯著負(fù)相關(guān)(odds ratio=-0.36,P=0.005),與同種幼苗密度呈顯著負(fù)相關(guān)(odds ratio=-0.35,P=0.011)。對于大于等于4a的幼苗,其存活率與土壤主成分分析的PC1呈顯著負(fù)相關(guān)(odds ratio=-1.01,P=0.029)。

    圖3 不同年齡結(jié)構(gòu)幼苗存活的優(yōu)勢比率Fig.3 Odds ratios of tree seedling survival for three different age classes圓圈代表每個(gè)參數(shù)的優(yōu)勢比率,實(shí)心點(diǎn)表示顯著相關(guān),空心點(diǎn)表示相關(guān)不顯著

    2.3 不同種子傳播類型的幼苗存活

    表3可以看出,種子風(fēng)力傳播的幼苗存活最優(yōu)模型是生物模型,而重力傳播的幼苗存活的最優(yōu)模型是生物+非生物模型。

    從圖4中可以看出,種子風(fēng)力傳播的幼苗存活率與同種幼苗密度呈顯著負(fù)相關(guān)(odds ratio= -0.31,P=0.006)。種子重力傳播的幼苗存活率與林冠開闊度呈顯著正相關(guān)(odds ratio=0.49,P=0.049),與草本蓋度(odds ratio=0.89,P=0.038)呈正相關(guān),與土壤PC3呈顯著負(fù)相關(guān)(odds ratio=-0.66,P=0.006),與異種大樹胸高斷面積呈顯著負(fù)相關(guān)(odds ratio=-0.60,P=0.001),與草本密度呈顯著負(fù)相關(guān)(odds ratio=-0.98,P=0.010)。

    2.4 物種水平的幼苗存活

    在物種水平上,由表3可以看出,簇毛槭(Acerbarbinerve)、色木槭(Acermono)、水曲柳(Fraxinusmandschurica)的最優(yōu)模型均為生物模型,假色槭(Acerpseudo-sieboldianum)的最優(yōu)模型為非生物模型,紫椴(Tiliaamurensis)的最優(yōu)模型為生物+非生物模型。

    而對以上5個(gè)樣地中最常見的幼苗存活的所有生物和非生物因素進(jìn)行分析,結(jié)果顯示:簇毛槭、假色槭、色木槭、水曲柳4個(gè)物種與文章中所測的所有生物因素及非生物因素均無顯著相關(guān)性(圖5)。而只有常見種之一的紫椴與土壤主成分PC2、PC3、異種大樹胸高斷面積、草本密度有相關(guān)性(圖5)。紫椴幼苗存活率與土壤主成分PC2(odds ratio=0.99,P=0.010)呈正相關(guān)。與土壤主成分PC3(odds ratio=-0.78,P=0.003)、異種大樹鄰居(odds ratio=-0.60,P=0.004)、草本密度(odds ratio=-1.11,P=0.007)均呈負(fù)相關(guān)。

    圖4 不同種子傳播方式幼苗存活的優(yōu)勢比率Fig.4 Odds ratios of tree seedling survival for two dispersalmode groups圓圈代表每個(gè)參數(shù)的優(yōu)勢比率,實(shí)心點(diǎn)表示顯著相關(guān),空心點(diǎn)表示相關(guān)不顯著

    圖5 物種水平幼苗存活的優(yōu)勢比率Fig.5 Odds ratios of seedling survival of two common species 圓圈代表每個(gè)參數(shù)的優(yōu)勢比率,實(shí)心點(diǎn)表示顯著相關(guān),空心點(diǎn)表示相關(guān)不顯著

    3 結(jié)論與討論

    為了探討影響溫帶闊葉紅松林幼苗存活的因素,從群落水平、不同年齡階段、不同種子傳播方式和物種水平4個(gè)方面檢驗(yàn)了生物因素和非生物環(huán)境對幼苗存活率的影響。

    3.1 影響幼苗存活的因素是多樣的

    幼苗存活受生物鄰居和非生物環(huán)境共同影響[20],但兩者對不同狀態(tài)下的幼苗影響的程度不同。本文的研究中群落水平上幼苗的存活受生物鄰居的影響要大于環(huán)境因子的影響(表3)。同種幼苗密度與幼苗存活呈顯著的負(fù)相關(guān)(圖2),許多熱帶、亞熱帶的研究也得到了相同的結(jié)論[27-29]。隨著幼苗的生長,影響幼苗存活的因子也在發(fā)生變化。4年生以下的幼苗存活受生物因子影響較大(圖3),而4年生及以上幼苗存活則受環(huán)境因子影響較大(圖3)。即隨著幼苗年齡的增長,生物鄰居對幼苗的影響程度逐漸變小,生境對其生長的影響程度反而增大。這可能由于幼苗對周圍環(huán)境中與生物鄰居之間競爭力增強(qiáng),對周圍環(huán)境中真菌感染[30],食草動(dòng)物和病原菌的抗性增強(qiáng),而不斷生長過程中,需要更多的養(yǎng)分、水分等,反而更依賴于生境環(huán)境。種子風(fēng)力傳播的幼苗存活更易受生物因子的影響(圖4),而重力傳播的物種則受到生物和非生物因素的雙重作用(圖4),具體原因還有待進(jìn)一步探討。不同物種的幼苗存活的影響因子明顯不同,紫椴幼苗同時(shí)受生物和非生物多種因素顯著影響,地上部分和地下部分同時(shí)對紫椴的生長起作用[31]。而水曲柳幼苗則基本不受任何因素影響,這與物種本身的生理特性有關(guān),水曲柳有更強(qiáng)的適應(yīng)能力,生長迅速,成林快,適生于東北東部山地的中生、潮濕和濕生立地[32]。

    3.2 負(fù)密度制約效應(yīng)在溫帶森林中被進(jìn)一步驗(yàn)證

    本研究發(fā)現(xiàn)群落水平總體幼苗存活率(圖2)、2—3年生幼苗存活率(圖3)和種子風(fēng)力傳播幼苗存活率(圖4)都與同種幼苗密度或同種大樹胸高斷面積呈顯著負(fù)相關(guān),說明在該森林群落中存在負(fù)密度制約效應(yīng)。這與白雪嬌等在長白山的研究、祝燕對負(fù)密度制約的研究結(jié)果相一致[17,20,33]。負(fù)密度制約假說主要描述由于資源競爭、有害生物侵害(比如病原微生物、食草動(dòng)物捕食)等,同種個(gè)體之間發(fā)生的相互損害行為;它主要強(qiáng)調(diào)同種個(gè)體之間的相互作用,解釋自然群落物種共存的機(jī)理。負(fù)密度制約假說在多種森林群落中得到驗(yàn)證,Carson等在2008綜合分析了53項(xiàng)有關(guān)負(fù)密度制約的檢驗(yàn),其中81.2%的研究在熱帶森林,3.8%的研究在亞熱帶森林,15%的研究在溫帶森林[34]。盡管生態(tài)學(xué)家們?nèi)匀粦岩韶?fù)密度制約效應(yīng)調(diào)節(jié)群落物種共存的重要性,但是目前的研究還沒有發(fā)現(xiàn)否定負(fù)密度制約假說的充分證據(jù)。本文的研究驗(yàn)證了負(fù)密度制約效應(yīng)在溫帶森林中存在。

    3.3 草本效應(yīng)影響幼苗的存活率

    以往對于大樣地中草本的研究主要集中在多樣性和季節(jié)動(dòng)態(tài)、形成格局及影響因素等方面[35-37],也有文章分析草本格局的影響因素中提到喬木層通過影響光照和降水及枯落物的含量而間接影響草本層的多樣性格局[37]。而至今在分析影響喬木幼苗存活率的因素中沒有提及草本植物的影響。本文將幼苗樣方中的草本密度和草本蓋度引入幼苗存活率的影響因素中來。分析結(jié)果顯示,草本植物對某些特定的幼苗存活率有影響,比如草本密度與紫椴幼苗存活率呈負(fù)相關(guān)(圖5),草本密度與重力傳播物種的幼苗存活率也存在負(fù)相關(guān)關(guān)系(圖4,b)。因?yàn)楸狙芯康臄?shù)據(jù)量有限,草本植物對某些幼苗存活率的影響沒有檢測出來。有待以后加大數(shù)據(jù)量,更為準(zhǔn)確的監(jiān)測草本植物對幼苗存活率的影響程度。

    [1] Wright S J, Muller-Landau H C, Calderon O, Hernandez A. Annual and spatial variation in seedfall and seedling recruitment in a neotropical forest. Ecology, 2005, 86(4): 848-860.

    [2] Mason W L, Edwards C, Hale S E. Survival and early seedling growth of conifers with different shade tolerance in a Sitka spruce spacing trial and relationship to understory light climate. Silva Fennica, 2004, 38(4): 357-370.

    [3] Poorter L, Rose S A. Light-dependent changes in the relationship between seed mass and seedling traits: a meta-analysis for rain forest tree species. Oecologia, 2005, 142(3): 378-387.

    [4] Silvera K, Skillman J B, Dalling J W. Seed germination, seedling growth and habitat partitioning in two morphotypes of the tropical pioneer tree Tremamicrantha in a seasonal forest in Panama. Journal of Tropical Ecology, 2003, 19(1): 27-34.

    [5] 彭閃江, 黃忠良, 彭少麟, 歐陽學(xué)軍, 徐國良. 植物天然更新過程中種子和幼苗死亡的影響因素. 廣西植物, 2004, 24(2): 113-121.

    [6] Condit R, Ashton P, Bunyavejchewin S, Dattaraja HS, Davies S, Esufali S, Ewango C, Foster R, Gunatilleke I, Gunatilleke CV, Hall P, Harms KE, Hart T, Hernandez C, Hubbell S, Itoh A, Kiratiprayoon S, Lafrankie J, de Lao SL, Makana JR, Noor MN, Kassim AR, Russo S, Sukumar R, Samper C, Suresh HS, Tan S, Thomas S, Valencia R, Vallejo M, Villa G, Zillio T. The importance of demographic niches to tree diversity.Science, 2006, 313(5783):98-101.

    [7] 劉慶. 林窗大小和位置對麗江云杉自然更新幼苗存活和生長的影響. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2004, 10(3): 281-285.

    [8] 盛海燕, 李偉成, 葛瀅. 明黨參幼苗存活與生長對光照強(qiáng)度的響應(yīng). 應(yīng)用生態(tài)學(xué)報(bào), 2006, 17(5): 783-788.

    [9] 王斌, 黃俞淞, 李先琨, 向悟生, 丁濤, 黃甫昭, 陸樹華, 韓文衡, 文淑均, 何蘭軍. 弄崗北熱帶喀斯特季節(jié)性雨林15ha監(jiān)測樣地的樹種組成與空間分布. 生物多樣性, 2014, 22(2): 141-156.

    [10] 賓粵, 葉萬輝, 曹洪麟, 黃忠良, 練琚愉. 鼎湖山南亞熱帶常綠闊葉林20公頃樣地幼苗的分布. 生物多樣性, 2011, 19(2): 127-133.

    [11] Lambers J H R, Clark J S. Effects of dispersal, shrubs, and density-dependent mortality on seed and seedling distributions in temperate forests. Canadian Journal of Forest Research, 2003, 33(5): 783-795.

    [12] Comita L S, Hubbell S P. Local neighborhood and species′shade tolerance influence survival in a diverse seedling bank. Ecology, 2009, 90(2): 328-334.

    [13] Queenborough S A, Burslem D F R P, Garwood N C, Valencia R. Taxonomic scale-dependence of habitat niche partitioning and bioticneighbourhood on survival of tropical tree seedlings. Proceedings of theRoyal Society B:Biological Sciences, 2009, 276(1767): 4197-4205.

    [14] Streng D R, Glitzenstein J S, Harcombe P A. Woody seedling dynamics in an east texas floodplain forest. Ecology Monograph, 1989, 59: 177-204.

    [15] 鞏合德, 楊國平, 魯志云, 劉玉洪, 曹敏. 哀牢山常綠闊葉林喬木樹種的幼苗組成及時(shí)空分布特征. 生物多樣性, 2011, 19(2): 151-157.

    [16] 汪殷華, 米湘成, 陳聲文, 李銘紅, 于明堅(jiān). 古田山常綠闊葉林主要樹種2002—2007年間更新動(dòng)態(tài). 生物多樣性, 2011, 19(2): 178-189.

    [17] Zhang J, Hao Z Q, Sun I F, Song B, Ye J, Li B H, Wang X G. Density dependence on tree survival in an old-growth temperate forest in northeastern China. Annals of Forest Science, 2009, 66(2): 204-204.

    [18] Yuan Z Q, Gazol A, Wang X G, Li B H, Lin F, Ye J, Bai X J, Hao Z Q. Scale specific determinants of tree diversity in an old growth temperate forest in China. Basic and Applied Ecology, 2011, 12(6): 488-495.

    [19] Wang X G, Wiegand T,Wolf A, Howe R, Davies S J, Hao Z Q. Spatial patterns of tree species richness in two temperate forests. Journal of Ecology, 2011, 99(6): 1382-1393.

    [20] Bai X J, Queenborough S A, Wang X G, Zhang J, Li B H,Yuan Z Q, Xing D L,Lin F, Ye J, Hao Z Q. Effects of local biotic neighbors and habitat heterogeneity on tree and shrub seedling survival in an old-growth temperate forest.Oecologia, 2012, 170: 755-765.

    [21] 張健, 李步杭, 白雪嬌, 原作強(qiáng), 王緒高, 葉吉, 郝占慶. 長白山闊葉紅松林喬木樹種幼苗組成及其年際動(dòng)態(tài). 生物多樣性, 2009, 17(4): 385-396.

    [22] 李文華, 鄧?yán)っ? 李飛. 長白山主要生態(tài)系統(tǒng)生物量生產(chǎn)量的研究. 森林生態(tài)系統(tǒng)研究, 1981, (2): 34-50.

    [23] 許廣山, 丁桂芳, 張玉華. 長白山北坡主要森林土壤有機(jī)質(zhì)及其特征的初步研究. 森林生態(tài)系統(tǒng)研究, 1980, (1): 215-220.

    [24] Condit R. Tropical Forest Census Plots: Methods and Results from BarroColorado Island, Panama and a Comparison with Other Plots. Berlin:Springer, 1998.

    [25] Burnham K P, Anderson D R. Model Selection and Multi-modelInference: APractical Information-theoretic Approach. New York:Springer, 2002.

    [26] Bates D, Maechler M, Bolker B. lme4: Linear mixed-effects models using S4 classes. R Package Version: 0. 999375-37edn, 2010.

    [27] Metz M R, Sousa W P, Valencia R. Widespread density-dependent seedling mortality promotes species coexistence in a highly diverse Amazonian rain forest. Ecology, 2010, 91(12): 3675-3685.

    [28] Queenborough S A,Burslem D F R P, Garwood NC, Valencia R. Neighborhood and community interactions determine the spatial pattern of tropical tree seedling survival. Ecology, 2007, 88(9): 2248-2258.

    [29] Webb C O, Gilbert G S, Donoghue M J. Phylodiversity-dependent seedling mortality, size structure, and disease in a bornean rain forest. Ecology, 2006, 87(7): S123-S131.

    [30] Masaki T, Nakashizuka T. Seedling demography of Swidacontroversa: effect of light and distance to conspecifics. Ecology, 2002, 83(12): 3497-3507.

    [31] 王軍邦, 王政權(quán), 胡秉民, 牛錚, 王長耀. 不同栽植方式下紫椴幼苗生物量分配及資源利用分析. 植物生態(tài)學(xué)報(bào), 2002, 26(6): 677-683.

    [32] 荊濤, 馬萬里, Kujansuu J, 羅菊春, 王廣發(fā), 孫波. 水曲柳萌芽更新的研究. 北京林業(yè)大學(xué)學(xué)報(bào), 2002, 24(4): 12-15.

    [33] 祝燕, 米湘成, 馬克平. 植物群落物種共存機(jī)制: 負(fù)密度制約假說. 生物多樣性, 2009, 17(6): 594-604.

    [34] Carson W P, Anderson J, Leigh E G, Schnitzer S A. Challengesassociated with testing and falsifying the Janzen-Connell hypothesis: a review and critique// Carson WP, Stefan AS. Tropical Forest Community Ecology.Oxford, UK: Wiley-Blackwell,2008:210-241.

    [35] 夏富才, 潘春芳, 趙秀海, 何海燕, 周海城. 長白山原始闊葉紅松林林下草本植物多樣性格局及其影響因素. 西北植物學(xué)報(bào), 2012, 32(2): 370-376.

    [36] 吳剛, 尹若波, 周永斌, 郝占慶, 梁秀英. 長白山紅松闊葉林林隙動(dòng)態(tài)變化對早春草本植物的影響. 生態(tài)學(xué)報(bào), 1999, 19(5): 659-663.

    [37] 郝占慶, 趙士洞, 陶大立. 長白山北坡闊葉紅松林草本植物物種多樣性及其季節(jié)動(dòng)態(tài). 生物多樣性, 1994, 2(3): 125-132.

    Effects of biotic and abiotic factors on tree seedling survival in a broad-leaved

    Korean pine (Pinuskoraiensis) mixed forest on Changbai Mountain

    XIAO Cui1, LIU Shuai1, HUANG Zhen1, FAN Ying1, WANG Junwei1, ZHAO Xiuhai1,*, TANG Jingyi2

    1TheKeyLaboratoryforSilvicultureandConservationofMinistryofEducation,BeijingForestryUniversity,Beijing100083,China2JilinProvincialBaiheForestryBureau,Jilin133600,China

    Seedling recruitment plays a key role in determining species composition and diversity. Exploring patterns of seedling dynamics can therefore provide crucial insights into the mechanisms that affect seedling recruitment. The broad-leaved Korean pine (Pinuskoraiensis) mixed forest is a temperate forest typical of northeast China. In 2007, we established a 1 hm2plot in a broad-leaved Korean pine mixed forest on Changbai Mountain, northeast China. We identified, tagged, and mapped each tree with diameter at breast height ≥1 cm. At the same time, we set up 100 seedling quadrats of 1m × 1m which were censused every July. In total, we collected data on 1288 adult trees and 965 seedlings in the 1 hm2plot between 2012 and 2013, and examined the importance of biotic and abiotic factors on seedling survival using a generalized mixed linear model. Biotic factors included presence of conspecific seedling neighbors and heterospecific seedling neighbors, presence of conspecific adult neighbors and heterospecific adult neighbors, and herb density and herb cover. Abiotic factors included soil organic matter, available soil phosphorus, available soil potassium, available soil nitrogen, total soil nitrogen, and canopy openness. Conspecific neighbors had a significant negative effect on seedling survival at the community level, indicating significant negative density-dependence effects. The effects of biotic and abiotic factors on seedlings differed with seedling age. Conspecific adult neighbors and canopy openness had significant positive effects on survival of, 1-year seedlings while heterospecific adult neighbors had a significant negative effect. In contrast, for 2-to 3-year-old seedlings, conspecific seedling neighbors and conspecific adult neighbors showed a significant negative effect on seedling survival, but abiotic factors had no significant effect. Moreover, for ≥4-year-old seedlings, soil PC axis 1 (associated with low organic matter available K, available N, and total N) showed a significant negative effect on survival. Factors that drove seedling survival varied with species dispersal-mode. For wind-dispersed species, conspecific seedling neighbors had a significant positive effect on seedling survival. For gravity-dispersed species, soil PC axis 3 (associated with high available N and total N, heterospecific adult neighbors, and herb density) had a significant negative effect on survival, whereas canopy openness showed a significant positive effect. The effects of biotic neighbors and habitat differed among species. For instance, soil PC axis 3, heterospecific adult neighbors, and herb density had significant negative effects on seedling survival ofTiliaamurensis, while soil PC axis 2 (associated with high organic matter and total N) showed a significant positive effect. Finally, herbs might influence seedling survival. Both herb density and herb cover impacted seedling survival by modifying local habitat. However, there was no significant effect of herbs on survival due to limited data. We conclude that both local biotic and abiotic factors play important roles in determining seedling recruitment in this temperate forest.

    broad-leaved Korean pine (Pinuskoraiensis) mixed forest on Changbai Mountain; generalized linear mixed-effect model; biotic factors; environmental factors; seedling survival

    國家“十二五”科技支撐課題(2012BAD22B0203, 2012BAC01B03); 北京市共建項(xiàng)目專項(xiàng)(天然林生物多樣性保護(hù)技術(shù)與楊樹抗逆機(jī)理研究2012)

    2014-03-04; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2014-12-04

    10.5846/stxb201403040371

    *通訊作者Corresponding author.E-mail: zhaoxh@bjfu.edu.cn

    肖翠, 劉帥, 黃珍, 樊瑩, 王均偉, 趙秀海, 唐景毅.長白山闊葉紅松林中影響喬木幼苗存活的關(guān)鍵因素.生態(tài)學(xué)報(bào),2015,35(19):6557-6565.

    Xiao C, Liu S, Huang Z, Fan Y, Wang J W, Zhao X H, Tang J Y.Effects of biotic and abiotic factors on tree seedling survival in a broad-leaved Korean pine (Pinuskoraiensis) mixed forest on Changbai Mountain.Acta Ecologica Sinica,2015,35(19):6557-6565.

    猜你喜歡
    紅松林樣方負(fù)相關(guān)
    藥王谷多種因素對漏蘆、野韭生長的影響
    昆明金殿國家森林公園林業(yè)有害生物調(diào)查研究
    沒有紅松的紅松林
    N-末端腦鈉肽前體與糖尿病及糖尿病相關(guān)并發(fā)癥呈負(fù)相關(guān)
    典型野生刺梨(Rosa roxburghii Tratt.)灌草叢植物多樣性研究
    紅松林采伐更新的確定方式
    云南九龍山自然保護(hù)區(qū)水青岡群落物種多樣性及群落特征
    東折棱河楓樺紅松林與柞樹紅松林結(jié)構(gòu)特征對比分析
    更 正
    翻譯心理與文本質(zhì)量的相關(guān)性探析
    考試周刊(2016年63期)2016-08-15 14:33:26
    日本 欧美在线| 欧美日韩黄片免| 亚洲av成人一区二区三| 中文在线观看免费www的网站| 久久久精品大字幕| 在线国产一区二区在线| 亚洲欧美日韩高清在线视频| 午夜福利成人在线免费观看| 成年女人毛片免费观看观看9| 91麻豆av在线| 精品久久久久久久久久免费视频| 亚洲av五月六月丁香网| 亚洲精品在线观看二区| 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频| 日韩欧美在线乱码| av中文乱码字幕在线| 国产欧美日韩精品亚洲av| 最近最新免费中文字幕在线| 国产高清videossex| 国产高清三级在线| 日本熟妇午夜| 大型黄色视频在线免费观看| 婷婷六月久久综合丁香| 精品久久久久久久末码| 色尼玛亚洲综合影院| 婷婷亚洲欧美| 国产高清激情床上av| 成人午夜高清在线视频| 亚洲第一电影网av| 国产精品久久久久久久电影 | 黄频高清免费视频| 成人鲁丝片一二三区免费| 9191精品国产免费久久| 欧美av亚洲av综合av国产av| 亚洲av熟女| 中出人妻视频一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 五月伊人婷婷丁香| 91久久精品国产一区二区成人 | 最新美女视频免费是黄的| 黄色日韩在线| 国产成年人精品一区二区| 一本久久中文字幕| 天堂av国产一区二区熟女人妻| 日韩欧美 国产精品| 99在线人妻在线中文字幕| 一本一本综合久久| 麻豆一二三区av精品| 午夜福利高清视频| 极品教师在线免费播放| 97超级碰碰碰精品色视频在线观看| 成人精品一区二区免费| 久久人人精品亚洲av| 舔av片在线| 18禁黄网站禁片午夜丰满| 国产av麻豆久久久久久久| 在线a可以看的网站| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 99久久综合精品五月天人人| av天堂中文字幕网| 18禁国产床啪视频网站| 久99久视频精品免费| 国产精品电影一区二区三区| 女人高潮潮喷娇喘18禁视频| 两个人视频免费观看高清| 99久久综合精品五月天人人| 国产精品98久久久久久宅男小说| 岛国视频午夜一区免费看| 黄片大片在线免费观看| 丁香六月欧美| 日韩免费av在线播放| 亚洲精品国产精品久久久不卡| 欧美日韩一级在线毛片| 无限看片的www在线观看| 日韩欧美精品v在线| 国产黄色小视频在线观看| 性欧美人与动物交配| 看黄色毛片网站| 18禁观看日本| 午夜亚洲福利在线播放| 国内久久婷婷六月综合欲色啪| av天堂中文字幕网| 久久香蕉国产精品| 少妇的丰满在线观看| 99国产精品一区二区三区| 国产一区在线观看成人免费| 黄片大片在线免费观看| 丁香六月欧美| 午夜福利18| 18禁美女被吸乳视频| 人妻久久中文字幕网| 亚洲av五月六月丁香网| www.999成人在线观看| 麻豆国产av国片精品| 国产精品野战在线观看| 精品人妻1区二区| 99国产综合亚洲精品| 性欧美人与动物交配| 日韩精品中文字幕看吧| 88av欧美| 狂野欧美白嫩少妇大欣赏| 黑人操中国人逼视频| 国产成人系列免费观看| 俺也久久电影网| 国内精品久久久久精免费| 男人舔女人的私密视频| 中文字幕熟女人妻在线| 女人被狂操c到高潮| 日韩国内少妇激情av| 午夜福利欧美成人| 久久久色成人| 一个人看视频在线观看www免费 | 国产精华一区二区三区| 国产亚洲精品av在线| 久久久国产欧美日韩av| 中国美女看黄片| 亚洲av成人精品一区久久| 悠悠久久av| or卡值多少钱| 精品99又大又爽又粗少妇毛片 | 午夜激情欧美在线| 欧美日韩精品网址| 精品乱码久久久久久99久播| 两人在一起打扑克的视频| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲精品不卡| 在线观看免费午夜福利视频| 激情在线观看视频在线高清| 男女床上黄色一级片免费看| 99久久精品国产亚洲精品| 一个人看视频在线观看www免费 | 桃红色精品国产亚洲av| 免费在线观看日本一区| 国产一区在线观看成人免费| 一级作爱视频免费观看| 少妇熟女aⅴ在线视频| 亚洲美女视频黄频| 无人区码免费观看不卡| 日韩国内少妇激情av| 男女做爰动态图高潮gif福利片| 嫁个100分男人电影在线观看| 成人一区二区视频在线观看| 午夜福利欧美成人| 麻豆成人午夜福利视频| 欧美乱色亚洲激情| 一进一出好大好爽视频| 国产精品美女特级片免费视频播放器 | 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品乱码久久久v下载方式 | 久久久成人免费电影| 制服丝袜大香蕉在线| 成熟少妇高潮喷水视频| 国内精品一区二区在线观看| 99久久精品热视频| 久99久视频精品免费| 一个人看视频在线观看www免费 | 真人一进一出gif抽搐免费| 中文字幕高清在线视频| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 亚洲成人免费电影在线观看| 精品免费久久久久久久清纯| 免费在线观看日本一区| 日韩欧美精品v在线| 日韩高清综合在线| 老司机深夜福利视频在线观看| 久久久成人免费电影| 日韩有码中文字幕| 黄色日韩在线| 香蕉丝袜av| 精品久久久久久久久久久久久| 我的老师免费观看完整版| 久久中文字幕一级| 国产人伦9x9x在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 最近在线观看免费完整版| 亚洲一区二区三区不卡视频| 色吧在线观看| 欧美黄色淫秽网站| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 精品国内亚洲2022精品成人| 最新中文字幕久久久久 | 麻豆成人av在线观看| 亚洲熟妇熟女久久| 熟女少妇亚洲综合色aaa.| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 丰满人妻熟妇乱又伦精品不卡| 国内少妇人妻偷人精品xxx网站 | 久久久久久九九精品二区国产| 18禁美女被吸乳视频| 亚洲人成网站高清观看| 男女下面进入的视频免费午夜| 老司机在亚洲福利影院| 亚洲精品粉嫩美女一区| 高清在线国产一区| 亚洲自拍偷在线| 男插女下体视频免费在线播放| 国产亚洲精品av在线| 国产欧美日韩一区二区三| 久久这里只有精品19| 日韩国内少妇激情av| 男女床上黄色一级片免费看| 亚洲精品乱码久久久v下载方式 | 国产高清有码在线观看视频| 久久久久久久久中文| 精品国产三级普通话版| 亚洲一区高清亚洲精品| 日本免费一区二区三区高清不卡| 午夜福利在线观看吧| 久久久久久久久中文| 日韩欧美一区二区三区在线观看| 亚洲一区高清亚洲精品| 男人舔女人的私密视频| 一本一本综合久久| 成人高潮视频无遮挡免费网站| 欧美午夜高清在线| 91麻豆av在线| 久久久久久九九精品二区国产| 久久久久久久精品吃奶| 中文资源天堂在线| 麻豆国产97在线/欧美| 亚洲男人的天堂狠狠| ponron亚洲| 老鸭窝网址在线观看| 少妇的逼水好多| 婷婷六月久久综合丁香| 99久久综合精品五月天人人| 亚洲精品在线观看二区| 国产真实乱freesex| 欧美日韩瑟瑟在线播放| 噜噜噜噜噜久久久久久91| 午夜福利在线在线| 亚洲成av人片在线播放无| 男人舔奶头视频| 欧美zozozo另类| 精品久久久久久久人妻蜜臀av| 嫩草影视91久久| 国模一区二区三区四区视频 | 欧美av亚洲av综合av国产av| 中文字幕精品亚洲无线码一区| 国产亚洲av高清不卡| 香蕉丝袜av| 啦啦啦观看免费观看视频高清| 午夜福利在线在线| 69av精品久久久久久| www.熟女人妻精品国产| 精品电影一区二区在线| 欧美日韩乱码在线| 国产成人精品久久二区二区免费| 日本精品一区二区三区蜜桃| 天堂影院成人在线观看| 男人舔女人下体高潮全视频| av片东京热男人的天堂| 日本免费一区二区三区高清不卡| 午夜精品在线福利| 熟女电影av网| а√天堂www在线а√下载| 欧美最黄视频在线播放免费| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 午夜福利在线在线| 黄色片一级片一级黄色片| 99视频精品全部免费 在线 | 麻豆久久精品国产亚洲av| 91久久精品国产一区二区成人 | 国产精品99久久99久久久不卡| 黄色视频,在线免费观看| 午夜精品一区二区三区免费看| 国内揄拍国产精品人妻在线| 中出人妻视频一区二区| 久久久国产欧美日韩av| 97超级碰碰碰精品色视频在线观看| 身体一侧抽搐| 久久九九热精品免费| 成年版毛片免费区| 欧美成人性av电影在线观看| 国产午夜精品论理片| 亚洲精品美女久久av网站| 精品久久久久久久久久久久久| 欧美av亚洲av综合av国产av| 国产69精品久久久久777片 | 国产精品久久电影中文字幕| 99re在线观看精品视频| 又爽又黄无遮挡网站| 久久久久久人人人人人| 99在线视频只有这里精品首页| 麻豆成人av在线观看| 最近最新中文字幕大全免费视频| 免费在线观看视频国产中文字幕亚洲| 床上黄色一级片| 最近最新中文字幕大全电影3| 亚洲国产看品久久| 老司机午夜福利在线观看视频| 久久欧美精品欧美久久欧美| 99热6这里只有精品| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 成年版毛片免费区| 免费搜索国产男女视频| 亚洲色图av天堂| 日本免费a在线| 国产精品野战在线观看| 亚洲精品国产精品久久久不卡| 久久热在线av| 美女高潮喷水抽搐中文字幕| 国产aⅴ精品一区二区三区波| 国产一区二区在线观看日韩 | 亚洲国产中文字幕在线视频| 一a级毛片在线观看| 首页视频小说图片口味搜索| 日本与韩国留学比较| 麻豆成人午夜福利视频| 国产高清视频在线播放一区| 在线观看免费午夜福利视频| 国内少妇人妻偷人精品xxx网站 | 在线观看一区二区三区| 精品日产1卡2卡| 亚洲精品一区av在线观看| 国产美女午夜福利| 亚洲av成人不卡在线观看播放网| 国产精品美女特级片免费视频播放器 | 久久中文字幕人妻熟女| 韩国av一区二区三区四区| 国产精品九九99| 夜夜躁狠狠躁天天躁| 在线免费观看的www视频| 99久久成人亚洲精品观看| 欧美在线黄色| 亚洲欧美精品综合久久99| 人人妻人人澡欧美一区二区| 国产69精品久久久久777片 | 最近最新中文字幕大全电影3| 久久精品国产亚洲av香蕉五月| 狂野欧美激情性xxxx| 亚洲乱码一区二区免费版| 无人区码免费观看不卡| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 亚洲成人中文字幕在线播放| 人妻久久中文字幕网| 一级a爱片免费观看的视频| 欧美3d第一页| 18禁黄网站禁片免费观看直播| 精品一区二区三区四区五区乱码| 小说图片视频综合网站| 欧美日韩乱码在线| 日日夜夜操网爽| 999久久久国产精品视频| 香蕉丝袜av| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久中文看片网| 老司机午夜福利在线观看视频| 午夜两性在线视频| 听说在线观看完整版免费高清| 日本在线视频免费播放| 99国产极品粉嫩在线观看| 免费搜索国产男女视频| 天堂av国产一区二区熟女人妻| 男女下面进入的视频免费午夜| 午夜激情欧美在线| 国产视频一区二区在线看| 99国产综合亚洲精品| 天堂网av新在线| 欧美在线黄色| 女警被强在线播放| 曰老女人黄片| 亚洲av电影在线进入| 级片在线观看| 在线视频色国产色| 欧美色视频一区免费| 久久久久国产精品人妻aⅴ院| 国产av不卡久久| 亚洲国产精品久久男人天堂| 久久久久久人人人人人| 久久久久国内视频| 亚洲va日本ⅴa欧美va伊人久久| 床上黄色一级片| 伦理电影免费视频| 亚洲欧美日韩无卡精品| 熟女少妇亚洲综合色aaa.| 精品国产三级普通话版| 九色国产91popny在线| 国产高清视频在线观看网站| 美女被艹到高潮喷水动态| 国产精品九九99| 亚洲中文日韩欧美视频| 淫秽高清视频在线观看| 日韩国内少妇激情av| 欧美zozozo另类| 国产精品1区2区在线观看.| 欧美xxxx黑人xx丫x性爽| 亚洲专区国产一区二区| 久久九九热精品免费| 悠悠久久av| 欧美乱妇无乱码| 国产精品99久久99久久久不卡| 天天添夜夜摸| 日日夜夜操网爽| 长腿黑丝高跟| 两个人看的免费小视频| 国产精品久久久久久精品电影| 久久草成人影院| 亚洲真实伦在线观看| 性色avwww在线观看| 欧美乱色亚洲激情| 麻豆国产97在线/欧美| 中文字幕最新亚洲高清| 美女免费视频网站| 很黄的视频免费| 国产99白浆流出| 色精品久久人妻99蜜桃| 国产精品亚洲美女久久久| 女生性感内裤真人,穿戴方法视频| 亚洲天堂国产精品一区在线| 国产精品久久久av美女十八| 精品一区二区三区视频在线 | 久久香蕉精品热| 日本a在线网址| netflix在线观看网站| 99热只有精品国产| 黄色丝袜av网址大全| 97人妻精品一区二区三区麻豆| 少妇丰满av| a在线观看视频网站| 国产精品99久久久久久久久| 亚洲无线在线观看| 我要搜黄色片| 中文资源天堂在线| 亚洲av电影在线进入| 麻豆一二三区av精品| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 国产真人三级小视频在线观看| xxxwww97欧美| 国内少妇人妻偷人精品xxx网站 | 免费看光身美女| 叶爱在线成人免费视频播放| 国产探花在线观看一区二区| 一进一出好大好爽视频| 国产激情久久老熟女| 亚洲av成人不卡在线观看播放网| 国产欧美日韩一区二区精品| 国产精品久久久久久久电影 | 国产精品精品国产色婷婷| 国产单亲对白刺激| 嫩草影院入口| 人妻丰满熟妇av一区二区三区| 国产成+人综合+亚洲专区| 天堂√8在线中文| 伦理电影免费视频| 91老司机精品| 精品人妻1区二区| 操出白浆在线播放| 成人特级黄色片久久久久久久| а√天堂www在线а√下载| 国产熟女xx| 国产一区二区三区在线臀色熟女| 成人欧美大片| 哪里可以看免费的av片| 熟女人妻精品中文字幕| 1000部很黄的大片| 成年人黄色毛片网站| 欧美日韩综合久久久久久 | 国产一区二区激情短视频| 国产激情欧美一区二区| 极品教师在线免费播放| av天堂中文字幕网| 狂野欧美白嫩少妇大欣赏| 91av网一区二区| 国产aⅴ精品一区二区三区波| 国产亚洲欧美在线一区二区| 老司机深夜福利视频在线观看| 这个男人来自地球电影免费观看| 久久精品人妻少妇| 校园春色视频在线观看| 天堂√8在线中文| 国产精品99久久久久久久久| 一本久久中文字幕| 女人高潮潮喷娇喘18禁视频| 国产av一区在线观看免费| 老汉色av国产亚洲站长工具| 又黄又粗又硬又大视频| 高清毛片免费观看视频网站| 伊人久久大香线蕉亚洲五| 精品日产1卡2卡| 亚洲18禁久久av| 亚洲无线观看免费| 国模一区二区三区四区视频 | 成熟少妇高潮喷水视频| 免费高清视频大片| 国产免费男女视频| 日本一二三区视频观看| 2021天堂中文幕一二区在线观| 热99re8久久精品国产| 国产乱人伦免费视频| 三级国产精品欧美在线观看 | av黄色大香蕉| 免费av毛片视频| 亚洲欧洲精品一区二区精品久久久| 国产精品av视频在线免费观看| www.熟女人妻精品国产| 99热只有精品国产| 久久久久久久久免费视频了| 国产免费av片在线观看野外av| 熟女少妇亚洲综合色aaa.| 国产91精品成人一区二区三区| 99久久综合精品五月天人人| 国产一区在线观看成人免费| 色综合欧美亚洲国产小说| 欧美av亚洲av综合av国产av| 欧美一级毛片孕妇| 国产真人三级小视频在线观看| 欧美另类亚洲清纯唯美| 1024手机看黄色片| 1000部很黄的大片| 可以在线观看的亚洲视频| 国产爱豆传媒在线观看| 美女扒开内裤让男人捅视频| 国产高清视频在线观看网站| 亚洲欧美日韩东京热| 午夜精品在线福利| av国产免费在线观看| 怎么达到女性高潮| 亚洲七黄色美女视频| 黄色 视频免费看| 精品久久蜜臀av无| 一二三四社区在线视频社区8| 免费看光身美女| 久9热在线精品视频| 日韩 欧美 亚洲 中文字幕| 日韩大尺度精品在线看网址| 欧美激情在线99| 99久久久亚洲精品蜜臀av| 99久久综合精品五月天人人| 国产黄a三级三级三级人| 激情在线观看视频在线高清| 18禁黄网站禁片午夜丰满| 久久这里只有精品中国| 天堂av国产一区二区熟女人妻| 母亲3免费完整高清在线观看| 成人国产一区最新在线观看| 国内久久婷婷六月综合欲色啪| 一个人免费在线观看的高清视频| 日本成人三级电影网站| aaaaa片日本免费| 亚洲,欧美精品.| 一本综合久久免费| 亚洲欧美日韩高清专用| avwww免费| 看黄色毛片网站| 十八禁人妻一区二区| 亚洲成人精品中文字幕电影| 人妻丰满熟妇av一区二区三区| 成人午夜高清在线视频| 大型黄色视频在线免费观看| 伦理电影免费视频| 少妇裸体淫交视频免费看高清| 999精品在线视频| 欧美另类亚洲清纯唯美| a级毛片a级免费在线| 熟女电影av网| 中出人妻视频一区二区| 国产三级黄色录像| 国语自产精品视频在线第100页| x7x7x7水蜜桃| 在线免费观看的www视频| 亚洲va日本ⅴa欧美va伊人久久| 成人av在线播放网站| 亚洲美女视频黄频| 国产亚洲精品一区二区www| 久久性视频一级片| 少妇丰满av| 色综合婷婷激情| 久久久久国内视频| 美女大奶头视频| av黄色大香蕉| 成年女人永久免费观看视频| 欧美绝顶高潮抽搐喷水| 操出白浆在线播放| 国产精品一及| 精品久久久久久久人妻蜜臀av| 一个人免费在线观看电影 | tocl精华| 男女下面进入的视频免费午夜| 国产亚洲精品综合一区在线观看| 免费观看的影片在线观看| 色噜噜av男人的天堂激情| 欧美高清成人免费视频www| 给我免费播放毛片高清在线观看| 在线观看一区二区三区| 久久精品国产综合久久久| 久久午夜亚洲精品久久| 精品国产三级普通话版| 国产亚洲精品一区二区www| 精品一区二区三区视频在线 | 桃红色精品国产亚洲av| 日韩 欧美 亚洲 中文字幕| 国产亚洲精品综合一区在线观看| 香蕉国产在线看| 久久精品国产亚洲av香蕉五月| 国产精品亚洲美女久久久| 精品熟女少妇八av免费久了| 国产私拍福利视频在线观看| 国产免费男女视频| 狂野欧美白嫩少妇大欣赏| 亚洲精品色激情综合| 激情在线观看视频在线高清| e午夜精品久久久久久久|