The retrieval of the Asian dust depolarization ratio in Korea with the correction of the polarization-dependent transmission

The linear particle depolarization ratios were retrieved from the observation with a multiwavelength Raman lidar at the Gwangju Institute of Science and Technology (GIST), Korea (35.11°N, 126.54°E). The measurements were carried out in spring (March to May) 2011. The transmission ratio measurements were performed to solve problems of the depolarization-dependent transmission at a receiver of the lidar and applied to correct the retrieved depolarization ratio of Asian dust at first time in Korea. The analyzed data from the GIST multiwavelength Raman lidar were classified into three categories according to the linear particle depolarization ratios, which are pure Asian dust on 21 March, the intermediate case which means Asian dust mixed with urban pollution on 13 May, and haze case on 10 April. The measured transmission ratios were applied to these cases respectively. We found that the transmission ratio is needed to be used to retrieve the accurate depolarization ratio of Asian dust and also would be useful to distinguish the mixed dust particles between intermediate case and haze. The particle depolarization ratios of pure Asian dust were approximately 0.25 at 532 nm and 0.14 at 532 nm for the intermediate case. The linear particle depolarization ratios of pure Asian dust observed with the GIST multiwavelength Raman lidar were compared to the linear particle depolarization ratios of Saharan dust observed in Morocco and Asian dust observed both in Japan and China.     

GRAPES混合同化方法在青藏高原区域的初步试验

针对青藏高原地区气象观测站点稀少和模式同化分析质量较低的问题,将GRAPES区域集合变分(En-3DVAR)混合同化系统应用于青藏高原地区,进行了单点理想试验和真实观测资料同化分析预报试验,分析评估青藏高原混合同化分析增量及预报误差的水平垂直结构特征及其合理性,并与中国东部平原地区进行比对。单点理想试验表明,En-3DVAR混合同化系统中背景误差协方差具备流依赖属性。真实资料混合同化试验结果表明,基于集合预报估计的分析增量具有流依赖的特征,特别是在高原地区和槽脊系统附近;青藏高原地区分析场的绝对误差总体低于3DVAR系统,其平均绝对误差在中层和高层低于平原地区,说明在青藏高原地区的改进效果略优于平原地区。需要关注的是,青藏高原地区En-3DVAR混合同化分析增量总体大于3DVAR,特别是近地面层u风分量分析增量明显偏大,这可能与青藏高原复杂地形有关。     

Parallel Computing of a Variational Data Assimilation Model for GPS／MET Observation Using the Ray-Tracing Method

The Spectral Statistical Interpolation (SSI) analysis system of NCEP is used to assimilate me-teorological data from the Global Positioning Satellite System (GPS/MET) refraction angles with thevariational technique. Verified by radiosonde, including GPS/ME% observations into the analysis makesan overall improvement to the analysis variables of temperature, winds, and water vapor. However, thevariational model with the ray-tracing method is quite expensive for numerical weather prediction andclimate research. For example, about 4 000 GPS/MET refraction a;~gles need to be assimilated to producean ideal global analysis. Just one iteration of minimization will take more than 24 hours CPU time onthe NCEP‘s Cray C90 computer. Although efforts have been taken to reduce the computational cost, it isstill prohibitive for operational data assimilation. In this paper, a parallel version of the three-dimensional variational data assimilation model of GPS/MET occultation measurement suitable for massive parallelprocessors architectures is developed. The divide-and-conquer strategy is used to achieve parallelism and isimplemented by message passing. The authors present the principles for the code‘s design and examine theperformance on the state-of-the-art parallel computers in China. The results show that this parallel modelscales favorably as the number of processors is increased. With the Memory-IO technique implemented bythe author, the wall clock time per iteration used for assimilating 1420 refraction angles is reduced from45 s to 12 s using 1420 processors. This suggests that the new parallelized code has the potential to beuseful in numerical weather prediction (NWP) and climate studies.     

Delaunay三角剖分法在降水量插值中的应用

