Non-Linear Numerical Model Predictions of Flow Over an Isolated Hill of Moderate Slope

Non-linear model simulations of atmospheric boundary-layer flow over the hill called Blashaval have been compared with observations and linear model predictions. Previous studies have shown that linear models can give good predictions of wind speed at the summit and on the upwind slopes of Blashaval. The non-linear model provided wind speed predictions of similar accuracy when compared with the mean observed values at these locations.The published experimental data showed that on the lee-slope the wind speeds at 8m were reduced to approximately 10% of their upstream value at the same height. This was associated with an 180° change in wind direction compared with the upstream flow, suggesting that flow separation had occurred. The non-linear model predictions of lee-slope wind speed, when used with high-resolution topography data, were significantly better than linear model predictions. However, the non-linear model predicted lee-slope wind speeds that were still stronger than observed. The non-linear model simulated flow separation more readily with a 1 1/2-order turbulence closure than with a first-order, mixing-length closure. The configuration of the non-linear model that gave best agreement with observations predicted an 8m lee-slope wind speed that was around 50% of the upstream value.     

Non-Linear, Microscale Modelling of the Flow Over Askervein Hill

The flow solver “3DWind” is used to explore new aspects of the Askervein hill flow case. Previous work has investigated sensitivities to the grid, the inflow boundary profile, the roughness and the turbulence model. Several different linear and non-linear numerical models have also been validated by means of the Askervein hill case. This analysis focuses on the flow sensitivity to the grid spacing, the incident wind direction and the vertical resolution of topographic data. The horizontal resolution is found to be fine enough to cause only minor differences compared to a grid where every second node is removed. The vertical resolution dependence is mainly attributed to the wall functions. Simulations are performed for wind directions 200°, 205°, 210° and 215° at the reference station. The smallest directional biases compared to experimental values along a line through the hilltop are found for the directions 200° and 205°. There are larger wind direction changes along this line through the hilltop in the 200° case than in the 215° case. Still the simulation results give less veering than found in the experimental results, and this is maybe caused by a slightly stable atmosphere. The sensitivity to the vertical resolution of the topographical data is found to be particularly high close to the ground at the top of the hill; this is where the speed-up is most important. Differences decrease with the height from the ground. At higher levels the speed-ups are smaller and caused by terrain formations with larger scales.     

An improved method for integrating the Mixed Spectral Finite Difference (MSFD) model equations

In their Mixed Spectral Finite Difference (MSFD) model for flow over complex terrain, Beljaars et al. (1987) solve a set of coupled, second-order ordinary differential equations (ODEs) for the first-order perturbations to the logarithmic velocity profile caused by nonuniform surface roughness and topography. To solve this set of ODEs, they employ a Forward Euler Shooting Method. It is demonstrated here that the shooting method is computationally unstable for this problem. An absolutely stable finite-difference method based on a block tridiagonal LU factorization of the finite-difference matrix is presented. The advantages of the present algorithm over the method used by Beljaars et al. are demonstrated both by theoretical argument and numerical experiment.     

A Mixed Spectral-Integration Model for Neutral Mean Wind Flow Over Hills

A linear model for neutral surface-layer flow over orography is presented. The Reynolds-Averaged Navier-Stokes and E– turbulence closure equations are expressed in a terrain-following coordinate system created from a simple analytical expression in the Fourier domain. The perturbation equations are solved spectrally horizontally and by numerical integration vertically. Non-dimensional solutions are stored in look-up tables for quick re-use. Model results are compared to measurements, as well as other authors’ flow models in three test cases. The model is implemented and tested in two-dimensional space; the equations for a full three-dimensional version are presented.     

Laboratory Study of Topographic Effects on the Near-surface Tornado Flow Field

