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

    Simulating Eastern-and Central-Pacific Type ENSO Using a Sim p le Coup led M odel

    2018-04-08 10:59:09XianghuiFANGandFeiZHENG
    Advances in Atmospheric Sciences 2018年6期

    Xianghui FANG and Fei ZHENG

    1Institute of Atmospheric Sciences,Fudan University,Shanghai 200438,China

    2International Center for Climate and Environment Science(ICCES),Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    3Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science and Technology,Nanjing 210044,China

    1. Introduction

    As the concept of central Pacific(CP)El Nio began to emerge(Fu and Fletcher,1985;Larkin and Harrison,2005;Ashok et al.,2007;Yu and Kao,2007;Kao and Yu,2009;Kim et al.,2009;Kug et al.,2009),comparisons between this new type and the traditional eastern Pacific(EP)type of ENSO(El Nio–Southern Oscillation)were intensively conducted,including the triggering factors(Vimont et al.,2001,2003;Kug et al.,2009;Lee and McPhaden,2010;Yu and Kim,2011;Ham et al.,2013),the occurrence frequency(Kug et al.,2009;Lee and M cPhaden,2010),the amplitude(e.g.,Zheng et al.,2014),the dynam ic mechanisms(Kug et al.,2009;Yu et al.,2010),the predictabilities(e.g.,Zheng and Yu,2017),and the impacts on the global climate(Kumar et al.,2006;Feng and Li,2011;Zhang et al.,2011).On the subject of the dynamic mechanisms,due to the fact that the anomalous warming center of EP ENSO is located in the EP,the mean thermocline(TH)is quite shallow and permits the perturbations on the subsurface to effectively influenc the sea surface temperature(SST)through upwelling processes,including both TH and Ekman(EK)feedback.In contrast,for CP ENSO,the major warming center is concentrated in the CP,where the TH is relatively deep.This means that the upwelling processes have less of an impact on the SST.However,because the zonal mean SST gradient is strongest in this region—due to the warm pool to the west and the cold tongue to the east—the anomalous zonal-current-related zonal advective(ZA)feedback plays the dominant role.The relative contributions of the TH,EK and ZA processes in the equatorial Pacific have also been widely diagnosed using observational datasets(Latif,1987;K leeman,1993;Kug et al.,2009).

    From another perspective,although ENSO is regarded as a kind of recurring event,with the positive phase named El Nio and the negative phase La Nia,many asymmetries exist between these two phases(Burgers and Stephenson,1999;Jinetal.,2003;Su et al.,2010;Sun and Zhang,2006).Forexample,El Nio is always stronger than La Nia;the anomalous cooling center during La Nia is located relatively farther west than during El Nio;and the decaying period of La Nia is longer than that of El Nio.Although many studies have investigated the differences between the two types of ENSO and the asymmetries between the two phases of ENSO,and have obtained certain results,it is still a great challenge for general circulation models(GCMs)to successfully simulate these ENSO characteristics.

    Because GCMs include many factors that can influenc ENSO events,it is not always easy to quantify the feedback contributions or to attribute the main defects to the correct causes.For this reason,simple models consisting of only the major processes that influenc ENSO are always good candidates for studying specific problems or mechanisms.In this study,a linear air–sea coupled model of the equatorial Pacific is utilized to investigate the importance of the major dynam ic mechanisms(e.g.,the TH,EK and ZA feedbacks)on the ENSO problems mentioned above.

    The rest of this paper is organized as follows:Section 2 gives a brief description of the data and the simple model used in this study.Section 3 describes the modification of the model and its ability to simulate the traditional EP type of ENSO.Section 4 tests the influenc of the ZA feedback on simulating the CP type of ENSO.In Section 5,using a simple nonlinear control scheme,the characteristics of the two types of El Nio are simultaneously simulated,as well as the asymmetries between the two phases of the ENSO events.The results are summarized and discussed in section 6.

    2. Datasets,methods and model

    2.1. Datasets and methods

    Observational and reanalysis datasets are used to analyze the developing phase of El Nio.The wind stress data are from NCEP Reanalysis-2(Kanamitsu et al.,2002);the SST data are from ERSST.v3b(Smith et al.,2008);and the monthly ocean temperature and oceanic circulation data are from GODAS(Behringer and Xue,2004).They are used to quantify the contributions of each oceanic feedback process on the development of the SST and calculate the TH depth(TCD—indicated by the 20°C isotherm depth).The analysis period is from 1980 to 2010.Anomalies presented in this study are calculated by removing the monthly mean climatology.Following Xiang et al.(2013), fiveCPElNio(1991/92,1994/95,2002/03,2004/05and 2009/10)and four EPElNio events(1982/83,1986/87,1997/98 and 2006/07)are chosen.The developing phase for each event is defined as from the month when the oceanic Nio index becomes larger than 0.5°C until the end of the calendar year.For example,the developing phase is from June to December for the 1991/92 El Nio event,and from May to December for the 1997/98 El Nio event.

    To measure the relative contributions of each oceanic feedback on the development of the SST,a straightforward method is used to conduct a budget analysis of the mixed layer temperature.The temperature equation is as follows:

    where overbars and primes represent the monthly climatology and anomaly,respectively.The variablesu,vandTare the zonal and meridional currents and oceanic temperature averaged over the mixed layer.The vertical velocity(w)is measured at the bottom of the mixed layer.Ris the residual term,which consists of thermodynamic processes and influences from variation other than that on the interannual time scale,amongst other factors.This method is the same as that used in Kug et al.(2009).

    Fig.1. (a)Observational linear contributions of ZA(black line and gray shading),TH(blue line and shading),EK(green line and shading),and the combination of TH and EK(red line and shading)feedbacks to the SST tendency along the equatorial Pacific.(b)Factors in the SST equations,Eqs.(2)and(3),as a function of longitude;here,α is in units of K(10 m)?1 month?1,β in K(0.1 Pa)?1 month?1,and γ in month?1.

    Figure 1a illustrates the relative contributions of each process[rhs of Eq.(1)]to the development of local SST tendencies[lhs of Eq.(1)]in the equatorial Pacific for the whole period(1980–2010),rather than for specific El Nio events,as calculated in Kug et al.(2009).It indicates that the ZA feedback(i.e.,the anomalous zonal current and mean zonal tem-perature gradient—black line in the figure)plays the dominant role in the CP region,where the largest mean zonal temperature gradient is located.In the western and eastern areas,the contribution from the ZA feedback is relatively weak.Instead,the TH feedback(i.e.,the mean vertical velocity and anomalous vertical temperature gradient—blue line)and the EK feedback(i.e.,the anomalous vertical velocity and mean vertical temperature gradient—green line)processes are dominant.Specifically,the TH feedback plays an important role in the central–eastern region,i.e.,110°–150°W,with a comparable contribution to the ZA feedback,and the EK feedback plays the dominant role in the eastern part of the Pacific,i.e.,80°–110°W.Because the TH and EK feedbacks are both related to vertical upwelling processes,their combined effects are also illustrated in Fig.1a(red line).It shows that the total vertical upwelling processes can explain the majority of the SST development over the central–eastern Pacific region,with the relative contribution exhibiting a near linear increase from the central to the eastern area.Meanwhile,these feedbacks have nearly no effect upon the far western Pacific(WP),where the mean TH is too deep and the zonal mean SST gradient is too small.The contributions from the other oceanic processes in Eq.(1)are very small,so they are not shown in the figure.Overall,this observational analysis indicates that the ZA feedback(a parabolic shape with maximal strength located in the CP)and the vertical upwelling processes(a near linear increase in strength from the CP to EP)play the dominant roles in the development of the SST in the CP and central–eastern Pacific,respectively.

    2.2. Model

    A linear air–sea coupled model of the equatorial Pacific(GMODEL V3.0)is used in this study(Burgers et al.,2002;Burgers and Van Oldenborgh,2003),which can be obtained at http://www.sciamachy-validation.org/research/CKO/gmodel.htm l.Its oceanic dynamic component is a wind-forced linear shallow water(“1.5 layer”)model of a baroclinic mode on a beta plane.

    The atmospheric model comprises simple bivariate linear regression patterns of observed wind stress anomalies to observed SST anomalies in the Nio3(5°S–5°N,90°–150°W)and Nio4(5°S–5°N,160°E–150°W)regions(Burgers and Van Oldenborgh,2003),based on the reality that these two fields are coupled quite well in the equatorial Pacific region.These regression patterns are then combined with noise and used during the run to calculate the likely atmospheric response to model SST anomalies.Details can be found in Burgers and Van Oldenborgh(2003)and Fang and Zheng(2014).

    Given the response of SST to TH depth and wind stress anomalies,the SST model(Burgers and Van Oldenborgh,2003)is valid only in the central and eastern equatorial Pacific,i.e.,

    where the factors α,β and γ determine the strength of the TH feedback(i.e.,the term proportional to the anomalous TCDh),zonal wind stress feedback(i.e.,the term proportional to the anomalous zonal wind stress τx),and the relaxation term(i.e.,the term proportional to the anomalous SSTT′),which vary considerably along the equator(Fig.1b).It can be seen,in the EP,that the SST tendency is linearly related to the TH depth anomaly,with the coefficient depending on longitude and becoming stronger in the EP(dotted line in Fig.1b).Comparing with the observation in Fig.1a,this distribution bears a strong resemblance to the vertical upwelling processes(red line),which means the TH feedback defined in the model actually contains both the TH and EK feedbacks.In the CP,where the ZA feedback dominates,a term proportional to the anomalous zonal wind stress,but not the zonal current,is used in the SST equation.As stated by Burgers and Van Oldenborgh(2003),the disadvantage of this setting is that the zonal velocity fluctuations that are not related to local zonal wind stress variations are not taken into account.The zonal wind stress feedback is not equivalent to the ZA feedback.Comparing with Fig.1a,the distribution of this coefficient(dashed line in Fig.1b)also bears a strong resemblance to the observational linear contribution of the ZA feedback to the SST development(black line).These comparisons indicate that the simple modeling of the major processes influencing the SST development is quite reasonable.In addition,the strength of the linear damping processes in the SST equation is also shown by the solid line in Fig.1b.

    In the simulations,external noise in the form of wind stress variations forces the oceanic shallow water model and induces the variation of the TH depth.The TH depth and zonal wind stress anomalies influenc the SST through the TH and wind stress feedbacks.The varied SST in turn impacts the zonal wind stress anomalies through the atmospheric model,and so on.Details of the model description and its performance when simulating ENSO can be found in related articles(Burgers et al.,2002;Burgers and Van Oldenborgh,2003;Philip and Van Oldenborgh,2010;Fang and Zheng,2014;Zhang et al.,2015;Zheng et al.,2015).In this study,each experiment conducted in GMODEL lasts for 100 years,and the analyzed period is the final 50 years.

    3. Modification of GMODEL and its performance

    3.1. Modifying GMODEL to GMODEL-ZA

    As mentioned above,many observational analyses have indicated that the ZA feedback is the dominant dynamic process for the development of CP ENSO.However,GMODEL does not directly represent this term,but instead utilizes a socalled wind stress feedback.To better depict the ZA feedback and measure its impact on the ENSO events,the wind stress feedback must be modified into the correct ZA feedback,which linearly correlates with the zonal current anomaly[i.e.,uin Eq.(3)],but not the zonal wind stress anomaly,i.e.,

    The feedback strength distribution remains due to its strong resemblance to the observations.To distinguish from the original GMODEL,the modified GMODEL is termed GMODEL-ZA.

    3.2. Performance of GMODEL-ZA when simulating ENSO events

    As a linear air–sea coupled model,the simulated ENSO pattern is quite unique.For this reason,and for achieving results that are more physically consistent,the method of combined empirical orthogonal function(CEOF)analysis is a straightforward way to exhibit the distribution of the simulated ENSO phenomenon.Figure 2 shows the leading CEOF patterns of the anomalous SST,zonal wind stress(Taux),and TCD of the 51st–100th model years simulated by GMODELZA,which is termed the CTRL run.Also,the observational composites of the anomalous SST,Taux and TCD during the developing phase of EP El Nio are illustrated in the righthand panels.It can be seen that the patterns of the three fields are strongly coupled together and exhibit typical EP ENSO patterns,i.e.,from the perspective of the positive phase,the major warming occurs in the EP,while the cooling part shows a lateral V-shape with the corner located in the west;the strong westerlies are mainly located in the CP and WP,while the weak easterlies are located in the EP and far WP;the TCD pattern is similar to that of the SST,with the TCD deepening in the EP while shallowing in the west.As the TCD feedback is much stronger in the EP than in the WP,the SST variation in the WP is relatively weak,although the local TCD variation is intense.The patterns of TCD and Taux are tightly linked by the Sverdrup balance,i.e.,the balance between the zonal wind stress and the TH depth tilt along the equator(Jin,1997).This indicates GMODEL-ZA can capture the major patterns of EP ENSO successfully,as well as the coupled relationships among the air–sea fields.

    Fig.2. Leading CEOF patterns of the 51st–100th model years simulated by the CTRL run(left-hand panels)and the observational composites of the developing phase of EP El Nio(right-hand panels)for(a,b)SST,(c,d)Taux and(e,f)TCD.The numbers at the top of(a)indicate the percentage of variance explained by the CEOF mode.Contour intervals for(a–f)are 0.4°C,0.4°C,0.01 N m?2,0.004 N m?2,10 m,and 10 m,respectively.Purple dashed lines in each panel are along the equator and dateline.

    Fig.3. Time–longitude diagrams of the monthly SST anomalies in the equatorial Pacific during the 51st–100th(left)and 91st–100th(right)model years simulated by the CTRL run.Purple dashed line in each panel is along the dateline.

    To depict the variation of the ENSO system,Fig.3a shows a time–longitude diagram for monthly SST anomalies along the equator,with the last 10 years(91–100)enlarged in Fig.3b.SST anomalies consistently propagate eastward,with the amplitude small in the west and gradually amplifying from the CP to the east under the combined effects of the ZA and TH feedbacks.This is the typical variation of the traditional EP type of ENSO,as indicated by both delayed oscillation and recharge oscillation.Calculating the power spectrum of the Nio3.4(5°S–5°N,120°–170°W)SST index,the major period of this CTRL run is 4.04 years,which is consistent with observations(Kao and Yu,2009).This experiment indicates that GMODEL-ZA can capture the main characteristics of EP ENSO reasonably well,reflecting that the standard TH feedback is much stronger than the ZA feedback in influencing the SST development.However,being a linear model,it can only simulate a unique pattern,which means it cannot simulate the CP ENSO or the asymmetries between El Nio and La Nia.

    4. Simulating CP ENSO in GMODEL-ZA

    Observational analyses have indicated that the occurrence and development of CP ENSO are mainly concentrated in the CP region,where the TCD is too deep to effectively impact the SST development.Instead,the ZA feedback related to the mean zonal SST gradient and anomalous zonal current plays the dominant role.As a result,if the TH feedback is not weakened in the model,even though the original warming center occurs in the CP,the SST in the EP will grow to be higher than that in the CP during the developing period and will make the event into an EP type eventually,due to the TH feedback being much stronger than the ZA feedback.Thus,to solely simulate the characteristics of CP ENSO,the TH feedback in GMODEL-ZA is switched off[i.e.the first term on the rhs of Eq.(3)is removed].This experiment is termed the CP ENSO run.

    Figure 4 shows the leading CEOF patterns of the anomalous SST,Taux and TCD of the51st–100thmodel years simulated in the CP ENSO run,along with the observational composites during the developing phase of CP El Nio illustrated in the right-hand panels.It indicates that the patterns are totally different from those of the CTRL run.From the perspective of the positive phase,the major warming center is situated in the CP while the weaker warming occurs in the eastern region;the corresponding westerly wind stress anomalies are also in the CP and WP,with a significant westward shift compared with the CTRL run,and the easterlies in the EP have a larger zonal stretch and amplitude.These results are consistent with the observations.The TCD anomalies still show a west-shallow ing–east-deepening pattern,but the positive TCD anomalies show a westward shift and are located in the CP and EP,while negative anomalies are limited to the west and have a narrower zonal stretch.The patterns of TCD and Taux are also tightly linked by the Sverdrup balance,with a consistently westward shift.This indicates that,after sw itching offthe TH feedback and only retaining the ZA feedback and the damping term,the air–sea interaction of the equatorial Pacific can induce the major patterns of CP ENSO reasonably well.

    Fig.4. As in Fig.2 but for the results simulated by the CP ENSO run.

    Figure 5a shows the time–longitude diagram for monthly SST anomalies along the equator,with the last 10 years(91–100)enlarged in Fig.5b.The major SST anomalies are mainly located in the CP,with an eastward propagation from the west to the CP,but no significant eastward propagation from the CP to the east.This manifests the characteristic local development of CP ENSO.It also shows that CP ENSO is much weaker than the EP type,due to the weaker strength of the ZA feedback compared with that of the TH feedback.The power spectrum of the Nio3.4 SST index shows that the major period of this CP ENSO run is concentrated in 2–3 years,with 2.7 years being the most significant.This is also consistent with observations(Kao and Yu,2009).It should be noted that,in the model,due to the TH(ZA)feedback mainly dominating the EP(CP),the Taux center is located in the CP and WP(more to the latter).Therefore,the balance relative to the Taux is maintained by the TCD anomalies between the WP and EP(CP).The system will take a longer(shorter)time to propagate the signal between the two regions via Kelvin and Ross by waves,producing a longer(shorter)period.The relative role of the ZA and TH feedbacks in shaping the characteristics of ENSO(including period and propagation orientation)was also investigated by Zhu et al.(2011).Their study mentioned that the more important the role played by the ZA feedback,the higher frequency oscillations would be observed,which is consistent with the results obtained here.

    This experiment indicates that GMODEL-ZA without TH feedback can capture the main characteristics of CP ENSO reasonably well,confirming the dominant role played by the ZA feedback.Nevertheless,because the model is linear,it can only simulate a unique pattern and cannot simulate EP ENSO or the asymmetries between El Nio and La Nia.

    5. Simulating the two types of ENSO and ENSO asymmetries in GMODEL-ZA

    In the last two sections,the EP and CP types of ENSO were simulated separately by the linear air–sea coupled model GMODEL-ZA,confirming the dominant role played by the TH and ZA feedbacks,respectively.In this section,the two types of ENSO are simulated simultaneously in an integrated setup.In addition,the asymmetries between the two phases of ENSO are taken into consideration.

    Fig.5. As in Fig.3 but for the results simulated by the CP ENSO run.

    Since GMODEL-ZA is a linear model,the positive and negative phases are totally symmetric,i.e.,the characteristics of El Nio and La Nia are the same.In the observations,however,the cooling center in La Nia is closer to the CP and always has a relatively weak amplitude,compared with that of El Nio.Physically,these two characteristics are tightly correlated,i.e.,when the cooling center is closer to the CP,the SST variation there will be influenced by the relatively weaker TH feedback than that in the EP type of ENSO,and the relatively weaker ZA feedback than that in the CP type of ENSO,which makes the amplitude of La Nia weaker than EP El Nio and stronger than CP El Nio,since the TH feedback is stronger than the ZA feedback.To sum up,during the developing period of ENSO,the sole EP(CP)El Nio will be dominated by strong(relatively weak)TH(ZA)feedback in the EP(CP),while La Nia will be influenced by both weaker-than-EP El Nio TH feedback and weaker-than-CP El Nio ZA feedback.

    To simulate the characteristics mentioned above,a simple nonlinear control scheme is added to GMODEL-ZA—namely,letting the model freely integrate at the beginning.When the Nio3 or Nio4 SST index is greater than 0.35°C,the model is sw itched into the El Nio developing scheme.Then,if the Nio3 index is greater than the Nio4 index,the model is sw itched into the EP ENSO scheme,i.e.,the ZA feedback is turned off;and if the Nio3 index is less than the Nio4 index,the model is sw itched into the CP ENSO run scheme,i.e.,the TH feedback is turned off.On the other hand,when the Nio3 and Nio4 SST indices are both less than?0.35°C(since La Nia is difficult to divide into two types),the model is sw itched into the La Nia developing scheme.This scheme includes both the weaker-than-EP El Nio TH feedback(i.e.,multiplying the TH feedback of the EP ENSO scheme by a factor of 0.5)and the weaker-than-CP El Nio ZA feedback(i.e.,multiplying the ZA feedback of the CP ENSO scheme by a factor of 0.5),since the observationally anomalous cooling center during La Nia is located west of the anomalous warming center during EP Nio but east of that during CP Nio.

    As the purpose of this experiment is to simulate the diversities in ENSO,Fig.6 shows the leading two CEOF patterns of the anomalous SST,Taux and TCD of the 51st–100th model years simulated in this run,which is termed the NLCTRL run.The following descriptions of the CEOF patterns are all from the perspective of the positive phase.For the first CEOF patterns,the three fields show an EP El Nio pattern and explain 40.31%of the variation,with the warming SST area located in the EP region;accordingly,the TCD deepens in the EP and shallows in the west,confirming the important role played by the TH feedback over the EP.This indicates that the primary SST variations exist in the EP,and the TCD variations exhibit an east–west seesaw pattern.In addition,the strong westerlies are mainly concentrated in the CP and WP,while the weak easterlies are located in the EP and far WP.For the second CEOF patterns(which explains 30.2%of the variation),the SST exhibits a distribution that is more like the CP type of ENSO,but with a slightly eastward shift;accordingly,the TCD deepens in the CP and shallows on both sides.Compared to the CP type of ENSO(Fig.4c),the TCD of this pattern deepens more significantly in the CP.This is because the TH feedback in this region plays a more important role in this experiment for supporting the develop-ment of La Nia,as compared to that in the sole CP ENSO simulation.In addition,Taux showsa pattern that is more like the CP type,with the westerly wind stress anomalies showing a significant westward shift compared with that of the EP ENSO,and the easterlies in the EP have a larger zonal stretch and amplitude.

    Fig.6. Leading two CEOF patterns of(a,b)SST,(c,d)Taux and(e,f)TCD of the 51st–100th model years simulated by the NLCTRL run.The numbers at the top of(a)and(b)indicate the percentage of variance explained by the CEOF modes.Contour intervals are 0.4,0.01,and 10 for SST,Taux and TCD,respectively.Purple dashed lines in each panel are along the equator and dateline.

    Figure 7 gives the time–longitude diagram for monthly SST anomalies along the equator.As expected,the SST variation patterns capture both the EP and CP El Nio types,with the CP type having a weaker amplitude.Meanwhile,the center and amplitude of some simulated La Nia events are also situated between the CP and EP types of El Nio,which bear a strong resemblance to the observations.The power spectrum of the Nio3.4 SST index shows that the major period of this NLCTRL run is concentrated in 2–5 years,with 4.04 and 2.7 years being the most significant,which is also consistent with the observations.This indicates that,through the simple nonlinear control method,many ENSO characteristics,including the CP and EP type of El Nio and the asymmetries between El Nia and La Nia,can be captured well using the simple linear air–sea coupled GMODEL-ZA.

    6. Summary and discussion

    It has been widely mentioned in observational analyses that the TH and ZA feedbacks play a dom inant role in the development of EP and CP ENSO,respectively.In this paper,a simple linear air–sea coupled model,GMODEL,is modified to correctly describe the ZA process,which is linearly related to the zonal current anomalies.This is then substituted for the original so-called wind stress feedback that is linearly related with the zonal wind stress anomalies in the model.The other settings of GMODEL,including the strength distribution of the TH and ZA feedbacks that bear a strong resemblance to the observations,are retained.The modified GMODEL is termed GMODEL-ZA.The results of the CTRL run mainly re fl ect the characteristics of the typical EP type of ENSO.This,on the one hand,demonstrates that the description of the air–sea interaction in the equatorial Pacific is reasonable;while on the other hand,it re fl ects the fact that the general purpose of previous ENSO models is to convincingly simulate the traditional EP type of ENSO,which is also referred to as canonical ENSO.

    Fig.7. As in Fig.3 but for the results simulated by the NLCTRL run.

    The TH feedback in GMODEL-ZA is switched off to simulate the sole CP type of ENSO.This is done because the TH feedback is much stronger than the ZA feedback in influencing the SST development,which will mean that even though the original warming center occurs in the CP region,the SST in the EP will grow to be higher than that in the CP during the developing period,and the event will eventually develop into an EP type.The results of this CP ENSO run indeed manifest the main characteristics of CP ENSO,both in the anomalous patterns and in the equatorial SST variations.These two experiments con firm the dominant role played by the TH and ZA feedbacks in the development of the EP and CP types of ENSO,respectively.

    However,because the model is linear,it can only simulate aunique ENSO pattern on the specific settings,and induce totally symmetric El Nio and La Nia events.To simulate the characteristics of the two types of El Nio and all the ENSO asymmetries,a simple nonlinear control scheme,which gives different settings under these three circumstances,is added to GMODEL-ZA.The result of this NLCTRL run indeed bears a strong resemblance to the observations,i.e.,it effectively exhibits the characteristics of both the EP and CP types of El Nio,and the center and amplitude of some simulated La Nia events are also situated between the CP and EP types of El Nio.This indicates that,through the simple nonlinear control method,many ENSO characteristics can be captured well using the simple linear air–sea coupled GMODEL-ZA,largely due to the realistic description of the two most important processes over the CP to EP regions,i.e.,the TH and ZA feedbacks,in the model.

    Analyzing the historical simulation results of CM IP5,Fang et al.(2015)indicated that most state-of-the-art models can accurately simulate the main characteristics of EP ENSO,which was regarded as the only ENSO type before the concept of the CP type emerged.This is also consistent with the verification that the extended-range of ENSO had met with encouraging results during the 1990s when the EP type of ENSO prevailed(Kirtman and Schopf,1998;Latif et al.,1998;Chen et al.,2004).However,severe biases exist in models when capturing realistic CP El Nio structures,i.e.,barely any model can simulate the significantly weaker warming anomalies in the EP(when compared to those in the CP).This is also a reason why the ENSO forecast skill during the recent period of 2002–11,when the CP type of ENSO began to occur more frequently,is relatively low compared to that in the 1980s and 1990s(Wang et al.,2010;Barnston et al.,2012;Xue et al.,2013;Zheng et al.,2016).Many observational analyses have demonstrated that the ZA feedback plays the dom inant role in the development of CP ENSO.The experiments conducted using GMODEL-ZA also confi rm that,using a reasonable depiction of the ZA feedback,the model can produce a reasonable simulation.Looking back at the performance of GCM s,nearly all show a common bias when constructing climatological SST distributions,i.e.,the simulated cold tongue stretches much farther west than observed.This leads to a significant equatorial cold SST bias in the central–eastern Pacific,and pushes the largest mean zonal temperature gradient westward,giving it a narrower zonal stretch than observed.Since the ZA process is the combination of the anomalous zonal current and the mean zonal temperature gradient,this model bias in the mean zonal temperature gradient will give an inaccurate depiction of the ZA feedback,and makes it difficult to simulate a realistic CP ENSO.

    Acknowledgements.Theauthorsw ish to thank the two anonymous reviewers for their very helpful comments and suggestions.This work was supported by a project funded by the China Postdoctoral Science Foundation(Grant No.2017M 610225),and the National Natural Science Foundation of China(Grant No.41576019).The author is grateful to Mu MU for his support and comments on the manuscript.The monthly ocean temperature and oceanic circulation data were obtained from http://www.cpc.ncep.noaa.gov/products/GODAS/.

    Ashok,K.,S.K.Behera,S.A.Rao,H.Y.Weng,and T.Yamagata,2007:El Nio Modoki and its possible teleconnection.J.Geophys Res.,112,C11007,https://doi.org/10.1029/2006 JC003798.

    Barnston,A.G.,M.K.Tippett,M.L.L’Heureux,S.H.Li,and D.G.Dew itt,2012:Skillof real-timeseasonal ENSO model predictions during 2002–11:Is our capability increasing?Bull.Amer.Meteor.Soc.,93,631–651,https://doi.org/10.1175/BAMS-D-11-00111.1.

    Behringer,D.,and Y.Xue,2004:Evaluation of the global ocean data assimilation system at NCEP:The Pacific Ocean.Preprints,Eighth Symp.on Integrated Observing and Assimilation Systems for Atmosphere,Oceans,and Land Surface,Seattle,WA,Washington State Convention and Trade Center,Amer.Meteor.Soc.,2.3.[Available online at http://ams.confex.com/ams/pdfpapers/70720.pdf.]

    Burgers,G.,and D.B.Stephenson,1999:The“normality”of El Nio.Geophys.Res.Lett.,26,1027–1030,https://doi.org/10.1029/1999GL900161.

    Burgers,G.,and G.J.Van Oldenborgh,2003:On the impact of local feedbacks in the central Pacific on the ENSO cycle.J.Climate,16,2396–2407,https://doi.org/10.1175/2766.1.

    Burgers,G.,M.A.Balmaseda,F.C.Vossepoel,G.J.Van Oldenborgh,and P.J.Van Leeuwen,2002:Balanced ocean-data assimilation near the equator.J.Phys.Oceanogr.,32,2509–2519,https://doi.org/10.1175/1520-0485-32.9.2509.

    Chen,D.K.,M.A.Cane,A.Kaplan,S.E.Zebiak,and D.J.Huang,2004:Predictability of El Nio over the past 148 years.Nature,428,733–736,https://doi.org/10.1038/nature02439.

    Fang,X.-H.,and F.Zheng,2014:Effect of decadal changes in air-sea interaction on the climate mean state over the tropical Pacific.Atmos.Oceanic Sci.Lett.,7,400–405,https://doi.org/10.3878/j.issn.1674-2834.14.0019.

    Fang,X.-H.,F.Zheng,and J.Zhu,2015:The cloud-radiative effect when simulating strength asymmetry in two types of El Nio events using CM IP5 models.J.Geophys.Res.,120,4357–4369,https://doi.org/10.1002/2014JC010683.

    Feng,J.,and J.P.Li,2011:influenc of El Nio Modoki on spring rainfall over south China.J.Geophys.Res.,116,D13102,https://doi.org/10.1029/2010JD015160.

    Fu,C.B.,and J.Fletcher,1985:Two patterns of equatorial warming associated with El Nio.Science Bulletin,30,1360–1364.

    Ham,Y.G.,J.S.Kug,J.Y.Park,and F.-F.Jin,2013:Sea surface temperature in the north tropical Atlantic as a trigger for El Nio/Southern Oscillation events.Nature Geoscience,6,112–116,https://doi.org/10.1038/ngeo1686.

    Jin,F.-F.,1997:An equatorial ocean recharge paradigm for ENSO.Part I:Conceptual model.J.Atmos.Sci.,54,811–829,https://doi.org/10.1175/1520-0469(1997)054<0811:AEORPF>2.0.CO;2.

    Jin,F.-F.,S.I.An,A.Timmermann,and J.X.Zhao,2003:Strong El Nio events and nonlinear dynamical heating.Geophys.Res.Lett.,30,1120,https://doi.org/10.1029/2002GL016356.

    Kanam itsu,M.,W.Ebisuzaki,J.Woollen,S.-K.Yang,J.J.Hnilo,M.Fiorino,and G.L.Potter,2002:NCEP-DEO AM IP-II Reanalysis(R-2).Bull.Amer.Meteor.Soc.,83,1631–1643,https://doi.org/10.1175/BAMS-83-11-1631.

    Kao,H.Y.,and J.-Y.Yu,2009:Contrasting eastern-Pacific and central-Pacific types of ENSO.J.Climate,22,615–632,https://doi.org/10.1175/2008JCLI2309.1.

    Kim,H.-M.,P.J.Webster,and J.A.Curry,2009:Impact of shifting patterns of Pacific Ocean warming on North Atlantic tropical cyclones.Science,325,77–80,https://doi.org/10.1126/science.1174062.

    Kirtman,B.P.,and P.S.Schopf,1998:Decadal Variability in ENSO Predictability and Prediction.J.Climate,11,2804–2822,https://doi.org/10.1175/1520-0442(1998)011<2804:DVIEPA>2.0.CO;2.

    K leeman,R.,1993:On the dependence of hindcast skill on ocean thermodynamics in a coupled ocean-atmosphere model.J.Climate,6,2012–2033,https://doi.org/10.1175/1520-0442(1993)006<2012:OTDOHS>2.0.CO;2.

    Kug,J.-S.,F.-F.Jin,and S.-I.An,2009:Two types of El Nio events:Cold tongue El Nio and warm pool El Nio.J.Climate,22,1499–1515,https://doi.org/10.1175/2008JCLI2624.1.

    Kumar,K.K.,B.Rajagopalan,M.Hoerling,G.Bates,and M.Cane,2006:Unraveling the mystery of Indian monsoon failure during El Nio.Science,314,115–119,https://doi.org/10.1126/science.1131152.

    Larkin,N.K.,and D.E.Harrison,2005:Global seasonal temperature and precipitation anomalies during El Nio autumn and winter.Geophys.Res.Lett.,32,L16705,https://doi.org/10.1029/2005GL022860.

    Latif,M.,1987:Tropical ocean circulation experiments.J.Phys.Oceanogr.,17,246–263,https://doi.org/10.1175/1520-0485(1987)017<0246:TOCE>2.0.CO;2.

    Latif,M.,and Coauthors,1998:A review of the predictability and prediction of ENSO.J.Geophy.Res.,103,14 375–14 393,https://doi.org/10.1029/97JC03413.

    Lee,T.,and M.J.McPhaden,2010:Increasing intensity of El Nio in the central-equatorial Pacific.Geophys.Res.Lett.,37,L14603,https://doi.org/10.1029/2010GL044007.

    Philip,S.Y.,and G.J.V.Van Oldenborgh,2010:Atmospheric properties of ENSO:Models versus observations.Climate Dyn.,34,1073–1091,https://doi.org/10.1007/s00382-009-0579-7.

    Smith,T.M.,R.W.Reynolds,T.C.Peterson,and J.Law rimore,2008:Improvements to NOAA’s historical merged landocean surface temperature analysis(1880–2006).J.Climate,21,2283–2296,https://doi.org/10.1175/2007JCLI2100.1.

    Su,J.Z.,R.H.Zhang,T.Li,X.Y.Rong,J.-S.Kug,and C.-C.Hong,2010:Causes of the El Nio and La Nia amplitude asymmetry in the equatorial eastern Pacific.J.Climate,23,605–617,https://doi.org/10.1175/2009JCLI2894.1.

    Sun,D.Z.,and T.Zhang,2006:A regulatory effect of ENSO on the time-mean thermal strati fi cation of the equatorial upper ocean.Geophys.Res.Lett.,33,L07710,https://doi.org/10.1029/2005GL025296.

    Vimont,D.J.,D.S.Battisti,and A.C.Hirst,2001:Footprinting:A seasonal connection between the tropics and mid-latitudes.Geophys.Res.Lett.,28,3923–3926,https://doi.org/10.1029/2001GL013435.

    Vimont,D.J.,J.M.Wallace,and D.S.Battisti,2003:The seasonal foot printing mechanism in the Pacific:Implications for ENSO.J.Climate,16,2668–2675,https://doi.org/10.1175/1520-0442(2003)016<2668:TSFM IT>2.0.CO;2.

    Wang,W.Q.,M.Y.Chen,and A.Kumar,2010:An assessment of the CFS real-time seasonal forecasts.Wea.Forecasting,25,950–969,https://doi.org/10.1175/2010WAF2222345.1.

    Xiang,B.Q.,B.Wang,and T.Li,2013:A new paradigm for the predominance of standing Central Pacific warming after the late 1990s.Climate Dyn.,41,327–340,https://doi.org/10.1007/s00382-012-1427-8.

    Xue,Y.,M.Y.Chen,A.Kumar,Z.-Z.Hu,and W.Q.Wang,2013:Prediction skill and bias of tropical Pacific sea surface temperatures in the NCEP Climate Forecast System Version 2.J.Climate,26,5358–5378,https://doi.org/10.1175/JCLI-D-12-00600.1.

    Yu,J.-Y.,and H.-Y.Kao,2007:Decadal changes of ENSO persistence barrier in SST and ocean heat content in dices:1958–2001.J.Geophys.Res.,112,125–138,https://doi.org/10.1029/2006JD007654.

    Yu,J.-Y.,and S.T.Kim,2011:Relationships between extratropical sea level pressure variations and the central Pacific and eastern Pacific types of ENSO.J.Climate,24,708–720.https://doi.org/10.1175/2010JCLI3688.1.

    Yu,J.-Y.,H.-Y.Kao,and T.Lee,2010:Subtropics-related interannual sea surface temperature variability in the central equatorial Pacific.J.Climate,23,2869–2884,https://doi.org/10.1175/2010JCLI3171.1.

    Zhang,W.J.,F.-F.Jin,J.P.Li,and H.-L.Ren,2011:Contrasting impacts of two-type El Nio over the western north Pacific during boreal autumn.J.Meteor.Soc.of Japan,89,563–569,https://doi.org/10.2151/jmsj.2011-510.

    Zhang,W.,Q.-L.Chen,and F.Zheng,2015:Bias corrections of the heat flux damping process to improve thesimulation of ENSO post-2000.SOLA,11,181–185,https://doi.org/10.2151/sola.2015-040.

    Zheng,F.,and J.-Y.Yu,2017:Contrasting the skills and biases of deterministic predictions for the two types of El Nio.Adv.Atmos.Sci.,34(12),1395–1403,https://doi.org/10.1007/s00376-017-6324-y.

    Zheng,F.,X.-H.Fang,J.-Y.Yu,and J.Zhu,2014:Asymmetry of the Bjerknes positive feedback between the two types of El Nio.Geophys.Res.Lett.,41,7651–7657,https://doi.org/10.1002/2014GL062125.

    Zheng,F.,W.Zhang,J.-Y.Yu,and Q.-L.Chen,2015:A possible bias of simulating the post-2000 changing ENSO.Science Bulletin,60(21),1850–1857,https://doi.org/10.1007/s11434-015-0912-y.

    Zheng,F.,X.-H.Fang,J.Zhu,J.-Y.Yu,and X.-C.Li,2016:Modulation of Bjerknes feedback on the decadal variations in ENSO predictability.Geophys.Res.Lett.,43,12 560–12 568,https://doi.org/10.1002/2016GL071636.

    Zhu,J.,G.Zhou,R.-H.Zhang,and Z.Sun,2011:On the role of oceanic entrainment temperature(Te)in decadal changes of El Nio/Southern Oscillation.Annales Geophysicae,29(3),529–540,https://doi.org/10.5194/angeo-29-529-2011.

    日韩av不卡免费在线播放| 一夜夜www| 黄片wwwwww| 日本午夜av视频| 国产精品人妻久久久影院| av国产免费在线观看| 我的老师免费观看完整版| 亚洲av日韩在线播放| 天堂中文最新版在线下载 | 大香蕉久久网| 久久人人爽人人片av| 亚洲欧美清纯卡通| 啦啦啦啦在线视频资源| 久久国内精品自在自线图片| 一级毛片久久久久久久久女| 大片免费播放器 马上看| 午夜福利视频精品| 精品久久久精品久久久| 国产免费视频播放在线视频 | 国产成人精品福利久久| 欧美最新免费一区二区三区| 午夜福利在线观看免费完整高清在| 国产乱人视频| 国产精品国产三级专区第一集| 深夜a级毛片| 亚洲熟女精品中文字幕| 久久精品夜夜夜夜夜久久蜜豆| av免费在线看不卡| 一级毛片电影观看| 国产精品女同一区二区软件| 久久99热这里只频精品6学生| 久久久久久国产a免费观看| 最近视频中文字幕2019在线8| 最新中文字幕久久久久| 97超视频在线观看视频| 一边亲一边摸免费视频| or卡值多少钱| 美女cb高潮喷水在线观看| 亚洲最大成人手机在线| 熟女人妻精品中文字幕| 日本欧美国产在线视频| 国产精品蜜桃在线观看| 天堂影院成人在线观看| 国产午夜精品论理片| 欧美最新免费一区二区三区| 亚洲精品视频女| 成人国产麻豆网| 99热网站在线观看| 1000部很黄的大片| 纵有疾风起免费观看全集完整版 | 我的老师免费观看完整版| 熟妇人妻久久中文字幕3abv| av专区在线播放| 国产成人精品一,二区| 成年版毛片免费区| 国产av码专区亚洲av| 91午夜精品亚洲一区二区三区| 成人毛片a级毛片在线播放| 日韩制服骚丝袜av| 又粗又硬又长又爽又黄的视频| 日本免费a在线| 春色校园在线视频观看| av女优亚洲男人天堂| 亚洲av电影不卡..在线观看| 中文字幕久久专区| 欧美高清成人免费视频www| 少妇的逼水好多| 亚洲欧美日韩东京热| 国产国拍精品亚洲av在线观看| 我的老师免费观看完整版| 午夜福利网站1000一区二区三区| 日本与韩国留学比较| 一级毛片 在线播放| 欧美性猛交╳xxx乱大交人| 成人无遮挡网站| 99热全是精品| 国产在视频线在精品| 日韩在线高清观看一区二区三区| 国产高清不卡午夜福利| 亚洲av国产av综合av卡| 精品久久久久久电影网| 亚洲av电影在线观看一区二区三区 | 国产在线一区二区三区精| 在线播放无遮挡| 一个人看视频在线观看www免费| av播播在线观看一区| 午夜亚洲福利在线播放| 免费看av在线观看网站| 大又大粗又爽又黄少妇毛片口| 综合色av麻豆| 国产成人精品婷婷| 国精品久久久久久国模美| 午夜福利视频1000在线观看| 亚洲精品中文字幕在线视频 | 日本黄大片高清| 精品久久久久久久久久久久久| 亚洲精品一区蜜桃| 晚上一个人看的免费电影| 国产在线男女| 乱系列少妇在线播放| 一区二区三区免费毛片| 久久久久精品性色| 韩国高清视频一区二区三区| 国产一区有黄有色的免费视频 | 国产乱人视频| av在线蜜桃| 男人爽女人下面视频在线观看| 一本一本综合久久| 国产色爽女视频免费观看| 国产久久久一区二区三区| 99久久中文字幕三级久久日本| 激情五月婷婷亚洲| 国产亚洲午夜精品一区二区久久 | 亚洲av.av天堂| 五月玫瑰六月丁香| 亚洲国产欧美人成| 大香蕉97超碰在线| 一级毛片我不卡| 秋霞在线观看毛片| 精品99又大又爽又粗少妇毛片| 国产在线一区二区三区精| 国产成人一区二区在线| 精品熟女少妇av免费看| 精华霜和精华液先用哪个| 久久久国产一区二区| 美女xxoo啪啪120秒动态图| 精品久久久精品久久久| 亚洲成人久久爱视频| 亚洲国产精品国产精品| 男女那种视频在线观看| 久久久久免费精品人妻一区二区| 日本熟妇午夜| 中文字幕av成人在线电影| 久久久久久伊人网av| 亚洲国产欧美在线一区| 国产一区亚洲一区在线观看| 久久久久久国产a免费观看| 国产探花在线观看一区二区| 亚洲内射少妇av| 久久97久久精品| 看非洲黑人一级黄片| 国产av在哪里看| 少妇人妻一区二区三区视频| 一级片'在线观看视频| 欧美激情在线99| 日韩 亚洲 欧美在线| 日韩一区二区三区影片| 日韩视频在线欧美| 搞女人的毛片| 天天躁夜夜躁狠狠久久av| 美女国产视频在线观看| 国产精品精品国产色婷婷| 在现免费观看毛片| 亚洲一区高清亚洲精品| 一级毛片aaaaaa免费看小| a级一级毛片免费在线观看| 国产熟女欧美一区二区| av播播在线观看一区| 日韩av免费高清视频| 亚洲精品自拍成人| 国内揄拍国产精品人妻在线| 高清av免费在线| 午夜福利视频精品| 中文字幕制服av| 亚洲无线观看免费| 亚洲av中文字字幕乱码综合| 极品教师在线视频| 亚洲精品一二三| 久久99蜜桃精品久久| 天美传媒精品一区二区| 精品久久久噜噜| 日产精品乱码卡一卡2卡三| 舔av片在线| 国产单亲对白刺激| 日本黄大片高清| 激情 狠狠 欧美| 晚上一个人看的免费电影| 九九爱精品视频在线观看| 亚洲经典国产精华液单| 国产亚洲最大av| 欧美成人精品欧美一级黄| 在线a可以看的网站| 亚洲国产欧美在线一区| 国产91av在线免费观看| 亚洲精品乱码久久久久久按摩| 偷拍熟女少妇极品色| 欧美成人精品欧美一级黄| 精品国产一区二区三区久久久樱花 | 亚洲自拍偷在线| 搡老乐熟女国产| 插阴视频在线观看视频| 尾随美女入室| 99久久中文字幕三级久久日本| 简卡轻食公司| 久久鲁丝午夜福利片| 亚洲国产精品专区欧美| 成年免费大片在线观看| 国产精品久久久久久久电影| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 黄片wwwwww| 18禁裸乳无遮挡免费网站照片| 国产淫语在线视频| 免费播放大片免费观看视频在线观看| 午夜爱爱视频在线播放| 成年版毛片免费区| 97在线视频观看| 国产单亲对白刺激| 成人毛片a级毛片在线播放| 欧美性猛交╳xxx乱大交人| 亚洲美女搞黄在线观看| 亚洲精品亚洲一区二区| 亚洲va在线va天堂va国产| 久久久精品免费免费高清| 久久精品熟女亚洲av麻豆精品 | 久久久久性生活片| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩卡通动漫| 亚洲精品久久午夜乱码| 99久国产av精品国产电影| 蜜桃亚洲精品一区二区三区| 色尼玛亚洲综合影院| 超碰97精品在线观看| 久久久久网色| 国产精品嫩草影院av在线观看| 国产精品一区www在线观看| 久久久久性生活片| 精品亚洲乱码少妇综合久久| 内地一区二区视频在线| 午夜激情福利司机影院| 男人舔奶头视频| 少妇高潮的动态图| 亚洲性久久影院| 91精品一卡2卡3卡4卡| 极品少妇高潮喷水抽搐| 国产免费又黄又爽又色| 亚洲精品视频女| 在线a可以看的网站| 人人妻人人澡人人爽人人夜夜 | 一级爰片在线观看| 观看免费一级毛片| 国产在视频线精品| 99久久中文字幕三级久久日本| 搡女人真爽免费视频火全软件| 日本猛色少妇xxxxx猛交久久| 色哟哟·www| 嫩草影院精品99| av在线天堂中文字幕| 毛片女人毛片| 亚洲国产精品国产精品| 亚洲无线观看免费| 午夜福利高清视频| 久久久久久伊人网av| 日日摸夜夜添夜夜爱| 午夜激情欧美在线| 成年女人在线观看亚洲视频 | 国产色婷婷99| 我要看日韩黄色一级片| 毛片女人毛片| 日本三级黄在线观看| 国产成人freesex在线| 男女啪啪激烈高潮av片| 亚洲精品色激情综合| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看成人毛片| 一级爰片在线观看| 精品不卡国产一区二区三区| 日韩欧美精品v在线| 97在线视频观看| 亚洲精品日本国产第一区| 午夜免费激情av| 人妻系列 视频| 精品国产三级普通话版| 天美传媒精品一区二区| 免费看av在线观看网站| 久久久久久久久中文| 欧美激情在线99| 久久久国产一区二区| 日本av手机在线免费观看| 一本一本综合久久| 久久精品夜色国产| 久久国内精品自在自线图片| 高清欧美精品videossex| 狂野欧美白嫩少妇大欣赏| 九九爱精品视频在线观看| 黄色欧美视频在线观看| 亚洲欧美日韩东京热| 中文天堂在线官网| 婷婷六月久久综合丁香| 青春草亚洲视频在线观看| 最近手机中文字幕大全| 国产精品一二三区在线看| 精品99又大又爽又粗少妇毛片| 亚洲在久久综合| 国产一区二区三区av在线| 日本一二三区视频观看| av天堂中文字幕网| 国产成人a∨麻豆精品| 在线播放无遮挡| 亚洲自偷自拍三级| 国产色爽女视频免费观看| eeuss影院久久| 少妇熟女欧美另类| 亚洲精品视频女| av女优亚洲男人天堂| 国产精品久久久久久久电影| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 免费观看性生交大片5| 国产又色又爽无遮挡免| 成人特级av手机在线观看| 尤物成人国产欧美一区二区三区| 色综合亚洲欧美另类图片| 午夜福利成人在线免费观看| av.在线天堂| 亚洲av男天堂| 亚洲欧美中文字幕日韩二区| kizo精华| 日韩精品有码人妻一区| 国产亚洲精品久久久com| 精品久久久久久久人妻蜜臀av| 少妇人妻一区二区三区视频| 日韩欧美精品免费久久| 成人一区二区视频在线观看| 亚洲成人精品中文字幕电影| 美女内射精品一级片tv| 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 一本久久精品| 97超碰精品成人国产| 日本wwww免费看| 人妻少妇偷人精品九色| 国产男女超爽视频在线观看| 丰满乱子伦码专区| 午夜久久久久精精品| 男人舔女人下体高潮全视频| 中文字幕免费在线视频6| 97人妻精品一区二区三区麻豆| 国产精品一二三区在线看| 国产黄色免费在线视频| 一级黄片播放器| 久久久久久九九精品二区国产| 毛片女人毛片| 1000部很黄的大片| 亚洲成人久久爱视频| 精品一区二区三区人妻视频| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 国产 亚洲一区二区三区 | 日韩 亚洲 欧美在线| 一区二区三区高清视频在线| 久久久久精品性色| 日韩av在线免费看完整版不卡| 国产人妻一区二区三区在| 国产一级毛片七仙女欲春2| 春色校园在线视频观看| 国产一级毛片七仙女欲春2| 国产亚洲精品久久久com| 国产精品久久久久久久电影| 国产精品熟女久久久久浪| 你懂的网址亚洲精品在线观看| 久久午夜福利片| 亚洲成色77777| 欧美成人a在线观看| 国产黄片美女视频| videossex国产| 国产男女超爽视频在线观看| 男插女下体视频免费在线播放| 欧美日韩综合久久久久久| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 国产亚洲精品久久久com| 国产午夜福利久久久久久| 插阴视频在线观看视频| 国产精品一区二区性色av| 老司机影院成人| 大话2 男鬼变身卡| 欧美xxxx性猛交bbbb| 日韩av不卡免费在线播放| 国产麻豆成人av免费视频| 久久久久久九九精品二区国产| 熟女人妻精品中文字幕| av在线观看视频网站免费| 人妻系列 视频| 国产黄频视频在线观看| 久久草成人影院| 久久国产乱子免费精品| 国产毛片a区久久久久| 偷拍熟女少妇极品色| 特级一级黄色大片| 亚洲aⅴ乱码一区二区在线播放| 女人久久www免费人成看片| 一个人观看的视频www高清免费观看| 大香蕉久久网| 一个人观看的视频www高清免费观看| 黄色配什么色好看| 午夜福利在线观看吧| 99视频精品全部免费 在线| 日本一二三区视频观看| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 亚洲乱码一区二区免费版| 亚洲精品中文字幕在线视频 | 成人av在线播放网站| 我的老师免费观看完整版| 精品人妻熟女av久视频| 2021天堂中文幕一二区在线观| 久久精品久久久久久久性| 久久国产乱子免费精品| 亚洲图色成人| 免费av毛片视频| 真实男女啪啪啪动态图| 精品久久久久久久久av| 国产伦理片在线播放av一区| 91久久精品国产一区二区成人| 久久鲁丝午夜福利片| 嫩草影院入口| 天堂中文最新版在线下载 | 亚洲精品在线美女| 日韩精品有码人妻一区| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 午夜福利一区二区在线看| 久久热在线av| 精品酒店卫生间| 国产免费一区二区三区四区乱码| 亚洲国产精品一区二区三区在线| 一二三四在线观看免费中文在| 汤姆久久久久久久影院中文字幕| 久久久久久久大尺度免费视频| 国产欧美日韩综合在线一区二区| 男女午夜视频在线观看| 2018国产大陆天天弄谢| 国产不卡av网站在线观看| 校园人妻丝袜中文字幕| 色网站视频免费| 色婷婷久久久亚洲欧美| 国产欧美亚洲国产| 亚洲欧美一区二区三区久久| 国产女主播在线喷水免费视频网站| 最近2019中文字幕mv第一页| 大香蕉久久成人网| 国产深夜福利视频在线观看| 中文精品一卡2卡3卡4更新| 在线观看免费视频网站a站| 久久影院123| 岛国毛片在线播放| 最近手机中文字幕大全| 1024视频免费在线观看| 中文字幕制服av| 色吧在线观看| 日本欧美视频一区| 精品第一国产精品| 国产成人精品在线电影| av女优亚洲男人天堂| 免费人妻精品一区二区三区视频| 飞空精品影院首页| 中国三级夫妇交换| 久久这里有精品视频免费| 久久综合国产亚洲精品| 国产高清不卡午夜福利| 下体分泌物呈黄色| 波多野结衣一区麻豆| 丝瓜视频免费看黄片| 久久久久视频综合| 熟女av电影| 亚洲色图综合在线观看| 久久人妻熟女aⅴ| 黄色 视频免费看| 欧美日韩精品成人综合77777| 亚洲av在线观看美女高潮| 自拍欧美九色日韩亚洲蝌蚪91| 80岁老熟妇乱子伦牲交| 日韩在线高清观看一区二区三区| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜制服| 秋霞在线观看毛片| 搡女人真爽免费视频火全软件| 久久国产精品男人的天堂亚洲| 国产人伦9x9x在线观看 | 亚洲av免费高清在线观看| freevideosex欧美| 狠狠精品人妻久久久久久综合| 久久久久视频综合| 国产毛片在线视频| 人妻人人澡人人爽人人| 最新的欧美精品一区二区| 18禁动态无遮挡网站| 日日撸夜夜添| 热99久久久久精品小说推荐| 亚洲成人手机| 边亲边吃奶的免费视频| 晚上一个人看的免费电影| 中国三级夫妇交换| 欧美激情 高清一区二区三区| 日日爽夜夜爽网站| 大话2 男鬼变身卡| 欧美亚洲日本最大视频资源| 色吧在线观看| 亚洲第一青青草原| 大片免费播放器 马上看| 欧美精品一区二区大全| 婷婷成人精品国产| 亚洲精品乱久久久久久| 欧美老熟妇乱子伦牲交| 亚洲欧美一区二区三区久久| 国产黄色免费在线视频| 亚洲伊人色综图| 国产精品欧美亚洲77777| www.av在线官网国产| 老司机亚洲免费影院| 91aial.com中文字幕在线观看| 丁香六月天网| 国产一区二区 视频在线| 免费高清在线观看日韩| 精品少妇一区二区三区视频日本电影 | 人体艺术视频欧美日本| 欧美成人午夜精品| 国产精品秋霞免费鲁丝片| 一区二区日韩欧美中文字幕| 免费观看无遮挡的男女| 永久免费av网站大全| 如何舔出高潮| 一边亲一边摸免费视频| 精品酒店卫生间| 欧美日韩av久久| 日韩伦理黄色片| 青春草视频在线免费观看| 日韩成人av中文字幕在线观看| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 亚洲天堂av无毛| 国产日韩欧美在线精品| 在线观看美女被高潮喷水网站| 精品国产露脸久久av麻豆| 免费黄色在线免费观看| 亚洲美女视频黄频| 国产精品国产av在线观看| 王馨瑶露胸无遮挡在线观看| 午夜日本视频在线| av在线app专区| 丰满迷人的少妇在线观看| 国产有黄有色有爽视频| 久久久精品区二区三区| 亚洲成国产人片在线观看| 在线精品无人区一区二区三| 女性被躁到高潮视频| 欧美人与性动交α欧美精品济南到 | 色网站视频免费| 一边摸一边做爽爽视频免费| 国产97色在线日韩免费| 人妻一区二区av| 春色校园在线视频观看| 成人黄色视频免费在线看| 国产日韩欧美在线精品| 亚洲,一卡二卡三卡| 亚洲精品aⅴ在线观看| 18禁动态无遮挡网站| 亚洲欧洲日产国产| 久久国产精品大桥未久av| 麻豆精品久久久久久蜜桃| 七月丁香在线播放| 久久av网站| 中文字幕精品免费在线观看视频| 国产精品无大码| 极品少妇高潮喷水抽搐| 天堂8中文在线网| 国产成人精品久久久久久| 久久久久久久大尺度免费视频| √禁漫天堂资源中文www| 免费高清在线观看日韩| 亚洲精品在线美女| 大片电影免费在线观看免费| videosex国产| 亚洲三区欧美一区| 一本—道久久a久久精品蜜桃钙片| 尾随美女入室| 国产老妇伦熟女老妇高清| 人人妻人人澡人人看| 1024香蕉在线观看| 不卡av一区二区三区| 婷婷色综合大香蕉| 两个人免费观看高清视频| 丝袜美足系列| 免费高清在线观看视频在线观看| 国产精品久久久av美女十八| 一区二区三区四区激情视频| 亚洲av欧美aⅴ国产| 精品国产超薄肉色丝袜足j| 亚洲av成人精品一二三区| 国产精品香港三级国产av潘金莲 | 免费在线观看完整版高清| 欧美人与性动交α欧美软件| 国产av精品麻豆| 美女视频免费永久观看网站| 亚洲欧洲日产国产| 精品国产露脸久久av麻豆| 飞空精品影院首页| 欧美少妇被猛烈插入视频| 不卡视频在线观看欧美| 国产黄色视频一区二区在线观看| 日本vs欧美在线观看视频| 在线精品无人区一区二区三| 国产 精品1| 亚洲综合色惰| 捣出白浆h1v1| 侵犯人妻中文字幕一二三四区| 少妇的丰满在线观看| 成人免费观看视频高清| 秋霞在线观看毛片| 一本色道久久久久久精品综合| 国产在线免费精品| 黄网站色视频无遮挡免费观看| av网站在线播放免费| 大陆偷拍与自拍|