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

    基于壓差的二氧化碳垂直采樣分析技術(shù)

    2019-10-23 12:06:44方雙喜馬千里郭敏銳
    中國環(huán)境科學(xué) 2019年10期
    關(guān)鍵詞:標(biāo)氣長管廓線

    梁 苗,方雙喜*,劉 毅,姚 波,王 勇,馬千里,易 游,郭敏銳

    基于壓差的二氧化碳垂直采樣分析技術(shù)

    梁 苗1,方雙喜1*,劉 毅2,姚 波1,王 勇2,馬千里3,易 游2,郭敏銳4

    (1.中國氣象局氣象探測中心,北京 100081;2.中國科學(xué)院大氣物理研究所,北京 100029;3.中國氣象局浙江臨安大氣本底站,浙江 臨安 311307;4.北京師范大學(xué)全球變化與地球系統(tǒng)科學(xué)研究院,北京 100875)

    本研究開發(fā)了一套基于直接采樣技術(shù)的二氧化碳(CO2)垂直廓線采樣分析系統(tǒng),分析近地面至25km高空CO2濃度的垂直分布.系統(tǒng)通過壓差實現(xiàn)垂直方向連續(xù)采樣,利用CRDS高精度分析技術(shù)對不同高度樣品進行定量分析,計算得到采樣區(qū)域CO2濃度廓線.于2018年6月13~14日,在內(nèi)蒙古錫林浩特國家氣候觀象臺利用平流層探空氣球平臺進行了觀測實驗.實驗室測試顯示,CO2分析準(zhǔn)確度優(yōu)于0.06×10-6,精度優(yōu)于0.08×10-6.外場實驗獲得區(qū)域近地面至25km高空CO2濃度的高分辨率垂直廓線,顯示CO2在不同高度的分層結(jié)構(gòu).考慮不同高度樣品擴散作用,系統(tǒng)垂直分辨率在10km高度以下優(yōu)于580m,在10~20km高度優(yōu)于3.3km.研究表明: 分析系統(tǒng)可搭載在合適探空平臺上進行CO2垂直觀測,獲取濃度廓線,可為傳輸模式提供數(shù)據(jù),并為碳衛(wèi)星遙感數(shù)據(jù)提供實測數(shù)據(jù)校驗.

    垂直采樣;二氧化碳;濃度廓線;垂直分辨率;平流層

    大氣二氧化碳(CO2)等溫室氣體濃度增長與全球變暖密切相關(guān),一般認為其濃度大幅增長是全球變暖的主因[1].由于排放源強、擴散條件、氣候因素等多種差異,大氣中溫室氣體濃度存在明顯的時空變化.目前全球范圍內(nèi)已建立了很多近地面網(wǎng)點用于大氣CO2濃度觀測,而垂直尺度的觀測相對較少[2-3].CO2濃度的垂直觀測可為認識CO2通量(源和匯)分布提供重要科學(xué)依據(jù)[4-5].目前國際上主要利用衛(wèi)星遙感技術(shù)和地基光譜觀測技術(shù)進行溫室氣體柱總量或廓線觀測,并利用模式反演出地面CO2通量分布.這一“自上而下”方法,對于定量監(jiān)測和評估人為活動及生態(tài)系統(tǒng)對CO2濃度和源匯總量的影響,并預(yù)測CO2變化趨勢至關(guān)重要[6-7].衛(wèi)星遙感技術(shù)可以提供CO2等溫室氣體的全球尺度柱總量,此類衛(wèi)星有中國碳衛(wèi)星(TanSat)、美國軌道碳觀測衛(wèi)星(Orbiting Carbon Observatory-2,OCO-2)[8]、日本溫室氣體觀測衛(wèi)星(Greenhouse Gases Observing Satellite, GOSAT)[9]、歐洲ENVISAT 衛(wèi)星上搭載的大氣吸收光譜掃描成像儀(Scanning Imaging Absorption Spectrometer for Atmospheric Cartography,SCIAMACHY)[10]等.此外,基于地基的遙感觀測也是獲取大氣溫室氣體廓線或總量的一種重要手段,國際上最具代表性的地基光譜儀器為全球總碳柱觀測網(wǎng)(TCCON)[11]廣泛使用的傅里葉變換光譜儀(FTS).

    利用原位觀測的CO2垂直廓線數(shù)據(jù),是驗證星載和地基光譜觀測反演結(jié)果的重要依據(jù)[12],同時還能減少區(qū)域碳源、匯反演模式的不確定度[11],提升區(qū)域尺度碳通量反演結(jié)果的準(zhǔn)確度.本研究利用基于大氣CO2廓線直接采樣的觀測技術(shù)(“氣芯”)[13],設(shè)計并組裝了一套下投式溫室氣體垂直廓線采樣分析系統(tǒng),可獲得地面直至下平流層(約25km)溫室氣體濃度廓線.利用此技術(shù)在內(nèi)蒙古錫林浩特國家氣候觀象臺進行了垂直觀測,初步得到了地面直至25km大氣CO2濃度垂直廓線.

    1 材料與方法

    1.1 系統(tǒng)構(gòu)造

    為了開展CO2垂直廓線探測,本文在前期相關(guān)研究的基礎(chǔ)上[7],研制了一套溫室氣體垂直廓線采樣分析系統(tǒng).系統(tǒng)由探空采樣系統(tǒng)和地面分析系統(tǒng)組成.

    1.1.1 探空采樣系統(tǒng) 包括盤狀長管采樣器、氣球探空平臺和地面控制模塊.其中長管采樣器由不銹鋼細長管、填充高氯酸鎂的不銹鋼干燥管(外徑10mm,長度8mm)、熔斷式自動關(guān)閉閥和截止閥構(gòu)成.長管由兩段30m長、外徑分別為6.35和3.18mm的管線(Restek,Inc.,美國)連接而成,管壁厚度分別為0.89和0.71mm,內(nèi)部容積約為565mL,重量約5kg,實現(xiàn)降落過程大氣樣品的采集和儲存.采樣管在使用之前持續(xù)通入干燥、高溫(100℃)高純氮氣(99.999%,氦普北分)12h,使管內(nèi)壁達到鈍化平衡狀態(tài).

    氣球探空平臺包括1000m3零壓氣球(填充氦氣)、切割器、降落傘和集成探空模塊.其中集成探空模塊由氣壓、氣溫和濕度傳感器、銥星定位系統(tǒng)(JT007,中國)、數(shù)據(jù)/信號收發(fā)裝置,以及數(shù)據(jù)處理單元構(gòu)成.集成探空模塊和長管采樣器共同包裹在塑料泡沫中作為載荷,保證實驗過程長管采樣器的溫度變化不超過20K,整體載荷重量約17kg.地面控制模塊包括數(shù)傳終端和控制單元.通過地面控制模塊與探空模塊間信號傳輸,實現(xiàn)采樣過程中環(huán)境大氣的大氣溫度、壓力、濕度數(shù)據(jù)及經(jīng)緯度、高度數(shù)據(jù)接收,以及載荷受控下降、定位和采樣管進樣口的自動閉合.

    1.1.2 地面分析平臺 如圖1,包括6口選擇進樣閥(Vici Inc.,美國)、CRDS分析主機(PicarroInc.,G2301,美國)、質(zhì)量流量控制器(MKS 246C,美國,精度0.1sccm)、氣壓表、三通球閥(Swagelok Inc.,美國)構(gòu)成,各組件通過不銹鋼管線(Swagelok,美國)連接,實現(xiàn)采樣器填充/沖洗、樣品分析和分析主機校準(zhǔn).CRDS分析主機單一標(biāo)氣CO2測量精度(1σ)達到0.05×10-6(0.5Hz).采用標(biāo)準(zhǔn)氣體外標(biāo)法對長管樣品的CO2和CO體積分數(shù)進行定值,CO2量值溯源至World Meteorological Organization (WMO)-X2007標(biāo)準(zhǔn);CO量值溯源至WMO-X2014A標(biāo)準(zhǔn).由于CRDS具有較好的穩(wěn)定性和線性,實驗時采用1瓶標(biāo)氣來定值[13].

    圖1 地面分析平臺結(jié)構(gòu)

    1.2 測量原理

    由于開始降落位置(即飛行最高處)采樣管內(nèi)會留存少量殘余氣體,為確定殘余氣體量,長管采樣器在升空前12h用高濃度CO標(biāo)氣(1971.2×10-9)作為填充氣(Fill Gas,FG)進行沖洗和填充,以區(qū)別于遠離排放源的平流層的CO低濃度大氣,從而準(zhǔn)確標(biāo)記最高點樣品開始分析的時間.氣球升空前,長管采樣器3.18mm端截止閥關(guān)閉,6.35mm端熔斷式開關(guān)閥開啟,采樣器隨氣球探空平臺上升,管內(nèi)FG隨環(huán)境氣壓減少而不斷排出.到達指定高度后(25~30km),地面控制模塊遙控剪切氣球,采樣管隨降落傘緩慢下降,在內(nèi)外壓差作用下進行被動式采樣,使不同高度的大氣依次填充在采樣管的不同位置.同時,數(shù)據(jù)采集和發(fā)送裝置完成傳感器探測數(shù)據(jù)和定位數(shù)據(jù)的實時傳輸,包括管壁溫度、載荷艙內(nèi)溫度、大氣壓力、經(jīng)緯度和高度等.載荷降落地面后,地面控制模塊遙控集成探空模塊中的繼電器,自動熔斷開關(guān)閥關(guān)閉,采樣完成.

    實驗人員通過地面控制系統(tǒng)實時獲取的定位數(shù)據(jù),迅速找回載荷,運送至地面分析平臺進行分析.長管采樣器中氣體樣品的分析過程中,管內(nèi)樣氣在標(biāo)準(zhǔn)氣體(濃度接近實際大氣,防止CO2擴散影響)的推動下,勻速緩慢進入CRDS主機進行分析,得到不同高度樣氣的CO2物質(zhì)的量分數(shù).采集過程中,氣體進入管子的流量小于235cm3/min,保持層流運動,氣體擴散主要為分子擴散,無其它混合作用.在采樣和等待回收的時間段內(nèi)(一般不超過4h),單個分子在擴散作用下沿管路的徑向運動距離(24h約擴散3.2m)與采樣管長度相比,可以忽略不計.說明采樣管可以對依次儲存的各段大氣樣品進行保存,從而對整個樣品的濃度分布進行記錄,并達到一定的垂直分辨率.

    數(shù)據(jù)處理時,需要將分析所得CO2物質(zhì)的量比匹配到相應(yīng)高度上.理想條件下,下降過程每一高度的管內(nèi)外壓力達到平衡.根據(jù)理想氣體狀態(tài)方程,可得時刻長管內(nèi)部的大氣樣品物質(zhì)的量:

    式中:nPT分別為對應(yīng)采樣高度采集所得的氣體物質(zhì)的量、大氣壓力、管壁溫度和管的體積.

    采樣管到達地面時管內(nèi)樣品物質(zhì)的量最大,為

    式中:PT分別為地面大氣壓和落地時刻的管壁溫度.

    分析時,通過質(zhì)量流量控制器精確控制,使長管內(nèi)樣品以恒定的物質(zhì)的量流量進入分析主機,即相同時間分析的樣品物質(zhì)的量一定,可得

    式中:n’t時間內(nèi)分析的樣品物質(zhì)的量分數(shù);t為分析時間,total為分析總時間.

    nn對應(yīng),即可將分析時間與采樣高度的環(huán)境氣壓對應(yīng),從而將CO2體積分數(shù)和采樣高度對應(yīng).

    1.3 測試方法

    1.3.1 單一標(biāo)氣填充測試 為檢驗長管采樣器對樣品分析精度和準(zhǔn)確度的影響,在3種條件下進行多次單一標(biāo)氣分析測試,條件1、2分別為單一標(biāo)氣通過長管直接分析和在長管中儲存12h后進行分析.條件3模擬飛行過程中內(nèi)部氣體排空和樣氣采集過程,長管先填充FG,管內(nèi)抽真空至壓力約47hPa,近似對應(yīng)20km高空氣壓.標(biāo)氣填充流量為100cm3/min (標(biāo)況下),為下降過程中的典型平均采集速率.當(dāng)管內(nèi)達到1atm時,停止填充,關(guān)閉進氣閥,通過分析平臺進行分析.常規(guī)條件下CRDS主機的分析精度為0.05×10-6.條件1~3下CRDS主機的分析精度分別為0.06×10-6、0.08×10-6、0.06×10-6.三種條件下CO2濃度值漂移(相對于標(biāo)氣參比濃度)分別為0.04×10-6、0.05×10-6、0.05×10-6. 3種條件下的氣體分析流量均為 40.0cm3/min(標(biāo)況下).

    1.3.2 分段標(biāo)氣填充測試 為檢驗長管采樣器在儲存和分析過程中,內(nèi)部不同濃度氣體的混合和擴散程度,采用間隔進樣的方式依次用高、低濃度CO2標(biāo)氣進行填充(濃度分別為375.79×10-6和454.54× 10-6).兩種標(biāo)氣首先通過多口進樣閥切換進入CRDS主機分析,確定分段測試“基線”,代表兩種標(biāo)氣在分析系統(tǒng)管路及主機腔體內(nèi)的混合和擴散程度.之后標(biāo)氣通過多口閥切換,以40cm3/min(標(biāo)況下)的流量間隔填充長管,每種標(biāo)氣填充5min,直至管內(nèi)壓力接近一個1atm,關(guān)閉進氣開關(guān)閥,將填充樣品進行分析,分析流量為40.0cm3/min(標(biāo)況下),分析結(jié)果與 “基線”進行對比.

    1.4 實驗過程

    依據(jù)前期氣象預(yù)報,選擇2018年6月13、14日滿足探空氣球升空條件(低層850hPa風(fēng)速小于6m/s)的時間,利用垂直廓線采樣分析系統(tǒng)進行了2次實驗,地點在內(nèi)蒙古錫林浩特國家氣候觀象臺(43°57'N, 116°07'E),實驗各步驟時刻見表1.6月13日,氣球于當(dāng)?shù)貢r間07:30開始釋放,約08:56上升至25.4km高度(25hPa),上升過程管內(nèi)FG在內(nèi)外壓差作用下不斷排出,管溫緩慢下降.在最高位置,氣球載荷平飄約1h的飛行,管內(nèi)氣體交換少,管溫下降變緩.平飄結(jié)束,地面控制單元遙控載荷降落,環(huán)境空氣連續(xù)進入管內(nèi).45min后載荷落地,進樣閥未能及時關(guān)閉,1h后由實驗人員找到采樣管后手動關(guān)閉.這段時間由于管溫稍有升高,導(dǎo)致部分近地面樣品排出.43min后采樣器送至分析平臺進行分析.6月14日,氣球到達最高點25.4km(25hPa)后下降,上升和下降過程與13號類似.即將落地(距地面50m)時地面控制模塊遙控控制探空模塊內(nèi)部繼電器,使自動熔斷開關(guān)閥關(guān)閉.等候分析的時間為1h17min.兩次實驗采樣、尋找及分析環(huán)節(jié)銜接順暢.兩次飛行路徑受氣象條件影響差別較大,如圖2所示.13日和14日的降落地點距釋放地點分別為71km和90km.

    表1 探空實驗以及樣品分析的關(guān)鍵時間

    大氣樣品分析前,分析系統(tǒng)的兩個三通球閥調(diào)節(jié)至 a 路(圖 1),通過標(biāo)定氣(CG)對分析主機進行校準(zhǔn),然后對長管采樣器不同高度層大氣樣品進行定量分析,并繪制CO2濃度廓線.分析時,三通球閥1 切換至b路,6口進樣閥切換高濃度CO推動氣(Push Gas,PG),對連接長管的管路進行沖洗,之后將長管6.35mm端與三通球閥1連接, 3.18mm端與三通球閥2連接;即6.35mm端與進樣閥連接,3.18mm端通過MFC與分析主機連接,高層樣氣先進行分析,近地面樣氣后分析,減少高層樣氣在分析過程中分辨率的損失.沖洗連接管是為了精確判斷采樣最低點的樣氣分析開始時間.之后將三通球閥2切換至b路,同時打開長管兩端的開關(guān)閥,采樣管內(nèi)樣品在PG的推動下進入主機進行分析.通過質(zhì)量流量控制器調(diào)節(jié)分析流量為40cm3/min(標(biāo)況下),以保證濃度廓線的垂直分辨率.通過減壓閥調(diào)節(jié)PG的進氣壓力略高于大氣壓(0.03atm),以減少壓力變化引起的氣體混合和擴散.分析時,FG先進入主機分析,之后是大氣樣品.當(dāng)樣氣全部分析完,PG進入主機并穩(wěn)定一段時間后,多口閥切換CG進入主機再次進行校準(zhǔn),分析結(jié)束.CG與PG/FG的CO2濃度分別為434.88×10-6和380.81×10-6,CO濃度分別為247.67×10-9和1971.27×10-9.

    2 結(jié)果與討論

    2.1 標(biāo)氣填充測試結(jié)果

    表2 分段標(biāo)氣填充測試結(jié)果

    單一標(biāo)氣填充,測試結(jié)果表明,該方法分析長管內(nèi)部儲存樣品的精度(濃度波動的標(biāo)準(zhǔn)偏差)優(yōu)于0.06×10-6(15min,1σ),與主機直接分析所得濃度的偏差優(yōu)于0.04×10-6.保存樣品12h后的分析精度優(yōu)于0.08×10-6,與主機直接分析所得濃度的偏差優(yōu)于0.05×10-6.長管抽真空后填充單一標(biāo)氣的分析精度優(yōu)于0.06×10-6,與主機直接分析所得濃度的偏差優(yōu)于0.05×10-6.說明長管可以準(zhǔn)確地保留樣品成分,準(zhǔn)確度優(yōu)于0.05×10-6,精度優(yōu)于0.08×10-6.

    圖3 分段填充測試結(jié)果

    圖4 采集樣品的分析結(jié)果

    分段填充測試結(jié)果如表2、圖3所示.2種標(biāo)氣在直接切換分析和間隔填充后分析兩種情況下,CO2濃度差值£0.06×10-6,表明長管采樣技術(shù)可以準(zhǔn)確地保留樣品成分.標(biāo)氣直接切換進入主機分析時,徑向混合量(兩種濃度標(biāo)氣混合部分的體積,即濃度過渡段時間×流量)為35.1mL,混合主要發(fā)生在主機之前進樣管路與主機腔體內(nèi).與“基線”相比,標(biāo)氣在長管及干燥管內(nèi)的混合造成較大的徑向混合量87.4mL.以上兩處混合作用會影響垂直分辨率.

    2.2 垂直探測結(jié)果

    圖4為2次實驗中主機分析所得原始數(shù)據(jù).具體分析步驟為: 沖洗氣體(使用PG)、FG、平流層樣氣、低層樣氣和PG依次進入分析系統(tǒng).最高點樣品的分析時間為FG與高層樣氣過渡階段的中間時刻,最低點樣品的分析時間為低層樣氣與PG過渡階段的中間時刻.FG在分析時突然升高,可能由于實驗之前管內(nèi)填充了12h的FG,管內(nèi)壁對FG造成了少量污染,影響了FG的CO濃度,但對CO2濃度沒有影響,2次實驗樣品含水率均低于0.002%,對CO2分析濃度影響可以忽略不計.圖5為數(shù)據(jù)處理后2次實驗的CO2濃度廓線,結(jié)果可反映大氣垂直結(jié)構(gòu)和氣團傳輸影響[16].近地面CO2濃度較低,源于夏季陸地生態(tài)系統(tǒng)吸收作用[17].隨高度增加,CO2濃度增高,并呈現(xiàn)一定波動(約0.3×10-6~3×10-6),在對流層上部達到最高值:13日在384.04hPa氣壓高度處達到409.2×10-6, 14日在347.04hPa氣壓高度處達到406.9×10-6.在平流層,由于上對流層與下平流層交換作用導(dǎo)致CO2濃度不斷降低[18].在廓線最頂端,CO2濃度降至最低, 13日為388.9×10-6(25hPa),14日為384.6×10-6(20hPa).13日實驗,當(dāng)采樣載荷下降到827hPa高度時,由于山脈等障礙物遮擋,定位信號無法直接傳輸至地面接收模塊,轉(zhuǎn)而使用銥星傳輸數(shù)據(jù),而銥星數(shù)傳模塊的數(shù)據(jù)采集頻率為1次/min,造成827hPa高度以下廓線數(shù)據(jù)只有兩個濃度值.14日的實驗在793hPa高度以下出現(xiàn)遮擋,僅有四個濃度值.

    2.3 垂直分辨率

    長管采樣分析技術(shù)的有效垂直分辨率與氣體流動過程中的擴散運動和CRDS主機光腔造成的拖尾作用有關(guān).拖尾效應(yīng)的產(chǎn)生原因如下:CRDS主機響應(yīng)0.5Hz,在110cm3/min(標(biāo)況下)流量下分析,即每3.7cm3(標(biāo)況下)樣氣體積測量一次.光腔體積為35cm3(標(biāo)況下),由于保持在1.87×104Pa壓力、45℃條件下,標(biāo)準(zhǔn)狀況下有效光腔體積(標(biāo)況下光腔內(nèi)氣體體積)近似為6cm3(標(biāo)況下).多次測量的氣體在光腔中發(fā)生混合,會產(chǎn)生拖尾效應(yīng).6cm3(標(biāo)況下)相當(dāng)于長管中采集氣體0.28m,近似于CO2分子擴散10min的距離.這一等效擴散距離遠小于分析之前樣氣在管內(nèi)分子擴散的距離.故在計算垂直分辨率時不考慮分析主機的拖尾效應(yīng).

    Karion等[19]的研究已證明,在長管采樣和分析條件下,管內(nèi)的氣體保持層流狀態(tài),氣體運動是分子擴散和泰勒分散的疊加,徑向分子擴散和泰勒分散共同導(dǎo)致氣體在管內(nèi)的徑向混合[20],可用于表征垂直分辨率.

    有效擴散系數(shù)為分子擴散和泰勒分散系數(shù)的加和[21]:

    而管內(nèi)氣體混合距離可以表示為:(考慮到氣體向兩個方向同時擴散和混合)

    式中:為落地后關(guān)閉采樣閥到分析結(jié)束的時間.

    垂直分辨率可以用壓力表示:

    式中:為采樣管長度;P為地面大氣壓.

    由式(6)可知,在主機腔體拖尾效應(yīng)可忽略、分析延遲時間一定的情況下,長管采樣分析的分辨率僅與管長和管徑有關(guān).對于同一采樣管,由于中高層(平流層底部)的壓力梯度遠小于近地面和中下層(對流層),采集到的中高層樣品量遠小于中下層樣品.為了更好地獲得平流層濃度分布,本研究采用兩段不同管徑連接組成采樣管,使得中高層大氣進入3.18mm一端,有利于獲得更高的平流層大氣CO2分辨率.

    圖6 實驗結(jié)果的垂直分辨率

    3 結(jié)語

    本實驗研發(fā)了下投式溫室氣體垂直廓線采樣分析系統(tǒng),通過兩次外場實驗,證明該系統(tǒng)可以直接采集近地面至平流層(0~25km)大氣樣品并進行高精度分析,分析精度優(yōu)于0.08×10-6,得到高分辨率CO2濃度廓線.

    [1] IPCC. Climate Change 2014:Synthesis report. Contribution of working groups I,II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [R]. IPCC,Geneva,Switzerland,2014.

    [2] Biraud S C, Torn M S, Smith J R, et al. A multi-year record of airborne CO2observations in the US Southern Great Plains [J]. Atmosphric Measurement Techniques, 2013,6:751-763.

    [3] GCOS: Systematic observation requirements for satellite-based data products for climate, Global Climate Observing System-154 [R]. avialable at: http://www.wmo.int/pages/prog/gcos/(last access:15O ctober 2016), 2011.

    [4] Sweeney C, Karion A, Wolter S, et al. Seasonal climatology of CO2across North America from aircraft measurements in the NOAA/ESRL Global Greenhouse Gas Reference Network [J]. Journal of Geophysical Research: Atmospheres, 2015,120(10):5155-5190.

    [5] Zhang H, Chen B, van der Laan-Luijkx, et al. Estimating Asian terrestrial carbon fluxes from CONTRAIL aircraft and surface CO2observations for the period [J]. Atmospheric Chemistry & Physics Discussions, 2013,3,27597-27639.

    [6] Gurney K R, Law R M, Denning A S, et al. Towards robust regional estimates of CO2sources and sinks using atmospheric transport models [J]. Nature, 2002,415:626.

    [7] Peters W, Miller J B, Whitaker J, et al. An ensemble data assimilation system to estimate CO2surface fluxes from atmospheric trace gas observations[J]. Journal of Geophysical Research: Atmospheres, 2005, 110(D24).

    [8] Crisp D, Atlas R M, Breon F M, et al. The Orbiting Carbon Observatory (OCO) mission[J]. Advances in Space Research, 2004, 34:700–709.

    [9] Butz A, Guerlet S, Hasekamp O, et al. Toward accurate CO2and CH4observations from GOSAT [J]. Geophysical Research Letters, 2011, 38:L14812.

    [10] Frankenberg C, Aben I, Bergamaschi P, et al. Global column-averaged methane mixing ratios from 2003 to 2009 as derived from SCIAMACHY: Trends and variability [J]. Journal of Geophysical Research: Biogeosciences, 2011,116:D04302.

    [11] Warneke T R, de Beek1, Buchwitz M, et al. Shipborne solar absorption measurements of CO2, CH4, N2O and CO and comparison with SCIAMACHY WFM-DOAS retrievals[J]. Atmospheric Chemistry & Physics, 2005,5(8):2029-2034.

    [12] Washenfelder R A, Toon G C, Blavier, J F, et al. Carbon dioxide column abundances at the Wisconsin Tall Tower site[J]. Journal of Geophysical Research: Atmospheres, 2006,111(D22),D22305.

    [13] Gerbig C, Lin J C, Wofsy S C, et al. Toward constraining regional-scale fluxes of CO2with atmospheric observations over a continent: 1. Observed spatial variability from airborne platforms [J]. Journal of Geophysical Research: Atmospheres, 2003,108(D24).

    [14] Tans P P, System and method for providing vertical profile measurements of atmospheric gases [P]. U.S. Patent 7597014, 2009.

    [15] Fnag S X, Zhou L X, Zang K P, et al. Measurement of atmospheric CO2mixing ratio by cavity ringdown spectroscopy (CRDS) at the 4background stations in China [J]. Acta Scientiae Circumstantiae, 2011,3:624-629.

    [16] Morguí J A and Rodó X. Physical atmospheric structure and tropospheric mixing information in vertical profiles of atmospheric CO2mixing ratios [J]. Journal of Geophysical Research: Atmospheres, 2010,115(D19).

    [17] Tans P P, Fung P I, Takahashi T. Observational contrains on the global atmospheric CO2budget[J]. Science, 1990,247:1431-1438.

    [18] Nakazawa T, Miyashita K, Aoki S, et al. Temporal and spatial variations of upper tropospheric and lower stratospheric carbon dioxide [J]. Tellus, 2010,43(2):106-117.

    [19] Karion A, Sweeney C, Tans P, et al. AirCore: an innovative atmospheric sampling system [J]. Journal of Atmospheric and Oceanic Technology, 2010,27(11):1839-1853.

    [20] Taylor G. Dispersion of soluble matter in solvent flowing slowly through a tube [J]. Proceedings of the Royal Society: A, 1953, 219(1137):186-203.

    [21] Aris R. On the dispersion of a solute in a fluid flowing through a tube [J]. Proceedings of the Royal Society: A, 1999,1:109-120.

    致謝:感謝法蘭克福歌德大學(xué)大氣與環(huán)境科學(xué)研究所(Institute for Atmospheric and Environmental Science, Goethe University Frankfurt)的Andreas Engel博士在分析系統(tǒng)和采樣器設(shè)計方面提供的技術(shù)指導(dǎo).感謝內(nèi)蒙古錫林浩特國家氣候觀象臺業(yè)務(wù)人員對外場實驗的支持.

    A balloon-borne sampler system based on the pressure gradient for vertical profile measurements of CO2.

    LIANG Miao1, FANG Shuang-xi1*, LIU Yi2, YAO Bo1, WANG Yong2, MA Qian-li3, YI You2, GUO Min-rui4

    (1.Meteorological Observation Centre, China Meteorological Administration, Beijing 100081, China;2.Key Laboratory for Middle Atmosphere and Global Environment Observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China;3.Zhejiang Lin’an Regional Atmosphere Watch Station, Lin’an 311307, China;4.College of Global Change and Earth System Science, Beijing Normal University, Beijing 100875, China)., 2019,39(10):4117~4124

    A balloon-borne sampling system for vertical profile of atmospheric CO2from the surface to 25 km was developed and field-tested. The sampler is deployed on a stratospheric balloon and relies on passively collecting air using the atmospheric pressure gradient during descent. The atmospheric sample through continuous sampling is measured by a cavity ring-down spectrometer (CRDS) for CO2upon recovery. Field campaigns were conducted for the first time on June 13thand 14th, 2018 at Xilinhaote Observatory. Measurements of CO2mole fractions in laboratory tests indicated a repeatability of 0.08×10-6and bias to better than 0.06×10-6(1σ) for CO2under various conditions. The vertical resolution with our configuration was determined to be 580m up to 10km and better than 2.0km up to 20km respectively according to the molecular diffusion of the sample. The system provides an in situ method of greenhouse gas measurements for validation of satellite data and estimation of regional flux.

    sampling system;CO2;vertical profiles;vertical resolution;stratosphere

    X830.2

    A

    1000-6923(2019)10-4117-08

    梁 苗(1987-),女,山西運城人,博士,高級工程師,主要研究方向為溫室氣體及相關(guān)微量成分.發(fā)表論文4篇.

    2019-01-24

    科技部國家重點研發(fā)計劃項目(2017YFB0504000);國家自然科學(xué)基金資助項目(41805129 & 41730103)

    * 責(zé)任作者, 研究員, fangsx@cma.gov.cn

    猜你喜歡
    標(biāo)氣長管廓線
    混合標(biāo)氣標(biāo)校溫室氣體數(shù)據(jù)的計算方法與實現(xiàn)
    廣東氣象(2024年1期)2024-03-13 13:49:18
    狀態(tài)監(jiān)測技術(shù)在壓縮天然氣長管拖車在線監(jiān)控上的智慧應(yīng)用研究
    氣體中總烴分析方法的探討
    低溫與特氣(2021年2期)2021-05-21 09:35:36
    不同降水強度下風(fēng)廓線雷達譜矩特征與測風(fēng)準(zhǔn)確性分析
    長管拖車運行風(fēng)險與事故分析和基于風(fēng)險的動態(tài)監(jiān)測
    動態(tài)氣體校準(zhǔn)儀參與大氣監(jiān)測設(shè)備量值溯源過程可行性初探
    同化Argo海洋廓線觀測對CAS-ESM-C的上層海洋溫度模擬的改進
    長管拖車安全應(yīng)用技術(shù)的分析與探究
    基于快速局域線性回歸的IRAS/FY-3B大氣溫濕廓線反演
    GRAPES模式頂外部背景廓線構(gòu)造方法初步研究
    日日干狠狠操夜夜爽| 国产成人福利小说| 大型黄色视频在线免费观看| 一个人看视频在线观看www免费 | 制服丝袜大香蕉在线| 在线观看一区二区三区| 三级男女做爰猛烈吃奶摸视频| 欧美一级a爱片免费观看看| 国产淫片久久久久久久久 | 久久久久性生活片| 91久久精品国产一区二区成人 | 两个人的视频大全免费| 熟妇人妻久久中文字幕3abv| 黄色丝袜av网址大全| 成人av在线播放网站| 精品国产乱码久久久久久男人| 亚洲人成伊人成综合网2020| 亚洲无线在线观看| 少妇人妻一区二区三区视频| 黑人巨大精品欧美一区二区mp4| 麻豆一二三区av精品| 人妻久久中文字幕网| 免费观看的影片在线观看| a级毛片a级免费在线| 国产一区在线观看成人免费| 999久久久精品免费观看国产| 亚洲国产日韩欧美精品在线观看 | 日本一本二区三区精品| 久久精品国产99精品国产亚洲性色| 桃红色精品国产亚洲av| 叶爱在线成人免费视频播放| 久久婷婷人人爽人人干人人爱| 国产三级中文精品| 三级男女做爰猛烈吃奶摸视频| 亚洲成av人片在线播放无| 亚洲avbb在线观看| 久久午夜亚洲精品久久| 国产精品精品国产色婷婷| 久久婷婷人人爽人人干人人爱| 国产精品综合久久久久久久免费| 搡老熟女国产l中国老女人| 亚洲精品乱码久久久v下载方式 | 十八禁人妻一区二区| 国产精品九九99| 搡老岳熟女国产| 夜夜看夜夜爽夜夜摸| 一区二区三区国产精品乱码| 99热这里只有精品一区 | 精品久久久久久久久久久久久| av女优亚洲男人天堂 | 国产激情偷乱视频一区二区| 久久草成人影院| 在线观看免费午夜福利视频| 国产伦人伦偷精品视频| 欧美黑人欧美精品刺激| 男人舔奶头视频| 婷婷精品国产亚洲av| 成人欧美大片| 99久久综合精品五月天人人| 国产精华一区二区三区| 制服丝袜大香蕉在线| 欧美黑人欧美精品刺激| 在线视频色国产色| 夜夜看夜夜爽夜夜摸| 国产成+人综合+亚洲专区| 最近最新中文字幕大全电影3| 很黄的视频免费| 一边摸一边抽搐一进一小说| 欧美性猛交╳xxx乱大交人| 国产欧美日韩一区二区三| 伊人久久大香线蕉亚洲五| 国产精品综合久久久久久久免费| 国产99白浆流出| 亚洲自偷自拍图片 自拍| 久久国产精品人妻蜜桃| 国产美女午夜福利| 亚洲欧美一区二区三区黑人| 久久久久久人人人人人| svipshipincom国产片| 欧美色视频一区免费| 91在线精品国自产拍蜜月 | 国产私拍福利视频在线观看| 日韩免费av在线播放| 一个人看的www免费观看视频| 人人妻人人看人人澡| 精品久久久久久久久久久久久| 97超级碰碰碰精品色视频在线观看| 亚洲五月婷婷丁香| 少妇裸体淫交视频免费看高清| 在线免费观看的www视频| 欧美午夜高清在线| 欧美日本亚洲视频在线播放| 欧美日韩黄片免| 在线十欧美十亚洲十日本专区| 精品国产乱子伦一区二区三区| 美女高潮的动态| 两性夫妻黄色片| 亚洲七黄色美女视频| 国产一区二区激情短视频| 亚洲欧美激情综合另类| 老鸭窝网址在线观看| 成在线人永久免费视频| 18禁国产床啪视频网站| www日本黄色视频网| 一进一出抽搐gif免费好疼| 一二三四在线观看免费中文在| 窝窝影院91人妻| 亚洲第一电影网av| e午夜精品久久久久久久| 国产美女午夜福利| 色播亚洲综合网| 国产精品av久久久久免费| 亚洲熟妇熟女久久| 在线永久观看黄色视频| 国产一区二区在线观看日韩 | 99热这里只有是精品50| svipshipincom国产片| 99久久综合精品五月天人人| 老司机午夜十八禁免费视频| 国产亚洲精品久久久com| 69av精品久久久久久| a级毛片a级免费在线| 国产av不卡久久| 香蕉av资源在线| 亚洲片人在线观看| 精品久久久久久,| 97碰自拍视频| 精品电影一区二区在线| 亚洲一区二区三区色噜噜| 99视频精品全部免费 在线 | 嫩草影院精品99| 亚洲中文av在线| 天堂动漫精品| 精品乱码久久久久久99久播| 男人舔女人下体高潮全视频| 国产男靠女视频免费网站| 免费在线观看影片大全网站| 一级黄色大片毛片| 国产精品国产高清国产av| 两性午夜刺激爽爽歪歪视频在线观看| 中文字幕av在线有码专区| 成人一区二区视频在线观看| 亚洲av成人精品一区久久| 亚洲天堂国产精品一区在线| 欧美日韩亚洲国产一区二区在线观看| 少妇的丰满在线观看| x7x7x7水蜜桃| 精品一区二区三区四区五区乱码| 成人性生交大片免费视频hd| 男人和女人高潮做爰伦理| 久久久久国内视频| 国产精品 国内视频| 18禁美女被吸乳视频| 国内毛片毛片毛片毛片毛片| 麻豆国产97在线/欧美| 久久久久久久久久黄片| 99久久精品热视频| 久久这里只有精品中国| 国产成人系列免费观看| 亚洲真实伦在线观看| 精品人妻1区二区| 亚洲 国产 在线| 午夜视频精品福利| 色哟哟哟哟哟哟| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久久电影 | 99国产精品一区二区三区| 无限看片的www在线观看| 国产精品99久久久久久久久| 国产亚洲精品一区二区www| 久久精品亚洲精品国产色婷小说| 成人永久免费在线观看视频| 超碰成人久久| 曰老女人黄片| 亚洲 欧美 日韩 在线 免费| 色哟哟哟哟哟哟| 婷婷精品国产亚洲av在线| 国产一区二区在线观看日韩 | 在线观看免费视频日本深夜| 国产99白浆流出| 神马国产精品三级电影在线观看| 在线视频色国产色| 亚洲片人在线观看| 亚洲国产欧洲综合997久久,| 成年女人永久免费观看视频| 在线观看美女被高潮喷水网站 | 亚洲精品乱码久久久v下载方式 | 国产精品一区二区免费欧美| 国产欧美日韩精品一区二区| 国产又色又爽无遮挡免费看| 男人舔女人下体高潮全视频| 久久天堂一区二区三区四区| 国产欧美日韩精品亚洲av| 嫩草影院精品99| 最近最新中文字幕大全免费视频| 欧美xxxx黑人xx丫x性爽| 久久午夜亚洲精品久久| 中文字幕精品亚洲无线码一区| 欧美黄色淫秽网站| 国产高清videossex| 国产伦精品一区二区三区四那| 亚洲无线观看免费| 激情在线观看视频在线高清| 成年女人看的毛片在线观看| 18禁黄网站禁片免费观看直播| 亚洲欧美日韩卡通动漫| 最近最新免费中文字幕在线| 精品熟女少妇八av免费久了| 最新在线观看一区二区三区| 蜜桃久久精品国产亚洲av| 真实男女啪啪啪动态图| 亚洲 欧美一区二区三区| 亚洲自偷自拍图片 自拍| 熟女电影av网| 欧美+亚洲+日韩+国产| 国产一区在线观看成人免费| 99精品欧美一区二区三区四区| 淫秽高清视频在线观看| 一个人观看的视频www高清免费观看 | 亚洲国产精品999在线| 丰满的人妻完整版| 日韩三级视频一区二区三区| 日本a在线网址| 青草久久国产| 国产97色在线日韩免费| 亚洲欧洲精品一区二区精品久久久| 美女被艹到高潮喷水动态| 成人国产综合亚洲| 国内精品美女久久久久久| 丁香欧美五月| 长腿黑丝高跟| 91在线精品国自产拍蜜月 | a在线观看视频网站| 无人区码免费观看不卡| 欧美成人免费av一区二区三区| 中文在线观看免费www的网站| 少妇的逼水好多| 婷婷精品国产亚洲av| 大型黄色视频在线免费观看| 少妇丰满av| 国产精品日韩av在线免费观看| 免费搜索国产男女视频| 午夜激情欧美在线| 日韩欧美国产在线观看| 欧美高清成人免费视频www| 日本成人三级电影网站| 亚洲av成人不卡在线观看播放网| 少妇丰满av| 亚洲国产欧美人成| 中亚洲国语对白在线视频| 国产麻豆成人av免费视频| 美女扒开内裤让男人捅视频| 色老头精品视频在线观看| 老熟妇仑乱视频hdxx| 亚洲国产看品久久| 国产一区二区在线av高清观看| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 久久久久精品国产欧美久久久| 午夜日韩欧美国产| 亚洲电影在线观看av| 亚洲欧美日韩无卡精品| 床上黄色一级片| 欧美xxxx黑人xx丫x性爽| 日本黄色片子视频| 国产探花在线观看一区二区| 欧美日韩精品网址| 99国产极品粉嫩在线观看| 91麻豆精品激情在线观看国产| 亚洲性夜色夜夜综合| av视频在线观看入口| 国产成人精品久久二区二区91| 一本久久中文字幕| www.www免费av| 成人三级做爰电影| 中文字幕av在线有码专区| 国产亚洲av高清不卡| 97人妻精品一区二区三区麻豆| 天堂√8在线中文| 国产精品久久久久久人妻精品电影| 久久这里只有精品19| 国产综合懂色| 少妇熟女aⅴ在线视频| 亚洲 欧美 日韩 在线 免费| 美女午夜性视频免费| 噜噜噜噜噜久久久久久91| 男人和女人高潮做爰伦理| 精品99又大又爽又粗少妇毛片 | 曰老女人黄片| 亚洲乱码一区二区免费版| 亚洲精品色激情综合| 亚洲五月天丁香| 久久久久国产一级毛片高清牌| 免费一级毛片在线播放高清视频| 久久香蕉精品热| 久久精品夜夜夜夜夜久久蜜豆| 嫩草影院精品99| 日韩有码中文字幕| 久久热在线av| 动漫黄色视频在线观看| 搡老岳熟女国产| 午夜福利免费观看在线| 1024手机看黄色片| 最新美女视频免费是黄的| 亚洲色图av天堂| 少妇人妻一区二区三区视频| a级毛片a级免费在线| 欧美日韩福利视频一区二区| 亚洲在线自拍视频| 亚洲人成网站在线播放欧美日韩| 99re在线观看精品视频| 最近视频中文字幕2019在线8| 香蕉国产在线看| 最近最新中文字幕大全电影3| 国产乱人伦免费视频| 好男人电影高清在线观看| 天堂√8在线中文| 亚洲成av人片免费观看| 制服人妻中文乱码| 一级黄色大片毛片| 国产伦精品一区二区三区视频9 | 一区福利在线观看| 国内毛片毛片毛片毛片毛片| 国内精品久久久久久久电影| 高潮久久久久久久久久久不卡| 国产精品亚洲一级av第二区| 热99在线观看视频| 午夜亚洲福利在线播放| 欧美色欧美亚洲另类二区| 日日夜夜操网爽| 麻豆久久精品国产亚洲av| 岛国视频午夜一区免费看| 国模一区二区三区四区视频 | 久久欧美精品欧美久久欧美| 老汉色av国产亚洲站长工具| 久久久久九九精品影院| 中文在线观看免费www的网站| 精品不卡国产一区二区三区| 在线永久观看黄色视频| 免费看十八禁软件| 日本三级黄在线观看| 亚洲午夜精品一区,二区,三区| 99re在线观看精品视频| av视频在线观看入口| 久久草成人影院| 成人高潮视频无遮挡免费网站| 亚洲成人中文字幕在线播放| 搡老妇女老女人老熟妇| 中出人妻视频一区二区| 999久久久精品免费观看国产| 一a级毛片在线观看| 精品国产乱码久久久久久男人| 非洲黑人性xxxx精品又粗又长| 啦啦啦观看免费观看视频高清| 精品人妻1区二区| 国产高清有码在线观看视频| 性欧美人与动物交配| 欧美高清成人免费视频www| 久久久久久人人人人人| www.999成人在线观看| 久久久久国内视频| 日本三级黄在线观看| www日本在线高清视频| 在线视频色国产色| 免费观看人在逋| 国产成人欧美在线观看| 亚洲无线观看免费| 国产野战对白在线观看| 久久中文字幕人妻熟女| xxxwww97欧美| 91在线观看av| 日本与韩国留学比较| 国产一区二区三区在线臀色熟女| 美女 人体艺术 gogo| 亚洲av片天天在线观看| 久久久久亚洲av毛片大全| 嫩草影院精品99| 欧美色欧美亚洲另类二区| 99久久99久久久精品蜜桃| 久久久久亚洲av毛片大全| 久久精品综合一区二区三区| 狂野欧美白嫩少妇大欣赏| 一区二区三区国产精品乱码| 伊人久久大香线蕉亚洲五| 九九热线精品视视频播放| 99久久国产精品久久久| 国产成人精品久久二区二区免费| 日韩精品青青久久久久久| 国产精品自产拍在线观看55亚洲| 中文字幕精品亚洲无线码一区| 在线观看舔阴道视频| 久久久国产欧美日韩av| 一夜夜www| 欧美另类亚洲清纯唯美| 婷婷精品国产亚洲av| 午夜免费观看网址| 成年人黄色毛片网站| 亚洲在线观看片| 国内少妇人妻偷人精品xxx网站 | 九色国产91popny在线| 亚洲一区二区三区色噜噜| 19禁男女啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 中文字幕av在线有码专区| 久久人人精品亚洲av| av视频在线观看入口| 国产精品久久电影中文字幕| 一级黄色大片毛片| 五月玫瑰六月丁香| 狂野欧美白嫩少妇大欣赏| 亚洲成av人片在线播放无| 日韩成人在线观看一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲欧洲精品一区二区精品久久久| 亚洲中文字幕日韩| 麻豆国产97在线/欧美| 99久久精品国产亚洲精品| 不卡av一区二区三区| 亚洲美女视频黄频| 国产亚洲av嫩草精品影院| 免费在线观看日本一区| 午夜免费成人在线视频| 狠狠狠狠99中文字幕| 国产综合懂色| 日本免费一区二区三区高清不卡| 听说在线观看完整版免费高清| 最新在线观看一区二区三区| 激情在线观看视频在线高清| 亚洲国产欧美网| 在线观看66精品国产| 黄片小视频在线播放| 精品一区二区三区视频在线 | 国产精品电影一区二区三区| 最好的美女福利视频网| 亚洲精华国产精华精| 啦啦啦韩国在线观看视频| 亚洲自拍偷在线| 99riav亚洲国产免费| 久久久久久久久中文| 国产不卡一卡二| 日本撒尿小便嘘嘘汇集6| 久久精品91蜜桃| 听说在线观看完整版免费高清| 高清在线国产一区| 天天躁狠狠躁夜夜躁狠狠躁| 久9热在线精品视频| 国产免费男女视频| 国产精品日韩av在线免费观看| 黄色丝袜av网址大全| 波多野结衣高清无吗| 桃红色精品国产亚洲av| 天天一区二区日本电影三级| 国产黄a三级三级三级人| 高清在线国产一区| 美女午夜性视频免费| 女警被强在线播放| 亚洲色图av天堂| 观看美女的网站| 欧美最黄视频在线播放免费| 九九久久精品国产亚洲av麻豆 | 不卡av一区二区三区| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 91字幕亚洲| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲无线在线观看| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 后天国语完整版免费观看| 欧美日韩精品网址| 最近视频中文字幕2019在线8| 中文在线观看免费www的网站| 国产精品久久久av美女十八| 无人区码免费观看不卡| 午夜福利在线观看免费完整高清在 | 国产成人福利小说| 亚洲国产精品久久男人天堂| 18禁观看日本| 18禁裸乳无遮挡免费网站照片| 国产毛片a区久久久久| 视频区欧美日本亚洲| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 亚洲欧美日韩高清专用| 老司机深夜福利视频在线观看| 亚洲av日韩精品久久久久久密| 精品久久久久久久久久久久久| 久久久国产欧美日韩av| 精品熟女少妇八av免费久了| 亚洲午夜精品一区,二区,三区| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 一区福利在线观看| 在线看三级毛片| ponron亚洲| 亚洲成av人片免费观看| 精品福利观看| 国内毛片毛片毛片毛片毛片| 国产av麻豆久久久久久久| 色综合亚洲欧美另类图片| 国产精品亚洲av一区麻豆| 午夜福利高清视频| 黄片大片在线免费观看| 国内精品一区二区在线观看| 淫秽高清视频在线观看| 一级毛片女人18水好多| 99久久精品热视频| 欧美中文日本在线观看视频| 国内少妇人妻偷人精品xxx网站 | 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看| a级毛片a级免费在线| 欧美最黄视频在线播放免费| 一二三四社区在线视频社区8| 中文字幕久久专区| 免费av毛片视频| 人妻久久中文字幕网| 男女午夜视频在线观看| 操出白浆在线播放| 动漫黄色视频在线观看| 欧美黑人巨大hd| av福利片在线观看| 午夜福利视频1000在线观看| 国产精品99久久99久久久不卡| 国产一区二区三区在线臀色熟女| 一本综合久久免费| 日本a在线网址| 成人高潮视频无遮挡免费网站| 久久中文字幕一级| 成人精品一区二区免费| 免费大片18禁| 欧美激情久久久久久爽电影| 中亚洲国语对白在线视频| 欧美日韩精品网址| 国模一区二区三区四区视频 | 亚洲av第一区精品v没综合| 色老头精品视频在线观看| 啪啪无遮挡十八禁网站| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区在线臀色熟女| 国产精品电影一区二区三区| 亚洲av免费在线观看| 色尼玛亚洲综合影院| 国产三级在线视频| 免费观看的影片在线观看| 免费看光身美女| 综合色av麻豆| 观看免费一级毛片| 一区福利在线观看| 国产精品亚洲av一区麻豆| 18禁观看日本| 午夜精品在线福利| 久久国产精品影院| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 久久精品国产清高在天天线| 国产男靠女视频免费网站| 久久久久九九精品影院| 欧美激情久久久久久爽电影| 91av网站免费观看| 午夜日韩欧美国产| 免费高清视频大片| 国产综合懂色| 久久这里只有精品中国| 亚洲精品456在线播放app | 亚洲av五月六月丁香网| 51午夜福利影视在线观看| 国产免费男女视频| 黄频高清免费视频| 一个人看视频在线观看www免费 | 天天躁日日操中文字幕| 日韩精品中文字幕看吧| 色播亚洲综合网| 国产精品1区2区在线观看.| 亚洲黑人精品在线| 国语自产精品视频在线第100页| 观看美女的网站| 一个人看视频在线观看www免费 | 久久中文字幕人妻熟女| 欧美乱码精品一区二区三区| 最近视频中文字幕2019在线8| 国产熟女xx| 99国产综合亚洲精品| 亚洲自偷自拍图片 自拍| 国产成人精品久久二区二区91| 岛国在线观看网站| 亚洲五月天丁香| 午夜福利在线观看免费完整高清在 | 国产黄a三级三级三级人| 日日摸夜夜添夜夜添小说| 综合色av麻豆| 亚洲一区二区三区不卡视频| 亚洲欧美一区二区三区黑人| 少妇熟女aⅴ在线视频| cao死你这个sao货| h日本视频在线播放| 中国美女看黄片| 亚洲熟女毛片儿| 精品欧美国产一区二区三| 久久久久久人人人人人| 国产成人一区二区三区免费视频网站| 此物有八面人人有两片| 夜夜夜夜夜久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 人妻久久中文字幕网| 精品国产三级普通话版| 国产高清videossex| 欧美黑人欧美精品刺激| 窝窝影院91人妻| 国产精品精品国产色婷婷| 少妇熟女aⅴ在线视频|