边界适应曲线网格在天文潮、风暴潮数值模拟中的应用

采用曲线坐标系下的准三维近岸流波流耦合数值模型,采用边界适应曲线计算网格和两层嵌套方法,设计一套简单实用的天文潮预报系统,并引入风场以及考虑天文潮和风暴潮之间非线性作用的开边界水位,对2003年渤海发生的一次温带风暴潮进行模拟,计算结果与实测潮位吻合较好。     

台州湾波浪条件数值模拟研究

依据台州湾外海的风浪资料建立数学模型,应用抛物形缓坡方程、Boussinesq方程及小风区成浪对工程区波浪条件进行计算,得到了工程区不同重现期波浪要素,并通过分析得到了工程区波浪年分频分级情况。通过计算与分析可知台州湾外海为台风大浪区;台州湾内由于受众多岛屿掩护,经过地形的折射、绕射和破碎影响,与外海相比波浪衰减达到50%左右;工程区NE向H1%波高约为6.5m,H4%波高超过1.0m的频率每年在20d以下,建港的波浪条件较好。     

The wind-forced response on a buoyant coastal current: Observations of the western Gulf of Maine plume

The freshwater plume in the western Gulf of Maine is being studied as part of an interdisciplinary investigation of the physical transport of a toxic alga. A field program was conducted in the springs of 1993 and 1994 to map the spatial and temporal patterns of salinity, currents and algal toxicity. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman drift dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to as little as 10 km.Using a simple slab model, we find qualitative agreement between the observed variations of plume width and those predicted by Ekman theory for short time scales of integration. Near surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate a dominant Ekman balance within the plume. A significant correlation between alongshore currents and winds suggests that interfacial drag may be important, although inclusion of a Raleigh drag term does not significantly improve the alongshore momentum balance.     

Modelling physical processes of haline stratification in ROFIs with application to the Rhine outflow region

A physical and numerical study is made of the processes governing the stratification and circulation in ROFIs (Regions of Freshwater Influence) where there is an important impact of wind and tides. Observations in the Rhine ROFI showed that the salinity field consists of a mean and a tidally oscillating part. The physical processes are first analysed using the analytical solutions from a one-dimensional two-layer model. A justification is given for the neglection of non-linear advective terms in the equations of momentum and salinity. The dimensionless forms of the solutions can be expressed in terms of a series of dimensionless numbers. It is shown in particular that stratification and cross-shore circulation largely depend on the balance between rotation and turbulent diffusion, which depends in turn on parameters such as the Ekman number, the bottom friction coefficient, the eddy viscosity ratio and the depth of the layer interface. Surface winds either enhance or destroy stratification depending on the wind angle. The response to wind forcing is discussed using classical Ekman theory. To verify the analytical theory numerical tests are performed with a point model including an advanced turbulence closure scheme. Differences arise due to the non-linear interaction between turbulence on the one hand and current shear and stratification on the other hand. It is shown in particular that the amplitude of the tidal forcing and the off-shore horizontal salinity gradient strongly affect the semi-diurnal and semi-monthly variation of stratification. The effect of the wind is found to be in good agreement with the analysis of the two-layer model. Finally, the numerical model is compared with existing observational data in the Rhine ROFI for October 1990.     

Variability of the western Galician upwelling system (NW Spain) during an intensively sampled annual cycle. An EOF analysis approach

The key features of the western Galician shelf hydrography and dynamics are analyzed on a solid statistical and experimental basis. The results allowed us to gather together information dispersed in previous oceanographic works of the region. Empirical orthogonal functions analysis and a canonical correlation analysis were applied to a high-resolution dataset collected from 47 surveys done on a weekly frequency from May 2001 to May 2002. The main results of these analyses are summarized bellow. Salinity, temperature and the meridional component of the residual current are correlated with the relevant local forcings (the meridional coastal wind component and the continental run-off) and with a remote forcing (the meridional temperature gradient at latitude 37°N). About 80% of the salinity and temperature total variability over the shelf, and 37% of the residual meridional current total variability are explained by two EOFs for each variable. Up to 22% of the temperature total variability and 14% of the residual meridional current total variability is devoted to the set up of cross-shore gradients of the thermohaline properties caused by the wind-induced Ekman transport. Up to 11% and 10%, respectively, is related to the variability of the meridional temperature gradient at the Western Iberian Winter Front. About 30% of the temperature total variability can be explained by the development and erosion of the seasonal thermocline and by the seasonal variability of the thermohaline properties of the central waters. This thermocline presented unexpected low salinity values due to the trapping during spring and summer of the high continental inputs from the River Miño recorded in 2001. The low salinity plumes can be traced on the Galician shelf during almost all the annual cycle; they tend to be extended throughout the entire water column under downwelling conditions and concentrate in the surface layer when upwelling favourable winds blow. Our evidences point to the meridional temperature gradient acting as an important controlling factor of the central waters thermohaline properties and in the development and decay of the Iberian Poleward Current.     

