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

    榆中縣生境質(zhì)量和生境退化的時空分異特征及其地形梯度效應

    2022-08-31 05:55:14屈雯謝保鵬陳英賀瑩馬麗婭
    甘肅農(nóng)業(yè)大學學報 2022年3期
    關鍵詞:榆中縣生境梯度

    屈雯,謝保鵬,陳英,賀瑩,馬麗婭

    (甘肅農(nóng)業(yè)大學管理學院,甘肅 蘭州 730070)

    生境質(zhì)量是生態(tài)系統(tǒng)為個體和種群提供適宜、持續(xù)性生存條件的能力[1],是區(qū)域生物多樣性及生態(tài)服務水平的重要反映[2-3]。近年來,隨著人為活動的不斷加強,使得土地利用方式、利用強度以及景觀格局發(fā)生了巨大的變化[4],導致區(qū)域生態(tài)系統(tǒng)結構和功能的變化,進而引發(fā)了生物多樣性銳減、生境質(zhì)量下降、生境退化等一系列生態(tài)環(huán)境問題[3]。地形作為自然環(huán)境中的重要影響因子,對地表物質(zhì)的遷移與能量的轉換有著顯著影響[2],研究表明地形因子是影響土地利用景觀格局變化的重要因素[3,5],同時也是影響生境系統(tǒng)物質(zhì)交換和能量循環(huán)的重要因素[4]。因此,從地形梯度效應角度出發(fā),對生境質(zhì)量及生境退化進行研究,將有助于準確的識別生境質(zhì)量的變化效應,對于挖掘區(qū)域生境質(zhì)量的地形信息圖譜具有重要的意義。目前,國內(nèi)外學者進行了大量關于生境質(zhì)量的評估研究,國外學者側重于研究生境斑塊對生物多樣性的影響[6-9],國內(nèi)學者側重于研究生境質(zhì)量的時空變化規(guī)律[10]??傮w而言,生境質(zhì)量的評估方法分為兩大類[1],一是實地調(diào)查法,通過對研究區(qū)進行實地調(diào)研,構建評價指標對生境質(zhì)量進行綜合評價;二是基于土地利用數(shù)據(jù),采用模型對研究區(qū)生境質(zhì)量進行評估,使用最普遍的是InVEST模型,此類研究尺度涉及全國[10]、?。?1-13]、市[14-15]縣[16-17]、流域[18-19]、城市群[20]、自然保護區(qū)[21-22]等,研究內(nèi)容多從生境質(zhì)量時空變化[4,23]、影響因素[24-25]、生境質(zhì)量與土地整治[26]、生境質(zhì)量與生態(tài)系統(tǒng)服務價值[27]等方面展開。

    位于甘肅省中部地區(qū)的榆中縣是東部季風區(qū)、西北干旱半干旱區(qū)、青藏高原區(qū)三大自然區(qū)的交匯區(qū),其地形地勢和土地利用方式具有明顯的過渡性特征[28],特殊的地理位置使其生境具有顯著的地域性特征。近年來,隨著蘭州市的不斷發(fā)展,其作為蘭州市“東擴”的主要承載區(qū),經(jīng)濟發(fā)展勢頭迅猛,人口不斷增加,資源與環(huán)境約束日益凸顯,因此,以榆中縣作為研究區(qū),基于土地利用數(shù)據(jù)、DEM數(shù)據(jù),借助InVEST模型評估榆中縣生境退化和生境質(zhì)量,以期為榆中縣生境保護、生態(tài)保護及國土空間規(guī)劃編制提供理論參考。

    1 材料與方法

    1.1 研究區(qū)概況

    榆中縣位于甘肅省中部地區(qū),是省會蘭州市的東部郊區(qū)(E 103°49′15″~104°34′40″,N 35°34′20″~36°26′30″),縣域海拔1 369~3 670 m,總面積為32 946.69 km2。地處黃土高原西部,同時也是西北干旱區(qū)、青藏高寒區(qū)和東部季風區(qū)的交匯地帶[28]??h域地形復雜,南部為高寒陰濕區(qū),馬銜山和興隆山從西北向東南橫列整個縣域的南部,該區(qū)降水充沛,植被茂盛,森林資源豐富,是榆中縣主要的生態(tài)保障區(qū);中部為興隆山山前沖積、洪積平原和宛川河谷地組成的川塬河谷平原區(qū),地勢平坦,土壤肥沃,交通便利,是榆中縣主要的工業(yè)和農(nóng)業(yè)生產(chǎn)區(qū);北部為干旱山區(qū),降水稀少,植被稀疏,是榆中縣水土流失最為嚴重的區(qū)域。

    1.2 數(shù)據(jù)來源

    (1)榆中縣DEM數(shù)據(jù)來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心,其分辨率為30 m。

    (2)榆中縣2008年和2018年土地利用柵格數(shù)據(jù)來源于榆中縣土地利用變更數(shù)據(jù)庫,分辨率為30 m。根據(jù)研究需要,按照《土地利用現(xiàn)狀分類GB/T 21010—2017》的分類依據(jù)[29]和土地利用/土地覆被遙感監(jiān)測數(shù)據(jù)分類標準[10],利用ArcGIS 10.2軟件平臺對變更數(shù)據(jù)中的地類進行進一步劃分,分為耕地、園地、林地、草地、水域、建設用地和未利用地7大類。

    1.3 研究方法

    1.3.1 基于In VEST模型的生境質(zhì)量和生境退化評估 In VEST模型中開展生境質(zhì)量評價首先要計算區(qū)域的生境退化度,其計算公式如下:

    式中:Dxj是生境退化指數(shù),R是脅迫因子的數(shù)量;Y r是脅迫因子r在土地利用柵格圖中的總個數(shù);ry是土地利用柵格單元上的脅迫因子個數(shù);ωr是脅迫因子的權重;βx是柵格x的可達性水平;Sjr是土地利用類型j對脅迫因子r的敏感性程度,取值范圍為0~1;irxy是脅迫因子的影響距離;d xy為柵格x與y之間的距離;d rmax是脅迫因子r的最大影響范圍。

    生境質(zhì)量的計算公式為:

    式中:Qxj為第j種土地利用類型中柵格單元x的生境質(zhì)量;H j為第j種土地利用類型的生境適宜;z為尺度常數(shù)是模型的默認參數(shù);k為半飽和常數(shù)。

    In VEST模型運行時需要輸入?yún)?shù),本研究在對InVEST模型手冊及相關研究[11-12,27-29]分析的基礎之上,參數(shù)設置具體參考劉春芳等[30]、賀瑩等[31]的研究,參數(shù)如表1~2所示。

    表1 研究區(qū)內(nèi)各威脅因子的參數(shù)Table 1 Parameters of threat factors in the study area

    1.3.2 圖譜分析 文章借助土地利用圖譜分析的相關原理[12,32-34],構建生境圖譜分析模型,以記錄和分析不同時空條件下生境質(zhì)量及生境退化的復合變化信息。計算公式如下:

    式中:A為整個研究期內(nèi)生境質(zhì)量或生境退化圖譜單元變化類型;A 1、A 2為研究初期和研究末期的生境質(zhì)量或生境退化圖譜柵格狀態(tài)編碼。

    1.3.3 地形梯度分析 地形條件在很大程度上會影響土地利用類型的分布[35,2],本研究引進復合地形因子,即地形位梯度[36],以綜合反映地形條件與生境質(zhì)量、生境退化之間的關系。計算公式如下:

    式中:T表示地形位梯度指數(shù),E表示研究區(qū)內(nèi)某一點的高程,S表示研究區(qū)某一點的坡度,Eˉ表示研究區(qū)的平均高程,Sˉ表示研究區(qū)的平均坡度。

    由于不同地形位梯度下生境質(zhì)量、生境退化面積存在較大的差異,本文為消除這種面積大小差異帶來的量綱影響[36],引進地形分布指數(shù),對榆中縣生境質(zhì)量及生境退化的地形位梯度效應進行研究。參照周潤芳等[37]對榆中縣地形地貌的研究,將地形位梯度分為5級,其計算公式如下:

    式中:P為分布指數(shù),Sie為第i級生境質(zhì)量或生境退化度在e級地形因子等級下的面積,Si為第i級生境質(zhì)量或生境退化度的總面積,Se為第e級地形因子所占的面積,S為研究區(qū)的總面積。

    表2 各地類對脅迫因子的敏感度參數(shù)Table 2 Sensitivity parameters of each category to stress factors

    2 結果與分析

    2.1 生境質(zhì)量和生境退化度的時空變化特征

    2.1.1 生境質(zhì)量時空變化特征 為了準確識別榆中縣生境質(zhì)量的變化規(guī)律,研究利用ArcGIS軟件中的自然斷點法將生境質(zhì)量分為低值區(qū)(1級)、較低值區(qū)(2級)、中值區(qū)(3級)、較高值區(qū)(4級)、高值區(qū)(5級)5個等級。

    由圖1可以看出,2008~2018年,榆中縣生境質(zhì)量的空間格局整體上較為穩(wěn)定,從南到北總體呈現(xiàn)出由西南向東北遞減的分布特點。生境質(zhì)量高值區(qū)域主要分布在縣域南部的興隆山、馬銜山沿線,區(qū)內(nèi)人為活動稀少,常年降水充足,植被茂密,土地利用多以草地、林地為主,生境質(zhì)量整體上優(yōu)于其他區(qū)域。生境質(zhì)量較高值區(qū)域主要分布在生境質(zhì)量高值區(qū)域的外圍,除此之外在貢井林場有小部分分布。生境質(zhì)量中值區(qū)域整體上分布零散,在園子岔鄉(xiāng)、貢井鄉(xiāng)、中連川鄉(xiāng)、清水驛鄉(xiāng)、韋營鄉(xiāng)等個別鄉(xiāng)鎮(zhèn)分布較為集中,除此之外在宛川河沿岸河谷兩側也有零星分布。生境質(zhì)量較低值區(qū)域位于北部山區(qū)和中部平原地區(qū),該區(qū)土地利用主要是耕地和其他草地、農(nóng)村居民點。低值區(qū)域位于縣域中部和北部,該區(qū)域分布大量的居民點用地、城鎮(zhèn)用地以及工礦用地,人類活動頻繁,因此生境質(zhì)量普遍偏低。由表3可以看出,2008年和2018年研究區(qū)平均生境質(zhì)量分別為0.228 9和0.226 6,生境質(zhì)量總體上略有下降。研究期間1等生境占比由46.29%變?yōu)?7.01%,結合圖2可以發(fā)現(xiàn),變化區(qū)域主要是由于受到建設用地需求增加的影響,人類對區(qū)域生態(tài)環(huán)境有所干預,破壞了原有的自然生境,使得生境質(zhì)量略有退化;2等生境占比由38.74%變?yōu)?8.06%,主要是由于研究期間耕地面積減少的影響;3等、4等、5等生境總體上保持相對穩(wěn)定的狀態(tài)。

    圖1 2008~2018年榆中縣生境質(zhì)量空間分布圖Figure 1 The spatial distribution of habitat quality in Yuzhong County from 2008 to 2018

    表3 2008~2018年榆中縣生境質(zhì)量變化情況Table 3 Habitat quality changes in Yuzhong County from 2008 to 2018

    圖2 榆中縣土地利用類型Figure 2 Land use types in Yuzhong County

    2.1.2 生境退化度時空變化特征 由圖3可以看出,2008~2018年,榆中縣生境退化度的空間格局整體上較為穩(wěn)定,并且表現(xiàn)出明顯的空間差異,具體表現(xiàn)為以宛川河為界,呈現(xiàn)“河北低而河南高”的分布特點。生境退化度最嚴重的區(qū)域為中部川區(qū)的和平鎮(zhèn)、定遠鎮(zhèn)等重點鄉(xiāng)鎮(zhèn),這些區(qū)域是榆中縣主要的城鎮(zhèn)及工農(nóng)業(yè)發(fā)展區(qū),如圖2所示,隨著縣域人為活動日益增強,對居住地、工礦用地以及交通運輸用地等建設用地的需求快速增加,使原有的自然生境遭到破壞,導致生境破碎度加劇、連通性變差,最終使得這些區(qū)域的生境質(zhì)量下降;整個研究期間縣域生境退化程度有所增加,平均生境退化度由2008年的0.007 370變?yōu)?018年的0.007 389。具體而言,研究期間生境退化度高值區(qū)域進一步擴散,2008年和2018年的占比分別為0.15%、0.49%,說明榆中縣中部川區(qū)存在生態(tài)退化的風險;生境退化度較低值區(qū)有所增加,占比由2008年的67.35%變?yōu)?018年的74.26%,主要分布在縣域北部石質(zhì)低山丘陵區(qū),主要是由于榆中縣近年來實施了一系列“退耕還林”工程和“生態(tài)北山”建設工程,加強了對北部山區(qū)生態(tài)環(huán)境的保護;生境退化較高值區(qū)域有所縮減,由2008年的9.20%變?yōu)?018年的7.54%,變化主要發(fā)生在縣域中部,轉化成了生境退化中值區(qū);生境退化中值區(qū)域也有所縮減,由2008年的23.30%變?yōu)?018年的7.54%,變化主要集中在縣域北部,主要轉化成了生境退化低值區(qū)域。

    2.2 生境質(zhì)量及生境退化的圖譜分析

    按照前文分析,將生境質(zhì)量5級區(qū)間的區(qū)間編碼設置為1、2、3、4、5,分別對應生境質(zhì)量低值、較低值、中值、較高值和高值5種狀態(tài)。由圖4可以看出,榆中縣生境質(zhì)量轉移圖譜共有25種類型,其中,有5類生境質(zhì)量圖譜單元未發(fā)生變化,其他20類生境質(zhì)量圖譜單元或多或少都發(fā)生了變化。整體而言,未發(fā)生變化的圖譜編碼為11(表示研究期初和期末的生境質(zhì)量狀態(tài)都為低值)、22、33、44、55總占比為71.45%,其中,低-低(編碼為11)的圖譜單元最為穩(wěn)定,占未變化單元的49.75%,主要分布在縣域北部和中部偏西的區(qū)域,和草地、未利用地的分布基本一致,中-中(編碼為33)的圖譜單元最不穩(wěn)定,占未變化單元的2%,分布比較零散,主要分布在縣域中東部區(qū)域。發(fā)生變化的圖譜占總圖譜的比例為28.55%,其中,“低-較低”和“較低-低”(編碼為12和21)圖譜變化最為明顯,由圖3可以看出,這兩類圖譜主要分布在榆中縣北部、東部,以及宛川河以南興隆山以北的區(qū)域,縣域北部和東部分布著大量的耕地,隨著人類耕地利用程度的不斷加深,使得原有的耕地景觀遭到破壞,生境質(zhì)量有所下降,宛川河以南興隆山以北的區(qū)域地勢平坦,土壤肥沃,是主要的生產(chǎn)生活場所,尤其是近幾年,縣域人口增加導致建設用地需求增長,從而影響了生境變化。由圖4可以看出,榆中縣生境退化轉移圖譜共有19種類型,其中,有4類生境退化圖譜單元未發(fā)生變化,其他15類生境退化圖譜單元或多或少都發(fā)生了變化。整體而言,未發(fā)生變化的圖譜編碼為22、33、44、55,總占比為81.10%,其中,較低-較低(編碼為22)的圖譜單元最為穩(wěn)定,占未變化單元的76.93%。發(fā)生變化的圖譜占總圖譜的比例為18.90%,其中,中-較低(編碼為32)圖譜變化最為明顯,集中分布在縣域北部,主要是“北山生態(tài)工程”的建設使得該區(qū)中等退化轉化成了較低程度的退化。

    圖4 2008~2018年榆中縣生境質(zhì)量及生境退化的圖譜空間分布Figure 4 The spatial distribution map of habitat quality and habitat degradation in Yuzhong County from 2008 to 2018

    表4 2008~2018年榆中縣生境退化變化情況Table 4 Habitat degradation changes in Yuzhong County from 2008 to 2018

    圖3 2008~2018年榆中縣生境退化度空間分布Figure 3 The spatial distribution of habitat degradation in Yuzhong County from 2008 to 2018

    2.3 生境質(zhì)量的地形梯度效應

    2.3.1榆中縣的地形梯度特征 研究上述公式(5)計算榆中縣地形位梯度,結果區(qū)間為0.22~1.09,采用自然斷點法將其分成5個等級,等級1表示地形位指數(shù)在0.22~0.42,等級2表示地形位指數(shù)在0.42~0.54,等級3表示地形位指數(shù)在0.54~0.64,等級4表示地形位指數(shù)在0.64~0.75,等級5表示地形位指數(shù)在0.75~1.09,等級越低表示地形越平坦。

    由圖5可以看出,榆中縣地形位指數(shù)呈現(xiàn)出“南部高、中部低、北部高”的分布格局,與縣域南高北低,中間低凹的地勢基本一致。榆中縣地形位以3級為主,占比為30.80%左右,其次是4級、2級、1級,占比分別為24.23%、16.94%、15.55%左右,5級地形位占比最小,為12.48%左右。由表5可知,榆中縣1等地形位和2等地形位的土地利用類型以耕地為主,占比達60%左右,3等、4等及5等地形位以未利用地為主,具體而言耕地、園地、水域、建設用地在1等地形位上占比最大,林地、草地在5等地形位上占比最大,未利用地在4等地形位上占比最大。隨著地形位指數(shù)的增加,耕地、園地、水域、建設用地的占比出現(xiàn)下降趨勢,相反,林地、草地、未利用地的占比出現(xiàn)了增加的趨勢。主要是由于人類活動強度最大的區(qū)域一般都是地形位低的位置,地形位指數(shù)值低表明區(qū)域自然環(huán)境越有利于人類生存,人類對土地的利用程度越強,相反,隨著地形位指數(shù)的增加,人類對自然環(huán)境的改造程度越小,因此土地利用多以林地、草地以及未利用地為主。

    表5 不同地形位梯度的土地利用特征Table 5 Land use characteristics of different terrain position gradients

    圖5 榆中縣地形位梯度分布情況Figure 5 Topographical position gradient distribution in Yuzhong County

    2.3.2 生境質(zhì)量及生境退化的地形位梯度分布特征地形要素是影響生物種群分布格局以及維持生物多樣性的一個重要因素[4],地形分異性會影響土地利用類型進而也會在很大程度上影響生境質(zhì)量及生境退化的空間分布格局。本文在計算了研究區(qū)地形位指數(shù)的基礎上,采用分布指數(shù)來描述生境質(zhì)量、生境退化在不同地形位梯度上的分布特征,以揭示綜合地形變化對生境質(zhì)量和生境退化的影響。如圖6、7所示,生境質(zhì)量和生境退化在地形位梯度上呈現(xiàn)出較大的差異。具體而言,低質(zhì)量生境在第3、第4級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,較低質(zhì)量生境在第1、第2、第3級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,中等質(zhì)量生境在第3、第4級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢。較高質(zhì)量生境及高質(zhì)量生境在第4、第5級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢。生境退化低值在第1級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,生境退化較低值在第3、第4、第5級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,生境退化中等區(qū)域在第1、第2、第3級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,較高生境退化和高生境退化第1級地形位梯度上的分布指數(shù)大于1,占絕對優(yōu)勢,因此,第3級地形位梯度是生境質(zhì)量及生境退化的混合梯度。

    圖6 生境質(zhì)量在地形梯度上的分布指數(shù)Figure 6 Distribution index of habitat quality on topographic gradient

    3 討論

    本研究應用InVEST模型和地形位指數(shù)對榆中縣生境質(zhì)量進行了比較充分的分析,其中生境質(zhì)量的時空分布結果與劉春芳等[11,28,30]的研究結果相近,說明本文的研究結果比較符合榆中縣的實際情況,應該能反映近年來榆中縣生境質(zhì)量及生境退化的變化,而且本研究首次從地形視角對榆中縣的生境質(zhì)量及生境退化的分布特點進行了分析,較好地揭示了研究區(qū)生境質(zhì)量的地形梯度效應,研究成果可為榆中縣生態(tài)環(huán)境保護、生物多樣性維護提供理論依據(jù)。但研究也存在以下的不足:(1)研究在對生境質(zhì)量進行評估時,應用的是傳統(tǒng)的評估模型,缺乏方法的創(chuàng)新性,而且模型計算中參數(shù)設置可能存在一定的主觀性;(2)生境質(zhì)量和生境退化會受到多種因素的影響,本文只是選取了較為重要的地形因素,對人為因素的研究還需進一步探討。

    圖7 生境退化度在地形梯度上的分布指數(shù)Figure 7 Distribution index of habitat degradation degree on topographic gradient

    榆中縣是西北地區(qū)重要的生態(tài)屏障區(qū),作為三大自然帶的交匯地帶具有明顯的地形分異特征,區(qū)域生境質(zhì)量及生境退化受到地形變化特征的影響??h域南部地勢高聳,分布大量的天然植被,林地、草地是該區(qū)的土地利用類型,是縣域生境質(zhì)量的高值區(qū);中部宛川河流域地勢低平,是近年來蘭州市東擴的主要承載區(qū),同時該區(qū)域是主要的城鄉(xiāng)建設用地增加區(qū),隨著人類活動日益加強,導致原有的土地利用結構、生態(tài)結構和功能遭到破壞,從而使得此處生境退化程度不斷加深,生境質(zhì)量不斷下降;北部地區(qū)系主要的黃土丘陵區(qū),區(qū)域自然環(huán)境較差,分布大量的草地和未利用地以及部分農(nóng)村居民點。因此,在南部山區(qū),應該加強對森林的保護,適當?shù)脑黾恿植莞采w度,同時有序控制人類活動,以保持其生態(tài)功能;在中部川區(qū),應該限制建設用地增長,同時要遏制生境退化,對用地結構進行適當調(diào)整,尤其是園地要盡可能的降低其利用強度;北部山區(qū),一方面要適當?shù)脑黾恿植莞采w率,另一方面要控制現(xiàn)有的居民點用地,盡可能使該區(qū)生態(tài)保持自然的狀態(tài),以使其進行必要的自我調(diào)節(jié)和恢復。

    4 結論

    本研究利用DEM數(shù)據(jù)、2008年和2018年土地利用柵格數(shù)據(jù),基于InVEST模型和地形位指數(shù),評估了榆中縣生境質(zhì)量時空的變化特征及其地形梯度效應。結果表明:1)近10 a來,榆中縣生境質(zhì)量整體上呈下降趨勢,生境退化度有小幅度的增加趨勢,但二者整體上變化幅度較小,同時生境質(zhì)量及生境退化未發(fā)生變化的圖譜單元占比遠遠大于發(fā)生變化的圖譜單元占比,說明榆中縣生境整體上存在惡化的風險,但這種風險比較小。

    2)研究期內(nèi)生境質(zhì)量總體上呈現(xiàn)出由西南向東北遞減的分布特點,生境退化度呈現(xiàn)出“河北低而河南高”的分布格局。

    3)榆中縣生境質(zhì)量及生境退化分布指數(shù)具有明顯的地形位梯度變化效應,二者都隨著地形位起伏發(fā)生變化,但生境質(zhì)量和生境退化的優(yōu)勢地形位梯度差異較大,高質(zhì)量生境在第4、5等地形位梯度上分布最多,占絕對優(yōu)勢,生境退化高值在1等地形位梯度上占絕對優(yōu)勢,二者的混合梯度都出現(xiàn)在第3級地形位梯度。

    猜你喜歡
    榆中縣生境梯度
    甘肅榆中縣掏9億元建公園 ,7年沒交付
    榆中縣關工委開展“義務理發(fā)”志愿服務活動
    下一代英才(2023年3期)2023-04-18 01:05:50
    一個改進的WYL型三項共軛梯度法
    棗樹適應干旱生境研究進展
    落葉果樹(2021年6期)2021-02-12 01:29:10
    一種自適應Dai-Liao共軛梯度法
    一類扭積形式的梯度近Ricci孤立子
    金秀瑤族長鼓舞傳承的文化生境
    體育科技(2016年2期)2016-02-28 17:06:07
    甘肅省榆中縣農(nóng)村居民社會經(jīng)濟狀況與高血壓關系調(diào)查
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    地溫梯度判定地熱異常的探討
    河南科技(2014年3期)2014-02-27 14:05:45
    一进一出抽搐动态| 国产精品久久久久久人妻精品电影| 50天的宝宝边吃奶边哭怎么回事| 日韩av在线大香蕉| 亚洲精品在线观看二区| 国产亚洲精品综合一区在线观看 | 激情在线观看视频在线高清| 精品久久久久久久久久免费视频| 99久久国产精品久久久| 一级,二级,三级黄色视频| 国产一区二区三区在线臀色熟女| 黄色女人牲交| 亚洲国产毛片av蜜桃av| 亚洲熟女毛片儿| 无限看片的www在线观看| 亚洲av片天天在线观看| 99精品在免费线老司机午夜| 精品午夜福利视频在线观看一区| e午夜精品久久久久久久| 禁无遮挡网站| cao死你这个sao货| 亚洲无线在线观看| 亚洲视频免费观看视频| 97人妻精品一区二区三区麻豆 | 午夜久久久在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲人成电影观看| 黑丝袜美女国产一区| 日本精品一区二区三区蜜桃| 黄色女人牲交| 无遮挡黄片免费观看| 国产一级毛片七仙女欲春2 | 丝袜美足系列| 久久午夜综合久久蜜桃| 夜夜夜夜夜久久久久| 久久香蕉激情| 亚洲无线在线观看| 亚洲男人天堂网一区| 午夜影院日韩av| 每晚都被弄得嗷嗷叫到高潮| 国产精品99久久99久久久不卡| 久久精品国产综合久久久| 亚洲色图综合在线观看| 国内精品久久久久久久电影| 亚洲色图 男人天堂 中文字幕| 免费女性裸体啪啪无遮挡网站| 欧美日韩中文字幕国产精品一区二区三区 | 日本vs欧美在线观看视频| 这个男人来自地球电影免费观看| 国产精品久久久人人做人人爽| 久久人人97超碰香蕉20202| 曰老女人黄片| 91精品国产国语对白视频| 非洲黑人性xxxx精品又粗又长| 国产精品精品国产色婷婷| 午夜福利,免费看| 村上凉子中文字幕在线| 亚洲精品国产区一区二| 99久久久亚洲精品蜜臀av| 一边摸一边抽搐一进一出视频| 50天的宝宝边吃奶边哭怎么回事| 国产高清videossex| www国产在线视频色| 亚洲一区高清亚洲精品| 两人在一起打扑克的视频| 极品人妻少妇av视频| 麻豆国产av国片精品| 在线观看日韩欧美| 香蕉国产在线看| 亚洲五月婷婷丁香| 高清黄色对白视频在线免费看| 亚洲性夜色夜夜综合| 欧美激情 高清一区二区三区| 午夜a级毛片| 午夜激情av网站| 国产乱人伦免费视频| 美女国产高潮福利片在线看| 无限看片的www在线观看| 国产精品乱码一区二三区的特点 | 亚洲激情在线av| 亚洲国产看品久久| 国产又爽黄色视频| 久热爱精品视频在线9| 亚洲va日本ⅴa欧美va伊人久久| 欧美 亚洲 国产 日韩一| 国产伦一二天堂av在线观看| 91在线观看av| 亚洲精品av麻豆狂野| 一区二区三区国产精品乱码| 一级毛片精品| 91精品国产国语对白视频| 国内精品久久久久精免费| 9色porny在线观看| 精品午夜福利视频在线观看一区| 欧美丝袜亚洲另类 | 黄色成人免费大全| 欧美乱色亚洲激情| 欧美绝顶高潮抽搐喷水| 亚洲最大成人中文| 久久久久久久久中文| 两个人看的免费小视频| 国产亚洲欧美精品永久| 成人三级做爰电影| 午夜久久久久精精品| 人妻久久中文字幕网| 老司机午夜福利在线观看视频| 成年女人毛片免费观看观看9| 成人国语在线视频| 狂野欧美激情性xxxx| 亚洲一区高清亚洲精品| 可以在线观看的亚洲视频| 精品熟女少妇八av免费久了| 久久精品91无色码中文字幕| 国产精品久久视频播放| 69av精品久久久久久| 亚洲av熟女| 成人免费观看视频高清| 亚洲精品国产精品久久久不卡| 19禁男女啪啪无遮挡网站| 十八禁网站免费在线| 天堂动漫精品| 长腿黑丝高跟| 日本精品一区二区三区蜜桃| 91精品三级在线观看| av视频免费观看在线观看| 啦啦啦免费观看视频1| av在线播放免费不卡| 少妇被粗大的猛进出69影院| 日本 av在线| 国产在线观看jvid| 嫩草影视91久久| 国产精品久久久久久亚洲av鲁大| 波多野结衣高清无吗| 一级a爱片免费观看的视频| 国产一区在线观看成人免费| 国产高清视频在线播放一区| 一进一出抽搐gif免费好疼| 欧美一区二区精品小视频在线| 97碰自拍视频| 亚洲成人久久性| 999精品在线视频| 女人被躁到高潮嗷嗷叫费观| 人人妻,人人澡人人爽秒播| 国产又色又爽无遮挡免费看| 国产欧美日韩精品亚洲av| 亚洲一码二码三码区别大吗| 午夜福利在线观看吧| 久久久国产成人免费| √禁漫天堂资源中文www| 黄片播放在线免费| 国产精品电影一区二区三区| 好看av亚洲va欧美ⅴa在| 久久久久国产精品人妻aⅴ院| 一区二区三区高清视频在线| www.自偷自拍.com| 免费少妇av软件| 亚洲男人天堂网一区| 少妇裸体淫交视频免费看高清 | 99久久综合精品五月天人人| 最新美女视频免费是黄的| 久久久久久久久中文| 精品卡一卡二卡四卡免费| 国产精品国产高清国产av| 亚洲 欧美 日韩 在线 免费| av中文乱码字幕在线| 免费在线观看完整版高清| 国产一区二区在线av高清观看| 满18在线观看网站| a在线观看视频网站| 国产精华一区二区三区| 91国产中文字幕| 一进一出抽搐动态| 精品久久蜜臀av无| 淫妇啪啪啪对白视频| 亚洲精品av麻豆狂野| www日本在线高清视频| 亚洲人成网站在线播放欧美日韩| 国产高清视频在线播放一区| 人人妻人人澡欧美一区二区 | 麻豆一二三区av精品| 国产精品 国内视频| 国产午夜福利久久久久久| 叶爱在线成人免费视频播放| 成人18禁在线播放| 久久亚洲精品不卡| 激情在线观看视频在线高清| 黑人巨大精品欧美一区二区蜜桃| 国产免费av片在线观看野外av| 日韩欧美免费精品| 波多野结衣一区麻豆| 午夜久久久在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品在线美女| 国产主播在线观看一区二区| 成人亚洲精品一区在线观看| 国产成人影院久久av| 日本vs欧美在线观看视频| 久久午夜综合久久蜜桃| 搞女人的毛片| 国产高清视频在线播放一区| 精品国产一区二区三区四区第35| 欧美激情极品国产一区二区三区| 色精品久久人妻99蜜桃| 国产极品粉嫩免费观看在线| 人人澡人人妻人| 欧美日韩瑟瑟在线播放| 男人舔女人的私密视频| 久久精品国产亚洲av香蕉五月| 我的亚洲天堂| 精品久久久久久久久久免费视频| 香蕉久久夜色| 国产99白浆流出| 亚洲国产看品久久| 亚洲精品美女久久av网站| 性色av乱码一区二区三区2| 午夜福利视频1000在线观看 | 一级片免费观看大全| 午夜福利一区二区在线看| 国产精品,欧美在线| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲av嫩草精品影院| 日韩欧美在线二视频| 午夜免费鲁丝| 老熟妇乱子伦视频在线观看| 黄色成人免费大全| 成人特级黄色片久久久久久久| 免费少妇av软件| 久热这里只有精品99| 88av欧美| 一二三四社区在线视频社区8| 国产精品日韩av在线免费观看 | 99久久久亚洲精品蜜臀av| 亚洲美女黄片视频| 国产亚洲精品久久久久久毛片| 国语自产精品视频在线第100页| 亚洲最大成人中文| 黄色a级毛片大全视频| 亚洲av成人一区二区三| 一级毛片精品| 国产成+人综合+亚洲专区| 精品福利观看| 精品国产超薄肉色丝袜足j| 精品国产亚洲在线| 一a级毛片在线观看| 18禁国产床啪视频网站| 在线天堂中文资源库| 成人手机av| 校园春色视频在线观看| 波多野结衣巨乳人妻| 成人国产综合亚洲| 又紧又爽又黄一区二区| 午夜影院日韩av| 亚洲,欧美精品.| 免费高清视频大片| 免费看美女性在线毛片视频| 国产精品久久久久久人妻精品电影| 成年人黄色毛片网站| 欧美日韩一级在线毛片| 国产精品秋霞免费鲁丝片| 国产高清激情床上av| 久久婷婷成人综合色麻豆| 亚洲 国产 在线| 长腿黑丝高跟| 欧美日韩亚洲综合一区二区三区_| 90打野战视频偷拍视频| 国产一区二区三区综合在线观看| 国产午夜精品久久久久久| 香蕉国产在线看| 欧美丝袜亚洲另类 | 亚洲视频免费观看视频| 神马国产精品三级电影在线观看 | 女性生殖器流出的白浆| 午夜福利视频1000在线观看 | 九色国产91popny在线| 欧美一级毛片孕妇| 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| 久久人妻av系列| 波多野结衣巨乳人妻| 国产精品乱码一区二三区的特点 | 热re99久久国产66热| 欧美日韩黄片免| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av香蕉五月| 精品国产一区二区三区四区第35| 久久中文看片网| 免费人成视频x8x8入口观看| 麻豆av在线久日| 亚洲成a人片在线一区二区| 国产欧美日韩一区二区精品| 女人精品久久久久毛片| 中出人妻视频一区二区| 波多野结衣巨乳人妻| 国产亚洲精品av在线| 在线观看66精品国产| 丝袜在线中文字幕| 亚洲av美国av| 国产一区二区激情短视频| 免费搜索国产男女视频| 亚洲五月婷婷丁香| 亚洲成a人片在线一区二区| 亚洲天堂国产精品一区在线| 日韩精品中文字幕看吧| 国产主播在线观看一区二区| 国产午夜精品久久久久久| 高清黄色对白视频在线免费看| 男人的好看免费观看在线视频 | 天天躁狠狠躁夜夜躁狠狠躁| 国产主播在线观看一区二区| 制服丝袜大香蕉在线| 十八禁人妻一区二区| 桃色一区二区三区在线观看| 夜夜夜夜夜久久久久| 99久久综合精品五月天人人| 欧美最黄视频在线播放免费| 黄网站色视频无遮挡免费观看| 婷婷六月久久综合丁香| av网站免费在线观看视频| 老司机午夜福利在线观看视频| 91精品国产国语对白视频| 乱人伦中国视频| 久久婷婷人人爽人人干人人爱 | 999久久久国产精品视频| 国产欧美日韩一区二区三区在线| 天天躁夜夜躁狠狠躁躁| 国产男靠女视频免费网站| 巨乳人妻的诱惑在线观看| 变态另类丝袜制服| 一个人观看的视频www高清免费观看 | 日本欧美视频一区| 精品人妻在线不人妻| 99精品欧美一区二区三区四区| 欧美黄色淫秽网站| 在线av久久热| 日韩欧美在线二视频| 久久久久久久久中文| 老熟妇乱子伦视频在线观看| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| av超薄肉色丝袜交足视频| www.熟女人妻精品国产| 美国免费a级毛片| 亚洲国产高清在线一区二区三 | 热99re8久久精品国产| 国产亚洲精品av在线| 国产高清视频在线播放一区| 一区福利在线观看| www.www免费av| 99久久综合精品五月天人人| 热re99久久国产66热| 曰老女人黄片| 亚洲av成人不卡在线观看播放网| 欧美老熟妇乱子伦牲交| 久久精品国产亚洲av高清一级| 久久九九热精品免费| 中文字幕av电影在线播放| 国产精品免费视频内射| 亚洲精品在线观看二区| 免费女性裸体啪啪无遮挡网站| 一级作爱视频免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 丁香六月欧美| 人人妻,人人澡人人爽秒播| 在线天堂中文资源库| 日本 欧美在线| x7x7x7水蜜桃| 亚洲情色 制服丝袜| 熟妇人妻久久中文字幕3abv| 久久亚洲精品不卡| 高清黄色对白视频在线免费看| 老司机午夜十八禁免费视频| 一级毛片高清免费大全| 精品久久久精品久久久| 一a级毛片在线观看| 国产精品av久久久久免费| 亚洲人成电影免费在线| xxx96com| 嫩草影视91久久| 变态另类丝袜制服| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| 精品高清国产在线一区| 又黄又粗又硬又大视频| 亚洲精品久久成人aⅴ小说| 免费观看人在逋| 热99re8久久精品国产| 亚洲色图综合在线观看| 欧美日韩精品网址| 日韩欧美国产一区二区入口| 欧美午夜高清在线| 满18在线观看网站| 欧美黄色淫秽网站| 欧美日本视频| 国产一区二区三区视频了| 好男人电影高清在线观看| 亚洲精品中文字幕一二三四区| 性少妇av在线| 看黄色毛片网站| 欧美一区二区精品小视频在线| 18禁国产床啪视频网站| 别揉我奶头~嗯~啊~动态视频| 9色porny在线观看| 日韩精品中文字幕看吧| 久久久国产欧美日韩av| 亚洲av电影在线进入| 操出白浆在线播放| 免费久久久久久久精品成人欧美视频| 国产成+人综合+亚洲专区| 久久人人精品亚洲av| 色综合站精品国产| 可以在线观看的亚洲视频| 国产精品国产高清国产av| 免费少妇av软件| 久久人妻熟女aⅴ| 亚洲人成伊人成综合网2020| 国产精品国产高清国产av| 精品国产一区二区三区四区第35| 在线观看免费日韩欧美大片| 黄色丝袜av网址大全| 精品高清国产在线一区| 大型黄色视频在线免费观看| 久久久久久久午夜电影| 亚洲成a人片在线一区二区| 人人澡人人妻人| 日韩中文字幕欧美一区二区| 日韩三级视频一区二区三区| 在线观看免费日韩欧美大片| 国产激情久久老熟女| 手机成人av网站| 国产又色又爽无遮挡免费看| 制服诱惑二区| 亚洲精品av麻豆狂野| 少妇 在线观看| АⅤ资源中文在线天堂| 在线观看日韩欧美| 黄色a级毛片大全视频| 波多野结衣av一区二区av| 午夜亚洲福利在线播放| 狂野欧美激情性xxxx| 9191精品国产免费久久| 国产精品香港三级国产av潘金莲| 日本免费a在线| 欧美另类亚洲清纯唯美| 韩国av一区二区三区四区| 国产精品免费一区二区三区在线| 日日干狠狠操夜夜爽| 免费久久久久久久精品成人欧美视频| 咕卡用的链子| 国产精品二区激情视频| 大香蕉久久成人网| 亚洲欧美精品综合久久99| 亚洲精品一区av在线观看| 日本黄色视频三级网站网址| 国产精品二区激情视频| 91av网站免费观看| 国产精品电影一区二区三区| 激情视频va一区二区三区| 亚洲久久久国产精品| 精品一区二区三区四区五区乱码| 香蕉国产在线看| 久久久久久久久中文| 中文亚洲av片在线观看爽| 亚洲第一青青草原| 十分钟在线观看高清视频www| 一本综合久久免费| 久久国产乱子伦精品免费另类| 久久精品人人爽人人爽视色| 欧美国产日韩亚洲一区| 亚洲精品一卡2卡三卡4卡5卡| 国产精品,欧美在线| 亚洲一区高清亚洲精品| 欧美性长视频在线观看| 在线视频色国产色| АⅤ资源中文在线天堂| 极品人妻少妇av视频| 亚洲av电影不卡..在线观看| av免费在线观看网站| 免费av毛片视频| 一个人观看的视频www高清免费观看 | 两个人视频免费观看高清| 国产亚洲欧美精品永久| 激情视频va一区二区三区| 男女下面插进去视频免费观看| 精品欧美一区二区三区在线| 久久精品亚洲熟妇少妇任你| 在线观看免费视频网站a站| 欧美色欧美亚洲另类二区 | 国产av在哪里看| 在线视频色国产色| 久99久视频精品免费| 丁香六月欧美| 成人国产一区最新在线观看| 国产欧美日韩精品亚洲av| 欧美中文综合在线视频| 国产乱人伦免费视频| 国产高清激情床上av| 亚洲五月天丁香| 午夜成年电影在线免费观看| 亚洲自偷自拍图片 自拍| 日韩免费av在线播放| 黑丝袜美女国产一区| a在线观看视频网站| 色综合欧美亚洲国产小说| 欧美激情久久久久久爽电影 | 亚洲av美国av| 欧美日韩亚洲国产一区二区在线观看| 国产精品98久久久久久宅男小说| 国产私拍福利视频在线观看| 夜夜夜夜夜久久久久| 成人国产综合亚洲| 精品一品国产午夜福利视频| 免费无遮挡裸体视频| 女性被躁到高潮视频| 免费av毛片视频| 欧美成人免费av一区二区三区| www日本在线高清视频| 色综合亚洲欧美另类图片| a在线观看视频网站| 老司机福利观看| 亚洲午夜精品一区,二区,三区| 一个人免费在线观看的高清视频| 欧美一级毛片孕妇| 国产午夜福利久久久久久| 国产精品九九99| 最新美女视频免费是黄的| 亚洲欧美激情在线| 亚洲成人国产一区在线观看| 午夜久久久在线观看| 欧美日韩黄片免| 国产精品 国内视频| 国产在线精品亚洲第一网站| 黄色 视频免费看| 亚洲五月婷婷丁香| 久久国产亚洲av麻豆专区| 精品人妻1区二区| 精品国产乱码久久久久久男人| 后天国语完整版免费观看| 国产亚洲精品一区二区www| 久久精品aⅴ一区二区三区四区| 亚洲五月天丁香| 亚洲国产精品合色在线| 99国产精品免费福利视频| 美女 人体艺术 gogo| 中文亚洲av片在线观看爽| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品在线观看二区| 中文字幕另类日韩欧美亚洲嫩草| 国产激情久久老熟女| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕人妻丝袜一区二区| 日韩有码中文字幕| 中文字幕高清在线视频| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 最好的美女福利视频网| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| 亚洲国产精品sss在线观看| 亚洲狠狠婷婷综合久久图片| 99久久精品国产亚洲精品| 黑人操中国人逼视频| 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| 久久九九热精品免费| 国产亚洲欧美精品永久| 色哟哟哟哟哟哟| 亚洲精品美女久久久久99蜜臀| 久久久久九九精品影院| 91在线观看av| 此物有八面人人有两片| 欧美激情 高清一区二区三区| 欧美丝袜亚洲另类 | 黄色a级毛片大全视频| 在线国产一区二区在线| 欧美亚洲日本最大视频资源| 俄罗斯特黄特色一大片| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 最近最新中文字幕大全免费视频| 国产精品久久电影中文字幕| 亚洲精品国产区一区二| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 美女午夜性视频免费| 国产av又大| 91九色精品人成在线观看| av在线天堂中文字幕| 淫妇啪啪啪对白视频| 国产三级在线视频| av片东京热男人的天堂| 又紧又爽又黄一区二区| 国产精品日韩av在线免费观看 | 成年女人毛片免费观看观看9| 欧美日韩精品网址| 男女做爰动态图高潮gif福利片 | 精品无人区乱码1区二区| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 亚洲精品国产一区二区精华液| 日本免费一区二区三区高清不卡 | 两性夫妻黄色片| 精品国产亚洲在线| 中文字幕另类日韩欧美亚洲嫩草| 国产高清videossex| 久久精品国产清高在天天线| 最新在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 一级毛片高清免费大全| 在线视频色国产色| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线观看二区| 又大又爽又粗| av电影中文网址| 欧美老熟妇乱子伦牲交|