Development of a High-Resolution Coastal Circulation Model for the Ocean Observatory in Lunenburg Bay

An advanced ocean observatory has been established in Lunenburg Bay of Nova Scotia, Canada as part of an interdisciplinary research project of marine environmental prediction. The development of a high-resolution coastal circulation model is one of important components of the observatory. The model horizontal resolution is 60m and the vertical resolution is about 1 m. The coastal circulation model is used to simulate the semi-diurnal tidal circulation and associated nonlinear dynamics with the M2 forcing specified at the model open boundaries. The model is also used to simulate the storm-induced circulation in the bay during Hurricane Juan in September 2003, with the model forcing to be the combination of tides and remotely generated waves specified at the model open boundaries and wind stress applied at the sea surface. The model results demonstrate strong interactions between the local wind stress, tidal forcing, and remotely generated waves during this period. Comparison of model results with the surface elevation and current observations demonstrates that the coastal circulation model has reasonable skills in simulating the tidal and storm-induced circulation in the bay.     

The effect of Coriolis-Stokes forcing on upper ocean circulation in a two-way coupled wave-current model

We investigated the Stokes drift-driven ocean currents and Stokes drift-induced wind energy input into the upper ocean using a two-way coupled wave-current modeling system that consists of the Princeton Ocean Model generalized coordinate system (POMgcs), Simulating WAves Nearshore (SWAN) wave model, and the Model Coupling Toolkit (MCT). The Coriolis-Stokes forcing (CSF) computed using the wave parameters from SWAN was incorporated with the momentum equation of POMgcs as the core coupling process. Experimental results in an idealized setting show that under the steady state, the scale of the speed of CSF-driven current was 0.001 m/s and the maximum reached 0.02 m/s. The Stokes drift-induced energy rate input into the model ocean was estimated to be 28.5 GW, taking 14% of the direct wind energy rate input. Considering the Stokes drift effects, the total mechanical energy rate input was increased by approximately 14%, which highlights the importance of CSF in modulating the upper ocean circulation. The actual run conducted in Taiwan Adjacent Sea (TAS) shows that: 1) CSF-based wave-current coupling has an impact on ocean surface currents, which is related to the activities of monsoon winds; 2) wave-current coupling plays a significant role in a place where strong eddies present and tends to intensify the eddy’s vorticity; 3) wave-current coupling affects the volume transport of the Taiwan Strait (TS) throughflow in a nontrivial degree, 3.75% on average.     

Modeling seasonal variations of ocean and sea ice circulation in the Beaufort and Chukchi Seas： A model-data fusion study

A 3.8-kin Coupled Ice-Ocean Model （C1OM） was implemented to successfully reproduce many observed phenomena in the Beaufort and Chukchi seas, including the Bering-inflow-originated coastal current that splits into three branches： Alaska Coastal Water （ACW） , Central Channel, and Herald Valley branches. Other modeled phenomena include the Beaufort Slope Current （BSC） , the Beaufort Gyre, the East Siberian Current （ ESC）, mesoscale eddies, seasonal landfast ice, sea ice ridging, shear, and deformation. Many of these downscaling processes can only be captured by using a high-resolution CIOM, nested in a global climate model. The seasonal cycles for sea ice concentration, thickness, velocity, and other variables are well reproduced with Solid validation by satellite measurements. The seasonal cycles for upper ocean dynamics and thermodynamics are also well reproduced, which include the formation of the cold saline layer due to the injection of salt during sea ice formation, the BSC, and the subsurface upwelling in winter that brings up warm, even more saline Atlantic Water along the shelfbreak and shelf along the Beaufort coast.     

A FVCOM-based unstructured grid wave, current, sediment transport model, I. Model description and validation

An effort was made to couple FVCOM (a three-dimensional (3D),unstructured grid,Finite Volume Coastal Ocean Model) and FVCOM-SWAVE (an unstructured grid,finite-volume surface wave model) for the study of nearshore ocean processes such as tides,circulation,storm surge,waves,sediment transport,and morphological evolution.The coupling between FVCOM and FVCOM-SWAVE was achieved through incorporating 3D radiation stress,wave-current-sediment-related bottom boundary layer,sea surface stress parameterizations,and morphology process.FVCOM also includes a 3D sediment transport module.With accurate fitting of irregular coastlines,the model provides a unique tool to study sediment dynamics in coastal ocean,estuaries,and wetlands where local geometries are characterized by inlets,islands,and intertidal marsh zones.The model was validated by two standard benchmark tests: 1) spectral waves approaching a mild sloping beach and 2) morphological changes of seabed in an idealized tidal inlet.In Test 1,model results were compared with both analytical solutions and laboratory experiments.A further comparison was also made with the structured grid Regional Ocean Model System (ROMS),which provides an insight into the performance of the two models with the same open boundary forcing.     

