水体(本文指湖泊、池塘和水库)为人类社会提供了基本的生态系统服务。水体的分布及面积变化在很大程度上受到气候变化和人类活动的影响。过去三十余年,随着中国的土地利用方式的转变,水体的面积、丰度和分布等也发生了明显变化。然而,由于水体提取的数据和方法的局限性,对中国的内陆水体特别是小水体(面积小于1 km2)的空间分布及变化依然缺乏较为全面的了解,成为我国水资源开发利用和可持续管理的瓶颈。因此,本研究将全国分成11个地区(即长江流域(YTR)、黄河流域(YR)、珠江流域(PR)、东南诸河地区(SER)、西南诸河地区(SWR)、青藏高原地区(TP)、西北诸河地区(NWR)、松花江流域(SHR)、辽河流域(LR)、海河流域(HR)、淮河流域(HuR)),利用1984年至2015年间所有Landsat遥感影像数据分析中国内陆水体(面积大于0.001 km2)的总体特征以及在过去两个时期(即1984-1999年和2000-2015年)的变化。本研究取得的结果如下:(1)中国大于0.001 km2的水体共有1 544 116个,总面积为132 662 km2;其中小水体(小于1 km2)数量为1 537 295个(占比99.5%),但面积只有20 362.69 km2(占比15.3%);中大型(大于1 km2)的水体数量为6 821个,是以往研究估计的1.2-1.7倍,水体面积为123 342 km2,比以往研究估计高出41-63%。(2)各个区域间水体的分布表现出较大的不均衡性,水体数量最多的是长江流域(538 462个),数量最少的是辽河流域(35 372个),水体总面积最大的是长江流域(37 297km2),总面积最小的为东南诸河区域(1 803 km2)。(3)中国水体的大小-丰度关系符合幂律分布的特征,帕累托系数为-0.835,与美国大陆水体的系数相近(-0.85),但高于全球水体的系数(-1.051)。在区域范围内,帕累托系数最高的是青藏高原(-0.66),表明在本区域中、大水体占的比重相对较高;最低的是珠江流域(-0.941),表明小水体在本区域占的比重较高。(4)中国水体岸线-面积的关系同样符合幂定律,岸线分形系数为1.254。区域范围内水体岸线分形系数最高的是珠江流域(1.321),表明水体岸线在珠江流域岸线形态复杂,受到人为干扰较大;岸线分形系数最低的是青藏高原区域(1.178),表明水体岸线在青藏高原较规整,在青藏高原受到人为干扰少。(5)在过去三十年间,全国共有5 197.09 km2的水体消失,约占全国水体总面积的3.92%;同时共有9 040.26 km2的非水体转化为水体,约占总面积的6.81%。全国范围内超过一半的水体在研究时期内保持稳定,其余水体均显示出不同程度的变化。(6)区域尺度上水体增加和减少同时存在。水体变化最剧烈的地区是长江流域(主要是由于人类活动)和青藏高原(主要是由于气候变化),前者水体表现为急剧的增加或减少,后者水体表现为急剧的增加。总体来看,本研究结果表明了(1)中国水体总体特征为水体在全国和地区上分布不均衡,差异较大,虽然中国水体的大小-丰度以及岸线长度-面积的关系符合幂律分布,但由于受到气候和土地利用方式的影响,分形特征在区域尺度上表现出较大的变异性。(2)中国水体的变化呈现出较为明显的区域特征,这主要与各区域间气候变化和土地利用方式紧密相关。本研究突出了利用长期卫星数据来研究大范围内水体动态特征的重要性,有利于在更小区域或复杂地形区验证水体分形理论的适用性,推动区域乃至全球尺度的生物地球化学过程,估计全球碳交易以及预测气候变化的轨迹,促进中国和世界其他地区水资源的保护及动态研究。