• <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
    亚洲成人免费电影在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美黄色片欧美黄色片| 妹子高潮喷水视频| 制服人妻中文乱码| 桃红色精品国产亚洲av| 日韩三级视频一区二区三区| 久久欧美精品欧美久久欧美| 国语自产精品视频在线第100页| 精品久久久久久成人av| av福利片在线| 欧美性猛交黑人性爽| 国产不卡一卡二| 日本成人三级电影网站| 在线播放国产精品三级| 很黄的视频免费| 宅男免费午夜| 嫩草影视91久久| 亚洲专区字幕在线| 怎么达到女性高潮| 成人亚洲精品av一区二区| 久久性视频一级片| 国产精品一区二区三区四区久久| avwww免费| 黄色 视频免费看| 韩国av一区二区三区四区| 国产精品永久免费网站| 欧美在线一区亚洲| 亚洲av成人av| 村上凉子中文字幕在线| 亚洲一区中文字幕在线| 黑人巨大精品欧美一区二区mp4| 亚洲成av人片在线播放无| www.999成人在线观看| 波多野结衣高清无吗| 91九色精品人成在线观看| 91麻豆av在线| 精品无人区乱码1区二区| 在线观看免费午夜福利视频| 欧美成人性av电影在线观看| 亚洲成人久久爱视频| 免费看十八禁软件| 五月伊人婷婷丁香| 国产亚洲精品久久久久5区| 亚洲专区国产一区二区| 亚洲一区二区三区不卡视频| 国产探花在线观看一区二区| 国产久久久一区二区三区| 狂野欧美激情性xxxx| 岛国视频午夜一区免费看| av天堂在线播放| 国产精品久久久久久亚洲av鲁大| 国产99白浆流出| 亚洲最大成人中文| 亚洲av电影不卡..在线观看| 村上凉子中文字幕在线| 一边摸一边做爽爽视频免费| 不卡av一区二区三区| 性欧美人与动物交配| 琪琪午夜伦伦电影理论片6080| 午夜视频精品福利| 成年版毛片免费区| 亚洲专区字幕在线| 国产99久久九九免费精品| 精品国产乱子伦一区二区三区| 久久人妻av系列| 成人18禁高潮啪啪吃奶动态图| 级片在线观看| 久久人妻av系列| 成年免费大片在线观看| 欧美性猛交╳xxx乱大交人| 露出奶头的视频| 久久久久久久精品吃奶| 青草久久国产| 国产视频内射| 亚洲一区高清亚洲精品| 无人区码免费观看不卡| 欧美色欧美亚洲另类二区| 亚洲成a人片在线一区二区| 给我免费播放毛片高清在线观看| 黄色视频不卡| 久久香蕉国产精品| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品成人综合色| 母亲3免费完整高清在线观看| 免费看美女性在线毛片视频| 色综合站精品国产| 搡老妇女老女人老熟妇| 精品欧美国产一区二区三| 两个人视频免费观看高清| 亚洲精品国产精品久久久不卡| 妹子高潮喷水视频| 狂野欧美激情性xxxx| 黄片大片在线免费观看| a级毛片在线看网站| 大型av网站在线播放| 在线观看美女被高潮喷水网站 | 男女下面进入的视频免费午夜| 国产又色又爽无遮挡免费看| 亚洲激情在线av| 亚洲av美国av| 久久亚洲精品不卡| 久久久国产欧美日韩av| 久热爱精品视频在线9| 一级毛片高清免费大全| 一级毛片精品| aaaaa片日本免费| 午夜精品在线福利| 叶爱在线成人免费视频播放| 亚洲av日韩精品久久久久久密| 久久精品影院6| 色综合亚洲欧美另类图片| 日韩成人在线观看一区二区三区| 欧美在线黄色| 国产精品亚洲一级av第二区| 亚洲欧美精品综合一区二区三区| 啦啦啦韩国在线观看视频| 欧美成人免费av一区二区三区| 午夜亚洲福利在线播放| 九九热线精品视视频播放| 最近最新中文字幕大全免费视频| 最近最新中文字幕大全电影3| 亚洲国产高清在线一区二区三| 亚洲av成人av| 国产视频内射| 精品不卡国产一区二区三区| 欧美性猛交黑人性爽| 又爽又黄无遮挡网站| 婷婷亚洲欧美| 国产av不卡久久| 国产成人av教育| 一区二区三区国产精品乱码| 国产免费男女视频| 老鸭窝网址在线观看| 搡老妇女老女人老熟妇| 午夜精品在线福利| 免费在线观看视频国产中文字幕亚洲| 国产精品一及| 亚洲成av人片免费观看| 国产伦一二天堂av在线观看| 美女黄网站色视频| 亚洲成a人片在线一区二区| 亚洲一区中文字幕在线| 国产av不卡久久| 夜夜夜夜夜久久久久| 亚洲国产精品合色在线| 变态另类丝袜制服| 成人亚洲精品av一区二区| 小说图片视频综合网站| or卡值多少钱| 欧美极品一区二区三区四区| 成人特级黄色片久久久久久久| 国产一区二区在线av高清观看| xxxwww97欧美| 三级男女做爰猛烈吃奶摸视频| www.999成人在线观看| xxxwww97欧美| 亚洲 欧美 日韩 在线 免费| 一夜夜www| 国产精品永久免费网站| 久久欧美精品欧美久久欧美| 中亚洲国语对白在线视频| 夜夜躁狠狠躁天天躁| 午夜福利在线在线| a级毛片在线看网站| 亚洲 欧美一区二区三区| 免费看a级黄色片| 亚洲精华国产精华精| 国产探花在线观看一区二区| 亚洲片人在线观看| 久久久精品国产亚洲av高清涩受| 在线观看66精品国产| 九色成人免费人妻av| 国产成人欧美在线观看| 日韩欧美 国产精品| 亚洲九九香蕉| xxx96com| 国产真人三级小视频在线观看| 91国产中文字幕| 亚洲国产欧美网| 久久精品夜夜夜夜夜久久蜜豆 | 日韩大尺度精品在线看网址| 国产真人三级小视频在线观看| 久久人妻福利社区极品人妻图片| 真人一进一出gif抽搐免费| 日本在线视频免费播放| 不卡一级毛片| 成人午夜高清在线视频| 成人精品一区二区免费| 国产私拍福利视频在线观看| 久久精品国产99精品国产亚洲性色| 久久国产精品影院| 亚洲性夜色夜夜综合| 欧洲精品卡2卡3卡4卡5卡区| 亚洲真实伦在线观看| 最近最新免费中文字幕在线| 国产午夜精品论理片| 1024香蕉在线观看| www国产在线视频色| 亚洲国产高清在线一区二区三| 亚洲性夜色夜夜综合| 级片在线观看| 亚洲成人中文字幕在线播放| 99久久国产精品久久久| 亚洲熟女毛片儿| 天天添夜夜摸| 黑人欧美特级aaaaaa片| 亚洲片人在线观看| 久久久国产成人精品二区| 亚洲中文字幕日韩| 两个人视频免费观看高清| 一级毛片女人18水好多| 婷婷六月久久综合丁香| videosex国产| 国产成人欧美在线观看| 亚洲专区中文字幕在线| 给我免费播放毛片高清在线观看| 亚洲色图av天堂| 变态另类丝袜制服| 午夜久久久久精精品| 人人妻人人澡欧美一区二区| 麻豆成人午夜福利视频| 在线观看免费视频日本深夜| 嫩草影院精品99| 中文字幕最新亚洲高清| 国产精品久久久av美女十八| 免费搜索国产男女视频| 精品高清国产在线一区| 88av欧美| 人妻久久中文字幕网| 少妇粗大呻吟视频| 这个男人来自地球电影免费观看| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 欧美久久黑人一区二区| 国产av一区二区精品久久| 国语自产精品视频在线第100页| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| 身体一侧抽搐| 1024视频免费在线观看| www日本在线高清视频| 亚洲av成人av| 观看免费一级毛片| 久久99热这里只有精品18| 亚洲欧美激情综合另类| 亚洲成人久久性| 欧美黑人巨大hd| 天堂av国产一区二区熟女人妻 | 99在线视频只有这里精品首页| 欧美大码av| 国产乱人伦免费视频| 日韩欧美精品v在线| 一个人免费在线观看电影 | 97超级碰碰碰精品色视频在线观看| 男女那种视频在线观看| 黄频高清免费视频| 在线观看免费午夜福利视频| 亚洲在线自拍视频| 久久国产乱子伦精品免费另类| or卡值多少钱| 一级作爱视频免费观看| 最近最新中文字幕大全电影3| 国产精品亚洲一级av第二区| 99热只有精品国产| 黄片大片在线免费观看| 欧美3d第一页| 亚洲国产欧美网| 日韩精品青青久久久久久| 美女黄网站色视频| 日本一二三区视频观看| 一区二区三区国产精品乱码| 久久久国产欧美日韩av| 日韩av在线大香蕉| 国产人伦9x9x在线观看| 黄色 视频免费看| 90打野战视频偷拍视频| 国产成人欧美在线观看| 欧美日韩精品网址| 丰满的人妻完整版| 亚洲成av人片免费观看| 看片在线看免费视频| 午夜福利18| 久久精品成人免费网站| 一级毛片精品| 日韩成人在线观看一区二区三区| 欧美一区二区国产精品久久精品 | 老熟妇仑乱视频hdxx| 久久精品国产亚洲av香蕉五月| 1024视频免费在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产三级中文精品| 一夜夜www| 亚洲欧洲精品一区二区精品久久久| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 国产精品九九99| 99久久综合精品五月天人人| 麻豆成人av在线观看| 波多野结衣高清作品| xxxwww97欧美| 黑人巨大精品欧美一区二区mp4| 国产一区二区激情短视频| 婷婷丁香在线五月| 亚洲 国产 在线| tocl精华| 亚洲精品在线美女| av超薄肉色丝袜交足视频| 欧美性长视频在线观看| 黄频高清免费视频| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| 啦啦啦观看免费观看视频高清| 91麻豆精品激情在线观看国产| 老司机福利观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费看日本二区| 欧美又色又爽又黄视频| 99久久国产精品久久久| 亚洲av成人不卡在线观看播放网| 国产成人av教育| 欧美一区二区精品小视频在线| 青草久久国产| 老汉色∧v一级毛片| 十八禁人妻一区二区| 国产黄片美女视频| 日韩av在线大香蕉| 国产精品99久久99久久久不卡| 亚洲av第一区精品v没综合| 黑人操中国人逼视频| 久久香蕉国产精品| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 国产精品久久久久久人妻精品电影| 精品日产1卡2卡| 亚洲熟女毛片儿| 精品国产乱码久久久久久男人| 成年女人毛片免费观看观看9| www日本在线高清视频| 在线永久观看黄色视频| 久久香蕉激情| 国产精品乱码一区二三区的特点| 丰满的人妻完整版| 久久久国产成人免费| 欧美精品啪啪一区二区三区| 亚洲精品色激情综合| 成年免费大片在线观看| 久久国产精品影院| 久久中文看片网| 女警被强在线播放| 国产伦在线观看视频一区| 国内揄拍国产精品人妻在线| cao死你这个sao货| 男男h啪啪无遮挡| 久久亚洲精品不卡| 国产成人av教育| 欧美不卡视频在线免费观看 | 免费在线观看亚洲国产| www日本黄色视频网| 午夜福利在线观看吧| 国内久久婷婷六月综合欲色啪| 全区人妻精品视频| 亚洲狠狠婷婷综合久久图片| 特大巨黑吊av在线直播| 亚洲一区二区三区色噜噜| 欧美zozozo另类| netflix在线观看网站| 婷婷六月久久综合丁香| 国产av又大| 日本熟妇午夜| 99久久精品热视频| 亚洲人成电影免费在线| 香蕉久久夜色| 色在线成人网| 久久精品成人免费网站| 窝窝影院91人妻| 日韩精品青青久久久久久| 国产成人精品久久二区二区91| 脱女人内裤的视频| 啦啦啦免费观看视频1| 1024手机看黄色片| 亚洲午夜精品一区,二区,三区| 国产午夜精品久久久久久| 美女扒开内裤让男人捅视频| www.自偷自拍.com| a在线观看视频网站| 亚洲欧美日韩东京热| 久久久水蜜桃国产精品网| 欧美国产日韩亚洲一区| 国产亚洲精品第一综合不卡| 精品欧美国产一区二区三| 777久久人妻少妇嫩草av网站| 亚洲av日韩精品久久久久久密| 欧美大码av| www日本黄色视频网| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 亚洲精华国产精华精| 欧美日韩国产亚洲二区| 中文字幕人妻丝袜一区二区| 女人被狂操c到高潮| 一区二区三区高清视频在线| 中国美女看黄片| 欧美成人午夜精品| 欧美乱妇无乱码| 12—13女人毛片做爰片一| 国产成人aa在线观看| 精品电影一区二区在线| 99riav亚洲国产免费| 一级毛片精品| 日韩大尺度精品在线看网址| 亚洲国产精品999在线| 国产一区二区在线av高清观看| 黑人欧美特级aaaaaa片| 黄片小视频在线播放| 亚洲第一电影网av| av欧美777| 成人国语在线视频| 亚洲中文字幕一区二区三区有码在线看 | 精品一区二区三区四区五区乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜免费激情av| 久久久久国产精品人妻aⅴ院| 国产欧美日韩精品亚洲av| 久久这里只有精品19| 变态另类成人亚洲欧美熟女| 淫秽高清视频在线观看| 一级毛片女人18水好多| 日韩精品免费视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线a可以看的网站| 男女视频在线观看网站免费 | 国产高清有码在线观看视频 | 日韩欧美精品v在线| 国产成人一区二区三区免费视频网站| 在线观看www视频免费| 欧美性猛交黑人性爽| 日韩欧美在线二视频| 成人国产一区最新在线观看| 青草久久国产| 超碰成人久久| 欧美性猛交黑人性爽| 久久热在线av| 国产免费av片在线观看野外av| 精品日产1卡2卡| 麻豆久久精品国产亚洲av| 国产亚洲精品久久久久5区| 精品国产乱子伦一区二区三区| 国产精品久久久久久精品电影| 精品日产1卡2卡| 色在线成人网| 亚洲人成伊人成综合网2020| 成人av一区二区三区在线看| 真人一进一出gif抽搐免费| 欧美日韩瑟瑟在线播放| 亚洲人成77777在线视频| 看黄色毛片网站| 亚洲精品色激情综合| 一本精品99久久精品77| 啦啦啦免费观看视频1| 成人亚洲精品av一区二区| 搡老岳熟女国产| 精品无人区乱码1区二区| 国产精品久久久久久久电影 | av国产免费在线观看| 狂野欧美激情性xxxx| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 亚洲美女黄片视频| 免费在线观看成人毛片| 又粗又爽又猛毛片免费看| 在线观看免费午夜福利视频| ponron亚洲| 午夜精品久久久久久毛片777| 久久中文字幕人妻熟女| 又爽又黄无遮挡网站| 最好的美女福利视频网| 三级毛片av免费| 国产精品综合久久久久久久免费| 毛片女人毛片| 亚洲成人免费电影在线观看| 国产精品 国内视频| 国产1区2区3区精品| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 国产精品久久久人人做人人爽| e午夜精品久久久久久久| av天堂在线播放| 给我免费播放毛片高清在线观看| av免费在线观看网站| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 搡老妇女老女人老熟妇| 欧美性长视频在线观看| 亚洲精品中文字幕一二三四区| 国产亚洲精品av在线| 99精品欧美一区二区三区四区| 亚洲男人天堂网一区| 国产片内射在线| 亚洲欧美日韩无卡精品| xxxwww97欧美| 久久欧美精品欧美久久欧美| 狂野欧美白嫩少妇大欣赏| ponron亚洲| 久9热在线精品视频| 国产午夜福利久久久久久| 亚洲国产精品成人综合色| 国产精品亚洲美女久久久| 亚洲性夜色夜夜综合| 亚洲成av人片在线播放无| 国产一区二区激情短视频| 精品一区二区三区av网在线观看| 欧美黑人欧美精品刺激| svipshipincom国产片| 国产一区在线观看成人免费| 99热6这里只有精品| 三级国产精品欧美在线观看 | 青草久久国产| 精品国产超薄肉色丝袜足j| 午夜免费激情av| 熟女电影av网| 午夜精品久久久久久毛片777| 久久婷婷成人综合色麻豆| 久久性视频一级片| 午夜视频精品福利| 极品教师在线免费播放| 中文字幕最新亚洲高清| 在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 黑人巨大精品欧美一区二区mp4| 狠狠狠狠99中文字幕| 中文资源天堂在线| 日韩精品青青久久久久久| 亚洲成人久久爱视频| 欧美三级亚洲精品| 俺也久久电影网| 超碰成人久久| 免费观看人在逋| 又黄又粗又硬又大视频| 欧美性猛交╳xxx乱大交人| 中文资源天堂在线| 国产精品av久久久久免费| 特级一级黄色大片| 亚洲七黄色美女视频| 国产精品一区二区三区四区免费观看 | 亚洲男人天堂网一区| 高潮久久久久久久久久久不卡| 亚洲精品美女久久久久99蜜臀| 久久热在线av| 国产免费av片在线观看野外av| e午夜精品久久久久久久| 妹子高潮喷水视频| 免费看十八禁软件| 免费在线观看视频国产中文字幕亚洲| 久久香蕉国产精品| 人人妻人人看人人澡| 成人午夜高清在线视频| 男人舔女人的私密视频| 伦理电影免费视频| 精品国产乱码久久久久久男人| 日本黄大片高清| 欧美另类亚洲清纯唯美| 90打野战视频偷拍视频| 亚洲中文字幕日韩| 少妇的丰满在线观看| 九色成人免费人妻av| 777久久人妻少妇嫩草av网站| 国内精品久久久久久久电影| 欧美黄色片欧美黄色片| 欧美精品啪啪一区二区三区| 精品熟女少妇八av免费久了| 黄色女人牲交| 一本综合久久免费| 麻豆成人午夜福利视频| 国产精品av视频在线免费观看| 9191精品国产免费久久| 日本一二三区视频观看| 99re在线观看精品视频| 中文字幕高清在线视频| 男插女下体视频免费在线播放| 精华霜和精华液先用哪个| 成人欧美大片| 这个男人来自地球电影免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 可以在线观看毛片的网站| 欧美成人一区二区免费高清观看 | 制服丝袜大香蕉在线| 99久久久亚洲精品蜜臀av| svipshipincom国产片| 久久中文字幕一级| 每晚都被弄得嗷嗷叫到高潮| 午夜成年电影在线免费观看| 国产精品av久久久久免费| 免费在线观看影片大全网站| 色尼玛亚洲综合影院| 少妇裸体淫交视频免费看高清 | 丝袜人妻中文字幕| 97人妻精品一区二区三区麻豆| 国产激情欧美一区二区| 亚洲精品久久国产高清桃花| 国产在线观看jvid| 天堂av国产一区二区熟女人妻 | 婷婷亚洲欧美| 久久久久性生活片| 长腿黑丝高跟| 三级男女做爰猛烈吃奶摸视频| 99久久国产精品久久久| 免费在线观看完整版高清| 国产伦一二天堂av在线观看| 日韩有码中文字幕| 在线免费观看的www视频| 白带黄色成豆腐渣| 成人国语在线视频|