一阶Buterworth递归式带通滤波器技术改进方案

把一阶Ｂｕｔｔｅｒｗｏｒｔｈ递归带通滤波器应用于滑动样本,输出序列出现了振幅和位相的相对误差,为了克服这种现象给制作中长期预报带来的困难,作者通过数值试验,分析了该滤波器的误差传递特性,并设计了“拟定初值的三步滤波方案”。大量试验表明,该方案可较为有效地抑制输出序列的误差和位相飘移。     

青藏高原区域性积雪增量序列及其变化特征

利用青藏高原83个气象台站的逐日积雪观测资料,充分考虑每次降雪过程所引起的积雪变化,建立了该地区1960-2007年度的区域性积雪增量序列,并讨论了其长期变化特征.结果表明:(1)积雪增量序列具有其他积雪参数指标的代表性,又较其他积雪参数指标序列有明显的统计学优点.(2)在48个年度里,1960-1966年度青藏高原区域性积雪增量处于负位相,1967-1998年度为多雪的正位相,1999-2007年度又为少雪的负位相,其中1998年度是个具有显著性的突变点.     

长期大气折射指数预报研究

本文对近地面大气折射指数气候变化统计预报进行了研究,建立了逐步回归多重分析时变模式。在模式中考虑了大气折射指数肘间序列的振动周期、前期预报因子(降水量),还特别考虑了时间序列的振动位相,并随着新数据的到来,适时地改变预报因子和预报方程的参数。本文用该模式对上海近地面6月上旬旬平均大气折射指数进行了长期预报试验,结果表明预报误差较小,预报效果稳定,具有较好的实用价值。这方面的工作可用来对无线电波的传播环境进行预测,从而对电子工程系统进行适时的调整控制。除此而外,该预报模式也可以应用于其它的预报领域。     

ＬＣ滤波器设计的优化解析公式

给出了改进的Ｃｈｅｂｙｓｈｅｖ和Ｂｕｔｔｅｒｗｏｒｔｈ低通原型ＬＣ滤波器的解析设计公式．此公式充分利用阶数取整带来的余量,在不改变滤波器阶数的前提下减小了通带波动,极大地方便了集成滤波器的实现．给出了集成有源高通滤波器设计实例,验证了改进的Ｃｈｅｂｙｓｈｅｖ和Ｂｕｔｔｅｒｗｏｒｔｈ低通原型ＬＣ滤波器解析设计公式的有效性．     

平流方程的数值研究

二维平流方程分别用二阶差分方案、四阶差分方案、一维三次样条函数方案及谱导数方案(假谱方案)进行长时间的数值积分,其中风场是给定的均匀旋转风场或形变风场,并分别以不同大小的圆锥状的物理量分布作为初值。数值试验的结果表明,谱导数方案的精度最高,没有位相误差;其次是一维三次样条函数方案;再其次是四阶差分方案;二阶差分方案的精度最差,畸变严重,位相误差大。 对上述四种方案给出用权重系数形式表示的统一的导数计算公式,看到高精度的导数计算公式是非局地性的。     

边界层方案对中国中东部地区降水日变化的影响

以华中区域中尺度业务模式WRF3D为平台,使用MYJ、ACM2边界层方案完成了2012年7月的批量敏感试验,并与高时空分辨率实况降水数据作对比,重点关注不同边界层方案对中国中东部地区降水日变化的影响。结果表明,两种方案下的降水预报均对我国中东部地区清晨和午后并存的降水双峰值有所表现,且两者对清晨降水峰值的预报差异不大,但预报的午后降水峰值差异明显,尤其表现在峰值幅度上,位相上也略有差异。总体而言,ACM2方案下的降水日变化特征更接近实况。进一步对高分辨率模式输出的诊断表明,午后降水峰值主要由隐式降水决定。在此基础上,从模式中隐、显式降水产生机制的角度考察了不同边界层方案下降水日变化差异的可能原因。     

基于WRF模式输出统计的逐时太阳总辐射预报初探

基于WRF(weather research and forecasting model)模式逐时输出结果,设计了逐时太阳总辐射的模式输出统计(model output statistics,MOS)预报流程。主要包括:对逐时观测序列进行低通滤波再除以天文辐射,对模式输出因子的筛选和降维,以及建立MOS预报方程,并对2009年1月、4月、8月和10月武汉站逐时太阳总辐射进行预报试验。结果表明,该方案在各月预报相对稳定,拟合和预报效果均较为理想,可使平均绝对百分比误差控制在20%～30%,相对均方根误差控制在30%～40%,相对模式直接预报辐射改进了50%左右。由此可见,通过对模式输出进行解释应用,可以有效提高辐射预报的准确率。此外,客观分析所得的气温、云量、露点、比湿、相对湿度、地面气压等13个模式输出因子可以作为各地区建立MOS辐射预报方程的参考因子。     

