• <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)分析
    亚洲人成网站在线播放欧美日韩| av天堂中文字幕网| 国产麻豆成人av免费视频| 亚洲av日韩精品久久久久久密| 日本与韩国留学比较| 99视频精品全部免费 在线| 天堂网av新在线| 成人二区视频| 国产精华一区二区三区| 夜夜爽天天搞| 床上黄色一级片| 在线观看一区二区三区| 全区人妻精品视频| 国产高清不卡午夜福利| 无人区码免费观看不卡| 亚洲va在线va天堂va国产| 亚州av有码| 日日夜夜操网爽| 免费观看人在逋| 小说图片视频综合网站| 欧美成人a在线观看| 女同久久另类99精品国产91| 日本 欧美在线| 麻豆成人av在线观看| 国产淫片久久久久久久久| 在线观看一区二区三区| 中文字幕av在线有码专区| 天美传媒精品一区二区| 啦啦啦韩国在线观看视频| 天天躁日日操中文字幕| 亚洲一级一片aⅴ在线观看| 无遮挡黄片免费观看| 在线天堂最新版资源| 亚洲中文日韩欧美视频| 国产精品久久电影中文字幕| 麻豆成人午夜福利视频| 99国产精品一区二区蜜桃av| 日日撸夜夜添| 午夜免费成人在线视频| 色哟哟哟哟哟哟| 又紧又爽又黄一区二区| 欧美高清成人免费视频www| 日韩欧美国产在线观看| 精品无人区乱码1区二区| 国产欧美日韩一区二区精品| 亚洲精品久久国产高清桃花| 欧美日韩中文字幕国产精品一区二区三区| 午夜福利在线观看免费完整高清在 | 精华霜和精华液先用哪个| 变态另类成人亚洲欧美熟女| 国产午夜精品论理片| 日韩一区二区视频免费看| 欧美另类亚洲清纯唯美| 51国产日韩欧美| 亚洲va在线va天堂va国产| av视频在线观看入口| 国产精品一区二区三区四区久久| 婷婷色综合大香蕉| 精品人妻偷拍中文字幕| 成人国产综合亚洲| 小蜜桃在线观看免费完整版高清| 内地一区二区视频在线| 亚洲最大成人手机在线| 免费看美女性在线毛片视频| 女生性感内裤真人,穿戴方法视频| 久久久午夜欧美精品| 麻豆成人av在线观看| 国产乱人伦免费视频| 免费人成在线观看视频色| 色视频www国产| 久久国内精品自在自线图片| 制服丝袜大香蕉在线| 嫩草影院精品99| 色吧在线观看| 日韩欧美精品v在线| 久久精品91蜜桃| 亚洲国产日韩欧美精品在线观看| 国产大屁股一区二区在线视频| 欧美色欧美亚洲另类二区| 一区二区三区免费毛片| 露出奶头的视频| 亚州av有码| 美女大奶头视频| 日韩欧美 国产精品| 动漫黄色视频在线观看| 国产黄片美女视频| 国产欧美日韩精品一区二区| 直男gayav资源| 一卡2卡三卡四卡精品乱码亚洲| 国内精品久久久久久久电影| 国产淫片久久久久久久久| 午夜福利在线观看吧| 免费搜索国产男女视频| 亚洲最大成人av| 国产高清视频在线观看网站| 免费av毛片视频| 日韩中字成人| 中文字幕精品亚洲无线码一区| 欧美激情在线99| 深爱激情五月婷婷| 床上黄色一级片| 亚洲一级一片aⅴ在线观看| 男人舔女人下体高潮全视频| 免费看a级黄色片| 国产爱豆传媒在线观看| 成人特级黄色片久久久久久久| 亚洲最大成人中文| 色综合婷婷激情| 亚洲人成网站高清观看| 最近最新中文字幕大全电影3| 蜜桃久久精品国产亚洲av| 在线看三级毛片| 国产欧美日韩一区二区精品| 日本一二三区视频观看| 午夜免费激情av| 97热精品久久久久久| 亚洲aⅴ乱码一区二区在线播放| 男女之事视频高清在线观看| 亚洲最大成人av| 精品久久久久久久久av| 亚洲一级一片aⅴ在线观看| 午夜a级毛片| 91av网一区二区| 日韩精品青青久久久久久| 精品一区二区三区视频在线观看免费| 国内精品一区二区在线观看| 亚洲中文日韩欧美视频| 免费一级毛片在线播放高清视频| 亚洲欧美精品综合久久99| 校园春色视频在线观看| 简卡轻食公司| 少妇人妻精品综合一区二区 | 中文字幕精品亚洲无线码一区| 欧美成人免费av一区二区三区| 亚洲av成人精品一区久久| 真实男女啪啪啪动态图| 99久久成人亚洲精品观看| 成人亚洲精品av一区二区| 国产精品伦人一区二区| xxxwww97欧美| 日本撒尿小便嘘嘘汇集6| bbb黄色大片| bbb黄色大片| 午夜激情欧美在线| 久9热在线精品视频| 好男人在线观看高清免费视频| 97热精品久久久久久| www.www免费av| 亚洲四区av| 国产视频内射| 亚洲av成人精品一区久久| 成人永久免费在线观看视频| 麻豆成人av在线观看| 久久国产乱子免费精品| 国产av不卡久久| 精品久久久久久久久av| 波多野结衣高清作品| 日韩欧美三级三区| 精品一区二区免费观看| avwww免费| 一级a爱片免费观看的视频| 亚洲成人精品中文字幕电影| 午夜免费成人在线视频| 成人国产一区最新在线观看| 国产av在哪里看| 国产 一区精品| 国产不卡一卡二| 91av网一区二区| 最新中文字幕久久久久| 亚洲av美国av| 日本免费a在线| 天天躁日日操中文字幕| 在线免费十八禁| 一卡2卡三卡四卡精品乱码亚洲| 色尼玛亚洲综合影院| 中文字幕久久专区| 少妇的逼水好多| 看十八女毛片水多多多| 日韩欧美 国产精品| 美女免费视频网站| 他把我摸到了高潮在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av第一区精品v没综合| 精品久久久久久久末码| 给我免费播放毛片高清在线观看| 免费看av在线观看网站| 亚洲美女搞黄在线观看 | 如何舔出高潮| 精品人妻一区二区三区麻豆 | 成人av在线播放网站| 99热这里只有是精品在线观看| 国产精品一及| 亚洲精品国产成人久久av| 天天躁日日操中文字幕| 成人无遮挡网站| 一区二区三区免费毛片| 午夜老司机福利剧场| www.www免费av| 男人的好看免费观看在线视频| 午夜福利在线观看吧| 精品久久久久久久末码| 日韩欧美精品免费久久| 亚洲七黄色美女视频| 校园人妻丝袜中文字幕| 婷婷亚洲欧美| 免费看光身美女| 成人特级黄色片久久久久久久| 国产精品国产三级国产av玫瑰| 中国美女看黄片| 高清在线国产一区| 白带黄色成豆腐渣| 久久人人精品亚洲av| 精品午夜福利在线看| 嫩草影院入口| 中文字幕精品亚洲无线码一区| 免费观看精品视频网站| 一进一出抽搐gif免费好疼| 欧美黑人欧美精品刺激| 丝袜美腿在线中文| 成人av在线播放网站| 国产 一区 欧美 日韩| 熟妇人妻久久中文字幕3abv| 成年女人看的毛片在线观看| 亚洲国产精品sss在线观看| 亚洲狠狠婷婷综合久久图片| 欧美xxxx性猛交bbbb| 直男gayav资源| 亚洲在线自拍视频| 少妇的逼水好多| 亚洲电影在线观看av| 亚洲成人中文字幕在线播放| 精品久久久久久久久av| 俺也久久电影网| 在线看三级毛片| 国产又黄又爽又无遮挡在线| 国产精品1区2区在线观看.| 一级a爱片免费观看的视频| 久久久久久大精品| 一进一出抽搐动态| 日日夜夜操网爽| 中文字幕av成人在线电影| 亚洲av成人av| 午夜免费激情av| 久久草成人影院| 国产欧美日韩精品一区二区| 亚洲最大成人中文| 久久久久久九九精品二区国产| 亚州av有码| 最后的刺客免费高清国语| 欧美另类亚洲清纯唯美| 欧美激情国产日韩精品一区| 免费av不卡在线播放| 国产高清视频在线播放一区| 少妇人妻一区二区三区视频| 热99在线观看视频| 色哟哟·www| 男女之事视频高清在线观看| 中国美女看黄片| 少妇人妻精品综合一区二区 | 国产精品一区www在线观看 | 91久久精品电影网| 深夜a级毛片| 美女黄网站色视频| av在线天堂中文字幕| 日韩欧美精品免费久久| 最后的刺客免费高清国语| 熟女人妻精品中文字幕| 亚洲av不卡在线观看| 十八禁国产超污无遮挡网站| a级一级毛片免费在线观看| 少妇熟女aⅴ在线视频| 亚洲一区二区三区色噜噜| 久久人人爽人人爽人人片va| 1024手机看黄色片| 九九爱精品视频在线观看| 日本-黄色视频高清免费观看| 国产真实伦视频高清在线观看 | 国产淫片久久久久久久久| 婷婷丁香在线五月| 日韩大尺度精品在线看网址| 欧美激情国产日韩精品一区| 亚洲成人免费电影在线观看| 男人舔奶头视频| 免费av毛片视频| 中文字幕熟女人妻在线| 亚洲精品久久国产高清桃花| 国产麻豆成人av免费视频| 99热这里只有精品一区| 免费一级毛片在线播放高清视频| 国产毛片a区久久久久| 国内少妇人妻偷人精品xxx网站| 少妇丰满av| 成人特级av手机在线观看| 一个人观看的视频www高清免费观看| 成人av在线播放网站| 欧洲精品卡2卡3卡4卡5卡区| 看黄色毛片网站| 国产毛片a区久久久久| 日韩中文字幕欧美一区二区| 亚洲欧美日韩无卡精品| 欧美黑人欧美精品刺激| 亚洲aⅴ乱码一区二区在线播放| bbb黄色大片| 老女人水多毛片| 国产探花极品一区二区| 91久久精品电影网| 婷婷六月久久综合丁香| 极品教师在线免费播放| 一本一本综合久久| 亚洲avbb在线观看| 国产激情偷乱视频一区二区| 国语自产精品视频在线第100页| 婷婷亚洲欧美| 国产极品精品免费视频能看的| 精品久久久久久久久亚洲 | 精品日产1卡2卡| 小蜜桃在线观看免费完整版高清| 国产熟女欧美一区二区| 97人妻精品一区二区三区麻豆| 国产黄a三级三级三级人| 中国美女看黄片| 91麻豆av在线| 亚洲午夜理论影院| 国产一区二区三区视频了| 国产高清激情床上av| 欧美一区二区亚洲| 亚洲成人久久爱视频| 乱码一卡2卡4卡精品| 国产毛片a区久久久久| av天堂在线播放| 亚洲内射少妇av| 亚洲欧美清纯卡通| 蜜桃久久精品国产亚洲av| 亚洲无线观看免费| 12—13女人毛片做爰片一| 香蕉av资源在线| 一区二区三区四区激情视频 | 国产老妇女一区| eeuss影院久久| 超碰av人人做人人爽久久| bbb黄色大片| 日韩大尺度精品在线看网址| 国产极品精品免费视频能看的| 国产一区二区三区在线臀色熟女| 蜜桃亚洲精品一区二区三区| 久久香蕉精品热| www日本黄色视频网| 国产精品精品国产色婷婷| 在线免费十八禁| 精品一区二区三区av网在线观看| .国产精品久久| 久久精品久久久久久噜噜老黄 | 国产精品,欧美在线| 色播亚洲综合网| 精品午夜福利视频在线观看一区| 尾随美女入室| 久久久午夜欧美精品| 麻豆国产97在线/欧美| 欧美色视频一区免费| 99热这里只有是精品50| 麻豆成人av在线观看| 国产精品av视频在线免费观看| 一区二区三区四区激情视频 | 亚洲一区二区三区色噜噜| 美女黄网站色视频| a在线观看视频网站| 在线播放国产精品三级| 看十八女毛片水多多多| 97人妻精品一区二区三区麻豆| 69人妻影院| 欧美日韩亚洲国产一区二区在线观看| 国产精品99久久久久久久久| 亚洲真实伦在线观看| 嫩草影院精品99| 国产麻豆成人av免费视频| 日韩一区二区视频免费看| 99久久久亚洲精品蜜臀av| 亚洲熟妇中文字幕五十中出| 嫁个100分男人电影在线观看| 高清在线国产一区| 国产精品98久久久久久宅男小说| 婷婷色综合大香蕉| 日本a在线网址| 午夜福利成人在线免费观看| 天天一区二区日本电影三级| 亚洲成人中文字幕在线播放| 狠狠狠狠99中文字幕| 99久国产av精品| 日韩,欧美,国产一区二区三区 | 久久欧美精品欧美久久欧美| 久久久久久大精品| 久久久久性生活片| 亚洲国产日韩欧美精品在线观看| 亚洲精品一区av在线观看| 美女xxoo啪啪120秒动态图| 日本黄色片子视频| 老司机深夜福利视频在线观看| 深夜a级毛片| 久久久久久久精品吃奶| 九九爱精品视频在线观看| 亚洲人成网站在线播放欧美日韩| 深夜精品福利| 人妻少妇偷人精品九色| 国产亚洲精品综合一区在线观看| 人人妻人人澡欧美一区二区| 91在线精品国自产拍蜜月| 亚洲欧美精品综合久久99| 色综合色国产| 99热网站在线观看| 毛片女人毛片| 精品国内亚洲2022精品成人| 亚洲精品一区av在线观看| 嫩草影院新地址| 国产久久久一区二区三区| 欧美区成人在线视频| 久久久久久久久久黄片| 色综合婷婷激情| 亚洲真实伦在线观看| 狠狠狠狠99中文字幕| 亚洲美女搞黄在线观看 | 日韩在线高清观看一区二区三区 | 波多野结衣高清作品| 午夜爱爱视频在线播放| 国产亚洲精品综合一区在线观看| 欧美日本亚洲视频在线播放| 亚洲七黄色美女视频| 亚洲人成网站高清观看| 97人妻精品一区二区三区麻豆| 免费观看的影片在线观看| 欧美人与善性xxx| 在线观看一区二区三区| 久久久色成人| netflix在线观看网站| 最新中文字幕久久久久| 国产淫片久久久久久久久| 国产成人一区二区在线| 他把我摸到了高潮在线观看| 三级毛片av免费| 欧美区成人在线视频| 亚洲av日韩精品久久久久久密| 欧美日本视频| 亚洲午夜理论影院| 欧美日韩精品成人综合77777| 97热精品久久久久久| 赤兔流量卡办理| 18禁黄网站禁片免费观看直播| 欧美xxxx黑人xx丫x性爽| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区 | 国产单亲对白刺激| 在线a可以看的网站| 日韩欧美免费精品| 九色国产91popny在线| 美女大奶头视频| 精品欧美国产一区二区三| 日韩一区二区视频免费看| 人人妻人人澡欧美一区二区| 一个人观看的视频www高清免费观看| 免费观看精品视频网站| 无遮挡黄片免费观看| 嫁个100分男人电影在线观看| 少妇的逼水好多| 久久亚洲精品不卡| 国产淫片久久久久久久久| 日韩,欧美,国产一区二区三区 | www日本黄色视频网| 美女cb高潮喷水在线观看| 搡老熟女国产l中国老女人| 一进一出好大好爽视频| 亚洲成人久久爱视频| 亚洲最大成人中文| 国产国拍精品亚洲av在线观看| 成人性生交大片免费视频hd| 伦理电影大哥的女人| 国产av不卡久久| 亚洲av不卡在线观看| av在线蜜桃| 一区二区三区四区激情视频 | 在线观看一区二区三区| 国产精品美女特级片免费视频播放器| 国产探花极品一区二区| 日韩欧美三级三区| 久久久国产成人精品二区| 中文亚洲av片在线观看爽| 亚洲狠狠婷婷综合久久图片| 色尼玛亚洲综合影院| 日韩一本色道免费dvd| 色播亚洲综合网| 国产亚洲精品综合一区在线观看| 黄色女人牲交| 国产亚洲欧美98| 动漫黄色视频在线观看| 天美传媒精品一区二区| 精品久久久噜噜| 国产一区二区在线观看日韩| 亚洲国产色片| 免费无遮挡裸体视频| 成年女人永久免费观看视频| 国产精品久久电影中文字幕| 国产爱豆传媒在线观看| 久久午夜福利片| 又黄又爽又免费观看的视频| 国产成人福利小说| 内地一区二区视频在线| 日韩人妻高清精品专区| 最近最新免费中文字幕在线| 亚洲中文日韩欧美视频| 一进一出好大好爽视频| 国产亚洲欧美98| 日本黄色视频三级网站网址| 观看美女的网站| 国产精品野战在线观看| 日日撸夜夜添| 91精品国产九色| ponron亚洲| 国产成人影院久久av| 亚洲最大成人中文| 少妇的逼好多水| 内地一区二区视频在线| 亚洲精品一区av在线观看| 久久国内精品自在自线图片| 国产毛片a区久久久久| 国产黄色小视频在线观看| 免费人成在线观看视频色| 美女黄网站色视频| 人妻久久中文字幕网| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看| 欧美人与善性xxx| 国产一区二区在线av高清观看| 日韩欧美在线二视频| 国产真实伦视频高清在线观看 | 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 久久国产精品人妻蜜桃| 欧美bdsm另类| 两人在一起打扑克的视频| 国产高潮美女av| 国产亚洲av嫩草精品影院| 亚洲av二区三区四区| 搡老妇女老女人老熟妇| 国产亚洲精品综合一区在线观看| 免费高清视频大片| 天天躁日日操中文字幕| 精品99又大又爽又粗少妇毛片 | 国产淫片久久久久久久久| 欧美高清性xxxxhd video| 亚洲国产精品sss在线观看| 桃色一区二区三区在线观看| 亚洲精品亚洲一区二区| 精品人妻熟女av久视频| 两个人的视频大全免费| 舔av片在线| 又爽又黄a免费视频| 久久亚洲精品不卡| 99久国产av精品| 亚洲美女黄片视频| 51国产日韩欧美| 成人国产综合亚洲| 久久人人爽人人爽人人片va| 久久国产精品人妻蜜桃| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月| 日本熟妇午夜| 免费看美女性在线毛片视频| 日日摸夜夜添夜夜添av毛片 | 亚洲图色成人| 欧美日韩精品成人综合77777| 久久久国产成人精品二区| 男人舔奶头视频| 午夜精品久久久久久毛片777| av.在线天堂| 18禁裸乳无遮挡免费网站照片| 黄色配什么色好看| 日本免费a在线| 国产成人一区二区在线| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 一卡2卡三卡四卡精品乱码亚洲| 国产一区二区在线观看日韩| 成熟少妇高潮喷水视频| 久久精品久久久久久噜噜老黄 | 人妻久久中文字幕网| 无遮挡黄片免费观看| 两人在一起打扑克的视频| 又粗又爽又猛毛片免费看| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 男人的好看免费观看在线视频| 在线观看66精品国产| 精品久久久久久久末码| 国语自产精品视频在线第100页| 91久久精品国产一区二区成人| 又黄又爽又免费观看的视频| а√天堂www在线а√下载| 国产精品99久久久久久久久| 亚洲av一区综合| 亚洲va日本ⅴa欧美va伊人久久| 国产伦人伦偷精品视频| 伊人久久精品亚洲午夜| 日日干狠狠操夜夜爽| 久久久久国产精品人妻aⅴ院| 国产精品久久视频播放| 国产极品精品免费视频能看的| 能在线免费观看的黄片| 日本五十路高清| 国产v大片淫在线免费观看| 在线观看一区二区三区| 久久国产精品人妻蜜桃| АⅤ资源中文在线天堂| 看片在线看免费视频| 日本免费a在线| 亚洲经典国产精华液单| 春色校园在线视频观看|