To study topographic effects on the near-surface tornado flow field, the Iowa State University tornado simulator was used to simulate a translating tornado passing over three different two-dimensional topographies: a ridge, an escarpment and a valley. The effect of the translation speed on maximum horizontal wind speeds is observed for translation speeds of 0.15 and 0.50 \(\hbox {m}\,\hbox {s}^{-1}\), with the lower value resulting in a larger maximum horizontal wind speed. The tornado translation over the three topographies with respect to flat terrain is assessed for changes in: (a) the maximum horizontal wind speeds in terms of the flow-amplification factor; (b) the maximum aerodynamic drag in terms of the tornado speed-up ratio; (c) the maximum duration of exposure at any location to high wind speeds of a specific range in terms of the exposure amplification factor. Results show that both the maximum wind amplification factor of 14%, as well as the maximum speed-up ratio of 14%, occur on the ridge. For all topographies, the increase in aerodynamic drag is observed to be maximized for low-rise buildings, which illustrates the importance of the vertical profiles of the horizontal wind speed near the ground. The maximum exposure amplification factors, estimated for the range of wind speeds corresponding to the EF2 (50–60 \(\hbox {m}\,\hbox {s}^{-1}\)) and EF3 (61–75 \(\hbox {m}\,\hbox {s}^{-1})\) scales, are 86 and 110% for the ridge, 4 and 60% for the escarpment and ? 6 and 47% for the valley, respectively.     

Spectral Characteristics and Correction of Long-Term Eddy-Covariance Measurements Over Two Mixed Hardwood Forests in Non-Flat Terrain

We present turbulence spectra and cospectra derived from long-term eddy-covariancemeasurements (nearly 40,000 hourly data over three to four years) and the transferfunctions of closed-path infrared gas analyzers over two mixed hardwood forests inthe mid-western U.S.A. The measurement heights ranged from 1.3 to 2.1 times themean tree height, and peak vegetation area index (VAI) was 3.5 to 4.7; the topographyat both sites deviates from ideal flat terrain. The analysis follows the approach ofKaimal et al. (Quart. J. Roy. Meteorol. Soc. 98, 563–589, 1972) whose results were based upon 15 hours of measurements atthree heights in the Kansas experiment over flatter and smoother terrain. Both thespectral and cospectral constants and stability functions for normalizing and collapsingspectra and cospectra in the inertial subrange were found to be different from those ofKaimal et al. In unstable conditions, we found that an appropriate stabilityfunction for the non-dimensional dissipation of turbulent kinetic energy is of the form ₍₎ = (1 - b^-)^-1/4 - c^-, where representsthe non-dimensional stability parameter. In stable conditions, a non-linear functionG_xy() = 1 + b_xy^c _xy (c_xy < 1) was found to benecessary to collapse cospectra in the inertial subrange. The empirical cospectralmodels of Kaimal et al. were modified to fit the somewhat more (neutraland unstable) or less (stable) sharply peaked scalar cospectra observed over forestsusing the appropriate cospectral constants and non-linear stability functions. Theempirical coefficients in the stability functions and in the cospectral models varywith measurement height and seasonal changes in VAI. The seasonal differencesare generally larger at the Morgan Monroe State Forest site (greater peak VAI) andcloser to the canopy.The characteristics of transfer functions of the closed-path infrared gas analysersthrough long-tubes for CO₂ and water vapour fluxes were studied empirically. This was done by fitting the ratio between normalized cospectra of CO₂ or watervapour fluxes and those of sensible heat to the transfer function of a first-order sensor.The characteristic time constant for CO₂ is much smaller than that for water vapour. The time constant for water vapour increases greatly with aging tubes. Three methods were used to estimate the flux attenuations and corrections; from June through August, the attenuations of CO₂ fluxes are about 3–4% during the daytime and 6–10% at night on average. For the daytime latent heat flux (Q_E), the attenuations are foundto vary from less than 10% for newer tubes to over 20% for aged tubes. Correctionsto Q_E led to increases in the ratio (Q_H + Q_E)/(Q^* - Q_G) by about 0.05 to0.19 (Q_H is sensible heat flux, Q^* is net radiation and Q_G is soil heat flux),and thus are expected to have an important impact on the assessment of energy balanceclosure.     

一个有限差分大气环流模式的并行效率

适合于大规模并行处理机（ＭＰＰ）的中国科学院大气物理研究所新一代九层大气环流模式（ＩＡＰ ９Ｌ ＡＧＣＭ）被发展成功。本文描述了并行模式的设计原则，给出了该模式的并行版本在我国几种主要的并行机上的并行执行效率数据。采用消息传递方式（ＭＰＩ）的区域分解方法用于实现程序的并行，结果表明，对于该模式，只有一维区域分解方法才能取得较好的并行可扩展性，即随着处理器的增加而能得到相应较好的并行加速比和并行效率。     

区域气候模式中侧边界地形缓冲区作用的数值试验

