Princeton Ocean Model Research Articles

Applying satellite sea surface temperature as Dirichlet-type surface thermal boundary condition in an ocean model

Following the Dirichlet-type boundary condition that specifies value at boundary of a system in mathematics and physics, this study suggests applying satellite sea surface temperature (SST) as Dirichlet-type surface thermal boundary condition (STBC) in ocean models. Numerical experiments with Dirichlet-type and Combined-type STBC with different configurations of satellite SST, for the period of January to April 2019 in the Northwest Pacific Ocean, were carried out based on the Princeton Ocean Model. The experiments were assessed using satellite and in situ observations of temperature, GOFS3.1 analysis and GLORYS12v1 reanalysis. Results show that applying satellite SST as Dirichlet-type STBC could constrain the modeled SST and near-surface temperature in upper 50 m depth well, which is better than the scheme that uses satellite SST as the relaxation terms of Combined-type STBC. The temperature section along the 137°E resulted from the Dirichlet-type STBC is comparable with GOFS3.1 analysis and better than GLORYS12v1 reanalysis in upper 50 m depth. These results suggest that applying high-accuracy satellite SST as Dirichlet-type STBC in ocean models has a promising prospect in numerical simulation.

Detecting the role of Stokes drift under typhoon condition by a fully coupled wave-current model

The impacts of Stokes drift and sea-state-dependent Langmuir turbulence (LT) on the three-dimensional ocean response to a tropical cyclone in the Bohai Sea are studied through two-way coupled wave-current simulations. The Stokes drift is calculated from the simulated wave spectrum of the wave model, Simulating Waves Nearshore (SWAN), and then input to the Princeton Ocean Model with the generalized coordinate system (POMgcs) to represent the Langmuir effect. The Langmuir circulation is included in the vertical mixing of the ocean model by adding the Stokes drift to the shear of the vertical mean current and by including LT enhancements to the Mellor-Yamada 2.5 turbulent closure submodel. Simulations are assessed through the case study of Typhoon Masta in 2005 with a set of diagnostic experiments that incorporated different terms of Stokes production (SP) respectively. It is shown that with the consideration of SP, a deeper mixed layer, an enhanced vertical mixing coefficient KMS, and a more accurate representation of the vertical temperature distribution could be derived. Moreover, the effect of LT in elevating the turbulence mixing is stronger than that of Coriolis Stokes force (CSF) and Craik-Leibovich vortex force (CLVF). LT has a greater influence on the vertical mixing during typhoon than that in normal weather.

Oil Spill Sensitivity Analysis of the Coastal Waters of Taiwan Using an Integrated Modelling Approach

Pollution caused by marine oil spills can lead to persistent ecological disasters and severe social and economic damages. Numerical simulations are useful and essential tools for accurate decision making during emergencies and planning response actions. In this study, we applied the Princeton Ocean Model (POM) to determine current data, including seawater velocity, salinity, and temperature, and we obtained the fate and trajectory of spilled oil using OpenOil. Several probable oil slicks around Taiwan were simulated over time (12 months) and space (four spill locations in the marine area of each coastal city or county) using the model. The percentage risk under the effect of an oil spill is estimated. The risk zone of the coastal waters of Taiwan was identified based on the frequency of simulated oil slicks hitting the coast and sensitive resources. This information not only helps authorities guide the preparation of effective plans to minimise the impacts of oil spill incidents but could also be used to improve regulations related to shipping and vessel navigation in regional seas.

Tidal impacts on the discharge division and freshwater transport in the Berau Delta, East Kalimantan, Indonesia

Flow division at bifurcations is important in controlling material carried from the terrestrial to the coastal zone in tidally influenced deltas. The present study aims to identify the impacts of tides on flow division and freshwater transport in the Berau Delta, East Kalimantan. The Princeton Ocean Model (POM) is applied to simulate the hydrodynamics with forcing from river discharges and tides. Tides at open seas and observed river discharge at upstream locations were used to set model boundaries. Model validation was accomplished by comparing measurements of tidal elevation and salinity time series with model results. The model results reproduce the observed in temporal variations of tidal elevation and salinity. Model results highlight that the tidal amplitude has large influence on discharge division at the tidal junction. Tides enhance equal subtidal flow distribution in the river junctions, with the influence of tide is 15% on subtidal flow distribution in the bifurcations. Freshwater discharge at the Berau Delta flows mainly to the north and middle outlet with 90% of the total discharge. Based on an analysis model results, the decomposition of freshwater transport reveals that advection and tidal pumping are major factor, which are fluctuate in fortnight cycle.

