太湖表面温度与表层水温差异及其影响因素的观测分析*

李雨竹,贾 磊,张 弥,2,肖 薇,2,张 圳,包恒鑫,葛 培,王 伟,2，吕 恒，万梓文

(1:南京信息工程大学气候与环境变化国际合作联合实验室大气环境中心,南京 210044) (2:南京信息工程大学江苏省农业气象重点实验室,南京 210044) (3:南京市江宁区气象局,南京 211199)

湖泊是气候变化的指示器[1],水温升高是湖泊对气候变暖最直接的响应[2-3]。原位观测和遥感反演显示全球湖泊水温以0.34℃/(10 a)的速率升高[4-5],增暖趋势在未来高排放情景下愈加显著[6-7]。湖泊升温改变了湖泊的物理、化学、生物和生态过程。首先,作为大气模型的下边界,湖面温度升高改变了湖泊与大气之间的温度和湿度梯度,影响湖气界面的感热和潜热交换[8]。作为湖体温度垂直分布的初始值[9],湖面温度升高改变了水温垂直分层,影响湖体内热量和动量的垂直交换[10-11]。其次,湖泊升温使得深水湖泊的热力层结增强、垂直混合减弱,深水区溶解氧下降和底栖物种代谢减速[12-13]。而且,湖泊升温促进CH4排放和蓝藻异常增殖。观测显示水温高时湖气界面CH4排放量是水温低时的5.2倍[14]。据估计,到21世纪末,北方湖泊CH4排放量最大增幅可达54%[15],蓝藻暴发频率将增加20%[16]。为评估气候变暖对湖泊物理、化学、生物和生态过程的影响,需要在不同时间尺度上观测湖泊温度。

目前观测湖泊温度的方法主要有原位观测和红外遥感反演,但两者获取的水温因测量深度不同存在本质差异。原位观测获取的是表层水温[5],即利用温度传感器在湖水表层一定深度处(约0.5 m)直接测量水温。而红外遥感反演获取的是湖泊表面温度[17],即基于卫星传感器的热红外数据通过辐射传输方程等算法反演得到湖泊与大气界面薄层(<500 μm)的温度。原位观测持续时间短、覆盖范围小[18-19],红外遥感反演能拓展湖面温度研究的时空尺度[20]。常用于反演湖面温度的卫星传感器有MODIS的Terra和Aqua、NOAA发布的AVHRR传感器、欧洲空间局设计的ATSR系列传感器、Landsat8的TIR传感器、Landsat7的ETM+传感器和Aqua卫星上的ASTER传感器[21],中国发射的风云(FY)系列卫星MERSI传感器和VIRR传感器[22]、环境资源卫星(HJ)的IRS传感器[23]也常被用于反演湖面温度。上述卫星传感器的参数见表 1。

红外遥感反演得到的湖面温度与原位观测得到的表层水温在物理过程、数值大小和影响因素上存在差异[24-25]。在湖气界面薄层内,湖面温度主要受分子热传导控制[26],而表层水温主要受水体湍流和太阳辐射共同控制[18]。由于物理过程不同,表面温度与表层水温的数值存在差异[27],称为“皮肤效应”(skin effect)。“皮肤效应”的大小通常以表面温度与表层水温的差值表示,根据数值正负,分为暖“皮肤效应”(表面温度高于表层水温)和冷“皮肤效应”(表面温度低于表层水温)[28-29]。“皮肤效应”主要受湖水吸收的太阳辐射、湖面与大气之间的热量交换和湖水湍流混合3个过程控制,存在明显的日变化。白天,“皮肤效应”主要取决于太阳辐射对上层水体的加热作用,通常表现为暖“皮肤效应”;夜间,“皮肤效应”由水气界面的感热、潜热交换和净长波辐射共同决定,表现为冷“皮肤效应”[27,29-30]。“皮肤效应”受气象要素影响显著,云量和风速分别通过改变太阳辐射强度和湍流混合强度影响“皮肤效应”[31-32]。