A THREE-LAYER OCEAN CIRCULATION MODEL APPLICATION TO NUMERICAL STUDIES ON THE NORTH PACIFIC OCEAN CIRCULATION

The formulation and justification of a three-layer baroclinic ocean model developed to simulate thegeneral circulation of the ocean are described in this paper.Test of the model in simulating the annualmean circulation patterns in the North Pacific under the prescribed atmospheric forcing,which consists ofthe climatological surface wind stress and sea surface heat flux,and comparison of the results withobservations showed that the model basically simulated the large scale features of the annual meancirculation patterns in the North Pacific Ocean such as those of the intensified western boundary currentsand the North Equatorial Currents and Undercurrents.But due to the coarse resolution of the model,some details of these currents were poorly reproduced.The seasonal variations of the North Pacific Oceancirculation driven by the seasonal mean sea surface wind stress was calculated,the different aspects of theseresults were analyzed and the main current(the intensified western boundary currents)transports we     

Impacts of wave and current on drag coefficient and wind stress over the tropical and northern Pacific

By taking into consideration the effects of ocean surface wave-induced Stokes drift velocity U,w and current velocityU,c on the drag coefficient,the spatial distributions of drag coefficient and wind stress in 2004 are computed over the tropical andnorthern Pacific using an empirical drag coefficient parameterization formula based on wave steepness and wind speed.The globalocean current field is generated from the Hybrid Coordinate Ocean Model (HYCOM) and the wave data are generated from Wave-watch Ⅲ (WW3).The spatial variability of the drag coefficient and wind stress is analyzed.Preliminary results indicate that theocean surface Stokes drift velocity and current velocity exert an important influence on the wind stress.The results also show thatconsideration of the effects of the ocean surface Stokes drift velocity and current velocity on the wind stress can significantly im-prove the modeling of ocean circulation and air-sea interaction processes.     

全球海潮模型最新进展及在中国沿海精度评估

介绍TPXO、FES、Chinatide、MIKE Global Tide、Utide等典型海潮模型，总结归纳其同化潮汐数据来源和最新的海洋地形数据，利用我国沿岸长期验潮站以外的26个中短期潮位观测站评估TPXO等海潮模型预报精度。结果表明，全球海潮模型对我国沿海M2分潮的预报精度普遍较低，且主导了几种海潮模型在中国海域的整体预报精度；相比MIKE Global Tide和TPXO7.2,TPXO8、TPXO_Yellow Sea 2010和TPXO_China&Ind模型在我国沿海的预报精度更高。     

Modeling the ocean circulation in the Bering Sea

With parameterized wave mixing, the circulation and the tidal current in the Bering Sea were simulated simultaneously using the three-dimensional Princeton Ocean Model. The simulated circulation pattern in the deep basin is relatively stable, cyclonic, and has little seasonal change. The Bering Slope Current between 200-1000 m isobaths was estimated to be 5 Sv in volume transport. The Kamchatka Current was estimated to be 20 Sv off the Kamchatka Peninsula. The Bering shelf circulations vary with season, driven mainly by wind. These features are consistent with historical esti- mates. A counter current was captured flowing southeastward approximately along the 200 m isobath of the Bering Slope, opposite to the northwestward Bering Slope Current, which needs to be validated by observations. An upwelling current is located in the shelf break （ 120-1000 m） area, which may imply the vertical advection of nutrients for supporting the Bering Sea Green Belt seasonal plankton blooms in the breakslope area. The Bering Slope Current is located in a downwelling area.     

Estimates of global M_2 internal tide energy fluxes using TOPEX/POSEIDON altimeter data

TOPEX/POSEIDON altimeter data from October 1992 to June 2002 are used to calculate the global barotropic M2 tidal currents using long-term tidal harmonic analysis. The tides calculated agree well with ADCP data obtained from the South China Sea (SCS). The maximum tide velocities along the semi-major axis and semi-minor axis can be computed from the tidal ellipse. The global distribution of M2 internal tide vertical energy flux from the sea bottom is calculated based on a linear internal wave generation model. The global vertical energy flux of M2 internal tide is 0.96 TW, with 0.36 TW in the Pacific, 0.31 TW in the Atlantic and 0.29 TW in the Indian Ocean, obtained in this study. The total horizontal energy flux of M2 internal tide radiating into the open ocean from the lateral boundaries is 0.13 TW, with 0.06 TW in the Pacific, 0.04TW in the Atlantic, and 0.03 TW in the Indian Ocean. The result shows that the principal lunar semi-diurnal tide M2 provides enough energy to maintain the large-scale thermohaline circulation of the ocean.     

