大气数值预报模式中高度地形追随坐标的设计研究：理论分析与理想试验

采用一种统一的地形追随坐标的形式,对Gal-Chen & Somerville(简称Gal.C.S坐标)、平缓坐标(smoothed level vertical coordinate,简称SLEVE坐标)等几种典型的高度地形追随坐标进行了气压梯度力计算误差影响和二维质量平流试验的理论分析,并与一种新提出的高度地形追随坐标——三角函数平缓坐标(简称COS坐标)进行比较.气压梯度力计算误差分析结果显示,与Gal.C.S坐标相比,单尺度平缓坐标(简称SLEVE1坐标)、双尺度平缓坐标(简称SLEVE2坐标)和COS坐标在减小气压梯度力计算误差上有不同程度的改进,SLEVE2坐标和COS坐标比其他两种坐标更具优势,衰减系数b和坐标转换的雅可比项对减小误差起决定性作用.二维质量平流试验也有类似的结果,与无地形的参考试验结果相比,COS坐标的质量输送计算误差最小,且经优化的COS坐标的质量输送计算误差几乎和参考计算误差完全重合,在4种坐标中最优.     

一种改进的平缓-混合地形追随坐标在GRAPES中尺度模式中的应用研究

平缓-混合地形追随坐标（T-F坐标）可以减小坐标面上的地形影响带来的各种计算误差。以余弦三角函数为基函数的平缓-混合坐标（COS坐标）高层坐标面水平,计算误差较小,但是低层坐标面之间的厚度较薄,增大了计算误差,给模式稳定性及模拟效果带来较大的影响。设计一种改进的COS坐标,使低层坐标面垂直分布更加均匀,应用于GRAPES-Meso模式进行理想试验和实际模拟试验。结果表明,改进的COS坐标相对COS坐标,中高层计算误差相当,低层地形作用衰减的垂直变化更加均匀,减小了计算误差,提高了计算稳定性;地形重力波试验结果显示,改进的COS坐标重力波破碎相对COS坐标有一定缓解,更接近解析值;批量模拟试验结果显示,改进的COS坐标各个层次上的月平均模拟偏差比单尺度双曲函数平缓-混合坐标（简称SLEVE1坐标）更小,均方根误差减小,距平相关系数增大。改进的COS坐标有效地解决了COS坐标的计算问题,提高了模式预报效果。      

一种新型高度地形追随坐标在GRAPES区域模式中的实现:理想试验与比较研究

将一种新的高度地形追随坐标(Klemp坐标)引入了GRAPES区域模式,并与传统追随坐标(Gal-Chen坐标)和平缓地形追随坐标(SLEVE,Smooth Level Vertical coordinate)进行了比较。对不同坐标下气压梯度力的计算误差通过理想静止大气试验进行了评估,结果表明:与Gal-Chen坐标和SLEVE坐标相比,Klemp坐标有效地减小了气压梯度力的计算误差。理想重力波模拟试验表明,Klemp坐标下对重力波的模拟相比其他两种坐标也更接近于解析解。模式进一步采用了Mahrer气压梯度计算方案减少了计算误差,并提高了模式的精度和稳定性。实际个例试验与理想试验的结论相似。     

数值预报模式中双三次曲面地形与水平气压梯度力计算

利用C3连续双三次曲面拟合了全球数值模式地形曲面;讨论构建了有复杂地形数值模式引入地形追随高度坐标((z)坐标)后,同时引入包含定常斜率、曲率和挠率的双三次曲面地形,又进一步讨论了双三次曲面地形模式大气的水平气压梯度力计算问题.结果表明,对(z)坐标模式大气的压、温、湿场,通过做经、纬向三次样条拟合,求得地形斜率“静力平衡”气压差,从而插值(反演)任一水平面(海平面)上的气压场,同时可以求得时变的参考大气,则计算水平气压梯度(力)的精度,完全依赖于插值(反演)对应的水平面(海平面)气压场的计算精度.并指出,理论上可按三次样条的曲率判断,做变量场(地形)的局域或单点平滑.     

不同平缓-混合坐标对一次高原准静止锋降水过程模拟影响对比分析