“皮肤效应”给湖泊温度响应气候变暖研究带来诸多挑战。首先,多数研究采用表层水温观测数据验证遥感反演得到的表面温度[33-34],虽然部分研究进行了表层水温至表面温度的转换[35-36],但由于湖泊地理位置、形态特征和气候条件存在差异,导致转换方程难以推广,且转换过程会带来约0.3 K的不确定性[18]。其次,为了量化湖泊变暖速率,通常采用多源遥感数据融合来延长时间序列[17]。由于卫星过境时间不同(表1),若不考虑“皮肤效应”的日变化,仅采用固定数值将遥感反演得到的表面温度转换为表层水温势必带来误差[37]。而且,研究发现若不考虑“皮肤效应”,仅用表层水温计算湖泊感热和潜热通量存在10%高估[38]。因此,量化湖泊“皮肤效应”及其影响因素,不仅有助于准确验证遥感反演结果和计算湖气热量交换,还可更加科学地评估湖泊变暖速率。为此,本文提出以下两个科学问题:(1)在不同的时间尺度和天气条件下,湖泊“皮肤效应”特征如何?(2)若不考虑“皮肤效应”,表面温度和表层水温表征的湖泊升温速率有多大差异?

表1 卫星传感器过境时太湖表面温度(Ts)与20 cm处水温(T20)和50 cm处水温(T50)的差值Tab.1 The difference between the surface temperature (Ts) and water temperature at 20 cm depth (T20) or at 50 cm depth (T50) when satellites pass over Lake Taihu

太湖中尺度通量网[39]拥有10余年(2010年6月至今)的多站点、高频且连续的表面温度和水温梯度观测数据,为回答以上两个科学问题提供了观测数据支撑。本文基于太湖中尺度通量网2011-2020年的辐射四分量、水温梯度和小气候观测数据,量化了不同时间尺度(日、季节与年际)上、不同天气条件(晴天大风、晴天小风、阴天大风和阴天小风)下太湖表面温度与表层水温的差异,以期为太湖物理、化学、生物和生态过程对气候变暖的响应研究提供参考。

1 数据与方法

1.1 研究区域

太湖位于亚热带季风气候区,1981-2010年平均气温为16.2℃,年总降水量为1122 mm,年平均风速为2.6 m/s[40]。太湖是我国第三大淡水湖,水域面积约2338 km2,平均水深为1.9 m[41]。太湖中尺度通量网(图1)始建于2010年6月,截止2020年底,该通量网一共有7个湖泊站点,按照建成顺序依次为梅梁湾站(MLW)、大浦口站(DPK)、避风港站(BFG)、小雷山站(XLS)、平台山站(PTS)、东太湖站(DTH)与新梅梁湾站(MLW2)。

