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

    BEPS模型中光合模塊機(jī)理性參數(shù)的適用性1)

    2018-06-15 01:44:06潘天石劉哲范文義
    關(guān)鍵詞:帽兒步長(zhǎng)通量

    潘天石 劉哲 范文義

    (東北林業(yè)大學(xué),哈爾濱,150040) (內(nèi)蒙古赤峰市松山區(qū)水利局) (東北林業(yè)大學(xué))

    碳循環(huán)研究最早開(kāi)始于20世紀(jì)70年代的IBP計(jì)劃,并于1975年完成了對(duì)全球凈初級(jí)生產(chǎn)力(NPP)的估算[1],但未解決由個(gè)體或群體結(jié)果推算大尺度生產(chǎn)力時(shí)的尺度轉(zhuǎn)換問(wèn)題?,F(xiàn)階段,一些碳循環(huán)機(jī)理模型在模擬大空間尺度碳收支方面被廣泛地使用,這些模型大多以植物生理過(guò)程模型為基礎(chǔ),對(duì)植物瞬時(shí)碳同化量做定量分析,并采用積分等不同的數(shù)學(xué)方法對(duì)不同時(shí)間步長(zhǎng)的碳積累量進(jìn)行模擬。

    BEPS由FOREST-BGC模型基礎(chǔ)上改善而來(lái),最早由Liu J et al.[2]利用1 km分辨率MODIS影像和插值得到的氣象數(shù)據(jù),在日步長(zhǎng)尺度上對(duì)加拿大陸地凈初級(jí)生產(chǎn)力進(jìn)行了模擬。該模型成功地解決了時(shí)間及空間由小尺度向大尺度轉(zhuǎn)換的問(wèn)題。在日步長(zhǎng)模型的基礎(chǔ)上,Ju W et al.[3]開(kāi)發(fā)了可以模擬每日NPP變化的小時(shí)步長(zhǎng)碳循環(huán)模型,并用該模型模擬了薩斯喀徹溫省1997—2004年際NPP。BEPS模型驅(qū)動(dòng)數(shù)據(jù),包括冠層聚集度指數(shù)、初始碳庫(kù)值、氣象數(shù)據(jù)(太陽(yáng)總輻射、溫度、濕度、降水、風(fēng)速)等。從模型所用的驅(qū)動(dòng)數(shù)據(jù)不難看出,氣象數(shù)據(jù)是模型運(yùn)行過(guò)程中一項(xiàng)重要參數(shù),這些參數(shù)直接或間接決定了植物光合作用的同化速率,而植物光合作用正是直接作用于植物碳吸收的生理過(guò)程。

    然而,許多學(xué)者未注意到在不同外界條件下,CO2同化量的理論模型存在機(jī)理參數(shù)上的差別,對(duì)最終模型模擬結(jié)果產(chǎn)生影響。Farquhar et al.[4]的研究表明,光合作用模型在不同生理?xiàng)l件下有3套機(jī)理性參數(shù),不同機(jī)理參數(shù)對(duì)CO2的同化速率的模擬結(jié)果有較大影響。該模型雖然指出了不同生理?xiàng)l件下的模型機(jī)理參數(shù)會(huì)有不同,但并未具體闡述外界條件變化對(duì)植物生理過(guò)程有何種具體影響。在應(yīng)用以宏觀氣象數(shù)據(jù)驅(qū)動(dòng)的BEPS模型進(jìn)行模擬時(shí),無(wú)法針對(duì)相應(yīng)的外界條件選擇最優(yōu)的參數(shù)。因此,分析不同參數(shù)在不同外界條件下的適用性,對(duì)提高模型模擬精度顯得十分重要。本研究以BEPS模型作為平臺(tái),將太陽(yáng)總輻射、氣溫、降水等宏觀外界條件視為模型模擬精度的主要因素進(jìn)行分析,通過(guò)實(shí)測(cè)數(shù)據(jù)驗(yàn)證,選取適合研究區(qū)不同外界條件的光合作用機(jī)理性參數(shù)。

    1 研究區(qū)概況

    帽兒山實(shí)驗(yàn)林場(chǎng)(127°30′~127° 34′ E,45° 20′~45°25′N(xiāo))位于黑龍江省東南部尚志市帽兒山鎮(zhèn)內(nèi),總面積約為26 620 hm2。東北林業(yè)大學(xué)帽兒山生態(tài)站位于老爺嶺施業(yè)區(qū)內(nèi),平均海拔400 m,平均坡度15°,主要土壤類(lèi)型為暗棕壤;氣候?qū)俅箨懶约撅L(fēng)氣候,夏季溫暖濕潤(rùn),冬季寒冷干燥,年平均氣溫3.1 ℃。該地區(qū)的植物區(qū)系為長(zhǎng)白植物區(qū)系,原始地帶性植被為以紅松為主的針闊混交林,經(jīng)過(guò)不同程度的干擾(采伐、經(jīng)營(yíng)、火燒、開(kāi)墾等)后,形成的東北東部山區(qū)典型天然次生林和人工林。主要樹(shù)種有樟子松(Pinussylvestrisvar.mongolica)、紅松(Pinuskoraiensis)、落葉松(Larixgmmelinii)、白樺(Betulaplatyphylla)、蒙古櫟(Quercusmongolica)、山楊(Populusdavidiana)等。

    2 研究方法

    2.1 數(shù)據(jù)來(lái)源

    本研究所用數(shù)據(jù)主要為逐日氣象數(shù)據(jù)、逐日葉面積指數(shù)(LAI)數(shù)據(jù)和逐日總初級(jí)生產(chǎn)力(GPP)數(shù)據(jù)。逐日氣象數(shù)據(jù)來(lái)自帽兒山生態(tài)實(shí)驗(yàn)站2011年測(cè)定數(shù)據(jù),包括氣溫、相對(duì)濕度、降水量、風(fēng)速、太陽(yáng)總輻射日值等數(shù)據(jù);逐日LAI數(shù)據(jù)通過(guò)大量的樣地實(shí)測(cè)數(shù)據(jù)進(jìn)行3次樣條插值方法提取[5];逐日GPP數(shù)據(jù)為帽兒山生態(tài)站通量塔的日步長(zhǎng)觀測(cè)數(shù)據(jù)。

    2.2 光合作用過(guò)程中電子傳遞速率的確定

    植物光合作用主要由3個(gè)過(guò)程組成,即光反應(yīng)過(guò)程、光合碳還原過(guò)程(PCR)和光合碳氧化過(guò)程(PCO)。Farquhar et al.[4]在相應(yīng)外界條件下,通過(guò)對(duì)這些過(guò)程中各物質(zhì)的產(chǎn)生和消耗進(jìn)行量化,從而對(duì)CO2同化速率做出定量描述,而B(niǎo)EPS光合模塊是通過(guò)對(duì)光合碳還原過(guò)程的描述來(lái)模擬CO2的同化量。

    光合碳還原過(guò)程中需要三磷酸腺苷(ATP)供能,還原型輔酶II(NADPH)作為還原劑,二者的作用與電子傳遞速率息息相關(guān)。不同條件下光合作用CO2同化速率公式及機(jī)理參數(shù)見(jiàn)表1。

    從CO2同化速率的量化過(guò)程可以看出,模擬光合速率的模型機(jī)理性參數(shù)并不唯一,3種對(duì)光合速率的描述具有顯著差異。用Farquhar模型模擬光合作用過(guò)程,必然涉及到參數(shù)的確定問(wèn)題。選用不同參數(shù)的模型對(duì)光和速率的模擬精度會(huì)有一定影響。

    表1不同生理?xiàng)l件下BEPS模型模擬CO2同化時(shí)所采用參數(shù)

    生理?xiàng)l件CO2同化速率參數(shù)受NADPH供給量限制A=JC-Γ?4C+8Γ?-Rd4和8受ATP供給量限制A=JC-Γ?4.5C+10.5Γ?-Rd4.5和10.5當(dāng)DHAP/PGA穿梭發(fā)生時(shí)A=JC-Γ?6C+18Γ?-Rd6和18

    注:表中J為供給量限制下的電子傳遞速率,Rd為光照條件下線(xiàn)粒體產(chǎn)生CO2量[6-9],C為葉內(nèi)CO2分壓[10-16];Γ*為CO2光補(bǔ)償點(diǎn)[17-21]。

    確定三套研究參數(shù)后,在模型程序中對(duì)應(yīng)模塊進(jìn)行參數(shù)修改,得到不同參數(shù)下的模擬結(jié)果,并以通量塔實(shí)測(cè)數(shù)據(jù)為參考進(jìn)行模擬精度對(duì)比驗(yàn)證,確定最優(yōu)參數(shù)組合,并結(jié)合氣象數(shù)據(jù),分析不同外界條件下參數(shù)的適用性。根據(jù)確定系數(shù)(R2)、均方根誤差(RMSE)、絕對(duì)誤差(Bias)評(píng)價(jià)機(jī)理參數(shù)的模擬結(jié)果。

    3 結(jié)果與分析

    3.1 帽兒山生態(tài)站全年GPP模擬

    由圖1可知,在本研究中,利用BEPS模型模擬帽兒山生態(tài)站單點(diǎn)全年每日GPP。運(yùn)行前先對(duì)模型中的生理參數(shù)進(jìn)行修改,以模擬不同條件下的光合速率。當(dāng)植物處于生長(zhǎng)季時(shí),除少數(shù)幾天外,選用參數(shù)4和8的模擬結(jié)果要比選用參數(shù)4.5和10.5的模擬結(jié)果以及選用參數(shù)6和18的模擬結(jié)果好,說(shuō)明選用參數(shù)4和8模擬的GPP是理想情況下的最大值。

    相較于其它兩套參數(shù)模擬的GPP,選用參數(shù)4.5和10.5的模擬結(jié)果更接近通量塔所測(cè)GPP,該套參數(shù)模擬的GPP描述的是光合速率的平均水平。而選用參數(shù)6和18的模擬結(jié)果在生長(zhǎng)季偏低,非生長(zhǎng)季的模擬結(jié)果又偏高,在實(shí)際模擬過(guò)程中并不能得到很好的結(jié)果,所以在后面的討論中將不再對(duì)這套參數(shù)進(jìn)行分析。

    3.2 選用不同生理參數(shù)模擬結(jié)果對(duì)比

    從圖2中可以看出,選用兩組不同參數(shù)模擬的年際GPP與通量塔實(shí)測(cè)年際GPP呈顯著的線(xiàn)性關(guān)系。

    圖1 選用同生理參數(shù)的BEPS日步長(zhǎng)模型模擬帽兒山生態(tài)站全年GPP與通量塔實(shí)測(cè)GPP的對(duì)比

    圖2選用不同生理參數(shù)下實(shí)測(cè)日GPP與模擬日GPP的線(xiàn)性回歸分析

    表2 選用不同參數(shù)模擬生長(zhǎng)季GPP精度對(duì)比

    由表2可知,在模擬生長(zhǎng)季GPP時(shí),選用參數(shù)4.5和10.5模擬結(jié)果優(yōu)于選用參數(shù)4和8的模擬結(jié)果。

    定義兩套參數(shù)模擬GPP的差異率為:(選用參數(shù)4和8模擬結(jié)果-選用參數(shù)4.5和10.5的模擬結(jié)果)/選用參數(shù)4.5和10.5的模擬結(jié)果。

    通過(guò)圖3可以看出,選用兩套不同生理參數(shù)在模擬生長(zhǎng)季GPP時(shí),模擬結(jié)果存在明顯差異,一般相差6%左右,最多相差14.08%。造成這種結(jié)果的原因是梅勒反應(yīng)中的電子傳輸速率大概為潛在電子傳遞速率的16%[4]。

    3.3 不同生理參數(shù)組合模擬結(jié)果差異對(duì)外界條件的響應(yīng)

    通過(guò)圖4可以看出隨著光合有效輻射(PAR)的增加,差值有明顯的增大趨勢(shì),二者間存在線(xiàn)性相關(guān)關(guān)系(R2=0.475 8),說(shuō)明隨PAR的增強(qiáng),選用兩套參數(shù)模擬結(jié)果差異性增強(qiáng);溫度對(duì)差值也有一定的影響,二者間存在線(xiàn)性關(guān)系(R2=0.619 7),差值隨溫度的增大而增大,隨溫度的升高,兩套參數(shù)模擬結(jié)果差異性增強(qiáng);而相對(duì)濕度和降水對(duì)差值的影響微乎其微,R2均小于0.1,相關(guān)關(guān)系極弱,其變化對(duì)模擬結(jié)果的差異沒(méi)有貢獻(xiàn)。結(jié)合表2的結(jié)果可知,選用參數(shù)4.5和10.5的結(jié)果要更接近于實(shí)測(cè)值,而隨著PAR的增強(qiáng)和溫度的升高,選用參數(shù)4和8的模擬結(jié)果與選用參數(shù)4.5和10.5的模擬結(jié)果差異加大,說(shuō)明選用參數(shù)4和8的模擬結(jié)果與真實(shí)值的偏差增大,選用參數(shù)4.5和10.5的模擬效果優(yōu)于選用參數(shù)4和8的模擬效果。

    圖3 選用不同生理參數(shù)模擬GPP時(shí)的差異率

    4 結(jié)論

    (1)通過(guò)對(duì)比選用不同生理參數(shù)的日步長(zhǎng)BEPS模型模擬GPP值與通量塔實(shí)測(cè)GPP值可以發(fā)現(xiàn),選用參數(shù)4和8的模擬結(jié)果絕大多數(shù)都大于真實(shí)值,可視為模型模擬每日GPP的上限;選用參數(shù)4.5和10.5的模擬結(jié)果雖然也大于真實(shí)值,但與真實(shí)值十分接近,代表了模型模擬每日GPP的平均水平,是利用模型模擬每日GPP的最佳參數(shù);選用參數(shù)6和18的模擬結(jié)果低估了生長(zhǎng)季GPP,又高估了非生長(zhǎng)季GPP,不適于作為機(jī)理參數(shù)模擬每日GPP。

    圖4 選用不同生理參數(shù)模擬GPP時(shí)的精度差異對(duì)外界條件的響應(yīng)

    (2)通過(guò)對(duì)比不同生理參數(shù)模擬GPP的差異,發(fā)現(xiàn)選用參數(shù)4和8的模擬結(jié)果與選用參數(shù)4.5和10.5的模擬結(jié)果最大相差14.08%,平均相差6%,選用不同參數(shù)時(shí)的模擬結(jié)果存在較大差異。

    (3)通過(guò)分析現(xiàn)有日步長(zhǎng)數(shù)據(jù)也得到了趨勢(shì)性的結(jié)論,即在模擬每日GPP時(shí),隨著PAR的增強(qiáng)和溫度的升高,選用參數(shù)4和8的模擬結(jié)果與真實(shí)值之間的差異會(huì)顯著增大,此時(shí)選用參數(shù)4.5和10.5的模擬結(jié)果會(huì)更接近于真實(shí)值。

    雖然滿(mǎn)足ATP消耗下電子傳遞速率計(jì)算的CO2同化速率更具有普適性,但并不是說(shuō)滿(mǎn)足NADPH消耗下電子傳遞速率計(jì)算的CO2同化速率是錯(cuò)誤的。本研究中所用數(shù)據(jù)為通量塔所測(cè)的日步長(zhǎng)數(shù)據(jù),缺少小時(shí)步長(zhǎng)的通量數(shù)據(jù),沒(méi)能通過(guò)對(duì)比半小時(shí)數(shù)據(jù)、日步長(zhǎng)數(shù)據(jù)和真實(shí)數(shù)據(jù)間的關(guān)系。因此,為了模擬結(jié)果更加科學(xué)合理,需要進(jìn)一步獲取準(zhǔn)確的小時(shí)步長(zhǎng)通量數(shù)據(jù)。

    參 考 文 獻(xiàn)

    [1] HELMUT L, WHITTAKER R H. Primary productivity of the biosphere[M]. New York: Springer,1975:305-328.

    [2] LIU J, CHEN J M, CIHLAR J, et al. Net primary productivity mapped for Canada at 1-km resolution[J]. Global Ecology & Biogeography,2002,11(2):115-129.

    [3] JU W, CHEN J M, BLACK T A, et al. Modelling multi-year coupled carbon and water fluxes in a boreal aspen forest[J]. Agricultural & Forest Meteorology,2006,140(1):136-151.

    [4] FARQUHAR G D, CAEMMERER S V. Modeling of photosynthetic response to environmental conditions[J]. Encyclopedia of Plant Physiology New,1982(5):550-584.

    [5] 盧偉,范文義.中小區(qū)域尺度時(shí)間序列林地LAI快速估測(cè)方法[J].農(nóng)業(yè)工程學(xué)報(bào),2016(5):188-193.

    [6] GRAHAM D, CHAPMAN E A. Interactions Between Photosynthesis and Respiration in Higher Plants[M]//Martin Gibbs, Erwin Latzko. Photosynthesis II. Berlin Heidelberg: Springer,1979:150-162.

    [7] MANGAT B S, LEVIN W B, BIDWELL R G S. Extent of dark respiration in illuminated leaves and its control by ATP levels[J]. Canadian Journal of Botany,2011,52(4):673-68.

    [8] PEISKER M, TICHA I, CATSKY J. Ontogenetic changes in the internal limitations to bean-leaf photosynthesis.7. Interpretation of the linear correlation between CO2compensation concentration and CO2evolution in darkness[J]. Photosynthetica,1981,15(2):161-168.

    [9] SOMERVILLE C R, OGREN W L. A phosphoglycolate phosphatase-deficient mutant of Arabidopsis[J]. Nature,1979,280:833-836.

    [10] DANIEL E. The temperature dependence of photoinhibition in leaves ofPhaseolusvulgaris(L.) Influence of CO2and O2concentrations[J]. Plant Science,1997,124(1):1-8.

    [11] VON C S, FARQUHAR G D. Some relationships between the biochemistry of photosynthesis and the ags exchange of leaves[J]. Planta,1981,153(4):376-387.

    [12] TENHUNEN J D, WEBER J A, FILIPEK L H, et al. Development of a photosynthesis model with an emphasis on ecological applications: III. Carbon dioxide and oxygen dependencies[J]. Oecologia,1979,41(2):145-162.

    [13] FARQUHAR G D, CAEMMERER S V, BERRY J A. A biochemical model of photosynthetic CO2assimilation in leaves of C3 species[J]. Planta,1980,149(1):78-90.

    [14] BAUWE H, APEL P, PEISKER M. Ribulose 1,5-bisphosphate carboxylase/oxygenase and CO2exchange characteristics in C3 and C3-C4 intermediate species checking mathematical models of carbon metabolism[J]. Photosynthetica,1980,14(3):550-556.

    [15] ALLEN J F. Cyclic, pseudocyclic and noncyclic photophosphorylation: new links in the chain[J]. Trends in Plant Science,2003,8(1):15-19.

    [16] SHAVIT N. Energy transduction in chloroplasts: structure and function of the ATPase complex[J]. Annual Review of Biochemistry,1980,49(1):111.

    [17] NOLAN W G, SMILLIE R M. Multi-temperature effects on Hill reaction activity of barley chloroplasts[J]. Biochim Biophys Acta,1976,440(3):461-475.

    [18] WEIS E. Short term acclimation of spinach to high temperatures: effect on chlorophyll fluorescence at 293 and 77 kelvin in intact leaves[J]. Plant Physiology,1984,74(2):402-407.

    [19] JOHN F A. Cyclic, pseudocyclic and noncyclic photophosphorylation: new links in the chain[J]. Trends in Plant Science,2003,8(1):15-19.

    [20] FINAZZI G, RAPPAPORT F, FURIA A, et al. Involvement of state transitions in the switch between linear and cyclic electron flow in chlamydomonas reinhardtii[J]. Embo Reports,2002,3(3):280-285.

    [21] ARMOND P A , SCHREIBER U, BJ?RKMAN O. Photosynthetic acclimation to temperature in the desert shrub, larrea divaricata: II. light-harvesting efficiency and electron transport 1[J]. Plant Physiology,1978,61(3):411-415.

    猜你喜歡
    帽兒步長(zhǎng)通量
    冬小麥田N2O通量研究
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    看山
    帽兒山景區(qū)民宿旅游發(fā)展策略研究
    海邊有棟小別墅
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    近山侏羅紀(jì)
    婦女之友(2014年9期)2014-09-26 09:41:25
    一種新型光伏系統(tǒng)MPPT變步長(zhǎng)滯環(huán)比較P&O法
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    又黄又爽又免费观看的视频| 99精品在免费线老司机午夜| 午夜免费激情av| 午夜影院日韩av| 熟女电影av网| 国产亚洲精品综合一区在线观看| 91在线精品国自产拍蜜月| 欧美激情久久久久久爽电影| 校园春色视频在线观看| 成人特级av手机在线观看| 在线观看av片永久免费下载| 国产av一区在线观看免费| 97热精品久久久久久| 18禁黄网站禁片免费观看直播| 午夜精品久久久久久毛片777| 欧美+亚洲+日韩+国产| 国产69精品久久久久777片| 国产私拍福利视频在线观看| 欧美高清成人免费视频www| 男女之事视频高清在线观看| 一级黄色大片毛片| 久久亚洲真实| 亚洲电影在线观看av| 免费搜索国产男女视频| 人妻夜夜爽99麻豆av| 中文字幕免费在线视频6| 亚洲自偷自拍三级| 少妇的逼水好多| 亚洲国产欧美人成| 一进一出抽搐gif免费好疼| 99热这里只有精品一区| 五月玫瑰六月丁香| 日韩,欧美,国产一区二区三区 | 国产高潮美女av| 国产成人a区在线观看| 国产精品嫩草影院av在线观看 | 最近最新免费中文字幕在线| 精品久久久久久久久亚洲 | 亚洲最大成人av| 久久精品夜夜夜夜夜久久蜜豆| 在线观看免费视频日本深夜| 久久久久久久亚洲中文字幕| 一进一出抽搐动态| 亚洲av一区综合| 天堂动漫精品| 自拍偷自拍亚洲精品老妇| 中文资源天堂在线| 国产在视频线在精品| 美女 人体艺术 gogo| 亚洲精品影视一区二区三区av| 亚洲精品日韩av片在线观看| 国产精品不卡视频一区二区| 天堂√8在线中文| 亚洲,欧美,日韩| 一边摸一边抽搐一进一小说| 男人和女人高潮做爰伦理| 日韩欧美在线乱码| 日韩亚洲欧美综合| 国内精品宾馆在线| 日本a在线网址| 美女免费视频网站| 精品不卡国产一区二区三区| 国产欧美日韩精品一区二区| 亚洲精品乱码久久久v下载方式| 日本欧美国产在线视频| 变态另类丝袜制服| 日韩欧美国产一区二区入口| 99久久精品一区二区三区| 国产中年淑女户外野战色| 日日干狠狠操夜夜爽| 香蕉av资源在线| 日韩欧美精品v在线| 国产高清不卡午夜福利| 日日夜夜操网爽| 亚洲男人的天堂狠狠| 乱系列少妇在线播放| 99久久无色码亚洲精品果冻| 中文字幕熟女人妻在线| 久久精品久久久久久噜噜老黄 | 国内精品一区二区在线观看| 国产精品电影一区二区三区| 久久久久久久久久久丰满 | 久久99热6这里只有精品| 免费搜索国产男女视频| 日本精品一区二区三区蜜桃| 成人三级黄色视频| 国产精品一区二区性色av| 国产黄片美女视频| 免费一级毛片在线播放高清视频| 久久久国产成人精品二区| 97热精品久久久久久| 有码 亚洲区| 91av网一区二区| 亚洲av二区三区四区| 中文字幕av在线有码专区| 夜夜夜夜夜久久久久| 久久天躁狠狠躁夜夜2o2o| 琪琪午夜伦伦电影理论片6080| 久久久午夜欧美精品| 国产毛片a区久久久久| 人妻丰满熟妇av一区二区三区| 久久草成人影院| 麻豆国产av国片精品| 91久久精品国产一区二区成人| 高清日韩中文字幕在线| 老熟妇仑乱视频hdxx| 小蜜桃在线观看免费完整版高清| 中出人妻视频一区二区| 最近中文字幕高清免费大全6 | 日本精品一区二区三区蜜桃| 成人av一区二区三区在线看| 91狼人影院| 亚洲在线观看片| 干丝袜人妻中文字幕| 嫁个100分男人电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 一区福利在线观看| 看十八女毛片水多多多| 久久人妻av系列| 啦啦啦啦在线视频资源| videossex国产| 亚洲最大成人av| 亚洲国产欧洲综合997久久,| 国产女主播在线喷水免费视频网站 | 一区二区三区免费毛片| 欧美精品国产亚洲| 少妇丰满av| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美日韩无卡精品| 悠悠久久av| 国产精品98久久久久久宅男小说| 亚洲,欧美,日韩| 又爽又黄a免费视频| 免费看日本二区| 久99久视频精品免费| 精华霜和精华液先用哪个| 亚洲一区高清亚洲精品| 亚洲欧美激情综合另类| 少妇猛男粗大的猛烈进出视频 | 中文亚洲av片在线观看爽| 国语自产精品视频在线第100页| 成人高潮视频无遮挡免费网站| 亚洲va日本ⅴa欧美va伊人久久| 国内精品一区二区在线观看| 色5月婷婷丁香| 日日夜夜操网爽| 国产伦精品一区二区三区四那| 国产欧美日韩精品一区二区| 国产aⅴ精品一区二区三区波| 黄色丝袜av网址大全| 国产aⅴ精品一区二区三区波| 亚洲avbb在线观看| 午夜福利成人在线免费观看| 久久精品国产亚洲av天美| 国产久久久一区二区三区| 亚洲avbb在线观看| 国产久久久一区二区三区| 久久草成人影院| 亚洲精品乱码久久久v下载方式| 国产精品永久免费网站| 亚洲不卡免费看| 成人亚洲精品av一区二区| 久久精品国产亚洲网站| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久爱视频| 精品一区二区三区视频在线| 精品久久久久久久末码| 色哟哟·www| 久久精品人妻少妇| 国产一级毛片七仙女欲春2| 欧美精品啪啪一区二区三区| 91在线观看av| 国产精品电影一区二区三区| 啦啦啦观看免费观看视频高清| 联通29元200g的流量卡| 国产精品久久电影中文字幕| 久久精品国产亚洲网站| 午夜精品久久久久久毛片777| 亚洲精华国产精华精| 日本 欧美在线| 嫩草影院入口| 最近最新免费中文字幕在线| a在线观看视频网站| 嫩草影院入口| 精品久久久久久久末码| videossex国产| 一区二区三区免费毛片| 欧美黑人巨大hd| 亚洲五月天丁香| 一边摸一边抽搐一进一小说| 日韩高清综合在线| 亚洲黑人精品在线| 一边摸一边抽搐一进一小说| 禁无遮挡网站| 美女 人体艺术 gogo| 国产v大片淫在线免费观看| 91麻豆av在线| 国产色爽女视频免费观看| 欧美潮喷喷水| 欧美日本视频| 亚洲经典国产精华液单| 夜夜看夜夜爽夜夜摸| 男人舔奶头视频| 亚洲成人中文字幕在线播放| 欧美bdsm另类| 色吧在线观看| 丰满乱子伦码专区| 亚洲欧美日韩卡通动漫| 免费观看的影片在线观看| 午夜激情欧美在线| 国产美女午夜福利| 色视频www国产| 女的被弄到高潮叫床怎么办 | 欧美一区二区精品小视频在线| 国产一区二区三区视频了| 亚洲精品一区av在线观看| 一区福利在线观看| 波多野结衣高清作品| 亚洲av第一区精品v没综合| 最近中文字幕高清免费大全6 | 桃红色精品国产亚洲av| 99精品久久久久人妻精品| 热99re8久久精品国产| 桃色一区二区三区在线观看| 国内精品久久久久久久电影| 久久99热6这里只有精品| 校园春色视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 婷婷丁香在线五月| 欧美一级a爱片免费观看看| 免费人成在线观看视频色| aaaaa片日本免费| 成人性生交大片免费视频hd| 直男gayav资源| 中国美白少妇内射xxxbb| 熟女电影av网| 国产午夜福利久久久久久| 一区二区三区激情视频| 白带黄色成豆腐渣| 日韩,欧美,国产一区二区三区 | 亚洲avbb在线观看| 亚洲av成人精品一区久久| 男女那种视频在线观看| 免费人成在线观看视频色| 欧美绝顶高潮抽搐喷水| 国产亚洲精品久久久久久毛片| 特级一级黄色大片| 最近在线观看免费完整版| 免费av不卡在线播放| 成人国产麻豆网| 露出奶头的视频| 国产主播在线观看一区二区| 日日撸夜夜添| 最近视频中文字幕2019在线8| 国产精品自产拍在线观看55亚洲| 在线观看午夜福利视频| 久久精品国产亚洲网站| 亚洲avbb在线观看| 99精品在免费线老司机午夜| 亚洲无线观看免费| 国产视频一区二区在线看| 国产高潮美女av| 亚洲精品成人久久久久久| 88av欧美| 久久久久久久精品吃奶| av女优亚洲男人天堂| 欧美丝袜亚洲另类 | ponron亚洲| 精品一区二区免费观看| 亚洲性夜色夜夜综合| 嫩草影院新地址| 一本一本综合久久| 欧美另类亚洲清纯唯美| 日韩一区二区视频免费看| 黄色日韩在线| 国产三级在线视频| 亚洲精品色激情综合| 日韩av在线大香蕉| 国产中年淑女户外野战色| 丰满乱子伦码专区| 亚洲性夜色夜夜综合| 一个人观看的视频www高清免费观看| 禁无遮挡网站| 亚洲男人的天堂狠狠| 亚洲无线观看免费| 一区二区三区高清视频在线| 亚洲美女视频黄频| 深爱激情五月婷婷| 午夜精品久久久久久毛片777| 色尼玛亚洲综合影院| 欧美人与善性xxx| 欧美一区二区精品小视频在线| 久久99热这里只有精品18| 日本免费a在线| av女优亚洲男人天堂| 中文在线观看免费www的网站| 欧美精品啪啪一区二区三区| 欧美+亚洲+日韩+国产| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| 国产真实乱freesex| 99热这里只有是精品50| 三级男女做爰猛烈吃奶摸视频| 国产男人的电影天堂91| 热99re8久久精品国产| 免费在线观看日本一区| 在线观看66精品国产| 一级a爱片免费观看的视频| 国产精品久久视频播放| 久久久久久久久久久丰满 | 真人做人爱边吃奶动态| 99视频精品全部免费 在线| 国产一区二区三区av在线 | 国内久久婷婷六月综合欲色啪| 国产欧美日韩精品一区二区| 夜夜看夜夜爽夜夜摸| 婷婷精品国产亚洲av在线| 亚洲一级一片aⅴ在线观看| 日日干狠狠操夜夜爽| 一进一出抽搐动态| 免费黄网站久久成人精品| 日日夜夜操网爽| 麻豆成人av在线观看| 亚洲精品在线观看二区| 深夜精品福利| 两个人视频免费观看高清| 精品一区二区三区视频在线| 日韩欧美国产在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产三级中文精品| 国产欧美日韩精品亚洲av| 最后的刺客免费高清国语| 成人二区视频| 赤兔流量卡办理| 久久久久久九九精品二区国产| 国产乱人伦免费视频| 国产免费av片在线观看野外av| 国产成人影院久久av| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕高清在线视频| 国产免费一级a男人的天堂| 高清日韩中文字幕在线| 免费搜索国产男女视频| 国产日本99.免费观看| 看片在线看免费视频| 在线天堂最新版资源| 少妇被粗大猛烈的视频| 99九九线精品视频在线观看视频| 亚洲成a人片在线一区二区| 亚洲中文日韩欧美视频| 一区二区三区高清视频在线| 长腿黑丝高跟| 午夜免费成人在线视频| 成年女人看的毛片在线观看| 午夜福利欧美成人| 舔av片在线| 国产精品嫩草影院av在线观看 | 窝窝影院91人妻| 91久久精品电影网| 亚洲精品粉嫩美女一区| 真人一进一出gif抽搐免费| 老女人水多毛片| 欧美最新免费一区二区三区| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 色视频www国产| 久久久久免费精品人妻一区二区| 亚洲第一区二区三区不卡| 香蕉av资源在线| 男插女下体视频免费在线播放| 久久久精品大字幕| 91久久精品电影网| 久久久久久国产a免费观看| 中文字幕精品亚洲无线码一区| 日韩在线高清观看一区二区三区 | 97人妻精品一区二区三区麻豆| 嫩草影院新地址| 如何舔出高潮| 日本成人三级电影网站| 在线免费观看的www视频| 给我免费播放毛片高清在线观看| 成人美女网站在线观看视频| 中文亚洲av片在线观看爽| 国产精品久久久久久精品电影| videossex国产| 日本一本二区三区精品| 日本-黄色视频高清免费观看| 精品午夜福利在线看| 一进一出抽搐动态| 午夜日韩欧美国产| 国产 一区精品| 亚洲精品影视一区二区三区av| av黄色大香蕉| 亚洲一区二区三区色噜噜| 成年女人永久免费观看视频| 在线看三级毛片| 黄色配什么色好看| 免费av毛片视频| 日本欧美国产在线视频| 国产精品电影一区二区三区| 又粗又爽又猛毛片免费看| 超碰av人人做人人爽久久| 人妻久久中文字幕网| 精品不卡国产一区二区三区| 国产男人的电影天堂91| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 久久精品久久久久久噜噜老黄 | 中文字幕熟女人妻在线| 桃色一区二区三区在线观看| h日本视频在线播放| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看 | 亚洲熟妇中文字幕五十中出| 午夜精品久久久久久毛片777| 韩国av一区二区三区四区| 色噜噜av男人的天堂激情| 毛片一级片免费看久久久久 | 在线a可以看的网站| 国产男人的电影天堂91| aaaaa片日本免费| 免费无遮挡裸体视频| 国产精品一区二区三区四区久久| 在线免费观看不下载黄p国产 | a级毛片a级免费在线| 国产精品久久久久久久久免| 白带黄色成豆腐渣| 色吧在线观看| 美女 人体艺术 gogo| 国产精品一区二区三区四区免费观看 | 国产伦精品一区二区三区四那| 亚洲国产精品sss在线观看| 最近最新免费中文字幕在线| 啦啦啦韩国在线观看视频| 日韩欧美精品v在线| 性插视频无遮挡在线免费观看| 免费在线观看日本一区| 无遮挡黄片免费观看| 嫩草影院入口| 成人精品一区二区免费| 久久精品久久久久久噜噜老黄 | 一边摸一边抽搐一进一小说| 日韩一区二区视频免费看| 日韩欧美精品v在线| 久久久色成人| 亚洲,欧美,日韩| 国产91精品成人一区二区三区| 可以在线观看的亚洲视频| 欧美日韩精品成人综合77777| 在线免费观看的www视频| 精品无人区乱码1区二区| 欧美激情久久久久久爽电影| 亚洲精品色激情综合| 麻豆av噜噜一区二区三区| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 成人鲁丝片一二三区免费| 变态另类成人亚洲欧美熟女| 国产高清视频在线观看网站| 日本撒尿小便嘘嘘汇集6| 又紧又爽又黄一区二区| 九九热线精品视视频播放| 欧美另类亚洲清纯唯美| 国产午夜福利久久久久久| 俄罗斯特黄特色一大片| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 男插女下体视频免费在线播放| 18+在线观看网站| 丰满的人妻完整版| 久久精品国产亚洲网站| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 国产综合懂色| 日韩欧美精品v在线| 国产亚洲精品久久久com| 欧美+日韩+精品| 国产男人的电影天堂91| 女的被弄到高潮叫床怎么办 | 97热精品久久久久久| 91在线观看av| 精品乱码久久久久久99久播| 性欧美人与动物交配| 日本-黄色视频高清免费观看| 校园人妻丝袜中文字幕| 99热精品在线国产| 久久精品夜夜夜夜夜久久蜜豆| bbb黄色大片| av福利片在线观看| 噜噜噜噜噜久久久久久91| 亚洲第一电影网av| 99国产极品粉嫩在线观看| 一本久久中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 十八禁网站免费在线| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费成人在线视频| 亚洲最大成人手机在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚州av有码| 一卡2卡三卡四卡精品乱码亚洲| 精品人妻偷拍中文字幕| 亚洲精品久久国产高清桃花| 在线播放无遮挡| 国产私拍福利视频在线观看| 午夜激情福利司机影院| 国产精品99久久久久久久久| 精品国产三级普通话版| 俺也久久电影网| 免费人成在线观看视频色| 亚洲精品在线观看二区| 嫩草影院新地址| 免费黄网站久久成人精品| 国产乱人伦免费视频| 国产成年人精品一区二区| 欧美丝袜亚洲另类 | 极品教师在线免费播放| 国产探花在线观看一区二区| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 日本五十路高清| 久久久国产成人免费| 一级黄色大片毛片| 日本爱情动作片www.在线观看 | 午夜久久久久精精品| av国产免费在线观看| bbb黄色大片| 动漫黄色视频在线观看| 嫩草影院入口| 午夜日韩欧美国产| 国产av不卡久久| 免费黄网站久久成人精品| 亚洲成人精品中文字幕电影| 精品一区二区三区视频在线| 亚洲成人免费电影在线观看| 欧美在线一区亚洲| 国产探花极品一区二区| 亚洲最大成人av| 日韩大尺度精品在线看网址| 久久久成人免费电影| 中出人妻视频一区二区| 午夜福利18| 精品无人区乱码1区二区| 国产一区二区三区视频了| www.www免费av| 老熟妇乱子伦视频在线观看| 国产精品女同一区二区软件 | 国产91精品成人一区二区三区| 蜜桃久久精品国产亚洲av| 人人妻人人澡欧美一区二区| 熟女人妻精品中文字幕| 国产极品精品免费视频能看的| 一边摸一边抽搐一进一小说| 看十八女毛片水多多多| 欧美在线一区亚洲| 国产精品国产高清国产av| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 一级毛片久久久久久久久女| 韩国av一区二区三区四区| a级毛片a级免费在线| 久久久久久久久中文| 有码 亚洲区| 免费在线观看日本一区| 国产av不卡久久| 大又大粗又爽又黄少妇毛片口| 天堂动漫精品| 成人美女网站在线观看视频| 国产精品免费一区二区三区在线| 天堂影院成人在线观看| 亚洲美女视频黄频| 日本在线视频免费播放| 精品国产三级普通话版| 亚洲精品乱码久久久v下载方式| 国产高清激情床上av| 精品久久久久久成人av| 国内揄拍国产精品人妻在线| 国产精品久久电影中文字幕| 国产不卡一卡二| 亚洲男人的天堂狠狠| 成人国产麻豆网| 欧美高清性xxxxhd video| 国产伦精品一区二区三区视频9| 亚洲国产精品久久男人天堂| 久久婷婷人人爽人人干人人爱| 亚洲图色成人| 日韩欧美一区二区三区在线观看| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 精品久久久久久久久久久久久| 日日撸夜夜添| 久久人妻av系列| 亚洲18禁久久av| 欧美潮喷喷水| 久久久久久大精品| a级一级毛片免费在线观看| 琪琪午夜伦伦电影理论片6080| 赤兔流量卡办理| 黄色视频,在线免费观看| 国产精品99久久久久久久久| 国产av一区在线观看免费| 国产蜜桃级精品一区二区三区| 99久久成人亚洲精品观看| 亚洲av美国av| 国产精品电影一区二区三区| 国国产精品蜜臀av免费| 一个人看视频在线观看www免费| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 亚洲国产色片| 久久久午夜欧美精品| 免费看av在线观看网站| 3wmmmm亚洲av在线观看| 欧美另类亚洲清纯唯美| 亚洲av.av天堂|