利用区域环境系统集成模式以1998年夏季为例,通过数值试验研究了侧边界地形缓冲区对区域气候模拟结果的影响。结果表明,设置侧边界地形缓冲区使模拟场和观测的均方根误差都有所减小,区域平均中尺度动能也减小,从而有利于维持区域模式长期积分的稳定性。设置地形缓冲区并不改变模拟降水分布格局,但对主要降水区各降水时段降水强度的模拟有不同程度的改进。     

Effects of Spectral Nudging on Oceanic States in a Coarse-Resolution Model

A 1-degree global model is used to investigate the skill of spectral nudging at coarse resolution by performing two numerical experiments, one with spectral nudging and the other without. In the spectral nudging experiment, the model temperature and salinity are nudged to an observed climatological monthly-mean field. The study compares the model mean state, as well as the interannual and decadal variability of oceanic quantities with observations, (e.g., sea surface height (SSH) and sea surface temperature (SST)). Spectral nudging is found to be effective in constraining model drift from the observed mean state of temperature and salinity in the global ocean, which has been reported in previous studies. The present study further shows that spectral nudging significantly improves the model skill of topostrophy (a measure of currents flowing along the topography) in water depth below 2000?m with no clear improvement elsewhere. Despite its known ability to damp oceanic variability at various time scales, spectral nudging can still represent the interannual and decadal variability of SSH and SST well, to a degree comparable to the other experiment.     

Large-Eddy Simulation of Atmospheric Boundary-Layer Flow Over Fluvial-Like Landscapes Using a Dynamic Roughness Model

A recently developed dynamic surface roughness model (Anderson and Meneveau, J Fluid Mech 679:288–314, 2011) for large-eddy simulation (LES) of atmospheric boundary-layer flow over multi-scale topographies is applied to boundary-layer flow over several types of fluvial-like landscapes. The landscapes are generated numerically with simulation of a modified Kardar–Parisi–Zhang equation (Passalacqua et al., Water Resour Res 42:WOD611, 2006). These surfaces possess the fractal-like channel network and anisotropic features often found in real terrains. The dynamic model is shown to lead to accurate flow predictions when the surface-height distributions exhibit power-law scaling (scale invariance) in the prevalent mean flow direction. In those cases, the LES provide accurate predictions (invariant to resolution) of mean velocity profiles. Conversely, some resolution dependence is found for applications in which the landscape’s streamwise spectra do not exhibit pure power-law scaling near wavenumbers corresponding to the LES grid resolution.     

湖北省地形特征对雷电参数的影响

根据2007—2019年湖北省地闪监测数据、地表覆盖数据以及数字地形高程数据,采用数理统计方法,分析研究了地表覆盖类型、海拔、坡向、坡度4类地形信息数据与闪电频数、正闪比例、地闪回击强度、地闪回击陡度等雷电参数的关系。结果表明:(1)建成区单位面积闪电频数和小幅值地闪比例较高,分别为87.1次/km^(2)和8.6%;(2)耕地区和水域区地闪回击陡度值较高;(3)湖北省地闪回击主要集中在海拔0~700 m、坡度0.5°~35°,地闪回击密度随海拔增加呈减少趋势,地闪回击强度平均值随海拔增加呈现先增加后减少的趋势,地闪回击陡度平均值随海拔增加而减少;(4)坡向朝南或朝东地闪回击密度较高,坡向朝南或朝北地闪回击强度平均值较高;(5)不同坡度区间和海拔区间雷电参数的变化规律较为一致。     

四川地形扰动对降水分布影响

引入一维加权平均的谱分析方法定量研究四川地形强迫对该区域降水分布的影响。结果表明:纬向地形和冬季降水谱峰锁相于同一波长（475.8 km）,呈共振关系,地形与其他季节降水呈漂移关系,这与经向和纬向上环流变动有关,即冬季纬向环流占主导,纬向地形触发的大气波动对冬季降水策动作用大;夏季降水是各种不同尺度系统相互作用的结果,地形是重要因素之一。经向和纬向地形特征尺度分别为296.8 km和475.8 km,反映了地形强迫的中尺度特征,且纬向地形谱峰比经向大1个数量级,纬向强迫更明显。夏季降水谱峰比冬季大2个数量级,降水系统纬向特征尺度比冬季小约150 km,说明夏季在纬向地形强迫下,降水系统尺度减小的同时其强度大大增加,这在一定程度上可以解释中尺度对流性降水在夏季偏多。四川夏季最大降水位于雅安地区,其地形扰动比四川整体扰动更明显,故产生的降水也更大。夏季降水和经向地形锁相于同一波长（37.1 km）,经向地形对雅安夏季强降水起关键作用。     

