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

    AMMI模型和GGE雙標(biāo)圖對寧夏不同氣候類型區(qū)苦蕎品種穩(wěn)產(chǎn)性適應(yīng)性分析

    2024-04-12 00:00:00穆蘭海??饲?/span>杜燕萍楊崇慶陳一鑫

    摘要:為準(zhǔn)確評價參試苦蕎新品種的豐產(chǎn)性、穩(wěn)產(chǎn)性和生態(tài)適應(yīng)性,利用AMMI模型和GGE雙標(biāo)圖對5個參試品種在寧夏南部山區(qū)的半干旱區(qū)、半干旱易旱區(qū)和干旱區(qū)的6縣(區(qū))多年多點試驗數(shù)據(jù)進(jìn)行了品種的豐產(chǎn)性、穩(wěn)產(chǎn)性、適應(yīng)性和試點代表性分析及綜合評價。通過主成分PCA與品種的穩(wěn)定性(Dg)和試點的穩(wěn)產(chǎn)性(De)分析,篩選出適宜在寧夏南部山區(qū)種植的豐產(chǎn)性穩(wěn)定性和適應(yīng)性好的品種為晉苦蕎2號和川蕎3號,晉苦蕎2號產(chǎn)量為1 671~2 593 kg·hm-2,比對照平均增產(chǎn)17.3%。川蕎3號的產(chǎn)量為1 474~2 513 kg·hm-2,與對照相比平均增產(chǎn)7.5%;試點辨別力代表性以半干旱易旱區(qū)的海原縣,中部干旱帶的同心縣及半干旱區(qū)的西吉縣比較理想;試驗經(jīng)歷了干旱年份、正常年份和豐水年份,參試品種在干旱年份較豐水年份產(chǎn)量減產(chǎn)28.0%~45.3%。

    關(guān)鍵詞:苦蕎品種;AMMI模型;GGE雙標(biāo)圖;豐產(chǎn)性;穩(wěn)定性;生態(tài)適應(yīng)性

    寧南山區(qū)8縣(區(qū))氣候類型區(qū)由南向北為陰濕區(qū)、半干旱區(qū)、半干旱易旱區(qū)和中部干旱帶,適宜苦蕎種植,是我國北方重要苦蕎種植區(qū),連同甜蕎每年種植面積在5.85萬hm2左右,面積和產(chǎn)量分別占全國的15.1%和6.7%[1]。篩選高產(chǎn)、穩(wěn)產(chǎn)和適應(yīng)性廣的苦蕎品種并在生產(chǎn)中推廣應(yīng)用,不斷推動苦蕎特色產(chǎn)業(yè)快速發(fā)展是育種者追求的目標(biāo)。要對選育或引進(jìn)苦蕎品種(系)進(jìn)行多年多點試驗資料統(tǒng)計分析和綜合評價豐產(chǎn)穩(wěn)產(chǎn)性及生態(tài)適應(yīng)性,并進(jìn)行品種間差異顯著性比較,由于品種(品系)基因型與試點環(huán)境間普遍存在互作[2],因而對參試品種的豐產(chǎn)穩(wěn)定性和適應(yīng)性評價具有重要意義。

    近年來,國內(nèi)外應(yīng)用 AMMI(Additive Main Effects and Multiplicative Interaction)和GGE(Genotypemain Effects and Genotype×Environ-mentinteraction)等模型來分析多點試驗,AMMI 模型能較客觀地分析基因型(G)與環(huán)境(E)的互作效應(yīng)(G×E)[1-3],而GGE雙標(biāo)圖能夠同時考慮基因型和環(huán)境互作對品種、試點和劃分生態(tài)區(qū)域進(jìn)行有效評價,是利用多點試驗數(shù)據(jù)進(jìn)行品種和試驗點評價的理想方法[4-5]。GGE雙標(biāo)圖法已在燕麥、花生、油菜、大豆和馬鈴薯等作物品種試驗中廣泛應(yīng)用[6-8]。但采用AMMI模型和GGE雙標(biāo)圖進(jìn)行蕎麥品種豐產(chǎn)性、穩(wěn)定性和試點代表性辨別力綜合性評價的文獻(xiàn)相對較少。

    本研究擬利用R語言支持下GGE-Biplo軟件進(jìn)行雙標(biāo)圖分析與AMMI模型評估5個苦蕎品種在6個試點的豐產(chǎn)穩(wěn)定性、適應(yīng)性和試點代表性,為科學(xué)布設(shè)試驗點和合理優(yōu)化品種種植區(qū)域,篩選適合各區(qū)域大面積推廣應(yīng)用的苦蕎品種提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗地概況

    2017-2019年在寧南山區(qū)6縣(區(qū))布設(shè)試驗點。其中原州區(qū)頭營、西吉縣馬建和彭陽縣城陽代表半干旱區(qū),海原縣樹臺代表半干旱區(qū)易旱區(qū),同心縣預(yù)旺和鹽池縣花馬池代表中部干旱帶,每個試點參試品種和田間設(shè)計相同。

    由表1和圖1可以看出,寧南山區(qū)半干旱區(qū)(原州區(qū)、西吉縣和彭陽縣)、半干旱區(qū)易旱區(qū)(海原縣)和中部干旱帶(同心縣、鹽池縣)多年年均降水量分別為416.6~450.0 mm、350.0 mm和215.1~255.6 mm,多年平均氣溫分別為6.2~7.2 ℃、7.1 ℃和8.6 ℃;苦蕎品種生育期降水量為280.9~391.9 mm、345.1 mm和210.6~251.2 mm,生育期平均氣溫分別為16.7~19.4 ℃、18.7 ℃和21.3~22.0 ℃。

    1.2 材料

    參試品種名稱及選育單位見表2。2021-2022年完成品種生產(chǎn)示范和品種審定。

    1.3 方法

    1.3.1 試驗設(shè)計

    試驗采用隨機(jī)區(qū)組設(shè)計,6個試驗點統(tǒng)一田間設(shè)計方案,小區(qū)面積 10 m2(長5 m×寬2 m),每個小區(qū)種植7行區(qū),平均行距30 cm,留苗密度控制在每小區(qū)1 050株左右,3次重復(fù),隨機(jī)排列。播種前整地結(jié)合旋耕統(tǒng)一基肥磷酸二銨肥料150 kg·hm-2。5月下旬播種,6月中旬定苗并調(diào)查基本苗,其他田間管理同大田。每小區(qū)收獲全部植株進(jìn)行脫粒,各試驗點品種考種樣取中間1 m行長植株,其中選代表性15株進(jìn)行考種,并計算樣段內(nèi)平均株高,主莖分枝數(shù)、節(jié)間數(shù)、單株粒重、粒數(shù)及千粒重等。

    1.3.2 方差分析和AMMI模型分析

    本研究采用DPS 9.5和Excel 2010軟件進(jìn)行數(shù)據(jù)分析。利用AMMI模型分析其品種及試驗點間的穩(wěn)定性和豐產(chǎn)性參數(shù)及綜合評價[9-12]。利用AMMI模型和R語言軟件包GGE作圖進(jìn)行產(chǎn)量數(shù)據(jù)方差分析和主成分特征值分析[13-14],計算品種的豐產(chǎn)性穩(wěn)定性參數(shù)(Dg)和試點穩(wěn)定性參數(shù)(De)。AMMI模型的主要特點是將方差分析和主成分分析有機(jī)結(jié)合在一起,公式如下:

    式中,yge表示在某試點e中某一基因型g的平均產(chǎn)量;μ為總體平均值;αg為基因型平均偏差(基因型平均值減去總的平均值);βe為環(huán)境的平均偏差(各個環(huán)境的平均值減去總的平均值);λn為第n個主成分特征值的平方根;γgn為第n個主成分的環(huán)境得分;δgn為第n個主成分的基因型主成分得分;θge為殘差。3次重復(fù),則誤差項為εger,它等于γgn平均值與r個重復(fù)的單個觀察值之間的偏差,并滿足可加性。

    1.3.3 品種穩(wěn)定性與適應(yīng)性分析

    采用GGE 運用R語言支持下的GGE-Biplo軟件包功能進(jìn)行雙標(biāo)圖分析[11-13],分別評價各品種的豐產(chǎn)性和穩(wěn)產(chǎn)性、各試點的區(qū)分力與代表性、品種生態(tài)適應(yīng)性等。

    穩(wěn)定性參數(shù)是品種或試驗點的交互效應(yīng)主成分值(Interaction Principal Component Axis,IPCA)在多維空間中圖標(biāo)與原點的歐氏距離[14],公式如下:

    式中,n表示達(dá)到顯著水平的IPCA個數(shù);Dg(e)用來度量品種基因型或在環(huán)境相對穩(wěn)定性條件下的參數(shù),即為基因型或環(huán)境在n個IPCA上的歐式距離。

    2 結(jié)果與分析

    2.1 苦蕎品種產(chǎn)量的方差和AMMI模型分析

    寧南山區(qū)苦蕎多年多點品種區(qū)域試驗中,苦蕎品種產(chǎn)量受制于基因(Genotype, G)、環(huán)境(Environment,E)以及基因與環(huán)境互作效應(yīng)(G×E)的影響[15-17]。由表3可知,環(huán)境及基因和交互作用變異占總平方和的比例從小到大排序為:環(huán)境(56.3%)>基因(24.0%)>基因與環(huán)境互作效應(yīng)(19.8%)。相同的苦蕎品種在不同氣候類型區(qū)種植,其產(chǎn)量表現(xiàn)差異顯著,其中處理和環(huán)境(地點)達(dá)到極顯著。由此也可以看出,品種對生態(tài)環(huán)境的適應(yīng)性至關(guān)重要。

    AMMI模型主成分分析中,基因型和環(huán)境互作效應(yīng)可以進(jìn)一步分解達(dá)到顯著水平的3個互作效應(yīng)主成分軸PCA1 (Interaction Principle Component Axis)、PCA2和PCA3的變異平方和分別占互作量的58.9%、25.4%和15.7%,且F值均達(dá)到極顯著。說明此AMMI模型可以較好地解釋品種和試點互作效應(yīng)。

    2.2 參試品種與試點GGE雙標(biāo)圖分析

    2.2.1 品種的豐產(chǎn)性穩(wěn)產(chǎn)性

    由表4和圖2可知,參試品種在所有試驗點的平均產(chǎn)量依次為KV1(晉苦蕎2號)2 188 kg·hm-2>KV2(川蕎3號)2 004 kg·hm-2>KV4(黔黑蕎5號)1 956 kg·hm-2>KV3(云蕎2號)1 944 kg·hm-2>KV5(固原當(dāng)?shù)乜嗍w)1 865 kg·hm-2;根據(jù)參試品種的互作效應(yīng)達(dá)到顯著的主成分PCA1、PCA2和PCA3可分別計算出各參試品種的穩(wěn)產(chǎn)性參數(shù)Dg值,其Dg值越小,表示品種的穩(wěn)定性越好。其品種的穩(wěn)定性參數(shù)Dg值依次為KV1(晉苦蕎2號)7.80lt;KV3(云蕎2號)11.80lt;KV2(川蕎3號)12.97lt;KV5(固原當(dāng)?shù)乜嗍w)14.76lt;KV4(黔黑蕎5號)17.57,苦蕎品種的豐產(chǎn)性穩(wěn)定性與品種在多年多點平均產(chǎn)量排名基本一致。因此,依照不同苦蕎麥品種在不同氣候資源類型和生產(chǎn)條件下對產(chǎn)量的影響程度看,新品種(系)選育應(yīng)重點考慮既有豐產(chǎn)性又具備穩(wěn)產(chǎn)性作為理想目標(biāo)。

    2.2.2 試點代表性與穩(wěn)定性

    參試地點互作主成分值(PCA)與試點辨別力的穩(wěn)定性(De),其值越大,表明試點的穩(wěn)定性較好。由表5可知,苦蕎品種的區(qū)試點穩(wěn)定性De值依次為E4(海原縣)18.32>E5(同心縣)14.65>E2(西吉縣)11.96>E3(彭陽縣)9.20>E6(鹽池縣)8.54>E1(原州區(qū))6.64。說明苦蕎品種在代表試點E2(或E3)、E4和E5的穩(wěn)定性和代表性比較好。各參試地點苦蕎的產(chǎn)量依次表現(xiàn)為E1gt;E3gt;E2gt;E4gt;E5gt;E6。

    另外,苦蕎品種的生產(chǎn)能力在半干旱區(qū)試點E1、E2和E3平均產(chǎn)量為2 139 kg·hm-2;半干旱易旱區(qū)的試點E4平均產(chǎn)量為1 987 kg·hm-2;中部干旱帶E5和E6平均產(chǎn)量為1 772 kg·hm-2左右,較半干旱區(qū)和半干旱易旱區(qū)分別減產(chǎn)17.2%和10.8%。

    2.2.3 試點代表性與辨別力

    利用R語言的GGE-Biplot軟件包對參試品種在6個試點進(jìn)行代表性和辨別力GGE雙坐標(biāo)圖分析(圖3)。圖中帶箭頭的直線是平均環(huán)境軸,試點線段和平均環(huán)境之間的夾角表示該試點的代表性[18-21],角度越小,代表性越強。如果一個試驗點和平均環(huán)境軸夾角為鈍角,則表示該試點不適合作為試驗點。各試點線段的長短表示該試點的辨別力。

    由圖3可知,通過試點代表性和辨別力GGE雙坐標(biāo)圖分析,試點E1、E2和E3環(huán)境軸線段均同處于半干旱區(qū),為了減少設(shè)置試點費用成本,半干旱區(qū)的苦蕎試點以選擇原州區(qū)或西吉縣比較好,半干旱易旱區(qū)的試點選擇E4(海原縣),而中部干旱帶的E5(同心縣)與E6(鹽池縣)雖然處于相同氣候環(huán)境,但試點E6卻與平均環(huán)境軸的夾角最小,而且試點E6亦是寧夏蕎麥雜糧作物規(guī)模化種植和產(chǎn)業(yè)化經(jīng)營的重要基地。因此中部干旱帶苦蕎試點選擇鹽池縣蕎麥品種試驗點比較理想。

    2.2.4 參試品種生態(tài)適應(yīng)性

    利用R語言中的GGE-Biplot軟件包對5個品種在6個試點數(shù)據(jù)進(jìn)行生態(tài)適應(yīng)性分析(圖4)可以看出,在GGE雙標(biāo)圖的橫坐標(biāo)品種基因與環(huán)境試點的互作成分(PCA1)與縱坐標(biāo)互作成分(PCA2)可以解釋78.19%和10.51%,使基因(G)與環(huán)境(E)互作變異信息達(dá)到88.7%。

    從參試品種在不同生態(tài)類型區(qū)適應(yīng)性定位看,以中心圓點為坐標(biāo)軸,圍繞縱橫坐標(biāo)距離原點將所有品種與試點形成散狀分布點,以最遠(yuǎn)的品種與試點連成一個多邊形,過原點形成多邊扇形圖[22-23],將雙標(biāo)圖內(nèi)的品種和試點分成若干個區(qū)域,處于多邊形頂端的品種則是各區(qū)域內(nèi)表現(xiàn)最好的品種。由圖4看出,第一扇區(qū)為品種KV3(云蕎2號)和KV2(川蕎3號)適合在試點E4(海原縣)、E5(同心縣)和E6(鹽池縣)推廣種植;第二扇區(qū)品種KV1(晉苦蕎2號)生態(tài)適應(yīng)性最好,適合在試點E1(原州區(qū))、E2(西吉縣)、E3(彭陽縣)地區(qū)進(jìn)行大面積種植。

    2.2.5 不同降水年份苦蕎品種產(chǎn)量差異及適應(yīng)性

    苦蕎品種在不同試點和不同降水年份條件下的適應(yīng)性差異明顯(表6),對品種×年份產(chǎn)量交互效應(yīng)差異及適應(yīng)性進(jìn)行了分析。由表6可知,品種KV1(晉苦蕎2號)在干旱年、正常年和豐水年的產(chǎn)量最高,分別為1 671~2 593 kg·hm-2,平均產(chǎn)量2 187.6 kg·hm-2,較KV5(CK)品種增產(chǎn)17.3%。產(chǎn)量排第二位的品種KV2(川蕎3號)產(chǎn)量為1 474~2 513 kg·hm-2,平均產(chǎn)量2 004.7 kg·hm-2,較對照品種增產(chǎn)7.5%。而品種KV3(云蕎2號)、KV4(黔黑蕎5號)產(chǎn)量分別為1 465~2 310 kg·hm-2和1 358~2 483 kg·hm-2,較CK品種增產(chǎn)4.2%~4.9%。所有參試品種干旱年份較豐水年份減產(chǎn)28.0%~45.3%。

    品種KV1(晉苦蕎2號)和KV2(川蕎3號)在正常年份和豐水年份的產(chǎn)量分別達(dá)到2 299~2 593 kg·hm-2和2 027~2 513 kg·hm-2,說明這兩個品種在E1~E6區(qū)可進(jìn)行大面積推廣種植。KV3(云蕎2號)為中晚熟品種,適宜在熱量條件較好的E4(海原縣)和E6(鹽池縣)豐水年份種植,KV4(黔黑蕎5號)僅適宜在半干旱區(qū)E1(原州區(qū))和E3(彭陽縣)水年份種植,干旱年份不宜種植,而KV5(固原苦蕎)可選擇性地在半干旱區(qū)種植。

    3 討論

    作物品種的穩(wěn)產(chǎn)性和適應(yīng)性是決定其推廣應(yīng)用價值的重要指標(biāo)[24],本文基于AMMI模型和GGE雙標(biāo)圖分析了寧南山區(qū)的半干旱區(qū)、半干旱區(qū)易旱區(qū)和中部干旱區(qū)的5個苦蕎品種在6個試點環(huán)境的豐產(chǎn)穩(wěn)定性,適應(yīng)性和試點代表性。試驗經(jīng)受了干旱年份、正常年份和豐水年份,AMMI綜合了主成分分析和方差分析的優(yōu)點是將基因與環(huán)境(G×E)效應(yīng)中選取達(dá)到顯著的主成分PCA1、PCA2和PCA3互作下的穩(wěn)定性值(Dg),一定程度提高了品種豐產(chǎn)性和穩(wěn)定性參數(shù)估計精度。Dg值越小,表示品種的穩(wěn)定性越好,參試品種Dg值依次為,晉苦蕎2號7.80、云蕎2號11.80、川蕎3號12.97、固原當(dāng)?shù)乜嗍w14.76和黔黑蕎5號17.57。品種的豐產(chǎn)性穩(wěn)定性與品種在多年多點平均產(chǎn)量排名基本一致。對參試地點互作主成分值(PCA)與試點辨別力的穩(wěn)定性(De)進(jìn)行了分析,其De值越大,表明試點的穩(wěn)定性較好。代表苦蕎品種的區(qū)試點穩(wěn)定性De值依次為,E4(海原縣)18.32>E5(同心縣)14.65>E2(西吉縣)11.96>E3(彭陽縣)9.20>E6(鹽池縣)8.54>E1(原州區(qū))6.64。依據(jù)苦蕎品種的穩(wěn)定性(Dg)與試點辨別力的穩(wěn)定性(De)綜合判斷,篩選出適宜在寧夏南部山區(qū)種植的品種為晉苦蕎2號和川蕎3號,其豐產(chǎn)性和穩(wěn)定性最好,云蕎2號和黔黑蕎5號品種的豐產(chǎn)性、穩(wěn)定性適應(yīng)性次之。試點辨別力和代表性以半干旱易旱區(qū)的海原縣,中部干旱帶的同心縣和半干旱區(qū)的西吉縣比較好,其次為原州區(qū)和鹽池縣。

    生產(chǎn)中推廣應(yīng)用的苦蕎品種,必需要考慮其品種在某個地區(qū)的豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性。本研究認(rèn)為,試點環(huán)境效應(yīng)(E)是影響起品種產(chǎn)量變化的主要原因,遠(yuǎn)高于基因基因(G)及基因與環(huán)境(G×E)互作效應(yīng),其研究結(jié)果與前人研究部分結(jié)論基本一致[25-26]。如5個苦蕎品種在寧南山區(qū)不同試點環(huán)境下進(jìn)行多年試驗,其環(huán)境效應(yīng)(E)占總變異平方和56.3%>基因(G)24.0%>基因與環(huán)境(G×E)互作效應(yīng)19.8%。由此看出,試點環(huán)境是影響品種適應(yīng)性的主要因素。GGE坐標(biāo)圖能夠充分考慮環(huán)境對基因型的影響,已被廣泛應(yīng)用于作物多年多點的試驗分析中。因此,AMMI方法能更準(zhǔn)確地解釋環(huán)境和基因的互作,對提高品種豐產(chǎn)性穩(wěn)定性判斷和綜合評價將起到積極作用。

    本研究參試品種經(jīng)歷了干旱年份、正常年份和豐水年份,因此利用AMMI模型和GGE雙標(biāo)圖對品種的豐產(chǎn)性、穩(wěn)產(chǎn)性、生態(tài)適應(yīng)性和試點代表性分析及綜合評價對研究該地區(qū)苦蕎品種對不同氣候的響應(yīng)有一定指導(dǎo)意義。但就苦蕎作物品種在不同氣候類型區(qū)生長期對氣溫、光照、降水量,耗水量與生長期水分滿足程度及品種的豐產(chǎn)性穩(wěn)定性的影響有待進(jìn)一步研究。另外,苦蕎品種在同一試點不同年份間也存在較大的差異,為了更準(zhǔn)確地反映品種實質(zhì),在推廣之前要進(jìn)行多年多點生產(chǎn)性能和適應(yīng)性試驗,然后根據(jù)品種的豐產(chǎn)穩(wěn)定性和適應(yīng)性,結(jié)合當(dāng)?shù)氐纳鷳B(tài)環(huán)境進(jìn)行品種的綜合評價,充分發(fā)揮氣候資源對生產(chǎn)潛力的貢獻(xiàn),才能實現(xiàn)高產(chǎn)高效。

    4 結(jié)論

    利用GGE雙坐標(biāo)圖和AMMI 模型依據(jù)苦蕎品種的穩(wěn)定性(Dg)與試點辨別力的穩(wěn)定性(De)綜合判斷,篩選出適宜在寧夏南部山區(qū)種植的品種為晉苦蕎2號和川蕎3號,其豐產(chǎn)性和穩(wěn)定性最好。試點辨別力和代表性以半干旱易旱區(qū)的海原縣,中部干旱帶的同心縣和半干旱區(qū)的西吉縣比較好,其次為原州區(qū)和鹽池縣。

    不同降水年份對參試品種產(chǎn)量差異及適應(yīng)性分析結(jié)果表明,晉苦蕎2號在經(jīng)歷干旱年、正常年和豐水年產(chǎn)量為1 671~2 593 kg·hm-2,較對照品種增產(chǎn)17.3%,其次為川蕎3號產(chǎn)量為1 474~2 513 kg·hm-2,較對照品種增產(chǎn)7.5%,而品種KV3(云蕎2號)和KV4(黔黑蕎5號)產(chǎn)量分別為1 465~2 310 kg·hm-2和1 358~2 483 kg·hm-2,較CK品種增產(chǎn)4.2%~4.9%,所有參試品種在干旱年份較豐水年份減產(chǎn)28.0%~45.3%。

    參考文獻(xiàn):

    [1]李旭霞,劉晶,師靜.寧夏地區(qū)蕎麥栽培現(xiàn)狀調(diào)查研究[J].安徽農(nóng)學(xué)通報,2023,29(10):20-22.

    [2]謝銳,路戰(zhàn)遠(yuǎn),金曉蕾,等.基于AMMI模型的內(nèi)蒙古中西部地區(qū)馬鈴薯品種(系)穩(wěn)定性和適應(yīng)性分析[J].北方農(nóng)業(yè)學(xué)報,2022,50(1):18-24.

    [3]黃大明,謝雄澤,白桂萍,等.AMMI模型和GGE雙標(biāo)圖在湖北省油菜區(qū)域試驗中的應(yīng)用[J].湖北農(nóng)業(yè)科學(xué),2018,57(12):24-29.

    [4]步清,魯月,郝德榮,等.AMMI模型和GGE雙標(biāo)圖在江蘇省糯玉米品種區(qū)域試驗中的應(yīng)用[J].分子植物育種,2022,20(22):7628-7636.

    [5]吳曹陽,梁詩涵,邱軍,等.基于連續(xù)12年國家苦蕎區(qū)域試驗的中國苦蕎品種選育現(xiàn)狀分析[J].中國農(nóng)業(yè)科學(xué),2020,53(19):3878-3894.

    [6]魏常敏,許衛(wèi)猛,邢永鋒,等.AMMI模型和GGE雙標(biāo)圖在黃淮海甜玉米中的應(yīng)用[J].分子植物育種,2021,19(17):5909-5916.

    [7]郭敏杰,鄧麗,任麗,等.基于R語言的AMMI和GGE雙標(biāo)圖在花生區(qū)試中的應(yīng)用[J].花生學(xué)報,2017,46(2):24-31.

    [8]袁金海,賈多英,劉續(xù)漢,等.GGE雙標(biāo)圖在馬鈴薯產(chǎn)量及試驗點評價分析中的應(yīng)用[J].中國種業(yè),2018(8):56-59.

    [9]李亞杰,李德明,范士杰,等.GGE雙標(biāo)圖在馬鈴薯品種適應(yīng)性及產(chǎn)量穩(wěn)定性分析中的應(yīng)用評價[J].蘭州大學(xué)學(xué)報(自然科學(xué)版),2016,52(5):617-622.

    [10]宋慧,劉金榮,王素英,等.GGE雙標(biāo)圖評價谷子‘豫谷18’的豐產(chǎn)穩(wěn)產(chǎn)性和適應(yīng)性[J].中國農(nóng)業(yè)大學(xué)學(xué)報,2020,25(1):29-38.

    [11]昝凱,陳亞光,徐淑霞,等.基于R語言的GGE雙標(biāo)圖在大豆區(qū)試中的應(yīng)用[J].大豆科技,2019(4):16-20.

    [12]唐啟義,馮明光.DPS數(shù)據(jù)處理系統(tǒng):實驗設(shè)計、統(tǒng)計分析及數(shù)據(jù)挖掘[M].北京:科學(xué)出版社,2007.

    [13]王育才,王化俊,李葆春,等.基于GGE雙標(biāo)圖對青稞區(qū)試點代表性及品系穩(wěn)定性的分析[J].草原與草坪,2018,38(5):22-27.

    [14]王迪.GGE雙標(biāo)圖分析黑龍江省甜瓜區(qū)域試驗品種適應(yīng)性和穩(wěn)定性[J].黑龍江農(nóng)業(yè)科學(xué),2021(10):64-68.

    [15]張志芬,任長忠,楊海順,等.基于SD-GGE雙標(biāo)圖進(jìn)行飼草燕麥品種歸類和區(qū)域劃分[J].中國農(nóng)業(yè)大學(xué)學(xué)報,2021,26(5):10-21.

    [16]張大愛,趙緒明,錢一萍,等.用GGE雙標(biāo)圖分析苦蕎品種的產(chǎn)量穩(wěn)定性及試驗地點相似性[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2015,43(7):101-108.

    [17]李琴琴,高乙萍,張志芬,等.甜蕎品種穩(wěn)定性和試驗地點相似性的GGE雙標(biāo)圖分析[J].干旱地區(qū)農(nóng)業(yè)研究,2013,31(1):67-70,121.

    [18]魏鵬程,陳杜,羅英艦,等.基于AMMI和GGE雙標(biāo)圖的貴州不同生態(tài)區(qū)春玉米豐產(chǎn)性穩(wěn)產(chǎn)性和試點辨別力評價[J].玉米科學(xué),2023,31(1):22-31.

    [19]劉志芳,王廣才,劉兵兵,等.基于AMMI模型和GGE雙標(biāo)圖的河北省冬小麥品系穩(wěn)定性及適應(yīng)性分析[J].東北農(nóng)業(yè)科學(xué),2023,48(3):26-32.

    [20]李建武,李高峰,文國宏,等.甘肅省多點聯(lián)合試驗馬鈴薯產(chǎn)量要素穩(wěn)定性及試點鑒別力分析[J].西北農(nóng)業(yè)學(xué)報,2022,31(11):1422-1434.

    [21]王星宇,程靜,高生,等.應(yīng)用AMMI模型和GGE雙標(biāo)圖評價裸燕麥品種在華北高寒區(qū)的適應(yīng)性[J].草業(yè)學(xué)報,2022,31(12):76-84.

    [22]肖繼兵,劉志,辛宗緒,等.基于GGE雙標(biāo)圖的谷子種質(zhì)資源耐旱性鑒定[J].華北農(nóng)學(xué)報,2022,37(5):87-96.

    [23]慕平,趙桂琴,柴繼寬.基于GGE-Biplot的甘肅省不同生態(tài)區(qū)燕麥生產(chǎn)性能及適應(yīng)性分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2015,23(6):705-712.

    [24]柴繼寬,慕平,趙桂琴.8個燕麥品種在甘肅的產(chǎn)量穩(wěn)定性及試點代表性研究[J].草地學(xué)報,2016,24(5):1100-1107.

    [25]崔同霞,張保軍,姚友旭,等.基于GGE雙標(biāo)圖法的甘肅省西片區(qū)域春小麥新品種的豐產(chǎn)性及適應(yīng)性評價[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報,2022,57(3):58-65.

    [26]岳海旺,魏建偉,卜俊周,等.河北省春播玉米品種產(chǎn)量和主要穗部性狀GGE雙標(biāo)圖分析[J].玉米科學(xué),2018,26(4):28-35.

    Analysis of Stable Yield Adaptability of Fagopyrum tataricum Varieties in Different Climate Regions of Ningxia Using AMMI Model and GGE Double Label Map

    MU Lanhai,CHANG Keqin,DU Yanping,YANG Chongqing,CHEN Yixin

    Abstract:In order to accurately evaluate the high-yield, stable yield, and ecological adaptability of the new varieties of Fagopyrum tataricum, the AMMI model and GGE double label map were used to analyze the mprehensively evaluate the high-yield, stable yield, and pilot representativeness of five Fagopyrum tataricum varieties in six counties (districts) of semi-arid areas, semi-arid prone areas, and central arid zones in southern Ningxia. By analyzed the stability (Dg) of principal component PCA and ariety, as well as the yield stability (De) of experimental points, suitable varieties with high yield stability and good adaptability for planting in the southern mountainous areas of Ningxia were selected as Jinkuqiao 2 and Chuanqiao 3.The yield of Jinkuqiao 2 was 1 671-2 593 kg·ha-1, with an average increase of 17.3% compared to the control variety. The yield of Chuanqiao 3 was 1 474-2 513 kg·ha-1, with an average increase of 7.5% compared to the control variety.The experimental sites were ideal for Haiyuan County, which represents the semi-arid and drought prone areas, Tongxin County in the central arid zone, and Xiji County in the semi-arid areas.The experiment went through drought years, normal years, and flood years and the tested varieties experienced a 28.0%—45.3% reduction in yield compared to flood years in drought years.

    Keywords:Fagopyrum tataricum varieties;AMMI model;GGE double label diagram; high yield; stable yield;ecological adaptability

    收稿日期:2023-08-23

    基金項目:國家重點研發(fā)計劃項目(2022YFD1602503);寧夏農(nóng)林科學(xué)院科技成果轉(zhuǎn)化院地合作項目(NNKZZCGZH-2021-04);寧夏科技攻關(guān)項目(YES-16-18)。

    第一作者:穆蘭海(1966-),男,高級農(nóng)藝師,從事蕎麥和燕麥作物新品種選育及高效節(jié)水種植技術(shù)研究與示范。E-mail:gymlhai@163.com。

    通信作者:??饲冢?965-),男,學(xué)士,研究員,從事蕎麥和燕麥作物新品種選育。E-mail:nxgychkq@163.com。

    国产伦精品一区二区三区视频9| av一本久久久久| 激情五月婷婷亚洲| 人人妻人人添人人爽欧美一区卜| 高清不卡的av网站| 久久午夜福利片| 热99国产精品久久久久久7| 少妇 在线观看| 久久精品国产亚洲网站| 免费看av在线观看网站| 97在线人人人人妻| 欧美97在线视频| 国产毛片在线视频| 三级经典国产精品| 你懂的网址亚洲精品在线观看| 日本黄大片高清| 午夜日本视频在线| 看非洲黑人一级黄片| 久久久久久久久大av| 少妇丰满av| 少妇人妻 视频| 国产精品国产三级国产av玫瑰| 日韩中文字幕视频在线看片| 亚洲精品国产成人久久av| 国产精品女同一区二区软件| 日韩欧美精品免费久久| 亚洲综合精品二区| 51国产日韩欧美| 香蕉精品网在线| 99久久人妻综合| 日韩视频在线欧美| 欧美人与善性xxx| h日本视频在线播放| 欧美三级亚洲精品| 日产精品乱码卡一卡2卡三| 中文欧美无线码| 如日韩欧美国产精品一区二区三区 | 久久97久久精品| 能在线免费看毛片的网站| a 毛片基地| 嘟嘟电影网在线观看| 黄色欧美视频在线观看| av专区在线播放| 欧美激情极品国产一区二区三区 | 99热6这里只有精品| 热re99久久精品国产66热6| 亚洲av在线观看美女高潮| 毛片一级片免费看久久久久| 最近手机中文字幕大全| 99热6这里只有精品| 欧美成人午夜免费资源| 久久精品熟女亚洲av麻豆精品| 十八禁网站网址无遮挡 | 亚洲国产最新在线播放| 狂野欧美白嫩少妇大欣赏| 女的被弄到高潮叫床怎么办| 久久久久久久国产电影| 大又大粗又爽又黄少妇毛片口| 国产日韩欧美在线精品| 久久亚洲国产成人精品v| 成年人免费黄色播放视频 | 国国产精品蜜臀av免费| 六月丁香七月| 51国产日韩欧美| 伊人久久国产一区二区| 久久亚洲国产成人精品v| 一级毛片黄色毛片免费观看视频| 国产精品福利在线免费观看| 亚洲精品乱久久久久久| 美女xxoo啪啪120秒动态图| 亚洲婷婷狠狠爱综合网| 少妇人妻精品综合一区二区| 中文字幕制服av| 国产成人精品一,二区| 亚洲国产精品999| 亚洲av二区三区四区| 伊人亚洲综合成人网| 欧美日韩精品成人综合77777| 久久久国产一区二区| 一级二级三级毛片免费看| 亚洲精品国产av成人精品| 丰满饥渴人妻一区二区三| av在线app专区| 中文资源天堂在线| 在线天堂最新版资源| 日产精品乱码卡一卡2卡三| 97超视频在线观看视频| 国产男人的电影天堂91| 97在线视频观看| 大片免费播放器 马上看| 国产中年淑女户外野战色| 亚洲国产精品国产精品| 国产午夜精品一二区理论片| 午夜福利,免费看| 久久鲁丝午夜福利片| 性高湖久久久久久久久免费观看| 成人国产av品久久久| 久久99热6这里只有精品| 中文资源天堂在线| av女优亚洲男人天堂| 自线自在国产av| 成人美女网站在线观看视频| 精品一品国产午夜福利视频| 性高湖久久久久久久久免费观看| 国产永久视频网站| 亚洲国产精品专区欧美| 欧美日韩视频高清一区二区三区二| 亚洲欧美清纯卡通| 嘟嘟电影网在线观看| 免费人成在线观看视频色| 国产在线视频一区二区| 人妻夜夜爽99麻豆av| 日韩精品有码人妻一区| 亚洲欧洲日产国产| 国精品久久久久久国模美| 日韩 亚洲 欧美在线| 97在线人人人人妻| 在线精品无人区一区二区三| 制服丝袜香蕉在线| 亚洲欧洲国产日韩| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品国产色婷婷电影| 午夜老司机福利剧场| 国产精品国产三级国产专区5o| 国产日韩欧美在线精品| 人人妻人人澡人人爽人人夜夜| 91aial.com中文字幕在线观看| 高清av免费在线| 美女内射精品一级片tv| 韩国av在线不卡| 国产色婷婷99| av天堂中文字幕网| 蜜臀久久99精品久久宅男| 亚洲伊人久久精品综合| 久久影院123| 亚洲四区av| 亚洲人成网站在线播| 两个人免费观看高清视频 | 婷婷色综合www| 中文字幕免费在线视频6| 国产精品不卡视频一区二区| 久久99一区二区三区| 亚洲伊人久久精品综合| 日本爱情动作片www.在线观看| 精品一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 日韩在线高清观看一区二区三区| 免费观看在线日韩| 久久精品国产自在天天线| 国产视频首页在线观看| 日韩亚洲欧美综合| 高清毛片免费看| 伦精品一区二区三区| a级毛片在线看网站| 亚洲综合色惰| 看非洲黑人一级黄片| 777米奇影视久久| 成人无遮挡网站| 国产精品蜜桃在线观看| 少妇被粗大的猛进出69影院 | 一二三四中文在线观看免费高清| 免费黄色在线免费观看| 亚洲av成人精品一区久久| 女性生殖器流出的白浆| 亚洲精品,欧美精品| 黑人猛操日本美女一级片| 日本欧美国产在线视频| av黄色大香蕉| 人人妻人人澡人人爽人人夜夜| 亚洲美女黄色视频免费看| 男女国产视频网站| av女优亚洲男人天堂| 成年av动漫网址| av有码第一页| 男女啪啪激烈高潮av片| 亚洲av福利一区| 性色avwww在线观看| 久久久久久久大尺度免费视频| 日韩在线高清观看一区二区三区| 久久久久久久久久久久大奶| 一本—道久久a久久精品蜜桃钙片| 六月丁香七月| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 国产真实伦视频高清在线观看| 成人影院久久| 黄色视频在线播放观看不卡| 精品一区二区三卡| 波野结衣二区三区在线| av在线观看视频网站免费| 亚洲精品乱久久久久久| 黄色怎么调成土黄色| 国产黄色视频一区二区在线观看| 欧美一级a爱片免费观看看| av.在线天堂| 少妇人妻久久综合中文| 一区二区三区四区激情视频| 欧美3d第一页| 久久亚洲国产成人精品v| 只有这里有精品99| 一级毛片电影观看| 丝袜脚勾引网站| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花| 亚洲欧美成人综合另类久久久| 高清不卡的av网站| 久热久热在线精品观看| 高清视频免费观看一区二区| 黄片无遮挡物在线观看| 精品久久国产蜜桃| 在线 av 中文字幕| 99热这里只有是精品50| 国产熟女欧美一区二区| 黑人高潮一二区| 免费观看的影片在线观看| 久久久a久久爽久久v久久| 成人亚洲欧美一区二区av| 搡女人真爽免费视频火全软件| 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| 久久av网站| av播播在线观看一区| 人人澡人人妻人| 色视频www国产| 91精品一卡2卡3卡4卡| 久久热精品热| 丰满乱子伦码专区| 51国产日韩欧美| 91精品一卡2卡3卡4卡| 一区二区三区精品91| 99九九线精品视频在线观看视频| 韩国高清视频一区二区三区| 尾随美女入室| 少妇人妻一区二区三区视频| 亚洲欧美日韩卡通动漫| 大香蕉97超碰在线| 亚洲图色成人| 一个人看视频在线观看www免费| 久久免费观看电影| 女人精品久久久久毛片| 99久久综合免费| 日本黄色片子视频| 赤兔流量卡办理| 美女国产视频在线观看| 欧美 日韩 精品 国产| 成人毛片a级毛片在线播放| 高清不卡的av网站| 人人妻人人爽人人添夜夜欢视频 | 嘟嘟电影网在线观看| 成人毛片a级毛片在线播放| 国产高清有码在线观看视频| 亚洲精品第二区| 日韩在线高清观看一区二区三区| 汤姆久久久久久久影院中文字幕| 久久ye,这里只有精品| 99热全是精品| 亚洲精品一二三| 自拍偷自拍亚洲精品老妇| 欧美性感艳星| 亚洲av不卡在线观看| 免费av不卡在线播放| 亚洲国产av新网站| 精品国产国语对白av| 午夜影院在线不卡| av网站免费在线观看视频| 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱| 久久国产乱子免费精品| 午夜福利视频精品| 国产亚洲av片在线观看秒播厂| 免费观看av网站的网址| 午夜久久久在线观看| 下体分泌物呈黄色| 91精品一卡2卡3卡4卡| 欧美日韩视频精品一区| 亚洲高清免费不卡视频| 青春草视频在线免费观看| 日本黄大片高清| 欧美另类一区| 国产亚洲一区二区精品| 免费看光身美女| 亚洲国产精品专区欧美| 女性被躁到高潮视频| 亚洲人成网站在线观看播放| 亚洲美女黄色视频免费看| 一二三四中文在线观看免费高清| 极品少妇高潮喷水抽搐| 伊人亚洲综合成人网| 99九九线精品视频在线观看视频| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99| 日韩三级伦理在线观看| 国产精品一区www在线观看| a级片在线免费高清观看视频| 日本午夜av视频| 大片免费播放器 马上看| 99热这里只有是精品50| 在线观看av片永久免费下载| 久久人人爽av亚洲精品天堂| 欧美日韩亚洲高清精品| 高清午夜精品一区二区三区| 在线观看三级黄色| 99热6这里只有精品| 黄色配什么色好看| www.色视频.com| 五月天丁香电影| 亚洲欧美成人精品一区二区| 成人黄色视频免费在线看| 精品国产国语对白av| 久久这里有精品视频免费| 黄片无遮挡物在线观看| 亚洲成人手机| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 能在线免费看毛片的网站| h视频一区二区三区| 在线观看免费日韩欧美大片 | 亚洲va在线va天堂va国产| 久热这里只有精品99| 街头女战士在线观看网站| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 人人妻人人看人人澡| 寂寞人妻少妇视频99o| 18禁裸乳无遮挡动漫免费视频| 黄色视频在线播放观看不卡| 久久精品国产亚洲av涩爱| 各种免费的搞黄视频| 春色校园在线视频观看| a级片在线免费高清观看视频| 国产欧美日韩综合在线一区二区 | 少妇熟女欧美另类| 大香蕉久久网| 国产男人的电影天堂91| 热re99久久精品国产66热6| 亚洲成人av在线免费| 夫妻午夜视频| 国产精品久久久久久精品古装| 夫妻性生交免费视频一级片| 男女边吃奶边做爰视频| 中文字幕亚洲精品专区| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 十八禁高潮呻吟视频 | 99热全是精品| 亚洲成人av在线免费| 青春草亚洲视频在线观看| 国产深夜福利视频在线观看| 欧美精品人与动牲交sv欧美| 精品少妇久久久久久888优播| tube8黄色片| 自拍偷自拍亚洲精品老妇| 午夜福利影视在线免费观看| 另类亚洲欧美激情| 人人妻人人澡人人看| 最近的中文字幕免费完整| 国产白丝娇喘喷水9色精品| 777米奇影视久久| 香蕉精品网在线| 国产免费一级a男人的天堂| 啦啦啦中文免费视频观看日本| 久久久国产欧美日韩av| 日日啪夜夜撸| 好男人视频免费观看在线| 久久婷婷青草| 国产欧美日韩综合在线一区二区 | 91精品一卡2卡3卡4卡| 亚洲av.av天堂| 水蜜桃什么品种好| 国产欧美亚洲国产| 精品久久国产蜜桃| 欧美成人午夜免费资源| 亚洲av国产av综合av卡| 一区二区三区乱码不卡18| 日日摸夜夜添夜夜添av毛片| 丰满人妻一区二区三区视频av| 久久精品久久久久久久性| 肉色欧美久久久久久久蜜桃| 成人漫画全彩无遮挡| 少妇精品久久久久久久| 成人黄色视频免费在线看| 精品一区二区三区视频在线| 大香蕉久久网| 国产视频首页在线观看| 在线观看人妻少妇| 在线观看免费日韩欧美大片 | 人妻系列 视频| 中文字幕人妻丝袜制服| 免费大片18禁| 久久久久久久亚洲中文字幕| 边亲边吃奶的免费视频| 日韩av免费高清视频| 内地一区二区视频在线| 女人久久www免费人成看片| 日本vs欧美在线观看视频 | av免费观看日本| 水蜜桃什么品种好| 国产男女超爽视频在线观看| 国产日韩一区二区三区精品不卡 | 性色av一级| 天堂中文最新版在线下载| 色视频在线一区二区三区| 久久精品熟女亚洲av麻豆精品| 人人妻人人澡人人看| 国产精品国产三级国产专区5o| 欧美三级亚洲精品| 亚洲欧美日韩卡通动漫| 另类亚洲欧美激情| 中国三级夫妇交换| 香蕉精品网在线| 91午夜精品亚洲一区二区三区| 精品久久久精品久久久| 国产伦在线观看视频一区| 国产亚洲91精品色在线| 自线自在国产av| 日本与韩国留学比较| 精品一区二区三区视频在线| 免费av不卡在线播放| 久久综合国产亚洲精品| 免费人成在线观看视频色| 精品久久久久久电影网| 亚洲怡红院男人天堂| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区 | 国产精品国产av在线观看| 亚洲欧洲国产日韩| 色吧在线观看| 日本91视频免费播放| 18禁裸乳无遮挡动漫免费视频| 亚洲图色成人| 又粗又硬又长又爽又黄的视频| 精品一品国产午夜福利视频| 成人18禁高潮啪啪吃奶动态图 | www.av在线官网国产| 国产精品欧美亚洲77777| 日本欧美国产在线视频| 曰老女人黄片| 成年女人在线观看亚洲视频| 大码成人一级视频| 亚洲成色77777| 亚洲av电影在线观看一区二区三区| 亚洲不卡免费看| 久久国产乱子免费精品| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费 | 欧美少妇被猛烈插入视频| 久久这里有精品视频免费| 国产亚洲午夜精品一区二区久久| 亚洲不卡免费看| 五月伊人婷婷丁香| 寂寞人妻少妇视频99o| 亚洲成人手机| 男女免费视频国产| 免费黄色在线免费观看| 久久精品国产自在天天线| av在线app专区| 九色成人免费人妻av| 熟女av电影| 久久99热这里只频精品6学生| 性高湖久久久久久久久免费观看| 日本欧美国产在线视频| 只有这里有精品99| 亚洲图色成人| 一级片'在线观看视频| 亚洲精品一二三| 亚洲欧美日韩卡通动漫| 亚洲欧美精品自产自拍| 天堂俺去俺来也www色官网| 99热国产这里只有精品6| 亚洲国产毛片av蜜桃av| 我要看日韩黄色一级片| 国产在线免费精品| 美女cb高潮喷水在线观看| 日韩一本色道免费dvd| 国产一级毛片在线| 成人国产av品久久久| 黄色配什么色好看| 大香蕉97超碰在线| 自拍偷自拍亚洲精品老妇| 亚洲不卡免费看| 欧美精品国产亚洲| 在线精品无人区一区二区三| 国语对白做爰xxxⅹ性视频网站| 国内少妇人妻偷人精品xxx网站| 成人二区视频| 97超视频在线观看视频| 熟女av电影| 六月丁香七月| 精品人妻偷拍中文字幕| 中文字幕亚洲精品专区| 欧美精品国产亚洲| 国产探花极品一区二区| 国产精品久久久久成人av| 男的添女的下面高潮视频| 美女视频免费永久观看网站| 日韩三级伦理在线观看| 亚洲av成人精品一区久久| 爱豆传媒免费全集在线观看| av在线老鸭窝| 精品久久久久久电影网| 蜜桃久久精品国产亚洲av| 中文欧美无线码| 妹子高潮喷水视频| 日韩一本色道免费dvd| 久久久久网色| av.在线天堂| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 国模一区二区三区四区视频| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频| 日韩熟女老妇一区二区性免费视频| 国产免费视频播放在线视频| av有码第一页| 18禁在线无遮挡免费观看视频| 老司机影院毛片| 欧美 亚洲 国产 日韩一| 中文天堂在线官网| 国产日韩欧美视频二区| 午夜日本视频在线| 在线观看www视频免费| 看非洲黑人一级黄片| av黄色大香蕉| 2018国产大陆天天弄谢| 欧美老熟妇乱子伦牲交| 校园人妻丝袜中文字幕| h日本视频在线播放| 99久久精品国产国产毛片| 中文天堂在线官网| 色视频www国产| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 国产亚洲91精品色在线| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 人妻 亚洲 视频| 欧美少妇被猛烈插入视频| 免费看不卡的av| 男女国产视频网站| 九九爱精品视频在线观看| 久久99蜜桃精品久久| 国产黄片视频在线免费观看| 人妻一区二区av| 观看美女的网站| 欧美区成人在线视频| 国产精品国产三级国产专区5o| 久久久久久久久大av| 中国美白少妇内射xxxbb| 啦啦啦中文免费视频观看日本| 亚洲精品aⅴ在线观看| 亚洲人成网站在线播| 黑人猛操日本美女一级片| av福利片在线观看| 青春草国产在线视频| av免费在线看不卡| 久久久久久久久久成人| 国产片特级美女逼逼视频| 日本午夜av视频| 成年女人在线观看亚洲视频| 精品少妇内射三级| 亚洲精品久久久久久婷婷小说| 伊人久久国产一区二区| 亚洲情色 制服丝袜| 永久网站在线| 亚洲精品国产av蜜桃| 一级毛片电影观看| 人体艺术视频欧美日本| 国产精品久久久久久精品古装| 午夜激情久久久久久久| 看十八女毛片水多多多| 99热全是精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩东京热| 男人爽女人下面视频在线观看| 日韩中字成人| 男人添女人高潮全过程视频| 在线观看免费高清a一片| 这个男人来自地球电影免费观看 | 国产高清不卡午夜福利| 性色av一级| 久久久国产一区二区| h日本视频在线播放| 国产视频内射| 黑人高潮一二区| 岛国毛片在线播放| 亚洲精品成人av观看孕妇| 各种免费的搞黄视频| 欧美激情极品国产一区二区三区 | 精品一区二区三区视频在线| 搡女人真爽免费视频火全软件| 在线播放无遮挡| 色婷婷久久久亚洲欧美| 亚洲av成人精品一二三区| 国产精品不卡视频一区二区| 在现免费观看毛片| 日本wwww免费看| 在线播放无遮挡| 五月天丁香电影| 99久久精品热视频| 日韩一区二区视频免费看| 日韩成人av中文字幕在线观看| 亚洲成人av在线免费| 一本一本综合久久| 天天操日日干夜夜撸| a级毛色黄片| 亚洲天堂av无毛| 高清毛片免费看| 一本色道久久久久久精品综合| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 日韩人妻高清精品专区| 免费在线观看成人毛片|