共查询到20条相似文献，搜索用时 46 毫秒
1.
《Ocean Modelling》2002,4(3-4):249-267
During the course of developing new numerical algorithms for a terrain-following ocean modeling system (TOMS), different numerical aspects have been evaluated through a comparison between two widely used community ocean models, the Princeton ocean model (POM) and the regional ocean modeling system (ROMS). While both models aim at modeling coastal to basin-scale problems using similar grids, their numerical algorithms, code structure, and parameterization options are very different. Sensitivity studies with an idealized channel flow and a steep seamount configuration demonstrate how different algorithms in the two models may affect numerical errors, the stability of the code and the computational efficiency. For example, new pressure gradient schemes using polynomial fits and new time stepping algorithms may reduce numerical errors and allow using longer time steps than standard schemes do. However, the new schemes may require more careful choices of time steps and the use of higher order advection schemes to maintain numerical stability. 相似文献
2.
Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution 总被引：5，自引：0，他引：5
In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper. 相似文献
3.
A new high-resolution Eulerian numerical method is proposed for modelling quasigeostrophic ocean dynamics in eddying regimes. The method is based on a novel, second-order non-dissipative and low-dispersive conservative advection scheme called CABARET. The properties of the new method are compared with those of several high-resolution Eulerian methods for linear advection and gas dynamics. Then, the CABARET method is applied to the classical model of the double-gyre ocean circulation and its performance is contrasted against that of the common vorticity-preserving Arakawa method. In turbulent regimes, the new method permits credible numerical simulations on much coarser computational grids. 相似文献
4.
渤海AVHRR多通道海冰密集度反演算法试验研究 总被引：2，自引：1，他引：1
为了得到更精确的渤海海冰密集度反演参数,采用辽东湾不同类型海冰ASD实测数据,在分析光谱特征的基础上,针对NOAA/AVHRR数据确定出合适海冰密集度反演算法阈值。继而,基于线性光谱混合模型的多通道反演算法进行了一系列算法试验。同时实现了基于LandSat5-TM数据的渤海海冰密集度场反演,并利用该结果与AVHRR单通道和多通道算法得到的海冰密集度反演结果进行比对分析。定量误差分析结果表明,当单通道算法或组合算法中包含1通道时,与Landsat5-TM反演结果的平均误差为正值,包含2通道且不包含1通道时,平均误差为负值;同时使用这两个通道较只包含其一的各种组合算法的平均误差明显偏小;在各种组合算法中,1245四个通道组合反演的海冰密集度结果误差最小,可应用于渤海AVHRR数据海冰密集度反演。 相似文献
5.
Novel laboratory experiments and numerical modelling have been performed to study the advection scales of suspended sediment in the swash zone. An experiment was designed specifically to measure only the sediment picked up seaward of the swash zone and during bore collapse. The advection scales and settling of this sediment were measured during the uprush along a rigid sediment-free beach face by a sediment trap located at varying cross-shore positions. Measurements were made using a number of repeated solitary broken waves or bores. Approximately 25% of the pre-suspended sediment picked up by the bores reaches the mid-swash zone (50% of the horizontal run-up distance), indicating the importance of the sediment advection in the lower swash zone. The pre-suspended sediment is sourced from a region seaward of the shoreline (still water line) which has a width of about 20% of the run-up distance. An Eulerian–Lagrangian numerical model is used to model the advection scales of the suspended sediment. The model resolves the hydrodynamics by solving the non-linear shallow water equations in an Eulerian framework and then solves the advection–diffusion equation for turbulence and suspended sediment in a Lagrangian framework. The model provides good estimates of the measured mass and distribution of sediment advected up the beach face. The results suggest that the correct modelling of turbulence generation prior to and during bore collapse and the advection of the turbulent kinetic energy into the lower swash is important in resolving the contribution of pre-suspended sediment to the net sediment transport in the swash zone. 相似文献
6.
Assimilation of image sequences in numerical models 总被引：1，自引：0，他引：1
By OLIVIER TITAUD ARTHUR VIDARD INNOCENT SOUOPGUI FRANÇOIS-XAVIER LE DIMET 《地球，A辑:动力气象学与海洋学》2010,62(1):30-47
Understanding and forecasting the evolution of geophysical fluids is a major scientific and societal challenge. Forecasting algorithms should take into account all the available information on the considered dynamic system. The variational data assimilation (VDA) technique combines all these informations in an optimality system (O.S.) in a consistent way to reconstruct the model inputs. VDA is currently used by the major meteorological centres. During the last two decades about 30 satellites were launched to improve the knowledge of the atmosphere and of the oceans. They continuously provide a huge amount of data that are still underused by numerical forecast systems. In particular, the dynamic evolution of certain meteorological or oceanic features (such as eddies, fronts, etc.) that the human vision may easily detect is not optimally taken into account in realistic applications of VDA. Image Assimilation in VDA framework can be performed using 'pseudo-observation' techniques: they provide apparent velocity fields, which are assimilated as classical observations. These measurements are obtained by certain external procedures, which are decoupled with the considered dynamic system. In this paper, we suggest a more consistent approach, which directly incorporates image sequences into the O.S. 相似文献
7.
《Ocean Modelling》2008,20(3):293-311
A diagnostic method for calculating physical and numerical mixing of tracers in ocean models is presented. The physical mixing is defined as the turbulent mean tracer variance decay rate. The numerical mixing due to discretisation errors of tracer advection schemes is shown to be the decay rate between the advected square of the tracer variance and the square of the advected tracer and can be easily implemented into any ocean model. The applicability of the method is demonstrated for four test cases: (i) a one-dimensional linear advection equation with periodic boundary conditions, (ii) a two-dimensional flat-bottom lock exchange test case without mixing, (iii) a two-dimensional marginal sea overflow study with mixing and entrainment and (iv) the DOME test case with a dense bottom current propagating down a broad linear slope. The method has a number of advantages over previously introduced estimates for numerical mixing. 相似文献
8.
Doppler radar radial wind observations are modelled in numerical weather prediction (NWP) within observation errors which consist of instrumental, modelling and representativeness errors. The systematic and random modelling errors can be reduced through a careful design of the observation operator (Part I). The impact of the random instrumental and representativeness errors can be decreased by optimizing the processing of the so-called super-observations (spatial averages of raw measurements; Part II).
The super-observation processing is experimentally optimized in this article by determining the optimal resolution for the super-observations for different NWP model resolutions. A 1-month experiment with the HIRLAM data assimilation and forecasting system is used for radial wind data monitoring and for generating observation minus background (OmB) differences. The OmB statistics indicate that the super-observation processing reduces the standard deviation of the radial wind speed OmB difference, while the mean vector wind OmB difference tends to increase. The optimal parameter settings correspond at a measurement range of 50 km (100 km) to an averaging area of 1.7 km2 (7.3 km2 ).
In conclusion, an accurate and computationally feasible observation operator for the Doppler radar radial wind observations is developed (Part I) and a super-observation processing system is optimized (Part II). 相似文献
The super-observation processing is experimentally optimized in this article by determining the optimal resolution for the super-observations for different NWP model resolutions. A 1-month experiment with the HIRLAM data assimilation and forecasting system is used for radial wind data monitoring and for generating observation minus background (OmB) differences. The OmB statistics indicate that the super-observation processing reduces the standard deviation of the radial wind speed OmB difference, while the mean vector wind OmB difference tends to increase. The optimal parameter settings correspond at a measurement range of 50 km (100 km) to an averaging area of 1.7 km
In conclusion, an accurate and computationally feasible observation operator for the Doppler radar radial wind observations is developed (Part I) and a super-observation processing system is optimized (Part II). 相似文献
9.
水下溢油数值模拟研究 总被引：1，自引：0，他引：1
基于Lagrange积分法和Lagrange粒子追踪法建立了一个水下溢油数值模型。该溢油模型由两个子模型组成:羽流动力模型和对流扩散模型,其中羽流动力模型用以模拟溢油的喷发阶段和浮力羽流阶段;对流扩散模型用以模拟溢油的对流扩散阶段。通过数值实验,结合实验室水槽实验和水下溢油现场实验的观测资料进行模型验证。实验结果表明,模拟结果与观测资料一致性较好,从而验证了本文溢油模型的合理性和准确性;羽流动力模型为对流扩散模型提供源,海流、海水的垂向密度结构和油滴的直径分布是影响溢油在对流扩散阶段运动和分布的主要因素。 相似文献
10.
Surface currents measured by high frequency （HF） radar arrays are assimilated into a regional ocean model over Qingdao coastal waters based on Kalman filter method. A series of numerical experiments are per- formed to evaluate the performance of the data assimilation schemes. In order to optimize the analysis pro- cedure in the traditional ensemble Kalman filter （ENKF）, a different analysis scheme called quasiensemble Kaman filter （QENKF） is proposed. The comparisons between the ENKF and the QENKF suggest that both them can improve the simulated error and the spatial structure. The estimations of the background error covariance （BEC） are also assessed by comparing three different methods： Monte Carlo method; Canadian quick covariance （CQC） method and data uncertainty engine （DUE） method. A significant reduction of the root-mean-square （RMS） errors between model results and the observations shows that the CQC method is able to better reproduce the error statistics for this coastal ocean model and the corresponding external forcing. In addition, the sensibility of the data assimilation system to the ensemble size is also analyzed by means of different scales of the ensemble size used in the experiments. It is found that given the balance of the computational cost and the forecasting accuracy, the ensemble size of 50 will be an appropriate choice in the Qingdao coastal waters. 相似文献
11.
An unrealistically high-salinity maximum is found to develop in a high-resolution model of the north and equatorial Atlantic below the shallow halocline in the Gulf of Guinea. The spurious water mass with salinities too high by as much as 1 psu is formed when the vertical advection is treated by the standard central-differencing advection scheme. The problem is considerably reduced either by increasing the vertical resolution of the numerical grid, or by switching to a higher-order upwind-weighted scheme for vertical advection. This note stresses the need for a careful consideration of vertical discretization even in typical high-resolution ocean general circulation models (OGCMs). Particular attention is needed for studying heat and salt budgets or transports of biogeochemical tracers. 相似文献
12.
一个高分辨率的长江口台风风暴潮数值预报模式及其应用 总被引：13，自引：1，他引：13
利用河口海岸海洋模式(ECOM-Si)建立了一个适用于长江口区风暴潮的数值预报模式.该模式采用对岸线有较好拟合能力的自然正交水平坐标系统和能分辨较复杂海底地形的垂直σ坐标系统.模式考虑了长江口径流量对风暴潮的影响,部分地考虑了天文潮和风暴潮非线性相互作用对风暴增水的影响.风暴潮预报的大气强迫场用模型气压场和模型风场.利用所建立的模式对长江口区台风风暴潮进行了8个个例模拟,模拟增水与实测增水的峰值相比较,平均绝对误差不足10cm.利用本研究建立的模式,就气象因子对风暴潮位的敏感性进行了数值试验.试验结果表明,台风中心气压降低(升高)20hPa可导致约100cm的风暴潮位升高(或降低).台风最大风速半径误差对台风增水的变化影响也较显著.试验还表明,长江径流量增加1倍(减半),可以造成风暴潮的平均增加25cm(减小13cm).天文潮位相变化对风暴增水的影响数值试验表明,当台风暴潮与天文潮在不同位相相互作用,可使风暴潮位最大增加达70cm或减小90cm. 相似文献
13.
14.
An Eta-Coordinate Version of the Princeton Ocean Model 总被引：1，自引：0，他引：1
The Princeton Ocean Model (POM) is widely used in both the research and real-time marine forecast communities. However, there
is a single heavily criticized feature and shortcoming of POM, i.e., the pressure gradient (PG) error associated with sloping
topography. To overcome this problem, in this paper we present an eta-coordinate POM (ECPOM). The eta-coordinate system was
originally adopted in atmospheric numerical models and has proved to be effective in removing the large errors associated
with the horizontal pressure gradient force and advection and diffusion terms along a steeply sloping topography. A familiar
idealized seamount problem is used to test the PG error in the ECPOM in order to isolate errors from other sources, as is
done in the literature. The model is set up so that, ideally, the ocean should be at rest forever; any development of circulation
is regarded as numerically artificial. Results show that in ECPOM, errors are dramatically reduced by two orders of magnitude
over a wide range of parameters. As a trade-off, the bottom boundary of the ECPOM is a stepwise instead of a continuous one,
as in the sigma-coordinate POM (SCPOM). Modification of the SCPOM to the ECPOM is straightforward and minimal. Compared to
the z-coordinate model, the ECPOM is numerically simpler since the eta-coordinate has a surface-following character.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
15.
采用质点跟踪方法对物质输运方程平流项数值格式的改进 总被引：1，自引：1，他引：0
用数值模式对河口海岸地区的物质输运进行计算时,平流项的数值格式必须要能对物质浓度锋面进行正确处理,以避免产生过多的数值耗散或频散。本文中设计了一种在网格内设置一些质点并对质点进行跟踪的格式计算平流项。结果表明,质点跟踪格式在一维情形下无频散和几乎没有耗散,在二维情形下无频散和在水深变化剧烈的地方基本避免了垂向数值耗散。与其他数值格式的耗散性和频散性相比,本文中设计的数值格式明显地提高了物质输运方程中平流项的计算精度,在河口海洋物质输运的计算中具有较大的应用价值。 相似文献
16.
This study provides a practical guide to the use of classical tidal prediction algorithms in coastal numerical forecasting models such as tide and tide-storm-surge models. Understanding tidal prediction parameter formulas and their limitations is key to successfully modifying and upgrading tidal prediction modules in order to increase the accuracy of perpetual interannual simulations and, in particular, storm-surge modeling studies for tide-dominated coastal environments. The algorithms for the fundamental prediction parameters, the five astronomical variables, used in tidal prediction are collated and tested. Comparisons between their estimation using different parameterizations shows that these methods yield essentially the same results for the period 1900–2099, revealing all are applicable for tidal forecasting simulation. Through experiments using a numerical model and a harmonic prediction program, the effects of nodal modulation correction and its update period on prediction accuracy and sensitivity are examined and discussed using a case study of the tidally-dominated coastal regime off the west coast of Korea. Results indicate that this correction needs updating within <30 days for accurate perpetual interannual tidal and mean sea-level predictions, and storm-surge model predictions requiring centimeter accuracy, for tidally-dominated coastal regimes. Otherwise, unacceptable systematic errors occur. 相似文献
17.
We apply the second-order moment (SOM) advection scheme of (Prather, M.J. 1986. Numerical advection by conservation of second-order moments. J. Geophys. Res. 91, 6671–6681.) to the simulation of the large-scale circulation of the Arctic Ocean with a coupled ocean–sea-ice model. Compared to three other advection schemes commonly employed in ocean simulations (centred differences, flux corrected transport, and multidimensional positive definite advection transport), the SOM method helps preserve the vertical structure of Arctic water masses. The depth, thickness and hydrographic properties of the Arctic Surface Water and the Arctic Atlantic Layer are better represented with SOM than with any of the other three advection algorithms. We also present a convenient method for calculating the implicit numerical diffusivity of upstream based schemes, such as the SOM method, and discuss three approaches for improving the monotonicity properties of the SOM algorithm. 相似文献
18.
A complication of finite-volume triangular C-grid methods is the numerical emergence of horizontal divergence errors that lead to grid-scale oscillations in vertical velocity. Nonlinear feedback via advection of momentum can lead to numerical instability in velocity modes via positive feedback with spurious vertical velocities induced by horizontal divergence truncation error. Existing strategies to mitigate divergence errors such as direct divergence averaging and increased diffusion do not completely mitigate horizontal vertical velocity oscillations. We present a novel elliptic filtering approach to mitigate this spurious error and more accurately represent vertical velocities via improved calculation of horizontal divergences. These results are applied to laminar curved channel flows, demonstrating the applicability of the method to reproduce secondary flow features. 相似文献
19.
在悬沙输运的数值模拟中,初始场的准确给定至关重要。目前诸多确定初始场的方案均存在一定的缺陷,初始场的准确性有待进一步提高。本文基于一个三维悬沙输运伴随同化模型,通过孪生实验和实际实验,对模型初始场进行了伴随法反演研究。在孪生实验中,首先验证了初始场的相对重要性;其次,探讨了初始场的反演结果对优化算法、初始猜测值、卫星遥感数据数量、同化时间窗口宽度和背景流场误差的敏感性;最后,比较了伴随法和插值法重构初始场的能力。孪生实验结果表明:最速下降法对初始场的优化反演效果要优于三种共轭梯度法和有限记忆BFGS法;初始场的反演效果对初始猜测值、卫星遥感数据数量和背景流场误差不敏感,而对同化窗口宽度较为敏感;与插值法相比,伴随法是重构模型初始场更有效的手段。实际实验中,在杭州湾海域同化典型的小潮时期和大潮时期的GOCI卫星遥感资料所得表层悬沙浓度数据,优化反演了初始场。实际实验结果表明:数据同化后,得到了更符合实际的最优初始场,表明伴随法是实现初始场优化反演的有效手段。该研究对进一步改进悬沙输运模型的初始化方案具有一定的参考价值,也对其他数值模型的初始化方案具有一定的借鉴价值。 相似文献
20.
本文使用HYCOM数值模式,根据两种海气通量数据集(COADS、ECMWF)和两种海气通量块体参数化方案(常数块体参数化方案和非常数块体参数化方案)的不同结合,构成4组数值实验,分别模拟了赤道及北太平洋的气候态海表温度.实验结果表明:1)在本文的实验中,非常数块体参数化方案优于常数块体参数化方案;在太平洋40°N-20°S区域内,采用前者得到的年平均海表温度比Pathfinder卫星资料高约0.21 °C,而采用后者得到的年平均海表温度比Pathfinder卫星资料高约0.63 °C.2)HYCOM数值模式很好的模拟了赤道及北太平洋的气候态海表温度变化及西太平洋暖池空间分布的月变化.特别是实验2(采用COADS数据集和非常数块体参数化方案),在太平洋40°N - 20°S区域内,冬春两季平均SST仅比Pathfinder卫星数据集高0.02 °C.3)不同的海气通量数据会对模式结果产生明显的影响.对比采用COADS数据集的实验2结果与采用ECMWF数据集的实验4结果可以发现,在模拟区域的西北部,实验2比实验4的年平均SST高约1 °C;在模拟区域的东南部,实验4比实验2的年平均SST高约1 °C.两者差的最大值出现在58°N、140°E附近及中国渤海,约为4 °C(实验2比实验4的年平均SST高约4 °C). 相似文献