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

    The Change Features of the West Boundary Bifurcation Line of the North Equatorial Current in the Pacific

    2015-04-01 01:57:18GUOJunruLIUYulongSONGJunBAOXianwenLIYanCHENShaoyangandYANGJinkun
    Journal of Ocean University of China 2015年6期

    GUO Junru, LIU Yulong, SONG Jun, 5), *, BAO Xianwen, LI Yan, CHEN Shaoyang, and YANG Jinkun

    ?

    The Change Features of the West Boundary Bifurcation Line of the North Equatorial Current in the Pacific

    GUO Junru1), 2), 4), LIU Yulong3), SONG Jun3), 5), *, BAO Xianwen1), 5), LI Yan3), CHEN Shaoyang3), and YANG Jinkun3)

    1),,266100,2),100194,3),300171,4),,,, ACT2602,5),,266100,

    The equatorial Current in the North Pacific (NEC) is an upper layer westward ocean current, which flows to the west boundary of the ocean, east of the Philippines, and bifurcates into the northerly Kuroshio and the main body of the southerly Mindanao current. Thus, NEC is both the south branch of the Subtropical Circulation and the north branch of the Tropical Circulation. The junction of the two branches extends to the west boundary to connect the bifurcation points forming the bifurcation line. The position of the North Pacific Equatorial Current bifurcation line of the surface determines the exchange between and the distribution of subtropical and tropical circulations, thus affecting the local or global climate. A new identification method to track the line and the bifurcation channel was used in this study, focusing on the climatological characteristics of the western boundary of the North Equatorial Current bifurcation line. The long-term average NEC west boundary bifurcation line shifts northwards with depth. In terms of seasonal variation, the average position of the western boundary of the bifurcation line is southernmost in June and northernmost in December, while in terms of interannual variation, from spring to winter in the years when ENSO is developing, the position of the west boundary bifurcation line of NEC is relatively to the north (south) in EI Ni?o (La Ni?a) years as compared to normal years.

    North Equatorial Current; Pacific Ocean; bifurcation line; climate change; ENSO

    1 Introduction

    NEC system connects the subtropical and tropical cy- clones in the Pacific and makes important contributions to the heat budget in the western Pacific warm pool (Qu, 1997), while the warm pool is a key region that affects the circulation of ENSO (Webster and Lukas, 1992). The western boundary current systems of NEC directly enters the global thermohaline circulation through the Indone- sian throughflow (Gordon, 1986) and plays a decisive role in the meridional transport of heat, salt and mass in the world oceans (Qu, 1998). The western Pacific region is an intersection for the middle and high latitude water masses (Fine, 1994), from which the concept of subtropical cell is derived: The net heating from the low latitude ocean is transported northwards from the surface to the NEC and then via Kuroshio to the higher latitudes; while the cold subsurface/deep water water flows southwards in the thermocline and then to the trop- ics through different passageways after being bifurcated by the NEC (Lu, 1998; Johnson and Mcphaden, 1999); and lastly turns upwards at the equator, constitute- ing the sallow meridional overturning circulation (Mc- Creary and Lu, 1994; Liu, 1994). The NEC links the sub- tropical and tropical regions, and is the boundary between subtropical and tropical circulations. At the western bound- ary, NEC bifurcates into two branches: NEC north branch passageway (flowing to Kuroshio) and NEC south branch passageway (flowing to the Mindanao current) (Xu and M-Rizzoli, 2013). There are different bifurcation points of NEC at different depths, the line connecting these bi- furcation points is the bifurcation line which determines the distribution of mass, heat and salt between the sub- tropical and tropical circulations. The interannual and decadal changes of the NEC bifurcation line can affect the tropical and subtropical heat distribution and thermal structure change, thus affecting the local or global climate (Song, 2011). Therefore, the change of the NEC bifurcation line plays a very important role in our re- search on the global climate change. Previous studies mostly focus on the change of the western boundary bifurcation points and less research is made on the change pattern of the bifurcation points, but this charac- teristic of bifurcation line can more directly and accu- rately identify the change of circulation systems (Wang and Liu, 2000).

    The position of the bifurcation point of NEC plays an important role in the Climate Change, so the research on the latitude of the NEC bifurcation point began as early as in the 1970s. According to the hydrological observations in 1934–1968, Nitani (1972) held that the surface of the NEC bifurcation point is about between 11?N and 14.5?N and the position of the bifurcation point moves north- wards with the increase of depth. Through analysis of the distribution of water masses, Toole(1988) estimated that the latitude of bifurcation was about 12?N. From the hydrological data from the Sino-American Joint Survey between September, 1987 and April, 1988, Toole(1990) defined the latitude where the transport flow func- tion is zero and concluded that the bifurcation latitude was near 13?N. Employing the model results, Qiu and Lukas (1996) defined the latitude where the average me- ridional speed is zero within the zone of two degrees of latitude along the western boundary and concluded that the bifurcation latitude was close to 13?N. According to the results of model calculation, Metzger and Hurlburt (1996) concluded that the average position of NEC moved in 15.6?N±0.6?N, shifting northwards with depth and reaching 18?N at the depth of 700m. On the basis of the climatic data of temperature and salinity, Qu(1999) defined the place where the average meridional speed is zero within two degrees of longitude from the coast as the position of bifurcation and concluded that near the surface the bifurcation line was at 13.5?N, at 500m in 18?N and the vertical average position between 0–500m was in about 15?N. By using the WOA (1994) historical observations, Qu(1999) concluded that in the upper thermocline, the NEC bifurcation point was at 15?N, moved northwards with depth and on the 27.2kgm?3isopycnic surface (about 800m), reached 20?N. By analysing the historical observations (WOA98), Qu and Lukas (2002) defined the latitude where the merid- ional transport within 5 degrees of longitude from the western boundary is zero as the latitude of bifurcation point and concluded that near the surface within 100m it was at 14?N, shifted northwards with depth and reached 23?N at 800m. By using the data of world ocean tempera- ture, salinity and dissolved oxygen in 1998 of WOD, Qu and Lukas (2003) found that the average position of the bifurcation latitude near surface (within 100 m) is 14.2?N, and if the effect of Ekman transport is added, it moves to 13.3?N while the bifurcation latitude shifts northwards with depth and at the depth of 800m, it is located in 20?N. Based on the calculated results of the model, Yaremchuk and Qu (2004) defined the latitude of the mean merid- ional speed change symbol and it was derived from cal- culation that the NEC average position is 14.3?±0.7?N. By employing the computed results of the high resolution OGCM model, Kim(2004) defined the latitude where the average meridional speed is zero within two degrees of longitude away from the coast of the Philip- pines as the bifurcation latitude and obtained the annual average bifurcation latitude of 15.5?N, which shifts northwards with depth, and the bifurcation latitude moves from 14.3?N at the surface to 16.6?N at 500 m. By use of the WOCE drifting buoy data in 1987–1998, Li(2005) defined that the position of the bifurcation point is between the trajectories of the buoy west of 130?E turning southwards and northwards. And in the case of the trajec- tories of the buoy turning southwards and northwards that cross west of 130?E, the latitude of the bifurcation point is defined as the latitude where the buoy trajectories cross and the surface bifurcation latitude is estimated to be between 11?N and 14.7?N. Wang and Hu (2006) used the satellite altimeter data (October, 1992–December, 2004) and concluded that the bifurcation latitude at the NEC western boundary surface is located at 13.4?N. Zhang(2008) used the satellite altimeter data of high spatial and temporal resolution to calculate the geostrophic current and defined the average latitude where the velocity gradient quickly increases and the me- ridional velocity is zero when reaching the shore as the criterion by which to determine the latitude of bifurcation. And the annual average bifurcation latitude was estimated to be 13.4?N.

    It can be seen that scientists (Zhai and Hu, 2012) have done a lot of researches on the flow field structure and bifurcation position of the western boundary current system, but most scholars located the bifurcation point in terms of the meridional velocity of the flow field, and as there are many factors affecting NEC and the change is complex, the difference in the results obtained from choosing to average over different width ranges is greater (Kim, 2004). Furthermore, even if the position of bifurcation point is determined, it is still impossible to reflect very well the bifurcation situation of the whole current system and obtain more detailed bifurcation fea- tures of the flow field. Hang and Wang (2001) and Wang and Huang (2005) proposed the method of quantitatively describing the barotrophic sea channel by wind stress, and that of quantitatively describing the baroclinic channel in recent years, by which to research into and analyse the property of the channel transport change with latitude. However, they did not make further research on the fea- tures of the bifurcation line.

    This paper compared the north Equatorial Current in the Pacific and tried to overcome the above mentioned shortcoming by using a new computational method,, identifying the position and change of bifurcation at each layer by calculating the isopycnal line, describing the NEC bifurcation line and making systematic analysis of its form and change pattern, analysing the plane distribu- tion of the bifurcation line on the two-dimensional plane, stratifying the data to obtain the vertical distribution of bifurcation lines and lastly deriving the climatic, seasonal and interannual changes of the bifurcation line though statistical analysis.

    2 Data and Methods

    2.1 Data Description

    This paper used the SODA (Simple Ocean Data As- similation) database version 2.0.2 (Carton and Giese, 2005), which is based on the GFDL (Geophysical Fluid Dynamics Laboratory) reanalysis data of the Parallel Ocean Program physics marine numerical model and data assimilation technology.

    This model includes the terrigenous freshwater flux with vertical resolution of 10m near the surface. The model’s wind stress forcing field is constructed by the COADS (Comprehensive Ocean-Atmosphere Data Set) and NECP (National Centers for Environmental Predic- tion) data (after 1992), of which the long-term trend of wind stress data not in harmony with the sea surface pressure trend observed has been eliminated. For the spe- cific description of these datasets involved in the research, refer to Carton and Giese (2008). The SODA data prod- ucts covers the 50 years from 1958 to 2007, with the global longitude range from 0.25? to 360.25? (from west to east) and the latitudinal range from ?75.25? to 89.25? (from south to north); the data have horizontal resolution 0.5?×0.5? and are vertically divided into 40 layers. The horizontal and vertical resolutions are high enough for studying the bifurcation of NEC.

    2.2 Method of Calculation and Data Processing

    The sea area selected in this study is in the North Pa- cific (100.25?E–120.25?W, 2.25?S–57.25?N). The paper used the data on the equipotential density surface layer by layer for analyzing the path. The specific method of cal- culation is as follows:

    The first step is to calculate the grid thermohaline data to get the potential temperature and density on thecoor- dinate (Except otherwise specified, the density or poten- tial density mentioned here after also represents the po- tential density). The layer-by-layer potential density av- eraged over the region is found and then the potential density layer is divided in equal intervals. In choosing the potential density layer, if it is too small, the low density water can only be reserved in the equatorial part of the west Pacific, which cannot display the path of water par- ticles; if it is too large, the bifurcation of western bound- ary trends to be blurred and disappears and so cannot be distinguished.

    The second step is to carry out the coordinate trans- formation, turning thecoordinate of raw data into the isopycnal surface coordinate for our analysis. First, the multiple density stratification is used to carry out the multi-layer division in the density stratification range and the algorithm of interpolation is employed to calculate this density range, and the meridional velocity, zonal ve- locity, temperature and salinity in isopycnal layer. Then, the isopycnal layer thickness is chosen according to the weight to represent the range. For the expanded and den- sified density layer, filter value taking is carried out ac- cording to the weight (the difference in the distance to the calculation point), and again the calculation is done to the required isopycnal surface, thus obtaining the variable required for the calculation under the density coordinate.

    The third step is to carry out the evaluation of integrals for the path with the nonlinear, high-order, one-step nu- merical algorithm of Runge-Kutta. The Runge-Kutta method is well-known as a multilevel algorithm, which is different from the integration method of ordinary differ- ential equations such as Adams-Bashorth, Adams-Monl- ton,(multiple methods) and it does not need to recal- culate the unknown function value of each order in the differential equation. It requires us to calculate the values offor a simple correction. The basic idea is to use the value of the intermediate step to replace the higher de- rivative.

    This paper uses for the first time the isopycnal surface stream line tracing and judging algorithm (Lagrangian particle tracking technique is used in this study), giving the accurate trace points of the NEC northern branch channel and southern branch channel and its bifurcation line on the calculated level by starting from the grid ve- locity field and tracking the stream lines within the chan- nels at different density layers. If the same calculation is done for the density layers, we can get the whole three-dimensional NEC western boundary current system. So this paper gives for the first time the precisely quanti- fied structure of the NEC western boundary channel and its bifurcation condition.

    3 Climatic Characteristics of the Pacific NEC Western Boundary Bifurcation Line

    This paper does not adopt the more limited method of bifurcation point, rather the integration path method to obtain the bifurcation line in the research. Fig.1 gives the long-term averaged stratified bifurcation line of the NEC western boundary. In the layer-by-layer change, with the increase of depth, the bifurcation line shifts northwards, which is consistent with the previous findings that in the sense of annual average the latitude of the western boundary bifurcation point moves northwards from the surface to deep layers. For the bifurcation line used in this paper, its distribution form of being high in the west and low in the east may be seen and the layer-by-layer bifur- cation lines tend to be parallel, do not overlap in the ver- tical direction and density changes little in the deeper layers, so the degree of separation between bifurcation lines increases on the deep-layer isopycnal surface. The latitudinal difference of a simple bifurcation line between the east and west ends of 130?E–170?E is about 1.5?.

    For the poleward shift of the NEC western boundary bifurcation line with depth, we approach this phenomenon from the oceanic circulation basis. In the vorticity equa- tion, the motion reaches a constant state and at the same time the exogenous process can be neglected (large-scale motion):

    Fig.1 Long-term averaged NEC western boundary stratification (23.4–26.4kgm?3) bifurcation lines.

    whereωis absolute vorticity;is horizontal flow veloc- ity;is density,is pressure.

    The relative vorticity of the large-scale motion is far smaller than the planetary vorticity:

    whereis the Coriolis parameter.

    The component form of baroclinic fluid may be ex- pressed as:

    hereis zonal coordinate with east positive;is merid- ional coordinate with north positive;is vertical coordi- nate with down positive.

    This is the thermogenic wind relationship, which con- structs the relationship between the variation of vertical flow velocity and the horizontal density (temperature), which is a very important relationship between flow ve- locity and density (temperature) in the ocean. Further- more, it is derived from the magnitude relation:

    In the meantime, for the vertical flow velocity shear and horizontal density gradients, we have:

    . (3-5)

    As the gradient of salinity is smaller, when the density (temperature) increases polewards, the increase of the current with depth quickly produces the flow westwards, as a result of which, the centre of the subtropical gyre (anticyclonic) moves polewards with the increase of depth (Liu, 2011). Therefore, we get the image of the shift polewards of the western boundary bifurcation line with depth due to the effect of the subtropical gyre.

    Figs.2 and 3 give the long-term averaged distribution of the NEC western boundary bifurcation lines in the sec- tion of 130?E under thecoordinate and potential density coordinate. Under the depth coordinate it may be seen that the depth of the first layer density 23.4kgm?3 in the whole section of 130?E over a span of 8?N–25?N between south and north is between 75m and 125m; on the depth coordinate are contained the information on the equal potential density line, isotherm and zonal velocity field. On the section of 130?E near the western boundary, the isotherm and the isopycnal line are basically parallel. At the same time, the upper temperature gradient in the low latitude area is greater than that in the high latitude area and the isopycnal line is far shallower. For instance, the isopycnal line of 26.4kgm?3is 210m at about 8?N, while at about 25?N it amounts to 555m. It can be seen from the figure that the western boundary bifurcation latitude in this section is not at the extremum centre of zonal veloc- ity, but is to the north a little. On the section of depth co- ordinate, the depth from the first layer to the last layer totals 500m in depth; the bifurcation line latitude varies for about 4?, with each 100m deepening the bifurcation latitude moves northwards for about 1 degree and the northward distance in the upper layer is slightly larger than that in the lower layer. In the density coordinates of Fig.3, the bifurcation line points shift in the density layer. As the upper layer is relatively thinner, the variation of bifurcation line points in the density coordinates is not obvious, while in the lower layer there is a distinct northward shift. This has been verified in Fig.2 in which in area of bifurcation line points, it may be clear that it is located at the largest bended location,, near the low-lying (deeper) place of the isopycnal surface.

    Fig.2 Long-term averaged 130?E sectional (depth coordinated) drawing, in which there are the NEC western boundary bifurcation lines (grey squares), isopycnal lines (solid lines) and equal zonal velocity lines (dotted lines).

    Fig.3 Long-term averaged 130?E sectional (depth coordinates) drawing, in which there are the NEC western boundary bifurcation lines (grey squares), isobaths lines (solid lines) and equal zonal transport lines (dotted lines).

    4 Seasonal Variation Features of the NEC Western Boundary Bifurcation Line

    Fig.4 gives the seasonal variations of the fine verti- cally-averaged NEC western boundary bifurcation line (23.4kgm?3–26.4kgm?3), derived on the basis of the channel algorithm. In winter (November, December and January), the bifurcation line is relatively straight and the difference in the latitude of bifurcation is not obvious in the east and west directions, which is related to the dis- tribution of NEC flow pattern in the winter, being about 15.7?N on the average. The bifurcation line is located at the higher latitudes throughout the year, being northern- most in November. Starting from February it obviously begins to move southwards gradually, showing inclining state of being high in the east and lower in the west. In June the westernmost end of the bifurcation line moves to the southernmost position in the year, in about 14.3?N on the average. This coincides with the findings of Qu and Lukas (2002) by use of the data on the west Pacific tem- perature, salinity, DO concentration,The research makes an integral calculation of the dynamic height and its corresponding circulation from surface to 1000 m, and in the light of the Ekman transport, concludes that the NEC bifurcation point begins to move southwards in January, reaches the southernmost end and then moves northwards, reaching the northernmost end in December. The time difference is mainly due to the different integra- tion methods used. This paper carries out the vertical in- tegration according to the density coordinate. And the calculated density layers are mostly within 100–600m, not including all the intervals of integration in 0–1000m. After June, the bifurcation line begins to shift northwards, the east-to-west dip is more obvious and the difference in latitude between the east and west ends of the line in au- tumn reaches the highest throughout the year, up to 2.8?, while such difference in winter when the bifurcation line is relatively straight is only 0.9?, and with the lapse of time, the bifurcation line moves northwards and the dip angle decreases gradually.

    Fig.4 Seasonal variations of the vertically averaged NEC western boundary bifurcation line (23.4–26.4kgm?3).

    Fig.5 is the seasonal variations of the NEC western boundary bifurcation lines layer by layer. In each sheet, there are 16 curves, representing the bifurcation lines in the 16 density layers. The distribution of the bifurcation lines on all the isopycnal surfaces in the whole area is relatively regular, and seasonally the latitude of bifurcation line moves northwards with depth and basically is parallel with the bifurcation line of the upper density layer. In winter, the bifurcation line is relatively higher in latitude and straight in shape while in summer and autumn it in- clines and at the west end of the line the wind direction changes and the wind stress magnitude is lower as com- pared with that in spring and winter. The location repre- sented by the first layer bifurcation line is not different from that indicated by the vertically averaged seasonal change,, it is moving northwards in winter and southwards in summer, being straight in winter and in- clining in autumn. The change in shape is displayed more obviously in the layer-by-layer bifurcation line. Among others, there exists a bifurcation line ridge layer by layer (the line bulges northwards), and the deviation of the ridge from the main axis of the bifurcation line may reach about 1.2?; it belongs to an obvious abnormal fluctuation. Such phenomenon begins to appear and develops in March and April; it speeds up in July, the intensity being lower at about 168?E, and in August propagates west- wards to 166?E, with a tendency to increase. In Septem- ber, the ridge moves to 163?E, further strengthening. By October, the ridge with the largest intensity has reached the location of 158?E; in November the follow-up fluc- tuation of October forms another ridge, which propagates westwards; and the former ridge moves westwards to about 155?E to form a double ridge with its intensity also dropping significantly. In December, the ridge propagates to the western boundary, which results in the up dip of the western boundary bifurcation point; and then it basically weakens and disappears. Figs.3–9 show the seasonal change in the NEC western boundary bifurcation line’s vertical dipping. A simple point-to-point method of inter- polation is used to treat the degree of diffusion of the layer-to-layer bifurcation line from bottom to surface along the longitudinal direction. The phenomenon with the maximum diffusion magnitude of the ridge in the whole calculated layer exceeding 4.4? appears in the spring and summer months of May and June, while the minimum magnitude of diffusion smaller than 2.1? occurs in the autumn months of September and October. It can be seen from the time evolution figure that the diffusion amplitude of the bifurcation line as a whole in January and February is smaller, the degree of diffusion is smaller, and the degree of diffusion in the western part of the line is even weaker than that in the eastern part. In March, the degree of diffusion of the bifurcation line in the vertical direction begins to increase and the amplitude fluctuation grows. By June, the vertical amplitude reaches the maxi- mum, generally exceeding 3.6? from west to east, which is the largest throughout the year. In July, the diffusion amplitude begins to drop and upon entering the autumn months of August, September and October, it contracts to the minimum level and part of the variation amplitude is as low as below 2.4?, which is the smallest throughout the year. Starting from November, the vertical variation am- plitude of the NEC western boundary bifurcation line again grows gradually, so we can say that the diffusion amplitude of the bifurcation line varies regularly with season.

    Fig.5 Seasonal variations of the layer-by-layer NEC western boundary bifurcation lines.

    Fig.7 gives the latitudes of the Pacific NEC western boundary bifurcation line and the depth values corre- sponding to it. In terms of the movement of the western boundary bifurcation line, as far as the 140m and 160m depth lines are concerned, the 160m depth line is located near the surface in June, which shows the depth of the bifurcation line is the largest in June. And in December, the 140m depth line is located far away from the surface, showing that the bifurcation line in December is the shal- lowest. Fluctuations are constantly seen on the bifurcation line to propagate from west to east. Starting from May there are two obvious wave crests located at the longi- tudes of 145?E and 170?E from top to bottom, constantly propagating westwards. By October the wave crest propa- gating from 170?E has developed most strongly while the 145?E wave crest has propagated up to the western boundary, thus causing the bifurcation line within 10 de- grees of longitude of the western boundary to swing from the southwest-northeast direction in May to the north- west-southeast direction in December and lowest latitu- dinal position to shift from the western boundary to about 140?E with a longitudinal difference of over 10 degrees.

    From the layer-by-layer angle with time as the variation axis, it may be clearly seen that the latitudinal change propagates from west to east and that the depth is the largest at the lowest latitude. Another phenomenon is that the lowest latitude of the western boundary bifurcation occurs around June while the latitude at 145?E which is the inflection point of the high value latitude propagation occurs around October, but at the bottom,, 26kgm?3, the lowest signal of the western boundary latitude occurs around July, the highest latitude is postponed to occur around November and possibly the signal that first changes propagates from the upper layer to the lower layer.

    5 Response of the Location of the NEC Western Boundary Bifurcation Line to the ENSO Event

    El Ni?o and La Ni?a years are chosen according to the Ni?o 3.4 exponent. The Ni?o 3.4 abnormal temperature maximum occurs in winter. The spring in the year when ENSO develops is usually chosen as the initial developing season and a multi-seasonal bifurcation line comparison is made with season as unit. Analysis is made according to four seasons of spring (March–May), summer (June–August), autumn (September–November) and winter (De- cember–February), which may be further divided into several time periods of the initial spring (March–May), developing summer (June–August), developing summer (June–August), developing autumn (September–Novem- ber), peak winter (December–February) and weakening spring (March–May) on the basis of the process of de- velopment-peak-weakening of ENSO and the variations of the bifurcation line are obtained from calculation as follows:

    Figs.8–15 give separately the NEC western boundary bifurcation lines in the spring (initial), summer (develop- ing), autumn (developing), and winter (peak) in the year when ENSO occurs and those in the next year’s spring (weakening), summer (weakening), autumn (weakening), and in the figures, the three lines indicate the El Ni?o year (Red Line), La Ni?a year (Blue Line), Long-Term Average (Black Line). It is usually believed that the spring of the year when ENSO occurs is the period when ENSO begins to develop, and in the spring of developing period, the upper bifurcation line first appears in the east to north in the form that it is higher (lower) of latitude in the El Ni?o (La Ni?a) year relative to normal years. With the change of seasons, up to the summer of the year when ENSO occurs, the three bifurcation line are completely separated from each other, the highest (lowest) latitude in the El Ni?o (La Ni?a) year. The bottom flow pattern is confused in the La Ni?a year and the situation in autumn is similar to that in summer, which will not be discussed in detail. The winter bifurcation lines in the year when ENSO occurs come back to the average state, the separa- tion degree at the east ends of the three bifurcation lines is reduced, the bifurcation lines begin to draw near, but their west ends still have a higher degree of separation, which is in consistent with the developing period.

    Fig.7 NEC western boundary bifurcating line and the variations at its corresponding depths (the shadow part is the latitude values at the two ends of the line).

    Fig.8 NEC western boundary bifurcation line in the initial spring: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    In the spring of the second year when ENSO occurs, the signed of the bifurcation line in the east by south (by north) further increases, propagating westwards to around 150?N, higher (lower) in the west in the El Ni?o year (La Ni?a year),which the bifurcation line in average years is located in the middle. The main intersection point of the summer bifurcation lines in the second year propagates up to the western boundary. In other parts except the western boundary, the El Ni?o (La Ni?a) year’s lower (higher) average lines are in the middle, and in the El Ni?o year there appear obvious fluctuation ridges of bi- furcation lines, the distribution patterns being basically constant in the upper and lower layers. In the second year’s autumn following the year when ENSO occurs, the bifurcation line transforms into the state where the El Ni?o (La Ni?a) year’s lower (higher) average line is in the middle, and in the meanwhile, the fluctuating peak formed in summer continues to propagate west wards, and the fluctuation in the El Ni?o year is stronger than that in the La Ni?a year. In the second year’s winter fol- lowing the year when ENSO occurs, the relative average value of the ENSO occurs, the relative average value of the NEC western boundary bifurcation line falls as a whole and in the following season returns to normal. By then, this ENSO event has ended.

    Fig.9 NEC western boundary bifurcation line in the developing summer: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.10 NEC western boundary bifurcation line in the developing autumn: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.11 NEC western boundary bifurcation line in the peak winter: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.12 NEC western boundary bifurcation line in the weakening spring: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.13 NEC western boundary bifurcation line in the weakening summer: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.14 NEC western boundary bifurcation line in the weakening autumn: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    Fig.15 NEC western boundary bifurcation line in the weakening winter: El Ni?o year (red line), La Ni?a year (blue line), long-term average (black line).

    6 Conclusions

    By use of the SODA data, this paper studies the cli- matic, interannual and seasonal variations of the NEC western boundary bifurcation line. It gives a visual dis- play of the channel bifurcation line proposed by this paper, makes an analysis of the position and its physical quantity change and obtains the following results with regularity:

    1) The long-term averaged NEC western boundary bi- furcation line shifts northwards with depth, and its main mechanism could be explained by the thermal wind rela- tionship in the ocean. The western boundary point of the bifurcation line is from 13.8?N of the 23.4kgm?3isopy- cnal surface to 17.7?N of the 26.4kgm?3isopycnal sur- face. It may be seen from the sections in different longi- tudes: with each 100 meters deepening, the latitudinal bifurcation moves northwards for about 1 degree and the northward shifting distance in the upper layer is slightly larger than that in the lower layer. The bifurcation line is located north of the centre with the maximum meridional velocity and near the low-lying (deeper) place of the isopycnal surface.

    2) The seasonal variation of the vertically averaged NEC western boundary bifurcation line is characterized by south-north moving and east-west dipping on the whole and the western boundary bifurcation point is a part of the whole movement of the bifurcation line. In terms of the zonal position of bifurcation line, it is north- ernmost in December with an average latitude of 15.7?N; In June, it is southernmost with an average latitude of about 14.3?N. The difference in the winter bifurcation latitude in the east and west directions is not obvious and the bifurcation line is relatively straight; The spring bi- furcation line begins to show a dipping state of being high in the east and low in the west, and the autumn latitudinal difference in the east and west directions is as large as over 2.8?.

    3) The seasonal variation of the NEC western boundary bifurcation line has a significant feature of the fluctuation propagating westwards, indicating that the bifurcation point is obviously modulated by the westward propagate- ing fluctuation (Rossby wave). Starting from July, the ridge occurring at 168?E (The bifurcation line bulges northwards) propagates westwards and is strengthened; it moves westwards in October to the place west of 158?E, with the amplitude reaching the maximum; In November it continues to propagate westwards but its intensity is obviously weakened; in December it propagates to the western boundary where it is basically weakened and disappears. The bifurcation lines of deeper layers in No- vember and December have a shape of three ridges and two troughs, but at the surface this distribution of forms is not distinct. In terms of the south-north movement of the NEC western boundary bifurcation line, the signal of latitude change in the lower layer lags behind that in the upper layer and the westward propagation of the bifurca- tion line ridge may cause the western boundary bifurca- tion point to swing between south and north. The bifurca- tion line in each layer has a northward jump by a large margin (The average latitude of bifurcation line in each layer has a northward jump of over 1.6?), which in time is in agreement with the first northward jump of the atmos- pheric subtropic high pressure system. In March the bi- furcation lines in all layers all withdraw to the south. Therefore the bifurcation lines may be divided into two clusters: being concentrated in the south in March–June and in the north in July–February. Corresponding to this, the average layer thickness of layers along the bifurcation line is also relatively small in March–June and larger in other months, which should be the result of vorticity conservation.

    4) The interannual variation of the NEC western boundary bifurcation line shows that in the developing and strengthening period of El Ni?o (La Ni?a), the bifur- cation line is to the north (to the south). Base on the proc- ess of development-peak-weakening of ENSO, several time periods are divided,, initial spring (March–May), developing summer (June–August), developing summer (June–August), developing autumn (September–Novem- ber), peak winter (December–February) and weakening spring (March–May). In the initial spring of El Ni?o, the bifurcation line is first in the east by north; in the devel- oping summer and autumn, the whole bifurcation line is to the north and reaches the maximum. In the peak winter, the bifurcation line restores the average state and the east part begins to go by south; in the weakening spring, the signal in the east by south further increases and propa- gates westwards. During the La Ni?a, the development of the bifurcation line is similar, but in opposite phase. Dur- ing the ENSO, there is an obvious feature of propagating westwards. The phenomenon that the maximum south-north shift of the western boundary bifurcation line oc- curs before the prime period of ENSO deserves further study in the future.

    Acknowledgements

    This work was supported by the National Natural Sci- ence Foundation of China (41206013, 41106004, 41376014, 41430963); Key Marine Science Foundation of the State Oceanic Administration of China for Young Scholar (2013203, 2012202, 2012223); POL Visiting Fellowship Program (Jun Song); the Public Science and Technology ResearchFundsProjectsofOcean(201205018,201005019); China Scholarship Council ([2008]3019, [2012]3013).

    Carton, J. A., and Giese, B. S., 2005. SODA: A reanalysis of ocean climate., 1: 1-30.

    Carton, J. A., and Giese, B. S., 2008. A reanalysis of ocean climate using simple ocean data assimilation (SODA)., 136: 2999-3017.

    Fine, R., Lukas, R., Bingham, F. M., Warner, M. J., and Gammon, R. H., 1994. The western equatorial Pacific is a water mass crossroads., 99: 25063-25080.

    Gordon, A. L., 1986. Interocean exchange of thermocline water., 91: 5037-5046.

    Huang, R. X., and Wang, Q., 2001. Interior communication from the subtropical to the tropical oceans., 31: 3538-3550.

    Johnson, G. C., and McPhaden, M. J., 1999. Interiod pycnocline flow from the subtropical to the equatorial Pacific Ocean., 29: 3073-3089.

    Kim, Y. Y., Qu, T., and Jensen, T, 2004. Seasonal and interannual variations of the North Equatorial Current bifurcation in a high-resolution OGCM., 109, C03040.

    Li, L. J., Liu, Q. Y., and Liu, W., 2005. Surface current speed andbifurcation of the North Equatorial Current in the Pacific Ocean., 35 (3): 370-374.

    Liu, Y. L., Wang, Q., Song, J., Zhu, X. D., Gong, X. Q., and Wu, F., 2011. Numerical study on the bifurcation of the North Equatorial Current., 10 (4): 305-313.

    Liu, Z., 1994. A simple model of the mass exchange between the subtropical and tropical ocean., 24: 1153-1165.

    Lu, P., McCreary, J. P., and Klinger, B. A., 1998. Meridional circulation cell and the source water of the Pacific equatorial undercurrent., 28:62-84.

    McCreary, J. P., and Lu, P., 1994. On the interaction between the subtropical and equatorial ocean circulation: The subtropical cell., 24: 466-497.

    Metzger, E. J., and Hurlburt, H. E., 1996. Coupled dynamics of the South China Sea, the Sulu Sea, and the Pacific Ocean., 101: 12331-12352.

    Nitani, H., 1972.. University of Washington Press, Washington, 129-163.

    Qiu, B., and Lukas, R., 1996. Seasonal and interannual variability of the North Equatorial Current, the Mindanao Current, and the Kuroshio along the Pacific western boundary., 101: 12315-12330.

    Qu, T. D., Meyers, J. G., and Godfrey, S., 1997. Upper ocean dynamics and its role in maintaining the annual mean western Pacific warm pool in a global GCM., 17: 711-724.

    Qu, T. D., Mitsudera, H., and Yamagata, T., 1998. On the western boundary currents in the Philippine Sea., 103 (4): 7537-7548.

    Qu, T. D., Mitsudera, H., and Yamagata, T., 1999. A climatology of the circulation and water mass distribution near the Phi- lippine coast., 29: 1488-1505.

    Qu, T., and Lukas, R., 2002. Depth distribution of the subtropical Gyre in the North Pacific., 58 (3): 525-529.

    Qu, T., and Lukas, R., 2003. The bifurcation of the North Equatorial Current in the Pacific., 33: 5-18.

    Song, J., Xue, H. J., Bao, X. W., and Wu, D. X., 2011. A spectral mixture model analysis of the Kuroshio variability and the water exchange between the Kuroshio and the East China Sea., 29 (2): 446-459.

    Toole, J. M., Zou, E., and Millard, R. C., 1988. On the circulation of the upper waters in the western equatorial Pacific Ocean., 35: 1451-1482.

    Toole, J., Millard, C., Wang, Z., and Pu, S., 1990. Observations of the Pacific North Equatorial Current bifurcation at the Philippine coast., 20: 307-318.

    Wang, Q., and Huang, R. X., 2005. Decadal variability of pycnocline flows from the subtropical to the equatorial Pacific., 35 (10): 1861-1875.

    Wang, Q., and Liu, Q. Y., 2000. The study on the Pacific Subtropical Cell., 1: 1-5.

    Wang, Q. Y., and Hu, D. X., 2006. Bifurcation of the North Equatorial Current derived from altimetry in the Pacific Ocean., 18 (5): 620-626.

    Webster, P. J., and Lukas, R., 1992. TOGA/COARE: The coupledocean-atmosphere response experiment., 73: 1377-1416

    Xu, D., and M-Rizzoli, P., 2013. Seasonal variation of the upper layer of the South China Sea and the Indonesian Seas: An ocean model study., 63: 103-130.

    Yaremchuk, M., and Qu, T. D., 2004. Seasonal variability of the large-scale currents near the coast of the Philippines., 34 (4): 844-855.

    Zhai, F. G., and Hu, D. X., 2012. Interannual variability of transport and bifurcation of the North Equatorial Current in the tropical North Pacific Ocean., 30 (1): 177-185.

    Zhang, X. D., Xiu, Y. R., Liu, J. F., Su, G., and Wang, Q. Y., 2008. The determination about the bifurcation latitude of the NEC and its variations., 25 (2): 33-41.

    (Edited by Xie Jun)

    DOI 10.1007/s11802-015-2594-0

    ISSN 1672-5182, 2015 14 (6): 957-968

    ? Ocean University of China, Science Press and Springer-Verlag Berlin Heidelberg 2015

    (February 6, 2014; revised March 12, 2014; accepted September 2, 2015)

    * Corresponding author. E-mail: thunder098@hotmail.com

    欧美绝顶高潮抽搐喷水| 欧美一级a爱片免费观看看 | 亚洲av片天天在线观看| 日日摸夜夜添夜夜添小说| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久久毛片| 每晚都被弄得嗷嗷叫到高潮| 99国产精品一区二区三区| 男男h啪啪无遮挡| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人av| 亚洲一区二区三区色噜噜| 日韩有码中文字幕| АⅤ资源中文在线天堂| 一本大道久久a久久精品| 怎么达到女性高潮| 欧美成人免费av一区二区三区| 美女 人体艺术 gogo| 极品教师在线免费播放| 国产主播在线观看一区二区| 他把我摸到了高潮在线观看| 少妇被粗大的猛进出69影院| 黄色视频,在线免费观看| 亚洲av成人av| 亚洲精品美女久久久久99蜜臀| 黄色女人牲交| 91老司机精品| 制服诱惑二区| а√天堂www在线а√下载| 美女免费视频网站| 久久久久久免费高清国产稀缺| 日本免费一区二区三区高清不卡| 国产av一区二区精品久久| 中文字幕熟女人妻在线| 亚洲成人久久爱视频| 午夜久久久久精精品| 国产免费男女视频| www国产在线视频色| 免费看日本二区| 久久久国产成人精品二区| 国产区一区二久久| 99久久精品热视频| 国产99白浆流出| 亚洲中文字幕日韩| 亚洲人成网站在线播放欧美日韩| 日韩精品中文字幕看吧| 欧美乱码精品一区二区三区| 日韩欧美在线二视频| 午夜精品在线福利| 又大又爽又粗| 18禁美女被吸乳视频| 亚洲av成人av| 久久精品综合一区二区三区| 成人三级黄色视频| 国产精品久久久久久亚洲av鲁大| 欧美不卡视频在线免费观看 | 一级毛片女人18水好多| 亚洲天堂国产精品一区在线| 中文字幕最新亚洲高清| 蜜桃久久精品国产亚洲av| 美女 人体艺术 gogo| 俺也久久电影网| 亚洲欧美精品综合久久99| 色综合婷婷激情| 国产亚洲av高清不卡| 国产av在哪里看| 亚洲成人国产一区在线观看| 淫秽高清视频在线观看| 亚洲精品国产精品久久久不卡| 亚洲国产精品999在线| 在线观看免费视频日本深夜| 99热这里只有精品一区 | 99精品在免费线老司机午夜| e午夜精品久久久久久久| 99在线视频只有这里精品首页| 亚洲人成电影免费在线| aaaaa片日本免费| 不卡av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 久久这里只有精品19| 欧美另类亚洲清纯唯美| www日本黄色视频网| 国产精品,欧美在线| 人成视频在线观看免费观看| 别揉我奶头~嗯~啊~动态视频| 国内久久婷婷六月综合欲色啪| 在线播放国产精品三级| x7x7x7水蜜桃| www国产在线视频色| 成人精品一区二区免费| 老司机靠b影院| 在线视频色国产色| 深夜精品福利| 成人三级黄色视频| 国产av不卡久久| 日韩大码丰满熟妇| 欧美一级a爱片免费观看看 | 午夜福利高清视频| 久久人人精品亚洲av| av欧美777| 成在线人永久免费视频| 中文字幕精品亚洲无线码一区| 亚洲国产中文字幕在线视频| 特级一级黄色大片| 亚洲国产精品999在线| 久久久精品大字幕| 99热这里只有是精品50| 国产在线精品亚洲第一网站| 美女扒开内裤让男人捅视频| 97超级碰碰碰精品色视频在线观看| 国产亚洲欧美在线一区二区| 伊人久久大香线蕉亚洲五| 99精品欧美一区二区三区四区| 1024视频免费在线观看| 香蕉国产在线看| 精品午夜福利视频在线观看一区| 正在播放国产对白刺激| 国产伦在线观看视频一区| 亚洲国产欧美网| 久久婷婷成人综合色麻豆| 国产91精品成人一区二区三区| 亚洲精品中文字幕一二三四区| 精品无人区乱码1区二区| xxx96com| 999久久久精品免费观看国产| www国产在线视频色| 91老司机精品| 一区二区三区高清视频在线| 国产午夜精品久久久久久| 亚洲色图 男人天堂 中文字幕| 久久伊人香网站| 成人欧美大片| 妹子高潮喷水视频| 男女床上黄色一级片免费看| 欧美色视频一区免费| 国产久久久一区二区三区| 亚洲中文日韩欧美视频| 99久久久亚洲精品蜜臀av| 久久久久久久久中文| 国产精品1区2区在线观看.| 欧美绝顶高潮抽搐喷水| 国产野战对白在线观看| 亚洲无线在线观看| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美精品综合一区二区三区| 18禁黄网站禁片免费观看直播| 90打野战视频偷拍视频| 美女 人体艺术 gogo| 老司机午夜福利在线观看视频| 久久久久久亚洲精品国产蜜桃av| 日韩有码中文字幕| 99久久国产精品久久久| 亚洲免费av在线视频| 日本三级黄在线观看| 在线十欧美十亚洲十日本专区| 2021天堂中文幕一二区在线观| 天堂√8在线中文| 一二三四社区在线视频社区8| 男女做爰动态图高潮gif福利片| 一个人免费在线观看电影 | 一个人观看的视频www高清免费观看 | 美女大奶头视频| 亚洲,欧美精品.| 精品欧美国产一区二区三| 最近最新中文字幕大全电影3| 欧美成人午夜精品| 亚洲精品色激情综合| 一本久久中文字幕| 欧美最黄视频在线播放免费| 久久婷婷成人综合色麻豆| 巨乳人妻的诱惑在线观看| 久久精品aⅴ一区二区三区四区| 搞女人的毛片| 精品一区二区三区av网在线观看| 伦理电影免费视频| bbb黄色大片| 天天添夜夜摸| 国产成人精品无人区| 欧美乱码精品一区二区三区| 精品久久久久久成人av| 中亚洲国语对白在线视频| 别揉我奶头~嗯~啊~动态视频| 男女之事视频高清在线观看| 性欧美人与动物交配| 精品久久久久久久久久久久久| 中文字幕久久专区| 男人舔奶头视频| 丝袜美腿诱惑在线| 热99re8久久精品国产| 亚洲中文日韩欧美视频| 18禁黄网站禁片免费观看直播| 欧美日韩乱码在线| 精品欧美一区二区三区在线| 国产乱人伦免费视频| 午夜a级毛片| 久久亚洲真实| 久久久久精品国产欧美久久久| 婷婷精品国产亚洲av| 欧美黑人精品巨大| tocl精华| 精品人妻1区二区| 国产精品久久久久久精品电影| 欧美成人一区二区免费高清观看 | 在线视频色国产色| 日本精品一区二区三区蜜桃| 国产一区二区在线av高清观看| 久久久水蜜桃国产精品网| 一本一本综合久久| 亚洲国产欧美网| 中文字幕人成人乱码亚洲影| 性色av乱码一区二区三区2| 中出人妻视频一区二区| 国产精品国产高清国产av| 全区人妻精品视频| 久久久久久国产a免费观看| 国产亚洲av嫩草精品影院| 亚洲国产精品久久男人天堂| 国产亚洲av高清不卡| 久久中文看片网| 香蕉国产在线看| 一级毛片高清免费大全| 国产精品乱码一区二三区的特点| 亚洲美女黄片视频| 欧美又色又爽又黄视频| 露出奶头的视频| 午夜福利成人在线免费观看| 欧美激情久久久久久爽电影| 亚洲av成人不卡在线观看播放网| 丝袜美腿诱惑在线| 亚洲欧美一区二区三区黑人| 狂野欧美激情性xxxx| 亚洲av片天天在线观看| 天天躁狠狠躁夜夜躁狠狠躁| av视频在线观看入口| xxx96com| 欧美一区二区精品小视频在线| 一个人免费在线观看的高清视频| 午夜福利成人在线免费观看| 亚洲一区二区三区色噜噜| 99国产综合亚洲精品| 制服丝袜大香蕉在线| 欧美成狂野欧美在线观看| 窝窝影院91人妻| 亚洲成人精品中文字幕电影| 国产成人影院久久av| 国产人伦9x9x在线观看| 999久久久国产精品视频| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| 久久久精品大字幕| 国产午夜精品论理片| 午夜福利18| 这个男人来自地球电影免费观看| 国产激情偷乱视频一区二区| 亚洲国产精品999在线| 国产又色又爽无遮挡免费看| 免费观看人在逋| 无人区码免费观看不卡| 在线视频色国产色| 精品免费久久久久久久清纯| 国产人伦9x9x在线观看| 久久久久久九九精品二区国产 | or卡值多少钱| 操出白浆在线播放| 国产激情偷乱视频一区二区| 嫩草影视91久久| 国产高清视频在线观看网站| 女人被狂操c到高潮| 亚洲18禁久久av| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲美女视频黄频| 亚洲成人久久爱视频| 日本熟妇午夜| 老汉色av国产亚洲站长工具| 久久香蕉国产精品| 美女 人体艺术 gogo| netflix在线观看网站| 亚洲一区中文字幕在线| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| 国产成人av教育| 50天的宝宝边吃奶边哭怎么回事| 国产成+人综合+亚洲专区| 麻豆久久精品国产亚洲av| 亚洲国产欧洲综合997久久,| 亚洲中文av在线| 日本一二三区视频观看| 亚洲最大成人中文| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 国产激情久久老熟女| 欧美一级毛片孕妇| 窝窝影院91人妻| 伊人久久大香线蕉亚洲五| 国产精品久久久久久精品电影| 亚洲人与动物交配视频| 不卡av一区二区三区| АⅤ资源中文在线天堂| 国产激情久久老熟女| 又黄又粗又硬又大视频| 国产精品一区二区三区四区久久| 两人在一起打扑克的视频| 久久久久免费精品人妻一区二区| 哪里可以看免费的av片| 首页视频小说图片口味搜索| 国产乱人伦免费视频| 香蕉av资源在线| 老汉色av国产亚洲站长工具| 久久伊人香网站| 又粗又爽又猛毛片免费看| 一a级毛片在线观看| 日本a在线网址| 非洲黑人性xxxx精品又粗又长| 国产精品免费视频内射| 极品教师在线免费播放| 亚洲成人久久爱视频| 国产激情偷乱视频一区二区| 欧美日韩瑟瑟在线播放| 精品人妻1区二区| 国产一区二区在线观看日韩 | 亚洲,欧美精品.| 丝袜人妻中文字幕| 欧美绝顶高潮抽搐喷水| 久久久久久久久久黄片| 国内精品久久久久久久电影| av免费在线观看网站| 19禁男女啪啪无遮挡网站| 高清毛片免费观看视频网站| 两个人看的免费小视频| 国产99久久九九免费精品| 欧美乱色亚洲激情| 好男人电影高清在线观看| 在线观看美女被高潮喷水网站 | 男男h啪啪无遮挡| 特级一级黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻1区二区| 国产亚洲精品久久久久5区| 俄罗斯特黄特色一大片| 最新在线观看一区二区三区| 亚洲国产精品合色在线| 日本成人三级电影网站| 久久久久亚洲av毛片大全| 欧美在线黄色| АⅤ资源中文在线天堂| 1024手机看黄色片| 久久久久国产精品人妻aⅴ院| 国产精品一及| 午夜老司机福利片| 亚洲精品在线美女| 久久久久久久精品吃奶| 亚洲国产精品合色在线| av超薄肉色丝袜交足视频| 免费无遮挡裸体视频| 一夜夜www| 麻豆久久精品国产亚洲av| 88av欧美| 久久国产精品人妻蜜桃| 两个人免费观看高清视频| 欧美3d第一页| 国产三级在线视频| 精品久久久久久久末码| 久久久国产成人精品二区| 变态另类成人亚洲欧美熟女| 99在线人妻在线中文字幕| 午夜成年电影在线免费观看| 最好的美女福利视频网| 99精品久久久久人妻精品| 久热爱精品视频在线9| 久久精品综合一区二区三区| 国产探花在线观看一区二区| 久久精品国产亚洲av香蕉五月| 99久久精品国产亚洲精品| 精品电影一区二区在线| 亚洲成av人片免费观看| 老司机深夜福利视频在线观看| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 久久这里只有精品中国| 亚洲精品国产一区二区精华液| 午夜两性在线视频| 亚洲精品久久国产高清桃花| 国产黄片美女视频| 一个人免费在线观看电影 | 国产69精品久久久久777片 | 久久久久久免费高清国产稀缺| 国产精品久久久久久人妻精品电影| 91在线观看av| www日本黄色视频网| 午夜老司机福利片| 我的老师免费观看完整版| 午夜福利在线观看吧| 99国产精品一区二区三区| 嫩草影院精品99| 亚洲av中文字字幕乱码综合| 久久久水蜜桃国产精品网| 又爽又黄无遮挡网站| 少妇的丰满在线观看| 91字幕亚洲| 国产欧美日韩一区二区精品| 999精品在线视频| a在线观看视频网站| 欧美性长视频在线观看| 午夜亚洲福利在线播放| 熟妇人妻久久中文字幕3abv| 淫秽高清视频在线观看| 久久 成人 亚洲| 哪里可以看免费的av片| √禁漫天堂资源中文www| 国产精品久久久av美女十八| 午夜激情福利司机影院| 男人舔女人的私密视频| 国产精品美女特级片免费视频播放器 | xxxwww97欧美| 亚洲av成人不卡在线观看播放网| 亚洲人成电影免费在线| 一进一出抽搐gif免费好疼| 狠狠狠狠99中文字幕| 一本大道久久a久久精品| av在线天堂中文字幕| 成人欧美大片| 哪里可以看免费的av片| 日日爽夜夜爽网站| 19禁男女啪啪无遮挡网站| 国产探花在线观看一区二区| 19禁男女啪啪无遮挡网站| 看免费av毛片| 久久精品影院6| 国产成人av激情在线播放| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3| 久久天堂一区二区三区四区| 久久人妻av系列| 精品乱码久久久久久99久播| 99热这里只有是精品50| 久久久久久久久久黄片| 日本一本二区三区精品| 欧美黄色淫秽网站| 动漫黄色视频在线观看| 香蕉国产在线看| 亚洲精品中文字幕在线视频| 成人手机av| 小说图片视频综合网站| 国产精品 国内视频| 变态另类成人亚洲欧美熟女| 亚洲五月婷婷丁香| 妹子高潮喷水视频| 后天国语完整版免费观看| 国产亚洲精品av在线| 亚洲精品粉嫩美女一区| 日日爽夜夜爽网站| 琪琪午夜伦伦电影理论片6080| 脱女人内裤的视频| 女生性感内裤真人,穿戴方法视频| 99热这里只有是精品50| 天堂动漫精品| 黄色a级毛片大全视频| 国产主播在线观看一区二区| www.熟女人妻精品国产| 精品乱码久久久久久99久播| 日韩大尺度精品在线看网址| 国产三级黄色录像| 中文字幕高清在线视频| 国产一区二区在线av高清观看| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 两性午夜刺激爽爽歪歪视频在线观看 | 国产欧美日韩一区二区精品| 亚洲专区中文字幕在线| 亚洲第一欧美日韩一区二区三区| 亚洲成人国产一区在线观看| 看黄色毛片网站| 黄频高清免费视频| 成人特级黄色片久久久久久久| 日本一区二区免费在线视频| 巨乳人妻的诱惑在线观看| 日韩欧美国产在线观看| 国产亚洲欧美在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 成人国语在线视频| 91麻豆精品激情在线观看国产| 一级毛片女人18水好多| 男女之事视频高清在线观看| 成人特级黄色片久久久久久久| 热99re8久久精品国产| 99re在线观看精品视频| 久久久精品大字幕| 精品国产乱码久久久久久男人| 久热爱精品视频在线9| 一级作爱视频免费观看| 亚洲成人中文字幕在线播放| 精品第一国产精品| 极品教师在线免费播放| 小说图片视频综合网站| 久久这里只有精品19| avwww免费| 国产探花在线观看一区二区| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸| x7x7x7水蜜桃| 国产又色又爽无遮挡免费看| 久久久久久久久免费视频了| 久久久久久久午夜电影| 亚洲色图av天堂| 午夜精品久久久久久毛片777| 男女午夜视频在线观看| 亚洲天堂国产精品一区在线| 亚洲av电影在线进入| 无限看片的www在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲 国产 在线| 99热这里只有精品一区 | 国产一区二区三区在线臀色熟女| 亚洲专区国产一区二区| 国产真实乱freesex| 午夜福利在线在线| 正在播放国产对白刺激| 久久久精品欧美日韩精品| av免费在线观看网站| 国产成人aa在线观看| 巨乳人妻的诱惑在线观看| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 欧美日韩黄片免| 久久精品国产99精品国产亚洲性色| 久久99热这里只有精品18| 波多野结衣巨乳人妻| 一个人观看的视频www高清免费观看 | 亚洲欧美日韩东京热| 一本久久中文字幕| 亚洲美女黄片视频| 可以免费在线观看a视频的电影网站| 国产免费av片在线观看野外av| 窝窝影院91人妻| 亚洲aⅴ乱码一区二区在线播放 | 国产精华一区二区三区| 日韩欧美精品v在线| 国产精品一区二区免费欧美| 国产精品av久久久久免费| 麻豆成人av在线观看| 亚洲欧美激情综合另类| 久久中文看片网| 国产亚洲欧美在线一区二区| 精品国内亚洲2022精品成人| 国产高清视频在线观看网站| 精品国产乱子伦一区二区三区| 国产精品九九99| 看片在线看免费视频| av有码第一页| 韩国av一区二区三区四区| 最近最新免费中文字幕在线| 久久欧美精品欧美久久欧美| 午夜a级毛片| 国产精品野战在线观看| 久久亚洲真实| 99国产极品粉嫩在线观看| 国产av一区在线观看免费| 九色成人免费人妻av| 岛国在线观看网站| bbb黄色大片| 少妇被粗大的猛进出69影院| 一级黄色大片毛片| 国产区一区二久久| 亚洲18禁久久av| 在线十欧美十亚洲十日本专区| 亚洲一卡2卡3卡4卡5卡精品中文| 99re在线观看精品视频| 波多野结衣巨乳人妻| 午夜福利在线观看吧| 免费看a级黄色片| 91老司机精品| 69av精品久久久久久| www.自偷自拍.com| 国产又黄又爽又无遮挡在线| 久久中文字幕人妻熟女| 夜夜看夜夜爽夜夜摸| 中出人妻视频一区二区| 亚洲国产欧美人成| 黑人巨大精品欧美一区二区mp4| x7x7x7水蜜桃| 美女高潮喷水抽搐中文字幕| 欧美日韩瑟瑟在线播放| 无人区码免费观看不卡| 国产91精品成人一区二区三区| 国产免费男女视频| 亚洲人成电影免费在线| 丰满的人妻完整版| 国产熟女xx| 1024手机看黄色片| 国产精品av视频在线免费观看| netflix在线观看网站| 最近最新免费中文字幕在线| 午夜影院日韩av| 国产99久久九九免费精品| 国产精品久久视频播放| 色在线成人网| 日日爽夜夜爽网站| 日韩高清综合在线| 精华霜和精华液先用哪个| 成人欧美大片| 丝袜人妻中文字幕| av国产免费在线观看| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 美女大奶头视频| 精品久久久久久久人妻蜜臀av| cao死你这个sao货| 999久久久国产精品视频|