• <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)測
    妹子高潮喷水视频| 黑丝袜美女国产一区| 亚洲国产精品一区三区| 中文精品一卡2卡3卡4更新| 最近最新中文字幕免费大全7| 亚洲伊人久久精品综合| 日本-黄色视频高清免费观看| 国产成人精品一,二区| 久久午夜福利片| 成人黄色视频免费在线看| 亚洲在久久综合| 老熟女久久久| 我要看黄色一级片免费的| 一级毛片 在线播放| 中文天堂在线官网| 国产成人午夜福利电影在线观看| 国产亚洲一区二区精品| 久久久久久久国产电影| 国产又爽黄色视频| 99国产综合亚洲精品| 国产精品一国产av| 中国国产av一级| 热99国产精品久久久久久7| 国产69精品久久久久777片| 久久久国产欧美日韩av| 黄色配什么色好看| 日韩av不卡免费在线播放| 国产精品秋霞免费鲁丝片| 国产视频首页在线观看| 最近手机中文字幕大全| 亚洲精品一区蜜桃| 亚洲,欧美,日韩| 寂寞人妻少妇视频99o| 哪个播放器可以免费观看大片| 国产成人免费观看mmmm| 中文字幕av电影在线播放| 国产 精品1| 中文字幕免费在线视频6| 日韩精品免费视频一区二区三区 | 午夜免费男女啪啪视频观看| 久久久久久久久久人人人人人人| 国产欧美日韩综合在线一区二区| 国产成人av激情在线播放| 高清欧美精品videossex| 青春草国产在线视频| 免费少妇av软件| 色网站视频免费| 狂野欧美激情性bbbbbb| 日韩制服丝袜自拍偷拍| 欧美激情极品国产一区二区三区 | 纯流量卡能插随身wifi吗| 亚洲国产精品999| 久久鲁丝午夜福利片| 亚洲av国产av综合av卡| 国产精品一区二区在线不卡| 久久精品夜色国产| 亚洲国产精品999| 1024视频免费在线观看| 伦理电影免费视频| 国产成人aa在线观看| 欧美97在线视频| 免费观看av网站的网址| 国产欧美另类精品又又久久亚洲欧美| 中文字幕精品免费在线观看视频 | 国产国拍精品亚洲av在线观看| 大香蕉97超碰在线| 久久婷婷青草| 只有这里有精品99| 国产高清国产精品国产三级| 亚洲成人一二三区av| 日韩成人伦理影院| av免费观看日本| 日韩成人av中文字幕在线观看| 99热全是精品| 亚洲国产日韩一区二区| www.色视频.com| 精品卡一卡二卡四卡免费| 免费黄网站久久成人精品| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄| 丝袜美足系列| 一个人免费看片子| av网站免费在线观看视频| 国产男女内射视频| 性色av一级| 亚洲,欧美精品.| 亚洲精品久久久久久婷婷小说| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 一级,二级,三级黄色视频| 91久久精品国产一区二区三区| 狂野欧美激情性bbbbbb| 捣出白浆h1v1| 日日撸夜夜添| 黄色视频在线播放观看不卡| 久久av网站| 99视频精品全部免费 在线| 午夜日本视频在线| 亚洲精品色激情综合| 亚洲成国产人片在线观看| 欧美日韩视频高清一区二区三区二| 日韩大片免费观看网站| 久久久久久人妻| 毛片一级片免费看久久久久| 精品亚洲乱码少妇综合久久| 日本黄大片高清| 国产成人精品婷婷| 国产亚洲一区二区精品| 亚洲成人av在线免费| 各种免费的搞黄视频| 国产免费一区二区三区四区乱码| 9热在线视频观看99| 性色av一级| 国产69精品久久久久777片| 丝瓜视频免费看黄片| 男人添女人高潮全过程视频| 国产精品久久久av美女十八| 国产乱人偷精品视频| 亚洲精品一区蜜桃| 精品亚洲成a人片在线观看| 亚洲国产看品久久| 一级爰片在线观看| 国产成人一区二区在线| 国产一区二区三区综合在线观看 | 日韩中文字幕视频在线看片| 中文字幕最新亚洲高清| 亚洲精品av麻豆狂野| 国产 精品1| 欧美精品高潮呻吟av久久| 少妇高潮的动态图| 观看av在线不卡| 综合色丁香网| 晚上一个人看的免费电影| 亚洲第一av免费看| 国产精品.久久久| 国产av精品麻豆| 国产一级毛片在线| 啦啦啦啦在线视频资源| 国产一区有黄有色的免费视频| 你懂的网址亚洲精品在线观看| 欧美亚洲 丝袜 人妻 在线| 欧美bdsm另类| 久久久久久久国产电影| 国产亚洲一区二区精品| 在线 av 中文字幕| 国产亚洲一区二区精品| 婷婷成人精品国产| 9色porny在线观看| 男男h啪啪无遮挡| 一本久久精品| 中文字幕亚洲精品专区| 国产爽快片一区二区三区| 午夜免费观看性视频| 亚洲精品美女久久av网站| av网站免费在线观看视频| 亚洲熟女精品中文字幕| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩一区二区视频在线观看视频在线| 人妻少妇偷人精品九色| 国产免费又黄又爽又色| 女的被弄到高潮叫床怎么办| 日本欧美视频一区| 欧美xxⅹ黑人| 啦啦啦在线观看免费高清www| 一区二区三区乱码不卡18| 久久精品国产亚洲av天美| 欧美日韩综合久久久久久| 免费观看a级毛片全部| 免费大片18禁| 欧美成人精品欧美一级黄| 熟女av电影| a级毛片在线看网站| 大片电影免费在线观看免费| 丝袜美足系列| 日本91视频免费播放| 国产成人精品福利久久| √禁漫天堂资源中文www| 在线 av 中文字幕| 人妻 亚洲 视频| 久久99一区二区三区| 欧美bdsm另类| 亚洲国产欧美在线一区| 亚洲av福利一区| 啦啦啦中文免费视频观看日本| 好男人视频免费观看在线| 国产日韩一区二区三区精品不卡| 久久99精品国语久久久| 国产福利在线免费观看视频| 欧美日韩视频精品一区| 亚洲国产精品专区欧美| 国国产精品蜜臀av免费| 国产在视频线精品| 99久国产av精品国产电影| 另类精品久久| 成年动漫av网址| 色94色欧美一区二区| 亚洲精品,欧美精品| 国产日韩一区二区三区精品不卡| 午夜老司机福利剧场| 国产亚洲一区二区精品| 国产精品蜜桃在线观看| 欧美老熟妇乱子伦牲交| 国产免费一级a男人的天堂| 国产一区亚洲一区在线观看| 亚洲国产精品专区欧美| 国产亚洲欧美精品永久| 七月丁香在线播放| 精品久久蜜臀av无| 超色免费av| 交换朋友夫妻互换小说| 精品福利永久在线观看| 日本av免费视频播放| 欧美精品一区二区免费开放| 亚洲精品自拍成人| 少妇人妻 视频| 大码成人一级视频| 大香蕉97超碰在线| 大香蕉久久网| 最黄视频免费看| 欧美精品av麻豆av| 韩国高清视频一区二区三区| 九草在线视频观看| 亚洲伊人久久精品综合| 性色avwww在线观看| 男女午夜视频在线观看 | 高清黄色对白视频在线免费看| h视频一区二区三区| 色哟哟·www| 午夜激情久久久久久久| 亚洲内射少妇av| 久久99热这里只频精品6学生| 成年av动漫网址| 国产在线免费精品| 大片电影免费在线观看免费| 久久久精品区二区三区| 亚洲精品aⅴ在线观看| 男女国产视频网站| 日韩制服丝袜自拍偷拍| 日本-黄色视频高清免费观看| 啦啦啦在线观看免费高清www| 亚洲精品第二区| 亚洲欧美成人精品一区二区| 亚洲美女视频黄频| 九色亚洲精品在线播放| 日韩熟女老妇一区二区性免费视频| 精品一区在线观看国产| 欧美成人午夜免费资源| 国产免费福利视频在线观看| 少妇人妻久久综合中文| 国产日韩欧美视频二区| 久久久久久久久久成人| 美女主播在线视频| 国产精品熟女久久久久浪| av在线app专区| 国产精品一区二区在线观看99| 99热网站在线观看| a 毛片基地| 中国美白少妇内射xxxbb| 亚洲色图 男人天堂 中文字幕 | 嫩草影院入口| 精品午夜福利在线看| 美女视频免费永久观看网站| 99精国产麻豆久久婷婷| 一二三四中文在线观看免费高清| 欧美 日韩 精品 国产| 久久国产精品大桥未久av| 最近最新中文字幕免费大全7| 欧美日韩亚洲高清精品| 成人毛片a级毛片在线播放| 国产视频首页在线观看| 日韩一区二区视频免费看| 亚洲熟女精品中文字幕| 亚洲欧美日韩另类电影网站| 人人妻人人澡人人爽人人夜夜| 精品久久国产蜜桃| 日本91视频免费播放| 男女边摸边吃奶| 日本欧美国产在线视频| 免费大片黄手机在线观看| 国产男女内射视频| 国产成人aa在线观看| 国产一区二区激情短视频 | 制服人妻中文乱码| 在线亚洲精品国产二区图片欧美| 狠狠婷婷综合久久久久久88av| 91在线精品国自产拍蜜月| 97在线人人人人妻| 黄片无遮挡物在线观看| 黄色视频在线播放观看不卡| 波野结衣二区三区在线| 多毛熟女@视频| 国产老妇伦熟女老妇高清| 美女中出高潮动态图| freevideosex欧美| 好男人视频免费观看在线| 秋霞伦理黄片| 欧美人与性动交α欧美精品济南到 | 看免费成人av毛片| 99久久中文字幕三级久久日本| 精品99又大又爽又粗少妇毛片| 亚洲欧美精品自产自拍| 日本爱情动作片www.在线观看| 国产一区有黄有色的免费视频| 母亲3免费完整高清在线观看 | freevideosex欧美| 99九九在线精品视频| 国产精品.久久久| 久久久国产一区二区| 超碰97精品在线观看| 黄色配什么色好看| 97在线视频观看| 中文字幕免费在线视频6| 人体艺术视频欧美日本| 午夜激情久久久久久久| 亚洲国产毛片av蜜桃av| 国产 一区精品| 美女内射精品一级片tv| 80岁老熟妇乱子伦牲交| 成人手机av| 欧美xxxx性猛交bbbb| 大香蕉久久网| 一级毛片电影观看| 日韩制服丝袜自拍偷拍| 街头女战士在线观看网站| 人妻少妇偷人精品九色| 精品酒店卫生间| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 日韩av在线免费看完整版不卡| 美女国产高潮福利片在线看| www.av在线官网国产| 亚洲av电影在线观看一区二区三区| 自线自在国产av| 高清黄色对白视频在线免费看| 丰满乱子伦码专区| 97在线人人人人妻| 国产精品国产三级专区第一集| 国产成人av激情在线播放| 99久久人妻综合| 欧美日韩精品成人综合77777| 国产亚洲精品久久久com| 亚洲国产精品999| 久久久久人妻精品一区果冻| 男女午夜视频在线观看 | 妹子高潮喷水视频| 亚洲一码二码三码区别大吗| 人妻人人澡人人爽人人| 欧美激情极品国产一区二区三区 | 18禁观看日本| 这个男人来自地球电影免费观看 | 国产精品 国内视频| 日日爽夜夜爽网站| av视频免费观看在线观看| 26uuu在线亚洲综合色| 五月天丁香电影| 伦理电影大哥的女人| 久久97久久精品| 成人综合一区亚洲| 国产麻豆69| 一边摸一边做爽爽视频免费| 亚洲中文av在线| 一边亲一边摸免费视频| 国产乱人偷精品视频| 黄片无遮挡物在线观看| av不卡在线播放| 全区人妻精品视频| 91精品伊人久久大香线蕉| 欧美变态另类bdsm刘玥| 十分钟在线观看高清视频www| 丁香六月天网| 久久99精品国语久久久| 97人妻天天添夜夜摸| 国产成人精品在线电影| 日本vs欧美在线观看视频| 中文字幕免费在线视频6| 一区二区三区精品91| 日韩免费高清中文字幕av| 亚洲久久久国产精品| 亚洲精品乱久久久久久| 97在线人人人人妻| 国产成人精品久久久久久| 中文精品一卡2卡3卡4更新| 搡老乐熟女国产| 免费黄色在线免费观看| 男女午夜视频在线观看 | 99re6热这里在线精品视频| 亚洲国产成人一精品久久久| 丁香六月天网| 国产成人一区二区在线| 成人手机av| 在线观看国产h片| 欧美3d第一页| 日本91视频免费播放| 久久久久国产网址| 亚洲精品视频女| 久久ye,这里只有精品| 老司机影院成人| 精品亚洲成国产av| 捣出白浆h1v1| 波多野结衣一区麻豆| 亚洲综合精品二区| 亚洲性久久影院| 亚洲精品,欧美精品| 中文字幕最新亚洲高清| 国产精品久久久久久精品电影小说| 午夜精品国产一区二区电影| 亚洲精品美女久久久久99蜜臀 | 久久久精品94久久精品| 青春草国产在线视频| 精品人妻熟女毛片av久久网站| 26uuu在线亚洲综合色| 在线观看国产h片| 久久久久国产精品人妻一区二区| 日韩电影二区| 香蕉丝袜av| 多毛熟女@视频| 欧美成人午夜精品| 天天影视国产精品| 永久免费av网站大全| 国产精品一区二区在线不卡| 中文字幕亚洲精品专区| 视频在线观看一区二区三区| av一本久久久久| 亚洲第一av免费看| 中文天堂在线官网| 免费久久久久久久精品成人欧美视频 | 国产亚洲精品第一综合不卡 | 一边亲一边摸免费视频| 精品国产国语对白av| 精品人妻熟女毛片av久久网站| 久久精品久久久久久噜噜老黄| 亚洲精品久久午夜乱码| 日韩免费高清中文字幕av| 少妇被粗大的猛进出69影院 | 久久久久国产网址| 中文字幕另类日韩欧美亚洲嫩草| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩另类电影网站| 午夜免费男女啪啪视频观看| 男人操女人黄网站| 亚洲国产最新在线播放| 黄网站色视频无遮挡免费观看| 日韩精品免费视频一区二区三区 | 满18在线观看网站| 2018国产大陆天天弄谢| 午夜免费鲁丝| 伊人久久国产一区二区| 国产精品无大码| 亚洲国产最新在线播放| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 久久久久久久国产电影| 中文字幕人妻丝袜制服| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区| 最近中文字幕2019免费版| 狂野欧美激情性xxxx在线观看| 韩国高清视频一区二区三区| 亚洲精品一二三| 午夜av观看不卡| 一区二区三区四区激情视频| 亚洲欧美日韩另类电影网站| 夫妻性生交免费视频一级片| 一级毛片电影观看| 99久久人妻综合| 国产免费一级a男人的天堂| 午夜免费鲁丝| 欧美人与善性xxx| 七月丁香在线播放| 激情视频va一区二区三区| 免费观看av网站的网址| 性色avwww在线观看| 精品一品国产午夜福利视频| 欧美成人午夜免费资源| 亚洲美女搞黄在线观看| 国产女主播在线喷水免费视频网站| 国产精品秋霞免费鲁丝片| 999精品在线视频| 国产有黄有色有爽视频| 色视频在线一区二区三区| 午夜影院在线不卡| 成人黄色视频免费在线看| 在线免费观看不下载黄p国产| 深夜精品福利| 国产精品人妻久久久影院| 高清毛片免费看| 天天操日日干夜夜撸| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 欧美国产精品一级二级三级| 九九爱精品视频在线观看| 边亲边吃奶的免费视频| www.色视频.com| 青春草国产在线视频| 赤兔流量卡办理| 色婷婷av一区二区三区视频| 国产一级毛片在线| 永久免费av网站大全| av天堂久久9| a级片在线免费高清观看视频| 一本—道久久a久久精品蜜桃钙片| 久久久久久久国产电影| 亚洲一区二区三区欧美精品| 国产精品99久久99久久久不卡 | 国产69精品久久久久777片| 日韩在线高清观看一区二区三区| 免费久久久久久久精品成人欧美视频 | 国精品久久久久久国模美| 国产色婷婷99| 少妇高潮的动态图| 黑人高潮一二区| 晚上一个人看的免费电影| 欧美xxⅹ黑人| 国产高清国产精品国产三级| 国产淫语在线视频| 午夜av观看不卡| 精品久久久精品久久久| 90打野战视频偷拍视频| 观看美女的网站| 国产淫语在线视频| 欧美精品一区二区免费开放| 狂野欧美激情性bbbbbb| 女人久久www免费人成看片| 国产成人精品婷婷| 狠狠精品人妻久久久久久综合| 国产成人精品在线电影| 国产精品欧美亚洲77777| 亚洲国产毛片av蜜桃av| 国产熟女午夜一区二区三区| 女性被躁到高潮视频| 美女脱内裤让男人舔精品视频| videossex国产| 欧美+日韩+精品| 免费久久久久久久精品成人欧美视频 | 亚洲欧美日韩卡通动漫| 欧美激情极品国产一区二区三区 | 国产免费一区二区三区四区乱码| 国产精品.久久久| 色视频在线一区二区三区| 男女无遮挡免费网站观看| 亚洲欧洲日产国产| 乱码一卡2卡4卡精品| 中文字幕人妻熟女乱码| 丝瓜视频免费看黄片| a级毛片在线看网站| 欧美成人午夜精品| 色吧在线观看| 又黄又爽又刺激的免费视频.| 观看av在线不卡| 欧美成人午夜免费资源| 一区二区三区精品91| www日本在线高清视频| 免费大片黄手机在线观看| 亚洲,一卡二卡三卡| 80岁老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 新久久久久国产一级毛片| 亚洲国产精品专区欧美| 晚上一个人看的免费电影| 国产精品 国内视频| 欧美人与性动交α欧美软件 | 激情视频va一区二区三区| 国产又爽黄色视频| 亚洲综合色惰| 18禁观看日本| 久久99一区二区三区| 九色亚洲精品在线播放| 欧美精品一区二区免费开放| 久久久久久人人人人人| 免费观看a级毛片全部| 18禁在线无遮挡免费观看视频| 国产男女超爽视频在线观看| 精品国产国语对白av| 久久久久精品人妻al黑| 精品午夜福利在线看| 久热久热在线精品观看| 最新的欧美精品一区二区| 婷婷色av中文字幕| 大香蕉久久成人网| 人人妻人人添人人爽欧美一区卜| 国产成人一区二区在线| 日本猛色少妇xxxxx猛交久久| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区日韩欧美中文字幕 | 两性夫妻黄色片 | 国产综合精华液| 亚洲人与动物交配视频| 亚洲精品,欧美精品| 侵犯人妻中文字幕一二三四区| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区国产| 成人手机av| 国产成人aa在线观看| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 欧美精品一区二区免费开放| 欧美成人午夜免费资源| 久久99精品国语久久久| 婷婷色av中文字幕| 2022亚洲国产成人精品| 精品一区二区三区视频在线| 精品人妻偷拍中文字幕| 好男人视频免费观看在线| 欧美精品一区二区大全| 少妇的逼好多水| 国产亚洲一区二区精品| 国产av码专区亚洲av| 嫩草影院入口| 亚洲精品美女久久久久99蜜臀 | 久久久国产一区二区| 亚洲精华国产精华液的使用体验| 久久ye,这里只有精品| 免费大片黄手机在线观看|