• <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
    氣候的未來掌握在我們手中
    日韩视频在线欧美| 一级黄片播放器| 成人二区视频| 十八禁网站网址无遮挡 | 色视频www国产| 一级爰片在线观看| 亚洲av二区三区四区| 日韩制服骚丝袜av| 最近中文字幕高清免费大全6| 精品人妻熟女av久视频| 成人美女网站在线观看视频| 亚洲中文av在线| 亚洲av国产av综合av卡| 黑人猛操日本美女一级片| 亚洲av免费高清在线观看| 中文字幕免费在线视频6| 国产精品人妻久久久久久| 我的老师免费观看完整版| 国产精品秋霞免费鲁丝片| 校园人妻丝袜中文字幕| 欧美成人精品欧美一级黄| 国产成人午夜福利电影在线观看| 亚洲真实伦在线观看| 午夜久久久在线观看| av不卡在线播放| 十八禁高潮呻吟视频 | 午夜老司机福利剧场| 黄色视频在线播放观看不卡| 美女福利国产在线| 亚洲不卡免费看| 国产成人一区二区在线| 免费高清在线观看视频在线观看| 国国产精品蜜臀av免费| 熟女人妻精品中文字幕| av专区在线播放| 国产免费福利视频在线观看| 久久久久视频综合| 久久久久网色| 99久久精品一区二区三区| 成年美女黄网站色视频大全免费 | 日日爽夜夜爽网站| 如何舔出高潮| 亚洲国产日韩一区二区| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 国产乱来视频区| 国产成人午夜福利电影在线观看| 免费高清在线观看视频在线观看| 最近中文字幕高清免费大全6| 成人二区视频| 久久av网站| 纵有疾风起免费观看全集完整版| 成人影院久久| 99热全是精品| 欧美老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 热re99久久国产66热| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 少妇丰满av| 91精品国产国语对白视频| 如何舔出高潮| 一级毛片黄色毛片免费观看视频| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的| 秋霞伦理黄片| 欧美亚洲 丝袜 人妻 在线| 日韩av免费高清视频| 曰老女人黄片| 久久久久久久亚洲中文字幕| 日韩人妻高清精品专区| 亚洲精品国产色婷婷电影| 人妻少妇偷人精品九色| 日韩精品免费视频一区二区三区 | 在线观看一区二区三区激情| 日韩欧美 国产精品| 日韩三级伦理在线观看| 国产中年淑女户外野战色| 亚洲人成网站在线观看播放| 欧美亚洲 丝袜 人妻 在线| 9色porny在线观看| 热re99久久精品国产66热6| a级片在线免费高清观看视频| 亚洲精品久久久久久婷婷小说| 插阴视频在线观看视频| 国产视频内射| 亚洲在久久综合| av在线观看视频网站免费| av天堂久久9| 国产精品99久久99久久久不卡 | 免费观看av网站的网址| 多毛熟女@视频| 亚洲欧美日韩另类电影网站| 久久亚洲国产成人精品v| 自拍偷自拍亚洲精品老妇| 九九爱精品视频在线观看| 日韩成人伦理影院| 日韩不卡一区二区三区视频在线| 老司机影院成人| 久热这里只有精品99| 边亲边吃奶的免费视频| 精品一区在线观看国产| 少妇熟女欧美另类| 午夜福利,免费看| 成人黄色视频免费在线看| 啦啦啦视频在线资源免费观看| 成年av动漫网址| 国产高清不卡午夜福利| 在线观看三级黄色| 国产69精品久久久久777片| 国产一区二区三区综合在线观看 | 日韩一本色道免费dvd| 日韩强制内射视频| 人妻系列 视频| 亚洲精品aⅴ在线观看| 大码成人一级视频| 亚洲国产精品成人久久小说| 欧美激情国产日韩精品一区| 亚洲综合色惰| 日韩 亚洲 欧美在线| 日韩精品免费视频一区二区三区 | 少妇的逼水好多| 只有这里有精品99| 美女内射精品一级片tv| 国产午夜精品一二区理论片| 五月伊人婷婷丁香| 好男人视频免费观看在线| 乱码一卡2卡4卡精品| 日本黄大片高清| 一本大道久久a久久精品| 精品亚洲成a人片在线观看| 国产又色又爽无遮挡免| 亚洲欧美精品专区久久| 夜夜爽夜夜爽视频| 亚洲国产精品一区二区三区在线| 亚洲国产精品一区三区| 蜜桃久久精品国产亚洲av| 黑丝袜美女国产一区| 亚洲精品色激情综合| 亚洲国产精品成人久久小说| 国产毛片在线视频| 一级毛片 在线播放| 简卡轻食公司| 国产在线视频一区二区| 久久人人爽人人爽人人片va| 国产亚洲一区二区精品| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频 | 纯流量卡能插随身wifi吗| 久久女婷五月综合色啪小说| 中国三级夫妇交换| 国产精品久久久久久精品古装| 看非洲黑人一级黄片| 亚洲内射少妇av| 热re99久久精品国产66热6| 人人妻人人澡人人爽人人夜夜| 成人国产av品久久久| 男人狂女人下面高潮的视频| 97超视频在线观看视频| 99九九线精品视频在线观看视频| 成人影院久久| 国语对白做爰xxxⅹ性视频网站| 99热这里只有是精品在线观看| 久久久久久久久久久丰满| 国产精品人妻久久久久久| 大陆偷拍与自拍| 日韩av不卡免费在线播放| 亚洲av福利一区| 男男h啪啪无遮挡| 热re99久久精品国产66热6| 久久久精品免费免费高清| 欧美日韩一区二区视频在线观看视频在线| 少妇人妻精品综合一区二区| 一区二区av电影网| 亚洲精品国产成人久久av| 久久精品夜色国产| 免费av不卡在线播放| 久久久久视频综合| 国产91av在线免费观看| 永久免费av网站大全| 久久久久久久精品精品| 久久久国产欧美日韩av| 日本av免费视频播放| 热re99久久国产66热| 成人毛片a级毛片在线播放| 国产亚洲午夜精品一区二区久久| 成年av动漫网址| 人妻夜夜爽99麻豆av| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| 亚洲一区二区三区欧美精品| 国语对白做爰xxxⅹ性视频网站| 国产成人午夜福利电影在线观看| 日韩不卡一区二区三区视频在线| 国产精品一二三区在线看| 亚洲经典国产精华液单| 国产精品欧美亚洲77777| 三上悠亚av全集在线观看 | 日日爽夜夜爽网站| 亚洲人成网站在线观看播放| 一本色道久久久久久精品综合| 熟女人妻精品中文字幕| 九九在线视频观看精品| 欧美性感艳星| 美女中出高潮动态图| 国产精品无大码| 国语对白做爰xxxⅹ性视频网站| 18+在线观看网站| 亚洲国产最新在线播放| 国产亚洲av片在线观看秒播厂| 亚洲一区二区三区欧美精品| 日本欧美国产在线视频| 欧美最新免费一区二区三区| av一本久久久久| 亚洲人成网站在线播| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区成人| 久久国产精品男人的天堂亚洲 | 亚洲精品aⅴ在线观看| 国产av国产精品国产| 如日韩欧美国产精品一区二区三区 | 亚洲精品色激情综合| 黄色视频在线播放观看不卡| 最近中文字幕2019免费版| 男女国产视频网站| 蜜桃久久精品国产亚洲av| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 久久鲁丝午夜福利片| 国产国拍精品亚洲av在线观看| 在线亚洲精品国产二区图片欧美 | 2018国产大陆天天弄谢| 精品人妻一区二区三区麻豆| 汤姆久久久久久久影院中文字幕| 老司机亚洲免费影院| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 久久精品国产鲁丝片午夜精品| 国产黄片视频在线免费观看| 亚洲美女搞黄在线观看| freevideosex欧美| 亚洲欧美中文字幕日韩二区| 人体艺术视频欧美日本| 91精品伊人久久大香线蕉| 久久久精品免费免费高清| 亚洲怡红院男人天堂| 在线免费观看不下载黄p国产| 亚洲av二区三区四区| 久久影院123| 免费黄频网站在线观看国产| 边亲边吃奶的免费视频| 五月伊人婷婷丁香| 国产视频内射| 日韩一本色道免费dvd| 精品久久久精品久久久| 亚洲在久久综合| 男女边吃奶边做爰视频| 久久99蜜桃精品久久| 最后的刺客免费高清国语| 日日摸夜夜添夜夜添av毛片| 午夜久久久在线观看| 亚洲国产精品国产精品| 五月伊人婷婷丁香| 一级av片app| 边亲边吃奶的免费视频| 国内精品宾馆在线| 夫妻性生交免费视频一级片| 91在线精品国自产拍蜜月| av专区在线播放| 七月丁香在线播放| 亚洲国产欧美日韩在线播放 | 青春草国产在线视频| 亚洲av日韩在线播放| 亚洲高清免费不卡视频| 99国产精品免费福利视频| 中文精品一卡2卡3卡4更新| 午夜精品国产一区二区电影| 亚洲真实伦在线观看| 亚洲国产精品成人久久小说| 日本黄大片高清| 国产日韩欧美亚洲二区| 久久久久久久久久久久大奶| 精品一区二区三区视频在线| 高清午夜精品一区二区三区| 久久国产亚洲av麻豆专区| kizo精华| 夫妻午夜视频| 日本欧美视频一区| av免费在线看不卡| 乱人伦中国视频| 69精品国产乱码久久久| 久久精品久久久久久久性| 我的老师免费观看完整版| 午夜老司机福利剧场| 黄色配什么色好看| 日本爱情动作片www.在线观看| 国产伦精品一区二区三区四那| 亚洲国产精品专区欧美| 老熟女久久久| 亚洲美女视频黄频| 丝瓜视频免费看黄片| 哪个播放器可以免费观看大片| 国产美女午夜福利| 日本免费在线观看一区| 好男人视频免费观看在线| 黄色配什么色好看| 国精品久久久久久国模美| 人妻一区二区av| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 另类精品久久| freevideosex欧美| 国产乱来视频区| 日韩人妻高清精品专区| 亚洲人与动物交配视频| 制服丝袜香蕉在线| 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 日本91视频免费播放| 国产精品成人在线| 秋霞在线观看毛片| 一区二区三区精品91| 欧美日韩av久久| 在线播放无遮挡| 插逼视频在线观看| 国产在视频线精品| 天堂中文最新版在线下载| 国产成人a∨麻豆精品| 国产成人免费观看mmmm| 国产69精品久久久久777片| 国产黄色免费在线视频| 欧美日韩av久久| 国产精品一区二区性色av| 国模一区二区三区四区视频| av福利片在线观看| 人妻人人澡人人爽人人| 你懂的网址亚洲精品在线观看| videossex国产| 国产精品.久久久| 人人妻人人添人人爽欧美一区卜| 高清视频免费观看一区二区| 偷拍熟女少妇极品色| 91精品一卡2卡3卡4卡| 国产精品三级大全| 日韩一区二区视频免费看| 国产精品人妻久久久影院| 久久精品国产鲁丝片午夜精品| 少妇人妻 视频| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 国产69精品久久久久777片| 免费人妻精品一区二区三区视频| 另类精品久久| 国产精品无大码| 在线观看www视频免费| 亚洲精品日韩在线中文字幕| 午夜激情福利司机影院| 自拍偷自拍亚洲精品老妇| 麻豆精品久久久久久蜜桃| 日产精品乱码卡一卡2卡三| 黄色视频在线播放观看不卡| 男女啪啪激烈高潮av片| 免费黄色在线免费观看| 汤姆久久久久久久影院中文字幕| 欧美3d第一页| 26uuu在线亚洲综合色| 国产精品久久久久久久久免| 观看av在线不卡| 男人狂女人下面高潮的视频| 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美 | 视频中文字幕在线观看| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 亚洲国产精品999| av天堂久久9| 久久毛片免费看一区二区三区| 男女国产视频网站| 亚洲,欧美,日韩| 丰满乱子伦码专区| 亚洲,欧美,日韩| 欧美日韩视频高清一区二区三区二| a级片在线免费高清观看视频| 成人18禁高潮啪啪吃奶动态图 | 女人久久www免费人成看片| 91久久精品国产一区二区成人| 久久久国产精品麻豆| 91精品国产国语对白视频| 亚洲经典国产精华液单| 春色校园在线视频观看| 国产成人精品久久久久久| 国内精品宾馆在线| 男女啪啪激烈高潮av片| 亚洲成人一二三区av| 日本午夜av视频| 嘟嘟电影网在线观看| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 最新中文字幕久久久久| 91成人精品电影| 亚洲自偷自拍三级| 午夜视频国产福利| 久久影院123| 一个人看视频在线观看www免费| 在线观看美女被高潮喷水网站| 国产成人freesex在线| 一二三四中文在线观看免费高清| 欧美最新免费一区二区三区| 国产亚洲精品久久久com| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 国产极品天堂在线| 免费看日本二区| 婷婷色综合大香蕉| 涩涩av久久男人的天堂| 精品午夜福利在线看| 久久国产精品男人的天堂亚洲 | 亚洲av男天堂| 男女无遮挡免费网站观看| 久久久久久久久久人人人人人人| 91久久精品国产一区二区成人| 夫妻性生交免费视频一级片| 热99国产精品久久久久久7| 在线看a的网站| 麻豆成人午夜福利视频| 自拍偷自拍亚洲精品老妇| 欧美成人精品欧美一级黄| 九九久久精品国产亚洲av麻豆| 国产乱人偷精品视频| xxx大片免费视频| 久久狼人影院| 9色porny在线观看| 亚洲经典国产精华液单| 老司机影院成人| 午夜日本视频在线| 国产欧美日韩一区二区三区在线 | 男男h啪啪无遮挡| 你懂的网址亚洲精品在线观看| 男人爽女人下面视频在线观看| 亚洲av国产av综合av卡| 18禁动态无遮挡网站| 亚洲美女视频黄频| 成人综合一区亚洲| 国产亚洲最大av| 高清视频免费观看一区二区| 好男人视频免费观看在线| 国产成人精品久久久久久| www.av在线官网国产| 日本爱情动作片www.在线观看| 亚洲av二区三区四区| 久久国内精品自在自线图片| 精品人妻偷拍中文字幕| 热re99久久精品国产66热6| av福利片在线| 99热这里只有是精品50| 久久影院123| 日本黄色日本黄色录像| 夜夜骑夜夜射夜夜干| 天天操日日干夜夜撸| 国产精品一区www在线观看| 国产精品国产三级专区第一集| 国产精品蜜桃在线观看| 一级毛片 在线播放| 日韩一区二区视频免费看| 欧美性感艳星| 久久久久久久久久人人人人人人| 人妻少妇偷人精品九色| 日本午夜av视频| 丁香六月天网| 成人影院久久| 国产成人freesex在线| 亚洲国产最新在线播放| 男女边摸边吃奶| 国产精品国产三级国产专区5o| 日韩中字成人| 精品久久久噜噜| 午夜久久久在线观看| 99国产精品免费福利视频| 精品一区二区免费观看| 我要看黄色一级片免费的| 肉色欧美久久久久久久蜜桃| 国产日韩欧美视频二区| 亚洲一区二区三区欧美精品| 狠狠精品人妻久久久久久综合| 黑人高潮一二区| 久久婷婷青草| 国产高清有码在线观看视频| 制服丝袜香蕉在线| 妹子高潮喷水视频| 少妇被粗大猛烈的视频| 老女人水多毛片| 午夜久久久在线观看| 不卡视频在线观看欧美| 日韩欧美精品免费久久| 精品人妻偷拍中文字幕| 国产精品蜜桃在线观看| 欧美精品国产亚洲| 十八禁网站网址无遮挡 | 寂寞人妻少妇视频99o| 秋霞在线观看毛片| 菩萨蛮人人尽说江南好唐韦庄| 日韩一区二区视频免费看| 国产高清有码在线观看视频| 男女免费视频国产| 欧美日韩av久久| 精品久久久久久久久亚洲| 久久久a久久爽久久v久久| 大片电影免费在线观看免费| 国产成人freesex在线| 国产亚洲精品久久久com| 一区二区三区精品91| 国产日韩欧美亚洲二区| 成人18禁高潮啪啪吃奶动态图 | 看免费成人av毛片| 91在线精品国自产拍蜜月| 日本黄色日本黄色录像| 日韩精品免费视频一区二区三区 | 久久青草综合色| 高清视频免费观看一区二区| 色视频在线一区二区三区| 久久毛片免费看一区二区三区| 有码 亚洲区| 久热这里只有精品99| 日韩三级伦理在线观看| 国产精品麻豆人妻色哟哟久久| 国产精品国产三级国产av玫瑰| 亚洲av不卡在线观看| 国产极品粉嫩免费观看在线 | 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 日韩亚洲欧美综合| 精品熟女少妇av免费看| 免费观看的影片在线观看| 三级国产精品欧美在线观看| 日韩亚洲欧美综合| 在线看a的网站| 亚洲欧美中文字幕日韩二区| a级毛片在线看网站| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 久久人人爽人人片av| 日日啪夜夜撸| 国产精品国产三级国产专区5o| freevideosex欧美| 十八禁高潮呻吟视频 | 国产一区二区三区综合在线观看 | 18禁在线无遮挡免费观看视频| 日韩av在线免费看完整版不卡| 精品99又大又爽又粗少妇毛片| 国产亚洲午夜精品一区二区久久| 亚洲av二区三区四区| 亚洲欧美日韩另类电影网站| 国产一区二区在线观看av| 午夜视频国产福利| 大片免费播放器 马上看| 亚洲av国产av综合av卡| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频| 一级毛片黄色毛片免费观看视频| 亚洲av综合色区一区| xxx大片免费视频| 日韩,欧美,国产一区二区三区| 在线 av 中文字幕| 人人妻人人澡人人看| 人人妻人人看人人澡| 一区二区三区乱码不卡18| 欧美日韩在线观看h| 夜夜爽夜夜爽视频| 久久国产亚洲av麻豆专区| 亚洲精品亚洲一区二区| 国产乱人偷精品视频| 偷拍熟女少妇极品色| 久久99一区二区三区| 晚上一个人看的免费电影| 亚洲欧洲日产国产| 亚洲av二区三区四区| 建设人人有责人人尽责人人享有的| 亚州av有码| 国产欧美亚洲国产| 97精品久久久久久久久久精品| 99久久精品热视频| 麻豆成人av视频| 日本欧美国产在线视频| 最新中文字幕久久久久| 日韩在线高清观看一区二区三区| 国语对白做爰xxxⅹ性视频网站| 如日韩欧美国产精品一区二区三区 | 免费黄色在线免费观看| 十分钟在线观看高清视频www | 日韩强制内射视频| 免费观看av网站的网址| 国产欧美日韩一区二区三区在线 | 麻豆精品久久久久久蜜桃| 亚洲自偷自拍三级| 精品酒店卫生间| 国产精品一区www在线观看| 国产精品秋霞免费鲁丝片| 人妻一区二区av| av女优亚洲男人天堂| 国产精品无大码| 精品人妻一区二区三区麻豆| 夫妻午夜视频| 夫妻性生交免费视频一级片| 亚洲av综合色区一区| 亚洲自偷自拍三级| 日日摸夜夜添夜夜添av毛片| 亚洲成人手机| 亚洲av国产av综合av卡| 日本欧美视频一区| 美女内射精品一级片tv| 欧美日韩视频精品一区| 赤兔流量卡办理| 久久国内精品自在自线图片| 亚洲欧洲国产日韩| 亚洲欧美成人精品一区二区|