对近期提出的太阳—QBO—天气之间关系的看法

联系太阳活动,考察了1875—1936年期间若干测站的地面气象资料,特别是想看一看这些历史资料是否与 van Loon 和 Labitzke(vLL)最近根据现代资料(1950年以后)所找出的太阳-QBO-天气之间的关系相一致。在把 vLL 的分析往前延伸时,主要问题是不知道 QBO 的位相。在该研究中,对 QBO 位相的数百万种可能的序列重复了 vLL 的计算。其结果揭示出在早期资料中再现 vLL 的结果时存在的一些问题。这表明,或者是 QBO 活动在过去是不同的,或者是vLL关于太阳-天气关系的研究结果对长序列资料是不稳定的。     

CMIP5模式对北极涛动的模拟评估

基于国际耦合模式比较计划第五阶段(CMIP5)历史试验(historical)的输出结果,评估了26个耦合模式对北极涛动(AO)的模拟能力。对各模式逐月AO指数序列的功率谱分析表明,有23个模式能够模拟出AO模态无显著变化周期的特征。这些模式也能够较好地再现冬季AO在海平面气压场上的主要分布特征,但均高估了AO模态的强度。对于伴随着冬季AO位相变化而出现的中高纬度偶极型的纬向平均纬向风异常,CMIP5中只有一些模式有较好的模拟表现,大多数模式对其中心位置和强度的模拟存在明显不足。对比之下,MPI-ESM-P对AO时空特征的模拟更接近观测结果。     

SST日变化对西太平洋暖池海表热通量季节内变化的影响(英文)

基于利用日最大太阳辐射、日平均海面风和日降水量近似计算海表温度（SST）日变化的振幅的方法,发展了一个计算SST日循环的参数化方案。利用周平均SST强迫美国国家大气研究中心（NCAR）的CCM3大气模式进行了有、无考虑SST日变化的比较试验。热带海洋与全球大气计划之耦合海气响应实验（TOGA COARE）的强化观测期间 IMET浮标的逐时海表观测资料不仅验证了该参数化方案的合理性,也表明了利用参数化方案对强迫场SST迭加日变化使CCM3较真实地模拟出西太平洋暖池海表热通量的季节内变化的位相结构。     

Two-Dimensional Variational Analysis of Near-Surface Moisture from Simulated Radar Refractivity-Related Phase Change Observations

Because they are most sensitive to atmospheric moisture content, radar refractivity observations can provide high-resolution information about the highly variable low-level moisture field. In this study, simulated radar refractivity-related phase-change data were created using a radar simulator from realistic highresolution model simulation data for a dryline case. These data were analyzed using the 2DVAR system developed specifically for the phase-change data. Two sets of experiments with the simulated observations were performed, one assuming a uniform target spacing of 250 m and one assuming nonuniform spacing between 250 m to 4 km. Several sources of observation error were considered, and their impacts were examined. They included errors due to ground target position uncertainty, typical random errors associated with radar measurements, and gross error due to phase wrapping. Without any additional information, the 2DVAR system was incapable of dealing with phase-wrapped data directly. When there was no phase wrapping in the data, the 2DVAR produced excellent analyses, even in the presence of both position uncertainty and random radar measurement errors. When a separate pre-processing step was applied to unwrap the phase-wrapped data, quality moisture analyses were again obtained, although the analyses were smoother due to the reduced effective resolution of the observations by interpolation and smoothing involved in the unwrapping procedure. The unwrapping procedure was effective even when significant differences existed between the analyzed state and the state at a reference time. The results affirm the promise of using radar refractivity phase-change measurements for near-surface moisture analysis.     

Autoregressive filtering versus linear detrending in estimation of fluxes by the eddy covariance method

