AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.     

Classification of Northeast China Cold Vortex Activity Paths in Early Summer Based on K-means Clustering and Their Climate Impact

The classification of the Northeast China Cold Vortex(NCCV)activity paths is an important way to analyze its characteristics in detail.Based on the daily precipitation data of the northeastern China(NEC)region,and the atmospheric circulation field and temperature field data of ERA-Interim for every six hours,the NCCV processes during the early summer(June)seasons from 1979 to 2018 were objectively identified.Then,the NCCV processes were classified using a machine learning method(k-means)according to the characteristic parameters of the activity path information.The rationality of the classification results was verified from two aspects,as follows:(1)the atmospheric circulation configuration of the NCCV on various paths;and(2)its influences on the climate conditions in the NEC.The obtained results showed that the activity paths of the NCCV could be divided into four types according to such characteristics as the generation origin,movement direction,and movement velocity of the NCCV.These included the generation-eastward movement type in the east of the Mongolia Plateau(eastward movement type or type A);generation-southeast longdistance movement type in the upstream of the Lena River(southeast long-distance movement type or type B);generationeastward less-movement type near Lake Baikal(eastward less-movement type or type C);and the generation-southward less-movement type in eastern Siberia(southward less-movement type or type D).There were obvious differences observed in the atmospheric circulation configuration and the climate impact of the NCCV on the four above-mentioned types of paths,which indicated that the classification results were reasonable.     

Phase Two of the Integrative Monsoon Frontal Rainfall Experiment (IMFRE-II) over the Middle and Lower Reaches of the Yangtze River in 2020

Phase Two of the Integrative Monsoon Frontal Rainfall Experiment(IMFRE-II)was conducted over the middle and lower reaches of the Yangtze River during the period 16 June to 19 July 2020.This paper provides a brief overview of the IMFRE-II field campaign,including the multiple ground-based remote sensors,aircraft probes,and their corresponding measurements during the 2020 mei-yu period,as well as how to use these numerous datasets to answer scientific questions.The highlights of IMFRE-II are:(1)to the best of our knowledge,IMFRE-II is the first field campaign in China to use ground-based,airborne,and spaceborne platforms to conduct comprehensive observations over the middle and lower reaches of the Yangtze River;and(2)seven aircraft flights were successfully carried out,and the spectra of ice particles,cloud droplets,and raindrops at different altitudes were obtained.These in-situ measurements will provide a“cloud truth”to validate the ground-based and satellite-retrieved cloud and precipitation properties and quantitatively estimate their retrieval uncertainties.They are also crucial for the development of a warm(and/or cold)rain conceptual model in order to better understand the cloud-to-rain conversion and accretion processes in mei-yu precipitation events.Through an integrative analysis of ground-based,aircraft,and satellite observations and model simulations,we can significantly improve our cloud and precipitation retrieval algorithms,investigate the microphysical properties of cloud and precipitation,understand in-depth the formation and dissipation mechanisms of mei-yu frontal systems,and improve cloud microphysics parameterization schemes and model simulations.     

CHARACTERISTICS OF TROPICAL CYCLONES IMPACTING A NUCLEAR POWER PLANT IN SOUTHERN FUJIAN AND ESTIMATE OF MINIMUM CENTRAL PRESSURE

Based on tropical cyclone (TC) data for the period 1949 to 2008 and following the Gumbe-I method, Pearson-Ⅲ method and determinacy method, this article estimates the possible minimum central pressure of TCs affecting southern Fujian where a nuclear power will be located. Results show that the observed minimum central pressure of TCs agrees well with the results determined with the methods above and there is little difference between them (the minimum central pressure is 867.4 hPa and 868.1 hPa, respectively, in a 1,000-yr return period). Established with the theory of atmospheric dynamics, the determinacy method yields a result of 867.28 hPa/1000 years, only a little smaller than the result of the probability method. Because of randomicity in parameter adjustment with the Pearson-Ⅲ method whereas the determinacy method is theoretically solid and its estimates are the smallest of the three methods, it is therefore reasonable, for security and conservative concerns, to adopt the result determined with the determinacy method as the possible maximum intensity of TC (with the central pressure being 867.28 hPa in a 1,000-yr return period).     

A Comparison of Two Bulk Microphysics Parameterizations for the Study of Aerosol Impacts on an Idealized Supercell

