The impact of open boundary forcing on forecasting the East Australian Current using ensemble data assimilation

We investigate the performance of an eddy resolving regional ocean forecasting system of the East Australian Current (EAC) for both ensemble optimal interpolation (EnOI) and ensemble Kalman filter (EnKF) with a focus on open boundary model nesting solutions. The performance of nesting into a global re-analysis; nesting into the system’s own analysis; and nesting into a free model is quantified in terms of forecast innovation error. Nesting in the global reanalysis is found to yield the best results. This is closely followed by the system that nests inside its own analysis, which seems to represent a viable practical option in the absence of a suitable analysis to nest within. Nesting into a global reanalysis without data assimilation and nesting into an unconstrained model were both found to be unable to constrain the mesoscale circulation at all times. We also find that for a specific interior area of the domain where the EAC separation takes place, there is a mixture of results for all the systems investigated here and that, whilst the application of EnKF generates the best results overall, there are still times when not even this method is able to constrain the circulation in this region with the available observations.     

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.     

ARGO data assimilation into the ocean dynamics model with high spatial resolution using Ensemble Optimal Interpolation (EnOI)

The article proposes parallel implementation of the Ensemble Optimal Interpolation (EnOI) data assimilation (DA) method in eddy-resolving World Ocean circulation model. The results of DA experiments in North Atlantic with ARGO drifters are compared with the multivariate optimal interpolation (MVOI) DA scheme. The sensitivity of the model error, i.e., the difference between the model and observations depending on the number of ensemble elements, is also assessed and presented. The effectiveness of this method over the MVOI scheme is confirmed. The model outputs with and without assimilation are also compared with independent sea surface temperature data from ARMOR 3d.     

Asynchronous data assimilation with the EnKF

This study revisits the problem of assimilation of asynchronous observations, or four-dimensional data assimilation, with the ensemble Kalman filter (EnKF). We show that for a system with perfect model and linear dynamics the ensemble Kalman smoother (EnKS) provides a simple and efficient solution for the problem: one just needs to use the ensemble observations (that is, the forecast observations for each ensemble member) from the time of observation during the update, for each assimilated observation. This recipe can be used for assimilating both past and future data; in the context of assimilating generic asynchronous observations we refer to it as the asynchronous EnKF. The asynchronous EnKF is essentially equivalent to the four-dimensional variational data assimilation (4D-Var). It requires only one forward integration of the system to obtain and store the data necessary for the analysis, and therefore is feasible for large-scale applications. Unlike 4D-Var, the asynchronous EnKF requires no tangent linear or adjoint model.     

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.     

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.     

The impact of ocean data assimilation on seasonal predictions based on the National Climate Center climate system model

An ensemble optimal interpolation (EnOI) data assimilation method is applied in the BCC_CSM1.1 to investigate the impact of ocean data assimilations on seasonal forecasts in an idealized twin experiment framework. Pseudo-observations of sea surface temperature (SST), sea surface height (SSH), sea surface salinity (SSS), temperature and salinity (T/S) profiles were first generated in a free model run. Then, a series of sensitivity tests initialized with predefined bias were conducted for a one-year period; this involved a free run (CTR) and seven assimilation runs. These tests allowed us to check the analysis field accuracy against the “truth”. As expected, data assimilation improved all investigated quantities; the joint assimilation of all variables gave more improved results than assimilating them separately. One-year predictions initialized from the seven runs and CTR were then conducted and compared. The forecasts initialized from joint assimilation of surface data produced comparable SST root mean square errors to that from assimilation of T/S profiles, but the assimilation of T/S profiles is crucial to reduce subsurface deficiencies. The ocean surface currents in the tropics were better predicted when initial conditions produced by assimilating T/S profiles, while surface data assimilation became more important at higher latitudes, particularly near the western boundary currents. The predictions of ocean heat content and mixed layer depth are significantly improved initialized from the joint assimilation of all the variables. Finally, a central Pacific El Ni?o was well predicted from the joint assimilation of surface data, indicating the importance of joint assimilation of SST, SSH, and SSS for ENSO predictions.     

A comparison between 3DVAR and EnOI techniques for satellite altimetry data assimilation

