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

    一類非線性發(fā)生率植被恢復(fù)模型的穩(wěn)定性分析

    2017-05-07 02:59:46紀(jì)德龍李硯趙新俊
    關(guān)鍵詞:覆蓋面積覆蓋度平衡點(diǎn)

    紀(jì)德龍,李硯,趙新俊

    (石河子大學(xué)理學(xué)院,新疆 石河子832003)

    植被是生態(tài)系統(tǒng)中最基本也是最重要的元素,植被恢復(fù)不僅可以修復(fù)地表和地下的生態(tài)系統(tǒng)[1],還對(duì)土壤理化性質(zhì)[2]有極其重要影響。沙漠邊緣區(qū)的環(huán)境惡劣,依靠人工保護(hù)措施是植被生態(tài)系統(tǒng)恢復(fù)的重要手段之一。目前,植被生態(tài)系統(tǒng)恢復(fù)的方式主要以人工保護(hù)措施為主,其中主要有圍欄封育、人工植被重建以及改善植被退化的影響因子等方式。圍欄封育可以改善土壤成分[3]和植被特征[4-5],但過(guò)長(zhǎng)的圍封年限會(huì)導(dǎo)致植被結(jié)構(gòu)和功能的退化[6]。影響人工植被重建的因素有很多,其中以植被蒸騰蒸發(fā)和土壤特性最為突出。通過(guò)對(duì)植被蒸騰蒸發(fā)[7]、土壤呼吸特性[8]的實(shí)驗(yàn)研究以及采用植被指標(biāo)對(duì)退化草地植被效果進(jìn)行綜合評(píng)價(jià)[9],為人工植被合理重建提供了一定的理論依據(jù)。在荒漠區(qū),影響植被生長(zhǎng)的關(guān)鍵因子地失調(diào)會(huì)導(dǎo)致植被生態(tài)系統(tǒng)的失衡,所以改善荒漠化驅(qū)動(dòng)因子是植被生態(tài)系統(tǒng)恢復(fù)的根本手段。在遏制植被退化過(guò)程中,考慮空間尺度效應(yīng)下的驅(qū)動(dòng)因子[10]、沙漠化驅(qū)動(dòng)因子的分類[11]及其動(dòng)力學(xué)模型[12],對(duì)遏制植被退化更具有針對(duì)性。隨著現(xiàn)代科技的發(fā)展,人們開(kāi)始利用RS、GIS技術(shù)和激光雷達(dá)圖像研究大規(guī)模植被重建下的時(shí)空演變和植被退化特征[13-14],這對(duì)植被恢復(fù)具有重要的意義。

    以上文獻(xiàn)主要圍繞人工保護(hù)措施的方式及其影響因素進(jìn)行研究,未能有效預(yù)測(cè)植被生態(tài)系統(tǒng)恢復(fù)狀況以及恢復(fù)所用時(shí)間。荒漠區(qū)資源匱乏這一特性對(duì)植被人工修復(fù)有很大的遏制作用,從而對(duì)人工保護(hù)措施力度也提出了苛刻的要求?;诖?,本文通過(guò)建立植被恢復(fù)模型定性分析人工保護(hù)措施力度對(duì)植被恢復(fù)的影響,并利用數(shù)值模擬驗(yàn)證所建模型的適用性,為失衡的植被生態(tài)系統(tǒng)更平穩(wěn)的恢復(fù)平衡提供參考依據(jù)。

    1 模型的建立

    由于植被蓋度大小與侵蝕大小是相互驅(qū)動(dòng)的,故本文基于非線性發(fā)生率的SIS傳染病模型,構(gòu)造加入人工保護(hù)系數(shù)的植被恢復(fù)模型,并對(duì)其平衡點(diǎn)的穩(wěn)定性進(jìn)行分析。傳染病模型的研究核心是對(duì)發(fā)生率的刻畫(huà),并且在經(jīng)典的傳染病傳播模型中,還常常設(shè)定發(fā)生率是雙線性的。Moghadas S M討論了DS(Differential Susceptibility)和DI(Differential Infectivity)模型的標(biāo)準(zhǔn)發(fā)生率,并分析了2種模型的穩(wěn)定性[17]。本文在Moghadas S M提出的DS和DI模型的標(biāo)準(zhǔn)發(fā)生率基礎(chǔ)之上,針對(duì)植被被侵蝕后,人工保護(hù)措施對(duì)侵蝕速率 (即植被覆蓋面積的動(dòng)態(tài)變化)的影響,研究一類具有非線性發(fā)生率的被侵蝕植被恢復(fù)模型。通過(guò)加入人工保護(hù)措施因素預(yù)測(cè)植被恢復(fù)的動(dòng)態(tài)變化,對(duì)資源匱乏的沙漠邊緣區(qū)的植被恢復(fù)具有深刻意義。

    基于SIS傳染病模型[18-19],構(gòu)造加入人工保護(hù)系數(shù)的模型,其具體系統(tǒng)模型為:

    其中:S表示t時(shí)刻健康且脆弱(即易被侵蝕)的植被覆蓋面積;I表示t時(shí)刻被侵蝕植被的覆蓋面積;r為內(nèi)稟自然增長(zhǎng)率;K為環(huán)境容納量;β為侵蝕率;α為植被被侵蝕后,人工保護(hù)的作用系數(shù);θ為自然恢復(fù)率;d為因植被被侵蝕后生長(zhǎng)環(huán)境改變所造成的死亡率。系統(tǒng)中所有的參數(shù)均為正數(shù)。

    為了證明和分析的簡(jiǎn)便,對(duì)系統(tǒng)模型(1)作變換dt=K(1+αS)dτ,則模型(1)變?yōu)?/p>

    基于系統(tǒng)模型(1)、(2)的實(shí)際生態(tài)學(xué)意義,本文只在區(qū)域G={(S,I)|S≥0,I≥0}內(nèi)討論問(wèn)題。

    2 平衡點(diǎn)的存在性

    引理1:模型(2)在G內(nèi)總存在無(wú)侵蝕平衡點(diǎn)E0(0,0)和E1(K,0),當(dāng)θα<β<ad-γ/K時(shí),模型(2)在G內(nèi)存在唯一的侵蝕平衡點(diǎn)E2(S*,I*),即植被生態(tài)系統(tǒng)始終被侵蝕,但植被生態(tài)系統(tǒng)不會(huì)崩潰,其中

    證明:分別令模型(2)中兩式為0,即

    可求得:

    故模型平衡點(diǎn)分別為E0(0,0)、E1(K,0)和E2(S*,I*)。

    引理1的生態(tài)學(xué)意義是:E0是無(wú)效的,即現(xiàn)實(shí)中不存在這中情況或?qū)е孪到y(tǒng)發(fā)散,E1和E2與現(xiàn)實(shí)中的植被被侵蝕的狀況相一致。

    3 平衡點(diǎn)的穩(wěn)定性分析

    下面分析模型及其平衡點(diǎn)的穩(wěn)定性[20-21]。由于模型(2)是模型(1)轉(zhuǎn)換而來(lái)的,所以這2個(gè)模型具有相同的平衡點(diǎn)和系統(tǒng)穩(wěn)定性。下面對(duì)模型(2)進(jìn)行分析。其中,定理2-5都是在定理1的基礎(chǔ)上展開(kāi)分析證明的,定理2和定理3分別分析證明了平衡點(diǎn)E0(0,0)和E1(K,0)的穩(wěn)定性,定理4和定理5分別證明了平衡點(diǎn)E2(S*,I*)的不同穩(wěn)定性,二者屬于一種遞進(jìn)關(guān)系。

    定理1:模型(2)的一切正解一定有界。

    證明:由模型(1),令V=S+I,則

    所以∨ε>0,T>0,當(dāng)t>T時(shí),

    因此,模型(1)的一切正解一定有界。

    故模型(2)的一切正解一定有界,證畢。

    定理2:無(wú)侵蝕平衡點(diǎn)E0(0,0)是不穩(wěn)定的。

    證明:模型(2)在E0(0,0)處的Jacobi矩陣為:

    則系統(tǒng)(2)所對(duì)應(yīng)的特征方程為

    顯然特征方程有2個(gè)正實(shí)根,故E0(0,0)是不穩(wěn)定的。

    定理3:當(dāng)β<θα?xí)r,無(wú)侵蝕平衡點(diǎn)E1(K,0)不僅是局部漸近穩(wěn)定的,還在區(qū)域G內(nèi)是全局漸近穩(wěn)定的。

    證明:(i)首先證明平衡點(diǎn)E1(K,0)的局部漸近穩(wěn)定。模型(2)在E1(K,0)處的Jacobi矩陣為:

    則模型(2)所對(duì)應(yīng)的特征方程為

    特征根λ1=-rK(1+αK)根為負(fù)實(shí)根;

    又因?yàn)棣拢鸡圈?,從而?<0,故無(wú)侵蝕平衡點(diǎn)E1(K,0)是局部漸近穩(wěn)定的。

    (ii)再證明平衡點(diǎn)E1(K,0)在G內(nèi)全局漸近穩(wěn)定。事實(shí)上,當(dāng)β<θα?xí)r,模型(2)在區(qū)域G內(nèi)存在唯一的無(wú)侵蝕平衡點(diǎn)E1(K,0),而且平衡點(diǎn)處在邊界上,故在區(qū)域G內(nèi)不存在閉軌線,且模型(2)從G內(nèi)出發(fā)的軌線均會(huì)一直停留在G內(nèi)。又由于區(qū)域G的有界性,對(duì)任給的一個(gè)初始條件,模型(2)滿足初值的解(S,J)最終都會(huì)趨于平衡點(diǎn)E1(K,0),再由平衡點(diǎn)E1(K,0)的局部漸近穩(wěn)定性,知可E1(K,0)是全局漸近穩(wěn)定的。

    由上面定理推導(dǎo)可知:在所研究區(qū)域內(nèi),無(wú)論初始被侵蝕植被生態(tài)系統(tǒng)有多么脆弱,最終被侵蝕植被生態(tài)系統(tǒng)都會(huì)轉(zhuǎn)化為健康的植被生態(tài)系統(tǒng)。

    定理4:當(dāng)θα<β<αd-r/K時(shí),被侵蝕平衡點(diǎn)E2(S*,I*)是局部漸近穩(wěn)定的。

    證明:模型(2)在E2(S*,I*)處的Jacobi矩陣為:

    其中:

    再由Hurwitz判別法可知:

    當(dāng)θα<β<αd-r/K時(shí),侵蝕平衡點(diǎn)E2(S*,I*)是局部漸近穩(wěn)定的。

    從上述定理可知,模型(2)在侵蝕平衡點(diǎn)E2(S*,I*)的鄰域內(nèi)是吸引的。若初始點(diǎn)(S0,I0)在其鄰域內(nèi),則模型(2)所描述的系統(tǒng)是穩(wěn)定的,即隨著時(shí)間推移,植被生態(tài)系統(tǒng)會(huì)達(dá)到一種平衡狀態(tài);若初始點(diǎn)不在(S0,I0)其鄰域內(nèi),模型(2)就有可能是不穩(wěn)定的,即隨著時(shí)間推移,植被生態(tài)系統(tǒng)就有可能趨于崩潰狀態(tài)。

    定理 5:當(dāng)θα<β<αd-r/K且Kα-1<0時(shí),侵蝕平衡點(diǎn)E2(S*,I*)在G內(nèi)是全局漸近穩(wěn)定的。

    根據(jù)定理4可知:一旦植被被侵蝕,侵蝕狀態(tài)就會(huì)一直存在,最終的植被易被侵蝕覆蓋面積S和被侵蝕覆蓋面積面積I將會(huì)穩(wěn)定在S*和I*。

    4 數(shù)值模擬及分析

    4.1 數(shù)據(jù)來(lái)源

    模型中涉及的參數(shù)是依據(jù)古爾班通古特沙漠南緣地區(qū)植被覆蓋面積變化以及植被各等級(jí)度之間實(shí)際覆蓋面積轉(zhuǎn)移的百分比[22]賦值的,這使得參數(shù)更具有實(shí)際效果,從而更加準(zhǔn)確地驗(yàn)證模型的適用性。

    表1是從2001年開(kāi)始到2005年,由于大型工程和居民活動(dòng)導(dǎo)致植被生態(tài)系統(tǒng)失衡后,該地區(qū)植被覆蓋面積的變化。

    表2是2004年國(guó)家對(duì)古爾班通古特沙漠南緣地區(qū)實(shí)施人工保護(hù)措施后,2005-2008年植被覆蓋面積的變化。

    表3和表4分別是表1和表2的覆蓋度等級(jí)之間的轉(zhuǎn)移矩陣,其中植被覆蓋度等級(jí)如下:a為覆蓋度 0-0.1;b為覆蓋度 0.1-0.2;c為覆蓋度 0.2-0.3;d為覆蓋度0.3-0.4;e為覆蓋度>0.4。

    表1 2001-2005年古爾班通古特沙漠南緣植被覆蓋面積變化 hm2Tab.1 Vegetation coverage area variation of the south rim of Gurbantunggut from 2001 to 2005

    表2 2005-2008年古爾班通古特沙漠南緣植被覆蓋面積變化 hm2Tab.2 Vegetation coverage area variation of the south rim of Gurbantunggut from 2005 to 2008

    表3 2001-2005年植被覆蓋度之間的轉(zhuǎn)移矩陣Tab.3 Transition matrix between vegetation coverage from 2001 to 2005 %

    表4 2005-2008年植被覆蓋度之間的轉(zhuǎn)移矩陣Tab.4 Transition matrix between vegetation coverage from 2005 to 2008 %

    依據(jù)表1和表2植被覆蓋度面積的變化情況,確定未被且易被侵蝕的植被覆蓋面積S為植被覆蓋度等級(jí)c、d的實(shí)際面積之和,被侵蝕植被覆蓋面積I為植被覆蓋度等級(jí)a、b的實(shí)際面積之和。

    4.2 參數(shù)的確定

    4.2.1 環(huán)境容納量的確定

    從表1可以看出,2001年覆蓋度等級(jí)為c的植被為82305 hm2,占所研究區(qū)域 90000 hm2的91.45%。

    由以上數(shù)據(jù)可知,植被幾乎未被侵蝕,植被生態(tài)系統(tǒng)正處于平衡狀態(tài)。故此時(shí)的植被覆蓋面積可近似看作已經(jīng)達(dá)到環(huán)境承載量,這具有一定的合理性。故環(huán)境承載量等于各等級(jí)植被實(shí)際蓋度面積之和,即K0=22400。由于作圖時(shí),取值不宜過(guò)大,故取K=22.4。

    4.2.2 內(nèi)稟增長(zhǎng)率r的確定

    在2001-2005年間,由于古爾班通古特沙漠南緣受到居民活動(dòng)和大型工程的影響,其自然內(nèi)稟增長(zhǎng)率也受到很大影響。自2004年起,國(guó)家在此地區(qū)開(kāi)始實(shí)施人工保護(hù)措施后,表2中圍欄內(nèi)覆蓋度等級(jí)為c、d的植被覆蓋面積不僅基數(shù)較大,而且增長(zhǎng)相對(duì)穩(wěn)定,故其轉(zhuǎn)入轉(zhuǎn)出的綜合情況所體現(xiàn)的自然內(nèi)稟增長(zhǎng)率更加接近實(shí)際(由于e的等級(jí)范圍過(guò)大,選取=0.5)。經(jīng)分析,內(nèi)稟增長(zhǎng)率r由下列公式確定:

    代入數(shù)據(jù),可求得r=0.2470。其中:λi→j表示覆蓋等級(jí)度i向覆蓋等級(jí)度j轉(zhuǎn)移的百分比,i(i=a,b,c,d,e)表示各覆蓋等級(jí)度的平均值。

    4.2.3 植被恢復(fù)率θ,被侵蝕率β,死亡率d、初值S0、I0和保護(hù)系數(shù)α的確定

    (1)通過(guò)表1和表2的覆蓋面積變化,以及表3和表4的植被等級(jí)度之間的面積轉(zhuǎn)移百分比可知:對(duì)于植被覆蓋等極度由b→c正向轉(zhuǎn)化73.74%,故由可得,θ=0.4910。

    (2)被侵蝕率β由c→b和d→b轉(zhuǎn)移的綜合情況確定,從而可得下列公式:

    采用表2和表4數(shù)據(jù),由上式解得β=0.0036。

    (3)對(duì)死亡率d的確定,通過(guò)覆蓋度等級(jí)為a,b的植被面積自身增長(zhǎng)率和r之間的差距來(lái)確定。結(jié)合4.2.2可知,d=0.0245。

    (4)對(duì)于初值的確定,通過(guò)表2數(shù)據(jù)并計(jì)算,可知2005年植被覆蓋度等級(jí)為c、d實(shí)際覆蓋面積之和為11000 hm2,植被覆蓋度等級(jí)為a、b的實(shí)際覆蓋面積之和為 7000 hm2,故令S0=11,I0=7。

    (5)通過(guò)模型1的結(jié)構(gòu)以及所確定的初值,并結(jié)合表3和表4的對(duì)比,可令人工保護(hù)系數(shù)α=0.8。

    4.3 模型求解

    將K=22.4,S0=11,I0=7,r=0.2470,θ=0.4910,β=0.0036,d=0.0245,α=0.8 代入模型(1),并利用軟件MATLAB7.1得到圖1。

    圖1 模型(1)的數(shù)值模擬圖Fig.1 Numerical simulation chart of model(1)

    由圖1可知:

    (1)經(jīng)過(guò)數(shù)年后,植被將逐漸趨于無(wú)侵蝕平衡點(diǎn)E1(K,0),即植被生態(tài)系統(tǒng)無(wú)侵蝕。

    (2)通過(guò)參數(shù)的選取可知β<θα,該圖符合定理3的結(jié)論。

    4.4 補(bǔ)充性檢驗(yàn)

    由于模型 (1)中所涉及參數(shù)均是由2001年至2008年古爾班通古特沙漠南緣植被覆蓋面積及覆蓋度等級(jí)面積轉(zhuǎn)移程度得到的,故圖1的曲線走勢(shì)與此地區(qū)2001年至2008年的植被覆蓋面積有很大關(guān)聯(lián),圖1中2005-2008年的植被覆蓋面積并不能很好的驗(yàn)證模型的可靠性。為了充分檢驗(yàn)?zāi)P?1)的預(yù)測(cè)功能,本文選取2010-2013年古爾班通古特沙漠南緣植被覆蓋面積實(shí)際值[23-24](通過(guò)文獻(xiàn)中數(shù)據(jù)及數(shù)學(xué)反演得到),并通過(guò)模型(1)計(jì)算得出2010-2013年此地區(qū)植被覆蓋面積預(yù)測(cè)值。具體數(shù)據(jù)如下:

    表5 近年古爾班通古特沙漠南緣植被覆蓋面積Tab.5 Vegetation coverage area of the south rim of Gurbantunggut from 2010 to 2013 hm2

    利用軟件SPSS對(duì)上述數(shù)據(jù)進(jìn)行配對(duì)t檢驗(yàn)分析,得到下列分析結(jié)果:

    表6 配對(duì)t檢驗(yàn)分析結(jié)果Tab 6 Pairingt-testing results

    由表6可知:Sig.=0.303>0.05,這說(shuō)明這兩組數(shù)據(jù)顯著性差異很小。由此可知,模型(1)的預(yù)測(cè)功能良好。

    5 結(jié)論

    由模型分析及數(shù)值模擬,并結(jié)合生態(tài)學(xué)原理,可得出以下結(jié)論:

    (1)本文通過(guò)微分方程幾何理論和生態(tài)學(xué)原理建立了植被恢復(fù)模型,該模型能有效預(yù)測(cè)植被恢復(fù)狀況,以便當(dāng)植被恢復(fù)狀況異常時(shí),能夠及時(shí)查明原因并采取應(yīng)對(duì)措施。

    (2)在植被恢復(fù)過(guò)程中,通過(guò)調(diào)控人工保護(hù)措施力度,即當(dāng)α≈β/θ時(shí),可以更好地節(jié)省植被恢復(fù)所需資源,這對(duì)資源匱乏的沙漠邊緣區(qū)的植被恢復(fù)具有重要意義。

    (3)當(dāng)β<θα?xí)r,根據(jù)植被生態(tài)系統(tǒng)所具有的特性,隨著時(shí)間推移,被侵蝕的植被生態(tài)系統(tǒng)將會(huì)逐漸恢復(fù),毗鄰已恢復(fù)植被生態(tài)系統(tǒng)地區(qū)的完全沙漠化地帶(植物極其稀少或不存在)的植被數(shù)量與種類也會(huì)趨向于增多的態(tài)勢(shì)。這表明:植被覆蓋面積的增大,生態(tài)系統(tǒng)的不斷修復(fù),土壤中微生物數(shù)量及種類也會(huì)不斷增多,進(jìn)而改善土壤環(huán)境,促進(jìn)植被生長(zhǎng),形成一種良性循環(huán)。

    (4)當(dāng)θα<β<αd-r/K且Kα-1<0時(shí),沙漠化面積會(huì)持續(xù)增大,侵蝕區(qū)域不斷向外拓展。從模型平衡點(diǎn)分析上看,沙漠化侵蝕持續(xù)存在,其侵蝕速度較為溫和,短時(shí)期內(nèi)生態(tài)系統(tǒng)不會(huì)出現(xiàn)崩潰狀態(tài)。這表明:隨著沙漠化面積越來(lái)越大,沙漠化侵蝕速度也會(huì)越來(lái)越快,這就形成了一種惡性循環(huán)。因此,應(yīng)加大保護(hù)力度,以改變參數(shù)之間的關(guān)系,以遏制沙漠化侵蝕。

    6 討論

    (1)大多數(shù)文獻(xiàn)通過(guò)植被生態(tài)因子的變化研究植被修復(fù)問(wèn)題,由于實(shí)驗(yàn)方法的不同和生態(tài)因子的側(cè)重點(diǎn)不同,其研究結(jié)果帶有一定的主觀性。本文所建立的模型,利用植被的種群特征及其數(shù)量關(guān)系,從宏觀層面上研究植被修復(fù)問(wèn)題,這使得研究結(jié)果更傾向于客觀化。

    (2)模型所涉及的參數(shù)由植被受保護(hù)前后若干年份的覆蓋面積變化來(lái)確定,所以,實(shí)施人工保護(hù)措施后,并不能在短時(shí)間內(nèi)預(yù)測(cè)未來(lái)植被恢復(fù)狀態(tài)。另外,此類模型并未考慮自然降水的異常改變,一旦降水量變化過(guò)大,模型所涉及的參數(shù)就會(huì)出現(xiàn)差異,這對(duì)預(yù)測(cè)結(jié)果有一定影響。

    (3)由于人工保護(hù)系數(shù)的改變,使得植被覆蓋面積的內(nèi)稟自然增長(zhǎng)率、恢復(fù)率、非自然死亡率以及被侵蝕率產(chǎn)生變化,從而對(duì)植被生態(tài)系統(tǒng)恢復(fù)時(shí)間的預(yù)測(cè)也會(huì)產(chǎn)生較大的影響。因此,有待對(duì)模型(1)中相關(guān)參數(shù)之間的相互影響進(jìn)行更深入的研究。

    參考文獻(xiàn):

    [1] 胡嬋娟,郭雷.植被恢復(fù)的生態(tài)效應(yīng)研究進(jìn)展[J].生態(tài)環(huán)境學(xué)報(bào),2012,21(9):1640-1646.Hu C J,Guo L.Advances in the research of ecological effects of vegetation restoration[J].Ecology and Environmental Sciences,2012,21(9):1640-1646.

    [2] 付標(biāo),齊雁冰,常慶瑞.不同植被重建管理方式對(duì)沙質(zhì)草地土壤及植被性質(zhì)的影響[J].草地學(xué)報(bào),2015,23(1):47-54.Fu B,Qi Y B,Chang Q R.Effect of different revegetation management methods on the soil and vegetation characteristics of degraded sandy grassland[J].Acta Agrestia Since,2015,23(1):47-54.

    [3] 楊新國(guó),宋乃平,李學(xué)斌,等.短期圍欄封育對(duì)荒漠草原沙化灰鈣土有機(jī)碳組分及物理穩(wěn)定性影響 [J].應(yīng)用生態(tài)學(xué)報(bào),2012,23(12):3325-3330.Yang X G,Song N P,Li X B,et al.Effects of short-term fencing on organic carbon fractions and physical stability of sandy sierozem in desert steppe of Northwest China[J].Chinese Journal of Applied Ecology,2012,23(12):3325-3330.

    [4] 安耕,王天河.圍欄封育改良荒漠化草地的效果[J].草業(yè)科學(xué),2011,28(5):874-876.An G,Wang T H.Effect of enclosure on restoration of degraded grassland[J].Pratacultural Science,2011,23(12):3325-3330.

    [5] 劉雪明,聶學(xué)敏.圍欄封育對(duì)高寒草地植被數(shù)量特征的影響[J].草業(yè)科學(xué),2012,29(1):112-116.Liu X M,Nie X M.Effects of enclosure on the quantitativecharacteristics of alpine vegetation[J].Pratacultural Science,2012,29(1):112-116.

    [6] 劉健,張克斌,程中秋,等.圍欄封育對(duì)沙化草地植被及土壤特性的影響[J].水土保持通報(bào),2011,31(4):180-184.Liu J,Zhang K B,Cheng Z Q,et al.Influences of fencing on vegetation and soil properties in sa ndy grassland[J].Bulletin of Soil and Water Conservation,2011,31(4):180-184.

    [7] 張志山,李新榮,王新平,等.沙漠人工植被區(qū)的蒸發(fā)蒸騰[J].生態(tài)學(xué)報(bào),2005,25(10):2484-2490.Zhang Z S,Li X R,Wang X P,et al.Evaporation and transpiration in revegetated desert area[J].Acta Ecologica Sinica,2005,25(10):2484-2490.

    [8] 高艷紅,劉立超,賈榮亮,等.沙坡頭人工植被演替過(guò)程的土壤呼吸特征[J].生態(tài)學(xué)報(bào),2012,32(8):2474-2482.Gao Y H,Liu L C,Jia R L,et al.Soil respiration patterns during restoration of vegetation in the Shapotou area Northern China[J].Acta Ecologica Sinica,2012,32(8):2474-2482.

    [9] 王雪佼,靳瑰麗,朱進(jìn)忠,等.放牧型退化草地植被重建效果的初步研究[J].新疆農(nóng)業(yè)科學(xué),2013,50(10):1893-1901.Wang X J,Jin G L,Zhu J Z,et al.Preliminary study on the effect of reconstructed grassland for grazing[J].Xinjiang Agricultural Sciences,2013,50(10):1893-1901.

    [10] 鄭偉,朱進(jìn)忠.新疆草地荒漠化過(guò)程及驅(qū)動(dòng)因素分析[J].草業(yè)科學(xué),2012,29(9):1340-1351.Zheng W,Zhu J Z.Analysis of desertification process and driving force factors in grassland ecosystem of Xinjiang[J].Pratacultural Science,2012,29(9):1340-1351.

    [11] 彭羽,米凱,卿鳳婷,等.影響植被退化生態(tài)因子的多尺度分析[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2015,23(增刊):11-19.Peng Y,Mi K,Qing F T,et al.Spatial scale effects of ecological factors on vegetation degradation[J].Journal of Basic Science and Engineering,2015,23(Supplement):11-19.

    [12] Jiao J,Cai S H,Li L M.Analysis of the Oasis Vegetation Degradation Dynamical model with Impulsive Irrigation in Arid Area[J].Mathematica Applicata,2016,29(2):409-417.

    [13] 秦偉,朱清科,左長(zhǎng)清,等.大規(guī)模植被重建背景下的黃土高原流域土地利用時(shí)空演變 [J].水土保持學(xué)報(bào),2014,28(5):43-50.Qin W,Zhu Q K,et al.Spatio-temporal changes of Land use/cover driven by large-scale vegetation restoration on watershed in the loess plateau[J].Journal of Soil and Water Conservation,2014,28(5):43-50.

    [14] 肖靜.基于激光雷達(dá)圖像分析的植被退化特征提取[J].激光雜志,2009,37(7):50-53.Xiao J.Vegetation degradation feature extraction based on the analysis of the laser radar image[J].Laser Journal,2009,37(7):50-53.

    [15] 杜艷可,徐瑞,段立江.一類具有標(biāo)準(zhǔn)發(fā)生率的SIS型傳染病模型的全局穩(wěn)定性[J].數(shù)學(xué)的實(shí)踐與認(rèn)知,2009,39(10):140-144.Du Y K,Xu R,Duan L J.Global stability of an SIS epidemic model with standard incidence Rate[J].Mathematics in Practice and Theory,2009,39(10):140-144.

    [16]朱慧強(qiáng).具有非線性發(fā)生率和非定常人口的傳染病傳播模型分析[D].長(zhǎng)沙:湖南大學(xué),2012,39-48.

    [17]Moghadas SM.Two core group models for sexual transmission of disease[J].Ecological modeling,2002,148(1):15-26.

    [18] 陳蘭蓀.數(shù)學(xué)生態(tài)學(xué)模型與研究方法[M].北京:科學(xué)出版社,1988:11-34.

    [19] 馬知恩,周義倉(cāng),王穩(wěn)定.傳染病動(dòng)力學(xué)的數(shù)學(xué)建模與研究[M].北京:科學(xué)出版社,2004:239-242.

    [20] 馬知恩,周義倉(cāng).常微分方程定性與穩(wěn)定性方法[M].北京:科學(xué)出版社,2001:59-64.

    [21] 張錦炎,馮貝葉.常微分方程幾何理論與分支問(wèn)題[M].北京:北京大學(xué)出版社,2000:31-37.

    [22]劉乾.古爾班通古特沙漠南緣植物群落多樣性與植被覆蓋變化分析[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2009:48-52.

    [23] 許玉鳳,楊井,陳亞寧,等.近32年來(lái)新疆地區(qū)植被覆蓋的時(shí)空變化[J].草業(yè)科學(xué),2015,32(5):702-709.Xu Y F,Yang J,Chen Y N,et al.The spatial-temporal variations of vegetation cover in Xinjiang from 1982 to 2013[J].Pratacultural Science,2015,32(5):702-709.

    [24] 許玉鳳,楊井,李衛(wèi)紅,等.1982-2013年新疆不同植被生長(zhǎng)時(shí)空變化[J].草業(yè)學(xué)報(bào),2016,25(1):47-63.Xu Y F,Yang J,Li W H,et al.Spatial-temporal change in different vegetation growth of Xinjiang from 1982 to 2013[J].Acta Prataculturae Sinica,2016,25(1):47-63.

    猜你喜歡
    覆蓋面積覆蓋度平衡點(diǎn)
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測(cè)
    基于NDVI的晉州市植被覆蓋信息提取
    金星上全是金子嗎
    低覆蓋度CO分子在Ni(110)面的吸附研究
    覆蓋面積和網(wǎng)圍結(jié)構(gòu)對(duì)水面蒸發(fā)抑制率的影響
    探尋中國(guó)蘋(píng)果產(chǎn)業(yè)的產(chǎn)銷平衡點(diǎn)
    電視庭審報(bào)道,如何找到媒體監(jiān)督與司法公正的平衡點(diǎn)
    在給專車(chē)服務(wù)正名之前最好找到Uber和出租車(chē)的平衡點(diǎn)
    行走在預(yù)設(shè)與生成的平衡點(diǎn)上共同演繹精彩政治課堂
    散文百家(2014年11期)2014-08-21 07:16:58
    有機(jī)硅表面活性劑Tech-408和Fairland2408對(duì)農(nóng)藥?kù)F滴在煙草葉片上覆蓋面積的影響
    bbb黄色大片| 老司机亚洲免费影院| 一级a爱视频在线免费观看| 纯流量卡能插随身wifi吗| 午夜福利一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 国产精品成人在线| 亚洲精品国产色婷婷电影| 一区二区三区乱码不卡18| 国产精品久久久久久精品古装| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 在线观看三级黄色| 亚洲精品日本国产第一区| 自线自在国产av| 91精品三级在线观看| 这个男人来自地球电影免费观看 | 欧美精品高潮呻吟av久久| 久久女婷五月综合色啪小说| 精品一区二区免费观看| 色播在线永久视频| 人人妻,人人澡人人爽秒播 | 色精品久久人妻99蜜桃| xxx大片免费视频| 精品少妇内射三级| 亚洲图色成人| 国产一区二区激情短视频 | 国产福利在线免费观看视频| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 久久久精品区二区三区| 久久久久精品人妻al黑| 久久精品久久久久久噜噜老黄| 亚洲婷婷狠狠爱综合网| 18禁动态无遮挡网站| 欧美人与善性xxx| 丝袜人妻中文字幕| 精品一区二区三区av网在线观看 | 亚洲一码二码三码区别大吗| 欧美少妇被猛烈插入视频| www日本在线高清视频| 国产亚洲av片在线观看秒播厂| 国产高清不卡午夜福利| 亚洲精品乱久久久久久| 欧美国产精品va在线观看不卡| 精品亚洲成国产av| 欧美在线一区亚洲| 国产精品 欧美亚洲| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 欧美激情高清一区二区三区 | 午夜福利一区二区在线看| 久久久亚洲精品成人影院| 久久久国产一区二区| 免费在线观看视频国产中文字幕亚洲 | 国产成人av激情在线播放| 秋霞在线观看毛片| 国产探花极品一区二区| 成人黄色视频免费在线看| 天堂8中文在线网| 激情视频va一区二区三区| 亚洲,欧美精品.| 免费高清在线观看日韩| 男女之事视频高清在线观看 | 亚洲少妇的诱惑av| 美女中出高潮动态图| 亚洲精品美女久久av网站| 国产精品一国产av| 亚洲情色 制服丝袜| 激情视频va一区二区三区| 国产成人一区二区在线| 亚洲av中文av极速乱| 操出白浆在线播放| 日韩大片免费观看网站| 日韩制服骚丝袜av| av片东京热男人的天堂| 国产男女内射视频| 中文字幕最新亚洲高清| 最近中文字幕高清免费大全6| 99久久综合免费| 性色av一级| 男女床上黄色一级片免费看| 日韩伦理黄色片| 黄色毛片三级朝国网站| 十八禁高潮呻吟视频| 最近中文字幕2019免费版| 一级毛片 在线播放| 无限看片的www在线观看| 波多野结衣av一区二区av| av.在线天堂| 可以免费在线观看a视频的电影网站 | 性高湖久久久久久久久免费观看| 各种免费的搞黄视频| 亚洲自偷自拍图片 自拍| 欧美av亚洲av综合av国产av | 韩国高清视频一区二区三区| 伦理电影大哥的女人| 极品人妻少妇av视频| 精品国产国语对白av| 少妇精品久久久久久久| 免费人妻精品一区二区三区视频| 国产免费视频播放在线视频| 国产色婷婷99| 波多野结衣一区麻豆| 国产99久久九九免费精品| 国产欧美日韩一区二区三区在线| 免费女性裸体啪啪无遮挡网站| 2018国产大陆天天弄谢| 蜜桃国产av成人99| 一级毛片黄色毛片免费观看视频| 日韩av免费高清视频| 最近2019中文字幕mv第一页| 国产色婷婷99| 秋霞伦理黄片| 青春草亚洲视频在线观看| 99久久精品国产亚洲精品| 日韩制服骚丝袜av| 国产精品 欧美亚洲| 久久精品久久精品一区二区三区| 美女国产高潮福利片在线看| 免费黄色在线免费观看| 女的被弄到高潮叫床怎么办| 精品亚洲成a人片在线观看| videos熟女内射| 久久久欧美国产精品| 国产一区二区激情短视频 | 婷婷色av中文字幕| 精品人妻在线不人妻| 99精国产麻豆久久婷婷| 亚洲av男天堂| 久久人人爽人人片av| 在线精品无人区一区二区三| 国产精品秋霞免费鲁丝片| 欧美人与性动交α欧美软件| 黄片小视频在线播放| 欧美日韩视频精品一区| 日本欧美国产在线视频| 国产一卡二卡三卡精品 | 秋霞在线观看毛片| 国产精品av久久久久免费| av网站在线播放免费| 免费黄色在线免费观看| 99国产综合亚洲精品| 色综合欧美亚洲国产小说| 国产精品国产三级国产专区5o| 不卡视频在线观看欧美| 最新在线观看一区二区三区 | 在线亚洲精品国产二区图片欧美| 欧美日韩综合久久久久久| 色综合欧美亚洲国产小说| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 狂野欧美激情性bbbbbb| 亚洲av日韩精品久久久久久密 | 夫妻午夜视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲av国产av综合av卡| av国产久精品久网站免费入址| 欧美97在线视频| 又大又爽又粗| 好男人视频免费观看在线| 亚洲图色成人| 亚洲精品,欧美精品| 大片免费播放器 马上看| 一级爰片在线观看| 69精品国产乱码久久久| 国产精品一区二区精品视频观看| 久久久久久久久久久免费av| 伦理电影免费视频| av线在线观看网站| 亚洲国产欧美一区二区综合| 欧美乱码精品一区二区三区| 精品人妻熟女毛片av久久网站| e午夜精品久久久久久久| 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频| 最新在线观看一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 只有这里有精品99| 午夜日韩欧美国产| 久久精品熟女亚洲av麻豆精品| 久久性视频一级片| 老司机亚洲免费影院| 欧美黑人欧美精品刺激| 亚洲综合色网址| 亚洲av综合色区一区| 天天添夜夜摸| 黑人欧美特级aaaaaa片| avwww免费| 亚洲熟女精品中文字幕| 亚洲欧美精品综合一区二区三区| 久久国产精品大桥未久av| 欧美日韩亚洲国产一区二区在线观看 | netflix在线观看网站| 人妻 亚洲 视频| 99香蕉大伊视频| 老司机影院成人| 欧美黑人欧美精品刺激| 只有这里有精品99| 亚洲国产精品成人久久小说| 久久精品熟女亚洲av麻豆精品| 日本午夜av视频| 精品久久久精品久久久| 十八禁网站网址无遮挡| 丰满乱子伦码专区| www.自偷自拍.com| 午夜av观看不卡| 老司机亚洲免费影院| 国产探花极品一区二区| 男女国产视频网站| 丁香六月天网| 在线观看免费午夜福利视频| 一边摸一边做爽爽视频免费| 亚洲av日韩在线播放| 国产成人啪精品午夜网站| 最近中文字幕高清免费大全6| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品成人av观看孕妇| av线在线观看网站| 伦理电影大哥的女人| 哪个播放器可以免费观看大片| 日韩伦理黄色片| 激情五月婷婷亚洲| xxxhd国产人妻xxx| 丝袜脚勾引网站| 日韩视频在线欧美| 欧美最新免费一区二区三区| 制服诱惑二区| 我的亚洲天堂| 一边摸一边做爽爽视频免费| 国产成人一区二区在线| 69精品国产乱码久久久| 国产1区2区3区精品| 国产欧美日韩一区二区三区在线| 国产在线免费精品| 街头女战士在线观看网站| 日韩精品有码人妻一区| 国产日韩一区二区三区精品不卡| 91精品伊人久久大香线蕉| 日韩制服丝袜自拍偷拍| 性高湖久久久久久久久免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人一二三区av| 天堂中文最新版在线下载| 国产一区二区 视频在线| 18禁动态无遮挡网站| videosex国产| 免费高清在线观看视频在线观看| 日韩av在线免费看完整版不卡| 免费久久久久久久精品成人欧美视频| 日本av手机在线免费观看| 黄色 视频免费看| 国产免费又黄又爽又色| 欧美精品人与动牲交sv欧美| 卡戴珊不雅视频在线播放| 精品视频人人做人人爽| 成人国产av品久久久| 人人妻人人爽人人添夜夜欢视频| 丝袜人妻中文字幕| 欧美日韩视频高清一区二区三区二| 亚洲精品成人av观看孕妇| 欧美变态另类bdsm刘玥| av网站免费在线观看视频| 欧美中文综合在线视频| 97人妻天天添夜夜摸| 国产精品国产av在线观看| 香蕉国产在线看| av女优亚洲男人天堂| 免费高清在线观看日韩| 久久韩国三级中文字幕| 看免费成人av毛片| 9色porny在线观看| 两个人免费观看高清视频| 国产精品亚洲av一区麻豆 | 日韩大片免费观看网站| 老熟女久久久| 亚洲美女视频黄频| 老司机靠b影院| 天天躁夜夜躁狠狠久久av| 热re99久久国产66热| 高清视频免费观看一区二区| 欧美精品av麻豆av| 久久毛片免费看一区二区三区| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 这个男人来自地球电影免费观看 | 青草久久国产| 亚洲欧美激情在线| 丝袜人妻中文字幕| 欧美精品亚洲一区二区| 满18在线观看网站| 中文字幕精品免费在线观看视频| 女性生殖器流出的白浆| a级毛片黄视频| 免费看不卡的av| 高清av免费在线| 精品亚洲成国产av| 精品一区二区免费观看| av线在线观看网站| 精品一区二区免费观看| 伊人久久国产一区二区| 亚洲欧洲国产日韩| 久热爱精品视频在线9| av.在线天堂| 王馨瑶露胸无遮挡在线观看| 午夜福利免费观看在线| 亚洲精品av麻豆狂野| 宅男免费午夜| 久久ye,这里只有精品| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码 | 日韩 亚洲 欧美在线| 免费高清在线观看日韩| 99九九在线精品视频| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 亚洲精华国产精华液的使用体验| 黑人猛操日本美女一级片| 免费少妇av软件| 最近最新中文字幕免费大全7| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 免费黄网站久久成人精品| 欧美精品一区二区免费开放| 精品少妇一区二区三区视频日本电影 | 日本vs欧美在线观看视频| 亚洲精品国产av蜜桃| 街头女战士在线观看网站| 男人爽女人下面视频在线观看| 国产伦人伦偷精品视频| 在现免费观看毛片| 精品免费久久久久久久清纯 | 少妇精品久久久久久久| 少妇人妻 视频| 国产一级毛片在线| 国产精品无大码| 青春草国产在线视频| 欧美精品人与动牲交sv欧美| 9热在线视频观看99| 久久人妻熟女aⅴ| 精品久久久精品久久久| 黄频高清免费视频| 老司机影院成人| 51午夜福利影视在线观看| 咕卡用的链子| 精品国产一区二区三区久久久樱花| 亚洲成人手机| 人体艺术视频欧美日本| 欧美人与性动交α欧美软件| 久久99一区二区三区| 最近手机中文字幕大全| 免费观看av网站的网址| 又黄又粗又硬又大视频| 黄频高清免费视频| 在现免费观看毛片| 美女大奶头黄色视频| 高清欧美精品videossex| 久久久久久久久久久免费av| 国产免费视频播放在线视频| 久久国产精品大桥未久av| 热re99久久国产66热| 90打野战视频偷拍视频| 亚洲欧美清纯卡通| 精品免费久久久久久久清纯 | 国产免费又黄又爽又色| 18禁裸乳无遮挡动漫免费视频| 丁香六月天网| 一本色道久久久久久精品综合| 极品少妇高潮喷水抽搐| 亚洲免费av在线视频| 日本91视频免费播放| 亚洲美女视频黄频| 十八禁高潮呻吟视频| e午夜精品久久久久久久| 中文字幕亚洲精品专区| 色精品久久人妻99蜜桃| netflix在线观看网站| 啦啦啦在线免费观看视频4| 久久久久久人妻| 日本av免费视频播放| 午夜激情av网站| 男女之事视频高清在线观看 | 亚洲av成人精品一二三区| 老司机亚洲免费影院| 亚洲成色77777| 中文字幕色久视频| 亚洲精品中文字幕在线视频| 你懂的网址亚洲精品在线观看| 亚洲,一卡二卡三卡| 精品一品国产午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 69精品国产乱码久久久| 久久久久久免费高清国产稀缺| 国产精品 国内视频| 老司机影院成人| 狠狠婷婷综合久久久久久88av| 午夜福利视频精品| 老汉色av国产亚洲站长工具| 久久午夜综合久久蜜桃| 新久久久久国产一级毛片| 香蕉国产在线看| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 超碰成人久久| 免费高清在线观看日韩| 久久精品久久精品一区二区三区| 日韩一区二区视频免费看| 中国三级夫妇交换| 精品亚洲乱码少妇综合久久| 日韩不卡一区二区三区视频在线| 麻豆精品久久久久久蜜桃| av网站免费在线观看视频| 老司机深夜福利视频在线观看 | 国产精品久久久久久精品古装| 中文天堂在线官网| 亚洲综合精品二区| 波野结衣二区三区在线| 免费观看性生交大片5| 亚洲欧美成人精品一区二区| 国产激情久久老熟女| 一级爰片在线观看| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 亚洲第一区二区三区不卡| 久久青草综合色| 亚洲,欧美精品.| 丰满乱子伦码专区| 午夜福利乱码中文字幕| 亚洲国产欧美在线一区| 熟妇人妻不卡中文字幕| 制服诱惑二区| 建设人人有责人人尽责人人享有的| 丝袜人妻中文字幕| 一边摸一边抽搐一进一出视频| 亚洲成人av在线免费| xxxhd国产人妻xxx| 黄片小视频在线播放| 欧美另类一区| 精品国产一区二区三区四区第35| 美女主播在线视频| 欧美成人午夜精品| 亚洲一区二区三区欧美精品| 老司机在亚洲福利影院| 久久精品人人爽人人爽视色| 精品第一国产精品| 成人午夜精彩视频在线观看| 国产成人免费观看mmmm| 自线自在国产av| 中文字幕人妻熟女乱码| 超碰成人久久| 午夜福利在线免费观看网站| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 欧美国产精品va在线观看不卡| 久久鲁丝午夜福利片| 黄片播放在线免费| 欧美黑人精品巨大| 中文精品一卡2卡3卡4更新| 国产在线一区二区三区精| 欧美精品一区二区免费开放| 在线免费观看不下载黄p国产| 十八禁网站网址无遮挡| 夫妻午夜视频| 国产免费又黄又爽又色| xxx大片免费视频| 精品第一国产精品| 精品少妇一区二区三区视频日本电影 | 久久精品aⅴ一区二区三区四区| 一本—道久久a久久精品蜜桃钙片| 精品人妻一区二区三区麻豆| 国产精品秋霞免费鲁丝片| 十八禁高潮呻吟视频| 18禁裸乳无遮挡动漫免费视频| 欧美在线黄色| 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| 久久久精品94久久精品| 另类亚洲欧美激情| 多毛熟女@视频| 熟女av电影| 日本爱情动作片www.在线观看| 亚洲视频免费观看视频| 黑人猛操日本美女一级片| 两个人免费观看高清视频| 一级毛片电影观看| 18禁国产床啪视频网站| 亚洲av电影在线观看一区二区三区| 国产精品国产av在线观看| 青春草亚洲视频在线观看| 亚洲在久久综合| 国产1区2区3区精品| av有码第一页| netflix在线观看网站| 精品国产乱码久久久久久男人| 一级毛片我不卡| 51午夜福利影视在线观看| 美女国产高潮福利片在线看| 国产精品久久久久久人妻精品电影 | 青草久久国产| a级毛片黄视频| 永久免费av网站大全| 午夜福利视频精品| 蜜桃在线观看..| 国产亚洲最大av| 成人午夜精彩视频在线观看| 免费久久久久久久精品成人欧美视频| 国产精品香港三级国产av潘金莲 | 校园人妻丝袜中文字幕| 国产成人精品福利久久| 一区二区三区乱码不卡18| 99国产综合亚洲精品| 在线 av 中文字幕| 满18在线观看网站| 国产精品嫩草影院av在线观看| 成人亚洲欧美一区二区av| 精品久久久精品久久久| 90打野战视频偷拍视频| 精品国产国语对白av| 午夜福利免费观看在线| 美国免费a级毛片| 国产福利在线免费观看视频| 久久精品亚洲熟妇少妇任你| 欧美在线一区亚洲| 国产 精品1| 十八禁人妻一区二区| 黄色视频在线播放观看不卡| 亚洲欧美一区二区三区黑人| 免费日韩欧美在线观看| 母亲3免费完整高清在线观看| 国产日韩欧美视频二区| 国产精品一国产av| av线在线观看网站| 色吧在线观看| 精品一品国产午夜福利视频| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠久久av| 亚洲精品国产区一区二| 国产极品粉嫩免费观看在线| 久久久国产欧美日韩av| 黑人猛操日本美女一级片| 欧美黑人欧美精品刺激| 中文字幕高清在线视频| 国产av码专区亚洲av| a 毛片基地| 亚洲色图 男人天堂 中文字幕| 亚洲精品一区蜜桃| 久久久久国产一级毛片高清牌| 久久婷婷青草| 黄色怎么调成土黄色| 天天添夜夜摸| 午夜激情av网站| 制服诱惑二区| 伊人久久国产一区二区| 日韩av免费高清视频| 亚洲男人天堂网一区| 一区二区三区四区激情视频| 老司机影院毛片| 亚洲精品日本国产第一区| 18在线观看网站| 女的被弄到高潮叫床怎么办| 国产麻豆69| 国产淫语在线视频| 久久人人爽av亚洲精品天堂| 久久久久久久久久久免费av| 免费在线观看视频国产中文字幕亚洲 | 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人| 99热全是精品| 久久精品久久久久久噜噜老黄| 丝袜喷水一区| 欧美国产精品一级二级三级| 国产探花极品一区二区| 午夜日韩欧美国产| 国产福利在线免费观看视频| 国产亚洲av高清不卡| 婷婷色综合大香蕉| 精品亚洲成a人片在线观看| 国产精品三级大全| 老汉色av国产亚洲站长工具| 人妻人人澡人人爽人人| 午夜福利网站1000一区二区三区| 国产 精品1| 欧美日韩成人在线一区二区| 亚洲国产成人一精品久久久| 操出白浆在线播放| 久久精品国产a三级三级三级| 国产日韩一区二区三区精品不卡| 天堂中文最新版在线下载| 久久久久久久国产电影| 激情五月婷婷亚洲| 久久人妻熟女aⅴ| 999久久久国产精品视频| 搡老乐熟女国产| 大话2 男鬼变身卡| 亚洲图色成人| 观看美女的网站| 欧美在线一区亚洲| av在线app专区| 久久精品久久精品一区二区三区| 日韩av在线免费看完整版不卡| 少妇精品久久久久久久| 久久 成人 亚洲| 丝袜人妻中文字幕| 国产xxxxx性猛交| 嫩草影院入口| 欧美日韩视频精品一区| 午夜福利一区二区在线看| 亚洲欧美清纯卡通| 日本一区二区免费在线视频| 国产精品久久久人人做人人爽| 老司机深夜福利视频在线观看 | 国产黄频视频在线观看|