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

    植被對土壤熱擴散特征的影響
    ——以長白山闊葉紅松林為例

    2015-02-06 03:52:34施婷婷鄭興波張麗波李輝東李偉莉高玉芳
    生態(tài)學(xué)報 2015年12期
    關(guān)鍵詞:深度

    施婷婷, 鄭興波, 張麗波, 楊 弘, 李輝東, 李偉莉, 高玉芳

    1 南京信息工程大學(xué), 應(yīng)用氣象學(xué)院, 生態(tài)氣象環(huán)境研究中心, 南京 210044 2 中國科學(xué)院沈陽應(yīng)用生態(tài)研究所, 沈陽 110016 3 沈陽市水利建筑勘測設(shè)計院, 沈陽 110015

    植被對土壤熱擴散特征的影響
    ——以長白山闊葉紅松林為例

    施婷婷1, 鄭興波2,*, 張麗波1, 楊 弘2, 李輝東2, 李偉莉3, 高玉芳1

    1 南京信息工程大學(xué), 應(yīng)用氣象學(xué)院, 生態(tài)氣象環(huán)境研究中心, 南京 210044 2 中國科學(xué)院沈陽應(yīng)用生態(tài)研究所, 沈陽 110016 3 沈陽市水利建筑勘測設(shè)計院, 沈陽 110015

    土壤溫度變化及熱傳遞是影響土壤和大氣水熱交換的重要過程,而植被是決定這種變化和影響的環(huán)境因子之一。通過比較林地與裸地土壤熱特性的差異分析植被對土壤熱擴散的影響。 研究內(nèi)容包括觀測2007年長白山闊葉紅松林區(qū)不同深度的土壤溫度以及同期土壤含水量和葉面積指數(shù)。分析不同深度林地與裸地土壤溫度年周期特征,根據(jù)熱傳導(dǎo)方程估算各土壤層溫度的阻尼深度和熱擴散率,并探討引起林地與裸地土壤熱特性差異的可能原因。結(jié)果表明,林地與裸地的土壤溫度有明顯的時空變化規(guī)律。隨著深度的增加,土壤溫度年周期的振幅逐漸減小、相位逐漸增大、平均值逐漸升高。林地土壤溫度年平均值低于裸地,表層年平均溫度相差約0.8℃;地表以下相同深度處,林地土壤溫度年周期的振幅約低于裸地2.6—2.9℃,相位約小于裸地0.2—0.24 rad(角速度),這表明林地土壤溫度極大值和極小值出現(xiàn)的時間比裸地滯后約11—14d。土壤溫度阻尼深度和熱擴散率隨深度的增加而逐漸增大,而在1.6—3.2m則略有降低的趨勢。林地與裸地土壤溫度和熱特性的時空特征和差異可能與土壤含水量和葉面積指數(shù)有關(guān)。

    土壤溫度; 熱擴散; 闊葉紅松林; 阻尼深度

    The influence of vegetation on soil thermal properties: a case study of

    土壤作為植物生長的媒介,為植物提供水分和礦物營養(yǎng)[1],土壤熱環(huán)境的變化,必然伴隨著土壤微生物和酶的活性改變,影響著土壤質(zhì)量,對植物的生長發(fā)育乃至群落的演替和生態(tài)系統(tǒng)的恢復(fù)均發(fā)揮著直接而深遠的影響。土壤溫度是計算地表與大氣之間輻射與能量交換的重要參數(shù),并且在一定程度上決定了地表的可用能量向能量收支各項的分配,包括顯熱、潛熱、土壤熱通量等[2]。傳統(tǒng)的土壤熱量傳輸過程主要考慮熱傳導(dǎo)過程,Bhumralkar[3]以正弦函數(shù)方程作為邊界條件,推導(dǎo)出地表以下給定深度的土壤溫度。Gao[4]則在此基礎(chǔ)上還考慮了土壤水分的垂直流動對土壤熱量傳遞的貢獻,以包含熱傳導(dǎo)項和熱對流項的土壤熱量傳輸方程來描述一個半干旱草地的土壤熱量交換過程。地表以下的土壤溫度及其垂直梯度將影響植物根部溫度、植被生長速率、初級生產(chǎn)力、分解以及蒸散速率[5],進而對土壤中的熱量傳遞和熱量儲存產(chǎn)生較強的驅(qū)動作用[6]。在植被覆蓋區(qū),遮蔭和植被的熱量傳輸過程較大程度的控制了植被與大氣之間的能量交換,因此可大幅降低土壤表面的溫度日變化[7]。自1980年代開始,植被對土壤熱量傳遞的影響逐漸被學(xué)者重視,并被嵌入陸地表面過程模式中,如Sellers等[8]的SiB2模型、Verseghy[9]的CLASS模型、Wang和Leuning[10]的CABLE模型、Dai和Zeng[11]的IAP94模型、Desborough和Pitman[12]的BASS模型等。植被類型對氣候和其它環(huán)境因子的反應(yīng)是綜合的,主要涉及熱量和降水[13-14]兩大要素。而植被對局地氣候狀況也有較強的反饋作用,如植被覆蓋下的氣溫和土壤溫度的振幅均小于裸地[15-16],植被覆蓋對土壤熱擴散特征影響的研究仍然較為缺乏。本文主要分析了林地與裸地不同深度的土壤溫度變化和熱擴散特征,并探究引起林地與裸地土壤熱特性差異的可能影響因子,不僅為理解植被覆蓋下的土壤熱通量的時空變化提供理論支持,也為基于過程的陸面過程模式的土壤熱擴散過程提供一定的數(shù)據(jù)支持。

    1 研究地區(qū)與方法

    1.1 研究地概況

    本研究在吉林省長白山自然保護區(qū)的中國科學(xué)院長白山森林生態(tài)系統(tǒng)定位站及其西邊約1km的闊葉紅松林內(nèi)1號標(biāo)準(zhǔn)地進行(42°24′N,128°6′E,海拔738m)。林內(nèi)主喬木為紅松(Pinuskoraiensis)、椴樹(Tiliaamurensis)、水曲柳(Fraxinusmandshurica)、色木槭(Acermono)、蒙古櫟(Quercusmongolica)等,林分為復(fù)層結(jié)構(gòu),平均株高約26m,立木株數(shù)約560株/hm2,郁閉度約0.8,下木覆蓋度40%[17]。闊葉紅松林帶土壤為白漿化的山地暗棕色森林土,坡降2%—4%,地表0—5cm為枯枝落葉層,5—11cm為深灰色或深灰棕色腐殖質(zhì)層,粘壤土,疏松,土壤顏色向下過渡明顯,根系分布很多,11—32cm為白漿化暗棕色粘土,根系分布較多,但向下明顯減少,32—105cm為為暗棕色粘土,核塊狀結(jié)構(gòu),基本沒有根系分布[18]。闊葉紅松林帶的年平均氣溫約為4.5℃,最冷月為1月(平均氣溫約-12.9℃),最熱月為7月(平均氣溫約為19.4℃),日平均溫度最高和最低值約為24.0℃和-21.8℃。

    在闊葉紅松林內(nèi)建有高約62m的為氣象觀測塔,在塔上32m高度安裝了光合有效輻射探頭(Li- 190Sb, LiCor Inc., USA),在2.5 m高度安裝了5個桿狀光合有效輻射探頭(LQS70- 10, APOGEE, USA),分別對冠層上方和下方的光合有效輻射進行測定。在氣象塔的地表以下5個深度(0、0.05、0.2、0.5、1m)設(shè)置兩套土壤溫度傳感器(105T和107T,Campbell,USA),并與數(shù)據(jù)采集器(CR23X,Campbell,USA)相連,每30分鐘記錄平均值,并計算日平均溫度。在長白山森林生態(tài)系統(tǒng)定位站內(nèi)設(shè)有標(biāo)準(zhǔn)氣象站自動觀測地溫,觀測系統(tǒng)為M520,根據(jù)地面氣象觀測規(guī)范,氣象站土壤溫度在面積為1m2左右的裸地下觀測,對土壤內(nèi)部5個深度(地表以下0、0.4、0.8、1.6、3.2m)進行日平均溫度的觀測和記錄。本文中數(shù)據(jù)采集時間為2007年1月1日至2007年12月31日。

    從5月1日至10月31日在闊葉紅松林內(nèi)對葉面積指數(shù)(LAI)進行人工觀測,采用葉面積分析儀(LAI- 2000, LiCor Inc., Lincoln, NE, USA)在林內(nèi)選取10個樣點進行觀測,觀測高度為1m。由于LAI是通過植被冠層的輻射傳輸模型來計算的,要求天空為散射輻射,因而選擇早晨日出前或陰天散射輻射較多時進行觀測。在生長季初末期(5月、9月和10月)葉片剛開始展葉或枯黃凋萎時,葉面積變化較大,每隔3d觀測1次。在生長旺盛期(6—8月)葉面積較為穩(wěn)定時,每隔5—10d觀測1次。

    1.2 研究方法

    對于水平均質(zhì)的土壤,設(shè)其熱量僅在垂直方向上交換,則土壤的熱傳導(dǎo)方程為[3]:

    (1)

    式中,Ts為土壤溫度(℃),t為時間(s),z為土壤深度(m),k為土壤的熱擴散率(m2/s),k=λ/c,λ為導(dǎo)熱率(W m-1K-1),c為土壤的容積熱容量(J m-3K-1))。給定土壤表面溫度變化的邊界條件:

    (2)

    (3)

    ΔTz=ΔT0·e-z/d

    (4)

    φz=φ0+z/d

    (5)

    由于闊葉紅松林的裸地和林內(nèi)各有五層不同深度的土壤溫度觀測值,根據(jù)相鄰兩個不同深度的土壤溫度資料可分別計算出該厚度土壤溫度年尺度的阻尼深度d和熱擴散率k。

    2 結(jié)果分析

    2.1 林地與裸地土壤溫度的年變化特征

    圖1 長白山闊葉紅松林裸地(a)與林內(nèi)(b)地下不同深度的土壤溫度Ts年變化的觀測值(點)與擬合曲線(線)

    長白山闊葉紅松林地區(qū)裸地的土壤溫度具有明顯的時空變化規(guī)律(圖1a)。在垂直方向上,冬季(12月—次年2月)隨著深度的增加,土壤溫度逐漸升高,3.2m深度處的土壤溫度比土壤表面溫度月平均值可高達22.3℃。夏季(6月—8月)隨著深度的增加,土壤溫度逐漸降低,3.2m處深度處的土壤溫度比土壤表面溫度月平均值可低至19.0℃。春秋季為溫度隨深度變化的過渡季。在年周期上,土壤表層溫度的年變化最明顯,冬季溫度較低,最冷月為1月,溫度為-15.1℃,日平均溫度最低為-23.7℃(1月13日)。隨著春季的來臨,氣溫逐漸升高,土壤溫度也隨之升高,其表層最熱月出現(xiàn)在7月,溫度為24.5℃,日平均溫度最高為33℃(7月26日)。當(dāng)進入秋季,隨著氣溫的降低,土壤溫度也逐漸降低,土壤表層溫度在年周期上呈單峰形。較深層土壤溫度的年周期與表層相似,但年最高溫與年最低溫的差異比表層小,且年最高溫和最低溫出現(xiàn)的日期遲于表層。0.4、0.8、1.6m和3.2m深度處的最高溫度及其出現(xiàn)日期分別為20.6℃(8月17日)、18.6℃(8月20日)、15.8℃(8月25日)、11.4℃(9月22日),而最低溫度及其出現(xiàn)日期分別為-2.9℃(1月29日)、-0.5℃(2月5日)、1.1℃(4月6日)、4℃(5月8日)。這表明,與土壤表層溫度年變化特征相比,隨著深度的增加,土壤溫度年最高值逐漸減小,而年最低值逐漸增大,且極大值和極小值出現(xiàn)的時間依次滯后。

    闊葉紅松林內(nèi)土壤溫度的時空變化規(guī)律與裸地相似但又有一定的差異。在年周期上,林地土壤溫度也呈單峰形變化,在空間上,與表層土壤溫度相比,隨著土壤深度的增加,土壤溫度在冬季逐漸升高,而在夏季逐漸降低,且極大值和極小值出現(xiàn)的時間依次滯后(圖1b)。林地1m處的土壤溫度月平均值與表層溫度差值可達-6.6℃(6月)至6.4℃(1月)。林地土壤表層溫度年平均值低于裸地,相差約0.8℃,且林地土壤表層溫度的年較差也明顯低于裸地,前者為21.8℃,后者為39.6℃,這表明地上植被和地表枯枝落葉層對土壤溫度的影響較大。

    2.2 林地與裸地土壤溫度平均值、振幅與相位的垂直廓線

    圖2 長白山闊葉紅松林裸地與林內(nèi)地下不同深度的土壤溫度Ts年周期的振幅ΔT、相位φ和平均值

    Table 1 The annual temperature amplitude ΔTz, phaseφz, average value Tzand determination coefficientsR2in different depthsz(m)of the bare land and broadleaved Korean pine forest, estimated with equation (3)

    地表條件Surfacecondition深度Depthz/mΔTz/℃φz/radTz/℃R2裸地Bareland019.85±0.271.87±0.016.23±0.190.940.411.2±0.122.15±0.016.97±0.080.960.89.44±0.112.36±0.017.1±0.080.951.67.15±0.072.69±0.017.39±0.050.963.23.15±0.033.43±0.017.39±0.020.96林內(nèi)Forestland010.91±0.112.09±0.015.59±0.080.960.0510.47±0.112.14±0.015.6±0.080.960.29.32±0.12.27±0.015.77±0.070.960.57.88±0.082.44±0.015.88±0.060.9616.17±0.062.65±0.016.06±0.040.97

    2.3 林地與裸地土壤溫度的阻尼深度與熱擴散率特征

    根據(jù)方程(3)可得出相應(yīng)的振幅ΔTz和相位φz,根據(jù)振幅和相位與阻尼深度d和熱擴散率k之間的關(guān)系,可分別由方程(4)(即振幅法)和方程(5)(即相位法)反推出d和k。結(jié)果表明,這兩種方法計算出的阻尼深度和熱擴散率在林地和裸地都表現(xiàn)出較明顯的垂直變化特征(表2)。從地表至1.6m深度的土壤層中,阻尼深度和熱擴散率隨著深度的增加而逐漸增大,而1.6—3.2m的土壤層內(nèi),阻尼深度和熱擴散率則略有降低的趨勢。在裸地,由振幅法得出的熱擴散率從0—0.4m土壤層的0.49×10-7m2/s逐漸升高至0.8—1.6m土壤層的8.25×10-7m2/s,而由相位法得出的熱擴散率由0—0.4m土壤層的2.06×10-7m2/s逐漸升高至0.8—1.6m土壤層的5.82×10-7m2/s。在闊葉紅松林內(nèi),由振幅法和相位法得出的熱擴散率較為接近,從表層土壤的1.04×10-7—1.46×10-7m2/s逐漸升高至1m處的4.18×10-7—5.64×10-7m2/s。由相位法和振幅法計算出的d和k沒有顯著差異,且在相同深度處,林地和裸地的d和k也無明顯的系統(tǒng)性差異性。

    表2 長白山闊葉紅松林裸地與林內(nèi)地下不同深度的土壤溫度在年尺度上的阻尼深度d和熱擴散率k、d和k分別由振幅法和相位法計算得出

    Table 2 The damping depthdand thermal diffusivitykat different depthszof the bare land and broadleaved Korean pine forest, estimated with amplitude- and phase-equation

    地表條件Surfacecondition深度Depthz/m由振幅法計算Estimatedwithamplitude-equationd/mk/(m2/s)由相位法計算Estimatedwithphase-equationd/mk/(m2/s)裸地Bareland 0—0.40.700.49×10-71.442.06×10-7 0.4—0.82.345.47×10-71.833.32×10-7 0.8—1.62.878.25×10-72.425.82×10-7 1.6—3.21.953.79×10-72.184.75×10-7林內(nèi)Forestland 0—0.051.211.46×10-71.021.04×10-70.05—0.21.291.66×10-71.161.35×10-7 0.2—0.51.783.17×10-71.743.03×10-7 0.5—12.044.18×10-72.385.64×10-7

    3 討論

    3.1 植被覆蓋對土壤溫度的影響

    由于植被對太陽輻射的反射和吸收作用、以及植被的蒸騰作用等,植被覆蓋對局地氣溫和地表輻射有較大的影響。與裸地相比,林地提高了夜間氣溫,降低了白晝氣溫,使氣溫日振幅變小[15],主要原因有兩方面,一是太陽輻射經(jīng)過林冠的反射與吸收,到達地表時輻射量已大幅降低,在晝間林內(nèi)輻射差額的正值有減小的效應(yīng),而在夜間林內(nèi)輻射差額的負值亦有減小的效應(yīng),因而氣溫日振幅減??;二是林冠降低風(fēng)速和湍流傳輸,阻隔水平方向和垂直方向的熱量與水汽交換,同時反射地面長波輻射,使得林內(nèi)溫度升降幅度趨緩,導(dǎo)致林地與裸地?zé)崃繝顩r差異較大。土壤內(nèi)部溫度變化的驅(qū)動力主要來自于土壤表層接受到的輻射和熱量。對于特定的植物種類,地表能接受到的太陽輻射主要受葉面積指數(shù)LAI的影響,而葉面積指數(shù)同時也反映了其蒸騰能力的大小,因而葉面積指數(shù)對土壤溫度的影響較大。裸地與林地不同深度處(0.4、0.5、0.8m和1m)土壤溫度的差異(裸地土壤溫度減去林地土壤溫度)都具有明顯的季節(jié)變化規(guī)律(圖3)。在冬季(11月至翌年2月),裸地土壤溫度比林地土壤溫度約低0.8—2.8℃。隨著春季的到來,植被開始展葉,葉面積指數(shù)從5月份開始逐漸增加(圖4),植被對土壤表面接受到的輻射和熱量的影響越來越大,林地土壤溫度的日平均值越來越低于裸地,至7月份差值達到最大,各深度處林地土壤溫度月平均值比裸地低約4.5—5℃。隨著秋季枝葉枯黃凋落,葉面積指數(shù)逐漸減小,植被對土壤表面輻射和熱量的影響也降低,林地和裸地土壤溫度的差值也逐漸減小,在10月份左右該差值由正值變?yōu)樨撝?,表明到了冬季,林地土壤溫度高于裸地?/p>

    圖3 裸地與林地在0.4、0.5、0.8m和1m深度處土壤溫度的差值,縱坐標(biāo)為裸地土壤溫度減去林地土壤溫度

    圖4 闊葉紅松林葉面積指數(shù)LAI的季節(jié)變化

    由于裸地與林地對太陽輻射的反射率不同,土壤溫度也有一定差異。植被覆蓋對不同季節(jié)土壤溫度的影響機制,已被廣泛應(yīng)用于農(nóng)業(yè)生產(chǎn)。冬季,積雪的土面反射率較大,達75%—95%,植被冠層的反射率則較小,約為5%—25%,用草覆蓋地面可減小反射率,有利于提高土壤溫度。如,在我國北方寒冷地區(qū),通過地面覆蓋,提高春季土壤溫度,活化土壤微生物和酶的活性,提高土壤質(zhì)量,延長作物的生育期[5]。

    3.2 土壤含水量對熱擴散率垂直梯度的的影響

    土壤是由3種物質(zhì)(固體、液體、氣體)組成的介質(zhì),熱量的傳遞就通過固體微粒、水和空氣進行,土壤導(dǎo)熱性能的高低取決于各種微粒相互接觸的形態(tài)、性質(zhì)以及土壤含水量的多少[20]。對于特定的研究地點,土壤容重和礦物組成變化較小,因而土壤熱擴散率主要受到土壤含水(含冰)量和溫度的影響。由于冰與水的導(dǎo)熱率有較大差異,因而土壤在凍結(jié)和融化狀態(tài)下的熱擴散率主導(dǎo)因子不同。在融化狀態(tài)下,土壤熱擴散率主要與土壤含水量有關(guān)。在凍結(jié)狀態(tài)下,土壤中總含水量變化不大,土壤熱擴散率主要與溫度和固液態(tài)水的比例有關(guān)[21]。因而,本研究選取5—10月的數(shù)據(jù)來討論土壤含水量對熱擴散率垂直梯度的影響(圖5),該期間內(nèi)土壤中無凍結(jié)現(xiàn)象,處于融化狀態(tài)。土壤的熱擴散率k隨著深度的增加而呈逐漸升高,利用SPSS軟件進行T檢驗,發(fā)現(xiàn)熱擴散率隨深度而升高的趨勢已達到顯著性水平(Sig.<0.05),土壤含水量隨深度而增加的趨勢也達到顯著性水平(Sig.<0.05)。在土壤的三相介質(zhì)中,固相含量相對變化較小,而水的熱容量(4.19×106J m-3K-1)與空氣的熱容量(0.0012×106J m-3K-1)有較大差異,水的導(dǎo)熱率(0.58 W m-1K-1)約為空氣導(dǎo)熱率(0.026 W m-1K-1)的22倍[22],因而土壤熱擴散率受其含水量的影響較大,Gao[4]的研究表明土壤熱擴散率隨著土壤含水量的升高而增加。這與本研究中土壤熱擴散率與土壤含水量呈正相關(guān)關(guān)系的結(jié)果相似。因而土壤物理學(xué)的研究中越來越重視土壤水分運移對熱量傳遞的影響,如Shao等[23]在分析土壤熱擴散方程時同時考慮了土壤中的熱傳導(dǎo)和熱對流,后者即為由于土壤水分的垂直運動而對土壤溫度產(chǎn)生的貢獻。

    圖5 闊葉紅松林內(nèi)地表以下0.05、0.2m和0.5m深度處土壤熱擴散率和土壤含水量5—10月的變化

    4 結(jié)論

    (1)植被覆蓋對土壤溫度的年周期具有較大的影響,長白山闊葉紅松林地區(qū)與與裸地的土壤溫度年變化都與正弦曲線相似,表層溫度的極大值和極小值分別出現(xiàn)在7月和1月,隨著深度的增加,極值出現(xiàn)的時間依次滯后。但林地土壤溫度的年平均值低于裸地,表層相差約0.8℃,0.4—1m的土壤層相差約1.1℃;同時,林地土壤溫度年較差也明顯小于裸地。

    (2)林地與裸地土壤溫度年周期均有明顯的垂直變化特征,同時也受到植被狀況的影響。隨著深度的增加,土壤溫度年周期的振幅逐漸減小、相位逐漸增大、平均值逐漸升高。而在地表以下相同深度處,林地土壤振幅比裸地低2.6—2.9℃,相位比裸地小0.2—0.24 rad,意味著極值出現(xiàn)的時間滯后約11—14d。

    (3)根據(jù)振幅法和相位法推算的土壤溫度的阻尼深度和熱擴散率較為接近,從地表至1.6m深度的土壤層中,阻尼深度和熱擴散率隨著深度的增加而逐漸增大,而1.6—3.2m的土壤層內(nèi),阻尼深度和熱擴散率則略有降低的趨勢。植被覆蓋對阻尼深度與熱擴散率無顯著影響。

    (4)表征植被狀況指標(biāo)之一的葉面積指數(shù)對林地與裸地土壤溫度差異具有較為顯著的影響,相同深度處林地與裸地土壤溫度的差異隨著葉面積指數(shù)的升高(降低)而增加(減小)。土壤熱擴散率隨著深度的增加而增加,可能與不同深度處土壤含水量不同有關(guān)。

    致謝:中國科學(xué)院沈陽應(yīng)用生態(tài)研究所關(guān)德新研究員、南京信息工程大學(xué)高志球研究員、澳大利亞聯(lián)邦科學(xué)與工業(yè)研究組織王應(yīng)平老師和美國農(nóng)業(yè)部林業(yè)局森林變化研究中心劉永強老師,對本文給予幫助,特此致謝。

    [1] 劉占鋒, 傅伯杰, 劉國華, 朱永官. 土壤質(zhì)量與土壤質(zhì)量指標(biāo)及其評價. 生態(tài)學(xué)報, 2006, 26(3): 901- 913.

    [2] Bateni S M, Entekhabi D. Relative efficiency of land surface energy balance components. Water Resources Research, 2012, 48(4), W04510, doi: 10.1029/2011WR011357.

    [3] Bhumralkar C M. Numerical experiments on the computation of ground surface temperature in an atmospheric general circulation model. Journal of Applied Meteorology, 1975, 14(7): 1246- 1258.

    [4] Gao Z Q. Determination of soil heat flux in a Tibetan short-grass prairie. Boundary-Layer Meteorology, 2005, 114(1): 165- 178.

    [5] 邵明安, 王全九, 黃明斌. 土壤物理學(xué). 北京: 高等教育出版社, 2006: 162- 166.

    [6] Hirota T, Pomeroy J W, Granger R J, Maule C P. An extension of the force-restore method to estimating soil temperature at depth and evaluation for frozen soils under snow. Journal of Geophysical Research, 2002, 107(D24): ACL 11- 1-ACL 11- 10, doi: 10.1029/2001JD001280.

    [7] Dickinson R E. The force-restore model for surface temperatures and its generalizations. Journal of Climate, 1988, 1(11): 1086- 1097.

    [8] Sellers P J, Tucker C J, Collatz G J, Los S O, Justice C O, Dazlich D A, Randall D A. A revised land surface parameterization (SiB2) for atmospheric GCMs. Part Ⅱ: The generation of global fields of terrestrial biophysical parameters from satellite data. Journal of Climate, 1996, 9(4): 706- 737.

    [9] Verseghy D L. CLASS—a Canadian land surface scheme for GCMs. I. Soil model. International Journal of Climatology, 1991, 11(2): 111- 133.

    [10] Wang Y P, Leuning R. A two-leaf model for canopy conductance, photosynthesis and partitioning of available energy I: Model description and comparison with a multi-layered model. Agricultural and Forest Meteorology, 1998, 91(1/2): 89- 111.

    [11] Dai Y J, Zeng Q C. A Land Surface Model (IAP94) for Climate Studies Part I: Formulation and Validation in Off-line Experiments. Advances in Atmospheric Sciences, 1997, 14(4): 433- 460.

    [12] Desborough C E, Pitman A J. The BASE land surface model. Global and Planetary Change, 1998, 19(1/4): 3- 18.

    [13] 方精云. 我國森林植被帶的生態(tài)氣候?qū)W分析. 生態(tài)學(xué)報, 1991, 11(4): 377- 387.

    [14] 周廣勝, 張新時. 中國氣候-植被關(guān)系初探. 植物生態(tài)學(xué)報, 1996, 20(2): 113- 119.

    [15] 吳家兵, 關(guān)德新, 代力民, 金明淑, 徐浩. 長白山闊葉紅松林夏季溫度特征研究. 生態(tài)學(xué)雜志, 2002, 21(5): 14- 17.

    [16] 孫金偉, 吳家兵, 關(guān)德新, 王安志, 袁鳳輝, 張新建. 森林與空曠地空氣溫濕度及土壤溫度的長期對比研究. 生態(tài)學(xué)雜志, 2011, 30(12): 2685- 2691.

    [17] 王戰(zhàn), 徐振邦, 李昕. 長白山北坡主要森林類型及其群落結(jié)構(gòu)特點(之一) // 中國科學(xué)院長白山森林生態(tài)系統(tǒng)定位站. 森林生態(tài)系統(tǒng)研究. 北京: 林業(yè)出版社, 1980, 1: 25- 42.

    [18] 程伯容, 許廣山, 丁桂芳, 張玉華. 長白山北坡自然保護區(qū)主要土壤類型及其基本特性 // 中國科學(xué)院長白山森林生態(tài)系統(tǒng)定位站. 森林生態(tài)系統(tǒng)研究. 北京: 林業(yè)出版社, 1981, 2: 196- 206.

    [19] 關(guān)德新, 吳家兵, 王安志, 金昌杰, 施婷婷, 韓士杰. 長白山紅松針闊葉混交林林冠層葉面積指數(shù)模擬分析. 應(yīng)用生態(tài)學(xué)報, 2007, 18(3): 499- 503.

    [20] 翁篤鳴, 陳萬隆, 沈覺成. 小氣候和農(nóng)田小氣候. 北京: 農(nóng)業(yè)出版社, 1981: 30- 30.

    [21] 李韌, 趙林, 丁永建, 焦克勤, 王銀學(xué), 喬永平, 杜二計, 劉廣岳, 孫琳嬋, 肖瑤. 青藏高原北部活動層土壤熱力特性的研究. 地球物理學(xué)報, 2010, 53(5): 1060- 1072.

    [22] Oke T R. Boundary Layer Climates. 2nd ed. New York: Routledge, 1987: 42- 48.

    [23] Shao M A, Horton R, Jaynes D B. Analytical solution for one-dimensional heat conduction-convection equation. Soil Science Society of America Journal, 1998, 62(1): 123- 128.

    broadleaved Korean pine forest in Changbai Mountain

    SHI Tingting1, ZHENG Xingbo2,*, ZHANG Libo1, YANG Hong2, LI Huidong2, LI Weili3, GAO Yufang1

    1InternationalCenterforEcology,MeteorologyandEnvironment,SchoolofAppliedMeteorology,NanjingUniversityofInformationScienceandTechnology,Nanjing210044,China2InstituteofAppliedEcology,ChineseAcademyofSciences,Shenyang110016,China3ShenyangHydraulicConstructionSurveyandDesignInstitute,Shenyang110015,China

    Vegetation is one of the determining factors for soil temperature variation and thermal transfer, which are important processes for the energy and water exchange on the land-atmosphere interface. This study was to understand the impacts of vegetation on soil thermal diffusivity through comparing observations on forested and bared lands. The observations were made in broadleaved Korean pine forest over Changbai Mountain during 2007. The annual cycle of soil temperature and thermal properties at five depths were measured, together with soil volumetric water content and leaf area index.. The soil thermal diffusivity and damping depth of the soil temperature fluctuations were estimated using the heat conduction equation and the sinusoidal-fluctuating temperature at the ground surface. The potential environmental factors for the thermal properties were explored. The results show that soil temperature in both bared land and forest had similar spatial and temporal variations but with different magnitudes. With increasing soil depth, the amplitude of annual temperature variation decreased gradually, while both its phase and annual average temperature itself increased. In comparison with the based soil, the annual average temperature of the surface soil layer in the forest was 0.8℃ lower, the annual amplitude of soil temperature was 2.6—2.9℃ lower, and the temperature phase was 0.2—0.24 rad lower, suggesting that the maximum and minimum temperatures of the forest lagged by 11—14 days behind those of the bare land. The damping depth and thermal diffusivity increased with the soil depth, with an exception for the 1.6—3.2m layers where slightly decrease with depth was observed. It is indicated that soil moisture and leaf area index are probably responsible for the different characteristics in soil temperature and thermal properties between the forest and bared land.

    soil temperature; thermal diffusivity; broadleaved korean pine forest; damping depth

    國家自然科學(xué)基金項目(41105112, 41105077); 南京信息工程大學(xué)生態(tài)氣象環(huán)境研究中心科研基金; 江蘇省自然科學(xué)基金面上項目(BK2012859); 江蘇高校優(yōu)勢學(xué)科建設(shè)工程(PAPD)

    2013- 08- 17;

    2014- 07- 02

    10.5846/stxb201308172098

    *通訊作者Corresponding author.E-mail: carolstt@gmail.com

    施婷婷, 鄭興波, 張麗波, 楊弘, 李輝東, 李偉莉, 高玉芳.植被對土壤熱擴散特征的影響——以長白山闊葉紅松林為例.生態(tài)學(xué)報,2015,35(12):3970- 3978.

    Shi T T, Zheng X B, Zhang L B, Yang H, Li H D, Li W L, Gao Y F.The influence of vegetation on soil thermal properties: a case study of broadleaved Korean pine forest in Changbai Mountain.Acta Ecologica Sinica,2015,35(12):3970- 3978.

    猜你喜歡
    深度
    深度理解不等關(guān)系
    四增四減 深度推進
    深度理解一元一次方程
    深度觀察
    深度觀察
    深度觀察
    深度觀察
    芻議深度報道的深度與“文”度
    新聞傳播(2016年10期)2016-09-26 12:14:59
    提升深度報道量與質(zhì)
    新聞傳播(2015年10期)2015-07-18 11:05:40
    微小提議 深度思考
    精品久久国产蜜桃| 国产日本99.免费观看| 18禁黄网站禁片午夜丰满| 亚洲四区av| 久久国内精品自在自线图片| 久9热在线精品视频| 99久久精品国产国产毛片| 色精品久久人妻99蜜桃| 人妻少妇偷人精品九色| 成人性生交大片免费视频hd| 欧美成人免费av一区二区三区| 美女被艹到高潮喷水动态| 久久久久国产精品人妻aⅴ院| 久久精品91蜜桃| 久久国产乱子免费精品| 色哟哟·www| 18禁裸乳无遮挡免费网站照片| 麻豆av噜噜一区二区三区| av黄色大香蕉| 搡老熟女国产l中国老女人| 一级黄色大片毛片| 日韩欧美 国产精品| 身体一侧抽搐| 精品久久久久久久久亚洲 | 欧美国产日韩亚洲一区| 欧美日韩综合久久久久久 | 97超级碰碰碰精品色视频在线观看| 欧美潮喷喷水| 男人舔奶头视频| 床上黄色一级片| 亚洲四区av| 亚洲天堂国产精品一区在线| 亚洲电影在线观看av| 亚洲国产日韩欧美精品在线观看| 亚洲精品在线观看二区| 一级毛片久久久久久久久女| 国产精品久久久久久久久免| 又黄又爽又免费观看的视频| 五月伊人婷婷丁香| 99久久无色码亚洲精品果冻| 久久久久精品国产欧美久久久| 伦理电影大哥的女人| 五月伊人婷婷丁香| 欧美成人一区二区免费高清观看| 久久人妻av系列| 精品不卡国产一区二区三区| 亚洲精品色激情综合| 性色avwww在线观看| 日韩欧美国产在线观看| 精品一区二区三区人妻视频| 成年女人永久免费观看视频| 亚洲人成伊人成综合网2020| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆| 床上黄色一级片| 国产精品一区www在线观看 | 最近最新免费中文字幕在线| 日本一本二区三区精品| 亚洲专区中文字幕在线| 深爱激情五月婷婷| 内地一区二区视频在线| 动漫黄色视频在线观看| 国产男人的电影天堂91| 欧美色欧美亚洲另类二区| 免费av观看视频| 午夜福利在线在线| 两个人视频免费观看高清| 校园春色视频在线观看| 三级男女做爰猛烈吃奶摸视频| 热99在线观看视频| 国产亚洲av嫩草精品影院| 国产精品一区二区免费欧美| 欧美+亚洲+日韩+国产| 九九热线精品视视频播放| 三级毛片av免费| 波多野结衣巨乳人妻| 中文字幕av成人在线电影| 男女边吃奶边做爰视频| 色视频www国产| 国产免费男女视频| 男女做爰动态图高潮gif福利片| 国产精品久久久久久av不卡| 国产又黄又爽又无遮挡在线| 男女啪啪激烈高潮av片| 日本a在线网址| 婷婷丁香在线五月| 日本成人三级电影网站| 亚洲内射少妇av| 能在线免费观看的黄片| 伦精品一区二区三区| 国产真实伦视频高清在线观看 | 免费观看在线日韩| 亚洲无线在线观看| 婷婷丁香在线五月| 成年女人看的毛片在线观看| 久久精品综合一区二区三区| 如何舔出高潮| 全区人妻精品视频| 色综合站精品国产| 亚洲av不卡在线观看| 国产乱人伦免费视频| 国产精品av视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 久久人妻av系列| 1000部很黄的大片| 国产精品女同一区二区软件 | 深夜精品福利| 在线免费观看不下载黄p国产 | 久久久久国内视频| 国产激情偷乱视频一区二区| av国产免费在线观看| 欧美成人a在线观看| 天堂影院成人在线观看| 91午夜精品亚洲一区二区三区 | 99视频精品全部免费 在线| 精品乱码久久久久久99久播| 亚洲人与动物交配视频| 内射极品少妇av片p| 国产一区二区在线观看日韩| 俺也久久电影网| 成人美女网站在线观看视频| 精品人妻1区二区| 制服丝袜大香蕉在线| 亚洲avbb在线观看| 12—13女人毛片做爰片一| 亚洲av日韩精品久久久久久密| 国内精品久久久久久久电影| 国产精品亚洲美女久久久| avwww免费| www.www免费av| 午夜精品一区二区三区免费看| 亚洲精品亚洲一区二区| 亚洲五月天丁香| 精品一区二区三区av网在线观看| 村上凉子中文字幕在线| 性色avwww在线观看| 可以在线观看毛片的网站| 97碰自拍视频| 久久久久国内视频| 日日摸夜夜添夜夜添av毛片 | 日本一本二区三区精品| 无遮挡黄片免费观看| 免费人成视频x8x8入口观看| 啦啦啦观看免费观看视频高清| 成人国产一区最新在线观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 91av网一区二区| 欧美激情在线99| 中文在线观看免费www的网站| 我的老师免费观看完整版| or卡值多少钱| 国产aⅴ精品一区二区三区波| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| 联通29元200g的流量卡| 欧美日本亚洲视频在线播放| 99riav亚洲国产免费| 国产精品电影一区二区三区| 久久亚洲精品不卡| 国产美女午夜福利| 精品福利观看| 亚洲无线在线观看| 亚洲天堂国产精品一区在线| 国产午夜精品久久久久久一区二区三区 | 最新中文字幕久久久久| 春色校园在线视频观看| 伦精品一区二区三区| 精品一区二区免费观看| 无人区码免费观看不卡| 日本成人三级电影网站| 一个人免费在线观看电影| 乱码一卡2卡4卡精品| av在线蜜桃| 国产精品精品国产色婷婷| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色视频三级网站网址| 一区二区三区四区激情视频 | 欧美极品一区二区三区四区| 久久久久久九九精品二区国产| 偷拍熟女少妇极品色| 亚洲第一电影网av| 国产三级中文精品| 18禁裸乳无遮挡免费网站照片| 亚洲中文日韩欧美视频| 国产91精品成人一区二区三区| 美女 人体艺术 gogo| 欧美不卡视频在线免费观看| 亚洲av日韩精品久久久久久密| 99热这里只有精品一区| 亚洲av中文av极速乱 | 又紧又爽又黄一区二区| 午夜亚洲福利在线播放| 天堂av国产一区二区熟女人妻| 搡老岳熟女国产| 午夜精品在线福利| 少妇高潮的动态图| 淫秽高清视频在线观看| 久久人妻av系列| 自拍偷自拍亚洲精品老妇| 最近最新免费中文字幕在线| 国产免费男女视频| 色综合站精品国产| 国产成人a区在线观看| av中文乱码字幕在线| 精品久久久久久久末码| 特级一级黄色大片| 女生性感内裤真人,穿戴方法视频| 极品教师在线视频| 欧美+亚洲+日韩+国产| 欧美成人a在线观看| 成人精品一区二区免费| 国产精品国产三级国产av玫瑰| xxxwww97欧美| 成人毛片a级毛片在线播放| 女人十人毛片免费观看3o分钟| 亚洲欧美日韩高清在线视频| 久久这里只有精品中国| 久久草成人影院| 国产一区二区三区视频了| 免费观看人在逋| 精品一区二区三区av网在线观看| 黄色女人牲交| 三级毛片av免费| 日本成人三级电影网站| 99国产精品一区二区蜜桃av| 免费av观看视频| 99久久中文字幕三级久久日本| 在线看三级毛片| 国产亚洲精品久久久久久毛片| 黄色视频,在线免费观看| 国产色爽女视频免费观看| 少妇人妻一区二区三区视频| 亚洲人与动物交配视频| av在线蜜桃| 亚洲国产精品成人综合色| 麻豆一二三区av精品| 久久久久久伊人网av| 国产日本99.免费观看| 一区二区三区激情视频| 欧美黑人巨大hd| 国产精品人妻久久久久久| 我的女老师完整版在线观看| 亚洲男人的天堂狠狠| 国产乱人视频| 亚洲av免费高清在线观看| 人妻制服诱惑在线中文字幕| 我要看日韩黄色一级片| 真人做人爱边吃奶动态| 国产男人的电影天堂91| 国产 一区 欧美 日韩| 日日啪夜夜撸| 国产真实乱freesex| 黄色日韩在线| 91在线精品国自产拍蜜月| 99热精品在线国产| 日本精品一区二区三区蜜桃| 在线播放国产精品三级| 亚洲精品一区av在线观看| 精品欧美国产一区二区三| 国产精品一区二区三区四区免费观看 | 国产伦人伦偷精品视频| 国产伦精品一区二区三区四那| 少妇丰满av| 国产精品美女特级片免费视频播放器| 欧美三级亚洲精品| 99riav亚洲国产免费| 日本熟妇午夜| 日本免费a在线| 亚洲av中文字字幕乱码综合| av视频在线观看入口| 久久精品人妻少妇| 免费人成视频x8x8入口观看| 高清日韩中文字幕在线| 国产精华一区二区三区| 日韩国内少妇激情av| 有码 亚洲区| 黄色丝袜av网址大全| 婷婷精品国产亚洲av在线| 亚洲欧美清纯卡通| 欧美激情久久久久久爽电影| 舔av片在线| 国产 一区 欧美 日韩| 男人舔女人下体高潮全视频| 熟女人妻精品中文字幕| 日本黄色视频三级网站网址| 最后的刺客免费高清国语| 日韩精品中文字幕看吧| 两个人视频免费观看高清| 色综合婷婷激情| 在线免费观看的www视频| 欧美成人免费av一区二区三区| 日日啪夜夜撸| 免费在线观看成人毛片| 亚洲av中文av极速乱 | 不卡一级毛片| 国产成人影院久久av| 成人国产综合亚洲| 岛国在线免费视频观看| 一个人看视频在线观看www免费| 久久久久久久久久黄片| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩av片在线观看| 欧美丝袜亚洲另类 | 久久草成人影院| 日日啪夜夜撸| 很黄的视频免费| 麻豆久久精品国产亚洲av| 国产精品久久久久久av不卡| 国产白丝娇喘喷水9色精品| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 亚洲熟妇熟女久久| 亚洲三级黄色毛片| 少妇的逼好多水| 九色国产91popny在线| 午夜a级毛片| 美女高潮的动态| 神马国产精品三级电影在线观看| 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 国语自产精品视频在线第100页| 超碰av人人做人人爽久久| 色综合婷婷激情| 国产一区二区激情短视频| 成人综合一区亚洲| 天堂√8在线中文| av在线老鸭窝| 国产高清视频在线观看网站| 亚洲三级黄色毛片| 亚洲精品粉嫩美女一区| 麻豆一二三区av精品| 床上黄色一级片| 啦啦啦观看免费观看视频高清| 日韩欧美三级三区| 久久国产精品人妻蜜桃| 在线免费观看的www视频| av黄色大香蕉| 在线免费观看的www视频| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 久久草成人影院| 日韩欧美精品v在线| 夜夜夜夜夜久久久久| 99国产极品粉嫩在线观看| 在线免费观看不下载黄p国产 | 欧美日韩亚洲国产一区二区在线观看| 黄片wwwwww| 精品午夜福利视频在线观看一区| 黄片wwwwww| 91麻豆av在线| 国产不卡一卡二| 免费看日本二区| 精品人妻偷拍中文字幕| 国产成人av教育| 一区二区三区高清视频在线| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 久久久国产成人免费| 亚洲中文字幕一区二区三区有码在线看| 成人无遮挡网站| 国产乱人伦免费视频| 少妇猛男粗大的猛烈进出视频 | 国产伦在线观看视频一区| 久久久久久大精品| 日韩在线高清观看一区二区三区 | 精品久久久噜噜| 精品午夜福利在线看| 成人美女网站在线观看视频| 精品久久久久久久久av| 成人精品一区二区免费| av在线天堂中文字幕| 18+在线观看网站| 亚洲最大成人av| 中文字幕高清在线视频| 欧美性猛交╳xxx乱大交人| 午夜日韩欧美国产| 露出奶头的视频| netflix在线观看网站| 黄色视频,在线免费观看| 色5月婷婷丁香| 国产精品精品国产色婷婷| 亚洲无线观看免费| 搡老妇女老女人老熟妇| 白带黄色成豆腐渣| 亚洲va在线va天堂va国产| 在线观看免费视频日本深夜| 久久久色成人| 日韩人妻高清精品专区| 免费在线观看成人毛片| 国产午夜精品久久久久久一区二区三区 | 国产精华一区二区三区| 两个人的视频大全免费| 好男人在线观看高清免费视频| 色哟哟哟哟哟哟| 欧美色视频一区免费| 1000部很黄的大片| 一级黄片播放器| 亚洲精品在线观看二区| 极品教师在线视频| 丰满的人妻完整版| 亚洲在线观看片| 国产国拍精品亚洲av在线观看| 国产精品电影一区二区三区| 舔av片在线| 欧美激情国产日韩精品一区| 欧美zozozo另类| 国产精品1区2区在线观看.| 久久久久国内视频| 两个人的视频大全免费| 日本撒尿小便嘘嘘汇集6| 人人妻人人澡欧美一区二区| 91麻豆av在线| 久久久久久久久大av| 久久久久久久亚洲中文字幕| 琪琪午夜伦伦电影理论片6080| 色综合婷婷激情| 在线观看美女被高潮喷水网站| 天天一区二区日本电影三级| 成人无遮挡网站| 欧美日韩黄片免| 99在线视频只有这里精品首页| 大型黄色视频在线免费观看| 国产精品伦人一区二区| 热99在线观看视频| 久久久精品大字幕| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 一夜夜www| 免费无遮挡裸体视频| 久久久久国内视频| 日本黄大片高清| 国产精品一区www在线观看 | av黄色大香蕉| 国产一区二区在线av高清观看| 亚洲精华国产精华液的使用体验 | 欧美性猛交黑人性爽| 91在线观看av| 男女视频在线观看网站免费| 欧美三级亚洲精品| 99热6这里只有精品| 美女大奶头视频| 中文字幕免费在线视频6| 亚洲av.av天堂| 欧美最黄视频在线播放免费| 国产精品自产拍在线观看55亚洲| 成年版毛片免费区| 真人一进一出gif抽搐免费| 在线观看av片永久免费下载| 男人狂女人下面高潮的视频| 国产免费av片在线观看野外av| 久久久精品大字幕| 国产精品精品国产色婷婷| 又黄又爽又免费观看的视频| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 尾随美女入室| 亚洲av熟女| 又紧又爽又黄一区二区| 毛片女人毛片| 国产av在哪里看| 免费看日本二区| 久久久久久伊人网av| 日日夜夜操网爽| 国产精品女同一区二区软件 | 国产欧美日韩一区二区精品| 午夜影院日韩av| 亚洲美女搞黄在线观看 | 好男人在线观看高清免费视频| 日韩欧美国产一区二区入口| 午夜福利成人在线免费观看| 香蕉av资源在线| 又爽又黄a免费视频| 免费不卡的大黄色大毛片视频在线观看 | 在线播放无遮挡| 久久久久久久精品吃奶| 亚洲中文字幕一区二区三区有码在线看| 在线免费观看不下载黄p国产 | 中文字幕久久专区| 我要搜黄色片| 中文亚洲av片在线观看爽| 国产真实乱freesex| 免费看a级黄色片| 中文字幕免费在线视频6| 亚洲精品色激情综合| netflix在线观看网站| 一个人观看的视频www高清免费观看| 日韩中文字幕欧美一区二区| 在线看三级毛片| 国产高清激情床上av| 亚洲在线自拍视频| 日韩欧美免费精品| 一级黄色大片毛片| 午夜激情福利司机影院| 男女做爰动态图高潮gif福利片| 99热精品在线国产| а√天堂www在线а√下载| 日韩中文字幕欧美一区二区| 99热这里只有是精品在线观看| 九九久久精品国产亚洲av麻豆| 国产精品电影一区二区三区| 一区二区三区四区激情视频 | 少妇高潮的动态图| 麻豆成人午夜福利视频| www日本黄色视频网| 久久人妻av系列| 我要看日韩黄色一级片| 日本成人三级电影网站| x7x7x7水蜜桃| 久久这里只有精品中国| 麻豆一二三区av精品| 日韩av在线大香蕉| 国内少妇人妻偷人精品xxx网站| 99久国产av精品| 亚洲自偷自拍三级| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利视频1000在线观看| 亚洲美女黄片视频| 最近最新中文字幕大全电影3| 舔av片在线| 精品不卡国产一区二区三区| 久久精品国产亚洲网站| 国产精品久久久久久精品电影| 日日夜夜操网爽| av视频在线观看入口| 我的老师免费观看完整版| 老司机午夜福利在线观看视频| 亚洲综合色惰| 18+在线观看网站| 欧美成人一区二区免费高清观看| 97碰自拍视频| 黄色丝袜av网址大全| 狂野欧美激情性xxxx在线观看| 亚洲色图av天堂| 国产美女午夜福利| 九九久久精品国产亚洲av麻豆| 波多野结衣巨乳人妻| 性插视频无遮挡在线免费观看| 男女啪啪激烈高潮av片| 久久精品国产鲁丝片午夜精品 | 婷婷亚洲欧美| 熟女人妻精品中文字幕| 婷婷六月久久综合丁香| 日韩欧美在线乱码| 俺也久久电影网| 村上凉子中文字幕在线| 亚洲在线自拍视频| 日日干狠狠操夜夜爽| 色哟哟哟哟哟哟| 亚洲av不卡在线观看| 国产一区二区在线av高清观看| 看十八女毛片水多多多| 一级黄片播放器| 免费电影在线观看免费观看| 久久久久免费精品人妻一区二区| 18禁裸乳无遮挡免费网站照片| 国产精品福利在线免费观看| 久久人妻av系列| 国产伦人伦偷精品视频| 亚洲成人久久性| 亚洲男人的天堂狠狠| 亚洲精品国产成人久久av| 精品不卡国产一区二区三区| 干丝袜人妻中文字幕| 国产精品一区二区三区四区久久| 国产毛片a区久久久久| 国产成人av教育| 久久久久久久久大av| 欧美人与善性xxx| 欧美色视频一区免费| 午夜a级毛片| 日韩欧美在线乱码| 美女cb高潮喷水在线观看| 美女大奶头视频| 国产av麻豆久久久久久久| 久久午夜福利片| 村上凉子中文字幕在线| 国内精品久久久久精免费| 两个人视频免费观看高清| 男女啪啪激烈高潮av片| 国产aⅴ精品一区二区三区波| 校园人妻丝袜中文字幕| 麻豆久久精品国产亚洲av| 内地一区二区视频在线| 国产精品精品国产色婷婷| 亚洲中文字幕一区二区三区有码在线看| 老司机福利观看| a级毛片免费高清观看在线播放| 国产亚洲精品久久久久久毛片| 免费看日本二区| 久久久久久国产a免费观看| 女同久久另类99精品国产91| 亚洲精品456在线播放app | 又黄又爽又免费观看的视频| 国产麻豆成人av免费视频| 亚洲综合色惰| 天天一区二区日本电影三级| 国产成人av教育| 久久久久久久久久黄片| 国产精品国产高清国产av| 久久精品国产亚洲网站| 国产男人的电影天堂91| 亚洲精品一卡2卡三卡4卡5卡| 欧美色欧美亚洲另类二区| 欧美精品国产亚洲| 亚洲av二区三区四区| xxxwww97欧美| 变态另类丝袜制服| 女同久久另类99精品国产91| 如何舔出高潮| 亚洲最大成人av| 麻豆久久精品国产亚洲av| 久久精品国产鲁丝片午夜精品 | 久久精品91蜜桃|