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

    Do European Blocking Events Precede North Atlantic Oscillation Events?

    2015-06-09 21:37:22YAOYaoandLUODehai
    Advances in Atmospheric Sciences 2015年8期

    YAO Yao and LUO Dehai

    Key Laboratory of Regional Climate-Environment for Temperate East Asia, Institute of Atmospheric Physics,Chinese Academy of Science,Beijing 100029

    Do European Blocking Events Precede North Atlantic Oscillation Events?

    YAO Yao and LUO Dehai?

    Key Laboratory of Regional Climate-Environment for Temperate East Asia, Institute of Atmospheric Physics,Chinese Academy of Science,Beijing 100029

    Using a two-dimensional blocking index,the cause and effect relationship between European blocking(EB)events and North Atlantic Oscillation(NAO)events is investigated.It is shown that the EB event frequency is enhanced over Northern(Southern)Europe for negative(positive)phases of the NAO.Enhanced EB events over Northern Europe precede the establishment of negative phase NAO(NAO-)events,while the enhanced frequency of EB events over Southern Europe lags positive phase NAO(NAO+)events.

    The physical explanation for why enhanced EB events over Northern(Southern)Europe lead(lag)NAO-(NAO+) events is also provided.It is found that the lead–lag relationship between EB events in different regions and the phase of NAO events can be explained in terms of the different latitudinal distribution of zonal wind associated with the different phases of NAO events.For NAO+events,the self-maintained eastward displacement of intensified midlatitude positive height anomalies owing to the intensified zonal wind can enhance the frequency of EB events over Southern Europe,thus supporting a standpoint that EB events over Southern Europe lag NAO+events.Over Northern Europe,EB events lead NAO-events because NAO-events arise from the self-maintained westward migration of intensified blocking anticyclones due to the weakened zonal wind in higher latitudes.

    European blocking,North Atlantic Oscillation,mutual relationship

    1.Introduction

    The North Atlantic Oscillation(NAO)is an important low frequency mode of climate variability in the Northern Hemisphere(NH),identi fied as a north–south dipole pattern confined to the Atlantic basin(Hurrell,1995).The basic dynamics of NAO events has been an important research direction (Feldstein,2003;Luo et al.,2007a,2007b;Luo and Cha, 2012)because of its importance in climate changes in the NH (Scaife et al.,2008;Sillmann et al.,2011).

    Many studies have demonstrated that negative(positive) phase NAO events,hereafter referred to as NAO-(NAO+) events,arise from the cyclonic(anticyclonic)wave breaking of synoptic-scale waves(Benedict et al.,2004;Franzke et al.,2004;Rivi`ere and Orlanski,2007;Woollings et al., 2008;Strong and Magnusdottir,2008),which are also suggested to correspond to the presence(absence)of blocked days over the North Atlantic(Luo et al.,2007a;Woollings et al.,2008).However,based upon a nonlinear multi-scale interaction model of NAO events developed by Luo et al. (2007a,2007b),Luo et al.(2008)suggested that synoptic wave breaking is not a cause of NAO onset,but instead is a concomitant result of the eddy-induced intensified NAO pattern.Other studies have established links between the phase of the NAO and blocking activity(Shabbar et al.,2001;Luo, 2005;Scherrer etal.,2006;Croci-Maspoli etal.,2007;Luo et al.,2007a,2014).Among these studies,an interesting result is that the negative(positive)phase of the NAO corresponds to an enhanced blocking frequency over Greenland(Europe). In terms of a wave breaking index,Woollings et al.(2008) suggested that NAO-(NAO+)events may arise from the succession of blocked(non-blocked)days over Greenland.Feldstein(2003)and Sung et al.(2011)attributed NAO-events to the existence of a blocking-type anomaly over Northern Europe.More recently,Davini et al.(2012b)found that the reduction of the Greenland blocking frequency leads to an eastward shift of the NAO pattern.Although these studies have revealed a close link between European blocking(EB) events and the phase of NAO events,the cause and effect relationship between EB events and NAO events remains unclear. This is a critical problem that we attempt to investigate in the present paper.

    The remainder of the paper is organized as follows:In section 2,we describe the data and blocking index used in the analysis.Section 3 depicts the geographical distributions of the winter-mean blocking frequency,as well as the instantaneous blocking(IB)frequencies associated with in-situ NAOevents and transition events.We examine the cause and effect relationship between EB events and the phase of the NAO in section 4.A physical explanation of the lead–lag relationship between EB events and the phase of the NAO is presented in section 5.A self-maintaining mechanism of the westward (eastward)migration of the NAO-(NAO+)dipole pattern is also proposed in this section.Finally,the key conclusions are summarized,and further discussion is provided,in section 6.

    2.Data and methodology

    2.1.Data description and definition of an NAO event

    This study uses the daily multi-level geopotential height and wind fields from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR)reanalysis at a horizontal resolution of 2.5?× 2.5?from 1950 to 2011.In this paper,we focus on winter, defined as the time interval from November to March.The daily anomaly,for which the seasonal cycle is removed,is defined as the deviation of the daily field from its climatological time mean.An NAO+(NAO-)event is defined if the normalized daily NAO index is greater(less)than or equal to +1.0(-1.0)standard deviation persisting for at least three consecutive days.The daily NAO index,based on the Rotated Empirical Orthogonal Function(REOF)of daily geopotential height fields at 500 hPa,is available from the National Oceanic and Atmospheric Administration/Climate Prediction Center(Kalnay et al.,1996).

    Following Luo et al.(2012a),NAO events are categorized into two types:in-situ(individual)NAO events and transition events.The in-situ NAO events are defined as those events with only one phase,either NAO+or NAO-.The transition events are defined as those events of one phase followed by the opposite phase,and must include both NAO+and NAO-events.A detailed definition of NAO transition events can be found in Luo et al.(2012a).

    Based on the definitions of these events,we detect 99 insitu NAO-events and 97 in-situ NAO+events,as well as 38 NAO+to NAO-and 30 NAO-to NAO+transition events in winter during the period from 1950 to 2011.In this paper, lag(0)represents the day with maximum amplitude during the NAO(NAO+and NAO-)life cycle.For transition NAO events,lag(0)represents the day with maximum amplitude during the NAO+of NAO+to NAO-transition events and the NAO-of NAO-to NAO+transition events,respectively.

    2.2.Two-dimensional blocking index

    In recent years,many two-dimensional(2D)blocking indices have been proposed to identify the spatial and temporal variations of blocking events(Barriopedro et al.,2006,2010; Diao et al.,2006;Scherrer et al.,2006;Davini et al.,2012a, 2012b)based upon the extension of one-dimensional blocking indices(Lejen¨as and ?kland,1983;Tibaldi and Molteni, 1990;Pelly and Hoskins,2003).The recent 2D blocking index developed by Davini et al.(2012a)has the advantage of differentiating high-latitude and low-latitude events of blocking flows in the Euro-Atlantic sector compared to other 2D indices(Diao et al.,2006).Thus,it is appropriate to use the 2D blocking index of Davini et al.(2012a,2012b)to perform our investigation because our focus is on the spatiotemporal variation of Euro-Atlantic blocking events.To detect IB events in the Euro-Atlantic region,analogous to Davini et al.(2012a)we define the meridional gradient reversal of the 500-hPa geopotential height at each grid point as

    where Z(λ0,φ0)is the daily 500-hPa geopotential height at the grid point(λ0,φ0),λ0(φ0)is the longitude(latitude)ranging from 0?to 360?W(30?N to 75?N),φN=φ0+15,and φS=φ0-15(Davinietal.,2012a,2012b).GHGNand GHGS represent the geopotential height gradients[units:m(?lat)-1] at 500 hPa over the north and south sides of(λ0,φ0).

    One can identify an IB event if conditions(1–4)are satisfied.Here,very low-latitude blocking events such as subtropicalridges in the Atlantic basin are excluded through constraint(4).According to the definition of a blocking event (Davini et al.,2012a,2012b),an IB event satisfies conditions (1–4).However,a blocking event is a continuous process of an IB event persisting for at least 5 days at a given grid point and has at least a large spatial scale[such as 5?(lat)×10?(lon),as described in Davini et al.(2012a)].Because the frequency of IB events exhibits a striking resemblance to that of blocking events,the IB frequency in the Euro-Atlantic sector is only calculated and used as a measure of the blocking occurrence in this paper.The sensitivity of this 2D blocking index defined in Eq.(1)to different data is examined in Davini et al.(2012a).

    3.Geographical distribution of Euro-Atlantic blocking events and its link with the phase of NAO events

    To understand the possible impact of NAO events on the spatial distribution of blocking events in the Euro-Atlantic sector,it is reasonable to use the 2D blocking index of Davini et al.(2012a,2012b)to calculate the IB frequency.Figure 1 shows the geographical distributions of the winter mean IB frequency as well as the spatial distributions of the composite IB frequencies associated with in-situ NAO and transition events during 1950–2011.Here,the winter mean IB frequency reflects the IB climatology in winter,which is measured as the percentage of blocked days with respect to the total days of winter during 1950–2011.The calculation is applied during the life cycle of an NAO event[from lag(-10) to lag(+10)for in-situ NAO events,but from lag(-10)to lag(+5)and lag(+5)to lag(+25)for NAO transition events]. It is evident that there are two frequency centers of blocking events in the Euro-Atlantic sector:Greenland and Central Europe(Fig.1a).However,the occurrence region of Euro-Atlantic blocking events is modulated by the variation of the NAO phase.For the positive phase of in-situ NAO events,the higher IB frequency is concentrated in one region from the central-eastern Atlantic to Southern Europe(Fig.1b),while the IB frequency for the negative phase has two maximum regions:Greenland and Northern Europe(Fig.1c).Thus,the action center of Euro-Atlantic blocking events is affected by the phase of in-situ NAO events.

    If an NAO event undergoes a transition from one phase to another,the geographical distribution of the IB frequency can exhibit a marked change from a southwest–northeast (southeast–northwest)pattern distribution to another in the southeast–northwest(southwest–northeast)direction(Figs. 1e–g).The results are evidently different from the previous results of Scherrer et al.(2006)and Croci-Maspoli et al. (2007),who only noted that the blocking frequency is en-hanced(reduced)over Europe for positive(negative)phases. However,here we detect a new result that there is an enhanced blocking frequency over Southern(Northern)Europe for positive(negative)phases.Thus,the phase of NAO events seems to play a crucial role in controlling the spatial distribution of EB events.

    As can be seen from the data in Table 1,the numbers of transition NAO events are smaller than those of in-situ NAO events.Statistical analysis indicates that the total number of days of in-situ NAO events accounts for about 60%of the total number of NAO event days.Actually,transition NAO events are a combination oftwo phases of in-situ NAOevents. The impact of NAO+(NAO-)episodes of NAO+to NAO-transition events on IB frequency is similar to in-situ NAO+(NAO-)events,and vice versa for NAO-to NAO+transition events.The regime transition of NAO events plays an important role in the interannual variability of the NAO(Luo et al., 2012a,2012b),but this is beyond the scope of the current paper.Thus,in this study,we focus only on analyzing in-situ NAO events.

    4.Cause and effect relationship between EB events and NAO events

    4.1.IB frequency indices and their relationship with NAO events

    To describe the longitudinal and latitudinal variations of Euro-Atlantic blocking events,the regions(60?–77.5?N, 10?W–30?E)and(35?–50?N,10?W–30?E)are defined as Northern Europe and Southern Europe,respectively.Correspondingly,the regions(60?–77.5?N,60?–10?W)and(35?–50?N,60?–10?W)are defined as the high-latitude North Atlantic and low-latitude North Atlantic,respectively.By averaging the IB frequencies in each grid within these four regions,we can construct time series of the winter IB frequency index in each region during 1950–2011 and show the normalized time series of the winter IB frequencies in the four regions in Fig.2.The IB frequency index in each region is defined to be high(low)if the normalized index of the IB frequency is greater(less)than or equal to+1.0(-1.0)the standard deviation.Based on the number of winters that correspond to high and low IB frequencies,the numbers of insitu NAO+,in-situ NAO-,NAO+to NAO-transition events, and NAO+to NAO-transition events can be calculated,the results of which are shown in Table 1.It can be seen that insitu NAO-(NAO+)events are more(less)frequent for high IB frequency winters compared to low IB frequency winters over Northern Europe.In Southern Europe,the frequency of in-situ NAO-(NAO+)events is higher for low(high)IB frequency winters compared to high(low)IB frequency winters. Thissuggests thatfrequent in-situ NAO-(NAO+)eventscorrespond to high(low)IB frequency over Northern Europe or low(high)IB frequency over Southern Europe.Similarly,the high frequency of in-situ NAO-(NAO+)events is seen for high(low)IB frequency winters in the high-latitude North Atlantic.In contrast,the in-situ NAO-(NAO+)events have a high frequency for low(high)IB frequency winters in the low-latitude North Atlantic.Thus,the difference in the blocking event frequency in different regions is closely related tothe phase of in-situ NAO events.

    Table 1.Event and year numbers of in-situ NAO events and transition events for high and low IB frequencies over Northern Europe, Southern Europe,the high-latitude North Atlantic,and low-latitude North Atlantic.Bold numbers represent the dominant events.

    4.2.Empirical Orthogonal Function patterns of the monthlymeangeopotentialheightanomaliesforhigh IB frequency winters in different regions

    In this subsection,Empirical Orthogonal Function(EOF) analysis is used to identify the leading modes of high IB index winters for different regions.

    Forhigh IBfrequency winters in Northern Europe,Southern Europe,the high-latitude North Atlantic,and low-latitude North Atlantic,we show the leading and second EOF components of monthly mean 500-hPa geopotential height fields and their corresponding principal component time series in Figs.3 and 4,respectively.It is clear that when the IB frequency in Northern Europe is higher(Fig.3a),the leading EOF(EOF1),which explains 26.3%of the total variance, shows a northwest–southeast oriented NAO-pattern,while the second EOF(EOF2),which explains 23.7%of the total variance,shows a Scandinavian blocking pattern(Cassou et al.,2004).However,the EOF1 and EOF2 correspond to NAO+and negative phase Eastern-Atlantic(EA-)patterns (Moore and Renfrew,2012)respectively forlow IBfrequency winters(less than or equal to-1.0 standard deviation for the IB frequency index)in Northern Europe(not shown).For a high IB frequency in Southern Europe,the EOF1(EOF2)exhibits an NAO+(EA-)pattern(Fig.3b).The EOF1 and EOF2 for the low IB frequency over Southern Europe are similar to those for high IB frequency over Northern Europe (not shown).For a high IB frequency in the high-latitude North Atlantic,the EOF1 and EOF2 describe the northwest–southeast and northeast–southwest oriented NAO-patterns (Fig.4a),in which the EOF2 pattern seemsto have a more evident eastward migration.However,for low IB frequency in the high-latitude North Atlantic,the EOF1 and EOF2 exhibit the NAO+and EA-patterns(not shown).This result is also seen for high IB frequency in the low-latitude North Atlantic (Fig.4b).When the IB frequency is low in the low-latitude North Atlantic,the EOF1 and EOF2 describe the northwest–southeast oriented NAO-pattern and Scandanavian blocking pattern,respectively(not shown),which is also detected for high IB frequency over Northern Europe and low IB frequency over Southern Europe.Thus,this suggests a strong linkage between changes in blocking frequencies in different regions and variations of monthly mean EOF patterns.Therefore,the connection between the phase of the NAO and the variation of blocking frequency can be seen in terms of the EOF analysis of the monthly mean height field data.This motivates us to further examine the lead–lag relationship between EB and NAO events from the daily data in the next subsection.

    4.3.Lead–lag relationship between the composite daily NAO and IB frequency indices in different regions

    To further reveal the lead–lag relationship between EB events and NAO events,we construct composites of the daily NAO indices and IB frequencies in different regions(Fig. 5).It can be seen that the daily NAO-index has a 2-day negative correlation delay of-0.89 with the IB frequency over Northern Europe(Fig.5a)from lag(-10)to lag(+10), which is statistically significant at the 99.9%confidence level (Glahn,1968).Further analysis indicates that 21(80.8%)of 26 in-situ NAO-events are lagged by EB events.Thus,EB events lead the onset of NAO-events.This result was also noted by Schwierz et al.(2004)and Sung et al.(2011).Sung et al.(2011)found that blocking ridges over Northern Europe contribute to the growth of NAO-events through a lowfrequency eddy interaction.Nevertheless,the IB frequency over Southern Europe tends to have a 3-day positive correlation delay of 0.81 with the composite daily NAO+index(Fig. 5b).This correlation is significant at the 99.9%confidence level.In addition,29(87.9%)of 33 in-situ NAO+events are found to precede EB events.Hence,EB events over Southern Europe lag NAO+Events.Thus,it seems reasonable to attribute an enhanced EB event frequency over Southern Eu-rope to the genesis of NAO+events.However,a noticeably different characteristic is found in the high-latitude North Atlantic and low-latitude North Atlantic,where the daily NAO index and IB frequency exhibit an in-phase correlation(Figs. 5c and d).The correlation coefficient is particularly high for the high IB frequencies over the high-latitude North Atlantic and low-latitude North Atlantic.For example,during the NAO-(NAO+)life cycle the coefficient of their negative(positive)correlation is-0.94(0.92).As a result,the relationship between the phase of the NAO and blocking events is different between the Atlantic basin and continental Europe. Thus,we can conclude from the results presented here that EB events over Northern(Southern)Europe may lead(lag) the genesis of NAO-(NAO+)events.The results are new, although the result of NAO-events lagging blocking events over Northern Europe was also noted by Sung et al.(2011).

    5.Physical mechanism of the lead–lag relationship between EB events and NAOevents

    5.1.Role of initial large-scale anomalies in the growth of the two phases of the NAO

    To reveal the physical mechanism of the lead(lag)correlation between EB events and NAO-(NAO+)events,we show the composite maps of daily 500-hPa geopotential height anomalies of in-situ NAO-(NAO+)events for high IB frequency over Northern(Southern)Europe from lag(-6) to lag(4)(during the life cycle of the NAO)in Fig.6. As can be seen from Fig.6a,an anticyclonic anomaly appears over Northern Europe at lag(-6).The anticyclonic anomaly grows owing to the high-frequency eddy forcing Luo et al.(2007a)and the cyclonic wave breaking(Benedict et al.,2004;Franzke et al.,2004;Rivi`ere and Orlanski,2007; Woollings et al.,2008).The amplified anticyclonic anomaly over Northern Europe undergoes a westward migration and extends toward the western North Atlantic,and then establishes an NAO-dipole pattern[at lag(0)]through combining with an amplified cyclonic anomaly to its south.Thus, blocking ridges over Northern Europe prior to the negative phase of the NAO may be considered as a precursor to NAO-events.Feldstein(2003)attributed the maintenance of the NAO anomaly to high-frequency eddy forcing,but the NAO decay to Ekmann dissipation;whereas,Luo et al.(2007a) suggested that the growth and decay of the NAO anomaly is due to positive synoptic-scale eddy forcing and the reversal of its sign.Sung et al.(2011)noted that NAO-events frequently accompany blocking ridges over the Barents Sea and are primarily forced by the strong positive low-frequency eddy forcing related to blocking ridgesoverNorthern Europe.

    In Fig.6b we can see that a cyclonic anomaly appears over the western North Atlantic at lag(-6).Thus,the cyclonic anomaly in the western North Atlantic is a precondition of the onset of NAO+events.The initial growth of the cyclonic anomaly results from the negative low-frequency eddy forcing over the western North Atlantic and Greenland (Sung et al.,2011).However,its subsequent growth(decay) during the entire life cycle is more likely to be due to highfrequency eddy forcing(Luo et al.,2007a)or due to the presence(absence)of the anticyclonic wave breaking of synopticscale eddies(Benedict et al.,2004;Woollings et al.,2008). Moreover,the anticyclonic anomaly in the low-latitude North Atlantic intensifies with the growth of the cyclonic anomaly in the high-latitude North Atlantic.This process will cause the establishment of an NAO+event.During the growing process,the NAO+dipole pattern can undergo eastward migration and extension toward continental Europe,which is consistent with the result of Luo and Cha(2012).Because the anticyclonic anomaly of the NAO+dipole pattern occupies a broader region in the midlatitudes,it can even extend into the Urals region during the mature stage of an NAO+event from lag(-1)to lag(+1)(not shown).In this case,it is inevitable that we see an enhanced frequency of EB events (in the form of an anticyclonic anomaly)over Southern Europe during NAO+growth.Consequently,itis concluded that EB events over Southern Europe lag NAO+events.

    It should be pointed out that EB events in Northern Europe in fact correspond to an anticyclonic anomaly over the region.When an EB event undergoes westward migration and extends toward the North Atlantic,it becomes an NAO-event.The observational study of Sung et al.(2011)indicates thatmostNAO-events resultfrom the westward migration of amplifying Northern European anticyclones.Thus,to some extent,the anticyclonic anomaly over Northern Europe(the North Atlantic)is the flow reflection of an EB(NAO-)event.

    5.2.Longitudinalmigrationofgeopotentialheight anomalies in different latitudinal bands during the life cycle of the NAO for its two phases

    To reflect the zonal migration of the geopotential height anomalies(especially positive anomalies)fordifferentphases of the NAO,Hovm¨oller diagrams of the composite daily geopotential height anomalies averaged over 60?–77.5?N (35?–50?N)for the negative(positive)phase are shown in Fig.7 for high IB frequencies over Northern Europe,Southern Europe,the high-latitude North Atlantic,and the lowlatitude North Atlantic.It is clear that for the negative phase of the NAO the positive anomaly over Northern Europe undergoes a noticeable westward shift after lag(-6)(Fig.7a), and continues its westward migration until lag(+6).The positive anomaly tends to reach its strongest amplitude near 40?W at lag(0).Such a westward migration is also seen in the high-latitude North Atlantic(Fig.7c),where only the distinct westward migration is detected before lag(0),and the strongest amplitude occurs at lag(0)and locates near 40?W.Although the westward shift of positive anomalies is somewhat weaker in the high-latitude North Atlantic than in Southern Europe,it is concluded that the genesis of NAO-events results from the westward migration of amplified anticyclonic anomalies that initiate over Northern Europe.

    On the other hand,for the positive phase of the NAO,it is found that the positive anomaly undergoes an eastward migration not only in the low-latitude North Atlantic(Fig.7d), but also in Southern Europe(Fig.7b).Such an eastwardshift persists during the entire lifespan of an NAO+event.In particular,the amplified positive anomaly can extend more eastward into inland continental Europe,leading to an enhanced EB event frequency over Southern Europe.Thus,it is concluded that EB events over Southern Europe lag NAO+events.The result that EB events over Northern(Southern) Europe lead(lag)NAO-(NAO+)events can be interpreted in terms of a self-maintaining mechanism of the longitudinal migration of the NAO dipole pattern,as proposed in the next subsection.

    5.3.Self-maintaining mechanism of the westward(eastward)migration of NAO-(NAO+)events

    To understand the physical cause of the westward(eastward)migration of positive anomalies associated with NAO-(NAO+)events,it seems reasonable to show the meridional profile of the time-mean zonal winds averaged in the Atlantic basin(60?–10?W)and continental Europe(10?W–30?E)during the entire lifespan of NAO events(Fig.8).It is interesting to see that the zonal wind over Northern Europe(50?–70?N)is relatively weak for the negative phase of the NAO relative to the positive phase(Fig.8b).The amplified positive anomaly inevitably undergoes a westward shift once the positive anomaly that initiates over Northern Europe intensi fies because of the strong eddy forcing.Because of the strengthening(weakening)of zonal wind in the low-latitude, 25?–40?N(high-latitude,45?–65?N)Atlantic basin associated with the growing NAO-anomaly(Fig.8a),the weak zonal wind in the high-latitude North Atlantic tends to render the westward shift of the NAO-dipole pattern.According to Jung et al.(2003),the strength of the mean zonal wind is a major cause of the zonal shift of the NAO dipole mode.Our result also holds if the basic mean zonal wind[de fined as the zonal wind prior to NAO onset(Yao and Luo,2014)]is used (not shown).However,Luo et al.(2010)found from both observed and theoretical aspects that the zonal shift of the NAO pattern is dominated by the meridional variation of the Atlantic jet associated with the NAO phase due to the different phase speed.This is a possible mechanism that can be used here to provide an explanation for the results of the present study.

    For the positive phase of the NAO,the zonal wind intensi fies and shifts northward due to the feedback of the growing NAO+pattern,such that the zonal wind signi ficantly intensifies in the midlatitudes centered at 50?N(Fig.8a).This result has been theoretically noted by Luo et al.(2007b;2010).Thestrengthening of the zonal wind in the low-latitude(0?–20?N) North Atlantic can favor the eastward migration of the anticyclonic anomaly of the NAO+dipole pattern and its further extension into Southern Europe.This leads to an enhanced frequency of EB events over Southern Europe.Hence,it is natural that EB events over Southern Europe lag NAO+events.

    The theoretical results of Luo et al.(2007b)based on weakly nonlinear multi-scale interaction model of NAO events have demonstrated that the strengthening(weakening) of zonal wind in the relatively low(high)latitudes of the Atlantic basin is a spontaneous consequence of NAO-genesis. The reason is that the zonal wind is intensified(weakened) in the middle–high latitudes for the positive(negative)phase of the NAO pattern due to the presence of the intensified NAO-(NAO+)anomaly(Luo et al.,2007b).The weakening (strengthening)ofzonalwind in the mid–high latitudes favors the westward(eastward)migration of the NAO-(NAO+) anomaly(Luo et al.,2010).Such a process is the so-called“self-maintaining”mechanism.

    6.Conclusions and discussion

    This paper examines the cause and effect relationship between EBeventsand the phase ofthe NAO.Itisfound thatthe frequency of EB events is enhanced over Northern(Southern)Europe for the negative(positive)phase of the NAO. This is a new result,different from previous findings reported by Scherrer et al.(2006),Woollings et al.(2008)and Davini et al.(2012b),who only noted an enhanced frequency of EB events for the positive phase of the NAO.Our observational analysis shows that EB events over Northern(Southern)Europe lead(lag)NAO-(NAO+)events.To some extent,the results derived from this study are a complement to the previous studies of Sung et al.(2011).They found that blocking ridges over Northern Europe frequently precede negative NAO events and suggested that the low-frequency eddy forcing associated with the blocking ridges contribute positively to the initial growth of negative NAO events.However,in this paper we further find that NAO+events precede enhanced EB events over Southern Europe.This is a new result,notreported previously.

    A self-maintaining mechanism of the westward(eastward)migration of NAO-(NAO+)events is proposed to explain why EB events over Northern(Southern)Europe lead (lag)NAO-(NAO+)events,based on the theoretical studies of Luo et al.(2007b,2010).This mechanism suggests that the westward(eastward)migration of NAO-(NAO+) events can be attributed to the weakening(strengthening)of the zonal wind in the middle—high latitudes owing to the feedback of the NAO-(NAO+)dipole pattern(Luo et al., 2007b).Thus,it is concluded that EB events over Northern (Southern)Europe leading(lagging)NAO-(NAO+)events during the NAO-(NAO+)life cycle is a self-maintaining phenomenon.

    Although the focus of the present paper is not on the physical mechanism of NAO genesis,the physical mechanism of the onset of NAO events can be better explained in terms of the nonlinear multi-scale interaction model of NAO events(Luo et al.,2007a,2007b).With the weakening(intensifying)of zonal winds in middle—high latitudes during the growth of an NAO-(NAO+)event,the amplified anticyclonic anomaly of the NAO-(NAO+)dipole pattern undergoes a marked westward(eastward)migration away from the source region.Such an opposite zonal migration between the NAO-and NAO+dipole patterns implies that there is a lead (lag)relationship between EB events over Northern(Southern)Europe and NAO-(NAO+)events because the anticyclonic anomaly is often identified as a blocking event.

    However,it must be pointed out that the blocking latitude location in continental Europe is significantly different between the positive and negative phases of the NAO.Whether or not the different latitude location of blocking events over Europe can exert a significant influence on extreme weather events is not addressed in the present paper,but further understanding of this problem is needed.

    Acknowledgements.The authors acknowledge the support from the National Natural Science Foundation of China(Grant Nos. 41375067 and 41430533)and we thank the two anonymous reviewers for their constructive comments for improving this paper.

    REFERENCES

    Barriopedro,D.,R.Garc′?a-Herrera,A.R.Lupo,and E.Hern′andez, 2006:A climatology of Northern Hemisphere blocking.J. Climate,19,1042–1063,doi:10.1175/JCLI3678.1.

    Barriopedro,D.,R.Garc′?a-Herrera,and R.M.Trigo,2010:Application of blocking diagnosis methods to General Circulation Models.Part I:A novel detection scheme.Climate Dyn.,35, 1373–1391,doi:10.1007/s00382-010-0767-5.

    Benedict J.J.,S.Lee,and S.B.Feldstein,2004:Synoptic view of the North Atlantic Oscillation.J.Atmos.Sci.,61,121–144.

    Cassou,C.,L.Terray,J.W.Hurrell,and C.Deser,2004:North Atlantic winter climate regimes:Spatial asymmetry,stationarity with time,and oceanic forcing.J.Climate,17,1055–1068.

    Croci-Maspoli,M.,C.Schwierz,and H.C.Davies,2007:Atmospheric blocking:Space–time links to the NAO and PNA.Climate Dyn.,29,713–725.

    Davini,P.,C.Cagnazzo,S.Gualdi,and A.Navarra,2012a:Bidimensional diagnostics,variability,and trends of Northern Hemisphere blocking.J.Climate,25,6496–6509.

    Davini,P.,C.Cagnazzo,R.Neale,and J.Tribbia,2012b:Coupling between Greenland blocking and the North Atlantic Oscillation pattern.Geophys.Res.Lett.,39,l14701,doi:10.1029/ 2012GL052315.

    Diao,Y.N.,J.P.Li,and D.H.Luo,2006:A new blocking index and its application:Blocking action in the Northern Hemisphere.J.Climate,19,4819–4839.

    Feldstein,S.B.,2003:The dynamics of NAO teleconnection pattern growth and decay.Quart.J.Roy.Meteor.Soc.,129,901–924.

    Franzke,C.,S.Lee,and S.B.Feldstein,2004:Isthe North Atlantic Oscillation a breaking wave?J.Atmos.Sci.,61,145–160.

    Glahn,H.R.,1968:Canonical correlation and its relationship to discriminant analysis and multiple regression.J.Atmos.Sci.,25,23–31.

    Hurrell,J.W.,1995:Decadal trends in the North Atlantic Oscillation:Regional temperatures and precipitation.Science,269, 676–679.

    Jung,T.,M.Hilmer,E.Ruprecht,S.Kleppek,S.K.Gulev,and O. Zolina,2003:Characteristics of the recent eastward shift of interannual NAO variability.J.Climate,16,3371–3382.

    Kalnay,E.,and Coauthors,1996:The NCEP/NCAR 40-Year Reanalysis Project.Bull.Amer.Meteor.Soc.,77,437–471.

    Lejen¨as,H.,and H.?kland,1983:Characteristics of northern hemisphere blocking as determined from a long time series of observational data.Tellus A,35A,350–362.

    Luo,D.H.,2005:Why is the North Atlantic block more frequent and long-lived during the negative NAO phase?Geophys.Res. Lett.,32,L20804,doi:10.1029/2005GL022927.

    Luo,D.H.,and J.Cha,2012:The North Atlantic Oscillation and the North Atlantic jet variability:Precursors to NAO regimesand transitions.J.Atmos.Sci.,69,3763–3787.

    Luo,D.H.,A.R.Lupo,and H.Wan,2007a:Dynamics of eddydriven low-frequency dipole modes.Part I:A simple model of North Atlantic Oscillations.J.Atmos.Sci.,64,3–38.

    Luo,D.H.,T.T.Gong,and Y.N.Diao,2007b:Dynamics of eddy-driven low-frequency dipole modes.Part III:Meridional displacement of westerly jet anomalies during two phases of NAO.J.Atmos.Sci.,64,3232–3248.

    Luo,D.H.,T.T.Gong,and Y.N.Diao,2008:Dynamics of eddydriven low-frequency dipole modes.Part IV:Planetary and synoptic wave-breaking processes during the NAO life cycle. J.Atmos.Sci.,65,737–765.

    Luo,D.H.,Z.H.Zhu,R.C.Ren,L.H.Zhong,and C.Z.Wang, 2010:Spatial pattern and zonal shift of the North Atlantic Oscillation.Part I:A dynamical interpretation.J.Atmos.Sci., 67,2805–2826.

    Luo,D.H.,J.Cha,and S.B.Feldstein,2012a:Weather regime transitions and the interannual variability of the North Atlantic Oscillation.Part I:A likely connection.J.Atmos.Sci., 69,2329–2346.

    Luo,D.H.,J.Cha,and S.B.Feldstein,2012b:Weather regime transitions and the interannual variability of the North Atlantic Oscillation.Part II:Dynamical processes.J.Atmos. Sci.,69,2347–2363.

    Luo,D.H.,Y.Yao,and A.G.Dai,2014:Decadal relationship between European blocking and North Atlantic oscillation during 1978–2011.Part I:Atlantic conditions.J.Atmos.Sci.,72, 1152–1173.

    Moore,G.W.K.,and I.A.Renfrew,2012:Cold European winters:Interplay between the NAO and the East Atlantic mode. Atmospheric Science Letters,13,1–8.

    Pelly,J.L.,and B.J.Hoskins,2003:A new perspective on blocking.J.Atmos.Sci.,60,743–755.

    Rivi`ere,G.,and I.Orlanski,2007:Characteristics of the Atlantic storm-track eddy activity and its relation with the North Atlantic Oscillation.J.Atmos.Sci.,64,241–266.

    Scaife,A.A.,C.K.Folland,L.V.Alexander,A.Moberg,and J. R.Knight,2008:European climate extremes and the North Atlantic oscillation.J.Climate,21,72–83,doi:10.1175/2007 JCLI1631.1.

    Scherrer,S.,M.Croci-Maspoli,C.Schwierz,and C.Appenzeller, 2006:Two-dimensional indices of atmospheric blocking and their statistical relationship with winter climate patterns in the Euro-Atlantic region.Inter.J.Climatol.,26,233–249.

    Schwierz,C.,M.Croci-Maspoli,and H.C.Davies,2004:Perspicacious indicators of atmospheric blocking.Geophys.Res. Lett.,31,L06125,doi:10.1029/2003GL019341.

    Shabbar,A.,J.P.Huang,and K.Higuchi,2001:The relationship between the wintertime North Atlantic Oscillation and blocking episodes in the North Atlantic.International Journal of Climatology,21,355–369,doi:10.1002/joc.612.

    Sillmann,J.,M.Croci-Maspoli,M.Kallache,and R.W.Katz, 2011:Extreme cold winter temperatures in Europe under the influence of North Atlantic atmospheric blocking.J.Climate, 24,5899–5913,doi:10.1175/2011JCLI4075.1.

    Sung,M.K.,G.H.Lim,J.S.Kug,and S.An,2011:A linkage between the North Atlantic Oscillation and its downstream development due to the existence of a blocking ridge.J.Geophys.Res.,116,doi:10.1029/2010JD015006.

    Strong,C.,and G.Magnusdottir,2008:Tropospheric Rossby wave breaking and the NAO/NAM.J.Atmos.Sci.,65,2861–2876.

    Tibaldi,S.,and F.Molteni,1990:On the operational predictability of blocking.Tellus A,42,343–365.

    Woollings,T.,B.Hoskins,M.Blackburn,and P.Berrisford,2008: A new Rossby wave–breaking interpretation of the North Atlantic Oscillation.J.Atmos.Sci.,65,609–626.

    Yao,Y.,and D.H.Luo,2014:Relationship between zonal position of the North Atlantic Oscillation and Euro-Atlantic blocking events and its possible effect on the weather over Europe.Science China Earth Sciences,57,2628–2636.

    :Yao,Y.,and D.H.Luo,2015:Do European blocking events precede North Atlantic Oscillation events?Adv. Atmos.Sci.,32(8),1106–1118,

    10.1007/s00376-015-4209-5.

    19 September 2014;revised 20 January 2015;accepted 2 February 2015)

    ?Corresponding author:LUO Dehai

    Email:ldh@mail.iap.ac.cn

    最近中文字幕高清免费大全6| 成年动漫av网址| 国产毛片在线视频| 久久久久久久久久久免费av| 黄频高清免费视频| 亚洲第一区二区三区不卡| 一区二区三区精品91| 在线观看免费午夜福利视频| 日韩大片免费观看网站| 咕卡用的链子| 91成人精品电影| 亚洲欧美成人精品一区二区| 91国产中文字幕| 精品久久蜜臀av无| 下体分泌物呈黄色| 考比视频在线观看| 老鸭窝网址在线观看| 色吧在线观看| 成年人午夜在线观看视频| 人人妻,人人澡人人爽秒播 | 精品一区二区三卡| 亚洲一码二码三码区别大吗| 午夜激情久久久久久久| 超碰成人久久| 日本欧美视频一区| 久久免费观看电影| 国产精品人妻久久久影院| 国产在视频线精品| 丝袜喷水一区| 天天操日日干夜夜撸| 少妇的丰满在线观看| 一区二区三区乱码不卡18| 国产深夜福利视频在线观看| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 一边摸一边抽搐一进一出视频| 欧美av亚洲av综合av国产av | 老司机影院成人| 国产亚洲av片在线观看秒播厂| 成人亚洲精品一区在线观看| 99热国产这里只有精品6| 久久久久精品国产欧美久久久 | av免费观看日本| 大香蕉久久成人网| 日本vs欧美在线观看视频| 一级毛片电影观看| 午夜久久久在线观看| 成人影院久久| 91老司机精品| 午夜91福利影院| 国产精品嫩草影院av在线观看| a级毛片黄视频| 亚洲欧美激情在线| 999精品在线视频| 精品国产一区二区三区四区第35| 国产成人一区二区在线| 成人免费观看视频高清| 成年美女黄网站色视频大全免费| 成人毛片60女人毛片免费| 精品国产一区二区久久| 男女边吃奶边做爰视频| 久久 成人 亚洲| 国产精品免费视频内射| 欧美亚洲日本最大视频资源| 美女午夜性视频免费| 日本黄色日本黄色录像| 亚洲国产欧美一区二区综合| 日韩制服骚丝袜av| 国产亚洲欧美精品永久| www.熟女人妻精品国产| 亚洲精品久久午夜乱码| 免费日韩欧美在线观看| 下体分泌物呈黄色| 男女午夜视频在线观看| 叶爱在线成人免费视频播放| 亚洲精品国产一区二区精华液| 国产老妇伦熟女老妇高清| 免费少妇av软件| 亚洲一码二码三码区别大吗| 9色porny在线观看| 亚洲精品在线美女| 看十八女毛片水多多多| 久久精品久久久久久噜噜老黄| 两个人看的免费小视频| 在线 av 中文字幕| 成人漫画全彩无遮挡| 菩萨蛮人人尽说江南好唐韦庄| 免费观看性生交大片5| 秋霞在线观看毛片| 国产又爽黄色视频| 美女脱内裤让男人舔精品视频| 黄片无遮挡物在线观看| 爱豆传媒免费全集在线观看| 丁香六月欧美| 建设人人有责人人尽责人人享有的| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 日韩不卡一区二区三区视频在线| av国产精品久久久久影院| 亚洲精品国产区一区二| 久久久久精品人妻al黑| 成人黄色视频免费在线看| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 97在线人人人人妻| 日本黄色日本黄色录像| 日本午夜av视频| 成人漫画全彩无遮挡| 国产精品99久久99久久久不卡 | 亚洲激情五月婷婷啪啪| e午夜精品久久久久久久| 精品卡一卡二卡四卡免费| 亚洲国产欧美日韩在线播放| 另类亚洲欧美激情| 99国产综合亚洲精品| 日本色播在线视频| 人人澡人人妻人| 日韩 亚洲 欧美在线| 国产色婷婷99| 午夜福利免费观看在线| 国产黄频视频在线观看| www.熟女人妻精品国产| 欧美日韩精品网址| 欧美日韩亚洲高清精品| 欧美老熟妇乱子伦牲交| 丝袜美腿诱惑在线| 亚洲熟女精品中文字幕| 国产片内射在线| 女人爽到高潮嗷嗷叫在线视频| 飞空精品影院首页| 亚洲,欧美精品.| 亚洲欧美一区二区三区黑人| 一级毛片 在线播放| 丁香六月欧美| 超碰97精品在线观看| 久久天堂一区二区三区四区| 国产精品av久久久久免费| 免费高清在线观看视频在线观看| av又黄又爽大尺度在线免费看| 亚洲av电影在线观看一区二区三区| 一级毛片我不卡| 老汉色∧v一级毛片| 午夜福利乱码中文字幕| 欧美日韩亚洲综合一区二区三区_| 啦啦啦在线观看免费高清www| 久久久久网色| 看免费成人av毛片| 中国三级夫妇交换| 一级a爱视频在线免费观看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱久久久久久| 亚洲精品第二区| 99国产精品免费福利视频| 国产一区二区在线观看av| 欧美 亚洲 国产 日韩一| 久热这里只有精品99| 街头女战士在线观看网站| 一二三四中文在线观看免费高清| 欧美精品高潮呻吟av久久| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 91成人精品电影| 叶爱在线成人免费视频播放| 91精品三级在线观看| 精品少妇久久久久久888优播| 国产免费视频播放在线视频| 少妇人妻精品综合一区二区| 伊人久久大香线蕉亚洲五| 91精品三级在线观看| 青草久久国产| 可以免费在线观看a视频的电影网站 | 99国产精品免费福利视频| av网站免费在线观看视频| 国产精品一区二区在线观看99| 欧美精品一区二区免费开放| 少妇的丰满在线观看| 男男h啪啪无遮挡| 在线观看www视频免费| 久久精品熟女亚洲av麻豆精品| 国产成人精品福利久久| 青草久久国产| 午夜免费鲁丝| 老司机靠b影院| 天天添夜夜摸| 黑人欧美特级aaaaaa片| 青春草视频在线免费观看| 国产精品国产三级专区第一集| 国产高清国产精品国产三级| 久久韩国三级中文字幕| 一本一本久久a久久精品综合妖精| 午夜福利免费观看在线| 91aial.com中文字幕在线观看| 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| 亚洲av在线观看美女高潮| 亚洲精华国产精华液的使用体验| 久久久久久久久久久久大奶| 国产精品 欧美亚洲| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜| 国产精品三级大全| 一个人免费看片子| 中文字幕人妻熟女乱码| www日本在线高清视频| kizo精华| av网站在线播放免费| 亚洲精品中文字幕在线视频| 男女下面插进去视频免费观看| 最近中文字幕高清免费大全6| 国产成人系列免费观看| 久久精品久久久久久噜噜老黄| 丝袜脚勾引网站| 女人久久www免费人成看片| 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 亚洲av综合色区一区| 女性被躁到高潮视频| 午夜福利视频在线观看免费| 免费观看人在逋| 国产男人的电影天堂91| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 久久精品人人爽人人爽视色| av福利片在线| 国产一区二区 视频在线| 国产黄色视频一区二区在线观看| 日韩电影二区| 在线看a的网站| 欧美中文综合在线视频| 黄色 视频免费看| 久久99一区二区三区| 久久青草综合色| 成年动漫av网址| 亚洲欧美成人综合另类久久久| 亚洲欧美一区二区三区国产| 黑丝袜美女国产一区| e午夜精品久久久久久久| 秋霞在线观看毛片| 国产探花极品一区二区| 男女高潮啪啪啪动态图| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人澡人人看| 9热在线视频观看99| 免费看不卡的av| 日本午夜av视频| 亚洲欧美一区二区三区久久| 精品卡一卡二卡四卡免费| 黄色毛片三级朝国网站| 亚洲第一av免费看| 午夜福利在线免费观看网站| 亚洲天堂av无毛| 国产av国产精品国产| 亚洲成人免费av在线播放| 国产精品av久久久久免费| 久久精品久久久久久噜噜老黄| 卡戴珊不雅视频在线播放| av在线app专区| 亚洲精品自拍成人| 我的亚洲天堂| 亚洲国产精品国产精品| 欧美日韩综合久久久久久| 一本色道久久久久久精品综合| 亚洲自偷自拍图片 自拍| 男女无遮挡免费网站观看| 男女下面插进去视频免费观看| 国产精品免费大片| 久热爱精品视频在线9| 欧美人与善性xxx| 一级片'在线观看视频| 18禁观看日本| 亚洲熟女精品中文字幕| 国产精品一区二区精品视频观看| 国产熟女欧美一区二区| 9191精品国产免费久久| 女性生殖器流出的白浆| 日日爽夜夜爽网站| 亚洲国产日韩一区二区| 青青草视频在线视频观看| 午夜免费观看性视频| 桃花免费在线播放| 三上悠亚av全集在线观看| 超色免费av| av女优亚洲男人天堂| 久久女婷五月综合色啪小说| 久久 成人 亚洲| 一本一本久久a久久精品综合妖精| 亚洲精品国产一区二区精华液| 高清在线视频一区二区三区| 国产精品秋霞免费鲁丝片| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区大全| 成年人午夜在线观看视频| 黑丝袜美女国产一区| 国产av国产精品国产| 久久亚洲国产成人精品v| 日韩av免费高清视频| 亚洲图色成人| 视频区图区小说| 一区二区三区激情视频| 欧美亚洲 丝袜 人妻 在线| 亚洲一区中文字幕在线| 国产女主播在线喷水免费视频网站| 99九九在线精品视频| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 另类亚洲欧美激情| 久久久久精品国产欧美久久久 | 亚洲欧洲日产国产| 久久婷婷青草| 欧美日韩亚洲国产一区二区在线观看 | 夜夜骑夜夜射夜夜干| 久久久国产精品麻豆| xxxhd国产人妻xxx| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 久久天堂一区二区三区四区| 美女扒开内裤让男人捅视频| 欧美97在线视频| 中国三级夫妇交换| 在线观看人妻少妇| 亚洲精品乱久久久久久| 超碰97精品在线观看| 99精国产麻豆久久婷婷| 在线观看三级黄色| 一边摸一边做爽爽视频免费| 19禁男女啪啪无遮挡网站| 最新在线观看一区二区三区 | 亚洲综合精品二区| 日韩精品免费视频一区二区三区| 国产精品人妻久久久影院| 亚洲少妇的诱惑av| 女人精品久久久久毛片| 男女午夜视频在线观看| 一区二区三区精品91| 国产人伦9x9x在线观看| 中文字幕人妻熟女乱码| 啦啦啦视频在线资源免费观看| 国产乱人偷精品视频| 男的添女的下面高潮视频| 看免费av毛片| 在线观看人妻少妇| 视频在线观看一区二区三区| 啦啦啦在线免费观看视频4| 免费少妇av软件| 免费不卡黄色视频| 在线观看人妻少妇| 久久久久国产一级毛片高清牌| 日本91视频免费播放| 在线观看免费午夜福利视频| 九九爱精品视频在线观看| 欧美人与性动交α欧美精品济南到| 热re99久久精品国产66热6| 18在线观看网站| 少妇人妻精品综合一区二区| 亚洲精品日韩在线中文字幕| 亚洲五月色婷婷综合| 欧美精品人与动牲交sv欧美| 777久久人妻少妇嫩草av网站| 欧美在线一区亚洲| 少妇人妻精品综合一区二区| 欧美亚洲 丝袜 人妻 在线| 秋霞在线观看毛片| 日本av免费视频播放| 青春草国产在线视频| 亚洲av中文av极速乱| 亚洲精品第二区| 亚洲av中文av极速乱| 亚洲五月色婷婷综合| 18在线观看网站| 国产乱来视频区| 免费人妻精品一区二区三区视频| 欧美亚洲日本最大视频资源| 精品国产乱码久久久久久小说| 中文字幕亚洲精品专区| 国产一区二区在线观看av| 人人妻,人人澡人人爽秒播 | 亚洲伊人久久精品综合| videos熟女内射| 日韩大片免费观看网站| 青春草视频在线免费观看| 熟女少妇亚洲综合色aaa.| h视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久久久精品古装| 街头女战士在线观看网站| 国产精品久久久久久精品古装| 男人操女人黄网站| 亚洲精品国产av蜜桃| 欧美日韩一级在线毛片| 亚洲图色成人| 中文字幕色久视频| 欧美黑人精品巨大| 久久热在线av| 日本黄色日本黄色录像| 精品卡一卡二卡四卡免费| 最近中文字幕2019免费版| 无限看片的www在线观看| 天天添夜夜摸| 看免费av毛片| 熟妇人妻不卡中文字幕| 天天躁夜夜躁狠狠躁躁| 欧美激情高清一区二区三区 | 国产精品久久久av美女十八| 久久天躁狠狠躁夜夜2o2o | 成年人午夜在线观看视频| 国产伦人伦偷精品视频| 女性被躁到高潮视频| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| tube8黄色片| 街头女战士在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 国产又爽黄色视频| 成人国产麻豆网| 一边摸一边做爽爽视频免费| 欧美xxⅹ黑人| 亚洲婷婷狠狠爱综合网| 国产av精品麻豆| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 亚洲欧美精品综合一区二区三区| 成人手机av| 中文字幕制服av| 精品一区二区免费观看| 精品国产一区二区三区四区第35| 免费av中文字幕在线| 亚洲国产精品999| 看免费av毛片| 美女脱内裤让男人舔精品视频| 色婷婷久久久亚洲欧美| 日韩av在线免费看完整版不卡| 亚洲熟女毛片儿| 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| 国产成人91sexporn| 999久久久国产精品视频| 黄片无遮挡物在线观看| 国产精品一区二区精品视频观看| 日韩熟女老妇一区二区性免费视频| 黑人欧美特级aaaaaa片| 久久久精品免费免费高清| 久久久精品国产亚洲av高清涩受| 在线免费观看不下载黄p国产| 成人亚洲欧美一区二区av| xxxhd国产人妻xxx| 中文字幕人妻丝袜一区二区 | 免费女性裸体啪啪无遮挡网站| av女优亚洲男人天堂| 欧美亚洲 丝袜 人妻 在线| 如何舔出高潮| 一本—道久久a久久精品蜜桃钙片| 国产精品香港三级国产av潘金莲 | 国产精品久久久久成人av| 久久av网站| 午夜福利视频精品| 高清黄色对白视频在线免费看| 老司机在亚洲福利影院| 成人黄色视频免费在线看| 亚洲av在线观看美女高潮| 在线观看三级黄色| 久久av网站| 中文字幕高清在线视频| 国产成人精品无人区| 国产乱来视频区| 久久久久精品久久久久真实原创| 国产免费又黄又爽又色| 精品第一国产精品| 免费高清在线观看日韩| 午夜激情av网站| 欧美日韩视频精品一区| 久久久国产欧美日韩av| 欧美少妇被猛烈插入视频| 亚洲精品美女久久久久99蜜臀 | 亚洲伊人色综图| 自拍欧美九色日韩亚洲蝌蚪91| 国产一卡二卡三卡精品 | 亚洲国产精品999| 久久韩国三级中文字幕| 欧美在线一区亚洲| 亚洲国产欧美日韩在线播放| 久久久久视频综合| 熟女av电影| 啦啦啦视频在线资源免费观看| 国产伦人伦偷精品视频| 久久久精品区二区三区| 国产成人免费观看mmmm| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 人成视频在线观看免费观看| 国产成人欧美| 成年女人毛片免费观看观看9 | 在线观看国产h片| 在线天堂中文资源库| 中文字幕精品免费在线观看视频| 这个男人来自地球电影免费观看 | 国产成人91sexporn| 精品亚洲乱码少妇综合久久| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 黄片无遮挡物在线观看| 国产极品粉嫩免费观看在线| 一区二区三区激情视频| 亚洲精品视频女| 美国免费a级毛片| 成年人免费黄色播放视频| 久久午夜综合久久蜜桃| 赤兔流量卡办理| 午夜激情久久久久久久| 国产色婷婷99| xxx大片免费视频| 国产日韩欧美在线精品| 成人三级做爰电影| 国产免费又黄又爽又色| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频在线观看免费| 亚洲成国产人片在线观看| 国产成人系列免费观看| 熟女少妇亚洲综合色aaa.| av视频免费观看在线观看| 蜜桃国产av成人99| 18禁裸乳无遮挡动漫免费视频| 欧美av亚洲av综合av国产av | 97在线人人人人妻| 18禁裸乳无遮挡动漫免费视频| 在线观看一区二区三区激情| 亚洲国产欧美在线一区| 搡老岳熟女国产| 一级片'在线观看视频| 久久久精品免费免费高清| √禁漫天堂资源中文www| 十八禁高潮呻吟视频| 99热网站在线观看| 日本91视频免费播放| 国产精品人妻久久久影院| 精品一品国产午夜福利视频| 久久久国产一区二区| 日韩一区二区视频免费看| 伦理电影免费视频| 欧美激情 高清一区二区三区| 日韩成人av中文字幕在线观看| 成年人免费黄色播放视频| 中文字幕高清在线视频| 最近中文字幕2019免费版| 中文乱码字字幕精品一区二区三区| 天天躁夜夜躁狠狠久久av| 97精品久久久久久久久久精品| 十八禁人妻一区二区| 悠悠久久av| 亚洲国产av影院在线观看| 国产日韩欧美亚洲二区| 日韩大片免费观看网站| 亚洲欧美色中文字幕在线| 亚洲国产中文字幕在线视频| 久久久欧美国产精品| 国产一区二区 视频在线| 日本黄色日本黄色录像| 国产1区2区3区精品| 夫妻性生交免费视频一级片| 婷婷色综合www| 国产成人精品久久二区二区91 | 久久精品aⅴ一区二区三区四区| 亚洲国产中文字幕在线视频| 久久毛片免费看一区二区三区| 天美传媒精品一区二区| 亚洲中文av在线| 欧美日韩国产mv在线观看视频| 制服丝袜香蕉在线| 黄色视频在线播放观看不卡| 久久性视频一级片| 纯流量卡能插随身wifi吗| 精品久久蜜臀av无| 精品酒店卫生间| 久热爱精品视频在线9| av片东京热男人的天堂| 亚洲久久久国产精品| a级毛片黄视频| 一区二区av电影网| 久久久欧美国产精品| 丁香六月欧美| 国产伦人伦偷精品视频| 秋霞在线观看毛片| 日本一区二区免费在线视频| 国产又爽黄色视频| 久久午夜综合久久蜜桃| 黄色视频不卡| 99久久综合免费| 欧美xxⅹ黑人| 国产精品一区二区精品视频观看| 两性夫妻黄色片| 老熟女久久久| 久久久久国产一级毛片高清牌| 国产精品国产av在线观看| 免费观看性生交大片5| 国产男女超爽视频在线观看| 免费日韩欧美在线观看| 亚洲精品久久久久久婷婷小说| 国产欧美日韩一区二区三区在线| 日韩大片免费观看网站| 亚洲在久久综合| 欧美国产精品va在线观看不卡| 日本av免费视频播放| 国产日韩欧美在线精品| 91精品伊人久久大香线蕉| 亚洲精品,欧美精品| 久久久国产欧美日韩av| 校园人妻丝袜中文字幕| 纵有疾风起免费观看全集完整版| 精品亚洲乱码少妇综合久久| 我的亚洲天堂| 一级,二级,三级黄色视频| 波多野结衣av一区二区av| 啦啦啦视频在线资源免费观看| 亚洲精品久久久久久婷婷小说| 久久国产精品男人的天堂亚洲| 国产成人a∨麻豆精品|