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

    食餌庇護所對斑塊環(huán)境下Leslie-Gower捕食系統(tǒng)的影響

    2018-06-07 03:07:30劉志廣
    生態(tài)學報 2018年8期
    關(guān)鍵詞:庇護所食餌時間尺度

    劉志廣

    河南大學應(yīng)用數(shù)學研究所數(shù)學與統(tǒng)計學院, 開封 475004

    捕食是許多生態(tài)系統(tǒng)常見而又非常重要的種間作用關(guān)系。數(shù)學模型是理解和分析捕食者-食餌系統(tǒng)動態(tài)行為的一個重要且有效的工具,這其中最早也最著名的是Lotka-Volterra 模型[1]。這個模型假設(shè)捕食者對食餌的捕食效率(即功能反應(yīng))隨食餌種群密度的增加而成直線增加,但是實驗數(shù)據(jù)表明這個假設(shè)對很多種群不太合理[2]。根據(jù)物種類型和它們所處的環(huán)境,Holling提出3種不同類型的功能反應(yīng)[3]。這些功能反應(yīng)函數(shù)都假設(shè)捕食者的捕食效率只依賴食餌的數(shù)量,因此被稱為食餌依賴的功能反應(yīng)。但越來越多的研究發(fā)現(xiàn),有些物種捕食者的捕食效率除了受到食餌數(shù)量的影響,還會受到捕食者數(shù)量的影響,因此提出各種捕食者依賴的功能反應(yīng)函數(shù),例如率依賴[4]和捕食者間相互干擾[5-6]的功能反應(yīng)函數(shù)。

    另一類有趣的捕食者-食餌模型是考慮捕食者的增長函數(shù)與捕食函數(shù)不同。假設(shè)捕食者的增長項依賴捕食者和食餌數(shù)量的比值,這就是所謂的Leslie-Gower捕食者-食餌模型[7]。通過數(shù)值模擬[8]或嚴格的數(shù)學證明[9-10],可以知道這個系統(tǒng)唯一的正平衡點是全局穩(wěn)定的?;贚eslie-Gower 框架,很多更加實際的具體因素被引入到這個模型或者其變種形式中,如其他形式的功能反應(yīng)[11-12]、擴散[13]、時滯效應(yīng)[14]、收獲效應(yīng)[15]、Allee效應(yīng)[16]和庇護所效應(yīng)[17]等。

    庇護所對捕食者-食餌系統(tǒng)的影響是數(shù)學生態(tài)學和理論生態(tài)學領(lǐng)域備受關(guān)注的課題。大部分研究關(guān)注食餌庇護所對系統(tǒng)穩(wěn)定性的影響。研究表明,食餌庇護所可以穩(wěn)定共存平衡點且防止食餌種群滅絕[18]。也有一些研究揭示在一定條件下可以破壞共存平衡點的穩(wěn)定性[19]。另一方面,時間和空間庇護所對雜食系統(tǒng)物種共存的影響也越來越受關(guān)注[20-22]。以往研究大多利用基于Lotka-Volterra框架的捕食模型。研究方法往往通過修改功能反應(yīng)函數(shù)來隱含的考慮食餌庇護所效應(yīng)。環(huán)境異質(zhì)性可以為食餌種群提供空間庇護所。斑塊環(huán)境下的種群相互作用模型是研究環(huán)境異質(zhì)性的一個有效方法[23-26]。此模型由兩部分組成,一部分用于描述局部的物種相互作用,另一部分用于描述斑塊間的遷移。受此啟發(fā),本研究提出一個兩斑塊環(huán)境的Leslie-Gower模型。其中一個斑塊表示開放的生境,另一個斑塊表示食餌的庇護所,從而就可以顯式地研究食餌庇護所的影響?;诖?本研究重點關(guān)注食餌庇護所對系統(tǒng)穩(wěn)定性和捕食者食餌密度變化的影響。

    1 捕食者-食餌的完全模型

    假設(shè)有兩個斑塊,斑塊 1 上面生活著捕食者和食餌。x1(t)和y(t)表示在斑塊1上t時刻的食餌和捕食者的種群密度,可以用以下微分方程來描述:

    (1)

    這就是所謂的Leslie-Gower捕食者-食餌模型[7- 11]。此模型假設(shè)捕食者的環(huán)境容納量與食餌數(shù)量成比例。它強調(diào)捕食者和食餌的增長率都存在一個上限的實際情況,而這是Lotka-Volterra模型沒有注意到的[10]。這個上限只有在適合的條件下才能達到:對于捕食者來說,是在單位捕食者占有的食餌數(shù)量最大時;對于食餌來說,是在單位食餌占有的捕食者最少時。其中,r1,a1,b1,r2和a2是模型參數(shù),都假設(shè)為正值。這些參數(shù)的生物學意義為:r1表示斑塊1上食餌的增長率,b1表示斑塊1上食餌的種內(nèi)競爭強度,a1是最大捕食率,r2表示斑塊1捕食者的增長率,a2和a1有類似的意義。

    另外,假設(shè)斑塊2為庇護所斑塊,只有食餌才可以遷移到上面,而捕食者不能遷移到上面。且假設(shè)斑塊2上面食餌的增長率為負,這是考慮到庇護所斑塊上面缺少資源和尋找配偶的機會[27]。假設(shè)遷移和生物相互作用發(fā)生在兩個不同的時間尺度上。遷移發(fā)生在快時間尺度,而種群的增長、死亡和捕食發(fā)生在慢時間尺度。本研究僅考慮非密度依賴的遷移方式。根據(jù)以上假設(shè),可以建立一個顯式含有庇護所的兩斑塊捕食者-食餌模型如下:

    (2)

    式中,參數(shù)ε是一個非常小且無量綱的正數(shù),它表示快慢時間尺度的比率。τ是快時間尺度。參數(shù)β表示食餌從斑塊1到斑塊2的遷移率,α的意義相反。參數(shù)d表示斑塊2上食餌的增長率,d>0表示在庇護所斑塊上食餌只能死亡。

    2 綜合模型

    基于已有的理論[23],可以將以上的完全系統(tǒng)簡化成只有兩個變量的綜合模型。這個方法的思想就是當有快慢兩個時間尺度的時候,認為在慢時間尺度系統(tǒng)的每一個時刻,快時間尺度系統(tǒng)都處于平衡態(tài)。因此可以將快時間尺度系統(tǒng)的平衡點代入到慢時間尺度的系統(tǒng)中,從而得到簡化的綜合系統(tǒng)。具體操作的方法是:第一步,忽略方程的慢時間尺度部分,從而研究在快時間尺度下的遷移模型。先來計算快時間尺度系統(tǒng)的平衡點。選取兩斑塊上物種的總密度作為綜合變量,即

    x1(τ)+x2(τ)=x(τ),y(τ)=y(τ)

    (3)

    在快時間尺度上,x(τ)和y(τ)是常數(shù)。遷移可以改變兩個斑塊上食餌的比例,但不能改變它們的總密度。對方程(2),令ε=0,快的時間系統(tǒng)的平衡點為

    (4)

    (I) 把模型(2)的前兩個方程相加,得到

    (5)

    (6)

    式中t=ετ是慢時間尺度。

    3 綜合模型正平衡點的穩(wěn)定性

    通過簡單計算,可以得到系統(tǒng)(6)唯一的內(nèi)部平衡點:

    (7)

    (8)

    因此,關(guān)于內(nèi)部平衡點穩(wěn)定性的主要結(jié)果表述如下:

    證明:首先構(gòu)造李雅普諾夫函數(shù):

    (9)

    顯然,V(x,y)當x>0,y>0是有定義且連續(xù)。這個函數(shù)滿足:

    (10)

    (11)

    沿著系統(tǒng)(6)的解,計算V(x,y)的導數(shù),并結(jié)合等式(8),得到

    4 食餌庇護所的影響

    4.1 食餌庇護所對捕食者密度的影響

    (12)

    4.2 食餌庇護所對食餌密度的影響

    改寫x*的表達式為如下形式:

    (13)

    通過計算x*關(guān)于u的導數(shù):

    (14)

    5 結(jié)果與討論

    在空間生態(tài)學的研究中,常見的數(shù)學模型可分為反應(yīng)擴散方程、基于個體模型、集合種群模型、斑塊模型和網(wǎng)格模型[24]。網(wǎng)格模型常用于考慮空間結(jié)構(gòu)和局部擴散對系統(tǒng)時空動態(tài)的影響。而斑塊模型常用于研究空間異質(zhì)性和全局擴散對生態(tài)系統(tǒng)時空動態(tài)的影響。本研究就是利用斑塊模型來考慮庇護所斑塊和普通斑塊間的異質(zhì)性對種群動態(tài)影響的。另外以往關(guān)于斑塊環(huán)境下捕食系統(tǒng)的研究大多假設(shè)遷移和局域相互作用發(fā)生在同一時間尺度,很少研究假設(shè)遷移發(fā)生在快時間尺度[25-26]。但很多現(xiàn)實情況是遷移和種群的增長、死亡及捕食發(fā)生在不同的時間尺度。比如,在水生環(huán)境下,浮游動物每天都表現(xiàn)出一種垂直遷移。兩個主要因子可以解釋這種現(xiàn)象:光和食物。在白天,光線增加了被魚類捕食的風險,這導致很多物種向下遷移到光線暗的地方從而減少捕食風險。晚上,這些物種在向上遷移去消費那些白天生產(chǎn)出的浮游植物[26]。在這種情況下,遷移相對于種群的增長、死亡和捕食來說是快時間尺度的。因為遷移每天都發(fā)生,而種群的增長、死亡和捕食發(fā)生在星期或者月的時間尺度上的。因此,假設(shè)遷移和種群過程不在同一個時間尺度往往更符合實際[23,25- 27]。另外,兩個時間尺度的好處在于可以利用奇異攝動理論將系統(tǒng)簡化[25],這樣數(shù)學上更方便分析。否則,包括兩個斑塊和遷移的模型非常復(fù)雜,很難得到一般的分析結(jié)果。

    大部分關(guān)于空間庇護所對捕食系統(tǒng)動態(tài)影響的研究都基于Lotka-Volterra模型框架,有的研究采用隱含空間庇護所的模型[19],有的研究采用顯式含有空間庇護所的模型[27]。本研究基于Leslie-Gower捕食模型系統(tǒng),而基于這個模型揭示食餌庇護所對系統(tǒng)動態(tài)的研究相對較少。因此本研究對全面認識和理解庇護所對不同類型生態(tài)系統(tǒng)動態(tài)的影響是非常有幫助的。最近,也有研究基于Leslie-Gower捕食系統(tǒng),采用空間隱含的方法研究食餌庇護所對系統(tǒng)動態(tài)的影響[17]。研究揭示食餌庇護所對系統(tǒng)的續(xù)存沒有影響,但是增加了食餌平衡密度;在一定條件下,捕食者平衡密度與庇護所的數(shù)量呈現(xiàn)單峰曲線。與之不同,本研究采用顯式含有空間庇護所的方法,這樣更能體現(xiàn)空間及空間異質(zhì)性所起到的作用。研究表明捕食者的平衡種群密度隨著躲進庇護所食餌數(shù)量增加而降低。這是因為被保護的食餌數(shù)量越多,捕食者的食物就越少,當然數(shù)量會越低。研究還發(fā)現(xiàn)在一定條件下,隨著躲進庇護所食餌數(shù)量的增加,食餌的平衡密度先增加后降低。這是因為我們假設(shè)庇護所上食餌不能生存。當躲進庇護所食餌數(shù)量比較小時,庇護所的保護作用大于捕食效應(yīng),所以食餌數(shù)量增加。但當被保護的食餌數(shù)量達到一定程度時,因庇護所上沒有食物,此時庇護所對躲在它上面的食餌的不利影響大于保護作用,因此食餌數(shù)量下降。這也就是說食餌種群存在一種利用庇護所的最佳策略。不能一味的躲進庇護所里面不出來,這實質(zhì)上是捕食和庇護作用的權(quán)衡。顯式含有和隱含空間庇護所的數(shù)學模型所得結(jié)論不一致,這意味著在研究庇護所對捕食系統(tǒng)種群動態(tài)影響時,空間結(jié)構(gòu)可能起著重要作用。到底哪種理論預(yù)測的結(jié)果更合理,還需結(jié)合以后的實驗工作來檢驗。

    本研究假設(shè)斑塊間食餌的遷移是一種非密度依賴的方式。它代表的是一種逃逸行為,但這里給出的是被動的擴散而已,并沒有體現(xiàn)出隨捕食壓力加重的逃逸行為。假設(shè)食餌遷移是一種密度依賴的方式會更加合理。比如假設(shè)食餌的遷移率依賴所在斑塊的捕食者密度[25]、食餌密度[23]或者同時依賴于捕食者和食餌密度[26]。另外,本研究為了分析方便,假設(shè)只有兩個空間斑塊,擴展到多斑塊的情形也是值得思考的問題。此時,庇護所的空間位置、大小及空間結(jié)構(gòu)對系統(tǒng)動態(tài)行為都會有影響[28]。

    參考文獻(References):

    [1] Lotka A J. Elements of Physical Biology. Baltimore: Williams and Wilkins Co., Inc., 1924.

    [2] Begon M, Townsend C R, Harper J L. Ecology: From Individuals to Ecosystems. 4th ed. Oxford: Blackwell Publishing, 2005.

    [3] Holling C S. Some characteristics of simple types of predation and parasitism. The Canadian Entomologist, 1959, 91(7): 385- 398.

    [4] Arditi R, Ginzburg L R. Coupling in predator-prey dynamics: ratio-dependence. Journal of Theoretical Biology, 1989, 139(3): 311- 326.

    [5] Beddington J R. Mutual interference between parasites or predators and its effect on searching efficiency. Journal of Animal Ecology, 1975, 44(1): 331- 340.

    [6] DeAngelis D L, Goldstein R A, O′Neill R V. A model for tropic interaction. Ecology, 1975, 56(4): 881- 892.

    [7] Leslie P H, Gower J C. The properties of a stochastic model for the predator-prey type of interaction between two species. Biometrika, 1960, 47(3/4): 219- 234.

    [8] Pielou EC. Mathematical Ecology. New York: John Wiley & Sons, 1977.

    [9] Hsu SB, Huang TW. Global stability for a class of predator-prey systems. SIAM Journal on Applied Mathematics, 1995, 55(3): 763- 783.

    [10] Korobeinikov A. A Lyapunov function for Leslie-Gower predator-prey models. Applied Mathematics Letters, 2001, 14(6): 697- 699.

    [11] Aziz-Alaoui M A, Okiye M D. Boundedness and global stability for a predator-prey model with modified Leslie-Gower and Holling-type II schemes. Applied Mathematics Letters, 2003, 16(7): 1069- 1075.

    [12] Huang J C, Ruan S G, Song J. Bifurcations in a predator-prey system of Leslie type with generalized Holling type III functional response. Journal of Differential Equations, 2014, 257(6): 1721- 1752.

    [13] Yang W S. Global asymptotical stability and persistent property for a diffusive predator-prey system with modified Leslie-Gower functional response. Nonlinear Analysis: Real World Applications, 2013, 14(3): 1323- 1330.

    [14] Nindjin A F, Aziz-Alaoui M A, Cadivel M. Analysis of a predator-prey model with modified Leslie-Gower and Holling-type II schemes with time delay. Nonlinear Analysis: Real World Applications, 2006, 7(5): 1104- 1118.

    [15] Gupta RP, Chandra P. Bifurcation analysis of modified Leslie-Gower predator-prey model with Michaelis-Menten type prey harvesting. Journal of Mathematical Analysis and Applications, 2013, 398(1): 278- 295.

    [16] Aguirre P, González-Olivares E, Sáez E. Three limit cycles in a Leslie-Gower predator-prey model with additive Allee effect. SIAM Journal on Applied Mathematics, 2009, 69(5): 1244- 1262.

    [17] Chen F D, Chen L J, Xie X D. On a Leslie-Gower predator-prey model incorporating a prey refuge. Nonlinear Analysis: Real World Applications, 2009, 10(5): 2905- 2908.

    [18] González-Olivares E, Ramos-Jiliberto R. Dynamic consequences of prey refuges in a simple model system: more prey, fewer predators and enhanced stability. Ecological Modelling, 2003, 166(1/2): 135- 146.

    [19] Ma Z H, Li W L, Zhao Y, Wang W T, Zhang H, Li Z Z. Effects of prey refuges on a predator-prey model with a class of functional responses: the role of refuges. Mathematical Biosciences, 2009, 218(2): 73- 79.

    [20] Finke D L, Denno R F. Spatial refuge from intraguild predation: implications for prey suppression and trophic cascades. Oecologia, 2006, 149(2): 265- 275.

    [21] Amarasekare P. Coexistence of intraguild predators and prey in resource-rich environments. Ecology, 2008, 89(10): 2786- 2797.

    [22] Liu Z G, Zhang F P. Species coexistence of communities with intraguild predation: the role of refuges used by the resource and the intraguild prey. Biosystems, 2003, 114(1): 25- 30.

    [23] Michalski J, Poggiale JC, Arditi R, Auger PM. Macroscopic dynamic effects of migrations in patchy predator-prey systems. Journal of Theoretical Biology, 1997, 185(4): 459- 474.

    [24] Briggs C J, Hoopes M F. Stabilizing effects in spatial parasitoid-host and predator-prey models: a review. Theoretical Population Biology, 2004, 65(3): 299- 315.

    [25] El Abdllaoui A, Auger P, Kooi B W, de laParra R B, Mchich R. Effects of density-dependent migrations on stability of a two-patch predator-prey model. Mathematical Biosciences, 2007, 210(1): 335- 354.

    [26] Mchich R, Auger P, Poggiale JC. Effect of predator density dependent dispersal of prey on stability of a predator-prey system. Mathematical Biosciences, 2007, 206(2): 343- 356.

    [27] Ma Z H, Wang S F, Li W D, Li Z Z. The effect of prey refuge in a patchy predator-prey system. Mathematical Biosciences, 2013, 243(1): 126- 130.

    [28] Mistro D C, Rodrigues L A D, Varriale M C. The role of spatial refuges in coupled map lattice model for host-parasitoid systems. Bulletin of Mathematical Biology, 2009, 71(8): 1934- 1953.

    猜你喜歡
    庇護所食餌時間尺度
    捕食-食餌系統(tǒng)在離散斑塊環(huán)境下強迫波的唯一性
    時間尺度上非完整系統(tǒng)的Noether準對稱性與守恒量
    “給我”庇護所
    一類具有修正的Leslie-Gower項的捕食-食餌模型的正解
    時間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    力學學報(2021年10期)2021-12-02 02:32:04
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    具有兩個食餌趨化項的一個Ronsenzwing-MacArthur捕食食餌模型的全局分歧
    How Archaeologist Discovered “Treasures of Apollo” Buried Below Ancient Greek Ruins
    一類帶有交叉擴散的捕食-食餌模型的正解
    大連市暴雨多時間尺度研究分析
    tocl精华| 动漫黄色视频在线观看| www国产在线视频色| 观看美女的网站| 国产又色又爽无遮挡免费看| 亚洲中文日韩欧美视频| 亚洲欧美一区二区三区黑人| 十八禁网站免费在线| 国语自产精品视频在线第100页| 99精品在免费线老司机午夜| 日本在线视频免费播放| 日韩欧美免费精品| 99re在线观看精品视频| 国产伦在线观看视频一区| 好看av亚洲va欧美ⅴa在| 亚洲欧美日韩东京热| 成人特级黄色片久久久久久久| 日韩欧美国产一区二区入口| 国产三级在线视频| 人妻夜夜爽99麻豆av| 午夜福利在线观看免费完整高清在 | 国产美女午夜福利| 中文字幕精品亚洲无线码一区| 狂野欧美激情性xxxx| 国产aⅴ精品一区二区三区波| 久久国产精品人妻蜜桃| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看亚洲国产| 日本三级黄在线观看| 美女被艹到高潮喷水动态| 在线十欧美十亚洲十日本专区| 日本成人三级电影网站| 老司机在亚洲福利影院| 午夜福利高清视频| av国产免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 岛国视频午夜一区免费看| 99在线人妻在线中文字幕| cao死你这个sao货| 狂野欧美白嫩少妇大欣赏| 国产视频内射| 高清在线国产一区| 亚洲欧美日韩高清在线视频| 少妇的逼水好多| 小说图片视频综合网站| 亚洲国产高清在线一区二区三| 亚洲中文字幕一区二区三区有码在线看 | 18禁黄网站禁片免费观看直播| 久久这里只有精品中国| 亚洲18禁久久av| 欧美日韩中文字幕国产精品一区二区三区| 99国产综合亚洲精品| 久久久水蜜桃国产精品网| 免费看a级黄色片| 国产精品1区2区在线观看.| 国产视频一区二区在线看| 熟妇人妻久久中文字幕3abv| 99久久99久久久精品蜜桃| 九九在线视频观看精品| 757午夜福利合集在线观看| 国产视频一区二区在线看| 国产精品一区二区三区四区久久| 99re在线观看精品视频| 中文亚洲av片在线观看爽| 深夜精品福利| 欧美zozozo另类| 黄色片一级片一级黄色片| 国产激情久久老熟女| 精品久久久久久久人妻蜜臀av| 国产精品99久久99久久久不卡| 免费看a级黄色片| 无遮挡黄片免费观看| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 男人和女人高潮做爰伦理| 欧美成人免费av一区二区三区| 在线观看舔阴道视频| h日本视频在线播放| av天堂在线播放| 丁香欧美五月| 久久这里只有精品中国| 老司机午夜十八禁免费视频| 欧美黑人巨大hd| 亚洲九九香蕉| 日韩av在线大香蕉| 国模一区二区三区四区视频 | 亚洲欧美日韩高清在线视频| 色综合婷婷激情| 天堂影院成人在线观看| 欧美最黄视频在线播放免费| 一个人看的www免费观看视频| 久久久国产欧美日韩av| 久久精品国产亚洲av香蕉五月| 男女床上黄色一级片免费看| 露出奶头的视频| 国产黄色小视频在线观看| 成人性生交大片免费视频hd| 少妇的丰满在线观看| 一级毛片高清免费大全| 亚洲国产高清在线一区二区三| 国产成人系列免费观看| 色综合站精品国产| 一边摸一边抽搐一进一小说| 色视频www国产| 中文字幕久久专区| 老司机在亚洲福利影院| 2021天堂中文幕一二区在线观| 男女床上黄色一级片免费看| 日本免费a在线| 99热6这里只有精品| 久久精品人妻少妇| av天堂中文字幕网| 久久精品国产清高在天天线| 亚洲激情在线av| 露出奶头的视频| 亚洲av中文字字幕乱码综合| 久久久国产成人免费| 日本 欧美在线| 美女午夜性视频免费| 欧美成人免费av一区二区三区| 18禁观看日本| 亚洲熟妇熟女久久| 亚洲成人免费电影在线观看| 丁香六月欧美| 成人无遮挡网站| 久久精品aⅴ一区二区三区四区| 高潮久久久久久久久久久不卡| 无限看片的www在线观看| 亚洲人与动物交配视频| 每晚都被弄得嗷嗷叫到高潮| 90打野战视频偷拍视频| 国产亚洲精品综合一区在线观看| 香蕉av资源在线| 亚洲成人免费电影在线观看| 老熟妇仑乱视频hdxx| 日本成人三级电影网站| 国产欧美日韩精品一区二区| 国产精品 国内视频| 日本 av在线| 91av网一区二区| 999精品在线视频| 国产激情久久老熟女| 美女黄网站色视频| 亚洲专区字幕在线| 国内久久婷婷六月综合欲色啪| 国产精品av视频在线免费观看| 少妇熟女aⅴ在线视频| 亚洲欧美日韩高清在线视频| 中出人妻视频一区二区| 噜噜噜噜噜久久久久久91| 搡老妇女老女人老熟妇| 亚洲人与动物交配视频| 日韩欧美国产在线观看| 熟女少妇亚洲综合色aaa.| 丁香六月欧美| 身体一侧抽搐| 日本免费一区二区三区高清不卡| 亚洲国产高清在线一区二区三| 小说图片视频综合网站| 亚洲成人久久性| 国产精品久久久人人做人人爽| 丰满的人妻完整版| 少妇裸体淫交视频免费看高清| 精品欧美国产一区二区三| 国内久久婷婷六月综合欲色啪| 欧美乱码精品一区二区三区| 一个人看的www免费观看视频| 大型黄色视频在线免费观看| 国产欧美日韩一区二区三| or卡值多少钱| 桃红色精品国产亚洲av| 99久久综合精品五月天人人| 最新美女视频免费是黄的| 人妻丰满熟妇av一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国产单亲对白刺激| 日本黄色片子视频| 亚洲午夜精品一区,二区,三区| 亚洲激情在线av| 在线观看免费午夜福利视频| 白带黄色成豆腐渣| 国产亚洲欧美在线一区二区| 亚洲男人的天堂狠狠| 亚洲精品在线美女| 欧美成人性av电影在线观看| 日韩大尺度精品在线看网址| 两个人的视频大全免费| 国产成人精品久久二区二区91| www.精华液| 午夜福利在线观看免费完整高清在 | 色视频www国产| 欧美xxxx黑人xx丫x性爽| 国模一区二区三区四区视频 | 99视频精品全部免费 在线 | 99久久成人亚洲精品观看| 中文字幕av在线有码专区| 国产一区二区三区视频了| 在线视频色国产色| 久久精品综合一区二区三区| 给我免费播放毛片高清在线观看| 级片在线观看| 手机成人av网站| 亚洲国产精品sss在线观看| 色综合站精品国产| 毛片女人毛片| 美女 人体艺术 gogo| 三级毛片av免费| 国产激情久久老熟女| 午夜日韩欧美国产| 日韩欧美国产一区二区入口| 床上黄色一级片| 老司机午夜福利在线观看视频| 两性夫妻黄色片| 女人被狂操c到高潮| 又黄又粗又硬又大视频| 国产av在哪里看| 90打野战视频偷拍视频| 熟女人妻精品中文字幕| 99热这里只有是精品50| 国产视频一区二区在线看| 国产欧美日韩一区二区三| 久久亚洲真实| 亚洲第一电影网av| 亚洲一区二区三区不卡视频| 日本五十路高清| 国产精品av视频在线免费观看| 亚洲,欧美精品.| 国产亚洲精品久久久久久毛片| 亚洲成av人片在线播放无| 精品人妻1区二区| 国产黄a三级三级三级人| 日韩欧美免费精品| www.www免费av| 亚洲人成电影免费在线| 一本精品99久久精品77| 超碰成人久久| 两个人的视频大全免费| 国产精品香港三级国产av潘金莲| 亚洲国产看品久久| 国产黄片美女视频| 午夜福利在线在线| 激情在线观看视频在线高清| 人人妻,人人澡人人爽秒播| 老汉色av国产亚洲站长工具| 女警被强在线播放| 在线国产一区二区在线| 人人妻人人澡欧美一区二区| 亚洲av电影在线进入| 国产高清有码在线观看视频| 在线播放国产精品三级| 在线观看免费视频日本深夜| 国产不卡一卡二| 俄罗斯特黄特色一大片| 在线观看日韩欧美| 91九色精品人成在线观看| 最好的美女福利视频网| 欧美zozozo另类| 欧美大码av| 精品国产超薄肉色丝袜足j| 亚洲av成人av| 一个人免费在线观看电影 | 成人特级黄色片久久久久久久| 十八禁网站免费在线| 一进一出抽搐动态| 久久精品亚洲精品国产色婷小说| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 久久精品国产99精品国产亚洲性色| 少妇熟女aⅴ在线视频| 在线观看一区二区三区| 高清毛片免费观看视频网站| 久久午夜综合久久蜜桃| 国产高潮美女av| 亚洲七黄色美女视频| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 日韩三级视频一区二区三区| 桃色一区二区三区在线观看| 一级黄色大片毛片| 亚洲av成人av| 三级毛片av免费| 国产精品国产高清国产av| 中文字幕人妻丝袜一区二区| e午夜精品久久久久久久| 久久国产精品影院| 亚洲性夜色夜夜综合| www国产在线视频色| 99热精品在线国产| 最好的美女福利视频网| 给我免费播放毛片高清在线观看| 好看av亚洲va欧美ⅴa在| 国产欧美日韩一区二区三| 免费在线观看亚洲国产| 在线a可以看的网站| 国产av不卡久久| xxxwww97欧美| 国语自产精品视频在线第100页| 国内揄拍国产精品人妻在线| 高清毛片免费观看视频网站| 国产又黄又爽又无遮挡在线| 日本黄大片高清| 免费电影在线观看免费观看| 麻豆国产av国片精品| 窝窝影院91人妻| 一区二区三区高清视频在线| 19禁男女啪啪无遮挡网站| 亚洲人成电影免费在线| 欧美日本视频| 最新中文字幕久久久久 | 国产精品永久免费网站| 亚洲中文字幕一区二区三区有码在线看 | 亚洲,欧美精品.| 久久性视频一级片| 欧美一级a爱片免费观看看| 熟女少妇亚洲综合色aaa.| 老司机午夜福利在线观看视频| 亚洲无线观看免费| 日韩中文字幕欧美一区二区| 成人特级av手机在线观看| 免费看日本二区| 小说图片视频综合网站| 久久精品综合一区二区三区| 国产成人福利小说| 国产欧美日韩精品一区二区| 欧美日韩国产亚洲二区| 岛国在线观看网站| 床上黄色一级片| 97超级碰碰碰精品色视频在线观看| 曰老女人黄片| ponron亚洲| 午夜视频精品福利| 亚洲国产高清在线一区二区三| 国产精品久久久av美女十八| 波多野结衣高清作品| 午夜福利在线观看吧| 国产av麻豆久久久久久久| 国产精品久久久av美女十八| 黄色片一级片一级黄色片| 久久热在线av| 午夜福利18| 1024手机看黄色片| 欧美不卡视频在线免费观看| 九色成人免费人妻av| 人妻久久中文字幕网| 一区二区三区激情视频| 岛国在线免费视频观看| av黄色大香蕉| 亚洲成人精品中文字幕电影| 十八禁人妻一区二区| 亚洲av片天天在线观看| 国产蜜桃级精品一区二区三区| 亚洲人与动物交配视频| 国产亚洲av高清不卡| 天天添夜夜摸| 午夜久久久久精精品| av中文乱码字幕在线| 一本久久中文字幕| 国产精品一区二区三区四区免费观看 | www.熟女人妻精品国产| 色在线成人网| 亚洲人与动物交配视频| 中文亚洲av片在线观看爽| 香蕉丝袜av| 亚洲 欧美一区二区三区| 免费看a级黄色片| 中国美女看黄片| 日韩大尺度精品在线看网址| 亚洲美女视频黄频| 性欧美人与动物交配| 国产视频一区二区在线看| 90打野战视频偷拍视频| 一本久久中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 成人三级做爰电影| 国产亚洲精品久久久久久毛片| 一个人看视频在线观看www免费 | xxxwww97欧美| 男人舔女人的私密视频| 国产亚洲av嫩草精品影院| 757午夜福利合集在线观看| 免费大片18禁| 91老司机精品| 最好的美女福利视频网| 这个男人来自地球电影免费观看| 国产探花在线观看一区二区| 在线观看66精品国产| 神马国产精品三级电影在线观看| 啪啪无遮挡十八禁网站| 欧美大码av| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区| 成人三级做爰电影| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 12—13女人毛片做爰片一| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| h日本视频在线播放| 啦啦啦韩国在线观看视频| 757午夜福利合集在线观看| 美女高潮喷水抽搐中文字幕| 欧美成人一区二区免费高清观看 | 91字幕亚洲| 国产黄色小视频在线观看| 亚洲欧美日韩高清在线视频| 日本与韩国留学比较| 观看美女的网站| 中文字幕久久专区| 一区福利在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品合色在线| av在线天堂中文字幕| 亚洲欧美日韩高清在线视频| 少妇丰满av| aaaaa片日本免费| 我的老师免费观看完整版| 国产成人精品无人区| 99国产精品99久久久久| 亚洲美女视频黄频| 夜夜夜夜夜久久久久| 亚洲 欧美 日韩 在线 免费| 啪啪无遮挡十八禁网站| 黑人操中国人逼视频| 久久精品影院6| 色av中文字幕| 亚洲 国产 在线| 性色avwww在线观看| 韩国av一区二区三区四区| 国产精品,欧美在线| 午夜成年电影在线免费观看| 午夜两性在线视频| 99国产极品粉嫩在线观看| 国产一区在线观看成人免费| xxxwww97欧美| 久久欧美精品欧美久久欧美| 99riav亚洲国产免费| 高潮久久久久久久久久久不卡| 色在线成人网| 夜夜夜夜夜久久久久| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 国产视频一区二区在线看| av视频在线观看入口| 亚洲在线观看片| 不卡一级毛片| 国产高清有码在线观看视频| xxxwww97欧美| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频| 精品久久久久久久毛片微露脸| 日韩国内少妇激情av| 国产成人精品久久二区二区91| 黄色成人免费大全| 麻豆成人午夜福利视频| 午夜福利18| 午夜成年电影在线免费观看| www.999成人在线观看| 人人妻人人看人人澡| 免费大片18禁| 麻豆成人av在线观看| 欧美色视频一区免费| 午夜福利在线在线| 美女扒开内裤让男人捅视频| 色综合亚洲欧美另类图片| 国内精品久久久久精免费| 淫秽高清视频在线观看| 欧美绝顶高潮抽搐喷水| 性色avwww在线观看| 一进一出抽搐gif免费好疼| 熟女人妻精品中文字幕| 91av网一区二区| 亚洲欧洲精品一区二区精品久久久| 精品乱码久久久久久99久播| 国产av麻豆久久久久久久| 免费看美女性在线毛片视频| 免费无遮挡裸体视频| 毛片女人毛片| 男女午夜视频在线观看| 高清在线国产一区| 国产一区二区激情短视频| 又黄又粗又硬又大视频| www.www免费av| 亚洲专区中文字幕在线| 久久久久性生活片| 婷婷六月久久综合丁香| 真实男女啪啪啪动态图| 久久香蕉国产精品| 搞女人的毛片| 美女cb高潮喷水在线观看 | 毛片女人毛片| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式 | 精品熟女少妇八av免费久了| 人人妻人人澡欧美一区二区| 国产高清videossex| 久久久水蜜桃国产精品网| 无限看片的www在线观看| 又大又爽又粗| 欧美+亚洲+日韩+国产| 精品久久久久久成人av| 男女床上黄色一级片免费看| 又黄又粗又硬又大视频| 给我免费播放毛片高清在线观看| 亚洲美女黄片视频| 精品久久久久久成人av| 女生性感内裤真人,穿戴方法视频| 韩国av一区二区三区四区| 精品久久久久久久末码| 97碰自拍视频| 国产av一区在线观看免费| 久久天堂一区二区三区四区| 亚洲成av人片在线播放无| 久久天躁狠狠躁夜夜2o2o| 人妻夜夜爽99麻豆av| 成年女人毛片免费观看观看9| 黄色女人牲交| 国产欧美日韩精品一区二区| 色播亚洲综合网| ponron亚洲| 校园春色视频在线观看| www.999成人在线观看| 免费在线观看视频国产中文字幕亚洲| 久久亚洲真实| 97碰自拍视频| av黄色大香蕉| 偷拍熟女少妇极品色| 神马国产精品三级电影在线观看| 色av中文字幕| 国产一区二区在线av高清观看| 久久精品亚洲精品国产色婷小说| 毛片女人毛片| 美女免费视频网站| 欧美最黄视频在线播放免费| 亚洲七黄色美女视频| 我的老师免费观看完整版| 特大巨黑吊av在线直播| av欧美777| 最新美女视频免费是黄的| 亚洲九九香蕉| 中文亚洲av片在线观看爽| 操出白浆在线播放| 女人高潮潮喷娇喘18禁视频| 中出人妻视频一区二区| 国产午夜福利久久久久久| 欧美黄色片欧美黄色片| 怎么达到女性高潮| 成人特级av手机在线观看| 色综合婷婷激情| 欧美成人性av电影在线观看| www日本在线高清视频| 91av网站免费观看| 色综合欧美亚洲国产小说| 啦啦啦观看免费观看视频高清| 国产视频一区二区在线看| 国产欧美日韩一区二区精品| 久久精品人妻少妇| 99久久综合精品五月天人人| 久久久成人免费电影| av黄色大香蕉| 国产免费av片在线观看野外av| 制服人妻中文乱码| 亚洲国产精品sss在线观看| 亚洲国产中文字幕在线视频| 美女黄网站色视频| 免费在线观看日本一区| 国产一区二区激情短视频| 欧美乱码精品一区二区三区| 日本免费a在线| 一个人看的www免费观看视频| 一边摸一边抽搐一进一小说| 久久国产精品影院| 可以在线观看的亚洲视频| 成人三级黄色视频| 国产精品永久免费网站| 亚洲熟妇熟女久久| 看免费av毛片| 免费大片18禁| 国产视频一区二区在线看| 可以在线观看的亚洲视频| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看| 久久精品国产清高在天天线| 亚洲 欧美一区二区三区| 制服丝袜大香蕉在线| 久久国产精品影院| 在线观看日韩欧美| 少妇人妻一区二区三区视频| 国产久久久一区二区三区| 久久天堂一区二区三区四区| 看免费av毛片| 少妇裸体淫交视频免费看高清| 熟女电影av网| 精品一区二区三区视频在线 | 日韩欧美国产一区二区入口| 欧美黄色淫秽网站| av国产免费在线观看| www.熟女人妻精品国产| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 91麻豆精品激情在线观看国产| 一区二区三区国产精品乱码| 91久久精品国产一区二区成人 | 又黄又粗又硬又大视频| 午夜免费激情av| 国产乱人视频| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜| 国产 一区 欧美 日韩| 午夜福利免费观看在线| 白带黄色成豆腐渣| 久99久视频精品免费| 国产一区在线观看成人免费|