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

    基于高斯函數(shù)的春玉米葉片功能期模型構(gòu)建與應用

    2019-07-23 08:41:00李姚姚范盼盼王春霞王克如謝瑞芝李少昆
    作物學報 2019年8期
    關(guān)鍵詞:葉位高斯葉片

    李姚姚 范盼盼 明 博 王春霞 王克如 侯 鵬 謝瑞芝 李少昆

    ?

    基于高斯函數(shù)的春玉米葉片功能期模型構(gòu)建與應用

    李姚姚 范盼盼 明 博 王春霞 王克如 侯 鵬 謝瑞芝*李少昆*

    中國農(nóng)業(yè)科學院作物科學研究所 / 農(nóng)業(yè)部作物生理生態(tài)重點實驗室, 北京 100081

    玉米; 葉片功能期; 高斯函數(shù); 生長模型

    玉米葉片是光合作用的主要場所, 是物質(zhì)生產(chǎn)和產(chǎn)量形成的基地[1-2]。前人將葉片生長發(fā)育分為展開、持續(xù)和衰老3個階段[3], 從葉片完全展開至1/2葉片變黃所持續(xù)的天數(shù)為葉片功能期[4], 是綠葉面積光合有效持續(xù)期, 可以作為評價葉片光合作用時間的關(guān)鍵指標, 對評估葉片光合生產(chǎn)能力具有重要作用。

    玉米各節(jié)位葉片生長有一定的規(guī)律性和順序性。前人認為玉米葉片展開順序按葉序自下而上依次展開, 玉米葉片衰老以基部葉為先, 而后由下部葉轉(zhuǎn)向上部葉再發(fā)展到中部葉[4]。玉米葉片展開與衰老順序的不同造成葉片功能期大小隨葉位上升呈先增加后減小的趨勢[5], 按其功能期大小將葉片劃分為葉組, 表現(xiàn)為中部葉組功能期最長(60~63 d), 上部葉組其次(45 d), 下部葉組再次(33~60 d), 基部葉組最短(18~32 d)[6]。品種遺傳特性[7-8], 缺氮[9-10]、增密[11-12]、干旱[13-16]、高溫[11,15]、冷害[13]、病蟲害[13]等逆境脅迫均會影響葉片功能期, 甚至造成早衰, 但對其變化規(guī)律及影響因素的研究多為定性描述。以數(shù)學模型量化表達已普遍應用于作物生長、發(fā)育及物質(zhì)生產(chǎn)[17-21]等相關(guān)研究, 關(guān)于整株葉面積隨生育進程變化的量化研究也有涉及[3,22-23], 但關(guān)于量化評估不同葉位葉片功能期的報道較少。2003年Lizaso等[3]曾提出用高斯函數(shù)對不同葉位功能期(GDD8)進行量化分析, 并建立葉片功能期模型, 2012年Kim等[24]借鑒Lizaso等[3]的方法, 將高斯模型嵌入MAIZISM模型中, 由于該模型構(gòu)建需要植株全部葉片的功能期, 費時耗力, 關(guān)于此函數(shù)模型的本地化研究和應用較少。本研究利用2015—2017年先玉335和鄭單958的葉片功能期數(shù)據(jù), 分析以高斯函數(shù)構(gòu)建的玉米葉片功能期模型的科學性和適用性, 探索簡化模型構(gòu)建、量化評估葉片功能期的技術(shù)與方法, 為相關(guān)研究提供借鑒。

    1 材料與方法

    1.1 試驗區(qū)概況

    中國農(nóng)業(yè)科學院作物科學研究所吉林公主嶺綜合試驗站(43°11′–44°9′N, 124°02′–125°18′E)位于吉林省中西部, 屬于北溫帶半濕潤大陸性季風氣候, 年平均降雨量594.8 mm, 全年無霜期144 d, 累積積溫(≥10℃) 2600~2800℃, 土壤類型為黑鈣土。

    1.2 試驗設計

    試驗于2015—2017年進行, 生育期內(nèi)氣溫變化見圖1。試驗小區(qū)2009—2017年長期不施氮, 土壤耕層含有機質(zhì)29.68 g kg–1、全氮0.99 g kg–1、速效磷28.16 mg kg–1、速效鉀240.28 g kg–1。選用品種先玉335 (XY335, PH6WC♀×PH4CV♂)和鄭單958 (ZD958, 鄭58♀×昌7-2♂), 種植密度為67,500株 hm–2, 行距為65 cm, 播前底肥配施過磷酸鈣(含P2O542.5 kg hm–2)、氯化鉀(含K2O 42.5 kg hm–2), 除草、噴藥等管理同大田生產(chǎn), 生育進程見表1。

    圖1 2015?2017年4月至9月最高、最低、平均氣溫圖

    灰色區(qū)域表示氣溫變化范圍, 黑色線表示平均氣溫。

    The gray area represents the temperature range, and the black line represents the average temperature.

    表1 2015?2017年生育進程

    1.3 測定項目與方法

    玉米出苗期選擇生長一致植株標記, 四葉期篩選出生長狀況穩(wěn)定的20個植株掛牌標記并標注葉片所在葉位, 九葉期確定長期觀察的5個植株, 記錄植株每片葉的展開與衰老時間。

    葉片完全展開的標準[25-26], 以葉環(huán)從上一展開葉基部露出, 兩葉葉環(huán)平齊, 葉耳不再包裹莖稈, 葉片平展并與莖稈形成一定夾角時為完全展開日期, 當綠葉面積下降至最大葉面積50%時為衰老日期[25-26], 以葉片完全展開日期與衰老日期內(nèi)的GDD8表示葉片功能期[3]。

    1.4 數(shù)據(jù)處理與分析

    采用Microsoft Excel 2010進行數(shù)據(jù)分析處理和作圖, 采用Curve Expert Professional 2.2.0進行曲線擬合。

    2 結(jié)果與分析

    2.1 玉米葉片展開時間與衰老時間

    如圖2所示, 以出苗為起點, 不同葉位葉片的展開時間隨著葉位的上升順序延遲, 上部葉位葉片展開時間比較接近; 葉片衰老時間則隨葉位上升先增大后減小, XY335葉片衰老最晚出現(xiàn)在18葉, ZD958衰老最晚出現(xiàn)在19葉, 衰老發(fā)生最晚的葉位均高于穗位葉, 品種間和年際間變化規(guī)律一致。

    2.2 葉片功能期模型的擬合與驗證

    圖2 不同葉位葉片展開與衰老時間

    圖A為先玉335葉片展開時間與衰老時間, 圖B為鄭單958葉片展開時間與衰老時間。實心圓表示出苗到葉片展開的天數(shù), 空心三角表示出苗到葉片衰老的天數(shù), 虛線表示穗位葉所在。

    Fig. A shows the leaf expanded and senescence time of Xianyu 335, and Fig. B shows the leaf expanded and senescence time of Zhengdan 958. The solid circle represents the days from sowing to leaf expanded, the hollow circle represents the days from sowing to leaf senescence, and the dashed line represents the position of ear leaf.

    圖3 2015年和2016年葉片功能期擬合圖

    虛線表示穗位葉所在, XY335表示先玉335, ZD958表示鄭單958。

    The dashed line represents the position of ear leaf. XY335: Xianyu 335; ZD958: Zhengdan 958.

    利用2017年實測數(shù)值分別對2個品種的擬合方程驗證顯示(圖4), XY335擬合兩者的歸一化均方根誤差NRMSE為8.2%, ZD958擬合兩者的歸一化均方根誤差NRMSE為6.3%, 表明利用高斯函數(shù)構(gòu)建的葉片功能期模型具有較好的預測性和適用性。

    Lizaso等[3]認為, 基于高斯函數(shù)的葉片功能期模型參數(shù)具有農(nóng)學生理意義, 其中+值代表最大功能期,值即最大功能期所對應葉位,值代表曲線的陡度。比較本研究2個品種的模型參數(shù)可以發(fā)現(xiàn), ZD958的最大功能期、曲線陡度大于XY335, 最大功能期所在葉位取整后均為16, 與實測值的吻合度較高, 品種間差異也可以通過模型參數(shù)值區(qū)分(表2)。

    圖4 模型擬合值與實測值比較

    XY335表示先玉335, ZD958表示鄭單958。XY335: Xianyu 335; ZD958: Zhengdan 958.

    表2 葉片功能期模型參數(shù)表

    、、代表高斯函數(shù)參數(shù),表示曲線的漸近線,表示曲線的振幅,+表示最大功能期,表示最大功能期所對應葉位,表示曲線的陡度。

    ,,represent the Gaussian function parameter,is the asymptote of the curve,is the amplitude of the curve,+is the longest longevity,is the nodal position of the most longevity, andis the width of the curve.

    2.3 葉片功能期模型構(gòu)建的簡化方法

    葉片功能期模型構(gòu)建過程中需整株葉片功能期數(shù)據(jù), 試驗數(shù)據(jù)測試工作量大。通過對已構(gòu)建的葉片功能期模型分析發(fā)現(xiàn), 對高斯函數(shù)求導可以獲得曲線變化的轉(zhuǎn)折點(表3)。一階導為0的點為功能期最大值所在點; 二階導為0的點為功能期變化速率最大的點; 三階導為0的點為功能期開始快速增大的點, 根據(jù)總?cè)~片數(shù)找到取值范圍內(nèi)的3個轉(zhuǎn)折點取整獲得轉(zhuǎn)折關(guān)鍵葉位, 而葉1、為模型的起始點和終止點, 影響著曲線的變幅。從理論上, 通過5個轉(zhuǎn)折葉位功能期可以擬合構(gòu)建玉米全株葉片功能期模型。

    表3 高斯模型特征分析

    “”表示頂部葉片所在葉位?!啊?represents the position of the most top leaf.

    利用本試驗獲取的1、4、9、16和20(22) 5個葉位實測葉片功能期構(gòu)建葉片功能期模型, 擬合曲線的2值均達到0.96。并用轉(zhuǎn)折葉位之外的其他葉葉片功能期對模型進行驗證, 歸一化均方根誤差NRMSE均為7.8%, 表明利用5個轉(zhuǎn)折葉位葉片功能期能夠擬合整株不同葉位葉片功能期。

    圖5 2015–2017年5個葉位葉片功能期擬合圖

    虛線表示穗位葉所在, XY335表示先玉335, ZD958表示鄭單958。

    The dashed line represents the position of ear leaf. XY335: Xianyu 335; ZD958: Zhengdan 958.

    圖6 模型擬合值與實測值比較

    XY335表示先玉335, ZD958表示鄭單958。

    XY335: Xianyu 335; ZD958: Zhengdan 958.

    2.4 葉片功能期模型的應用

    圖7 葉片分組示意圖

    S1表示緩慢增長葉組, S2表示快速增長葉組, S3表示緩慢下降葉組。

    S1: a slow-growing leaf group; S2: a fast-growing leaf group; S3: a slow-declining leaf group.

    表4 葉片分組區(qū)間

    S1表示緩慢增長葉組, S2表示快速增長葉組, S3表示緩慢下降葉組。

    S1: a slow-growing leaf group; S2: a fast-growing leaf group; S3: a slow-declining leaf group.

    3 討論

    為了量化植株的生產(chǎn)能力, 前人對葉片生長發(fā)育過程做了大量的技術(shù)和方法研究。利用WOFOST模型[27-28]、DSSAT模型[29-30]、APSIM模型[31-32]等作物生長模型對葉片伸長過程、葉面積增長過程、葉面積指數(shù)動態(tài)變化及產(chǎn)量形成建立數(shù)學模型, 對不同空間尺度和環(huán)境條件下進行驗證和適用性研究。Lizaso等[3]提出用高斯函數(shù)描述葉片功能期的方法, 并聯(lián)系CRESE-MAIZE模型更準確地分析葉片生長發(fā)育過程, 但因其僅使用了1個品種、1個年份的數(shù)據(jù), 受分析數(shù)據(jù)量限制, 研究不夠深入; 2012年Kim等[24]借鑒Lizaso等[3]的方法, 將高斯函數(shù)應用于MAIZISM模型中, 對該方法進行了拓展。本文用2015—2017年3年2個品種的數(shù)據(jù)分析, 發(fā)現(xiàn)這種方法在年際間穩(wěn)定, 品種間有區(qū)別度, 并賦予了函數(shù)參數(shù)的生理學意義, 證實了高斯函數(shù)對玉米葉片功能期模擬的可行性, 為該方法的應用提供了科學的佐證, 為玉米生產(chǎn)能力的量化分析提供了思路。

    玉米不同節(jié)位葉片的生產(chǎn)能力不同, 對物質(zhì)生產(chǎn)及產(chǎn)量形成的貢獻存在主次之分[33]。前人[34-35]按葉片生長速度、供生長中心將玉米葉片分為根葉(1~6)、莖葉(7~11)、穗葉(12~16)和粒葉(17~頂葉) 4組, 也有人[6,36-37]按功能期天數(shù)長短將玉米葉片分為基部、下部、中部和上部葉4組, 與上述按功能劃分的區(qū)間基本一致。本研究以XY335和ZD958為試驗材料, 建立了不同葉位葉片功能期(GDD8)模型, 并通過模型特征參數(shù)從玉米整株葉片功能期變化動態(tài)切入將葉片分為緩慢增長(S1)、快速增長(S2)、緩慢下降(S3) 3組, 第4葉是S1與S2的分界點, 第16葉是S2與S3的分界點。三葉期是玉米的離乳期, 前期生長所需養(yǎng)分由籽粒胚乳提供[38-39], 所以1~3葉功能期較短且增速緩慢, 4葉以上功能期開始快速增長, 至16葉功能期達最大, 16葉以上營養(yǎng)物質(zhì)迅速向籽粒轉(zhuǎn)運, 功能期縮短。XY335的最大功能期所在葉位為穗上第2葉, 與范盼盼等[40]研究XY335穗上葉光合生產(chǎn)能力大于穗位葉的結(jié)果一致, ZD958的最大功能期所在葉位與穗位葉重合, 傳統(tǒng)上認為穗位葉功能期最長在品種間表現(xiàn)不一致。本研究提出的基于葉片功能期分組為解析葉片生長發(fā)育規(guī)律提供了參考。

    本研究是在長期氮肥定位試驗地的不施氮處理下進行, 玉米生育期縮短, 雖然能夠完整表現(xiàn)不同葉位葉片功能期, 仍與正常管理的植株生產(chǎn)存在一定差異, 但本研究的數(shù)據(jù)分析表明, 高斯函數(shù)準確地模擬了試驗條件下葉片功能期的變化特征, 證明該模型的科學性。因此可根據(jù)本研究結(jié)果, 拓展到不同生態(tài)條件、不同管理措施的葉片功能期模擬, 并利用高斯函數(shù)的特征參數(shù)進行玉米功能期不同影響因素的量化分析, 拓展了玉米生長研究的分析角度, 為玉米生產(chǎn)能力量化評估提供了借鑒。

    4 結(jié)論

    用2個品種3個年份的研究數(shù)據(jù)實證了利用高斯函數(shù)構(gòu)建玉米葉片功能期模型的科學性和可靠性,并利用高斯函數(shù)的特征參數(shù)提出了利用第1、第4、第9、第16和頂部葉片5個葉位構(gòu)建先玉335和鄭單958整株葉片功能期模型的方法, 簡化了模型的構(gòu)建過程。還提出了基于高斯函數(shù)構(gòu)建葉片功能期模型的葉片分組, 探索了葉片功能期模型的應用, 為葉片功能期的量化分析提供了理論和方法。

    [1] 陳國平. 玉米的干物質(zhì)生產(chǎn)與分配(綜述). 玉米科學, 1994, 2(1): 48–53. Chen G P. Dry matter production and distribution of maize (review)., 1994, 2(1): 48–53 (in Chinese).

    [2] 董樹亭, 高榮岐, 胡昌浩, 王群瑛, 王空軍. 玉米花粒期群體光合性能與高產(chǎn)潛力研究. 作物學報, 1997, 23: 318–325. Dong S T, Gao R Q, Hu C H, Wang Q Y, Wang K J. Study of canopy photosynthesis property and high yield potential after anthesis in maize., 1997, 23: 318–325 (in Chinese with English abstract).

    [3] Lizaso J I, Batchelor W D, Westgate M E. A leaf area model to simulate cultivar-specific expansion and senescence of maize leaves., 2003, 80: 1–17.

    [4] 羅瑤年, 張建華, 李霞. 玉米葉片的衰老. 玉米科學, 1992, 0(1): 40–43. Luo Y N, Zhang J H, Li X. Senescence of maize leaves., 1992, 0(1): 40–43 (in Chinese).

    [5] 王志剛. 超高產(chǎn)春玉米根冠結(jié)構(gòu)、功能特性及農(nóng)藝節(jié)水補償機制研究. 內(nèi)蒙古農(nóng)業(yè)大學博士學位論文, 內(nèi)蒙古呼和浩特, 2009. Wang Z G. Structural and Functional Properties of Canopy and Root of Super High Yield Spring Maize & Agronomic Water Saving Compensatory Mechanism. PhD Dissertation of Inner Mongolla Agricultural University, Hohhot, Inner Mongolla, China, 2009 (in Chinese with English abstract).

    [6] 張福喜, 蔡瑞國, 張敏, 李彥生, 韓金玲, 周印富, 王文頗. 春玉米葉片生長動態(tài)變化及其對種植密度的響應. 河北科技師范學院學報, 2009, 23(2): 15–21. Zhang F X, Cai R G, Zhang M, Li Y S, Han J L, Zhou Y F, Wang W P. Dynamic change of leaf growth in spring maize and its response to planting density., 2009, 23(2): 15–21 (in Chinese with English abstract).

    [7] 李文娟. 玉米(L.)氮素營養(yǎng)生理及其與葉片衰老的關(guān)系. 中國農(nóng)業(yè)科學院博士學位論文, 北京, 2012. Li W J. Nitrogen Physiology and Its Relation to Leaf Senescence of Maize (L.). PhD Dissertation of Chinese Academy of Agricultural Sciences, Beijing, China, 2012 (in Chinese with English abstract).

    [8] 張子山, 李耕, 高輝遠, 劉鵬, 楊程, 孟祥龍, 孟慶偉. 玉米持綠與早衰品種葉片衰老過程中光化學活性的變化. 作物學報, 2013, 39: 93–100. Zhang Z S, Li G, Gao H Y, Liu P, Yang C, Meng X L, Meng Q L. Changes of photochemistry activity during senescence of leaves in stay green and quick-leaf-senescence inbred lines of maize., 2013, 39: 93–100 (in Chinese with English abstract).

    [9] Wolfe D W, Henderson D W, Hsiao T C, Alvino A. Interactive water and nitrogen effects on senescence of maize: I. Leaf area duration, nitrogen distribution, and yield., 1962, 80: 859–864.

    [10] Trápani N, Hall A J. Effects of leaf position and nitrogen supply on the expansion of leaves of field grown sunflower (L.)., 1996, 184: 331–340.

    [11] 朱玉芹, 楊雙, 蔡鑫茹, 李萬良, 李景和. 玉米葉片早衰及子粒敗育. 吉林農(nóng)業(yè)科學, 2008, 33(3): 9–10. Zhu Y Q, Yang S, Cai X R, Li W L, Li J H. Studies on leaves early ageing and kernels abortion of maize.,2008, 33(3): 9–10 (in Chinese with English abstract).

    [12] 李軍虎, 張翠綿, 杜義英, 霍立勇, 滕慧穎, 趙會薇. 化控條件下密度對夏玉米產(chǎn)量及產(chǎn)量性狀的影響. 玉米科學, 2014, 22(3): 72–76. Li J H, Zhang C M, Du Y Y, Huo L Y, Teng H Y, Zhao H W. Effects of density on yield and yield traits of summer maize under chemical control., 2014, 22(3): 72–76 (in Chinese with English abstract).

    [13] 杜桂娟, 張哲, 楊姝, 馬鳳江, 劉洋. 淺談玉米早衰及其調(diào)控. 遼寧農(nóng)業(yè)科學, 2012, (2): 58–61.Du G J, Zhang Z, Yang S, Ma F J, Liu Y. Discussion on premature senescence of maize and its prevention technology., 2012, (2): 58–61 (in Chinese with English abstract).

    [14] 劉小芳. 玉米根系吸水調(diào)控機制. 中國科學院博士學位論文, 北京, 2013. Liu X F. Corn Root Water Absorption Control Mechanism. PhD Dissertation of Chinese Academy of Sciences, Beijing, China, 2013 (in Chinese with English abstract).

    [15] 劉全永, 李林峰, 李娜, 寇艷玲, 邵瑞鑫, 楊青華. 玉米早衰影響因素及其調(diào)控措施. 河南農(nóng)業(yè)科學, 2015, 44(5): 18–21. Liu Q Y, Li L F, Li N, Kou Y L, Shao R X, Yang Q H. Factors affecting premature senescence of maize and its controlling measures.,2015, 44(5): 18–21 (in Chinese with English abstract).

    [16] 熊炳霖, 王鑫月, 陳道鉗, 王仕穩(wěn), 殷俐娜, 鄧西平. 苗期玉米葉片碳氮平衡與干旱誘導的葉片衰老之關(guān)系. 西北植物學報, 2016, 36: 534–541. Xiong B L, Wang X Y, Chen D Q, Wang S W, Yin L N, Deng X P. Carbon/nitrogen balance associate with drought-induced leaf senescence in maize () seedling., 2016, 36: 534–541 (in Chinese with English abstract).

    [17] Pereira A R. Simulation of ecophysiological processes of growth in several annual crops., 1992, 57: 312–314.

    [18] 謝云, Kiniry J R. 國外作物生長模型發(fā)展綜述. 作物學報, 2002, 28: 190–195. Xie Y, Kiniry J R. A review on the development of crop modeling and its application., 2002, 28: 190–195 (in Chinese with English abstract).

    [19] 林忠輝, 莫興國, 項月琴. 作物生長模型研究綜述. 作物學報, 2003, 29: 750–758. Lin Z H, Mo X G, Xiang Y Q. Research advances on crop growth models., 2003, 29: 750–758 (in Chinese with English abstract).

    [20] 孫寧, 馮利平. 利用冬小麥作物生長模型對產(chǎn)量氣候風險的評估. 農(nóng)業(yè)工程學報, 2005, 21(2): 106–110. Sun N, Feng L P. Assessing the climatic risk to crop yield of winter wheat using crop growth models., 2005, 21(2): 106–110 (in Chinese with English abstract).

    [21] Wei X, Holman I, Conway D, Lin E, Li Y. A crop model cross calibration for use in regional climate impacts studies., 2008, 213: 365–380.

    [22] Birch C J, Hammer G L, Rickert K G. Improved methods for predicting individual leaf area and leaf senescence in maizeL.., 1998, 49: 249–262.

    [23] Birch C J, Vos J, Putten P E L. Plant development and leaf area production in contrasting cultivars of maize grown in a cool temperate environment in the field., 2003, 19: 173–188.

    [24] Kim S H, Yang Y, Timlin D J, Fleisher D H, Dathe A, Reddy V R, Staver K. Modeling temperature responses of leaf growth, deve-lopment, and biomass in maize with MAIZSIM., 2012, 104: 1523–1537.

    [25] 李少昆. 北方春玉米田間種植手冊. 北京: 中國農(nóng)業(yè)出版社, 2011. pp 112–113. Li S K. North Spring Corn Field Planting Manual. Beijing: China Agriculture Press, 2011. pp 112–113 (in Chinese).

    [26] Abendroth L J, Elmore R W, Boyer M J, Marlay S K. Corn Growth and Development. Iowa: Iowa State University Extension, 2011. pp 6–7.

    [27] Ittersum M K V, Leffelaar P A, Keulen H V, Kropff M J, Bastiaans L, Goudriaan J. On approaches and applications of the Wageningen crop models., 2003, 18: 201–234.

    [28] 鄔定榮, 歐陽竹, 趙小敏, 于強, 羅毅. 作物生長模型WOFOST在華北平原的適用性研究. 植物生態(tài)學報, 2003, 27: 594–602. Wu D R, Ou-Yang Z, Zhao X M, Yu Q, Luo Y. The applicability research of WOFOST model in north China plant., 2003, 27: 594–602 (in Chinese with English abstract).

    [29] Thorp K R, Hunsaker D J, French A N. Assimilating leaf area index estimates from remote sensing into the simulations of a cropping systems model., 2010, 53: 251–262.

    [30] Lizaso J I, Boote K J, Jones J W, Porter C H, Echarte L, Westgate M E, Sonohat G. CSM-IXIM: A new maize simulation model for DSSAT version 4.5., 2011, 103: 766–779.

    [31] Asseng S, Keating B A, Fillery I R P, Gregory P J, Bowden J W, Turner N C, Palta J A, Abrecht D G. Performance of the APSIM-wheat model in western Australia., 1998, 57: 163–179.

    [32] Asseng S, Keulen H V, Stol W. Performance and application of the APSIM N wheat model in the Netherlands., 2000, 12: 37–54.

    [33] 王忠孝. 山東玉米. 北京: 中國農(nóng)業(yè)出版社, 1999. pp 87–88. Wang Z X. Shandong Corn. Beijing: China Agriculture Press, 1999. pp 87–88 (in Chinese).

    [34] 任其云. 玉米不同生長階段“生長中心”和“供長中心葉”的探討. 植物學報, 1975, 17: 66–70. Ren Q Y. Studies on “Growth Center” and “Growth Center Leaf” at different stages of development in the corn plant., 1975, 17: 66–70 (in Chinese with English abstract).

    [35] 文宗群, 苑華, 陳慶惠. Sc704玉米葉片生長、根系發(fā)育及穗分化觀察. 新疆農(nóng)業(yè)科學, 1989, (2): 13–14. Wen Z Q, Fan H, Chen Q H. Sc704 observation of leaf growth, root development and ear differentiation of maize., 1989, (2): 13–14 (in Chinese).

    [36] 譚鳳梧. 玉米葉片生長發(fā)育規(guī)律及其功能的研究. 沈陽農(nóng)業(yè)大學學報, 1989, 20(1): 30–34. Tan F W. Studies on maize leaf’s growth and development rhythm and its function., 1989, 20(1): 30–34 (in Chinese).

    [37] 張建華, 李霞. 礦質(zhì)元素對玉米三葉期幼苗生長的影響. 中國土壤與肥料, 1992, (1): 32–35. Zhang J H, Li X. Effect of mineral elements on growth of maize seedlings at three-leaf stage., 1992, (1): 32–35 (in Chinese).

    [38] 郭慶法, 王慶成, 汪黎明. 中國玉米栽培學. 上海: 上??茖W技術(shù)出版社, 2004. pp 17–167. Guo Q F, Wang Q C, Wang L M. Chinese Maize Cultivation. Shanghai: Shanghai Scientific and Technical Publishers, 2004. pp 17–167 (in Chinese).

    [39] 王洋, 齊曉寧, 邵金鋒, 劉兆永, 李澤興. 遮陰對不同品種玉米葉片形態(tài)和功能的影響. 土壤與作物, 2009, 25: 450–453. Wang Y, Qi X N, Zhao J F, Liu Z Y, Li Z X. Effects of shading at the whole growth stage on the leaf morphology and function of different maize varieties., 2009, 25: 450–453 (in Chinese with English abstract).

    [40] 范盼盼, 謝瑞芝, 明博, 李姚姚, 王克如, 侯鵬, 李少昆. 基于不同葉位受光條件的玉米冠層光合生產(chǎn)能力分析. 玉米科學, 2017, 25(5): 68–72.Fan P P, Xie R X, Ming B, Li Y Y, Wang K R, Li S K. Analysis of photosynthetic productivity of maize canopy based on light conditions of different leaf positions., 2017, 25(5): 68–72 (in Chinese with English abstract).

    Establishment and application of spring maize leaf longevity model based on Gaussian function

    LI Yao-Yao, FAN Pan-Pan, MING Bo, WANG Chun-Xia, WANG Ke-Ru, HOU Peng, XIE Rui-Zhi*, and LI Shao-Kun*

    Institute of Crop Sciences, Chinese Academy of Agricultural Sciences / Key Laboratory of Crop Physiology and Ecology, Beijing 100081, China

    maize; leaf longevity; Gaussian function; growth model

    2018-12-11;

    2019-04-15;

    2019-04-22.

    10.3724/SP.J.1006.2019.83082

    謝瑞芝, E-mail: xieruizhi@caas.cn, Tel: 010-82105791; 李少昆, E-mail: lishaokun@caas.cn, Tel: 010-82108891

    E-mail: 981457781@qq.com, Tel: 010-82105791

    本研究由國家重點研發(fā)計劃項目(2017YFD0300302), 國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設專項(CARS-02-20)和中國農(nóng)業(yè)科學院科技創(chuàng)新工程項目資助。

    This study was supported by the National Key Research and Development Program of China (2017YFD0300302), the China Agriculture Research System (CARS-02-20), and the Science and Technology Innovation Program of Chinese Academy of Agricultural Sciences.

    URL: http://kns.cnki.net/kcms/detail/11.1809.S.20190419.1016.006.html

    猜你喜歡
    葉位高斯葉片
    小高斯的大發(fā)現(xiàn)
    云煙87不同葉位煙葉烤后單葉重及等級質(zhì)量分析
    月季葉片“雕刻師”
    干旱脅迫下辣椒幼苗光合與熒光參數(shù)測定的最佳葉片
    不同葉位桑葉茶的營養(yǎng)活性成分分析
    兩種噴涂方法在TRT葉片上的對比應用及研究
    天才數(shù)學家——高斯
    云煙85不同葉位和不同部位烤后煙葉碳氮含量和碳氮比的差異性研究
    基于CFD/CSD耦合的葉輪機葉片失速顫振計算
    航空學報(2015年4期)2015-05-07 06:43:30
    有限域上高斯正規(guī)基的一個注記
    日本成人三级电影网站| 夜夜爽天天搞| 国产精品一区二区三区四区免费观看| 能在线免费观看的黄片| 亚洲无线在线观看| 少妇丰满av| 美女脱内裤让男人舔精品视频 | 尤物成人国产欧美一区二区三区| 岛国在线免费视频观看| 国产黄a三级三级三级人| 欧美日韩在线观看h| 大香蕉久久网| 成人欧美大片| 国产精品一二三区在线看| 色播亚洲综合网| 国产成人精品婷婷| 国产av在哪里看| 国产精品一区二区三区四区免费观看| 观看免费一级毛片| 在线免费观看不下载黄p国产| 天美传媒精品一区二区| 国产私拍福利视频在线观看| 日本-黄色视频高清免费观看| 久久久久久久亚洲中文字幕| 女人被狂操c到高潮| 97热精品久久久久久| 国产黄a三级三级三级人| www日本黄色视频网| 悠悠久久av| 特级一级黄色大片| 免费看光身美女| 中文字幕制服av| 搡老妇女老女人老熟妇| 久久亚洲精品不卡| 精品一区二区三区人妻视频| 熟妇人妻久久中文字幕3abv| 日韩大尺度精品在线看网址| 一级黄片播放器| 国产精品久久视频播放| 欧美性感艳星| 哪里可以看免费的av片| 日本色播在线视频| 国产午夜福利久久久久久| 人人妻人人看人人澡| 久久精品国产鲁丝片午夜精品| 久久九九热精品免费| 校园人妻丝袜中文字幕| 日韩欧美三级三区| 成人三级黄色视频| 两个人的视频大全免费| 免费看光身美女| 可以在线观看毛片的网站| 日本撒尿小便嘘嘘汇集6| 免费人成视频x8x8入口观看| 国产精品一区二区三区四区久久| 中文欧美无线码| 亚洲人成网站在线播放欧美日韩| 成人永久免费在线观看视频| 久久久久久久亚洲中文字幕| 如何舔出高潮| 可以在线观看的亚洲视频| 亚洲高清免费不卡视频| 国产精品精品国产色婷婷| 国产精品,欧美在线| 亚洲av中文av极速乱| 国产精品.久久久| 人人妻人人看人人澡| 亚洲第一电影网av| 99热6这里只有精品| 午夜福利视频1000在线观看| 久久6这里有精品| 岛国毛片在线播放| 又爽又黄a免费视频| 99热6这里只有精品| 日韩亚洲欧美综合| 久久午夜福利片| 搡老妇女老女人老熟妇| 成人特级av手机在线观看| 99久久精品热视频| 欧美日韩在线观看h| 欧美又色又爽又黄视频| .国产精品久久| 不卡视频在线观看欧美| 精品国内亚洲2022精品成人| 欧美成人免费av一区二区三区| 亚洲欧美中文字幕日韩二区| 69人妻影院| 91久久精品国产一区二区三区| 国产精品一及| 亚洲自拍偷在线| 国产不卡一卡二| 99久久中文字幕三级久久日本| 久久久久久九九精品二区国产| 午夜福利高清视频| 18禁黄网站禁片免费观看直播| 午夜福利在线在线| 黄片无遮挡物在线观看| 91aial.com中文字幕在线观看| 亚洲自拍偷在线| 日日啪夜夜撸| 三级国产精品欧美在线观看| 日韩成人伦理影院| 国产精品电影一区二区三区| 成人二区视频| 亚洲精品456在线播放app| 人妻系列 视频| 国产成人影院久久av| 给我免费播放毛片高清在线观看| 亚洲五月天丁香| 欧美性感艳星| 国产精品免费一区二区三区在线| 久久国产乱子免费精品| av.在线天堂| 99视频精品全部免费 在线| 白带黄色成豆腐渣| 国产伦一二天堂av在线观看| 欧美最黄视频在线播放免费| 免费观看精品视频网站| 久久久久久久久中文| 小蜜桃在线观看免费完整版高清| 99久久成人亚洲精品观看| 99久久人妻综合| 啦啦啦韩国在线观看视频| 亚洲天堂国产精品一区在线| 亚洲图色成人| 日韩 亚洲 欧美在线| 国产亚洲av嫩草精品影院| 亚洲aⅴ乱码一区二区在线播放| 99在线人妻在线中文字幕| 男人舔女人下体高潮全视频| 亚洲一区二区三区色噜噜| 色哟哟哟哟哟哟| 亚洲精品自拍成人| 精品人妻偷拍中文字幕| 一级黄色大片毛片| 成人永久免费在线观看视频| 欧美日韩国产亚洲二区| 国产色婷婷99| 人妻久久中文字幕网| 欧美3d第一页| 99久国产av精品| 一级av片app| 狂野欧美白嫩少妇大欣赏| 男人和女人高潮做爰伦理| 日韩亚洲欧美综合| 国产精品久久久久久精品电影小说 | 此物有八面人人有两片| 国产精品一区二区性色av| 国产亚洲精品久久久久久毛片| 亚洲在久久综合| 91狼人影院| 老熟妇乱子伦视频在线观看| 夫妻性生交免费视频一级片| 亚洲av二区三区四区| 亚洲国产高清在线一区二区三| 人妻制服诱惑在线中文字幕| 91狼人影院| 免费看光身美女| 在线天堂最新版资源| 国产熟女欧美一区二区| av天堂中文字幕网| 亚洲电影在线观看av| 三级国产精品欧美在线观看| 男的添女的下面高潮视频| 乱码一卡2卡4卡精品| 亚洲一区高清亚洲精品| 国产v大片淫在线免费观看| 中文字幕免费在线视频6| 别揉我奶头 嗯啊视频| 赤兔流量卡办理| 免费看光身美女| 两性午夜刺激爽爽歪歪视频在线观看| 日本熟妇午夜| 国内精品一区二区在线观看| 变态另类成人亚洲欧美熟女| 网址你懂的国产日韩在线| 婷婷亚洲欧美| 天堂av国产一区二区熟女人妻| 国产一级毛片在线| 免费人成视频x8x8入口观看| 色播亚洲综合网| 国产精华一区二区三区| 精品无人区乱码1区二区| 成人二区视频| 国产大屁股一区二区在线视频| 老女人水多毛片| 国产色婷婷99| 成人美女网站在线观看视频| 免费一级毛片在线播放高清视频| 国产精品不卡视频一区二区| 欧美+日韩+精品| 九草在线视频观看| 亚洲成人av在线免费| 欧美日韩综合久久久久久| av.在线天堂| 国内少妇人妻偷人精品xxx网站| 色哟哟哟哟哟哟| 国产一区二区三区av在线 | 毛片一级片免费看久久久久| 99久久精品国产国产毛片| 午夜福利视频1000在线观看| 精品一区二区三区视频在线| 国产视频内射| 午夜福利成人在线免费观看| 观看免费一级毛片| 麻豆成人av视频| 黄色日韩在线| 久久热精品热| 久久久久久久久大av| a级毛片a级免费在线| 男女视频在线观看网站免费| 能在线免费观看的黄片| 日韩大尺度精品在线看网址| 成人毛片60女人毛片免费| 亚州av有码| 久久韩国三级中文字幕| 亚洲人成网站在线播| 国产亚洲精品av在线| 中文亚洲av片在线观看爽| 亚洲精品乱码久久久久久按摩| 91麻豆精品激情在线观看国产| 日韩强制内射视频| 搡老妇女老女人老熟妇| 亚洲国产精品成人久久小说 | 小蜜桃在线观看免费完整版高清| 三级经典国产精品| 欧美一区二区精品小视频在线| 成人三级黄色视频| av又黄又爽大尺度在线免费看 | av.在线天堂| 久久99蜜桃精品久久| 美女国产视频在线观看| 深夜精品福利| 嫩草影院精品99| 国产成人a∨麻豆精品| 桃色一区二区三区在线观看| 午夜久久久久精精品| 亚洲欧美日韩高清在线视频| kizo精华| 成年版毛片免费区| 两个人的视频大全免费| 久久精品国产自在天天线| 久久久久久久午夜电影| av在线播放精品| 男人的好看免费观看在线视频| 亚洲成人久久爱视频| 桃色一区二区三区在线观看| 老女人水多毛片| 老师上课跳d突然被开到最大视频| 日本三级黄在线观看| 免费av观看视频| 少妇高潮的动态图| 成年版毛片免费区| 日韩一本色道免费dvd| 日韩欧美一区二区三区在线观看| 91久久精品国产一区二区三区| 狠狠狠狠99中文字幕| 国内精品宾馆在线| 男人狂女人下面高潮的视频| 99国产极品粉嫩在线观看| 国产亚洲精品久久久久久毛片| 亚洲乱码一区二区免费版| 美女xxoo啪啪120秒动态图| 亚洲国产精品成人久久小说 | 欧美不卡视频在线免费观看| 国产女主播在线喷水免费视频网站 | 国产真实伦视频高清在线观看| 国产成年人精品一区二区| 搡老妇女老女人老熟妇| 亚洲av免费高清在线观看| 嘟嘟电影网在线观看| 可以在线观看毛片的网站| 成人高潮视频无遮挡免费网站| 两个人的视频大全免费| 男女做爰动态图高潮gif福利片| 97在线视频观看| 中国美女看黄片| 久久99热6这里只有精品| 国产亚洲av片在线观看秒播厂 | 国产一级毛片七仙女欲春2| 18禁黄网站禁片免费观看直播| 亚洲自偷自拍三级| 日日干狠狠操夜夜爽| 简卡轻食公司| 99久久成人亚洲精品观看| 欧美日本视频| 日韩精品青青久久久久久| 国产三级在线视频| 国产精品免费一区二区三区在线| 免费观看人在逋| 国产精品久久久久久精品电影小说 | 成人综合一区亚洲| 三级经典国产精品| 亚洲欧美日韩无卡精品| 国产精品国产三级国产av玫瑰| 午夜亚洲福利在线播放| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久丰满| 国产免费男女视频| 91午夜精品亚洲一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 97热精品久久久久久| 婷婷色av中文字幕| 国产av麻豆久久久久久久| 亚洲人成网站在线播放欧美日韩| 亚洲欧美精品专区久久| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 一卡2卡三卡四卡精品乱码亚洲| 国产精品人妻久久久久久| 国产精品久久久久久久电影| 97人妻精品一区二区三区麻豆| 国产亚洲91精品色在线| 欧美成人a在线观看| 女的被弄到高潮叫床怎么办| 国产探花极品一区二区| 神马国产精品三级电影在线观看| 亚洲精品色激情综合| 日本欧美国产在线视频| 中文字幕久久专区| 亚洲一区二区三区色噜噜| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲成人av在线免费| 欧美bdsm另类| 日本爱情动作片www.在线观看| 非洲黑人性xxxx精品又粗又长| 联通29元200g的流量卡| 成年女人看的毛片在线观看| 99久久精品热视频| 欧美在线一区亚洲| 蜜桃亚洲精品一区二区三区| 中文资源天堂在线| 日本免费一区二区三区高清不卡| 色吧在线观看| 一个人观看的视频www高清免费观看| 少妇人妻精品综合一区二区 | 欧美不卡视频在线免费观看| 91aial.com中文字幕在线观看| 国产一区二区三区在线臀色熟女| 最近2019中文字幕mv第一页| 国产成人精品一,二区 | 丝袜美腿在线中文| 国产视频内射| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 国产成人91sexporn| 一个人观看的视频www高清免费观看| 欧美3d第一页| 亚洲性久久影院| 2021天堂中文幕一二区在线观| 一级毛片电影观看 | www.色视频.com| 91久久精品国产一区二区成人| 成年免费大片在线观看| 国产亚洲精品av在线| 最好的美女福利视频网| 99九九线精品视频在线观看视频| 国产黄色视频一区二区在线观看 | 99热这里只有是精品在线观看| 日韩大尺度精品在线看网址| 国产极品精品免费视频能看的| 搡女人真爽免费视频火全软件| 丝袜美腿在线中文| 成年av动漫网址| 99久久九九国产精品国产免费| 黄片wwwwww| 日本av手机在线免费观看| 久久人人爽人人片av| 日韩在线高清观看一区二区三区| 91午夜精品亚洲一区二区三区| 国产成人91sexporn| a级毛片免费高清观看在线播放| 人妻少妇偷人精品九色| 国内精品久久久久精免费| 亚洲成人久久性| 国产国拍精品亚洲av在线观看| 日韩亚洲欧美综合| 99热精品在线国产| 神马国产精品三级电影在线观看| 男女边吃奶边做爰视频| 午夜激情福利司机影院| 免费一级毛片在线播放高清视频| 久久欧美精品欧美久久欧美| 中文字幕熟女人妻在线| a级毛色黄片| 国产成人一区二区在线| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 一级黄片播放器| 久久久久久久久久成人| 3wmmmm亚洲av在线观看| 美女内射精品一级片tv| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 男女下面进入的视频免费午夜| 国产精品久久视频播放| 中文在线观看免费www的网站| 三级男女做爰猛烈吃奶摸视频| 久久久久久伊人网av| 最后的刺客免费高清国语| 乱码一卡2卡4卡精品| 亚洲电影在线观看av| 波多野结衣高清作品| 亚洲国产色片| 欧美高清成人免费视频www| av在线观看视频网站免费| 一个人免费在线观看电影| 人妻久久中文字幕网| 夫妻性生交免费视频一级片| 国产精品免费一区二区三区在线| 久久九九热精品免费| 免费一级毛片在线播放高清视频| 日产精品乱码卡一卡2卡三| 嫩草影院新地址| 国产精品永久免费网站| 国产精品久久久久久久电影| 有码 亚洲区| 色视频www国产| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 少妇人妻精品综合一区二区 | 一区二区三区高清视频在线| 日韩强制内射视频| 久久精品国产99精品国产亚洲性色| 你懂的网址亚洲精品在线观看 | 日本免费a在线| 午夜福利成人在线免费观看| 在线天堂最新版资源| 精品久久国产蜜桃| 欧美3d第一页| 三级男女做爰猛烈吃奶摸视频| 长腿黑丝高跟| 久久精品久久久久久久性| av福利片在线观看| 欧美性猛交╳xxx乱大交人| 国内精品久久久久精免费| 26uuu在线亚洲综合色| 男插女下体视频免费在线播放| 人人妻人人看人人澡| 日韩三级伦理在线观看| 青春草亚洲视频在线观看| 成年av动漫网址| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 直男gayav资源| 国产精品综合久久久久久久免费| 免费观看人在逋| 中文字幕人妻熟人妻熟丝袜美| 日韩av不卡免费在线播放| 欧美zozozo另类| 黄色欧美视频在线观看| 精品人妻视频免费看| 97热精品久久久久久| 成人午夜精彩视频在线观看| 一区福利在线观看| 成人av在线播放网站| 久久鲁丝午夜福利片| 国产精品一区二区三区四区久久| 麻豆av噜噜一区二区三区| 国产老妇女一区| 深爱激情五月婷婷| 级片在线观看| 国产极品天堂在线| 精品久久久久久久末码| 国产成人午夜福利电影在线观看| 99久久精品热视频| 青春草国产在线视频 | 自拍偷自拍亚洲精品老妇| 一级毛片我不卡| 亚洲,欧美,日韩| 亚洲av成人av| 精品久久久久久成人av| 一区二区三区高清视频在线| 69av精品久久久久久| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 亚洲人与动物交配视频| 国产白丝娇喘喷水9色精品| 成年女人永久免费观看视频| 亚洲精品乱码久久久v下载方式| 精品人妻偷拍中文字幕| 色播亚洲综合网| 精品久久久久久久人妻蜜臀av| 九九久久精品国产亚洲av麻豆| 免费在线观看成人毛片| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久亚洲中文字幕| 村上凉子中文字幕在线| 在现免费观看毛片| 久久久久久久久久久免费av| 97在线视频观看| 99久久成人亚洲精品观看| 亚洲成人精品中文字幕电影| 国产淫片久久久久久久久| 国产黄片美女视频| 我要看日韩黄色一级片| 欧美日韩精品成人综合77777| 好男人在线观看高清免费视频| 不卡视频在线观看欧美| 国产精品三级大全| 亚洲国产精品国产精品| 国产三级在线视频| 亚洲成人久久爱视频| 在线观看66精品国产| 黄色视频,在线免费观看| 久久久久久久亚洲中文字幕| 26uuu在线亚洲综合色| 日本撒尿小便嘘嘘汇集6| 寂寞人妻少妇视频99o| 亚洲五月天丁香| 欧美+日韩+精品| 精品久久国产蜜桃| 高清在线视频一区二区三区 | a级一级毛片免费在线观看| 久久精品国产清高在天天线| 女人被狂操c到高潮| 日韩一本色道免费dvd| 日本一二三区视频观看| 日韩在线高清观看一区二区三区| 免费看a级黄色片| 一本久久精品| kizo精华| 婷婷色综合大香蕉| 成人av在线播放网站| 在线a可以看的网站| 午夜亚洲福利在线播放| 国产成人午夜福利电影在线观看| 国产成人精品婷婷| 十八禁国产超污无遮挡网站| av免费在线看不卡| 中文字幕久久专区| 色综合站精品国产| 国产在线男女| 少妇的逼好多水| 最近最新中文字幕大全电影3| 欧美bdsm另类| 女人被狂操c到高潮| 亚洲精品日韩在线中文字幕 | 国产精品一区二区三区四区免费观看| 舔av片在线| 麻豆一二三区av精品| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲av天美| 国产真实伦视频高清在线观看| 亚洲四区av| 插逼视频在线观看| 少妇猛男粗大的猛烈进出视频 | 尾随美女入室| 亚洲欧洲日产国产| 99久久精品一区二区三区| 长腿黑丝高跟| 天堂网av新在线| 免费不卡的大黄色大毛片视频在线观看 | 久久久久免费精品人妻一区二区| 99在线人妻在线中文字幕| 99热这里只有是精品50| 国产精品一及| 亚洲图色成人| 国内揄拍国产精品人妻在线| 极品教师在线视频| 亚洲欧美精品专区久久| 韩国av在线不卡| 久久欧美精品欧美久久欧美| 简卡轻食公司| 三级国产精品欧美在线观看| 久久久久九九精品影院| 一个人看的www免费观看视频| 两个人视频免费观看高清| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| 国产午夜精品一二区理论片| 有码 亚洲区| 成人一区二区视频在线观看| 国产高清视频在线观看网站| 久久人妻av系列| 久久中文看片网| 久久精品国产鲁丝片午夜精品| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕制服av| 婷婷色综合大香蕉| 成熟少妇高潮喷水视频| 亚洲国产精品sss在线观看| 麻豆乱淫一区二区| 久久欧美精品欧美久久欧美| 最近的中文字幕免费完整| 亚洲真实伦在线观看| 嘟嘟电影网在线观看| 日韩精品有码人妻一区| 日本三级黄在线观看| www.色视频.com| 天堂影院成人在线观看| 国产成人91sexporn| 又爽又黄a免费视频| av免费在线看不卡| 久久久a久久爽久久v久久| 久久中文看片网| 久久人妻av系列| 深夜精品福利| 日韩视频在线欧美| 色哟哟·www| 日韩成人伦理影院| 91在线精品国自产拍蜜月| 国产成人freesex在线| 国产高清有码在线观看视频| 国产精品国产高清国产av| 男女视频在线观看网站免费| 人人妻人人澡人人爽人人夜夜 | 中出人妻视频一区二区| 成人三级黄色视频| 国产精品久久久久久久久免| eeuss影院久久| 国产久久久一区二区三区| www.色视频.com| 99热这里只有是精品50| 国产av不卡久久|