月平均环流的长期数值预报

本文对当前月平均环流的长期数值预报状况做了总结,用大气环流模式做月平均500hPa高度预报,相关系数约0.4。对两类集合预报,即蒙特卡罗预报与落后平均预报的优缺点进行了比较。最后讨论了月平均环流预报的可预报性,指出有可能做6—11个月的预报,但预报准确率因受气候噪声影响而有一定限度。     

大气初始场对短期气候数值预测的影响

利用中国科学院大气物理研究所发展的九层大气环流模式(简称IAP 9L2°×2.5°-AGCM)分析了大气初始场对短期气候数值预测的影响,分别从实际观测海温、海冰的外强迫和气候态海温、海冰的外强迫出发,进行两组集合回报试验,每组包含3个试验,分别将实时NCEP-Ⅱ资料和对NCEP-Ⅱ资料经5 d平滑、11 d平滑后的资料作为大气初始场,进行17 a(1988—2004年)集合回报试验,采用相关分析方法对试验结果进行对比分析。结果表明,比较相关系数定量检验出大气初始场对热带地区可预报性影响较小,而对中、高纬度地区影响很大。通过对6个试验中对应气象要素在对流层各层距平时间相关系数以及17 a空间异常相关系数均值比较分析发现对东亚(中国)地区夏季气候而言,NCEP-Ⅱ资料经5 d平滑后生成的大气初始场对应回报试验结果相对最好。     

我国大气动力学和数值天气预报研究工作的进展

我国自解放以来大气动力学和数值天气预报研究工作蓬勃发展,取得了不少可喜的成果。在祖国各地已初步形成了一支从事动力学和数值预报工作的力量,研究与大气预报实践密切相关的各种动力学问题,尤其是对我国天气影响很大的问题,同时也注意到基础理论问题的研究。我国数值天气预报的理论及数学方法的研究工作也有许多成果,数     

Comparison of prediction performance using statistical postprocessing methods

As the 2018 Winter Olympics are to be held in Pyeongchang, both general weather information on Pyeongchang and specific weather information on this region, which can affect game operation and athletic performance, are required. An ensemble prediction system has been applied to provide more accurate weather information, but it has bias and dispersion due to the limitations and uncertainty of its model. In this study, homogeneous and nonhomogeneous regression models as well as Bayesian model averaging (BMA) were used to reduce the bias and dispersion existing in ensemble prediction and to provide probabilistic forecast. Prior to applying the prediction methods, reliability of the ensemble forecasts was tested by using a rank histogram and a residualquantile-quantile plot to identify the ensemble forecasts and the corresponding verifications. The ensemble forecasts had a consistent positive bias, indicating over-forecasting, and were under-dispersed. To correct such biases, statistical post-processing methods were applied using fixed and sliding windows. The prediction skills of methods were compared by using the mean absolute error, root mean square error, continuous ranked probability score, and continuous ranked probability skill score. Under the fixed window, BMA exhibited better prediction skill than the other methods in most observation station. Under the sliding window, on the other hand, homogeneous and non-homogeneous regression models with positive regression coefficients exhibited better prediction skill than BMA. In particular, the homogeneous regression model with positive regression coefficients exhibited the best prediction skill.     

Estimating the propagation velocity of atmospheric fronts from surface wind observations

Abstract

Surface mesonet winds recorded at 10‐min intervals are used to estimate the propagation velocities of atmospheric fronts in East Coast winter storms during the Canadian Atlantic Storms Program (CASP). The frontal motion is modelled locally as the translation of a line across which there is an abrupt shift in wind direction. The mesonet is used to detect the propagation velocity of the windshift line.

Frontal velocities estimated using mesonet winds for all cases in which fronts passed through the mesonet (two cold fronts and three warm fronts) are in close agreement with those deduced from synoptic charts. Recommendations are given for using the method as a research tool to estimate frontal motions in oceanographic studies of wind‐driven circulation.     

Modification of fronts by the Alps: Simulations and numerical experimentation

Summary An isentropic limited area model is used to simulate and investigate the frontal passages of 3 May and 8 October 1987. It is demonstrated that a southward outbreak of air with high potential vorticity on 3 May 1987 affected the propagation of the front to the north of the Alps and the related formation of an orographic jet. Moreover, the outbreak plays a crucial role in the genesis of a lee cyclone. On the other hand, no such outbreak occurred on 8 October 1987 and it is shown that the propagation of the front near the ground was hardly affected by the flow at upper levels.With 10 Figures     

Effects of mesoscale sea-surface temperature fronts on the marine atmospheric boundary layer

