999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

毛竹林土壤呼吸研究進展

2012-01-26 04:26:19雷海清傅懋毅李正才
浙江林業科技 2012年2期
關鍵詞:模型研究

雷海清,傅懋毅,李正才,童 龍

(1. 中國林業科學院亞熱帶林業研究所,浙江 富陽 311400;2. 浙江省亞熱帶作物研究所,浙江 溫州 325005)

毛竹林土壤呼吸研究進展

雷海清1,2,傅懋毅1*,李正才1,童 龍1

(1. 中國林業科學院亞熱帶林業研究所,浙江 富陽 311400;2. 浙江省亞熱帶作物研究所,浙江 溫州 325005)

介紹和分析了毛竹(Phyllostachys heterocycla cv. pubescens)林土壤呼吸的研究方法、影響因子、時空變化、干擾響應、組分區分、模型構建等方面取得的研究成果,提出了今后毛竹林土壤呼吸研究應該加強:①土壤溫度對毛竹林土壤呼吸影響機制研究;②對毛竹林土壤呼吸的時空變化進行長期的定位研究;③繼續優化數學模型;④采用更多樣化的方法來研究土壤呼吸各組分。

土壤呼吸;毛竹;研究;進展

土壤呼吸是土壤碳庫的主要輸出途徑和大氣CO2的重要來源,土壤以呼吸的形式釋放到大氣中的CO2量是化石燃料釋放量的10倍以上[1],其微小的變化都可能導致大氣CO2濃度的變化[2]。從19世紀20年代開始,對土壤呼吸的研究就已經開始了,早期的研究主要是為了描述土壤的代謝特征。20世紀90年代以來,對土壤呼吸的研究主要受全球氣候變化的驅動。對大氣CO2數據的分析表明,雖然人類活動釋放的碳大部分可能被陸地生物圈吸收[3],但大氣CO2濃度的不斷升高加劇了溫室效應,導致全球變暖,全球變暖又會大大刺激呼吸作用,導致更多的CO2釋放到大氣中。這些研究激發了研究者對土壤呼吸對于全球氣候變化影響的極大關注,使得越來越多的研究者開始關注土壤呼吸的研究。森林生態系統作為陸地生態系統的重要組成部份,儲存了陸地生態系統有機碳地上部分的80%,地下部分的40%[4],研究林地土壤碳循環對于探討全球CO2變化及其影響具有重要意義。

全球竹林總面積約2 200萬hm2,毛竹(Phyllostachys heterocycla cv. pubescens)是其中面積最大、分布最廣的竹種,具有豐富的經濟和生態價值。毛竹林在調控地球系統的大氣CO2濃度和氣候動態方面起著十分關鍵的作用。對毛竹林土壤呼吸的研究吸引了來自土壤學、生態學、森林培育(竹林培育)等多學科專家的共同興趣。為此筆者期望通過對毛竹林土壤呼吸的研究方法、影響因子、時空變化、干擾響應、組分區分、模型構建等方面研究成果的分析以及對未來發展方向的展望,為毛竹林土壤呼吸研究提供參考。

1 土壤呼吸研究概況

1.1 土壤呼吸的定義

土壤呼吸是通過根呼吸、微生物對凋落物和土壤有機質的分解、動物呼吸,從土壤中釋放CO2的生態系統過程。它是陸地生態系統和大氣生態系統之間碳轉移的主要途徑之一,也是全球碳循環中最大的通量之一,在生態系統的碳平衡中扮演重要角色。土壤呼吸通量的微小變化將會對大氣CO2濃度變化產生較大的影響[5]。土壤呼吸通常是通過測定土壤表面的CO2通量得到,在沒有較大淋溶和沉積的情況下,土壤表面的CO2通量等于土壤呼吸速率[6]。

1.2 土壤呼吸的意義