Delaunay三角剖分方法在空间分析中具有重要地位,文中简要介绍了Delaunay三角网特性和常用的3类算法,并对随机增长法实现过程进行了详细阐述。根据三角分片线性插值原理,求得插值系数,实现对任意点的三角分片线性插值。利用2008年中国2200个观测站的08时24 h降水量资料,对全中国范围及划分的8个区域内相应的0.28125°×0.28125°降水量格点场,使用交叉检验方法,对比分析了三角分片线性插值和反距离权重法的估值准确率。结果表明：在各区域,三角分片线性插值法的均方根误差偏小;在站点较密集的区域,均方根误差、平均绝对误差比较中,三角分片线性插值都有一定的优势;在平均误差对比中,三角分片线性插值优势明显,在全中国范围交叉检验中,三角分片线性插值法对应的年平均误差是0.005 mm,而反距离权重法为-0.107 mm,对其可能的原因进行了分析,证明了Delaunay三角剖分法的合理性。同时,从图形上展示了降水量的Delaunay三角网的三维结构图和三角分片线性插值后的格点场, 在直观上,Delaunay三角剖分后得到降水分布和实况保持一致,并有较好的视觉效果;通过三角分片线性插值得到的格点场降水量分布图,克服了反距离权重法的固有缺陷,使获得的降水量格点场趋于合理,提高了插值精度。最后,探讨了Delaunay三角网在气象领域的应用前景。     

Multi-species chemical data assimilation with the Danish Eulerian hemispheric model: system description and verification

Satellite retrievals of atmospheric composition provide a wealth of data on a global scale. These complement results from atmospheric chemistry-transport models (CTMs), and can be combined using data assimilation. We present two assimilation schemes coupled to the Danish Eulerian Hemispheric Model (DEHM), a three-dimensional, off-line CTM with full photochemistry: a variant on the ensemble Kalman filter and the three-dimensional variational scheme. The aim of this paper is to describe the two schemes and present an initial assessment of their impacts on model skill. Retrievals of multiple atmospheric trace gases are assimilated, namely: NO₂ tropospheric column densities, CH₄ total column densities, and partial column concentrations of O₃, CO and CH₄; these data are retrieved from four satellite sensors. Data for each species are assimilated independently of one another, and other species are only adjusted indirectly via the model’s chemistry and dynamics. Assimilation results are compared with measurements from surface monitoring stations and other satellite retrievals, and preliminary validation results are presented.Reference simulations (without assimilation) grossly underestimate surface CO concentrations, and both assimilation schemes eliminate this large and systematic model bias. The assimilation improves the spatial correlation of modelled CO with surface observations, and improves the spatial correlation between forecasts and retrievals for CO, NO₂ and O₃. Results for CH₄ show a loss of skill due to a mismatch in model bias between two assimilated CH₄ data-sets. Finally, we discuss differences in methodology and results between this paper and a recent study on multi-species chemical data assimilation. Joint optimisation of initial conditions and emission rates offers a promising direction for improving modelled boundary-layer concentrations.     

中尺度区域模式层顶高度对高空风场数值模拟的影响分析

本文基于中尺度区域模式WRF,开展模式层顶高度变化对高空气象要素,特别是高空风场数值模拟影响的研究。通过设计模式顶高45、5 hPa两个试验,同化来源于NOAA-15、NOAA-18、NOAA-19和METOP-2的AMSU-A辐射计高通道数据,表明提高模式层顶能够使卫星更多的高通道样本数量进入同化系统,达到减小背景场误差,同时减小高于层顶通道辐射能量对低层通道影响的目的,一定程度上改进了同化效果,从而改善高空气象要素,特别是风场的模拟效果,与观测值的均方根误差减小了约0.4~0.5 m·s^-1。     

充分利用数据库资料提高雷雨预报准确率

2000年7月17日午后及晚上双流机场产生了两次强雨天气过程，机场预报员成功地对上述过程做出准确以预报；本文就当日的天气过程进行了较为详细地分析，进一步总结出民航气象数据库WAFS（世界区域预报系统的简称，以下同）产品500hPa高空风既流场形势预预告图对近两年来双流机场多起夏季雷雨天气过程成功预测所起的作用，为双流机场夏季雷雨的预报提供了一个着眼点。     

