徐福敏,李 婷
(河海大学港口海岸与近海工程学院,江苏 南京 210098)
海洋升降流研究对于气候变化[1-2]、海洋环境、渔业养殖和社会经济均具有重大意义[3]。上升流可以把海底丰富的营养物质带到海洋上层,是渔场形成的主要条件[4-5]。海岸河口水域上升流的产生、发展与风有密切关系。笔者根据升降流理论,利用中分辨率成像光谱仪(MODIS)数据反演加拿大西北部高纬度Mackenzie河口水域海洋表面温度;结合当地实测风速、风向数据综合分析上升流产生条件、发展、强度,以及它们随纬度、风向及风速的变化关系,并分析上升流异常现象产生的原因,以期为当地社区生态、渔业和海洋环境提供科学依据。
MODIS是美国国家航空航天局(NASA)新一代对地观测系统(EOS)中搭载在TERRA和AQUA卫星上的重要传感器[6-7],它的多波段数据可以同时提供反映地表面状况、云边界、云特性、海洋水色、地表温度、云顶温度和大气温度等信息。笔者基于MODIS数据(第31,32波段),采用多通道分裂窗算法反演海洋表面温度。原始MODIS数据(来源于NASA的MODIS网站)及常用的算法模式见文献[8-10]。
海岸河口水域上升流基本运动方程为[11]
对于与海岸成β角的直线风系生成的上升流,边界条件为
式中:u,v,w——垂直于岸线方向、平行于岸线方向、沿水深方向的水质点运动速度;ζ——自由水位;Az,Al——垂直和水平湍流黏滞系数;f——科氏力;ρ——海水密度;τx——垂直于海岸的风应力;τy——平行于海岸的风应力;L——直线风系宽度。
设xOz平面流函数ψ(x,z)为文献[11]中流函数的实部:
式中:ω——地转角速度;φx,φy——与离岸距离、水深有关的量在x和y方向上的投影[11];φ——地理纬度。
平面流速度为
由式(3)可知,流函数由两部分组成,一部分由τx决定,另一部分由τy决定。对于平行于海岸的直线风系(τx=0),其流函数为式(3)的特殊情况。因此,非平行于海岸的直线风系与平行于海岸的直线风系导致的近岸流动相似。
对于与岸线交角为β的风系,由式(3)可得其上升流垂直质量运输为
当φx(x,z),φy(x,z)取最大负值时M有最大值,此时将τx=τsinβ,τy=τcosβ代入式(5)得
式中τ为风应力。
根据式(6)绘出M随β的变化曲线,如图1所示。由图1可知,β=21.5°时上升流最强烈;较强上升流产生的角度区域为0°~40°;由于式(6)分母有f项,因此纬度越低,上升流越强。
Mackenzie河口水域位于加拿大西北部,北纬52°13'~69°09'、西经108°11'~134°05'之间。Mackenzie河的众多南部支流(主要有19条)向北流入Mackenzie河口,继而向北进入Beaufort海。Mackenzie河口靠近北极,地理位置独特,冰期较长,全年大部分时间被海冰覆盖。每年6月初海冰开始融化,随温度逐渐升高,海冰融化范围变大,开冰季节一般为7—11月。夏秋季Mackenzie河将南部高温淡水(15℃以上)带入低温的Mackenzie河口和Beaufort海,一薄层低密度、高温淡水像羽流一样浮于低温海水之上,同时向外海扩散。此时Mackenzie河口水域海洋表面温度普遍较高(一般高于10℃),而低密度、高温淡水之下却依然是低温海水。若该水域发生上升流,上升流会使下层低温海水上翻,上升流产生水域的海洋表面温度便会低于周围海洋表面温度。
图1 β对上升流的影响Fig.1 Impact of angle between wind direction and coastline(β)on upwelling
本文采用数据为National Climate Data and Information Archive提供的Pelly气象站(-135.43°W,69.62°N)的观测风资料,Pelly气象站地理位置见图2(a)。由风速、风向观测数据可知,Mackenzie河口主要风向为东南风与西北风,其中东南风最常见(30%),常见风速为4~6 m/s,月峰值达14 m/s;西北风发生频率相对较小(20%~25%),但风速较大,常见风速为6~8 m/s,月峰值达26 m/s[12-13]。本文根据2007年8—9月Mackenzie河口海洋表面温度的MODIS数据和此期间Pelly气象站观测风资料研究该河口水域上升流特征,重点分析上升流的产生条件、发展、强弱,以及它们随纬度、风向及风速的变化关系。
根据预处理后的MODIS数据反演得到Mackenzie河口2007年8—9月海洋表面温度(TSST)分布。图2为典型日(8月22日、23日、29日和31日)Mackenzie河口海洋表面温度分布,图3为8月18日至9月1日Pelly气象站观测的风速、风向(风向采用笛卡尔坐标系,0°为西风,逆时针旋转)。
图2 Mackenzie河口海洋表面温度分布Fig.2 SST distribution of Mackenzie Estuary
一般来说,低纬度海域上升流现象较强烈,高纬度海域上升流现象较弱。然而图2表明位于高纬度的Mackenzie河口夏季存在明显、强烈的上升流现象。
图3 2007年8月18日至9月1日(第230~244天)Pelly气象站风速、风向Fig.3 Wind speed and direction at Pelly Station from August 18 to September 1,2007
图2(a)(b)分别为第234~235天Mackenzie河口海洋表面温度分布,可以看出在近岸处有上升流产生,下层低温海水上翻造成局部海洋表面温度降低,上升流产生范围在近海水域至20 m等深线附近。图3表明在上升流产生的前2 d(第232~234天)主要风向为东向风(160°~200°),风速约为7 m/s,风吹刮时间长于2 d。风向在第235天时由东南转为偏南(100°左右)。由此可见,东向风长时间吹刮会产生上升流,当风向偏离东风时上升流现象将减弱。
图2(c)为第241天海洋表面温度分布,此时冷水上翻现象相当明显。图3显示,在此上升流现象产生之前东风为主的风吹刮时间长于4 d,风速较大(8 m/s左右),风向接近于正东风(180°),此时上升流范围超过20 m等深线。风向在第243天时已转向西北(-20°~-60°),风向变化较大,风速较小(2 m/s左右),图2(d)显示此时上升流现象已消失。
由此可见,Mackenzie河口水域在东向风长时间吹刮时会产生上升流,正东风时最强,其他风向不能产生上升流。由于Mackenzie河口岸线自然角度为200°,可知Mackenzie河口160°~200°的风向适于产生较强上升流,二者相一致。同时风速大小也影响上升流强度和范围,东风持续风速越大,上升流强度越强、影响范围越大。
Mackenzie河口海洋动力条件复杂,河口径流、近海海冰、波浪、水流、内波等对上升流的存在均有影响。本文研究中第217天(2007-08-05)Mackenzie河口海洋表面温度分布(图4)显示无上升流产生,但其对应的第213~218天Pelly气象站风速、风向表明,在此之前东风已经持续长时间(图5),这显然与理论及第3节的资料分析结果相悖。
图4 2007年8月5日Mackenzie河口海洋表面温度分布Fig.4 SST distribution of Mackenzie Estuary in August 5,2007
图5 2007年8月1—6日Pelly气象站风速、风向Fig.5 Wind speed and direction at Pelly Station from August 1 to 6,2007
海冰数据显示2007年8月5日Mackenzie河口与外海相接处仍有大面积海冰未完全融化,东风吹过的海面仍有部分海冰存在,因此表层温度较高的低盐水难以产生离岸流动,从而阻碍底层低温海水的上升流动,不能形成上升流。虽然此时南部Mackenzie河的径流量仍然较大,高温、低密度淡水覆盖在海水上,但是近岸海冰的封锁使Mackenzie河口开敞水域面积受限制,因此长时间的持续东向风无法形成上升流,可见上升流的形成和当地海域综合物理海洋环境紧密相关。
基于升降流理论,利用MODIS数据反演海洋表面温度,结合相应的实测风速、风向探讨了Mackenzie河口夏秋季(2007年8—9月)上升流与风场的关系以及上升流产生的条件、强弱、范围等。结果表明:
a.上升流与纬度、风向、风速等因素有关。Mackenzie河口虽处于高纬度区域,但持续东向风使该水域产生上升流。当东向风吹刮时间长于半天,即有上升流产生。一般情况下,正东风(180°左右)所致上升流最强,东向偏南(160°~180°)与东向偏北(180°~200°)风也能产生显著上升流。吹刮时间不足或其他风向均不能产生上升流。
b.持续东风吹刮条件下,风速越大产生的上升流越强,底层冷水上翻现象越明显。上升流产生的范围一般为沿岸到20 m等深线附近,较强东向风引起的上升流会越过20 m等深线处到达30 m等深线附近。
c.由于Mackenzie河口水域复杂的海洋环境条件,很多海洋动力因素影响到上升流的产生和强度。例如当近海的冰未完全融化导致河口开敞水域面积受限制,此时即使有长时间持续的东向风也无法产生上升流。总之,需要多因素综合研究海岸河口水域的上升流现象。
[1]顾润源,刘晓东,白桦.厄尔尼诺研究的进展和现状[J].宁夏大学学报,1990,11(3):83-91.(GU Runyuan,LIU Xiaodong,BAI Hua.Advances and present conditions in El Nino[J].Journal of Ningxia University,1990,11(3):83-91.(in Chinese))
[2]张月丛,孟宪锋.厄尔尼诺和拉尼娜的成因及其对中国气候的影响[J].承德民族师专学报,2005,25(2):59-60.(ZHANG Yuecong,MENGXianfeng.El Nino and La Nina weather to China the causes and the influence of climate[J].Journal of Chengde Teachers’College for Nationalities,2005,25(2):59-60.(in Chinese))
[3]吕新刚.黄东海上升流机制数值研究[D].青岛:中国科学院海洋研究所,2010.
[4]赵保仁,李徽翡,杨玉玲.长江口海区上升流现象的数值模拟[J].海洋科学集刊,2003,45:64-76.(ZHAO Baoren,LI Huifei,YANG Yuling.Numerical simulation of upwelling in the Changjiang River Mouth area[J].Studia Marina Sinica,2003,45:64-76.(in Chinese))
[5]赵保仁,任广法,曹德明,等.长江口上升流海区的生态环境特征[J].海洋与湖沼,2001,32(3):327-333.(ZHAO Baoren,REN Guangfa,CAO Deming,et al.Characteristics of the ecological environment in upwelling area adjacent to the Changjiang River estuary[J].Oceanologia Et Limnologia Sinica,2001,32(3):327-333.(in Chinese))
[6]张春桂,陈家金,谢怡芳,等.利用MODIS多通道数据反演近海海表温度[J].气象,2008,34(3):30-36.(ZHANG Chungui,CHEN Jiajin,XIE Yifang,et al.Retrieval of sea surface temperature in coastal region by using MODISmulti-channel data[J].Meteorological Monthly,2008,34(3):30-36.(in Chinese))
[7]张春桂,张星,曾银东,等.台湾海峡海表面温度的遥感反演及精度检验[J].海洋学报,2008,30(2):153-160.(ZHANG Chungui,ZHANG Xing,ZENG Yindong,et al.Retrieval and validation of sea surface temperature in the Taiwan Strait using MODISdata[J].Acta Oceanologica Sinica,2008,30(2):153-160.(in Chinese))
[8]刘良明,周军元.MODIS数据的海洋表面温度反演[J].地理空间信息,2006,4(2):7-9.(LIU Liangming,ZHOU Junyuan.Using MODISimagery to map sea surface temperature[J].Geospatial Information,2006,4(2):7-9.(in Chinese))
[9]郑嘉淦,李继龙,杨文波.利用MODIS遥感数据反演东海海域海表温度的研究[J].海洋渔业,2006,28(2):141-146.(ZHENG Jiagan,LI Jilong,YANG Wenbo.Retrieval of SST in the east China sea with MODISremote sensing data[J].Marine Fisheries,2006,28(2):141-146.(in Chinese))
[10]刘玉洁,杨忠东,王长耀,等.MODIS遥感信息处理原理与算法[M].北京:科学出版社,2001:293-294.
[11]叶安乐,李凤岐.物理海洋学[M].青岛:青岛海洋大学出版社,1992:252-280.
[12]FISSEL D B,MARKO J R,Álvarez M S.The changing met-ocean and ice conditions in the Beaufort Sea:implications for Offshore Oil and Gas[C]//GHUNG JS.Annual meeting of the International Society of Offshore and Polar Engineers(ISOPE).Osaka:[s.n.],2009:1-8.
[13]HUDAK D R,YOUNG M C.Storm climatology of the southern Beaufort Sea[J].Atmosphere-Ocean,2001,40(2):145-148.