A two-dimensional response to a tropical storm on the Gulf of Mexico shelf

Surface current data from drifting buoys and remotely sensed wind data recorded over the continental shelf in the northeastern Gulf of Mexico during the passage of tropical storm Josephine in October 1996 are examined. Drifter data show the existence of a strong surface jet (velocities reaching 1 m s⁻¹) that moves up the west Florida shelf and westward along the Louisiana–Texas shelf, and lasts for nearly 1 week. The coastal jet occurs during an intense synoptic scale wind event where wind speeds reach 15 m s⁻¹. A simple force balance and statistical analysis are performed to assess the role of strong wind forcing. The primary balance shows an Ekman-type current. The role of local acceleration is greatest when winds are directed along bathymetry. A simple two-dimensional strongly forced shelf response model developed from the linear steady-state momentum equations also indicates larger along-shore currents due to both Ekman-type forcing by cross-shore winds and a cross-shore pressure gradient arising from conservation of mass. Model parameters fit empirically are within 15% of theoretical values. The simple model explains 30% and 46% of the variance in the observed along-shore and cross-shore surface currents, respectively.     

Variability of river plumes off Northwest Iberia in response to wind events

The Western Iberian Buoyant Plume (WIBP) is a low-salinity lens formed by river discharge and continental run-off extending along the shelf off Northwest Iberia. The variability of this structure is evaluated with a numerical model forced by real meteorological data and climatologic river discharge during late 2002, when conditions were those of a typical autumn. The direction and intensity of the wind-induced Ekman transport, but also the previous conditions and the duration of the event are found to determine plume behavior. We have identified three characteristic situations: a) confinement of the plume to the coast during downwelling — southerly-winds, b) expansion of the plume during the declining phase of the downwelling event by relaxation of the wind, and c) expansion of the plume by upwelling — northerly-winds. The short time scale of the response of the plume (1–3 h) adds timing between wind events and the phase of the tide as an additional source of variability. In all cases the Iberian Poleward Current (IPC), a saltier and warmer poleward current flowing over the slope, responds as well to wind changes. Furthermore, our simulations illustrate how topography and differences in the river discharge induce local differences in dynamics. Comparisons to available observations show a reasonable model skill. Differences between wind measurements and wind forcing applied to the model appear to be a major source of uncertainty in model results.     

海洋平台设计中的风浪参数选取标准

以我国渤海某区1970年至1993年风暴过程的后报资料为基础,提出了泊松二维冈贝尔逻辑分布,并将其用于风暴过程中伴生的风速与波高的联合概率计算,给出了美国石油协会(API)行业标准有关环境条件三种设计标准的具体做法.基于新的统计模式,同时提出了联合概率法设计标准.比较了多种标准所得设计参数的差异,给出了适合工程所在海区的环境荷载参数.计算结果表明,新的统计模式适用于遭受风暴影响的海洋工程结构设计,特别是边际油田的开发建设.     

Circulation features in the Japan/East Sea related to statistically obtained wind patterns in the warm season

A response of the circulation in the Japan/East Sea (JES) to different kinds of wind forcing is studied, with the emphasis on the warm season, using a primitive equation oceanic model. Wind forcing is based on typical patterns obtained from complex empirical orthogonal functions of 1°-gridded NCEP/NCAR 6 h winds for 1998–2005. These patterns are distinguished by a prevailing wind direction. Northwestern wind and strong cyclonic (C) curl prevail in winter, while a variety of patterns occur in the warm season, differing in the wind direction and curl. Three model runs are performed to examine the circulation in response to a prevailing C wind stress curl or an alternating C and anticyclonic (AC) curl or a strong C curl in the warm season. The simulated features are consistent with the observational evidence, in particular with thermal fronts and frequent eddy locations derived from multi-year infrared satellite imagery. The simulated C circulation intensifies and the subarctic region extends southward with the strengthening of a summer C wind stress curl over the JES. Variability of Subarctic Front (SF) in the western JES (between 130°E and 133°E) is strongly affected by summer wind stress curl. Forcing by an AC curl tends to shift SF northward, while SF shifts to the south under the forcing by a C curl, reaching the southern Ulleung Basin in the case of the strong C curl. In the northwestern JES (off Peter the Great Bay, Russia, and North Korea), the SF northwestern branch (NWSF) is simulated. It is a known feature in autumn and early winter and can also occur in the warm season. The simulation results suggest an AC wind stress curl as the forcing of the formation of the NWSF in the warm season. The Siberia Seamount and sharply bending coastline near Peter the Great Bay facilitate partial separation of the Primorye (Liman) Current from the coast. The wind stress curl can be an additional forcing of the Tsushima Warm Current (TWC) branching off the Korea Strait to the East Korea Warm Current (EKWC) and the offshore branch (OB). In the warm season, the simulated TWC bifurcation occurs farther north, the EKWC is strong, and the OB is weak under the forcing of the AC wind stress curl. The EKWC is weak and the OB is strong under the forcing of the strong C wind stress curl.     

