• <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
    12—13女人毛片做爰片一| 国产又色又爽无遮挡免费看| 国产精品自产拍在线观看55亚洲| 一级作爱视频免费观看| 久久国产精品人妻蜜桃| 日日夜夜操网爽| 欧美日韩av久久| 99国产极品粉嫩在线观看| 国产精品一区二区三区四区久久 | 国产精品久久久av美女十八| 丝袜美腿诱惑在线| 热99re8久久精品国产| 午夜免费鲁丝| 免费在线观看影片大全网站| 日韩欧美一区视频在线观看| 亚洲欧美一区二区三区黑人| 一级片'在线观看视频| 国产亚洲精品一区二区www| 999久久久国产精品视频| 午夜福利,免费看| 久久天躁狠狠躁夜夜2o2o| 97超级碰碰碰精品色视频在线观看| 一级毛片女人18水好多| 午夜视频精品福利| 免费高清视频大片| 精品熟女少妇八av免费久了| 高潮久久久久久久久久久不卡| 国产精品久久久久成人av| 国产深夜福利视频在线观看| 亚洲性夜色夜夜综合| 在线观看免费午夜福利视频| 脱女人内裤的视频| 久久人人97超碰香蕉20202| 亚洲久久久国产精品| 最近最新免费中文字幕在线| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 久久国产精品男人的天堂亚洲| 日韩视频一区二区在线观看| 国产成人影院久久av| 熟女少妇亚洲综合色aaa.| 亚洲自偷自拍图片 自拍| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 美女 人体艺术 gogo| xxxhd国产人妻xxx| 中文字幕精品免费在线观看视频| 国产深夜福利视频在线观看| 国产伦人伦偷精品视频| 人妻久久中文字幕网| 搡老乐熟女国产| 亚洲av成人一区二区三| 久久人人97超碰香蕉20202| av天堂在线播放| 琪琪午夜伦伦电影理论片6080| 在线观看午夜福利视频| 黄片播放在线免费| 欧美日韩黄片免| 亚洲,欧美精品.| 久久天躁狠狠躁夜夜2o2o| 性少妇av在线| 国产aⅴ精品一区二区三区波| 精品久久久精品久久久| 麻豆成人av在线观看| 亚洲精品中文字幕在线视频| 久久欧美精品欧美久久欧美| 女人高潮潮喷娇喘18禁视频| 男女下面进入的视频免费午夜 | 日韩大码丰满熟妇| 无人区码免费观看不卡| 成年女人毛片免费观看观看9| 国产精品久久电影中文字幕| 在线播放国产精品三级| 黄色a级毛片大全视频| 日韩免费高清中文字幕av| 香蕉久久夜色| 一边摸一边抽搐一进一出视频| 777久久人妻少妇嫩草av网站| ponron亚洲| 亚洲国产毛片av蜜桃av| 国产精品免费一区二区三区在线| 久久精品影院6| 国产精品久久久久久人妻精品电影| 日本vs欧美在线观看视频| 亚洲精品久久午夜乱码| 岛国在线观看网站| 久热爱精品视频在线9| 午夜精品国产一区二区电影| 亚洲精品一区av在线观看| 午夜免费成人在线视频| 男人舔女人下体高潮全视频| a级片在线免费高清观看视频| 久久久久久大精品| 高清黄色对白视频在线免费看| 大香蕉久久成人网| 精品无人区乱码1区二区| 狂野欧美激情性xxxx| 午夜久久久在线观看| 中文字幕av电影在线播放| 亚洲午夜精品一区,二区,三区| 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 91麻豆av在线| 不卡一级毛片| 在线天堂中文资源库| 欧美日韩福利视频一区二区| 99国产精品99久久久久| 亚洲精品成人av观看孕妇| 悠悠久久av| 久久国产乱子伦精品免费另类| avwww免费| 亚洲av电影在线进入| 这个男人来自地球电影免费观看| 国产成人欧美在线观看| 亚洲av片天天在线观看| 亚洲一区二区三区色噜噜 | 99精国产麻豆久久婷婷| 999久久久精品免费观看国产| 亚洲专区中文字幕在线| 亚洲国产看品久久| 免费高清视频大片| 精品福利永久在线观看| 免费在线观看亚洲国产| 欧美大码av| 九色亚洲精品在线播放| 久久狼人影院| 久久精品人人爽人人爽视色| 久久久久久久精品吃奶| 91国产中文字幕| 热99re8久久精品国产| 国产片内射在线| 久久国产亚洲av麻豆专区| 国产免费男女视频| 久久婷婷成人综合色麻豆| 80岁老熟妇乱子伦牲交| 亚洲av电影在线进入| 两性夫妻黄色片| 十分钟在线观看高清视频www| 9色porny在线观看| 亚洲av五月六月丁香网| 国产欧美日韩精品亚洲av| 国产国语露脸激情在线看| 又黄又粗又硬又大视频| 天堂动漫精品| 操美女的视频在线观看| 侵犯人妻中文字幕一二三四区| 亚洲,欧美精品.| 亚洲一区二区三区欧美精品| 97超级碰碰碰精品色视频在线观看| 亚洲成a人片在线一区二区| 久久久久久久久久久久大奶| 视频区图区小说| 夜夜躁狠狠躁天天躁| 免费日韩欧美在线观看| 丰满饥渴人妻一区二区三| 亚洲国产精品999在线| 久久青草综合色| 亚洲自偷自拍图片 自拍| 精品一区二区三区视频在线观看免费 | av超薄肉色丝袜交足视频| 日韩一卡2卡3卡4卡2021年| www.精华液| 好看av亚洲va欧美ⅴa在| 丝袜在线中文字幕| 男女午夜视频在线观看| 两个人免费观看高清视频| 久久香蕉国产精品| 狠狠狠狠99中文字幕| 18禁裸乳无遮挡免费网站照片 | 免费看a级黄色片| 久久久久久大精品| 三级毛片av免费| 免费在线观看日本一区| 757午夜福利合集在线观看| 精品一区二区三区av网在线观看| 夜夜看夜夜爽夜夜摸 | 最近最新免费中文字幕在线| 国产av一区在线观看免费| 水蜜桃什么品种好| 国产xxxxx性猛交| 亚洲七黄色美女视频| 成人亚洲精品av一区二区 | 午夜福利,免费看| 国产精品久久久久久人妻精品电影| 免费在线观看亚洲国产| 黄色毛片三级朝国网站| 久久香蕉精品热| 日本 av在线| 成在线人永久免费视频| 国产熟女午夜一区二区三区| 国产人伦9x9x在线观看| 精品福利观看| 乱人伦中国视频| 亚洲熟妇熟女久久| 最近最新免费中文字幕在线| 免费不卡黄色视频| 国产精品 欧美亚洲| netflix在线观看网站| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三| 高清黄色对白视频在线免费看| 性少妇av在线| 日韩大尺度精品在线看网址 | 夜夜看夜夜爽夜夜摸 | 成人黄色视频免费在线看| 国产成人av激情在线播放| 在线观看免费视频日本深夜| 国产av一区二区精品久久| 成人免费观看视频高清| 青草久久国产| 亚洲男人的天堂狠狠| 久久精品亚洲精品国产色婷小说| 亚洲精品在线观看二区| 午夜福利在线免费观看网站| 国产精华一区二区三区| 中文字幕色久视频| 99国产精品99久久久久| 国产一区二区三区视频了| 另类亚洲欧美激情| 满18在线观看网站| 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 九色亚洲精品在线播放| 久久国产精品男人的天堂亚洲| 日本a在线网址| 少妇的丰满在线观看| 色播在线永久视频| 波多野结衣高清无吗| 久久中文字幕一级| 搡老乐熟女国产| 国产精品一区二区精品视频观看| 国产欧美日韩一区二区三| 国产欧美日韩综合在线一区二区| 一区二区三区激情视频| 免费高清在线观看日韩| 成年女人毛片免费观看观看9| 美女国产高潮福利片在线看| 亚洲自拍偷在线| 国产一区在线观看成人免费| 午夜福利免费观看在线| 亚洲一卡2卡3卡4卡5卡精品中文| 精品第一国产精品| 亚洲色图av天堂| 日本 av在线| 中亚洲国语对白在线视频| 成人18禁高潮啪啪吃奶动态图| 日韩精品中文字幕看吧| 久久精品国产清高在天天线| 久久久久国内视频| 国产av一区在线观看免费| 黄色 视频免费看| 亚洲精品美女久久久久99蜜臀| 久久九九热精品免费| 国产高清国产精品国产三级| 亚洲精品国产区一区二| 天堂√8在线中文| 在线观看午夜福利视频| svipshipincom国产片| 久久久久国产精品人妻aⅴ院| 韩国精品一区二区三区| 亚洲精品粉嫩美女一区| 亚洲一卡2卡3卡4卡5卡精品中文| 日日爽夜夜爽网站| 国产深夜福利视频在线观看| av在线天堂中文字幕 | 国产精品久久电影中文字幕| 中文字幕人妻丝袜制服| 精品免费久久久久久久清纯| 中文字幕人妻丝袜制服| 婷婷丁香在线五月| 国产在线精品亚洲第一网站| 制服人妻中文乱码| 午夜两性在线视频| 精品乱码久久久久久99久播| 99在线视频只有这里精品首页| 亚洲精品av麻豆狂野| 变态另类成人亚洲欧美熟女 | 97超级碰碰碰精品色视频在线观看| 国产在线精品亚洲第一网站| 少妇粗大呻吟视频| 亚洲av成人av| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费 | 少妇粗大呻吟视频| 国产精品成人在线| 国产精品一区二区精品视频观看| 满18在线观看网站| 欧美性长视频在线观看| 在线永久观看黄色视频| 国产精品美女特级片免费视频播放器 | 中文字幕人妻丝袜制服| 日本三级黄在线观看| 大码成人一级视频| 黑人猛操日本美女一级片| www.熟女人妻精品国产| 日日爽夜夜爽网站| 欧美日韩精品网址| 男人舔女人的私密视频| 午夜福利欧美成人| 国产精品二区激情视频| 19禁男女啪啪无遮挡网站| 亚洲久久久国产精品| 黑丝袜美女国产一区| 精品福利永久在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产色婷婷电影| 久久精品亚洲精品国产色婷小说| 欧美 亚洲 国产 日韩一| 男人舔女人的私密视频| 久久中文看片网| 精品一区二区三区视频在线观看免费 | 香蕉久久夜色| 亚洲国产看品久久| 18禁裸乳无遮挡免费网站照片 | 婷婷六月久久综合丁香| 国产成年人精品一区二区 | 久久精品亚洲av国产电影网| 久久久久国产一级毛片高清牌| 国产视频一区二区在线看| 久久中文看片网| 18禁裸乳无遮挡免费网站照片 | 亚洲av成人不卡在线观看播放网| 91成年电影在线观看| 婷婷丁香在线五月| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 日韩免费高清中文字幕av| 成年人免费黄色播放视频| 精品久久久久久成人av| 日本黄色视频三级网站网址| 亚洲自偷自拍图片 自拍| 精品福利观看| av视频免费观看在线观看| 成人手机av| 欧美大码av| 人人妻人人澡人人看| 精品卡一卡二卡四卡免费| www.自偷自拍.com| 国产精品av久久久久免费| 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 国产高清激情床上av| 桃红色精品国产亚洲av| 日韩大尺度精品在线看网址 | 黄色成人免费大全| 长腿黑丝高跟| 国产区一区二久久| 在线观看免费视频日本深夜| av视频免费观看在线观看| 精品熟女少妇八av免费久了| 又黄又粗又硬又大视频| 80岁老熟妇乱子伦牲交| 夜夜爽天天搞| 亚洲欧洲精品一区二区精品久久久| 日韩精品中文字幕看吧| 91成人精品电影| 久久久久久大精品| 久久草成人影院| 99国产精品一区二区三区| 美女大奶头视频| 91av网站免费观看| 亚洲人成电影免费在线| 亚洲自偷自拍图片 自拍| 日韩欧美一区二区三区在线观看| 欧美日韩视频精品一区| 黑人操中国人逼视频| 亚洲国产看品久久| 久久久久国内视频| xxxhd国产人妻xxx| 国产精品综合久久久久久久免费 | 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 在线观看免费日韩欧美大片| 黄色毛片三级朝国网站| 午夜福利欧美成人| 99精国产麻豆久久婷婷| 91麻豆av在线| 欧美亚洲日本最大视频资源| 天堂√8在线中文| av片东京热男人的天堂| 少妇粗大呻吟视频| 成年人黄色毛片网站| 亚洲国产精品合色在线| 亚洲av成人av| 法律面前人人平等表现在哪些方面| 侵犯人妻中文字幕一二三四区| videosex国产| 夜夜夜夜夜久久久久| 老司机靠b影院| 国产精品av久久久久免费| 啦啦啦在线免费观看视频4| 久久狼人影院| 精品国产一区二区久久| 亚洲国产精品合色在线| 国产精品免费一区二区三区在线| 很黄的视频免费| 一级a爱视频在线免费观看| 视频区欧美日本亚洲| 美女福利国产在线| 亚洲九九香蕉| 99国产精品一区二区三区| 狠狠狠狠99中文字幕| 日韩精品中文字幕看吧| 视频在线观看一区二区三区| 一边摸一边做爽爽视频免费| 真人一进一出gif抽搐免费| 亚洲成人久久性| 久99久视频精品免费| 国产精品综合久久久久久久免费 | 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 精品一品国产午夜福利视频| av天堂在线播放| 国内毛片毛片毛片毛片毛片| 50天的宝宝边吃奶边哭怎么回事| 久久天堂一区二区三区四区| 午夜免费鲁丝| 91成年电影在线观看| 亚洲欧美一区二区三区久久| 精品免费久久久久久久清纯| 国产精品98久久久久久宅男小说| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 91国产中文字幕| 亚洲成人久久性| 99久久国产精品久久久| 欧美日本亚洲视频在线播放| 高清av免费在线| av欧美777| 日韩欧美在线二视频| 国产欧美日韩一区二区三| 人妻久久中文字幕网| 久久久久久免费高清国产稀缺| 亚洲一区二区三区色噜噜 | 日本免费a在线| 我的亚洲天堂| 日本免费一区二区三区高清不卡 | 日本a在线网址| 亚洲精品av麻豆狂野| 亚洲成人久久性| 99热只有精品国产| 99久久精品国产亚洲精品| 国产亚洲精品久久久久久毛片| 777久久人妻少妇嫩草av网站| 两个人免费观看高清视频| 精品久久久久久电影网| 亚洲中文字幕日韩| 一个人观看的视频www高清免费观看 | 悠悠久久av| 久久精品影院6| 精品欧美一区二区三区在线| 男女午夜视频在线观看| 亚洲一区二区三区欧美精品| 亚洲国产精品合色在线| 少妇裸体淫交视频免费看高清 | 精品高清国产在线一区| 啦啦啦免费观看视频1| 免费不卡黄色视频| 一进一出抽搐gif免费好疼 | 超碰成人久久| 国产精品影院久久| 女人高潮潮喷娇喘18禁视频| 18美女黄网站色大片免费观看| 级片在线观看| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 好看av亚洲va欧美ⅴa在| 成人亚洲精品一区在线观看| 亚洲精品久久午夜乱码| www.999成人在线观看| 高潮久久久久久久久久久不卡| 久久中文看片网| 亚洲精品国产色婷婷电影| 一区福利在线观看| 欧美成人免费av一区二区三区| 午夜精品国产一区二区电影| 99香蕉大伊视频| 一级a爱视频在线免费观看| 人成视频在线观看免费观看| 9热在线视频观看99| 夜夜爽天天搞| 日韩免费高清中文字幕av| 十八禁网站免费在线| xxx96com| 悠悠久久av| 国产精品偷伦视频观看了| 欧美精品啪啪一区二区三区| 五月开心婷婷网| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品999在线| 女警被强在线播放| 亚洲精品在线美女| 国产亚洲精品久久久久5区| 黄色怎么调成土黄色| 国产在线观看jvid| 午夜精品在线福利| 女同久久另类99精品国产91| 成在线人永久免费视频| 国产精品一区二区免费欧美| 一a级毛片在线观看| 国产成人一区二区三区免费视频网站| 一级a爱视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 色综合婷婷激情| 美女高潮喷水抽搐中文字幕| 成人亚洲精品av一区二区 | 中国美女看黄片| 黄色成人免费大全| 久久久国产成人精品二区 | 热99re8久久精品国产| 午夜精品国产一区二区电影| 日韩有码中文字幕| 18禁观看日本| 操出白浆在线播放| 国产精品免费视频内射| 天堂√8在线中文| 最新美女视频免费是黄的| 88av欧美| 99精品在免费线老司机午夜| 亚洲久久久国产精品| 中文字幕色久视频| 国产欧美日韩一区二区三区在线| 99香蕉大伊视频| 激情在线观看视频在线高清| 国产精品九九99| 两个人看的免费小视频| 日本 av在线| 亚洲国产中文字幕在线视频| 另类亚洲欧美激情| 美女扒开内裤让男人捅视频| 成人三级黄色视频| 首页视频小说图片口味搜索| 国产不卡一卡二| 法律面前人人平等表现在哪些方面| 国产又色又爽无遮挡免费看| 午夜久久久在线观看| 午夜a级毛片| 久久久国产一区二区| 国产野战对白在线观看| 欧美丝袜亚洲另类 | 在线天堂中文资源库| 色综合婷婷激情| 黄色视频,在线免费观看| 一级毛片精品| 精品久久久久久久毛片微露脸| 99久久99久久久精品蜜桃| 国产一区二区三区在线臀色熟女 | 五月开心婷婷网| 国产精品一区二区免费欧美| 80岁老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看| 美女大奶头视频| 成年女人毛片免费观看观看9| 在线观看一区二区三区激情| 国产97色在线日韩免费| 亚洲成人精品中文字幕电影 | 激情视频va一区二区三区| 亚洲免费av在线视频| 99热只有精品国产| 国产免费男女视频| 一区在线观看完整版| 美女午夜性视频免费| 亚洲 欧美一区二区三区| 51午夜福利影视在线观看| 丝袜在线中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲中文av在线| 手机成人av网站| 99香蕉大伊视频| 国产主播在线观看一区二区| 久久久久久久久中文| 国产黄色免费在线视频| 国产野战对白在线观看| 黄色a级毛片大全视频| 欧美成人性av电影在线观看| 老汉色∧v一级毛片| 成人亚洲精品av一区二区 | 超碰成人久久| 午夜精品久久久久久毛片777| 国产精品香港三级国产av潘金莲| 麻豆一二三区av精品| 18美女黄网站色大片免费观看| 午夜免费成人在线视频| 国产精品一区二区免费欧美| а√天堂www在线а√下载| 午夜视频精品福利| svipshipincom国产片| 黑丝袜美女国产一区| 亚洲 欧美 日韩 在线 免费| 首页视频小说图片口味搜索| 天堂√8在线中文| 国产乱人伦免费视频| aaaaa片日本免费| 人人妻人人爽人人添夜夜欢视频| 国产1区2区3区精品| 久久中文字幕人妻熟女| 亚洲av五月六月丁香网| 色精品久久人妻99蜜桃| 欧美日韩精品网址| 高清毛片免费观看视频网站 | 国产免费av片在线观看野外av| 夜夜躁狠狠躁天天躁| 国产亚洲欧美在线一区二区| 国产99久久九九免费精品| 女生性感内裤真人,穿戴方法视频| 男男h啪啪无遮挡| 国产一区二区三区在线臀色熟女 | 亚洲精品一二三| 色尼玛亚洲综合影院| 国产主播在线观看一区二区| 久久国产精品人妻蜜桃| 欧美日本中文国产一区发布|