Temperature structure and wind-induced water level anomalies in the Caspian Sea: A study on coastal upwelling

In this study, the index of coastal upwelling was examined based on wind-induced water level anomaly using the Princeton ocean model (POM) in the Caspian Sea. The data of wind fields and atmospheric fluxes were obtained from European Centre for Medium-Range Weather Forecasts (ECMWF) database and considering freshwater inflows in the Volga, Kura, and Ural Rivers. This model was implemented for 10 years, from 2008 to 2017. Findings indicated that winds, currents, water levels, and sea surface temperatures are highly correlated in upwelling events. Ekman transport and surface geostrophic current (southward) occur due to wind field and water level anomaly (about 10 cm, between the coastal and offshore areas) in the eastern shores of the middle basin, respectively. Partial and full coastal upwelling mainly occur in the cold and warm months, respectively, in the middle basin (upwelling center). So, in June, July, and August, on the eastern shores of the middle and southern basins, a greater number of water level declines occurred, and the upwelling index increased compared to other months. Topography and the weak wind field effects can be the mechanism causing the occurrence of partial coastal upwelling. Also, the findings implied that seawater temperature anomalies are higher in the windy state in both June and January compared to calm situations. The intensity of sea current in the windy state in the warm months is less than in the cold months due to the weak counter-clockwise circulation.

Hindcast of oil spill pollution in the East China Sea in January 2018

In this study, an integrated model is developed for studying the Sanchi oil spill event, which occurred in the East China Sea in January 2018. The results of the Advanced Research Weather Research and Forecasting model (WRF-ARW) as well as the Princeton Ocean Model (POM) are used for meteorological forecasting and the hydrodynamic simulations, respectively. These data are adopted as inputs for the OpenOil, a sub-module of OpenDrift, for the oil spill model. Some reference experiments are examined for short-term hindcast. The satellite image is used to validate the numerical result. The oil slicks of the satellite image and the numerical result are of similar shapes. Quantitatively, the simulated oil slick and that from the satellite image are located closely and have similar dimensions of 56 km by 34 km and 54 km by 29 km, respectively. It is found that the accurate results can be obtained by the proposed integrated model with the high-frequency (hourly) and high-spatial-resolution data as inputs, and the wind drift factor has to be added. The long-term 1-month simulation showed that most of the oil particles would move to the northeast of the sinking location and be trapped by the Kuroshio current.

An operational improvement of A-4DEnVar and its application to the estimation of the spatially varying bottom friction coefficients of the M2 constituent in the Bohai and Yellow seas

The analytical four-dimensional ensemble variational (A-4DEnVar) data assimilation scheme inherits the advantages of the conventional four-dimensional variational (4D-Var) data assimilation scheme and removes the adjoint model. However, compatible operational improvements such as the reduction of the computational costs and the localization method should be considered when it is used in realistic systems. In this paper, the computational complexity of calculating the inverse of background error covariance (the B matrix) is decreased by a precondition transform method, i.e., introducing a new state variable whose product with the B matrix is the original state variable to be optimized in the cost function. Furthermore, an independent point (IP) scheme is combined to construct an implicit localization method and further decreases the computational cost. Based on the Princeton Ocean Model with the generalized coordinate system (POMgcs), the operational improved A-4DEnVar is applied to optimize the spatially varying bottom friction coefficients (BFCs) of the M2 constituent in the Bohai and Yellow seas. A twin experiment with idealized observations is designed to validate the effectiveness of the proposed method. In practical experiments, with no more than 10 IPs, the algorithm can assimilate observations from the National Astronomical Observatory (NAO) dataset and obtain a good simulation. The experimental performances increase with the increase of either the IPs or observations, which indicates the efficacy of the proposed algorithm.

Variational parameter estimation in a two-equation turbulence model: A case study with a 3D primitive-equation ocean model

