城市面源SO2允许排放总量的网格分配方法

为了弥补A—P值法的不足，本文以A值法和SSIM模式为基础。提出了一种合理分配城市面源SO2允许排放总量的方法。该方法将A值法计算出的面源SO2允许排放总量作为一个约束条件放入SSIM（Source Strength Inversed Model）模式，调整该模式中各组源强的分担率，直至模式反演出的城市控制区各网格内的面源允许排放量之和与约束条件相等为止，此时模式输出的网格面源SO2允许排放量就是城市控制区内面源SO2允许排放总量的分配结果。应用该方法对石家庄市面源SO2允许排放总量进行网格分配，并采用美国工业源长期浓度复合模式（ISC3LT）对石家庄市SO2总量控制效果进行检验。结果表明：采用该方法分配面源SO2允许排放总量，可保证城市总量控制区内各网格上的SO2浓度值均满足国家标准。     

Footprint estimation for scalar flux measurements in the atmospheric surface layer

The flux footprint is the contribution, per unit emission, of each element of a surface area source to the vertical scalar flux measured at height z _m ; it is equal to the vertical flux from a unit surface point source. The dependence of the flux footprint on crosswind location is shown to be identical to the crosswind concentration distribution for a unit surface point source; an analytic dispersion model is used to estimate the crosswind-integrated flux footprint. Based on the analytic dispersion model, a normalized crosswind-integrated footprint is proposed that principally depends on the single variable z/z _m , where z is a measure of vertical dispersion from a surface source. The explicit dependence of the crosswind-integrated flux footprint on downwind distance, thermal stability and surface roughness is contained in the dependence of z on these variables. By also calculating the flux footprint with a Lagrangian stochastic dispersion model, it is shown that the normalized flux footprint is insensitive to the analytic model assumption of a self-similar vertical concentration profile.The National Center for Atmospheric Research is funded by the National Science Foundation.     

基于敏感源分析的动态大气污染排放方案模拟

选取2014年11月6—11日一次典型污染过程,以北京城区为重点关注的目标区域,基于印痕分析技术判别对该区域影响较大的敏感源区,设计重点区域减排试验和敏感源区逐日动态减排试验,利用区域化学传输模式WRF-Chem进行模拟对比。结果显示:上述两种方案对源强较高的减排当地PM_2.5浓度降低均有明显的改善作用,且在传输作用下会辐射影响到下游地区;但仅就目标区域而言,敏感源区减排方案的减排效率要远远高于重点区域减排方案。为了验证基于敏感源分析的动态减排方法的适用性,进一步开展了不同季节以及不同背景浓度的个例模拟。结果表明:基于敏感源区分析结果制定动态逐日减排措施,可降低削减成本、提高减排效率,以达到最具经济环境效益的减排效果。     

Temporal variation and source identification of black carbon at Lin’an and Longfengshan regional background stations in China

Black carbon (BC) is a component of fine particulate matter (PM_2.5), associated with climate, weather, air quality, and people’s health. However, studies on temporal variation of atmospheric BC concentration at background stations in China and its source area identification are lacking. In this paper, we use 2-yr BC observations from two background stations, Lin’an (LAN) and Longfengshan (LFS), to perform the investigation. The results show that the mean diurnal variation of BC has two significant peaks at LAN while different characteristics are found in the BC variation at LFS, which are probably caused by the difference in emission source contributions. Seasonal variation of monthly BC shows double peaks at LAN but a single peak at LFS. The annual mean concentrations of BC at LAN and LFS decrease by 1.63 and 0.26 μg m^–3 from 2009 to 2010, respectively. The annual background concentration of BC at LAN is twice higher than that at LFS. The major source of the LAN BC is industrial emission while the source of the LFS BC is residential emission. Based on transport climatology on a 7-day timescale, LAN and LFS stations are sensitive to surface emissions respectively in belt or approximately circular area, which are dominated by summer monsoon or colder land air flows in Northwest China. In addition, we statistically analyze the BC source regions by using BC observation and FLEXible PARTicle dispersion model (FLEXPART) simulation. In summer, the source regions of BC are distributed in the northwest and south of LAN and the southwest of LFS. Low BC concentration is closely related to air mass from the sea. In winter, the source regions of BC are concentrated in the west and south of LAN and the northeast of the threshold area of s_tot at LFS. The cold air mass in the northwest plays an important role in the purification of atmospheric BC. On a yearly scale, sources of BC are approximately from five provinces in the northwest/southeast of LAN and the west of LFS. These findings are helpful in reducing BC emission and controlling air pollution.     

