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

    全球PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空格局

    2021-12-01 09:50:18張亮林潘竟虎
    中國(guó)環(huán)境科學(xué) 2021年11期
    關(guān)鍵詞:人口區(qū)域濃度

    張亮林,潘竟虎

    全球PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空格局

    張亮林,潘竟虎*

    (西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,甘肅 蘭州 730070)

    基于PM2.5遙感數(shù)據(jù)和人口格網(wǎng)數(shù)據(jù),利用污染物人口暴露風(fēng)險(xiǎn)模型、Theil-Sen Media和Mann-Kendall等方法,分析了2000~2016年全球PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空分布特征,并識(shí)別出暴露高風(fēng)險(xiǎn)區(qū)域.結(jié)果表明,PM2.5遙感數(shù)據(jù)和人口格網(wǎng)數(shù)據(jù)可以客觀地評(píng)價(jià)暴露風(fēng)險(xiǎn)程度.全球PM2.5平均濃度在各大洲差異顯著,PM2.5污染的高值區(qū)域主要分布在東亞、南亞和東南亞.PM2.5質(zhì)量濃度的多年平均值從高到低分別是亞洲14.7μg/m3、非洲8.1μg/m3、歐洲8.03μg/m3、南美洲5.69μg/m3、北美洲4.41μg/m3和大洋洲1.27μg/m3.2000~2016年,全球PM2.5人口暴露風(fēng)險(xiǎn)在宏觀尺度上呈逐漸減少的趨勢(shì),而在區(qū)域內(nèi)則呈現(xiàn)出差異性.空間上,全球PM2.5人口暴露風(fēng)險(xiǎn)各大洲從高到低依次為亞洲5.94、非洲0.62、歐洲0.45、南美洲0.32、北美洲0.27和大洋洲0.01.時(shí)間上,2000~2016年,亞洲和非洲PM2.5人口暴露風(fēng)險(xiǎn)呈增長(zhǎng)趨勢(shì),歐洲和北美洲呈減少趨勢(shì),大洋洲和南美洲變化幅度較小.

    PM2.5;人口暴露風(fēng)險(xiǎn);時(shí)空變化;全球

    空氣污染暴露(Air Pollution Exposure)是指?jìng)€(gè)體居民與空氣污染物直接接觸使其暴露在空氣污染中的狀態(tài)或過(guò)程[1].研究表明,細(xì)顆粒物(PM2.5)是最致命的空氣污染形式,也是造成全球非傳染性疾病的主要原因[2].2017年全球疾病負(fù)擔(dān)(Global Burden of Disease, GBD)研究表明[3],環(huán)境顆粒物(PM)污染導(dǎo)致全球約294萬(wàn)人死亡,其中約35.47%發(fā)生在東南亞、東亞和大洋洲,是造成死亡負(fù)擔(dān)的第八大主要風(fēng)險(xiǎn)因素.PM2.5空氣污染對(duì)健康的影響并沒(méi)有在所有公民之間平等分享,在全球范圍內(nèi),與PM2.5空氣污染相關(guān)的健康風(fēng)險(xiǎn)在中低收入國(guó)家更高,其中91%的過(guò)早死亡是由包括PM2.5在內(nèi)的空氣污染造成的[4].研究發(fā)現(xiàn)[5],在北美、歐洲、亞洲、非洲和大洋洲,較貧窮地區(qū)的PM2.5空氣污染水平往往較高.隨著經(jīng)濟(jì)全球化和區(qū)域一體化快速發(fā)展,全球PM2.5污染監(jiān)測(cè)與治理難度加大.許多國(guó)家受自身經(jīng)濟(jì)狀況制約,尚未建立PM2.5地面監(jiān)測(cè)站點(diǎn),缺乏持續(xù)的PM2.5濃度監(jiān)測(cè)數(shù)據(jù).部分國(guó)家受戰(zhàn)亂影響,政局不穩(wěn),人口普查無(wú)法實(shí)現(xiàn),數(shù)據(jù)的缺失給全球尺度的PM2.5人口暴露研究帶來(lái)了挑戰(zhàn).長(zhǎng)時(shí)間序列的PM2.5人口暴露風(fēng)險(xiǎn)定量化研究,有助于剖析全球PM2.5暴露風(fēng)險(xiǎn)變化,為PM2.5及暴露風(fēng)險(xiǎn)的健康效應(yīng)影響、環(huán)境管理、流行病學(xué)的研究提供可靠、詳盡的輔助基礎(chǔ)數(shù)據(jù)和參考[6],是對(duì)行政單元上暴露風(fēng)險(xiǎn)分析研究的很好補(bǔ)充.因此,基于遙感數(shù)據(jù)定量分析全球PM2.5人口暴露風(fēng)險(xiǎn),有助于區(qū)域聯(lián)動(dòng),進(jìn)行大氣污染治理和暴露風(fēng)險(xiǎn)規(guī)避. 此外,由于PM2.5易于遠(yuǎn)距離乃至跨境傳輸,發(fā)達(dá)國(guó)家借助貿(mào)易使欠發(fā)達(dá)國(guó)家增加了經(jīng)濟(jì)產(chǎn)值,但同時(shí)也可能通過(guò)空氣污染造成不可忽視的經(jīng)濟(jì)和健康損失[7],國(guó)際貿(mào)易中隱含的PM2.5轉(zhuǎn)移問(wèn)題越來(lái)越受到學(xué)者的關(guān)注[8-9],分析全球PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空格局,可為化解全球和區(qū)域生態(tài)環(huán)境風(fēng)險(xiǎn)、推進(jìn)國(guó)際貿(mào)易的健康發(fā)展提供參考.

    已有部分學(xué)者開(kāi)展了PM2.5暴露風(fēng)險(xiǎn)的嘗試性探索研究,但總體上目前PM2.5人口暴露風(fēng)險(xiǎn)研究案例較少,從內(nèi)容來(lái)看主要集中在流行病學(xué)領(lǐng)域[10-12].研究尺度上,大多數(shù)研究集中在某個(gè)特定的國(guó)家(尤其是中國(guó)和印度兩個(gè)人口最多的國(guó)家)[6,13-16],或?qū)⒚總€(gè)國(guó)家作為單一研究目標(biāo)[17].個(gè)別學(xué)者對(duì)全球PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空格局開(kāi)展了研究.Han等研究發(fā)現(xiàn)[18],2010年全球陸地面積的11.0′106km2(約占全球陸地面積的8%)暴露于PM2.5污染(PM2.5質(zhì)量濃度>35μg/m3),自2000年以來(lái)增加了4.3′106km2;2010年,全球19.4億人暴露于PM2.5污染.Roberts等[19]利用哥白尼大氣監(jiān)測(cè)系統(tǒng)(CAMS)調(diào)查了2016~2019年全球人口PM2.5暴露,發(fā)現(xiàn)143個(gè)國(guó)家的平均人口加權(quán)PM2.5總表面濃度超過(guò)10μg/m3的世衛(wèi)組織建議值,每年有6720萬(wàn)人暴露在PM2.5水平分類(lèi)為“危險(xiǎn)”(PM2.5質(zhì)量濃度>250.5μg/m3)的環(huán)境.Li等[20]發(fā)現(xiàn)1998年至2016年全球大陸PM2.5暴露人口占全球總?cè)丝诘?6.0%.傳統(tǒng)的PM2.5人口暴露風(fēng)險(xiǎn)研究多以行政區(qū)為研究單元,雖然在行政單元上可以考慮人口性別、教育程度和收入等社會(huì)經(jīng)濟(jì)因素,但這類(lèi)數(shù)據(jù)在用于估算空氣污染物的實(shí)時(shí)暴露風(fēng)險(xiǎn)時(shí),存在一定的局限性,因?yàn)樗鼈冎皇菍⒚總€(gè)人口普查區(qū)塊的人口視為一個(gè)同質(zhì)實(shí)體,無(wú)法深入研究污染人口暴露的空間異質(zhì)性.遙感因其宏觀、動(dòng)態(tài)、客觀、周期性強(qiáng)等優(yōu)勢(shì),為在不同時(shí)空尺度上開(kāi)展環(huán)境監(jiān)測(cè)和污染分布及人口暴露風(fēng)險(xiǎn)評(píng)估提供了良好的數(shù)據(jù)源[21].基于遙感的暴露風(fēng)險(xiǎn)分析能夠更直觀地識(shí)別區(qū)域內(nèi)部差異,從而為提出差別化、針對(duì)性的應(yīng)對(duì)措施提供了依據(jù)[22].

    基于以上原因,本文利用2000~2016年全球PM2.5與人口分布遙感反演數(shù)據(jù),采用暴露風(fēng)險(xiǎn)模型、趨勢(shì)分析和穩(wěn)定性分析等方法,在像元尺度上研究全球PM2.5人口暴露風(fēng)險(xiǎn)的時(shí)空格局演變特征,并準(zhǔn)確識(shí)別出PM2.5人口高暴露風(fēng)險(xiǎn)區(qū)域,分析PM2.5人口暴露的穩(wěn)定性,探尋暴露風(fēng)險(xiǎn)的時(shí)空演變規(guī)律,以期全面系統(tǒng)地揭示PM2.5人口暴露狀況,為減緩大氣污染、改善人居環(huán)境提供科學(xué)參考.相較于前人研究,本研究以遙感為數(shù)據(jù)源,結(jié)果可落實(shí)到空間柵格上,分辨率高,可視化效果好;采用長(zhǎng)時(shí)間序列衛(wèi)星產(chǎn)品,保證了結(jié)果的連續(xù)性、客觀性和可比較性.

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來(lái)源

    PM2.5質(zhì)量濃度遙感數(shù)據(jù)來(lái)自于MODIS/MISR和SeaWiFS氣溶膠光學(xué)厚度(AOD)聯(lián)合反演得到的全球逐年P(guān)M2.5格網(wǎng)產(chǎn)品(http://earthdata.nasa.gov),時(shí)間序列為2000~2016年,空間分辨率為0.01o×0.01o.人口格網(wǎng)數(shù)據(jù)為L(zhǎng)andScan全球人口格網(wǎng)數(shù)據(jù),獲取自https://landscan.ornl.gov,空間分辨率約1km.LandScan數(shù)據(jù)通過(guò)收集各國(guó)高精度的人口普查數(shù)據(jù)、道路、坡度、土地覆蓋、夜間燈光和城市密度數(shù)據(jù),利用GIS技術(shù)和高分辨率遙感影像模擬人口空間分布,并利用建筑區(qū)、居民點(diǎn)等人口相關(guān)指示因子進(jìn)行模擬結(jié)果的校核和驗(yàn)證,確保了數(shù)據(jù)精度.世界境界線矢量底圖來(lái)自于自然資源部標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站,根據(jù)審圖號(hào)為GS(2016)1665號(hào)的世界地圖繪制,底圖無(wú)修改.

    1.2 研究方法

    1.2.1 PM2.5人口暴露風(fēng)險(xiǎn)計(jì)算 PM2.5人口暴露風(fēng)險(xiǎn)指數(shù)可以評(píng)估空間子單元內(nèi)的暴露狀況[23]:

    式中:是格網(wǎng)號(hào);R表示格網(wǎng)內(nèi)的PM2.5人口暴露風(fēng)險(xiǎn)指數(shù)值;POP是格網(wǎng)內(nèi)的人口數(shù);C是格網(wǎng)內(nèi)的PM2.5濃度值;是研究區(qū)內(nèi)的格網(wǎng)數(shù)之和,這里指全球格網(wǎng)總數(shù)之和.

    1.2.2 趨勢(shì)分析 Theil-Sen Median趨勢(shì)分析可以模擬每個(gè)柵格的變化趨勢(shì),能客觀地反映長(zhǎng)時(shí)間序列PM2.5人口暴露風(fēng)險(xiǎn)的演化趨勢(shì).計(jì)算公式為[24]:

    式中:S是擬合方程的斜率;R是第年的暴露風(fēng)險(xiǎn)指數(shù)值;R是第年的暴露風(fēng)險(xiǎn)指數(shù)值.當(dāng)S>0,表示PM2.5人口暴露風(fēng)險(xiǎn)呈增長(zhǎng)趨勢(shì);反之,表示PM2.5人口暴露風(fēng)險(xiǎn)呈減小趨勢(shì).

    Mann-Kendall檢驗(yàn)是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)的方法,用來(lái)判斷趨勢(shì)的顯著情況.計(jì)算公式為[25]:

    設(shè)定時(shí)間序列{R},=2000, 2001, …, 2016,定義統(tǒng)計(jì)量為:

    其中,

    其中,

    式中:為統(tǒng)計(jì)變量;()是的方差;RR分別表示像元第和年的PM2.5人口暴露風(fēng)險(xiǎn)值;表示時(shí)間序列的長(zhǎng)度;sgn()為符號(hào)函數(shù).當(dāng)|Z|>1-α/2時(shí),表明在水平上存在顯著變化.本文判斷在=0.05置信水平上暴露風(fēng)險(xiǎn)在時(shí)間序列上的顯著性[26].

    1.2.3 遙感數(shù)據(jù)的精度檢驗(yàn) 為全面客觀地評(píng)價(jià)本文使用的PM2.5數(shù)據(jù)精度,采用遙感預(yù)測(cè)值和地面監(jiān)測(cè)值之間的2、和指標(biāo)評(píng)價(jià)PM2.5遙感數(shù)據(jù)的精度.具體含義和計(jì)算方法如下:

    式中:PM2.5為遙感數(shù)值;PM2.5為地面站點(diǎn)監(jiān)測(cè)值;為建模數(shù)據(jù)集記錄總數(shù).

    2 結(jié)果與討論

    2.1 全球PM2.5和人口遙感數(shù)據(jù)的精度檢驗(yàn)

    通過(guò)海量文獻(xiàn)搜索的方式獲得了全球210個(gè)PM2.5地面環(huán)境監(jiān)測(cè)站點(diǎn)數(shù)據(jù),經(jīng)交叉驗(yàn)證顯示驗(yàn)證點(diǎn)處的估算值和實(shí)測(cè)值一致性較好(2=0.81)[27].但是,該產(chǎn)品估算時(shí)是以歐美地區(qū)的實(shí)測(cè)數(shù)據(jù)構(gòu)建模型的,對(duì)于其他地區(qū)的適用性并沒(méi)有進(jìn)行說(shuō)明,加之這些科研團(tuán)隊(duì)和學(xué)者沒(méi)有公布更加詳細(xì)的驗(yàn)證數(shù)據(jù),如和等,該數(shù)據(jù)在全球的精度有待驗(yàn)證.PM2.5數(shù)據(jù)是本文的基礎(chǔ)數(shù)據(jù),為確保PM2.5遙感數(shù)據(jù)的可靠性,分別選取了中國(guó)、歐盟和加拿大為驗(yàn)證區(qū)域,對(duì)2016年三個(gè)區(qū)域的年平均PM2.5進(jìn)行匹配驗(yàn)證.選擇驗(yàn)證點(diǎn)時(shí),以空間分布均勻?yàn)榛驹瓌t,從而保證PM2.5遙感數(shù)據(jù)的檢驗(yàn)結(jié)果客觀可靠.從選取的29個(gè)國(guó)家共2899個(gè)PM2.5監(jiān)測(cè)站點(diǎn)中,篩選出有效的驗(yàn)證數(shù)據(jù)442組,其空間分布如圖1所示,對(duì)其進(jìn)行精度檢驗(yàn).

    2016年中國(guó)共有監(jiān)測(cè)站點(diǎn)1542個(gè),本文在空間上均勻且隨機(jī)地選取了驗(yàn)證點(diǎn)230個(gè),將遙感數(shù)據(jù)與地面實(shí)測(cè)數(shù)據(jù)進(jìn)行匹配驗(yàn)證,在<0.01水平上,2=0.7264,MPE=8.09%,RMSE=10.387μg/m3.歐盟境內(nèi)共有PM2.5監(jiān)測(cè)站點(diǎn)1327個(gè),從中均勻隨機(jī)地選取180個(gè)驗(yàn)證點(diǎn),2=0.798,MPE=1.57%, RMSE=2.685μg/m3.加拿大共有30個(gè)監(jiān)測(cè)站點(diǎn),主要分布在南部,從中選取驗(yàn)證點(diǎn)12個(gè),經(jīng)驗(yàn)證,2= 0.7408,MPE=0.876%,RMSE=1.144μg/m3.將三個(gè)區(qū)域的PM2.5驗(yàn)證結(jié)果進(jìn)行對(duì)比發(fā)現(xiàn),歐盟的2最高,加拿大次之,中國(guó)的2最低;加拿大的MPE和RMSE最優(yōu),歐盟次之,中國(guó)最差.總體來(lái)看,本文所使用的的PM2.5遙感數(shù)據(jù)精度滿足研究需求,精度較好.

    為驗(yàn)證LandScan全球人口格網(wǎng)數(shù)據(jù)的精度,選取世界銀行公布的全球各個(gè)國(guó)家人口總和數(shù)據(jù)對(duì)其進(jìn)行精度檢驗(yàn).使用GIS空間分區(qū)統(tǒng)計(jì)工具,統(tǒng)計(jì)了每個(gè)國(guó)家的人口總量,并和世界銀行公布的人口總和數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證.本文從世界銀行數(shù)據(jù)庫(kù)下載了2016年全球每個(gè)國(guó)家的人口總和,使用ArcMap10.4的區(qū)域統(tǒng)計(jì)工具,統(tǒng)計(jì)了LandScan的各個(gè)國(guó)家人口總和,通過(guò)與世界銀行統(tǒng)計(jì)數(shù)據(jù)進(jìn)行匹配驗(yàn)證,發(fā)現(xiàn)精度良好,<0.001水平下,2=0.99,平均偏差僅為2.43%,可以滿足本研究需要.

    2.2 全球PM2.5和人口空間分布格局

    從圖2可知,2000~2016年全球PM2.5平均濃度在空間上呈亞洲最高、非洲較高、歐洲較低和其他大洲極低的分布狀態(tài),分別為亞洲14.7μg/m3、非洲8.1μg/m3、歐洲8.03μg/m3、南美洲5.69μg/m3、北美洲4.41μg/m3、大洋洲1.27μg/m3.PM2.5現(xiàn)狀值(2016年)分別為亞洲16.6μg/m3、非洲9.15μg/m3、歐洲8.9μg/m3、南美洲5.66μg/m3、北美洲3.6μg/m3、大洋洲1.34μg/m3.可以直觀地發(fā)現(xiàn)亞洲的PM2.5污染程度遠(yuǎn)超其他大洲,尤其以印度和中國(guó)最為嚴(yán)重.中國(guó)的PM2.5濃度高值區(qū)分布在東部,17年間的平均濃度約30μg/m3,超過(guò)WHO的第三級(jí)健康標(biāo)準(zhǔn)約15μg/m3.PM2.5濃度的高值區(qū)主要分布在東北平原中部和以北京-上海-常德為頂點(diǎn)所圍成的三角形區(qū)域內(nèi),這些地區(qū)地形平坦,工業(yè)發(fā)達(dá),產(chǎn)業(yè)集群化程度高,人口密集,大氣污染較嚴(yán)重.東北平原中部的PM2.5濃度整體變化趨勢(shì)為增加,而以北京-上海-常德為頂點(diǎn)所圍成的三角形區(qū)域的變化趨勢(shì)為先增加再減少,2008年,PM2.5濃度最高且污染范圍最廣,2008年以后PM2.5濃度明顯下降.南亞的PM2.5高值區(qū)主要分布在印度和孟加拉國(guó).印度的PM2.5呈現(xiàn)出逐步增加的趨勢(shì),多年均值為33.36μg/m3,到2016年達(dá)到最大值41.05μg/m3.東南亞地區(qū)的PM2.5也呈現(xiàn)出不斷增大的趨勢(shì),其中印度尼西亞和泰國(guó)的17年平均值分別為7.36μg/m3、27.76μg/m3,印度尼西亞的PM2.5增長(zhǎng)尤為劇烈,泰國(guó)的PM2.5多年均值和現(xiàn)狀值均極高.非洲的PM2.5在空間上呈現(xiàn)出中部高、南北低,西部高東部低的分布狀態(tài),高值區(qū)主要分布在幾內(nèi)亞灣附近.非洲的PM2.5多年變化趨勢(shì)為緩慢增加,增加的區(qū)域主要位于中部和南部區(qū)域.

    北美洲PM2.5年均值維持在較低水平,17年P(guān)M2.5平均濃度為4.4μg/m3,其中美國(guó)為6.9μg/m3,加拿大為2.73μg/m3,多年P(guān)M2.5變化呈現(xiàn)出逐步減少的趨勢(shì).南美洲PM2.5年均值濃度為5.69μg/m3,略高于北美.歐洲地區(qū)PM2.5變化則呈現(xiàn)出先增加再減少的趨勢(shì).大洋洲PM2.5常年維持在極低水平,平均值為1.34μg/m3,現(xiàn)狀值為10.6μg/m3.

    從2000~2016年均全球人口空間分布可知(圖3),人口在空間上呈亞洲最高、非洲和歐洲次之、南北美洲較低和大洋洲最低的分布狀態(tài).全球人口持續(xù)增長(zhǎng)的國(guó)家主要分布在非洲,持續(xù)減少的國(guó)家主要分布在歐洲.人口基數(shù)最大的是亞洲,龐大的人口消耗著大量的資源和能源,資源與能源加工生產(chǎn)造成了大量的PM2.5排放,嚴(yán)重的PM2.5污染又反作用于人類(lèi),影響著人類(lèi)的身體健康.非洲人口增長(zhǎng)速度極快.

    圖3 2000~2016年全球人口空間分布

    2.3 全球PM2.5人口暴露風(fēng)險(xiǎn)分布格局

    將PM2.5人口暴露風(fēng)險(xiǎn)分為7個(gè)級(jí)別:極安全、安全、較安全、一般、較危險(xiǎn)、危險(xiǎn)和極危險(xiǎn)[28],詳細(xì)劃分標(biāo)準(zhǔn)如表1所示.

    表1 PM2.5人口暴露風(fēng)險(xiǎn)級(jí)別劃分

    2.3.1 總體分布格局 全球年均PM2.5人口暴露風(fēng)險(xiǎn)呈逐年減少的趨勢(shì),由于全球不同地區(qū)的工業(yè)化和消費(fèi)水平差異顯著,PM2.5人口暴露風(fēng)險(xiǎn)呈現(xiàn)出明顯的空間差異(圖4).PM2.5人口暴露風(fēng)險(xiǎn)的高值區(qū)分布在人口密度大且經(jīng)濟(jì)快速增長(zhǎng)的地區(qū)和國(guó)家,包括中國(guó)、印度、東南亞各國(guó)、非洲尼羅河流域和幾內(nèi)亞灣沿岸、西歐等地,PM2.5人口暴露風(fēng)險(xiǎn)極低值區(qū)主要分布在自然條件惡劣,不適宜人類(lèi)生產(chǎn)發(fā)展的高原、高寒、沙漠、熱帶雨林等地區(qū),如青藏高原、格陵蘭島、東西伯利亞平原、撒哈拉大沙漠、亞馬遜河等.全球PM2.5人口暴露風(fēng)險(xiǎn)在空間上呈亞洲最高(5.94)、非洲次高(0.62)、歐洲、南美洲和北美洲次之(分別為0.45、0.32和0.27),大洋洲最低(0.01)的分布狀態(tài).以2016年為例,PM2.5人口高暴露風(fēng)險(xiǎn)(暴露風(fēng)險(xiǎn)等級(jí)為高風(fēng)險(xiǎn)與極高風(fēng)險(xiǎn))狀態(tài)下的全球各國(guó)家人口總量在空間上呈亞洲和美洲高,歐洲、非洲和大洋洲低的分布狀態(tài).高暴露風(fēng)險(xiǎn)下的人口總量31億的國(guó)家包括中國(guó)、印度、美國(guó)、印度尼西亞、巴西、巴基斯坦、尼日利亞和孟加拉國(guó),暴露人口總量35000萬(wàn)且<1億人的國(guó)家包括俄羅斯、墨西哥、日本和埃及.2016年全球高暴露風(fēng)險(xiǎn)狀態(tài)下的人口密度3200人的國(guó)家包括印度、日本、韓國(guó)、孟加拉國(guó)、荷蘭、盧旺達(dá)和布隆迪,人口密度3100且<200人km2的國(guó)家有13個(gè),包括英國(guó)、德國(guó)、瑞士、比利時(shí)、意大利、巴基斯坦、尼泊爾、朝鮮、越南、菲律賓、多米尼加、薩爾瓦多和尼日利亞.

    為便于分析,將2000~2016年每隔4a統(tǒng)計(jì),PM2.5濃度均值從10~45μg/m3,每隔5μg/m3分級(jí),計(jì)算超過(guò)各級(jí)限值的人口累計(jì)百分比(圖5),結(jié)果發(fā)現(xiàn)超過(guò)各級(jí)限制的暴露人口累計(jì)百分比在2000~2008年有明顯增大的趨勢(shì),2012年顯著回落,2016年又有所回升.2000年全球有73.82%的人口生活在PM2.5濃度年均值為10μg/m3以上的環(huán)境中,有12.59%的人口所生活的環(huán)境PM2.5濃度平均值在35μg/m3以上,還有6.26%的人口生活在PM2.5濃度平均值在45μg/m3以上.2016年,全球仍有75.72%的人口生活在PM2.5濃度年均值為10μg/m3以上的環(huán)境中,27.51%的人口所生活的環(huán)境PM2.5濃度平均值在35μg/m3以上,18.8%的人口暴露在PM2.5濃度平均值在45μg/m3以上.

    2.3.2 分區(qū)域空間分布特征 從圖6可以發(fā)現(xiàn)亞洲的PM2.5人口暴露風(fēng)險(xiǎn)遠(yuǎn)遠(yuǎn)超過(guò)其他大洲,17年的暴露風(fēng)險(xiǎn)均值為5.94,屬于極危險(xiǎn)級(jí)別.中國(guó)PM2.5人口暴露風(fēng)險(xiǎn)高的區(qū)域主要分布在東部地區(qū),平均值為3.21,達(dá)到了較危險(xiǎn)水平.中國(guó)PM2.5人口暴露風(fēng)險(xiǎn)變化趨勢(shì)為先增加后減少,PM2.5濃度最高且污染范圍最廣是2008年,2008年之后暴露風(fēng)險(xiǎn)逐漸降低.日本年均暴露風(fēng)險(xiǎn)值為4.48,屬于危險(xiǎn)等級(jí),高暴露風(fēng)險(xiǎn)區(qū)主要分布在東部沿海地區(qū).南亞的印度、孟加拉國(guó)和巴基斯坦東部區(qū)域是PM2.5人口暴露極危險(xiǎn)區(qū),且呈逐年增加的趨勢(shì),2016年的風(fēng)險(xiǎn)值達(dá)到最大.這三個(gè)國(guó)家的17年均PM2.5人口暴露風(fēng)險(xiǎn)分別為24.03、79.03和8.75.東南亞地區(qū)是全球范圍內(nèi)又一PM2.5人口暴露風(fēng)險(xiǎn)高值區(qū),越南、泰國(guó)、柬埔寨、緬甸、印度尼西亞暴露風(fēng)險(xiǎn)均值分別為8.44、5.89、2.03、2.25、2.76.2000~2016年?yáng)|南亞地區(qū)的PM2.5人口暴露風(fēng)險(xiǎn)呈現(xiàn)出不斷增加的趨勢(shì),其中中南半島的增長(zhǎng)幅度最大.中國(guó)和印度人口眾多,研究時(shí)段內(nèi)都處于經(jīng)濟(jì)飛速發(fā)展的階段,形成了印度北部、東部和中國(guó)華北、華東等高暴露風(fēng)險(xiǎn)地區(qū).日本重工業(yè)主要布局于東部沿海地區(qū),其PM2.5人口暴露風(fēng)險(xiǎn)處于極高風(fēng)險(xiǎn)級(jí)別,但17年間的變化幅度不大.

    圖5 全球PM2.5暴露水平超過(guò)各極限值的人口百分比

    圖6 亞洲PM2.5人口暴露風(fēng)險(xiǎn)空間分布

    非洲的PM2.5人口暴露風(fēng)險(xiǎn)17年均值為0.62,暴露風(fēng)險(xiǎn)等級(jí)為較安全,在空間上呈現(xiàn)出南高北低的空間分布格局,PM2.5人口暴露風(fēng)險(xiǎn)高值區(qū)主要集中在5個(gè)區(qū)域,分別為幾內(nèi)亞灣沿岸北部、埃及境內(nèi)的尼羅河沿岸和三角洲、埃塞俄比亞、環(huán)維多利亞湖和南非的比勒陀利亞地區(qū)(圖7).這5個(gè)區(qū)域是非洲經(jīng)濟(jì)較繁榮的地區(qū),人類(lèi)活動(dòng)頻繁,工業(yè)發(fā)展速度快,形成了高暴露風(fēng)險(xiǎn)區(qū).撒哈拉沙漠暴露風(fēng)險(xiǎn)等級(jí)為極安全,但這個(gè)地區(qū)對(duì)于人口個(gè)體來(lái)說(shuō)是非常危險(xiǎn)的.埃及的PM2.5人口暴露風(fēng)險(xiǎn)均值為1.08,但尼羅河流域的PM2.5人口暴露風(fēng)險(xiǎn)等級(jí)為極危險(xiǎn).幾內(nèi)亞灣是世界石油儲(chǔ)量最豐富的地區(qū)之一,尼日利亞以石油開(kāi)采和加工為支柱的產(chǎn)業(yè)經(jīng)濟(jì)造成了嚴(yán)重的大氣污染,埃及的尼羅河流域、埃塞俄比亞和維多利亞湖周?chē)丝诒姸?采掘業(yè)和加工業(yè)形成了一定規(guī)模,但污染治理無(wú)法跟進(jìn),人口暴露風(fēng)險(xiǎn)很高.

    圖7 非洲PM2.5人口暴露風(fēng)險(xiǎn)空間分布

    圖8 歐洲PM2.5人口暴露風(fēng)險(xiǎn)空間分布

    歐洲PM2.5人口暴露風(fēng)險(xiǎn)均值為0.45,僅次于亞洲和非洲.比利時(shí)、意大利、德國(guó)、英國(guó)、波蘭、法國(guó)和烏克蘭的PM2.5人口暴露風(fēng)險(xiǎn)較高,年均風(fēng)險(xiǎn)值分別為5.35、3.63、3.42、2.81、2.76、1.58和1.29,PM2.5人口暴露風(fēng)險(xiǎn)的變化趨勢(shì)為先增加再減少,2000~ 2006年逐漸增加,2006年之后開(kāi)始緩慢減少,其中PM2.5人口暴露風(fēng)險(xiǎn)高值區(qū)是以各國(guó)的首都為中心呈點(diǎn)狀分布(圖8),這與鄭軍關(guān)于2006~2014年歐盟28國(guó)PM2.5呈下降趨勢(shì)的研究結(jié)論一致[29].歐洲的工業(yè)和制造業(yè)發(fā)達(dá),且總體已進(jìn)入后工業(yè)化時(shí)代,雖然個(gè)別地區(qū)工業(yè)型城市的暴露風(fēng)險(xiǎn)仍較嚴(yán)重,但整體變化趨勢(shì)是減少的.

    南美洲人口總量少,密度低,人口主要分布在西北和東南沿海一帶的區(qū)域,故PM2.5人口暴露風(fēng)險(xiǎn)低,以極安全和安全為主.暴露級(jí)別為危險(xiǎn)與極危險(xiǎn)的區(qū)域只零星分布在亞馬孫河口三角洲、巴西東部和、阿根廷的內(nèi)格羅河與科羅拉河流域.北美洲PM2.5人口暴露風(fēng)險(xiǎn)水平較低,以極安全和安全級(jí)別為主.北美工業(yè)化程度整體較高,工業(yè)污染少,人口密度小,墨西哥、美國(guó)和加拿大的PM2.5人口暴露風(fēng)險(xiǎn)均值分別為1.27、0.41和0.03,現(xiàn)狀暴露風(fēng)險(xiǎn)值為0.98、0.3和0.02,多年暴露均值和現(xiàn)狀暴露值均較低.北美洲高暴露風(fēng)險(xiǎn)區(qū)以特大城市呈團(tuán)簇狀分布,2000~2016年,PM2.5人口暴露風(fēng)險(xiǎn)變化呈現(xiàn)出逐步減少的趨勢(shì).大洋洲PM2.5人口暴露風(fēng)險(xiǎn)維持在極安全和安全級(jí)別.近年來(lái),在澳大利亞中部、南部和新西蘭的暴露風(fēng)險(xiǎn)在逐漸增加,沙塵引起的PM2.5濃度較高.此外,該地區(qū)桉樹(shù)發(fā)生火災(zāi)后釋放大量的可吸收顆粒物,使得PM2.5人口暴露風(fēng)險(xiǎn)增加[30].

    2.4 PM2.5人口暴露時(shí)空變化特征

    2.4.1 暴露風(fēng)險(xiǎn)線性變化趨勢(shì) 利用Theil-Sen Median和Mann-Kendall方法刻畫(huà)PM2.5人口暴露在時(shí)間序列上的顯著性變化.當(dāng)>0時(shí),表明PM2.5人口暴露在時(shí)間序列上呈增加趨勢(shì).當(dāng)<0時(shí),PM2.5人口暴露則呈減少趨勢(shì).對(duì)檢驗(yàn)結(jié)果的值分別進(jìn)行分級(jí),將值的絕對(duì)值結(jié)果按照顯著性分為極顯著變化(||>2.58)、顯著變化(1.96<||£2.58)、弱顯著變化(1.65<||£1.96)和無(wú)顯著變化(0<||£1.65)4個(gè)等級(jí)[26,31].將值與||值重分類(lèi)后再進(jìn)行同像元疊加,得到全球PM2.5人口暴露線性變化趨勢(shì)(圖9).

    圖9 全球PM2.5人口暴露線性變化趨勢(shì)

    由圖9可見(jiàn)2000~2016年全球PM2.5人口暴露風(fēng)險(xiǎn)整體上呈增加趨勢(shì),極顯著增加的區(qū)域主要分布在亞洲和非洲,顯著增加的區(qū)域分布在極顯著增加的周?chē)?極顯著減少的區(qū)域主要分布在西歐和美國(guó).亞洲極顯著增長(zhǎng)幅度大的區(qū)域主要為南亞、中國(guó)東北、青藏高原、帕米爾高原、里海地區(qū)和東南亞其中中南半島增長(zhǎng)幅度最大.非洲的變化特征以極顯著增加為主,主要包括剛果、贊比亞、安哥拉、坦桑尼亞、贊比亞、南非、馬達(dá)加斯加島、尼日爾南部、尼日利亞北部、布基納法索、乍得西部和阿爾及利亞中部地區(qū)等,顯著增加與弱顯著增加的區(qū)域環(huán)繞在極顯著增加區(qū)域的周?chē)?非洲極顯著減少的區(qū)域分布在幾內(nèi)亞灣北部、西撒哈拉、埃及和利比亞.

    北美洲的變化趨勢(shì)呈兩極化.美國(guó)以極顯著減少為主要變化特征,而加拿大西部和中部地區(qū)呈現(xiàn)出大范圍的極顯著增加區(qū).南美洲變化特征以增加為主,減少的區(qū)域極少,變化趨勢(shì)為極顯著增加的區(qū)域包括亞馬孫河口三角洲、巴西東部、阿根廷的內(nèi)格羅河與科羅拉河流域.歐洲基本上均為顯著減少區(qū).澳大利亞中部和西南部出現(xiàn)了極顯著增加區(qū).

    2.4.2 PM2.5人口暴露風(fēng)險(xiǎn)時(shí)間變化 由圖10可知,各大洲變化差異顯著.亞洲呈增長(zhǎng)趨勢(shì),其暴露風(fēng)險(xiǎn)值明顯高于其他大洲,最低值在2000年,暴露風(fēng)險(xiǎn)值為3.55,最高值出現(xiàn)在2008年,暴露值為3.87.非洲的暴露風(fēng)險(xiǎn)值也呈增加趨勢(shì),最低值為2009年的0.33,最高值為0.45.歐洲和北美洲呈減少趨勢(shì),歐洲和北美洲的暴露風(fēng)險(xiǎn)最高值均出現(xiàn)在2000年,分別為0.39和0.24,最低值分別為2008年的0.23和2016年的0.13.大洋洲和南美洲暴露風(fēng)險(xiǎn)值較低且變化也較小,尤其大洋洲的暴露風(fēng)險(xiǎn)值始終維持在極低的狀態(tài).

    圖10 2000~2016年全球年均PM2.5人口暴露風(fēng)險(xiǎn)時(shí)間變化

    2.5 討論

    傳統(tǒng)的PM2.5人口暴露風(fēng)險(xiǎn)是基于行政單元進(jìn)行研究,時(shí)空連續(xù)性方面存在缺陷.研究依賴于人口普查數(shù)據(jù)和環(huán)境監(jiān)測(cè)站點(diǎn)的PM2.5數(shù)據(jù),是基于行政單元進(jìn)行分析,忽視了其內(nèi)部的空間異質(zhì)性,以行政單元為主體研究PM2.5人口暴露風(fēng)險(xiǎn),割裂了空間單元,不利于宏觀上的污染物防治[32-33].然而,對(duì)于環(huán)境監(jiān)控網(wǎng)絡(luò)來(lái)說(shuō),要獲取空間連續(xù)的監(jiān)測(cè)值,在全球尺度是難以實(shí)現(xiàn)的,而對(duì)于人口普查數(shù)據(jù)來(lái)說(shuō),可信度雖高,但10年或更長(zhǎng)時(shí)間更新一次的時(shí)間分辨率不僅無(wú)法滿足實(shí)時(shí)變化的暴露風(fēng)險(xiǎn)研究,還不能與PM2.5數(shù)據(jù)相匹配.該研究以遙感手段,基于柵格量化了PM2.5人口暴露風(fēng)險(xiǎn).分析了宏觀尺度、長(zhǎng)時(shí)間序列的PM2.5濃度和PM2.5人口暴露風(fēng)險(xiǎn)的變化.從全球宏觀視角分析了PM2.5人口暴露現(xiàn)狀,識(shí)別高暴露風(fēng)險(xiǎn)區(qū),剖析高暴露污染的時(shí)空變化規(guī)律,探究了各國(guó)高暴露的原因,為降低暴露風(fēng)險(xiǎn)和跨國(guó)跨區(qū)域大氣污染聯(lián)動(dòng)治理提供時(shí)空連續(xù)的基礎(chǔ)數(shù)據(jù).許多發(fā)展中國(guó)家受經(jīng)濟(jì)條件和國(guó)情需要,PM2.5監(jiān)測(cè)歷史較短,例如,中國(guó)2012年大氣環(huán)境監(jiān)測(cè)系統(tǒng)才將PM2.5的濃度值納入監(jiān)測(cè)指標(biāo)體系,此前全國(guó)PM2.5監(jiān)測(cè)處于空白,且無(wú)公開(kāi)歷史數(shù)據(jù)[34].2000~2016年長(zhǎng)時(shí)間序列的PM2.5人口暴露分析,彌補(bǔ)了部分國(guó)家缺失的數(shù)據(jù).該研究采用了目前可獲取的更高空間分辨率的PM2.5和人口分布柵格數(shù)據(jù),基于柵格尺度,利用污染物人口暴露風(fēng)險(xiǎn)法、Theil-Sen Median和Mann-Kendall等方法,客觀地刻畫(huà)PM2.5、人口、PM2.5人口暴露風(fēng)險(xiǎn)的時(shí)空格局,準(zhǔn)確掌握高暴露風(fēng)險(xiǎn)區(qū)域,還可以直觀地測(cè)度各空間柵格在空間單元整體PM2.5人口暴露中的貢獻(xiàn)大小.

    本文PM2.5人口暴露風(fēng)險(xiǎn)指數(shù)是結(jié)合PM2.5空間分布與人口空間分布的綜合指數(shù),其風(fēng)險(xiǎn)程度由二者共同決定,撒哈拉與塔克拉瑪干等沙漠地區(qū),有時(shí)受大風(fēng)沙塵天氣影響PM2.5濃度較高,但自然條件惡劣,人跡罕至,幾乎無(wú)人類(lèi)活動(dòng)和經(jīng)濟(jì)活動(dòng),故暴露風(fēng)險(xiǎn)級(jí)別是安全的,但這些區(qū)域的瞬時(shí)暴露對(duì)于人口個(gè)體來(lái)說(shuō)依然是非常危險(xiǎn)的.受數(shù)據(jù)獲取的限制,缺少日、月和季節(jié)平均的PM2.5和人口數(shù)據(jù),因此無(wú)法分析PM2.5人口暴露的日、月和季節(jié)差異.此外,可獲取的人口柵格數(shù)據(jù)缺少屬性,故無(wú)法區(qū)分不同年齡、不同職業(yè)的人口暴露風(fēng)險(xiǎn)差異.

    3 結(jié)論

    3.1 全球PM2.5時(shí)空格局具有顯著差異.PM2.5污染最嚴(yán)重的區(qū)域是印度和中國(guó)東部,污染持續(xù)增長(zhǎng)的區(qū)域主要分布在印度和東南亞.人口在空間上呈亞洲最高、非洲歐洲次之、南北美洲較低和大洋洲最低的分布狀態(tài).全球人口持續(xù)增長(zhǎng)的國(guó)家主要分布在非洲,持續(xù)減少的國(guó)家主要分布在歐洲.

    3.2 全球PM2.5人口暴露風(fēng)險(xiǎn)存在顯著的時(shí)空差異.全球PM2.5人口暴露風(fēng)險(xiǎn)在宏觀尺度上呈逐漸減少的趨勢(shì),空間上呈亞洲最高、非洲次高、歐洲較低和其他大洲極低的分布狀態(tài).PM2.5人口暴露風(fēng)險(xiǎn)的高值區(qū)分布在中國(guó)、印度、東南亞等經(jīng)濟(jì)高速發(fā)展且人口密集的區(qū)域,而PM2.5人口暴露風(fēng)險(xiǎn)極低值區(qū)主要分布在高山、高寒、沙漠、熱帶雨林等人跡罕至,開(kāi)發(fā)程度極低的區(qū)域.極顯著增加的區(qū)域主要有亞洲、非洲、南美洲和加拿大西部,極顯著減少的區(qū)域主要有美國(guó)和西歐.

    [1] Ott W R. Concepts of human exposure to air pollution [J]. Environment International, 1982,7(3):179-196.

    [2] Gakidou E, Afshin A, Abajobir A A, et al. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990~2016: a systematic analysis for the Global Burden of Disease Study 2016 [J]. The Lancet, 2017,390(10100):1345–1422.

    [3] Stanaway J D, Afshin A, Gakidou E, et al. Global, regional, and national comparative risk assessment of 84 behavioral, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990~2017: a systematic analysis for the Global Burden of Disease Study 2017 [J]. Lancet, 392(2018):1923– 1994,10.1016/s0140-6736(18)32225-6.

    [4] WHO. Ambient (outdoor) air quality and health. World Health Organization; 2018-05-02. http://www.who.int/mediacentre/factsheets/ fs313/en/.Accessed 9Sept 2019.

    [5] Hajat A, Allison M, Diez-Roux A V, et al. Long-term exposure to air pollution and markers of inflammation, coagulation, and endothelial activation: A repeat-measures analysis in the multi-ethnic study of atherosclerosis (MESA) [J]. Epidemiology, 2015,26(3):310-320.

    [6] Zhang L L, Pan J H. Estimation of PM2.5mass concentrations in Beijing–Tianjin–Hebei Region based on GWR and spatial downscaling method[J]. Journal of the Indian Society of Remote Sensing, 2020,48(12):1691-1703.

    [7] Nansai K, Tohno S, Chatani S, et al. Affluent countries inflict inequitable mortality and economic loss on Asia via PM2.5emissions [J]. Environment International, 134:105238.

    [8] 慕 航,何 超,阮秋明,等.“一帶一路”沿線國(guó)家PM2.5污染與人口暴露風(fēng)險(xiǎn)的時(shí)空分布特征[J]. 環(huán)境科學(xué)學(xué)報(bào), 2021,41(6):2229-2240.

    Mu H, He C, Ruan Q M, et al. Spatiotemporal distribution characteristics and population exposure risks to PM2.5in countries along the Belt and Road [J]. Acta Scientiae Circumstantiae, 2021,41(6):2229-2240.

    [9] Zhang Q, Jiang X, Tong D, et al. Transboundary health impacts of transported global air pollution and international trade [J]. Nature, 543(7647):705-709.

    [10] Hamra G B, Guha N, Cohen A, et al. Outdoor particulate matter exposure and lung cancer: a systematic review and meta-analysis [J]. Environmental Health Perspect, 2014,122(9):906–911.

    [11] Pope C A, Brook R D, Burnett R T, et al. How is cardiovascular disease mortality risk affected by duration and intensity of fine particulate matter exposure? An integration of the epidemiologic evidence [J]. Air Quality, Atmosphere & Health, 2010,4(1):5-14.

    [12] Stafoggia M, Cesaroni G, Peters A. Long-term exposure to ambient air pollution and incidence of cerebrovascular events: Results from 11European cohorts within the ESCAPE project [J]. Environ. Health Perspect, 2014,122(9):919-925.

    [13] Saraswat A, Kandlikar M, Brauer M, et al. PM2.5population exposure in New Delhi using a probabilistic simulation framework [J]. Environmental Science & Technology, 2016,50(6):3174-3183.

    [14] Han L, Zhou W, Pickett S, et al. Risks and causes of population exposure to cumulative fine particulate (PM2.5) pollution in China [J]. Earth's Future, 2019,7(6):615–622.

    [15] Saini P, Sharma M. Cause and age-specific premature mortality attributable to PM2.5exposure: An analysis for million-plus Indian cities [J]. Science of The Total Environment, 2019,710:135230.

    [16] Wu W, Yao M, Yang X, et al. Mortality burden attributable to long- term ambient PM2.5exposure in China: using novel exposure-response functions with multiple exposure windows [J]. Atmospheric Environment, 2020,246(10):118098.

    [17] Yang B Y, Qianmin Z, Howard S W, et al. Global association between ambient air pollution and blood pressure: A systematic review and meta-analysis [J]. Environmental Pollution, 2018,235:576-588.

    [18] Han L, Zhou W, Li W, et al. Global population exposed to fine particulate pollution by population increase and pollution expansion. Air Quality Atmosphere & Health, 2017,10:1221–1226.

    [19] Roberts G, Wooster M J. Global impact of landscape fire emissions on surface level PM2.5concentrations, air quality exposure and population mortality [J]. Atmospheric Environment, 2021,(4):118210.

    [20] Li J M, Han X L, Jin M J, et al. Globally analysing spatiotemporal trends of anthropogenic PM2.5concentration and population's PM2.5exposure from 1998 to 2016 [J]. Environment International, 2019,128:46-62.

    [21] Zhang L L, Pan J H. Estimation of PM2.5mass concentrations in Beijing–Tianjin–Hebei Region based on geographically weighted regression and spatial downscaling method [J]. Indian Society of Remote Sensing. 2020,48(12):1691–1703.

    [22] 張亮林,潘竟虎.中國(guó)PM2.5人口暴露風(fēng)險(xiǎn)時(shí)空格局[J]. 中國(guó)環(huán)境科學(xué), 2020,40(1):1-12.

    Zhang L L, Pan J H. Spatial-temporal pattern of population exposure risk to PM2.5in China [J].China Environmental Science, 2020,40(1):1-12.

    [23] 鄒 濱,彭 芬,焦利民,等.高分辨率人口空氣污染暴露GIS空間區(qū)劃研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2013,38(3):334-338.

    Zou B, Peng F, Jiao L M, et al. GIS aided spatial zoning of high- resolution population exposure to air pollution [J]. Geomatics and Information Science of Wuhan University, 2013,38(3):334-338.

    [24] Sen P K. Estimates of the regressions coefficient based on Kendall’s tau [J]. Journal of the American Statistical Association, 1968,63(324): 1379-1389.

    [25] Kendall M G. Rank correlation methods [J]. British Journal of Psychology, 1955,25(1):86-91.

    [26] 張永年,潘竟虎.基于DMSP/OLS數(shù)據(jù)的中國(guó)碳排放時(shí)空模擬與分異格局[J]. 中國(guó)環(huán)境科學(xué), 2019,39(4):1436-1446.

    Zhang Y N, Pan J H. Spatio-temporal simulation and differentiation pattern of carbon emissions in China based on DMSP/OLS nighttime light data [J]. China Environmental Science, 2019,39(4):1436-1446.

    [27] Zhang Y L, Cao F. Fine particulate matter PM2.5in China at a city level [J]. Scientific Reports, 2015,5(1):14884.

    [28] 別 同,韓立建,田淑芳,等.城市化對(duì)空氣污染人群暴露貢獻(xiàn)的定量方法研究[J]. 生態(tài)學(xué)報(bào), 2018,38(13):4570-4583.

    Bie T, Han L J, Tian S F, et al. Method for quantifying the contribution of urbanization on population exposure to air pollution [J]. Acta Ecologica Sinica, 2018,38(13):4570-4583.

    [29] 鄭 軍,袁 鈺.歐洲大氣污染防治經(jīng)驗(yàn)及對(duì)我國(guó)的啟示[J]. 環(huán)境保護(hù), 2017,45(18):64-66.

    [30] Kolbe A, Gilchrist K L. An extreme bushfire smoke pollution event: Health impacts and public health challenges [J]. New South Wales Public Health Bulletin, 2009,20(2):19-23.

    [31] 徐建華.現(xiàn)代地理學(xué)中的數(shù)學(xué)方法[M]. 北京:高等教育出版社, 2016:133-136.

    Xu J H. Mathematical methods in contemporary geography [M]. Beijing: Higher Education Press, 2016:133-136.

    [32] 李名升,任曉霞,于 洋,等.中國(guó)大陸城市PM2.5污染時(shí)空分布規(guī)律[J]. 中國(guó)環(huán)境科學(xué), 2016,36(3):641-65.

    Li M S, Ren X X, Yu Y, et al. Spatiotemporal pattern of ground-level fine particulate matter (PM2.5) pollution in mainland China [J]. China Environmental Science,2016,36(3):641-650.

    [33] 黃曉軍,祁明月,李艷雨,等.關(guān)中地區(qū)PM2.5時(shí)空演化及人口暴露風(fēng)險(xiǎn)[J]. 環(huán)境科學(xué), 2020,41(12):5245-5255.

    Huang X J, Qi M Y, Li Y Y, et al. Spatio-temporal Evolution and Population Exposure Risk to PM2.5in the Guanzhong Area [J]. Environmental Science, 2020,41(12):5245-5255.

    [34] 何映月,張松林.2000~2016年中國(guó)PM2.5時(shí)空格局[J]. 中國(guó)環(huán)境科學(xué), 2020,40(8):3284-3293.

    He Y Y, Zhong S L. Spatiotemporal pattern of PM2.5from 2000 to 2016 in China [J]. China Environmental Sciencece, 2020,40(8):3284-3293.

    Spatial-temporal pattern of population exposure risk to PM2.5in Global.

    ZHANG Liang-lin, PAN Jing-hu*

    (College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, China)., 2021,41(11):5391~5404

    Based on PM2.5remote sensing data and population grid data, the spatial-temporal distribution characteristics of global population exposure risk to PM2.5from 2000 to 2016 were analyzed by using exposure risk model. Theil-Sen Median, Mann-Kendall, and the high-risk areas were accurately identified. The results show that PM2.5remote sensing data and population grid data had good accuracy. China, European Union and Canada were selected to verify the PM2.5mass concentrations and population grid data with good accuracy. The global average PM2.5mass concentrations varies significantly among different continents, and the high-value PM2.5pollution regions are mainly distributed in East Asia, South Asia and Southeast Asia. The annual mean PM2.5mass concentrations ranged from 14.7μg/m3in Asia, 8.1μg/m3in Africa, 8.03μg/m3in Europe, 5.69μg/m3in South America, 4.41μg/m3in North America and 1.27μg/m3in Oceania, respectively. The population exposure risk to PM2.5in the world showed a gradually decreasing trend in the macro scale, while it showed a different trend in the region. In terms of spatial scale, the population exposure risks to PM2.5in all continents rank from high to low in Asia, 0.62 in Africa, 0.45 in Europe, 0.32 in South America, 0.27 in North America and 0.01 in Oceania. Time series, global population exposure risk to PM2.5is significantly different from 2000 to 2016 years. Asia and Africa showed an increasing trend, Europe and North America showed a decreasing trend, Oceania and South America showed a small range of change.

    PM2.5;population exposure risk;spatial-temporal evolution;Global

    X511

    A

    1000-6923(2021)11-5391-14

    張亮林(1993-),男,甘肅慶陽(yáng)人,西北師范大學(xué)碩士研究生,主要從事大氣環(huán)境遙感方面的研究.發(fā)表論文6篇.

    2021-04-02

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41661025,42071216)

    * 責(zé)任作者, 教授, panjh@nwnu.edu.cn

    猜你喜歡
    人口區(qū)域濃度
    拍照濃度有待降低
    《世界人口日》
    人口轉(zhuǎn)型為何在加速 精讀
    神醫(yī)的煩惱——濃度與配比
    人口最少的國(guó)家
    1723 萬(wàn)人,我國(guó)人口數(shù)據(jù)下滑引關(guān)注
    關(guān)于四色猜想
    分區(qū)域
    對(duì)改良的三種最小抑菌濃度測(cè)試方法的探討
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    国产av不卡久久| 最近视频中文字幕2019在线8| 男女边吃奶边做爰视频| 欧美成人一区二区免费高清观看| 狂野欧美白嫩少妇大欣赏| 一级av片app| 国产精品麻豆人妻色哟哟久久 | 国产一级毛片在线| 久久这里有精品视频免费| 国产伦在线观看视频一区| 午夜福利在线观看免费完整高清在| 自拍偷自拍亚洲精品老妇| 久久久成人免费电影| 中国国产av一级| 午夜福利在线观看吧| 国产精品一区www在线观看| 老司机福利观看| 欧美变态另类bdsm刘玥| 成人毛片60女人毛片免费| 精品久久久久久久久亚洲| 国产亚洲91精品色在线| 久久国内精品自在自线图片| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看| 国产免费一级a男人的天堂| 日本午夜av视频| АⅤ资源中文在线天堂| 国语对白做爰xxxⅹ性视频网站| 国产黄色视频一区二区在线观看 | 亚洲av免费在线观看| 免费观看人在逋| 精品人妻熟女av久视频| 日韩视频在线欧美| 在线观看66精品国产| 神马国产精品三级电影在线观看| 国产成人精品婷婷| 在现免费观看毛片| 国产在视频线在精品| 伦理电影大哥的女人| 精品人妻视频免费看| 有码 亚洲区| 欧美性猛交黑人性爽| 99在线人妻在线中文字幕| 亚洲久久久久久中文字幕| av在线播放精品| 观看免费一级毛片| 最近2019中文字幕mv第一页| 亚洲四区av| 国产真实乱freesex| 成人鲁丝片一二三区免费| 99热这里只有精品一区| 精品久久久噜噜| 久久亚洲国产成人精品v| .国产精品久久| 国产爱豆传媒在线观看| 日韩中字成人| 国产精品一区二区三区四区久久| 观看美女的网站| 日本欧美国产在线视频| 一级av片app| 国产成人a区在线观看| 久久久久久久久久久免费av| 黄片wwwwww| 亚洲av熟女| 狂野欧美白嫩少妇大欣赏| 91aial.com中文字幕在线观看| 欧美一区二区国产精品久久精品| 亚洲人与动物交配视频| 色5月婷婷丁香| 神马国产精品三级电影在线观看| 久久综合国产亚洲精品| 中文字幕免费在线视频6| 黄色配什么色好看| 高清日韩中文字幕在线| 最近手机中文字幕大全| 国产成人a区在线观看| 1000部很黄的大片| 中国国产av一级| 啦啦啦啦在线视频资源| 黄色日韩在线| 床上黄色一级片| 久99久视频精品免费| 秋霞伦理黄片| 蜜臀久久99精品久久宅男| 久久综合国产亚洲精品| 精品人妻一区二区三区麻豆| 欧美日韩一区二区视频在线观看视频在线 | 日韩欧美 国产精品| 国产黄a三级三级三级人| 三级毛片av免费| 啦啦啦韩国在线观看视频| 午夜福利高清视频| 色吧在线观看| 久99久视频精品免费| 中文字幕av在线有码专区| 欧美精品国产亚洲| 中文字幕精品亚洲无线码一区| 亚洲成av人片在线播放无| 日产精品乱码卡一卡2卡三| 中国国产av一级| 美女高潮的动态| 人妻制服诱惑在线中文字幕| 99热网站在线观看| 精品一区二区三区视频在线| 久久久久久久久中文| 日本五十路高清| 天美传媒精品一区二区| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩一区二区视频在线观看视频在线 | 天堂网av新在线| 免费观看在线日韩| 黄片无遮挡物在线观看| 内射极品少妇av片p| 最近中文字幕高清免费大全6| 国产精品三级大全| 成人午夜高清在线视频| 久99久视频精品免费| 韩国高清视频一区二区三区| 欧美成人精品欧美一级黄| 国产在视频线在精品| 久久久a久久爽久久v久久| 黄色日韩在线| 国产精品国产三级专区第一集| 国产精品av视频在线免费观看| 男人舔女人下体高潮全视频| 一夜夜www| 岛国在线免费视频观看| 精品一区二区三区视频在线| 久久精品久久久久久久性| 一级黄片播放器| 成年av动漫网址| 少妇的逼水好多| 晚上一个人看的免费电影| ponron亚洲| 亚洲精品久久久久久婷婷小说 | 久久久亚洲精品成人影院| 欧美区成人在线视频| 如何舔出高潮| 国产精品国产三级国产专区5o | 国产成人精品一,二区| 高清在线视频一区二区三区 | 国内少妇人妻偷人精品xxx网站| 成人三级黄色视频| 中文天堂在线官网| 欧美3d第一页| 午夜福利成人在线免费观看| 日韩强制内射视频| 97在线视频观看| 久久久久久久午夜电影| 国产精品嫩草影院av在线观看| 国产一级毛片在线| 九九爱精品视频在线观看| 国产在线男女| 色综合色国产| 天天躁夜夜躁狠狠久久av| 久久精品影院6| 乱人视频在线观看| 在线免费十八禁| 乱系列少妇在线播放| 99热精品在线国产| 99久久九九国产精品国产免费| 自拍偷自拍亚洲精品老妇| 亚洲av成人av| 国产真实乱freesex| 国产 一区 欧美 日韩| 黄色配什么色好看| 免费观看a级毛片全部| 国产成人午夜福利电影在线观看| 国产午夜精品久久久久久一区二区三区| 国产高潮美女av| 亚洲精品日韩在线中文字幕| av免费在线看不卡| 男女视频在线观看网站免费| 日韩一区二区三区影片| 天堂影院成人在线观看| 在线播放国产精品三级| 免费看av在线观看网站| 天堂网av新在线| 国产精品久久久久久久电影| 久久久久精品久久久久真实原创| 免费看美女性在线毛片视频| 三级毛片av免费| 成人欧美大片| av在线蜜桃| 免费电影在线观看免费观看| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美精品专区久久| 男女边吃奶边做爰视频| 国产精品熟女久久久久浪| 老女人水多毛片| 三级经典国产精品| 丰满乱子伦码专区| 国产淫片久久久久久久久| 国产成人午夜福利电影在线观看| 99热网站在线观看| 中文字幕精品亚洲无线码一区| 国产黄片美女视频| 五月玫瑰六月丁香| 亚洲国产成人一精品久久久| 纵有疾风起免费观看全集完整版 | 亚洲国产最新在线播放| 国产伦一二天堂av在线观看| 日韩欧美三级三区| 亚洲久久久久久中文字幕| ponron亚洲| 色播亚洲综合网| 国产高清三级在线| 国产不卡一卡二| 亚洲欧美中文字幕日韩二区| 99久久九九国产精品国产免费| 床上黄色一级片| 欧美极品一区二区三区四区| 日韩av不卡免费在线播放| 国产成人a∨麻豆精品| 国产激情偷乱视频一区二区| 亚洲国产精品专区欧美| 久久久国产成人精品二区| 国产一级毛片七仙女欲春2| 国产精品女同一区二区软件| 一级爰片在线观看| 免费观看a级毛片全部| 1024手机看黄色片| 亚洲av二区三区四区| 久久精品综合一区二区三区| 国产av一区在线观看免费| 汤姆久久久久久久影院中文字幕 | 国产麻豆成人av免费视频| 国产精品1区2区在线观看.| 综合色av麻豆| 久久精品熟女亚洲av麻豆精品 | 国产精品人妻久久久影院| 日本黄色片子视频| 天天躁夜夜躁狠狠久久av| 亚洲婷婷狠狠爱综合网| 伊人久久精品亚洲午夜| 淫秽高清视频在线观看| 大又大粗又爽又黄少妇毛片口| 精华霜和精华液先用哪个| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看| 国产精品爽爽va在线观看网站| 国产又黄又爽又无遮挡在线| 色噜噜av男人的天堂激情| 国产黄色视频一区二区在线观看 | 久久精品夜色国产| 午夜福利在线观看免费完整高清在| 精品免费久久久久久久清纯| 国产成人午夜福利电影在线观看| 亚洲熟妇中文字幕五十中出| 国产久久久一区二区三区| 亚洲精品乱码久久久久久按摩| 精品久久久久久久末码| 久久亚洲国产成人精品v| 亚洲熟妇中文字幕五十中出| 又爽又黄a免费视频| 日韩欧美精品v在线| 久久久色成人| 精品国产一区二区三区久久久樱花 | 一区二区三区乱码不卡18| 熟妇人妻久久中文字幕3abv| 精品久久国产蜜桃| 国产日韩欧美在线精品| 最近最新中文字幕免费大全7| 国产精品久久久久久精品电影小说 | 中文字幕亚洲精品专区| 夜夜爽夜夜爽视频| 午夜福利高清视频| 中国国产av一级| 亚洲国产成人一精品久久久| 免费av不卡在线播放| 久久久久久伊人网av| 一本久久精品| 国产成人a区在线观看| 边亲边吃奶的免费视频| АⅤ资源中文在线天堂| 日本黄色片子视频| 三级经典国产精品| 久久人人爽人人爽人人片va| 久久人人爽人人片av| 国产成人freesex在线| 日本wwww免费看| 三级毛片av免费| 国产午夜福利久久久久久| 男女视频在线观看网站免费| 成人高潮视频无遮挡免费网站| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 亚洲国产精品sss在线观看| 水蜜桃什么品种好| 边亲边吃奶的免费视频| 欧美色视频一区免费| 精品不卡国产一区二区三区| 亚洲av熟女| 欧美日韩国产亚洲二区| 丰满乱子伦码专区| 国产精品久久电影中文字幕| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩无卡精品| 国产视频内射| 只有这里有精品99| 好男人在线观看高清免费视频| av天堂中文字幕网| 超碰av人人做人人爽久久| 国产伦精品一区二区三区四那| 七月丁香在线播放| videos熟女内射| 国产一区二区在线观看日韩| 久热久热在线精品观看| 日韩欧美精品v在线| 非洲黑人性xxxx精品又粗又长| 1000部很黄的大片| 久久久国产成人免费| 日本黄色视频三级网站网址| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 69人妻影院| 亚洲成av人片在线播放无| 麻豆国产97在线/欧美| 男女那种视频在线观看| 亚洲伊人久久精品综合 | 欧美不卡视频在线免费观看| 国产黄片视频在线免费观看| 国产探花极品一区二区| 日本av手机在线免费观看| 人人妻人人看人人澡| 变态另类丝袜制服| 99热精品在线国产| 亚洲欧洲日产国产| 免费一级毛片在线播放高清视频| 又粗又硬又长又爽又黄的视频| 观看免费一级毛片| 国产精品人妻久久久影院| 亚洲天堂国产精品一区在线| 亚洲综合色惰| 真实男女啪啪啪动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品爽爽va在线观看网站| 中文字幕久久专区| 国产亚洲午夜精品一区二区久久 | 别揉我奶头 嗯啊视频| av播播在线观看一区| 久久人人爽人人片av| 久久精品人妻少妇| 久久精品国产亚洲网站| 国产精品国产三级专区第一集| 欧美xxxx性猛交bbbb| 在线观看av片永久免费下载| 国产美女午夜福利| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 国产爱豆传媒在线观看| 最近的中文字幕免费完整| 边亲边吃奶的免费视频| 热99re8久久精品国产| 国产乱人视频| a级毛色黄片| 久久欧美精品欧美久久欧美| 国产亚洲一区二区精品| 精华霜和精华液先用哪个| 91精品国产九色| 哪个播放器可以免费观看大片| 色视频www国产| 亚洲精品久久久久久婷婷小说 | 国产精品一区二区性色av| 一本一本综合久久| 日本免费一区二区三区高清不卡| 国产麻豆成人av免费视频| 国产乱人视频| 久久99热6这里只有精品| 2021天堂中文幕一二区在线观| 精品人妻熟女av久视频| 国产av在哪里看| 国产亚洲一区二区精品| 欧美日韩精品成人综合77777| 99久国产av精品| 成人三级黄色视频| 亚洲美女搞黄在线观看| 国产精品久久久久久精品电影| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 午夜精品一区二区三区免费看| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| 九九热线精品视视频播放| 亚洲精华国产精华液的使用体验| 国产黄色小视频在线观看| 不卡视频在线观看欧美| 国产精品久久久久久av不卡| 小说图片视频综合网站| 日韩亚洲欧美综合| 精品人妻偷拍中文字幕| 久久精品国产亚洲网站| 日本午夜av视频| 99热这里只有精品一区| 99久久中文字幕三级久久日本| 国产高清三级在线| 成年女人永久免费观看视频| 色5月婷婷丁香| 久久99热这里只频精品6学生 | 噜噜噜噜噜久久久久久91| 国产又色又爽无遮挡免| 午夜亚洲福利在线播放| 国产精品一区二区三区四区免费观看| 性色avwww在线观看| 国产av码专区亚洲av| 亚洲欧美清纯卡通| 国产爱豆传媒在线观看| 国产69精品久久久久777片| 欧美日韩在线观看h| 一级av片app| 最新中文字幕久久久久| av黄色大香蕉| 国产成人福利小说| 中国美白少妇内射xxxbb| 免费看美女性在线毛片视频| 看十八女毛片水多多多| 亚洲va在线va天堂va国产| 身体一侧抽搐| 99热全是精品| 久久精品久久久久久噜噜老黄 | 久久久久久久午夜电影| 精品少妇黑人巨大在线播放 | 欧美精品一区二区大全| 免费在线观看成人毛片| av免费观看日本| 中文字幕免费在线视频6| 免费av毛片视频| 亚洲av成人精品一区久久| 国产精品三级大全| 91精品一卡2卡3卡4卡| 成年女人永久免费观看视频| av在线老鸭窝| 午夜a级毛片| 久久久久性生活片| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 又爽又黄a免费视频| 哪个播放器可以免费观看大片| 最近2019中文字幕mv第一页| 三级男女做爰猛烈吃奶摸视频| 乱人视频在线观看| 联通29元200g的流量卡| 中文精品一卡2卡3卡4更新| 麻豆一二三区av精品| 最近视频中文字幕2019在线8| 日韩,欧美,国产一区二区三区 | 亚洲精品aⅴ在线观看| 亚洲成人久久爱视频| 国产极品天堂在线| 免费黄网站久久成人精品| 久久久久久久久大av| 国产色婷婷99| 免费看美女性在线毛片视频| 亚洲欧美成人精品一区二区| 九九在线视频观看精品| 一区二区三区免费毛片| 欧美精品一区二区大全| 国产一级毛片在线| 国产私拍福利视频在线观看| 成人午夜精彩视频在线观看| 亚洲av熟女| 高清午夜精品一区二区三区| 99热6这里只有精品| 国产一区二区在线观看日韩| 亚洲av成人精品一区久久| 淫秽高清视频在线观看| 国产精品.久久久| 国内精品宾馆在线| 日本与韩国留学比较| 亚洲精品乱码久久久v下载方式| 美女黄网站色视频| 欧美高清成人免费视频www| 69人妻影院| 高清午夜精品一区二区三区| 少妇熟女欧美另类| 热99在线观看视频| 成人欧美大片| av在线播放精品| av.在线天堂| 免费在线观看成人毛片| 婷婷六月久久综合丁香| 国产伦精品一区二区三区四那| 日本av手机在线免费观看| 欧美区成人在线视频| 国产麻豆成人av免费视频| eeuss影院久久| 一级爰片在线观看| 一级毛片久久久久久久久女| 精品国产一区二区三区久久久樱花 | 建设人人有责人人尽责人人享有的 | 久久6这里有精品| 天堂中文最新版在线下载 | 综合色丁香网| 小说图片视频综合网站| 秋霞在线观看毛片| 精品国产三级普通话版| 日韩一区二区视频免费看| 少妇裸体淫交视频免费看高清| 人体艺术视频欧美日本| 日韩一区二区三区影片| 久久久久久久久大av| 精品久久久久久久久久久久久| 日韩大片免费观看网站 | 亚洲国产欧洲综合997久久,| 国产亚洲精品av在线| 一区二区三区免费毛片| 男人舔奶头视频| 日本三级黄在线观看| 天堂√8在线中文| 国产精品不卡视频一区二区| 午夜老司机福利剧场| 日韩一本色道免费dvd| 欧美另类亚洲清纯唯美| 99久久中文字幕三级久久日本| 亚洲精华国产精华液的使用体验| 成人无遮挡网站| 天天躁夜夜躁狠狠久久av| 亚洲人成网站在线播| 天堂√8在线中文| 少妇熟女aⅴ在线视频| 国产精品嫩草影院av在线观看| 亚洲欧美日韩无卡精品| 一级爰片在线观看| 最近中文字幕高清免费大全6| 国产一级毛片七仙女欲春2| 干丝袜人妻中文字幕| 国产精品一及| 亚洲性久久影院| 日本一二三区视频观看| 亚洲欧美日韩无卡精品| 欧美一级a爱片免费观看看| 校园人妻丝袜中文字幕| 人妻夜夜爽99麻豆av| 欧美成人免费av一区二区三区| АⅤ资源中文在线天堂| 99久久九九国产精品国产免费| 能在线免费观看的黄片| 久久这里只有精品中国| 免费av观看视频| 亚洲精品影视一区二区三区av| 日韩av在线免费看完整版不卡| 99久久人妻综合| 激情 狠狠 欧美| 插逼视频在线观看| 女人被狂操c到高潮| 亚洲国产欧美在线一区| www.av在线官网国产| 久99久视频精品免费| 搡女人真爽免费视频火全软件| 少妇的逼水好多| 成人无遮挡网站| 久久99蜜桃精品久久| 亚洲在线自拍视频| 69人妻影院| 欧美一区二区精品小视频在线| 亚洲av不卡在线观看| 一边亲一边摸免费视频| 久久久久久久久久久免费av| 秋霞在线观看毛片| 国产精品国产三级国产av玫瑰| 国产免费男女视频| 成人av在线播放网站| 婷婷六月久久综合丁香| 中文资源天堂在线| 在现免费观看毛片| 久久精品久久久久久噜噜老黄 | 毛片一级片免费看久久久久| 久久久午夜欧美精品| 国产在视频线在精品| 我要搜黄色片| 一区二区三区高清视频在线| 日本三级黄在线观看| 青春草亚洲视频在线观看| 国产熟女欧美一区二区| 欧美3d第一页| 变态另类丝袜制服| 老司机影院毛片| 美女高潮的动态| 99在线人妻在线中文字幕| 亚洲经典国产精华液单| 久久99热6这里只有精品| 亚洲国产欧洲综合997久久,| 久久久久久久久久久免费av| 国产视频内射| АⅤ资源中文在线天堂| 老司机影院毛片| 女人被狂操c到高潮| 国产黄a三级三级三级人| av福利片在线观看| 超碰97精品在线观看| 久久人人爽人人爽人人片va| 波野结衣二区三区在线| 美女高潮的动态| 床上黄色一级片| 欧美日韩综合久久久久久| 国产精品伦人一区二区| 少妇的逼好多水| 欧美日韩综合久久久久久| 蜜桃亚洲精品一区二区三区| 久久久久久久午夜电影| 免费观看在线日韩| 亚洲人成网站在线播| 99九九线精品视频在线观看视频| 1000部很黄的大片| 尤物成人国产欧美一区二区三区| 搡老妇女老女人老熟妇| 亚洲激情五月婷婷啪啪| 亚洲人成网站高清观看| 国产精品久久久久久av不卡| 黄片无遮挡物在线观看| 久久这里有精品视频免费| 男人的好看免费观看在线视频| 看非洲黑人一级黄片| 午夜免费男女啪啪视频观看| 男插女下体视频免费在线播放|