GRAPES_Meso模式预报存在南风偏大、虚假降水偏多等问题,且在大地形下游地区异常明显。平缓-混合坐标可以有效减小气压梯度力计算误差以及平流输送误差,而这两种误差与风场和水汽场预报密切相关。基于GRAPES_Meso模式选择四种平缓-混合坐标对一次典型的高原东部准静止锋降水过程进行模拟分析。模拟结果表明,较弱的天气形势演变下,平缓-混合坐标的改进效果比较明显,可以有效缓解风场预报偏差、虚假降水、虚假天气系统等问题,个例模拟的结果与实况更接近。     

非静力中尺度高分辨率模式模拟中的垂直坐标影响研究

近几年来,随着高性能计算机技术的高速发展,甚高分辨率的中尺度数值预报模式业务应用已成为可能,与之相关的一系列模式技术的新问题也随之提了出来,垂直坐标系的影响就是其中之一。文中借助于美国新一代数值预报模式WRF(WeatherReseachandForecast),比较了非静力中尺度模式高分辨率模拟应用的垂直坐标影响问题。研究表明,当选用几何高度(z)和气压(p)来构造地形追随坐标时,低层两坐标引起的误差基本一样,中高层高度地形追随坐标引起的误差小于气压地形追随坐标;而且分辨率越高差异愈大。高分辨率模拟结果也表明,这种差异趋势是存在的;此外,本文对天气过程的预报要素场进行了相关的分析。     

地形追随坐标系下埃特尔位涡的计算

随着数值天气预报模式的日益普及以及位涡在不同天气、气候时、空尺度上的广泛应用,对采用模式输出资料计算的埃特尔（Ertel）位涡的精度提出了更高要求。将气象要素从中尺度模式普遍采用的地形追随坐标系向定义位涡的z或p坐标系的插值过程（插值法）,是导致埃特尔位涡计算误差的主要原因。文中利用地形追随坐标系与z（p）坐标系的变换关系,导出了z坐标系定义的埃特尔位涡在地形追随坐标系中的表达形式,实现了直接在模式格点上计算埃特尔位涡（直接法）,在提高计算精度的同时,保留了其在z坐标系中的物理意义,方便了分析与应用。为了进一步分析直接法在减少计算误差方面的作用,使用WRF模式对2016年7月发生在华北地区复杂地形背景下的一次黄淮气旋爆发过程进行了模拟,利用模式输出资料对此次天气过程中强降水时段插值法计算误差的分布进行了分析。结果表明,利用直接法计算位涡可有效减少插值法引入的计算误差,中低层均方根误差可达0.5 PVU,中高层可达0.3 PVU。将其应用到中小尺度对流性天气精细化结构的分析及气候统计研究中,可以有效减少插值法误差对结果产生的不良影响,提高计算和分析准确度。      

一个基于Eta垂直坐标的新WRF动力框架及其数值试验

针对我国陡峭地形数值预报难题,本文在国际先进WRF(Weather Research and Forecasting)模式动力框架中引入阶梯地形垂直坐标,以期为改进复杂地形区域数值天气预报提供模式发展可选方案.设计气柱质量变换方法,实现阶梯地形和追随地形两种垂直坐标下动力方程组的形式一致,从而简化方程组离散及程序实现的...     

Pressure gradient errors in a covariant method of implementing the σ-coordinate: idealized experiments and geometric analysis

本文针对经典σ坐标的气压梯度误差(PGF误差),采用多种地形展开理想试验,对比经典σ坐标的经典方案和协变方案的PGF误差。结果表明:计算空间中,协变方案始终能减小经典方案的误差,地形越陡,效果越明显。然而,几何分析和理想试验均表明:协变方案仅能减小计算空间的误差,不能减小物理空间的误差;相比经典方案,正交地形追随坐标能同时减小计算空间和物理空间的误差。     

分块地形坐标大气环流模式框架的计算稳定性及数值试验

用数值试验方法讨论了分块地形η坐标大气环流模式动力框架的计算稳定性; 采用实际初始场, 该框架能够长期稳定地积分.在积分前100d内总有效能量变化很小, 总质量守恒、各物理量场的数值和波形都得到很好保持, 这表明该框架稳定性好.在分析了该框架的积分结果后可知其合乎物理事实, 与IAP 2.0的动力框架相比, 能更好地反映出大地形的环流形势, 这说明有大地形时采用分块地形η坐标要较σ坐标好.另外, 还对青藏高原等大地形的纯机械强迫作了数值试验, 得到了地形强迫平均槽脊的分布特征, 给出了地形纯机械强迫造成的准定常行星波的图像. 这些工作表明采用分块地形η坐标处理大地形是合理可行的.     

