A hybrid model for numerical wave forecasting and its implementation——Ⅰ. The wind wave model

The authors make an endeavor to explain why a new hybrid wave model is here proposed when several such models have already been in operation and the so- called third generation wave modej is proving attractive. This part of the paper is devoted to the wind wave model. Both deep and shallow water models have been developed, the former being actually a special case of the latter when water depth is great. The deep water model is exceptionally simple in form. Significant wave height is the only prognostic variable. In comparison with the usual methods to compute the energy input and dissipations empirically or by "tuning", the proposed model has the merit that the effects of all source terms are combined into one term which is computed through empirical growth relations for significant waves, these relations being, relatively speaking, easier and more reliable to obtain than those for the source terms in the spectral energy balance equation. The discrete part of the model and the implementation of the mode     

LAGFD-WAM numerical wave model-Ⅱ. Characteristics inlaid scheme and its application

In this paper the parameterizational approach of nonlinear source function and the implicit scheme of the model are discussed in detail. The matching problem is solved between time and space steps using the characteristics inlaid scheme with very strong physical meaning. The computational comparison in typical winds shows some improvements to the WAM model. That the hindcast results of the model for typhoon cases are in good agreement with real data illustrates its applicability to wave forecast and engineering study.     

A hybrid model for numerical wave forecasting and its implementation——Ⅱ. The discrete part and implementation of the model

In the first part of the present paper we have explained why we manage to formulate another wave prediction model when so many of them, including the so-called third generation model, have already been in use. The wind-wave part of the proposed model has also been given. Now we proceed to discuss the swell part,the implementation of the model as a prediction method,mumerical experiments done with ideal wind fields and hindcasts made in the Bohai Sea,in the neighboring seas adjacent to China and in the Northwest Pacific.     

Quadtree grid numerical model of nearshore wave–current interaction

