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

    遼寧長??h附近海域顆粒有機碳和溶解有機碳交換能力的數(shù)值模擬*I模式的建立和驗證

    2015-03-20 02:37:25郭俊如夏穎穎鮑獻文劉玉龍陳新平姚志剛袁澤軼
    關(guān)鍵詞:長??h黃海潮汐

    郭俊如, 宋 軍, 夏穎穎, 鮑獻文,4, 劉玉龍, 陳新平, 姚志剛,4, 袁澤軼

    (1.中國海洋大學海洋環(huán)境學院,山東 青島 266100;2.國家海洋減災中心,北京 100194;3. 國家海洋信息中心,天津 300171;4.中國海洋大學物理海洋教育部重點實驗室,山東 青島 266100)

    ?

    遼寧長海縣附近海域顆粒有機碳和溶解有機碳交換能力的數(shù)值模擬*I模式的建立和驗證

    郭俊如1,2,4, 宋 軍1,3,4, 夏穎穎1, 鮑獻文1,4, 劉玉龍3, 陳新平2, 姚志剛1,4, 袁澤軼3

    (1.中國海洋大學海洋環(huán)境學院,山東 青島 266100;2.國家海洋減災中心,北京 100194;3. 國家海洋信息中心,天津 300171;4.中國海洋大學物理海洋教育部重點實驗室,山東 青島 266100)

    基于無結(jié)構(gòu)網(wǎng)格和有限體積法的海洋模型FVCOM(Finite Volume Coastal Ocean Model),并結(jié)合2個連續(xù)測流站、2個水位站的觀測資料,選取對養(yǎng)殖潛力影響最大的顆粒有機碳(POC)和溶解有機碳(DOC)作為本研究海區(qū)營養(yǎng)物質(zhì)的表征變量,依據(jù)不同季節(jié)和影響因子共設計了18種計算方案,對長??h附近海域營養(yǎng)物質(zhì)的交換能力分別進行了診斷和計算,對各個影響因子在營養(yǎng)物質(zhì)交換過程中的重要性進行了初步評估。研究得出,北黃海環(huán)流對該海域物質(zhì)交換能力的促進作用強于風生流的作用。

    北黃海; 長??h; 水交換能力; POC; DOC; FVCOM; 數(shù)值模擬

    長??h位于遼東半島東側(cè)的黃海北部,東與朝鮮半島隔海相望,西南與山東省廟島全島相對,西部和北部海域與普蘭店等相毗鄰,共轄112個島嶼,5處群礁以及51個明礁,東西跨度達到百余公里。長??h的經(jīng)濟支柱主要以水產(chǎn)養(yǎng)殖為主,全縣共有適合浮筏養(yǎng)殖淺海水面1.6萬hm2,適合底播增殖貝、藻類的海底面積9.7萬hm2,適合魚類人工流放增殖的海域面積13萬hm2。對長??h附近海域海水營養(yǎng)物質(zhì)交換問題的研究具有非?,F(xiàn)實的經(jīng)濟和社會意義。

    長??h附近海域的海水交換及其營養(yǎng)物質(zhì)輸運受到多種因素的制約和影響[1-3]。首先研究海域內(nèi)島嶼星羅密布,岸線水深復雜多變,潮汐潮流特征顯著,從而導致該海區(qū)內(nèi)溫鹽和環(huán)流結(jié)構(gòu)異常復雜[4-8];其次,風場等大氣強迫對研究海域內(nèi)水動力環(huán)境的演變過程也有著十分顯著的影響[9-14],這一點可以從長??h所處的北黃海海域溫鹽及環(huán)流場的季節(jié)演變中得到證實[2,15-18],從而可以預見到研究海域內(nèi)營養(yǎng)物質(zhì)的交換也應當存在比較明顯的季節(jié)特征;第三,在近海海域營養(yǎng)物質(zhì)的分布和輸運的研究中,徑流輸入也是一個應當考慮的因素,一方面徑流所攜入海的各種營養(yǎng)鹽等提供了近海營養(yǎng)物質(zhì)變化的“源”,另一方面徑流自身也影響和改變了周邊的水動力環(huán)境場,進而影響了營養(yǎng)物質(zhì)的分布和輸運;此外,長??h附近海域所處的北黃海背景溫鹽場以及環(huán)流結(jié)構(gòu)的變化也對研究海域內(nèi)的營養(yǎng)物質(zhì)的輸運有著重要影響,因此,對長??h附近海域營養(yǎng)物質(zhì)交換問題的研究中應當綜合考慮上述所有影響因子。這其中一個比較重要的問題即是對交換過程中各影響因子作用及重要性的研究,這有助于我們加深對該海區(qū)營養(yǎng)物質(zhì)交換過程的理解和認識,亦即本文的主要研究目標。

    1 研究方法

    基于上述研究目標,本文采用了當前國際先進的非結(jié)構(gòu)有限體積模型FVCOM(Finite Volume Coastal Ocean Model)[3],并結(jié)合現(xiàn)場觀測,針對影響長??h營養(yǎng)物質(zhì)輸運的各因子設計了一系列的敏感性試驗,對研究海區(qū)內(nèi)營養(yǎng)物質(zhì)輸運過程進行了研究。特別值得一提的是,F(xiàn)VCOM模型基于非結(jié)構(gòu)的三角網(wǎng)格,從而在岸線地形復雜的島嶼附近具有較高的空間分辨率,可以較好的模擬復雜地理特征下的水動力學過程,因而特別適合本文的研究海域-位于北黃海的長海縣附近海域。

    1.1 模型配置

    研究海域為長海縣全海域(見圖1),模型計算區(qū)域的范圍為31.8°N~41.0°N,117.5°E~127.0°E(見圖2),包括整個渤黃海海域,最大網(wǎng)格分辨率3′(約5000m)。為盡可能的模擬出長??h海域的海流狀況,對長海縣海域進行了加密,最小水平網(wǎng)格分辨率達到20m,平均網(wǎng)格分辨率小于100m,垂向采用sigma坐標,共分6層。網(wǎng)格的分辨率幾乎達到了模擬真實地形的程度[3]。模型的底摩擦系數(shù)基于張繼才與呂咸青[19]的同化研究結(jié)果。外海開邊界的溫度、鹽度與流通量來自經(jīng)過驗證的大區(qū)ROMS模型的結(jié)果[3,24],而開邊界處用于預報諧振潮潮位與潮流強迫條件的調(diào)和常數(shù)則來自于NAO99 全球潮汐模型所提供的8 個分潮(M2、S2、N2、K2、K1、O1、P1、Q1)的數(shù)據(jù)(http://www.miz.nao.ac.jp/staffs/nao99/index_En.html)。海表熱通量條件來自于空間分辨率為2(°)×2(°)、時間分辨率為每日一次,來源于美國國家環(huán)境預報中心數(shù)據(jù)(http://nomads.ncep.noaa.gov/txt_descriptions/servers.shtml)。同樣為每日一次,空間分辨率為1/4(°)×1/4(°)的海表面風應力場數(shù)據(jù)則來源于 the Blended Sea Winds(http://www.ncdc.noaa.gov/oa/rsad/seawinds.html)。模型岸邊界條件由如下公式給出:

    圖1 長??h海域行政區(qū)域劃分圖(綠線)與驗潮站、測流站站位分布圖(a)以及水深等值線圖(b)Fig.1 The stations of water lever(W)、current(S) observation(a) and Topography(b)

    ((a)FVCOM模型的計算區(qū)域(31.8°N~41°N,117.5°E~127°E)與水深;(b)長??h海域網(wǎng)格配置。(a)Calculation area and depth of fvcom; (b)The gricl configure of Changhai country area.)

    圖2 模型區(qū)域及網(wǎng)格配置
    Fig.2 Model area and its mesh-system configure

    (1)

    (2)

    其中:C為物質(zhì)的濃度;D為水深(平均水深與水位的和);u、v和w分別為x、y和z方向的流速分量;W為源匯項;AH與KH為物質(zhì)的水平擴散系數(shù)與垂向擴散系數(shù)。需要說明的是,由于長??h所轄海域無陸地邊界,海域附近的鴨綠江等主要主要徑流攜帶大量淡水和有機物質(zhì)通量注意該海域的周邊海域(外海),所以在以長??h海域為研究目標時,徑流作為影響該海域外海的條件引入,而受水質(zhì)監(jiān)測資料的限制,在考慮該海域與外海物質(zhì)交換時,外海海域的有機物質(zhì)濃度只能給出平均值。所以,本文中徑流的影響主要是指徑流影響水動力環(huán)境而進一步對物質(zhì)交換造成的影響,而徑流的物質(zhì)通量的影響已經(jīng)包括在所有列及的模型實驗中。

    1.2 模型驗證

    為了獲得長海縣海域近期的水文動力特征以對模型進行驗證。本研究開展了該海域冬季(2007年1月) 2個站位1個月的連續(xù)水位觀測和2個站位海流大、中、小潮的25h連續(xù)觀測[25](站位分布如圖1所示,這里的中潮觀測選取的是大潮日和小潮日中間的一天,即2007年1月6日)。圖3~6為模型控制實驗的計算結(jié)果與實測水位、流速資料的對比圖,模式結(jié)果與實測資料的誤差已達到了此次研究的需要。其中水位的計算結(jié)果與實測值幾乎吻合,流速與流向結(jié)果的誤差小于10%。于是,基于這一相對可靠的數(shù)值模型,作者可以更進一步的加入各種因素并研究其對作者所關(guān)心問題的影響程度。

    圖3 W1驗潮站處的計算水位(藍線)與實測水位(紅線)對比圖Fig.3 Validation of the result of water level (blue line) from ocean model with observation (red line) in the station of W1

    圖4 W2驗潮站處的計算水位(藍線)與實測水位(紅線)對比圖

    2 敏感性數(shù)值試驗設計

    基于經(jīng)過觀測數(shù)據(jù)驗證的潮汐模型[3],主要針對冬、夏2個典型季節(jié)下潮汐潮流、環(huán)流、風和徑流等不同影響因子的作用,共設計了18種不同強迫條件下的數(shù)值模擬實驗,分別包括冬、夏季純潮驅(qū)動試驗(方案1,11)以及冬、夏季不同強迫條件組合(方案2~10,方案12~18)。夏季計算時間段為2006年6~8月,冬季計算時間段為2006年12月~2007年2月,具體設計方案見表1。

    通過對這些數(shù)值計算的結(jié)果進行診斷分析,對各影響因子在研究海區(qū)內(nèi)水交換及營養(yǎng)物質(zhì)交換過程中的作用將在本文和本文后續(xù)文章中進行較全面的評估。本文基于溶解有機物濃度與顆粒有機物濃度的調(diào)查觀測資料,并結(jié)合經(jīng)典的對流擴散模型[20-23], 計算了各時間段內(nèi)海區(qū)有機物質(zhì)總量在不考慮生物消耗與生產(chǎn)情況下的絕對變化。

    圖5 S1測流站大潮底層流速(左)流向(右)模擬結(jié)果(虛線)與實測(實線)對比

    圖6 S2測流站小潮表層流速(左)流向(右)模擬結(jié)果(虛線)與實測(實線)對比

    夏季Sumer冬季W(wǎng)inter夏季純潮正壓模型(方案11)冬季純潮正壓模型(方案1)潮+夏季溫鹽(方案12)潮+冬季溫鹽(方案2)潮+夏季溫鹽+夏季徑流(方案13)潮+冬季溫鹽+冬季徑流(方案3)潮+夏季溫鹽+夏季徑流+夏季風(方案14)潮+冬季溫鹽+冬季徑流+冬季風(方案4)潮+夏季溫鹽+夏季雙倍流量徑流+夏季風(方案15)潮+冬季溫鹽+冬季雙倍流量徑流+冬季風(方案5)潮+夏季溫鹽+夏季徑流+夏季雙倍風速風(方案16)潮+冬季溫鹽+冬季徑流+冬季雙倍風速風(方案6)潮+夏季溫鹽+夏季徑流+夏季氣候態(tài)風(方案17)潮+冬季溫鹽+冬季徑流+冬季氣候態(tài)風(方案7)潮+夏季風(方案18)潮+冬季風(方案8)潮+冬季氣候態(tài)風(方案9)潮+冬季溫鹽+冬季徑流+冬季1.5倍風速風(方案10)

    (3)

    其中:Mt表示時間段t內(nèi)所研究物質(zhì)的量的變化大??;Q為某時刻研究區(qū)域物質(zhì)的總量;He為三維計算單元e的垂向高度;Se為三維計算單元e的網(wǎng)格面積;Ce為三維計算單元e內(nèi)的物質(zhì)濃度;t1為時間段t的結(jié)束時刻;t0為時間段t的開始時刻;本文中Mt的測算單位為噸?;谑?的計算結(jié)果,作者將在下文中對不同時間段內(nèi)各種影響因素對研究海域內(nèi)營養(yǎng)物質(zhì)交換所造成的影響進行討論。

    3 潮汐潮流對長??h海域水交換與營養(yǎng)物質(zhì)交換的影響

    長??h海域潮流運動呈現(xiàn)明顯的旋轉(zhuǎn)流特性,潮汐類型表現(xiàn)為正規(guī)半日潮。大區(qū)的潮流規(guī)律為漲潮時海流自黃海東部流入,落潮時由渤海自渤海海峽流出。長??h及其各鄉(xiāng)海域在此潮流背景下,不同潮時(大中小潮,其中本文中潮定義為大潮日與小潮日中間一天的潮汐)營養(yǎng)物質(zhì)(顆粒有機碳與溶解有機碳)的交換情況對比如圖7~9所示。

    圖7 冬季長海縣海域各站在不同潮時下顆粒有機碳(POC)交換量Fig.7 Different exchange strength of POC in different kinds of winter tide days in country Changhai

    3.1 冬季交換情況

    (1)代表顆粒有機物的顆粒有機碳(POC)的交換情況(見圖7)。

    由圖7可以看出,長??h大部分海域在中潮時POC交換能力最強,其次是大潮時,小潮時POC交換能力最小。但在獐子島與海洋鄉(xiāng),POC交換能力在大潮時最大,次之為中潮與小潮。從POC交換能力的垂直結(jié)構(gòu)上來看,由于積分計算是根據(jù)觀測時所采用的六點法分層(即表面、0.2H、0.4H、0.6H、0.8H、底層)所分出的5個垂向平分的水層進行的,所以各層的交換量具有直觀的對比性。僅就潮汐作用下的水交換量來看,除海洋鄉(xiāng)外,其他區(qū)域皆為表層的POC交換能力最強,海洋鄉(xiāng)POC交換能力表現(xiàn)為底層較強。各海域具體的交換能力可以參考圖7中的數(shù)值。

    (2) 代表溶解有機物的溶解有機碳(DOC)的交換情況(見圖8)。

    圖8 冬季長海縣海域各站在不同潮時下溶解有機碳(DOC)交換量Fig.8 Different exchange strength of DOC in different kinds of winter tide days in country Changhai

    由圖8可以看出,與POC的情況不同,DOC的交換能力在獐子島與海洋鄉(xiāng)表現(xiàn)為中潮時最強,次之為大潮時,小潮時最小,而在大小長山則表現(xiàn)為大潮時最大,中潮與小潮次之。廣鹿鄉(xiāng)海域則表現(xiàn)為小潮時交換能力最強。各海域具體的交換潛力可以參考圖8中的數(shù)值。值得注意的是,廣鹿鄉(xiāng)在大潮的時候DOC的交換潛力是正的,即增加的,而在中潮和小潮時卻是減小的。

    3.2 夏季交換情況

    圖9為基于方案14的長??h海域數(shù)值模擬計算結(jié)果。從圖中可以看出,與冬季的交換情況不同,夏季長??h海域POC與DOC不同潮時下交換能力差別顯著。大潮時交換量明顯較大,其次為中潮,小潮時交換量最小。從垂向分布結(jié)構(gòu)來看,可以看出無論是POC還是DOC,均表現(xiàn)為中間3層的交換量較表層、底層交換量高。

    圖9 夏季不同潮時下長??h海域POC(左圖)和DOC(右圖)交換量對比Fig.9 Different exchange strength of POC and DOC in different kinds of summer tide days in country Changhai

    4 環(huán)流對長海縣海域水交換與營養(yǎng)物質(zhì)交換的影響

    環(huán)流一般包括密度流、局地風生流和外海傳入流,本研究所涉及風和風應力都是指局地風。為了方便細化對比,在此局地風生流對比時提到的環(huán)流指的是無局地風影響的環(huán)流。北黃海的環(huán)流主要包括黃海暖流及黃海沿岸流,其中黃海暖流是外海水輸入的主要來源。雖然北黃海主要環(huán)流流速僅有最大潮流的十分之一,但考慮到其流向終年比較穩(wěn)定,與周期往復的潮流相比,北黃海的環(huán)流對研究海區(qū)水體以及物質(zhì)輸運有著更加顯著的影響。此外,風應力對研究海域的環(huán)流,尤其是表層環(huán)流也有著重要影響,具體表現(xiàn)為冬季表層流多南向流,而夏季則多偏北流,表層海流具有明顯的風海流性質(zhì)。由于黃海環(huán)流的流速有著冬季強、夏季弱的特點,因此本文中重點對冬季大潮期間黃海北部環(huán)流背景場與風應力作用對研究海區(qū)營養(yǎng)物質(zhì)交換情況的影響進行分析。

    圖10描述的是長??h海域各層在僅有潮汐作用的正壓潮汐模型(方案1)、不加風與徑流的斜壓潮汐模型(方案2)、不加風加徑流的斜壓潮汐模型(方案3)以及加風加徑流的斜壓潮汐模型(方案4)3種情況作用下,營養(yǎng)物質(zhì)的交換情況的對比。

    圖10 長??h海域冬季潮汐、黃海北部環(huán)流(無局地風影響)和局地風生流對POC(左圖)和DOC(右圖)交換量的影響對比Fig.10 Different exchange strength of POC and DOC are influenced by tide, circulation and wind in country Changhai

    從圖中可以看出,風對表層的營養(yǎng)物質(zhì)輸運起著決定性的作用,并且可能會使表層輸運計算結(jié)果的數(shù)值反號,即表層原本在潮汐作用下會增加的量由于風應力作用反而會減小。但風應力對下層的水層,特別是水深較淺的情況下,由于質(zhì)量守恒會產(chǎn)生反向的補償流,一般會加強僅有潮汐作用情況下的營養(yǎng)物質(zhì)輸運程度。

    此外可以看出,北黃海環(huán)流對研究海區(qū)內(nèi)營養(yǎng)物質(zhì)的輸運一般呈加強的趨勢,其增加的量甚至可以達到一倍以上。這也與上面的理論分析結(jié)果一致。

    值得說明的是,雖然在風應力的持續(xù)作用下會對流場特別是表層流場有著顯著的影響,但其對整個水層的影響存在明顯的垂直變化,因此整體而言,風應力對研究海區(qū)內(nèi)營養(yǎng)物質(zhì)輸運的影響較北黃海背景場環(huán)流偏弱,尤其是考慮到二者間的非線性作用時,北黃海環(huán)流背景場對營養(yǎng)物質(zhì)輸運的影響會更加顯著。

    5 徑流輸入對長??h海域水交換與營養(yǎng)物質(zhì)交換的影響

    一般認為海區(qū)周邊的徑流輸入會對海域內(nèi)余流的流速流向產(chǎn)生一定影響,此外,徑流入海所輸入的營養(yǎng)物質(zhì)對研究海域內(nèi)營養(yǎng)物質(zhì)的改變也有顯著的影響。因此本文考慮了渤海和北黃海海域周邊的黃河、海河、灤河、遼河、愛河、大洋河及鴨綠江等主要徑流輸入對長??h海域營養(yǎng)物質(zhì)輸運情況的影響,并與無徑流輸入的控制試驗進行了對比。

    圖11分別給出了長??h海域各層在僅有潮汐作用的正壓潮汐模型(方案1),不加風與徑流的斜壓潮汐模型(方案2),以及不加風但加徑流的斜壓潮汐模型(方案3)3種情況作用下營養(yǎng)物質(zhì)的交換情況的對比。由上節(jié)的分析與圖11可以看出,徑流(主要是鴨綠江)對長海縣海域營養(yǎng)物質(zhì)的輸運有一定的影響,但影響不占主導地位,不超過5%。將徑流的流量在平均水平上人為加大一倍之后,其影響也不超過10%。

    圖11 長??h海域夏季徑流對POC(左圖)和DOC(右圖)交換量的影響Fig.11 Different exchange strength of POC(Left) and DOC(Right) are influenced by river in country Changhai

    6 結(jié)論

    本文應用無結(jié)構(gòu)網(wǎng)格和有限體積法的海洋模型FVCOM(Finite Volume Coastal Ocean Model),結(jié)合4個連續(xù)測流站、3個水位站以及1個坐底ADCP的觀測資料,選取顆粒有機物(POM)和溶解有機物(DOM)作為本研究海區(qū)營養(yǎng)物質(zhì)的表征變量,依據(jù)不同季節(jié)和影響因子共設計了18種計算方案,對潮汐潮流、環(huán)流、風和徑流在營養(yǎng)物質(zhì)交換過程中的重要性進行了評估。研究結(jié)論如下:

    (1)長??h海域冬、夏季不同潮時下顆粒有機物、溶解有機物的交換能力不同。夏季長海縣海域,POC與DOC在不同潮時下交換能力差別顯著,大潮時交換量明顯較大,其次為中潮,小潮時交換能力最小。從垂向分布結(jié)構(gòu)來看,無論是POC還是DOC,均表現(xiàn)出中層水的交換量最大,表層水、底層水的交換量次之的特點;

    (2)北黃海環(huán)流對研究海區(qū)內(nèi)營養(yǎng)物質(zhì)的輸運一般呈加強的趨勢,其增加的量甚至可以達到一倍以上;

    (3)整體而言,風應力對該海區(qū)內(nèi)營養(yǎng)物質(zhì)輸運的影響小于北黃海環(huán)流的影響;

    (4)徑流(主要是鴨綠江)對長??h海域營養(yǎng)物質(zhì)的輸運有一定的影響,但影響在各因子總的影響中所占的比重小于5%。

    [1] 鮑獻文, 宋軍, 姚志剛, 等. 北黃海潮流、余流垂直結(jié)構(gòu)及其季節(jié)變化 [J]. 中國海洋大學學報:自然科學版, 2010, 40(11): 11-18.

    [2] 趙俊生, 耿世江. 北黃海內(nèi)潮對潮流垂直結(jié)構(gòu)的影響 [J]. 海洋學報, 1990, 12(6): 677-692.

    [3] 宋軍, 水交換模型的理論方法及應用研究 [D].青島: 中國海洋大學海洋環(huán)境學院, 2010.

    [4] 宋新, 林霄沛, 王悅. 夏季黃海冷水團的多年際變化及原因淺析 [J]. 廣東海洋大學學報, 2009(3): 59-63.

    [5] 江蓓潔, 鮑獻文, 吳德星, 等. 北黃海冷水團溫、鹽多年變化特征及影響因素 [J]. 海洋學報: 中文版, 2007(4): 1-10.

    [6] 張以懇, 楊玉玲. 夏季北黃海冷水團多年變化特征分析 [J]. 海洋預報, 1996(4): 15-21.

    [7] 蘇紀蘭, 黃大吉. 黃海冷水團的環(huán)流結(jié)構(gòu) [J]. 海洋與湖沼, 1995(S1): 1-7.

    [8] 赫崇本, 汪圓祥, 雷宗友, 等. 黃海冷水團的形成及其性質(zhì)的初步探討 [J]. 海洋與湖沼, 1959(1): 1-14.

    [9] 王輝武, 于非, 呂連港, 等. 冬季黃海暖流區(qū)的空間變化和年際變化特征 [J]. 海洋科學進展, 2009(2): 140-148.

    [10] 許一, 于非, 張志欣, 等. 冬季黃海暖流的診斷計算 [J]. 海洋科學進展, 2005(4): 398-407.

    [11] 湯毓祥, 鄒娥梅, Lie Heung-Jae. 冬至初春黃海暖流的路徑和起源 [J]. 海洋學報: 中文版, 2001(1): 1-12.

    [12] 藍淑芳. 黃海暖流水的調(diào)查研究 [J]. 海洋科學, 1993(1): 38-40.

    [13] 樂肯堂. 冬季黃海暖流水的起源 [J]. 海洋學報: 中文版, 1992(2): 9-19.

    [14] 李繁華, 孔祥德. 關(guān)于黃海暖流某些特征的分析 [J]. 海洋湖沼通報, 1987(4): 1-6.

    [15] 修日晨, 李繁華. 渤海及北黃海潮流場的基本特征 [J]. 海洋科學, 1989(5): 1-7.

    [16] 戚建華, 蘇育嵩. 黃海潮生陸架鋒的數(shù)值模擬研究 [J]. 海洋與湖沼, 1998, 29(3): 247-254.

    [17] 周玲, 高煥臣. 渤, 黃海沿岸潮災性質(zhì)分析 [J]. 海洋預報, 1993, 10(1): 37-39.

    [18] 王以嬌, 張延廷. 北黃海天文潮與風暴潮耦合水位的數(shù)值預報方法 [J]. 黃渤海海洋, 1993, 11(1): 1-8.

    [19] 張繼才, 呂咸青.渤、黃、東海二維潮汐模式底摩擦系數(shù)的反演研究 [J].計算力學學報, 2007, 24(4): 430-435.

    [20] 董禮先, 蘇紀蘭.象山港水交換數(shù)值研究 I. 對流擴散型的水交換模式 [J]. 海洋與湖沼, 1999, 30(4): 410-415.

    [21] Sandery, Kampf. Winter-spring flushing of bass strait, South-Eastern Australia: a numerical modelling study [J]. Estuarine, Coastal and Shelf Science, 2005, 63(1-2): 23-31.

    [22] Otsmann M, Astok V, Suursaar ü. A model for water exchange between the Baltic Sea and the Gulf of Riga [J]. Nordic Hydrology, 1997, 28(4-5): 351-364.

    [23] Krauss W, Brügge B. Wind-produced water exchange between the deep basins of the Baltic Sea [J]. Journal of Physical Oceanography, 1991, 21(3): 373-384.

    [24] Song J, Xue H J, Bao X W, et al. A spectral mixture model analysis of the Kuroshio variability and the water exchange between the Kuroshio and the East China Sea [J]. Chinese Journal of Oceanology and Limnology, 2011, 29(2): 446-459

    [25] Song Jun, Guo Junru, Qiao Lulu, et al. Vertical structure and seasonal variation of the tide currents and residual currents in the North Yellow Sea [J]. Acta Oceanologica Sinica, 2011, 33(5): 14-22.

    責任編輯 龐 旻

    Numerical Simulation of Exchange of POC and DOC in the Sea Waters Adjacent to Changhai County, Liaoning Province, China Part I: Model Development and Verification

    GUO Jun-Ru1,2,4, SONG Jun1, 3,4, XIA Ying-Ying1, BAO Xian-Wen1,4, LIU Yu-Long3, CHEN Xin-Ping2, YAO Zhi-Gang1,4, YUAN Ze-Yi3

    (1 College of Physical and Environmental Oceanography, Ocean University of China, Qingdao 266100, China;2 National Marine Hazard Mitigation Service, SOA, Beijing 100194, China;3 National Marine Data and Information Service, SOA, Tianjin 300171, China; 4 The Key Caboratory of Physical Oceanography, Ministry of Education, Ocean University of China, Qingdao 266100, China.)

    By making use of FVCOM(An Unstructured Grid, Finite-Volume Coastal Ocean Model) and being based on the consecutively observed data from four hydrometric stations, three waterlevel instruments and one ADCP in this sea area, choosing POC and DOC which have the most significant impact on the potential of mariculture as study variables to represent the seeds of nutrients, we calculated and diagnostically analyzed the exchange capacity of nutrients County Chinghai in eighteen conditions based on the proportions of different factors and seasons. The importance assessment of each factor during the nutrients exchange is provided. This study indicated that the material exchange in North Yellow Sea was strongly influenced by baroclinic flow and wind-driven current, and the former had more effect.

    Yellow Sea; County Chinghai; water exchange ability; POC; DOC; FVCOM; numorical simulation

    國家自然科學基金項目(41206013;41106004;41376014);國家海洋局青年科學基金重點項目(2012202;2013203; 2012223);教育部物理海洋重點實驗室開放基金項目;海洋公益性行業(yè)科研專項(201205018;201005019);國家建設高水平大學公派研究生項目(留金出[2008]3019,[2012]3013)資助

    2013-11-05;

    2013-12-07

    郭俊如(1986-),女,博士生。

    ** 通訊作者: E-mail:thunder098@hotmail.com

    P731.2

    A

    1672-5174(2015)04-018-09

    10.16441/j.cnki.hdxb.20130410

    猜你喜歡
    長??h黃海潮汐
    人民日報:電纜通海島 供電更可靠
    潮汐與戰(zhàn)爭(上)
    東方濕地 黃海明珠
    華人時刊(2020年21期)2021-01-14 01:32:28
    大愛小和愛
    黃海簡介
    絕美海灘
    潮汐式灌溉控制系統(tǒng)的設計及應用
    電子制作(2017年9期)2017-04-17 03:00:56
    三角恒等變換去哪兒了
    干法紙的潮汐
    生活用紙(2016年6期)2017-01-19 07:36:25
    長??h工會服務職工辦實事系列活動掠影
    国产人妻一区二区三区在| 欧美成人精品欧美一级黄| 欧美性猛交╳xxx乱大交人| 亚洲一级一片aⅴ在线观看| 国产精品三级大全| 日韩大片免费观看网站| 夫妻午夜视频| 国产av不卡久久| 街头女战士在线观看网站| 全区人妻精品视频| 亚洲色图av天堂| 亚洲成色77777| 国产白丝娇喘喷水9色精品| 欧美激情久久久久久爽电影| 欧美不卡视频在线免费观看| 国产精品一区二区在线观看99 | 男人舔奶头视频| 亚洲乱码一区二区免费版| 男女下面进入的视频免费午夜| 久久久久久久午夜电影| 一级毛片我不卡| 中文字幕亚洲精品专区| 九九爱精品视频在线观看| 哪个播放器可以免费观看大片| 欧美不卡视频在线免费观看| 国产伦理片在线播放av一区| 哪个播放器可以免费观看大片| 精品欧美国产一区二区三| 久99久视频精品免费| 97精品久久久久久久久久精品| 午夜激情福利司机影院| 国产白丝娇喘喷水9色精品| 大陆偷拍与自拍| 亚洲精品aⅴ在线观看| 亚洲综合色惰| 一区二区三区乱码不卡18| 午夜福利成人在线免费观看| 成年版毛片免费区| 国产一区二区在线观看日韩| 精品久久久久久久末码| 欧美极品一区二区三区四区| 波野结衣二区三区在线| 日韩伦理黄色片| 乱人视频在线观看| 欧美xxxx性猛交bbbb| 国产一区有黄有色的免费视频 | 精品一区二区免费观看| 免费观看无遮挡的男女| 啦啦啦啦在线视频资源| 午夜免费男女啪啪视频观看| 免费看a级黄色片| 成人特级av手机在线观看| 好男人在线观看高清免费视频| 极品少妇高潮喷水抽搐| 国产精品1区2区在线观看.| 五月伊人婷婷丁香| 国产综合精华液| 99热6这里只有精品| 啦啦啦韩国在线观看视频| 日本-黄色视频高清免费观看| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 国产69精品久久久久777片| 亚洲av中文字字幕乱码综合| 婷婷色麻豆天堂久久| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 最近2019中文字幕mv第一页| 一级爰片在线观看| 欧美三级亚洲精品| 国产不卡一卡二| 久久久久精品久久久久真实原创| 一级毛片黄色毛片免费观看视频| 久久亚洲国产成人精品v| 亚洲精品一二三| 国产免费又黄又爽又色| 晚上一个人看的免费电影| 久99久视频精品免费| 久热久热在线精品观看| 天堂影院成人在线观看| 午夜免费男女啪啪视频观看| 白带黄色成豆腐渣| 精品一区二区三区人妻视频| 自拍偷自拍亚洲精品老妇| 亚洲欧美成人综合另类久久久| 美女黄网站色视频| 亚洲精品国产成人久久av| 久久精品久久精品一区二区三区| 国产亚洲av嫩草精品影院| 亚洲久久久久久中文字幕| 黄色配什么色好看| 国产av不卡久久| 亚洲精品国产av成人精品| 熟妇人妻不卡中文字幕| 成人亚洲精品一区在线观看 | 亚洲国产高清在线一区二区三| 精品一区二区免费观看| 午夜精品在线福利| 亚洲人成网站高清观看| 纵有疾风起免费观看全集完整版 | 十八禁国产超污无遮挡网站| 日韩av不卡免费在线播放| 国产精品三级大全| 国产亚洲一区二区精品| 国产一区二区三区综合在线观看 | 哪个播放器可以免费观看大片| 69av精品久久久久久| 久久精品国产鲁丝片午夜精品| 欧美日韩综合久久久久久| 色吧在线观看| 国产老妇伦熟女老妇高清| 秋霞伦理黄片| 深夜a级毛片| 日日摸夜夜添夜夜爱| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 夫妻午夜视频| 久久草成人影院| 精品酒店卫生间| 久久99热这里只有精品18| 日韩欧美精品免费久久| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人 | 免费少妇av软件| 丝瓜视频免费看黄片| 久久久a久久爽久久v久久| 午夜久久久久精精品| 国产亚洲精品av在线| av在线蜜桃| 最新中文字幕久久久久| 亚洲国产欧美人成| 久久午夜福利片| 久久精品久久久久久久性| 日韩视频在线欧美| 白带黄色成豆腐渣| 久久99蜜桃精品久久| 国产av不卡久久| 大香蕉97超碰在线| 国产一区有黄有色的免费视频 | 亚洲欧美清纯卡通| 亚洲丝袜综合中文字幕| 国产精品福利在线免费观看| 好男人视频免费观看在线| 国产精品麻豆人妻色哟哟久久 | 伦理电影大哥的女人| 亚洲熟女精品中文字幕| 国产精品久久视频播放| 夜夜爽夜夜爽视频| 日韩人妻高清精品专区| 舔av片在线| 熟妇人妻不卡中文字幕| 在线天堂最新版资源| 美女国产视频在线观看| 最后的刺客免费高清国语| 成人av在线播放网站| 欧美高清成人免费视频www| 亚洲欧美精品专区久久| 亚洲婷婷狠狠爱综合网| 中文字幕人妻熟人妻熟丝袜美| 午夜福利高清视频| 久久99热这里只频精品6学生| 欧美xxxx性猛交bbbb| 91狼人影院| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区三区| 女人被狂操c到高潮| 白带黄色成豆腐渣| 久久久久久伊人网av| 亚洲国产色片| 精品国内亚洲2022精品成人| 永久免费av网站大全| 听说在线观看完整版免费高清| 久久国内精品自在自线图片| 国产在视频线在精品| videos熟女内射| 青青草视频在线视频观看| 成人性生交大片免费视频hd| 国产成人精品一,二区| 久久久久久久久久成人| 国产精品一区二区三区四区免费观看| 午夜福利在线在线| 99久久精品热视频| 国产精品av视频在线免费观看| 欧美激情在线99| 国产一区有黄有色的免费视频 | 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产| 亚洲不卡免费看| av天堂中文字幕网| 久久久久久久久大av| 97热精品久久久久久| 国产一区有黄有色的免费视频 | 久久久精品免费免费高清| 亚洲国产高清在线一区二区三| 26uuu在线亚洲综合色| 18禁在线播放成人免费| 狠狠精品人妻久久久久久综合| 夜夜看夜夜爽夜夜摸| 亚洲熟妇中文字幕五十中出| 久热久热在线精品观看| 亚洲av二区三区四区| 色哟哟·www| 人妻系列 视频| 黄色日韩在线| 久久精品综合一区二区三区| av线在线观看网站| 色综合站精品国产| 国产精品久久久久久久久免| 中文天堂在线官网| 成人漫画全彩无遮挡| videos熟女内射| 1000部很黄的大片| 天天躁日日操中文字幕| 蜜桃亚洲精品一区二区三区| 国产高清不卡午夜福利| 偷拍熟女少妇极品色| 高清毛片免费看| ponron亚洲| 亚洲aⅴ乱码一区二区在线播放| 欧美潮喷喷水| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 久久久久久久久久人人人人人人| 舔av片在线| 丝袜喷水一区| 亚洲av国产av综合av卡| 一本久久精品| 亚洲国产日韩欧美精品在线观看| 久久99热这里只有精品18| 午夜福利视频1000在线观看| 精品酒店卫生间| 美女黄网站色视频| 特大巨黑吊av在线直播| av.在线天堂| 国产精品日韩av在线免费观看| 在现免费观看毛片| 亚洲人成网站高清观看| 伦理电影大哥的女人| 五月伊人婷婷丁香| 色综合站精品国产| 好男人在线观看高清免费视频| 联通29元200g的流量卡| 亚洲欧美日韩无卡精品| 中文乱码字字幕精品一区二区三区 | 国产精品麻豆人妻色哟哟久久 | 久久草成人影院| 黄色日韩在线| av在线天堂中文字幕| 精品99又大又爽又粗少妇毛片| 国产爱豆传媒在线观看| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 秋霞伦理黄片| 全区人妻精品视频| 人人妻人人澡欧美一区二区| 午夜免费男女啪啪视频观看| 亚洲成色77777| 只有这里有精品99| 丰满乱子伦码专区| 色哟哟·www| 水蜜桃什么品种好| 欧美高清性xxxxhd video| 国产免费又黄又爽又色| 欧美日韩一区二区视频在线观看视频在线 | www.色视频.com| 亚洲欧美日韩卡通动漫| 99久国产av精品| 中文字幕av在线有码专区| 国产av码专区亚洲av| 国产不卡一卡二| 日本免费a在线| 男人舔奶头视频| 少妇被粗大猛烈的视频| 久久久亚洲精品成人影院| 国产欧美另类精品又又久久亚洲欧美| 精品欧美国产一区二区三| 欧美日韩亚洲高清精品| 久久久精品免费免费高清| 欧美日韩综合久久久久久| 亚洲av中文av极速乱| 亚洲国产欧美人成| 搡老妇女老女人老熟妇| 日日干狠狠操夜夜爽| 久久精品夜色国产| 成人亚洲精品av一区二区| 亚洲三级黄色毛片| 欧美最新免费一区二区三区| 超碰av人人做人人爽久久| 日韩欧美国产在线观看| 欧美bdsm另类| 人妻系列 视频| 永久免费av网站大全| 最近中文字幕高清免费大全6| 看黄色毛片网站| 久久久久精品久久久久真实原创| 亚洲一区高清亚洲精品| 国产一级毛片七仙女欲春2| 午夜福利视频精品| 真实男女啪啪啪动态图| freevideosex欧美| 婷婷色综合www| 人人妻人人澡人人爽人人夜夜 | 日韩欧美精品免费久久| 国产精品一区二区三区四区免费观看| 亚洲18禁久久av| 欧美丝袜亚洲另类| 人妻少妇偷人精品九色| 97超碰精品成人国产| 久久热精品热| 在线 av 中文字幕| 91狼人影院| av免费在线看不卡| 亚洲av免费高清在线观看| 成人无遮挡网站| 噜噜噜噜噜久久久久久91| 禁无遮挡网站| 亚洲天堂国产精品一区在线| 不卡视频在线观看欧美| av天堂中文字幕网| 国产精品人妻久久久影院| 免费观看的影片在线观看| 国产色婷婷99| 亚洲电影在线观看av| 少妇丰满av| 日韩欧美精品免费久久| 国产在视频线在精品| 亚洲第一区二区三区不卡| 亚洲av.av天堂| 亚洲精品国产av成人精品| 午夜久久久久精精品| 成人毛片60女人毛片免费| 国产亚洲91精品色在线| 国产午夜精品论理片| 欧美日韩亚洲高清精品| 黄色配什么色好看| 色哟哟·www| 成人毛片60女人毛片免费| 精品久久国产蜜桃| 99久国产av精品| 熟女电影av网| 免费观看的影片在线观看| 亚洲经典国产精华液单| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 国产成年人精品一区二区| 晚上一个人看的免费电影| 欧美日韩在线观看h| 能在线免费观看的黄片| 中文字幕av成人在线电影| 少妇裸体淫交视频免费看高清| 97超碰精品成人国产| 久久精品夜夜夜夜夜久久蜜豆| a级毛色黄片| 亚洲av中文av极速乱| 欧美成人一区二区免费高清观看| 国产精品一区二区三区四区免费观看| 一级毛片久久久久久久久女| 国产伦在线观看视频一区| 亚洲自偷自拍三级| 亚洲在线观看片| 亚洲欧洲国产日韩| 欧美成人午夜免费资源| 神马国产精品三级电影在线观看| 日本黄色片子视频| 白带黄色成豆腐渣| 午夜福利在线观看吧| 亚洲色图av天堂| 欧美一区二区亚洲| 免费看光身美女| 高清av免费在线| 精品一区二区三卡| 美女黄网站色视频| 中文字幕久久专区| 国产成人免费观看mmmm| 国产亚洲av片在线观看秒播厂 | 亚洲国产欧美在线一区| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 免费av不卡在线播放| 美女脱内裤让男人舔精品视频| 欧美三级亚洲精品| 男女国产视频网站| 舔av片在线| 建设人人有责人人尽责人人享有的 | 亚洲精品久久久久久婷婷小说| 一区二区三区免费毛片| 99re6热这里在线精品视频| 噜噜噜噜噜久久久久久91| 青春草视频在线免费观看| 一级a做视频免费观看| 日韩欧美一区视频在线观看 | 一级av片app| 街头女战士在线观看网站| 人妻少妇偷人精品九色| 国产黄色免费在线视频| 成人av在线播放网站| 国产av码专区亚洲av| 欧美激情在线99| 久久亚洲国产成人精品v| 高清av免费在线| 美女被艹到高潮喷水动态| 色哟哟·www| 亚洲在线观看片| 午夜激情久久久久久久| 尤物成人国产欧美一区二区三区| 身体一侧抽搐| 国产精品熟女久久久久浪| 麻豆久久精品国产亚洲av| 一本一本综合久久| av.在线天堂| 一本一本综合久久| 亚洲国产高清在线一区二区三| 毛片女人毛片| 在线a可以看的网站| videossex国产| 美女脱内裤让男人舔精品视频| 一个人观看的视频www高清免费观看| 成年版毛片免费区| 亚洲国产精品成人综合色| 国产精品一区二区在线观看99 | 日韩欧美一区视频在线观看 | 成人二区视频| 精品国产三级普通话版| 成人亚洲精品av一区二区| 欧美bdsm另类| 特大巨黑吊av在线直播| 色综合亚洲欧美另类图片| 蜜桃亚洲精品一区二区三区| 亚洲国产av新网站| 午夜久久久久精精品| 亚洲欧美成人综合另类久久久| 国产乱人偷精品视频| 久久久久精品久久久久真实原创| 久久精品国产鲁丝片午夜精品| 国产一区二区在线观看日韩| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 免费看日本二区| av国产免费在线观看| 又爽又黄无遮挡网站| 日本色播在线视频| 国产一级毛片在线| 免费av毛片视频| 晚上一个人看的免费电影| 亚洲成人精品中文字幕电影| 亚洲精品乱码久久久v下载方式| 亚洲av中文av极速乱| 女人被狂操c到高潮| 国产免费福利视频在线观看| 男人和女人高潮做爰伦理| 日韩欧美 国产精品| 午夜爱爱视频在线播放| 亚洲乱码一区二区免费版| a级毛色黄片| 97人妻精品一区二区三区麻豆| 国产精品久久久久久精品电影小说 | 色视频www国产| 一区二区三区高清视频在线| 麻豆av噜噜一区二区三区| 国产一区二区在线观看日韩| 国产成人精品一,二区| 亚洲av日韩在线播放| 久久久久性生活片| 国产av在哪里看| 乱系列少妇在线播放| 人妻制服诱惑在线中文字幕| 色综合色国产| 国产伦精品一区二区三区视频9| 亚洲国产色片| 在线天堂最新版资源| 日本免费a在线| 老女人水多毛片| 亚洲av男天堂| 午夜视频国产福利| 国产精品女同一区二区软件| 久久鲁丝午夜福利片| 国产老妇伦熟女老妇高清| 大陆偷拍与自拍| 91av网一区二区| 在线免费观看的www视频| 久久久久精品性色| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 在线观看免费高清a一片| 国产v大片淫在线免费观看| 男的添女的下面高潮视频| 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 在线观看人妻少妇| 91久久精品国产一区二区三区| 国产伦理片在线播放av一区| 能在线免费看毛片的网站| 国产老妇女一区| 亚洲成人久久爱视频| 亚洲精品日韩av片在线观看| 国产91av在线免费观看| 亚洲国产欧美在线一区| 国产色爽女视频免费观看| 有码 亚洲区| 看免费成人av毛片| 欧美日韩国产mv在线观看视频 | videos熟女内射| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久久成人| 精品久久久久久久末码| 欧美区成人在线视频| 三级毛片av免费| 黄片wwwwww| 高清毛片免费看| 男女边吃奶边做爰视频| 联通29元200g的流量卡| a级毛片免费高清观看在线播放| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 欧美xxⅹ黑人| 婷婷色av中文字幕| 看非洲黑人一级黄片| 国产色婷婷99| 国产av不卡久久| a级一级毛片免费在线观看| 最近中文字幕高清免费大全6| 免费无遮挡裸体视频| 乱系列少妇在线播放| 看非洲黑人一级黄片| 精品久久久久久久久久久久久| a级毛片免费高清观看在线播放| 国产又色又爽无遮挡免| 卡戴珊不雅视频在线播放| 人妻夜夜爽99麻豆av| 在线免费观看不下载黄p国产| 久久精品久久久久久噜噜老黄| 亚洲图色成人| 免费av毛片视频| 精品人妻一区二区三区麻豆| 亚洲欧美成人综合另类久久久| 一级片'在线观看视频| 狂野欧美白嫩少妇大欣赏| av.在线天堂| 婷婷色综合www| 亚洲成人av在线免费| 非洲黑人性xxxx精品又粗又长| 国产精品无大码| 丝瓜视频免费看黄片| 看非洲黑人一级黄片| 中文字幕人妻熟人妻熟丝袜美| 天堂√8在线中文| 最近最新中文字幕大全电影3| 亚洲国产欧美在线一区| 黄色配什么色好看| 水蜜桃什么品种好| 国产真实伦视频高清在线观看| 国产色爽女视频免费观看| 精品午夜福利在线看| 国产探花极品一区二区| 国产精品久久久久久精品电影小说 | 男插女下体视频免费在线播放| 两个人的视频大全免费| 亚洲久久久久久中文字幕| 国产成人福利小说| 国产探花极品一区二区| 能在线免费看毛片的网站| av在线蜜桃| 成年av动漫网址| 男女国产视频网站| 国产探花极品一区二区| 丰满少妇做爰视频| 国内少妇人妻偷人精品xxx网站| 欧美3d第一页| 自拍偷自拍亚洲精品老妇| 如何舔出高潮| 亚洲精品色激情综合| 成人午夜精彩视频在线观看| 少妇裸体淫交视频免费看高清| .国产精品久久| 成年免费大片在线观看| 欧美精品国产亚洲| 97超碰精品成人国产| 毛片女人毛片| 免费看a级黄色片| 亚洲精品中文字幕在线视频 | 性插视频无遮挡在线免费观看| 国产女主播在线喷水免费视频网站 | 三级男女做爰猛烈吃奶摸视频| 一级毛片我不卡| 街头女战士在线观看网站| 乱码一卡2卡4卡精品| 国产高清国产精品国产三级 | 免费看不卡的av| 一夜夜www| 日韩av不卡免费在线播放| 如何舔出高潮| 国产成人精品婷婷| 波多野结衣巨乳人妻| 国产av国产精品国产| 国产黄色免费在线视频| 国产色爽女视频免费观看| 一级毛片电影观看| 国产有黄有色有爽视频| 在现免费观看毛片| 国产午夜精品久久久久久一区二区三区| 国产成人a区在线观看| 亚洲18禁久久av| 亚洲成色77777| 亚洲成人一二三区av| 免费播放大片免费观看视频在线观看| 夜夜爽夜夜爽视频| 久久午夜福利片| 国产乱人视频| 免费黄色在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人妻系列 视频| 亚洲精品,欧美精品| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄|