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

    海洋溫鹽度資料多變量同化研究進展

    2013-08-15 00:46:26張春玲李宏許建平王振峰
    海洋預報 2013年1期
    關鍵詞:變分鹽度約束

    張春玲,李宏,許建平,2,王振峰

    (1.國家海洋局第二海洋研究所,浙江杭州 310012;2.衛(wèi)星海洋環(huán)境動力學國家重點實驗室,浙江杭州 310012;3.中國海洋大學海洋環(huán)境學院,山東青島 266003;4.東海艦隊司令部海洋水文氣象中心,浙江寧波 312122)

    1 引言

    至少在Argo計劃出現(xiàn)以前,海洋鹽度觀測資料比較匱乏,海溫觀測資料相對豐富些;另一方面,在平均狀態(tài)下,由溫度引起的密度變化要遠大于由鹽度引起的密度變化[1-2],這使得人們在研究一般海洋環(huán)流特征時,通常忽略鹽度的影響。所以早期的海洋資料同化主要考慮單個變量即海溫資料的同化問題,而鹽度和流場在同化中保持不變,僅僅是通過動力學模式來調(diào)整。

    近些年來,人們研究發(fā)現(xiàn)鹽度變化對密度場和海洋環(huán)流的影響是不可忽略的,單變量(海溫)的數(shù)據(jù)同化有時會嚴重地惡化密度場,導致模式計算的流場比沒有同化觀測資料時還要差[3],所以有關鹽度的變化對全球海洋的影響成了人們非常關注的問題。實際上,對于一致的模式動力學來說,溫度和鹽度應該同時被監(jiān)測,尤其在西太平洋海域,這一點顯得更為重要[4]。鹽度垂直結(jié)構對于西太平洋熱量再分布具有重大的影響,鹽度變化對海水層化的改變可以影響西太平洋地區(qū)上層海洋的熱量和動量收支[5-6];在太平洋暖池,由鹽度控制下的熱傳輸代表了影響此區(qū)域熱預算的一個重要過程;并且,近些年的研究表明,只修正溫度的單變量同化系統(tǒng)是有缺陷的,如,1996年美國國家環(huán)境預報中心西太平洋赤道海區(qū)的業(yè)務化海洋分析結(jié)果中,海面動力高度就有5—10 cm的誤差,這很可能是由于缺乏鹽度實測資料而忽略鹽度變化引起的[7];Troccoli等[8]發(fā)現(xiàn)當只同化溫度時,要得到一個好的溫度同化場就必須正常地調(diào)整鹽度廓線。這些研究都表明了在海洋資料同化中糾正鹽度的重要性。

    為了彌補鹽度觀測資料的不足,解決同化中只有溫度沒有鹽度的問題,需要把鹽度和溫度或海面高度關聯(lián)起來,比如可以引入溫、鹽度之間的約束關系。因此,許多學者針對一些目前被廣泛應用的資料同化方法,例如最優(yōu)插值(OI)、集合卡爾曼濾波(EnKF)、集合最優(yōu)插值(EnOI)和三維變分(3D-VAR)等,提出了不同的溫、鹽度多變量調(diào)整方案,力求在同化過程中,能夠保證溫、鹽度場均能得到調(diào)整。

    2 用于多變量分析的幾種主要的海洋資料同化方法

    目前的資料同化方法根據(jù)其理論可分為兩類,一類是基于統(tǒng)計估計理論的,如最優(yōu)插值、卡爾曼濾波、集合卡爾曼濾波、集合最有插值等;另外一類是基于變分方法的,如三維變分,強約束四維變分和弱約束四維變分等。這里簡單回顧廣泛用于多變量分析的海洋資料同化方法的基本原理,主要包括:最優(yōu)插值、集合卡爾曼濾波、集合最優(yōu)插值及三維變分等四種。

    2.1 最優(yōu)插值法(OI)

    Gandin[9]首次獨立地推出了多元OI方程組,并把它們應用于在前蘇聯(lián)進行的客觀分析。在OI方法中,分析場是背景場與由權重矩陣加權的修正量之和,采用最小二乘方法求得最佳線性無偏估計方程中的最優(yōu)權重矩陣。此方法的基本假設是,對于一個模型變量,在確定它的增量時,只有幾個觀測值是重要的。因此,OI易于編碼,并且計算量相對較小,這是它的主要優(yōu)點。但基于這一假設,也使得OI分析結(jié)果并非全局最優(yōu),分析在空間上不協(xié)調(diào)[10],并且,OI所用協(xié)方差矩陣是固定的,不隨時間變化,這就限制了它不能將動力模式和觀測信息很好地融合在一起;另外,OI是針對線性系統(tǒng)發(fā)展起來的,難以處理觀測算子非線性的情況;無法確保大小尺度分析的一致性。

    2.2 集合卡爾曼濾波(EnKF)

    1994年Evensen[11]在Kalman濾波方法[12]的基礎上,發(fā)展了基于蒙特卡羅算法的集合卡爾曼濾波方法。它結(jié)合了Kalman濾波和集合預報的優(yōu)點,即用有限的集合樣本來估算誤差協(xié)方差矩陣的不確定性,這樣計算量明顯減小。在系統(tǒng)為線性,且樣本數(shù)量趨于無窮時,EnKF和KF是等價的[13]。由于其概念簡單,不需要作線性假設,不需要求解模式的切線性及其伴隨,適合于并行計算等優(yōu)點,是一個目前比較流行的方法。但由于EnKF是通過選取有限的樣本來構造背景場誤差協(xié)方差,這勢必使得樣本集合離散度不夠(樣本量有限),產(chǎn)生樣本誤差問題;在實際操作中,系統(tǒng)的非線性性,及通常利用擾動觀測法獲取樣本初值,使得這種樣本誤差問題更為明顯;另外,EnKF計算量依然很大。

    2.3 集合最優(yōu)插值(EnOI))

    Evensen[14-16]將集合思想吸收到最優(yōu)插值同化技術中,提出了集合最優(yōu)插值數(shù)據(jù)同化法:格點的分析值在一個固定的模式向量樣本集合(如長時間序列的模式積分)空間內(nèi)進行計算,模式統(tǒng)計誤差不隨時間變化,從而減少計算量。但在誤差計算時仍沿用EnKF集合預報的方式,以獲得較傳統(tǒng)OI方法更優(yōu)的分析值。EnOI能夠保持準動力一致性,避免假設均勻和各向同性等[17],并且EnOI的計算量要比EnKF小得多。但由于模式誤差不隨模式積分時間改變,EnOI較EnKF得到一個次優(yōu)解。我國有學者[18-19]將此方法與其他同化方法進行了對比同化試驗,EnOI是一種嶄新的數(shù)據(jù)同化方法,近年來逐漸被國內(nèi)外學者廣泛研究應用[20-22]。

    2.4 三維變分(3D-VAR)

    三維變分基于極大似然估計的理論基礎,通過求解一個目標函數(shù)(也稱代價函數(shù))的極小值(一般利用目標函數(shù)的梯度求其極小值),產(chǎn)生一個分析時刻的綜合考慮背景場和觀測值的大氣或海洋真實狀態(tài)的最大似然估計,并且給出背景場和觀測場各自相應的精度。3D-VAR進行的是三維空間的全局分析,避免了分析不是全局最優(yōu)的問題;也可以處理觀測算子是非線性的情況,這樣可以同化各種不同來源的觀測資料[23-24]。但此方法是在某一時刻進行的分析,前一時刻的同化結(jié)果可作為后一時刻模式運行的初始場。但在使用時,無法用后面時刻的資料來訂正前面的結(jié)果,同化的解在時間上不連續(xù)。

    3 不同的溫、鹽度多變量調(diào)整方案

    單純的溫度資料同化只能夠調(diào)整溫度場,而鹽度場得不到任何訂正,即便是溫、鹽度兩種觀測資料同時進入同化系統(tǒng)中,獨立地訂正溫度場、鹽度場,而不考慮溫、鹽度之間的約束關系,也往往會帶來虛假信息。那么在缺乏鹽度直接觀測的情況下,如何來估計和調(diào)整鹽度場成為海洋學家們備受關注的問題。目前,國內(nèi)外學者提出了一些能夠同時調(diào)整溫度和鹽度的方案,并將其應用于上述海洋資料同化方法中。

    3.1 基于溫-鹽(T-S)關系估計鹽度場

    海水溫度和鹽度是反映海水物理性質(zhì)的兩個基本變量,某一地區(qū)的水團特性通常對應某種特定的溫鹽關系。Stommel[25]基于從氣候態(tài)數(shù)據(jù)集得到的溫度與鹽度之間的統(tǒng)計相關性,即溫-鹽(T-S)關系,首次提出了由溫度觀測來估計鹽度的思想。后來,許多學者[26-27]改進了這種思想方法,以擴展其應用范圍。

    這種方法的基本思想是,當溫度分析場在單變量溫度資料同化中獲取后,通過局地T-S關系,計算出對應的鹽度分析場,而當有鹽度資料時,再同化鹽度資料,進而修正T-S關系。此方案基于兩個基本假設:一是僅考慮海水的溫度和鹽度隨深度的變化關系,并且溫、鹽度垂直剖面數(shù)據(jù)的物理(或統(tǒng)計)相關性(T-S關系)在一定時間和區(qū)域內(nèi)是保持不變的,Troccoli等[28-29]的研究指出,在熱帶海域,這種T-S關系能夠保持2到4周;二是假定海洋混合層深度對溫度和鹽度來說是一致的。第一個基本假設是該方法的核心。

    構造一種較好的T-S關系是此方法的關鍵問題,Reynolds等[30]曾指出,如果我們能夠選取一種表現(xiàn)形式來有效刻畫海表鹽度觀測信息,那么我們可以借助這一信息來改善上層海洋500m的鹽度垂直剖面。最初有學者從氣候態(tài)溫、鹽度資料通過函數(shù)擬合來構造T-S關系,但氣候態(tài)資料一般是長時期的平均態(tài)資料,不能捕捉到季節(jié)變化特征,且不能充分構建出密度場的變化特性[7],通常情況下,這是一種可行但比較粗略的方法。因此,更為可靠的是利用時間和空間上更為連續(xù)的溫、鹽度場來構造這一關系,Troccoli等[31]建議利用模式預報溫、鹽度場,而韓桂軍等則是利用就近的CTD溫、鹽度觀測資料。韓桂軍等[32]借鑒Troccoli的鹽度調(diào)整方案,利用多重網(wǎng)格三維變分海洋數(shù)據(jù)同化方法,開發(fā)了中國近海及鄰近海域海面高度和三維溫鹽流的23年海洋再分析產(chǎn)品。

    然而,T-S關系并不是一個很精確的假設,而且在海洋非等熵動力過程存在的地方(比如混合層,河流入??诘葏^(qū)域),溫度和鹽度垂直剖面會變得高度不相關。此時,利用T-S關系根據(jù)溫度場來構造鹽度場就不太可取,對海洋混合層而言,額外的海表鹽度觀測值的加入能夠彌補T-S約束在此時的不適用性。在這些特殊區(qū)域,對鹽度一般不作調(diào)整[7]。

    3.2 利用耦合溫-鹽EOF方法構造鹽度場

    通過溫-鹽關系來調(diào)整鹽度場,雖然簡單易行,但在某些海區(qū)T-S關系具有很大的不確定性。Maes[1]呼吁要關注T-S關系隨時間的變化,并建議用一種新的途徑來估計上層鹽度剖面數(shù)據(jù):將以溫-鹽函數(shù)關系為基礎的經(jīng)驗正交函數(shù)(Empirical Orthogonal Functions,EOFs)與海表鹽度及海表動力高度數(shù)據(jù)相結(jié)合來估計溫鹽廓線。通過用原始的獨立數(shù)據(jù)集進行誤差估計表明,與利用傳統(tǒng)的溫-鹽關系調(diào)整鹽度場的方法相比,該方法在一定程度上是成功的。此方法是通過把觀測到的海面變率、海表溫度和海表鹽度從諸多模態(tài)中分離出來實現(xiàn)的,充足的溫、鹽度觀測資料是這種方法的基礎。

    在熱帶太平洋海域,海洋觀測資料較其他地區(qū)相對充足,Maes等[33]在此海區(qū),使用歷史CTD資料,先將一段時期的溫、鹽度剖面資料按區(qū)域進行劃分,在選取的每個小區(qū)域中,構造耦合的溫-鹽EOF模態(tài),并選取幾個主要的模態(tài)作為重新構造溫、鹽度變化場的基函數(shù),再現(xiàn)溫、鹽度變率,通過求權重函數(shù)最小值來確定系數(shù),并在權重函數(shù)中加入鹽度約束項,以保證鹽度值在合理的范圍之內(nèi);并通過增添高度計數(shù)據(jù)拓展了上述方法[29];Fujii等[34]也用同樣的方法,借助于變分同化方法重構了日本東部海域的溫、鹽度分布場;另外,Bellucci等[35]借鑒Maes提出的耦合溫-鹽EOF模態(tài),并針對鹽度同化的困難,利用耦合溫-鹽來構建背景場誤差協(xié)方差,同時考慮背景場協(xié)方差的時空結(jié)構對鹽度的影響,發(fā)展了一個多變量降秩最優(yōu)插值同化系統(tǒng)。

    利用耦合溫-鹽EOF方法來構造鹽度場可以將不同的海洋觀測信息加入到權重函數(shù)中作為約束條件,從而提高精確度,并能根據(jù)我們重點關心的問題來選擇對應的EOF信號。但耦合溫-鹽EOF方法要求溫、鹽度剖面數(shù)據(jù)越多的時候,其反映的信號越全面真實,由于海洋觀測資料空間分布的不均勻,此方法難以推廣到全球海洋。

    3.3 借助海面動力高度訂正鹽度場

    海面動力高度也是海洋中一個很重要的變量,其能捕捉到海洋中固定參考層上的所有斜壓過程,這些斜壓過程與溫度、鹽度的變化有著密切聯(lián)系。因此,在海洋資料同化中,海面高度信息對鹽度場的訂正有著積極的意義。Vossepoel等[36-37]提出了一種基于海面動力高度信息的鹽度調(diào)整方案,并將這一鹽度訂正方案融入到三維變分同化系統(tǒng)中,在常規(guī)三維變分代價函數(shù)的基礎上加入了與海面動力高度相關的溫、鹽度訂正項,以此來強迫海面動力高度對溫、鹽度的一致性約束。該方案先利用溫度觀測資料和氣候態(tài)T-S關系獲取鹽度估計值,然后根據(jù)由溫、鹽度觀測值計算的動力高度觀測值與由溫度觀測值和鹽度估計值計算的動力高度估計值之差來計算鹽度訂正值。在此基礎上,閆長香等[38]考慮了背景誤差協(xié)方差矩陣的垂向相關,利用海面動力高度觀測數(shù)據(jù)來估計溫、鹽度垂向剖面。目前,海面動力高度觀測資料較為豐富:繼20世紀90年代成功實施了Topex/Poseidon衛(wèi)星高度計觀測計劃后,Jason-1海面高度觀測衛(wèi)星(Jason-1 Mission)也在2001年12月發(fā)射升空,這個觀測系統(tǒng)可以提供幾乎全球覆蓋(范圍為66°N-66°S)、時間分辨率為10 d的海面高度觀測資料。因此利用這一方法來構造溫、鹽度場,可以彌補某些區(qū)域資料的不足,也可作為訂正溫、鹽度資料的一個參考方法。

    3.4 基于平衡約束的溫、鹽度多變量調(diào)整方案

    平衡約束是另外一種可以訂正鹽度的方案。在有溫、鹽度兩個模式變量的最優(yōu)插值或三維變分同化方案中,背景場誤差協(xié)方差是一個至關重要的量,其決定了不同變量之間的空間分布信息。如果不考慮溫、鹽度之間的任何約束,那么背景誤差協(xié)方差矩陣就是一個塊對角矩陣,而實際海洋中溫鹽關系是存在的。Derber等[39]提出了一種線性平衡約束方案,假設溫度及鹽度的背景場與真實場的偏差之間存在線性約束關系,引入一個約束溫度和鹽度關系的線性平衡算子,從而保證溫、鹽度場之間具有相關性。此時,僅同化溫度或鹽度資料時,可以將溫度和鹽度間的信息相互傳遞,使得溫、鹽度場都能同時得到調(diào)整。在此基礎上,Ricci等[2]沿用Troccoli等[7]的方式,采用了一種代表局地背景場溫、鹽度約束關系的方案:假設在任何分析格點,鹽度可表示為一個背景場溫度的函數(shù),并且此函數(shù)是可微分的,鹽度擾動量可以根據(jù)溫度擾動量來計算。并鑒于溫-鹽關系在海洋混合層難以得到保持,引進一個控制系數(shù)來控制不同垂直深度上的溫度和鹽度擾動量之間的關系,使得在混合層等一些特殊區(qū)域,對鹽度不作調(diào)整。

    線性約束方案的一個優(yōu)點是,單一的觀測資料進入同化系統(tǒng)時能夠?qū)?、鹽度資料進行同步調(diào)整。然而在海洋中,平衡約束可以是非線性的,甚至是強非線性的。

    韓桂軍等[31]人的研究表明,一些溫-鹽關系具有高階多項式(如四階)的形式,并提出了一個把溫-鹽關系作為弱約束的3D-VAR同化方案,利用歷史CTD資料,通過多項式擬合來構造局地溫-鹽關系,給出了包含溫、鹽度兩個分析變量的垂向一維代價函數(shù),其最后一項包括非線性的溫鹽約束關系,并通過代價函數(shù)最小化同時求得溫度和鹽度的最優(yōu)分析場。基于三維變分中Derber[38]提出的線性平衡約束形式,朱江[23]針對那些作用于模式變量(而不是增量變量)上的約束,導出背景場很好的滿足平衡關系和背景場不必滿足平衡約束兩種情況下的非線性形式,提出了在三維變分同化中考慮非線性平衡約束的方案;并利用該方案,進行了從海面動力高度來估計溫、鹽度的資料同化數(shù)值試驗。結(jié)果表明,采用了非線性平衡約束可以提高同化的質(zhì)量;如果背景場的平衡關系滿足不好的話,第二種形式將產(chǎn)生較第一種形式更好的結(jié)果。朱江等[40]介紹了新完成的一個海洋資料三維變分同化系統(tǒng)OVALS(Ocean Variational Analysis System)的設計方案,并利用熱帶太平洋的大洋環(huán)流模式進行了實際海洋溫度、鹽度和衛(wèi)星高度計資料21a的同化實驗,對此系統(tǒng)的性能進行了檢驗。結(jié)果表明,經(jīng)過同化后,溫度和鹽度場較非同化試驗產(chǎn)生了顯著的改善。OVALS的設計方案考慮了背景誤差的垂向相關和非線性的溫-鹽關系,通過同化高度計資料來直接調(diào)整模式的溫度和鹽度場。就當前的觀測系統(tǒng)而言,OVALS可以說是一個比較全面的同化系統(tǒng),可以對多種觀測資料進行同化,并對多變量進行分析。

    4 小結(jié)與討論

    隨著研究的加深,這些不同的鹽度調(diào)整方案將對大型海洋數(shù)值模式及同化系統(tǒng)的建立起到積極的推動作用。但任何溫鹽約束調(diào)整方案都無法完全解決由單變量同化所帶來的一系列問題,必須依靠直接同化溫度和鹽度觀測資料來解決這些問題,特別是次表層鹽度觀測資料的同化,而鹽度觀測在時空分布上的嚴重不足是阻礙直接對其同化的一個重要因素。

    2000年底正式啟動國際Argo計劃,其觀測目標是能取得世界大洋中精確度分別為±0.005℃和±0.01的海水溫度和鹽度資料[41]。由3000個自動剖面觀測浮標(簡稱Argo浮標)組成的全球Argo實時海洋觀測網(wǎng)也已于2007年10月末初步建成。截止到2010年12月底,國際Argo計劃在全球海洋0—2000m深度范圍內(nèi)累計獲得的溫、鹽度剖面已達70余萬條,并正以每年10萬條剖面以上的速度在增加。Argo剖面浮標觀測資料以其水平分布廣,觀測深度深、層次密以及數(shù)據(jù)量大和測量精度高的優(yōu)勢,在很大程度上改善了海洋鹽度觀測資料匱乏的現(xiàn)狀。由2011年5月召開的第六屆全國海洋資料同化研討會的情況來看,隨著Argo資料的逐漸積累,鹽度資料相對匱乏的問題已經(jīng)得到改善,使直接進行海洋溫鹽同步同化成為可能。

    但由于Argo資料空間分布不均勻,而且缺少表層(0 m層)觀測,一般海洋模式用戶在直接使用時存在較大困難,從而使得豐富的Argo溫、鹽度觀測資料不能充分發(fā)揮其優(yōu)勢。而2005年正式開始實施的高分辨率海面溫度計劃(high-resolution sea surface temperature Pilot Project,簡稱GHRSST-PP),針對業(yè)務化海洋學和氣候研究與預測,可提供全球覆蓋的高時空分辨率的海面溫度觀測資料,其空間分辨率為10 km(部分地區(qū)可以達到2 km),時間分辨率可達到6 h[42]??v觀其余如衛(wèi)星遙感資料的應用過程,其獲得普遍使用主要是解決了網(wǎng)格化和等時間間隔數(shù)據(jù)集的制作問題,從而便于模式引用。因此,將Argo資料與衛(wèi)星遙感SST數(shù)據(jù)融合,盡快推出適用于海洋模式的Argo數(shù)據(jù)產(chǎn)品或網(wǎng)格數(shù)據(jù)集,是進一步擴大和推進Argo數(shù)據(jù)應用的重中之重。

    [1]Maes C.A note on the vertical scales of temperature and salinity and their signature in dynamic height in the western Pacific Ocean:Implications for data assimilation[J].J.Geophys.Res,1999,104:15575-15585.

    [2]Ricci S A,Weaver J,Vialard,Rogel P.Incorporating State dependent temperature salinity constraints in the background error covariance of variational ocean data assimilation.Mon[J].Wea.Rev.2005,133:317-338.

    [3]Cooper N S.The effect of salinity in tropical ocean models[J].J.Phys.Oceanogr,1988,18:1333-1347.

    [4]Woodgate R A.Can we assimilate temperature data alone into a full equation of state model[J].Ocean Model,1997,114:4-5.

    [5]Vialard J,Delecluse P.An OGCM study for the TOGA Decade.Part I:Role of salinity in the physics of the western Pacific fresh pool[J].J Phys.Oceanogr,1998a,28:1071-1088.

    [6]Vialard J,Delecluse P.An OGCM study for the TOGA decade.Part II:Barrier layer formation and variability[J].J.Phys.Oceanogr,1998b,28:1089-1106.

    [7]Troccoli A,Coauthors.Salinity adjustments in the presence of temperature data assimilation[J].Mon.Wea.Rev,2002,130:89-102.

    [8]Ji M,Reynolds R W,Behringer D.Use of TOPEX/Poseidon sea level data for ocean analysis and ENSO prediction:Some early results[J].J Clim 2000,13:216-231.

    [9]Gandin L S.Objective analysis of meteorological fields[M].Gidromet:Leningrad,1963.1-242.

    [10]Kalnay E.Atmospheric Modeling,Data Assimilation and Predictability[M].Cambridge:Cambridge University Press,2003.1-343.

    [11]Evensen G.Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics[J].J Geophys Res,1994,99:10143-10162.

    [12]Kalman R E,Bucy R.New results in linear filtering and predication theory[J].Journal of Basic Engineering,1961,83:95-108.

    [13]Houtekamer P L,Mitchell H L.Ensemble Kalman filtering[J].Meteorol.Soc,2005,131:3269-3289.

    [14]Evensen G.The Ensemble Kalman Filter:theoretical formulation and practical implementation[J].Ocean Dyn.,2003,53:343-367.

    [15]Oke P,Brassington G B,et al.The bluelink ocean data assimilation system(BODAS)[J].Ocean Model,2008,21:46-70.

    [16]Counillon F,Bertino L.Ensemble Optimal Interpolation:multivariate properties in the Gulf of Mexico.Tellus-A,2009,61:296-308.

    [17]Oke P,Brassington G B,et al.Ocean data assimilation:a case for ensemble optimal interpolation[J].Australian Meteorological and Oceanographic Journal,2010,59:67-76.

    [18]Fu W,Zhu J,Yan C.A comparison between 3DVAR and EnOI techniques for satellite altimetry data assimilation[J].Ocean Model,2009,26:206-16.

    [19]Wan L,Bertino L,Zhu J.Assimilating Altimetry Data into a HYCOM Model of the Pacific:Ensemble Optimal Interpolation versus Ensemble Kalman Filter[J].Journal of Atmospheric and Oceanic Technology,2010,27:753-765.

    [20]Bertino L,Lisaeter K A.The TOPAZ monitoring and prediction system for the Atlantic and Arctic Oceans[J].Operat.Oceanogr,2008,2:15-18.

    [21]Xie J,Zhu J.Ensemble optimal interpolation schemes for assimilating Argo profiles into a hybrid coordinate ocean model[J].Ocean Model,2010,33:283-298.

    [22]Yan Changxiang,Zhu Jiang,Xie Jiping.An ocean reanalysis system for the joining area of Asia and Indian-Pacific ocean[J].Atmospheric and Oceanic Science Letters,2010,3(2):81-86

    [23]朱江,閆長香.三維變分資料同化中的非線性平衡約束[J].中國科學D輯:地球科學,2005,35(12):1187-1192.

    [24]高山,王凡,李明.中尺度渦的高度計資料同化模擬[J].中國科學D輯:地球科學,2007,37(12):1669-1678.

    [25]Stommel H M.Note on the use of T-S correlation for dynamic height anomaly computation[J].J Mar Res,1947,6:85-92.

    [26]Flierl G R.Correcting expendable bathythermograph(XBT)data for salinity effects to compute dynamic heights in Gulf Stream rings[J].Deep-Sea Res I,1978,25:129-134.

    [27]Kessler W S,Taft B A.Dynamic heights and zonal geostrophic transports in the central tropical Pacific during 1979-84[J].J.Phys.Oceanogr,1987,17:97-122.

    [28]Emery W J,Dewar L J.Mean temperature-salinity,salinity-depth,and temperature-depth curves in the North Atlantic and North Pacific[J].Progress in Oceanography,1982,16:219-305.

    [29]Maes C,Behringer D.Using Satellite-derived Sea Level and Temperature Profiles for Determining the Salinity Variability:A NewApproach[J].J.Geophys Res,2000b,105:8537-8547.

    [30]Reynolds R W,Ji M,Leetmaa A.Use of salinity to improve ocean modeling[J].Phys Chem Earth,1998,23:545-555.

    [31]Han G,Zhu J,Zhou G.Salinity estimation using the T-S relation in the context of variational data assimilation[J].J Geophys Res,2004,109,C03018,doi:10.1029/2003 JC 001781.

    [32]韓桂軍,李威,張學峰,等.中國近海及鄰近海域海洋再分析技術報告[R].天津:國家海洋信息中心,2009,1-6.

    [33]Maes C,Behringer D,Reynolds R W.Retrospective analysis ofthe salinity variability in the western Tropical Pacific Ocean using an indirect minimization approach[J].J.Atmos Ocean Technol,2000a,17:512—524.

    [34]Fujii Y,Kamachi M.A reconstruction of observed profiles in the Sea East of Japan using vertical coupled temperature-salinity EOF modes[J].J Oceanogr,2003,59:173-186.

    [35]Bellucci A,Masina S,DiPietro P,Navarra A.Using temperature-salinity relations in a global ocean implementation of a multivariate data assimilation scheme[J].Mon.Weather Rev.2007,135:3785-3807.

    [36]Vossepoel F,Reynolds R W,Miller L.The use of sea level observations to estimate salinity variability in the tropical Pacific[J].J AtmosOceanicTechnol,1999,16:1400-1414.

    [37]Vossepoel F,Behringer D.Impact of sea level assimilation on salinity variability in the western equatorial Pacific[J].J Phys Oceanogr,2000,30:1706-1721.

    [38]Yan Changxiang,Jiang Zhu,Rongfeng Li,Guangqing Zhou.Roles of vertical correlations of background error and T-S relations in estimation of temperature and salinity profiles from sea surface dynamic height[J].J.Geophys.Res.,2004,109,C08010,doi:10.1029/2003JC002224.

    [39]Derber J,Bouttier F.A reformulation of the background error covariance in ECMWF global data assimilation[J].Tellus,1999,51:195-221.

    [40]朱江,周廣慶,閆長香,等.一個三維變分海洋資料同化系統(tǒng)的設計和初步應用[J].中國科學D輯:地球科學,2007,37(2):261-271.

    [41]許建平,劉增宏.中國Argo大洋觀測網(wǎng)試驗[M].北京:氣象出版社,2007:4-5.

    [42]GHRSST Porject Office.The GHRSST-PP Development and Im-plementation Plan[R].UK:Met Office,2003.1-6.

    猜你喜歡
    變分鹽度約束
    “碳中和”約束下的路徑選擇
    逆擬變分不等式問題的相關研究
    約束離散KP方程族的完全Virasoro對稱
    求解變分不等式的一種雙投影算法
    關于一個約束變分問題的注記
    一個擾動變分不等式的可解性
    鹽度和pH對細角螺耗氧率和排氨率的影響
    鹽度脅迫對入侵生物福壽螺的急性毒性效應
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    適當放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    久久久久久伊人网av| 夜夜夜夜夜久久久久| 直男gayav资源| av专区在线播放| 国产精品99久久久久久久久| 欧美一区二区精品小视频在线| eeuss影院久久| 欧美性感艳星| 一级av片app| 日日夜夜操网爽| 波野结衣二区三区在线| 久久精品夜夜夜夜夜久久蜜豆| 超碰av人人做人人爽久久| 久久久国产成人免费| 老女人水多毛片| 婷婷亚洲欧美| 在线观看一区二区三区| 一个人观看的视频www高清免费观看| 男人狂女人下面高潮的视频| 日本熟妇午夜| 午夜a级毛片| 69人妻影院| 国产中年淑女户外野战色| 成人综合一区亚洲| 久久久久久九九精品二区国产| 亚洲国产欧洲综合997久久,| 亚洲国产色片| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区色噜噜| 丰满的人妻完整版| 69人妻影院| 91精品国产九色| 亚洲无线观看免费| 美女免费视频网站| 亚洲精品影视一区二区三区av| 国产视频内射| 欧美激情国产日韩精品一区| av国产免费在线观看| 国内毛片毛片毛片毛片毛片| 特大巨黑吊av在线直播| 色综合色国产| 色噜噜av男人的天堂激情| 日本免费a在线| 在线看三级毛片| 精品久久久久久久人妻蜜臀av| av在线蜜桃| 我的老师免费观看完整版| 欧美另类亚洲清纯唯美| 亚洲国产色片| 午夜激情欧美在线| 日韩 亚洲 欧美在线| bbb黄色大片| x7x7x7水蜜桃| 亚洲av不卡在线观看| 亚洲欧美精品综合久久99| 99久久成人亚洲精品观看| 国产高清三级在线| 真人做人爱边吃奶动态| 日韩国内少妇激情av| 亚洲专区国产一区二区| 色哟哟哟哟哟哟| 久久人人精品亚洲av| 国内揄拍国产精品人妻在线| 琪琪午夜伦伦电影理论片6080| 成人av一区二区三区在线看| 精品人妻偷拍中文字幕| 亚洲精品成人久久久久久| 在线免费观看的www视频| 欧美性感艳星| 美女cb高潮喷水在线观看| 夜夜夜夜夜久久久久| 两人在一起打扑克的视频| 又爽又黄无遮挡网站| 干丝袜人妻中文字幕| 日本撒尿小便嘘嘘汇集6| 久久精品国产99精品国产亚洲性色| 亚洲天堂国产精品一区在线| 亚洲三级黄色毛片| 亚洲欧美日韩东京热| 永久网站在线| 成人三级黄色视频| 日日撸夜夜添| 黄色视频,在线免费观看| 两个人视频免费观看高清| 午夜精品在线福利| 91精品国产九色| 午夜免费男女啪啪视频观看 | 看黄色毛片网站| 国产精品野战在线观看| 美女 人体艺术 gogo| 中文字幕久久专区| 久久久久国产精品人妻aⅴ院| 亚洲精品一区av在线观看| 日本精品一区二区三区蜜桃| 最近最新免费中文字幕在线| 亚洲一区高清亚洲精品| 国产av麻豆久久久久久久| 全区人妻精品视频| av在线观看视频网站免费| 3wmmmm亚洲av在线观看| 色哟哟·www| 欧美黑人巨大hd| 热99re8久久精品国产| 亚洲国产精品久久男人天堂| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 中文字幕av成人在线电影| 草草在线视频免费看| 天堂√8在线中文| 色精品久久人妻99蜜桃| 国产单亲对白刺激| 看片在线看免费视频| 又紧又爽又黄一区二区| 无遮挡黄片免费观看| 美女高潮喷水抽搐中文字幕| 欧美日韩乱码在线| 久久久久久久精品吃奶| 国产私拍福利视频在线观看| 久久久精品欧美日韩精品| 1024手机看黄色片| 精品久久久久久成人av| 亚洲午夜理论影院| 免费在线观看影片大全网站| 一进一出好大好爽视频| 久久久久久久精品吃奶| 日本撒尿小便嘘嘘汇集6| 色尼玛亚洲综合影院| 中文资源天堂在线| 在线观看免费视频日本深夜| 久久人人精品亚洲av| 88av欧美| 丝袜美腿在线中文| 国产精品美女特级片免费视频播放器| 日本在线视频免费播放| netflix在线观看网站| 永久网站在线| 男女下面进入的视频免费午夜| 精品久久久久久久久av| 亚洲七黄色美女视频| 熟女人妻精品中文字幕| 亚洲在线自拍视频| 亚洲经典国产精华液单| av在线观看视频网站免费| 国产黄色小视频在线观看| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在 | 午夜福利成人在线免费观看| 日韩中文字幕欧美一区二区| 我要搜黄色片| av在线老鸭窝| 欧美日本亚洲视频在线播放| 午夜精品在线福利| 特级一级黄色大片| 国产高清激情床上av| 欧美色视频一区免费| 国产在线男女| 他把我摸到了高潮在线观看| 日本黄色片子视频| 亚洲av二区三区四区| 五月玫瑰六月丁香| 特大巨黑吊av在线直播| 亚洲中文字幕一区二区三区有码在线看| 麻豆国产97在线/欧美| 久久精品国产99精品国产亚洲性色| 三级男女做爰猛烈吃奶摸视频| 国产高清视频在线观看网站| 亚洲人成网站在线播| 成年女人看的毛片在线观看| 搡老妇女老女人老熟妇| 国产高清三级在线| 精品久久久久久,| 午夜福利在线观看吧| 欧美最黄视频在线播放免费| 久久精品夜夜夜夜夜久久蜜豆| 国产精品嫩草影院av在线观看 | 亚洲熟妇熟女久久| 国产久久久一区二区三区| 国产精品一区二区性色av| 免费看美女性在线毛片视频| 国产成人a区在线观看| 人妻夜夜爽99麻豆av| 88av欧美| 亚洲 国产 在线| 亚洲真实伦在线观看| 不卡视频在线观看欧美| 91午夜精品亚洲一区二区三区 | 日本撒尿小便嘘嘘汇集6| 亚洲av中文字字幕乱码综合| 悠悠久久av| 免费电影在线观看免费观看| 欧美黑人巨大hd| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产日韩欧美精品在线观看| 国产三级中文精品| 久久久久久大精品| 热99re8久久精品国产| 午夜视频国产福利| 国产免费一级a男人的天堂| 一进一出好大好爽视频| 美女高潮的动态| 全区人妻精品视频| 国产伦人伦偷精品视频| 成人亚洲精品av一区二区| 中文字幕高清在线视频| 成年女人毛片免费观看观看9| 日日摸夜夜添夜夜添av毛片 | a级毛片免费高清观看在线播放| 一个人免费在线观看电影| 三级毛片av免费| a级一级毛片免费在线观看| 全区人妻精品视频| 我的女老师完整版在线观看| 毛片一级片免费看久久久久 | 看十八女毛片水多多多| av国产免费在线观看| 最近中文字幕高清免费大全6 | 三级男女做爰猛烈吃奶摸视频| 精品久久久久久,| 我要搜黄色片| 91精品国产九色| 久久久久国产精品人妻aⅴ院| 午夜福利高清视频| 国产亚洲精品综合一区在线观看| 永久网站在线| 一级a爱片免费观看的视频| 久久人人爽人人爽人人片va| 99国产精品一区二区蜜桃av| 国产v大片淫在线免费观看| 国产探花极品一区二区| 99国产极品粉嫩在线观看| 我要搜黄色片| 午夜福利成人在线免费观看| 韩国av一区二区三区四区| 欧美成人免费av一区二区三区| 18+在线观看网站| 999久久久精品免费观看国产| 国产男人的电影天堂91| 99精品久久久久人妻精品| 91狼人影院| 草草在线视频免费看| av女优亚洲男人天堂| 久久久久免费精品人妻一区二区| 欧美又色又爽又黄视频| 毛片一级片免费看久久久久 | 日本熟妇午夜| 五月伊人婷婷丁香| 俺也久久电影网| 真实男女啪啪啪动态图| 亚洲 国产 在线| 99久久无色码亚洲精品果冻| 亚洲欧美日韩高清在线视频| 嫩草影视91久久| 美女xxoo啪啪120秒动态图| 亚洲人与动物交配视频| 精品国内亚洲2022精品成人| 久久精品91蜜桃| 在线免费观看的www视频| 国产精品嫩草影院av在线观看 | 国产午夜福利久久久久久| 亚洲av成人av| 99精品久久久久人妻精品| a级一级毛片免费在线观看| 神马国产精品三级电影在线观看| 最近视频中文字幕2019在线8| 啦啦啦啦在线视频资源| 日韩在线高清观看一区二区三区 | 色播亚洲综合网| av黄色大香蕉| 国产午夜福利久久久久久| 亚洲aⅴ乱码一区二区在线播放| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 亚洲男人的天堂狠狠| 久久欧美精品欧美久久欧美| 日韩 亚洲 欧美在线| 夜夜看夜夜爽夜夜摸| 亚洲av熟女| 亚洲成av人片在线播放无| 88av欧美| 又爽又黄a免费视频| 一级a爱片免费观看的视频| 欧美xxxx性猛交bbbb| 亚洲人成网站在线播放欧美日韩| 丝袜美腿在线中文| 又黄又爽又免费观看的视频| 天美传媒精品一区二区| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 国产色爽女视频免费观看| 午夜a级毛片| 啦啦啦啦在线视频资源| 国产三级在线视频| 中文字幕精品亚洲无线码一区| 国产熟女欧美一区二区| 夜夜爽天天搞| 国产黄片美女视频| 我要搜黄色片| 伦精品一区二区三区| 欧美一区二区国产精品久久精品| 成人美女网站在线观看视频| 日本精品一区二区三区蜜桃| 精品久久久久久久久久免费视频| 我的女老师完整版在线观看| 老师上课跳d突然被开到最大视频| 黄片wwwwww| 99热6这里只有精品| 永久网站在线| 久久久色成人| 日本-黄色视频高清免费观看| 69av精品久久久久久| 欧美激情国产日韩精品一区| 校园人妻丝袜中文字幕| 国产精品亚洲美女久久久| 国产一区二区激情短视频| av在线老鸭窝| 五月玫瑰六月丁香| 午夜精品一区二区三区免费看| 高清日韩中文字幕在线| 乱码一卡2卡4卡精品| 日本色播在线视频| 亚洲欧美精品综合久久99| 色在线成人网| 两人在一起打扑克的视频| 色av中文字幕| 久久久久免费精品人妻一区二区| 国产一区二区三区av在线 | 一个人看的www免费观看视频| 亚洲内射少妇av| 极品教师在线视频| 国产高清视频在线播放一区| av视频在线观看入口| 国产精品电影一区二区三区| 国产一区二区在线观看日韩| 简卡轻食公司| 黄色日韩在线| 综合色av麻豆| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 嫩草影院入口| 啦啦啦韩国在线观看视频| 无遮挡黄片免费观看| 国产精品久久久久久av不卡| 成人一区二区视频在线观看| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 中国美白少妇内射xxxbb| 国产精品国产三级国产av玫瑰| 国产69精品久久久久777片| 国产老妇女一区| 欧美性猛交黑人性爽| 日韩欧美在线二视频| 欧美又色又爽又黄视频| 中国美女看黄片| 中文字幕av在线有码专区| 在线观看一区二区三区| 能在线免费观看的黄片| 成人国产综合亚洲| 亚洲国产欧美人成| 99精品久久久久人妻精品| 成年女人永久免费观看视频| bbb黄色大片| 欧美日韩国产亚洲二区| 日韩精品有码人妻一区| 国产欧美日韩一区二区精品| 欧美日韩国产亚洲二区| 久久中文看片网| 国产亚洲欧美98| 日本-黄色视频高清免费观看| 黄色视频,在线免费观看| 老师上课跳d突然被开到最大视频| 亚洲真实伦在线观看| 男插女下体视频免费在线播放| 国产精品亚洲美女久久久| 亚洲18禁久久av| 中文字幕免费在线视频6| 日韩一本色道免费dvd| 久久午夜福利片| 久久久久免费精品人妻一区二区| 免费电影在线观看免费观看| 老司机午夜福利在线观看视频| 国产高清视频在线播放一区| 一区二区三区免费毛片| 91在线观看av| 亚洲av美国av| 女的被弄到高潮叫床怎么办 | 夜夜爽天天搞| 中文资源天堂在线| 五月玫瑰六月丁香| 欧美国产日韩亚洲一区| 国产一区二区三区av在线 | 免费av毛片视频| 亚洲一区高清亚洲精品| 少妇丰满av| 亚洲图色成人| 五月玫瑰六月丁香| 少妇的逼水好多| 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添av毛片 | 国产熟女欧美一区二区| 99在线视频只有这里精品首页| 两个人的视频大全免费| 久久精品国产亚洲网站| 免费看a级黄色片| 夜夜夜夜夜久久久久| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 国产高清视频在线观看网站| 可以在线观看毛片的网站| 丰满乱子伦码专区| 国产成人一区二区在线| 99久久无色码亚洲精品果冻| av.在线天堂| 亚洲国产色片| av在线观看视频网站免费| 国产单亲对白刺激| 国产免费av片在线观看野外av| av黄色大香蕉| 少妇的逼好多水| 午夜免费男女啪啪视频观看 | 欧美成人免费av一区二区三区| 国内揄拍国产精品人妻在线| 久久久久久久亚洲中文字幕| 国产男人的电影天堂91| a级一级毛片免费在线观看| 我要看日韩黄色一级片| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看吧| 日本-黄色视频高清免费观看| 国产精品一区www在线观看 | 成熟少妇高潮喷水视频| 美女黄网站色视频| 小蜜桃在线观看免费完整版高清| 日本爱情动作片www.在线观看 | 99久久中文字幕三级久久日本| 淫妇啪啪啪对白视频| 精品久久国产蜜桃| 伦精品一区二区三区| 欧美在线一区亚洲| 在线天堂最新版资源| 亚洲欧美日韩高清专用| 国产精品久久久久久亚洲av鲁大| 成熟少妇高潮喷水视频| 97超视频在线观看视频| 亚洲三级黄色毛片| 欧美一区二区亚洲| 国产高清不卡午夜福利| 精品一区二区三区视频在线| 午夜福利18| 亚洲不卡免费看| 丰满人妻一区二区三区视频av| 久久国产精品人妻蜜桃| 日韩强制内射视频| 看十八女毛片水多多多| 最近在线观看免费完整版| 亚洲av中文av极速乱 | 日韩欧美三级三区| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 国产精品1区2区在线观看.| 男人和女人高潮做爰伦理| 欧美一区二区精品小视频在线| 色视频www国产| 国产高清激情床上av| 久久久精品大字幕| 我要看日韩黄色一级片| 国语自产精品视频在线第100页| 99久久精品国产国产毛片| 国产爱豆传媒在线观看| 国产精品伦人一区二区| 免费av毛片视频| 日日夜夜操网爽| 日韩欧美免费精品| 久久香蕉精品热| 1024手机看黄色片| 亚洲国产精品sss在线观看| 久久精品国产鲁丝片午夜精品 | 国产精品乱码一区二三区的特点| 日韩高清综合在线| 日韩人妻高清精品专区| 蜜桃亚洲精品一区二区三区| 在线免费观看的www视频| 赤兔流量卡办理| 成人午夜高清在线视频| 一区二区三区激情视频| 一级黄片播放器| 午夜福利欧美成人| 精品一区二区三区av网在线观看| av在线天堂中文字幕| 国语自产精品视频在线第100页| 亚洲欧美日韩高清专用| 成人三级黄色视频| 老师上课跳d突然被开到最大视频| 国内久久婷婷六月综合欲色啪| 国产精品无大码| 亚洲一区高清亚洲精品| 2021天堂中文幕一二区在线观| 欧美色欧美亚洲另类二区| 亚洲精品在线观看二区| 欧美xxxx性猛交bbbb| 日韩一本色道免费dvd| 熟女人妻精品中文字幕| 男人舔奶头视频| 熟女人妻精品中文字幕| 国内少妇人妻偷人精品xxx网站| 亚洲国产色片| 成人av一区二区三区在线看| 真人一进一出gif抽搐免费| 国内精品美女久久久久久| 一本一本综合久久| 12—13女人毛片做爰片一| 午夜爱爱视频在线播放| 成人毛片a级毛片在线播放| 简卡轻食公司| 日韩一本色道免费dvd| aaaaa片日本免费| 日韩欧美精品v在线| 如何舔出高潮| 久久这里只有精品中国| 99久久成人亚洲精品观看| 国产精品,欧美在线| 亚洲人成网站高清观看| 少妇的逼好多水| 国产精品不卡视频一区二区| 男女之事视频高清在线观看| 成年版毛片免费区| 日韩欧美精品免费久久| av.在线天堂| 人人妻人人澡欧美一区二区| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 99riav亚洲国产免费| 99视频精品全部免费 在线| 国产亚洲av嫩草精品影院| 精品不卡国产一区二区三区| a在线观看视频网站| 亚洲第一区二区三区不卡| 在现免费观看毛片| 亚洲中文字幕一区二区三区有码在线看| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 久久精品国产亚洲网站| 村上凉子中文字幕在线| 婷婷丁香在线五月| 久久这里只有精品中国| 中文在线观看免费www的网站| 久久中文看片网| 免费观看精品视频网站| 色综合站精品国产| 国产亚洲精品av在线| 搡老妇女老女人老熟妇| h日本视频在线播放| 亚洲欧美清纯卡通| 免费在线观看影片大全网站| 51国产日韩欧美| 日本五十路高清| 一级av片app| 久久精品国产亚洲av香蕉五月| 久久久久久久亚洲中文字幕| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添av毛片 | 一本精品99久久精品77| av在线天堂中文字幕| 精品一区二区三区视频在线观看免费| 极品教师在线视频| 成年女人永久免费观看视频| 一级黄色大片毛片| 女人被狂操c到高潮| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 午夜激情福利司机影院| 精品乱码久久久久久99久播| 国产精品乱码一区二三区的特点| 能在线免费观看的黄片| 丰满的人妻完整版| 特级一级黄色大片| 午夜爱爱视频在线播放| 免费av毛片视频| 特级一级黄色大片| 乱系列少妇在线播放| 搡女人真爽免费视频火全软件 | 亚洲色图av天堂| 国内精品久久久久久久电影| 美女高潮喷水抽搐中文字幕| 日韩精品中文字幕看吧| 中文字幕精品亚洲无线码一区| 我的女老师完整版在线观看| 国产精品乱码一区二三区的特点| bbb黄色大片| 丰满人妻一区二区三区视频av| 少妇熟女aⅴ在线视频| 欧美xxxx性猛交bbbb| 亚洲一级一片aⅴ在线观看| 亚洲最大成人中文| 中文字幕精品亚洲无线码一区| 在线免费观看的www视频| 又黄又爽又免费观看的视频| 国内精品久久久久久久电影| 精品午夜福利在线看| 国语自产精品视频在线第100页| 亚洲成a人片在线一区二区| 天堂动漫精品| .国产精品久久| 成年免费大片在线观看| 久久精品影院6| 亚洲专区国产一区二区| 国产精品伦人一区二区| 久久精品综合一区二区三区| 国产精品亚洲一级av第二区| 特大巨黑吊av在线直播| or卡值多少钱| 亚洲第一区二区三区不卡| 免费黄网站久久成人精品| 听说在线观看完整版免费高清| 亚洲中文字幕一区二区三区有码在线看| 99热这里只有精品一区|