Two conceptually different assimilation schemes, three dimensional variational (3DVAR) assimilation and Ensemble Optimum Interpolation (EnOI) are compared in the context of satellite altimetric data assimilation. Similarities and differences of the two schemes are briefly discussed and their impacts on the model simulation are investigated.With a tropical Pacific ocean model, two assimilation experiments of sea level anomaly (SLA) data from TOPEX/Poseidon are performed for 5 years from 1997 to 2001. Annual mean states of temperature and salinity fields are compared with analysis data and some independent observations. It is found that EnOI generally produces moderate improvements on both temperature and salinity fields, while changes induced by 3DVAR assimilation are strong and vary remarkably in different areas. For instance, 3DVAR tends to excessively modify the temperature field along the thermocline depth and even deteriorate the simulation, but it is more effective than EnOI below the thermocline depth. However, for the salinity field 3DVAR outperforms EnOI nearly for almost the whole layer. As the difference relative to the WOA01 analysis is compared, it is apparently reduced to below 0.3 psu in most areas in the 3DVAR experiment. On the other hand, the pattern of difference in the EnOI experiment resembles that of the simulation and the magnitude is only diminished to some extent. One advantage of EnOI is that it yields more consistent improvements even in areas where there are large model errors. It is more reliable than 3DVAR in such a sense. It is also revealed that the T–S relation plays a very important role in altimetric data assimilation. Further, the distinct performance of the two schemes can be partly accounted for by their inherent assumptions and settings.     

The representer method, the ensemble Kalman filter and the ensemble Kalman smoother: A comparison study using a nonlinear reduced gravity ocean model

This paper compares contending advanced data assimilation algorithms using the same dynamical model and measurements. Assimilation experiments use the ensemble Kalman filter (EnKF), the ensemble Kalman smoother (EnKS) and the representer method involving a nonlinear model and synthetic measurements of a mesoscale eddy. Twin model experiments provide the “truth” and assimilated state. The difference between truth and assimilation state is a mispositioning of an eddy in the initial state affected by a temporal shift. The systems are constructed to represent the dynamics, error covariances and data density as similarly as possible, though because of the differing assumptions in the system derivations subtle differences do occur. The results reflect some of these differences in the tangent linear assumption made in the representer adjoint and the temporal covariance of the EnKF, which does not correct initial condition errors. These differences are assessed through the accuracy of each method as a function of measurement density. Results indicate that these methods are comparably accurate for sufficiently dense measurement networks; and each is able to correct the position of a purposefully misplaced mesoscale eddy. As measurement density is decreased, the EnKS and the representer method retain accuracy longer than the EnKF. While the representer method is more accurate than the sequential methods within the time period covered by the observations (particularly during the first part of the assimilation time), the representer method is less accurate during later times and during the forecast time period for sparse networks as the tangent linear assumption becomes less accurate. Furthermore, the representer method proves to be significantly more costly (2–4 times) than the EnKS and EnKF even with only a few outer iterations of the iterated indirect representer method.     

Predicting the East Australian Current

Results are presented from an ensemble prediction study (EPS) of the East Australian Current (EAC) with a specific focus on the examination of the role of dynamical instabilities and flow dependent growing errors. The region where the EAC separates from the coast, is characterized by significant mesoscale eddy variability, meandering and is dominated by nonlinear dynamics thereby representing a severe challenge for operational forecasting. Using analyses from OceanMAPS, the Australian operational ocean forecast system, we explore the structures of flow dependent forecast errors over 7 days and examine the role of dynamical instabilities. Forecast ensemble perturbations are generated using the method of bred vectors allowing the identification of those perturbations to a given initial state that grow most rapidly. We consider a 6 month period spanning the Austral summer that corresponds to the season of maximum eddy variability. We find that the bred vector (BV) structures occur in areas of instability where forecast errors are large and in particular in regions associated with the Tasman Front and EAC extension. We also find that very few BVs are required to identify these regions of large forecast error and on that basis we expect that even a small BV ensemble would prove useful for adaptive sampling and targeted observations. The results presented also suggest that it may be beneficial to supplement the static background error covariances typically used in operational ocean data assimilation systems with flow dependent background errors calculated using a relatively cheap EPS.     

海洋资料同化对气候季节-年际预测技巧及初始场的影响试验

