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.     

基于ROMS模式利用集合最优插值同化沿轨海表面高度异常数据的研究

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.     

利用热带太平洋TAO浮标阵列评估全球高分辨率浪-潮-流耦合模式同化结果

In order to evaluate the assimilation results from a global high resolution ocean model, the buoy observations from tropical atmosphere ocean(TAO) during August 2014 to July 2015 are employed. The horizontal resolution of wave-tide-circulation coupled ocean model developed by The First Institute of Oceanography(FIOCOM model) is 0.1°×0.1°, and ensemble adjustment Kalman filter is used to assimilate the sea surface temperature(SST), sea level anomaly(SLA) and Argo temperature/salinity profiles. The simulation results with and without data assimilation are examined. First, the overall statistic errors of model results are analyzed. The scatter diagrams of model simulations versus observations and corresponding error probability density distribution show that the errors of all the observed variables, including the temperature, isotherm depth of 20°C(D20), salinity and two horizontal component of velocity are reduced to some extent with a maximum improvement of 54% after assimilation. Second, time-averaged variables are used to investigate the horizontal and vertical structures of the model results. Owing to the data assimilation, the biases of the time-averaged distribution are reduced more than70% for the temperature and D20 especially in the eastern Pacific. The obvious improvement of D20 which represents the upper mixed layer depth indicates that the structure of the temperature after the data assimilation becomes more close to the reality and the vertical structure of the upper ocean becomes more reasonable. At last,the physical processes of time series are compared with observations. The time evolution processes of all variables after the data assimilation are more consistent with the observations. The temperature bias and RMSE of D20 are reduced by 76% and 56% respectively with the data assimilation. More events during this period are also reproduced after the data assimilation. Under the condition of strong 2014/2016 El Ni?o, the Equatorial Undercurrent(EUC) from the TAO is gradually increased during August to November in 2014, and followed by a decreasing process. Since the improvement of the structure in the upper ocean, these events of the EUC can be clearly found in the assimilation results. In conclusion, the data assimilation in this global high resolution model has successfully reduced the model biases and improved the structures of the upper ocean, and the physical processes in reality can be well produced.     

Generation of high resolution sea surface temperature using multi-satellite data for operational oceanography

In the present article, we introduce a high resolution sea surface temperature(SST) product generated daily by Korea Institute of Ocean Science and Technology(KIOST). The SST product is comprised of four sets of data including eight-hour and daily average SST data of 1 km resolution, and is based on the four infrared(IR) satellite SST data acquired by advanced very high resolution radiometer(AVHRR), Moderate Resolution Imaging Spectroradiometer(MODIS), Multifunctional Transport Satellites-2(MTSAT-2) Imager and Meteorological Imager(MI), two microwave radiometer SSTs acquired by Advanced Microwave Scanning Radiometer 2(AMSR2), and Wind SAT with in-situ temperature data. These input satellite and in-situ SST data are merged by using the optimal interpolation(OI) algorithm. The root-mean-square-errors(RMSEs) of satellite and in-situ data are used as a weighting value in the OI algorithm. As a pilot product, four SST data sets were generated daily from January to December 2013. In the comparison between the SSTs measured by moored buoys and the daily mean KIOST SSTs, the estimated RMSE was 0.71°C and the bias value was –0.08°C. The largest RMSE and bias were 0.86 and –0.26°C respectively, observed at a buoy site in the boundary region of warm and cold waters with increased physical variability in the Sea of Japan/East Sea. Other site near the coasts shows a lower RMSE value of 0.60°C than those at the open waters. To investigate the spatial distributions of SST, the Group for High Resolution Sea Surface Temperature(GHRSST) product was used in the comparison of temperature gradients, and it was shown that the KIOST SST product represents well the water mass structures around the Korean Peninsula. The KIOST SST product generated from both satellite and buoy data is expected to make substantial contribution to the Korea Operational Oceanographic System(KOOS) as an input parameter for data assimilation.     

The impact of mean dynamic topography on a sea-level anomaly assimilation in the South China Sea based on an eddy-resolving model

The sea-level anomaly (SLA) from a satellite altimeter has a high accuracy and can be used to improve ocean state estimation by assimilation techniques. However, the lack of an accurate mean dynamic topography (MDT) is still a bothersome issue in an ocean data assimilation. The previous studies showed that the errors in MDT have significant impacts on assimilation results, especially on the time-mean components of ocean states and on the time variant parts of states via nonlinear ocean dynamics. The temporal-spatial differences of three MDTs and their impacts on the SLA analysis are focused on in the South China Sea (SCS). The theoretical analysis shows that even for linear models, the errors in MDT have impacts on the SLA analysis using a sequential data assimilation scheme. Assimilation experiments, based on EnOI scheme and HYCOM, with three MDTs from July 2003 to June 2004 also show that the SLA assimilation is very sensitive to the choice of different MDTs in the SCS with obvious differences between the experimental results and observations in the centre of the SCS and in the vicinity of the Philippine Islands. A new MDT for assimilation of SLA data in the SCS was proposed. The results from the assimilation experiment with this new MDT show a marked reduction (increase) in the RMSEs (correlation coefficient) between the experimental and observed SLA. Furthermore, the subsurface temperature field is also improved with this new MDT in the SCS.     

