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

    氣候變化情景下黃淮海區(qū)域熱量資源及夏玉米溫度適宜度

    2016-06-22 02:52:06李萌申雙和呂厚荃韓艷褚榮浩沙修竹
    大氣科學(xué)學(xué)報(bào) 2016年3期
    關(guān)鍵詞:黃淮海積溫夏玉米

    李萌,申雙和*,呂厚荃,韓艷,褚榮浩,沙修竹

    ① 南京信息工程大學(xué) 應(yīng)用氣象學(xué)院,江蘇 南京 210044;② 南京信息工程大學(xué) 氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室/氣候與環(huán)境變化國際合作聯(lián)合實(shí)驗(yàn)室/氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心/江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210044;

    氣候變化情景下黃淮海區(qū)域熱量資源及夏玉米溫度適宜度

    李萌①②,申雙和①②*,呂厚荃③,韓艷④,褚榮浩①,沙修竹⑤

    ① 南京信息工程大學(xué) 應(yīng)用氣象學(xué)院,江蘇 南京 210044;② 南京信息工程大學(xué) 氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室/氣候與環(huán)境變化國際合作聯(lián)合實(shí)驗(yàn)室/氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心/江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210044;

    ③ 中國氣象局 國家氣象中心,北京 100081;

    ④ 河南大學(xué) 環(huán)境與規(guī)劃學(xué)院,河南 開封 475004;

    ⑤ 河南省人工影響天氣領(lǐng)導(dǎo)小組辦公室,河南 鄭州 450003

    2015-06-11收稿,2016-03-08接受

    公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201106020);干旱氣象科學(xué)研究資助項(xiàng)目(GYHY201506001)

    摘要利用黃淮海區(qū)域90個(gè)站點(diǎn)1971—2000年逐日氣象資料以及國家氣候中心發(fā)布的未來氣候變化情景(A1B)下區(qū)域氣候模式(RegCM3)模擬的黃淮海區(qū)域1951—2070年0.25°×0.25°格點(diǎn)氣象資料,結(jié)合夏玉米主要生育期對溫度的需求,構(gòu)建了黃淮海區(qū)域的溫度適宜度和變異系數(shù)模型,并對1951—2070年黃淮海區(qū)域熱量資源、夏玉米主要生育期的溫度適宜度及其變異系數(shù)的時(shí)空變化特征進(jìn)行分析。結(jié)果表明:1)黃淮海區(qū)域≥10 ℃積溫和80%保證率下日平均溫度≥10 ℃的初日均呈現(xiàn)由北向南依次增加的趨勢,且隨時(shí)間推移,分別呈增加和提前趨勢。2)黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度隨時(shí)間整體呈逐漸上升的變化趨勢、其變異系數(shù)隨時(shí)間呈降—升—降的變化趨勢;出苗—抽雄期的溫度適宜度隨時(shí)間呈先降后升的變化趨勢、其變異系數(shù)呈降—升—降—升的變化趨勢;抽雄—成熟期的溫度適宜度空間上呈現(xiàn)2010年前北低南高、未來情景下中部低四周高的分布趨勢,時(shí)間上呈2010年前穩(wěn)定、未來情景下先降后升的變化趨勢,其變異系數(shù)呈相反變化趨勢;3)黃淮海區(qū)域夏玉米溫度適宜度及其變異系數(shù)從播種—出苗期—出苗—抽雄期—抽雄—成熟期均呈反相位的變化關(guān)系。

    關(guān)鍵詞

    夏玉米

    熱量資源

    溫度適宜度

    變異系數(shù)

    IPCC第五次評估報(bào)告指出,19世紀(jì)以來,全球平均地表溫度(GMST)均呈現(xiàn)上升趨勢,1901—2012年,全球陸地溫度上升了0.89 ℃,1951—2012年,全球海洋溫度上升了0.72 ℃。在IPCC排放情景特別報(bào)告(SRES)中指出,1990—2100年,全球平均地表溫度將上升1.4~5.8 ℃(姚潔等,2010;李忠賢等,2011;杜美芳等,2015)。全球變暖將使未來高溫極端事件的頻率、持續(xù)時(shí)間、量級大大增加(秦大河等,2014),這將給人類生存、社會經(jīng)濟(jì)、生態(tài)系統(tǒng)等方面造成極大的影響(Rosenzweig and Parry,1994;Darwin et al.,1995;Tao et al.,2016)。

    玉米是C4作物,其光合生理特征使其在日益增長的大氣CO2濃度中受益,卻可能會遭致增溫帶來的生育期縮短、產(chǎn)量降低等不利影響(金之慶和葛道闊,1996)。黃淮海區(qū)域夏玉米種植面積居全國首位,但其生長季的氣候條件處于次適宜狀態(tài),導(dǎo)致其產(chǎn)量不穩(wěn)(崔讀昌等,1984)。為定量分析吉林省農(nóng)作物生長季熱量資源對主要農(nóng)作物各生長發(fā)育期的滿足程度,馬樹慶(1996)建立了適宜度函數(shù),之后該函數(shù)普遍用于夏玉米、冬小麥和水稻適宜度研究。代立芹等(2010)對河北省冬小麥降水適宜度及時(shí)空特征進(jìn)行分析,得出降水不足是限制河北省冬小麥生產(chǎn)最重要的因子。李秀芬等(2013)應(yīng)用WOFOST(WOrld FOod STudies)模型模擬玉米生物量,并建立了東北地區(qū)玉米生長發(fā)育期氣象條件適宜度指數(shù)。目前關(guān)于適宜度的研究大多集中在某種作物,對于作物種植風(fēng)險(xiǎn)的研究相對較少,且時(shí)間序列較短。

    本文在前人研究的基礎(chǔ)上,基于區(qū)域氣候模式RegCM3輸出的未來氣候情景(SRES A1B情景,Special Report on Emissions Scenarios the A1 storyline and scenario family featured by a Balance across all sources)下的氣候資料,結(jié)合黃淮海區(qū)域夏玉米生長發(fā)育參數(shù)和夏玉米對溫度的需求指標(biāo),建立了黃淮海區(qū)域夏玉米溫度適宜度模型,針對未來情景氣候變化過程中黃淮海區(qū)域熱量資源的變化情況、夏玉米溫度適宜度及其變異系數(shù)的時(shí)空變化特征進(jìn)行分析,為該區(qū)域應(yīng)對氣候變化、促進(jìn)高產(chǎn)穩(wěn)產(chǎn)、合理利用農(nóng)業(yè)資源和調(diào)整農(nóng)業(yè)結(jié)構(gòu)提供了理論依據(jù)。

    1資料與方法

    1.1研究區(qū)概況

    黃淮海區(qū)域包括黃河、淮河、海河中下游流域的北京、天津、山東、河南全部、河北大部以及江蘇安徽兩省北部。該區(qū)域?qū)贉貛О霛駶櫄夂?無霜期170~220 d,多年平均年降水量500~1 000 mm左右,平原區(qū)一般為500~800 mm,年降水量的60%集中在6—8月;山丘區(qū)一般為1 000 mm上下,年降水約70%~80%集中于6—9月;夏玉米全生育期需水量300~400 mm,需水高峰期一般為8月;玉米播種面積占全國31.98%,產(chǎn)量占全國32.92%,占糧食總產(chǎn)量25.95%(李運(yùn)朝等,2004)。

    1.2資料來源

    本研究數(shù)據(jù)來源于國家氣候中心RegCM3模式模擬A1B情景下1951—2070年0.25°×0.25°格點(diǎn)氣象資料日值,包括日平均氣溫、日最高和最低氣溫、日降水量、日平均風(fēng)速、日平均相對濕度、日總輻射、日凈輻射;以及國家信息中心發(fā)布的1971—2000年黃淮海區(qū)域90個(gè)站點(diǎn)(圖1)逐日氣象資料,包括日平均氣溫、日最高和最低氣溫、日降水量、日平均風(fēng)速、日平均相對濕度。

    圖1 黃淮海區(qū)域氣象站點(diǎn)分布Fig.1 Distribution of meteorological stations in the Huang-Huai-Hai region

    1.3數(shù)據(jù)訂正

    利用Matlab讀取RegCM3模式模擬黃淮海區(qū)域1951—2070年的逐日平均溫度、日最高和最低溫度、日降水量、日平均風(fēng)速、日相對濕度;利用雙線性插值法將RegCM3模式數(shù)據(jù)進(jìn)行插值得到1971—2000年逐日站點(diǎn)的模式數(shù)據(jù);用1971—2000年歷史數(shù)據(jù)與模式數(shù)據(jù)的插值,對黃淮海區(qū)域1951—2070年0.25°×0.25°的模式數(shù)據(jù)進(jìn)行訂正(冶明珠等,2012)。

    1.4研究方法

    將1951—2070年的數(shù)據(jù)資料分成4個(gè)時(shí)段進(jìn)行分析,分別為1951—1980、1981—2010、2011—2040、2041—2070。

    1.4.1夏玉米生育期的確定

    采用5 d滑動平均法確定逐年穩(wěn)定通過10 ℃的初日作為夏玉米的適播期(孫衛(wèi)國,2008;孫孟梅等,1998);夏玉米生長的終止期采用積溫劃定。根據(jù)黃淮海區(qū)域主要農(nóng)業(yè)氣象站點(diǎn)夏玉米各生育期的起止日期,計(jì)算各站點(diǎn)夏玉米各生育期所需積溫,之后將整個(gè)黃淮海區(qū)域取平均值,為了計(jì)算方便,在不影響整體結(jié)果的情況下,對其進(jìn)行取整,具體各生育期積溫指標(biāo)參見表1。

    表1黃淮海區(qū)域夏玉米各生育期積溫指標(biāo)

    Table 1Accumulated temperature index of summer maize in each growth period in the Huang-Huai-Hai region

    ℃·d

    適播期與夏玉米生長終止期的間隔天數(shù)為夏玉米生長季長度。生長季內(nèi)≥10 ℃的活動積溫反映了夏玉米全生育期對熱量的要求。對≥10 ℃初日、生長季長度及熱量資源按80%保證率(陳群等,2014)取值。指定概率條件下界限值的計(jì)算公式(顏宏等,2002)如下:

    x=(1-a)xj+axj+1。

    (1)

    其中:j為序號,j=p(n+1),j取整數(shù),p為指定概率(又稱保證率,下同),n為記錄個(gè)數(shù);xj做由大到小的排序;a為中間變量,a=p(n+1)-j。

    1.4.2溫度適宜度函數(shù)的確定

    溫度適宜度計(jì)算公式(馬樹慶,1996)如下:

    (2)

    其中:T為不同生育期內(nèi)日平均氣溫,T0、T1、T2分別為夏玉米不同生育期內(nèi)生長發(fā)育所需的最適、下限和上限溫度(表2)(信乃詮等,1999)。

    表2黃淮海區(qū)域夏玉米各生育期三基點(diǎn)溫度

    Table 2Three fundamental temperature points in each growth period of summer maize in the Huang-Huai-Hai region

    1.4.3變異系數(shù)

    變異系數(shù)CV(Coefficient of Variation)是根據(jù)變量平均值的離散程度大小,來客觀地對其賦權(quán)(羅良清和魏和清,2011)。本文以黃淮海區(qū)域夏玉米不同發(fā)育期的溫度適宜度的變異系數(shù)來表示其穩(wěn)定情況。計(jì)算公式如下:

    (3)

    2結(jié)果與分析

    發(fā)育期內(nèi)≥10 ℃積溫的變化反應(yīng)了夏玉米對熱量資源的需求,在農(nóng)業(yè)氣象上,本文以穩(wěn)定通過10 ℃的初日為夏玉米適播期,根據(jù)積溫條件依次確定1951—2070年夏玉米的播種、出苗、抽雄和成熟日期。由式(2)、(3)分別計(jì)算得到1951—2070年4個(gè)不同時(shí)段夏玉米不同生育階段的溫度適宜度及其變異系數(shù)。

    2.1黃淮海區(qū)域熱量資源變化特征

    2.1.1黃淮海區(qū)域夏玉米生育期內(nèi)≥10 ℃積溫

    由圖2可以看出,黃淮海區(qū)域≥10 ℃積溫呈現(xiàn)由北向南依次增加的趨勢,且隨著時(shí)間的推移,整個(gè)黃淮海區(qū)域≥10 ℃積溫呈增加趨勢。

    1951—1980年(簡稱P1),除河北中部、北京北部和河南西部小部分區(qū)域,其他區(qū)域≥10 ℃積溫均大于3 600 ℃·d;相對于P1,1981—2010年(簡稱P2)安徽北部和河南東南部≥10 ℃積溫在4 600 ℃·d以上的區(qū)域面積有所擴(kuò)大;相對于P2,2011—2040年(簡稱P3)黃淮海區(qū)域除河北西部和北部、北京北部、河南西部、山東中部和東部,大部分區(qū)域≥10 ℃積溫增加到4 600 ℃·d以上;相對于P3,2041—2070年(簡稱P4)黃淮海區(qū)域≥10 ℃積溫大部分已達(dá)4 600 ℃·d以上,其中河南西南部和北部的小部分區(qū)域、山東西部、安徽北部增加到5 600 ℃·d以上。

    2.1.2黃淮海區(qū)域80%保證率下日平均氣溫≥10 ℃初日

    圖2 黃淮海區(qū)域夏玉米生育期內(nèi)≥10 ℃積溫  a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.2 The accumulated temperature of 10 ℃ in the maize growth season in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    圖3 黃淮海區(qū)域80%保證率下日平均氣溫≥10 ℃初日(儒略日)a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.3 The first day daily average temperature≥10 ℃ with a guaranteed rate of 80% in the Huang-Huai-Hai region(Julian days):(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    圖4 黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.4 Temperature suitability of summer maize in the stage of seeding to emergence in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    如圖3所示,黃淮海區(qū)域80%保證率下日平均溫度≥10 ℃的初日呈現(xiàn)由北向南逐漸提前的趨勢,隨時(shí)間推移,黃淮海區(qū)域80%保證率下日平均溫度≥10 ℃的初日整體呈現(xiàn)提前趨勢。P1黃淮海區(qū)域≥10 ℃初日大部分在3月中下旬及以后,河北北部≥10 ℃初日出現(xiàn)最晚,達(dá)4月中上旬。相對于P1,P2河北東南部、山東中部≥10 ℃初日提前了18 d左右,達(dá)3月中下旬。相對于P2,P3黃淮海區(qū)域≥10 ℃初日變化不大。相對于P3,P4黃淮海區(qū)域≥10 ℃初日有所提前,在3月中下旬范圍的區(qū)域面積有所擴(kuò)大,包括安徽北部、江蘇北部、河南中部和北部、山東全部、河北南部、北京南部和天津全部,河南南部和安徽西北部≥10 ℃初日達(dá)2月下旬或3月上旬。

    2.2黃淮海區(qū)域夏玉米不同生育期的溫度適宜度及變異系數(shù)

    2.2.1播種—出苗期

    從圖4可以看出,黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度隨時(shí)間整體呈逐漸上升的趨勢。P1黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度大部分在0.31~0.40之間。其中河北、天津、山東東部較低,在0.3以下。安徽北部局部地區(qū)溫度適宜度較高,達(dá)0.41~0.50;相對于P1,P2黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度變化不大,其中江蘇東部和安徽北部局部地區(qū)下降了0.1左右,山東西部局部地區(qū)上升了0.1左右;相對于P2,P3黃淮海區(qū)域的溫度適宜度有所上升,大部分在0.3以上,其中河北中西部、山東中北部達(dá)0.4以上,說明在該時(shí)段內(nèi)黃淮海區(qū)域播種—出苗期的溫度條件更適宜于夏玉米的生長;相對于P3,P4年黃淮海區(qū)域的溫度適宜度繼續(xù)上升,大部分區(qū)域達(dá)到了0.4以上,根據(jù)上文得出的該時(shí)間段內(nèi)積溫增加的結(jié)論,由于夏玉米從播種到出苗所需要的積溫是固定值,若該時(shí)間段內(nèi)的積溫增加,將更利于夏玉米該階段的生長,夏玉米該階段所需的生長時(shí)間將會縮短,溫度適宜度將會升高。

    圖5為黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度變異系數(shù),可以看出:P1,黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度變異系數(shù)大部分在0.26~0.50之間,河南中東部和山東東南部局部地區(qū)變異系數(shù)較大,說明該時(shí)段內(nèi),上述區(qū)域的溫度適宜度有所波動,應(yīng)注意密切關(guān)注溫度變化,以及時(shí)采取應(yīng)對措施;相對于P1,P2黃淮海大部分區(qū)域的溫度適宜度變異系數(shù)有所下降,河北南部、山東中西部、河南大部、江蘇和安徽北部下降了0.1左右,其中山東東南部小部分區(qū)域下降了0.35左右,說明該時(shí)間段內(nèi)黃淮海區(qū)域的溫度適宜度相對穩(wěn)定;相對于P2,2011—2070年黃淮海區(qū)域的溫度適宜度變異系數(shù)持續(xù)上升,幅度在0.15左右。

    2.2.2出苗—抽雄期

    從圖6可以看出,黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度整體高于播種—出苗期,且時(shí)間上呈先降后升的分布趨勢。P1,黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度大部分在0.61~0.70之間,其中河北中西部、河南西部、山東中東部、江蘇東北部溫度適宜度較低,在0.6以下;相對于P1,P2黃淮海區(qū)域除河北中部、河南西南部、安徽北部局部地區(qū),其它區(qū)域的溫度適宜度均下降0.1左右;相對于P2,P3整個(gè)黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度有明顯上升,幅度達(dá)0.1~0.2,大部分區(qū)

    圖5 黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度變異系數(shù)a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.5 CV of temperature suitability of summer maize in the stage of seeding to emergence in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    圖6 黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.6 Temperature suitability of summer maize in the stage of emergence to heading in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    域在0.6以上,其中河北南部、河南、山東北部達(dá)0.7以上;相對于P3,P4黃淮海區(qū)域的溫度適宜度變化不大,主要表現(xiàn)為河北中部、山東西部有0.1左右的上升。

    圖7為黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度變異系數(shù),由圖可知:P1,黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度變異系數(shù)大部分在0.2以下,其中河北中部、河南西南部、山東東南部、安徽北部和江蘇東北局部較高,達(dá)0.21~0.30,河北中部小部分區(qū)域較小,在0.15以下;相對于P1,P2黃淮海大部分區(qū)域的溫度適宜度變異系數(shù)有所下降,河南西南部和安徽北部下降幅度較大,達(dá)0.05~0.10,說明該時(shí)段內(nèi),上述區(qū)域溫度適宜度波動較小;相對于P2,P3黃淮海區(qū)域內(nèi)河南南部、安徽北部和江蘇西北部溫度適宜度變異系數(shù)上升0.05~0.10,未來情景下,上述區(qū)域應(yīng)根據(jù)具體溫度變化,合理調(diào)整種植制度,河北中部、山東東南部、江蘇東北部下降0.05左右,其它區(qū)域變化不大;相對于P3,P4黃淮海區(qū)域的溫度適宜度變異系數(shù)有0.05左右的下降,除河南西部、安徽西北局部,大部分區(qū)域在0.2以下,說明該時(shí)間段的夏玉米溫度適宜度相對穩(wěn)定。

    2.2.3抽雄—成熟期

    從圖8可以看出,黃淮海區(qū)域夏玉米抽雄—成熟期的溫度適宜度整體高于出苗—抽雄期,時(shí)間上呈現(xiàn)2010年前較穩(wěn)定、未來情景下逐漸下降趨勢,空間上呈現(xiàn)2010年之前北低南高、未來情景下中部低四周高的分布趨勢。1951—2010年,黃淮海區(qū)域除河北北部局部區(qū)域,其他大部分區(qū)域夏玉米抽雄—成熟期的溫度適宜度在0.85以上,說明該階段溫度條件能夠基本滿足夏玉米抽雄—成熟期的溫度要求;相對于P2,在未來情景下,P3河北、北京南部、天津大部、河南東北部、山東西北部溫度適宜度

    圖7 黃淮海區(qū)域夏玉米出苗—抽雄期的溫度適宜度變異系數(shù)a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.7 CV of temperature suitability of summer maize in the stage of emergence to heading in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    圖8 黃淮海區(qū)域夏玉米抽雄—成熟期的溫度適宜度a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.8 Temperature suitability of summer maize in the stage of heading to maturity in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    有所下降,在0.15~0.20,其他區(qū)域基本在0.85以上;P4,黃淮海區(qū)域整體呈下降趨勢,其中河北南部、北京、天津、河南大部、山東中西部、安徽、江蘇北部下降幅度較明顯,達(dá)0.15~0.30,尤以河北南部、山東西北部最為明顯。

    圖9為黃淮海區(qū)域夏玉米抽雄—成熟期的溫度適宜度變異系數(shù),結(jié)合圖8可以看出,該時(shí)間段夏玉米抽雄—成熟期的溫度適宜度變異系數(shù)時(shí)空變化時(shí)間上呈現(xiàn)2010年前較穩(wěn)定、未來情景下逐漸上升的趨勢。1951—2010年,黃淮海區(qū)域夏玉米抽雄—成熟期的溫度適宜度變異系數(shù)基本在0.2以下,說明該時(shí)間段溫度適宜度波動極小。在未來情景下,2011—2070年,黃淮海區(qū)域除河北南部、天津西南部、北京大部、河南東北部和山東西北部溫度適宜度變異系數(shù)持續(xù)上升,應(yīng)密切注意未來情景下的溫度條件,以應(yīng)對夏玉米抽雄—成熟期的溫度需求。其它區(qū)域的溫度適宜度變異系數(shù)基本維持在0.2以下。

    圖9 黃淮海區(qū)域夏玉米抽雄—成熟期的溫度適宜度變異系數(shù)a.1951—1980年;b.1981—2010年;c.2011—2040年;d.2041—2070年Fig.9 CV of temperature suitability of summer maize in the stage of heading to maturity in the Huang-Huai-Hai region:(a)1951—1980;(b)1981—2010;(c)2011—2040;(d)2041—2070

    3結(jié)論和討論

    黃淮海區(qū)域干濕季節(jié)分明,熱量資源對于夏玉米的生長發(fā)育尤為重要。本文在氣候變化的背景下針對性地研究了黃淮海區(qū)域熱量資源以及夏玉米溫度適宜度及其變異,為提高氣候資源利用效率、促進(jìn)夏玉米合理種植提供了科學(xué)依據(jù)。

    1)黃淮海區(qū)域≥10 ℃積溫呈現(xiàn)由北向南依次增加的趨勢,且隨時(shí)間推移,整個(gè)黃淮海區(qū)域≥10 ℃積溫呈增加趨勢。

    2)黃淮海區(qū)域80%保證率下日平均溫度≥10 ℃的初日呈現(xiàn)由北向南逐漸提前的趨勢,且隨時(shí)間推移,整體呈提前趨勢。

    3)黃淮海區(qū)域夏玉米播種—出苗期的溫度適宜度隨時(shí)間整體逐漸上升的變化趨勢,該生育期的溫度適宜度變異系數(shù)隨時(shí)間呈降—升—降的變化趨勢;出苗—抽雄期的溫度適宜度時(shí)間上呈先降后升的變化趨勢,該生育期的溫度適宜度變異系數(shù)隨時(shí)間呈降—升—降的變化趨勢;抽雄—成熟期空間上呈現(xiàn)2010年前北低南高、未來情景下中部低四周高的分布趨勢,時(shí)間上呈2010年前較穩(wěn)定、未來情景下逐漸下降的變化趨勢,該生育期的溫度適宜度變異系數(shù)隨時(shí)間呈2010年前較穩(wěn)定、未來情景下逐漸上升的變化趨勢。

    4)黃淮海區(qū)域夏玉米溫度適宜度從播種—出苗期—出苗—抽雄期—抽雄—成熟期呈依次上升趨勢,其變異系數(shù)呈依次下降趨勢,二者呈反相位的變化關(guān)系。

    目前關(guān)于氣候變化對作物影響方面的研究很多,主要側(cè)重于對作物播種期和生育期的變化、產(chǎn)量以及作物種植結(jié)構(gòu)變化的影響方面,而且大多側(cè)重于對歷史事實(shí)的檢測,而對未來黃淮海區(qū)域夏玉米溫度適宜度狀況及其變異情況尚缺乏系統(tǒng)研究。因此,研究未來情景下黃淮海區(qū)域夏玉米溫度適宜度時(shí)空演變趨勢及其波動情況,對充分利用氣候資源,指導(dǎo)黃淮海區(qū)域夏玉米生產(chǎn)具有重要意義和價(jià)值。在氣候變暖背景下,夏玉米的品種和耕作制度將會有一定變化,夏玉米的生長發(fā)育過程也將發(fā)生變化。目前關(guān)于夏玉米適宜度的研究主要在區(qū)域氣候模式數(shù)據(jù)的基礎(chǔ)上,并結(jié)合現(xiàn)有品種的生理指標(biāo)進(jìn)行的,因此通過結(jié)合田間試驗(yàn)和未來品種夏玉米的預(yù)測指標(biāo)進(jìn)行夏玉米適宜度的研究是未來研究的方向。

    參考文獻(xiàn)(References)

    陳群,于歡,侯雯嘉,等.2014.氣候變暖對黃淮海區(qū)域冬小麥生育進(jìn)程與產(chǎn)量的影響[J].麥類作物學(xué)報(bào),34(10):1363-1372.Chen Q,Yu H,Hou W J,et al.2014.Impacts of climate warming on growth development process and yield of winter wheat in Huang-huai-hai region of China[J].J Triticeae Crops,34(10):1363-1372.(in Chinese).

    崔讀昌,劉洪順,閔瑾如,等.1984.中國主要農(nóng)作物氣候資源圖集[M].北京:氣象出版社.Cui D C,Liu H S,Min J R,et al.1984.Atlas of climate resources for main crops in China[M].Beijing:China Meteorological Press.(in Chinese).

    代立芹,李春強(qiáng),姚樹然,等.2010.氣候變暖背景下河北省冬小麥凍害變化分析[J].中國農(nóng)業(yè)氣象,31(3):467-471.Dai L Q,Li C Q,Yao S R,et al.2010.Variation analysis of freezing injury on winter wheat under climate warming in Hebei province[J].Chinese Journal of Agrometeorology,31(3):467-471.(in Chinese).

    Darwin R,Tsigas M,Lewandrowski J,et al.1995.World agriculture and climate change:Economic adaptations[R].United States Department of Agriculture,Economic Research Service,Agricultural Economics Reports.

    杜美芳,徐海明,鄧潔淳.2015.基于CMIP5資料的熱帶太平洋和印度洋海溫變化與降水變化的關(guān)系及其成因分析[J].大氣科學(xué)學(xué)報(bào),38(6):742-752.Du M F,Xu H M,Deng J C.2015.Relationship between sea surface temperature variation and precipitation variation over tropical Pacific and Indian Ocean and its cause analysis based on CMIP5 data[J].Trans Atmos Sci,38(6):742-752.(in Chinese).

    金之慶,葛道闊.1996.評價(jià)全球氣候變化對我國玉米生產(chǎn)的可能影響[J].作物學(xué)報(bào),22(5):513-524.Jin Z Q,Ge D K.1996.Assessing the potential impacts of global climate change on maize production in China[J].Acta Agron Sin,22(5):513-524.(in Chinese).

    李秀芬,馬樹慶,宮麗娟,等.2013.基于WOFOST的東北地區(qū)玉米生育期氣象條件適宜度評價(jià)[J].中國農(nóng)業(yè)氣象,34(1):43-49.Li X F,Ma S Q,Gong L J,et al.2013.Evaluation of meteorological suitability degree during maize growth period based on wofost in northeast China[J].Chinese Journal of Agrometeorology,34(1):43-49.(in Chinese).

    李運(yùn)朝,王元東,崔彥宏,等.2004.玉米抗旱性鑒定研究進(jìn)展[J].玉米科學(xué),12(1):63-68.Li Y C,Wang Y D,Cui Y H,et al.2004.Approach of study on maize drought-resistance identification[J].J Maize Sci,12(1):63-68.(in Chinese).

    李忠賢,陳海山,曾剛,等.2011.海溫強(qiáng)迫下的東亞夏季大氣環(huán)流潛在可預(yù)報(bào)性特征[J].大氣科學(xué)學(xué)報(bào),34(3):281-287.Li Z X,Chen H S,Zeng G,et al.2011.Potential predictability of East Asian summer atmospheric circulation forced by sea surface temperature[J].Trans Atmos Sci,34(3):281-287.(in Chinese)

    羅良清,魏和清.2011.統(tǒng)計(jì)學(xué)[M].北京:中國財(cái)政經(jīng)濟(jì)出版社.Luo L Q,Wei H Q.2011.Statistics[M].Beijing:China Financial Economic Publishing House.(in Chinese).

    馬樹慶.1996.吉林省農(nóng)業(yè)氣候研究[M].北京:氣象出版社.Ma S Q.1996.Agricultural climate research in Jilin province[M].Beijing:China Meteorological Press.(in Chinese).

    秦大河,Thomas S,沈永平,等.2014.IPCC第五次評估報(bào)告第一工作組報(bào)告的亮點(diǎn)結(jié)論[J].氣候變化研究進(jìn)展,10(1):1-6.Qin D H,Thomas S,Shen Y P,et al.2014.Highlights of the IPCC working groupⅠfifth assessment report[J].Adv Climate Change Res,10(1):1-6.(in Chinese).

    Rosenzweig C,Iglesias A,Yang X B,et al.2001.Climate change and extreme weather events:Implications for food production,plant diseases,and pests[J].Global Change & Human Health,2(2):90-104.

    Rosenzweig C,Parry M L.1994.Potential impact of climate change on world food supply[J].Nature,367(6459):133-138.

    孫孟梅,姜麗霞,于榮環(huán),等.1998.玉米生育期熱量指標(biāo)及其不同品種栽培北界的研究[J].中國農(nóng)業(yè)氣象,19(4):8-16.(in Chinese).Sun M M,Jiang L X,Yu R H,et al.1998.Study on heat index at growing stages of corn and the planting border of different varieties in the north[J].Chinese Journal of Agrometeorology,19(4):8-16.(in Chinese).

    孫衛(wèi)國.2008.氣候資源學(xué)[M].北京:氣象出版社.Sun W G.2008.Science of climatic resources[M].Beijing:China Meteorological Press.(in Chinese).

    Tao S,Shen S H,Li Y,et al.2016.Projected crop production under regional climate change using scenario data and modeling:Sensitivity to chosen sowing date and cultivar[J].Sustainability,8(3):214.

    信乃詮,崔讀昌,高亮之,等.1999.中國農(nóng)業(yè)氣象學(xué)[M].北京:中國農(nóng)業(yè)出版社.Xin N Q,Cui D C,Gao L Z,et al.1999.Agricultural meteorology of China[M].Beijing:China Agriculture Press.(in Chinese).

    顏宏,沈國權(quán),毛耀順.2002.中華人民共和國氣候圖集[M].北京:氣象出版社.Yan H,Shen G Q,Mao Y S.2002.The climatic atlas of the People’s Republic of China[M].Beijing:China Meteorological Press.(in Chinese).

    姚潔,陳海山,朱偉軍.2010.北半球陸面過程對全球變暖響應(yīng)特征的初步分析[J].大氣科學(xué)學(xué)報(bào),33(2):220-226.Yao J,Chen H S,Zhu W J.2010.Prelimilary analysis on response features of land surface processes over Northern Hemisphere to global warming[J].Tans Atmos Sci,33(2):220-226.(in Chinese).

    冶明珠,郭建平,袁彬,等.2012.氣候變化背景下東北地區(qū)熱量資源及玉米溫度適宜度[J].應(yīng)用生態(tài)學(xué)報(bào),23(10):2786-2794.Ye M Z,Guo J P,Yuan B,et al.2012.Thermal resources and maize temperature suitability in northeast China under climate change[J].Chinese Journal of Applied Ecology,23(10):2786-2794.(in Chinese).

    Based on the daily meteorological data of 90 stations for the period 1971 to 2000 and daily 0.25°×0.25° gridded meteorological data under the A1B climate scenario(1951—2070) extracted from the regional climate model RegCM3 released by the National Climate Center,a temperature suitability and coefficient of variation(CV) model for the Huang-Huai-Hai region was built,combing with the growth demand of summer maize for temperature.The spatial and temporal variations of thermal resources,summer maize temperature suitability and its CV in the Huang-Huai-Hai region for the period 1951—2070 were analyzed.The Huang-Huai-Hai area has distinct dry and wet seasons.The heat resources are particularly important for the growth of summer maize.The analysis of this region in this paper,against the background of climate change,is aimed at improving the efficiency of the utilization of climate resources,and providing a scientific basis for summer maize planting.The results showed that:

    (1)The accumulated temperature of 10 ℃ with a guaranteed rate of 80% in the Huang-Huai-Hai region showed a gradually increasing trend from north to south,increasing over time.(2)The first day of 10 ℃ with a guaranteed rate of 80% in the Huang-Huai-Hai region showed a gradually increasing trend from north to south,advancing over time.(3)In the stage of seeding to emergence,the temperature suitability showed a trend of rising firstly and then falling,with the CV falling firstly and then rising over time.In the stage of emergence to heading,a falling-rising-falling trend was apparent,with the CV showing a falling-rising-falling-rising trend over time.In the stage of heading to maturity,spatially,the temperature suitability was higher in the south than the north before 2010,but showed a future pattern of being higher in the center than the surroundings.Temporally,it was stable before 2010,but showed a falling followed by rising pattern for the future,which was opposite to the pattern for the CV.(4)The temperature suitability and CV fluctuated in positive phase with time from seeding to maturity.

    To date,most research on summer maize suitability has mainly been based on regional climate model data,combined with physiological indexes of existing varieties,which mainly focus on the change in seeding time,growth period,yield,and the effect of crop planting structure by analyzing existing dates.Crucially,there remains a lack of systematic research on the variation of summer maize temperature suitability in the future.Therefore,combining field results with future projections,such as the present work in the Huang-Huai-Hai region,may be the most appropriate research direction moving forwards.This study of the spatiotemporal trends of summer maize temperature suitability and its variation under a future climate scenario in the Huang-Huai-Hai region should help to make full use of climate resources,which may also guide the production of summer maize in the region and provide important information for decision-makers.

    summer maize;thermal resources;temperature suitability;coefficient of variation

    (責(zé)任編輯:孫寧)

    Thermal resources and summer maize temperature suitability in the Huang-Huai-Hai region under future climate change

    LI Meng1,2,SHEN Shuanghe1,2,LYU Houquan3,HAN Yan4,CHU Ronghao1,SHA Xiuzhu5

    1CollegeofAppliedMeteorology,NanjingUniversityofInformationScience&Technology,Nanjing210044,China;2KeyLaboratoryofMeteorologicalDisasters,MinistryofEducation(KLME)/JointInternationalResearchLaboratoryofClimateandEnvironmentChange(ILCEC)/CollaborativeInnovationCenteronForecastandEvaluationofMeteorologicalDisasters(CIC-FEMD)/JiangsuKeyLaboratoryofAgriculturalMeteorology,NanjingUniversityofInformationScience&Technology,Nanjing210044,China;3NationalMeteorologicalCenterofChinaMeteorologicalAdministration,Beijing100081,China;4CollegeofEnvironmentandPlanning,HenanUniversity,Kaifeng475004,China;5WeatherModificationLeadingGroupOfficeofHenanProvince,Zhengzhou450003,China

    doi:10.13878/j.cnki.dqkxxb.20150611001

    引用格式:李萌,申雙和,呂厚荃,等.2016.氣候變化情景下黃淮海區(qū)域熱量資源及夏玉米溫度適宜度[J].大氣科學(xué)學(xué)報(bào),39(3):391-399.

    LiM,ShenSH,LyuHQ,etal.2016.ThermalresourcesandsummermaizetemperaturesuitabilityintheHuang-Huai-Hairegionunderfutureclimatechange[J].TransAtmosSci,39(3):391-399.doi:10.13878/j.cnki.dqkxxb.20150611001.(inChinese).

    *聯(lián)系人,E-mail:yqzhr@nuist.edu.cn

    猜你喜歡
    黃淮海積溫夏玉米
    2022 年國審黃淮海夏播玉米新品種
    黃淮海北部地區(qū)夏玉米穩(wěn)產(chǎn)高產(chǎn)的播期優(yōu)選
    湟源縣30年0℃活動積溫變化特征分析
    小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
    石河子地區(qū)近50年積溫變化特征分析
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    黃淮海地區(qū)高蛋白夏大豆高產(chǎn)栽培技術(shù)
    大豆科技(2016年6期)2016-03-11 15:30:22
    石河子市近51a≥0℃積溫變化特征分析
    伊春地區(qū)近40 a積溫變化的動態(tài)分析
    黑龍江氣象(2012年3期)2012-12-25 09:38:22
    天津市2011年良種補(bǔ)貼推介品種目錄
    在线观看一区二区三区| av国产免费在线观看| 欧美日韩综合久久久久久| 蜜桃在线观看..| 涩涩av久久男人的天堂| 亚洲三级黄色毛片| 国产精品99久久久久久久久| 国产午夜精品久久久久久一区二区三区| 男女下面进入的视频免费午夜| 十分钟在线观看高清视频www | 国产男人的电影天堂91| 免费观看性生交大片5| 亚洲婷婷狠狠爱综合网| 毛片女人毛片| 国产精品蜜桃在线观看| 日韩电影二区| 色综合色国产| 亚洲精品久久久久久婷婷小说| 看非洲黑人一级黄片| 精品人妻偷拍中文字幕| 国产精品人妻久久久影院| 97超视频在线观看视频| 视频中文字幕在线观看| 天堂8中文在线网| 毛片一级片免费看久久久久| 麻豆精品久久久久久蜜桃| 99热这里只有精品一区| 80岁老熟妇乱子伦牲交| 伦理电影大哥的女人| 大香蕉久久网| 精品久久久久久电影网| 夫妻午夜视频| 午夜福利视频精品| 尤物成人国产欧美一区二区三区| 色5月婷婷丁香| 亚洲欧美精品自产自拍| 免费高清在线观看视频在线观看| 亚洲国产精品专区欧美| 亚洲欧美清纯卡通| 特大巨黑吊av在线直播| av视频免费观看在线观看| 国产精品一区二区性色av| 日本av免费视频播放| 亚洲av男天堂| 欧美少妇被猛烈插入视频| 国产成人a∨麻豆精品| 国产片特级美女逼逼视频| 久久久久久久亚洲中文字幕| 日韩,欧美,国产一区二区三区| av.在线天堂| 高清黄色对白视频在线免费看 | 搡女人真爽免费视频火全软件| 日韩中文字幕视频在线看片 | 最后的刺客免费高清国语| 久久久久久久大尺度免费视频| 国产日韩欧美在线精品| 国产精品99久久久久久久久| 国产高潮美女av| av播播在线观看一区| www.色视频.com| 丰满迷人的少妇在线观看| 综合色丁香网| 97精品久久久久久久久久精品| 国产精品蜜桃在线观看| 男人添女人高潮全过程视频| 男人舔奶头视频| 高清日韩中文字幕在线| 激情 狠狠 欧美| 国产黄片视频在线免费观看| 欧美 日韩 精品 国产| h视频一区二区三区| 精品酒店卫生间| 久久精品国产自在天天线| 国产男女内射视频| 亚洲国产最新在线播放| 亚洲精品国产av成人精品| 国产无遮挡羞羞视频在线观看| 亚洲图色成人| 99热6这里只有精品| 国产欧美日韩精品一区二区| 国产av精品麻豆| 亚洲av成人精品一二三区| 国产精品精品国产色婷婷| 建设人人有责人人尽责人人享有的 | 交换朋友夫妻互换小说| 国产精品秋霞免费鲁丝片| 久久久色成人| 成人亚洲精品一区在线观看 | 亚洲丝袜综合中文字幕| 久久亚洲国产成人精品v| 日本av手机在线免费观看| 天堂中文最新版在线下载| 国产人妻一区二区三区在| 草草在线视频免费看| h视频一区二区三区| 噜噜噜噜噜久久久久久91| 人妻一区二区av| 最近2019中文字幕mv第一页| 在线观看国产h片| 欧美区成人在线视频| 亚洲欧美清纯卡通| 久久久久国产网址| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 久久 成人 亚洲| 亚洲精品日本国产第一区| 午夜日本视频在线| 热99国产精品久久久久久7| 国产精品成人在线| 国产在线男女| 亚洲欧美一区二区三区国产| 日韩三级伦理在线观看| 十八禁网站网址无遮挡 | 女人久久www免费人成看片| 亚洲欧美日韩另类电影网站 | 只有这里有精品99| 中文字幕久久专区| 激情五月婷婷亚洲| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 久久久久久久亚洲中文字幕| 啦啦啦视频在线资源免费观看| 国产av码专区亚洲av| 欧美丝袜亚洲另类| 亚洲欧美精品专区久久| 国产高清三级在线| 日韩免费高清中文字幕av| 日本vs欧美在线观看视频 | 国产 一区精品| 久久精品国产亚洲av涩爱| 人妻系列 视频| 男女边吃奶边做爰视频| 在线播放无遮挡| 国产老妇伦熟女老妇高清| 国产一区二区在线观看日韩| 免费看光身美女| 久久久久国产精品人妻一区二区| 国产欧美另类精品又又久久亚洲欧美| 国产熟女欧美一区二区| 校园人妻丝袜中文字幕| 国内少妇人妻偷人精品xxx网站| 国产爽快片一区二区三区| 大片电影免费在线观看免费| 亚洲国产高清在线一区二区三| 亚洲国产精品国产精品| 国产高清有码在线观看视频| 中文天堂在线官网| 涩涩av久久男人的天堂| 国产av国产精品国产| 日本wwww免费看| 看免费成人av毛片| 国产精品嫩草影院av在线观看| 91久久精品国产一区二区三区| 午夜福利在线在线| 久久亚洲国产成人精品v| 国产在线一区二区三区精| 久久午夜福利片| 日韩人妻高清精品专区| 国产伦精品一区二区三区四那| 在线免费十八禁| 亚洲,一卡二卡三卡| 日本vs欧美在线观看视频 | 99久国产av精品国产电影| 久久精品国产自在天天线| 欧美人与善性xxx| 久久久久久久国产电影| 欧美日韩国产mv在线观看视频 | 亚洲内射少妇av| 亚洲欧美一区二区三区黑人 | 亚洲人成网站在线观看播放| 国产精品女同一区二区软件| 人人妻人人添人人爽欧美一区卜 | 国产综合精华液| 精品一区二区免费观看| 国产伦精品一区二区三区视频9| 国产精品无大码| 黄片无遮挡物在线观看| 欧美bdsm另类| 少妇被粗大猛烈的视频| 日韩大片免费观看网站| 国产一区二区在线观看日韩| 色网站视频免费| 黑丝袜美女国产一区| 男女免费视频国产| 亚洲第一区二区三区不卡| 亚洲精品亚洲一区二区| 在线观看av片永久免费下载| 日韩av在线免费看完整版不卡| 亚洲av欧美aⅴ国产| 男人舔奶头视频| 高清午夜精品一区二区三区| 少妇高潮的动态图| 看十八女毛片水多多多| 色综合色国产| 精品国产一区二区三区久久久樱花 | 夜夜骑夜夜射夜夜干| 久久97久久精品| 色吧在线观看| 青春草亚洲视频在线观看| 中国国产av一级| 校园人妻丝袜中文字幕| 一级毛片久久久久久久久女| 亚洲精品国产色婷婷电影| 午夜激情久久久久久久| 久久ye,这里只有精品| 久久综合国产亚洲精品| 成人无遮挡网站| 国产人妻一区二区三区在| 干丝袜人妻中文字幕| 精品亚洲乱码少妇综合久久| 国产91av在线免费观看| 亚洲三级黄色毛片| 国产免费又黄又爽又色| 日本欧美视频一区| 99热全是精品| 晚上一个人看的免费电影| av在线app专区| 欧美97在线视频| 亚洲国产精品999| 国产成人午夜福利电影在线观看| 人人妻人人爽人人添夜夜欢视频 | 噜噜噜噜噜久久久久久91| 久久99热这里只频精品6学生| 观看美女的网站| 欧美成人一区二区免费高清观看| 好男人视频免费观看在线| 色视频在线一区二区三区| 黑人猛操日本美女一级片| 色视频在线一区二区三区| 美国免费a级毛片| 狂野欧美激情性bbbbbb| 自线自在国产av| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品aⅴ一区二区三区四区| 成年动漫av网址| 一本一本久久a久久精品综合妖精| 性高湖久久久久久久久免费观看| 成年美女黄网站色视频大全免费| 免费观看a级毛片全部| 又大又黄又爽视频免费| 狂野欧美激情性xxxx| 婷婷丁香在线五月| 亚洲欧美日韩高清在线视频 | 999久久久国产精品视频| 日本五十路高清| 少妇被粗大的猛进出69影院| 国产一区二区三区综合在线观看| 国产伦人伦偷精品视频| 亚洲一区二区三区欧美精品| 青草久久国产| 久9热在线精品视频| 亚洲国产av新网站| 丰满人妻熟妇乱又伦精品不卡| 国产成人免费无遮挡视频| 夜夜骑夜夜射夜夜干| 啦啦啦啦在线视频资源| 悠悠久久av| 999久久久国产精品视频| 欧美av亚洲av综合av国产av| 在线观看免费视频网站a站| 宅男免费午夜| 久久鲁丝午夜福利片| 99国产精品一区二区蜜桃av | √禁漫天堂资源中文www| 精品一区二区三区av网在线观看 | 久久青草综合色| 人人妻人人添人人爽欧美一区卜| 人妻人人澡人人爽人人| 久久久久国产精品人妻一区二区| 水蜜桃什么品种好| 精品少妇一区二区三区视频日本电影| 国产精品久久久久久精品古装| 精品亚洲成国产av| 亚洲男人天堂网一区| 18禁国产床啪视频网站| 亚洲成av片中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 九草在线视频观看| 巨乳人妻的诱惑在线观看| 美女中出高潮动态图| 亚洲欧美精品自产自拍| 在线亚洲精品国产二区图片欧美| 亚洲国产精品成人久久小说| 亚洲国产精品成人久久小说| 欧美成人精品欧美一级黄| 午夜福利乱码中文字幕| 国产免费现黄频在线看| 狠狠婷婷综合久久久久久88av| 久久久久久久久久久久大奶| 欧美日韩亚洲综合一区二区三区_| 久久天躁狠狠躁夜夜2o2o | 脱女人内裤的视频| 一级黄色大片毛片| 一个人免费看片子| av有码第一页| 十八禁人妻一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产色婷婷电影| 日日夜夜操网爽| 久久人人爽人人片av| 国产在线观看jvid| 大香蕉久久网| 久久99精品国语久久久| 亚洲av成人精品一二三区| 国产精品三级大全| av天堂久久9| 一本久久精品| 欧美日韩视频高清一区二区三区二| 高清欧美精品videossex| 欧美成人精品欧美一级黄| 国产精品久久久av美女十八| 亚洲成av片中文字幕在线观看| 久久99精品国语久久久| 大片电影免费在线观看免费| 精品一区二区三区av网在线观看 | 国语对白做爰xxxⅹ性视频网站| a级毛片黄视频| 少妇精品久久久久久久| 国产精品久久久人人做人人爽| 久久人妻熟女aⅴ| 制服诱惑二区| 91麻豆av在线| 丰满少妇做爰视频| 少妇 在线观看| 国产精品久久久久久精品古装| 午夜免费观看性视频| av网站免费在线观看视频| e午夜精品久久久久久久| 91精品伊人久久大香线蕉| 久久精品久久久久久久性| 日韩视频在线欧美| 黑丝袜美女国产一区| 欧美成狂野欧美在线观看| 激情五月婷婷亚洲| 久久久久网色| 涩涩av久久男人的天堂| 视频区欧美日本亚洲| 午夜av观看不卡| xxx大片免费视频| netflix在线观看网站| 亚洲国产成人一精品久久久| 国产精品二区激情视频| 欧美精品人与动牲交sv欧美| 大片电影免费在线观看免费| 亚洲欧美一区二区三区久久| 老熟女久久久| 国产男女超爽视频在线观看| 永久免费av网站大全| 中文字幕色久视频| 国产1区2区3区精品| 欧美成人精品欧美一级黄| 亚洲精品美女久久久久99蜜臀 | 午夜免费鲁丝| 国产熟女欧美一区二区| 日本色播在线视频| 天堂8中文在线网| 另类精品久久| 日韩精品免费视频一区二区三区| 日本a在线网址| 中文精品一卡2卡3卡4更新| 一本综合久久免费| 亚洲精品在线美女| 亚洲成av片中文字幕在线观看| 久9热在线精品视频| 人人澡人人妻人| 午夜视频精品福利| 国产一区二区激情短视频 | 久久女婷五月综合色啪小说| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 国产亚洲av高清不卡| 久久亚洲精品不卡| av网站免费在线观看视频| 一级片'在线观看视频| 欧美日韩av久久| 久久亚洲精品不卡| 国产女主播在线喷水免费视频网站| av又黄又爽大尺度在线免费看| www日本在线高清视频| 国产精品麻豆人妻色哟哟久久| 人成视频在线观看免费观看| 中国国产av一级| 亚洲专区中文字幕在线| 国产在线免费精品| 少妇人妻久久综合中文| 精品福利观看| 国产xxxxx性猛交| 伊人亚洲综合成人网| 免费一级毛片在线播放高清视频 | 久久人人爽人人片av| 女性生殖器流出的白浆| 精品久久蜜臀av无| 老司机亚洲免费影院| 深夜精品福利| 女警被强在线播放| 在线 av 中文字幕| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频 | 好男人视频免费观看在线| 两个人免费观看高清视频| 日韩人妻精品一区2区三区| 男男h啪啪无遮挡| 久久毛片免费看一区二区三区| 极品人妻少妇av视频| 午夜两性在线视频| 亚洲欧美精品自产自拍| 国产精品一区二区免费欧美 | 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品在线美女| 满18在线观看网站| 人人妻人人添人人爽欧美一区卜| 日韩精品免费视频一区二区三区| 亚洲欧洲国产日韩| 亚洲专区国产一区二区| 亚洲男人天堂网一区| 久久久久视频综合| 九色亚洲精品在线播放| 美女脱内裤让男人舔精品视频| 久久人人爽av亚洲精品天堂| 国产成人系列免费观看| 男人操女人黄网站| 欧美日韩av久久| 美女高潮到喷水免费观看| 中文精品一卡2卡3卡4更新| 在线亚洲精品国产二区图片欧美| cao死你这个sao货| 好男人视频免费观看在线| www.999成人在线观看| 免费久久久久久久精品成人欧美视频| 成人国产一区最新在线观看 | 国产精品久久久久成人av| 一本一本久久a久久精品综合妖精| 精品一区在线观看国产| 国产精品麻豆人妻色哟哟久久| 99香蕉大伊视频| 欧美性长视频在线观看| 在线天堂中文资源库| 亚洲av电影在线观看一区二区三区| 日日爽夜夜爽网站| 国产一区有黄有色的免费视频| 亚洲美女黄色视频免费看| 高清欧美精品videossex| 91九色精品人成在线观看| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 久久综合国产亚洲精品| h视频一区二区三区| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 99精国产麻豆久久婷婷| 看免费av毛片| 国产精品国产av在线观看| 超碰成人久久| 国产午夜精品一二区理论片| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 最近最新中文字幕大全免费视频 | 久久精品人人爽人人爽视色| 国产精品秋霞免费鲁丝片| 妹子高潮喷水视频| 中国美女看黄片| a级毛片黄视频| 大码成人一级视频| 一级黄片播放器| 天堂中文最新版在线下载| 国产1区2区3区精品| 久久亚洲国产成人精品v| 曰老女人黄片| 尾随美女入室| 亚洲午夜精品一区,二区,三区| 国产精品二区激情视频| 七月丁香在线播放| 久久久久久久大尺度免费视频| 成人手机av| 99热全是精品| 一级片免费观看大全| 欧美性长视频在线观看| 国产成人av教育| 少妇粗大呻吟视频| 中文欧美无线码| 又紧又爽又黄一区二区| 久久久亚洲精品成人影院| 精品一品国产午夜福利视频| 在线av久久热| 丁香六月天网| 少妇人妻久久综合中文| 国产三级黄色录像| 欧美xxⅹ黑人| 狠狠精品人妻久久久久久综合| 性色av一级| 建设人人有责人人尽责人人享有的| av在线老鸭窝| 国产一区二区三区av在线| 欧美大码av| 在线亚洲精品国产二区图片欧美| 国产一区二区 视频在线| 精品亚洲乱码少妇综合久久| 一级毛片黄色毛片免费观看视频| 久久久精品94久久精品| 97在线人人人人妻| 欧美在线一区亚洲| 咕卡用的链子| 成年人免费黄色播放视频| 精品一区在线观看国产| 美女午夜性视频免费| 欧美亚洲 丝袜 人妻 在线| 最黄视频免费看| 国产精品一国产av| 视频区图区小说| 国产精品99久久99久久久不卡| 亚洲av日韩精品久久久久久密 | 国产日韩欧美亚洲二区| 又大又爽又粗| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 丰满少妇做爰视频| 99精国产麻豆久久婷婷| 午夜福利视频精品| 在线天堂中文资源库| 50天的宝宝边吃奶边哭怎么回事| 日韩 欧美 亚洲 中文字幕| 七月丁香在线播放| 午夜免费成人在线视频| 成人国产av品久久久| 免费日韩欧美在线观看| 亚洲一区二区三区欧美精品| bbb黄色大片| 日本a在线网址| 国产成人免费无遮挡视频| 国产精品一二三区在线看| 国产欧美日韩精品亚洲av| 亚洲国产欧美在线一区| 悠悠久久av| 免费看十八禁软件| 欧美 亚洲 国产 日韩一| 大香蕉久久网| 亚洲情色 制服丝袜| 国产精品麻豆人妻色哟哟久久| bbb黄色大片| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av | 9191精品国产免费久久| 亚洲成av片中文字幕在线观看| 在线av久久热| 国产在线一区二区三区精| 美女主播在线视频| 97人妻天天添夜夜摸| 在线看a的网站| 丁香六月欧美| 看十八女毛片水多多多| 国产99久久九九免费精品| 深夜精品福利| 成人影院久久| 18禁国产床啪视频网站| av网站在线播放免费| 国产又爽黄色视频| 成人国产av品久久久| av网站免费在线观看视频| 不卡av一区二区三区| 久久中文字幕一级| 老司机深夜福利视频在线观看 | videosex国产| 啦啦啦中文免费视频观看日本| 久久毛片免费看一区二区三区| √禁漫天堂资源中文www| 80岁老熟妇乱子伦牲交| 国产精品 欧美亚洲| 99久久人妻综合| 亚洲欧洲国产日韩| www.自偷自拍.com| 免费在线观看黄色视频的| 少妇粗大呻吟视频| 观看av在线不卡| 三上悠亚av全集在线观看| 久久久久久久久久久久大奶| 天天躁日日躁夜夜躁夜夜| 少妇的丰满在线观看| 国产日韩欧美亚洲二区| 亚洲欧洲国产日韩| 中文字幕制服av| 在线观看免费午夜福利视频| 999精品在线视频| 最黄视频免费看| 亚洲av综合色区一区| 午夜av观看不卡| 午夜视频精品福利| 欧美日韩福利视频一区二区| 好男人电影高清在线观看| 亚洲国产av影院在线观看| 久久久国产欧美日韩av| 最新在线观看一区二区三区 | 99国产精品免费福利视频| 亚洲欧美日韩高清在线视频 | 久久久精品94久久精品| 免费看不卡的av| 后天国语完整版免费观看| 国产一级毛片在线| 99热国产这里只有精品6| 韩国精品一区二区三区| 在线av久久热| 999久久久国产精品视频| 亚洲国产精品一区二区三区在线| 9热在线视频观看99| 麻豆av在线久日| 国产亚洲一区二区精品| 午夜av观看不卡| 亚洲av在线观看美女高潮| 性少妇av在线| 少妇精品久久久久久久| 少妇人妻 视频| 新久久久久国产一级毛片| 一本色道久久久久久精品综合| 精品亚洲成国产av| 欧美黄色淫秽网站|