A three-dimensional and complete adjoint model of the Princeton Ocean Model with a generalized coordinate system (POMgcs) is developed to construct the 4D-Variational data assimilation (4D-Var) algorithm in this study. Uncertain parameters in the Mellor-Yamada 2.5 turbulence submodel (MY-2.5) which is enclosed in POMgcs, are tentatively estimated via the 4D-VAR algorithm within a biased model framework. Here, the control variables in the biased model are set to two uncertain wave-affected parameters (wave energy factor α and Charnock coefficient β ) in the MY-2.5 turbulence model, which play a crucial role in modulating the heat content distribution in the upper coastal sea. First of all, the ocean temperature and salinity in a typical coastal sea, Bohai Sea, are simulated by the model to validate the rationally of the MY-2.5 parameterization scheme for both constructing the “truth model” and generating the “pseudo-observations” in the data assimilation studies. Then, after thoroughly testing the ability of the 4D-Var to optimize the initial state fields of the POMgcs model, a series of parameter estimation experiments are carried out to investigate whether and to what degree the parameters embedded in high-order turbulence models can be significantly optimized. Results of parameter estimation with perfect initial fields show the two estimated parameters in the MY-2.5 submodel can successfully converge to the “truth” value. The local minimum of the cost function can be effectively and efficiently jumped out once two kinds of optimization algorithms, LBFGS and LMBM, are jointly used. In addition, the estimated parameter will converge to the optimal value rather than the truth one to compensate for the initial field error when the state-parameter are estimated simultaneously. Further, the performance of the parameter estimation is also deeply discussed when the observation samples are noised. Finally, prescribing the initial field and parameter as error source, a forecasting experiment for sea temperature is performed. The experiment results indicate that assimilating “pseudo-observations” to the model based on 4D-Var can significantly improve the sea temperature simulation. Moreover, adjusting the initial field and parameter leads to a better result than the only initial field, and this conclusion is more evident at the surface than in the deeper ocean.

Power Maximisation of Wind Energy Using Wind Speed Sensors on Stewart Island.

This paper evaluates the feasibility of using wind power for power supply to coastal communities isolated from the main supply grid. The case study is Stewart Island, where the cost of electricity provided by a central diesel power station is higher than the grid network in New Zealand. The Princeton Ocean Model (POM) conducted by MetOcean Solutions Limited (MSL) is used to find Foveaux as an optimized site for generating wind power. Global Wind Atlas is used to plot the wind rose of current wind patterns in New Zealand. In the next step, wind speed data from each site are imported from the NASA database to WRPLOT view software and Homer Pro to find wind frequency distribution and output power in the area. The maximum annual power can be seen in WSW (32,299 kW hours/year), SW (20,111 kW hours/year) and W (15,622 kW hour/year) directions, respectively.

Numerical Study of Circulation and Seasonal Variability in the Southwestern Yellow Sea

A nested-grid ocean circulation modelling system (NGMS-swYS) is used for examining the impact of tides and winds on the three-dimensional (3D) circulation, hydrography and seasonal variability over the southwestern Yellow Sea (swYS). The modelling system is based on the Princeton Ocean Model (POM) and uses a nested-grid setup, with a fine-resolution (~2.7 km) inner model nested inside a coarse-resolution (~9.0 km) outer model. The domain of the outer model covers the China Seas and adjacent deep ocean waters. The domain of the fine-resolution inner model covers the swYS and adjacent waters. The NGMS-swYS is driven by a suite of external forcings, including the atmospheric forcing, tides, freshwater discharge and currents specified at the lateral open boundaries. A comparison of model results with observations and previous numerical studies demonstrates the satisfactory performance of the NGMS-swYS in simulating tides, seasonal mean circulation and distribution of temperature and salinity. Five additional numerical experiments were conducted using NGMS-swYS with different combinations of external forcing. Analysis of model results demonstrates that the monthly mean circulation over the swYS is affected significantly by tides and winds, with large seasonal variability. The northward Subei Shoal Current occurred in both winter and summer months in 2015, with persistent strong southeastward mean currents induced by tides along the 50 m isobath. Model results also demonstrated that strong wind-induced currents occurred with large sea surface cooling during Typhoon Chan-Hom.

Coastal upwelling by wind-driven forcing in the Caspian Sea: A numerical analysis

In this study, wind-driven coastal upwelling in the Caspian Sea was investigated using a developed three-dimensional hydrodynamic numerical model based on the Princeton Ocean Model (POM). The model was forced with wind fields and atmospheric fluxes from the ECMWF database and it considers freshwater inflows from the Volga, Kura and Ural Rivers. This model was implemented for 10 years (2008–2018). Findings indicated that the upwelling in the Caspian Sea was due to effects of wind and bottom topography, often occurring from May to September. In June and July, in the eastern part of the middle and sometimes southern basins, up to 3°C water temperature difference occurs between coastal and offshore areas. The vertical temperature gradient in the middle basin was larger than that in the southern basin. Upwelling in August in the eastern coasts of the middle basin within 25 km of coast from the depth of 15 m to the surface was shown, which was due to the effects of wind and bottom topography. In the middle basin, the highest vertical velocities caused by upwelling in June, July and August were 12, 13.82, and 10.36 m/day, respectively.

