渤海海冰和海洋的动力相互作用模拟试验

利用冰-海洋动力耦合模式, 对渤海海冰和海洋之间的相互作用进行了研究和模拟试验.耦合模式中的海洋模式采用普林斯顿海洋模式(POM), 冰模式采用国家海洋局海洋环境预报中心(NMEFC)的渤海海冰业务化数值预报模式.通过模拟研究, 分析了M2分潮对海冰质点运动轨迹、冰外缘线位置、瞬时冰速和冰厚分布的影响; 同时探讨了海冰对冰下海洋的动力作用.     

冻土过程参数化方案与中尺度大气模式的耦合

将发展的冻土过程参数化方案与耦合了NCARLSM陆面过程模式的中尺度大气模式MM5相耦合,对包括我国北方,蒙古国,东西伯利亚的高纬度地区进行了模拟.模拟时间选为北半球季节冻土生成的10月份.模拟结果表明,模式能够较好地模拟出该区域的冻土分布,并且对大气模式中海平面气压场、气温场的模拟改进显著.     

典型中尺度数值预报模式参数化方案的综述与展望

中尺度天气现象对人类的生产生活有重要影响,中尺度数值预报模式是进行数值天气预报的主要工具。受模式分辨率和对天气现象物理机制认识不足的限制,许多复杂天气过程只能用参数化方案来隐式表达,对参数化方案进行研究有助于推动中尺度数值预报模式的模拟和预报效果的不断优化。在探讨国内外典型中尺度数值预报模式特性的基础上,综述了积云对流参数化方案、云微物理参数化方案、边界层参数化方案、陆面过程参数化方案和辐射传输过程参数化方案的研究进展,比较了代表性参数化方案的理论基础和应用场景,发现对独立或组合参数化方案效果的对比试验和改进试验是该领域主要的研究范式。认为参数化方案未来将走向深入探讨各种天气现象的影响因素和物理机制,多要素、多尺度和多方案耦合,更加重视对“灰色区域”的模拟,应用选择更加多元化,并与机器学习技术融合的发展模式。     

青藏高原微地形对极值风速的影响及其参数化模拟方法研究

利用青藏高原铁路沿线1961-2006年7个气象台站和2003年9月-2004年9月7个野外观测点风资料, 结合地理信息系统分析了青藏高原腹地微地形对极值风速的影响. 通过地形因子的参数化处理, 建立了极值风速随海拔和地形参数变化的拟合模型以精确推算复杂地形的极值风速, 并利用临时观测点风资料对其进行了检验. 结果表明, 利用地理信息系统和地形参数化处理方法研究青藏高原微地形对极值风速的影响具有可操作性.     

基于光滑质点流体动力学的海冰热力-动力数值模式

考虑海冰热力因素对其厚度、密集度的影响,在光滑质点流体动力学(SPH)基础上发展了一个海冰热力-动力数值模式。该模式既解决了传统欧拉有限差分法和质点网格法存在的数值扩散问题,同时弥补了光滑质点动力学海冰动力模式未考虑热力因素的不足,具有精确模拟冰缘线运动、计算精度高等优点。首先介绍了光滑质点流体动力学的基本原理,并对海冰生消的热力因素进行了分析,将影响冰厚和密集度的热力因素引入到光滑质点流体动力学的海冰动力模式中,得到该热力-动力模式的控制方程。应用该数值模式对渤海海冰进行了48 h数值模拟,得到了海冰厚度和速度矢量的分布规律;对JZ20-2海域的海冰厚度、冰内温度场分布以及热力因素的变化特性进行了讨论。数值模拟结果表明,该数值模式能够很好地适用于渤海海冰数值模拟,是一种有效的海冰数值模拟方法。     

全球冰-海洋耦合模式的海冰模拟

海冰是全球气候系统的重要分量 ,与大气和海洋的相互作用 ,直接影响大气环流和海洋环流 ,对气候及其变化具有重要影响。文中依据冰、海洋间的热力、动力耦合相互作用 ,改进冰海洋热力耦合方案 ,利用由中国科学院大气物理研究所的 30层海洋模式和基于Flato空化流体流变学的海冰动力模式和Hibler表面热收支平衡的零层海冰热力模式 ,建立全球冰海洋耦合模式。利用大气月平均气候资料 ,利用冰海洋耦合模式对全球海冰的分布及其季节性变化、海冰漂移进行了耦合模拟和分析。模拟的南半球海冰分布及季节变化与实际分析资料非常接近 ,比 2 0层冰海洋耦合模式的结果有显著改进。北半球海冰范围偏小 ,但季节变化的量值与实际相当一致。模拟的海冰速度场反映了南、北半球海冰漂移的主要特征 ,如北极的穿极漂流和南大洋的绕极环流等。对海冰密集度的分析表明 ,模拟结果得以改进原因在于改进的冰海洋热力耦合方案增强了融冰期冰海洋耦合系统海洋热通量增加—密集度减小—能量收支增加的正反馈机制。     

