Two-Dimensional Mathematical Model of Tidal Current and Sediment for Oujiang Estnary and Wenzhou Bay

A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the WenzhouBay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport arepresented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with theirinitial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of thecritical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested.The model verification computation shows that, the tidal levels computed from the model are in good agreement with the fielddata at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measuredalong the totally 52 synchronously observed verticals distributed over 8 cross sections. The computed tidal water throughputsthrough the Huangda‘ao cross section are close to the measured data. And the computed values of bed deformation fromYangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992.The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of thereclamation project on the Wenzhou shoal are predicted by means of this model.     

数字地球与海洋科学的发展

本文从数字地球的概念出发，根据海洋科学的发展要求，详细论述了数字地球与海洋科学的发展关键。并提出了海洋科学的新的发展战略和配套技术的开发研究理论，从而为我国海洋事业的发展和数字地球的建立有一个协调和统一的认识观点。     

Simulation of Formation and Spreading of Salinity Minimum Associated with NPIW Using a High-Resolution Model

A series of numerical experiments were conducted with a high-resolution (eddy-permitting) North Pacific model to simulate the formation and spreading of the salinity minimum associated with the North Pacific Intermediate Water (NPIW). It was found that two factors are required to simulate a realistic configuration of the salinity minimum: a realistic wind stress field and small-scale disturbances. The NCEP reanalyzed wind stress data lead to better results than the Hellerman and Rosenstein wind stress data, due to the closer location of the simulated Oyashio and Kuroshio at the western boundary. Small-scale disturbances formed by relaxing computational diffusivity included in the advection scheme promote the large-scale isopycnal mixing between the Oyashio and Kuroshio waters, simulating a realistic configuration of the salinity minimum. A detailed analysis of the Oyashio water transport was carried out on the final three-year data of the experiment with reduced computational diffusivity. Simulated transport of the Kuroshio Extension in the intermediate layer is generally smaller than the observed value, while those of the Oyashio and the flow at the subarctic front are comparable to the observed levels. In the Oyashio-Kuroshio interfrontal zone the zonally integrated southward transport of the Oyashio water (140–155°E) is borne by the eddy activity, though the time-mean flow reveals the existence of a coastal Oyashio intrusion. In the eastern part (155°E–180°) the zonally integrated transport of the Oyashio water indicates a southward peak at the southern edge of the Kuroshio Extension, which corresponds to the branching of the recirculating flow from the Kuroshio Extension. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Three‐dimensional simulation of tsunami run‐up around conical island

At the circular Babi Island in the Flores tsunami (1992) and pear shaped island in the Okushiri event (1993), unexpectedly large tsunami run‐up heights in the lee of conic islands were observed. The flume and basin physical model studies were conducted in the Coastal Hydraulic Laboratory, Engineering Research and Development Center, U.S. Army Corps of Engineers to provide a better understanding of the physical phenomena and verify numerical models used in predicting tsunami wave run‐up on beaches, islands, and vertical walls. Reasonably accurate comparison of run‐up height of solitary waves on a circular island has been obtained between laboratory experimental results and two‐dimensional computation model results. In this study we apply three‐dimensional RANS model to simulate wave run‐up on conical island. In the run‐up computation we obtain that 3D calculations are in very good comparison with laboratory and 2D numerical results. A close examination of the three‐dimensional velocity distribution around conical island to compare with depth‐integrated model is performed. It is shown that the velocity distribution along the vertical coordinate is not uniform: and velocity field is weaker in the bottom layer and higher on the sea surface. The maximum difference (about 40%) appears at the time when solitary wave reached the circular island.