边界层方案对华北低层O3垂直分布模拟的影响

利用WRF-Chem模式,采用3种边界层参数化方案 (YSU, MYJ和ACM2),针对1个晴空、静稳日 (2013年8月26日20:00—27日20:00(北京时)) 进行模拟,着重分析不同边界层参数化方案对夜间残留层形成及日出前后O₃浓度垂直分布形式的模拟效果,并与固城站地面及垂直同步观测资料进行对比。结果表明:3种边界层参数化方案均能够模拟出温度及风速的区域分布形式以及风温垂直结构的变化特征;相比之下,MYJ方案模拟的夜间边界层高度较YSU方案和ACM2方案明显偏高,该对比结果可能是导致近地面污染物浓度模拟差异的重要原因;在夜间稳定层结至日出后稳定状态打破的边界层结构演变过程中,采用YSU方案和ACM2方案模拟的温度和风速垂直扩线形式与观测结果更为接近;同样采用非局地闭合的YSU方案和同时考虑局地和非局地闭合的ACM2方案,对于边界层高度内O₃浓度垂直分布形式的模拟效果具有明显优势。     

Evaluation of the Weather Research and Forecast/Urban Model Over Greater Paris

Meteorological modelling in the planetary boundary layer (PBL) over Greater Paris is performed using the Weather Research and Forecast (WRF) numerical model. The simulated meteorological fields are evaluated by comparison with mean diurnal observational data or mean vertical profiles of temperature, wind speed, humidity and boundary-layer height from 6 to 27 May 2005. Different PBL schemes, which parametrize the atmospheric turbulence in the PBL using different turbulence closure schemes, may be used in the WRF model. The sensitivity of the results to four PBL schemes (two non-local closure schemes and two local closure schemes) is estimated. Uncertainties in the PBL schemes are compared to the influence of the urban canopy model (UCM) and the updated Coordination of Information on the Environment (CORINE) land-use data. Using the UCM and the CORINE land-use data produces more realistic modelled meteorological fields. The wind speed, which is overestimated in the simulations without the UCM, is improved below 1,000 m height. Furthermore, the modelled PBL heights during nighttime are strongly modified, with an increase that may be as high as 200 %. At night, the impact of changing the PBL scheme is lower than the impact of using the UCM and the CORINE land-use data.     

复杂地形上气象场对空气质量数值模拟 结果影响的研究

利用WRF模式三种边界层参数化方案(YSU、MYJ、ACM2)产生的气象场分别驱动多尺度空气质量模式CMAQ,对兰州市西固区冬季2005年1月27日至2月2日期间SO₂和NO₂浓度进行了数值模拟,将模拟结果与同期监测的污染物浓度进行对比分析,结果表明:WRF模式不同边界层参数化方案模拟输出的气象场驱动CMAQ模式所模拟的SO₂和NO₂浓度均可以反映出污染物的时空变化特征,CMAQ模式具有模拟复杂下垫面高分辨率污染物输送特征的能力;WRF模式的边界层参数化方案选为局地与非局地闭合方案(ACM2)时,模拟的气象场驱动CMAQ模式得到的空气污染物浓度分布特征最优,这主要是由于ACM2的湍流输送机制较为合理,模拟的边界层低层气象场更接近实际,从而可以较好地模拟污染物的输送特征;当CMAQ模式的垂直混合方案与WRF模式的湍流输送方案一致时(均采用ACM2方案),模式间的兼容性好.     

Sensitivity analysis of planetary boundary layer schemes using the WRF model in Northern Colombia during 2016 dry season

This study conducted meteorological simulations in northern Colombia by analyzing different planetary boundary layer (PBL) schemes available in the numerical Weather Research and Forecasting (WRF) model. The study area included three nested domains with horizontal resolutions of 18 km, 6 km, and 2 km, with 38 vertical levels. The evolution and structure of the PBL were analyzed during the driest months (March, April, and May 2016) and in regions with the highest particulate matter concentrations. Sensitivity analysis of the WRF model was performed with two local and two non-local PBL schemes. The results were validated using observations of the surface air temperature, relative humidity, and surface wind speed collected from three meteorological stations in the area. The PBL heights were experimentally determined using radiosonde data provided by a station located in the center of the study area. Variations in PBL heights were estimated using linear regression methods and minimization of statistical errors for the bulk Richardson number, as well as analysis of vertical temperature and wind profiles. The WRF model reliably reproduced the daily values and diurnal cycles of temperature, relative humidity, and wind speed within the PBL and accounted for the influence of topography and sea breezes. Horizontal heat advection dominates the upwelling of air masses when sea breezes are active. The onshore wind direction starts to change from east to northwest, implying a decay in the land breeze regime. All schemes overestimate the mixing height and tend to underestimate surface air temperature values at night. All show wetter conditions and underestimate wind speed. Although the non-local Yonsei University (YSU) scheme shows the best performance, it also shows the largest sources of errors when determining the behavior of the surface layer during stable conditions. Relative humidity and wind speed estimates provided by the local Mellor‐Yamada‐Nakanishi‐Niino (MYNN) scheme were closer to those recorded at the meteorological stations.     

