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

    Modifications in aroma characteristics of ‘Merlot’ dry red wines aged in American,French and Slovakian oak barrels with different toasting degrees

    2024-02-20 07:14:46HaochengLuBinhaoChengYibinLanChangqingDuanFeiHe

    Haocheng Lu,Binhao Cheng,Yibin Lan,Changqing Duan,Fei He,*

    a Center for Viticulture and Enology, College of Food Science and Nutritional Engineering, China Agricultural University, Beijing 100083, China

    b Key Laboratory of Viticulture and Enology, Ministry of Agriculture and Rural Affairs, Beijing 100083, China

    Keywords:Wine Aroma Oak Barrel Oak species Toasting degrees Geographic origin Aging

    ABSTRACT Aging in oak barrels is widely used in enology which could bring flavor changes and aromatic complexity to wines.In the present study,the aroma compounds were analyzed from the ‘Merlot’ dry red wines,which were fermented in two types of fermenters (stainless steel tank and rotated oak barrel) and aged in six types of oak barrels (three geographic origins × two toasting degrees) for different time (0,3,6 and 9 months,respectively).Results showed that 30 volatiles were associated with barrels and increased during oak aging.The fermenters could influence the intensities of the toast,leathery,smoky,fruity,floral and caramel aromas.The concentration of whisky lactone,eugenol,cis-isoeugenol,and the intensities of the toast and spicy aromas were highest in the wines aged in American oak and were lowest in the wines aged in French oak barrels.The concentrations of guaiacol,syringol,trans-isoeugenol,furfural alcohol,vanilla,cis-whisky lactone enabled the medium toasting barrels to be distinguished from the light toasting ones.The compounds originating from the barrels could be used to distinguish the types of different barrels,but the other general grape-derived and fermentation-derived volatiles could not. The fermenters,oak species and toasting degrees of the barrels all had significant effects on the aroma profiles of the aged ‘Merlot’ dry red wines,but the influence of the geographic origin was not obvious.

    1. Introduction

    Aroma performance is a key sensory quality of the wine that directly influences the preference of consumers.Generally,wine aroma compounds could be derived from grapes,fermentation,and aging[1].Wine aging could be divided into two successive periods:barrel (or tank) aging and bottle aging,each of which had specific effects on the wine aromas.Barrel aging had been used for more than 200 years[2],and the corresponding effects on wine included not only the extraction of certain compounds from the wood,such as lactones,furanic compounds,phenolic aldehydes and volatile phenols,but also the formation of new volatile compounds[3],as well as the transformation,oxidation[4],hydrolysis-esterification[5-7],and sorption[8]of some compounds originating from the grapes and fermentation.Many researchers had demonstrated many factors which influenced the quantities of the aroma compounds associated with barrel aging,including the botanic oak species[9-11],the geographic origin of the oak[12-14],the wood grain tightness[11,15],the wood seasoning and toasting process[9],the ratio of the barrel surface to its volume[16],and the contact time between wines and barrels[4,14,17].

    Nowadays,France and America are the major sources of oaks for barrels.In Europe,the most popular oak species employed for cooperage areQuercusrobur(Q.pedunculata) andQ.sessilis(Q.petraea,andQ.sessiliflora),whileQ.albais the most widely distributed American oak species.The most significant feature that could differentiate French oak from American oak wascis-whisky lactone (also named oak lactone) and the ratios ofcis-/trans-whisky lactone[14,18,19]were higher in the latter.With the development of the wine industry,the limited production and the relatively high price of the French oaks in particular,oaks from other European countries had attracted increased interest.These countries include Spain,Hungary and Slovakia,which had the same oak species as France.Besides,some Russian oaks (Q.sessilis) low in lactones were also used.Such oaks could be referred to as East European oaks.Andrei Prida[14]found that the composition of East European oaks was in between that of French and American oaks,and could be distinguished from French oaks by 2-phenylethanol,eugenol and vanilla.

    The degrees of toasting during barrel-making were another factor that influenced the aging aromas of the wines.The concentration of many oak aroma compounds,such as furanic compounds,phenolic aldehydes and volatile phenols could be increased by toasting[20].But as toasting continues,the concentrations of some compounds,such as vanillin,syringaldehyde and whisky lactones decrease through volatilization and thermal degradation[20].

    Just as the differences between oak barrels could influence the aroma compounds in the wines,wine composition itself could also influence the extraction of volatile compounds from oak woods.Higher alcohol levels could increase the extraction of whisky lactones and vanilla[9,21].Other influencing factors included the SO2concentration[22],ellagitannins[21],polyphenolic compounds[21].

    In China,extensive studies have been performed regarding the techniques to improve wine quality.Most of the studies focus on the viticulture or the fermentation stage.There is much space for investigating the choice of oak barrels during wine aging.Although many previous studies have investigated the influence of oak barrels on wine composition.There were few reports in the Yanqing region with a typical continental monsoon climate.In the present research,the ‘Merlot’ wines fermented in tanks or oak barrels were aged for 0,3,6 and 9 months in new AmericanQ.alba,French and SlovakiaQ.petraeabarrels,of light or medium toasting,to follow the evolution patterns of the aroma compounds during aging and to determine the characteristic compounds those were obtained from the different barrels.

    2. Materials and methods

    2.1 Samples

    In 2012,‘Merlot’ grapes were harvested at commercial maturity with a total acidity of 6 to 7 g/L and total soluble solids of 24-25 °Brix from a vineyard in Yanqing County (40°45’34’’N,115°97’2’’E),Beijing.Vines in this winery were planted in the year of 2000 and grown on their own roots under a north-south row orientation with a modified vertical shoot positioning (M-VSP) trellis system.The harvested grapes were immediately transported to the winery,destemmed and crushed,and then sulfur dioxide (60 mg/L) was added to the must.

    Two types of fermenters were used in this study: Type 1: stainless steel tanks (TF) of 20 hL capacity.Type 2: rotated oak barrels (BF)made from French (Q.petraea) oaks with a volume of 20 hL,which had been used for seven years.All fermenters were located in the same fermentation environment.The fermenters were filled to about 75%-80% capacity.

    After loading in the fermenters,20 g/hL of the commercial yeast L2323 (LALLEMAND,France) was added after rehydration to each fermenter to start the alcoholic fermentation.The fermentation temperature was kept at 20-25 oC with full punching down three times daily.When the wine reached the end of the alcoholic fermentation(the reducing sugar was below 4 g/L),the wine was pressed off the skins and seeds.The malolactic fermentation was then started by inoculation with the commercial malolactic bacteria LALVIN 31(LALLEMAND,France).After the malolactic fermentation,the wines were settled in the fermenter for five days and then the free-run wine was racked off the gross lees and transferred into oak barrels of different types.

    The oak woods used for new barrels were harvested from Missouri in the USA (Q.alba,AO),France (Q.petraea,FO) and Slovakia (Q.petraea,SO),respectively.The barrels made from these oaks were toasted at two levels using traditional methods.The starting temperature of the light toasting program was 180 oC for 10 min,and then 200 oC for 20 min.For medium toasting,the starting temperature was 180 oC which was held for 10 min,then raised to 220 oC for 35 min.All the oak barrels were the same Bordeaux style and had the same volume of 225 L.

    The two wines from the two different fermenters were put into the six types of new barrels with a duplicate for each type.And also put into stainless tanks as the control.The wines aged in light toasting barrels,medium toasting barrels,and stainless tanks were coded as LT,MT,and CK,respectively.There were three factors investigated in our study: 1) the fermenter type,which was coded with TF and BF;2) the barrel geographic origin,which was coded with AO,FO,and SO;3) the toast levels,which was coded with MT and LT.The CK was represented for wines aging in stainless tanks.In the following analysis,we used ‘XX-XX-XX’ to respent for “the fermenter typethe barrel geographic origin-the toast level” for the fourteen wines:TF-AO-MT,TF-AO-LT,TF-SO-MT,TF-SO-LT,TF-FO-MT,TF-FO-LT,TF-CK,BF-AO-MT,BF-AO-LT,BF-SO-MT,BF-SOLT,BF-FO-MT,BF-FO-LT and BF-CK.They were aged for nine months and were sampled four times: 0,3,6,and 9 months.For each wine,there were two biological replicates × two technical replicates.These samples were stored in freezers at -20 oC until the aroma compounds of these samples were analyzed.

    2.2 Physicochemical parameter detection

    The enological parameters of the samples,including alcohol level,residual sugar,total acidity and pH were analyzed based on the methods given by GB/T 15038-2006 (the National Standard of the People’s Republic of China for wine,2006).In brief,the alcohol content was determined by the distillation/densimetry method.The residual sugar was detected by titration with Fehling’s reagent.Wine pH was determined using a pH meter (Sartorius PB-10,Germany).Total acidity of wines was titrated with 0.05 mol/L NaOH to an endpoint of pH 8.2 and was expressed as tartaric acid equivalents.

    2.3 Headspace solid-phase extraction

    The majority of aromatic compounds from the wine samples were extracted by headspace solid-phase microextraction (HS-SPME)coupled with gas chromatography-mass spectrometry (GC-MS)[23-25].5 mL wine mixed with 10 μL of 2.004 g/L 4-methyl-2-pentanol(internal standard) and 1 g NaCl was prepared in 20 mL vial capped with a PTFE-silicon septum.Then the vial containing the sample was moved to the heating equipment.After equilibrating at 40 °C for 30 min,the SPME fiber was inserted into the headspace of the vial to extract the volatile compounds at 40 °C for 30 min.The heating/stirring equipment was stirred at 500 r/min.The SPME fiber was thermally desorbed by being inserted into the GC injector for 8 min.

    2.4 Liquid-liquid extraction

    Liquid-liquid extraction was only used to extract compounds mainly extracted from the oak barrels.The 25 mL of wine was put into a 50 mL centrifuge tube,5 g/L ammonium sulfate was added and the tube was shaken until the ammonium sulfate was completely dissolved.The internal standards in this part wereo-vanillin(2.002 mg/L),3,4-dimethylphenol (1.027 mg/L) andγ-caprolactone(2.043 mg/L).Every 25 mL sample had 40 μL of mixed internal standards added.Then,5 mL of dichloromethane was added to the solution.The solution was shaken to uniformly mix the sample and then centrifuged 15 min at 8 000 r/min.The organic liquid layer was decanted to another tube.Then a further 5 mL of dichloromethane was added to the inorganic solution,shaken and centrifuged,before being added to the first extract.The combined dichloromethane extract was dried with anhydrous sodium sulphate,and transferred to another tube.That solution was concentrated to a final volume of 1 mL under an N2stream,prior to being analyzed by GC-MS.

    2.5 GC-MS Conditions and volatile compounds analysis

    An Agilent 6890 GC coupled to an Agilent 5975 mass spectrometer was used.The capillary column was 60 m × 0.25 mm and 0.25 μm film thickness (J &W Scientific,Folsom,CA,USA).The flow rate of the helium carrier gas was 1 mL/min.In order to inject the samples,the SPME fiber was placed at GC inlet,for 25 min,which was operated in the splitless mode.The oven temperature program was: 1 min at 50 oC then raised to 220 oC at 3 oC/min and held at this temperature for 5 min.The mass detector was operated in the full scan mode (m/z30-350) with electron ionization (EI) mode at 70 eV.Peak areas were found by setting the Chemstation software (Agilent Technologies) to autotune mode.Compound identification was based on retention indices of reference standards combined with mass spectra matching in the standard NIST 11 library.When reference standards were not available,tentative identification was made based on mass spectra matching in the standard NIST 11 library,or based on the comparison of the retention indices sourced from the NIST Standard Reference Database (NIST,Chemistry Web Book)[23].

    For the quantitative analysis,a synthetic matrix was synthesized according to the basic enological parameters of the sample wines.It was made from distilled water with 14% (V/V) ethanol (HPLC quality),5 g/L tartaric acid,and the pH was adjusted to 3.8 with NaOH solution.All standards were dissolved in ethanol (HPLC quality) and mixed.Then the combined solution was diluted by the model wine to 15 levels.For the general volatiles,4-methyl-2-pentanol was added into the combined solution as the internal standard.For oak-related volatiles,o-vanillin was used as the internal standard for furans and lactones,3,4-dimethylphenol was used as the internal standard for volatile phenols andγ-caprolactone was used as the internal for phenolic aldehydes.The extraction and analytical procedures were identical for these samples.The regression coefficients of the calibration curves were all above 99%.If there weren’t standards for some volatile compounds,the concentrations were calculated from the calibration curves of compounds with the closest chemical structures and the number of carbon atoms.The detailed quantification information is in Supplementary Table 1.

    Table 1 Physicochemical parameters of the original wines fermented in the two fermenters before aging.

    2.6 Odor activity values (OAVs) and aroma series

    Odor activity values (OAVs) were used to represent the sensory impact of the aroma compounds.The OAVs were calculated by dividing the concentration of the compounds by their sensory thresholds[24,26,27].If the compounds were present at abovethreshold levels (OAV >1),these compounds were recognized to have the potential to actively contribute to the aroma of the wines[28,29].If the OAVs were lower than 1 but higher than 0.1,these compounds could also potentially contribute to the wine aromas by additive effects[28].But if their OAVs were lower than 0.1,their effects on wine aromas could be negligible[30].

    The active aroma compounds were divided into thirteen series with similar odor descriptors.From this,it should be possible to predict the overall aroma profiles of the wines based on their numerical GC-MS analysis data.The thirteen series were: (1) fruity,(2) floral,(3) herbaceous,(4) nutty,(5) caramel,(6) earthy,(7)chemical,(8) fatty,(9) toast,(10) oaky,(11) spicy,(12) smoky and(13) leathery.They were modeled from the literature and the widely employed wine aroma wheel[31].

    Some compounds were allocated to more than one series[32,33],because of the great complexity of olfactory perception.Using the OAVs of each compound in a given series,it was possible to model the total intensities of each aroma series and hence the series contribution to the overall aroma.This would also allow the correlation of the quantitative information obtained by the chemical analysis with the sensory properties of the samples.

    2.7 Statistical analysis

    Statistical software SPSS 19.0 was used to perform the factor analysis and identify general compounds and compounds associated with barrel aging.Before the analysis,all variables were meancentered and scaled to unit variance.This multivariate analysis method reduced one set of variables into a lesser number of new variables which could represent the original variables and found the associations among the original variables.To better interpret these factors,the principal components were used to extract factors from all the original variables after varimax rotation of the data[21,34,35].The factor values for the samples drawn from the wines during aging were calculated to facilitate an objective comparison during the wines development.The value of each factor was calculated by multiplying the coefficients of volatile compounds in the factor loading matrix by their corresponding concentrations.The Figures were prepared by using Origin 9.0 software.

    3. Results and discussion

    3.1 Physicochemical parameters

    The enological parameters of the original wines fermented in the two fermenters (TF and BF) before aging were shown in Table 1.The total acidity and pH in the original wines showed no difference indicating the similar maturity of grapes fermenting in the two fermenters,although the tank fermented wines had a little lower alcohol concentration than the barrel fermented wines.

    3.2 Effects of aging in different barrels on compounds associated with barrels

    As shown in Supplementary Table 1,33 volatiles that resulted from oak aging were detected by liquid-liquid extraction,including 17 volatile phenols,5 phenolic aldehydes,9 furan,derivatives and 2 lactones.There were only 17 compounds of OAVs >0.1 associated with oak barrels and their concentrations were shown in Supplementary Table 2.Most of the selected compounds had increased concentrations in the first three months of the aging period and slowed down in the last six months,except for ethyl phenols and vinyl phenols.When selecting the influencing factors using the factor analysis,it had been suggested to explain at least 80% of the total variances[21].We found that five factors,which included all of the 17 volatiles mentioned above,had higher eigen values than 1,and explained 88.8% of the total variance 1.Using varimax rotation,the coefficients of all variables in the five factors were calculated and were shown in Table 2.To simplify the presentation of the result,only coefficients >0.50 were listed.

    Table 2 Factor loading matrix after varimax rotation of volatile compounds associated with oak barrels.

    Factor 1 explained 29.2% of the total variance of all the positive variables.Most of the compounds incorporated in factor 1 originated from greenwood,except for vanillin and furfural,which came from the degradation of lignins and carbohydrates during the toasting process.As shown in Fig.1a,there was little difference between the wines fermented in the stainless steel tank and oak barrel.This might be because the oak barrel has been used for seven years and a previous study showed that the amount of extractable compounds in the oaks of the oak fermenter could be expected to be minimal after 8 years’utilization[36].The factor values of all the wines increased quickly in the first three months due to the rapidly increased concentrations of eugenol,trans-eugenol,cis-andtrans-whisky lactone.The increased concentration ofcis-whisky lactone in AO was most obvious.After aging for three months,the concentration ofcis-whisky lactone was 195.81 μg/L in TF-AO-MT,266.38 μg/L in TF-AO-LT,256.47 μg/L in BF-AO-MT and 294.06 μg/L in BF-AO-LT,which was in agreement with the reports from Ortega-Heras[21]and Pérez-Prieto[37].Over the next 6 months,the factor values of wines fermented in tanks increased,which was contrary to the barrelfermented wines.For the duration of aging,the factor values of the barrel-fermented wines were higher than the tank-fermented wines.Furfural was the principal compound giving these results.However,at the end of aging (9 month),the factor values of the barrel-fermented wines were no more higher than the tank-fermented wines because of the decreased concentration of furfural.This might because that the furfural could be reacted to furfural alcohol[38],2-furanmethanethiol[39],and adducts with (+)-catechin[40].

    Factor 2 included five compounds originating in the toasting process,although the 4-ethylphenol,which increased during aging,was likely to come mainly from decarboxylated cinnamic acids due to the action ofBrettanomycesandDekkerayeast[21,41].When its concentration was greater than its threshold (620 μg/L)[42],it had a negative sensory effect.Theo-cresol and 4-propylguaiacol,originating from the thermal dehydration during toasting[43,44],were positively correlated with factor 2.Acetovanilone was also a product of the toasting process[45],and had a negative association with ethyl vanillate.The values of factor 2 (Fig.1b) were stable or declined slightly during the first three months and were less than zero,because the extraction of acetovanilone and ethyl vanillate was more rapid than volatile phenols[46],and the concentration of 4-ethylphenol was still low.But after the first three months,the values of these wines increased quickly,approximately linearly.Non-barrel origin 4-ethylphenol,and other volatile phenols were responsible for this phenomenon.

    As was shown in Table 2,factor 3 included four variables,all of which had positive contributions.Thetrans-whisky lactone,with its wood and coconut aromas,andtrans-isoeugenol,contributing clove and spicy aromas[28],which existed in unseasoned wood,were the markers of the oak aged wines.Vanillin originated from the degradation of lignins during the aging period and 2-furfury alcohol was transformed from furfural.In factor 3,the factor values of the tank fermented wines were higher than those of the wines fermented in the barrels (Fig.1c).For the latter wine,the factor values increased slowly during the first 6 months but more rapidly for the last 3 months.In contrast,the factor values of the wines fermented in the tank increased quickly during the first 3 months and then slowly for the last 6 months.Furfuryl alcohol,which was biologically formed from furfural but not extracted from the wood[47,48],was the cause of the difference.For example,after aging for 3 months,the concentration of furfuryl alcohol was 8,027.96 μg/L in TF-AOMT,nearly one hundred times higher than that in the original wine(84.57 μg/L).At the same time,the concentration of furfuryl alcohol in BF-AO-MT was only 2154.36 μg/L.But after aging for 9 months,the concentration of furfuryl alcohol in BF-AO-MT wines was augmented to 10,814.78 μg/L and there was 10,527.54 μg/L in TF-AO-MT wines.As shown by factor 1,the timing of furfural transformation to the corresponding alcohol was seemingly dependent on the fermenters used for the base.

    The variance explained by factor 4 (14.2%) was similar to that explained by factor 3.It included three compounds,all of which came from the toasting process.The values of this factor (Fig.1d) were consistent in most of the wines.The values were stable during 0 to 3 months,increased rapidly in the next 3 months.Then the treatments showed two distinct trends.The values kept increasing in BF-CK,TFCK,TF-AO-MT,TF-AO-LT,TF-SO-MT,BF-FO-LT while declined quickly during the last 3 months in the last treatments.Guaiacol was the main contributor to these changes.For example,TF-SO-LT wines had the highest factor value when aging 6 months with 816.61 μg/L guaiacol,but there was only 75 μg/L guaiacol in BF-FL-LT after 6 months of aging,which had the lowest factor value.The declining trend of guaiacol was also found by Chira et al.that the maximum extraction of guaiacol occurred during the first nine months and declined then at the end of the aging point[43].Two possible reasons might lead to the decline of guaiacol: (i) the guaiacol in the wine and oak barrels was a dynamic balance that could be adsorbed on the surface of the oak barrel;(ii) during aging,the guaiacol could form polymeric pigments with anthocyanins which led to its decrease[49].

    Factor 5 comprised two volatile phenols: 4-vinyl guaiacol and 4-vinyl phenol.Hydroxycinnamate decarboxylase fromBrettanomycesconverted ferulic acid,p-coumaric acid,caffeic acid,the cinnamic acids of grapes and wine into vinylphenol,which in turn was converted into ethylphenols by vinylphenol reductase[50].So the concentration of ethylphrnol had a significant effect on the wine quality.The values of wines (Fig.1e) fermented in the barrels were lower than the tank fermented ones,most of which stayed stable during the aging period.The difference of the barrel fermented and tank fermented wines might be due to the action ofBrettanomyces,which was normally considered as unwanted spoilage yeast,although in all wines the concentration of vinylphenols and ethylphenols were below their threshold.

    The principal component analysis of all compounds associated with oak barrels was used to analyze the influence of the species and the geographic origins of oak barrels (Fig.2a and Fig.2b).Factor 1 and factor 2 were examined further because they explained 47.2% of the total variance.The positive correlation of eugenol,cis-andtrans-isoeugenol,cis-whisky lactone,in factor 1 (Fig.2) showed that American oak barrels could be distinguished from European barrels,which was in agreement with the findings of Prida[14]and others who worked with different species[36,51].A previous research reported that furanic aldehydes had the potential for distinguishing the oak samples from different species[52],and furfural was found to have a positive impact on factor 1 in the present study.That might result from the different concentrations of carbohydrate in the green wood of the different oak species.As shown in Fig.2,the samples from French and Slovakian barrels could not be separated by factor 1 and factor 2.This lack of discrimination in the geographic origins was,at least in part,due to the fact that the same oak species were used,and the side of the forest where the trees came from,the part of the tree used for the staves,and the cutting and seasoning of the staves were very similar[20].However,the stepwise discriminant analysis did distinguish these two geographic origins with the discriminant rates greater than 90%.The structure of the stepwise discriminant function matrix revealed the compounds which could distinguish Slovakia and French barrels (Table 3).The concentrations ofcis-whisky lactone,eugenol andcis-isoeugenol were higher than in the American barrels than those in the Slovakian barrels and were lost in French barrels,which was again in agreement with Prida[14],partially in agreement with Pizarro[52],who considered that the French barrels had the highestcis-eugenol.Thetrans-whisky lactone levels were highest in the wines aged in the Slovakian oak,followed by the French barrels and then the American barrels.Prida[14]suggested that eugenol and vanilla could be used as markers to differentiate the French and the East European oaks,which was in agreement with our results.Neither the geographic origins nor the toasting levels affected the concentrations of vanilla and ethyl vanillate.

    Fig.2 Results of PCA for the barrels of different geographic origins (a,b) and toasting degrees (c,d).a: Loading plot for the wines aged in the barrels of the three geographic origins with the wines sampled at 3,6 and 9 months.b: Scores plot for the volatile compounds associated with the oak barrels.c: Loading plot for the wines aged in the barrels of the two toasting degrees sampled at,3,6,and 9 months.d: Scores plot for the volatile compounds associated with the oak barrels.AO,the wines aged in the America oak barrels;FO,the wines aged in the French oak barrels;SO,the wines aged in the Slovakian oak barrels.MT,the wines aged in the medium toasting oak barrels;LT,the wines aged in the light toasting oak barrels.

    Table 3 Structural matrix of the stepwise discriminant function.

    Considering that the compounds in factors 2 and 4 arose from the toasting process,these two factors were used for a PCA analysis(Fig.2c and Fig.2d).At the same time,the stepwise discriminant function was used to indicate the results (Table 3).Six compounds,including guaiacol,syringol,trans-isoeugenol,furfural alcohol,vanilla,andcis-whisky lactone,distinguished the light toasting barrels from the medium toasting ones (Fig.2d).It was well known that toasting oak altered the concentrations of many compounds[20,32].Volatile phenols and phenolic aldehyde were formed by the thermal dehydration of lignins[20,21].Furfural could be transformed into furfural alcohol during aging originating from the thermal dehydration of carbohydrates[21,52],and part of whisky lactone arose from the lipid oxidation during the toasting process[52,53].For this reason,these compounds could indicate the difference in toasting.It was found that the concentration oftrans-isoeugenol increased in the wines aged in the medium toasting barrels.

    The concentrations of the above aroma compounds were analyzed and then the sum of OAVs were used to describe the special aroma characteristics of various barrels (Supplementary Table 3).The main aroma descriptors were nutty,roasted,oaky,spicy,smoky and leathery (Supplementary Table 1).Furfural was the only compound that had a nutty aroma.The nutty aroma of wines from medium toasting barrels was always stronger than in the wines from light toasting ones because of the source of furfural[21,52].The nutty aroma had the lowest intensity among all of the aroma series,just from 0 in TF-CK to 0.14 in TF-AO-MT.The roasted aroma,which included the contributions from furfural and furfuryl alcohol,was higher in the wines fermented in the tanks and the wines aged in medium toasting barrels.During the first 3 months of aging,this aroma type increased then stayed steady or increased slowly over the next 3 months,and had no consistent trend in the last 3 months.The oak aroma was associated with eugenol,cis-andtrans-whisky lactone,vanilla and ethyl vanillate,so the wines from the American barrels were highest on this character.For example,the OAV of oak in BF-AO-LT was 23.63,much higher than those of 9.33 in BF-SO-LT and 10.15 in BFFO-LT.The oak aroma of wines from the light toasting American and Slovakian barrels was higher than the medium toasting ones because of the concentration of whisky lactone.Though the three volatile phenols with smoky aromas all came from toasting,the intensity of the smoky aroma did not reflect the level of toasting.The greatest increase of this aroma took place during 3 to 6 months,but decreased or did not vary during the last 3 months.The spicy notes were mainly due to volatile phenols and phenolic aldehydes,and were in keeping with factor 4.Ethyl phenol had leathery aroma,and was identified with factor 5.Some compounds associated with barrels had caramel aroma but were of low intensity,so they were combined together and were analyzed as below.

    3.3 Effects of aging in different barrels on general compounds

    Apart from compounds closely associated with barrels,35 general volatiles were selected with OAVs higher than 0.1.Although most of these compounds came from grapes and fermentation,some of them,such as ethyl lactate,ethyl 2-hydroxy-4-methylpentanoate,isoamyl lactate and 3-methyl-1-pentanol were generated during aging.The concentrations of most higher alcohols,esters and fatty acids decreased during the first 3 aging months.But most of the terpenes,norisoprenoids and even ethyl 2-hydroxy-4-methylpentanoate and ethyl benzoate stayed stable throughout aging.

    Factor analysis was used to reveal the similarities between the wines and found that 90.9% of the total variance could be explained by the six factors that included all of the 35 compounds.It was also found that even factor 6 had an eigenvalue above 1.The coefficients of the variables,Table 4,were higher than 0.50 after varimax rotation of the six factors.

    Table 4 Factor loading matrix after varimax rotation of the primary and the secondary compounds.

    Factor 1 explained 34% of the total variance and was associated with nearly half of the total number of the fermentation-related variables,mainly including higher alcohols and ethyl esters (Table 4).From Fig.3a,the wines which were fermented in the tanks could be distinguished easily from the barrel fermented ones and had higher concentrations than the latter during aging.The decrease in values in wines during the first 3 months of aging might be due to ester hydrolysis.For example,the concentrations of ethyl acetate were 110 359.49 μg/L in TF and 116 179.11 μg/L in BF,respectively.After aging in the barrel for 3 months,the concentration of ethyl acetate ranged from 67 085.96 μg/L to 98 565.5 μg/L in all of the wines.Most of the esters contribute floral and fruity notes to the aroma of the wines[1,28],and were biochemically synthesized more than formed as a result of reaching equilibrium concentrations from acid and alcohol esterification during fermentation[54,55],and gradually hydrolyzed to their corresponding alcohols and acids to approach to an equilibriumduring aging[53].Also in this period,the concentration of higher alcohols decreased to produce their corresponding carbonyl compounds and other compounds.The concentration of propyl alcohol and isobutanol decreased quickly.There were 71 149.36 μg/L of these compounds in BF before aging.Three months later,there were only 37 438.75 μg/L of them in BF-AO-MT,which was nearly half of the original wine.From 3 to 9 months of aging,the factor values of each wine were stable because the concentration of higher alcohol augmented with the concentration of esters declining.

    Fig.3 Changes in the values factors 1 to 6 for the compounds associated with grapes and fermentation in the wines aged in different oak barrels for 9 months.In the legend,‘XX-XX-XX’ represented “the fermenter type– the barrel geographic origin -the toast level” for the fourteen wines.BF,the wines fermented in the oak barrels;TF,the wines fermented in the stainless steel tanks;AO,the wines aged in the America oak barrels;FO,the wines aged in the French oak barrels;SO,the wines aged in the Slovakian oak barrels;MT,the wines aged in the medium toasting oak barrels;LT,the wines aged in the light toasting oak barrels.

    Factor 2 included C6-alcohol,esters and norisoprenoids and the factor values were below 0 for all of the wines,which originated from grape maturation or fermentation process.From Fig.3b,the wines from the tank fermenters could be separated from the oak barrel fermented wines.And the variation trends of the values for the wines that came from different fermenters were distinctive.The factor values of the tank fermented wines declined in the first 6 months and had different variation trends in the last 3 months,because the concentration of 1-hexanol decreased and the concentration of acetoin increased (Supplementary Table 2).Furthermore,1-hexanol had a positive correlation with factor 2 and was the single compound with an OAV higher than 0.1 in C6-alcohols,which came from the grapes and usually had vegetal and herbaceous characteristics,causing negative effects[55,56].It might transform to the corresponding esters to reduce the vegetal aroma during the aging period.However,the factor values of the wines fermented in barrels increased in the first 3 months but fell at the 6 months,caused by the concentration changes of isoamyl acetate andβ-damascenone.The latter compound,with its floral and fruit aromas and low threshold[28],was one of the most important norisoprenoids and played an essential role in the varietal character of the wines.It had a positive correlation with factor 2.Isoamyl acetate was also an important compound that had obvious banana characters and could be used to distinguish the wines of different aging time as reported by Castro-Vázquez[4].

    All variables,principally fatty acids and esters,had positive correlations in factor 3.The production of fatty acids was determined by grapes and fermentation conditions[28].Previous research indicated that yeasts produced more fatty acids if the fermentation conditions were anaerobic[57].According to some reports[4,6],fatty acids,which had fatty,fruity,cheese and rancid characteristics[28],decreased during the aging period,because of sorption by the barrels.This could explain the reason why the factor values of the wines decreased during the first 3 months.These factors also included ethyl lactate,ethyl 2-hydroxy-4-methylpentanoate and isoamyl lactate,which were mainly formed during the aging period and led to the increase of the factor values during 3 to 9 months.

    From the values of factors 1 to 6,it could be found that the concentration of the general aroma compounds fluctuated with the changes of the oak species,the geographic origins and the toasting degrees of barrels,but there were no consistent trends in each independent variable.To further evaluate the influence of the oak species and the geographic origins on these general compounds,the above factors were used to do a PCA analysis.But this did not distinguish the wines from the two oak species and the three different geographic origins.Then,the discriminant function was used for that purpose but had the same outcome.The results showed no general compounds that could be used as the markers to distinguish the American,French and Slovakian oaks in the present study.When these statistical analysis methods were used to evaluate the influence of toasting degrees on the general compounds,the same results were obtained.

    The sum of OAVs were also used to describe the special aroma characteristics of various barrels.The main aroma series associated with general compounds were fruity,floral,herbaceous,caramel,earthy,chemical and fatty (Supplementary Table 3).The intensity of the fruity aroma decreased in the first 3 months because of the dehydration of esters[53],but stayed stable from 3 to 9 months because of the increase of 1-propanol,and ethyl ester.The variation in floral and caramel intensities was almost the same,and the wines from different fermenters were inseparable because they all includedβ-damascenone and phenylacetaldehyde at concentrations greater than their thresholds and masked the effects of other compounds.The influence of the oak species,the geographic origins and the barrel toasting levels on the fruity,floral and caramel aroma series were not distinct.There was no regular pattern for all of the wines during the whole period when the herbaceous aroma was analyzed,which related to C6-alcohols,higher alcohols and carbonyl compounds.Another special aroma series associated with general compounds was the earthy aroma,which correlated only with 1-octen-3-ol,whose organoleptic property was described as the mushroom.The wines fermented in different fermenters could be distinguished before aging,but these wines changed individually during aging beyond 3 months.At the end of the aging period,the earthy aroma of the tank fermented wines augmented evidently compared to the beginning of the barrel aging.The concentration of 1-octen-3-ol fluctuated from 4.95 μg/L to 5.63 μg/L after 9 months of aging,which was higher than 3.72 μg/L in the tank fermented original wine (Supplementary Table 2).Compounds that contributed to fatty aroma series were not only associated with fatty acid,but also some esters and carbonyl compounds.The fatty aroma notes of the barrel fermented wines aged in the American oak barrels increased more rapidly during the last 3 months.But the influence of the barrels with their various parameters of the oak species,the geographic origins and the toasting degrees were not obvious on the aroma series associated with general compounds during the whole aging period.

    3.4 Aroma profiles of all wines at the end of the aging period

    After 9 months,at the end of the trial aging period,the sum OAVs of all thirteen aroma series were used to construct the aroma profiles of the various trial wines,as shown in Fig.4.The floral,caramel and fruity aromas dominated the profiles,accounting for 27.1%-32.7%,25.2%-30.1% and 12.6%-17.3%,respectively.The wines aged in the medium toasting French barrels had the highest floral and caramel aromas,and the lowest was in the wines aged in the American barrels.The wines fermented in the tanks had the highest floral aroma.While the fatty and chemical aroma notes also played important roles in all wine profiles,the wines were not obviously similar.Another four aroma series to make important contributions,but much less than the first three aroma series,were smoky,leathery,spicy and oaky,which were the aromas associated with the barrels.The oak aromas of the wines from most of the light toasting barrels were stronger than those of the medium toasting ones,except for the wines from the barrel fermenters and aged in the Slovakian barrels.The toasting degrees and the fermenters did not influence the spicy aroma.The intensity of the smoky aroma was higher and similar for the wines fermented in the tanks.The leathery aroma was similar to the smoky aroma,except for the wines from the barrel fermenters and aged in the light toasting Slovakian barrels.In all,each wine had its own organoleptic features after 9 months’ aging,according to the types of the fermenters used,the oak species,the geographic origins and the degrees of the barrel toasting.

    Fig.4 The spider plot for the different categories of aroma types associated with the grapes and fermentation process for the wines aged in different barrels for 9 months.In the legend,‘XX-XX-XX’ represented “the fermenter type-the barrel geographic origin-the toast level” for the fourteen wines.BF,the wines fermented in the oak barrels;TF,the wines fermented in the stainless steel tanks;AO,the wines aged in the America oak barrels;FO,the wines aged in the French oak barrels;SO,the wines aged in the Slovakian oak barrels;MT,the wines aged in the medium toasting oak barrels;LT,the wines aged in the light toasting oak barrels.

    4. Conclusion

    After aging,the ‘Merlot’ red wines fermented in the tanks had higher toasting notes,and more leathery and smoky aromas but lower fruity,floral and caramel aromas than those fermented in the barrels.After 3 months of barrel aging,the concentrations of most esters,higher alcohols and fatty acids decreased noticeably.In contrast,some compounds,such as whisky lactone,vanilla,eugenol and so on had been quickly extracted from the oak barrels in the same period.However,the rates of extraction slowed during the following 6 months.General aroma compounds,such as higher alcohols and esters,could not be recognized as the markers to distinguish the oak species,the geographic origins and the toasting degrees of the barrel oaks,though their concentrations could be impacted.The wines from the American oak barrels had the highest total amounts of whisky lactone,eugenol,cis-isoeugenol,as well as the highest toast and spicy aromas.The high concentration oftrans-whisky lactone could be used as a marker of the Slovakian barrels.The difference between the American and European barrels was greater than that between the French and Slovakian ones,reflecting that the influence of oak species was greater than the importance of the geographic origins.The higher concentrations of guaiacol,syringol,trans-isoeugenol,furfural alcohol,vanilla,ciswhisky lactone and the total whisky lactone could distinguish the medium toasting degree barrels from the light toasting ones.The former could bring more toast aroma,but a less smoky aroma to the wines.

    Conflict of interest statement

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    The authors gratefully acknowledge the financial support received from The Key Project of R&D Program of Ningxia Hui Autonomous Region,China (2022BBF01003) and China Agriculture Research System of MOF and MARA (CARS-29).

    Appendix A.Supplementary data

    Supplementary data associated with this article can be found,in the online version,at http://doi.org/10.26599/FSHW.2022.9250032.

    999精品在线视频| 侵犯人妻中文字幕一二三四区| 一进一出抽搐动态| 男女床上黄色一级片免费看| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 天堂动漫精品| 久久国产亚洲av麻豆专区| 久久中文字幕人妻熟女| 久热爱精品视频在线9| 99久久精品国产亚洲精品| 日韩成人在线观看一区二区三区| 国产色视频综合| 国产欧美日韩精品亚洲av| 免费在线观看影片大全网站| 国产单亲对白刺激| 男女之事视频高清在线观看| 日韩精品青青久久久久久| 久久国产精品人妻蜜桃| 久久青草综合色| 男女下面插进去视频免费观看| 国产av在哪里看| 免费在线观看视频国产中文字幕亚洲| 欧美色欧美亚洲另类二区 | 午夜激情av网站| 99久久国产精品久久久| 丝袜美腿诱惑在线| 国产一级毛片七仙女欲春2 | 国产欧美日韩一区二区三区在线| 亚洲 欧美一区二区三区| 国产国语露脸激情在线看| 国产av精品麻豆| 精品久久久久久,| 女性被躁到高潮视频| 国产精品久久电影中文字幕| 久久精品成人免费网站| 91麻豆精品激情在线观看国产| 国产午夜福利久久久久久| 成人18禁高潮啪啪吃奶动态图| 亚洲精品久久国产高清桃花| 1024香蕉在线观看| 黄网站色视频无遮挡免费观看| 亚洲av电影不卡..在线观看| 久久久久久久精品吃奶| av在线播放免费不卡| 久久久久国产精品人妻aⅴ院| 99久久综合精品五月天人人| 亚洲五月天丁香| av欧美777| av电影中文网址| 多毛熟女@视频| 国产欧美日韩一区二区精品| 一个人免费在线观看的高清视频| 999久久久国产精品视频| 黑人操中国人逼视频| 自拍欧美九色日韩亚洲蝌蚪91| 91av网站免费观看| 国产在线观看jvid| 成年版毛片免费区| 一卡2卡三卡四卡精品乱码亚洲| cao死你这个sao货| 免费高清在线观看日韩| 日本精品一区二区三区蜜桃| 欧美午夜高清在线| 男女做爰动态图高潮gif福利片 | 岛国视频午夜一区免费看| 亚洲国产精品999在线| 亚洲中文字幕一区二区三区有码在线看 | 免费看a级黄色片| 在线观看www视频免费| 成人国产综合亚洲| 老司机深夜福利视频在线观看| 在线国产一区二区在线| 精品国产国语对白av| 免费在线观看日本一区| 免费观看精品视频网站| 在线国产一区二区在线| 精品久久久久久成人av| 亚洲色图 男人天堂 中文字幕| 搡老妇女老女人老熟妇| 波多野结衣av一区二区av| 女同久久另类99精品国产91| 精品一区二区三区av网在线观看| 亚洲色图综合在线观看| 丁香六月欧美| 亚洲午夜理论影院| 9191精品国产免费久久| 真人做人爱边吃奶动态| 成人欧美大片| 欧美最黄视频在线播放免费| 久久青草综合色| 啦啦啦韩国在线观看视频| 中文字幕高清在线视频| 午夜日韩欧美国产| 成人欧美大片| 亚洲色图综合在线观看| 99国产精品99久久久久| 色综合婷婷激情| 精品高清国产在线一区| 成人精品一区二区免费| 久久国产亚洲av麻豆专区| 夜夜爽天天搞| 在线观看日韩欧美| 视频区欧美日本亚洲| 一区二区日韩欧美中文字幕| 亚洲一区二区三区不卡视频| 丰满人妻熟妇乱又伦精品不卡| 日本撒尿小便嘘嘘汇集6| bbb黄色大片| 亚洲欧美激情综合另类| 首页视频小说图片口味搜索| 成熟少妇高潮喷水视频| 在线播放国产精品三级| 久久人妻av系列| 国产高清videossex| 午夜久久久久精精品| 精品日产1卡2卡| 国产精品亚洲av一区麻豆| 日韩一卡2卡3卡4卡2021年| 亚洲五月天丁香| av超薄肉色丝袜交足视频| cao死你这个sao货| 欧美日韩乱码在线| 国产精品精品国产色婷婷| 琪琪午夜伦伦电影理论片6080| 久久久久久大精品| 69av精品久久久久久| 男女下面插进去视频免费观看| 男女之事视频高清在线观看| 亚洲一码二码三码区别大吗| 日韩大码丰满熟妇| 给我免费播放毛片高清在线观看| 国产熟女午夜一区二区三区| 女警被强在线播放| 久久久国产成人精品二区| 校园春色视频在线观看| 在线播放国产精品三级| 国产精品一区二区免费欧美| 国产精品一区二区在线不卡| 久久影院123| 一级毛片女人18水好多| 黄色毛片三级朝国网站| 色播在线永久视频| 女性被躁到高潮视频| av免费在线观看网站| 成人av一区二区三区在线看| 国产精品影院久久| 婷婷六月久久综合丁香| 久久中文看片网| 一二三四社区在线视频社区8| 国产午夜精品久久久久久| 国产一区二区三区在线臀色熟女| 亚洲黑人精品在线| 男人操女人黄网站| 免费不卡黄色视频| 欧美另类亚洲清纯唯美| 国产av一区在线观看免费| 久久人人爽av亚洲精品天堂| 国产成人免费无遮挡视频| 亚洲国产欧美网| 亚洲午夜理论影院| 9热在线视频观看99| 热re99久久国产66热| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品久久成人aⅴ小说| 两个人看的免费小视频| 色综合婷婷激情| 777久久人妻少妇嫩草av网站| 久热爱精品视频在线9| 国产一区二区三区视频了| 免费在线观看影片大全网站| 国产三级在线视频| 亚洲激情在线av| 神马国产精品三级电影在线观看 | 国产成人精品久久二区二区91| 亚洲精品中文字幕一二三四区| 男女下面插进去视频免费观看| 亚洲五月婷婷丁香| 亚洲人成伊人成综合网2020| 欧美丝袜亚洲另类 | 大型黄色视频在线免费观看| 精品电影一区二区在线| 日本 av在线| 国产精品免费一区二区三区在线| 亚洲av日韩精品久久久久久密| 最近最新免费中文字幕在线| 久9热在线精品视频| 老司机福利观看| 国产成人av激情在线播放| 国产精品乱码一区二三区的特点 | 亚洲九九香蕉| 亚洲人成电影观看| 久久这里只有精品19| 黄色丝袜av网址大全| 精品久久久精品久久久| 午夜福利高清视频| 日本欧美视频一区| 免费高清视频大片| 亚洲av电影在线进入| 18禁美女被吸乳视频| 禁无遮挡网站| 国产麻豆69| 国产在线精品亚洲第一网站| 午夜日韩欧美国产| 日韩国内少妇激情av| 欧美亚洲日本最大视频资源| 欧美av亚洲av综合av国产av| 欧美激情 高清一区二区三区| 久久久久久久午夜电影| 国产成人精品久久二区二区91| 好男人在线观看高清免费视频 | 非洲黑人性xxxx精品又粗又长| 亚洲国产精品久久男人天堂| 人妻丰满熟妇av一区二区三区| 午夜福利影视在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 麻豆一二三区av精品| 一边摸一边抽搐一进一出视频| 久久久久久免费高清国产稀缺| 国产三级在线视频| 变态另类丝袜制服| 黄片大片在线免费观看| 久久久久国产一级毛片高清牌| 一级毛片精品| 亚洲欧美激情在线| 最好的美女福利视频网| 免费av毛片视频| 久久久久久久午夜电影| 大香蕉久久成人网| 日本在线视频免费播放| 国产亚洲精品av在线| 精品国产超薄肉色丝袜足j| 天天躁狠狠躁夜夜躁狠狠躁| 国产午夜精品久久久久久| 免费看十八禁软件| 久久草成人影院| 91精品国产国语对白视频| 亚洲av五月六月丁香网| 女人高潮潮喷娇喘18禁视频| 亚洲狠狠婷婷综合久久图片| 国产极品粉嫩免费观看在线| 正在播放国产对白刺激| 一边摸一边抽搐一进一出视频| 乱人伦中国视频| АⅤ资源中文在线天堂| 人成视频在线观看免费观看| 精品熟女少妇八av免费久了| 伦理电影免费视频| 久久 成人 亚洲| 精品国内亚洲2022精品成人| 午夜福利影视在线免费观看| 大型黄色视频在线免费观看| 宅男免费午夜| 欧美在线一区亚洲| 桃色一区二区三区在线观看| 久久人人爽av亚洲精品天堂| 日韩国内少妇激情av| 精品电影一区二区在线| 丰满的人妻完整版| 国产成人精品久久二区二区91| 国产97色在线日韩免费| 99精品在免费线老司机午夜| 中文亚洲av片在线观看爽| 国产精品香港三级国产av潘金莲| 国产av精品麻豆| 亚洲色图 男人天堂 中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品av在线| 久久人人97超碰香蕉20202| 国产精品久久久久久精品电影 | 国产精品一区二区三区四区久久 | 97人妻精品一区二区三区麻豆 | 亚洲av五月六月丁香网| 自线自在国产av| 国产av一区二区精品久久| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 精品久久久久久成人av| 韩国av一区二区三区四区| 亚洲伊人色综图| 久久久水蜜桃国产精品网| 18禁裸乳无遮挡免费网站照片 | 精品久久久久久久毛片微露脸| 久久香蕉激情| 91老司机精品| 亚洲一区中文字幕在线| 精品久久久久久,| 亚洲精品美女久久久久99蜜臀| 亚洲国产欧美一区二区综合| 国产精品亚洲av一区麻豆| 成人精品一区二区免费| 97碰自拍视频| 久久久久久大精品| 国内精品久久久久久久电影| 天堂影院成人在线观看| 一二三四在线观看免费中文在| 欧美av亚洲av综合av国产av| 黑丝袜美女国产一区| 亚洲成a人片在线一区二区| 丁香六月欧美| www.自偷自拍.com| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一小说| 级片在线观看| 中文亚洲av片在线观看爽| 国产精品久久久久久精品电影 | 亚洲七黄色美女视频| 十分钟在线观看高清视频www| 韩国av一区二区三区四区| 国产精品影院久久| 久久久久久国产a免费观看| 国产极品粉嫩免费观看在线| 脱女人内裤的视频| 亚洲成人久久性| 久久精品影院6| 18禁黄网站禁片午夜丰满| 久久 成人 亚洲| 大型黄色视频在线免费观看| 色综合站精品国产| 香蕉国产在线看| 午夜免费鲁丝| 亚洲三区欧美一区| 99热只有精品国产| 欧美成人性av电影在线观看| 夜夜夜夜夜久久久久| 亚洲色图综合在线观看| 国产午夜福利久久久久久| 一级a爱片免费观看的视频| 国产精品精品国产色婷婷| 操出白浆在线播放| 窝窝影院91人妻| 纯流量卡能插随身wifi吗| АⅤ资源中文在线天堂| 色播亚洲综合网| 精品一区二区三区四区五区乱码| 日韩大尺度精品在线看网址 | 91成年电影在线观看| 大香蕉久久成人网| 国产aⅴ精品一区二区三区波| av免费在线观看网站| 久久天堂一区二区三区四区| 99国产精品一区二区蜜桃av| 成人三级做爰电影| 超碰成人久久| 99国产极品粉嫩在线观看| 欧美成人性av电影在线观看| 老司机福利观看| 成人永久免费在线观看视频| 波多野结衣av一区二区av| 亚洲精品中文字幕一二三四区| 狂野欧美激情性xxxx| 精品久久久精品久久久| 亚洲精品一区av在线观看| 久久人妻福利社区极品人妻图片| 国产极品粉嫩免费观看在线| 波多野结衣巨乳人妻| 国产野战对白在线观看| 51午夜福利影视在线观看| 午夜久久久久精精品| 香蕉久久夜色| 在线视频色国产色| 色播亚洲综合网| 最新美女视频免费是黄的| 欧美成人午夜精品| 在线观看一区二区三区| 日本 欧美在线| 精品欧美一区二区三区在线| 亚洲色图av天堂| 宅男免费午夜| 麻豆成人av在线观看| 久久欧美精品欧美久久欧美| 丝袜美腿诱惑在线| 欧美另类亚洲清纯唯美| 国产精品国产高清国产av| 日韩欧美一区二区三区在线观看| www.www免费av| 男男h啪啪无遮挡| 久久国产精品男人的天堂亚洲| 99riav亚洲国产免费| 亚洲av成人一区二区三| 美女免费视频网站| tocl精华| 精品国产国语对白av| 国产精品香港三级国产av潘金莲| 精品欧美国产一区二区三| 热re99久久国产66热| 男女午夜视频在线观看| 国产高清有码在线观看视频 | 国产精品98久久久久久宅男小说| 午夜福利一区二区在线看| 国产av又大| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费观看网址| 精品一区二区三区四区五区乱码| 欧美性长视频在线观看| 99热只有精品国产| 在线av久久热| 久久精品亚洲熟妇少妇任你| 1024香蕉在线观看| 亚洲国产精品成人综合色| 成人亚洲精品av一区二区| 免费无遮挡裸体视频| 国产精品久久久人人做人人爽| 国产av又大| 熟妇人妻久久中文字幕3abv| 国产精品,欧美在线| 给我免费播放毛片高清在线观看| 国产成人精品久久二区二区免费| 午夜精品国产一区二区电影| 久久久精品欧美日韩精品| 国产激情久久老熟女| 欧美一区二区精品小视频在线| 免费看美女性在线毛片视频| 亚洲国产中文字幕在线视频| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 亚洲中文字幕一区二区三区有码在线看 | 又黄又粗又硬又大视频| 美女免费视频网站| 欧美乱妇无乱码| 色播在线永久视频| 亚洲激情在线av| 69精品国产乱码久久久| 搞女人的毛片| 欧美激情 高清一区二区三区| 很黄的视频免费| 亚洲自偷自拍图片 自拍| 999久久久国产精品视频| 亚洲,欧美精品.| 一二三四社区在线视频社区8| 欧美av亚洲av综合av国产av| 亚洲成人久久性| 在线观看免费视频日本深夜| 亚洲精品中文字幕在线视频| 人人妻,人人澡人人爽秒播| 亚洲熟妇熟女久久| 99精品久久久久人妻精品| 一边摸一边抽搐一进一出视频| 日本精品一区二区三区蜜桃| 一进一出抽搐动态| 欧美另类亚洲清纯唯美| 色综合婷婷激情| 精品一品国产午夜福利视频| 亚洲人成电影免费在线| 侵犯人妻中文字幕一二三四区| 国产成人精品在线电影| 久久久久精品国产欧美久久久| 国产三级在线视频| 制服诱惑二区| 欧美精品啪啪一区二区三区| av有码第一页| 中文字幕最新亚洲高清| 亚洲美女黄片视频| 国产色视频综合| 欧美性长视频在线观看| 国产午夜精品久久久久久| 正在播放国产对白刺激| 啦啦啦免费观看视频1| 免费在线观看黄色视频的| 9热在线视频观看99| 黑丝袜美女国产一区| 午夜免费鲁丝| av中文乱码字幕在线| 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 亚洲中文av在线| 两个人看的免费小视频| 亚洲国产毛片av蜜桃av| 欧美日韩一级在线毛片| 亚洲精品一卡2卡三卡4卡5卡| or卡值多少钱| 我的亚洲天堂| 欧美在线黄色| 欧美色欧美亚洲另类二区 | 国产精品 国内视频| 韩国av一区二区三区四区| 天堂动漫精品| svipshipincom国产片| 人人妻,人人澡人人爽秒播| 亚洲黑人精品在线| 亚洲专区字幕在线| 久久香蕉国产精品| 国产激情久久老熟女| 午夜福利18| 久久久国产成人免费| 亚洲最大成人中文| 亚洲精品粉嫩美女一区| 亚洲avbb在线观看| 亚洲成人精品中文字幕电影| 麻豆国产av国片精品| 国产精品亚洲一级av第二区| 黄频高清免费视频| 波多野结衣巨乳人妻| 亚洲精品美女久久久久99蜜臀| 九色国产91popny在线| 国产又色又爽无遮挡免费看| 日本一区二区免费在线视频| 久久国产精品影院| 欧美日韩乱码在线| 亚洲国产欧美网| 91精品三级在线观看| 9色porny在线观看| 久热这里只有精品99| 在线av久久热| 97碰自拍视频| 99riav亚洲国产免费| 免费无遮挡裸体视频| 欧美激情高清一区二区三区| 久久青草综合色| 国产欧美日韩一区二区三| 女人精品久久久久毛片| 免费在线观看影片大全网站| 母亲3免费完整高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 90打野战视频偷拍视频| 日日夜夜操网爽| 亚洲久久久国产精品| 欧美大码av| 丰满的人妻完整版| 欧美另类亚洲清纯唯美| 色在线成人网| 国产麻豆成人av免费视频| 麻豆久久精品国产亚洲av| 18禁裸乳无遮挡免费网站照片 | 日本一区二区免费在线视频| 亚洲精品国产区一区二| 久久久久国内视频| 精品卡一卡二卡四卡免费| 国产一区二区激情短视频| 人妻丰满熟妇av一区二区三区| 91麻豆精品激情在线观看国产| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 亚洲第一欧美日韩一区二区三区| 午夜影院日韩av| 操出白浆在线播放| 亚洲精品中文字幕在线视频| 99riav亚洲国产免费| 亚洲av五月六月丁香网| 国产精品精品国产色婷婷| 亚洲色图综合在线观看| 岛国视频午夜一区免费看| www.熟女人妻精品国产| 国产精品久久久久久精品电影 | 国产欧美日韩综合在线一区二区| 丝袜在线中文字幕| 男人操女人黄网站| 欧美成人午夜精品| 一区二区三区精品91| 在线观看www视频免费| 成人手机av| 激情视频va一区二区三区| 桃红色精品国产亚洲av| 亚洲熟女毛片儿| 午夜精品在线福利| 999精品在线视频| 一区二区三区精品91| 999久久久精品免费观看国产| 两性夫妻黄色片| 999精品在线视频| 久99久视频精品免费| 精品国产超薄肉色丝袜足j| 成人手机av| 亚洲国产欧美日韩在线播放| 在线观看免费午夜福利视频| 大型黄色视频在线免费观看| 欧美日韩一级在线毛片| 自线自在国产av| 别揉我奶头~嗯~啊~动态视频| 丝袜在线中文字幕| 国产成人av教育| 国产av在哪里看| 久久这里只有精品19| 日韩精品青青久久久久久| 老司机午夜十八禁免费视频| 亚洲自拍偷在线| 精品欧美一区二区三区在线| 丁香欧美五月| 日韩一卡2卡3卡4卡2021年| 人人妻人人澡欧美一区二区 | 神马国产精品三级电影在线观看 | av在线播放免费不卡| 亚洲男人的天堂狠狠| 日韩欧美一区二区三区在线观看| 国产精品久久电影中文字幕| 国产精品久久久人人做人人爽| 久久人人爽av亚洲精品天堂| 亚洲av成人av| 国产激情欧美一区二区| 色综合站精品国产| 咕卡用的链子| 欧美日韩乱码在线| 9热在线视频观看99| 久久久久九九精品影院| 在线观看免费视频日本深夜| 如日韩欧美国产精品一区二区三区| 悠悠久久av| 非洲黑人性xxxx精品又粗又长| 国产三级在线视频| 精品日产1卡2卡| 18禁国产床啪视频网站| 黄色a级毛片大全视频| 亚洲熟妇熟女久久| 91九色精品人成在线观看| 丝袜美腿诱惑在线| 欧美激情久久久久久爽电影 | 又黄又爽又免费观看的视频| 免费看十八禁软件| 老汉色∧v一级毛片| 操出白浆在线播放| 久久人人精品亚洲av| 首页视频小说图片口味搜索| 欧美在线黄色| 高清在线国产一区|