大氣CO2濃度上升引起的溫室效應及其所帶來的一系列生態環境變化已成為生態學、環境科學等領域的熱點問題。碳循環是生態系統物質循環最重要的過程之一,在全球環境變化的背景下,尤其是隨著《京都議定書》的簽署和生效實施后,其在各個尺度上的循環過程及其特征研究更是倍受關注。盡管土壤是一個巨大的碳庫(相當于大氣中CO2貯量的2倍,陸地植被碳貯量的3倍)[4,7],但是土壤每年以呼吸的形式釋放到大氣中的CO2量是化石燃料釋放量的10倍以上[1],其微小的變化都可能導致大氣CO2濃度的變化[2]。

1.3 土壤呼吸的研究歷史

土壤呼吸的研究已有一百余年的歷史,可以追溯到Wollny[8]、Boussingualt & Lery[9]和Moller[10]的文章。20世紀土壤呼吸的研究可以大致分為4個主要時期:

在20世紀初的幾十年,土壤呼吸的研究主要是利用農業土壤在實驗室中進行的,被用于評價土壤肥力和土壤中生物活動的一個指標[11]。在那個時期,產生了一些較為原始的測量土壤呼吸的方法,如堿吸收法[12]、原位測量法[13]等。含水量被認為是土壤呼吸的影響因子之一[14]。對于土壤呼吸的原理,Turpin研究認為土壤釋放的CO2主要源自細菌的分解作用[15]。

從20世紀30年代末到50年代初這段時期,土壤呼吸的研究相對來說不那么活躍,這可能是由于這個時期世界范圍內的社會動蕩造成的。

從50年代末到80年代,土壤呼吸的研究活動又重新開始進行。這一時期在技術上取得了一個重要的進步就是紅外氣體分析儀(IRGA)被用于土壤呼吸的測量,同時有關土壤呼吸調控因子的結果也層出不窮,許多研究證明了土壤呼吸與溫度呈指數相關[16~17],土壤呼吸的組分被清楚地分為自養和異樣呼吸[18~19]兩大類。

20世紀90年代以來,對土壤呼吸的研究主要受全球變化的驅動,測量技術的進步——便攜式紅外氣體分析儀的廣泛應用,使得土壤呼吸研究更加活躍[20]。研究明確了土壤呼吸受溫度、濕度、底物供應等多因子的交互作用[6,21~22],土壤呼吸的時空變化特征[23~24],土壤呼吸對干擾的響應[22~25],土壤呼吸各組分的區分[26~27]以及模型集成與分析[28~29]。

2 毛竹林土壤呼吸的研究進展

2.1 毛竹林土壤呼吸的研究意義

毛竹是我國南方重要的森林資源,隸屬禾本科竹亞科剛竹屬,毛竹林是地球陸地上重要的森林植被類型。由于毛竹特殊的無性繁殖方式,每年有大量新筍長大成竹向周圍擴展。在全球森林面積急劇下降的今天,竹林面積卻以每年3%左右的速度在遞增,這意味著竹林將是一個不斷增大的碳匯[30]。毛竹林土壤呼吸作為毛竹林生態系統碳循環的一個主要組成部分,與毛竹林生態系統的許多組分都有聯系,與營養過程(如分解作用和礦化作用)也密切相關。而且,毛竹林與其他生態類型的土壤呼吸一起,在調控地球系統的大氣CO2濃度和氣候動態、區域及全球尺度的碳循環上都有相當重要的作用。因此,研究毛竹林生態系統碳儲量變化,無疑會對深入研究中國乃至全球森林生態系統碳平衡問題起重要作用。

中國竹林面積以毛竹為主,毛竹大約占總面積的70%,其它竹類面積只占約30%。通過對包括6次森林資源清查在內的歷史竹林數據統計分析,過去50 a中國的竹林面積隨著全國森林面積的增長而增長,其中以近15 a增長最快,未來的50 a內中國竹林面積還將隨著森林面積的繼續增長而增長,竹林碳儲量仍然會保持增長,擁有穩定的碳匯功能。在用竹林現存生物量計算竹林碳儲量時,已經剔出了常年的竹材采伐量。如果考慮到在采伐竹材利用中,絕大部分還會存留較長時間,竹林的實際固碳能力將大于其碳儲量。目前中國竹林碳儲量占森林總碳儲量的比例已超過11%,成為一個舉足輕重的生態系統碳庫單元。隨著經營水平的不斷提高,中國竹林立竹密度還會增大,相應的竹林碳密度也將增大,這意味著未來竹林碳儲量在整個森林碳儲量中所占的比例還可能增大[31]。