This paper describes an adaptive quadtree-based 2DH wave–current interaction model for evaluating nearly horizontal wave-induced currents in the surf-zone. The model accounts for wave breaking, shoaling, refraction, diffraction, wave–current interaction, set-up and set-down, mixing processes (turbulent diffusion), bottom frictional effects, and movement of the land–water interface at the shoreline. The wave period- and depth-averaged governing equations, which conserve mass, momentum, energy and wave action, are discretised explicitly by means of an Adams–Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Grid adaptation is achieved through seeding points distributed according to flow criteria (e.g. local current gradients). The model is verified for nearshore circulation at a sinusoidal beach and nearshore currents at a multi-cusped beach. Reasonable agreement is obtained with experimental data from da Silva Lima [da Silva Lima, S.S.L., 1981. Wave-induced Nearshore Currents. PhD Thesis, Department of Civil Engineering, University of Liverpool] and Borthwick et al. [Borthwick, A.G.L., Foote, Y.L.M., Ridehalgh, A., 1997. Nearshore measurements at a cusped beach in the UK Coastal Research Facility, Coastal Dynamics '97, Plymouth, 953–962]. The modelling approach presented herein should be useful in simulating nearshore processes in complicated natural coastal domains. Of particular value is the local grid enrichment capability, which permits refined modelling of important localised flow behaviour such as rip currents and surf-zone circulation systems.     

Theoretical wind wave frequency spectra in deep water——Ⅰ. Form of spectrum

In this part ot the paper theoretical wind-wave spectra nave been derived oy (I) expressing the spectrum in series composed of exponential terms; (2) assuming that the spectrum satisfies a high order linear ordinary differential equation; (3) introducing proper parameters in the spectrum; and (4) making use of some known charateristics of wind-wave spectrum, for instance, the law governing the equilibrium range. The spectrum obtained contains the zero order moment of the spectrum m0, the peak frequency ω0 and the ratio R =ω/ω0 (ω being the mean zero-crossing frequency) as parameters. The shape of the nondimensional spectrum S(ω) = ω0S(ω)/m0(ω=ω/ω0) changes with R and theoretically reduces to a Dirac delta function δ(ω-1) when R = 1. A spectrum of simplified form is given for practical uses, in which R is replaced by a peakness factor P=S(1).     

LAGFD-WAM海浪数值模式——Ⅰ:基本物理模型

LAGFD-WAM型海浪数值模式是一种第三代海浪数值模式,本文详细地介绍了该模式的物理模型,给出了能谱平衡方程、复杂特征线方程、各种源函数,以及物理空间和相空间的边界条件。     

Dissipation source function and an improvement to LAGFD-WAM model

In this paper, a new theoretical expression of dissipation term is presented on the basis of statistical model of breaking wave, which is an improvement to LAGFD-WAM wave model. The computational results in three typical wind fields show a good improvement to LAGFD-WAM model and a better accuracy in comparison with the observed data in the South China Sea.     

CHGS method for numerical forecasting of typhoon waves——Ⅰ. Spectrum of waves in growing phase

Owing to the fact that the wind speed and direction of typhoon vary rapidly with time and space in typhoon fetch; the nearer to the typhoon eye the greater the wind velocity, and the shorter the wind fetch the smaller the wind time,as a result,the more difficult for the wind wave to fully grow. Hence.in typhoon wave numerical calculation it is impossible to use the model for a fully grown wave spectrum. Lately, the author et at. presented a CHGS method for numerical forecasting of typhoon waves, where a model for the growing wave spectrum was set up (see Eq. (2) in the text). The model involves a parameter indicating the growing degree of wind wave, i. e. ,the mean wave age β. When βvalue is small, the wave energy is chiefly concentrated near the peak frequency, so that the spectral peak gets high and steep; with the increase of β the spectral shape gradually gets lower and gentler; when β=Ⅰ, the wave fully grows, the growing spectrum becomes a fully grown P-M spectrum. The model also shows a spect     

A study on the numerical prediction method for the vertical thermal structure in the Bohai Sea and the Huanghai Sea-Ⅰ. One-dimensional numerical prediction model

In this paper,on the basis of the heat conduction equation without consideration of the advection and turbulence effects,one-dimensional model for describing surface sea temperature (T_s),bottom sea temperature (T_h) and the thickness of the upper homogeneous layer (h) is developed in terms of the dimensionless temperature θ_T and depth η and self-simulation function θ_T-f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T_s,T_h,h,the thickness and gradient of thermocline are satisfactory to some extent.     

A numerical model of mesoscale frontal instabilities and plankton dynamics — I. Model formulation and initial experiments

Previous observational and modelling studies of open ocean frontal regions have found large vertical velocities associated with instabilities on the frontal jet. A combined physical/ecosystem numerical model is used to investigate the impact of jet instability and the associated vertical motions on the local ecosystem. The evolution of the instability of a mesoscale frontal jet gives rise to vertical transport of nutrients into the euphotic zone and subduction of biota out of the euphotic zone. The upwelling of nutrients stimulates increases in primary production, with resulting increases in phytoplankton stocks. The reaction of the ecosystem is found to be dependent on the physical characteristics of the front, but the increase in primary production can be locally of the order of 100%, and of the order of 10% when averaged over the frontal region. The action of upwelling and subduction introduces spatial heterogeneity in primary production and plankton biomass. The heterogeneity is at a variety of length scales, from the order of a few kilometres for thin filaments and up to 50 km for coherent features. With increases in new production occurring over several degrees of latitude, frontal dynamics may make a significant contribution to the strength of the biological pump.     

A hybrid Lagrangian-Eulerian numerical model for sea-ice dynamics

A hybrid Lagrangian-Eulerian(HLE) method is developed for sea ice dynamics,which combines the high computational efficiency of finite difference method(FDM) with the high numerical accuracy of smoothed particle hydrodynamics(SPH).In this HLE model,the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations.These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function.The FDM is used to determine the ice velocities at Eulerian grid nodes,and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also.The thicknesses and concentrations of ice particles are determined based on their new locations.With the HLE numerical model,the ice ridging process in a rectangular basin is simulated,and the simulated results are validated with the analytical solution.This method is also applied to the simulation of sea ice dynamics in a vortex wind field.At last,this HLE model is applied to the Bohai Sea,and the simulated concentration,thickness and velocity match the satellite images and the field observed data well.     

Three-dimensional numerical simulation of obfique wave action on vertical walls

A three-dimensional numerical model based on the potential theory was developed to study the oblique wave action on vertical walls. A source term inside the domain was used to generate incident waves and outgoing waves were dissipated by sponge layers and transmitted by radiation boundaries. The finite difference method was used to solve the governing equations and boundary conditions in the regular transformed domain in σ-coordinate. Satisfactory agreements between the numerical predictions and experimental results of wave force were obtained. It is concluded that the maximum wave force acting on the vertical walls is induced by the obliquely incident waves rather than the normally incident waves.     

A model for prediction SST in limited region-Ⅱ. The model''s physical equation

-The model's physical equation is used to parameterize some subgrid-scale processes and physical processes in the present numerical model. The transmission and attenuation of the solar energy in the ocean are considered. A simple diagnostic equation for the cloud fractions k proposed on the basis of the humidities at the surface and the mid-troposphere. The parameterized formulae of both entrainment and Ekman pumping are improved.In the numerical integration, the treatment on damping the inertial oscillations is emphasized. The initialization and the objective analysis of the data which are necessary for the operational prediction will be presented in another paper.Results of SST prediction and some numerical experiments are given here. The model is computationally stable and successful in modelling the behaviors of the drift current and the mixed layer physics, and the AMD (absolute mean deviations) ≤1. 2℃ , RC (correlation coefficients ) ≥85% for 3-day forecasting.     

A coupled ice-ocean model for the Bohai Sea Ⅰ. Study on model and parameter

According to the earlier international studies on the coupled ice-ocean model and the hydrology, meteorology, and ice features in the Bohai Sea, a coupled ice-ocean model is developed based on the National Marine Environment Forecast Center's (NMEFC) numerical forecasting ice model of the Bohai Sea and the Princeton ocean model (POM).In the coupled model, the transfer of momentum and heat between ocean and ice is two-way, and the change of ice thickness and concentration depends on heat budget not only at the surface and bottom of ice, but also at the surface of open water between ices. The dynamic and thermodynamic coupling process is expatiated emphatically. Some thermodynamic parameters are discussed as well.     

The application of wave refraction-diffraction model with friction

In this paper a refraction-diffraction model with friction is used to compute wave characteristics in a region near a certain port. Comparing with the results from refraction model, and with the data observed during a typhoon in 1985, it is found that the characteristics from the refraction-diffraction model with friction are reasonable, and that the results are in rather good agreement with observations. Thus it can be concluded that the model is effective for computing coastal wave characteristics over complicated bottom topography.     

A coupled ice-ocean model for the Bohai Sea Ⅰ. Study on model and parameter

According to the earlier international studies on the coupled ice-ocean model and the hydrology, meteorology, and ice features in the Bohai Sea, a coupled ice-ocean model is developed based on the National Marine Environment Forecast Center‘ s (NMEFC) numerical forecasting ice model of the Bohai Sea and the Princeton ocean model (POM). In the coupled model, the transfer of momentum and heat between ocean and ice is two-way, and the change of ice thickness and concentration depends on heat budget not only at the surface and bottom of ice, but also at the surface of open water between ices. The dynamic and thermodynamic coupling process is expatiated emphatically. Some thermodynamic parameters are discussed as well.     

Study on short-range numerical forecasting of ocean current in the East China Sea——Ⅰ Basic problems of ocean current forecasting and structure of the models

Ocean current forecasting is still in explorative stage of study. In the study, we face some problems that have not been met before. The solving of these problems has become fundamental premise for realizing the ocean current forecasting. In the present paper are discussed in depth the physical essence for such basic problems as the predictability of ocean current, the predictable currents, the dynamical basis for studying respectively the tidal current and circulation, the necessity of boundary model, the models on regions with different scales and their link. The foundations and plans to solve the problems are demonstrated. Finally a set of operational numerical forecasting system for ocean current is proposed.     

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅰ. Meander simulation and linear growth rate

Ａｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｏｃｅａｎｇｅｎｅｒａｌｃｉｒｃｕｌａｔｉｏｎｍｏｄｅｌｆｏｒｍｅｓｏｓｃａｌｅｅｄｄｉｅｓ──ⅠＭｅａｎｄｅｒｓｉｍｕｌａｔｉｏｎａｎｄｌｉｎｅａｒｇｒｏｗｔｈｒａｔｅ￥ＷａｎｇＪｉａａｎｄＭｏｔｏＩｋｅｄａ（Ｒ...