江苏海域波浪分布特征研究

以1991—2011年风场数据为驱动,运用第3代波浪模型SWAN对江苏海域的波浪分布特征进行了模拟研究。首先运用基于浅水方程的Delft3D-flow模块对江苏海域的潮波进行模拟,给波浪模型提供一个准确的边界条件;其次运用SWAN模型对江苏海域的波浪进行模拟并与观测值进行比较;最后对江苏海域的波浪分布特征进行研究。结果表明:江苏海域年平均波高由南往北依次递减,年平均最大值为1.5 m左右,强浪向为NE向,四季中冬季的平均波浪最大,春秋次之,夏季最小。     

Ekman transport along the Galician Coast (NW, Spain) calculated from QuikSCAT winds

Ekman transport is studied close to the Galician coast by means of wind data provided by the QuikSCAT satellite from November 1999 to October 2005. Three different coastal zones are identified, western coast from Miño River to Cape Finisterre, middle coast from Cape Finisterre to Cape Ortegal and northern coast, from Cape Ortegal to Cape Peñas. In addition to existence of long-term variations, the periodicity of the transport signal is characterized by an annual component (365 days), a seasonal fluctuation (50–80 days) and a time scale related to passing storms (15–20 days). Although the periodicity of the signal is similar at the three zones due to external meteorological forcing, the Ekman transport is modulated by the presence of the coast, in such a way that seasonal patterns vary in intensity and direction along the coast. Thus, the spring–summer pattern is characterized by high transport at the western coast, pointing seaward perpendicular to the shore-line. The same orientation is observed at the middle coast although with a lower magnitude. Finally, Ekman transport at the northern coast points landward and oblique to the shore-line. The different transport orientations are shown to be responsible for the upwelling probability variation along the coast. On the other hand, the autumn–winter pattern does not show a clear trend with important inter-annual differences showing the high variability of Ekman transport for this period.     

Lagrangian data in a high-resolution numerical simulation of the North Atlantic: I. Comparison with in situ drifter data

A model/data comparison was performed between simulated drifters from a high-resolution numerical simulation of the North Atlantic and a data set from in situ surface drifters. The comparison makes use of pseudo-Eulerian statistics such as mean velocity and eddy kinetic energy, and Lagrangian statistics such as integral time scales. The space and time distribution of the two data sets differ in the sense that the in situ drifters were released inhomogeneously in space and time while the simulated drifters were homogeneously seeded at the same time over a regular 1° grid. Despite this difference, the total data distributions computed over the complete data sets show some similarities that are mostly related to the large-scale pattern of Ekman divergence/convergence.Comparisons of eddy kinetic energy and root mean square velocity indicate that the numerical model underestimates the eddy kinetic energy in the Gulf Stream extension and in the ocean interior. In addition, the model Lagrangian time scales are longer in the interior than the in situ time scales by approximately a factor of 2. It is suggested that this is primarily due to the lack of high-frequency winds in the model forcing, which causes an underestimation of the directly forced eddy variability. Regarding the mean flow, the comparison has been performed both qualitatively and quantitatively using James' statistical test. The results indicate that over most of the domain, the differences between model and in situ estimates are not significant. However, some areas of significant differences exist, close to high-energy regions, notably around the Gulf Stream path, which in the model lies slightly north of the observed path, although its strength and structure are well represented overall. Mean currents close to the buffer zones, primarily the Azores Current, also exhibit significant differences between model results and in situ estimates. Possibilities for model improvement are discussed in terms of forcings, buffer zone implementations, turbulence and mixed layer parameterizations, in light of our model/data comparison.     

Model test investigation of a spar floating wind turbine

