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

    氣候變暖和覆膜對(duì)新疆不同熟性棉花種植區(qū)劃的影響

    2019-02-25 02:07:46胡莉婷潘學(xué)標(biāo)馬雪晴王香茹張恒恒
    關(guān)鍵詞:新疆區(qū)域

    胡莉婷,胡 琦,潘學(xué)標(biāo),馬雪晴,徐 琳,王香茹,張恒恒

    ?

    氣候變暖和覆膜對(duì)新疆不同熟性棉花種植區(qū)劃的影響

    胡莉婷1,2,胡 琦1※,潘學(xué)標(biāo)1,馬雪晴1,徐 琳1,王香茹2,張恒恒2

    (1. 中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100193; 2. 中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所/棉花生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,安陽(yáng) 455000)

    利用新疆地區(qū)及外延200 km范圍內(nèi)173個(gè)氣象站點(diǎn)1960-2015年逐日氣溫資料和ANUSPLIN插值軟件,綜合分析了氣候變暖背景下新疆≥10 ℃的有效積溫、無霜凍期和7月平均溫度的時(shí)空變化特征,并結(jié)合地膜增溫效應(yīng)機(jī)制,探討了氣候變暖和覆膜對(duì)不同熟性棉區(qū)種植界限及可種植面積的影響。結(jié)果表明:新疆熱量資源總體呈增加的趨勢(shì),空間分布呈平原和盆地高(多)于山地。氣候變暖背景下,中熟棉區(qū)、早中熟棉區(qū)的種植面積增加,早熟棉區(qū)的可種植面積變化不明顯,特早熟棉區(qū)和不適宜區(qū)均減少。地膜增溫機(jī)制對(duì)中熟、特早熟棉區(qū)和不宜棉區(qū)的種植界限基本無影響,但對(duì)準(zhǔn)噶爾盆地區(qū)域各棉區(qū)的種植界限影響顯著,其中準(zhǔn)噶爾盆地內(nèi)早中熟棉區(qū)東擴(kuò)65 km,早熟棉區(qū)東擴(kuò)范圍在0~300 km;地膜增溫機(jī)制下早中熟Ⅱ葉塔次亞棉區(qū)可種植區(qū)域增加了5.47×104km2,早熟和特早熟棉區(qū)分別減少1.4%和1.6%,但對(duì)中熟棉區(qū)種植面積基本無影響。研究能為新疆不同熟性棉花的區(qū)劃和高產(chǎn)提供科學(xué)依據(jù)。

    氣候變暖;棉花;區(qū)劃;地膜;新疆;種植界限

    0 引 言

    熱量資源是農(nóng)業(yè)氣候資源的重要組成部分,其變化直接影響農(nóng)作物布局、種植制度和農(nóng)業(yè)生產(chǎn)[1-2]。IPCC第五次報(bào)告指出,近百年來全球海陸表面平均溫度呈線性上升趨勢(shì),且2003-2012年平均溫度較1850-1990年上升了0.78 ℃,氣候變暖現(xiàn)象毋庸置疑[3]。近年來,研究學(xué)者越來越重視有關(guān)氣候變暖對(duì)農(nóng)業(yè)生產(chǎn)的影響研究,有研究表明,氣候變暖,農(nóng)業(yè)熱量資源增加,對(duì)農(nóng)業(yè)生產(chǎn)的影響突出表現(xiàn)為作物生長(zhǎng)季延長(zhǎng)、農(nóng)作物種植界限北移[4-5]、適宜種植喜溫作物的范圍擴(kuò)大等[6-7]。

    新疆是中國(guó)優(yōu)質(zhì)棉生產(chǎn)基地,目前新疆棉花總產(chǎn)量占全國(guó)總產(chǎn)量的60%以上[8]。研究表明,隨著氣候變暖,新疆年均氣溫升高、≥10 ℃積溫增加、無霜凍期延長(zhǎng),對(duì)新疆棉花的生長(zhǎng)發(fā)育產(chǎn)生一定的影響[9-10]。受熱量資源變化的影響,全疆宜棉區(qū)面積增加,次宜棉區(qū)和不宜棉區(qū)減少[11-12];南疆中熟和早中熟棉區(qū)的面積擴(kuò)大,特早熟和不宜棉區(qū)減少[13]。在實(shí)際生產(chǎn)中,新疆棉花品種繁雜,且部分地區(qū)出現(xiàn)了品種熟性選擇不合理、盲目擴(kuò)大種植面積的問題,嚴(yán)重影響棉花品質(zhì)和產(chǎn)量[14-15]。因此,研究近幾十年來氣候變化對(duì)不同熟性棉花種植區(qū)劃的影響具有較大的積極和現(xiàn)實(shí)意義。覆膜是新疆植棉的關(guān)鍵技術(shù),至1994年,地膜植棉面積已達(dá)全疆總播種面積的96.67%[16]。有研究表明,地膜覆蓋不僅增加土壤溫度、保墑,還會(huì)對(duì)棉株生長(zhǎng)中所需的氣積溫存有補(bǔ)償?shù)淖饔肹17-20]。全球氣候變暖和地膜植棉技術(shù)都會(huì)影響棉花生長(zhǎng)發(fā)育中所需的熱量條件,進(jìn)而可能會(huì)影響棉花的種植區(qū)域。然而目前,大多研究者僅對(duì)全球變暖背景下新疆棉花適宜性分區(qū)進(jìn)行研究[11-12,21],缺乏對(duì)氣候變暖背景下不同熟性棉花種植區(qū)劃變化的研究以及定量分析地膜增溫機(jī)制對(duì)不同熟性棉花種植區(qū)域的影響。本文擬基于ANUSPLIN氣象插值軟件和80%保證率,利用新疆及其鄰近國(guó)內(nèi)外173個(gè)觀測(cè)站點(diǎn)的氣候資料,分析氣候變暖下新疆不同熟性棉花種植區(qū)劃變化特征,并結(jié)合地膜覆蓋增溫效應(yīng),進(jìn)一步研究不同熟性棉花種植界限變化特征,為新疆地區(qū)棉花品種熟性選擇和合理布局提供科學(xué)依據(jù),促進(jìn)新疆棉花高產(chǎn)穩(wěn)產(chǎn)的發(fā)展。

    1 資料來源及處理方法

    1.1 研究區(qū)域概況及數(shù)據(jù)來源

    新疆(34°15′~49°10′N,73°20′~96°25′E)地處亞歐大陸的中部,遠(yuǎn)離海洋,具有典型的大陸性氣候特點(diǎn),氣溫變化劇烈,日照充足,降水稀少。新疆地貌形態(tài)獨(dú)特復(fù)雜,基本可概括為“三山夾兩盆”,天山橫亙于新疆中部,習(xí)慣上稱天山以北為北疆,以南為南疆,哈密、吐魯番盆地為東疆(圖1)。新疆豐富的光熱資源、廣袤的土地資源以及山地降水和三大山脈的積雪冰川融水所形成的較穩(wěn)定的河流徑流灌溉資源,為棉花種植提供了得天獨(dú)厚的生態(tài)環(huán)境。

    注:由于本文研究區(qū)域僅為新疆區(qū)域,故其他地區(qū)的海拔高度沒有顯示。

    研究的氣象數(shù)據(jù)包括1960—2015年新疆維吾爾自治區(qū)內(nèi)及其外延200 km范圍的國(guó)內(nèi)外總共173個(gè)氣象站點(diǎn)逐日氣溫資料(平均溫度、最高溫度和最低溫度),其中國(guó)內(nèi)站點(diǎn)88個(gè),分布于新疆、甘肅和青海,國(guó)外站點(diǎn)85個(gè),集中在俄羅斯、蒙古、哈薩克斯坦、烏茲別克斯坦、吉爾吉斯斯坦和塔吉克斯坦等國(guó)家,數(shù)據(jù)來源于中國(guó)氣象數(shù)據(jù)網(wǎng)[22]和美國(guó)國(guó)家海洋和大氣管理局(national oceanic and atmospheric administration,NOAA)全球日值數(shù)據(jù)集 (global surface summary of the day,GSOD)[23](圖1)。國(guó)外站點(diǎn)僅為空間插值時(shí)使用,其目的是提高氣象要素空間插值的準(zhǔn)確性,降低對(duì)研究區(qū)域邊緣插值的誤差[24]。

    數(shù)字高程模型(digital elevation mode,DEM)數(shù)據(jù)來自美國(guó)地質(zhì)調(diào)查局(United States Geological Survey,USGS)[25],原始空間分辨率(約900 m)的數(shù)據(jù)經(jīng)過重采樣得到8 km×8 km的柵格數(shù)據(jù)集。

    1.2 研究方法

    1.2.1 熱量資源計(jì)算方法

    有效積溫為一段時(shí)間內(nèi)日平均氣溫與生物學(xué)最低溫度之差的總和,≥10 ℃有效積溫(AT10)是棉花生長(zhǎng)重要的熱量指標(biāo)之一。本文中≥10 ℃有效積溫由棉花生長(zhǎng)季內(nèi)日有效溫度(即日平均溫度減去10 ℃的差值)相加而得,即

    式中和分別為棉花生長(zhǎng)季內(nèi)逐日平均氣溫≥10 ℃的起始日期和終止日期,采用5日滑動(dòng)平均法確定;T為≥10 ℃的日平均氣溫,℃。

    以日最低溫度≤0 作為霜日出現(xiàn)的指標(biāo),終霜日至初霜日之間的持續(xù)日數(shù)即為無霜凍期,初終霜日采用5日滑動(dòng)平均法確定。

    分別使用氣候傾向率及檢驗(yàn)[26]研究分析新疆地區(qū)1960—2015年≥10 ℃的有效積溫、無霜凍期和7月平均溫度等熱量資源及其顯著水平的時(shí)空變化特征。

    1.2.2 保證率計(jì)算

    保證率指大于等于或小于等于某要素值出現(xiàn)的可能性或幾率。在進(jìn)行農(nóng)業(yè)氣候分析或提供農(nóng)業(yè)氣候建議時(shí),不能只利用氣候要素的多年平均值來分析問題,至少要有80%以上的保證程度,現(xiàn)常用80%保證率為指標(biāo)值。

    1.2.3 有效氣積溫補(bǔ)償值推算

    肖明等[20,27]在棉田地膜覆蓋增溫效應(yīng)的研究中指出,地膜覆蓋增溫效應(yīng)主要在開花前的苗期和蕾期,并給出了根據(jù)苗期和蕾期的氣溫推算地膜覆蓋下有效氣積溫補(bǔ)償值的計(jì)算方法:

    1)根據(jù)≥10 ℃的初日、苗期和蕾期所需的積溫值估算出新疆地區(qū)各年份露地棉苗期、蕾期天數(shù)N1和N2;

    2)根據(jù)露地棉苗期、蕾期的逐日平均氣溫及天數(shù)采用試算法估算地膜棉的苗期和蕾期天數(shù)N1和N2,計(jì)算公式如下:

    式中為地積溫對(duì)地膜棉有效氣積溫的補(bǔ)償系數(shù),單位(℃·d) /(℃·d)可省略,苗期=0.843,蕾期=0.207;NN分別為露地棉和地膜棉的某生育期天數(shù),d;TT分別為露地棉和地膜棉某生育期內(nèi)的日平均氣溫,℃。

    3)由地膜棉的苗期和蕾期天數(shù)N1和N2估算地膜棉苗期和蕾期的有效氣積溫補(bǔ)償值DATM1和DATM2,公式如下:

    式中DATM為某一生育期內(nèi)地膜棉的有效氣積溫補(bǔ)償值,℃·d。

    4)由式(4)估算苗期和蕾期有效氣積溫總補(bǔ)償值ΣDATM

    ΣDATM=DATM1+DATM2(4)

    趙黎等[28]將此計(jì)算過程添加到棉花模型GOSSYM中進(jìn)而建立新的氣象文件,并用該模型對(duì)北疆3個(gè)試點(diǎn)(石河子148團(tuán)、奎屯、瑪納斯)2個(gè)棉花品種的生育期、產(chǎn)量及主要性狀進(jìn)行模擬研究。3個(gè)試點(diǎn)地膜棉苗期天數(shù)模擬值與實(shí)測(cè)值一致;蕾期天數(shù)模擬值與實(shí)測(cè)值擬合率均在90%以上,相對(duì)誤差均在誤差允許范圍內(nèi),其中1997年石河子148團(tuán)地膜棉蕾期天數(shù)模擬值與實(shí)測(cè)值相對(duì)誤差為3.3%,1999年瑪納斯蕾期天數(shù)相對(duì)誤差為3.8%;3個(gè)地區(qū)地膜棉苗期和蕾期天數(shù)模擬值與實(shí)測(cè)值擬合程度較好。說明此計(jì)算過程能夠在不同氣候條件、不同品種熟性下較理想地估算地膜棉地有效氣積溫補(bǔ)償值。

    1.2.4 區(qū)劃指標(biāo)

    ≥10 ℃積溫是評(píng)價(jià)一個(gè)地區(qū)植棉有利程度的重要指標(biāo),其與霜前花產(chǎn)量呈正相關(guān)關(guān)系,該積溫值越高表明棉花生長(zhǎng)期越長(zhǎng),霜前花產(chǎn)量越高;無霜凍期是表明棉花利用熱量條件的限制因子;而棉花花期需要一定的熱量強(qiáng)度來滿足優(yōu)質(zhì)纖維生長(zhǎng)的熱量需求,即用7月平均溫度作為輔助指標(biāo)[29-31]。為此,綜合前人研究結(jié)果,本文選用的區(qū)劃指標(biāo)為≥10 ℃的有效積溫、無霜凍期和7月平均溫度,具體劃分標(biāo)準(zhǔn)見表1[29-31]。

    表1 新疆棉區(qū)分區(qū)氣候指標(biāo)

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

    氣候要素的空間柵格數(shù)據(jù)利用ANUSPLIN軟件插值來獲取。基于ANUSPLIN平臺(tái),以經(jīng)度和緯度為獨(dú)立自變量,以DEM數(shù)據(jù)為協(xié)變量,對(duì)新疆全區(qū)及外延200 km范圍內(nèi)的國(guó)內(nèi)外173個(gè)氣象站點(diǎn)56 a逐日氣溫資料進(jìn)行插值運(yùn)算,生成新疆地區(qū)1960-2015年逐日氣溫8 km×8 km分辨率的時(shí)間序列柵格數(shù)據(jù)集。ANUSPLIN是基于薄盤樣條理論編寫的針對(duì)氣候數(shù)據(jù)曲面擬合的函數(shù)統(tǒng)計(jì)模型,是運(yùn)用普通薄盤和局部薄盤樣條函數(shù)作為理論的插值方法,能同時(shí)進(jìn)行多個(gè)表面的空間插值,適用于時(shí)間序列的氣象數(shù)據(jù);ANUSPLIN除引入普通的樣條自變量外,還允許引入海拔等協(xié)變量[32-33]。該插值方法的總體精度水平優(yōu)于反距離加權(quán)插值法和克里格插值法,且對(duì)于中國(guó)山區(qū)氣溫的插值具有較高精度[34]。柵格數(shù)據(jù)處理均利用Matlab2010軟件實(shí)現(xiàn)。

    相關(guān)研究計(jì)算指出,1960-2015年間,1989年為中國(guó)溫帶地區(qū)氣溫突變年份[24]。故本文將1960-2015年分為2個(gè)時(shí)間段(1960-1989年、1990-2015年)研究全球氣候變暖對(duì)新疆不同熟性棉花種植區(qū)域的影響。

    2 結(jié)果與分析

    2.1 熱量資源空間分布特征

    圖2為80%保證率下≥10 ℃積溫、無霜凍期和7月平均溫度及氣候傾向率的空間分布特征。

    圖2 1960-2015年新疆地區(qū)80%保證率下農(nóng)業(yè)熱量資源指標(biāo)、氣候傾向率及其氣候傾向率顯著水平的空間分布

    1960-2015年新疆地區(qū)≥10 ℃積溫、無霜凍期和7月平均溫度的空間分布如圖2a~圖2c所示。新疆地區(qū)近56 a≥10 ℃積溫平均為2 700 ℃·d,無霜凍期平均為201.3 d,7月平均溫度為20.0 ℃。新疆各地?zé)崃抠Y源的變化與地勢(shì)密切相關(guān),有較明顯的區(qū)域性差異,全疆≥10 ℃的積溫、無霜凍期及最熱月7月平均溫度空間分布基本一致,即積溫和7月平均溫度南疆高、北疆低,平原和盆地高、山地低;無霜凍期南疆長(zhǎng)于北疆,平原和盆地長(zhǎng)于山地。塔里木盆地東部、吐魯番盆地和哈密盆地地區(qū)≥10 ℃積溫均大于4 500 ℃·d,無霜凍期均大于200 d,7月平均溫度均大于29 ℃,是新疆地區(qū)熱量資源最豐富的地區(qū);且以該3個(gè)地區(qū)為中心,積溫、無霜凍期和7月平均溫度的空間分布趨勢(shì)基本一致,即由中心向盆地周邊依次減小。準(zhǔn)噶爾盆地和新疆東北角鴨子泉地區(qū)積溫均位于3 500 ℃·d以上,無霜凍期均>180 d,7月平均溫度均>24.6 ℃,其熱量資源僅次于上述3個(gè)熱量資源最豐富的地區(qū)。由于地勢(shì)海拔較高的原因,阿爾泰山、天山和昆侖山地區(qū)熱量資源較少,積溫均低于3 200 ℃·d,無霜凍期均小于165 d,7月平均溫度均小于23 ℃。

    1960-2015年新疆地區(qū)≥10 ℃積溫、無霜凍期和7月平均溫度的氣候傾向率及其顯著水平空間分布如圖2d~圖2i所示。全疆近56 a≥10 ℃積溫氣候傾向率為64.7 ℃·d/10 a,無霜凍期氣候傾向率為3.3 d/10 a,7月平均溫度為0.2 ℃/10 a。結(jié)合表2得,近56 a來新疆積溫基本呈現(xiàn)顯著升高趨勢(shì),積溫顯著升高的區(qū)域占全疆面積的93.6%,新疆南部、吐魯番盆地、哈密盆地及準(zhǔn)噶爾盆地地區(qū)升高幅度最為明顯,積溫氣候傾向率均大于90.0 ℃·d/10 a,其中哈密盆地及新疆東南地區(qū)庫(kù)姆塔格沙漠的增速最大可達(dá)到195.0 ℃·d/10 a;新疆山區(qū)積溫氣候傾向率基本為負(fù)值,56 a來積溫顯著降低。無霜凍期呈顯著增加趨勢(shì)的面積占全疆84.6%,表現(xiàn)為新疆南部地區(qū)增速明顯,最高可達(dá)15.2 d/10 a;山區(qū)增速均大于3.0 d/10 a,且無霜凍期氣候傾向率高于盆地及平原地區(qū)。7月平均溫度呈顯著增高趨勢(shì)的面積占69.8%,其中0.03%的格點(diǎn)的7月平均溫度顯著降低,集中在塔里木盆地中的和田河流域;新疆南部和東部地區(qū)的增速最大,最大可達(dá)0.9 ℃/10a;56 a來準(zhǔn)噶爾盆地和和塔里木盆地地區(qū)的7月平均溫度基本無變化。3個(gè)農(nóng)業(yè)資源熱量指標(biāo)在塔里木盆地中西部均無顯著變化。

    表2 1960-2015年新疆地區(qū)≥10 ℃積溫、無霜凍期和7月平均溫度氣候傾向率不同顯著水平面積占比

    2.2 氣候變暖背景下不同熟性棉花種植區(qū)劃

    2.2.1 氣候變暖背景下棉花區(qū)劃變化特征

    根據(jù)以上區(qū)劃指標(biāo)和方法,分別得出80%保證率下全疆1960-1989年、1990-2015年2個(gè)時(shí)段不同熟性棉花種植區(qū)劃結(jié)果(圖3)。對(duì)比分析了1960-1989年(圖3a)和1990-2015年(圖3b)80%保證率下新疆地區(qū)不同熟性棉花種植區(qū)域的變化。氣候變暖背景下中熟棉區(qū)的種植面積明顯增加;早中熟棉區(qū)的可種植區(qū)域呈增大趨勢(shì);早熟棉區(qū)的可種植面積變化趨勢(shì)不明顯;特早熟棉區(qū)和不適宜區(qū)的可種植區(qū)域均減小。

    注:1989年為中國(guó)溫帶地區(qū)氣溫突變年份[24]。

    由圖3和表3可得:

    1)中熟棉區(qū):中熟棉區(qū)是全疆最小的棉區(qū),1960-1989年該區(qū)僅分布于吐魯番盆地和哈密盆地的部分地區(qū),可種植面積僅1.49×104km2;1990-2015年時(shí)中熟棉區(qū)不僅在吐魯番盆地、哈密盆地的分布面積明顯擴(kuò)大,而且還分布于塔里木盆地東部羅布泊地區(qū),其可種植總面積為5.31×104km2;與1960-1989年相比,中熟棉區(qū)面積增加3.82×104km2,占全疆總面積的比率增加2.3%。

    2)早中熟棉區(qū):早中熟棉區(qū)分布較廣,是新疆的最大棉區(qū),由于該區(qū)經(jīng)緯度跨度大,生態(tài)環(huán)境差異顯著,故又分為2個(gè)次亞區(qū):葉塔次亞區(qū)和塔哈次亞區(qū),葉塔次亞區(qū)分布面積顯著大于塔哈次亞區(qū)。在1960-1989年內(nèi)葉塔次亞棉區(qū)可種植面積為44.82×104km2,主要分布于吐魯番盆地、哈密盆地以及塔里木盆地的大部分地區(qū),此外還零星分布于新疆東北部海拔500 m以下的地區(qū);1990-2015年內(nèi)葉塔次亞區(qū)增至47.48×104km2,與1960-1989年相比,該棉區(qū)除在吐魯番盆地、哈密盆地和新疆東北部的分布面積明顯擴(kuò)大外,另在準(zhǔn)噶爾盆地西部有大面積出現(xiàn),其占地面積增加2.66×104km2;此外在兩時(shí)段內(nèi),該棉區(qū)位于塔里木盆地的分布區(qū)域存在不同,與1960-1989年相比,該棉區(qū)在塔里木盆地的中西部和東部的分布減少,但在南部的分布明顯增加。1960-1989年內(nèi),塔哈次亞區(qū)分布于準(zhǔn)噶爾盆地的中西部、吐魯番盆地和哈密盆地邊緣、塔里木盆地的北部、東部邊緣及盆地南部和新疆東北部海拔500m以上的地區(qū),其可種植面積14.61×104km2,1990-2015年內(nèi)該棉區(qū)可種植面積為14.77×104km2,雖在準(zhǔn)噶爾盆地、吐魯番盆地、哈密盆地及塔里木盆地南部的面積明顯減少,但在塔里木盆地中東部、北部以及新疆東北部海拔500 m以上的地區(qū)的面積增加,故兩時(shí)段內(nèi)可種植面積變化較小。

    3)早熟棉區(qū):該棉區(qū)主要分布于北疆和東疆。北疆片位于準(zhǔn)噶爾盆地、塔城區(qū)附近以及新疆東北部海拔500~1 000 m區(qū)域內(nèi),東疆片位于吐魯番盆地和哈密盆地邊緣500~1 000 m區(qū)域內(nèi),此外還零星分布在天山南坡的和靜縣附近及塔里木盆地邊緣。兩段時(shí)間內(nèi),該棉區(qū)可種植面積基本穩(wěn)定在6.97×104km2;在南疆的分布區(qū)域基本無變化,北疆片在準(zhǔn)噶爾盆地南部分布的面積減少,但在準(zhǔn)噶爾盆地的西部及新疆東北部海拔500~1 000 m區(qū)域內(nèi)分布的面積增加。

    4)特早熟棉區(qū):該棉區(qū)主要分布于準(zhǔn)噶爾盆地、吐魯番盆地、哈密盆地邊緣以及新疆東北部的海拔1 000~1 500 m區(qū)域,此外還零星分布于天山南坡和靜縣附近海拔1 000~1 500 m區(qū)域以及塔里木盆地邊緣。1960-1989年時(shí)該棉區(qū)可種植面積為8.13×104km2,而1990-2015年時(shí)面積為7.64×104km2,減少了0.49×104km2。

    5)不適宜棉區(qū):阿爾泰山、天山及昆侖山脈周邊由于地勢(shì)海拔等因素導(dǎo)致農(nóng)業(yè)熱量資源不足,不適宜種植棉花。兩時(shí)段內(nèi),不適宜種植區(qū)面積減少6.14×104km2,棉花適宜種植區(qū)域增加,其主要是由于氣候變暖背景下北疆西部伊犁河流域和準(zhǔn)噶爾盆地東部海拔1 000 m以上區(qū)域熱量條件已滿足棉花生長(zhǎng)的需求,故該區(qū)域棉花可種植面積增加。

    表3 1960-2015年新疆各棉區(qū)可種植區(qū)域變化

    2.2.2 不同熟性棉花區(qū)劃的熱量限制因素

    基于區(qū)劃指標(biāo)以及全疆各地區(qū)的平均氣溫,確定全疆各地區(qū)棉花種植區(qū)域的熱量資源限制因素,圖4表示新疆不同熟性棉區(qū)熱量限制因素的空間分布特征。由圖4可知,山區(qū)附近各地?zé)崃抠Y源限定因素較為一致,而盆地和平原等海拔較低的地區(qū)不同熟性棉花種植的限制因素較為復(fù)雜。阿爾泰山、天山和昆侖山脈處由于地勢(shì)海拔較高,氣溫較低,受≥10 ℃積溫、無霜凍期和7月平均溫度3因素共同限制而不適宜種植棉花;除山脈外,三大山區(qū)附近不宜種植棉花的主要限制因素為≥10 ℃積溫和7月平均溫度,且兩熱量因素也是塔里木盆地處絕大部分的早中熟葉塔次亞棉區(qū)的熱量限制因素;全疆受≥10 ℃積溫和無霜凍期兩因素共同限制的地區(qū)較少,主要零星分布于天山和昆侖山區(qū)的不適宜棉區(qū)內(nèi);塔里木盆地處早中熟塔哈次亞棉區(qū)以及早中熟葉塔次亞棉區(qū)西部、中部和東部的部分區(qū)域主要受7月平均溫度的限制,此外特早熟棉區(qū)的絕大多數(shù)區(qū)域也受此熱量因素的限制;準(zhǔn)噶爾盆地、塔城區(qū)附近、新疆東北部、吐魯番盆地、哈密盆地及塔里木盆地區(qū)域內(nèi)的早中熟塔哈次亞棉區(qū)、早熟棉區(qū)和部分特早熟棉區(qū)主要受≥10 ℃積溫限制。

    2.3 地膜增溫機(jī)制下棉花區(qū)劃變化

    地膜植棉技術(shù)于20世紀(jì)90年代初基本應(yīng)用于整個(gè)新疆棉區(qū),本文基于棉區(qū)區(qū)劃指標(biāo)、地膜覆蓋對(duì)有效氣溫的補(bǔ)償值[20,27,35],定性分析1990—2015年時(shí)80%保證率下新疆地區(qū)不同熟性棉花在考慮地膜增溫效應(yīng)前后,其種植界限和可種植區(qū)域的變化特征。

    由圖3b、圖5a和表4可知,地膜覆蓋下氣積溫的補(bǔ)償效應(yīng)機(jī)制對(duì)不同熟性棉區(qū)的可種植區(qū)域面積和種植界限均存在一定影響,且對(duì)其種植界限影響顯著??紤]地膜增溫效應(yīng)前后,中熟棉區(qū)的可種植面積基本穩(wěn)定在5.31×104km2,其種植界限及區(qū)域無變化;早中熟葉塔次亞棉區(qū)的可種植區(qū)域有明顯的增大趨勢(shì);早中熟塔哈次亞棉區(qū)、早熟棉區(qū)、特早熟棉區(qū)的可種植區(qū)域和不適宜區(qū)面積均減小,其中不適宜區(qū)變化不明顯,可種植區(qū)域占全疆棉區(qū)的比率僅減少0.1%。

    圖4 1960—2015年新疆地區(qū)棉區(qū)熱量限制因素的空間分布

    圖5 1990-2015年地膜覆蓋下新疆棉花可種植區(qū)域及種植界限及可種植區(qū)域變化特征

    表4 1990-2015年地膜覆蓋下新疆各棉區(qū)可種植區(qū)域變化

    早中熟棉區(qū):早中熟兩次亞棉區(qū)在考慮地膜增溫效應(yīng)前后,其區(qū)域界限和可種植區(qū)域均發(fā)生一定的變化。從圖5b和表4可知,考慮地膜增溫效應(yīng)前后早中熟葉塔次亞棉區(qū)的四大分布區(qū)(準(zhǔn)噶爾盆地、吐魯番盆地及哈密盆地、塔里木盆地和新疆東北部)的區(qū)域界限均有向四周擴(kuò)的變化趨勢(shì),其中準(zhǔn)噶爾盆地區(qū)域東擴(kuò)趨勢(shì)最大,準(zhǔn)噶爾盆地區(qū)域南擴(kuò)趨勢(shì)次之,平均擴(kuò)大65和30 km;哈密盆地區(qū)域東擴(kuò)趨勢(shì)明顯,變化范圍在10~90 km;新疆東北部的北擴(kuò)范圍在10~40 km;塔里木盆地區(qū)域南擴(kuò)最大為20 km;考慮地膜增溫效應(yīng)后,葉塔次亞棉區(qū)的可種植區(qū)域面積由47.48×104km2增至52.95×104km2,增加了3.3%。從圖5c和表4中可知,考慮地膜增溫效應(yīng)后,塔哈次亞棉區(qū)在準(zhǔn)噶爾盆地的區(qū)域界限向東擴(kuò)的趨勢(shì)較為明顯,平均擴(kuò)大70 km;雖在新疆西北部的塔城區(qū)附近原早熟棉區(qū)現(xiàn)可種植塔哈次亞棉區(qū),且在新疆東北部地區(qū)處塔哈次亞棉區(qū)的可種植界限稍微向西北方向擴(kuò)大,但吐魯番、哈密地區(qū)南部處(41°N附近)原塔哈次亞棉區(qū)的部分區(qū)域現(xiàn)可種植葉塔次亞棉區(qū),致使塔哈次亞棉區(qū)可種植區(qū)域略減小,可種植面積由14.77×104km2縮小至14.45×104km2。

    早熟棉區(qū):由圖5d和表4可知,該棉區(qū)在準(zhǔn)噶爾盆地區(qū)域的種植界限有明顯的東擴(kuò)趨勢(shì),變化范圍在0~300 km,但該區(qū)域內(nèi)東擴(kuò)大的面積較?。辉撁迏^(qū)原可種植區(qū)域(塔城區(qū)附近、新疆東北部海拔500~1 000 m區(qū)域內(nèi)、哈密盆地和塔里木盆地邊緣)的部分區(qū)域現(xiàn)可用來種植早中熟塔哈次亞棉區(qū),且該些區(qū)域原種植界限無明顯平移外擴(kuò)趨勢(shì),致使該棉區(qū)可種植區(qū)域面積明顯減少,減少2.32×104km2,減少了1.4%。

    特早熟棉區(qū):由圖3b、圖5a和表4可知,考慮地膜增溫效應(yīng)后,特早熟棉區(qū)的種植界限基本無變化,而原種植區(qū)域的較大部分現(xiàn)可以種植早熟棉花,因而特早熟棉區(qū)的可種植面積減小,可種植面積由7.64×104km2減少至4.98×104km2,減少了1.6%。

    3 討 論

    前人研究表明,全球氣候變暖背景下,新疆北疆和南疆熱量資源總體各呈現(xiàn)“平原和盆地多,山區(qū)少”的空間分布格局[11,13,36];絕大部分地區(qū)10 ℃初(終)日提早(推遲),終(初)霜凍日提前(推遲),≥10 ℃積溫持續(xù)天數(shù)和無霜凍期均延長(zhǎng),積溫增加,棉花適宜種植區(qū)擴(kuò)大[11,16,36-37],本文研究結(jié)果與前人研究結(jié)論類似。

    文中研究表明氣候變暖,熱量資源得到改善,使新疆地區(qū)中熟棉區(qū)和早中熟棉區(qū)的可種植區(qū)域均不同程度的增加,不適宜種植棉花的區(qū)域面積減小,此與李景林等[11-12]認(rèn)為北疆宜棉區(qū)大幅擴(kuò)張,次宜棉區(qū)、風(fēng)險(xiǎn)棉區(qū)和不宜棉區(qū)不同程度的減小的結(jié)果類似;本文研究得出南疆早中熟棉區(qū)可種植面積增大、早熟棉區(qū)可種植區(qū)域變化不明顯及特早熟棉區(qū)減少等特征,與張山清等[13]得到的南疆棉花區(qū)劃變化較為接近。本文結(jié)果表明,與北疆相比,南疆宜棉區(qū)中不同熟性棉區(qū)的主要熱量資源限制因子除了受≥10 ℃積溫限制外,還受7月平均溫度的限制,這表明在南疆各地制定和規(guī)劃棉花區(qū)劃時(shí),除要考慮積溫外,還必須高度重視7月平均溫度這一指標(biāo)。

    地膜增溫效應(yīng)機(jī)制下早熟和特早熟棉區(qū)可種植區(qū)域減少,中熟棉區(qū)和早中熟棉區(qū)的種植界限變化趨勢(shì)主要為向四周外擴(kuò),其可種植區(qū)域均不同程度的增加,這與本文及前人研究得出的氣候變暖背景下(1990-2015年與1960-1989年比較)棉區(qū)的種植界限和可種植面積的變化趨勢(shì)類似[11-13,21],但其主要影響熱量需求相對(duì)較低的特早熟棉區(qū)和早熟棉區(qū),而對(duì)熱量需求最高的中熟棉區(qū)的種植界限及可種植區(qū)域均無影響,原因可能為地膜增溫效應(yīng)僅作用于棉花生育前期的苗期和蕾期,對(duì)有效氣積溫的補(bǔ)償值有一定的限度[20,27],不足以使其他棉區(qū)的熱量資源增加至滿足中熟棉花生長(zhǎng)發(fā)育的需求。

    本文中棉田地膜增溫公式全部引用肖明等[20,27]的研究結(jié)果,該公式采用試算法由露地棉苗期、蕾期的逐日平均氣溫及天數(shù)來估算地膜棉的苗期和蕾期天數(shù),進(jìn)而求得地膜棉的有效氣積溫補(bǔ)償值。該計(jì)算過程以土壤溫度與氣溫的相關(guān)關(guān)系為核心,考慮了氣溫變化對(duì)地溫的影響,且在不同年份和不同試點(diǎn)均得到較好的驗(yàn)證,可操作性較強(qiáng)[28]。但未來仍需用試驗(yàn)數(shù)據(jù)對(duì)該公式進(jìn)行驗(yàn)證已確保其時(shí)效性和準(zhǔn)確性。

    本文關(guān)于氣候變化和覆膜下新疆棉花可種植區(qū)域變化特征的研究中,對(duì)不同熟性棉花的區(qū)劃為理論計(jì)算的氣候區(qū)劃。在區(qū)劃指標(biāo)中僅考慮了熱量資源因素,而一個(gè)地區(qū)實(shí)際的種植制度不僅取決于熱量資源,同時(shí)受水資源、土壤條件、當(dāng)?shù)氐慕?jīng)濟(jì)水平、社會(huì)效益等因素綜合影響,故未來可綜合考慮水資源、土壤條件、社會(huì)經(jīng)濟(jì)等因素,對(duì)棉花種植區(qū)劃進(jìn)行更深入地研究,為新疆棉花區(qū)劃提供科學(xué)依據(jù),實(shí)現(xiàn)棉花高產(chǎn)穩(wěn)產(chǎn)。

    4 結(jié) 論

    本文基于柵格計(jì)算了1960-2015年新疆地區(qū)的≥10 ℃積溫、無霜凍期和7月平均溫度,研究了氣候變化背景下新疆地區(qū)農(nóng)業(yè)熱量資源的時(shí)空變化特征,比較和探討了氣候變暖和覆膜影響下,不同熟性棉花可種植區(qū)的變化特征以及棉花區(qū)劃的熱量資源限制因素,得出以下主要結(jié)論:

    1)1960-2015年新疆地區(qū)總體表現(xiàn)為≥10 ℃積溫、無霜凍期和7月平均溫度分別以64.7 ℃·d/10a、3.3 d/10 a和0.2 ℃/10 a的傾向率呈顯著增加趨勢(shì);全疆熱量資源有較為明顯的區(qū)域性差異,≥10 ℃積溫、無霜凍期和7月平均溫度的變化與地勢(shì)海拔密切相關(guān),總體呈現(xiàn)南疆高(多)于北疆,平原和盆地高(多)于山地。

    2)與1960-1989年相比1990-2015年時(shí)不同熟性棉區(qū)的可種植區(qū)存在不同程度的變化趨勢(shì),中熟棉區(qū)的種植面積明顯增加3.82×104km2,主要增加區(qū)域位于塔里木盆地東部地區(qū);早中熟棉區(qū)的可種植區(qū)域有增大趨勢(shì),早中熟葉塔次亞棉區(qū)可種植區(qū)的增加區(qū)域主要分布于準(zhǔn)噶爾盆地東部、吐魯番盆地南部、哈密盆地南部以及塔里木盆地西部;早熟棉區(qū)的可種植面積變化趨勢(shì)不明顯,特早熟棉區(qū)的可種植面積和不適宜區(qū)的面積均減少。

    3)地膜增溫機(jī)制對(duì)中熟、特早熟棉區(qū)和不宜棉區(qū)的種植界限基本無影響,但對(duì)準(zhǔn)噶爾盆地區(qū)域棉區(qū)的種植界限影響顯著,準(zhǔn)噶爾盆地區(qū)域內(nèi)早中熟葉塔次亞棉區(qū)、塔哈次亞棉區(qū)平均分別東擴(kuò)65和70 km,早熟棉區(qū)東擴(kuò)范圍在0~300 km;地膜增溫機(jī)制下對(duì)中熟和不宜棉區(qū)的種植面積基本無影響,但早中熟葉塔次亞棉區(qū)可種植面積增加5.47×104km2,增加3.3%,特早熟棉區(qū)面積減小1.6%。

    4)全疆棉花種植的熱量資源限制因素的空間分布特征為,不宜棉區(qū)大部分受≥10 ℃積溫和7月平均溫度限制,而三大山脈處不宜棉區(qū)還受無霜凍期限制;北疆宜棉區(qū)中不同熟性棉區(qū)的主要熱量資源限制因子為≥10 ℃積溫,南疆主要為≥10 ℃積溫和7月平均溫度兩者共同限制或7月平均溫度單個(gè)因素限制,大部分早中熟塔哈次亞棉區(qū)、早熟棉區(qū)和部分特早熟棉區(qū)受≥10 ℃積溫限制。

    [1] 范文波,江煜,吳普特,等. 新疆石河子墾區(qū)50年氣候變化對(duì)棉花種植的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2011,29(6):244-248.

    Fan Wenbo, Jiang Yu, Wu Pute, et al. Impacts of climate change on planting proportion of cotton in Shihezi area in the recent 50 years[J].Agricultural Research in the Arid Areas, 2011, 29(6): 244-248. (in Chinese with English abstract)

    [2] Chi Zaixiang, Mo Jianguo, Kang Xueliang, et al. Meticulous regionalization of climate suitability about spring potato planting in western Guizhou based on GIS[J]. Meteorological and Environmental Research, 2012, 3(11): 65-70.

    [3] IPCC. Climate Change 2013: The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel On Climate Change[M]. Cambridge: Cambridge University Press, 2014.

    [4] 劉德祥,董安祥,梁東升,等. 氣候變暖對(duì)西北干旱區(qū)農(nóng)作物種植結(jié)構(gòu)的影響[J]. 中國(guó)沙漠,2007,27(5):831-836.

    Liu Dexiang, Dong Anxiang, Liang Dongsheng, et al. Affect of climate warming on crops planting structure in arid zone of northwestern China[J]. Journal of Desert Research, 2007, 27(5): 831-836. (in Chinese with English abstract)

    [5] 李克南,楊曉光,慕臣英,等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅷ:氣候變化對(duì)中國(guó)冬小麥冬春性品種種植界限的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué),2013,46(8):1583-1594.

    Li Kenan, Yang Xiaoguang, Mu Chenying, et al. The possible effects of global warming on cropping systems in China Ⅷ: The effects of climate change on planting boundaries of different winte-spring varieties of winter wheat in China[J]. Scientia Agricultura Sinica, 2013, 46(8): 1583-1594. (in Chinese with English abstract)

    [6] 趙俊芳,郭建平,馬玉平,等. 氣候變化背景下我國(guó)農(nóng)業(yè)熱量資源的變化趨勢(shì)及適應(yīng)對(duì)策[J]. 應(yīng)用生態(tài)學(xué)報(bào),2010,21(11):2922-2930.

    Zhao Junfang, Guo Jianping, Ma Yuping, et al. Change trends of China agricultural thermal resources under climate change and related adaptation countermeasures[J]. Chinese Journal of Applied Ecology, 2010, 21(11): 2922-2930. (in Chinese with English abstract)

    [7] 王潤(rùn)元,張強(qiáng),劉宏誼,等. 氣候變暖對(duì)河西走廊棉花生長(zhǎng)的影響[J]. 氣候變化研究進(jìn)展,2006,2(1):4-6.

    Wang Runyuan, Zhang Qiang, Liu Hongyi, et al. Impact of climate warming on cotton growth in the Hexi corridor area[J].Climate Change Research, 2006, 2(1): 4-6. (in Chinese with English abstract)

    [8] 黃璐,宋玉蘭. 新疆棉花生產(chǎn)效率發(fā)展現(xiàn)狀分析[J]. 山西農(nóng)業(yè)科學(xué),2017,45(6):1020-1023.

    Huang Lu, Song Yulan. Analysis on development status of cotton production efficiency in Xinjiang[J]. Journal of Shanxi Agricultural Sciences, 2017, 45(6): 1020-1023. (in Chinese with English abstract)

    [9] 曹占洲,毛煒嶧,李迎春,等. 近49年新疆棉區(qū)≥10 ℃終日和初霜期的變化及對(duì)棉花生長(zhǎng)的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào),2011,27(8):355-361.

    Cao Zhanzhou, Mao Weiyi, Li Yingchun, et al. The change of ≥10 ℃ terminalday and the first frost date in Xinjiang cotton-growing grea and affect of cotton growth in recent 49 years[J]. Chinese Agricultural Science Bulletin, 2011, 27(8): 355-361. (in Chinese with English abstract)

    [10] 只娟,張山清,王榮曉,等. 近40年天山北坡經(jīng)濟(jì)帶熱量資源時(shí)空變化特征研究[J]. 干旱地區(qū)農(nóng)業(yè)研究,2015,33(3):254-259.

    Zhi Juan, Zhang Shanqing, Wang Rongxiao, et al. Research on spatial-temporal variation of heat resources in economic zone of northern slope of Tianshan Mountain for recent 40 years[J]. Agricultural Research in the Arid Areas, 2015, 33(3): 254-259. (in Chinese with English abstract)

    [11] 李景林,普宗朝,張山清,等. 近52年北疆氣候變化對(duì)棉花種植氣候適宜性分區(qū)的影響[J]. 棉花學(xué)報(bào),2015,27(1):22-30.

    Li Jinglin, Pu Zongchao, Zhang Shanqing, et al. Impact of climate change on cotton-planting zoning in North Xinjiang during 1961-2012[J]. Cotton Science, 2015, 27(1): 22-30. (in Chinese with English abstract)

    [12] 普宗朝,張山清,賓建華,等. 氣候變暖對(duì)新疆烏昌地區(qū)棉花種植區(qū)劃的影響[J]. 氣候變化研究進(jìn)展,2012,8(4):257-264.

    Pu Zongchao, Zhang Shanqing, Bin Jianhua, et al. Impact of global warming on cotton-planting zoning in the Urumqi-Changji region of Xinjiang[J]. Climate Change Research, 2012, 8(4): 257-264. (in Chinese with English abstract)

    [13] 張山清,普宗朝,李景林,等. 氣候變暖背景下南疆棉花種植區(qū)劃的變化[J]. 中國(guó)農(nóng)業(yè)氣象,2015,36(5):594-601.

    Zhang Shanqing, Pu Zongchao, Li Jinglin, et al. Cotton-planting zoning change in southern Xinjiang under the background of global warming[J]. Chinese Journal ofAgrometeorology, 2015, 36(5): 594-601. (in Chinese with English abstract)

    [14] 姚源松. 新疆棉花區(qū)劃新論[J]. 中國(guó)棉花,2001,28(2):2-5.

    Yao Yuansong. New discussion on cotton-planting areas division in Xinjiang[J]. China Cotton, 2001, 28(2): 2-5. (in Chinese with English abstract)

    [15] 黃滋康,崔讀昌. 中國(guó)棉花生態(tài)區(qū)劃[J]. 棉花學(xué)報(bào),2002,14(3):185-190.

    Huang Zikang, Cui Duchang. Ecological regionalization of cotton production in China[J]. Cotton Science, 2002, 14(3): 185-190. (in Chinese with English abstract)

    [16] 戴良佐. 北疆植棉和地膜棉的起始及推廣[J]. 新疆地方志,1995(4):42-43.

    Dai Liangzuo. The beginning and extension of cotton planting and film cotton in Northern Xinjiang[J]. Local Chronicles of Xinjiang, 1995(4): 42-43. (in Chinese with English abstract)

    [17] 胡明芳,田長(zhǎng)彥. 新疆棉田地膜覆蓋耕層土壤溫度效應(yīng)研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2003,11(3):128-130.

    Hu Mingfang, Tian Changyan. The temperature effects of membrane coverage on cultivated soil cotton field in Xinjiang[J]. Chinese Journal of Eco-Agriculture, 2003, 11(3): 128-130. (in Chinese with English abstract)

    [18] 畢繼業(yè),王秀芬,朱道林. 地膜覆蓋對(duì)農(nóng)作物產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2008,24(11):172-175.

    Bi Jiye, Wang Xiufen, Zhu Daolin. Effect of plastic-film mulch on crop yield[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(11): 172-175. (in Chinese with English abstract)

    [19] 祖米來提·吐爾干,林濤,嚴(yán)昌榮,等. 地膜覆蓋時(shí)間對(duì)新疆棉田水熱及棉花耗水和產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(11):113-120.

    Zumilaiti Tuergan, Lin Tao, Yan Changrong, et al. Effect of plastic film mulching duration on soil temperature and moisture in field and cotton water comsumption and yield in Xinjiang[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(11): 113-120. (in Chinese with English abstract)

    [20] 肖明,鐘俊平,趙黎. 棉田土壤溫度與氣溫的關(guān)系及膜地增溫效應(yīng)對(duì)有效氣積溫的補(bǔ)償作用的研究[J]. 新疆農(nóng)業(yè)大學(xué)學(xué)報(bào),1998(4):5-9.

    Xiao Ming, Zhong Jungping, Zhao Li. Studies on the correlation between cotton soil temperature and air temperature, and the compensating action of soil temperature increase in the mulched cotton field on effective air accumulation temperature[J]. Journal of Xinjiang Agricultural University, 1998(4): 5-9. (in Chinese with English abstract)

    [21] 李迎春,謝國(guó)輝,李新建,等. 基于GIS的新疆阿克蘇地區(qū)棉花區(qū)劃[J]. 中國(guó)農(nóng)業(yè)氣象,2007,28(4):450-452.

    Li Yingchun, Xie Guohui, Li Xinjian, et al. Regionalization of cotton cultivation in Akesu region of Xinjiang based on GIS[J]. Chinese Journal ofAgrometeorology, 2007, 28(4): 450-452. (in Chinese with English abstract)

    [22] 中國(guó)氣象局. 氣象數(shù)據(jù)[EB/OL].2016-08-30 [2017-06-02]. http://data.cma.cn/.

    [23] National Oceanic and Atmospheric Administration(NOAA). Global summary of the day [EB/OL]. 2017-01-01 [2017-06-02]. https://www7.ncdc.noaa.gov/CDO/cdoselect.cmd/.

    [24] 張煦庭,潘學(xué)標(biāo),徐琳,等. 中國(guó)溫帶地區(qū)不同界限溫度下農(nóng)業(yè)熱量資源的時(shí)空演變[J]. 資源科學(xué),2017,39(11):2104-2115.

    Zhang Xuting, Pan Xuebiao, Xu Lin, et al. Spatio- temporal variation of agricultural thermal resources at different critical temperatures in China’s temperate zone[J]. Resources Science, 2017, 39(11): 2104-2115. (in Chinese with English abstract)

    [25] Danielson J J, Gesch D B. Global multi-resolution terrain elevation data 2010 (Gmted2010)[EB/OL].2010-09-01 [2017-06-02]. https://topotools.cr.usgs.gov/gmted_viewer/.

    [26] 魏鳳英. 現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M]. 北京:氣象出版社,2007.

    [27] 肖明,鐘俊平,李保成,等. 棉田地膜覆蓋增溫效應(yīng)和規(guī)律及對(duì)GOSSYM模型的修改[J]. 新疆農(nóng)業(yè)大學(xué)學(xué)報(bào),1997,20(4):32-39.

    Xiao Ming, Zhong Junping, Li Baocheng, et al. The effect and pattern of soil temperature increase in cotton field mulched with plastic film and the modification on the GOSSYM simulation model[J]. Journal of Xinjiang Agricultural University, 1997, 20(4): 32-39. (in Chinese with English abstract)

    [28] 趙黎,鐘俊平,趙富強(qiáng),等. 棉花生長(zhǎng)模擬模型GOSSYM及其在北疆棉區(qū)的應(yīng)用[J]. 新疆農(nóng)業(yè)大學(xué)學(xué)報(bào),1999,22(1):43-50.

    Zhao Li, Zhong Junping, Zhao Fuqiang, et al. GOSSYM and applications under the condition of north Xinjiang cotton region[J]. Journal of Xinjiang Agricultural University, 1999, 22(1): 43-50. (in Chinese with English abstract)

    [29] 徐德源. 從農(nóng)業(yè)氣候看新疆細(xì)絨棉的發(fā)展地區(qū)[J]. 新疆農(nóng)業(yè)科學(xué),1963(3):92-97.

    [30] 徐德源. 新疆棉花生態(tài)氣候區(qū)劃[J]. 新疆農(nóng)業(yè)科學(xué),1987(6):4-7.

    [31] 中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所. 中國(guó)棉花品種志:1978-2007[M]. 北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,2009.

    [32] Hutchinson M F, Xu T B. Anusplin Version 4. 2 User Guide[M]. Canberra:The Australian National University, 2004.

    [33] 姜曉劍,湯亮,劉小軍,等. 中國(guó)主要稻作區(qū)水稻生產(chǎn)氣候資源的時(shí)空特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(7):238-245.

    Jiang Xiaojian, Tang Liang, Liu Xiaojun, et al. Spatial and temporal characteristics of rice production climatic resources in main growing regions of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 238-245. (in Chinese with English abstract)

    [34] 譚劍波,李愛農(nóng),雷光斌. 青藏高原東南緣氣象要素Anusplin和Cokriging空間插值對(duì)比分析[J]. 高原氣象,2016,35(4):875-886.

    Tan Jianbo, Li Ainong, Lei Guangbin. Contrast on ANUSPLIN and Cokriging meteorological spatial interpolation in southeastern margin of Qinghai-Xizang Plateau[J]. Plateau Meteorology, 2016, 35(4): 875-886. (in Chinese with English abstract)

    [35] 劉德章,馮利平. 地膜覆蓋棉田土壤溫度的變化動(dòng)態(tài)及其模擬研究[J]. 棉花學(xué)報(bào),1994,6(2):99-102.

    Liu Dezhang, Feng Liping. Studies of variation trends of soil temperatures and its models on film-mulching cotton[J]. Cotton Science, 1994, 6(2): 99-102. (in Chinese with English abstract)

    [36] 李迎春,謝國(guó)輝,王潤(rùn)元,等. 北疆棉區(qū)棉花生長(zhǎng)期氣候變化特征及其對(duì)棉花發(fā)育的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2011,29(2):253-258.

    Li Yingchun, Xie Guohui, Wang Runyuan, et al. Climatic change during Apr to Oct in recent 47 years and its effects on growing period of cotton in Beijiang cotton planting region[J]. Agricultural Research in the Arid Areas, 2011, 29(2): 253-258. (in Chinese with English abstract)

    [37] 張山清,普宗朝,李景林,等. 氣候變暖背景下新疆無霜凍期時(shí)空變化分析[J]. 資源科學(xué),2013,35(9):182-190.

    Zhang Shanqing, Pu Zongchao, Li Jinglin, et al. The impact of global warming on frost-free periods from 1961 to 2010 in Xinjiang[J]. Resources Science, 2013, 35(9): 182-190. (in Chinese with English abstract)

    Effects of global warming and plastic mulching on cotton-planting zoning with different matures in Xinjiang

    Hu Liting1,2, Hu Qi1※, Pan Xuebiao1, Ma Xueqing1, Xu Lin1, Wang Xiangru2, Zhang Hengheng2

    (1.100193; 2.455000)

    Xinjiang is a main producing area of cotton in China. This study investigated the effects of global warming and plastic mulching on cotton-planting zoning of different matures in Xinjiang. The climatic data were from 173 meteorological stations during 1960-2015 in temperate zone around Xinjiang within a range of 200 km. The digital elevation mode data were from the United States Geological Survey, which had the resolution about 900 m. We analyzed the spatio-temporal variation of accumulated temperature not less than 10 ℃, frost-free period and mean temperature in July with the aid of ANUSPLIN interpolation software. The cotton-growing area was classified based on these 3 parameters. In addition, the change of planting boundaries and area of different cotton matures were studied with climate change and the compensation effect of plastic mulching. The main results included: 1) the accumulated temperature not less than 10 ℃, frost-free period and mean temperature in July showed an increasing trend in 1960-2015. Regional differences on thermal resources were obvious in Xinjiang. The accumulated temperature not less than 10 ℃, frost-free period and mean temperature in July were closely related to the terrain elevation. The thermal resources wasricher in the southern Xinjiang than those the northern Xinjiang and richer in plain areas than in mountain areas. 2) under global warming, the ratio of planting areas of middle mature and early-middle mature cotton all increased.The area of middle mature cotton region increased by 3.82×104km2, and the main increase areas were located in the eastern part of tarim basin. The increased area of early-middle mature cotton was widely distributed in the eastern Junggar Basin, the southern Turpan Basin, the southern Hami Basin and the western Tarim Basin. But the proportion of planting area of early maturing cotton did not change significantly and exceptional early-mature and unsuitable areas decreased. 3) the accumulated temperature not less than 10oC and mean temperature in July were the main limiting factors in unsuitable planting areas, besides the areas near Tianshan, Altai and Kunlun Mountains were also restricted by frost-free period. The main thermal limiting factors in early-middle mature cotton areas (Taha), early mature cotton areas and most regions of exceptional early-mature were limited by the accumulated temperature not less than 10 ℃. 4) cotton plastic mulching had no significant influence on the planting boundaries of middle mature, early mature and unsuitable planting areas. But it had a significant influence on the planting boundaries of cotton areas in the area of the Junggar basin. With plastic mulching, theplanting boundary of early-middle mature cotton extended 65 km to the east approximately, and that of early mature cotton areas extended 0-300 km to the east. With the plastic mulching, the planting area of early-middle mature cotton (Yeta) increased by 5.47×104km2, some of which was from the former Taha cotton area. In addition, the planting areas of early mature and exceptional early-mature cotton decreased by 1.4% and 1.6% respectively. This study provide an effective method for cotton-planting zoning of different matures in Xinjiang.

    global warming; cotton; zoning; plastic mulching; Xinjiang; planting boundaries

    10.11975/j.issn.1002-6819.2019.02.012

    S162.5+4

    A

    1002-6819(2019)-02-0090-10

    2018-08-17

    2019-01-01

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2017YFD0101600);公益性行業(yè)科研專項(xiàng)(GYHY201206022)

    胡莉婷,博士生,主要從事氣候變化影響和農(nóng)業(yè)作物模型研究。Email:huliting@cau.edu.cn

    胡琦,博士,講師,主要從事氣候變化影響評(píng)價(jià)研究。Email:s10020292@cau.edu.cn

    胡莉婷,胡 琦,潘學(xué)標(biāo),馬雪晴,徐 琳,王香茹,張恒恒. 氣候變暖和覆膜對(duì)新疆不同熟性棉花種植區(qū)劃的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(2):90-99. doi:10.11975/j.issn.1002-6819.2019.02.012 http://www.tcsae.org

    Hu Liting, Hu Qi, Pan Xuebiao, Ma Xueqing, Xu Lin, Wang Xiangru, Zhang Hengheng. Effects of global warming and plastic mulching on cotton-planting zoning with different matures in Xinjiang[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(2): 90-99. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.02.012 http://www.tcsae.org

    猜你喜歡
    新疆區(qū)域
    走進(jìn)新疆
    國(guó)畫家(2022年2期)2022-04-13 09:07:46
    永久基本農(nóng)田集中區(qū)域“禁廢”
    在新疆(四首)
    分割區(qū)域
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    新疆多怪
    絲綢之路(2014年9期)2015-01-22 04:24:46
    區(qū)域
    民生周刊(2012年10期)2012-10-14 09:06:46
    新疆對(duì)外開放山峰
    91字幕亚洲| 久久久久久国产a免费观看| 琪琪午夜伦伦电影理论片6080| 岛国在线免费视频观看| 欧美黄色淫秽网站| 啦啦啦韩国在线观看视频| 婷婷精品国产亚洲av| 村上凉子中文字幕在线| 天天添夜夜摸| 在线观看日韩欧美| 国产精品野战在线观看| 99久久久亚洲精品蜜臀av| 1024香蕉在线观看| 精品久久久久久久末码| 亚洲精品乱码久久久v下载方式 | 色噜噜av男人的天堂激情| 亚洲狠狠婷婷综合久久图片| aaaaa片日本免费| 两个人看的免费小视频| 小说图片视频综合网站| 黄片小视频在线播放| 欧美成人性av电影在线观看| 国产激情久久老熟女| 不卡av一区二区三区| 97碰自拍视频| 偷拍熟女少妇极品色| 18禁黄网站禁片午夜丰满| 99热这里只有是精品50| 午夜视频精品福利| 亚洲中文av在线| 精品一区二区三区视频在线 | 在线国产一区二区在线| 男女午夜视频在线观看| 俄罗斯特黄特色一大片| 国产日本99.免费观看| 可以在线观看毛片的网站| 一区二区三区激情视频| 亚洲精品美女久久久久99蜜臀| a级毛片在线看网站| 色av中文字幕| av天堂中文字幕网| 岛国视频午夜一区免费看| 日韩高清综合在线| 精品久久久久久久久久免费视频| 国产视频一区二区在线看| 日本黄色视频三级网站网址| 最新中文字幕久久久久 | 夜夜夜夜夜久久久久| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 精品国产乱码久久久久久男人| 国产日本99.免费观看| 9191精品国产免费久久| 男人舔奶头视频| 视频区欧美日本亚洲| 国产精品美女特级片免费视频播放器 | 一区二区三区高清视频在线| 国产一区二区激情短视频| 欧美色视频一区免费| 脱女人内裤的视频| 国产av在哪里看| 国产精品自产拍在线观看55亚洲| 99久久综合精品五月天人人| 黄色成人免费大全| 五月玫瑰六月丁香| 午夜激情欧美在线| 免费在线观看视频国产中文字幕亚洲| 99re在线观看精品视频| 黄片大片在线免费观看| 啦啦啦韩国在线观看视频| 12—13女人毛片做爰片一| 哪里可以看免费的av片| 国产成人av教育| 少妇的逼水好多| 成人av一区二区三区在线看| 天天添夜夜摸| 色综合欧美亚洲国产小说| 欧美乱色亚洲激情| 在线观看一区二区三区| 禁无遮挡网站| 午夜日韩欧美国产| 欧美av亚洲av综合av国产av| 欧美黑人巨大hd| 久久伊人香网站| 三级男女做爰猛烈吃奶摸视频| 久久天堂一区二区三区四区| 成人av在线播放网站| 一级毛片精品| 亚洲国产精品sss在线观看| 观看美女的网站| 狂野欧美激情性xxxx| 国产免费男女视频| 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 国产在线精品亚洲第一网站| 成年女人看的毛片在线观看| 国产成年人精品一区二区| 1024手机看黄色片| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| 久久中文字幕人妻熟女| 人妻久久中文字幕网| 免费无遮挡裸体视频| 久久精品aⅴ一区二区三区四区| 欧美xxxx黑人xx丫x性爽| 久久午夜综合久久蜜桃| 国产乱人视频| 91av网站免费观看| 欧美中文综合在线视频| 麻豆久久精品国产亚洲av| 国产三级黄色录像| 国产精品亚洲av一区麻豆| 国内精品一区二区在线观看| 久久中文看片网| 国产真实乱freesex| 亚洲中文日韩欧美视频| 最新中文字幕久久久久 | 国产毛片a区久久久久| 国内揄拍国产精品人妻在线| 国产淫片久久久久久久久 | 97超级碰碰碰精品色视频在线观看| 国产伦在线观看视频一区| 欧美在线黄色| 嫁个100分男人电影在线观看| 久久久精品大字幕| 999久久久国产精品视频| 一个人免费在线观看的高清视频| 成熟少妇高潮喷水视频| 老汉色∧v一级毛片| 成人特级av手机在线观看| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看| 国产精品av视频在线免费观看| 久久亚洲真实| 91麻豆av在线| 久久久久久久久免费视频了| av视频在线观看入口| 色精品久久人妻99蜜桃| 精品无人区乱码1区二区| 国产av一区在线观看免费| 中文资源天堂在线| 在线a可以看的网站| 国产成人系列免费观看| 亚洲av成人av| 男女视频在线观看网站免费| 亚洲一区高清亚洲精品| 一区二区三区国产精品乱码| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 国产欧美日韩一区二区三| 日韩高清综合在线| 精品久久蜜臀av无| 99久久99久久久精品蜜桃| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区三| 欧美zozozo另类| 99久久精品国产亚洲精品| 最近最新中文字幕大全电影3| 窝窝影院91人妻| 两个人视频免费观看高清| 91在线精品国自产拍蜜月 | 天堂动漫精品| 亚洲欧美一区二区三区黑人| 久久久久九九精品影院| 可以在线观看的亚洲视频| 久久伊人香网站| 熟妇人妻久久中文字幕3abv| 久久精品91蜜桃| 久久国产乱子伦精品免费另类| 一区福利在线观看| 亚洲av电影在线进入| 两性夫妻黄色片| 一本精品99久久精品77| 精品福利观看| 国产真人三级小视频在线观看| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩高清专用| 欧美成人免费av一区二区三区| 欧美高清成人免费视频www| 日本黄色视频三级网站网址| 精华霜和精华液先用哪个| 我要搜黄色片| 日本黄色片子视频| 日韩欧美在线乱码| 成人av一区二区三区在线看| 岛国视频午夜一区免费看| 午夜福利视频1000在线观看| 校园春色视频在线观看| 久久久久国产精品人妻aⅴ院| 美女大奶头视频| 18禁美女被吸乳视频| 熟女人妻精品中文字幕| 婷婷丁香在线五月| 特级一级黄色大片| 老汉色∧v一级毛片| 国产野战对白在线观看| 久久香蕉精品热| 成人国产一区最新在线观看| 午夜福利在线在线| 在线视频色国产色| 90打野战视频偷拍视频| 欧美成人性av电影在线观看| 在线看三级毛片| 国产高清三级在线| 亚洲av五月六月丁香网| 少妇丰满av| 国产精品自产拍在线观看55亚洲| 久久香蕉国产精品| 母亲3免费完整高清在线观看| 美女 人体艺术 gogo| 怎么达到女性高潮| 国产精品99久久99久久久不卡| 亚洲真实伦在线观看| 国产精品日韩av在线免费观看| 久久久色成人| 日本在线视频免费播放| 国产精品乱码一区二三区的特点| 精品国产亚洲在线| 国内揄拍国产精品人妻在线| 99国产综合亚洲精品| 国产亚洲欧美98| 国产成人精品久久二区二区免费| 欧美不卡视频在线免费观看| 五月伊人婷婷丁香| 久久久精品欧美日韩精品| 午夜精品久久久久久毛片777| 精品久久久久久久久久久久久| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看| 亚洲精品美女久久av网站| 男人的好看免费观看在线视频| 亚洲第一电影网av| 色老头精品视频在线观看| 亚洲自拍偷在线| 岛国在线免费视频观看| 制服人妻中文乱码| 欧美三级亚洲精品| 国产亚洲欧美在线一区二区| 色播亚洲综合网| 琪琪午夜伦伦电影理论片6080| 麻豆国产av国片精品| 在线播放国产精品三级| 动漫黄色视频在线观看| 中文字幕人成人乱码亚洲影| 欧美成人性av电影在线观看| 国产成人啪精品午夜网站| 熟妇人妻久久中文字幕3abv| 国产三级中文精品| 久久久久久久久久黄片| 欧美一区二区国产精品久久精品| 免费观看的影片在线观看| 亚洲av电影在线进入| 热99在线观看视频| 听说在线观看完整版免费高清| 亚洲aⅴ乱码一区二区在线播放| 最新在线观看一区二区三区| 成人精品一区二区免费| 国产真实乱freesex| 国产v大片淫在线免费观看| 中文字幕高清在线视频| 亚洲国产看品久久| 91久久精品国产一区二区成人 | 欧美av亚洲av综合av国产av| 国产欧美日韩一区二区三| 97热精品久久久久久| 三级毛片av免费| 亚洲人成网站在线观看播放| 午夜日本视频在线| 中文字幕熟女人妻在线| 久久精品国产亚洲网站| 三级经典国产精品| 中文精品一卡2卡3卡4更新| 桃色一区二区三区在线观看| 国产69精品久久久久777片| 一级毛片aaaaaa免费看小| 一区二区三区高清视频在线| 一边亲一边摸免费视频| 中文精品一卡2卡3卡4更新| 91在线精品国自产拍蜜月| 国产精品久久电影中文字幕| 亚洲无线观看免费| 午夜a级毛片| 亚洲欧美日韩卡通动漫| 一卡2卡三卡四卡精品乱码亚洲| 国产精品.久久久| 黑人高潮一二区| 水蜜桃什么品种好| 日韩在线高清观看一区二区三区| ponron亚洲| 美女xxoo啪啪120秒动态图| 亚洲精品一区蜜桃| 国产精品综合久久久久久久免费| 99九九线精品视频在线观看视频| 国产午夜精品一二区理论片| 又粗又硬又长又爽又黄的视频| 如何舔出高潮| 亚洲精品aⅴ在线观看| 高清av免费在线| 国产黄片美女视频| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| www.av在线官网国产| 国内揄拍国产精品人妻在线| 校园人妻丝袜中文字幕| 成人亚洲精品av一区二区| 国产乱人视频| 精品久久久久久久末码| 欧美变态另类bdsm刘玥| 亚洲欧美日韩无卡精品| 嫩草影院入口| 国产精品不卡视频一区二区| av国产免费在线观看| 欧美zozozo另类| 波多野结衣巨乳人妻| 久久久成人免费电影| 99久国产av精品国产电影| 麻豆一二三区av精品| 91aial.com中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 国产成人a区在线观看| 亚洲国产精品sss在线观看| 久久6这里有精品| 一级av片app| 成人毛片a级毛片在线播放| 久久99热这里只频精品6学生 | 久久久久久久国产电影| 国产熟女欧美一区二区| 色哟哟·www| 欧美3d第一页| 看十八女毛片水多多多| 亚洲精品色激情综合| av女优亚洲男人天堂| 久久久久国产网址| 美女国产视频在线观看| av女优亚洲男人天堂| 久久久久久久久久成人| 亚洲av男天堂| 久久精品夜色国产| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 三级国产精品片| 波野结衣二区三区在线| 97在线视频观看| 天堂网av新在线| 亚洲精品色激情综合| 最近最新中文字幕大全电影3| 亚洲伊人久久精品综合 | 中文欧美无线码| 三级毛片av免费| 天天躁夜夜躁狠狠久久av| 成人性生交大片免费视频hd| 午夜福利高清视频| 少妇熟女欧美另类| 欧美日本亚洲视频在线播放| 国产v大片淫在线免费观看| 国产 一区精品| 久久精品久久久久久久性| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 成人毛片a级毛片在线播放| 在现免费观看毛片| 国内少妇人妻偷人精品xxx网站| 成人性生交大片免费视频hd| 久久欧美精品欧美久久欧美| 久久精品久久久久久久性| 亚洲乱码一区二区免费版| 99热全是精品| 成人三级黄色视频| 亚洲av电影不卡..在线观看| 久久久成人免费电影| 国产精品永久免费网站| 免费不卡的大黄色大毛片视频在线观看 | 高清在线视频一区二区三区 | 久久草成人影院| 美女大奶头视频| 如何舔出高潮| 日韩制服骚丝袜av| 久久人妻av系列| 免费av毛片视频| 老女人水多毛片| 日韩大片免费观看网站 | 乱码一卡2卡4卡精品| 日韩一本色道免费dvd| 欧美xxxx性猛交bbbb| 日韩av在线免费看完整版不卡| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 麻豆av噜噜一区二区三区| av播播在线观看一区| 干丝袜人妻中文字幕| 村上凉子中文字幕在线| 中文字幕免费在线视频6| 在线观看一区二区三区| 级片在线观看| 长腿黑丝高跟| 国产精品不卡视频一区二区| 岛国在线免费视频观看| 色网站视频免费| 日韩高清综合在线| 亚洲精品aⅴ在线观看| 最近2019中文字幕mv第一页| 2021天堂中文幕一二区在线观| 国产亚洲一区二区精品| 在线观看一区二区三区| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 久久精品熟女亚洲av麻豆精品 | 日韩成人av中文字幕在线观看| 国产免费男女视频| 亚洲高清免费不卡视频| 亚洲中文字幕一区二区三区有码在线看| 国产淫片久久久久久久久| 久久精品人妻少妇| 日本黄色片子视频| 国内揄拍国产精品人妻在线| 日本一本二区三区精品| 色哟哟·www| 亚洲精华国产精华液的使用体验| 免费观看人在逋| 99久久九九国产精品国产免费| av国产久精品久网站免费入址| 国产熟女欧美一区二区| 亚洲成人久久爱视频| 国产av在哪里看| 极品教师在线视频| 婷婷色麻豆天堂久久 | 又爽又黄a免费视频| 精品久久久久久久久亚洲| 久久久国产成人精品二区| 国产不卡一卡二| 日韩欧美国产在线观看| 男女边吃奶边做爰视频| 精品少妇黑人巨大在线播放 | 91久久精品电影网| 日本黄色片子视频| 亚洲怡红院男人天堂| 国产一区有黄有色的免费视频 | 日日摸夜夜添夜夜爱| 欧美日韩综合久久久久久| 精品欧美国产一区二区三| 亚洲激情五月婷婷啪啪| 女的被弄到高潮叫床怎么办| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 97在线视频观看| 亚洲,欧美,日韩| 卡戴珊不雅视频在线播放| 精品国产一区二区三区久久久樱花 | 又爽又黄无遮挡网站| 亚洲精品,欧美精品| 国产老妇伦熟女老妇高清| 免费看日本二区| 看免费成人av毛片| 熟女电影av网| av在线老鸭窝| 99久久精品热视频| 内地一区二区视频在线| 草草在线视频免费看| 久久久国产成人免费| 欧美区成人在线视频| 日韩av在线大香蕉| 国产伦一二天堂av在线观看| 欧美性猛交黑人性爽| 99久久精品一区二区三区| 又粗又硬又长又爽又黄的视频| 国产 一区 欧美 日韩| 亚洲美女视频黄频| 99久国产av精品| 中文欧美无线码| 国产乱来视频区| 日韩精品青青久久久久久| 成人av在线播放网站| 国产精品一区二区性色av| 午夜免费激情av| 亚洲经典国产精华液单| 日韩亚洲欧美综合| 2021天堂中文幕一二区在线观| 国产精品久久电影中文字幕| 一边亲一边摸免费视频| 日本熟妇午夜| 毛片女人毛片| 村上凉子中文字幕在线| 午夜久久久久精精品| 国产成人精品婷婷| 国产伦在线观看视频一区| 久久午夜福利片| 22中文网久久字幕| 欧美极品一区二区三区四区| АⅤ资源中文在线天堂| 欧美高清成人免费视频www| 色播亚洲综合网| 日本-黄色视频高清免费观看| 欧美色视频一区免费| 日韩欧美精品v在线| 国产大屁股一区二区在线视频| 男女国产视频网站| 国产精品1区2区在线观看.| 成人性生交大片免费视频hd| 3wmmmm亚洲av在线观看| 联通29元200g的流量卡| 国产高清有码在线观看视频| 国产免费一级a男人的天堂| 婷婷色av中文字幕| 国产欧美另类精品又又久久亚洲欧美| 日韩成人伦理影院| 欧美xxxx黑人xx丫x性爽| 在线观看av片永久免费下载| 亚洲久久久久久中文字幕| 亚洲内射少妇av| 人妻制服诱惑在线中文字幕| 久久久国产成人免费| 精华霜和精华液先用哪个| 级片在线观看| 乱人视频在线观看| 国产高清不卡午夜福利| 午夜激情福利司机影院| 精品久久久久久电影网 | 亚州av有码| 亚洲无线观看免费| 春色校园在线视频观看| 亚洲欧美精品专区久久| 真实男女啪啪啪动态图| 深爱激情五月婷婷| 九九在线视频观看精品| 日本黄色片子视频| 爱豆传媒免费全集在线观看| 麻豆成人av视频| 午夜福利视频1000在线观看| 国产精品无大码| 99热这里只有是精品50| 日本爱情动作片www.在线观看| 中文亚洲av片在线观看爽| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 小蜜桃在线观看免费完整版高清| 好男人视频免费观看在线| 午夜福利高清视频| 亚洲精品乱码久久久v下载方式| 国产亚洲午夜精品一区二区久久 | 国产伦在线观看视频一区| 最近手机中文字幕大全| 村上凉子中文字幕在线| 精品人妻一区二区三区麻豆| 舔av片在线| 亚洲中文字幕日韩| 亚洲一区高清亚洲精品| 亚洲av福利一区| 久久久久国产网址| 精品人妻熟女av久视频| 久久久成人免费电影| 亚洲国产日韩欧美精品在线观看| 一级毛片我不卡| 国产精品福利在线免费观看| 亚洲av成人精品一区久久| 国产黄片视频在线免费观看| 91久久精品电影网| av线在线观看网站| 青春草亚洲视频在线观看| 国产精品久久久久久久电影| 国产一级毛片在线| 国产一区二区在线av高清观看| 亚洲经典国产精华液单| 麻豆乱淫一区二区| 黑人高潮一二区| 九草在线视频观看| 51国产日韩欧美| 亚洲国产精品专区欧美| 国产精品不卡视频一区二区| 久久国内精品自在自线图片| eeuss影院久久| 天堂av国产一区二区熟女人妻| 精品免费久久久久久久清纯| ponron亚洲| 内地一区二区视频在线| 最近视频中文字幕2019在线8| 日本黄色视频三级网站网址| 特大巨黑吊av在线直播| 久久久精品94久久精品| 熟女电影av网| 特大巨黑吊av在线直播| 男女国产视频网站| 日本免费在线观看一区| 亚洲av不卡在线观看| 九九热线精品视视频播放| 夜夜看夜夜爽夜夜摸| 久久精品国产99精品国产亚洲性色| 高清av免费在线| av福利片在线观看| 亚洲一区高清亚洲精品| 床上黄色一级片| 色网站视频免费| 嘟嘟电影网在线观看| 久久久久久国产a免费观看| 人体艺术视频欧美日本| 如何舔出高潮| 综合色av麻豆| 亚洲中文字幕一区二区三区有码在线看| 欧美另类亚洲清纯唯美| 成人特级av手机在线观看| 久久精品91蜜桃| 精品99又大又爽又粗少妇毛片| 亚洲精品国产成人久久av| 日本-黄色视频高清免费观看| 国产91av在线免费观看| 国产精品精品国产色婷婷| 国内精品美女久久久久久| 别揉我奶头 嗯啊视频| 在线观看美女被高潮喷水网站| 国产高潮美女av| 亚洲国产精品国产精品| 嫩草影院新地址| 午夜日本视频在线| 天堂网av新在线| 午夜福利在线观看吧| 日本色播在线视频|