充分利用数据库资料提高雷雨预报准确率

2000年7月17日午后及晚上双流机场产生了两次强雷 雨天气过程，机场 预报员成功地对上述过程做出准确预报；本文就当日的天气过程进行了较为详细的分析，进 一步总结出民航气象数据库WAFS(世界区域预报系统的简称，以下同)产品500hPa高空风既流 场形势预告图对近两年来双流机场多起夏季雷雨天气过程成功预测所起的作用，为双流机场 夏季雷雨的预报提供了一个着眼点。     

Atmosphere data assimilation system for the Siberian region with the WRF-ARW model and three-dimensional variational analysis WRF 3D-Var

The system of the cyclic assimilation of data on atmospheric conditions used in the West Siberian Administration for Hydrometeorology and Environmental Monitoring is described. It is based on the WRF-ARW mesoscale atmospheric model and on the WRF 3D-Var system of the three-dimensional variational analysis of data. The system is verified when the first approximation data (6-hour forecast) and WRF-ARW forecasts with the lead time up to 24 hours are compared with the observational data. The problems of assimilation of observations from the AMSU-A and AIRS satellite instruments are considered. The effect of using AMSU-A and AIRS for the analysis in the Novosibirsk region is estimated. The experiments demonstrated that the cyclic data assimilation system operates successfully. The AMSU-A observations improve the quality of analyses and forecasts in winter. In summer the impact of satellite observations on the forecast skill scores is ambiguous. Good short-term forecasts are provided by the initial conditions obtained using the system of detailing of the NCEP large-scale analysis.     

Intercomparison of the performance of MM5/WRF with and without satellite data assimilation in short-range forecast applications over the Indian region

The present study is conducted to verify the short-range forecasts from mesoscale model version5 (MM5)/weather research and forecasting (WRF) model over the Indian region and to examine the impact of assimilation of quick scatterometer (QSCAT) near surface winds, spectral sensor microwave imager (SSM/I) wind speed and total precipitable water (TPW) on the forecasts by these models using their three-dimensional variational (3D-Var) data assimilation scheme for a 1-month period during July 2006. The control (without satellite data assimilation) as well as 3D-Var sensitivity experiments (with assimilating satellite data) using MM5/WRF were made for 48 h starting daily at 0000 UTC July 2006. The control run is analyzed for the intercomparison of MM5/WRF short-range forecasts and is also used as a baseline for assessing the MM5/WRF 3D-Var satellite data sensitivity experiments. As compared to the observation, the MM5 (WRF) control simulations strengthened (weakened) the cross equatorial flow over southern Arabian sea near peninsular India. The forecasts from MM5 and WRF showed a warm and moist bias at lower and upper levels with a cold bias at the middle level, which shows that the convective schemes of these models may be too active during the simulation. The forecast errors in predicted wind, temperature and humidity at different levels are lesser in WRF as compared to MM5, except the temperature prediction at lower level. The rainfall pattern and prediction skill from day 1 and day 2 forecasts by WRF is superior to MM5. The spatial distribution of forecast impact for wind, temperature, and humidity from 1-month assimilation experiments during July 2006 demonstrated that on average, for 24 and 48-h forecasts, the satellite data improved the MM5/WRF initial condition, so that model errors in predicted meteorological fields got reduced. Among the experiments, MM5/WRF wind speed prediction is most benefited from QSCAT surface wind and SSM/I TPW assimilation while temperature and humidity prediction is mostly improved due to latter. The largest improvement in MM5/WRF rainfall prediction is due to the assimilation of SSM/I TPW. The assimilation of SSM/I wind speed alone in MM5/WRF degraded the humidity and rainfall prediction. In summary the assimilation of satellite data showed similar impact on MM5/WRF prediction; largest improvement due to SSM/I TPW and degradation due to SSM/I wind speed.     

Fast data assimilation using a nonlinear Kalman filter and a model surrogate: An application to the Columbia River estuary