集合卡尔曼滤波(Ensemble Kalman filter, EnKF)是一种国内外广泛使用的海洋资料同化方案, 用集合成员的状态集合表征模式的背景误差协方差, 结合观测误差协方差, 计算卡尔曼增益矩阵, 有效地将观测信息添加到模式初始场中。由于季节、年际预测很大程度上受到初始场的影响, 因此资料同化可以提高模式的预测性能。本文在NUIST-CFS1.0预测系统逐日SST nudging的初始化方案上, 利用EnKF在每个月末将全场(full field)海表温度(sea surface temperature, SST)、温盐廓线(in-situ temperature and salinity profiles, T-S profiles)以及卫星观测海平面高度异常(sea level anomalies, SLA)观测资料同化到模式初始场中, 对比分析了无海洋资料同化以及加入同化后初始场的区别、加入海洋资料同化后模式提前1~24个月预测性能的差异以及对于厄尔尼诺-南方涛动(El Niño-southern oscillation, ENSO)预测技巧的影响。结果表明, 加入海洋资料同化能有效地改进初始场, 并且呈现随深度增加初始场改进越显著的特征。加入同化后, 对全球SST、次表层海水温度的平均预测技巧均有一定的提高, 也表现出随深度增加预测技巧改进越明显的特征。但加入海洋资料同化后, 模式对ENSO的预测技巧有所下降, 可能是由于模式误差的存在, 使得同化后的预测初始场从接近观测的状态又逐渐恢复到与模式动力相匹配的状态, 加剧了赤道太平洋冷舌偏西、中东部偏暖的气候平均态漂移。     

Impact of SARAL/AltiKa-Derived Sea Level Anomaly in a Data Assimilative Ocean Prediction System for the Indian Ocean

The impact of SARAL/AltiKa derived sea level anomaly (SLA) has been studied by assimilating it along with Jason-2 and Cryosat-2 SLA in the Princeton Ocean model (POM) using ensemble optimal interpolation (EnOI) technique. For isolating the extra benefit brought by SARAKL/Altika, a parallel run with assimilation of only Jason-2 and Cryosat-2 SLA has also been conducted. The importance of SARAL SLA in a data assimilative ocean prediction system has been evaluated with special emphasis on the improvement in thermocline depth, depth of the 20° isotherm, subsurface temperature and currents. Comparison with RAMA buoy has shown a positive impact of up to 13% for 20°C isotherm and up to 17% for thermocline depth after assimilating SARAL SLA. An overall improvement in temperature profile is also observed when compared with analogous profiles from RAMA buoys and Argo floats. Improvement in zonal currents away from the equator has also been noticed.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part III - Observation impact and observation sensitivity in the California Current System

The critical role played by observations during ocean data assimilation was explored when the Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation system was applied sequentially to the California Current circulation. The adjoint of the 4D-Var gain matrix was used to quantify the impact of individual observations and observation platforms on different aspects of the 4D-Var circulation estimates during both analysis and subsequent forecast cycles. In this study we focus on the alongshore and cross-shore transport of the California Current System associated with wind-induced coastal upwelling along the central California coast. The majority of the observations available during any given analysis cycle are from satellite platforms in the form of SST and SSH, and on average these data exert the largest controlling influence on the analysis increments and forecast skill of coastal transport. However, subsurface in situ observations from Argo floats, CTDs, XBTs and tagged marine mammals often have a considerable impact on analyses and forecasts of coastal transport, even though these observations represent a relatively small fraction of the available data at any particular time.During 4D-Var the observations are used to correct for uncertainties in the model control variables, namely the initial conditions, surface forcing, and open boundary conditions. It is found that correcting for uncertainties in both the initial conditions and surface forcing has the largest impact on the analysis increments in alongshore transport, while the cross-shore transport is controlled mainly by the surface forcing. The memory of the circulation associated with the control variable increments was also explored in relation to 7 day forecasts of the coastal circulation. Despite the importance of correcting for surface forcing uncertainties during analysis cycles, the coastal transport during forecast cycles initialized from the analyses has less memory of the surface forcing corrections, and is controlled primarily by the analysis initial conditions.Using the adjoint of the entire 4D-Var system we have also explored the sensitivity of the coastal transport to changes in the observations and the observation array. A single integration of the adjoint of 4D-Var can be used to predict the change that occurs when observations from different platforms are omitted from the 4D-Var analysis. Thus observing system experiments can be performed for each data assimilation cycle at a fraction of the computational cost that would be required to repeat the 4D-Var analyses when observations are withheld. This is the third part of a three part series describing the ROMS 4D-Var systems.     

海洋三维温盐流数值模拟研究的有关进展和问题

就海洋三维温盐流数值模拟使用的海洋模式和数据同化方法、在中尺度数值预报和再分析中的应用,以及所需支撑条件三方面,简述了国内外研究有关进展和问题。表述了开展大范围分辨中尺度乃至次中尺度涡的高分辨率海洋三维温盐流数值模拟正在研究的有关问题,扼要说明提供相匹配的高性能计算机模拟平台的必要性。初步探讨制约该研究快速发展的有关问题。     