Several floating wind turbine designs whose hull designs reflect those used in offshore petroleum industry have emerged as leading candidates for the future development of offshore wind farms. This article presents the research findings from a model basin test program that investigated the dynamic response of a 1:50 scale model OC3 spar floating wind turbine concept designed for a water depth of 200 m. In this study the rotor was allowed to rotate freely with the wind speed and this approach eliminated some of the undesirable effects of controlling wind turbine rotational speed that were observed in earlier studies. The quality of the wind field developed by an array of fans was investigated as to its uniformity and turbulence intensity. Additional calibration tests were performed to characterize various components that included establishing the baseline wind turbine tower frequencies, stiffness of the delta type mooring system and free decay response behaviour. The assembled system was then studied under a sequence of wind and irregular wave scenarios to reveal the nature of the coupled response behaviour. The wind loads were found to have an obvious influence on the surge, heave and pitch behaviour of the spar wind turbine system. It was observed from the experimental measurements that bending moment at the top of the support tower is dominated by the 1P oscillation component and somewhat influenced by the incoming wave. Further it was determined that the axial rotor thrust and tower-top shear force have similar dynamic characteristics both dominated by tower’s first mode of vibration under wind-only condition while dominated by the incident wave field when experiencing wind-wave loading. The tensions measured in the mooring lines resulting from either wave or wind-wave excitations were influenced by the surge/pitch and heave couplings and the wind loads were found to have a clear influence on the dynamic responses of the mooring system.     

Two parametric tropical cyclone models for storm surge modeling

In this paper,the two parametric tropical cyclone models for storm surge modeling are further developed.The analytical expressions of tangential and radial velocity distribution are derived from the governing momentum equations,based on the general symmetric pressure distribution proposed by Holland and Fujita.On the basis of the data of several tropical cyclones that occurred in East China Ocean,the shape parameter in pressure model is estimated.Finally,the Fred cyclone(typhoon 199417)is calculated,and comparisons of measured and calculated air pressures and wind speed are presented.     

Variability in hydrographic conditions to the east and northwest of South Georgia, 1996–2001

Six years of high-resolution hydrographic data from the eastern and northwestern sides of South Georgia (southwest Atlantic) are used to study the changing circulation and water mass properties of the region. One year of data from these locations was used previously to describe the oceanographic conditions at those times; using the much greater volume of data now available, we identify which features appear temporally robust and which are transient, and begin addressing topics relating to the forcing of the inter-annual variability and the potential consequences for the local ecosystem. Waters on the shelf and those over the adjacent deep ocean invariably have different hydrographic properties, though the transition between them can be abrupt or gradual. The onshelf/offshelf differences vary greatly from year to year, due to the combined influences of local and remote processes. There are several instances of strong physical coupling between the eastern and northwestern sides of South Georgia; this offers potential for distinguishing physically-induced ecosystems changes separately from biologically-induced ecosystems changes. On the northeast side of the shelf, close to Cumberland Bay, there is evidence of an often intense, but variable, cyclonic circulation that is the result of interaction with the local bathymetry. This may act as a retention mechanism, and enhance local productivity. Two examples of extreme cold anomalies are present in the series of measurements. One of these (in 2000/2001) affected a limited area at the eastern side of the region surveyed and was due to an intrusion of the Southern Antarctic Circumpolar Current Front. The other (in early 1998) was due to the combined effects of the passage of a large-scale ocean anomaly that had its origins upstream in the Pacific Ocean, and strong air/sea interaction. Both of these were associated with the strong 1997/1998 El Niño event. Whilst previous studies have observed a link between El Niño forcing and ocean response around South Georgia with a temporal lag of around 3 years, we observe a much more rapid response to the extremely strong 1997/1998 El Niño event. This indicates that the ocean and ecosystem around South Georgia are more immediately susceptible to extreme instances of remote climatic forcing than had been supposed.     

半潜式钻井服务支持平台随机风载荷特性研究

风载荷是海洋平台设计载荷之一,直接关系到平台稳性,确定平台在不同工况下的风荷载对于平台安全设计具有重要的工程意义.该文基于Fluent软件结合自编UDF程序,考虑Davenport、NPD以及API三种典型风谱的影响,利用简谐波叠加法将风谱由频域转换为时域内的随机脉动风速,引入雷诺平均法求解NS方程结合分离涡(DES)湍流模型对半潜式钻井服务支持平台在自存海况下的风力和风倾力矩开展了数值研究,并将数值模拟得到的最大风力及风倾力矩与ABS、DNV以及CCS的计算结果进行了对比验证.计算结果表明:平台受到的风力和风倾力矩与风谱自身特性、平台倾斜角及风向角等因素密不可分;同一工况下采用Davenport与NPD风谱计算时平台受到的平均风力(矩)较为接近;NPD风谱作用时平台受到的随机最大风力(矩)最大;采用API与NPD风谱计算时,各工况下最小风力(矩)随风向角的变化趋势、计算结果均基本一致.     

