A coupled ice-ocean model for the Bohai Sea Ⅱ. Case study

The coupled ice-ocean model for the Bohai Sea is used for simulating the freezing, melting, and variation of ice cover and the heat balance at the sea-ice, air-ice, and air-sea interfaces of the Bohai Sea during the entire winter in 1998-1999 and 2000-2001. The coupled model is forced by real time numerical weather prediction fields. The results show that the thermodynamic effects of atmosphere and ocean are very important for the evolvement of ice in the Bohai Sea, especially in the period of ice freezing and melting. Ocean heat flux plays a key role in the thermodynamic coupling. The simulation also presents the different thermodynamic features in the ice covered region and the marginal ice zone. Ice thickness, heat budget at the interface, and surface sea temperature, etc. between the two representative points are discussed.     

A coupled ice-ocean model for the Bohai Sea Ⅱ. Case study

The coupled ice-ocean model for the Bohai Sea is used for simulating the freezing, melting, and variation of ice cover and the heat balance at the sea-ice, air-ice, and air-sea interfaces of the Bohai Sea during the entire winter in 1998～1999 and 2000～2001. The coupled model is forced by real time numerical weather prediction fields. The results show that the thermodynamic effects of atmosphere and ocean are very important for the evolvement of ice in the Bohai Sea, especially in the period of ice freezing and melting. Ocean heat flux plays a key role in the thermodynamic coupling. The simulation also presents the different thermodynamic features in the ice covered region and the marginal ice zone. Ice thickness, heat budget at the interface, and surface sea temperature, etc. between the two representative points are discussed.     

Application of an ice-ocean coupled model to Bohai Sea ice simulation

The HAMSOM(Hamburg Shelf Ocean Model),a high-resolution regional ice-ocean coupled model,was applied to investigate the seasonal evolution of Bohai Sea ice for winter 2015/2016.HAMSOM was initialized with monthly climatological temperature and salinity data from WOA13 and driven by hourly meteorological data obtained from the NCEP above the sea surface and tides at the open boundary.The ice model used here is a modifi ed Hibler-type dynamic-thermodynamic sea ice model based upon viscous-plastic rheology.The ice extent,concentration,area,thickness,length of ice season as well as the distance between the top of Liaodong Bay(North China)and the outer ice edge line were simulated and compared with the observed data.Three types of modeling experiments were carried out to investigate the eff ects of wind,tide,and both wind and tide on Bohai Sea ice.The results show that wind,as both a dynamic and a thermodynamic factor,has a signifi cant impact on the ice thickness,ice area,and ice-freezing and ice breakup dates as well as the ice velocity,while tides are a dynamic factor that infl uences only the ice velocity.During the severe ice period,the wind speed intensity increased by 25%,the average ice thickness thickened by approximately 4.0 cm in Liaodong Bay,approximately 2.1 cm in Bohai Bay and approximately 2.5 cm in Laizhou Bay,and the total ice coverage area and total ice actual area increased by about 2×104 km 2 and 1.4×104 km 2,respectively.While the tidal amplitude intensity increased by 25%,the average ice velocity increased by approximately 0.1 m/s.     

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.     

A coupled ice-ocean model for the Bohai Sea Ⅱ. Case study

1 Introduction As mentioned in the previous paper (Su et al., 2004), investigation and application of the coupled ice- ocean model play an important role in the polar and climate research. In recent years, great attention is paid to the research on coupled models of ice margin zone and regional seas and its application in marine engineering and forecasting. For example, simulation study on seasonal evolution rule of sea ice and upper layer ocean during a winter in the Labrador Sea was conduct…     

Observations and modeling of the ice-ocean conditions in the coastal Chukchi and Beaufort Seas