MEAN SEA LEVEL AND SEA SURFACE VARIABILITY OF NQRTHWEST PACIFIC OCEAN AND EASTERN CHINA SEAS FROM GEOSAT ALTIMETRY

Collinear analysis technique is widely used for determining sea surface variability with Geosat altimeterdata from its Exact Repeat Mission(ERM).But most of the researches have been only on global scaleor in oceans deeper than 2000 m.In shallow shelf waters this method is hampered by the inaccuracy ofocean tide data supplied with Geosat Geophysical Data Records(GDRs).This work uses a modified collinearanalysis technique characterized by simultaneous separation of mean sea level and ocean tide with theleast squares method,to compute sea surface variability in the Northwest Pacific Ocean and eastern ChinaSeas.The mean sea level map obtained contains not ouly bathymetric but also dynamic features such asamphidromes,indicating considerable improvement over previous works.Our sea surface variability mapsshow clearly the main current system,the well-known Zhejiang coastal upwelling,and a northern East Chi-na Sea meso-scale eddy in good agreement with satellite sea surface temperature(SST)observation and his-to     

Impacts of SIS and CICE as Sea Ice Components in BCC_CSM on the Simulation of the Arctic Climate

Sea ice is a sensitive indicator of climate change and an important component of climate system models. The Los Alamos Sea Ice Model 5.0(CICE5.0) was introduced to the Beijing Climate Center Climate System Model(BCC_CSM) as a new alternative to the Sea Ice Simulator(SIS). The principal purpose of this paper is to analyze the impacts of these two sea ice components on simulations of basic Arctic sea ice, atmosphere, and ocean states. Two sets of experiments were conducted with the same configurations except for the sea ice component used, i.e., SIS and CICE. The distributions of sea ice concentration and thickness reproduced by the CICE simulations in both March and September were closer to actual observations than those reproduced by SIS simulations, which presented a very thin sea ice cover in September. Changes in sea ice conditions also brought about corresponding modifications to the atmosphere and ocean circulation. CICE simulations showed higher agreement with the reference datasets than did SIS simulations for surface air temperature, sea level pressure, and sea surface temperature in most parts of the Arctic Ocean. More importantly, compared with simulations with SIS, BCC_CSM with CICE revealed stronger Atlantic meridional overturning circulation(AMOC), which is more consistent with actual observations. Thus, CICE shows better performance than SIS in BCC_ CSM. However, both components demonstrate a number of common weaknesses, such as overestimation of the sea ice cover in winter, especially in the Nordic Sea and the Sea of Okhotsk. Additional studies and improvements are necessary to develop these components further.     

基于SWOT分析的厦门海洋文化产业发展研究

海洋文化产业作为新的产业形态是目前一些临海国家或区域发展的经济热点。拥有丰富海洋文化资源且以良好的海洋经济发展作为依托,厦门海洋文化产业发展初显成效,但过于依赖滨海旅游业、海洋文化产业人才匮乏等问题制约了厦门海洋文化产业的发展。研究认为厦门海洋文化产业今后的发展应从打造闽台海洋文化人文旅游线路、发展厦门新兴海洋文化产业、提高厦门海洋文化产业附加值、培养本土海洋文化产业人才等4个方面需求突破。     

Estimates of global M 2 internal tide energy fluxes using TOPEX/POSEIDON altimeter data

TOPEX/POSEIDON altimeter data from October 1992 to June 2002 are used to calculate the global barotropic M ₂ tidal currents using long-term tidal harmonic analysis. The tides calculated agree well with ADCP data obtained from the South China Sea (SCS). The maximum tide velocities along the semi-major axis and semi-minor axis can be computed from the tidal ellipse. The global distribution of M ₂ internal tide vertical energy flux from the sea bottom is calculated based on a linear internal wave generation model. The global vertical energy flux of M ₂ internal tide is 0.96 TW, with 0.36 TW in the Pacific, 0.31 TW in the Atlantic and 0.29 TW in the Indian Ocean, obtained in this study. The total horizontal energy flux of M ₂ internal tide radiating into the open ocean from the lateral boundaries is 0.13 TW, with 0.06 TW in the Pacific, 0.04TW in the Atlantic, and 0.03 TW in the Indian Ocean. The result shows that the principal lunar semi-diurnal tide M ₂ provides enough energy to maintain the large-scale thermohaline circulation of the ocean. Supported by the National Basic Research Program of China (973 Program, No. 2005CB422303), the International Cooperation Program (No. 2004DFB02700), and the National Natural Science Foundation of China (No. 40552002). The TOPEX/POSEIDON data are provided by Physical Oceanography Distributed Active Archive Center (PO DACC)     