Idealized supercell storms are simulated with two aerosol-aware bulk microphysics schemes(BMSs),the Thompson and the Chen-Liu-Reisner(CLR),using the Weather Research and Forecast(WRF)model.The objective of this study is to investigate the parameterizations of aerosol effects on cloud and precipitation characteristics and assess the necessity of introducing aerosols into a weather prediction model at fine grid resolution.The results show that aerosols play a decisive role in the composition of clouds in terms of the mixing ratios and number concentrations of liquid and ice hydrometeors in an intense supercell storm.The storm consists of a large amount of cloud water and snow in the polluted environment,but a large amount of rainwater and graupel instead in the clean environment.The total precipitation and rain intensity are suppressed in the CLR scheme more than in the Thompson scheme in the first three hours of storm simulations.The critical processes explaining the differences are the auto-conversion rate in the warm-rain process at the beginning of storm intensification and the low-level cooling induced by large ice hydrometeors.The cloud condensation nuclei(CCN)activation and auto-conversion processes of the two schemes exhibit considerable differences,indicating the inherent uncertainty of the parameterized aerosol effects among different BMSs.Beyond the aerosol effects,the fall speed characteristics of graupel in the two schemes play an important role in the storm dynamics and precipitation via low-level cooling.The rapid intensification of storms simulated with the Thompson scheme is attributed to the production of hail-like graupel.     

THE IMPACTS OF MADDEN-JULIAN OSCILLATION ON SPRING RAINFALL IN EAST CHINA

Phase composite analyses are conducted to investigate the possible effect of the Madden–Julian oscillation(MJO)on the spring rainfall anomalies in East China by using the Real-time Multivariate MJO(RMM)index from Australian Meteorological Bureau.The results show that the rainfall anomalies over the mid-and lower-valley of Yangtze River are positive when the MJO shifts eastward to the mid-and eastern-Indian Ocean,and anomalous precipitation over South China are positive when the MJO moves further eastward to the maritime continent,whereas spring rainfall anomalies over East China are negative in the other MJO episodes.The MJO impacts on the precipitation over East China result from the changes in large-scale atmospheric circulation as well as vorticity and water vapor transportation in the mid-and lower-troposphere.     

Hydro-climatic Characteristics of Yarlung Zangbo River Basin since the Last Glacial Maximum

Global climate changes significantly impact the water condition of big rivers in glacierized high mountains. However,there is a lack of studies on hydrological changes within river basins caused by climate changes over a geological timescale due to the impossibility of direct observations. In this study, we examine the hydro-climatic variation of the Yarlung Zangbo River Basin in the Tibet Plateau since the Last Glacial Maximum(LGM) by combining δ18 O proxy records in Indian and Omani caves with the simulated Indian summer monsoon, surface temperature, precipitation, evapotranspiration and runoff via the Community Climate System Model and the reconstructed glacier coverage via the Parallel Ice Sheet Model. The mean river runoff was kept at a low level of 145 billion cubic meters per year until an abrupt increase at a rate of 8.7 million cubic meters per year in the B?lling-Aller?d interval(BA). The annual runoff reached a maximum of 250 billion cubic meters in the early Holocene and then reduced to the current value of 180 billion cubic meters at a rate of 6.4 million cubic meters per year. The low runoff in the LGM and Heinrich Stadial 1(HS1) is likely attributed to such a small contribution of precipitation to runoff and the large glacier cover. The percentage of precipitation to runoff was only 20%during the LGM and HS1. Comparison of glacier area among different periods indicates that the fastest deglaciation occurred during the late HS1, when nearly 60% of glacier area disappeared in the middle reach, 50% in the upper reach,and 30% in the lower reach. The rapid deglaciation and increasing runoff between the late HS1 and BA may have accelerated widespread ice-dam breaches and led to extreme outburst flood events. Combining local geological proxy records and regional simulations could be a useful approach for the study of paleo-hydrologic variations in big river basins.     

Verification of Subseasonal-to-Seasonal Forecasts for Major Stratospheric Sudden Warmings in Northern Winter from 1998/99 to 2012/13

This study reports verification results of hindcast data of four systems in the subseasonal-to-seasonal(S2 S) prediction project for major stratospheric sudden warmings(MSSWs) in northern winter from 1998/99 to 2012/13. This report deals with average features across all MSSWs, and possible differences between two MSSW types(vortex displacement and split types). Results for the average features show that stratospheric forecast verifications, when further averaged among the four systems, are judge...     

