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

    黃海綠潮分布年際變化特征分析

    2016-12-22 01:53:53郭偉趙亮李秀梅
    海洋學(xué)報(bào) 2016年12期
    關(guān)鍵詞:綠潮覆蓋面積年際

    郭偉,趙亮*,李秀梅

    (1. 天津科技大學(xué) 海洋與環(huán)境學(xué)院,天津 300457)

    ?

    黃海綠潮分布年際變化特征分析

    郭偉1,趙亮1*,李秀梅1

    (1. 天津科技大學(xué) 海洋與環(huán)境學(xué)院,天津 300457)

    2008—2015年,連續(xù)8 年在黃海海域暴發(fā)大規(guī)模綠潮,但因暴發(fā)時(shí)間、規(guī)模及漂移路徑的不同,對(duì)沿海地區(qū)造成的環(huán)境影響和經(jīng)濟(jì)損失大不相同。本文利用EOS/MODIS衛(wèi)星的多通道資料,采用NDVI算法獲取綠潮信息,給出了2007年以來5—8月所有無云或少云晴朗天氣下黃海海面綠潮的分布情況。黃海綠潮覆蓋面積變化每年呈現(xiàn)單峰值分布,MODIS衛(wèi)星在5月中旬至6月初首次發(fā)現(xiàn)綠潮,隨后30~40 d內(nèi)達(dá)到覆蓋面積峰值,7月快速消退,8月上旬完全消失。2007年綠潮出現(xiàn)首年覆蓋面積極低,發(fā)現(xiàn)綠潮時(shí)面積均未超過50 km2;2008年和2009年綠潮覆蓋面積峰值分別為3 110 km2和4 075 km2,自此綠潮災(zāi)害成為新的海洋環(huán)境事件;2010—2012年綠潮暴發(fā)規(guī)模異常低值,各年覆蓋面積峰值均未超過1 800 km2,但從2013年開始綠潮規(guī)模逐漸攀升,到2015年綠潮覆蓋面積峰值達(dá)5 629 km2,持續(xù)可達(dá)98 d,覆蓋面積和持續(xù)時(shí)間為歷年之最。綠潮漂移路徑可歸納分為3類:2008年和2011年綠潮主體先北向漂移越過34°30′N后,西北向垂直于岸線漂移,主要對(duì)連云港、日照和青島造成較大影響;2009年和2012年綠潮主體先北向漂移越過35°N后,東北向平行于岸線漂移,故只有少許綠潮上岸;2010年、2013—2015年綠潮主體北向漂移至近岸后東北向沿岸漂移,對(duì)日照、青島和榮成沿海造成大面積影響。所有年份綠潮影響范圍均限于南黃海內(nèi),東側(cè)邊界最遠(yuǎn)未越過124.2°E。

    綠潮;覆蓋面積;漂移路徑;年際變化;黃海

    1 引言

    大型海洋綠藻過度增殖引起的一種新的海洋生態(tài)災(zāi)害現(xiàn)象,被稱為綠潮(Green Tide)。綠潮已成為世界范圍性的河口、潟湖和內(nèi)灣等海域的一個(gè)普遍現(xiàn)象。自20世紀(jì)70年代以來,綠潮在世界范圍內(nèi)的發(fā)生頻率和影響規(guī)模呈現(xiàn)上升趨勢(shì),美國、加拿大、丹麥、荷蘭、法國、意大利、日本、韓國和菲律賓等國家均有綠潮災(zāi)害報(bào)道,其中以法國布列塔尼沿海的情況最為嚴(yán)重[1—4]。我國自2008年開始連續(xù)8年在黃海海域暴發(fā)大規(guī)模綠潮,對(duì)水產(chǎn)養(yǎng)殖、濱海旅游、海上交通運(yùn)輸及海上賽事等相關(guān)產(chǎn)業(yè)影響嚴(yán)重,給江蘇及山東沿海城市帶來巨大的經(jīng)濟(jì)損失和社會(huì)影響。因此,開展關(guān)于黃海綠潮的溯源、分布變化規(guī)律和發(fā)生機(jī)理的研究,對(duì)區(qū)域經(jīng)濟(jì)穩(wěn)定發(fā)展具有重要意義。

    不同學(xué)者基于衛(wèi)星遙感和數(shù)值模擬分析了黃海綠潮分布特征,遙感反演發(fā)現(xiàn)黃海綠潮多呈條帶狀分布[5],其分布受到風(fēng)和海流的共同影響。李曰嵩等[6]利用海洋模型FVCOM對(duì)綠潮早期發(fā)生地的蘇北海域進(jìn)行水動(dòng)力模擬,發(fā)現(xiàn)潮流的作用導(dǎo)致早期綠潮聚集,呈現(xiàn)與輻射沙洲相一致放射狀的條帶狀或斑塊狀分布;張?zhí)K平等[7]和衣立等[8]從黃海水文氣象環(huán)境方面分別對(duì)2008年和2009年綠潮聚集和定向移動(dòng)進(jìn)行分析,指出綠潮聚集可能受到輻合風(fēng)場的直接驅(qū)動(dòng),綠潮聚集區(qū)的移動(dòng)與盛行風(fēng)向一致;喬方利等[9]基于MASNUM海浪-潮流-環(huán)流耦合數(shù)值模式,發(fā)現(xiàn)風(fēng)場驅(qū)動(dòng)下的海洋表層流場是綠潮漂移路徑變化的主要原因;Keesing等[10]基于MODIS圖像給出了2007—2009年黃海綠潮的覆蓋面積,指出風(fēng)場是影響黃海綠潮漂移路徑的主要原因。黃娟等[11]給出2008—2013年黃海綠潮首次出現(xiàn)、最大覆蓋及分布面積的變化,指出綠潮分布和漂移路徑存在明顯的年際變化特征。目前國內(nèi)外學(xué)者的研究主要聚焦于黃海綠潮的溯源、監(jiān)測方法和分布特征等方面,對(duì)覆蓋面積和漂移路徑的研究多局限在特定年份,對(duì)綠潮分布特征的年際變化研究較少。

    在綠潮的暴發(fā)和持續(xù)階段,使用MODIS衛(wèi)星數(shù)據(jù)可以較為準(zhǔn)確的獲取綠潮分布及覆蓋面積[4—5,8,10],這彌補(bǔ)了傳統(tǒng)的調(diào)查船耗時(shí)費(fèi)力,覆蓋不全面的缺陷。本文基于MODIS衛(wèi)星數(shù)據(jù),重點(diǎn)分析2007—2015年黃海綠潮生消過程、覆蓋面積、漂移路徑的年際變化特征,以期為黃海綠潮的預(yù)測、預(yù)警和防控提供科學(xué)研究基礎(chǔ)。

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

    數(shù)據(jù)來源于NASA數(shù)據(jù)共享網(wǎng)站提供的EOS衛(wèi)星AUQA和TEREA的中等分辨率成像光譜儀(MODIS)資料(https://ladsweb.nascom.nasa.gov/data/search.html)。MODIS是新一代“圖譜合一”的光學(xué)遙感儀器,具有電磁波譜范圍在0.4~14.5 μm之間的36個(gè)離散波段通道。MODIS地面最高分辨率為250 m,能夠覆蓋整個(gè)黃海,且每天過境兩次,是一種可確保連續(xù)觀測的數(shù)據(jù)源。文獻(xiàn)報(bào)道,綠潮自2007年開始出現(xiàn)[5,12],每年4—7月為綠潮暴發(fā)和持續(xù)期,8—10月為綠潮消亡期[13]。早期漂浮綠藻和消亡期懸浮綠藻可利用船舶監(jiān)測,而只有當(dāng)綠潮聚集到一定規(guī)模才可被MODIS發(fā)現(xiàn),故本文關(guān)注時(shí)間為2007—2015年每年的5—8月,觀測海域范圍:31.5°~37.5°N,119°~124.5°E。本文選擇下載無云或少云條件下空間分辨率為250 m的MODIS 1B級(jí)HDF文件,使用其中的紅外波段(波段1:620~670 nm)和近紅外波段(波段2:841~876 nm)數(shù)據(jù),利用歸一化植被指數(shù)(NDVI)提取綠潮信息[7—8,10]:

    NDVI=(NIR-R)/(NIR+R),

    (1)

    式中,NIR和R分別為近紅外波段和紅外波段的反射率。綠潮水體近紅外波段反射率高,紅外波段反射率較低,而正常海水在紅外和近紅外波段的反射率相當(dāng),因此綠潮聚集區(qū)的NDVI有較高的值。此外利用紅外、近紅外和紅外波段合成的RGB假彩色影像圖上,海水常呈現(xiàn)深藍(lán)色或黑色,綠潮水體常呈現(xiàn)翠綠色,陸地植被則呈現(xiàn)黃綠色,三者具有較明顯差異[14]。利用NDVI值和假彩色影像能夠明顯區(qū)分綠潮水體及正常海水和陸地植被。

    利用IDL語言編程批量計(jì)算NDVI,通過設(shè)置合理的NDVI閾值T提取綠潮信息,獲取綠潮二值圖。將綠潮二值圖疊加到RGB假彩色合成圖進(jìn)行目視解譯[15—16],進(jìn)一步驗(yàn)證綠潮信息提取的正確性,剔除由于島嶼、薄云等造成的誤判,最后將確認(rèn)后的結(jié)果轉(zhuǎn)成矢量圖輸入到ArcGIS軟件,統(tǒng)計(jì)研究區(qū)NDVI>T的像元個(gè)數(shù)N,計(jì)算綠潮覆蓋面積S[10]:

    (2)

    式中,r為MODIS衛(wèi)星空間分辨率250 m。

    3 結(jié)果與討論

    利用MODIS衛(wèi)星數(shù)據(jù)反演了2007—2015年9年內(nèi)所有無云或少云天氣下黃海綠潮信息,基于反演結(jié)果分析黃海綠潮生消過程的變化特征,并對(duì)綠潮的覆蓋面積、分布特征及漂移路徑進(jìn)行統(tǒng)計(jì)分析,給出黃海綠潮年際差異及變化趨勢(shì)。

    3.1 黃海綠潮生消過程年際變化特征

    選擇2006—2015年每年5—8月所有無云或少云晴朗天氣的MODIS衛(wèi)星圖像,通過目視解譯方法鑒別挑選這些天中存在綠潮的圖像(表1)。統(tǒng)計(jì)發(fā)現(xiàn),2006年5—8月無云晴朗天總數(shù)為14 d,時(shí)間上較為均勻分布在5—8月,均未發(fā)現(xiàn)綠潮,故2006年黃海并未出現(xiàn)綠潮,而2005年也未發(fā)現(xiàn)綠潮存在。因此,2007年首次通過MODIS衛(wèi)星監(jiān)測到黃海綠潮的發(fā)生。

    表1 2006—2015年5—8月無云及綠潮天數(shù)統(tǒng)計(jì)表

    2007年黃海綠潮發(fā)生規(guī)模較小,僅在6月17日、7月16日和7月18日發(fā)現(xiàn)零星條帶狀分布,且覆蓋面積均不足50 km2(圖1),很難用于分析覆蓋面積和漂移路徑的演變規(guī)律,故不納入年際變化的分析中。2008年黃海發(fā)生大規(guī)模綠潮,一度影響青島奧帆賽的正常舉辦,并對(duì)青島、日照和海陽等沿海城市造成巨大的經(jīng)濟(jì)損失和社會(huì)環(huán)境影響[17],因而常常認(rèn)定2008年為黃海綠潮大規(guī)模暴發(fā)的第一年。統(tǒng)計(jì)2008—2015年5—8月黃海無云或少云條件下存在綠潮天數(shù)列于表1,各年無云晴朗總天數(shù)在12~24 d之間,MODIS發(fā)現(xiàn)綠潮天數(shù)為9~16 d,其中在規(guī)模性綠潮形成和持續(xù)階段(6—7月)的每個(gè)月發(fā)現(xiàn)綠潮天數(shù)都在3 d以上,這樣的綠潮信息量可以支持我們的研究。

    圖1 2007年黃海綠潮分布圖Fig.1 Distributions of the Green Tide in the Yellow Sea in 2007

    對(duì)歷年綠潮分布特征(包括首次發(fā)現(xiàn)時(shí)間、最大面積、增值時(shí)間等)進(jìn)行統(tǒng)計(jì)(表2),發(fā)現(xiàn)除2012年綠潮出現(xiàn)時(shí)間較早,其他年份首次出現(xiàn)時(shí)間集中在5月中旬到6月初。2008年和2009年發(fā)現(xiàn)綠潮時(shí)間與Keesing等[10]的研究結(jié)果完全一致,衣立等[8]指出2009年黃海西部降水明顯增多的時(shí)間較2008年晚15 d左右,是2009年綠潮暴發(fā)比2008年晚15 d左右的重要原因;2010年發(fā)現(xiàn)綠潮時(shí)間最晚,但比喬方利等[9]研究結(jié)果早4 d,且分布范圍要??;2011年和2013年首次發(fā)現(xiàn)時(shí)間與黃娟等[11]的研究結(jié)果基本一致。

    表2 2008—2015年黃海綠潮發(fā)生情況統(tǒng)計(jì)表

    圖2 2008—2015年MODIS衛(wèi)星首次發(fā)現(xiàn)黃海綠潮分布Fig.2 Distribution of the Green Tide appeared first time in the Yellow Sea from 2008 to 2015 by MODIS

    反演得到綠潮首次出現(xiàn)時(shí)的分布顯示,各年早期綠潮多以條帶或斑塊狀分布,分布位置主要集中在3片區(qū)域(圖2)。區(qū)域A:蘇北淺灘輻射沙脊群。2009年和2012年5月首次在該區(qū)域發(fā)現(xiàn)條帶狀綠潮,其分布形狀和輻射沙脊的地形類似。區(qū)域B:中心距連云港110 km海域(33.7°~35°N,120.5°~121.5°E),2008年、2010年、2011年、2014年及2015年首次發(fā)現(xiàn)黃海綠潮區(qū)域。Keesing等[10]通過對(duì)比2008年和2009年綠潮漂移的時(shí)間和位置,認(rèn)為2008年的發(fā)源地很可能和2009年一樣,東南風(fēng)驅(qū)動(dòng)使得2008年早期綠潮從A區(qū)漂向B區(qū)。對(duì)比發(fā)現(xiàn),2009年5月26日出現(xiàn)在A區(qū)的綠潮,在6月14日運(yùn)動(dòng)到B區(qū),2012年5月6日出現(xiàn)在A區(qū)的綠潮,在5月16日運(yùn)動(dòng)到B區(qū)。首先在A區(qū)出現(xiàn)的綠潮都會(huì)流經(jīng)B區(qū),推測首先在B區(qū)出現(xiàn)的綠潮在衛(wèi)星可監(jiān)測前也流經(jīng)A區(qū),隨后受夏季東南風(fēng)、蘇北沿岸流及長江沖淡水等共同作用漂移至B區(qū)。故推測A區(qū)和B區(qū)的綠潮可能都來自A區(qū)最南端的如東沿?;蛘哂筛系钠鹪吹仄贫鴣怼^(qū)域C:中心在如東以東230 km海域(31.6°~33.5°N,123.2°~124.2°E),是2013年和2014年首次發(fā)現(xiàn)黃海綠潮區(qū)域。2013年5月13日在該片海域內(nèi)發(fā)現(xiàn)綠潮,5月21日衛(wèi)星顯示該片海域已無綠潮,而在區(qū)域A和B處出現(xiàn)覆蓋面積為276 km2的條帶狀綠潮,分布形狀和輻射沙脊較為一致;2014年5月12日在區(qū)域C和B同時(shí)出現(xiàn)綠潮,5月22日,區(qū)域C綠潮消失,區(qū)域B綠潮則迅速生長北移。南黃海5月上旬主要以東南風(fēng)為主,C區(qū)的綠潮不可能由A區(qū)綠潮漂移所致;而在短短的10 d時(shí)間內(nèi),C區(qū)綠潮也無法漂移至A區(qū)和B區(qū)。推測C區(qū)出現(xiàn)綠潮與A區(qū)和B區(qū)不同,可能是由更南的起源地漂移而來,而C區(qū)5月上旬海溫在10~13℃,低于綠潮最適生長溫度15~25℃[18],不適合綠潮生長,故很快消失了,區(qū)域C綠潮的來源和消衰還有待于進(jìn)一步確認(rèn)。

    黃海綠潮自每年6月進(jìn)入快速增殖期,一般30~40 d時(shí)間內(nèi)達(dá)到本年覆蓋面積峰值,到8月中旬基本消失(表2)。黃海綠潮整個(gè)生消過程持續(xù)2~3個(gè)月時(shí)間,2009年至2011年綠潮生消過程時(shí)間較短,在70~80 d左右,其余5年在90~100 d左右。2008—2015年黃海綠潮無論出現(xiàn)時(shí)間早晚,最后發(fā)現(xiàn)綠潮消亡時(shí)間均集中在8月上半月,且當(dāng)天綠潮以零星斑塊狀出現(xiàn),覆蓋面積在10~40 km2之間。每年黃海北部8月中上旬海表溫度在26~28℃,超出綠潮的適宜生長溫度,與每年綠潮的消亡時(shí)間較為一致,推測黃海綠潮的消亡受溫度影響顯著,海表溫度迅速上升超出綠潮生長的適宜溫度是綠潮快速消亡的重要原因之一。

    3.2 黃海綠潮面積及分布年際變化特征

    圖3給出了2007—2015年9 年內(nèi)所有無云或少云晴朗天氣下綠潮覆蓋面積隨時(shí)間的變化。除2010年綠潮覆蓋面積變化較為特殊,呈現(xiàn)雙峰值分布外,其他年份均呈單峰值分布:5月中旬到6月初首次發(fā)現(xiàn)后,6月經(jīng)歷快速增殖,在30~40 d時(shí)間內(nèi)達(dá)到面積峰值,7月份快速消退,到8月上旬完全消失。

    圖3 2007—2015年黃海無云或少云晴朗天氣綠潮覆蓋面積統(tǒng)計(jì)圖Fig.3 Variation of the coverage area of the Green Tides from 2007 to 2015

    比較各年綠潮覆蓋面積極大值及各年覆蓋面積整體變化趨勢(shì)發(fā)現(xiàn):2007年暴發(fā)首年綠潮覆蓋面積最低;2008年青島海域暴發(fā)高密度、大面積綠潮,覆蓋面積極大值達(dá)3 110 km2,年平均覆蓋面積達(dá)1 266 km2,對(duì)沿海社會(huì)經(jīng)濟(jì)造成嚴(yán)重影響;2009年同樣迎來黃海綠潮大暴發(fā),其最大覆蓋面積達(dá)4 075 km2,超過2008年,但因綠潮漂移路徑和岸線平行東北向漂移,未對(duì)沿海海域造成大的影響;2010-2012年,黃海綠潮發(fā)生規(guī)模處于同一水平的低值,各年的峰值點(diǎn)均未超過1 800 km2,而年平均覆蓋面積在350~450 km2之間;2013年6月中旬,黃海綠潮最大覆蓋面積達(dá)3 618 km2,年平均覆蓋面積達(dá)1 250 km2,與2008年和2009年綠潮覆蓋規(guī)模相當(dāng);2014年綠潮覆蓋規(guī)模和2013年相比偏小,但影響范圍并未減?。?015綠潮覆蓋面積最大值達(dá)5 629 km2,年平均覆蓋面積達(dá)2 700 km2,是其他年份的2~8倍,持續(xù)時(shí)間達(dá)98 d之久,覆蓋面積、持續(xù)時(shí)間都達(dá)到歷史最高水平。

    圖4給出了2008—2015年覆蓋面積達(dá)極大值時(shí)的黃海綠潮分布,發(fā)現(xiàn)各年綠潮覆蓋面積極值出現(xiàn)時(shí)間、影響范圍及面積大小都存在差異。黃海綠潮覆蓋面積極大值在每年6月中旬至7月初出現(xiàn),覆蓋區(qū)域海表平均溫度在18~23℃,符合綠潮生長的最適溫度范圍。因各年綠潮到達(dá)峰值點(diǎn)時(shí)間和黃海海表動(dòng)力的差異,導(dǎo)致每年綠潮最大覆蓋面積的地理分布存在明顯變化,但綠潮主體均未越過124.2°E以東及江蘇大豐縣以南,部分年份綠潮峰值點(diǎn)的時(shí)間集中在6月中旬,綠潮主體還未漂移至當(dāng)年最北端,故未對(duì)沿岸產(chǎn)生影響。2008年、2013年和2014年綠潮達(dá)最大覆蓋面積時(shí)分布較為密集,綠潮主體分布在山東乳山以西、江蘇濱海以北海域,2008年和2013年綠潮沖襲上岸對(duì)日照—青島—海陽等山東沿海城市造成嚴(yán)重影響,而2014年在綠潮覆蓋面積達(dá)極大值時(shí)并未對(duì)沿海造成很大影響;2009年綠潮達(dá)最大覆蓋面積時(shí)主體偏東,覆蓋區(qū)域在青島以東、濱海以北海域,因7月初至8月初黃海北部以西-西南風(fēng)為主,導(dǎo)致綠潮與岸線平行東北向漂至124.2°E,并未對(duì)山東沿海造成影響;2010年和2015年綠潮達(dá)最大覆蓋面積時(shí)影響面積較其他年份偏大,綠潮南邊界在江蘇大豐沿海,東邊界在山東高角,2010年綠潮覆蓋面積峰值點(diǎn)密度較低,在連云港—日照—青島等沿海發(fā)現(xiàn)斑塊狀綠潮,造成的影響較2008年和2013小很多,2015年覆蓋面積在達(dá)到極大值后繼續(xù)北移,在山東沿岸大面積登陸;2011年和2012年最大覆蓋面積相對(duì)其他年份小,分布在連云港以北、海州灣外海海域,2011年綠潮主體北端在青島以南70 km海域,而2012年綠潮北端已在海陽沿岸登陸。

    圖4 2008—2015年黃海綠潮最大覆蓋面積分布Fig.4 Distribution of the maximum coverage area of the Green Tide in the Yellow Sea from 2008 to 2015

    圖5 2007—2015年黃海綠潮覆蓋面積月平均年際變化圖Fig.5 Variation of the monthly mean Green Tide coverage area in the Yellow Sea from 2007 to 2015

    因海上多云天氣的限制,每年獲取黃海綠潮清晰的衛(wèi)星圖像數(shù)量不同,且選用日期無固定規(guī)律,為進(jìn)行比較研究,引入覆蓋面積的月平均值。圖5給出了2007—2015年5月至8月黃海綠潮覆蓋面積月平均年際變化。發(fā)現(xiàn),因2009和2010年綠潮均在5月底6月初出現(xiàn),發(fā)生時(shí)間較晚,故這兩年綠潮的覆蓋面積是7月最高,6月次之;而其他年份5月中旬首次出現(xiàn),6月經(jīng)歷快速增殖,在30~40 d時(shí)間內(nèi)達(dá)到面積峰值,因此覆蓋面積6月最高,7月次之,到8月基本消失。

    3.3 黃海綠潮漂移路徑年際變化特征

    本文將綠潮分布區(qū)域的幾何中心確定為綠潮主體位置。2008—2015年黃海綠潮主體以偏北向漂移為主,且絕大部分年份會(huì)沖襲上岸,對(duì)山東和江蘇沿海城市造成不同程度的影響。但綠潮漂移路徑在具體方向和時(shí)間上存在明顯的年際變化(圖6),主要可以分為以下3種情況。

    圖6 2008-2015年黃海綠潮主體漂移路徑年際變化Fig.6 Trajectories of the Green Tide in the Yellow Sea from 2007 to 2015

    2008年和2011年綠潮主體先北向漂移越過34°30′N后,西北向垂直于岸線漂移,主要對(duì)連云港、日照和青島造成較大影響。2008年自5月中旬發(fā)現(xiàn)綠潮開始北向漂移,至5月底開始向西北方向垂直于岸線漂移,6月25日后,綠潮大規(guī)模北上入侵青島海岸,岸灘上堆積的綠潮厚度平均0.3~0.5 m[9],當(dāng)?shù)卣塾?jì)清運(yùn)綠潮超過1 000 000 t;2011年綠潮自6月初出現(xiàn)后先北向漂移,至6月中旬開始西北向漂移,7月中旬零星綠潮上岸,但相比2008年,綠潮主體并未大規(guī)模登陸,隨后綠潮主體沿岸線西南向漂移,在連云港至青島沿岸均有少許登陸。

    2009年和2012年綠潮主體較其他年份而言整體偏東,綠潮在南黃海輻射沙脊處出現(xiàn)后先西北向漂移越過34°30′N ,隨后東北向平行于岸線漂移,故只有少許綠潮上岸。2009年綠潮自5月末發(fā)現(xiàn)后先西北向漂移,6月下旬漂出輻射沙脊后開始平行于岸線東北向漂移,少許綠潮在海陽-乳山登陸,這與高松等[4]研究一致,認(rèn)為造成2009年和2008年綠潮漂移路徑不同的主要原因是綠潮發(fā)生海域的季風(fēng)差異及其引起的海表面環(huán)流差異共同作用的結(jié)果;2012年綠潮漂移路徑與2009年類似,在5月初出現(xiàn)后北向漂移,5月中旬漂出輻射沙脊后向東北向漂移,7月上旬也僅有少量綠潮在青島-海陽沿岸登陸,但其影響規(guī)模要比2009年小很多。

    2010、2013、2014和2015年綠潮主體西北向漂移至近岸后東北向沿岸漂移,對(duì)江蘇和山東沿海造成大面積影響。2010年6月初發(fā)現(xiàn)綠潮后先北向漂移,6月中下旬西北向漂移,至6月底在日照-青島登陸,隨后綠潮沿海岸線擴(kuò)展蔓延,南至連云港沿岸,北至海陽沿岸均有綠潮登陸;2013年先西北向漂移,6月下旬綠潮主體在連云港-日照-青島-海陽沿海登陸,隨后綠潮主體沿岸線東北向漂移,6月底對(duì)青島-乳山沿岸造成較大影響,7月下旬綠潮東邊界漂移至124°E;2014年自5月中旬發(fā)現(xiàn)綠潮開始呈擴(kuò)散式生長,并北向漂移,至6月中旬綠潮主體南端擴(kuò)散生長至大豐沿岸,大面積影響鹽城沿海,北端擴(kuò)散至海陽近海,隨后綠潮主體東北向沿岸漂移至榮成近岸,部分綠潮西北向漂移并在連云港大面積登陸;2015年先西北向漂移,5月底綠潮在大豐-濱海沿岸大面積登陸,隨后綠潮北向漂移至近岸后沿岸東北向漂移,6月中旬綠潮主體北端漂移至乳山近海,7月初綠潮已在日照-青島-榮成大面積登陸,影響范圍幾乎遍布南黃海。

    究其影響機(jī)制,國內(nèi)學(xué)者主要認(rèn)為綠潮所在海域的海面風(fēng)與表層環(huán)流的年際變化是綠潮漂移路徑年際變化的主要因素之一。Bao等[19]通過在南黃海輻射沙脊群釋放漂流瓶和Argos浮標(biāo)模擬綠潮的漂移路徑,認(rèn)為綠潮在偏南風(fēng)的作用下才能擺脫輻射沙脊地形限制,當(dāng)綠潮越過34.5°N后,漂移路徑變得復(fù)雜,但主體以偏北向漂移為主,這與本文的研究基本一致。

    4 總結(jié)

    黃海綠潮分布變化規(guī)律對(duì)于理解這一新的生態(tài)災(zāi)害過程,評(píng)估其經(jīng)濟(jì)和社會(huì)損失,開展預(yù)測、預(yù)警和防控有重要的借鑒意義。本文通過分析近9年黃海綠潮MODIS衛(wèi)星數(shù)據(jù)資料,對(duì)衛(wèi)星首次發(fā)現(xiàn)綠潮的時(shí)間、持續(xù)和消亡時(shí)間、覆蓋面積、分布以及綠潮的漂移路徑等特征進(jìn)行分析,發(fā)現(xiàn):

    (1)2007年為黃海綠潮暴發(fā)首年,至2015年已連續(xù)9年在南黃海暴發(fā)綠潮災(zāi)害。

    (2)不同年份間,黃海綠潮具有相似的發(fā)生發(fā)展及消亡過程。綠潮整個(gè)生消過程持續(xù)2~3個(gè)月,5月出現(xiàn),6月快速增殖,7月下旬快速消退,到8月中旬完全消失。

    (3)黃海綠潮覆蓋面積呈現(xiàn)波動(dòng)變化。經(jīng)歷2008年和2009年綠潮大暴發(fā)后,2010—2012年綠潮暴發(fā)規(guī)模異常低值,2013—2015年綠潮規(guī)模呈增長趨勢(shì),2015年綠潮暴發(fā)規(guī)模為歷年之最。

    (4)綠潮漂移路徑存在明顯的年際差異,各年綠潮北向漂移超出34.5°N后,漂移路徑變得復(fù)雜,大體可分為西北向垂直于岸線漂移、東北向平行于岸線漂移及繼續(xù)北向漂移3類。所有年份綠潮影響范圍均限于南黃海內(nèi),東側(cè)邊界最遠(yuǎn)未越過124.2°E。

    本文對(duì)于黃海綠潮覆蓋面積和漂移路徑年際變化特征方面得到一些認(rèn)識(shí),可為黃海綠潮年際變化研究提供一定的基礎(chǔ)資料。但在水文特征、風(fēng)場和營養(yǎng)鹽等因素等對(duì)綠潮生物量、分布及漂移路徑的影響、綠潮發(fā)生機(jī)理等方面還有待進(jìn)一步深入研究。

    [1] Charlier R H,Morand P,F(xiàn)inkl C W,et al. Green tides on the Brittany coasts[C]//IEEE US/EU Baltic International Symposium.Klaipeda,Lithuania:IEEE,2006.

    [2] Yabe T,Ishii Y,Amano Y,et al. Green tide formed by free-floatingUlvaspp. at Yatsu tidal flat,Japan[J]. Limnology,2009,10(3):239-245.

    [3] 唐啟升,張曉雯,葉乃好,等.綠潮研究現(xiàn)狀與問題[J]. 中國科學(xué)基金,2010(1):5-9.

    Tang Qisheng,Zhang Xiaowen,Ye Naihao,et al. Review on the research progress on marine green tide[J]. Bulletin of National Natural Science Foundation of China,2010(1):5-9.

    [4]高松,黃娟,白濤,等.2008年與2009年黃海綠潮漂移路徑分析[J].海洋科學(xué),2014,38(2):86-90.

    Gao Song,Huang Juan,Bai Tao,et al.Analysis on drifting path of Green Tides in the Yellow Sea in 2008 and 2009[J].Marine Sciences,2014,38(2):86-90.

    [5] 李大秋,賀雙顏,楊倩,等.青島海域滸苔來源與外海分布特征研究[J].環(huán)境保護(hù),2008(16):45-46.

    Li Daqiu,He Shuangyan,Yang Qian,et al.The origin and distribution of theEnteromorphaproliferain Qingdao coastal waters[J]. Environmental Protection,2008(16):45-46.

    [6] 李曰嵩,肖文軍,楊紅,等.2012年黃海綠潮藻早期發(fā)生和聚集動(dòng)力學(xué)成因分析[J].海洋環(huán)境科學(xué),2015,34(2):268-273.

    Li Yuesong,Xiao Wenjun,Yang Hong,et al. Dynamic analyses of early development and gather of green macroalgae in 2012[J]. Marine Environmental Science,2015,34(2):268-273.

    [7] 張?zhí)K平,劉應(yīng)辰,張廣泉,等.基于遙感資料的2008年黃海綠潮滸苔水文氣象條件分析[J]. 中國海洋大學(xué)學(xué)報(bào),2009,39(5):870-876.

    Zhang Suping,Liu Yingchen,Zhang Guangquan,et al. Analysis on the hydro-moteorological conditions from remote sensing data for the 2008 algal blooming in the Yellow Sea[J]. Periodical of Ocean University of China,2009,39(5):870-876.

    [8] 衣立,張?zhí)K平,殷玉齊.2009年黃海綠潮滸苔爆發(fā)與漂移的水文氣象環(huán)境[J].中國海洋大學(xué)學(xué)報(bào),2010,40(10):15-23.

    Yi Li,Zhang Suping,Yin Yuqi. Influence of environmental hydro-meteorological conditions toEnteromorphaproliferablooms in Yellow Sea,2009[J]. Periodical of Ocean University of China,2010,40(10):15-23.

    [9]喬方利,王關(guān)鎖,呂新剛,等.2008與2010年黃海滸苔漂移輸運(yùn)特征對(duì)比[J].科學(xué)通報(bào),2011,56(18):1470-1476.

    Qiao Fangli,Wang Guansuo,Lü Xin’gang,et al. Drift characteristics of green macroalgae in the Yellow Sea in 2008 and 2010[J]. Chinese Science Bulletin,2011,56(18):1470-1476.

    [10] Keesing J K,Liu Dongyan,F(xiàn)earns P,et al. Inter-and intra-annual patterns ofUlvaproliferagreen tides in the Yellow Sea during 2007-2009,their origin and relationship to the expansion of coastal seaweed aquaculture in China[J]. Marine Pollution Bulletin,2011,62(6):1169-1182.

    [11] 黃娟,吳玲娟,高松,等.黃海綠潮分布年際變化分析[J]. 激光生物學(xué)報(bào),2014,23(6):572-578.

    Huang Juan,Wu Lingjuan,Gao Song,et al. Analysis on the interannual distribution variation of Green Tide in Yellow Sea[J]. Acta Laser Biology Sinica,2014,23(6):572-578.

    [12] Jiang Peng,Wang Jinfeng,Cui Yulin,et al. Molecular phylogenetic analysis of attached Ulvaceae species and free-floatingEnteromorphafrom Qingdao coasts in 2007[J]. Chinese Journal of Oceanology and Limnology,2008,26(3):276-279.

    [13] 劉峰,逄少軍.黃海滸苔綠潮及其溯源研究進(jìn)展[J].海洋科學(xué)進(jìn)展,2012,30(3):441-449.

    Liu Feng,Pang Shaojun. Research advances on Green Tides in the Yellow Sea[J]. Advances in Marine Science,2012,30(3):441-449.

    [14] 吳孟泉,郭浩,張安定,等.2008年-2012年山東半島海域滸苔時(shí)空分布特征研究[J].光譜學(xué)與光譜分析,2014,34(5):1312-1318.

    Wu Mengquan,Guo Hao,Zhang Anding,et al.Research on the characteristics ofUlvaProliferain Shandong Peninsula during 2008-2012 based on MODIS data[J]. Spectroscopy and Spectral Analysis,2014,34(5):1312-1318.

    [15] 王玉麗,馬震.應(yīng)用ENVI軟件目視解譯TM影像土地利用分類[J].現(xiàn)代測繪,2011,34(1):11-13.

    Wang Yuli,Ma Zhen. Visual interpretation TM image land use classification by applied the software of ENVI[J]. Modern Surveying and Mapping,2011,34(1):11-13.

    [16] 蔡曉晴,崔廷偉,鄭榮兒,等.靜止海洋水色衛(wèi)星(GOCI)綠潮探測算法對(duì)比研究[J]. 遙感信息,2014,29(5):44-50.

    Cai Xiaoqing,Cui Tingwei,Zheng Rong’er,et al. Comparison of algorithms for green macro-algae bloom detection based on Geostationary Ocean Color Imager[J]. Remote Sensing Information,2014,29(5):44-50.

    [17] 李德萍,楊育強(qiáng),董海鷹,等. 2008年青島海域滸苔大爆發(fā)天氣特征及成因分析[J]. 中國海洋大學(xué)學(xué)報(bào),2009,39(6):1165-1170.

    Li Deping,Yang Yuqiang,Dong Haiying,et al. Cause analysis and synoptic characteristics for outbreaks of enteromorpha around Qingdao in 2008[J]. Periodical of Ocean University of China,2009,39(6):1165-1170.

    [18] 吳洪喜,徐愛光,吳美寧.滸苔實(shí)驗(yàn)生態(tài)的初步研究[J]. 浙江海洋學(xué)院學(xué)報(bào)(自然科學(xué)版),2000 ,19(3):230-234.

    Wu Hongxi,Xu Aiguang,Wu Meining. Preliminary study on experimental ecology ofEnteromorphaprolifera(Miill.)[J]. Journal of Zhejiang Ocean University (Natural Science),2000,19(3):230-234.

    [19] Bao Min,Guan Weibing,Yang Yang,et al. Drifting trajectories of green algae in the western Yellow Sea during the spring and summer of 2012[J]. Estuarine,Coastal & Shelf Science,2015,163:9-16.

    The interannual variation of Green Tide in the Yellow Sea

    Guo Wei1,Zhao Liang1,Li Xiumei1

    (1.CollegeofMarineandEnvironmentalSciences,TianjinUniversityofScienceandTechnology,Tianjin300457,China)

    From 2008 to 2015,large-scale floating Green Tides have consecutively occurred in the Yellow Sea of China,which have great impacts on the marine environment and economic. Based on the EOS/MODIS remote sensing ocean color data,theNDVIalgorithm was adopted to inverse the distribution of Green Tides during May to August. The Green Tides were firstly found from MODIS images in 2007 with the coverage area less than 50 km2in the Yellow Sea. Seasonally,the Green Tides first appeared during the middle of May to early June,then its coverage grew rapidly and continuously and reached a peak in 30-40 days. In July,the Green Tides depleted gradually and disappeared in August. The maximum coverage areas of the Green Tides were 3 110 km2and 4 075 km2in 2008 and 2009 respectively. The Green Tides coverages were in low phase from 2010 to 2012 and the peak areas in these past three years were less than 1 800 km2. Since 2013,the scales of Green Tides increased rapidly,and the peak coverage reached 5 629 km2in 2015. (The Green Tides covered the regions ranged from 124°E in the east,Lianyungang of Jiangsu Province in the west,Rudong of Jiangsu Province in the south and 37°N in the north). The trajectories of the Green Tides could be divided into three categories:(a) moved northward to 34°30′N,then turned to northwest toward the coast in 2008,2010 and 2011; (b) moved northward to 35°N,then drifted northeast parallel to the coast in 2009 and 2012; (c) moved northward to the coast of Shandong province,then turned along the coast in 2013,2014 and 2015.

    Green Tide; coverage area; trajectories; interannual variation; Yellow Sea

    2016-04-18;

    2016-05-16。

    中國科學(xué)院戰(zhàn)略先導(dǎo)專項(xiàng)(XDA11020305);國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFA0601301);國家自然科學(xué)基金(41276016,41376006,41376112)。

    郭偉(1990—),男,江蘇省南通市人,從事淺海動(dòng)力學(xué)研究。E-mail:guowei72@tust.edu.cn

    *通信作者:趙亮(1975—),男,教授,從事海洋生態(tài)系統(tǒng)動(dòng)力學(xué)研究。E-mail:zhaoliang@ouc.edu.cn

    10.3969/j.issn.0253-4193.2016.12.004

    X55

    A

    0253-4193(2016)12-0036-10

    郭偉,趙亮,李秀梅. 黃海綠潮分布年際變化特征分析[J].海洋學(xué)報(bào),2016,38(12):36—45,

    Guo Wei,Zhao Liang,Li Xiumei. The interannual variation of Green Tide in the Yellow Sea[J]. Haiyang Xuebao,2016,38(12):36—45,doi:10.3969/j.issn.0253-4193.2016.12.004

    猜你喜歡
    綠潮覆蓋面積年際
    HY-1衛(wèi)星CZI影像在黃海綠潮監(jiān)測的應(yīng)用
    綠潮監(jiān)測行業(yè)標(biāo)準(zhǔn)通過專家審查
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    金星上全是金子嗎
    覆蓋面積和網(wǎng)圍結(jié)構(gòu)對(duì)水面蒸發(fā)抑制率的影響
    人民長江(2019年3期)2019-10-20 04:50:19
    基于GOCI數(shù)據(jù)的2015年黃海海域綠潮監(jiān)測及時(shí)空變化分析
    測繪工程(2018年1期)2018-01-09 06:03:31
    基于MODIS數(shù)據(jù)的2016年黃海綠潮災(zāi)害動(dòng)態(tài)監(jiān)測研究
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    基于M-K法對(duì)圖們江干流含沙量年際變化的分析
    青春草国产在线视频| 五月玫瑰六月丁香| 大又大粗又爽又黄少妇毛片口| 天堂av国产一区二区熟女人妻| 国产亚洲一区二区精品| 校园人妻丝袜中文字幕| 亚洲av男天堂| 只有这里有精品99| 国产高潮美女av| 网址你懂的国产日韩在线| 成人高潮视频无遮挡免费网站| 精品久久久久久成人av| 午夜爱爱视频在线播放| 91久久精品国产一区二区成人| 国产黄a三级三级三级人| 国产精品精品国产色婷婷| 青青草视频在线视频观看| 嫩草影院新地址| 91在线精品国自产拍蜜月| 少妇裸体淫交视频免费看高清| 欧美bdsm另类| 99re6热这里在线精品视频| 极品少妇高潮喷水抽搐| 高清毛片免费看| 久久精品国产亚洲网站| 男人和女人高潮做爰伦理| 一级毛片我不卡| 91在线精品国自产拍蜜月| 色哟哟·www| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜免费激情av| 久久97久久精品| 国产精品爽爽va在线观看网站| 日韩欧美三级三区| 日韩精品有码人妻一区| 日韩欧美精品v在线| 久久热精品热| 国产成人精品久久久久久| 老司机影院毛片| www.av在线官网国产| 国产伦在线观看视频一区| av专区在线播放| 成人一区二区视频在线观看| 国产成人a区在线观看| 精品久久久久久成人av| 如何舔出高潮| 日本黄色片子视频| 人妻少妇偷人精品九色| 免费观看的影片在线观看| 亚洲婷婷狠狠爱综合网| 色播亚洲综合网| 欧美日韩综合久久久久久| 亚洲不卡免费看| 国产黄片视频在线免费观看| 国产黄片视频在线免费观看| 免费观看无遮挡的男女| 久久久精品94久久精品| 五月玫瑰六月丁香| 欧美精品国产亚洲| 亚洲欧洲日产国产| 日韩精品有码人妻一区| 黄片wwwwww| 日本-黄色视频高清免费观看| 波野结衣二区三区在线| 欧美激情久久久久久爽电影| 国产午夜精品一二区理论片| 日日啪夜夜爽| 午夜精品一区二区三区免费看| 国产又色又爽无遮挡免| 少妇高潮的动态图| 日韩欧美 国产精品| 国产av码专区亚洲av| 国产午夜精品论理片| 亚洲av中文字字幕乱码综合| 小蜜桃在线观看免费完整版高清| 男女那种视频在线观看| videos熟女内射| 国产乱人偷精品视频| 久久精品熟女亚洲av麻豆精品 | 日韩人妻高清精品专区| 国产精品av视频在线免费观看| 97超视频在线观看视频| av在线播放精品| 成年人午夜在线观看视频 | 人妻制服诱惑在线中文字幕| 国产免费一级a男人的天堂| 亚洲国产精品sss在线观看| 777米奇影视久久| 一区二区三区四区激情视频| 亚洲国产最新在线播放| 亚洲一区高清亚洲精品| 大又大粗又爽又黄少妇毛片口| 国产男女超爽视频在线观看| 2022亚洲国产成人精品| 久久久久久久午夜电影| 国产伦在线观看视频一区| 美女脱内裤让男人舔精品视频| 亚洲乱码一区二区免费版| 国产精品久久久久久精品电影| 国产亚洲最大av| av卡一久久| 97精品久久久久久久久久精品| 中文字幕免费在线视频6| 成人综合一区亚洲| 久久久久久久久久成人| 欧美潮喷喷水| 天堂av国产一区二区熟女人妻| av在线亚洲专区| 欧美+日韩+精品| 三级国产精品欧美在线观看| 国产av在哪里看| 免费看a级黄色片| 婷婷六月久久综合丁香| 久久精品综合一区二区三区| 亚洲欧洲日产国产| 国产成人福利小说| 在线观看一区二区三区| 久久久成人免费电影| 国产一区二区在线观看日韩| 美女主播在线视频| 男人舔女人下体高潮全视频| 日本与韩国留学比较| 日韩一本色道免费dvd| 欧美精品一区二区大全| 永久免费av网站大全| 人体艺术视频欧美日本| 欧美不卡视频在线免费观看| 天天一区二区日本电影三级| 亚洲人成网站在线观看播放| 小蜜桃在线观看免费完整版高清| 黄色欧美视频在线观看| 亚洲性久久影院| 久久久久久久久久久免费av| 国产精品人妻久久久久久| 午夜免费观看性视频| 成人一区二区视频在线观看| 国产乱人偷精品视频| 精品久久久噜噜| 国产成人一区二区在线| 国产有黄有色有爽视频| 日韩精品青青久久久久久| 黄片无遮挡物在线观看| 最新中文字幕久久久久| 亚洲内射少妇av| 偷拍熟女少妇极品色| 美女内射精品一级片tv| 简卡轻食公司| 极品教师在线视频| av黄色大香蕉| 久久久欧美国产精品| 内射极品少妇av片p| 少妇熟女aⅴ在线视频| 国产伦理片在线播放av一区| 男的添女的下面高潮视频| 国产精品综合久久久久久久免费| 国产乱人视频| 一级毛片 在线播放| 久久99热这里只频精品6学生| 久久久色成人| 麻豆av噜噜一区二区三区| 国产久久久一区二区三区| 国产精品麻豆人妻色哟哟久久 | 亚洲电影在线观看av| 国产一级毛片七仙女欲春2| 日韩欧美精品v在线| 亚洲精品成人久久久久久| 少妇被粗大猛烈的视频| 少妇的逼好多水| 欧美日韩视频高清一区二区三区二| 成人无遮挡网站| 男插女下体视频免费在线播放| 国产黄色视频一区二区在线观看| a级毛色黄片| 中文字幕亚洲精品专区| 国产欧美日韩精品一区二区| 高清午夜精品一区二区三区| 精品久久久久久久久亚洲| 久久久精品欧美日韩精品| 亚洲国产精品成人久久小说| 欧美极品一区二区三区四区| 黄色欧美视频在线观看| 九草在线视频观看| 久久热精品热| 亚洲国产精品专区欧美| 午夜亚洲福利在线播放| 久热久热在线精品观看| 日本一本二区三区精品| 色网站视频免费| 国产精品av视频在线免费观看| 免费观看无遮挡的男女| 三级毛片av免费| 狠狠精品人妻久久久久久综合| 亚洲精品乱码久久久v下载方式| 高清毛片免费看| 日韩视频在线欧美| 国产 一区 欧美 日韩| 精品久久久久久久人妻蜜臀av| 尤物成人国产欧美一区二区三区| 久久久久久久国产电影| 久久精品国产亚洲av天美| 色播亚洲综合网| 18+在线观看网站| 99热6这里只有精品| 老女人水多毛片| 国产亚洲av片在线观看秒播厂 | 人妻制服诱惑在线中文字幕| 内射极品少妇av片p| 国产高清不卡午夜福利| 1000部很黄的大片| 男的添女的下面高潮视频| 亚洲国产精品专区欧美| 我的女老师完整版在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲经典国产精华液单| 国产精品.久久久| 亚洲欧美日韩东京热| 亚洲第一区二区三区不卡| 1000部很黄的大片| 91午夜精品亚洲一区二区三区| 国产极品天堂在线| av免费观看日本| 亚州av有码| 女人十人毛片免费观看3o分钟| 国产精品蜜桃在线观看| 亚洲无线观看免费| 大话2 男鬼变身卡| 久久99热这里只有精品18| 国产精品精品国产色婷婷| 校园人妻丝袜中文字幕| 欧美成人a在线观看| 男人狂女人下面高潮的视频| 熟妇人妻久久中文字幕3abv| 久久精品久久久久久噜噜老黄| 午夜福利在线观看吧| 亚洲不卡免费看| www.av在线官网国产| 国产白丝娇喘喷水9色精品| 国产黄频视频在线观看| 亚洲av福利一区| 国产淫片久久久久久久久| 久久这里有精品视频免费| 久久久亚洲精品成人影院| 久久久久久国产a免费观看| 久久久久精品性色| 日韩一区二区视频免费看| 22中文网久久字幕| 26uuu在线亚洲综合色| 我的女老师完整版在线观看| 久久久色成人| 青春草视频在线免费观看| 亚洲激情五月婷婷啪啪| 日韩强制内射视频| 青春草国产在线视频| 亚洲精华国产精华液的使用体验| 十八禁网站网址无遮挡 | 色综合色国产| 亚洲国产精品专区欧美| 午夜福利视频1000在线观看| 亚洲国产精品国产精品| 精品一区二区三区人妻视频| 在线免费观看的www视频| 成人av在线播放网站| 欧美日本视频| 26uuu在线亚洲综合色| 免费观看的影片在线观看| 欧美一级a爱片免费观看看| 日韩欧美 国产精品| 少妇的逼好多水| 亚洲精品国产av成人精品| 免费观看精品视频网站| 国内少妇人妻偷人精品xxx网站| 日本色播在线视频| 午夜福利视频精品| 久久精品国产亚洲av天美| 欧美区成人在线视频| xxx大片免费视频| 五月天丁香电影| 色5月婷婷丁香| 免费播放大片免费观看视频在线观看| 国产免费视频播放在线视频 | 欧美日韩一区二区视频在线观看视频在线 | 51国产日韩欧美| 成人鲁丝片一二三区免费| 成人午夜精彩视频在线观看| 久久久久久九九精品二区国产| 婷婷六月久久综合丁香| 欧美性猛交╳xxx乱大交人| 男插女下体视频免费在线播放| 国产91av在线免费观看| 免费大片黄手机在线观看| 亚洲av中文av极速乱| 日韩av免费高清视频| 99久久精品国产国产毛片| 亚洲久久久久久中文字幕| 亚洲精华国产精华液的使用体验| 国产精品综合久久久久久久免费| 日韩欧美三级三区| 在线观看免费高清a一片| 免费黄网站久久成人精品| 欧美日韩一区二区视频在线观看视频在线 | 赤兔流量卡办理| 91午夜精品亚洲一区二区三区| 一级毛片 在线播放| 国内精品宾馆在线| 看黄色毛片网站| 只有这里有精品99| 色播亚洲综合网| 午夜免费观看性视频| 男女边摸边吃奶| 2018国产大陆天天弄谢| 亚洲精品第二区| 日韩国内少妇激情av| 亚洲av男天堂| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 色吧在线观看| 在线免费十八禁| 午夜免费男女啪啪视频观看| 少妇人妻一区二区三区视频| 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩无卡精品| 美女内射精品一级片tv| 伊人久久国产一区二区| 91久久精品电影网| 日本av手机在线免费观看| av线在线观看网站| 成人毛片60女人毛片免费| 超碰av人人做人人爽久久| 亚洲av在线观看美女高潮| 国产高潮美女av| 久久久精品欧美日韩精品| 久久精品国产亚洲av天美| 欧美日韩一区二区视频在线观看视频在线 | 国产中年淑女户外野战色| 国产综合懂色| 国产极品天堂在线| 一本久久精品| 午夜免费观看性视频| 日产精品乱码卡一卡2卡三| 日本免费a在线| 日韩一区二区三区影片| 久久这里有精品视频免费| 纵有疾风起免费观看全集完整版 | 色播亚洲综合网| 黄色欧美视频在线观看| 午夜福利在线观看免费完整高清在| 亚洲精品成人久久久久久| 好男人视频免费观看在线| 国产精品一二三区在线看| 国产午夜精品一二区理论片| 亚洲真实伦在线观看| 91久久精品国产一区二区三区| 最近的中文字幕免费完整| 国产成人精品一,二区| 免费高清在线观看视频在线观看| 午夜福利在线观看吧| 国产亚洲av片在线观看秒播厂 | 亚洲精品成人久久久久久| 少妇裸体淫交视频免费看高清| 建设人人有责人人尽责人人享有的 | 国产精品一区二区三区四区免费观看| 亚洲国产精品专区欧美| 亚洲欧美清纯卡通| 精品一区二区三卡| 亚洲精品亚洲一区二区| 精品一区在线观看国产| 日本黄色片子视频| 人妻夜夜爽99麻豆av| 青青草视频在线视频观看| 久久久a久久爽久久v久久| 国产黄色免费在线视频| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 国产不卡一卡二| 人妻少妇偷人精品九色| 精品国产一区二区三区久久久樱花 | 久久韩国三级中文字幕| 欧美潮喷喷水| 少妇被粗大猛烈的视频| 亚洲精品影视一区二区三区av| 中文字幕免费在线视频6| 三级经典国产精品| 男女下面进入的视频免费午夜| 男女视频在线观看网站免费| 日本一本二区三区精品| 午夜福利高清视频| 淫秽高清视频在线观看| 韩国高清视频一区二区三区| 黄色欧美视频在线观看| 身体一侧抽搐| 日本爱情动作片www.在线观看| 丝瓜视频免费看黄片| 中文天堂在线官网| 青青草视频在线视频观看| 精品久久久久久成人av| 日本-黄色视频高清免费观看| 国产一区二区亚洲精品在线观看| 如何舔出高潮| 边亲边吃奶的免费视频| 一二三四中文在线观看免费高清| 国产精品久久久久久久久免| 白带黄色成豆腐渣| 热99在线观看视频| 亚洲精品国产av成人精品| 黄色欧美视频在线观看| 纵有疾风起免费观看全集完整版 | 一个人观看的视频www高清免费观看| 亚洲不卡免费看| 中文字幕av成人在线电影| 插阴视频在线观看视频| 99re6热这里在线精品视频| 亚洲欧洲日产国产| 午夜福利成人在线免费观看| 免费观看av网站的网址| 超碰av人人做人人爽久久| 观看免费一级毛片| 在线观看av片永久免费下载| 亚洲aⅴ乱码一区二区在线播放| 高清在线视频一区二区三区| 午夜日本视频在线| 日本av手机在线免费观看| videos熟女内射| 亚洲精品久久久久久婷婷小说| 一二三四中文在线观看免费高清| 亚洲av不卡在线观看| av.在线天堂| 校园人妻丝袜中文字幕| 九九爱精品视频在线观看| 夜夜看夜夜爽夜夜摸| 丰满乱子伦码专区| 国产精品熟女久久久久浪| 99久久精品国产国产毛片| 中文字幕av在线有码专区| 亚洲怡红院男人天堂| 欧美丝袜亚洲另类| 黄色一级大片看看| 久久久精品94久久精品| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 国产亚洲av嫩草精品影院| 精品久久久久久久人妻蜜臀av| 欧美日韩在线观看h| 男人和女人高潮做爰伦理| 国产又色又爽无遮挡免| 久久国产乱子免费精品| 欧美成人一区二区免费高清观看| 亚洲av成人精品一区久久| 国产黄色免费在线视频| 校园人妻丝袜中文字幕| 夜夜爽夜夜爽视频| 婷婷六月久久综合丁香| 免费黄色在线免费观看| .国产精品久久| 亚洲最大成人av| 老女人水多毛片| 亚洲熟女精品中文字幕| 18禁裸乳无遮挡免费网站照片| 天堂av国产一区二区熟女人妻| 国产极品天堂在线| 久久99热这里只频精品6学生| 精品亚洲乱码少妇综合久久| 国产免费福利视频在线观看| 亚洲av一区综合| av专区在线播放| 美女黄网站色视频| 99热这里只有精品一区| 少妇丰满av| 国产精品蜜桃在线观看| 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 麻豆国产97在线/欧美| 高清毛片免费看| 2021天堂中文幕一二区在线观| 五月天丁香电影| 欧美区成人在线视频| 嫩草影院新地址| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 国产精品伦人一区二区| ponron亚洲| 亚洲四区av| 亚洲色图av天堂| 少妇的逼好多水| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 青春草国产在线视频| 国内少妇人妻偷人精品xxx网站| 深爱激情五月婷婷| 免费观看精品视频网站| 美女黄网站色视频| 亚洲精品456在线播放app| 精品国产一区二区三区久久久樱花 | 国产爱豆传媒在线观看| 插阴视频在线观看视频| 看免费成人av毛片| 亚洲国产欧美人成| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 色综合站精品国产| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 午夜免费激情av| 欧美激情国产日韩精品一区| 床上黄色一级片| 久热久热在线精品观看| 身体一侧抽搐| 午夜免费激情av| 永久网站在线| 国产亚洲av嫩草精品影院| 国产精品.久久久| 国产免费视频播放在线视频 | 亚洲精品国产成人久久av| 少妇的逼水好多| 久久久国产一区二区| 亚洲在线观看片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 成人亚洲欧美一区二区av| 少妇裸体淫交视频免费看高清| 尾随美女入室| 国产精品精品国产色婷婷| 中文字幕av在线有码专区| 色吧在线观看| 成年av动漫网址| 欧美日韩在线观看h| 午夜激情福利司机影院| 日韩精品青青久久久久久| 国产黄色视频一区二区在线观看| 国产精品1区2区在线观看.| 一二三四中文在线观看免费高清| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 秋霞伦理黄片| 一级二级三级毛片免费看| 夫妻性生交免费视频一级片| 亚洲精品成人av观看孕妇| 久久99热这里只有精品18| 亚洲人成网站在线观看播放| 国产精品久久视频播放| 神马国产精品三级电影在线观看| 亚洲图色成人| av一本久久久久| 国产高清不卡午夜福利| a级毛色黄片| 国产精品国产三级国产av玫瑰| 国产伦一二天堂av在线观看| 最近手机中文字幕大全| 亚洲av国产av综合av卡| 亚洲精品国产av蜜桃| 一级a做视频免费观看| 国产白丝娇喘喷水9色精品| 国产麻豆成人av免费视频| 亚洲怡红院男人天堂| kizo精华| 国产亚洲av嫩草精品影院| 亚洲熟女精品中文字幕| 大又大粗又爽又黄少妇毛片口| 日日摸夜夜添夜夜爱| 久久久亚洲精品成人影院| 精品亚洲乱码少妇综合久久| av在线天堂中文字幕| 天堂中文最新版在线下载 | 久久久久久久大尺度免费视频| 床上黄色一级片| 人人妻人人澡人人爽人人夜夜 | 纵有疾风起免费观看全集完整版 | 久久99蜜桃精品久久| 一级爰片在线观看| 欧美激情在线99| 久久久亚洲精品成人影院| 啦啦啦啦在线视频资源| 亚洲精品成人久久久久久| 亚洲精品第二区| 91午夜精品亚洲一区二区三区| av国产久精品久网站免费入址| av天堂中文字幕网| 嫩草影院入口| 国产精品.久久久| 成人午夜高清在线视频| 91狼人影院| av在线天堂中文字幕| 中文字幕av在线有码专区| 亚洲欧洲日产国产| 成年女人在线观看亚洲视频 | 国产亚洲91精品色在线| 亚洲av在线观看美女高潮| 久久久久久伊人网av| 亚洲av不卡在线观看| 亚洲av男天堂| 亚洲人成网站高清观看| 国产黄a三级三级三级人| 欧美97在线视频| 国产av不卡久久| 可以在线观看毛片的网站| 噜噜噜噜噜久久久久久91| 亚洲,欧美,日韩| 国产黄片视频在线免费观看| 看十八女毛片水多多多| 高清欧美精品videossex| 91午夜精品亚洲一区二区三区| 精品一区在线观看国产| 国产一区二区三区综合在线观看 | 国产伦一二天堂av在线观看| 国产成人a∨麻豆精品| 天堂中文最新版在线下载 | 色网站视频免费| 性色avwww在线观看| 男女边摸边吃奶| 极品教师在线视频| 亚洲精品,欧美精品| 好男人在线观看高清免费视频| 久久久久久九九精品二区国产|