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

    生態(tài)安全條件下土地利用格局優(yōu)化
    ——以皇甫川流域?yàn)槔?/h1>
    2014-10-11 10:42:58李曉兵
    生態(tài)學(xué)報(bào) 2014年12期
    關(guān)鍵詞:皇甫蓋度模數(shù)

    喻 鋒,李曉兵,王 宏

    (1.北京師范大學(xué)地表過(guò)程與資源生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,資源學(xué)院,北京 100875;2.國(guó)土資源部信息中心,國(guó)土資源部國(guó)土資源戰(zhàn)略研究重點(diǎn)實(shí)驗(yàn)室,北京 100812)

    生態(tài)安全條件下土地利用格局優(yōu)化
    ——以皇甫川流域?yàn)槔?/p>

    喻 鋒1,2,李曉兵1,*,王 宏1

    (1.北京師范大學(xué)地表過(guò)程與資源生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,資源學(xué)院,北京 100875;2.國(guó)土資源部信息中心,國(guó)土資源部國(guó)土資源戰(zhàn)略研究重點(diǎn)實(shí)驗(yàn)室,北京 100812)

    在人類(lèi)活動(dòng)導(dǎo)致的生態(tài)環(huán)境問(wèn)題中,土地利用格局變化對(duì)區(qū)域生態(tài)安全起著決定性作用。選擇我國(guó)乃至在世界上都罕見(jiàn)的多沙、粗沙、強(qiáng)烈水土流失的黃河皇甫川流域作為研究區(qū)域,針對(duì)土壤侵蝕和生態(tài)用水這兩大影響流域生態(tài)安全的關(guān)鍵問(wèn)題,開(kāi)展土地利用變化模擬、土壤侵蝕3S監(jiān)測(cè)、適宜植被蓋度估算、土壤侵蝕情景模擬、土地利用格局優(yōu)化等方面的綜合研究。結(jié)果表明:(1) 土地利用格局變化是導(dǎo)致土壤侵蝕變化的重要因素。從1987年到2015年,建筑用地、林地和灌叢面積增加,水域、沙地和裸砒砂巖面積減少,導(dǎo)致流域年平均土壤侵蝕模數(shù)由16160.72 t km-2a-1降至9710.72 t km-2a-1。(2) 土地利用類(lèi)型與土壤侵蝕模數(shù)密切相關(guān)并存在一定規(guī)律。6種不同土地利用類(lèi)型多年平均土壤侵蝕模數(shù)的大小順序?yàn)椋郝闩皫r>沙地>耕地>草地>林地>灌叢,表明灌、林措施是流域植被恢復(fù)和土壤侵蝕減少的首選。(3) 土地利用結(jié)構(gòu)和空間格局優(yōu)化有助于提升區(qū)域生態(tài)安全水平。在優(yōu)化格局下,2007年和2015年土壤侵蝕模數(shù)比優(yōu)化前分別減少5469.57 t km-2a-1和5432.77 t km-2a-1,優(yōu)化后2015年土壤侵蝕模數(shù)僅為4277.95 t km-2a-1,低于流域自然侵蝕臨界值或允許土壤侵蝕模數(shù)5300 t km-2a-1,優(yōu)化土地利用格局成為加強(qiáng)流域生態(tài)環(huán)境建設(shè)的重要途徑。

    生態(tài)安全; 土地利用; 土壤侵蝕; 適宜植被蓋度; 格局優(yōu)化

    隨著全球性環(huán)境問(wèn)題的日益顯著,生態(tài)安全問(wèn)題越來(lái)越受到人類(lèi)的普遍關(guān)注。近十多年來(lái),國(guó)內(nèi)外生態(tài)安全研究主要集中在基因工程生物的生態(tài)風(fēng)險(xiǎn)與生態(tài)安全、生態(tài)入侵、生態(tài)系統(tǒng)健康評(píng)價(jià)、生態(tài)安全評(píng)價(jià)方法及模型和生態(tài)安全預(yù)警等方面[1]。區(qū)域生態(tài)安全狀況受自然因素和人類(lèi)活動(dòng)雙重影響,而在人類(lèi)活動(dòng)導(dǎo)致的生態(tài)環(huán)境問(wèn)題中,土地利用變化過(guò)程又對(duì)區(qū)域生態(tài)安全起著決定性作用[2]。從土地利用或土地生態(tài)系統(tǒng)視角來(lái)探討區(qū)域生態(tài)安全問(wèn)題,成為土地利用/覆蓋變化研究和生態(tài)安全研究的熱點(diǎn)。

    土地利用格局優(yōu)化,即不同類(lèi)型土地利用單元在時(shí)間和空間尺度上的優(yōu)化配置,可以協(xié)調(diào)生態(tài)效益、經(jīng)濟(jì)效益和社會(huì)效益三者之間的關(guān)系。長(zhǎng)期以來(lái),土地利用格局研究側(cè)重于對(duì)土地利用類(lèi)型和數(shù)量的優(yōu)化配置以實(shí)現(xiàn)經(jīng)濟(jì)效益的最大化,即追求生態(tài)安全中的經(jīng)濟(jì)生態(tài)安全目標(biāo)。近年來(lái),生態(tài)因素已成為土地利用優(yōu)化配置的重要約束條件和優(yōu)化目標(biāo),生態(tài)安全的理念也開(kāi)始在土地利用結(jié)構(gòu)優(yōu)化中得以體現(xiàn):Herrmann等[3]應(yīng)用系統(tǒng)工程方法,從土壤肥力、地下水質(zhì)量、地表水、群落生境和景觀五方面選取指標(biāo)作為生態(tài)約束條件,開(kāi)展鄉(xiāng)村土地利用結(jié)構(gòu)優(yōu)化設(shè)計(jì);Makowski等[4]以歐共體農(nóng)用土地資源面臨的最主要的污染問(wèn)題為導(dǎo)向,以氮流失量最小為規(guī)劃目標(biāo),建立了歐共體農(nóng)業(yè)土地利用結(jié)構(gòu)優(yōu)化模型;徐學(xué)選等[5]應(yīng)用線性規(guī)劃模型,以土壤侵蝕量作為生態(tài)約束條件,探討了黃土丘陵區(qū)生態(tài)建設(shè)中農(nóng)林牧土地結(jié)構(gòu)的優(yōu)化模式;林彰平[6]提出了生態(tài)安全條件下土地利用模式優(yōu)化研究的概念框架,并采用灰色線性規(guī)劃模型,探討了東北農(nóng)牧交錯(cuò)帶土地利用結(jié)構(gòu)優(yōu)化方法;劉艷芳等[7]對(duì)基于綠當(dāng)量的最佳森林覆蓋率標(biāo)準(zhǔn)的生態(tài)優(yōu)化方法進(jìn)行了探討,對(duì)土地利用優(yōu)化中關(guān)于生態(tài)標(biāo)準(zhǔn)的量化研究有著重要的啟發(fā)意義;蘇偉等[8]建立了土地利用格局優(yōu)化模擬模型LUOS,基于GLP方法和CA模型,得出了中國(guó)北方農(nóng)牧交錯(cuò)帶在6種情景下的土地利用優(yōu)化方案;呂紅峰等[9]對(duì)呼倫貝爾沙地北部海拉爾沙帶的土地利用進(jìn)行了生態(tài)安全評(píng)價(jià)和優(yōu)化,建立了海拉爾沙帶生態(tài)安全條件下的土地利用結(jié)構(gòu)優(yōu)化方案;粱偉恒等[10]運(yùn)用主成分分析法、聚類(lèi)分析法和灰色線性規(guī)劃法,提出了重慶市開(kāi)縣的土地利用優(yōu)化模式。

    為此,本文選擇我國(guó)乃至在世界上都罕見(jiàn)的多沙、粗沙、強(qiáng)烈水土流失的黃河皇甫川流域作為研究區(qū)域,針對(duì)土壤侵蝕和生態(tài)用水兩大影響流域生態(tài)安全的關(guān)鍵問(wèn)題,以土地利用變化模擬、土壤侵蝕3S監(jiān)測(cè)、適宜植被蓋度估算、土壤侵蝕情景模擬、土地利用格局優(yōu)化等為重點(diǎn),開(kāi)展生態(tài)安全條件下土地利用格局優(yōu)化研究,有利于為合理開(kāi)發(fā)國(guó)土空間、科學(xué)配置“三生”(生產(chǎn)、生活和生態(tài))空間提供科學(xué)依據(jù),對(duì)大力推進(jìn)生態(tài)文明建設(shè)、實(shí)現(xiàn)區(qū)域可持續(xù)發(fā)展具有重要戰(zhàn)略意義。

    1 研究區(qū)

    皇甫川流域位于E110.3°—111.2°,N39.2°—39.9°之間,地跨鄂爾多斯高原與黃土高原,流域全長(zhǎng)約125 km,總面積約為3240 km2?;矢Υ饔?qū)贉貛О敫珊荡箨懶詺夂颍昶骄鶜鉁?.2—7.2 ℃;年平均降雨量在379—420 mm之間,主要集中在夏季,6—8月份的降雨量占總降水量的61%。流域內(nèi)砒砂巖大面積裸露,原生植被破壞殆盡,再加上較大的地勢(shì)高差和較強(qiáng)的暴雨沖刷導(dǎo)致水土流失十分嚴(yán)重,年平均向黃河輸送約0.5 億t泥沙,深刻影響著當(dāng)?shù)氐纳鷳B(tài)環(huán)境和農(nóng)牧業(yè)生產(chǎn)條件。土壤母質(zhì)類(lèi)型為砒砂巖土類(lèi)、黃土類(lèi)、風(fēng)沙土類(lèi)、沖積土類(lèi),分別占流域面積的39.8%、41.3%、17.0%、1.9%,土壤類(lèi)型相應(yīng)為栗鈣土、黃綿土、風(fēng)沙土、草甸土和潮土,其中大部分屬于易侵蝕類(lèi)型,這也決定了該區(qū)域的土壤侵蝕比較強(qiáng)烈。

    皇甫川流域?qū)儆谵r(nóng)牧交錯(cuò)帶,主要的土地利用類(lèi)型為耕地、草地和林地,形成了以草地利用為主、耕地鑲嵌、小片林地星散分布、侵蝕溝網(wǎng)嵌套的土地利用格局。其中,流域中西部和北部人口較為稀少,又有大面積連片分布的天然草地,適合于開(kāi)展草地畜牧業(yè);而東南部人口相對(duì)稠密,具備了開(kāi)展農(nóng)業(yè)生產(chǎn)的勞動(dòng)力資源,加之地形破碎、土地開(kāi)墾程度高,經(jīng)濟(jì)發(fā)展以農(nóng)業(yè)為主。耕地主要集中在流域內(nèi)沙圪堵鄉(xiāng)、納林鄉(xiāng)與府谷縣境內(nèi),其中坡耕地占到耕地總面積的46%,但大多數(shù)坡耕地土壤含水量差、淋溶作用強(qiáng)、土壤肥力低、平均單產(chǎn)低。

    2 方法

    2.1 土地利用變化模擬

    2.1.1 數(shù)量結(jié)構(gòu)——Markov模型

    作為一種基于柵格的空間概率模型,在假設(shè)當(dāng)前社會(huì)經(jīng)濟(jì)模式發(fā)展不變情況下,通過(guò)計(jì)算土地利用變化的轉(zhuǎn)移概率矩陣,Markov(馬爾柯夫)模型可以預(yù)測(cè)特定情景下土地利用數(shù)量的未來(lái)變化趨勢(shì)。傳統(tǒng)的Markov模型可表示為:

    P(n)=P(n-1)P=P(0)Pn

    (1)

    式中,P(0)和P(n)分別表示初始時(shí)刻和任何時(shí)刻土地利用的狀態(tài)概率向量;P為土地利用類(lèi)型之間相互轉(zhuǎn)化的轉(zhuǎn)移概率矩陣,也可以把它看成是多個(gè)時(shí)段的平均轉(zhuǎn)移概率向量,即:

    (2)

    式中,P1…Pn為各個(gè)時(shí)段的轉(zhuǎn)移概率矩陣,n1…nn為各個(gè)時(shí)段的步長(zhǎng)。

    本文通過(guò)計(jì)算1987—1995年、1995—2000年、2000—2007年3個(gè)時(shí)段的平均轉(zhuǎn)移概率矩陣,預(yù)測(cè)未來(lái)直至2015年皇甫川流域土地利用的變化趨勢(shì)。

    2.1.2 空間格局——Bayes模型

    Bayes(貝葉斯)模型能彌補(bǔ)Markov模型不能準(zhǔn)確刻畫(huà)未來(lái)土地利用空間格局的缺陷。基于Bayes最大似然概率原則,該模型考慮某一柵格轉(zhuǎn)化為不同土地利用類(lèi)型的可能性,從而最終確定土地利用的未來(lái)空間格局。其計(jì)算方法為:

    (3)

    本文采用5×5的柵格窗口計(jì)算每一柵格轉(zhuǎn)變?yōu)椴煌恋乩妙?lèi)型的概率,基于“耕地—草地—灌叢—林地—水域—建筑用地—裸砒砂巖—沙地”的順序,按照概率的高低依次確定各土地利用類(lèi)型柵格的空間位置,直至滿足Markov模型預(yù)測(cè)后的數(shù)量要求為止,最終確定皇甫川流域2015年土地利用空間格局。

    2.2土壤侵蝕變化模擬[11]

    黃土高原地區(qū)地表破碎且降雨集中,水蝕起著關(guān)鍵作用;此外,重力侵蝕與水蝕相互復(fù)合,也是造成當(dāng)?shù)販衔g強(qiáng)烈的重要?jiǎng)恿?。研究成果表明[12- 13]:水蝕和重力侵蝕約占皇甫川流域土壤侵蝕總量的83%;其中,當(dāng)坡度小于35°時(shí),重力侵蝕可以忽略不計(jì),而當(dāng)坡度大于35°時(shí),重力侵蝕才顯著增大,約為水蝕量的30%。因此,基于計(jì)算土壤水蝕的USLE(通用土壤侵蝕方程),皇甫川流域土壤侵蝕模型的表達(dá)式為:

    (4)

    式中,E為年均土壤侵蝕量(t hm-2a-1);R為降雨侵蝕力因子(MJ mm hm-2h-1a-1);K為土壤可蝕性因子(t hm2h hm-2MJ-1mm-1);L為坡長(zhǎng)因子;S為坡度因子;C為地表覆蓋因子;P為水土保持措施因子;θ為地表坡度。

    上述6個(gè)參數(shù)分別有不同的計(jì)算方法(表1)。在預(yù)測(cè)未來(lái)土地利用情景下的土壤侵蝕時(shí),參數(shù)C是在未來(lái)時(shí)點(diǎn)土地利用格局和現(xiàn)狀時(shí)點(diǎn)植被蓋度狀況的基礎(chǔ)上得到的,參數(shù)R則采用流域的多年平均值。

    表1 皇甫川流域土壤侵蝕模型參數(shù)計(jì)算[13- 15]

    2.3適宜植被蓋度計(jì)算[16- 17]

    皇甫川流域?qū)儆诟珊蛋敫珊档貐^(qū),水分是植物生長(zhǎng)的主要限制因子。經(jīng)過(guò)長(zhǎng)期的水土流失綜合治理,流域植被蓋度有明顯增加,但隨之植被需水量也日益增加。因此,植被與水的矛盾,即水分生態(tài)安全成為影響流域生態(tài)安全的重要因素。研究流域的水分生態(tài)安全,應(yīng)考慮目前植被的正常生長(zhǎng)能否維持或能維持多久,不同林、灌、草等土地利用配置模式的水分豐欠差異,以及主要植物種的種群密度是否合理等。適宜植被蓋度是指在適宜的狀態(tài)下,即植物群落的水分輸入和輸出達(dá)到平衡時(shí),群落中主要植被的蓋度。研究適宜植被蓋度是解決水分生態(tài)安全問(wèn)題的重要突破口,有助于優(yōu)化皇甫川流域的能量流動(dòng)和物質(zhì)循環(huán),緩解水土流失。

    皇甫川流域喬灌木和半灌木的適宜植被蓋度的計(jì)算公式為:

    冠幅面積=π×(冠幅長(zhǎng)×冠幅寬)/4

    單位面積植物年耗水量=蒸騰強(qiáng)度×葉干重×蒸騰時(shí)數(shù)/1000

    單株年耗水量=單位面積年耗水量×冠幅面積

    單株?duì)I養(yǎng)面積=單株年耗水量/單位土壤面積的有效水分

    群落適宜蓋度=冠幅面積/營(yíng)養(yǎng)面積×100%

    (5)

    皇甫川流域小半灌木和草本植物的計(jì)算公式為:

    單位面積內(nèi)株(叢)面積=單位面積×蓋度

    單位面積植物年耗水量=蒸騰強(qiáng)度×單位面積內(nèi)葉干重×蒸騰時(shí)數(shù)/1000

    營(yíng)養(yǎng)面積=單位面積年耗水量×單位面積/單位土壤面積的有效水分

    群落適宜蓋度=株(叢)面積/營(yíng)養(yǎng)面積×100%

    (6)

    式中,冠幅面積是當(dāng)植物冠幅的垂直投影區(qū)為橢圓形時(shí)所對(duì)應(yīng)的面積。

    單位土壤面積上的有效水分值按兩種假設(shè)條件計(jì)算:(1) 無(wú)效降水量及蒸發(fā)量占降水量的10%,其余降水量均可被植物利用,無(wú)地下水補(bǔ)給;(2) 無(wú)效降水量及蒸發(fā)量占降水量的10%,年降水量的20%作為徑流流失,無(wú)地下水補(bǔ)給。

    2.4 土壤侵蝕情景模擬

    將皇甫川流域土地利用的可能情景簡(jiǎn)化為全林地、全灌叢、全草地、全耕地這4種單一類(lèi)型覆蓋格局,分別模擬其在植被覆蓋現(xiàn)狀和水分平衡(即各植被類(lèi)型蓋度為其適宜植被蓋度)條件下所對(duì)應(yīng)的土壤侵蝕,從而為下一步基于土壤侵蝕量進(jìn)行土地利用格局空間優(yōu)化提供定量依據(jù)。主要步驟包括:

    (1) 從現(xiàn)有植被類(lèi)型與蓋度出發(fā),分別計(jì)算出各柵格單元在全林地、全灌叢、全草地、全耕地4種情況下的土壤侵蝕。

    (2) 根據(jù)其現(xiàn)有水分條件和植被的實(shí)際蒸騰量,分別推算各柵格單元作為林地、灌叢、草地和耕地時(shí)各自的適宜蓋度(采用2.3的計(jì)算結(jié)果)。

    (3) 根據(jù)適宜植被蓋度的計(jì)算結(jié)果,分別計(jì)算出水分平衡條件下不同土地利用格局情景下(全林地、全灌叢、全草地、全耕地)所對(duì)應(yīng)的土壤侵蝕模數(shù)。

    2.5 土地利用格局優(yōu)化

    2.5.1 土地利用類(lèi)型總量控制

    土地利用類(lèi)型的總量控制主要考慮四個(gè)方面的約束條件:一是強(qiáng)制約束性因素,即要滿足強(qiáng)制性邊界條件要求,如濕地、森林公園、地質(zhì)公園等自然保護(hù)區(qū)(本文不存在),以及湖泊和水庫(kù)等不允許或不能進(jìn)行覆被調(diào)整的地區(qū);二是政策導(dǎo)向性因素,如坡度大于25°的陡坡耕地實(shí)施退耕還林還草及城鎮(zhèn)擴(kuò)展等規(guī)劃建設(shè)用地需求等;三是經(jīng)濟(jì)協(xié)調(diào)性因素,即要滿足糧食安全、畜牧業(yè)和林果業(yè)發(fā)展、增加農(nóng)民純收入等的要求,既確保較為穩(wěn)定的耕地面積,同時(shí)適當(dāng)提高林地、草地的比例;四是生態(tài)保障性因素,即滿足生態(tài)優(yōu)化和生物多樣性要求,如考慮林地、灌叢、草地的適宜比例,同時(shí)減少裸砒砂巖和沙地的面積。

    2.5.2 土地利用類(lèi)型空間優(yōu)化

    空間格局優(yōu)化主要依據(jù)圖像中各柵格單元對(duì)不同覆被類(lèi)型的歸屬概率而進(jìn)行。歸屬概率是指某一柵格,從有利于生態(tài)安全(本文主要考察減少土壤侵蝕)的角度出發(fā),轉(zhuǎn)化為不同植被類(lèi)型的可能性,其計(jì)算方法如下[18]:

    (7)

    土地利用類(lèi)型空間優(yōu)化的具體步驟如下:

    (1) 首先維持現(xiàn)有的水域空間格局不發(fā)生變化,同時(shí)滿足對(duì)建筑用地空間分配的預(yù)期要求。

    (2) 從流域現(xiàn)有植被覆蓋與土地利用類(lèi)型出發(fā),推算水分供給條件,以此作為進(jìn)行土地覆蓋優(yōu)化布局的關(guān)鍵約束因子之一。

    (3) 對(duì)于25°以上的坡地,根據(jù)其現(xiàn)有水分條件和植被實(shí)際蒸騰量,分別推算柵格作為林地、灌叢、草地時(shí)各自的覆蓋度,并且確保其蓋度不超過(guò)各自類(lèi)型的適宜蓋度。從減少土壤侵蝕模數(shù)的角度,分別計(jì)算出柵格被賦于上述3種類(lèi)型的歸屬概率。

    (4) 同理,對(duì)于25°以下的坡地,推算柵格為林地、灌叢、草地、耕地時(shí)各自覆蓋度,并分別計(jì)算出柵格被賦于上述四種類(lèi)型的歸屬概率。

    (5) 按照不同土地利用類(lèi)型的賦值順序(耕地-草地-灌叢-林地),在適宜蓋度范圍內(nèi)按各柵格的歸屬概率從大到小進(jìn)行屬性分配,直至滿足各土地利用類(lèi)型的總量控制條件為止。

    (6) 將剩余的柵格賦值為沙地和裸砒砂巖,考慮其原有空間位置及土壤侵蝕模數(shù)大小以滿足總量控制所確定的面積。

    (7) 考察優(yōu)化后的土地利用格局所對(duì)應(yīng)的土壤侵蝕模數(shù)和植被蓋度,討論是否滿足生態(tài)安全條件。

    3 結(jié)果及分析

    3.1 土地利用變化

    從1987年到2007年(表2和圖1),流域土地利用發(fā)生了顯著變化:灌叢、林地和建筑用地面積均呈現(xiàn)持續(xù)增長(zhǎng),所占比例分別增長(zhǎng)了11.96、4.93和0.12個(gè)百分點(diǎn),而主要位于流域中部的裸砒砂巖和沙地面積則持續(xù)減少,到2007年所占比例分別下降至13.11%和1.37%,表明經(jīng)過(guò)近20年的水土保持綜合治理,流域生態(tài)恢復(fù)的效果明顯;城鎮(zhèn)化和經(jīng)濟(jì)發(fā)展也推動(dòng)了建筑用地的劇烈擴(kuò)張;水域面積由8.05%減至2.06%,減少了194.08 km2,對(duì)流域未來(lái)植被恢復(fù)以及農(nóng)業(yè)生產(chǎn)等極為不利。此外,由于不合理的土地利用行為不斷加劇,流域西北部草地面積大幅減少,轉(zhuǎn)變?yōu)樯车鼗蛘咴斐膳皫r裸露。從2007年到2015年(表2和圖1),建筑用地和耕地的面積仍將以26.67%和14.89%的速度保持增長(zhǎng),而水域面積則繼續(xù)縮減,進(jìn)一步加深人-水、畜-水矛盾,成為阻礙當(dāng)?shù)亟?jīng)濟(jì)發(fā)展、社會(huì)安定和生態(tài)恢復(fù)的瓶頸性因素;沙地和裸砒砂巖面積存在反彈趨勢(shì),流域生態(tài)環(huán)境可能面臨再次惡化的局面。

    從土地利用類(lèi)型間變化來(lái)看(表3和圖1),建筑用地?cái)U(kuò)張主要是通過(guò)占用灌叢和草地而來(lái),而對(duì)耕地的占用相對(duì)較少,說(shuō)明當(dāng)?shù)刂匾暩乇Wo(hù)的同時(shí)把建筑用地的需求轉(zhuǎn)移到灌叢和草地上面;林地面積增加159.73 km2,灌叢對(duì)其貢獻(xiàn)最大,其次是裸砒砂巖和草地,反映出在流域綜合治理中發(fā)展林果業(yè)成效不錯(cuò),其中灌叢和草地上主要栽種經(jīng)濟(jì)林,裸砒砂巖則以生態(tài)林為主;減少的水域主要轉(zhuǎn)變?yōu)楦亍⒉莸睾吐闩皫r,表明流域降水條件并不理想,同時(shí)大面積的植被恢復(fù)也消耗了大量的水分;而對(duì)裸砒砂巖的治理基本上灌、草措施為主,同時(shí)加強(qiáng)耕作措施的維護(hù)。此外,沙地的治理效果也非常明顯,其中有144.37 km2的沙地被治理成草地。

    表2 皇甫川流域土地利用變化模擬與預(yù)測(cè)

    3.2 土壤侵蝕變化

    1987年、1995年、2000年、2007年和2015年的年平均土壤侵蝕模數(shù)分別為16160.72、14939.42、13943.32、11354.03 t/km2和9710.72 t/km2。在此基礎(chǔ)上,根據(jù)《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)(SLl90—2007)》[19]進(jìn)行土壤侵蝕強(qiáng)度劃分。

    從1987年到2000年(表4和圖2),微度侵蝕和輕度侵蝕的面積分別由1066.15 km2和558.85 km2增加到了1546.16 km2和704.85 km2,中度侵蝕一直到極劇烈侵蝕的面積都呈減少態(tài)勢(shì),流域南部地區(qū)劇烈侵蝕和極劇烈侵蝕的面積顯著減少,均反映出水土流失綜合治理取得良好成效。從2000年到2015年(表4和圖2),微度侵蝕面積持續(xù)減少,輕度侵蝕面積比例在22%上下波動(dòng),雖然兩者面積總和在2015年將保持在2016.58 km2,但流域總體生態(tài)環(huán)境仍比較脆弱:中度、強(qiáng)度和劇烈侵蝕面積均有所回升,三者之和將超過(guò)總面積的四分之一。

    圖1 皇甫川流域1987—2015年土地利用格局Fig.1 Land use patterns from 1987 to 2015 in Huangfuchuan watershed1987、1995、2000和2007分別為當(dāng)年的土地利用現(xiàn)狀格局,2015(F)為2015年土地利用預(yù)測(cè)格局

    土地利用類(lèi)型Landusetype林地Woodland灌叢Shrub草地Grassland耕地Croplan沙地Sand裸砒砂巖Softrock建筑用地Constructionland水域Water林地Woodland0.214540.367010.126430.076700.009660.198670.000370.00662灌叢Shrub0.084140.496980.177100.130970.027960.070690.000860.01130草地Grassland0.034110.282860.478480.042610.044560.101580.000270.01553耕地Cropland0.095450.387170.099510.335730.006530.057260.000580.01777沙地Sand0.023410.338100.391700.031360.173970.037810.001420.00223裸砒砂巖Softrock0.055490.212700.270590.029680.003630.413240.000410.01426建筑用地Constructionland0.018250.087020.079850.015850.007750.082170.695490.01362水域Water0.040720.198240.199640.097720.009510.098940.000480.35475

    表4 皇甫川流域土壤侵蝕變化情況

    圖2 皇甫川流域1987—2015年土壤侵蝕格局Fig.2 Soil erosion patterns from 1987 to 2015 in Huangfuchuan watershed1987、1995、2000和2007分別為當(dāng)年的土壤侵蝕現(xiàn)狀格局,2015(F)為2015年土壤侵蝕預(yù)測(cè)格局

    從不同土地利用類(lèi)型的年平均土壤侵蝕模數(shù)來(lái)看(表5),從1987年到2007年,除沙地的侵蝕模數(shù)增長(zhǎng)80.42%外,其余均呈下降態(tài)勢(shì),其中裸砒砂巖減少了15467.27 t/km2,其次為草地和灌叢,分別減少了1250.44 t/km2和1047.50 t/km2,林地則減少了996.21 t/km2。到2015年,雖然沙地和裸砒砂巖的受蝕程度明顯減弱,但耕地和草地的侵蝕模數(shù)將出現(xiàn)大幅反彈,林地和灌叢也將回復(fù)到2000年的水平。此外, 6種不同土地利用類(lèi)型多年平均土壤侵蝕模數(shù)的大小順序?yàn)椋郝闩皫r>沙地>耕地>草地>林地>灌叢,表明灌、林措施是流域植被恢復(fù)和土壤侵蝕減少的首選,草地減蝕的效果也比較明顯,而沙地尤其是分布面積較廣的裸砒砂巖則是水土流失綜合治理的難題和關(guān)鍵。

    表5 皇甫川流域不同土地利用類(lèi)型年平均土壤侵蝕模數(shù)

    3.3 適宜植被蓋度

    經(jīng)初步計(jì)算,皇甫川流域不同植被類(lèi)型的適宜植被蓋度為[6]:楊樹(shù)37%—48%,油松36%—47%,沙柳46%—59%,檸條43%—55%,沙棘61%—79%,百里香68%—87%,沙打旺50%—65%,羊柴為54%—69%。根據(jù)本研究的具體目標(biāo),選定研究區(qū)4種主要的土地利用類(lèi)型的植被適宜蓋度為:林地42%,灌叢55%,草地77.5%,耕地40%,其中耕地的適宜植被蓋度是通過(guò)與研究區(qū)其它植被類(lèi)型生態(tài)用水的比較來(lái)進(jìn)行估算的[7]。

    3.4 土壤侵蝕情景

    從土壤侵蝕情景模擬結(jié)果來(lái)看(表6),隨著流域植被蓋度的增加(2007年流域植被蓋度約為40%,低于林、灌、草、耕各單一類(lèi)型的適宜植被蓋度),單一土地利用覆蓋類(lèi)型下流域的土壤侵蝕模數(shù)顯著減小,其中草地和灌叢的提升空間較大,主要是因?yàn)槠洮F(xiàn)有蓋度與適宜蓋度之間差異明顯,且面積比例較大的緣故;其次為耕地,而林地通過(guò)增加植被蓋度來(lái)減少土壤侵蝕的效果最差。雖然水分平衡條件下現(xiàn)狀格局比土地利用現(xiàn)狀格局的土壤侵蝕減少了近3900 t km-2a-1,但仍遠(yuǎn)高于流域自然侵蝕臨界值或土壤允許侵蝕量5300 t km-2a-1,表明只提高植被蓋度,而不從數(shù)量和空間位置上對(duì)土地利用進(jìn)行調(diào)整,無(wú)法滿足實(shí)現(xiàn)控制水土流失的生態(tài)安全目標(biāo)。

    3.5 土地利用優(yōu)化

    2007年優(yōu)化前后(表7和圖3),皇甫川流域土地利用格局由“灌—草—裸巖”型轉(zhuǎn)變?yōu)閮?yōu)化格局下的“灌—草—耕”型,裸砒砂巖和沙地面積均減少了45%以上,草地和耕地則呈增加態(tài)勢(shì)。從空間上看,流域中部的沙地基本上被草地所代替,中南部裸砒砂巖則大部分被治理成灌叢;東部灌叢變得相對(duì)破碎,而草地在其間隙中得到了較好的恢復(fù)和擴(kuò)展;耕地則主要是依靠河道向兩側(cè)擴(kuò)張。從效果上看,優(yōu)化后流域土壤侵蝕模數(shù)由11354.03 t km-2a-1減至5884.46 t km-2a-1,考慮到生態(tài)恢復(fù)和重建周期長(zhǎng)、難度大的特點(diǎn),優(yōu)化效果良好,但仍高出流域允許土壤侵蝕模數(shù)。

    2015年優(yōu)化前后(表7和圖3),皇甫川流域土地利用格局將由預(yù)測(cè)格局下的“灌—草—裸巖”型轉(zhuǎn)變?yōu)閮?yōu)化格局下的“灌—草—林”型,裸砒砂巖和沙地面積減幅均在35%左右,林地則有近40%的快速增長(zhǎng)。從空間上看,草地在流域北部和東側(cè)面積有所增加;除在原有中西部地區(qū)有所擴(kuò)張外,城鎮(zhèn)居民點(diǎn)周?chē)_(kāi)始出現(xiàn)分布較為集中的林地;河道兩側(cè)的耕地則繼續(xù)減少。對(duì)比2015年優(yōu)化前后土地利用格局的土壤侵蝕模數(shù)可以發(fā)現(xiàn),在加強(qiáng)林、灌、草植被建設(shè)(包括盡量增加其面積比例和滿足其適宜植被蓋度),適當(dāng)縮減耕地面積,優(yōu)化土地空間布局的基礎(chǔ)上,土壤侵蝕模數(shù)降為4277.95 t km-2a-1,低于流域允許土壤侵蝕模數(shù),表明土地利用格局優(yōu)化明顯改善了生態(tài)環(huán)境狀況,初步實(shí)現(xiàn)了流域生態(tài)安全的目標(biāo)。

    表6 皇甫川流域不同土地利用格局土壤侵蝕模數(shù)

    表7 優(yōu)化前后皇甫川流域土地利用格局及其土壤侵蝕模數(shù)

    4 結(jié)論

    基于土地利用/覆蓋變化與適宜植被蓋度,本文對(duì)皇甫川流域不同土地利用情景下的土壤侵蝕進(jìn)行了定量評(píng)價(jià),進(jìn)而提出包括總量控制和空間配置在內(nèi)的土地利用格局優(yōu)化模式,較好地實(shí)現(xiàn)了以減少土壤侵蝕、增加植被覆蓋為目標(biāo)的區(qū)域生態(tài)安全。主要結(jié)論如下:

    (1) 土地利用結(jié)構(gòu)與植被覆蓋變化是導(dǎo)致流域土壤侵蝕變化的重要因素。從1987年到2015年,隨著建筑用地、林地和灌叢面積持續(xù)增加,水域、沙地和裸砒砂巖面積不斷減少,年平均土壤侵蝕模數(shù)由16160.72 t km-2a-1降至9710.72 t km-2a-1;在水分平衡條件下,若使各主要植被類(lèi)型達(dá)到其適宜植被蓋度,2007年土壤侵蝕模數(shù)將由11354.03 t km-2a-1降至7476.71 t km-2a-1。兩者仍遠(yuǎn)高于流域自然侵蝕臨界值或土壤允許侵蝕量5300 t km-2a-1,表明單純調(diào)整土地利用數(shù)量結(jié)構(gòu)或提高植被覆蓋水平,而不從數(shù)量、空間、覆蓋狀況同步對(duì)土地利用進(jìn)行調(diào)整,不能完全滿足實(shí)現(xiàn)遏制流域水土流失的生態(tài)安全目標(biāo)。

    圖3 皇甫川流域1987年和2015年土地利用和土壤侵蝕優(yōu)化格局Fig.3 Optimized patterns of land use and soil erosion in 1987 and 2015 in Huangfuchuan watershed2007(O)和2015(O)分別為2007年和2015年的土地利用和土壤侵蝕優(yōu)化格局

    (2) 土地利用類(lèi)型與土壤侵蝕密切相關(guān)并存在一定規(guī)律。皇甫川流域6種不同土地利用類(lèi)型多年平均土壤侵蝕模數(shù)的大小順序?yàn)椋郝闩皫r>沙地>耕地>草地>林地>灌叢,表明灌、林措施是流域植被恢復(fù)和減少土壤侵蝕的首選,草地減蝕的效果也比較明顯,而沙地尤其是分布面積較廣的裸砒砂巖則是土壤侵蝕綜合治理的難題和關(guān)鍵。

    (3) 土地利用格局優(yōu)化有助于改善生態(tài)系統(tǒng)功能與服務(wù),實(shí)現(xiàn)流域生態(tài)安全目標(biāo)。在優(yōu)化格局下,2007年和2015年皇甫川流域土壤侵蝕模數(shù)比優(yōu)化前分別減少了5469.57 t km-2a-1和5432.77 t km-2a-1;2015年優(yōu)化后,流域土地利用格局由“灌—草—裸巖”型轉(zhuǎn)變“灌—草—林”型,土壤侵蝕模數(shù)僅為4277.95 t km-2a-1,已經(jīng)實(shí)現(xiàn)了低于流域允許土壤侵蝕模數(shù)的目標(biāo)。

    5 討論

    (1) 在土壤侵蝕模擬過(guò)程中,由于缺乏相關(guān)數(shù)據(jù)和野外實(shí)驗(yàn)的支持,某些參數(shù)的選取采用了前人在黃土高原其它地區(qū)的研究成果,雖然模擬精度仍有待提高,但對(duì)土壤侵蝕的總體趨勢(shì)變化不會(huì)產(chǎn)生大的影響;而在預(yù)測(cè)過(guò)程中,由于無(wú)法獲取流域未來(lái)時(shí)點(diǎn)下的植被蓋度狀況,2015年土壤侵蝕預(yù)測(cè)結(jié)果會(huì)有所偏差。

    (2) 在適宜植被蓋度估算中,文本僅對(duì)單一植被類(lèi)型的適宜蓋度進(jìn)行了模擬,考慮到皇甫川流域存在的植被混合分布情況,下一步應(yīng)該對(duì)由不同植被種類(lèi)組成的植被群落適宜蓋度進(jìn)行分析和模擬。

    (3) 在土地利用格局空間優(yōu)化過(guò)程中,不同土地利用類(lèi)型的分配順序會(huì)對(duì)土地單元空間位置的確定產(chǎn)生影響,進(jìn)而影響到土壤侵蝕模數(shù)的計(jì)算。本文采用的“耕地—草地—灌叢—林地”的分配順序,基本上滿足了流域生態(tài)安全目標(biāo)實(shí)現(xiàn)過(guò)程中的“先生存條件安全、后生態(tài)環(huán)境安全”的要求,但就不同分配順序?qū)?yōu)化過(guò)程所產(chǎn)生的影響及其現(xiàn)實(shí)意義仍可作進(jìn)一步的探討。

    [1] Yu F, Li X B, Wang H, Yu H J.Land use change and eco-security assessment of Huangfuchuan watershed.Acta Geographica Sinca, 2006, 61(6): 645- 653.

    [2] Shi P J, Gong P, Li X B, Chen J, Qi Y, Pan Y Z.The methods and practice of study on the land use and land cover change.Beijing: Science Press, 2000.

    [3] Herrmann S, Osinski E.Planning sustainable land use in rural areas at different spatial using levels GIS and modeling tools.Landscape and Urban Planning, 1999, 46(1/3): 93- 101.

    [4] Makowski D, Hendrix E M T, Van Ittersum M K, Rossing W A H.A framework to study nearly optimal solutions of linear programming models developed for agricultural land use exploration.Ecological Modelling, 2000, 131(1): 65- 77.

    [5] Xu X X, Zhang S B, Wang S Q.Optimizing land use structure by linear programming in the loess hilly region for constructing eco-agriculture.Agricultural Research in the Arid Areas, 2001, 19(2): 94- 99.

    [6] Lin Z P, Liu X N.A case study on land use pattern under ecological security in ecotone between agriculture and animal husbandry in Northeastern China.Chinese Journal of Ecology, 2002, 21(6): 15- 19.

    [7] Liu Y F, Ming D P, Yang J Y.Optimization of land use structure based on ecological green equivalent.Geomatics and Information Science of Wuhan University, 2002, 27(5): 493- 498,515- 515.

    [8] Su W, Chen Y H, Wu Y F, Li J, Zhang J S.Study on land use pattern optimization based on ecosecurity---a case study of agriculture and animal husbandry interlaced zone in Northern China.Progress in Natural Science, 2006, 16(2): 207- 214.

    [9] Lü H F, Pan D H, Liu Z, Wang J A.Optimization of land use structure of desert regions at ecological security level---a case study of typical region in Hulunbeier sandy land.Journal of Hulunbeier College, 2006, 14(4): 6- 8.

    [10] Liang W H, Liao H P, Yang W, Wei H B, Pan Z.Optimization of land use pattern based on ecological security in the mountainous region in Southeast China——a case study of Kaixian County, Chongqing.Journal of Southwest China Normal University (Natural Science Edition), 2012, 37(5): 119- 126.

    [11] Yu F, Li X B, Chen Y H, Wang H, Yang M C.Land use change and soil erosion evaluation in Huangfuchuan watershed.Acta Ecologica Sinica, 2006, 26(6): 1947- 1956.

    [12] Shi P J, Liu B Y, Zhang K L, Jin Z P.Soil erosion process and model studies.Resources Science, 1999, 21(5): 9- 18.

    [13] Jin Z P, Shi P J, Hou F C.System mode and conservation pattern of soil erosion in Huangfuchuan watershed.Beijing: China Ocean Press, 1992.

    [14] You S C, Li W Q.Estimation of soil erosion supported by GIS---a case study in Guanji township, Taihe, Jiangxi.Journal of Natural Resources, 1999, 14(1): 62- 68.

    [15] Williams J R, Renard K G, Dyke P T.EPIC---a new method for assessing erosion′s effect on soil productivity.Journal of Soil and Water Conservation, 1983, 38(6): 381- 383.

    [16] Yang J, Gao Q Z, Li G Q, Jin Z P.A study on the water ecology of some dominant plants in Huangfuchuan Basin.Acta Ecologica Sinica, 2004, 24(11): 2387- 2394.

    [17] Yang J, Song B Y, Piao S J, Tong C, Gao Q Z.Experimental study on ecological use of water of a small catchment in huangfuchuan area.Journal of Natural Resources, 2003, 18(5): 513- 521.

    [18] Chen Y H, Li X B, Shi P J, Gong A D, Dou W.Land cover pattern optimization for local ecological security using remotely sensed data.International Journal of Remote Sensing, 2006, 27(10): 2003- 2010.

    [19] The Ministry of Water Resources of PRC.SLl90—2007 Industry Standard of the People′s Republic of China, Standards for Classification and Gradation of Soil Erosion.Beijing: China Water & Power Press, 2008.

    參考文獻(xiàn):

    [1] 喻鋒, 李曉兵, 王宏, 余弘婧.皇甫川流域土地利用變化與生態(tài)安全評(píng)價(jià).地理學(xué)報(bào), 2006, 61(6): 645- 653.

    [2] 史培軍, 宮鵬, 李曉兵,陳晉,齊曄,潘耀忠.土地利用/覆蓋變化研究的方法與實(shí)踐.北京: 科學(xué)出版社, 2000.

    [5] 徐學(xué)選, 張世彪, 王栓全.黃土丘陵區(qū)生態(tài)建設(shè)中農(nóng)林牧土地結(jié)構(gòu)優(yōu)化模式探討.干旱地區(qū)農(nóng)業(yè)研究, 2001, 19(2): 94- 99.

    [6] 林彰平, 劉湘南.東北農(nóng)牧交錯(cuò)帶土地利用生態(tài)安全模式案例研究.生態(tài)學(xué)雜志, 2002, 21(6): 15- 19.

    [7] 劉艷芳, 明冬萍, 楊建宇.基于生態(tài)綠當(dāng)量的土地利用結(jié)構(gòu)優(yōu)化.武漢大學(xué)學(xué)報(bào): 信息科學(xué)版, 2002, 27(5): 493- 498,515- 515.

    [8] 蘇偉, 陳云浩, 武永峰, 李京, 張瑾水.生態(tài)安全條件下的土地利用格局優(yōu)化模擬研究——以中國(guó)北方農(nóng)牧交錯(cuò)帶為例.自然科學(xué)進(jìn)展, 2006, 16(2): 207- 214.

    [9] 呂紅峰, 潘東華, 劉珍, 王靜愛(ài).基于生態(tài)安全條件的沙區(qū)土地利用結(jié)構(gòu)優(yōu)化——以呼倫貝爾沙地典型樣區(qū)為例.呼倫貝爾學(xué)院學(xué)報(bào), 2006, 14(4): 6- 8.

    [10] 梁偉恒, 廖和平, 楊偉, 魏洪斌, 潘卓.基于生態(tài)安全的西南山地丘陵區(qū)土地利用優(yōu)化模式研究——以重慶市開(kāi)縣為例.西南師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2012, 37(5): 119- 126.

    [11] 喻鋒, 李曉兵, 陳云浩, 王宏, 楊明川.皇甫川流域土地利用變化與土壤侵蝕評(píng)價(jià).生態(tài)學(xué)報(bào), 2006, 26(6): 1947- 1956.

    [12] 史培軍, 劉寶元, 張科利, 金爭(zhēng)平.土壤侵蝕過(guò)程與模型研究.資源科學(xué), 1999, 21(5): 9- 18.

    [13] 金爭(zhēng)平, 史培軍, 候福昌.黃河皇甫川流域土壤侵蝕系統(tǒng)模型和治理模式.北京: 海洋出版社, 1992.

    [14] 游松財(cái), 李文卿.GIS支持下的土壤侵蝕量估算——以江西省泰和縣灌溪鄉(xiāng)為例.自然資源學(xué)報(bào), 1999, 14(1): 62- 68.

    [16] 楊劼, 高清竹, 李國(guó)強(qiáng), 金爭(zhēng)平.皇甫川流域幾種主要植物水分生態(tài)特征.生態(tài)學(xué)報(bào), 2004, 24(11): 2387- 2394.

    [17] 楊劼, 宋炳煜, 樸順姬, 仝川, 高清竹.皇甫川丘陵溝壑區(qū)小流域生態(tài)用水實(shí)驗(yàn)研究.自然資源學(xué)報(bào), 2003, 18(5): 513- 521.

    [19] 中華人民共和國(guó)水利部.SLl90—2007 中華人民共和國(guó)行業(yè)標(biāo)準(zhǔn), 土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn).北京: 中國(guó)水利水電出版社, 2008.

    Optimizationoflandusepatternbasedoneco-securityacasestudyinthehuangfuchuanwatershed

    YU Feng1,2, LI Xiaobing1,*, WANG Hong1

    1StateKeyLaboratoryofEarthSurfaceProcessesandResourceEcology,CollegeofResourcesScienceandTechnology,BeijingNormalUniversity,Beijing100875,China2InformationCenterofMinistryofLandandResources,Beijing100812,China

    In terms of ecological and environmental problems caused by human activities, changes in the pattern of land use play a decisive role in regional eco-security.In recent years, ecological considerations have been an important limiting factor and optimization goal, with the concept of eco-security reflected in land use pattern optimization.Supported by 3S technology, we chose the Huangfuchuan watershed of the Yellow River basin as the study area—a rare, sandy, gritty, intense erosion area of the world, where soil erosion and ecological water are important factors determining regional eco-security.Here we carried out comprehensive research that involved the monitoring and forecasting of change in land use, soil erosion simulation, estimation of appropriate vegetation coverage, and optimization of land use.The results were as follows.(1) Over the last 10 years, there has been a general change in land use in the Huangfuchuan watershed, with the proportion of construction land, woodland, cropland and shrub increasing gradually over time (although the rate of increase has slowed), and the proportion of water, grassland, sand and bare rock decreasing, most markedly in water resources.During this course of development, the watershed pattern of land use has been under tremendous pressure from the conflict between rapid urbanization, economic development and the conservation and rebuilding of the ecological environment.This scenario looks set to continue for some time in the future, and moreover, the distinct decline in water area would be the main obstacle in the above conflict.(2) Change in the pattern of land use was an important factor in the course of soil erosion, the annual modulus of soil erosion decreased rapidly from 16160.72 t/km2in 1987 to 9710.72 t/km2in 2015.To a degree, soil erosion was closely related to type of land use.According to the annual modulus of soil erosion in the same year, different types of land use impacted on soil erosion with the following order of magnitude: bare rock > sand > cropland > grassland > woodland > shrub, indicating that woodland and shrub were the preferred measures for vegetation restoration and a reduction in soil erosion.Furthermore, the alteration of sand and bare rock would be critical to the comprehensive control of soil erosion, a fact that could not be ignored.(3) The effect of soil and water conservation has been significant over the last 10 years; however, owing to the high proportion of sand and bare rock and the proliferating cropland in the watershed, the annual modulus of soil erosion remained higher than the permissible degree of soil erosion, namely 5300 t/km2.Land use pattern optimization can help to improve ecosystem services in the Huangfuchuan watershed, and achieve the goal of promoting eco-security.After optimization, the annual modulus of soil erosion could decrease by 5469.57 t/km2and 5432.77 t/km2in 2007 and 2015, with 4277.95 t/km2in 2015, which is considerably below the permissible degree of soil erosion.Generally, our paper has emphasized the impact of different types and patterns of land use on soil erosion, highlighted the adjustment or optimization in pattern of land use, and suggested improvements in the scientific assessment of soil erosion for the simulation and forecasting of land use change, which has broad future applications.

    eco-security; land use change; soil erosion evaluation; appropriate vegetation coverage; pattern optimization

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目(2014CB138803);國(guó)家自然科學(xué)基金資助項(xiàng)目(41030535,31000229)

    2013- 07- 01;

    2014- 02- 24

    10.5846/stxb201307011807

    *通訊作者Corresponding author.E-mail: xbli@bnu.edu.cn

    喻鋒,李曉兵,王宏.生態(tài)安全條件下土地利用格局優(yōu)化——以皇甫川流域?yàn)槔?生態(tài)學(xué)報(bào),2014,34(12):3198- 3210.

    Yu F, Li X B, Wang H.Optimization of land use pattern based on eco-security: a case study in the huangfuchuan watershed.Acta Ecologica Sinica,2014,34(12):3198- 3210.

    猜你喜歡
    皇甫蓋度模數(shù)
    Absorption spectra and enhanced Kerr nonlinearity in a four-level system
    基于單片機(jī)和模數(shù)化設(shè)計(jì)的低壓側(cè)電壓監(jiān)視與保護(hù)裝置
    能源工程(2021年2期)2021-07-21 08:40:02
    黃山市近30 a植被蓋度時(shí)空變化遙感分析
    模數(shù)化設(shè)計(jì)方法在景觀鋪裝設(shè)計(jì)中的應(yīng)用
    綠色科技(2020年11期)2020-08-01 02:23:58
    黃土高原地區(qū)植被蓋度對(duì)產(chǎn)流產(chǎn)沙的影響
    《柳青在皇甫》《柳青言論集》出版
    趙晶、皇甫舟楠設(shè)計(jì)作品
    基于LID模式的城區(qū)排澇模數(shù)探析
    一種新型的RSA密碼體制模數(shù)分解算法
    坡位與植被蓋度對(duì)楓香造林成效的影響

    超碰97精品在线观看| 国产欧美日韩精品亚洲av| 少妇精品久久久久久久| 免费少妇av软件| 欧美激情 高清一区二区三区| 91精品三级在线观看| 亚洲av日韩在线播放| 69精品国产乱码久久久| 2018国产大陆天天弄谢| 精品少妇内射三级| a 毛片基地| 亚洲一区中文字幕在线| 欧美人与性动交α欧美软件| www日本在线高清视频| 国产精品免费视频内射| 国产一级毛片在线| 一区福利在线观看| 天天操日日干夜夜撸| 日本wwww免费看| 看免费av毛片| 国产日韩欧美视频二区| 欧美午夜高清在线| 久久久久久久久久久久大奶| 香蕉国产在线看| 1024视频免费在线观看| 久久久国产精品麻豆| 久久精品人人爽人人爽视色| 免费观看av网站的网址| 国产精品二区激情视频| 成人影院久久| 免费不卡黄色视频| 视频区图区小说| 久久久久久亚洲精品国产蜜桃av| 免费av中文字幕在线| 日日爽夜夜爽网站| 啦啦啦中文免费视频观看日本| 日韩中文字幕欧美一区二区| 久久亚洲精品不卡| 日本a在线网址| 十分钟在线观看高清视频www| 99国产精品99久久久久| 精品亚洲乱码少妇综合久久| 女性生殖器流出的白浆| 日本vs欧美在线观看视频| 91精品国产国语对白视频| 飞空精品影院首页| 女人久久www免费人成看片| 一区在线观看完整版| 丰满迷人的少妇在线观看| 亚洲欧美精品综合一区二区三区| 久久天堂一区二区三区四区| 又大又爽又粗| 9色porny在线观看| 18禁观看日本| 美女脱内裤让男人舔精品视频| 日本wwww免费看| 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 国产在线视频一区二区| 日本av免费视频播放| 制服人妻中文乱码| 国产视频一区二区在线看| 男人爽女人下面视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 热99久久久久精品小说推荐| 丝袜喷水一区| 男女无遮挡免费网站观看| 国产真人三级小视频在线观看| 日韩电影二区| 亚洲成人国产一区在线观看| 天天躁夜夜躁狠狠躁躁| 女人爽到高潮嗷嗷叫在线视频| 俄罗斯特黄特色一大片| 91av网站免费观看| 国产成人精品在线电影| 男女下面插进去视频免费观看| 久久久久国内视频| 欧美成狂野欧美在线观看| 欧美老熟妇乱子伦牲交| 日韩电影二区| 国产又色又爽无遮挡免| 亚洲精品一二三| 国产成人精品在线电影| 99国产精品一区二区蜜桃av | 我要看黄色一级片免费的| 色老头精品视频在线观看| 久久久久国产精品人妻一区二区| 多毛熟女@视频| 日韩 欧美 亚洲 中文字幕| 啦啦啦免费观看视频1| 巨乳人妻的诱惑在线观看| 精品国产一区二区三区四区第35| 男女之事视频高清在线观看| 日韩大码丰满熟妇| 成年人黄色毛片网站| 最黄视频免费看| 一区在线观看完整版| 精品国产乱码久久久久久小说| 在线永久观看黄色视频| 男女国产视频网站| 成人黄色视频免费在线看| 后天国语完整版免费观看| 曰老女人黄片| 99九九在线精品视频| 亚洲欧美精品综合一区二区三区| 亚洲专区中文字幕在线| 国产精品久久久av美女十八| 欧美黑人欧美精品刺激| 伊人亚洲综合成人网| 一进一出抽搐动态| 老司机深夜福利视频在线观看 | 69精品国产乱码久久久| 亚洲欧美清纯卡通| 他把我摸到了高潮在线观看 | 法律面前人人平等表现在哪些方面 | 成人18禁高潮啪啪吃奶动态图| 久久久精品区二区三区| 19禁男女啪啪无遮挡网站| 操出白浆在线播放| 热99国产精品久久久久久7| av在线老鸭窝| 日韩人妻精品一区2区三区| 国产日韩一区二区三区精品不卡| 18禁观看日本| 一本一本久久a久久精品综合妖精| 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 欧美一级毛片孕妇| 丰满迷人的少妇在线观看| 日韩 亚洲 欧美在线| 久久狼人影院| 欧美黄色片欧美黄色片| 免费看十八禁软件| svipshipincom国产片| 高潮久久久久久久久久久不卡| 一级毛片精品| 真人做人爱边吃奶动态| 大片免费播放器 马上看| 精品少妇一区二区三区视频日本电影| 亚洲少妇的诱惑av| 人人妻人人爽人人添夜夜欢视频| 久久免费观看电影| 18禁裸乳无遮挡动漫免费视频| 黄片播放在线免费| 十八禁高潮呻吟视频| 精品免费久久久久久久清纯 | 男女边摸边吃奶| 亚洲精品国产色婷婷电影| 超碰97精品在线观看| 90打野战视频偷拍视频| 首页视频小说图片口味搜索| 亚洲欧美日韩另类电影网站| 国产一区二区激情短视频 | tube8黄色片| av网站免费在线观看视频| 丰满迷人的少妇在线观看| 成年av动漫网址| 久久中文字幕一级| 伦理电影免费视频| 精品人妻1区二区| 亚洲七黄色美女视频| 欧美少妇被猛烈插入视频| 9热在线视频观看99| 国产日韩欧美在线精品| 黄色视频,在线免费观看| 久久精品亚洲av国产电影网| 色婷婷av一区二区三区视频| 两个人看的免费小视频| 欧美精品啪啪一区二区三区 | 亚洲av电影在线进入| 波多野结衣一区麻豆| 人人澡人人妻人| 另类精品久久| 成年美女黄网站色视频大全免费| 精品国产乱子伦一区二区三区 | 多毛熟女@视频| 美女视频免费永久观看网站| 丝袜美腿诱惑在线| 99热国产这里只有精品6| 国产一区二区 视频在线| 啦啦啦中文免费视频观看日本| 老司机靠b影院| 欧美在线一区亚洲| 操美女的视频在线观看| 日韩免费高清中文字幕av| 无遮挡黄片免费观看| 精品乱码久久久久久99久播| 国产精品av久久久久免费| 日韩大码丰满熟妇| 夜夜夜夜夜久久久久| 国产成人免费无遮挡视频| 999久久久精品免费观看国产| 久久精品国产亚洲av香蕉五月 | 80岁老熟妇乱子伦牲交| 免费女性裸体啪啪无遮挡网站| 伊人亚洲综合成人网| 国产高清videossex| 亚洲三区欧美一区| 午夜福利一区二区在线看| 国产精品一区二区在线观看99| 女警被强在线播放| 一级毛片女人18水好多| 黄网站色视频无遮挡免费观看| 91麻豆精品激情在线观看国产 | 日本wwww免费看| 少妇猛男粗大的猛烈进出视频| 久久久久国产精品人妻一区二区| 1024香蕉在线观看| 亚洲一区二区三区欧美精品| 久久午夜综合久久蜜桃| 亚洲精品成人av观看孕妇| 1024香蕉在线观看| 男女免费视频国产| 一级,二级,三级黄色视频| 中文字幕最新亚洲高清| 精品福利观看| 亚洲精品久久成人aⅴ小说| 精品福利永久在线观看| 久久ye,这里只有精品| 999精品在线视频| 国产成人精品在线电影| 欧美黄色淫秽网站| 曰老女人黄片| 亚洲成av片中文字幕在线观看| 黑丝袜美女国产一区| av线在线观看网站| 日本精品一区二区三区蜜桃| 国产亚洲精品久久久久5区| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 欧美精品一区二区免费开放| 国产精品一区二区免费欧美 | 亚洲av日韩精品久久久久久密| 老鸭窝网址在线观看| 久久热在线av| 欧美成人午夜精品| 黄片小视频在线播放| 久久性视频一级片| 成人18禁高潮啪啪吃奶动态图| 国产高清国产精品国产三级| 黄片小视频在线播放| 国产欧美亚洲国产| 极品人妻少妇av视频| 亚洲成人手机| 在线观看一区二区三区激情| 1024视频免费在线观看| 精品久久久久久电影网| netflix在线观看网站| 美女视频免费永久观看网站| 黄色视频,在线免费观看| 在线av久久热| 麻豆国产av国片精品| 亚洲熟女精品中文字幕| 亚洲精品av麻豆狂野| 日韩欧美国产一区二区入口| 男女免费视频国产| av网站在线播放免费| 色老头精品视频在线观看| 国产淫语在线视频| 亚洲精品久久久久久婷婷小说| 久久午夜综合久久蜜桃| 秋霞在线观看毛片| 日本精品一区二区三区蜜桃| 精品国产一区二区久久| 咕卡用的链子| 午夜免费成人在线视频| 中文字幕另类日韩欧美亚洲嫩草| av网站在线播放免费| 建设人人有责人人尽责人人享有的| 亚洲,欧美精品.| 视频区图区小说| 久久久久网色| 91精品三级在线观看| 久久久国产成人免费| 天天躁日日躁夜夜躁夜夜| 悠悠久久av| 丝袜美腿诱惑在线| 伦理电影免费视频| 韩国高清视频一区二区三区| 亚洲av电影在线进入| www.精华液| 大片免费播放器 马上看| 高清欧美精品videossex| 精品亚洲成国产av| 啪啪无遮挡十八禁网站| 亚洲欧洲精品一区二区精品久久久| 女人被躁到高潮嗷嗷叫费观| 免费观看人在逋| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 久久精品人人爽人人爽视色| 久久久久国产一级毛片高清牌| 啦啦啦中文免费视频观看日本| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 51午夜福利影视在线观看| 亚洲av电影在线观看一区二区三区| 亚洲中文字幕日韩| 大香蕉久久网| 十分钟在线观看高清视频www| 亚洲午夜精品一区,二区,三区| 国产91精品成人一区二区三区 | 亚洲国产av新网站| 国产免费视频播放在线视频| 国产在线观看jvid| 伊人亚洲综合成人网| xxxhd国产人妻xxx| 久久久国产欧美日韩av| 1024视频免费在线观看| 成人影院久久| av线在线观看网站| 最近最新免费中文字幕在线| 男女高潮啪啪啪动态图| 动漫黄色视频在线观看| 一个人免费在线观看的高清视频 | 黑人巨大精品欧美一区二区mp4| 在线精品无人区一区二区三| 成人国语在线视频| 欧美激情久久久久久爽电影 | 欧美成人午夜精品| 五月开心婷婷网| 精品一区二区三区四区五区乱码| 免费日韩欧美在线观看| 日韩视频在线欧美| 黄片大片在线免费观看| 亚洲色图综合在线观看| av在线老鸭窝| 国产黄色免费在线视频| 正在播放国产对白刺激| 妹子高潮喷水视频| 天天添夜夜摸| 丰满人妻熟妇乱又伦精品不卡| 中文字幕精品免费在线观看视频| 国产欧美日韩一区二区三 | 亚洲黑人精品在线| 黄色怎么调成土黄色| 久久久国产欧美日韩av| 国产精品成人在线| 欧美 亚洲 国产 日韩一| 国产精品二区激情视频| 操美女的视频在线观看| 97精品久久久久久久久久精品| 国产成人精品在线电影| 免费av中文字幕在线| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 丰满少妇做爰视频| 99re6热这里在线精品视频| 777米奇影视久久| 丝袜在线中文字幕| 国产av精品麻豆| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 国产精品久久久久久人妻精品电影 | 69av精品久久久久久 | 啦啦啦视频在线资源免费观看| 午夜免费观看性视频| 欧美日韩一级在线毛片| 成年人午夜在线观看视频| 蜜桃在线观看..| 桃花免费在线播放| 国产精品国产三级国产专区5o| 精品国产一区二区三区久久久樱花| 国产高清videossex| 国产野战对白在线观看| 久久久久久久久久久久大奶| 亚洲 欧美一区二区三区| 国产xxxxx性猛交| 精品熟女少妇八av免费久了| 黄色视频在线播放观看不卡| 亚洲精品国产精品久久久不卡| 狂野欧美激情性xxxx| 老司机靠b影院| 最近最新免费中文字幕在线| 搡老熟女国产l中国老女人| 久久国产精品男人的天堂亚洲| 亚洲av欧美aⅴ国产| 精品久久久精品久久久| 一二三四社区在线视频社区8| 一级片'在线观看视频| 亚洲全国av大片| 老鸭窝网址在线观看| 亚洲国产欧美网| 亚洲成av片中文字幕在线观看| 成年女人毛片免费观看观看9 | 交换朋友夫妻互换小说| 久久 成人 亚洲| 1024香蕉在线观看| 成人国语在线视频| 欧美日韩黄片免| a级毛片在线看网站| 人妻人人澡人人爽人人| 91麻豆精品激情在线观看国产 | 亚洲专区字幕在线| 中文字幕人妻丝袜一区二区| 免费久久久久久久精品成人欧美视频| 亚洲全国av大片| 成年人黄色毛片网站| 欧美精品亚洲一区二区| 黄色视频,在线免费观看| 黄色a级毛片大全视频| 韩国精品一区二区三区| av超薄肉色丝袜交足视频| 亚洲欧美一区二区三区久久| 人妻一区二区av| www.熟女人妻精品国产| 免费黄频网站在线观看国产| 两性夫妻黄色片| a 毛片基地| 国内毛片毛片毛片毛片毛片| 91老司机精品| av欧美777| 亚洲免费av在线视频| 婷婷色av中文字幕| 久久人妻福利社区极品人妻图片| 天堂8中文在线网| h视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 91精品国产国语对白视频| 久久天堂一区二区三区四区| 大码成人一级视频| 久久精品国产a三级三级三级| 欧美精品一区二区免费开放| 侵犯人妻中文字幕一二三四区| 满18在线观看网站| 在线观看舔阴道视频| 国产精品九九99| 精品少妇久久久久久888优播| 少妇被粗大的猛进出69影院| 欧美一级毛片孕妇| 久久国产精品影院| 亚洲欧美精品自产自拍| 91精品国产国语对白视频| 欧美黑人欧美精品刺激| 一级毛片电影观看| 少妇精品久久久久久久| 69精品国产乱码久久久| 亚洲va日本ⅴa欧美va伊人久久 | 日韩欧美免费精品| 精品第一国产精品| 巨乳人妻的诱惑在线观看| 国产麻豆69| 久久久久久久国产电影| 亚洲人成77777在线视频| 久久久久久久大尺度免费视频| 美女大奶头黄色视频| 欧美另类一区| 日韩视频在线欧美| 亚洲免费av在线视频| 亚洲成人免费电影在线观看| 天天添夜夜摸| 免费在线观看日本一区| 亚洲va日本ⅴa欧美va伊人久久 | 国产成人欧美| 亚洲精华国产精华精| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 黄片大片在线免费观看| 无限看片的www在线观看| 日韩熟女老妇一区二区性免费视频| 人人妻人人添人人爽欧美一区卜| 巨乳人妻的诱惑在线观看| 桃花免费在线播放| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| a级毛片在线看网站| 女性被躁到高潮视频| 91av网站免费观看| 男男h啪啪无遮挡| 久久精品人人爽人人爽视色| 精品国产乱子伦一区二区三区 | 免费久久久久久久精品成人欧美视频| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 天堂俺去俺来也www色官网| 亚洲激情五月婷婷啪啪| 老司机在亚洲福利影院| 久久人妻熟女aⅴ| 热re99久久国产66热| 国产伦理片在线播放av一区| 久久久久精品国产欧美久久久 | 亚洲国产欧美一区二区综合| 又黄又粗又硬又大视频| 亚洲精品美女久久av网站| 黑人欧美特级aaaaaa片| 1024香蕉在线观看| 国产精品 欧美亚洲| 一个人免费看片子| 国产成人免费观看mmmm| 亚洲伊人色综图| 精品国产一区二区久久| 欧美 日韩 精品 国产| 亚洲色图 男人天堂 中文字幕| 建设人人有责人人尽责人人享有的| 亚洲精品日韩在线中文字幕| 欧美精品亚洲一区二区| 国产福利在线免费观看视频| 高清欧美精品videossex| 肉色欧美久久久久久久蜜桃| 丰满少妇做爰视频| 国产精品av久久久久免费| 一本大道久久a久久精品| 亚洲七黄色美女视频| 久久精品亚洲熟妇少妇任你| 男女之事视频高清在线观看| 热99国产精品久久久久久7| 国产亚洲精品久久久久5区| 天天躁日日躁夜夜躁夜夜| 国产三级黄色录像| 久久久久久久精品精品| 一二三四社区在线视频社区8| 免费在线观看日本一区| 19禁男女啪啪无遮挡网站| 脱女人内裤的视频| 五月开心婷婷网| 亚洲成人手机| 美女脱内裤让男人舔精品视频| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女 | 一区二区三区激情视频| av超薄肉色丝袜交足视频| a 毛片基地| 啪啪无遮挡十八禁网站| 亚洲精品一卡2卡三卡4卡5卡 | 久久ye,这里只有精品| 国产又色又爽无遮挡免| 一级毛片电影观看| 久久精品aⅴ一区二区三区四区| 国产精品 欧美亚洲| 午夜福利一区二区在线看| 高清在线国产一区| 亚洲视频免费观看视频| 正在播放国产对白刺激| 欧美激情 高清一区二区三区| 国产人伦9x9x在线观看| 国产不卡av网站在线观看| 青春草亚洲视频在线观看| 亚洲第一青青草原| 精品少妇久久久久久888优播| 19禁男女啪啪无遮挡网站| 美女福利国产在线| 免费高清在线观看日韩| netflix在线观看网站| 国产成人精品在线电影| 国产亚洲精品一区二区www | 欧美亚洲日本最大视频资源| 一区在线观看完整版| 丝瓜视频免费看黄片| 亚洲国产av影院在线观看| 新久久久久国产一级毛片| 中文字幕色久视频| 久久精品国产a三级三级三级| 亚洲欧美成人综合另类久久久| 韩国精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 亚洲欧美色中文字幕在线| 成人三级做爰电影| 精品人妻在线不人妻| 欧美少妇被猛烈插入视频| 欧美激情 高清一区二区三区| 国产片内射在线| 久久亚洲国产成人精品v| 青青草视频在线视频观看| 亚洲精品国产精品久久久不卡| 精品国产国语对白av| 久久中文字幕一级| 99国产精品一区二区三区| 久久免费观看电影| 亚洲精品在线美女| 免费在线观看黄色视频的| 国产欧美日韩一区二区三区在线| 黄色片一级片一级黄色片| 无遮挡黄片免费观看| 亚洲精品国产一区二区精华液| 亚洲成人国产一区在线观看| 精品少妇久久久久久888优播| 捣出白浆h1v1| 王馨瑶露胸无遮挡在线观看| 精品国产一区二区三区久久久樱花| 午夜日韩欧美国产| 在线看a的网站| 一级黄色大片毛片| 麻豆国产av国片精品| 一级黄色大片毛片| 欧美日韩视频精品一区| 自线自在国产av| 成年女人毛片免费观看观看9 | 亚洲av欧美aⅴ国产| 日韩一区二区三区影片| 成年人黄色毛片网站| 亚洲全国av大片| 一本—道久久a久久精品蜜桃钙片| av在线老鸭窝| 男女国产视频网站| 日韩三级视频一区二区三区| av电影中文网址| 久久免费观看电影| 老司机靠b影院| e午夜精品久久久久久久| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区三区| 久久久国产精品麻豆| 免费日韩欧美在线观看| 免费高清在线观看视频在线观看| 国产一级毛片在线| 欧美黄色淫秽网站| 他把我摸到了高潮在线观看 | 国产亚洲精品一区二区www | 丝袜喷水一区| www.熟女人妻精品国产| 国产在线免费精品| 欧美国产精品va在线观看不卡|