南京三千公尺高空之风向与天气之预测

晚近日本籐原笑平(Fujiwhora)博士于地球物理杂志发表「根据三千公尺高空等压线,以预测天气之一例证」一文,谓日本最近用三千公尺高空之等压线,作每日天气之预测,已得相当成就。氏之经验法则谓自九月以迄五月,日本太平洋沿岸,三千公尺高空之等压线,来自西南者,可形去致雨,而来自西北者,则可期晴明。此种倾向颇为显著。     

The Distribution and Uncertainty Quantification of Wind Profile in the Stochastic General Ekman Momentum Approximation Model

The general Ekman momentum approximation boundary-layer model(GEM) can be effectively used to describe the physical processes of the boundary layer. However, eddy viscosity, which is an approximated value, can lead to uncertainty in the solutions. In this paper, stochastic eddy viscosity is taken into consideration in the GEM, and generalized polynomial chaos is used to quantify the uncertainty. The goal of uncertainty quantification is to investigate the effects of uncertainty in the eddy viscosity on the model and to subsequently provide reliable distribution of simulation results. The performances of the stochastic eddy viscosity and generalized polynomial chaos method are validated based on three different types of eddy viscosities, and the results are compared based on the Monte Carlo method. The results indicate that the generalized polynomial chaos method can be accurately and efficiently used in uncertainty quantification for the GEM with stochastic eddy viscosity.     

台风“巨爵”近海急剧加强的特征及成因分析

源自东风波动的0915号台风"巨爵"在自东向西靠近广东沿海时出现了近海急剧加强的异常现象。利用热带气旋定位资料、NCEP全球同化分析资料、卫星云图以及海上平台自动站探测资料等对"巨爵"近海加强特征和成因进行了诊断分析。结果表明:"巨爵"形成源自于西太副高南侧的东风波动,副高演变及引起的"北高南低"和"东高西低"的形势与"巨爵"环流增强发展密切相关;弱的垂直风切变、高空明显的辐散出流以及台风中心正涡度增大并向对流层中上层拓展有利于"巨爵"强度增强;低层弱冷空气南下,一方面从增加气压梯度和温度梯度来加强台风环流和暖心结构,另一方面通过触发外围对流云团发展并往台风中心输送进而对"巨爵"近海加强产生重要的触发作用;孟加拉湾西南季风气流、105°~110°E附近越赤道气流和西太副高西侧东南气流3支气流汇合为"巨爵"近海加强提供了充沛水汽条件。     

温室黄瓜生育期模拟模型的研究

根据黄瓜（Cucumis sativus）发育的光温反应过程,建立了基于生理发育时间的黄瓜生育期模拟模型,并利用不同品种的试验资料对模型进行了检验。结果表明,黄瓜从播种到第一批瓜成熟需积累生理发育时间30d,播种—幼苗、幼苗—伸蔓、伸蔓—开花、开花—结瓜、结瓜—成熟所需的生理发育时间分别为3、11、9、5、2d。发芽期、幼苗期、伸蔓期、开花期、结瓜期等各生育期持续时间的模拟值与实际观测值的回归估计标准误差（root mean squared error;RMSE）分别为0、2．6、1．7、0．8、2．1d,从播种到第一次采收的模拟值与观测值的RMSE为1．4d。而用有效积温法对发芽期、幼苗期、伸蔓期、开花期、结瓜期等各生育期持续时间的模拟值与实际观测值的RMSE分别为0．7、10．0、5．7、2．4、2．5d,从播种到第一次采收的模拟值与观测值的RMSE为22．0d。本模型比用有效积温法能更准确地预测黄瓜各个生育期的起止日期和黄瓜的收获期。     

T106数值预报产品在赣南温度预报中的释用

利用1995年7～9月、1996年7～8月24～144h、1000～500hPa的T106格点资料与赣南17个县（市）的最高、最低、平均温度作相关分析;用逐步回归方法建立了分县逐日滚动预报方程。并用1996年9月份T106格点资料对方法进行检验,同时对预报方法及T106产品在温度预报中的释用能力作了初步分析。     

雹云判别和适宜作业时机分析

人工防雹作业中关键是雹云判别和适时作业。雹云判别用4因子组合表征雹云特征，还能判别雹落区；在cb云路径上确定查强区进行雹云判别，以决定适时作业的时间。     

RG13型雨量传感器测量结果的不确定度评定