The Chukchi and Beaufort Seas include several important hydrological features： inflow of the Pacific water, Alaska coast current （ ACC ）, the seasonal to perennial sea ice cover, and landfast ice ＇along the Alaskan coast. The dynamics of this coupled ice-ocean system is important for both regional scale oceanography and large-scale global climate change research. A mumber of moorings were deployed in the area by JAMSTEC since 1992, and the data revealed highly variable characteristics of the hydrological environment. A regional high-resolution coupled ice-ocean model of the Chukchi and Beaufort Seas was established to simulate the ice-ocean environment and unique seasonal landfast ice in the coastal Beaufort Sea. The model results reproduced the Beaufort gyre and the ACC. The depthaveraged annual mean ocean currents along the Beaufort Sea coast and shelf hreak compared well with data from four moored ADCPs, but the simulated velocity had smaller standard deviations, which indicate small-scale eddies were frequent in the region. The model resuits captured the sea,real variations of sea ice area as compared with remote sensing data, and the simulated sea ice velocity showed an ahnost stationary area along the Beaufort Sea coast that was similar to the observed landfast ice extent. It is the combined effects of the weak oceanic current near the coast, a prevailing wind with an onshore component, the opposite direction of the ocean current, and the blocking hy the coastline that make the Beaufort Sea coastal areas prone to the formation of landfast ice.     

Simulation of Typhoon Muifa using a mesoscale coupled atmosphere–ocean model

A mesoscale coupled atmosphere–ocean model has been developed based on the GRAPES(Global and Regional Assimilation and Prediction System) regional typhoon model(GRAPES_TYM) and ECOM-si(estuary, coast and ocean model(semi-implicit)). Coupling between the typhoon and ocean models was conducted by exchanging wind stress, heat, moisture fluxes, and sea surface temperatures(SSTs) using the coupler OASIS3.0. Numerical prediction experiments were run with and without coupling for the case of Typhoon Muifa in the western North Pacific. To investigate the impact of using more accurate SST information on the simulation of the track and the intensity of Typhoon Muifa, experiments were also conducted using increased SST resolution in the initial condition field of the control test. The results indicate that increasing SST resolution in the initial condition field somewhat improved the intensity forecast, and use of the coupled model improved the intensity forecast significantly, with mean absolute errors in maximum wind speed within 48 and 72 h reduced by 32% and 20%, respectively. Use of the coupled model also resulted in less pronounced over-prediction of the intensity of Typhoon Muifa by the GRAPES_TYM. Moreover, the effects of using the coupled model on the intensity varied throughout the different stages of the development of Muifa owing to changes in the oceanic mixed layer depth. The coupled model had pronounced effects during the later stage of Muifa but had no obvious effects during the earlier stage. The SSTs predicted by the coupled model decreased by about 5–6°C at most after the typhoon passed, in agreement with satellite data. Furthermore, based on analysis on the sea surface heat flux, wet static energy of the boundary layer, atmospheric temperature, and precipitation forecasted by the coupled model and the control test, the simulation results of this coupled atmosphere–ocean model can be considered to reasonably reflect the primary mechanisms underlying the interactions between tropical cyclones and oceans.     

Numerical Simulation of Bohai Oil Spill in the Winter Sea Ice Period

The Bohai Sea is a seasonal icy sea area that has the lowest latitude of any sea experiencing icing in the northern hemisphere, and simulation studies on oil spills during its sea ice period are the key to analyzing winter oil spill accidents. This study applied the three-dimensional free surface to establish a high-resolution hydrodynamic model and simulate tidal distributions in the Bohai Sea. Then, the oil spill model of the open sea area and thermodynamic model were combined to establish a numerical model for the Bohai oil spill during the winter sea ice period. The hydrodynamic model and sea ice growth and melting model were verified, and the parameters were adjusted based on the measured values, which indicate that the numerical model established in this paper is of high accuracy,stability and ubiquity. Finally, after checking the calculations repeatedly, the diffusion coefficient for the Bohai Sea was determined to be 1.0×10~(–7 )m~2/s. It is better that the comprehensive weathering attenuation coefficient is lower than that of a non-winter oil spill, with 1.3×10~(–7 )m~2/s being the most appropriate coefficient. This study can provide the reliable technical support for the operational safety and reduction in losses caused by winter oil spill accidents for the petroleum industry.     

A comparison of two global ocean-ice coupled models with different horizontal resolutions