图1 太湖中尺度通量网观测平台:左图为太湖中尺度通量网站点分布,右图为避风港站(BFG)(四源净辐射传感器(上)和水温梯度观测系统(下)实景,改自Lee等[39]Fig.1 The Taihu Eddy Flux Network and instrument platform at the BFG site

1.2 观测数据与质量控制

由于BFG站观测时间长(始于2011年12月15日)、数据质量优[40-42],本文以BFG站2012-2020年0.5 h的20 cm水温(T20)、50 cm水温(T50)、辐射四分量和小气候观测资料为基础,选择MLW站(建于2010年6月)2011年观测数据为补充,并以距离BFG站最近的东山站(站号:58358)2011-2020年的气温观测数据(http://data.cma.cn)为参考。BFG站2014年水温梯度数据缺测最少(<6%)[40],故选取BFG站2013年12月1日-2014年11月30日一整年的观测数据,来分析太湖各层水温的日变化、季节变化以及不同天气条件下的“皮肤效应”。

BFG站小气候系统架设高度为8.5 m,由温湿度传感器(型号:HMP155A;公司:Vaisala)和风速风向计(型号:05103;公司:R.M.Young)组成,分别用于观测空气温度、湿度、风速和风向。水温梯度观测系统(型号:109-L;公司:Campbell Scientific)用于观测水面以下20、50、100和150 cm的水温以及底泥温度。四源净辐射传感器(型号:CNR4;公司:Kipp &Zonen B.V)架设高度为2.5 m,用于测量向下短波辐射、向上短波辐射、向下长波辐射和向上长波辐射[42-43]。参考Zhang等[40]的方法对BFG站的水温和辐射数据进行以下质量控制和缺测插补。首先,根据BFG站的气候条件设定各变量的阈值,水温变化范围为-5～45℃,向下短波辐射(K↓)变化范围为0～1200 W/m2,向上短波辐射(K↑)变化范围为0～200 W/m2,向上长波辐射(L↑)变化范围为300～600 W/m2,向下长波辐射(L↓)变化范围为150～550 W/m2,剔除超出阈值范围的异常点。其次,运用五点滑动平均两倍标准差方法剔除异常值。数据插补分两步,若缺测时长短于1 h,进行线性内插;若缺测时长超过1 h,采用邻近站点的相同变量观测值插补,即将BFG站不同年份、不同季节的原始数据与其它站点的原始数据进行回归,选择相关系数最高的站点,利用该回归方程对BFG站缺失数据进行插补。质量控制和缺测插补完成后,再运用五点滑动平均两倍标准差方法剔除异常值,并对剔除后的数据进行插补,直至没有空缺。BFG站2013年12月1日-2014年11月30日水温和辐射数据插补后的时间序列如图2所示。

图2 BFG站2013年12月1日-2014年11月30日20 cm处水温T20(a)、50 cm处水温T50(b)、向上短波辐射K↑(c)、向下短波辐射K↓(d)、向上长波辐射L↑(e)和向下长波辐射L↓(f)原始数据(蓝色)与插补数据(红色)的时间序列Fig.2 The time series of raw data (blue) and gap-filled data(red) of T20(a),T50(b),K↑(c),K↓(d),L↑(e) and L↓(f) from December 1, 2013 to November 30, 2014 at the BFG site

1.3 计算方法

1.3.1 湖面温度的计算 基于四源净辐射传感器观测的向上长波辐射L↑和向下长波辐射L↓,根据Stefan-Boltzmann定律计算湖泊表面温度Ts,即:

(1)

式中,ε为湖面比辐射率,取值为0.97[44];σ为Stefan-Boltzmann常数,为5.6696×10-8W/(m2·K4)。

从观测设备维护和数据交叉验证两个方面保证湖面温度数据的可靠性。首先,每月对太湖中尺度通量网各站点的四源净辐射传感器进行镜头擦拭与水平校准,每年将生产商校准过的同款四源净辐射传感器安装在各站点平行观测10天,对比显示测量偏差小于1.0 W/m2。其次,对比BFG站建站第一年红外测温传感器(型号:SI-411;公司:Apogee Instruments Inc.)观测结果与式(1)计算的湖面温度,两者的相关系数和一致性系数均超过0.99,平均绝对误差为0.1 K。依据表层水温观测规范[45]和太湖已有观测试验[46],结合本研究水温梯度观测实际,以湖面温度(Ts)和20、50 cm深度处水温(T20、T50)的差值Ts-T20和Ts-T50来量化太湖水温的“皮肤效应”。

1.3.2 天气状况的确定 本文采用晴空指数(Kt)来区分晴天和阴天。晴空指数指某一太阳高度角下地表水平面上接受的太阳总辐射S(W/m2)与大气上界水平面上的天文辐射Se(W/m2)的比值[47]。

(2)

Se=Ssc[1+0.033cos(360td/365)]sinβ

(3)

式中,Ssc为太阳常数,取1367 W/m2;td为日序;β为太阳高度角。定义Kt> 0.6时为晴天,Kt< 0.3时为阴天[47]。

基于BFG 0.5 h风速观测数据,采用分位数法定义大风和小风。风速高于75%分位数时(Ws> 5.8 m/s)为大风,低于25%分位数时(Ws<3.2 m/s)为小风。

2 结果与讨论

2.1 太湖水温多尺度时间变化特征

2.1.1 太湖表面温度和表层水温的日变化 太湖BFG站四季和年平均的Ts、T20和T50皆呈单峰型日变化(图3)。从年平均日变化来看(图3e),夜间湖面向上发射的长波辐射大于大气逆辐射,长波辐射冷却使得与大气直接接触的表面温度迅速降低,低于表层水温,表现为冷“皮肤效应”,7:00-8:00最强,此时Ts分别比T20和T50低0.50和0.34℃。白天,太阳辐射加热使湖面迅速升温,且升温速率快于表层水温,使得表面温度逐渐高于表层水温,表现为暖“皮肤效应”,15:00-16:00最强,此时Ts分别比T20和T50高1.54和1.95℃。此外,由于表层水温白天升温和夜晚降温速率均慢于表面温度,其日变化曲线滞后于表面温度,且随着深度增加,日变化幅度减小,滞后愈加明显。

图3 2013年12月-2014年11月BFG站表面温度(Ts)、20 cm处水温(T20)、50 cm处水温(T50)和气温(Ta)的平均日变化(a:春季;b:夏季;c:秋季;d:冬季;e:年平均)(图中阴影表示各卫星传感器过境太湖的时刻)Fig.3 Diurnal variations in the surface temperature (Ts), the water temperature at 20 cm depth (T20), the water temperature at 50 cm depth (T50) and the air temperature (Ta) at the BFG site from December 2013 to November 2014 (a: spring; b: summer; c: autumn; d: winter; e: annual mean) (The passing over time of each satellites are denoted by the gray bars)

研究卫星在太湖过境时刻的“皮肤效应”,不仅可以量化表层水温验证遥感反演表面温度结果时的误差,还可以为多源遥感产品融合构建太湖水温长时间序列提供参考。常用于表面温度反演的卫星热红外传感器在太湖过境时的Ts与T20、Ts与T50的差值如表 1所示。以年平均日变化为例,AVHRR日间过境时(14:00),暖“皮肤效应”显著,Ts-T20和Ts-T50分别为1.38和1.73℃;凌晨过境时(2:00)呈冷“皮肤效应”,强度明显降低,Ts-T20和Ts-T50分别为-0.42和-0.32℃。MODIS-Aqua凌晨过境(01:30)和午后(13:30)过境时,太湖水温呈现与AVHRR类似的“皮肤效应”昼夜对比特征。01:30时Ts-T20和Ts-T50分别为-0.40和-0.29℃;13:30时Ts-T20和Ts-T50分别为1.35和1.64℃。其他卫星多在10:00-10:30过境,此时太湖水温暖“皮肤效应”数值约为峰值的一半。如Landsat的TM/ETM+/TIRS传感器、环境资源卫星(HJ)的IRS传感器、FY-3的MERSI与VIRR传感器以及ENVISAT-AATSR过境时,Ts-T20和Ts-T50在0.58～0.75℃之间。

太湖“皮肤效应”强度在四季之间存在差异。暖“皮肤效应”在春季最强,Ts-T20和Ts-T50最大值分别为1.54和1.90℃;夏季最弱,最大值分别为1.07和1.44℃。秋冬两季介于春夏之间,秋季Ts-T20和Ts-T50最大值分别为1.30和1.81℃,冬季分别为1.59和1.77℃。可见,由于“皮肤效应”的存在,表层水温观测值无法直接验证卫星午后过境反演得到的太湖湖面温度,特别是在春季,其偏差最大接近2℃。在夜间,这种验证带来的偏差小于0.5℃。

2.1.2 太湖表面温度和表层水温的季节变化 2013年12月-2014年11月BFG站Ts、T20和T50月均值的季节变化如图4所示。太湖水温呈夏高冬低的季节变化特征,水温在春季开始升高,7月达到峰值,月均值接近30℃,其后水温下降,1月降至最低,月均值约为6℃。

图4 2013年12月-2014年11月BFG站表面温度(Ts)、20 cm处水温(T20)、50 cm处水温(T50)和气温(Ta)的季节变化箱式图Fig.4 Box plot of seasonal variation in surface temperature (Ts), water temperature at 20 cm depth (T20), water temperature at 50 cm depth (T50) and the air temperature (Ta) at the BFG site from December 2013 to November 2014

相较于日尺度,季节尺度上太湖表面温度与表层水温之间差异更小,这与太湖的浅水特性有关[44-45]。太湖水温在3和5月存在微弱的热力分层,Ts-T20和Ts-T50分别为0.65和0.71℃。除1月外,冬季Ts-T50均小于0.4℃。虽然1月气温全年最低(5.8℃)、太阳辐射较弱(110 W/m2),但该月风速(4.1 m/s)低于年平均风速(4.9 m/s),低风速抑制了湖水的动力混合,可能导致该月呈现微弱的“皮肤效应”,Ts-T20和Ts-T50分别为0.53和0.67℃。

2.1.3 太湖表面温度和表层水温的年际变化特征 太湖表面温度和表层水温在年际尺度上均呈上升趋势(图5)。10年中,Ts、T20和T50年均最低值均出现在2012年,分别为17.5、17.4和17.4℃,最高值均出现在2020年,分别为18.8、18.5和18.4℃。10年间,太湖水温升温速率随深度增加而递减。Ts升温速率最快,为0.14℃/a(P<0.01),T20和T50的升温速率相当,均为0.12℃/a(P<0.05)。对比东山站同期2 m气温可见,10年间,气温数值均低于BFG站3层水温,但上升趋势与水温一致,升温速率与Ts接近,为0.14℃/a(P<0.01)。此外,从10年间温度波动来看,太湖表面温度的标准差与气温相当,均为0.43℃,高于T20和T50的标准差(0.36℃)。可见,太湖表面温度能快速响应气象条件变化,升温速率与气温相当,但比表层水温升温速率快0.02℃/a。

图5 BFG站表面温度(Ts)、20 cm处水温(T20)、50 cm处水温(T50)和DS站气温(Ta)的年际变化(**表示年际变化趋势在0.01水平上显著,*表示年际变化趋势在0.05水平上显著)Fig.5 Interannual variations in surface temperature (Ts), water temperature at 20 cm depth (T20), water temperature at 50 cm depth (T50) at the BFG site and air temperature at the DS site

2.2 不同天气条件下的太湖水温“皮肤效应”

2.2.1 晴天与阴天 由图6可见,“皮肤效应”数值随风速增大而减小,除风速大于11 m/s外,晴天时Ts-T20和Ts-T50均较阴天时更大。晴天时,当风速小于0.5 m/s时,Ts-T20和Ts-T50数值均超过3℃,随着风速增大,Ts-T20和Ts-T50迅速减小。当风速达到2.5 m/s后,Ts-T20和Ts-T50随风速增加趋于平稳。当风速超过8 m/s 时,Ts-T20与Ts-T50的数值趋于零,此时各层湖水均匀混合,这与Yang等[31]研究结果一致。阴天时,Ts-T20与Ts-T50随风速增大而减小的趋势与晴天结果一致,但Ts-T20与Ts-T50的数值比晴天时小0.15～0.28℃。可见,晴天与阴天条件下,太湖水温的“皮肤效应”随风速的变化幅度高达3℃,与Prats等[35]的研究结果相近。而风速相同时,晴天比阴天的“皮肤效应”最多高0.5℃。

图6 BFG站晴天与阴天条件下Ts-T20(a)、Ts-T50(b)随风速的变化(▽和○分别代表晴天和阴天条件下的区间平均值,区间大小为1.5 m/s,对应的误差线为区间内观测值的标准差)Fig.6 Changes of Ts-T20 (a) and Ts-T50 (b) with wind speed under clear-sky and cloudy conditions at the BFG site (The triangle and the circle are the means of the intervals under clear-sky and cloudy conditions, respectively. The width of the interval is 1.5 m/s, and the corresponding error bar is standard deviation)

2.2.2 大风与小风 大风条件下,随着晴空指数Kt增加,Ts-T20和Ts-T50呈现先增后减的变化特征(图7a)。当Kt<0.5时,Ts-T20和Ts-T50随Kt增加而增大;当Kt接近于0.5时,达到峰值;当Kt>0.5时,温差随Kt增加而下降。小风条件下,Ts-T20和Ts-T50随晴空指数的变化趋势与大风时结果一致,但数值更大(图7b)。大风和小风时Ts-T20和Ts-T50的差异在Kt接近0.5时最大,分别为1.19和1.85℃。随着Kt偏离0.5,差异减小。当Kt接近于1时,两者无明显差异。Kt接近于0时,“皮肤效应”在大风条件下虽然接近于0,但在小风条件下依然能观测到约1℃的正“皮肤效应”。

图7 BFG站点大风与小风条件下Ts-T20(a)、Ts-T50(b)随晴空指数的变化(□和△分别代表小风和大风条件下的区间平均值,区间大小为0.2,对应的误差线为区间内观测值的标准差)Fig.7 Changes of Ts-T20(a) and Ts-T50(b) with clearness index under strong wind and weak wind conditions at the BFG site(The square and the triangle are the means of the intervals under the conditions of weak wind and strong wind, respectively. The width of the interval is 0.2, and the corresponding error bar is standard deviation)

如图8所示,太湖水温“皮肤效应”随着太阳辐射增强而增大,且增大趋势在小风天气下更加明显。当K↓由50 W/m2增至950 W/m2时,大风条件下,Ts-T20和Ts-T50均增加了1℃;小风条件下,Ts-T20和Ts-T50分别增加了1.9和2.3℃。当K↓=50 W/m2时,大风和小风条件下的“皮肤效应”差异不足0.4℃,但两者之间的差异随着K↓增强而增大;当K↓=950 W/m2时,大风和小风条件下的“皮肤效应”差异达到1.9℃。

图8 BFG站点大风与小风条件下Ts-T20(a)、Ts-T50(b)随太阳辐射的变化(□和△分别代表小风和大风条件下的区间平均值,区间大小为100 W/m2,对应的误差线为区间内观测值的标准差)Fig.8 Changes of Ts-T20(a) and Ts-T50(b) with solar radiation under strong wind and weak wind conditions at the BFG site(The square and the triangle are the means of the intervals under the conditions of weak wind and strong wind, respectively. The size of the interval is 100 W/m2, and the corresponding error bar is standard deviation)

晴天小风时Ts-T20和Ts-T50最大(图9),均值分别为1.23和1.64℃,中位数分别为0.95和1.39℃。阴天大风时Ts-T20和Ts-T50最小,均值分别为0.32和0.39℃,中位数分别为0.14和0.19℃。可见在晴天小风条件下,太湖“皮肤效应”最显著,这与以往研究结果一致[17,26-27,31,35,48]。而较弱的太阳辐射与强风速带来的湖水混合,使太湖观测不到热力分层。阴天小风下Ts-T20和Ts-T50均值分别为1.04和1.49℃,与晴天小风时结果最为接近,但强于晴天大风下的结果(0.58和0.77℃)。结合图7、8和9可知,相较于太阳辐射,风速对于太湖水温“皮肤效应”的影响更大,说明对于太湖这样的浅水湖泊而言,风的扰动是影响湖水热力分层的主要因素。Yang等[31]指出太湖水体在日尺度上的对流主要由风扰动引起,风可以驱动湖水的湍流运动,使水体上下混合均匀,从而减弱“皮肤效应”。鉴于风速对“皮肤效应”的重要性,诸多研究提出了利用风速计算“皮肤效应”的参数化方案[19,28,49-50]。但这些参数方案均源自海洋观测,会严重低估太湖水温的“皮肤效应”,在低风速时尤为明显(图10)。为此,根据太湖观测数据,本文给出了晴天条件下适用于太湖的“皮肤效应”的风速参数化方案为:ΔT= 0.4885+6.359exp(-1.173Ws),阴天条件下适用于太湖的“皮肤效应”的风速参数化方案为:ΔT= 0.1639+3.498exp(-0.7639Ws)(图10)。Ts-T50随风速变化的拟合结果与Ts-T20类似。

图9 4种天气条件下太湖表面温度与表层水温差异(Ts-T20和Ts-T50)的小提琴图Fig.9 The violin plot of the difference between the surface temperature and the bulk temperature of Lake Taihu (Ts-T20 and Ts-T50) under four different weather conditions

图10 BFG站点晴天(a)和阴天(b)条件下Ts-T20随风速变化的拟合曲线Fig.10 Fitting curves of Ts-T20 with wind speed at the BFG site under clear-sky (a) and cloudy (b) conditions

基于太湖中尺度通量网10年的表面温度和表层水温观测数据,得到太湖水温“皮肤效应”白天为正、夜间为负,这一昼夜变化特征与以往研究结果一致[18,27,35],但不同内陆水体间的“皮肤效应”数值存在差异。如法国Bimont Reservoir在上午10:00时水下55 cm处“皮肤效应”约为0.34℃[35],低于本文结果(0.86℃)。美国Lake Tahoe[24-25]、Lake Ramsey[51]的“皮肤效应”常年为负值,主要是因为表层水温观测深度较浅(2～5 cm),该薄层内水温主要受分子热传导控制,使得该观测深度处水温常年高于表面温度。作为典型的大型浅水湖泊,太湖观测的暖“皮肤效应”在午后接近2℃,可以推断深水湖泊暖“皮肤效应”会更强,利用表层水温观测值直接验证卫星反演的深水湖泊表面温度带来的误差可能更大。

太湖表面温度和表层水温的升温趋势与气温一致,尤其是表面温度(图5),这与以往研究结果一致[6,47,52-55]。太湖升温速率明显高于Sharma等[5](0.3℃/(10 a))和Schneider等[3](0.13℃/(10 a))得到的全球湖泊的平均升温速率,这一差异可能与研究时段不同有关。若将研究时段固定在2010-2019年,张运林等[46]研究发现50 cm处太湖水温升温速率为1.79℃/(10 a),高于本研究的结果(1.13℃/(10 a)),这一差异可能与水温观测站点的选取有关。张运林等[46]采用太湖栈桥水域每日2:00、14:00与20:00 3个时刻50 cm处水温观测数据,该观测点离岸约50 m,易受周边人为活动等的影响。本研究的BFG站离岸最近距离超过4 km,观测数据受周边环境影响较小。对于大型湖泊,采用近岸观测数据可能会高估湖泊升温速率,有必要在湖泊中心开展观测,以客观地表征湖泊水体的升温速率。表面温度和表层水温都可以用来量化湖泊升温速率[3,5-6,46,52-55]。太湖近10年观测数据显示其表面温度上升速率比50 cm处水温上升速率快0.02℃/a,若不考虑“皮肤效应”,直接使用遥感反演的表面温度表征的太湖升温速率会比传统的表层水温观测结果快0.02℃/a。

3 结论

基于太湖2011-2020年辐射四分量、水温梯度和小气候观测数据,量化了不同时间尺度(日、季节与年际)上、不同天气条件(晴天大风、晴天小风、阴天大风和阴天小风)下太湖表面温度与表层水温的差异,结论如下:

1)太湖四季和年平均的Ts、T20和T50皆呈现单峰型日变化。暖“皮肤效应”在15:00-16:00最强,Ts-T20和Ts-T50分别为1.54和1.95℃;冷“皮肤效应”在7:00-8:00最强,Ts-T20和Ts-T50分别为-0.50和-0.34℃。“皮肤效应”强度春季最强,Ts-T20和Ts-T50日最大值为1.54～1.90℃;夏季最弱,Ts-T20和Ts-T50日最大值为1.07～1.44℃。用表层水温观测值直接验证卫星午后过境反演得到的太湖湖面温度,其偏差接近2℃,尤其在春季。在夜间,验证偏差小于0.5℃。

2)太湖Ts、T20和T50均呈夏高冬低的季节变化特征,三层水温数值差异不显著。太湖表面温度与表层水温的在年际尺度上均呈上升趋势,其中Ts上升最快,为0.14℃/a(P<0.01),与同期气温上升速率相当。T20和T50的上升速率均为0.12℃/a(P<0.05)。直接使用遥感反演的表面温度表征的太湖升温速率会比传统的表层水温观测结果快0.02℃/a。

3)晴天小风时暖“皮肤效应”最显著,Ts-T20和Ts-T50均值分别为1.23和1.64℃;阴天大风时“皮肤效应”最不显著,Ts-T20和Ts-T50均值分别为0.32和0.39℃。相较于太阳辐射,风速对太湖水温“皮肤效应”的影响更大。现有的参数方案会低估太湖水温的“皮肤效应”,在低风速时尤为明显。基于10年观测数据建立了适用于太湖晴天和阴天条件下“皮肤效应”风速参数化方案。