Seasonal variability of zonal heat advection in the mixed layer of the tropical Pacific

Zonal heat advection （ZHA） plays an important role in the variability of the thermal structure in the tropical Pacific Ocean, especially in the western Pacific warm pool （WPWP）. Using the Simple Ocean Data Assimilation （SODA） Version 2.02/4 for the period 1958-2007, this paper presents a detailed analysis of the climatological and seasonal ZHA in the tropical Pacific Ocean. Climatologically, ZHA shows a zonal- band spatial pattern associated with equatorial currents and contributes to forming the irregular eastern boundary of the WPWP （EBWP）. Seasonal variation of ZHA with a positive peak from February to July is most prominent in the Nifio3.4 region, where the EBWP is located. The physical mechanism of the seasonal cycle in this region is examined. The mean advection of anomalous temperature, anomalous advection of mean temperature and eddy advection account for 31%, 51%, and 18% of the total seasonal variations, respectively. This suggests that seasonal changes of the South Equatorial Current induced by variability of the trade winds are the dominant contributor to the anomalous advection of mean temperature and hence, the seasonality of ZHA. Heat budget analysis shows that ZHA and surface heat flux make comparable contributions to the seasonal heat variation in the Nifio3.4 region, and that ZHA cools the upper ocean throughout the calendar year except in late boreal spring. The connection between ZHA and EBWP is further explored and a statistical relationship between EBWP, ZHA and surface heat flux is established based on least squares fitting.     

�й��غ���ͬ����ģ�͵���б���ɷ���

???????????????????6???????????CSR4.0??EOT11a??FES2004??GOT4.7??NAO99b??TPXO7.2?????й????????????????????????????????????????б????Ч????????????????RMS???????RSS????????????????????б????????10^-8??rad??????????10^-7??rad?????????й??????????????????????????????????????????????     

Fast and slow responses of the North Pacific mode water and Subtropical Countercurrent to global warming

Six coupled general circulation models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) are employed for examining the full evolution of the North Pacific mode water and Subtropical Countercurrent (STCC) under global warming over 400 years following the Representative Concentration Pathways (RCP) 4.5. The mode water and STCC first show a sharp weakening trend when the radiative forcing increases, but then reverse to a slow strengthening trend of smaller magnitude after the radiative forcing is stablized. As the radiative forcing increases during the 21st century, the ocean warming is surface-intensified and decreases with depth, strengthening the upper ocean’s stratification and becoming unfavorable for the mode water formation. Moving southward in the subtropical gyre, the shrinking mode water decelerates the STCC to the south. After the radiative forcing is stabilized in the 2070s, the subsequent warming is greater at the subsurface than at the sea surface, destabilizing the upper ocean and becoming favorable for the mode water formation. As a result, the mode water and STCC recover gradually after the radiative forcing is stabilized.     

��ͬģ���º������ɶԸ߾���GPS������Ӱ��

???й?????????IGS????,??????FES2004????????NAO99b?????????????????????????????????λ?????????????????????????????????λ????????cm??????????????U????????????????????????3~4????????????????????????????????mm????????????????????????С????????????????GPS????????????????10^-8??????????????????????????????????С??mm?????????????????????????????1 cm??     

中国沿海不同海潮模型的倾斜负荷分析

采用积分格林函数方法及6种全球海潮模型（CSR4.0,EOT11a,FES2004,GOT4.7,NAO99b和TPXO7.2）和中国近海潮汐资料,计算了我国沿海大地控制点上的海潮倾斜负荷效应。通过标准差、均方根RMS及和方根RSS等综合分析表明,海潮倾斜负荷普遍为10^-8rad量级,最大达10^-7rad量级。在中国沿海区域,各模型差异较大,应针对不同区域采用更适合本区的模型计算。     

利用PPP反演海潮负荷位移参数

选取高纬地区的3个CORS站以及中纬地区的SHAO测站,利用PPP浮点解计算测站未经海潮负荷位移改正的坐标时间序列,采用傅里叶变换方法提取分潮的周期信息,进而反演其中8个主要分潮的海潮负荷位移参数,将其与5个全球海潮模型DTU10、EOT11a、GOT4.8、TPXO.7.2、FES2004计算得到的海洋负荷位移参数比较。结果表明,PPP反演得到的8个主要分潮的海潮负荷位移参数与海潮模型计算结果具有较好的一致性,两者均方根差异为mm级,说明利用PPP反演海潮负荷位移参数是可行的。