Quality control methods for KOOS operational sea surface temperature products

Sea surface temperature SST obtained from the initial version of the Korea Operational Oceanographic System(KOOS) SST satellite have low accuracy during summer and daytime. This is attributed to the diurnal warming effect. Error estimation of SST data must be carried out to use the real-time forecasting numerical model of the KOOS. This study suggests two quality control methods for the KOOS SST system. To minimize the diurnal warming effect, SSTs of areas where wind speed is higher than 5 m/s were used. Depending on the wind threshold value, KOOS SST data for August 2014 were reduced by 0.15°C. Errors in SST data are considered to be a combination of random, sampling, and bias errors. To estimate bias error, the standard deviation of bias between KOOS SSTs and climatology SSTs were used. KOOS SST data yielded an analysis error standard deviation value similar to OSTIA and NOAA NCDC(OISST) data. The KOOS SST shows lower random and sampling errors with increasing number of observations using six satellite datasets. In further studies, the proposed quality control methods for the KOOS SST system will be applied through more long-term case studies and comparisons with other SST systems.     

Ocean temperature responses to Typhoon Mstsa in the East China Sea

The MASNUM wave-tide-circulation coupled model, with 21 layers in the vertical and （1/8） °horizontal resolution, was employed to investigate the oceanic responses to Typhoon Mstsa which traversed the East China Sea （ECS） during the period of 4 - 6 August, 2005. Numerical experiment results are analyzed and compared with observation. The responses of the sea surface temperature （SST）, in a focused area of （27° -29°N, 121° - 124°E）, include heating and cooling stages. The heating is mainly due to warm Kuroshio water transportation and downwelling due to the water accumulation. In the cooling stage, the amplitude of the simulated cold wake （ -3℃ ）, located on the right side of this typhoon track, is compared quite well with that of the satellite observed SST data. The wave-induced mixing（Bv） plays a key role for the SST cooling. Bv still plays a leading role, which accounts for 36%, for the ocean temperature drop in the upper ocean of 0 - 40 m, while the upwelling is responsible for 84% of the cooling for the lower layer of 40 - 70 m. The mixed layer depth （MLD） increased quickly from 28 to 50 m in the typhoon period. However, the simulated MLD without the wave-induced vertical mixing, evolution from 13 to 32 m, was seriously underestimated. The surface wave is too important to be ignored for the ocean responses to a typhoon.     

Evaluating the contribution of satellite measurements to the re-construction of three-dimensional ocean temperature fields in combination with Argo profiles

Assimilation systems absorb both satellite measurements and Argo observations. This assimilation is essential to diagnose and evaluate the contribution from each type of data to the reconstructed analysis, allowing for better configuration of assimilation parameters. To achieve this, two comparative reconstruction schemes were designed under the optimal interpolation framework. Using a static scheme, an in situ-only field of ocean temperature was derived by correcting climatology with only Argo ...     

东海对台风云娜的海洋响应

Many typhoons pass through the East China Sea(ECS) and the oceanic responses to typhoons on the ECS shelf are very energetic. However, these responses are not well studied because of the complicated background oceanic environment. The sea surface temperature(SST) response to a severe Typhoon Rananim in August 2004 on the ECS shelf was observed by the merged cloud-penetrating microwave and infrared SST data. The observed SST response shows an extensive SST cooling with a maximum cooling of 3°C on the ECS shelf and the SST cooling lags the typhoon by about one day. A numerical model is designed to simulate the oceanic responses to Rananim.The numerical model reasonably simulates the observed SST response and thereby provides a more comprehensive investigation on the oceanic temperature and current responses. The simulation shows that Rananim deepens the ocean mix layer by more than 10 m on the ECS shelf and causes a cooling in the whole mixed layer. Both upwelling and entrainment are responsible for the cooling. Rananim significantly deforms the background Taiwan Warm Current on the ECS shelf and generates strong Ekman current at the surface. After the typhoon disappears, the surface current rotates clockwise and vertically, the current is featured by near inertial oscillation with upward propagating phase.     

沿海观测站海表温度质量控制和半日特征分析

