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

    2017年遼東灣夏季潮汐鋒位置變化的分析*

    2020-02-06 06:00:48張廣躍肖勁根張海彥李志成
    海洋與湖沼 2020年1期

    張廣躍 魏 皓 肖勁根 張海彥 李志成

    (天津大學海洋科學與技術(shù)學院 天津 300072)

    潮汐鋒是海洋鋒的一種,指在增溫季節(jié)由潮混合與層化共同作用形成的近岸混合區(qū)與遠岸層化區(qū)之間的過渡海域(Simpsonet al,1974)。鋒區(qū)的輻聚、輻散作用是鯷魚卵子和仔稚魚呈斑塊分布的原因(Weiet al,2003;魏皓等,2007;周鋒等,2008),潮汐鋒還會對葉綠素a和浮游生物的分布產(chǎn)生影響(Liuet al,2003;李洪波等,2004;Sunet al,2010;劉光興等,2015;Choiet al,2017),韋欽勝等(2011)進而提出了鋒區(qū)生態(tài)系統(tǒng)的觀點。因此,由表層風摩擦作用與底層潮混合作用而使上下混合層在淺海相互貫通而形成的潮汐鋒,對動力海洋學和海洋生態(tài)學研究都具有重要意義(Simpsonet al,1978;趙保仁,1987;趙保仁等,1993;魏皓等,1993;Liuet al,2003;唐啟升等,2005)。此外,渤海底部低氧區(qū)的范圍與鋒面所圍冷水團基本一致(張華等,2016),研究渤海夏季潮汐鋒變化,對于探討低氧形成及變化機制同樣具有重要意義。

    趙保仁(1985)計算了黃海和東海北部Simpson-Hunter層化參數(shù)KSH=log(H/Umax3)(Simpsonet al,1974),發(fā)現(xiàn)KSH=1.8等值線走向與整個底層黃海冷水團的邊緣一致,其中H為水深,Umax為潮流最大特征流速。魏皓等(1993)根據(jù)Stigebrandt公式(Stigebrandt,1981)對夏季南黃海底層潮汐鋒位置進行了較為準確的計算,認為底層潮流和風對底層鋒的位置起決定性作用。Bi等(1993)對黃海西部潮汐鋒的形成和演化過程進行了研究并探討了潮混合、熱通量和風對潮汐鋒的影響,認為鋒區(qū)和深海層化區(qū)的溫度結(jié)構(gòu)受到潮混合作用和平流的共同影響,海流和鋒面會對風應(yīng)力的變化作出滯后響應(yīng),且鋒面的強度和位置還會受到Ekman輸運的影響。楊殿榮等(1994)在Stigebrandt公式基礎(chǔ)上加入浮力平流項,進一步提高了南黃海海洋鋒的預(yù)報質(zhì)量。周鋒等(2008)綜合調(diào)查資料和MODIS/Terra衛(wèi)星遙感海表面溫度(sea surface temperature,SST)資料,并利用HAMSOM(Hamburg Shelf Ocean Model)模式計算Simpson-Hunter參數(shù)分布,對南黃海西北部夏季潮鋒的不連續(xù)分布現(xiàn)象做了分析和討論,認為潮流較弱、地形相對平緩和層化作用較強是潮鋒不連續(xù)分布的主要成因,而夏季風場對鋒面分布的影響基本可以忽略。Lü等(2010)基于POM(Princeton Ocean Model)模式建立了波浪-潮汐-環(huán)流耦合數(shù)值模型,用以研究潮汐混合對夏季黃海環(huán)流造成的影響,結(jié)果顯示黃海冷水團被典型的潮汐鋒所包圍,而較大的斜壓梯度使鋒區(qū)成為了上升流最活躍的區(qū)域,同時黃海冷斑的形成主要是由鋒區(qū)上升流和潮混合造成的。

    近年來國內(nèi)學者對于潮汐鋒的研究大都集中在黃、東海及南海等區(qū)域(朱鳳芹等,2014;程雪麗等,2017;王謙等,2019),針對渤海潮汐鋒的研究比較少。趙保仁等(2001)對渤海夏季無風狀態(tài)下潮汐鋒的分布特征進行過較為系統(tǒng)的研究,結(jié)果表明渤海海峽、遼東灣東部近海、遼東灣西部六股河口近海及渤海灣口的南、北近岸水域分布著潮汐鋒。劉浩等(2007)利用POM模式對渤海潮汐鋒及層化結(jié)構(gòu)的季節(jié)變化進行了研究,通過與Loder參數(shù)(KL=log(h4/U3),h、U分別為水深和潮流速度)計算結(jié)果對比,認為采用表底溫差指標研究渤海的層化結(jié)構(gòu)特征和潮汐鋒是可靠的,同時模擬結(jié)果表明渤海層化形成于4月份,之后隨著太陽輻射的增強以及垂向混合的減弱,層化區(qū)面積以及穩(wěn)定性都進一步增強并在7月份達到頂峰,隨著層化區(qū)域范圍的不斷擴大和層化強度的增加,潮汐鋒也隨之向淺水區(qū)移動。

    綜上所述,前人對潮汐鋒的研究,大都集中在潮汐鋒的季節(jié)變化、潮汐鋒的形成機制等方面,對夏季期間潮汐鋒位置的變化討論很少。根據(jù)宋貴生等人2017、2018年秦皇島外海低氧調(diào)查結(jié)果得知,夏季潮汐鋒恰是影響秦皇島外海低氧維持和分布的重要原因(個人交流)。因此本文綜合利用現(xiàn)場觀測數(shù)據(jù)和ROMS(Regional Ocean Modeling System)模擬結(jié)果,并基于Stigebrandt公式對2017年夏季遼東灣潮汐鋒位置進行診斷計算,探討風、凈熱通量和浮力平流對潮汐鋒位置的影響,為研究遼東灣生物聚集和低氧分布機制提供依據(jù)。

    1 數(shù)據(jù)和方法

    1.1 數(shù)據(jù)介紹

    現(xiàn)場觀測數(shù)據(jù)來自國家重點研發(fā)計劃“海洋環(huán)境安全保障”重點專項于2017年7月11—20日和2017年8月11—20日兩個時間段內(nèi)對遼東灣東部海域(120.5°—121.5°E,39.3°—40.3°N)(圖)進行的逐月水文調(diào)查。在7月份由于天氣等原因?qū)е挛磳1、G1、G4站位進行采樣。溫度、鹽度、密度等觀測數(shù)據(jù),主要利用加拿大RBR公司生產(chǎn)的CTD(Conductivity-Temperature-Depth profiler)獲得,溫度精度 0.002°C,電導(dǎo)率精度0.003mS/cm,設(shè)定采樣頻率6Hz。利用Matlab-2014a對CTD獲得的數(shù)據(jù)資料進行質(zhì)量控制處理、剔除異常值并進行0.2m平均后進行繪圖分析。

    1.2 模式配置

    本文使用的數(shù)值模式ROMS,是一個采用自由表面、靜力假定、σ坐標、曲線正交網(wǎng)格的三維原始方程模式(Shchepetkinet al,2005)。目前該模式已被廣泛應(yīng)用于中國近海水動力模擬研究中(Bianet al,2013;Zhouet al,2017)。本文ROMS模式區(qū)域覆蓋了渤、黃海及東海部分區(qū)域(117.01°—131.66°E,29.04°—42.09°N),地形采用成均館大學提供的分辨率為1’×1’的地形數(shù)據(jù),模式最小水深為5m。模式水平分辨率約為1/24°×1/24°,且在渤海海區(qū)進行了加密(分辨率約為3km),垂向分為30層,在表底進行了加密。模式垂向采用MY2.5(Mellor-Yamada 2.5)湍混合參數(shù)化方案(Melloret al,1982),二次底拖曳系數(shù)為0.0015。風場、海面熱通量等大氣強迫數(shù)據(jù)來自歐洲中期天氣預(yù)報中心(European Centre for Medium-range Weather Forecast,ECMWF,https://www.ecmwf.int/)的 ERAInterim再分析數(shù)據(jù)資料。模式有兩個開邊界,南側(cè)開邊界位于杭州灣南側(cè),東側(cè)開邊界位于朝鮮海峽附近,開邊界處的溫度、鹽度、水位及流速數(shù)據(jù)來自HYCOM(Hybrid Coordinate Ocean Model)。開邊界處添加M2、S2、N2、K2、K1、O1、P1和Q1,共8個分潮,潮汐數(shù)據(jù)來自于OTPS(OSU Tidal Prediction Software)。模式初始場的溫度及鹽度來自WOA2013(World Ocean Atlas 2013)一月份的氣候態(tài)數(shù)據(jù),水位及流速的初始值設(shè)為零。模式共包括9條流入渤、黃、東海的主要河流,分別為遼河、灤河、海河、黃河、鴨綠江、漢江、淮河、長江和錢塘江。河流流量數(shù)據(jù)來源于中國河流泥沙公報及相關(guān)文獻(Zhang,1996;Liuet al,2009;Tonget al,2015)。模式運行時間為2006年1月1日—2017年12月31日,本文使用2017年的模擬結(jié)果對遼東灣潮汐鋒位置的變化進行研究。從模式的渤海海表溫度和底溫月平均序列來看(圖略),模式在研究年份已達到穩(wěn)定。

    1.3 模式校驗

    圖2 夏季(6—8月)渤海垂向平均流場結(jié)構(gòu)Fig.2 Vertical structure of circulations of Bohai Sea in summer average(June-August)

    由圖2所示夏季渤海垂向平均流場結(jié)構(gòu)可知,夏季期間渤海大部分區(qū)域深度平均流速小于10cm/s,渤海中部存在反氣旋式環(huán)流,而遼東灣存在氣旋式環(huán)流,這與Zhou等(2017)文章中模擬的渤海流場結(jié)構(gòu)基本一致。夏季渤海海峽南部存在出流區(qū),北部存在入流區(qū),整體為“北進南出”的海流結(jié)構(gòu),這與黃大吉等(1998)、魏澤勛等(2003)、冀承振等(2019)研究結(jié)果比較一致。

    選取圖1所示的D、E、F斷面觀測站點得到的溫度數(shù)據(jù)與ROMS模擬的溫度結(jié)果進行垂向?qū)Ρ?圖3)。雖然模式的整體溫度偏低,但月份之間的溫度變化趨勢模擬較好。由2017年遼東灣調(diào)查數(shù)據(jù)(圖3中柱條)可以看出,7月份D、E、F斷面實測溫躍層大部分位于6—8m水深處,在D3站點東側(cè)、E3站點附近以及F4站點附近水體充分混合,潮汐鋒鋒線處水深為18—26m(圖3a、3c和3e)。8月份實測D、E、F斷面水溫整體升高約1°C,E斷面溫躍層強度減弱、近岸混合區(qū)范圍有所擴大,D、F斷面溫躍層變動不太明顯,8月份實測D3站位西側(cè)、E3站位西側(cè)、F3站位附近存在潮汐鋒,鋒線處水深為25—28m,鋒線位置向深水區(qū)移動(圖3b、3d和3f)。通過模擬結(jié)果(圖3背景顏色)可以看出,7月份D、E、F斷面躍層基本位于12.5m水深處,底層冷水團范圍較大,邊緣位于25m等深線附近,鋒面位于冷水團的邊緣(圖3a、3c和3e)。8月份D、E、F斷面模擬的溫度較7月份高約0.8°C,這與實測溫度變化趨勢一致,溫躍層深度變化不大,但冷水團的范圍向深水區(qū)收縮,其邊緣位于28m等深線附近,潮汐鋒的位置也隨之向深水區(qū)移動(圖3b、3d和3f)。

    選取遼東灣東部D、E、F斷面存在溫度躍層的站點,對站點位置的溫躍層強度和溫躍層深度進行了數(shù)據(jù)校驗。7、8月份溫躍層強度變化不大,模式和觀測的相關(guān)指數(shù)R2為0.84,表明模式能夠較好的模擬出躍層的強度(圖4a)。7、8月份的溫躍層深度沒有明顯變化,模式和觀測的相關(guān)指數(shù)R2為0.86,因此模式能夠較好地模擬出躍層的位置(圖4b)。

    圖3 2017年7、8月份遼東灣東部D(a、b)、E(c、d)、F(e、f)斷面溫度分布Fig.3 Temperature distribution of Section D(a,b),E(c,d),and F(e,f)in the eastern Liaodong Bay in July and August 2017

    圖4 2017年7(藍點)、8(紅點)月份遼東灣東部觀測站點溫度躍層相關(guān)模擬結(jié)果與觀測數(shù)據(jù)Fig.4 Simulation and observation of the thermocline at the sampling stations of Liaodong Bay in July(blue point)and August(red point)2017

    綜上所敘,模型能夠較好的模擬出渤海夏季流場結(jié)構(gòu)、底部冷水團位置、溫躍層強度和溫躍層深度,較為準確的反映出近海混合區(qū)與遠海層化區(qū)的位置,在一定程度上反映出實測潮汐鋒的位置變化,可以使用模擬結(jié)果進行潮汐鋒的進一步分析。

    1.4 潮汐鋒位置的估算方法

    參照劉浩等(2007)對渤海潮汐鋒的研究并結(jié)合觀測數(shù)據(jù)和模擬結(jié)果,本文將表底溫差0.2°C等值線作為實際潮汐鋒位置。由于Simpson-Hunter層化參數(shù)KSH=log(H/Umax3)(Simpsonet al,1974)僅考慮了水深和潮流對潮汐鋒位置的影響,未考慮其他影響潮汐鋒的要素,不能夠準確反映出潮汐鋒位置的月度變化。本文參考魏皓等(1993)和楊殿榮等(1994)等引入浮力平流效應(yīng)并改進的Stigebrandt公式來計算潮汐鋒的位置:

    其中,Hfront為潮汐鋒所處的水深,A為底層潮流對鋒的作用項,B為表層風場對鋒的作用項,C代表凈熱通量對鋒的作用項,M代表浮力平流對鋒的作用項。

    A的表達式如下:

    式(2)中,mob為下層挾卷參量,可取為1.12(Stigebrandt,1981;魏皓等,1993),而為底層潮流流速,海底阻力系數(shù)取值與ROMS模式設(shè)置一致。

    B的表達式如下:

    式(3)中,mow上層挾卷參量,可取為1.25(Garwood,1977;朱建榮等,1993;魏皓等,1993),而Uw為表層風速,海水平均密度空氣密度海面阻力系數(shù)。

    C的表達式如下:

    式(4)中,Q為凈熱通量,9.8g=m/s2,熱膨脹系數(shù)系統(tǒng)比熱。

    M的表達式如下:

    式(5)中,ρ1、ρ2分別為上層海水密度,h1、h2分別為上、下層海水的厚度,u1、v1分別為上層海水東分量流速和北分量流速,u2、v2分別為下層海水東分量流速和北分量流速。

    為了減小噪音和誤差,參照楊殿榮等(1994)對M進行處理,將式(5)做如下轉(zhuǎn)換:

    對近岸水體而言,M可表示為:

    當不考慮M項時,Stigebrandt公式為:

    2 結(jié)果分析

    2.1 遼東灣水團的變化

    潮汐鋒作為水團的邊界,是近岸混合區(qū)與遠海層化區(qū)的過渡帶。由遼東灣底層溫度和表底溫差的分布變化,可知2017年6、7、8月份,遼東灣整體處于升溫階段(圖5)。6月份遼東灣底層水溫相對較低,大部分區(qū)域底層水溫低于17°C,僅在遼東灣北部底層水溫為18°C左右(圖5a)。7月份底層溫度整體升高,低于17°C的冷水范圍向西南方向收縮至40°N,40°N以北的大部分區(qū)域水溫在20°C以上(圖5b)。8月份底層溫度進一步升高,底層冷水的范圍進一步縮小,大部分區(qū)域溫度升高至22°C以上(圖5c)。總的來說冷水團的范圍不斷向深水區(qū)縮小,底層溫度在6—7月升溫幅度較大,平均升幅在5°C左右,7—8月升溫幅度較小,平均升幅在2°C左右。

    圖5 2017年夏季遼東灣底層溫度(a、b、c)和表底溫差(d、e、f)月度分布Fig.5 Monthly distribution of bottom temperature(a,b and c)and temperature difference between the surface and bottom(d,e,and f)in Liaodong Bay during summer 2017

    通過2017年6、7、8月份表底溫差(圖5d、5e和5f)分布變化可知,6月份表底溫差最大值不超過4°C,說明6月份層化較弱,而小于0.2°C溫差等值線的區(qū)域主要為遼東灣沿岸海域,混合區(qū)的范圍較小。7月份表底溫差最大值7°C,說明7月份遼東灣層化明顯加強,而0.2°C溫差等值線在遼東灣北部位置基本保持不變,但在東、西兩岸,0.2°C溫差等值線略微向深水區(qū)移動。8月份,0.2°C溫差等值線明顯向深水區(qū)移動,遼東灣東北部0.2°C溫差等值線位置變化最大,向深水區(qū)移動大約18km。整體來看,遼東灣近岸混合區(qū)范圍不斷向深水區(qū)擴展,遼東灣中部及西南部層化不斷加強,但層化區(qū)的范圍在減小。通過計算Simpson-Hunter層化參數(shù)得到的潮汐鋒位置(圖5洋紅線)在6、7、8月份基本保持不變,這說明僅考慮潮流和地形因素計算潮汐鋒位置的Simpson-Hunter層化參數(shù)并不能夠很好的反映遼東灣夏季期間潮汐鋒的位置變化,尤其是在8月份。

    2.2 潮汐鋒位置的估算

    圖6 2017年夏季遼東灣底層潮流作用項A(a、b、c)、表層風場作用項B(d、e、f)、凈熱通量C(g、h、i)和浮力平流效應(yīng)項M(j、k、l)月度分布Fig.6 Monthly distribution of bottom tidal mixing term A(a,b,and c),wind friction term B(d,e,and f),net heat fluxes term C(g,h,and i)and buoyancy advection term M(j,k,and l)in Liaodong Bay during summer 2017

    由式(1)可知,表層風場、底層潮流和凈熱通量對潮汐鋒的位置有著重要作用,而浮力平流對潮汐鋒的位置起到一定的調(diào)節(jié)作用。如圖6所示,2017年6、7、8月份底層潮流作用項A基本保持不變,數(shù)值量級為10-6m3/s3,整體呈現(xiàn)出遼東灣東岸高西岸低的分布態(tài)勢(圖6a、6b和6c)。2017年6、7、8月份遼東灣風場作用項B,數(shù)值量級為10-7—10-6m3/s3(圖6d、6e和6f)。6月份風場作用項B高值區(qū)位于遼東灣中部和南部,量值約為0.7×10-7m3/s3(圖6d)。7月份風場作用項B高值區(qū)位于遼東灣西北部,量值約為0.71×10-7m3/s3(圖6e)。8月份風場作用項B明顯減小,平均約為0.18×10-7m3/s3(圖6f)。整體來說,6、7月份表層風場作用變化不大,6月份遼東灣中南部B項值略高,而7月份遼東灣北部B項值略高,8月份B項值明顯減小。2017年6、7、8月份遼東灣凈熱通量項C,整體呈現(xiàn)出減小的變化趨勢,且均為自南向北逐漸降低的分布態(tài)勢(圖6g、6h和6i)。6、7月份,遼東灣絕大部分區(qū)域凈熱通量作用項C大于1×10-7m3/s3,7月份較6月份C項略微減小約0.1×10-7m3/s3(圖6g、圖6h)。8月份遼東灣凈熱通量作用項C均小于0.8×10-7m3/s3,8月份較 7月份C項大幅減少約 0.8×10-7m3/s3(圖6i)。2017年6、7、8月份遼東灣浮力平流效應(yīng)項M(圖6j、圖6k和圖6l)均大于0。6、7月份遼東灣大部分區(qū)域M項數(shù)值約為1.4×10-7m3/s3(圖6j、圖6k),且在遼東灣北部和東、西兩岸離岸較近的地方,存在M項的極大值區(qū),8月份遼東灣大部分區(qū)域M項數(shù)值小于1×10-7m3/s3(圖6l),遼東灣北部和東、西兩岸存在的M項極大值區(qū)向深水區(qū)移動,尤以遼東灣東北部移動最為明顯。

    通過圖7對潮汐鋒位置的診斷可知,2017年6、7月份,遼東灣潮汐鋒位置變化不大,鋒線主要位于遼東灣北部淺海和東、西沿岸,呈“幾”字形分布。7月份遼東灣東、西兩岸鋒線略微向深水區(qū)移動,8月份潮汐鋒位置向深水區(qū)大幅移動,尤其以遼東灣東北部鋒線向深水區(qū)移動最為明顯,最大移動距離約為20km。通過圖7加入平流項(紅線)和未加入平流項(藍線)得到的潮汐鋒線的對比可知,2017年6、7月份兩種情況下計算得到的潮汐鋒線與模擬結(jié)果都比較吻合,但在8月份,未加入平流項計算得到的潮汐鋒位置與模擬結(jié)果相比,過于偏向深水區(qū),計算結(jié)果誤差較大。而加入平流項得到的潮汐鋒位置與模擬結(jié)果較為接近,最大偏差距離不超過12km,能夠較為準確的計算出潮汐鋒的位置變化情況。這可能是因為2017年6、7月份遼東灣浮力平流對潮汐鋒的作用不夠明顯,8月份在凈熱通量和風的作用都大幅減少的情況下,浮力平流的作用使計算得到的潮汐鋒位置向淺水區(qū)偏移的效果更加明顯。

    圖7 2017年夏季遼東灣,式(1)得到的鋒線(紅線)和式(8)得到的鋒線(藍線)與模型鋒線(黑線)位置對比Fig.7 Comparison in the position of the tidal front obtained by equation(1)(red line),equation(8)(blue line),and ROMS(black line)in Liaodong Bay during summer 2017

    3 討論

    3.1 風場和流場的變化

    由圖8可知,2017年6、7月份渤海地區(qū)主要以南風為主,風速較大,遼東灣平均風速為5m/s左右,且7月份遼東灣北部風速要大于6月份。8月份,整個渤海風場轉(zhuǎn)向,主要以東風為主,風速約減小為2.5m/s,遼東灣風場大小為東北低西南高的分布態(tài)勢。Bi等(1993)對黃海西部34oN斷面潮汐峰進行過理論研究,其結(jié)果顯示南風誘發(fā)的Ekman輸送離岸,可使該區(qū)域的潮汐鋒位置向深水區(qū)移動,北風誘發(fā)的Ekman輸送向岸,可使該區(qū)域的潮汐鋒位置向淺海移動。依據(jù)以上理論結(jié)果,6月份和7月份,遼東灣盛行南風,此時Ekman輸送基本沿著鋒線,對遼東灣灣口的潮汐鋒影響較小。而遼東灣西側(cè)Ekman輸送為離岸,遼東灣東側(cè)Ekman輸送為向岸,因此遼東灣西側(cè)潮汐鋒向深水移動,而遼東灣東側(cè)潮汐鋒向淺水移動,但圖7a和圖7b中遼東灣東側(cè)和西側(cè)潮汐鋒變動范圍不顯著,8月份由于風速轉(zhuǎn)向和減弱,Ekman效應(yīng)更加不顯著。同時由式(1)可知,當?shù)讓映绷骰颈3植蛔儠r,表層風場增大會使潮汐鋒向深水區(qū)移動,風場減小會使潮汐鋒向淺水區(qū)移動,而凈熱通量項C作為分母,當凈熱通量增加時,會使潮汐鋒位置向淺水區(qū)移動,凈熱通量降低會使潮汐鋒向深水區(qū)移動。所以結(jié)合圖6可知,當?shù)讓映绷髯饔庙桝基本保持不變時,2017年7月份遼東灣凈熱通量作用項C較6月份略微減小,風速7月份較6月份略微增大,風和凈熱通量的共同作用使得遼東灣東、西兩岸7月份鋒線位置較6月份稍微向深水區(qū)移動。而8月份凈熱通量項C較7月份大幅降低約50%,表層風場作用項B也大幅減小70%,而潮汐鋒位置較7月份整體向深水區(qū)移動約8km。由式(1)可知,底層潮流作用A和表層風場作用B同時作為分子項,由于A項數(shù)值較B項數(shù)值要大很多(圖6),在A項基本不變的情況下,B項的變化對于潮汐鋒的影響并不敏感,而凈熱通量作用項C作為分母,對潮汐鋒的位置變化有著更為顯著的影響,因此可知8月份凈熱通量減弱是造成遼東灣潮汐鋒向深水區(qū)移動的主要原因。

    圖8 2017年夏季渤海表層風場(a、b、c)月度分布和表底溫差與環(huán)流(d、e、f)的月度分布Fig.8 Monthly distribution of sea surface wind field(a,b,c),surface-bottom temperature difference with circulation(d,e,f)in Bohai Sea during summer 2017

    而由圖8d、8e和8f可知,渤海中部夏季整體為反氣旋環(huán)流,而在遼東灣中部存在一個逆時針流環(huán),逆時針流環(huán)的強度、范圍和位置會隨時間有所變化,這一模擬結(jié)果與畢聰聰(2013)、Zhou等(2017)等利用ROMS模式模擬的渤海環(huán)流結(jié)果基本一致。2017年夏季渤海表底溫差的最大值基本上位于逆時針流環(huán)的中心位置和順時針與逆時針流環(huán)交接的地帶。遼東灣范圍內(nèi)浮力平流項M(圖6j、6k和6l)的極大值區(qū)基本分布在逆時針流環(huán)邊緣,根據(jù)楊殿榮等(1994)的結(jié)論,當M>0時,這種熱鹽環(huán)流導(dǎo)致的浮力平流使得潮汐鋒位置向淺水區(qū)偏移,尤其在8月份浮力平流起到了較好的修正效果。

    3.2 凈熱通量的變化分析

    由3.1分析可知,8月份遼東灣潮汐鋒位置向深水區(qū)移動,這主要是由凈熱通量減弱造成的。凈熱通量的變化主要由海表凈短波輻射、海表凈長波輻射、潛熱通量和感熱通量所決定,各熱收支項夏季變化如圖9所示,2017年6、7、8月份遼東灣海表凈短波輻射逐漸降低,其中7—8月海表凈短波輻射降幅約60.0W/m2。2017年6、7、8月份遼東灣海表凈長波輻射均為負值,屬于失熱項,其絕對值先減少后略微增加,而潛熱通量在6、7月份基本保持不變,8月份潛熱通量變化幅度約為30.0W/m2。2017年6、7、8月,3個月之間的感熱通量變化幅度基本保持在10.0W/m2左右,其中6、7月份感熱通量為正值,8月份為負值。

    通過表1可以得出,2017年7月份遼東灣海面凈熱通量較6月份略微減小10.0W/m2,但8月份海面凈熱通量較7月份大幅減少約109.0W/m2,可見這是導(dǎo)致潮汐鋒在8月份向深水區(qū)大幅移動的主要因素。8月份凈熱通量的降低主要是由于海表凈短波輻射和潛熱通量的減少造成,其中海表凈短波輻射的減少是導(dǎo)致凈熱通量減少的最主要因素,約占凈熱通量降幅的55%。

    圖9 2017年夏季遼東灣海表凈短波輻射(a、b、c)、海表凈長波輻射(d、e、f)、潛熱通量(g、h、i)和感熱通量(j、k、l)月度分布Fig.9 Monthly distribution of net short-wave radiation(a,b,and c),net long-wave radiation(d,e,and f),latent heat fluxes(g,h,and i)and sensible heat fluxes(j,k,and l)in Liaodong Bay during summer 2017

    表1 2017年夏季遼東灣海域凈熱通量及熱收支主要因子項變化幅度Tab.1 The changes of net heat flux and heat budget in Liaodong Bay area during summer 2017

    4 結(jié)論

    通過2017年7、8月份遼東灣東部海域現(xiàn)場實測CTD資料,對ROMS模擬結(jié)果進行了驗證,結(jié)果表明模式可以用于遼東灣潮汐鋒位置變化的分析?;赗OMS模擬的溫度場,選取表底0.2°C溫差等值線與觀測數(shù)據(jù)得到的潮汐鋒的位置進行對比,發(fā)現(xiàn)二者變化趨勢比較一致,因而表底0.2°C溫差等值線可作為遼東灣潮汐鋒位置的代用指標。

    利用ROMS模擬結(jié)果,基于Stigebrandt公式分兩種情況(加入平流項和不加入平流項)對潮汐鋒的位置進行計算,得到如下結(jié)論:2017年夏季遼東灣潮汐鋒位于遼東灣北部和東、西沿岸淺水區(qū)域,呈“幾”字形分布。6月份到8月份,雖然遼東灣平均風速由較強南風逐漸減弱并轉(zhuǎn)為弱東風,但是Ekman作用總體上不顯著。6、7月份鋒線位置離岸較近,位置變化不大,浮力平流作用效果相對不明顯。8月份遼東灣表層風場和海表凈熱通量均減弱,潮汐鋒由沿岸淺水區(qū)向離岸深水區(qū)收縮。8月份海表凈熱通量的減弱可能是鋒線向深水區(qū)收縮的主要影響因素,進一步研究發(fā)現(xiàn)海表凈短波輻射的減少是引起8月份海表凈熱通量減弱的主要原因??紤]浮力平流后,與未加浮力平流的計算結(jié)果相比,鋒線位置向淺水區(qū)偏移,與ROMS模擬的鋒線結(jié)果更加接近,可見浮力平流對8月份遼東灣潮汐鋒位置的準確估算起到了較好的修正作用。

    人妻制服诱惑在线中文字幕| 久久这里只有精品中国| 日本五十路高清| 国内少妇人妻偷人精品xxx网站| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久久久免费av| 伦精品一区二区三区| 中文资源天堂在线| 国产成人午夜福利电影在线观看| 中国美白少妇内射xxxbb| 久久久久国产网址| 夜夜看夜夜爽夜夜摸| 九九在线视频观看精品| 最近手机中文字幕大全| 1024手机看黄色片| 三级毛片av免费| 美女cb高潮喷水在线观看| 国产老妇伦熟女老妇高清| 97热精品久久久久久| 久久99热这里只频精品6学生 | 男女视频在线观看网站免费| 91在线精品国自产拍蜜月| 搡老妇女老女人老熟妇| 青春草亚洲视频在线观看| 男女边吃奶边做爰视频| 国产一区二区在线av高清观看| 有码 亚洲区| 丰满人妻一区二区三区视频av| 国产不卡一卡二| 99热网站在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久久国产成人免费| 欧美成人a在线观看| 欧美xxxx黑人xx丫x性爽| 一夜夜www| АⅤ资源中文在线天堂| 色综合亚洲欧美另类图片| 菩萨蛮人人尽说江南好唐韦庄 | 变态另类丝袜制服| 五月玫瑰六月丁香| 国产高清国产精品国产三级 | 又粗又爽又猛毛片免费看| 高清av免费在线| 丰满少妇做爰视频| 天堂影院成人在线观看| 永久网站在线| 亚洲欧洲日产国产| 欧美丝袜亚洲另类| 少妇熟女aⅴ在线视频| 中文欧美无线码| 99九九线精品视频在线观看视频| 老司机福利观看| 天天一区二区日本电影三级| 亚洲综合色惰| 久久久国产成人免费| 精品国内亚洲2022精品成人| 国产精品国产三级专区第一集| 色综合站精品国产| 熟女人妻精品中文字幕| 国产伦一二天堂av在线观看| 午夜a级毛片| 嫩草影院入口| 啦啦啦观看免费观看视频高清| 亚洲精品aⅴ在线观看| АⅤ资源中文在线天堂| 青春草国产在线视频| 久久久久免费精品人妻一区二区| 欧美bdsm另类| 一级毛片我不卡| 最近最新中文字幕大全电影3| 99视频精品全部免费 在线| 国产精品久久久久久精品电影| 国产精品久久视频播放| 中文字幕久久专区| 国产极品天堂在线| 日本免费a在线| 国产精品伦人一区二区| 亚洲成人av在线免费| 亚洲av男天堂| 久99久视频精品免费| 国产在线男女| 久久久午夜欧美精品| 麻豆国产97在线/欧美| 91久久精品国产一区二区成人| 1024手机看黄色片| 日本黄色片子视频| 久热久热在线精品观看| 别揉我奶头 嗯啊视频| 久久精品国产鲁丝片午夜精品| 亚洲美女搞黄在线观看| 免费看日本二区| 乱码一卡2卡4卡精品| 亚洲人成网站在线观看播放| 91狼人影院| 日韩成人伦理影院| 在线天堂最新版资源| 亚洲成人中文字幕在线播放| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站高清观看| 免费看a级黄色片| 日本免费a在线| 中文字幕制服av| 六月丁香七月| 久久人妻av系列| 精品一区二区三区人妻视频| 免费观看人在逋| 狂野欧美白嫩少妇大欣赏| 韩国高清视频一区二区三区| 波多野结衣高清无吗| 午夜福利在线观看免费完整高清在| 欧美又色又爽又黄视频| 亚洲精品aⅴ在线观看| 精华霜和精华液先用哪个| 99久久精品热视频| 久久精品综合一区二区三区| 国产精品精品国产色婷婷| 午夜老司机福利剧场| a级毛色黄片| 成人亚洲精品av一区二区| 99热精品在线国产| 不卡视频在线观看欧美| 三级国产精品欧美在线观看| 少妇丰满av| 黑人高潮一二区| 国产精品一区二区在线观看99 | 国产欧美日韩精品一区二区| 国产高清不卡午夜福利| 午夜a级毛片| 国产成人a∨麻豆精品| 午夜福利成人在线免费观看| 少妇的逼水好多| 亚洲精品aⅴ在线观看| 亚洲18禁久久av| 插阴视频在线观看视频| 女的被弄到高潮叫床怎么办| 日韩欧美在线乱码| 日本免费在线观看一区| 夫妻性生交免费视频一级片| 久久久久精品久久久久真实原创| or卡值多少钱| 51国产日韩欧美| 欧美成人午夜免费资源| 亚洲性久久影院| 三级男女做爰猛烈吃奶摸视频| 亚洲成人精品中文字幕电影| 变态另类丝袜制服| 亚洲aⅴ乱码一区二区在线播放| 日本一二三区视频观看| 成人性生交大片免费视频hd| 嘟嘟电影网在线观看| 成人综合一区亚洲| 99久久九九国产精品国产免费| 又粗又爽又猛毛片免费看| 国产精品国产三级国产av玫瑰| 综合色丁香网| 最后的刺客免费高清国语| 人妻系列 视频| 国产高清三级在线| 伊人久久精品亚洲午夜| 波多野结衣高清无吗| 国产精品女同一区二区软件| 国产色爽女视频免费观看| 男人舔奶头视频| 国产黄片视频在线免费观看| 大又大粗又爽又黄少妇毛片口| 只有这里有精品99| 国内精品宾馆在线| 少妇高潮的动态图| 一级毛片电影观看 | 人妻夜夜爽99麻豆av| 深爱激情五月婷婷| 男女那种视频在线观看| 深夜a级毛片| 午夜老司机福利剧场| 麻豆一二三区av精品| 国产不卡一卡二| 日本欧美国产在线视频| 精品酒店卫生间| 久久精品影院6| 啦啦啦啦在线视频资源| 日韩欧美 国产精品| 久久99热这里只频精品6学生 | 在线观看一区二区三区| 亚洲国产色片| 人人妻人人澡人人爽人人夜夜 | av视频在线观看入口| 又黄又爽又刺激的免费视频.| 国产三级中文精品| 亚洲国产精品成人综合色| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 秋霞在线观看毛片| 成人亚洲欧美一区二区av| 日韩中字成人| 免费搜索国产男女视频| 亚洲av电影不卡..在线观看| 欧美+日韩+精品| 欧美日韩在线观看h| 免费大片18禁| 青春草国产在线视频| 亚洲成av人片在线播放无| 91午夜精品亚洲一区二区三区| 日日摸夜夜添夜夜爱| 中文字幕久久专区| 精品一区二区三区人妻视频| 免费观看a级毛片全部| 国产精品人妻久久久久久| 国产精品乱码一区二三区的特点| av福利片在线观看| 搡女人真爽免费视频火全软件| 一个人免费在线观看电影| 亚洲国产精品成人综合色| 久久草成人影院| 国产日韩欧美在线精品| 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 欧美精品国产亚洲| 欧美性感艳星| 亚洲成人久久爱视频| 啦啦啦观看免费观看视频高清| 亚洲色图av天堂| 日产精品乱码卡一卡2卡三| 久久婷婷人人爽人人干人人爱| 色播亚洲综合网| 中国美白少妇内射xxxbb| 亚洲国产精品专区欧美| 亚洲av不卡在线观看| 男的添女的下面高潮视频| or卡值多少钱| 久久久色成人| 色噜噜av男人的天堂激情| 国产精品一区二区三区四区久久| 欧美成人一区二区免费高清观看| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 日韩精品青青久久久久久| 长腿黑丝高跟| 免费观看人在逋| 亚洲不卡免费看| 全区人妻精品视频| 国产精品一区二区性色av| 欧美97在线视频| 亚洲在久久综合| 亚洲av.av天堂| 简卡轻食公司| 两性午夜刺激爽爽歪歪视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 成人鲁丝片一二三区免费| 人妻制服诱惑在线中文字幕| 九九久久精品国产亚洲av麻豆| 色播亚洲综合网| 日本一本二区三区精品| 久久久久久久久久久免费av| 亚洲欧美成人综合另类久久久 | 精品少妇黑人巨大在线播放 | 国产精品久久久久久av不卡| 内地一区二区视频在线| 欧美不卡视频在线免费观看| 久久久成人免费电影| av国产久精品久网站免费入址| 亚洲人成网站在线播| 国产又色又爽无遮挡免| 亚洲自拍偷在线| 91在线精品国自产拍蜜月| 成人亚洲精品av一区二区| 国产精品国产高清国产av| 久久这里只有精品中国| 亚洲五月天丁香| 男女那种视频在线观看| 精品一区二区三区人妻视频| 亚洲av福利一区| 乱系列少妇在线播放| 久久久a久久爽久久v久久| 亚洲真实伦在线观看| 亚洲五月天丁香| 中文字幕亚洲精品专区| h日本视频在线播放| 免费不卡的大黄色大毛片视频在线观看 | 成人毛片a级毛片在线播放| 色综合站精品国产| 国产一级毛片在线| 麻豆成人av视频| 国产三级在线视频| 久久久精品94久久精品| 国产人妻一区二区三区在| 国产日韩欧美在线精品| 久久久午夜欧美精品| 69人妻影院| 建设人人有责人人尽责人人享有的 | 久久精品综合一区二区三区| 午夜激情欧美在线| av专区在线播放| 国产精品无大码| 欧美成人a在线观看| 熟妇人妻久久中文字幕3abv| 99九九线精品视频在线观看视频| kizo精华| 欧美性猛交黑人性爽| 一级爰片在线观看| 国产精品,欧美在线| 欧美极品一区二区三区四区| 日韩中字成人| 亚洲av不卡在线观看| 国产三级在线视频| 日本黄色片子视频| 村上凉子中文字幕在线| 亚洲成人久久爱视频| 国产一级毛片七仙女欲春2| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品合色在线| 国产高清有码在线观看视频| 一个人看视频在线观看www免费| 亚洲精品国产成人久久av| 精品久久久噜噜| 99久久精品热视频| 在线播放国产精品三级| 免费av毛片视频| 中文字幕免费在线视频6| 欧美日本亚洲视频在线播放| 麻豆精品久久久久久蜜桃| 日韩人妻高清精品专区| 久久精品夜色国产| 波多野结衣高清无吗| 亚洲欧美中文字幕日韩二区| 国产免费一级a男人的天堂| 性色avwww在线观看| 久久久久久久久久久免费av| av卡一久久| 黄色一级大片看看| 一个人观看的视频www高清免费观看| 国产精品99久久久久久久久| 天美传媒精品一区二区| 国产免费男女视频| 国产黄片视频在线免费观看| 国产精品乱码一区二三区的特点| av视频在线观看入口| 搞女人的毛片| 日本爱情动作片www.在线观看| 91狼人影院| 我要看日韩黄色一级片| 91精品一卡2卡3卡4卡| 麻豆国产97在线/欧美| 国产免费又黄又爽又色| 中文精品一卡2卡3卡4更新| 国产色爽女视频免费观看| 国产成人91sexporn| 免费黄网站久久成人精品| 黄片wwwwww| 欧美成人精品欧美一级黄| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的 | 寂寞人妻少妇视频99o| 久久久精品94久久精品| 午夜激情欧美在线| 色5月婷婷丁香| 最新中文字幕久久久久| 好男人视频免费观看在线| 国产大屁股一区二区在线视频| 国产av码专区亚洲av| av又黄又爽大尺度在线免费看 | 亚洲欧美成人精品一区二区| 久久这里只有精品中国| 人妻少妇偷人精品九色| 亚洲最大成人中文| 国产久久久一区二区三区| 一区二区三区高清视频在线| 两个人视频免费观看高清| 欧美97在线视频| 99热这里只有是精品在线观看| 99热这里只有是精品50| 国内精品美女久久久久久| 国产成人精品婷婷| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| 精品久久久久久电影网 | 少妇被粗大猛烈的视频| 性插视频无遮挡在线免费观看| 精品人妻一区二区三区麻豆| 午夜福利在线观看免费完整高清在| 亚洲欧洲日产国产| 国产精品伦人一区二区| 精品熟女少妇av免费看| 中文字幕熟女人妻在线| 极品教师在线视频| 日本黄大片高清| 综合色丁香网| 啦啦啦观看免费观看视频高清| 少妇被粗大猛烈的视频| 99久久人妻综合| 免费av毛片视频| 观看免费一级毛片| 国产成人福利小说| 欧美日本视频| 校园人妻丝袜中文字幕| 免费看光身美女| 两个人视频免费观看高清| 99久久九九国产精品国产免费| 欧美成人a在线观看| 亚洲成人av在线免费| 亚洲精品国产av成人精品| 国产极品精品免费视频能看的| 中文乱码字字幕精品一区二区三区 | 麻豆久久精品国产亚洲av| 成人av在线播放网站| 91午夜精品亚洲一区二区三区| 99热全是精品| eeuss影院久久| 亚洲图色成人| 日韩欧美三级三区| 国产激情偷乱视频一区二区| 听说在线观看完整版免费高清| 亚洲无线观看免费| av在线播放精品| 色综合色国产| 美女cb高潮喷水在线观看| 成人毛片60女人毛片免费| 最近2019中文字幕mv第一页| 欧美日韩在线观看h| 午夜福利成人在线免费观看| av免费在线看不卡| 日韩一区二区视频免费看| 噜噜噜噜噜久久久久久91| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| 国产精品久久久久久久电影| 国产黄色视频一区二区在线观看 | 欧美高清性xxxxhd video| 美女国产视频在线观看| 国产精品电影一区二区三区| 国产精品女同一区二区软件| 网址你懂的国产日韩在线| 一区二区三区乱码不卡18| 精品熟女少妇av免费看| 热99re8久久精品国产| av播播在线观看一区| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| av在线老鸭窝| 综合色av麻豆| 黄色配什么色好看| 男人的好看免费观看在线视频| 国产女主播在线喷水免费视频网站 | 成年免费大片在线观看| 精品国产三级普通话版| 国产91av在线免费观看| 国产亚洲5aaaaa淫片| 亚洲欧美日韩高清专用| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放| 不卡视频在线观看欧美| 亚洲av成人av| 晚上一个人看的免费电影| 人人妻人人澡欧美一区二区| 国产精品一区www在线观看| h日本视频在线播放| 国国产精品蜜臀av免费| 国产久久久一区二区三区| 美女大奶头视频| 在线免费观看不下载黄p国产| 1000部很黄的大片| 亚洲综合精品二区| 亚洲经典国产精华液单| 国产黄片视频在线免费观看| 看免费成人av毛片| 国产亚洲午夜精品一区二区久久 | 国产高清三级在线| 噜噜噜噜噜久久久久久91| 久久久久国产网址| 国产亚洲av嫩草精品影院| 又粗又爽又猛毛片免费看| 久久草成人影院| 亚洲av.av天堂| 国产伦精品一区二区三区四那| 在线观看美女被高潮喷水网站| 国产白丝娇喘喷水9色精品| 2021少妇久久久久久久久久久| 亚洲精品亚洲一区二区| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出| 高清视频免费观看一区二区 | 亚洲天堂国产精品一区在线| a级毛片免费高清观看在线播放| 国产精品乱码一区二三区的特点| 国产精华一区二区三区| 久久99精品国语久久久| 亚洲国产精品国产精品| 欧美一区二区国产精品久久精品| 日本熟妇午夜| 午夜精品在线福利| 91aial.com中文字幕在线观看| 亚洲国产精品久久男人天堂| 国产一区二区在线av高清观看| 内射极品少妇av片p| 性色avwww在线观看| 亚洲精品自拍成人| 亚洲不卡免费看| 99在线人妻在线中文字幕| 91狼人影院| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 国产精品久久电影中文字幕| 国产成人福利小说| 黄色一级大片看看| av女优亚洲男人天堂| 亚洲乱码一区二区免费版| 国产乱人偷精品视频| 精品久久国产蜜桃| 国产精品一区二区三区四区免费观看| 国内精品宾馆在线| 深夜a级毛片| 真实男女啪啪啪动态图| 美女xxoo啪啪120秒动态图| 亚洲成人av在线免费| 少妇人妻精品综合一区二区| 亚洲怡红院男人天堂| 久久精品熟女亚洲av麻豆精品 | 久久久久九九精品影院| 日本五十路高清| 亚洲国产精品国产精品| 综合色av麻豆| 精品不卡国产一区二区三区| 人妻系列 视频| 美女国产视频在线观看| 在线播放无遮挡| 97热精品久久久久久| 午夜日本视频在线| 高清在线视频一区二区三区 | 国产免费又黄又爽又色| 日本爱情动作片www.在线观看| 午夜福利在线观看免费完整高清在| 亚洲一区高清亚洲精品| 亚洲av男天堂| 国产午夜精品论理片| 一区二区三区乱码不卡18| 国产亚洲5aaaaa淫片| 精品人妻一区二区三区麻豆| 听说在线观看完整版免费高清| 精品国产一区二区三区久久久樱花 | 亚洲丝袜综合中文字幕| 国产免费视频播放在线视频 | 久久精品综合一区二区三区| 亚洲经典国产精华液单| 非洲黑人性xxxx精品又粗又长| 精品久久久久久电影网 | 天天躁日日操中文字幕| 少妇熟女欧美另类| 国产精品久久久久久精品电影| 免费av观看视频| 日韩高清综合在线| 亚洲成人精品中文字幕电影| 天堂网av新在线| 两个人视频免费观看高清| 欧美精品一区二区大全| 亚洲不卡免费看| 久久久精品欧美日韩精品| 在线天堂最新版资源| 成人鲁丝片一二三区免费| 亚洲18禁久久av| 国内精品一区二区在线观看| 亚洲成人久久爱视频| 国产伦精品一区二区三区视频9| 看黄色毛片网站| 日韩欧美精品免费久久| 亚洲自拍偷在线| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩卡通动漫| 久久热精品热| 精品人妻一区二区三区麻豆| 免费观看人在逋| 国产女主播在线喷水免费视频网站 | 丝袜喷水一区| 69人妻影院| 精品人妻偷拍中文字幕| 精品久久久久久久久久久久久| 欧美性感艳星| 国产免费一级a男人的天堂| 亚洲久久久久久中文字幕| 欧美性猛交黑人性爽| 中文亚洲av片在线观看爽| 干丝袜人妻中文字幕| 国产av不卡久久| 亚洲成色77777| 韩国高清视频一区二区三区| 亚洲熟妇中文字幕五十中出| 看十八女毛片水多多多| 三级国产精品片| 亚洲av福利一区| 网址你懂的国产日韩在线| 人妻系列 视频| 亚洲国产欧美在线一区| 看十八女毛片水多多多| 青春草国产在线视频| 亚洲av成人av| av福利片在线观看| 女人被狂操c到高潮| 国产一区二区在线av高清观看| 三级男女做爰猛烈吃奶摸视频| 一级黄片播放器| 亚洲av成人精品一二三区| 欧美97在线视频| 久久99蜜桃精品久久| 一区二区三区四区激情视频| 成年版毛片免费区| 免费观看的影片在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产人妻一区二区三区在| 国产精品综合久久久久久久免费| 麻豆乱淫一区二区| 在线观看一区二区三区| av黄色大香蕉| 成年女人永久免费观看视频| 国产午夜精品论理片| 亚洲高清免费不卡视频| 亚洲av日韩在线播放| 麻豆乱淫一区二区|