2.2 毛竹林土壤呼吸的研究方法

目前對毛竹林土壤呼吸研究的測定方法主要集中在兩種方法上,一種是采用封閉式動態氣室法[30,32~39],還有一種是采用封閉式靜態氣室法[40]。對于采用其他方法來研究土壤呼吸,目前還未見報道。

2.3 毛竹林土壤呼吸的影響因子

目前研究已經明確的影響土壤呼吸的因子有溫度、濕度、氧氣、氮、土壤質地和土壤pH等。姜艷[30]、范少輝[32]、王燕[33]、李正才[34]、王超[35]、李雅紅[36]、劉勝[37]、劉源月[38]、肖復明[39]等人研究了土壤溫度對毛竹林土壤呼吸的影響,結果表明毛竹林土壤總呼吸及組分呼吸均與空氣溫度、土壤溫度變化相類似。李雅紅[36]、劉源月[38]對土壤濕度對土壤呼吸的影響做了初步研究,李雅紅研究結果表明土壤體積含水量與土壤呼吸不存在顯著關系,土壤體積含水量僅可以解釋土壤呼吸變化的14.19%[36];劉源月研究結果表明,土壤水分與土壤呼吸呈負相關[38]。對于土壤氧氣、土壤氮、土壤質地和土壤pH值對毛竹林土壤呼吸的研究,目前還未見報道。

2.4 毛竹林土壤呼吸的時空變化

土壤呼吸的時間變化包括日變化、季節變化、年際變化。姜艷[30]、范少輝[32]、王燕[33]、李正才[34]、王超[35]、李雅紅[36]、劉勝[37]、劉源月[38]、肖復明[39]、葉耿平[40]都研究了毛竹林土壤呼吸的日變化,結果基本一致,毛竹林土壤呼吸在中午至半下午時達到高峰,然后下午和整個夜間都隨溫度的降低而下降。姜艷[30]、范少輝[32]、王燕[33]、王超[35]、李雅紅[36]、劉勝[37]、劉源月[38]、肖復明[39]、葉耿平[40]還對季節變化進行了研究,結果表明土壤呼吸通常在夏季最高而冬季最低。

對毛竹林土壤呼吸的空間變異研究開展得比較少,姜艷[30]對海拔進行了研究,結果表明,隨著海拔升高,土壤總呼吸、根系呼吸及微生物呼吸均減小。其他如緯度、經度、地形和演替年齡等對毛竹林土壤呼吸的影響均未見報道。

2.5 毛竹林土壤呼吸對干擾的響應

人類活動產生的各種干擾也會影響土壤呼吸。葉耿平[40]研究表明,整個生長季,集約經營毛竹的土壤呼吸速率高于粗放經營毛竹林。劉勝[37]對雪災對土壤呼吸的影響做了初步研究,結果表明重度受災林分土壤呼吸速率分別是中度受災林分的1.65倍,輕度受災林分的1.22倍。其他干擾如CO2濃度升高、降雨變化、底物的增減、氮沉降等對毛竹林土壤呼吸的影響都未見報道。

2.6 毛竹林土壤呼吸組分

