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

    結(jié)合衛(wèi)星遙感技術(shù)的太湖藍藻水華形成風場特征

    2016-04-16 07:12:23李亞春謝小萍朱小莉景元書江蘇省氣象局江蘇南京0008常州市氣象局江蘇常州0008南京信息工程大學(xué)應(yīng)用氣象學(xué)院江蘇南京0044
    中國環(huán)境科學(xué) 2016年2期
    關(guān)鍵詞:數(shù)值模擬

    李亞春,謝小萍,杭 鑫,朱小莉,黃 珊,景元書(.江蘇省氣象局,江蘇 南京 0008;.常州市氣象局,江蘇 常州 0008;.南京信息工程大學(xué)應(yīng)用氣象學(xué)院,江蘇 南京 0044)

    ?

    結(jié)合衛(wèi)星遙感技術(shù)的太湖藍藻水華形成風場特征

    李亞春1*,謝小萍1,杭 鑫1,朱小莉2,黃 珊1,景元書3(1.江蘇省氣象局,江蘇 南京 210008;2.常州市氣象局,江蘇 常州 210008;3.南京信息工程大學(xué)應(yīng)用氣象學(xué)院,江蘇 南京 210044)

    摘要:為進一步了解太湖藍藻水華形成和分布與近地面風場的關(guān)系,利用太湖湖面及周邊地區(qū)2003~2013年氣象與衛(wèi)星觀測數(shù)據(jù)分析、并應(yīng)用WRF3.5.1數(shù)值模型模擬,發(fā)現(xiàn)太湖藍藻水華主要出現(xiàn)在衛(wèi)星觀測時刻前6h平均風速為0.5~3.4m/s的區(qū)間,占比達94.7%;藍藻水華面積總體上隨風速增大而減小,大范圍藍藻水華主要出現(xiàn)在前6h平均風速≤2m/s的情形下,占比達89%;風向則主要影響藍藻水華在太湖的空間分布格局.結(jié)果表明局地風場對于太湖藍藻水華的形成、輸移和分布具有重要作用.

    關(guān)鍵詞:藍藻水華;近地面風場;衛(wèi)星遙感;數(shù)值模擬

    * 責任作者, 研究員級高級工程師, jsqxlyc@163.com

    氣候變化和水質(zhì)富營養(yǎng)化被認為是近幾十年來全球藍藻水華頻發(fā)的兩個重要原因[1-3],也是近年來太湖藍藻水華增多趨重的主要原因[4-5].藍藻水華的形成除了其本身的生理特點以外,氣象條件、營養(yǎng)鹽和其他環(huán)境因素也都起著關(guān)鍵作用,但任何一種單因子都難以驅(qū)動水華形成[6-8].當湖泊水質(zhì)富營養(yǎng)化程度難以得到控制的情況下,磷、氮濃度可能不再是藻類生長和水華形成的限制因子,其他環(huán)境因素可能取而代之,其中氣象條件顯得尤為重要[9].

    氣象和水文條件一般通過影響湖泊水體的分層、熱量傳遞交換以及光照、營養(yǎng)鹽的可利用性等,直接或間接地影響藻類個數(shù)、群落、分布和生命周期等[10].現(xiàn)有研究已證實,氣溫、風、光照、降水等氣象因子都會對太湖藍藻生長和水華形成產(chǎn)生重要影響[11-13].其中,風的影響相對復(fù)雜,一方面,其產(chǎn)生的風浪和湖流改變了藍藻的位置,使其積聚形成水華,另一方面則因風浪擾動會促使大量的營養(yǎng)鹽從底泥中釋放,導(dǎo)致湖泊水體重新分層與混合,增加了藻類可利用的營養(yǎng)鹽[14-16].相關(guān)研究[17-19]表明了風或風浪對太湖水體中膠體態(tài)營養(yǎng)鹽、浮游植物和藻類存在重要的影響, 而Kanoshina等[20]認為風引起的水流可能決定了藍藻水華的空間分布,Zhou等[21]則認為風的持續(xù)時間及其動力性是影響太湖藍藻水華時空間變化的關(guān)鍵因子.有研究表明,通常情況下微風有利于藍藻水華的生長和漂浮,孫小靜等[17]發(fā)現(xiàn)風速<4m/s的小風浪有利于藍藻生長或漂浮,王銘瑋等[22]則認為小于2.2m/s的低風速有利于淀山湖的藍藻水華形成.Wu等[23-24]利用MODIS衛(wèi)星圖像結(jié)合實地抽樣調(diào)查研究發(fā)現(xiàn),太湖藍藻水華面積與風速成負相關(guān),較小風速有利于藍藻水華形成,但短時間強風(30分鐘平均風速超過6m/s)引起的強烈混合作用則會增加大面積藍藻水華形成的機會,并且太湖地區(qū)四月和十月大風頻繁可能正是藍藻水華大面積發(fā)生的重要原因.王成林等[25-26]應(yīng)用數(shù)值模擬方法研究了適宜太湖藍藻水華形成的風場輻散特征及其形成機制,認為太湖區(qū)域存在特有的輻散場,藍藻水華在此輻散風場的驅(qū)動下會改變遷移方向.王文蘭等[27]利用WRF V2模式分析了近地面風場變化對太湖藍藻暴發(fā)的影響,認為藍藻水華面積及輸移方向?qū)孛骘L場的響應(yīng)相當迅速.由此可見,風對藍藻水華的形成、輸移和分布的影響很大,但到目前為止還沒有完全弄清其影響機制,而防控藍藻水華越來越需要能夠?qū)ζ浒l(fā)生發(fā)展和位置變化作出精準預(yù)測,風向、風速就應(yīng)當作為必須考慮的環(huán)境驅(qū)動因子[28].此外,現(xiàn)有的研究大多或基于實驗數(shù)據(jù),或基于短時間、間斷性的實測數(shù)據(jù),有些甚至僅僅是個例分析,而太湖是一個大型淺水湖泊,藍藻水華面積常常達數(shù)百平方公里,因此有些結(jié)果難以反映實際情況,而且不同的結(jié)果之間有時會存在明顯差異.

    為此,本文利用太湖湖面及周邊地區(qū)2003~ 2013年連續(xù)高頻的氣象觀測數(shù)據(jù)及同步衛(wèi)星監(jiān)測數(shù)據(jù),從氣象學(xué)角度深入分析太湖藍藻水華形成、輸移與分布的風場特征,明確藍藻水華、尤其是大面積藍藻水華形成的適宜風速參數(shù),并對太湖區(qū)域風場進行高分辨率數(shù)值模擬,分析近地面風場變化對藍藻水華的影響,以期為進一步研究藍藻水華形成機理和輸移過程提供理論依據(jù),為太湖藍藻水華的預(yù)測、預(yù)警和防控提供技術(shù)支持.

    1 資料與方法

    1.1 資料

    氣象資料來源于江蘇省氣象局,主要選用了太湖湖區(qū)及周邊區(qū)域的5個氣象基本站和23個區(qū)域自動站的觀測資料(圖1),基本站觀測資料包括宜興、無錫、蘇州、吳中和東山站2003~2013年逐時風向、風速觀測資料,區(qū)域自動站觀測資料為2009~2013年時間間隔為10min的風向、風速觀測數(shù)據(jù).

    圖1 氣象數(shù)據(jù)采樣點位置示意Fig.1 Location of the sampling meteorological site

    衛(wèi)星數(shù)據(jù)來源于國家衛(wèi)星氣象中心和江蘇省氣象局,選用2003~2013年具有較高時空分辨率的Aqua/Terra衛(wèi)星(EOS/MODIS)和 FY-3衛(wèi)星(MERSI)觀測的影像數(shù)據(jù).

    1.2 太湖藍藻的遙感解譯方法

    衛(wèi)星遙感在監(jiān)測大型湖泊水質(zhì)和大面積藍藻水華方面具有獨特優(yōu)勢.藍藻水華的遙感解譯方法很多,包括單波段、波段差值、波段比值、歸一化植被指數(shù)(NDVI)及藻類指數(shù)模型等[29-31].本研究選用目前業(yè)務(wù)上常用的NDVI方法,解譯得到2003~2013年累計805幅太湖藍藻水華面積≥1km2的衛(wèi)星遙感影像,解讀出藍藻水華頻次、面積、時間及位置信息等.東太湖和東部沿岸區(qū)水質(zhì)較好,很少發(fā)生藍藻聚集現(xiàn)象,在NDVI提取過藍藻水華過程中以掩膜技術(shù)進行處理.

    為考察大范圍藍藻水華形成與風的關(guān)系,定義單次藍藻水華面積超過太湖水面面積的20% 或468km2以上的為大范圍藍藻水華,得到大范圍藍藻水華衛(wèi)星遙感影像樣本62個,除2003年以外的其他年份均有發(fā)生.

    1.3 風的統(tǒng)計分析方法

    太湖藍藻水華的形成與輸移對風的響應(yīng)較為迅速,但同時也是前期氣象條件影響的累積,本研究統(tǒng)計了衛(wèi)星觀測時刻前1~6h的風觀測數(shù)據(jù),考察各時間段平均風速、最大風速、最小風速和風向與藍藻水華之間的關(guān)系.風的統(tǒng)計依據(jù)《地面氣象觀測規(guī)范》[32]規(guī)定的方法.

    1.4 風場的數(shù)值模擬方法

    為全面、詳細地了解近地層風場變化對藍藻水華的影響,利用數(shù)值模式WRF 3.5.1對太湖區(qū)域的風場進行高分辨率模擬,分析近地面風場變化對藍藻水華的強度、輸移及分布范圍的影響.模擬時間為2010年9月3日08:00~9月8日08:00,期間太湖藍藻水華經(jīng)歷了較明顯的生消過程.WRF 3.5.1數(shù)值模式采用三層嵌套結(jié)構(gòu),格距分別為30、10和3.3km,格點數(shù)分別為232×172、166×121和109×79.模擬方案的初始場和邊界氣象資料取自美國國家環(huán)境預(yù)報中心(NCEP)提供的空間分辨率為1°×1°,時間間隔為6h的FNL(Final Operational Global Analysis)全球分析資料,云微物理參數(shù)化方案采用Purdue Lin,長波輻射方案為RRTM,短波輻射方案為Dudhia,陸面過程為Noah陸面模式,積云參數(shù)化方案為Grell-Devenyi ensemble scheme,邊界層方案為Mellor-Yamada-Janjic.

    2 結(jié)果與討論

    2.1 風速對藍藻水華頻次的影響

    2003~2013年,衛(wèi)星解譯得到太湖出現(xiàn)藍藻水華共計805次,統(tǒng)計衛(wèi)星觀測時刻前6h的平均風速,得到觀測時刻前6h的平均風速分布情況(圖2).衛(wèi)星觀測到各次藍藻水華前6h的平均風速均介于0.2~4.8m/s之間,最小值0.2m/s,最大值4.8m/s,平均值1.8m/s.

    圖2 2003~2013年衛(wèi)星觀測時刻前6h平均風速分布Fig.2 Average wind speed in 6hours before satellites observation time during 2003~2013

    表1 衛(wèi)星觀測時刻前1~6h平均風速范圍及藍藻水華頻次總體分布Table 1 Average wind speed in 1~6h before satellites observation time and cyanobacteria bloom frequency

    為分析不同風速對藍藻水華頻次的影響,按照1.0m/s的風速間隔進行分段統(tǒng)計,考慮到在風速0.5m/s以下和3.5m/s以上藍藻水華頻次均較少,故將其分別設(shè)置為上下限.表1為衛(wèi)星監(jiān)測到藍藻水華前6h的最大、最小風速和不同時段不同風速區(qū)間的藍藻水華頻次分布情況.由表1可知,衛(wèi)星監(jiān)測到藍藻水華前6h的風速一般都超過0.5m/s,<0.5m/s的藍藻水華頻次占比均不超過2%,而風速≥3.5m/s的頻次占比也僅為3.4%~ 10.4%,大部分藍藻水華都出現(xiàn)在平均風速為0.5~3.4m/s的區(qū)間內(nèi),占比從1h最小為89.3%增至6h時最大為94.7%,其中1.5~2.4m/s區(qū)間占比最高,各時段占比均超過40%.表明太湖絕大多數(shù)藍藻水華都出現(xiàn)在0.5~3.4m/s的風速范圍內(nèi),在前6h平均風速小于0.5或大于3.5m/s的情況下,衛(wèi)星均較少觀測到藍藻水華的出現(xiàn).

    2.2 風速對藍藻水華面積的影響

    風速對藍藻水華的面積具有顯著的影響,在風速較小的情況下,藻類會順著風向漂移,易形成藍藻水華,而較大的風速則會對藍藻水華產(chǎn)生一定程度的抑制作用.圖3為衛(wèi)星觀測到的太湖藍藻水華面積與對應(yīng)前6h平均風速分布圖,由圖3看出,太湖藍藻水華面積大多在600km2以下,占比達95%,總體上隨著藍藻水華面積的增大,相對應(yīng)的風速趨于減小.

    圖3 2003~2013年太湖藍藻水華面積及對應(yīng)6h平均風速分布Fig.3 Cyanobacteria bloom area and average wind speed in 6h before satellites observation time during 2003~2013

    圖4 2003~2013年藍藻水華不同面積區(qū)間及對應(yīng)6h平均風速Fig.4 Cyanobacteria bloom area ranges and corresponding average wind speed in 6h during 2003~2013

    為進一步分析太湖藍藻水華面積與風速之間的關(guān)系,根據(jù)每個區(qū)間之間樣本數(shù)等分原則,將面積分成0~13、14~32、33~63、64~103、104~152、153~241、242~427及428以上km2共8個區(qū)間,分別統(tǒng)計區(qū)間內(nèi)前6h平均最小風速、最大風速及80%概率區(qū)間對應(yīng)的風速,繪制成相應(yīng)的風速箱線圖(圖4).由圖可以看出,各面積區(qū)間對應(yīng)的前6h平均最大風速隨面積增大而明顯下降,平均最小風速則隨面積增大而緩慢下降,風速波動范圍隨面積增大呈明顯的收窄趨勢,風速范圍由最大的1.1~3.6縮減至0.7~2.1m/s.

    2.3 風速對大范圍藍藻水華的影響

    表2 大范圍藍藻水華前1~6h平均風速范圍及藍藻水華頻次總體分布Table 2 Average wind speed in 1~6h before large scalecyanobacteria bloom observed by satellites andcyanobacteria bloom frequency

    2003~2013年衛(wèi)星觀測到大范圍藍藻水華樣本數(shù)共計62個,除2003年以外的其他各年均有發(fā)生.統(tǒng)計大范圍藍藻水華前太湖沿岸氣象站1~6h的平均風速分布情況,分析大范圍藍藻水華與風速的關(guān)系(表2),表明:衛(wèi)星觀測到大范圍藍藻水華前1~6h各時間段平均風速最小均小于0.5m/s,最大平均風速均小于2.9m/s;各時間段藍藻水華主要出現(xiàn)在平均風速為1.1~2m/s區(qū)間,頻次占比均超過50%,隨著統(tǒng)計的時間段延長,0~1m/s區(qū)間的藍藻水華頻次占比由1h的6%穩(wěn)定增至5h的35%,而2.1~3m/s區(qū)間的頻次占比則由30%穩(wěn)定降至11%,0~2m/s區(qū)間的頻次占比由68%增至89%.由此得到:大范圍藍藻水華主要出現(xiàn)在前6h平均風速小于2m/s的情況下,在風速小于0.5m/s或大于3.0m/s的情況下,衛(wèi)星幾乎觀測不到大面積藍藻水華的聚集.

    2.4 風向?qū)λ{藻水華的影響

    圖5 各季節(jié)太湖藍藻水華頻次分布及風向玫瑰Fig.5 Cyanobacteria bloom frequency and corresponding wind direction rose chart in the four seasons during 2003~2013

    風向主要影響藍藻水華的輸移方向,是影響太湖藍藻水華在全湖分布格局變化的重要驅(qū)動因素.圖5為2003~2013年各季節(jié)藍藻水華頻次分布及太湖區(qū)域風向玫瑰圖.由圖5可看出,在春季(3~5月)藍藻水華復(fù)蘇期,隨溫度上升以及光輻射加強,冬季累積在底泥表面的藍藻“種源”開始復(fù)蘇,并隨氣溫升高快速生長,形成藍藻水華.由于春季太湖區(qū)域盛行東南東風向,聚集成群的藍藻水華在風的作用下向西北部方向漂移,從而在太湖西北部沿岸區(qū)形成相對高頻區(qū).在藍藻水華活躍的夏季(6~8月),氣象條件適宜,藍藻水華頻次明顯增加,在主導(dǎo)風向仍為東南東風的作用下,藍藻水華向西北部漂移,與西部沿岸區(qū)、梅梁湖和竺山湖區(qū)域自身生長的藍藻共同形成水華,成為太湖藍藻水華的頻發(fā)重發(fā)區(qū).秋季(9~11月)仍處于藍藻水華活躍期,但主導(dǎo)風向已轉(zhuǎn)為北風,在相對高頻率的北~西北風的作用下,大量藍藻水華向南部漂移集聚,頻發(fā)重發(fā)區(qū)域向湖心區(qū)和南部沿岸區(qū)轉(zhuǎn)移.進入冬季(12月~次年2月),隨著氣溫的下降,藍藻水華逐漸進入休眠期,前期形成的大量藍藻水華不可能在短時間內(nèi)迅速消亡,在主導(dǎo)風向為NW的風作用下,藍藻水華高頻區(qū)域繼續(xù)向湖心區(qū)轉(zhuǎn)移.

    2.5 數(shù)值模擬近地面風場對藍藻水華的影響

    由于湖水和陸地之間存在熱力差異,對于太湖這樣的大型湖泊,主要受湖陸風的影響,湖面氣象觀測站點稀少及觀測資料缺乏,難以全面了解湖面風場的變化情況.利用WRF 3.5.1對太湖區(qū)域的風場進行高分辨率的數(shù)值模擬,分析近地面風場對藍藻水華輸移及分布范圍的影響,以期全面、詳細地了解風對藍藻水華的影響.

    以2010年9月4~7日的藍藻水華過程作為個例(表3).2010年9月4~7日,太湖湖區(qū)及沿岸地區(qū)日平均氣溫在29℃左右,日平均風速<3m/s,氣象條件適宜,衛(wèi)星均觀測到了藍藻水華,但面積和范圍不同,4日和6日藍藻水華面積分別達789km2和795km2,而7日僅有6km2.

    對照藍藻水華強度分布圖和近地面風場數(shù)值模擬結(jié)果(圖6),可以看到,在9月4日13:00時的風場模擬圖上,梅梁湖、西部沿岸區(qū)及湖心區(qū)北部區(qū)域為低風速區(qū),風速為2~3m/s,部分區(qū)域甚至更低,無明顯的主導(dǎo)風向,由湖內(nèi)向湖岸風速逐漸增大,風向由東北部偏西北向——湖心區(qū)中部偏東向——西北部沿岸偏東南向變化,受此影響,在湖心區(qū)、西部沿岸區(qū)和梅梁湖等區(qū)域觀測到大范圍藍藻水華,在貢湖和竺山湖大部、西部近岸等區(qū)域由于風速較大,藍藻水華范圍或強度明顯偏小;在9月5 日13:00時的風場模擬圖上,太湖湖面風速增大,僅在南部沿岸區(qū)和湖心區(qū)南部有2~3m/s的低風速區(qū),而這一區(qū)域由于水質(zhì)較好藍藻水華出現(xiàn)頻次較低,在主導(dǎo)風向為東北的風場影響下,僅在西南部沿岸區(qū)和湖心區(qū)觀測到小面積藍藻水華;9月6日11:00時的風場模擬圖顯示,在太湖西南部大片區(qū)域風速僅為1~2m/s,且湖面風向主要為東或東北向,在此風場作用下,觀測到西部沿岸區(qū)和湖心區(qū)西部出現(xiàn)大面積藍藻水華,而包括竺山湖、梅梁湖、貢湖、湖心區(qū)北部及西北部近岸區(qū)幾乎未出現(xiàn)藍藻水華;9月7日11:00時的風場模擬結(jié)果顯示,太湖湖面風速明顯增大,大部分區(qū)域風速都在4m/s以上,僅在南部沿岸區(qū)和東部沿岸區(qū)的很小區(qū)域風速較小,前一天的大面積藍藻水華已不見蹤影.由此可見:藍藻水華的面積和分布范圍對近地面風場的變化響應(yīng)非常迅速,較小的風速有利于藍藻顆粒的上浮聚集形成水華,較大的風速使藍藻顆粒不易上浮,并在風浪和環(huán)流的共同作用下于水體中混合無法集聚形成水華.

    表3 2010年9月4~7日的藍藻水華過程Table 3 A developing process of cyanobacteria bloom from 4th to 7th Sep. in 2010

    圖6 2010年9月4~7日藍藻水華分布及近地面風場數(shù)值模擬對照Fig.6 Cyanobacteria bloom distribution and corresponding surface wind simulation chart from Sep. 4th to 7th in 2010

    2.6 討論

    由于衛(wèi)星觀測到的主要是湖泊表面藍藻水華形成的現(xiàn)象,難以反映水華形成前藍藻生物量及其水中垂直分布,因而本文主要探討藍藻水華形成的風場特征及風對藍藻水華輸移和分布的影響.通常認為微風或無風有利于藍藻水華的形成,這種說法不夠準確.氣象上規(guī)定的“微風”風速范圍為3.4~5.4m/s,“無風”為0~0.2m/s,兩者相差甚遠.根據(jù)分析結(jié)果,風速<0.5m/s時藍藻水華的頻次占比均不超過2%,而大范圍藍藻水華前1~6h最小風速為0.4~0.5m/s,表明水華的形成需要一定風力及其產(chǎn)生的風浪和湖流的擾動,真正風平浪靜時藍藻反而不易聚集形成水華.目前報道的大部分風速閾值都處于“微風”范圍,如Cao等[33]確定為3.1m/s,而孫小靜等[17]、王成林等[34]為<4m/s,小型湖泊則更小,如淀山湖為<2.2m/s[22].造成差異的原因很多,數(shù)據(jù)來自于實驗還是實地觀測、觀測位置和取樣不同,其結(jié)果就會不一樣,此外,風速統(tǒng)計采用的方法不同也會影響結(jié)果.統(tǒng)計風速的時間長短對閾值的影響很大,總體上統(tǒng)計時間越長越容易平滑風的波動性,閾值也趨小,如大范圍藍藻水華發(fā)生前1h的平均風速主要位于1.1~2.9m/s區(qū)間,占比高達94%,而前6h平均風速則主要位于0~2.0m/s,占比89%,因此對于各種風速閾值需要進行具體分析.

    藍藻水華的形成是包括氣象、水文等眾多環(huán)境因子綜合影響的結(jié)果,任何單一因子都難以驅(qū)動藍藻的生長和水華形成[10].本文的研究沒有考慮其他環(huán)境因子,因此分析結(jié)果難免會受到影響,需要進一步研究其他環(huán)境因子及其協(xié)同作用對藍藻水華形成的影響,以更好地理解藍藻水華的發(fā)生機理,為防控藍藻水華提供理論依據(jù).

    3 結(jié)論

    3.1 藍藻水華的形成需要一定風力及其產(chǎn)生的風浪和湖流的擾動,太湖藍藻水華主要出現(xiàn)于前6h平均風速為0.5~3.4m/s的區(qū)間,衛(wèi)星觀測到藍藻水華時刻前6h平均風速在此區(qū)間占比達94.7%.

    3.2 藍藻水華面積總體上隨風速增大而減小,大范圍藍藻水華主要出現(xiàn)在前1h平均風速1.1~ 2.9m/s區(qū)間,占比達94%,而6h平均風速為0~ 2.0m/s,占比89%.

    3.3 風向主要影響藍藻水華的分布格局.在春季復(fù)蘇期和夏季活躍期,主導(dǎo)風向為東南東,藍藻水華主要向西北部方向漂移,與西部沿岸區(qū)、梅梁湖和竺山湖區(qū)域自身生長的藍藻共同形成藍藻水華,成為藍藻水華頻發(fā)重發(fā)區(qū);在秋季活躍期,主導(dǎo)風向轉(zhuǎn)為北風,大量藍藻水華向太湖南部漂移集聚.12月份藍藻水華逐漸進入衰敗期,前期形成的大量藍藻水華在主導(dǎo)風向為西北的作用下,繼續(xù)向太湖湖心區(qū)轉(zhuǎn)移.

    參考文獻:

    [1] O’Neila J M, Davisb T W, Burfordb M A, et al. The rise of harmful cyanobacteria blooms: The potential roles of eutrophication and climate change [J]. Harmful Algae, 2012,14:313–334.

    [2] Te S H, Chen E Y, Gin K Y. Comparison of quantitative PCR and droplet digital PCR multiplex assays for two genera of bloomforming cyanobacteria, cylindrospermopsis and Microcystis [J]. Appl. Environ. Microbiol., 2015,81(15):5203-5211.

    [3] Anna Rigosi, Cayelan C Carey, Bas W Ibelings, et al. The interaction between climate warming and eutrophication to promote cyanobacteria is dependent on trophic state and varies among taxa [J]. Limnology and Oceanography, 2014,59(1):99-114.

    [4] 商兆堂,任 健,秦銘榮.氣候變化與太湖藍藻暴發(fā)的關(guān)系 [J].生態(tài)學(xué)雜志, 2010,29(1):55-56.

    [5] 王成林,潘維玉,韓月琪,等.全球氣候變化對太湖藍藻水華發(fā)展演變的影響 [J]. 中國環(huán)境科學(xué), 2010,30(6):822-828.

    [6] 孔繁翔,高 光.大型淺水富營養(yǎng)化湖泊中藍藻水華形成機理的思考 [J]. 生態(tài)學(xué)報, 2005,25(3):589-595.

    [7] Bista D, Heckathorn S A, Bridgeman T, et al. Interactive effects of temperature, nitrogen, and zooplankton on growth and protein and carbohydrate content of cyanobacteria from Western Lake Erie [J]. Journal of Water Resource and Protection, 2014,6(12): 1139-1153.

    [8] 趙巧華,杭蓉蓉,王 玲,等.夏季太湖水體勢能異常的空間分異特征及其機制 [J]. 中國環(huán)境科學(xué), 2015,35(1):227-235.

    [9] 謝國清,李 蒙,魯韋坤,等.滇池藍藻水華光譜特征、遙感識別及暴發(fā)氣象條件 [J]. 湖泊科學(xué), 2010,22(3):327-336.

    [10] 孔繁翔,宋立榮,等.藍藻水華形成過程及其環(huán)境特征研究 [M].北京:科學(xué)出版社, 2011.

    [11] Zhang Min, Duan Hongtao, Shi Xiaoli, et al. Contributions of meteorology to the phenology of cyanobacterial blooms: implications for future climate change [J]. Wat. Res., 2012,46(2):442-452.

    [12] Reichwaldt E S, Ghadouani A. Effects of rainfall patterns on toxic cyanobacterial blooms in a changing climate: between simplistic scenarios and complex dynamics [J]. Wat. Res., 2011,46(5):1372-1393.

    [13] 巫 娟,陳雪初,孔海南,等.光照度對水華魚腥藻細胞比重與藻絲長度的影響研究 [J]. 中國環(huán)境科學(xué), 2012,32(5):875-879.

    [14] 秦伯強,胡維平,高 光,等.太湖沉積物懸浮的動力機制及內(nèi)源釋放的概念性模式 [J]. 科學(xué)通報, 2003,48(17):1822-1831.

    [15] 范成新,張 路,秦伯強,等.風浪作用下太湖懸浮態(tài)顆粒物中磷的動態(tài)釋放估算 [J]. 中國科學(xué)(D輯), 2003,33(8):760-768.

    [16] 朱廣偉,秦伯強,高 光,等.風浪擾動引起大型淺水湖泊內(nèi)源磷暴發(fā)性釋放的直接證據(jù) [J]. 科學(xué)通報, 2005,50(1):66-71.

    [17] 孫小靜,秦伯強,朱廣偉,等.風浪對太湖水體中膠體態(tài)營養(yǎng)鹽和浮游植物的影響 [J]. 環(huán)境科學(xué), 2007,28(3):506-511.

    [18] 尤本勝,王同成,范成新,等.風浪作用下太湖草型湖區(qū)水體N、P動態(tài)負荷模擬 [J]. 中國環(huán)境科學(xué), 2008,28(1):33-38.

    [19] 朱永春,蔡啟銘.風場對藻類在太湖中遷移影響的動力學(xué)研究[J]. 湖泊科學(xué), 1997,9(2):152-158.

    [20] Kanoshina I, Lips U, Leppanen L M. The influence of weather conditions (temperature and wind) on cyanobacterial bloom development in the Gulf of Finland (Baltic Sea) [J]. Harmful Algae, 2003,(2):29-41.

    [21] Zhou Jian, Qin Boqiang, Casenave C, et al. Effects of wind wave turbulence on the phytoplankton community composition in large, shallow Lake Taihu [J]. Environmental Science and Pollution Research, 2015,22(16):12737-12746.

    [22] 王銘瑋,徐啟新,車 越,等.淀山湖藍藻水華暴發(fā)的氣象水文因素探討 [J]. 華東師范大學(xué)學(xué)報(自然科學(xué)版), 2011,2011(1):21-31.

    [23] Wu Tingfeng, Qin Boqiang, Zhu Guangwei, et al. Dynamics of cyanobacteria bloom formation during short-term hydrodynamic fluctuation in a large shallow, eutrophic, and wind-exposed Lake Taihu, China [J]. Environmental Science and Pollution Research, 2013,20(12):8546-8556.

    [24] Wu Tingfeng, Qin Boqiang, Justin D, et al. The influence of changes in wind patterns on the areal extension of surface cyanobacterial blooms in a large shallow lake in China [J]. Science of The Total Environment, 2015,(518–519):24-30.

    [25] 王成林,陳黎明,潘維玉,等.適宜太湖藍藻水華形成的風場輻散特征及其形成機制 [J]. 中國環(huán)境科學(xué), 2010,30(9):1168-1176.

    [26] 王成林,黃 娟,錢 新.高溫微風條件下太湖流域風場時空特征分析 [J]. 湖泊科學(xué), 2011,23(1):122-128.

    [27] 王文蘭,曾明劍,任 健.近地面風場變化對太湖藍藻暴發(fā)影響的數(shù)值研究 [J]. 氣象科學(xué), 2011,31(6):718-725.

    [28] 王長友,于 洋,孫運坤,等.基于ELCOM-CAEDYM模型的太湖藍藻水華早期預(yù)測探討 [J]. 中國環(huán)境科學(xué), 2013,33(3): 491-502.

    [29] 馬榮華,孔維娟,段洪濤,等.基于MODIS影像估測太湖藍藻暴發(fā)期藻藍素含量 [J]. 中國環(huán)境科學(xué), 2009,29(3):254-260.

    [30] 韓秀珍,吳朝陽,鄭 偉,等.基于水面實測光譜的太湖藍藻衛(wèi)星遙感研究[J].應(yīng)用氣象學(xué)報, 2010,21(6):724-731.

    [31] Hu C M, Lee Z P, Ma R H, et al. Moderate Resolution Imaging Spectroradiometer (MODIS) observations of cyanobacteria blooms in Taihu Lake, China [J]. Journal of Geophysical Research, 2010,115(C04002):1-20.

    [32] 中國氣象局.地面氣象觀測規(guī)范 [M]. 北京:氣象出版社, 2003.

    [33] Cao H S, Kong F X, Luo L C, et al. Effects of wind-induced waves on vertical photoplankton distribution and surface blooms of Microcystis aeruginosain Lake Taihu [J]. Journal of Freshwater Ecology, 2006,21(2):231-238.

    [34] 王成林.太湖藍藻水華及藻源性“湖泛”形成機制的氣象因素研究[學(xué)位論文].南京:南京大學(xué), 2008.

    Analysis of wind field features causing cyanobacteria bloom in Taihu Lake combined with remote sensing methods.

    LI Ya-chun1*, XIE Xiao-ping1, HANG Xing1, ZHU Xiao-li2, HUANG Shan1, JING Yuan-shu2(1. Meteorological Bureau of Jiangsu Province, Nanjing 210008, China;2.Meteorological Bureau of Changzhou City, Changzhou 213022, China;3.Institute of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210044, China). China Environmental Science, 2016,36(2):525~533

    Abstract:On the basis of the meteorological observation data and remote sensing data over the Taihu Lake and its surrounding areas during 2003 to 2013, combined with the WRF3.5.1 simulations, the relationships between the cyanobacteria blooms in Taihu Lake and surface wind speeds were analyzed. The results show that cyanobacteria blooms mainly occurred in the period of the weak wind speed about 0.5~3.4m/s, about six hours earlier than the MODIS satellite overpass time in Taihu Lake, which accounts for 94.7% of all cyanobacteria blooms events. As the area of cyanobacteria blooms decreasing with the wind speed strengthening, a large bloom of cyanobacteria bloom mainly occurred 6 hours ahead with the average wind speed less than 2m/s, accounting for 89%. In addition, the wind direction mainly affects the spatial patterns of cyanobacteria blooms in Taihu Lake. All these results suggest that the wind patterns play an important role on the formation, drifting and distribution of cyanobacteria blooms.

    Key words:cyanobacteria bloom;surface wind speeds;remote sensing;numerical simulation

    作者簡介:李亞春(1966-),男,江蘇武進人,研究員級高級工程師,碩士,主要從事遙感應(yīng)用和應(yīng)用氣象研究.發(fā)表論文20余篇.

    基金項目:江蘇省基礎(chǔ)研究計劃太湖專項(BK2007745);江蘇省科技支撐計劃項目(BE2011840);江蘇省氣象局重點項目(KZ201403)

    收稿日期:2015-07-15

    中圖分類號:X87

    文獻標識碼:A

    文章編號:1000-6923(2016)02-0525-09

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風場信息的風機輪轂處風速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    日韩欧美三级三区| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 欧美在线一区亚洲| 精品一区二区三卡| 9色porny在线观看| 狠狠婷婷综合久久久久久88av| 久热爱精品视频在线9| 国产1区2区3区精品| 精品国产一区二区久久| 99国产综合亚洲精品| 国产蜜桃级精品一区二区三区 | 国产亚洲精品第一综合不卡| 亚洲色图av天堂| 欧美久久黑人一区二区| 国产精品秋霞免费鲁丝片| 啪啪无遮挡十八禁网站| 午夜福利乱码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕最新亚洲高清| a在线观看视频网站| 久久天堂一区二区三区四区| 老司机靠b影院| 正在播放国产对白刺激| 日本撒尿小便嘘嘘汇集6| 久久久久国内视频| 老司机影院毛片| 日本精品一区二区三区蜜桃| 99re在线观看精品视频| av免费在线观看网站| 色播在线永久视频| 国产成人免费无遮挡视频| 欧美中文综合在线视频| 国产麻豆69| 成人手机av| 亚洲精品国产精品久久久不卡| 巨乳人妻的诱惑在线观看| 老熟女久久久| 97人妻天天添夜夜摸| 欧美不卡视频在线免费观看 | 国产不卡一卡二| 亚洲午夜精品一区,二区,三区| 日日爽夜夜爽网站| 亚洲 欧美一区二区三区| 正在播放国产对白刺激| 很黄的视频免费| 老鸭窝网址在线观看| 在线天堂中文资源库| 免费少妇av软件| 动漫黄色视频在线观看| 777久久人妻少妇嫩草av网站| 精品第一国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 国产单亲对白刺激| 无人区码免费观看不卡| 又大又爽又粗| 国产精品久久视频播放| 中文字幕另类日韩欧美亚洲嫩草| 一二三四社区在线视频社区8| 宅男免费午夜| 亚洲视频免费观看视频| 日韩 欧美 亚洲 中文字幕| 韩国精品一区二区三区| 国产高清视频在线播放一区| 天堂俺去俺来也www色官网| 手机成人av网站| 欧美性长视频在线观看| 国产1区2区3区精品| 国产在线精品亚洲第一网站| 正在播放国产对白刺激| 亚洲欧洲精品一区二区精品久久久| 国产精品亚洲av一区麻豆| 女人精品久久久久毛片| 亚洲熟女精品中文字幕| 亚洲一区二区三区欧美精品| 欧美日韩亚洲综合一区二区三区_| 久久久久精品人妻al黑| 精品一区二区三卡| 大型黄色视频在线免费观看| 变态另类成人亚洲欧美熟女 | 国产欧美日韩一区二区三| 大型黄色视频在线免费观看| 热re99久久精品国产66热6| 激情在线观看视频在线高清 | 在线天堂中文资源库| 久久精品aⅴ一区二区三区四区| 日本五十路高清| 动漫黄色视频在线观看| 又紧又爽又黄一区二区| 久久精品亚洲熟妇少妇任你| 欧美日韩国产mv在线观看视频| 日日爽夜夜爽网站| 一级毛片女人18水好多| 国产高清国产精品国产三级| 夫妻午夜视频| 欧洲精品卡2卡3卡4卡5卡区| 夜夜夜夜夜久久久久| 午夜老司机福利片| 99热国产这里只有精品6| av福利片在线| 超色免费av| 高清毛片免费观看视频网站 | av中文乱码字幕在线| 亚洲 国产 在线| 国产在线一区二区三区精| 国产无遮挡羞羞视频在线观看| 又紧又爽又黄一区二区| 午夜激情av网站| 免费在线观看视频国产中文字幕亚洲| 亚洲 欧美一区二区三区| 美女福利国产在线| 精品高清国产在线一区| 法律面前人人平等表现在哪些方面| 两性夫妻黄色片| 国产欧美日韩综合在线一区二区| 国产av又大| 18在线观看网站| 久久人妻福利社区极品人妻图片| 午夜福利在线免费观看网站| 日韩熟女老妇一区二区性免费视频| 国产精品久久久人人做人人爽| 久久久久久人人人人人| 久久久久久久国产电影| 亚洲全国av大片| 黄网站色视频无遮挡免费观看| 国产亚洲av高清不卡| av免费在线观看网站| 激情视频va一区二区三区| 国产aⅴ精品一区二区三区波| 777米奇影视久久| 国产欧美日韩一区二区精品| 亚洲黑人精品在线| 国产一卡二卡三卡精品| 久久婷婷成人综合色麻豆| 国产成人av激情在线播放| 国产精品 欧美亚洲| 亚洲av第一区精品v没综合| 黄色 视频免费看| 中文亚洲av片在线观看爽 | 高清视频免费观看一区二区| 精品人妻在线不人妻| 少妇粗大呻吟视频| 国产精品国产高清国产av | 91字幕亚洲| 精品少妇一区二区三区视频日本电影| 美女福利国产在线| 一级片'在线观看视频| 国产成人欧美在线观看 | 热99国产精品久久久久久7| 91大片在线观看| 免费久久久久久久精品成人欧美视频| 狂野欧美激情性xxxx| 久久久久国内视频| 老鸭窝网址在线观看| 嫁个100分男人电影在线观看| 99国产综合亚洲精品| 老司机影院毛片| 久久中文字幕人妻熟女| 中文字幕人妻丝袜一区二区| 久久久久精品人妻al黑| 午夜免费观看网址| 午夜福利乱码中文字幕| 热99久久久久精品小说推荐| 在线观看免费视频网站a站| 亚洲精品在线观看二区| 天天躁夜夜躁狠狠躁躁| 国产亚洲精品久久久久久毛片 | 欧美日韩中文字幕国产精品一区二区三区 | 老司机在亚洲福利影院| 国产精品乱码一区二三区的特点 | 日韩制服丝袜自拍偷拍| 99国产精品一区二区三区| 又黄又粗又硬又大视频| 亚洲 欧美一区二区三区| 国产成人免费无遮挡视频| 国产精品偷伦视频观看了| 欧美日韩亚洲综合一区二区三区_| 一进一出好大好爽视频| 无遮挡黄片免费观看| 99re在线观看精品视频| 国产精品二区激情视频| 天堂动漫精品| 亚洲av美国av| 欧美午夜高清在线| 亚洲精品在线美女| 国产高清视频在线播放一区| 三上悠亚av全集在线观看| 欧美+亚洲+日韩+国产| 丝瓜视频免费看黄片| 天天躁夜夜躁狠狠躁躁| 涩涩av久久男人的天堂| 久久久久久久国产电影| 国产野战对白在线观看| 亚洲第一av免费看| 婷婷精品国产亚洲av在线 | 久久久国产精品麻豆| 亚洲精品美女久久av网站| 人人妻人人澡人人爽人人夜夜| 国内毛片毛片毛片毛片毛片| 国产精品 国内视频| 三级毛片av免费| 婷婷丁香在线五月| 久久热在线av| 高潮久久久久久久久久久不卡| 日日摸夜夜添夜夜添小说| 欧美不卡视频在线免费观看 | www.999成人在线观看| 欧美在线黄色| 啦啦啦 在线观看视频| 日韩欧美免费精品| 视频区欧美日本亚洲| 亚洲精品美女久久久久99蜜臀| 另类亚洲欧美激情| 精品视频人人做人人爽| 欧美黄色淫秽网站| 女人精品久久久久毛片| 老司机午夜福利在线观看视频| 老司机午夜福利在线观看视频| 免费av中文字幕在线| 校园春色视频在线观看| 亚洲在线自拍视频| 亚洲中文日韩欧美视频| 黑人操中国人逼视频| 欧美精品啪啪一区二区三区| av电影中文网址| 国产精品99久久99久久久不卡| 激情在线观看视频在线高清 | 成在线人永久免费视频| 国产成人精品在线电影| 国产男女超爽视频在线观看| 日韩人妻精品一区2区三区| 女警被强在线播放| 中出人妻视频一区二区| 国产激情欧美一区二区| 国产成+人综合+亚洲专区| 欧美日韩黄片免| 村上凉子中文字幕在线| 国产精品香港三级国产av潘金莲| 亚洲成人国产一区在线观看| 国产免费现黄频在线看| 少妇粗大呻吟视频| 国产精品国产av在线观看| 国产精品免费一区二区三区在线 | 黄色怎么调成土黄色| 欧洲精品卡2卡3卡4卡5卡区| 黑人欧美特级aaaaaa片| 欧美乱码精品一区二区三区| 亚洲欧美激情在线| 乱人伦中国视频| 1024香蕉在线观看| 国产欧美日韩一区二区三| 亚洲av成人av| 女人久久www免费人成看片| 丝袜美腿诱惑在线| 久久九九热精品免费| 国产国语露脸激情在线看| 99久久综合精品五月天人人| 美女扒开内裤让男人捅视频| 中文字幕高清在线视频| 老司机在亚洲福利影院| 午夜视频精品福利| 午夜福利欧美成人| videosex国产| videosex国产| 亚洲一区二区三区欧美精品| 欧美精品一区二区免费开放| 91麻豆av在线| 国产1区2区3区精品| 老司机在亚洲福利影院| 国产精品久久久久久精品古装| av天堂在线播放| 美女高潮喷水抽搐中文字幕| 天堂中文最新版在线下载| 久久精品国产亚洲av高清一级| 99精品久久久久人妻精品| 50天的宝宝边吃奶边哭怎么回事| 老司机在亚洲福利影院| 午夜精品国产一区二区电影| 国产不卡一卡二| 成人亚洲精品一区在线观看| 叶爱在线成人免费视频播放| 黄色视频,在线免费观看| 精品国产一区二区三区久久久樱花| 免费观看a级毛片全部| 欧美日韩成人在线一区二区| 久久午夜亚洲精品久久| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 12—13女人毛片做爰片一| 欧美精品高潮呻吟av久久| 久久久精品免费免费高清| 欧美日韩亚洲综合一区二区三区_| 一进一出抽搐动态| 国产亚洲欧美98| 国产高清国产精品国产三级| 777久久人妻少妇嫩草av网站| 国产免费av片在线观看野外av| 露出奶头的视频| 久99久视频精品免费| 国产高清国产精品国产三级| 激情在线观看视频在线高清 | 黄频高清免费视频| 久久久国产成人精品二区 | 丁香六月欧美| 亚洲三区欧美一区| 精品一区二区三卡| 免费观看人在逋| 一级作爱视频免费观看| 精品欧美一区二区三区在线| 亚洲自偷自拍图片 自拍| 亚洲精品国产区一区二| 80岁老熟妇乱子伦牲交| 俄罗斯特黄特色一大片| 十分钟在线观看高清视频www| 99久久99久久久精品蜜桃| 新久久久久国产一级毛片| 亚洲av美国av| 亚洲七黄色美女视频| 久久午夜综合久久蜜桃| 欧美乱码精品一区二区三区| 午夜免费成人在线视频| 国产伦人伦偷精品视频| 免费不卡黄色视频| 99国产精品一区二区三区| 午夜福利视频在线观看免费| 大香蕉久久网| 中文欧美无线码| 99久久99久久久精品蜜桃| 亚洲专区中文字幕在线| 国产一卡二卡三卡精品| 欧美在线一区亚洲| 成人18禁在线播放| 午夜两性在线视频| 亚洲视频免费观看视频| 青草久久国产| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区mp4| 午夜福利乱码中文字幕| 满18在线观看网站| 操美女的视频在线观看| 色综合婷婷激情| 亚洲成人免费电影在线观看| 亚洲精品久久成人aⅴ小说| 欧美精品人与动牲交sv欧美| 伦理电影免费视频| 国产亚洲一区二区精品| 亚洲第一青青草原| 人成视频在线观看免费观看| e午夜精品久久久久久久| 天天躁日日躁夜夜躁夜夜| 国产日韩欧美亚洲二区| 亚洲国产精品一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 精品电影一区二区在线| 久久中文字幕一级| 成人三级做爰电影| 搡老熟女国产l中国老女人| 激情在线观看视频在线高清 | 两性夫妻黄色片| 欧美 亚洲 国产 日韩一| 免费少妇av软件| 俄罗斯特黄特色一大片| 国产区一区二久久| 亚洲成人手机| 韩国精品一区二区三区| 99精国产麻豆久久婷婷| 午夜老司机福利片| 中文字幕精品免费在线观看视频| 亚洲精品久久成人aⅴ小说| 男人的好看免费观看在线视频 | 精品国产美女av久久久久小说| 老鸭窝网址在线观看| 淫妇啪啪啪对白视频| 国产无遮挡羞羞视频在线观看| 欧美激情极品国产一区二区三区| 久久精品成人免费网站| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 两个人免费观看高清视频| 日韩大码丰满熟妇| 久久国产精品人妻蜜桃| 精品福利观看| 日韩欧美一区视频在线观看| 高清毛片免费观看视频网站 | 欧美激情高清一区二区三区| 国产单亲对白刺激| 欧美不卡视频在线免费观看 | 精品高清国产在线一区| 色在线成人网| 狠狠狠狠99中文字幕| 亚洲精品国产一区二区精华液| 在线观看免费视频日本深夜| 日韩三级视频一区二区三区| 一进一出抽搐动态| 王馨瑶露胸无遮挡在线观看| 极品少妇高潮喷水抽搐| 成在线人永久免费视频| 丰满迷人的少妇在线观看| 在线观看66精品国产| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人| 欧美国产精品一级二级三级| 亚洲欧美日韩另类电影网站| 亚洲成av片中文字幕在线观看| 视频区图区小说| 国产精品乱码一区二三区的特点 | 最新美女视频免费是黄的| 日本黄色视频三级网站网址 | 国产片内射在线| 国产无遮挡羞羞视频在线观看| 久久久久久久午夜电影 | 高清毛片免费观看视频网站 | 两性夫妻黄色片| 亚洲在线自拍视频| 亚洲 国产 在线| 国产亚洲精品久久久久5区| 大片电影免费在线观看免费| 久久久久国产精品人妻aⅴ院 | 国产单亲对白刺激| 亚洲欧美激情在线| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 精品福利观看| 高清视频免费观看一区二区| 波多野结衣一区麻豆| 天天躁日日躁夜夜躁夜夜| 精品人妻在线不人妻| 十八禁人妻一区二区| 午夜日韩欧美国产| 一进一出好大好爽视频| 欧美日韩一级在线毛片| 国产精品乱码一区二三区的特点 | 正在播放国产对白刺激| 欧美日韩黄片免| 色94色欧美一区二区| 校园春色视频在线观看| а√天堂www在线а√下载 | 9热在线视频观看99| 国产精品香港三级国产av潘金莲| 青草久久国产| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 男人的好看免费观看在线视频 | 中文字幕另类日韩欧美亚洲嫩草| 成人影院久久| 好男人电影高清在线观看| 性色av乱码一区二区三区2| 一a级毛片在线观看| 亚洲精品美女久久久久99蜜臀| 国产成人免费观看mmmm| 精品一区二区三区视频在线观看免费 | 色综合欧美亚洲国产小说| 飞空精品影院首页| 国产成人精品久久二区二区91| 亚洲色图综合在线观看| 亚洲欧美色中文字幕在线| 青草久久国产| 亚洲色图av天堂| 国产免费现黄频在线看| 中文字幕最新亚洲高清| 午夜影院日韩av| 一区在线观看完整版| 老熟妇乱子伦视频在线观看| 国产亚洲欧美在线一区二区| 久久午夜综合久久蜜桃| 午夜影院日韩av| 欧美日韩一级在线毛片| 曰老女人黄片| 成人免费观看视频高清| 欧美一级毛片孕妇| 成人特级黄色片久久久久久久| 国产亚洲精品久久久久5区| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 中文字幕人妻熟女乱码| 欧美日韩精品网址| 狠狠婷婷综合久久久久久88av| 黄色怎么调成土黄色| 日韩 欧美 亚洲 中文字幕| 精品久久久久久电影网| 国产精品久久久久久精品古装| 亚洲av第一区精品v没综合| 亚洲国产欧美网| 国产一区二区激情短视频| av天堂久久9| 在线观看一区二区三区激情| 午夜精品久久久久久毛片777| 又黄又爽又免费观看的视频| 国产成人精品在线电影| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 三级毛片av免费| 女人被躁到高潮嗷嗷叫费观| 又大又爽又粗| 欧美 亚洲 国产 日韩一| 91成年电影在线观看| 欧美大码av| 国产不卡一卡二| 欧美丝袜亚洲另类 | 国产精品欧美亚洲77777| 男女午夜视频在线观看| 亚洲熟女精品中文字幕| 麻豆乱淫一区二区| 国产伦人伦偷精品视频| 男男h啪啪无遮挡| www.999成人在线观看| 欧美大码av| 日韩 欧美 亚洲 中文字幕| 国产97色在线日韩免费| 好男人电影高清在线观看| 麻豆国产av国片精品| 亚洲精品一卡2卡三卡4卡5卡| 满18在线观看网站| 老司机深夜福利视频在线观看| 人人妻人人澡人人爽人人夜夜| 看免费av毛片| 人妻一区二区av| 正在播放国产对白刺激| 在线观看免费视频网站a站| 国产单亲对白刺激| 新久久久久国产一级毛片| 精品国产一区二区久久| 12—13女人毛片做爰片一| 欧美色视频一区免费| 精品国产超薄肉色丝袜足j| 亚洲一区二区三区欧美精品| 国产精品1区2区在线观看. | 精品午夜福利视频在线观看一区| 91成人精品电影| 亚洲成a人片在线一区二区| 久久精品国产a三级三级三级| bbb黄色大片| 亚洲 国产 在线| 国产熟女午夜一区二区三区| 亚洲熟女毛片儿| 国产精华一区二区三区| 亚洲欧美日韩另类电影网站| av网站免费在线观看视频| 一边摸一边抽搐一进一出视频| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 少妇裸体淫交视频免费看高清 | 亚洲中文av在线| 国产成人精品久久二区二区免费| 亚洲中文av在线| 热re99久久国产66热| 大陆偷拍与自拍| 久久精品国产综合久久久| 亚洲av电影在线进入| 色播在线永久视频| 亚洲人成77777在线视频| 久久国产亚洲av麻豆专区| 色94色欧美一区二区| 一区二区三区国产精品乱码| 国产精品.久久久| 成人精品一区二区免费| 女人久久www免费人成看片| 丝袜美足系列| 黑丝袜美女国产一区| 香蕉久久夜色| а√天堂www在线а√下载 | 狠狠狠狠99中文字幕| 手机成人av网站| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 在线观看午夜福利视频| 国产亚洲一区二区精品| 国产精华一区二区三区| 欧美 日韩 精品 国产| 在线观看日韩欧美| 亚洲伊人色综图| 女人精品久久久久毛片| 欧美精品高潮呻吟av久久| av网站在线播放免费| 国产高清激情床上av| 国产精品 欧美亚洲| 热99久久久久精品小说推荐| 51午夜福利影视在线观看| 在线观看舔阴道视频| 狠狠婷婷综合久久久久久88av| 亚洲色图综合在线观看| 国产一区二区三区综合在线观看| 亚洲精品av麻豆狂野| 黄色视频,在线免费观看| 又黄又粗又硬又大视频| 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 欧美乱妇无乱码| 国产高清视频在线播放一区| aaaaa片日本免费| 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 亚洲视频免费观看视频| 亚洲欧美日韩高清在线视频| 一a级毛片在线观看| 久久久精品国产亚洲av高清涩受| 极品少妇高潮喷水抽搐| 久久久国产成人精品二区 | 欧美在线黄色| av一本久久久久| 欧美日本中文国产一区发布| 9191精品国产免费久久| 亚洲av第一区精品v没综合| 亚洲五月色婷婷综合| 免费观看a级毛片全部| 国产精品影院久久| 色综合欧美亚洲国产小说| 亚洲片人在线观看| 久久久久久久久免费视频了| videosex国产| 757午夜福利合集在线观看| 这个男人来自地球电影免费观看| 每晚都被弄得嗷嗷叫到高潮| 一级作爱视频免费观看| 男女午夜视频在线观看|