将一个由水动力、粘性泥沙和水质生态模型搭建的综合模拟系统应用于长江口及邻近海域,采用枯季(11-4月)平均水文气象条件设置并驱动模型。将模拟结果与2003年3月观测资料及其他历史资料进行比对发现,模拟结果总体上很好地再现了研究海区枯季温度、盐度、悬浮泥沙、无机氮、无机磷等生态水质要素的水平和垂直分布趋势,说明模型可反映研究海区的水动力、泥沙和水质生态的关键过程。模拟结果中,DIN的分布趋势呈现出“近岸河口高、外海低”的特点,长江口和钱塘江口附近含量最高;DIP的高浓度区不仅仅在长江口及杭州湾,在浙江近海其它区域的浓度值也较高;浮游植物量最大的区域位于舟山以东海域。另外,3个条件响应数值试验表明,枯季长江口及邻近海区无机氮和无机磷均处于富余状态,尤其是无机氮,即使在陆源输入减半的情况下,仍然富余,没有被浮游植物充分利用;研究区域的高质量浓度悬浮泥沙通过限制光照强度控制着浮游植物的生长。
Abstract
An integrated model system that combines POM,a viscous sediment model and a water quality ecological model was applied to Changjiang River Estuary and its adjacent areas. The computational domain covered an area of 29°00'~33° 30' N, 120°00'~ 125°10' E with a mesh grid which has a finest horizontal resolution of 998 m and average of 2 789 m. In vertical, the grid was divided into 11 sigma layers. The model system was initialized and driven by averaged hydrometeorological conditions of the dry season from November to April. The runoff of the Changjiang Riveris given by the perennial average from December to April,with an average discharge of 15 200 m3/s and that for the Qiantang River was set as 925 m3/s. At the open boundary, eight tidal components were calculated, that is Q1,O,P1, K1,N2,M2,S2 and K2 and the amplitudes and phase lags were interpolated from the results of OTPS. The sediment mass concentration from the upper stream of the Changjiang and Qiantang Riversis 0. 547 kg/m3 and 0. 2 kg/m3,respectively. The biochemical elements for the upper streams of the two rivers were DIN 100 imol/L,DIP 0.7 imol/L, detrital organic matter of N 1. 0 imol/L,detrital organic matter of P 0. 08 imol/L, phytoplankton N 1. 2 imol/L, zooplankton N 0. 4 imol/L, and dissolved oxygen 8. 62 mg/L. On the off-sea open boundary, all the scalar fields use the approach given by Thatcher-Harleman boundary conditions at the hypothesis that the scalar field values outside the domain all remain the initial set values, and the time lag is set as 2 h.
According to the observed data in March 2003,high DIN concentration occurred near the Changjiang Estuary and Hangzhou bay and the value near the surface is 65 imol/L while that of the bottom layer is a little bit higher with the highest value of 80 imol/L near the bottom. For the modeling result, the highest DIN concentration zone matches with the observation data only the value is greater than 50 imol/L at the surface. According to the observed data in March 2003,the highest surface DIP concentration is approximately 0. 8 imol/L and that at the bottom approaches 1 imol/L. In the modeling results, the high concentration zones of DIP are not only in the Changjiang River Estuary and Hangzhou bay,but also in other areas off the Zhejiang Province and DIP concentration near the surface is 0. 75 imol/L, a little bit lower than that near the bottom, which matches well with the observed data. Chlorophyll a is usually used to represent the phytoplankton. The observed data of chlorophyll a in March 2003 indicates that there's little difference between the mass concentration near surface and bottom while the highest chlorophyll a mass concentration zone is located to the north of Changjiang River Estuary and comparatively high mass concentration zone at the open sea near the Zhoushan Islands. According to the observed data in November 2002,the zone of highest chlorophyll a mass concentration near surface is located near 30°~31°N,123°E and the zone of highest chlorophyll a mass concentration near the surface is distinctly different from that near the bottom. The modeling result of chlorophyll a doesn't match with the observed data, which may be caused by the transparency of the seawater. Because chlorophyll a was measured with fluorescence method, it's affected badly by the suspended sediment. Comparisons above show that the modeling results match well with the distributions of temperature, salinity, suspended sediment, inorganic nitrogen and inorganic phosphorus both horizontally and vertically, which indicates that the model system has simulated the key processes of hydrodynamics, siltation and water quality ecology in the computational domain.
Finally,three numerical experiments were run to check the modeling response to the change of initial conditions, which led to some helpful conclusions. For example, the high concentration zone of chlorophyll a is in the Changjiang River Estuary with a much higher value if the impact of suspended sediment on the light intensity is not considered. This indicates the great suppression on the growth of phytoplankton by sediment. If only the imput of DIN and detrital organic matter of nitrogen from the two rivers into the area was decreased by 50%,the modeling results of phytoplankton distribution in the surface layer is almost the same as that from standard operation test. But if only the imput of DIP and detrital organic matter of phosphorus from the two rivers into the area was decreased by 50% while keep the DIN value in the standard operation test unchanged, the resulted peak chlorophyll a mass concentration changed little. These three numerical experiments indicate that during dry season, the nutrients of inorganic nitrogen and phosphorous in the Changjiang River Estuary and its adjacent area are surplus especially the former. Even half of the input of the terrigenous inorganic nutrients is still surplus. The factor that controls the growth of phytoplankton is not nutrients but the high mass concentration silt that has great impacts on the light intensity and thus on the growth of phytoplankton.
关键词
长江口 /
三维模型 /
水质 /
生态
Key words
Changjiang River Estuary /
3D model /
ecosystem /
water quality
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] HU Yi-feng, GUO Peng-jun.Study on red tide and paralytic shellfish poisoning in Zhoushan area since 2001[J]. Marine Environmental Science ,2008 , 27(supp.1) :54-56.
胡益峰,郭朋军.2001年以来舟山海域赤潮及麻痹性贝毒发生情况的研究[J].海洋环境科学,2008 ,27(增刊1):54-56.
[2] ZHOU Ming-jiang, YAN Tian, ZOU Jing-zhong.Preliminary analysis of the characteristics of red tide areas in Changjiang River estuary and its adjacent sea[J]. Chinese Journal of Applied Ecology ,2003,14(7):1 031-1 038.
周明江,颜天,邹景忠.长江口邻近海域赤潮发生区基本特征初探[J].应用生态学报, 2003,14(7):1 031-1 038.
[3] LU Sai-ying, GE Ren-feng, LIU Li-hui.Seasonal variation and physical transportation of nutrients in continental areas of East China Sea[J]. Acta Oceanoligica Sinica,1996, 18(5):41-51.
陆赛英,葛人峰,刘丽慧.东海陆架水域营养盐的季节变化和物理输运的规律[J].海洋学报,1996, 18(5):41-51.
[4] PINCKNEY J L, PAERL H W, TESTER P, et al.The role of nutrient loading and eutrophication in estuarine ecology[J]. Environmental Health Perspectives, 2001 , 109(Suppl S): 699-706.
[5] CLOERNJ E.Turbidity as a control on phytoplankton biomass and productivity in estuaries[J]. Continental Shelf Research, 1987,7(11/12):1 367-1 381.
[6] CLOERN J E.The relative importance of light and nutrient limitation of phytoplankton growth: A simple index of coastal ecosystem sensitivity to nutrient enrichment[J]. Aquatic Ecology, 1999,33(1): 3-15.
[7] DONEY SC, GLOVER DM, NAJJAR R B.A new coupled, one-dimensional biological-physical model for the upper ocean: Applications to the JGOFS Bermuda Atlantic time-series study (BATS) site[J]. Deep-Sea Research II, 1996,43(2/3) :591-624.
[8] HOOD R R, BATES N R, CAPONE D G, et al.Modeling the effect of nitrogen fixation on carbon and nitrogen fluxes at BATS[J]. Deep-Sea Research II,2001 ,48:1609-1 648.
[9] MARRA J, HO C.Initiation of the spring bloom in the northeast Atlantic (47-degrees-N, 20 degrees W)-A numerical-simulation[J]. Deep-Sea Research II,1993 ,40( 1/2):55-73.
[10] LIMA I D,DONEY S C. A three-dimensional, multinutrient and size-structured ecosystem model for the North Atlantic[J]. Global Biogeochemical Cycles, 2004,18(3): GB3019, doi: 10.1029/2003G B002146.
[11] SKOGEN M, SVENDSEN D,BERNSTEN J, et al.Modeling the primary production in the North-Sea using a coupled three-dimensional physical-chemical-biological ocean model[J]. Estuarine Coastal and Shelf Science 1995 ,41(5): 545-565.
[12] MCCLAIN C, ARRIGO R,TAI K S, et al.Observations and simulations of physical and biological processes at ocean weather station P 1 951-1 980[J]. Journal of Geophysical Research O-ceans,1996, 101(C2):3697-3 713.
[13] OSCHLIES A, GARCON V.An eddy-permit ting coupled physical-biological model of the North Atlantic-1. Sensitivity to advection numerics and mixed layer physics[J]. Global Biogeo-chemical Cycles, 1999, 13(1):135-160.
[14] HOODR R, KOHLER K E, MCCREARY J P, et al.A four dimensional validation of a coupled physical-biological model of the Arabian Sea[J]. Deep-Sea Research II, 2003, 50( 22-26): 2917-2 945.
[15] MOOREJ K, DONEY SC, KLEYPASJ A, et al.An intermediate complexity marine ecosystem model for the global domain[J]. Deep Sea Research II, 2002 ,49(1-3) :403-462.
[16] CHEN C, JI R, SCHWAB D, et al.A coupled biological and physical model study of the ecosystem in Lake Michigan Part I: A 1-D experiment[J]. Ecological Modeling, 2002, 152(2/3): 145-168.
[17] ZHENG L,CHEN C,ZHANG F Y.Development of water quality model in the Satilla River Estuary, Georgia[J]. Ecological Modeling ,2004,178(3/4):457-482.
[18] GUAN Wei-bing, WONG Lai-ah, XU Dong feng. Modeling nitrogen and phosphorus cycles and dissolved oxygen in the Zhu-jiang Estuary part I. Model development[J]. Acta Oceanologica Sinica , 2003 ,25(1):52-60.
管卫兵,王丽娅,许东峰.珠江河口氮和磷循环及溶解氧的数值模拟——I.模式建立[J].海洋学报,2003,25(1):52-60.
[19] GUAN Wei-bing, WONG Lai-ah, XU Dong-feng.Modeling nitrogen and phosphorus cycles and dissolved oxygen in the Zhu-jiang Estuary part II. Model results[J]. Acta Oceanologica Sinica,2003,25(1 ):61-68.
管卫兵,王丽娅,许东峰.珠江河口氮和磷循环及溶解氧的数值模拟——I.模拟结果[J].海洋学报,2003,25(1):61-68.
[20] SHEN Yong-ming,ZHENG Yong-hong.3D water quality dynamic model on the trans portation and transformation of pollution in coastal area[J]. Progress in Natural Science, 2004, 14(6) : 694-699.
沈永明,郑永红.近岸海域污染物迁移转化的三维水质动力学模型[J].自然科学进展,2004, 14(6) :694-699.
[21] ZHANG Su-xiang.Modeling research on the 3D hydrodynamics and water quality in Yangtze Estuary[D]. Nanjing: Hohai University , 2007.
张素香. 长江口三维水动力与水质模拟研究[D].南京:河海大学,2007.
[22] ZHOU Ji-fu, WANG Tao, LI Jia-chun, et al.Influence of run-off and tide on sediment transport in the Yangtze Estuary[J]. Journal of Hydrodynamics: Ser. A,1999,14(1) :90-100.
周济福,王涛,李家春,等.径流与潮流对长江口泥沙输运的影响[J].水动力学研究与进展:A辑,1999, 14(1) :90-100.
[23] SHEN Huan-ting, HE Song-lin, PAN Ding-an, et al.A study of turbidity maximum in the Changjiang Estuary[J]. Acta Geo-graphica Sinica, 1992 ,47(5):472-479.
沈焕庭,贺松林,潘定安,等.长江河口最大浑浊带研究[J].地理学报,1992 ,47(5) :472-479.
[24] China Bay Records Compiling Committee. Chinese bay records: The fourteenth section[M]. Beijing: Ocean Press, 1992: 133-134.
中国海湾志编纂委员会.中国海湾志,第十四分册[M].北京:海洋出版社,1992:133-134.
[25] ZHU Jian-rong, SHEN Huan-ting.The expansion mechanism of Changjiang River diluted water[M]. Shanghai: Press of East China Normal University,1997:1-17.
朱建荣,沈焕庭.长江冲淡水扩展机制[M].上海:华东师范大学出版社,1997:1-17.
[26] BLUMBERG A F, MELLOR G R.A description of a three dimensional coastal ocean circulation model[C]// HEAPS N S. Three-dimensional coastal ocean models, coastal and estuarine Sci,Vol. 4. Washington DC: American Geophysical Union, 1987:1-19.
[27] WU Xiao yan, GUAN Wei-bing. 3D variable boundary numerical modeling of tidal current in the Xiangshan Bay[J]. Journal of Marine Sciences , 2009 ,27(2):23-31.
吴晓燕,管卫兵.象山港内三维动边界潮流的数值模拟[J].海洋学研究, 2009 ,27(2):23-31.
[28] GUAN Wei-bing, KOT SC, WOLANSKI E.3-D fluid-mud dynamics in the Jiaojiang Estuary, China[J]. Estuary, Coastal and Shelf Science, 2005, 65(4) :747-762.
[29] LI Xiao-hong, DONG Li-xian.The low-frequency sea level fluctuations in the Yellow Sea and the Bohai Sea during a winter storm[J]. Oceanologia et Limnologia Sinica, 2011,42(4): 467-473.
李晓红,董礼先.冬季寒潮期间黄渤海水位低频波动研究[J].海洋与湖沼,2011,42(4):467-473.
[30] EGBERT G D, EROFEEVA S Y.Efficient inverse modeling of barotropic ocean tides[J]. Journal of Atmospheric and Oceanic Technology , 2002 , 19(2) :183-204.
[31] FLATHER R A,HEAPS N S.Tidal computations for More-cambe Bay[J]. Geophys J R Astron Soc, 1975,42(2):489-517.
[32] The Ministry of Water Resource of the People's Republic of China. Bulletin of Chinese rivers and sediments 2009[M]. Beijing: China Water Power Press , 2009.
中华人民共和国水利部编.中国河流泥沙公报2009[M].北京: 中国水利水电出版社,2009.
[33] SHEN Huan-ting.Material flux of Yangtze Estuary[M]. Beijing: Ocean Press,2001.
沈煥庭. 长江河口物质通量[M].北京:海洋出版社,2001.
[34] THATCHER M L, HARLEMAN DR F.A mathematical mod-el for the prediction of unsteady salinity intrusion in estuaries[M]. Cambridge: Sea Grant Project Office, 1972.
[35] WONG L A,CHEN J C,XUE H, et al. A model study of the circulation in the Pearl River Estuary ( PRE) andits adjacent coastal waters: l. Simulations and comparison with observations [J]. Journal of Geophysical Research, 2003, 108: 3156,dio: 10.1029/2002JC001451.
[36] Ocean Atlas Commission.Ocean atlas of Bohai Sea,Yellow Sea and East China Sea ( hydrology)[G]. Bejing: Ocean Press, 1992:430-431.
海洋图集编委会.渤海黄海东海海洋图集(水文)[G].北京:海洋出版社,1992:430-431.
[37] PAN Chong, WANG Hui-qun, GUAN Wei-bing, et al.Study on real-time prediction of oil spill in Changjiang River Estuary and its adjacent sea areas[J]. Journal of Marine Sciences, 2011, 29(3) :176-186.
潘冲,王惠群,管卫兵,等.长江口及邻近海域溢油实时预测研究[J].海洋学研究,2011,29(3):176-186.
基金
国家重点基础研究发展规划资助项目(2010CB428704); 国家海洋局第二海洋研究所基本科研业务专项资金项目(JT1007,JG1106)