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

    不同種源杉木樹輪α纖維素δ13C對年氣候因子的響應(yīng)

    2020-04-21 03:59:26劉小英段愛國張建國張雄清朱安明
    林業(yè)科學(xué)研究 2020年2期
    關(guān)鍵詞:樹輪氣候因子最低氣溫

    劉小英,段愛國,2*,張建國,2,張雄清,朱安明

    (1. 中國林業(yè)科學(xué)研究院林業(yè)研究所,國家林業(yè)和草原局林木培育重點實驗室,北京 100091;2. 南京林業(yè)大學(xué)南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,江蘇 南京 210000;3. 中國林業(yè)科學(xué)研究院林業(yè)科技信息研究所,北京 100091)

    造林樹種種源選擇已成為人類應(yīng)對氣候的重要途徑,探討不同地理種源主要造林樹種樹輪指標對氣候變化的響應(yīng)及評價未來氣候變化對林木種源生長的影響已成為亟待解決的重要科學(xué)問題。林木種源是指取得種子或其他繁殖材料的地理來源或原產(chǎn)地[1]。地理種源的變異是樹木長期氣候環(huán)境適應(yīng)性進化的產(chǎn)物,從遺傳角度看,不同地理種源經(jīng)過自然選擇基因型組成會發(fā)生變化,形成固有的遺傳基礎(chǔ)[2-4],因此,在引種栽培時往往采取“就近原則”,選擇同一種源或同一種源區(qū)的幼苗進行栽植[5]。然而,有研究者通過對不同種源山毛櫸(Fagus sylvaticaL.)應(yīng)對氣候變化的響應(yīng)情況進行研究,發(fā)現(xiàn)不同種源山毛櫸對異地(栽植地)的氣候具有不同的響應(yīng)特征,而最適宜種源并不滿足“就近原則”[6-7]。

    樹輪穩(wěn)定碳同位素(δ13C)技術(shù)被廣泛應(yīng)用于生理生態(tài)、氣候變化、水文等研究領(lǐng)域,現(xiàn)已逐步成為了解生物與其生存環(huán)境相互關(guān)系的工具;而樹木是如何應(yīng)對氣候變化依舊被氣候?qū)W家關(guān)注,與樹輪寬度、木材密度相似,樹輪穩(wěn)定碳同位素指標對不同尺度上氣候變化的響應(yīng)存在很大的差異[8]。目前,針對氣候響應(yīng)研究的氣候因子主要有平均氣溫[9]、降水[10]、相對濕度[11-12]、最高氣溫、最低氣溫[13]、最小相對濕度及日照時數(shù)[6]。有研究表明,樹輪穩(wěn)定碳同位素比樹輪寬度指標對氣候變化的響應(yīng)更靈敏。在種源層面上,樹木種源對環(huán)境的影響具有差異性,研究發(fā)現(xiàn),不同種源的白刺(Nitraria tangutorumBobr.)對不同梯度的水分響應(yīng)存在較大的區(qū)別[14];不同種源的落葉松(Larix gmelinii(Rupr.) Kuzen.)對原產(chǎn)地環(huán)境條件的響應(yīng)也出現(xiàn)差異[15]。不同種源樹木樹輪的早材、晚材、整輪寬度對年平均氣溫和年最高氣溫較敏感[16]。杉木(Cunninghamia lancealata(Lamb.)Hook.)是我國重要的鄉(xiāng)土針葉用材樹種[17],其生長呈現(xiàn)明顯的地理緯度變異,針對杉木樹輪指標在種源層面上響應(yīng)氣候變化的研究相對較少,而將高靈敏度的樹輪穩(wěn)定碳同位素指標引入不同種源樹木生長對氣候變化響應(yīng)的研究國內(nèi)外未見報道。

    鑒于此,本研究選擇設(shè)置于廣西柳州的34年生杉木種源試驗林中的30個代表性地理種源,探討其樹輪α纖維素中穩(wěn)定碳同位素與年際氣候因子間的關(guān)系,以期篩選出關(guān)鍵氣候因子,為未來氣候變化條件下種源層面上區(qū)域杉木的培育提供科學(xué)支撐。

    1 研究區(qū)自然概況

    杉木種源試驗林位于廣西柳州六峰山林場,由中國林科院林業(yè)研究所主持,全國杉木地理種源試驗協(xié)作組于20世紀70年代末開展的第2次杉木全分布區(qū)種源試驗時期布設(shè)[18]。依據(jù)杉木產(chǎn)區(qū)的劃分,其屬于杉木的中心產(chǎn)區(qū)。試驗林于1981年春季,采用1年實生苗造林。種源試驗林所處林地較平坦,隨機區(qū)組設(shè)計。造林行間距2 m×2 m,每個小區(qū)4株,重復(fù)10次,進行了1次疏伐,至2013年秋季,杉木種源試驗林保留密度為1 261株·hm-2,平均胸徑20.14 cm,平均樹高13.45 m。試驗林地處桂中丘陵地區(qū) (23°42′ N,109°50′ E),年平均氣溫21.1℃,年降水量1 418.5 mm,年蒸發(fā)量1 969.1 mm,相對濕度76%,全年日照時數(shù)1 823.7 h。

    2 研究材料與方法

    2.1 樹輪資料

    結(jié)合種源林商業(yè)采伐,對所有種源保留單株在20 cm高度處截取了圓盤,考慮地理分布的代表性及圓盤數(shù)量大于10個,選取貴州六枝(GZl)、湖南雙牌(HNs)以及福建永安(FJy)等30個地理種源各11個單株,合計330個單株進行分析,其中,中心產(chǎn)區(qū)10個種源,一般產(chǎn)區(qū)5個種源,邊緣產(chǎn)區(qū)15個種源。平均胸徑20.33 cm,平均樹高13.57 m。應(yīng)用LignoStation年輪分析儀對不同種源各單株圓盤進行掃描,分別以樹木的南、北2個方向進行測量,應(yīng)用年輪分析儀自帶軟件進行年輪界線的劃分,并量測樹輪自髓心往外各年輪寬度,后借助交叉定年軟件COFECHA程序?qū)λ卸ê媚甑男蛄凶鳈z驗。

    2.2 α-纖維素提取

    將樹輪圓盤放在干凈的玻璃板上,選擇向陽的方位[19],在顯微鏡下用大頭針標注經(jīng)交叉定年的樹輪。考慮到δ13C值會受到幼齡效應(yīng)的影響[20],在取樣時去除前6年的樹輪。利用干凈的刻刀進行樹輪的逐年剝離,將同一種源的11個圓盤同一年的年輪樣品剝?nèi)『蠡旌显谝黄?,放入高速離心球磨儀將樣品粉碎至過60目篩(300 μm)。按照有機溶劑(苯:無水乙醇=2∶1,丙酮180 mL)萃取、漂白(次氯酸鈉∶乙酸=3∶2)及堿化(10%與17%的氫氧化鈉依次處理)等步驟提取α-纖維素[21-22]。用真空冷凍干燥機(-45℃)進行凍干處理,采用快速核酸提取儀進行α-纖維素粉末的混合。

    2.3 碳穩(wěn)定同位素(δ13C)的測量

    采用儀器Thermo Finnigan DELTA plus XP氣象同位素質(zhì)譜儀測量杉木全木以及α-纖維素的δ13C的值,該過程在植被與變化國家重點實驗室完成。用百萬分之一天平(Starorius)稱量粉末α-纖維素(0.20±0.05)mg,用錫囊將樣品包裹,卷好密封后上機測試,得到13C/12C值(相對VPDB標準),即δ13C。整個碳同位素(δ13C)測量精度≤0.15‰。

    式中:Rsample表示樣品中的重輕同位素峰度之比,Rstandard是國際通用標準的重輕同位素峰度之比[23]。

    2.4 氣候數(shù)據(jù)來源

    圖 1 廣西來賓氣象站月平均溫度與總降水量變化(1985-2013)Fig. 1 The change trend of the mean temperature and year total precipitation of each month of the climate station of Laibin, Guangxi (1985-2013)

    本實驗采用廣西來賓氣象站點的數(shù)據(jù)(圖1),該站點距離種源試驗林場100 km,中間無高山阻隔。所用氣象資料包括1985—2013年的年平均氣溫(T)、年最高氣溫(Ta)、年最低氣溫(Ti)、年降水量(P),年平均相對濕度(RH)、年最小相對濕度(RHi)和日照時數(shù)(S)(氣候資料來源于國家氣象信息中心來賓氣象站點的數(shù)據(jù))。經(jīng)檢驗,來賓站點氣象數(shù)據(jù)可靠,無明顯突變,該站點氣候數(shù)據(jù)代表了當?shù)貧夂蜃兓厔荨?/p>

    3 結(jié)果與分析

    3.1 不同種源杉木δ13C的統(tǒng)計特征

    對30個種源交叉定年后,選用去除大氣CO2濃度的年表作為后續(xù)分析年表,年表統(tǒng)計特征包括各種源樹輪δ13C的最小值、最大值、平均值、標準差、一階自相關(guān)系數(shù)、二階自相關(guān)系數(shù)及平均敏感度。

    從杉木樹輪α-纖維素δ13C的測定結(jié)果可知:1987—2013年杉木樹輪δ13C的變幅為-20.997‰~-25.910‰(表1)。δ13C序列隨時間呈波動變化,所有種源的δ13C平均值達-24.022‰,標準差變幅為0.24‰~0.47‰??傮w上,不同種源樹輪樣品δ13C均值比較離散,為基于樹輪穩(wěn)定碳同位素探討不同種源對氣候變化的響應(yīng)提供了必要條件。

    通過對不同地理種源杉木年輪 α-纖維素δ13C序列進行自相關(guān)分析(表1)發(fā)現(xiàn):不同地理種源杉木的一階、二階自相關(guān)系數(shù)較大,一階自相關(guān)系數(shù)為 0.181~0.866,而二階自相關(guān)系數(shù)為0.044~0.749,且大部分種源的一階和二階自相關(guān)系數(shù)超過0.4。由此可以說明,大部分種源杉木樹輪α-纖維素中δ13C值在響應(yīng)氣候變化時存在明顯“滯后效應(yīng)”,故選擇當年、前1年、前2年的氣候因子進行相關(guān)分析。

    3.2 不同種源杉木樹輪δ13C值與氣候因子相關(guān)性分析

    3.2.1 杉木樹輪δ13C與年氣象因子的相關(guān)分析從圖2可看出:大部分杉木種源樹輪α-纖維素中δ13C與年平均氣溫、年最高氣溫、年最低氣溫呈顯著正相關(guān),與年平均相對濕度、年最小相對濕度呈顯著負相關(guān),而與年降水、日照時數(shù)及濕潤度指數(shù)不顯著相關(guān)。這是因為樹輪α纖維素中δ13C的累積主要通過植物的光合作用進行,而溫度的升高,有利于光合作用過程中酶活性的提高,能更有效的吸收大氣中的CO2,有效促進光合作用,從而使樹輪穩(wěn)定碳同位素可以更靈敏的響應(yīng)氣候的變化。杉木生長需要的濕度相對較大,當植物葉片處于干燥的環(huán)境時,水分的減少導(dǎo)致氣孔導(dǎo)度降低,限制了葉片中CO2的供給,從而樹輪α纖維素中δ13C增大,致使樹輪穩(wěn)定碳同位素可以更靈敏的響應(yīng)氣候的變化。

    不同產(chǎn)區(qū)杉木樹輪的δ13C對氣候因子的響應(yīng)靈敏程度存在差異。圖2表明:對于湖南安化、福建南平、福建崇安、廣西三江、廣東樂昌、廣西融水、湖南祁陽、江西遂川、湖南雙牌、廣東始興等杉木種源而言,不同種源樹輪的δ13C與當年的最小相對濕度和前2年的平均溫度、最低氣溫、最高氣溫及平均相對濕度的相關(guān)系數(shù)較高,其中,與前2年的平均相對濕度呈顯著負相關(guān)(p<0.05),與當年的最小相對濕度呈極顯著負相關(guān)(p<0.01)。對于四川汶川、浙江龍泉、江西修水、四川古藺、云南西疇等杉木種源而言,不同種源樹輪δ13C與當年的最低氣溫、前1年的最小相對濕度、平均氣溫及前2年的平均相對濕度、最高氣溫的相關(guān)系數(shù)較高,其中,與前1年最小相對濕度呈極顯著負相關(guān)(p<0.01),與當年的平均相對濕度呈顯著負相關(guān)(p<0.05);而廣東紫金、貴州六枝、湖南新寧、江蘇句容、廣西博白、福建永安、湖北谷城、陜西南鄭、云南屏邊、湖北羅田、河南商城、安徽東至、貴州麻江、浙江麗水、福建武平的杉木種源,不同種源樹輪δ13C值與當年的最小相對濕度、平均相對濕度和前2年的平均氣溫、最低氣溫、最高氣溫的相關(guān)系數(shù)較大,其中,與當年的最小相對濕度和平均相對濕度呈極顯著負相關(guān)(p<0.01),且與當年的最小相對濕度相關(guān)系數(shù)較大。不同產(chǎn)區(qū)原產(chǎn)地的氣候、地貌、土壤等立地條件均不相同,由于受到不同因素的影響,不同杉木種源的遺傳基礎(chǔ)出現(xiàn)差異,當將其引種到同一氣候條件下進行栽培時,這種差異性導(dǎo)致了不同種源對栽植地氣候因子的差異性響應(yīng)。

    表 1 1987-2013年不同種源杉木樹輪穩(wěn)定碳同位素(δ13C)序列的基本統(tǒng)計特征Table 1 Basic statistical characteristics of stable carbon isotope sequences of different tree provenances in China from 1987 to 2013

    不同地理種源杉木樹輪的δ13C與當年降水呈不顯著相關(guān)(圖2)。種源FJw與當年最低氣溫呈現(xiàn)出負相關(guān)性,而大多數(shù)種源與當年、前1年及前2年最低氣溫呈顯著或極顯著正相關(guān),相關(guān)系數(shù)最大達0.677(p<0.01)。對最小相對濕度的相關(guān)性分析中發(fā)現(xiàn),除種源Ync相關(guān)性不顯著外,其余種源均與當年、前1年及前2年最小相對濕度呈顯著或極顯著負相關(guān),相關(guān)系數(shù)最大達-0.744(p<0.01)。與日照時數(shù)的相關(guān)性分析中發(fā)現(xiàn),僅少數(shù)種源與當年、前1年及前2年呈顯著正相關(guān),相關(guān)系數(shù)最大達0.457(p<0.01)。不同地理種源杉木樹輪δ13C與平均氣溫的相關(guān)性分析中發(fā)現(xiàn)(圖2),大部分種源與當年、前1年及前2年呈顯著或極顯著正相關(guān),而種源GXb與平均氣溫呈顯著負相關(guān),相關(guān)系數(shù)可達-0.418(p<0.05)。多數(shù)種源與當年、前1年及前2年平均相對濕度呈顯著或極顯著負相關(guān),而種源GDz、GXs及Ynb與前1年呈顯著或極顯著正相關(guān)。不同地理種源杉木樹輪δ13C與前2年最高氣溫的相關(guān)性分析中發(fā)現(xiàn)(圖2),種源GXs與其呈現(xiàn)極顯著正相關(guān),相關(guān)系 數(shù) 0.569(p<0.01) , 種 源 SXn、 FJw、 HBt、HNq、HNs、GDs、GZm、Ynb、SCc與最高氣溫呈現(xiàn)出顯著的正相關(guān)關(guān)系,相關(guān)系數(shù)可達0.447(p<0.05)。

    圖 2 杉木樹輪δ13C與年際氣候因子相關(guān)分析Fig. 2 Correlation analysis between δ13C and annual climatic factors of fir tree rings

    3.2.2 不同地理種源杉木樹輪δ13C與當年氣候因子相關(guān)性隨經(jīng)緯度的變化 從杉木地理種源樹輪δ13C與當年氣候因子相關(guān)性隨經(jīng)度的變化趨勢(圖3)可知:杉木種源樹輪δ13C與年降水的相關(guān)性受經(jīng)度影響較?。╮=0.028),與年平均氣溫、最低氣溫、最高氣溫、日照時數(shù)的相關(guān)性隨經(jīng)度的增加總體上呈先增強后減弱的變化規(guī)律,且在110~115°E間正相關(guān)性最強。杉木種源樹輪δ13C與平均相對濕度及最小相對濕度的負相關(guān)性隨經(jīng)度的增大總體呈先增強后減弱的趨勢,且在110~115° E間呈現(xiàn)最強負相關(guān)性,δ13C與最小相對濕度的負相關(guān)性隨經(jīng)度的這種先增強后減弱的變動趨勢可由二次多項式描述,且達顯著水平(r=-0.462)。從穩(wěn)定碳同位素與各年際氣候因子的相關(guān)系數(shù)隨緯度的變化趨勢圖(圖3)可以發(fā)現(xiàn),杉木種源樹輪δ13C與年降水量的相關(guān)系數(shù)受種源所處緯度的影響較大,隨緯度的增加呈先上升后下降的顯著變化趨勢(r=0.434),而不同種源樹輪δ13C與年平均氣溫、最低氣溫、最高氣溫、平均相對濕度、日照時數(shù)、最小相對濕度的相關(guān)關(guān)系則受緯度影響較小。

    3.3 不同種源杉木樹輪α-纖維素中δ13C與氣候因子的多元回歸分析

    圖 3 δ13C與當年各年際氣候因子的相關(guān)系數(shù)隨著經(jīng)(緯)度的變化散點圖Fig. 3 Scatter diagram of correlation coefficients between the δ13C value and inter-annual climate factors along with longitude and latitude of provenances

    為了更深入地了解氣候因子對不同種源杉木樹輪δ13C的影響情況,以樹輪α纖維素δ13C為因變量,自變量選取相關(guān)性較強的當年最小相對濕度(X1)及前2年的平均氣溫(X2)、最低氣溫(X3)、最高氣溫(X4)和平均相對濕度(X5),回歸方法選擇逐步回歸法,置信度<0.05,輸出結(jié)果見表2。

    從表中2中α-纖維素δ13C序列的逐步回歸分析可知:在剔除了對α-纖維素中δ13C序列影響不顯著的因子后,HNh、FJc、GXs、GDc、GXr、HNq、JXs、HNs、GDs、ZJq、JXx、GDz、GZl、HNx、JCr、FJy、HBc、SXn、Ynb、HBt、HeN、Ahd、GZm、ZJs、FJw等25個種源樹輪α-纖維素δ13C序列受氣候因子當年最小相對濕度影響顯著,其中,部分種源同時受到前2年平均相對濕度或前2年最低氣溫的顯著影響,R2處于0.176~0.627之間,說明以當年最小相對濕度為主的氣候因子對杉木不同地理種源樹輪α-纖維素中δ13C序列的貢獻量為17.6%~62.7%;而福建南平(FJn)和四川汶川(SCc)種源均分別與前2年的最低氣溫和平均相對濕度呈顯著正線性及負線性相關(guān),廣西博白(GXb)種源與前2年平均氣溫呈顯著負線性相關(guān),云南西疇(Ync)、四川古藺(SCl)種源與前2年平均相對濕度具顯著負線性相關(guān)關(guān)系。

    表 2 δ13C與年值氣候因子多元回歸分析模型Table 2 Multiple regression analysis model between δ13C values and yearly climate factors

    4 討論

    不同種源杉木(1987—2013)樹輪α纖維素中穩(wěn)定碳同位素(δ13C)的變幅為-20.997‰~-25.910‰,且表現(xiàn)出較強的一階自相關(guān),相關(guān)系數(shù)為0.181~0.866,表明廣西柳州不同種源杉木樹輪δ13C的年際變化相對穩(wěn)定,能較好的呈現(xiàn)環(huán)境的變化,這與靳翔等[9]對岷江冷杉樹輪δ13C變幅及年際變化的研究結(jié)果相一致。

    本研究發(fā)現(xiàn),大部分杉木種源δ13C與年平均氣溫、最高氣溫、最低氣溫呈顯著正相關(guān),與最小相對濕度、平均相對濕度呈極顯著負相關(guān),而與年日照總時數(shù)呈不顯著正相關(guān)。杉木不同種源樹輪δ13C與氣候因子的這種正負相關(guān)性與對亞熱帶樹種猴面包樹[24]及亞高山暗針葉林[25]的研究相似。

    研究還發(fā)現(xiàn),不同種源的杉木對氣候因子響應(yīng)的靈敏度不同,可見,影響樹輪中δ13C值的因素并不單一,是多個因子的綜合效應(yīng)[26]。樹輪α纖維素中δ13C值與葉片氣孔導(dǎo)度及環(huán)境溫度相關(guān)[13],而葉片是通過感知環(huán)境中水分多少來調(diào)節(jié)氣孔導(dǎo)度,當葉片處于相對干燥的環(huán)境中,水分的減少導(dǎo)致氣孔導(dǎo)度降低,限制了葉片CO2的供給,最終導(dǎo)致樹輪中δ13C值增大[27]。溫度是影響CO2分餾的主要條件,當溫度低于平均溫度,且二者之間差值較大時,將會影響樹輪中碳的固定[28]。通過對杉木樹輪中δ13C值與當年、前1年、前2年氣候因子相關(guān)性的分析可以看出,17個種源對前2年平均氣溫、15個種源與前2年最低氣溫、21個種源與前2年最高氣溫以及11個種源與前2年的平均相對濕度相關(guān)系數(shù)較大,而14個種源與當年最小相對濕度呈顯著負相關(guān)性,除了廣東紫金、貴州六枝和湖南新寧的3個種源樹輪δ13C值對氣候因子的響應(yīng)是即時效應(yīng)外,其余的種源對氣候因子的響應(yīng)均存滯后效應(yīng)。由此可以看出,各氣候因子對杉木樹輪中穩(wěn)定碳同位素的影響是長期積累的結(jié)果,這與在低海拔地區(qū)樹輪穩(wěn)定碳同位素對年氣候因子的響應(yīng)出現(xiàn)明顯滯后效應(yīng)的結(jié)論一致[29]。

    通過對相關(guān)系數(shù)較高的氣候因子進行逐步回歸分析,發(fā)現(xiàn)除了福建南平(FJn)種源與前2年的最低氣溫和平均相對濕度呈線性關(guān)系,廣西博白(GXb)種源與前2年平均氣溫呈線性關(guān)系,云南西疇(Ync)、四川古藺(SCl)種源與前2年平均相對濕度呈線性關(guān)系外,其余26個種源杉木樹輪α纖維素中δ13C值隨當年最小相對濕度的增大而出現(xiàn)減小的趨勢。因此,當年最小相對濕度可以作為影響杉木樹輪α纖維素中δ13C的主導(dǎo)氣候因子。在濕度較大的環(huán)境中,杉木往往生長較好,其有機物的分餾作用越強,導(dǎo)致樹輪中δ13C減小[28]。

    5 結(jié)論

    中亞熱帶柳州地區(qū)不同地理種源杉木樹輪α-纖維素中δ13C對同一氣候因子的響應(yīng)方向基本一致,大部分種源樹輪δ13C與平均相對濕度和最小相對濕度呈負相關(guān),與平均氣溫、最高氣溫和最低氣溫呈顯著正相關(guān),與日照時數(shù)、降水和濕潤度指數(shù)的相關(guān)性不大,且有90%的杉木種源樹輪δ13C對氣候的響應(yīng)存在明顯滯后效應(yīng)。不同地理種源杉木樹輪α-纖維素中δ13C對最小相對濕度這一氣候因子的響應(yīng)較靈敏,從生物學(xué)基礎(chǔ)上闡明了樹輪穩(wěn)定碳同位素的固定、分餾受到生長環(huán)境中最小相對濕度的影響較大,而對年降水量多少的依賴性較小。研究結(jié)果有效的填補了種源層面上樹輪穩(wěn)定碳同位素對氣候響應(yīng)研究的不足,也為未來氣候變化條件下杉木人工林健康培育提供了依據(jù)。

    猜你喜歡
    樹輪氣候因子最低氣溫
    基于氣候變量的苧麻產(chǎn)量SSA-BP預(yù)測模型
    樹輪研究方法與進展
    天山南北坡樹輪穩(wěn)定碳同位素對氣候的響應(yīng)差異
    北辰地區(qū)日最低氣溫特征及影響因素研究
    天津科技(2019年3期)2019-03-30 07:17:14
    67年來濟南最高和最低氣溫變化特征
    樹木年輪對生態(tài)環(huán)境變化的響應(yīng)及應(yīng)用研究
    數(shù)學(xué)選擇題練習(xí)
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    阿爾泰山中部樹輪寬度年表特征及其氣候響應(yīng)分析
    久久亚洲真实| 精品人妻1区二区| 999精品在线视频| 欧美成人午夜精品| 老汉色∧v一级毛片| 日韩成人在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | av在线播放免费不卡| 草草在线视频免费看| 欧美黄色片欧美黄色片| 国产一区二区三区在线臀色熟女| 日日干狠狠操夜夜爽| 成人手机av| 18禁黄网站禁片午夜丰满| 午夜老司机福利片| 国产精品久久久久久人妻精品电影| 久久中文看片网| 久久久水蜜桃国产精品网| 国产高清视频在线观看网站| 欧美性猛交黑人性爽| 国产精品 国内视频| 亚洲五月婷婷丁香| 好男人在线观看高清免费视频| 久久性视频一级片| 桃色一区二区三区在线观看| 99久久国产精品久久久| 97超级碰碰碰精品色视频在线观看| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 亚洲午夜理论影院| 欧美性长视频在线观看| 亚洲熟女毛片儿| 久久久久久九九精品二区国产 | 99在线视频只有这里精品首页| 性欧美人与动物交配| 正在播放国产对白刺激| 色综合婷婷激情| 亚洲 欧美 日韩 在线 免费| 亚洲成人精品中文字幕电影| 天堂√8在线中文| 男女视频在线观看网站免费 | 久久国产乱子伦精品免费另类| 天堂√8在线中文| 大型av网站在线播放| 国产成人一区二区三区免费视频网站| 高潮久久久久久久久久久不卡| 国产亚洲精品av在线| 久久国产精品影院| 伊人久久大香线蕉亚洲五| 他把我摸到了高潮在线观看| 美女高潮喷水抽搐中文字幕| 美女高潮喷水抽搐中文字幕| 一夜夜www| 一进一出抽搐动态| 在线观看免费午夜福利视频| 亚洲免费av在线视频| 色综合欧美亚洲国产小说| 国产精品永久免费网站| 国产又黄又爽又无遮挡在线| 热99re8久久精品国产| 欧美成人免费av一区二区三区| 成人国语在线视频| 真人做人爱边吃奶动态| 又紧又爽又黄一区二区| 亚洲专区中文字幕在线| 亚洲成人中文字幕在线播放| 日韩三级视频一区二区三区| 精品国产乱码久久久久久男人| 国产黄色小视频在线观看| 国产亚洲欧美在线一区二区| 国产精品,欧美在线| 91字幕亚洲| 亚洲专区字幕在线| 日日爽夜夜爽网站| 可以免费在线观看a视频的电影网站| 亚洲五月婷婷丁香| 亚洲av成人av| 国产精品香港三级国产av潘金莲| 国产高清激情床上av| 亚洲熟妇熟女久久| 丝袜美腿诱惑在线| 很黄的视频免费| 国产探花在线观看一区二区| 夜夜看夜夜爽夜夜摸| 色噜噜av男人的天堂激情| 亚洲熟妇熟女久久| 免费人成视频x8x8入口观看| 午夜福利视频1000在线观看| 国产高清有码在线观看视频 | 好男人在线观看高清免费视频| 性欧美人与动物交配| 欧美极品一区二区三区四区| 国产av麻豆久久久久久久| 亚洲欧美精品综合一区二区三区| 日本黄色视频三级网站网址| 一级a爱片免费观看的视频| 亚洲av成人av| 欧美人与性动交α欧美精品济南到| 国产又黄又爽又无遮挡在线| 亚洲一区二区三区不卡视频| 午夜福利高清视频| 国产精品影院久久| 动漫黄色视频在线观看| 久久天堂一区二区三区四区| 免费人成视频x8x8入口观看| 十八禁人妻一区二区| 校园春色视频在线观看| 精品高清国产在线一区| 国产一区二区三区视频了| 日韩成人在线观看一区二区三区| 国产精品九九99| 成人手机av| 成人18禁高潮啪啪吃奶动态图| 在线十欧美十亚洲十日本专区| 午夜福利18| 久久久国产精品麻豆| 禁无遮挡网站| www.精华液| 一本综合久久免费| www.自偷自拍.com| 午夜免费观看网址| 久久久久久国产a免费观看| 久久性视频一级片| 欧美大码av| 男男h啪啪无遮挡| 精品电影一区二区在线| 19禁男女啪啪无遮挡网站| 亚洲欧美一区二区三区黑人| 制服诱惑二区| 黄色丝袜av网址大全| 国产1区2区3区精品| 精品国产亚洲在线| 国产又色又爽无遮挡免费看| 亚洲无线在线观看| 精品免费久久久久久久清纯| av片东京热男人的天堂| 成人三级做爰电影| 国产亚洲精品av在线| 欧美最黄视频在线播放免费| 国产精品免费视频内射| 午夜福利18| 国产精品,欧美在线| 九九热线精品视视频播放| 欧美在线一区亚洲| 十八禁人妻一区二区| 91麻豆av在线| 国内精品久久久久精免费| 欧美性长视频在线观看| 国产精品乱码一区二三区的特点| 亚洲精品一卡2卡三卡4卡5卡| www.自偷自拍.com| 1024香蕉在线观看| 在线观看美女被高潮喷水网站 | 全区人妻精品视频| 国产精品久久久久久久电影 | 国产精品日韩av在线免费观看| 欧美成人性av电影在线观看| 免费看十八禁软件| 操出白浆在线播放| 国产激情久久老熟女| 国产1区2区3区精品| 亚洲最大成人中文| 久久精品91蜜桃| 欧美又色又爽又黄视频| av在线播放免费不卡| 亚洲精品粉嫩美女一区| 欧美日韩乱码在线| 俺也久久电影网| 三级毛片av免费| 人妻丰满熟妇av一区二区三区| 亚洲 欧美 日韩 在线 免费| 亚洲成人中文字幕在线播放| avwww免费| 午夜免费观看网址| 亚洲一码二码三码区别大吗| 精品日产1卡2卡| 岛国在线观看网站| 日韩精品免费视频一区二区三区| 日本a在线网址| 精品一区二区三区四区五区乱码| 中文在线观看免费www的网站 | 黄片小视频在线播放| 男女视频在线观看网站免费 | 十八禁网站免费在线| 午夜免费激情av| 亚洲人成网站高清观看| 不卡av一区二区三区| 曰老女人黄片| 不卡一级毛片| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品999在线| 日韩欧美三级三区| av视频在线观看入口| 精品午夜福利视频在线观看一区| 啦啦啦观看免费观看视频高清| 在线十欧美十亚洲十日本专区| 国产99白浆流出| 国产一区二区在线观看日韩 | 999精品在线视频| 国产爱豆传媒在线观看 | 久久 成人 亚洲| 夜夜看夜夜爽夜夜摸| 国产欧美日韩精品亚洲av| 国产真实乱freesex| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 国产成人系列免费观看| 老汉色av国产亚洲站长工具| 国产精品久久久久久人妻精品电影| 日本三级黄在线观看| 久久精品91蜜桃| 在线免费观看的www视频| 好男人在线观看高清免费视频| 99热这里只有是精品50| 曰老女人黄片| 女人高潮潮喷娇喘18禁视频| 99国产综合亚洲精品| 桃色一区二区三区在线观看| 亚洲av中文字字幕乱码综合| 亚洲专区字幕在线| 亚洲国产精品999在线| av视频在线观看入口| 国产av一区在线观看免费| 女同久久另类99精品国产91| 亚洲av成人av| 国产成人av教育| 午夜日韩欧美国产| 亚洲乱码一区二区免费版| 舔av片在线| 两人在一起打扑克的视频| 无限看片的www在线观看| 波多野结衣高清无吗| 亚洲人成77777在线视频| 在线观看日韩欧美| 亚洲国产欧美网| 亚洲国产高清在线一区二区三| 午夜免费观看网址| 在线观看免费日韩欧美大片| 日韩欧美国产在线观看| 欧美成人午夜精品| 黄色女人牲交| 日韩高清综合在线| 国内精品久久久久精免费| 黄色 视频免费看| 床上黄色一级片| 1024手机看黄色片| 级片在线观看| 又黄又爽又免费观看的视频| 久久国产乱子伦精品免费另类| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 国产片内射在线| 欧美又色又爽又黄视频| 色精品久久人妻99蜜桃| 国产在线精品亚洲第一网站| 亚洲专区字幕在线| 久久国产精品人妻蜜桃| 亚洲午夜精品一区,二区,三区| 最好的美女福利视频网| 久久久国产精品麻豆| 午夜老司机福利片| 国产av麻豆久久久久久久| 一区二区三区高清视频在线| 色综合站精品国产| 久久久久久亚洲精品国产蜜桃av| netflix在线观看网站| 97人妻精品一区二区三区麻豆| 亚洲一码二码三码区别大吗| 黄色丝袜av网址大全| 国产99白浆流出| 老汉色av国产亚洲站长工具| 无人区码免费观看不卡| 黄色a级毛片大全视频| 欧美日韩黄片免| 长腿黑丝高跟| 欧美黑人巨大hd| 一卡2卡三卡四卡精品乱码亚洲| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 亚洲无线在线观看| 国产日本99.免费观看| 成人欧美大片| 老熟妇仑乱视频hdxx| 老汉色av国产亚洲站长工具| 制服诱惑二区| 国内久久婷婷六月综合欲色啪| 黑人巨大精品欧美一区二区mp4| 97人妻精品一区二区三区麻豆| 91av网站免费观看| 最新在线观看一区二区三区| 免费观看精品视频网站| www日本在线高清视频| 大型黄色视频在线免费观看| 国产精品一区二区三区四区久久| 国产av又大| 国产探花在线观看一区二区| 12—13女人毛片做爰片一| 亚洲激情在线av| 久久久久国内视频| 国产av一区在线观看免费| 黑人巨大精品欧美一区二区mp4| av在线播放免费不卡| 老鸭窝网址在线观看| 亚洲男人的天堂狠狠| 岛国在线观看网站| 亚洲乱码一区二区免费版| 视频区欧美日本亚洲| av福利片在线| 亚洲av中文字字幕乱码综合| 成人国产一区最新在线观看| 19禁男女啪啪无遮挡网站| 99久久国产精品久久久| 天堂av国产一区二区熟女人妻 | 色老头精品视频在线观看| 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 怎么达到女性高潮| 999久久久精品免费观看国产| 在线免费观看的www视频| 国产亚洲av嫩草精品影院| 两个人看的免费小视频| 12—13女人毛片做爰片一| 一个人免费在线观看的高清视频| 色av中文字幕| 欧美激情久久久久久爽电影| 精品久久久久久,| or卡值多少钱| 老司机午夜福利在线观看视频| 搞女人的毛片| 男人舔女人下体高潮全视频| 亚洲欧美一区二区三区黑人| www.熟女人妻精品国产| 亚洲人成网站在线播放欧美日韩| www日本黄色视频网| 国产高清有码在线观看视频 | 国产探花在线观看一区二区| 免费在线观看影片大全网站| 床上黄色一级片| 国内精品久久久久久久电影| 精品高清国产在线一区| 18禁观看日本| 国产三级黄色录像| 精品熟女少妇八av免费久了| 波多野结衣巨乳人妻| 人人妻,人人澡人人爽秒播| 叶爱在线成人免费视频播放| 最近最新免费中文字幕在线| 日本在线视频免费播放| 国语自产精品视频在线第100页| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 国产人伦9x9x在线观看| 亚洲电影在线观看av| 亚洲精品色激情综合| 国产成人一区二区三区免费视频网站| 亚洲av中文字字幕乱码综合| 国产精品影院久久| 国产成人精品无人区| 一本精品99久久精品77| 全区人妻精品视频| 国产精品一区二区三区四区久久| 好男人在线观看高清免费视频| 啦啦啦观看免费观看视频高清| 91字幕亚洲| 亚洲性夜色夜夜综合| 欧美成人免费av一区二区三区| 一边摸一边做爽爽视频免费| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 亚洲国产日韩欧美精品在线观看 | 国内精品久久久久久久电影| 搞女人的毛片| 亚洲激情在线av| 久久香蕉激情| 又黄又粗又硬又大视频| 日韩欧美一区二区三区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人操中国人逼视频| 久久精品91蜜桃| 99热只有精品国产| aaaaa片日本免费| 色老头精品视频在线观看| 久久天堂一区二区三区四区| 久99久视频精品免费| 欧美绝顶高潮抽搐喷水| 国产精品av视频在线免费观看| 桃红色精品国产亚洲av| 国产探花在线观看一区二区| 久久人人精品亚洲av| 一级片免费观看大全| 中文字幕久久专区| 搡老岳熟女国产| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 无遮挡黄片免费观看| 欧美黑人巨大hd| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 久久中文字幕人妻熟女| 怎么达到女性高潮| 亚洲国产精品sss在线观看| 久久久久久久精品吃奶| 九色国产91popny在线| 日韩欧美国产一区二区入口| 日日干狠狠操夜夜爽| 国产高清videossex| 老熟妇仑乱视频hdxx| 久久精品国产亚洲av高清一级| 伦理电影免费视频| 波多野结衣高清无吗| 欧美成人性av电影在线观看| 国产aⅴ精品一区二区三区波| 午夜福利高清视频| 中亚洲国语对白在线视频| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 亚洲精品久久国产高清桃花| 久久 成人 亚洲| av在线播放免费不卡| 1024手机看黄色片| 欧美乱码精品一区二区三区| 国产主播在线观看一区二区| 欧美色视频一区免费| 韩国av一区二区三区四区| 18禁国产床啪视频网站| 欧美色欧美亚洲另类二区| 757午夜福利合集在线观看| 国产精品,欧美在线| 在线看三级毛片| 亚洲 欧美 日韩 在线 免费| 一个人免费在线观看的高清视频| 精品久久久久久成人av| 99久久久亚洲精品蜜臀av| 精品国产乱码久久久久久男人| 亚洲第一电影网av| 啦啦啦免费观看视频1| 在线观看免费日韩欧美大片| 香蕉久久夜色| 国产欧美日韩一区二区三| 日本一区二区免费在线视频| 国产精品乱码一区二三区的特点| 亚洲av美国av| 97超级碰碰碰精品色视频在线观看| 国产精品98久久久久久宅男小说| 女生性感内裤真人,穿戴方法视频| 欧美中文日本在线观看视频| 悠悠久久av| 欧美zozozo另类| 亚洲av成人精品一区久久| 欧美av亚洲av综合av国产av| 中亚洲国语对白在线视频| 看片在线看免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 成人18禁在线播放| 无遮挡黄片免费观看| 久久久精品欧美日韩精品| 99在线人妻在线中文字幕| 蜜桃久久精品国产亚洲av| 国内少妇人妻偷人精品xxx网站 | 十八禁网站免费在线| 国产激情久久老熟女| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 欧美一级a爱片免费观看看 | 最新在线观看一区二区三区| 久久精品91蜜桃| 最近最新中文字幕大全免费视频| 国产激情偷乱视频一区二区| 国产免费av片在线观看野外av| 国产日本99.免费观看| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 亚洲av成人一区二区三| 国产91精品成人一区二区三区| 99久久精品国产亚洲精品| 一边摸一边做爽爽视频免费| 国产精品 国内视频| 国产亚洲精品一区二区www| 五月玫瑰六月丁香| 色av中文字幕| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 丁香六月欧美| 久久中文看片网| 国产成人啪精品午夜网站| 日本五十路高清| 欧美乱妇无乱码| 婷婷亚洲欧美| 91成年电影在线观看| 久久久久精品国产欧美久久久| 午夜福利欧美成人| 免费在线观看亚洲国产| 我的老师免费观看完整版| 久久性视频一级片| 国语自产精品视频在线第100页| 性色av乱码一区二区三区2| 一级毛片精品| 国产欧美日韩一区二区三| 亚洲欧美精品综合久久99| 日韩欧美精品v在线| 非洲黑人性xxxx精品又粗又长| 一区二区三区激情视频| 大型av网站在线播放| 亚洲专区字幕在线| 亚洲精品国产精品久久久不卡| √禁漫天堂资源中文www| 久久国产乱子伦精品免费另类| 日韩有码中文字幕| 亚洲成人免费电影在线观看| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 九九热线精品视视频播放| 精品一区二区三区av网在线观看| 美女大奶头视频| 禁无遮挡网站| 国产精品99久久99久久久不卡| 亚洲九九香蕉| 成人欧美大片| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 波多野结衣高清无吗| 午夜精品久久久久久毛片777| 变态另类成人亚洲欧美熟女| 亚洲九九香蕉| 看片在线看免费视频| videosex国产| 50天的宝宝边吃奶边哭怎么回事| 男女做爰动态图高潮gif福利片| 欧美中文日本在线观看视频| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 亚洲色图av天堂| 人成视频在线观看免费观看| 国产精品亚洲美女久久久| 最近最新中文字幕大全免费视频| 天天一区二区日本电影三级| 午夜亚洲福利在线播放| 久久精品亚洲精品国产色婷小说| 91成年电影在线观看| 脱女人内裤的视频| 日韩欧美三级三区| 伦理电影免费视频| 日韩精品免费视频一区二区三区| 欧美丝袜亚洲另类 | 好男人在线观看高清免费视频| 变态另类成人亚洲欧美熟女| 一个人观看的视频www高清免费观看 | 757午夜福利合集在线观看| 国产精品自产拍在线观看55亚洲| 日韩欧美在线乱码| 91成年电影在线观看| 亚洲一区二区三区不卡视频| 精品免费久久久久久久清纯| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲综合一区二区三区_| 在线观看舔阴道视频| 久久中文字幕一级| 久久精品亚洲精品国产色婷小说| 精品国产乱码久久久久久男人| 天天一区二区日本电影三级| 亚洲黑人精品在线| 精品少妇一区二区三区视频日本电影| 在线观看免费日韩欧美大片| 久久婷婷成人综合色麻豆| 美女 人体艺术 gogo| 久久国产乱子伦精品免费另类| 精品久久蜜臀av无| 日韩 欧美 亚洲 中文字幕| 国产成人啪精品午夜网站| 久久精品aⅴ一区二区三区四区| 最近最新中文字幕大全电影3| 色综合亚洲欧美另类图片| 亚洲欧洲精品一区二区精品久久久| 亚洲激情在线av| 亚洲在线自拍视频| 国产精品九九99| 日韩欧美在线乱码| 制服丝袜大香蕉在线| 午夜激情av网站| 亚洲av熟女| 国产单亲对白刺激| 日韩精品免费视频一区二区三区| 黄色丝袜av网址大全| 19禁男女啪啪无遮挡网站| 一本综合久久免费| 国产精品久久视频播放| 亚洲精品中文字幕在线视频| 久久久久性生活片| 亚洲成人中文字幕在线播放| 91大片在线观看| 久久精品综合一区二区三区| 亚洲自偷自拍图片 自拍| 免费搜索国产男女视频| 身体一侧抽搐| 久久久国产成人精品二区| 一进一出好大好爽视频| 久久 成人 亚洲| 性色av乱码一区二区三区2| 国产精品亚洲一级av第二区| 日韩有码中文字幕| 久久国产乱子伦精品免费另类| 国产精品 欧美亚洲| 19禁男女啪啪无遮挡网站| 99国产精品一区二区三区| 禁无遮挡网站| 欧美又色又爽又黄视频| 丰满人妻熟妇乱又伦精品不卡| 天堂√8在线中文| 欧美乱色亚洲激情| a在线观看视频网站| www.自偷自拍.com| 国产激情久久老熟女| 中文字幕人妻丝袜一区二区| 黑人操中国人逼视频|