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

    全國(guó)縣域復(fù)種指數(shù)時(shí)空變化特征研究

    2021-12-08 09:29:23淮賀舉李奇峰王小慧胡海棠李存軍
    江蘇農(nóng)業(yè)科學(xué) 2021年22期
    關(guān)鍵詞:重心時(shí)空變化縣域

    淮賀舉 李奇峰 王小慧 胡海棠 李存軍

    摘要:復(fù)種指數(shù)是衡量耕地資源集約化利用程度的基礎(chǔ)性指標(biāo),提高復(fù)種指數(shù)是糧食持續(xù)增產(chǎn)的重要途徑之一。為分析全國(guó)復(fù)種指數(shù)變化特征、空間差異、重心遷移,系統(tǒng)收集耕地面積、農(nóng)作物總播種面積等數(shù)據(jù),基于GIS 空間展示與分析,應(yīng)用重心分析法、全局空間自相關(guān)、局部空間自相關(guān)等方法,研究全國(guó)復(fù)種指數(shù)時(shí)空變化特征。結(jié)果顯示:(1)1998年,全國(guó)復(fù)種指數(shù)最高,達(dá)到165%,目前保持在123%左右。(2)從空間上看,復(fù)種指數(shù)呈現(xiàn)從北到南逐漸增加的趨勢(shì);從時(shí)間上看,復(fù)種指數(shù)大多數(shù)地區(qū)呈現(xiàn)“上升—降低—上升”的趨勢(shì)。1995年,復(fù)種指數(shù)在201%~300%之間的區(qū)域范圍最大,之后2005年和2015年不斷減少,同時(shí)有向西北方向移動(dòng)的趨勢(shì)。1985—2015年,全國(guó)復(fù)種指數(shù)上升區(qū)域主要集中在東北、黃淮海、西南、四川盆地農(nóng)作區(qū),下降區(qū)域主要集中在北部低中高原、長(zhǎng)江中下游(東南部)、江南丘陵、華南農(nóng)作區(qū)。(3)1985—2005年空間上的聚集態(tài)勢(shì)增加,2005—2015年降低。與1985年相比,2015年全國(guó)復(fù)種指數(shù)重心向西南移動(dòng)33 km,其中復(fù)種指數(shù)<100%的重心向東南移動(dòng)433 km,復(fù)種指數(shù)100%~200%的重心向東北移動(dòng)104 km,復(fù)種指數(shù)>200%的重心向東北移動(dòng) 319 km。我國(guó)復(fù)種指數(shù)存在空間上的自相關(guān)性,南方地區(qū)復(fù)種指數(shù)降低、重心向西南遷移,是我國(guó)復(fù)種指數(shù)變化的最重要特征。

    關(guān)鍵詞:復(fù)種指數(shù);時(shí)空變化;縣域;空間相關(guān)性;重心

    中圖分類號(hào): S126;F326.11? 文獻(xiàn)標(biāo)志碼: A

    文章編號(hào):1002-1302(2021)22-0076-08

    收稿日期:2020-12-31

    基金項(xiàng)目:國(guó)家重點(diǎn)研發(fā)計(jì)劃(編號(hào):2016YFD0300201);北京市農(nóng)林科學(xué)院青年科研基金(編號(hào):QNJJ202126)。

    作者簡(jiǎn)介:淮賀舉(1983—),男,河北邢臺(tái)人,博士,助理研究員,主要從事農(nóng)業(yè)區(qū)域發(fā)展與農(nóng)業(yè)信息化研究。E-mail:huaihj@nercita.org.cn。

    通信作者:李存軍,博士,研究員,主要從事農(nóng)業(yè)數(shù)字生態(tài)研究。E-mail:licj@nercita.org.com。

    糧食安全是關(guān)系國(guó)民經(jīng)濟(jì)發(fā)展、社會(huì)穩(wěn)定和國(guó)家自立的全局性重大戰(zhàn)略問題。在國(guó)內(nèi)外形勢(shì)復(fù)雜變化的當(dāng)下,保障糧食安全的戰(zhàn)略地位十分重要。耕地面積、復(fù)種指數(shù)、單位面積產(chǎn)量是決定我國(guó)糧食產(chǎn)量的三大關(guān)鍵因素。隨著城市化、工業(yè)化進(jìn)程的加快,我國(guó)大量耕地被占用,非農(nóng)化、非糧化問題日益嚴(yán)重。如何在保證耕地面積不減少、糧食單產(chǎn)逐步提高的基礎(chǔ)上,提高復(fù)種指數(shù),擴(kuò)大農(nóng)作物播種面積,是我國(guó)未來實(shí)現(xiàn)糧食持續(xù)增產(chǎn)的重要途徑之一。

    多熟種植是我國(guó)重要的農(nóng)作制度之一,對(duì)于提高耕地生產(chǎn)能力和確保糧食安全具有重要作用[1]。復(fù)種指數(shù)是耕作制度研究中衡量耕地資源集約化利用程度的基礎(chǔ)性指標(biāo),也是宏觀評(píng)價(jià)耕地資源利用基本狀況的重要技術(shù)指標(biāo)。復(fù)種指數(shù)分為潛在復(fù)種指數(shù)和實(shí)際復(fù)種指數(shù)[2],潛在復(fù)種指數(shù)常利用農(nóng)業(yè)氣象法估算充分利用自然資源時(shí)能達(dá)到的最大復(fù)種指數(shù),實(shí)際復(fù)種指數(shù)一般基于不同行政單元播種面積和耕地面積的統(tǒng)計(jì)數(shù)據(jù)計(jì)算。黃國(guó)勤分析了南方復(fù)種指數(shù)的演變與特點(diǎn),并通過回歸分析法研究了耕地復(fù)種指數(shù)與糧食總產(chǎn)的關(guān)系,提出了提高復(fù)種指數(shù)的途徑與措施[3]。劉巽浩對(duì)不同區(qū)域的熟制和水熱資源進(jìn)行了實(shí)際調(diào)查,估算了全國(guó)耕地復(fù)種指數(shù)[4]。段紅平分析了1949—1998年湖南省復(fù)種指數(shù)的歷史變化,認(rèn)為提高復(fù)種指數(shù)是必然趨勢(shì),并提出發(fā)展適合大面積種植、高效益的冬種作物等措施[5]。楊曉光等采用公認(rèn)的農(nóng)業(yè)氣候指標(biāo)計(jì)算方法,使用ArcGIS繪出了1950—1980年和1981—2007年2個(gè)時(shí)段全國(guó)種植制度北界圖,認(rèn)為一年一熟帶、一年兩熟帶、一年三熟帶都不同程度地向北移動(dòng)[6]。張志國(guó)以市為單元計(jì)算了河南省復(fù)種指數(shù)和潛力復(fù)種指數(shù),建立了河南省復(fù)種指數(shù)和糧食產(chǎn)量之間的函數(shù)關(guān)系,認(rèn)為耕地復(fù)種指數(shù)和糧食產(chǎn)量顯著相關(guān),且復(fù)種指數(shù)是糧食產(chǎn)量變化的原因,但河南省再依靠提高復(fù)種指數(shù)來增加糧食產(chǎn)量的余地不大[7]。楊丹等基于空間數(shù)據(jù)分析(ESDA)與GIS技術(shù),探討江西省縣域尺度復(fù)種指數(shù)的空間格局和演變規(guī)律,認(rèn)為江西省的耕地面積、糧食播種面積、耕地復(fù)種指數(shù)均呈現(xiàn)下降的態(tài)勢(shì),但局部空間集聚顯著,有強(qiáng)烈的空間相關(guān)性[8]。何堅(jiān)堅(jiān)等運(yùn)用區(qū)域差異度量模型、空間自相關(guān)模型對(duì)2005—2014年中原經(jīng)濟(jì)區(qū)復(fù)種指數(shù)的時(shí)空差異和驅(qū)動(dòng)力進(jìn)行研究,認(rèn)為中原經(jīng)濟(jì)區(qū)復(fù)種指數(shù)整體呈現(xiàn)波動(dòng)下降的趨勢(shì),自然環(huán)境因素對(duì)耕地復(fù)種指數(shù)及復(fù)種潛力指數(shù)均起著重要的作用[9]。

    總體上看,基于省市統(tǒng)計(jì)數(shù)據(jù)進(jìn)行全國(guó)和省區(qū)等尺度復(fù)種指數(shù)研究較多,縣域單元尺度的全國(guó)性研究不足,利用重心模型、空間相關(guān)分析等指標(biāo)進(jìn)行時(shí)間維度變化的研究較少。本研究系統(tǒng)收集整理全國(guó)縣域耕地面積、農(nóng)作物播種面積等數(shù)據(jù),分析我國(guó)復(fù)種指數(shù)時(shí)空變化特征,研究我國(guó)縣域復(fù)種指數(shù)的空間差異特征,以期為未來我國(guó)糧食安全研究提供理論支撐。

    1 資料與方法

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

    全國(guó)1985—2017年耕地面積,農(nóng)作物、糧食作物、蔬菜播種面積數(shù)據(jù)來源于國(guó)家統(tǒng)計(jì)局。1985、1995、2005、2015年不同縣域的耕地面積、農(nóng)作物播種面積數(shù)據(jù)來源于農(nóng)業(yè)農(nóng)村部發(fā)布的農(nóng)業(yè)分縣統(tǒng)計(jì)數(shù)據(jù),并結(jié)合縣級(jí)和省級(jí)農(nóng)業(yè)統(tǒng)計(jì)數(shù)據(jù)進(jìn)行校正。縣域行政區(qū)數(shù)字地圖來源于國(guó)家基礎(chǔ)地理信息中心,中國(guó)農(nóng)作制分區(qū)圖(圖1)根據(jù)《中國(guó)農(nóng)作制》[10]繪制。

    1.2 研究方法

    1.2.1 復(fù)種指數(shù)

    復(fù)種指數(shù)反映某一種植制度對(duì)耕地的利用程度,計(jì)算公式[4]如下:

    MCI=ASAC×100%。

    式中:MCI為復(fù)種指數(shù);AS為農(nóng)作物總播種面積,hm2;AC為耕地總面積,hm2。

    1.2.2 復(fù)種指數(shù)變化

    復(fù)種指數(shù)變化程度(K)是指某一空間單元復(fù)種指數(shù)從i時(shí)期到j(luò)時(shí)期發(fā)生的變化。

    K=MCIj-MCIi。

    1.2.3 重心分析法

    借鑒力學(xué)原理,引入?yún)^(qū)域重心的概念來確定區(qū)域復(fù)種指數(shù)重心點(diǎn)在不同年份移動(dòng)的方向和距離,用來對(duì)復(fù)種指數(shù)動(dòng)態(tài)演化過程進(jìn)行分析。

    xj=∑ni=1(Qij×xi)/∑ni=1Qij;

    yj=∑ni=1(Qij×yi)/∑ni=1Qij;

    dm=(xk+m-xk)2+(yk+m-yk)2。

    式中,xj和yj分別為復(fù)種指數(shù)重心的經(jīng)度與緯度;Qij表示復(fù)種指數(shù)在第i個(gè)縣域第j年的數(shù)據(jù);dm表示復(fù)種指數(shù)在不同年份間移動(dòng)距離;yk+m和yk分別表示復(fù)種指數(shù)在第k+m年和k年的緯度坐標(biāo);xk+m和xk分別表示復(fù)種指數(shù)在第k+m年和k年的經(jīng)度坐標(biāo)。P(xj,yj) 代表了研究復(fù)種指數(shù)第j年重心的地理坐標(biāo)。將上述公式在Excel中運(yùn)算之后,將P(xj,yj)在ArcGIS 10.2中進(jìn)行空間顯示。

    1.2.4 全局空間自相關(guān) 全局自相關(guān)是復(fù)種指數(shù)間整體的相關(guān)性,主要指標(biāo)有Morans I指數(shù)、Gearys c系數(shù)和Getis-Ord G指數(shù),本研究采用Morans I指數(shù),計(jì)算公式為:

    Morans I=∑ni=1∑nj=1Wij(yi-Y)(yj-Y)S2∑ni=1∑nj=1Wij;

    S2=1n∑ni=1(yi-Y)2;

    Y=1n∑ni=1yi。

    其中,n為樣本量,即空間單元個(gè)數(shù);yi、yj是空間單元i和j的復(fù)種指數(shù);Wij為二進(jìn)制的鄰接空間權(quán)值矩陣,即空間權(quán)重矩陣,用以度量空間單元間的依賴程度。在構(gòu)造空間權(quán)重矩陣時(shí),常用的是相鄰原則中的ROOK 標(biāo)準(zhǔn),即有共同邊界的原則。則有:

    Wij=f(x)=1,當(dāng)?shù)貐^(qū)i和地區(qū)j相鄰

    0,當(dāng)?shù)貐^(qū)i和地區(qū)j不相鄰。

    Morans I 指數(shù)的取值范圍是 [-1,1],取值大于0時(shí)表明空間存在正相關(guān),小于0時(shí)表明空間存在,等于0時(shí)表明呈獨(dú)立隨機(jī)分布。

    1.2.5 局部空間自相關(guān) Morans I 反映的是空間整體上的聚集或者是離散態(tài)勢(shì),并不能反映局部單元的具體情況。而冷熱點(diǎn)分析則解決了這一問題,它反映的是局部空間的關(guān)聯(lián)程度,可以標(biāo)識(shí)出變量的空間高值區(qū)和低值區(qū),其計(jì)算公式為:

    Gi(d)2=∑nj=1(yi-Y)(yj-Y)S2∑ni=1∑nj=1Wij;

    Gi(d)2=∑nj=1Wij(d)Xj∑nj=1Xj。

    式中:Xj為基本空間單元j的復(fù)種指數(shù),Wij(d)是以空間距離d度量的權(quán)重。結(jié)果為正值的時(shí)候表明周圍屬于熱點(diǎn)區(qū)域,高值聚集;結(jié)果為負(fù)值的時(shí)候表明周圍為冷點(diǎn)區(qū)域,低值聚集。

    1.3 數(shù)據(jù)的統(tǒng)計(jì)分析和空間表達(dá)

    基于Excel軟件對(duì)基礎(chǔ)數(shù)據(jù)進(jìn)行處理,采用ArcGIS 10.0軟件進(jìn)行空間分析與表達(dá)。

    2 結(jié)果與分析

    2.1 全國(guó)復(fù)種指數(shù)及農(nóng)作物播種面積變化分析

    1985年以來,除1996年統(tǒng)計(jì)方法進(jìn)一步完善引起的耕地面積增長(zhǎng)外,我國(guó)耕地面積整體呈減少趨勢(shì)(圖2)。農(nóng)作物播種面積整體呈現(xiàn)波動(dòng)性增長(zhǎng)的趨勢(shì),從1985年的14.36×107 hm2持續(xù)增長(zhǎng)到目前的16.63×107 hm2左右。其中1992—1994和2002—2004年出現(xiàn)波動(dòng)性下降。糧食作物播種面積自1985—1998年呈現(xiàn)波動(dòng)性增長(zhǎng)趨勢(shì),1999年播種面積增加到11.38×107 hm2,到2003年僅有9.94×107 hm2;之后在國(guó)家政策的支持下,穩(wěn)步增長(zhǎng)到2017年的11.80×107 hm2。蔬菜面積呈現(xiàn)整體直線上升的趨勢(shì),從1985年的0.48×107 hm2增加到2017年的2.00×107 hm2,在一定程度上支撐了農(nóng)作物總播種面積的持續(xù)增長(zhǎng)趨勢(shì)。1985—1995年,我國(guó)復(fù)種指數(shù)呈現(xiàn)整體上升趨勢(shì);1999年由于耕地面積統(tǒng)計(jì)方法的變化,復(fù)種指數(shù)急劇下降;2000—2007年略有下降;2008年明顯提升;之后2008—2013年略有提升;2013至今比較平穩(wěn),復(fù)種指數(shù)保持在123%左右。

    2.2 全國(guó)縣域復(fù)種指數(shù)時(shí)空特征研究

    我國(guó)復(fù)種指數(shù)從空間上看呈現(xiàn)從北到南逐漸增加的趨勢(shì),從時(shí)間上看大多數(shù)地區(qū)呈現(xiàn)“上升—降低—上升”的趨勢(shì)。1995年,復(fù)種指數(shù)在201%~300%之間的區(qū)域范圍最大,之后2005年和2015年不斷減少,同時(shí)有向西北方向移動(dòng)的趨勢(shì)(圖3)。

    1985年,復(fù)種指數(shù)≤100%的有389個(gè)縣域,占全國(guó)的18.12%,主要分布在東北農(nóng)作區(qū)興安嶺亞區(qū)、北部低中高原農(nóng)作區(qū)的內(nèi)蒙古高原北部亞區(qū)、西北農(nóng)作區(qū)、青藏高原農(nóng)作區(qū);復(fù)種指數(shù)在101%~200%的有1 345個(gè)縣域,主要分布在黃淮海、西南、東北(南部)、江南丘陵農(nóng)作區(qū),以及長(zhǎng)江中下游、四川盆地農(nóng)作區(qū)局部,占全國(guó)的62.65%;復(fù)種指數(shù)在201%~300%的有394個(gè)縣域,主要分布在長(zhǎng)江中下游、華南(大部)、四川盆地(局部)農(nóng)作區(qū),占18.35%;復(fù)種指數(shù)>300%的有19個(gè)縣域,呈零星點(diǎn)狀分布。

    1995年,復(fù)種指數(shù)≤100%的有359個(gè)縣域,占全國(guó)15.87%,主要分布東北農(nóng)作區(qū)的興安嶺亞區(qū)、北部低中高原農(nóng)作區(qū)的內(nèi)蒙古高原北部亞區(qū)、西北農(nóng)作區(qū)的河套河西亞區(qū)和北疆亞區(qū)、青藏高原農(nóng)作區(qū);復(fù)種指數(shù)在101%~200%的有1 147個(gè)縣域,主要分布在黃淮海、西南(西部)、東北(南部)農(nóng)作區(qū),占50.71%;復(fù)種指數(shù)在201%~300%的有721個(gè)縣域,主要分布在長(zhǎng)江中下游、江南丘陵、華南、四川盆地農(nóng)作區(qū),占31.87%;復(fù)種指數(shù)>300%的有35個(gè)縣域,呈零星點(diǎn)狀分布。

    2005年,復(fù)種指數(shù)≤100%的有551個(gè)縣域,占21.75%,主要分布東北農(nóng)作區(qū)的興安嶺亞區(qū)、北部低中高原農(nóng)作區(qū)、西北農(nóng)作區(qū)的河套河西亞區(qū)和北疆亞區(qū)、青藏高原農(nóng)作區(qū);復(fù)種指數(shù)在101%~200%的有1 307個(gè)縣域,主要分布在黃淮海、西南(西部)、東北(南部)、西南(大部)農(nóng)作區(qū)以及長(zhǎng)江中下游農(nóng)作區(qū)北部,占51.60%;復(fù)種指數(shù)在201%~300%的有599個(gè)縣域,主要分布在長(zhǎng)江中下游農(nóng)作區(qū)的兩湖平原亞區(qū)、江南丘陵、華南、四川盆地農(nóng)作區(qū),占23.65%;復(fù)種指數(shù)>300%的有76個(gè)縣域,呈零星點(diǎn)狀分布。

    2015年,復(fù)種指數(shù)≤100%的有565個(gè)縣域,占22.31%,主要分布東北農(nóng)作區(qū)的興安嶺丘陵山地亞區(qū)、北部低中高原農(nóng)作區(qū)大部、青藏高原農(nóng)作區(qū);復(fù)種指數(shù)在101%~200%的有1 248個(gè)縣域,主要分布在黃淮海、西南(南部)、東北(南部)、西南(大部)、長(zhǎng)江中下游(北部)、江南丘陵、華南(局部)農(nóng)作區(qū),占49.29%;復(fù)種指數(shù)在201%~300%的636個(gè)縣域,主要分布在長(zhǎng)江中下游、江南丘陵(北部)、四川盆地農(nóng)作區(qū),占25.12%;復(fù)種指數(shù)>300%的有83個(gè)縣域,呈零星點(diǎn)狀分布(表1)。

    2.2 全國(guó)縣域復(fù)種指數(shù)變化研究

    1985—1995年,全國(guó)復(fù)種指數(shù)整體呈現(xiàn)上升趨勢(shì),其中復(fù)種指數(shù)變化在1%~20%的縣域有867個(gè),僅占全國(guó)的41.36%,復(fù)種指數(shù)變化>20%的縣域有691個(gè),僅占全國(guó)的32.97%。復(fù)種指數(shù)降低的縣域主要分布在北部低中高原、西北、青藏高原農(nóng)作區(qū),以及長(zhǎng)江中下游農(nóng)作區(qū)東部的局部地區(qū),其中復(fù)種指數(shù)變化≤-20%的縣域有123個(gè),僅占全國(guó)的5.84%;復(fù)種指數(shù)變化在-21%~0%之間的縣域有415個(gè),僅占全國(guó)的19.80%。

    1995—2005年,復(fù)種指數(shù)上升的區(qū)域主要分布在東北、黃淮海、 西南(南部)農(nóng)作區(qū),復(fù)種指數(shù)降低

    的縣域主要分布在長(zhǎng)江中下游、江南丘陵、華南、青藏高原農(nóng)作區(qū),其中復(fù)種指數(shù)變化集中于1%~20%的縣域有542個(gè),僅占全國(guó)的24.22%;復(fù)種指數(shù)變化>20%的縣域有456個(gè),僅占全國(guó)的20.38%;復(fù)種指數(shù)變化≤-20%的縣域有728個(gè),僅占全國(guó)的32.53%;復(fù)種指數(shù)變化在-21%~0之間的縣域有812個(gè),占全國(guó)的36.28%。

    2005—2015年,全國(guó)復(fù)種指數(shù)變化空間差異較大,整體上升區(qū)域主要集中在東北、黃淮海、長(zhǎng)江中下游(北部)農(nóng)作區(qū),整體下降區(qū)域主要集中在長(zhǎng)江中下游(南部)、江南丘陵、華南農(nóng)作區(qū)。其中復(fù)種指數(shù)變化在1%~20%的縣域有478個(gè),僅占全國(guó)的20.71%;復(fù)種指數(shù)變化>20%的縣域有657個(gè),僅占全國(guó)的28.47%;復(fù)種指數(shù)變化≤-20%的縣域有709個(gè),僅占全國(guó)的30.72%;復(fù)種指數(shù)變化在 -21%~0之間的縣域有464個(gè),占全國(guó)的20.10%。

    1985—2015年,全國(guó)復(fù)種指數(shù)上升區(qū)域主要集中在東北、黃淮海、西南、四川盆地農(nóng)作區(qū),整體下降區(qū)域主要集中在北部低中高原、長(zhǎng)江中下游(東南部)、江南丘陵、華南農(nóng)作區(qū)。其中復(fù)種指數(shù)變化在1%~20%的縣域有381個(gè),僅占全國(guó)的18.05%;復(fù)種指數(shù)變化>20%的縣域有782個(gè),僅占全國(guó)的37.04%;復(fù)種指數(shù)變化≤-20%的縣域有617個(gè),僅占全國(guó)的29.23%,復(fù)種指數(shù)變化在 -21%~0之間的縣域有331個(gè),占全國(guó)的15.68%(圖4、表2)。

    2.3 全國(guó)縣域復(fù)種指數(shù)空間差異分析

    2.3.1 總體空間差異

    由表3可見,1985、1995、2005、2015年的全局Morans I分別為0.356 3、0.484 0、0.601 2、0.544 1,在4個(gè)年份數(shù)據(jù)中Z值>2.58,且都為正值,P值皆小于0.01,表示此數(shù)據(jù)并不是隨機(jī)生成的,在99.9%的置信度情況下成顯著正相關(guān),說明我國(guó)復(fù)種指數(shù)存在空間上的自相關(guān)性。4個(gè)年份中1985年復(fù)種指數(shù)與與空間的相關(guān)性最小;1985—2005年的空間相關(guān)性呈現(xiàn)增長(zhǎng)趨勢(shì),表明我國(guó)縣域的復(fù)種指數(shù)在空間上的聚集態(tài)勢(shì)增加;2005—2015年相關(guān)性降低;與1985年相比,2005年的空間相關(guān)性達(dá)到最大值,表明復(fù)種指數(shù)在空間上的聚集態(tài)勢(shì)達(dá)到1985—2005年近20年來的最高狀態(tài)。

    2.3.2 局部空間差異

    為了更加直觀地表示我國(guó)復(fù)種指數(shù)在局部空間的聚集和擴(kuò)散狀態(tài),對(duì)1985、1995、2005、2015年4個(gè)時(shí)期全國(guó)縣域復(fù)種指數(shù)的冷熱點(diǎn)進(jìn)行分析,生成全國(guó)地級(jí)尺度的冷熱點(diǎn)空間分布圖(圖5),熱點(diǎn)區(qū)域表示復(fù)種指數(shù)高值區(qū)域聚集,冷點(diǎn)區(qū)域表示復(fù)種指數(shù)低值區(qū)域聚集。

    由圖5可見,1985年我國(guó)復(fù)種指數(shù)的熱點(diǎn)區(qū)主要分布在長(zhǎng)江中下游(南部)、江南丘陵(中部)、華南農(nóng)作區(qū),冷點(diǎn)地區(qū)分布在東北、北部低中高原、西北、青藏高原農(nóng)作區(qū)。1995年我國(guó)復(fù)種指數(shù)的熱點(diǎn)區(qū)域進(jìn)一步擴(kuò)大,江南丘陵、華南、西南(東部)農(nóng)作區(qū)均有不同程度增加,呈現(xiàn)出明顯的空間聚集態(tài)勢(shì)。2005年,冷熱點(diǎn)聚集程度和分布均有所下降,其中熱點(diǎn)區(qū)域增加了四川盆地農(nóng)作區(qū),減少了沿海地區(qū)農(nóng)作區(qū);冷點(diǎn)區(qū)域中東北、西北、青藏高原農(nóng)作區(qū)空間聚集程度有所降低。2015年復(fù)種指數(shù)的空間聚集程度進(jìn)一步降低,熱點(diǎn)區(qū)域向西移動(dòng),主要分布在長(zhǎng)江中下游(西部)、西南(東部、南部)、四川盆地農(nóng)作區(qū),江南丘陵(中部)、華南農(nóng)作區(qū)不再是熱點(diǎn)區(qū)域;冷點(diǎn)地區(qū)分布在東北、北部低中高原、西北(東部)農(nóng)作區(qū),西北(西部)、青藏高原農(nóng)作區(qū)不再是冷點(diǎn)區(qū)域。

    2.4 全國(guó)復(fù)種指數(shù)重心遷移

    我國(guó)復(fù)種指數(shù)重心1985年位于湖北省樊城區(qū),1995年向南移動(dòng)至南漳縣中部,2005年向北移動(dòng)至南漳縣北部,2015年繼續(xù)北上移動(dòng)至谷城縣;2015年與1985年相比,重心向西南移動(dòng)33 km。復(fù)種指數(shù)≤100%的重心在1985年位于內(nèi)蒙古鄂托克旗,1995年北移至杭錦旗,2005年向東南移動(dòng)至內(nèi)蒙南部與山西北部交界神木縣最北部, 2015年繼續(xù)向東南移動(dòng)至陜西省祈府區(qū);2015年與1985年相比,重心向東南移動(dòng)433 km。復(fù)種指數(shù)101%~200%的重心在1985年位于湖北省與河南省交界的丹江口,1995年北移至河南省洛寧縣,2005向東南移動(dòng)至南召縣,2015年在南召縣內(nèi)又稍微向西移動(dòng);2015年與1985年相比,重心向東北移動(dòng)104 km。復(fù)種指數(shù)>200%重心在1985年位于湖南最東部的瀏陽市,1995年向西南移動(dòng)至湘鄉(xiāng)市,2005年向西北移動(dòng)至新化縣,2015年繼續(xù)向西北移動(dòng)至湖南北部的慈利縣;2015年與1985年相比,重心向東北移動(dòng)319 km。

    3 討論與結(jié)論

    復(fù)種是在1塊地上1年內(nèi)種收1季以上作物的種植方式。復(fù)種指數(shù)作為農(nóng)作制中衡量耕地集約化利用程度的重要指標(biāo),既是農(nóng)學(xué)學(xué)科的研究范疇,也是地理學(xué)科的研究范疇,具有多學(xué)科交叉的特點(diǎn)[2]?;谛姓卧y(tǒng)計(jì)數(shù)據(jù)的復(fù)種指數(shù)研究方法具有數(shù)據(jù)容易獲取、計(jì)算方便等優(yōu)勢(shì),在國(guó)家、區(qū)域、省等不同尺度范圍得到了廣泛應(yīng)用。近年來,隨著對(duì)地觀測(cè)技術(shù)的快速發(fā)展,利用遙感技術(shù)監(jiān)測(cè)從生長(zhǎng)到成熟、衰落整個(gè)作物周期的活動(dòng)過程[11-18]受到越來越多的研究人員關(guān)注,復(fù)種指數(shù)監(jiān)測(cè)精度達(dá)到了地塊尺度,可以更好地反映復(fù)種指數(shù)的真實(shí)情況。鑒于統(tǒng)計(jì)單元內(nèi)部的空間異質(zhì)性,統(tǒng)計(jì)方法不能準(zhǔn)確描述種植制度的空間特征,在數(shù)據(jù)獲取上也存在一定的滯后性[19],加上遙感技術(shù)在監(jiān)測(cè)精度、工作量等方面的不足,筆者認(rèn)為融合遙感監(jiān)測(cè)技術(shù)與統(tǒng)計(jì)方法,進(jìn)一步提高復(fù)種指數(shù)監(jiān)測(cè)的精準(zhǔn)性、及時(shí)性是未來發(fā)展的重點(diǎn)。

    本研究基于縣域統(tǒng)計(jì)數(shù)據(jù)對(duì)實(shí)際復(fù)種指數(shù)的時(shí)空特征進(jìn)行了分析,認(rèn)為從時(shí)間上大多數(shù)地區(qū)呈現(xiàn)“上升—降低—上升”的趨勢(shì),這與李卓等對(duì)黃淮海地區(qū)耕地復(fù)種指數(shù)的研究結(jié)論[20]一致。研究認(rèn)為1995年我國(guó)復(fù)種指數(shù)在201%~300%的區(qū)域范圍最大,之后2005年和2015年不斷減少,這與蔣敏等研究的趨勢(shì)[7-8,21]一致。

    本研究結(jié)果表明,我國(guó)復(fù)種指數(shù)在空間上存在較強(qiáng)的自相關(guān)性,具有空間上聚集的特征,熱點(diǎn)從東南部向西南部轉(zhuǎn)移,冷點(diǎn)在北方地區(qū)范圍縮減,同時(shí)復(fù)種指數(shù)重心發(fā)生變化,這可能與我國(guó)氣候變暖而導(dǎo)致的種植制度北界北移[6]相關(guān),也可能與當(dāng)?shù)氐纳鐣?huì)經(jīng)濟(jì)條件、農(nóng)業(yè)生產(chǎn)條件、農(nóng)業(yè)種植結(jié)構(gòu)等有關(guān),需要進(jìn)一步研究以明確。

    參考文獻(xiàn):

    [1]黃國(guó)勤,孫丹平. 中國(guó)多熟種植的發(fā)展現(xiàn)狀與研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào),2017,33(3):35-43.

    [2]吳文斌,余強(qiáng)毅,陸 苗,等. 耕地復(fù)種指數(shù)研究的關(guān)鍵科學(xué)問題[J]. 中國(guó)農(nóng)業(yè)科學(xué),2018,51(9):1681-1694.

    [3]黃國(guó)勤. 論提高我國(guó)南方耕地復(fù)種指數(shù)[J]. 自然資源,1995,17(1):30-37.

    [4]劉巽浩. 論我國(guó)耕地種植指數(shù)(復(fù)種)的潛力[J]. 作物雜志,1997(3):1-3.

    [5]段紅平. 湖南省耕作制度50年(1949—1998年)演變分析[J]. 耕作與栽培,2001(3):1-4.

    [6]楊曉光,劉志娟,陳 阜. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅰ.氣候變暖對(duì)中國(guó)種植制度北界和糧食產(chǎn)量可能影響的分析[J]. 中國(guó)農(nóng)業(yè)科學(xué),2010,43(2):329-336.

    [7]張志國(guó). 河南省復(fù)種指數(shù)的時(shí)空特征及其與糧食產(chǎn)量關(guān)系研究[J]. 湖北農(nóng)業(yè)科學(xué),2011,50(17):3653-3656.

    [8]楊 丹,葉長(zhǎng)盛. 基于縣域單元的江西省復(fù)種指數(shù)時(shí)空格局分析[J]. 農(nóng)業(yè)現(xiàn)代化研究,2013,34(6):722-727.

    [9]何堅(jiān)堅(jiān),龐 博,張鵬巖,等. 區(qū)域耕地復(fù)種指數(shù)時(shí)空差異測(cè)算及可挖掘潛力分析——以中原經(jīng)濟(jì)區(qū)為例[J]. 江蘇農(nóng)業(yè)科學(xué),2018,46(16):249-255.

    [10]劉巽浩,陳 阜. 中國(guó)農(nóng)作制[M]. 北京:中國(guó)農(nóng)業(yè)出版社,2005.

    [11]Panigrahy S,Manjunath K R,Ray S S. Deriving cropping system performance indices using remote sensing data and GIS[J]. International Journal of Remote Sensing,2005,26(12):2595-2606.

    [12]Frolking S,Yeluripati J B,Douglas E. New district-level maps of rice cropping in India:A foundation for scientific input into policy assessment[J]. Field Crop Research,2006,98(2):164-177.

    [13]Sakamoto T,van Nguyen N,Ohno H,et al. Spatio-temporal distribution of rice phenology and cropping systems in the Mekong Delta with special reference to the seasonal water flow of the Mekong and Bassac rivers[J]. Remote Sensing of Environment,2006,100(1):1-16.

    [14]彭代亮,黃敬峰,金輝民. 基于MODIS-NDVI的浙江省耕地復(fù)種指數(shù)監(jiān)測(cè)[J]. 中國(guó)農(nóng)業(yè)科學(xué),2006,39(7):1352-1357.

    [15]Canisius F,Turral H,Molden D. Fourier analysis of historical NOAA time series data to estimate bimodal agriculture[J]. International Journal of Remote Sensing,2007,28(24):5503-5522.

    [16]唐鵬欽,姚艷敏,吳文斌,等. 基于遙感技術(shù)的耕地復(fù)種指數(shù)研究進(jìn)展[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃,2010,31(2):21-27.

    [17]楊 婷,趙文利,王哲怡,等. 基于遙感影像NDVI數(shù)據(jù)的中國(guó)種植制度分布變化[J]. 中國(guó)農(nóng)業(yè)科學(xué),2015,48(10):1915-1925.

    [18]Yan H,Xiao X,Huang H,et al. Multiple cropping intensity in China derived from agro-meteorological observation and MODIS data[J]. Chinese Geographical Science,2014,24(2):205-219.

    [19]閆慧敏,曹明奎,劉紀(jì)遠(yuǎn),等. 基于多時(shí)相遙感信息的中國(guó)農(nóng)業(yè)種植制度空間格局研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2005,21(4):85-90.

    [20]李 卓,劉淑亮,孫然好,等. 黃淮海地區(qū)耕地復(fù)種指數(shù)的時(shí)空格局演變[J]. 生態(tài)學(xué)報(bào),2018,38(12):4454-4460.

    [21]蔣 敏,李秀彬,辛良杰,等. 南方水稻復(fù)種指數(shù)變化對(duì)國(guó)家糧食產(chǎn)能的影響及其政策啟示[J]. 地理學(xué)報(bào),2019,74(1):32-43.

    猜你喜歡
    重心時(shí)空變化縣域
    縣域消防專項(xiàng)規(guī)劃研究
    山東縣域GDP排名出爐
    商周刊(2017年6期)2017-08-22 03:42:46
    英漢句子結(jié)構(gòu)差異淺析
    人間(2016年27期)2016-11-11 16:35:51
    網(wǎng)球比賽中運(yùn)動(dòng)員的身體重心變化及步伐特征研究
    護(hù)理專業(yè)的學(xué)生必須要掌握的物理知識(shí)
    基于主成分分析法及熵值法結(jié)合的耕地集約利用時(shí)空變化
    商(2016年29期)2016-10-29 11:10:35
    江蘇經(jīng)濟(jì)發(fā)展區(qū)域差異時(shí)空變化研究
    南海聲速主躍層的時(shí)空分布
    縣域就診率為何差了40%
    一項(xiàng)縣域中招政策的回顧與反思
    久久午夜福利片| 日本wwww免费看| 内地一区二区视频在线| 精品99又大又爽又粗少妇毛片| 国产又爽黄色视频| 欧美激情国产日韩精品一区| 一级片免费观看大全| 亚洲av男天堂| 少妇猛男粗大的猛烈进出视频| 国产片特级美女逼逼视频| 纯流量卡能插随身wifi吗| 久久亚洲国产成人精品v| 亚洲av日韩在线播放| 国产麻豆69| 啦啦啦中文免费视频观看日本| 国产精品一二三区在线看| 在线观看免费视频网站a站| 国产有黄有色有爽视频| 日韩中文字幕视频在线看片| 成年人午夜在线观看视频| 久久精品国产鲁丝片午夜精品| 亚洲,欧美精品.| 日韩一本色道免费dvd| 国产男人的电影天堂91| 九色亚洲精品在线播放| av视频免费观看在线观看| 亚洲一码二码三码区别大吗| 亚洲欧洲日产国产| 免费黄色在线免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产xxxxx性猛交| 黑丝袜美女国产一区| 免费播放大片免费观看视频在线观看| videossex国产| 精品99又大又爽又粗少妇毛片| 国产熟女午夜一区二区三区| 日本av手机在线免费观看| 国产精品久久久久久精品电影小说| 丝袜喷水一区| 永久免费av网站大全| 色吧在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av在线观看美女高潮| 少妇高潮的动态图| 99热网站在线观看| 一二三四中文在线观看免费高清| 久久精品国产鲁丝片午夜精品| 久久女婷五月综合色啪小说| 爱豆传媒免费全集在线观看| 国产精品久久久久久精品电影小说| 日韩视频在线欧美| 亚洲一区二区三区欧美精品| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩另类电影网站| 免费黄频网站在线观看国产| 亚洲国产精品999| 欧美最新免费一区二区三区| 成人国语在线视频| 亚洲综合精品二区| 久久精品aⅴ一区二区三区四区 | 亚洲精品久久久久久婷婷小说| 亚洲伊人久久精品综合| 777米奇影视久久| 黄片播放在线免费| 热re99久久精品国产66热6| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 日日摸夜夜添夜夜爱| 亚洲精品456在线播放app| 欧美激情国产日韩精品一区| 精品久久久久久电影网| 黑人欧美特级aaaaaa片| 人妻人人澡人人爽人人| 欧美丝袜亚洲另类| 久久精品国产亚洲av涩爱| 国产精品成人在线| 伦理电影免费视频| 校园人妻丝袜中文字幕| 咕卡用的链子| 免费黄网站久久成人精品| 女性生殖器流出的白浆| av又黄又爽大尺度在线免费看| 国产精品熟女久久久久浪| videos熟女内射| 看免费成人av毛片| 寂寞人妻少妇视频99o| 亚洲国产av新网站| 欧美人与性动交α欧美软件 | 一边亲一边摸免费视频| 免费看av在线观看网站| 在线天堂中文资源库| 欧美bdsm另类| 日本欧美国产在线视频| 日韩在线高清观看一区二区三区| 国产一区二区三区av在线| 国产综合精华液| 美女视频免费永久观看网站| 日本免费在线观看一区| 久久影院123| 国产日韩欧美视频二区| 亚洲三级黄色毛片| 九草在线视频观看| 日本wwww免费看| 日韩免费高清中文字幕av| 乱码一卡2卡4卡精品| 亚洲精品美女久久av网站| 久久国产精品大桥未久av| 超色免费av| 9色porny在线观看| 18禁动态无遮挡网站| 久久久国产精品麻豆| 十分钟在线观看高清视频www| 99久久精品国产国产毛片| 精品国产一区二区三区四区第35| 中文欧美无线码| 亚洲欧美日韩卡通动漫| av在线老鸭窝| 久久99精品国语久久久| 亚洲高清免费不卡视频| 精品一区二区三卡| 69精品国产乱码久久久| 国内精品宾馆在线| 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 亚洲,一卡二卡三卡| 国产在线一区二区三区精| 国产精品久久久久久久电影| 国产精品三级大全| 亚洲精品国产av蜜桃| 全区人妻精品视频| 2018国产大陆天天弄谢| 久久久久久久大尺度免费视频| 波野结衣二区三区在线| 精品一区在线观看国产| 国产欧美日韩一区二区三区在线| 久久国产精品大桥未久av| 中文精品一卡2卡3卡4更新| 狂野欧美激情性bbbbbb| 波野结衣二区三区在线| 欧美精品av麻豆av| av在线app专区| 久久久久久久久久久久大奶| 日日摸夜夜添夜夜爱| 极品少妇高潮喷水抽搐| 免费高清在线观看视频在线观看| 一本久久精品| 熟女av电影| 国产欧美日韩综合在线一区二区| 曰老女人黄片| 国产成人精品福利久久| 日韩av免费高清视频| 下体分泌物呈黄色| 亚洲国产欧美日韩在线播放| 欧美精品一区二区免费开放| 精品久久久久久电影网| 2018国产大陆天天弄谢| 99久国产av精品国产电影| 韩国av在线不卡| 老熟女久久久| av黄色大香蕉| 久久毛片免费看一区二区三区| 国产精品.久久久| 精品国产露脸久久av麻豆| 91aial.com中文字幕在线观看| 日本vs欧美在线观看视频| 欧美日韩视频高清一区二区三区二| 国产黄色视频一区二区在线观看| 成人18禁高潮啪啪吃奶动态图| 国产精品蜜桃在线观看| 男的添女的下面高潮视频| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 亚洲精品自拍成人| 成人无遮挡网站| 精品一品国产午夜福利视频| 精品一区二区免费观看| 亚洲欧美一区二区三区国产| 黑丝袜美女国产一区| 国国产精品蜜臀av免费| 伊人久久国产一区二区| av播播在线观看一区| 美女内射精品一级片tv| 18禁国产床啪视频网站| 一本色道久久久久久精品综合| 日韩 亚洲 欧美在线| 国产欧美亚洲国产| 欧美国产精品一级二级三级| 亚洲少妇的诱惑av| 日本av免费视频播放| 99九九在线精品视频| 在线观看免费视频网站a站| 99香蕉大伊视频| 在线 av 中文字幕| videosex国产| 久久影院123| 十八禁高潮呻吟视频| 视频在线观看一区二区三区| 亚洲成人手机| 国产精品人妻久久久久久| 国产麻豆69| 侵犯人妻中文字幕一二三四区| 成人亚洲精品一区在线观看| 精品第一国产精品| 国产亚洲av片在线观看秒播厂| 最新的欧美精品一区二区| 卡戴珊不雅视频在线播放| 亚洲丝袜综合中文字幕| 亚洲一级一片aⅴ在线观看| 国产亚洲一区二区精品| 高清黄色对白视频在线免费看| 久久精品久久精品一区二区三区| 80岁老熟妇乱子伦牲交| 亚洲一级一片aⅴ在线观看| 在线观看www视频免费| 99视频精品全部免费 在线| 制服人妻中文乱码| 日韩一区二区视频免费看| 日韩中字成人| a级毛色黄片| 肉色欧美久久久久久久蜜桃| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 欧美97在线视频| 少妇被粗大猛烈的视频| 啦啦啦啦在线视频资源| 天天躁夜夜躁狠狠久久av| 久久韩国三级中文字幕| 人人妻人人爽人人添夜夜欢视频| 国产精品国产三级国产av玫瑰| 亚洲精品456在线播放app| 在线看a的网站| av国产精品久久久久影院| 亚洲精品久久成人aⅴ小说| 亚洲欧洲国产日韩| 黄网站色视频无遮挡免费观看| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美日韩在线播放| 波野结衣二区三区在线| 日日撸夜夜添| 日本av手机在线免费观看| 纯流量卡能插随身wifi吗| 国产精品.久久久| 精品一区二区三区四区五区乱码 | av一本久久久久| 51国产日韩欧美| 亚洲,欧美,日韩| 少妇被粗大猛烈的视频| xxxhd国产人妻xxx| 久久久久久人妻| 亚洲精品视频女| 91精品伊人久久大香线蕉| 又黄又爽又刺激的免费视频.| 成人亚洲欧美一区二区av| 视频区图区小说| 国产1区2区3区精品| 久久这里有精品视频免费| 亚洲综合色惰| 狠狠精品人妻久久久久久综合| 这个男人来自地球电影免费观看 | 少妇猛男粗大的猛烈进出视频| 亚洲欧美成人精品一区二区| 日韩大片免费观看网站| 十八禁高潮呻吟视频| 国产一区二区三区av在线| 人妻少妇偷人精品九色| 欧美xxⅹ黑人| 丝袜美足系列| 黄色配什么色好看| 在线观看美女被高潮喷水网站| 国产色婷婷99| 大码成人一级视频| 一个人免费看片子| 中文欧美无线码| 熟女人妻精品中文字幕| videosex国产| 日韩中文字幕视频在线看片| 男人爽女人下面视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 男女边吃奶边做爰视频| 亚洲国产精品国产精品| 美女主播在线视频| 国产成人欧美| 日本色播在线视频| 性色av一级| 男的添女的下面高潮视频| 亚洲精品国产av成人精品| 九草在线视频观看| 亚洲,一卡二卡三卡| 欧美老熟妇乱子伦牲交| www.熟女人妻精品国产 | 人人澡人人妻人| 久久久久久久亚洲中文字幕| 汤姆久久久久久久影院中文字幕| 人人妻人人添人人爽欧美一区卜| 黑人巨大精品欧美一区二区蜜桃 | 成年av动漫网址| 欧美日韩成人在线一区二区| 人妻系列 视频| 国产女主播在线喷水免费视频网站| 久久午夜福利片| 黄色一级大片看看| xxxhd国产人妻xxx| 欧美变态另类bdsm刘玥| 国产国语露脸激情在线看| 91精品国产国语对白视频| 婷婷色av中文字幕| 久久久国产精品麻豆| 亚洲国产精品一区三区| 国产成人午夜福利电影在线观看| 国产片内射在线| 青春草视频在线免费观看| 国产成人91sexporn| 如何舔出高潮| 国产黄频视频在线观看| 三级国产精品片| 男女边摸边吃奶| 色婷婷久久久亚洲欧美| 国产一区二区三区综合在线观看 | 男女免费视频国产| 亚洲一级一片aⅴ在线观看| 亚洲五月色婷婷综合| 高清黄色对白视频在线免费看| 亚洲综合精品二区| av.在线天堂| 亚洲国产最新在线播放| 少妇的逼水好多| 日本与韩国留学比较| 色视频在线一区二区三区| 亚洲欧美色中文字幕在线| 亚洲精品aⅴ在线观看| 两个人看的免费小视频| 乱码一卡2卡4卡精品| 日韩人妻精品一区2区三区| 亚洲色图综合在线观看| 日韩成人伦理影院| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 亚洲性久久影院| 七月丁香在线播放| 老司机亚洲免费影院| a级片在线免费高清观看视频| 九色成人免费人妻av| 夜夜爽夜夜爽视频| 久久久国产一区二区| 中文字幕亚洲精品专区| 国产福利在线免费观看视频| 青春草亚洲视频在线观看| 天天躁夜夜躁狠狠躁躁| 高清黄色对白视频在线免费看| 成人18禁高潮啪啪吃奶动态图| 青春草视频在线免费观看| 亚洲国产欧美在线一区| 韩国精品一区二区三区 | 亚洲综合精品二区| 另类亚洲欧美激情| 国产精品免费大片| 成人午夜精彩视频在线观看| 搡老乐熟女国产| 国产xxxxx性猛交| 在线观看三级黄色| 久久这里有精品视频免费| 国产精品免费大片| kizo精华| 免费大片18禁| 久久久久久久久久成人| 一区二区三区四区激情视频| 亚洲中文av在线| 久久青草综合色| 天天影视国产精品| 免费高清在线观看视频在线观看| 欧美日韩国产mv在线观看视频| 夫妻性生交免费视频一级片| 99九九在线精品视频| 色网站视频免费| 国语对白做爰xxxⅹ性视频网站| 精品国产露脸久久av麻豆| 一本久久精品| 纵有疾风起免费观看全集完整版| 97精品久久久久久久久久精品| 亚洲精品自拍成人| 欧美 亚洲 国产 日韩一| 女人久久www免费人成看片| 欧美国产精品一级二级三级| 美女中出高潮动态图| 中文精品一卡2卡3卡4更新| 免费人成在线观看视频色| 又粗又硬又长又爽又黄的视频| 久久国产精品男人的天堂亚洲 | 91精品三级在线观看| 免费观看性生交大片5| 国产成人精品福利久久| 一级毛片我不卡| 国产av码专区亚洲av| 亚洲,一卡二卡三卡| a级片在线免费高清观看视频| 国产欧美另类精品又又久久亚洲欧美| 精品人妻一区二区三区麻豆| 777米奇影视久久| 久久婷婷青草| 亚洲av电影在线观看一区二区三区| 亚洲精品一二三| 青春草亚洲视频在线观看| 亚洲成av片中文字幕在线观看 | 亚洲四区av| 天天操日日干夜夜撸| 国产黄色视频一区二区在线观看| 在线观看免费日韩欧美大片| 久热久热在线精品观看| 精品一区二区三区四区五区乱码 | 亚洲欧洲国产日韩| 国产高清三级在线| 一级毛片电影观看| 十八禁网站网址无遮挡| 久久午夜福利片| 亚洲第一区二区三区不卡| 女性被躁到高潮视频| 天天操日日干夜夜撸| 久久人人爽人人片av| 91国产中文字幕| 久久久久久人人人人人| 桃花免费在线播放| 在线观看人妻少妇| 日韩制服骚丝袜av| 如何舔出高潮| 人成视频在线观看免费观看| 美女内射精品一级片tv| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 高清在线视频一区二区三区| av免费观看日本| 麻豆乱淫一区二区| 国产成人精品无人区| 久久精品国产综合久久久 | 国产精品久久久久久精品电影小说| 视频在线观看一区二区三区| 女的被弄到高潮叫床怎么办| 日韩大片免费观看网站| 亚洲图色成人| av福利片在线| 美女视频免费永久观看网站| 日韩一本色道免费dvd| 国产精品一二三区在线看| 看非洲黑人一级黄片| 乱码一卡2卡4卡精品| 午夜久久久在线观看| 免费人成在线观看视频色| 成人18禁高潮啪啪吃奶动态图| 九草在线视频观看| 大香蕉97超碰在线| 欧美日韩成人在线一区二区| 极品人妻少妇av视频| 韩国av在线不卡| av免费观看日本| 久久精品国产a三级三级三级| 99香蕉大伊视频| 内地一区二区视频在线| 久久人人爽人人爽人人片va| 宅男免费午夜| 啦啦啦中文免费视频观看日本| 国产av精品麻豆| av不卡在线播放| 亚洲精品日韩在线中文字幕| 成年av动漫网址| 26uuu在线亚洲综合色| 欧美亚洲日本最大视频资源| 亚洲一区二区三区欧美精品| 欧美国产精品va在线观看不卡| 国产毛片在线视频| 国产精品三级大全| 少妇被粗大的猛进出69影院 | 日韩电影二区| 国产综合精华液| 成人亚洲欧美一区二区av| 国产精品一区www在线观看| 一区二区三区精品91| 九九爱精品视频在线观看| 爱豆传媒免费全集在线观看| 尾随美女入室| 国产免费一区二区三区四区乱码| 亚洲精品aⅴ在线观看| 超碰97精品在线观看| 国产一区二区在线观看av| 久久精品熟女亚洲av麻豆精品| 久久热在线av| 亚洲人成77777在线视频| 国产精品久久久久久精品电影小说| 中文欧美无线码| 极品少妇高潮喷水抽搐| 中文字幕最新亚洲高清| 夫妻性生交免费视频一级片| 51国产日韩欧美| 26uuu在线亚洲综合色| 日韩视频在线欧美| 最近中文字幕高清免费大全6| 中文字幕人妻丝袜制服| 80岁老熟妇乱子伦牲交| 啦啦啦视频在线资源免费观看| 亚洲欧美一区二区三区黑人 | 巨乳人妻的诱惑在线观看| 久久久国产一区二区| 美女中出高潮动态图| 国产精品欧美亚洲77777| 成年动漫av网址| 男女边吃奶边做爰视频| 五月开心婷婷网| 亚洲国产精品成人久久小说| 日本色播在线视频| 免费少妇av软件| 精品第一国产精品| 一级毛片黄色毛片免费观看视频| 99热6这里只有精品| 亚洲欧美成人综合另类久久久| 夜夜爽夜夜爽视频| 久久精品久久精品一区二区三区| 国产男女超爽视频在线观看| av电影中文网址| 全区人妻精品视频| 99久久精品国产国产毛片| 制服丝袜香蕉在线| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人 | 亚洲精品美女久久av网站| 免费日韩欧美在线观看| 五月伊人婷婷丁香| 久久精品熟女亚洲av麻豆精品| 91在线精品国自产拍蜜月| 亚洲第一区二区三区不卡| 成人二区视频| 大片电影免费在线观看免费| 国产一区二区三区综合在线观看 | 欧美日韩视频精品一区| 国产成人精品一,二区| 黄片播放在线免费| 寂寞人妻少妇视频99o| 丝袜在线中文字幕| 亚洲内射少妇av| 看十八女毛片水多多多| 免费黄色在线免费观看| 亚洲av.av天堂| 亚洲精华国产精华液的使用体验| 人妻人人澡人人爽人人| 亚洲美女黄色视频免费看| 天天影视国产精品| 亚洲精品av麻豆狂野| 精品福利永久在线观看| 黑人巨大精品欧美一区二区蜜桃 | 国产成人精品在线电影| 国产日韩欧美视频二区| 亚洲成国产人片在线观看| 精品少妇黑人巨大在线播放| 在线 av 中文字幕| 亚洲欧洲国产日韩| 亚洲精华国产精华液的使用体验| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 精品一品国产午夜福利视频| 亚洲高清免费不卡视频| 蜜臀久久99精品久久宅男| tube8黄色片| 久久精品久久精品一区二区三区| 男的添女的下面高潮视频| av线在线观看网站| 久久热在线av| 最近最新中文字幕大全免费视频 | 色哟哟·www| 久久精品人人爽人人爽视色| 人妻系列 视频| 精品久久蜜臀av无| 女人久久www免费人成看片| 在线免费观看不下载黄p国产| 亚洲av.av天堂| 久久精品国产亚洲av涩爱| 亚洲国产精品成人久久小说| 视频区图区小说| 18禁国产床啪视频网站| 欧美国产精品va在线观看不卡| 26uuu在线亚洲综合色| 国产又色又爽无遮挡免| 欧美 亚洲 国产 日韩一| 欧美bdsm另类| 18禁国产床啪视频网站| 在现免费观看毛片| 一区二区av电影网| 伦精品一区二区三区| 国产精品久久久久久精品电影小说| 日韩制服骚丝袜av| 99国产精品免费福利视频| 亚洲精品日本国产第一区| 在线观看免费日韩欧美大片| 高清毛片免费看| 五月玫瑰六月丁香| 热99久久久久精品小说推荐| 少妇猛男粗大的猛烈进出视频| 免费女性裸体啪啪无遮挡网站| 人妻少妇偷人精品九色| 街头女战士在线观看网站| 午夜影院在线不卡| 人人澡人人妻人| 母亲3免费完整高清在线观看 | 赤兔流量卡办理| 国语对白做爰xxxⅹ性视频网站| 哪个播放器可以免费观看大片| h视频一区二区三区| 99九九在线精品视频| 免费观看a级毛片全部| 久久毛片免费看一区二区三区| 高清欧美精品videossex| 又大又黄又爽视频免费| 亚洲性久久影院| 水蜜桃什么品种好|