The application of autoregressive running mean filtering (RMF) and linear detrending (LDT) in the estimation of turbulent fluxes by the eddy covariance method is analysed. The systematic, as well as the random, errors of the fluxes arising from filtering and/or limited observation time effects are described. To observe negligible systematic errors in fluxes, the RMF has to be applied with moderately long time constants. However, the obtained flux values are subject to increased random errors during periods of non-stationarity and the method leads to systematic overestimation of variances. These shortcomings are not inherent in the LDT approach, which is recommended for use. But the systematic errors of fluxes due to LDT are not negligible under certain experimental conditions and have to be accounted for. The corrections are important because the relatively small errors in short-period fluxes can translate to significant errors in long-period averages. The corrections depend on the turbulence time scales, which should be preferably estimated as ensemble mean variables for a particular site.     

附加相位延迟到弯曲角反演过程中的误差传递模型

为了研究全球卫星导航定位系统无线电掩星探测中的误差传递过程,文中讨论了掩星探测中附加相位延迟到弯曲角的反演过程中,附加相位延迟随机误差对反演弯曲角误差的影响。在对全球卫星导航定位系统和低轨道卫星进行精密定轨的情况下,忽略轨道误差对反演结果的影响,将轨道参数视为常数,考虑附加相位延迟△L存在随机误差δ△L时对反演弯曲角α带来的弯曲角误差δα。利用EGOPS仿真软件,采用数值模拟的方法模拟了无电离层影响的大气附加相位延迟数据,利用CHAMP真实轨道数据,MSIS-90大气模型和全三维射线追踪器进行前向建模,从中选取了一个掩星事件对误差传递模型进行了验证。结果表明:附加相位延迟的变化会给弯曲角的反演带来相应的误差:当随机误差为5 cm左右时,反演的弯曲角误差在±0.2 mrad范围内;当随机误差为1 cm左右时,反演的弯曲角误差在±4μrad范围内;当随机误差为1 mm左右时,反演的弯曲角误差在±40μrad范围内。最后从模拟的掩星事件中随机抽取了20个掩星事件进行统计分析,给出了20个掩星事件存在1 cm附加相位延迟测量误差引起弯曲角误差的均值误差,其误差特性虽略大于个例研究结果,但与模型完全吻合,反演的弯曲角误差也在±40μrad范围内,进一步验证了模型的准确性。     

Displaced-Beam Small Aperture Scintillometer Test. Part I: The Wintex Data-Set

The friction velocity (u_*) and the sensible heat flux density (H) determined with a displaced-beam small aperture scintillometer (DBSAS) and a hot-film eddy correlation system are compared. Random errors in the DBSAS are relatively small, compared to scatter found with two eddy-correlation systems. Assuming that the hot-film system yields the true fluxes, theDBSAS appears to overestimate u_* when u_* is less than 0.2 m s^-1 and to underestimate u_* at high wind speeds. This implies that the DBSAS measurements of theinner scale length of turbulence, l₀, a direct measure for the dissipation rate of kinetic turbulent energy, are biased. Possible causes for these results are discussedin detail. A correction procedure is presented to account for effects of random noise and of so-called inactive turbulence or sensor vibrations. The errors in u_* cause errors in the DBSAS measurements of the structure parameter of temperature C_T ². The derived H appears to be less sensitive to errors in l₀ and C_T ², because errors in these quantities tend to cancel out.     

Summer temperature in the eastern part of southern South America: its variability in the twentieth century and a teleconnection with Oceania

The 1907–2001 summer-to-summer surface air temperature variability in the eastern part of southern South America (SSA, partly including Patagonia) is analysed. Based on records from instruments located next to the Atlantic Ocean (36°S–55°S), we define indices for the interannual and interdecadal timescales. The main interdecadal mode reflects the late-1970s cold-to-warm climate shift in the region and a warm-to-cold transition during early 1930s. Although it has been in phase with the Pacific Decadal Oscillation (PDO) index since the 1960s, they diverged in the preceding decades. The main interannual variability index exhibits high spectral power at ~3.4 years and is representative of temperature variability in a broad area in the southern half of the continent. Eleven-years running correlation coefficients between this index and December-to-February (DJF) Niño3.4 show significant decadal fluctuations, out-of-phase with the running correlation with a DJF index of the Southern Annular Mode. The main interannual variability index is associated with a barotropic wavetrain-like pattern extending over the South Pacific from Oceania to SSA. During warm (cold) summers in SSA, significant anticyclonic (cyclonic) anomalies tend to predominate over eastern Australia, to the north of the Ross Sea, and to the east of SSA, whereas anomalous cyclonic (anticyclonic) circulation is observed over New Zealand and west of SSA. This teleconnection links warm (cold) SSA anomalies with dry (wet) summers in eastern Australia. The covariability seems to be influenced by the characteristics of tropical forcing; indeed, a disruption has been observed since late 1970s, presumably due to the PDO warm phase.     