发展方程计算稳定性的比较分析

对线性与非线性发展方程的几种差分格式,以一维线性和非线性平流方程为例,对线性与非线性发展方程差分格式的计算稳定性进行了比较分析,揭示了差分格式结构和初值形式与非线性计算稳定性的关系。     

Stably Stratified Boundary-Layer Flow Over Low Hills: A Comparison of Model Results and Field Data

A linear theory of airflow over low hills andthe MSFD–STAB model are comparedwith published observations made at Cooper'sRidge, north–west of Goulburn inNew South Wales, Australia. The MSFD–STABmodel results show good agreementwith field data in weak to moderate stablystratified flow cases. The originallinear theory overpredicts the wind speed-upratio even in weakly stablecases. After some modification, the agreementis greatly improved.     

Optimization of Upstream Profiles in Modelled Flow Over Complex Terrain

This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.     

“7.15”宜昌大暴雨的地形影响特征

湖北地形复杂,受地形影响产生的局地强降水造成的灾害很多。2010年7月15日宜昌地区的强降水是在满足大尺度降水条件下受地形作用加强的局地暴雨。本文利用实况高空、地面、加密自动站、雷达以及LAPS再分析资料对背景场、低层流场、雷达回波等几个方面进行了分析。得出：冷空气南下和副热带高压的加强促进低层风场的调整,在大尺度降水条件下,地形对低层风场的辐合作用触发了局地强降水的发生。此次过程中地形对降水的触发主要有两方面的作用：一方面为地形迎风坡抬升触发作用,另一方面为地形对近地层流场的影响造成的辐合触发。     

非线性迭代时间积分方案在IAP AGCM-II中的实施及其模拟结果

将非线性迭代时间积分方案用于大气环流模式 IAP AGCM-II中, 并用 Rossby-Hauwitz波进行了检验, 结果表明, 该方案能长期稳定的积分。此外, 还将采用该方案的模式(其中物理过程不变)和原 IAP AGCM-II模式作了模拟结果的比较, 以考察其效果; 发现采用该方案模式对降水尤其是东亚地区降水的模拟能力较原模式有所提高。     

南半球500hPa月平均场谱物理量的气候特征

本文利用南半球500hPa多年（1972年5月至1987年12月）逐月高度场资料,在地转近似条件下,计算了500hPa高度场的各谱物理量,包括月平均环流场、纬向动能谱、经向动能谱、扰动角动量输送以及扰动动能向平均动能转换率。分析了它们的气候特征,发现谱物理量前3波的分布基本上决定了其总体分布,优势域的位置和强度随季节的变化较显著。     

在均匀经纬度网格中计算高纬纬向差商的问题

针对在高纬地区计算纬向差商时,常须把时间步长取得很小,或采取Fourier滤波的问题,提出了可以避免这些困难的差分格式。这种格式可以根据任何给定的、稳定的差分格式构造。在形式上它和后者相同,但要求调整网格距,使之可以取和中纬地区同样大的时间步长,并满足计算稳定性条件。在计算中它很有用,也很方便,几乎没有额外的运算量。而且,如果所根据的格式具有质量或能量守恒性质,新格式也同样具有。     

Performance Evaluation of the GRAPES Model in Wind Simulations Over South China

In the present study, the performance of the GRAPES model in wind simulation over south China was assessed. The simulations were evaluated by using surface observations and two sounding stations in south China. The results show that the GRAPES model could provide a reliable simulation of the distribution and diurnal variation of the wind. It showed a generally overestimated southerly wind speed especially over the Pearl River Delta region and the south of Jiangxi Province as well as the coastal region over south China. GRAPES also exhibited a large number of stations with the opposite surface wind directions over the east of Guangxi and the south of Jiangxi during the nocturnal- to-morning period, as well as an overall overestimation of surface wind over the coastal regions during the afternoon. Although GRAPES could simulate the general evolutional characteristics of vertical wind profile, it underestimated wind speed above 900 hPa and overestimated wind speed below 900 hPa. Though the parameterization scheme of gravity wave drag proved to be an effective method to alleviate the systematic deviation of wind simulation, GRAPES still exhibited large errors in wind simulation, especially in the lower and upper troposphere.