ECOM模式在丁字湾的应用

应用河口、陆架和海洋模式 (ECOM模型 ) ,引入干湿网格法模拟潮滩涨落的改进 ,并建立丁字湾及近岸海域的三维变动边界潮流模型。该模型考虑了湾口拦门沙、湾内水道和人工围海等地形特点。计算结果与实测值比较符合良好 ,较好地刻画出丁字湾 M2 分潮潮流场的时空分布特点。     

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.     

3D gravity modelling for Anyongbok Seamount in the East Sea

A seamount chain with an approximately WNW trend is observed in the northeastern Ulleung Basin. It has been argued that these seamounts, including two islands called Ulleung and Dok islands, were formed by a hotspot process or by ridge related volcanism. Many geological and geophysical studies have been done for all the seamounts and islands in the chain except Anyongbok Seamount, which is close to the proposed spreading ridge. We first report morphological characteristics, sediment distribution patterns, and the crustal thickness of Anyongbok Seamount using multibeam bathymetry data, seismic reflection profiles, and 3D gravity modeling. The morphology of Anyongbok Seamount shows a cone shaped feature and is characterized by the development of many flank cones and flank rift zones. The estimated surface volume is about 60 km³, and implies that the seamount is smaller than the other seamounts in the chain. No sediments have been observed on the seamount except the lower slope, which is covered by more than 1,000 m of strata. The crustal structure obtained from a 3D gravity modeling (GFR = 3.11, SD 3.82 = mGal) suggests that the seamount was formed around the boundary of the Ulleung Plateau and the Ulleung Basin, and the estimated crustal thickness is about 20 km, which is a little thicker than other nearby seamounts distributed along the northeastern boundary of the Ulleung Basin. This significant crustal thickness also implies that Anyongbok Seamount might not be related to ridge volcanism.     

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.