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

    Processes Leading to Second-Year Cooling of the 2010–12 La Ni?na Event, Diagnosed Using GODAS

    2015-02-24 03:39:54FENGLichengZHANGRongHuaWANGZhangguiandCHENXingrong
    Advances in Atmospheric Sciences 2015年3期

    FENG Licheng,ZHANG Rong-Hua,WANG Zhanggui,and CHEN Xingrong

    1National Marine Environmental Forecasting Center,State Oceanic Administration,Beijing100081

    2Key Laboratory of Ocean Circulation and Waves,Institute of Oceanology, Chinese Academy of Sciences,Qingdao266071

    3Earth System Science Interdisciplinary Center(ESSIC),University of Maryland, College Park,Maryland,USA,20740

    Processes Leading to Second-Year Cooling of the 2010–12 La Ni?na Event, Diagnosed Using GODAS

    FENG Licheng?1,ZHANG Rong-Hua2,3,WANG Zhanggui1,and CHEN Xingrong1

    1National Marine Environmental Forecasting Center,State Oceanic Administration,Beijing100081

    2Key Laboratory of Ocean Circulation and Waves,Institute of Oceanology, Chinese Academy of Sciences,Qingdao266071

    3Earth System Science Interdisciplinary Center(ESSIC),University of Maryland, College Park,Maryland,USA,20740

    Isopycnal analyses were performed on the Global Ocean Data Assimilation System(GODAS)to determine the oceanic processes leading to so-called second-year cooling of the La Ni?na event.In 2010–12,a horseshoe-like pattern was seen, connecting negative temperature anomalies off and on the Equator,with a dominant infuence from the South Pacifc.During the 2010 La Ni?na event,warm waters piled up at subsurface depths in the western tropical Pacifc.Beginning in early 2011, these warm subsurface anomalies propagated along the Equator toward the eastern basin,acting to reverse the sign of sea surface temperature(SST)anomalies(SSTAs)there and initiate a warm SSTA.However,throughout early 2011,pronounced negative anomalies persisted off the Equator in the subsurface depths of the South Pacifc.As isopycnal surfaces outcropped in the central equatorial Pacifc,negative anomalies from the subsurface spread upward along with mean circulation pathways,naturally initializing a cold SSTA.In the summer,a cold SSTA reappeared in the central basin,which subsequently strengthened due to the off-equatorial effects mostly in the South Pacifc.These SSTAs acted to initiate local coupled air–sea interactions,generating atmospheric–oceanic anomalies that developed and evolved with the second-year cooling in the fall of 2011.However,the cooling tendency in mid-2012 did not develop into another La Ni?na event,since the cold anomalies in the South Pacifc were not strong enough.An analysis of the 2007–09 La Ni?na event revealed similar processes to the 2010–12 La Ni?na event.

    La Ni?na,second-year cooling,off-equatorial effects,isopycnal analyses,circulation pathways,GODAS

    1. Introduction

    The El Ni?no–Southern Oscillation(ENSO)is the leading mode of interannual variability in the tropical Pacifc climate system,signifcantly impacting global weather and climate. In the past several decades,extensive studies have led to substantial progress in understanding,modeling and predicting El Ni?no events(e.g.,McCreary and Anderson,1984; Cane and Zebiak,1985;Zebiak and Cane,1987;Philander, 1992;Wang et al.,2011,2013).The delayed oscillator mechanism has been proposed to explain ENSO dynamics and its interannual oscillation within the tropical Pacifc climate system(Battisti and Hirst,1989).This theory emphasizes equatorial wave processes(Rossby wave and its refection along the low-latitude western boundaryinto a Kelvin wave).Another is the recharge/discharge mechanism(Jin,1997), which focuses on water exchange in the ocean on and off the Equator.As implied by these theories,the ENSO can be a cyclic oscillation between El Ni?no and La Ni?na conditions within the tropical Pacifc climate system.

    However,as observed,the ENSO also exhibits signifcant variability from one event instance to another.For example,multi-year cooling events can be seen during ENSO cycles from historical SST data(e.g.,Hu et al.,2014).During 2010–12,the tropical Pacifc had a persistent La Ni?na condition,with a second-year sea surface cooling that occurred in the fall of 2011.Further,many coupled models have failed to predict the Ni?no 3.4 sea surface temperature(SST)cooling when initialized from early-to mid-2011.Yet,one intermediate coupled model—an integrated climate model(ICM) operated at the Earth System Science Interdisciplinary Center(ESSIC),University of Maryland(UMD),the so-called ESSIC ICM(Zhang et al.,2003,2005)—gave a successful forecast of the 2011 negative SSTAs with a lead time of oneyear or so(Zhang et al.,2013,2014aZhang,R.-H,L.C.Feng,and Z.G.Wang,2014:Role of atmospheric wind forcing in the second-year cooling of the 2010–12 La Ni?na event.Atmos. Sci.Lett.,submitted.).This presents a challenge to the ENSO prediction community and indicates an urgent need to understand processes leading to the secondyear cooling.

    Previously,ICM-based experiments were carried out to examine the roles played by the temperature of subsurface water entrained into the mixed layer,and wind forcing (Zhang et al.,2013,2014aZhang,R.-H,L.C.Feng,and Z.G.Wang,2014:Role of atmospheric wind forcing in the second-year cooling of the 2010–12 La Ni?na event.Atmos. Sci.Lett.,submitted.).The reappearance of a negative SSTA in the central equatorial Pacifc in early summer of 2011 was closely related to off-equatorial thermal anomalies in the South Pacifc.However,the three-dimensional structure and evolution of these have not been illustrated,as theoceanicprocessesresponsibleforthesecond-yearcooling duringthe2010–12LaNi?na eventare still poorlyunderstood. The causes of the occurrence of a multi-year La Ni?na in general,and the 2011–12 La Ni?na event in particular,are not fully understood(Hu et al.,2014).

    In this paper,we examine the oceanic processes responsible for the second-year cooling of the 2010–12 La Ni?na event using reanalysis data,with a focus on the roles played by off-equatorial subsurface anomalies in the South Pacifc.To better represent pathways,isopycnal analyses were performed using three-dimensional temperature and salinity felds(Zhang and Rothstein,2000).Since subsurface temperature anomalies tend to propagate along density surfaces,an isopycnal analysis can better characterize the threedimensional structure and time evolution in a natural and physical way,therefore enabling us to trace pathways consistently throughout the basin.Our major fnding was that a distinct pathway of off-equatorial temperature anomalies occurred along the South Equatorial Current(SEC),clearly associated with the onset of second-year cooling during the 2010–12 La Ni?na event.Through examining the subsurface temperature evolution on isopycnals,connections were more clearly illustrated between thermal anomalies at the subsurface and surface,and off and on the Equator,leading to an improved understanding of ENSO variability.Additionally,re-evaluating the historical ENSO evolution showed that another multi-year cooling case occurred in the tropical Pacifc in 2007–09.The similarities and differences of these two events were analyzed to describe the nature of these strikingly different ENSO evolutions associated with various forcings and feedbacks within the Pacifc climate system.

    The remainder of the paper is organized as follows.We introduce the data and methodologyused in this work in section 2.The results are presented in section 3,followed by a summary and discussion in section 4.

    2. Data and methodology

    Monthly-mean data for currents,sea surface height,temperature and salinity came from the Global Ocean Data Assimilation System(GODAS)(Behringer and Xue,2004),operational at the National Centers for Environmental Prediction(NCEP).GODAS has a horizontal resolution of 1°× 1/3°in the zonal and meridional directions;it has 40 levels in the vertical,with a 10 m resolutionin the upper 200 m.We used the GODAS data coveringthe periodfromJanuary1980 throughDecember2012.Additionally,surface winds at 10 m height were from the NCEP–NCAR(National Center for Atmospheric Research)Reanalysis(Kalnay et al.,1996),with a longitudinal and latitudinal resolution of 1.904°×1.875°on a T62 Gaussian grid(192×94).

    Long-term climatological felds were formed from the period 1980–2012,including monthly-mean current vectors. Interannualanomaliesfortemperature,windstress andothers werethencalculatedrelativetotheir climatologicalfelds.Finally,isopycnal surfaces were estimated using monthly temperatureandsalinity data.Thetemperatureanomaliesat level depths were interpolated to constant density surfaces by using a cubic spline.Climatological current vectors on isopycnal surfaces were formed in the same way.In this study, interannual anomaly felds on isopycnal surfaces were used to investigate the roles played by anomalous temperature advection in the 2010–12 and 2007–09 La Ni?na events.

    3. Results

    3.1.SST evolution

    Figure 1 illustrates the horizontal distributions of SSTAs and surface wind anomalies for selected time intervals in 2011.In January,there was a La Ni?na state over the tropical Pacifc.Consequently,negative SSTAs prevailed in the central and eastern tropical Pacifc with the maxima exceeding-2°C between 150°and 170°W along the Equator.Surface easterly winds were stronger than normal over the western central equatorial Pacifc and southeasterly wind anomalies dominated off the Equator in the South Pacifc(Fig.1a). Thereafter,the cold SSTA diminished and the SSTA became normal in the eastern tropical Pacifc domain.Simultaneously,the wind stress anomalies weakened in most regions (Fig.1b).This warming process persisted during the following months and peaked in June,when a neutral SST state prevailed throughout the Equator except for a weak negative anomaly along 160°W.At this time,easterly wind anomalies weakened in the central tropical Pacifc(Fig.1d).In August, negative SSTA strengthened in the central equatorial Pacifc (Fig.1f),and this cooling tendency persisted during the following months(Figs.1g–h).

    Themechanismofformationofthe coldSSTA in the central eastern equatorial Pacifc during mid–late 2011 has not been fully explained.Some possible factors,such as wind forcing or a subsurface thermal anomaly,may play an important role.Note that southeasterly wind became stronger in the tropical South Pacifc(Fig.1e),forcing the cold waters located in the South Pacifc to move to the equatorial band (Fig.1e)and leadingto the negativeSSTA.However,the cur-rent driven by anomalous wind was not enough to produce strong and persistent negative SSTAs in the central equatorial Pacifc,especially after the wind anomalies changed direction during September and October(Figs.1g–h).Otherprocesses,such as subsurface effects,are required to fully understand the cause of the second-year cooling.

    3.2.Subsurface temperature anomaly pathway

    The climatological Bernoulli function(B)was calculated on the isopycnal surfaces to study the mean fow pattern.According to Cox and Bryan(1984),Bcan be written as

    whereσ=ρ-1000 is an isopycnal surface,andρis density with units of kg m-3;ρ0is mean density;gis the acceleration due to gravity,andηis dynamic height.Brepresents geostrophic streamlines that measure the geostrophic fow away from the Equator;thus,it can be used to illustrate fow paths on isopycnal surfaces.

    Figures 2a and 2c show the mean depth distributions of the 23.4 and 25.2 isopycnal surfaces.These two isopycnal surfaces had similar patterns in the tropical Pacifc.On the Equator,the thermocline was deep in the west and shallow in the east.The deepest regions on the isopycnal surfaces were located around 15°N and 5°S,respectively,in the western central Pacifc,with a relatively shallow band between 6°and 10°N.The isopycnal surfaces shoaled eastward along the Equator and reached minima in the far-eastern Pacifc.The 23.4 isopycnal surface intersected with the sea surface(i.e., outcropped)in the central and eastern basin on and south of the Equator(Fig.2a).

    Pathways along which off-equatorial waters move onto the Equator have been examined in many studies(e.g.Zhang et al.,1999;Zhang and Busalacchi,1999;Wang et al.,2007). However,mostpreviousanalysesfocusedontheeffectsinthe North Pacifc,with fewer studies in the South Pacifc.Chang et al.(2001)pointed out the potential importance of south tropical Pacifc variability in the decadal modulation of the ENSO.Luoetal.(2003)investigatedtheoriginofthedecadal ENSO-like variation.Luo et al.(2005)carried out 49-year simulations,and found that decadal variability of temperatureandsalinityalongtheEquatororiginatesfromsubsurface spiciness anomalies in the South Pacifc.

    FromFigs.2b andd,onecan see clear pathways originating from the southeastern tropical Pacifc:water carried by the South Equator Current(SEC)extending northwestward to south of the equatorial band and then transported by the strong Equator Undercurrent(EUC)onto the Equator.The South Pacifc water pathways intersect with the surface in the eastern equatorial and Southeast Pacifc domain(Fig.2b).

    Figure 3 gives subsurface temperature anomalies evaluated on the 25.2 isopycnal surface(see Fig.2c for its depth information)at some selected time periods in 2011;the vertical distribution of temperature anomalies in the upper ocean along the Equator is presented in Fig.4.During the 2010–11 La Ni?na event,there was a buildup of warm waters in the western Pacifc Ocean due to stronger than normal easterly winds in the central basin,characterizedby positive thermal anomalies in the upper ocean.For example,in January 2011,a large positive anomaly was observed in the western central tropical Pacifc and a negative anomaly was locatedin the central eastern tropical Pacifc regions.These two anomaly bands with opposite signs intersected along 160°W with a sharp temperature front(Figs.3a and 4a).Beginning in early 2011,accompanied by the seasonal strengthening of the EUC,warm waters in the western Pacifc expanded eastward across the Equator;cold anomalies in the central eastern equatorial Pacifc diminished and reversed to above normal(Figs.3b and 4b).This warming tendency peakedin April(Fig.3c),when positive temperature anomalies occupied almost the whole equatorial Pacifc except for near 150°W.Temperature anomalies reached more than 2°C in the far-eastern equatorial Pacifc.In the meantime,cold waters retreated to northeastern and southeastern regions off the Equator.As seen from the vertical section along the Equa-tor(Fig.4c),cold waters shrank back dramatically,and were confned to a narrow region of the central Pacifc.

    In May,positive anomalies along the Equator were seen to have two separate western and eastern bands(Fig.3d), with below-normal temperature anomalies amplifed in the regions of 140°–160°W(Fig.4d).Subsequently,the neg-ative anomalies dominated over the central Pacifc in June (Fig.3e),forming a horseshoe-like thermal anomaly pattern connecting large negative thermal anomalies on and off the Equator.Comparing Figs.3e and 3d,the EUC decelerated in the far-eastern equatorial Pacifc in June(Yu et al.,1997), but the off-equatorialcold anomalies strengthenedin the central South Pacifc.These changes were in favor of cold water advection to the equatorial regions through the well-defned SouthPacifcwaterpathway(Fig.2d),andthenextendedinto the equatorial region to combine with the negative anomalies located north of the Equator.In July,the EUC weakened further,and was even replaced by the SEC in the eastern Pacifc on the 25.2 isopycnal surface.At this time,cold anomalies were transported by SEC from the Southeastern Pacifc, and amplifed on and off the central equatorial Pacifc.This cooling tendency persisted in the following months.Positive anomalies along the Equator disappeared gradually,and cold anomalies dominated over the whole equatorial band(Figs. 3g and h).The vertical sections along the Equator displayed the same behavior(Figs.4e–h).

    3.3.Phase relationships between subsurface and surface temperature anomalies

    As analyzed above,the subsurface thermal anomalies at the Equator exhibited similar evolution to the SSTAs,but with a 2 month phase lead time:negative sea temperatureanomalies(Fig.3)re-strengthened at subsurface depths in June,while those in Fig.1 re-strengthened at the sea surface in August.This indicates the existence of close links between subsurface temperature anomalies and the SSTAs. During boreal spring,positive SSTAs in the far-eastern equatorial Pacifc(Figs.1b and c)cannot be explained by surface temperature advection,and they are likely to originate from the outcrop of subsurface warm anomalies(Figs.4b and c). This process can be described as follows.During the previous La Ni?na event,warm waters piled up in the western Pacifc Ocean due to stronger than normal easterly winds in the central basin.As the EUC became seasonally strengthened, the subsurface warm water was transported from the western Pacifc to the central and eastern Pacifc across the Equator (Figs.3b–c).Since the thermocline shoaled eastward(Figs. 2a and c),the warm water was exposed to the sea surface in the eastern Pacifc,acting to generate positive SSTAs(Figs. 1b and c;Figs.4b and c).

    As for the sea surface coolingin the fall of 2011,it can be traced to the subsurface anomalies.Beginning in mid-2011, subsurfacecoldanomalieslocatedinthesoutheasterntropical PacifcwerecontinuallyadvectednorthwestwardbytheSEC, to the south of the equatorial band,and then transported by the EUC to the Equator,where they were accumulated(Figs. 3e–h).But how did the subsurface cold water in the central Pacifc affect the sea surface?Since there was no systematic surface wind stress curl(fgures not shown),the related Ekman pumpingwas not a major factor infuencingthe outcropping of subsurface cold water,so the upwelling can only be driven by oceanic processes.Figure 5 presents the tempera-ture anomalies,and the horizontal and vertical velocity felds on the 23.4 and 25.2 isopycnals.The convergence pattern of the horizontal currents agreed reasonably well with the vertical velocity feld.For example,the convergence center was located on the Equator near 110°W,where the EUC met the SEC,giving rise to a strong upwelling(Fig.5b).

    InJune,small coldanomalieswere accompaniedbyweak upwelling in the central equatorial Pacifc(Fig.5a).With time,both cold anomalies and vertical velocity strengthened in the central equatorial Pacifc on the 25.2 isopycnal surface (Figs.5b and c).For example,in June the cold anomalies were confned between 130°W and 150°W along the Equator,but it dominated the eastern central Pacifc in July.These changes were induced by the weakened EUC and strengthened SEC,which favored the accumulation of cold water at the Equator.Figures 5e–h indicate that the vertical current in the upperlayerwas strongerthanthat at the lowerlayer(Figs. 5a–d),and the cold anomalies appeared later than that on the subsurface layer,which confrmed that the cold water originated from the subsurface.As discussed above,there was a clear pathway along which subsurface cold water was transported to the sea surface.Firstly,the subsurface cold water located in the southeastern tropical Pacifc was advected by the SEC south of the Equator.Subsequently,the EUC transported it to the equatorial Pacifc,where the EUC met the SEC and induced upwelling.Finally,under the effects of the EUC and SEC,the cold water spread upward and westward to the sea surface.

    In the fall of 2011(Figs.1g and h),negative SSTAs dominated in the central and eastern equatorial Pacifc basin.The negative SSTAs in the east affected winds to the west,which in turn affected the thermocline and SST in the east.This essentiallyinvolvedinteractionsamonganomaliesofSST,wind and the thermocline,forming a coupling loop and leading to the second-year cooling during 2010–12.

    3.4.Evolution during the 2012 decay phase

    Figure 6 gives the horizontal distributions of SSTAs and surface wind anomalies at some selected time intervals in 2012.From February onwards,the cold SSTA diminished and the SSTA became normal in the eastern tropical Pacifc domain(Fig.6a).This warming process persisted during the following months,and the SSTAs in the central and eastern tropical Pacifc rose above normal(Fig.6d),except in the far-eastern Pacifc.Figure 7 illustrates the subsurface temperature anomalies evaluated on the 25.2 isopycnal surface at some selected time periods in 2012;the vertical distribution of temperature anomalies in the upper ocean along the Equator is presented in Fig.8.Beginning in early 2012,accompanied by the seasonal strengthening of the EUC,warm waters in the western Pacifc expanded eastward across the Equator (Figs.7a and 8a).In May,with the seasonal maximum EUC, warm anomalies occupied the whole central eastern equatorial Pacifc(Figs.7band8b).Negativeanomaliesre-emerged twice(Figs.7c and e;Figs.8c and e)in the central equatorial Pacifc,since the EUC decelerated from June onwards. However,these coolingprocessesdidnotpersist anddevelop, perhaps because the cold anomalies in the South Pacifc weretoo weak to provideenoughcold water(Figs.7d–f;Figs.8d–f).Finally,the SSTAs did not return to the La Ni?na state,as happened during 2011.

    3.5.Evolution during the 2008 La Ni?na event

    Figure 9 gives the horizontal distributions of the SSTAs and surface wind anomalies at some selected time intervalsin 2008.In January,a La Ni?na state occupied the tropical Pacifc:negative SSTAs prevailed in the central and eastern tropical Pacifc with the maxima exceeding-2.5°C,located at 170°W along the Equator(Fig.9a).Thereafter,the cold SSTA diminished and the SSTA increased above normal in the far-eastern tropical Pacifc domain(Fig.9b).This warm-ingprocesspersistedduringthefollowingmonthsandpeaked in August(Fig.9d).In September,the negative SSTA restrengthened in the central equatorial Pacifc(Fig.9e),and this cooling tendency persisted during the following months (Fig.9f).

    Figure 10 illustrates the subsurface temperature anomalies evaluated on the 25.4 isopycnal surface at some selected time periods in 2008;the vertical distribution of temperature anomalies in the upper ocean along the Equator is presented in Fig.11.Beginning in early 2008,accompanied by the seasonal strengthening of the EUC,warm waters in the western Pacifc expanded eastward across the Equator(Fig.10b). This warming tendency peaked in mid-2008(Figs.10d and 11d),when positive temperature anomalies occupied almost the whole equatorial Pacifc.There was a 1–2 month lead time into the SSTAs.Compared with the warming process in 2011(Figs.3c and 4c),it lagged by about 2 months,possibly attributable to stronger negative anomalies in the eastern tropical Pacifc.Beginning in August,the subsurface cold anomalies located in the southeastern tropical Pacifc were continually advected northwestward by the SEC to the south of the equatorial band,and then transported by the EUC to the Equator,where they accumulated(Figs.10e–h).The cold anomalies were then transported by a vertical current to the sea surface and induced negative SSTAs.From September,cold water re-strengthened in the central-equatorial Pacifc(Fig.9e),and this cooling tendency persisted and extendedeastwardduringthefollowingmonths(Figs.9h);consequently,the double-trough La Ni?na developed.

    4. Summary and discussion

    The reanalysis products from GODAS were used to produce isopycnalsurfaces to better illustrate and understandthe processes leading to the second-year cooling of the 2010–12 La Ni?na event.We found anomaly patterns originating at depth from the southeastern tropical Pacifc that could be responsible for generating and sustaining negative SSTAs in the central equatorial Pacifc.

    A sequence of events leading to the La Ni?na conditions in the fall of 2011 was described.During the 2010 La Ni?na event,warm waters piled up at subsurface depths in the western tropical Pacifc.Beginning in early 2011,and accompanied by a strongEUC,subsurfacewarm waters in the western Pacifc transmitted eastward along the Equator.Positive temperature anomalies occupied the equatorial Pacifc in April, and cold waters retreated to northeastern and southeastern off-equatorial Pacifc regions.Since the thermocline shoaled along the Equator and was close to the surface in the eastern Pacifc,subsurface warm waters were directly exposed to the sea surface in the eastern Pacifc,and induced a warm SSTA. Normal SST conditions appeared in the central and eastern equatorial Pacifc in mid-2011.

    In August a negative SSTA reappeared in the central Pacifc.We hypothesized that this anomaly came from the subsurface cold waters off the Equator through the Southern Pacifc pathway.Based on the GODAS analyses,the processes were described as follows:Cold anomalies located in southeastern tropical Pacifc region were advected continually by the SEC northwestward to the south of the equatorial band, and then by the EUC northeastward to the Equator.With time,the EUC weakened and the SEC strengthened in the eastern equatorial Pacifc,inducing cold waters that accumulated in the central tropical Pacifc and then tended to spread upwardwith theconvergenceofhorizontalcurrentsandeventually outcropped to the surface.These subsurface-induced SSTAs actedtoinitiatelocalcoupledair–seainteractionsgenerating atmospheric–oceanic anomalies that developed and evolved with the second-year cooling in the fall of 2011.

    Further study of the 2012 processes indicated that the cooling tendency did not develop into another La Ni?na event, since the cold anomalies in the South Pacifc were not strong enough.An analysis around the 2007–09 La Ni?na event revealed similar evolution processes with around a 2—month phase lag,compared to the 2010–12 La Ni?na event.

    These analyses provide an observational basis for an understanding of the processes involved.The results can be used to explain the ways in which coupled models predict the second-year cooling case,and offer guidance for historical analysesforothermulti-yearcoolingevents.Furthersupporting modeling studies are needed to quantify the role played by off-equatorial subsurface anomalies in triggering La Ni?na events in the tropical Pacifc.Here,we discussed the effect of interannual variability on the multi-year cooling.The effect of modulation of decadal to interdecadal timescale variability on the multi-year cooling,such as tropical Pacifc decadal variability(Choi et al.,2013),requires further study.

    Acknowledgements.This work has benefted a great deal from Prof.A.J.BUSALACCHI’s support.This research was jointly supported by National Natural Science Foundation of China (Grant No.40906014),the Ocean Public Welfare Scientifc Research Project(Grant No.201205018-2),the National Key Basic Research Program of China(Grant No.2010CB950302),and the China Scholarship Council(CSC).ZHANG is supported partly by the National Science Foundation(NSF)(Grant No.ATM-0727668),NASA(Grant No.NNX08AI74G),and the National Oceanic and Atmospheric Administration(NOAA)(Grant No. NA08OAR4310885).

    REFERENCES

    Battisti,D.S.,and A.C.Hirst,1989:Interannual variability in the tropical atmosphere-ocean system:Infuence of the basis state,ocean geometry and nonlinearity.J.Atmos.Sci.,46, 1687–1712.

    Behringer,D.W.,and Y.Xue,2004:Evaluation of the global ocean data assimilation system at NCEP:The Pacifc Ocean. Preprints,Eighth Symp.on Integrated Observing and Assimilation Systems for Atmosphere,Oceans,and Land Surface, Seattle,WA,Amer.Meteor.Soc.,2.3.[Available online athttps://ams.confex.com/ams/84Annual/techprogram/paper 70720.htm.]

    Cane,M.A.,and S.E.Zebiak,1985:A theory for El Ni?no and the Southern Oscillation.Science,228,1085–1087.

    Chang,P.,B.S.Giese,L.Ji,H.F.Seidel,and F.Wang,2001: Decadal change in the South Tropical Pacifc in a global assimilation analysis.Geophys.Res.Lett.,28,3461–3464.

    Choi,J.,S.I.An.,S.W.Yeh,and J.Y.Yu,2013:ENSO-like and ENSO-induced tropical Pacifc decadal variability in CGCMs.J.Climate,26,1485–1501.

    Cox,M.D.,and K.Bryan,1984:A numerical model of the ventilated thermocline,J.Phys.Oceanogr.,14,674–687.

    Hu,Z.-Z.,A.Kumar,Y.Xue,and B.Jha,2014:Why were some La Ni?na followed by another La Ni?na?Climate Dyn.,42, 1029–1042,doi:10.1007/s00382-013-1917-3.

    Jin,F.-F.,1997:An equatorial ocean recharge paradigm for ENSO.Part I:Conceptual model.J.Atmos.Sci.,54,811–829.

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

    Luo,J.-J.,S.Masson,S.K.Behera,P.Delecluse,S.Gualdi, A.Navarra,and T.Yamagata,2003:South Pacifc Origin of the decadal ENSO-like variation as simulated by a coupled GCM.Geophys.Res.Lett.,30(24),2250,doi:10.1029/ 2003GL018649.

    Luo,Y.-Y.,L.M.Rothstein,R.-H.Zhang,and A.J.Busalacchi,2005:On the connection between South Pacifc subtropical spiciness anomalies and decadal equatorial variability in an ocean general circulation model.J.Geophys.Res.,110, C10002,doi:10.1029/2004JC002655.

    McCreary,J.P.Jr.,and D.L.T.Anderson,1984:A simple model of El Ni?no and the Southern Oscillation.Mon.Wea.Rev.,112, 934–946.

    Philander,S.G.H,1992:Ocean-atmosphere interactions in the tropics:A review of recent theories and models.J.Appl.Meteoro.,31,938–945.

    Wang,X.,C.Y.Li,and W.Zhou,2007:Interdecadal mode and its propagating characteristics of SSTA in the South Pacifc.Meteor.Atmos.Phys.,98,115–124,doi:10.1007/s00703-006-0235-2.

    Wang,X.,C.Z.Wang,W.Zhou,D.X.Wang,and J.Song,2011: Teleconnected infuence of NorthAtlanticsea surface temperature on the El Ni?no onset.Climate Dyn.,37,663–676,doi: 10.1007/s00382-010-0833-z.

    Wang,X.,C.Z.Wang,W.Zhou,L.Liu,and D.X.Wang,2013: Remote infuence of North Atlantic SST on the equatorial westerly wind anomalies in the western Pacifc for initiating an El Ni?no event:An Atmospheric General Circulation Model Study.Atmos.Sci.Lett.,14,107–111.

    Yu,Z.-J.,P.S.Schopf,and J.P.McCreary Jr.,1997:On the annual cycle of upper-ocean circulation in the eastern equatorial Pacifc.J.Phys.Oceanogr.,27 309–324.

    Zebiak,S.E.,and M.A.Cane,1987:A model El Ni?no-Southern Oscillation.Mon.Wea.Rev.,115,2262–2278.

    Zhang,R.-H.,and A.J.Busalacchi,1999:A possible link between off-equatorial warm anomalies propagating along the NECC path and the onset of the 1997–98 El Ni?no.Geophys.Res. Lett.,26(18),2873–2876.

    Zhang,R.-H.,and L.M.Rothstein,2000:Role of off-equatorial subsurface anomalies in initiating the 1991-1992 El Ni?no as revealed by the National Centers for Environmental Prediction ocean reanalysis data.J.Geophys.Res.,105(C3),6327–6339.

    Zhang,R.-H.,L.M.Rothstein,A.J.Busalacchi,and X.Z.Liang, 1999:The onset of the 1991-92 El Ni?no event in the tropical Pacifc Ocean:The NECC subsurface pathway.Geophys. Res.Lett.,26(7),847–850.

    Zhang,R.-H.,S.E.Zebiak,R.Kleeman,and N.Keenlyside, 2003:A new intermediate coupled model for El Ni?no simulation and prediction.Geophys.Res.Lett.,30(19),doi: 10.1029/2003GL018010.

    Zhang,R.-H.,S.E.Zebiak,R.Kleeman,and N.Keenlyside,2005: Retrospective El Nino forecast using an improved intermediate coupled model.Mon.Wea.Rev.,133,2777–2802.

    Zhang,R.-H.,F.Zheng,J.Zhu,and Z.G.Wang,2013:A successful real–time forecast of the 2010–11 La Ni?na event.Sci. Rep.,3,1108,doi:10.1038/srep01108.

    :Feng,L.C.,R.-H.Zhang,Z.G.Wang,and X.R.Chen,2015:Processes leading to the second-year cooling of the 2010–12 La Ni?na event,diagnosed using GODAS.Adv.Atmos.Sci.,32(3),424–438,

    10.1007/s00376-014-4012-8.

    (Received 7 February 2014;revised 12 July 2014;accepted 22 July 2014)

    ?Corresponding author:FENG Licheng

    Email:fenglich@nmefc.gov.cn

    午夜福利在线免费观看网站| 777米奇影视久久| 又大又爽又粗| 成在线人永久免费视频| 熟女av电影| 高清视频免费观看一区二区| 国产精品久久久久久精品电影小说| 免费av中文字幕在线| 亚洲欧洲国产日韩| 首页视频小说图片口味搜索 | 女性生殖器流出的白浆| 欧美老熟妇乱子伦牲交| 久久久精品免费免费高清| 午夜福利一区二区在线看| 两性夫妻黄色片| 18禁国产床啪视频网站| 一区福利在线观看| 乱人伦中国视频| 欧美日韩黄片免| 亚洲av综合色区一区| 国产精品一二三区在线看| 日日夜夜操网爽| 国产1区2区3区精品| 9色porny在线观看| 亚洲精品国产色婷婷电影| 一边亲一边摸免费视频| 18禁裸乳无遮挡动漫免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷色麻豆天堂久久| 麻豆av在线久日| 国产精品 国内视频| 久久女婷五月综合色啪小说| 亚洲精品一二三| 欧美精品高潮呻吟av久久| 国产日韩欧美视频二区| 又黄又粗又硬又大视频| 日日摸夜夜添夜夜爱| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品自拍成人| 大话2 男鬼变身卡| 国产精品国产三级专区第一集| 青草久久国产| 99精国产麻豆久久婷婷| 国产真人三级小视频在线观看| 日韩制服丝袜自拍偷拍| 亚洲 国产 在线| 悠悠久久av| av在线app专区| 亚洲国产精品国产精品| 嫩草影视91久久| 一个人免费看片子| 久热这里只有精品99| 精品人妻1区二区| 亚洲精品在线美女| 中文字幕制服av| 777久久人妻少妇嫩草av网站| 欧美中文综合在线视频| 高清av免费在线| 亚洲国产欧美在线一区| av一本久久久久| 成人手机av| 麻豆av在线久日| 老司机午夜十八禁免费视频| 激情五月婷婷亚洲| 成年人午夜在线观看视频| 丝袜在线中文字幕| 日韩 欧美 亚洲 中文字幕| 久久鲁丝午夜福利片| 国产一卡二卡三卡精品| 老司机影院毛片| 中文字幕人妻熟女乱码| 久久ye,这里只有精品| 国语对白做爰xxxⅹ性视频网站| 老司机靠b影院| 亚洲av成人不卡在线观看播放网 | 一级a爱视频在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 久久狼人影院| 精品久久久久久久毛片微露脸 | 狠狠婷婷综合久久久久久88av| 亚洲一区二区三区欧美精品| 成人国产av品久久久| 国产欧美亚洲国产| 久久av网站| 亚洲精品自拍成人| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品一区三区| 久久久久久久久免费视频了| 欧美日韩视频高清一区二区三区二| 国产成人a∨麻豆精品| 久久精品国产综合久久久| 日韩中文字幕欧美一区二区 | 免费高清在线观看日韩| 国产一区二区激情短视频 | 日本五十路高清| 国产亚洲av片在线观看秒播厂| 久久久国产欧美日韩av| 久久av网站| 99精国产麻豆久久婷婷| 人体艺术视频欧美日本| 欧美精品人与动牲交sv欧美| 最黄视频免费看| 一本—道久久a久久精品蜜桃钙片| 激情五月婷婷亚洲| 精品卡一卡二卡四卡免费| 亚洲欧美中文字幕日韩二区| 777久久人妻少妇嫩草av网站| 欧美97在线视频| 超碰97精品在线观看| 国产高清视频在线播放一区 | 天堂俺去俺来也www色官网| 亚洲av综合色区一区| 韩国精品一区二区三区| 欧美老熟妇乱子伦牲交| 肉色欧美久久久久久久蜜桃| 香蕉国产在线看| 亚洲国产欧美在线一区| 色播在线永久视频| 亚洲中文字幕日韩| 国产xxxxx性猛交| 国产欧美亚洲国产| 久久精品亚洲av国产电影网| 一级毛片 在线播放| 亚洲色图综合在线观看| 水蜜桃什么品种好| 涩涩av久久男人的天堂| 国产伦理片在线播放av一区| 精品福利永久在线观看| 一级毛片 在线播放| 免费一级毛片在线播放高清视频 | 日韩精品免费视频一区二区三区| 精品国产一区二区三区久久久樱花| 一区在线观看完整版| 成人18禁高潮啪啪吃奶动态图| 一级a爱视频在线免费观看| 日韩欧美一区视频在线观看| 亚洲国产日韩一区二区| 日韩av在线免费看完整版不卡| 免费不卡黄色视频| 国产激情久久老熟女| 久久久久久久久免费视频了| 婷婷色综合www| 亚洲精品日本国产第一区| 无限看片的www在线观看| 亚洲中文字幕日韩| 精品人妻一区二区三区麻豆| 亚洲精品久久成人aⅴ小说| 久久精品国产亚洲av涩爱| 国产成人精品在线电影| 天堂中文最新版在线下载| 人人妻人人澡人人看| 丰满迷人的少妇在线观看| 欧美日韩av久久| 亚洲国产精品国产精品| 日本黄色日本黄色录像| 日韩人妻精品一区2区三区| 亚洲免费av在线视频| 久久精品久久精品一区二区三区| 国产精品久久久久成人av| 亚洲欧美日韩另类电影网站| 精品亚洲成国产av| 一二三四社区在线视频社区8| 久久久亚洲精品成人影院| 亚洲成国产人片在线观看| 青青草视频在线视频观看| a级毛片黄视频| 无遮挡黄片免费观看| 五月开心婷婷网| 欧美日韩成人在线一区二区| 亚洲国产欧美在线一区| 亚洲精品av麻豆狂野| 亚洲中文字幕日韩| 欧美97在线视频| 岛国毛片在线播放| 国产成人影院久久av| 精品久久蜜臀av无| 9191精品国产免费久久| 欧美亚洲日本最大视频资源| avwww免费| 精品国产乱码久久久久久男人| 久9热在线精品视频| 国产一区二区 视频在线| 日韩 亚洲 欧美在线| 在线观看一区二区三区激情| 亚洲成色77777| 观看av在线不卡| 伦理电影免费视频| 久久国产精品影院| 搡老岳熟女国产| 免费观看人在逋| 国产欧美日韩综合在线一区二区| 亚洲av电影在线观看一区二区三区| 亚洲精品日本国产第一区| 久久久久久久大尺度免费视频| 日韩欧美一区视频在线观看| 99国产精品一区二区蜜桃av | 久久99精品国语久久久| 中文字幕人妻熟女乱码| 精品一区二区三区四区五区乱码 | 亚洲午夜精品一区,二区,三区| 三上悠亚av全集在线观看| 国产免费又黄又爽又色| 美女高潮到喷水免费观看| 母亲3免费完整高清在线观看| 人人妻人人澡人人看| 91老司机精品| 超碰成人久久| 亚洲熟女毛片儿| 午夜视频精品福利| 久久久欧美国产精品| 天天添夜夜摸| 波多野结衣一区麻豆| av天堂久久9| 国产精品99久久99久久久不卡| 国产一区有黄有色的免费视频| 国产在线视频一区二区| 亚洲一区二区三区欧美精品| 欧美亚洲 丝袜 人妻 在线| av一本久久久久| 国产视频一区二区在线看| 免费少妇av软件| 99精国产麻豆久久婷婷| 午夜老司机福利片| 亚洲成人手机| av线在线观看网站| 成人亚洲精品一区在线观看| 日韩制服骚丝袜av| 欧美日韩国产mv在线观看视频| 午夜福利视频在线观看免费| 免费观看人在逋| 亚洲国产欧美在线一区| 9191精品国产免费久久| 乱人伦中国视频| 九草在线视频观看| 18在线观看网站| 中文精品一卡2卡3卡4更新| 国产精品欧美亚洲77777| 丰满迷人的少妇在线观看| 日韩 亚洲 欧美在线| 亚洲精品一区蜜桃| 久久久欧美国产精品| 丝袜美腿诱惑在线| 成人黄色视频免费在线看| 美女扒开内裤让男人捅视频| 精品少妇久久久久久888优播| 五月天丁香电影| 乱人伦中国视频| 九草在线视频观看| 国产精品偷伦视频观看了| 国产一区二区在线观看av| 免费在线观看黄色视频的| 午夜激情久久久久久久| 丰满饥渴人妻一区二区三| 人人妻人人澡人人爽人人夜夜| xxx大片免费视频| 五月天丁香电影| 在线观看免费视频网站a站| 欧美在线一区亚洲| 夫妻午夜视频| 欧美成人午夜精品| 人人妻人人爽人人添夜夜欢视频| 国产精品久久久久久人妻精品电影 | 91九色精品人成在线观看| 波多野结衣av一区二区av| 51午夜福利影视在线观看| 亚洲精品国产一区二区精华液| 精品福利观看| 18禁国产床啪视频网站| 欧美日韩视频精品一区| 啦啦啦中文免费视频观看日本| 自线自在国产av| 少妇裸体淫交视频免费看高清 | 午夜影院在线不卡| 国产欧美亚洲国产| 久久性视频一级片| 一级毛片电影观看| 超碰成人久久| 国产在线观看jvid| 中文字幕人妻丝袜制服| 爱豆传媒免费全集在线观看| 五月天丁香电影| av视频免费观看在线观看| 婷婷色麻豆天堂久久| 国产成人av激情在线播放| 91字幕亚洲| 国产有黄有色有爽视频| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 国产欧美日韩一区二区三 | 视频区欧美日本亚洲| 久久久久视频综合| 2021少妇久久久久久久久久久| av在线老鸭窝| 亚洲男人天堂网一区| 国产在线观看jvid| 亚洲国产精品一区二区三区在线| 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 一区福利在线观看| 免费高清在线观看视频在线观看| 青春草亚洲视频在线观看| 亚洲专区国产一区二区| 18禁国产床啪视频网站| 精品福利永久在线观看| 精品第一国产精品| 美女视频免费永久观看网站| a级片在线免费高清观看视频| 久久免费观看电影| 少妇人妻 视频| 久久久久精品国产欧美久久久 | 男人操女人黄网站| 大码成人一级视频| 中文乱码字字幕精品一区二区三区| 久久毛片免费看一区二区三区| 又大又黄又爽视频免费| 老司机影院成人| 两个人免费观看高清视频| 欧美在线一区亚洲| 制服人妻中文乱码| 侵犯人妻中文字幕一二三四区| 精品视频人人做人人爽| 久久天躁狠狠躁夜夜2o2o | 亚洲少妇的诱惑av| 女性生殖器流出的白浆| 免费在线观看黄色视频的| 99久久99久久久精品蜜桃| 黑人巨大精品欧美一区二区蜜桃| 亚洲成色77777| 一本—道久久a久久精品蜜桃钙片| 高清欧美精品videossex| 久热这里只有精品99| 男女边摸边吃奶| 日韩av免费高清视频| 国产精品国产三级国产专区5o| h视频一区二区三区| 丝袜人妻中文字幕| 大片免费播放器 马上看| 另类亚洲欧美激情| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 18禁观看日本| 精品福利永久在线观看| 国产片内射在线| 黄色 视频免费看| 狠狠婷婷综合久久久久久88av| 亚洲欧洲日产国产| 狂野欧美激情性bbbbbb| 久久久久精品人妻al黑| 亚洲成国产人片在线观看| 建设人人有责人人尽责人人享有的| 国产精品.久久久| 成人三级做爰电影| 欧美亚洲日本最大视频资源| 肉色欧美久久久久久久蜜桃| 精品第一国产精品| 精品亚洲成国产av| 国产黄频视频在线观看| 一本久久精品| 国产91精品成人一区二区三区 | 赤兔流量卡办理| 操出白浆在线播放| 少妇 在线观看| 国产有黄有色有爽视频| 国产男女内射视频| 日韩精品免费视频一区二区三区| 男女高潮啪啪啪动态图| 午夜福利,免费看| 黄片播放在线免费| 国产精品九九99| 在线观看免费高清a一片| 国产在线视频一区二区| 亚洲中文日韩欧美视频| 成在线人永久免费视频| 日本欧美国产在线视频| a级毛片在线看网站| 老司机在亚洲福利影院| 十八禁人妻一区二区| 99热全是精品| 成人黄色视频免费在线看| 悠悠久久av| 操出白浆在线播放| 国产日韩欧美亚洲二区| 成年女人毛片免费观看观看9 | 少妇人妻久久综合中文| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9 | 国产精品亚洲av一区麻豆| 波多野结衣av一区二区av| 精品久久蜜臀av无| 国产av一区二区精品久久| 亚洲欧洲精品一区二区精品久久久| 国产一卡二卡三卡精品| 亚洲精品久久午夜乱码| 亚洲中文av在线| 午夜日韩欧美国产| 美女脱内裤让男人舔精品视频| 中文字幕精品免费在线观看视频| 男的添女的下面高潮视频| 少妇 在线观看| 嫁个100分男人电影在线观看 | 欧美黄色淫秽网站| 一区福利在线观看| 丰满迷人的少妇在线观看| 在线精品无人区一区二区三| 国产成人a∨麻豆精品| 大陆偷拍与自拍| 免费在线观看视频国产中文字幕亚洲 | 人成视频在线观看免费观看| 国产亚洲av高清不卡| 成人手机av| 久久久久久久大尺度免费视频| 亚洲欧洲国产日韩| 亚洲欧美色中文字幕在线| 国产男女内射视频| 亚洲男人天堂网一区| 久久精品久久精品一区二区三区| 国产男女超爽视频在线观看| av电影中文网址| 成人免费观看视频高清| 黑人欧美特级aaaaaa片| 巨乳人妻的诱惑在线观看| 观看av在线不卡| 国产精品免费视频内射| 视频区欧美日本亚洲| √禁漫天堂资源中文www| 无遮挡黄片免费观看| 久久精品久久久久久噜噜老黄| 欧美精品亚洲一区二区| 国产av精品麻豆| 亚洲欧美日韩高清在线视频 | 免费一级毛片在线播放高清视频 | 男女边吃奶边做爰视频| videos熟女内射| 国产欧美日韩综合在线一区二区| 欧美成人午夜精品| 99国产精品一区二区蜜桃av | 国产一区有黄有色的免费视频| 99热国产这里只有精品6| 久久精品久久久久久久性| 男女边摸边吃奶| 国产老妇伦熟女老妇高清| 一级,二级,三级黄色视频| 亚洲精品在线美女| 老司机午夜十八禁免费视频| 日韩精品免费视频一区二区三区| 黄片播放在线免费| 成在线人永久免费视频| 午夜激情av网站| 超色免费av| 天天躁夜夜躁狠狠躁躁| 91麻豆av在线| 国产又色又爽无遮挡免| 久久国产亚洲av麻豆专区| 中文字幕另类日韩欧美亚洲嫩草| 97精品久久久久久久久久精品| 亚洲熟女精品中文字幕| 国产日韩欧美亚洲二区| 黑人猛操日本美女一级片| www.999成人在线观看| 热re99久久国产66热| 亚洲七黄色美女视频| 亚洲九九香蕉| 日本av免费视频播放| 狂野欧美激情性bbbbbb| 精品免费久久久久久久清纯 | 乱人伦中国视频| 又大又黄又爽视频免费| 老司机靠b影院| 欧美精品一区二区免费开放| 少妇 在线观看| 母亲3免费完整高清在线观看| 女性被躁到高潮视频| 免费观看av网站的网址| 丰满少妇做爰视频| 天天添夜夜摸| 亚洲成人手机| 国产一区二区激情短视频 | 精品一区二区三卡| 成年女人毛片免费观看观看9 | 国产亚洲欧美精品永久| 亚洲成av片中文字幕在线观看| 99久久99久久久精品蜜桃| 国产女主播在线喷水免费视频网站| a级毛片黄视频| 如日韩欧美国产精品一区二区三区| 日韩制服丝袜自拍偷拍| 午夜久久久在线观看| 老汉色av国产亚洲站长工具| 欧美 亚洲 国产 日韩一| 免费看十八禁软件| 亚洲精品国产av蜜桃| 中文字幕色久视频| 亚洲男人天堂网一区| 91九色精品人成在线观看| 久久免费观看电影| 成年av动漫网址| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久久久久婷婷小说| 91字幕亚洲| 各种免费的搞黄视频| 18禁国产床啪视频网站| 精品国产国语对白av| a级毛片在线看网站| 中文字幕制服av| 啦啦啦在线观看免费高清www| 高清视频免费观看一区二区| 亚洲精品国产av成人精品| 最新的欧美精品一区二区| 黑人巨大精品欧美一区二区蜜桃| 国产成人啪精品午夜网站| 久久鲁丝午夜福利片| 国产精品熟女久久久久浪| 中文欧美无线码| 亚洲成人免费av在线播放| 中国国产av一级| 大陆偷拍与自拍| 亚洲五月色婷婷综合| 久久久精品国产亚洲av高清涩受| 亚洲精品久久成人aⅴ小说| 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 丝袜脚勾引网站| 亚洲一卡2卡3卡4卡5卡精品中文| 男人添女人高潮全过程视频| 亚洲色图综合在线观看| 欧美精品一区二区大全| 国产爽快片一区二区三区| www.精华液| 日韩,欧美,国产一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产福利在线免费观看视频| 欧美亚洲日本最大视频资源| 99久久精品国产亚洲精品| 国产亚洲精品第一综合不卡| 两个人看的免费小视频| 国产有黄有色有爽视频| 妹子高潮喷水视频| 狂野欧美激情性xxxx| 久久av网站| 我要看黄色一级片免费的| 大香蕉久久成人网| 国产亚洲av片在线观看秒播厂| 丰满迷人的少妇在线观看| 精品国产乱码久久久久久男人| 国产亚洲一区二区精品| 十分钟在线观看高清视频www| 久久av网站| 欧美性长视频在线观看| 美女主播在线视频| 欧美黑人精品巨大| 99久久综合免费| 国产成人影院久久av| 国产亚洲av片在线观看秒播厂| 黄片播放在线免费| 在线看a的网站| 国产欧美日韩一区二区三区在线| 熟女少妇亚洲综合色aaa.| 看十八女毛片水多多多| 亚洲精品成人av观看孕妇| 一级毛片 在线播放| 高清视频免费观看一区二区| www日本在线高清视频| 亚洲一卡2卡3卡4卡5卡精品中文| h视频一区二区三区| 色播在线永久视频| 在线观看免费日韩欧美大片| 免费少妇av软件| 啦啦啦在线免费观看视频4| 日韩人妻精品一区2区三区| 午夜精品国产一区二区电影| 激情五月婷婷亚洲| 欧美激情高清一区二区三区| 另类精品久久| 熟女少妇亚洲综合色aaa.| 国产免费现黄频在线看| 亚洲熟女毛片儿| 国产精品亚洲av一区麻豆| 欧美日本中文国产一区发布| 黄色视频在线播放观看不卡| 精品国产一区二区久久| 午夜福利,免费看| 成年人午夜在线观看视频| 少妇人妻 视频| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美软件| 精品熟女少妇八av免费久了| av视频免费观看在线观看| 午夜福利影视在线免费观看| 久9热在线精品视频| 80岁老熟妇乱子伦牲交| 一本综合久久免费| 男人操女人黄网站| 国产av一区二区精品久久| 一本久久精品| 最黄视频免费看| 1024香蕉在线观看| 中文字幕精品免费在线观看视频| 18在线观看网站| 狠狠精品人妻久久久久久综合| e午夜精品久久久久久久| 在线观看一区二区三区激情| 亚洲精品国产一区二区精华液| 亚洲视频免费观看视频| 久久久久久久精品精品| 视频区欧美日本亚洲| 又黄又粗又硬又大视频| 日韩av不卡免费在线播放| 久久久久久免费高清国产稀缺| 精品国产一区二区久久| 成年动漫av网址| 黄色怎么调成土黄色| 久久中文字幕一级| 亚洲精品一卡2卡三卡4卡5卡 | 久久99一区二区三区| 久久精品久久久久久久性|