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

    武漢市大氣干濕沉降中溶解態(tài)黑碳的季節(jié)分布特征

    2021-02-05 07:50:20白冰冰李思思郭孟如孔少飛
    地球化學 2021年1期
    關鍵詞:大氣

    白冰冰, 洪 軍*, 李思思, 劉 佳, 郭孟如, 黃 穎, 孔少飛

    (中國地質大學(武漢) 環(huán)境學院, 湖北 武漢 430074)

    0 引 言

    黑碳(black carbon, BC)是有機質(生物質及化石燃料等)不完全燃燒的產物[1], 其物理與化學性質都十分穩(wěn)定, 是一種惰性極高的物質, 在全球碳循環(huán)中扮演著非常重要的角色[2-3]。黑碳排放以后, 一部分進入大氣, 另一部分進入土壤。進入大氣和土壤中的黑碳經過大氣干濕沉降和地表徑流等遷移途徑, 最終沉積在河流、土壤、湖泊和海洋等環(huán)境中[4-5]。黑碳因其吸光特性, 在大氣輻射升溫作用中扮演重要角色[6], 黑碳的多孔性和高比表面積使其對有機物有較高的吸附能力, 對環(huán)境中污染物的遷移和降解有重要影響[7]。

    溶解態(tài)黑碳(dissolve black carbon, DBC), 即黑碳的水溶性組分, 其操作上定義為通過固相萃取(solid-phase extraction, SPE)可吸附的黑碳[8], 可以通過浸出和地表徑流從黑碳中釋放出來[9]。溶解態(tài)黑碳在分子水平上為多環(huán)芳族碳的總和[10], 目前普遍認為溶解態(tài)黑碳的分子結構為被部分親水性官能團(羧基、羥基等)取代的濃縮稠環(huán)結構(5~8 個苯環(huán))[8,11,12]。溶解態(tài)黑碳主要來源于土壤黑碳的降解、大氣黑碳的氧化降解及大氣中極細顆粒黑碳的沉降[13-15]。各類淡水生態(tài)系統(tǒng)[10,11,16,17]及海水[8,18]中均檢測到溶解態(tài)黑碳的輸出, 之前的研究顯示, 全球河流輸送到全球海洋的溶解態(tài)黑碳通量每年約(26.5±1.8)百萬噸[9], 這一過程是迄今為止最大的海洋溶解態(tài)黑碳來源。而近期研究顯示, 大氣沉降物也是海洋溶解態(tài)黑碳的重要來源之一, 大氣氣溶膠中的溶解態(tài)黑碳沉降在全球海洋中的通量每年約(1.8±0.8)百萬噸[19]。對河流集水區(qū)溶解態(tài)黑碳來源的研究顯示, 氣溶膠黑碳對溶解態(tài)黑碳河流通量的貢獻為5%~18%[20], 擬合模型表明, 黑碳氣溶膠沉積速率對溶解態(tài)黑碳濃度的影響是土壤黑碳儲量的兩倍[21]。溶解態(tài)黑碳的轉化途徑主要是光降解和生物降解[18,22], 已有研究證明, 河流和海洋表層水體中的溶解態(tài)黑碳在太陽輻射下極易發(fā)生光降解[18,22]。溶解態(tài)黑碳在環(huán)境水體中的遷移轉化會對水質和水生態(tài)環(huán)境產生負面影響。全球超過10 億人依賴于冰川作為水源[23], 而當前研究已顯示, 冰川中含有溶解態(tài)黑碳和其他污染物[24-25], 會使人類的飲用水安全存在潛在危險。溶解態(tài)黑碳光降解過程會產生許多小分子有機物(包括脂肪烴、含C==O 鍵的有機物和濃縮稠環(huán)程度較低的有機物)和活性氧基團(如1O2和O2?等)[18,26,27]。有機物能為微生物的生長提供碳源, 并在消毒過程中產生鹵代甲烷等消毒副產物[28]。

    目前, 對溶解態(tài)黑碳的研究多集中在海洋、河流、冰川以及人工制備的生物質黑碳中, 主要研究其時空分布[17-18], 化學結構特征[10-12]及光化學反應過程[18,22], 而對于大氣干濕沉降中溶解態(tài)黑碳的研究較少。近年來已有針對海洋[19]和森林[29]大氣沉降中溶解態(tài)黑碳的研究, 但對于城市大氣沉降中的溶解態(tài)黑碳及其影響因素還鮮有研究。大氣干濕沉降是大氣中黑碳向陸地及海洋遷移的重要途徑, 研究城市大氣干濕沉降物中溶解態(tài)黑碳的含量及其影響因素對全面理解溶解態(tài)黑碳、完善溶解態(tài)黑碳在全球碳循環(huán)中的作用以及溶解態(tài)黑碳在環(huán)境介質中的遷移轉化規(guī)律具有重要的意義。武漢市區(qū)黑碳氣溶膠的主要來源是燃料燃燒源和以汽車尾氣排放為代表的人為源, 其本底濃度與國內外城市相比處于中等水平, 具有顯著的季節(jié)變化特征[30-31]。武漢市是中國水域面積最大的城市, 黑碳氣溶膠沉降后受納的水體豐沛, 在水體中釋放的溶解態(tài)黑碳會對區(qū)域的水生態(tài)系統(tǒng)造成一定的影響。本研究擬針對武漢市城區(qū)大氣干濕沉降中的溶解態(tài)黑碳定量分析, 獲得溶解態(tài)黑碳濃度與沉降通量, 探討武漢市大氣干濕沉降中溶解態(tài)黑碳的季節(jié)分布特征及其影響因素。

    1 實驗部分

    1.1 儀器、材料與試劑

    研究中主要使用的試劑和化學品有1,2,3-苯三羧酸(1,2,3-B3CA, ≥98%)、1,2,4-苯三羧酸(1,2,4- B3CA, ≥98%)、1,3,5-苯三羧酸(1,3,5-B3CA, ≥ 98%)、1,2,4,5-苯四羧酸(1,2,4,5-B4CA, ≥98%)、苯五羧酸(B5CA, ≥98%)、苯六羧酸(B6CA, ≥98%)、四丁基溴化銨(≥99%)、無水乙酸鈉(99.99%)、甲醇(高效液相色譜級)、鹽酸(36%~38%)和硝酸(65%)。使用的儀器為Waters1500 高效液相色譜儀(美國Waters)、Waters2998 二極管陣列檢測器(美國 Waters)、ZR-3901 型全自動降水采樣器(青島眾瑞, 中國)、固相萃取裝置(上海那艾, 中國)、氮吹儀(崢嶸儀器, 中國)、pH 計(HACH, 美國)、總有機碳分析儀(Hitachi, 日本)。使用的耗材為固相萃取小柱 Bond Elut-PPL (1 g, 6 mL) (Agilent 公司, 美國)和玻璃纖維濾膜(GF/F, 0.7 μm)。

    1.2 實驗方法

    溶解態(tài)黑碳測定方法參照Dittmar[8]利用高效液相色譜(high performance liquid chromatography, HPLC)分析溶解態(tài)黑碳的苯多羧酸法, 其原理是溶解態(tài)黑碳的稠環(huán)結構在濃硝酸高溫高壓條件下氧化消解, 生成單一苯環(huán)的苯多羧酸(benzenepolycarboxylic acids, BPCAs), 通過色譜分析獲得的BPCAs 含量可轉化為溶解態(tài)黑碳含量。苯多羧酸法具有所需樣品量少、根據(jù)BPCAs 含量可獲得一定結構信息等特點, 是目前溶解態(tài)黑碳定量的主流方法。

    1.2.1 樣品采集

    總沉降樣品(干沉降及濕沉降樣品的混合物)的采集主要參照《環(huán)境空氣 降塵的測定 重量法》(GB/T 15265—94), 濕沉降樣品采集參照《大氣降水采樣和分析方法》(GB13580.1—13580.13)。采樣點位于武漢市洪山區(qū)中國地質大學空氣質量自動監(jiān)測站(30°31′13.6″N, 114°23′54.2″E), 屬于典型的城市環(huán)境大氣觀測點。采樣高度距地面約8 m, 總沉降樣品使用6 個直徑15 cm、高30 cm 的玻璃筒于鐵架上露天采集(圖1), 濕沉降樣品使用全自動降水采樣器采集。采樣結束后, 將水樣用20 μm濾膜過濾, ?18 ℃冷凍保存。2018 年3 月至2019 年2 月武漢市氣象數(shù)據(jù)及洪山區(qū)中國地質大學空氣質量監(jiān)測站空氣質量數(shù)據(jù)分別來源于湖北省氣象局(http: //hb.cma. gov.cn/)和武漢市生態(tài)環(huán)境局(http: //hbj.wh.gov.cn/)。

    圖1 總沉降采樣裝置示意圖 Fig.1 Schematic of the total deposition sampling device

    1.2.2 預處理

    樣品預處理步驟如下。

    測定溶解性有機碳 使用總有機碳分析儀分析溶解性有機碳(dissolve organic carbon, DOC)的濃度。

    過濾酸化 根據(jù)每月樣品溶解性有機碳濃度高低將1~2 L 水樣過濾酸化, 樣品解凍后立即通過0.7 μm 玻璃纖維濾膜(預先在馬弗爐內450 ℃下灼燒4 h)過濾水樣, 并用鹽酸酸化至pH = 2。

    固相萃取 使用安捷倫Bond Elut PPL (1 g)固相萃取小柱通過固相萃取(SPE)從水樣中分離溶解性有機物(dissolve organic matter, DOM) (以下將固相可提取的DOM 稱為SPE-DOM), 使用固相萃取儀批量處理樣品并控制流速, 在使用前立即用甲醇(10 mL, HPLC 級)沖洗柱子, 然后將經過濾酸化的樣品以小于15 mL/min 的流速通過萃取小柱, 在萃取之后用20 mL 0.01 mol/L HCl 脫鹽, 將小柱在氮氣流或空氣流下干燥, 并將被填料吸附的SPE-DOM用8 mL 甲醇洗脫, 定容至10 mL。將10 mL 洗脫物儲存在酸洗過的聚乙烯瓶中, 冷藏保存。

    確定提取率 使用氮吹儀將1 mL 甲醇提取物蒸發(fā)至干燥, 將其重新溶解在0.1 倍固相萃取所用水樣量中, 測定該溶液溶解性有機碳濃度, 并計算其占原始樣品溶解性有機碳濃度的比例, 以確定每個樣品溶解性有機碳提取率。本研究中溶解性有機碳提取率均值為60%, 與使用該方法的其他研究[8]一致。

    硝酸氧化 將1~2 mL 甲醇提取物轉移到2 mL玻璃安瓿瓶中, 在氮氣流下蒸發(fā)至干, 并在瓶中加入0.5 mL HNO3(65%), 然后將安瓿瓶使用酒精噴燈密封, 置于100 mL 不銹鋼高壓消解罐中, 在烘箱170 ℃下消解9 h。消解完成后使用氮吹儀50 ℃氮氣流下蒸發(fā)硝酸, 將處理后的樣品用500 μL 流動相A 溶解, 并用0.22 μm 濾膜過濾, 轉移至進樣小瓶等待進樣。

    1.2.3 液相色譜分析

    使用配備有二級管陣列檢測器的HPLC 檢測6種BPCAs 標準物質(苯三羧酸 (1,2,3-B3CA、1,2,4- B3CA、1,3,5-B3CA)、苯四羧酸(1,2,4,5-B4CA)、苯五羧酸(B5CA)、苯六羧酸(B6CA))配制成的混標溶液(5~80 μmol/L), 優(yōu)化色譜條件以獲得最佳分離效果。對比文獻資料的檢測條件[8], 本研究色譜條件的柱溫由16 ℃變?yōu)?0 ℃, 進樣量由4 μL 調至10 μL, 流動相混合梯度中流動相B 的初始值由6%變?yōu)?%, 可獲得較好的分離效果(表1 和表2)。使用梯度混標溶液確定檢測限并制作標準曲線, 然后將預處理好的樣品在相同色譜條件下上機測試。

    1.2.4 質量控制和質量保證

    為保證實驗結果的準確性, 在進行樣品測定的同時進行空白樣和平行樣的測定, 使用BPCAs 混合標準溶液計算了整個分析過程的 BPCAs 回收率,各BPCAs 回收率在85%~93%之間(N=8)。分析使用的標準曲線, 要求至少5 個濃度梯度, 相關系數(shù)R2必須大于0.99, 各BPCAs 標準曲線及檢測限如表3所示。每分析10 個樣品, 需校準標準曲線, 要求分析值與標準值相差不大于5%, 平行樣與樣品的測定結果偏差小于10%。

    表1 色譜條件參數(shù) Table 1 Settings of HPLC parameters

    表2 流動相混合梯度 Table 2 Gradient of mobile phase

    表3 各BPCAs 標準曲線及檢測限(信噪比s/n=3) Table 3 Standard curve and detection limit (s/n=3) of each BPCAs

    1.2.5 溶解態(tài)黑碳濃度和沉降通量計算

    總沉降樣品每月蒸發(fā)量差異較大, 其采樣體積不具有可比性, 因此總沉降樣品體積取對應月份濕沉降體積的值進行計算。不考慮收集桶內二次起塵和液體可能發(fā)生的物理、化學和生物過程, 總沉降、濕沉降各BPCAs 及溶解態(tài)黑碳的濃度、沉降通量計算過程如下。

    首先使用標準曲線函數(shù)將HPLC 分析所得的各BPCAs 積分面積轉換為濃度cBPCAs; 濕沉降消解產生的各 BPCAs 濃度c1和總沉降消解產生的各BPCAs 濃度c2由式(1)、式(2)計算可得。

    式中, 0.5 為進樣前溶解硝酸消解產物所用流動相A 的量, 單位mL;a為從SPE-DOM 中取出氮吹的用量占SPE-DOM 總量比例;b為固相萃取后測得溶解性有機碳濃度除以初始溶解性有機碳濃度所得固相提取率;l為固相萃取用水量。

    式中,V1、V2分別為該月樣品采集的總沉降、濕沉降總體積;A1、A2分別為總沉降面積和濕沉降面積,A1= 0.09 m2,A2= 0.07 m2。

    溶解態(tài)黑碳濃度(以碳計)使用Dittmar[8]根據(jù)溶解態(tài)黑碳分子結構創(chuàng)建的計算公式計算。

    總沉降及濕沉降溶解態(tài)黑碳日沉降通量N計算公式如下, 式中t為每月采樣時間(d)。

    2 結果與討論

    2.1 大氣干濕沉降中溶解態(tài)黑碳含量

    2018 年3 月至2019 年2 月總沉降和濕沉降樣品中溶解態(tài)黑碳及各BPCAs 的濃度、溶解態(tài)黑碳沉降通量結果如表 4 所示。

    表4 總沉降及濕沉降溶解態(tài)黑碳、溶解性有機碳以及各BPCAs 濃度 Table 4 Concentrations of dissolved black carbon, dissolved organic carbon, and each BPCAs in total deposition and wet deposition

    (續(xù)表4)

    總沉降樣品溶解態(tài)黑碳濃度范圍為 0.12~ 0.83 mg/L, 平均值為0.40 mg/L; 濕沉降樣品溶解態(tài)黑碳濃度范圍為0.04~0.18 mg/L, 平均值為0.10 mg/L, 濕沉降占總沉降的平均比例為29.4%, 說明總沉降溶解態(tài)黑碳含量中干沉降的貢獻更大, 約占71%。最新研究[29]顯示, 在美國東南部長葉松林中, 濕沉降中溶解態(tài)黑碳濃度均值為0.024 mg/L。與森林降雨相比, 武漢市降雨中的溶解態(tài)黑碳濃度較高, 這是由于森林地區(qū)黑碳排放源較少, 且樹木對氣溶膠中的黑碳有一定的吸附作用, 黑碳濃度較低, 而武漢市區(qū)黑碳氣溶膠的主要來源是化石燃料燃燒源和以汽車尾氣排放為代表的流動源[30], 黑碳本底濃度較高。我國長江下游溶解態(tài)黑碳濃度均值為0.04 mg/L[32], 與武漢大氣干濕沉降中的溶解態(tài)黑碳濃度相比較低, 且質量平衡計算表明, 長江的大部分溶解態(tài)黑碳(78%~85%)來自生物質燃燒, 只有15%~22%來自化石燃料燃燒[32], 這表明武漢市長江中的溶解態(tài)黑碳大部分來源于土壤, 少量來源于大氣沉降。

    總沉降溶解態(tài)黑碳每月的日均沉降通量范圍為0.26~1.75 mg·m?2·d?1, 年通量為269 mg·m?2·a?1, 濕沉降溶解態(tài)黑碳每月的日均沉降通量范圍為0.08~ 0.29 mg·m?2·d?1, 年通量為65 mg·m?2·a?1, 干沉降(總沉降減濕沉降)溶解態(tài)黑碳每月的日均沉降通量范圍為0.15~1.51 mg·m?2·d?1, 年通量為204 mg·m?2·a?1。根據(jù)溶解態(tài)黑碳沉降通量平均值粗略估計出武漢市每年通過干濕沉降直接進入湖水中的溶解態(tài)黑碳達600 t, 其中通過降雨進入湖水中的溶解態(tài)黑碳約占1/4。我國黃海和東海[19]干沉降中水溶性黑碳通量分別為(0.042±0.033) mg·m?2·d?1和(0.10±0.043) mg·m?2·d?1, 低于武漢市干沉降溶解態(tài)黑碳通量均值(0.56 mg·m?2·d?1), 黃海和東海的主要河流排放分別來自黃河和長江, 水溶性黑碳的大氣沉積量約為黃河和長江流域河流排放量的40%。通過上述對大氣沉降物中溶解態(tài)黑碳含量的定量分析以及與前人的溶解態(tài)黑碳研究對比顯示, 武漢市大氣沉降物中溶解態(tài)黑碳濃度和沉降通量顯著高于森林、海洋等人類活動較少的區(qū)域。大氣沉降物并非環(huán)境水體中溶解態(tài)黑碳最主要的來源, 但其對水體中溶解態(tài)黑碳的含量也有一定的貢獻。

    由于BPCAs 的羧基數(shù)目與苯環(huán)的共邊數(shù)目相對應, 所以單個BPCAs 分子的相對比例可指示溶解態(tài)黑碳分子大小[33]??偝两蹬c濕沉降中溶解態(tài)黑碳衍生的各BPCAs 比例無明顯差異, 總沉降及濕沉降中溶解態(tài)黑碳衍生的各BPCAs 占總BPCAs 的比例分別為B3CA(38.63%)、B4CA(38.06%)、B5CA(14.56%)和B6CA(8.76%), 該比例與已報道的河流溶解態(tài)黑碳研究結果差異較大。潮間帶系統(tǒng)河流[33]中該比例為B3CA(15%)、B4CA (43%)、B5CA (31%)和B6CA (11%), 森林河流[17]該比例為 B3CA(13%)、B4CA (41%)、B5CA (36%)和B6CA (10%), 草原河流[16]該比例為B3CA(8%)、B4CA (38%)、B5CA (43%)和B6CA (11%)。河流中溶解態(tài)黑碳衍生的各BPCA 比例具有B4CA 和B5CA 相對比例較高的統(tǒng)一特征。與河流對比而言, 總沉降中溶解態(tài)黑碳的B5CA 和B6CA 相對比例偏小, 說明總沉降中的溶解態(tài)黑碳分子具有較少的苯環(huán)數(shù)目和相對分子質量??缮藤彽某鞘袎m埃標準物質中溶出的溶解態(tài)黑碳衍生BPCAs 的特征含有大量的B3CA+B4CA (86%)和少量B5CA+B6CA (14%)[34], 與本研究中各BPCAs 的相對比例類似。這表明從人為的大氣沉積顆粒物中提取的溶解態(tài)黑碳分子具有較少的縮聚結構。土壤中木炭溶出的溶解態(tài)黑碳的研究發(fā)現(xiàn), 在土壤中老化了10 年的木炭的水溶性部分中具有與河流相似的單個BPCAs 比例, 而新產生木炭的水溶性部分很少, 且僅含有 B3CA(88%)和少量B4CA(12%)[35]。這表明新產生的木炭中釋放出的可溶性部分較少且分子較小, 在土壤中經歷生物或非生物氧化過程后, 增加了含氧官能團的數(shù)量和總體極性, 使其可釋放出更多較大分子的可溶性組分。

    本研究與河流、土壤及城市塵埃研究中溶解態(tài)黑碳衍生的各BPCAs 相對比例的對比可發(fā)現(xiàn), 河流中溶解態(tài)黑碳衍生的各BPCAs 相對比例與土壤中老化多年的黑碳溶出的溶解態(tài)黑碳具有相似特征(B4CA 和B5CA 相對比例較高, 總比例達(77±3)%), 可以推測河流溶解態(tài)黑碳主要來源于土壤中長時間老化的黑碳。而大氣沉降物、城市塵埃和新產生木炭溶出的溶解態(tài)黑碳衍生的 BPCAs 中 B3CA 和B4CA 的相對比例較高, 這可能是未經長時間老化的黑碳溶出溶解態(tài)黑碳的統(tǒng)一特征。但 B3CA+ B4CA 總比例大氣沉降物(77%)<城市塵埃(86%)<新產生木炭(約 100%), 這些樣品中黑碳的燃燒源材料、生成條件、暴露環(huán)境和氧化時間都有一定的差異, 從而導致溶解態(tài)黑碳衍生的各BPCAs 相對比例的差異。

    2.2 溶解態(tài)黑碳與溶解性有機碳相關性分析

    溶解態(tài)黑碳是溶解性有機碳中具有濃縮稠環(huán)結構特征的一部分, 前人研究顯示, 全球河流中溶解態(tài)黑碳約占溶解性有機碳的10%[8,29], 海洋中溶解態(tài)黑碳約占溶解性有機碳的2%[13,19], 海洋氣溶膠中溶解態(tài)黑碳約占溶解性有機碳的2.7%[19]。美國東南部森林降雨中溶解態(tài)黑碳占溶解性有機碳的1.6%[29]。本研究中總沉降溶解態(tài)黑碳占溶解性有機碳的平均比例為9.26%, 濕沉降的比例為6.51%, 該比例均高于海洋氣溶膠和森林降雨。本研究中采樣點位于武漢市市區(qū), 與海洋和森林大氣中黑碳的來源不同, 武漢市大氣中黑碳的來源主要為人為源, 如化石燃料燃燒源和機動車尾氣等[30]。本研究中采樣頻率為每月1 次, 樣品在放置過程中可能會出現(xiàn)微生物降解或光降解反應, 對溶解態(tài)黑碳占溶解性有機碳的比例造成一定影響。對干濕沉降DBC/DOC的比值與氣象條件進行相關性分析的結果顯示, 干沉降DBC/DOC 比值與風速在0.05 顯著性水平上顯著負相關(相關系數(shù)R2= ?0.69, 顯著性), 濕沉降DBC/DOC 比值與溫度在0.05 顯著性水平上顯著負相關(R2= ?0.69), 說明武漢市干沉降和濕沉降中DBC/DOC 比值也受氣象條件影響。

    前人研究中不少關于河流溶解態(tài)黑碳的研究顯示, 溶解態(tài)黑碳濃度與溶解性有機碳濃度具有良好的線性相關關系[9-10], 對本研究的總沉降和濕沉降樣品進行了溶解態(tài)黑碳和溶解性有機碳的相關性分析(圖2), 結果顯示, 總沉降和濕沉降溶解態(tài)黑碳濃度與溶解性有機碳濃度均有十分顯著的正相關關系, 總沉降溶解態(tài)黑碳與溶解性有機碳相關系數(shù)R2= 0.94, 擬合函數(shù)見式(5)。式(5)與樣本量較大的世界河流溶解態(tài)黑碳研究[9]中的擬合函數(shù)十分相近, 見式(6)。濕沉降溶解態(tài)黑碳濃度與溶解性有機碳濃度相關系數(shù)R2= 0.90, 擬合函數(shù)見式(7)。

    由于溶解性有機碳為常見的環(huán)境和水文監(jiān)測指標, 溶解性有機碳與溶解態(tài)黑碳良好的線性關系有助于利用樣品溶解性有機碳含量指示溶解態(tài)黑碳的大致濃度。

    2.3 干濕沉降溶解態(tài)黑碳季節(jié)分布特征

    圖2 總沉降、濕沉降溶解態(tài)黑碳與溶解性有機碳散點分布圖 Fig.2 Scatter distributions of dissolved black carbon and dissolved organic carbon in total deposition, wet deposition

    圖3 總沉降及濕沉降溶解態(tài)黑碳沉降通量、降雨量 和風向季節(jié)變化 Fig.3 Seasonal changes in dissolved black carbon deposition flux, rainfall, and wind direction in total deposition and wet deposition

    為去除降雨對干濕沉降濃度的影響, 以沉降通量表示干濕沉降溶解態(tài)黑碳的季節(jié)變化。武漢市屬北亞熱帶季風性(濕潤)氣候, 具有常年雨量豐沛、熱量充足和四季分明等特點, 溶解態(tài)黑碳沉降通量隨季節(jié)變化而呈現(xiàn)不同的分布特征(圖3)??偝两等芙鈶B(tài)黑碳沉降通量季節(jié)變化為春季(1.01 mg·m?2·d?1)>冬季(0.84 mg·m?2·d?1)>夏季(0.79 mg·m?2·d?1)>秋季(0.31 mg·m?2·d?1), 濕沉降溶解態(tài)黑碳沉降通量季節(jié)變化為冬季(0.25 mg·m?2·d?1)>春季(0.23 mg·m?2·d?1)> 與O3的非均相化學反應使黑碳表面生成酮、內酯和酸酐等含氧物種, 且光照能夠促進O2與黑碳的非均相化學反應而使黑碳表面氧化[36], 因此夏季大氣中的黑碳氧化程度較高, 這可能是夏季總沉降樣品溶解態(tài)黑碳衍生的BPCAs 中B3CA 相對比例較低的重要原因之一。武漢市黑碳氣溶膠濃度主要來源為人為源[30], 受人類活動影響不同季節(jié)大氣中黑碳的排放源有所不同, 如夏季由于暑假導致采樣點附近機動車排放黑碳減少(本研究中采樣點位于校園內)。除此之外, 使用后向軌跡法模擬氣團軌跡發(fā)現(xiàn), 夏季氣團軌跡來源主要為南方, 其他季節(jié)氣團軌跡來源 夏季(0.17 mg·m?2·d?1)>秋季(0.11 mg·m?2·d?1)。秋季降雨量最少, 總沉降和濕沉降溶解態(tài)黑碳沉降通量均為最少。總沉降溶解態(tài)黑碳中各BPCAs 相對比例隨季節(jié)變化也出現(xiàn)不同特征(如圖4), 其中夏季變化幅度最大, 夏季總沉降樣品中 B4CA、B5CA 和B6CA 相對比例都比其他季節(jié)高, 而B3CA 相對比例最低。這說明了夏季總沉降樣品中的溶解態(tài)黑碳具有較少的縮聚結構和相對分子質量。夏季與其他季節(jié)相比最顯著的特征是溫度高、太陽輻射強度大, 大氣顆粒物中, 可吸入顆粒物(PM10)、細顆粒物(PM2.5)濃度低, O3濃度較高(可參考圖5)。大氣氧化性氣體與黑碳的非均相化學反應研究中發(fā)現(xiàn), 黑碳主要為北方。大氣中黑碳的來源和傳輸路徑的差異可能也是夏季總沉降樣品溶解態(tài)黑碳衍生的BPCAs中B3CA 相對比例較低的原因之一。

    圖4 總沉降各BPCAs 相對比例季節(jié)變化 Fig.4 Seasonal changes in BPCAs in the total deposition

    2.4 干濕沉降溶解態(tài)黑碳季節(jié)分布影響因素

    使用 SPSS 軟件對溶解態(tài)黑碳沉降通量、各BPCAs 相對比例與降雨量、溫度和風速等氣象因素進行了相關性分析, 結果顯示, 濕沉降溶解態(tài)黑碳沉降通量與降雨量在0.01 顯著性水平上呈顯著正相關關系(R2= 0.78), 說明濕沉降溶解態(tài)黑碳主要受降雨量影響。干沉降溶解態(tài)黑碳(總沉降減濕沉降)與降雨量的相關性很弱(R2= 0.16), 而與風速在0.01 顯著性水平上呈顯著正相關關系(R2= 0.74), 說明風速可能是影響干沉降溶解態(tài)黑碳的重要因素。總沉降溶 解態(tài)黑碳與降雨量相關性較弱(R2= 0.52), 這是因為降雨量較大, 降雨時間較長時, 濕沉降溶解態(tài)黑碳較多, 但由于空氣相對清潔, 干沉降中溶解態(tài)黑碳比例將會減少, 如圖3 中5 月降雨量高于4 月, 但5月干沉降溶解態(tài)黑碳較少, 所以5 月總沉降溶解態(tài)黑碳低于4 月。除此之外總沉降中溶解態(tài)黑碳通量還受風速、風向和空氣質量等多種因素共同影響, 如圖3 中3 月和4 月降雨量和風速相差不大, 但4月總沉降溶解態(tài)黑碳通量遠高于3 月, 3 月主風向為東北風, 而4 月為南風, 使用后向軌跡法模擬氣團軌跡(圖6)發(fā)現(xiàn), 3 月氣團軌跡為東北-海陸混合型(NE-mixed), 4 月氣團軌跡為南部-大陸型(S-land), 推測是由于污染物來源及運動軌跡不同導致了溶解態(tài)黑碳沉降通量的差異。

    圖5 溶解態(tài)黑碳沉降通量與各氣象因素(a)、空氣質量指標(b)的季節(jié)變化 Fig.5 Seasonal changes in deposition fluxes of dissolved black carbon, meteorological factors (a), and air quality indicators (b)

    圖6 采樣點3 月(a)及4 月(b)典型氣團軌跡圖 Fig.6 Typical air mass trajectory at sampling points in March (a) and April (b)

    對各BPCAs 相對比例和氣象因素的相關性分析發(fā)現(xiàn), B3CA 相對比例與溫度(R2= ?0.61)、能見度(R2= ?0.62)和日照時數(shù)(R2= ?0.58)等氣象條件在0.05 顯著性水平上呈顯著負相關關系。這證明了總沉降溶解態(tài)黑碳中B3CA 相對比例隨太陽輻射強度增大而減小, 而太陽輻射強度與大氣黑碳老化關系密切[32], 這也側面證明了溶解態(tài)黑碳衍生的B3CA相對比例與黑碳氧化程度具有一定的相關關系。

    圖5 為溶解態(tài)黑碳總沉降通量、濕沉降通量及各氣象因素與空氣質量指標季節(jié)變化趨勢。氣象因素中降雨量、濕度、風速、總云量與總沉降、濕沉降溶解態(tài)黑碳通量的變化趨勢基本一致, 而溫度、日照時數(shù)和總沉降、濕沉降溶解態(tài)黑碳通量的變化趨勢相反, 溫度高、日照長的季節(jié)總沉降、濕沉降溶解態(tài)黑碳通量較少。之前武漢市黑碳氣溶膠的研究中顯示, 黑碳氣溶膠濃度季節(jié)變化為冬季>秋季>春季>夏季, 且黑碳濃度與PM2.5、CO 濃度有顯著的正相關關系, 與O3濃度有顯著負相關[30-31], 據(jù)此相關關系推測空氣質量指標中PM2.5、CO 和O3等指標可能通過影響黑碳氣溶膠濃度而間接對大氣沉降中溶解態(tài)黑碳通量造成一定影響??諝赓|量指標中各指標與總沉降、濕沉降溶解態(tài)黑碳通量的變化趨勢都不太一致, 這是因為秋季降雨量最低, 濕沉降對大氣污染物的去除率較低, 因此秋季PM10、PM2.5、CO、SO2和 NO2空氣質量分指數(shù)(individual air quality index, IAQI)相對較高, 而總沉降、濕沉降溶解態(tài)黑碳主要受降雨量影響沉降通量較低, 除秋季外其他季節(jié)PM10、PM2.5、CO 與總沉降、濕沉降溶解態(tài)黑碳通量的變化趨勢基本一致。O3主要受太陽輻射影響, 與總沉降溶解態(tài)黑碳趨勢相反。以上分析說明, 降雨量、風速和溫度等氣象條件是武漢市大氣沉降物中溶解態(tài)黑碳沉降通量季節(jié)變化的主要影響因素, 但氣象條件變化較小時, 空氣質量(PM2.5、CO 和O3)對大氣沉降溶解態(tài)黑碳沉降通量的季節(jié)變化也有重要影響。

    3 結 論

    (1) 武漢市洪山區(qū)中國地質大學空氣質量監(jiān)測站2018 年3 月至2019 年2 月總沉降溶解態(tài)黑碳月濃度范圍為0.12~0.83 mg/L, 平均值為0.40 mg/L, 與世界河流溶解態(tài)黑碳均值十分接近; 濕沉降溶解態(tài)黑碳濃度范圍為0.04~0.18 mg/L, 平均值為0.10 mg/L, 濕沉降占總沉降平均比例為29.4%。

    (2) 武漢市總沉降及濕沉降溶解態(tài)黑碳年沉降通量分別為269 mg·m?2·a?1和 65 mg·m?2·a?1。溶解態(tài)黑碳沉降通量隨季節(jié)變化而呈現(xiàn)不同的分布特征: 總沉降溶解態(tài)黑碳沉降通量季節(jié)變化為春季>冬季>夏季>秋季, 濕沉降溶解態(tài)黑碳沉降通量季節(jié)變化為冬季>春季>夏季>秋季。

    (3) 相關性分析表明, 濕沉降溶解態(tài)黑碳沉降通量的季節(jié)變化主要受降雨量影響, 干沉降溶解態(tài)黑碳沉降通量的季節(jié)變化主要受風速影響, 而總沉降溶解態(tài)黑碳沉降通量的季節(jié)變化除受降雨量影響外還受風速、風向和空氣質量等因素共同影響。

    (4) 總沉降中溶解態(tài)黑碳占溶解性有機碳的9.26%, 濕沉降中溶解態(tài)黑碳占溶解性有機碳的6.51%, 且總沉降和濕沉降中溶解態(tài)黑碳與溶解性有機碳均有顯著的正相關關系。總沉降溶解態(tài)黑碳中各BPCAs 相對比例隨季節(jié)變化出現(xiàn)不同特征, 相關性分析發(fā)現(xiàn), B3CA 相對比例與溫度、能見度和日照時數(shù)等, 以及太陽輻射強度相關的氣象條件均呈顯著負相關關系。

    兩位審稿專家和編輯部對本次研究提出了非常寶貴和細致的修改意見, 筆者在此表示衷心地感謝!

    猜你喜歡
    大氣
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    首次發(fā)現(xiàn)系外行星大氣中存在CO2
    科學(2022年5期)2022-12-29 09:48:56
    宏偉大氣,氣勢與細膩兼?zhèn)?Vivid Audio Giya G3 S2
    太赫茲大氣臨邊探測儀遙感中高層大氣風仿真
    有“心氣”才大氣
    如何“看清”大氣中的二氧化碳
    學生天地(2020年18期)2020-08-25 09:29:24
    大氣穩(wěn)健的美式之風Polk Audio Signature系列
    稚拙率真 圓融大氣
    中國篆刻(2017年3期)2017-05-17 06:20:46
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    新農業(yè)(2016年18期)2016-08-16 03:28:27
    母亲3免费完整高清在线观看| 国产亚洲精品久久久久久毛片| 日本五十路高清| 国产野战对白在线观看| 久久久久久久午夜电影| 黄片大片在线免费观看| 色尼玛亚洲综合影院| 天天一区二区日本电影三级| 午夜福利高清视频| 国产黄a三级三级三级人| 欧美黑人巨大hd| www国产在线视频色| 午夜日韩欧美国产| 亚洲精品亚洲一区二区| 欧美一区二区国产精品久久精品| 精品无人区乱码1区二区| 成人午夜高清在线视频| 午夜老司机福利剧场| 欧美另类亚洲清纯唯美| 国语自产精品视频在线第100页| 搡老岳熟女国产| 亚洲熟妇中文字幕五十中出| 美女高潮喷水抽搐中文字幕| 91久久精品国产一区二区成人 | 成人国产一区最新在线观看| 色在线成人网| 热99在线观看视频| 色综合站精品国产| 高清在线国产一区| 男女做爰动态图高潮gif福利片| 97超视频在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 国内揄拍国产精品人妻在线| 色综合亚洲欧美另类图片| 国产高清视频在线观看网站| 露出奶头的视频| 亚洲aⅴ乱码一区二区在线播放| 免费观看的影片在线观看| 操出白浆在线播放| 亚洲av不卡在线观看| 国产又黄又爽又无遮挡在线| 亚洲在线观看片| 日日夜夜操网爽| 亚洲国产精品久久男人天堂| 久久伊人香网站| 一卡2卡三卡四卡精品乱码亚洲| xxx96com| 99久久99久久久精品蜜桃| 一级毛片女人18水好多| 2021天堂中文幕一二区在线观| 国产精品亚洲一级av第二区| 欧美高清成人免费视频www| 丰满人妻熟妇乱又伦精品不卡| 午夜久久久久精精品| 国产精品野战在线观看| 熟妇人妻久久中文字幕3abv| 久久香蕉精品热| 99久久久亚洲精品蜜臀av| 亚洲成人久久性| 五月伊人婷婷丁香| 精品乱码久久久久久99久播| 中文字幕熟女人妻在线| 欧美成人a在线观看| 日韩亚洲欧美综合| 国产精品 欧美亚洲| 亚洲最大成人中文| 久久精品国产亚洲av香蕉五月| 国产一区二区在线av高清观看| 国产一区二区在线观看日韩 | 看片在线看免费视频| 成人特级黄色片久久久久久久| 国产精品99久久99久久久不卡| 久久国产乱子伦精品免费另类| 国产成人福利小说| 99久久无色码亚洲精品果冻| 日本 欧美在线| 国产精品久久久久久久久免 | 国产亚洲欧美98| 亚洲在线观看片| 黄色成人免费大全| 亚洲精品日韩av片在线观看 | 午夜激情欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 9191精品国产免费久久| 国产一区二区三区视频了| 精品人妻一区二区三区麻豆 | 国模一区二区三区四区视频| 亚洲欧美日韩高清专用| 国产成人啪精品午夜网站| 亚洲av成人精品一区久久| 天堂av国产一区二区熟女人妻| 特级一级黄色大片| 熟女电影av网| 国产男靠女视频免费网站| 久久久久久久久大av| 免费看美女性在线毛片视频| 午夜免费激情av| 国产亚洲欧美98| xxx96com| www日本在线高清视频| 国产免费一级a男人的天堂| 日韩成人在线观看一区二区三区| 国产91精品成人一区二区三区| 亚洲av一区综合| 欧美一级毛片孕妇| 国产一级毛片七仙女欲春2| 变态另类丝袜制服| 亚洲精品456在线播放app | 9191精品国产免费久久| 婷婷六月久久综合丁香| 国产淫片久久久久久久久 | 国语自产精品视频在线第100页| 国产精品一及| 在线看三级毛片| 亚洲va日本ⅴa欧美va伊人久久| av欧美777| 18美女黄网站色大片免费观看| 黑人欧美特级aaaaaa片| 一区二区三区高清视频在线| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 三级毛片av免费| 国产免费男女视频| 深夜精品福利| 久久精品国产99精品国产亚洲性色| 色老头精品视频在线观看| 国产精品久久久久久精品电影| 日本一二三区视频观看| 日日夜夜操网爽| 亚洲精品一区av在线观看| а√天堂www在线а√下载| 色综合欧美亚洲国产小说| 精品不卡国产一区二区三区| 九色成人免费人妻av| 免费在线观看成人毛片| 亚洲va日本ⅴa欧美va伊人久久| 黄色女人牲交| 久久久色成人| 日韩国内少妇激情av| 久久99热这里只有精品18| 免费观看的影片在线观看| 一区二区三区国产精品乱码| 啦啦啦免费观看视频1| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 午夜免费观看网址| 美女高潮的动态| 国产黄片美女视频| 综合色av麻豆| 中文字幕人妻熟人妻熟丝袜美 | 亚洲精品456在线播放app | 搞女人的毛片| 免费高清视频大片| 色综合欧美亚洲国产小说| 亚洲国产高清在线一区二区三| 国产一区二区亚洲精品在线观看| 亚洲久久久久久中文字幕| 亚洲精品在线美女| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久久久久久久| 亚洲成人久久性| 免费大片18禁| 中文字幕久久专区| 午夜a级毛片| 好男人在线观看高清免费视频| 欧美乱码精品一区二区三区| 99久久精品热视频| 夜夜夜夜夜久久久久| 手机成人av网站| 男女床上黄色一级片免费看| 欧美另类亚洲清纯唯美| 精品人妻偷拍中文字幕| 性色avwww在线观看| a级一级毛片免费在线观看| 精品久久久久久,| 日韩欧美精品免费久久 | 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| 真实男女啪啪啪动态图| 国产午夜精品久久久久久一区二区三区 | 最近视频中文字幕2019在线8| 国产精品久久电影中文字幕| 国产爱豆传媒在线观看| www.www免费av| 中文字幕高清在线视频| 高清日韩中文字幕在线| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 乱人视频在线观看| 天堂网av新在线| 亚洲自拍偷在线| 高清毛片免费观看视频网站| 国产精品99久久久久久久久| 变态另类丝袜制服| 亚洲av电影不卡..在线观看| 一区二区三区高清视频在线| 岛国在线免费视频观看| 日本三级黄在线观看| 国产激情欧美一区二区| 午夜精品在线福利| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| h日本视频在线播放| 色噜噜av男人的天堂激情| 精品久久久久久久毛片微露脸| 久久久色成人| 日本三级黄在线观看| 国产伦精品一区二区三区四那| 国产精品亚洲一级av第二区| 最新中文字幕久久久久| 超碰av人人做人人爽久久 | 国产高清有码在线观看视频| aaaaa片日本免费| 18+在线观看网站| 亚洲av成人不卡在线观看播放网| 免费高清视频大片| 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 黄色成人免费大全| 一进一出好大好爽视频| 亚洲欧美精品综合久久99| av在线蜜桃| 亚洲欧美日韩卡通动漫| 国产欧美日韩精品一区二区| 免费电影在线观看免费观看| 精品日产1卡2卡| 99久久久亚洲精品蜜臀av| 久久精品综合一区二区三区| 色精品久久人妻99蜜桃| 国内精品一区二区在线观看| 亚洲真实伦在线观看| 九色成人免费人妻av| 亚洲成a人片在线一区二区| 欧美乱色亚洲激情| 久久亚洲精品不卡| 最近最新免费中文字幕在线| 国产色爽女视频免费观看| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 久久久色成人| a级一级毛片免费在线观看| 激情在线观看视频在线高清| 婷婷精品国产亚洲av在线| 天天添夜夜摸| 最新在线观看一区二区三区| 少妇高潮的动态图| 亚洲av成人不卡在线观看播放网| 岛国在线免费视频观看| 午夜精品久久久久久毛片777| a级一级毛片免费在线观看| 欧美+亚洲+日韩+国产| 久久香蕉国产精品| 一区二区三区免费毛片| 成人18禁在线播放| 亚洲国产精品999在线| 99精品久久久久人妻精品| 一个人看的www免费观看视频| 一个人免费在线观看的高清视频| 窝窝影院91人妻| 丰满的人妻完整版| 国产成人啪精品午夜网站| av中文乱码字幕在线| 久久九九热精品免费| 手机成人av网站| 欧美性猛交╳xxx乱大交人| 制服人妻中文乱码| 热99在线观看视频| 亚洲avbb在线观看| 亚洲国产精品sss在线观看| 99在线视频只有这里精品首页| 黄色成人免费大全| 一个人观看的视频www高清免费观看| 日韩欧美在线乱码| 丰满乱子伦码专区| 国产欧美日韩精品一区二区| 天堂√8在线中文| 国内久久婷婷六月综合欲色啪| 日本与韩国留学比较| 一本一本综合久久| 精品国产超薄肉色丝袜足j| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 日日干狠狠操夜夜爽| 99久久99久久久精品蜜桃| 黄色女人牲交| 精品熟女少妇八av免费久了| 国产亚洲欧美在线一区二区| 亚洲成人久久爱视频| 久久精品国产亚洲av涩爱 | 成年免费大片在线观看| 欧美乱码精品一区二区三区| eeuss影院久久| 国产精品美女特级片免费视频播放器| 国内精品美女久久久久久| 亚洲五月天丁香| 日本熟妇午夜| 日本a在线网址| 在线国产一区二区在线| 成人特级av手机在线观看| 久久久久久久亚洲中文字幕 | 成人无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 18禁在线播放成人免费| 亚洲无线观看免费| 国产精品自产拍在线观看55亚洲| 小蜜桃在线观看免费完整版高清| 麻豆国产97在线/欧美| xxx96com| 男女床上黄色一级片免费看| 国产又黄又爽又无遮挡在线| 热99在线观看视频| 久久伊人香网站| av国产免费在线观看| 免费在线观看成人毛片| 女人高潮潮喷娇喘18禁视频| 99久久精品国产亚洲精品| 国内少妇人妻偷人精品xxx网站| 男人的好看免费观看在线视频| 黄色日韩在线| 成人无遮挡网站| 国产乱人伦免费视频| 69人妻影院| 搡老熟女国产l中国老女人| 成人特级av手机在线观看| 国产探花在线观看一区二区| 国产黄a三级三级三级人| 日本一本二区三区精品| 国产精品一区二区三区四区久久| 日韩大尺度精品在线看网址| 免费高清视频大片| 窝窝影院91人妻| 国产精品一区二区三区四区久久| 真实男女啪啪啪动态图| 久久久久久久久大av| svipshipincom国产片| 国产一区二区在线观看日韩 | 日韩欧美国产一区二区入口| 免费在线观看成人毛片| 国产野战对白在线观看| 狂野欧美激情性xxxx| 午夜影院日韩av| 免费在线观看亚洲国产| 亚洲乱码一区二区免费版| 在线免费观看不下载黄p国产 | 国产99白浆流出| a在线观看视频网站| 婷婷精品国产亚洲av| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 国产伦人伦偷精品视频| 免费看光身美女| 男女之事视频高清在线观看| 国产探花极品一区二区| 叶爱在线成人免费视频播放| 久久精品亚洲精品国产色婷小说| а√天堂www在线а√下载| 两性午夜刺激爽爽歪歪视频在线观看| 成人午夜高清在线视频| 国产精品一区二区三区四区免费观看 | 两个人看的免费小视频| 一二三四社区在线视频社区8| 午夜精品一区二区三区免费看| 国产黄a三级三级三级人| e午夜精品久久久久久久| 老司机在亚洲福利影院| 亚洲午夜理论影院| 九九久久精品国产亚洲av麻豆| 午夜激情欧美在线| 精品一区二区三区视频在线观看免费| 国产伦人伦偷精品视频| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 欧美丝袜亚洲另类 | 亚洲av免费在线观看| 亚洲av电影不卡..在线观看| 午夜免费成人在线视频| 成年人黄色毛片网站| 三级男女做爰猛烈吃奶摸视频| 岛国在线免费视频观看| 日韩大尺度精品在线看网址| 精品国产三级普通话版| 欧美3d第一页| 亚洲精品亚洲一区二区| 国模一区二区三区四区视频| 日韩精品中文字幕看吧| 久久久精品欧美日韩精品| 精品一区二区三区av网在线观看| 亚洲乱码一区二区免费版| 国产精品嫩草影院av在线观看 | 国产精品久久久人人做人人爽| 51国产日韩欧美| 在线天堂最新版资源| 日韩欧美免费精品| 天堂av国产一区二区熟女人妻| avwww免费| 亚洲成人久久爱视频| h日本视频在线播放| 日韩欧美免费精品| 国产精品国产高清国产av| 亚洲aⅴ乱码一区二区在线播放| 婷婷亚洲欧美| 嫁个100分男人电影在线观看| 动漫黄色视频在线观看| 国产美女午夜福利| 精品久久久久久,| 色综合欧美亚洲国产小说| 久久精品夜夜夜夜夜久久蜜豆| 丰满人妻熟妇乱又伦精品不卡| 91麻豆精品激情在线观看国产| 国产色爽女视频免费观看| 日韩欧美三级三区| 成人亚洲精品av一区二区| 麻豆成人av在线观看| 国产成+人综合+亚洲专区| 国产成年人精品一区二区| 久久久国产成人免费| 此物有八面人人有两片| 国产成年人精品一区二区| 免费人成在线观看视频色| 久久香蕉国产精品| 国产亚洲精品一区二区www| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 久久精品91无色码中文字幕| 不卡一级毛片| 美女 人体艺术 gogo| av天堂在线播放| xxxwww97欧美| 亚洲自拍偷在线| 午夜福利成人在线免费观看| 美女黄网站色视频| 看片在线看免费视频| 国产精品爽爽va在线观看网站| 91九色精品人成在线观看| 一夜夜www| 精品福利观看| 亚洲成av人片免费观看| 好男人电影高清在线观看| 免费在线观看亚洲国产| 我的老师免费观看完整版| 亚洲无线在线观看| 在线观看美女被高潮喷水网站 | 国产一区二区激情短视频| 日本撒尿小便嘘嘘汇集6| 亚洲av五月六月丁香网| 久久精品综合一区二区三区| 两个人看的免费小视频| 少妇丰满av| 在线看三级毛片| 国产免费一级a男人的天堂| 亚洲精华国产精华精| 亚洲欧美日韩东京热| 亚洲狠狠婷婷综合久久图片| 成人一区二区视频在线观看| 国产91精品成人一区二区三区| 亚洲男人的天堂狠狠| 亚洲国产精品sss在线观看| 国内精品久久久久久久电影| 成人av一区二区三区在线看| 美女cb高潮喷水在线观看| 在线观看66精品国产| 欧美激情在线99| 欧美+日韩+精品| 叶爱在线成人免费视频播放| 看黄色毛片网站| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清在线视频| 久久中文看片网| 国产av一区在线观看免费| 国产精品香港三级国产av潘金莲| 亚洲精品456在线播放app | 免费观看精品视频网站| 国产91精品成人一区二区三区| 精品久久久久久久久久久久久| 热99在线观看视频| 欧美国产日韩亚洲一区| 国产精品久久久久久精品电影| 亚洲国产精品久久男人天堂| 国产真实伦视频高清在线观看 | 黄色片一级片一级黄色片| 精品久久久久久,| 精品久久久久久久久久久久久| 麻豆国产97在线/欧美| 搞女人的毛片| 青草久久国产| 日本五十路高清| 亚洲国产中文字幕在线视频| 亚洲内射少妇av| 又黄又粗又硬又大视频| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av在线| 日韩欧美在线二视频| 一二三四社区在线视频社区8| 国产午夜精品论理片| 一区二区三区免费毛片| 欧美+亚洲+日韩+国产| 欧美极品一区二区三区四区| 亚洲av成人av| 特级一级黄色大片| 毛片女人毛片| 男女午夜视频在线观看| 欧美大码av| 嫩草影院精品99| 在线观看66精品国产| 18禁黄网站禁片午夜丰满| 日韩 欧美 亚洲 中文字幕| 亚洲人成网站高清观看| 欧美性感艳星| 日韩中文字幕欧美一区二区| 国产视频内射| 高清日韩中文字幕在线| 国产熟女xx| 亚洲精品久久国产高清桃花| 日韩欧美国产一区二区入口| 国内久久婷婷六月综合欲色啪| 有码 亚洲区| 国内精品一区二区在线观看| 一级作爱视频免费观看| 一a级毛片在线观看| 亚洲人成网站高清观看| 久久精品国产自在天天线| 国产成人a区在线观看| 18禁在线播放成人免费| 国产视频一区二区在线看| 精品无人区乱码1区二区| 婷婷六月久久综合丁香| 午夜免费男女啪啪视频观看 | 又粗又爽又猛毛片免费看| 亚洲avbb在线观看| 国产一区二区亚洲精品在线观看| 精品福利观看| 中文字幕精品亚洲无线码一区| 99久久精品热视频| 在线国产一区二区在线| 性色avwww在线观看| 最好的美女福利视频网| 国产色爽女视频免费观看| 香蕉丝袜av| 亚洲电影在线观看av| 欧美丝袜亚洲另类 | 一a级毛片在线观看| 国产精品电影一区二区三区| 免费人成视频x8x8入口观看| 好男人电影高清在线观看| 国产精品 欧美亚洲| 90打野战视频偷拍视频| 一级黄色大片毛片| av片东京热男人的天堂| 亚洲精品色激情综合| 日韩国内少妇激情av| 99在线视频只有这里精品首页| eeuss影院久久| 国产伦精品一区二区三区四那| 母亲3免费完整高清在线观看| 日本免费一区二区三区高清不卡| 亚洲人成网站在线播| 国产三级黄色录像| 久久国产精品影院| 桃色一区二区三区在线观看| 午夜免费观看网址| 国产免费男女视频| 99精品久久久久人妻精品| 午夜a级毛片| 午夜免费男女啪啪视频观看 | 2021天堂中文幕一二区在线观| 久久人妻av系列| 午夜亚洲福利在线播放| 悠悠久久av| 亚洲久久久久久中文字幕| 成人特级av手机在线观看| 亚洲人成网站在线播放欧美日韩| 夜夜夜夜夜久久久久| 国产成年人精品一区二区| 又黄又爽又免费观看的视频| 美女 人体艺术 gogo| 精品乱码久久久久久99久播| 免费看光身美女| 亚洲人与动物交配视频| 宅男免费午夜| www日本在线高清视频| 国产成人av教育| 宅男免费午夜| www日本在线高清视频| bbb黄色大片| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩福利视频一区二区| 日韩欧美一区二区三区在线观看| 欧美+日韩+精品| 欧美日韩一级在线毛片| 黄色女人牲交| 国产精品99久久久久久久久| 欧美日韩福利视频一区二区| 久久精品国产清高在天天线| 国产免费男女视频| 国产美女午夜福利| 美女黄网站色视频| 国产毛片a区久久久久| 午夜激情福利司机影院| 久久中文看片网| 国产精品国产高清国产av| 亚洲欧美激情综合另类| 波多野结衣高清作品| 黄色成人免费大全| 久久精品国产99精品国产亚洲性色| 亚洲精品粉嫩美女一区| 在线观看免费视频日本深夜| 国产精品久久久久久久久免 | 色精品久久人妻99蜜桃| 精品久久久久久,| 看片在线看免费视频| 青草久久国产| 国产高清激情床上av|