姜艷[30]等把毛竹林土壤呼吸組分界定為土壤根系呼吸、土壤微生物呼吸及土壤動物呼吸,采用壕溝法排除根系、電棒+衛生球法排除土壤動物的方法得出是土壤異養呼吸(微生物呼吸、動物呼吸)日貢獻率為42.85% ~ 53.43%(其中微生物呼吸為35.31%~45.52%,動物呼吸為6.92%~8.42%)。范少輝[32]等同樣也采用壕溝法研究不同組分,其中異養呼吸、自養呼吸和凋落物呼吸分別占毛竹林地土壤年釋放CO2量的59.5%、28.3%和12.2%。王燕[33]等采用根切除法測定無根土壤呼吸速率。肖復明[39]等同樣用壕溝法測定了土壤根系呼吸、礦質土壤呼吸、凋落物呼吸。以上研究結果基本都采用壕溝法來區分土壤呼吸各組分,屬于實驗處理法中的一種,除了實驗處理法外,還可以用同位素法、推理和模型法來研究土壤呼吸的各組分,但目前為止未見有人在毛竹林土壤呼吸的研究中采用這些技術。

2.7 毛竹林土壤呼吸模型的構建

對毛竹林土壤呼吸的模型主要集中在溫度—呼吸模型上,姜艷[30]、王超[35]、葉耿平[40]、李正才[34]都用vant’Hoff模型擬合了土壤呼吸和溫度的關系,結果表明毛竹林土壤呼吸速率與土壤溫度之間均具有指數正相關關系。姜艷[30]還同時采用了Lloyd & Taylor模型,結果表明兩種模型均能很好地表達土壤呼吸與土壤溫度的相關性,相比vant’ Hoff模型,Lloyd & Taylor方程擬合能力更好,擬合度R2稍有增大。李正才[34]還用直線方程、二次方程、三次方程對土壤呼吸隨氣溫變化的趨勢進行了擬合,結果表明直線方程和vant’ Hoff模型擬合的效果較好,而二次方程、三次方程模擬的效果較差。

對于土壤呼吸與土壤含水量的關系模型,姜艷[30]用線性模型模擬兩者之間的關系,結果在不同海拔、不同土壤深度含水量兩者相關關系不盡相同。王超[35]用回歸分析計算了土壤呼吸速率與土壤含水量之間的關系,結果表明土壤呼吸速率與土壤含水量變化呈一定的負相關關系。葉耿平[40]對土壤CO2排放量與土壤含水量進行相關性分析,結果表明集約經營毛竹林和常綠闊葉林CO2排放通量與土壤含水量呈指數關系,粗放經營毛竹林CO2排放通量與含水量呈對數關系。

對于底物—呼吸模型研究得比較少,只有王超[35]用回歸分析計算了土壤呼吸與凋落物的關系,結果表明土壤呼吸速率與凋落物量呈顯著的正相關關系

3 展望

由于毛竹林特殊的無性繁殖方式,地下埋藏著豐富的竹鞭,每年有大量新筍長大成竹向周圍擴展,使毛竹林土壤呼吸的研究有一定的難度。今后在開展毛竹林土壤呼吸研究的過程中,筆者建議從以下幾方面加強:

(1)繼續加強土壤溫度對毛竹林土壤呼吸影響機制的研究,深入開展土壤濕度對土壤呼吸影響的研究,同時,對于影響土壤呼吸的其他因子,如土壤氧氣、土壤氮、土壤質地和土壤pH值等,也應開展進一步研究。

(2)毛竹林土壤呼吸的時空變化長期定位研究是必要的,從中可獲得土壤呼吸的年際變化。在許多生態系統中都觀察到了土壤呼吸具有顯著的年際變化,如草原、山毛櫸林、溫帶混交林、美國黃松林和人工林[40~44]。特別是毛竹林的生長有明顯的大小年,植物由于對年際間氣候變化和干擾產生響應而引起的生理變化也可能影響土壤呼吸的年際變化??臻g變化上,對于緯度、經度、地形和演替年齡等對毛竹林土壤呼吸的影響也需做進一步的研究。同時毛竹林是一個人類干擾較大的生態系統,因人類活動(包括多種經營措施,如施肥、覆蓋等)對土壤呼吸的影響相對明顯的。