Effects of terrain-following vertical coordinates on high-resolution NWP simulations

With increasing resolution in numerical weather prediction(NWP)models,the model topography can be described with finer resolution and includes steeper slopes.Consequently,negative effects of the traditional terrain-following vertical coordinate on high-resolution numerical simulations become more distinct due to larger errors in the pressure gradient force(PGF)calculation and associated distortions of the gravity wave along the coordinate surface.A series of numerical experiments have been conducted in this study,including idealized test cases of gravity wave simulation over a complex mountain,error analysis of the PGP estimation over a real topography,and a suite of real-data test cases.The GRAPES-Meso model is utilized with four different coordinates,i.e.,the traditional terrain-following vertical coordinate proposed by Gal-Chen and Somerville(hereinafter referred to as the Gal.C.S coordinate),the one-scale smoothed level(SLEVE1),the two-scale smoothed level(SLEVE2),and the COSINE(COS)coordinates.The results of the gravity wave simulation indicate that the GRAPES-Meso model generally can reproduce the mountain-induced gravity waves,which are consistent with the analytic solution.However,the shapes,vertical structures,and intensities of the waves are better simulated with the SLEVE2 coordinate than with the other three coordinates.The model with the COS coordinate also performs well,except at lower levels where it is not as effective as the SLEVE2 coordinate in suppressing the PGF errors.In contrast,the gravity waves simulated in both the Gal.C.S and SLEVE1 coordinates are relatively distorted.The estimated PGF errors in a rest atmosphere over the real complex topography are much smaller(even disappear at the middle and upper levels)in the GRAPES-Meso model using the SLEVE2 and COS coordinates than those using the Gal.C.S and SLEVE1 coordinates.The results of the real-data test cases conducted over a one-month period suggest that the three modified vertical coordinates(SLEVE1,SLEVE2,and COS coordinates)give better results than the traditional Gal.C.S coordinate in terms of forecasting bias and root mean square error,and forecasting anomaly correlation coefficients.In conclusion,the SLEVE2 coordinate is proved to be the best option for the GRAPES-Meso model.     

Large-Eddy Simulations of Atmospheric Flows Over Complex Terrain Using the Immersed-Boundary Method in the Weather Research and Forecasting Model

Atmospheric flow over complex terrain, particularly recirculation flows, greatly influences wind-turbine siting, forest-fire behaviour, and trace-gas and pollutant dispersion. However, there is a large uncertainty in the simulation of flow over complex topography, which is attributable to the type of turbulence model, the subgrid-scale (SGS) turbulence parametrization, terrain-following coordinates, and numerical errors in finite-difference methods. Here, we upgrade the large-eddy simulation module within the Weather Research and Forecasting model by incorporating the immersed-boundary method into the module to improve simulations of the flow and recirculation over complex terrain. Simulations over the Bolund Hill indicate improved mean absolute speed-up errors with respect to previous studies, as well an improved simulation of the recirculation zone behind the escarpment of the hill. With regard to the SGS parametrization, the Lagrangian-averaged scale-dependent Smagorinsky model performs better than the classic Smagorinsky model in reproducing both velocity and turbulent kinetic energy. A finer grid resolution also improves the strength of the recirculation in flow simulations, with a higher horizontal grid resolution improving simulations just behind the escarpment, and a higher vertical grid resolution improving results on the lee side of the hill. Our modelling approach has broad applications for the simulation of atmospheric flows over complex topography.     

Evaluation of the Global and Regional Assimilation and Prediction System for Predicting Sea Fog over the South China Sea