Sea surface temperature(SST) data obtained from coastal stations in Jiangsu, China during 2010–2014 are quality controlled before analysis of their characteristic semidiurnal and seasonal cycles, including the correlation with the variation of the tide. Quality control of data includes the validation of extreme values and checking of hourly values based on temporally adjacent data points, with 0.15°C/h considered a suitable threshold for detecting abnormal values. The diurnal variation amplitude of the SST data is greater in spring and summer than in autumn and winter. The diurnal variation of SST has bimodal structure on most days, i.e., SST has a significant semidiurnal cycle. Moreover, the semidiurnal cycle of SST is negatively correlated with the tidal data from March to August, but positively correlated with the tidal data from October to January. Little correlation is detected in the remaining months because of the weak coastal–offshore SST gradients. The quality control and understanding of coastal SST data are particularly relevant with regard to the validation of indirect measurements such as satellitederived data.     

Satellite-derived SST validation based on in-situ data during summer in the East China Sea and western North Pacific

Satellite-derived sea surface temperature (SST) is validated based on in-situ data from the East China Sea (ECS) and western North Pacific where most typhoons, which make landfall on the Korean peninsula, are formed and pass. While forecasting typhoons in terms of intensity and track, coupled ocean-typhoon models are significantly influenced by initial ocean condition. Potentially, satellite-derived SST is a very useful dataset to obtain initial ocean field because of its wide spatial coverage and high temporal resolution. In this study, satellite-derived SST from various sources such as Tropical Rainfall Measuring Mission Microwave Imager (TMI), Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) and New Generation Sea Surface Temperature for Open Ocean (NGSST-O) datasets from merged SSTs were compared with in-situ observation data using an indirect method which is using near surface temperature for validation of satellite derived SST. In-situ observation data included shipboard measurements such as Expendable Bathythermograph (XBT), and Conductivity, Temperature, Depth (CTD), and Argo buoy data. This study shows that in-situ data can be used for microwave derived SST validation because homogeneous features of seawater prevail at water depths of 2 m to 10 m under favorable wind conditions during the summer season in the East China Sea. As a result of validation, root-mean-square errors (RMSEs) are shown to be 0.55 °C between microwave SST and XBT/CTD data mostly under weak wind conditions, and 0.7 °C between XBT/CTD measurement and NGSST-O data. Microwave SST RMSE of 0.55 °C is a potentially valuable data source for general application. Change of SST before and after typhoon passing may imply strength of ocean mixing due to upwelling and turbulent mixing driven by the typhoon. Based on SST change, ocean mixing, driven by Typhoon Nari, was examined. Satellite-derived SST reveals a significant SST drop around the track immediately following the passing of Typhoon Nari in October, 2007.     

4DVAR Data Assimilation with the Regional Ocean Modeling System (ROMS): Impact on the Water Mass Distributions in the Yellow Sea

In this study, we evaluate the performance of the recently developed incremental strong constraint 4-dimensional variational (4DVAR) data assimilation applied to the Yellow Sea (YS) using the Regional Ocean Modeling System (ROMS). Two assimilation experiments are compared: assimilating remote-sensed sea surface temperature (SST) and both the SST and in-situ profiles measured by shipboard CTD casts into a regional ocean modeling from January to December of 2011. By comparing the two assimilation experiments against a free-run without data assimilation, we investigate how the assimilation affects the hydrographic structures in the YS. Results indicate that the SST assimilation notably improves the model behavior at the surface when compared to the non-assimilative free-run. The SST assimilation also has an impact on the subsurface water structure in the eastern YS; however, the improvement is seasonally dependent, that is, the correction becomes more effective in winter than in summer. This is due to a strong stratification in summer that prevents the assimilation of SST from affecting the subsurface temperature. A significant improvement to the subsurface temperature is made when the in-situ profiles of temperature and salinity are assimilated, forming a tongue-shaped YS bottom cold water from the YS toward the southwestern seas of Jeju Island.     

Assimilation of sea surface temperature in the Northwest Pacific Ocean and its marginal seas using the ensemble Kalman filter