A numerical modelling study is presented focusing on the effects of mesoscale sea-surface temperature (SST) variability on surface fluxes and the marine atmospheric boundary-layer structure. A basic scenario is examined having two regions of SST anomaly with alternating warm/cold or cold/warm water regions. Conditions upstream from the anomaly region have SST values equal to the ambient atmosphere temperature, creating an upstream neutrally stratified boundary layer. Downstream from the anomaly region the SST is also set to the ambient atmosphere value. When the warm anomaly is upstream from the cold anomaly, the downstream boundary layer exhibits a more complex structure because of convective forcing and mixed layer deepening upstream from the cold anomaly. An internal boundary layer forms over the cold anomaly in this case, generating two distinct layers over the downstream region. When the cold anomaly is upstream from the warm anomaly, mixing over the warm anomaly quickly destroys the shallow cold layer, yielding a more uniform downstream boundary-layer vertical structure compared with the warm-to- cold case. Analysis of the momentum budget indicates that turbulent momentum flux divergence dominates the velocity field tendency, with pressure forcing accounting for only about 20% of the changes in momentum. Parameterization of surface fluxes and boundary-layer structure at these scales would be very difficult because of their dependence on subgrid-scale SST spatial order. Simulations of similar flow over smaller scale fronts (<5 km) suggest that small-scale SST variability might be parameterized in mesoscale models by relating the effective heat flux to the strength of the SST variance.     

A numerical investigation of the JASIN atmospheric boundary layer

A numerical model of the cloudy marine boundary layer is described and used to investigate the role of entrainment instability on the developing boundary layer. In general, previous studies have been limited to boundary layers capped by convectively stable inversions or have described only cumulus fields. Here we extend a stratus-capped boundary-layer model to consider the transition to a convectively unstable cloud layer capped by an inversion across which there is a rapid decrease in wet-bulb or equivalent potential temperature. In this case, the inversion is very active and the entrainment rate is determined by the active instability at the interface, in contrast to the mean turbulent motion within the boundary layer.The model is used to interpret the observed boundary layer from the JASIN experiment. Cool, dry air is modified by prolonged passage over increasingly warmer ocean which leads to the development of a convectively unstable cloud layer.     

Generalized inversion of a global numerical weather prediction model

Summary We construct the generalized inverse of a global numerical weather prediction (NWP) model, in order to prepare initial conditions for the model at time t=0 hrs. The inverse finds a weighted, least-squares best-fit to the dynamics for –24<t<0, to the previous initial condition att=–24, and to data att=–24,t=–18,t=–12 andt=0. That is, the inverse is a weak-constraint, four-dimensional variational assimilation scheme. The best-fit is found by solving the nonlinear Euler-Lagrange (EL) equations which determine the local extrema of a penalty functional. The latter is quadratic in the dynamical, initial and data residuals. The EL equations are solved using iterated representer expansions. The technique yields optimal conditioning of the very large minimization problem, which has 10⁹ hydrodynamical and thermodynamical variables defined on a 4-dimensional, space-time grid.In addition to introducing the inverse NWP model, we demonstrate it on a medium-sized problem, namely, a study of the impact of reprocessed cloud track wind observations (RCTWO) from the 1990 Tropical Cyclone Motion Experiment (TCM-90). The impact is assessed in terms of the improvement of forecasts in the South China Sea att=+48 hours. The calculation shows that the computations are manageable, the iteration scheme converges, and that the RCTWO have a beneficial impact.     

A high resolution numerical study of Gulf of Mexico fronts and eddies

Summary The Gulf of Mexico (GOM) circulation is simulated using the DieCAST ocean model, with a horizontal resolution of 1/12° and 20 vertical layers. The results compare well with observations of both large and small scale features, including Loop Current frontal occlusions associated with frontal eddies. The simulation is carried out without any data assimilation. The frontal eddies tend to be spaced at about 90° intervals around the Loop Current, leading to a Loop Current head shaped like a square with rounded corners. The pattern rotates as the eddies circle the Loop, and frontal eddies elongate as they squeeze through the Florida Strait. Major warm core eddies separate regularly from the Loop Current and propagate to the western GOM. Old warm core eddies in the western Gulf dissipate through bottom drag effects, which also generate cyclonic parasitic eddies. Newly arrived warm core eddies merge with old ones in the western GOM. Recently separated elongated Loop Current eddies can rotate and reattach temporarily to the Loop Current. The barotropic flow component develops eddies between the main separated warm core eddy and the Loop Current due to eastward dispersion, as the main eddy itself propagates westward into the Gulf.With 10 Figures     

Surface heat fluxes from a numerical weather prediction system

Twenty-two months (July 1983-April 1985) of surface heat fluxes predicted at day 1 from a numerical weather prediction system have been processed. Monthly means and monthly standard deviations of available surface short-wave, long-wave, latent and sensible heat fluxes as well as annual means have been computed. The global mean of the annual net sea-surface heat flux is about 40 Wm^–2 and is therefore far from equilibrium. When used to force an oceanic model, these fluxes would tend to warm the ocean and would produce an unrealistic transport of heat by the oceanic general circulation. They therefore need to be corrected. This correction appears feasible because the main difference between these fluxes and long-term climatologies appears largely independent of the month and the latitude. This suggests that the errors have a systematic origin. The corrected fluxes allow both the reproduction of a realistic seasonal migration of the zero net heat-flux line and the reproduction of the annual meridional heat transport in the different oceans, within the range of previous estimates.     