Comparative performance of HWRF model coupled with POM and HYCOM for tropical cyclones over North Indian Ocean

The demonstration of air-sea interaction in the high-resolution numerical weather prediction (NWP) model is the necessary condition to predict the extreme event like tropical cyclone. It is already established that the coupled mesoscale models are capable of forecasting the track and intensity of tropical cyclones with reasonable skill over North Indian Ocean (NIO) basin and over other ocean basins as well. The Hurricane Weather Research and Forecasting (HWRF) model is consequently operational by India Meteorological Department (IMD) for the forecasting tropical cyclones like many other National Meteorological and Hydrological Services. Two different coupled versions of HWRF modeling systems with feature based Princeton Ocean Model (POM) and eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) are employed for the real-time forecasting of tropical cyclones over NIO. In this study, the comparison of forecast skills between two coupled variants in predicting track and intensity for all eight cyclones during 2019 has been carried out. The results show that the track prediction of all cyclones by HWRF model coupled with HYCOM is comparatively better than that with POM. But the intensity prediction skills of both coupled models do not portray any obvious superiority for all forecast hours. The absolute intensity errors of HWRF model with POM coupling are large compared to HYCOM during initial phase up to 60 hours of forecast except for Vayu cyclone. But, the errors in intensity exhibited by POM are a bit smaller compared to HYCOM for longer lead period. However, the standard deviation is very large for both models and numbers of forecasts verified are also less. In case of two cyclones of 2020, the ability of coupled variant with HYCOM is superior to POM in forecasting the tracks of cyclones. Considering the skills of both the models for intensity prediction, it is clear that the use of coupled HCOM model for short-lived weaker storms is advantageous but the performance of the model is not consistently maintained with longer lead time for intense and long-lived cyclone. A brief assessment of surface fluxes and Sea Surface Temperature reveals that HYCOM portray more realistic ocean condition in turn surface fluxes compared to POM and therefore the coupled variant of HWRF with HYCOM is superior to the variant with POM while forecasting cyclones over deep Ocean.

Development of a Tsunami-tide Interaction Model with Simulation of Tsunami Wave Propagation to Korea Coasts

Kim, K.O.; Jung, K.T., and Ha, T.M., 2021. Development of a tsunami-tide interaction model with simulation of tsunami wave propagation to Korea coasts. In: Lee, J.L.; Suh, K.-S.; Lee, B.; Shin, S., and Lee, J. (eds.), Crisis and Integrated Management for Coastal and Marine Safety. Journal of Coastal Research, Special Issue No. 114, pp. 91–95. Coconut Creek (Florida), ISSN 0749-0208. Simulation studies of historic and prognostic tsunamis have been performed with a focus on tsunami-tide interactions. A modeling system based on the 3D POM (Princeton Ocean Model) has been developed that is composed of a global tsunami model and a regional multinesting tide/tsunami/tsunami-tide interaction model covering the Northwest Pacific region. Propagations of two historic tsunamis, the 2010 Chilean and 2011 Tohoku tsunamis, were simulated to estimate the effects of tsunami-tide interactions by comparing the results from tsunami-only and tsunami simulations with tide simulations. Calculations show that tsunami waves are significantly transformed after traveling a long distance and subsequently affected by the tidal phase. As prognostic tsunami propagation, the Nankai Trough earthquake-induced tsunami wave in Japan was investigated. A series of calculations were performed with a systematic shift of the tsunami occurrence time. The effects of tsunami-tide interactions were examined in terms of local tidal phase changes, focusing on the changes in the first tsunami heights and occurrence times arriving at the selected tidal stations in Korea. The tsunami-tide interaction mostly results in a positive change in the tsunami first peak heights with a maximum change near low tide probably due to the interaction between the tsunami-induced current and tidal current. The water depth increase near high tide advances the arrival of the tsunami first peaks, whereas the water depth decrease near low tide delays the arrival of the tsunami first peaks.

Numerical Study of the Thermal Structure and Circulation in a Large and Deep Dimictic Lake Over Tibetan Plateau