Satellite-borne sea surface temperature (SST) data were assimilated with the ensemble Kalman filter (EnKF) in a Northwest Pacific Ocean circulation model to examine the effect of data assimilation. The model domain included the northwestern part of the Pacific Ocean and its marginal seas, such as the Yellow Sea and East/Japan Sea. The performance of the data assimilation was evaluated by comparing the simulated ocean state with that observed. Spatially averaged root-mean-squared errors in the SST and sea surface height (SSH) decreased by 0.44 °C and 4 cm, respectively, by the assimilation. The results of the numerical experiments substantiated the effectiveness of the SST assimilation via the EnKF for all marginal seas, as well as the Kuroshio region. The benefit of the data assimilation depended on the characteristics of each marginal sea. The variation of the SST in the East/Japan Sea and the Kuroshio extension (KE) region were improved 34% and those in the Yellow Sea 12.5%. The variation of the SSH was improved approximately 36% in the KE region. This large improvement was achieved in the deep-water regions because assimilation of SST data corrected the separation point of the western boundary currents, such as the Kuroshio and the East Korea Warm Current, and the associated horizontal surface currents. The SST assimilation via the EnKF also improved the subsurface temperature profiles. The effectiveness of SST assimilation was seasonally dependent, with the improvement being relatively larger in winter than in summer, which was related to the seasonal variation of the vertical mixing and stratification in the ocean surface layer.     

Research and development of the New Generation Sea Surface Temperature for Open Ccean (NGSST-O) product and its demonstration operation

Real-time generation and distribution of the New Generation Sea Surface Temperature for Open Ocean (NGSST-O) product began in September 2003 as a demonstration operation of the Global Ocean Data Assimilation Experiment (GODAE) High-Resolution Sea Surface Temperature Pilot Project. Satellite sea surface temperature (SST) observations from infrared radiometers (AVHRR, MODIS) and a microwave radiometer (AMSR-E) are objectively merged to generate the NGSST-O product, which is a quality-controlled, cloud-free, high-spatial-resolution (0.05° gridded), wide-coverage (13–63° N, 116–166° E), daily SST digital map. The NGSST-O demonstration operation system has been developed in cooperation with the Japanese Space Agency (JAXA) and has produced six years of continuous data without gaps. Comparison to in situ SSTs measured by drifting buoys indicates that the root mean-square error of NGSST-O has been kept at approximately 0.9°C.     

Spatial classification of La Niña events

Thirty La Niña events have been selected from monthly mean sea surface temperature (SST) data of the Hadley Centre Global Sea Ice and Sea Surface Temperature (HadISST) and Centennial Observation-Based Estimates (COBE SST2) datasets from 1870 to 2013 based on a criterion of –0.5°С for a cold anomaly in the Niño 3.4 region (5° N–5° S, 170°–120° W) and its minimum duration of 5 months. The selected events are classified by hierarchical clustering analysis according to two characteristics: geographic coordinates and SST anomalies during the mature phase of La Niña. The objective classification method identifies two types of La Niña differing by the evolution of negative SST anomalies in the equatorial Pacific and by the Southern Oscillation Index.     

中国近海现场海洋观测系统设计评估

中国科学院正在发展一个在中国近海(包括黄海、东海和南海)现场海洋观测系统。观测系统包括3个沿岸观测站点、4个近海离岸浮标和由观测船只按固定航线做的船舶观测断面。观测站点、浮标和断面的位置已经预先确定,这个计划在2008-2011实施。利用基于卡尔曼理论的样本集合方法对这样一个能够监测大尺度的季节和年季变率的观测系统设计进行了评估。根据卡尔曼滤波理论,用集合样本的方法能够给出经过同化这个观测系统位置的观测资料后能够减少多少分析误差和分析场的不确定性。用2个来自不同模式、不同分辨率的模式的结果作为集合样本来计算静态的背景误差协方差,这2套样本分别是来自分辨率是0.5°×0.5°的模式同化结果和高分辨0.125°×0.125°的模式结果。由这2个不同资料得到的结果是一致的。发现来自3个近岸和4个离岸浮标得到的观测能够有效地减少SST在渤海、黄海、东海和南海中部的分析误差。然而在越南东部和台湾东部海域,分析误差减少的百分比相对要小。最后,给出了中国近海最优的观测位置序列设计。     

基于卫星遥感海温数据的南海SST预报误差订正

尝试利用卫星遥感高分辨率海表温度资料GHRSST (Group for High Resolution Sea Surface Temperature) 与海表温度(sea surface temperature, SST)数值预报产品之间的误差, 建立一种南海SST模式预报订正方法。首先, 利用南海的Argo浮标上层海温数据对GHRSST 海温数据进行验证, 结果表明两者之间均方根误差约为0.3℃, 相关系数为0.98, GHRSST 海温数据可用于南海业务化数值预报SST的订正。预报订正后的SST与Argo浮标海温数据相比, 24h、48h和72h的均方根误差均由0.8℃左右下降到0.5℃以内。与GHRSST 海温数据相比, 南海北部海域(110°E—121°E, 13°N—23°N)订正后的24h、48h和72h的SST预报空间误差均显著减小, 在冷空气影响南海期间或中尺度涡存在的过程中, SST预报订正效果也较为显著。因此, 该方法可考虑在南海业务化SST数值预报系统中应用。