A fast and adjoint-free nonlinear data assimilation (DA) system was developed to simulate 3D baroclinic circulation in estuaries, leveraging two recently developed technologies: (1) a nonlinear model surrogate that executes forward simulation three orders of magnitude faster than a forward numerical circulation code and (2) a nonlinear extension to the reduced-dimension Kalman filter that estimates the state of the model surrogate. The noise sources in the Kalman filter were calibrated using empirical cross-validation and accounted for errors in model and model forcing.The DA system was applied to assimilate in situ measurements of water levels, salinities, and temperatures in simulations of the Columbia River estuary. To validate the DA results, we used a combination of cross-validation studies, process-oriented studies, and tests of statistical and dynamical consistency. The validation studies showed that DA improved the representation of several important processes in the estuary, including nonlinear tidal propagation, salinity intrusion, estuarine residual circulation, heat balance, and response of the estuary to coastal winds.     

Fast data assimilation using a nonlinear Kalman filter and a model surrogate: An application to the Columbia River estuary

A fast and adjoint-free nonlinear data assimilation (DA) system was developed to simulate 3D baroclinic circulation in estuaries, leveraging two recently developed technologies: (1) a nonlinear model surrogate that executes forward simulation three orders of magnitude faster than a forward numerical circulation code and (2) a nonlinear extension to the reduced-dimension Kalman filter that estimates the state of the model surrogate. The noise sources in the Kalman filter were calibrated using empirical cross-validation and accounted for errors in model and model forcing.The DA system was applied to assimilate in situ measurements of water levels, salinities, and temperatures in simulations of the Columbia River estuary. To validate the DA results, we used a combination of cross-validation studies, process-oriented studies, and tests of statistical and dynamical consistency. The validation studies showed that DA improved the representation of several important processes in the estuary, including nonlinear tidal propagation, salinity intrusion, estuarine residual circulation, heat balance, and response of the estuary to coastal winds.     

Identification of the optimal integration of the hydrologic model with meteorology data using the advanced GIS object-oriented programming

An object-oriented framework called GIS and Hydrologic Information System Modeling Object (GHISMO) that is tightly coupled with a prototype geographic information system (GIS), is presented in this research. The proposed GHISMO framework is used to simulate a hydrologic system on the base of the storage-release concept and with multiple weather data such as gage and NEXRAD data. Both quantitative and qualitative results of simulation for the Pleasant Run Creek and Little Buck Creek watersheds in Indiana (USA) demonstrate the method prospects.     

基于雷达资料四维变分同化及云模式的中尺度对流系统数值临近预报试验

基于雷达资料快速刷新四维变分同化(RR4DVar)初始化的三维数值云模式,利用京津冀6部新一代多普勒天气雷达和区域自动气象站观测资料,针对2013年7月4日出现在京津冀平原地区的中尺度对流系统(MCS),开展了数值临近预报试验。研究结果表明,充分考虑雷达观测信息的对流尺度数值临近预报具有很大的优势,但也存在不足:(1)模式能够较好地把握中尺度对流系统的组织发展和移动演变特征,对风暴回波带的走向和尺度特征有较好的预报,但对强回波的强度和位置预报存在一定偏差;(2)模式预报可以反映风暴系统的中小尺度扰动特征,对风暴冷池和出流边界(阵风锋)的发展变化均有较为合理的预报;(3)模式对强降水中心和雨带位置的预报有很大优势,能较好地预报弱降水雨带的分布形势和雨量,但对强降水落区的预报偏大;(4)模式对风暴造成的对流性强降水的预报准确率较高,对0.5—10 mm阈值的降水范围预报偏差比较合理,对10 mm以上降水范围的预报偏大,但是对弱降水风暴的弱回波较强回波的预报性能要好;(5)由于三维数值云模式对京津冀复杂地形的处理不够完善,对山前风场预报偏差较大,造成对山前风暴的发展演变和山前降水的预报偏差较大。     

IAP94陆面过程模式在淮河流域的验证试验

