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

    斯里蘭卡穹頂區(qū)的形成和演變機(jī)制分析

    2019-05-30 01:10:52王衛(wèi)強(qiáng)

    王衛(wèi)強(qiáng)

    摘要 通過(guò)SODA再分析資料和AVISO觀測(cè)資料研究了斯里蘭卡穹頂區(qū)(SLD)的遷移和消散機(jī)制.斯里蘭卡穹頂區(qū)是孟加拉灣西南部的一個(gè)氣旋渦旋,主要出現(xiàn)在西南季風(fēng)(5—9月)期間,與西南季風(fēng)海流侵入孟加拉灣同時(shí)存在.正風(fēng)應(yīng)力旋度引起的Ekman抽吸是形成SLD的主要原因.回歸分析結(jié)果表明SLD區(qū)域的風(fēng)應(yīng)力旋度與Ekman抽吸存在較強(qiáng)的正相關(guān)(r2=0.93,p>0.5).此外,結(jié)果表明SLD在發(fā)展過(guò)程中的移動(dòng)主要受正風(fēng)應(yīng)力旋度移動(dòng)的影響,SLD的消減與該正風(fēng)應(yīng)力旋度減弱和西傳的暖Rossby波有關(guān),而冷Rossby波的傳播有益于SLD的發(fā)展.在SLD消減時(shí)期,孟加拉灣渦旋(BBD)獨(dú)立發(fā)展并進(jìn)一步與SLD融合,回歸分析發(fā)現(xiàn)BBD區(qū)域的Ekman抽吸與當(dāng)?shù)仫L(fēng)應(yīng)力旋度的關(guān)系密切(r2=0.76,p>0.5),這表明了BBD在形成階段由局地的風(fēng)應(yīng)力主導(dǎo).9月之后,風(fēng)應(yīng)力旋度減弱,BBD和SLD開(kāi)始了合并過(guò)程.動(dòng)力方面,EKE分析顯示SLD衰退的同時(shí),BBD的EKE大幅增加;熱力方面,10—11月時(shí),由Ekman抽吸引起的SLD和BBD次表層冷水匯合,清晰地表明了二者之間的熱動(dòng)力學(xué)聯(lián)系.

    關(guān)鍵詞 斯里蘭卡穹頂區(qū);風(fēng)應(yīng)力旋度;Ekman抽吸;Rossby波;孟加拉灣渦旋

    中圖分類(lèi)號(hào) P444

    文獻(xiàn)標(biāo)志碼 A

    0 引言

    斯里蘭卡穹頂區(qū)(SLD)是孟加拉灣西南部的一個(gè)氣旋渦旋,主要出現(xiàn)在西南季風(fēng)(5—9月)期間,和西南季風(fēng)海流侵入孟加拉灣同時(shí)存在.前人研究通過(guò)分析氣候態(tài)海表溫度發(fā)現(xiàn),SLD位于斯里蘭卡東部,在5月生成,9月開(kāi)始衰減,盛期出現(xiàn)在7月.對(duì)SLD形成與衰退的機(jī)制研究表明,在西南季風(fēng)期間斯里蘭卡東部存在一個(gè)正的風(fēng)應(yīng)力旋度,從而引起局地的Ekman抽吸形成SLD;之后風(fēng)應(yīng)力旋度減弱,雖然正的風(fēng)應(yīng)力旋度一直持續(xù)到12月,但是西傳的暖Rossby波會(huì)引起穹頂區(qū)的下沉從而導(dǎo)致SLD的衰退.另外,在SLD衰減過(guò)程中,在斯里蘭卡東北部形成一個(gè)新的氣旋式渦旋,即孟加拉灣渦旋(BBD).有觀點(diǎn)認(rèn)為BBD是SLD氣旋式渦旋信號(hào)傳播導(dǎo)致,而其他的觀點(diǎn)認(rèn)為BBD是獨(dú)立發(fā)展形成的,并與SLD相互影響.本研究試圖通過(guò)分析簡(jiǎn)單海洋數(shù)據(jù)同化再分析2.0.2版本的20年月平均氣候態(tài)數(shù)據(jù),來(lái)描述SLD的遷移和消散機(jī)制.

    本文研究發(fā)現(xiàn),在西南季風(fēng)期間,斯里蘭卡東側(cè)存在較復(fù)雜的風(fēng)應(yīng)力旋度分布,首先是位于斯里蘭卡東側(cè)的正的風(fēng)應(yīng)力旋度(SLDPWC),與SLDPWC相關(guān)的Ekman抽吸機(jī)制對(duì)SLD的形成給出了合理的解釋,回歸分析也表明了SLDPWC和SLD的上升流之間存在著很強(qiáng)的正相關(guān)關(guān)系(r2=0.93,p >0.5).該氣旋式環(huán)流與正風(fēng)應(yīng)力旋度一起在5月形成,6—8月逐漸增強(qiáng)并向東北方向傳播,清楚表明了SLD在形成過(guò)程中和正風(fēng)應(yīng)力旋度的一致性.關(guān)于SLD的衰退機(jī)制,由Hovmoller 分析顯示,東邊界反射的暖Rossby波在8月到達(dá)斯里蘭卡東海岸,使得SLD減弱;另外,局地正風(fēng)應(yīng)力旋度的減弱也是SLD消減的因素之一.

    其次,除了SLDPWC之外,SLDPWC北部還有另外2個(gè)風(fēng)應(yīng)力渦旋(正文圖2).正的風(fēng)應(yīng)力渦旋位于印度東海岸(東印度正風(fēng)應(yīng)力旋渦,EIPWC),在EIPWC和SLDPWC之間存在一個(gè)負(fù)的風(fēng)應(yīng)力旋度(負(fù)風(fēng)應(yīng)力旋度,NWC).本文研究發(fā)現(xiàn),SLD的衰退與BBD的發(fā)展與這3個(gè)風(fēng)應(yīng)力旋度的變化密切相關(guān).首先,在EIPWC作用下BBD開(kāi)始發(fā)展,回歸分析也表明EIPWC和BBD的上升流之間存在著很強(qiáng)的正相關(guān)關(guān)系(r2=0.76,p >0.5) (正文圖9).隨著NWC的減弱,正的風(fēng)應(yīng)力旋度逐漸統(tǒng)治了孟加拉灣西南部,這也是最終SLD和BBD合并的關(guān)鍵因子,進(jìn)一步的EKE分析也證實(shí)了兩者之間的這種動(dòng)力學(xué)聯(lián)系.除此之外,次表層與SLD和BBD相關(guān)的冷中心之間的合并過(guò)程,也清晰地表明了二者之間的熱動(dòng)力學(xué)聯(lián)系.

    Mechanisms of the Sri Lanka Dome and its evolutionary aspects

    K.B.S.S.J EKANAYAKA1 WANG Weiqiang1

    1 State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology,

    Chinese Academy of Sciences,Guangzhou 510301

    AbstractThis study focuses on the formation,migration and dissipation mechanisms of the Sri Lanka Dome (SLD) and its combination with Bay of Bengal Dome (BBD) using Simple Ocean Data Assimilation reanalysis data and AVISO data.The SLD is a cyclonic eddy in the southwest Bay of Bengal,which occurs during May through September coinciding with the intrusion of southwest monsoon current into the Bay of Bengal.The results show that,Ekman pumping due to positive wind stress curl east off Sri Lanka is the governing mechanism of the SLD formation and migration during developing stage.Besides,a positive impact of cold Rossby waves propagation contribute to the SLD developing as well.Further,weakening positive wind stress curl and the eastern boundary reflected warm Rossby wave effect are important factors for dissipation of SLD.During the decaying stage of SLD,the BBD evolves independently because of the Ekman Pumping driven by local positive wind curl.The regression analysis proves that the local wind stress curl has a strong positive relationship with upwelling within the BBD.After September,weak positive wind stress curl dominates in the southwest BOB,the BBD starts to merge with the SLD.There are two kinds of links between SLD and BBD.The first is a dynamic interaction revealed by eddy kinetic energy analysis and the second is a thermodynamic interaction described by combination of subsurface cold cores of the SLD and the BBD.

    Key wordsSri Lanka Dome;wind stress curl;Ekman Pumping;Rossby;Bay of Bengal Dome

    0 Introduction

    Due to the transition pattern of the Inter-Tropical Convergence Zone (ITCZ) over the Bay of Bengal (BoB) region,there are two annually reversing monsoon currents,which can be observed as Northeast Monsoon Current (NMC) and Southwest Monsoon Current (SMC),driven by Northeast Monsoon (NEM) winds (December to April) and Southwest Monsoon (SWM) winds (May to September) respectively (Vos et al.[1]).The Sri Lanka Dome (SLD) circulation occurs during the SWM period coinciding with the intrusion of SMC eastward flow into BoB,which is one of the most important mesoscale cyclonic circulation in the Bay of Bengal.

    After passing the Sri Lankan island,part of the SMC curves around its southeastern boundary and intrude into the BoB,while the rest of SMC flow continues eastward.A part of this aforementioned BoB intrusion flows along a cyclonic path near the east coast of Sri Lanka and contributes to the SLD formation (Vinayachandran and Yamagata[2];Vinayachandran et al.[3];Schott and McCreary[4]).As a result,the western flank southward coastal current of the SLD flows against the prevailing local winds along off eastern coast of Sri Lanka (Shetye et al.[5]).The southernand eastern flanks of the SLD are closed by eastward SMC flow andBoB intrusion of SMC respectively (Vinayachandran and Yamagata[2]).

    From the climatological sea surface temperature analysis,Vinayachandran and Yamagata[2] mentioned that,the SLD starts from May and decays in September.Further,they noticed the well-developed stage of the SLD is occurring in east of Sri Lanka during July.McCreary et al.[6],Vinayachandran and Yamagata[2] and Burns et al.[7] demonstrated a positive wind stress curl occurs in the east of Sri Lanka and generates an Ekman upwelling in this region during SWM.Furthermore,Vos et al.[1] ran an idealized model to examine the hypothesis that,the SLD is formed through the interaction between the SMC and the topography (Island Wake effect).They concluded that the island wake effect is a supportive mechanism for the SLD formation.

    Although wind stress curl near the east coast of Sri Lanka remains positive until December,the arrival of warm Rossby waves propagating from the east,results in decay of the SLD due to the downwelling effect (Vinayachandran and Yamagata[2];Burns et al.[7]).Vinayachandran and Yamagata[2] reported that,these warm Rossby waves originate from the eastern boundary of BoB during April-June and reach the east coast of Sri Lanka during September.In addition to the negative impact of warm Rossby waves,Burns et al.[7] mentioned the northern BoB low salinity water invasion is also a factor for the SLD dissipation.During the decay stage of SLD,another cyclonic eddy was observed,which located northeastward to Sri Lanka named as Bay of Bengal Dome (BBD).McCreary et al.[6] showed that,the SLD cyclonic ocean circulation gradually migrates to northward during the ending of SWM and becomes the BBD.As opposed to the above evidence,Vinayachandran and Yamagata[2] showed that,the BBD independently evolves during September and devoid of further dynamic interaction with SLD.

    All aforementioned studies have well addressed the formation of the SLD,but their focus deprived of explanation about the migration of the SLD throughout the whole evolution.The current study focuses on the formation and migration of the SLD and its relevant mechanisms.Furthermore,the decay stage of the SLD seems bit controversial owing to the several demonstrations given by previous studies.The present study will attempt to explain the decay stage of the SLD and its related interaction with BBD.

    1 Data and ethodology

    The Simple Ocean Data Assimilation reanalysis version 2.0.2. is used in this study.This ocean model is based on Parallel Ocean Program physics with an average 0.25°×0.4°×40-level resolution.Observations include virtually all available ocean station data,moored temperature,hydrographic profile data and salinity time series,surface temperature and salinity observations of various types,and nighttime infrared satellite SST data.The processed data are in a monthly-composites form and they were mapped onto a uniform 0.5°×0.5°×40-level grid (Carton and Giese[8]). wind stress (N/m2),sea surface height (m),ocean temperature (including SST) (°C),salinity,horizontal ocean velocities (cm·s-1) and vertical velocity(cm·s-1) were used for the analysis.Monthly data were used for climatological and time series data analysis.All done by using NCAR Command Language (NCL) 6.3.0.

    The wind stress curl (×τ) was calculated using τy and τx component of the wind stress data according to the following equation (Kraus and Businger[9]):

    Where x and y are the longitude and latitude respectively.

    To understand the dynamic link between SLD and adjacent mesoscale eddies,the Sea Surface Height Anomaly (SSHA) and Eddy Kinetic Energy (EKE) are plotted from the Archiving Validation and Interpretation of Satellite Oceanographic (AVISO) SSHA and Geostrophic Current data.These homogeneous,inter-calibrated and highly accurate long time series AVISO data processed from various altimetric missions (Topex/Poseidon,ERS-1/2,Jason-1,Envisat and OSTM/Jason-2)(IPRC[10]).The EKE is calculated by using AVISO geostrophic current data according to the following equation:

    Where Ug and Vg are the geostrophic velocities(Jia et al.,[11]).

    Fig.1 Monthly ocean current at 35 m (vectors) from SODA overlaid on SSHA (shaded) from AVISO March to

    November.The black box (7°-9°N,83°-86°E) in June is selected to represent the region during the SLD peak

    2 Results

    As a cyclonic eddy,SLD can be distinguished with the negative SSHA due to its divergent upwelling formation.In this study,the occurrence of the cyclonic ocean circulation is studied with climatological ocean current data at 35 m overlaid with the SSHA data.The SSHA data show positive values east of Sri Lanka until May due to the anticyclonic ocean circulation (Fig.1).In May,there is a small negative circular SSHA patch and corresponding cyclonic ocean circulation,which appears first with -2 to -4 cm SSHA near the southeastern coast of Sri Lanka.From June to August,negative SSHA patch develop and become stronger more than -10 cm.In particular,it gradually move northeastward and results in the decay of SLD.During September, the SLD is vanished.Instead,a widened negative SSHA core can be seen between 10°N and 13°N resulting in the formation of BBD near the northeast coast of Sri Lanka.

    2.1 Developing stage of the Sri Lanka Dome

    Previous studies described the Ekman pumping as the major factor for the formation of SLD (McCreary et al.[6];Vinayachandran and Yamagata[2];Burns et al.[7]).They showed occurrence of positive wind stress curl generating the SLD.Although,the impact of this positive wind stress curl on the SLD migration is yet to be analyzed.

    Fig.2 Monthly wind stress curl (unit:10-8 N·m-3) from May to October (SODA).Three parallelly

    northwestward wind stress curls in June are marked by EIPWC (east Indian positive wind stress curl),NWC

    (negative wind curl) and SLDPWC (Sri Lanka Dome related positive wind stress curl)

    First,in resemblance with previous studies,F(xiàn)igure 2 also captured the strong positive wind stress curl,which can be observed during the boreal summer in the east part of Sri Lanka,with the evolution of the SLD.The patch of positive wind curl develops from the southern coast of Sri Lanka during May with a distinctive northeastward projection.From May to August,this strong positive wind stress curl gradually strengthens and further extends its projection northeastward till August (hereafter this Sri Lanka Dome related positive wind stress curl will be mentioned as SLDPWC).In the same manner,the SLD migrates northeastward via expanding its periphery from the monthly ocean current and SSHA in Figure 1.

    Further,to understand the relationship between local wind stress curl and SLD,a detailed analysis of wind stress curl is carried out over the region of 7°-9°N and 83°-86°E,which is indicated as SLD region in Jun in Figure 1.The results illustrate that,the well-developed wind stress curl over the SLD region develops from May,strengthens till August,and decay after August (Fig.3).At the same time,vertical velocity at 35 m over the SLD region shows similar variations with the wind stress curl,indicating that positive wind curl induced Ekman upwelling in the SLD region is dominant.Correspondingly,the patch of positive wind stress curl gradually retract to the southern coast of Sri Lanka at the end of SWM in October (Fig.2).Afterward,the SLD cyclonic ocean circulation deflects northwestward and merge with BBD (Fig.1).

    In addition to the SLDPWC,there are two other wind stress curls located northward to the SLDPWC,which run parallelly northeastward during boreal summer (Fig.2).The positive wind stress curl located near east coast of India (hereafter mentioned as east Indian positive wind stress curl,EIPWC).The other is a negative one,which can be seen between the SLDPWC and EIPWC (hereafter mentioned as NWC).This NWC results in the northeast coastal current of Sri Lanka flowing northward and separating from the southward western flank flow of the SLD (Fig.1).After August,this northward flow gradually decays with the weakening of this NWC.

    2.2 Dissipation of the Sri Lanka Dome

    Although the positive wind stress curl is still remaining in east of Sri Lanka in September,SLD starts to decay (Vinayachandran and Yamagata[2];Burns et al.[7]),in which the warm Rossby wave is a major factor contributing to the dissipation of the SLD.In detail,the eastward travelling equatorial Kelvin waves,generated by zonally oscillated seasonal equatorial winds,result in the generation of upwelling (cold) and downwelling (warm) Rossby waves after hitting on the eastern boundary of the BoB (Wyrtki[15];Rao et al.[16];Sreenivas et al.[17]).As a result,the warm Rossby waves start in April and reach the east coast of Sri Lanka after September (Vinayachandran and Yamagata[2];Burns et al.[7]),in this study,the downwelling Rossby waves,which start during March and reach the east coast of Sri Lanka during August.In addition to the warm Rossby wave effect,the SLDPWC weakening is also an important factor for dissipation of SLD.

    Fig.4 Time-longitude diagrams of climatological average,meridional velocity over the region of

    7°N to 8°N (left) and zonal velocity over the region of 8°N to 9°N (right)〖TS)〗

    Moreover,the impact of remote forcing Rossby wave signals on the SLD is studied using Hovmoller plots of the zonal velocity and the meridional velocity in this study (Fig.4).Here,meridional velocity over zonal belt of 7°N to 8°N and zonal velocity over zonal belt of 8°N to 9°N are selected,since these two zonal belts cover western and northern flank of SLD.The warm Rossby waves,depicted by meridional and zonal velocity in Figure 4,start from the east during March and gradually reach the east coast of Sri Lanka during August.The downwelling Rossby signals favor strengthening the northward currents and eastward currents.Thus arrival of this warm meridional velocity Rossby signals result in weakening southward western flank flow of SLD? and weakening westward northern flank flow of SLD during August (Fig.4).As a result,the anticyclonic ocean circulation and corresponding downwelling condition start to develop in this region.In addition,as a supportive mechanism to SLD development,it also can be seen from upwelling Rossby waves propagation (Fig.4),which start from January and reach the east coast of Sri Lanka during April through July.It results in strengthen the westward northern flank flow of SLD and the SLD cyclonic circulation.

    In addition to the warm Rossby waves,weakening of the positive wind stress curl is also an important factor for SLD dissipation.Figure 2 and Figure 3 clearly that,the positive wind stress curl,directly forcing SLD,gradually decays after August.Correspondingly,the upwelling of SLD also gradually weakens according to vertical velocity variation shown in Figure 3 because of the strong positive relationship between local directly forcing wind stress curl and SLD.Therefore,the combined effect of the weakening of local positive wind stress curl and the arrival of warm Rossby waves results in decaying of the SLD after August.

    2.3 SLD and BBD

    After August,the SLD moves to northeastward and merges with BBD.So far,the interaction between SLD and BBD is still controversial due to two different viewpoints.First,McCreary et al.[6] showed that,the SLD cyclonic ocean circulation gradually migrates northward during the ending of SWM and becomes the BBD during which the dynamic link between SLD and BBD is not described.As opposed to the above illustration,Vinayachandran and Yamagata[2] showed that,the BBD independently evolves from the SLD during September and devoid of further dynamic interaction with the SLD.

    In contrast,the present study shows that,the BBD evolves independently initiated by Ekman pumping resulted from local positive wind stress curl and subsequently merges with the SLD accompany with shrinkage of negative wind stress curl sandwiched between two positive wind curl.Correspondingly,the subsurface temperature clearly shows evolvement of the BBD and combination of the BBD and the SLD (Fig.5c).As a result,the upwelled cold water masses of the SLD and the BBD merge to be a large cold water gyre in northeast of Sri Lanka.Vinayachandran and Yamagata[2] also stated that,the lateral convergence is the largest contributor for the heat changes in the BBD during the formation stage.It indicates that the cold water mass of decaying SLD plays an important role in the thermodynamics of BBD.Therefore,the present study proves both dynamic and thermodynamic links between the SLD and the BBD,further investigate in detail below.

    As mentioned above,the BBD evolves independently in the formation stage.Moreover,the Ocean currents and SSHA in Figure 5b indicate that,the formation of the BBD also to westward projection of the SLD during September.This BBD formation occurs in northeast of Sri Lanka simultaneous the eastward propagation of the EIPWC (Figs.5a and 5b).A regression analysis is carried out for the Bay of Bengal Dome box (black square marked in Figure 5b) (10°-13°N,81°-83°E) to understand the relationship between local wind stress curl and the BBD formation (Fig.6).The results reveal that,the local wind stress curl has a strong positive relationship with the upwelling of the BBD (r2:0.76,p>0.5),which means the BBD forms as a result of the EIPWC driven Ekman Pumping.In order to confirm this independent evolution of the BBD,a detailed analysis of EKE is carried out from June to December (Fig.7).Since the analysis on EKE variations during the SLD decaying,for convenience,a new Sri Lanka Dome box (8°-12°N,83°-87°E) (red square in Figure 5b) is defined,which is different from the region during the SLD peak (Fig.1).The Figure 7 illustrates that,the EKE of SLD box gradually decreases after June,while the EKE of? the BBD box gradually increases until November.Most importantly,a considerably higher value of EKE (~120 cm2·s-1) can be seen from the BBD box prior to the SLD interaction with the BBD during August confirms the BBD circulation starts to occur independently from the SLD.

    Fig.5 Monthly climatological characteristics,(a) wind stress curl,(b) ocean current vectors (35 m depth) overlaid

    on SSHA (shaded),and (c) temperature at 35 m depth from September to October.Red square is for the SLD box

    Fig.6 Relationship between monthly wind stress curl and

    monthly 35 m depth velocity (cm·s-1) in

    BBD box (10°-13°N,81°-83°E) (SODA,1988-2008),

    with black line for linear regression

    With the degradation of the NWC,the positive wind stress curl dominates over southwest of the BoB and favors the combination of the weakening SLD and the newly formed BBD (Figs.5a and 5b).The spatial temperature distributions and ocean currents clearly indicate the merging of these two separate during October (Figs.5b and 5c).Furthermore,the EKE of the BBD,which remains constant (-120 cm2·s-1) during August to October,suddenly increases to the maximum level (-160 cm2·s-1) in November,whereas the EKE of the SLD gets the lowest (-10 cm2·s-1)(Fig.7).These two simultaneous EKE variations are associated with the SLD completion and the BBD peak in November,for which the SLD is decaying in November and afterward the major part of the SLD EKE moves to the BBD and favor BBD? (Fig.5b).It clearly shows that,the upwelled cold water masses of these two eddies merge owing to the interaction of their two cyclonic ocean circulations (Fig.5c).This cold water masses interaction further reveals the thermodynamic link between the SLD and the BBD relevant to their dynamic interaction.Eventually,the merging of the SLD and the BBD creates a large cold water gyre in southwest of the BoB.

    Fig.7 Monthly variation of the Eddy Kinetic Energy

    (EKE) of the Sri Lanka Dome box (red square in Fig.5b)

    (red color) and the Bay of Bengal Dome box (black

    square in Fig.5b) (blue color) from June to December based

    on monthly composites of AVISO data during 1993-2008

    3 Discussions and Conclusions

    As studied above,the SLDPWC related Ekman Pumping mechanism offers a reasonable explanation for the formation of the SLD.The regression analysis proves that,there is a strong positive relationship between the SLDPWC and the upwelling of the SLD (r2:0.93,p>0.5) (Fig.3).The cyclonic circulation,which starts from May,gradually strengthens during June to August and moves northeast accompany with northeast propagation of the SLDPWC.The present study clearly demonstrates the consistency between the SLD migration and positive wind stress curl propagation during SLD evolution.Further,Ekman upwelling of SLD results in the upward movement of subsurface high salinity cold water east of Sri Lanka during boreal summer.

    The role of eastern boundary reflected Rossby wave is importantin the SLD evolution.Cold Rossby waves propagation,observed in time-longitude diagrams of zonal velocity,favors the SLD developing.Similarly,the reflected warm Rossby waves reach the east coast of Sri Lanka during August and play a major role in dissipation of the SLD (Fig.4).In addition,the present study shows that,the weakening of SLDPWC is also an important factor for the dissipation mechanism because of the strong positive relationship between SLDPWC and the SLD.

    In addition to the SLDPWC,there are two other wind stress curls located northward to the SLDPWC,which run parallelly northeastward during boreal summer (Fig.2).One is the positive wind stress curl located near east coast of India (EIPWC),and the other is a negative wind curl (NWC) which can be seen between the related wind curl and EIPWC.The present study shows that,evolution of the SLD and the BBD are closely related to variations of these three wind curls.Regarding to the BBD,it first evolves independently from the SLD due to the EIPWC driven Ekman Pumping,for which the regression analysis confirms that the EIPWC has a strong positive relationship with the upwelling of the BBD (r2:0.76,p>0.5) (Fig.6).Although these three wind curls weaken after September,the BBD develops and peaks accompany with disappearance of NWC.As a result,weak positive wind curl is dominant over southwest of BoB in November,which favors combination of SLD and BBD.As a dynamic interaction,the EKE transfers to the BBD when SLD decaying (Fig.3).Simultaneous,also as a thermodynamic interaction,the upwelled cold water mass of the SLD and the BBD merge revealed by subsurface temperature variations (Fig.5c).It illustrates the close relationship between the SLD and the BBD,for which they have not only a thermodynamic interaction,but also a dynamic interaction.

    Acknowledgement

    We would like toexpress our special gratitude to China-Sri Lanka Joint Center for Education and Research,Chinese Academy of Sciences and University of Ruhuna (Sri Lanka) for their guidance and the administrative support.This study was funded by the Strategic Priority Research Program of Chinese Academy of Sciences under contract No.XDA20060502;the National Natural Science Foundation of China under contract Nos.41676013,41521005,41731173;the National Key Research and Development Program of China under contract No.2016YFC1401401;the Independent Research Project Program of State Key Laboratory of Tropical Oceanography under contract No.LTOZZ1702;and the CAS/SAFEA International Partnership Program for Creative Research Teams.Further we would like to extend our gratitude to NCAR,as all the plots that are included in this paper are created by using NCAR Command Language (NCL).

    References

    [1] De Vos A,Pattiaratchi C B,Wijeratne E M S.Surface circulation and upwelling patterns around Sri Lanka[J].Biogeosciences,2014,11(20):5909-5930

    [2] Vinayachandran P N,Yamagata T.Monsoon response of the sea around Sri Lanka:generation of thermal domes and anticyclonic vortices[J].J Phys Oceanogr,1998,28(10):1946-1960

    [3] Vinayachandran P N,Masumoto Y,Mikawa T,et al.Intrusion of the Southwest Monsoon Current into the Bay of Bengal[J].J Geophys Res-Oceans,1999,1041(C5):11077-11085

    [4] Schott F A,McCreary Jr J P.The monsoon circulation of the Indian Ocean[J].Prog Oceanogr,2001,51(1):1-123

    [5] Shetye S R,Shenoi S S C,Gouveia A D,et al.Wind-driven coastal upwelling along the western boundary of the Bay of Bengal during the southwest monsoon[J].Cont Shelf Res,1991,11(11):1397-1408

    [6] McCreary Jr J P,Han W,Shankar D,et al.Dynamics of the East India Coastal Current 2.Numerical solutions[J].J Geophys Res,1996,101(C6):13993-14010

    [7] Burns J M,Subrahmanyam B,Murty V S N.On the dynamics of the Sri Lanka Dome in the Bay of Bengal[J].J Geophys Res-Oceans,2017,122(9):7737-7750

    [8] Carton J A,Giese B S.A reanalysis of ocean climate using simple ocean data assimilation (SODA)[J].Mon Weather Rev,2008,136(8):2999-3017

    [9] Kraus E B,Businger J A.Atmosphere-ocean interaction[M].London:Oxford University Press,1994:362

    [10] Asia-Pacific Data-Research Center of the IPRC.AVISO maps of sea level anomalies & geostrophic velocity anomalies (MSLA)[EB/OL].[2018-01-26].http://apdrc.soest.hawaii.edu/datadoc/aviso_upd_msla_mon.php

    [11] Jia F,Wu L X,Qiu B.Seasonal modulation of eddy kinetic energy and its formation mechanism in the Southeast Indian Ocean[J].J Phys Oceanogr,2011,41(4):657-665

    [12] Vinayachandran P N,Chauhan P,Mohan M,et al.Biological response of the sea around Sri Lanka to summer monsoon[J].Geophys Res Lett,2004,31(1):1-4

    [13] Jensen T G.Arabian Sea and Bay of Bengal exchange of salt and tracers in an ocean model[J].Geophys Res Lett,2001,28(20):3967-3970

    [14] Wilson E A,Riser S C.An assessment of the seasonal salinity budget for the upper Bay of Bengal[J].J Phys Oceanogr,2016,46(5):1361-1376

    [15] Wyrtki K.An equatorial jet in the Indian Ocean[J].Science,1973,181(4096):262-264

    [16] Rao R R,Kumar M S G,Ravichandran M,et al.Interannual variability of Kelvin wave propagation in the wave guides of the equatorial Indian Ocean,the coastal Bay of Bengal and the southeastern Arabian Sea during 1993-2006[J].Deep Sea Res Pt I,2009,57(1):1-13

    [17] Sreenivas P,Gnanaseelan C,Prasad K V S R.Influence of El Nio and Indian Ocean Dipole on sea level variability in the Bay of Bengal[J].Global Planet Change,2011,80-81:215-225

    观看av在线不卡| 精品一区二区三卡| 人体艺术视频欧美日本| 成人三级做爰电影| av国产久精品久网站免费入址| 男女边摸边吃奶| 国精品久久久久久国模美| 18禁观看日本| 亚洲人成电影免费在线| 国产91精品成人一区二区三区 | 久久国产精品男人的天堂亚洲| 日韩 亚洲 欧美在线| 国产福利在线免费观看视频| 男女午夜视频在线观看| 又粗又硬又长又爽又黄的视频| 欧美日韩黄片免| 91麻豆精品激情在线观看国产 | 成人午夜精彩视频在线观看| 亚洲国产最新在线播放| 纵有疾风起免费观看全集完整版| 亚洲欧美中文字幕日韩二区| 9色porny在线观看| av福利片在线| 免费高清在线观看视频在线观看| 日韩一区二区三区影片| 人人妻人人爽人人添夜夜欢视频| 亚洲五月婷婷丁香| 国产精品一区二区在线观看99| 巨乳人妻的诱惑在线观看| 精品一区二区三区av网在线观看 | 欧美另类一区| 国产一区二区三区av在线| 中文字幕人妻丝袜制服| 成年动漫av网址| 婷婷色av中文字幕| 男人操女人黄网站| 亚洲激情五月婷婷啪啪| 久久精品久久久久久噜噜老黄| 欧美日韩亚洲国产一区二区在线观看 | a级毛片在线看网站| 精品第一国产精品| 亚洲七黄色美女视频| 国产精品免费视频内射| www.熟女人妻精品国产| 美女主播在线视频| 一区二区av电影网| 日韩大片免费观看网站| 免费不卡黄色视频| 午夜免费男女啪啪视频观看| 18禁国产床啪视频网站| 日韩制服骚丝袜av| 成年动漫av网址| 国产激情久久老熟女| 精品免费久久久久久久清纯 | 精品少妇一区二区三区视频日本电影| 国产福利在线免费观看视频| 自线自在国产av| 狠狠精品人妻久久久久久综合| 国产无遮挡羞羞视频在线观看| 欧美成人午夜精品| 欧美久久黑人一区二区| 青草久久国产| 蜜桃国产av成人99| 亚洲国产精品一区三区| 中文字幕人妻熟女乱码| 人人妻人人澡人人爽人人夜夜| 国产麻豆69| 免费高清在线观看日韩| 丝袜人妻中文字幕| 母亲3免费完整高清在线观看| 国产又爽黄色视频| 国语对白做爰xxxⅹ性视频网站| 十八禁人妻一区二区| 亚洲av美国av| 一级黄色大片毛片| www.熟女人妻精品国产| 亚洲国产av影院在线观看| 青草久久国产| 别揉我奶头~嗯~啊~动态视频 | 老司机影院毛片| 99国产精品99久久久久| 国产黄色免费在线视频| 成人黄色视频免费在线看| 国产爽快片一区二区三区| 丝瓜视频免费看黄片| 免费看十八禁软件| 国产日韩欧美在线精品| 日韩精品免费视频一区二区三区| 久久99热这里只频精品6学生| 中文字幕高清在线视频| 国产成人免费无遮挡视频| 久久久精品94久久精品| 制服人妻中文乱码| 亚洲成国产人片在线观看| 精品少妇久久久久久888优播| 亚洲精品av麻豆狂野| 热re99久久国产66热| 国产免费视频播放在线视频| 十八禁人妻一区二区| 国产成人av教育| 亚洲国产欧美在线一区| 亚洲 欧美一区二区三区| 精品少妇内射三级| 日韩精品免费视频一区二区三区| 国产激情久久老熟女| 久久亚洲精品不卡| 最近最新中文字幕大全免费视频 | 久久精品亚洲熟妇少妇任你| 国产深夜福利视频在线观看| 天堂俺去俺来也www色官网| 日本91视频免费播放| 人妻一区二区av| 中文字幕人妻丝袜制服| 在线观看免费午夜福利视频| 成人三级做爰电影| 91精品伊人久久大香线蕉| a 毛片基地| 涩涩av久久男人的天堂| 欧美激情 高清一区二区三区| 在线看a的网站| cao死你这个sao货| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜日韩欧美国产| 美女午夜性视频免费| cao死你这个sao货| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 久久99一区二区三区| 黑人猛操日本美女一级片| 午夜两性在线视频| www.999成人在线观看| 两性夫妻黄色片| 日本vs欧美在线观看视频| 欧美精品啪啪一区二区三区 | 丝袜在线中文字幕| 大片免费播放器 马上看| 少妇人妻 视频| 在线天堂中文资源库| 欧美变态另类bdsm刘玥| 高清不卡的av网站| 国产成人影院久久av| 久久精品久久久久久久性| 亚洲欧美一区二区三区国产| 超碰97精品在线观看| 满18在线观看网站| 搡老乐熟女国产| 免费av中文字幕在线| 国产精品 欧美亚洲| 国产成人av激情在线播放| 午夜福利一区二区在线看| 久久久久网色| 满18在线观看网站| 亚洲精品久久成人aⅴ小说| 激情视频va一区二区三区| 一级毛片女人18水好多 | 亚洲欧美色中文字幕在线| 首页视频小说图片口味搜索 | 国产成人精品无人区| 丝袜在线中文字幕| 亚洲,欧美,日韩| 亚洲情色 制服丝袜| 热re99久久精品国产66热6| 9热在线视频观看99| 午夜免费成人在线视频| 亚洲久久久国产精品| 激情视频va一区二区三区| 男女边摸边吃奶| 91麻豆精品激情在线观看国产 | 国产麻豆69| 少妇粗大呻吟视频| 亚洲欧美精品综合一区二区三区| 国产av国产精品国产| 久久av网站| 久久久精品免费免费高清| netflix在线观看网站| 亚洲一区二区三区欧美精品| 欧美xxⅹ黑人| 亚洲精品久久久久久婷婷小说| 日韩制服丝袜自拍偷拍| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲自偷自拍图片 自拍| 熟女少妇亚洲综合色aaa.| 亚洲,一卡二卡三卡| 亚洲精品日本国产第一区| 亚洲免费av在线视频| 久久女婷五月综合色啪小说| 欧美日韩亚洲高清精品| 爱豆传媒免费全集在线观看| 免费人妻精品一区二区三区视频| 美女大奶头黄色视频| 亚洲欧美中文字幕日韩二区| 精品久久久精品久久久| 少妇 在线观看| 日本av免费视频播放| 另类亚洲欧美激情| 日本一区二区免费在线视频| 午夜福利免费观看在线| 国产淫语在线视频| 97精品久久久久久久久久精品| 国产成人欧美| 日本91视频免费播放| 超碰97精品在线观看| 午夜久久久在线观看| 久久人人爽av亚洲精品天堂| 国产亚洲欧美在线一区二区| 欧美精品人与动牲交sv欧美| 亚洲精品一二三| 亚洲综合色网址| 波多野结衣一区麻豆| 亚洲欧洲日产国产| 国产高清不卡午夜福利| 嫩草影视91久久| 国产精品久久久久成人av| 一级a爱视频在线免费观看| 成人国产一区最新在线观看 | 日韩 欧美 亚洲 中文字幕| 色网站视频免费| 制服诱惑二区| 日韩电影二区| 亚洲欧美一区二区三区黑人| 亚洲九九香蕉| 亚洲精品国产色婷婷电影| 国产不卡av网站在线观看| 精品少妇内射三级| 中文字幕av电影在线播放| 国产视频一区二区在线看| www日本在线高清视频| videos熟女内射| 国产精品国产三级专区第一集| 日韩制服骚丝袜av| 好男人电影高清在线观看| 美女视频免费永久观看网站| 亚洲精品国产av蜜桃| 精品人妻在线不人妻| 久久精品人人爽人人爽视色| 久久精品久久精品一区二区三区| 国产在线免费精品| 精品国产一区二区三区四区第35| 婷婷色av中文字幕| 男女无遮挡免费网站观看| 午夜福利一区二区在线看| 国产国语露脸激情在线看| 亚洲人成电影观看| 欧美日韩亚洲高清精品| 午夜免费观看性视频| 亚洲国产精品一区二区三区在线| 免费高清在线观看视频在线观看| 色婷婷av一区二区三区视频| 免费久久久久久久精品成人欧美视频| 天天影视国产精品| 熟女少妇亚洲综合色aaa.| 亚洲精品在线美女| 亚洲精品在线美女| av在线老鸭窝| 成人亚洲精品一区在线观看| 国产免费现黄频在线看| 永久免费av网站大全| 欧美精品高潮呻吟av久久| 久久精品久久久久久久性| 婷婷色av中文字幕| 深夜精品福利| 人人妻人人爽人人添夜夜欢视频| 国产在线视频一区二区| 国产精品久久久久久人妻精品电影 | av在线app专区| 久久久久久亚洲精品国产蜜桃av| 国产色视频综合| 午夜久久久在线观看| 亚洲国产精品一区三区| 午夜日韩欧美国产| 精品少妇一区二区三区视频日本电影| 久久影院123| 在线看a的网站| 亚洲国产欧美一区二区综合| 最新在线观看一区二区三区 | 中文字幕精品免费在线观看视频| 国产一级毛片在线| 亚洲成人手机| 高清不卡的av网站| 国产av国产精品国产| 手机成人av网站| 啦啦啦视频在线资源免费观看| 99热全是精品| 精品亚洲成a人片在线观看| 国产成人一区二区在线| 国产成人精品无人区| 亚洲伊人久久精品综合| 一级毛片我不卡| 久久狼人影院| 在线观看www视频免费| 欧美黄色片欧美黄色片| 亚洲av美国av| bbb黄色大片| 不卡av一区二区三区| 国语对白做爰xxxⅹ性视频网站| 成年动漫av网址| 欧美成狂野欧美在线观看| 国产精品三级大全| 亚洲图色成人| 国产成人精品无人区| 人妻人人澡人人爽人人| 宅男免费午夜| 精品久久蜜臀av无| 老司机靠b影院| 你懂的网址亚洲精品在线观看| 大码成人一级视频| 国产黄色视频一区二区在线观看| 在线天堂中文资源库| 欧美日韩成人在线一区二区| 我要看黄色一级片免费的| 久久久久国产精品人妻一区二区| 18禁国产床啪视频网站| 七月丁香在线播放| 又紧又爽又黄一区二区| 国产女主播在线喷水免费视频网站| 欧美久久黑人一区二区| 女人精品久久久久毛片| 一本大道久久a久久精品| 啦啦啦 在线观看视频| 亚洲国产av影院在线观看| 在线亚洲精品国产二区图片欧美| 精品免费久久久久久久清纯 | 亚洲国产欧美一区二区综合| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 国产日韩欧美在线精品| 国产日韩欧美在线精品| 咕卡用的链子| 精品少妇一区二区三区视频日本电影| 日本五十路高清| videos熟女内射| 免费在线观看完整版高清| 九草在线视频观看| 高清黄色对白视频在线免费看| 国产国语露脸激情在线看| 各种免费的搞黄视频| 大片免费播放器 马上看| 中文字幕最新亚洲高清| 9热在线视频观看99| 国产在线观看jvid| 男人操女人黄网站| 婷婷丁香在线五月| 91精品伊人久久大香线蕉| 9191精品国产免费久久| 午夜免费成人在线视频| 大型av网站在线播放| 夫妻午夜视频| 欧美 日韩 精品 国产| 少妇的丰满在线观看| 欧美另类一区| 999久久久国产精品视频| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 亚洲国产av影院在线观看| 亚洲av国产av综合av卡| 在现免费观看毛片| 美女高潮到喷水免费观看| 国产亚洲av高清不卡| 色网站视频免费| av天堂久久9| 91成人精品电影| www日本在线高清视频| 50天的宝宝边吃奶边哭怎么回事| 后天国语完整版免费观看| 不卡av一区二区三区| 国产成人一区二区在线| av网站在线播放免费| 一级a爱视频在线免费观看| 在线观看免费高清a一片| 美女午夜性视频免费| 亚洲视频免费观看视频| 亚洲国产av影院在线观看| 欧美成人午夜精品| 2018国产大陆天天弄谢| 久久人人爽人人片av| 香蕉丝袜av| 视频区图区小说| 久久久久国产精品人妻一区二区| 久久中文字幕一级| 99香蕉大伊视频| 午夜激情av网站| 啦啦啦视频在线资源免费观看| 欧美变态另类bdsm刘玥| 超碰成人久久| 久久毛片免费看一区二区三区| 观看av在线不卡| 亚洲,一卡二卡三卡| 久热这里只有精品99| 国产成人一区二区三区免费视频网站 | 精品少妇久久久久久888优播| 在线观看一区二区三区激情| 久久天躁狠狠躁夜夜2o2o | 久久久久久久久免费视频了| 国产视频一区二区在线看| 国产野战对白在线观看| 久久久久久久大尺度免费视频| 99国产综合亚洲精品| 视频在线观看一区二区三区| 男女国产视频网站| 十八禁高潮呻吟视频| 日韩一卡2卡3卡4卡2021年| 在线观看免费午夜福利视频| 久热爱精品视频在线9| 老司机靠b影院| 欧美 日韩 精品 国产| 成人亚洲欧美一区二区av| 亚洲av欧美aⅴ国产| 无遮挡黄片免费观看| 久久人妻熟女aⅴ| 国产亚洲欧美在线一区二区| 国产极品粉嫩免费观看在线| 亚洲一卡2卡3卡4卡5卡精品中文| 日日爽夜夜爽网站| 丁香六月天网| 精品国产乱码久久久久久男人| 久久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| 人人妻,人人澡人人爽秒播 | 精品一区二区三区四区五区乱码 | 久久久久精品人妻al黑| 国产精品久久久久久精品古装| 在线 av 中文字幕| 精品国产一区二区三区四区第35| 日韩一本色道免费dvd| 精品福利永久在线观看| 中文字幕人妻丝袜制服| 国产精品亚洲av一区麻豆| 国产精品成人在线| 久久影院123| 日韩一区二区三区影片| 国产不卡av网站在线观看| 这个男人来自地球电影免费观看| 中国美女看黄片| 午夜福利影视在线免费观看| 一区二区三区乱码不卡18| 亚洲视频免费观看视频| 久久亚洲国产成人精品v| 美女中出高潮动态图| 熟女少妇亚洲综合色aaa.| 一本—道久久a久久精品蜜桃钙片| 亚洲美女黄色视频免费看| 视频区欧美日本亚洲| 99热网站在线观看| 亚洲综合色网址| 高潮久久久久久久久久久不卡| 亚洲av综合色区一区| 亚洲伊人色综图| 99香蕉大伊视频| 99精国产麻豆久久婷婷| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 一区二区三区精品91| 久久久久国产精品人妻一区二区| 黄网站色视频无遮挡免费观看| 亚洲欧美日韩高清在线视频 | 精品一区二区三区四区五区乱码 | 啦啦啦中文免费视频观看日本| 国产亚洲精品第一综合不卡| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 午夜激情av网站| 亚洲成色77777| 亚洲一卡2卡3卡4卡5卡精品中文| 97人妻天天添夜夜摸| 老汉色av国产亚洲站长工具| 视频区欧美日本亚洲| 亚洲成国产人片在线观看| 久久久久久免费高清国产稀缺| 视频在线观看一区二区三区| 91老司机精品| 亚洲精品中文字幕在线视频| 亚洲av欧美aⅴ国产| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 久久人人爽人人片av| 欧美黄色片欧美黄色片| av又黄又爽大尺度在线免费看| 国语对白做爰xxxⅹ性视频网站| 飞空精品影院首页| 久久精品国产亚洲av高清一级| 777米奇影视久久| 超碰成人久久| 在线观看免费视频网站a站| 香蕉国产在线看| 女人爽到高潮嗷嗷叫在线视频| 久久毛片免费看一区二区三区| 亚洲美女黄色视频免费看| 国产精品免费视频内射| 美女高潮到喷水免费观看| 精品久久蜜臀av无| 国产真人三级小视频在线观看| 国产精品欧美亚洲77777| 免费日韩欧美在线观看| 天堂中文最新版在线下载| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲欧美在线一区二区| 久久九九热精品免费| 高清不卡的av网站| 国产一卡二卡三卡精品| 日韩av在线免费看完整版不卡| av天堂久久9| 亚洲国产成人一精品久久久| 如日韩欧美国产精品一区二区三区| 久久久久国产精品人妻一区二区| 久久性视频一级片| 亚洲午夜精品一区,二区,三区| 一级,二级,三级黄色视频| 亚洲精品日本国产第一区| 亚洲精品美女久久久久99蜜臀 | 午夜av观看不卡| 欧美黄色淫秽网站| 国产一级毛片在线| 国产三级黄色录像| 免费av中文字幕在线| 国产成人欧美| 丝袜人妻中文字幕| 一本大道久久a久久精品| 久久99一区二区三区| 欧美日韩视频高清一区二区三区二| 欧美日韩视频高清一区二区三区二| 国产片特级美女逼逼视频| 精品人妻在线不人妻| 99香蕉大伊视频| 99九九在线精品视频| 久久国产精品男人的天堂亚洲| 国产精品国产三级专区第一集| 国产精品 欧美亚洲| 飞空精品影院首页| 18在线观看网站| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 免费观看a级毛片全部| 建设人人有责人人尽责人人享有的| 五月开心婷婷网| 波多野结衣一区麻豆| 亚洲精品av麻豆狂野| 国产淫语在线视频| 日韩人妻精品一区2区三区| 亚洲欧洲日产国产| 91国产中文字幕| 男女国产视频网站| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品中文字幕在线视频| 男女之事视频高清在线观看 | videos熟女内射| 日本色播在线视频| netflix在线观看网站| 午夜福利视频在线观看免费| 日本午夜av视频| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 啦啦啦在线免费观看视频4| 精品国产超薄肉色丝袜足j| 80岁老熟妇乱子伦牲交| 男的添女的下面高潮视频| 亚洲精品国产色婷婷电影| 国产精品久久久人人做人人爽| 日韩制服骚丝袜av| 丝袜在线中文字幕| 亚洲欧美一区二区三区黑人| 50天的宝宝边吃奶边哭怎么回事| av天堂在线播放| 99热网站在线观看| 亚洲三区欧美一区| 在线看a的网站| 欧美亚洲 丝袜 人妻 在线| 亚洲精品久久午夜乱码| 欧美人与善性xxx| xxxhd国产人妻xxx| 午夜福利,免费看| 亚洲中文日韩欧美视频| 久久久精品免费免费高清| 久久综合国产亚洲精品| 国产成人系列免费观看| 免费高清在线观看视频在线观看| 国产欧美日韩一区二区三区在线| 亚洲人成电影观看| 久热爱精品视频在线9| 成人手机av| 亚洲精品日韩在线中文字幕| 99热全是精品| 天天躁夜夜躁狠狠久久av| 亚洲精品成人av观看孕妇| 大香蕉久久网| 久久久国产一区二区| 啦啦啦啦在线视频资源| 一级毛片黄色毛片免费观看视频| 国产精品一区二区免费欧美 | 大陆偷拍与自拍| av片东京热男人的天堂| 老司机靠b影院| 男女之事视频高清在线观看 | 亚洲精品中文字幕在线视频| cao死你这个sao货| 久久国产亚洲av麻豆专区| 巨乳人妻的诱惑在线观看| 操美女的视频在线观看| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av涩爱| 精品久久久精品久久久| 天天添夜夜摸| 一级a爱视频在线免费观看| 日韩一本色道免费dvd| 久久人人爽人人片av| 亚洲av日韩精品久久久久久密 | e午夜精品久久久久久久| 我要看黄色一级片免费的| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 天堂俺去俺来也www色官网| 亚洲久久久国产精品| av网站在线播放免费| 久久国产精品影院| 国产成人一区二区三区免费视频网站 |