全球冰—海洋耦合模式的海冰模拟

海冰是全球气候系统的重要分量,与大气和海洋的相互作用,直接影响大气环流和海洋环流,对气假及其变化及具有重要影响,文中依据冰、海洋间的热力、动力耦合相互作用,改进冰－海洋热力耦合方案,利用由中国科学院大气物理研究所的３０层海洋模式和基于Ｆｌａｔｏ空化流体流变学的海冰动力模式和Ｈｉｂｌｅｒ表面热收支平衡的零层少冰热力模式,建立全球冰－海洋耦合模式,利用大月平均气候资料,模拟的南半球少冰分布及海冰的分布及其季节性变化、海冰漂移进行了耦合模拟和分析,模拟的南半球海冰分布及季节变化与实际分析资料非常接近,比２０层冰－海洋耦合模式的结果有显著改进,北半球海冰范围偏小,但季节变化的量值与实际相当一致,模拟的海冰速度场反映了南、北半球海冰漂移的主要特征,如北极的穿极漂流和南大洋的绕极环流等,对海冰密集度的分析表明,模拟结果得以     

两种数值模式资料的平均海平面气压和地面风速在中国区域的评估

利用中国气象局国家级自动站(2 421个)的观测资料, 分别对2012年的ECMWF(欧洲中期数值预报中心)和JMA(日本气象厅)数值模式资料的平均海平面气压和地面风速在中国地区的适用性进行了对比研究.结果表明:两种数值模式资料均能在一定程度上反映观测资料所具有的时空分布特征, 东部地区的适用性均要高于西部地区. 对于平均海平面气压, 在西南地区JMA比ECMWF资料更接近实际观测; 而在其他地区, 两种数值模式资料都较接近实际观测, 冬季的结果比夏季好. 对于地面风速, 这两种数值模式资料各具优势.在中国东南部地区, JMA相对于ECMWF的地面风速资料更接近实际观测值, 而在中国西部地区, 则相反.就8个时次的年变化而言, ECMWF资料的年变化趋势与观测资料更为一致, 而JMA资料的地面风速大小与观测资料更为接近.     

陆气耦合模型在实时暴雨洪水预报中的应用

采用加拿大区域性中尺度大气模式MC2(Canadian Mesoscale Compressible Community)和新安江模型单向耦合模型系统,对2005年7月4～15日发生在淮河流域的一场暴雨洪水,进行了实时预报.采用王家坝以上流域的实测降水和王家坝断面的实测洪水资料,对MC2预报降水的时空分布和陆气耦合模型预报的洪水过程进行了分析.结果表明,MC2对该场强降水过程具有很好的预报能力,陆气耦合模型有效地增长了洪水预报的预见期,具有很好的应用前景.     

基于Polar WRF模型的南极冰盖气温、风速和气压数值模拟北大核心CSCD

在全球变暖的背景下,南极已成为全球气候变化研究的热点,然而其区域内的观测站点稀疏且缺乏较长的时间序列,限制了人们对南极气候变化机制的分析与理解。Polar WRF作为目前最先进的极地区域气候模型之一,有力弥补了观测资料不足的缺陷,然而模式存在误差,在应用之前有必要对其定量评估。本文利用Polar WRF3.9.1对2004-2013年南极冰盖2m气温、10m风速和地表气压进行了数值模拟,并与28个气象站数据进行了对比分析,结果表明:模式对气温的模拟值在东南极沿岸偏低,在内陆偏高,在南极半岛既存在冷偏差也存在暖偏差,而对风速和气压的模拟整体呈高估。而从均方根误差和平均绝对误差的空间分布来看,模式对气温和气压的模拟结果在东南极沿岸的精度高于内陆和南极半岛,而风速则在内陆的精度要高于沿岸地区。但总体来说模拟效果较为理想,在2004-2013年间气温、风速、气压的模拟值的变化趋势与实测值的变化趋势相同。模式模拟的年平均2m气温和近地面气压在所有站点都通过了α=0.1的显著性检验,季节误差和月误差整体较小,且所有月份的相关系数都分别大于0.90与0.79。模式对10m风速的模拟精度要略低,部分沿岸站点的年平均误差超过了7.5m·s^(-1),但整体而言其在四季和各个月份的相关性均大于0.5且误差小于4.5m·s^(-1)。虽然Polar WRF作为天气模式,但在模拟长时间尺度的气候方面仍然表现较好。     