北京城市大气环境污染机理与调控原理

该文主要介绍了科技部国家重点基础研究发展规划项目“首都北京及周边地区大气、水、土环境污染机理及调控原理”大气分项目的研究成果。项目分别于2001年和2003年重点开展了BECAPEX科学试验 (Beijing City Air Pollution Experiment)。BECAPEX试验同步进行城市边界层气象与大气化学观测, 通过卫星遥感、地面观测, 即城市空间和地面以及点与面结合的技术途径, 以揭示北京城市污染“空气穹隆”大气化学结构特征及其变化规律, 为城市环境大气动力-化学模式提供基本科学参数, 给出城市边界层大气物理化学过程综合模型, 为提高城市环境大气物理-化学过程耦合模式的准确性和可靠性提供科学依据。该项目揭示了北京城区和城近郊区城市边界层结构与湍流特征, 城市大气污染垂直结构特征; 发现了城市大气污染空间结构多尺度特征, 其中包括大气污染源影响和城市热岛多尺度特征; 揭示了城市大气重污染过程周边源影响域, 以及北京及周边地区气溶胶影响域和区域气候响应; 提出了北京市典型污染源排放清单; 发展了城市气象模式系统, 包括冠层模式、街谷环流和热力结构以及城市高大建筑群周围风环境数值模拟; 发展了空气质量模式技术, 包括二次气溶胶模拟试验、北京地区SO₂污染的长期模拟及不同类型排放源影响的计算与评估、影响北京地区的沙尘暴输送模拟、区域化学输送模式中NO_x和O₃源示踪法, 城市尺度的大气污染CAPPS模式及统计模型的应用、大气污染及紫外辐射数值预报模式和CMAQ-MOS空气质量预报方法; 改进了美国公共多尺度空气质量预报模式, 建立了CMAQ-MOS区域空气质量动力-统计模型预报模式, 以及发展的源同化技术, 突破了当前空气质量模式技术“瓶颈”, 使模式预报准确率明显提高。     

Transfer Matrices From the AES‐LRT Model

Abstract

The Canadian Atmospheric Environment Service (A ES) Long‐ Range Transport (LRT) model has been used in the Canada‐United States Memorandum of Intent programme to compute transfer matrices in order to quantify the source‐receptor relationships between emission regions and selected receptor sites. Four‐day backward trajectories were computed from the selected sites for the year 1978 and were started from the 925‐mb level(～600 m). The Lagrangian concentration / deposition model computed sulphur concentrations and depositions for 9 receptor sites using an emissions inventory divided into 15 Canadian and25 United States emission regions.

The 40×9 source‐receptor matrices show that the greatest impact on a receptor site usually results from an emission region close to the site although the regions giving the greatest impact for air concentration and wet deposition are not necessarily identical. In addition, the matrices show the impacts of all the emission regions on all of the receptor sites ranked by the magnitudes of the matrix elements. The per cent contribution from each emission region at each receptor site is shown as well as the overall per cent contributions from both Canada and the United States.

These matrices are an attempt to quantify source‐receptor relationships in Canada and the United States for assessing emission control strategies. The uncertainties associated with the matrices are being studied.     

Seasonal persistence and propagation of intraseasonal patterns over the Indian monsoon region

