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

    未來氣候變化情景下東北地區(qū)中晚熟春玉米機械粒收氣候適宜區(qū)分析*

    2023-08-24 02:44:42蘇正娥劉志娟楊婉蓉祝光欣史登宇楊曉光
    中國農(nóng)業(yè)氣象 2023年8期
    關(guān)鍵詞:吉林省籽粒遼寧省

    蘇正娥,劉志娟,楊婉蓉,祝光欣,史登宇,楊曉光

    未來氣候變化情景下東北地區(qū)中晚熟春玉米機械粒收氣候適宜區(qū)分析*

    蘇正娥,劉志娟**,楊婉蓉,祝光欣,史登宇,楊曉光

    (中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100193)

    玉米生產(chǎn)全過程機械化是玉米生產(chǎn)的發(fā)展方向,當(dāng)前機械粒收是中國轉(zhuǎn)變玉米生產(chǎn)方式的重點。籽粒含水率是影響機械粒收的主要因素,其主要受玉米農(nóng)藝性狀、生育后期生態(tài)氣象因子以及栽培管理措施等因素影響。東北地區(qū)作為中國玉米主產(chǎn)區(qū)之一,也是對氣候變化最敏感的區(qū)域之一,因此,明確未來氣候條件下該地區(qū)玉米機械粒收的氣候適宜性及適宜區(qū),可為未來進(jìn)一步提高玉米生產(chǎn)的機械化水平,實現(xiàn)玉米生產(chǎn)全過程機械化提供科學(xué)依據(jù)。本文以RCP4.5和RCP8.5兩個溫室氣體濃度情景,分別代表中等排放情景和高排放情景,結(jié)合兩種情景下預(yù)估的2021?2060年氣候數(shù)據(jù)和春玉米生育期數(shù)據(jù),利用籽粒含水率預(yù)測模型確定未來40a當(dāng)前廣泛種植的春玉米中晚熟品種機械粒收各等級氣候適宜區(qū)北界,并基于此得出基于多品種的機械粒收氣候適宜區(qū)。結(jié)果表明:未來40a不同品種的氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)表現(xiàn)為由西南向東北延伸的分布特征,不同品種間各級適宜區(qū)北界由西南向吉林省東南部山區(qū)的長白、東崗和松江一帶以及高緯度的黑龍江北部變動。RCP4.5(RCP8.5)氣候情景下,50%和80%春玉米中晚熟品種的最適宜區(qū)均位于遼寧省南部(遼寧省西部和南部),最北分別可達(dá)42.0°N和41.0°N(42.8°N和41.9°N),最東可達(dá)123.4°E和123.3°E(124.5°E和123.3°E),面積分別占研究區(qū)域面積的5.9%和1.8%(11.2%和5.9%);適宜區(qū)分別位于遼寧省北部和中部(遼寧省北部和吉林省西南部),最北可至43.0°N和42.7°N(44.8°N和42.9°N),最東可至124.7°E和124.4°E(124.7°E和124.5°E),該區(qū)域面積占比分別為8.3%和8.9%(4.7%和6.6%)。未來40a,與RCP4.5氣候情景相比,RCP8.5氣候情景下春玉米中晚熟品種機械粒收氣候適宜性更強,但兩個氣候情景下氣候最適宜區(qū)和適宜區(qū)范圍均較小。未來春玉米機械粒收氣候適宜區(qū)在品種間差異較大,未來可選育種植生育期更短、脫水速率更快的適宜機械粒收品種,以提高機械粒收的質(zhì)量和效率。

    春玉米;籽粒脫水;機械粒收;氣候適宜區(qū);氣候變化;東北地區(qū)

    玉米是世界主要糧食作物之一,中國是世界第二大玉米生產(chǎn)國,近5a(2017-2021年)玉米播種面積維持在4100萬hm2以上,總產(chǎn)量穩(wěn)定在2.5億t以上,播種面積和產(chǎn)量均占中國糧食生產(chǎn)的35%以上[1]。當(dāng)前,生產(chǎn)全過程機械化是中國玉米生產(chǎn)的主要發(fā)展方向,其中利用機械進(jìn)行玉米籽粒收獲(簡稱機械粒收)是玉米機械化生產(chǎn)的重要環(huán)節(jié)。機械粒收不僅可大大節(jié)約勞動力成本,還可降低晾曬、脫粒過程中的籽粒霉?fàn)€與損失,是玉米機械收獲的發(fā)展方向和今后玉米生產(chǎn)方式轉(zhuǎn)變的重點[2]。

    相較于機械穗收和人工收獲,機械粒收每公頃分別可節(jié)約成本375~600元和1050~1650元,若以單產(chǎn)6t·hm?2計算,每生產(chǎn)1t玉米成本分別降低62.5~100.0元和175~275元[3]。美國、德國等國在20世紀(jì)50年代玉米收獲以機械穗收為主,到70年代則全面采用田間機械粒收,60-90年代集中開展了玉米機械粒收技術(shù)相關(guān)基礎(chǔ)理論及品種選育的研究,為機械粒收的全面推廣提供了支撐[2,4?6]。目前中國玉米收獲以人工收獲和機械穗收為主,而機械粒收仍處于起步階段,相關(guān)研究多集中于機械粒收質(zhì)量及其影響因素分析和適宜品種篩選[7?13]。機械粒收質(zhì)量主要指標(biāo)包括籽粒破碎率、雜質(zhì)率和損失率。其中破碎率主要受籽粒含水率和籽粒力學(xué)強度影響,雜質(zhì)主要成分為穗軸,雜質(zhì)率與穗軸含水率及力學(xué)強度密切相關(guān)[14]。研究發(fā)現(xiàn),籽粒含水率是影響機械粒收的主要因素,與籽粒破碎率、損失率和雜質(zhì)率呈顯著相關(guān),即籽粒含水率越高,籽粒破碎率和雜質(zhì)率越高,田間損失率越低[15?18]。由于收獲時籽粒含水率是籽粒破碎率的主要決定因素[19],大量學(xué)者對玉米籽粒脫水及籽粒含水率的影響因素進(jìn)行了深入探究,認(rèn)為收獲時籽粒含水率主要由生理成熟前后籽粒脫水速率控制,而脫水速率主要受玉米農(nóng)藝性狀、生育后期生態(tài)氣象因子以及栽培管理措施等因素影響[20]。此外,研究發(fā)現(xiàn),玉米籽粒含水率與授粉后積溫呈極顯著非線性關(guān)系[21]?;诓煌贩N的含水率預(yù)測模型,有學(xué)者進(jìn)一步估算了含水率降到機械粒收理想含水率時所需的活動積溫,確定了不同地區(qū)適宜機械粒收的時間[22?25]。因此通過選育脫水速率快、生理成熟期含水率低的品種和適當(dāng)推遲收獲時間可降低收獲時籽粒含水率,從而降低機械粒收籽粒破碎率,提高機械粒收質(zhì)量[12]。

    當(dāng)前關(guān)于中國機械粒收的研究為國內(nèi)玉米機械粒收的推廣與應(yīng)用奠定了科學(xué)基礎(chǔ),為提高機械粒收質(zhì)量提供了理論保障。但基于區(qū)域氣象條件,評估機械粒收生產(chǎn)技術(shù)適應(yīng)性及適宜區(qū)的研究還鮮見報道,特別是氣候變化背景下未來機械粒收氣候適宜性的研究還未見到。東北地區(qū)作為中國玉米主產(chǎn)區(qū)之一,也是對氣候變化最敏感的區(qū)域之一,明確氣候變化背景下機械粒收在該地區(qū)的適宜性對未來進(jìn)一步提高玉米生產(chǎn)的機械化水平,實現(xiàn)玉米生產(chǎn)全程機械化至關(guān)重要。

    本文基于2021?2060年氣象數(shù)據(jù)和春玉米生育期數(shù)據(jù),利用籽粒含水率預(yù)測模型,預(yù)測初霜日前春玉米的籽粒含水率,確定未來40a不同春玉米品種機械粒收的各等級適宜區(qū)北界以及基于多品種的機械粒收氣候適宜區(qū),以期為東北地區(qū)合理推廣與應(yīng)用春玉米機械粒收技術(shù),提高春玉米生產(chǎn)效率提供參考。

    1 資料與方法

    1.1 研究區(qū)域概況

    以東北三?。ê邶埥?、吉林和遼寧?。┐河衩诐撛诜N植區(qū)為研究區(qū)域。春玉米潛在種植區(qū)是根據(jù)80%保證率下≥10℃活動積溫高于2100℃·d確定的,包括黑龍江省中南部、吉林省和遼寧省。該地區(qū)地勢平坦、土壤肥沃,氣候為溫帶濕潤半濕潤大陸性季風(fēng)氣候,光水資源優(yōu)越,唯熱量條件略顯不足,年均氣溫為?3.8~11.3℃,喜溫作物溫度生長期內(nèi)≥10℃活動積溫為1680~3850℃·d,日照時數(shù)806~1545h,降水量為298~880mm[26?27]。研究區(qū)域內(nèi)包含67個氣象觀測站點和45個農(nóng)業(yè)氣象觀測站點,分別用于收集研究區(qū)內(nèi)的氣象數(shù)據(jù)和作物生產(chǎn)數(shù)據(jù)。研究區(qū)和所有觀測站點分布見圖1。

    圖1 研究區(qū)域(東北三省春玉米潛在種植區(qū))氣象站點和農(nóng)業(yè)氣象站點分布

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

    潛在種植區(qū)內(nèi)春玉米實際種植面積提取自2000?2015年全國玉米收獲面積高精度(1km)遙感數(shù)據(jù)集[28],利用ArcGIS10.6軟件對逐年全國玉米收獲面積進(jìn)行疊加并基于研究區(qū)域矢量圖層提取。

    政府間氣候變化專門委員會(IPCC)為了對未來氣候進(jìn)行評估,以2100年輻射強迫水平達(dá)到4.5W·m?2和8.5W·m?2分別代表中等排放情景(RCP4.5)和高排放情景(RCP8.5)。采用的未來氣候情景數(shù)據(jù)來源于英國氣象局開發(fā)的區(qū)域氣候模擬系統(tǒng)PRECIS在RCP4.5和RCP8.5情景下模擬的未來時段(2021?2060年)逐日氣象數(shù)據(jù),模式時間積分步長為5min,水平輸出分辨率在旋轉(zhuǎn)坐標(biāo)下的經(jīng)緯度各為0.44°。該數(shù)據(jù)集氣象要素包括日平均溫度、最高氣溫、最低氣溫、降水量、太陽總輻射和平均風(fēng)速,對上述氣象要素輸出結(jié)果進(jìn)行訂正,并插值為空間分辨率0.5°×0.5°的網(wǎng)格數(shù)據(jù),由中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所提供,基于此數(shù)據(jù)集提取研究區(qū)域內(nèi)67個氣象站的未來氣候數(shù)據(jù)。

    春玉米生育期數(shù)據(jù)來源于研究區(qū)域內(nèi)45個農(nóng)業(yè)氣象觀測站,包括2010?2015年春玉米抽雄期和成熟期。利用反距離加權(quán)法插值獲得各氣象站點逐年生育期,并假設(shè)未來40a春玉米生育期不變,取2010?2015年各氣象站6a平均值作為未來40a春玉米的平均生育期。

    1.3 機械粒收氣候適宜區(qū)的判別方法

    1.3.1 計算籽粒含水率

    玉米籽粒脫水過程不僅受生長發(fā)育控制,還與環(huán)境因素有關(guān)。因此,在玉米籽粒形成過程中籽粒含水率(MC,%)由生理控制的趨勢含水率(Trend_ MC,%)和由環(huán)境因素控制的氣象含水率(Climate_ MC,%)組成,三者關(guān)系為[29]

    其中,趨勢含水率用Logistic Power含水率動態(tài)預(yù)測模型模擬[21],即

    式中,a、b、c為模型參數(shù),a為籽粒初始含水率(%),設(shè)定為90;b和c為由品種特性決定的品種參數(shù)。采用東北地區(qū)主栽的多個中晚熟春玉米品種進(jìn)行機械粒收氣候適宜性分析,相較于單一品種,可降低品種差異對機械粒收氣候適宜性的影響,使結(jié)果更具有區(qū)域代表性。本文基于黃兆福等[25]確定的東北地區(qū)29個主栽中晚熟春玉米品種參數(shù),模擬當(dāng)?shù)?9個春玉米主栽品種的趨勢含水率;TT為授粉? 計算籽粒含水率期間≥0℃積溫(℃·d)。

    根據(jù)玉米籽粒脫水特征,春玉米脫水可分為兩個階段,第一階段為生理成熟前與灌漿相關(guān)的“生理脫水”,第二階段為生理成熟后籽粒水分的散失過程,兩個階段籽粒脫水的主要環(huán)境影響因子存在差異[30?31]。氣象含水率采用高尚等[29]基于田間試驗確定的回歸模型進(jìn)行模擬,分為生理成熟前(Climate_MC1)和生理成熟后(Climate_MC2),預(yù)測模型分別為

    式中,x1、x2和x3分別為授粉?生理成熟期的日平均氣溫(℃)、平均風(fēng)速(m·s?1)和平均蒸發(fā)量(mm);x4和x5分別為生理成熟期?測定含水率日期的日平均氣溫(℃)和平均空氣相對濕度,其中相對濕度采用小數(shù)形式。日蒸發(fā)量用參考作物蒸散量表示,采用Penman-Monteith公式計算[32],即

    式中,ET0為日參考作物蒸散量(mm);Rn為到達(dá)作物表面的凈輻射(MJ·m?2·d?1);G為土壤熱通量密度(MJ·m?2·d?1),取值為0;T為作物冠層2m高處的空氣溫度(℃);U2為作物冠層2m高處的風(fēng)速(m·s?1);ed為飽和水汽壓(kPa);eα為實際水汽壓(kPa);Δ為水汽壓對溫度的斜率(kPa·℃?1);γ為干濕球常數(shù)。其中,Rn、G、Δ和U2通過氣象臺站觀測資料計算求得。

    1.3.2 劃分機械粒收氣候適宜區(qū)的標(biāo)準(zhǔn)

    以秋季日最低氣溫低于0℃的第一天確定各氣象站未來40a(2021?2060年)逐年初霜日[33],基于以上籽粒含水率計算方法,分別計算各年份29個中晚熟品種初霜日前一天的春玉米籽粒含水率。在某一站點,若某一品種某年初霜日前籽粒含水率可降至25%(即MC≤25),則認(rèn)為該地區(qū)在該年份對這一品種可進(jìn)行機械粒收[25]。在此基礎(chǔ)上考慮累積概率,進(jìn)行單一品種機械粒收適宜區(qū)的劃分,即將未來40a內(nèi)有0~20%、20%~50%、50%~80%和80%以上年份可進(jìn)行機械粒收的區(qū)域分別劃分為該品種的機械粒收氣候不適宜區(qū)、次適宜區(qū)、適宜區(qū)和最適宜區(qū)。

    對于各個地區(qū),如果29個品種均屬于同一等級適宜區(qū),則定義該地區(qū)為春玉米中晚熟品種的機械粒收該等級適宜區(qū),如某一站點,所有品種均為最適宜區(qū),則該地區(qū)為春玉米中晚熟品種的機械粒收最適宜區(qū)。如果不同品種所屬適宜區(qū)等級不一致,則定義為敏感帶。最適宜區(qū)與適宜區(qū)的敏感帶為敏感帶一,即該地區(qū)對于一部分品種屬于最適宜區(qū),對于另一部分品種則屬于適宜區(qū)。類似地,最適宜區(qū)、適宜區(qū)與次適宜區(qū)的敏感帶為敏感帶二,次適宜區(qū)與不適宜區(qū)的敏感帶為敏感帶三。本文確定了50%與80%品種保證率的機械粒收氣候適宜區(qū),即得出的機械粒收氣候適宜區(qū)分別對于50%與80%品種是滿足的,生產(chǎn)者可基于自身的風(fēng)險承擔(dān)程度選擇是否選用機械進(jìn)行春玉米中晚熟品種的籽粒收獲。

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

    采用R語言計算初霜日和籽粒含水率,采用ArcGIS 10.6制作氣候適宜區(qū)空間分布示意圖。

    2 結(jié)果與分析

    2.1 未來氣候情景下單一品種機械粒收氣候適宜區(qū)分析

    由圖2可見,RCP4.5氣候情景下,東北地區(qū)29個中晚熟春玉米品種機械粒收的氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)均表現(xiàn)出相同特點即依次由西南向東北延伸,適宜性整體較差,氣候最適宜區(qū)和適宜區(qū)主要位于遼寧省和吉林省。不同春玉米品種間,機械粒收氣候適宜區(qū)空間差異較大,其中迪卡159、ND865、澤玉501和東單507機械粒收氣候適宜性最差,僅大連最適宜機械粒收,適宜區(qū)位于遼寧省最南部。金珠58、德育919、澤玉8911和京農(nóng)科728等品種機械粒收氣候適宜性較高,氣候不適宜區(qū)主要位于海拔較高的吉林省東南部、黑龍江省東南部和緯度較高的黑龍江省北部地區(qū),氣候最適宜區(qū)與適宜區(qū)位于遼寧省、吉林省西部和黑龍江省西南部。圖3顯示,RCP8.5氣候情景下機械粒收氣候適宜區(qū)空間分布特征與RCP4.5一致,但相較于RCP4.5,RCP8.5氣候情景下最適宜區(qū)與適宜區(qū)面積更大,機械粒收氣候適宜性更高。

    2.2 未來氣候情景下機械粒收氣候適宜區(qū)北界及其敏感帶分析

    圖4為RCP4.5氣候情景下未來40a(2021?2060年)東北地區(qū)不同春玉米中晚熟品種機械粒收氣候適宜區(qū)北界變動范圍及敏感帶分布。其中,圖4a、b和c分別為最適宜區(qū)、適宜區(qū)和次適宜區(qū)北界變動范圍。由圖可知,機械粒收氣候最適宜區(qū)北界范圍主要為遼寧省除清原、桓仁外的所有區(qū)域和吉林省西部,即黑龍江省泰來、吉林省前郭爾羅斯、四平、遼寧省開原、章黨、本溪和寬甸一帶及其以西地區(qū)。機械粒收氣候適宜區(qū)北界變動范圍位于吉林省除南部外的所有區(qū)域和吉林省西部以及黑龍江省西南部和東部的富錦、寶清和虎林一帶,即適宜區(qū)變動范圍南界為遼寧省興城、錦州、營口和莊河一帶以北,北界為黑龍江省的齊齊哈爾、安達(dá)、哈爾濱、吉林省的三岔河、長春、清原和桓仁一帶。機械粒收氣候次適宜區(qū)北界變動范圍位于遼寧省東北部、吉林省除東南部外的所有區(qū)域和黑龍江省,變動范圍的南界位于遼寧省開原、章黨、本溪、寬甸一帶以南,最東可至吉林省東南部山地地區(qū),即臨江、靖宇、敦化、牡丹江和延吉一帶。圖4d為RCP4.5氣候情景下機械粒收氣候適宜區(qū)北界變動范圍敏感帶。由圖可知,遼寧省興城、錦州、營口和莊河一帶以北,彰武、沈陽和丹東一帶及其以南為敏感帶一(最適宜區(qū)與適宜區(qū)敏感帶),即對部分品種該地區(qū)為最適宜區(qū),對于部分為適宜區(qū)。黑龍江省泰來、吉林省前郭爾羅斯、四平一帶及其以西和遼寧省的寬甸、本溪、章黨和開原一帶為敏感帶二(最適宜區(qū)、適宜區(qū)與次適宜區(qū)的敏感帶),即對于部分品種為最適宜區(qū),部分為適宜區(qū),部分為次適宜區(qū)。黑龍江省齊齊哈爾、安達(dá)、哈爾濱、吉林省三岔河、長春、遼寧省清原和桓仁一帶為敏感帶三(次適宜區(qū)與不適宜區(qū)的敏感帶),即對于部分品種該帶為次適宜區(qū),另一部分為不適宜區(qū)。

    圖5為RCP8.5氣候情景下未來40a(2021?2060年)東北地區(qū)春玉米中晚熟品種機械粒收氣候適宜區(qū)北界變動范圍及敏感帶分布。圖5a、b和c分別為最適宜區(qū)、適宜區(qū)和次適宜區(qū)北界變動范圍,圖5d為氣候適宜區(qū)北界變動范圍敏感帶。由圖可知,RCP8.5氣候情景下,不同品種最適宜區(qū)北界變動范圍為黑龍江省泰來、吉林省前郭爾羅斯、長春、章黨和寬甸一帶及其以西地區(qū)。適宜區(qū)北界范圍為朝陽、阜新、沈陽、本溪和岫巖一帶及其以北,黑龍江齊齊哈爾、安達(dá)、綏化、哈爾濱、吉林省梅河口、遼寧省桓仁一帶及其以西地區(qū)以及黑龍江省佳木斯、富錦、寶清、虎林和綏芬河一帶。次適宜區(qū)北界變動范圍與RCP4.5氣候情景下一致。朝陽、阜新、彰武、沈陽和岫巖一帶為敏感帶一。敏感帶二與RCP4.5氣候情景下一致,敏感帶三相較于RCP4.5氣候情景增加了吉林省梅河口和黑龍江省佳木斯和綏芬河等地。

    注:對任一品種,最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)分別指未來40a內(nèi)有80%以上、50%~80%、20%~50%和0~20%年份可進(jìn)行機械粒收的區(qū)域,該品種在某年可進(jìn)行機械粒收的判斷指標(biāo)為:該年初霜日前玉米籽粒含水率可降至25%(即MC≤25%)。下同。

    Note: For each variety, the very suitable, suitable, moderately suitable and marginally suitable zones refer to the regions where maize kernel can be harvested mechanically with more than 80%, 50%?80%, 20%?50% and 0?20% of the future 40 years . The judgment index for mechanical kernel harvesting of this variety in a certain year is that the grain moisture content can be reduced to 25% (i.e. MC≤25%) before the first frost day. The same as below.

    圖3 RCP8.5氣候情景下未來(2021?2060年)中晚熟春玉米品種(29個)機械粒收氣候適宜區(qū)

    2.3 未來氣候情景下基于多品種的機械粒收氣候適宜區(qū)綜合劃分

    圖6為RCP4.5和RCP8.5氣候情景下未來40a(2021?2060年)東北地區(qū)50%與80%中晚熟春玉米品種的機械粒收氣候適宜區(qū)分布及面積占比。表1為各氣候適宜區(qū)的經(jīng)緯度范圍和包含的實際春玉米種植面積及占比。結(jié)合圖6和表1可知,RCP4.5氣候情景下,50%春玉米品種的機械粒收氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)面積分別占研究區(qū)域面積的5.9%、8.3%、18.9%和66.8%。其中,最適宜區(qū)位于遼寧省南部,即綏中、興城、錦州、黑山、鞍山和莊河一帶及其以南,最北可達(dá)北緯42.0°,最東可達(dá)東經(jīng)123.4°。適宜區(qū)位于遼寧省北部與中部,即葉柏壽、朝陽、阜新、彰武、沈陽、本溪、岫巖和丹東一帶以及吉林省雙遼地區(qū),最北可達(dá)北緯43.0°,最東可至東經(jīng)124.7°。次適宜區(qū)位于黑龍江省西南部和吉林省西北部以及遼寧省東部,即黑龍江省泰來、安達(dá)、哈爾濱、吉林省長春和四平一帶及其以西以及遼寧省的開原、章黨和寬甸一帶,最北可達(dá)北緯47.8°,最東可至東經(jīng)127.1°。不適宜區(qū)位于黑龍江省除西南部外的所有區(qū)域以及吉林省東部、南部和遼寧省清原、桓仁等地。80%春玉米品種機械粒收氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)面積占研究區(qū)域面積分別為1.8%、8.9%、10.9%和78.5%。最適宜區(qū)位于遼寧省南部的營口、熊岳、瓦房店和莊河一帶及其以南地區(qū),最北可達(dá)北緯41.0°,最東可至東經(jīng)123.3°。適宜區(qū)位于遼寧省西部和中部,即阜新、黑山、沈陽、鞍山、丹東一帶及其以西地區(qū),最北可達(dá)北緯42.7°,最東可至東經(jīng)124.4°。次適宜區(qū)位于吉林省西部、遼寧省中部地區(qū),即吉林省白城、前郭爾羅斯、雙遼和四平一帶及其以西以及遼寧開原、彰武、本溪和岫巖一帶,最北可達(dá)北緯46.2°,最東可至東經(jīng)125.2°。不適宜區(qū)位于黑龍江省、吉林省除西北外的所有區(qū)域以及遼寧省東部地區(qū)。

    圖4 RCP4.5氣候情景下未來(2021?2060年)東北地區(qū)中晚熟春玉米品種機械粒收氣候適宜區(qū)北界范圍及敏感帶分布

    注:圖a、b和c分別為29個中晚熟春玉米品種氣候最適宜區(qū)、適宜區(qū)和次適宜區(qū)北界范圍,圖d顯示各適宜區(qū)的敏感帶,最適宜區(qū)與適宜區(qū)的敏感帶為敏感帶一,即該地區(qū)對于一部分品種屬于最適宜區(qū),對于另一部分品種則屬于適宜區(qū)。類似地,最適宜區(qū)、適宜區(qū)與次適宜區(qū)的敏感帶為敏感帶二,次適宜區(qū)與不適宜區(qū)的敏感帶為敏感帶三。下同。

    Note: Fig.a, Fig.b and Fig.c are respectively the ranges of northern boundaries of very suitable, suitable and moderately suitable zone for 29 middle-late maturing spring maize varieties, Fig.d shows the sensitive regions of suitable zones. The sensitive region of the very suitable and suitable zones is the sensitive region Ⅰ, that is, the region belongs to very suitable zone for some varieties and suitable zone for others. Similarly, the sensitive region of very suitable, suitable and moderately suitable zone is the sensitive zone II, and the sensitive region of moderately suitable zone and marginally suitable zone is the sensitive zone III. The same as below.

    RCP8.5氣候情景下,50%品種機械粒收最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)面積分別占研究區(qū)域面積的11.2%、4.7%、18.0%和66.2%。最適宜區(qū)位于遼寧省中部、西部和南部,即彰武、沈陽、鞍山、莊河和丹東一帶及其以西以南地區(qū),最北可達(dá)北緯42.8°,最東可至東經(jīng)124.5°。適宜區(qū)位于吉林省雙遼、四平、遼寧省本溪和岫巖一帶,最北可達(dá)北緯44.8°,最東可至東經(jīng)124.72°。次適宜區(qū)位于黑龍江省泰來、哈爾濱、吉林省三岔河、長春一帶及其以西地區(qū)以及遼寧省開原、章黨和寬甸一帶,最北可達(dá)北緯47.7°,最東可至東經(jīng)127.2°。不適宜區(qū)位于黑龍江省齊齊哈爾、安達(dá)、尚志、吉林省梅河口、遼寧省清原、桓仁一帶及其以東以北地區(qū)。80%春玉米品種機械粒收氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)面積分別占研究區(qū)域面積的5.9%、6.6%、12.8%和74.7%。最適宜區(qū)位于遼寧省南部的葉柏壽、錦州、黑山、鞍山和莊河一帶及其以南地區(qū),最北可達(dá)41.9°,最東可至東經(jīng)123.3°。適宜區(qū)位于遼寧省朝陽、阜新、彰武、沈陽、岫巖和丹東一帶,最北可達(dá)北緯42.9°,最東可至東經(jīng)124.5°。次適宜區(qū)為黑龍江省泰來、吉林省前郭爾羅斯、長春、四平一帶及其以西,以及遼寧省開原和本溪一帶,最北可達(dá)北緯46.7°,最東可至東經(jīng)125.3°。不適宜區(qū)位于遼寧省東北部、吉林省中部、東部和黑龍江省。

    圖5 RCP8.5氣候情景下未來(2021?2060年)東北地區(qū)中晚熟春玉米品種機械粒收氣候適宜區(qū)北界范圍及敏感帶分布

    圖6 RCP4.5(a)和RCP8.5(b)氣候情景下未來(2021?2060)東北地區(qū)50%(1)與80%(2)中晚熟春玉米品種機械粒收氣候適宜區(qū)劃分及其面積占比

    注:餅狀圖為各級適宜區(qū)面積占研究區(qū)域即春玉米潛在種植區(qū)面積的比例。

    Note:The pie chart shows the proportion of the suitable areas at all levels to the potential planting area of spring maize in the study area.

    表1 多品種綜合機械粒收氣候適宜區(qū)的經(jīng)緯度范圍及其目前實際種植面積和占比

    注:“面積”為各級適宜區(qū)內(nèi)包含的春玉米實際種植面積,“占比”為各級適宜區(qū)內(nèi)春玉米實際種植面積占整個研究區(qū)域內(nèi)春玉米實際總種植面積的比例。

    Note: "Area" in the table refers to the actual planting area of spring maize included in the appropriate areas, and "Percentage" refers to the proportion of the actual planting area of spring maize in the appropriate areas to the actual total planting area of spring maize in the whole study area.

    RCP4.5氣候情景下,50%春玉米中晚熟品種的機械粒收氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)分別占研究區(qū)域內(nèi)春玉米總種植面積的7.9%、11.2%、24.9%和55.9%,80%品種的各級氣候適宜區(qū)占比分別為2.2%、11.9%、13.7%和72.3%。RCP8.5氣候情景下,50%春玉米中晚熟品種的各級機械粒收氣候適宜區(qū)占比分別為15.3%、6.1%、24.1%和54.6%,80%春玉米品種各級適宜區(qū)占比分別為7.7%、8.8%、16.9%和66.6%。在RCP4.5氣候情景下,最適宜區(qū)與適宜區(qū)內(nèi)包含的春玉米種植面積占研究區(qū)域總種植面積的比例在50%與80%品種保證率下分別為19.1%和14.0%,RCP8.5氣候情景下分別為21.4%和16.6%,說明RCP8.5氣候情景下春玉米機械粒收氣候適宜性較RCP4.5情景強。在兩個氣候情景下,50%與80%品種保證率各級適宜區(qū)經(jīng)緯度范圍以及各區(qū)域面積占比存在較大差異,因此未來可選育種植生育期更短、脫水速率更快的適宜機械粒收品種,以提高機械粒收的質(zhì)量和效率。

    3 結(jié)論與討論

    3.1 討論

    近50a東北地區(qū)呈現(xiàn)暖干化趨勢,年平均氣溫和春玉米生長季內(nèi)平均氣溫每10a分別升高0.38℃和0.14℃,降水量分別減少8.5mm和8.0mm[26, 34]。IPCC第六次評估報告指出,未來幾十年氣候?qū)⒊掷m(xù)變暖[35]。對于東北地區(qū),未來氣候變暖使春玉米種植北界大范圍北移,中晚熟品種種植范圍北擴(kuò)至黑龍江省的最北部地區(qū)[36]。溫度和降水是影響玉米籽粒脫水最主要的因素,未來氣候持續(xù)變暖將進(jìn)一步影響東北地區(qū)春玉米機械粒收的氣候適宜性。因此,本研究基于未來氣候情景數(shù)據(jù)和當(dāng)?shù)刂髟灾型硎齑河衩灼贩N,分析未來氣候變化下機械粒收的氣候適宜區(qū),可為當(dāng)?shù)匚磥砗侠硗茝V應(yīng)用機械粒收技術(shù)提供參考,同時,基于東北地區(qū)未來情景數(shù)據(jù)與春玉米生育期數(shù)據(jù),通過預(yù)測不同春玉米品種初霜日前的籽粒含水率,確定了29個中晚熟春玉米品種2021-2060年機械粒收的氣候適宜區(qū),結(jié)果表明,不同春玉米中晚熟品種的機械粒收氣候適宜區(qū)空間分布差異較大,說明在相同的脫水環(huán)境下,東北地區(qū)春玉米機械粒收適宜性受春玉米品種脫水特性影響較大,這與王克如等[20]的結(jié)論一致。基于當(dāng)前遼河流域地區(qū)的主栽品種,在兩個氣候情景下未來40a春玉米機械粒收最適宜區(qū)和適宜區(qū)北界均主要位于遼寧省,總體氣候適宜性較低,說明東北地區(qū)近年多選種生育期較長、收獲前可用于籽粒脫水的時間較短的春玉米品種,這與黃兆福等[25]的研究結(jié)果一致,當(dāng)前遼河流域主栽春玉米品種含水率降至適宜機械粒收的積溫需求高于國審適宜機械粒收品種,在當(dāng)?shù)爻R?guī)收獲期無法實現(xiàn)高質(zhì)量機械粒收[34, 37]。因此,未來可通過選育和種植脫水速率快、生育期相對較短的春玉米品種,尤其是對于各個敏感帶,選種合適的品種可提高機械粒收的適宜性,利于機械粒收的高質(zhì)量應(yīng)用。

    玉米籽粒脫水速率在品種間存在由基因型控制的可遺傳差異,與苞葉薄厚和松緊、果穗大小和長短、穗軸粗細(xì)以及籽粒長度等植株性狀相關(guān),同時受空氣濕度、溫度和降水量等氣象因素以及種植密度、水肥管理和播期等栽培措施的影響[20]。本文采用趨勢含水率與氣象含水率分離的方式,主要基于氣象條件估算初霜日前春玉米的籽粒含水率,從而確定東北地區(qū)機械粒收氣候適宜區(qū),其中氣象含水率由前人在黃淮海地區(qū)確定的玉米籽粒含水率與氣象因子的回歸關(guān)系模型進(jìn)行初步計算,不同品種與氣象因子間相互作用的差異性以及栽培管理措施的影響可能導(dǎo)致估算的氣象含水率存在不確定性,因此,未來需融合更多的試驗數(shù)據(jù)進(jìn)一步確定不同生態(tài)區(qū)春玉米與氣象條件的定量關(guān)系,進(jìn)行更準(zhǔn)確的氣象含水率計算。采用積溫?含水率預(yù)測模型(Logistic power模型)模擬玉米籽粒趨勢含水率,可對趨勢含水率進(jìn)行較好的預(yù)測。前人研究表明,通過調(diào)整播期、增大玉米種植密度和控制氮肥量等可以在熱量限制區(qū)域加快玉米干物質(zhì)累積速率,減小生理成熟時玉米籽粒含水率[38?40]。統(tǒng)計模型無法較好地反映不同管理措施下籽粒的脫水差異,未來可通過結(jié)合數(shù)理模型與作物生長模型模擬的方式,在確定玉米籽粒脫水生理過程的基礎(chǔ)上預(yù)測不同管理措施下的玉米籽粒含水率,進(jìn)一步探索機械粒收的適宜性。此外,重點分析東北地區(qū)中晚熟品種未來機械粒收氣候適宜性,對于東北地區(qū)中南部地區(qū),采用當(dāng)前多個品種確定未來機械粒收適宜區(qū)具有較好的代表性,但對于黑龍江省等因未來氣候變暖、熱量增加而新增的中晚熟種植區(qū),熱量條件不能滿足籽粒脫水的熱量要求,可能導(dǎo)致低估了這些區(qū)域的機械粒收適宜性。例如德美亞系列品種具有抗性強、適應(yīng)性廣、耐密植、脫水快、產(chǎn)量高、品質(zhì)好并適宜機械化收獲的特點,引進(jìn)后在黑龍江省、吉林省等地得到迅速推廣[41]。未來可通過在該地區(qū)開展系列試驗研究工作,明確其籽粒脫水速率與氣象因子的關(guān)系,進(jìn)一步明確黑龍江省等地的機械粒收氣候適宜性并進(jìn)行氣候適宜分區(qū),為機械粒收在黑龍江省的推廣應(yīng)用提供參考。

    3.2 結(jié)論

    未來40a不同春玉米品種的氣候最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)均表現(xiàn)為由西南向東北延伸的空間分布特征,對于氣候適宜性較強的品種,不適宜區(qū)主要位于高海拔的吉林省東部山區(qū)(長白、東崗和松江一帶)和高緯度的黑龍江省北部地區(qū)。兩個氣候情景下,未來40a最適宜區(qū)北界范圍主要為遼寧省除清原、桓仁外的所有區(qū)域和吉林省西部;適宜區(qū)北界變動范圍位于吉林省除南部外的所有區(qū)域和吉林省西部以及黑龍江省西南部;次適宜區(qū)北界變動范圍位于遼寧省東北部、吉林省除東南部外的所有區(qū)域和黑龍江省。RCP4.5氣候情景下,50%和80%春玉米品種的最適宜區(qū)位于遼寧省南部,最北分別可達(dá)北緯42.0°和41.0°,最東可達(dá)東經(jīng)123.4°和123.3°,包含的春玉米種植面積分別占研究區(qū)域內(nèi)春玉米總種植面積的7.9%和2.2%;適宜區(qū)位于遼寧省北部與中部,最北可至北緯43.0°和42.7°,最東可至東經(jīng)124.7°和124.4°,包含的春玉米種植面積占比分別為11.2%和11.9%。RCP8.5氣候情景下,未來40a內(nèi)50%和80%品種的最適宜區(qū)遼寧省西部和南部,最北分別可達(dá)北緯42.8°和41.9°,最東可至東經(jīng)124.5°和123.3°,包含的春玉米種植面積占比分別為15.3%和7.7%;適宜區(qū)位于遼寧省北部和吉林省西南部,最北分別可達(dá)北緯44.8°和42.9°,最東可至東經(jīng)124.7°和124.5°,包含春玉米種植面積占比分別為6.1%和8.8%。兩個氣候情景下,氣候不適宜區(qū)和次適宜區(qū)位于黑龍江省和吉林省以及遼寧省東北部。

    [1] 中華人民共和國國家統(tǒng)計局.中國統(tǒng)計年鑒2021[M].北京:中國統(tǒng)計出版社,2022.

    National Bureau of Statistics of China.China statistical yearbook of 2021[M].Beijing:China Statistics Press,2022. (in Chinese)

    [2] 李少昆,王克如,謝瑞芝,等.機械粒收推動玉米生產(chǎn)方式轉(zhuǎn)型[J].中國農(nóng)業(yè)科學(xué),2018,51(10):1842-1844.

    Li S K,Wang K R,Xie R Z,et al.Grain mechanical harvesting technology promotes the transformation of maize production mode[J].Scientia Agricultura Sinica,2018, 51(10):1842-1844.(in Chinese)

    [3] 李少昆,王克如,謝瑞芝,等.實施密植高產(chǎn)機械化生產(chǎn)實現(xiàn)玉米高產(chǎn)高效協(xié)同[J].作物雜志,2016(4):1-6.

    Li S K,Wang K R,Xie R Z,et al.Implementing higher population and full mechanization technologies to achieve high yield and high efficiency in maize production[J]. Crops,2016(4):1-6.(in Chinese)

    [4] Waelti H.Physical properties and morphological characteristics of maize and their influence on threshing injury of kernels[D].Ames:Iowa State University,1967.

    [5] Waelti H,Buchele W F.Factors affecting corn kernel damage combine cylinders[J].Transactions of the ASAE, 1969,12(1):55-59.

    [6] Dutta P K.Effects of grain moisture,drying methods, and variety on breakage susceptibility of shelled corns as measured by the Wisconsin Breakage Tester[D].Ames: Iowa State University,1986.

    [7] 王克如,孔令杰,袁建華,等.江蘇沿海地區(qū)夏玉米機械粒收質(zhì)量與品種篩選研究[J].玉米科學(xué),2018,26(5):110-116.

    Wang K R,Kong L J,Yuan J H,et al.Grain quality of maize mechanical kernel harvest and varieties screen in Jiangsu coastal area[J].Journal of Maize Science,2018,26(5): 110- 116.(in Chinese)

    [8] 王克如,李少昆,王延波,等.遼寧中部適宜機械粒收玉米品種的篩選[J].作物雜志,2018b,(3):97-102.

    Wang K R,Li S K,Wang Y B,et al.Screening maize varieties suitable for mechanical harvesting grain in the central Liaoning province[J].Crops,2018b,(3):97-102.(in Chinese)

    [9] 李洪,閻曉光,王國梁,等.山西省春玉米機械粒收品種篩選與大田示范[J].玉米科學(xué),2018,26(5):117-122.

    Li H,Yan X G,Wang G L,et al.Cultivars selection of high yield spring maize suitable for mechanical kernel harvest and filed demonstration in Shanxi[J].Journal of Maize Science,2018,26(5):117-122.(in Chinese)

    [10] 李少昆,王克如,王立春,等.吉林玉米機械粒收質(zhì)量影響因素研究及品種篩選[J].玉米科學(xué),2018,26(4):55-62.

    Li S K,Wang K R,Wang L C,et al.Study on the influencing factors of maize grain quality of mechanically harvesting and corresponding cultivar selection in Jilin province[J]. Journal of Maize Science,2018,26(4):55- 62.(in Chinese)

    [11] 丁德芳,柴宗文,徐喜俊,等.河西灌區(qū)玉米密植高產(chǎn)機械粒收品種篩選[J].干旱地區(qū)農(nóng)業(yè)研究,2019,37(5):144-148.

    Ding D F,Chai Z W,Xu X J,et al.The selection of high yield maize cultivars suitable for dense planting and machine harvest in Hexi irrigation district[J].Agricultural Research in Arid Areas,2019,37(5):144-148.(in Chinese)

    [12] 孔凡磊,趙波,詹小旭,等.四川省夏玉米機械粒收適宜品種篩選與影響因素分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 2020,28(6): 835-842.

    Kong F L,Zhao B,Zhan X X,et al.Varieties screening of mechanical grain harvest and analysis of influencing factors of summer maize in Sichuan province[J].Chinese Journal of Eco-Agriculture(Chinese and English),2020, 28(6):835-842.(in Chinese)

    [13] 董朋飛,王群,王振華,等.河南省夏玉米品種篩選與機械粒收質(zhì)量研究[J].玉米科學(xué),2020,28(2):128-134.

    Dong P F,Wang Q,Wang Z H,et al.Grain quality of maize mechanical kernel harvest and varieties screen in Henan province[J].Journal of Maize Sciences,2020,28(2):128-134. (in Chinese)

    [14] 趙波,吳雅薇,李小龍,等.玉米強弱勢粒間機械脫粒破碎率的差異[J].中國生態(tài)農(nóng)業(yè)學(xué)報(中英文),2020,28(6): 843-851.

    Zhao B,Wu Y W, Li X L,et al.Differences in mechanical threshing broken rate between superior and inferior maize grains[J].Chinese Journal of Eco-Agriculture(Chinese and English),2020,28(6):843-851.(in Chinese)

    [15] 閻曉光,李洪,董紅芬,等.山西省春玉米機械粒收質(zhì)量及其影響因素分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報(中英文),2021, 29(10):1704-1711.

    Yan X G,Li H,Dong H F,et al.Analysis of mechanical grain harvesting quality and its influential factors of spring maize in Shanxi province[J].Chinese Journal of Eco-Agriculture, 2021,29(10):1704-1711.(in Chinese)

    [16] 謝瑞芝,雷曉鵬,王克如,等.黃淮海夏玉米子粒機械收獲研究初報[J].作物雜志,2014(2):76-79.

    Xie R Z,Lei X P,Wang K R,et al.Research on corn mechanically harvesting grain quality in Huanghuaihai plain[J].Crops,2014(2): 76-79.(in Chinese)

    [17] 李少昆.我國玉米機械粒收質(zhì)量影響因素及粒收技術(shù)的發(fā)展方向[J].石河子大學(xué)學(xué)報(自然科學(xué)版),2017,35(3): 265-272.

    Li S K.Factors affecting the quality of maize grain mechanical harvest and the development trend of grain harvest technology[J].Journal of Shihezi University (Natural Science Edition),2017,35(3):265-272.(in Chinese)

    [18] 李璐璐,雷曉鵬,謝瑞芝,等.夏玉米機械粒收質(zhì)量影響因素分析[J].中國農(nóng)業(yè)科學(xué),2017,50(11):2044-2051.

    Li L L,Lei X P,Xie R Z,et al.Analysis of influential factors on mechanical grain harvest quality of summer maize[J]. Scientia Agricultura Sinica,2017,50(11):2044-2051.(in Chinese)

    [19] 柴宗文,王克如,郭銀巧,等.玉米機械粒收質(zhì)量現(xiàn)狀及其與含水率的關(guān)系[J].中國農(nóng)業(yè)科學(xué),2017,50(11):2036- 2043.

    Chai Z W,Wang K R,Guo Y Q,et al.Current status of maize mechanical grain harvesting and its relationship with grain moisture content[J].Scientia Agricultura Sinica,2017,50 (11):2036-2043.(in Chinese)

    [20] 王克如,李少昆.玉米籽粒脫水速率影響因素分析[J].中國農(nóng)業(yè)科學(xué),2017,50(11):2027-2035.

    Wang K R,Li S K.Analysis of factors affecting corn kernel dehydration rate[J].Scientia Agricultura Sinica,2017,50(11): 2027-2035.(in Chinese)

    [21] 李璐璐,明博,高尚,等.夏玉米籽粒脫水特性及與灌漿特性的關(guān)系[J].中國農(nóng)業(yè)科學(xué),2018,51(10):1878-1889.

    Li L L,Ming B,Gao S,et al.Study on grain dehydration characters of summer maize and its relationship with grain filling[J].Scientia Agricultura Sinica,2018,51(10):1878-1889. (in Chinese)

    [22] 李璐璐,明博,謝瑞芝,等.黃淮海夏玉米品種脫水類型與機械粒收時間的確立[J].作物學(xué)報,2018,44(12):1764-1773.

    Li L L,Ming B,Xie R Z,et al.Grain dehydration types and establishment of mechanical grain harvesting time for summer maize in the Yellow-Huai-Hai river plain[J].Acta Agronomica Sinica,2018,44(12):1764-1773.(in Chinese)

    [23] 郭亞南,明博,王克如,等.北京地區(qū)玉米子粒脫水特性與適宜收獲期研究[J].玉米科學(xué),2019,27(5):130-136.

    Guo Y N,Ming B,Wang K R,et al.Maize grain dehydration characteristics and suitable harvest time in Beijing region [J].Journal of Maize Science,2019,27(5):130-136.(in Chinese)

    [24] 李紅燕.寧夏玉米籽粒脫水特性與機械粒收時期預(yù)測[D].石河子:石河子大學(xué),2019.

    Li H Y.Dehydration characteristics of maize kernels in Ningxia and prediction of mechanical grain harvesting period[D].Shihezi: Shihezi University,2019.(in Chinese)

    [25] 黃兆福,明博,王克如,等.遼河流域玉米籽粒脫水特點及適宜收獲期分析[J].作物學(xué)報,2019,45(6):922-931.

    Huang Z F,Ming B,Wang K R,et al.Characteristics of maize grain dehydration and prediction of suitable harvest period in the Liao River Basin[J].Acta Agronomica Sinica, 2019,45(6):922-931.(in Chinese)

    [26] 劉志娟,楊曉光,王文峰,等.氣候變化背景下我國東北三省農(nóng)業(yè)氣候資源變化特征[J].應(yīng)用生態(tài)學(xué)報,2009,20(9): 2199-2206.

    Liu Z J,Yang X G,Wang W F,et al.Characteristics of agricultural climate resources in three provinces of Northeast China under global climate change[J].Chinese Journal of Applied Ecology,2009,20(9):2199-2206.(in Chinese)

    [27] 李克南,楊曉光,劉志娟,等.全球氣候變化對中國種植制度可能影響分析Ⅲ.中國北方地區(qū)氣候資源變化特征及其對種植制度界限的可能影響[J].中國農(nóng)業(yè)科學(xué),2010, 43(10):2088-2097.

    Li K N,Yang X G,Liu Z J,et al.Analysis of the potential influence of global climate change on cropping systems in China Ⅲ.The change characteristics of climatic resources in northern China and its potential influence on cropping system[J].Scientia Agricultura Sinica,2010,43(10):2088- 2097.(in Chinese)

    [28] Luo Y C,Zhang Z,Li Z Y,et al.Data for:Identifying the spatiotemporal changes of annual harvesting areas for three staple crops in China by integrating multi-data sources[J]. Mendeley Data,2020,V2.doi:10.17632/jbs44b2hrk.2.

    [29] 高尚,明博,李璐璐,等.黃淮海夏玉米籽粒脫水與氣象因子的關(guān)系[J].作物學(xué)報,2018,44(12):1755-1763.

    Gao S,Ming B,Li L L,et al.Relationship between grain dehydration and meteorological factors in the Yellow- Huai-Hai rivers summer maize[J].Acta Agronomica Sinica, 2018,44(12):1755-1763.(in Chinese)

    [30] Brooking Ian R.Maize ear moisture during grain-filling, and its relation to physiological maturity and grain- drying[J].Elsevier,1990,23(1), 55-68.

    [31] 馮東升,高樹仁,楊克軍.解析玉米籽粒脫水的動力問題[J].中國種業(yè),2017(10):35-36.

    Feng D S,Gao S R,Yang K J.Analysis of the dynamic problem of corn kernel dehydration[J].China Seed Industry, 2017(10):35-36. (in Chinese)

    [32] Allen R G,Pereira L S,Raes D,et al.Crop Evapotranspiration: guidelines for computing crop water requirements[J].FAO Irrigation and Drainage Paper 56.1998.

    [33] 寧曉菊,張麗君,楊群濤,等.1951 年以來中國無霜期的變化趨勢[J].地理學(xué)報,2015,70(11):1811-1822.

    Ning X J,Zhang L J,Yang Q T,et al.Trends in the frost-free period in China from 1951 to 2012[J].Acta Geographica Sinica,2015,70(11):1811-1822.(in Chinese)

    [34] Liu Z J,Yang X G,Chen F,et al.The effects of past climate change on the northern limits of maize planting in Northeast China[J].Climatic Change,2013,117(4):891-902.

    [35] Allen M R,Johansen T G,Poloczanska E,et al.Global warming of 1.5°C.An IPCC special report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways,in the context of strengthening the global response to the threat of climate change,sustainable development,and efforts to eradicate poverty[R].Intergovernmental Panel on Climate Chang(IPCC),2018.

    [36] 劉志娟,楊曉光,王文峰,等.全球氣候變暖對中國種植制度可能影響Ⅳ.未來氣候變暖對東北地區(qū)春玉米種植北界的可能影響[J].中國農(nóng)業(yè)科學(xué),2010,43(11):2280-2291.

    Liu Z J,Yang X G,Wang W F,et al.The possible effects of global warming on cropping systems in China Ⅳ.The possible impact of future climatic warming on the northern limits of spring maize in three provinces of northeast China[J].Scientia Agricultura Sinica,2010,43(11):2280-2291. (in Chinese)

    [37] 淮賀舉,李存軍,李奇峰,等.東北地區(qū)春玉米物候期時空動態(tài)變化特征[J].玉米科學(xué),2019,27(1):81-86.

    Huai H J,Li C J,Li Q F,et al.Changes of spring maize phenology from 1981 to 2012 in the Northeast Plain of China[J].Journal of Maize Sciences,2019,27(1):81-86.(in Chinese)

    [38] 付江鵬,賀正,賈彪,等.水肥一體化施氮水平對玉米籽粒灌漿和脫水過程的影響[J].中國農(nóng)業(yè)氣象,2019,40(12): 772-782.

    Fu J P,He Z,Jia B,et al.Effect of integrated fertilization level of water and fertilizer on maize grain filling and dehydration process[J].Chinese Journal of Agrometeorology, 2019,40(12):772-782.(in Chinese)

    [39] 鄒俊麗,崔兆韻,徐祎,等.黃淮海地區(qū)不同播期夏玉米籽粒脫水過程和產(chǎn)量形成的比較[J].中國農(nóng)業(yè)氣象, 2022,43(2):112-123.

    Zou J L,Cui Z Y,Xu Y,et al.Comparison of grain dehydration process and yield formation for summer maize under different sowing dates in Huanghuaihai Plain[J]. Chinese Journal of Agrometeorology,2022,43(2):112-123. (in Chinese)

    [40] 劉夢,葛均筑,吳錫冬,等.華北熱資源限制區(qū)實現(xiàn)夏玉米高產(chǎn)高效和機械粒收的適宜種植密度和施氮量[J].植物營養(yǎng)與肥料學(xué)報,2022,28(8):1444-1456.

    Liu M,Ge J Z,Wu X D,et al.Optimizing planting density and nitrogen application of summer maize to achieve high yield and efficiency,and mechanical grain harvesting in the thermal resource-limited area of North China[J].Journal of Plant Nutrition and Fertilizers,2022,28(8):1444-1456.(in Chinese)

    [41] 寧毅,安英輝,張志輝.德美亞系列玉米品種栽培技術(shù)[J].種子世界,2016(5):46-47.

    Ning Y,An Y H,Zhang Z H.Cultivation techniques of Demeiya series corn varieties[J].Seed World,2016(5):46- 47.(in Chinese)

    Climate Suitable Zones for Mechanical Kernel Harvesting of Middle-late Maturing Spring Maize in Northeast China under Future Climate Scenarios

    SU Zheng-e, LIU Zhi-juan, YANG Wan-rong, ZHU Guang-xin, SHI Deng-yu, YANG Xiao-guang

    (College of Resources and Environment Sciences, China Agricultural University, Beijing 100193, China)

    Fully mechanized maize production has been trending in China, and using agricultural machinery to harvest maize grains is the focus of the transformation of maize production methods. Grain moisture content is the main factor affecting mechanical kernel harvesting, which is mainly affected by the agronomic characteristics of maize, ecological meteorological factors at the late growth stage, cultivation and management practice and other factors. Northeast China (NEC) is one of the main maize production regions and is the most sensitive regions to climate change in China. Clarifying the climate suitability and suitable zones for mechanical maize kernel harvesting in NEC under future climate conditions can provide scientific basis for improving the mechanization level of maize production and realizing the fully mechanized maize production in the future. In this study, based on the predicted climate data of 2021?2060 under two greenhouse gas concentration scenarios of RCP4.5 and RCP8.5 (the radiation forcing levels will reach 4.5W·m?2and 8.5W·m?2by 2100, which representing medium and high emission scenario, respectively) and the growth period data of spring maize, the northern boundaries of climate suitable zones for mechanical kernel harvesting of currently widely planted 29 middle-late spring maize varieties were determined using the kernel moisture content prediction models, and further the multi-variety based climate suitable zones were analyzed for the future 40 years. The results showed that climate very suitable, suitable, moderately suitable and marginally suitable zones for mechanical kernel harvesting will extend from southwest to northeast in NEC. And among different varieties, the northern boundaries of climate suitable zones fluctuate from the southwest to the southeastern mountainous areas of Jilin province (Changbai, Donggang and Songjiang areas) and the northern high latitude areas of Heilongjiang province. Under RCP4.5 (RCP8.5) climate scenario, the very suitable zones for 50% and 80% of middle-late maturing spring maize varieties located in the southern part of Liaoning province (the western and southern parts of Liaoning province). The northernmost can reach 42.0°N and 41.0°N (42.8°N and 41.9°N) latitude, and the easternmost can reach 123.4°E and 123.3°E (124.5°E and 123.3°E) longitude, respectively. The areas account for 5.9% and 1.8% (11.2% and 5.9%) of the potential growing area of spring maize in NEC, respectively. Suitable zones located in the northern and central parts of Liaoning province (the northern part of Liaoning province and southwestern part of Jilin province), with the northernmost reaching 43.0°N and 42.7°N (44.8°N and 42.9°N) latitude, the easternmost reaching 124.7°E and 124.4°E (124.7°E and 124.5°E) longitude, and the areas accounting for 8.3% and 8.9% (4.7% and 6.6%) of potential growing area of spring maize, respectively. In the future 40 years, compared with the RCP4.5 climate scenario, the middle-late maturing varieties of spring maize under RCP8.5 climate scenario have stronger climatic suitability for mechanical kernel harvesting, but the climate very suitable and suitable zones under the two climate scenarios are both small. The climate suitable zones for mechanical kernel harvesting will differ greatly between varieties, so varieties with short growth periods and fast dehydration rates which are suitable for mechanical kernel harvesting should be selected to improve the quality and efficiency of mechanical kernel harvesting in the future.

    Spring maize; Kernel dehydration; Mechanical kernel harvesting; Climate suitable zones; Climate change; Northeast China

    10.3969/j.issn.1000-6362.2023.08.001

    蘇正娥,劉志娟,楊婉蓉,等.未來氣候變化情景下東北地區(qū)中晚熟春玉米機械粒收氣候適宜區(qū)分析[J].中國農(nóng)業(yè)氣象,2023, 44(8):649-663

    2022?09?23

    國家重點研發(fā)計劃項目(2019YFA0607402)

    劉志娟,副教授,博士,博士生導(dǎo)師,主要從事氣候變化對農(nóng)業(yè)影響與適應(yīng)研究,E-mail:zhijuanliu@cau.edu.cn

    蘇正娥,E-mail:suzhenge_cau@163.com

    猜你喜歡
    吉林省籽粒遼寧省
    吉林省教育廳新年賀詞
    籽粒莧的飼用價值和高產(chǎn)栽培技術(shù)
    籽粒莧的特性和種植技術(shù)
    吉林省“十四五”食品安全信息化建設(shè)的幾點思考
    吉林省梅河口老年大學(xué)之歌
    吉林省完成1.4萬公里農(nóng)村公路“暢返不暢”整治
    石油瀝青(2019年6期)2019-02-13 04:24:34
    玉米機械脫粒籽粒含水量與破碎率的相關(guān)研究
    遼寧省各級工會收看十九大開幕會
    遼寧省環(huán)境保護(hù)廳與遼寧省氣象局進(jìn)一步深化合作
    遼寧省自然環(huán)境分類探索
    亚洲国产毛片av蜜桃av| 91午夜精品亚洲一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 一区二区三区乱码不卡18| 22中文网久久字幕| 又黄又爽又刺激的免费视频.| 简卡轻食公司| 欧美性感艳星| 午夜日本视频在线| 人妻制服诱惑在线中文字幕| 新久久久久国产一级毛片| 在线 av 中文字幕| 亚洲欧美日韩卡通动漫| 亚洲精品国产色婷婷电影| 成人手机av| 国产av码专区亚洲av| 国产色婷婷99| 久久久久久人妻| 美女国产高潮福利片在线看| 亚洲色图综合在线观看| 国产一区亚洲一区在线观看| 国产欧美日韩综合在线一区二区| 国产精品国产三级专区第一集| 你懂的网址亚洲精品在线观看| 免费看光身美女| 18+在线观看网站| 国产精品麻豆人妻色哟哟久久| 麻豆乱淫一区二区| 日本欧美国产在线视频| 王馨瑶露胸无遮挡在线观看| 国产亚洲最大av| 亚洲综合色惰| 哪个播放器可以免费观看大片| 国产精品国产三级专区第一集| 成人毛片60女人毛片免费| 女人久久www免费人成看片| 国产不卡av网站在线观看| 97在线人人人人妻| 国产爽快片一区二区三区| 丰满饥渴人妻一区二区三| 亚洲欧美中文字幕日韩二区| 国产精品熟女久久久久浪| 国产高清有码在线观看视频| 有码 亚洲区| 国产一区二区在线观看av| 日本色播在线视频| 亚洲精品自拍成人| 亚洲国产精品一区二区三区在线| 高清视频免费观看一区二区| 成人午夜精彩视频在线观看| 久久精品国产亚洲av天美| 99热6这里只有精品| 夜夜骑夜夜射夜夜干| 国产成人a∨麻豆精品| 成人18禁高潮啪啪吃奶动态图 | 精品熟女少妇av免费看| 纯流量卡能插随身wifi吗| 亚洲国产色片| 99九九在线精品视频| 男女边吃奶边做爰视频| 日本av免费视频播放| 插逼视频在线观看| 不卡视频在线观看欧美| 亚洲人成网站在线观看播放| 黑人巨大精品欧美一区二区蜜桃 | 亚洲,一卡二卡三卡| 美女国产视频在线观看| 欧美亚洲 丝袜 人妻 在线| 波野结衣二区三区在线| 午夜福利,免费看| 日韩一区二区视频免费看| 亚洲人与动物交配视频| 国产免费一区二区三区四区乱码| 亚洲欧洲日产国产| 日本vs欧美在线观看视频| 一级,二级,三级黄色视频| 欧美三级亚洲精品| 久热这里只有精品99| 精品视频人人做人人爽| 在线亚洲精品国产二区图片欧美 | 欧美最新免费一区二区三区| 国产av码专区亚洲av| 精品亚洲成国产av| 国产精品久久久久久精品古装| 校园人妻丝袜中文字幕| 欧美日韩精品成人综合77777| 国产精品 国内视频| 久久国产精品男人的天堂亚洲 | 777米奇影视久久| 免费高清在线观看视频在线观看| 日韩av在线免费看完整版不卡| 狠狠婷婷综合久久久久久88av| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 精品视频人人做人人爽| 青春草视频在线免费观看| 国产在线一区二区三区精| 欧美 日韩 精品 国产| 天堂俺去俺来也www色官网| 国产成人aa在线观看| 精品久久久精品久久久| 国模一区二区三区四区视频| 肉色欧美久久久久久久蜜桃| 国产淫语在线视频| 国产亚洲av片在线观看秒播厂| 18禁观看日本| av线在线观看网站| 国产深夜福利视频在线观看| 亚洲精品456在线播放app| 少妇 在线观看| 亚洲精品第二区| 日韩免费高清中文字幕av| 在线观看美女被高潮喷水网站| 日本与韩国留学比较| 欧美bdsm另类| 国产 一区精品| 国产精品成人在线| 日韩av不卡免费在线播放| 制服丝袜香蕉在线| 三级国产精品欧美在线观看| 精品熟女少妇av免费看| 99视频精品全部免费 在线| 色94色欧美一区二区| 黄色配什么色好看| 久久久国产精品麻豆| 91精品三级在线观看| 伦精品一区二区三区| 97超碰精品成人国产| 亚洲丝袜综合中文字幕| 久热这里只有精品99| 青春草视频在线免费观看| 中文字幕久久专区| av专区在线播放| 亚洲性久久影院| 欧美+日韩+精品| 99国产精品免费福利视频| 免费看光身美女| 国产色爽女视频免费观看| 五月伊人婷婷丁香| 亚洲精品久久久久久婷婷小说| 岛国毛片在线播放| 国产精品.久久久| 国产色爽女视频免费观看| 国产在线一区二区三区精| 高清在线视频一区二区三区| 色5月婷婷丁香| 人妻 亚洲 视频| 大又大粗又爽又黄少妇毛片口| 男女无遮挡免费网站观看| 国产淫语在线视频| 99久久精品一区二区三区| 人妻 亚洲 视频| 成人综合一区亚洲| 国产精品人妻久久久久久| 两个人免费观看高清视频| 另类亚洲欧美激情| 久久久精品免费免费高清| 亚洲精品久久成人aⅴ小说 | 尾随美女入室| 亚洲精华国产精华液的使用体验| 亚洲欧美成人精品一区二区| 国产精品不卡视频一区二区| 国产毛片在线视频| 九草在线视频观看| 91精品三级在线观看| 久久鲁丝午夜福利片| 看免费成人av毛片| 大香蕉久久网| 看十八女毛片水多多多| 黑丝袜美女国产一区| 欧美精品人与动牲交sv欧美| 久久久午夜欧美精品| 美女大奶头黄色视频| 三级国产精品欧美在线观看| 伊人亚洲综合成人网| 日韩av不卡免费在线播放| 久久久午夜欧美精品| 街头女战士在线观看网站| 欧美日韩在线观看h| 最新中文字幕久久久久| 九九久久精品国产亚洲av麻豆| 亚洲国产av影院在线观看| 免费大片黄手机在线观看| 欧美变态另类bdsm刘玥| 国产精品国产av在线观看| 伊人久久精品亚洲午夜| 国产在线免费精品| 夫妻性生交免费视频一级片| 寂寞人妻少妇视频99o| 母亲3免费完整高清在线观看 | 亚洲欧洲国产日韩| 亚洲美女视频黄频| videos熟女内射| 亚洲精华国产精华液的使用体验| 久久久久久久久久久免费av| 超色免费av| 国产黄频视频在线观看| 久久人人爽人人爽人人片va| av在线播放精品| 大香蕉97超碰在线| 啦啦啦啦在线视频资源| 在线观看一区二区三区激情| 日韩,欧美,国产一区二区三区| 边亲边吃奶的免费视频| 亚洲怡红院男人天堂| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 飞空精品影院首页| 大又大粗又爽又黄少妇毛片口| 亚洲无线观看免费| 欧美亚洲 丝袜 人妻 在线| 国产高清有码在线观看视频| 欧美97在线视频| 亚洲精品久久久久久婷婷小说| 国产av国产精品国产| 一级爰片在线观看| 欧美日韩精品成人综合77777| 欧美日韩国产mv在线观看视频| 成人综合一区亚洲| 婷婷色综合大香蕉| 好男人视频免费观看在线| 午夜免费观看性视频| 成人手机av| 大片免费播放器 马上看| 国产色婷婷99| 亚洲天堂av无毛| 亚洲高清免费不卡视频| 欧美一级a爱片免费观看看| 这个男人来自地球电影免费观看 | 久久国产精品大桥未久av| 2018国产大陆天天弄谢| 久久久久久久国产电影| 纵有疾风起免费观看全集完整版| 免费av中文字幕在线| 国产精品一区www在线观看| 蜜桃国产av成人99| 成人影院久久| 啦啦啦啦在线视频资源| 午夜视频国产福利| 精品熟女少妇av免费看| 久久久国产一区二区| 日韩亚洲欧美综合| 日本欧美视频一区| av网站免费在线观看视频| 亚洲综合色惰| 我要看黄色一级片免费的| 中文字幕久久专区| 永久网站在线| 欧美另类一区| 建设人人有责人人尽责人人享有的| 91精品三级在线观看| 69精品国产乱码久久久| 国产av国产精品国产| tube8黄色片| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 精品亚洲成a人片在线观看| 久久99一区二区三区| 国国产精品蜜臀av免费| 欧美精品亚洲一区二区| 欧美+日韩+精品| 97精品久久久久久久久久精品| 午夜老司机福利剧场| 成人综合一区亚洲| 蜜臀久久99精品久久宅男| 欧美国产精品一级二级三级| 老司机影院成人| 高清视频免费观看一区二区| 五月伊人婷婷丁香| 热re99久久国产66热| 麻豆精品久久久久久蜜桃| 性高湖久久久久久久久免费观看| www.色视频.com| 最近中文字幕高清免费大全6| 最后的刺客免费高清国语| 久久人人爽av亚洲精品天堂| 王馨瑶露胸无遮挡在线观看| 欧美丝袜亚洲另类| 欧美激情 高清一区二区三区| 制服丝袜香蕉在线| 国产精品无大码| 午夜日本视频在线| 免费久久久久久久精品成人欧美视频 | 熟妇人妻不卡中文字幕| 精品99又大又爽又粗少妇毛片| 三级国产精品片| 免费日韩欧美在线观看| 丝袜在线中文字幕| 丝袜喷水一区| 久久鲁丝午夜福利片| 女人久久www免费人成看片| 人妻夜夜爽99麻豆av| 老熟女久久久| 久久久国产欧美日韩av| 18+在线观看网站| 99re6热这里在线精品视频| 国产免费一级a男人的天堂| 国产精品国产三级专区第一集| 最黄视频免费看| 精品国产一区二区久久| 18禁在线播放成人免费| 人妻人人澡人人爽人人| 啦啦啦啦在线视频资源| 久久久久久久久久久丰满| 热99久久久久精品小说推荐| 国产精品一区二区三区四区免费观看| 九草在线视频观看| 中文天堂在线官网| 国产 一区精品| 在线天堂最新版资源| 久久99蜜桃精品久久| 如日韩欧美国产精品一区二区三区 | 丝袜在线中文字幕| 高清毛片免费看| 久久久久久人妻| 人妻一区二区av| 欧美日韩精品成人综合77777| 国产成人a∨麻豆精品| 日本91视频免费播放| 国产欧美亚洲国产| av一本久久久久| 飞空精品影院首页| 亚洲av日韩在线播放| 国产免费一区二区三区四区乱码| 午夜视频国产福利| 精品国产一区二区久久| 欧美亚洲日本最大视频资源| 永久网站在线| 久久久久人妻精品一区果冻| 精品一区二区三卡| 欧美日韩精品成人综合77777| 老司机影院成人| 国产有黄有色有爽视频| 久久影院123| 亚洲av不卡在线观看| 欧美激情国产日韩精品一区| 男女高潮啪啪啪动态图| 丝袜在线中文字幕| 在线看a的网站| 性色avwww在线观看| 精品人妻在线不人妻| 欧美精品一区二区免费开放| 国产精品久久久久成人av| 国产av精品麻豆| 黑人猛操日本美女一级片| 中文字幕免费在线视频6| 九九在线视频观看精品| 国精品久久久久久国模美| 在线观看一区二区三区激情| 哪个播放器可以免费观看大片| 好男人视频免费观看在线| 岛国毛片在线播放| 多毛熟女@视频| 蜜桃久久精品国产亚洲av| 欧美日韩视频精品一区| 五月天丁香电影| 国产在线视频一区二区| 又大又黄又爽视频免费| 国产精品一区www在线观看| 99re6热这里在线精品视频| 一区二区三区精品91| 国产老妇伦熟女老妇高清| 亚洲在久久综合| 国产男人的电影天堂91| 免费黄网站久久成人精品| 欧美精品一区二区大全| 国产成人a∨麻豆精品| 国产精品麻豆人妻色哟哟久久| 日产精品乱码卡一卡2卡三| 免费观看无遮挡的男女| 性高湖久久久久久久久免费观看| 人人澡人人妻人| 在线观看人妻少妇| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 欧美日韩一区二区视频在线观看视频在线| 毛片一级片免费看久久久久| kizo精华| 亚州av有码| 老女人水多毛片| 亚洲av中文av极速乱| 99久久综合免费| 一区二区日韩欧美中文字幕 | 伦理电影大哥的女人| 午夜激情久久久久久久| 又黄又爽又刺激的免费视频.| 国产精品国产三级专区第一集| 日韩免费高清中文字幕av| 欧美日韩综合久久久久久| 中文字幕最新亚洲高清| 又黄又爽又刺激的免费视频.| 成人毛片60女人毛片免费| 久久午夜综合久久蜜桃| 久久久精品区二区三区| 国产 一区精品| 老司机亚洲免费影院| 国产av码专区亚洲av| 欧美性感艳星| 一级爰片在线观看| 亚洲怡红院男人天堂| 久久精品国产亚洲av天美| 亚洲av中文av极速乱| 婷婷色综合www| 制服丝袜香蕉在线| 国产成人a∨麻豆精品| 欧美老熟妇乱子伦牲交| 人人妻人人添人人爽欧美一区卜| 国产色婷婷99| 菩萨蛮人人尽说江南好唐韦庄| 新久久久久国产一级毛片| 亚洲av国产av综合av卡| 大香蕉97超碰在线| 成年人免费黄色播放视频| 国产av国产精品国产| 亚洲精品久久成人aⅴ小说 | 久久ye,这里只有精品| 久热这里只有精品99| 精品国产一区二区久久| 日本黄色片子视频| 在线观看免费日韩欧美大片 | 丝瓜视频免费看黄片| 久久久久精品久久久久真实原创| 亚洲精品第二区| 成年美女黄网站色视频大全免费 | 欧美日韩一区二区视频在线观看视频在线| 亚洲欧洲精品一区二区精品久久久 | 色视频在线一区二区三区| av在线老鸭窝| 日韩不卡一区二区三区视频在线| videos熟女内射| 久久精品国产自在天天线| 人妻系列 视频| 国产熟女午夜一区二区三区 | 两个人免费观看高清视频| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 久久久久精品久久久久真实原创| 久久久a久久爽久久v久久| 日本av免费视频播放| 十分钟在线观看高清视频www| av不卡在线播放| 日韩成人伦理影院| 我要看黄色一级片免费的| 国产精品一区二区三区四区免费观看| 伦理电影免费视频| 国产日韩欧美亚洲二区| 亚洲国产欧美在线一区| 亚洲国产精品国产精品| 日产精品乱码卡一卡2卡三| 飞空精品影院首页| 日韩一本色道免费dvd| 国产成人免费无遮挡视频| 97超碰精品成人国产| av又黄又爽大尺度在线免费看| 简卡轻食公司| 欧美日韩视频精品一区| 久久99热6这里只有精品| 26uuu在线亚洲综合色| 51国产日韩欧美| 人妻夜夜爽99麻豆av| 久久精品国产亚洲网站| 午夜福利网站1000一区二区三区| 我要看黄色一级片免费的| 日韩av不卡免费在线播放| 99久久精品国产国产毛片| 少妇的逼水好多| 夫妻午夜视频| 欧美97在线视频| 亚洲av电影在线观看一区二区三区| 一级爰片在线观看| 久久影院123| 男女无遮挡免费网站观看| 丰满迷人的少妇在线观看| 日本欧美视频一区| 日韩人妻高清精品专区| 国产视频内射| 亚洲av欧美aⅴ国产| 国产视频内射| 日韩人妻高清精品专区| 国产精品熟女久久久久浪| 日韩av不卡免费在线播放| 日韩免费高清中文字幕av| 伦精品一区二区三区| 成人影院久久| 国产精品人妻久久久久久| 一二三四中文在线观看免费高清| 亚洲国产日韩一区二区| 99久国产av精品国产电影| 国产高清不卡午夜福利| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 精品国产国语对白av| 国产精品99久久久久久久久| 能在线免费看毛片的网站| 亚洲成色77777| 亚洲国产色片| 天堂中文最新版在线下载| 一本久久精品| 亚洲av中文av极速乱| 亚洲人与动物交配视频| 免费观看av网站的网址| 国产亚洲精品久久久com| 国产有黄有色有爽视频| 亚洲欧美日韩另类电影网站| 交换朋友夫妻互换小说| 久久久久久久精品精品| 永久网站在线| 色婷婷久久久亚洲欧美| 91aial.com中文字幕在线观看| 少妇 在线观看| 国产欧美亚洲国产| 国产片特级美女逼逼视频| 亚洲欧美精品自产自拍| 国产精品国产三级国产专区5o| 最新的欧美精品一区二区| 欧美日韩av久久| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 天天躁夜夜躁狠狠久久av| 国产 精品1| 久久久久久久久久久久大奶| 男人爽女人下面视频在线观看| 免费看光身美女| 日日爽夜夜爽网站| 国产白丝娇喘喷水9色精品| 丝袜脚勾引网站| 我的老师免费观看完整版| 亚洲精品久久午夜乱码| 男女免费视频国产| 亚洲国产精品成人久久小说| 欧美bdsm另类| 久久久久久久大尺度免费视频| 啦啦啦在线观看免费高清www| videos熟女内射| 一区二区三区精品91| 国产亚洲一区二区精品| 69精品国产乱码久久久| 中文字幕av电影在线播放| 天天操日日干夜夜撸| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 韩国高清视频一区二区三区| 2022亚洲国产成人精品| 丝袜在线中文字幕| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 国产精品人妻久久久影院| 国产色婷婷99| 免费观看无遮挡的男女| 国产欧美亚洲国产| 亚洲情色 制服丝袜| 久久久久国产网址| 校园人妻丝袜中文字幕| 永久免费av网站大全| 男人爽女人下面视频在线观看| 精品国产乱码久久久久久小说| 国产无遮挡羞羞视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 美女中出高潮动态图| 午夜激情久久久久久久| 女性生殖器流出的白浆| av国产久精品久网站免费入址| 亚洲激情五月婷婷啪啪| 一级黄片播放器| 在线看a的网站| 免费黄色在线免费观看| 亚洲av成人精品一二三区| 亚洲少妇的诱惑av| 制服诱惑二区| 久久婷婷青草| av天堂久久9| 欧美日韩一区二区视频在线观看视频在线| 一本一本综合久久| 国产极品粉嫩免费观看在线 | 91精品三级在线观看| 久久国内精品自在自线图片| 国产精品久久久久久久电影| 亚洲欧美清纯卡通| 日本色播在线视频| 欧美 日韩 精品 国产| 男女边吃奶边做爰视频| 国产午夜精品一二区理论片| 一本久久精品| 大话2 男鬼变身卡| 建设人人有责人人尽责人人享有的| 精品久久久久久久久亚洲| 国产精品蜜桃在线观看| 男女国产视频网站| 国产精品99久久久久久久久| 美女国产高潮福利片在线看| 中文天堂在线官网| av黄色大香蕉| 国产午夜精品久久久久久一区二区三区| 18禁动态无遮挡网站| 日韩不卡一区二区三区视频在线| 乱人伦中国视频| av免费在线看不卡| 国产精品久久久久久精品古装| 有码 亚洲区| 尾随美女入室| 亚洲精品日本国产第一区| 两个人的视频大全免费| av线在线观看网站| av又黄又爽大尺度在线免费看| 性色av一级| 国产极品粉嫩免费观看在线 | a级片在线免费高清观看视频| 久久午夜综合久久蜜桃| 满18在线观看网站| 国产黄色视频一区二区在线观看| 99久久中文字幕三级久久日本| 精品久久久噜噜| 亚洲婷婷狠狠爱综合网| 高清午夜精品一区二区三区| 在线亚洲精品国产二区图片欧美 |