(3)繼續優化數學模型。在經典經驗模型的基礎上,探討溫度—呼吸模型、濕度—呼吸模型、底物—呼吸模型和多因子模型,并根據不同的氣候帶提出各自的經驗參數。

(4)采用更多樣化的方法來研究土壤呼吸各組分??梢試L試使用同位素法,利用同位素自然豐度的差異或者在同位素標記實驗中,給植物或生態系統添加痕量同位素以形成不同的同位素源值來區分土壤呼吸的組分。

[1] Raich J W , Potter C S. Global patterns of carbon dioxide emissions from soils[J]. Glob Biogeochem Cyc,1995, 9(1):23-36.

[2] 劉紹輝,方精云. 土壤呼吸的影響因素及全球尺度下溫度的影響[J]. 生態學報,1997,17(5):469-476.

[3] Tans P P, Fung I Y, Takahashi T. Observational constraints on the global atmospheric CO2budget[J]. Science,1990, 247(49):1 431-1 438.

[4] 康博文,劉建軍,黨坤良,等. 秦嶺火地塘林區油松林土壤碳循環研究[J]. 應用生態學報,2006,17(5):759-764.

[5] Raich J W, Tufekeioglu A. Vegetation and soil respiration: Correlation and Controls[J]. Biogeochemistry,2000, 48(1):71-90.

[6] Raich J W, Sehlesinger W H. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate[J]. Tellus, 1992, 44(2):81-99

[7] 陳述悅,李俊,陸佩玲,等. 華北平原麥田土壤呼吸特征[J]. 應用生態學報,2004,15(9):152-156.

[8] Wollny E. Undersuchungen uber den Einfluss der physikalischen Eigenschaften des Bodens auf dessen Gehalt an freier Kohlensaure[J]. ForGebiete-Agr Phys, 1831(4):1-28.

[9] Boussingault J B, Levy B. Me moire sur la composition del’aire confine dans la terre vegetable[J]. Ann Chem Phys, 1853(37):5-50.

[10] Moller. Uber die freie Kohlensaure im Boden[J]. For Gebiete-Agr Physiol, 1879(2):329-338.

[11] Russell E J, Appleyard A. The atmosphere of the soil its composition and the causes of variation[J]. J Agr Sci, 1915(7):1-48.

[12] Lundegardh H. Carbon dioxide evolution of soil and crop growth[J]. Soil Sci, 1927(23):417-453.

[13] Humfeld H. A method for measuring carbon dioxide evolution from soil[J]. Soil Sci, 1930,30(1):1-11.

[14] Greaves J R, Carter E G. Influence of moisture on the bacterial activities of the soil[J]. Soil Scie, 1920(10):361-387.

[15] Turpin H W. The carbon dioxide of the soil air[J]. Cornell Univ Agr Exp Stat Mem, 1920(32):319-362.

[16] Wiant H v. Influence of temperature on the rate of soil respiration[J]. J For, 1967(65):489-490.

[17] Kucera C L, Kirkham. Soil respiration studies in tallgrass prairie in Missouri[J]. Ecology, 1971(52):912-915.

[18] Coleman D C. Compartmental analysis of “total soil respiration”, an exploratory study[J]. Oikos, 1973(24):361-366.

[19] Edwards N T , Sollins. Continuous measurement of carbon dioxide evolution from partitioned forest floor components[J]. Ecology, 1973(54):406-412.

[20] Norman J M, Garcia R, Verma S B. Soil surface CO2fluxes and the carbon budget of a grassland[J]. J Geophys Res Atmos, 1992, 97(D17):18 845-18 853.

[21] Vanhala P. Seasonal variation in the soil respiration rate in coniferous forest soils[J]. Soil Biol Biochem, 2002(34):1 375-1379.

[22] Christopher W H, John M. B, Philip A F, et al. Increased rainfall variability and reduced rainfall amount decreases soil CO2flux in a grassland ecosystem[J].Glob Chang Biol, 2005, 11(2):322-334.

[23] Nay S M, Bormann B T. Soil carbon changes: Comparing flux monitoring and mass balance in a box lysimeter experiment[J]. Soil Sci Soc Am J, 2000(64):943-948.