AbstractA three‐dimensional (3‐D) hydrodynamic model based on the Princeton Ocean Model (POM) was applied to simulate the thermal structure and circulation of Lake Nam Co (LNC), the third largest lake over the Tibetan Plateau (TP), during May–December 2013. Compared with a spatially distributed set of one‐dimensional thermal diffusion lake models, POM better reproduced the observed seasonal evolution of the horizontal distribution of lake surface temperature and the vertical thermal structure. A heat budget analysis confirmed that the lateral heat exchange made significant contributions to the horizontal variability of lake temperature. The model results showed that LNC was thermally stratified in summer, had a weak inverse stratification since mid‐December, and was fully turned over during late spring and autumn. During both overturning phases, the modeled “thermal bar” was developed as a result of the density‐driven convection in response to the radiative heating (surface cooling) during spring (autumn). The 3‐D model results showed that the monthly mean circulation featured a predominant mid‐lake cyclonic gyre throughout the ice‐free period; upwelling along the western coast and strong coastal currents occurred in all months except in July–August. Model sensitivity experiments confirmed that the lake circulation was primarily driven by the barotropic dynamics of the prevailing southwesterly wind, while the baroclinic process made a secondary contribution. The results pointed out the necessity to resolve lateral processes when modeling large TP lakes.

The Cross‐Shore Component in the Vertical Structure of Wave‐Induced Currents and Resulting Offshore Transport

AbstractA comprehensive field survey was conducted, and a numerical model based on the Princeton Ocean Model was established to investigate the vertical structure of current induced by surface waves and identify its effect on sand loss from a beach. Previous studies have examined wave‐induced longshore transport, whereas this study focused on cross‐shore transport caused by waves. The introduction of vertically dependent radiation stress revealed that wave‐induced vertical gyres existed with onshore and offshore flows in the upper and lower layers, respectively, besides the well‐known wave‐induced longshore flow. It was found that the vertical gyres were caused by the onshore wave force (a gradient of radiation stress) in the upper layer and the offshore wave force usually existed in the lower layer, as well as the horizontal imbalance of local wave force. The locations of wave‐induced vertical gyres were found to correspond to the location of sediment erosion, and thus the lower part of the vertical gyres drove the offshore transport of eroded sediment. Moreover, the vertical gyres were located in areas with a highly variable topographic gradient (e.g., at a bar), and their driven offshore currents (undertow/return flow) extended far beyond the surf zone. The topographic slope had a profound effect on the vertical structure. Numerical experiments showed that an extended breakwater from the eastern headland of the bay could prevent sediment from being eroded. Furthermore, the shear front at the rim of a horizontal gyre could retain suspended sediment in the center of the gyre.

Revisiting the Ocean’s Nonisostatic Response to 5-Day Atmospheric Loading: New Results Based on Global Bottom Pressure Records and Numerical Modeling

AbstractWe use bottom pressure records from 59 sites of the global tsunami warning system to examine the nonisostatic response of the World Ocean to surface air pressure forcing within the 4–6-day band. It is within this narrow “5-day” band that sea level fluctuations strongly depart from the isostatic inverted barometer response. Numerical simulations of the observed bottom pressures were conducted using a two-dimensional Princeton Ocean Model forced at the upper boundary by two versions of the air pressure loading: (i) an analytical version having the form of the westward propagating, 5-day Rossby–Haurwitz air pressure mode; and (ii) an observational version based on a 16-yr record of global-scale atmospheric reanalysis data with a spatial resolution of 2.5°. Simulations from the two models—consisting of barotropic standing waves of millibar amplitudes and near uniform phases in the Pacific, Atlantic, and Indian Oceans—are in close agreement and closely reproduce the observed bottom pressures. The marked similarity of the outputs from the two models and the ability of both models to accurately reproduce the seafloor pressure records indicate a pronounced dynamic response of the World Ocean to nonstationary air pressure fields resembling the theoretical Rossby–Haurwitz air pressure mode.

3D Numerical Simulation of Kuroshio-induced Wake Near Green Island, Taiwan

