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

    氣溶膠對利用近紅外波段反演CO2濃度精度的影響研究

    2016-04-14 10:50:02王舒鵬張興贏王維和陳林方莉國家衛(wèi)星氣象中心北京0008中國科學(xué)院遙感與數(shù)字地球研究所北京000
    關(guān)鍵詞:氣溶膠

    王舒鵬 張興贏 王維和 陳林 方莉( 國家衛(wèi)星氣象中心,北京 0008; 中國科學(xué)院遙感與數(shù)字地球研究所,北京 000)

    ?

    氣溶膠對利用近紅外波段反演CO2濃度精度的影響研究

    王舒鵬1張興贏1王維和1陳林1方莉2
    (1 國家衛(wèi)星氣象中心,北京 100081;2 中國科學(xué)院遙感與數(shù)字地球研究所,北京 100101)

    摘要:近紅外波段CO2濃度反演誤差主要與氣溶膠散射作用引起的吸收光程改變量難以精確估計有關(guān)。此外,不同的氣溶膠模式及下墊面反射率條件下吸收光程改變量也不盡相同。針對不同氣溶膠模式條件下,忽略氣溶膠散射對CO2

    濃度反演精度的影響進(jìn)行分析。研究發(fā)現(xiàn),對鄉(xiāng)村型、對流層型和海洋型氣溶膠來說,忽略氣溶膠散射效應(yīng),在暗地表時(反射率<0.1)會導(dǎo)致CO2濃度反演結(jié)果低估;當(dāng)?shù)乇矸瓷渎食^0.1后,CO2濃度反演結(jié)果均為高估,且隨地表反射率的增加,反演誤差增加。對城市型氣溶膠模式來說,忽略氣溶膠影響均會導(dǎo)致CO2濃度反演結(jié)果低估,隨地表反射率增加,反演誤差遞減。在典型觀測幾何下,地表反射率(1.6μm)和氣溶膠光學(xué)厚度(0.55μm)在0.1~0.3范圍內(nèi)時,忽略氣溶膠影響,城市型、海洋型、鄉(xiāng)村型以及對流層氣溶膠分別會引入-0.1%~-0.5%、0.22%~1.92%、0.09%~1.46%及0.02%~0.45%的反演誤差。在利用GOSAT觀測光譜對XCO2(干燥空氣下CO2平均混合比)進(jìn)行反演時發(fā)現(xiàn),相比GOSAT發(fā)布結(jié)果(389.814ppm,1ppm=10-6),基于地基實測氣溶膠特性數(shù)據(jù)的XCO2反演結(jié)果(390.95ppm)與地基觀測結(jié)果(390.737ppm)有更好的一致性。

    關(guān)鍵詞:XCO2反演,近紅外波段,GOSAT,氣溶膠

    第一作者:王舒鵬(1984—),Email:wangsp@cma.gov.cn

    資助信息:國家高技術(shù)研究發(fā)展計劃(2011AA12A104);國家自然科學(xué)基金(41306185)

    0 引言

    基于近紅外高光譜衛(wèi)星遙感的CO2監(jiān)測,由于其高時空分辨率、全球覆蓋、低成本以及對近地面CO2濃度變化的高敏感性,逐漸成為監(jiān)測與人類活動相關(guān)的CO2源匯信息的重要手段[1-2]。2002年搭載于ENVISAT上的SCIAMACHY[3],2009年發(fā)射的近紅外波段的高光譜探測儀GOSAT[4]和2014年美國發(fā)射的OCO-2[2]都可以通過觀測近紅外波段的太陽光反射光譜來反演近地面CO2含量。即將發(fā)射的下一代近紅外高光譜溫室氣體載荷,包括歐洲的CarbonSat[5]以及中國的TanSat等多個載荷[6-7],都將為更好地獲取近地面溫室氣體的源匯信息提供支持。

    影響近紅外太陽反射光譜反演CO2濃度的要素眾多[8]。但最大的不確定性來自大氣中氣溶膠的散射作用所帶來的光程改變的不確定性[9]。氣溶膠的散射作用可能增加光程,也可能減小光程。在下墊面反射率較低時,傳感器接收到的光子通過大氣散射直接進(jìn)入傳感器而不經(jīng)過地表反射,此時的光程是減小的。相反,在下墊面反射率較高時,到達(dá)傳感器的光子由于氣溶膠或云與地表的多次散射作用,光程增加[10]。GOSAT 官方發(fā)布的CO2反演產(chǎn)品的地基驗證結(jié)果表明,忽略氣溶膠和云的散射作用對光程的改變,即使在晴空條件下,都可能對CO2反演精度帶來較大的誤差[11-12]。研究發(fā)現(xiàn),撒哈拉沙漠的沙塵型氣溶膠可以對SCIAMACHY的CO2總量反演誤差達(dá)到10%[13],遠(yuǎn)超過在區(qū)域尺度上(約1000km×1000km)監(jiān)測碳源匯的需求精度(<1%)[14]。以往很多研究基于單次散射假設(shè)對氣溶膠散射作用引起的CO2濃度反演誤差進(jìn)行定量化分析。也有研究在氣溶膠多次散射假設(shè)下,對氣溶膠散射引起的CO2濃度反演誤差進(jìn)行分析[10,15],但并沒有考慮不同下墊面反射對反演精度的影響??紤]到地表反射率是決定有效光程的重要參數(shù),Aben等[16]在更大范圍的地表反射率條件下,分析了不同氣溶膠光學(xué)厚度和薄卷云存在條件下的散射誤差。楊東旭等[17]通過不同氣溶膠模式對1.6μm波段強(qiáng)弱吸收帶的輻射光譜比的影響來評價氣溶膠散射對反演的影響。

    考慮到不同氣溶膠模式和下墊面反射率條件下,氣溶膠散射效應(yīng)并不相同,本研究重點在不同地表反射率和不同氣溶膠模式條件下,對忽略氣溶膠散射對XCO2反演精度的影響進(jìn)行分析。此外,本文在實測氣溶膠模式和幾種典型氣溶膠模式假設(shè)下,分別利用GOSAT觀測光譜分別對XCO2進(jìn)行反演,并對不同氣溶膠假設(shè)引起的反演誤差進(jìn)行分析。

    1 正向輻射傳輸模擬及反演

    本文所采用的正向輻射傳輸模型是SCIATRAN 3.2[18],該輻射傳輸模式基于離散坐標(biāo)法對輻射傳輸方程進(jìn)行求解,為簡便起見,在求解中忽略偏振貢獻(xiàn)。在模擬傳感器的卷積效應(yīng)時,采用Boxcar型響應(yīng)函數(shù),傳感器各波段半波寬(HWMF)與GOSAT對應(yīng)波段參數(shù)相同,1.6μm波段和2.05μm波段分別約為0.03512nm和0.054nm。

    在計算吸收氣體的消光系數(shù)時,為了提高計算精度,采用逐線積分(LBL)的方式計算,采樣點間隔為0.001nm。在輻射傳輸計算時,將整個大氣分為20層。在所選波段內(nèi)存在的吸收氣體為CO2和H2O。這兩種吸收氣體每一層的消光系數(shù)都通過該層所包含的各子層的消光系數(shù)相加得到。每一子層內(nèi)近似認(rèn)為壓強(qiáng)恒定。太陽吸收光譜來自于Kurucz[19]的太陽吸收線表。吸收氣體的吸收截面通過HITRAN2008數(shù)據(jù)庫[20]中的吸收線參數(shù)計算獲得。圖1給出在GOSAT光譜分辨率下,1.6μm波段使用HITRAN2012[21]和HITRAN2008的光譜差異。雖然相比2008版,2012版中CO2的吸收線條數(shù)增加了許多,但是兩者的總線強(qiáng)幾乎相同,不同HITRAN數(shù)據(jù)庫版本得到的透過率相差很小,只有在少數(shù)CO2強(qiáng)吸收的地方存在一些波動,卷積后兩者差在整個波段上幾乎為0,可以認(rèn)為2008版和2012版在該波段中對CO2的敏感性相同。

    圖1 使用HITRAN2008和HITRAN2012光譜庫得到的模擬光譜差異(1.6μm波段)Fig.1 Difference between the spectrum at 1.6μm simulated using HITRAN2008 and HITRAN2012

    在模擬過程中同時考慮大氣分子和氣溶膠粒子的散射影響??紤]到在實際CO2反演中,有云存在的條件下反演結(jié)果可信度往往較低,假定模擬場景為無云情況。通過瑞利散射假設(shè)對大氣分子散射進(jìn)行模擬[22],其對應(yīng)的光學(xué)厚度是波長和大氣分子數(shù)密度的函數(shù)。在模擬過程中選擇的不同氣溶膠模式來自于LOWTRAN氣溶膠數(shù)據(jù)庫[23]。不同的氣溶膠模式為:海洋型氣溶膠(能見度:23km,相對濕度:70%),城市型氣溶膠(能見度:10km,相對濕度:70%),鄉(xiāng)村型氣溶膠(能見度:23km,相對濕度:70%)和一般對流層氣溶膠(能見度:23km,相對濕度:70%)。氣溶膠光學(xué)厚度(AOD)的變化范圍為0.05~1.0(0.55μm)。氣溶膠散射特性模擬基于球形粒子假設(shè),通過Mie散射對其散射相函數(shù)進(jìn)行模擬。假定地表為朗伯體,兩個CO2吸收波段地表反射率相同。利用迭代的非線性最小二乘法視線光譜匹配[24],對待反演參數(shù)進(jìn)行反演。表1給出正向模擬時不同場景的參數(shù)變化范圍。

    表1 不同氣溶膠模式的單次散射反照率(SSA)和不對稱因子(g)(1.536μm)Table 1 The SSA and asymmetry factor of different aerosol models(at 1.536μm)

    2 忽略氣溶膠散射對反演精度影響分析

    大氣氣溶膠對CO2反演的主要影響體現(xiàn)在其散射作用會改變到達(dá)傳感器的入瞳光的光程。雖然在反演波段內(nèi)氣溶膠的散射特性相對一致,該光程改變量的估計誤差仍會直接影響CO2反演精度。氣溶膠的光學(xué)特性(AOD)、微物理特性(SSA和不對稱因子)、氣溶膠粒子的垂直分布、下墊面反射率以及觀測幾何都會影響由散射作用引起的光程改變量,甚至光程改變的符號。

    傳感器接收到的光強(qiáng)Itot可以表達(dá)為[10]:

    其中,Idir代表直射太陽光,即沒有經(jīng)過大氣散射、僅經(jīng)過一次地表反射后直接進(jìn)入傳感器的光強(qiáng);Ibsct代表大氣后向散射,即僅經(jīng)過大氣中的分子和氣溶膠后向散射進(jìn)入傳感器的光強(qiáng);Isrfl代表大氣散射光經(jīng)地表反射直接進(jìn)入傳感器的光強(qiáng);Irsct代表直射太陽光經(jīng)地表反射后又經(jīng)過一次大氣散射進(jìn)入傳感器的光強(qiáng);Imsct代表經(jīng)過大氣和地表多次散射進(jìn)入傳感器的光強(qiáng)。若沒有氣溶膠散射存在,傳感器件接收到的光強(qiáng)以直射光項Idir為主,大氣散射項可以忽略。在氣溶膠散射存在的情況下,大氣散射項不能忽略,與直射光項Idir相比,大氣后向散射項Ibsct引起吸收光程縮短,而前向散射項(Isrfl,Irsct,Imsct)引起吸收光程增加。

    考慮到大氣散射對短波紅外波段光程的改變是氣溶膠散射與下墊面反射聯(lián)合作用的結(jié)果[25]。圖2給出對應(yīng)于LOWTRAN的不同氣溶膠模式(鄉(xiāng)村型、對流層型、海洋型和城市型氣溶膠)在不同AOD(0.55μm)下,太陽天頂角為30°、觀測天頂角為0°時,如忽略氣溶膠散射,CO2柱濃度反演精度隨地表反射率的變化。各相關(guān)要素的變化范圍見表2。

    圖2 鄉(xiāng)村型(a)、對流層型(b)、海洋型(c)和城市型(d)氣溶膠模式在不同AOD(0.55μm)情況下,忽略氣溶膠存在時,XCO2反演相對誤差(真值為385ppm)隨地表反射率的變化Fig.2 Variation of the XCO2retrieval error introduced by ignoring aerosol scattering with the surface reflectance for different AOD(0.55μm)under(a)rural,(b)tropospheric,(c)marine and(d)urban aerosol model

    表2 氣溶膠對CO2反演精度影響分析中各要素的變化范圍Table 2 The variation range of parameters used in the study on CO2retrieval accuracy affected by aerosols

    在相同AOD和地表反射率情況下,忽略海洋型和城市型氣溶膠散射的反演誤差明顯大于忽略鄉(xiāng)村型和對流層型氣溶膠的誤差(AOD>0.5時,反演誤差可能超過15%)。在海洋型和城市型氣溶膠模式下忽略氣溶膠散射影響進(jìn)行反演時,為了保證反演結(jié)果的合理性,需要對模擬光譜進(jìn)行多項式校正,減小由于氣溶膠寬通道信號引入的反演誤差。

    對鄉(xiāng)村型、對流層型和海洋型氣溶膠模式,在相同AOD下,當(dāng)?shù)乇矸瓷渎瘦^低(<0.1)時,忽略氣溶膠散射影響會導(dǎo)致XCO2反演結(jié)果低估;當(dāng)?shù)乇矸瓷渎食^0.1時,XCO2反演結(jié)果高估且隨地表反射率的增加反演誤差逐漸增加;在相同地表反射率條件下,隨AOD增加,反演誤差增加。這是由于在暗地表時,傳感器接收到的光強(qiáng)主要來自于大氣后向散射項Ibsct,相比忽略氣溶膠散射時的主要貢獻(xiàn)項Idir,吸收光程縮短,因此反演的CO2含量偏低。隨著地表反射率的增加,大氣前向散射逐漸成為傳感器入瞳光的主要貢獻(xiàn)。此時前向散射項(Isrfl,Irsct,Imsct)的貢獻(xiàn)增加,后向散射項(Ibsct)貢獻(xiàn)減弱,前向散射項對應(yīng)的吸收光程要大于直射光項,因此反演的CO2含量增加。城市型氣溶膠模式與上面三種氣溶膠模式不同:在相同AOD情況下,忽略氣溶膠散射影響,均導(dǎo)致CO2反演結(jié)果低估;隨地表反射率增加,誤差遞減,即暗地表反演誤差最大。在相同地表反射率條件下,隨AOD增加,反演誤差增加。

    當(dāng)?shù)乇矸瓷渎剩?.6μm)和AOD(0.55μm)在0.1~0.3時,CO2柱總量反演誤差對城市型氣溶膠為-0.1%~-0.5%,海洋型氣溶膠為0.22%~1.92%,鄉(xiāng)村型氣溶膠為0.09%~1.46%,對流層氣溶膠為0.02%~0.45%。本文的結(jié)果與Aben等[16]得到的結(jié)論基本一致。Aben等[16]發(fā)現(xiàn)強(qiáng)吸收的城市氣溶膠會引起反演結(jié)果偏低,而其他氣溶膠模式除在暗地表會出現(xiàn)反演結(jié)果偏低外,當(dāng)?shù)乇矸瓷渎食^0.1時都會出現(xiàn)反演結(jié)果偏高。本文得到的海洋型氣溶膠和城市型氣溶膠的反演誤差范圍小于Aben等[16]的結(jié)果(分別引入0.1%~2.8%和-1.1%~0的反演誤差)。這是由于本文在這兩種氣溶膠模式下利用二階多項式對寬通道信號進(jìn)行校正,壓制了反演誤差所致。Mao等[10]通過強(qiáng)弱吸收帶的輻射光譜比的變化來研究不同氣溶膠模式對XCO2反演精度的影響,發(fā)現(xiàn)海洋型氣溶膠對XCO2反演精度影響最大(~15ppm),吸收性城市型氣溶膠影響最?。ā?ppm)。與本文結(jié)論差異的原因在于,基于強(qiáng)弱吸收帶的輻射光譜比的校正可以大幅度削弱地表反射率或氣溶膠引入的寬通道信號所導(dǎo)致的誤差,但對由氣溶膠散射引起的吸收光程變化造成的誤差削弱效果較差,即對城市吸收型氣溶膠(SSA為0.52)引入的誤差校正最好,而對在該波段散射性較強(qiáng)的海洋型氣溶膠(SSA為0.98)引入的誤差校正效果較差。該差異也說明,無論反演誤差是由于氣溶膠或地表吸收還是由于散射引起的光程改變所引入,基于二次多項式擬合寬通道信號的校正結(jié)果都較好。如果反演中無法獲取氣溶膠特性信息,可以考慮用該校正方法提高反演精度。

    3 反演實例分析

    2013年8月在敦煌定標(biāo)場(40.19°N,94.31°E)開展了針對GOSAT L2級CO2產(chǎn)品的精度評價試驗。為了進(jìn)一步說明氣溶膠模式的不確定性對CO2濃度反演的影響,在地基實際氣溶膠參數(shù)和以上不同氣溶膠模式假設(shè)條件下分別利用GOSAT觀測光譜對XCO2進(jìn)行反演。

    考慮到CO2反演需要無云場景,采用GOSAT L2 級 CAI 云檢測產(chǎn)品對無云場景進(jìn)行辨識。經(jīng)過比較,在地基同步觀測時間內(nèi)(2013年8月19—25日),僅有8月24日的GOSAT對應(yīng)點觀測光譜數(shù)據(jù)滿足良好的無云場景條件。表3給出觀測站位置信息及同步觀測的環(huán)境參數(shù)信息。選擇與地基觀測站點位置最為接近的曝光點數(shù)據(jù)進(jìn)行反演,利用地基觀測數(shù)據(jù)得到反演時刻附近的環(huán)境及大氣參數(shù)。圖3給出2013年8 月24日拉曼偏振大氣探測激光雷達(dá)在測得的GOSAT過境前后0.532μm大氣消光系數(shù)廓線。圖4和圖5給出0.532μm消光系數(shù)時空演化圖和0.532μm偏振系數(shù)時空演化圖??梢钥闯鲈囼灝?dāng)日,敦煌地區(qū)大氣潔凈,上午9:00—夜間20:00,0.532μm的大氣光學(xué)厚度在0.1以下,且整日比較穩(wěn)定。GOSAT過境時(15:30前后),大氣中氣溶膠粒子濃度較上午有所增長。整日大氣偏振系數(shù)大多在0.15以下,大氣中氣溶膠主要為規(guī)則的細(xì)粒子。14:00以后3km以下空氣中有較薄的沙塵傳輸。由消光系數(shù)廓線可以看出,當(dāng)日敦煌地區(qū)大氣邊界層約為3km,峰值在1.5km處。圖6為2013年8月24日10:00—19:00 10km高度內(nèi)大氣粒子有效半徑時空演化圖,結(jié)果顯示GOSAT過境前后,地面到3km之間粒子有效半徑主要集中在1~1.5μm。圖7給出基于地基CE318觀測數(shù)據(jù)反演得到的GOSAT過境前后氣溶膠粒子體積譜分布。

    圖3 2013年8月24日GOSAT過境前后敦煌站0.532μm大氣消光系數(shù)廓線Fig.3 The atmospheric extinction profiles at 0.532μm measured at Dunhuang near GOSAT overpass on August 24th,2013

    圖4 2013年8月24日敦煌站0.532μm消光系數(shù)時空演化圖Fig.4 The space-time diagram of atmospheric extinction at 0.532μm measured at Dunhuang on August 24th,2013

    圖5 2013年8月24日敦煌站0.532μm偏振系數(shù)時空演化圖Fig.5 The space-time diagram of aerosol polarization coefficient at 0.532μm measured at Dunhuang on August 24th,2013

    圖6 2013年8月24日敦煌站粒子有效半徑時空演化圖Fig.6 The space-time diagram of aerosol effective radius measured at Dunhuang on August 24th,2013

    圖7 基于CE318觀測數(shù)據(jù)反演得到的GOSAT過境前后氣溶膠粒子體積譜分布Fig.7 The aerosol volume size distribution during the GOSAT overpass retrieved from CE-318 measurements

    表3 觀測站位置信息及同步觀測的環(huán)境參數(shù)Table 3 The location of the ground observation station and simultaneous environmental parameters

    根據(jù)表3所列的環(huán)境參數(shù)先驗值,反演XCO2。表4給出反演結(jié)果、GOSAT L2產(chǎn)品和地基觀測結(jié)果的比較。在反演中除了XCO2,還對地表反射率,大氣垂直柱水汽含量、溫度和壓強(qiáng)廓線修正量進(jìn)行反演。根據(jù)反演結(jié)果,地表反射率和先驗觀測值差異明顯,在1.6μm和2.0μm波段分別為0.289和0.291。這是由于地表反射率實際觀測值是在垂直方向上觀測獲得的,沒有考慮地表的方向反射特性。

    表4 XCO2反演結(jié)果比較Table 4 Comparison of XCO2results

    為了進(jìn)一步說明不同氣溶膠模式對XCO2反演結(jié)果的影響,圖8給出在實測AOD情況下,不同LOWTRAN氣溶膠模式和實測氣溶膠模式情況下XCO2反演結(jié)果。此時觀測幾何及其他環(huán)境參數(shù)均采用地基實測數(shù)據(jù)。相同AOD,不同氣溶膠模式下,城市型和海洋型氣溶膠引入的XCO2反演誤差最大,分別為-1.97ppm和1.58ppm;鄉(xiāng)村型氣溶膠引入的反演誤差最?。?.48ppm)。這可能是由于鄉(xiāng)村型氣溶膠與實際氣溶膠模式的差異相對較小,特別是SSA差異為0.01,而海洋型氣溶膠和城市型氣溶膠的SSA與實際偏差分別為0.12和-0.31。此時城市型氣溶膠模式假設(shè)下XCO2反演結(jié)果偏低與上文的分析結(jié)論是一致的。

    圖8 不同氣溶膠模式和實際觀測條件下XCO2反演結(jié)果Fig.8 The XCO2retrievals with different assumed aerosol types and the ground-based measured aerosol model.

    4 結(jié)論

    本文主要研究了不同氣溶膠模式條件下,忽略氣溶膠散射對近地面CO2濃度反演精度的影響。

    對鄉(xiāng)村型、對流層型和海洋型氣溶膠模式來說,在相同AOD條件下,在低地表反射率(<0.1)時,忽略氣溶膠散射影響會導(dǎo)致XCO2反演結(jié)果低估;當(dāng)?shù)乇矸瓷渎食^0.1時,XCO2反演結(jié)果高估,且隨地表反射率的增加,反演誤差逐漸增加。而對城市型氣溶膠模式來說,在相同AOD條件下,忽略氣溶膠影響,均會導(dǎo)致XCO2的反演結(jié)果低估;且隨地表反射率增加,誤差遞減。

    忽略海洋型和城市型氣溶膠散射影響,XCO2反演誤差明顯大于其他兩種氣溶膠(AOD>0.5條件下反演誤差可能超過15%)。采用二次多項式擬合寬通道信號校正觀測光譜后可以取得較好的反演結(jié)果,在反演中如果無法精確獲取氣溶膠特性信息可以采用該校正方法提高反演精度。

    在不同氣溶膠模式假設(shè)條件下利用GOSAT觀測光譜對XCO2進(jìn)行反演。通過反演結(jié)果對比發(fā)現(xiàn),基于地基實測氣溶膠特性數(shù)據(jù)的XCO2反演結(jié)果與地基觀測結(jié)果有更好的一致性。此外,在實測AOD、不同氣溶膠模式下,城市型和海洋型氣溶膠假設(shè)下引入的XCO2反演誤差最大,分別為-1.97ppm和1.58ppm;而鄉(xiāng)村型氣溶膠引入的反演誤差最?。?.48ppm),這主要是由于鄉(xiāng)村型氣溶膠與實際氣溶膠模式的SSA差異相對較小所致。

    參考文獻(xiàn)

    [1]IPCC.Climate Change 2007:The physical science basis.Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change.Cambridge,United Kingdom and New Tork,NY:Cambridge University Press,2007.

    [2]Boland S,B?sch H,Brown L,et al.The need for atmospheric carbon dioxide measurements from space:Contributions from a Rapid Reflight of the Orbiting Carbon Observatory.Technical Report,OCO Science Team,NASA Jet Propulsion Laboratory,2009.

    [3]Schneising O,Buchwitz M,Bovensmann H,et al.Three years of SCIAMACHY carbon dioxide and methane colμmn-averaged dry air mole fraction measurements.Proceedings ENVISAT Symposium.2007.

    [4]Yokota T,Yoshida Y,Eguchi N,et al.Global concentrations of CO2and CH4retrieved from GOSAT:First preliminary results.Sola,2009,5:160-163.

    [5]Bovensmann H,Buchwitz M,Burrows J P,et al.A remote sensing technique for global monitoring of power plant CO2emissions from space and related applications.Atmos Meas Tech,2010,3:781-811.

    [6]劉毅,呂達(dá)仁,陳洪濱,等.衛(wèi)星遙感大氣CO2的技術(shù)與方法進(jìn)展綜述.遙感技術(shù)與應(yīng)用,2011,26(2):247-254.

    [7]施海亮,熊偉,羅海燕,等.新型超光譜大氣CO2遙感探測技術(shù).光電工程,2013(8):36-41.

    [8]鄒銘敏,陳良富,陶金花,等.短波紅外通道CO2觀測的溫度敏感性分析.紅外與毫米波學(xué)報,2012,31(5):455-462.

    [9]Oshchepkov S,Bril A,Yokota T,et al.Effects of atmospheric light scattering on spectroscopic observations of greenhouse gases from space:Validation of PPDF-based CO2retrievals from GOSAT.Journal of Geophysical Research,2012,117:D12305-D12305.

    [10]Mao J,Kawa S R.Sensitivity studies for space-based measurement of atmospheric total column carbon dioxide by reflected sunlight.Applied Optics,2004,43(4):914-927.

    [11]Reuter M,Buchwitz M,Schneising O,et al.A method for improved SCIAMACHY CO2retrieval in the presence of optically thin clouds.Atmos Meas Tech,2010,3:209–232,.

    [12]Qu Y,Zhang C M,Wang D Y,et al.Comparison of atmospheric CO2observed by GOSAT and two ground stations in China,International Journal of Remote Sensing,2012,34:3938-3946.

    [13]Houweling S,Hartmann W,Aben I,et al.Evidence of systematic errors in SCIAMACHY-observed CO2due to aerosols.Atmos Chem Phys,2005,5:3003–3013.

    [14]Patra P K,Maksyutov S,Sasano Y,et al.An evaluation of CO2observations with solar occultationFTS for inclined-orbit satellite sensor for surface source inversions.J Geophys Res,2003,108.

    [15]Buchwitz M,Burrows J P.Retrieval of CH4,CO,and CO2total column amounts from SCIAMACHY near-infrared nadir spectra:retrieval algorithm and first results.Proceedings of SPIE 5235,remote sensing of clouds and the atmosphere Ⅷ,2004.

    [16]Aben I,Hasekamp O,Hartmann W.Uncertainties in the spacebased measurements of CO2columns due to scattering in the Earth’s atmosphere.J Quant Spectrosc Radiat.Transfer,2007,104:450–459.

    [17]Yang D X,Liu Y,Cai Z N.Simulations of aerosol optical properties to top of atmospheric reflected sunlight in the near infrared CO2weak absorption band.Atmos Oceanic Sci Lett,2013,6:60-64.

    [18]Rozanov V V,Rozanov A V,Kokhanovsky A A,et al.Radiative transfer through terrestrial atmosphere and ocean:software package SCIATRAN.Journal of Quantitative Spectroscopy and Radiative Transfer,2014,133:13-71.

    [19]Kurucz R L.New atlases for solar flux,irradiance,central intensity,and limb intensity.Memorie della Società Astronomica Italiana Supplement,2005,8:189.

    [20]Rothman L S,Gordon I E,Barbe A,et al.The HITRAN 2008 molecular spectroscopic database.Journal of Quantitative Spectroscopy and Radiative Transfer,2009,110:533-572.

    [21]Rothman L S,Gordon I E,Babikov Y,et al.The HITRAN2012 molecular spectroscopic database.Journal of Quantitative Spectroscopy and Radiative Transfer,2013,130:4-50.

    [22]Bodhaine B A,Wood N B,Dutton E G,et al.On Rayleigh optical depth calculations.Journal of Atmospheric and Oceanic Technology,1999,16:1854-1861

    [23]Shettle E P,and R W Fenn.Models of the atmospheric aerosols and their optical properties,in AGARD Conference Proceedings No.183,ADA028-615,1976.

    [24]Crisp D,Bosch H,Brown L,et al.OCO(Orbiting Carbon Observatory)-2 Level 2 full physics retrieval algorithm theoretical basis,Tech Rep OCO D-65488,NASA Jet Propulsion Laboratory,California Institute of Technology,Pasadena,CA,version 1.0 Rev 4.available at:http://disc.sci.gsfc.nasa.gov/acdisc/documentation/OCO-2_L2_FP_ATBD_v1_rev4_Nov10.pdf(last access:January 2012),2010.13042,13045.

    [25]Dubuisson P,Borde R,Schmechtig,C,et al.Surface pressure estimates from satellite data in the Oxygen A-band:Applications to the MOS sensor over land.J Geophys Res,2001,106:27277–27286.

    The lnfluence of Aerosol Particles on Retrieval Accuracy of CO2Concentration Using Near Infrared Bands

    Wang Shupeng1,Zhang Xingying1,Wang Weihe1,Chen Lin1,Fang Li2
    (1 National Satellite Meteorological Centre,Beijing 100081 2 Institute of Remote Sensing and Digital Earth,Chinese Academy of Sciences,Beijing 100101)

    Abstract:The main difficulty in retrieving CO2concentration from near-infrared observations comes from the uncertainty of light path change brought about by scattering of the atmospheric aerosols.The scattering effect described by different aerosol models and surface reflectances may lead to different light path changes.This paper aims to study the influence of neglecting aerosol scattering on CO2retrieval accuracy for different aerosol models.It is found that for rural,tropospheric and marine aerosol models,the neglecting of aerosol scattering results in underestimation of CO2retrievals at low surface reflectances(<0.1),but overestimation at surface reflectance over 0.1 and the increasing retrieval errors with increasing surface reflectances.For the urban aerosol model,the neglecting of aerosol scattering results in underestimation of CO2retrievals at any surface reflection,and the retrieval errors decrease with increasing surface reflectances.For typical observation geometry,a surface reflectance of 0.1-0.3 at 1.6μm and aerosol optical depth(AOD)of 0.1-0.3 at 0.55μm,the neglecting of aerosol scattering leads to retrieval errors of(-0.1%)-(-0.5%),0.22%-1.92%,0.09%-1.46% and 0.02%-0.45% for urban,marine,rural and tropospheric aerosol models respectively.XCO2(the dry air column averaged mixing ratio of CO2) was retrieved from GOSAT observations with in-situ measurements of aerosol properties as input parameters.Comparison with ground-based XCO2measurements shows that the retrieved XCO2(390.95ppm)is more consistent with the measured XCO2(390.737ppm)as compared to that of GOSAT L2 product(389.814ppm).

    Keywords:retrieval of XCO2,near infrared bands,GOSAT,aerosols

    收稿日期:2015年8月31日;修回日期:2015年12月29日

    DOI:10.3969/j.issn.2095-1973.2016.01.016

    猜你喜歡
    氣溶膠
    基于AERONET的多源多尺度衛(wèi)星氣溶膠產(chǎn)品真實性檢驗
    基于飛機(jī)觀測的四川盆地9月氣溶膠粒子譜分析
    氣溶膠傳播之謎
    CF-901型放射性氣溶膠取樣泵計算公式修正
    氣溶膠中210Po測定的不確定度評定
    氣溶膠指數(shù)與臭氧總量的相關(guān)性初步研究
    Aerosol absorption optical depth of fine-mode mineral dust in eastern China
    四川盆地秋季氣溶膠與云的相關(guān)分析
    氣溶膠對雷暴云起電以及閃電發(fā)生率影響的數(shù)值模擬
    氣溶膠科學(xué)
    中文字幕最新亚洲高清| 欧美激情国产日韩精品一区| 精品人妻熟女毛片av久久网站| 久久久久久人妻| 日韩亚洲欧美综合| 在线 av 中文字幕| 亚洲人与动物交配视频| 秋霞在线观看毛片| 国产精品国产三级国产专区5o| 在现免费观看毛片| 欧美精品亚洲一区二区| 一级毛片我不卡| 建设人人有责人人尽责人人享有的| 高清毛片免费看| 欧美激情 高清一区二区三区| 亚洲成色77777| 亚洲天堂av无毛| 18禁在线播放成人免费| 婷婷色av中文字幕| 国产成人一区二区在线| 免费观看a级毛片全部| 99九九线精品视频在线观看视频| 国产精品一区www在线观看| 丝袜喷水一区| 亚洲人成77777在线视频| 久久精品久久久久久噜噜老黄| 男女边吃奶边做爰视频| 少妇 在线观看| 成人综合一区亚洲| 精品人妻偷拍中文字幕| 搡老乐熟女国产| 欧美精品亚洲一区二区| 尾随美女入室| 成人手机av| 国产高清国产精品国产三级| 欧美丝袜亚洲另类| 国产在线免费精品| 99热国产这里只有精品6| 黄色怎么调成土黄色| 国产日韩欧美视频二区| 精品人妻熟女毛片av久久网站| 日韩av不卡免费在线播放| 精品久久久久久久久av| 九色亚洲精品在线播放| 高清视频免费观看一区二区| 欧美bdsm另类| 亚洲成人手机| 精品亚洲成a人片在线观看| 看十八女毛片水多多多| 视频区图区小说| 国产av国产精品国产| 免费黄色在线免费观看| av黄色大香蕉| 在线观看人妻少妇| 午夜久久久在线观看| 高清av免费在线| 国产综合精华液| 久久国产亚洲av麻豆专区| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久久免| 高清欧美精品videossex| 国产精品99久久99久久久不卡 | 麻豆精品久久久久久蜜桃| 青春草国产在线视频| 亚洲精品乱码久久久v下载方式| 日本91视频免费播放| 一个人免费看片子| 免费看不卡的av| 性高湖久久久久久久久免费观看| 久久精品国产鲁丝片午夜精品| 国产精品久久久久成人av| 精品人妻熟女av久视频| 卡戴珊不雅视频在线播放| 91精品国产国语对白视频| 欧美性感艳星| av国产久精品久网站免费入址| 久久国产亚洲av麻豆专区| 国产片特级美女逼逼视频| 美女视频免费永久观看网站| 国产乱人偷精品视频| 精品一区二区三卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人国语在线视频| 色哟哟·www| 亚洲国产精品一区三区| 日韩欧美精品免费久久| av卡一久久| 人妻系列 视频| 如何舔出高潮| 午夜av观看不卡| 一区在线观看完整版| 国产精品国产三级国产av玫瑰| 99热国产这里只有精品6| 国产亚洲欧美精品永久| 久久青草综合色| 99久国产av精品国产电影| 日本黄色片子视频| 一级爰片在线观看| 日韩成人伦理影院| 欧美精品一区二区大全| 国产一区亚洲一区在线观看| 高清av免费在线| 精品一区在线观看国产| 91在线精品国自产拍蜜月| 美女福利国产在线| 亚洲精品国产av成人精品| 美女cb高潮喷水在线观看| 亚洲精品中文字幕在线视频| 成年女人在线观看亚洲视频| 伦理电影大哥的女人| 亚洲av中文av极速乱| 亚洲av男天堂| av天堂久久9| 亚洲精品av麻豆狂野| 老司机影院成人| 国产乱来视频区| 欧美日韩亚洲高清精品| 亚洲少妇的诱惑av| 熟女av电影| 亚洲第一区二区三区不卡| 亚洲色图 男人天堂 中文字幕 | 三上悠亚av全集在线观看| 中文字幕制服av| 国产精品久久久久久久久免| 国产69精品久久久久777片| 久久精品夜色国产| 丝袜在线中文字幕| 黄色欧美视频在线观看| 久久久欧美国产精品| 日日摸夜夜添夜夜添av毛片| 日日摸夜夜添夜夜添av毛片| 亚洲精品一二三| 99国产精品免费福利视频| 精品少妇黑人巨大在线播放| 精品人妻熟女毛片av久久网站| videos熟女内射| 国产高清三级在线| 成人综合一区亚洲| 成人午夜精彩视频在线观看| 欧美老熟妇乱子伦牲交| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单| 五月开心婷婷网| 久久久久久久久久久免费av| 亚洲精品美女久久av网站| 亚洲久久久国产精品| 18+在线观看网站| 午夜视频国产福利| 国产午夜精品一二区理论片| 久热这里只有精品99| 国产av一区二区精品久久| 久久狼人影院| tube8黄色片| 精品亚洲成a人片在线观看| 日日啪夜夜爽| 精品亚洲成国产av| 亚洲欧美成人精品一区二区| 日本黄色日本黄色录像| 看十八女毛片水多多多| 99精国产麻豆久久婷婷| 亚洲经典国产精华液单| 欧美xxxx性猛交bbbb| 亚洲精品乱码久久久v下载方式| 国产乱人偷精品视频| 精品少妇黑人巨大在线播放| 亚洲精品,欧美精品| 精品久久久久久电影网| 日韩欧美一区视频在线观看| 成年美女黄网站色视频大全免费 | 一级毛片电影观看| 制服人妻中文乱码| 成人无遮挡网站| 高清毛片免费看| 久久久国产精品麻豆| 99热国产这里只有精品6| 卡戴珊不雅视频在线播放| 欧美亚洲日本最大视频资源| 欧美成人精品欧美一级黄| 国产又色又爽无遮挡免| 母亲3免费完整高清在线观看 | 亚洲人与动物交配视频| 少妇精品久久久久久久| 伦理电影大哥的女人| 精品国产一区二区久久| 性色avwww在线观看| av黄色大香蕉| 亚洲av二区三区四区| 久久久久久伊人网av| 国产免费一区二区三区四区乱码| 另类精品久久| av专区在线播放| 人人妻人人澡人人爽人人夜夜| 婷婷色麻豆天堂久久| 黑人欧美特级aaaaaa片| 亚洲人与动物交配视频| 国产淫语在线视频| 欧美精品一区二区大全| 十分钟在线观看高清视频www| 黄片播放在线免费| 国精品久久久久久国模美| 免费看光身美女| 久久精品人人爽人人爽视色| 只有这里有精品99| √禁漫天堂资源中文www| 国内精品宾馆在线| 亚洲av成人精品一区久久| 色网站视频免费| 亚洲欧美成人综合另类久久久| 成人漫画全彩无遮挡| 久久人人爽人人片av| 午夜精品国产一区二区电影| 人妻制服诱惑在线中文字幕| 精品国产乱码久久久久久小说| 男男h啪啪无遮挡| 老女人水多毛片| 男的添女的下面高潮视频| 欧美国产精品一级二级三级| 日韩一区二区三区影片| av视频免费观看在线观看| 日本av免费视频播放| 国产黄片视频在线免费观看| 免费看av在线观看网站| 99国产精品免费福利视频| 亚洲欧美一区二区三区国产| 亚洲精品中文字幕在线视频| 亚洲在久久综合| 日韩不卡一区二区三区视频在线| 亚洲av电影在线观看一区二区三区| 99久久精品国产国产毛片| 成人综合一区亚洲| 国产极品天堂在线| 亚洲精品av麻豆狂野| 亚洲精品一区蜜桃| 日产精品乱码卡一卡2卡三| 一本久久精品| 成年人免费黄色播放视频| 18禁裸乳无遮挡动漫免费视频| 精品少妇久久久久久888优播| 成年美女黄网站色视频大全免费 | 超色免费av| 大码成人一级视频| 国产精品偷伦视频观看了| 桃花免费在线播放| 狂野欧美激情性xxxx在线观看| 欧美丝袜亚洲另类| 国产精品久久久久久久电影| 免费观看无遮挡的男女| 91精品一卡2卡3卡4卡| 婷婷色综合大香蕉| 国产精品99久久99久久久不卡 | 欧美另类一区| 久久久久久久亚洲中文字幕| 欧美日韩精品成人综合77777| 一本久久精品| 欧美激情极品国产一区二区三区 | 国产免费一区二区三区四区乱码| 热99久久久久精品小说推荐| 一本色道久久久久久精品综合| 久久精品夜色国产| 在线看a的网站| 校园人妻丝袜中文字幕| 日韩欧美一区视频在线观看| 一级毛片 在线播放| 国产69精品久久久久777片| 免费观看在线日韩| 一二三四中文在线观看免费高清| 80岁老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 久久99热这里只频精品6学生| 99热这里只有是精品在线观看| 九草在线视频观看| 黑人欧美特级aaaaaa片| 精品久久久噜噜| 国产av码专区亚洲av| 久久久久视频综合| 久热久热在线精品观看| av有码第一页| 少妇的逼好多水| 狠狠精品人妻久久久久久综合| 男女免费视频国产| av卡一久久| 97超碰精品成人国产| 亚洲av男天堂| 中文字幕制服av| 男女边摸边吃奶| 最后的刺客免费高清国语| 国产伦理片在线播放av一区| 国国产精品蜜臀av免费| 日本欧美视频一区| 久久精品久久久久久噜噜老黄| 免费日韩欧美在线观看| 天天影视国产精品| 亚洲国产av影院在线观看| 99热这里只有是精品在线观看| 亚洲经典国产精华液单| 天美传媒精品一区二区| 亚洲精华国产精华液的使用体验| 最近最新中文字幕免费大全7| 麻豆精品久久久久久蜜桃| 综合色丁香网| 久热这里只有精品99| 极品人妻少妇av视频| 亚洲欧美中文字幕日韩二区| 精品熟女少妇av免费看| 美女内射精品一级片tv| 少妇高潮的动态图| √禁漫天堂资源中文www| 精品亚洲成a人片在线观看| 如日韩欧美国产精品一区二区三区 | 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕 | 精品午夜福利在线看| 亚洲美女搞黄在线观看| 少妇丰满av| 久久久精品区二区三区| 中文欧美无线码| 97精品久久久久久久久久精品| 色网站视频免费| 久久 成人 亚洲| 国产伦理片在线播放av一区| 黑人欧美特级aaaaaa片| 国产成人一区二区在线| 激情五月婷婷亚洲| 久久久精品免费免费高清| 精品午夜福利在线看| 精品亚洲成国产av| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 91国产中文字幕| 久久久久国产精品人妻一区二区| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 99国产综合亚洲精品| 热re99久久国产66热| 日韩,欧美,国产一区二区三区| 多毛熟女@视频| 一级爰片在线观看| 啦啦啦在线观看免费高清www| 国产深夜福利视频在线观看| 少妇人妻 视频| 91精品国产九色| 久久99蜜桃精品久久| 99热全是精品| 国产精品成人在线| 美女视频免费永久观看网站| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 91午夜精品亚洲一区二区三区| 夜夜骑夜夜射夜夜干| 午夜免费男女啪啪视频观看| 国产国语露脸激情在线看| 九色亚洲精品在线播放| 成人亚洲精品一区在线观看| av天堂久久9| 亚洲av成人精品一区久久| 五月伊人婷婷丁香| 一本久久精品| 亚洲,欧美,日韩| 日韩三级伦理在线观看| 18禁动态无遮挡网站| 人妻夜夜爽99麻豆av| 婷婷色综合www| 黑丝袜美女国产一区| 精品卡一卡二卡四卡免费| 五月开心婷婷网| 黄色毛片三级朝国网站| 成人免费观看视频高清| videossex国产| 免费高清在线观看日韩| 中文字幕精品免费在线观看视频 | 亚洲av成人精品一二三区| 亚洲人成网站在线播| av免费观看日本| 欧美 日韩 精品 国产| 欧美日韩在线观看h| 青春草国产在线视频| 最新的欧美精品一区二区| 狂野欧美激情性bbbbbb| 18在线观看网站| 久久精品国产亚洲网站| 日本av免费视频播放| 成人二区视频| 日韩不卡一区二区三区视频在线| 欧美成人午夜免费资源| 久久久久精品久久久久真实原创| 久久久久精品性色| 国产男女超爽视频在线观看| 欧美性感艳星| av在线app专区| 18在线观看网站| 亚洲精品久久成人aⅴ小说 | 国产精品国产三级国产av玫瑰| 国产亚洲最大av| 妹子高潮喷水视频| 永久网站在线| 免费观看性生交大片5| 国产高清三级在线| 国产一区二区三区av在线| 一本大道久久a久久精品| 国产伦理片在线播放av一区| 国产精品无大码| 卡戴珊不雅视频在线播放| 成人午夜精彩视频在线观看| 七月丁香在线播放| 丰满饥渴人妻一区二区三| 午夜福利在线观看免费完整高清在| 9色porny在线观看| 国产一区有黄有色的免费视频| 大香蕉久久成人网| 一本一本综合久久| 亚洲经典国产精华液单| 欧美性感艳星| 精品一品国产午夜福利视频| 少妇人妻 视频| 建设人人有责人人尽责人人享有的| 天堂中文最新版在线下载| 五月天丁香电影| 亚洲av欧美aⅴ国产| 日韩欧美精品免费久久| 五月开心婷婷网| 青春草亚洲视频在线观看| 免费观看无遮挡的男女| 中国美白少妇内射xxxbb| 国产一区亚洲一区在线观看| 一级爰片在线观看| 精品国产乱码久久久久久小说| 亚洲精品成人av观看孕妇| 永久免费av网站大全| 男人爽女人下面视频在线观看| 99国产精品免费福利视频| 国产视频首页在线观看| 久久国产精品男人的天堂亚洲 | 汤姆久久久久久久影院中文字幕| 少妇猛男粗大的猛烈进出视频| 免费看不卡的av| 青青草视频在线视频观看| 国产在线免费精品| 国产在视频线精品| 日韩av不卡免费在线播放| 在线免费观看不下载黄p国产| 午夜日本视频在线| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 亚洲图色成人| 一级片'在线观看视频| 80岁老熟妇乱子伦牲交| av在线播放精品| 久久久久久久久久人人人人人人| 久久国内精品自在自线图片| 亚洲精品日韩在线中文字幕| 午夜精品国产一区二区电影| 五月玫瑰六月丁香| 尾随美女入室| 一区二区三区免费毛片| 国产伦理片在线播放av一区| 久久综合国产亚洲精品| 人人妻人人澡人人爽人人夜夜| 丝袜美足系列| 国产有黄有色有爽视频| 狂野欧美激情性bbbbbb| 精品国产乱码久久久久久小说| 自线自在国产av| 亚洲精品色激情综合| 午夜视频国产福利| 美女xxoo啪啪120秒动态图| 久久久国产精品麻豆| 一本大道久久a久久精品| 日韩一本色道免费dvd| 亚洲国产欧美日韩在线播放| 国产成人免费观看mmmm| 精品人妻在线不人妻| 日本av手机在线免费观看| 国产男女超爽视频在线观看| 九色亚洲精品在线播放| 国产一区二区在线观看日韩| 久久久久久久亚洲中文字幕| 丰满迷人的少妇在线观看| 麻豆成人av视频| 国产精品久久久久久av不卡| 成人毛片60女人毛片免费| 久久精品久久久久久久性| 日韩制服骚丝袜av| h视频一区二区三区| 天美传媒精品一区二区| 精品熟女少妇av免费看| 日本黄色片子视频| 亚洲伊人久久精品综合| 国产av一区二区精品久久| 一区二区三区四区激情视频| 制服丝袜香蕉在线| 亚洲av综合色区一区| 国产成人免费观看mmmm| 免费看不卡的av| av卡一久久| 中文天堂在线官网| 在线亚洲精品国产二区图片欧美 | 欧美变态另类bdsm刘玥| 中文字幕av电影在线播放| 亚洲欧美精品自产自拍| 亚洲av综合色区一区| 桃花免费在线播放| 欧美激情 高清一区二区三区| 国产av国产精品国产| 日本免费在线观看一区| 麻豆精品久久久久久蜜桃| 草草在线视频免费看| 一级毛片aaaaaa免费看小| 午夜免费鲁丝| 久久人妻熟女aⅴ| 久久99热6这里只有精品| 九九久久精品国产亚洲av麻豆| 天天影视国产精品| 午夜免费鲁丝| 国产免费福利视频在线观看| 欧美亚洲日本最大视频资源| 久久久久视频综合| 亚洲在久久综合| 国产精品久久久久成人av| 国产黄频视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产男人的电影天堂91| 日本黄色日本黄色录像| 在线 av 中文字幕| 日本午夜av视频| 久久久久久久亚洲中文字幕| 黄色怎么调成土黄色| 国产精品一区二区在线观看99| 中文乱码字字幕精品一区二区三区| 亚洲精品自拍成人| 日日啪夜夜爽| 亚洲av电影在线观看一区二区三区| 亚洲中文av在线| 午夜激情av网站| 国产国语露脸激情在线看| 亚洲丝袜综合中文字幕| 高清不卡的av网站| 精品人妻偷拍中文字幕| 亚洲精品国产av蜜桃| 久久人人爽人人爽人人片va| 黄色一级大片看看| 简卡轻食公司| 亚洲欧洲国产日韩| 久久久精品免费免费高清| 热re99久久国产66热| 久久免费观看电影| 久久精品国产a三级三级三级| 免费久久久久久久精品成人欧美视频 | 丝瓜视频免费看黄片| 一边亲一边摸免费视频| 伦理电影大哥的女人| 熟女人妻精品中文字幕| 日韩精品有码人妻一区| 在线观看免费视频网站a站| 久久国产亚洲av麻豆专区| 久久免费观看电影| 97在线人人人人妻| 国产综合精华液| 亚洲欧美精品自产自拍| 在线亚洲精品国产二区图片欧美 | 99国产精品免费福利视频| 一本一本综合久久| 精品国产乱码久久久久久小说| 欧美一级a爱片免费观看看| 免费观看性生交大片5| 亚洲av欧美aⅴ国产| 久久精品国产亚洲av天美| 国产片内射在线| 国产有黄有色有爽视频| 欧美亚洲日本最大视频资源| 国产国拍精品亚洲av在线观看| 亚洲国产日韩一区二区| 亚洲欧美一区二区三区国产| 久久久久久久久大av| 日韩成人伦理影院| 大香蕉97超碰在线| 3wmmmm亚洲av在线观看| 考比视频在线观看| 日韩免费高清中文字幕av| 中文欧美无线码| 美女中出高潮动态图| 国精品久久久久久国模美| 一级,二级,三级黄色视频| 精品亚洲成国产av| 久久精品久久精品一区二区三区| 免费观看在线日韩| 三级国产精品片| 国产一区亚洲一区在线观看| 岛国毛片在线播放| 99热全是精品| 精品午夜福利在线看| 女人久久www免费人成看片| 高清欧美精品videossex| 亚洲av国产av综合av卡| 晚上一个人看的免费电影| 国产精品国产三级国产av玫瑰| 日韩电影二区| 国产综合精华液| 不卡视频在线观看欧美| 18在线观看网站| √禁漫天堂资源中文www| av有码第一页| 国产精品一区www在线观看| 99热国产这里只有精品6| 国产精品麻豆人妻色哟哟久久| 国产一区二区在线观看日韩| av线在线观看网站| 国产男人的电影天堂91| 精品人妻熟女毛片av久久网站| 午夜免费鲁丝| 肉色欧美久久久久久久蜜桃| 国产精品国产三级国产专区5o| 少妇被粗大猛烈的视频| 国产亚洲精品第一综合不卡 | 日本色播在线视频| 亚洲av在线观看美女高潮|