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

    冬季黑潮延伸體區(qū)域渦旋尺度海面溫度季節(jié)內變化特征?

    2020-01-07 16:14:58楊小繪賈英來
    關鍵詞:黑潮邊界層渦旋

    楊小繪,賈英來

    (中國海洋大學海洋與大氣學院 海洋氣象學系,山東 青島 266100)

    黑潮延伸體區(qū)域是中緯度海洋-大氣相互作用的關鍵區(qū)域[1-2]。Wang等[3]發(fā)現(xiàn)黑潮延伸體區(qū)域大尺度的海面溫度異常存在季節(jié)內變化特征,并且通過海-氣相互作用驅動北太平洋海盆尺度的風應力旋度,構成耦合模態(tài),該季節(jié)內振蕩信號在夏季最強。黑潮延伸體區(qū)域海洋渦旋活躍,渦旋通過海氣相互作用影響其上層大氣。而海氣相互作用的關鍵是海-氣界面間的熱通量交換,熱通量和SST間的關系在不同空間尺度(比如大尺度和渦旋尺度)上并不相同,Li等[4]用懲罰回歸分析(Penalized spectral regression,PSR)得出湍流熱通量和SST異常的關系在渦旋尺度(<1 000 km)上同位相變化且關系密切,是渦旋尺度SST強迫大氣;而在海盆尺度上,熱通量和SST呈負相關[5-6],這時為大氣強迫海洋。Smirnov等[7]和Putrasahan等[8]通過利用數(shù)值模式也得到類似結論。

    前人發(fā)現(xiàn)海洋渦旋可以引起海氣界面熱量通量以及邊界層大氣中風速等要素的響應[9-10],該響應在冬季最強[11-12],其主要通過垂直混合機制影響大氣[13-15]。通過渦旋合成的方法[16-18]發(fā)現(xiàn)湍流熱通量異常在渦旋中心附近更強,在渦旋外部最小[17]。在暖渦上空風速增大,冷渦上空風速減小[18]。并且渦旋對大氣的影響存在季節(jié)變化,冬春最大、夏秋較弱[11]。但是,前人研究海洋渦旋上空的海-氣相互作用過程時,多通過渦旋合成的方法來進行,并未探討其時間變化特征。黑潮延伸體區(qū)域的流軸彎曲(Meander)存在3~60 d的季節(jié)內變化信號[19],該季節(jié)內變化信號是否也存在于渦旋尺度SST的時間變化之中?為探討黑潮延伸體區(qū)域渦旋尺度SST的時間變化特征,本文利用觀測資料和再分析資料對大尺度和渦旋尺度SST的時間變化進行了分析,比較了兩者在季節(jié)內變化上的區(qū)別,并進一步分析了和冬季渦旋尺度SST季節(jié)內變化相關的大氣邊界層的響應特征。

    1 資料和方法

    本文所用SST 數(shù)據(jù)有歐洲中期天氣預報中心ECMWF(European center for medium-range weather forecasts)的再分析資料ERA-interim SST數(shù)據(jù)(水平分辨率為0.75°×0.75°,時間間隔為6 h)和OISST(NOAA Optimum Interpolation Sea Surface Temperature,OISST)觀測數(shù)據(jù)(水平分辨率為0.25°×0.25°的逐日數(shù)據(jù))。取兩者共同的時間段(2002年1月~2017年間)進行研究。因ERA數(shù)據(jù)不僅包含SST,而且包含大氣邊界層的各變量,本文中主要選取ERA的數(shù)據(jù)集進行分析。本文中冬季定義為每年的11月1日~次年3月31日,夏季定義為每年的5月1日~9月30日。根據(jù)黑潮延伸體區(qū)域海洋渦旋和鋒面的特征尺度[20],采用5°×5°二維空間濾波(2D Boxcar)方法分離渦旋尺度的SST異常和大尺度的SST異常信號。利用小波功率譜分析方法研究時間變化顯著周期,采用Lanczos帶通濾波器對資料進行帶通濾波以得到相應的頻段。為對季節(jié)內時間變化信號進行研究,所有數(shù)據(jù)均進行去(2002—2017年)氣候態(tài)月平均處理后再進行相應的后續(xù)處理。

    采用ERA-interim再分析數(shù)據(jù)的海表面潛熱通量(Surface Latent Heat Flux,SLHF)、海表面感熱通量(Surface Sensible Heat Flux,SSHF)、蒸發(fā)(Evaporation,E)、10 m風場(10 metre speed,SPE)、2 m氣溫(2 metre temperature,2mT)等數(shù)據(jù)分析大氣邊界層對SST信號的響應。采用與SST相同的方法對大氣邊界層各變量進行渦旋尺度信號和大尺度信號的分離。為考察渦旋尺度SST異常信號引起的大氣邊界層響應,采用二維模態(tài)相關(Pattern correlation)的方法進行分析兩者在空間分布形態(tài)上的一致性。關于渦旋引起的大氣的響應,以前的文獻主要通過渦旋合成分析[9-10,16-19],但是渦旋合成后的場是多渦旋平均的結果,并不能探討渦旋及大氣響應兩者關系的時間變化。模態(tài)相關,即兩個變量之間的線性相關的Pearson積矩系數(shù)(Pearson product-moment coefficient)用于衡量兩個變量各自模態(tài)之間的相關性。其主要適用于分析兩個異常場空間分布形態(tài)的異同。若兩個異常場(比如渦旋SST異常場和SLHF異常場)各異常中心的位置和形狀近似,則相關系數(shù)高。利用模態(tài)相關可以研究渦旋SST異常信號引起的大氣異常響應。若兩者相關系數(shù)大,則說明渦旋SST異常能夠引起大氣變量在空間同位相的響應。模態(tài)相關系數(shù)的時間變化還可以顯示大氣對SST異常響應的時間變化。在本文中,我們將利用模態(tài)相關分析揭示渦旋SST異常及大氣響應的季節(jié)內變化特征。

    2 渦旋尺度海面溫度季節(jié)內振蕩的時空分布特征

    為說明黑潮-親潮延伸體區(qū)域大尺度海溫與渦旋尺度海溫在季節(jié)內變化特征上的不同,圖1給出了SSTA(Sea Surface Temperature Anomaly)場90 d滑動標準差在冬季、夏季的空間分布。兩種數(shù)據(jù)的結果都表明,大尺度SST的季節(jié)內信號具有夏季強、冬季弱的特點(見圖1(a)、(c)及(e)、(g)),其中大尺度SST季節(jié)內變化信號的空間分布和Wang 等[3]的結果一致。渦旋尺度SST的季節(jié)內變化信號則是冬季強、夏季弱(見圖1(b)、(d)及(e)、(f))。大尺度和渦旋尺度SST季節(jié)內信號標準差的空間分布也不相同,渦旋尺度信號主要分布在黑潮-親潮海洋鋒面區(qū)域,在日本沿岸最強,并沿黑潮、親潮鋒面向東延伸。在兩個數(shù)據(jù)中冬季渦旋尺度SST信號的標準差可達1 ℃左右,尤其是在黑潮、親潮鋒面區(qū)域,是SST的主要變化信號。OISST數(shù)據(jù)中渦旋尺度季節(jié)內信號的標準差更大可能是因為其分辨率高的原因。

    為進一步考察渦旋尺度季節(jié)內信號的變化周期,選取黑潮海洋鋒面上的一點(146.75°E,36.25°N),對其SST時間序列進行小波分析(見圖2)。大尺度SST的顯著周期以20~40 d為主,而渦旋尺度SST的顯著周期則在40 d以上,兩者的小波譜值均存在明顯的季節(jié)變化。為進一步考察黑潮-親潮延伸體區(qū)SST的季節(jié)內變化周期,分0~40 d和40~100 d兩個季節(jié)變化顯著的頻段進行研究,圖3給出了兩個頻段內大尺度和渦旋尺度SST信號在夏季和冬季期間小波譜值最大值的空間分布,以考察該頻段周期信號的特征。對于大尺度SST信號,0~40 d頻段的小波譜值主要在日本沿岸較強,而40~100 d頻段的信號還多分布在黑潮延伸體下游。大尺度SST信號的0~40 d變化和40~100 d信號均在夏季最明顯(見圖3),這和Wang 等[3]關于該區(qū)域大尺度SST季節(jié)內變化的結論一致。對于渦旋尺度SST變化信號,40~100 d頻段的小波譜值(約15~20 ℃2)明顯強于0~40 d頻段的(6 ℃2以內)(見圖4),并且渦旋尺度SST的季節(jié)內變化在冬季明顯強于夏季。冬季渦旋尺度SST的40~100 d振蕩的顯著信號主要分布在黑潮、親潮海洋鋒面處 (見圖4),這可能和黑潮流軸彎曲(Meander)中的季節(jié)內信號有關[19]。

    3 冬季渦旋尺度海面溫度的季節(jié)內變化信號對邊界層大氣的影響

    由上述可知,大尺度SST與渦旋尺度SST的季節(jié)內振蕩信號并不一致,在冬季,黑潮-親潮延伸體區(qū)域SST變化主要是渦旋尺度SST的40~100 d變化造成的。根據(jù)前人研究[9-10],渦旋尺度SST變化可以引起海氣界面熱量通量以及邊界層大氣中風速等要素的響應,那么,渦旋尺度SST的季節(jié)內變化能否引起大氣邊界層各要素的相應變化呢?為回答這個問題,本文首先考察了黑潮鋒面一點(146.75°E,36.25°N)上的SLHF的小波功率譜分布(見圖2(e)、(f)、(g)、(h))。

    (填色部分為海面溫度90 d滑動平均的標準差,等值線為相應季節(jié)的海面溫度多年平均氣候態(tài)值(單位:℃)。圖中的藍色點為下文分析中選取的點(146.75°E,36.25°N)。The shaded parts are 90-day running mean standard deviation of SST and the contours are climatological distribution of SST,(Unit:℃.) The blue point (146.75°E,36.25°N) indicate the token point for the following analysis.)
    圖1 海面溫度標準差圖
    Fig.1 The average standard deviation distribution of SST

    ((a、c、e、g)中黑色等值線圈住的填色區(qū)域為通過95%信度水平的小波功率譜值,黑色虛線為影響錐的部分曲線,黑色曲線內為可信區(qū)域。(b、d、f、h)中藍色實線為小波全譜,紅色虛線為95%信度水平的標準紅噪音譜。The dotted lines in (a,c,e,g) are part of conial and the conial area is the credible part;The solid line is the period that passes the 95% confidence level.Solid lines in (b,d,f,h) represent the wavelet full spectrum of the oscillation,and the dashed lines indicate the standard red-noise spectra at the 95% confidence level.)
    圖2 點(146.75°E,36.25°N)的小波功率譜以及小波全譜圖
    Fig.2 Wavelet spectrum and wavelet full spectrum on the point (146.75°E,36.25°N)

    (填色部分為大尺度海面溫度的小波全譜的最大譜值,單位:℃2,等值線部分為相應季節(jié)的多年海面溫度平均值,單位:℃。The shaded parts are maximum wavelet full spectrum power of SST (Unit:℃2) and the contours are climatological distribution of corresponding SST (Unit:℃).)
    圖3 大尺度海面溫度的小波全譜最大譜值平面分布圖
    Fig.3 Distribution of the large scale SST maximum wavelet power

    (圖中的實線方框為下文分析所選取的關鍵區(qū)域,范圍為(142°E~160°E,33°N~44°N)。The rectangle indicate the key region for the following analysis:(142°E~160°E,33°N~44°N).)
    圖4 同圖3,但填色部分為渦旋尺度海面溫度
    Fig.4 As in Fig.3,but for mesoscale SST

    由圖可見,大尺度SLHF及渦旋尺度SLHF均存在季節(jié)內振蕩,大尺度SLHF的顯著周期以40 d以上為主,峰值在60~70 d。渦旋尺度SLHF的顯著周期也具有60~100 d的顯著周期,兩者的小波譜值均存在明顯的季節(jié)變化。其中渦旋尺度SLHF與渦旋尺度SST的功率譜結構相似。

    為進一步探討渦旋尺度SST和SLHF在季節(jié)內變化上的相關性,將冬季渦旋尺度SST和大氣邊界層各要素的異常進行20~100 d帶通濾波以提取季節(jié)內變化信號,并對其進行區(qū)域平均后(142°E~160°E,33°N~44°N)計算了超前滯后相關系數(shù)(見圖5)。為和大尺度信號進行對比,同時給出了大尺度SST和大氣邊界層各要素間的相關關系。圖5明確表明,大尺度SST與大氣邊界層各要素間的相關性比渦旋尺度信號間的相關性弱得多。大氣邊界層各變量和SST相關系數(shù)低且呈現(xiàn)SST超前或滯后的相關(大氣強迫海洋)。而渦旋尺度SST信號則和SLHF、SSHF、BLH、E,TSA以及SPE都是同期相關,并且渦旋尺度SST信號和SLHF、SSHF、BLH、E,TSA以及SPE這幾個變量間的同期相關系數(shù)高達0.8。兩者渦旋尺度信號的變化一致性與前人用渦旋追蹤合成方法得出的結果一致[9-10,16-19]。這說明冬季渦旋尺度SST信號能夠引起大氣邊界層各要素的同期響應,冬季,黑潮延伸體區(qū)域,在渦旋尺度上,是海洋SST驅動大氣。

    (黑色實線為渦旋尺度SST與渦旋尺度大氣要素的超前滯后相關系數(shù),黑色虛線為大尺度SST與大尺度大氣要素的超前滯后相關系數(shù)。The solid line and;Dashed line indicate the mesoscale and large scale,respectively.)
    圖5 冬季(142°E~160°E,33°N~44°N)區(qū)域內平均的海面溫度與大氣要素的超前滯后相關
    Fig.5 Lead lag correlation between area-mean SST and atmospheric elements in the key region (142°E~160°E,33°N~44°N) during winter

    上文發(fā)現(xiàn)冬季渦旋尺度SST存在季節(jié)內變化信號,當渦旋尺度SST發(fā)生變化時,大氣要素產生與之一致的響應。前人關于渦旋引起大氣的響應主要通過渦旋合成分析[9-10,16-19],但是渦旋合成后的場是多渦旋平均的結果,這并不能探討大氣對渦旋的響應是否隨著渦旋SST信號的時間而產生相應的變化。那么渦旋尺度SST的季節(jié)內信號能否引起與大氣響應的季節(jié)內變化?

    為回答這個問題,本文使用二維模態(tài)相關分析來研究這兩者之間的時間變化關系。冬季(142°E~160°E,33°N~44°N)的渦旋場與大氣邊界層各要素響應的時間變化關系如圖6藍色實線所示。由圖可見,對與渦旋尺度SST同位相時間相關密切的熱通量、蒸發(fā)和海氣溫差項,模態(tài)相關系數(shù)同樣較強,高達0.8以上,可見兩者在空間分布上也非常一致。而邊界層高度和海面風場與渦旋尺度SST的關系則較差。這說明渦旋尺度SST的季節(jié)內變化信號能夠引起大氣邊界層各要素,尤其是熱通量的相應的季節(jié)內變化。同時,兩者的模態(tài)相關系數(shù)的時間變化也說明在有的年份,兩者的關系較好,在有的時間段,兩者的關系較差。那么大氣邊界層對渦旋尺度SST異常的響應強度與什么因素有關系呢?為此,統(tǒng)計了上述區(qū)域較強渦旋(即暖渦中心值>0.6 ℃與冷渦中心值<-0.6 ℃)總個數(shù)的時間變化(見圖6黑色實線)??梢钥吹捷^強渦旋的個數(shù)與大氣對SST異常響應密切相關,強渦旋多的時候大氣邊界層異常與SST異常在空間分布上更一致,這時候由大氣邊界層更多地受渦旋SST的影響,而渦旋少的時候大氣邊界層對渦旋SST異常的響應較弱。圖6紅色實線給出的是點(146.75°E,36.25°N)渦旋尺度SST變化在40~100 d頻段內的小波功率譜值時間變化序列,可以看出,渦旋SST季節(jié)內變化信號明顯的時間段模態(tài)相關系數(shù)較大,強渦個數(shù)較多;而渦旋SST季節(jié)內變化信號不明顯的時間段,強渦個數(shù)較小,模態(tài)相關系數(shù)也較小。這說明渦旋SST季節(jié)內變化表現(xiàn)為強渦個數(shù)的變化,而強渦的增多也加強了大氣邊界層的響應。Shan和Dong[21]利用WRF(Weather research and forecasting model)模式理想試驗的結果揭示了SST異常信號的增強會加強大氣的響應,這和我們的結果是一致的。為更直觀地進行驗證,以SLHF為例,選取了6個不同的時間做了分析。參見圖6(a)中的3個紅點(代表模態(tài)相關系數(shù)大值點)和3個綠點(代表小值點)。圖7給出了這6個時間的渦旋SST場以及SLHF場的空間分布。由圖可見,在渦旋SST較弱的時間,大氣的渦旋尺度異常信號仍然存在但是分布位置和形態(tài)與渦旋SST并不一致,較弱的渦旋不易對大氣邊界層產生影響。而較強渦旋上空則大多會出現(xiàn)海表面熱通量的響應。并且強渦旋大多沿著黑潮、親潮鋒面分布,這與圖1(d)渦旋尺度SST季節(jié)內信號沿鋒面分布也是一致的??傊?,較強渦旋少的時候模態(tài)相關系數(shù)低,較強渦旋多的時候模態(tài)相關系數(shù)高。在該區(qū)域較強渦旋信號個數(shù)較多的時間段,SST異常的季節(jié)內變化更明顯,大氣與渦旋尺度SST信號的相關程度更大。

    (藍色縱軸代表關鍵區(qū)內渦旋尺度海面溫度與大氣各要素模態(tài)相關系數(shù)90 d滑動平均的大小,黑色縱軸代表關鍵區(qū)內較強渦旋的個數(shù),紅色縱軸代表點(146.75°E,36.25°N)40~100 d小波譜的平均譜值。(a)中的三個紅點代表模態(tài)相關系數(shù)的高值點,三個綠點代表模態(tài)相關系數(shù)的低值點。The blue solid lines indicate 90 d running mean pattern correlation curves between mesoscale SST and atmospheric elements in the key region.Black lines indicate the number of more stronger eddies (the centers maximum of warm eddies more than 0.6 ℃ and the centers minimum of cold eddies less than -0.6 ℃) 90 d running mean curves in the key region (142°E~160°E,33°N~44°N).Red lines indicate the SST mean power between 40~100 d on the point (146.75°E,36.25°N)90 d running mean curves.In (a) the three red points are the high points and the three green points are the low points in the temporal variation of pattern correlation coefficients.)
    圖6 冬季模態(tài)相關系數(shù)、較強渦旋個數(shù)以及40~100 d的小波譜譜值的時間序列
    Fig.6 Times series of pattern correlation,the number of stronger eddies and the mean wavelet power between 40~100 d in winter

    (填色部分為渦旋尺度海面溫度異常,單位:℃;紅藍等值線為渦旋尺度海表面潛熱通量異常(海表面潛熱通量正值代表海洋從大氣中吸收熱量),單位:W/m2;黑色等值線為冬季的海面溫度平均值,單位:℃。(a,b,c)為模態(tài)相關系數(shù)的高值點,(d,e,f)為模態(tài)相關系數(shù)的低值點。The shaded parts are mesoscale SST variability and the contours are mesoscale SLHF variability (the positive value of SLHF means downward into the ocean).Distributions of the high points of pattern correlation coefficients are present in (a,b,c) and low points in (d,e,f).)
    圖7 渦旋尺度海面溫度和海表面潛熱通量異常的分布圖
    Fig.7 Distribution of mesoscale SST variability and mesoscale SLHF variability

    冬季黑潮延伸體區(qū)域盛行西北風,把歐亞大陸的干冷空氣吹到溫暖潮濕的海洋上,這種情況下,黑潮延伸體區(qū)域會有強的向上的湍流熱通量釋放,從而會對上空大氣造成重要影響。這種釋放是通過海洋的中尺度結構海洋鋒面和海洋渦旋來完成的,所以由海洋渦旋引起的海面溫度異常有重要影響。為更進一步驗證黑潮-親潮延伸體區(qū)渦旋尺度SST的季節(jié)內變化信號會引起SLHF的響應,同樣分0~40 d和40~100 d兩個季節(jié)內變化顯著的頻段進行研究。圖8和9給出了這兩個頻段內大尺度和渦旋尺度SLHF信號在夏季和冬季期間小波譜值的最大值分布,以考察相應頻段周期的強弱。與大尺度SST不一致,大尺度SLHF信號40~100 d頻段的小波譜值明顯強于0~40 d頻段,并且大尺度SLHF信號的40~100 d變化在冬季最明顯,信號大值區(qū)的分布位置主要在黑潮流軸以南。對于渦旋尺度SLHF變化信號,40~100 d 頻段的小波譜值明顯強于0~40 d頻段,在黑潮、親潮海洋鋒面處,冬季也明顯存在40~100 d振蕩的顯著信號分布(見圖9),與渦旋尺度SST的變化一致。這進一步說明渦旋尺度SST的季節(jié)內變化能夠引起SLHF的相應變化,使SLHF在黑潮、親潮海洋鋒面處出現(xiàn)相似的季節(jié)內變化周期。從SLHF對渦旋尺度SST的響應,可以進一步證實渦旋尺度SST對大氣要素的影響較強,冬季渦旋尺度SST的季節(jié)內振蕩周期以40~100 d周期為主,并且會引起大氣邊界層的響應,激發(fā)海氣界面湍流熱通量、海面氣溫、邊界層高度等同位相的季節(jié)內振蕩。

    圖8 同圖3,但填色部分為大尺度海表面潛熱通量(單位:(W/m2)2)Fig.8 As in Fig.3,but shaded for large scale SLHF (Unit:(W/m2)2)

    4 結論與討論

    本文利用ERA_interim再分析資料以及OISST高分辨率海面溫度衛(wèi)星觀測數(shù)據(jù),以及小波分析、二維模態(tài)相關(Pattern correlation)分析等方法系統(tǒng)地分析了黑潮延伸體區(qū)域渦旋尺度SST信號的季節(jié)內變化特征及其對邊界層大氣的影響。結果表明,渦旋尺度SST的季節(jié)內變化和大尺度的信號在空間分布和季節(jié)變化上均存在顯著差別,渦旋尺度SST季節(jié)內振蕩信號在冬季最強,夏季最弱,而大尺度信號則在夏季最強。渦旋尺度SST季節(jié)內變化信號主要分布在黑潮-親潮海洋鋒面區(qū)域,在日本沿岸振幅最強,并沿黑潮、親潮鋒面向東延伸,其標準差高達1 ℃,是該海域冬季SST的主要變化信號。渦旋尺度SST的季節(jié)內振蕩周期以40~100 d周期為主,并且會引起大氣邊界層的響應,激發(fā)海氣界面湍流熱通量、海面氣溫、邊界層高度等同位相的季節(jié)內振蕩。強渦旋少的時候渦旋尺度SST信號和海氣界面湍流熱通量間模態(tài)相關系數(shù)低,強渦旋多的時候兩者模態(tài)相關系數(shù)高。在黑潮延伸體上游區(qū)域強渦個數(shù)較多的時間段,SST異常的季節(jié)內變化更明顯,邊界層大氣與渦旋尺度SST信號的相關程度更大。

    圖9 同圖3,但填色部分為渦旋尺度海表面潛熱通量(單位:(W/m2)2)Fig.9 As in Fig.3,but shaded for mesoscale SLHF (Uint:(W/m2)2)

    渦旋尺度SST的季節(jié)內變化信號和黑潮延伸體區(qū)域的鋒面變化以及該區(qū)域的海洋渦旋有著密切的關系,那么,該信號是否也存在于海平面高度場中?這是我們今后要研究的問題。另外,該SST變化信號是否能夠通過影響大氣邊界層的變化,尤其是和大氣之間的熱量通量而進一步影響上層大氣也值得研究。

    猜你喜歡
    黑潮邊界層渦旋
    譜松弛動力降尺度方法及其在呂宋海峽黑潮模擬中的應用
    海洋通報(2022年6期)2023-01-07 02:48:56
    基于PM算法的渦旋電磁波引信超分辨測向方法
    基于HIFiRE-2超燃發(fā)動機內流道的激波邊界層干擾分析
    光渦旋方程解的存在性研究
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計數(shù)據(jù)的黑潮變異特征
    海洋科學(2016年1期)2016-07-06 09:42:34
    東海冷渦對黑潮熱輸運年代際躍變的響應
    一類具有邊界層性質的二次奇攝動邊值問題
    變截面復雜渦旋型線的加工幾何與力學仿真
    非特征邊界的MHD方程的邊界層
    亚洲免费av在线视频| 俄罗斯特黄特色一大片| 一本一本久久a久久精品综合妖精| 久热这里只有精品99| 搡老岳熟女国产| 国产高清激情床上av| 久久天堂一区二区三区四区| av片东京热男人的天堂| 叶爱在线成人免费视频播放| 成年动漫av网址| 国产精品国产av在线观看| 日韩精品免费视频一区二区三区| 国产淫语在线视频| 精品国产一区二区三区四区第35| 亚洲精品一卡2卡三卡4卡5卡| 老司机在亚洲福利影院| 欧美 亚洲 国产 日韩一| 老司机深夜福利视频在线观看| 午夜91福利影院| 亚洲专区中文字幕在线| 久久久国产精品麻豆| 午夜精品久久久久久毛片777| 久久人妻av系列| av免费在线观看网站| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| av又黄又爽大尺度在线免费看| 久久国产精品人妻蜜桃| 少妇猛男粗大的猛烈进出视频| 午夜两性在线视频| 欧美黑人精品巨大| 欧美日韩亚洲综合一区二区三区_| 免费观看人在逋| 国产成人av教育| 亚洲国产av新网站| 久久青草综合色| 超碰97精品在线观看| kizo精华| 国产精品偷伦视频观看了| 多毛熟女@视频| 一边摸一边抽搐一进一出视频| 亚洲国产毛片av蜜桃av| 国产成人精品无人区| 首页视频小说图片口味搜索| 老司机亚洲免费影院| 一个人免费看片子| 亚洲一码二码三码区别大吗| 18禁裸乳无遮挡动漫免费视频| 一区二区三区乱码不卡18| 韩国精品一区二区三区| 国产精品.久久久| 久久久国产成人免费| 精品一区二区三区av网在线观看 | 人妻 亚洲 视频| 在线看a的网站| 久久久久国产一级毛片高清牌| av电影中文网址| 十八禁高潮呻吟视频| 精品久久久久久久毛片微露脸| 99热国产这里只有精品6| 国产成人av激情在线播放| 国产日韩欧美视频二区| 十八禁网站网址无遮挡| 热99re8久久精品国产| 午夜福利视频精品| 欧美日韩精品网址| 亚洲国产看品久久| a在线观看视频网站| 在线播放国产精品三级| 曰老女人黄片| 欧美变态另类bdsm刘玥| 欧美成狂野欧美在线观看| 久久国产精品人妻蜜桃| 国产深夜福利视频在线观看| 国产av精品麻豆| 国产精品秋霞免费鲁丝片| 国产免费视频播放在线视频| 新久久久久国产一级毛片| 久久国产精品大桥未久av| 两人在一起打扑克的视频| 精品国产国语对白av| 亚洲少妇的诱惑av| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品男人的天堂亚洲| 在线看a的网站| 国产男女内射视频| 国产一区二区三区视频了| 捣出白浆h1v1| 日本vs欧美在线观看视频| 一本综合久久免费| 午夜福利一区二区在线看| 久久久久久久久久久久大奶| 天堂8中文在线网| 少妇粗大呻吟视频| 黑丝袜美女国产一区| 日本撒尿小便嘘嘘汇集6| 美国免费a级毛片| 日韩成人在线观看一区二区三区| 久久亚洲真实| 高潮久久久久久久久久久不卡| 美女视频免费永久观看网站| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 在线十欧美十亚洲十日本专区| 男女无遮挡免费网站观看| 亚洲国产欧美一区二区综合| 午夜免费鲁丝| 亚洲国产毛片av蜜桃av| 十八禁人妻一区二区| 免费在线观看日本一区| 波多野结衣av一区二区av| 精品免费久久久久久久清纯 | 人人妻人人澡人人看| 精品国产亚洲在线| 亚洲精品一卡2卡三卡4卡5卡| 国产日韩欧美视频二区| 免费在线观看影片大全网站| 国产又色又爽无遮挡免费看| 老汉色av国产亚洲站长工具| 亚洲欧美一区二区三区黑人| 视频在线观看一区二区三区| 欧美精品一区二区免费开放| 一区二区av电影网| 欧美中文综合在线视频| 脱女人内裤的视频| 我要看黄色一级片免费的| 在线观看www视频免费| 午夜福利视频在线观看免费| 久久久久精品人妻al黑| 亚洲七黄色美女视频| 男女下面插进去视频免费观看| 丰满饥渴人妻一区二区三| 久久精品亚洲熟妇少妇任你| 日韩一区二区三区影片| 国产在线一区二区三区精| 国产精品国产av在线观看| 极品教师在线免费播放| 1024视频免费在线观看| 久久性视频一级片| 亚洲国产中文字幕在线视频| 久久久久久久久免费视频了| 男男h啪啪无遮挡| av免费在线观看网站| 乱人伦中国视频| 久久精品aⅴ一区二区三区四区| 亚洲av电影在线进入| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美一区二区综合| 国产99久久九九免费精品| 精品人妻1区二区| 男女免费视频国产| 窝窝影院91人妻| 女人高潮潮喷娇喘18禁视频| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频 | 久久精品亚洲av国产电影网| 一区二区三区乱码不卡18| 午夜福利乱码中文字幕| 精品国产乱码久久久久久小说| 国产在视频线精品| 国产伦理片在线播放av一区| 国产精品亚洲av一区麻豆| 久热爱精品视频在线9| 欧美日本中文国产一区发布| 极品人妻少妇av视频| 757午夜福利合集在线观看| 18禁国产床啪视频网站| 十分钟在线观看高清视频www| 黄色视频在线播放观看不卡| 国产免费现黄频在线看| 国产精品一区二区在线不卡| 丝袜喷水一区| 男人操女人黄网站| 日韩免费高清中文字幕av| 在线观看免费视频网站a站| 国产精品偷伦视频观看了| 亚洲国产欧美网| 亚洲国产看品久久| 国产欧美日韩综合在线一区二区| av免费在线观看网站| 亚洲精品美女久久av网站| 国产精品美女特级片免费视频播放器 | 动漫黄色视频在线观看| 国产精品自产拍在线观看55亚洲 | 国产老妇伦熟女老妇高清| a级片在线免费高清观看视频| 99久久国产精品久久久| 老熟妇仑乱视频hdxx| 黑丝袜美女国产一区| 高清在线国产一区| 亚洲全国av大片| 久久99一区二区三区| 淫妇啪啪啪对白视频| 亚洲 国产 在线| 人成视频在线观看免费观看| 一级毛片精品| 久久久久久人人人人人| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 亚洲第一av免费看| 人成视频在线观看免费观看| 精品亚洲成a人片在线观看| 色综合欧美亚洲国产小说| 菩萨蛮人人尽说江南好唐韦庄| 一夜夜www| 国产精品免费一区二区三区在线 | 女人爽到高潮嗷嗷叫在线视频| 最近最新中文字幕大全电影3 | 99九九在线精品视频| 午夜福利影视在线免费观看| av欧美777| 欧美另类亚洲清纯唯美| 久久久国产欧美日韩av| 在线av久久热| 国产精品影院久久| 国产不卡一卡二| 麻豆av在线久日| 搡老乐熟女国产| 丝袜在线中文字幕| 一个人免费在线观看的高清视频| 国产单亲对白刺激| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 国产欧美日韩综合在线一区二区| av电影中文网址| 一区福利在线观看| 最新的欧美精品一区二区| 久久青草综合色| 亚洲美女黄片视频| 人成视频在线观看免费观看| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 日本a在线网址| 久久久精品94久久精品| 久久久久久人人人人人| 国产精品av久久久久免费| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩精品网址| 十八禁网站免费在线| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看舔阴道视频| 国产又爽黄色视频| 桃红色精品国产亚洲av| 黄色毛片三级朝国网站| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 亚洲精华国产精华精| 久久99热这里只频精品6学生| 日本五十路高清| 夜夜爽天天搞| 国产精品电影一区二区三区 | videosex国产| 亚洲精品国产一区二区精华液| 久久久精品区二区三区| 亚洲性夜色夜夜综合| 欧美中文综合在线视频| 午夜视频精品福利| 亚洲国产精品一区二区三区在线| 天堂8中文在线网| 免费看a级黄色片| 亚洲国产欧美一区二区综合| 国产av国产精品国产| 日韩成人在线观看一区二区三区| cao死你这个sao货| 日日摸夜夜添夜夜添小说| 国产真人三级小视频在线观看| 欧美日韩成人在线一区二区| 精品国产国语对白av| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人看| 在线 av 中文字幕| 19禁男女啪啪无遮挡网站| 可以免费在线观看a视频的电影网站| 欧美成人午夜精品| 一个人免费在线观看的高清视频| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 久久久久久久久久久久大奶| 精品少妇久久久久久888优播| 91成人精品电影| 纵有疾风起免费观看全集完整版| 亚洲第一欧美日韩一区二区三区 | av天堂在线播放| 国产欧美亚洲国产| 久9热在线精品视频| 最近最新中文字幕大全电影3 | 国产精品偷伦视频观看了| 亚洲国产中文字幕在线视频| 五月天丁香电影| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 久久午夜综合久久蜜桃| 国产视频一区二区在线看| 国产主播在线观看一区二区| a级毛片在线看网站| 色视频在线一区二区三区| 天天影视国产精品| 精品久久久久久电影网| 在线亚洲精品国产二区图片欧美| 中文欧美无线码| 成年人午夜在线观看视频| 久久久国产精品麻豆| 亚洲情色 制服丝袜| 午夜福利欧美成人| 成年女人毛片免费观看观看9 | 日韩精品免费视频一区二区三区| 狂野欧美激情性xxxx| 自线自在国产av| 欧美日韩视频精品一区| 黄网站色视频无遮挡免费观看| 男女下面插进去视频免费观看| 黑人操中国人逼视频| 国产欧美亚洲国产| 日日夜夜操网爽| 亚洲精品一二三| 亚洲国产中文字幕在线视频| 精品视频人人做人人爽| 一级片'在线观看视频| 国产野战对白在线观看| 久久这里只有精品19| 欧美日韩黄片免| 最新的欧美精品一区二区| 亚洲午夜精品一区,二区,三区| 久久中文字幕一级| 国产在视频线精品| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www | 日韩中文字幕视频在线看片| 国产在线一区二区三区精| 99精品久久久久人妻精品| 日韩一区二区三区影片| 国产黄色免费在线视频| 国产精品免费视频内射| 欧美精品人与动牲交sv欧美| 欧美日韩精品网址| 亚洲精华国产精华精| a级毛片黄视频| 叶爱在线成人免费视频播放| 国产在视频线精品| 国产伦人伦偷精品视频| 我要看黄色一级片免费的| 亚洲人成77777在线视频| 国产av国产精品国产| 亚洲久久久国产精品| 国产精品一区二区在线不卡| 久久天堂一区二区三区四区| 久久人妻av系列| 久热这里只有精品99| 日韩中文字幕视频在线看片| 国产xxxxx性猛交| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜精品国产一区二区电影| 日韩有码中文字幕| 黑人巨大精品欧美一区二区蜜桃| 一区二区av电影网| av线在线观看网站| 精品乱码久久久久久99久播| 精品第一国产精品| 午夜福利,免费看| 久久人妻av系列| 欧美激情 高清一区二区三区| 欧美国产精品一级二级三级| 国产av一区二区精品久久| 久久久国产一区二区| 国产一区有黄有色的免费视频| 性少妇av在线| 91成人精品电影| 成人国产av品久久久| 亚洲天堂av无毛| 999久久久精品免费观看国产| 亚洲人成电影观看| 动漫黄色视频在线观看| 波多野结衣一区麻豆| 超碰97精品在线观看| 国产精品成人在线| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 久久久久久久大尺度免费视频| 久久香蕉激情| 欧美+亚洲+日韩+国产| 久久香蕉激情| 国产精品影院久久| 久久久久视频综合| avwww免费| 国产激情久久老熟女| 日韩精品免费视频一区二区三区| 在线播放国产精品三级| 亚洲,欧美精品.| 久久影院123| 欧美午夜高清在线| 黄片播放在线免费| 亚洲精品国产精品久久久不卡| 黄色毛片三级朝国网站| 久久中文字幕人妻熟女| 国产老妇伦熟女老妇高清| 男女下面插进去视频免费观看| 下体分泌物呈黄色| 中国美女看黄片| 亚洲中文字幕日韩| 亚洲天堂av无毛| 狠狠精品人妻久久久久久综合| 免费在线观看黄色视频的| 999久久久国产精品视频| 国产成人一区二区三区免费视频网站| 新久久久久国产一级毛片| 夜夜夜夜夜久久久久| 国产野战对白在线观看| 2018国产大陆天天弄谢| 少妇被粗大的猛进出69影院| 久久热在线av| 久久久水蜜桃国产精品网| 窝窝影院91人妻| 国产精品.久久久| 丰满少妇做爰视频| 人成视频在线观看免费观看| 亚洲三区欧美一区| 99香蕉大伊视频| 久久人人97超碰香蕉20202| 国产成人一区二区三区免费视频网站| tube8黄色片| 老司机福利观看| 另类亚洲欧美激情| 国产99久久九九免费精品| 90打野战视频偷拍视频| 老司机深夜福利视频在线观看| 亚洲七黄色美女视频| 精品福利观看| 青草久久国产| 老熟女久久久| 高清av免费在线| 80岁老熟妇乱子伦牲交| 亚洲精品久久午夜乱码| 亚洲精品国产精品久久久不卡| 亚洲午夜理论影院| 狠狠狠狠99中文字幕| 欧美精品亚洲一区二区| 欧美成狂野欧美在线观看| 国产成人欧美| 国产麻豆69| 无遮挡黄片免费观看| 99国产精品免费福利视频| 国产日韩一区二区三区精品不卡| 午夜福利,免费看| 啦啦啦视频在线资源免费观看| 免费高清在线观看日韩| 老司机亚洲免费影院| 母亲3免费完整高清在线观看| 亚洲欧美日韩另类电影网站| 日本av免费视频播放| 久久久久久久大尺度免费视频| 免费不卡黄色视频| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 国产淫语在线视频| 黑人欧美特级aaaaaa片| 五月开心婷婷网| 亚洲专区字幕在线| 亚洲精华国产精华精| 精品国产超薄肉色丝袜足j| av福利片在线| 丝袜喷水一区| 高清在线国产一区| 国产午夜精品久久久久久| 午夜福利影视在线免费观看| 久久精品成人免费网站| 三上悠亚av全集在线观看| 亚洲免费av在线视频| 在线观看免费视频网站a站| 两个人看的免费小视频| 大陆偷拍与自拍| a级毛片黄视频| 如日韩欧美国产精品一区二区三区| 久久人人97超碰香蕉20202| 精品欧美一区二区三区在线| 亚洲男人天堂网一区| 亚洲午夜理论影院| 久久天躁狠狠躁夜夜2o2o| 国产高清视频在线播放一区| 国产97色在线日韩免费| 国产一卡二卡三卡精品| 精品国产国语对白av| 日韩人妻精品一区2区三区| 欧美亚洲日本最大视频资源| 日韩大片免费观看网站| 欧美激情 高清一区二区三区| 大片免费播放器 马上看| 男女下面插进去视频免费观看| 999久久久精品免费观看国产| av片东京热男人的天堂| kizo精华| 热re99久久国产66热| av网站免费在线观看视频| 超碰97精品在线观看| 免费人妻精品一区二区三区视频| 女人高潮潮喷娇喘18禁视频| 免费av中文字幕在线| 日本黄色日本黄色录像| 亚洲自偷自拍图片 自拍| 99riav亚洲国产免费| 嫩草影视91久久| www.精华液| 国产在线视频一区二区| 老司机深夜福利视频在线观看| 黄片大片在线免费观看| av天堂在线播放| 久久久久精品国产欧美久久久| 国产免费现黄频在线看| 九色亚洲精品在线播放| 国产国语露脸激情在线看| 天天影视国产精品| 看免费av毛片| 欧美在线一区亚洲| 日韩一区二区三区影片| 成年人黄色毛片网站| 狠狠精品人妻久久久久久综合| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| 一本色道久久久久久精品综合| 国产精品国产av在线观看| 久久国产亚洲av麻豆专区| 性色av乱码一区二区三区2| 另类精品久久| 欧美人与性动交α欧美精品济南到| 三上悠亚av全集在线观看| 一区二区av电影网| 大型黄色视频在线免费观看| 亚洲天堂av无毛| 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 亚洲av成人不卡在线观看播放网| 18禁裸乳无遮挡动漫免费视频| 在线观看免费视频网站a站| 国产高清videossex| 国产一卡二卡三卡精品| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三| 久久影院123| 啦啦啦 在线观看视频| 在线观看舔阴道视频| 69av精品久久久久久 | av有码第一页| 欧美大码av| 热99久久久久精品小说推荐| 亚洲精品国产色婷婷电影| 人妻一区二区av| 国精品久久久久久国模美| 亚洲人成77777在线视频| 欧美 亚洲 国产 日韩一| 69av精品久久久久久 | 亚洲黑人精品在线| 一二三四在线观看免费中文在| 亚洲黑人精品在线| 亚洲va日本ⅴa欧美va伊人久久| 天堂中文最新版在线下载| 欧美激情久久久久久爽电影 | 中国美女看黄片| 激情在线观看视频在线高清 | 亚洲五月婷婷丁香| 不卡一级毛片| 无人区码免费观看不卡 | 在线av久久热| h视频一区二区三区| 亚洲精品国产色婷婷电影| 高清视频免费观看一区二区| 91精品国产国语对白视频| 国产有黄有色有爽视频| 国产精品国产高清国产av | 久久久久视频综合| 亚洲视频免费观看视频| 久久国产精品人妻蜜桃| 国产亚洲一区二区精品| 大香蕉久久网| 久久人人爽av亚洲精品天堂| 国产精品一区二区在线观看99| 欧美 亚洲 国产 日韩一| 电影成人av| 热99国产精品久久久久久7| 亚洲黑人精品在线| 亚洲精品在线美女| 波多野结衣一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲免费av在线视频| 久久毛片免费看一区二区三区| 脱女人内裤的视频| 亚洲性夜色夜夜综合| 欧美激情久久久久久爽电影 | 丰满饥渴人妻一区二区三| 成人黄色视频免费在线看| 91成人精品电影| 正在播放国产对白刺激| 成人亚洲精品一区在线观看| 成年版毛片免费区| 大香蕉久久成人网| 久久久久久人人人人人| 国产区一区二久久| 两个人免费观看高清视频| 久久国产亚洲av麻豆专区| 天天操日日干夜夜撸| 亚洲精品国产区一区二| 成人国产av品久久久| 成人手机av| 亚洲精品国产区一区二| 亚洲国产av影院在线观看| 一区福利在线观看| 99精品在免费线老司机午夜| 国产亚洲精品久久久久5区| 变态另类成人亚洲欧美熟女 |