基于热力学过程的渤海海冰生消模型

为了对渤海海冰生消问题进行模拟,基于三维自由水面垂向分层动网格的Euler-Lagrangian模式,采用VC（Vertex-Centered）方式的非结构化有限体积方法离散三维浅水方程组,模拟渤海的水流运动过程。以此水动力场为背景,引入热力学过程的影响,对渤海海域的冬季结冰过程进行模拟,建立了海冰生消模型。其中热力学参数主要包括大气温度、相对湿度、风场特征、太阳辐射、感热通量系数以及潜热通量系数。以渤海2011/2012年常冰年的结冰过程为例,将各个热力学参数对冰情影响的敏感性进行分析,对模型进行了验证。结果表明,感热与潜热系数是对冰情影响最敏感的一项,即当感热和潜热通量系数仅减小了0.000 2后,其海冰的最大厚度就减小了15 cm。最后,通过连续实测的水动力数据与2009/2010年冬季的典型海冰灾害过程再一次对模型进行了验证,充分说明了模型具有较高的精确性、稳定性与实用性。     

A sensitivity study of the WRF model in offshore wind modeling over the Baltic Sea

Accurate wind modeling is important for wind resources assessment and wind power forecasting.To improve the WRF model configuration for the offshore wind modeling over the Baltic Sea,this study per-formed a sensitivity study of the WRF model to multiple model configurations,including domain setup,grid resolution,sea surface temperature,land surface data,and atmosphere-wave coupling.The simu-lated offshore wind was evaluated against LiDAR observations under different wind directions,atmo-spheric stabilities,and sea status.Generally,the simulated wind profiles matched observations,despite systematic underestimations.Strengthening the forcing from the reanalysis data through reducing the number of nested domains played the largest role in improving wind modeling.Atmosphere-wave cou-pling further improved the simulated wind,especially under the growing and mature sea conditions.Increasing the vertical resolution,and updating the sea surface temperature and the land surface infor-mation only had a slight impact,mainly visible during very stable conditions.Increasing the horizontal resolution also only had a slight impact,most visible during unstable conditions.Our study can help to improve the wind resources assessment and wind power forecasting over the Baltic Sea.     

SWAN浅水波浪模式在渤海的应用研究——Phillips线性增长比例系数的改进

将国际上通用的SWAN波浪数值模式应用到渤海,以检验其适应性。通过对由风引起的波浪的增长和耗散作用源函数的实验研究,发现SWAN中所采用的Phillips线性增长部分的比例系数和波浪成长有较大的关系,此比例系数随摩擦速度变化(原模式取此系数为一常数)。基于模式试验,引入了一个新的Phillips线性增长项公式,并且利用3个风过程对模式的改进进行了验证。3个风过程的计算结果和实测值比较表明,计算结果在3个过程当中都得到了明显的改进。     

渤海层化结构及潮汐锋面季节变化的数值研究

采用三维斜压海流模式(POM)模拟了渤海海温的季节变化,以海表与海底温差ΔT作为判别依据,发现3月份前整个渤海的表、底温差小于0.5℃,说明渤海处于充分混合状态;进入4月份以后,莱州湾、渤海中部以及渤海海峡的局部水域出现超过2℃的表、底温差,意味着垂向层化结构开始形成;层化区域面积随着海表热通量增大的趋势可一直持续到8月,9月后由于日照量逐渐减小,季节性温越层逐渐消失,11月以后渤海又恢复到充分混合的状态.     

Impact of modified physics in limited area model forecasts

A number of physical factors have been introduced to improve limited area model forecasts. The factors include land surface fluxes, shallow convection and radiation. The model including these additional physical factors (modified physics) is run for five cases of monsoon depression which made landfall over the Indian coast, and the results are compared with those of the control run. The forecasts are verified by computing the root mean square and mean errors. The differences in these skill scores between the two model runs are tested for their statistical significance. It is found that the modified physics has a statistically significant effect on the model skill with the maximum impact on the mean sea level pressure and the temperature. Detailed analyses of mean sea level pressure, wind, rainfall and temperature further confirm that the modified physics has maximum impact on mean sea level pressure and temperature and marginal impact on wind and rainfall. Furthermore, analyses of some model parameters related to physics at a grid point for one case of depression were done. The results show that the inclusion of the land surface physics, shallow convection and radiative processes have produced a better precipitation forecast over the grid point.     

