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

    賀蘭山東麓‘赤霞珠’品質(zhì)形成氣象條件與評級方法研究*

    2018-03-07 03:38:03馬力文李劍萍韓穎娟李萬春

    馬力文, 李劍萍, 韓穎娟, 李萬春

    ?

    賀蘭山東麓‘赤霞珠’品質(zhì)形成氣象條件與評級方法研究*

    馬力文, 李劍萍**, 韓穎娟, 李萬春

    (中國氣象局旱區(qū)特色農(nóng)業(yè)氣象災(zāi)害監(jiān)測預(yù)警與風(fēng)險管理重點實驗室/寧夏氣象防災(zāi)減災(zāi)重點實驗室 銀川 750002)

    葡萄品質(zhì)很大程度上決定了葡萄酒品質(zhì), 其與氣象條件關(guān)系密切。在正常管理和農(nóng)業(yè)技術(shù)水平下,衡量葡萄生長氣象條件是否適宜品質(zhì)形成的方法稱為釀酒葡萄品質(zhì)氣象評價技術(shù)。研究釀酒葡萄品質(zhì)形成的氣象評級技術(shù)和氣象指標(biāo), 可為葡萄酒商業(yè)評級、年份酒鑒定和窖藏提供氣候參考依據(jù)。利用賀蘭山東麓2003—2011年多點采樣化驗的101份‘赤霞珠’總糖、總酸、糖酸比、pH和單寧資料, 選取其中的51份進(jìn)行相關(guān)普查, 確定了有生物學(xué)意義的顯著氣象因子, 以各因子的決定系數(shù)分別建立了5項品質(zhì)指標(biāo)的權(quán)重模型, 通過單項品質(zhì)指標(biāo)對綜合品質(zhì)的貢獻(xiàn)構(gòu)建了綜合品質(zhì)氣象評價模型。參照前人對葡萄酒品質(zhì)、釀酒葡萄品質(zhì)的研究和葡萄酒和葡萄品質(zhì)標(biāo)準(zhǔn), 根據(jù)釀造葡萄酒所需的總糖、總酸、糖酸比、pH和單寧, 確定‘赤霞珠’各項品質(zhì)指標(biāo)的閾值, 利用模型反推出與葡萄酒品質(zhì)相對應(yīng)的氣象因子閾值和‘赤霞珠’品質(zhì)氣象分類標(biāo)準(zhǔn)?;卮鷻z驗顯示5項品質(zhì)指標(biāo)均通過了0.001的檢驗和檢驗, 與實測值接近且變化趨勢一致??偺?、總酸、糖酸比和單寧的模擬效果較好,≥0.59且RMSE較小, 但pH誤差相對較大。用50份未參與建模的樣本估算了‘赤霞珠’品質(zhì)等級??偺?、總酸與實測接近; 個別樣本的糖酸比、單寧誤差相對較大, 但變化趨勢與實測一致; 樣本pH均在適宜范圍內(nèi)且較穩(wěn)定。從綜合評分和等級來看, 28個樣本與實況等級相同, 18個樣本誤差在1級內(nèi), 僅4個樣本相差2級, 準(zhǔn)確反映出原料的質(zhì)量?!嘞贾椤焚|(zhì)氣象評價指標(biāo)和模型為評價賀蘭山東麓釀酒葡萄品質(zhì)提供了一套可行的方法。

    賀蘭山東麓; 赤霞珠; 釀酒葡萄; 氣象品質(zhì); 氣象因子; 品質(zhì)等級; 模型評價

    賀蘭山東麓光熱充足, 干燥少雨, 是業(yè)界認(rèn)可的釀酒葡萄世界最佳產(chǎn)區(qū), 葡萄酒屢獲國內(nèi)外大獎, 2003年獲“葡萄酒國家地理標(biāo)志產(chǎn)品”保護(hù)區(qū)認(rèn)證[1]。業(yè)界普遍認(rèn)為優(yōu)質(zhì)葡萄酒主要取決于品種、土壤和氣象條件[2], 氣象條件是決定釀造原料品質(zhì)的關(guān)鍵[3]。氣象條件好的年份葡萄酒的濃縮度高, 單寧、色素、酚類物質(zhì)、礦物質(zhì)、維生素等成分含量高[4]。因此, 開展釀酒葡萄品質(zhì)氣象評價可為葡萄酒評級、提高附加值和市場認(rèn)可度提供重要參考依據(jù)。

    研究釀酒葡萄品質(zhì), 首先要明確決定葡萄酒品質(zhì)的主要釀造指標(biāo)。李記明[5]提出葡萄酒品質(zhì)決定于葡萄總糖、總酸、糖酸比、pH和單寧含量??偺菦Q定了酒度和風(fēng)味, 170 g?L-1以上才能釀出優(yōu)質(zhì)葡萄酒[6]。總糖與成熟度有關(guān), 王華等[7]用水熱系數(shù)反映葡萄成熟度。Vr?i?等[8]、Neumann等[9]認(rèn)為氣溫升高導(dǎo)致葡萄總糖增加, 總酸下降。毛如志等[10]認(rèn)為氣溫高是云南‘美樂’葡萄總糖、單寧比山東高的原因, 但Webb等[11-12]研究表明澳大利亞‘赤霞珠’和‘黑比諾’品質(zhì)隨氣溫升高而下降??偹崤c口感有關(guān), 李華[13]確定釀酒葡萄的適宜酸度為6~10 g?L-1, 過少使葡萄酒乏味、少筋、平淡, 過多酸澀。翟衡等[14]認(rèn)為釀酒葡萄總酸決定于光合呼吸、漿果膨大的稀釋及糖、醛、酸的代謝轉(zhuǎn)化。釀酒葡萄糖酸比相對穩(wěn)定。劉玉蘭等[15-16]認(rèn)為糖酸比與采收前35 d水熱系數(shù)關(guān)系密切, 32左右較合適, 過高、過低不易釀造出優(yōu)質(zhì)葡萄酒。陳曉前等[17]認(rèn)為pH是衡量葡萄酒品質(zhì)的指標(biāo); 李記明等[6]給出釀造優(yōu)質(zhì)葡萄酒的pH為3.0~3.6, 過高會增加微生物活性, 降低花色苷的顯色比例和游離SO2的有效量, 縮短陳釀潛力。張曉煜等[18]認(rèn)為釀酒葡萄pH與開花至成熟期最高氣溫關(guān)系密切。岳俊波等[19]提出葡萄酒中的單寧能沉淀蛋白質(zhì), 提高結(jié)構(gòu)感, 穩(wěn)定色素, 抗氧化, 抗自由基, 抗菌, 防止還原味。葡萄品質(zhì)與土壤有關(guān), 周濤等[20]認(rèn)為土壤質(zhì)地影響賀蘭山東麓釀酒葡萄品質(zhì)。趙磊等[21]認(rèn)為賀蘭山東麓‘赤霞珠’和‘美樂’香氣物質(zhì)比河北沙城高的原因是兩地氣候、土壤的差異。王秀琴等[22]認(rèn)為各地氣候、土壤差異與葡萄果實花青素、糖、pH、可滴定酸、果皮和種子多酚及葡萄酒的醇含量、蘋果酸、總酚、顏色強度、顏色穩(wěn)定性及口感均有關(guān)系。Pircalabu等[23]、Jones等[24]認(rèn)為氣候變率導(dǎo)致葡萄風(fēng)味年際間差異明顯。此外, 還有一些成分影響葡萄酒品質(zhì)和風(fēng)味。如Cheng等[25]提出新疆≥35 ℃高溫日數(shù)偏多、降水偏少可增加釀酒葡萄花青素濃度。張娟等[26]、周淑珍等[27]認(rèn)為氣象條件影響葡萄酒總酚含量。楊曉帆等[28]認(rèn)為積溫與葡萄香氣物質(zhì)的積累有關(guān)。房玉林等[29]試驗得到水分虧缺對釀酒葡萄生長和品質(zhì)均有影響。一些國家嘗試用氣候評價葡萄年份酒。如Webb等[11-12]建立了優(yōu)質(zhì)葡萄酒品種的品質(zhì)氣候敏感性估測模型, Gouveia等[30]用葡萄牙春夏季最高氣溫、降水量和霜日建立了年份酒氣候評價方法, 張曉煜等[18,31]采用調(diào)和權(quán)重法構(gòu)建了賀蘭山東麓釀酒葡萄的綜合品質(zhì)權(quán)重模型。

    綜上所述, 國內(nèi)外對葡萄品質(zhì)與氣象、土壤養(yǎng)分的關(guān)系研究較多, 個別國家建立了評價葡萄年份酒的氣候模型, 但國內(nèi)尚未建立釀酒葡萄品質(zhì)的氣象評價模型。葡萄酒質(zhì)量70%決定于原料, 30%決定于工藝和設(shè)備。本文在歸納、借鑒前人研究的基礎(chǔ)上, 把不同時段的溫度、積溫、最高氣溫、最低氣溫、降水量、日照時數(shù)等氣象要素分別與總糖、總酸、糖酸比、pH和單寧進(jìn)行相關(guān)普查, 篩選出對品質(zhì)影響生物學(xué)意義明顯的顯著氣象因子, 根據(jù)其貢獻(xiàn)權(quán)重合成最佳模擬方程。參考釀酒葡萄品質(zhì)評價相關(guān)標(biāo)準(zhǔn)[32-35]及李記明等[5-6]的研究, 結(jié)合釀酒師的經(jīng)驗, 分別制定了‘赤霞珠’總糖、總酸、糖酸比、pH和單寧等級及關(guān)鍵氣象因子的指標(biāo)范圍, 構(gòu)建了‘赤霞珠’綜合品質(zhì)的氣象評價模型和指標(biāo), 可根據(jù)氣象條件評價釀造原料應(yīng)形成的品質(zhì)級別, 為保障生產(chǎn)優(yōu)質(zhì)葡萄酒、商業(yè)評級提供參考依據(jù)。

    1 材料與方法

    1.1 資料來源

    所用數(shù)據(jù)資料為2003—2011年賀蘭山東麓7個市縣的葡萄基地和酒莊的101份‘赤霞珠’化驗和生育進(jìn)程觀測結(jié)果。所涉及的葡萄生育期包括開花期、轉(zhuǎn)色期和采收期, 品質(zhì)化驗指標(biāo)有可溶性固形物、總糖、總酸、糖酸比、pH及單寧。采用斐林試劑滴定法測定總糖含量, NaOH滴定法測定總酸(以酒石酸計), 硫酸苯酚法測定單寧。選擇其中51份樣品建立模型, 另50份樣品用于檢驗效果。氣象資料為與化驗?zāi)攴莺筒蓸拥攸c對應(yīng)的自動氣象站逐日氣溫、降水量資料, 日照選用樣點所在市縣氣象局大監(jiān)站的資料(表1)。

    1.2 ‘赤霞珠’品質(zhì)指標(biāo)與氣象因子的膨化相關(guān)普查

    采用‘赤霞珠’可溶性固形物、總糖、總酸、糖酸比、pH及單寧化驗結(jié)果, 分別與生長期間不同階段平均氣溫、最高和最低氣溫、≥10 ℃積溫、累計降水量、日均降水量、累計日照時數(shù)等因子進(jìn)行膨化相關(guān)普查; 時段分為開花至采收、轉(zhuǎn)色至采收以及距離采收10~40 d內(nèi)間隔5 d, 以擴大因子的篩選范圍, 粗選出通過顯著性檢驗的氣象因子。

    1.3 ‘赤霞珠’品質(zhì)要素與氣象因子的關(guān)系模型

    分別研究可溶性固形物、總糖、總酸、糖酸比、pH及單寧與相關(guān)極顯著的氣象因子的線性、非線性關(guān)系, 選擇通過檢驗且對品質(zhì)形成具有生物學(xué)意義的氣象因子, 建立不同品質(zhì)指標(biāo)的關(guān)系模型, 根據(jù)入選氣象因子復(fù)相關(guān)系數(shù)的權(quán)重建立評價釀酒葡萄單項品質(zhì)要素的綜合權(quán)重模擬方程。

    1.4 釀酒葡萄單項品質(zhì)要素等級判別指標(biāo)及綜合品質(zhì)氣象等級模型

    參考4項葡萄酒和葡萄品質(zhì)評價標(biāo)準(zhǔn)、國際植物遺傳委員會(BPGR)葡萄品質(zhì)分級標(biāo)準(zhǔn)[32-35]及相關(guān)文獻(xiàn)[5-6], 結(jié)合玉泉營、立蘭酒莊釀酒師的分級經(jīng)驗, 分別確定評價釀酒葡萄總糖、總酸、糖酸比、pH和單寧優(yōu)劣的分級指標(biāo), 每項指標(biāo)均分為5級(特優(yōu))、4級(優(yōu))、3級(良)、2級(中)和1級(差)5個等級標(biāo)準(zhǔn), 反推顯著氣象因子的等級指標(biāo)范圍??紤]到以上5項品質(zhì)要素對釀酒葡萄綜合品質(zhì)等級的貢獻(xiàn), 以每個品質(zhì)要素氣象因子綜合模擬方程的2作為權(quán)重, 集成得到綜合品質(zhì)氣象等級模型。按照上述等級劃分方法, 將綜合品質(zhì)等級也劃分為5個級別。評價時將相關(guān)氣象因子代入每項品質(zhì)氣象模擬模型, 分別得到總糖、總酸、糖酸比、pH和單寧的估測值和相應(yīng)等級, 代入綜合品質(zhì)等級模型, 可評價釀酒葡萄綜合品質(zhì)等級。

    表1 ‘赤霞珠’葡萄采樣年份、地點及對應(yīng)的氣象站點位置

    1.5 品質(zhì)氣象評價模型的效果檢驗

    為檢驗該方法的適宜性, 選用2003—2011年期間永寧、青銅峽和蘆花臺葡萄基地的‘赤霞珠’化驗樣本, 將參與模擬的51份樣本進(jìn)行回代檢驗, 其余未參與模擬的50份樣本進(jìn)行評估檢驗。對估算值和實測值均按照分類閾值進(jìn)行品質(zhì)等級歸類, 與實測值歸類結(jié)果進(jìn)行比較, 計算單項品質(zhì)指標(biāo)評分和綜合評分的準(zhǔn)確率。

    2 結(jié)果與分析

    2.1 ‘赤霞珠’品質(zhì)與氣象因子的膨化相關(guān)普查

    膨化相關(guān)普查表明(表2), 可溶性固形物、總糖含量及糖酸比與開花至采收期的平均氣溫、最高氣溫極顯著負(fù)相關(guān), 與≥10 ℃積溫和日照正相關(guān); 轉(zhuǎn)色至采收期間上述顯著因子的相關(guān)更好, 表明轉(zhuǎn)色后常出現(xiàn)不利氣象條件, 影響總糖和可溶性固形物的累積及糖酸轉(zhuǎn)化。總酸含量與開花至采收期間平均氣溫正相關(guān), 與日照負(fù)相關(guān), 均通過0.001的檢驗, 而轉(zhuǎn)色至采收期間相關(guān)性弱一些, 說明總酸主要在結(jié)實前中期積累, 轉(zhuǎn)色后的氣象條件主要影響糖酸的轉(zhuǎn)化。pH僅與轉(zhuǎn)色至采收日照時數(shù)的正相關(guān)和采收前30 d降水量的負(fù)相關(guān)達(dá)到0.005極顯著水平。單寧與轉(zhuǎn)色至采收期間的積溫和日照正相關(guān), 與轉(zhuǎn)色后的平均氣溫、最低氣溫負(fù)相關(guān), 與采收前10~30 d的最低氣溫、降水量負(fù)相關(guān)更顯著。表明轉(zhuǎn)色后氣溫越低, 日照越少, 單寧越高; 轉(zhuǎn)色至采收間隔越長, 期間積溫和日照越多, 代表晴天多, 氣溫日較差大, 單寧越多; 期間降水越多, 單寧越少。

    表2 ‘赤霞珠’葡萄品質(zhì)與不同生育期氣象因子的相關(guān)分析(n=51)

    ****、***、**和*表示通過0.001、0.005、0.01和0.05的顯著性檢驗。****, ***, ** and * indicate passing= 0.001, 0.005, 0.01 and 0.05 significance test.

    2.2 ‘赤霞珠’品質(zhì)與氣象因子的關(guān)系

    2.2.1 總糖含量與氣象因子的關(guān)系

    選擇與總糖相關(guān)極顯著的氣象因子建立線性或非線性模型(表3)??偺桥c開花至采收期間日照時數(shù)呈極顯著正相關(guān), 日照越多, 糖分積累也越多。與≥10 ℃積溫呈二次曲線關(guān)系, 適時采收含糖量最高, 采收過早或過遲, 都會影響釀酒葡萄的含糖量??偺桥c轉(zhuǎn)色至采收期間日照時數(shù)呈極顯著的線性關(guān)系, 與最低氣溫、≥10 ℃積溫均呈二次曲線關(guān)系。采收前30 d的最低氣溫越高, 表明采收越早, 成熟度越差, 含糖量越低。采收前30 d降水量平均每增多10 mm, 總糖降低9.7 g×L-1, 與劉玉蘭等[16]、修德仁等[36]得出采收前1個月降水量每增加10 mm, 含糖量降低8.2 g×L-1的結(jié)果相近。

    參考國際植物遺傳委員會(BPGR)葡萄品質(zhì)分級標(biāo)準(zhǔn)[32-35], 結(jié)合Smith[37]對‘赤霞珠’等4個品種的分級指標(biāo)和李記明等[5-6]給出的國際上釀造優(yōu)質(zhì)、特優(yōu)級葡萄酒對總糖的要求, 將‘赤霞珠’總糖分為特優(yōu)(≥210 g×L-1)、優(yōu)(195.1~210 g×L-1)、良(170.1~195 g×L-1)、中(155.1~170 g×L-1)和差(≤155 g×L-1)5個等級, 根據(jù)總糖與氣象因子的關(guān)系反推不同氣象因子的閾值范圍(表3)??梢? 釀造優(yōu)質(zhì)葡萄酒要求開花至采收期間平均氣溫在22.6 ℃以下, ≥10 ℃積溫≥2 189 ℃·d, 日照≥867 h。其中轉(zhuǎn)色至采收平均氣溫、平均最低氣溫、≥10 ℃積溫和日照分別為≤21.7 ℃、≤16.6 ℃、≥789.6 ℃·d 和≥303 h, 采收前30 d最低氣溫、降水量需≤16.1 ℃和≤49.1 mm。生產(chǎn)上可根據(jù)表3的方程和指標(biāo)估算釀酒葡萄總糖含量, 達(dá)到優(yōu)質(zhì)級別后開始采摘。

    表3 ‘赤霞珠’葡萄總糖含量(y)與不同生育期氣象因子(x)的關(guān)系及不同總糖等級的氣象因子范圍

    =51,0.05=4.038,0.01=7.20,20.005=0.151,20.001=0.201.

    2.2.2 總酸含量與氣象因子的關(guān)系

    葡萄總酸含量與開花至采收期間的平均氣溫呈二次曲線關(guān)系, 氣溫高導(dǎo)致酸分解, 總酸含量降低??偹崤c轉(zhuǎn)色至采收期間的日均降水量呈線性正相關(guān), 轉(zhuǎn)色后降水越多, 總酸越多??偹崤c開花至采收、轉(zhuǎn)色至采收累積日照時數(shù)均呈對數(shù)關(guān)系, 日照越多, 總酸含量越低。采收前10 d最高氣溫與總酸含量呈二次曲線關(guān)系, 氣溫越高, 表示采摘越早, 糖酸轉(zhuǎn)化不充分, 總酸含量越高(表4)。

    按照葡萄酒質(zhì)量要求, 將‘赤霞珠’總酸含量分為特優(yōu)(8.0~9.4 g×L-1)、優(yōu)(9.5~10.5 g×L-1)、良(10.6~12.0 g×L-1)、中(12.1~13.5 g×L-1)和差(>13.5 g×L-1)5個等級, 根據(jù)各因子的關(guān)系方程反推不同總酸水平的氣象因子閾值, 發(fā)現(xiàn)達(dá)‘赤霞珠’優(yōu)質(zhì)葡萄酒的總酸含量要求開花至采收期間平均氣溫在22.2 ℃以下, 累積日照930 h以上。轉(zhuǎn)色至采收期間≥10 ℃積溫946 ℃·d以上, 日均降水量1.4 mm以下, 累積日照357.5 h以上, 采收前10 d平均最高氣溫在25.7 ℃以下。

    表4 ‘赤霞珠’葡萄總酸含量(y)與氣象因子(x)的關(guān)系及不同總酸等級的氣象因子范圍

    =51,0.05=4.038,0.01=7.20,20.05=0.076,20.01=0.125 7,20.005=0.151,20.001=0.201.

    2.2.3 糖酸比與氣象因子的關(guān)系

    隨著開花至采收期間平均氣溫升高, 糖酸比下降, 平均氣溫每升高1 ℃, 糖酸比下降4.3。采收越遲, 開花至采收的平均氣溫越低, ≥10 ℃積溫和累積日照越多, 糖酸比越高?!?0 ℃積溫每增加100 ℃×d, 糖酸比升高1.2; 日照每增加100 h, 糖酸比升高2.8。轉(zhuǎn)色至采收期間的氣象因子與糖酸比的關(guān)系類似。采收前30 d的降水量與糖酸比存在負(fù)指數(shù)關(guān)系, 隨著降水增多, 糖酸比下降。采收前30 d平均最低氣溫與糖酸比負(fù)相關(guān), 低溫可促進(jìn)糖酸轉(zhuǎn)化, 過早采摘糖酸比偏低。

    王華等[38]用糖酸比衡量葡萄的成熟度, 分為特優(yōu)(27~30)、優(yōu)(24~26.9或30.1~33.0)、良(19.0~23.9或33.1~36.0)、中(14.0~18.9或36.1~40.0)和差(<14.0或>40.0)5個等級。根據(jù)表5的方程反推不同氣象因子的閾值, 可見優(yōu)質(zhì)葡萄酒的理想糖酸比需開花至采收期間平均氣溫19.1~21.2 ℃, ≥10 ℃積溫2 654 ℃·d以上, 日照1 068 h以上。轉(zhuǎn)色至采收期間平均氣溫15.9~19.4 ℃, ≥10 ℃積溫1 238.9~ 1 953.2 ℃·d, 日照501.1~813.6 h。采收前30 d降水量21.8 mm以下, 平均最低氣溫5.6~12.0 ℃。歷年樣本大部分處于優(yōu)質(zhì)指標(biāo)范圍內(nèi), 表明賀蘭山東麓‘赤霞珠’適合形成釀造優(yōu)質(zhì)葡萄酒所需的糖酸比。

    2.2.4 pH與氣象因子的關(guān)系

    釀酒葡萄的pH與開花、轉(zhuǎn)色至采收的積溫為正相關(guān), 與開花至采收期的日照時數(shù)、采收前30 d降水量均為負(fù)相關(guān), 與轉(zhuǎn)色至采收的日照時數(shù)呈二次曲線關(guān)系(表6)。pH隨采收前30 d降水量的增大而減小。采集樣本的pH為2.9~3.6, 與李記明等[6]給出的適宜釀造優(yōu)質(zhì)葡萄酒的pH基本相同, 表明賀蘭山東麓釀酒葡萄的pH適宜釀造優(yōu)質(zhì)酒, 氣象條件對pH的影響有限。

    優(yōu)質(zhì)釀酒葡萄pH為3.0~3.6, 按照0.1的間距分為特優(yōu)(3.2~3.4)、優(yōu)(3.41~3.5或3.1~3.19)、良(3.51~3.6或3.0~3.09)、中(3.61~3.7或2.9~2.99)和差(>3.7或<2.9)5個等級。由于化驗樣本的pH均在2.9~3.6, 大部分處于優(yōu)質(zhì)級別, 有些氣象因子無法反推分類等級。按照釀造優(yōu)質(zhì)葡萄酒的pH分類等級, 需要開花至采收期間≥10 ℃積溫1671~2 671 ℃·d, 日均日照在12.5 h×d-1以下, 轉(zhuǎn)色至采收期間≥10 ℃積溫500~1 360.8 ℃·d, 日照時數(shù)203.5~525.0 h, 且采收前30 d降水量少于79.4 mm。

    2.2.5 單寧含量與氣象因子的關(guān)系

    單寧含量與開花至采收和轉(zhuǎn)色至采收期間平均最低氣溫均呈二次曲線關(guān)系, 最適值分別為15.2 ℃和13.3 ℃, 轉(zhuǎn)色后的影響更大。轉(zhuǎn)色至采收日照每增加100 h, 單寧增2.6 mmol×kg-1。采收前30 d的平均最低氣溫、20 d的降水量均與單寧含量有關(guān), 前20 d的降水量增加, 單寧值減少, 當(dāng)降水量由30 mm增加到40 mm時, 單寧減少1.2 mmol×kg-1。采收前10 d的氣溫日較差在10~13 ℃時單寧含量較高, 過小、過大單寧含量均會下降(表7)。

    表5 ‘赤霞珠’葡萄糖酸比(y)與氣象因子(x)的關(guān)系及不同糖酸比等級的氣象因子范圍

    =51,0.05=4.038,0.01=7.20,20.05=0.076,20.01=0.125 7,20.005=0.150 5,20.001=0.200 7.

    表6 ‘赤霞珠’葡萄pH(y)與氣象因子(x)的關(guān)系及不同pH等級的氣象因子范圍

    =51,0.05=4.038,0.01=7.20,20.05=0.076,20.01=0.125 7,20.005=0.150 5,20.001=0.200 7.

    按照葡萄酒釀造質(zhì)量對單寧的要求, 分為特優(yōu)(>30.1 mmol×kg-1)、優(yōu)(25.1~30 mmol×kg-1)、良(20.1~25 mmol×kg-1)、中(15.1~20 mmol×kg-1)和差(≤15.0 mmol×kg-1)5個等級。釀造優(yōu)質(zhì)葡萄酒需開花、轉(zhuǎn)色至采收期間平均日最低氣溫分別為15.0~17.3 ℃和13.5~17.1 ℃, 轉(zhuǎn)色至采收日照430.6~633.1 h, 采收前30 d平均最低氣溫9.6~17.4 ℃, 采收前20 d降水量不超過59.5 mm, 采收前10 d平均氣溫日較差在9.9 ℃以上。

    2.3 ‘赤霞珠’品質(zhì)的氣象因子綜合模擬模型

    根據(jù)2.2中釀酒葡萄不同品質(zhì)要素與氣象因子的關(guān)系模型, 以顯著性強且相互獨立的因子按照決 定系數(shù)的權(quán)重分別構(gòu)建了5項品質(zhì)要素的綜合評價模型(表8), 復(fù)相關(guān)系數(shù)均有較大幅度提升, 各地可根據(jù)氣象實況推算‘赤霞珠’應(yīng)表現(xiàn)出的各項品質(zhì)。

    2.4 ‘赤霞珠’品質(zhì)等級評價綜合判別模型與指標(biāo)

    參考李記明[5]、張曉煜等[31]的研究, 結(jié)合釀酒師的建議, 按照各項品質(zhì)要素對葡萄酒品質(zhì)的貢獻(xiàn)加權(quán):

    p=0.54Sc+0.10Ac+0.21S/A+0.10Tan+ 0.05Ph(1)

    式中:p為‘赤霞珠’綜合氣象品質(zhì);ScAcS/ATan和pH分別為總糖、總酸、糖酸比、單寧和pH的品質(zhì)級別, 均分為特優(yōu)(5級)、優(yōu)(4級)、良(3級)、中(2級)和差(1級)5個級別, 分別按照表3至表6歸類, 得到‘赤霞珠’品質(zhì)綜合氣象等級判別指標(biāo)(表9)。

    表7 ‘赤霞珠’葡萄單寧含量(y)與氣象因子(x)的關(guān)系及不同單寧含量等級的氣象因子范圍

    =51,0.05=4.038,0.01=7.20,20.05=0.076,20.01=0.125 7,20.005=0.150 5,20.001=0.200 7.

    表8 ‘赤霞珠’葡萄各項品質(zhì)要素與氣象因子的最優(yōu)關(guān)系模型

    =51。fh、dfh: 開花—采收日照時數(shù)(h)和日均日照時數(shù)(h×d-1);fh: 開花—采收平均氣溫(℃); ATsch: 轉(zhuǎn)色—采收≥10 ℃積溫(℃·d);nch: 轉(zhuǎn)色—采收平均最低氣溫(℃);ch: 轉(zhuǎn)色—采收日照時數(shù)(h);n30: 采收前30 d平均最低氣溫(℃);30、20: 采收前30 d、20 d的降水量(mm);m10、d10: 采收前10 d平均最高氣溫和平均氣溫日較差(℃)。fhanddfhare sun duration and average daily sunshine hours from flower to harvest (h×d-1);fhis average temperature from flower to harvest (℃); ATschis accumulated temperature ≥10 ℃from pigmentation to harvest (℃·d);nchis average minimum temperature from pigmentation to harvest (℃);chis sun duration from pigmentation to picking (h);n30is mean min temperature 30 days before picking (℃);30and20are precipitation 30 days and 20 days before picking (mm);m10andd10are mean max temperature and daily temperature difference 10 days before picking (℃).

    表9 ‘赤霞珠’葡萄品質(zhì)綜合氣象等級判別指標(biāo)

    2.5 估算模型的效果檢驗

    用建模的51個樣本回代檢驗(圖1), 5項品質(zhì)要素氣象模型回代均通過了0.001的檢驗和檢驗, 其中總糖、總酸、糖酸比和單寧的模擬效果較好,在0.59以上, RMSE較小。pH模擬離散度相對較大, 與樣本的pH變化不大, 氣象因子的影響相對有限有關(guān)。

    圖1 ‘赤霞珠’葡萄5項品質(zhì)要素的氣象因子回代模擬效果

    為檢驗?zāi)P偷脑u估效果, 選用2003—2011年永寧、青銅峽和蘆花臺葡萄基地未參與建模的‘赤霞珠’化驗樣本50份, 按照最優(yōu)模型估算了品質(zhì)要素值, 樣本按照總糖由小到大排列(圖2)。結(jié)果表明, 總糖、總酸含量與實測結(jié)果接近, 變化趨勢一致。糖酸比除2005年、2011年的樣本估算誤差較大外, 其他年份基本準(zhǔn)確, 單寧個別樣本的估算誤差接近10 mmol×kg-1, 誤差相對較大, 但變化趨勢一致。pH估算值較穩(wěn)定, 沒有反映出樣本的差異。總的來看, ‘赤霞珠’各品質(zhì)要素的實測值與估算值同比增減, 變化趨勢較為一致。

    圖2 用其他未參與建模的‘赤霞珠’葡萄樣本對5項品質(zhì)要素的模擬檢驗

    按照‘赤霞珠’單項品質(zhì)分類閾值, 分別將估算的每項品質(zhì)要素按照5個等級歸類, 實測值和估算值均按照式(1)計算樣本的綜合得分值, 再根據(jù)表9進(jìn)行綜合品質(zhì)等級分類。未參與建模的50個樣本中, 有28個樣本與實況等級相同; 有11個樣本將3級估算為4級; 有7個樣本將5級估算為4級; 僅有4個樣本將3級估算為5級(圖3)。

    3 討論

    3.1 光熱條件主要影響釀酒葡萄的成熟度, 模型可用于預(yù)測適宜采收期

    釀酒葡萄的品質(zhì)與葡萄成熟度有關(guān)[8], 充分成熟的葡萄含糖量高, 糖酸轉(zhuǎn)化充分, pH適度, 葡萄籽和果皮內(nèi)芳香類物質(zhì)豐富, 單寧較高。從模擬結(jié)果來看, 開花至采收的日數(shù)與總糖、糖酸比、單寧均為極顯著正相關(guān), 與總酸極顯著負(fù)相關(guān)??偺?、糖酸比和單寧均隨著采收的推遲而增加, 總酸則相反, 表明結(jié)實期越長, 葡萄成熟度越好, 品質(zhì)也越高。王華等[38]認(rèn)為溫度和光照是影響葡萄各種有效成分的主要因素, 無論是開花至采收、轉(zhuǎn)色至采收還是采收前一個月的積溫、累積日照與上述物質(zhì)的含量均顯著相關(guān), 反映了氣象條件與成熟度的關(guān)系, 也說明后期光熱資源不足影響葡萄成熟度。由于不同地域的氣象條件差異很大, 且高溫、低溫年果實的生長日數(shù)差距很大, 種植戶難以隨時測定諸多品質(zhì)要素, 往往憑經(jīng)驗決定采收期。采用積溫、累積日照等氣象因子估算各類物質(zhì)的含量能客觀反映葡萄的成熟度, 可為生產(chǎn)上確定適宜采收期提供客觀依據(jù)。

    3.2 著色后期的氣象條件對釀酒葡萄品質(zhì)的影響更大

    光熱條件主要影響葡萄總糖、總酸的積累和轉(zhuǎn)化, 進(jìn)而影響糖酸比和pH?!嘞贾椤偺?、糖酸比與轉(zhuǎn)色至采收期間及采收前30 d達(dá)到極顯著的氣象因子多于開花至采收期, 表明后期氣象條件的影響更大, 也說明后期光熱不足和水分過多更影響總糖積累和轉(zhuǎn)化??偹崤c開花至采收期間氣象因子的相關(guān)系數(shù)比轉(zhuǎn)色后的大, 說明總酸的積累與結(jié)實前中期的氣象條件關(guān)系密切。pH與轉(zhuǎn)色后的日照及采收前30 d的降水量通過檢驗, 表明轉(zhuǎn)色后熱量不足、陰雨天氣影響pH。單寧與轉(zhuǎn)色后通過顯著性檢驗的氣象因子多, 表明轉(zhuǎn)色后溫度下降、降水偏多影響單寧含量。采收越早, 單寧越少, 后期光熱越充足, 降水越少, 單寧含量越高。這與鄧浩亮等[39]、李濤等[40]得出的果實著色至成熟期適度水分脅迫可增加花青苷、還原糖、單寧、總酚含量, 且有效抑制可滴定酸積累的結(jié)果相符。

    3.3 釀酒葡萄品質(zhì)氣象等級劃分參考了釀造標(biāo)準(zhǔn), 檢驗效果較好

    品質(zhì)氣象等級的劃分參考了國外葡萄質(zhì)量等級劃分標(biāo)準(zhǔn)[41]、山葡萄劃分的5級標(biāo)準(zhǔn)及4項頒布標(biāo)準(zhǔn), 依據(jù)葡萄釀造品質(zhì)對總糖、總酸、糖酸比、pH和單寧含量的要求劃分為5級, 并按此等級確定了氣象因子的閾值和品質(zhì)氣象等級。因此, 利用氣象模型估算的結(jié)果可與化驗結(jié)果相對應(yīng), 能客觀反映氣象條件應(yīng)形成的品級, 而不是根據(jù)品質(zhì)波動區(qū)間劃分的相對優(yōu)劣等級, 對評價和預(yù)估葡萄酒質(zhì)量有一定幫助。

    圖3 未參試的‘赤霞珠’葡萄樣本的綜合氣象品質(zhì)評分估算檢驗

    樣本回代效果表明, 總糖、總酸、糖酸比和單寧的回代效果準(zhǔn)確, pH與氣象因子的關(guān)系較弱, 說明賀蘭山東麓的氣象條件不是限制釀酒葡萄pH變化的主要因素。用未參與建模的樣本檢驗, 總糖、總酸含量與實測結(jié)果接近, 糖酸比基本準(zhǔn)確, 單寧誤差相對較大, pH沒有反映出樣本的差異。綜合等級分值與實況相同的樣本占56%, 相差1個級別以內(nèi)的樣本占36%, 2個級別誤差為8%??傮w來看, 該方法和指標(biāo)能客觀反映‘赤霞珠’的品質(zhì)等級, 對綜合品質(zhì)等級的評定效果較好。

    3.4 模型的局限與不確定性

    葡萄酒的品質(zhì)除了原料品質(zhì)外, 還與釀造技術(shù)關(guān)系很大, 氣象條件只是影響釀酒葡萄品質(zhì)的一個重要方面, 化驗的品質(zhì)要素僅上述幾項, 特別是缺少芳香烴類等眾多品質(zhì)要素的化驗資料, 綜合氣象品質(zhì)還難以完全反映葡萄酒的色、香、味、柔度等口感。因此, 本方法可作為評價年份酒優(yōu)劣的客觀參考依據(jù), 不能將其與葡萄酒品質(zhì)完全劃等號。其次, 采用氣象因子模擬的先決條件是氣象因子在取樣樣本的正常范圍內(nèi)波動, 若氣象因子超出模擬范圍, 模擬的準(zhǔn)確性會受到影響。第三, 有些統(tǒng)計關(guān)系中采用二次曲線模擬上升段或下降段, 但實況氣象資料超過了模擬范圍, 如采收前30 d的降水量建模時最大才到60 mm, 歷史上同期降水量超過100 mm, 應(yīng)用時要注意因子的值域范圍不應(yīng)超過建模樣本的值域區(qū)間, 否則可能得出錯誤的結(jié)論。另外, 采用的品質(zhì)化驗資料僅限于寧夏范圍, 結(jié)果有一定的局限性, 應(yīng)用到我國其他產(chǎn)區(qū)時要謹(jǐn)慎, 要做本地化驗證。今后可參照本方法建立基于全國各地采樣化驗資料的品質(zhì)要素及氣象品質(zhì)模型, 也可根據(jù)釀造要求增加其他必要的品質(zhì)成分指標(biāo), 加大氣候差異, 擴大模型的適用性。還可以仿照此方法建立其他葡萄品種的品質(zhì)氣象評價模型, 形成釀酒葡萄氣候品質(zhì)評價與鑒定的系列化、客觀化和定量化。

    4 結(jié)論

    葡萄酒的品質(zhì)主要取決于優(yōu)質(zhì)釀酒原料, 與品種、土壤和氣象條件均有關(guān)系, 在正常水肥管理條件下, 氣象條件是造成品質(zhì)年際間波動的主因。本文研究了‘赤霞珠’5項品質(zhì)要素與氣象因子的關(guān)系, 構(gòu)建了最優(yōu)單項品質(zhì)要素的氣象評價模型, 選擇生物學(xué)意義明顯的因子, 分別構(gòu)建了5項品質(zhì)要素的綜合加權(quán)模型和氣象因子等級閾值; 借鑒前人研究、標(biāo)準(zhǔn)和釀酒師的意見, 構(gòu)建了‘赤霞珠’綜合品質(zhì)等級評價權(quán)重模型和等級分類標(biāo)準(zhǔn), 用非參與建模的樣本檢驗, 顯示模型較準(zhǔn)確地反映了賀蘭山東麓‘赤霞珠’品質(zhì)的綜合等級。其方法、等級分類指標(biāo)可為年份酒評級提供參考依據(jù), 有助于打造我國優(yōu)質(zhì)釀酒葡萄產(chǎn)地, 實施高端品牌戰(zhàn)略, 提升品牌知名度。

    致謝: 立蘭酒莊首席釀酒師邵青松先生對本文提出了寶貴意見, 在此表示感謝!

    [1] 李玉鼎, 李欣, 張光弟, 等. 寧夏葡萄酒產(chǎn)業(yè)持續(xù)發(fā)展的思考與技術(shù)創(chuàng)新[J]. 中外葡萄與葡萄酒, 2010, (3): 74-77 LI Y D, LI X, ZHANG G D, et al. Thoughts and technology innovation on sustainable development of wine industry in Ningxia[J]. Sino-overseas Grapevine & Wine, 2010, (3): 74-77

    [2] 李玉鼎, 張軍翔, 王戰(zhàn)斗, 等. 寧夏賀蘭山東麓葡萄年份酒與氣候[J]. 中外葡萄與葡萄酒, 2004, (2): 54-57 LI Y D, ZHANG J X, WANG Z D, et al. Vintage wine and weather in the area of Ningxia at the eastern foot of Helan Mountain[J]. Sino-overseas Grapevine & Wine, 2004, (2): 54–57

    [3] 張軍翔, 李玉鼎, 王戰(zhàn)斗, 等. 氣象因子對葡萄酒質(zhì)量影響的研究[J]. 山西果樹, 2004, (2): 3-5 ZHANG J X, LI Y D, WANG Z D, et al. The study of influence of climatic factor on the quality of wine[J]. Shanxi Fruits, 2004, (2): 3–5

    [4] 趙清華. 年份氣候和采摘期對瓊瑤漿葡萄酒品質(zhì)的影響[J]. 山西食品工業(yè), 2004, (3): 23-24 ZHAO Q H. Influence of year climate and picking period on the quality of the wine gewurztraminer[J]. Shanxi Food Industry, 2004, (3): 23–24

    [5] 李記明. 關(guān)于葡萄品質(zhì)的評價指標(biāo)[J]. 中外葡萄與葡萄酒, 1999, (1): 54-57 LI J M. Evaluation index for grape quality[J]. Sino-overseas Grapevine & Wine, 1999, (1): 54-57

    [6] 李記明, 李華. 不同地區(qū)釀酒葡萄成熟度與葡萄酒質(zhì)量的研究[J]. 西北農(nóng)業(yè)學(xué)報, 1996, 5(4): 71–74 LI J M, LI H. Studies on winegrape maturity and wine quality in different ecological zones[J]. Acta Agriculturae Boreali-occidentalis Sinica, 1996, 5(4): 71–74

    [7] 王華, 李記明, 陶士衡, 等. 有效積溫與降水量的比值(T/P)與釀酒葡萄成熟度關(guān)系的研究[J]. 葡萄栽培與釀酒, 1994, (2): 8-9 WANG H, LI J M, TAO S H, et al. Study on the relationship between wine grape maturity and the ratio of effective accumulated temperature and precipitation (T/P)[J]. Viticulture & Enology, 1994, (2): 8-9

    [8] VR?I? S, ?U?TAR V, PULKO B, et al. Trends in climate parameters affecting winegrape ripening in northeastern Slovenia[J]. Climate Research, 2014, 58(3): 257-266

    [9] NEUMANN P A, MATZARAKIS A. Potential climate change impacts on winegrape must density and titratable acidity in Southwest Germany[J]. Climate Research, 2014, 59(2): 161–172

    [10] 毛如志, 張國濤, 邵建輝, 等. 低海拔和高海拔產(chǎn)區(qū)氣象因子對‘美樂’葡萄漿果品質(zhì)和代謝組的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2016, 24(4): 506–516 MAO R Z, ZHANG G T, SHAO J H, et al. Response of ‘Merlot’ grape berry quality and metabolome to meteorological factors at both low and high altitudes[J]. Chinese Journal of Eco-Agriculture, 2016, 24(4): 506–516

    [11] WEBB L B, WHETTON P H, BARLOW E W R. Climate change and winegrape quality in Australia[J]. Climate Research, 2008, 36(2): 99-111

    [12] WEBB L B, WHETTON P H, BARLOW E W R. Modelling the relationship between climate, winegrape price and winegrape quality in Australia[J]. Climate Research, 2008, 36(2): 89–98

    [13] 李華. 葡萄酒品嘗學(xué)[M]. 北京: 中國青年出版社, 1992 LI H. Wine Tasting[M]. Beijing: China Youth Publishing House, 1992

    [14] 翟衡, 杜金華管雪強, 等. 釀酒葡萄栽培及加工技術(shù)[M]. 北京: 中國農(nóng)業(yè)出版社, 2001: 199-238 ZHAI H, DU J H, GUAN X Q, et al. Wine Grape Cultivation and Processing Technology[M]. Beijing: China Agricultural Publishing House, 2001: 199–238

    [15] 劉玉蘭. 釀酒葡萄的品質(zhì)指標(biāo)分析及其與氣象條件關(guān)系的試驗研究[D]. 南京: 南京信息工程大學(xué), 2006: 1-43 LIU Y L. The analysis of quality indexes of wine grape and the experimental study on meteorological condition’s impact on some qualities[D]. Nanjing: Nanjing University of Information Science & Technology, 2006: 1-43

    [16] 劉玉蘭, 鄭有飛, 張曉煜. 氣象條件對釀酒葡萄品質(zhì)影響的研究進(jìn)展[J]. 中外葡萄與葡萄酒, 2008, (1): 28-29 LIU Y L, ZHENG Y F, ZHANG X Y. Research progress of effect of meteorological condition on wine grape quality[J]. Sino-overseas Grapevine & Wine, 2008, (1): 28-29

    [17] 陳曉前, 張鐵. 試述pH值與葡萄酒的關(guān)系[J]. 釀酒, 2002, 29(2): 62-63 CHEN X Q, ZHANG T. Relationship between pH value and the wine [J]. Liquor Making, 2002, 29(2): 62–63

    [18] 張曉煜, 劉玉蘭, 張磊, 等. 氣象條件對釀酒葡萄若干品質(zhì)因子的影響[J]. 中國農(nóng)業(yè)氣象, 2007, 28(3): 326-330 ZHANG X Y, LIU Y L, ZHANG L, et al. Influence of meteorological conditions on some quality factors of wine grape[J]. Chinese Journal of Agrometeorology, 2007, 28(3): 326-330

    [19] 岳俊波, 侯勇. 淺談單寧在紅葡萄酒中的作用[J]. 釀酒, 2000, 27(6): 51 YUE J B, HOU Y. Discussion on the function of tannin in red wine[J]. Liquor-Making, 2000, 27(6): 51

    [20] 周濤, 張富國, 白國勝, 等. 風(fēng)沙土土壤的磷素狀況及施磷對釀酒葡萄品質(zhì)的影響[J]. 中國農(nóng)業(yè)科學(xué), 2002, 35(2): 169–173 ZHOU T, ZHANG F G, BAI G S, et al. Phosphorus status of sand soil and effects of phosphorus fertilization on quality of wine-grapes[J]. Scientia Agricultura Sinica, 2002, 35(2): 169–173

    [21] 趙磊, 吳玉文, 段長青, 等. 寧夏賀蘭山地區(qū)與河北沙城地區(qū)赤霞珠和美樂葡萄酒香氣物質(zhì)的比較[J]. 中外葡萄與葡萄酒, 2009, (11): 15-19 ZHAO L, WU Y W, DUAN C Q, et al. Comparison of aromatic compounds of Cabernet Sauvignon wine and Merlot wine originated from Helan Mountain of Ningxia and Shacheng of Hebei[J]. Sino-overseas Grapevine & Wine, 2009, (11): 15–19

    [22] 王秀琴, 陳小波, 戰(zhàn)吉成, 等. 生態(tài)因素對釀酒葡萄和葡萄酒品質(zhì)的影響[J]. 食品科學(xué), 2006, 27(12): 791–797 WANG X Q, CHEN X B, ZHAN J C, et al. Effects of ecological factors on quality of winegrape and wine[J]. Food Science, 2006, 27(12): 791–797

    [23] PIRCALABU L A, TUDORACHE A, SERDINESCU A, et al. The quantitatively and qualitatively potential of the wine grape harvest in relation with the wine climate milessima[J]. Bulletin of University of Agricultural Sciences, 2015, 72(1): 165–169

    [24] JONES G V, REID R, VILKS A. Climate, grapes, and wine: Structure and suitability in a variable and changing climate[M]//DOUGHERTY P H. The Geography of Wine: Regions, Terroir and Techniques. Dordrecht: Springer, 2012: 109–133

    [25] CHENG G, HE Y N, YUE T X, et al. Effects of climatic conditions and soil properties on Cabernet Sauvignon berry growth and anthocyanin profiles[J]. Molecules, 2014, 19(9): 13683–13703

    [26] 張娟, 王曉宇, 田呈瑞, 等. 基于酚類物質(zhì)的釀酒紅葡萄品種特性分析[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(7): 1370–1382 ZHANG J, WANG X Y, TIAN C R, et al. Analysis of phenolic compounds in red grape varieties[J]. Scientia Agricultura Sinica, 2015, 48(7): 1370–1382

    [27] 周淑珍, 張軍翔, 曹景麗. 不同氧化程度和年份葡萄酒總酚含量研究[J]. 中外葡萄與葡萄酒, 2008, (1): 49–51 ZHOU S Z, ZHANG J X, CAO J L. Study on wine total phenol content for different oxidation degree and age[J]. Sino-overseas Grapevine & Wine, 2008, (1): 49–51

    [28] 楊曉帆, 高媛, 韓梅梅, 等. 云南高原區(qū)釀酒葡萄果實香氣物質(zhì)的積累規(guī)律[J]. 中國農(nóng)業(yè)科學(xué), 2014, 47(12): 2405–2416 YANG X F, GAO Y, HAN M M, et al. Accumulation characteristics of volatile compounds in wine grape berries grown in high altitude regions of Yunnan[J]. Scientia Agricultura Sinica, 2014, 47(12): 2405–2416

    [29] 房玉林, 孫偉, 萬力, 等. 調(diào)虧灌溉對釀酒葡萄生長及果實品質(zhì)的影響[J]. 中國農(nóng)業(yè)科學(xué), 2013, 46(13): 2730–2738 FANG Y L, SUN W, WAN L, et al. Effects of regulated deficit irrigation (RDI) on wine grape growth and fruit quality[J]. Scientia Agricultura Sinica, 2013, 46(13): 2730–2738

    [30] Gouveia C M, Gani é A, Liberato M L R. Climate variability and wine quality over Portuguese regions[C]. Geophysical Research Abstracts, 2015, 17: 12361

    [31] 張曉煜, 亢艷莉, 袁海燕, 等. 釀酒葡萄品質(zhì)評價及其對氣象條件的響應(yīng)[J]. 生態(tài)學(xué)報, 2007, 27(2): 740–745 ZHANG X Y, KANG Y L, YUAN H Y, et al. The quality evaluation of wine grape and its respond to weather condition[J]. Acta Ecologica Sinca, 2007, 27(2): 740–745

    [32] 中國輕工業(yè)聯(lián)合會. GB 15037—2006 葡萄酒[S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2008 China Light Industry Federation. GB 15037—2006 Wine[S]. Beijing: China Standards Press, 2008

    [33] 全國原產(chǎn)地域產(chǎn)品標(biāo)準(zhǔn)化工作組. GB/T 19504—2008 地理標(biāo)志產(chǎn)品賀蘭山東麓葡萄酒[S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2008 National Work Group on Standardization of Geographical Origin Products. GB/T 19504—2008 Product of Geographical Indication — Wine in Helan Mountain East Region[S]. Beijing: China Standards Press, 2008

    [34] 全國原產(chǎn)地域產(chǎn)品標(biāo)準(zhǔn)化工作組. GB/T 19049—2008 地理標(biāo)志產(chǎn)品昌黎葡萄酒[S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2009 National Work Group on Standardization of Geographical Origin Products. GB/T 19049—2008 Product of Geographical Indication — Changli Wine[S]. Beijing: China Standards Press, 2009

    [35] 王同坤, 項殿芳, 齊永順, 等. DB13/T 912—2007 釀酒葡萄質(zhì)量標(biāo)準(zhǔn)[S]. 石家莊: 河北省技術(shù)監(jiān)督局, 2007 WANG T K, XIANG D F, QI Y S, et al. DB13/T 912—2007 Quality Standard for Wine Grape[S]. Shijiazhuang: Hebei Provincial Administration of Quality and Technical Supervision, 2007

    [36] 修德仁, 周潤生, 晁無疾, 等. 干紅葡萄酒用品種氣候區(qū)域化指標(biāo)分析及基地選擇[J]. 葡萄栽培與釀酒, 1997, (3): 22–26 XIU D R, ZHOU R S, CHAO W J, et al. Index analysis of climate regionalization for dry red wine varieties and the varieties base selection[J]. Viticulture & Enology, 1997, (3): 22–26

    [37] SMITH P. Assessment of relationships between grape chemical composition and grape allocation grade for Cabernet Sauvignon, Shiraz and Chardonnay[J]. Australian and New Zealand Grapegrower and Winemaker, 2015, (620): 30–32

    [38] 王華, 李華, 李記明, 等. 釀酒葡萄成熟度與主要氣象因子關(guān)系研究[J]. 西北農(nóng)業(yè)學(xué)報, 1993, 2(4): 71–74 WANG H, LI H, LI J M, et al. Studies on relationship between the main climatologic factors and the maturity of grapeyine[J]. Acta Agriculturae Boreali-occidentalis Sinica, 1993, 2(4): 71–74

    [39] 鄧浩亮, 孔維萍, 張恒嘉, 等. 不同生育期調(diào)虧灌溉對釀酒葡萄耗水及果實品質(zhì)的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2016, 24(9): 1196–1205 DENG H L, KONG W P, ZHANG H J, et al. Effect of regulated deficit irrigation at different growth stages on water consumption and fruit quality of wine grape[J]. Chinese Journal of Eco-Agriculture, 2016, 24(9): 1196–1205

    [40] 李濤, 張建豐, 張江輝, 等. 以灌溉定額為參數(shù)的葡萄果實含糖量和果型指數(shù)的數(shù)學(xué)模型[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2010, 18(2): 348–351 LI T, ZHANG J F, ZHANG J H, et al. Fuzzy model for grape sugar content and fruit shape index with irrigation quota as parameter[J]. Chinese Journal of Eco-Agriculture, 2010, 18(2): 348–351

    [41] 劉裕嚴(yán). 評價葡萄品質(zhì)的新方法[J]. 葡萄栽培與釀酒, 1985, (3): 39–41 LIU Y Y. A new method for evaluating grape quality[J]. Viticulture & Enology, 1985, (3): 39–41

    Meteorological conditions and rating method of quality formation of ‘Cabernet Sauvignon’ grape in eastern foothills of Helan Mountain*

    MA Liwen, LI Jianping**, HAN Yingjuan, LI Wanchun

    (Key Laboratory for Meteorological Disaster Monitoring, Early Warning and Risk Management of Characteristic Agriculture in Arid Regions, China Meteorological Administration / Ningxia Key Lab of Meteorological Disaster Prevention and Reduction, Yinchuan 750002, China)

    The quality of grape wine is largely determined by the quality of grapes, which is closely related to meteorological conditions during grape growing season. Under normal management conditions and agricultural technology, the method of evaluation of meteorological conditions for fruit growth is called quality meteorological evaluation of wine grape. Meteorological condition rating techniques and meteorological indexes research on wine grape quality formation provide the base lines of climatic conditions for rating commercial wine, vintage identification and hoarding. In this paper, five indexes were tested for total sugar, total acid, sugar-acid ratio, pH and tannin content of 101 ‘Cabernet Sauvignon’ fruit samples collected at various sites in eastern Helan Mountain during 2003–2011. Test data for 51 samples were selected to analyze for correlations with meteorological data for fruit growth. Based on the analysis of biological effects of meteorological factors on quality, meteorological factors affecting each quality index were determined for significant level. Based on the determination coefficient of each meteorological factor to a quality index, five weight models of quality indexes with significant meteorological factors were constructed. Thus the quality indexes of ‘Cabernet Sauvignon’ grape were reckoned with the meteorological conditions for fruit growth period. Using determination coefficients of meteorological factors as weight, an integrated meteorological simulation model with 5 quality indicators was constructed and a comprehensive quality meteorological evaluation model was established by using the weight of individual quality indicators. Using related research results coupled with the winemaker suggestion, a comprehensive meteorological scoring model for wine grape quality was constructed based on the contribution of each quality index and wine quality. In relation to previous studies and the criteria of quality of grape and wine, 5 quality indicators and quality grades were classified. The threshold meteorological factors corresponding to the wine quality and classification criteria for grape climatic conditions of ‘Cabernet Sauvignon’ were also determined. It implied that based on the requirements of total sugar, total acid, sugar-acid ratio, pH and tannin content needed for brewing high quality wine, the 5 quality indicators and comprehensive meteorological grades of ‘Cabernet Sauvignon’ were divided into five grades. This included excellent (grade 5), very good (grade 4), good (grade 3), medium (grade 2) and poor (grade 1). Then based on the threshold of each quality grade and the comprehensive meteorological simulation equation of the above 5 quality indicators, the threshold of the meteorological factors affecting each quality grade was deduced. This was in turn used as the meteorological factor classification index for evaluating the climate condition for quality of ‘Cabernet Sauvignon’ grape. The effects of simulation and grape classification of the 5 quality indicators and ‘Cabernet Sauvignon’ grape in the foothills of Helan Mountain were tested using the 51 samples selected to establish the model. The results showed that all the models of the 5 quality indicators passed the 0.001-test and-test, indicating that the estimated values were within the measured values and that the variation trends were consistent. Among these, the simulation effects of total sugar content, total acid content, sugar-acid ratio and tannin content were better (≥ 0.59 and small RMSE) than the effects of pH (where there was a relative large simulation dispersion). To predict the effects of the evaluation model and indexes by using the other 50 samples that were not included in the model development and meteorological data during fruit development, the effects of the grading were tested for ‘Cabernet Sauvignon’ grape quality in the eastern foothills of Helan Mountain. The results indicated that the contents of total sugar and total acid were close to the measured values. In addition to a relative large error for few samples of sugar-acid ratio and tannin content, the error for the other samples was relatively small and the overall trend was consistent with the measured values. pH was relatively stable and all the samples were within the appropriate range for brewing high-quality wine, and differences in yearly meteorological conditions were difficult to determine, indicating that the local climatic conditions were conducive for the pH needed for high quality wine brewing. From the comprehensive score and grade, 28 samples were the same as the actual level and only 4 samples differed at 2 levels. The rest of the samples were within level 1, which more accurately reflected the quality of raw materials. The meteorological evaluation indexes and models of ‘Cabernet Sauvignon’ provided a feasible tool for evaluating the quality of wine grape in the eastern foothills of Helan Mountain.

    Eastern foothill of Helan Mountain; Cabernet Sauvignon; Wine grape; Meteorological quality; Meteorological factor; Quality grade; Evaluation model

    , E-mail: lijp_111@163.com

    10.13930/j.cnki.cjea.170579

    P49

    A

    1671-3990(2018)03-0453-14

    李劍萍, 主要研究方向為農(nóng)業(yè)氣候資源與生態(tài)。E-mail: lijp_111@163.com 馬力文, 主要從事特色作物農(nóng)業(yè)氣象、農(nóng)作物氣候品質(zhì)認(rèn)證等方面的研究。E-mail: ma_liwen@163.com

    2017-06-24

    2017-10-26

    * This study was supported by the Science and Technology Support Program of Ningxia Hui Autonomous Region (2013ZYN245), the Key Technology Integration Project of China Meteorological Administration (CMATG2007M44) and the Natural Science Foundation of Ningxia (ZN16202).

    * 寧夏回族自治區(qū)科技支撐計劃項目(2013ZYN245)、中國氣象局關(guān)鍵技術(shù)集成項目(CMATG2007M44)和寧夏自然科學(xué)基金項目 (ZN16202)資助

    Jun. 24, 2017; accepted Oct. 26, 2017

    馬力文, 李劍萍, 韓穎娟, 李萬春. 賀蘭山東麓‘赤霞珠’品質(zhì)形成氣象條件與評級方法研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2018, 26(3): 453-466

    MAL W, LIJ P, HANY J, LIW C. Meteorological conditions and ratingmethod of quality formation of ‘Cabernet Sauvignon’ grape in eastern foothills of Helan Mountain[J]. Chinese Journal of Eco-Agriculture, 2018, 26(3): 453-466

    中文字幕精品免费在线观看视频| 亚洲国产中文字幕在线视频| 在线av久久热| 熟女少妇亚洲综合色aaa.| 操美女的视频在线观看| 精品久久蜜臀av无| 在线播放国产精品三级| 少妇猛男粗大的猛烈进出视频| 欧美日韩国产mv在线观看视频| 成年人午夜在线观看视频| 国产亚洲av高清不卡| 欧美亚洲日本最大视频资源| 在线观看舔阴道视频| 国产精品国产高清国产av | 一级片免费观看大全| 亚洲精品自拍成人| 中文字幕人妻熟女乱码| 欧美一级毛片孕妇| 国产99久久九九免费精品| 日韩免费高清中文字幕av| 麻豆成人av在线观看| 激情视频va一区二区三区| 国产亚洲精品一区二区www | 不卡一级毛片| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 丝袜喷水一区| 国产在线免费精品| 日韩大码丰满熟妇| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区免费欧美| 午夜91福利影院| 水蜜桃什么品种好| 在线观看人妻少妇| 高清视频免费观看一区二区| av超薄肉色丝袜交足视频| 午夜精品久久久久久毛片777| 成人国产av品久久久| 热99国产精品久久久久久7| 欧美日韩av久久| 老司机靠b影院| 搡老岳熟女国产| 日韩欧美三级三区| 欧美精品一区二区大全| 国产精品美女特级片免费视频播放器 | 精品国产国语对白av| 午夜激情av网站| 久久久精品区二区三区| 成年动漫av网址| 美女高潮喷水抽搐中文字幕| 俄罗斯特黄特色一大片| 欧美精品一区二区免费开放| 国产精品 欧美亚洲| 美女福利国产在线| 国产无遮挡羞羞视频在线观看| 黄网站色视频无遮挡免费观看| 国产男靠女视频免费网站| 国产免费av片在线观看野外av| 国产三级黄色录像| 黄片大片在线免费观看| 亚洲免费av在线视频| 伦理电影免费视频| 757午夜福利合集在线观看| 超碰成人久久| 国产精品99久久99久久久不卡| 亚洲精品在线观看二区| 久久久久精品人妻al黑| 男女免费视频国产| 国产麻豆69| 黄色丝袜av网址大全| 99久久99久久久精品蜜桃| 侵犯人妻中文字幕一二三四区| 大型黄色视频在线免费观看| 美女午夜性视频免费| 五月开心婷婷网| 免费女性裸体啪啪无遮挡网站| 男人操女人黄网站| 无人区码免费观看不卡 | 少妇 在线观看| 久久久久久久大尺度免费视频| 久久国产精品男人的天堂亚洲| 欧美日本中文国产一区发布| 色精品久久人妻99蜜桃| 操美女的视频在线观看| 我的亚洲天堂| 亚洲精品一卡2卡三卡4卡5卡| 久久人人97超碰香蕉20202| 男女无遮挡免费网站观看| 亚洲中文av在线| 夜夜爽天天搞| 桃花免费在线播放| 色视频在线一区二区三区| 欧美精品啪啪一区二区三区| 18禁观看日本| a级片在线免费高清观看视频| videosex国产| 99国产精品99久久久久| 日本黄色视频三级网站网址 | www.自偷自拍.com| 久久久精品免费免费高清| 岛国毛片在线播放| 91国产中文字幕| 久久久久国产一级毛片高清牌| 一级毛片女人18水好多| 老熟妇乱子伦视频在线观看| 乱人伦中国视频| 久久久欧美国产精品| 黑丝袜美女国产一区| 啪啪无遮挡十八禁网站| 99久久精品国产亚洲精品| 亚洲欧洲精品一区二区精品久久久| 欧美乱码精品一区二区三区| 久久久久久人人人人人| www.999成人在线观看| cao死你这个sao货| 国产精品美女特级片免费视频播放器 | 超色免费av| av网站免费在线观看视频| 男女边摸边吃奶| 99香蕉大伊视频| 亚洲第一欧美日韩一区二区三区 | 激情在线观看视频在线高清 | 欧美在线黄色| 成人亚洲精品一区在线观看| 欧美日韩一级在线毛片| 亚洲 欧美一区二区三区| 一边摸一边抽搐一进一出视频| 好男人电影高清在线观看| 亚洲,欧美精品.| 久久人妻熟女aⅴ| 国产麻豆69| 99精品在免费线老司机午夜| 少妇猛男粗大的猛烈进出视频| 精品少妇内射三级| 久久精品国产综合久久久| 久久精品国产综合久久久| 老汉色∧v一级毛片| 自线自在国产av| 91精品三级在线观看| 天堂中文最新版在线下载| 一进一出抽搐动态| 2018国产大陆天天弄谢| 国产老妇伦熟女老妇高清| 精品亚洲乱码少妇综合久久| 人妻一区二区av| 亚洲精品久久成人aⅴ小说| 国产av一区二区精品久久| 天天添夜夜摸| 免费观看a级毛片全部| 国产片内射在线| 亚洲自偷自拍图片 自拍| 自线自在国产av| 99国产精品免费福利视频| 国产福利在线免费观看视频| 丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 18禁美女被吸乳视频| 99国产精品一区二区蜜桃av | 国产精品 欧美亚洲| 欧美精品av麻豆av| 国产亚洲午夜精品一区二区久久| 老司机福利观看| 久久久精品94久久精品| 老司机福利观看| 国产精品 欧美亚洲| 色视频在线一区二区三区| 一级毛片精品| 王馨瑶露胸无遮挡在线观看| 久久 成人 亚洲| 窝窝影院91人妻| 亚洲自偷自拍图片 自拍| 老熟女久久久| 天天添夜夜摸| 狂野欧美激情性xxxx| 黄色怎么调成土黄色| 黑人操中国人逼视频| 在线十欧美十亚洲十日本专区| 精品一区二区三区视频在线观看免费 | 极品少妇高潮喷水抽搐| 午夜福利视频在线观看免费| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区三| 精品乱码久久久久久99久播| 他把我摸到了高潮在线观看 | 在线看a的网站| 男女之事视频高清在线观看| 亚洲精品自拍成人| 欧美精品一区二区免费开放| 中文字幕另类日韩欧美亚洲嫩草| 久久久国产成人免费| 亚洲精品国产一区二区精华液| 一本大道久久a久久精品| 三上悠亚av全集在线观看| 精品乱码久久久久久99久播| 日本撒尿小便嘘嘘汇集6| 看免费av毛片| 亚洲专区中文字幕在线| 黑丝袜美女国产一区| 一区二区av电影网| 王馨瑶露胸无遮挡在线观看| 人人妻,人人澡人人爽秒播| 19禁男女啪啪无遮挡网站| 91国产中文字幕| 最新美女视频免费是黄的| 午夜激情av网站| 久久久欧美国产精品| 亚洲av第一区精品v没综合| 久久免费观看电影| 欧美中文综合在线视频| 国产精品1区2区在线观看. | 老司机深夜福利视频在线观看| 丝袜在线中文字幕| 成年人免费黄色播放视频| 叶爱在线成人免费视频播放| 精品久久久久久电影网| 中文字幕最新亚洲高清| 亚洲精品国产区一区二| 亚洲成av片中文字幕在线观看| 亚洲五月婷婷丁香| 中文字幕最新亚洲高清| 女人被躁到高潮嗷嗷叫费观| 亚洲av美国av| 人成视频在线观看免费观看| 中文字幕最新亚洲高清| 国产日韩欧美视频二区| 国产在线观看jvid| 免费看十八禁软件| 丁香六月欧美| 亚洲欧美色中文字幕在线| 午夜福利免费观看在线| 日韩中文字幕视频在线看片| 曰老女人黄片| 一进一出好大好爽视频| 麻豆国产av国片精品| 人人澡人人妻人| 国产亚洲一区二区精品| 亚洲久久久国产精品| 成人黄色视频免费在线看| 久久99热这里只频精品6学生| 妹子高潮喷水视频| 中文字幕人妻丝袜一区二区| 国产色视频综合| 亚洲精品乱久久久久久| 看免费av毛片| a级毛片黄视频| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 在线播放国产精品三级| 又紧又爽又黄一区二区| 日本精品一区二区三区蜜桃| 欧美成狂野欧美在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费电影在线观看| 久久久久精品人妻al黑| 极品人妻少妇av视频| netflix在线观看网站| 国产野战对白在线观看| 高清黄色对白视频在线免费看| 亚洲精品国产一区二区精华液| 国产精品一区二区精品视频观看| 午夜福利一区二区在线看| 久久免费观看电影| 久久99热这里只频精品6学生| 欧美+亚洲+日韩+国产| 成人亚洲精品一区在线观看| 久久精品国产99精品国产亚洲性色 | 嫁个100分男人电影在线观看| 高清在线国产一区| 久久婷婷成人综合色麻豆| 欧美激情极品国产一区二区三区| 国产成人一区二区三区免费视频网站| 色综合婷婷激情| av天堂在线播放| 日韩 欧美 亚洲 中文字幕| 十八禁人妻一区二区| 两个人看的免费小视频| 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| 黑丝袜美女国产一区| 18在线观看网站| 欧美日本中文国产一区发布| 啦啦啦中文免费视频观看日本| 国产精品免费视频内射| 91成人精品电影| 一二三四社区在线视频社区8| 少妇被粗大的猛进出69影院| 免费在线观看完整版高清| 亚洲七黄色美女视频| 国产成人一区二区三区免费视频网站| 怎么达到女性高潮| 黄网站色视频无遮挡免费观看| av在线播放免费不卡| 精品人妻熟女毛片av久久网站| 午夜福利免费观看在线| 一级片免费观看大全| 欧美日韩视频精品一区| 欧美激情极品国产一区二区三区| 国产福利在线免费观看视频| 极品人妻少妇av视频| 丰满迷人的少妇在线观看| 老司机亚洲免费影院| 日韩一区二区三区影片| 99久久人妻综合| 欧美午夜高清在线| 免费高清在线观看日韩| 汤姆久久久久久久影院中文字幕| 国产精品亚洲av一区麻豆| 色视频在线一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲av成人不卡在线观看播放网| 午夜福利一区二区在线看| 久久人妻熟女aⅴ| 国产熟女午夜一区二区三区| 男女高潮啪啪啪动态图| 国产成+人综合+亚洲专区| 色婷婷久久久亚洲欧美| 老司机午夜十八禁免费视频| e午夜精品久久久久久久| 高清在线国产一区| 无限看片的www在线观看| 国产一区二区三区视频了| 精品国产一区二区三区四区第35| 黄色视频,在线免费观看| 精品福利永久在线观看| 99热国产这里只有精品6| 高清在线国产一区| 亚洲av片天天在线观看| 老司机亚洲免费影院| 精品亚洲乱码少妇综合久久| 岛国毛片在线播放| 欧美成狂野欧美在线观看| av网站在线播放免费| 国产精品久久电影中文字幕 | 日韩免费av在线播放| avwww免费| 国产精品av久久久久免费| 免费在线观看黄色视频的| 国产精品一区二区免费欧美| 一级a爱视频在线免费观看| 99久久人妻综合| 大香蕉久久网| 黄片大片在线免费观看| 黄色丝袜av网址大全| av不卡在线播放| 欧美+亚洲+日韩+国产| 国产精品二区激情视频| 久久亚洲真实| 免费不卡黄色视频| 精品少妇内射三级| 国产av一区二区精品久久| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 欧美日韩av久久| 免费少妇av软件| 久久人妻福利社区极品人妻图片| 制服诱惑二区| 黑人欧美特级aaaaaa片| 成年版毛片免费区| 波多野结衣av一区二区av| 欧美日本中文国产一区发布| 后天国语完整版免费观看| 男女高潮啪啪啪动态图| 9色porny在线观看| 欧美激情久久久久久爽电影 | 亚洲国产精品一区二区三区在线| 日韩精品免费视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 欧美日本中文国产一区发布| 免费观看a级毛片全部| 我的亚洲天堂| 日韩一区二区三区影片| 久久精品aⅴ一区二区三区四区| 午夜福利,免费看| 最新的欧美精品一区二区| 99久久国产精品久久久| 国产成人免费无遮挡视频| 国产亚洲一区二区精品| 国产精品久久久人人做人人爽| 亚洲国产精品一区二区三区在线| 精品久久蜜臀av无| 久久青草综合色| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 啪啪无遮挡十八禁网站| 老司机午夜福利在线观看视频 | 久久久精品区二区三区| 性色av乱码一区二区三区2| 亚洲欧美日韩高清在线视频 | 老司机福利观看| 999久久久国产精品视频| 久久久久国内视频| 搡老熟女国产l中国老女人| 久久天堂一区二区三区四区| 1024视频免费在线观看| 一级毛片电影观看| 一级毛片女人18水好多| 亚洲精品中文字幕在线视频| 美女高潮喷水抽搐中文字幕| 国产黄频视频在线观看| 国产在线精品亚洲第一网站| 黄网站色视频无遮挡免费观看| 国产不卡av网站在线观看| 日韩人妻精品一区2区三区| 日韩中文字幕视频在线看片| 在线 av 中文字幕| 国产精品麻豆人妻色哟哟久久| 国产男女超爽视频在线观看| 视频区图区小说| 国产精品一区二区在线不卡| 黄色成人免费大全| 极品少妇高潮喷水抽搐| 久久中文字幕人妻熟女| 最新在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 日韩视频在线欧美| 我要看黄色一级片免费的| 手机成人av网站| 国产在线观看jvid| 人人妻,人人澡人人爽秒播| 国产成人系列免费观看| 国产精品av久久久久免费| 国产单亲对白刺激| 午夜激情av网站| 男人舔女人的私密视频| 国产三级黄色录像| 无人区码免费观看不卡 | 色尼玛亚洲综合影院| 欧美精品啪啪一区二区三区| 好男人电影高清在线观看| 天天添夜夜摸| videosex国产| 69精品国产乱码久久久| 精品午夜福利视频在线观看一区 | 国产色视频综合| 久久青草综合色| 纯流量卡能插随身wifi吗| 大香蕉久久网| 日本一区二区免费在线视频| 亚洲情色 制服丝袜| 不卡一级毛片| 国产真人三级小视频在线观看| 国产伦人伦偷精品视频| www.精华液| 亚洲色图av天堂| 亚洲第一av免费看| 交换朋友夫妻互换小说| 可以免费在线观看a视频的电影网站| 女性被躁到高潮视频| 婷婷成人精品国产| 十八禁网站免费在线| 欧美日韩视频精品一区| 国产精品一区二区在线不卡| 国产成人欧美| 男女床上黄色一级片免费看| bbb黄色大片| 亚洲欧美一区二区三区久久| 午夜两性在线视频| 日韩 欧美 亚洲 中文字幕| 精品高清国产在线一区| 国精品久久久久久国模美| 757午夜福利合集在线观看| 夜夜骑夜夜射夜夜干| 婷婷丁香在线五月| 国产av国产精品国产| 女人被躁到高潮嗷嗷叫费观| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 国产精品98久久久久久宅男小说| 国产精品免费一区二区三区在线 | 国产xxxxx性猛交| 麻豆av在线久日| 国产高清国产精品国产三级| 麻豆乱淫一区二区| 精品少妇黑人巨大在线播放| 成人手机av| 免费一级毛片在线播放高清视频 | 亚洲欧美精品综合一区二区三区| 黄片小视频在线播放| 日韩大码丰满熟妇| 我的亚洲天堂| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 免费在线观看黄色视频的| 99久久99久久久精品蜜桃| 精品国产超薄肉色丝袜足j| 精品少妇黑人巨大在线播放| 操美女的视频在线观看| 99国产精品一区二区三区| 久久亚洲精品不卡| 午夜福利免费观看在线| 如日韩欧美国产精品一区二区三区| 亚洲av成人一区二区三| 岛国毛片在线播放| 日日夜夜操网爽| 在线天堂中文资源库| 免费在线观看日本一区| 精品乱码久久久久久99久播| 91大片在线观看| 亚洲成人国产一区在线观看| 国产人伦9x9x在线观看| 久热这里只有精品99| √禁漫天堂资源中文www| 巨乳人妻的诱惑在线观看| av又黄又爽大尺度在线免费看| 亚洲国产毛片av蜜桃av| 狂野欧美激情性xxxx| 天天躁狠狠躁夜夜躁狠狠躁| 一本—道久久a久久精品蜜桃钙片| 女人高潮潮喷娇喘18禁视频| 极品人妻少妇av视频| 91成人精品电影| 性少妇av在线| 丁香六月天网| 最新在线观看一区二区三区| videos熟女内射| 亚洲免费av在线视频| 少妇精品久久久久久久| e午夜精品久久久久久久| 新久久久久国产一级毛片| 日本wwww免费看| 精品人妻1区二区| 夜夜夜夜夜久久久久| 日韩成人在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 97在线人人人人妻| 久久久久久亚洲精品国产蜜桃av| 波多野结衣av一区二区av| 亚洲第一欧美日韩一区二区三区 | av视频免费观看在线观看| 久久久精品94久久精品| 亚洲视频免费观看视频| 十分钟在线观看高清视频www| 一个人免费看片子| 国产区一区二久久| 欧美精品高潮呻吟av久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久狼人影院| 国产精品 欧美亚洲| 久久中文看片网| 黄色丝袜av网址大全| 日韩欧美三级三区| 精品少妇久久久久久888优播| 天堂动漫精品| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| 久久午夜综合久久蜜桃| 国产欧美日韩精品亚洲av| 亚洲情色 制服丝袜| 19禁男女啪啪无遮挡网站| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 老司机深夜福利视频在线观看| 亚洲九九香蕉| 久久精品亚洲熟妇少妇任你| 老鸭窝网址在线观看| 成年版毛片免费区| 欧美亚洲 丝袜 人妻 在线| 亚洲精品在线美女| 最黄视频免费看| 十八禁高潮呻吟视频| 欧美一级毛片孕妇| 精品视频人人做人人爽| 国产成人影院久久av| 欧美亚洲日本最大视频资源| 国产成人av教育| 女人高潮潮喷娇喘18禁视频| 亚洲第一欧美日韩一区二区三区 | 免费看十八禁软件| 老汉色∧v一级毛片| 黄色片一级片一级黄色片| 在线观看人妻少妇| 国产亚洲欧美精品永久| 在线永久观看黄色视频| 最新美女视频免费是黄的| 成年人免费黄色播放视频| 我的亚洲天堂| 亚洲av日韩精品久久久久久密| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 久久精品亚洲熟妇少妇任你| 中亚洲国语对白在线视频| 一区二区三区激情视频| 免费观看a级毛片全部| av天堂在线播放| 亚洲avbb在线观看| 男女床上黄色一级片免费看| 99国产精品99久久久久| 少妇裸体淫交视频免费看高清 | 丝袜在线中文字幕| 99久久精品国产亚洲精品| 另类精品久久| 性少妇av在线| 久久精品熟女亚洲av麻豆精品| 国产人伦9x9x在线观看| 欧美激情极品国产一区二区三区| 精品高清国产在线一区| 亚洲国产精品一区二区三区在线| 久久人人97超碰香蕉20202| 亚洲精品国产一区二区精华液| 露出奶头的视频| 一本—道久久a久久精品蜜桃钙片| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产单亲对白刺激| 男女无遮挡免费网站观看| 国产片内射在线| 亚洲成av片中文字幕在线观看| 丝袜美足系列| 51午夜福利影视在线观看| 麻豆乱淫一区二区| 黄色片一级片一级黄色片| 亚洲第一av免费看| 丰满迷人的少妇在线观看| 性高湖久久久久久久久免费观看| 热99国产精品久久久久久7|