[24] Autsin A T, Osvaldo E S. Carbon and nitrogen dynamics across a natural precipitation gradient in Patagonia,Argentina[J]. J Veg Sci, 2002, 13(3):351-360.

[25] Zhou X, Sherry R A, Yuan An, et al. Main and interactive effects of warming, clipping, and doubled precipitation on soil CO2efflux in a grassland ecosystem[J]. Glob Biogeocheml Cyc, 2006, 20,GB1003, doi: 10.1029/2005 GB002526.

[26] Vose J M, Ryan M G. Seanonal respiration of foliage, fine roots, and woody tissue in relation to growth, tissue N, and photosynthesis[J]. Glob Chang Biol, 2002, 8(2):182-193.

[27] Cheng W, Johnson D W. Elevated CO2, rhizosphere processes,and soil organic matter decomposition[J].Plant Soil, 1998, 202(2):167-174.

[28] Rainch J W, Potter C S, Bhagawati D. Interannual variability in global soil respiration,1980-94[J]. Glob Chang Biol, 2002, 8(8):800-812.

[29] Frank A B, Liebig M A, Hanson J D. Soil carbon dioxide fluxes in northern semiarid grasslands[J]. Soil Biol Biochem, 2002, 34(9):1 235 -1 241.

[30] 姜艷. 毛竹林土攘呼吸及其三個生物學過程的時空格局變化研究[D]. 北京:中國林業科學研究院,2010.

[31] 陳先剛,張一平,張小全,等. 過去50年中國竹林碳儲量變化[J]. 生態學報,2008,28(11):5 218-5 227.

[32] 范少輝,肖復明,汪思龍,等. 湖南會同林區毛竹林地的土壤呼吸[J]. 生態學報,2009,29(11):5 971-5 977.

[33] 王燕. 江西大崗山毛竹林生態系統碳平衡研究[D]. 北京:中國林業科學研究院,2008.

[34] 李正才. 土地利用變化對土壤有機碳的影響[D]. 北京:中國林業科學研究院,2006.

[35] 王超,楊智杰,等. 萬木林保護區毛竹林土壤呼吸特征及影響因素[J]. 應用生態學報,2011,22(5):1 212-1 218.

[36] 李雅紅,江洪,原煥英,等. 西天目山毛竹林土壤呼吸特征及其影響因子[J]. 生態學報,2010,30(17):4 590-4 597.

[37] 劉勝. 雪災對毛竹林土壤呼吸和土壤微生物量的影響[D]. 南京:南京林業大學,2010.

[38] 劉源月,江洪,邱忠平,等. 亞熱帶典型森林生態系統土壤呼吸[J]. 西南交通大學學報,2009,44(4):590-594.

[39] 肖復明. 毛竹林生態系統碳平衡特征的研究[D]. 北京:中國林業科學研究院,2007.

[40] 葉耿平,劉娟,姜培坤,等. 集約經營措施對毛竹林生長季土壤呼吸的影響[J]. 浙江農林大學學報,2011,28(1):18-25.

[41] Epron D, Nouvellon Y, Roupsard O, et al. Spatial and temporal variations of soil respiration in a Eucalyptus plantation in Congo[J]. For Ecol Manag, 2004, 202(1~3):149-160.

[42] Wickland K P, Striegl R G, Mast M A, et al. Interannual variation of soil respiration in two New England forests[J]. Glob Biogeochem Cycs, 2001, 15(2):337-350.

[43] Irvine J, Law B E. Contrasting soil respiration in young and old-growth ponderosa pine forests[J]. Glob Chang Biol, 2002, 8(12):11 83 -1 194.

[44] King J S, Hanson P J, Bernhardt E, et al. A multiyear synthesis of soil respiration responses to elevated atmospheric CO2from four forest face experiments[J]. Glob Chang Biol, 2004, 10(6):1 027-1 042.