渤海古沙漠沉积特征分析

根据世界沙漠的分布特点和沙漠的发育条件,作者认为晚更新世末期渤海气候干冷,雨量稀少,西北风强盛,松散沉积物丰富,具备发育沙漠的条件.文章从以下方面证明渤海古沙漠的存在:(1)渤海约有四分之一区域为残留砂沉积;(2)在全新世盖层下的更新世砂质沉积具有高倾角的前积纹层;(3)渤海东南缘陆地和岛屿有古渤海沙漠衍生黄土分布.综合分析认为末次冰期时,渤海中部和辽东湾古沙漠发育程度好,渤海湾次之,莱州湾未发育沙漠.     

Performance of nested WRF model in typhoon simulations over West Pacific and South China Sea

Forecasting skill of weather research and forecasting (WRF) model in simulating typhoons over the West Pacific and South China Sea with different trajectories has been studied in terms of track direction and intensity. Four distinct types of typhoons are chosen for this study in such a way that one of them turns toward left during its motion and had landfall, while the second took a right turn before landfall. The third typhoon followed almost a straight line path during its course of motion, while the fourth typhoon tracked toward the coast and just before landfall, ceased its motion and travelled in reverse direction. WRF model has been nested in one way with a coarse resolution of 9?km and a fine resolution of 3?km for this study, and the experiments are performed with National Center for Environmental Prediction-Global Forecasting System (NCEP-GFS) analyses and forecast fields. The model has been integrated up to 96?h and the simulation results are compared with observed and analyzed fields. The results show that the WRF model could satisfactorily simulate the typhoons in terms of time and location of landfall, mean sea-level pressure, maximum wind speed, etc. Results also show that the sensitivity of model resolution is less in predicting the track, while the fine-resolution model component predicted slightly better in terms of central pressure drop and maximum wind. In the case of typhoon motion speed, the coarse-resolution component of the model predicted the landfall time ahead of the actual, whereas the finer one produced either very close to the best track or lagging little behind the best track though the difference in forecast between the model components is minimal. The general tendency of track error forecast is that it increases almost linearly up to 48?h of model simulations and then it diverges quickly. The results also show that the salient features of typhoons such as warm central core, radial increase of wind speed, etc. are simulated well by both the coarse and fine domains of the WRF model.     

Computer model of wind, waves, and longshore currents during a coastal storm

A mathematical model has been developed to forecast or hindcast wind, waves, and longshore currents during the passage of a coastal storm. Storm intensity is a function of the barometric pressure gradient which is modeled by rotating an inverted normal curve around the center of an ellipse. The length and orientation of the major and minor axes of the ellipse control the size and shape of the storm. The path of the storm is determined by a sequence of storm positions for the hindcast mode, and by interpolated positions assuming constant speed and direction for the forecast mode. The site location, shoreline orientation, and nearshore bottom slope provide input data for the shore position. The geostrophic wind speed and direction at the shore site are computed from the latitude and barometric pressure gradient. The geostrophic wind is converted into surface wind speed and direction by applying corrections for frictional effects over land and sea. The surface wind speed and direction, effective fetch, and wind duration are used to compute wave period, breaker height, and breaker angle at the shore site. The longshore current velocity is computed as a function of wave period, breaker height and angle, and nearshore slope. The model was tested by comparing observed data for several coastal locations with predicted values for wind speed, wave period and height, and longshore current velocity. Forecasts were made for actual storms and for hypothetical circular and elliptical storms.     

Operational wave prediction of extreme storms in Northern Europe

The disastrous effects of numerous winter storms on the marine environment in the North Sea and the Baltic Sea during the last decade show that wind waves generated by strong winds actually represent natural hazards and require high quality wave forecast systems as warning tools to avoid losses due to the impact of rough seas. Hence, the operational wave forecast system running at the German Weather Service including a regional wave model for the North Sea and the Baltic Sea is checked extensively whether it provides reasonable wave forecasts, especially for periods of extraordinary high sea states during winter storms. For two selected extreme storm events that induced serious damage in the area of interest, comprehensive comparisons between wave measurements and wave model forecast data are accomplished. Spectral data as well as integrated parameters are considered, and the final outcome of the corresponding comparisons and statistical analysis is encouraging. Over and above the capability to provide good short-term forecast results, the regional wave model is able to predict extreme events as severe winter storms connected with extraordinary high waves already about 2 days in advance. Therefore, it represents an appropriate warning tool for offshore activities and coastal environment.