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

    基于GOCI數(shù)據(jù)的霾天氣氣溶膠輻射強(qiáng)迫的日內(nèi)變化——以長江三角洲為例

    2019-02-27 09:12:12李雅雯張海龍金雙根
    中國環(huán)境科學(xué) 2019年2期
    關(guān)鍵詞:反照率氣溶膠反演

    李雅雯,陳 健*,張海龍,金雙根

    ?

    基于GOCI數(shù)據(jù)的霾天氣氣溶膠輻射強(qiáng)迫的日內(nèi)變化——以長江三角洲為例

    李雅雯1,陳 健1*,張海龍2,金雙根1

    (1.南京信息工程大學(xué)遙感與測繪工程學(xué)院,江蘇 南京 210044;2.南京信息工程大學(xué)海洋科學(xué)學(xué)院,江蘇 南京 210044)

    以長江三角洲(以下簡稱長三角)為研究區(qū),基于2017年10~12月內(nèi)共16d 94景的晴空靜止衛(wèi)星 GOCI L1B圖像和6S模型,采用深藍(lán)算法反演長三角地區(qū)氣溶膠光學(xué)厚度(AOD),并利用2個實(shí)測站數(shù)據(jù)進(jìn)行驗(yàn)證.再通過反演得到的AOD結(jié)果進(jìn)一步模擬計算該區(qū)域地表和大氣層頂?shù)臍馊苣z直接輻射效應(yīng)(ADRE),并結(jié)合典型的霧霾天氣進(jìn)行分析.結(jié)果表明:利用GOCI數(shù)據(jù)反演的AOD具有較高擬合精度,北辰樓站點(diǎn)擬合相關(guān)性2為0.68,太湖站點(diǎn)為0.67.空間上,由于氣溶膠存在制冷效應(yīng),AOD和地表面、大氣層頂氣溶膠直接輻射強(qiáng)度均呈現(xiàn)出顯著的線性關(guān)系.時間上,由于氣溶膠擴(kuò)散和風(fēng)向等因素,早上和下午輻射效應(yīng)強(qiáng)度較高,中午較低.其中10:00和11:00影像可以較好地模擬日均ADRE的特征,利用ADRE日內(nèi)變化成功捕捉到一次霧霾爆發(fā)并消失的現(xiàn)象,對氣象部門監(jiān)測近地表氣溫和分析灰霾的形成等具有重要意義.

    灰霾;氣溶膠;直接輻射效應(yīng);長江三角洲;GOCI

    氣溶膠光學(xué)厚度(AOD)是氣溶膠最重要的光學(xué)特性之一,可以用于PM2.5[1]等顆粒物質(zhì)濃度估測以及氣溶膠氣候效應(yīng)研究等.氣溶膠直接輻射效應(yīng)(ADRE)是指氣溶膠組分含量的變化所產(chǎn)生的地氣系統(tǒng)能量平衡的變化[2],單位為W/m2.IPCC第五次會議報告中強(qiáng)調(diào)了氣溶膠的氣候效應(yīng)研究的重要性,而氣溶膠的日內(nèi)變化對于氣溶膠的氣候效應(yīng)評估目前仍存在較大的不確定性.

    目前,已有少數(shù)專家學(xué)者開展了ADRE的日內(nèi)變化研究.Arola等[3]利用大量的站點(diǎn)數(shù)據(jù)分析AOD的日變化對大氣層頂(TOA)日均ADRE的影響,發(fā)現(xiàn)單個站點(diǎn)的ADRE估計上出現(xiàn)明顯的偏移;當(dāng)對所有站點(diǎn)取平均時,AOD日變化對日均ADRE的影響很小,相對誤差主要在5%以內(nèi).Wang等[4]研究合肥區(qū)域發(fā)現(xiàn)氣溶膠日變化在氣溶膠濃度較高時確實(shí)會影響日均輻射強(qiáng)迫,并分析了3種氣溶膠光學(xué)性質(zhì),結(jié)果發(fā)現(xiàn)當(dāng)SSA和ASY參數(shù)范圍分別為0.85~ 0.99和0.61~0.73之間時,24-h平均輻射強(qiáng)迫的絕對誤差和相對誤差分別為0.3±4.2W/m2(?7.6~15.6W/ m2)和11.8%±6.6% W/m2(0.1%~28.5% W/m2). Kassianov等[5]發(fā)現(xiàn)即使在AOD日變特征強(qiáng)烈的情況下,僅需氣溶膠光學(xué)特性的日均值就能精確估測ADRE的日均值,并進(jìn)一步分析了氣溶膠光學(xué)特性的日變化對估測日均ADRE的影響.近年來,ADRE的衛(wèi)星遙感評估研究雖然已取得重要進(jìn)展,但針對ADRE的日變化特征的研究才剛剛起步,所采用的數(shù)據(jù)例如MODIS時相較少,限制了氣溶膠直接輻射效應(yīng)日變化特征的研究.

    長江三角洲地區(qū)作為我國東部沿海重要的發(fā)展區(qū)域,其經(jīng)濟(jì)和人口迅猛增長、工業(yè)化程度高度飽和[6],是我國污染問題最嚴(yán)重的區(qū)域之一.該地區(qū)氣溶膠組分非常復(fù)雜,因而針對ADRE的日內(nèi)變化的研究還很少.本文以長三角為研究區(qū),利用高時空分辨率的靜止衛(wèi)星GOCI反演得到的AOD數(shù)據(jù),結(jié)合地表反照率和地面觀測站數(shù)據(jù)進(jìn)行高污染地區(qū)日內(nèi)ADRE的研究.研究結(jié)果將對長三角區(qū)域環(huán)境的監(jiān)測、分析、治理等提供理論支持,并為政府制定相應(yīng)的環(huán)境法規(guī)提供參考.

    1 研究區(qū)與研究數(shù)據(jù)

    1.1 研究區(qū)概況

    以長江三角洲地區(qū)(28.04°~33.41°N,118.33°~ 122.43°E)為研究區(qū),包含江蘇的南京、揚(yáng)州、泰州、南通、鎮(zhèn)江、常州、無錫、蘇州8個城市,上海市,浙江的湖州、嘉興、杭州、紹興、寧波、舟山6個城市.長江三角洲地區(qū)北部以低海拔的平原為主,南部多為高山丘陵.地表覆蓋類型多樣,主要包括林地、草地、耕地、水域和建設(shè)用地等.

    長江三角洲地區(qū)是我國的主要經(jīng)濟(jì)區(qū)之一,能源的巨大消耗和機(jī)動車尾氣的大量排放,導(dǎo)致區(qū)域污染逐步加重,大氣中的各種氣體、細(xì)顆粒物質(zhì)和痕量有毒氣體構(gòu)成復(fù)合污染[7],造成區(qū)域能見度下降、霧霾天氣增多,危害居民和當(dāng)?shù)氐纳鷳B(tài)環(huán)境.

    1.2 研究數(shù)據(jù)

    1.2.1 GOCI數(shù)據(jù) 2009年,韓國發(fā)射了第一顆主要應(yīng)用于海洋監(jiān)測和氣象服務(wù)研究的地球靜止海洋氣象衛(wèi)星,GOCI為其搭載的星上水色傳感器,包含從可見光到近紅外共8個波段.其最明顯的優(yōu)勢是時間分辨率高達(dá)1h,每天可以提供目標(biāo)區(qū)域8個時相的觀測數(shù)據(jù),實(shí)時監(jiān)測朝鮮、俄羅斯、韓國、日本以及中國區(qū)域內(nèi)約為2500km′2500km的區(qū)域范圍[8].

    本文選取GOCI 2017年10~12月共計94景的晴空(云量覆蓋小于10%)圖像為研究數(shù)據(jù),其L1B數(shù)據(jù)與GOCI數(shù)據(jù)處理軟件可從韓國海洋衛(wèi)星中心網(wǎng)站下載.該官網(wǎng)上2017年的數(shù)據(jù)每天共有8景,個別天數(shù)由于網(wǎng)站限制只有7景.由于現(xiàn)有的AOD的反演算法主要適用于晴空圖像,因此篩選出共16d的晴空數(shù)據(jù)進(jìn)行研究.后經(jīng)試驗(yàn)發(fā)現(xiàn),GOCI的第一景和最后一景由于BRDF的影響,和本文假定地面為朗伯體的情況不符,因此只選擇其中的6景數(shù)據(jù)進(jìn)行反演.

    1.2.2 MODIS數(shù)據(jù) MODIS是搭載在Terra和Aqua衛(wèi)星上的傳感器,掃描寬度為2330km.該傳感器有可見光到熱紅外共36個光譜波段,空間分辨率有250m單日合成等5種形式[9],可以同時提供氣溶膠、云特性、植被指數(shù)、海洋特性數(shù)據(jù)等信息.MODIS具有較高的時空覆蓋度,且數(shù)據(jù)獲取非常容易,因而在我國得到了廣泛的應(yīng)用.

    GOCI的藍(lán)光波段和MODIS對應(yīng)波段范圍并不一致,本文利用MOD02輻亮度數(shù)據(jù)對GOCI數(shù)據(jù)進(jìn)行了交叉輻射校正.基于MOD09A1產(chǎn)品構(gòu)建地表反射率庫,以消除短期內(nèi)地表反射率的變化.為提高ADRE計算精度,選擇MCD43A3地表反照率產(chǎn)品代入6S(太陽光譜中衛(wèi)星信號的二次模擬)模型中進(jìn)行模擬.這2種均為SIN正弦曲線投影Grid格式的MODIS數(shù)據(jù)[10].

    1.2.3 CE-318數(shù)據(jù) NASA和法國國家科學(xué)研究中心構(gòu)建的AERONET(AErosol RObotic NETwork)是使用全自動型CE-318太陽分光光度計測量、精度較高的地基氣溶膠遙感監(jiān)測網(wǎng)絡(luò),可以提供長期的、連續(xù)的、免費(fèi)的AOD及其輻射特性的監(jiān)測,其AOD觀測誤差約為0.01~0.02[11].它可以自動追蹤太陽測量太陽直射輻射,具有容易攜帶、安裝方便的優(yōu)點(diǎn)[12].CE-318太陽光度計不僅廣泛應(yīng)用對于大氣氣溶膠和氣象氣候的監(jiān)測研究,還可以進(jìn)行大氣光學(xué)特性相關(guān)參數(shù)的推算.

    本文采用的驗(yàn)證數(shù)據(jù)為南京信息工程大學(xué)北辰樓站點(diǎn)(32.20°N,118.71°E)和AERONET太湖站(31.42°N,120.21°E)的地面實(shí)測AOD數(shù)據(jù),時間范圍和GOCI對應(yīng)一致.驗(yàn)證結(jié)果時需要將原始實(shí)測AOD值轉(zhuǎn)為550nm處的AOD[13].

    2 研究方法

    2.1 計算原理

    假設(shè)地表為朗伯體且大氣水平方向均一的條件下,入瞳處的輻射亮度TOA可表示為[14]:

    式中:v為衛(wèi)星天頂角的余弦;s為太陽天頂角的余弦;0為衛(wèi)星接收方向的路徑輻射;為相對方位角;0是大氣層頂?shù)奶柟廨椛渫棵芏?是大氣透過率;為地表反射率;為大氣整層向下路徑的半球反射率.、0和(s)(v)這3個參數(shù)含有大氣的信息,而包含地表反射率的信息.因此,要先剔除地表反射率的影響實(shí)現(xiàn)地氣解耦,進(jìn)一步反演得到AOD.

    ADRE的模擬是通過計算地表和大氣層分別在氣溶膠存在和不存在2種情況下短波輻射通量的差值來實(shí)現(xiàn).即在其他參數(shù)不變時,一次取反演的AOD值,一次取0值.公式為:

    式中:ΔTOA和ΔSUR分別為大氣層和地面處的ADRE;(CLR)TOA、(CLRSUR為無氣溶膠時大氣層和地面處的凈輻射通量密度;(AERO)TOA、(AERO)SUR為有氣溶膠時大氣層頂和地面處的凈輻射通量密度;↑、↓分別為上行和下行的輻射通量值,W/m2;為地表反照率[15].ADRE值的正負(fù)取決于氣溶膠粒子對太陽輻射吸收和其自身的反射能力等因素.上述AOD和ADRE模擬過程均是IDL編程并調(diào)用6S來實(shí)現(xiàn).

    2.2 6S模型反演流程

    6S可模擬傳感器在整個0.25~4um短波光譜通道內(nèi)接收到的輻射量.參考程晨[16]對長三角氣溶膠粒子組分的研究,將氣溶膠設(shè)為大陸型氣溶膠;通過GOCI原始數(shù)據(jù)獲得經(jīng)緯度、AOD、太陽天頂角和方位角、衛(wèi)星方位角和天頂角等逐像元信息;地面設(shè)為朗伯體;反照率為MCD43產(chǎn)品,其余為6S內(nèi)置參數(shù).技術(shù)流程如圖1所示.

    圖1 ADRE模擬流程

    3 結(jié)果與討論

    3.1 氣溶膠直接輻射效應(yīng)結(jié)果

    經(jīng)反演、計算后,得到研究區(qū)16d共94景AOD 和ADRE結(jié)果,其中12月24日典型重度霧霾天氣下的長三角地區(qū)的AOD、地表和大氣層頂?shù)腁DRE結(jié)果如下圖2~4所示,圖5為該日平均AOD和地表處ADRE的線性回歸分析.

    圖2 2017年12月24日氣溶膠光學(xué)厚度

    a.9:00;b.10:00;c.11:00;d.12:00;e.13:00;f.14:00

    a.9:00;b.10:00;c.11:00;d.12:00;e.13:00;f.14:00

    圖4 2017年12月24日ΔFTOA結(jié)果

    a.9:00;b.10:00;c.11:00;d.12:00;e.13:00;f.14:00

    圖5 2017年12月24日日均AOD和ΔFSUR的回歸分析

    圖3中,ΔSUR絕對值越大表明氣溶膠阻擋到達(dá)地表的太陽光越多,即氣溶膠濃度越大.圖4中, ΔTOA為負(fù)值說明由于氣溶膠的影響,更多的太陽光被反射回了地球外部空間,對地表產(chǎn)生了降溫效應(yīng).圖5可以表明AOD和其對應(yīng)的ΔSUR有非常顯著的線性相關(guān)關(guān)系,這與Wang等[4]的研究結(jié)果一致.ΔTOA的結(jié)果和Zhuang等[17]得出的晴空條件下日均ADRE為-6.9W/m2相似,驗(yàn)證了本文結(jié)果的可靠性.參照Li等[18]的研究,當(dāng)水和其他陸地表面反照率小于0.16時,AOD值和對應(yīng)的ΔTOA之間存在明顯的負(fù)相關(guān)關(guān)系,即ΔTOA隨著AOD的增加值越來越負(fù);反照率接近0.16時,盡管氣溶膠對大氣-地表輻射傳輸系統(tǒng)有微弱的影響,它以冷卻表面的代價來提高大氣層的溫度;當(dāng)反照率大于0.16時,AOD和ΔTOA之間有明顯的正相關(guān)關(guān)系.本文研究區(qū)內(nèi)地表的反照率值98%以上都小于0.16,因此AOD和ΔTOA之間存在負(fù)相關(guān)關(guān)系,和上述描述一致,并且與Gong等[19]研究中國典型城市和郊區(qū)得出的AOD和ΔTOA呈現(xiàn)明顯的負(fù)相關(guān)關(guān)系的結(jié)論一致,再一次驗(yàn)證了本文的結(jié)論.即圖2,4中,當(dāng)AOD值越高時,ΔTOA值越小即絕對值越大.綜上所述,由于氣溶膠的制冷效應(yīng),ΔTOA和ΔSUR隨著AOD值的增加而減小,即氣溶膠直接輻射強(qiáng)度增大.

    從ADRE日內(nèi)變化可以看出,ΔTOA變化特征不明顯,可能和地表反射率、6S中大陸型氣溶膠類型的設(shè)定有關(guān).

    3.2 精度評價

    本文選用研究區(qū)內(nèi)2017年10 ~12月共16d北辰樓和太湖2個站點(diǎn)衛(wèi)星過境時間段±0.5h內(nèi)的實(shí)測AOD數(shù)據(jù)與反演結(jié)果進(jìn)行時空匹配并進(jìn)行精度驗(yàn)證.由于站點(diǎn)數(shù)據(jù)不完整,最終北辰樓站點(diǎn)共有77組數(shù)據(jù),太湖站有51組數(shù)據(jù).結(jié)果如圖6所示,北辰樓站點(diǎn)的2=0.68,太湖站點(diǎn)2=0.67,且二者統(tǒng)計檢驗(yàn)值<0.001,均呈現(xiàn)了較高的擬合精度.

    進(jìn)一步分析誤差的來源,交叉輻射校正、地表反射率庫的精度,6S中大陸型氣溶膠和朗伯體地表的假定均對反演結(jié)果造成一定的誤差.由于缺乏相應(yīng)的實(shí)測資料,6S中這些參數(shù)的假定仍然是目前較為通用的且為大家所接受的一種方式.此外,地基觀測資料和衛(wèi)星輻射觀測本身的誤差都會對反演結(jié)果造成一定的影響.由于長三角地區(qū)ADRE是基于該地區(qū)AOD反演結(jié)果的二次模擬產(chǎn)品,且周杰等[13]、王文君等[15]的研究已驗(yàn)證過該算法的可靠性,因此AOD結(jié)果精度良好一定程度上表明ADRE模擬結(jié)果的可靠性.

    3.3 不同地物對ADRE的影響分析

    本文進(jìn)一步選取不受氣流影響的典型霾天12月20日并對該日6個時刻不同地物類型的ΔTOA值進(jìn)行統(tǒng)計分析,結(jié)果如表1所示.就水體而言,其吸收較強(qiáng)而反射較弱,從而反照率較低,因此6景圖像中均是水體的最小值和平均值最低.相反地,高反射率地表類型如城市的ADRE值較高,最大值均為正值,與Chen[20]等的研究結(jié)果一致.水體對應(yīng)的ΔTOA值一天內(nèi)變化最大,可見地表反照率對ΔTOA有顯著影響,這與劉建軍[21]、Zhuang等[17]得出的結(jié)論一致.森林和農(nóng)作物均屬于植被類型,因此同一時刻二者差異不大.

    表1 2017年12月20日6景不同地物類型ΔFTOA值統(tǒng)計表(W/m2)

    3.4 ADRE的日內(nèi)變化規(guī)律

    為探究ADRE的日內(nèi)變化特征,本文計算了20日6景圖像的平均值,并比較分析每個時刻和平均值的差異,結(jié)果如下圖所示.本文中距平指每個時刻的結(jié)果和當(dāng)日6個時刻均值之差.

    圖7 2017年12月20日ΔFSUR均值分布

    從表2得知,每一景平均值和標(biāo)準(zhǔn)偏差相比最值有較大偏差.MODIS過境時間約為當(dāng)?shù)?0:30左右,而10:00和11:00影像的平均值和標(biāo)準(zhǔn)偏差最小,即這個時間段可以較好的模擬日均ADRE值,與Kassianov等[5]認(rèn)為的MODIS氣溶膠產(chǎn)品可以合理的估計日均氣溶膠光學(xué)特性和ADRE的結(jié)果一致,證實(shí)了本文反演結(jié)果的可靠性.Arola[3]也認(rèn)為在沒有AOD日變化數(shù)據(jù)的基礎(chǔ)上,可近似利用MODIS的AOD數(shù)據(jù)計算ADRE來表征當(dāng)天的平均情況.

    從圖3、4可以看出,ADRE值的變化趨勢表征了當(dāng)日長三角地區(qū)氣溶膠受氣流影響的分布狀況,并成功捕捉到了一次該地區(qū)冬季霧霾爆發(fā)并消失的過程.從時間分布上來看,09:00影像受到云區(qū)的影響存在部分缺值,之后該云區(qū)消失.該影像上長三角西北部地區(qū)出現(xiàn)了較廣范圍的濃度極高的氣溶膠,并逐漸向東南方向移動.10:00左右,氣溶膠范圍擴(kuò)散到了整個長三角區(qū)域,范圍很大.11:00左右,氣溶膠濃度明顯降低,這可能是受到風(fēng)向及氣溶膠沉降的影響.早上08:00,09:00和下午17:00,18:00是上、下班高峰期,氣溶膠濃度相對較高且氣溶膠輻射效應(yīng)強(qiáng)度也高;中午時分熱力強(qiáng)迫和垂直混合作用的加強(qiáng)導(dǎo)致氣溶膠容易擴(kuò)散且濃度降低.此時邊界層高度升高,直到下班高峰期這種現(xiàn)象才會消失,氣溶膠濃度也將再次升高,這對氣溶膠的時間分布趨勢也有一定的影響.12:00左右,氣溶膠逐漸擴(kuò)散到江蘇南部、上海等區(qū)域,范圍和濃度均進(jìn)一步減小,最終該氣溶膠團(tuán)往南部擴(kuò)散并消失.張悅等[22]通過設(shè)置3個不同氣溶膠濃度的試驗(yàn)區(qū)分其對輻射效應(yīng)的影響,發(fā)現(xiàn)不論是氣溶膠直接、半直接還是間接輻射效應(yīng)均使得污染區(qū)短波輻射減少、2m以內(nèi)氣溫下降、大氣邊界層高度降低,不利于水汽與污染物的擴(kuò)散,從而造成空氣污染進(jìn)一步加重;霧結(jié)構(gòu)進(jìn)一步穩(wěn)定,并延長了霧的持續(xù)時間,發(fā)展高度更高.可見ADRE與大氣污染息息相關(guān),能進(jìn)一步加重大氣污染程度,并利于霧的發(fā)生和發(fā)展.綜上所述,由于太陽輻射是地球熱力系統(tǒng)的直接來源,因而局地氣溫、風(fēng)速和邊界層高度等氣象要素都受到氣溶膠輻射效應(yīng)的影響,分析ADRE的日內(nèi)變化對氣象部門監(jiān)測近地表氣溫、分析灰霾的形成到消失過程、研究邊界層高度等氣象要素具有重要的意義,并對政府部門制定相關(guān)政策提供幫助.

    圖8 2017年12月20日9:00~14:00 ΔFSUR距平分布

    a.9:00;b.10:00;c.11:00;d.12:00;e.13:00;f.14:00

    表2 2017年12月20日ΔFSUR距平統(tǒng)計結(jié)果

    4 結(jié)論

    4.1 2017年10~12月長江三角洲地區(qū)GOCI數(shù)據(jù)反演的AOD值和地面實(shí)測資料具有很好的擬合精度,其中北辰樓2=0.68,太湖2=0.67.

    4.2 分析ADRE時空分布和典型的灰霾天氣發(fā)現(xiàn),氣溶膠存在明顯的制冷效應(yīng).AOD和其對應(yīng)的大氣層頂及地表ADRE強(qiáng)度有顯著的線性關(guān)系,高值和低值分布的區(qū)域大致相同.

    4.3 分析不同地物類型的ΔTOA值發(fā)現(xiàn),由于水體的地表反照率比較低,因而其ΔTOA最低且一天內(nèi)變化最大,可見地表反照率對ΔTOA有顯著影響.相反地,高反射率地表類型的ΔTOA較高,甚至出現(xiàn)了正值,這與前人的結(jié)果一致.

    4.4 GOCI 10:00和11:00影像可以較好地模擬日均ADRE.從時間分布上來看,早上和下午的氣溶膠輻射效應(yīng)強(qiáng)度最大、中午最小,這是受到風(fēng)向以及氣溶膠沉降的影響.ADRE的日變特征能成功捕捉一次霧霾爆發(fā)并消失的過程,對氣象部門監(jiān)測近地表氣溫、分析霧的形成和消失的過程及政府部門制定相關(guān)政策提供幫助.

    [1] Shao P, Xin J Y, An J L, et al. The empirical relationship between PM2.5and AOD in Nanjing of the Yangtze River Delta [J]. Atmospheric Pollution Research, 2017,8:233-243.

    [2] 宋 姚,麻金繼.利用MODIS數(shù)據(jù)計算北京地區(qū)氣溶膠直接輻射效應(yīng)[J]. 大氣與環(huán)境光學(xué)學(xué)報, 2012,7(6):1673-6141.

    [3] Arola A, Eck T F, Huttunen J, et al. Influence of observed diurnal cycles of aerosol optical depth on aerosol direct radiative effect [J]. Atmospheric Chemistry and Physics, 2013,13:7895–7901.

    [4] Wang Z, Liu D, Wang Y, et al. Diurnal aerosol variations do affect daily averaged radiative forcing under heavy aerosol loading observed in Hefei, China [J]. Atmospheric Measurement Techniques, 2015,8: 2901–2907.

    [5] Kassianov E, Barnard J, Pekour M, et al. Do Diurnal aerosol changes affect daily average radiative forcing? [J]. Geophysical Research Letters, 2013,40:3265–3269.

    [6] 段 婧,毛節(jié)泰.長江三角洲大氣氣溶膠光學(xué)厚度分布和變化趨勢研究[J]. 環(huán)境科學(xué)學(xué)報, 2007,27(40):537-543.

    [7] 翟 華,朱 彬,趙雪婷,等.長江三角洲初冬一次重污染天氣成因分析[J]. 中國環(huán)境科學(xué), 2018,38(11):4001-4009.

    [8] 李冠男,王 林,王 祥,等.靜止水色衛(wèi)星GOCI及其應(yīng)用進(jìn)展[J]. 海洋環(huán)境科學(xué), 2014,33(6):966-971.

    [9] 王中挺,厲 青,王 橋,等.利用深藍(lán)算法從HJ-1數(shù)據(jù)反演陸地氣溶膠[J]. 遙感學(xué)報, 2012,16(3):596-610.

    [10] 張 璐.基于國產(chǎn)遙感衛(wèi)星數(shù)據(jù)的北京市氣溶膠光學(xué)厚度反演研究[D]. 上海:華東師范大學(xué), 2016.

    [11] 張志薇,王宏斌,張 鐳,等.中國地區(qū)3個AERONET站點(diǎn)氣溶膠直接輻射強(qiáng)迫分析[J]. 中國科學(xué)院大學(xué)學(xué)報, 2014,31(3):297-305.

    [12] 李 禮,余家燕,楊 燦,等.CE-318太陽光度計在大氣環(huán)境監(jiān)測中的應(yīng)用[J]. 環(huán)境科學(xué)與管理, 2012,37(2):107-110.

    [13] 周 杰,陳 健,張海龍,等.基于MODIS數(shù)據(jù)的長三角地區(qū)氣溶膠直接輻射效應(yīng)研究[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2018,34(1):54-63.

    [14] Kaufman Y J, Tanri D, Remer L A, et al. Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer [J]. Journal of Geophysical Research, 1997,102(14):17051-17067.

    [15] 王文君.基于衛(wèi)星和地面觀測的長三角地區(qū)氣溶膠直接輻射效應(yīng)研究[D]. 南京:南京信息工程大學(xué), 2016.

    [16] 程 晨.基于HJ-1衛(wèi)星的南京氣溶膠光學(xué)厚度反演和時空變化分析[D]. 南京:南京信息工程大學(xué), 2013.

    [17] Zhuang B L, Wang T J, Li S et al. Optical properties and radiative forcing of urban aerosols in Nanjing, China [J]. Atmospheric Environment, 2014,83:43–52.

    [18] Li Z Q, Lee K H, Wang Y S, et al. First observation-based estimates of cloud-free aerosol radiative forcing across China [J]. Journal of Geophysical Research -Atmospheres, 2010,115(7):1029-1037.

    [19] Gong C S, Xin J Y, Wang S G, et al. Anthropogenic aerosol optical and radiative properties in the typical urban/suburban regions in China [J]. Atmospheric Research, 2017,197:177-187.

    [20] Chen J, Zhou J, Wang W J, et al. Satellite observations of the aerosol direct radiative effect during a pollution episode in Nanjing, China, using high-resolution HJ-1CCD imagery [J]. International Journal of Remote Sensing, 2017,38(16):4535-4552.

    [21] 劉建軍.長三角太湖地區(qū)云和氣溶膠輻射特性的地基遙感研究 [D]. 南京:南京信息工程大學(xué), 2012.

    [22] 張 悅,藩曙先,李 皓,等.氣溶膠輻射效應(yīng)在華東地區(qū)一次霧霾過程中的作用[J]. 氣象學(xué)報, 2016,74(3):465-478.

    致謝:本研究采用了韓國GOCI官網(wǎng)、Aeronet官網(wǎng)、南信大北辰樓地面觀測站基地提供的數(shù)據(jù),在此表示衷心的感謝!同時感謝審稿人對本文提出的寶貴的意見!

    Diurnal variations of aerosol radiative effect under haze weather condition using GOCI data-A case study of Yangtze River Delta.

    LI Ya-wen1, CHEN Jian1*, ZHANG Hai-long2, JIN Shuang-gen1

    (1.School of Remote Sensing & Geomatics Engineering, Nanjing University of Information Science & Technology, Nanjing 210044, China;2.School of Marine Science, Nanjing University of Information Science & Technology, Nanjing 210044, China)., 2019,39(2):497~505

    Based on the Second Simulation of the Satellite Signal in the Solar Spectrum (6S) model and the Geostationary Ocean Color Imager (GOCI) L1B data of the clear-sky geostationary satellites with a total of 94 images observed in 16 days selected from October to December in 2017, the Aerosol Optical Depth (AOD) in the Yangtze River Delta was inversed with the deep blue algorithm, which was validated by in-suit observation data at two stations. The Aerosol Direct Radiation Effect (ADRE) on the surface and the top layer of the atmosphere were further calculated and analyzed under typical haze weather. The results showed that the AOD retrieved from GOCI data had high fitting precision. In addition, the fitting correlation2at Beichen building site and Taihu site were 0.68 and 0.67, respectively. Due to the refrigeration effect of aerosols, the AOD had a significant linear relationship with the direct radiant intensity of aerosols on the surface and at the top of the atmosphere spatially. Due to aerosol diffusion and wind direction, the intensity of radiation effects was high in the morning and afternoon and low at noon temporally. The images observed at 10:00 LT and 11:00 LT could better simulate the characteristics of the daily average ADRE. The process of an outburst and disappearance of the haze was successfully captured by using the ADRE diurnal variations, which was greatly significant to measure the near-surface temperature and analyze the formation of haze for the meteorological agency.

    haze;aerosol;direct radiation effect;the Yangtze River Delta;GOCI

    X513

    A

    1000-6923(2019)02-0497-09

    李雅雯(1995-),女,江蘇連云港人.南京信息工程大學(xué)碩士研究生,主要從事定量遙感方面的研究.發(fā)表論文3篇.

    2018-07-04

    國家重點(diǎn)研發(fā)計劃(2018YFC1506404)

    * 責(zé)任作者, 副教授, chjnjnu@163.com

    猜你喜歡
    反照率氣溶膠反演
    基于藍(lán)天空反照率的氣溶膠輻射強(qiáng)迫模擬
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    薩吾爾山木斯島冰川反照率時空變化特征研究
    冰川凍土(2022年6期)2022-02-12 08:31:06
    長江三角洲地區(qū)大氣氣溶膠柱單次散射反照率特性研究
    氣溶膠傳播之謎
    氣溶膠中210Po測定的不確定度評定
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    四川盆地秋季氣溶膠與云的相關(guān)分析
    大氣氣溶膠成核監(jiān)測
    美女午夜性视频免费| 制服人妻中文乱码| 在线免费观看的www视频| 国产一区二区三区视频了| 亚洲av中文字字幕乱码综合| 免费在线观看成人毛片| 日韩欧美免费精品| 2021天堂中文幕一二区在线观| 好看av亚洲va欧美ⅴa在| 色av中文字幕| 99热6这里只有精品| 精品欧美国产一区二区三| 麻豆av在线久日| 国产午夜精品久久久久久| 免费av毛片视频| 少妇人妻一区二区三区视频| 国产爱豆传媒在线观看 | 久久香蕉国产精品| 婷婷六月久久综合丁香| 亚洲美女黄片视频| 日韩欧美 国产精品| 啦啦啦韩国在线观看视频| 可以免费在线观看a视频的电影网站| 男插女下体视频免费在线播放| 舔av片在线| 精品国产美女av久久久久小说| 五月伊人婷婷丁香| 国产99白浆流出| 老熟妇乱子伦视频在线观看| 国产免费男女视频| 日韩大尺度精品在线看网址| 免费看日本二区| 成人av在线播放网站| 最近最新免费中文字幕在线| 舔av片在线| 亚洲狠狠婷婷综合久久图片| 国产三级在线视频| 国产一级毛片七仙女欲春2| 午夜免费激情av| 中文亚洲av片在线观看爽| 男女床上黄色一级片免费看| 国产一区二区在线观看日韩 | 日本五十路高清| 小说图片视频综合网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久免费视频了| 欧美中文综合在线视频| 又爽又黄无遮挡网站| 特大巨黑吊av在线直播| 国产蜜桃级精品一区二区三区| 久久久久久久午夜电影| 久久欧美精品欧美久久欧美| 在线观看66精品国产| 又黄又粗又硬又大视频| 亚洲人成网站高清观看| 亚洲第一电影网av| 午夜老司机福利片| 中国美女看黄片| 久久中文字幕人妻熟女| 波多野结衣高清作品| 草草在线视频免费看| 中文字幕人成人乱码亚洲影| 午夜福利在线在线| 亚洲中文字幕日韩| 欧美性长视频在线观看| 欧美乱码精品一区二区三区| 精品国产乱子伦一区二区三区| 欧美成人免费av一区二区三区| 久久草成人影院| 别揉我奶头~嗯~啊~动态视频| 人成视频在线观看免费观看| 久久久国产精品麻豆| 免费高清视频大片| 亚洲av美国av| 波多野结衣高清作品| 床上黄色一级片| 国产精品免费一区二区三区在线| 欧美黄色片欧美黄色片| 久久中文看片网| 日本精品一区二区三区蜜桃| 精品乱码久久久久久99久播| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 啦啦啦免费观看视频1| 97人妻精品一区二区三区麻豆| 2021天堂中文幕一二区在线观| www.自偷自拍.com| 国产精品久久久久久亚洲av鲁大| 嫩草影视91久久| 亚洲全国av大片| 在线观看www视频免费| 日本黄色视频三级网站网址| 国产精品一区二区精品视频观看| 激情在线观看视频在线高清| 少妇熟女aⅴ在线视频| 日本熟妇午夜| 国产午夜精品论理片| 一本大道久久a久久精品| netflix在线观看网站| 日韩欧美在线乱码| 午夜日韩欧美国产| 久久久久久久久中文| 久久久久久久久中文| 日韩欧美在线乱码| 日韩欧美在线乱码| 国产精品亚洲美女久久久| 欧美 亚洲 国产 日韩一| 欧美精品啪啪一区二区三区| 亚洲成av人片在线播放无| 精品高清国产在线一区| 波多野结衣巨乳人妻| 中文字幕高清在线视频| 波多野结衣巨乳人妻| 韩国av一区二区三区四区| 香蕉久久夜色| 欧美高清成人免费视频www| 欧美 亚洲 国产 日韩一| 国产成+人综合+亚洲专区| 搞女人的毛片| 又紧又爽又黄一区二区| 日韩欧美 国产精品| 大型黄色视频在线免费观看| 日韩欧美精品v在线| 亚洲一区二区三区色噜噜| 国产爱豆传媒在线观看 | 久久亚洲精品不卡| 久久中文看片网| 一夜夜www| 这个男人来自地球电影免费观看| 中国美女看黄片| 久久亚洲真实| 国产精品久久视频播放| 大型av网站在线播放| 亚洲国产中文字幕在线视频| 成年女人毛片免费观看观看9| 色尼玛亚洲综合影院| 黄色丝袜av网址大全| 怎么达到女性高潮| 搡老熟女国产l中国老女人| 韩国av一区二区三区四区| 亚洲最大成人中文| 欧美不卡视频在线免费观看 | 欧美人与性动交α欧美精品济南到| 成年免费大片在线观看| 亚洲第一欧美日韩一区二区三区| 成熟少妇高潮喷水视频| www.www免费av| 精品国产乱码久久久久久男人| 51午夜福利影视在线观看| 69av精品久久久久久| 亚洲黑人精品在线| 亚洲黑人精品在线| 亚洲人与动物交配视频| 成人18禁在线播放| 最新美女视频免费是黄的| 婷婷丁香在线五月| 可以免费在线观看a视频的电影网站| 亚洲欧美激情综合另类| 精品无人区乱码1区二区| 中文字幕精品亚洲无线码一区| 1024手机看黄色片| 国产精品 国内视频| 中亚洲国语对白在线视频| 久99久视频精品免费| 国产亚洲av嫩草精品影院| 国产日本99.免费观看| 黑人欧美特级aaaaaa片| 亚洲欧美日韩东京热| 欧美日本视频| 亚洲国产欧美网| 亚洲欧美日韩东京热| 又黄又爽又免费观看的视频| 欧美日韩瑟瑟在线播放| 国产伦在线观看视频一区| 好看av亚洲va欧美ⅴa在| 中文资源天堂在线| 国产精品久久久久久久电影 | 老汉色av国产亚洲站长工具| 国产精品爽爽va在线观看网站| 夜夜爽天天搞| 别揉我奶头~嗯~啊~动态视频| 99久久综合精品五月天人人| 脱女人内裤的视频| 婷婷精品国产亚洲av| 伊人久久大香线蕉亚洲五| 床上黄色一级片| 又大又爽又粗| 久久久久精品国产欧美久久久| 日本成人三级电影网站| 国产精品久久久人人做人人爽| 一级片免费观看大全| 久久伊人香网站| 黑人巨大精品欧美一区二区mp4| 黄色女人牲交| 嫁个100分男人电影在线观看| 观看免费一级毛片| 国产亚洲av高清不卡| 在线免费观看的www视频| 国产亚洲av嫩草精品影院| 中文字幕久久专区| 老熟妇仑乱视频hdxx| 麻豆国产97在线/欧美 | 亚洲avbb在线观看| 欧美日韩国产亚洲二区| www.自偷自拍.com| 亚洲精品美女久久久久99蜜臀| 亚洲九九香蕉| 国产视频内射| 脱女人内裤的视频| 亚洲专区国产一区二区| 黄色女人牲交| 亚洲黑人精品在线| 女人被狂操c到高潮| 一本大道久久a久久精品| 神马国产精品三级电影在线观看 | 观看免费一级毛片| 在线观看免费视频日本深夜| 丁香欧美五月| 精品免费久久久久久久清纯| 非洲黑人性xxxx精品又粗又长| 国产精品永久免费网站| 亚洲精品一卡2卡三卡4卡5卡| 久久精品aⅴ一区二区三区四区| 久久久精品欧美日韩精品| 国产成人av教育| 久久 成人 亚洲| 每晚都被弄得嗷嗷叫到高潮| 久久精品影院6| 亚洲精品久久国产高清桃花| 伦理电影免费视频| 成年版毛片免费区| av在线播放免费不卡| 亚洲精品中文字幕一二三四区| x7x7x7水蜜桃| 黄色a级毛片大全视频| 一边摸一边做爽爽视频免费| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 欧美一级毛片孕妇| 香蕉久久夜色| 国产三级中文精品| 亚洲成人中文字幕在线播放| 91麻豆精品激情在线观看国产| 国产精品九九99| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品av麻豆狂野| 成人精品一区二区免费| 日韩精品中文字幕看吧| 成人高潮视频无遮挡免费网站| 全区人妻精品视频| 大型黄色视频在线免费观看| 丁香欧美五月| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 搡老妇女老女人老熟妇| 欧美成人一区二区免费高清观看 | 色精品久久人妻99蜜桃| 亚洲免费av在线视频| 色av中文字幕| 两人在一起打扑克的视频| 亚洲第一电影网av| 一区二区三区国产精品乱码| 成人永久免费在线观看视频| 亚洲中文字幕日韩| 亚洲真实伦在线观看| 国产免费av片在线观看野外av| АⅤ资源中文在线天堂| 日本a在线网址| 桃红色精品国产亚洲av| 老司机福利观看| 看免费av毛片| 亚洲一区二区三区不卡视频| 国产精品1区2区在线观看.| 精品乱码久久久久久99久播| 精品熟女少妇八av免费久了| 50天的宝宝边吃奶边哭怎么回事| 999久久久国产精品视频| 美女免费视频网站| 国产三级黄色录像| 狠狠狠狠99中文字幕| 国产欧美日韩一区二区精品| 亚洲一区高清亚洲精品| 夜夜躁狠狠躁天天躁| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲 欧美一区二区三区| 岛国在线观看网站| 亚洲成人久久性| 美女大奶头视频| 国产一区二区激情短视频| 一级片免费观看大全| 美女 人体艺术 gogo| 成熟少妇高潮喷水视频| 一级毛片精品| 老司机深夜福利视频在线观看| 国产精品99久久99久久久不卡| 亚洲av成人不卡在线观看播放网| 国产主播在线观看一区二区| 国产精品自产拍在线观看55亚洲| 欧美乱妇无乱码| 欧美中文综合在线视频| 一本综合久久免费| 欧美日韩黄片免| 国产免费av片在线观看野外av| 黄频高清免费视频| 三级毛片av免费| 国产成人av激情在线播放| 国内精品一区二区在线观看| 午夜免费成人在线视频| 一二三四社区在线视频社区8| 中国美女看黄片| 成人精品一区二区免费| 在线视频色国产色| 欧美成人性av电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 最好的美女福利视频网| av在线天堂中文字幕| 特大巨黑吊av在线直播| 全区人妻精品视频| 少妇被粗大的猛进出69影院| 一夜夜www| 国产日本99.免费观看| 啦啦啦韩国在线观看视频| 欧美成人午夜精品| 欧美在线黄色| 免费观看人在逋| 国内精品久久久久久久电影| 免费观看精品视频网站| 亚洲精品一卡2卡三卡4卡5卡| 久久欧美精品欧美久久欧美| 欧美黑人欧美精品刺激| 老鸭窝网址在线观看| e午夜精品久久久久久久| 精品国产乱子伦一区二区三区| 黄色成人免费大全| 美女免费视频网站| 国产一区二区激情短视频| 男插女下体视频免费在线播放| 人妻久久中文字幕网| 一区二区三区国产精品乱码| 一区二区三区激情视频| 巨乳人妻的诱惑在线观看| 欧美一区二区国产精品久久精品 | 国产区一区二久久| 看片在线看免费视频| 国产一级毛片七仙女欲春2| 欧美黑人精品巨大| 国内毛片毛片毛片毛片毛片| 久99久视频精品免费| 免费在线观看日本一区| 亚洲电影在线观看av| 他把我摸到了高潮在线观看| 国产亚洲精品久久久久5区| x7x7x7水蜜桃| 国产激情欧美一区二区| 久热爱精品视频在线9| 欧美黑人精品巨大| 色播亚洲综合网| 国产v大片淫在线免费观看| 特级一级黄色大片| 久久久久亚洲av毛片大全| 国产精品一区二区三区四区久久| 国产精品美女特级片免费视频播放器 | 国产一区二区在线av高清观看| 成人手机av| 观看免费一级毛片| 国产视频一区二区在线看| 他把我摸到了高潮在线观看| 欧美黑人巨大hd| 欧美另类亚洲清纯唯美| 国产精品乱码一区二三区的特点| 神马国产精品三级电影在线观看 | 久久人妻av系列| 久久香蕉精品热| 国产精品一及| 午夜福利在线观看吧| 夜夜躁狠狠躁天天躁| 老司机靠b影院| 国产精品1区2区在线观看.| 亚洲精品一区av在线观看| 国产精品日韩av在线免费观看| 别揉我奶头~嗯~啊~动态视频| 久久久精品大字幕| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久久久99蜜臀| 午夜免费激情av| 观看免费一级毛片| 一区二区三区激情视频| 深夜精品福利| 国产欧美日韩一区二区精品| 精品久久久久久久久久免费视频| 美女免费视频网站| 久久久国产成人免费| 不卡av一区二区三区| 久久中文看片网| 久9热在线精品视频| 欧美久久黑人一区二区| 97碰自拍视频| 日本一二三区视频观看| 国产精品国产高清国产av| АⅤ资源中文在线天堂| 久久久久国产一级毛片高清牌| 久久精品91蜜桃| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 欧美一区二区精品小视频在线| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频| 久久婷婷人人爽人人干人人爱| 这个男人来自地球电影免费观看| 国产aⅴ精品一区二区三区波| 亚洲电影在线观看av| 国产高清视频在线观看网站| 别揉我奶头~嗯~啊~动态视频| 一级毛片高清免费大全| 精品人妻1区二区| 性色av乱码一区二区三区2| 成人三级做爰电影| 一二三四在线观看免费中文在| 香蕉av资源在线| 成人av在线播放网站| 欧美绝顶高潮抽搐喷水| 亚洲片人在线观看| 青草久久国产| 欧美日本视频| 男人舔女人下体高潮全视频| 久久久久九九精品影院| 久久 成人 亚洲| www日本黄色视频网| 免费人成视频x8x8入口观看| 中文字幕最新亚洲高清| 午夜视频精品福利| 亚洲av第一区精品v没综合| 国产av一区在线观看免费| 老汉色∧v一级毛片| 欧美一区二区国产精品久久精品 | 男女那种视频在线观看| 国产爱豆传媒在线观看 | 免费一级毛片在线播放高清视频| 悠悠久久av| av在线天堂中文字幕| 国产成人啪精品午夜网站| 精品国产亚洲在线| 淫秽高清视频在线观看| 可以免费在线观看a视频的电影网站| 精品福利观看| 十八禁网站免费在线| 亚洲人成伊人成综合网2020| 国内精品一区二区在线观看| 亚洲成人国产一区在线观看| 日本三级黄在线观看| 777久久人妻少妇嫩草av网站| 宅男免费午夜| 成人亚洲精品av一区二区| 熟女电影av网| 黄色视频不卡| 国产亚洲欧美98| 欧美不卡视频在线免费观看 | 少妇的丰满在线观看| 国产精品免费视频内射| 日韩欧美三级三区| 麻豆国产av国片精品| 看黄色毛片网站| 国产在线观看jvid| a级毛片在线看网站| 国产又黄又爽又无遮挡在线| 999精品在线视频| 久久久久久久午夜电影| aaaaa片日本免费| 欧美精品啪啪一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲精品久久国产高清桃花| 亚洲欧美日韩高清在线视频| 亚洲av成人不卡在线观看播放网| 免费人成视频x8x8入口观看| 天堂av国产一区二区熟女人妻 | 午夜免费观看网址| 国产欧美日韩一区二区精品| 看黄色毛片网站| 亚洲国产精品sss在线观看| 久久亚洲真实| 男女那种视频在线观看| 男人舔女人的私密视频| 制服丝袜大香蕉在线| √禁漫天堂资源中文www| 国产精品久久久久久久电影 | 精品无人区乱码1区二区| 欧美日韩国产亚洲二区| 久久精品夜夜夜夜夜久久蜜豆 | 免费在线观看成人毛片| 日本熟妇午夜| 国产高清视频在线观看网站| 很黄的视频免费| 激情在线观看视频在线高清| 国产精品久久久人人做人人爽| 麻豆成人午夜福利视频| 好看av亚洲va欧美ⅴa在| 亚洲片人在线观看| 国产av又大| 亚洲av五月六月丁香网| 99国产综合亚洲精品| 波多野结衣高清作品| 俺也久久电影网| 欧美高清成人免费视频www| 亚洲成人中文字幕在线播放| 精品久久久久久久末码| 日日爽夜夜爽网站| 国产精品影院久久| 国产高清有码在线观看视频 | 看免费av毛片| 亚洲乱码一区二区免费版| 亚洲av熟女| 男女那种视频在线观看| 老汉色∧v一级毛片| 九九热线精品视视频播放| 亚洲av电影不卡..在线观看| 亚洲精品中文字幕在线视频| 久久久精品欧美日韩精品| 午夜成年电影在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 久久热在线av| 国产男靠女视频免费网站| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 国产精品久久久av美女十八| 午夜免费激情av| 女人爽到高潮嗷嗷叫在线视频| 午夜免费观看网址| 国产亚洲欧美98| 国产一级毛片七仙女欲春2| 亚洲欧美一区二区三区黑人| 国产探花在线观看一区二区| av在线播放免费不卡| 免费在线观看黄色视频的| 欧美精品亚洲一区二区| 天天添夜夜摸| 亚洲在线自拍视频| 日韩成人在线观看一区二区三区| 午夜免费激情av| 在线十欧美十亚洲十日本专区| 亚洲精品美女久久av网站| 国内精品久久久久精免费| 久久婷婷人人爽人人干人人爱| 亚洲成a人片在线一区二区| 国产精品综合久久久久久久免费| 久久午夜亚洲精品久久| 亚洲欧美精品综合一区二区三区| 国产伦人伦偷精品视频| 午夜日韩欧美国产| АⅤ资源中文在线天堂| 精品第一国产精品| 午夜免费观看网址| √禁漫天堂资源中文www| 午夜福利在线在线| 可以在线观看的亚洲视频| 久久久久久免费高清国产稀缺| 男女视频在线观看网站免费 | 欧美黑人欧美精品刺激| 美女扒开内裤让男人捅视频| 午夜福利18| 午夜老司机福利片| 人妻夜夜爽99麻豆av| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 淫秽高清视频在线观看| 午夜精品一区二区三区免费看| 丰满人妻一区二区三区视频av | 欧美乱码精品一区二区三区| 久久精品国产亚洲av香蕉五月| 人妻久久中文字幕网| 亚洲av成人不卡在线观看播放网| 搡老妇女老女人老熟妇| 亚洲av成人一区二区三| 久久久久久久久久黄片| 国产高清videossex| 国产av一区二区精品久久| 精品日产1卡2卡| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 国产真人三级小视频在线观看| 国产单亲对白刺激| 男女下面进入的视频免费午夜| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 午夜久久久久精精品| www.999成人在线观看| 中文字幕久久专区| 久久精品综合一区二区三区| 国产av麻豆久久久久久久| or卡值多少钱| 亚洲成av人片在线播放无| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 国产精品av久久久久免费| 91九色精品人成在线观看| 九九热线精品视视频播放| 精品高清国产在线一区| 九九热线精品视视频播放| 老司机在亚洲福利影院| 亚洲狠狠婷婷综合久久图片| 麻豆一二三区av精品| 国产一区二区激情短视频| 一区二区三区国产精品乱码| 欧美日韩中文字幕国产精品一区二区三区| 国产探花在线观看一区二区| 九色成人免费人妻av| 看免费av毛片| 国产午夜精品久久久久久| 国产精品野战在线观看| 最近视频中文字幕2019在线8| 欧美乱码精品一区二区三区| 一夜夜www| 精品福利观看| 成年人黄色毛片网站| 身体一侧抽搐| 国产成人一区二区三区免费视频网站|