A global eddy-permitting ocean-ice coupled model with a horizontal resolution of 0.25 by 0.25 is established on the basis of Modular Ocean Model version 4 (MOM4) and Sea Ice Simulator (SIS). Simulation results are compared with those of an intermediate resolution ocean-ice coupled model with a horizontal resolution of about 1 by 1 . The results show that the simulated ocean temperature, ocean current and sea ice concentration from the eddy-permitting model are better than those from the intermediate resolution model. However, both the two models have the common problem of ocean general circulation models (OGCMs) that the majority of the simulated summer sea surface temperature (SST) is too warm while the majority of the simulated subsurface summer temperature is too cold. Further numerical experiments show that this problem can be alleviated by incorporating the non-breaking surface wave-induced vertical mixing into the vertical mixing scheme for both eddy-permitting and intermediate resolution models.     

A numerical study of the South China Sea deep circulation and its relation to the Luzon Strait transport

A fine-resolution MOM code is used to study the South China Sea basin-scale circulationand its relation to the mass transport through the Luzon Strait. The model domain includes the South China Sea, part of the East China Sea, and part of the Philippine Sea so that the currents in the vicinity of the Luzon Strait are free to evolve. In addition, all channels between the South China Sea and the Indonesian seas are closed so that the focus is on the Luzon Strait transport. The model is driven by specified Philippine Sea currents and by surface heat and salt flux conditions. For simplicity, no wind-stress is applied at the surface.The simulated Luzon Strait transport and the South China Sea circulation feature a sandwich vertical structure from the surface to the bottom. The Philippine Sea water is simulated to enter the South China Sea at the surface and in the deep ocean and is carried to the southern basin by western boundary currents. At the intermediate depth, the net Luzon Strait transport is out of t     

The estimate of sea ice resources quantity in the Bohai Sea based on NOAA/AVHRR data

The research on sea ice resources is the academic base of sea ice exploitation in the Bohai Sea. According to the ice-water spectrum differences and the correlation between ice thickness and albedo, this paper comes up with a sea ice thickness inversion model based on the NOAA/AVHRR data. And then a sea ice resources quantity (SIQ) time series of Bohai Sea is established from 1987 to 2009. The results indicate that the average error of inversion sea ice thickness is below 30%. The maximum sea ice resources quantity is about 6 × 10 9 m 3 and the minimum is 1.3 × 10 9 m 3 . And a preliminary analysis has been made on the errors of the estimate of sea ice resources quantity (SIQ).     

基于高光谱遥感的渤海海冰厚度半经验模型

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness(SEMSIT) that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands(spectral dimension information) are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.     

A PARAMETRIC MODEL FOR THERMAL STRUCTURE FEATURES OF THE OCEAN UPPER LAYER

Based on the non-dimensional general function for the thermal structure features presented by S. A. Kitaigorodsky et al. (1963, 1965)[10,11]. this paper tries to establish semi-empirical and semi-theoretical models bfor the thickness of the upper homogeneous layer of the ocean, thermocline intensity and lower oundary depth of thermocline by consecutive observations from 159 stations in the Bohai Sea,Huanghai Sea and East China Sea in the warm half of the years from 1957 to 1964 and the heat budget at the sea surface computed with the simplified computing formulae proposed by Wang (1983)[22]. This model indicates the main factors forming the thermal structure features in the upper layer of the ocean and their function. With the model, one can directly use the sea surface temperature, air temperature and wind speed to compute the thermal feature in the upper layer of the ocean.     

OSTIA数据在中国近海业务化环流模型中的同化应用

