DB32T-大氣二氧化碳和甲烷本底濃度篩分技術(shù)規(guī)范報(bào)批稿

Q/LB.□XXXXX-XXXXDB32/TXXXX—XXXX目次TOC\o"1-1"\h前言 II引言 III1范圍 12規(guī)范性引用文件 13術(shù)語和定義 14資料收集和處理 25本底篩分方法 56篩分方法選擇 67篩分流程 7附錄A（資料性）本底篩分流程 8附錄B（資料性）觀測(cè)站點(diǎn)信息 9參考文獻(xiàn) 10前言本文件按照GB/T1.1—2020《標(biāo)準(zhǔn)化工作導(dǎo)則第1部分：標(biāo)準(zhǔn)化文件的結(jié)構(gòu)和起草規(guī)則》的規(guī)定起草。請(qǐng)注意本文件的某些內(nèi)容可能涉及專利。本文件的發(fā)布機(jī)構(gòu)不承擔(dān)識(shí)別專利的責(zé)任。本文件由江蘇省氣象局提出。本文件由江蘇省氣象標(biāo)準(zhǔn)化技術(shù)委員會(huì)歸口。本文件起草單位：江蘇省氣候中心、江蘇省氣象探測(cè)中心。本文件主要起草人：陳燕、惠品宏、徐家平、徐萌、沈瑱、李崇志、韋芬芬、周晨、喬賀、徐進(jìn)。引言實(shí)現(xiàn)碳達(dá)峰碳中和是一場(chǎng)廣泛而深刻的經(jīng)濟(jì)社會(huì)系統(tǒng)性變革，掌握溫室氣體本底濃度及其變化情況，有助于客觀、準(zhǔn)確地測(cè)算區(qū)域溫室氣體源匯總量及變化，是實(shí)現(xiàn)碳達(dá)峰碳中和的重要基礎(chǔ)。為規(guī)范我省區(qū)域大氣二氧化碳和甲烷本底濃度篩分，特制定本標(biāo)準(zhǔn)。大氣二氧化碳和甲烷本底濃度篩分技術(shù)規(guī)范范圍本文件規(guī)定了大氣二氧化碳、甲烷的本底濃度篩分所需的資料、篩分方法、篩分方案選擇、篩分流程的要求。本文件適用于氣象、生態(tài)、環(huán)境等行業(yè)在受人為或自然源影響的區(qū)域開展二氧化碳、甲烷的本底濃度篩分和數(shù)據(jù)處理,其他行業(yè)可參考本文件。規(guī)范性引用文件下列文件中的內(nèi)容通過文中的規(guī)范性引用而構(gòu)成本文件必不可少的條款。其中，注日期的引用文件，僅該日期對(duì)應(yīng)的版本適用于本文件；不注日期的引用文件，其最新版本（包括所有的修改單）適用于本文件。GB/T3840制定地方大氣污染物排放標(biāo)準(zhǔn)的技術(shù)方法GB/T35221地面氣象觀測(cè)規(guī)范總則GB/T35227地面氣象觀測(cè)規(guī)范風(fēng)向和風(fēng)速Q(mào)X/T118地面氣象觀測(cè)資料質(zhì)量控制QX/T174大氣成分站選址要求術(shù)語和定義下列術(shù)語和定義適用于本文件。

溫室氣體greenhousegas大氣中能夠吸收紅外輻射的氣體成分，主要包括二氧化碳（CO2）、甲烷（CH4）、氧化亞氮（N2O）、水汽（H2O）、氫氟碳化物（HFCs）、全氟碳化物（PFCs）、六氟化硫（SF6）和臭氧（O3）等。[來源：QX/T125-2011，3.1，有修改]

二氧化碳carbondioxide分子式為CO2，在大氣中的滯留時(shí)間可達(dá)幾十年或上百年，是影響地球輻射平衡的主要溫室氣體。源主要是化石燃料和生物質(zhì)的燃燒、土地利用變化以及工業(yè)過程排放，匯主要是陸地生態(tài)系統(tǒng)和海洋的吸收。[來源：QX/T125-2011，4.2，有修改]

甲烷methane分子式為CH4，在大氣中的滯留時(shí)間約為12年。源主要是農(nóng)業(yè)（如牲畜和稻田）、垃圾處理過程、化石燃料生產(chǎn)和使用（煤炭、天然氣/石油開采）、生物燃料和生物質(zhì)燃燒的排放、濕地和其他自然過程等，匯主要是羥基自由基氧化反應(yīng)。[來源：QX/T125-2011，4.3，有修改]

