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

    基于氣候適宜度的河南省夏玉米產(chǎn)量預(yù)報研究

    2015-06-27 08:53:26李樹巖余衛(wèi)東
    關(guān)鍵詞:省轄市夏玉米氣候

    李樹巖, 余衛(wèi)東

    (中國氣象局河南省農(nóng)業(yè)氣象保障與應(yīng)用技術(shù)重點實驗室,河南省氣象科學(xué)研究所,河南 鄭州 450003)

    ?

    基于氣候適宜度的河南省夏玉米產(chǎn)量預(yù)報研究

    李樹巖, 余衛(wèi)東

    (中國氣象局河南省農(nóng)業(yè)氣象保障與應(yīng)用技術(shù)重點實驗室,河南省氣象科學(xué)研究所,河南 鄭州 450003)

    從夏玉米不同發(fā)育階段所需的光、溫、水等生物學(xué)特性出發(fā),構(gòu)建了河南省夏玉米氣候適宜度模型,計算了1980—2007年全省13個主產(chǎn)省轄市夏玉米生長季綜合氣候適宜度;根據(jù)各階段氣候適宜度與夏玉米產(chǎn)量的相關(guān)關(guān)系,建立了以旬為步長的夏玉米產(chǎn)量預(yù)報模型,并利用2008—2012年數(shù)據(jù)進(jìn)行預(yù)報檢驗。結(jié)果表明,氣候適宜度指數(shù)與夏玉米產(chǎn)量呈顯著正相關(guān),氣候適宜度變化能較為客觀地反映夏玉米產(chǎn)量水平及其動態(tài)變化。1980—2007年各省轄市模型回代檢驗準(zhǔn)確率為87.5%~94.6%,全省平均為91.8%,2008—2012年各省轄市預(yù)報準(zhǔn)確率逐旬平均為92.3%~98.4%,全省平均為96.0%,能基本滿足農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)需求。

    氣候適宜度;產(chǎn)量預(yù)報;預(yù)報模型;夏玉米

    夏玉米是河南省重要的糧食作物,其產(chǎn)量波動將直接影響全年糧食生產(chǎn)的穩(wěn)定性。夏玉米產(chǎn)量豐歉受氣候、土壤、品種和農(nóng)業(yè)生產(chǎn)技術(shù)水平等多種因素影響,其中溫度、降水、光照等氣象條件貫穿整個生長過程,是最主要的環(huán)境因素,因此實時跟蹤,定量評價各生育階段氣象條件對作物的影響,是農(nóng)業(yè)氣象系列化服務(wù)的重要內(nèi)容[1]。隨著社會經(jīng)濟的發(fā)展,農(nóng)業(yè)生產(chǎn)對農(nóng)業(yè)氣象產(chǎn)量預(yù)報的需求日益增加,準(zhǔn)確預(yù)報夏玉米產(chǎn)量變化,為農(nóng)業(yè)生產(chǎn)提供預(yù)測信息,對保障糧食安全和農(nóng)業(yè)發(fā)展具有重要意義。20世紀(jì)70年代末就開始了作物產(chǎn)量預(yù)報的研究和服務(wù),取得了許多研究成果,建立了多種預(yù)報模式[2-5]。各模式多以統(tǒng)計方法為主,綜合了氣象、農(nóng)學(xué)、經(jīng)濟、遙感等多學(xué)科技術(shù)和方法[6],取得了較好的預(yù)報效果,隨著作物模擬模型技術(shù)的發(fā)展,基于作物模型的產(chǎn)量預(yù)報方法也開始逐步應(yīng)用,但模型中部分所需實時數(shù)據(jù)在業(yè)務(wù)中仍難以獲得,參數(shù)的調(diào)整也較為困難[7,8]。近年來,氣候適宜度理論不斷發(fā)展和完善[9,10],通過構(gòu)建光、溫、水綜合的氣候適宜度模型,開展作物產(chǎn)量預(yù)報的方法研究日趨成熟,魏瑞江等[11]通過計算河北省8個市夏玉米生育期內(nèi)逐旬氣候適宜度,建立了8個市夏玉米不同時段產(chǎn)量預(yù)報模型,試報準(zhǔn)確率滿足了業(yè)務(wù)服務(wù)需求。易雪等[12]將氣候適宜度指數(shù)在早稻產(chǎn)量動態(tài)預(yù)報上進(jìn)行了應(yīng)用。李曼華等[13]結(jié)合山東省小麥生理特性構(gòu)建適宜度模型,進(jìn)行冬小麥產(chǎn)量預(yù)報,單產(chǎn)預(yù)報準(zhǔn)確率在95%以上。劉偉昌等[14]建立了河南省冬小麥氣候適宜度模型,姜麗霞等[15]通過對玉米生長過程的干物質(zhì)累積曲線進(jìn)行逐時段的溫度和水分訂正,建立了黑龍江省玉米氣候產(chǎn)量的預(yù)測模式。但不論何種作物及產(chǎn)量預(yù)報方法,都有一定的區(qū)域針對性,不能將其他地區(qū)建立的預(yù)報模型照辦挪用,需根據(jù)本區(qū)域的氣候特點,綜合考慮光、溫、水氣象要素來構(gòu)建產(chǎn)量預(yù)報模型。夏玉米的生長發(fā)育和產(chǎn)量形成是一個動態(tài)變化的過程,需根據(jù)氣象條件變化進(jìn)行實時動態(tài)預(yù)報。因此,在前人研究的基礎(chǔ)上,本研究利用河南省夏玉米13個主產(chǎn)省轄市的產(chǎn)量資料和氣象資料,改進(jìn)氣候適宜度計算方法,構(gòu)建基于氣候適宜度的河南省夏玉米產(chǎn)量預(yù)報模型,實現(xiàn)逐旬滾動產(chǎn)量預(yù)報,全面、動態(tài)、客觀地揭示了夏玉米生產(chǎn)與氣象條件的關(guān)系,為糧食豐歉年景預(yù)測提供信息服務(wù)。

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

    1.1 資料選取

    根據(jù)河南省夏玉米播種面積及總產(chǎn),確定安陽、焦作、開封、洛陽、新鄉(xiāng)、鄭州、漯河、許昌、平頂山、周口、商丘、南陽和駐馬店13個省轄市作為主產(chǎn)區(qū)代表市,這13個省轄市夏玉米播種面積占全省的95%以上,可以代表全省夏玉米生產(chǎn)情況。氣象資料為上述13個省轄市1980—2012年逐旬平均氣溫、旬降水量和旬日照時數(shù)。產(chǎn)量資料為相應(yīng)年份的夏玉米單產(chǎn)資料,來源于河南省統(tǒng)計年鑒。

    1.2 氣象產(chǎn)量提取

    一般來說,農(nóng)作物產(chǎn)量可分為3個部分,即趨勢產(chǎn)量、氣象產(chǎn)量和隨機“噪聲”,表示為:

    Y=Yt+Yw+ε

    (1)

    式中:Y為作物歷年單產(chǎn);Yt為趨勢產(chǎn)量;Yw為氣象產(chǎn)量;ε為隨機“噪聲”,一般忽略不計,故(1)式可簡化為:

    Y=Yt+Yw

    (2)

    為消除各地產(chǎn)量水平的差異,重點考慮氣象條件對產(chǎn)量波動的影響,還需計算相對氣象產(chǎn)量Ywa,如式(3)

    Ywa=(Y-Yt)/Yt

    (3)

    本研究采用三次多項式和直線滑動平均2種方法來進(jìn)行趨勢產(chǎn)量的模擬,并計算氣象產(chǎn)量。利用趨勢產(chǎn)量與實際產(chǎn)量的擬合系數(shù),進(jìn)行2種模擬方法比較。

    1.2.1 三次多項式模擬 三次多項式模擬的商丘市夏玉米趨勢產(chǎn)量見圖1。趨勢產(chǎn)量序列方程為:

    商丘:Y=-0.041X3+4.586X2+19.49X+2 949

    式中Y為趨勢產(chǎn)量,X是從1980—2010年的年代序數(shù),根據(jù)方程可求得歷年的趨勢產(chǎn)量,同理建立其他省轄市的模擬方程。

    圖1 三次多項式模擬趨勢產(chǎn)量Fig.1 The method to simulate the tendency yield by cubic polynomial model

    1.2.2 直線滑動平均模擬 該方法將線性回歸模型與滑動平均相結(jié)合,將玉米產(chǎn)量的時間序列在某個階段內(nèi)的變化看做線性函數(shù),呈一直線。隨著階段的連續(xù)滑動,直線不斷變換位置,后延滑動,從而反映產(chǎn)量歷史演變趨勢變化。依次求取各階段內(nèi)的直線回歸模型,而各時間點上各直線滑動回歸模擬值的平均即為其趨勢產(chǎn)量[16]。設(shè)某階段的線性趨勢方程為;

    yi=ai+bit

    (4)

    式中:i為方程個數(shù),i=n-K+1,K為滑動步長,n為樣本序列個數(shù);t是時間序號。計算每個方程在t點上的函數(shù)值yi(t),每個t點上分別有q個函數(shù)值,q的多少與n,K有關(guān)。本研究中n=30,K取11。然后求算每個t點上q個函數(shù)的平均值:

    (5)

    1.2.3 2種模擬方法比較 用2種方法計算的實際產(chǎn)量與趨勢產(chǎn)量做擬合效果分析(表1),結(jié)果表明,二者的擬合效果差異不大,直線滑動平均方法擬合的相關(guān)系數(shù)全省平均為0.842,三次多項式方法為0.802,直線滑動平均方法稍好??紤]計算的簡便性,本研究采用三次多項式方法計算氣象產(chǎn)量。

    表1 直線滑動平均和三次多項式方法擬合系數(shù)比較Table 1 The fitting coefficients by the linear moving average model compared with cubic polynomial model

    2 夏玉米氣候適宜度模型建立

    2.1 溫度適宜度模型

    玉米喜溫,整個生育期均需要較高的溫度,不同生育階段適宜溫度范圍不同,過低過高都會對生長產(chǎn)生抑制作用。根據(jù)馬樹慶等[17]的研究成果,建立溫度適宜度模型:

    (6)

    式中:F(tij)省轄市為各點第j年第i旬溫度適宜度;tij為各點第j年第i旬平均氣溫,tL,tH,t0分別為夏玉米第i旬所需的旬平均最低氣溫、旬平均最高氣溫和旬平均適宜氣溫。tL,tH,t0的確定分別參照河南省夏玉米產(chǎn)區(qū)各市的旬平均最低氣溫、旬平均最高氣溫和旬平均氣溫的歷年平均值以及夏玉米各發(fā)育期對溫度的需求, 見表2。

    2.2 降水適宜度模型

    夏玉米生長在高溫和蒸發(fā)量大的夏季,一生需水量多,且不同生育階段對水分要求不同。一般抽雄前10 d到抽雄后20 d為需水臨界期,降水不足將嚴(yán)重影響夏玉米產(chǎn)量;降水量過多亦會發(fā)生玉米澇漬及連陰雨災(zāi)害,尤其在玉米生長中后期,直接影響開花授粉及籽粒灌漿,造成減產(chǎn)[18]。因此根據(jù)夏玉米各生育階段需水量是否滿足[19]計算降水適宜度。

    (7)

    表2 氣候適宜度模型參數(shù)Table 2 The parameters of climate suitability model

    注:E-表示1個月的前10 d,上旬;M-表示1個月中間10 d,中旬;L-表示1個月最后10 d,下旬;下同。

    Note: E- is the first ten-day of a month. M- is the middle ten-day of a month. L- is last ten-day of a month. The same as below.

    式中:△Wij為農(nóng)田水分盈虧值,ΔWij=Rij-Wij,Wij=Kc·Eij,Rij為第j年第i旬降水量,Wij為農(nóng)作物理論需水量,Kc為夏玉米作物系數(shù),Eij為該時段潛在蒸散量用彭曼公式計算所得。

    2.3 日照適宜度模型

    玉米喜光怕陰,充足的日照能促進(jìn)玉米高產(chǎn),夏玉米對光照最敏感的時段是雌穗分化期和開花吐絲期,如果此時光照不足,使玉米植株正常發(fā)育受阻或花絲、花粉活力降低造成空稈或結(jié)實不良。日照適宜度模型為:

    (8)

    2.4 綜合適宜度模型

    2.4.1 逐旬綜合適宜度模型 旬綜合適宜度為光、溫、水各氣象要素適宜度加權(quán)求和所得。

    F(Cij)省轄市=btiF(tij)省轄市+bRiF(Rij)省轄市+
    bsiF(Sij)省轄市

    (9)

    式中:F(Cij)省轄市為第j年第i旬綜合適宜度;bti,bRi,bsi分別為各旬溫度,降水和日照的適宜度權(quán)重系數(shù);bti,bRi,bsi的確定方法為將每旬各氣象要素適宜度與相應(yīng)年份的氣象產(chǎn)量做相關(guān)分析。用各要素相關(guān)系數(shù)除以3個要素的相關(guān)系數(shù)絕對值之和。因此bti+bRi+bsi=1,確定方法如下。

    (10)

    式中:αti,αRi,αsi為分別為各旬溫度、降水和日照時數(shù)適宜度與相對氣象產(chǎn)量的相關(guān)系數(shù);αci為各氣象要素相關(guān)系數(shù)絕對值之和。

    2.4.2 全生育期綜合適宜度模型 河南夏玉米全生育期共12旬,將各旬綜合氣候適宜度與相應(yīng)年份氣象產(chǎn)量做相關(guān)分析,用各旬相關(guān)系數(shù)除以12旬相關(guān)系數(shù)絕對值之和,計算各旬氣候適宜度權(quán)重系數(shù),然后累加計算全生育期氣候適宜度,如公式(11)

    (11)

    (12)

    F(Cj)省轄市為各省轄市第j年夏玉米全生育期綜合適宜度,Ki為各旬適宜度權(quán)重系數(shù),αi為各旬綜合氣候適宜度與氣象產(chǎn)量的相關(guān)系數(shù),n為夏玉米全生育期總旬?dāng)?shù),本研究n=12。

    2.4.3 全省氣候適宜度模型 全省氣候適宜度由各省轄市氣候適宜度加權(quán)求和所得,見公式(13)~(17)。

    他們在一起時總愛咬文嚼字,說些富有象征意義的暗語?!安祭侍睂λ馕吨∪?,對愛拉則是向人求助的女人。“你的文學(xué)午餐”是他用來指她不貞的說法,有時以風(fēng)趣的口吻說出,有時則一本正經(jīng)。(2014:204)

    (13)

    (14)

    (15)

    (16)

    (17)

    F(tij)省,F(xiàn)(Rij)省,F(xiàn)(Sij)省分別為第j年第i旬全省溫度、降水、日照適宜度。Qi為各省轄市適宜度影響權(quán)重系數(shù),由各省轄市夏玉米播種面積除以13個省轄市播種面積之和所得。權(quán)重系數(shù)計算結(jié)果見表3,F(xiàn)(Cij)省為第j年第i旬全省綜合適宜度,F(xiàn)(Cj)省為第j年全省夏玉米全生育期綜合氣候適宜度。

    表3 各省轄市種植面積權(quán)重系數(shù)Table 3 The weight coefficients in every city according to plant area

    2.5 氣候適宜度模型檢驗

    將13個省轄市1980—2007年氣候適宜度指數(shù)與對應(yīng)相對氣象產(chǎn)量做相關(guān)分析,如表4所示,除焦作市相關(guān)系數(shù)通過0.1顯著性檢驗外,其他省轄市均通過0.01的極顯著性檢驗,說明本研究建立的夏玉米氣候適宜度模型能較好地反映河南省氣候適宜性動態(tài)變化和產(chǎn)量增減趨勢,可以用來評價夏玉米生長適宜情況。結(jié)合各地的農(nóng)業(yè)基礎(chǔ)條件可知,焦作、新鄉(xiāng)等豫北地區(qū)灌溉條件較好,遇到干旱的年份可及時灌溉,一定程度上掩蓋了氣象條件對生產(chǎn)的影響,因此氣象產(chǎn)量與氣候適宜度的相關(guān)性較低。

    2.6 夏玉米產(chǎn)量動態(tài)預(yù)報模型

    河南夏玉米一般6月上旬播種,擬從7月中旬開始進(jìn)行逐旬滾動的預(yù)報模型,因此,計算6月上旬至相應(yīng)起報旬的累積適宜度指數(shù),建立各起報旬累積適宜度指數(shù)與相對氣象產(chǎn)量的回歸方程,進(jìn)行夏玉米產(chǎn)量滾動預(yù)報,如表5所示。x為相應(yīng)起報時刻的累積適宜度指數(shù),Y為預(yù)報所得相對氣象產(chǎn)量。對預(yù)報模型進(jìn)行F檢驗可知,除開封和周口7月中旬的預(yù)報方程未通過檢驗外,其他省轄市各旬產(chǎn)量預(yù)報方程具有顯著的統(tǒng)計意義,可用來預(yù)報夏玉米產(chǎn)量。

    表4 氣候適宜度指數(shù)與氣象產(chǎn)量的相關(guān)性檢驗Table 4 The correlation test between climate suitability index and meteorological yield

    表5 各起報時刻夏玉米產(chǎn)量預(yù)報擬合方程Table 5 The fitting equation of summer maize yield forecast in different forecasting date beginning

    注:*和**分別代表0.05和0.01的顯著水平。

    Note: * and ** mean the significance at 0.05 and 0.01 level respectively.

    3 夏玉米產(chǎn)量預(yù)報模型的檢驗

    計算夏玉米生育期內(nèi)綜合氣候適宜度,自7月中旬開始,進(jìn)行夏玉米氣象產(chǎn)量逐旬滾動預(yù)報。利用三次多項式擬合方法計算各省轄市趨勢產(chǎn)量,加上預(yù)報所得的相對氣象產(chǎn)量,獲得各地及全省的預(yù)報單產(chǎn)。預(yù)報產(chǎn)量與實際產(chǎn)量相比較,計算預(yù)報準(zhǔn)確率。

    3.1 單產(chǎn)預(yù)報檢驗

    3.1.1 單產(chǎn)預(yù)報回代檢驗 1980—2007年回代檢驗準(zhǔn)確率如表6所示,各省轄市旬平均回代檢驗準(zhǔn)確率在87.5%~94.6%,除駐馬店、南陽外,其他省轄市均在90%以上,其中安陽、焦作最高,回代檢驗準(zhǔn)確率分別為94.6%和94.0%。全省各旬平均回代檢驗準(zhǔn)確率在91.4%~92.0%(平均為91.8%)。

    統(tǒng)計各省轄市回代準(zhǔn)確率在90%以上的年份,計算年份所占比例,如表7所示,全省約七成的年份預(yù)報準(zhǔn)確率在90%以上。各省轄市之間差異較大,其中駐馬店僅46.9%的年份回代檢驗準(zhǔn)確率在90%以上,其次開封、南陽和周口市也較低,回代檢驗準(zhǔn)確率在90%的年份小于70%,安陽最高為86.2%。各旬平均在71.6%~73.3%(全省平均為72.7%)

    表6 1980—2007年回代檢驗準(zhǔn)確率Table 6 The accuracy rate of return test from 1980 to 2007

    表7 1980—2007年回代檢驗準(zhǔn)確率大于90%的年份比例Table 7 The proportion of years when the accuracy rate of return test is greater than 90% from 1980 to 2007 %

    3.1.2 單產(chǎn)試報檢驗 利用預(yù)報模型在2008—2012年產(chǎn)量預(yù)報中進(jìn)行應(yīng)用,預(yù)報準(zhǔn)確率如表8所示,各省轄市預(yù)報準(zhǔn)確率在92.3%~98.4%,除開封和駐馬店外,其他省轄市旬平均預(yù)報準(zhǔn)確率均>95%,鄭州最高為98.4%,全省平均為96.0%,基本能滿足業(yè)務(wù)服務(wù)的需求。

    3.2 產(chǎn)量趨勢預(yù)報檢驗

    根據(jù)當(dāng)年產(chǎn)量與上一年單產(chǎn)的增減情況進(jìn)行趨勢產(chǎn)量預(yù)報,如表9所示。結(jié)果表明,1980—2007年趨勢預(yù)報的回代準(zhǔn)確率為67%,除平頂山,許昌,周口3個省轄市較小外,其他省轄市均大于60%。2008—2012年趨勢預(yù)報準(zhǔn)確率平均約50%,5年中大部分省轄市只報對2年的增減趨勢,表明趨勢預(yù)報的準(zhǔn)確率還較低。

    4 結(jié)論與討論

    (1)從夏玉米不同發(fā)育階段上限溫度、最適溫度、下限溫度、需水量、需光性等生物學(xué)特性出發(fā),參考前人研究方法,建立了河南省夏玉米氣候適宜度模型,并根據(jù)河南夏玉米生長特點對模型參數(shù)進(jìn)行修訂。通過氣候適宜度與對應(yīng)氣象產(chǎn)量相關(guān)分析,2者呈極顯著的正相關(guān),表明氣候適宜度指數(shù)能較為客觀地反映夏玉米產(chǎn)量動態(tài)變化,可作為產(chǎn)量預(yù)報的驅(qū)動因子。

    表8 2008—2012年產(chǎn)量預(yù)報準(zhǔn)確率Table 8 The accuracy rate of yield forecasting from 2008 to 2012 %

    表9 產(chǎn)量趨勢預(yù)報檢驗Table 9 The test of tendency yield forecasting

    (2)計算各省轄市氣候適宜度指數(shù),自7月中旬起建立產(chǎn)量預(yù)報回歸模型,實現(xiàn)以旬為時間步長的產(chǎn)量動態(tài)預(yù)報。經(jīng)檢驗,各省轄市旬平均回代檢驗準(zhǔn)確率為87.5%~94.6%,全省平均91.8%;各省轄市預(yù)報準(zhǔn)確率為92.3%~98.4%,全省平均96.0%,能基本滿足業(yè)務(wù)服務(wù)需求。

    (3)應(yīng)用基于氣候適宜度的產(chǎn)量預(yù)報模型能達(dá)到較高的預(yù)報準(zhǔn)確率,但模型未考慮到重大災(zāi)害天氣對產(chǎn)量的影響,個別年份模型預(yù)報準(zhǔn)確率很低,這是基于統(tǒng)計方法的產(chǎn)量預(yù)報模型較為共性的問題[20,21],如2003年的周口、駐馬店、商丘、南陽等地,大范圍水澇災(zāi)害,當(dāng)年實際單產(chǎn)僅相當(dāng)于正常年份的20%~50%,預(yù)報準(zhǔn)確率較差。因此預(yù)報災(zāi)害對產(chǎn)量損失的影響的方法,還有待于進(jìn)一步探討。

    (4)用當(dāng)年產(chǎn)量與上一年相比的增減進(jìn)行產(chǎn)量趨勢預(yù)報,趨勢回代及預(yù)報檢驗準(zhǔn)確率約50%~85%,準(zhǔn)確率較低。表明基于適宜度方法的產(chǎn)量預(yù)報模型對產(chǎn)量增減的預(yù)報不敏感,預(yù)報精度有待于提高。

    [1] 陳懷亮,張雪芬. 玉米生產(chǎn)農(nóng)業(yè)氣象服務(wù)指南[M].北京:氣象出版社,1999.

    [2] 王建林,楊霏云,宋迎波. 西北地區(qū)玉米產(chǎn)量動態(tài)業(yè)務(wù)預(yù)報方法探討[J]. 應(yīng)用氣象學(xué)報,2004, 15(1) : 51-57.

    [3] 鄭昌玲, 楊霏云, 王建林, 等. 早稻產(chǎn)量動態(tài)預(yù)報模型[J].中國農(nóng)業(yè)氣象,2007, 28(4) : 412-416.

    [4] 王建林,趙四強. 全國棉花產(chǎn)量預(yù)報模式[J]. 氣象,1990,16(5):26-30.

    [5] 宋迎波,王建林,陳 暉,等. 中國油菜產(chǎn)量動態(tài)預(yù)報方法研究[J]. 氣象,2008,34(3):93-99.

    [6] 錢 拴,王建林. 農(nóng)業(yè)氣象作物產(chǎn)量預(yù)報的特點與思考[J]. 氣象科技,2003 (5): 33-38.

    [7] 黃晚華,薛昌穎,李忠輝. 基于作物生長模擬模型的產(chǎn)量預(yù)報方法研究進(jìn)展[J]. 中國農(nóng)業(yè)氣象,2009, 30(増1) : 140-143.

    [8] 帥細(xì)強,王石立,馬玉平,等. 基于水稻生長模型的氣象影響評價和產(chǎn)量動態(tài)預(yù)測[J]. 應(yīng)用氣象學(xué)報,2008, 19(1):71-81.

    [9] 趙 峰,千懷遂,焦士興. 農(nóng)作物氣候適宜度模型研究[J]. 資源科學(xué),2003,25(6):75-82.

    [10]冶明珠,郭建平,蔣躍林,等. 氣候變化對農(nóng)作物氣候適宜度影響研究進(jìn)展[J]. 安徽農(nóng)業(yè)科學(xué),2011,39(15):9104-9105,9134.

    [11]魏瑞江,宋迎波,王 鑫. 基于氣候適宜度的玉米產(chǎn)量動態(tài)預(yù)報方法[J]. 應(yīng)用氣象學(xué)報,2009,20(5):622-627.

    [12]易 雪,王健林,宋應(yīng)波. 氣候適宜度指數(shù)在早稻產(chǎn)量動態(tài)預(yù)報上的應(yīng)用[J]. 氣象,2010, 36(6):85-89.

    [13]李曼華,薛曉萍,李鴻怡. 基于氣候適宜度指數(shù)的山東省冬小麥產(chǎn)量動態(tài)預(yù)報[J]. 中國農(nóng)學(xué)通報,2012,28(12):291-295.

    [14]劉偉昌,陳懷亮,余衛(wèi)東,等.基于氣候適宜度指數(shù)的冬小麥動態(tài)產(chǎn)量預(yù)報技術(shù)研究[J]. 氣象與環(huán)境科學(xué),2008,31(2):21-24.

    [15]姜麗霞,王育光,孫孟梅,等. 黑龍江省玉米產(chǎn)量預(yù)報模式的研究[J]. 中國農(nóng)業(yè)氣象,2004, 25(1) : 13-16.

    [16]中國氣象局.QX/T 107—2009 小麥干旱災(zāi)害等級[S].北京:中國標(biāo)準(zhǔn)出版社,2007.

    [17]馬樹慶. 吉林省農(nóng)業(yè)氣候研究[M]. 北京:氣象出版社,1994:33.

    [18]孫景生,肖俊夫,段愛旺,等.夏玉米耗水規(guī)律及水分脅迫對其生長發(fā)育和產(chǎn)量的影響[J]. 玉米科學(xué),1999,7(2):45-48.

    [19]張 強,楊賢為,黃朝迎. 近30年氣候變化對黃土高原地區(qū)玉米生產(chǎn)潛力的影響[J]. 中國農(nóng)業(yè)氣象,1995,16(6) :19-24.

    [20]鄭昌玲,王建林,宋迎波,等. 大豆產(chǎn)量動態(tài)預(yù)報模型研究[J]. 大豆科學(xué),2008, 27(6):943-948.

    [21]游 超,蔡元剛,張玉芳. 基于氣象適宜指數(shù)的四川盆地水稻氣象產(chǎn)量動態(tài)預(yù)報技術(shù)研究[J]. 高原山地氣象研究,2011,31(1):51-55.

    (責(zé)任編輯:朱秀英)

    Research of summer maize yield forecasting based on climate suitability in Henan

    LI Shuyan, YU Weidong

    (China Meteorological Administration, Key Laboratory of Agrometeorological Safeguard and Applied Technique in Henan Province, Henan Institute of Meteorological Sciences, Zhengzhou 450003, China)

    The climate suitability model of summer maize in Henan was established based on the biological characteristics such as sunlight, temperature and water in different growth stages. The climate suitability of summer maize growing season in 13 cities in Henan was calculated from 1980 to 2007. According to the relationship between climate suitability and summer maize yield in different growth stages, the yield forecasting models by step of 10 days were established. At last the forecasting accuracy rates were tested using the data from 2008 to 2012. The results showed, there was a significant positive correlation between climate suitability and summer maize yield. The climate suitability could respond the level of summer maize yield and its dynamic change. From 1980 to 2007, the return test accuracy rates in every city were from 87.5% to 94.6%, with an average of 91.8%. From 2008 to 2012, the forecasting accuracy rates in every city were from 92.3 to 98.4%, with an average of 96.0%. It could meet basic needs for agro-meteorological operation and service.

    climate suitability; yield forecast; forecast model; summer maize

    1000-2340(2015)01-0027-08

    2014-01-17

    公益性行業(yè)(氣象)科研專項(GYHY201206022)

    李樹巖(1979-),女,河北唐山人,高級工程師,主要從事氣候資源利用與農(nóng)業(yè)減災(zāi)研究。

    S 513

    A

    猜你喜歡
    省轄市夏玉米氣候
    小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
    2018年河南省各省轄市及直管縣(市)專利申請量統(tǒng)計表(3月)
    河南科技(2018年12期)2018-09-10 05:12:39
    瞧,氣候大不同
    氣候變暖會怎樣?
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    2017年河南省各省轄市及直管縣(市)專利申請量統(tǒng)計表(2月)
    河南科技(2017年4期)2017-06-06 11:59:01
    2016年03月河南省各省轄市專利申請量統(tǒng)計表
    河南科技(2016年6期)2016-08-13 08:18:26
    立冬
    省轄市國民經(jīng)濟主要指標(biāo)(2013年)
    江蘇年鑒(2014年0期)2014-03-11 17:10:14
    氣候的未來掌握在我們手中
    99热这里只有是精品50| 天天躁日日操中文字幕| 2021天堂中文幕一二区在线观| 欧美潮喷喷水| 欧美精品啪啪一区二区三区| 亚洲成人久久性| 国产精品久久久久久亚洲av鲁大| 黄色一级大片看看| 亚洲久久久久久中文字幕| 欧美性猛交黑人性爽| 久久久成人免费电影| 婷婷精品国产亚洲av在线| 久久精品91蜜桃| 亚洲片人在线观看| 夜夜爽天天搞| 99热这里只有是精品在线观看 | 久久久久精品国产欧美久久久| 亚洲五月天丁香| 99riav亚洲国产免费| 国产精华一区二区三区| 一本一本综合久久| 欧美激情在线99| 国产精品日韩av在线免费观看| 深夜精品福利| 久久99热这里只有精品18| 亚洲自偷自拍三级| 偷拍熟女少妇极品色| 18美女黄网站色大片免费观看| 能在线免费观看的黄片| 欧美另类亚洲清纯唯美| 赤兔流量卡办理| 国产伦一二天堂av在线观看| 一本综合久久免费| 婷婷亚洲欧美| 亚洲,欧美,日韩| 国产精品野战在线观看| 亚洲第一欧美日韩一区二区三区| 在现免费观看毛片| 国产真实乱freesex| 国产精品久久视频播放| 亚洲精华国产精华精| 精品日产1卡2卡| 国产精品电影一区二区三区| 可以在线观看毛片的网站| 女人被狂操c到高潮| 99精品久久久久人妻精品| 我要看日韩黄色一级片| 丰满乱子伦码专区| www.999成人在线观看| 国产av不卡久久| 好男人在线观看高清免费视频| 老女人水多毛片| 亚洲色图av天堂| 1000部很黄的大片| 久久欧美精品欧美久久欧美| 亚洲第一区二区三区不卡| 免费无遮挡裸体视频| 国产在视频线在精品| 日本黄色片子视频| 直男gayav资源| 天天一区二区日本电影三级| 全区人妻精品视频| 757午夜福利合集在线观看| 亚洲三级黄色毛片| 男人狂女人下面高潮的视频| 色av中文字幕| 老女人水多毛片| 高潮久久久久久久久久久不卡| 内射极品少妇av片p| 精华霜和精华液先用哪个| 制服丝袜大香蕉在线| 亚洲人成伊人成综合网2020| 观看免费一级毛片| 一区福利在线观看| 亚洲精品456在线播放app | 丰满人妻一区二区三区视频av| 51午夜福利影视在线观看| 午夜精品久久久久久毛片777| 性欧美人与动物交配| 精品人妻熟女av久视频| 狠狠狠狠99中文字幕| 国产av不卡久久| 国产精品精品国产色婷婷| 一区二区三区四区激情视频 | 欧美在线一区亚洲| 亚州av有码| 国产精品av视频在线免费观看| 一级a爱片免费观看的视频| 国产高清激情床上av| 亚洲国产欧洲综合997久久,| 99国产综合亚洲精品| 免费观看的影片在线观看| 午夜激情福利司机影院| 老司机午夜十八禁免费视频| 欧美中文日本在线观看视频| 亚洲国产精品合色在线| 男人舔奶头视频| 亚洲成人中文字幕在线播放| av天堂在线播放| 成人特级黄色片久久久久久久| 99久久九九国产精品国产免费| 亚洲欧美日韩东京热| 亚洲最大成人中文| 亚洲一区二区三区色噜噜| 少妇的逼好多水| 欧美成狂野欧美在线观看| 精品99又大又爽又粗少妇毛片 | 精品午夜福利在线看| 18禁黄网站禁片免费观看直播| 国产不卡一卡二| 国产精品,欧美在线| 欧美一区二区国产精品久久精品| 欧美色欧美亚洲另类二区| 成人欧美大片| 女同久久另类99精品国产91| 波野结衣二区三区在线| 看黄色毛片网站| av在线老鸭窝| 欧美潮喷喷水| 91狼人影院| 三级国产精品欧美在线观看| 综合色av麻豆| bbb黄色大片| 最近最新中文字幕大全电影3| 成人精品一区二区免费| 亚洲无线观看免费| 黄色配什么色好看| 成人av一区二区三区在线看| 欧美又色又爽又黄视频| 国产午夜精品久久久久久一区二区三区 | 国产亚洲欧美98| 欧美成人免费av一区二区三区| 窝窝影院91人妻| 亚洲人成网站在线播放欧美日韩| 国产免费av片在线观看野外av| avwww免费| 亚洲av成人av| 亚洲自偷自拍三级| 国产一区二区激情短视频| 女生性感内裤真人,穿戴方法视频| 舔av片在线| 69av精品久久久久久| 老司机午夜十八禁免费视频| 精品不卡国产一区二区三区| or卡值多少钱| 尤物成人国产欧美一区二区三区| 国产精品影院久久| 欧美激情在线99| 又爽又黄无遮挡网站| 国产精品一区二区免费欧美| 最近最新中文字幕大全电影3| 男女那种视频在线观看| 国产 一区 欧美 日韩| 淫妇啪啪啪对白视频| 国产毛片a区久久久久| 亚洲经典国产精华液单 | or卡值多少钱| 91午夜精品亚洲一区二区三区 | 亚洲成人免费电影在线观看| 国产男靠女视频免费网站| 99国产精品一区二区蜜桃av| 国产大屁股一区二区在线视频| 永久网站在线| 亚洲国产精品成人综合色| 老女人水多毛片| 嫩草影视91久久| 51午夜福利影视在线观看| 国产大屁股一区二区在线视频| 在线观看av片永久免费下载| 观看美女的网站| 午夜久久久久精精品| 国产乱人视频| 国产精华一区二区三区| 麻豆一二三区av精品| 少妇人妻一区二区三区视频| 能在线免费观看的黄片| 成年版毛片免费区| a级毛片免费高清观看在线播放| 成人亚洲精品av一区二区| 日韩中文字幕欧美一区二区| 国内毛片毛片毛片毛片毛片| 午夜精品一区二区三区免费看| 国产午夜精品久久久久久一区二区三区 | 一区二区三区激情视频| 国产黄a三级三级三级人| 女人被狂操c到高潮| 成人特级黄色片久久久久久久| 噜噜噜噜噜久久久久久91| 日韩欧美在线二视频| 国产高清视频在线播放一区| 啪啪无遮挡十八禁网站| 可以在线观看毛片的网站| av国产免费在线观看| 久久久久九九精品影院| 中文资源天堂在线| 99在线人妻在线中文字幕| 精品午夜福利视频在线观看一区| 久久久久久九九精品二区国产| av视频在线观看入口| 高潮久久久久久久久久久不卡| 少妇人妻一区二区三区视频| 国产麻豆成人av免费视频| 亚洲一区二区三区不卡视频| 长腿黑丝高跟| 日韩欧美国产在线观看| 国产真实伦视频高清在线观看 | 桃色一区二区三区在线观看| 国产av麻豆久久久久久久| 全区人妻精品视频| 国产免费av片在线观看野外av| av女优亚洲男人天堂| 国产精品99久久久久久久久| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 亚洲专区中文字幕在线| 久久九九热精品免费| 99国产精品一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 人人妻人人澡欧美一区二区| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 少妇的逼好多水| 欧美另类亚洲清纯唯美| 午夜影院日韩av| 男人和女人高潮做爰伦理| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品国产亚洲| 午夜久久久久精精品| 日本与韩国留学比较| 性色avwww在线观看| 欧美+日韩+精品| 欧美黄色片欧美黄色片| 精品99又大又爽又粗少妇毛片 | 嫩草影院精品99| 国内久久婷婷六月综合欲色啪| 又粗又爽又猛毛片免费看| 天天躁日日操中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产av不卡久久| 欧美高清性xxxxhd video| 久久欧美精品欧美久久欧美| 嫁个100分男人电影在线观看| 91狼人影院| 亚洲国产欧洲综合997久久,| 黄色丝袜av网址大全| 国产三级中文精品| 日韩亚洲欧美综合| 国产精品久久久久久久电影| 久久久久久久精品吃奶| 亚洲,欧美,日韩| 九色成人免费人妻av| 国内久久婷婷六月综合欲色啪| 久久人人精品亚洲av| 久久九九热精品免费| 国产精品一区二区三区四区免费观看 | 亚洲成人久久性| 欧美成人a在线观看| 91九色精品人成在线观看| 亚洲国产欧洲综合997久久,| 久久人妻av系列| 国产一区二区在线观看日韩| 欧美中文日本在线观看视频| 国产av在哪里看| 日韩亚洲欧美综合| 亚洲 欧美 日韩 在线 免费| 午夜老司机福利剧场| 亚洲一区二区三区不卡视频| 国产一区二区三区在线臀色熟女| a级一级毛片免费在线观看| 久久精品国产亚洲av香蕉五月| 最新在线观看一区二区三区| 欧美+亚洲+日韩+国产| 禁无遮挡网站| 亚洲中文字幕日韩| 亚洲精品久久国产高清桃花| 精品福利观看| 亚洲aⅴ乱码一区二区在线播放| 小蜜桃在线观看免费完整版高清| 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 99久久精品热视频| 亚洲国产精品久久男人天堂| 精华霜和精华液先用哪个| 一级黄片播放器| 亚洲电影在线观看av| 久久亚洲真实| 精品福利观看| 国产精品久久久久久人妻精品电影| 亚洲无线观看免费| 国产精品精品国产色婷婷| 成人欧美大片| 有码 亚洲区| 久久人妻av系列| 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区| 91av网一区二区| 免费在线观看亚洲国产| 久久精品人妻少妇| 国产成人aa在线观看| 亚洲国产高清在线一区二区三| 毛片一级片免费看久久久久 | 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 国产黄片美女视频| 日本 欧美在线| 久久久久免费精品人妻一区二区| 成年女人永久免费观看视频| 久久久精品大字幕| www.999成人在线观看| 国产在视频线在精品| 中文字幕人成人乱码亚洲影| 老司机福利观看| 日韩成人在线观看一区二区三区| 久久99热6这里只有精品| 久久精品91蜜桃| 中亚洲国语对白在线视频| 757午夜福利合集在线观看| 亚洲狠狠婷婷综合久久图片| 小蜜桃在线观看免费完整版高清| 国产成人福利小说| 国产伦精品一区二区三区四那| 日本精品一区二区三区蜜桃| 亚洲人成网站高清观看| 精品乱码久久久久久99久播| 日本一本二区三区精品| 波多野结衣高清作品| 老女人水多毛片| 国产免费av片在线观看野外av| 成人亚洲精品av一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲av免费高清在线观看| 小说图片视频综合网站| 欧美成人一区二区免费高清观看| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 夜夜爽天天搞| 亚洲无线在线观看| 少妇的逼水好多| 国产一区二区亚洲精品在线观看| 亚洲中文日韩欧美视频| 舔av片在线| 午夜日韩欧美国产| 国产一区二区三区在线臀色熟女| 国产熟女xx| 亚洲性夜色夜夜综合| 欧美高清成人免费视频www| 欧美性猛交黑人性爽| 看片在线看免费视频| 天堂动漫精品| 丁香六月欧美| 亚洲狠狠婷婷综合久久图片| 色5月婷婷丁香| 亚洲最大成人av| 免费av观看视频| 成年女人永久免费观看视频| 18+在线观看网站| 一个人免费在线观看的高清视频| 欧美高清成人免费视频www| a级毛片a级免费在线| 美女被艹到高潮喷水动态| www.999成人在线观看| 午夜精品一区二区三区免费看| 亚洲自偷自拍三级| 99精品在免费线老司机午夜| 色哟哟·www| 天堂av国产一区二区熟女人妻| 国产一区二区激情短视频| 日日摸夜夜添夜夜添小说| 亚洲国产欧美人成| 色综合婷婷激情| 中文字幕熟女人妻在线| 内射极品少妇av片p| 精品免费久久久久久久清纯| 白带黄色成豆腐渣| 色精品久久人妻99蜜桃| 夜夜看夜夜爽夜夜摸| 99riav亚洲国产免费| 夜夜看夜夜爽夜夜摸| 国产亚洲av嫩草精品影院| 亚洲第一区二区三区不卡| 久久午夜亚洲精品久久| 亚洲成人久久性| 日韩有码中文字幕| 波多野结衣高清无吗| 在线免费观看不下载黄p国产 | 最近视频中文字幕2019在线8| 91在线观看av| 日日干狠狠操夜夜爽| 欧美不卡视频在线免费观看| 色精品久久人妻99蜜桃| 永久网站在线| 51国产日韩欧美| 嫩草影视91久久| 97超级碰碰碰精品色视频在线观看| 欧美日本亚洲视频在线播放| 亚洲av电影不卡..在线观看| 精品一区二区免费观看| 国内精品久久久久精免费| 欧美精品国产亚洲| 午夜免费男女啪啪视频观看 | 精品久久久久久久久av| 国产精品不卡视频一区二区 | 小说图片视频综合网站| 精品欧美国产一区二区三| 三级国产精品欧美在线观看| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 国产精品久久久久久久久免 | 欧美中文日本在线观看视频| 精品国产三级普通话版| 一本精品99久久精品77| 国产一区二区在线av高清观看| 国产av麻豆久久久久久久| 综合色av麻豆| 欧美日韩瑟瑟在线播放| 一区二区三区免费毛片| 99久久久亚洲精品蜜臀av| 在现免费观看毛片| 欧美不卡视频在线免费观看| 听说在线观看完整版免费高清| 成人欧美大片| 国产欧美日韩精品亚洲av| 日日摸夜夜添夜夜添小说| 舔av片在线| 伊人久久精品亚洲午夜| 一边摸一边抽搐一进一小说| 一进一出抽搐gif免费好疼| 精品无人区乱码1区二区| 日韩成人在线观看一区二区三区| 国产精品永久免费网站| 男女床上黄色一级片免费看| 免费人成在线观看视频色| 久久精品综合一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 色精品久久人妻99蜜桃| 12—13女人毛片做爰片一| x7x7x7水蜜桃| 男女下面进入的视频免费午夜| 天堂影院成人在线观看| 国产三级黄色录像| 伊人久久精品亚洲午夜| 国产成人影院久久av| 国产高潮美女av| 97超级碰碰碰精品色视频在线观看| 黄色丝袜av网址大全| 免费电影在线观看免费观看| 99热精品在线国产| av在线天堂中文字幕| 伊人久久精品亚洲午夜| 久久久久九九精品影院| 中文字幕av成人在线电影| aaaaa片日本免费| 婷婷色综合大香蕉| av在线老鸭窝| 嫩草影院精品99| 少妇高潮的动态图| 一本一本综合久久| 国产野战对白在线观看| 日本五十路高清| 白带黄色成豆腐渣| 久久久精品欧美日韩精品| 真实男女啪啪啪动态图| 国产在线男女| 欧美性感艳星| 久久精品国产亚洲av涩爱 | .国产精品久久| 99久久无色码亚洲精品果冻| 淫秽高清视频在线观看| 亚洲人成网站高清观看| 亚洲激情在线av| 中文在线观看免费www的网站| 97超视频在线观看视频| 欧美一级a爱片免费观看看| 亚洲自偷自拍三级| 中文字幕高清在线视频| 一二三四社区在线视频社区8| 看片在线看免费视频| 制服丝袜大香蕉在线| 一个人看的www免费观看视频| 午夜福利视频1000在线观看| 亚洲精品亚洲一区二区| 婷婷丁香在线五月| 黄片小视频在线播放| av在线老鸭窝| 亚洲最大成人av| 欧美精品啪啪一区二区三区| 十八禁人妻一区二区| 午夜福利高清视频| 老司机午夜福利在线观看视频| 久久这里只有精品中国| 免费在线观看日本一区| 最近视频中文字幕2019在线8| 99久久九九国产精品国产免费| 久久性视频一级片| 午夜福利在线在线| 欧美另类亚洲清纯唯美| 十八禁网站免费在线| 亚洲一区二区三区不卡视频| 18美女黄网站色大片免费观看| 身体一侧抽搐| 日本黄色片子视频| 香蕉av资源在线| 可以在线观看毛片的网站| 美女xxoo啪啪120秒动态图 | 午夜精品久久久久久毛片777| 欧美日本亚洲视频在线播放| 高清在线国产一区| 国产精品日韩av在线免费观看| 少妇熟女aⅴ在线视频| 一本久久中文字幕| 午夜福利在线在线| 国产在视频线在精品| 久久久久久国产a免费观看| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩高清专用| 精品不卡国产一区二区三区| 99国产精品一区二区蜜桃av| 欧美激情国产日韩精品一区| 有码 亚洲区| 九色国产91popny在线| 男女下面进入的视频免费午夜| 国产精品伦人一区二区| 亚洲,欧美精品.| 亚洲午夜理论影院| 亚洲精品久久国产高清桃花| 波多野结衣高清无吗| 亚洲国产精品久久男人天堂| 久久精品综合一区二区三区| 看十八女毛片水多多多| 成人性生交大片免费视频hd| 一a级毛片在线观看| 日韩大尺度精品在线看网址| а√天堂www在线а√下载| 精品国产亚洲在线| 女人被狂操c到高潮| 免费电影在线观看免费观看| 欧美日本视频| 午夜免费激情av| 日本五十路高清| 久久久久精品国产欧美久久久| 噜噜噜噜噜久久久久久91| 波野结衣二区三区在线| 热99re8久久精品国产| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 久久天躁狠狠躁夜夜2o2o| 亚洲电影在线观看av| 成人性生交大片免费视频hd| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 嫩草影视91久久| 色视频www国产| 婷婷六月久久综合丁香| 在线国产一区二区在线| 中国美女看黄片| 国产三级中文精品| 99久久成人亚洲精品观看| 露出奶头的视频| 在线a可以看的网站| 国产视频内射| 亚洲av美国av| 国产探花极品一区二区| 成人av在线播放网站| 少妇被粗大猛烈的视频| 国产激情偷乱视频一区二区| 久久亚洲真实| 国产精品久久久久久亚洲av鲁大| 午夜久久久久精精品| 美女xxoo啪啪120秒动态图 | 久久99热6这里只有精品| 欧美黄色片欧美黄色片| 久久精品国产亚洲av涩爱 | 日韩中字成人| 亚洲av中文字字幕乱码综合| 亚洲精品影视一区二区三区av| 神马国产精品三级电影在线观看| 精品久久久久久久久久久久久| 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 欧美黄色淫秽网站| 亚洲熟妇中文字幕五十中出| 亚洲av五月六月丁香网| 免费在线观看成人毛片| 国产精品一及| 欧美激情国产日韩精品一区| 成人性生交大片免费视频hd| 免费av毛片视频| 国产午夜精品论理片| 久久久国产成人免费| 一区福利在线观看| 18+在线观看网站| 美女cb高潮喷水在线观看| 国产一区二区在线观看日韩| 波多野结衣巨乳人妻| 日韩欧美精品免费久久 | 日本一本二区三区精品| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 又爽又黄a免费视频| 免费大片18禁| av在线蜜桃| 内地一区二区视频在线| 黄色视频,在线免费观看| 国产一区二区三区视频了| 欧美另类亚洲清纯唯美| 欧美性猛交╳xxx乱大交人| 两个人的视频大全免费| 日本三级黄在线观看| 欧美性猛交╳xxx乱大交人| 99久久久亚洲精品蜜臀av| 亚洲精品乱码久久久v下载方式| 国产精品女同一区二区软件 | 女生性感内裤真人,穿戴方法视频|