The prediction of sea surface temperature(SST) is an essential task for an operational ocean circulation model. A sea surface heat flux, an initial temperature field, and boundary conditions directly affect the accuracy of a SST simulation. Here two quick and convenient data assimilation methods are employed to improve the SST simulation in the domain of the Bohai Sea, the Yellow Sea and the East China Sea(BYECS). One is based on a surface net heat flux correction, named as Qcorrection(QC), which nudges the flux correction to the model equation; the other is ensemble optimal interpolation(En OI), which optimizes the model initial field. Based on such two methods, the SST data obtained from the operational SST and sea ice analysis(OSTIA) system are assimilated into an operational circulation model for the coastal seas of China. The results of the simulated SST based on four experiments, in 2011, have been analyzed. By comparing with the OSTIA SST, the domain averaged root mean square error(RMSE) of the four experiments is 1.74, 1.16, 1.30 and 0.91°C, respectively; the improvements of assimilation experiments Exps 2, 3 and 4 are about 33.3%, 25.3%, and 47.7%, respectively.Although both two methods are effective in assimilating the SST, the En OI shows more advantages than the QC,and the best result is achieved when the two methods are combined. Comparing with the observational data from coastal buoy stations, show that assimilating the high-resolution satellite SST products can effectively improve the SST prediction skill in coastal regions.     

渤海模型冰运动特征的现场观测研究

Sea ice drift is mainly controlled by ocean currents, local wind, and internal ice stress. Information on sea ice motion, especially in situ synchronous observation of an ice velocity, a current velocity, and a wind speed, is of great significance to identify ice drift characteristics. A sea ice substitute, the so-called "modelled ice", which is made by polypropylene material with a density similar to Bohai Sea ice, is used to complete a free drift experiment in the open sea. The trajectories of isolated modelled ice, currents and wind in the Bohai Sea during non-frozen and frozen periods are obtained. The results show that the currents play a major role while the wind plays a minor role in the free drift of isolated modelled ice when the wind is mild in the Bohai Sea. The modelled ice drift is significantly affected by the ocean current and wind based on the ice–current–wind relationship established by a multiple linear regression. The modelled ice velocity calculated by the multiple linear regression is close to that of the in situ observation, the magnitude of the error between the calculated and observed ice velocities is less than12.05%, and the velocity direction error is less than 6.21°. Thus, the ice velocity can be estimated based on the observed current velocity and wind speed when the in situ observed ice velocity is missing. And the modelled ice of same thickness with a smaller density is more sensitive to the current velocity and the wind speed changes. In addition, the modelled ice drift characteristics are shown to be close to those of the real sea ice, which indicates that the modelled ice can be used as a good substitute of real ice for in situ observation of the free ice drift in the open sea, which helps solve time availability, safety and logistics problems related to in situ observation on real ice.     

NUMERICAL EXPERIMENTS OF SEA ICE IN THE BOHAI SEA

A rather complete sea ice model is given, which deals with not only thermodynamic and dynamic processes commonly used in previous models of sea ice but also a melting process of ice driven into warmer waters. A series of numerical experiments have been carried out in order to search after a mechanism of the growth and decay of sea ice in the Bohai sea, and the principal result shows that the melting process of sea ice driven into the warmer waters must be taken into consideration when the ice condition in such a partially frozen sea as the Bohai Sea is calculated.     

The long-term trend of Bohai Sea ice in different emission scenarios

Based on a coupled ocean-sea ice model, this study investigates how changes in the mean state of the atmosphere in different CO_2 emission scenarios(RCP 8.5, 6.0, 4.5 and 2.6) may affect the sea ice in the Bohai Sea, China,especially in the Liaodong Bay, the largest bay in the Bohai Sea. In the RCP 8.5 scenario, an abrupt change of the atmospheric state happens around 2070. Due to the abrupt change, wintertime sea ice of the Liaodong Bay can be divided into 3 periods: a mild decreasing period(2021–2060), in which the sea ice severity weakens at a nearconstant rate; a rapid decreasing period(2061–2080), in which the sea ice severity drops dramatically; and a stabilized period(2081–2100). During 2021–2060, the dates of first ice are approximately unchanged, suggesting that the onset of sea ice is probably determined by a cold-air event and is not sensitive to the mean state of the atmosphere. The mean and maximum sea ice thickness in the Liaodong Bay is relatively stable before 2060, and then drops rapidly in the following decade. Different from the RCP 8.5 scenario, atmospheric state changes smoothly in the RCP 6.0, 4.5 and 2.6 scenarios. In the RCP 6.0 scenario, the sea ice severity in the Bohai Sea weakens with time to the end of the twenty-first century. In the RCP 4.5 scenario, the sea ice severity weakens with time until reaching a stable state around the 2070 s. In the RCP 2.6 scenario, the sea ice severity weakens until the2040 s, stabilizes from then, and starts intensifying after the 2080 s. The sea ice condition in the other bays of the Bohai Sea is also discussed under the four CO_2 emissions scenarios. Among atmospheric factors, air temperature is the leading one for the decline of the sea ice extent. Specific humidity also plays an important role in the four scenarios. The surface downward shortwave/longwave radiation and meridional wind only matter in certain scenarios, while effects from the zonal wind and precipitation are negligible.     