濃度concentration特定氣體摩爾數(shù)與同體積干空氣內(nèi)全部氣體分子摩爾數(shù)之比。

本底濃度backgroundconcentration遠(yuǎn)離局地排放源、濾除了局地自然環(huán)境和人為活動(dòng)直接影響、基本混合均勻的大氣組成成分濃度。

局地濃度localconcentration受局地排放源和吸收匯影響的非本底大氣組成成分濃度。

梯度理查森數(shù)gradientRichardsonnumber用兩個(gè)高度上的平均風(fēng)速和平均溫度代替理查森數(shù)表達(dá)式中風(fēng)速垂直切變項(xiàng)和溫度垂直梯度，來確定大氣穩(wěn)定度。

穩(wěn)健局部回歸robustlocalregression在一定時(shí)間窗口內(nèi)對(duì)觀測(cè)值進(jìn)行逐步逼近回歸擬合的統(tǒng)計(jì)學(xué)穩(wěn)健估計(jì)方法。資料收集和處理資料種類溫室氣體觀測(cè)站點(diǎn)信息：溫室氣體觀測(cè)站點(diǎn)的經(jīng)緯度、海拔高度、進(jìn)氣口高度、下墊面類型、周圍環(huán)境（表B.1）。溫室氣體觀測(cè)數(shù)據(jù)：溫室氣體觀測(cè)站點(diǎn)建站以來的CO2、CH4濃度觀測(cè)數(shù)據(jù)，CO2濃度單位為10-6mol/mol，CH4濃度單位為10-9mol/mol，下同。氣象觀測(cè)站點(diǎn)信息：氣象觀測(cè)站點(diǎn)的經(jīng)緯度、海拔高度、風(fēng)觀測(cè)高度、下墊面類型、周圍環(huán)境（表B.2）。氣象觀測(cè)數(shù)據(jù)：溫室氣體觀測(cè)站點(diǎn)同址、同期的風(fēng)速、風(fēng)向、氣溫、氣壓、云量資料。如果無同址觀測(cè)，需選擇距離近、觀測(cè)環(huán)境一致性好的氣象站點(diǎn)同期觀測(cè)數(shù)據(jù)。資料收集溫室氣體觀測(cè)站點(diǎn)信息：來源于國家測(cè)繪地理信息主管部門、地方測(cè)繪部門、溫室氣體觀測(cè)部門的地面和遙感觀測(cè)，觀測(cè)站點(diǎn)選址符合QX/T174的規(guī)定。溫室氣體濃度觀測(cè)數(shù)據(jù)：來自溫室氣體觀測(cè)部門，包括建站以來經(jīng)過質(zhì)量控制后的逐小時(shí)的CO2、CH4濃度數(shù)據(jù)，符合時(shí)間長度可以是月、季、年、年際尺度。氣象觀測(cè)站點(diǎn)信息：來源于國家測(cè)繪地理信息主管部門、地方測(cè)繪部門、氣象部門。氣象觀測(cè)數(shù)據(jù)：來自氣象部門，包括經(jīng)過質(zhì)量控制后的地面及大氣邊界層高度以內(nèi)的逐小時(shí)風(fēng)速、風(fēng)向、氣溫、氣壓數(shù)據(jù)以及逐小時(shí)云量數(shù)據(jù)，符合GB/T35221、GB/T35227及QX/T118的規(guī)定。數(shù)據(jù)處理溫室氣體濃度觀測(cè)數(shù)據(jù)數(shù)據(jù)完整率用于本底篩分的CO2和CH4濃度觀測(cè)數(shù)據(jù)的時(shí)間分辨率為逐小時(shí)，總數(shù)據(jù)及逐月數(shù)據(jù)完整率均需達(dá)到95%以上。計(jì)算數(shù)據(jù)完整率時(shí)不統(tǒng)計(jì)由于儀器故障、停機(jī)、系統(tǒng)維護(hù)等造成的缺測(cè)。逐月逐風(fēng)向溫室氣體平均濃度基于觀測(cè)時(shí)段內(nèi)CO2或CH4濃度觀測(cè)數(shù)據(jù)以及風(fēng)向數(shù)據(jù)，逐月逐風(fēng)向計(jì)算CO2或CH4平均濃度，十六風(fēng)向所對(duì)應(yīng)的方位、符號(hào)、中心角度、角度范圍以及靜風(fēng)標(biāo)準(zhǔn)見表1，計(jì)算方法見式（1）： Cm,n=i=1Nm,nC式中：Cm,n——特定月份特定風(fēng)向CO2或CH4平均濃度，m表示1~12月，n表示1~16Cm,ni——特定月份特定風(fēng)向逐時(shí)CO2或Nm,n——特定月份特定風(fēng)向參與統(tǒng)計(jì)的CO2或CH4觀測(cè)數(shù)風(fēng)向符號(hào)與度數(shù)對(duì)照表序號(hào)方位符號(hào)中心角度/°角度范圍/°1北N0348.76-11.252北東北NNE22.511.26-33.753東北NE4533.76-56.254東東北ENE67.556.26-78.755東E9078.76-101.256東東南ESE112.5101.26-123.757東南SE135123.76-146.258南東南SSE157.5146.26-168.759南S180168.76-191.2510南西南SSW202.5191.26-213.7511西南SW225213.76-236.2512西西南WSW247.5236.26-258.7513西W270258.76-281.2514西西北WNW292.5281.26-303.7515西北NW315303.76-326.2516北西北NNW337.5326.26-348.75注：風(fēng)速小于或等于0.2m/s按靜風(fēng)處理，不計(jì)入十六風(fēng)向。逐月溫室氣體十六風(fēng)向平均濃度基于4.3.1.2計(jì)算的逐月逐風(fēng)向CO2或CH4平均濃度，逐月計(jì)算十六風(fēng)向平均濃度，計(jì)算方法見式（2）： Cm,1?16=n=116式中：Cm,1?16——特定月份CO2或CH4十六風(fēng)向平均濃度逐月溫室氣體十六風(fēng)向濃度標(biāo)準(zhǔn)差基于4.3.1.2計(jì)算的逐月逐風(fēng)向CO2或CH4平均濃度，逐月計(jì)算十六風(fēng)向濃度標(biāo)準(zhǔn)差，計(jì)算方法見式（3）： σm,1?16=n=116Cm,n式中：σm,1?16——特定月份CO2或CH4十六風(fēng)向濃度標(biāo)準(zhǔn)差逐月溫室氣體平均濃度基于觀測(cè)時(shí)段內(nèi)CO2或CH4濃度觀測(cè)數(shù)據(jù)，逐月計(jì)算CO2或CH4平均濃度，計(jì)算方法見式（4）： Cm=i=1NmC式中：Cm——特定月份參與統(tǒng)計(jì)的CO2或CH4平均濃度Cmi——特定月份參與統(tǒng)計(jì)的逐時(shí)CO2或CHNm——特定月份參與統(tǒng)計(jì)的CO2或CH4觀測(cè)數(shù)逐月溫室氣體濃度標(biāo)準(zhǔn)差基于觀測(cè)時(shí)段內(nèi)CO2或CH4濃度觀測(cè)數(shù)據(jù)，逐月計(jì)算CO2或CH4濃度標(biāo)準(zhǔn)差，計(jì)算方法見式（5）： σm=i=1NmC式中：σm——特定月份參與統(tǒng)計(jì)的CO2或CH4濃度標(biāo)準(zhǔn)差大氣穩(wěn)定度梯度理查森數(shù)利用地面和大氣邊界層以內(nèi)氣溫、氣壓、風(fēng)速、風(fēng)向計(jì)算各時(shí)刻梯度理查森數(shù)Ri，并根據(jù)表2對(duì)大氣穩(wěn)定度進(jìn)行等級(jí)劃分，A、B、C、D、E和F分別表示強(qiáng)不穩(wěn)定、不穩(wěn)定、弱不穩(wěn)定、中性、較穩(wěn)定和穩(wěn)定。梯度理查森數(shù)計(jì)算方法見式（6）和式（7）： Ri=g?Z1?Z θ=T?(P0P式中：g——重力加速度，取9.8m/s2；Z1——低層Z2——高層θ1——低層θ2——高層V1——低層V2——高層θa——低層、高層位溫θ——位溫（K）；T——?dú)鉁兀↘）；P——?dú)鈮海╤Pa）；P0——標(biāo)準(zhǔn)氣壓，取1000hPakd——常數(shù)，取0.286基于梯度理查森數(shù)的大氣穩(wěn)定度等級(jí)劃分穩(wěn)定度等級(jí)Ri范圍ARi<-3.433B-3.433≤Ri<-1.375C-1.375≤Ri<-0.233D-0.233≤Ri<0.145E0.145≤Ri<0.636F0.636≤Ri帕斯奎爾-特納爾穩(wěn)定度按照GB/T3840中的方法，利用逐時(shí)云量、風(fēng)速數(shù)據(jù)，采用帕斯奎爾-特納爾穩(wěn)定度分類法來劃分大氣穩(wěn)定度等級(jí)，首先計(jì)算出每個(gè)時(shí)刻的太陽高度角，結(jié)合云量數(shù)據(jù)確定太陽輻射等級(jí)；隨后結(jié)合太陽輻射等級(jí)與地面風(fēng)速確定大氣穩(wěn)定度等級(jí)，可分為強(qiáng)不穩(wěn)定、不穩(wěn)定、弱不穩(wěn)定、中性、較穩(wěn)定和穩(wěn)定六級(jí)，分別由A、B、C、D、E和F表示。本底篩分方法穩(wěn)健局部回歸篩分法基于穩(wěn)健局部回歸算法對(duì)CO2或CH4有效觀測(cè)數(shù)據(jù)進(jìn)行本底篩分，定義觀測(cè)的濃度值為式（8），假設(shè)Cobs是連續(xù)的，可采用局部線性回歸模型對(duì)Cobs進(jìn)行估算，因此在任何給定時(shí)刻t0附近足夠小的鄰域內(nèi)，可以認(rèn)為Cobs(ti)為一條直線，在t0附近的數(shù)可以應(yīng)用最小二乘法求得?？紤]到大氣CO2和CH4濃度的季節(jié)變化將帶寬設(shè)置為60天（1440小時(shí)） Cobsti=Cbgt式中：Cobsti——CO2或CH4Cbgti——CO2或CH4C?tiEi綜合篩分法穩(wěn)健局部回歸篩分采用5.1穩(wěn)健局部回歸篩分法對(duì)CO2或CH4有效觀測(cè)數(shù)據(jù)進(jìn)行本底篩分。大氣穩(wěn)定度判定梯度理查森數(shù)判定經(jīng)過穩(wěn)健局部回歸篩分后，對(duì)剩余CO2或CH4濃度觀測(cè)值進(jìn)行大氣穩(wěn)定度判定，利用CO2或CH4觀測(cè)同期氣象數(shù)據(jù)，根據(jù)4.3.2.1計(jì)算梯度理查森數(shù)，并進(jìn)行大氣穩(wěn)定度分級(jí)，若穩(wěn)定度等級(jí)為D、E、F，表示高低層混合不充分，受局地排放源影響較大，則將相應(yīng)時(shí)刻的CO2或CH4觀測(cè)數(shù)據(jù)篩分為局地值。帕斯奎爾-特納爾穩(wěn)定度判定若氣象數(shù)據(jù)無法滿足梯度理查森數(shù)計(jì)算，根據(jù)4.3.2.2計(jì)算帕斯奎爾-特納爾穩(wěn)定度并進(jìn)行分級(jí)，若穩(wěn)定度等級(jí)為D、E、F，則將相應(yīng)時(shí)刻的CO2或CH4觀測(cè)數(shù)據(jù)篩分為局地值?；旌蠒r(shí)段判定若氣象數(shù)據(jù)無法滿足帕斯奎爾-特納爾穩(wěn)定度以及梯度理查森數(shù)計(jì)算，則將每日18:00～次日8:00（北京時(shí)間）的CO2或CH4觀測(cè)數(shù)據(jù)篩分為局地值。風(fēng)向與濃度判定經(jīng)過穩(wěn)健局部回歸篩分及大氣穩(wěn)定度判定后，對(duì)剩余CO2或CH4濃度觀測(cè)值進(jìn)一步基于風(fēng)向與濃度進(jìn)行判定。結(jié)合CO2或CH4觀測(cè)同期逐時(shí)風(fēng)向數(shù)據(jù)，根據(jù)4.3.1.2計(jì)算CO2或CH4的逐月逐風(fēng)向平均濃度Cm,n，根據(jù)4.3.1.3計(jì)算逐月十六風(fēng)向平均濃度Cm,1?16，根據(jù)4.3.1.4計(jì)算逐月十六風(fēng)向濃度標(biāo)準(zhǔn)差σm,1?16，對(duì)于m月n風(fēng)向，若Cm,n?Cm,1?16>2σ風(fēng)速判定經(jīng)過穩(wěn)健局部回歸篩分、大氣穩(wěn)定度判定及風(fēng)向與濃度判定后，對(duì)剩余CO2或CH4濃度觀測(cè)值進(jìn)一步基于風(fēng)速閾值進(jìn)行判定。結(jié)合CO2或CH4觀測(cè)同期逐時(shí)風(fēng)速數(shù)據(jù)，將風(fēng)速小于1.5m/s（一級(jí)風(fēng)）的時(shí)刻所對(duì)應(yīng)的CO2或CH4觀測(cè)濃度值篩分為局地值。異常值判定經(jīng)過穩(wěn)健局部回歸篩分、大氣穩(wěn)定度判定、風(fēng)向與濃度判定及風(fēng)速判定后，針對(duì)剩余CO2或CH4濃度觀測(cè)值，根據(jù)4.3.1.5逐月計(jì)算平均值Cm，根據(jù)4.3.1.6逐月計(jì)算標(biāo)準(zhǔn)差σm，對(duì)于某時(shí)刻，若Cmi?Cm>3σm，則將該時(shí)刻數(shù)據(jù)篩分方法選擇優(yōu)先采用綜合篩分法，若無氣象數(shù)據(jù)，則采用穩(wěn)健局部回歸篩分法。采用綜合篩分法時(shí)，基于所能收集到的氣象數(shù)據(jù)種類選擇大氣穩(wěn)定度判定方法，具體如下：a)有CO2或CH4濃度觀測(cè)數(shù)據(jù)同期地面風(fēng)速和風(fēng)向資料，大氣穩(wěn)定度判定采用混合時(shí)段判定；b)有CO2或CH4濃度觀測(cè)數(shù)據(jù)同期地面風(fēng)速、風(fēng)向、云量資料，大氣穩(wěn)定度判定采用帕斯奎爾-特納爾穩(wěn)定度判定；c)有CO2或CH4濃度觀測(cè)數(shù)據(jù)同期地面風(fēng)速、風(fēng)向、氣溫、氣壓資料，同期大氣邊界層以內(nèi)風(fēng)速、氣溫、氣壓資料，大氣穩(wěn)定度判定采用梯度理查森數(shù)判定。篩分流程按照附錄A中圖A.1所示的流程進(jìn)行資料收集、數(shù)據(jù)處理、本底篩分、篩分選擇。