The space-time evolution of convection over the monsoon region containing the Indian subcontinent, the Indian Ocean and the West Pacific has been studied. A multi-channel singular spectrum analysis of the daily outgoing longwave radiation has yielded two intraseasonal oscillatory patterns and two large-scale standing patterns as the most dominant modes of intraseasonal variability. The oscillatory modes vary on time scales of about 45 and 28 days and their average cycles of variability are shown to correspond to the life cycles of active and break periods of monsoon rainfall over India. During an active (break) cycle, a convection (dry) anomaly zone first appears in the equatorial Indian Ocean, subsequently expands to cover the Indian subcontinent and finally contracts to disappear in the northern part of India. Some eastward and northward movements are found to be associated with both oscillatory modes, while westward movement may also be associated with the 28-day mode. The oscillatory modes are shown to have a large spatial scale extending to the West Pacific. One of the standing modes has anomalies of uniform sign covering the entire region and is related to El Niño and southern oscillation (ENSO) pattern. The other standing mode has a dipole structure in the equatorial Indian Ocean associated with large-scale anomalies over India with the same sign as those over the western part of the dipole. These two standing modes persist throughout the monsoon season, each maintaining its respective pattern. The seasonal mean monsoon is mainly determined by the two standing patterns, without much contribution from the oscillatory modes. The relative role of the standing patterns (ENSO mode and dipole mode) seems to be important in determining the seasonal mean during certain years.     

Nonlinear Dynamics in Atmospheric Chemistry Rate Equations

Nonlinear terms in the gas-phase rate equations of atmospheric trace constituents give rise to unexpected oscillations in the concentrations. For a simplified chemical scheme of the troposphere that contains only the generic reaction paths, the underlying dynamical structure is investigated. It is shown that the source strengths of CO and NO are the controlling parameters of the system. A linear stability analysis reveals that the steady state solutions lose stability due to the occurrence of two Hopf bifurcations. Furthermore, it appears that the dynamical behaviour of the oscillatory solutions is dominated by O₃ and CO. Based on the two-variable system (CO–O₃) it is shown that the oscillatory solution involves an autocatalytic ozone production phase which is followed by a phase in which CO is oxidised quickly. A simple expression is presented from which the period of the oscillation can be obtained. The implications for the present troposphere are unclear, since other hydrocarbons are present and transport is taking place. Nevertheless, the system nicely shows the general nonlinear mechanisms that operate in the tropospheric chemistry equations.     

RMAPS_Chem V1.0系统SO2排放清单优化效果评估

RMAPS_Chem V1.0系统是基于WRF_Chem模式建立的服务于华北区域雾霾等污染预报业务的模式系统,该研究着重针对系统中污染排放清单不确定性带来的SO₂浓度预报偏差较大问题,采用EnKF源反演和误差统计订正相结合的方法对排放清单进行了改进,形成了一套优化后的华北区域SO₂排放清单。通过输入初始清单和优化清单对2017年10月进行模拟,并与华北地区616个地面环境监测站观测值进行对比,结果表明:EnKF源反演结合误差统计订正的排放清单优化方法适用于SO₂排放清单的改进,有效降低了清单系统性偏差,针对主要区域及重点城市的检验显示模拟结果接近观测值;排放清单优化后模拟误差显著降低,如河北南部、山东西部至北京一带模式预报均方根误差与归一化平均绝对误差明显下降,区域内站点模拟误差呈正态分布特征,误差分布范围、最大概率出现范围均明显变窄,且最大误差概率明显上升。     

Turbulent Dispersion of Non-uniformly Emitted Passive Tracers in the Convective Boundary Layer

The impact of spatially non-uniform emissions on the turbulence dispersion of passive tracers in the convective boundary layer is studied by means of large-eddy simulation. We explicitly calculated the different terms of the budget equations for the concentrations, fluxes and variances, and used sub-domain averaging where each sub-domain is the typical size of a large-scale model grid cell. We found that the concentration profiles in the sub-domain where the emission takes place are lightly affected by the size of the emission release. This effect becomes more relevant in the downwind sub-domain. Although sub-domain averaged fluxes are not affected by the emission source size, concentration variances are dramatically increased when the emission shrinks. This increase originates from the mixing of highly concentrated air parcels with those of low concentrations. We also found that the concentration variance at the surface is driven neither by the position of the emission source nor the strength of the shear forcing but solely by the emission variance.     

基于奇异谱分析的南方涛动指数短期气候预测试验

The Southern Oscillation Index (SOI) time series is analyzed by means of the singular spectrum analysis (SSA) method with 60-month window length. Two major oscillatory pairs are found in the series whose pe riods are quasi-four and quasi-two years respectively. The auto-regressive model, which is developed on the basis of the Maximum Entropy Spectrum Analy sis, is fitted to each of the 9 leading components including the oscillatory pairs. The prediction of SOI with the 36-month lead is obtained from the reconstruction of these extrapolated series. Correlation coefficient between predicted series and 5 months running mean of observed series is up to 0.8. The model can successfully predict the peak and duration of the strong ENSO event from 1997 to 1998. It＇s also shown that the proper choice of reconstructed components is the key to improve the model prediction.     