Optimizing the Determination of Roughness Parameters for Model Urban Canopies

We present an objective optimization procedure to determine the roughness parameters for very rough boundary-layer flow over model urban canopies. For neutral stratification the mean velocity profile above a model urban canopy is described by the logarithmic law together with the set of roughness parameters of displacement height d, roughness length \(z_0\), and friction velocity \(u_*\). Traditionally, values of these roughness parameters are obtained by fitting the logarithmic law through (all) the data points comprising the velocity profile. The new procedure generates unique velocity profiles from subsets or combinations of the data points of the original velocity profile, after which all possible profiles are examined. Each of the generated profiles is fitted to the logarithmic law for a sequence of values of d, with the representative value of d obtained from the minima of the summed least-squares errors for all the generated profiles. The representative values of \(z_0\) and \(u_*\) are identified by the peak in the bivariate histogram of \(z_0\) and \(u_*\). The methodology has been verified against laboratory datasets of flow above model urban canopies.     

Validation of synoptic circulation patterns simulated by the Canadian climate centre general circulation model for western north America: Research note

Abstract

In order to assess the ability of a GCM to simulate regional to synoptic scale atmospheric structures, a correlation‐based computer‐assisted gridded map typing procedure is used to compare daily pressure (MSL) and geopotential height fields (500 hPa) from a GCM simulation of the present climate to a decade of NMC analyses. The model is able to reproduce the entire range of synoptic circulation types. However, statistically significant differences in the seasonal frequencies and variances of the main circulation types are evident. These differences, which are most pronounced in the winter (at 500 hPa) and in spring and autumn in the MSL fields, are consistent with subtle errors in the predicted fields at the hemispheric scale. The lack of agreement between the NMC climatology and the “control” simulation precludes extension of this approach to investigation of climate change impacts in western north America, and to more meteorologically dynamic extra‐tropical regions. The map‐typing procedure is shown to be an appropriate GCM synoptic‐scale validation tool that permits direct comparison of GCM output and observed fields. As such, it has the potential to elucidate the regional‐scale impacts of global climatic change through established synoptic circulation environment relationships.     

基于Parsivel天气现象仪资料的郑州“7·20”罕见特大暴雨微物理特征分析

基于2016-2018年ECMWF模式温度预报和浙江省72个国家基本站观测资料,根据温度日变化特征,采用K-近邻(KNN)回归算法进行误差订正,改进浙江省172 h精细化温度预报。在KNN回归算法中,将模式起报时刻的温度视作“背景”,由模式预报减去起报时刻温度消除“背景”影响,得到温度日变化曲线,通过温度日变化曲线构建差异指标,选取历史相似个例。根据历史相似个例的误差特征,对温度预报进行订正,得到改进的温度预报。检验结果表明,KNN方案的温度预报平均绝对误差较ECMWF和30 d滑动平均误差订正方案(OCF)的分别减小26.2%和5.2%;日最高和最低温度预报误差绝对值小于2℃,准确率较ECMWF的分别提高14.8%和4.3%,较OCF的分别提高3.0%和1.3%。KNN方案对地形复杂地区的温度预报改进效果更为明显,对冷空气活动和夏季高温等天气过程预报改善效果也较稳定。     

Near-surface wind estimates using statistics from a planetary boundary-layer model

This paper shows the possibilities of a procedure for estimating near-surface wind statistics, by means of the numerical integration of a simple boundary-layer model with a second-order turbulent closure. Standard and easily available synoptic data are used as initial and boundary conditions. The development of this methodology is impelled by increasing requirements of a quick and precise knowledge of the wind characteristics in many regions of South America, which confronts the serious limitation of a reduced number of extended observational series, scattered over a vast continent. In order to evaluate the methodology, near-surface wind statistics from observed data at two locations are compared with model output statistics. Relative errors are about 0.2 for daily mean velocities and about 0.1 for weekly mean velocities, observed and computed time series being highly correlated in both cases. Calculated frequency distribution of wind directions is in good agreement with the observed one, and the absolute mean error in the daily mean wind direction is about 20 deg. Even though a wide variety of large-scale synoptic situations has been indirectly considered through boundary conditions, basic model output statistics resemble fairly well those observed at different levels between the surface and 100 m.