国内民航机场主要使用的雨量观测设备为芬兰维萨拉公司生产的RG13型雨量传感器,为保证雨量测量数据的真实可靠,对其测量结果的不确定度分析很有必要。根据自动气象站现场校准方法,分别进行大雨强和小雨强的重复测试,并依据JJF1059.1-2012测量不确定度的评定与表示要求,进行A类不确定度评定。分析测量过程中的B类不确定度来源,进行B类评定,最终给出扩展不确定度。结果表明：在小雨强下,测量不确定度为U95=0.17mm,包含因子k=2。在大雨强下,测量不确定度为U95=0.16mm,包含因子k=2。该研究完善了雨量传感器的现场校准工作流程,对雨量传感器测量结果的可信度评定具有参考价值。     

西风带南支槽对云南天气的影响

利用1980-2008年近30年逐日地面、高空观测资料,统计了影响云南的西风带南支槽个例,并分析了南支槽的时空分布和对云南降水的影响。结果表明,影响云南的南支槽平均每年出现18.76次,11月-次年5月平均每月出现次数相当(6-10月西风带北撤,转换为孟加拉湾槽);约5.88%的南支槽过程对云南产生大到暴雨天气,54.83%的南支槽产生小到中雨,另有17.28%的南支槽产生冰雹天气过程。南支槽的进退与西风带环流形势、副热带高压位置和高原大地形等关系密切,南支槽位置、水汽输送、湿度锋区、低空急流和冷空气强弱等条件的不同决定了降水的强弱或是否有强对流天气出现。     

Hadoop在气象数据密集型处理领域中的应用

气象资料的统计分析计算属于数据密集型计算,目前的处理方式多为单机处理,对大量数据的处理比较慢,难以应对日益增长的数据,对气象资料的研究形成一定的制约。针对数据密集型气象数据的处理,尝试应用Hadoop的MapReduce思想提高计算效率;对Hadoop在处理大量小文件组成的气象数据时的低效率,提出对原始文件进行预处理,将多个小文件整合成能直接用于计算的大文件。试验证明,该方法解决了Hadoop处理大量小文件时的低效率问题,通过与Oracle入库检索的比较,应用Hadoop处理数据密集型气象资料具有实际意义。     

梨树始花期预报

基于丰县1984—2009年气温、日照和降水量资料,分析越冬期不同时段的气象因子对梨树始花期的影响。结果表明:①后期因子影响大,前期相对较小;②温度、日照因子影响大,降水因子影响较小。越冬期温度水平高,光照充足,梨树始花期就偏早,反之始花期就偏晚;在上年12月至当年1月,除12月下旬外,最高气温偏高,影响梨树休眠期的需冷量积累,开花就越晚,最高气温偏低,需冷量的积累充足,有利于解除梨树休眠,开花就越早。采用逐步回归方法,利用SPSS软件,建立了梨树始花期预报模型,模型拟合效果较好,预报准确率满足业务要求,及时为政府和果农提供准确的始花期预报。     

降水对雾中能见度参数化的影响

雾的能见度是雾滴数浓度和含水量的函数,可以利用雾滴数浓度和含水量进行参数化得到能见度拟合值。但是当雾中有降水存在时,降水会影响雾滴数浓度和含水量,进而影响能见度参数化。2009年1月19日至3月5日在湖北恩施雷达站进行了包括雾微物理特征、雨滴谱及雨滴末速度、能见度以及基本气象要素的观测。利用不同的能见度参数化公式得到能见度拟合值,并与实测值进行对比分析,考虑雾中降水影响时的能见度拟合值要比不考虑降水时的能见度偏差大,因此在有降水伴随雾产生时,降水的效应可忽略不计。同时,利用40％的雾滴谱观测数据建立了良好的模型用于拟合当地的能见度,并利用Gultepe公式得到相应的参数,对当地的能见度进行了计算和预报。     

广义帕雷托分布在重庆暴雨强降水研究中的应用

引进广义帕雷托分布(GPD),借助于现代L-矩估计方法,模拟重庆地区极端降水事件,推算一定重现期的极端降水量分位数。模拟试验表明,基于超门限峰值法(POT)的GPD不但计算简便,而且基本不受原始序列样本量的影响,具有全部取值域的高精度稳定拟合(包括高端厚尾部),与GEV模拟结果相比,GPD具有更高精度和稳定性,更为实用。