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

何晏春; 郜永祺; 王会军; JOHANNESSENMOla; 于雷

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

1.
北京大学物理学院大气与海洋科学系, 北京 100871;中国科学院大气物理研究所竺可桢南森国际研究中心, 北京 100029
2.
中国科学院大气物理研究所竺可桢南森国际研究中心, 北京 100029;南森环境与遥感中心, 挪威卑尔根 5006;皮耶克尼斯气候研究中心, 挪威卑尔根 5007
3.
中国科学院大气物理研究所竺可桢南森国际研究中心, 北京 100029
4.
南森环境与遥感中心, 挪威卑尔根 5006

基金项目: 全球变化研究国家重大科学研究计划(2010CB951800);挪威研究理事会(Norwegian Research Council)"NorClim"项目。

计量
- 文章访问数: 1735
- HTML全文浏览量: 34
- PDF下载量: 1515
- 被引次数: 0
出版历程
- 收稿日期: 2011-08-31
- 修回日期: 2011-10-11

Transport of nuclear leakage from Fukushima Nuclear Power Plant in the North Pacific

1.
Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing, 100871;Nansen-Zhu International Research Center, Beijing, 100029
2.
Nansen-Zhu International Research Center, Beijing, 100029;Nansen Environmental and Remote Sensing Center, Bergen 5006, Norway;Bjerknes Center for Climate Research, Bergen 5007, Norway
3.
Nansen-Zhu International Research Center, Beijing, 100029
4.
Nansen Environmental and Remote Sensing Center, Bergen 5006, Norway

摘要

摘要: 使用全球版本的迈阿密等密度海洋环流模式对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以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。
- 核泄漏 /
- 海洋环流模式 /
- 示踪物 /
- 北太平洋副热带涡旋
Abstract: The massive nuclear leakage into ocean from Fukushima Nuclear Power Plant was observed on March 25^th,2011. The transport of leaked radioactive pollutant from the Fukushima Nuclear Power Plant was simulated using a global version of Miami Isopycnic Coordinate Ocean Model (MICOM). The simulation was performed by assuming constant and continuous leakage for 20 days (scenario 1) and for one year (scenario 2) starting from March 25th,2011 and was integrated for 20 years. Taking into account of the annual and inter-annul variability of the ocean circulation,the model for each scenario was forced with four different time periods of NCEP/NCAR daily atmospheric forcing,namely,repeated 2010,1991-2011,1971-1991 and 1951-1971. There is no remarkable difference of transport pathways (both on surface and subsurface layers),transport timescales and vertical penetration for the nuclear waste among different ensemble members. An ensemble mean for each scenario was obtained in order to reduce the uncertainty in the simulations. The results of the ensembles indicate that the nuclear pollutant for both scenarios transports eastward to eastern Pacific and thereafter follows a southwest pathway towards the western Pacific. It takes about 10 to 15 years to reach the coast of East Asia. The transport pathway is governed by the ocean current system of the North Pacific Subtropical Gyre. The transport at the subsurface (100-400 m depth) is about 5 years faster than that at the surface when the nuclear signal reaches the coast of the East Asia. The time evolution of maximum tracer concentration shows that the surface and subsurface maximum concentrations of leaked nuclear pollutant at the end of the simulation are only 0.1% and 1% of the values at the first year. Most of nuclear pollutant stays in the North Pacific throughout the integration with almost 100% of the released nuclear pollutants remaining in the North Pacific for first 10 years and 86±1.5% at the end of the integration. The penetration depth of the nuclear pollutant is less than 1 000 m with relative high concentration between the surface and 600 m depth. The simulations also suggest that the ocean circulation and mixing instead of the source function dominate the variability of the maximum concentration. It should be emphasized that a realistic source function is required and atmospheric fallout and role of ocean ecology should also be taken into account,in order to get a more reliable assessment of possible impact of the radioactive leakage on the ocean environment.
- nuclear leakage /
- Ocean General Circulation Model /
- ocean tracers /
- North Pacific Subtropical Gyre

HTML全文

参考文献(17)

唐森铭,商照荣. 近海辐射环境与生物多样性保护[J]. 核安全,2009(002):1-10.

