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

    紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)

    2017-10-13 12:32:53董伊晨劉艷紅
    生態(tài)學(xué)報 2017年17期
    關(guān)鍵詞:紅松林苗齡紅松

    董伊晨, 劉艷紅

    北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083

    紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)

    董伊晨, 劉艷紅*

    北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083

    通過分析不同苗齡紅松幼苗在控制溫度和光照條件下葉性狀的變化差異,以及各控制因素交互作用對紅松幼苗葉性狀的影響程度和紅松葉性狀的線性相關(guān)關(guān)系,探討紅松幼苗葉性狀響應(yīng)氣候環(huán)境因素變化的特征和所采取的物質(zhì)分配上的變化策略。實(shí)驗(yàn)采取人工氣候箱培養(yǎng)紅松幼苗,設(shè)置3組溫度為10、20、30℃,每組溫度下設(shè)置有效光照輻射4組:45.3、90.5、181.1、271.6(μmol m-2s-1),每組觀測苗齡1—4a。觀測性狀為:比葉重(LMA)、單位面積光合能力(Aarea)、單位質(zhì)量葉氮含量(Nmass)和單位質(zhì)量葉磷含量(Pmass)。結(jié)果顯示:(1)紅松幼苗的LMA、Aarea、Nmass、Pmass受到溫度、光照和苗齡的影響,各控制組葉性狀間差異顯著(P<0.05);(2)溫度、光照和苗齡對紅松幼苗葉性狀的交互作用效應(yīng)顯著(P<0.05),影響程度依次為:苗齡>溫度>溫度×光照×苗齡>光照>光照×苗齡>溫度×光照>溫度×苗齡;(3)LMA與溫度和苗齡呈正相關(guān)性,Aarea與苗齡呈正相關(guān)性,與溫度呈負(fù)相關(guān)性,Nmass與溫度、光照、苗齡、LMA和Aarea均呈負(fù)相關(guān)性,Pmass與溫度、Nmass呈正相關(guān)性而與苗齡和LMA呈負(fù)相關(guān)性。通過結(jié)果分析可以得出以下結(jié)論:(1)紅松幼苗葉性狀對溫度、光照和苗齡3個因素均有明顯響應(yīng)但響應(yīng)程度不一,影響性狀差異的主要因素是溫度和苗齡,光照在其溫度和苗齡二者的作用基礎(chǔ)上對性狀產(chǎn)生影響;(2)紅松幼苗葉性狀之間的線性關(guān)系實(shí)際反映的是對溫度、光照和苗齡響應(yīng)的結(jié)果,葉性狀間的密切聯(lián)系反映出葉片在響應(yīng)環(huán)境因素時所采取的分配策略的調(diào)整,其目的是維持植物體內(nèi)同化作用和異化作用的平衡,確保紅松幼苗穩(wěn)定生長。

    紅松幼苗; 控制實(shí)驗(yàn); 比葉重; 光合能力; 葉氮磷含量; 因素交互作用

    Abstract: In the present study, changes in leaf traits ofPinuskoraiensis(Korean pine) seedlings in response to seedling age, temperature, and light conditions were investigated, as well as the correlation of these factors with Korean pine seedling leaf traits. Korean pine seedlings were cultivated in artificial climate chambers with three levels of temperature settings (10, 20, and 30 ℃) and four levels of effective light radiation applied to each temperature level (45.3, 90.5, 181.1, and 271.6 μmol m-2s-1); the observed seedlings in each group were aged between 1 and 4 years. The observed leaf mass traits were per area (LMA), area-based photosynthetic capacity (Aarea), mass-based nitrogen concentration (Nmass), and mass-based phosphorus concentration (Pmass). Our results showed that (1) LMA,Aarea,Nmass, andPmassof Korean pine seedlings were affected by temperature, light, and seedling age. The differences in leaf traits among different controlled groups were statistically significant (P<0.05). (2) Temperature, light, and seedling age in combination affected Korean pine seedling leaf traits. The order of significance of effects was seedling age > temperature > temperature × light × seedling age > light > light × seedling age > temperature × light > temperature × seedling age. (3) LMA was positively correlated with temperature and seedling age;Aareawas positively correlated with seedling age but negatively correlated with temperature;Nmasswas negatively correlated with temperature, light, seedling age, LMA, andAarea; andPmasswas positively correlated with temperature andNmassbut was negatively correlated with seedling age and LMA. Our results suggest the following: (1) Korean pine seedling leaf traits responded to all three studied factors (temperature, light, and seedling age) to differing degrees. The main factors that caused differences among traits were temperature and seedling age, while the effect of light was influenced by these two factors. (2) The correlation between Korean pine seedling leaf traits reflected their response to temperature, light, and seedling age. The close relationship between leaf traits is consistent with the allocation strategies that the leaves use in response to environmental factors, which maintain the balance between assimilation and respiration to ensure stable growth in Korean pine seedlings.

    KeyWords: Korean pine seedlings; controlled experiments; leaf mass per area; photosynthetic capacity; leaf nitrogen and phosphorus concentrations; factor interaction

    紅松針闊混交林是中國東北東部山區(qū)的典型森林群落類型,紅松作為主要優(yōu)勢樹種在該地區(qū)森林演替更新、群落結(jié)構(gòu)動態(tài)和生物多樣性組成等方面起到關(guān)鍵作用,因此其生態(tài)特性一直是人們關(guān)注和研究的對象[1-4]。此外,紅松林群落及其物種多樣性對氣候動態(tài)響應(yīng)關(guān)系等方面的研究又使其與全球氣候變化、多樣性保護(hù)等熱點(diǎn)問題[5-6]緊密聯(lián)系在一起。近年來,在對前輩研究不斷深入和補(bǔ)充的基礎(chǔ)上,利用植物功能性狀(plant functional traits)和葉功能性狀(leaf functional traits)評價紅松群落生態(tài)學(xué)特征和氣候適應(yīng)性的方法,對紅松群落數(shù)量化理論和模型研究起到推動作用[7-8]。通過對紅松葉性狀的長期監(jiān)測不僅可以掌握紅松林群落動態(tài)變化,而且可以對紅松生長的環(huán)境影響程度和物質(zhì)分配策略進(jìn)行評價[9- 10],為紅松響應(yīng)環(huán)境變化方面的研究提供理論依據(jù)和方法。

    植物功能性狀和葉功能性狀研究的主要目的是闡述植物在進(jìn)化過程中所形成的環(huán)境適應(yīng)機(jī)制[11]。植物性狀在環(huán)境干擾評價、量化解釋生態(tài)系統(tǒng)功能和多樣性以及易于調(diào)查獲取等方面的優(yōu)勢廣泛應(yīng)用于生態(tài)學(xué)不同尺度的研究[9,12-17]。其中葉經(jīng)濟(jì)型譜(leaf economics spectrum)已形成較完整的理論體系[18- 19],在其基礎(chǔ)上建立的全球植物性狀數(shù)據(jù)庫[20- 21]也廣受關(guān)注。此外,植物葉性狀與植物代謝和植物相關(guān)生長等生態(tài)學(xué)理論聯(lián)系緊密,通過模型擬合的方法可以進(jìn)一步深化植物生長與氣候環(huán)境適應(yīng)關(guān)系的研究[22-28]。

    目前氣候環(huán)境因素影響植物葉性狀的研究多以全球、區(qū)域或群落尺度為主,而對特定種植物葉性狀的研究卻較鮮見[21],盡管如此,個體物種研究依然可以通過其性狀的生理特征闡述植物的環(huán)境適應(yīng)性和生態(tài)系統(tǒng)過程和功能[11,17]。因此研究紅松葉性狀無論是對紅松樹種還是紅松林群落意義都很明顯:首先通過葉性狀描述紅松對氣候環(huán)境的適應(yīng)性過程符合目前植物葉性狀的主要研究方向;其次通過葉性狀間的線性關(guān)系論述紅松內(nèi)部組織的物質(zhì)分配策略符合葉經(jīng)濟(jì)型譜的主要研究內(nèi)容[18];紅松及其群落作為東北亞地區(qū)典型森林群落類型,是世界葉經(jīng)濟(jì)型譜中重要的組成部分[20,29]。

    氣候環(huán)境因素對植物生長的影響是綜合性的[30-32],溫度和光照是其中重要的影響因素。溫度限制植物個體生長、生物產(chǎn)量和物種分布范圍[33-36],植物葉性狀對溫度變化的響應(yīng)機(jī)制反映了其生長適應(yīng)性,其形態(tài)、物候及生長型都是長期適應(yīng)進(jìn)化的結(jié)果[27,37]。光照強(qiáng)度影響植物光合作用中光響應(yīng)過程,進(jìn)而影響植物的生長和生物量積累[38-40]。溫度和光照對紅松生長分布及更新的意義明顯,紅松葉性狀對溫度和光照變化的響應(yīng)特征在紅松林群落結(jié)構(gòu)和物種多樣性研究中作用突出,尤其是光照對天然紅松林更新的影響,作為紅松林群落演替過程的重點(diǎn)問題一直備受關(guān)注[2,41-47]。

    基于上述原因,本文旨在通過分析紅松幼苗葉性狀在氣候環(huán)境因素影響下的響應(yīng)特征來判斷:(1)能否通過葉性狀的變化差異判斷紅松受到環(huán)境因素影響;(2)能否通過葉性狀的響應(yīng)差異判斷環(huán)境因素間影響程度是否相同;(3)能否通過葉性狀相關(guān)性判斷環(huán)境因素的影響是否改變紅松物質(zhì)分配的策略。此外,我國東北地區(qū)的紅松闊葉林所面臨現(xiàn)實(shí)困境依然嚴(yán)峻:如紅松林分布減少、紅松林?jǐn)?shù)量萎縮等風(fēng)險[3]。因此進(jìn)行紅松個體性狀研究除可以為紅松群落功能性狀提供數(shù)據(jù)支持,也為紅松的保護(hù)和栽培提供實(shí)驗(yàn)依據(jù)。

    1 實(shí)驗(yàn)材料與方法

    1.1 實(shí)驗(yàn)材料

    紅松幼苗產(chǎn)自黑龍江小興安嶺南坡。先將紅松幼苗在北京林業(yè)大學(xué)氣象站實(shí)驗(yàn)地盆栽培養(yǎng),其栽培用土為產(chǎn)苗當(dāng)?shù)赝寥馈_x擇苗齡在1—4a的健壯紅松幼苗作為實(shí)驗(yàn)樣本。

    1.2 實(shí)驗(yàn)控制條件

    實(shí)驗(yàn)時間為2014年3—9月間,將紅松幼苗放置在人工氣候箱(RXZ- 280C型)中進(jìn)行溫度和光照處理。溫度、光照、樹齡三者對紅松幼苗是否存在交互作用是本實(shí)驗(yàn)關(guān)注的問題之一,為有效判斷這些因素的功效作用,實(shí)驗(yàn)采取的控制條件設(shè)置為:(1)根據(jù)東北地區(qū)紅松針闊混交林在生長季時的溫度范圍[3,33]設(shè)置10、20、30℃,3個溫度處理組,視10—20℃為紅松生長季溫度,將30℃視為極端溫度的臨界值,低于10℃可能因生理低溫的影響不作為考慮溫度[48-49]。(2)根據(jù)天然紅松林下更新的遮蔭影響[45- 47],將光照范圍設(shè)定在實(shí)際光有效輻射40—300μmol m-2s-1之間,以此測定紅松幼苗在微光環(huán)境下的葉性狀變化。根據(jù)每個溫度組分別設(shè)置4個光照控制組:2200、4400、8800、13200lx(氣候箱實(shí)際光照度最大值22000lx),文中將光強(qiáng)換算為光有效輻射(換算公式:1μE-1·s-1=48.6lx)[50],后文均以光輻射值45.3、90.5、181.1、271.6 μmol m-2s-1表示光照條件。(3)按照實(shí)驗(yàn)條件的設(shè)置,每個溫度下均有4個不同光照的氣候箱,每個氣候箱中放置1—4a苗齡的紅松幼苗,每苗齡5株重復(fù),全部實(shí)驗(yàn)組共計(jì)48組(即48水平處理組)。每日光照時長為14h,箱內(nèi)空氣濕度(RH%)保持在(60±5)%。

    1.3 葉性狀測量方法

    在紅松幼苗生長旺盛的8—9月進(jìn)行性狀測定??紤]到針葉壽命長和往年生葉片的差異,故選擇生長情況較一致的當(dāng)年生葉片進(jìn)行性狀測量,以排除選擇干擾。葉性狀測定方法:(1)比葉重LMA(g/m2),取針葉3簇,視紅松葉為5等分圓柱體,視等分體外側(cè)弧狀面為光輻射面,葉面積即為光輻射面積總和。利用尺量葉長,掃描葉橫截面(Image-Pro Plus 6.0計(jì)算弧長),葉面積=弧長×葉長×5(未計(jì)楔狀面);烘干葉片稱取葉干重, LMA=葉干重/葉面積,LMA計(jì)算曾通過經(jīng)驗(yàn)?zāi)P陀?jì)算法比對[16,51]。(2)單位面積葉光合能力Aarea(μmolC m-2s-1),利用LI6400便攜式光合測定儀(LiCor, Lincoln, NE,USA)和針狀葉室(needle chamber, IRGA)測定紅松幼苗凈光合速率,每株取葉3簇。測定環(huán)境與處理環(huán)境相同,測定前記錄測定葉寬,測定時手動輸入:測量面積=葉寬×葉室長(若葉不及葉室長則輸入葉面積),選擇最大瞬時光合速率(10:00—11:00)作為單位面積光合能力Aarea[16,52]。(3)單位質(zhì)量葉氮含量Nmass(%)利用凱式定氮法測定。(4)單位質(zhì)量葉磷含量Pmass(%)利用鉬銻抗火焰光度計(jì)法測定。

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

    紅松幼苗葉性狀數(shù)據(jù)利用SPSS 18.0進(jìn)行多因素方差分析(GLM Multivariate),得到各組均值(mean±SE)、Duncan一致性子集檢驗(yàn)各組差異顯著值以及多因素功效檢驗(yàn)解釋值。利用相關(guān)性分析(Bivariate, pearson correlation)、一元線性回歸(Linear Regression, least-squares)分析性狀與因素關(guān)系、性狀間關(guān)系等,其中葉性狀指標(biāo)均進(jìn)行對數(shù)轉(zhuǎn)化(log10-transformed)[18-19],以降低方差非齊性可能造成的數(shù)據(jù)離散和判定概率系數(shù)P值偏小的干擾。利用Origin 8.0和Excel 2007制作圖表。

    2 結(jié)果分析

    2.1 不同苗齡紅松幼苗葉性狀對溫度、光照條件變化的響應(yīng)

    紅松幼苗葉性狀指標(biāo):LMA、Aarea、Nmass和Pmass對溫度、光照、苗齡3因素48水平(溫度3水平、光照4水平、苗齡4水平)的響應(yīng)呈現(xiàn)不同的變化趨勢(表1—表4)。實(shí)驗(yàn)結(jié)果表明溫度、光照、苗齡的變化對紅松幼苗葉性狀均產(chǎn)生明顯影響。

    表1 紅松幼苗比葉重(LMA)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    表2 紅松幼苗單位面積光合能力(Aarea)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    (1)紅松幼苗LMA的變化(表1):在不同溫度或光照組條件下,紅松幼苗LMA隨著苗齡的增長而上升,各苗齡組間差異顯著(P<0.05)。在固定苗齡組比較不同光照條件中LMA均值隨溫度的升高而呈現(xiàn)上升的趨勢,各溫度組間差異顯著(P<0.05)。相同溫度和苗齡組,紅松幼苗LMA受光輻射強(qiáng)度影響的組間差異顯著(P<0.05)。

    (2)紅松幼苗Aarea變化趨勢(表2):在不同的光照和苗齡組,溫度由10℃上升到20℃時Aarea呈現(xiàn)隨之升高的變化,而溫度持續(xù)升高到30℃時Aarea出現(xiàn)明顯下降,溫度組間差異顯著(P<0.05)。在相同的溫度和光照組中,隨苗齡的增長Aarea呈現(xiàn)上升的變化趨勢,各苗齡組差異顯著(P<0.05)。而Aarea隨光照有效輻射的升高,同一苗齡組有上升的趨勢,各光照組間差異顯著(P<0.05)。

    (3)紅松幼苗Nmass的變化趨勢(表3):在不同的溫度和光照組,Nmass隨苗齡的增長而下降,各苗齡間差異顯著(P<0.05)。在相同的溫度和苗齡組Nmass隨光有效輻射的增高而下降,各光照組間差異顯著(P<0.05)。Nmass在同一苗齡組隨溫度上升而減少,各溫度組間差異顯著(P<0.05)。

    (4)紅松幼苗Pmass的變化趨勢(表4):在不同的溫度和光照組中,Pmass隨苗齡的增長呈現(xiàn)小幅降低,各苗齡組間差異顯著(P<0.05)。在相同的光照和苗齡組中,Pmass隨溫度的升高而出現(xiàn)波動變化,但仍呈現(xiàn)下降的趨勢,各溫度組間差異顯著(P<0.05)。在同一的溫度組中Pmass隨光有效輻射的增強(qiáng)出現(xiàn)不同的變化,但在同一的苗齡組中隨光照變化的趨勢不明顯,各光照組差異顯著(P<0.05)。

    表3 紅松幼苗單位質(zhì)量葉氮含量(Nmass)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    表4 紅松幼苗單位質(zhì)量葉磷含量(Pmass)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    2.2 溫度、光照和苗齡對紅松幼苗葉性狀影響的交互作用功效分析

    紅松幼苗葉性狀對溫度、光照和苗齡3個因素的響應(yīng)變化趨勢不同,通過交互作用分析和功效檢驗(yàn)對三者的影響程度做進(jìn)一步分析。表5給出的是溫度、光照和苗齡3個因素對紅松幼苗葉性狀的交互作用和功效分析,3個單項(xiàng)(溫度、光照、苗齡)表示因素交互作用的主效應(yīng),3個二維項(xiàng)(溫度×光照、溫度×苗齡、光照×苗齡)表示的是兩兩因素的一階交互效應(yīng),1個三維項(xiàng)(溫度×光照×苗齡)則表示的是3因素的二階交互效應(yīng)。Wilks′ Lambda值為基于組間無差別的零級假設(shè)判別系數(shù),系數(shù)越接近于0則說明因素功效越有差異,檢驗(yàn)概率P<0.05表明各因素的影響差異顯著。

    主效應(yīng)判別中溫度和苗齡的系數(shù)最小(0.003、0.035),低于光照的系數(shù)(0.163),說明溫度和苗齡的獨(dú)立作用效率高于光照。二維交互項(xiàng)判別中溫度×光照系數(shù)最小(0.115),其次是光照×苗齡(0.136),溫度×苗齡系數(shù)最大(0.223),說明溫度×光照和光照×苗齡的交互作用效率高于溫度×苗齡,也即是說溫度和苗齡二者獨(dú)立作用效率高于二者的交互作用,光照×溫度、光照×苗齡的交互作用效率高于光照的獨(dú)立作用。三維交互項(xiàng)溫度×光照×苗齡系數(shù)(0.087)高于其他所有二維項(xiàng),說明溫度×光照×苗齡三者交互作用效率高于任意兩兩交互作用。通過交互作用判別可以看出光照的獨(dú)立作用程度受到溫度和苗齡的制約,同時光照也使溫度和苗齡的交互作用效率有所降低。

    效應(yīng)歸因解釋(Partial Eta-Squared)給出的因素解釋概率η2值指出溫度、光照和苗齡3因素對紅松幼苗葉性狀影響程度的強(qiáng)弱關(guān)系。各因素交互作用對葉性狀的解釋率分別為:苗齡對LMA的主效應(yīng)解釋概率最高(0.861)。溫度對Aarea的主效應(yīng)解釋率最高(0.990),其次為苗齡(0.885)、溫度×光照(0.703)、光照(0.629),說明紅松幼苗光合能力對環(huán)境因素的響應(yīng)過程比較復(fù)雜。光照和苗齡對Nmass的主效應(yīng)解釋率較高(0.630、0.794),其次是光照×苗齡(0.604),說明其二者的交互效應(yīng)對紅松幼苗葉性狀的影響既有獨(dú)立性又有相互促進(jìn)性。苗齡對Pmass的主效應(yīng)解釋率最高(0.792),其次為溫度(0.645)和溫度×光照×苗齡(0.613)??偨Y(jié)表5功效解釋率依次為:苗齡>溫度>溫度×光照×苗齡>光照>光照×苗齡>溫度×光照>溫度×苗齡。

    表5 各影響因素功效檢驗(yàn)和對葉性狀的交互作用解釋

    Wilks′ Lambda值為各因素差異性檢驗(yàn),F檢驗(yàn)概率P<0.05說明差異顯著,效應(yīng)估計(jì)η2值表示各因素功效解釋度;*因素為溫度,光照,苗齡,溫度×光照,溫度×苗齡,光照×苗齡,溫度×光照×苗齡

    通過分析比較紅松幼苗各葉性狀與溫度、光照和苗齡的交互作用和功效解釋率,發(fā)現(xiàn)各因素間的主效應(yīng)、二維交互效應(yīng)和三維交互效應(yīng)都較顯著(P<0.05),且對葉性狀的影響過程較為復(fù)雜。說明溫度、光照和苗齡3因素對紅松幼苗葉性狀的影響程度不同,因素間的獨(dú)立性和依賴性在交互過程中由于分組水平的變化而產(chǎn)生差異,因此在描述葉性狀的差異變化上顯得龐雜,然而這也是客觀現(xiàn)實(shí)的體現(xiàn)。利用因素的解釋度來分析溫度、光照和苗齡對紅松幼苗葉性狀影響程度,可以使結(jié)果更加清晰。

    2.3 紅松幼苗葉性狀間及與溫度、光照及苗齡等因素的相關(guān)關(guān)系

    在肯定溫度、光照和苗齡對紅松幼苗葉性狀影響程度存在差別的基礎(chǔ)上,通過相關(guān)性分析對各葉性狀與因素之間的密切程度,以及葉性狀間的相互關(guān)聯(lián)進(jìn)行分析。表6為紅松幼苗葉性狀與溫度、光照、苗齡的相關(guān)關(guān)系,以及葉性狀之間的相關(guān)關(guān)系。采用Pearson相關(guān)性檢驗(yàn)和雙尾顯著性t-檢驗(yàn),概率分為P<0.01、P<0.05表示顯著性。其中LMA與溫度和苗齡存在顯著正相關(guān)關(guān)系;Aarea與苗齡存在顯著正相關(guān)關(guān)系,與溫度成顯著負(fù)相關(guān)關(guān)系;Nmass與溫度、光照、苗齡、LMA和Aarea均成顯著負(fù)相關(guān)關(guān)系;Pmass與溫度、Nmass成顯著正相關(guān)關(guān)系而與苗齡和LMA則呈顯著負(fù)相關(guān)關(guān)系。

    表6 紅松幼苗葉性狀與溫度、光照間的相關(guān)性

    各葉性狀對數(shù)值采用Pearson Correlation 相關(guān)性檢驗(yàn)(2- tailed),**表示P<0.01顯著,*表示P<0.05顯著

    3 討論

    3.1 紅松幼苗葉性狀與溫度、光照和苗齡的響應(yīng)關(guān)系

    通過實(shí)驗(yàn)結(jié)果分析可知,紅松幼苗葉性狀對溫度、光照和苗齡響應(yīng)差異顯著,溫度、光照和苗齡既有獨(dú)立作用也有交互作用,且影響過程復(fù)雜,光照的作用是在溫度和苗齡的作用基礎(chǔ)上體現(xiàn)的。紅松幼苗葉性狀與溫度、光照和苗齡相關(guān)性表現(xiàn)為:溫度和苗齡對4種葉性狀都具有顯著相關(guān)性,而光照除與Nmass相關(guān)外與其他性狀沒有顯著相關(guān)性,說明紅松葉性狀更多表現(xiàn)為響應(yīng)溫度和苗齡的變化,響應(yīng)關(guān)系明顯。這些關(guān)系反映出:一方面苗齡對Aarea變化功效顯著,兩者顯著正相關(guān),即苗齡越高紅松的生物積累越多,另一方面溫度與LMA、Pmass呈正顯著相關(guān)性,說明紅松幼苗因?yàn)闇囟鹊纳仙訌?qiáng)了呼吸和代謝的作用,LMA上升是葉片呼吸造成水分散失加劇,導(dǎo)致葉面積收縮相對減少的結(jié)果。而LMA和苗齡的顯著正相關(guān)性又說明紅松幼苗生長和抵御環(huán)境干擾的能力也隨苗齡增強(qiáng)。磷素在植物細(xì)胞核質(zhì)和線粒體為補(bǔ)充轉(zhuǎn)錄RNA提供構(gòu)件[53-54],在群落研究中屬于生境限制因素[19,29],部分研究也利用磷素在膜質(zhì)結(jié)構(gòu)的多少來判斷植物速生與否[21,26,55]。對于紅松樹種而言,其本身非速生樹種,Pmass更多反映的是線粒體的代謝水平,它與LMA都反映紅松幼苗對溫度升高的響應(yīng)措施:隨著溫度逐步升高代謝增強(qiáng)導(dǎo)致Pmass下降,水分物質(zhì)流失加劇,對紅松生長帶來不利影響。而Pmass與苗齡呈負(fù)顯著相關(guān)性,即在溫度持續(xù)升高條件下紅松幼苗降低磷素代謝的水平隨年齡遞增而降低,這樣的措施抑制了部分水分物質(zhì)散失量,光合作用所積累的物質(zhì)消耗相對減少,可以避免紅松生長時因環(huán)境干擾所帶來不利影響。

    紅松幼苗的生長需光量的研究國內(nèi)已有不少報道[2,45- 47,56],涉及到紅松林下更新和幼苗生長適應(yīng)性等問題,對紅松幼苗需光或需庇蔭的討論較多[46-47,56]。從本實(shí)驗(yàn)的數(shù)據(jù)結(jié)果上看:首先功效檢驗(yàn)溫度×光照、光照×苗齡、溫度×光照×苗齡交互作用效率要比光照單獨(dú)作用顯著,其次紅松幼苗葉性狀中的除Nmass與光照顯著負(fù)相關(guān)外,LMA、Aarea、Pmass與光照均無顯著相關(guān)性,說明光照對紅松幼苗的生長獨(dú)立影響不明顯,而是與溫度和苗齡共同影響。此外在本實(shí)驗(yàn)中對紅松幼苗Aarea的功效解釋較高的因素有溫度、光照、苗齡、溫度×光照,也可以看出紅松光合作用受到多種因素限制,在已知的研究中也體現(xiàn)出這一特點(diǎn)[57-58]。不僅如此,本實(shí)驗(yàn)中Aarea在有效光輻射45.3,90.5,181.1,271.6(μmol m-2s-1)4種條件下數(shù)值差異顯著且沒有低于補(bǔ)償值范圍,同化依然高于呼吸,說明即使在這種光照偏低的條件下紅松幼苗依然可以進(jìn)行光合作用。在前人的研究中通過紅松幼苗的光響應(yīng)曲線判斷其在林下光適應(yīng)的大致范圍從20到900不等,全光可接近2000 μmol m-2s-1,說明紅松幼苗的光適應(yīng)范圍較寬[56,58-59]。

    綜上所述,溫度和光照對不同苗齡的紅松幼苗葉性狀影響差異顯著,葉性狀表現(xiàn)是溫度、光照和苗齡的獨(dú)立作用和交互作用共同影響的結(jié)果。溫度和苗齡影響功效更顯著,光照在其溫度和苗齡二者的影響基礎(chǔ)上產(chǎn)生作用。

    3.2 紅松幼苗葉性狀間線性相關(guān)關(guān)系

    通過相關(guān)性分析得知紅松幼苗葉性狀間存在顯著相關(guān)性,說明在響應(yīng)環(huán)境因素影響時葉性狀之間存在密切聯(lián)系。表7為紅松幼苗葉性狀logLMA、logAarea、logNmass和logPmass的一元線性擬合結(jié)果。其中l(wèi)ogLMA與logNmass、logPmass為負(fù)線型相關(guān),解釋概率分別為35%和30%,即紅松幼苗LMA與Nmass和Pmass變化趨勢相反。紅松幼苗比葉重在溫度、光照影響下隨苗齡上升而上升,而葉片中氮磷元素含量卻隨之減少,這說明紅松幼苗在維持光合呼吸過程中消耗的元素增多,而對氮磷元素吸收的能力因?yàn)榄h(huán)境因素的影響而降低[45]。logNmass和logPmass的正線性相關(guān)說明二者的變化趨勢是類似的(解釋概率12%),原因可能是在響應(yīng)外界干擾時,紅松幼苗需要同時調(diào)節(jié)光合和呼吸水平[58- 59]。除此以外,logNmass與logAarea之間呈顯著負(fù)線型關(guān)系,氮素含量的下降并沒有使紅松幼苗光合能力降低,說明說明光合作用在氮含量降低的情況下仍能維持。同時logNmass與logAarea之間的關(guān)系也并非簡單的數(shù)量上相關(guān)(解釋概率2%),還與氮參與光合作用的形式、機(jī)制以及利用效率有關(guān)[60],氮磷元素含量在光合作用和葉形態(tài)上相互關(guān)聯(lián)[61- 62],這其中氮素提供光合作用所需原件,磷素分別為膜系統(tǒng)和分解代謝提供所需原件[63-64],它們含量的降低主要是溫度上升導(dǎo)致LMA上升,水分和元素散失加劇的影響。紅松幼苗比葉重和光合能力在溫度和光照條件影響下對氮磷含量降低的表現(xiàn)反映的是其維持光合、呼吸之間的平衡的結(jié)果。

    表7 紅松幼苗葉性狀線性回歸分析

    根據(jù)性狀相關(guān)關(guān)系的顯著性進(jìn)行線性回歸分析,表中各項(xiàng)數(shù)值基于最小二乘法(least-squares regression)的結(jié)果

    紅松維持光合和呼吸平衡的過程中LMA和Aarea間雖無顯著線性相關(guān)性,但二者均與溫度、苗齡顯著相關(guān),存在間接關(guān)聯(lián)。LMA和Aarea都反映植物干物質(zhì)積累與葉生長關(guān)系[18-19,55],光合作用又受到胞間CO2濃度、水分利用能力和氣孔等限制[60,65]。紅松幼苗LMA的變化是由溫度的升高引起的,葉片水分在高溫下散失加劇,導(dǎo)致葉表面積收縮干物質(zhì)比重上升,而葉片通過降低氣體流量和水分散失來應(yīng)對這種干擾,導(dǎo)致光合能力下降,當(dāng)生長年限增加時紅松對這種干擾的抵抗能力都有所加強(qiáng)。所以盡管LMA與Aarea之間不存在線性相關(guān)性,但對相同影響因素所做出的響應(yīng)都是降低物質(zhì)水分消耗,維持紅松幼苗正常生長的過程,具有一定的相互影響。

    綜上所述,紅松幼苗葉性狀間的線性關(guān)系是對溫度、光照和苗齡響應(yīng)的結(jié)果。葉性狀間緊密聯(lián)系反應(yīng)的是紅松為響應(yīng)環(huán)境因素影響所采取的物質(zhì)分配策略的變化。在這個過程中同一因素對紅松葉性狀的影響方式和程度都有差別,同一性狀所響應(yīng)的因素也不相同,其目的都是維持植物體內(nèi)同化和異化作用的平衡,即保障紅松幼苗穩(wěn)定生長。

    4 結(jié)論

    通過上述討論可以明確紅松幼苗葉性狀對溫度、光照和苗齡的響應(yīng)本質(zhì)是其在不同條件下維持自身生長所采取的措施。葉性狀與3種因素間的關(guān)系,以及葉性狀之間的變化關(guān)系是其在調(diào)節(jié)同化和異化作用平衡所采取的相應(yīng)策略。因此總結(jié)數(shù)據(jù)分析和討論的結(jié)果可以得出的結(jié)論為:

    (1)紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)差異顯著,葉性狀差異是溫度、光照和苗齡的獨(dú)立作用和交互作用共同影響的結(jié)果。溫度和苗齡影響功效更顯著,光照在其溫度和苗齡二者的影響基礎(chǔ)上產(chǎn)生作用。

    (2)紅松幼苗葉性狀間的線性關(guān)系是對溫度、光照和苗齡響應(yīng)的結(jié)果。葉性狀間緊密聯(lián)系反應(yīng)的是紅松為響應(yīng)環(huán)境因素影響所采取的物質(zhì)分配策略的變化,其目的是維持植物體內(nèi)同化和異化作用的平衡保障紅松幼苗穩(wěn)定生長。

    除此以外,從實(shí)驗(yàn)中不同苗齡紅松幼苗葉性狀做出的響應(yīng)過程看氣候環(huán)境因素的作用是綜合的。如何對氣候環(huán)境因素分別進(jìn)行定量劃分范圍和歸因解釋是一個難點(diǎn),而進(jìn)行人為控制實(shí)驗(yàn)正是要將不同影響因素進(jìn)行逐一驗(yàn)證,尋找準(zhǔn)確的解釋方法。本實(shí)驗(yàn)選取的溫度、光照是重要的氣候因素,對紅松生長發(fā)育意義明顯,而苗齡(或樹齡)既是性狀的限定因素,又可以指示紅松林群落生物積累,是聯(lián)系個體紅松生長和種群分布、更新的橋梁。對幾種因素間的獨(dú)立、交互作用效果分別進(jìn)比較,可以較為準(zhǔn)確地判斷引起性狀變化的原因,為今后更有針對性地開展紅松的性狀研究提供了依據(jù)。

    [1] 朱濟(jì)凡, 劉慎諤, 王戰(zhàn), 馮宗煒, 劉同生. 小興安嶺紅松針闊葉混交林. 林業(yè)科學(xué), 1958, 4(4): 355- 369.

    [2] 李俊清. 闊葉紅松林中紅松的分布格局及其動態(tài). 東北林業(yè)大學(xué)學(xué)報, 1986, 14(1): 33- 38.

    [3] 王業(yè)蘧. 闊葉紅松林. 哈爾濱: 東北林業(yè)大學(xué)出版社, 1995.

    [4] 于振良, 趙士洞, 董立榮, 張?jiān)? 闊葉紅松林的研究綜述. 吉林林學(xué)院學(xué)報, 1996, 12(4): 235- 238.

    [5] 鄧慧平, 吳正方, 周道瑋. 全球氣候變化對小興安嶺闊葉紅松林影響的動態(tài)模擬研究. 應(yīng)用生態(tài)學(xué)報, 2000, 11(1): 43- 46.

    [6] 吳曉莆, 朱彪, 趙淑清, 樸世龍, 方精云. 東北地區(qū)闊葉紅松林的群落結(jié)構(gòu)及其物種多樣性比較. 生物多樣性, 2004, 12(1): 174- 181.

    [7] 胡耀升, 么旭陽, 劉艷紅. 長白山不同演替階段森林植物功能性狀及其與地形因子間的關(guān)系. 生態(tài)學(xué)報, 2014, 34(20): 5915- 5924.

    [8] 么旭陽, 胡耀升, 劉艷紅. 長白山闊葉紅松林不同群落類型的植物功能性狀與功能多樣性. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 2014, 42(3): 77- 84.

    [9] McIntyre S, Lavorel S, Landsberg J, Forbes T D A. Disturbance response in vegetation-towards a global perspective on functional traits. Journal of Vegetation Science, 1999, 10(5): 621- 630.

    [10] Wright I J, Reich P B, Cornelissen J H C, Falster D S, Groom P K, Hikosaka K, Lee W, Lusk C H, Niinemets ü, Oleksyn J, Osada N, Poorter H, Warton D I, Westoby M. Modulation of leaf economic traits and trait relationships by climate. Global Ecology and Biogeography, 2005, 14(5): 411- 421.

    [11] 孟婷婷, 倪健, 王國宏. 植物功能性狀與環(huán)境和生態(tài)系統(tǒng)功能. 植物生態(tài)學(xué)報, 2007, 31(1): 150- 165.

    [12] 毛偉, 李玉霖, 張銅會, 趙學(xué)勇, 黃迎新, 宋琳琳. 不同尺度生態(tài)學(xué)中植物葉性狀研究概述. 中國沙漠, 2012, 32(1): 33- 41.

    [13] Díaz S, Cabido M C, Casanoves F. Plant functional traits and environmental filters at a regional scale. Journal of Vegetation Science, 1998, 9(1): 113- 122.

    [14] Díaz S, Cabido Marcelo C, Zak M, Carretero E M, Araníbar J. Plant functional traits, ecosystem structure and land-use history along a climatic gradient in central-western Argentina. Journal of Vegetation Science, 1999, 10(5): 651- 660.

    [15] Díaz S, Lavorel S, Mcintyre S, Falczuk V, Casanoves F, Milchunas D G., Skarpe C, Rusch G, Sternberg M, Noy-Meir I, Landsberg J, Zhang W, Clark H, Campbell B. Plant trait responses to grazing-a global synthesis. Global Change Biology, 2007, 13(2): 313- 341.

    [16] Cornelissen J H C, Lavorel S, Garnier E, Díaz S, Buchmann N, Gurvich D E, Reich P B, Steege H T, Morgan H D, van der Heijden M G A, Pausas J G, Poorter H. A Handbook of Protocols for Standardised and Easy Measurement of Plant Functional Traits Worldwide. Australian Journal of Botany, 2003, 51(4): 335- 380.

    [17] 肖玉, 謝高地, 安凱, 魯春霞. 基于功能性狀的生態(tài)系統(tǒng)服務(wù)研究框架. 植物生態(tài)學(xué)報, 2012, 36(4): 353- 362.

    [18] Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H C, Diemer M, Flexas J, Garnier E, Groom P K, Gulias J, Hikosaka K, Lamont B B, Lee T, Lee W, Lusk C, Midgley J J, Marie-Laure N, Niinemets ü, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov V I, Roumet C, Thomas S C, Tjoelker M G, Veneklaas E J, Villar R. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821- 827.

    [19] Wright I J, Reich P B, Cornelissen J H C, Falster D S, Garnier E, Hikosaka K, Lamont B B, Lee W, Oleksyn J, Osada N, Poorter H, Villar R, Warton D I, Westoby M. Assessing the generality of global leaf trait relationships. New Phytologist, 2005, 166(2): 485- 496.

    [20] Kattge J, Díaz S, Lavorel S, Prentice I C, Leadley P, B?nisch G, Garnier E, Westoby M, Reich P B, Wright I J, Cornelissen J H C, Violle C, Harrison S P, van Bodegom P M, Reichstein M, Enquist B J, Soudzilovskaia N A, Ackerly D D, Anand M, Atkin O, Bahn M, Baker T R, Baldocchi D, Bekker R, Blanco C C, Blonder B, Bond W J, Bradstock R, Bunker D E, Casanoves F, Cavender-Bares J, Chambers J Q, ChapinF S, Chave J, Coomes D, Cornwell W K, Craine J M, Dobrin B H, Duarte L, Durka W, Elser J, Esser G, Estiarte M, Fagan W F, Fang J, Fernández-Méndez F, Fidelis A, Finegan B, Flores O, Ford H, Frank D, Freschet G T, Fyllas N M, Gallagher R V, Green W A, Gutierrez A G, Hickler T, Higgins S I, Hodgson J G, Jalili A, Jansen S, Joly C A, Kerkhoff A J, Kirkup D, Kitajima K, Kleyer M, Klotz S, Knops J M H, Kramer K, Kühn I, Kurokawa H, Laughlin D, Lee T D, Leishman M, Lens F, Lenz T, Lewis S L, Lloyd J, Llusià J, Louault F, Ma S, Mahecha M D, Manning P, Massad T, Medlyn B E, Messier J, Moles A T, Müller S C, Nadrowski K, Naeem S, Niinemets ü, N?llert S, Nüske A, Ogaya R, Oleksyn J, Onipchenko V G, Onoda Y, Ordoez J, Overbeck G, Ozinga W A, Patio S, Paula S, Pausas J G, Peuelas J, Phillips O L, Pillar V, Poorter H, Poorter L, Poschlod P, Prinzing A, Proulx R, Rammig A, Reinsch S, Reu B, Sack L, Salgado-Negret B, Sardans J, Shiodera S, Shipley B, Siefert A, Sosinski E, Soussana F J, Swaine E, Swenson N, Thompson K, Thornton P, Waldram M, Weiher E, White M, White S, Wright S J, Yguel B, Zaehle S, Zanne A E, Wirth C. TRY-a global database of plant traits. Global Change Biology, 2011, 17(9): 2905- 2935.

    [21] 陳瑩婷, 許振柱. 植物葉經(jīng)濟(jì)譜的研究進(jìn)展. 植物生態(tài)學(xué)報, 2014, 38 (10): 1135- 1153.

    [22] Swenson N G, Enquist B J. Ecological and evolutionary determinants of a key plant functional trait: wood density and its community-wide variation across latitude and elevation. American Journal of Botany, 2007, 94(3): 451- 459.

    [23] van Kleunen M, Weber E, Fischer M. A meta-analysis of trait differences between invasive and non-invasive plant species. Ecology Letters, 2010, 13(2): 235- 245.

    [24] Sack L, Scoffoni C, Grace P J, Poorter H, Mason C M, Mendez-Alonzo R, Donovan L A. How do leaf veins influence the worldwide leaf economic spectrum? Review and synthesis. Journal of Experimental Botany, 2013, 64(13): 4053- 4080.

    [25] Mason C M, McGaughey S E., Donovan L A. Ontogeny strongly and differentially alters leaf economic and other key traits in three diverse Helianthus species. Journal of Experimental Botany, 2013, 64(13): 4089- 4099.

    [26] Price C A, Wright I J, Ackerly D D, Niinemets ü, Reich P B, Veneklaas E J. Are leaf functional traits ‘invariant’ with plant size and what is ‘invariance’ anyway? Functional Ecology, 2014, 28(6): 1330- 1343.

    [27] Atkin O K, Bloomfield K J, Reich P B, Tjoelker M G, Asner G P, Bonal D, B?nisch G, Bradford M G, Cernusak L A, Cosio E G, Creek D, Crous K Y, Domingues T F, Dukes J S, Egerton J J G, Evans J R, Farquhar G D, Fyllas N M, Gauthier P P G, Gloor E, Gimeno T E, Griffin K L, Guerrieri R, Heskel M A, Huntingford C, Ishida F Y, Kattge J, Lambers H, Liddell M J, Lloyd J, Lusk C H, Martin R E, Maksimov A P, Maximov T C, Malhi Y, Medlyn B E, Meir P, Mercado L M, Mirotchnick N, Ng D, Niinemets ü, O′Sullivan O S, Phillips O L, Poorter L, Poot P, Prentice I C, Salinas N, Rowland L M, Ryan MG, Sitch S, Slot M, Smith N G, Turnbull M H, VanderWel M C, Valladares F, Veneklaas EJ, Weerasinghe L K, Wirth C, Wright I J, Wythers K R, Xiang J, Xiang S, Zaragoza-Castells J. Global variability in leaf respiration in relation to climate, plant functional types and leaf traits. New Phytologist, 2015, 206(2): 614- 636.

    [28] 祝介東, 孟婷婷, 倪健, 蘇宏新, 謝宗強(qiáng), 張守仁, 鄭元潤, 肖春旺. 不同氣候帶間成熟林植物葉性狀間異速生長關(guān)系隨功能型的變異. 植物生態(tài)學(xué)報, 2011, 35(7): 687- 698.

    [29] Chen Y H, Han W X, Tang L Y, Tang Z Y, Fang J Y. Leaf nitrogen and phosphorus concentrations of woody plants differ in responses to climate, soil and plant growth form. Ecography, 2013, 36(2): 178- 184.

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

    [31] 周廣勝, 張新時, 高素華, 白克智, 延曉東, 鄭元潤. 中國植被對全球變化反應(yīng)的研究. 植物學(xué)報, 1997, 39(9): 879- 888.

    [32] 黎磊, 陳家寬. 氣候變化對野生植物的影響及保護(hù)對策. 生物多樣性, 2014, 22(5): 549- 563.

    [33] 衛(wèi)民, 王輝民, 王其冬, 劉允芬, 賀慶棠, 袁嘉祖, 邵海榮, 宋從和. 氣候變化對我國紅松林的影響. 地理研究, 1995, 14(1): 17- 26.

    [34] 郭泉水, 閻洪, 徐德應(yīng), 王兵. 氣候變化對我國紅松林地理分布影響的研究. 生態(tài)學(xué)報, 1998, 18(5): 484- 488.

    [35] 郝占慶, 代力民, 賀紅士, David J.Malandnoff, 邵國凡. 氣候變暖對長白山主要樹種的潛在影響. 應(yīng)用生態(tài)學(xué)報.2001, 12(5):653- 658.

    [36] 冷文芳, 賀紅士, 布仁倉, 胡遠(yuǎn)滿. 氣候變化條件下東北森林主要建群種的空間分布. 生態(tài)學(xué)報, 2006, 26(12): 4257- 4266.

    [37] 馮秋紅, 史作民, 董莉莉. 植物功能性狀對環(huán)境的響應(yīng)及其應(yīng)用. 林業(yè)科學(xué), 2008, 44(4): 125- 131.

    [38] 陳悅, 許大全. 植物光合速率對光強(qiáng)從飽和到有限轉(zhuǎn)變響應(yīng)方式的物種依賴性. 植物生理與分子生物學(xué)學(xué)報, 2007, 33(6): 538- 546.

    [39] 王凱, 朱教君, 于立忠, 孫一榮, 張金鑫. 光環(huán)境對胡桃楸幼苗生長與光合作用的影響. 應(yīng)用生態(tài)學(xué)報, 2010, 21(4): 821- 826.

    [40] 馮玉龍, 曹坤芳, 馮志立, 馬玲. 四種熱帶雨林樹種幼苗比葉重, 光合特性和暗呼吸對生長光環(huán)境的適應(yīng). 生態(tài)學(xué)報, 2002, 22(6): 901- 910.

    [41] 劉慎諤. 關(guān)于大小興安嶺的森林更新問題. 林業(yè)科學(xué), 1957, 3(3): 263- 280.

    [42] 王戰(zhàn). 對于小興安嶺紅松林更新和主伐方式的意見. 林業(yè)科學(xué), 1957, 3(3): 281- 288.

    [43] 李俊清, 朱春全, 柴一新, 張東力, 劉傳照, 姚成濱. 闊葉紅松林的營養(yǎng)結(jié)構(gòu)與動態(tài)特性. 吉林林學(xué)院學(xué)報, 1989, 5(2): 1- 16.

    [44] 李俊清, 柴一新, 張東力. 人工闊葉紅松林的結(jié)構(gòu)與生產(chǎn)力. 林業(yè)科學(xué), 1990, 26(1): 1- 8.

    [45] 李俊清, 柴一新, 朱春全, 王業(yè)蘧. 不同光強(qiáng)下紅松幼樹光合作用和營養(yǎng)物質(zhì)含量的季節(jié)模式. 生態(tài)學(xué)雜志, 1991, 10(5): 1- 5.

    [46] 劉傳照, 李俊清, 金奎剛. 林下光照條件與紅松幼樹生長的相關(guān)性研究. 東北林業(yè)大學(xué)學(xué)報, 1991, 19(3): 103- 108.

    [47] 詹鴻振, 劉傳照, 徐桂林, 程文俊, 李茹秀, 李景文. 不同庇蔭的光照特點(diǎn)及對紅松苗的影響. 東北林業(yè)大學(xué)學(xué)報, 1986, 14(2): 25- 33.

    [48] 李晶, 閻秀峰, 祖元剛. 低溫脅迫下紅松幼苗活性氧的產(chǎn)生及保護(hù)酶的變化. 植物學(xué)報, 2000, 42(2): 148- 152.

    [49] 閻秀峰, 李晶, 祖元剛. 干旱脅迫對紅松幼苗保護(hù)酶活性及脂質(zhì)過氧化作用的影響. 生態(tài)學(xué)報, 1999, 19(6): 850- 854.

    [50] 聶修和, 聶宜茂, 聶俊華, 姜黎霞. 光合有效輻射測量原理及其單位換算. 山東農(nóng)業(yè)大學(xué)學(xué)報, 1992, 23(3): 247- 258.

    [51] Vile D, Garnier é, Shipley B, Laurent G R, Marie-Laure N, Roumet C, Lavorel S, DíAz S, Hodgson J G, Lloret F, Midgley G F, Poorter H, Rutherford M C, Wilson P J, Wright I J. Specific leaf area and dry matter content estimate thickness in laminar leaves. Annals of Botany, 2005, 96(6): 1129- 1136.

    [52] 許大全. 光合作用測定及研究中一些值得注意的問題. 植物生理學(xué)通訊, 2006, 42(6): 1163- 1167.

    [53] López-Bucio J, Cruz-Ramírez A, Herrera-Estrella L. The role of nutrient availability in regulating root architecture. Current Opinion in Plant Biology, 2003, 6(3): 280- 287

    [54] Reid K E, Olsson N, Schlosser J, Peng F, Lund S T. An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biology, 2006, 6: 27.

    [55] Santiago L S, Wright S J. Leaf functional traits of tropical forest plants in relation to growth form. Functional Ecology, 2007, 21(1): 19- 27.

    [56] 孫一榮, 朱教君, 于立忠, 閆巧玲, 王凱. 不同光環(huán)境對紅松幼苗光合生理特征的影響. 生態(tài)學(xué)雜志, 2009, 28 (5): 850- 857.

    [57] 張雷明, 曹沛雨, 朱亞平, 李慶康, 張軍輝, 王曉凌, 戴冠華, 李金功. 長白山闊葉紅松林生態(tài)系統(tǒng)光能利用率的動態(tài)變化及其主控因子. 植物生態(tài)學(xué)報, 2015, 39(12): 1156- 1165.

    [58] 韓麗冬, 牟長城, 張軍輝. 透光撫育對長白山闊葉紅松林冠下紅松光合作用的影響. 東北林業(yè)大學(xué)學(xué)報, 2016, 44(4): 38- 40.

    [59] 張彌, 吳家兵, 關(guān)德新, 施婷婷, 陳鵬獅, 紀(jì)瑞鵬. 長白山闊葉紅松林主要樹種光合作用的光響應(yīng)曲線. 應(yīng)用生態(tài)學(xué)報, 2006, 17(9): 1575- 1578.

    [60] 鄭淑霞, 上官周平. 不同功能型植物光合特性及其與葉氮含量、比葉重的關(guān)系. 生態(tài)學(xué)報, 2007, 27(l): 171- 181.

    [61] Garnier E, Gobin O, Poorter H. Nitrogen productivity depends on photosynthetic nitrogen use efficiency and on nitrogen allocation within the plant. Annals of Botany, 1995, 76(6): 667- 672.

    [62] Wright I J, Reich P B, Westoby M. Strategy shifts in leaf physiology, structure and nutrient content between species of high- and low-rainfall and high- and low-nutrient habitats. Functional Ecology, 2001, 15(4): 423- 434.

    [63] Elser J J, Urabe J. The stoichiometry of consumer-driven nutrient recycling: theory, observations, and consequences. Ecology, 1999, 80(3): 735- 751.

    [64] Elser J J, Sterner R W, Gorokhova E, Fagan W F, Markow T A, Cotner J B, Harrison J F, Hobbie S E, Odell G M, Weider L J. Biological stoichiometry from genes to ecosystems. Ecology Letters, 2000, 3(6): 540- 550.

    [65] 許大全. 光合作用氣孔限制分析中的一些問題. 植物生理學(xué)通訊, 1997, 33(4): 241- 244.

    ChangesintheresponseofleaftraitsinPinuskoraiensis(Koreanpine)seedlingsofdifferentagestocontrolledtemperaturesandlightconditions

    DONG Yichen, LIU Yanhong*

    TheKeyLaboratoryforForestResourcesandEcosystemProcessesofBeijing,BeijingForestryUniversity,Beijing100083,China

    國家“十二五”科技支撐計(jì)劃研究(2012BAC01B03)

    2015- 11- 12;

    2016- 08- 24

    10.5846/stxb201511122293

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

    董伊晨, 劉艷紅.紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng).生態(tài)學(xué)報,2017,37(17):5662- 5672.

    Dong Y C, Liu Y H.Changes in the response of leaf traits inPinuskoraiensis(Korean pine) seedlings of different ages to controlled temperatures and light conditions.Acta Ecologica Sinica,2017,37(17):5662- 5672.

    猜你喜歡
    紅松林苗齡紅松
    繞口令
    沒有紅松的紅松林
    不同苗齡大花序桉造林效果分析
    不同接穗苗齡對小果型西瓜嫁接育苗效果的影響
    紅松林采伐更新的確定方式
    東折棱河楓樺紅松林與柞樹紅松林結(jié)構(gòu)特征對比分析
    苗齡、光照強(qiáng)度和施氮量對水稻營養(yǎng)物質(zhì)的影響及其與抗褐飛虱的關(guān)系
    選擇紅松寶就是選擇財(cái)富
    選擇紅松寶就是選擇財(cái)富
    高見的詩
    久久久久久大精品| 日韩av在线大香蕉| 美女午夜性视频免费| 日韩欧美国产一区二区入口| 久久久久久久午夜电影 | 精品少妇一区二区三区视频日本电影| 国产一区在线观看成人免费| 男女之事视频高清在线观看| 亚洲成国产人片在线观看| 老司机深夜福利视频在线观看| 久久婷婷成人综合色麻豆| 悠悠久久av| 嫩草影视91久久| 最好的美女福利视频网| 叶爱在线成人免费视频播放| 午夜福利在线观看吧| 国产成人欧美在线观看| 超碰成人久久| 日本欧美视频一区| 亚洲国产精品一区二区三区在线| 一级片'在线观看视频| 国产精品 欧美亚洲| 久久人妻av系列| 老司机在亚洲福利影院| 男女做爰动态图高潮gif福利片 | 真人一进一出gif抽搐免费| 三级毛片av免费| 精品国产乱子伦一区二区三区| 亚洲专区国产一区二区| 亚洲激情在线av| 男女下面插进去视频免费观看| 性少妇av在线| 成人永久免费在线观看视频| 国产精品久久久av美女十八| 久久精品国产亚洲av高清一级| 一进一出抽搐动态| 欧美日韩国产mv在线观看视频| 国产欧美日韩一区二区三区在线| www.999成人在线观看| 亚洲一区二区三区欧美精品| 校园春色视频在线观看| 脱女人内裤的视频| 成人特级黄色片久久久久久久| 免费在线观看视频国产中文字幕亚洲| 最好的美女福利视频网| 亚洲av日韩精品久久久久久密| 一区二区三区精品91| 成人国语在线视频| 一区二区三区国产精品乱码| 色哟哟哟哟哟哟| 久久狼人影院| 精品国产乱码久久久久久男人| 欧美乱妇无乱码| 自拍欧美九色日韩亚洲蝌蚪91| 日日爽夜夜爽网站| 国产一区二区三区视频了| 亚洲美女黄片视频| 久久人妻熟女aⅴ| 成人黄色视频免费在线看| 老司机福利观看| 啦啦啦在线免费观看视频4| 国产乱人伦免费视频| 韩国精品一区二区三区| 啦啦啦 在线观看视频| 久久亚洲真实| www.精华液| 午夜免费鲁丝| 麻豆成人av在线观看| 久久亚洲精品不卡| 亚洲国产精品一区二区三区在线| 性欧美人与动物交配| 在线播放国产精品三级| 国产精品98久久久久久宅男小说| 国产一区在线观看成人免费| 少妇裸体淫交视频免费看高清 | 亚洲精品粉嫩美女一区| 日韩一卡2卡3卡4卡2021年| 久久热在线av| 999久久久国产精品视频| 国产蜜桃级精品一区二区三区| 成在线人永久免费视频| 精品日产1卡2卡| 黄色怎么调成土黄色| 亚洲av五月六月丁香网| 久久中文看片网| 一个人观看的视频www高清免费观看 | www日本在线高清视频| 纯流量卡能插随身wifi吗| av视频免费观看在线观看| 亚洲七黄色美女视频| 日本wwww免费看| 青草久久国产| 午夜久久久在线观看| 9191精品国产免费久久| 久久久久亚洲av毛片大全| 十八禁网站免费在线| 亚洲五月婷婷丁香| 91大片在线观看| 97人妻天天添夜夜摸| 国产精品影院久久| 国产伦人伦偷精品视频| aaaaa片日本免费| 免费在线观看日本一区| 热re99久久精品国产66热6| 一区二区日韩欧美中文字幕| 老鸭窝网址在线观看| 成年女人毛片免费观看观看9| 亚洲国产欧美网| 老汉色∧v一级毛片| 18禁国产床啪视频网站| 大香蕉久久成人网| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| a级片在线免费高清观看视频| a级毛片在线看网站| 91麻豆av在线| 亚洲自偷自拍图片 自拍| 久久欧美精品欧美久久欧美| 成在线人永久免费视频| 精品久久蜜臀av无| 亚洲av片天天在线观看| 妹子高潮喷水视频| 亚洲成av片中文字幕在线观看| 成人黄色视频免费在线看| 最近最新中文字幕大全电影3 | 高清毛片免费观看视频网站 | 亚洲 欧美一区二区三区| 一级a爱片免费观看的视频| 好看av亚洲va欧美ⅴa在| 日日摸夜夜添夜夜添小说| 无人区码免费观看不卡| 国产亚洲欧美98| aaaaa片日本免费| 人妻久久中文字幕网| 免费av中文字幕在线| 1024视频免费在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品久久久久5区| 美女高潮喷水抽搐中文字幕| 99精国产麻豆久久婷婷| 在线观看免费视频网站a站| 中文字幕人妻丝袜制服| 久久久久久亚洲精品国产蜜桃av| 丰满饥渴人妻一区二区三| 国产一卡二卡三卡精品| bbb黄色大片| 精品国内亚洲2022精品成人| 一进一出抽搐gif免费好疼 | 在线观看午夜福利视频| 亚洲在线自拍视频| 亚洲精品成人av观看孕妇| 80岁老熟妇乱子伦牲交| 亚洲人成77777在线视频| 久久久久九九精品影院| 欧美乱码精品一区二区三区| 99国产极品粉嫩在线观看| 国产精品98久久久久久宅男小说| 婷婷精品国产亚洲av在线| 美女福利国产在线| 日韩免费高清中文字幕av| 男女高潮啪啪啪动态图| 日韩一卡2卡3卡4卡2021年| 操出白浆在线播放| 久久这里只有精品19| 97超级碰碰碰精品色视频在线观看| 日本撒尿小便嘘嘘汇集6| 99在线视频只有这里精品首页| 男人操女人黄网站| 国产极品粉嫩免费观看在线| 黄片大片在线免费观看| 亚洲黑人精品在线| 男女之事视频高清在线观看| 操出白浆在线播放| 欧美乱码精品一区二区三区| 久久国产精品男人的天堂亚洲| 欧美激情高清一区二区三区| 欧美成人免费av一区二区三区| 成人永久免费在线观看视频| 国产欧美日韩一区二区三| 老鸭窝网址在线观看| 欧美另类亚洲清纯唯美| 国产精品1区2区在线观看.| 欧美在线黄色| 国产精品美女特级片免费视频播放器 | 欧美日韩乱码在线| 成人影院久久| 日韩免费av在线播放| 欧美成人性av电影在线观看| 搡老乐熟女国产| 91字幕亚洲| 99精国产麻豆久久婷婷| 久久久国产一区二区| 91在线观看av| 日本欧美视频一区| 国产精华一区二区三区| 人妻丰满熟妇av一区二区三区| 人人妻,人人澡人人爽秒播| 久久中文字幕一级| 午夜成年电影在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲av熟女| 亚洲第一欧美日韩一区二区三区| 激情视频va一区二区三区| 在线永久观看黄色视频| 又大又爽又粗| 国产精品秋霞免费鲁丝片| 免费av中文字幕在线| 真人一进一出gif抽搐免费| 高清毛片免费观看视频网站 | 国产主播在线观看一区二区| 乱人伦中国视频| 涩涩av久久男人的天堂| ponron亚洲| 嫩草影视91久久| 嫩草影视91久久| 久久人人精品亚洲av| 精品人妻在线不人妻| 亚洲视频免费观看视频| 精品乱码久久久久久99久播| 日韩视频一区二区在线观看| 大型av网站在线播放| 大香蕉久久成人网| 夜夜爽天天搞| 国产精品久久视频播放| 人成视频在线观看免费观看| 国产成人影院久久av| 欧美在线一区亚洲| 欧美乱妇无乱码| 女人精品久久久久毛片| 亚洲欧美精品综合久久99| 国产成人av激情在线播放| 亚洲 欧美一区二区三区| 欧美日韩瑟瑟在线播放| 激情视频va一区二区三区| 亚洲精品一区av在线观看| 久久精品aⅴ一区二区三区四区| 一二三四社区在线视频社区8| 国产成人av教育| 中文字幕高清在线视频| 成人三级黄色视频| 香蕉国产在线看| 国产成人精品久久二区二区免费| 亚洲国产精品sss在线观看 | 性欧美人与动物交配| 亚洲精品国产区一区二| 久久人妻熟女aⅴ| 精品久久久久久成人av| 婷婷精品国产亚洲av在线| 国产成人影院久久av| 中文字幕av电影在线播放| 免费看a级黄色片| 久久99一区二区三区| 嫁个100分男人电影在线观看| 波多野结衣一区麻豆| 久久精品国产亚洲av高清一级| 国产有黄有色有爽视频| 亚洲黑人精品在线| 成人黄色视频免费在线看| 欧美丝袜亚洲另类 | 欧美成人免费av一区二区三区| 一本综合久久免费| 亚洲欧美激情在线| 国产成人系列免费观看| www.www免费av| 中文字幕av电影在线播放| 久久国产乱子伦精品免费另类| 母亲3免费完整高清在线观看| 9色porny在线观看| av天堂在线播放| 99精国产麻豆久久婷婷| 国产成人啪精品午夜网站| 亚洲熟女毛片儿| 色综合站精品国产| 午夜成年电影在线免费观看| 日本免费一区二区三区高清不卡 | 亚洲精华国产精华精| 人人妻人人添人人爽欧美一区卜| 欧美日本亚洲视频在线播放| 在线观看免费日韩欧美大片| 在线观看舔阴道视频| 亚洲国产看品久久| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美一区二区三区久久| avwww免费| 中文欧美无线码| 香蕉久久夜色| 久久天堂一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 12—13女人毛片做爰片一| 叶爱在线成人免费视频播放| 我的亚洲天堂| 在线视频色国产色| 欧美日韩一级在线毛片| 欧美日韩精品网址| 国产成年人精品一区二区 | 丰满迷人的少妇在线观看| 国产精品影院久久| 午夜免费激情av| 美女福利国产在线| 亚洲第一青青草原| 免费看a级黄色片| 久久久久精品国产欧美久久久| 丁香欧美五月| 99在线视频只有这里精品首页| 国产成+人综合+亚洲专区| 国产高清激情床上av| 高清毛片免费观看视频网站 | 色综合欧美亚洲国产小说| 九色亚洲精品在线播放| 在线永久观看黄色视频| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 免费高清在线观看日韩| 一进一出抽搐动态| 日韩免费av在线播放| 啪啪无遮挡十八禁网站| 日日摸夜夜添夜夜添小说| 国产xxxxx性猛交| 大香蕉久久成人网| 亚洲自偷自拍图片 自拍| 另类亚洲欧美激情| 国产亚洲欧美精品永久| 成人国产一区最新在线观看| 黄频高清免费视频| 亚洲avbb在线观看| 日本黄色日本黄色录像| 嫩草影院精品99| 大香蕉久久成人网| 亚洲午夜理论影院| 99热只有精品国产| 色综合欧美亚洲国产小说| 国产成人欧美| 精品福利永久在线观看| netflix在线观看网站| 国产精品免费一区二区三区在线| 国产欧美日韩综合在线一区二区| 动漫黄色视频在线观看| 18禁美女被吸乳视频| 久久午夜综合久久蜜桃| 亚洲 欧美 日韩 在线 免费| 五月开心婷婷网| 亚洲国产中文字幕在线视频| 91字幕亚洲| 久久精品91无色码中文字幕| 一区二区三区国产精品乱码| 午夜福利一区二区在线看| 热99re8久久精品国产| 亚洲专区国产一区二区| 男女之事视频高清在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲成av片中文字幕在线观看| 久久久国产一区二区| 纯流量卡能插随身wifi吗| 国产精品久久视频播放| 欧美 亚洲 国产 日韩一| 国产午夜精品久久久久久| 免费久久久久久久精品成人欧美视频| 久久久久亚洲av毛片大全| 国产黄色免费在线视频| 香蕉久久夜色| 亚洲av美国av| 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 免费观看人在逋| 淫妇啪啪啪对白视频| 精品一品国产午夜福利视频| 桃色一区二区三区在线观看| 一进一出抽搐gif免费好疼 | 亚洲精品国产一区二区精华液| 久热爱精品视频在线9| 亚洲精品国产精品久久久不卡| 中文字幕色久视频| 久久香蕉精品热| 夜夜爽天天搞| 不卡av一区二区三区| 琪琪午夜伦伦电影理论片6080| 欧洲精品卡2卡3卡4卡5卡区| 欧美在线黄色| 夜夜夜夜夜久久久久| 不卡一级毛片| 视频区图区小说| 天天添夜夜摸| 久久国产精品人妻蜜桃| 美女午夜性视频免费| 老司机午夜十八禁免费视频| 丰满饥渴人妻一区二区三| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 午夜福利在线免费观看网站| 怎么达到女性高潮| 国产xxxxx性猛交| 一进一出抽搐动态| 少妇粗大呻吟视频| 黄色片一级片一级黄色片| 一二三四在线观看免费中文在| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| 天堂动漫精品| 中文字幕最新亚洲高清| 亚洲久久久国产精品| 夫妻午夜视频| 天天影视国产精品| 久久久国产一区二区| 久久久国产欧美日韩av| 久久久国产一区二区| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 99久久人妻综合| tocl精华| 亚洲精品一二三| 亚洲黑人精品在线| 他把我摸到了高潮在线观看| 午夜福利在线免费观看网站| 大香蕉久久成人网| 怎么达到女性高潮| 女同久久另类99精品国产91| 国产亚洲精品第一综合不卡| 一a级毛片在线观看| 黄色视频不卡| 十八禁人妻一区二区| xxx96com| 亚洲精品一卡2卡三卡4卡5卡| 大型黄色视频在线免费观看| 久久精品影院6| 欧美激情高清一区二区三区| 免费不卡黄色视频| 欧美人与性动交α欧美软件| 免费少妇av软件| 夜夜爽天天搞| 亚洲情色 制服丝袜| 色哟哟哟哟哟哟| a级片在线免费高清观看视频| 真人一进一出gif抽搐免费| 精品国产国语对白av| 91成年电影在线观看| 窝窝影院91人妻| 亚洲视频免费观看视频| 国产成人啪精品午夜网站| 久久婷婷成人综合色麻豆| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 亚洲男人天堂网一区| 69精品国产乱码久久久| 成人三级黄色视频| 色综合站精品国产| 免费在线观看亚洲国产| 丝袜美足系列| 久热这里只有精品99| 视频在线观看一区二区三区| 久久久国产一区二区| 亚洲中文字幕日韩| 久久精品影院6| 国产精品美女特级片免费视频播放器 | 午夜福利在线免费观看网站| 日日爽夜夜爽网站| 日本精品一区二区三区蜜桃| 国产高清videossex| 国产真人三级小视频在线观看| 国产精华一区二区三区| 九色亚洲精品在线播放| 午夜福利在线观看吧| 九色亚洲精品在线播放| 亚洲精华国产精华精| 欧美日韩av久久| 国产亚洲欧美在线一区二区| 美女福利国产在线| 午夜福利在线免费观看网站| 欧美 亚洲 国产 日韩一| 日本免费a在线| 欧美乱色亚洲激情| 免费观看人在逋| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 国产亚洲精品综合一区在线观看 | 亚洲av成人av| 国产高清国产精品国产三级| 国产又色又爽无遮挡免费看| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看一区二区三区激情| 久久久久久亚洲精品国产蜜桃av| 无遮挡黄片免费观看| 性少妇av在线| 国产99白浆流出| 嫁个100分男人电影在线观看| 在线观看66精品国产| 麻豆av在线久日| 国产精品99久久99久久久不卡| 欧美在线一区亚洲| 成年版毛片免费区| 日韩人妻精品一区2区三区| 国产乱人伦免费视频| 男女下面插进去视频免费观看| 国产精品国产高清国产av| 日本免费a在线| 亚洲va日本ⅴa欧美va伊人久久| 天天添夜夜摸| 亚洲一区二区三区不卡视频| www.www免费av| 亚洲精品一卡2卡三卡4卡5卡| 免费看十八禁软件| 亚洲中文字幕日韩| tocl精华| 嫩草影视91久久| 久久这里只有精品19| 精品一品国产午夜福利视频| 国产激情欧美一区二区| 黄频高清免费视频| 国产99白浆流出| 女人精品久久久久毛片| 国产一卡二卡三卡精品| 91在线观看av| 丝袜人妻中文字幕| 嫩草影视91久久| 欧美中文综合在线视频| 国产一区二区激情短视频| www日本在线高清视频| 欧美乱妇无乱码| 午夜影院日韩av| 国产高清国产精品国产三级| 日韩视频一区二区在线观看| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 一进一出好大好爽视频| 亚洲国产欧美网| 中亚洲国语对白在线视频| 亚洲专区中文字幕在线| 99精品久久久久人妻精品| 一进一出抽搐gif免费好疼 | 精品福利观看| 啪啪无遮挡十八禁网站| 亚洲精品一卡2卡三卡4卡5卡| 这个男人来自地球电影免费观看| 九色亚洲精品在线播放| 中亚洲国语对白在线视频| 久久香蕉精品热| 黄片大片在线免费观看| 精品一区二区三区视频在线观看免费 | 两个人免费观看高清视频| 精品国产乱子伦一区二区三区| 久久久久久久午夜电影 | 黄色视频,在线免费观看| 老司机深夜福利视频在线观看| 成年人黄色毛片网站| 大型黄色视频在线免费观看| 多毛熟女@视频| 日韩精品免费视频一区二区三区| 一级毛片精品| 国产一区二区三区在线臀色熟女 | 国产一区二区三区综合在线观看| 女性被躁到高潮视频| 一级,二级,三级黄色视频| 免费在线观看影片大全网站| 在线av久久热| 美女午夜性视频免费| 日日爽夜夜爽网站| 精品国产乱码久久久久久男人| 国产免费男女视频| 亚洲欧美一区二区三区黑人| 欧美日韩亚洲高清精品| 夜夜夜夜夜久久久久| 美女午夜性视频免费| xxx96com| 国产伦人伦偷精品视频| 亚洲精品成人av观看孕妇| 亚洲人成电影观看| 村上凉子中文字幕在线| 多毛熟女@视频| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久成人av| 18禁美女被吸乳视频| 国产高清视频在线播放一区| 亚洲欧美激情综合另类| 黄色女人牲交| 看免费av毛片| 天天添夜夜摸| 精品国产超薄肉色丝袜足j| 精品无人区乱码1区二区| 久久国产精品男人的天堂亚洲| 一级毛片精品| 亚洲九九香蕉| 亚洲午夜精品一区,二区,三区| av福利片在线| 久久人人精品亚洲av| 免费av中文字幕在线| 久久中文字幕人妻熟女| 手机成人av网站| 极品人妻少妇av视频| 美女大奶头视频| 一级毛片精品| 在线看a的网站| 精品一区二区三区视频在线观看免费 | 欧美日韩乱码在线| 亚洲专区国产一区二区| 首页视频小说图片口味搜索| 他把我摸到了高潮在线观看| 欧美人与性动交α欧美软件| 这个男人来自地球电影免费观看| 国产野战对白在线观看| 夫妻午夜视频| 人人妻人人爽人人添夜夜欢视频| 日日爽夜夜爽网站| 欧美成人性av电影在线观看| 天堂俺去俺来也www色官网| 免费在线观看亚洲国产| 精品一品国产午夜福利视频| 亚洲精品国产色婷婷电影| 亚洲第一青青草原| 国产免费男女视频| 超色免费av| 12—13女人毛片做爰片一| 国产片内射在线| 亚洲七黄色美女视频| 久久久久久人人人人人| 中文字幕人妻熟女乱码| 男女午夜视频在线观看|