Advance of Soil Respiration in Phyllostachys heterocycla cv. pubescens Stands

LEI Hai-qing1,2,FU Mao-Yi1*,LI Zheng-Cai1,TONG Long1
(1. Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Fuyang 311400, China;
2. Zhejiang Subtropica Crop Institute, Wenzhou 325005, China)

The paper reviewed research progress on soil respiration in Phyllostachys heterocycla cv. pubescens stands, including research methods, affect factors, spatial-temporal variation, response to disturbing, component division, model building. Some fields of study on soil respiration in Ph. heterocycla cv. pubescens stands need to be strengthened, such as effect of soil temperature, long-term positioning study on spatial-temporal variation, improvement of mathematic models, respiration component division with different methods.

soil respiration; Phyllostachys heterocycla cv. pubescens stands; research; advance

S718.5

A

1001-3776(2012)02-0071-05

2011-12-23;

2012-02-20

雷海清(1977-),女,云南曲靖人,工程師,博士研究生,從事森林資源可持續經營和管理研究;*通訊作者。

猜你喜歡
模型研究
一半模型
FMS與YBT相關性的實證研究
2020年國內翻譯研究述評
遼代千人邑研究述論
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
視錯覺在平面設計中的應用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
EMA伺服控制系統研究
新版C-NCAP側面碰撞假人損傷研究
3D打印中的模型分割與打包
主站蜘蛛池模板: 亚洲—日韩aV在线| 亚洲欧美在线精品一区二区| 久久精品国产91久久综合麻豆自制| 青草视频网站在线观看| 国产日本一区二区三区| 狠狠色婷婷丁香综合久久韩国| 成年人免费国产视频| 日本午夜三级| 五月婷婷综合色| 欧美色亚洲| 午夜福利无码一区二区| 国产极品美女在线播放| 伊人91视频| 国产网友愉拍精品| 高清视频一区| 欧美人与牲动交a欧美精品| 国产成人三级在线观看视频| 毛片最新网址| 91精品久久久无码中文字幕vr| 蜜芽一区二区国产精品| 欧美日本不卡| 天天综合网色| av大片在线无码免费| 国产日韩久久久久无码精品| 精品三级在线| 一级毛片不卡片免费观看| 这里只有精品在线播放| 国产浮力第一页永久地址| 欧美成人怡春院在线激情| 国产人人射| 免费一级毛片在线观看| 啪啪啪亚洲无码| 亚洲一区二区视频在线观看| 精品国产aⅴ一区二区三区| 午夜毛片免费观看视频 | 亚洲免费黄色网| 性69交片免费看| 最新精品久久精品| 国产精品手机在线播放| 欧美日韩中文字幕二区三区| 国产玖玖视频| 成年免费在线观看| 国产香蕉在线视频| 99久久精品免费看国产电影| 久久亚洲美女精品国产精品| 伊人久综合| 在线观看欧美国产| 国内精品久久久久久久久久影视 | 亚洲无码在线午夜电影| 亚洲视频二| 色综合天天娱乐综合网| 日韩人妻精品一区| 国产资源站| 综合社区亚洲熟妇p| 国产精品性| 丁香六月综合网| 亚洲精品天堂在线观看| 国内精品久久人妻无码大片高| 五月天天天色| 亚洲欧美不卡中文字幕| 亚洲国产精品成人久久综合影院| 亚洲网综合| 美女潮喷出白浆在线观看视频| 亚洲成年网站在线观看| 久久精品日日躁夜夜躁欧美| 色135综合网| 亚洲色图欧美激情| 日韩av无码DVD| 亚洲自偷自拍另类小说| 97视频在线精品国自产拍| 综合色天天| 欧美综合区自拍亚洲综合天堂| 秋霞国产在线| 青青操国产| 国产精品主播| 欧美色图久久| JIZZ亚洲国产| 国产最新无码专区在线| 国产精品jizz在线观看软件| 日韩中文字幕亚洲无线码| www.精品视频| 国产在线无码av完整版在线观看|