新疆不同降水观测资料的比较及其差异的可能原因

利用新疆地区89个气象站和38个水文站1961~2005年的降水资料,对比分析了这两类降水观测资料在表征新疆地区降水的时间演变特征上所存在的差异,进而探究了造成这种差异的原因,并在此基础上研究了增暖背景下新疆地区降水与地形(包括海拔高度、坡度和坡向)的关系。研究发现:1)从区域平均来看,水文站观测的年降水量(221.4 mm)明显高于气象站的观测值(152.1 mm),其差值为69.3 mm,但是二者具有相同的时空变化特征;2)气象站与水文站降水均值之间的差异,主要受海拔高度的影响,坡度和坡向的影响次之。新疆地区地形地貌复杂,降水时空变化较大,观测数据分布不均匀均可导致对降水的估算存在较大差异,因此多源降水数据相互融合是客观估算该地区降水量及其变化的一个途径。     

Vorticity waves in a shallow basin

Observed oscillatory current patterns in the southern basin of Lake Michigan, with a distinctive peak in the energy spectrum at a period of about 90 h, are simulated using a linear potential vorticity conservation model. Solutions of the forced vorticity equation in a paraboloidal basin show rotational, oscillatory motions tuned to the low-frequency topographic modes that are very similar to the observed flow patterns. Topography-controlled vorticity waves are excited most effectively by wind episodes with frequency nearly in resonance with the topographic modes. Bottom resistance has no significant effect on the frequency equation; it simply decays the waves slowly in the open lake and more quickly near the coast. Flow patterns of both the gravest free vorticity wave and the corresponding forced wave consist of two opposite circulation cells separated by a null streamline through the center of the basin and rotating cyclonically near the free wave and atmospheric forcing frequencies, respectively. Interactions between the forced and free waves result in an apparent rotational pattern with a frequency the median of the two. A combination of elliptic—paraboloidal basin and shorter period forced modes can approximate the observed Lake Michigan response. Doppler shift, due to the persistence of cyclonic vorticity in the flow field, is also determined to be a factor in shifting the elliptical basin mode to a higher frequency.     

INTEX-B NO2排放源在中国地区的应用与检验

将INTEX-B排放源应用到空气质量模型Model3-CMAQ中,对中国地区对流层NO₂的浓度分布进行了数值模拟,并与OMI卫星对流层NO₂资料进行了对比和验证。结果表明：将INTEX-B排放源应用到Model3-CMAQ模式,模拟的NO₂浓度在中国地区的分布、季节变化规律与卫星资料所得结果一致。敏感性试验表明,工业及电厂排放对NO₂的浓度贡献最大,而交通排放的贡献相对较小,两种排放均主要集中在京津、长江三角洲等经济发达地区。     

INTEX-B NO_2排放源在中国地区的应用与检验

将INTEX-B排放源应用到空气质量模型Model3-CMAQ中,对中国地区对流层NO2的浓度分布进行了数值模拟,并与OMI卫星对流层NO2资料进行了对比和验证。结果表明：将INTEX-B排放源应用到Model3-CMAQ模式,模拟的NO2浓度在中国地区的分布、季节变化规律与卫星资料所得结果一致。敏感性试验表明,工业及电厂排放对NO2的浓度贡献最大,而交通排放的贡献相对较小,两种排放均主要集中在京津、长江三角洲等经济发达地区。     

亚热带森林非甲烷碳氢化合物的研究——Ⅰ.季节变化(钢瓶采样)

对中国亚热带森林地区的非甲烷碳氢化合物（ＮＭＨＣ）进行了一年多的钢瓶采样，并利用气相色谱仪对空气样品进行了分析，得到了大量ＮＭＨＣ的实测资料以及它们的季节变化规律和相互之间的关系。研究表明，ＮＭＨＣ的主要来源是当地植物的排放。对异戊二烯及其影响因子进行了详细研究。     