（資料性）

本底篩分流程本底篩分流程圖

（資料性）

觀測(cè)站點(diǎn)信息溫室氣體觀測(cè)站點(diǎn)信息表站名地址經(jīng)度（度）緯度（度）海拔高度（米）填寫日期觀測(cè)站下墊面類型進(jìn)氣口高度（米）觀測(cè)場(chǎng)周圍環(huán)境觀測(cè)站點(diǎn)十六方位照片北北東北東北東東北東東東南東南南東南南南西南西南西西南西西西北西北北西北氣象觀測(cè)站點(diǎn)信息表站名地址經(jīng)度（度）緯度（度）海拔高度（米）填寫日期觀測(cè)站下墊面類型風(fēng)觀測(cè)高度（米）觀測(cè)場(chǎng)周圍環(huán)境觀測(cè)站點(diǎn)十六方位照片北北東北東北東東北東東東南東南南東南南南西南西南西西南西西西北西北北西北參考文獻(xiàn)[1]QX/T125-2011溫室氣體本底濃度觀測(cè)術(shù)語[2]20thWMO/IAEAMeetingonCarbonDioxide.OtherGreenhouseGasesandRelatedMeasurementTechniques(GGMT-2019).[3]IrwinJS.Estimatingplumedispersion-arecommendedgeneralizedscheme.FourthSymp.Turbulence.Diffusion,AirPollut.Reno,NevadaAmMetSoc,1979.62-69.[4]O'DohertyS,SimmondsP,CunnoldD,etal.InsituchloroformmeasurementsatAdvancedGlobalAtmosphericGasesExperimentatmosphericresearchstationsfrom1994to1998.JournalofGeophysicalResearch:Atmospheres,2001,106(D17):20429-20444.[5]PuJJ,XuHH,HeJ,etal.EstimationofregionalbackgroundconcentrationofCO2atLin'anStationinYangtzeRiverDelta,China.AtmosphericEnvironment,2014,94:402-408.[6]RuckstuhlAF,JacobsonMP,FieldRW,etal.Baselinesubtractionusingrobustlocalregressionestimation.JournalofQuantitativeSpectroscopy&RadiativeTransfer,2001,68:179-193.[7]RuckstuhlAF,HenneS,ReimannS,etal.Robustextractionofbaselinesignalofatmospherictracespeciesusinglocalregression.AtmosphericMeasurementTechniques,2012,5:2613–2624.[8]ThoningKW,TansPP,KomhyrWD.AtmosphericcarbondioxideatMaunaL