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

    Predictability of the Summer East Asian Upper-Tropospheric Westerly Jet in ENSEMBLES Multi-Model Forecasts

    2015-06-09 21:24:03LIChaofanandLINZhongda
    Advances in Atmospheric Sciences 2015年12期

    LI Chaofanand LIN Zhongda

    1Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2State Key Laboratory of Numerical Modelling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    Predictability of the Summer East Asian Upper-Tropospheric Westerly Jet in ENSEMBLES Multi-Model Forecasts

    LI Chaofan1and LIN Zhongda?2

    1Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2State Key Laboratory of Numerical Modelling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    The interannual variation of the East Asian upper-tropospheric westerly jet(EAJ)significantly affects East Asian climate in summer.Identifying its performance in model prediction may provide us another viewpoint,from the perspective of uppertropospheric circulation,to understand the predictability of summer climate anomalies in East Asia.This study presents a comprehensive assessment of year-to-year variability of the EAJ based on retrospective seasonal forecasts,initiated from 1 May,in the f i ve state-of-the-art coupled models from ENSEMBLES during 1960–2005.It is found that the coupled models show certain capability in describing the interannual meridional displacement of the EAJ,which reflects the models’performance in the first leading empirical orthogonal function(EOF)mode.This capability is mainly shown over the region south of the EAJ axis.Additionally,the models generally capture well the main features of atmospheric circulation and SST anomalies related to the interannual meridional displacement of the EAJ.Further analysis suggests that the predicted warm SST anomalies in the concurrent summer over the tropical eastern Pacific and northern Indian Ocean are the two main sources of the potential prediction skill of the southward shift of the EAJ.In contrast,the models are powerless in describing the variation over the region north of the EAJ axis,associated with the meridional displacement,and interannual intensity change of the EAJ,the second leading EOF mode,meaning it still remains a challenge to better predict the EAJ and, subsequently,summer climate in East Asia,using current coupled models.

    East Asian westerly jet,seasonal prediction,coupled model,meridional displacement

    1.Introduction

    The East Asian upper-tropospheric westerly jet(EAJ), one of the most important components in the Asian monsoon system,is intimately related to the East Asian climate anomalies in boreal summer(Lin,2013;Lu et al.,2013).The interannual variations in both the EAJ location and intensity correspond closely to the rainfall and surface air temperature anomalies over East Asia,especially for the meridional displacement of the EAJ location(Lau et al.,2000;Lu,2004; Wang et al.,2013).The poleward displacement of the EAJ tends to be linked to def i cient rainfall along the East Asian subtropical rainy belt(Lau et al.,2000)and more frequent high temperature extremes in Southeast China(Wang et al., 2013).Thus,interannual variation of the EAJ is expected to have some implications for the predictability of summer climate anomalies over East Asia(Liang and Wang,1998; Zhang and Guo,2005).As the current status of seasonal prediction for East Asian summer rainfall remains a challenge,it is imperative to comprehensively assess the prediction of the summer EAJ.

    The interannual variation of the EAJ is characterized by meridional displacement of the EAJ location and intensity change,which are the first and second leading EOF modes of the EAJ,respectively(Lin and Lu,2005).The atmospheric circulation anomalies related to the meridional shift of the EAJ location involve two wave trains:the Pacific–Japan(PJ) teleconnection(Nitta,1987;Kosaka and Nakamura,2006;Lu and Lin,2009)or the East Asia–Pacific(EAP)teleconnection(Huang and Sun,1992;Huang,2004)in the meridional direction over the East Asia–western North Pacific(WNP) region,and the“Silk Road”teleconnection in the zonal direction along the upper-tropospheric Asian westerly jet(Lu et al.,2002;Enomoto et al.,2003;Enomoto,2004;Ding and Wang,2005).In addition,the transient eddy forcing also demonstrates signifi cant impact on the meridional displacement of the EAJ,especially its northward progression (Xiang and Yang,2012).On the other hand,Lu(2004)also revealed a signifi cant meridional connection of the EAJ in-tensity change with convective activities over the Philippine Sea,though it is weaker than that related to the meridional displacement of the EAJ location.

    Despite all the above atmospheric forcing,the interannual meridional displacement of the EAJ location is also infl uenced by the tropical SST anomalies.Lin(2010)found that,associated with the southward shift of the EAJ in July and August,a warm SST anomaly in the tropical eastern Pacific persistsfrom the preceding spring to the concurrentsummer.Additionally,Qu and Huang(2012)revealed that the SST anomalies over the tropical Indian Ocean are also signi ficantly related to the meridional displacement of the EAJ. When SST warms in the tropical Indian Ocean,especially the northern Indian Ocean(NIO),the EAJ shifts southward.

    However,atmospheric general circulation models (AGCMs)are de fi cient in describing the SST-induced interannual variation of the EAJ.Lu et al.(2006)investigated the external and internal atmospheric variability over the WNP and East Asia based on the simulations forced by the prescribed observed SST in the Met Of fi ce AGCM,known as HadAM3.They indicated that the interannual anomalies of the EAJ,including both meridional displacement of the EAJ location and its intensity change,are dominated by atmospheric internal variability,different from the lowertropospheric circulation anomaly over the WNP,in which external variability plays an important role.The strong effect of the internal variability suggests a low potential predictability of the EAJ variation in the AGCM.In addition,current AGCMs are unable to properly simulate the Asian summer monsoon rainfall because the ocean–atmosphere coupling contributes signi ficantly to the climate anomalies over the WNP and East Asia(Wang et al.,2004,2005).This will further increase the possibility of AGCM de fi ciencies in the prediction of the summer EAJ,in terms of the meridional teleconnection between the interannual variation of the EAJ and convective activities over the WNP and East Asia.As a consequence,ocean–atmosphere coupled general circulation models(CGCMs)are more appropriate for the prediction of the interannual variation of the summer EAJ.

    Some studies have analyzed the potential predictability of summer climate and circulation anomalies in the Northern Hemisphere in CGCMs(Chowdary et al.,2010;Lee et al., 2011;Kosaka et al.,2012;Li et al.,2012,2014).Their results imply that atmosphere–ocean interactions,induced by the tropical Indian and Pacific oceans,contribute remarkably to the prediction reliability.Lee et al.(2011),for example, examined the predictability of summer upper-tropospheric circulation in the Northern Hemisphere.They pointed out that the seasonal forecast skill for the 200 hPa geopotential height derived mainly from the CGCM’s ability to predict the first two EOF modes in the Northern Hemisphere,for the 25-year period of 1981–2005.The first two EOF modes are characteristics of zonally uniform anomalies of uppertropospheric geopotential height over the tropics(EOF1), due to the prolonged impact of El Nin~o–Southern Oscillation(ENSO),and over the midlatitudes(EOF2)related to the developing ENSO.The zonally uniform anomalies of geopotential height at 200 hPa suggest a potential predictability of upper-tropospheric zonal winds in the context of geostrophic wind balance.Nevertheless,the current status of the seasonal predictability of the EAJ variation by CGCMs remains unclear and has not yet been well documented in the literature.

    Recently,the ENSEMBLES project(ensembles-based predictions of climate changes and their impacts),advocated by the European Union(van der Linden and Mitchell,2009), provided seasonal multi-model prediction products of f i ve leading CGCMs developed in Europe for the period 1960–2005.Based on the retrospective forecast(hindcast)of ENSEMBLES,Li et al.(2012)investigated the prediction of the summer WNP subtropical high and suggested that ENSEMBLES shows considerable capability in describing the interannual variation of the summer lower-tropospheric circulation and precipitation anomalies over the WNP.Moreover, related to the WNP subtropical high,the meridional teleconnection,in the lower troposphere,can be captured by these coupled models.

    This study attempts to comprehensively assess the predictability of the interannual variation of the EAJ,based on the products of the ENSEMBLES forecast system.Assessment of the EAJ will hopefully provide us another viewpoint, from the perspective of upper-tropospheric circulation,on the seasonal predictability of the summer monsoon and rainfall over East Asia.The rest of this paper is organized as follows.Section 2 describes the model hindcasts and observed datasets used in this study.Section 3 presents the comprehensive assessment of the EAJ,including its climatology,interannual variability,temporal correlation coeff icients,two leading EOF modes,circulation and tropical SST anomalies associated with the meridional displacement of the EAJ location.The predicted decadal change and potential prediction skill of the EAJ location and its implications are discussed in section 4.Conclusions are summarized in section 5.

    2.Model hindcasts and observational datasets

    The ENSEMBLES project was a seasonal-to-annual multi-model project developed by the European Union(van der Linden and Mitchell,2009).It comprised fi ve fully coupled atmosphere–ocean models,including the UK Met Offi ce(UKMO),the Me′te′o-France(MF),the European Centre for Medium-Range Weather Forecasts(ECMWF),the Leibniz Institute of Marine Sciences at Kiel University(IFMGEOMAR)and the Euro-Mediterranean Center for Climate Change(CMCC-INGV).The atmosphere and ocean were initialized using realistic estimates of their observed states.All of these models included major radiative forcing and had no flux adjustments.An ensemble of nine initial conditions was run for each model.Further details on the ENSEMBLES multi-model project,the models and the initial condition perturbations can be found in Weisheimer et al.(2009) and Doblas-Reyes et al.(2010).

    A 46-year hindcast from 1960 to 2005 was carried out in the fi ve models.They were initialized on 1 May for eachyear and run for seven months.The hindcast results for June, July and August are analyzed in this study.Additionally,the multi-model-ensemble(MME)results are calculated through a simple composite by applying equal weight to all the f i ve models.

    The observational datasets used here for validating the model simulation include:(1)National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR)monthly reanalysis data from 1960 to 2005 (Kalnay et al.,1996);(2)National Oceanic and Atmospheric Administration Extended Reconstructed monthly mean SST V3 data from 1960 to 2005(Smith and Reynolds,2004);and (3)monthly precipitation data obtained from the Global Precipitation Climatology Project during 1979–2005(Adler et al.,2003).

    3.Results

    In this section,we show the performance of the prediction of the EAJ in the f i ve models in ENSEMBLES.The predicted climatology and interannual standard deviation are presented in subsection 3.1,the prediction skill of the interannual variation of the EAJ in subsection 3.2,and the circulation and SST anomalies associated with the meridional shift of the EAJ location in subsections 3.3 and 3.4.

    3.1.Climatology and interannual variability

    Figure 1 shows the climatology of JJA(June–July–August)-mean 200 hPa zonal wind.In the observation,the upper-tropospheric jet stream is basically zonally oriented over East Asia.The axis of the EAJ,with the maximum 200 hPa zonal wind in each longitude over East Asia,is located at approximately 40?N(Fig.1a).The above features of the EAJ are successfully captured by the f i ve models and the MME prediction(Figs.1b–g).The location and intensity of the summer-mean EAJ described by the models are closely consistent with those shown in the observation,with a pattern correlation coeff i cient ranging from 0.96 to 0.98.The MME prediction performs generally better than the individual models,as it can reduce the noise presented in the individual forecasts.Relatively,the intensity of the JJA-mean zonal wind over East Asia is a little stronger in the ECMWF and IFM-GEOMAR models and weaker in the CMCC-INGV model.The axis of the EAJ is located a little south in the MF model than that in the observation.

    The interannual variation of upper-tropospheric zonal wind,during 1960–2005 in the observation,mostly appears in the midlatitudes between 20?N and 60?N,with two maxima residing to the south and north sides of the EAJ axis, respectively(Fig.2a).The maximum interannual standard deviation is greater than 3.5 m s?1both to the north and south sides of the EAJ axis.In terms of the model predictions,though it is weaker than observed,the maximum standard deviation,greater than 2 m s?1in the f i ve models and 1.5 m s?1in their MME prediction,is captured over the region south of the EAJ axis(Figs.2b–g).The result suggests these current coupled models show certain capability in capturing the interannual variation of the EAJ,particularly when it shifts southward.On the other hand,the standard deviation to the north of the EAJ axis is much weaker than that in the observation.This weaker variability is related to the ensemble mean of different initial conditions for the model prediction,and it further suggests that many synoptic perturbations contribute to the variations of the EAJ,especially to the north of the EAJ axis(e.g.,Dole and Black,1990;Wu et al.,2006;Xiang and Yang,2012).This is also the reason that the interannual variability of the MME prediction is a little weaker than that of the fi ve individual models.Among the fi ve models,the MF model(Fig.2e)predicts a weaker standard deviation over East Asia than the other four models.

    3.2.Prediction skill of the EAJ interannual variation

    Figure 3 illustrates the spatial distribution of the temporal correlation coeff i cient(TCC)between the model predictions and observation for 200 hPa zonal wind during 1960–2005. The MME prediction performs generally the best(Fig.3a). Good predictions of the upper-tropospheric zonal wind are found mainly over the south of the EAJ axis,extending northeastward from India to Japan,where the TCC skill over most of these regions is signifi cant at the 95%confidence level in the MME and all f i ve models except for the MF model. Consistent with the performance of the interannual variability (Fig.2e),the MF model exhibits relatively weaker capability in describing the variation of 200 hPa zonal wind over East Asia and the good prediction skill is unable to extend northeastward to North China and Japan(Fig.3d).Compared with that south of the EAJ axis,the skill over the region north of the EAJ axis is lower,though certain skill is also found over the north of the EAJ axis,especially over the east of Lake Baikal,in the MME,ECMWF,IFM-GEOMAR,and UKMO predictions.It is interesting to note that good predictions are located on the two sides of the EAJ axis,not in the axis region itself.This may suggest that the coupled models tend to better predict the meridional displacement of the EAJ location rather than its intensity change.

    Next,we investigate how well the models hindcasts describe the dominant modes of the year-to-year variability of the summer 200 hPa zonal wind over East Asia.To identify the major modes of the EAJ,we apply EOF analysis to the JJA 200 hPa zonal wind over the region(27.5?–55?N,120?–150?E),with strong interannual variability(Fig.2),following Lin and Lu(2005).

    Figure 4 shows the spatial distribution of the first leading EOF mode(EOF1)of the JJA 200 hPa zonal wind for the observation,the MME and the predictions of the f i ve models. EOF1,which accounts for 50%of the totalvariance in the observation,is characterized by the meridional variation of the EAJ location,with positive(negative)anomalies to the south (north)of the EAJ axis(Fig.4a).The models show good capability in describing the spatial distribution of EOF1(Figs. 4b–g).Variation of the EAJ meridional displacement for the model prediction averaged by different ensemble members accounts for a much larger variance,especially for the MMEprediction that comprises more members.This suggests that these coupled models are apt at reproducing the distribution pattern for this variation.Nevertheless,the models slightly overestimate the variation of the positive anomalies over the southern extent and underestimate the negative anomalies over the northern extent.This corresponds well to the interannual variability exhibited by the models,as stronger variability is found to the south than the north side of the EAJ axis(Figs.2b–g).

    The models are also capable of predicting the temporal evolution corresponding to EOF1(Table 1).The correlation coeff i cient of the associated principal component time series (PC1)between the MME prediction and the observation is 0.4,which is signifi cant at the 99.5%confidence level.Similarly,most of the f i ve models,except for the MF model,show similarly good capabilities.The IFM-GEOMAR model performs the best among the f i ve models and the prediction correlation coeff i cient reaches 0.5.The prediction correlation coeff i cient of 0.28 for the ECMWF and UKMO models exceeds the 94%confidence level,but is not greater than the 95%confidence level.The good performance of the coupled models in terms of EOF1 offers certain realizable predictability of the interannual meridional variation of the EAJ location.

    Table 1.Correlation coeff icients of the first two principal components(PC1 and PC2)between the observation and model predictions.The values in bold and those underlined represent statistical significance exceeding the 99%and 95%confidence levels,respectively.

    The second dominant mode(EOF2)features the changes of the EAJ intensity(Fig.5).It is associated with signifi cant positive values along the EAJ axis,accounting for 21.8%of the total variance in the observation.In general,a few models,including the ECMWF and IFM-GEOAR models,and the MME,tend to properly reproduce the spatial distribution of the EOF2 mode,but with a smaller explainable variance (Figs.5b–d).The other three models show similar explainable variance to the observation,but the positive anomalies extend more northward than the observation(Figs.5e–g), which contribute significantly to the northward distribution of the positive anomalies in the MME prediction(Fig.5b).In addition,the models are powerless in describing the interannual variation of the associated PC time series for the EOF2 mode(PC2).The correlation coeff icients are around zero forthe predictions of the MME and f i ve models(Table 1).Thus, the current coupled models have diff iculty in capturing the interannual intensity change of the EAJ.

    3.3.Circulation anomalies related to the meridional shift of the EAJ location

    Lu(2004)def i ned an index to depict the meridional shift of the EAJ location by the difference of 200 hPa zonal winds averaged over 120?–150?E between 10?to the south of the EAJ axis and 10?to the north.Similarly,we def i ne a JJA-mean EAJ location index(EAJLI)as the difference of 200 hPa zonal winds averaged between the regions(30?–40?N, 120?–150?E)and(40?–50?N,120?–150?E).Indeed,this EAJLI depicts the first EOF leading mode of interannual variations of the EAJ(Lin and Lu,2005).The EAJ moves southward when the EAJLI is positive,and it shifts northward when the EAJLI is negative.Time series of the EAJLI in the observation,MME,and f i ve models are shown in Fig.6.The correlation coeff i cient between the EAJLI and PC1 is 0.92 in the observation,0.96 in the MME,and larger than 0.9 in all of the f i ve models for the period 1960–2005.In the following analysis,the EAJLI,instead of PC1,is used to depict meridional shift of the EAJ location because it is more convenient to discuss the relative contribution from 200 hPa zonal winds on the southern(30?–40?N,120?–150?E)and northern(40?–50?N,120?–150?E)sides,respectively.

    The prediction skill of the EAJLI is presented in Table 2.The correlation coeff i cient of the EAJLI between the MME and observation is 0.35 during 1960–2005,which is signifi cant at the 95%confidence level.Similarly,signifi cant correlation coeff icients are also found in the IFM-GEOMAR(0.31)and UKMO(0.34)models,at the 95%con fi dence level,and in the CMCC-INGV(0.27)model,at the 90%confi dence level.The skill of the EAJLI is mainly attributed to that of the southern side.The correlation coef fi cient of the southern side is 0.51 in the MME,0.59 in the IFM-GEOMAR model,0.38 in the UKMO model,and 0.44 in the CMCCINGV model,all signi ficant at the 99%con fi dence level.For the other two models,the correlation coeff i cient of the southern side is 0.33 in the ECMWF model,signifi cant at the 95% confi dence level,and 0.22 in the MF model,though the correlation coeff icients of the corresponding EAJLI for these two models are only 0.19 and 0.07,respectively.However,the prediction skill of the northern side is poor,with a maximum correlation coeff i cient of 0.17 and minimum of?0.1between the observation and the f i ve individual models and their MME.

    Table 2.Correlation coeff icients of the EAJ location index(EAJLI) and its two sides between the observation and model predictions. The values in bold and those underlined represent statistical significance exceeding the 99%and 95%confidence levels,respectively.

    Further analysis shows that the main features of the circulation anomalies in the upper and lower troposphere,related to the EAJLI,are well predicted by the models,as illustrated in Fig.7.In the observation,associated with the EAJLI are a westerly anomaly to the south of 40?N and an easterly anomaly to the north over East Asia in the upper troposphere (Fig.7a),consistent with a southward shift tendency of the EAJ.In addition,the westerly and easterly wind anomalies over East Asia extend westward into West Asia so that the Asian westerly jet together moves southward.In the lower troposphere,a meridional wave train is observed,with a cyclonic anomaly over East Asia and an anticyclonic anomaly over the western North Pacific(WNP),i.e.,the PJ or EAP pattern.This meridional connection between the EAJ and the lower-tropospheric circulation anomalies over the East Asia–WNP region is consistent with the result previously identif i edby Lu and Lin(2009).Associated with the meridional connection,rainfall is enhanced over East Asia and suppressed over the WNP,Northeast China,and southern Russia(Fig. 7e).

    The circulation anomalies in the upper and lower troposphere and rainfall anomalies related to the EAJLI,discussed above in the observation,are successfully captured by the model predictions.In the MME prediction,associated withthe EAJLI,the Asian westerly jet moves southward in the upper troposphere(Fig.7b)and the meridional connection is established with a cyclonic anomaly over East Asia and an anticyclonic anomaly over the WNP in the lower troposphere(Fig.7d).Meanwhile,rainfall increases over East Asia and decreases over the WNP(Fig.7f).The successful reproducibility of the enhanced rainfall over East Asia further presents an implication for the important contribution of the EAJ displacement to the East Asian summer prediction.However,the relationship between the EAJ and tropical anomalies from the East Indian Ocean to the western Pacific, in the MME prediction,is significantly strengthened,compared with that in the observation.Related to the EAJLI,in the MME prediction,there are signifi cant westerly anomalies over the tropical regions in the upper troposphere(Fig.7b), and westerly anomalies over the tropical western Pacific and easterly anomalies over the East Indian Ocean in the lower troposphere(Fig.7d),concurrent with increased rainfall in the tropical western Pacific and decreased rainfall over the Maritime Continent(Fig.7f).This enhanced EAJLI–tropics connection is due to the overestimated effect of the tropical SST in the central-eastern Pacific and Indian Ocean in the model predictions,which is discussed in the next subsection.

    3.4.Tropical SST-EAJ location relationship

    The next step is to investigate what is responsible for the prediction skill of the EAJLI.To explore the possible impact of external forcing,we regress JJA-mean SST anomalies onto the EAJLI in the observation(Fig.8a)and the MME prediction(Fig.8b).As shown in Fig.8a,in the observation,the southward shift of the EAJ is significantly correlated with simultaneous warm SST anomalies in the tropical eastern Pacific and Indian Ocean,especial the northern Indian Ocean (NIO).

    The tropical SST anomalies depicted above are successfully reproduced in the MME prediction(Fig.8b),including the warm SST anomalies in the tropical eastern Pacific and NIO.This result therefore suggests a possible role of the tropical SST anomalies in the prediction of the EAJ location. However,the signals are much stronger in the model predictions than in the observation,arising from the average of different ensemble members in the model predictions.Here,we use two SST indices,the Ni~no3 index and NIO index(NIOI), to represent the interannual variation of SST in the tropical eastern Pacific and NIO,respectively.The Ni~no3 index is def i ned as the JJA-mean SST anomalies averaged over the region(5?S–5?N,90?–150?W)and the NIOI over the region (0?–20?N,40?–100?E).The correlation coeff i cient between the EAJLI and Ni~no3 index increases from 0.26 in the observation to 0.73 in the MME,and it ranges from 0.35 to 0.66 in the f i ve individual models(Table 3).Similarly,the correlation coeff i cient between the EAJLI and NIOI rises from 0.44 in the observation to 0.67 in the MME,and it changes from 0.52 to 0.63 in all of the models except the UKMO(0.44)and CMCC-INGV(0.39)models(Table 3).

    Table 3.Correlation coeff icients between the EAJLI and SST anomalies averaged over the Ni~no3 and NIO region.

    To confi rm the crucial role of the tropical SST anomalies of the eastern Pacific and NIO in the prediction of the EAJ location,we build a multiple linear regression model based on the predicted JJA Ni~no3 index and NIOI,to statistically predict the EAJLI.Accordingly,the statistical prediction skill of the EAJLI is due to the linear contribution of the dynamically predicted tropical eastern Pacific and NIO SST anomalies in the MME prediction.The cross-validated time series of the statistically predicted EAJLI are presented in Fig.9.Its correlation coeff i cient with the observation is 0.32,close to the value of 0.35 between the MME prediction and the observation during 1960–2005.Moreover,the cross-validated EAJLI agrees well with the MME-predicted EAJLI,with a correlation coeff i cient of 0.72 between them,indicating that about half of the predicted variance of the EAJLI is contributed linearly by the predicted tropical SST anomalies in the eastern Pacific and NIO in the MME.

    The SST anomalies associated with the summertime EAJLI were also revealed by Lin(2010),particularly in July and August.Lin(2010)proposed that the tropical eastern Pacific warm SST anomaly links to a southward shift of the EAJ through a meridionalteleconnection overthe WNP–EastAsia region,which is associated with increased precipitation over the tropical western Pacific and decreased precipitation over the subtropical WNP.The rainfall anomalies are partially indicated by the EAJLI-related rainfall in the MME prediction, as shown in Fig.7f.On the other hand,Shaman and Tziperman(2007)additionally highlighted the role of the westward propagating Rossby wave,induced by the tropical eastern Pacific warm SST-related divergence forcing over the equatorial Pacific in the uppertroposphere.Thisinteractswith the North Africa–Asia westerly jet and manifests as positive vorticity anomalies along the westerly jet,causing the southward shift of the EAJ.

    Furthermore,Qu and Huang(2012)also found a positive correlation between the warm NIO SST anomalies and southward shift of the EAJ.They proposed that when the tropical Indian Ocean SST is higher than normal,the anomalous tropical convection forces a Kelvin wave wedge penetrating into the equatorial western Pacific,leading to a decrease in precipitation near the Philippines(Xie et al.,2009).Combined with the climatological easterly shear over the subtropical WNP (Lim and Chang,1983,1986;Wang and Xie,1996;Lu,2004; Lin,2010),the EAP teleconnection(Huang and Sun,1992)is induced along the East Asian coast,which accelerates westerly on the southern f l ank of the EAJ and decelerates westerly on the northern f l ank in the upper troposphere(Kosaka and Nakamura,2006;Lu and Lin,2009).Thus,the EAJ shifts southward.

    Furthermore,the anomalous SSTs show a much stronger teleconnection with the 200 hPa zonal winds on the southern side of the EAJ axis than on its northern side.The correlation coeff i cient between the southern side of the EAJ axis and the Ni~no3(NIO)SST anomalies is 0.50(0.47)in the observation.This close teleconnection is reproduced by the model; the corresponding correlation coeff icients are both 0.73 in the MME prediction,presenting more predictable signals of the wind anomalies over the region.However,the correlation coeff i cient between the northern side of the EAJ axis and the Ni~no3(NIO)SST anomalies is only 0.07(?0.20) in the observation and 0.37(0.25)in the MME prediction. The stronger relationship to the southern side of the EAJ axis of zonal winds at 200 hPa with the tropical SST anomalies, leads to higher prediction skill of the southern side compared to the northern side(Table 2).The lower prediction skill on the northern side is probably due to the effect of synoptic transient activities(Xiang and Yang,2012),which have less predictability.

    4.Discussion

    4.1.Predicted decadal change of the EAJ location

    Some previous studies have reported a decadal tendency of the southward shift of the EAJ after the late 1970s(Yu andZhou,2007;Zhang and Huang,2011).This decadal shift in the EAJ is also revealed in the observation and MME prediction of the present study(Fig.6),suggesting a possible predicted decadal change of the EAJ location in ENSEMBLES. To reveal the predictability of decadal change of the EAJ location,we isolate the decadal component of the EAJLI,with its period longer than 8 years,based on a Fourier harmonics fi lter,and the residual is referred to as the interannual component.

    Table 4 shows the correlation coeff icients of the decadal and interannual components of the EAJLI,separately,between the observation and model simulations.The correlation coeff i cient of the EAJLI between the observation and MME,for the decadal component,is 0.47.Similarly,the decadal component is also predicted in the f i ve model predictions,with the correlation coeff icients ranging from 0.26 to 0.40.The potential skill of the decadal variation of the EAJLI is attributed significantly to summer SST anomalies in the tropical central-eastern Pacific and Indian Ocean,with warm tropical SST anomalies related to the southward shift of the EAJ,similar to that shown in Fig.8.Moreover,these tropical SST anomalies can be traced back to an El Ni~no–like SST pattern in the previous winter(figure not shown).The positive SST anomalies may induce atmospheric warming in the tropics(Yulaeva and Wallace,1994;Kumar and Hoerling, 2003;Huang and Huang,2009),which increases the meridional temperature gradient in the subtropics to the south of the EAJ axis and results in a southward shift of the EAJ.

    On the other hand,the potential prediction skill for the interannual component remains stable,compared with the total EAJLI.The correlation coeff i cient of the interannual component of the EAJLI between the MME prediction and the observation,as shown in Table 4,is 0.32,close to the value of 0.35 for the total EAJLI(Table 2),signifi cant at the 95% confi dence level.For the f i ve models,the prediction skill of the interannual component,compared with the total EAJLI, is also stable,with a correlation coeff i cient with the observation of 0.35(0.34)in the UKMO model,0.30(0.31)in the IFM-GEOMAR model,0.27(0.27)in the CMCC-INGV model,0.13(0.19)in the ECMWF model,and 0.0(0.07)in the MF model,for the interannual component(total)of the EAJLI.We also regress circulation and SST anomalies upon the interannual component of the EAJLI,and obtain similar results to those in Figs.7 and 8,except for weakened warm summer SST anomalies in the tropical eastern Pacific in the observation,suggesting a close relationship of the tropical SST anomalies in the eastern Pacific and the EAJ location on decadal timescales,as described in the above paragraph.

    Table 4.Correlation coeff icients of the decadal and interannual component of the EAJLI between the observation and model predictions.

    4.2.Potential prediction skill of the EAJ location

    The interannual variability of midlatitude anomalies is, generally,considered unpredictable(Lu et al.,2006;Wang, 2008;Kosaka et al.,2012).Lu et al.(2006),for example, found that the anomaly of the EAJ—both the meridional shift of the EAJ location and its intensity change—is dominated by atmospheric internal variability,different from the lowertropospheric circulation anomaly over the WNP in which the external variability plays an important role.

    In the present study,we f i nd a correlation coeff i cient of 0.35 ofthe EAJLIduring 1960–2005 between the observation and MME prediction in the previous section.It is signifi cant at the 95%confidence level and better improved compared to that in an AGCM simulation(Lu et al.,2006).However,the potential prediction skill of the EAJLI is still lower than the WNP anomalies,with a correlation coeff i cient of 0.68 for the WNP summer monsoon index(Li et al.,2012)between the MME predictions and the observation.Moreover,the meridional shift of the EAJ location is successfully captured as the EOF1 mode in the f i ve models and their MME in ENSEMBLES,the same as in the observation,while it is represented by the EOF2 mode in the SST-forced AGCM simulation in HadAM3(Lu et al.,2006).The improved potential skill of the EAJ prediction may help increase the prediction skill of climate anomalies in East Asia,since they are closely tied on the interannual timescale(Liang and Wang,1998;Lau et al., 2000;Lu,2004).

    5.Summary

    This study investigates the predictability of the EAJ in summer based on 46-yr predictions initiated on 1 May in the f i ve models of ENSEMBLES during 1960–2005 and NCEP/NCAR reanalysis data(taken as the observation).The results show that the f i ve models can capture the climatological location and intensity of the summer-mean EAJ.In addition,the models also predict strong interannual standard deviation of upper-tropospheric zonal winds over East Asia to the south of the EAJ axis,though its magnitude is weaker than in the observation.

    Moreover,the prediction skill for the first two EOF leading modes,i.e.,the meridional displacement of the EAJ location and intensity change of the EAJ,respectively,is evaluated.The spatial pattern of the EOF1 mode is successfully captured by the f i ve models and their MME.The temporal correlation coeff i cient between the observation and MME is 0.4 for the PC1 and 0.35 for the EAJ location index(EAJLI). Further examination on the two sides of the EAJLI shows that the prediction skill of the EAJLI can be attributed to that in the southern side.For the EOF2 mode,only the ECMWF and IFM-GEOMAR models,among the f i ve,reproduce the spatial pattern of interannual intensity change of the EAJ,but none of them,including their MME,skillfully predict its temporal variation.

    Finally,we compare the circulation and SST anomalies associated with the EAJLI between the observation andmodel predictions.The results show that the models can adequately capture the meridional shift of the Asian westerly jet in the upper troposphere and a meridional wave train over East Asia and the WNP in the lower troposphere,associated with the EAJLI.In addition,the warm SST anomalies in the tropical eastern Pacific and NIO related to the southward shift of the EAJ are also successfully captured,though their teleconnections are stronger in the model predictions than in the observation.Based on a statistical multiple linear regression model built on two model-predicted SST indices,representing the interannual variation of SST in the tropical eastern Pacific and NIO,respectively,we reproduce about half of the year-to-year variance of the model-predicted EAJLI,suggesting that the prediction skill of the EAJLI is contributed to significantly by the linear impact of the model-predicted tropical SST anomalies in the above two regions.The possible mechanisms linking the tropical SST anomalies to the interannual variation of the EAJ location are also discussed.

    Acknowledgements.We thank two anonymous reviewers for their valuable comments.We also thank Dr.LU Riyu from the Institute of Atmospheric Physics,Chinese Academy of Sciences, for discussion during preparation of the manuscript.This research was supported by the National Natural Science Foundation of China (Grant Nos.41375086,41320104007 and 41305067).

    REFERENCES

    Adler,R.F.,and Coauthors,2003:The version-2 global precipitation climatology project(GPCP)monthly precipitation analysis(1979–Present).Journal of Hydrometeorology,4,1147–1167.

    Chowdary,J.S.,S.P.Xie,J.-Y.Lee,Y.Kosaka,and B.Wang, 2010:Predictability of summer Northwest Pacific climate in 11 coupled model hindcasts:Local and remote forcing.J. Geophys.Res.,115,D22121,doi:10.1029/2010JD014595.

    Ding,Q.H.,and B.Wang,2005:Circumglobal teleconnection in the Northern Hemisphere summer.J.Climate,18,3483–3505,doi:10.1175/JCLI3473.1.

    Doblas-Reyes,F.J.,A.Weisheimer,T.N.Palmer,J.M.Murphy, and D.Smith,2010:Forecast quality assessment of the ENSEMBLES seasonal-to-decadal Stream 2 hindcasts.ECMWF Technical Memorandum,No.621,ECMWF,Reading,UK,45 pp.

    Dole,R.M.,and R.X.Black,1990:Life cycles of persistent anomalies.Part II:The development of persistent negative height anomalies over the North Pacific Ocean.Mon.Wea. Rev.,118,824–846.

    Enomoto,T.,2004:Interannual variability of the Bonin high associated with the propagation of Rossby waves along the Asian jet.J.Meteor.Soc.Japan,82,1019–1034.

    Enomoto,T.,B.J.Hoskins,and Y.Matsuda,2003:The formation mechanism of the Bonin high in August.Quart.J.Roy. Meteor.Soc.,129,157–178.

    Huang,G.,2004:An index measuring the interannual variation of the East Asian summer monsoon—The EAP index.Adv. Atmos.Sci.,21,41–52,doi:10.1007/BF02915679.

    Huang,P.,and R.H.Huang,2009:Southern-Northern Hemispheres symmetric and asymmetric effect of El Ni~no events on general circulation and analysis on its mechanism.Chinese J.Atmos.Sci.,33,1–15,doi:10.3878/j.issn.1006-9895.2009.01.01.(in Chinese).

    Huang,R.H.,and F.Y.Sun,1992:Impacts of the tropical western Pacific on the East Asian summer monsoon.J.Meteor.Soc. Japan,70,243–256.

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

    Kosaka,Y.,and H.Nakamura,2006:Structure and dynamics of the summertime Pacific-Japan teleconnection pattern.Quart. J.Roy.Meteor.Soc.,132,2009–2030.

    Kosaka,Y.,J.S.Chowdary,S.P.Xie,Y.-M.Min,and J.-Y.Lee, 2012:Limitations of seasonal predictability for summer climate over East Asia and the Northwestern Pacific.J.Climate, 25,7574–7589,doi:10.1175/JCLI-D-12-00009.1.

    Kumar,A.,and M.P.Hoerling,2003:The nature and causes for the delayed atmospheric response to El Ni~no.J.Climate,16, 1391–1403.

    Lau,K.-M.,K.-M.Kim,and S.Yang,2000:Dynamical and boundary forcing characteristics of regional components of the Asian summer monsoon.J.Climate,13,2461–2482.

    Lee,J.-Y.,B.Wang,Q.Ding,K.-J.Ha,J.-B.Ahn,A.Kumar,B. Stern,and O.Alves,2011:How predictable is the northern hemisphere summer upper-tropospheric circulation?Climate Dyn.,37,1189–1203,doi:10.1007/s00382-010-0909-9.

    Li,C.F.,R.Y.Lu,and B.W.Dong,2012:Predictability of the western North Pacific summer climate demonstrated by the coupled models of ENSEMBLES.Climate Dyn.,39,329–346,doi:10.1007/s00382-011-1274-z.

    Li,C.F.,R.Y.Lu,and B.W.Dong,2014:Predictability of the western North Pacific summer climate associated with different ENSO phases by ENSEMBLES multi-model seasonal forecasts.Climate Dyn.,43,1829–1845,doi:10.1007/ s00382-013-2010-7.

    Liang,X.Z.,and W.-C.Wang,1998:Associations between China monsoon rainfall and tropospheric jets.Quart.J.Roy.Meteor. Soc.,124,2597–2623,doi:10.1002/qj.49712455204.

    Lim,H.,and C.P.Chang,1983:Dynamics of teleconnections and walker circulations forced by equatorial heating.J.Atmos.Sci.,40,1897–1915.

    Lim,H.,and C.P.Chang,1986:Generation of internal-and external-mode motions from internal heating:Effects of vertical shear and damping.J.Atmos.Sci.,43,948–960.

    Lin,Z.D.,2010:Relationship between meridional displacement of the monthly East Asian jet stream in the summer and sea surface temperature in the tropical central and eastern Pacific.Atmos.Oceanic Sci.Lett.,3,40–44.

    Lin,Z.D.,2013:Impacts of two types of northward jumps of the East Asian upper-tropospheric jet stream in midsummer on rainfall in eastern China.Adv.Atmos.Sci.,30,1224–1234, doi:10.1007/s00376-012-2105-9.

    Lin,Z.D.,and R.Y.Lu,2005:Interannual meridional displacement of the East Asian upper-tropospheric jet stream in summer.Adv.Atmos.Sci.,22,199–211,doi:10.1007/BF 02918509.

    Lu,R.Y.,2004:Associations among the components of the East Asian summer monsoon system in the meridional direction.J.Meteor.Soc.Japan,82,155–165.

    Lu,R.Y.,and Z.D.Lin,2009:Role of subtropical precipitation anomalies in maintaining the summertime meridional teleconnection over the western North Pacific and East Asia.J. Climate,22,2058–2072.

    Lu,R.Y.,J.-H.Oh,and B.-J.Klim,2002:A teleconnection pattern in upper-level meridional wind over the North African and Eurasian continent in summer.Tellus,54A,44–55.

    Lu,R.Y.,Y.Li,and B.W.Dong,2006:External and internal summer atmospheric variability in the western North Pacific and East Asia.J.Meteor.Soc.Japan,84,447–462.

    Lu,R.Y.,Z.D.Lin,and Y.C.Zhang,2013:Variability of the East Asian upper-tropospheric jet in summer and its impacts on the East Asian monsoon.Chinese J.Atmos.Sci.,37,331–340, doi:10.3878/j.issn.1006-9895.2012.12310.(in Chinese).

    Nitta,T.,1987:Convective activities in the tropical western Pacific and their impact on the Northern Hemisphere summer circulation.J.Meteor.Soc.Japan,65,373–390.

    Qu,X.,and G.Huang,2012:Impacts of tropical Indian Ocean SST on the meridional displacement of East Asian jet in borealsummer.Inter.J.Climatol.,32,2073–2080,doi:10.1002/ joc.2378.

    Shaman,J.,and E.Tziperman,2007:Summertime ENSO-North African-Asian jet teleconnection and implications for the Indian monsoons.Geophys.Res.Lett.,34,L11702,doi: 10.1029/2006GL029143.

    Smith,T.M.,and R.W.Reynolds,2004:Improved extended reconstruction of SST(1854–1997).J.Climate,17,2466–2477.

    van der Linden,P.,and J.F.B.Mitchell,2009:ENSEMBLES: Climate change and its impact:summary of research and results from ENSEMBLES project.Met Off i ce Hadley Centre, FitzRoy Road,Exeter EX1 3PB,UK,160 pp.

    Wang,B.,2008:Thrusts and prospects on understanding and predicting Asian monsoon climate.Acta Meteorologica Sinica, 22,383–403.

    Wang,B.,and X.S.Xie,1996:Low-frequency equatorial waves in vertically sheared zonal f l ow.Part I:Stable waves.J.Atmos. Sci.,53,449–467.

    Wang,B.,I.-S.Kang,and J.-Y.Lee,2004:Ensemble simulations of Asian-Australian monsoon variability by 11 AGCMs.J. Climate,17,803–818.

    Wang,B.,Q.H.Ding,X.H.Fu,I.-S.Kang,K.Jin,J.Shukla,and F.Doblas-Reyes,2005:Fundamental challenge in simulation and prediction of summer monsoon rainfall.Geophys.Res. Lett.,32,L15711,doi:10.1029/2005GL022734.

    Wang,W.W.,W.Zhou,X.Wang,S.K.Fong,and K.C. Leong,2013:Summer high temperature extremes in Southeast China associated with the East Asian jet stream and circumglobal teleconnection.J.Geophys.Res.,118,8306–8319, doi:10.1002/jgrd.50633.

    Weisheimer,A.,F.J.Doblas-Reyes,T.N.Palmer,A.Alessandri,A.Arribas,M.D′equ′e,N.Keenlyside,M.MacVean, A.Navarra,and P.Rogel,2009:ENSEMBLES:A new multi-model ensemble for seasonal-to-annual predictions—Skill and progress beyond DEMETER in forecasting tropical Pacific SSTs.Geophys.Res.Lett.,36,L21711,doi: 10.1029/2009GL040896.

    Wu,W.J.,J.H.He,H.-S.Chung,C.-H.CHo,and R.Y.Lu,2006: The relationship between the East Asian up-tropospheric jet stream in summer and climatic characteristics of synopticscale disturbance.Climatic and Environmental Research,11, 525–534.(in Chinese).

    Xiang,Y.,and X.Q.Yang,2012:The effect of transient eddy on interannual meridional displacement of summer East Asian subtropical jet.Adv.Atmos.Sci.,29,484–492,doi:10.1007/ s00376-011-1113-5.

    Xie,S.P.,K.M.Hu,J.Hafner,H.Tokinaga,Y.Du,G.Huang, and T.Sampe,2009:Indian Ocean capacitor effect on Indo-Western Pacific climate during the summer following El Ni~no.J.Climate,22,730–747.

    Yu,R.C.,and T.J.Zhou,2007:Seasonality and three-dimensional structure of interdecadal change in the East Asian monsoon.J.Climate,20,5344–5355,doi:10.1175/2007JCLI1559.1.

    Yulaeva,E.,and J.M.Wallace,1994:The signature of ENSO in global temperature and precipitation fields derived from the microwave sounding unit.J.Climate,7,1719–1736.

    Zhang,Y.C.,and L.L.Guo,2005:The deviation of the East Asian subtropical westerly jet and the simulation of seasonal variation of the East China rain bands.Chinese Science Bulletin, 50,1394–1399.(in Chinese).

    Zhang,Y.C.,and D.Q.Huang,2011:Has the East Asian westerly jet experienced a poleward displacement in recent decades?Adv.Atmos.Sci.,28,1259–1265,doi:10.1007/s00376-011-9185-9.

    :Li,C.F.,and Z.D.Lin,2015:Predictability of the summer East Asian upper-tropospheric westerly jet in ENSEMBLES multi-model forecasts.Adv.Atmos.Sci.,32(12),1669–1682,

    10.1007/s00376-015-5057-z.

    16 February 2015;revised 12 May 2015;accepted 10 June 2015)?

    LIN Zhongda Email:zdlin@mail.iap.ac.cn

    韩国高清视频一区二区三区| 成人免费观看视频高清| 亚洲欧洲国产日韩| 久久久久国产网址| 美女xxoo啪啪120秒动态图| 男女啪啪激烈高潮av片| a级片在线免费高清观看视频| 国产69精品久久久久777片| 日韩欧美精品免费久久| 午夜免费鲁丝| 伦精品一区二区三区| 国产无遮挡羞羞视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利网站1000一区二区三区| 国产成人午夜福利电影在线观看| 最近的中文字幕免费完整| 日日爽夜夜爽网站| 青春草亚洲视频在线观看| 国产极品粉嫩免费观看在线 | 亚洲欧美日韩另类电影网站| 久久综合国产亚洲精品| 少妇人妻 视频| 汤姆久久久久久久影院中文字幕| 色94色欧美一区二区| 久久久国产欧美日韩av| 国产精品嫩草影院av在线观看| 日韩不卡一区二区三区视频在线| 久久毛片免费看一区二区三区| av黄色大香蕉| 亚洲国产成人一精品久久久| 男人爽女人下面视频在线观看| av国产久精品久网站免费入址| 国产精品一区二区在线不卡| 只有这里有精品99| 久久久久久久久久人人人人人人| a级片在线免费高清观看视频| 丝袜在线中文字幕| 一本大道久久a久久精品| 9色porny在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产男女超爽视频在线观看| 欧美日韩黄片免| 一区二区av电影网| 91成年电影在线观看| 少妇精品久久久久久久| 欧美日韩精品网址| 久久青草综合色| av片东京热男人的天堂| 国产av精品麻豆| 国产av一区二区精品久久| 欧美 日韩 精品 国产| 搡老熟女国产l中国老女人| 国产极品粉嫩免费观看在线| 国产在线视频一区二区| 精品亚洲乱码少妇综合久久| 悠悠久久av| 一区二区三区激情视频| 成年人午夜在线观看视频| 国产精品98久久久久久宅男小说| 一本久久精品| av免费在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| 久久中文字幕一级| 亚洲九九香蕉| 国产在线观看jvid| 亚洲一区中文字幕在线| 久久精品国产综合久久久| 久久人妻av系列| 伦理电影免费视频| 一区在线观看完整版| 蜜桃国产av成人99| svipshipincom国产片| 欧美精品一区二区免费开放| 国产伦理片在线播放av一区| 黑人巨大精品欧美一区二区蜜桃| 色播在线永久视频| 亚洲情色 制服丝袜| 国产区一区二久久| 亚洲一码二码三码区别大吗| 免费在线观看完整版高清| 丝袜喷水一区| 在线观看免费视频日本深夜| 亚洲av片天天在线观看| 国产97色在线日韩免费| 亚洲精品中文字幕一二三四区 | 热99re8久久精品国产| 后天国语完整版免费观看| 最新的欧美精品一区二区| 中文欧美无线码| 男女边摸边吃奶| 在线永久观看黄色视频| 在线永久观看黄色视频| 亚洲精品美女久久久久99蜜臀| 99re在线观看精品视频| 国产福利在线免费观看视频| 性高湖久久久久久久久免费观看| 久久精品91无色码中文字幕| 一边摸一边抽搐一进一小说 | 男人舔女人的私密视频| 日本五十路高清| videos熟女内射| 久久久精品区二区三区| 亚洲欧美精品综合一区二区三区| 黑人操中国人逼视频| 夫妻午夜视频| 国产成人欧美在线观看 | 精品免费久久久久久久清纯 | 亚洲国产看品久久| 纯流量卡能插随身wifi吗| 欧美国产精品va在线观看不卡| 91老司机精品| 在线播放国产精品三级| 麻豆av在线久日| 中文欧美无线码| 女警被强在线播放| aaaaa片日本免费| 日韩中文字幕视频在线看片| 亚洲欧洲精品一区二区精品久久久| cao死你这个sao货| 考比视频在线观看| 日本精品一区二区三区蜜桃| 一进一出好大好爽视频| 久久久精品免费免费高清| av有码第一页| 一级,二级,三级黄色视频| 十八禁网站免费在线| 日韩免费av在线播放| 国产精品久久久久久精品电影小说| 大片免费播放器 马上看| 我的亚洲天堂| 2018国产大陆天天弄谢| 成人国语在线视频| 99热国产这里只有精品6| 高清毛片免费观看视频网站 | 成人手机av| 精品国产乱码久久久久久男人| 捣出白浆h1v1| 99国产综合亚洲精品| 五月开心婷婷网| 亚洲av日韩精品久久久久久密| 国产精品av久久久久免费| 男女免费视频国产| 一区在线观看完整版| 国产免费福利视频在线观看| 男女免费视频国产| 啪啪无遮挡十八禁网站| 国产免费福利视频在线观看| 91成年电影在线观看| 国产成人精品在线电影| 国产av精品麻豆| 中亚洲国语对白在线视频| 久久人妻福利社区极品人妻图片| 国产男女超爽视频在线观看| 精品少妇黑人巨大在线播放| 国产精品九九99| 岛国毛片在线播放| 精品免费久久久久久久清纯 | 一区福利在线观看| 国产真人三级小视频在线观看| 纯流量卡能插随身wifi吗| 精品视频人人做人人爽| 性色av乱码一区二区三区2| 黄频高清免费视频| 国产国语露脸激情在线看| 久久午夜综合久久蜜桃| 精品福利观看| av天堂在线播放| 十八禁网站免费在线| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器 | 欧美日本中文国产一区发布| 人人妻人人澡人人爽人人夜夜| 久久亚洲真实| 中文字幕人妻熟女乱码| 精品久久蜜臀av无| 国产aⅴ精品一区二区三区波| 99九九在线精品视频| av电影中文网址| 久久久久久久大尺度免费视频| 日韩欧美国产一区二区入口| 18禁国产床啪视频网站| 黑丝袜美女国产一区| 国产精品成人在线| 亚洲欧美一区二区三区久久| 精品久久蜜臀av无| 欧美久久黑人一区二区| 怎么达到女性高潮| 黑人巨大精品欧美一区二区mp4| 欧美日韩av久久| 午夜福利视频精品| 曰老女人黄片| 国产精品99久久99久久久不卡| 又紧又爽又黄一区二区| 黄网站色视频无遮挡免费观看| 国产高清视频在线播放一区| 老司机靠b影院| 女人高潮潮喷娇喘18禁视频| 无遮挡黄片免费观看| 丝袜喷水一区| 精品视频人人做人人爽| 日韩欧美一区视频在线观看| 亚洲第一青青草原| 高清视频免费观看一区二区| 日韩欧美三级三区| 国产欧美日韩精品亚洲av| 亚洲精品美女久久久久99蜜臀| 少妇裸体淫交视频免费看高清 | 国产片内射在线| 欧美久久黑人一区二区| 亚洲专区字幕在线| 亚洲熟女毛片儿| 99久久人妻综合| 韩国精品一区二区三区| 成人av一区二区三区在线看| 黄色 视频免费看| 在线观看一区二区三区激情| 欧美老熟妇乱子伦牲交| 亚洲视频免费观看视频| 丝袜美腿诱惑在线| 国产精品成人在线| 精品熟女少妇八av免费久了| 精品福利永久在线观看| 十八禁人妻一区二区| 一级片'在线观看视频| 水蜜桃什么品种好| 免费在线观看影片大全网站| 精品国产亚洲在线| 日韩免费av在线播放| 午夜福利视频在线观看免费| 久久亚洲精品不卡| www.999成人在线观看| 久久久水蜜桃国产精品网| 国产精品二区激情视频| 99国产精品免费福利视频| 19禁男女啪啪无遮挡网站| 国产有黄有色有爽视频| 国产免费福利视频在线观看| 国产亚洲欧美在线一区二区| 国产欧美亚洲国产| 女人久久www免费人成看片| 丝袜人妻中文字幕| 精品欧美一区二区三区在线| 国产99久久九九免费精品| 黄色视频在线播放观看不卡| 无遮挡黄片免费观看| 日韩视频在线欧美| 国产深夜福利视频在线观看| 欧美 日韩 精品 国产| 亚洲成人免费av在线播放| 免费黄频网站在线观看国产| 国产精品久久久av美女十八| 亚洲熟妇熟女久久| 欧美日韩亚洲综合一区二区三区_| 一进一出好大好爽视频| 亚洲精品中文字幕一二三四区 | 欧美精品av麻豆av| 嫩草影视91久久| avwww免费| 欧美在线一区亚洲| 久9热在线精品视频| 亚洲欧美日韩高清在线视频 | 国产精品免费一区二区三区在线 | 亚洲专区国产一区二区| 丁香六月天网| 亚洲第一av免费看| 天堂动漫精品| 极品教师在线免费播放| 亚洲av国产av综合av卡| 精品国产一区二区久久| 9色porny在线观看| 精品久久久精品久久久| 狠狠精品人妻久久久久久综合| 欧美午夜高清在线| 亚洲国产欧美一区二区综合| 三级毛片av免费| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡| 亚洲国产毛片av蜜桃av| cao死你这个sao货| 国产国语露脸激情在线看| 99久久精品国产亚洲精品| 首页视频小说图片口味搜索| 欧美国产精品一级二级三级| 丰满人妻熟妇乱又伦精品不卡| 99九九在线精品视频| 国产1区2区3区精品| 99国产精品99久久久久| 亚洲va日本ⅴa欧美va伊人久久| 在线观看免费视频日本深夜| 中文字幕制服av| 精品国产超薄肉色丝袜足j| 韩国精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 99国产综合亚洲精品| 亚洲人成电影观看| 日韩中文字幕视频在线看片| 精品国内亚洲2022精品成人 | 免费看十八禁软件| av片东京热男人的天堂| 国产在线精品亚洲第一网站| 我要看黄色一级片免费的| 妹子高潮喷水视频| 久久亚洲精品不卡| 香蕉丝袜av| 国产欧美日韩综合在线一区二区| 国产亚洲午夜精品一区二区久久| 天天躁日日躁夜夜躁夜夜| 在线亚洲精品国产二区图片欧美| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 美女扒开内裤让男人捅视频| 女同久久另类99精品国产91| 好男人电影高清在线观看| 日本黄色视频三级网站网址 | 国产精品国产高清国产av | 黄色视频在线播放观看不卡| 一边摸一边做爽爽视频免费| 久久久国产欧美日韩av| 精品人妻熟女毛片av久久网站| 狠狠狠狠99中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色丝袜av网址大全| 老司机在亚洲福利影院| 妹子高潮喷水视频| 国产精品二区激情视频| 国产精品偷伦视频观看了| 国产精品免费大片| 青青草视频在线视频观看| 激情在线观看视频在线高清 | 97人妻天天添夜夜摸| 无遮挡黄片免费观看| 日韩熟女老妇一区二区性免费视频| 最近最新中文字幕大全免费视频| 精品国产一区二区久久| 色视频在线一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 色在线成人网| 大香蕉久久网| 操美女的视频在线观看| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 人妻一区二区av| 久久久国产成人免费| 久久国产精品人妻蜜桃| 国产亚洲一区二区精品| 欧美黄色片欧美黄色片| 欧美乱妇无乱码| 国产成人av激情在线播放| 久久婷婷成人综合色麻豆| 91麻豆精品激情在线观看国产 | 精品第一国产精品| 欧美精品亚洲一区二区| 在线永久观看黄色视频| 最近最新中文字幕大全免费视频| 国产一区二区在线观看av| 午夜免费鲁丝| 精品国内亚洲2022精品成人 | 日韩大片免费观看网站| 99久久99久久久精品蜜桃| 三上悠亚av全集在线观看| 午夜福利视频在线观看免费| 亚洲av电影在线进入| 国产单亲对白刺激| 岛国在线观看网站| 午夜福利,免费看| 黄片大片在线免费观看| 天堂8中文在线网| 日本黄色视频三级网站网址 | 国产精品香港三级国产av潘金莲| 亚洲精品中文字幕一二三四区 | 国产男靠女视频免费网站| 人成视频在线观看免费观看| 真人做人爱边吃奶动态| 一级a爱视频在线免费观看| 美女午夜性视频免费| 午夜福利在线观看吧| 中文字幕最新亚洲高清| 老司机靠b影院| 精品久久久久久久毛片微露脸| 无人区码免费观看不卡 | 一二三四在线观看免费中文在| tube8黄色片| 交换朋友夫妻互换小说| 欧美中文综合在线视频| 在线观看免费视频日本深夜| 少妇被粗大的猛进出69影院| 中亚洲国语对白在线视频| 最近最新中文字幕大全电影3 | 成人18禁在线播放| 一区二区三区乱码不卡18| kizo精华| 日本wwww免费看| 国产精品久久久久成人av| 深夜精品福利| 亚洲免费av在线视频| 少妇 在线观看| 国产在线视频一区二区| 国产精品亚洲av一区麻豆| 手机成人av网站| 国产伦理片在线播放av一区| 十分钟在线观看高清视频www| 日韩免费高清中文字幕av| 1024香蕉在线观看| 日本一区二区免费在线视频| 女性被躁到高潮视频| 欧美日本中文国产一区发布| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩精品网址| 五月开心婷婷网| 国产精品免费大片| 女人高潮潮喷娇喘18禁视频| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 757午夜福利合集在线观看| 他把我摸到了高潮在线观看 | 美女国产高潮福利片在线看| 丝袜喷水一区| 亚洲精品中文字幕一二三四区 | 在线观看www视频免费| 国产精品九九99| 最黄视频免费看| 黄色成人免费大全| 99精品欧美一区二区三区四区| 操美女的视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | xxxhd国产人妻xxx| 十八禁高潮呻吟视频| 人成视频在线观看免费观看| 黄片播放在线免费| av天堂在线播放| 最近最新免费中文字幕在线| 精品少妇久久久久久888优播| 自线自在国产av| 精品免费久久久久久久清纯 | 操美女的视频在线观看| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 国产单亲对白刺激| 亚洲三区欧美一区| 91成年电影在线观看| 夫妻午夜视频| 男人舔女人的私密视频| 男女下面插进去视频免费观看| 50天的宝宝边吃奶边哭怎么回事| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 精品卡一卡二卡四卡免费| 久久人妻熟女aⅴ| 国产精品二区激情视频| 91国产中文字幕| 色视频在线一区二区三区| 国产精品亚洲一级av第二区| 90打野战视频偷拍视频| 午夜久久久在线观看| 老熟女久久久| 久久精品成人免费网站| 热99国产精品久久久久久7| 亚洲精品美女久久久久99蜜臀| 国产精品免费大片| 天天影视国产精品| bbb黄色大片| 黑人巨大精品欧美一区二区mp4| 免费观看av网站的网址| 国产99久久九九免费精品| 久久亚洲精品不卡| 国产麻豆69| 麻豆乱淫一区二区| 国产免费福利视频在线观看| 中文字幕高清在线视频| 精品一区二区三卡| 亚洲人成电影免费在线| 一级a爱视频在线免费观看| 亚洲国产精品一区二区三区在线| 免费一级毛片在线播放高清视频 | 久久久久久久大尺度免费视频| 欧美在线一区亚洲| 香蕉久久夜色| 99在线人妻在线中文字幕 | 成人国语在线视频| 久久婷婷成人综合色麻豆| 国产免费福利视频在线观看| 丁香六月欧美| 欧美老熟妇乱子伦牲交| 国产精品久久电影中文字幕 | 免费在线观看视频国产中文字幕亚洲| 80岁老熟妇乱子伦牲交| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 女人久久www免费人成看片| 国产片内射在线| 99re在线观看精品视频| 久久中文字幕一级| 久久久久久亚洲精品国产蜜桃av| 在线观看免费高清a一片| 久久久精品免费免费高清| 午夜老司机福利片| 日韩中文字幕欧美一区二区| 国产有黄有色有爽视频| 成人影院久久| 80岁老熟妇乱子伦牲交| 夜夜夜夜夜久久久久| 国产精品久久久久久精品古装| 高清黄色对白视频在线免费看| 窝窝影院91人妻| 日韩精品免费视频一区二区三区| 欧美黄色淫秽网站| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品无人区| 大片免费播放器 马上看| 国产精品免费一区二区三区在线 | 国产免费现黄频在线看| 最近最新中文字幕大全电影3 | 免费女性裸体啪啪无遮挡网站| 脱女人内裤的视频| 国产免费av片在线观看野外av| 亚洲成a人片在线一区二区| 另类精品久久| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| h视频一区二区三区| 99热网站在线观看| 成人特级黄色片久久久久久久 | 国产麻豆69| 国产精品久久久久久精品古装| 夜夜爽天天搞| 久久精品国产综合久久久| 天堂动漫精品| 国产精品欧美亚洲77777| 国产男靠女视频免费网站| 中文字幕人妻熟女乱码| 亚洲熟女毛片儿| 夜夜爽天天搞| 一本大道久久a久久精品| 大香蕉久久网| 淫妇啪啪啪对白视频| 久久毛片免费看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 天天躁夜夜躁狠狠躁躁| 国产精品偷伦视频观看了| 国产精品亚洲一级av第二区| 最近最新中文字幕大全免费视频| 精品高清国产在线一区| 高清视频免费观看一区二区| svipshipincom国产片| 国产在线观看jvid| 在线av久久热| 国产在线一区二区三区精| 一区在线观看完整版| 日韩一区二区三区影片| a级毛片在线看网站| 在线十欧美十亚洲十日本专区| 久久国产精品大桥未久av| 热re99久久国产66热| 19禁男女啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 久久久欧美国产精品| 国产精品久久电影中文字幕 | 搡老熟女国产l中国老女人| 性高湖久久久久久久久免费观看| 真人做人爱边吃奶动态| 国产男女超爽视频在线观看| 高清欧美精品videossex| 亚洲国产精品一区二区三区在线| 欧美日韩av久久| 国产成人精品无人区| 水蜜桃什么品种好| 十八禁人妻一区二区| 黄网站色视频无遮挡免费观看| 久久性视频一级片| 女同久久另类99精品国产91| 国产黄色免费在线视频| 黄色成人免费大全| 这个男人来自地球电影免费观看| 国产视频一区二区在线看| 国产黄色免费在线视频| 九色亚洲精品在线播放| 日韩大码丰满熟妇| 欧美成人午夜精品| tocl精华| 久久人人97超碰香蕉20202| 国产成人一区二区三区免费视频网站| 亚洲国产成人一精品久久久| 免费在线观看黄色视频的| 乱人伦中国视频| 国产成人av教育| 人人妻,人人澡人人爽秒播| 亚洲成a人片在线一区二区| 亚洲av成人不卡在线观看播放网| 99热网站在线观看| 国产一区有黄有色的免费视频| 在线av久久热| 日韩 欧美 亚洲 中文字幕| 亚洲欧洲日产国产| 国产av一区二区精品久久| 成人18禁高潮啪啪吃奶动态图| 久久天堂一区二区三区四区| 伦理电影免费视频| 两个人免费观看高清视频| 黑人猛操日本美女一级片| 欧美老熟妇乱子伦牲交| 国产精品免费大片| 另类亚洲欧美激情| 看免费av毛片| 国产av精品麻豆| 亚洲va日本ⅴa欧美va伊人久久| 亚洲黑人精品在线| 欧美日韩精品网址| 欧美亚洲 丝袜 人妻 在线| 麻豆成人av在线观看| 免费观看a级毛片全部| 欧美+亚洲+日韩+国产| 免费不卡黄色视频| 搡老岳熟女国产| 正在播放国产对白刺激| 99国产精品免费福利视频|