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

    基于遙感與GIS的臨淄區(qū)土地利用變化與生態(tài)環(huán)境效應研究

    2020-06-15 06:38齊學蕾
    綠色科技 2020年6期
    關鍵詞:臨淄區(qū)貢獻率林地

    齊學蕾,常 軍

    (山東師范大學,山東 濟南 250300)

    1 引言

    土地利用和土地覆被變化可以直觀反映出一個地區(qū)土地資源的利用方式[1]。全球氣候變化和人類活動的加速導致了土地利用和土地覆被發(fā)生變化,從而引起了區(qū)域社會、經(jīng)濟和環(huán)境的變化,也引發(fā)了更多自然現(xiàn)象和生態(tài)過程的變化,如荒漠化、水環(huán)境變化以及植被覆蓋變化等。因此,雖然土地利用與覆被變化(LUCC)研究已經(jīng)過去20年,但仍然是全球環(huán)境變化研究的核心領域[2]。土地利用變化及其生態(tài)環(huán)境響應是地理研究的重點,且土地利用變化和生態(tài)環(huán)境響應與我國的總體發(fā)展格局密切相關[3]。因而,對區(qū)域土地利用的外部功能隨時間變化進行研究可以從區(qū)域發(fā)展的功能角度分析土地利用變化和區(qū)域生態(tài)環(huán)境響應的關系。以山東省淄博市臨淄區(qū)為例,利用2008年、2012年以及2018年三期土地利用數(shù)據(jù),探討了臨淄區(qū)10年間土地利用覆被變化及其生態(tài)環(huán)境效應的關系,以期為區(qū)域生態(tài)環(huán)境規(guī)劃提供參考。

    2 研究區(qū)概況

    淄博市位于魯中與魯北的交接地帶,臨淄區(qū)位于淄博市的東北部, 36°37′51″~37°00′30″N, 118°06′27″~118°29′30″E,地勢南高北低,并且向東北傾斜,最大海拔高度420 m,最低海拔高度31 m。臨淄區(qū)由南向北逐漸變緩,依次分布著丘陵、平原、平地、凹地等地貌單元,地質(zhì)特征是:古老的結晶基地上覆蓋有較新的沉積巖,是雙層結構,褶皺構造發(fā)育不完全,主要為斷裂構造。臨淄區(qū)屬于溫帶季風氣候,多年平均氣溫14.5 ℃,年極端最高氣溫39.4 ℃,出現(xiàn)在7月,年極端最低氣溫-16.6℃,出現(xiàn)在12月[4]。水資源主要以地下水為主,依靠降雨補充,多年平均降水量為650.2 mm,主要集中在夏季,全年最大日降水為50.3 mm。全區(qū)土地總面積663.68 km2,耕地面積約占70%。臨淄區(qū)境內(nèi)礦產(chǎn)資源非常豐富,主要有石灰石、鐵、鈷、黏土、煤以及各類稀有金屬等。臨淄區(qū)轄區(qū)內(nèi)有7個鎮(zhèn)、5個街道、414個行政村和66個社區(qū),總人口約61.74萬人,共有30個民族,以漢族為主。2016年,臨淄區(qū)的GDP總值為867.2億元,產(chǎn)業(yè)結構以第二產(chǎn)業(yè)為主。并且,臨淄區(qū)的糧食播種面積呈下降趨勢,工業(yè)及工業(yè)生產(chǎn)總值穩(wěn)步增長,產(chǎn)業(yè)結構在不斷調(diào)整中。

    2 數(shù)據(jù)與方法

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

    本研究中涉及的數(shù)據(jù)主要來自于3個部分:衛(wèi)星遙感數(shù)據(jù)、統(tǒng)計數(shù)據(jù)和實地調(diào)查數(shù)據(jù)。衛(wèi)星遙感數(shù)據(jù)主要使用了2008年、2012年和2018年共三期遙感影像,影像數(shù)據(jù)獲取自地理空間數(shù)據(jù)云(http://www.gscloud.cn)。統(tǒng)計數(shù)據(jù)主要來自《臨淄區(qū)統(tǒng)計年鑒》及相關學者公開發(fā)表的研究。實地調(diào)查數(shù)據(jù)來自2008年至2018年共5次對臨淄區(qū)范圍內(nèi)實地調(diào)查。

    2.2 土地利用現(xiàn)狀分類

    本研究參照《中國土地利用現(xiàn)狀調(diào)查技術規(guī)程》和中科院土地分類系統(tǒng)簡歷土地利用分類體系,并且從地理學與景觀生態(tài)學角度出發(fā),考慮到研究的自然地理狀況和經(jīng)濟發(fā)展狀況,在本研究中將土地利用類型的分類確定為耕地、林地、草地、水域、建設用地和其他用地六類。

    2.3 研究方法

    應用Erdas9.2進行波段融合和地理投影,ENVI5.3進行影響校正,用監(jiān)督分類方法進行土地利用分類,再用Arcgis分析數(shù)據(jù),獲取2008年、2012年和2018年臨淄區(qū)的土地利用/覆被類型數(shù)據(jù)并進行分析,借助Arcgis和Fragstats軟件進行生態(tài)環(huán)境效應各項指標的計算和繪圖。

    生態(tài)環(huán)境效應各項指標具體計算方法為:

    (1)生態(tài)環(huán)境質(zhì)量指數(shù):

    (1)

    式(1)中,LUi和Ci為研究區(qū)內(nèi)t時期第i類土地利用類型的面積及其生態(tài)環(huán)境質(zhì)量指數(shù);TA為區(qū)域總面積;n為區(qū)域內(nèi)土地利用類型數(shù)量[5]。

    (2)生態(tài)敏感性:

    (2)

    式(2)中,Ci表示空間單元對第i項因子的敏感指數(shù),m為第i個因子的指標個數(shù),Sk為第k個指標的敏感性空間分布。

    (3)生態(tài)貢獻率:

    LEI=(LEt+1-LEt)LA/TA

    (3)

    式(3)中,LEt+1、LEt分別為某種土地利用變化類型所反映的變化末期和初期土地利用變化類型所具有的生態(tài)環(huán)境質(zhì)量指數(shù);LA為該變化類型的面積;TA為區(qū)域總面積[6]。

    (4)生態(tài)風險指數(shù)。采用景觀指數(shù)法計算土地利用/覆被變化的生態(tài)風險[7],具體計算方法如下:

    Ri=Si×Fi

    (4)

    式(4)中,Ri為景觀格局指數(shù),Si為景觀干擾度指數(shù),F(xiàn)i為景觀脆弱度指數(shù)。

    Si=a×Ci+b×Ni+c×Di

    (5)

    式(5)中,Ci為景觀破碎度指數(shù),Ni為景觀分離度指數(shù),Di為景觀分維數(shù),a、b、c為對應景觀指數(shù)的權重,分別為0.5,0.3,0.2。

    (6)

    式(6)中,ERIi為風險小區(qū)的景觀風險指數(shù),Aki為第k個風險小區(qū)內(nèi)景觀類型i的面積,Ak為第k個小區(qū)的總面積,Ri為景觀類型i的景觀損失度指數(shù)。

    3 結果分析

    3.1 土地利用與覆被變化

    3.1.1 數(shù)量和結構變化

    2008年的土地利用類型為:草地面積5536827.179 m2,耕地面積478530122 m2,林地面積16072220.54 m2,水域面積2433781.045 m2,建設用地面積171883845.95 m2,其他用地8588815.77 m2,各個土地利用類型的面積大小排序依次為:耕地>建設用地>林地>其他用地>草地>水域,耕地的面積最大,水域的面積最小。

    2012年的土地利用類型為:草地10931915.19 m2,耕地457257708.6 m2,林地17605514.34 m2,水域2103624.918 m2,建設用地186893904.37 m2,其他用地7853634.31 m2,各個土地利用類型的面積大小排序依次為:耕地>建設用地>林地>草地>其他用地>水域,耕地的面積最大,水域的面積最小。

    2018年的土地利用類型為:草地13007616.93 m2,耕地439933796.4 m2,林地12937590.01 m2,水域1708786.1 m2,建設用地207540089.74 m2,其他用地7917890.44 m2,各個土地利用類型的面積大小排序依次為:耕地>建設用地>草地>林地>其他用地>水域,耕地的面積最大,水域的面積最小。

    在2008~2012年期間,臨淄區(qū)的土地利用的數(shù)量變化表現(xiàn)為:草地增加了5395088.008 m2,耕地減少了21272413.35 m2,林地增加了1533293.799 m2,水域減少了330156.1267 m2,建設用地增加了15010058.43 m2,其他用地減少了735181.47 m2。2008~2018年期間,草地增加了7470789.75 m2,耕地減少了-3134630.5 m2,林地減少了3134630.5 m2,水域減少了724994.94 m2,建設用地增加了35656243.79 m2,其他用地減少了670925.34 m2。

    綜上分析,臨淄區(qū)在2008年到2018年這期間的土地利用類型的變化規(guī)律明顯,具體表現(xiàn)為:草地面積持續(xù)增加,耕地面積持續(xù)減少,林地面積先增加后較少且減少量顯著大于增加量,水域面積持續(xù)減少,建設用地面積持續(xù)增加,其他用地面積先較少后增加且減少量大于增加量(圖1、表1)。

    2008~2012年,耕地向草地轉移了529.09 m2,向建設用地轉移了1479.30 m2,向林地轉移了164.13 m2;建設用地向耕地轉移了47.51 m2;林地向草地轉移了8.64 m2,向建設用地轉移了21.6 m2;其他用地向草地轉移了12.96 m2,向建設用地轉移了53.99 m2,向林地轉移了19.44 m2;水域向草地轉移了30.23 m2。

    在2012年,草地的增加面積有529.09 m2來自耕地,8.64 m2來自林地,12.96 m2來自其他用地,30.23 m2來自水域;耕地的增加面積有47.51 m2來自建設用地;建設用地的增加面積有1479.30 m2來自耕地,21.60 m2來自林地,53.99 m2來自其他用地;林地的變化面積有164.13 m2來自耕地,19.44 m2來自其他用地。

    2012~2018年,草地向耕地轉移了43.19 m2,向建設用地轉移了88.54 m2,向其他用地轉移了82.06 m2;耕地向草地轉移了51.83 m2,向建設用地轉移了2019.19 m2,向其他用地轉移了241.87 m2;建設用地有269.95 m2變成了耕地;林地有302.34 m2變成了草地,114.46 m2變成了耕地,51.83 m2變成了建設用地,10.80 m2變成了其他用地;其他用地中有177.08 m2變成了耕地,155.49 m2變成了建設用地;水域的面積中有36.71變成了建設用地。

    2018年,草地面積有51.83 m2來自于耕地,302.34 m2來自林地;耕地的面積中有43.19 m2來自于草地,269.95 m2來自于建設用地,114.46 m2來自于林地,177.08 m2來自于其他用地;建設用地的面積中有88.54 m2來自于草地,2019.19 m2來自于耕地,51.83 m2來自于林地,155.49 m2來自于其他用地,36.71 m2來自于水域;其他用地的面積中,82.06 m2來自于草地,241.87 m2來自于耕地,10.80 m2來自于林地(表1)。

    圖1 2008~2018年臨淄區(qū)土地利用現(xiàn)狀

    表1 2008~2018年臨淄區(qū)土地利用變化轉移矩陣 m2

    3.1.2 重心變化

    在2008年,臨淄區(qū)的草地比例為0.81%,耕地比重為70.06%,建設用地比重為23.12%,林地比重為2.35%,其他用地比重為3.30%,水域比重為0.36%。此時的土地利用類型以耕地為主,建設用地次之,草地、林地、其他用地和水域都極少。

    在2012年,草地的比重為1.60%,比2008年增加了0.79%;耕地的比重為66.98%,比2008年減少了3.08%;建設用地的比重為25.35%,比2008年增加了2.23%;林地的比重為2.58%,比2008年增加了0.53%;其他用地的比重為3.18%,比2008年減少了0.12%;水域的比重為0.31%,比2008年減少了0.05%??梢钥闯觯瑥?008年到2012年,臨淄區(qū)土地利用重心由耕地緩慢移向建設用地。

    在2018年,草地的比重為1.90%,比2012年增加了0.3%,比2008年增加了1.09%;耕地的比重為64.61%,比2012年下降了2.58%,相較于2008年下降了5.65%;建設用地的比重為28.36%,比2012年增加了3.01%,比2008年增加了5.24%;林地的比重為1.89%,比2012年下降了0.68%,比2008年降低了0.46%;其他用地的比重為3.18%,和2012年相同,比2008年下降了0.12%;水域的比重為0.25%,比2012年降低了0.06%,比2008年降低了0.11%。

    總體來看,從2008年到2018年,這10年間,臨淄區(qū)的土地利用類型中,耕地、林地、其他用地和水域的比重下降,草地和建設用地的比重增加,土地利用和覆被的重心由耕地、林地、水域及其他用地類型逐漸轉向建設用地和草地類型(表2)。

    3.2 生態(tài)環(huán)境效應分析

    3.2.1 生態(tài)環(huán)境質(zhì)量

    臨淄區(qū)2008年的生態(tài)環(huán)境質(zhì)量指數(shù)為0.24713118,2012年的生態(tài)環(huán)境質(zhì)量指數(shù)為0.242766019,2018年的生態(tài)環(huán)境質(zhì)量指數(shù)為0.23242998。2008年到2012年生態(tài)環(huán)境質(zhì)量指數(shù)的年變化率為-0.00109129,2012年到2018年的生態(tài)環(huán)境質(zhì)量指數(shù)的年變化率為—0.001722672。

    表2 臨淄區(qū)土地利用重心變化 %

    2008年到2018年期間的土地利用變化導致臨淄區(qū)的生態(tài)環(huán)境質(zhì)量逐年下降,可能是由于耕地逐年減少、建設用地逐年增加,土地破碎化逐年程度增加而造成的(圖2)。

    圖2 臨淄區(qū)2008~2018年生態(tài)環(huán)境質(zhì)量指數(shù)變化情況

    3.2.2 生態(tài)敏感

    臨淄區(qū)草地的生態(tài)敏感指數(shù)在2008年為0.04821,敏感程度為低敏感;在2012年敏感指數(shù)降低為0.4609,敏感程度為低敏感;在2018年,敏感性指數(shù)為0.04438,敏感程度依舊為低敏感。建設用地的敏感性指數(shù)在2008年為0.01750,2012年為0.01899,2018年為0.02108,敏感在3個時間點都是不敏感。水域的生態(tài)敏感性指數(shù)在2008年為0.00024,2012年為0.00021,2018年為0.00017,敏感都是高敏感,林地的生態(tài)敏感性指數(shù)在2008年為0.00057,2012年為0.00116,2018年為0.00080,敏感程度都是中敏感。其他用地的生態(tài)敏感程度在2008年為0.00088,2012年為0.00081,2018年0.00130,敏感程度都是不敏感。耕地的生態(tài)敏感性指數(shù)在2008年為0.00164,2012年為0.00180,2018年為0.00132,敏感程度都是高敏感。

    總體分析來看,雖然臨淄區(qū)的6個土地利用類型的生態(tài)敏感性指數(shù)在2008年到2018年有不同程度的變化,但是生態(tài)敏感沒有變化。其中,耕地和水域時高敏感區(qū)域,林地為中敏感,草地低敏感,建設用地和其他用地都是不敏感(圖3)。

    圖3 臨淄區(qū)2008年~2018年不同土地利用類型生態(tài)敏感分布變化

    3.2.3 生態(tài)貢獻率

    2008~2012年,造成臨淄區(qū)的生態(tài)環(huán)境質(zhì)量惡化的土地利用變化驅動因子有4個:耕地轉換為建設用地、建設用地轉換為其他用地、水域轉換為草地,林地轉換為草地,其中,耕地轉換為建設用地的生態(tài)貢獻率最大,占總貢獻率的98.5392%。造成臨淄區(qū)的生態(tài)環(huán)境質(zhì)量改善的驅動因子有7個,分別是:耕地轉換為林地、耕地轉換為草地、建設用地轉換為草地、建設用地轉換為耕地、其他用地轉換為建設用地、其他用地轉換為林地、林地轉換為建設用地,其中,耕地轉換為林地和耕地轉換為草地兩種土地利用變化的生態(tài)貢獻率最大,分別為25.9869%和48.7411%??傮w來看,土地利用變化的生態(tài)環(huán)境惡化貢獻率為0.000000637,改善的生態(tài)環(huán)境貢獻率為0.0000002544 ,惡化的貢獻率大于了改善的貢獻率,所以,臨淄區(qū)2008年到2012年的土地利用變化造成臨淄區(qū)生態(tài)環(huán)境惡化。

    2012年到2018年,造成臨淄區(qū)的生態(tài)環(huán)境質(zhì)量惡化的土地利用變化的驅動因子有7個,分別為:草地轉換為耕地,占總生態(tài)貢獻率的3.4008%;草地轉換為建設用地,占總生態(tài)貢獻率的19.6078;草地轉換為其他用地,占總生態(tài)貢獻率的18.7789%;林地轉換為草地,占總生態(tài)貢獻率的17.1104%;林地轉換為耕地,占總生態(tài)貢獻率的15.4897%;林地轉換為建設用地,占總生態(tài)貢獻率的14.4110%;水域轉換為建設用地,占總生態(tài)貢獻率的11.2014%。造成臨淄區(qū)生態(tài)環(huán)境改善的土地利用變化的驅動因子有4個,分別是:耕地轉換為草地,占總生態(tài)貢獻率的5.7932%;建設用地轉換為耕地,占總生態(tài)貢獻率的54.6881%;其他用地轉換為耕地,占總生態(tài)貢獻率的37.7310%;其他用地轉換為建設用地,占總生態(tài)貢獻率的1.7877%??傮w來看,土地利用變化的生態(tài)環(huán)境惡化的貢獻率為0.0000002977,生態(tài)環(huán)境改善的貢獻率為0.0000002097,惡化的貢獻率大于了改善的貢獻率,故而2012年到2018年臨淄區(qū)的土地利用變化造成了臨淄區(qū)生態(tài)環(huán)境質(zhì)量的惡化。

    綜上,對臨淄區(qū)生態(tài)環(huán)境質(zhì)量影響的土地利用類型變化的驅動因子中,2008~2012年,耕地轉換為建設用地是生態(tài)環(huán)境質(zhì)量惡化的主要因子,耕地轉換為草地是生態(tài)環(huán)境質(zhì)量改善的主要因子;2012~2018年,草地轉換為建設用地是生態(tài)環(huán)境質(zhì)量惡化的主要因子,建設用地轉換為耕地是生態(tài)環(huán)境質(zhì)量改善的主要因子(表3)。

    表3 2008~2018 年臨淄區(qū)土地利用變化的生態(tài)貢獻率

    3.2.4 生態(tài)風險

    2008年,臨淄區(qū)低風險的面積為16072221m2,較低風險的面積為5536827 m2,中風險的面積為480960485 m2,較高風險的面積為131224049 m2,高風險的面積為49254203 m2。各個風險等級的面積大小排序依次為:中風險>較高風險>高風險>低風險>較低風險。

    2012年,臨淄區(qū)低風險的面積為17413670 m2,較低風險的面積為11331221 m2,中風險的面積為458178866 m2,較高風險的面積為124443930 m2,高風險的面積為71680011 m2,高風險的面積為75322666 m2。各個風險等級的面積大小排序依次為:中風險>較高風險>高風險>低風險>較低風險。

    2018年,臨淄區(qū)低風險的面積為12937590 m2,較低風險的面積為13007617 m2,中風險的面積為441639227 m2,較高風險的面積為140141654 m2。各個風險等級的面積大小排序依次為:中風險>較高風險>高風險>較低風險>低風險。

    綜合來看,從2008~2018年,中生態(tài)風險的面積不斷下降,較高生態(tài)風險的面積逐漸增加,高生態(tài)風險的面積逐漸增加,低風險的面積逐漸降低,較低風險的面積逐漸增加。2008~2018年期間的土地利用類型變化導致臨淄區(qū)中高風險、較高風險和低風險區(qū)域的面積不斷增加,而中風險和低風險的面積不斷減少,臨淄區(qū)的區(qū)域生態(tài)風險水平增加(表4)。

    表4 2008~2018年臨淄區(qū)各風險等級面積 m2

    4 討論與結論

    臨淄區(qū)2008~2018年的主要土地利用類型為耕地和建設用地,草地和建設用地的面積增加,耕地、林地、水域和其他用地的面積減小。并且,耕地減少的面積最大,建設用地增加的面積最多。2008年到2012年,草地和建設用地增加的面積主要來自于耕地。2012年到2018年,草地增加的面積來自于耕地和林地,來自林地的部分更多;建設用地增加的面積主要來自于耕地。草地面積的增加可能與“退耕還草”的生態(tài)環(huán)境治理政策有關,建設用地的增加則是社會經(jīng)濟飛速發(fā)展和人口迅速增加的結果。

    臨淄區(qū)2008年到2018年的土地利用類型中,耕地和水域為生態(tài)高敏感類型,林地為中敏感,草地低敏感,建設用地和其他用地不敏感。在這十年內(nèi),生態(tài)環(huán)境質(zhì)量指數(shù)逐年下降,土地利用的變化對生態(tài)環(huán)境質(zhì)量造成了消極影響。對臨淄區(qū)生態(tài)環(huán)境質(zhì)量影響的土地利用類型變化的驅動因子中,2008~2012年,耕地轉換為建設用地是生態(tài)環(huán)境質(zhì)量惡化的主要因子,耕地轉換為草地是生態(tài)環(huán)境質(zhì)量改善的主要因子;2012~2018年,草地轉換為建設用地是生態(tài)環(huán)境質(zhì)量惡化的主要因子,建設用地轉換為耕地是生態(tài)環(huán)境質(zhì)量改善的主要因子。以上的規(guī)律說明,建設用地的增加是導致生態(tài)環(huán)境質(zhì)量下降的主要原因,而草地、耕地面積的增加則有助于改善區(qū)域的生態(tài)環(huán)境質(zhì)量。

    2008~2018年期間的土地利用類型變化導致臨淄區(qū)中高風險、較高風險和低風險區(qū)域的面積不斷增加,而中風險和低風險的面積不斷減少。說明臨淄區(qū)的整體區(qū)域生態(tài)風險在不斷增加,這可能與建設用地增加、耕地林地減少,生態(tài)環(huán)境質(zhì)量下降,生態(tài)功能穩(wěn)定性降低有關。

    猜你喜歡
    臨淄區(qū)貢獻率林地
    Summer Is Coming
    平衡車
    明水縣林地資源現(xiàn)狀及動態(tài)變化淺析
    眼睛
    淺談明水縣林地資源動態(tài)變化狀況
    14.8%
    高等教育體制改革的貢獻率研究:1998—2011年的數(shù)據(jù)觀察
    我和板栗捉迷藏
    淺談林地的特性及其重要地位
    午夜福利在线观看吧| 757午夜福利合集在线观看| 不卡一级毛片| 午夜福利,免费看| 久久精品国产亚洲av高清一级| 三级毛片av免费| 男女高潮啪啪啪动态图| 亚洲色图av天堂| 丁香六月天网| 国产免费福利视频在线观看| 亚洲成国产人片在线观看| 99久久99久久久精品蜜桃| 12—13女人毛片做爰片一| 18禁国产床啪视频网站| 亚洲精华国产精华精| 自拍欧美九色日韩亚洲蝌蚪91| 99精品在免费线老司机午夜| 亚洲全国av大片| 十八禁网站网址无遮挡| 别揉我奶头~嗯~啊~动态视频| 成人18禁高潮啪啪吃奶动态图| 成年人午夜在线观看视频| 丝袜喷水一区| 性高湖久久久久久久久免费观看| 在线观看免费午夜福利视频| 窝窝影院91人妻| 变态另类成人亚洲欧美熟女 | 后天国语完整版免费观看| 夫妻午夜视频| 亚洲欧美激情在线| 高清视频免费观看一区二区| 久久av网站| 久久国产精品大桥未久av| 亚洲色图 男人天堂 中文字幕| 成人影院久久| 亚洲中文字幕日韩| 亚洲精品美女久久久久99蜜臀| 99riav亚洲国产免费| 黑人猛操日本美女一级片| a级毛片黄视频| 欧美在线黄色| 日韩大码丰满熟妇| 看免费av毛片| 成人国产一区最新在线观看| 欧美日韩亚洲综合一区二区三区_| 一区二区av电影网| 女人爽到高潮嗷嗷叫在线视频| 男女午夜视频在线观看| 91av网站免费观看| 国产精品久久电影中文字幕 | 亚洲成人免费电影在线观看| 老汉色av国产亚洲站长工具| 国产精品秋霞免费鲁丝片| 啦啦啦视频在线资源免费观看| 黄色 视频免费看| 十八禁网站网址无遮挡| 亚洲人成伊人成综合网2020| av不卡在线播放| 丁香六月天网| 欧美精品一区二区大全| 国产精品一区二区在线不卡| 亚洲国产欧美日韩在线播放| 成在线人永久免费视频| 十分钟在线观看高清视频www| 精品久久久精品久久久| 天天影视国产精品| 亚洲国产欧美一区二区综合| av一本久久久久| 国产精品99久久99久久久不卡| 国产主播在线观看一区二区| 欧美精品一区二区免费开放| 精品国内亚洲2022精品成人 | 午夜福利视频在线观看免费| 欧美在线一区亚洲| www.熟女人妻精品国产| 亚洲精品成人av观看孕妇| 黄片小视频在线播放| 免费人妻精品一区二区三区视频| 久久 成人 亚洲| 大型黄色视频在线免费观看| 国产有黄有色有爽视频| 国产亚洲午夜精品一区二区久久| 婷婷成人精品国产| 啦啦啦 在线观看视频| 在线观看免费午夜福利视频| 夜夜骑夜夜射夜夜干| 国产在线观看jvid| 99riav亚洲国产免费| 丰满迷人的少妇在线观看| 精品一区二区三区视频在线观看免费 | 淫妇啪啪啪对白视频| 亚洲精品国产区一区二| 免费高清在线观看日韩| 国产精品免费大片| 女同久久另类99精品国产91| 国产一区二区激情短视频| 日韩一卡2卡3卡4卡2021年| 在线观看舔阴道视频| 亚洲成人免费av在线播放| 欧美激情高清一区二区三区| 高清在线国产一区| 女人被躁到高潮嗷嗷叫费观| 日韩欧美免费精品| 搡老乐熟女国产| 精品国产超薄肉色丝袜足j| 91精品三级在线观看| 欧美日韩福利视频一区二区| 免费在线观看视频国产中文字幕亚洲| 国产免费av片在线观看野外av| 汤姆久久久久久久影院中文字幕| 久久 成人 亚洲| 精品亚洲乱码少妇综合久久| 久久久久久久久久久久大奶| 久久青草综合色| 大型黄色视频在线免费观看| 啪啪无遮挡十八禁网站| 国产99久久九九免费精品| 男人舔女人的私密视频| 久久久久久久精品吃奶| 黄色丝袜av网址大全| 在线观看www视频免费| 欧美日本中文国产一区发布| 国产区一区二久久| 在线观看免费视频日本深夜| 国产精品九九99| 最新在线观看一区二区三区| 亚洲精品成人av观看孕妇| 十八禁网站网址无遮挡| 婷婷丁香在线五月| 不卡一级毛片| 亚洲七黄色美女视频| svipshipincom国产片| 水蜜桃什么品种好| 搡老岳熟女国产| 91成年电影在线观看| 18禁黄网站禁片午夜丰满| 制服人妻中文乱码| 午夜激情av网站| 久久久精品免费免费高清| 久久99一区二区三区| 99热网站在线观看| 91精品国产国语对白视频| 建设人人有责人人尽责人人享有的| av欧美777| 亚洲中文av在线| 亚洲精品自拍成人| 人妻久久中文字幕网| 国产精品久久久久成人av| 亚洲全国av大片| 人人妻人人爽人人添夜夜欢视频| av网站免费在线观看视频| 国产成人欧美在线观看 | 久久精品亚洲精品国产色婷小说| 亚洲国产欧美日韩在线播放| 国产又色又爽无遮挡免费看| 国产成人精品在线电影| 别揉我奶头~嗯~啊~动态视频| 久久99热这里只频精品6学生| 老司机在亚洲福利影院| 欧美+亚洲+日韩+国产| 日韩三级视频一区二区三区| 男人操女人黄网站| 搡老乐熟女国产| 久热这里只有精品99| 精品国产一区二区三区久久久樱花| 在线观看免费午夜福利视频| 黄色a级毛片大全视频| 亚洲色图av天堂| 久久狼人影院| 国产成人精品久久二区二区91| 中文字幕色久视频| 在线观看免费高清a一片| 精品国产亚洲在线| 一区福利在线观看| 男女之事视频高清在线观看| 操美女的视频在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲精品中文字幕在线视频| 国产一卡二卡三卡精品| a级毛片在线看网站| 人妻久久中文字幕网| av不卡在线播放| 一个人免费看片子| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 少妇猛男粗大的猛烈进出视频| 天天躁夜夜躁狠狠躁躁| 五月开心婷婷网| 99久久国产精品久久久| 午夜福利视频精品| 国产一卡二卡三卡精品| 欧美黄色淫秽网站| 欧美成人午夜精品| videosex国产| 欧美国产精品一级二级三级| 国产亚洲精品第一综合不卡| tocl精华| 777久久人妻少妇嫩草av网站| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看 | 国产av一区二区精品久久| 在线观看免费午夜福利视频| 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 法律面前人人平等表现在哪些方面| 欧美日韩亚洲高清精品| 国产福利在线免费观看视频| 色综合婷婷激情| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| 中文亚洲av片在线观看爽 | 亚洲免费av在线视频| 麻豆成人av在线观看| 午夜福利在线免费观看网站| 黄色a级毛片大全视频| 国产成人系列免费观看| 99久久国产精品久久久| 超碰97精品在线观看| 精品免费久久久久久久清纯 | 飞空精品影院首页| 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 久久中文看片网| 首页视频小说图片口味搜索| 日韩免费av在线播放| 男女下面插进去视频免费观看| 国产av国产精品国产| 久久精品aⅴ一区二区三区四区| 在线观看免费午夜福利视频| 激情视频va一区二区三区| 香蕉丝袜av| 黄片播放在线免费| av欧美777| 国产成人一区二区三区免费视频网站| 成人特级黄色片久久久久久久 | 亚洲精品av麻豆狂野| 精品视频人人做人人爽| 桃花免费在线播放| 在线亚洲精品国产二区图片欧美| 亚洲av国产av综合av卡| 国产在视频线精品| 亚洲精品美女久久av网站| 精品午夜福利视频在线观看一区 | 1024视频免费在线观看| 免费日韩欧美在线观看| 精品久久久久久久毛片微露脸| 最新美女视频免费是黄的| 精品亚洲成a人片在线观看| 在线观看免费日韩欧美大片| 国产精品98久久久久久宅男小说| 亚洲精品中文字幕一二三四区 | 大型av网站在线播放| 窝窝影院91人妻| 国产精品国产av在线观看| 亚洲国产av影院在线观看| 日本黄色视频三级网站网址 | 美女扒开内裤让男人捅视频| 国产一区有黄有色的免费视频| 99香蕉大伊视频| 国产精品一区二区精品视频观看| 久久中文字幕人妻熟女| 91精品国产国语对白视频| 国产不卡av网站在线观看| 在线av久久热| av线在线观看网站| 精品熟女少妇八av免费久了| 国产精品久久久久久精品古装| cao死你这个sao货| 99九九在线精品视频| 99精品欧美一区二区三区四区| 国产亚洲欧美精品永久| 在线观看人妻少妇| 波多野结衣一区麻豆| 国产精品久久久久久精品电影小说| 91成年电影在线观看| 91老司机精品| 久久精品国产亚洲av香蕉五月 | 免费日韩欧美在线观看| 中文字幕人妻丝袜制服| 在线观看舔阴道视频| 久久久久精品人妻al黑| 丰满迷人的少妇在线观看| av在线播放免费不卡| 操美女的视频在线观看| 韩国精品一区二区三区| 精品熟女少妇八av免费久了| 丝袜在线中文字幕| 中文欧美无线码| 9热在线视频观看99| 精品国产一区二区三区久久久樱花| 亚洲精品久久午夜乱码| 宅男免费午夜| 正在播放国产对白刺激| 久热这里只有精品99| 久久久久视频综合| 久久久国产欧美日韩av| 久久久久久亚洲精品国产蜜桃av| av天堂在线播放| 久久久国产成人免费| 国产日韩欧美视频二区| 18在线观看网站| 9热在线视频观看99| 精品国产乱码久久久久久男人| 国精品久久久久久国模美| 午夜福利,免费看| 欧美国产精品一级二级三级| 久热这里只有精品99| 日韩中文字幕视频在线看片| 日韩成人在线观看一区二区三区| 十分钟在线观看高清视频www| 大型黄色视频在线免费观看| 国产国语露脸激情在线看| 999久久久精品免费观看国产| 电影成人av| 国产精品久久久久成人av| 免费高清在线观看日韩| 午夜91福利影院| 国产午夜精品久久久久久| 人人澡人人妻人| 国产精品国产av在线观看| 亚洲人成电影免费在线| 纯流量卡能插随身wifi吗| 在线观看免费日韩欧美大片| 在线观看www视频免费| 中国美女看黄片| 老司机在亚洲福利影院| 一级片免费观看大全| 欧美 亚洲 国产 日韩一| 岛国毛片在线播放| 一级毛片精品| 亚洲成av片中文字幕在线观看| 国产高清国产精品国产三级| kizo精华| 久久热在线av| 日韩欧美一区二区三区在线观看 | 久久影院123| 考比视频在线观看| 法律面前人人平等表现在哪些方面| 十八禁高潮呻吟视频| 日本一区二区免费在线视频| 在线亚洲精品国产二区图片欧美| 亚洲熟女精品中文字幕| 亚洲国产成人一精品久久久| 人人妻人人澡人人看| 成人av一区二区三区在线看| 国产精品美女特级片免费视频播放器 | 亚洲av第一区精品v没综合| 精品久久久精品久久久| 肉色欧美久久久久久久蜜桃| 热99国产精品久久久久久7| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美网| 国产免费现黄频在线看| 精品高清国产在线一区| 色94色欧美一区二区| 啦啦啦视频在线资源免费观看| 欧美日本中文国产一区发布| 99国产精品99久久久久| 99精国产麻豆久久婷婷| 亚洲自偷自拍图片 自拍| 久久中文字幕一级| 无限看片的www在线观看| 精品乱码久久久久久99久播| 成人黄色视频免费在线看| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月 | 国产黄色免费在线视频| 国产精品1区2区在线观看. | 亚洲九九香蕉| 精品久久久精品久久久| 亚洲av日韩在线播放| www日本在线高清视频| 免费女性裸体啪啪无遮挡网站| 一区二区三区乱码不卡18| 国产精品av久久久久免费| 99久久99久久久精品蜜桃| 亚洲人成77777在线视频| 国产成人精品久久二区二区免费| 99精国产麻豆久久婷婷| 老司机亚洲免费影院| 18禁美女被吸乳视频| 99国产精品一区二区三区| 亚洲精品中文字幕一二三四区 | 国产av又大| 精品一区二区三区视频在线观看免费 | 自线自在国产av| 午夜老司机福利片| 丝瓜视频免费看黄片| 精品久久久久久电影网| 日韩成人在线观看一区二区三区| 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 人人澡人人妻人| 日本黄色日本黄色录像| 男男h啪啪无遮挡| 亚洲av日韩精品久久久久久密| 色播在线永久视频| 男女之事视频高清在线观看| 欧美精品亚洲一区二区| 桃红色精品国产亚洲av| 悠悠久久av| 欧美日韩福利视频一区二区| 亚洲少妇的诱惑av| 一级片'在线观看视频| 久热这里只有精品99| 搡老乐熟女国产| 深夜精品福利| 在线观看免费高清a一片| 免费高清在线观看日韩| 岛国在线观看网站| 汤姆久久久久久久影院中文字幕| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成伊人成综合网2020| 国产xxxxx性猛交| 亚洲精品国产一区二区精华液| 国产主播在线观看一区二区| 91老司机精品| 一个人免费在线观看的高清视频| 搡老熟女国产l中国老女人| 18禁美女被吸乳视频| 欧美日韩av久久| 久久免费观看电影| 超碰97精品在线观看| 日韩有码中文字幕| 无遮挡黄片免费观看| 2018国产大陆天天弄谢| 亚洲,欧美精品.| 精品少妇久久久久久888优播| 国产精品.久久久| 欧美日韩国产mv在线观看视频| 免费观看av网站的网址| 黄色成人免费大全| 亚洲色图 男人天堂 中文字幕| 亚洲欧洲日产国产| 国产精品亚洲av一区麻豆| 国产成人欧美| 久久精品亚洲av国产电影网| 变态另类成人亚洲欧美熟女 | 午夜久久久在线观看| 亚洲美女黄片视频| 两人在一起打扑克的视频| 日韩制服丝袜自拍偷拍| 久久99一区二区三区| 2018国产大陆天天弄谢| 国产日韩欧美亚洲二区| 欧美精品一区二区大全| 男女高潮啪啪啪动态图| 美女午夜性视频免费| 国产成人免费观看mmmm| 午夜精品国产一区二区电影| 超色免费av| 久久国产精品影院| 大型av网站在线播放| 超碰成人久久| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 深夜精品福利| 久久午夜综合久久蜜桃| 人妻 亚洲 视频| 大片电影免费在线观看免费| 免费看a级黄色片| 国产成人一区二区三区免费视频网站| 国产成人欧美| 波多野结衣一区麻豆| 久久久精品区二区三区| a在线观看视频网站| 亚洲欧美日韩另类电影网站| 男女下面插进去视频免费观看| 不卡av一区二区三区| 日韩欧美免费精品| 一本久久精品| 日韩视频一区二区在线观看| 9191精品国产免费久久| 日韩熟女老妇一区二区性免费视频| 看免费av毛片| 久久中文看片网| 色尼玛亚洲综合影院| 国产不卡一卡二| 一本—道久久a久久精品蜜桃钙片| 日本五十路高清| 无人区码免费观看不卡 | 欧美人与性动交α欧美精品济南到| 日本黄色视频三级网站网址 | 欧美日韩一级在线毛片| 丝袜在线中文字幕| 国产免费视频播放在线视频| 国产在线免费精品| 脱女人内裤的视频| 精品少妇黑人巨大在线播放| 国产一区二区 视频在线| 蜜桃在线观看..| 丰满人妻熟妇乱又伦精品不卡| 国产淫语在线视频| 丰满人妻熟妇乱又伦精品不卡| 国产欧美亚洲国产| 在线 av 中文字幕| 成人黄色视频免费在线看| 91成人精品电影| 国产高清videossex| 女人久久www免费人成看片| 免费日韩欧美在线观看| tocl精华| 在线观看一区二区三区激情| 黑人欧美特级aaaaaa片| 黑丝袜美女国产一区| 欧美日韩av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品高潮呻吟av久久| 搡老乐熟女国产| 人人妻人人澡人人爽人人夜夜| 黄色视频在线播放观看不卡| 操出白浆在线播放| 精品一品国产午夜福利视频| 丁香六月欧美| 老鸭窝网址在线观看| 丝袜在线中文字幕| 色综合欧美亚洲国产小说| 亚洲成av片中文字幕在线观看| 日韩免费高清中文字幕av| 黄片大片在线免费观看| 色在线成人网| 91老司机精品| 51午夜福利影视在线观看| 妹子高潮喷水视频| 成年女人毛片免费观看观看9 | 天天添夜夜摸| 欧美成人午夜精品| 国产区一区二久久| 人人妻人人爽人人添夜夜欢视频| 高清欧美精品videossex| 久久性视频一级片| 国产精品麻豆人妻色哟哟久久| 久久影院123| 国产精品av久久久久免费| 麻豆乱淫一区二区| 久久久久久亚洲精品国产蜜桃av| 久久av网站| 一本大道久久a久久精品| 一区二区三区精品91| 国产男靠女视频免费网站| 制服人妻中文乱码| 91av网站免费观看| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 十八禁网站网址无遮挡| 国产精品国产高清国产av | 国产成人精品久久二区二区免费| 视频在线观看一区二区三区| 精品人妻1区二区| 中文字幕最新亚洲高清| 香蕉丝袜av| 欧美日韩黄片免| 狂野欧美激情性xxxx| 99re6热这里在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 久久九九热精品免费| 国产成人系列免费观看| 国产99久久九九免费精品| 人成视频在线观看免费观看| 欧美+亚洲+日韩+国产| 国产片内射在线| 成人免费观看视频高清| 国产亚洲av高清不卡| 日韩大片免费观看网站| 天堂动漫精品| 手机成人av网站| 黄色片一级片一级黄色片| 少妇被粗大的猛进出69影院| 亚洲午夜理论影院| 成人手机av| 99re在线观看精品视频| 少妇猛男粗大的猛烈进出视频| 麻豆av在线久日| 国产亚洲精品第一综合不卡| 我的亚洲天堂| 啪啪无遮挡十八禁网站| 99国产极品粉嫩在线观看| 国产在线精品亚洲第一网站| 久久99一区二区三区| 久久人人97超碰香蕉20202| 中文字幕最新亚洲高清| 男女之事视频高清在线观看| 中文字幕制服av| 国产片内射在线| 国产亚洲精品久久久久5区| 成人av一区二区三区在线看| 成人18禁在线播放| 日韩欧美免费精品| 午夜视频精品福利| 国产精品 欧美亚洲| 黄色视频,在线免费观看| 超碰97精品在线观看| 国产男女内射视频| 成在线人永久免费视频| 婷婷成人精品国产| 国产精品欧美亚洲77777| 亚洲九九香蕉| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片女人18水好多| 久久精品91无色码中文字幕| 丝袜喷水一区| 久久 成人 亚洲| 亚洲人成电影观看| 国产真人三级小视频在线观看| 精品少妇内射三级| 欧美精品一区二区免费开放| av欧美777| 国产亚洲精品一区二区www | 一级毛片精品| 丝袜人妻中文字幕| 久久国产精品男人的天堂亚洲| 亚洲成人手机| 制服诱惑二区| 久久毛片免费看一区二区三区| 超碰成人久久|