Evolving the subspace of the three-dimensional multiscale ocean variability: Massachusetts Bay

A data and dynamics driven approach to estimate, decompose, organize and analyze the evolving three-dimensional variability of ocean fields is outlined. Variability refers here to the statistics of the differences between ocean states and a reference state. In general, these statistics evolve in time and space. For a first endeavor, the variability subspace defined by the dominant eigendecomposition of a normalized form of the variability covariance is evolved. A multiscale methodology for its initialization and forecast is outlined. It combines data and primitive equation dynamics within a Monte-Carlo approach.The methodology is applied to part of a multidisciplinary experiment that occurred in Massachusetts Bay in late summer and early fall of 1998. For a 4-day time period, the three-dimensional and multivariate properties of the variability standard deviations and dominant eigenvectors are studied. Two variability patterns are discussed in detail. One relates to a displacement of the Gulf of Maine coastal current offshore from Cape Ann, with the creation of adjacent mesoscale recirculation cells. The other relates to a Bay-wide coastal upwelling mode from Barnstable Harbor to Gloucester in response to strong southerly winds. Snapshots and tendencies of physical fields and trajectories of simulated Lagrangian drifters are employed to diagnose and illustrate the use of the dominant variability covariance. The variability subspace is shown to guide the dynamical analysis of the physical fields. For the stratified conditions, it is found that strong wind events can alter the structures of the buoyancy flow and that circulation features are more variable than previously described, on multiple scales. In several locations, the factors estimated to be important include some or all of the atmospheric and surface pressure forcings, and associated Ekman transports and downwelling/upwelling processes, the Coriolis force, the pressure force, inertia and mixing.     

系泊船舶动力学特性的计算机仿真研究

以三阶操纵运动方程为基础，引入定常的风力、潮流作用力和二阶波浪力，建立了系泊系统三自由度的运动微分方程。在此数学模型的基础上．建立了系泊系统的多自由度的计算机仿真模型。在潮流作用、潮流和风作用以及风浪流联合作用等三种情形下，对系泊船舶的动力学响应分别进行了仿真研究。研究表明，系泊系统的动力学行为具有强烈的非线性特征。对于单点系泊船舶而言，在定常的风浪流作用下．顶风顶浪顶流状态并不一定是最为危险的工况。     

Navigating Institutional Constraints in the Process of Computational Environmental Modeling

Abstract

The use of computational models has proliferated in the last few decades due to the practical benefits that they offer in terms of understanding complex systems. However, the process of developing a computational model is not only shaped by scientific and practical concerns. Modelers must also navigate the institutional structures in which the models are constructed and implemented. In the process, they must often draw upon existing or create new social relationships that can, themselves, alter the institutional structures they are attempting to navigate. In this article, I use an ethnographic approach to examine three examples in which the process of building and implementing a computational model was constrained by institutional factors and the strategies the modelers used to navigate them. The research was conducted with computational modelers at the Chesapeake Bay Program who are involved in constructing the Chesapeake Bay modeling system for nutrient management in the watershed. Modelers in this context had to: build relationships with the broader scientific community to reinforce the “believability” of the model, draw upon institutional relationships in order to work around limitations on data-sharing between organizations, and manage differing incentive structures to motivate and coordinate the research needed to complete the model     

Numerical simulations of seasonal and interannual variability of the Black Sea thermohaline circulation

The Black Sea general circulation is simulated by a primitive equation model with active free surface. The forcing is seasonally variable and is based on available climatic data. The model reproduces the main features of the Black Sea circulation, including the river discharge effects on the mean sea level and the Bosphorus outflow. Model results show that the simulated sea surface elevation increases in spring over the whole sea, reaching a maximum in the Danube delta area. In the same region, a minimum is observed in winter. The amplitude of the seasonal oscillations (about 8–12 cm over the whole basin) is of the same order of magnitude as the maximum horizontal variations (about 15–18 cm between the coastal areas and the basin interior). This strong signal formed mostly by river discharges, along with the seasonal variability in the other forcing functions and the local dynamics creates a well-pronounced interannual variability. The performance of the model in simulating the seasonal and interannual variability is critically analyzed, with a special attention on the cold intermediate water formation and the circulation in the upper 150 m. The simulations demonstrate that the source of intermediate waters is on the shelf, and that the water mass in the core of cold intermediate layer changes with time as a response to the periodic forcing at sea surface. This type of variability is characterized by pronounced interannual changes, proving that important differences could exist between water mass structure in different years, even when using identical atmospheric forcings each year.