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

    西北干旱地區(qū)氣候變化及其對(duì)草地生產(chǎn)潛力的影響

    2020-06-08 11:13:02張娟李劍萍王譽(yù)陶張翼井樂李建平
    生態(tài)科學(xué) 2020年3期
    關(guān)鍵詞:生產(chǎn)潛力降水量寧夏

    張娟, 李劍萍, 王譽(yù)陶, 張翼, 井樂, 李建平, 3, *

    西北干旱地區(qū)氣候變化及其對(duì)草地生產(chǎn)潛力的影響

    張娟1, 李劍萍2, 王譽(yù)陶1, 張翼1, 井樂1, 李建平1, 3, *

    1. 寧夏大學(xué)農(nóng)學(xué)院, 寧夏 銀川 750021 2. 寧夏防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室, 寧夏 銀川 750021 3. 西北退化生態(tài)系統(tǒng)恢復(fù)與重建國家重點(diǎn)實(shí)驗(yàn)室培育基地, 寧夏 銀川 750021

    依據(jù)寧夏各縣區(qū)24個(gè)氣象站37年(1981—2017年)間的逐月氣溫、降水資料, 采用Miami模型和Thornthwaite Memorial模型計(jì)算寧夏草地生產(chǎn)潛力, 分析其時(shí)空變化。結(jié)果表明: 從時(shí)間變化來看, 寧夏平均氣溫以0.47 ℃·(10a)-1的速度持續(xù)上升, 降水量呈先下降后上升的趨勢(shì), 從2010—2017年年降雨量增加了40 mm; 從空間變化來看, 寧夏氣候從西北向東南由“暖干”向“冷濕”逐漸變化; 寧夏草地氣候生產(chǎn)潛力從西北向東南逐漸遞增; 草地氣候生產(chǎn)潛力變化趨勢(shì)與降水變化趨勢(shì)一致, 各時(shí)間段氣候生產(chǎn)潛力與降水呈顯著正相關(guān)(<0.01), 與溫度呈負(fù)相關(guān)關(guān)系。因此, 降水是寧夏草地氣候生產(chǎn)潛力的限制因素, 草地生產(chǎn)潛力較大區(qū)域在寧夏中南部, 而以溫度升高為特征的氣候變化, 對(duì)寧夏草地氣候生產(chǎn)潛力影響較小。

    草地生產(chǎn)潛力; 氣候變化; 寧夏; Miami模型

    0 前言

    草地是全球陸地生態(tài)系統(tǒng)的重要組成部分, 也是最常見的植被類型之一, 在食物安全、生態(tài)、碳平衡和全球氣候變化中發(fā)揮著關(guān)鍵作用。我國草地面積為4億公頃, 占全國土地面積的40%[1-2]。草地對(duì)全球變化較為敏感[3], 近年來, 氣候變暖、干旱和極端氣候事件使草原植物物候發(fā)生改變, 初級(jí)生產(chǎn)力和草原載畜能力降低[4]。

    凈初級(jí)生產(chǎn)力是指綠色植物在單位面積、單位時(shí)間內(nèi)光合作用所積累的有機(jī)物總量[5], 反映植物群落的固碳能力, 是植被生長狀況和生態(tài)系統(tǒng)健康的重要指標(biāo)[6], 也是研究氣候變化對(duì)陸地生態(tài)系統(tǒng)影響的重要參數(shù)[7]。草地地上凈初級(jí)生產(chǎn)力與年平均溫度和降水密切相關(guān)[8]。IPCC第五次評(píng)估報(bào)告指出未來幾十年全球變暖將會(huì)持續(xù)下去, 對(duì)草地生產(chǎn)力產(chǎn)生顯著影響[9]。草地氣候生產(chǎn)潛力是草地凈初級(jí)生產(chǎn)力的一個(gè)方面, 研究草地生產(chǎn)潛力的時(shí)空分布及其對(duì)氣候變化的響應(yīng)對(duì)于合理利用草地資源、提高草地生產(chǎn)力、保障草地生態(tài)安全和可持續(xù)發(fā)展具有重要的指導(dǎo)意義[10-11]。

    草地是中國干旱半干旱區(qū)主要的植被類型, 也是我國北方地區(qū)最重要的生態(tài)屏障[12]。寧夏回族自治區(qū)(以下簡稱“寧夏”)位于我國北方農(nóng)牧交錯(cuò)區(qū)和生態(tài)脆弱帶, 屬于中溫帶干旱半干旱氣候區(qū), 天然草地面積為301.4萬公頃, 草原類型以北部荒漠草原類和南部溫性草原類為主, 全區(qū)降水稀少, 天然草地生產(chǎn)力的高低主要取決于當(dāng)?shù)貧夂颦h(huán)境及其變化[13]。本研究以西北干旱地區(qū)草地為研究區(qū), 以寧夏24個(gè)氣象站點(diǎn)37年間(1981—2017年)的逐月氣溫、降水資料為基礎(chǔ), 運(yùn)用Miami模型、Thornthwaite Memorial模型和克里格空間插值(Kriging) 法等分析草地生產(chǎn)潛力的時(shí)空變化特征, 以期揭示氣候變化對(duì)草地生產(chǎn)力影響的時(shí)空規(guī)律及其對(duì)氣候變化的響應(yīng), 為草地可持續(xù)發(fā)展、生產(chǎn)潛力動(dòng)態(tài)檢測(cè)以及草—畜平衡布局提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    寧夏位于中國西北部(35°14'—39°23'N, 104°17'—107°39'E), 總面積為6.64萬平方千米, 是我國北方典型的農(nóng)牧交錯(cuò)區(qū), 生態(tài)條件脆弱, 水土流失、草地退化、沙漠化等問題嚴(yán)重。地貌類型多樣, 北部為干旱剝蝕、風(fēng)蝕地貌, 南部為流水侵蝕地貌。研究區(qū)屬溫帶大陸性干旱半干旱氣候, 年平均降水量270—320 mm, 多集中于7—9月, 年平均氣溫為7—9.2 ℃, 無霜期120—218 d, 年日照時(shí)數(shù)2710—3124 h。植被以草原為主, 主要植被類型為溫性荒漠草原、溫性草原、溫性草原化荒漠、溫性山地草甸等(圖1)。

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

    氣象資料來源于寧夏各縣區(qū)24個(gè)氣象站點(diǎn)的37年間(1981—2017年)的逐月氣溫、降水資料。

    注: 來源于《寧夏草地資源》附圖。

    Figure 1 Distribution pattern of grassland types in Ningxia

    1.3 計(jì)算方法

    1.3.1 Miami模型

    該模型根據(jù)氣溫、降水量兩個(gè)氣候因子分別計(jì)算溫度生產(chǎn)潛力和降水生產(chǎn)潛力。

    Y=30000/(1+1.315–0.119T) (1)

    Y=30000/(1––0.000664R) (2)

    式中:為年平均氣溫(℃);為年降水量(mm);YY為根據(jù)年平均氣溫和年降水量計(jì)算的草地生產(chǎn)潛力, 稱為溫度生產(chǎn)潛力和降水生產(chǎn)潛力(kg·hm-2·a-1)。

    1.3.2 Thornthwaite Memorial模型

    該模型依據(jù)蒸散量計(jì)算草地生產(chǎn)潛力, 其中蒸散量是氣溫和降水的函數(shù)。

    Y=30000×[1––0.0009695(E–20)] (3)

    式中:Y為根據(jù)蒸散量計(jì)算的草地生產(chǎn)潛力(kg·hm–2·a–1);為實(shí)際年平均蒸散量(mm), 計(jì)算公式如下:

    =1.05×[1+(1.05/)2]–0.5

    式中:為年降水量(mm);為年平均最大蒸散量(mm), 計(jì)算公式如下:

    =300+25+0.053

    式中:為年平均氣溫(℃)。

    1.3.3 草地氣候生產(chǎn)潛力

    根據(jù)Liebig定律取(1)—(3)中的最小值作為寧夏草地的氣候生產(chǎn)潛力(kg·hm-2·a-1)。

    =min(Y,Y,Y)

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

    利用Origin 8.0對(duì)氣象資料進(jìn)行處理, 采用ArcGIS 10.2中的克里格法空間插值模塊對(duì)24個(gè)氣象站點(diǎn)的點(diǎn)數(shù)據(jù)進(jìn)行分時(shí)段(1981—1990年、1991—2000年、2001—2010年、2011—2017年)處理, 獲得寧夏年平均氣溫、年降水量和草地氣候生產(chǎn)潛力的空間分布特征圖, 用R 3.5.1對(duì)氣溫、降水、蒸散量與草地生產(chǎn)潛力進(jìn)行Pearson相關(guān)性分析。

    2 結(jié)果與分析

    2.1 寧夏氣候變化

    2.1.1 年平均氣溫、年降水量的年際變化

    寧夏37年間的年平均氣溫呈明顯的上升趨勢(shì)(圖2), 每10年增加0.47 ℃(<0.01), 氣溫平均值為8.40 ℃。1996年之前, 氣溫普遍位于平均值以下, 1996年之后, 氣溫普遍位于平均值以上, 這與全球氣候變暖的趨勢(shì)一致。最低值出現(xiàn)在1984年, 為6.76 ℃, 最高值出現(xiàn)在2013年, 為9.41 ℃。

    圖2 寧夏年平均氣溫變化曲線

    Figure 2 Ningxia annual average temperature curve

    寧夏37年間的降水量變化趨勢(shì)不明顯(圖3), 圍繞平均值上下波動(dòng), 波動(dòng)幅度較大, 降水的年際變化大。降水量平均值為281.69 mm, 最高值出現(xiàn)在1990年, 為371.91 mm, 最低值出現(xiàn)在1982年, 為172.30 mm。

    2.1.2 年平均氣溫的空間分布

    寧夏年均氣溫從西北向東南逐漸降低(圖4), 在1981—2017年的四個(gè)時(shí)間段內(nèi), 各縣區(qū)的氣溫均有所升高, 且升高幅度逐漸減小。1981—1990年(圖4a)、1991—2000年(圖4b)、2001—2010年(圖4c)、2011—2017年(圖4d)四個(gè)時(shí)段內(nèi)寧夏各縣區(qū)平均氣溫分別為7.71 ℃、8.33 ℃、8.79 ℃、8.92 ℃。

    1981—2017年間, 寧夏南北溫差增大, 溫差由8.84 ℃上升到9.25 ℃。最高氣溫出現(xiàn)在中部干旱帶的西北部以及北部引黃灌區(qū)的西部, 最低氣溫出現(xiàn)在南部山地丘陵區(qū), 六盤山海拔較其他地區(qū)高, 是研究區(qū)氣溫最低點(diǎn)。

    圖3 寧夏年降水量變化曲線

    Figure 3 Ningxia annual precipitation curve

    圖4 寧夏年平均氣溫空間分布

    Figure 4 Spatial distribution of annual average temperature in Ningxia

    1981—1990年, 寧夏最高年均氣溫較低, 不超過10.50 ℃。氣溫在9.01 ℃—10.50 ℃之間的地區(qū)包括青銅峽市、中寧縣北部、紅寺堡區(qū)中西部和西夏區(qū)、賀蘭縣、平羅縣、大武口區(qū)等地的大部分范圍, 1991—2000年, 范圍擴(kuò)大至引黃灌區(qū)和中部干旱帶大部分地區(qū), 2001—2010年, 范圍繼續(xù)擴(kuò)大, 到2011—2017年, 范圍基本覆蓋了引黃灌區(qū)和中部干旱帶。2001—2010年, 年平均氣溫首先在青銅峽市中南部、中寧縣北部和利通區(qū)西北部突破10.50 ℃, 此后范圍逐步擴(kuò)大至紅寺堡區(qū)中西部、中寧縣中北部、利通區(qū)大部, 此外, 在永寧縣、西夏區(qū)和金鳳區(qū)的部分地區(qū)的年平均氣溫也突破10.50 ℃。南部山區(qū)六盤山附近氣溫最低, 1981—1990年, 該地區(qū)的最高氣溫不超過1.50 ℃, 2011—2017年, 該地最低氣溫超過1.50 ℃。

    2.1.3 年降水量的空間分布

    與氣溫變化梯度方向相反, 降水量的空間分布呈現(xiàn)由西北向東南遞增的趨勢(shì)(圖5)。各縣區(qū)的降水量總體穩(wěn)定, 年際變化不明顯。1981—1990年(圖5a)、1991—2000年(圖5b)、2001—2010年(圖5c)、2011—2017年(圖5d)四個(gè)時(shí)間段內(nèi)各縣區(qū)的平均降水量分別為285.38 mm、270.31 mm、270.00 mm、309.38 mm。整體上, 從1981—2010年, 各縣區(qū)降水量沿降水梯度稍有減少, 2010年后降水量增加。自治區(qū)降水量空間變化范圍大致為160—700 mm, 北部引黃灌區(qū)和中部干旱帶的西北部地區(qū)降水量不超過220 mm。1991—2000年(圖5b)平羅縣西部和大武口東南部地區(qū)降水量低于160 mm。涇源縣降雨量最多, 年降水量總體上大于640 mm。

    圖5 寧夏年降水量空間分布

    Figure 5 Spatial distribution of annual precipitation in Ningxia

    2.2 草地氣候生產(chǎn)潛力計(jì)算結(jié)果

    根據(jù)Liebig定律, 將氣溫、降水、蒸散量所決定的草地生產(chǎn)潛力的最小值作為當(dāng)?shù)氐臍夂蛏a(chǎn)潛力值。從表1中可以看出, 由氣溫決定的生產(chǎn)潛力值高于由降水和蒸散量決定的生產(chǎn)潛力值(六盤山除外), 由蒸散量決定的生產(chǎn)潛力值略高于由降水決定的生產(chǎn)潛力值, 因此降水是寧夏各地區(qū)草地生產(chǎn)潛力的限制因素。六盤山的氣溫明顯低于其他地區(qū), 降水量遠(yuǎn)遠(yuǎn)高于其他地區(qū), 因此氣溫是六盤山草地生產(chǎn)潛力的主要限制因素。各地區(qū)草地生產(chǎn)潛力有明顯差異, 從引黃灌區(qū)到南部山區(qū)的草地生產(chǎn)潛力有逐漸升高的趨勢(shì)。

    從氣溫、降水量、蒸散量與草地生產(chǎn)潛力相關(guān)分析可以看出, 氣溫()和溫度生產(chǎn)潛力()、降水量()和降水生產(chǎn)潛力()、蒸散量()和由蒸散量決定的生產(chǎn)潛力()均呈正相關(guān)關(guān)系(<0.01)(圖6), 由、、三者中的最小值構(gòu)成的氣候生產(chǎn)潛力()與降水量呈顯著正相關(guān)(<0.01), 充分說明降水量是寧夏草地氣候生產(chǎn)潛力的決定因素。蒸散量由氣溫和降水量共同影響, 由蒸散量決定的草地生產(chǎn)潛力實(shí)際上是由氣溫和降水量共同決定的, 氣溫升高, 降水量減少, 草地生產(chǎn)潛力降低。

    表1 寧夏草地生產(chǎn)潛力

    注:為氣溫,為降水量,為蒸散量,Y、Y、Y為根據(jù)氣溫、降水量、蒸散量計(jì)算的草地生產(chǎn)潛力,為草地氣候生產(chǎn)潛力。

    2.3 草地氣候生產(chǎn)潛力分析

    2.3.1 草地氣候生產(chǎn)潛力的年際變化

    草地氣候生產(chǎn)潛力自1980年以來呈緩慢上升趨勢(shì)(圖7), 圍繞平均值上下波動(dòng), 波動(dòng)幅度大, 草地氣候生產(chǎn)潛力年際變化大。最高值出現(xiàn)在1985年, 為6223.99 kg·hm-2·a-1, 最低值出現(xiàn)在1982年, 為3094.87 kg·hm-2·a-1。草地氣候生產(chǎn)潛力年際變化特征與年降水量年際變化特征一致, 充分說明降水是寧夏草地氣候生產(chǎn)潛力的限制因素。

    注: T為氣溫, R為降水量, E為蒸散量, YT、YR、YE為根據(jù)氣溫、降水量、蒸散量計(jì)算的草地生產(chǎn)潛力, Y為草地氣候生產(chǎn)潛力。

    Figure 6 Correlation analysis of temperature, precipitation, evaporation and potential grassland productivity

    2.3.2 草地氣候生產(chǎn)潛力的空間分布

    寧夏草地氣候生產(chǎn)潛力的空間分布特征與年降水量的空間分布特征一致。1981—2017年各區(qū)域草地氣候生產(chǎn)潛力由西北向東南逐級(jí)遞增(圖8)。西北地區(qū)氣候生產(chǎn)潛力范圍在3000 kg·hm-2·a-1—4000 kg·hm-2·a-1之間, 1991—2000年西北地區(qū)大武口區(qū)出現(xiàn)最低值, 為2800.96 kg·hm-2·a-1(圖8b)。涇源縣南部的草地氣候生產(chǎn)潛力最高, 最低值在9000 kg·hm-2·a-1以上, 最高可達(dá)到10515.70 kg·hm-2·a-1(圖8d)。

    圖7 寧夏草地氣候生產(chǎn)潛力變化曲線

    Figure 7 Curve of climate production potential change in Ningxia grasslan

    1981—2017年, 各縣區(qū)的草地氣候生產(chǎn)潛力先減少后增加。1981—1990年(圖8a), 各縣區(qū)的氣候生產(chǎn)潛力平均值為4819.01 kg·hm-2·a-1, 1991—2000年, 平均值下降至4714.78 kg·hm-2·a-1(圖8b), 2001—2010年(圖8c)達(dá)到最低值, 為4710.94 kg·hm-2·a-1, 2011—2017年(圖8d)增加至5000 kg·hm-2·a-1以上。

    1981—1990年(圖9a)、1991—2000年(圖9b)、2001—2010年(圖9c)、2011—2017年(圖9d)四個(gè)時(shí)間段內(nèi), 降水量、蒸散量和草地氣候生產(chǎn)潛力均呈顯著正相關(guān)(<0.01), 從空間分布來看, 南部山區(qū)氣溫低, 降水量高, 蒸散量大, 草地氣候生產(chǎn)潛力高。

    3 結(jié)論和討論

    從時(shí)間變化來看, 研究區(qū)年平均氣溫呈顯著上升趨勢(shì)(<0.01), 增幅為0.47 ℃·(10a)-1, 與全球氣候變暖現(xiàn)象一致; 降水量呈緩慢上升的趨勢(shì), 未達(dá)到顯著水平。其中, 1981—2010年, 降水量呈下降趨勢(shì),2010年以后, 降水量有所增加, 這與內(nèi)蒙古中部草原的氣候變化特征相同[14]。從空間變化來看, 研究區(qū)氣候條件從西北向東南由“暖干”逐漸向“冷濕”轉(zhuǎn)變。馬甜[15]等人認(rèn)為, 寧夏中部干旱帶氣候在1961—2008年間呈現(xiàn)“暖干化”發(fā)展趨勢(shì), 羅瑞敏[16]等人認(rèn)為云霧山草原1992—2011年間降水量呈下降趨勢(shì), 這與本研究中1981—2010年間降水量變化趨勢(shì)吻合。寧夏自治區(qū)位于農(nóng)牧交錯(cuò)帶, 生態(tài)條件脆弱, 長期不合理的放牧方式導(dǎo)致草地被破壞, 生物量減少, 地皮裸露沙化, 干旱加劇, 這可能是2010年前寧夏各區(qū)域年降水量減少的原因, 從2000年以來, 自治區(qū)實(shí)施封山禁牧政策, 10年來, 森林覆蓋率大幅增加, 沙化面積減少, 水土流失得到遏制, 下墊面的改變促進(jìn)降水量增加。氣溫與降水量、蒸散量、草地氣候生產(chǎn)潛力呈顯著負(fù)相關(guān)(<0.01)。氣溫上升, 溫度生產(chǎn)潛力增加, 但氣溫上升引起蒸散量增加, 從而導(dǎo)致氣候生產(chǎn)潛力下降[16]。

    圖8 寧夏草地氣候生產(chǎn)潛力空間分布

    Figure 8 Spatial distribution of potential climatic grassland productivity in Ningxia

    注: T1、T2、T3、T4為氣溫, R1、R2、R3、R4為降水量, E1、E2、E3、E4為蒸散量, YT1、YT2、YT3、YT4、YR1、YR2、YR3、YR4、YE1、YE2、YE3、YE4為根據(jù)氣溫、降水量、蒸散量計(jì)算的草地生產(chǎn)潛力, Y1、Y2、Y3、Y4為草地氣候生產(chǎn)潛力。

    Figure 9 Correlation analysis of temperature, precipitation, evaporation and potential grassland productivity (four time periods)

    在干旱半干旱區(qū), 草地生態(tài)系統(tǒng)對(duì)氣候變化極為敏感, 溫度和降水決定植被生長的水熱條件, 從而影響草地氣候生產(chǎn)潛力[17]。溫度升高促進(jìn)植物生長, 但高溫又會(huì)增加干燥化[18], 降低草地氣候生產(chǎn)潛力, 而增加降雨量可緩解干燥化, 使草地氣候生產(chǎn)潛力提高。因此, 不同研究區(qū)域特有的氣候特征決定著該區(qū)域的草地生產(chǎn)潛力, 如我國東部地區(qū)降水豐沛, 草地生產(chǎn)潛力主要受制于溫度, 而西北干旱地區(qū)受制于降水[19]。在本研究中, 寧夏草地氣候生產(chǎn)潛力與降水量呈正相關(guān)(<0.01), 二者變化特征一致, 降水是當(dāng)?shù)貧夂蛏a(chǎn)潛力的決定因子。李廣[20]、孫慧蘭[21]等的研究也證實(shí)了降水是制約半干旱區(qū)草地生產(chǎn)潛力的主控因素。龔相文等人研究發(fā)現(xiàn), 屬干旱半干旱氣候的科爾沁地區(qū)草地氣候生產(chǎn)潛力與降水量呈極顯著正相關(guān), 水分是當(dāng)?shù)刂脖簧L的主要限制因子[22], 這和本研究結(jié)果一致, 即降水是寧夏草地氣候生產(chǎn)潛力的制約因素, 溫度變化對(duì)該區(qū)草地生產(chǎn)潛力影響較小。

    研究區(qū)草地生產(chǎn)潛力空間分布上從東南向西北遞減, 與我國生物生產(chǎn)量的空間分布規(guī)律一致[23]。寧夏草地氣候生產(chǎn)潛力范圍為3167.25 kg·hm-2·a-1—10212.89 kg·hm-2·a-1, 寧夏各縣區(qū)平均值為4869.46 kg·hm-2·a-1, 北部引黃灌區(qū)氣候生產(chǎn)潛力平均值為3332.65 kg·hm-2·a-1, 中部干旱帶氣候生產(chǎn)潛力平均值為4321.04 kg·hm-2·a-1, 南部山區(qū)氣候生產(chǎn)潛力平均值為7832.58 kg·hm-2·a-1。生產(chǎn)潛力的限制因素隨空間分布而不同[24], 引黃灌區(qū)和中部干旱帶地勢(shì)平坦, 海拔低, 降水少, 氣候生產(chǎn)潛力小, 南部山區(qū)海拔高, 地形復(fù)雜, 多地形雨, 氣候生產(chǎn)潛力大。寧夏各縣區(qū)草地氣候生產(chǎn)潛力總體的限制因素是降水, 六盤山因海拔較高, 氣溫較低, 降水豐富, 限制因素為氣溫。

    采用不同的模型進(jìn)行估算獲得的草地生產(chǎn)潛力值不同。本研究采用Thornthwaite Memorial 模型和Miami模型估算草地生產(chǎn)潛力, 兩者相比, 前者考慮的氣候因素更全面, 計(jì)算結(jié)果偏差較小[25]。但兩種模型都沒有考慮到土壤特性、草地利用狀況等因素對(duì)草地生產(chǎn)潛力的影響, 計(jì)算結(jié)果與實(shí)際情況的擬合程度不夠好。朱玉果[26]等人采用CASA模型對(duì)寧夏草地凈初級(jí)生產(chǎn)力進(jìn)行估算, 結(jié)果與當(dāng)?shù)夭莸厣a(chǎn)力實(shí)測(cè)值擬合較好。周平[27]等人對(duì)8種模型進(jìn)行分析后發(fā)現(xiàn), 周廣勝模型和LPA模型對(duì)草地生產(chǎn)力的模擬結(jié)果較好。下一步研究將對(duì)估算模型進(jìn)行優(yōu)化選擇, 使研究結(jié)果更加精準(zhǔn)。

    [1] LIU Yangyang, YANG Yue, WANG Qian, et al. Evaluating the responses of net primary productivity and carbon use efficiency of global grassland to climate variability along an aridity gradient[J]. Science of the Total Environment, 2019, 652: 671–682.

    [2] WANG Zhaoqi, ZHANG Yanzhen, YANG Yue, et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai-Tibet Plateau, in China[J]. Ecological Informatics, 2016, 33: 32–44.

    [3] ZHANG Ying, ZHANG Chaobin, WANG Zhaoqi, et al. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012[J]. Science of The Total Environment, 2016, 563/564: 210–220.

    [4] ZHANG Haiyan, FAN Jiangwen, WANG Junbang, et al. Spatial and temporal variability of grassland yield and its response to climate change and anthropogenic activities on the Tibetan Plateau from 1988 to 2013[J]. Ecological Indicators, 2018, 95: 141–151.

    [5] 任璇, 鄭江華, 穆晨, 等. 新疆15年草地NPP動(dòng)態(tài)變化與氣象因子的相關(guān)性研究[J]. 生態(tài)科學(xué), 2017, 36(3): 43– 51.

    [6] GANG Chengcheng, ZHAO Wei, ZHAO Ting, et al. The impacts of land conversion and management measures on the grassland net primary productivity over the Loess Plateau, Northern China[J]. Science of The Total Environ-ment, 2018, 645: 827–836.

    [7] 劉鳳, 曾永年. 近16年青海高原植被NPP時(shí)空格局變化及氣候與人為因素的影響[J]. 生態(tài)學(xué)報(bào), 2019, 39(5): 1528–1540.

    [8] GUO Liang, CHENG Jimin, LUEDELING E, et al. Critical climate periods for grassland productivity on China’s Loess Plateau[J]. Agricultural and Forest Meteorology, 2017, 233: 101–109.

    [9] 耿元波, 王松, 胡雪荻. 高寒草甸草原凈初級(jí)生產(chǎn)力對(duì)氣候變化響應(yīng)的模擬[J]. 草業(yè)學(xué)報(bào), 2018, 27(1): 1–13.

    [10] 羅瑞敏, 郭梁, 程積民. 近20 a云霧山草地生產(chǎn)力對(duì)氣候變化的響應(yīng)[J]. 干旱區(qū)研究, 2018, 35(1): 77–84.

    [11] 王松, 耿元波, 母悅. 典型草原凈初級(jí)生產(chǎn)力對(duì)氣候變化響應(yīng)的模擬[J]. 草業(yè)學(xué)報(bào), 2016, 25(12): 4–13.

    [12] HAN Dongmei, WANG Guoqiang, LIU Tiangxi, et al. Hydroclimatic response of evapotranspiration partitioning to prolonged droughts in semiarid grassland[J]. Journal of Hydrology, 2018, 563: 766–777.

    [13] 李艷春. 寧夏干旱區(qū)氣候承載力分布特征分析[J]. 干旱區(qū)資源與環(huán)境, 2010, 24(8): 96–99.

    [14] 高浩, 潘學(xué)標(biāo), 符瑜. 氣候變化對(duì)內(nèi)蒙古中部草原氣候生產(chǎn)潛力的影響[J]. 中國農(nóng)業(yè)氣象, 2009, 30(3): 277–282, 288.

    [15] 馬甜, 王俊波, 張治華, 等. 寧夏中部干旱帶天然草地氣候生產(chǎn)潛力研究[J]. 草地學(xué)報(bào), 2013, 21(2): 236–242.

    [16] 羅瑞敏, 郭梁, 程積民. 近20a云霧山草地生產(chǎn)力對(duì)氣候變化的響應(yīng)[J]. 干旱區(qū)研究, 2018, 35(1): 77–84.

    [17] CHEN Tao, BAO Anming, Guli Jiapaer, et al. Disentan-gling the relative impacts of climate change and human activities on arid and semiarid grasslands in Central Asia during 1982-2015[J]. Science of the Total Environment, 2019, 653: 1311–1325.

    [18] 陳楠楠. 2000–2014年氣候變化背景下錫林郭勒盟草地氣候生產(chǎn)力的時(shí)空動(dòng)態(tài)分析[D]. 呼和浩特: 內(nèi)蒙古大學(xué), 2016: 15.

    [19] 劉黎明, 張鳳榮, 趙英偉. 我國草地資源生產(chǎn)潛力分析及其可持續(xù)利用對(duì)策[J]. 中國人口·資源與環(huán)境, 2002, 12(4): 100–105.

    [20] 孫慧蘭, 李衛(wèi)紅, 徐遠(yuǎn)杰. 新疆伊犁河流域牧草氣候生產(chǎn)潛力的時(shí)空變化特征分析[J]. 草業(yè)學(xué)報(bào), 2010, 19(6): 55–61.

    [21] 楊正禮, 楊改河. 中國高寒草地生產(chǎn)潛力與載畜量研究[J]. 資源科學(xué), 2000, 22(4): 72–77.

    [22] 龔相文, 李玉強(qiáng), 王旭洋, 等. 科爾沁地區(qū)植被凈初級(jí)生產(chǎn)力時(shí)空動(dòng)態(tài)特征[J]. 草業(yè)科學(xué), 2018, 35(10): 2306– 2317.

    [23] 陳國南. 用邁阿密模型測(cè)算我國生物生產(chǎn)量的初步嘗試[J]. 自然資源學(xué)報(bào), 1987, 2(3): 270–278.

    [24] 張偉科, 劉玉杰, 封志明. 內(nèi)蒙古天然草地生產(chǎn)潛力及其限制性研究[J]. 草地學(xué)報(bào), 2008, 16(6): 572–579.

    [25] 張憲洲. 我國自然植被凈第一性生產(chǎn)力的估算與分布[J]. 自然資源, 1993, 7(1): 15.

    [26] 朱玉果, 杜靈通, 謝應(yīng)忠, 等. 2000—2015 年寧夏草地凈初級(jí)生產(chǎn)力時(shí)空特征及其氣候響應(yīng)[J]. 生態(tài)學(xué)報(bào), 2019, 39(2): 518–529.

    [27] 周平, 武威, 王瑞, 等. 不同草地NPP估算模型對(duì)中國草地的模擬計(jì)算分析[J]. 草業(yè)科學(xué), 2018, 35(10): 2381– 2388.

    Climate change in arid regions of Northwest China and its impact on potential grassland productivity

    ZHANG Juan1, LI Jianping2,WANG Yutao1, ZHANG Yi1, JING Le1, LI Jianping1,3,*

    1. School of Agriculture, Ningxia University, Ningxia Yinchuan 750021, China 2. Ningxia Key Laboratory of Disaster Prevention and Reduction, Yinchuan 750021, China 3. National Key Laboratory Breeding Base for Restoration and Reconstruction of Degraded Ecosystems in Northwest China, Yinchuan 750021, China

    Based on the monthly temperature and precipitation data of 24 meteorological stations in Ningxia during the past 37 years (1981-2017), the Miami model and the Thornthwaite Memorial model were used to calculate the potential grassland productivity in Ningxia,and analyzed its spatial and temporal changes. The results showed that from the perspective of time changes, the average temperature in Ningxia continued to rise at a rate of 0.47°C every 10 years,and the precipitation first decreased and then increased; the annual precipitation increased by 40 mm from 2010 to 2017. From northwest to southeast in Ningxia, the phenology was changed from "warm and dry" to "cold and wet". Meanwhile, the potential climatic productivity of grassland in Ningxia gradually increased from northwest to southeast. The trend of potentialclimatic productivity of grassland was consistent with the trend of precipitation. There was a significant positive correlation between potential climatic productivity and precipitation (<0.01)each time period, and a negative correlation between potential climatic productivity and temperature. Therefore, precipitation was the limited factor for the potential climatic productivity of grassland in Ningxia, the largest potential grassland productivity was in the south-central part of Ningxia, and the climate change characterized by rising temperature has a weak impact on the potential climatic productivity of Ningxia grassland.

    potential grassland productivity; climate change; Ningxia; Miami model

    10.14108/j.cnki.1008-8873.2020.03.024

    S812

    A

    1008-8873(2020)03-182-11

    2019-06-23;

    2019-09-08

    國家自然科學(xué)基金(31660143)資助; 寧夏大學(xué)草學(xué)一流學(xué)科建設(shè)項(xiàng)目(NXYLXK2017A01)

    張娟(1993—), 女, 山西太原人, 碩士, 主要從事草地生態(tài)與模型構(gòu)建研究, E-mail: 18435121531@163.com

    李建平(1982—), 男, 陜西延安人, 博士, 副教授, 碩導(dǎo), 主要從事草地生態(tài)與模型構(gòu)建研究, E-mail: lijianpingsas@163.com

    張娟, 李劍萍, 王譽(yù)陶,等. 西北干旱地區(qū)氣候變化及其對(duì)草地生產(chǎn)潛力的影響[J]. 生態(tài)科學(xué), 2020, 39(3): 182–192.

    ZHANG Juan, LI Jianping, WANG Yutao, et al. Climate change in arid regions of Northwest China and its impact on potential grassland productivity [J]. Ecological Science, 2020, 39(3): 182–192.

    猜你喜歡
    生產(chǎn)潛力降水量寧夏
    繪制和閱讀降水量柱狀圖
    1978-2017年河南省氣候生產(chǎn)潛力時(shí)空演變及趨勢(shì)分析
    基于DSSAT模型和水分滿足率的的農(nóng)業(yè)耕種潛力時(shí)空評(píng)估新方式研究
    寧夏
    降水量是怎么算出來的
    啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
    1961~2015年貴州省氣候生產(chǎn)潛力估算及趨勢(shì)分析*
    1988—2017年呼和浩特市降水演變特征分析
    某化工企業(yè)清潔生產(chǎn)潛力分析
    寧夏煤電博弈
    能源(2017年5期)2017-07-06 09:25:50
    寧夏
    国产三级在线视频| 搡女人真爽免费视频火全软件 | 最近在线观看免费完整版| 欧美极品一区二区三区四区| 久久精品国产亚洲av涩爱 | 国产精品自产拍在线观看55亚洲| xxxwww97欧美| netflix在线观看网站| 一进一出抽搐gif免费好疼| 久久中文看片网| 色播亚洲综合网| 国产成人欧美在线观看| 搞女人的毛片| 最新中文字幕久久久久| 亚洲av第一区精品v没综合| 国内精品久久久久久久电影| 中文资源天堂在线| 亚洲精品一区av在线观看| 床上黄色一级片| 欧美成狂野欧美在线观看| 高清在线国产一区| 欧美日韩瑟瑟在线播放| 热99在线观看视频| 精品久久久久久久末码| 给我免费播放毛片高清在线观看| 久久久久精品国产欧美久久久| 中文资源天堂在线| 久久精品国产自在天天线| 偷拍熟女少妇极品色| 很黄的视频免费| 国产激情偷乱视频一区二区| 人人妻,人人澡人人爽秒播| 女生性感内裤真人,穿戴方法视频| 给我免费播放毛片高清在线观看| 99热这里只有是精品50| 全区人妻精品视频| 免费高清视频大片| 午夜激情欧美在线| 色综合婷婷激情| 日韩人妻高清精品专区| 国产精品精品国产色婷婷| 黄色片一级片一级黄色片| 可以在线观看毛片的网站| 国内精品美女久久久久久| 丰满人妻一区二区三区视频av | 天天一区二区日本电影三级| 国产高清激情床上av| 69av精品久久久久久| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 亚洲精品成人久久久久久| 18禁美女被吸乳视频| 最新在线观看一区二区三区| 99久久九九国产精品国产免费| eeuss影院久久| 国产高清有码在线观看视频| 国内久久婷婷六月综合欲色啪| 白带黄色成豆腐渣| 美女高潮喷水抽搐中文字幕| 可以在线观看的亚洲视频| 人妻丰满熟妇av一区二区三区| 欧美日韩黄片免| 美女高潮喷水抽搐中文字幕| 非洲黑人性xxxx精品又粗又长| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久亚洲av鲁大| 69人妻影院| 99久国产av精品| 天天一区二区日本电影三级| 神马国产精品三级电影在线观看| 男人舔奶头视频| 成熟少妇高潮喷水视频| 啦啦啦韩国在线观看视频| 欧美最新免费一区二区三区 | 成人无遮挡网站| 一边摸一边抽搐一进一小说| 欧美日韩福利视频一区二区| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 免费在线观看日本一区| 少妇的逼水好多| 天天一区二区日本电影三级| 欧美三级亚洲精品| 午夜免费成人在线视频| 男人舔女人下体高潮全视频| 一二三四社区在线视频社区8| 国产黄色小视频在线观看| 嫩草影院入口| 亚洲在线观看片| 两个人的视频大全免费| 偷拍熟女少妇极品色| 蜜桃亚洲精品一区二区三区| 久久香蕉精品热| 中文字幕人妻丝袜一区二区| 最新在线观看一区二区三区| 久久国产精品影院| 国产一区二区激情短视频| 非洲黑人性xxxx精品又粗又长| 91九色精品人成在线观看| 亚洲国产精品999在线| 亚洲欧美日韩卡通动漫| 国产精品一区二区免费欧美| av中文乱码字幕在线| 我要搜黄色片| 精品人妻一区二区三区麻豆 | 天天躁日日操中文字幕| 午夜精品一区二区三区免费看| 一进一出抽搐gif免费好疼| 午夜福利免费观看在线| 亚洲av电影不卡..在线观看| 国产精品99久久99久久久不卡| 久久精品国产清高在天天线| 最近在线观看免费完整版| 人人妻人人看人人澡| 男女午夜视频在线观看| 久久精品91无色码中文字幕| 午夜精品一区二区三区免费看| 精品电影一区二区在线| 亚洲七黄色美女视频| 法律面前人人平等表现在哪些方面| 婷婷精品国产亚洲av| 色吧在线观看| 日本黄色片子视频| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 国产精品久久久久久久久免 | 国内少妇人妻偷人精品xxx网站| 久久久国产精品麻豆| 黄色片一级片一级黄色片| 国产免费男女视频| 亚洲18禁久久av| 热99re8久久精品国产| 老司机深夜福利视频在线观看| 亚洲片人在线观看| 操出白浆在线播放| 国产毛片a区久久久久| 久久久成人免费电影| 欧美乱妇无乱码| 免费一级毛片在线播放高清视频| 3wmmmm亚洲av在线观看| 最近在线观看免费完整版| 9191精品国产免费久久| 国产不卡一卡二| 最后的刺客免费高清国语| 国产成人av激情在线播放| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 国产伦精品一区二区三区四那| 日本与韩国留学比较| 国产精品久久久久久人妻精品电影| 午夜久久久久精精品| 亚洲成人久久性| 嫩草影院精品99| 免费av毛片视频| 国产亚洲av嫩草精品影院| 一级毛片高清免费大全| 草草在线视频免费看| 欧美成狂野欧美在线观看| 在线国产一区二区在线| 宅男免费午夜| 99热6这里只有精品| 欧美精品啪啪一区二区三区| 国产成人福利小说| 亚洲男人的天堂狠狠| 精品一区二区三区人妻视频| 岛国在线免费视频观看| 国产激情欧美一区二区| www日本黄色视频网| 久久性视频一级片| 国内精品一区二区在线观看| 亚洲精品粉嫩美女一区| 欧美中文综合在线视频| 欧美+亚洲+日韩+国产| 女人十人毛片免费观看3o分钟| 久久人妻av系列| 大型黄色视频在线免费观看| 国产精品女同一区二区软件 | 97超级碰碰碰精品色视频在线观看| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 国产成人欧美在线观看| 欧美一区二区精品小视频在线| 我要搜黄色片| 精品国内亚洲2022精品成人| av天堂在线播放| 啪啪无遮挡十八禁网站| 色综合婷婷激情| 男女做爰动态图高潮gif福利片| 精品一区二区三区视频在线 | 给我免费播放毛片高清在线观看| 色老头精品视频在线观看| 亚洲精品色激情综合| 国内精品久久久久久久电影| 国产一区二区三区视频了| 少妇裸体淫交视频免费看高清| 精品人妻一区二区三区麻豆 | 精品一区二区三区视频在线 | 亚洲av免费高清在线观看| 9191精品国产免费久久| 老汉色av国产亚洲站长工具| 日韩欧美在线二视频| 在线观看av片永久免费下载| 久久久久国内视频| 岛国在线观看网站| 亚洲熟妇中文字幕五十中出| www.熟女人妻精品国产| 91av网一区二区| 国产亚洲av嫩草精品影院| 国产精品99久久久久久久久| 亚洲av美国av| 日韩高清综合在线| 成年女人毛片免费观看观看9| 真人做人爱边吃奶动态| 国产一级毛片七仙女欲春2| 91麻豆精品激情在线观看国产| 午夜亚洲福利在线播放| 97超级碰碰碰精品色视频在线观看| 少妇人妻精品综合一区二区 | 中亚洲国语对白在线视频| 国产激情欧美一区二区| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 我的老师免费观看完整版| av视频在线观看入口| av黄色大香蕉| 欧美bdsm另类| 国产日本99.免费观看| 女同久久另类99精品国产91| 18禁美女被吸乳视频| 一级毛片女人18水好多| 久久精品国产清高在天天线| 亚洲人成网站在线播| ponron亚洲| 91久久精品国产一区二区成人 | 亚洲av美国av| 成人18禁在线播放| 97人妻精品一区二区三区麻豆| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| 免费大片18禁| 国产伦精品一区二区三区四那| 99久久精品热视频| 亚洲第一欧美日韩一区二区三区| 国产精品影院久久| 亚洲成人久久爱视频| 欧美日韩精品网址| 91在线精品国自产拍蜜月 | 亚洲七黄色美女视频| 丁香六月欧美| 亚洲18禁久久av| 亚洲avbb在线观看| 在线视频色国产色| 免费在线观看亚洲国产| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 国产精品1区2区在线观看.| 国产精品爽爽va在线观看网站| 亚洲天堂国产精品一区在线| 婷婷六月久久综合丁香| 欧美成人免费av一区二区三区| 国内少妇人妻偷人精品xxx网站| 久久精品91无色码中文字幕| 国产免费男女视频| 99国产综合亚洲精品| 亚洲在线自拍视频| 在线免费观看的www视频| 免费看日本二区| 每晚都被弄得嗷嗷叫到高潮| 男插女下体视频免费在线播放| 日韩国内少妇激情av| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟人妻熟丝袜美 | 两个人看的免费小视频| 久久99热这里只有精品18| e午夜精品久久久久久久| 国产一级毛片七仙女欲春2| 丝袜美腿在线中文| 又黄又爽又免费观看的视频| 免费大片18禁| 欧洲精品卡2卡3卡4卡5卡区| 99精品在免费线老司机午夜| 身体一侧抽搐| 国产精品久久视频播放| 国产精品国产高清国产av| 综合色av麻豆| 丁香六月欧美| 亚洲激情在线av| 十八禁人妻一区二区| 啦啦啦韩国在线观看视频| 久久99热这里只有精品18| 欧美日韩一级在线毛片| 搡老熟女国产l中国老女人| 久久精品91蜜桃| 成人特级av手机在线观看| 欧美中文综合在线视频| 少妇熟女aⅴ在线视频| 久久婷婷人人爽人人干人人爱| 欧美激情在线99| 真实男女啪啪啪动态图| 久久精品国产自在天天线| 桃色一区二区三区在线观看| 激情在线观看视频在线高清| 美女免费视频网站| 丝袜美腿在线中文| 脱女人内裤的视频| 一个人看的www免费观看视频| 女人高潮潮喷娇喘18禁视频| 欧美成人a在线观看| 久久伊人香网站| 欧美区成人在线视频| 哪里可以看免费的av片| av片东京热男人的天堂| 国产精品98久久久久久宅男小说| 美女黄网站色视频| 狂野欧美白嫩少妇大欣赏| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| x7x7x7水蜜桃| 一个人看的www免费观看视频| 亚洲av成人av| 老司机午夜福利在线观看视频| 亚洲专区国产一区二区| 九色成人免费人妻av| 欧美色欧美亚洲另类二区| 精品福利观看| 一个人免费在线观看电影| 18+在线观看网站| 成年女人永久免费观看视频| 天堂av国产一区二区熟女人妻| 亚洲熟妇中文字幕五十中出| 中文字幕久久专区| 麻豆国产97在线/欧美| 高潮久久久久久久久久久不卡| 午夜免费观看网址| 久久久久精品国产欧美久久久| 国内揄拍国产精品人妻在线| svipshipincom国产片| 国产不卡一卡二| 亚洲欧美精品综合久久99| 亚洲欧美日韩无卡精品| 蜜桃亚洲精品一区二区三区| 午夜福利欧美成人| 校园春色视频在线观看| av欧美777| 午夜福利成人在线免费观看| 国产一区二区在线观看日韩 | 久久久久久九九精品二区国产| 禁无遮挡网站| 悠悠久久av| 最近在线观看免费完整版| 色哟哟哟哟哟哟| 桃红色精品国产亚洲av| 身体一侧抽搐| 国产美女午夜福利| av天堂在线播放| x7x7x7水蜜桃| 亚洲自拍偷在线| 久久久久性生活片| 无遮挡黄片免费观看| 欧美3d第一页| 3wmmmm亚洲av在线观看| 亚洲国产精品久久男人天堂| 久久精品国产清高在天天线| 国产成人a区在线观看| 18禁国产床啪视频网站| 午夜精品一区二区三区免费看| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 国产欧美日韩精品亚洲av| 亚洲av成人精品一区久久| 国产精品,欧美在线| 欧美大码av| 亚洲,欧美精品.| 精品久久久久久成人av| 精品人妻1区二区| 哪里可以看免费的av片| 久久婷婷人人爽人人干人人爱| 男插女下体视频免费在线播放| 一本综合久久免费| 亚洲av美国av| 久久久久精品国产欧美久久久| av天堂在线播放| 精品熟女少妇八av免费久了| 网址你懂的国产日韩在线| 亚洲精品美女久久久久99蜜臀| 日本黄色片子视频| 又紧又爽又黄一区二区| 午夜福利免费观看在线| 级片在线观看| 日本五十路高清| 99riav亚洲国产免费| 搞女人的毛片| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品青青久久久久久| 在线观看午夜福利视频| 国产亚洲精品综合一区在线观看| 国产伦精品一区二区三区视频9 | 校园春色视频在线观看| 欧美又色又爽又黄视频| 1024手机看黄色片| 在线天堂最新版资源| 欧美成人a在线观看| 亚洲最大成人中文| or卡值多少钱| 欧美一级a爱片免费观看看| 久久婷婷人人爽人人干人人爱| 又爽又黄无遮挡网站| 日韩欧美 国产精品| 老鸭窝网址在线观看| 波多野结衣巨乳人妻| 亚洲av免费在线观看| 国产精品久久久久久久久免 | 网址你懂的国产日韩在线| 国产精品女同一区二区软件 | 精品欧美国产一区二区三| h日本视频在线播放| 深爱激情五月婷婷| 丁香六月欧美| 美女黄网站色视频| 91麻豆av在线| 日韩人妻高清精品专区| 99riav亚洲国产免费| 亚洲最大成人中文| 国产亚洲欧美在线一区二区| 国产 一区 欧美 日韩| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 国产一区二区三区视频了| 国产精品久久久久久人妻精品电影| 欧美+日韩+精品| 露出奶头的视频| 国产精品一区二区免费欧美| 最新中文字幕久久久久| 岛国在线免费视频观看| 国产激情偷乱视频一区二区| 天堂动漫精品| 黄色丝袜av网址大全| 午夜福利18| 日韩亚洲欧美综合| 国产精品98久久久久久宅男小说| av欧美777| a在线观看视频网站| 国产视频一区二区在线看| 亚洲无线在线观看| 国产精品亚洲av一区麻豆| 国产av一区在线观看免费| 亚洲一区二区三区色噜噜| 国产精品女同一区二区软件 | 日本 av在线| 日本撒尿小便嘘嘘汇集6| 国产午夜福利久久久久久| 亚洲国产精品成人综合色| 香蕉久久夜色| 小蜜桃在线观看免费完整版高清| 久久中文看片网| 狂野欧美激情性xxxx| 国产爱豆传媒在线观看| 男女视频在线观看网站免费| 在线播放无遮挡| 有码 亚洲区| 高潮久久久久久久久久久不卡| 一区福利在线观看| 麻豆成人午夜福利视频| 国语自产精品视频在线第100页| 国产精品爽爽va在线观看网站| 亚洲在线观看片| 美女黄网站色视频| 日本三级黄在线观看| 美女大奶头视频| 91麻豆av在线| 国产真人三级小视频在线观看| 国产成人欧美在线观看| 精品人妻偷拍中文字幕| 国产成人影院久久av| 757午夜福利合集在线观看| 九九热线精品视视频播放| 色尼玛亚洲综合影院| 色噜噜av男人的天堂激情| 国产v大片淫在线免费观看| 午夜福利18| 成人鲁丝片一二三区免费| 国产真实伦视频高清在线观看 | 美女高潮喷水抽搐中文字幕| 久久久久精品国产欧美久久久| 国产一区二区在线观看日韩 | 岛国在线观看网站| 日本a在线网址| 99久久99久久久精品蜜桃| 亚洲熟妇熟女久久| 观看免费一级毛片| 亚洲片人在线观看| 精品久久久久久成人av| 免费看a级黄色片| tocl精华| 国产乱人伦免费视频| 成人无遮挡网站| 午夜两性在线视频| 午夜老司机福利剧场| 欧美区成人在线视频| 国产亚洲av嫩草精品影院| 国产一区二区三区视频了| 免费观看的影片在线观看| 久久香蕉精品热| 久久性视频一级片| 色哟哟哟哟哟哟| 欧美色视频一区免费| 国产一区二区三区视频了| 国产成人aa在线观看| 日本五十路高清| 日本撒尿小便嘘嘘汇集6| 又粗又爽又猛毛片免费看| 操出白浆在线播放| 欧美日韩黄片免| 午夜免费激情av| 日本黄色视频三级网站网址| 岛国在线免费视频观看| 欧美绝顶高潮抽搐喷水| 在线观看美女被高潮喷水网站 | 国产精品 欧美亚洲| 精品国产美女av久久久久小说| svipshipincom国产片| 91字幕亚洲| 99精品久久久久人妻精品| 亚洲成人久久性| 免费在线观看亚洲国产| 色在线成人网| 国产亚洲欧美98| 亚洲成人精品中文字幕电影| 夜夜夜夜夜久久久久| 国产三级中文精品| 午夜激情福利司机影院| 又粗又爽又猛毛片免费看| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 国产97色在线日韩免费| 亚洲美女黄片视频| 首页视频小说图片口味搜索| 91av网一区二区| 亚洲不卡免费看| 欧美黄色片欧美黄色片| 国产伦精品一区二区三区四那| 香蕉丝袜av| 一个人看的www免费观看视频| 午夜老司机福利剧场| 国产精品99久久99久久久不卡| 亚洲精华国产精华精| 黄色日韩在线| 国模一区二区三区四区视频| 亚洲国产高清在线一区二区三| 日韩亚洲欧美综合| 嫩草影视91久久| 亚洲欧美激情综合另类| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 国产三级中文精品| 午夜激情福利司机影院| 婷婷精品国产亚洲av在线| 99久久精品国产亚洲精品| 免费看a级黄色片| 亚洲第一欧美日韩一区二区三区| 一a级毛片在线观看| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产| 观看免费一级毛片| 欧美黄色片欧美黄色片| 日本一本二区三区精品| 天堂√8在线中文| 免费在线观看日本一区| 中文字幕人妻丝袜一区二区| 中文字幕高清在线视频| 两个人视频免费观看高清| 午夜两性在线视频| 亚洲人与动物交配视频| 亚洲av成人不卡在线观看播放网| 很黄的视频免费| 丁香欧美五月| 在线播放国产精品三级| 国产精品一区二区三区四区久久| 亚洲av中文字字幕乱码综合| 极品教师在线免费播放| 国产精品国产高清国产av| 国产成人影院久久av| 久久久久久人人人人人| 桃色一区二区三区在线观看| 一区福利在线观看| 两个人看的免费小视频| a级毛片a级免费在线| 757午夜福利合集在线观看| 亚洲精品久久国产高清桃花| 我的老师免费观看完整版| 久久久久性生活片| 热99在线观看视频| 亚洲在线观看片| 听说在线观看完整版免费高清| 亚洲av五月六月丁香网| 一边摸一边抽搐一进一小说| 好男人电影高清在线观看| 高清在线国产一区| 一夜夜www| 嫩草影视91久久| 亚洲人成网站高清观看| 狠狠狠狠99中文字幕| 日韩欧美国产一区二区入口| 亚洲在线自拍视频| 99久久无色码亚洲精品果冻| 国产三级在线视频| 成人欧美大片| 午夜老司机福利剧场| 啦啦啦免费观看视频1| 亚洲国产欧美网| 给我免费播放毛片高清在线观看| 国内揄拍国产精品人妻在线| 一本综合久久免费| 国产精品美女特级片免费视频播放器|