利用ＩＡＰ９４陆面过程模式对淮河流域不同季节（春、夏、秋）、不同下垫面（旱地、森林、水田）进行了单点数值模拟试验，结果表明ＩＡＰ９４不但能较好地模拟出观测的地气间各种能量通量以及植被的冠层温度，同时还很好地模拟出整层土壤含水量的变化趋势，从而证实了ＩＡＰ９４陆面过程模式具有正确描述东亚半湿润季风区陆气相互作用的能力。另外比较结果还揭示了水田表面强烈蒸发作用的重要性，需要在陆过程模式中予以重视。     

A simple model of dry convective helical vortices (with applications to the atmospheric dust devil)

An asymptotic solution of inviscid Boussinesq equations for a ‘dry convective Rankine vortex’ with prescribed buoyant forcing is given. The obtained vortex solution demonstrates monotonic growth with height of the vortex core radius, which becomes infinite at a certain critical altitude, and the corresponding attenuation of the vertical vorticity. This idealized vortex is then embedded in a convectively unstable boundary layer; the resulting approximate vortex solution has been applied to determine the maximum rotational wind speed and diameter of dry convective dust-devil-like vortices.     

A method of correlating wind data between two stations with application to the Alberta oil sands*

Abstract

Two years of hourly wind data from two stations, Fort McMurray and Mildred Lake, in the Alberta Oil Sands area are examined by means of contingency tables of simultaneous occurrences of various wind‐speed class and wind‐direction combinations. Statistical calculations with the tabulated data show that in 1974–75 the same wind‐speed class occurred simultaneously 49% of the time, and the same or next adjacent 8‐point wind direction occurred simultaneously 75% of the time. A matrix of relative frequencies of simultaneous occurrences is calculated to express this correlation between the two stations.

The 1974–75 period is shown to have statistically significant differences from a longer‐term (1963–75) set of data at Fort McMurray, pointing out the need for a method, described herein, of generating wind data that Will represent climatological conditions as well as account for the differences between two locations. The resulting synthetic wind data can be used as input to the Climatological Dispersion Model(Cdm) in order to compute climatological values of ground‐level pollutant concentration.     

On theoretical wind speed and temperature profiles over the sea with applications to data from sable island,Nova Scotia

Abstract

In October 1985, the Boundary‐Layer Research Division of the Atmospheric Environment Service conducted an experiment on Sable Island, Nova Scotia, where 10‐m wind measurements were made at a number of locations. Wind data were also collected at 4 levels on one of the 10‐m masts and at 6 levels on a 26‐m mast, both located on the South Beach. Other data used in the present study consisted of air temperature measured at 9 m and sea temperature measured at the beach.

The theory for wind speed and temperature profiles over the sea is reviewed. A method of deriving over‐sea profile parameters (u_*, θ_*, Z₀, L) from wind data at one level and the air‐sea temperature difference is described. The method is limited to applications either over homogeneous open ocean or, provided measurements are taken above the internal boundary layer generated by the change of roughness at the coastline, over a flat beach (without coastal orography). The heights at which the method is applied must be within the surface layer which must not have any discontinuities in wind speed or temperature in the vertical, such as are often associated with inversion layers. An application to data collected at beach sites in onshore flow during the October 1985 experiment is illustrated.

Once the above parameters are obtained, theoretical wind profiles may be computed and compared with observed profiles. In order to make a proper comparison it is essential to account for internal boundary layers generated at the shoreline by the step‐change in surface roughness. Only the data measured above the internal boundary layer are representative of over‐sea conditions and may, therefore, be used for verifying the theoretical profiles. The agreement between calculated and measured data is generally very good. One complication, however, is a slight upstream‐blockage effect due to a 7‐m high dune located about 140 m downwind of the 26‐m mast. Estimates of the magnitude of this effect partially account for small discrepancies in the results at the 26‐ and 10‐m mast locations.

An estimation of the most probable errors in the calculated parameters, based on assumed measurement errors, is included in the computer program. Results suggest that small measurement errors can explain the above discrepancies.