中国大陆黑碳气溶胶排放清单

The detailed high-resolution emission inventory of black carbon (BC) from China in the year 2000 was calculated. The latest fuel consumption data including fossil and biomass fuels, and socio-economic statistics used were obtained from government agencies, mostly at the county level, and some new emission factors (Efs) from local measurements were also used. National and regional summaries of emissions were presented at 0.2°×0.2°resolution. Total BC emission was 1499.4 Gg in 2000, mainly due to the burning of coal and biomass. The BC emission estimated here is higher than those in previous studies, mainly because the emissions from coal burning by rural industries and residences were previously underestimated. More BC aerosols were emitted from eastern China than western China. A strong seasonal dependence was observed for BC emissions, with peaks in January and December and low emissions in July and August; and this seasonality is mainly due to patterns in residential heating and the open burning of crop residues.     

The Interannual Variation of Transboundary Contributions from Chinese Emissions of PM_(2.5) to South Korea

In recent years,several studies pointed out that anthropogenic emission sources in China which significantly contribute to the PM_2.5mass burden was an important cause of particulate pollution in South Korea.However,most studies generally focused upon a single pollution event.It is rare to see comprehensive research that captures those features prone to interannual variations concerning the transboundary pollutant contribution in South Korea using a unified method.In this paper,we establish the emission inventories covering East Asia in 2010,2015,and 2017,and then conduct the source apportionment by applying a coupled regional air quality model called the Integrated Source Apportionment Module(ISAM).Comparison of simulated and observed PM_2.5mass concentration at 165 CNEMC(China National Environmental Monitoring Center)sites suggests that the PM_2.5concentrations are well represented by the modeling system.The model is used to quantitatively investigate the contribution from emission sources in China to PM_2.5concentrations over South Korea and those features found to be prone to interannual variations are then discussed.The results show that the average annual contribution of PM_2.5has dropped significantly from 28.0%in 2010 to 15.7%in 2017,which strongly suggests that China has achieved remarkable results in the treatment of atmospheric particulates.     

北京城区二氧化碳时空分布及湍流谱特征

利用北京325 m气象塔上安装的7层CO₂涡动相关系统在2014年12月到2015年11月的观测资料,分析了北京城区不同高度上CO₂浓度、通量时空分布及湍流谱的特征。结果表明:城市CO₂浓度日变化除了冬季都呈现双峰型,冬季由于人为碳源排放的大幅增加,双峰型不明显。每层的CO₂浓度、通量都有明显的季节变化:冬季最高,春末、夏季最低。CO₂浓度整体随高度的增加而降低。北京城区是CO₂源,CO₂通量的日变化不如CO₂浓度日变化规律明显。CO₂通量在47 m以下为负,47 m以上为正。通量在140 m以下随高度的增加而增加;140m以上随高度的增加而减少。根据对CO₂时空分布的分析可知:边界层CO₂浓度、通量强烈受到碳源、下垫面植被、大气稳定度、环境温度和天气过程等因素的影响。各变量谱与Kaimal等的研究结果接近:归一化速度谱和CO₂谱在惯性子区有－2/3的斜率,在低频区与稳定度参数（Z/L）有一定的关系。这说明复杂地形的城市下垫面的湍流谱结构与平坦地形相比没有太大的实质性差异。     

改进的互功率谱相位法在雷声声源定位中的应用

互相关函数法是一种常用的时延估计方法,但在低信噪比和闪电存在多分叉的情况下,该方法的效果不理想。针对这一问题,该文选取广州地区2009年8月24日19:08:04和2012年5月5日08:01:00的两次闪电过程的雷声记录,结合高速摄像资料,比较了互相关函数法和改进的互功率谱相位法在雷声声源定位中的应用效果。首先,利用互相关函数法和改进的互功率谱相位法计算声源点的方向信息,定位结果与高速摄像照片对比具有较好的一致性,证明了定位方法的可靠性,对比结果表明:在低信噪比环境和多分叉闪电情况下,改进的互功率谱相位法能够获得比互相关函数法更好的定位效果。最后给出了综合利用改进的互功率谱相位法计算得到的方向信息和声源点与观测点的距离信息得到的雷声声源三维定位结果。