In this research, Kuroshio-induced vortex shedding in the lee of Green Island was numerically studied using the three-dimensional Princeton Ocean Model. The temporal and spatial characteristics of the vortex street were analyzed at different tidal levels, wind fields and water depths. The nu-merical results were validated by comparing with in-situ measurements from the Green Island tidal station, shipboard acoustic doppler current profiler (ADCP), conductivi-ty-temperature-depth (CTD) profiler, and moderate resolution imaging spectroradiometer (MODIS) satellite image datasets. Also, the three-dimensional results of velocity, temperature, and salinity fields were compared well with observational data. Furthermore, the influence depths of the wind effect and vor-tex shedding were found to be approximately 100 and 300 meters, respectively. The aspect ratio and dimensionless width were evaluated to be 3.64 and 1.46, respectively. In addition, the temporal scales (Strouhal number), and propagation speeds were analyzed at various water depths. Finally, the elongation of the Kuroshio-induced vortex street was studied under wind and tide effects at various water depths.

A Nowcast/Forecast System for Japan’s Coasts Using Daily Assimilation of Remote Sensing and In Situ Data

We have developed an ocean state nowcast/forecast system (JCOPE-T DA) that targets the coastal waters around Japan and assimilates daily remote sensing and in situ data. The ocean model component is developed based on the Princeton Ocean Model with a generalized sigma coordinate and calculates oceanic conditions with a 1/36-degree (2–3 km) resolution and an hourly result output interval. To effectively represent oceanic phenomena with a spatial scale smaller than 100 km, we adopted a data assimilation scheme that explicitly separates larger and smaller horizontal scales from satellite sea surface temperature data. Our model is updated daily through data assimilation using the latest available remote-sensing data. Here we validate the data assimilation products of JCOPE-T DA using various kinds of in situ observational data. This validation proves that the JCOPE-T DA model output outperforms those of a previous version of JCOPE-T, which is based on nudging the values of temperature and salinity toward those provided by a different coarse grid data-assimilated model JCOPE2M. Parameter sensitivity experiments show that the selection of horizontal scale separation parameters considerably affects the representation of sea surface temperature. Additional experiments demonstrate that the assimilation of daily-updated satellite sea surface temperature data actually improves the model’s efficiency in representing typhoon-induced disturbances of sea surface temperature on a time scale of a few days. Assimilation of additional in situ data, such as temperature/salinity/ocean current information, further improves the model’s ability to represent the ocean currents near the coast accurately.

Comparison of ADCP measurements to Kuroshio current flow simulations for ocean current turbines

Numerical simulations are compared to field measurements of ocean currents in the Tokara Strait in this work. The Kuroshio Current flows through this strait and an ongoing ocean current turbine device is under development to harness its energy. The usefulness for engineering design input of modeling flow conditions for a large domain in space and time motivates the use of the numerical simulations. Confidence in the accuracy of the simulations can be provided from this comparison to Acoustic Doppler Current Profiler (ADCP) observations during the summer of 2018 at four locations. Numerical simulations of ocean currents which overlap in time with the field observations are presented in an attempt to compare the data on a time-domain basis. The simulations were produced using a Princeton Ocean Model based JCOPE-T-Tokara500 model. The analysis describes the capability of the numerical model to match the flow profile throughout the water column in the time domain and on a statistical basis using histograms and rose diagrams. While instances of speed peaks in the measurement data possibly representing internal waves did not readily appear in the simulation, overall the analysis supports the continued use of simulation current flow for project design input.

Improving the Estimation of Temperature and Salinity by Assimilation of Observed Sound Speed Profiles

AbstractSound Speed Profile (SSP) affecting underwater acoustics is closely related to the temperature and the salinity fields. It is of great value to obtain the temperature and the salinity information through the high-precision sound speed profiles. In this paper, a data assimilation scheme by introducing sound speed profiles as a new constraint is proposed within the framework of 3DVAR data assimilation (referenced as SSP-constraint 3DVAR (SSPC-3DVAR) ), which aims at improving the analysis accuracy of initial fields of the temperature and salinity in coastal sea areas. In order to validate the performance of the new assimilation scheme, ideal experiments are firstly carried out to show the advantages of the new proposed SSPC-3DVAR. Then the temperature, the salinity, and the SSP observations from field experiments in a coastal area are assimilated into the Princeton Ocean Model to validate the performance of short-time forecasts, adopting the SSPC-3DVAR scheme. Results show that it is efficient to improve the estimate accuracy by as much as 14.6% (11.1%) for the temperature (salinity), compared with the standard 3DVAR. It demonstrates that the proposed SSPC-3DVAR approach works better in practice than the standard 3DVAR and will primarily benefit from variously and widely distributed observations in the future.