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

    模擬氮沉降下杜仲幼苗的光譜響應(yīng)

    2014-07-02 01:45:41哀建國杜江麗金松恒宋新章
    浙江林業(yè)科技 2014年6期
    關(guān)鍵詞:中氮杜仲植被指數(shù)

    哀建國,杜江麗,金松恒,宋新章

    (1.浙江農(nóng)林大學(xué) 天目學(xué)院,浙江 臨安 311300;2.浙江農(nóng)林大學(xué) 林業(yè)與生物技術(shù)學(xué)院,浙江 臨安 311300)

    模擬氮沉降下杜仲幼苗的光譜響應(yīng)

    哀建國1,杜江麗2,金松恒1,宋新章2

    (1.浙江農(nóng)林大學(xué) 天目學(xué)院,浙江 臨安 311300;2.浙江農(nóng)林大學(xué) 林業(yè)與生物技術(shù)學(xué)院,浙江 臨安 311300)

    設(shè)置對(duì)照(0 kg·hm-2·a-1)、低氮(30 kg·hm-2·a-1)、中氮(60 kg·hm-2·a-1)3組模擬大棚試驗(yàn),研究了杜仲(Eucommia ulmoides)幼苗在不同氮沉降處理下的光譜響應(yīng)特征。結(jié)果表明:杜仲幼苗在低氮處理下長勢(shì)良好,其光譜反射率總體處于最低水平,紅邊位置出現(xiàn)“紅移”趨勢(shì),中氮使得杜仲葉片在可見光及近紅外波段的反射率普遍升高,并發(fā)生“藍(lán)移”現(xiàn)象,杜仲生長受到抑制,后期抑制則逐漸消除;氮沉降處理間的差異顯著性光譜區(qū)域隨著時(shí)間的推移逐漸增多,到10月植物生長后期差異顯著性開始降低;反高斯模型參數(shù)Rs、R0、λ0、λp和σ可作為反映杜仲受氮沉降影響整體水平的指標(biāo),其中 R0最為顯著;與兩個(gè)歸一化植被指數(shù) NDVI和gNDVI相比,比值植被指數(shù)RVI700和RVI550表現(xiàn)出更為有效的指示脅迫作用。

    杜仲;模擬氮沉降;光譜反射率;紅邊

    近幾十年來,農(nóng)業(yè)和工業(yè)活動(dòng)的密集化,例如礦物燃料的燃燒、含氮化肥的生產(chǎn)和使用及畜牧業(yè)等人類活動(dòng)向大氣中排放了大量的含氮化合物[1]。目前中國南方已經(jīng)出現(xiàn)了大面積的氮沉降危害區(qū),而我國也已成為繼歐、美之后的第三大氮沉降集中區(qū)[1~2]。現(xiàn)有研究表明,氮沉降會(huì)對(duì)植物生產(chǎn)力[3~4]、光合生理[5]、營養(yǎng)元素含量[6]等產(chǎn)生影響,但以上研究多以農(nóng)作物、經(jīng)濟(jì)作物和部分森林樹種或整個(gè)森林生態(tài)系統(tǒng)為對(duì)象,對(duì)藥用植物特別是木本藥用植物則研究較少。

    杜仲(Eucommia ulmoides)為杜仲科(Eucommiaceae)杜仲屬(Eucommia)多年生落葉喬木,為我國特有珍貴中藥材和工業(yè)提膠原料樹種。杜仲在國內(nèi)栽培區(qū)域大致位于華東、華中和西南暖溫帶氣候區(qū)內(nèi)[7],其分布區(qū)大體上和長江流域相吻合,而此區(qū)域亦是我國氮沉降分布最集中地區(qū)。

    目前,高光譜遙感由于具有很高的光譜分辨率,能夠提供更為豐富的地物信息,已被廣泛用于監(jiān)測(cè)和研究環(huán)境變化引起的植物功能的變化。當(dāng)植物處于病蟲害、干旱、營養(yǎng)不良以及各類大氣污染物,諸如氮沉降、臭氧、重金屬等環(huán)境脅迫時(shí),植物葉片的一些生理特性會(huì)隨環(huán)境因子的脅迫而發(fā)生改變,進(jìn)而將引起葉片光譜曲線形態(tài)的變化,因此可以利用遙感數(shù)據(jù)(光譜反射率)來獲取植物對(duì)環(huán)境脅迫的響應(yīng)信息[8~9]。本文利用高光譜遙感技術(shù)研究模擬氮沉降對(duì)杜仲幼苗葉片反射光譜的影響,以期為實(shí)現(xiàn)浙江高氮沉降區(qū)杜仲的光譜實(shí)時(shí)監(jiān)測(cè)提供必要的基礎(chǔ)數(shù)據(jù)和技術(shù)支持。

    1 試驗(yàn)地概況

    試驗(yàn)地位于浙江省臨安市浙江農(nóng)林大學(xué)天目學(xué)院實(shí)驗(yàn)大棚,該地區(qū)年平均氣溫 15.8℃,7月為最熱月,歷年平均為28.1℃。年降水量1 426 mm,6月雨量最多,12月最少,屬中緯度北亞熱帶季風(fēng)氣候。試驗(yàn)地處于近自然狀態(tài),沒有人為經(jīng)營干擾。

    2 研究方法

    2.1 試驗(yàn)設(shè)計(jì)

    2010年4月初選取81株生長健壯、高度相近的一年生杜仲實(shí)生苗,移栽到高20 cm,內(nèi)徑18 cm的塑料花盆中,栽培土壤由黃紅壤和肥料土按7:1的比例混合而成,并置于實(shí)驗(yàn)大棚內(nèi)培養(yǎng)。實(shí)驗(yàn)共設(shè)9個(gè)樣方,每個(gè)樣方面積為1 m×1m,均勻放置9盆杜仲幼苗,樣方間距為0.5 m。實(shí)驗(yàn)分為對(duì)照(ck:0 kg·hm-2·a-1)、低氮(LN:30 kg·hm-2·a-1)、中氮(MN:60 kg·hm-2·a-1)3個(gè)處理組,9株為一重復(fù),每個(gè)處理分別有3個(gè)重復(fù)。當(dāng)受試杜仲苗過了緩苗期后,分別將每個(gè)處理所需的NH4NO3溶解在50 mL水中(相當(dāng)于浙江省年新增降水0.6 mm),自2013年5月開始,在杜仲苗生長發(fā)育期內(nèi)(5-10月),每月一次定期用噴霧器均勻噴灑受試杜仲苗。對(duì)照組則噴灑等量自來水。之后對(duì)苗木進(jìn)行常規(guī)栽培管理。氮沉降處理后每隔25 d測(cè)定反射光譜。

    2.2 反射光譜數(shù)據(jù)的采集

    葉片光譜反射率采用美國Unispec-SC型單通道光纖光譜儀(波段范圍:310 ~ 1 100 nm;光譜分辨率:3 nm;采樣間隔:1 nm)進(jìn)行采集。光譜反射率的采集選擇在晴朗無云或少云的天氣進(jìn)行,采集時(shí)間為10:00-14:00(太陽高度角大于45°)。實(shí)驗(yàn)時(shí)從每個(gè)樣方中隨機(jī)選取3株杜仲苗,在每株的中上部隨機(jī)選取光照充分、葉色均勻無病斑的成熟葉片6片,每葉1個(gè)數(shù)據(jù)采集點(diǎn)。每次測(cè)量前需用標(biāo)準(zhǔn)白板校對(duì)調(diào)零。采集得的葉片光譜反射率數(shù)據(jù)先用Multispec5.1導(dǎo)出,利用Excel進(jìn)行初步處理并運(yùn)算,最后再用Origin 8.0進(jìn)行整理分析并作圖。

    2.3 數(shù)據(jù)處理

    2.3.1 光譜反射率 對(duì)18組原始光譜反射率數(shù)據(jù)求平均值得到平均光譜反射率,并對(duì)受試杜仲苗在低氮和中氮處理下400 ~ 800nm的平均光譜反射率做差,進(jìn)行差異顯著性分析。

    2.3.2 植被指數(shù) 本實(shí)驗(yàn)中,選取了4個(gè)常用的植被指數(shù)(表 1)分析其在不同氮沉降水平的模擬實(shí)驗(yàn)作用下對(duì)受試杜仲樣本葉的區(qū)分度。其中,歸一化植被指數(shù)(NDVI)是應(yīng)用最為廣泛的一種植被指數(shù),常被用于植物脅迫的檢驗(yàn)[10]。在后續(xù)計(jì)算中,一律采用760 ~ 900 nm和630 ~ 690 nm所在波段范圍內(nèi)光譜反射率的平均值來表示NDVI 算式中的近紅外和紅光波段的值。

    表1 植被指數(shù)Table 1 Vegetation index

    2.3.3 反高斯模型 反高斯模型(inverted-Gaussian,IG)由Bonham-Carter于1988年首次提取,能夠模擬植被在670 ~ 780 nm范圍內(nèi)的光譜反射率曲線,用于定量分析植被的紅邊特性,其表達(dá)式為:

    式中,R0為紅光區(qū)光譜反射率最小值,Rs為近紅外區(qū)域最大光譜反射率值;λp為紅邊波長位置(REP),λ0為紅光區(qū)光譜反射率最小值所對(duì)應(yīng)的波長;σ為高斯方程偏差參數(shù),單位為nm,也等于λp-λ0。

    本實(shí)驗(yàn)中,R0和Rs分別選定為波長670 ~ 685 nm和780 ~ 795 nm范圍的光譜反射率平均值,并利用最小二乘法來估算685 ~ 780 nm波段內(nèi)另外兩個(gè)參數(shù)λ0和σ[14]:

    利用參數(shù)R0和Rs在B(λ)和λ間通過最小二乘法進(jìn)行線性擬合,并對(duì)紅邊參數(shù)模型的未知系數(shù)α0和α1進(jìn)行估算,則IG模型的參數(shù)λ0和σ分別為:

    3 結(jié)果與分析

    3.1 氮沉降對(duì)杜仲葉片光譜反射率的影響

    圖1 模擬氮沉降處理對(duì)杜仲平均光譜反射率的影響Figure 1 Mean spectral reflectance of seedlings under different nitrogen deposition

    反射率曲線是宏觀上觀察植物受脅迫狀態(tài)的最佳指標(biāo),反射率越高,則表明植物吸收的光越少,受脅迫越大。圖1為受試杜仲苗在5個(gè)月實(shí)驗(yàn)觀測(cè)中不同氮沉降處理下的平均光譜反射率比較。實(shí)驗(yàn)組和對(duì)照組葉片的原始反射光譜變化趨勢(shì)是一致的,都具有450 nm和650 nm處兩個(gè)明顯的吸收谷、550 nm處綠光反射峰、680 ~740 nm反射率急劇增高的紅邊區(qū)以及780 ~ 1 100 nm的近紅外高原平臺(tái),但在不同處理間差異顯著。實(shí)驗(yàn)初期和中期,在可見光范圍內(nèi)(400 ~ 700 nm),中氮組杜仲葉片平均光譜反射率最高,低氮組最低,說明中氮處理對(duì)受試杜仲苗造成了一定的脅迫,而低氮處理則起到了“施肥效應(yīng)”。后期,低氮和中氮實(shí)驗(yàn)組兩條光譜反射曲線幾乎重合,其平均光譜反射率都低于對(duì)照組,這種變化可能是杜仲苗對(duì)中濃度的氮沉降脅迫產(chǎn)生了一定耐受性的緣故。

    3.2 低氮和中氮處理的差異顯著性分析

    從圖1可知,中氮(MN)組的杜仲葉片光譜反射率在實(shí)驗(yàn)前期和中期都要高于低氮(NM)處理。為了進(jìn)一步說明受試杜仲苗對(duì)兩種濃度模擬氮沉降的響應(yīng)程度,計(jì)算二者平均光譜反射率間的差值(LN-MN),并通過Excel分析其差異大小,用Origin 8.0作圖比較,結(jié)果如圖2。

    圖2 低氮和中氮處理下杜仲葉片的平均光譜反射率差值(LN-MN)Figure 2 Difference of mean spectral reflectance of seedlings between LM and MN treatment

    由圖2可以看出,在550 nm為中心的綠光反射區(qū)、700 nm附近的紅邊以及近紅外部分,低氮和中氮處理的平均光譜反射率間差異比較顯著。此外,5-9月,在可見光部分低氮處理的光譜反射率幾乎都低于中氮處理,并且隨著時(shí)間的推移差值逐漸增大,到2013年9月差值達(dá)到最大。這種差值的差異顯著性(黑色部分)與差值大小變化基本一致,即差異顯著性的光譜區(qū)間隨著時(shí)間的推移也在逐漸增多,到10月份植物生長后期開始降低。上述結(jié)果表明,在杜仲生長的初期和中期,中氮處理使杜仲生長受到了一定程度的脅迫,而低氮處理長勢(shì)良好,差異顯著;后期,兩種處理間的差異不顯著。

    3.3 模擬氮沉降處理對(duì)杜仲植被指數(shù)的比較

    對(duì)整個(gè)實(shí)驗(yàn)期間杜仲的光譜反射率平均值進(jìn)行植被指數(shù)計(jì)算,結(jié)果如表2。由表2可知,兩種濃度模擬氮沉降處理下的杜仲植被指數(shù)變化特征基本一致,各個(gè)指數(shù)LN處理的值基本上都高于MN和ck的值,表明在低氮處理下受試杜仲苗生長狀態(tài)更佳,中氮處理則在一定程度上抑制了杜仲幼苗的生長。進(jìn)一步分析比較發(fā)現(xiàn),實(shí)驗(yàn)組和對(duì)照組的RVI700和RVI550差值比gNDVI和NDVI的差值更大,說明利用特定波長位置上的光譜反射率所構(gòu)建的植被指數(shù),對(duì)有效區(qū)分不同氮沉降處理水平下的杜仲樣本更加有效。

    3.4 模擬氮沉降對(duì)杜仲反高斯模型紅邊參數(shù)的影響

    近年來隨著高光譜遙感的興起而發(fā)展起來的光譜數(shù)據(jù)分析技術(shù)中,植被“紅邊”位移現(xiàn)象是研究最多、成效最顯著的成果之一?!凹t邊”是綠色植物在紅光范圍(680 ~ 750 nm)之間反射率增高最快的點(diǎn),即紅光范圍內(nèi)一階導(dǎo)數(shù)最大值所對(duì)應(yīng)的光譜位置。

    表2 不同處理時(shí)間的杜仲植被指數(shù)Table 2 Vegetation index of seedlings at different time

    本實(shí)驗(yàn)中,利用反高斯(IG) 模型來模擬 杜仲在氮沉降處理下的光譜反射 率的紅邊位置曲線,并對(duì)模型中的參數(shù)進(jìn)行定量分析,進(jìn)而篩選出相對(duì)明顯的參數(shù)作為反應(yīng)氮沉降脅迫的指標(biāo)。IG模型中共有5個(gè)參數(shù),分別為Rs、R0、λ0、λp和σ0,經(jīng)計(jì)算,5個(gè)參數(shù)同處理水平間具有較好的相關(guān)性,結(jié)果如圖3。

    Rs是近紅外區(qū)域肩反射率值,結(jié)果顯示,低氮組Rs值總體上最大(7月和10月除外),表明低氮處理會(huì)降低杜仲在近紅外區(qū)域肩吸收強(qiáng)度,中氮和對(duì)照組則無明顯規(guī)律,不同處理的Rs值總體上差異不明顯。R0值位于葉綠素a的吸收峰值區(qū),其值受葉綠素含量的影響。前期和中期,3組氮沉降處理的R0值呈LN > ck > MN,后期則對(duì)照組最大。說明在杜仲生長季低濃度的氮沉降有助于葉片中葉綠素的形成,這也能解釋實(shí)驗(yàn)中低氮處理下的杜仲葉片在多數(shù)情況下比對(duì)照組高的現(xiàn)象。

    圖3 氮沉降處理下反高斯模型參數(shù)比較Figure 3 Comparison of IG model parameters under different nitrogen deposition treatments(mean±SE)

    λ0和λp分別指的是紅光區(qū)域葉綠素吸收谷和紅邊所對(duì)應(yīng)的波長位置,二者間的關(guān)系為:λp= λ0+σ0中期,低氮組的λ0值相對(duì)較大,中氮組和對(duì)照組差異不顯著,初期和后期則無明顯規(guī)律。5-8月,低氮組的λp值要高于對(duì)照和中氮組,中氮最低,表明低氮處理會(huì)使杜仲產(chǎn)生“紅移”,促進(jìn)了杜仲生長,而中氮處理在這一時(shí)期則產(chǎn)生了一定的抑制作用。此外,6-8月,中氮組杜仲紅邊位置發(fā)生了“藍(lán)移”,后期又逐漸向長波方向移動(dòng),推測(cè)可能因?yàn)樵缙谥袧舛鹊牡两狄种屏硕胖偕L,但在時(shí)間足夠長的情況下,杜仲會(huì)對(duì)逆境脅迫表現(xiàn)出一定的適應(yīng)性和承受力。隨著杜仲葉片長大,葉面積指數(shù)增大,同時(shí)葉綠素含量增多,中濃度的氮沉降所造成的脅迫在一定程度上就會(huì)被減弱甚至消除。σ是反高斯模型的標(biāo)準(zhǔn)差系數(shù),對(duì)應(yīng)紅光吸收谷的寬度,其值受λ0和λp影響,無明顯規(guī)律。

    4 結(jié)論與討論

    (1)研究結(jié)果表明,與對(duì)照組相比,杜仲幼苗在低氮處理下長勢(shì)良好,其光譜反射率總體處于最低水平,表明低濃度的氮沉降促進(jìn)了杜仲生長,在整個(gè)實(shí)驗(yàn)期間都起到了一定的“施肥效應(yīng)”;中濃度的氮沉降使得杜仲葉片在可見光及近紅外波段的反射率普遍升高,杜仲生長受到抑制,這種抑制作用在后期逐漸消除,且與低氮組的光譜反射曲線幾乎重合。氮沉降處理間的差異顯著,但到10月植物生長后期開始降低。

    (2)與NDVI和gNDVI相比,植被指數(shù)RVI700和RVI550表現(xiàn)出更為有效的指示脅迫作用。本次研究中考慮到利用特定波長位置上的光譜反射率所構(gòu)建的植被指數(shù),對(duì)有效區(qū)分不同氮沉降處理水平下的杜仲樣本比較寬波長范圍的傳統(tǒng)植被指數(shù)NDVI更加有效,因此,可以考慮用RVI700等特定波長位置上有顯著指示作用的指數(shù)代替NDVI指數(shù)作為杜仲受氮沉降脅迫的指標(biāo)。

    (3)反高斯模型參數(shù)Rs、R0、λ0、λp和σ在6個(gè)時(shí)期均可作為反映杜仲受氮沉降脅迫整體水平的指標(biāo),其中參數(shù) R0最為顯著。在杜仲生長季,低氮組杜仲的R0最低,表明在該時(shí)期低濃度的氮沉降有助于葉片中葉綠素的形成,這也能解釋實(shí)驗(yàn)中低氮處理下的杜仲葉片在多數(shù)情況下比對(duì)照組高的現(xiàn)象。

    (4)在實(shí)驗(yàn)初期和中期,與對(duì)照組相比,低氮組杜仲光譜反射率紅邊位置總體上呈“紅移”趨勢(shì),中氮組則發(fā)生了“藍(lán)移”,并且表現(xiàn)出時(shí)間上的累積性,反映出低氮處理對(duì)杜仲生長的促進(jìn)和中氮處理的脅迫。但在后期,中氮組也發(fā)生了“紅移”現(xiàn)象,分析可能是因?yàn)槎胖僭谏L后期表現(xiàn)出了一定的適應(yīng)性。

    [1] Galloway J N, Aber J D, Erisman J W, et al.The nitrogen cascade[J].Bioscience, 2003, 53(4): 341-356.

    [2] Townsend AR, Braswell B H, Holland E A, et al.Spatial and temporal patterns in terrestrial carbon storage due to deposition of fossil fuel nitrogen[J].Ecol Appl, 1996, 6(3): 804-814.

    [3] Magill A H, Aber J D, Berntson G M, et al.Long-term nitrogen additions and nitrogen saturation in two temperate forests[J].Ecosystems, 2000(3): 238-253.

    [4] Sonnleitner M A, Günthardt-Goerg M S, Bucher-Wallin I K, et al.Influence of soil type on the effects of elevated atmospheric CO2and N deposition on the water balanceand growth of a young spruce and beech forest[J].Water, Air, Soil Pollut, 2001(126): 271-290.

    [5] Nakaji T, Takenaga S, Kuroha M, et al.Photosynthetic response of Pinus densiflora seedlings to high nitrogen load[J].Environ Sci, 2002, 9(4): 269-282.

    [6] Bauer G A, Bazzaz F A, Minocha R, et al.Effects of chronic N additions on tissue chemistry, photosynthetic capacity, and carbon sequestration potential of a red pine (Pinus resinosa) stand in the NE United States[J].For Ecol Manage, 2004(196): 173-186.

    [7] 湯詩杰,李和平,賀善安.杜仲研究的現(xiàn)狀與展望[J].林業(yè)科技開發(fā),2007,21(2):8-12.

    [8] 潘瑜春,趙春江.地理信息技術(shù)在精準(zhǔn)農(nóng)業(yè)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào),2003,19(4):1-5.

    [9] McKenzie R L, Aucamp P J, Bais A F, et al.Changes in biologically-active ultraviolet radiation reaching the Earth's surface[J].Photochem Photobiol Sci,2007, 6(3): 218-231.

    [10] Yang Z, RaoM N, Elliott N C, et al.Using ground-based multispectral radiometry to detect stress in wheat caused by greenbug (Homoptera: Aphididae) infestation[J].Comput Electr Agric,2005(47): 121-135.

    [11] Rouse Jr J W, Haas R H, Schell J A, et al.Monitoring vegetation systems in the Great Plains with ERTS[J].NASA Spec Publ, 1974(351): 309-317.

    [12] Gitelson A A, Kaufman Y J, Merzlyak M N.Use of a green channel in remote sensing of global vegetation from EOS-MODIS[J].RS Environ, 1996, 58(3): 289-298.

    [13] Lichtenthaler H K, Gitelson A, Lang M.Non-destructive determination of chlorophyll content of leaves of a green and an aurea mutant of tobacco by reflectance measurements[J].J Plant Physiol, 1996, 148(3): 483-493.

    [14] Miller J R, Hare E W, Wu J.Quantitative characterization of the vegetation red edge reflectance.An inverted-Gaussian reflectance model[J].Remote Sens, 1990, 11(10): 1 755-1 773.

    Spectral Response of Eucommia ulmoides Seedlings to Simulated Nitrogen Deposition

    AI Jian-guo1,DU Jiang-li2,JIN Song-heng1,SONG Xin-zhang2
    (1.Tianmu College of Zhejiang A & F University, Lin’an 311300, China; 2.School of Forestry and Biotechnology, Zhejiang A & F University,Lin’an 311300, China)

    Experiments were conducted on spectral response characteristics of Eucommia ulmoides seedlings under different simulated nitrogen deposition like control (CK: 0 kg/ha·a-1), low nitrogen (LN: 30 kg/ha·a-1), medium nitrogen (MN: 60 kg/ha·a-1).The result showed that the seedlings grew well under low nitrogen treatment, with the lowest spectral reflectance, and the red edge position presented a trend of “red shift”.While the spectral reflectance in visible light and near infrared band of seedlings under medium nitrogen deposition was generally higher, and the“blue shift” phenomenon occurred, consequently, the growth of seedlings was restrained but grew well later.The spectral regions of significant difference among different treatments increased with time, while in October, differences reduced.The inverted-Gaussian model parameters such as Rs, R0, λ0, λpand σ could be indicators of spectral response of seedlings under simulated nitrogen deposition, among which, R0 performed the best.Compared with NDVI and gNDVI, RVI700and RVI550could better indicate stress of nitrogen.

    Eucornmia ulmoides; simulated nitrogen deposition; spectral reflectance; red edge

    S718.43

    A

    1001-3776(2014)06-0072-06

    2014-04-25;

    2014-09-20

    浙江省自然科學(xué)基金資助項(xiàng)目(Y3110393)

    哀建國(1968-),男,浙江海寧人,副教授,博士,從事植被生態(tài)學(xué)與植物資源開發(fā)利用研究。

    猜你喜歡
    中氮杜仲植被指數(shù)
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    HPLC法同時(shí)測(cè)定杜仲-當(dāng)歸藥對(duì)中5種成分
    中成藥(2018年8期)2018-08-29 01:28:14
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    略陽杜仲
    清水江水體和沉積物中氮、磷的分布及變化趨勢(shì)
    蒸餾滴定法測(cè)定高溫合金中氮的含量
    殃及池魚
    周末
    西藏濕地生態(tài)系統(tǒng)中氮循環(huán)微生物數(shù)量和多樣性研究
    西藏科技(2015年4期)2015-09-26 12:12:59
    主要植被指數(shù)在生態(tài)環(huán)評(píng)中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    日本a在线网址| 午夜免费鲁丝| 日韩一卡2卡3卡4卡2021年| 中出人妻视频一区二区| 啦啦啦在线免费观看视频4| 一级a爱视频在线免费观看| 久久欧美精品欧美久久欧美| 午夜福利,免费看| 精品免费久久久久久久清纯| 国产区一区二久久| 中文字幕高清在线视频| 国产高清视频在线播放一区| 在线永久观看黄色视频| 国产有黄有色有爽视频| 久久久久国内视频| 视频在线观看一区二区三区| 十八禁人妻一区二区| 久久久久久久久久久久大奶| 一进一出好大好爽视频| 免费av毛片视频| 亚洲视频免费观看视频| 亚洲少妇的诱惑av| 巨乳人妻的诱惑在线观看| 亚洲国产看品久久| 另类亚洲欧美激情| 免费在线观看完整版高清| www.999成人在线观看| 在线天堂中文资源库| 亚洲成av片中文字幕在线观看| 成年版毛片免费区| 黑人操中国人逼视频| 欧美中文综合在线视频| 午夜成年电影在线免费观看| 欧美一级毛片孕妇| 日韩欧美在线二视频| 亚洲国产中文字幕在线视频| 国产人伦9x9x在线观看| 黑人欧美特级aaaaaa片| 9色porny在线观看| 黄色怎么调成土黄色| 中文字幕av电影在线播放| 黑人操中国人逼视频| 欧美乱色亚洲激情| 国产乱人伦免费视频| 婷婷精品国产亚洲av在线| 久久久水蜜桃国产精品网| 18美女黄网站色大片免费观看| 国产精品亚洲av一区麻豆| 看片在线看免费视频| 热re99久久精品国产66热6| 欧美黑人精品巨大| 麻豆av在线久日| 老司机福利观看| 韩国av一区二区三区四区| 99久久久亚洲精品蜜臀av| 免费在线观看完整版高清| 母亲3免费完整高清在线观看| 欧美黄色淫秽网站| 不卡av一区二区三区| 国产深夜福利视频在线观看| 99精国产麻豆久久婷婷| 亚洲免费av在线视频| 国产亚洲精品综合一区在线观看 | 如日韩欧美国产精品一区二区三区| www.精华液| 欧美日韩瑟瑟在线播放| 国产主播在线观看一区二区| 在线观看午夜福利视频| 精品第一国产精品| 一区二区三区激情视频| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 窝窝影院91人妻| 欧美亚洲日本最大视频资源| 琪琪午夜伦伦电影理论片6080| 男人的好看免费观看在线视频 | 午夜福利一区二区在线看| 欧美成人性av电影在线观看| 免费av中文字幕在线| 欧美日韩精品网址| 亚洲第一青青草原| 99在线视频只有这里精品首页| 免费观看精品视频网站| 少妇粗大呻吟视频| 欧美黄色片欧美黄色片| 在线观看免费午夜福利视频| 在线观看免费高清a一片| 女人被躁到高潮嗷嗷叫费观| 亚洲av成人av| 日本撒尿小便嘘嘘汇集6| 国产精品野战在线观看 | 久久久久久亚洲精品国产蜜桃av| 热re99久久国产66热| 国产深夜福利视频在线观看| 亚洲精品中文字幕一二三四区| 最新在线观看一区二区三区| 757午夜福利合集在线观看| 久久国产精品男人的天堂亚洲| 51午夜福利影视在线观看| 首页视频小说图片口味搜索| 国产精品乱码一区二三区的特点 | 日本免费a在线| 国产成人欧美在线观看| 一二三四在线观看免费中文在| 日韩欧美一区视频在线观看| 成人手机av| 欧美久久黑人一区二区| 妹子高潮喷水视频| 91成年电影在线观看| 国产麻豆69| 在线国产一区二区在线| 波多野结衣av一区二区av| 少妇裸体淫交视频免费看高清 | 夜夜看夜夜爽夜夜摸 | 91麻豆av在线| 欧美成人午夜精品| 欧美日韩中文字幕国产精品一区二区三区 | 久久狼人影院| 一边摸一边抽搐一进一小说| 色综合站精品国产| 好看av亚洲va欧美ⅴa在| 女性被躁到高潮视频| 动漫黄色视频在线观看| 色哟哟哟哟哟哟| 精品久久久久久,| 欧美黄色片欧美黄色片| 亚洲av第一区精品v没综合| 如日韩欧美国产精品一区二区三区| www国产在线视频色| 亚洲精品久久午夜乱码| 亚洲色图av天堂| 久久久久亚洲av毛片大全| 亚洲精品美女久久久久99蜜臀| 人成视频在线观看免费观看| 在线观看免费午夜福利视频| 两个人看的免费小视频| 亚洲午夜理论影院| av免费在线观看网站| 悠悠久久av| 80岁老熟妇乱子伦牲交| 欧美日韩中文字幕国产精品一区二区三区 | 欧美性长视频在线观看| 激情视频va一区二区三区| 天天影视国产精品| 美女高潮到喷水免费观看| 久久久久久久久久久久大奶| 国产亚洲欧美98| 国产精品一区二区免费欧美| 免费日韩欧美在线观看| 欧美av亚洲av综合av国产av| 大香蕉久久成人网| 欧美乱码精品一区二区三区| 国产精品国产av在线观看| 日韩大码丰满熟妇| 在线观看舔阴道视频| 亚洲第一av免费看| 香蕉国产在线看| 丁香六月欧美| 日本 av在线| 大香蕉久久成人网| 亚洲专区字幕在线| 欧美黑人精品巨大| 三级毛片av免费| 超碰成人久久| 在线永久观看黄色视频| 很黄的视频免费| 免费在线观看视频国产中文字幕亚洲| 欧美激情 高清一区二区三区| 99精品久久久久人妻精品| 在线国产一区二区在线| 一进一出好大好爽视频| 亚洲成a人片在线一区二区| 日韩精品中文字幕看吧| 国产精品九九99| 性欧美人与动物交配| 老熟妇乱子伦视频在线观看| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 老司机在亚洲福利影院| videosex国产| 国产人伦9x9x在线观看| av天堂在线播放| 满18在线观看网站| 黄片播放在线免费| 在线视频色国产色| 制服诱惑二区| 日韩免费高清中文字幕av| 一区二区日韩欧美中文字幕| 热99re8久久精品国产| 日韩国内少妇激情av| 国产av一区二区精品久久| 亚洲色图 男人天堂 中文字幕| 中文字幕高清在线视频| 久久精品国产清高在天天线| 老司机午夜十八禁免费视频| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看 | 搡老熟女国产l中国老女人| 日韩精品免费视频一区二区三区| 日韩有码中文字幕| 九色亚洲精品在线播放| 黄色视频,在线免费观看| 老司机福利观看| 动漫黄色视频在线观看| 中文字幕人妻丝袜制服| 高潮久久久久久久久久久不卡| 久久欧美精品欧美久久欧美| 超色免费av| 国产成人精品久久二区二区免费| 婷婷六月久久综合丁香| 精品一区二区三区四区五区乱码| 久久香蕉精品热| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 亚洲专区中文字幕在线| 一边摸一边抽搐一进一小说| 国产免费男女视频| 久久香蕉国产精品| 国产片内射在线| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 最近最新中文字幕大全电影3 | 丁香六月欧美| 欧美午夜高清在线| 在线av久久热| 久久精品成人免费网站| 看免费av毛片| 两个人免费观看高清视频| 天堂中文最新版在线下载| 夜夜夜夜夜久久久久| 欧美国产精品va在线观看不卡| 亚洲全国av大片| 香蕉丝袜av| 国产片内射在线| 乱人伦中国视频| 中国美女看黄片| 18美女黄网站色大片免费观看| 免费一级毛片在线播放高清视频 | 国产伦人伦偷精品视频| 最新美女视频免费是黄的| 神马国产精品三级电影在线观看 | 国产精品99久久99久久久不卡| 国产黄a三级三级三级人| 人妻久久中文字幕网| 亚洲中文字幕日韩| 国产又色又爽无遮挡免费看| 日本免费a在线| 宅男免费午夜| 久久亚洲真实| 亚洲avbb在线观看| 欧美人与性动交α欧美软件| 一级a爱片免费观看的视频| 国产精品二区激情视频| 香蕉丝袜av| 成人黄色视频免费在线看| 成人三级做爰电影| 欧美 亚洲 国产 日韩一| 天天躁狠狠躁夜夜躁狠狠躁| 露出奶头的视频| 国产精品秋霞免费鲁丝片| 麻豆一二三区av精品| 纯流量卡能插随身wifi吗| 99在线人妻在线中文字幕| 在线观看免费午夜福利视频| 久久青草综合色| 欧美一区二区精品小视频在线| 午夜福利,免费看| 真人一进一出gif抽搐免费| 丰满迷人的少妇在线观看| 午夜激情av网站| 久久久久久久久免费视频了| 水蜜桃什么品种好| 在线观看一区二区三区| 一级片免费观看大全| 免费久久久久久久精品成人欧美视频| 女性被躁到高潮视频| 国产xxxxx性猛交| 欧美黄色片欧美黄色片| 桃红色精品国产亚洲av| 国产在线精品亚洲第一网站| 日韩中文字幕欧美一区二区| 他把我摸到了高潮在线观看| av网站免费在线观看视频| 成熟少妇高潮喷水视频| 女人高潮潮喷娇喘18禁视频| 香蕉国产在线看| 一本大道久久a久久精品| 亚洲精品美女久久av网站| 大码成人一级视频| 最好的美女福利视频网| 国产国语露脸激情在线看| 亚洲五月色婷婷综合| 俄罗斯特黄特色一大片| 91精品国产国语对白视频| 国产精品久久久人人做人人爽| 最新在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 日本 av在线| 国产极品粉嫩免费观看在线| 在线观看66精品国产| 黑丝袜美女国产一区| 成人亚洲精品av一区二区 | 日本a在线网址| 久久影院123| 999久久久精品免费观看国产| 欧美激情高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 精品国产美女av久久久久小说| 欧美中文日本在线观看视频| 亚洲伊人色综图| 制服诱惑二区| 动漫黄色视频在线观看| 这个男人来自地球电影免费观看| 水蜜桃什么品种好| 男女下面插进去视频免费观看| 色播在线永久视频| 女性生殖器流出的白浆| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 国产精品 欧美亚洲| 久久人人精品亚洲av| 国产欧美日韩一区二区精品| 夜夜看夜夜爽夜夜摸 | 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 成人影院久久| 一进一出好大好爽视频| 精品日产1卡2卡| www.自偷自拍.com| 99热国产这里只有精品6| 国产亚洲av高清不卡| 欧美日韩乱码在线| 日本vs欧美在线观看视频| 黄色 视频免费看| 中国美女看黄片| 国产伦一二天堂av在线观看| 天天添夜夜摸| 可以在线观看毛片的网站| 18禁国产床啪视频网站| 成人18禁在线播放| 男女做爰动态图高潮gif福利片 | 男女高潮啪啪啪动态图| 两个人免费观看高清视频| 精品久久蜜臀av无| 在线视频色国产色| www国产在线视频色| 国产精品影院久久| 热99re8久久精品国产| 丝袜在线中文字幕| 久9热在线精品视频| 一边摸一边做爽爽视频免费| 日韩人妻精品一区2区三区| 老司机深夜福利视频在线观看| 女性生殖器流出的白浆| 麻豆一二三区av精品| 十八禁网站免费在线| www.熟女人妻精品国产| 日韩精品免费视频一区二区三区| 高清av免费在线| 999久久久国产精品视频| 这个男人来自地球电影免费观看| 在线观看66精品国产| 日韩av在线大香蕉| 久久久久国产精品人妻aⅴ院| 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 天堂√8在线中文| 黄片播放在线免费| 亚洲av成人av| 国产激情欧美一区二区| 99在线人妻在线中文字幕| 黑人欧美特级aaaaaa片| 中文字幕人妻丝袜一区二区| 精品国产乱子伦一区二区三区| 亚洲激情在线av| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 国产激情久久老熟女| 亚洲国产精品sss在线观看 | 亚洲国产欧美日韩在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 国产欧美日韩一区二区三区在线| 亚洲七黄色美女视频| 在线观看舔阴道视频| 88av欧美| 黄色怎么调成土黄色| 91麻豆精品激情在线观看国产 | 中文字幕人妻丝袜一区二区| 日韩大尺度精品在线看网址 | 在线观看午夜福利视频| 亚洲欧美激情在线| 嫁个100分男人电影在线观看| 亚洲欧美精品综合一区二区三区| 亚洲av熟女| 国产精品乱码一区二三区的特点 | 夜夜夜夜夜久久久久| av中文乱码字幕在线| 免费搜索国产男女视频| 在线观看舔阴道视频| 大陆偷拍与自拍| 在线观看舔阴道视频| 午夜福利在线免费观看网站| 国产熟女xx| 黑丝袜美女国产一区| 在线国产一区二区在线| 免费高清在线观看日韩| 亚洲国产欧美网| 淫秽高清视频在线观看| 在线观看舔阴道视频| 俄罗斯特黄特色一大片| 国产精品乱码一区二三区的特点 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲全国av大片| 18禁黄网站禁片午夜丰满| 色哟哟哟哟哟哟| 69精品国产乱码久久久| 91字幕亚洲| 日韩有码中文字幕| 欧美+亚洲+日韩+国产| 黄色片一级片一级黄色片| 午夜精品国产一区二区电影| 精品久久久久久电影网| 国产野战对白在线观看| 精品久久蜜臀av无| 一区福利在线观看| 天堂影院成人在线观看| 99精品欧美一区二区三区四区| 精品一区二区三区av网在线观看| 天堂√8在线中文| 国产激情欧美一区二区| 亚洲狠狠婷婷综合久久图片| 女人爽到高潮嗷嗷叫在线视频| 中文字幕最新亚洲高清| 免费av毛片视频| 欧美一级毛片孕妇| 大型黄色视频在线免费观看| 日韩高清综合在线| 亚洲va日本ⅴa欧美va伊人久久| av在线播放免费不卡| 黄色成人免费大全| 男女床上黄色一级片免费看| 三级毛片av免费| 精品福利观看| 国产午夜精品久久久久久| 国产激情久久老熟女| 亚洲少妇的诱惑av| 国产精品久久久久久人妻精品电影| 9色porny在线观看| 成人三级黄色视频| 国产三级黄色录像| 精品国内亚洲2022精品成人| 国产亚洲精品久久久久久毛片| 99国产极品粉嫩在线观看| 90打野战视频偷拍视频| 久久热在线av| 1024香蕉在线观看| 少妇的丰满在线观看| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 777久久人妻少妇嫩草av网站| 极品教师在线免费播放| xxxhd国产人妻xxx| 欧美日本亚洲视频在线播放| 女人被狂操c到高潮| 欧美大码av| 亚洲成人免费av在线播放| 国产熟女午夜一区二区三区| 又紧又爽又黄一区二区| 国产亚洲精品久久久久5区| 国产成人av激情在线播放| 亚洲男人天堂网一区| 久久国产精品人妻蜜桃| 欧美一级毛片孕妇| 自线自在国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 精品少妇一区二区三区视频日本电影| 免费观看人在逋| 一级毛片精品| 欧美老熟妇乱子伦牲交| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费 | 久久精品亚洲熟妇少妇任你| 久久精品影院6| 熟女少妇亚洲综合色aaa.| 国产av在哪里看| 人人妻人人澡人人看| 999久久久国产精品视频| 99热国产这里只有精品6| 三上悠亚av全集在线观看| 两性夫妻黄色片| 手机成人av网站| 91av网站免费观看| 黑人操中国人逼视频| 久久久久久久午夜电影 | 在线免费观看的www视频| 多毛熟女@视频| 丁香欧美五月| 女人被狂操c到高潮| 国产精品影院久久| 日本一区二区免费在线视频| 久久精品91蜜桃| 国产成人影院久久av| av免费在线观看网站| 亚洲一码二码三码区别大吗| 国产成+人综合+亚洲专区| 国产亚洲欧美在线一区二区| 欧美最黄视频在线播放免费 | 久久精品亚洲av国产电影网| 国产精品亚洲一级av第二区| 免费看十八禁软件| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 亚洲精品一区av在线观看| 日韩欧美在线二视频| 亚洲熟妇熟女久久| 欧美国产精品va在线观看不卡| 国产成人影院久久av| 欧美+亚洲+日韩+国产| 一级片'在线观看视频| 欧美黄色片欧美黄色片| 精品久久久精品久久久| 久久久久国内视频| 久久久久亚洲av毛片大全| 亚洲avbb在线观看| 在线天堂中文资源库| 两性午夜刺激爽爽歪歪视频在线观看 | 正在播放国产对白刺激| 视频区欧美日本亚洲| 美女国产高潮福利片在线看| 久久久国产成人免费| 在线观看一区二区三区激情| 最新美女视频免费是黄的| 久久久精品欧美日韩精品| 国产乱人伦免费视频| 99精品久久久久人妻精品| 国产又爽黄色视频| 免费观看精品视频网站| 亚洲精品中文字幕一二三四区| 亚洲精品在线美女| 色婷婷av一区二区三区视频| 亚洲精品一二三| 国产亚洲精品综合一区在线观看 | 激情视频va一区二区三区| 欧美人与性动交α欧美精品济南到| www.自偷自拍.com| 成人18禁在线播放| 亚洲av日韩精品久久久久久密| 亚洲一区中文字幕在线| 国产av又大| 美女大奶头视频| 夜夜躁狠狠躁天天躁| 日本免费a在线| 精品国产亚洲在线| 狂野欧美激情性xxxx| 欧美黄色片欧美黄色片| 中文字幕精品免费在线观看视频| 中文欧美无线码| 日韩成人在线观看一区二区三区| 久久国产精品人妻蜜桃| 欧美激情高清一区二区三区| 热re99久久精品国产66热6| 欧美中文日本在线观看视频| 天堂动漫精品| 热99国产精品久久久久久7| 成人黄色视频免费在线看| 波多野结衣高清无吗| 首页视频小说图片口味搜索| 精品国产一区二区三区四区第35| 欧美黄色淫秽网站| 国产三级在线视频| 久久精品国产99精品国产亚洲性色 | 视频区图区小说| 久久久久亚洲av毛片大全| 性欧美人与动物交配| 欧美av亚洲av综合av国产av| 国产精品久久视频播放| 中文字幕精品免费在线观看视频| 国产精品av久久久久免费| 午夜免费鲁丝| 久久国产精品男人的天堂亚洲| 免费在线观看黄色视频的| 亚洲少妇的诱惑av| 久久久国产成人精品二区 | 久久精品亚洲av国产电影网| 夜夜看夜夜爽夜夜摸 | 淫妇啪啪啪对白视频| 国产黄a三级三级三级人| 亚洲av片天天在线观看| 51午夜福利影视在线观看| 手机成人av网站| 日韩高清综合在线| 麻豆久久精品国产亚洲av | 欧美一区二区精品小视频在线| 日韩一卡2卡3卡4卡2021年| 在线观看一区二区三区激情| 国产激情久久老熟女| 日韩欧美三级三区| 国产精品电影一区二区三区| 大陆偷拍与自拍| 久久久久国产精品人妻aⅴ院| 麻豆av在线久日| 夜夜看夜夜爽夜夜摸 | 日日摸夜夜添夜夜添小说| 亚洲九九香蕉| 老司机在亚洲福利影院| 成年人黄色毛片网站| www.自偷自拍.com| 久久久国产一区二区| 国产精品二区激情视频| 看黄色毛片网站| 午夜亚洲福利在线播放| 侵犯人妻中文字幕一二三四区| 国产国语露脸激情在线看| 国产不卡一卡二| 我的亚洲天堂|