Intercomparison of Planetary Boundary-Layer Parametrizations in the WRF Model for a Single Day from CASES-99

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes—the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes—the Mellor–Yamada–Janjić (MYJ), quasi-normal scale elimination (QNSE), and Bougeault–Lacarrére (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.     

Influence of Vertical Eddy Diffusivity Parameterization on Daily and Monthly Mean Concentrations of O3 and NOy

Two parameterization schemes for vertical eddy diffusivity were utilized to investigate their impacts on both the daily and monthly mean concentrations of ozone and NOy, which are the major fractions of the sum of all reactive nitrogen species, i.e., NOy=NO＋NO2＋HNO3＋PAN. Simulations indicate that great changes in the vertical diffusivity usually occur within the planetary boundary layer （PBL）. Daily and monthly mean concentrations of NOy are much more sensitive to changes in the vertical diffusivity than those of ozone and ozone and NOy levels only at or in （relatively） clean sites and areas, where long-range transport plays a crucial role, display roughly equivalent sensitivity. The results strongly suggest that a widely-accepted parameterization scheme be selected and the refinement of the model＇s vertical resolution in the PBL be required, even for regional and long-term studies, and ozone only being examined in an effort to judge the model＇s performance be unreliable, and NOy be included for model evaluations.     

A 2-D Non-local Closure Model for Atmospheric Boundary Layer Simulations

Ａ　２－Ｄ　Ｎｏｎ－ｌｏｃａｌ　Ｃｌｏｓｕｒｅ　Ｍｏｄｅｌ　ｆｏｒ　Ａｔｍｏｓｐｈｅｒｉｃ　Ｂｏｕｎｄａｒｙ　Ｌａｙｅｒ　ＳｉｍｕｌａｔｉｏｎｓＪｉａｎｇＷｅｉｍｅｉａｎｄＷａｎｇＸｕｅｍｅｉ（ＧｕｅｓｔＲｅｓｅａｒｃｈｅｒｉｎＬＡＰＣ，Ｂｅｉｊｉｎ...     

尺度自适应大气边界层参数化改进及其对一次海雾的数值模拟研究

大气边界层湍流运动是地球大气运动最重要的能量输送过程之一.当数值模式分辨率接近活跃含能湍涡长度尺度时,湍流运动被部分解析,被称为"灰色区域",传统的边界层方案不适合此时模式湍流问题的描述.为了提高模式边界层方案在包括"灰色区域"的不同网格尺度上的描述能力,适应不同分辨率模式的需要,在雷诺平均湍流理论基础上,修正Mell...     

AN AIR POLLUTION PREDICTION TECHNIQUE FOR URBAN DISTRICTS BASED ON MESO-SCALE NUMERICAL MODEL

1 INTRODUCTION Since the 1930’s, urban air pollution has caused increasing concerns across the world[1]. Various pollution prediction techniques have been put into use in the 1970’s and significant results have been achieved[2]. Air quality prediction has been operational over recent years in many parts of China owing to cooperation between meteorological and environmental protection departments[3 – 9]. The prediction is conducted via either numerical model, potential or stochastic app…     

WRF模式多种边界层参数化方案对四川盆地不同量级降水影响的数值试验

利用中尺度模式WRF三种边界层参数化方案(MYJ、YSU和ACM2),对2012年四川盆地夏季连续40天逐日降水量进行数值试验,并检验评估了不同边界层参数化方案下模式对分级降水量和边界层结构的模拟能力,分析了各参数化方案对降水量模拟差异的可能原因。结果表明:三种边界层参数化方案对较小量级(小雨和中雨)降水量的模拟,24 h时效优于48 h,ACM2方案效果较好;对较大量级(大雨和暴雨)降水的模拟,48 h时效优于24 h,YSU方案模拟效果较好。对比分析温江站加密探空观测与模式模拟的大气边界层结构表明,ACM2方案对小量级降水时边界层结构的模拟较为准确,而YSU方案更适合于温江站大量级降水时边界层结构的模拟。不同边界层参数化方案对各量级降水量模拟差异的可能原因是边界层湍流混合强度的不同,MYJ方案湍流混合作用较弱,导致底层大量水汽积聚,不稳定性强,容易产生虚假降水,因此对各量级降水模拟能力均有限;YSU方案具有强烈的垂直混合强度,有利于局地水汽的向上输送,更易达到大量级降水发生发展的条件,适用于盆地较大量级降水的模拟;ACM2方案在保证足够湍流混合强度的同时,在较稳定条件下会关闭非局地输送,不致于产生过强降水,适合盆地较小量级降水的数值模拟     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

春季黄海海雾WRF参数化方案敏感性研究

利用2005—2011年10次春季黄海海雾个例开展WRF模式参数化方案敏感性研究。结果表明:边界层方案对WRF模式雾区模拟结果起决定作用,而微物理方案影响较小,它主要影响海雾浓度和高度。边界层与微物理方案的最佳组合为YSU与Lin方案,最差为Mellor-Yamada与WSM5方案;Mellor-Yamada和QNSE方案模拟的近海面湍流过强,导致边界层过高,不利于海雾的发展与维持;而MYNN与YSU方案刻画的湍流强度与边界层高度合适,有利于海雾发展与维持。MYNN方案虽与YSU方案相当,但在大多数海雾个例中,后者明显优于前者,而在有些个例中却刚好相反。因此对于某一具体海雾个例而言,所用边界层方案仍需在它们之中选择最优者。这些信息可为黄海海雾WRF模式边界层与微物理方案的选择与改进提供参考。     

A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

北京交通环境污染诊断模拟初步研究

建立了能够用于评估大型城市交通环境污染与气象条件关系的数值模拟技术,即三维边界层光化学数值模式.其中,高分辨率的精细城市边界层气象场、细致的城市地理与环境参数、空气污染物输送扩散与大气化学模拟以及交通废气污染源的分布等诸多因子是实施模拟研究的几个重要环节.从模拟结果看出,各污染物在地面附近浓度相对高值区位于排放源附近和下风方.在边界层顶附近,污染物浓度分布较均匀,分布趋势和气流一致.虽然仅仅考虑了城区污染物本底浓度和主要交通干道机动车的排放,但交通排放对城区污染物浓度的贡献是较大的.     

Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes,Boundary Conditions and Spin-Up

We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about \(2\) K) and substantially underestimates it at night (about \(4\) K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( \(\approx \) 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.     

Numerical modeling of pollutant transport and chemistry during a high‐ozone event in northern Taiwan

An air pollution prediction model system (APOPS) is developed and applied to northern Taiwan with complex terrain and local thermal circulations. It consists of a nonhydrostatic mesoscale meteoro‐logical model system (MMPMS) and a gas/aerosol transport and air quality model (GATAM). The basic processes relevant to modeling the urban air pollution problems such as meteorology, dispersion, chemistry and deposition are solved at the same time on practically the same grid. The APOPS was tested on a high‐ozone event in northern Taiwan on 16 November 1998. Comparison with observed surface winds shows able to predict local flow patterns such as sea/land breezes and mountain‐valley wind in this high air pollution episode. The predicted surface concentrations of ozone and other pollutants are compared with measured values, and a fairly good agreement with the mean normalized biases of −6%, −11%, for one‐day simulation and for daytime, respectively, is obtained for ozone. Thus, it is confirmed that the APOPS can be utilized to predict urban air quality in complex terrain area.     

The role of precursor emissions on ground level ozone concentration during summer season in Poland

Three online coupled chemical transport model simulations were analyzed for three summer months of 2015 in Poland. One of them was run with default emission inventory, the other two with NO_x and VOC emissions reduced by 30%, respectively. Obtained ozone concentrations were evaluated with data from air quality measurement stations and ozone sensitivity to precursor emissions was estimated by ozone concentration differences between simulations and with the use of indicator ratios. They were calculated based on modeled mixing ratios of ozone, total reactive nitrogen and its components, nitric acid and hydrogen peroxide. The results show that the model overestimates ozone concentrations with the largest errors in the morning and evening, which is primarily related to the way vertical mixing is resolved by the model. Better model performance for ozone is achieved in rural than urban environment, as PBL and mixing mechanisms play more significant role in urban areas. Modeled ozone shows mixed sensitivity to precursor concentrations, similarly to other European regions, but indicator ratios have different values than are found in literature, particularly H₂O₂/HNO₃ is larger than in southern Europe. However, indicator ratios often differ between locations and transition values need to be established individually for a given region.     

气溶胶辐射效应对边界层结构及夹卷特征影响的观测分析

2017年12月22日至2018年1月18日利用无人机携带温、湿和颗粒物浓度探测仪对南京地区灰霾污染条件下大气边界层垂直结构开展加密观测。通过比较不同灰霾污染条件下温度、湿度和PM_2.5（直径小于2.5微米的颗粒物）浓度的垂直结构差异,结合地面热通量、2米空气温度、相对湿度、风速、风向及主要大气污染物（如臭氧和PM_2.5）浓度,定量评估了气溶胶辐射效应对边界层和夹卷过程的影响。分析表明,灰霾或气溶胶削弱到达地表太阳辐射,减小地表感热通量,延迟边界层发展,增加近地层大气稳定度,降低边界层高度,并加重灰霾污染。灰霾污染物在混合层顶处累积,导致PM_2.5浓度最大变化出现在边界层顶部而不是近地层。气溶胶辐射效应对夹卷特征及其特征参数有重要影响。灰霾浓度升高时,夹卷区厚度增加;无量纲化夹卷速度随对流理查逊数的变化不再符合负1次方幂函数关系,与大涡模拟结果一致。本研究进一步指出,为提高重霾污染条件下天气和空气质量数值预报水平,必须考虑气溶胶辐射效应对边界层和夹卷参数化的影响。     

基于旋翼无人机的大气边界层环境气象垂直观测及订正方法的研究

基于旋翼无人机开展大气边界层观测可为气象要素和大气污染物垂直结构的研究提供具有高时空解析能力的新方法,有助于深入理解低层大气物理化学变化机制。本文详述了旋翼无人机在开展大气边界层环境气象垂直观测实验的应用及优势。基于自主研发的旋翼无人机环境气象观测平台,通过开展传感器在无人机上不同的搭载位置,以及旋翼无人机与探空气球、高塔的对比观测实验,明确了旋翼无人机对气象环境观测的影响及合理的搭载方式。研究进一步在湖北重污染天气条件下开展了0—1000 m的大气边界层垂直观测,并研发了基于旋翼无人机姿态数据的大气边界层气象要素及污染物垂直观测的订正方法。结果表明:实验获取了2—10 m垂直分辨率的高质量大气廓线数据,可精细捕捉大气边界层及其逆温层高度和污染物浓度等要素的垂直变化特征。本文旨在为无人机观测的科研应用提供一种技术可行且数据可靠的观测手段。     

Investigation of Turbulence Parametrization Schemes with Reference to the Atmospheric Boundary Layer Over the Aegean Sea During Etesian Winds

The spatial structure of the marine atmospheric boundary layer (MABL) over the Aegean Sea is investigated using the Weather Research and Forecasting (WRF) mesoscale model. Two ‘first-order’ non-local and five ‘1.5-order’ local planetary boundary-layer (PBL) parametrization schemes are used. The predictions from the WRF model are evaluated against airborne observations obtained by the UK Facility for Airborne Atmospheric Measurements BAe-14 research aircraft during the Aegean-GAME field campaign. Statistical analysis shows good agreement between measurements and simulations especially at low altitude. Despite the differences between the predicted and measured wind speeds, they reach an agreement index of 0.76. The simulated wind-speed fields close to the surface differ substantially among the schemes (maximum values range from 13 to \(18\hbox { m s}^{-1}\) at 150-m height), but the differences become marginal at higher levels. In contrast, all schemes show similar spatial variation patterns in potential temperature fields. A warmer (1–2 K) and drier (2–3\(\hbox { g kg}^{-1})\) layer than is observed, is predicted by almost all schemes under stable conditions (eastern Aegean Sea), whereas a cooler (up to 2 K) and moister (1–2\(\hbox { g kg}^{-1})\) layer is simulated under near-neutral to nearly unstable conditions (western Aegean Sea). Almost all schemes reproduce the vertical structure of the PBL and the shallow MABL (up to 300 m) well, including the low-level jet in the eastern Aegean Sea, with non-local schemes being closer to observations. The simulated PBL depths diverge (up to 500 m) due to the different criteria applied by the schemes for their calculation. Under stable conditions, the observed MABL depth corresponds to the height above the sea surface where the simulated eddy viscosity reaches a minimum; under neutral to slightly unstable conditions this is close to the top of the simulated entrainment layer. The observed sensible heat fluxes vary from ?40 to \(25\hbox { W m}^{-2}\), while the simulated fluxes range from ?40 to \(40\hbox { W m}^{-2}\); however, all of the schemes’ predictions are close to the observations under unstable conditions. Finally, all schemes overestimate the friction velocity, although the simulated range (from 0.2 to \(0.5\hbox { m s}^{-1})\) is narrower than that observed (from 0.1 to \(0.7\hbox { m s}^{-1})\).