• <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é)
    蜜臀久久99精品久久宅男| 亚洲精品aⅴ在线观看| 亚洲经典国产精华液单| 国产精品一区二区在线观看99| 亚洲国产精品999| 少妇被粗大的猛进出69影院 | 国产亚洲5aaaaa淫片| 国产亚洲5aaaaa淫片| 自拍偷自拍亚洲精品老妇| 又大又黄又爽视频免费| 欧美精品一区二区大全| 黑人高潮一二区| 精品卡一卡二卡四卡免费| a级毛色黄片| 免费少妇av软件| 国产精品一区二区在线观看99| 国产熟女午夜一区二区三区 | 丝袜脚勾引网站| 亚洲精品乱码久久久v下载方式| 欧美高清成人免费视频www| 制服丝袜香蕉在线| 自拍欧美九色日韩亚洲蝌蚪91 | 天天操日日干夜夜撸| 亚洲国产精品专区欧美| av国产久精品久网站免费入址| 亚洲精品456在线播放app| 男女无遮挡免费网站观看| 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| 国产男女内射视频| 水蜜桃什么品种好| 亚洲人成网站在线观看播放| 婷婷色麻豆天堂久久| 又粗又硬又长又爽又黄的视频| 最后的刺客免费高清国语| 三上悠亚av全集在线观看 | 最后的刺客免费高清国语| 狠狠精品人妻久久久久久综合| 国产成人一区二区在线| 亚洲一区二区三区欧美精品| 国产精品一区二区在线不卡| 99热6这里只有精品| 国产淫语在线视频| a级毛片在线看网站| 亚洲国产欧美在线一区| 国产亚洲精品久久久com| 丝瓜视频免费看黄片| 五月开心婷婷网| 久久国产乱子免费精品| 在线观看美女被高潮喷水网站| 亚洲精品第二区| 精品一品国产午夜福利视频| 国内少妇人妻偷人精品xxx网站| 国产一区亚洲一区在线观看| 国产淫片久久久久久久久| tube8黄色片| 99精国产麻豆久久婷婷| 一本色道久久久久久精品综合| 国模一区二区三区四区视频| 亚洲av成人精品一区久久| 国产av精品麻豆| 最新中文字幕久久久久| 午夜91福利影院| 高清毛片免费看| 性高湖久久久久久久久免费观看| 成人国产av品久久久| 少妇人妻精品综合一区二区| 午夜视频国产福利| 亚洲激情五月婷婷啪啪| 偷拍熟女少妇极品色| 国产精品一二三区在线看| 久久毛片免费看一区二区三区| 亚洲精品乱码久久久久久按摩| 高清av免费在线| 国产精品嫩草影院av在线观看| 国产色爽女视频免费观看| 熟女电影av网| 亚洲内射少妇av| 老司机影院成人| 在线观看一区二区三区激情| 在线天堂最新版资源| 国产精品国产三级国产av玫瑰| 18禁在线无遮挡免费观看视频| 午夜精品国产一区二区电影| 五月开心婷婷网| 91久久精品电影网| 欧美精品一区二区免费开放| 亚洲av.av天堂| 国产国拍精品亚洲av在线观看| 丰满迷人的少妇在线观看| 高清不卡的av网站| 欧美精品亚洲一区二区| 超碰97精品在线观看| 久久久久久人妻| 乱码一卡2卡4卡精品| av专区在线播放| 十八禁网站网址无遮挡 | 18禁在线无遮挡免费观看视频| 99久国产av精品国产电影| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久精品古装| 又爽又黄a免费视频| 男女无遮挡免费网站观看| 国产精品99久久久久久久久| 少妇人妻精品综合一区二区| 高清在线视频一区二区三区| 亚洲,欧美,日韩| 久久国产乱子免费精品| 春色校园在线视频观看| 日韩一区二区视频免费看| 久久久久久久久久久久大奶| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满| 国产精品久久久久久久电影| 国产精品免费大片| 少妇熟女欧美另类| 精品少妇黑人巨大在线播放| 99久久中文字幕三级久久日本| 婷婷色综合www| 免费少妇av软件| 国产永久视频网站| 久久午夜综合久久蜜桃| 中文字幕免费在线视频6| 久久女婷五月综合色啪小说| 精品久久久久久久久av| 久久99热6这里只有精品| 在线精品无人区一区二区三| 少妇人妻 视频| 一级毛片我不卡| 亚洲国产毛片av蜜桃av| 街头女战士在线观看网站| 色视频在线一区二区三区| 嘟嘟电影网在线观看| 黄色配什么色好看| 五月玫瑰六月丁香| 免费黄色在线免费观看| 不卡视频在线观看欧美| 十分钟在线观看高清视频www | 性高湖久久久久久久久免费观看| 国产成人精品婷婷| 少妇人妻精品综合一区二区| 精品少妇黑人巨大在线播放| 国产一区亚洲一区在线观看| 免费看av在线观看网站| 中文字幕人妻丝袜制服| 久久综合国产亚洲精品| 亚洲人成网站在线播| 老司机影院毛片| 99久久中文字幕三级久久日本| 久久久久网色| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 国产精品熟女久久久久浪| 一级毛片aaaaaa免费看小| 亚洲精品亚洲一区二区| 91久久精品国产一区二区成人| 在线亚洲精品国产二区图片欧美 | 少妇人妻久久综合中文| 在线 av 中文字幕| 内射极品少妇av片p| 26uuu在线亚洲综合色| 国产免费一区二区三区四区乱码| a级一级毛片免费在线观看| 黄片无遮挡物在线观看| 黄色毛片三级朝国网站 | 多毛熟女@视频| 另类精品久久| 九九在线视频观看精品| 在线观看国产h片| 91精品国产九色| av天堂中文字幕网| 夫妻午夜视频| 熟妇人妻不卡中文字幕| 亚洲一区二区三区欧美精品| 大片免费播放器 马上看| 国产高清三级在线| 2022亚洲国产成人精品| av国产久精品久网站免费入址| 青春草视频在线免费观看| 久久99精品国语久久久| 亚洲婷婷狠狠爱综合网| 日韩欧美一区视频在线观看 | 久久久国产欧美日韩av| 免费高清在线观看视频在线观看| 国产在视频线精品| 久久av网站| 亚洲欧洲国产日韩| 最近2019中文字幕mv第一页| 三级经典国产精品| 男女边吃奶边做爰视频| 亚洲av综合色区一区| 国产黄片美女视频| 97在线视频观看| 嘟嘟电影网在线观看| 国产av码专区亚洲av| 久久久精品94久久精品| 乱人伦中国视频| 人妻夜夜爽99麻豆av| 亚洲情色 制服丝袜| 亚洲国产精品国产精品| 久久久午夜欧美精品| 亚洲色图综合在线观看| 国产一区有黄有色的免费视频| 国产黄色免费在线视频| 丰满少妇做爰视频| 色5月婷婷丁香| 亚洲欧美成人综合另类久久久| 日韩亚洲欧美综合| www.av在线官网国产| 色吧在线观看| 国产精品久久久久久精品古装| 日本免费在线观看一区| 久久 成人 亚洲| 精品一品国产午夜福利视频| av免费在线看不卡| 国产精品伦人一区二区| 美女大奶头黄色视频| 国产av国产精品国产| 精品视频人人做人人爽| 亚洲综合色惰| 久热久热在线精品观看| 日韩人妻高清精品专区| 午夜福利影视在线免费观看| 国产精品.久久久| 国产男女内射视频| 久久久a久久爽久久v久久| 亚洲av国产av综合av卡| 精品久久久精品久久久| 欧美激情极品国产一区二区三区 | 最新中文字幕久久久久| 日本黄色片子视频| 91久久精品电影网| 性色avwww在线观看| 人妻 亚洲 视频| 又大又黄又爽视频免费| 亚洲成人一二三区av| 在线看a的网站| 久久 成人 亚洲| 赤兔流量卡办理| 精品国产一区二区久久| 一级av片app| 成人毛片60女人毛片免费| 91在线精品国自产拍蜜月| 精品一区二区三卡| 精品国产一区二区久久| 91精品国产九色| av网站免费在线观看视频| 高清不卡的av网站| 内射极品少妇av片p| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 亚洲欧美精品专区久久| 亚洲欧美日韩另类电影网站| 久久人妻熟女aⅴ| 嫩草影院新地址| 乱人伦中国视频| 日韩人妻高清精品专区| 一本久久精品| 人人妻人人爽人人添夜夜欢视频 | 18禁动态无遮挡网站| 日韩av不卡免费在线播放| 最近中文字幕高清免费大全6| 久久久久久人妻| 天堂8中文在线网| 欧美日韩国产mv在线观看视频| 99热6这里只有精品| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 久久6这里有精品| 日本91视频免费播放| 蜜臀久久99精品久久宅男| 亚洲高清免费不卡视频| 久久精品国产a三级三级三级| 欧美日韩综合久久久久久| 五月伊人婷婷丁香| 一二三四中文在线观看免费高清| 国产深夜福利视频在线观看| 岛国毛片在线播放| 欧美xxⅹ黑人| 搡老乐熟女国产| 国产爽快片一区二区三区| 精品国产一区二区三区久久久樱花| 久久久国产一区二区| 在线观看www视频免费| 九九在线视频观看精品| 欧美 日韩 精品 国产| 麻豆乱淫一区二区| 国产午夜精品一二区理论片| 一级,二级,三级黄色视频| a级毛片免费高清观看在线播放| 777米奇影视久久| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 亚洲一级一片aⅴ在线观看| 少妇熟女欧美另类| 性高湖久久久久久久久免费观看| 亚洲av男天堂| 久久久久国产精品人妻一区二区| 夫妻性生交免费视频一级片| av免费在线看不卡| 久久人人爽人人片av| 亚洲欧美日韩东京热| 免费久久久久久久精品成人欧美视频 | 一区在线观看完整版| 精品国产国语对白av| 乱码一卡2卡4卡精品| 国产有黄有色有爽视频| 精品久久国产蜜桃| 久久精品夜色国产| 成年人免费黄色播放视频 | 麻豆精品久久久久久蜜桃| 成人亚洲欧美一区二区av| 精品久久国产蜜桃| 日本黄色片子视频| 色94色欧美一区二区| 亚洲av福利一区| 国产精品99久久99久久久不卡 | 大香蕉久久网| 日本vs欧美在线观看视频 | 2022亚洲国产成人精品| 99久久综合免费| 少妇的逼水好多| 18禁裸乳无遮挡动漫免费视频| 亚洲av综合色区一区| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 在线观看免费视频网站a站| 国产午夜精品久久久久久一区二区三区| 91精品国产国语对白视频| 国产av一区二区精品久久| 99热网站在线观看| 免费在线观看成人毛片| 五月开心婷婷网| 日本猛色少妇xxxxx猛交久久| 免费在线观看成人毛片| 国产白丝娇喘喷水9色精品| 视频区图区小说| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 在线观看免费日韩欧美大片 | 如何舔出高潮| 欧美精品人与动牲交sv欧美| 91精品伊人久久大香线蕉| 久久精品熟女亚洲av麻豆精品| 成人黄色视频免费在线看| 亚洲精品日本国产第一区| .国产精品久久| 爱豆传媒免费全集在线观看| 99热全是精品| 少妇高潮的动态图| 亚洲av电影在线观看一区二区三区| 黄色欧美视频在线观看| 国产精品福利在线免费观看| 午夜91福利影院| 国产伦在线观看视频一区| 色哟哟·www| av在线app专区| freevideosex欧美| 嫩草影院入口| 日韩视频在线欧美| 久久久久精品久久久久真实原创| 日本欧美国产在线视频| 免费黄色在线免费观看| 蜜桃在线观看..| 亚洲av.av天堂| 中文字幕免费在线视频6| 欧美精品一区二区免费开放| 在线 av 中文字幕| 女人久久www免费人成看片| 日本猛色少妇xxxxx猛交久久| 日本91视频免费播放| 一区二区av电影网| 嘟嘟电影网在线观看| 亚洲av中文av极速乱| 午夜福利影视在线免费观看| 性高湖久久久久久久久免费观看| 久久人人爽人人爽人人片va| 日韩亚洲欧美综合| 日韩,欧美,国产一区二区三区| 日本vs欧美在线观看视频 | 久久久久久久久久久久大奶| 久久国产亚洲av麻豆专区| 大又大粗又爽又黄少妇毛片口| 成人黄色视频免费在线看| 亚洲精品日本国产第一区| www.av在线官网国产| 国产精品欧美亚洲77777| 精品亚洲乱码少妇综合久久| 亚洲一区二区三区欧美精品| 插阴视频在线观看视频| 97超碰精品成人国产| 777米奇影视久久| 高清毛片免费看| 国产精品成人在线| 亚洲欧美清纯卡通| av视频免费观看在线观看| 日本vs欧美在线观看视频 | 久久99蜜桃精品久久| 欧美丝袜亚洲另类| 观看av在线不卡| 中文字幕亚洲精品专区| 精品国产一区二区久久| 午夜影院在线不卡| 一级毛片电影观看| 午夜免费男女啪啪视频观看| 黄色怎么调成土黄色| 黄色一级大片看看| 一本久久精品| 22中文网久久字幕| 国产精品国产三级国产专区5o| 国产一区二区三区av在线| 在线亚洲精品国产二区图片欧美 | 性色avwww在线观看| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 日产精品乱码卡一卡2卡三| 国产乱来视频区| 成人黄色视频免费在线看| 国产91av在线免费观看| 日日撸夜夜添| h视频一区二区三区| 欧美另类一区| 国产成人freesex在线| 久久韩国三级中文字幕| 成年美女黄网站色视频大全免费 | 国产伦精品一区二区三区四那| 国产高清国产精品国产三级| 久久久久国产网址| 一级毛片我不卡| 亚洲国产精品一区二区三区在线| 秋霞在线观看毛片| av在线观看视频网站免费| 晚上一个人看的免费电影| 亚洲怡红院男人天堂| 亚洲图色成人| 国产 一区精品| 欧美3d第一页| 免费观看的影片在线观看| 国产av一区二区精品久久| 我的女老师完整版在线观看| 亚洲精品中文字幕在线视频 | 在线观看av片永久免费下载| 午夜老司机福利剧场| 国产精品人妻久久久久久| 桃花免费在线播放| 中文字幕免费在线视频6| 久久久午夜欧美精品| 亚洲四区av| 亚洲成人av在线免费| 人人妻人人添人人爽欧美一区卜| 成人综合一区亚洲| 少妇人妻久久综合中文| 菩萨蛮人人尽说江南好唐韦庄| 一级av片app| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久噜噜| 日韩强制内射视频| 成人免费观看视频高清| 中文乱码字字幕精品一区二区三区| 最新的欧美精品一区二区| 18禁动态无遮挡网站| 99久久中文字幕三级久久日本| 午夜免费观看性视频| 国产日韩一区二区三区精品不卡 | 日本-黄色视频高清免费观看| 男人爽女人下面视频在线观看| 免费高清在线观看视频在线观看| 国产欧美日韩一区二区三区在线 | 老司机影院成人| 国产极品天堂在线| 久久久久久人妻| 亚洲国产欧美在线一区| 亚洲欧美成人精品一区二区| 久久 成人 亚洲| 天堂中文最新版在线下载| 狂野欧美激情性bbbbbb| 最新中文字幕久久久久| 久久狼人影院| 人人澡人人妻人| 色94色欧美一区二区| 国产片特级美女逼逼视频| 国产精品福利在线免费观看| 极品少妇高潮喷水抽搐| 国产亚洲最大av| 亚洲国产日韩一区二区| 高清欧美精品videossex| 婷婷色麻豆天堂久久| 亚洲人成网站在线播| 又大又黄又爽视频免费| 如何舔出高潮| 国产精品福利在线免费观看| 免费黄网站久久成人精品| 亚洲国产欧美在线一区| 在线精品无人区一区二区三| 中文字幕精品免费在线观看视频 | 日本黄大片高清| 日本黄色片子视频| 一级毛片 在线播放| 国产精品国产av在线观看| 91精品国产国语对白视频| 三级国产精品欧美在线观看| 国产一级毛片在线| 美女国产视频在线观看| 亚洲精华国产精华液的使用体验| 国产成人a∨麻豆精品| 亚洲欧美精品自产自拍| 一级黄片播放器| 国产精品秋霞免费鲁丝片| h视频一区二区三区| 久久97久久精品| 在线看a的网站| 国产精品久久久久久精品古装| 亚洲一级一片aⅴ在线观看| 亚洲欧美成人精品一区二区| 日韩欧美精品免费久久| 男女边摸边吃奶| 亚洲精品国产av蜜桃| 日韩av免费高清视频| 看免费成人av毛片| 久久久久精品性色| 69精品国产乱码久久久| 中文字幕久久专区| 亚洲,一卡二卡三卡| 永久免费av网站大全| 超碰97精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久成人| 99久国产av精品国产电影| 午夜av观看不卡| 观看免费一级毛片| 观看美女的网站| 永久网站在线| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 国产成人freesex在线| 五月天丁香电影| 美女主播在线视频| h视频一区二区三区| 啦啦啦视频在线资源免费观看| 欧美精品一区二区大全| 欧美成人午夜免费资源| 久久ye,这里只有精品| av有码第一页| 免费黄色在线免费观看| 七月丁香在线播放| 欧美+日韩+精品| 视频中文字幕在线观看| 亚洲第一区二区三区不卡| 亚洲精品第二区| 久久久久久伊人网av| 中国美白少妇内射xxxbb| 婷婷色综合www| 午夜av观看不卡| 美女中出高潮动态图| 日本爱情动作片www.在线观看| 国产精品一区二区性色av| 午夜精品国产一区二区电影| 一级毛片电影观看| 插逼视频在线观看| 一级av片app| 日韩av免费高清视频| 国产白丝娇喘喷水9色精品| 免费黄色在线免费观看| 十分钟在线观看高清视频www | 丰满乱子伦码专区| 人体艺术视频欧美日本| 亚洲av综合色区一区| 大片电影免费在线观看免费| 女性生殖器流出的白浆| 久久人人爽人人片av| 丰满人妻一区二区三区视频av| 99re6热这里在线精品视频| 热99国产精品久久久久久7| 少妇被粗大猛烈的视频| h视频一区二区三区| 王馨瑶露胸无遮挡在线观看| 这个男人来自地球电影免费观看 | 永久免费av网站大全| 九九爱精品视频在线观看| 在线观看国产h片| 制服丝袜香蕉在线| 熟妇人妻不卡中文字幕| 国产日韩一区二区三区精品不卡 | 日本av免费视频播放| 国产淫语在线视频| 精品国产国语对白av| 欧美精品一区二区免费开放| 成年女人在线观看亚洲视频| 精品国产国语对白av| 国产一区二区三区av在线| 精品一品国产午夜福利视频| 亚洲国产最新在线播放| 国产精品成人在线| 亚洲丝袜综合中文字幕| 一级毛片我不卡| 五月伊人婷婷丁香| 免费看不卡的av| 国产精品嫩草影院av在线观看| 少妇人妻 视频| 伦精品一区二区三区| 黄色视频在线播放观看不卡| 久久久久久久久大av| 不卡视频在线观看欧美| 51国产日韩欧美| 女人精品久久久久毛片| 欧美人与善性xxx| 国产精品99久久99久久久不卡 | 夜夜骑夜夜射夜夜干| 亚洲美女黄色视频免费看| 欧美97在线视频| 男男h啪啪无遮挡| 欧美精品一区二区免费开放| av专区在线播放| 一本大道久久a久久精品| 国内揄拍国产精品人妻在线|