NIES H,HARMS I H,KARCHER M J,et al. Anthropogenic radioactivity in the Arctic Oceanreview of the results from the joint German project[J]. The Science of the total Environment,1999,237-238:181-191,doi: 10.1016/S0048-9697(99)00134-5.

GERLAND S,LIND B,DOWDALL M,et al. 99Tc in seawater in the West Spitsbergen Current and adjacent areas[J]. Journal of Environmental Radioactivity,2003,69(1/2): 119-127,doi: 10.1016/S0265931X(03)000900.

GAO Y,DRANGE H,Bentsen M,et al. Simulating transport of nonChernobyl ¹³⁷Cs and ⁹⁰Sr in the North AtlanticArctic region[J]. Journal of environmental radioactivity,2004,71(1):1-16,doi: 10.1016/S0265-931X(03)00108-5.

GAO Y,DRANGE H,BENTSEN M,et al. Tracerderived transit time of the waters in the eastern Nordic Seas[J]. Tellus B,2005,57(4): 332 340,doi: 10.1111/j.1600-0889.2005.00153.x.

GAO Y,DRANGE H,JOHANNESSEN O M,et al. Sources and pathways of ⁹⁰Sr in the North AtlanticArctic region: present day and global warming[J]. Journal of Environmental Radioactivity,2009,100(5): 375-395,doi: 10.1016/j.jenvrad.2009.01.003.

JOHANNESSEN O M,VOLKOV V A,PETTERSSON L M,et al. Radioactivity and Pollution in the Nordic Seas and Arctic Region: Observations,Modeling,and Simulations . Springer,2010.

乔方利,王关锁,赵伟,等. 2011年3月日本福岛核泄漏物质输运扩散路径的情景模拟和预测[J].科学通报,2011,56(12): 964-971,doi: 10.1007/s11434-011-4513-0.

BLECK R,ROOTH C,HU D,et al. Salinitydriven thermocline transients in a windand thermohalineforced isopycnic coordinate model of the North Atlantic[J]. Journal of Physical Oceanography,1992,22(12):1486-1505,doi: 10.1175/1520-0485(1992)022<1486:SDTTIA>2.0.CO;2.

ENDOH T,MITSUDERA H,XIE SP,et al. Thermohaline structure in the subarctic North Pacific simulated in a general circulation model[J]. Journal of Physical Oceanography,Citeseer,2004,34(2): 360-371,doi: 10.1175/1520-0485(2004)034<0360:TSITSN>2.0.CO;2.

WANG D,WANG J,WU L,et al. Relative importance of wind and buoyancy forcing for interdecadal regime shifts in the Pacific Ocean. Science in China Series D: Earth Sciences,Springer,2003,46(5): 417-427,doi: 10.1360/03yd9037.

WANG D,WANG J,LIXIN W,et al. Regime shifts in the North Pacific simulated by a COADS-driven Isopycnal model[J]. Advances in Atmospheric Sciences,2003,20(5): 743-754,doi: 10.1007/BF02915399.

DUTAY J,BULLISTER J L,DONEY S C,et al. Evaluation of ocean model ventilation with CFC-11: comparison of 13 global ocean models[J]. Ocean Modelling,2002,4(2): 89-120,doi: 10.1016/S1463-5003(01)00013-0.

ORRE S,GAO Y,DRANGE H,et al. Diagnosing ocean tracer transport from Sellafield and Dounreay by equivalent diffusion and age[J]. Advances in Atmospheric Sciences,2008,25(5): 805-814,doi: 10.1007/s00376-008-0805-y

KALNAY E,KANAMITSU M,KISTLER R,et al. The NCEP/NCAR 40year reanalysis project[J]. Bulletin of the American Meteorological Society,1996,77(3):437,doi: 10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

MONTEREY G,LEVITUS S,1997.Seasonal Variability of Mixed Layer Depth for the World Ocean . NOAA Atlas NESDIS 14,U.S. Gov. Printing Office,Wash D C.

LIU Q,HU H. A subsurface pathway for low potential vorticity transport from the central North Pacific toward Taiwan Island[J]. Geophysical Research Letters,2007,34(12):L12710,doi: 10.1029/2007GL029510.

施引文献

资源附件(0)

访问统计