,
,
,
,
,
,
2. 中国科学院青藏高原地球科学卓越创新中心 北京 100101;
3. 中国科学院大学 北京 100049;
4. 兰州大学 资源环境学院 兰州 730000;
5. 清华大学 地球系统科学系 北京 100080;
6. 郑州大学 地球科学与技术学院 郑州 450001;
7. 暨南大学 生态学系 广州 510632
, , , , , ,
2. cas center for excellence in tibetan plateau earth sciences, beijing 100101, china;
3. university of chinese academy of sciences, beijing 100049, china;
4. college of earth and environmental sciences, lanzhou university, lanzhou 730000, china;
5. department of earth system sciences, tsinghua university, beijing 100080, china;
6. school of geoscience and technology, zhengzhou university, zhengzhou 450001, china;
7. department of ecology, jinan university, guangzhou 510632, china
青藏高原分布着众多的冰川、湖泊, 是大江大河的发源地, 也是南亚、中亚水资源的重要补给保障[1], 被称为"亚洲水塔" [2]。青藏高原的湖泊众多, 其中面积大于1 km2的湖泊数量1 171个, 总面积约46 500 km2, 占青藏高原土地面积的1.8%, 超过全国湖泊总面积的50%[3]①。青藏高原地区98个气象站点观测数据表明, 1982-2012年平均气温升高幅度达到1.9℃[4], 是全球平均升温幅度的2倍。与此同时, 青藏高原降水量也呈现增加趋势[5], 但风速呈减弱趋势[6], 潜在蒸发量有所下降[7]。作为"亚洲水塔"的重要组成部分, 青藏高原广泛分布的湖泊是地表水汇聚和蒸发过程的重要环节, 并且通过湖泊面积、水量及其理化性质和生态条件对气候变化呈现敏感响应。
① 2014年数据，未统计已经作为矿产开采的盐湖。
1 湖泊数量、面积与水量变化 1.1 湖泊数量与面积变化随着全球变暖, 青藏高原湖泊呈现明显的数量增加和面积增大趋势, 80%以上的湖泊在扩张。王苏民和窦鸿身[8]利用20世纪60-80年代地形图研究发现, 青藏高原大于1 km2的湖泊有1 081个(总面积4.5×104 km2)。zhang等[9]利用landsat影像数据, 对青藏高原近40年来的湖泊数量和面积变化及其与气温和降水关系进行了详细研究。结果发现, 1976年面积大于1 km2的湖泊数量为1 080个, 总面积为4.0×104±766.5 km2; 1995年减少到930个, 面积下降了5.6%;到2000年又增加到1 174个, 总面积为4.1×104±443.7 km2; 至2018年数量达到1 424个, 总面积高达5.0×104±791.4 km2。20世纪70年代-2018年, 青藏高原湖泊总面积增长了25.4%, 但趋势不是均一的, 具有20世纪90年代中期的低值, 2000-2010年的快速增长, 2010-2016年的缓慢增加, 而近2年又出现扩张加剧的现象(图 1)。湖泊数量的增加主要是由于湖泊面积的扩张, 即湖泊面积从小于1 km2增加到大于1 km2, 而完全新出现(即从无到有)的湖泊很少, 它们对增加的湖泊数量贡献很小[10]。
进入21世纪以来, 卫星雷达与激光测高技术不断发展并产出了一批青藏高原地区的湖泊高程数据。 2003-2009年, 青藏高原约200个有可利用的icesat激光测高数据的湖泊表明[11], 湖泊平均的水位变化率为每年0.14 m, 包括152个(占湖泊个数的76 %)水位升高的湖泊(平均变化率为每年上升0.21 m)和48个(占湖泊个数的24 %)水位下降的湖泊(每年下降0.08 m) (图 2)。面积较大的色林错显示了水位快速升高, 根据色林错的水量平衡关系估算, 该湖1979-2017年水位累计上升14 m [12], 面积则由1667 km2增加到2 389 km2, 超过纳木错的面积(2 026 km2), 而成为西藏目前最大的湖泊[13]。在空间分布上, 青藏高原中-北部内流区湖泊水位显示明显升高, 而在南部的雅鲁藏布江流域, 湖泊水位以下降为主。
结合湖泊的水位变化及面积数据, 可对湖泊的水量变化进行估算。 yang等[14]利用"航天飞机雷达地形任务数字高程模型"(shuttle radar topography mission, digital elevation model, srtm dem)高程数据和landsat影像数据获取的湖泊面积, 建立了青藏高原地区面积大于50 km2的湖泊面积与水量变化关系, 并通过不同大小的湖泊实地测深结果验证了方法的可靠性, 据此估算了114个封闭湖泊1976- 1990年、1990-2000年、2000-2005年和2005- 2013年4个时期湖泊水量变化, 发现1976-2013年, 大于50 km2的全部湖泊水量共增加了1 026.4亿立方米, 并且在2000-2005年期间增速最大。qiao等[15]进一步将研究对象扩大到大于10 km2的315个湖泊, 发现这些湖泊在1976-2013年期间的水量共增加了1 171.1亿立方米, 尽管进一步考虑了201个面积介于10-50 km2的中小湖泊水量的增加, 这些湖泊的水量增量仅为大于50 km2的大中型湖泊水量增量的14.1%, 说明青藏高原大中型湖泊水量增加控制了其总体变化的态势。
2 湖水基本物理化学性质和浮游生物变化 2.1 湖水基本物理化学性质变化水温和盐度是湖泊对气候变化响应的重要物理化学指标。对纳木错水温与气温变化的观测和模拟研究发现, 1979-2012年湖水表层夏季平均水温升高率为(0.52℃±0.25℃)/10 a, 湖水温跃层(湖水上部具有水温季节变化的薄层与下部水温相对稳定的厚层之间出现水温急剧下降的层)的分层开始日期以(4.20±2.02) d/10 a的速率提前, 而分层的持续时间以(6.00±3.54) d/10 a的速率递增, 因子相关分析表明, 气温升高和湖水接受的长波辐射是造成湖泊水温上升、温跃层提前和持续时间延长的主要原因[16]。对色林错最近40年的盐度变化比较显示, 湖水矿化度已经从1979年的18.5 g/l下降到2017年的12.4 g/l, 而大量的以降水为主的淡水补给则是造成湖泊盐度下降的主要原因[13]。
湖泊透明度不仅是湖泊水体的基本物理化学指标, 也对湖泊水生生态环境具有重要影响。对青藏高原24个具有不同透明度变化的湖泊研究表明, 湖泊水色与其透明度具有很好的相关性[17]。影响湖泊水色的3个要素是湖水中的浮游植物叶绿素、有机黄色物质和无机悬浮物。青藏高原地区气候环境恶劣, 大多数湖泊基本不受人类活动影响, 湖泊生态系统生产力与内外源有机质输入微弱。作者调查实测数据表明, 该地区大部分湖泊叶绿素a浓度仅为0.1-5 μg/l, 有机质含量和悬浮物浓度仅为0.1-10 mg/l。因此, 影响青藏高原湖泊水色, 进而决定其透明度变化的主要因素是溶解性物质及其浓度。利用水色-透明度反演模型[18]重建的2000-2017年青藏高湖泊透明度变化结果表明, 在所研究的152个湖泊中, 有91个表现出明显的透明度增加(增加率为0.09 m/a), 而其他61个湖泊呈现透明度下降(下降率为−0.04 m/a), 这说明整个青藏高原的湖泊透明度在2000-2017年具有增加的态势(图 3), 而湖泊透明度增加与降水的变化具有显著的正相关关系。
|
| 图 3 青藏高原面积大于50 km2的湖泊2000-2017年的透明度(以sd深度计算)变化 蓝色圆圈表示湖泊sd深度增加, 红色圆圈表示湖泊sd深度下降; 圆圈的大小对应于变化率的绝对值 |
湖泊的浮游生物对湖水的温度和盐度变化具有明显响应。根据2012-2015年夏季对青藏高原49个湖泊的调查结果, 从浮游植物的种类组成可以将湖泊聚为4个大类:第一类湖泊浮游植物生物量相对高(1.3-4.2 mg/l), 优势种类以个体较大的硅藻为主, 具有水体盐度较低(< 5.5 g/l)和水温相对较高(>17℃)的特征; 第二类湖泊浮游植物生物量很低(0.06-0.09 mg/l), 优势种类以小型硅藻为主, 具有盐度相对较高(20-43 g/l)和水温相对较高(>16℃)的特征; 第三和第四类湖泊的浮游植物生物量中等(0.02-0.43 mg/l), 其盐度(0.22-6 g/l)和水温(13℃-15℃)也相对适中, 但第三类以皮克(pg)级浮游植物为主要优势种类, 第四类以鞭毛藻类为主要优势类群。不同盐度类型的湖泊浮游动物种类组成差别比较大:枝角类分布的盐度范围为0- 27.5 g/l; 桡足类分布的盐度范围为0-46 g/l; 无甲类分布在盐度范围为25-76 g/l的湖泊中。盐度超过100 g/l的湖泊中则未发现浮游动物分布。
青藏高原的湖泊水温上升, 对食物链各营养级间的相互作用产生了影响。在只有浮游动物和浮游植物2个营养级的湖泊中, 水温的上升促进了食物链的传递效率, 更加有利于浮游动物发展, 出现浮游动物与浮游植物生物量比值较高和浮游植物生物量与总磷比值较低; 相反, 在具有鱼类、浮游动物和浮游植物3个营养级的湖泊中, 水温的上升促进鱼类发展, 加剧了鱼类对浮游动物的捕食压力, 造成浮游动物与浮游植物生物量比值较低和浮游植物生物量与总磷比值较高。湖泊盐度升高使得浮游植物生态系统结构简单化, 浮游动物由小型枝角类、轮虫和桡足类等逐渐演替为大型滤食性枝角类占据优势; 随着盐度继续上升, 大型滤食性枝角类又被无甲类所取缔。因此, 高盐度湖泊中的敞水区往往有利于大型滤食性枝角类(西藏溞)或无甲类(卤虫)生存。湖泊盐度变化对浮游动物多样性的影响强度取决于食物链的长度, 2个营养级湖泊受影响程度比3个营养级湖泊更为强烈。在盐度为3-5 g/l或25-28 g/l的湖泊中, 浮游动物种类组成对盐度变化最为敏感[19]。
3 区域气候变化时空差异对湖泊主要补给的影响 3.1 湖泊主要补给的区域差异对气候变化响应近期青藏高原湖泊扩张和水量增加的原因主要是降水和冰川融水增加以及蒸发减少等, 而冰川融水与蒸发变化则主要受气温变化影响。气温持续升高使得冰川、冻土等加速融化[20], 但气温升高对青藏高原不同地区湖泊变化的影响各不相同。在青藏高原西北部, 冰川融水增多是湖泊扩张和水量增加的主导因素[21]。在羌塘地区东南部, 冰川融水对湖泊扩张也具有明显影响[22];而对纳木错水量变化的定量分析表明, 冰川融水对湖泊水量增量的贡献率为52.9%[23]。在青藏高原腹地, 依布茶卡和色林错水位的上升均与流域内积雪面积变化显著相关[24], 而冻土退化释放的水量对湖泊变化也有重要影响[25]。相较于流域内无冰川、冻土及积雪分布的湖泊, 气温主要通过影响蒸发从而改变湖泊的水量平衡[26, 27]。在青藏高原西南部的玛旁雍错流域, 1974-2003年冰川面积从107.92 km2减少到100.39 km2, 但区域气候暖干化使得年降水量减少、蒸发量增大, 冰川加速融化带来的补给并未能使湖泊发生明显扩张[28]。
对于青藏高原大多数地区的湖泊, 区域降水变化是造成湖泊扩张的主要原因。通过对青藏高原大于20 km2的109个内流封闭湖泊面积遥感估算和部分湖泊水位调查, 发现这些湖泊1976-1999年的变化各不相同, 而1999-2010年出现了明显的扩张。尽管潜在蒸发减少和冰川融水增加可能有助于湖泊的扩张, 但统计发现这些湖泊扩张的主要原因还是区域降水增加[29]。通过对大于50 km2的114个封闭湖泊的水量年变化速率与气象因子的相关分析, 发现1990年以前, 低温通过抑制融水发生使得湖泊水量趋向负平衡; 1990-2000年, 气温升高对湖泊水量的增加影响有限; 2000-2005年, 湖泊水量呈现显著增加的主要因素是较高的降水; 2005-2013年, 强烈蒸发对湖泊水量增加的抑制作用超过了由高温带来的融水增加的影响[14]。
3.2 湖泊变化与西风-印度季风环流及地表水循环的关系基于湖泊补给过程研究的湖泊群变化能够敏感地反映区域地表水循环的特征。在气候条件相似的青藏高原东北部, 那曲地区和可可西里地区的湖泊总体呈扩张趋势, 而黄河源区的湖泊则总体呈萎缩状态[30]。1999-2010年, 青藏高原大于20 km2的109个内流封闭湖泊均呈现明显的扩张, 但具有北部扩张明显、南部扩张微弱的空间分布特征, 其原因可能与该时期印度季风区降水减少和西风区降水增加有关[31]。进一步对其中典型湖泊水位监测和cryosat卫星测高数据分析, 结合大气降水和重力卫星grace质量变化, 认为青藏高原中北部和东北部湖泊快速扩张主要受夏季降水显著增加的影响, 而青藏高原西北部湖泊快速扩张更多与冰川消融和春节积雪增加的变化有关[32]。这一结论也得到对青藏高原西北部34个湖泊面积和水位变化研究结果的支持[20]。
通过对青藏高原大于50 km2的114个封闭湖泊在20世纪70年代-2013年的面积与水量变化研究, 发现湖泊水量随时间变化呈现3种类型且有明显的区域分异:①羌塘高原东南部和西北部以及阿里地区的湖泊变化主要受降水增加的影响, 呈现由平稳或缓升向急升转变的趋势; ②羌塘高原中部和昆仑山北部的湖泊水量受降水和温度升高导致的冰川融水增加的共同影响, 呈现先下降然后快速上升的趋势; ③而藏南湖泊的水量变化具有波动下降态势, 与温度、降水的关系并不明显(图 4) [14]。因此, 在降水主导的青藏高原湖泊变化的整体态势中, 湖泊的时空变化特征可能敏感地反映了西风和印度季风降水的差异, 而由于部分地区冰川融水的明显影响, 湖泊变化并不完全受降水变化的控制, 这说明西风和印度季风环流对地表水循环的影响过程和机制更为复杂。
青藏高原的湖泊变化受降水、蒸发以及气温增加引起的冰川融水变化影响, 而湖泊蒸发又与表层水温和风速具有密切联系[12]。这些因素在不同地区的影响程度具有很大差异[15]。对纳木错的表层水温监测与模拟表明, 湖泊表层水温随着气温升高和长波辐射增加呈现明显的增加趋势[16]。在考虑上述要素的湖泊水量平衡分析中, 通过湖泊水量变化, 反推了引起湖泊水量变化的气象要素变化[33]。结果表明, 内陆封闭湖泊面积从1995年的24 930 km2增加到2015年的33 741 km2(增幅达35%)的过程中, 反推的降水增幅达21%±7%, 与"全球降水气候中心"(global precipitation climatology centre, gpcc)的结果[34]具有高度一致性(图 5)。进一步根据"跨领域影响模型比较项目"(the inter-sectoral impact model intercomparison project, isimip) [35]预测的降水分析: 2016-2025年, 气候变化速率可能与现在近似的情况下, 湖泊面积将继续增加4 000 km2; 而2026-2035年, 由于气候的温暖湿润程度更加强烈, 湖泊可能出现更强的扩张。
4.2 青藏高原湖泊变化研究需要关注的问题
青藏高原的湖泊变化对气候变化具有非常敏感的响应, 不仅通过不同类型的补给变化反映气候变化及各要素的影响过程, 而且在湖泊理化性质和生态条件上产生一系列连锁响应。因此, 开展更加精细的湖泊水量变化与水量平衡、湖泊盐度和水生生态系统变化研究, 有助于深入理解气候变化对湖泊变化的影响程度, 从而准确评判未来气候变化条件下的湖泊变化趋势。
(1)宏观尺度的湖泊水量赋存和水量平衡研究。青藏高原的湖泊有多少水量赋存, 气候变化对其赋存条件和程度具有哪些影响, 是湖泊研究回答其在"亚洲水塔"中的作用的最根本问题。过去几十年遥感技术与数据的发展, 使得利用卫星遥感数据已经较为容易和快捷地获取湖泊的面积指标, 但由于湖盆形状和岸线坡度的差异, 湖泊面积相同并不等于水量赋存接近, 不同湖泊之间的面积变化幅度与其水量变化幅度也不一致。湖泊水量变化是气候变化带来的实质结果, 因此需要根据湖盆数字高程模型(dem)对每一个湖泊建立符合其本身特点的面积和水量关系, 才能通过遥感大数据获取整个高原各个湖泊不同时期的水量赋存, 分析其补给过程的时空变化, 进而准确评价整个青藏高原湖泊水量变化与气候变化的关系。
(2)湖泊主要理化性质与典型湖泊生态系统的系统调查与分析。青藏高原的湖泊水体质量如何, 气候变化对其有何影响, 是湖泊研究回答其在"亚洲水塔"中作为可利用水资源的核心问题。从大空间尺度来看, 获取代表湖泊主要理化性质和影响湖泊生态系统的要素是回答上述问题的关键。湖泊透明度与盐度变化对湖泊生态系统具有重要的影响。遥感影像提供了长时期、大范围的湖泊水色光谱数据, 由于湖泊水色敏感地反映湖泊透明度变化, 而水色与盐度之间具有较好的相关性, 因此具有利用遥感影像的湖泊水色数据重建湖泊透明度与盐度变化的基础。获取较高精度时间序列的湖泊透明度和盐度的变化, 对于理解湖泊生态系统的气候变化响应具有重要意义。
(3)大尺度气候变化对湖泊变化的影响及区域水循环过程。西风和印度季风是影响青藏高原地区的两大环流系统。近年来, 西风和印度季风作用区出现了明显的降水变化, 从而导致不同地区的湖泊变化具有相应的时空分异, 并使得这些湖泊与大气间的水分和热量交换发生明显改变, 进而影响了区域水循环过程。同时, 青藏高原地区的湖泊除了受降水影响外, 还与流域内的冰川融水、冻土退化等具有紧密联系。通过获取年际尺度的大范围的湖泊水量、透明度与盐度变化, 与流域尺度的气温、降水、蒸发、冰川融水、冻土变化等要素进行大数据分析, 才能发现它们在不同时间和区域上的联系, 从而评价西风和印度季风作用下的气候要素变化如何影响湖泊变化以及区域水循环过程。
| [1] |
pritchard h d. asia's shrinking glaciers protect large populations from drought stress. nature, 2019, 569: 649-654. |
| [2] |
immerzeel w w, van beek l p h, bierkens m f p. climate change will affect the asian water towers. science, 2010, 328: 1382-1385. |
| [3] |
wan w, long d, hong y, et al. data descriptor:a lake data set for the tibetan plateau from the 1960s, 2005, and 2014. scientific data, 2016, 3: 160039. |
| [4] |
chen b x, zhang x z, tao j, et al. the impact of climate change and anthropogenic activities on alpine grassland over the qinghai-tibet plateau. agricultural and forest meteorology, 2014, 189-190: 11-18. |
| [5] |
yang k, wu h, qin j, et al. recent climate changes over the tibetan plateau and their impacts on energy and water cycle:a review. global and planetary change, 2014, 112: 79-91. |
| [6] |
lin c, yang k, qin j, et al. observed coherent trends of surface and upper-air wind speed over china since 1960. journal of climatology, 2013, 26: 2891-2903. |
| [7] |
zhang y, liu c, tang y, et al. trends in pan evaporation and reference and actual evapotranspiration across the tibetan plateau. journal of geophysical research-atmosphere, 2007, 112: d12110. |
| [8] |
王苏民, 窦鸿身. 中国湖泊志. 北京: 科学出版社, 1998.
|
| [9] |
zhang g, luo w, chen w, et al. a robust but variable lake expansion on the tibetan plateau. science bulletin, 2019, 64: 1306-1309. |
| [10] |
zhang g, yao t, chen w, et al. regional differences of lake evolution across china during 1960s-2015 and its natural and anthropogenic causes. remote sensing of environment, 2019, 221: 386-404. |
| [11] |
zhang g, yao t, xie h, et al. increased mass over the tibetan plateau:from lakes or glaciers?. geophysical research letters, 2013, 40(10): 2125-2130. |
| [12] |
guo y, zhang y, ma n, et al. long-term changes in evaporation over siling co lake on the tibetan and its impact on recent rapid lake expansion. atmospheric research, 2019, 216: 141-150. |
| [13] |
zhu l p, wang j b, ju j t, et al. climatic and lake environmental changes in the serling co region of tibet over a variety of timescales. science bulletin, 2019, 64: 422-424. |
| [14] |
yang r m, zhu l p, wang j b, et al. spatiotemporal variations in volume of closed lakes on the tibetan plateau and their climatic responses from 1976 to 2013. climatic change, 2017, 140(3-4): 621-633. |
| [15] |
qiao b j, zhu l p, yang r m. temporal-spatial differences in lake water storage changes and their links to climate change throughout the tibetan plateau. remote sensing of environment, 2019, 222: 232-243. |
| [16] |
huang l, wang j b, zhu l p, et al. the warming of large lakes on the tibetan plateau:evidence from a lake model simulation of nam co, china, during 1979-2012. journal of geophysical research-atmosphere, 2017, 122: 13095-13107. |
| [17] |
刘翀, 朱立平, 王君波, 等. 基于modis的青藏高原湖泊透明度遥感反演. |
| [18] |
lee z p, shang s l, hu c m, et al. secchi disk depth:a new theory and mechanistic model for underwater visibility. remote sensing of environment, 2015, 169: 139-149. |
| [19] |
lin q q, xu l, hou j z, et al. responses of trophic structure and zooplankton community to salinity and temperature in tibetan lakes:implication for the effect of climate warming. water research, 2017, 124: 618-629. |
| [20] |
kang s, xu y, you q, et al. review of climate and cryospheric change in the tibetan plateau. environmental research letters, 2010, 5: 015101. |
| [21] |
qiao b j, zhu l p. difference and cause analysis of change in lakes of different supply types in the northwestern tibetan plateau. hydrological processes, 2017, 31(15): 2752-2763. |
| [22] |
song c, sheng y. contrasting evolution patterns between glacier-fed and non-glacier-fed lakes in the tanggula mountains and climate cause analysis. climatic change, 2016, 135(3-4): 493-507. |
| [23] |
zhu l p, xie m p, wu y h. quantitative analysis of lake area variations and the influence factors from 1971 to 2004 in the nam co basin of the tibetan plateau. |
| [24] |
zhang g q, xie h j, yao t d, et al. snow cover dynamics of four lake basins over tibetan plateau using time series modis data (2001-2010). |
| [25] |
li y, liao j, guo h, et al. patterns and potential drivers of dramatic changes in tibetan lakes, 1972-2010. plos one, 2014, 9: e111890. |
| [26] |
lei y, yao t, bird b w, et al. coherent lake growth on the central tibetan plateau since the 1970s:characterization and attribution. journal of hydrology, 2013, 483: 61-67. |
| [27] |
zhou j, wang l, zhang y, et al. exploring the water storage changes in the largest lake (selinco) over the tibetan plateau during 2003-2012 from a basin-wide hydrological modeling. water resource research, 2015, 51: 8060-8086. |
| [28] |
叶庆华, 姚檀栋, 郑红星, 等. 西藏玛旁雍错流域冰川与湖泊变及其对气候变化的响应. 地理研究, 2008, 27(5): 1178-1190. |
| [29] |
lei y, yang k, wang b, et al. response of inland lake dynamics over the tibetan plateau to climate change. climatic change, 2014, 125: 281-290. |
| [30] |
万玮, 肖鹏峰, 冯学智, 等. 近30年来青藏高原羌塘地区东南部湖泊变化遥感分析. |
| [31] |
yao t d, thompson l g, yang w, et al. different glacier status with atmospheric circulations in tibetan plateau and surroundings. nature climatic chang, 2012, 2(9): 663-667. |
| [32] |
lei y, yao t, yang k, et al. lake seasonality across the tibetan plateau and their varying relationship with regional mass changes and local hydrology. geophysical research letters, 2017, 44(2): 892-900. |
| [33] |
yang k, lu h, yue s, et al. quantifying recent precipitation change and predicting lake expansion in the inner tibetan plateau. climatic change, 2018, 147: 149-163. |
| [34] |
schneider u, becker a, finger p, et al. gpcc's new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle. theoretical apply of climatology, 2014, 115: 15-40. |
| [35] |
hempel s, frieler k, warszawski l, et al. a trend-preserving bias correction-the isimip approach. |
| [36] |
budyko m i. the heat balance of the earth's surface. washington, dc: national weather service, u. s. department of commercial, 1958: 144-155.
|