Efficient identification of ocean thermodynamics in a physical/biogeochemical ocean model with an iterative Importance Sampling method

Efficient identification of parameters in numerical models remains a computationally demanding problem. Here we present an iterative Importance Sampling approach and demonstrate its application to estimating parameters that control the heat uptake efficiency of a physical/biogeochemical ocean model coupled to a simple atmosphere. The algorithm has similarities to a previously-developed ensemble Kalman filtering (EnKF) method applied to similar problems, but is more flexible and powerful in the case of nonlinear models and non-Gaussian uncertainties. The method is somewhat more computationally demanding than the EnKF but may be preferred in cases where the approximations that the EnKF relies upon are unsound. Our results suggest that the three-dimensional structure of ocean tracer fields may act as a useful constraint on ocean mixing and consequently the heat uptake of the climate system under anthropogenic forcing.     

EnKF集合同化下黄海海雾数值确定性预报初始场构造方法的探究

在黄海海雾的数值模拟中,EnKF（ensemble Kalman filter）是一种优于3DVAR（three-dimensional variational）的数据同化方法。研究发现,对EnKF初始场集合体采取常用的集合平均所产生的确定性预报初始场,会出现初始场中海雾在预报开始后就迅速消失以及接下来海雾难以生成的异常现象。通过详细的海雾个例研究,清晰地揭示并解释了此现象,指出这是集合平均造成初始场中云水与温度湿度之间存在不协调关系所导致的后果,并提出了一种择优加权平均方法来取代常用的集合平均。研究结果表明,海雾确定性预报采用择优加权平均所构建的初始场,可以消除上述异常现象,显著改进海雾模拟效果。     

A stochastic operational forecasting system of the Black Sea: Technique and validation

In this article, we present the latest version of an ensemble forecasting system of the hydrodynamics of the Black Sea, based on the GHER model. The system includes the Weakly Constrained Ensembles algorithm to generate random, but physically balanced perturbations to initialize members of the ensemble. On top of initial conditions, the ensemble accounts also for uncertainty on the atmospheric forcing fields, and on some scalar parameters such as river flows or model diffusion coefficients. The forecasting system also includes the Ocean Assimilation Kit, a sequential data assimilation package implementing the SEEK and Ensemble Kalman filters. A novel aspect of the forecasting system is that not only our best estimate of the future ocean state is provided, but also the associated error estimated from the ensemble of models. The primary goal of this paper is to quantitatively show that the ensemble variability is a good estimation of the model error, regardless of the magnitude of the forecast errors themselves. In order for this estimation to be meaningful, the model itself should also be well validated. Therefore, we describe the model validation against general circulation patterns. Some particular aspects critical for the Black Sea circulation are validated as well: the mixed layer depth and the shelfopen sea exchanges. The model forecasts are also compared with observed sea surface temperature, and errors are compared to those of another operational model as well.     

Corrections to ocean surface forcing in the California Current System using 4D variational data assimilation

The option for surface forcing correction, recently developed in the 4D-variational (4DVAR) data assimilation systems of the Regional Ocean Model System (ROMS), is presented. Assimilation of remotely-sensed (satellite sea surface height anomaly and sea surface temperature) and in situ (from mechanical and expendable bathythermographs, Argo floats and CTD profiles) oceanic observations has been applied in a realistic, high resolution configuration of the California Current System (CCS) to sequentially correct model initial conditions and surface forcing, using the Incremental Strong constraint version of ROMS-4DVAR (ROMS-IS4DVAR). Results from both twin and real data experiments are presented where it is demonstrated that ROMS-IS4DVAR always reduces the difference between the model and the observations that are assimilated. However, without corrections to the surface forcing, the assimilation of surface data can degrade the temperature structure at depth. When using surface forcing adjustment in ROMS-IS4DVAR the system does not degrade the temperature structure at depth, because differences between the model and surface observations can be reduced through corrections to surface forcing rather than to temperature at depth. However, corrections to surface forcing can generate abnormal spatial and temporal variability in the structure of the wind stress or surface heat flux fields if not properly constrained. This behavior can be partially controlled via the choice of decorrelation length scales that are assumed for the forcing errors. Abnormal forcing corrections may also arise due to the effects of model error which are not accounted for in IS4DVAR. In particular, data assimilation tends to weaken the alongshore wind stress in an attempt to reduce the rate of coastal upwelling, which seems to be too strong due to other sources of error. However, corrections to wind stress and surface heat flux improve systematically the ocean state analyses. Trends in the correction of surface heat fluxes indicate that, given the ocean model used and its potential limitations, the heat flux data from the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) used to impose surface conditions in the model are generally too low except in spring-summer, in the upwelling region, where they are too high. Comparisons with independent data provide confidence in the resulting forecast ocean circulation on timescales ～14 days, with less than 1.5 °C, 0.3 psu, and 9 cm RMS error in temperature, salinity and sea surface height anomaly, respectively, compared to observations.     

