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

    基于動態(tài)模態(tài)分解的長江口海表溫度時空分布特征重構(gòu)研究

    2022-02-21 04:37:12陳嘉星江迪張霄宇
    浙江大學學報(理學版) 2022年1期
    關(guān)鍵詞:觀測站長江口典型

    陳嘉星,江迪,張霄宇,3,4*

    (1.浙江大學地球科學學院,浙江 杭州 310027;2.浙江大學信息與電子工程學院,浙江 杭州 310027;3.浙江大學海南研究院,海南 三亞 572000;4.浙江大學海洋研究院,浙江 舟山 316000)

    海洋表面溫度(sea surface temperature,SST)是海洋-大氣系統(tǒng)物質(zhì)和能量交換過程的重要參數(shù)之一,是海水熱量傳遞的直接表現(xiàn)形式,SST的變化可影響海洋-大氣系統(tǒng)之間的熱量交換,進而影響區(qū)域氣候[1-2]。

    長江口海域陸地、海洋、大氣及人類活動等過程的相互作用強烈,水動力過程復雜,長江口海洋表面溫度的分布特征與變化趨勢影響并指示海水運動,是該海域海洋復雜內(nèi)在動力機制的綜合表現(xiàn)。此外,SST與表層葉綠素a濃度分布、海洋養(yǎng)殖、海洋生物的分布等密切相關(guān)[3]。研究長江口海洋表面溫度的時空演變規(guī)律對于理解時空耦合系統(tǒng)的深層動力機制、淺海環(huán)流及物質(zhì)輸運等過程極為重要。

    關(guān)于長江口海洋表面溫度的時空分布特征,前人已有不少研究。周曉英等[4]利用功率譜分析、MK突變檢驗分析了水文站的月均海溫數(shù)據(jù)和徑流量數(shù)據(jù),研究發(fā)現(xiàn)長江口海溫的最高值出現(xiàn)在8月,最低值出現(xiàn)在2月,季風和長江徑流對其有顯著影響。孫凡等[5]基于1982—2016年的OISST數(shù)據(jù)及其他氣象數(shù)據(jù),采用回歸分析和滯后相關(guān)系數(shù)法進行了研究,結(jié)果表明,長江口海域海表溫度的長期升溫趨勢主要受該區(qū)域徑向熱輸送及長江平流熱輸送增強的影響。買佳陽等[6]基于2000—2012年的中等分辨率成像光譜儀(moderate-resolution imaging spectroradiometer,MODIS-Terra)LIB數(shù) 據(jù),通 過優(yōu)化劈窗算法,建立了適合長江口水域的海表溫度反演模型,并進行了海表溫度反演,分析發(fā)現(xiàn),長江口海表溫度主要受太陽輻射影響,且口內(nèi)口外的海表溫度具有不同的季節(jié)性差異。何全軍等[7]基于2015—2017年我國周邊海域的衛(wèi)星數(shù)據(jù)和現(xiàn)場數(shù)據(jù),采用多元回歸模型,開發(fā)了適用于FY-3C/VIRR數(shù)據(jù)的區(qū)域SST反演算法。AI等[8]利用渤海地區(qū)的晴天紅外遙感數(shù)據(jù)和浮標實測數(shù)據(jù),構(gòu)建了一種適用于渤海地區(qū)的SST神經(jīng)網(wǎng)絡(luò)反演模型,并將反演結(jié)果與MODIS的SST產(chǎn)品進行了精度評估,結(jié)果表明,深度學習可應(yīng)用于SST反演。以上研究或偏向于水文資料分析,受水文站分布限制,分析僅局限于某幾個固定點位,無法擴展至整個海域;或注重海溫反演模型的建立,其反演模型是基于實測站位建立的,但研究選取的實測站位是否可代表整個海域有待商榷,同時也缺乏對長江口海溫時空演變模式的模態(tài)分析。

    20世紀以來,遙感技術(shù)迅速發(fā)展,借助遙感手段可獲得長時序、全覆蓋的SST數(shù)據(jù)。但因云層、霧、霾等影響,常常出現(xiàn)數(shù)據(jù)缺失問題,為獲得具有連續(xù)時間序列的全覆蓋衛(wèi)星遙感資料,國內(nèi)外已有多種數(shù)據(jù)插值方法,如克里金插值法、最優(yōu)插值法、期望最大化法、奇異譜方法、本征模態(tài)分解法、卡爾曼濾波法等。其中大多受限于數(shù)據(jù)的原始狀況,或依賴相關(guān)的參數(shù)和先驗值等信息,在適用性以及計算效率上均受限制。經(jīng)驗正交函數(shù)分解法(data interpolating empirical orthogonal functions,DINEOF)是目前較為有效的遙感數(shù)據(jù)重構(gòu)方法,具有無需先驗值、自適應(yīng)、時空相關(guān),以及適用于大面積缺測數(shù)據(jù)重構(gòu)等傳統(tǒng)重構(gòu)方法不具備的優(yōu)勢[9-10]。

    動態(tài)模態(tài)分解(dynamic mode decomposition,DMD)是Schmid于2008年提出的一種數(shù)據(jù)分析方法,其本質(zhì)是將高維非線性系統(tǒng)矩陣轉(zhuǎn)化為低維線性系統(tǒng)的過程,已被廣泛應(yīng)用于股票價格預測、圖像壓縮處理、物理流場分析等數(shù)據(jù)處理研究[11-13]。當數(shù)據(jù)的缺失率過高時,DINEOF具有很好的重構(gòu)效果,在能夠獲得足夠時序數(shù)據(jù)的情況下,數(shù)據(jù)驅(qū)動的算法能夠更好地解決動態(tài)系統(tǒng)的采樣問題。海洋是一個典型的高維動態(tài)系統(tǒng),尤其在近岸區(qū)域,受潮汐、風浪等因素影響,數(shù)據(jù)隨時間動態(tài)變化劇烈。因此,可將其視為由許多不同頻率的組分組成,DMD算法求得的模態(tài)和特征值分別描述了系統(tǒng)的空間特性和時間特性,符合相關(guān)海洋表面溫度研究的要求。DMD按照頻率對系統(tǒng)進行排序,提取系統(tǒng)特征頻率,觀察不同頻率的流場結(jié)構(gòu)對流場的貢獻,同時DMD模態(tài)特征值可預測流場[14]。在能夠獲取足夠時序數(shù)據(jù)的情況下,DMD算法可以更好地解決動態(tài)系統(tǒng)的采樣問題,從而實現(xiàn)對海洋表面溫度的預測。

    在實際應(yīng)用中,由于數(shù)據(jù)維度高、樣本數(shù)量大、數(shù)據(jù)本身存在噪聲等問題,DMD算法需針對具體應(yīng)用場景做相應(yīng)改進??紤]海洋學的研究特性,為節(jié)約成本、提高效率并為實際站位布設(shè)提供參考,本文擬結(jié)合正交三角(orthogonal right triangular,QR)分解算法[15]通過典型觀測站位還原SST時空分布情況,進行精度校驗,并與單一DMD算法進行比較。

    本文采用2003年1月至2019年12月MODIS L3級月平均SST遙感數(shù)據(jù),主要開展兩方面工作:(1)用DMD算法獲得SST基本模態(tài),并結(jié)合QR分解算法獲得典型觀測站位;(2)分別利用DMD算法、DMD結(jié)合QR分解算法對SST數(shù)據(jù)進行重構(gòu),并評價2種算法的重構(gòu)精度,對誤差原因進行初步分析。

    1 研究區(qū)域概況

    如圖1所示,研究區(qū)域位于118°E~128°E,26°N~34°N,為長江口及其鄰近海域,區(qū)域內(nèi)主要入海河流為長江和錢塘江等。長江口區(qū)域位于東亞季風區(qū),夏季盛行西南風,冬季盛行東北風。水文結(jié)構(gòu)一方面受近岸的長江沖淡水和沿岸流影響:北有低溫、低鹽、高沙的黃海沿岸流沿蘇北南下,西有低鹽、中沙的長江沖淡水輸入,其影響范圍存在季節(jié)性差異;另一方面受外海的臺灣暖流和黑潮次表層水的相互作用,兩者從臺灣島附近侵入,沿約200 m等深線北上,控制陸架外緣地區(qū),但存在明顯的季節(jié)性差異[16-17]。此外,長江沖淡水與臺灣暖流混合水團的鹽度梯度明顯,主要影響長江口外羽狀鋒區(qū)域[18-19]。除季風與環(huán)流相互作用外,長江口區(qū)域還受潮汐、風浪等因素影響,導致長江口區(qū)域的水動力環(huán)境較為復雜,進而影響SST的時空分布。

    圖1 研究區(qū)域地理位置示意[16-17]Fig.1 Geographical map of the study area[16-17]

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

    2.1 MODIS/SST

    MODIS中分辨率成像光譜儀是一種“圖譜合一”的光學儀器,可獲取從可見光到紅外共36個波段的觀測數(shù)據(jù),其最大空間分辨率為250 m。MODIS可提供包含陸地、海洋、云、大氣、溫度等特征信息的數(shù)據(jù)產(chǎn)品,廣泛用于陸地、海洋及大氣環(huán)境監(jiān)測等領(lǐng)域[20]。

    考慮數(shù)據(jù)的時間連續(xù)性、適中的空間分辨率以及獲取數(shù)據(jù)的便捷性,選用2003—2019年MODISAqua遙感器的L3級月平均SST產(chǎn)品,數(shù)據(jù)從https://oceandata.sci.gsfc.nasa.gov/獲取???間分辨率為4.6 km×4.6 km,每幀影像大小均為240×192,時間幀為204幀,各參數(shù)的觀測像元總數(shù)均為9 400 320個,且均采用NetCDF4格式。NetCDF4是一種科學數(shù)據(jù)存儲格式文件,能分層存儲數(shù)據(jù),且可用C、Fortran和MATLAB等語言進行操作。

    前處理主要包括掩膜制作、數(shù)據(jù)缺失率統(tǒng)計及空白數(shù)據(jù)填補、去噪處理等。

    首先,獲取全球海岸線shp文件(網(wǎng)址:https://www.ngdc.noaa.gov/mgg/shorelines/data/gshhs),導出研究區(qū)域的shp文件,并將導出的shp文件在ARCMAP中進行完善。然后,在MATLAB中將shp文件轉(zhuǎn)化為矩陣格式掩膜,令陸地部分為0,海洋部分為1,為方便處理,將陸地中的湖泊和河流賦值為0。掩膜處理完成后,統(tǒng)計海洋數(shù)據(jù),數(shù)據(jù)缺失率為1.57%。因數(shù)據(jù)缺失率較低,考慮程序的運行效率,用每幀的平均值填補空白數(shù)據(jù)。最后,對填補完成的數(shù)據(jù)進行去噪處理。

    2.2 DMD算法原理

    將海洋視為由不同頻率的組分組成的高維非線性系統(tǒng),由DMD算法得到的模態(tài)描述其空間特性。下面介紹DMD算法的原理。

    首先,從數(shù)值仿真或?qū)嶒炛惺占疦個時刻的快照,形成數(shù)據(jù)矩陣{ψ0,ψ1,…,ψN},第i個時刻的快照為第i列,假設(shè)上述數(shù)據(jù)矩陣在時間上是等距的,時間步長為Δt。利用快照序列構(gòu)建2個時間相鄰數(shù)據(jù)矩陣:

    利用連續(xù)快照之間的線性關(guān)系,通過矩陣A將式(1)中的2個快照矩陣相連接。矩陣A為高維系統(tǒng)的系統(tǒng)矩陣,包含大量動態(tài)關(guān)系,反映系統(tǒng)的動態(tài)特征。由于矩陣A的維數(shù)高,需要對其進行降階處理。動態(tài)模態(tài)分解提供了矩陣A的低階表示方法:

    DMD通過對上述快照矩陣進行數(shù)學變換,提取主導特征值及主要模態(tài)。本征正交分解(proper orthogonal decomposition,POD)又稱主成分分析,將高維系統(tǒng)表示為少量POD模與其對應(yīng)系數(shù)的線性組合,從而實現(xiàn)對高維系統(tǒng)降維。本文擬采用奇異值分解與POD相結(jié)合的方法,通過奇異值分解對高階算子進行相似變換,此方法具有較好的數(shù)值穩(wěn)定性。

    利用X的POD模,尋找秩為r的矩陣F,近似表示高維矩陣A:

    其中,U為通過對X進行奇異值分解所獲得的POD模,U*為酉矩陣,有

    其中,Σ為對角矩陣,對角線元素為r個奇異值,對應(yīng)特征值σ1,σ1,…,σr,同時奇異值分解得到的酉矩陣U和V滿 足:

    F的求解過程可看作將X和Y之間的Frobenius范數(shù)最小化,將式(3)和(4)代入,得到

    求解式(6),得到

    可用F近似表示A,記F的第i個特征值為μi,特征向量為ωi,則第i個DMD模態(tài)為

    由以上定義可知,模態(tài)是所研究的系統(tǒng)中具有某種全局特征的結(jié)構(gòu),本文研究的全局特征是SST的“頻率”,即SST的變化率。另外,若特征向量為實數(shù),則對應(yīng)的特征向量隨特征值的正負變化分別呈指數(shù)級增長或衰退。

    2.3 DMD結(jié)合QR分解算法確定典型觀測站位

    DMD算法僅能解決理想情況下流場基本模態(tài)的獲得及預測問題。在實際應(yīng)用中,由于數(shù)據(jù)維度高、樣本數(shù)量大、數(shù)據(jù)本身存在噪聲等,需根據(jù)具體應(yīng)用場景對DMD算法做相應(yīng)改進。QR分解算法用以求解矩陣特征值,可提高求解速度。用QR分解算法選擇典型觀測站位的思路是,對于給定的典型觀測站位數(shù)k,通過尋找最佳置換矩陣最大化|detAk|(det為矩陣的行列式值,Ak為QR分解算法中的上三角矩陣),尋找穩(wěn)定、準確的典型觀測站位[21]??紤]海洋學的研究特性,本文擬結(jié)合QR分解算法計算最佳模態(tài)數(shù)和典型觀測站位數(shù),利用典型觀測站位還原海溫數(shù)據(jù)。

    2.4 數(shù)據(jù)處理流程

    數(shù)據(jù)處理流程如下:

    (1)下載原始數(shù)據(jù),構(gòu)造矩陣。用制作好的掩膜,剔除每個時間幀原始數(shù)據(jù)中的陸地和島嶼,只保留海洋部分的數(shù)據(jù),計算平均值,用平均值替代其中的NaN(not a number),最后展開為列向量,有效數(shù)據(jù)共n個。對數(shù)據(jù)進行去噪處理,同時以p次間隔(p=204,時間間隔設(shè)為一個月)的觀測數(shù)據(jù)構(gòu)建矩陣Xp×n,時間設(shè)定為起始時間t1到結(jié)束時間tq,且滿足tq=tq-1+Δt。

    (2)抽取前80%的連續(xù)數(shù)據(jù)作為訓練數(shù)據(jù),數(shù)據(jù)量為m(m=163,時間間隔為一個月),其余的20%作為預測數(shù)據(jù),數(shù)據(jù)量為q(q=41,時間間隔為一個月)。

    (3)由上述訓練數(shù)據(jù)矩陣Xm×n構(gòu)造子矩陣。令子矩陣的時間間隔為Δt,得到Xm-11,Xm2:

    (4)求特征值和特征向量。經(jīng)步驟(3)后得到2個子矩陣,對X1m-1進行奇異值分解,Σ、U、V*分別為奇異值、左奇異向量、右奇異向量,同時結(jié)合矩陣計算低階矩陣F的特征值及其特征向量。

    (5)將得到的特征值和對應(yīng)的特征向量組成矩陣,根據(jù)QR分解算法,確定最佳模態(tài)數(shù)和最佳典型觀測站位數(shù),并將其代入FDMD中,重構(gòu)歷史數(shù)據(jù)。

    3 結(jié)果與討論

    3.1 長江口長時序SST的DMD模態(tài)分解

    利 用2003—2019年 的MODIS L3級 月 平 均SST產(chǎn)品,由DMD得到SST的主要時空演變模態(tài),用這一典型非線性動力學現(xiàn)象說明DMD在海洋系統(tǒng)中的應(yīng)用過程。訓練數(shù)據(jù)的時間段為2003年1月至2016年7月,預測數(shù)據(jù)的時間段為2016年8月至2019年12月。

    將DMD結(jié)合QR分解算法尋找典型觀測站位的方法(以下簡稱典型觀測站位還原法)應(yīng)用于上述訓練數(shù)據(jù),得到的重構(gòu)誤差與典型觀測站位數(shù)的關(guān)系如圖2所示。

    圖2 典型觀測站位效果對比Fig.2 Comparison of effect of typical observation stations

    其中,重構(gòu)誤差用均方根誤差(root mean square error,RMSE)表征,p表示典型觀測站位數(shù),r表示DMD模態(tài)數(shù)。由圖2可知,均方根誤差隨DMD模態(tài)數(shù)的增加而降低,且若DMD模態(tài)數(shù)不變,通過增加典型觀測站位數(shù)能有效降低重構(gòu)誤差。DMD算法得到的主要模態(tài)如圖3所示(因篇幅所限,此處只給出前4個模態(tài))。

    圖3 用DMD算法分解長江口海域SST數(shù)據(jù)得到的前4個主要模態(tài)Fig.3 The first four main modes decomposed by DMD algorithm of SST data in the Yangtze River Estuary

    由圖3可知,模態(tài)1在整個海域表現(xiàn)為負相關(guān),從空間分布看,蘇北沿岸、長江口及杭州灣和浙閩沿岸區(qū)域是SST變化率絕對值較大的區(qū)域,變小很快,且隨著長江沖淡水的外推呈舌狀遞增狀態(tài),而SST變化率絕對值較小值出現(xiàn)在沖繩島西北海域的大洋區(qū)域,變小很慢。推測出現(xiàn)這種分布的原因可能由近岸區(qū)域受陸地徑流輸入的淡水與海水混合所致,而低值區(qū)受黑潮和臺灣暖流影響,海水性質(zhì)比較均一。季節(jié)特征表現(xiàn)為由秋入冬。

    模態(tài)2中SST變化率與近岸區(qū)域及南側(cè)海域呈正相關(guān)、與中北側(cè)海域呈負相關(guān),表明2個區(qū)域呈現(xiàn)相反的變化趨勢。這可能反映地形對SST的影響,蘇北沿岸、長江口及杭州灣、浙閩沿岸區(qū)域水淺、比熱容小,易受大陸氣溫和沿岸流的影響,SST變化率最大值高于0.02℃,呈現(xiàn)SST變大、變快的特征;南側(cè)海域受黑潮和臺灣暖流的影響,表層SST變化不大;中北側(cè)由于遠離陸地,受環(huán)流影響小,SST變化率呈負值。季節(jié)特征表現(xiàn)為由春入夏。

    模態(tài)3中SST變化率呈正、負2種分布情況,蘇北沿岸、長江口及杭州灣、浙閩沿岸區(qū)域及北側(cè)海域為正值區(qū)域,南側(cè)海域為負值區(qū)域。這與第1模態(tài)所反映的受沖淡水、環(huán)流的影響相同。但季節(jié)特征表現(xiàn)為由春入夏。

    模態(tài)4則表現(xiàn)為近岸部分和150 m等深線以外的大洋區(qū)域為SST變化率正值區(qū),大致為被30 m和150 m等深線所夾區(qū)域為SST變化率負值區(qū)。這與第2模態(tài)反映的地形影響一致。但季節(jié)特征表現(xiàn)為由秋入冬。

    3.2 DMD結(jié)合QR分解算法獲取典型觀測站位

    上述結(jié)果表明,當DMD模態(tài)為25、典型觀測站位為25時,重構(gòu)誤差已經(jīng)相當小。因此,考慮運行速度和重構(gòu)精度,本文選取25個DMD模態(tài)、25個典型觀測站位對長江口海域的海表溫度進行稀疏還原,結(jié)果如圖4所示。

    圖4 2016年8月長江口海域SST數(shù)據(jù)稀疏還原結(jié)果Fig.4 The results of sparse restoration of SST data in Yangtze River Estuary in August 2016

    由圖4可知,相比單純采用DMD算法,DMD結(jié)合QR分解算法的還原結(jié)果能很好地代表整個海域,海溫分布與真值更接近。因此,若不考慮成本和運行速度,單一追求重構(gòu)精度,用DMD算法能得到更高精度的還原結(jié)果。若綜合考慮成本、運行速度、重構(gòu)精度,則DMD結(jié)合QR分解算法的重構(gòu)表現(xiàn)更佳。

    3.3 精度分析

    選取2016年8月相同區(qū)域的原始數(shù)據(jù)與重構(gòu)數(shù)據(jù)的像元值進行精度評價,評價結(jié)果如圖5所示。從圖5(a)和(b)的對比中可知,2種方法的原始值與重構(gòu)值集中分布在直線y=x附近,精度指標均較理想。其中,DMD結(jié)合QR分解算法的RMSE為0.005 0。據(jù)統(tǒng)計結(jié)果,測試數(shù)據(jù)集的歷史數(shù)據(jù)與DMD結(jié)合QR分解算法重構(gòu)數(shù)據(jù)之間的RMSE在0.005 0~0.011 2,平均值為0.007 6。因此,雖然在精度上DMD算法的最佳還原結(jié)果比DMD結(jié)合QR分解算法的還原結(jié)果表現(xiàn)更佳,但考慮運行速度和成本,仍可認為DMD結(jié)合QR分解算法的結(jié)果較好。

    另外,由圖5可知,無論是DMD算法還是DMD結(jié)合QR分解算法重構(gòu)得到的結(jié)果均具有相當高的精度。為進一步分析重構(gòu)效果,隨機選取2個典型觀測站位,對其原始值、DMD最佳還原值、典型觀測站位還原值做對比分析。

    由圖6可知,DMD最佳還原值和典型觀測站位還原值與原始值耦合良好,說明DMD算法和DMD結(jié)合QR分解算法均較合理可靠。因此,可以認為選取的25個典型觀測站位能較好地代表整個研究海域海洋表面溫度的分布特征。這與圖5的DMD算法和DMD結(jié)合QR分解算法的結(jié)果高度吻合。

    為更好地評價2種方法的重構(gòu)效果,選擇預測數(shù) 據(jù) 集 中 的2017年4月、2017年7月、2017年10月及2018年1月(分別代表春、夏、秋、冬4個季節(jié))進行重構(gòu),得到的結(jié)果如圖7所示。

    由圖7可知,從空間分布上看,2種方法的重構(gòu)結(jié)果均與原始值的趨勢保持一致,誤差較小,這與圖5和圖6得到的結(jié)論一致。

    圖5 2016年8月的SST重構(gòu)精度Fig.5 The accuracy of SST reconstruction in August 2016

    圖6 某2個典型觀測站位重構(gòu)前后對比Fig.6 Comparison of reconstruction of two typical observation stations

    圖7 長江口海域典型年份四季SST數(shù)據(jù)稀疏還原結(jié)果Fig.7 Sparse restoration results of SST data in four seasons of a particular year in the Yangtze River Estuary

    4 結(jié)論

    基于DMD算法的基本理論,將QR分解算法應(yīng)用于海洋學觀測,分析了長江口海域的海溫數(shù)據(jù)。結(jié)果表明,DMD結(jié)合QR分解算法對長江口海溫數(shù)據(jù)進行稀疏還原,能有效去除噪聲,填補缺失數(shù)據(jù)。增加DMD模態(tài)數(shù)能有效降低重構(gòu)誤差,當模態(tài)數(shù)一定時,可通過增加典型觀測站位數(shù)降低重構(gòu)誤差。當采用25個模態(tài)、25個典型觀測站位時,重構(gòu)數(shù)據(jù)的RMSE為0.005 0。

    DMD算法分解得到的前4個模態(tài)中,模態(tài)1和模態(tài)3可能反映了長江沖淡水及環(huán)流的影響,模態(tài)1的季節(jié)特征為由秋入冬,模態(tài)3則為由春入夏;模態(tài)2和模態(tài)4可能代表海底地形的影響,模態(tài)2季節(jié)特征為由春入夏,模態(tài)4則為由秋入冬。

    在數(shù)據(jù)充足的情況下,DMD算法能很好地解決動態(tài)系統(tǒng)的采樣問題,對時序數(shù)據(jù)的重構(gòu)和分析具有應(yīng)用價值。單純采用DMD算法重構(gòu),精度相對較高,但采用DMD結(jié)合QR分解算法重構(gòu)效率較高,并且可在稀疏站位條件下重構(gòu)海洋溫度場,這對海洋學調(diào)查尤為重要。

    猜你喜歡
    觀測站長江口典型
    GPS導航對抗數(shù)據(jù)質(zhì)量特征實例分析
    用最典型的事寫最有特點的人
    四川省甘孜州:航拍四川稻城高海拔宇宙線觀測站
    科學導報(2022年39期)2022-07-04 23:15:46
    多項式求值題的典型解法
    典型胰島素瘤1例報道
    去中心化時差頻差直接定位方法
    航空學報(2017年5期)2017-11-20 03:44:46
    長江口鳳鱭繁殖群體的年齡結(jié)構(gòu)和生長特性
    長江口及其鄰近水域仔稚魚種類組成及分布特征
    風場對夏季長江口外東北海域低鹽水團的作用
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    人妻系列 视频| 国产毛片a区久久久久| 国产精品综合久久久久久久免费| 日日撸夜夜添| 久久国内精品自在自线图片| av线在线观看网站| 色5月婷婷丁香| 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| 免费观看性生交大片5| 韩国av在线不卡| 久久这里有精品视频免费| 欧美日韩一区二区视频在线观看视频在线 | 少妇熟女欧美另类| 男人和女人高潮做爰伦理| 精品人妻视频免费看| 亚洲电影在线观看av| 国产极品精品免费视频能看的| 久久精品91蜜桃| 久久久久精品久久久久真实原创| 精品免费久久久久久久清纯| 中文字幕制服av| 天堂av国产一区二区熟女人妻| av在线蜜桃| 亚洲av成人精品一二三区| 天美传媒精品一区二区| 免费看日本二区| 91精品一卡2卡3卡4卡| 中文字幕av成人在线电影| 日韩欧美三级三区| 亚洲国产色片| 成人午夜精彩视频在线观看| 2022亚洲国产成人精品| 亚洲av一区综合| 日韩一区二区三区影片| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 毛片一级片免费看久久久久| 午夜久久久久精精品| 黄色日韩在线| 亚州av有码| 午夜老司机福利剧场| 日本三级黄在线观看| 国产免费男女视频| 亚洲欧洲国产日韩| 亚洲精品影视一区二区三区av| 亚洲国产成人一精品久久久| 色综合色国产| av免费在线看不卡| 国产片特级美女逼逼视频| 午夜精品一区二区三区免费看| 日本黄色视频三级网站网址| 久久精品影院6| 精品一区二区三区视频在线| kizo精华| 成人毛片60女人毛片免费| 亚洲国产欧美在线一区| 久久久国产成人免费| 啦啦啦观看免费观看视频高清| 国产免费一级a男人的天堂| 亚洲自偷自拍三级| 热99re8久久精品国产| 成人午夜精彩视频在线观看| 日韩欧美精品v在线| 99九九线精品视频在线观看视频| 特大巨黑吊av在线直播| 插逼视频在线观看| 亚洲国产欧洲综合997久久,| 亚洲国产精品成人久久小说| 亚洲无线观看免费| 亚洲欧洲国产日韩| 日韩高清综合在线| 久久这里有精品视频免费| 日韩av不卡免费在线播放| 欧美一区二区国产精品久久精品| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 少妇被粗大猛烈的视频| 国产亚洲精品久久久com| 黄色欧美视频在线观看| av播播在线观看一区| 日本av手机在线免费观看| 中文字幕亚洲精品专区| av.在线天堂| 国产亚洲91精品色在线| 99热这里只有是精品在线观看| 国产成人一区二区在线| 午夜a级毛片| 亚洲高清免费不卡视频| 亚洲在线观看片| 国产午夜福利久久久久久| 亚洲国产欧美在线一区| 桃色一区二区三区在线观看| 波野结衣二区三区在线| 又粗又硬又长又爽又黄的视频| 久99久视频精品免费| 亚洲精品成人久久久久久| 亚洲精华国产精华液的使用体验| 久久精品夜色国产| 天堂影院成人在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品日韩av片在线观看| 超碰av人人做人人爽久久| 热99re8久久精品国产| 可以在线观看毛片的网站| 国产熟女欧美一区二区| 美女黄网站色视频| 2021少妇久久久久久久久久久| 男女那种视频在线观看| 国产精品久久久久久久电影| 免费av毛片视频| 一本久久精品| 天堂网av新在线| 我的老师免费观看完整版| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 春色校园在线视频观看| 国产淫语在线视频| 村上凉子中文字幕在线| 欧美成人a在线观看| 九九爱精品视频在线观看| 国产精品一区二区三区四区免费观看| 国产一区亚洲一区在线观看| 狂野欧美白嫩少妇大欣赏| 久久精品熟女亚洲av麻豆精品 | 日韩强制内射视频| 免费在线观看成人毛片| www.av在线官网国产| 在线观看一区二区三区| 男人和女人高潮做爰伦理| 联通29元200g的流量卡| 亚洲自拍偷在线| 成人特级av手机在线观看| 秋霞在线观看毛片| 国产精品.久久久| 久久婷婷人人爽人人干人人爱| 少妇猛男粗大的猛烈进出视频 | 干丝袜人妻中文字幕| 亚洲在久久综合| 久久精品国产自在天天线| 国产精品久久视频播放| 欧美潮喷喷水| 哪个播放器可以免费观看大片| 日韩av在线大香蕉| 久久久久久伊人网av| av国产久精品久网站免费入址| 最近最新中文字幕大全电影3| 国产免费福利视频在线观看| 中文字幕熟女人妻在线| 亚洲精品乱码久久久v下载方式| 国产老妇女一区| 国产高清视频在线观看网站| 免费看av在线观看网站| 国产熟女欧美一区二区| 一级二级三级毛片免费看| 99热这里只有是精品在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲欧洲国产日韩| 精品人妻熟女av久视频| 岛国毛片在线播放| 高清av免费在线| 99热全是精品| 国产高清不卡午夜福利| 九九热线精品视视频播放| 一级av片app| 亚洲成人av在线免费| 精品午夜福利在线看| av女优亚洲男人天堂| 国产亚洲精品av在线| 哪个播放器可以免费观看大片| 亚洲国产最新在线播放| 成人一区二区视频在线观看| av国产免费在线观看| 久久午夜福利片| 国产色婷婷99| www日本黄色视频网| 一级毛片我不卡| 成人一区二区视频在线观看| 午夜a级毛片| 日本一二三区视频观看| 久久久久免费精品人妻一区二区| 国产成人91sexporn| 久99久视频精品免费| 一个人看视频在线观看www免费| 99久久精品热视频| 建设人人有责人人尽责人人享有的 | 国产成人午夜福利电影在线观看| 日韩制服骚丝袜av| 在线免费观看不下载黄p国产| 内射极品少妇av片p| 精品久久久久久电影网 | av免费观看日本| 日韩欧美精品免费久久| 亚洲丝袜综合中文字幕| 亚洲av成人精品一二三区| 久久久久网色| 亚洲经典国产精华液单| 91久久精品国产一区二区三区| 91精品国产九色| 人体艺术视频欧美日本| 嫩草影院新地址| 亚洲精品乱码久久久v下载方式| 日韩欧美 国产精品| 一个人免费在线观看电影| 高清在线视频一区二区三区 | 99久久精品国产国产毛片| 国产亚洲精品av在线| 色综合亚洲欧美另类图片| 亚洲欧美日韩高清专用| 国语自产精品视频在线第100页| 高清av免费在线| 欧美成人精品欧美一级黄| 国产精品久久久久久久电影| av视频在线观看入口| 成人综合一区亚洲| 最新中文字幕久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲色图av天堂| 国产黄a三级三级三级人| 特大巨黑吊av在线直播| 2022亚洲国产成人精品| 99久久精品国产国产毛片| 亚洲,欧美,日韩| 国产精品一及| 国产一级毛片在线| 精品人妻偷拍中文字幕| 一级毛片我不卡| 久久久久免费精品人妻一区二区| 少妇的逼好多水| 嘟嘟电影网在线观看| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图| 国产av一区在线观看免费| 日产精品乱码卡一卡2卡三| 欧美成人a在线观看| 国产午夜精品一二区理论片| 在线免费观看的www视频| a级毛色黄片| 一级黄色大片毛片| 一个人免费在线观看电影| 欧美一级a爱片免费观看看| 婷婷色综合大香蕉| 少妇人妻精品综合一区二区| 老司机影院毛片| 综合色av麻豆| 欧美97在线视频| 亚洲成人中文字幕在线播放| 最近的中文字幕免费完整| 国产精品1区2区在线观看.| 久久久久久久午夜电影| 十八禁国产超污无遮挡网站| 在线观看一区二区三区| 免费看av在线观看网站| 久久欧美精品欧美久久欧美| 亚洲av免费高清在线观看| 麻豆成人av视频| 在线免费观看的www视频| 国产v大片淫在线免费观看| 日韩欧美 国产精品| 国模一区二区三区四区视频| www.av在线官网国产| 国产视频首页在线观看| 大又大粗又爽又黄少妇毛片口| 男人舔女人下体高潮全视频| 久久草成人影院| 日韩欧美 国产精品| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区| 91久久精品国产一区二区三区| 国产高清三级在线| 色吧在线观看| 如何舔出高潮| 亚洲欧美精品专区久久| 午夜激情欧美在线| 免费观看在线日韩| av国产久精品久网站免费入址| 特级一级黄色大片| 长腿黑丝高跟| 丝袜喷水一区| 搡女人真爽免费视频火全软件| a级毛色黄片| 91久久精品电影网| 午夜精品在线福利| 中文亚洲av片在线观看爽| av在线播放精品| 国产精品不卡视频一区二区| 国产av不卡久久| 别揉我奶头 嗯啊视频| 欧美日本亚洲视频在线播放| 在线免费观看不下载黄p国产| 国产一级毛片七仙女欲春2| 亚洲五月天丁香| 搞女人的毛片| 听说在线观看完整版免费高清| 天堂网av新在线| 亚洲国产精品专区欧美| 国产单亲对白刺激| 在线观看av片永久免费下载| 久久久久网色| 成人毛片60女人毛片免费| 18禁在线无遮挡免费观看视频| 国产伦一二天堂av在线观看| 国产高清国产精品国产三级 | 国产成人福利小说| 大又大粗又爽又黄少妇毛片口| 亚洲综合色惰| 日本三级黄在线观看| 91aial.com中文字幕在线观看| 人妻少妇偷人精品九色| 人妻夜夜爽99麻豆av| 久久久久九九精品影院| 精品国产一区二区三区久久久樱花 | 午夜激情福利司机影院| 一区二区三区高清视频在线| 亚洲国产最新在线播放| 99视频精品全部免费 在线| 色播亚洲综合网| 伦精品一区二区三区| 国产精品爽爽va在线观看网站| 国产免费又黄又爽又色| 天堂√8在线中文| videos熟女内射| 26uuu在线亚洲综合色| 色视频www国产| 日韩一本色道免费dvd| 久久国内精品自在自线图片| 夜夜看夜夜爽夜夜摸| 日日干狠狠操夜夜爽| 久久久久免费精品人妻一区二区| 偷拍熟女少妇极品色| 午夜视频国产福利| 国国产精品蜜臀av免费| 99久久精品一区二区三区| 国产一区亚洲一区在线观看| 日韩视频在线欧美| 久久人人爽人人爽人人片va| 午夜精品一区二区三区免费看| 丰满人妻一区二区三区视频av| 高清毛片免费看| av国产免费在线观看| 一级av片app| 亚洲自拍偷在线| 麻豆av噜噜一区二区三区| 亚洲激情五月婷婷啪啪| 人妻系列 视频| 国内少妇人妻偷人精品xxx网站| 亚洲最大成人手机在线| 国产免费福利视频在线观看| 麻豆成人午夜福利视频| 久久久亚洲精品成人影院| 国产精品嫩草影院av在线观看| 欧美区成人在线视频| 欧美成人a在线观看| 看十八女毛片水多多多| 嫩草影院入口| av天堂中文字幕网| 男人的好看免费观看在线视频| av福利片在线观看| 欧美成人一区二区免费高清观看| 三级男女做爰猛烈吃奶摸视频| 毛片一级片免费看久久久久| .国产精品久久| 免费搜索国产男女视频| 黄色配什么色好看| 色播亚洲综合网| 亚洲色图av天堂| 一个人观看的视频www高清免费观看| 不卡视频在线观看欧美| av视频在线观看入口| 黄色配什么色好看| 久久午夜福利片| 亚洲成人av在线免费| 日韩欧美精品免费久久| 国内精品美女久久久久久| 国产精品熟女久久久久浪| 日本黄大片高清| 日韩欧美精品免费久久| 亚洲av日韩在线播放| 久久久久网色| 婷婷色麻豆天堂久久 | 男人舔女人下体高潮全视频| 直男gayav资源| 亚洲美女视频黄频| 久久精品国产鲁丝片午夜精品| 搡女人真爽免费视频火全软件| 色尼玛亚洲综合影院| 亚洲av.av天堂| 天堂影院成人在线观看| 亚洲av成人av| 99热精品在线国产| 成人特级av手机在线观看| 久久久久久久久中文| 男女国产视频网站| 国产又色又爽无遮挡免| 久久综合国产亚洲精品| 伦理电影大哥的女人| 久久6这里有精品| 国产毛片a区久久久久| 欧美zozozo另类| 色综合色国产| 成人无遮挡网站| 欧美区成人在线视频| 久久国内精品自在自线图片| 日韩人妻高清精品专区| 人妻系列 视频| 少妇人妻一区二区三区视频| 欧美成人免费av一区二区三区| 精品国产三级普通话版| 天堂影院成人在线观看| АⅤ资源中文在线天堂| 欧美变态另类bdsm刘玥| 一级二级三级毛片免费看| 亚洲欧美日韩东京热| 性插视频无遮挡在线免费观看| 最近中文字幕2019免费版| 亚洲综合精品二区| 91久久精品电影网| 亚洲一区高清亚洲精品| 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 高清日韩中文字幕在线| 国产欧美日韩精品一区二区| or卡值多少钱| 黄色一级大片看看| 能在线免费看毛片的网站| 大话2 男鬼变身卡| 麻豆乱淫一区二区| 国产 一区精品| 一卡2卡三卡四卡精品乱码亚洲| 久久精品夜夜夜夜夜久久蜜豆| 亚洲图色成人| 国产精品人妻久久久影院| 国产精品麻豆人妻色哟哟久久 | 水蜜桃什么品种好| 亚洲自拍偷在线| 蜜桃久久精品国产亚洲av| 色综合亚洲欧美另类图片| 久久久精品大字幕| 丰满乱子伦码专区| 内射极品少妇av片p| 欧美成人a在线观看| 精品熟女少妇av免费看| 国产美女午夜福利| 精品久久久久久久久亚洲| 国产乱人偷精品视频| 国产亚洲91精品色在线| 亚洲性久久影院| 色播亚洲综合网| 国产熟女欧美一区二区| 三级毛片av免费| 国产三级在线视频| 亚洲av不卡在线观看| 日韩制服骚丝袜av| 亚洲丝袜综合中文字幕| 国产一区二区在线av高清观看| 国产大屁股一区二区在线视频| 婷婷色综合大香蕉| 51国产日韩欧美| 国产高潮美女av| 免费搜索国产男女视频| 中国美白少妇内射xxxbb| 97超视频在线观看视频| 久热久热在线精品观看| 午夜福利高清视频| 国模一区二区三区四区视频| 看免费成人av毛片| 91久久精品国产一区二区三区| 天堂网av新在线| 搡女人真爽免费视频火全软件| 男人狂女人下面高潮的视频| 欧美一区二区精品小视频在线| 九九爱精品视频在线观看| 亚洲精品乱码久久久v下载方式| 欧美又色又爽又黄视频| 99热全是精品| 午夜日本视频在线| 精品人妻熟女av久视频| 国产成人a区在线观看| 精品一区二区免费观看| 免费看av在线观看网站| 91aial.com中文字幕在线观看| 桃色一区二区三区在线观看| 男人舔奶头视频| 51国产日韩欧美| 91午夜精品亚洲一区二区三区| 国产精品久久久久久精品电影小说 | 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| 男女那种视频在线观看| 高清日韩中文字幕在线| 99久久成人亚洲精品观看| 久久久成人免费电影| 中文字幕免费在线视频6| 男人舔女人下体高潮全视频| 女的被弄到高潮叫床怎么办| 亚洲自拍偷在线| 国产v大片淫在线免费观看| 国产一区二区三区av在线| 欧美丝袜亚洲另类| 国产免费一级a男人的天堂| 99久国产av精品国产电影| 日本-黄色视频高清免费观看| 国语对白做爰xxxⅹ性视频网站| 老司机影院毛片| 久久韩国三级中文字幕| 亚洲经典国产精华液单| 色播亚洲综合网| 在线天堂最新版资源| 嫩草影院入口| 精品不卡国产一区二区三区| 亚洲国产精品成人综合色| 免费看光身美女| 亚洲av熟女| 白带黄色成豆腐渣| 少妇裸体淫交视频免费看高清| 久久久久久久久久成人| 日韩成人av中文字幕在线观看| 日韩国内少妇激情av| 国产高清有码在线观看视频| 少妇熟女aⅴ在线视频| 69人妻影院| 黄色配什么色好看| 岛国在线免费视频观看| 久久鲁丝午夜福利片| 69av精品久久久久久| 日韩,欧美,国产一区二区三区 | 麻豆成人av视频| 欧美一区二区亚洲| av在线亚洲专区| 超碰av人人做人人爽久久| 久久久久久久久久成人| 两个人视频免费观看高清| 网址你懂的国产日韩在线| 少妇猛男粗大的猛烈进出视频 | 一级黄片播放器| 日韩av在线免费看完整版不卡| 亚洲国产色片| 22中文网久久字幕| 亚洲成人久久爱视频| 免费看a级黄色片| 又爽又黄无遮挡网站| 哪个播放器可以免费观看大片| 午夜精品国产一区二区电影 | 成人二区视频| 看非洲黑人一级黄片| 成人亚洲精品av一区二区| 欧美bdsm另类| 极品教师在线视频| 国产不卡一卡二| 天堂影院成人在线观看| 爱豆传媒免费全集在线观看| 久久久久久久亚洲中文字幕| 国产成年人精品一区二区| 久久精品熟女亚洲av麻豆精品 | 亚洲欧洲国产日韩| 精品久久久久久久久久久久久| 我的女老师完整版在线观看| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 亚洲性久久影院| 毛片一级片免费看久久久久| 老司机影院毛片| 久久久色成人| 亚洲精品久久久久久婷婷小说 | 国产精品1区2区在线观看.| 一级毛片我不卡| 边亲边吃奶的免费视频| 三级毛片av免费| 亚洲欧美日韩无卡精品| 亚洲国产成人一精品久久久| 亚洲综合色惰| 国产色婷婷99| 噜噜噜噜噜久久久久久91| 欧美人与善性xxx| 欧美成人a在线观看| 免费观看在线日韩| 久久精品夜色国产| 深夜a级毛片| 国产av一区在线观看免费| 亚洲18禁久久av| 免费一级毛片在线播放高清视频| 久久综合国产亚洲精品| 99久久九九国产精品国产免费| 嫩草影院入口| 亚洲第一区二区三区不卡| 午夜激情欧美在线| 国产真实伦视频高清在线观看| 国产久久久一区二区三区| 国产精品,欧美在线| 大又大粗又爽又黄少妇毛片口| 激情 狠狠 欧美| 欧美激情久久久久久爽电影| 六月丁香七月| 两个人的视频大全免费| 国产精品一二三区在线看| 婷婷六月久久综合丁香| 夜夜爽夜夜爽视频| 日韩欧美三级三区| 日韩欧美国产在线观看| 国产高清视频在线观看网站| 欧美日本视频| 青春草亚洲视频在线观看| 亚洲精华国产精华液的使用体验| 欧美性猛交╳xxx乱大交人| 亚洲av二区三区四区| 老司机影院成人| av播播在线观看一区| 99久久人妻综合| 国产一区二区在线av高清观看| 国产免费男女视频| 久久午夜福利片| 欧美一区二区精品小视频在线| 亚洲精品456在线播放app|