In the South China Sea, sea fog brings severe disasters every year, but forecasters have yet to implement an effective seafog forecast. To address this issue, we test a liquid-water-content-only(LWC-only) operational sea-fog prediction method based on a regional mesoscale numerical model with a horizontal resolution of about 3 km, the Global and Regional Assimilation and Prediction System(GRAPES), hereafter GRAPES-3 km. GRAPES-3 km models the LWC over the sea, from which we infer the visibility that is then used to identify fog. We test the GRAPES-3 km here against measurements in 2016 and 2017 from coastal-station observations, as well as from buoy data, data from the Integrated Observation Platform for Marine Meteorology, and retrieved fog and cloud patterns from Himawari-8 satellite data. For two cases that we examine in detail, the forecast region of sea fog overlaps well with the multi-observational data within 72 h. Considering forecasting for0–24 h, GRAPES-3 km has a 2-year-average equitable threat score(ETS) of 0.20 and a Heidke skill score(HSS) of 0.335,which is about 5.6%(ETS) and 6.4%(HSS) better than our previous method(GRAPES-MOS). Moreover, the stations near the particularly foggy region around the Leizhou Peninsula have relatively high forecast scores compared to other sea areas.Overall, the results show that GRAPES-3 km can roughly predict the formation, evolution, and dissipation of sea fog on the southern China coast.     

Model evaluation experiments in the North Atlantic Basin: simulations in nonlinear terrain-following coordinates

A primitive equation ocean circulation model in nonlinear terrain-following coordinates is applied to a decadal-length simulation of the circulation in the North Atlantic Ocean. In addition to the stretched sigma coordinate, novel features of the model include the utilization of a weakly dissipative, third-order scheme for tracer advection, and a conservative and constancy-preserving time-stepping algorithm. The objectives of the study are to assess the quality of the new terrain-following model in the limit of realistic basin-scale simulations, and to compare the results obtained with it against those of other North Atlantic models used in recent multi-model comparison studies.The new model is able to reproduce many features of both the wind-driven and thermohaline circulation, and to do so within error bounds comparable with prior model simulations (e.g., CME and DYNAMO). Quantitative comparison with comparable results obtained with the Miami Isopycnic Coordinate Model (MICOM) show our terrain-following solutions are of similar overall quality when viewed against known measures of merit including meridional overturning and heat flux, Florida Straits and Gulf Stream transport, seasonal cycling of temperature and salinity, and upper ocean currents and tracer fields in the eastern North Atlantic Basin. Sensitivity studies confirm that the nonlinear vertical coordinate contributes significantly to model fidelity, and that the global inventories and spatial structure of the tracer fields are affected in important ways by the choice of lateral advection scheme.     

Numerical simulations of stratiform boundary-layer clouds on the meso-γ-scale. Part I: The influence of terrain height differences

A higher order closure mesoscale model is used to study the influence of terrain height differences on the meso--scale on stratiform boundary-layer clouds. The model is hydrostatic, has a terrain-following coordinate system and a sub-grid scale condensation scheme. It also has a radiation parameterisation for shortwave and longwave radiation in order to calculate radiative cooling/heating. The simulations show that the cloud base height variations induced by the terrain can be much larger than motivated by terrain height variations alone. It is also shown how this behavior is dependent on upstream boundary-layer conditions and/or changes in the turbulence field. Other features studied include the wave in the lee of a ridge/hill and the associated lifting of the cloud base. The results are compared with some simpler physical models, and limitations in those models are demonstrated.     

复杂地形上稳定边界层二维流场的数值模拟

建立了一个二维局地尺度数值模式，对水平扩散采用了水平平滑处理方法，提出了一种在地形追随坐标系中具有较高计算精度的差分近似方法，此方法兼具了地形追随坐标系和直角坐标系的优点，特别适用于地形特征尺度较小和垂直分层要求较细的情况。数值试验表明：该模式具有模拟局地尺度环流的能力，且计算机内存和ＣＰＵ时间都较节省，算法稳定。利用所建立的模式对经过重庆市中区的某一实际地形下的流场进行了数值模拟。     

对流云降水过程中地形作用的数值模拟

该文利用地形追随坐标三维云模式，对湖北的一次对流降水过程进行了平坦地面和理想斜坡地形的模拟试验。结果发现平坦地面模拟的回波强度、最大降水强度和累积降水都与实况有相当的差距。理想地形的对流过程虽然不能模拟出对流群中多单体的生消过程，但最大降水强度和累积降水都比平坦地面的大，更接近实测结果。模拟试验表明，在地形坡度较大时，会产生较强的上升气流，从而使系统对流发展旺盛，产生较大的降水和较强的回波。所以在模拟山区和丘陵地带的对流系统时，采用地形坐标并考虑地形将会改善模式的模拟结果。