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

    國(guó)際全球海洋環(huán)流預(yù)報(bào)系統(tǒng)的現(xiàn)狀與展望

    2020-07-29 09:41:42韓鵬李宇航揭曉蒙
    海洋預(yù)報(bào) 2020年3期
    關(guān)鍵詞:環(huán)流動(dòng)力數(shù)值

    韓鵬,李宇航,揭曉蒙

    (中國(guó)21世紀(jì)議程管理中心,北京100038)

    1 引言

    海洋動(dòng)力環(huán)境預(yù)報(bào)是基于對(duì)海洋過(guò)去和當(dāng)前狀態(tài)及其演變規(guī)律的認(rèn)識(shí),運(yùn)用數(shù)值模式、資料同化等手段,對(duì)多種時(shí)空尺度上的海洋現(xiàn)象和海洋狀況(海溫、鹽度、海流、海浪、潮汐等)進(jìn)行分析和預(yù)測(cè),是業(yè)務(wù)化海洋學(xué)研究中的關(guān)鍵技術(shù)之一[1-2]。海洋動(dòng)力環(huán)境預(yù)報(bào)技術(shù)主要包括動(dòng)力環(huán)境數(shù)值預(yù)報(bào)模式和多源觀測(cè)資料同化技術(shù)。其中,動(dòng)力環(huán)境數(shù)值預(yù)報(bào)模式主要包括全球和區(qū)域海洋環(huán)流模式、波浪模式、風(fēng)暴潮模式和海冰模式等;主要的多源觀測(cè)資料同化技術(shù)按具體同化方法區(qū)分有三維變分同化、四維變分同化和卡爾曼濾波等。近年來(lái),隨著衛(wèi)星遙感、高性能計(jì)算等技術(shù)的深入應(yīng)用,國(guó)內(nèi)外的海洋動(dòng)力環(huán)境預(yù)報(bào)技術(shù)不斷發(fā)展,預(yù)報(bào)技巧逐年提升[3-5]。海洋動(dòng)力環(huán)境預(yù)報(bào)技巧的提高主要得益于數(shù)值預(yù)報(bào)模式動(dòng)力框架和次網(wǎng)格物理過(guò)程的改進(jìn)、高性能計(jì)算技術(shù)的廣泛應(yīng)用、海洋資料同化方案的優(yōu)化和改進(jìn)等。在各類(lèi)海洋動(dòng)力環(huán)境數(shù)值模式中最核心的是全球海洋環(huán)流數(shù)值模式,它的分析和預(yù)報(bào)產(chǎn)品可為區(qū)域海洋環(huán)流模式、生態(tài)模式、波浪模式和風(fēng)暴潮模式提供動(dòng)力背景、邊界條件或者強(qiáng)迫條件等。因此,全球海洋環(huán)流數(shù)值預(yù)報(bào)模式及其配套同化技術(shù)是最重要的研究與應(yīng)用方向。

    國(guó)際上一直非常重視海洋數(shù)值預(yù)報(bào)模式和資料同化技術(shù)的發(fā)展。20世紀(jì)90年代,全球海洋環(huán)流預(yù)報(bào)系統(tǒng)已經(jīng)在英國(guó)氣象局(Met Office)的海洋同化模式預(yù)測(cè)模型(Forecasting Ocean Assimilation,Model,F(xiàn)OAM)[6]、美國(guó)海軍混合坐標(biāo)海洋模式(The HYbrid Coordinate Ocean Model,HYCOM)系統(tǒng)[7]和歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Fore-casts,ECMWF)ORAS1模式系統(tǒng)建立并運(yùn)行。2008年前后,國(guó)際上主要機(jī)構(gòu)已普遍實(shí)現(xiàn)了以數(shù)值預(yù)報(bào)為核心的預(yù)報(bào)變革,形成了完善的全球、區(qū)域海洋數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)。2011年,歐美主要預(yù)報(bào)機(jī)構(gòu)已普遍實(shí)現(xiàn)渦分辨數(shù)值預(yù)報(bào)系統(tǒng)的業(yè)務(wù)應(yīng)用。目前,全球海洋環(huán)流模式水平分辨率已經(jīng)從5°(約500 km)提高了到0.1°(約10 km),海洋環(huán)流數(shù)值預(yù)報(bào)模式系統(tǒng)正逐步朝著更高分辨率(分辨率高、網(wǎng)格精細(xì))、更多過(guò)程(方程組方程數(shù)目增加)和更快計(jì)算速度的方向發(fā)展。

    本文以國(guó)際上主要海洋預(yù)報(bào)業(yè)務(wù)中心的全球海洋動(dòng)力環(huán)境預(yù)報(bào)系統(tǒng)為主線,重點(diǎn)介紹國(guó)際上環(huán)流數(shù)值預(yù)報(bào)領(lǐng)域的主要進(jìn)展,并與國(guó)內(nèi)相關(guān)領(lǐng)域的研究進(jìn)展進(jìn)行比較,為我國(guó)未來(lái)海洋動(dòng)力環(huán)境預(yù)報(bào)系統(tǒng)的研發(fā)和應(yīng)用提供參考和借鑒。

    2 全球海洋環(huán)流預(yù)報(bào)模式

    海洋數(shù)值模式主要基于描述動(dòng)力和物理過(guò)程的偏微分方程,通過(guò)求解復(fù)雜的納維斯托克斯(Navier-Stokes)方程對(duì)海洋系統(tǒng)進(jìn)行建模和模擬。海洋環(huán)流模式、海浪模式、潮汐模式、風(fēng)暴潮模式等是納維斯托克斯方程不同的簡(jiǎn)化。

    美國(guó)海軍(Naval Oceanographic Office,NAVOCEANO)和美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction,NCEP)預(yù)報(bào)系統(tǒng)所用的全球海洋環(huán)流模式都是基于全球深度-等位密度面HYCOM開(kāi)發(fā)的,模式的分辨率設(shè)置基本一致,分辨率水平1/12°,垂直30層[7-9]。但是,兩個(gè)預(yù)報(bào)系統(tǒng)使用的海表強(qiáng)迫完全不一樣,前者采用美國(guó)海軍天氣預(yù)報(bào)模式輸出的海表通量作為海表強(qiáng)迫,后者采用NCEP的3 h業(yè)務(wù)化全球風(fēng)場(chǎng)的動(dòng)量、輻射和降水通量。未來(lái)將進(jìn)一步提高水平和垂直分辨率,水平分辨率提高到1/25°以上,垂直分辨率提高到41層,增加潮汐強(qiáng)迫,并將全球HYCOM與第三代海浪模式(WAVEWATCH III)、海冰模式以及大氣環(huán)流模式耦合。

    法國(guó)、英國(guó)和澳大利亞的全球海洋預(yù)報(bào)系統(tǒng)均基于NEMO(Nucleus for European Modelling of the Ocean)海洋模式開(kāi)發(fā),但在具體設(shè)置上有所不同。法國(guó)麥卡托國(guó)際海洋中心作為全球海洋預(yù)報(bào)系統(tǒng)的領(lǐng)先機(jī)構(gòu),全球海洋環(huán)流分析與預(yù)報(bào)系統(tǒng)體系完善(見(jiàn)圖1),承接了歐盟天氣尺度海洋預(yù)報(bào)業(yè)務(wù),其水平分辨率為1/12°、垂直方向?yàn)?0層[10];英國(guó)全球預(yù)報(bào)系統(tǒng)的水平分辨率為1/4°,垂向?yàn)?5層[6];澳大利亞全球預(yù)報(bào)系統(tǒng)則采用準(zhǔn)全球(75°S~75°N)的NEMO模式,水平分辨率為1/10°,垂直方向?yàn)?1層[11]。另外,法國(guó)的預(yù)報(bào)系統(tǒng)耦合的海冰模式是LIM(Louvain-la-Neuve Sea Ice Model),英國(guó)的預(yù)報(bào)系統(tǒng)耦合的海冰模式是CICE(Community Ice CodE and the Los Alamos Sea Ice Model),澳大利亞的預(yù)報(bào)系統(tǒng)目前則沒(méi)有耦合任何海冰模式。此外,這3個(gè)國(guó)家的預(yù)報(bào)系統(tǒng)采用的海表強(qiáng)迫場(chǎng)也不同,法國(guó)預(yù)報(bào)系統(tǒng)采用ECMWF發(fā)布的海表強(qiáng)迫場(chǎng),英國(guó)和澳大利亞則各自采用本國(guó)的中期天氣數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)輸出的海表強(qiáng)迫作為模式的上邊界條件。

    日本極為注重?cái)?shù)值預(yù)報(bào)系統(tǒng)的自主開(kāi)發(fā)和應(yīng)用,堅(jiān)持獨(dú)立發(fā)展數(shù)值模式與相關(guān)技術(shù)在實(shí)際應(yīng)用中的迭代優(yōu)化。日本氣象廳的全球海洋預(yù)報(bào)系統(tǒng)基于該國(guó)自主研制的Z坐標(biāo)海洋環(huán)流模式MRI.COM version 4。模式采用三極網(wǎng)格坐標(biāo),其水平分辨率相對(duì)較低,為1°×0.5°,垂向方向?yàn)?0層[12-14],采用了觀測(cè)的海冰密集度數(shù)據(jù)作為邊界條件。模式在海表采用了日本氣象廳發(fā)布的JRA55再分析數(shù)據(jù)進(jìn)行大氣強(qiáng)迫,并使用了水動(dòng)力模型(CaMa-Flood)計(jì)算的真實(shí)河流徑流量。

    圖1 法國(guó)海洋動(dòng)力預(yù)報(bào)系統(tǒng)及業(yè)務(wù)產(chǎn)品(引自:法國(guó)麥卡托(Mecator)海洋中心的2019年度國(guó)家報(bào)告)

    3 海洋資料同化技術(shù)

    數(shù)據(jù)同化方法可以結(jié)合觀測(cè)資料和數(shù)值模式,得到更合理的模擬結(jié)果或提高預(yù)測(cè)的精度。通常數(shù)值預(yù)報(bào)的誤差主要包括3方面來(lái)源:一是動(dòng)力模式控制方程組中包含有各種近似,或者物理過(guò)程與機(jī)制考慮不完全;二是數(shù)值模式初始場(chǎng)的質(zhì)量不高,模式參數(shù)不確定;三是模式分辨率低,不能解析次網(wǎng)格過(guò)程,而物理過(guò)程參數(shù)化方案等近似也存在誤差。因此,通過(guò)數(shù)據(jù)同化方法將觀測(cè)信息不斷地引入到數(shù)值模式,獲取可靠的數(shù)值模擬結(jié)果,或者通過(guò)改善數(shù)值模式的初始場(chǎng),提高預(yù)報(bào)和預(yù)測(cè)的準(zhǔn)確度是必不可少的業(yè)務(wù)化運(yùn)行方案。目前,國(guó)際上主要的海洋預(yù)報(bào)業(yè)務(wù)機(jī)構(gòu)通常采用三維變分同化(3DVAR)、集合最優(yōu)插值(ENOI)或者集合卡爾曼濾波(ENKF)等方法進(jìn)行資料同化。

    美國(guó)海軍全球海洋預(yù)報(bào)系統(tǒng)(Global Ocean Forecast System,GOFS)采用美國(guó)海軍耦合海洋數(shù)據(jù)同化三維變分同化(Navy Coupled Ocean Data Assimilation-Three Dimensional Variational Data Assimilation,NCODA-3DVar)[15-16]方案,同化的觀測(cè)數(shù)據(jù)包括衛(wèi)星遙感高度計(jì)(Jason-2、Altika、Cryosat-2)、衛(wèi) 星 遙 感 海 表 溫 度(NOAA-18、NOAA-19、METOP-A、METOP-B、GOES-13、GOES-15、MSG、MTSAT-2、COMS-1、WindSAT、NPP-VIIRS)、固定和漂流浮標(biāo)、現(xiàn)場(chǎng)剖面觀測(cè)(Argo、CTD、XBT、海洋Gliders)、表層溫度和鹽度船測(cè)資料、多波段微波輻射掃描儀(Spatial Sensor Microwave/Imager,SSM/I)和SSMIS海冰密集度、表層漂流計(jì)等。NCEP實(shí)時(shí)海洋預(yù)報(bào)系統(tǒng)(Real Time Ocean Forecast System,RTOFS)同化系統(tǒng)是在美國(guó)海軍NCODA-3DVAR同化方法基礎(chǔ)上建立的。主要更新和改進(jìn)包括:NCEP的RTOFS資料同化是利用混合坐標(biāo)層24 h預(yù)報(bào)作為第一猜測(cè)值,而NAVOCEANO則利用日平均場(chǎng)插值到z坐標(biāo)層上;NCEP直接利用整個(gè)水柱同化動(dòng)力高度(Altimeter Dynamic Topography,ADT),而海軍利用海面高度異常(Sea Level Anomaly,SLA)和1 000 m深的海水溫度鹽度廓線[17-18]。

    圖2 美國(guó)海軍耦合海洋資料同化方案

    法國(guó)麥卡托國(guó)際海洋中心的全球海洋預(yù)報(bào)系統(tǒng)的同化系統(tǒng)(SAM-2)主要基于降階卡曼濾波方法,是在減秩卡爾曼濾波(Singular Evolutive Extended Kalman Filter,SEEK)基礎(chǔ)上發(fā)展起來(lái)的[8,19-22]。對(duì)于預(yù)報(bào)系統(tǒng)的物理部分,SAM-2同化系統(tǒng)目前同化的數(shù)據(jù)包括海面高度計(jì)數(shù)據(jù)、海表面溫度數(shù)據(jù)、溫度與鹽度的垂直廓線數(shù)據(jù)以及海冰密集度數(shù)據(jù),其中對(duì)于背景誤差協(xié)方差矩陣采用多元三維誤差模式。該系統(tǒng)包括自適應(yīng)誤差與區(qū)域算法。自2007年4月,SAM-2同化系統(tǒng)已經(jīng)在所有的業(yè)務(wù)化系統(tǒng)中采用。該中心下一步將基于新版本NEMO3.6更新全球1/12°預(yù)報(bào)系統(tǒng),根據(jù)4D同化方案改善預(yù)報(bào)精度,同化新海表溫度觀測(cè)資料(高分辨L3 SST衛(wèi)星觀測(cè))。

    Met Office的業(yè)務(wù)系統(tǒng)采用三維變分同化方案(Nucleus for European Modelling of the Ocean/Variational Data Assimilation,NEMOVAR)。該同化方案采用了多元增量公式,是由歐洲科學(xué)計(jì)算研究 中 心(Centre Européen de Recherche et de Formation Avancée en Calcul Scientifique,CERFACS)、ECMWF、MetOffice和法國(guó)計(jì)算機(jī)科學(xué)與自動(dòng)化研究所(INRIA/LJK)等研究機(jī)構(gòu)共同研發(fā)的。NEMOVAR可以同化多種觀測(cè)資料,包括表層觀測(cè)資料(海表溫度、衛(wèi)星遙感海表高度)和次表層觀測(cè)數(shù)據(jù)(XBTs、Argo、CTDs、錨定浮標(biāo)、gliders),特別是可以合理地分辨出海洋中尺度渦旋的空間結(jié)構(gòu)。NEMOVAR還包含了一個(gè)擴(kuò)散算子,根據(jù)空間相關(guān)性傳播觀測(cè)信息,多元平衡算子考慮了海洋變量之間的協(xié)方差。Waters等[23-24]在Met Office預(yù)報(bào)系統(tǒng)中繼續(xù)發(fā)展了NEMOVAR,包括誤差協(xié)方差參數(shù)化的發(fā)展、高度計(jì)資料變分偏差校正方案[25-27]、衛(wèi)星遙感海表面溫度(Sea Surface Temperature,SST)偏差校正方案以及海冰密集度同化模塊。

    澳大利亞業(yè)務(wù)預(yù)報(bào)系統(tǒng)采用的數(shù)據(jù)同化方案為BODA[11,28],同化觀測(cè)資料包括Argo廓線資料、海水電導(dǎo)率、溫鹽深測(cè)量結(jié)果以及Jason-1 IGDR、Jason-2 IGDR、Envisat衛(wèi)星高度計(jì)和 NAVOCEANO L2P海表溫度觀測(cè)資料。

    日本氣象廳的海洋資料同化系統(tǒng)是由Fujii等[29-30]設(shè)計(jì)的3DVAR-IAU發(fā)展而來(lái)的。該同化系統(tǒng)采用非線性下降方法(POPULAR)和非線性約束,能夠高效地處理多變量(溫度和鹽度)耦合EOF垂直模態(tài)分解和背景誤差協(xié)方差矩陣中的非對(duì)角元素。該系統(tǒng)還在同化分析過(guò)程中采用了非線性和非高斯型約束,有效地改善了同化分析場(chǎng)對(duì)強(qiáng)海洋鋒空間結(jié)構(gòu)特征的刻畫(huà)。對(duì)于海冰密集度數(shù)據(jù)(MGDSST),該系統(tǒng)采用了松弛逼近和3DVAR同化技術(shù),并利用優(yōu)化之后的大氣環(huán)流場(chǎng)作為海表強(qiáng)迫進(jìn)行同化。為進(jìn)一步改善海冰同化,還開(kāi)發(fā)了基于海冰流變學(xué)的海冰伴隨模式,并進(jìn)行了一系列敏感性數(shù)值試驗(yàn)[31]。

    4 全球海洋環(huán)流預(yù)報(bào)的國(guó)際合作計(jì)劃

    與全球天氣預(yù)報(bào)業(yè)務(wù)類(lèi)似,海洋動(dòng)力環(huán)境的觀測(cè)、模式和同化技術(shù)的研發(fā)和改進(jìn)是一項(xiàng)長(zhǎng)期的、艱巨而又復(fù)雜的任務(wù),需要多學(xué)科合作以及大量人力和研究經(jīng)費(fèi)的投入,為此有必要開(kāi)展國(guó)際合作。在世界氣象組織(World Meteorological Organization,WMO)的推動(dòng)下,20世紀(jì)80年代開(kāi)展了與數(shù)值天氣預(yù)報(bào)相關(guān)的國(guó)際合作計(jì)劃,并極大地推動(dòng)了天氣預(yù)報(bào)技術(shù)的進(jìn)步和預(yù)報(bào)技巧的提高。受到數(shù)值天氣預(yù)報(bào)業(yè)務(wù)系統(tǒng)研發(fā)經(jīng)驗(yàn)的啟發(fā),國(guó)際上從20世紀(jì)90年代末就開(kāi)始在海洋動(dòng)力環(huán)境預(yù)報(bào)領(lǐng)域開(kāi)展國(guó)際合作。

    在1996—1997年期間,全球海洋數(shù)據(jù)同化試驗(yàn)(GlobalOceanDataAssimilationExperiment,GODAE)的倡議得到了各國(guó)科學(xué)家的廣泛響應(yīng),于是開(kāi)始了持續(xù)10 a(1998—2008)的GODAE國(guó)際計(jì)劃[2]。該計(jì)劃旨在:(1)將最先進(jìn)的海洋環(huán)流模式和資料同化技術(shù)應(yīng)用于短期海洋環(huán)流預(yù)報(bào)和氣候預(yù)測(cè),以及區(qū)域海洋預(yù)報(bào)系統(tǒng)的邊界條件;(2)發(fā)展海洋再分析數(shù)據(jù)集,用以理解海洋動(dòng)力過(guò)程及其可預(yù)報(bào)性,并應(yīng)用于海洋觀測(cè)系統(tǒng)的設(shè)計(jì)并提高觀測(cè)系統(tǒng)的有效性。2008年,GODAE國(guó)際合作計(jì)劃圓滿地完成了預(yù)定的各項(xiàng)研究任務(wù),極大地促進(jìn)了世界各國(guó)海洋資料同化技術(shù)和海洋預(yù)報(bào)系統(tǒng)的發(fā)展。在GODAE計(jì)劃的10 a時(shí)間里,一些發(fā)達(dá)國(guó)家業(yè)務(wù)預(yù)報(bào)中心基于水平分辨率為10 km左右的海洋模式建立了業(yè)務(wù)預(yù)測(cè)系統(tǒng),可以同化現(xiàn)有各種海表和上層海洋的觀測(cè)數(shù)據(jù),使得海洋環(huán)境預(yù)報(bào)進(jìn)入了一個(gè)可以分辨和預(yù)報(bào)中尺度海洋渦旋的時(shí)代。

    GODAE國(guó)際合作計(jì)劃結(jié)束后,各國(guó)科學(xué)家又推動(dòng)了一個(gè)新的國(guó)際合作計(jì)劃GODAE OceanView(GOV)計(jì)劃(2009—2018)[3]。該計(jì)劃可以看作是GODAE計(jì)劃的進(jìn)一步延續(xù),但是更加致力于推動(dòng)和加強(qiáng)科研成果向業(yè)務(wù)化運(yùn)行和實(shí)際應(yīng)用方向轉(zhuǎn)化,目前有包括中國(guó)在內(nèi)的15家海洋預(yù)報(bào)業(yè)務(wù)機(jī)構(gòu)成員(見(jiàn)表1)。該計(jì)劃旨在提升海洋環(huán)境監(jiān)測(cè)和預(yù)報(bào)能力,發(fā)展和完善4種能力:(1)衛(wèi)星高度計(jì)觀測(cè)海表中尺度動(dòng)力特征的能力;(2)ARGO浮標(biāo)系統(tǒng)觀測(cè)海洋溫度和鹽度剖面的能力;(3)高分辨率海洋模式模擬和預(yù)報(bào)中尺度渦旋的能力;(4)數(shù)據(jù)同化技術(shù)能夠有效地同化上述觀測(cè)數(shù)據(jù),并可以為動(dòng)力模式提供一個(gè)精確的初始場(chǎng)用于預(yù)報(bào)。目前,GOV計(jì)劃已經(jīng)順利完成,使得海洋動(dòng)力環(huán)境預(yù)報(bào)技巧得到了顯著提高。

    作為GOV計(jì)劃的后續(xù),實(shí)際上也是GODAE計(jì)劃的第三個(gè)階段,目前GOV正在過(guò)渡到Ocean Predict(2019至今,引自:https://www.godae-oceanview.org/about/gov-op-transition/),作為一個(gè)國(guó)際科學(xué)論壇和活動(dòng)平臺(tái),它將充分整合GODAE和GOV計(jì)劃的研究成果,同時(shí)還加強(qiáng)與海洋領(lǐng)域其他大型國(guó)際合作計(jì)劃(例如全球觀測(cè)系統(tǒng)計(jì)劃(Global Ocean Observation System,GOOS))的合作。更為重要的是,Ocean Predict計(jì)劃將會(huì)促進(jìn)科學(xué)研究和業(yè)務(wù)化運(yùn)行兩者之間的合作,并向其成員與合作伙伴提供專(zhuān)業(yè)知識(shí)和合作機(jī)會(huì),進(jìn)而形成一個(gè)從觀測(cè)-綜合信息分析-預(yù)報(bào)預(yù)測(cè)-科學(xué)評(píng)估-最終用戶的完整鏈條。

    5 發(fā)展趨勢(shì)

    隨著全球海洋環(huán)流預(yù)報(bào)技術(shù)的不斷完善,以及海洋觀測(cè)技術(shù)和高性能計(jì)算機(jī)的快速發(fā)展,未來(lái)有關(guān)技術(shù)主要呈現(xiàn)5方面發(fā)展趨勢(shì):

    表1 GODAE OceanView下的國(guó)際業(yè)務(wù)化海洋學(xué)機(jī)構(gòu)的預(yù)報(bào)系統(tǒng)

    (1)針對(duì)關(guān)鍵海區(qū)開(kāi)展強(qiáng)化觀測(cè)試驗(yàn)。理論研究、觀測(cè)體系的建設(shè)將更加緊密地與數(shù)值模式發(fā)展相結(jié)合,圍繞改進(jìn)數(shù)值模式的次網(wǎng)格物理參數(shù)化過(guò)程和提升模式的數(shù)值預(yù)報(bào)技巧,通過(guò)海洋環(huán)流的可預(yù)報(bào)性研究和目標(biāo)觀測(cè)研究,尋找影響海洋環(huán)流預(yù)報(bào)技巧的關(guān)鍵海區(qū),并在關(guān)鍵海區(qū)針對(duì)中小尺度動(dòng)力過(guò)程等開(kāi)展長(zhǎng)期連續(xù)的現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)。

    (2)進(jìn)一步改進(jìn)海洋模式的動(dòng)力框架并引進(jìn)多圈層相互作用。海洋數(shù)值預(yù)報(bào)模式的發(fā)展重點(diǎn)是繼續(xù)提高模式的水平和垂直分辨率、實(shí)現(xiàn)海-陸-氣-冰全耦合、物理模式與生態(tài)模式的耦合、發(fā)展自適應(yīng)非結(jié)構(gòu)網(wǎng)格、大規(guī)模并行計(jì)算技術(shù)等。數(shù)據(jù)同化發(fā)展的重點(diǎn)是實(shí)現(xiàn)變分、集合濾波及兩者混合等先進(jìn)數(shù)據(jù)同化方法和衛(wèi)星遙感觀測(cè)數(shù)據(jù)的同化應(yīng)用。

    (3)加強(qiáng)高性能計(jì)算和大數(shù)據(jù)技術(shù)的應(yīng)用。隨著海洋數(shù)值預(yù)報(bào)系統(tǒng)的計(jì)算需求越來(lái)越高,將海洋數(shù)值預(yù)報(bào)與云計(jì)算技術(shù)有機(jī)結(jié)合,突破海洋數(shù)值預(yù)報(bào)服務(wù)云和海洋數(shù)值預(yù)報(bào)平臺(tái)云構(gòu)建技術(shù),為海洋數(shù)值預(yù)報(bào)研究用戶和業(yè)務(wù)用戶提供便捷、高效、低廉的支撐環(huán)境,成為海洋數(shù)值預(yù)報(bào)的一個(gè)重要發(fā)展趨勢(shì)。

    (4)提高海洋環(huán)流預(yù)報(bào)時(shí)效。未來(lái)海洋環(huán)流的預(yù)報(bào)時(shí)效會(huì)不斷增加,目前業(yè)務(wù)化的海洋環(huán)流預(yù)報(bào)時(shí)效主要局限于短期(7 d)或者中期(7~30 d),未來(lái)預(yù)報(bào)時(shí)效將拓展至季節(jié)-年代際時(shí)間尺度,預(yù)報(bào)變量將會(huì)拓展到海冰和海表氣象要素等。預(yù)報(bào)時(shí)效的增加主要來(lái)源于數(shù)值模式性能的改進(jìn)、觀測(cè)數(shù)據(jù)的增加和資料同化技術(shù)的進(jìn)步。

    (5)發(fā)展精細(xì)化預(yù)報(bào)技術(shù)。區(qū)域海洋環(huán)流模式水平分辨率將進(jìn)一步提高,未來(lái)有可能達(dá)到1 km量級(jí),并與高精度的海浪和風(fēng)暴潮模式耦合,針對(duì)重點(diǎn)港口、海灣、重大基礎(chǔ)設(shè)施、關(guān)鍵經(jīng)濟(jì)目標(biāo)和典型人口密集區(qū)等沿海重點(diǎn)保障目標(biāo),開(kāi)展潮汐(流)、風(fēng)暴潮、海浪、海嘯、海冰、海溫等要素的預(yù)報(bào),形成對(duì)近岸工程環(huán)境的精細(xì)化預(yù)報(bào)與服務(wù)保障能力。

    6 未來(lái)發(fā)展建議

    近年來(lái),國(guó)際上數(shù)值模式技術(shù)框架已經(jīng)基本穩(wěn)定,未來(lái)發(fā)展的主要趨勢(shì)是在現(xiàn)有框架基礎(chǔ)上不斷完善物理過(guò)程、提高分辨率以及改進(jìn)同化技術(shù)。從海洋動(dòng)力環(huán)境預(yù)報(bào)效果來(lái)看,SST的預(yù)報(bào)誤差降低非常有限,均方根誤差近十年保持在0.6~1.0℃左右,近年來(lái)甚至出現(xiàn)略升[32];在表面鹽度、海面高度,及三維溫、鹽、流的長(zhǎng)期誤差檢驗(yàn)方面,受觀測(cè)時(shí)空積累限制相關(guān)研究有待進(jìn)一步加強(qiáng)。與海洋環(huán)流預(yù)報(bào)相似,其它業(yè)務(wù)化預(yù)報(bào)精度改進(jìn)也緩下步伐。以ECMWF預(yù)報(bào)產(chǎn)品為例,1992—2017年,海浪3 d預(yù)報(bào)誤差從42%下降到30%,25 a間約降12%;海浪5 d預(yù)報(bào),從50%下降到42%,25 a間約降8%[33]。同樣情況也出現(xiàn)在天氣數(shù)值預(yù)報(bào)領(lǐng)域,1982—2000年,其在北半球的500 hPa位勢(shì)高度3 d預(yù)報(bào)技巧從86%增長(zhǎng)到92%,18 a間增長(zhǎng)6%;5 d預(yù)報(bào)技巧同期增長(zhǎng)13%。2000—2015年,在北半球的3 d預(yù)報(bào)技巧從92%增長(zhǎng)到98%,增長(zhǎng)了6%;5 d預(yù)報(bào)同期增長(zhǎng)6%左右[34]。當(dāng)前海洋環(huán)境數(shù)值預(yù)報(bào)運(yùn)維復(fù)雜度逐漸增大[35],預(yù)報(bào)誤差降速明顯趨緩、數(shù)值模式技術(shù)總體框架已經(jīng)基本穩(wěn)定,結(jié)合近年來(lái)觀測(cè)手段的發(fā)展和大數(shù)據(jù)分析等新技術(shù)的應(yīng)用,對(duì)我國(guó)未來(lái)海洋動(dòng)力環(huán)境數(shù)值預(yù)報(bào)技術(shù)發(fā)展提出如下建議:

    (1)加強(qiáng)全球海洋觀測(cè)基礎(chǔ)能力建設(shè)

    發(fā)展全球海洋觀測(cè)及探測(cè)技術(shù),加強(qiáng)國(guó)際海洋領(lǐng)域的合作,完善海洋監(jiān)測(cè)/探測(cè)裝備體系,補(bǔ)齊海洋動(dòng)力環(huán)境監(jiān)測(cè)應(yīng)用傳感器技術(shù)不足。利用衛(wèi)星、飛機(jī)、船舶、浮/潛標(biāo)、無(wú)人航行器、海底觀測(cè)網(wǎng)和岸基觀測(cè)站等全球空天/水面/水下平臺(tái),建立立體連續(xù)觀測(cè)系統(tǒng)。對(duì)海洋多尺度現(xiàn)象和生態(tài)環(huán)境等進(jìn)行多參數(shù)連續(xù)實(shí)時(shí)/準(zhǔn)實(shí)時(shí)細(xì)化觀測(cè),實(shí)現(xiàn)海洋動(dòng)力環(huán)境要素的時(shí)空連續(xù)探測(cè)和信息感知。

    (2)加強(qiáng)海洋觀測(cè)數(shù)據(jù)共享

    實(shí)現(xiàn)信息融合、標(biāo)準(zhǔn)統(tǒng)一和數(shù)據(jù)共享,采用云計(jì)算、云存儲(chǔ)等技術(shù),對(duì)信息服務(wù)存儲(chǔ)系統(tǒng)軟件等進(jìn)行統(tǒng)一設(shè)計(jì)與集約化管理。支撐集約化的計(jì)算服務(wù)、存儲(chǔ)服務(wù)和應(yīng)用服務(wù),開(kāi)展國(guó)內(nèi)外海洋各類(lèi)監(jiān)測(cè)預(yù)警信息的快速處理、高效存儲(chǔ)管理和綜合分析應(yīng)用。采用按需推送、智能調(diào)取和信息交互等多種信息服務(wù)方式,整合和集成現(xiàn)有海洋信息資源,構(gòu)建共享共用的海洋信息共享平臺(tái),支持海洋信息資源共享。

    (3)加強(qiáng)核心關(guān)鍵技術(shù)攻關(guān)

    以現(xiàn)有國(guó)內(nèi)外先進(jìn)數(shù)值技術(shù)為基礎(chǔ),集合業(yè)務(wù)部門(mén)和教育科研機(jī)構(gòu)的優(yōu)勢(shì)力量,研發(fā)具有完全或高度自主知識(shí)產(chǎn)權(quán)并實(shí)現(xiàn)應(yīng)用的數(shù)值模式動(dòng)力框架、次網(wǎng)格物理過(guò)程參數(shù)化方案體系、高精度數(shù)值離散技術(shù)、精細(xì)化網(wǎng)格環(huán)境數(shù)據(jù)同化技術(shù)等,在不同尺度水平與垂直混合過(guò)程研究、高精度平流計(jì)算方案、四維變分同化和集合同化技術(shù)等細(xì)分領(lǐng)域加緊跨領(lǐng)域合作研究。

    (4)加強(qiáng)大數(shù)據(jù)、人工智能和高性能計(jì)算等技術(shù)應(yīng)用。隨著科學(xué)與技術(shù)的發(fā)展,各種觀測(cè)、模擬手段的突破,數(shù)據(jù)種類(lèi)在不斷增加、數(shù)據(jù)質(zhì)量在不斷提高。海洋數(shù)據(jù)已經(jīng)具備了大數(shù)據(jù)的典型特征(4V),即數(shù)據(jù)量大(Volume)、獲取數(shù)據(jù)速度快(Velocity)、數(shù)據(jù)類(lèi)型多樣(Variety)和數(shù)據(jù)真實(shí)性強(qiáng)(Veracity)。為此,建議在大數(shù)據(jù)與云計(jì)算技術(shù)方面,發(fā)展云架構(gòu)下海洋數(shù)值預(yù)報(bào)服務(wù)系統(tǒng)、數(shù)值模式自動(dòng)優(yōu)化分析系統(tǒng)、海洋環(huán)境大數(shù)據(jù)存儲(chǔ)云、海洋數(shù)值預(yù)報(bào)計(jì)算云、海洋大數(shù)據(jù)與人工智能預(yù)報(bào)技術(shù)。

    猜你喜歡
    環(huán)流動(dòng)力數(shù)值
    用固定數(shù)值計(jì)算
    內(nèi)環(huán)流控溫技術(shù)應(yīng)用實(shí)踐與發(fā)展前景
    學(xué)習(xí)動(dòng)力不足如何自給自足
    數(shù)值大小比較“招招鮮”
    熱鹽環(huán)流方程全局弱解的存在性
    胖胖一家和瘦瘦一家(10)
    謎底大揭秘
    動(dòng)力船
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    兩相坐標(biāo)系下MMC環(huán)流抑制策略
    激情视频va一区二区三区| 国产无遮挡羞羞视频在线观看| 久久人妻av系列| 国产aⅴ精品一区二区三区波| 校园春色视频在线观看| 国产精华一区二区三区| 午夜福利在线观看吧| 1024视频免费在线观看| 免费看十八禁软件| 中文字幕人妻熟女乱码| 久久久国产成人免费| 90打野战视频偷拍视频| ponron亚洲| 成年女人毛片免费观看观看9| 久久午夜亚洲精品久久| 久久久久九九精品影院| 国产成人精品在线电影| 国产成年人精品一区二区 | 亚洲国产毛片av蜜桃av| 美女国产高潮福利片在线看| 在线国产一区二区在线| 久久久久国产一级毛片高清牌| 欧美日韩福利视频一区二区| 国产成人精品久久二区二区91| 国产又色又爽无遮挡免费看| 成人国产一区最新在线观看| 丝袜美足系列| 校园春色视频在线观看| 国产精品秋霞免费鲁丝片| 色精品久久人妻99蜜桃| 欧美老熟妇乱子伦牲交| 国产三级在线视频| 精品福利观看| 色尼玛亚洲综合影院| 一级作爱视频免费观看| 免费搜索国产男女视频| 国产av一区二区精品久久| 午夜精品国产一区二区电影| 精品人妻1区二区| 在线观看66精品国产| 一区二区三区激情视频| 国产乱人伦免费视频| 黄网站色视频无遮挡免费观看| 日韩精品青青久久久久久| 亚洲欧美激情在线| 手机成人av网站| 国内久久婷婷六月综合欲色啪| 女警被强在线播放| 一夜夜www| www.自偷自拍.com| 超碰97精品在线观看| 级片在线观看| 一区二区三区精品91| 亚洲视频免费观看视频| 亚洲美女黄片视频| 少妇被粗大的猛进出69影院| 欧美在线黄色| www.自偷自拍.com| 在线观看免费视频日本深夜| 好看av亚洲va欧美ⅴa在| 男女高潮啪啪啪动态图| 无人区码免费观看不卡| 久久影院123| 欧美色视频一区免费| cao死你这个sao货| 国产黄a三级三级三级人| 欧美在线黄色| 少妇粗大呻吟视频| 国产熟女xx| 亚洲精品国产色婷婷电影| 黄色片一级片一级黄色片| 亚洲七黄色美女视频| 天堂动漫精品| av在线播放免费不卡| 亚洲成人久久性| 久久香蕉精品热| 欧美乱色亚洲激情| 丝袜美足系列| 国产在线精品亚洲第一网站| 色综合欧美亚洲国产小说| 久99久视频精品免费| 丝袜人妻中文字幕| 色精品久久人妻99蜜桃| 亚洲九九香蕉| 久久中文看片网| 丰满的人妻完整版| 国产主播在线观看一区二区| 男人舔女人下体高潮全视频| 久久久久亚洲av毛片大全| 亚洲精品在线美女| 九色亚洲精品在线播放| 激情在线观看视频在线高清| 午夜91福利影院| 亚洲一区二区三区不卡视频| www.熟女人妻精品国产| 亚洲欧美精品综合久久99| 国产亚洲欧美精品永久| 成人亚洲精品av一区二区 | 超碰成人久久| netflix在线观看网站| 看黄色毛片网站| 在线观看舔阴道视频| 久久人人爽av亚洲精品天堂| 中文字幕高清在线视频| 搡老熟女国产l中国老女人| 热re99久久国产66热| 久久国产亚洲av麻豆专区| 久久影院123| 精品午夜福利视频在线观看一区| 电影成人av| 日韩欧美免费精品| 啦啦啦 在线观看视频| 国产三级黄色录像| 丰满的人妻完整版| 一级a爱片免费观看的视频| 夜夜看夜夜爽夜夜摸 | 色精品久久人妻99蜜桃| 国产午夜精品久久久久久| 精品久久久久久,| 亚洲成人免费av在线播放| 激情视频va一区二区三区| 九色亚洲精品在线播放| 热99国产精品久久久久久7| 国产激情久久老熟女| 大码成人一级视频| 18禁黄网站禁片午夜丰满| 亚洲精品国产精品久久久不卡| 午夜福利一区二区在线看| 巨乳人妻的诱惑在线观看| 亚洲专区国产一区二区| 欧美 亚洲 国产 日韩一| 午夜免费鲁丝| 日韩大尺度精品在线看网址 | 女人高潮潮喷娇喘18禁视频| 最近最新中文字幕大全电影3 | 中文字幕人妻丝袜制服| 嫩草影视91久久| 国产激情欧美一区二区| 亚洲人成电影免费在线| 国产精品 国内视频| 国产黄色免费在线视频| 日韩av在线大香蕉| 99国产精品一区二区蜜桃av| 在线观看日韩欧美| 黄色片一级片一级黄色片| 搡老乐熟女国产| 人人澡人人妻人| 天堂影院成人在线观看| 久久久国产成人精品二区 | 欧美激情久久久久久爽电影 | 亚洲精品久久午夜乱码| 校园春色视频在线观看| 欧美一级毛片孕妇| 国产片内射在线| 中出人妻视频一区二区| 久久草成人影院| 欧美另类亚洲清纯唯美| 亚洲免费av在线视频| 十八禁人妻一区二区| 亚洲精品久久成人aⅴ小说| 人妻久久中文字幕网| 多毛熟女@视频| 级片在线观看| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区精品| 久久久久久久精品吃奶| 1024香蕉在线观看| 免费少妇av软件| 伊人久久大香线蕉亚洲五| 久久国产精品影院| 韩国av一区二区三区四区| 窝窝影院91人妻| 搡老乐熟女国产| 国产亚洲欧美精品永久| 久久久国产成人精品二区 | 国产黄色免费在线视频| 99国产综合亚洲精品| 久久香蕉国产精品| 大陆偷拍与自拍| 露出奶头的视频| 长腿黑丝高跟| 黄色视频不卡| 欧美黄色片欧美黄色片| 国产男靠女视频免费网站| 亚洲精品在线观看二区| 欧美日韩亚洲高清精品| av天堂在线播放| 好男人电影高清在线观看| 亚洲全国av大片| 涩涩av久久男人的天堂| 国产精品香港三级国产av潘金莲| 亚洲伊人色综图| 又紧又爽又黄一区二区| 变态另类成人亚洲欧美熟女 | 18禁黄网站禁片午夜丰满| 黄片大片在线免费观看| 午夜免费激情av| 国产精华一区二区三区| 国产成人精品无人区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久大精品| 俄罗斯特黄特色一大片| 国产成人精品在线电影| 免费女性裸体啪啪无遮挡网站| 大型av网站在线播放| 欧美日本中文国产一区发布| 黄片小视频在线播放| 在线免费观看的www视频| www.色视频.com| 国内揄拍国产精品人妻在线| 美女免费视频网站| 精品久久久久久,| 免费看a级黄色片| 精品一区二区三区人妻视频| .国产精品久久| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 变态另类成人亚洲欧美熟女| 亚洲专区中文字幕在线| 人人妻,人人澡人人爽秒播| 国产91精品成人一区二区三区| 日本 av在线| 国产伦精品一区二区三区视频9| 男女床上黄色一级片免费看| 国产一区二区激情短视频| 听说在线观看完整版免费高清| 丰满人妻一区二区三区视频av| 五月玫瑰六月丁香| 久久久久久国产a免费观看| 精品久久国产蜜桃| 一个人免费在线观看的高清视频| 国产精品日韩av在线免费观看| 变态另类丝袜制服| 最近视频中文字幕2019在线8| 制服丝袜大香蕉在线| 男女之事视频高清在线观看| 亚洲欧美日韩卡通动漫| av视频在线观看入口| 日韩av在线大香蕉| 18禁裸乳无遮挡免费网站照片| 亚洲国产色片| 午夜激情欧美在线| 欧美极品一区二区三区四区| 高清在线国产一区| 亚洲第一电影网av| 欧美日韩亚洲国产一区二区在线观看| 国产单亲对白刺激| 国产久久久一区二区三区| 高清在线国产一区| 91九色精品人成在线观看| 天堂av国产一区二区熟女人妻| 免费av观看视频| 欧美日本视频| 免费看光身美女| 国产精品乱码一区二三区的特点| 国产成+人综合+亚洲专区| 少妇熟女aⅴ在线视频| 成人高潮视频无遮挡免费网站| 国产探花极品一区二区| 伊人久久精品亚洲午夜| h日本视频在线播放| 夜夜看夜夜爽夜夜摸| 亚洲精品一区av在线观看| 久久欧美精品欧美久久欧美| 色哟哟哟哟哟哟| 永久网站在线| 色哟哟哟哟哟哟| 91字幕亚洲| av在线观看视频网站免费| 婷婷色综合大香蕉| 香蕉av资源在线| 成人特级av手机在线观看| 一个人免费在线观看的高清视频| 亚洲国产精品成人综合色| 搡老妇女老女人老熟妇| www.www免费av| 久9热在线精品视频| 亚洲一区高清亚洲精品| 窝窝影院91人妻| 九九在线视频观看精品| 精品一区二区三区人妻视频| 神马国产精品三级电影在线观看| eeuss影院久久| 国产精品三级大全| 久久性视频一级片| 精品人妻一区二区三区麻豆 | 日本 欧美在线| 精品人妻熟女av久视频| 国产高清三级在线| 国产精品不卡视频一区二区 | 少妇丰满av| 一区二区三区四区激情视频 | 免费人成视频x8x8入口观看| 岛国在线免费视频观看| 老鸭窝网址在线观看| 搡老妇女老女人老熟妇| 日韩亚洲欧美综合| 色视频www国产| 欧美高清成人免费视频www| 国产精品国产高清国产av| 成人特级黄色片久久久久久久| 午夜福利18| 日韩欧美三级三区| 极品教师在线免费播放| 午夜福利视频1000在线观看| 如何舔出高潮| 色5月婷婷丁香| 国产精品一区二区性色av| 美女大奶头视频| 久久中文看片网| 国产av不卡久久| 一级毛片久久久久久久久女| 91在线精品国自产拍蜜月| 亚洲黑人精品在线| 简卡轻食公司| 首页视频小说图片口味搜索| 男人狂女人下面高潮的视频| 国产av一区在线观看免费| 久久久久久久久中文| netflix在线观看网站| 午夜免费男女啪啪视频观看 | 日本与韩国留学比较| 亚洲成人中文字幕在线播放| 最新中文字幕久久久久| 97超视频在线观看视频| 色综合站精品国产| 午夜亚洲福利在线播放| 91在线精品国自产拍蜜月| 国产黄a三级三级三级人| 88av欧美| 亚洲精品一区av在线观看| 亚洲国产精品sss在线观看| 男女床上黄色一级片免费看| 欧美一区二区亚洲| 国产亚洲精品久久久com| 非洲黑人性xxxx精品又粗又长| 如何舔出高潮| 中文字幕av在线有码专区| 免费无遮挡裸体视频| 在线观看av片永久免费下载| 免费电影在线观看免费观看| 51午夜福利影视在线观看| 禁无遮挡网站| АⅤ资源中文在线天堂| av视频在线观看入口| 美女 人体艺术 gogo| 国产日本99.免费观看| 欧美高清性xxxxhd video| 韩国av一区二区三区四区| 亚洲中文字幕日韩| 欧美bdsm另类| 欧美性猛交黑人性爽| 国产精品日韩av在线免费观看| 18美女黄网站色大片免费观看| 国产日本99.免费观看| 欧美国产日韩亚洲一区| 91在线观看av| 久久香蕉精品热| 小说图片视频综合网站| 亚洲国产精品sss在线观看| 男女床上黄色一级片免费看| 国产一区二区三区在线臀色熟女| 午夜老司机福利剧场| 亚洲精品粉嫩美女一区| 久久午夜福利片| 成人鲁丝片一二三区免费| 国产大屁股一区二区在线视频| 国产av麻豆久久久久久久| 亚洲专区国产一区二区| 色在线成人网| www.999成人在线观看| 色精品久久人妻99蜜桃| 可以在线观看的亚洲视频| 欧美高清性xxxxhd video| 国产伦在线观看视频一区| 国产精品综合久久久久久久免费| 深夜精品福利| 性插视频无遮挡在线免费观看| 国产视频内射| 亚洲av不卡在线观看| 日韩欧美在线乱码| 欧美日韩瑟瑟在线播放| 少妇人妻一区二区三区视频| 听说在线观看完整版免费高清| 午夜老司机福利剧场| 国产精品永久免费网站| 欧美不卡视频在线免费观看| 美女大奶头视频| 精品久久久久久久久av| 色精品久久人妻99蜜桃| 欧美中文日本在线观看视频| 直男gayav资源| 亚洲精品在线美女| 一级a爱片免费观看的视频| 欧美一区二区国产精品久久精品| 色精品久久人妻99蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 搡老熟女国产l中国老女人| 国产乱人视频| 偷拍熟女少妇极品色| 久久99热6这里只有精品| 丁香六月欧美| 永久网站在线| 欧美午夜高清在线| 久久久久九九精品影院| 日本在线视频免费播放| 一个人看视频在线观看www免费| 国产精品综合久久久久久久免费| 欧美黑人欧美精品刺激| 亚洲精品久久国产高清桃花| 蜜桃亚洲精品一区二区三区| 欧美潮喷喷水| 国产免费一级a男人的天堂| 热99在线观看视频| 床上黄色一级片| 午夜精品久久久久久毛片777| 国产日本99.免费观看| 免费一级毛片在线播放高清视频| 女同久久另类99精品国产91| 午夜a级毛片| 久久久久国内视频| 色精品久久人妻99蜜桃| 美女高潮的动态| 国产伦精品一区二区三区四那| 99视频精品全部免费 在线| 久久精品国产清高在天天线| 一区二区三区免费毛片| 亚洲天堂国产精品一区在线| 三级男女做爰猛烈吃奶摸视频| 欧美黑人欧美精品刺激| 99精品久久久久人妻精品| 久久久久久久亚洲中文字幕 | 国产av在哪里看| 少妇裸体淫交视频免费看高清| 国产日本99.免费观看| 波多野结衣巨乳人妻| 国产精品综合久久久久久久免费| 久久久久久大精品| 久久久国产成人精品二区| 老司机深夜福利视频在线观看| 日韩欧美三级三区| 三级毛片av免费| a级一级毛片免费在线观看| 嫩草影视91久久| 18禁在线播放成人免费| 97超级碰碰碰精品色视频在线观看| 在线十欧美十亚洲十日本专区| 国产精品,欧美在线| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 色在线成人网| 欧美性感艳星| 日本黄色片子视频| 韩国av一区二区三区四区| 亚洲真实伦在线观看| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 国产白丝娇喘喷水9色精品| 香蕉av资源在线| 老鸭窝网址在线观看| 国产精品美女特级片免费视频播放器| 9191精品国产免费久久| 国产精品久久久久久亚洲av鲁大| 国产av不卡久久| 免费在线观看亚洲国产| 在线免费观看的www视频| 少妇人妻精品综合一区二区 | 成年女人毛片免费观看观看9| 在线国产一区二区在线| 国产野战对白在线观看| 国产亚洲欧美98| 国产欧美日韩一区二区精品| 国产精华一区二区三区| 亚洲人成网站在线播| .国产精品久久| 国产精品av视频在线免费观看| 国模一区二区三区四区视频| 乱码一卡2卡4卡精品| 亚洲电影在线观看av| 人妻夜夜爽99麻豆av| 国产欧美日韩一区二区三| 1024手机看黄色片| 91九色精品人成在线观看| 少妇丰满av| www.www免费av| 亚洲av中文字字幕乱码综合| 男人舔奶头视频| 国产成人福利小说| 欧美色欧美亚洲另类二区| 黄色一级大片看看| 国产成年人精品一区二区| 麻豆国产av国片精品| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 我要看日韩黄色一级片| 成人亚洲精品av一区二区| 国产探花极品一区二区| 18禁黄网站禁片免费观看直播| 日本一二三区视频观看| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 国产国拍精品亚洲av在线观看| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 国语自产精品视频在线第100页| 在线观看一区二区三区| 国产亚洲精品久久久com| 欧美日韩亚洲国产一区二区在线观看| 国产毛片a区久久久久| 久久精品影院6| 麻豆成人午夜福利视频| 欧美国产日韩亚洲一区| 小蜜桃在线观看免费完整版高清| 少妇的逼水好多| 久久久成人免费电影| 美女 人体艺术 gogo| 亚洲国产精品合色在线| 欧美日本亚洲视频在线播放| 日本 欧美在线| 小蜜桃在线观看免费完整版高清| ponron亚洲| 免费大片18禁| 欧美中文日本在线观看视频| av福利片在线观看| 88av欧美| 中文字幕免费在线视频6| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 别揉我奶头 嗯啊视频| 在线观看av片永久免费下载| 高清毛片免费观看视频网站| 成人国产综合亚洲| 噜噜噜噜噜久久久久久91| 夜夜躁狠狠躁天天躁| 国产免费av片在线观看野外av| 欧美黄色淫秽网站| 中文字幕av成人在线电影| 男人的好看免费观看在线视频| 伦理电影大哥的女人| 99国产综合亚洲精品| 免费在线观看日本一区| 我的女老师完整版在线观看| 亚洲avbb在线观看| 亚洲成人精品中文字幕电影| 18禁黄网站禁片午夜丰满| 亚洲第一区二区三区不卡| www.www免费av| 亚洲av五月六月丁香网| 观看美女的网站| 18+在线观看网站| 色噜噜av男人的天堂激情| 亚洲精品影视一区二区三区av| 亚洲精品乱码久久久v下载方式| 麻豆久久精品国产亚洲av| 亚洲最大成人av| 国产在视频线在精品| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放| 蜜桃久久精品国产亚洲av| 欧美又色又爽又黄视频| 国产精品久久久久久人妻精品电影| 特大巨黑吊av在线直播| 婷婷亚洲欧美| 日韩欧美三级三区| 美女高潮的动态| 国产午夜福利久久久久久| 一夜夜www| 脱女人内裤的视频| 日韩免费av在线播放| 国产精品久久久久久人妻精品电影| 亚洲av成人精品一区久久| 97碰自拍视频| 亚洲欧美日韩卡通动漫| 波多野结衣巨乳人妻| 日韩人妻高清精品专区| 夜夜爽天天搞| 赤兔流量卡办理| 日韩欧美精品v在线| 一本一本综合久久| 国产免费男女视频| 国产高清激情床上av| 国内精品美女久久久久久| 一区二区三区高清视频在线| 亚洲无线在线观看| 精品久久久久久久久久久久久| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 日韩大尺度精品在线看网址| 丰满的人妻完整版| 亚洲欧美日韩高清专用| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆| 51国产日韩欧美| 成年女人看的毛片在线观看| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式| 好男人电影高清在线观看| 久久性视频一级片| 亚洲在线自拍视频| 日本与韩国留学比较| 久久午夜亚洲精品久久| 少妇的逼水好多| 亚洲av二区三区四区| 18禁在线播放成人免费| 日日夜夜操网爽| 精品午夜福利在线看| 国产单亲对白刺激| 亚洲人成网站在线播放欧美日韩| av在线天堂中文字幕| 美女高潮的动态| 1000部很黄的大片| av女优亚洲男人天堂| 99视频精品全部免费 在线| 我要看日韩黄色一级片| 免费一级毛片在线播放高清视频|