摘要
污染源同化反演方法是利用观测数据对污染源进行"自上而下"约束更新,可为大气污染源清单提供一种快速校验和更新的方法,还可为"自下而上"的污染源清单补充更多时空分布信息。综述了质量平衡法、集合卡尔曼滤波法和四维变分法等同化反演方法的基本原理及在污染源清单反演优化中的应用方法;介绍了观测数据、初始污染源清单的获取和其不确定性的计算方法以及污染源与观测数据之间关系的获取方法;总结了污染源同化反演方法优化在不同尺度污染源清单反演中的应用,以及在国内外的应用情况。采用集合卡尔曼滤波法或四维变分法并利用伴随模式或去耦合直接法确定污染源与观测数据的关系,可在一定程度上避免模式中的非线性问题,但提高污染源同化反演效果仍需在观测数据质量和模式模拟性能方面进一步改进。
Data assimilation is to use the observed data to update the pollution sources from top to bottom, which could provide a fast way to update and verify the emission inventories, and also add further information about the temporal variation and some sources that could not easily be quantified by traditional bottom-up methods. The basic principles of several assimilation inversion methods, including the mass balance,ensemble Calman filter and fourdimensional variation methods were reviewed, and their application methods in emission inventory inversion and optimization introduced. The acquisition method of observation data and the primary emission inventory, the calculation method of uncertainty, and the acquisition method of the relationship between pollution sources and observation data were also introduced. The application of pollution source assimilation inversion methods at different scales of pollution source inventory inversion, both at home and abroad, were summarized. The ensemble Kalman filter or the 4D variation method combined with adjoint model or decoupling direct method was used to determine the relationship between pollution sources and observation data, able to avoid the nonlinearities in the model to a certain extent. The more accurate observations and more precise air quality model were needed in the future to improve the effect of pollution sources inversion.
引文
[1]STREETS D G,BOND T C,CARMICHAEL G R,et al.An inventory of gaseous and primary aerosol emissions in Asia in the year 2000[J].Journal of Geophysical Research Atmospheres,2003,108(21):1-23.
[2]ZHANG Q,STREETS D G,CARMICHAEL G R,et al.Asian emissions in 2006 for the NASA INTEX-B mission[J].Atmospheric Chemistry&Physics Discussions,2009,9(14):5131-5153.
[3]雷宇.中国人为源颗粒物及关键化学组分的排放与控制研究[D].北京:清华大学,2008.
[4]BOCQUET M,ELBERN H,ESKES H,et al.Data assimilation in atmospheric chemistry models:current status and future prospects for coupled chemistry meteorology models[J].Atmospheric Chemistry&Physics,2015,15(10):5325-5358.
[5]STREETS D G,CANTY T,CARMICHAEL G R,et al.Emissions estimation from satellite retrievals:a review of current capability[J].Atmospheric Environment,2013,77(3):1011-1042.
[6]白晓平,李红,方栋,等.资料同化在空气质量预报中的应用[J].地球科学进展,2007,22(1):66-73.BAI X P,LI H,FANG D,et al.Applications of data assimilation in air quality prediction[J].Advances in Earth Science,2007,22(1):66-73.
[7]朱江,汪萍.集合卡尔曼平滑和集合卡尔曼滤波在污染源反演中的应用[J].大气科学,2006,30(5):871-882.ZHU J,WANG P.Ensemble Kalman smoother and ensemble Kalman filter approaches to the joint air quality state and emission estimation problem[J].Chinese Journal of Atmospheric Sciences,2006,30(5):871-882.
[8]徐祥德,丁国安.城市大气污染源追踪“广义判识”与Nudging源同化集成技术途径[J].中国工程科学,2012,14(10):10-22.XU X D,DING G A.Generalized method of urban air pollutant source tracking and integration approach of Nudging source assimilation[J].Engineering Science,2012,14(10):10-22.
[9]TANG W,COHAN D S,LAMSAL L N,et al.Inverse modeling of Texas NOxemissions using space-based and ground-based NO2observations[J].Atmospheric Chemistry&Physics,2013,13(21):11005-11018.
[10]TANG X,ZHU J,WANG Z F,et al.Inversion of CO emissions over Beijing and its surrounding areas with ensemble Kalman filter[J].Atmospheric Environment,2013,81(4):676-686.
[11]朱国富.数值天气预报中分析同化基本方法的历史发展脉络和评述[J].气象,2015,41(8):986-996.ZHU G F.Remarks on development of basic methods of atmospheric data assimilation for numerical weather rediction[J].Meteorological Monthly,2015,41(8):986-996.
[12]RODGERS C D.Inverse methods for atmospheric sounding[M]//Inverse methods for atmospheric sounding.Maryland:University of Maryland,2000.
[13]朱国富.理解大气资料同化的根本概念[J].气象,2015,41(4):456-463.ZHU G F.Understanding of the fundamental concept of atmospheric data assimilation[J].Meteorological Monthly,2015,41(4):456-463.
[14]MARTIN R V,JACOB D J,CHANCE K,et al.Global inventory of nitrogen oxide emissions constrained by space-based observations of NO2,columns[J].Journal of Geophysical Research Atmospheres,2003,108(Suppl 17):4537.
[15]BOERSMA K F,JACOB D J,BUCSLA E J,et al.Validation of OMI tropospheric NO2,observations during INTEX-B and application to constrain NOxmath container loading mathjax,emissions over the Eastern United States and Mexico[J].Atmospheric Environment,2008,42(19):4480-4497.
[16]PALMER P I,JACOB D J,FIORE A M,et al.Mapping isoprene emissions over North America using formaldehyde column observations from space[J].Physics&Atmospheric Science,2003,108(6):315-323.
[17]GHUDE S D,KULKARNI S H,JENA C,et al.Application of satellite observations for identifying regions of dominant sources of nitrogen oxides over the Indian subcontinent[J].Journal of Geophysical Research Atmospheres,2013,118(2):1075-1089.
[18]LIN J T,MCELROY M B,BOERSMA K F.Constraint of anthropogenic NOxemissions in China from different sectors:a new methodology using multiple satellite retrievals[J].Atmospheric Chemistry&Physics Discussions,2010,10(1):63-78.
[19]HOOGHIEMSTRA P B,KROL M C,MEIRINK J F,et al.Optimizing global CO emission estimates using a four-dimensional variational data assimilation system and surface network observations[J].Atmospheric Chemistry&Physics,2011,11(10):4705-4723.
[20]BERGAMASCHI P,FRANKENBERG C,MEIRINK J F,et al.Inverse modeling of global and regional CH4emissions using SCIAMACHY satellite retrievals[J].Journal of Geophysical Research Atmospheres,2009,114(Suppl 22):2191-2196.
[21]COURTIER P,THPAUT J N,HOLLINGSWORTH A.A strategy for operational implementation of 4D-var,using an incremental approach[J].Quarterly Journal of the Royal Meteorological Society,1994,120:1367-1387.
[22]ELBERN H,STRUNK A,SCHMIDT H,et al.Emission rate and chemical state estimation by 4-dimensional variational inversion[J].Atmospheric Chemistry&Physics,2007,7(1):3749-3769.
[23]MIYAZAKI K,ESKES H J,SUDO K.Global NOxemission estimates derived from an assimilation of OMI tropospheric NO2columns[J].Atmospheric Chemistry&Physics,2012,11(5):31523-31583.
[24]HOUTEKAMER P L,MITCHELL H L.A sequential ensemble Kalman filter for atmospheric data assimilation[J].Monthly Weather Review,2001,129(1):123-137.
[25]EVENSEN G.Sequential data assimilation with a nonlinear quasigeostrophic model using Monte Carlo methods to forecast error statistics[J].Journal of Geophysical Research Oceans,1994,99(5):10143-10162.
[26]唐晓,朱江,王自发,等.基于集合卡尔曼滤波的区域臭氧资料同化试验[J].环境科学学报,2013,33(3):796-805.TANG X,ZHU J,WANG Z F,et al.Regional ozone data assimilation experiment based on ensemble Kalman filter[J].Acta Scientiae Circumstantiae,2013,33(3):796-805.
[27]MIJLING B,van DER A R J.Using daily satellite observations to estimate emissions of short-lived air pollutants on a mesoscopic scale[J].Journal of Geophysical Research Atmospheres,2012,117(17):1-20.
[28]吴钲,谢旻,高阳华,等.利用集合均方根卡尔曼滤波反演重庆地区SO2源排放[J].环境科学研究,2018,31(1):25-33.WU Z,XIE M,GAO Y H,et al.Inversion of SO2emissions over Chongqing with ensemble square root Kalman filter[J].Research of Environmental Scineces,2018,31(1):25-33.
[29]XU X,XIE L,CHENG X,et al.Application of an adaptive nudging scheme in air quality forecasting in China[J].Journal of Applied Meteorology and Climatology,2008,47(8):2105-2114.
[30]程兴宏,徐祥德,安兴琴,等.2013年1月华北地区重霾污染过程SO2和NOx的CMAQ源同化模拟研究[J].环境科学学报,2016,36(2):638-648.CHENG X H,XU X D,AN X Q,et al.Inverse modeling of SO2and NOxemissions using an adaptive nudging scheme implemented in CMAQ model in North China during heavy haze episodes in January 2013[J].Acta Scientiae Circumstantiae,2016,36(2):638-648.
[31]李明君,陈东升,程水源,等.遗传算法用于CMAQ模式污染源清单优化的研究[J].北京工业大学学报,2011(12):1862-1868.LI M J,CHEN D S,CHENG S Y,et al.Optimizing emission inventory of CMAQ by means of genetic algorithms[J].Journal of Beijing University of Technology,2011(12):1862-1868.
[32]冯帆,王自发,唐晓.一个基于打靶法的大气污染源反演自适应算法[J].大气科学,2016,40(4):719-729.FENG F,WANG Z F,TANG X.Development of an adaptive algorithm based on the shooting method and its application in the problem of estimating air pollutant emissions[J].Chinese Journal of Atmospheric Sciences,2016,40(4):719-729.
[33]XU W Y,ZHAO C S,RAN L,et al.A new approach to estimate pollutant emissions based on trajectory modeling and its application in the North China Plain[J].Atmospheric Environment,2011(11):75-83.
[34]ESKES H J,BOERSMA K F.Averaging kernels for DOAS totalcolumn satellite retrievals[J].Atmospheric Chemistry&Physics,2003,3(5):895-910.
[35]BOERSMA K F,ESKES H J,BRINKSMA E J.Error analysis for tropospheric NO2retrieval from space[J].Journal of Geophysical Research Atmospheres,2004,109(4):D04311-1-D04311-20.
[36]US EPA.Technical assistance document for the chemiluminescence measurement of nitrogen dioxide:EPA-600/4-75-003[R].Washington DC:Environmental Monitoring and Support Laboratory,US EPA,1975.
[37]WILLIAMS E J,FEHSENFELD F C,JOBSON B T,et al.Comparison of ultraviolet absorbance,chemiluminescence,and DOAS instruments for ambient ozone monitoring[J].Environmental Science&Technology,2006,40(18):5755-5762.
[38]HITZENBERGER R,BERNER A,GALAMBOS Z,et al.Intercomparison of methods to measure the mass concentration of the atmospheric aerosol during INTERCOMP2000:influence of instrumentation and size cuts[J].Atmospheric Environment,2004,38:6467-6476.
[39]KONOVALOV I B,BEEKMANN M,RICHTER A,et al.Inverse modelling of the spatial distribution of NOxemissions on a continental scale using satellite data[J].Atmospheric Chemistry&Physics Discussions,2005,5(6):1747-1770.
[40]GU D,WANG Y,SMELTZER C,et al.Anthropogenic emissions of NOxover China:reconciling the difference of inverse modeling results using GOME-2 and OMI measurements[J].Journal of Geophysical Research Atmospheres,2014,119(12):7732-7740.
[41]OLIVIER J G J,BERDOWSKI J J M.Global emission sources and sinks[J].Computer Aided Optimum Design of StructuresⅧ,2001,2698:622.
[42]National Academy of Science.EMEP-CORINAIR emission inventory guidebook[M].San Francisco:ARRB Group Limited,2007.
[43]HEALD C L,JACOB D J,JONES D B A,et al.Comparative inverse analysis of satellite(MOPITT)and aircraft(TRACE-P)observations to estimate Asian sources of carbon monoxide[J].Journal of Geophysical Research Atmospheres,2004,109(23):D23306-1-D23306-17.
[44]GU D,WANG Y,YIN R,et al.Inverse modelling of NOx emissions over eastern China:uncertainties due to chemical nonlinearity[J].Atmospheric Measurement Techniques,2016,9(10):5193-5201.
[45]GEDDES J A,MARTIN R V.Global deposition of total reactive nitrogen oxides from 1996 to 2014 constrained with satellite observations of NO2columns[J].Atmospheric Chemistry&Physics,2017,17(16):1-44.
[46]HUNEEUS N,CHEVALLIER F,BOUCHER O.Estimating aerosol emissions by assimilating observed aerosol optical depth in a global aerosol model[J].Atmospheric Chemistry&Physics,2012,12(10):4585-4606.
[47]MIYAZAKI K,ESKES H J,SUDO K,et al.Simultaneous assimilation of satellite NO2,O3,CO,and HNO3data for the analysis of tropospheric chemical composition and emissions[J].Atmospheric Chemistry&Physics&Discussions,2012,264:1017-1023.
[48]RUSSELL A R,VALIN L C,BUCSELA E J,et al.Space-based constraints on spatial and temporal patterns of NOxemissions in California,2005-2008[J].Environmental Science&Technology,2010,44(9):3608-3615.
[49]YOSHIDA Y,DUNCAN B N,RETSCHER C,et al.The impact of the 2005 Gulf hurricanes on pollution emissions as inferred from ozone monitoring instrument(OMI)nitrogen dioxide[J].Atmospheric Environment,2010,44(11):1443-1448.
[50]SOURI A H,CHON Y,JEON W,et al.Constraining NOx,emissions using satellite NO2,measurements during 2013DISCOVER-AQ Texas campaign[J].Atmospheric Environment,2016,131:371-381.
[51]苏芳,邵敏,蔡旭辉,等.利用逆向轨迹反演模式估算北京地区甲烷源强[J].环境科学学报,2002,22(5):586-591.SU F,SHAO M,CAI X H,et al.Estimation of methane emissions in Beijing area using backward trajectory inversion model[J].Acta Scientiae Circumstantiae,2002,22(5):586-591.
[52]孟凯,程兴宏,徐祥德,等.基于CMAQ源同化反演方法的京津冀局地污染源动态变化特征模拟研究[J].环境科学学报,2017,37(1):52-60.MENG K,CHENG X H,XU X D,et al.Spatial-temporal variations of pollutant emission sources inversed by adaptive nudging scheme over Beijing-Tianjin-Hebei region based on the CMAQ model[J].Acta Scientiae Circumstantiae,2017,37(1):52-60.
[53]KOURTIDIS K,GEORGOULIAS A K,MIJLING B,et al.A new method for deriving trace gas emission inventories from satellite observations:the case of SO2,over China[J].Science of the Total Environment,2018,612:923-930.
[54]LIN J,MCELROY M B.Detection from space of a reduction in anthropogenic emissions of nitrogen oxides during the Chinese economic downturn[J].Atmospheric Chemistry and Physics,2011,11(15):8171-8188.
[55]LI M,KLIMONT Z,ZHANG Q,et al.Comparison and evaluation of anthropogenic emissions of SO2and NOx,over China[J].Atmospheric Chemistry&Physics,2018,18(5):1-28.
[56]COOPER M,MARTIN R V,PADMANMBHAN A,et al.Comparing mass balance and adjoint methods for inverse modeling of nitrogen dioxide columns for global nitrogen oxide emissions[J].Journal of Geophysical Research,2017,122:4718-4734.
[57]JIANG Z,WORDEN J R,JONES D B A,et al.Constraints on Asian ozone using Aura TES,OMI and Terra MOPITT[J].Atmospheric Chemistry&Physics,2015,14(13):99-112.