2011年3月日本福岛核电站核泄漏在海洋中的传输

使用全球版本的迈阿密等密度海洋环流模式对2011年3月日本福岛核电站泄漏在海洋中的传输以及扩散进行了数值模拟。数值模式中核废料(示踪物)排放情景采取等通量连续排放,排放时间从3月25日开始,分别持续20 d以及1 a,两种情形分别积分20 a。为了减少海洋环流年际变化带来的数值模拟的的不确定性,20 a的模式积分分别用2010年、1991-2011年、1971-1991年以及1951-1971年4个不同时段的NCEP/NCAR逐日再分析资料作为大气强迫场,因此每种排放情形包含4个数值试验。模拟结果的分析表明,不同核废料排放情景及其在不同时段大气资料对海洋模式的驱动下,模拟的示踪物总体的传输扩散路径(包括表层以及次表层)、传输速率以及垂直扩展的范围没有显著的差异。集合平均数值模拟的结果显示:在两种排放情景下,日本福岛核泄漏在海洋的传输路径受北太平洋副热带涡旋洋流系统主导,其传输路径首先主要向东,到达东太平洋后,再向南向西扩散至西太平洋,可能在10~15 a左右影响到我国东部沿海海域,且海洋次表层的传输信号比表层信号早5 a左右。通过进一步分析模式积分过程中最大示踪物浓度随时间变化发现,在积分第20 a(2031年3月),海洋表层和次表层浓度的最高值分别只有模式积分第一年浓度的0.1%和1%。在积分的20 a里,排放的核废料主要滞留在北太平洋海域(超过86%±1.5%的核废料在积分结束时,滞留在北太平洋),而在积分的前10 a(2021年之前),几乎所有的核废料滞留在北太平洋;在核废料的垂直分布上,主要集中在海洋表层至600 m的深度,在积分的20 a时间里,没有核废料信号扩散至1 000 m以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。     

2011年3月日本福岛核电站核泄漏在海洋中的传输

使用全球版本的迈阿密等密度海洋环流模式对2011年3月日本福岛核电站泄漏在海洋中的传输以及扩散进行了数值模拟。数值模式中核废料(示踪物)排放情景采取等通量连续排放,排放时间从3月25日开始,分别持续20 d以及1 a,两种情形分别积分20 a。为了减少海洋环流年际变化带来的数值模拟的的不确定性,20 a的模式积分分别用2010年、1991-2011年、1971-1991年以及1951-1971年4个不同时段的NCEP/NCAR逐日再分析资料作为大气强迫场,因此每种排放情形包含4个数值试验。模拟结果的分析表明,不同核废料排放情景及其在不同时段大气资料对海洋模式的驱动下,模拟的示踪物总体的传输扩散路径(包括表层以及次表层)、传输速率以及垂直扩展的范围没有显著的差异。集合平均数值模拟的结果显示:在两种排放情景下,日本福岛核泄漏在海洋的传输路径受北太平洋副热带涡旋洋流系统主导,其传输路径首先主要向东,到达东太平洋后,再向南向西扩散至西太平洋,可能在10~15 a左右影响到我国东部沿海海域,且海洋次表层的传输信号比表层信号早5 a左右。通过进一步分析模式积分过程中最大示踪物浓度随时间变化发现,在积分第20 a(2031年3月),海洋表层和次表层浓度的最高值分别只有模式积分第一年浓度的0.1%和1%。在积分的20 a里,排放的核废料主要滞留在北太平洋海域(超过86%±1.5%的核废料在积分结束时,滞留在北太平洋),而在积分的前10 a(2021年之前),几乎所有的核废料滞留在北太平洋;在核废料的垂直分布上,主要集中在海洋表层至600 m的深度,在积分的20 a时间里,没有核废料信号扩散至1 000 m以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。