论大气边界层的数值模拟方法

本文采用分层模式和连续模式模拟了一维、稳定的大气边界层中一些主要的气象变量随时间的变化规律,讨论了几种气象条件及下垫面条件下计算结果之间的差异,并且还讨论了常值湍流热通量层存在的条件及其高度随时间的变化,从而给出了两种模式适用的范围.     

数值天气预报中的不一致性问题综述

近年来随着数值天气预报的不断发展,数值天气预报中出现了一个新的值得注意的问题:预报不一致性,即连续多次预报中,前后两次预报所做出的预报结果之间差异较大的现象。本文较系统地概述了这一问题,总结了预报不一致性的概念、定量分析方法、产生的可能原因及其与预报误差的关系等方面的研究成果。目前预报不一致性的定量分析方法还较少,其适用性和客观性等还需进一步研究。在预报不一致性产生的原因方面,现有的研究仅限于理论上的简要分析,初始误差和模式误差在预报不一致性形成过程中的具体作用等尚不清楚。最后阐述了在预报不一致性研究中存在的一些问题及可能的研究方向。     

Initialization of cloud water in a numerical weather prediction model

Summary In recent years many studies have shown the importance of treating condensation processes in a consistent manner in numerical weather prediction models. Among emerging improvements is the explicit treatment of cloud water, and in some cases precipitating water. An unresolved problem then is how to initialize the cloud water, especially since this quantity is not treated in the most commonly used analysis schemes.In this study, a method for initializing the cloud water in a numerical weather prediction (NWP) model will be presented and tested. The implications for the model's spin-up are investigated. Information from an earlier run (first guess fields) is used, together with satellite data. If necessary, humidity enhancement is performed where clouds are indicated by those sources. The results indicate that initialization of the cloud water field by itself does not have a large effect on the spin-up of precipitation and clouds. A much larger effect is obtained when the humidity field is enhanced. The spin-up time for precipitation is then reduced from 12 to 6 hours, while for cloud cover it is reduced to only 1–2 hours. The method is computationally very efficient, and is particularly useful over data-sparse areas, such as the oceans.An investigation of the different terms in the cloud water tendency equation is done and the results interpreted in terms of spin-up of cloud parameters. These tests confirm that the cloud water field only accounts for a small part of the spin-up effect. These also show that the production of cloud water per time step increases throughout the simulation.With 9 Figures     

动力初值预报试验

本文提出一种动力初值方案,这一方案考虑了不同纬度大气适应过程的不同特征:在低纬度,主要调整气压场和温度场,而对流场只作适当调整;在中高纬度,则主要调整流场,而对气压场和温度场只作适当调整。此方案计算收敛很快,只需相当于做2至3小时预报的时间。用此初值做48小时预报,得到了较好的预报结果。     

Analysis of results of numerical experiments with an atmospheric general circulation model with allowance for peculiarities of the underlying surface in megalopolises

A hypothesis is put forward that underlying surface peculiarities over large cities should be taken into account in the atmospheric circulation modeling, climate change studies, and in the numerical weather forecast. It is demonstrated based on the analysis of satellite images that the scales of changes in the underlying surface properties amount to hundreds of kilometers and that their influence can be studied with a large-scale atmospheric general circulation model. The influence of changes in the underlying surface characteristics in the cities on the surface temperature is analyzed. It is shown that the model atmosphere is sensitive to changes in surface properties in the cities. It is also shown that the surface air temperature approaches the observed climate over most land part of the urban territories, where peculiarities of the urban territories are taken into account. The model of the city should be changed in the areas, where such a result was not obtained. According to the experimental results, the model of the city in the mid- and high latitudes should differ from the model in the southern latitudes, but it should be taken into consideration for all urbanized areas of the world.     

Retrieval of atmospheric and oceanic parameters and the relevant numerical calculation

It is well known that retrieval of parameters is usually ill-posed and highly nonlinear, so parameter retrieval problems are very difficult. There are still many important theoretical issues under research, although great success has been achieved in data assimilation in meteorology and oceanography. This paper reviews the recent research on parameter retrieval, especially that of the authors. First, some concepts and issues of parameter retrieval are introduced and the state-of-the-art parameter retrieval technology in meteorology and oceanography is reviewed briefly, and then atmospheric and oceanic parameters are retrieved using the variational data assimilation method combined with the regularization techniques in four examples： retrieval of the vertical eddy diffusion coefficient; of the turbulivity of the atmospheric boundary layer; of wind from Doppler radar data, and of the physical process parameters. Model parameter retrieval with global and local observations is also introduced.