• <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
    寧夏
    久久婷婷成人综合色麻豆| 成年女人毛片免费观看观看9| 人人妻,人人澡人人爽秒播| 国产一级毛片七仙女欲春2 | 女同久久另类99精品国产91| 国产成人啪精品午夜网站| 亚洲精品久久国产高清桃花| 国产免费av片在线观看野外av| 一进一出好大好爽视频| 久久亚洲真实| 亚洲欧美日韩高清在线视频| 91av网站免费观看| 亚洲美女黄片视频| 亚洲国产欧美网| 国产麻豆成人av免费视频| 国产不卡一卡二| 欧美日本视频| av片东京热男人的天堂| 亚洲自拍偷在线| 亚洲国产欧美网| 国产精品久久电影中文字幕| 男人的好看免费观看在线视频 | 一个人免费在线观看的高清视频| 久久狼人影院| 亚洲国产高清在线一区二区三 | 嫩草影视91久久| 国产亚洲欧美精品永久| 天堂√8在线中文| 人人妻人人澡欧美一区二区| 欧美日韩精品网址| 欧美一级毛片孕妇| av超薄肉色丝袜交足视频| 人成视频在线观看免费观看| 欧美激情高清一区二区三区| av有码第一页| 99国产综合亚洲精品| 69av精品久久久久久| av超薄肉色丝袜交足视频| 亚洲欧美一区二区三区黑人| 国产高清有码在线观看视频 | 热re99久久国产66热| xxxwww97欧美| 99热只有精品国产| 日韩欧美三级三区| 90打野战视频偷拍视频| 亚洲激情在线av| 1024手机看黄色片| 精品一区二区三区四区五区乱码| 免费一级毛片在线播放高清视频| 91成年电影在线观看| 村上凉子中文字幕在线| 高清在线国产一区| or卡值多少钱| 777久久人妻少妇嫩草av网站| 给我免费播放毛片高清在线观看| 黄色成人免费大全| 午夜福利视频1000在线观看| 亚洲精品在线观看二区| 日日摸夜夜添夜夜添小说| 久久人妻福利社区极品人妻图片| 欧美在线黄色| 亚洲一区二区三区不卡视频| 校园春色视频在线观看| 99国产精品一区二区蜜桃av| 午夜影院日韩av| 老司机在亚洲福利影院| 黄色丝袜av网址大全| 巨乳人妻的诱惑在线观看| 麻豆一二三区av精品| 久久久精品国产亚洲av高清涩受| 日本免费a在线| 国产高清有码在线观看视频 | 波多野结衣高清无吗| 一本大道久久a久久精品| 久久精品亚洲精品国产色婷小说| 亚洲av五月六月丁香网| 亚洲国产欧美网| 午夜激情av网站| 久久香蕉精品热| 日韩欧美国产一区二区入口| 一区福利在线观看| 黄色视频不卡| 两个人视频免费观看高清| 国产三级黄色录像| 十八禁网站免费在线| 两个人视频免费观看高清| 两性夫妻黄色片| 欧美黄色淫秽网站| 高清在线国产一区| 国产aⅴ精品一区二区三区波| 日韩大尺度精品在线看网址| 国产熟女xx| 在线观看66精品国产| 999精品在线视频| 啦啦啦免费观看视频1| 久久精品夜夜夜夜夜久久蜜豆 | 国产成人影院久久av| 黄频高清免费视频| 午夜a级毛片| 国语自产精品视频在线第100页| 欧美性长视频在线观看| 精品欧美一区二区三区在线| 69av精品久久久久久| 一级片免费观看大全| 免费在线观看完整版高清| 国产精品,欧美在线| 久久中文看片网| 日日干狠狠操夜夜爽| 精品国产乱码久久久久久男人| 两性夫妻黄色片| 校园春色视频在线观看| 人妻丰满熟妇av一区二区三区| netflix在线观看网站| 热re99久久国产66热| 国产精品免费视频内射| 日韩免费av在线播放| 美女大奶头视频| 19禁男女啪啪无遮挡网站| 国产精品 国内视频| 1024手机看黄色片| www.精华液| 亚洲人成网站高清观看| 色综合欧美亚洲国产小说| 亚洲av电影不卡..在线观看| 在线播放国产精品三级| 淫秽高清视频在线观看| 国产97色在线日韩免费| 99热只有精品国产| 日韩精品中文字幕看吧| 国产高清videossex| 欧美人与性动交α欧美精品济南到| 一区福利在线观看| 白带黄色成豆腐渣| 久久久国产成人免费| 中文字幕人妻丝袜一区二区| 免费高清视频大片| 深夜精品福利| 老汉色av国产亚洲站长工具| 久久久久久久久久黄片| 久久精品aⅴ一区二区三区四区| 日本五十路高清| 少妇粗大呻吟视频| 国产精品亚洲一级av第二区| 老熟妇乱子伦视频在线观看| 日韩精品中文字幕看吧| 妹子高潮喷水视频| 日韩成人在线观看一区二区三区| 宅男免费午夜| 亚洲狠狠婷婷综合久久图片| 午夜福利视频1000在线观看| 亚洲欧美精品综合久久99| xxx96com| 欧美不卡视频在线免费观看 | 岛国在线观看网站| 伊人久久大香线蕉亚洲五| 婷婷精品国产亚洲av在线| 丝袜在线中文字幕| 一区二区三区高清视频在线| 国产又色又爽无遮挡免费看| 日本精品一区二区三区蜜桃| 久久久水蜜桃国产精品网| 国产精品av久久久久免费| 俺也久久电影网| 亚洲欧美一区二区三区黑人| 制服诱惑二区| 久久草成人影院| 啦啦啦观看免费观看视频高清| 国产熟女xx| 国产精品亚洲一级av第二区| 国产一区在线观看成人免费| 欧美性猛交╳xxx乱大交人| 中亚洲国语对白在线视频| 性欧美人与动物交配| 国产成人精品久久二区二区免费| 少妇 在线观看| av天堂在线播放| 欧美激情 高清一区二区三区| 51午夜福利影视在线观看| 欧美乱妇无乱码| 精品久久蜜臀av无| 在线国产一区二区在线| 久久久久国产精品人妻aⅴ院| 桃色一区二区三区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 在线观看免费午夜福利视频| 日韩欧美在线二视频| 久9热在线精品视频| 欧美日韩福利视频一区二区| www.www免费av| 在线天堂中文资源库| 亚洲国产精品成人综合色| 亚洲精品美女久久久久99蜜臀| 成人国产综合亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产精品999在线| 欧美黑人欧美精品刺激| 51午夜福利影视在线观看| 91麻豆av在线| 欧美性长视频在线观看| av中文乱码字幕在线| 国产精品免费视频内射| 色尼玛亚洲综合影院| 夜夜夜夜夜久久久久| 欧美黄色淫秽网站| 中文在线观看免费www的网站 | 亚洲一区高清亚洲精品| 久久久国产成人免费| 久久久国产成人免费| 亚洲男人的天堂狠狠| 亚洲av电影不卡..在线观看| 高潮久久久久久久久久久不卡| 黄片播放在线免费| 国产精品1区2区在线观看.| av中文乱码字幕在线| 亚洲成人精品中文字幕电影| 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 精品久久久久久久人妻蜜臀av| 国产亚洲精品一区二区www| 在线观看www视频免费| 国内少妇人妻偷人精品xxx网站 | 亚洲男人的天堂狠狠| 欧美绝顶高潮抽搐喷水| 一区二区三区高清视频在线| 亚洲av美国av| 又大又爽又粗| 法律面前人人平等表现在哪些方面| 91老司机精品| 少妇熟女aⅴ在线视频| 悠悠久久av| 亚洲自拍偷在线| 亚洲中文av在线| 欧美乱妇无乱码| 久久人妻av系列| 在线观看免费视频日本深夜| 国产一区在线观看成人免费| 90打野战视频偷拍视频| 久久人妻福利社区极品人妻图片| 91国产中文字幕| 久久午夜综合久久蜜桃| 一级黄色大片毛片| 国产真人三级小视频在线观看| 国产精品一区二区三区四区久久 | 欧美黄色淫秽网站| 无限看片的www在线观看| 成人国产一区最新在线观看| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 老司机午夜十八禁免费视频| 亚洲中文日韩欧美视频| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 禁无遮挡网站| 免费在线观看日本一区| 亚洲电影在线观看av| 国产免费av片在线观看野外av| 久久精品国产亚洲av高清一级| 成人免费观看视频高清| 大型av网站在线播放| 日韩欧美免费精品| 久久九九热精品免费| av福利片在线| 成人三级做爰电影| 91麻豆av在线| 老司机深夜福利视频在线观看| 亚洲一区高清亚洲精品| 黑人操中国人逼视频| 欧美 亚洲 国产 日韩一| 亚洲 欧美 日韩 在线 免费| 伊人久久大香线蕉亚洲五| 亚洲男人天堂网一区| 午夜福利视频1000在线观看| 国产成人av激情在线播放| 妹子高潮喷水视频| 99精品欧美一区二区三区四区| 欧美人与性动交α欧美精品济南到| 中文字幕高清在线视频| 老熟妇仑乱视频hdxx| 欧美激情高清一区二区三区| 中文字幕精品亚洲无线码一区 | 日韩免费av在线播放| 在线天堂中文资源库| 女警被强在线播放| 男男h啪啪无遮挡| 亚洲av成人不卡在线观看播放网| 成人免费观看视频高清| 18禁黄网站禁片午夜丰满| 99国产精品99久久久久| 亚洲全国av大片| 黑丝袜美女国产一区| 三级毛片av免费| www.www免费av| 少妇的丰满在线观看| 国产99久久九九免费精品| 日韩精品中文字幕看吧| 老熟妇仑乱视频hdxx| 欧美激情 高清一区二区三区| 日本免费a在线| 欧美成狂野欧美在线观看| 久久久国产成人精品二区| 国产久久久一区二区三区| 91九色精品人成在线观看| 日本免费一区二区三区高清不卡| 久久久久久久久久黄片| 午夜福利在线观看吧| 美女扒开内裤让男人捅视频| 色综合婷婷激情| 亚洲三区欧美一区| 国产91精品成人一区二区三区| 日韩欧美免费精品| 国产成人av教育| 99精品在免费线老司机午夜| 757午夜福利合集在线观看| 午夜免费成人在线视频| 午夜福利欧美成人| 精品久久久久久久毛片微露脸| 一边摸一边抽搐一进一小说| 观看免费一级毛片| 一区二区三区高清视频在线| 成在线人永久免费视频| 久久草成人影院| 侵犯人妻中文字幕一二三四区| 亚洲国产欧洲综合997久久, | 午夜两性在线视频| 欧美激情极品国产一区二区三区| 国产极品粉嫩免费观看在线| 人人妻人人看人人澡| 波多野结衣高清无吗| 搡老妇女老女人老熟妇| 人妻久久中文字幕网| 国产熟女xx| 精品人妻1区二区| 91av网站免费观看| 熟妇人妻久久中文字幕3abv| 欧美+亚洲+日韩+国产| 听说在线观看完整版免费高清| 国产色视频综合| 成人一区二区视频在线观看| 国语自产精品视频在线第100页| 成人三级黄色视频| 免费高清在线观看日韩| x7x7x7水蜜桃| 不卡av一区二区三区| 一区二区日韩欧美中文字幕| 中国美女看黄片| 午夜久久久久精精品| 亚洲电影在线观看av| 国产精品日韩av在线免费观看| 夜夜夜夜夜久久久久| 不卡av一区二区三区| 免费电影在线观看免费观看| a在线观看视频网站| 真人做人爱边吃奶动态| 99久久久亚洲精品蜜臀av| 成人午夜高清在线视频 | 天天躁夜夜躁狠狠躁躁| 免费女性裸体啪啪无遮挡网站| 精品欧美一区二区三区在线| 亚洲成人免费电影在线观看| 18禁国产床啪视频网站| 一级作爱视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 国产又色又爽无遮挡免费看| 亚洲成av人片免费观看| 一本综合久久免费| 欧美黄色片欧美黄色片| 麻豆成人午夜福利视频| 欧美精品啪啪一区二区三区| 99热只有精品国产| 99riav亚洲国产免费| 中出人妻视频一区二区| 中文亚洲av片在线观看爽| 亚洲国产毛片av蜜桃av| 精品国产乱码久久久久久男人| 一级作爱视频免费观看| 成人欧美大片| 国产一区二区三区视频了| 男女做爰动态图高潮gif福利片| 欧美黑人巨大hd| 欧美丝袜亚洲另类 | 色综合婷婷激情| 一a级毛片在线观看| 嫁个100分男人电影在线观看| 91九色精品人成在线观看| 久久久久久久久久黄片| 首页视频小说图片口味搜索| 国产爱豆传媒在线观看 | 午夜精品在线福利| 一本大道久久a久久精品| 免费电影在线观看免费观看| 国产区一区二久久| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品999在线| 最近在线观看免费完整版| 9191精品国产免费久久| 午夜福利欧美成人| 免费在线观看黄色视频的| 少妇熟女aⅴ在线视频| 十八禁网站免费在线| 国产视频内射| 熟妇人妻久久中文字幕3abv| 久热这里只有精品99| 在线观看免费午夜福利视频| 老熟妇乱子伦视频在线观看| 成人亚洲精品一区在线观看| 91麻豆av在线| 国产99久久九九免费精品| 亚洲精品中文字幕在线视频| 露出奶头的视频| 久久这里只有精品19| 欧美激情久久久久久爽电影| 日韩 欧美 亚洲 中文字幕| 嫩草影院精品99| 狂野欧美激情性xxxx| 白带黄色成豆腐渣| 他把我摸到了高潮在线观看| 国产欧美日韩一区二区精品| 怎么达到女性高潮| 午夜视频精品福利| 免费观看人在逋| 女人被狂操c到高潮| 欧美精品亚洲一区二区| 99热只有精品国产| 精品一区二区三区四区五区乱码| 欧美性长视频在线观看| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 国产免费男女视频| 日韩欧美三级三区| www.精华液| 波多野结衣av一区二区av| 18禁黄网站禁片午夜丰满| 亚洲 欧美 日韩 在线 免费| 亚洲aⅴ乱码一区二区在线播放 | 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 啦啦啦观看免费观看视频高清| 宅男免费午夜| 亚洲成a人片在线一区二区| 精品少妇一区二区三区视频日本电影| 亚洲五月天丁香| 国产私拍福利视频在线观看| 久99久视频精品免费| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| aaaaa片日本免费| 亚洲av成人一区二区三| 午夜福利高清视频| 视频在线观看一区二区三区| 超碰成人久久| 国产精品爽爽va在线观看网站 | 免费无遮挡裸体视频| 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址| 男女床上黄色一级片免费看| 一本精品99久久精品77| 亚洲专区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 黄网站色视频无遮挡免费观看| 欧美黑人巨大hd| 女性被躁到高潮视频| 色综合站精品国产| 最新美女视频免费是黄的| 亚洲自拍偷在线| 欧美日韩福利视频一区二区| 欧美日韩精品网址| 在线免费观看的www视频| 高清在线国产一区| 99re在线观看精品视频| 最新在线观看一区二区三区| 久久久久久九九精品二区国产 | 亚洲精品av麻豆狂野| 国产av一区二区精品久久| 夜夜夜夜夜久久久久| 日韩免费av在线播放| 国产成人精品无人区| 韩国av一区二区三区四区| 欧美激情 高清一区二区三区| 日本在线视频免费播放| 少妇粗大呻吟视频| 亚洲色图 男人天堂 中文字幕| 不卡一级毛片| 免费观看人在逋| 欧美乱码精品一区二区三区| 国产av不卡久久| xxxwww97欧美| 亚洲精品国产区一区二| 久久久国产精品麻豆| 亚洲第一电影网av| 久久香蕉精品热| 一个人免费在线观看的高清视频| 国产av一区在线观看免费| 国产亚洲精品第一综合不卡| 这个男人来自地球电影免费观看| 中文资源天堂在线| 亚洲第一av免费看| 久久九九热精品免费| 天天一区二区日本电影三级| 欧美黄色片欧美黄色片| 国产久久久一区二区三区| 婷婷精品国产亚洲av在线| 18禁裸乳无遮挡免费网站照片 | 国产精品综合久久久久久久免费| 中国美女看黄片| 久久久久免费精品人妻一区二区 | 亚洲av熟女| 婷婷丁香在线五月| 国产亚洲精品久久久久久毛片| 深夜精品福利| 欧美黄色淫秽网站| 国产爱豆传媒在线观看 | 久久香蕉激情| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全免费视频| 日韩国内少妇激情av| 欧美性猛交╳xxx乱大交人| 国产亚洲精品综合一区在线观看 | 中文字幕人成人乱码亚洲影| 女人被狂操c到高潮| 国产野战对白在线观看| bbb黄色大片| 午夜福利在线在线| 午夜老司机福利片| 神马国产精品三级电影在线观看 | av片东京热男人的天堂| 久久久久久九九精品二区国产 | 欧美成人性av电影在线观看| 精品国产国语对白av| 亚洲无线在线观看| √禁漫天堂资源中文www| 久久精品国产清高在天天线| 国内精品久久久久精免费| 欧美黄色片欧美黄色片| 人妻丰满熟妇av一区二区三区| 岛国视频午夜一区免费看| 欧美一级a爱片免费观看看 | 国产不卡一卡二| 女人被狂操c到高潮| 久久国产精品男人的天堂亚洲| 狂野欧美激情性xxxx| 国产亚洲精品第一综合不卡| 操出白浆在线播放| 男人舔女人的私密视频| 国产精华一区二区三区| 亚洲精品国产区一区二| a在线观看视频网站| 老汉色∧v一级毛片| 夜夜爽天天搞| 99久久99久久久精品蜜桃| 国产精品九九99| 日日夜夜操网爽| АⅤ资源中文在线天堂| 成人18禁高潮啪啪吃奶动态图| 午夜福利欧美成人| 18禁黄网站禁片午夜丰满| av天堂在线播放| 免费在线观看亚洲国产| 国产真人三级小视频在线观看| 婷婷丁香在线五月| 日韩欧美一区视频在线观看| 这个男人来自地球电影免费观看| 在线观看66精品国产| 午夜福利在线观看吧| 成人国产综合亚洲| av视频在线观看入口| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区黑人| 亚洲国产中文字幕在线视频| 国产精品香港三级国产av潘金莲| 很黄的视频免费| 男人舔女人的私密视频| 后天国语完整版免费观看| 精品一区二区三区视频在线观看免费| 丝袜人妻中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 国产真实乱freesex| 亚洲av中文字字幕乱码综合 | 国产精品香港三级国产av潘金莲| 久久草成人影院| 99久久久亚洲精品蜜臀av| 精品久久久久久久末码| 女人爽到高潮嗷嗷叫在线视频| 90打野战视频偷拍视频| 亚洲中文日韩欧美视频| 极品教师在线免费播放| www日本在线高清视频| 亚洲av日韩精品久久久久久密| 91麻豆av在线| 亚洲国产看品久久| 中出人妻视频一区二区| 国内少妇人妻偷人精品xxx网站 | 国产一级毛片七仙女欲春2 | 每晚都被弄得嗷嗷叫到高潮| 18禁黄网站禁片午夜丰满| 国产av在哪里看| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 午夜久久久久精精品| 日韩三级视频一区二区三区| 成人欧美大片| 黄频高清免费视频| 免费无遮挡裸体视频| cao死你这个sao货| 亚洲成人免费电影在线观看| 久久亚洲真实| 少妇粗大呻吟视频| 精品欧美国产一区二区三| 亚洲成人免费电影在线观看| 久久久国产成人免费| 97碰自拍视频| 18禁黄网站禁片午夜丰满| 亚洲真实伦在线观看| 精品午夜福利视频在线观看一区|