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

    An Incursion of O ff-EquatorialSubsurface Cold Water and Its Role in Triggering the“Double Dip”La Ni?na Event of 2011

    2015-04-20 05:59:10ZHENGFeiFENGLishaandZHUJiang
    Advances in Atmospheric Sciences 2015年6期

    ZHENG FeiFENG Lishaand ZHU Jiang

    1International Center for Climate and Environment Science,Institute ofAtmospheric Physics, Chinese Academy of Sciences,Beijing100029

    2College of Atmospheric Sciences,Chengdu University of Information Technology,Chengdu610225

    An Incursion of O ff-EquatorialSubsurface Cold Water and Its Role in Triggering the“Double Dip”La Ni?na Event of 2011

    ZHENG Fei?1,FENG Lisha1,2,and ZHU Jiang1

    1International Center for Climate and Environment Science,Institute ofAtmospheric Physics, Chinese Academy of Sciences,Beijing100029

    2College of Atmospheric Sciences,Chengdu University of Information Technology,Chengdu610225

    Based on Global Ocean Data Assim ilation System(GODAS)and NCEP reanalysis data,atmospheric and oceanic processes possibly responsible for the onsetof the 2011/12 La Ni?na event,which followed the 2010/11 La Ni?na even—referred to as a“double dip”La Ni?na—are investigated.The key mechanisms involved in activating the 2011/12 La Ni?na are illustrated by these datasets.Results show that neutral conditions were already evident in the equatorial eastern Pacific during the decaying phase of the 2010/11 La Ni?na.However,isothermalanalyses show obviously cold water still persisting at the surface and atsubsurface depths in off-equatorial regions throughoutearly 2011,being mostpronounced in the tropicalSouth Pacific.The negative SST anomalies in the tropical South Pacific acted to strengthen a southern w ind across the equator.The subsurface cold water in the tropical South Pacific then spread northward and broke into the equatorial region at the thermocline depth.This incursion process ofoff-equatorialsubsurface cold watersuccessfully interrupted the eastern propagation of warm water along the equator,which had previously accumulated atsubsurface depths in the warm pool during the 2010/11 La Ni?na event.Furthermore,the incursion process strengthened as a resultof the off-equatorialeffects,mostly in the tropical South Pacific.The negative SST anomalies then reappeared in the centralbasin in summer 2011,and acted to trigger local coupled air–sea interactions to produce atmospheric–oceanic anomalies thatdeveloped and evolved w ith the second cooling in the fallof 2011.

    “double dip”La Ni?na,subsurface cold water,incursion,off-equatorialprocess

    1.Introduction

    ENSO is the strongestsignalofair–sea interaction,which periodically develops in the tropical Pacific and exerts a surprisingly significant influence on short-term global climate change(e.g.,Huang and Wu,1989;Neelin et al., 1998;Alexanderand Scott,2002;Toniazzo and Scaife,2006; Philippon et al.,2012).The significance of ENSO has been documented by many researchers(e.g.,Wang and Picaut, 2004;M cPhaden et al.,2006;Wang et al.,2012).Recently, ENSO has been receiving increased attention due to the appearance of a new type of El Ni?no,the Central-Pacific(CP) El Ni?no(e.g.,Kao and Yu,2009;Yu etal.,2010),which has maximum SST anomalies centered near the dateline.

    As many efforts have been made to develop and enforce the performance of ENSO models(e.g.,Jietal.,1996;Latif etal.,1998;Kleeman and Moore,1999;Zhang etal.,2003; Zheng et al.,2006a,2006b;2007;Zheng and Zhu,2008; Jin etal.,2008;Zheng etal.,2009;Zheng and Zhu,2010a, 2010b;Barnston etal.,2012),theirskills in predicting ENSO have gradually improved.However,the real-time ENSO prediction skillof models in the pastdecade has been somewhat lower than the less advanced models of the 1980s and 1990s (Barnston etal.,2012).Additionally,Barnston etal.(2012) suggested an apparent retrogression in skillexists for ENSO prediction in the 2002–11 study period compared w ith that in the 1981–2010 period.The 9-yr sliding correlation for the hindcasts over the 1981–2010 period has an average score of 0.65 ata 6-month lead time,but the prediction skill of correlation decreases to 0.42 for the 2002–11 period.There are two possible explanations for this abnormal phenomenon in the pastdecade:one is due to the more frequentoccurrences of CP El Ni?no(Lee and M cPhaden,2010;M cPhaden et al., 2011;Xiang etal.,2012);and the other is the appearance ofa specialkind of La Ni?na(Zhang etal.,2013;Hu etal.,2014), which follows a previous La Ni?na event(e.g.,the 2011/12 event),and is referred to as“double dip”La Ni?na hereafter in this paper.

    M cPhaden and Zhang(2009)argued that there is a ten-dency for some La Nin?a events to persistand re-emerge for a second cooling year after 1980.However,most ENSO models are unable to successfully predict this kind of second cooling eventuntilnow(Barnston etal.,2012).Specifically, among allENSO prediction models(14 dynam icaland 8 statistical)listed on the InternationalResearch Institute for Climate and Society(IRI)website athttp://iri.columbia.edu/climate/ENSO/currentinfo/SST table.htm l,most failed to forecastthe Nin?o3.4 SST anomalies of the 2011/12 event,starting from early-m id 2011,and only one model[an intermediate coupled model(ICM)used by Zhang et al.(2013)]gave a successful prediction of 2011 cold SST conditions over the tropical Pacific one year in advance.This presents a critical challenge to the ENSO prediction community(Barnston et al.,2012)and suggests an imperative demand to understand the physical processes leading to the second-year cooling in 2011.Zhang et al.(2013)demonstrated that thermocline feedback is a crucial mechanism affecting the interannual SST evolution of the“double dip”in 2011,which is specifically represented by the relationship between the oceanic entrainment temperature(Te)and sea level(SL)in the ICM. At the same time,the persistentand sufficiently strong easterly w ind anomaly over the western-central equatorial Pacific,which operates to sustain the large negative subsurface thermal variability off and on the equator to significantly affect the SST evolution during 2010–12 associated w ith the thermocline feedback,has been demonstrated to be another key process by Zhang etal.(2013).Hu etal.(2014)argued there can be two preconditions responsible for the continuation of a La Ni?na event:one,which is necessary but insufficienton its own,is that the previous La Ni?na mustbe a strong event;and the other is whether there are eastward propagating downwelling warm equatorial Kelvin waves during the decaying phase of the fi rst La Ni?na.Moreover,it has been suggested that the subtropical–tropical teleconnection is important and cannot be ignored for the occurrence of ENSO events(Vimontetal.,2001;Wang and Yang,2014).Vimont etal.(2001)proposed a seasonal footprinting mechanism to explain how m idlatitude atmospheric variability could lead to tropical ENSO variability via atmosphere–ocean coupling in the subtropics.Wang and Yang(2014)investigated the impactof subtropical Pacific sea surface temperature anomalies(SSTAs)on the equatorialocean,and proposed a w ind–evaporation–SST mechanism as playing the dominant role in equatorialSST changes.However,allpreviousstudies on the 2011/12 La Ni?na event lack a detailed description of its onsetand evolution due to the qualification of observations and oceanic 3D monitoring,which is crucial in order to understand the mechanism causing the second-year cooling during the long-lived 2010–12 La Ni?na event.

    In this work,to elucidate the trigger mechanism of the 2011/12 La Ni?na event,we exam ine the atmospheric and oceanic processes related to the evolution of the 2011/12 La Ni?na in detail,through analyzing the 3D reanalysis fields from the GlobalOcean Data Assimilation System(GODAS) and its calculated diagnostic variables.We focus on the roles played by off-equatorialsubsurface temperature anomalies in the tropical South and North Pacific through exam ining the subsurface temperature evolution of isotherms.In order to observe the pathways consistently throughout the basin from the off-equatorial to the equatorial region,isothermal analyses are performed using the 3D temperature fi elds.Seeing as the subsurface temperature anomalies tend to propagate along isothermalsurfaces,an isothermalanalysisisnecessary and can better characterize the temporalevolution and the 3D structure of the upper ocean.The major finding in this work is that the enhanced southeasterly w ind anomalies over the tropical Pacific from March to June in 2011,associated w ith the negative SST anomalies in the tropical South Pacific that were leftbehind by the 2010 moderate La Ni?na event,played a key role in spreading the subsurface cold water over the tropical South Pacific northward.The subsurface cold water then broke into the equatorial region at thermocline depth to activate the onsetof the second cooling in summer 2011.In addition,analyzing the temperature anomalies at isothermal surfaces makes itclear thatthe incursive pathway of cold subsurface waterwas along the thermocline depth and originated from off-equatorial regions,mostly in the tropical South Pacific.

    The remainderof the paper is organized as follows.Section 2 briefl y introduces the datasetand analysis method used in this work.Section 3 describes the onset and evolution of the 2011/12 La Ni?na event in detail.Section 4 is concerned w ith the incursion process from the extratropical subsurface cold water into the equatorial region.And finally,a discussion and conclusions are presented in section 5.

    2.Datasetand analysis method

    The monthly mean data for temperature and SL are from GODAS(Behringer and Xue,2004),which are operationally provided by the National Centers for Environmental Prediction(NCEP).The horizontal resolution in the zonal and meridional direction of GODAS is 1?× (1/3)?and the products have 40 levels w ith a 10 m resolution in the upper 200 m.The monthly mean atmospheric reanalysis data,including sea level pressure(SLP)and w ind stress,are from the NECP–DOE(National Centers for Environmental Prediction–Departmentof Energy)Reanalysis II(Kanam itsu etal.,2002),w ith a resolution of2.5?×2.5?and a T62 Gaussian grid(192×94),respectively.Both the GODAS products and atmospheric NCEP reanalysis datasets cover the period from January 1982 to December2011,and the long-term climatological fields are generated based on the period 1982–2004.

    It is well known that the dynam ic and thermodynam ic processes at subsurface layers(e.g.,the m ixing and entrainment)contribute to interannual SST variations(e.g.,Zebiak and Cane,1987;Wang and M cPhaden,2000).Moreover,the thermocline and its variations play an important role in affecting ENSO evolution,butobservation systems for subsurface structure have been very scarce in the past(Wang etal., 2012).The available 3D GODAS products provide a uniqueopportunity to investigate the 3D structure of the upperocean and ENSO evolution in more detail.Furthermore,it is also worth mentioning that the 3D structure of the temperature field in the GODAS products—especially the distributions of the subsurface temperature,which play a critical role in ENSO evolution—is highly correlated w ith thatof the Argo datasetover the tropicalPacific.In this study,due to the poor quality of the salinity fields in GODAS,instead of the isopycnal analyses,isothermal analyses are performed using the 3D temperature fields to better characterize the temporalevolution and 3D structure over the tropical Pacific during the 2011/12“double dip”La Ni?na event.The thermocline depth is defined as the depth where the 20?C isotherm is located, and the climatologicaland interannually varying temperature at the thermocline depth are also calculated.

    3.Onsetand evolution of the 2011/12 La Ni?na

    Various GODAS and NCEP reanalysis data involving the atmospheric and oceanic processes responsible for the onset and evolution of the 2011/12 La Ni?na eventare investigated in this section.Figure 1 shows the temporalevolutions of the interannual variations in SST,SL,zonal w ind stress(Taux), meridional w ind stress(Tauy),and thermocline depth in the equatorial Pacific in 2011.The developmentof the 2011/12 La Ni?na event is quite clear at the equator,and there is a close relationship among these anomalous fields.For the SST anomalies,there is a shiftof the cold–warm–cold conditions over the equatorial central-eastern Pacifi c during the whole of 2011(Fig.1a).A second cooling at the sea surface occurs in m id-to-late 2011,follow ing the fi rst cooling in the fallof 2010(notshown).The 2010 moderate La Ni?na event seems to end in spring 2011(Zhang etal.,2013).Then,the equatorialeastern Pacific experiences neutralSST conditions and w itnesses this weak La Ni?na event(i.e.,“double dip”La Ni?na event)in August 2011.The weak La Ni?na conditions gradually develop to moderate intensity during November 2011.The key feature here is the appearance of the second cooling at the sea surface in the boreal fallof2011.Corresponding to the characteristics of the SST anomalies,the anomalous SL and thermocline depth show sim ilar variation patterns(negative–positive–negative)over the central to eastern equatorial Pacific(Figs.1b and e).Meanwhile,the combined data presented in Figs.1c and d,which depict the easterly and southerly w ind anomalies,respectively,demonstrate a predom inance of the anomalous southeasterly w inds over the whole of the equatorial region.In addition,the depth of the thermocline(Fig.1e),represented by the 20?C isotherm depth,and the SST(Fig.1a),are highly correlated,w ith a cold SST associated w ith a shallow thermocline,and vice versa.Moreover,from February to May 2011,the positive oceanic heat anomalies in the equatorial western Pacific ex-cite a downwelling oceanic Kelvin wave,which are depicted clearly by the thermocline depth anomalies(Fig.1e).Associated w ith the eastward propagating Kelvin wave,a pronounced warm ing process exists in the equatorial centraleastern Pacific during early–m id 2011(Fig.1a).However, after the Kelvin wave successfully reaches the eastern boundary of the tropical Pacific,the warm ing processdoes notcontinue afterward and is replaced by a cooling process in the whole equatorial Pacific.

    Fig.1.Temporal evolutions of interannualanomalies along the equatorial Pacific(averaged between 2?S and 2?N)in 2011 for(a)SST,(b)SL,(c)Taux,(d)Tauy,and(e)thermocline depth.The contour interval is 0.5?C in(a),5 cm in(b),0.4 dyn cm-2in(c)and(d),and 9 m in(e).

    Figure 2 displays the horizontaldistributions of the interannual SST anomalies over the equatorial Pacific region in 2011.The vigorous period of the 2010 La Ni?na persists to January 2011,so the negative SST anomalies prevail extensively in the tropical Pacific,w ith a m inimum of less than -2?C along the equator between 150?and 170?W(Fig.2a). From March to June,the region of negative SST anomalies shrinks gradually in the equatorialeastern Pacific,while a small-scale warm SST signal appearing in March gradually expands its area and increases its strength during the next three months(Figs.2b–e).As a result,neutral SST conditions emerge in the eastern equatorial Pacific in m id-2011.During this period,negative SST anomalies persist over off-equatorial regions,being most pronounced in the tropical South Pacific.Thereafter,the domain of positive SST anomalies gradually dim inishes and is displaced by negative ones.In particular,the negative SST anomalies reappear in the equatorial central Pacific in August(Fig.2g),and this cooling condition tends to persistand strengthen in the follow ing months(Fig.2h),leading to a“double dip”La Ni?na event.The condition shift of the SST anomalies(i.e., negative–positive–negative)in the central-eastern Pacific is the most direct response to this kind of continuous La Ni?na event.

    In contrast to the variations of temperature at the surface layer,those in the subsurface layers are clearer.Figure 3 presents the horizontal distributions of interannual temperature anomalies at the thermocline depth,associated w ith the current anomalies.In January,the cold and warm water at the subsurface layers over the tropical Pacific nearly separate at 160?W.The cold subsurface water is located in the eastern tropical Pacific,while the warm subsurface water is located in the western tropical Pacific(Fig.3a).The positive subsurface temperature anomalies,which previously accumulated atsubsurface depthsin the tropicalwestern Pacific (Fig.3a),startto propagate along the equator toward the eastern basin in February and March(Fig.3b).Simultaneously, pronounced negative anomalies persist in the subsurface layers of the off-equatorial regions,and warm water occupies almost the entire equatorial Pacific.Note that the cold water appears unexpectedly in 145?–155?W and 2?N–2?S and interrupts the eastern propagation of the warm water along the equator in April(Fig.3c).Then,the region of warm water shrinks gradually in the eastern equatorial Pacifi c and the area of equatorial cold water expands gradually in the nextfew months(Figs.3d–h).Finally,the cold water is seen to totally occupy the tropicaleastern Pacific again(Figs.3d and h).The pattern of the subsurface temperature anomalies in October is similar to that in January,only w ith a weakeramplitude.Obviously,the temperature anomalies at the thermocline depth(Fig.3)can reflect the evolution of the second cooling eventbetter than the SST anomalies(Fig.2).For the subsurface current field,a pathway originating from the tropical South Pacific can be clearly seen:the cold water carried by the southeastern currents is eventually transported to the equator(Fig.3).

    Fig.2.Horizontal distributions of the interannual SST anomalies over the equatorial Pacific region in some selected months of 2011.The contour interval is 1?C.

    Fig.3.The same as in Fig.2,but for the temperature(shading)and ocean current(vectors)anomalies at the thermocline depth.The units forocean current is cm s-1.

    Due to the strong interaction between ocean and atmosphere over the tropical Pacific,atmospheric components also play significant roles in the development of the“double dip”La Ni?na event.Here,the NCEP reanalysis data are used to describe the atmospheric processes closely related to this event.Figure 4 shows the horizontal distributions of the SLP anomalies(shading)and w ind stress anomalies(vectors).Consistentw ith the negative SST anomalies,the center of the positive SLP anomalies exists for several months in the tropical central-eastern Pacific during the 2011 La Ni?na event,acting to strengthen a southern meridionalw ind across the equator in early–m id 2011,especially from February to May(Figs.4b–d).The gradientof the SLP leads to an equatorward transportof cold water,and acts to produce the crossequatorialsouthern w ind.The cross-equatorialsouthern w ind can make the cold surface water in the tropical South Pacific expand northward to the equator and m ix w ith the cold surface water in the tropical North Pacifi c,off-equator.At the same time,the southeasterly w ind can cause the stratification variation of the thermocline to furtherspread the cold subsurface water to the equator from the off-equatorial regions.

    During the decaying phase of the 2011/12 La Ni?na,the region of negative SST anomalies gradually shrinks in the equatorialeastern Pacific,but these cooling processes do not persistand develop,perhaps because the off-equatorial cold anomalies(i.e.,South and North Pacific)are too weak to provide enough cold water,and the meridionalw inds also cannot supportan equatorward current in 2012(fi gures notshown). Finally,the SSTA does not return to the La Ni?na state,as happened during 2011.

    Fig.4.The same as in Fig.2,but for the sea levelpressure anomalies(shading)and the w ind stress anomalies (vectors).

    Figure 5 shows the vertical distribution of the interannual temperature anomalies in the upper oceans(i.e.,from the surface to 300 m)over the equatorial region.In January 2011,positive temperature anomalies can be observed at the surface in the tropical western Pacific,while negative temperature anomalies are located in the tropical central-eastern Pacific regions at the surface layer.However,positive temperature anomalies and negative temperature anomalies exist in almostequalmeasure in the subsurface layers,w ith a sharp temperature frontalong 160?W.Meanwhile,the thermocline over the equatorial western Pacific is relatively deeper than thatin the equatorialeastern Pacific(Fig.5a).Thereafter,the positive subsurface temperature anomalies,which had previously accumulated at subsurface depths over the tropical western Pacific,start to propagate toward the eastern basin along the equator,acting to dim inish the negative temperature anomalies there,reverse the sign of SST anomalies,and initiate warm conditions(Fig.5b).In April,the negative temperature anomalies in the equatorial central-eastern Pacific almost totally disappear except for a small surface region between the dateline and 140?W,and the whole surface equatorial Pacific is almost completely occupied by neutralwarm SST conditions(Fig.5c).In May,a striking feature is evident in the form of a cold water mass unpredictably appearing and persisting at the thermocline depth between 140?and 150?W,when the warm wateroccupies almost the entire equatorial Pacific both at the surface and in the subsurface layers(Fig.5d).In the follow ing months,the range of this cold water mass gradually expands and eventually reoccupies the whole equatorialcentral-eastern Pacific(Figs.5e–h). Based on these analysis results,we can conclude that this incursion of cold water into the equatorial region took place at the thermocline depth,and the incursion process was the key mechanism preventing the developmentof a warm eventand causing a second-year cooling during the 2010–12 La Ni?na event.However,it is still difficult to distinguish whether the source of the cold water was from the south or the north,as wellas whatkind of physicalprocess was able to trigger the incursion of cold water.These questions are furtherexam ined in the follow ing section.

    4.The incursion process from the offequatorial subsurface water into the equatorial region

    The incursion process is clearly displayed in Figs.6a–d by show ing the horizontal distributions of the interannual temperature anomalies at the thermocline depth in a specifically smaller region(10?S–10?N,165?W–120?W)and a shorter period(from March to June).In March,the negative subsurface anomalies still prevailover the off-equatorial region,while the equatorial Pacific is occupied by warm water in the subsurface layers(Fig.6a).Note that the cold subsurface water,which is located in the tropical South Pacific,starts to spread northward and break into the equatorial region between 145?and 155?W,successfully interrupting the eastern propagation of warm water along the equator in April(Fig.6b).In May,the incursion of the cold subsurface water from the tropical South Pacific into the equatorial region strengthens,and the cold water tends to move across the equator to the north(Fig.6c).Subsequently,the crossequatorial cold subsurface water from the tropical South Pacific mixes w ith that located in the tropicalNorth Pacific,and the off-equatorial negative temperature anomalies converge over the equator(Fig.6d).Furthermore,the cold waters in the tropicalSouth and North Pacific(the two selected specific regions)are previously accumulated at the subsurface depths. Consistent w ith the surface w inds,the subsurface currents,combined w ith the persistent cold water at the subsurface depths over the tropical South Pacific,are active enough to generate a cold water pathway to the equatorial region.The area of the junction gradually expandsin the next few months and the cold water can again be seen to completely occupy the tropical eastern Pacific(Fig.3).As depicted above,the cold water in this incursion process is mostly from the offequatorial regions(i.e.,the tropical South Pacific).

    Fig.5.Vertical distributions of the interannual temperature anomalies in the upper oceans of the equatorial region(averaged between 2?S and 2?N)in 2011 for(a)January,(b)March,(c)April,(d)May,(e) June,(f)July,(g)August,and(h)October.The green line represents the thermocline depth diagnosed by the 20?C isotherm.

    Fig.6.Horizontal distributions of the interannual temperature(shading)and ocean current(vectors)anomalies at the thermocline depth over the central to eastern equatorial Pacific(10?S–10?N,165?–120?W)in 2011 for(a)March,(b)April,(c)May,and(d)June.The contour interval is 2?C,and the units for ocean current is cm s-1.The bold black curve represents the zero-line.

    Fig.7.Temperature tendency in 2011 over the incursion region (150?–140?W,2?S–0?)at the thermocline depth.The black line represents the total temperature tendency;the red line represents the effect of climatological current w ith temperature anomalies on the temperature tendency;and the blue line represents the effectof anomalous currentw ith temperature anomalies on the temperature tendency.

    Consistent w ith the demonstrations in Huang et al. (2010),the climatological currentw ith anomalous temperature plays a more important role in affecting the tropical temperature tendency than the anomalous current w ith temperature anomalies.In 2011,the climatological current associated w ith temperature anomalies provides a cooling effecton modulating the temperature tendency over the tropical South Pacific at the thermocline depth(Fig.7).Especially during the key MAM-tim ing(March–April–May tim ing)for the onsetof“double dip”La Nina,the cooling effectis obvious and significant in affecting the temperature tendency.In contrast, over the whole of2011,the anomalouscurrenthas little effect on enhancing the temperature tendency.

    In order to investigate the 3D structure of the pathway along which the off-equatorial cold water moves into the equatorial region,we compute an oblique section from(2?N, 150.5?W)to(8?S,141.5?W).The vertical distributions of the interannual temperature anomalies in the upper oceans (5–300 m)from(2?N,150.5?W)to(8?S,141.5?W)are illustrated in Fig.8.It is clear that the negative temperature anomalies during the incursion process almost break into the equatorial region at the thermocline depth(Figs.8a–d), which has already been discussed based on Figs.3 and 6. In March,in the subsurface layers,there are negative temperature anomalies over the south of the equator and positive temperature anomalies over the equatorial region(Fig.8a).A northward position of the temperature zero-line in the subsurface depths occurs in April,and the positive temperature anomalies at the equator dissipate at that time(Fig.8b).In May,the negative subsurface temperature anomalies in the tropical South Pacific start to break into the equatorial region along the thermocline via horizontal advection(Fig.8c),and then mix w ith the negative subsurface temperature anomalies located in the tropicalNorth Pacific region(Fig.8d).

    As described in section 3,consistent w ith the negative SST anomalies over the tropical South Pacific,the positive SLP anomalies over that area act to produce a southern meridional w ind across the equator in early–m id 2011.As the interannual variations of the upper ocean over the tropical Pacific are controlled by the surface w ind(e.g.,Wang and M cPhaden,2000),a potential trigger for the activation of the incursion process of the subsurface cold water can be found in the form of the enhanced southern meridional w ind throughoutearly–mid 2011.Zhang etal.(2013)argued that interannualw ind forcing effects were important to SST evolution during 2011.Figure 9 shows the zonal–time sections for anomalous w ind stress changes during 2011 to illustrate the key factor triggering the 2011/12 La Ni?na event. In Fig.9a,the southeasterly w ind reappears over the tropical South Pacific(100?–110?W)around March,and spreads to almost the entire tropical South Pacific(160?E–80?W)from April.Moreover,the w ind is crossing the equator(Fig.4). On the other hand,southerly anomalies exist over the tropical North Pacific(Fig.9b),which are not conducive to the invasion of the subsurface cold water from the north of the equator.All these phenomena could have contributed to the incursion process of the off-equatorial subsurface cold water having mostly arisen from the tropical South Pacific,indicating significant roles played by the w ind in activating the“double dip”La Ni?na event.It is also worth noting that the subsurface cold water in the North Pacific is also significant. However,the surface cold water over that region is weaker than that in the South Pacific,and there are no corresponding northeasterly w inds to cause the stratification variation of the upperoceansto furtherspread the cold subsurface water from the north to the equator.

    5.Discussion and conclusions

    Fig.8.Verticaldistributions of the interannual temperature anomalies in the upperoceans from(2?N,150.5?W) to(8?S,141.5?W)in 2011 for(a)March,(b)April,(c)May,and(d)June.The contour interval is 1?C.The green line represents the thermocline depth diagnosed by the 20?C isotherm.The bold black curve represents the zero-line,and the verticalorange line represents the equator.

    The atmospheric and oceanic processes that could have been responsible for the onset of the La Ni?na event of 2011 have been investigated,and an incursion process of subsurface cold water from off-equatorial regions into the equatorial region has been illustrated by GODAS 3D temperature productsand NCEP reanalysisatmospheric datasets.Furthermore,ithas been demonstrated that this incursion process of cold subsurface waterwas the key process in precipitating the“double dip”La Ni?na.During 2010–12,a prolonged La Ni?na eventarose in the tropicalPacific,w ith a second cooling atthe sea surface follow ing an initialcooling in the fallof2010.After the moderate La Ni?na event in 2010,positive subsurface anomalies,which had previously accumulated at subsurface depths over the tropicalwestern Pacific,started to propagate along the equator toward the eastern basin,acting to reverse the sign of the SST anomalies there and initiate warm conditions.Then,in m id-2011,near normal SST conditions were present in the equatorial eastern Pacific.However,throughoutearly 2011,obviously negative anomalies persisted at the surface and at subsurface depths in off-equatorial regions, being most pronounced in the tropical South Pacific.The negative SST anomalies in the tropical South Pacific acted to strengthen a southern meridional w ind across the equator.The subsurface negative temperature anomalies in the tropical South Pacific then spread northward and broke into the equatorial region at the thermocline depth.This incursion process ofoff-equatorialsubsurface cold watersuccessfully interrupted the eastern propagation of the warm water along the equator,and was strengthened as a resultof the offequatorial effects,mostly in the tropical South Pacific.The cold SST anomalies then reappeared in the central basin in summer 2011,and acted to trigger local coupled air–sea interactions producing atmospheric–oceanic anomalies thatdeveloped and evolved w ith the second cooling in the fall of 2011.

    The evolution of the 2011“double dip”La Ni?na eventdescribed in this study is unprecedented,and depended on the obvious negative temperature anomalies thatpersisted at the surface and at subsurface depths over off-equatorial regions, as wellas the triggering mechanism of the atmosphere.During the decaying phase of the 2010 La Ni?na event,the surface and subsurface cold water persisted over subtropical regions,and the subsurface cold water maintained a stronger negative value and retreated slower than that at the surface (Feng etal.,2015).Moreover,different from the mechanism of subtropical–tropical teleconnection,the off-equator cold anomaly in early–m id 2011 may nothave been conducive to the equatorial recharge process,and as a result could have favored the persistence of a cold ocean subsurface temperature anomaly and prevented the transition from La Ni?na to El Ni?no(Hu et al.,2014).Two issues have been demonstrated as key conditions responsible for the appearance of“double dip”La Ni?na.One is subsurface cold water stillexisting in off-equatorial regions during the decaying phase of the previous La Ni?na event.Accumulated cold water in the off-equator subsurface layers provides a precondition for the generation ofa“double dip”La Ni?na event.The other is persistent negative SST anomalies in the tropical South Pacific acting to produce and strengthen southern w inds across the equator.The enhanced southern w ind then plays a triggering role thatencourages subsurface cold water in the South Pacific to spread northward and break into the equatorial region at the thermocline depth.

    Fig.9.Zonal–time sections of the w ind stress anomalies in 2011:(a)south of the equator(averaged between 5?S and the equator);(b)north of the equator(averaged between the equator and 5?N).

    The fi ndings of our work can also partly explain why most models,especially coupled atmosphere–ocean general circulation models,failed to predict the“double dip”La Ni?na event that started from early–mid 2011.The failure possibly resulted from the initialization in such models.In this paper,we indicate that the persistent southeasterly w inds and prescribed subsurface cold conditions w ith an equatorward incursion would have been the key factors for predicting the 2011/12 La Ni?na.However,for most ENSO models,their initialization methods are based on“ocean-only”assim ilation schemes(i.e.,only oceanic observations are assim ilated into the model,and the initial atmospheric states are keptunchanged).Moreover,due to the lim itations of the quality and quantity of ocean currentobservations,it is very hard to assimilate the observed current directly(Zheng and Zhu,2010a).Thus,providing more accurate atmospheric and oceanic subsurface information in the initial conditions ofmodels based on more advanced coupled data assim ilation methods(e.g.,Zheng and Zhu,2008,2010a)would clearly benefi t seasonal-to-interannual predictions of“double dip”La Ni?na events.

    However,these conclusions have only been drawn based on one“double dip”La Ni?na case.There are other“double dip”La Ni?na cases that have happened over the tropical Pacific(e.g.,in 2008/09 and 2000/01),and so the similarities and differences among these events need to be analyzed to better describe the nature of these strikingly different ENSO evolutions associated w ith various physical processes w ithin the Pacific climate system(Hu etal.,2014).For instance—as exemplified by the fact thata sim ilar incursion of subsurface cold water from the North Pacific happened in 2008,while there was no obvious incursion process from the off-equatorial regions in 2000—whetherornot the persistent subsurface cold watercan break into the equatorial region can differ from one case to the next.Furthermore,the differences in the location and orientation of the incursion process vary and depend on the atmospheric surface conditions(i.e.,easterly and equatorward w inds).

    Moreover,due to the lack of 3D oceanic data in the past, modeling studies are also needed to validate the roles played by the off-equatorialsurface and subsurface anomalies in activating“double dip”La Ni?na events over the tropical Pacific.For example,approximately 15“double dip”La Ni?na cases could be simulated by a long-term(i.e.,500 years) pre-industrial control run in a fully coupled atmosphere–ocean model(Zheng,2014),which could help analyze the key mechanisms favoring the occurrence of“double dip”La Ni?na events through providing more than 10 samples of the second cooling case.

    Acknow ledgements.We thank the two anonymous reviewers for their valuable comments.This work was supported by the National Program for Supportof Top-notch Young Professionals,the National Basic Research Program of China(Grant No. 2012CB955202),the Chinese Academy Sciences’Project“Western Pacific Ocean System:Structure,Dynam ics and Consequences”(WPOS;Grant No.XDA10010405),and the National Natural Science Foundation of China(GrantNo.41176014).

    REFERENCES

    A lexander,M.,and J.Scott,2002:The influence of ENSO on airsea interaction in the Atlantic.Geophys.Res.Lett.,29(14), 1701,46-1–46-4,doi:10.1029/2001GL014347.

    Barnston,A.G.,M.K.Tippett,M.L.L’Heureux,S.Li,and D.G. DeWitt,2012:Skill of real-time seasonal ENSO model predictions during 2002–11,is our capability increasing?Bull. Amer.Meteor.Soc.,93,631–651,doi:10.1175/BAMS-D-11-00111.1.

    Behringer,D.,and Y.Xue,2004:Evaluation of the global ocean data assim ilation system at NCEP:The Pacific Ocean. Preprints,Eighth Symposium on Integrated Observing and Assimilation Systems for Atmosphere,Oceans,and Land Surface,Washington State Convention and Trade Center,Seattle, Washington,Amer.Meteor.Soc.

    Feng,L.,R.-H.Zhang,Z.Wang,and X.Chen,2015:Processes leading to the second-year cooling of the 2010–12 La Ni?na event,diagnosed using GODAS.Adv.Atmos.Sci.,doi: 10.1007/s00376-014-4012-8.

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

    Huang,B.,Y.Xue,X.Zhang,A.Kumar,and M.J.M cPhaden, 2010:The NCEP GODAS ocean analysis of the tropical Pacific mixed layer heatbudgeton seasonal to interannual time scales.J.Climate,23,4901–4925.

    Huang,R.H.,and Y.F.Wu,1989:The influence of ENSO on the summer climate change in China and its mechanisms.Adv. Atmos.Sci.,6(1),21–32,doi:10.1007/BF02656915.

    Ji,M.,A.Leetmaa,and V.E.Kousky,1996:Coupled model forecasts of ENSO during the 1980 and 1990s atthe NationalMeteorological Center for environmental prediction.J.Climate, 9,3105–3120.

    Jin,E.K.,and Coauthors,2008:Currentstatus ofENSO prediction skill in coupled ocean–atmosphere models.Climate Dyn.,31, 647–664,doi:10.1007/s00382-008-0397-3.

    Kao,H.-Y.,and J.-Y.Yu.2009:Contrasting eastern-Pacific and centralPacific types of ElNi?no.J.Climate,22,615–632,doi: 10.1175/2008JCLI2309.1.

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

    K leeman,R.,and A.M.Moore,1999:A new method for determining the reliability of dynamical ENSO predictions.Mon. Wea.Rev.,127,694–705.

    Latif,M.,and Coauthors,1998:A review of the predictability and prediction of ENSO.J.Geophys.Res.,103(C7),143 75–14 393.

    Lee,T.,and M.J.M cPhaden,2010:Increasing intensity of El Ni?no in the central-equatorialPacific.Geophys.Res.Lett.,37, L14603,doi:10.1029/2010GL044007.

    M cPhaden,M.J.,S.E.Zebiak,and M.H.Glantz,2006:ENSO as an integrating concept in earth science.Science,314,1740–1745.

    M cPhaden,M.J.,and X.B.Zhang,2009:Asymmetry in zonal phase propagation of ENSO sea surface temperature anomalies.Geophys.Res.Lett.,36,L13703,doi:10.1029/2009GL 038774.

    McPhaden,M.J.,T.Lee,and D.M cClurg,2011:El Ni?no and its relationship to changing background conditions in the tropical Pacific Ocean.Geophys.Res.Lett.,38,L15709,doi:10.1029/ 2011GL048275.

    Neelin,J.D.,D.S.Battisti,A.C.Hirst,F.-F.Jin,Y.Wakata,T. Yamagata,and S.E.Zebiak,1998:ENSO theory.J.Geophys. Res.,103(C7),14 262–14 290.

    Philippon,N.,M.Rouault,Y.Richard,and A.Favre,2012:The influence of ENSO on w inter rainfall in South Africa.Int.J. Climatol.,32,2333–2347,doi:10.1002/joc.3403.

    Toniazzo,T.,and A.A.Scaife,2006:The influence of ENSO on w inter North Atlantic climate.Geophys.Res.Lett.,33, L24704,doi:10.1029/2006GL027881.

    Vimont,D.,J.M.Wallace,and D.S.Battisti,2001:Footprinting: A seasonalconnection between the mid-latitudes and tropics. Geophys.Res.Lett.,28,3923–3926.

    Wang,C.,and J.Picaut,2004:Understanding ENSO physics—A review.Earth’s Climate:The Ocean-Atmosphere Interaction,C.Wang,S.-P.Xie,and J.Carton,American GeophysicalUnion,Washington,D.C.,21–48.

    Wang,C.,C.Deser,J.-Y.Yu,P.DiNezio,and A.Clement,2012: El Ni?no and Southern Oscillation(ENSO):A review.Coral Reefs of the Eastern Pacific,P.W.Glynn,Spring,3–19.

    Wang,L.,and H.-J.Yang,2014:The role of atmospheric teleconnection in the subtropical thermal forcing on the equatorial Pacific.Adv.Atmos.Sci.,31(4),985–994,doi: 10.1007/s00376-013-3173-1.

    Wang,W.,and M.J.M cPhaden,2000:The surface-layer heatbalance in the equatorialPacific Ocean.Part II:Interannualvariability.J.Phys.Oceanogr.,30,2989–3008.

    Xiang,B.Q.,B.Wang,and T.Li,2012:A new paradigm for the predom inance of standing Central Pacific Warm ing after the late 1990s.Climate Dyn.,41,327–340,doi:10.1007/s00382-012-1427-8.

    Yu,J.-Y.,H.-Y.Kao,and T.Lee,2010:Subtropics-related interannualsea surface temperature variability in the equatorialcentral Pacific.J.Climate,23,2869–2884,doi:10.1175/2010 JCLI3171.1.

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

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

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

    Zheng,F.,2014:Variability simulated by a coupled general circulation model:ECHAM 5/MPI-OM.Atmos.Oceanic Sci.Lett., 7,471–475,doi:10.3878/j.issn.1674-2834.14.0014.

    Zheng,F.,and J.Zhu,2008:Balanced multivariate model errors of an intermediate coupled model for ensemble Kalman fi lter data assimilation.J.Geophys.Res.,113,C07002,doi: 10.1029/2007JC004621.

    Zheng,F.,and J.Zhu,2010a:Coupled assim ilation for an intermediated coupled ENSO prediction model.Ocean Dyn.,60, 1061–1073,doi:10.1007/s10236-010-0307-1.

    Zheng,F.,and J.Zhu,2010b:Spring predictability barrier of ENSO events from the perspective of an ensemble prediction system.Global Planet Change,72,108–117,doi:10.1016/ j.gloplacha.2010.01.021.

    Zheng,F.,J.Zhu,R.-H.Zhang,and G.-Q.Zhou,2006a:Ensemble hindcasts of SST anomalies in the tropicalPacific using an intermediate coupled model.Geophys.Res.Lett.,33,L19604,doi:10.1029/2006GL026994.

    Zheng,F.,J.Zhu,R.-H.Zhang,and G.-Q.Zhou,2006b:Improved ENSO forecasts by assim ilating sea surface temperature observations into an intermediate coupled model.Adv.Atmos. Sci.,23(4),615–624,doi:10.1007/s00376-006-0615-z.

    Zheng,F.,J.Zhu,and R.-H.Zhang,2007:The impactof altimetry data on ENSO ensemble initializations and predictions.Geophys.Res.Lett.,34,L13611,doi:10.1029/2007GL030451.

    Zheng,F.,J.Zhu,H.Wang,and R.-H.Zhang,2009:Ensemble hindcasts of ENSO events over the past 120 years using a large numberofensembles.Adv.Atmos.Sci.,26(2),359–372, doi:10.1007/s00376-009-0359-7.

    :Zheng,F.,L.S.Feng,and J.Zhu,2015:An incursion of off-equatorial subsurface cold water and its role in triggering the“double dip”La Ni?na eventof2011.Adv.Atmos.Sci.,32(6),731–742,

    10.1007/s00376-014-4080-9.

    (Received 15 April2014;revised 19 October2014;accepted 22 October2014)

    ?Corresponding author:ZHENG Fei

    Email:zhengfei@mail.iap.ac.cn

    ?Institute of Atm ospheric Physics/Chinese Academ y of Sciences,and Science Press and Springer-Verlag Berlin Heidelberg 2015

    av福利片在线观看| 搡女人真爽免费视频火全软件 | 在线观看午夜福利视频| 欧美极品一区二区三区四区| 一级毛片女人18水好多| 亚洲成人久久性| 国产免费av片在线观看野外av| 一夜夜www| x7x7x7水蜜桃| 91在线观看av| 午夜福利在线在线| 90打野战视频偷拍视频| 欧美日韩福利视频一区二区| 成人欧美大片| 久久精品人妻少妇| 国产黄a三级三级三级人| 女人十人毛片免费观看3o分钟| 久久久久性生活片| 国产成+人综合+亚洲专区| 国产三级中文精品| av在线天堂中文字幕| 日韩有码中文字幕| 黑人欧美特级aaaaaa片| 欧美高清成人免费视频www| 欧美乱码精品一区二区三区| 国产成人av激情在线播放| 国产97色在线日韩免费| 老司机午夜十八禁免费视频| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看成人毛片| 国产日本99.免费观看| 美女 人体艺术 gogo| 免费在线观看成人毛片| 一级毛片高清免费大全| 12—13女人毛片做爰片一| 深夜精品福利| 午夜精品久久久久久毛片777| 亚洲人成网站高清观看| 欧美一区二区精品小视频在线| 99在线人妻在线中文字幕| 欧美成狂野欧美在线观看| 久久久久精品国产欧美久久久| 亚洲av二区三区四区| 亚洲美女视频黄频| 午夜福利在线观看吧| 欧美极品一区二区三区四区| 在线观看免费视频日本深夜| 久久精品国产清高在天天线| 免费看十八禁软件| 欧美大码av| av福利片在线观看| 久久精品影院6| 欧美成狂野欧美在线观看| 黑人欧美特级aaaaaa片| 大型黄色视频在线免费观看| 偷拍熟女少妇极品色| 国产麻豆成人av免费视频| 国产69精品久久久久777片| 欧美中文日本在线观看视频| 18禁裸乳无遮挡免费网站照片| 亚洲七黄色美女视频| 亚洲熟妇中文字幕五十中出| 久久中文看片网| 丰满的人妻完整版| 五月玫瑰六月丁香| 床上黄色一级片| 亚洲熟妇熟女久久| 99国产综合亚洲精品| 一个人观看的视频www高清免费观看| 在线观看日韩欧美| 观看免费一级毛片| 两个人的视频大全免费| 午夜福利免费观看在线| 久久草成人影院| 天天一区二区日本电影三级| 久久99热这里只有精品18| 午夜老司机福利剧场| 国产色爽女视频免费观看| 全区人妻精品视频| 久久久久久久午夜电影| 国产日本99.免费观看| 成人永久免费在线观看视频| 久久99热这里只有精品18| 最后的刺客免费高清国语| 精品久久久久久久久久久久久| 亚洲专区国产一区二区| 99在线人妻在线中文字幕| 亚洲成人久久爱视频| 真人做人爱边吃奶动态| 欧美在线黄色| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 又粗又爽又猛毛片免费看| 露出奶头的视频| 国产日本99.免费观看| 国产免费男女视频| 搡女人真爽免费视频火全软件 | 窝窝影院91人妻| 久久精品人妻少妇| 日韩av在线大香蕉| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 久久久久久久久久黄片| 在线免费观看不下载黄p国产 | 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 国产精品女同一区二区软件 | 99热6这里只有精品| 欧美日韩国产亚洲二区| 97人妻精品一区二区三区麻豆| 欧美成人性av电影在线观看| 国产aⅴ精品一区二区三区波| 午夜福利在线观看免费完整高清在 | 久久精品国产清高在天天线| 欧美日韩福利视频一区二区| 久久久久免费精品人妻一区二区| 美女免费视频网站| 最近最新免费中文字幕在线| 免费一级毛片在线播放高清视频| 国产黄片美女视频| 成年版毛片免费区| 熟女电影av网| 国产黄片美女视频| 三级国产精品欧美在线观看| 男女床上黄色一级片免费看| 两人在一起打扑克的视频| 久久精品国产亚洲av涩爱 | 一个人观看的视频www高清免费观看| 中文在线观看免费www的网站| 精品一区二区三区人妻视频| av在线天堂中文字幕| 校园春色视频在线观看| 亚洲 欧美 日韩 在线 免费| 嫩草影院精品99| 国产精品久久久久久久久免 | 国内少妇人妻偷人精品xxx网站| 噜噜噜噜噜久久久久久91| 久久精品国产99精品国产亚洲性色| 日韩精品中文字幕看吧| 国产欧美日韩精品亚洲av| 午夜激情欧美在线| 免费av不卡在线播放| 亚洲18禁久久av| aaaaa片日本免费| 又黄又爽又免费观看的视频| 长腿黑丝高跟| 99久久综合精品五月天人人| 麻豆久久精品国产亚洲av| 天天添夜夜摸| 真实男女啪啪啪动态图| 国产av在哪里看| 国内毛片毛片毛片毛片毛片| 欧美乱妇无乱码| 特大巨黑吊av在线直播| 又爽又黄无遮挡网站| 欧美一级毛片孕妇| 中文字幕久久专区| 欧美高清成人免费视频www| 精品一区二区三区视频在线观看免费| 午夜两性在线视频| 丁香欧美五月| 欧美午夜高清在线| 搡老妇女老女人老熟妇| 女警被强在线播放| 久久久久国产精品人妻aⅴ院| 久久精品国产清高在天天线| aaaaa片日本免费| 国产乱人伦免费视频| 日本精品一区二区三区蜜桃| 中出人妻视频一区二区| 老熟妇仑乱视频hdxx| 偷拍熟女少妇极品色| 国产精品永久免费网站| 亚洲精品一区av在线观看| 母亲3免费完整高清在线观看| 亚洲在线观看片| 波多野结衣高清无吗| 欧美在线一区亚洲| 9191精品国产免费久久| 国产91精品成人一区二区三区| 在线观看免费午夜福利视频| 岛国在线观看网站| 久久久国产精品麻豆| 特大巨黑吊av在线直播| 欧美+亚洲+日韩+国产| 欧美在线一区亚洲| 看免费av毛片| 国产精品av视频在线免费观看| 日韩亚洲欧美综合| 日本熟妇午夜| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 亚洲国产高清在线一区二区三| 亚洲狠狠婷婷综合久久图片| 两个人看的免费小视频| 国产精品香港三级国产av潘金莲| 免费在线观看成人毛片| 精品一区二区三区人妻视频| 嫩草影院入口| 国产成+人综合+亚洲专区| 精品乱码久久久久久99久播| www国产在线视频色| 中文亚洲av片在线观看爽| 九九热线精品视视频播放| 欧美日韩亚洲国产一区二区在线观看| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| netflix在线观看网站| 黄色女人牲交| 中出人妻视频一区二区| 日韩欧美国产一区二区入口| 国产精华一区二区三区| 久久久久国产精品人妻aⅴ院| 亚洲av成人av| 国产av一区在线观看免费| 亚洲 国产 在线| 最近最新免费中文字幕在线| 青草久久国产| 欧美黄色淫秽网站| 免费av观看视频| 国内少妇人妻偷人精品xxx网站| 91在线精品国自产拍蜜月 | 黄色女人牲交| aaaaa片日本免费| 久久香蕉精品热| 香蕉丝袜av| 国产精品一区二区免费欧美| 91av网一区二区| 麻豆成人午夜福利视频| 亚洲av五月六月丁香网| 制服丝袜大香蕉在线| 久久久精品大字幕| 国产成人系列免费观看| 中文字幕精品亚洲无线码一区| av福利片在线观看| 两个人视频免费观看高清| 国产精品电影一区二区三区| 久久精品国产清高在天天线| 老司机福利观看| 在线观看av片永久免费下载| 亚洲成av人片免费观看| 黄色日韩在线| 一个人免费在线观看电影| 夜夜躁狠狠躁天天躁| 国产一区二区在线观看日韩 | 中文亚洲av片在线观看爽| 免费观看人在逋| 精品国内亚洲2022精品成人| 日韩欧美三级三区| 精品午夜福利视频在线观看一区| 亚洲 国产 在线| 99riav亚洲国产免费| 久久久久免费精品人妻一区二区| 亚洲熟妇中文字幕五十中出| 最新美女视频免费是黄的| 少妇人妻一区二区三区视频| 精品熟女少妇八av免费久了| 在线天堂最新版资源| 国产伦人伦偷精品视频| 午夜日韩欧美国产| 国产黄a三级三级三级人| 在线播放国产精品三级| av欧美777| 非洲黑人性xxxx精品又粗又长| 日本成人三级电影网站| 丝袜美腿在线中文| 中文资源天堂在线| 亚洲av电影在线进入| 男女那种视频在线观看| 黄色女人牲交| 国产精品久久电影中文字幕| 国产真人三级小视频在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 国产毛片a区久久久久| 久久精品综合一区二区三区| 精品一区二区三区av网在线观看| 色综合亚洲欧美另类图片| 国产精品亚洲av一区麻豆| 免费在线观看影片大全网站| 亚洲精华国产精华精| 无限看片的www在线观看| 女生性感内裤真人,穿戴方法视频| 日韩欧美国产在线观看| 国产免费av片在线观看野外av| 狂野欧美激情性xxxx| 亚洲精品亚洲一区二区| 操出白浆在线播放| 欧美一区二区精品小视频在线| 日韩欧美精品免费久久 | www日本黄色视频网| 无人区码免费观看不卡| 日本免费一区二区三区高清不卡| 最近在线观看免费完整版| 欧美日韩综合久久久久久 | 国产高清激情床上av| 国产精品亚洲av一区麻豆| 日韩人妻高清精品专区| 一级作爱视频免费观看| 人妻夜夜爽99麻豆av| 亚洲精品久久国产高清桃花| 国产乱人伦免费视频| 天美传媒精品一区二区| 欧美绝顶高潮抽搐喷水| 观看美女的网站| a级一级毛片免费在线观看| 日韩欧美三级三区| 亚洲人成网站高清观看| 亚洲激情在线av| 亚洲国产中文字幕在线视频| 午夜福利欧美成人| 国产v大片淫在线免费观看| 亚洲真实伦在线观看| av女优亚洲男人天堂| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| 国产精品99久久99久久久不卡| 91在线精品国自产拍蜜月 | 成人无遮挡网站| 成人国产综合亚洲| 久久6这里有精品| 香蕉久久夜色| 日本免费a在线| 韩国av一区二区三区四区| 久久久久久久午夜电影| 久久久久久久久中文| 国产成人福利小说| 免费人成在线观看视频色| 岛国视频午夜一区免费看| 尤物成人国产欧美一区二区三区| 婷婷六月久久综合丁香| 国产精品永久免费网站| 午夜免费观看网址| 真实男女啪啪啪动态图| 男女床上黄色一级片免费看| 看免费av毛片| 国产一区二区亚洲精品在线观看| 亚洲av中文字字幕乱码综合| 久久久久久大精品| 午夜免费男女啪啪视频观看 | 又黄又粗又硬又大视频| 叶爱在线成人免费视频播放| 久久精品亚洲精品国产色婷小说| 制服丝袜大香蕉在线| 欧美乱码精品一区二区三区| 欧美日韩乱码在线| 日本熟妇午夜| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 搡老妇女老女人老熟妇| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 国产精品99久久99久久久不卡| 亚洲黑人精品在线| av福利片在线观看| 桃红色精品国产亚洲av| 人妻夜夜爽99麻豆av| 久久午夜亚洲精品久久| 老熟妇仑乱视频hdxx| 欧美zozozo另类| 国产精品99久久久久久久久| av视频在线观看入口| 国产探花极品一区二区| 日韩av在线大香蕉| 日本成人三级电影网站| 久久久久国内视频| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 亚洲专区国产一区二区| 日本一二三区视频观看| 免费在线观看影片大全网站| 无人区码免费观看不卡| 午夜影院日韩av| 色视频www国产| 成年女人看的毛片在线观看| 麻豆国产av国片精品| 不卡一级毛片| 99精品久久久久人妻精品| 老熟妇仑乱视频hdxx| 亚洲成人久久爱视频| 伊人久久大香线蕉亚洲五| 亚洲成人久久爱视频| 又粗又爽又猛毛片免费看| 一级黄片播放器| 国产精品久久久久久亚洲av鲁大| 久久久国产成人免费| h日本视频在线播放| 99热只有精品国产| 国内少妇人妻偷人精品xxx网站| 色吧在线观看| 在线观看午夜福利视频| 免费av观看视频| 床上黄色一级片| 可以在线观看的亚洲视频| 国产高清三级在线| 成人三级黄色视频| 一本一本综合久久| 99精品在免费线老司机午夜| 成人av一区二区三区在线看| 99精品欧美一区二区三区四区| 久9热在线精品视频| 亚洲国产欧美网| 波多野结衣高清无吗| 女生性感内裤真人,穿戴方法视频| 国产精品亚洲av一区麻豆| 又黄又粗又硬又大视频| 亚洲美女黄片视频| 极品教师在线免费播放| 久久性视频一级片| 国产探花在线观看一区二区| 精品一区二区三区视频在线观看免费| 中文字幕人妻丝袜一区二区| 亚洲精品影视一区二区三区av| 国产熟女xx| 天堂√8在线中文| 国产97色在线日韩免费| 日韩欧美 国产精品| 欧美色视频一区免费| 99视频精品全部免费 在线| 99国产精品一区二区蜜桃av| 亚洲午夜理论影院| 哪里可以看免费的av片| 国产aⅴ精品一区二区三区波| 草草在线视频免费看| 国产精品久久电影中文字幕| 国产野战对白在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区色噜噜| 高清毛片免费观看视频网站| 老熟妇乱子伦视频在线观看| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩高清在线视频| 午夜a级毛片| 麻豆成人午夜福利视频| 亚洲电影在线观看av| 成人无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 18禁国产床啪视频网站| 亚洲精品乱码久久久v下载方式 | 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 久久久久久国产a免费观看| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 亚洲真实伦在线观看| 久久精品国产亚洲av香蕉五月| 免费观看的影片在线观看| 韩国av一区二区三区四区| 久久婷婷人人爽人人干人人爱| 欧美乱妇无乱码| 亚洲午夜理论影院| 久久精品国产亚洲av香蕉五月| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡| or卡值多少钱| 一级毛片高清免费大全| 少妇的逼好多水| 99热这里只有是精品50| 色综合亚洲欧美另类图片| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 中出人妻视频一区二区| 少妇人妻一区二区三区视频| 久久性视频一级片| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| 精品久久久久久成人av| 国产高潮美女av| 免费看美女性在线毛片视频| 国产精品久久电影中文字幕| 精品国产超薄肉色丝袜足j| 19禁男女啪啪无遮挡网站| 国产精品亚洲av一区麻豆| 黄色成人免费大全| 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 我的老师免费观看完整版| 免费观看人在逋| 精品一区二区三区av网在线观看| 少妇人妻一区二区三区视频| 精品无人区乱码1区二区| 日韩欧美 国产精品| 99精品久久久久人妻精品| 国产爱豆传媒在线观看| 亚洲成人久久爱视频| 好男人电影高清在线观看| 成人性生交大片免费视频hd| 中文亚洲av片在线观看爽| 天天躁日日操中文字幕| 99国产精品一区二区三区| 欧美性感艳星| 美女被艹到高潮喷水动态| 香蕉av资源在线| 国产精品 欧美亚洲| 亚洲成人中文字幕在线播放| 两人在一起打扑克的视频| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| 久久久久免费精品人妻一区二区| 国产私拍福利视频在线观看| 国产视频内射| av在线蜜桃| 欧美成人a在线观看| 亚洲午夜理论影院| 18禁国产床啪视频网站| 国产69精品久久久久777片| 老汉色∧v一级毛片| 91在线精品国自产拍蜜月 | 深爱激情五月婷婷| 欧美乱色亚洲激情| 成人午夜高清在线视频| 国产黄片美女视频| 性色avwww在线观看| 高清在线国产一区| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 18禁黄网站禁片免费观看直播| 色av中文字幕| av欧美777| 午夜久久久久精精品| 深夜精品福利| 日本三级黄在线观看| 国产蜜桃级精品一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲av成人不卡在线观看播放网| 午夜免费激情av| 亚洲色图av天堂| 免费看美女性在线毛片视频| 国产精品电影一区二区三区| 草草在线视频免费看| 成熟少妇高潮喷水视频| av天堂在线播放| 好男人在线观看高清免费视频| 精品不卡国产一区二区三区| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 男女床上黄色一级片免费看| 久久伊人香网站| 精华霜和精华液先用哪个| 无人区码免费观看不卡| 一进一出抽搐gif免费好疼| 亚洲欧美日韩高清专用| 嫩草影视91久久| 欧美成人免费av一区二区三区| 国产伦精品一区二区三区视频9 | 麻豆久久精品国产亚洲av| 在线观看免费午夜福利视频| 国内少妇人妻偷人精品xxx网站| 成人性生交大片免费视频hd| 国内少妇人妻偷人精品xxx网站| 久9热在线精品视频| 12—13女人毛片做爰片一| 禁无遮挡网站| 久久国产精品影院| 欧美在线一区亚洲| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站 | 国产成人aa在线观看| 精品国产超薄肉色丝袜足j| 国模一区二区三区四区视频| 动漫黄色视频在线观看| bbb黄色大片| 色在线成人网| 欧美中文综合在线视频| 免费电影在线观看免费观看| 成人欧美大片| 在线观看美女被高潮喷水网站 | 成人永久免费在线观看视频| 色播亚洲综合网| 久久国产精品影院| 给我免费播放毛片高清在线观看| 国产极品精品免费视频能看的| 男插女下体视频免费在线播放| 在线天堂最新版资源| 国产精品久久视频播放| 精品不卡国产一区二区三区| 一边摸一边抽搐一进一小说| av福利片在线观看| 午夜福利视频1000在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品久久视频播放| 麻豆成人av在线观看| 在线观看午夜福利视频| 女同久久另类99精品国产91| 午夜免费激情av| 亚洲中文字幕一区二区三区有码在线看| 18禁裸乳无遮挡免费网站照片| 成人特级黄色片久久久久久久| 青草久久国产| 亚洲真实伦在线观看| 91麻豆av在线| 欧美午夜高清在线| 99久久精品一区二区三区| 两个人视频免费观看高清| 亚洲 欧美 日韩 在线 免费| 亚洲欧美激情综合另类| 久久精品国产亚洲av香蕉五月| 国产精品野战在线观看| 日本免费a在线| 亚洲中文字幕日韩| 老司机午夜十八禁免费视频| 久久这里只有精品中国| 丁香六月欧美| av视频在线观看入口| 亚洲专区国产一区二区| 久久久色成人| 老熟妇仑乱视频hdxx| 两人在一起打扑克的视频| 18禁黄网站禁片免费观看直播| 蜜桃久久精品国产亚洲av| 亚洲乱码一区二区免费版| 成熟少妇高潮喷水视频| 热99在线观看视频|