Simulation of the seasonal thermal structure in the Bohai Sea

The seasonal thermal structure in the Bohai Sea are examined with a three-dimensional boroclinic primitive equation model for shelf sea.The evolution of the seasonal thermal stratification is well simulated.The stratification appears early in April,first in the area off Qinhuangdao and it is well developed in the middle of May.It intensifies with synoptic and neap-spring fluctuations throughout the summer and reaches its maximum in the middle of July.Eventually,it is destroyed at the end of September.There are cold water belts between well-mixed and stratified regions.They are loGated on the mixed side of tidal fronts,and coincide with the isolines for a temperature difference of 1-2℃ between surface and bottom.The sea surface temperature (SST) distribution shows local maxima at the head of three bays and to the south of Qinhuangdao during the summer.The Bohai Sea responds to the variability in the atmospheric forcing and in tides with the synoptic and neap-spring variations of SST,as well as in the stratification and in variable positions of tidal fronts.     

基于MODIS热红外数据的渤海海冰厚度反演

Level ice thickness distribution pattern in the Bohai Sea in the winter of 2009–2010 was investigated in this paper using MODIS night-time thermal infrared imagery.The cloud cover in the imagery was masked out manually.Level ice thickness was calculated using MODIS ice surface temperature and an ice surface heat balance equation.Weather forcing data was from the European Centre for Medium-Range Weather Forecasts(ECMWF) analyses.The retrieved ice thickness agreed reasonable well with in situ observations from two off-shore oil platforms.The overall bias and the root mean square error of the MODIS ice thickness are –1.4 cm and 3.9 cm,respectively.The MODIS results under cold conditions(air temperature –10°C) also agree with the estimated ice growth from Lebedev and Zubov models.The MODIS ice thickness is sensitive to the changes of the sea ice and air temperature,in particular when the sea ice is relatively thin.It is less sensitive to the wind speed.Our method is feasible for the Bohai Sea operational ice thickness analyses during cold freezing seasons.     

Effect of Langmuir circulation on upper ocean mixing in the South China Sea

Effect of Langmuir circulation (LC) on upper ocean mixing is investigated by a two-way wave-current coupled model. Themodel is coupled of the ocean circulationmodel ROMS (regional ocean modeling system) to the surface wave model SWAN (simulating waves nearshore) via the model-coupling toolkit. The LC already certified its importance by many one-dimensional (1D) research andmechanismanalysis work. This work focuses on inducing LC’s effect in a three-dimensional (3-D) model and applying it to real field modeling. In ROMS, theMellor-Yamada turbulence closuremixing scheme is modified by including LC’s effect. The SWAN imports bathymetry, free surface and current information fromthe ROMS while exports significant wave parameters to the ROMS for Stokes wave computing every 6 s. This coupled model is applied to the South China Sea (SCS) during September 2008 cruise. The results show that LC increasing turbulence and deepening mixed layer depth (MLD) at order of O (10 m) in most of the areas, especially in the north part of SCS where most of our measurements operated. The coupled model further includes wave breaking which will bringsmore energy into water. When LC works together with wave breaking,more energy is transferred into deep layer and accelerates the MLD deepening. In the north part of the SCS, their effects aremore obvious. This is consistent with big wind event in the area of the Zhujiang River Delta. The shallow water depth as another reasonmakes themeasy to influence the oceanmixing as well.