• <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.

    亚洲久久久久久中文字幕| 亚洲人成网站高清观看| 欧美最新免费一区二区三区| 日本与韩国留学比较| 一进一出好大好爽视频| 欧美黑人欧美精品刺激| videossex国产| 很黄的视频免费| 国产成人a区在线观看| 波野结衣二区三区在线| 久久婷婷人人爽人人干人人爱| 国内毛片毛片毛片毛片毛片| 国产淫片久久久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 亚洲午夜理论影院| 国产爱豆传媒在线观看| 在线a可以看的网站| 日本免费一区二区三区高清不卡| 99精品在免费线老司机午夜| or卡值多少钱| 日韩欧美精品v在线| 老熟妇仑乱视频hdxx| eeuss影院久久| 亚洲成人免费电影在线观看| 2021天堂中文幕一二区在线观| 欧美最新免费一区二区三区| 亚洲欧美激情综合另类| 999久久久精品免费观看国产| 成人无遮挡网站| 99视频精品全部免费 在线| 免费大片18禁| 美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 中国美女看黄片| 99国产精品一区二区蜜桃av| 禁无遮挡网站| 男人的好看免费观看在线视频| 欧美丝袜亚洲另类 | 网址你懂的国产日韩在线| 国产精品久久久久久久电影| 美女免费视频网站| a级毛片a级免费在线| videossex国产| 两人在一起打扑克的视频| 久久久精品大字幕| 亚洲电影在线观看av| 欧美bdsm另类| 校园人妻丝袜中文字幕| 精品欧美国产一区二区三| 欧美极品一区二区三区四区| 伊人久久精品亚洲午夜| 色在线成人网| 国产精品女同一区二区软件 | 精品人妻1区二区| 97超视频在线观看视频| 亚洲精品久久国产高清桃花| 一进一出抽搐gif免费好疼| 国产精品久久久久久亚洲av鲁大| 九色成人免费人妻av| 欧美黑人巨大hd| 国产美女午夜福利| 美女被艹到高潮喷水动态| 在线天堂最新版资源| 日韩中字成人| 国产精品爽爽va在线观看网站| 久久国内精品自在自线图片| 久久国产精品人妻蜜桃| 在线天堂最新版资源| 日韩精品有码人妻一区| 一区二区三区免费毛片| 国内久久婷婷六月综合欲色啪| 欧美日本视频| 三级国产精品欧美在线观看| 国产精品一区二区三区四区久久| 国内精品宾馆在线| 人妻丰满熟妇av一区二区三区| 欧美不卡视频在线免费观看| 日韩欧美国产一区二区入口| 欧美性猛交黑人性爽| 男插女下体视频免费在线播放| 午夜福利成人在线免费观看| 男人狂女人下面高潮的视频| 日本色播在线视频| а√天堂www在线а√下载| 色精品久久人妻99蜜桃| 国产一区二区三区视频了| 人妻丰满熟妇av一区二区三区| 男人和女人高潮做爰伦理| 男女做爰动态图高潮gif福利片| 免费看光身美女| 又爽又黄无遮挡网站| 国产精品久久久久久av不卡| 精品久久久久久久久亚洲 | 日日夜夜操网爽| 69av精品久久久久久| 尾随美女入室| 毛片一级片免费看久久久久 | 日韩欧美三级三区| 中文字幕av在线有码专区| 亚洲五月天丁香| 亚洲真实伦在线观看| 精品一区二区三区av网在线观看| 午夜激情欧美在线| 真人一进一出gif抽搐免费| 一个人免费在线观看电影| 亚洲无线观看免费| 男女做爰动态图高潮gif福利片| 特级一级黄色大片| 波多野结衣高清无吗| 国国产精品蜜臀av免费| 日本撒尿小便嘘嘘汇集6| 国产精品人妻久久久久久| 免费看a级黄色片| 日本免费a在线| 性欧美人与动物交配| 亚洲,欧美,日韩| 黄片wwwwww| 夜夜爽天天搞| 日韩大尺度精品在线看网址| 婷婷精品国产亚洲av在线| 舔av片在线| 啦啦啦啦在线视频资源| 久久香蕉精品热| 男女下面进入的视频免费午夜| 免费不卡的大黄色大毛片视频在线观看 | 日本 欧美在线| 国产精品久久久久久av不卡| 欧美成人性av电影在线观看| 亚洲精品久久国产高清桃花| 毛片女人毛片| 99国产精品一区二区蜜桃av| 啪啪无遮挡十八禁网站| 欧美成人性av电影在线观看| 夜夜爽天天搞| 精品日产1卡2卡| 日本熟妇午夜| 久久中文看片网| 日韩亚洲欧美综合| 免费av观看视频| 国产免费一级a男人的天堂| 久久人人爽人人爽人人片va| 97超级碰碰碰精品色视频在线观看| 欧美潮喷喷水| 2021天堂中文幕一二区在线观| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 在线观看av片永久免费下载| 在线观看美女被高潮喷水网站| 中文字幕免费在线视频6| 婷婷精品国产亚洲av| 国产精品久久久久久精品电影| 亚洲va在线va天堂va国产| 国产成人a区在线观看| 看黄色毛片网站| 一区二区三区激情视频| 日日摸夜夜添夜夜添av毛片 | 91在线观看av| 成年版毛片免费区| 欧美三级亚洲精品| av中文乱码字幕在线| 日韩欧美一区二区三区在线观看| 久久久久久大精品| 美女高潮的动态| 天天一区二区日本电影三级| 五月玫瑰六月丁香| 99热6这里只有精品| 中文在线观看免费www的网站| 久久久久久久久久成人| 内射极品少妇av片p| 变态另类成人亚洲欧美熟女| 国产一级毛片七仙女欲春2| 精品人妻1区二区| 国产伦精品一区二区三区视频9| 亚洲国产欧洲综合997久久,| 18禁在线播放成人免费| 精品人妻熟女av久视频| 国产免费男女视频| 国产男人的电影天堂91| 麻豆av噜噜一区二区三区| 欧美激情久久久久久爽电影| 免费高清视频大片| 国产单亲对白刺激| 我的老师免费观看完整版| 露出奶头的视频| 欧美成人免费av一区二区三区| 在线a可以看的网站| 精品久久久久久久久av| 精品日产1卡2卡| av.在线天堂| 久久久久久久午夜电影| 日韩欧美在线二视频| 麻豆久久精品国产亚洲av| or卡值多少钱| 久久久久久久亚洲中文字幕| 日韩欧美精品v在线| 国产精品久久视频播放| 在线播放无遮挡| 国产精品久久久久久精品电影| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| 一区福利在线观看| 午夜激情欧美在线| 女同久久另类99精品国产91| 免费看美女性在线毛片视频| 国产大屁股一区二区在线视频| 91久久精品电影网| 久久亚洲真实| 亚洲成av人片在线播放无| 国产高清激情床上av| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 免费不卡的大黄色大毛片视频在线观看 | 一本一本综合久久| 日本熟妇午夜| 嫩草影院精品99| 亚洲不卡免费看| 999久久久精品免费观看国产| 美女xxoo啪啪120秒动态图| 日韩欧美在线二视频| 久久热精品热| 午夜免费激情av| av在线亚洲专区| 国产免费av片在线观看野外av| 成人av在线播放网站| 少妇熟女aⅴ在线视频| 在线免费观看的www视频| 国内精品久久久久久久电影| 熟女电影av网| av视频在线观看入口| 国内精品久久久久精免费| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 亚洲欧美日韩东京热| 亚洲综合色惰| 国产单亲对白刺激| 亚洲精品亚洲一区二区| videossex国产| 国产在线男女| 亚州av有码| 亚洲精品456在线播放app | 最新在线观看一区二区三区| 熟女电影av网| 久久久久久久精品吃奶| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产免费一级a男人的天堂| 一区二区三区四区激情视频 | 亚洲精品乱码久久久v下载方式| 12—13女人毛片做爰片一| 岛国在线免费视频观看| 欧美日本视频| 亚洲国产欧美人成| 中文字幕久久专区| 日本 av在线| 狂野欧美白嫩少妇大欣赏| h日本视频在线播放| or卡值多少钱| 亚洲av成人精品一区久久| 国产精品福利在线免费观看| 熟女人妻精品中文字幕| 免费av不卡在线播放| 黄色视频,在线免费观看| 级片在线观看| 日韩欧美国产一区二区入口| 麻豆成人午夜福利视频| 麻豆成人av在线观看| www.色视频.com| 亚洲熟妇熟女久久| avwww免费| 国产精品人妻久久久影院| 麻豆国产av国片精品| 国产大屁股一区二区在线视频| 动漫黄色视频在线观看| 亚洲熟妇熟女久久| 日本黄色片子视频| 国产激情偷乱视频一区二区| 少妇熟女aⅴ在线视频| 久久婷婷人人爽人人干人人爱| 91久久精品国产一区二区成人| 欧美丝袜亚洲另类 | 国产精品av视频在线免费观看| 中文字幕精品亚洲无线码一区| av福利片在线观看| а√天堂www在线а√下载| 国产伦精品一区二区三区四那| 久久婷婷人人爽人人干人人爱| 国产视频内射| 国产精品久久久久久亚洲av鲁大| 午夜爱爱视频在线播放| 一级黄色大片毛片| 欧美黑人巨大hd| 国产午夜精品论理片| 精品午夜福利视频在线观看一区| 久久久久久久午夜电影| 久9热在线精品视频| 又黄又爽又刺激的免费视频.| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 久久久久性生活片| 国产精品一区www在线观看 | 国产精品国产高清国产av| 亚洲美女黄片视频| 国产黄a三级三级三级人| 热99re8久久精品国产| 中文资源天堂在线| 欧美bdsm另类| 成人性生交大片免费视频hd| 日本黄大片高清| 在线观看午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 免费看美女性在线毛片视频| 在线播放无遮挡| 亚洲熟妇中文字幕五十中出| 久久久久久久亚洲中文字幕| 精品乱码久久久久久99久播| 乱码一卡2卡4卡精品| 国产主播在线观看一区二区| 老司机福利观看| 最好的美女福利视频网| 国产中年淑女户外野战色| 亚洲国产高清在线一区二区三| 女人十人毛片免费观看3o分钟| 国产亚洲精品综合一区在线观看| 成年女人永久免费观看视频| 精品久久久久久久久久免费视频| 99热网站在线观看| 亚洲不卡免费看| 色噜噜av男人的天堂激情| 成年免费大片在线观看| 亚洲精品一区av在线观看| 精品人妻熟女av久视频| 久久久久久九九精品二区国产| 国产亚洲精品久久久com| 欧美最新免费一区二区三区| 嫩草影院入口| 黄色日韩在线| 毛片一级片免费看久久久久 | 欧美精品国产亚洲| 直男gayav资源| 久久精品国产自在天天线| 国产69精品久久久久777片| 22中文网久久字幕| 在线免费十八禁| 亚洲人成网站在线播放欧美日韩| 最后的刺客免费高清国语| 伦理电影大哥的女人| 日本免费a在线| 91狼人影院| 国产一区二区激情短视频| 国产成人a区在线观看| 高清日韩中文字幕在线| 亚洲男人的天堂狠狠| 国产精品无大码| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费| 久久国内精品自在自线图片| 成年女人看的毛片在线观看| 国模一区二区三区四区视频| 亚洲专区国产一区二区| 日韩欧美精品免费久久| 国产视频一区二区在线看| 久久亚洲真实| 精品福利观看| 啪啪无遮挡十八禁网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲aⅴ乱码一区二区在线播放| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 91在线观看av| 国内毛片毛片毛片毛片毛片| 亚洲性久久影院| 极品教师在线免费播放| 亚洲,欧美,日韩| 国产精品女同一区二区软件 | www.www免费av| .国产精品久久| 日韩av在线大香蕉| 亚洲av二区三区四区| 成人午夜高清在线视频| 免费无遮挡裸体视频| 69人妻影院| 欧美+日韩+精品| 男女之事视频高清在线观看| 欧美激情国产日韩精品一区| 日本免费a在线| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放| 老熟妇乱子伦视频在线观看| 欧美成人a在线观看| 在线观看66精品国产| 毛片一级片免费看久久久久 | 精品免费久久久久久久清纯| 伦精品一区二区三区| 深爱激情五月婷婷| 国产探花极品一区二区| 欧美日韩乱码在线| 少妇熟女aⅴ在线视频| 一区二区三区四区激情视频 | 国产午夜精品论理片| 国产精品女同一区二区软件 | 男插女下体视频免费在线播放| 亚洲中文字幕日韩| 在现免费观看毛片| 精品人妻偷拍中文字幕| 97超级碰碰碰精品色视频在线观看| av中文乱码字幕在线| 麻豆av噜噜一区二区三区| 午夜福利18| 日本黄色片子视频| eeuss影院久久| 别揉我奶头 嗯啊视频| 日韩欧美一区二区三区在线观看| av在线老鸭窝| 少妇的逼好多水| 麻豆一二三区av精品| 琪琪午夜伦伦电影理论片6080| 久久久久久久久中文| 嫩草影院精品99| 伦精品一区二区三区| 久久久久久久精品吃奶| 欧美日韩综合久久久久久 | 精品福利观看| 欧美激情在线99| 日韩欧美在线乱码| 2021天堂中文幕一二区在线观| 日本五十路高清| 亚洲黑人精品在线| 国产伦精品一区二区三区四那| 久久久成人免费电影| 欧美色欧美亚洲另类二区| 国产精品98久久久久久宅男小说| 91久久精品国产一区二区三区| 国产又黄又爽又无遮挡在线| 午夜影院日韩av| 久9热在线精品视频| 色综合色国产| 天堂动漫精品| 成人亚洲精品av一区二区| 亚洲精品亚洲一区二区| 精品一区二区三区视频在线| 在线观看舔阴道视频| 在线免费观看的www视频| 久久草成人影院| 色在线成人网| 97超级碰碰碰精品色视频在线观看| 成人av一区二区三区在线看| 久久6这里有精品| 亚洲无线观看免费| 国产男人的电影天堂91| 欧美日韩精品成人综合77777| 久久天躁狠狠躁夜夜2o2o| 99热这里只有是精品50| 男人舔女人下体高潮全视频| 美女 人体艺术 gogo| 国产亚洲精品久久久久久毛片| 亚洲av熟女| 午夜精品一区二区三区免费看| ponron亚洲| 高清毛片免费观看视频网站| 国产探花极品一区二区| 亚洲无线在线观看| 久久久久久九九精品二区国产| 女生性感内裤真人,穿戴方法视频| 免费观看人在逋| 亚洲在线自拍视频| 国产一区二区三区av在线 | 免费在线观看影片大全网站| 性插视频无遮挡在线免费观看| 一夜夜www| 又黄又爽又免费观看的视频| 99久久精品国产国产毛片| 91午夜精品亚洲一区二区三区 | 国产成年人精品一区二区| 亚洲av中文字字幕乱码综合| 欧美zozozo另类| 全区人妻精品视频| 无人区码免费观看不卡| 最新中文字幕久久久久| 亚洲在线观看片| 国产精品精品国产色婷婷| 欧美不卡视频在线免费观看| 99久久精品国产国产毛片| 日日啪夜夜撸| 欧美国产日韩亚洲一区| 色哟哟·www| 亚洲,欧美,日韩| 成人无遮挡网站| 国产女主播在线喷水免费视频网站 | 中文字幕免费在线视频6| 亚洲精品色激情综合| 人人妻人人看人人澡| 91久久精品国产一区二区成人| 国产一区二区三区视频了| 亚洲人成伊人成综合网2020| 久久精品夜夜夜夜夜久久蜜豆| 国产在视频线在精品| 成人特级av手机在线观看| 亚洲午夜理论影院| 国产精品一区二区三区四区久久| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 99riav亚洲国产免费| 久久久久久九九精品二区国产| 99久国产av精品| 一区福利在线观看| 国产伦精品一区二区三区四那| 日韩强制内射视频| 变态另类成人亚洲欧美熟女| 国产亚洲精品久久久com| 欧美日本亚洲视频在线播放| 如何舔出高潮| 变态另类丝袜制服| 国产在线男女| 欧美成人a在线观看| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 天天一区二区日本电影三级| 欧美色视频一区免费| 日本成人三级电影网站| 真人一进一出gif抽搐免费| 国产真实伦视频高清在线观看 | 小说图片视频综合网站| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 99久国产av精品| 天堂动漫精品| 日韩欧美在线乱码| 国模一区二区三区四区视频| 我的老师免费观看完整版| 91麻豆av在线| 91麻豆精品激情在线观看国产| 国产亚洲精品综合一区在线观看| 国产黄色小视频在线观看| 国产一级毛片七仙女欲春2| 国产亚洲精品久久久com| 国产精品无大码| 九九爱精品视频在线观看| 欧美一区二区国产精品久久精品| 久久精品国产亚洲av涩爱 | 男女做爰动态图高潮gif福利片| 99精品在免费线老司机午夜| 国产成年人精品一区二区| 尾随美女入室| 偷拍熟女少妇极品色| av福利片在线观看| 1024手机看黄色片| 97人妻精品一区二区三区麻豆| 韩国av在线不卡| 高清日韩中文字幕在线| 欧美潮喷喷水| 日本 欧美在线| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲欧美98| 午夜影院日韩av| 成人美女网站在线观看视频| 国国产精品蜜臀av免费| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 一级a爱片免费观看的视频| 亚洲不卡免费看| 中文字幕熟女人妻在线| 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 欧美日韩综合久久久久久 | 欧美高清成人免费视频www| 天天躁日日操中文字幕| 亚洲中文字幕日韩| 国产v大片淫在线免费观看| 草草在线视频免费看| 国产亚洲精品av在线| 日韩欧美免费精品| 99久久久亚洲精品蜜臀av| 夜夜夜夜夜久久久久| 五月玫瑰六月丁香| 午夜免费男女啪啪视频观看 | 中国美女看黄片| 亚洲,欧美,日韩| av女优亚洲男人天堂| 又粗又爽又猛毛片免费看| 一本精品99久久精品77| 麻豆国产97在线/欧美| av天堂中文字幕网| videossex国产| 97碰自拍视频| 欧美日韩亚洲国产一区二区在线观看| 国产av不卡久久| 久久精品国产亚洲av涩爱 | 身体一侧抽搐| 国产美女午夜福利| 午夜福利欧美成人| 国产色爽女视频免费观看| 久久久久久久久大av| 少妇的逼水好多| 禁无遮挡网站| 嫩草影视91久久| 国产精品一区二区免费欧美| 亚洲狠狠婷婷综合久久图片| 欧美日韩黄片免| a在线观看视频网站| 五月伊人婷婷丁香| 午夜免费男女啪啪视频观看 | 国产 一区精品| 亚洲图色成人| 免费大片18禁| 久久久久免费精品人妻一区二区| 毛片女人毛片| 男女做爰动态图高潮gif福利片| 啦啦啦观看免费观看视频高清| 亚洲精华国产精华精| 国内精品美女久久久久久| 久久精品国产亚洲av天美| 国产黄a三级三级三级人| 国产视频一区二区在线看| 一本一本综合久久|